- 1Division of Crop Improvement, ICAR-Indian Grassland and Fodder Research Institute, Jhansi, India
- 2Department of Chemistry and Biochemistry, Chaudhary Sarwan Kumar Himachal Pradesh Agriculture University, Palampur, HP, India
- 3Division of Seed Technology, Indian Council of Agricultural Research-Indian Grassland and Fodder Research Institute, Jhansi, India
- 4Germplasm Evaluation Division, National Bureau of Plant Genetic Resources, New Delhi, India
- 5Division of Social Science, Indian Council of Agricultural Research-Indian Grassland and Fodder Research Institute, Jhansi, India
The study was undertaken to characterize the total phenolics, flavonoids, essential oils, quinones, tannins and antioxidant activity of 15 samples of wild Murraya koenigii (L.) Spreng. (MK) leaves obtained from different locations of Himachal Pradesh at various growth stages. The results indicated a significant variation in total phenolic content which ranged from [(170.09 ± 4.59 to 303.57 ± 7.94) in pre-flowering, (266.48 ± 7.49 to 450.01 ± 11.78) in the flowering stage, and (212.72 ± 5.37 to 363.85 ± 9.79) in fruiting stage], expressed as mg tannic acid equivalents (TAE)/g. The total flavonoid content ranged from [(15.17 ± 0.36 to 33.40 ± 0.81) in pre-flowering, (25.16 ± 0.67 to 58.17 ± 1.52) in flowering stage, and (17.54 ± 0.42 to 37.34 ± 0.97) in fruiting stage], expressed as mg catechin equivalent (CE)/g. Total tannin content ranged from [(75.75 ± 1.69 to 143 ± 3.74) in pre-flowering, (116 ± 3.26 to 207 ± 5.42) in the flowering stage, and (47 ± 1.18 to 156 ± 4.05) in fruiting stage], expressed as mg TAE/g. The essential oil content ranged from (0.64 ± 0.01 to 0.89 ± 0.02%) in pre-flowering, (0.85 ± 0.02 to 1 ± 0.02%) in flowering stage, and (0.54 ± 0.01 to 0.7 ± 0.01%) in fruiting stage. Quinones ranged from [(2.05 ± 0.05 to 2.97 ± 0.07) in pre-flowering, (3.07 ± 0.07 to 4.95 ± 0.13) in flowering stage, and (1.02 ± 0.02 to 1.96 ± 0.04) in fruiting stage], expressed as mM/min/g tissue. Antioxidant activity ranged from [(4.01 ± 0.09 to 7.42 ± 0.17) in pre-flowering, (8.08 ± 0.19 to 13.60 ± 0.35) in flowering stage, and (3.11 ± 0.06 to 6.37 ± 0.15) in fruiting stage], expressed as μg/ml. Data was subjected to multivariate analysis using principal component analysis (PCA), hierarchical clustering analysis (HCA). This was used for elucidating the intricate relationships between the phytochemical properties. All evaluated phytochemical parameters significantly increased during the growth transition from pre-flowering to the flowering stage, followed by their gradual decrease during the fruiting stage. The present study can serve as rationale for commercializing MK for aromatic and phytopharmaceutical industries.
Introduction
Earth is a plant-oriented planet and plants are adding value to earth’s diversity and are fundamental to all life (Grygorieva et al., 2021). They have been used for food, medicine or aesthetic purposes for centuries. Using novel plant-based natural products for preventing, treating diseases and enhancing health has recently gained significant scientific attention (Binici et al., 2021). Traditional plant-based medicines have a long history since the ancient civilization where plant-based materials were used as the primary ingredient for drug synthesis (Chandra and Nagani, 2009). An ethnomedicinal plant Murraya koenigii (L.) Spreng. or Curry-leaf tree displays distinct bioactivities and medicinal properties. These properties can enable its use in the food, cosmetic, and pharmaceutical industries (Verma et al., 2013; Rajendran et al., 2014; Suthar et al., 2022).
The Murraya koenigii (L.) Spreng. (MK) leaves are moderately bitter, slightly pungent and mildly acidic in nature. These leaves are nutritionally rich owing to the numerous phytochemicals, bioactives, essential oils and nutritionally relevant biomolecules (Abeysinghe et al., 2021). The MK leaves are used as appetizers, digestives, analgesics, and anthelmintics in Indian cooking (Bhandari, 2012; Desai et al., 2012). Modern research has now highlighted the potential of MK leaves as a potent antidiabetic (Yadav et al., 2002; Adebajo et al., 2005; Kesari et al., 2005), antifungal (Murugan et al., 2013), nephroprotective (Mahipal and Pawar, 2017), nitric oxide scavenging (Gul et al., 2012), antioxidant (Bhatt et al., 2020), anti-inflammatory, antibacterial (Rautela et al., 2018), anti-cancer (Satyavarapu et al., 2020), immuno-modulatory (Gurjar and Pal, 2021), antihypertensive (An et al., 2020) anti-obesity (Kumar et al., 2020) and lipid lowering agent (Sandamali et al., 2020).These unique medicinal properties of MK leaves are accredited to their diverse metabolites. These include carbazole alkaloids, flavonoids, terpenoids, carbohydrates, phenolics, carotenoids, vitamins, essential oils, tannins, quinones, nicotinic acid, etc. (Balakrishnan et al., 2020).
A unique integration between environmental factors (temperature, photoperiod), soil type, and growth stage can alter the plant phenotype and its phytochemical composition, consequently changing the depot of bioactive compounds present in the MK leaves (Verma et al., 2012; Jimoh et al., 2019a). Also, variation in the germination response of MK to soil type suggests that mineral composition, water relation, and soil texture determines the plant growth, primarily the quantity and diversity of phytochemicals produced by the plant (Yang et al., 2018). Generally, for medicinal plants like MK, metabolic changes during the three growth stages (pre-flowering, flowering, and fruiting) and those caused by the environment can redirect the induction, attenuation, synthesis and transport of various bioactive phytochemicals (phenolics, tannins, antioxidant activity, essential oils, quinones, flavonoids etc.), consequently changing the complete phytochemical profile (Verma and Shukla, 2015). These variations can be viewed as a strategy imparting the plant an ability to adapt and thrive in the changing environments and for establishing a stable ecological relationship between the plant and other microorganisms (Metlen et al., 2009).
Multivariate data analysis (MVDA) is an important statistical approach to efficiently manage complex data (Granato et al., 2018). MVDA delineates the simultaneous analysis of one or more outcome variables (dependent variables) against two or more input variables (independent variables; Hidalgo and Goodman, 2013). It helps to determine the variable dimensions, which is the intricate relationships between several variables contemplated together. Pertinent usage of MVDA can appropriately assess the complex relationships between the phytochemical properties of the assessed sample. It also helps to decipher the similitudes and dissimilarities between multiple sample types, their biochemical characteristics, growth stages, region of cultivation or for gauging the samples in a two-/three-dimensional factor-plane, based on their distinct characteristics (Keithley et al., 2009). Techniques for pattern recognition like hierarchical cluster analysis (HCA) and principal component analysis (PCA) are primarily exploratory data analysis tools that are used for enhancing our comprehension of the data structure and can assess the distribution and association between phytochemicals and their place of cultivation (Zielinski et al., 2018).
Even though the phytochemistry of MK leaves have been investigated in the past, the information on this crop is very scanty and dispersed, especially for those popularly grown in the Indian subcontinent. There is however no report on the phytochemical composition of MK leaves at three different growth stages (pre-flowering, flowering, and fruiting) obtained from different elevations and the use of MVDA tools for deciphering the phytochemical diversity. Thus, the current investigation was performed to characterize the total phenolics, flavonoids, essential oils, quinones, tannins and antioxidant activity of 15 wild samples of MK leaves obtained from different locations of Himachal Pradesh at various growth stages using MVDA.
Materials and methods
Sample collection and preparation
MK fresh leaves were harvested at three different stages viz. pre-flowering, flowering and fruiting stage representing 15 different locations of Kangra and Mandi Districts of Himachal Pradesh. The locations include Pudva, Dhramn, Dhira, Mlan, Malghota, Tang, Nagri, Sagned, Dharampur, Lad Badhol, Neri, Madi, Drang, Kunnu, and Kotropi. The areas explored lies between (31°47′–32°50′ N and 76°25′–76°56′ E and 689–1,339 m altitudinal range). Standard practice and procedures for collecting the samples were followed. Mostly, the local plants were the target of the collection. Random and bulk sampling methods were followed as per the population/quantity of available germplasm material. Passport data, viz. location of the site, i.e., village, tehsil, district, latitude, longitude, and altitude were recorded at the different locations (Table 1; Figure 1). After harvesting, the leaf samples were cleaned and divided in two parts. One part was wrapped in aluminum foil, kept in airtight bags and then stored at −20°C. These samples were further used for essential oil yield quantification. The other part was kept for drying in an oven at 60°C–70°C till constant weight. These temperatures do not cause any significant loss of plant phytochemicals. About 10 g of oven dried samples were ground into fine powder followed by manual sieving by using a sieve size of 250 mm and were further utilized for estimation of total phenols, flavonoids, quinones, tannins, and antioxidant activity.
Total phenolic content
The phenolic content in MK leaf samples was estimated according to the Folin–Ciocalteu method as described by Julkunen-Tiitto (1985) using a tannic acid standard. Briefly, 200 mg of dried powder in triplicate were taken in test tubes containing 10 ml of 70% acetone and extracted at 37°C for 2 h using a shaking water bath. Subsequently, the tubes were cooled to room temperature and centrifuged at 1,000 g for 20 min. The resulting supernatant obtained was filtered, collected in fresh tubes and used as an extract for estimating total phenols. Hundred μl of sample extract in triplicates was then taken in different test tubes containing 900 μl distilled water. Concurrently blank was prepared in separate tubes by adding 1 ml distilled water. For standard preparation, tannic acid with different concentrations (0.005, 0.01, 0.015, 0.020, 0.025, 0.030, 0.035, 0.040, 0.045, and 0.05 mg) was added to a separate set of test tubes and volume was made to 1 ml with distilled water. Further, 0.5 ml of Folin–Ciocalteu reagent (1N) and 2.5 ml of 20 per cent Na2CO3 (w/v) was added to each of the three sets of test tubes (sample, blank and standards). The test tubes were wrapped in aluminum foil and kept at incubation for 40 min at room temperature. Folin–Ciocalteu reagent is a mixture of molybdates and tungstates and works by the oxidation–reduction reaction. The mixture of molybdates-heteropolyphosphotungsates is reduced by phenolic compounds forming a blue colored chromogen under basic condition by Na2CO3. The absorbance of the resulting solution (sample, standard) was measured at 725 nm against a blank by using a UV–visible spectrophotometer. The total phenolic content was computed using a standard curve made from tannic acid and the results are depicted as mg tannic acid equivalents (TAE) per gram dry weight basis (DWB) of MK leaves.
Total flavonoid content
The total flavonoids content were determined based on the procedure proposed by (Swain and Hillis, 1959) with slight modification. One-gram dried powdered MK leaf samples in triplicate were ground with 10 ml of methanol and the whole content was transferred to the test tubes. The tubes were then agitated constantly for 12 h in a shaking water bath maintained at 37°C. Thereafter, the content was centrifuged for 20 min at 3,500 g. The filtrate was collected in a beaker and dried with the help of a rotary evaporator. The dried extract obtained was re-dissolving in 1 ml of methanol. Briefly, 250 μl of methanol extract was taken in triplicate in test tubes and then 2,250 μl distilled water was added to it. Simultaneously, blank was prepared by adding 2,500 μl distilled water in separate tubes. For the preparation of standard, catechin, 5, 10, 15, 20, 25, 30, and 35 μg was added to a separate set of tubes and the volume was made to 2,500 μl with distilled water. After that 150 μl of 15%, NaNO3 solution was added to the reaction mixture (sample, blank and standards) and then incubated for 6 min at ambient temperature. After incubation, 150 μl of 10% AlCl3 was added and then 2,000 μl of 4% NaOH and 200 μl of distilled water were added. The solution was mixed well and the absorbance of the reaction mixture (sample, standard) was recorded at 510 nm against blank with the help of a UV–visible spectrophotometer. The total phenolic content was calculated with reference to the catechin standard curve and results are expressed in mg catechin equivalent (CE) per gram DWB of MK leaf powder.
In vitro antioxidant activity assays
This assay determined the potential of MK crude extracts to neutralize the free radical, 2,2-diphenyl-1-picrylhydrazyl (DPPH). DPPH is a stable free radical which has an unpaired valence electron at one atom of nitrogen bridge. The antioxidant activity was determined according to a method by Sharma and Bhat (2009) with slight modifications. Each extract in triplicates and catechin standard (20, 40, 60, 80, 100, 120, and 140 μl) was taken in a separate set of test tubes, subsequently, 3.0 ml of methanol and 1.0 ml of DPPH was added to it. A combination of 3.0 ml of methanol and 1.0 ml of DPPH in a separate test tube was used as a control. The three sets of a tube (sample, standard, and control) were kept in dark at 30°C for 30 min. The absorbance of each solution (sample, standard and control) was measures at 517 nm using methanol as blank. Finally, the percentage of DPPH free radical scavenging activity (per cent inhibition) of MK extracts was calculated according to Eq (1)
IC50 value (to measure the antioxidant amount needed to reduce the free radical DPPH concentration by 50%) for each solution was computed from the regression line generated from the graph of per cent DPPH inhibition against the concentration of each solution using Eq (2)
Essential oil yield
Essential oil yield was determined through a method by Yadegarinia et al. (2006) using the Clevenger apparatus. Eighty-gram fresh leaf of each sample in triplicate was placed in a 1,000 ml round bottom flask and subjected to hydro-distillation for 6 h. The oil collected was dried over a hot water bath and preserved at 4°C for further analysis. Essential oil yield (%) was determined using Eq (3):
Quinones
Quinones were estimated through a method by Mahadevan and Sridhar (1986). Briefly, 1.0 g sample in triplicate were ground with chilled 10 ml 0.1 M sodium phosphate buffer (pH 6.6) using pestle and mortar. The mixture was then centrifuged at 1,000 g for 30 min at 4°C. The supernatant was collected in fresh tubes and was used as an extract. Accurately 1.5 ml extract in triplicate was taken in a test tube. Then 3.0 ml phosphate buffer and 3.0 ml standard catechol were added to it. The mixture was vortexed gently and incubated in a water bath for 1 h at 60°C. Nearly 2.0 ml sample was withdrawn in a test tube in duplicates. Four ml of Trichloroacetic acid (without ascorbic acid) was added to one test tube and 4.0 ml of trichloroacetic acid (containing ascorbic acid) to another. The precipitate was filtered and the absorbance of each solution was measured at 400 nm against a blank. The results were expressed as mM/min/g tissue DWB.
Total tannin content
The total tannin content of MK samples was estimated by the Folin–Denis method (Pansera et al., 2003) with slight modifications. Briefly, 200 mg of dried powder in triplicate were taken in test tubes containing 10 ml of 70% acetone and extracted at 37°C for 2 h using a shaking water bath. Subsequently, the tubes were cooled to room temperature and centrifuged at 1,000 g for 20 min. The resulting supernatant obtained was filtered, collected in fresh tubes and used as an extract for estimating total tannin content. Briefly, 400 μl of extract in triplicate was taken in test tubes and 400 μl of Folin–Denis reagent was added to it. After 3 min, 400 μl of 8% Na2CO3 solution was added to these test tubes. The solution was mixed properly and then incubated for 1 h. The solution was centrifuged for 5 min at 1,000 g and the absorbance of each solution was taken at 725 nm against blank. For the preparation of standard, tannic acid (12.5, 25, 50, 100, 500, and 1,000 μg) was used and the results are expressed as mg tannic acid equivalents (TAE) per gram DWB of MK leaves.
Statistics
ArcGIS 10.1 platform (ESRI Inc., United States) was used for preparing the map for the study area and sample location sites. Unless stated otherwise all the analysis was performed as triplicates. Statistical significance was set at a 95% confidence level. Hierarchical cluster analysis (HCA) was done by Ward clustering algorithm using squared Euclidean distance by a trial version of IBM SPSS. The correlation between total phenolics, flavonoids, essential oils, quinones, tannins and antioxidant activities were assessed using Pearson’s correlation coefficient. The correlation analysis was performed using Jamovi version 1.2.27 at a 5% level of significance. The differences between various phytochemical attributes were statistically assessed by one-way analysis of variance (ANOVA). The differences among means were evaluated by post-hoc test (Duncan’s multiple range test) at (p < 0.05) level using a trial version of IBM SPSS. The data were subjected to principal component analysis (PCA) for assessing the phytochemical constituents in MK leaves that determine the observed variation and visualize the tentative relationships among these components using Jamovi version 1.2.27. To maximize the loading of a parameter in the component axis and for interpreting the pattern of a specific parameter, the orthogonal varimax rotation was applied to PCA. The heatmap for relative composition and Pearson’s correlation coefficient was constructed using a web interface (MetaboAnalyst; Chong et al., 2019).
Results and discussion
Multivariate data analysis
To enhance our understanding of the data organization and structure, its classification and distribution, PCA and HCA were used for visualizing the phytochemical composition of wild MK leaves collected from 15 different locations of Himachal Pradesh, India, at different growth stages, viz. pre-flowering, flowering, and fruiting stages.
Hierarchical clustering analysis
Hierarchical clustering analysis (HCA) is a statistical exploratory analysis technique that is used for differentiating sample groups or individual variables using similarity or density measurements. HCA was used to explore the sample organization into individual groups, and among groups, delineating their hierarchy. The results of the HCA are represented through a dendrogram, which shows the association between samples in form of a tree. In a typical dendrogram, the samples are joined through branches with corresponding similarity values. Similarities present in the traits of these samples are determined by an inverse function of their inter-distance. Various matrices are used for assessing the distance between two samples in a specific multivariate space, most common being Mahalanobis distance and Euclidean distance. In the present study, Ward’s linkage was used as the linkage criterion, which is based on the optimal values of a target function(Granato et al., 2015) and Euclidean distance was used as the adequate linkage benchmark and sample distance metric. Also, the samples were agglomeratively clustered, where each sample was initially regarded as a cluster and later, the cluster pairs were merged.
The HCA was performed at each stage of MK leaf development (viz. pre-flowering, flowering and fruiting). The leaf samples were grouped into four clusters for each of the development stages. The mean value of the phytochemical composition for each cluster for every growth stage is indicated in Table 2 and is depicted by Box-and-whisker plots in Supplementary Figures S1–S3. The hierarchical clusters along with the scaled relative phytochemical composition is depicted through heat maps in Figures 2A–F.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/963150/fpls-13-963150-HTML/image_m/fpls-13-963150-t002.jpg)
Table 2. Mean values of phytochemicals in MK leaves organized into different clusters at each growth stage.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/963150/fpls-13-963150-HTML/image_m/fpls-13-963150-g002.jpg)
Figure 2. The assessed samples were organized into four different clusters, depicting their hierarchy and presented as a dendrogram and scaled phytochemical attributes are represented through a heat map. (A) represents the hierarchical cluster integrated with a heat map representing the scaled phytochemical attributes for each sample at the pre-flowering stage. The data is scaled between +2 and −2; (B) represents the hierarchical cluster integrated with a heat map representing the mean scaled value of each cluster at the pre-flowering stage. The data is scaled between +1 and –1; (C) represents the hierarchical cluster integrated with a heat map representing the scaled phytochemical attributes for each sample at the flowering stage. The data is scaled between +1 and −1; (D) represents the hierarchical cluster integrated with a heat map representing the mean scaled value of each cluster at the flowering stage. The data is scaled between +1 and −1; (E) represents the hierarchical cluster integrated with a heat map representing the scaled phytochemical attributes for each sample at the fruiting stage. The data is scaled between +2 and −2; (F) represents the hierarchical cluster integrated with a heat map representing the mean scaled value of each cluster at the fruiting stage. The data is scaled between +1 and −1.
For the pre-flowering stage, cluster I was characterized by samples having lower phenol content (180.32 ± 16.56 mg TAE/g) and tannin (76.25 ± 0.66 mg TAE/g) content. Cluster II was characterized by samples with higher content of essential oils (0.82% ± 0.061%), antioxidant activity (6.3 ± 1.20 μg/ml) and flavonoids (25.62 ± 6.58 mg CE/g). Cluster III is indicated by samples with higher content of quinones (2.58 ± 0.35 mM/min/g tissue). Cluster IV is characterized by germplasm with higher content of phenol content (294.5 ± 6.45 mg TAE/g), tannins (137.5 ± 4.50 mg TAE/g) and lower content of essential oils (0.73 ± 0.11%). For the flowering stage, cluster I was represented by samples having lower antioxidant activity. Cluster II was characterized by samples having higher content of phenol content (450.01 ± 9.04 mg TAE/g), flavonoids (58.17 ± 1.04 mg CE/g), antioxidant activity (13.6 ± 0.75 μg/ml), essential oil (1% ± 0%), tannins (207 ± 2.14 mg TAE/g), and quinones (4.95 ± 0.89 mM/min/g tissue). Cluster III was characterized by samples having lower content of phenols (287.86 ± 12.93 mg TAE/g), flavonoids (30.21 ± 2.94 mg CE/g), antioxidant activity (9.09 ± 0.54 μg/ml), essential oil (0.9% ± 0.05%), tannins (119.9 ± 2.43 mg TAE/g), and quinones (3.65 ± 0.56 mM/min/g tissue). For fruiting stage, cluster I was characterized by samples having a higher content of flavonoids (29.32 ± 2.43 mg CE/g), essential oils (0.67 ± 0.03%), and tannins (142.33 ± 14.57 mg TAE/g). The cluster II was displayed by samples having a lower antioxidant activity (4.46 ± 0.77 μg/ml). Cluster III was depicted by germplasms having lower content of phenols (229.75 ± 14.33 mg TAE/g), tannins (54.19 ± 10.76 mg TAE/g), essential oil (0.56% ± 0.02%), and quinones (1.41 ± 0.22 mM/min/g tissue). Cluster IV was indicated by germplasms with higher content of quinones (1.68 ± 0.27 mM/min/g tissue), antioxidant activity (5.39 ± 0.53 μg/ml) and lower content of flavonoids (24.23 ± 6.0 mg CE/g).
Principal component analysis
The phytochemical composition of MK leaves for pre-flowering, flowering and fruiting stages were scaled and subjected to factor analysis using principal component analysis (PCA). This was done to distinguish the traits that determine the observed variability in the composition and visualize the potential associations among these traits. PCA algorithms investigate the direction having the maximum variation within a multidimensional data space. This is based on the postulation that the maximum variability (depicted by a high value of variance) equates to more information. To discern the accuracy of this direction, the data matrix was mean-centered column-wise, ascertaining that the axes rotate about the data centroid. The rotating axis explaining the highest variance delineates the first (lowest-order) principal component (PC1). The second axis resides in the direction with the maximum variance noncorrelated (orthogonal) among all directions with respect to PC1. Subsequently, the second PC (PC2) explains all possible information that is not explained by PC1 and so on for other PCs. The PCs substantiate the same information obtained from the assessed parameters but exhibits the additional advantage of being uncorrelated reciprocally, circumventing any copious information between the parameters (Granato et al., 2018). PCs can be regarded as a linear amalgamation of original variables, with each multiplied with a weight coefficient, known as loading. Through a geometrical perspective, loadings portray the cosine values of angle between PCs and the original variables. These loading values range between −1 and + 1, insinuating their role in defining a specific PC. Greater is the cosine value, lesser is the rotational angle between PC and the variable, nearer are the two directions, higher is the contribution of that variable to the given PC. PCA has a distinct ability to display the lengthy and complex information present in the data through bi and tridimensional plots (biplots, score plots, and loading plots). Each PC accounts for a part of the variation in the data is determined through Eigenvalue (Dobriban, 2017). Every Eigenvalue accords to an Eigenvector, determining the variance among the associated PCs. Factors with Eigenvalues > 1.0 were seen in agreement with Kaiser’s criteria (Kaiser, 1960).
Before subjecting the data to component analysis, Kaiser-Meyer-Olkin (KMO) Measure of Sampling Adequacy (MSA) and Bartlett’s test for sphericity was used as the main assumption checks. Bartlett’s test juxtaposes the identity matrix to the observed correlation matrix (Tobias and Carlson, 1969). It helps to identify the possible redundancy among variables that can be reiterated by a small number of factors. This test validates whether or not the variables in the data are correlated in such a way that they can be potentially outlined by a smaller set of factors. The null hypothesis for this test is that the variables are uncorrelated or orthogonal. The alternate hypothesis states that the variables are sufficiently correlated to where the correlation matrix significantly diverges from the identity matrix. Thus if Bartlett’s test is significant at (p < 0.5), the correlation matrix is significantly divergent from the null hypothesis, indicating that the data reduction through PCA can meaningfully compress the data. The values of Bartlett’s test for our data were p < 0.001 for phenols, flavonoids, antioxidant activity, essential oils, tannins, and quinones (Supplementary Table S1).
The KMO is another statistical test that determines the suitability of data for factor analysis through PCA. It assesses the variance proportion of the variables that could result from the underlying factors. Lesser is the proportion of variation, more is the suitability of the data for PCA (Lestari, 2021). Values ranging between 0.9 and 1.0 are contemplated as excellent for subjecting the data to PCA. The values between 0.7 and 0.9 are regarded as good and those between 0.5 and 0.7 are adequate. If the values are lower than 0.5, the PCA results will be unable to provide any meaningful interpretation. The KMO values in the present study ranged from 0.665 to 0.850 (Supplementary Table S2), showing the high appropriateness of the data for PCA.
Following the determination of the representative PCs on the basis of explained variance and sample grouping/differentiation, the loadings were assessed to elucidate the fundamental association between the data structure. The loading can be contemplated as linear combination coefficients of the initial variables from which the PCs are synthesized. These factor loadings (FLs) are the variable coordinates divided by the square root of the eigenvalue associated with the specific component. FLs having positive values indicate that the factor will lie higher on the positive axis of that PC. In the present study, quinones (FL1: 0.911), essential oils (FL: 0.882), antioxidant activity (FL1: 0.816), tannins (FL2: 0.910), phenols (FL2: 0.809), flavonoids (FL3: 0.902; Supplementary Table S3).
PCA revealed that the first two components explained 99.1% variability in the data, in terms of phenols, flavonoids, antioxidant activity, essential oils, tannins, and quinones. PC1 described 93.9% of the total variation, in which the dominant parameters include quinones, essential oils and antioxidant activity. PC2 explained 5.2% of the total variation in which the chief parameters include phenols and tannins. These determining factors showed high values of FLs and lower uniqueness. Based on the calculated FL values, the 2D and 3D PCA score plots showing the distribution of phytochemicals in ML leaves at three different growth stages is indicated in Figure 3.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/963150/fpls-13-963150-HTML/image_m/fpls-13-963150-g003.jpg)
Figure 3. Represents the PCA (Principal component analysis); (A) two-dimensional and (B) three-dimensional PCA component loading plot for all three developmental stages (pre-flowering, flowering, and fruiting). The samples lying in each coordinate are clearly indicated; (C) indicated the three-dimensional PCA score plot; (D) indicates two-dimensional component loading plot for each of the assessed phytochemical parameters.
Phytochemical composition of MK leaves and its relation with growth stage
Plants synthesize and accumulate a diversity of high and low molecular weight secondary metabolites which exhibit potent bioactivities (Quideau et al., 2011; Kasote et al., 2015). These natural bioactive metabolites belong to distinct chemical groups like flavonoids, phenols, quinones, tannins, diterpenes, isoflavones, anthocyanins (Jimoh et al., 2019b). These phytochemicals display unique overlapping roles as hydrogen donors, singlet oxygen quenchers, reducing agents, plant defense agents against biotic and abiotic stresses, radiation absorbers, allelopathy, pollinator attractors, transition metal chelators (Harborne and Williams, 2000; Demidchik, 2015; do Nascimento et al., 2018). The concentration of these show variation in different plant parts, growth phases, planting seasons (Lagnika et al., 2016). A unique combination of externally provided mineral resources from the soil (place of growth) and developmental stage delineate the carbon-nutrient balance of the plant. This determines the biosynthesis and translocation of these metabolites in the plants (Mahmood et al., 2012). It is therefore critical to assess the phytochemical composition of MK leaves at different growth stages to accurately ascertain the pharmacological potential of the plant. Our results indicated that the growth stages significantly influence the phytochemistry of the plant. Our data indicated that all evaluated phytochemical parameters in MK leaves significantly increased during the growth transition from pre-flowering to flowering stage, followed by a gradual decrease during the fruiting stage (Tanko et al., 2005). This could possibly result from the variation in expression of phenylpropanoid pathway mediating enzymes during different growth stages (Kariñho-Betancourt et al., 2019) and differences in root absorption capacity, soil nutrient availability, climatic conditions, photosynthesis rate, etc. (Briat et al., 2020; Tomar et al., 2021). The mean phytochemical composition of MK leaves is represented through Box and Whisker Plots at all the three developmental stages in Figure 4, by bar graphs in Supplementary Figure S4 and through comprehensive heat maps in Figure 5. Also the phytochemical composition is indicated in Table 3.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/963150/fpls-13-963150-HTML/image_m/fpls-13-963150-g004.jpg)
Figure 4. Represents the phytochemical profile of MK (Murraya koenigii (L.) Spreng.) leaves through Box and Whisker Plots at the three developmental stages (pre-flowering, flowering, and fruiting) for (A) phenols; (B) flavonoids; (C) antioxidant activity; (D) tannins; (E) essential oils; and (F) quinones.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/963150/fpls-13-963150-HTML/image_m/fpls-13-963150-g005.jpg)
Figure 5. Represents the phytochemical composition of MK (Murraya koenigii (L.) Spreng.) leaves through a heat map. (A1) Indicates the scaled mean phytochemical profile of MK leaves at three development stages. The data is scaled between +1 and –1. (A2) Represents the scaled photochemical profile of MK leaves for each sample at each of the three development stages. The data is scaled between +3 and −3; (B) Indicates the Pearson’s r correlation values through a heat map. The Pearson’s r (+1 to −1) is indicated through different color intensities.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/963150/fpls-13-963150-HTML/image_m/fpls-13-963150-t003.jpg)
Table 3. Values of phytochemicals in MK leaves organized into different clusters at each growth stage.
For the pre-flowering stage, the total phenols in MK leaves ranged from 170.09 ± 4.59 mg TAE/g in samples collected from Drang to 303.57 ± 7.94 mg TAE/g in samples obtained from Dharampur. For the flowering stage, the total phenols ranged from 266.48 ± 7.49 mg TAE/g in Sagned samples to 450.01 ± 11.78 mg TAE/g in Dharampur samples. In the fruiting stage, the content ranged from 212.72 ± 5.37 mg TAE/g in Mlan samples to 363.85 ± 9.79 mg TAE/gm in Madi samples. The results corroborated with those of (Sivakumar and Meera, 2013) who indicated the total phenol content of 111.6 ± 3.85 to 532.8 ± 2.81 mg GAE/g, similarly (Kaur et al., 2016) indicated a content of 119.6 mg GAE/g in MK leaves. Also for the pre-flowering stage, the total flavonoid content MK leaves ranged from 15.17 ± 0.36 mg CE/g in Malghota samples to 33.40 ± 0.81 mg CE/g in Dhramn samples. For the flowering stage, the content ranged from 25.16 ± 0.67 mg CE/g in Madi samples to 58.17 ± 1.52 mg CE/g in Dharampur samples. For the fruiting stage, the content ranged from 17.54 ± 0.42 mg CE/g in Lad badhol samples to 37.34 ± 0.97 mg CE/g in Dharampur samples. (Sindhu and Arora, 2012) also indicated a similar flavonoid content of 43.58 ± 1.89 mg CE/g.
This increase in the phenolic and flavonoid content from pre-flowering to flowering stage and their subsequent decrease at fruiting stage was also observed by (Kasera et al., 2018) in medicinal plants namely Withania coagulans, Dipcadi erythraeum, Corbichonia Decumbens and Arisaema tortuosum. Vlaisavljević et al. (2017) stated that variations in phenolics and flavonoids during transitions from one growth phase to another is fundamentally related to the changes in the biosynthetic pathways at the molecular level. These changes alter the synthesis, allocation and disintegration of these metabolites during developmental stages (Pretti et al., 2018). Further, the phenolic content can decrease as a result of the breakdown of secondary metabolites into simple intermediate metabolites like lignin (Shao et al., 2014). These findings relate to an earlier study that indicated that the phenolic and flavonoid contents were highest in the flowering stage and substantially decreased to half at the fruiting stages (Pinhatti et al., 2010; Ben Farhat et al., 2015; Oszmiański et al., 2018; Farhadi et al., 2020). For the pre-flowering stage, the tannins in MK leaves ranged from 75.75 ± 1.69 mg TAE/g in Dhira samples to 143 ± 3.74 in Dharampur samples. For the flowering stage, the tannin content ranged from 116 ± 3.26 mg TAE/g in Sagned samples to 207 ± 5.42 mg TAE/g in Dharampur samples. For the fruiting stage, the tannin content ranged from 47 ± 1.18 mg TAE/g in Mlan samples to 156 ± 4.05 mg TAE/g in Nagri samples. A similar content of 206.05 ± 7.50 mg TAE/g was reported by (Uraku and Nwankwo, 2015) and 0.86 ± 0.02 g/100 g by (Igara et al., 2016).
High flavonoid and tannin content in flowering stages was also confirmed in Ziziphora clinopodioides by (Ding et al., 2014), Nigella sativa by (Zribi et al., 2014). An increase in flavonoid concentration at the flowering stage could result could be attributable to the fact that flavonoids determine flower and fruit coloration, flower aroma, pigmentation, ultra-violet protection (Nishihara and Nakatsuka, 2011). More flavonoids and tannins are likely to be synthesized by the plant with the onset of flowering. An increase in biosynthesis of flavonoids and tannins during flowering could be a defense mechanism against a wide variety of pests that could attack the plants and could serve as chemical signals for attracting potential pollinators, oviposition stimulants, phytoalexins, feeding attractants, deterrents (Iwashina, 2003).
Phenylpropanoids like phenols, flavonoids, tannins are synthesized from the shikimate pathway. The first committed step for protein and phenolic synthesis is catalyzed by phenylalanine ammonia-lyase, using phenylalanine (an aromatic amino acid) as a precursor (Yadav et al., 2020). The protein competition model states that the competition between proteins and phenolic synthesis for the common limiting phenylalanine, leads to a process-level tradeoff between the phenolic versus phenolic synthesis rates, establishing an inverse relation between phenolic and protein allocation (Jones and Hartley, 1999). Thus, when protein synthesis rates are high, the phenolic synthesis must be low and vice versa. Protein synthesis is increased during the pre-flowering and fruiting stage of plant development, while it is comparably less during the flowering stage (Frenkel et al., 1968). Thus a logical explanation could indicate that the high rates of protein synthesis during pre-flowering and fruiting stages could channelize the phenylalanine and carbon sources toward proteins synthesized, compromising the phenylpropanoid biosynthesis. Similarly lower protein synthesis during the flowering stages could enable the carbon flow toward the biosynthesis of phenols and flavonoids. Studies also indicate that these compounds have important roles in adaptation and providing defense against pathogens, herbivores and other environmental stresses during the critical flowering stages (Akula and Ravishankar, 2011).
For the pre-flowering stage, the essential oils in MK leaves ranged from 0.64 ± 0.01% in samples obtained from Mlan to 0.89% ± 0.02% in Dhramn samples. For the flowering stage, the content ranged from 0.85 ± 0.02% in Malghota samples to 1 ± 0.02% Dharampur samples. For the fruiting stage, the content ranged from 0.54 ± 0.01% in Lad badhol samples to 0.7% ± 0.01% in Dharampur samples. A similar essential oil yield of 0.24% was reported by (Rana et al., 2004) and 0.12%–0.18% by (Verma et al., 2012). Castelo et al. (2012) stated that the essential oil production was enhanced during the flowering period to attract pollinators such as bees and other insects.
Quinones are oxidized derivatives of aromatic compounds readily made from reactive aromatic compounds with electron-donating substituents such as phenols and catechols. Quinones are a group of compounds that are majorly found in the cell membranes of living organisms. They have a hydrophilic head and non-polar isoprenoid side chain, giving these molecules lipid-soluble properties. Quinones primarily function as proton and electron carriers in respiratory and photosynthetic electron transport chains, with additional antioxidant functions (Nowicka and Kruk, 2010). For the pre-flowering stage, the quinone in MK leaves ranged from 2.05 ± 0.05 mM/min/g tissue in Madi samples to 2.97 ± 0.07 mM/min/g tissue in Tang samples. For the flowering stage, the quinone content ranged from 3.07 ± 0.07 mM/min/g tissue in Kunnu samples to 4.95 ± 0.13 mM/min/g tissue in Dharampur samples. For the fruiting stage, the content ranged from 1.02 ± 0.02 mM/min/g tissue in Pudva samples to 1.96 ± 0.04 mM/min/g tissue in Neri samples.
For the pre-flowering stage, the antioxidant activity in MK leaves ranged from 4.01 ± 0.09 μg/ml in Kunnu samples to 7.42 ± 0.17 μg/ml in Malghota samples. For the flowering stage, the content ranged from 8.08 ± 0.19 μg/ml in Dhramn samples to 13.60 ± 0.35 μg/ml in Dharampur samples. For the fruiting stage, the content ranged from 3.11 ± 0.06 μg/ml in Dhira samples to 6.37 ± 0.15 μg/ml in Lad badhol. The high amount of antioxidant activities during flowering could be attributable to the higher content of phenols, flavonoids, quinones, tannins and essential oils during these growth stages. Studies indicate that phenols, flavonoids, tannins, and quinones have a distinct property to scavenge reactive oxygen species (ROS), which encompass nonradical oxygen species (HOCl, 1O2, H2O2, NO•, HO•, O2–•), radical oxygen species and oxidatively generated free radicals like ROO• and RO• like those created through oligonucleic acids (RNA and DNA), proteins, low-density lipoproteins (LDLs; Quideau et al., 2011). These compounds also display their antioxidant activities by chelating metal ions like iron(III)/copper(II) and iron(II)/copper(I) which mediate the conversion of O2–• and H2O2 into highly reactive HO• via Haber-Weiss/Fenton-type reactions (Andjelković et al., 2006). They also inhibit the activities of enzymes generating O2–• like protein kinase C and xanthine oxidase (Pietta, 2000). These compounds exhibit their antioxidant activities by two main mechanisms. First is their ability to donate a hydrogen atom to free radicals, making the phenolic compounds as free radicals. The second mechanism is through single-electron transfer to a free radical forming a stable radical cation.
Correlation analysis
The correlation coefficient was studied among phenols, flavonoids, antioxidant activity, essential oils, tannins, and quinones. Pearson’s correlation test was used for determining significant correlations among these parameters. Pearson’s r > 0 at p < 0.05 indicate significant positive correlation, while r < 0 at p < 0.05 indicate significant negative correlations. The values of Pearson’s r are indicated through a heat map in Figure 5, and their precise values are indicated in Supplementary Table S4. Antioxidant activity showed significant positive correlation with phenol content (r = 0.536; p < 0.001), essential oils (r = 0.707; p < 0.001), tannins (r = 0.564; p < 0.001), quinones (r = 0.820; p < 0.001), and flavonoids (r = 0.540; p < 0.001). As already mentioned, these highly significant positive correlations arise due to the ability of these compounds to reduce cellular oxygen concentration, intercept singlet oxygen (1O2) production, block the reaction chain initiation by scavenging hydroxyl radicals, chelate metal ions, terminate free radicals. These antioxidant activities rely on the unique structure and hydroxyl group arrangement in these compounds (Shahidi and Ambigaipalan, 2015). It is now clearly indicated by various studies that these phytochemicals have lower reduction potential than that of a free radical (oxidized species), thus they can instantaneously transfer their hydrogen atom to the free radical (Neudörffer et al., 2006; Quideau et al., 2011). Phenols showed significant positive correlations with flavonoids (r = 0.734; p < 0.001), tannins (r = 0.857; p < 0.001), quinones (r = 0.535; p < 0.001). This could be because all these compounds originate from the same polyketide and/or shikimate-derived phenylpropanoid pathway(s). Also, essential oils were positively correlated with phenols (r = 0.513; p < 0.001) and flavonoids (r = 0.470; p = 0.001) since the pathways of terpene (primary component of essential oils) biosynthesis (mevalonic acid and methylerythritol phosphate pathway) and phenol biosynthesis (malonic acid pathways and shikimic acid pathway) are interconnected through pyruvate and acetyl CoA (Mijts and Schmidt-Dannert, 2003; Mahajan et al., 2020). Thus, the increase in the synthesis of phenols, upregulates the terpene biosynthesis, subsequently increasing the essential oil content. Essential oils also indicated a significant positive correlation with antioxidant activity (r = 0.707; p < 0.001). This could be because of the strong peroxide decomposition and free radicals scavenging activities of terpenes present in the MK leaf essential oils (Bhavaniramya et al., 2019; Liang et al., 2021).
Conclusion
The present study indicated the evidence that geographical locations and plant developmental stages primarily determine the total phenolics, flavonoids, essential oils, quinones, tannins and antioxidant activity in MK leaves. All evaluated phytochemical parameters significantly increased during the growth transition from pre-flowering to the flowering stage, followed by their gradual decrease during the fruiting stage. MVDA tools like PCA and HCA were used for visualizing the photochemical composition of MK leaves at different growth stages. HCA classified the samples into four different clusters based on their phytochemistry at three different developmental stages. PCA indicated that quinones, essential oils and antioxidant activity primarily determined the variability in the data. The study suggests that for utilizing the maximum phytochemical constituents, biological potency and industrial utility, MK leaves should be harvested at the flowering time. These phytochemicals can be used as raw materials in pharmaceuticals, cosmetics, insecticides, dyes, fragrances, nutraceuticals, etc. This study also establishes MVDA as an effective statistical tool to meaningfully decipher the phytochemical diversity in MK leaves. The information can significantly contribute toward the use of MK leaves aromatic and phytopharmaceutical industries.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author.
Author contributions
RV and NS designed the work. RV collected and analyzed the samples. RV and AR executed parts of the experiments in the laboratory. MT, RB, and DD statistically analyzed the research data. RV and MT drafted and revised the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the Department of Chemistry and Biochemistry, Chaudhary Sarwan Kumar Himachal Pradesh Agriculture University, Palampur, HP (India).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2022.963150/full#supplementary-material
References
Abeysinghe, D. T., Kumara, K. A. H., Kaushalya, K. A. D., Chandrika, U. G., and Alwis, D. D. D. H. (2021). Phytochemical screening, total polyphenol, flavonoid content, in vitro antioxidant and antibacterial activities of Sri Lankan varieties of Murraya koenigii and Micromelum minutum leaves. Heliyon 7:e07449. doi: 10.1016/j.heliyon.2021.e07449
Adebajo, A. C., Olayiwola, G., Eugen Verspohl, J., Iwalewa, E. O., Omisore, N. O. A., Bergenthal, D., et al. (2005). Evaluation of the Ethnomedical claims of Murraya koenigii. Pharm. Biol. 42, 610–620. doi: 10.1080/13880200490902518
Akula, R., and Ravishankar, G. A. (2011). Influence of abiotic stress signals on secondary metabolites in plants. Plant Signal. Behav. 6, 1720–1731. doi: 10.4161/psb.6.11.17613
An, Y., Guo, Y., Wu, H., Huang, Y., Xie, D., and Qin, J. (2020). Conceivable protective role of Murraya koenigii leaf extract loaded poly(D, L-lactic-co-glycolic acid)-gold-nanoparticles on the gestational diabetes mellitus of rats induced by Streptozotocin. Sci. Adv. Mater. 12, 87–92. doi: 10.1166/sam.2020.3569
Andjelković, M., Van Camp, J., De Meulenaer, B., Depaemelaere, G., Socaciu, C., Verloo, M., et al. (2006). Iron-chelation properties of phenolic acids bearing catechol and galloyl groups. Food Chem. 98, 23–31. doi: 10.1016/j.foodchem.2005.05.044
Balakrishnan, R., Vijayraja, D., Jo, S.-H., Ganesan, P., Su-Kim, I., and Choi, D.-K. (2020). Medicinal profile, phytochemistry, and pharmacological activities of Murraya koenigii and its primary bioactive compounds. Antioxidants 9:101. doi: 10.3390/antiox9020101
Ben Farhat, M., Chaouch-Hamada, R., Sotomayor, J. A., Landoulsi, A., and Jordán, M. J. (2015). Antioxidant properties and evaluation of phytochemical composition of Salvia verbenaca L. extracts at different developmental stages. Plant Foods Hum. Nutr. 70, 15–20. doi: 10.1007/s11130-015-0466-9
Bhandari, P. R. (2012). Curry leaf (Murraya koenigii) or cure leaf: review of its curative properties. J. Med. Nutr. Nutraceut. 1:92. doi: 10.4103/2278-019X.101295
Bhatt, S., Singh, B., and Gupta, M. (2020). Antioxidant and prebiotic potential of Murraya koenigii and Brassica oleracea var. botrytis leaves as food ingredient. J. Agric. Food Res. 2:100069. doi: 10.1016/j.jafr.2020.100069
Bhavaniramya, S., Vishnupriya, S., Al-Aboody, M. S., Vijayakumar, R., and Baskaran, D. (2019). Role of essential oils in food safety: antimicrobial and antioxidant applications. Grain Oil Sci. Technol. 2, 49–55. doi: 10.1016/j.gaost.2019.03.001
Binici, H., Şat, İ., and Aoudeh, E. (2021). The effect of different drying methods on nutritional composition and antioxidantactivity of purslane (Portulaca oleracea). Turk. J. Agric. For. 45, 680–689. doi: 10.3906/tar-2012-60
Briat, J.-F., Gojon, A., Plassard, C., Rouached, H., and Lemaire, G. (2020). Reappraisal of the central role of soil nutrient availability in nutrient management in light of recent advances in plant nutrition at crop and molecular levels. Eur. J. Agron. 116:126069. doi: 10.1016/j.eja.2020.126069
Castelo, A. V. M., Menezzi, C. H. S., and Resck, I. S. (2012). Seasonal variation in the yield and the chemical composition of essential oils from two Brazilian native arbustive species. J. Appl. Sci. 12, 753–760. doi: 10.3923/jas.2012.753.760
Chandra, S., and Nagani, K. (2009). In vitro models for antioxidant activity evaluation and some medicinal plants possessing antioxidant properties: an overview. AJMR 3, 981–996. doi: 10.5897/AJMR.9000401
Chong, J., Yamamoto, M., and Xia, J. (2019). MetaboAnalystR 2.0: from raw spectra to biological insights. Metabolites 9:57. doi: 10.3390/metabo9030057
Demidchik, V. (2015). Mechanisms of oxidative stress in plants: from classical chemistry to cell biology. Environ. Exp. Bot. 109, 212–228. doi: 10.1016/j.envexpbot.2014.06.021
Desai, S. N., Patel, D. K., Devkar, R. V., Patel, P. V., and Ramachandran, A. V. (2012). Hepatoprotective potential of polyphenol rich extract of Murraya koenigii L.: an in vivo study. Food Chem. Toxicol. 50, 310–314. doi: 10.1016/j.fct.2011.10.063
Ding, W., Yang, T., Liu, F., and Tian, S. (2014). Effect of different growth stages of Ziziphora clinopodioides lam. On its chemical composition. Pharmacogn. Mag. 10, S1–S5. doi: 10.4103/0973-1296.127329
do Nascimento, K. S., Gasparotto Sattler, J. A., Lauer Macedo, L. F., Serna González, C. V., Pereira de Melo, I. L., da Silva Araújo, E., et al. (2018). Phenolic compounds, antioxidant capacity and physicochemical properties of Brazilian Apis mellifera honeys. LWT 91, 85–94. doi: 10.1016/j.lwt.2018.01.016
Dobriban, E. (2017). Sharp detection in PCA under correlations: all eigenvalues matter. Ann. Stat. 45, 1810–1833. doi: 10.1214/16-AOS1514
Farhadi, N., Babaei, K., Farsaraei, S., Moghaddam, M., and Ghasemi Pirbalouti, A. (2020). Changes in essential oil compositions, total phenol, flavonoids and antioxidant capacity of Achillea millefolium at different growth stages. Ind. Crop Prod. 152:112570. doi: 10.1016/j.indcrop.2020.112570
Frenkel, C., Klein, I., and Dilley, D. R. (1968). Protein synthesis in relation to ripening of pome fruits 1. Plant Physiol. 43, 1146–1153. doi: 10.1104/pp.43.7.1146
Granato, D., Karnopp, A. R., and Ruth, S. M. (2015). Characterization and comparison of phenolic composition, antioxidant capacity and instrumental taste profile of juices from different botanical origins. J. Sci. Food Agric. 95, 1997–2006. doi: 10.1002/jsfa.6910
Granato, D., Santos, J. S., Escher, G. B., Ferreira, B. L., and Maggio, R. M. (2018). Use of principal component analysis (PCA) and hierarchical cluster analysis (HCA) for multivariate association between bioactive compounds and functional properties in foods: a critical perspective. Trends Food Sci. Technol. 72, 83–90. doi: 10.1016/j.tifs.2017.12.006
Grygorieva, O., Klymenko, S., Kuklina, A., Vinogradova, Y., Vergun, O., Sedlackova, V., et al. (2021). Evaluation of Lonicera caerulea L. genotypes based on morphological characteristics offruits germplasm collection. Turk. J. Agric. For. 45, 850–860. doi: 10.3906/tar-2002-14
Gul, M. Z., Attuluri, V., Qureshi, I. A., and Ghazi, I. A. (2012). Antioxidant and α-glucosidase inhibitory activities of Murraya koenigii leaf extracts. Pharm. J. 4, 65–72. doi: 10.5530/pj.2012.32.12
Gurjar, V. K., and Pal, D. (2021). “Natural compounds extracted from medicinal plants and their Immunomodulatory activities,” in Bioactive Natural Products for Pharmaceutical Applications Advanced Structured Materials. eds. D. Pal and A. K. Nayak (Cham: Springer International Publishing), 197–261.
Harborne, J. B., and Williams, C. A. (2000). Advances in flavonoid research since 1992. Phytochemistry 55, 481–504. doi: 10.1016/S0031-9422(00)00235-1
Hidalgo, B., and Goodman, M. (2013). Multivariate or multivariable regression? Am. J. Public Health 103, 39–40. doi: 10.2105/AJPH.2012.300897
Igara, C., Omoboyowa, D., Ahuchaogu, A., Orji, N., and Ndukwe, M. (2016). Phytochemical and nutritional profile of Murraya Koenigii (Linn) Spreng leaf. J. Pharmacogn. Phytochem. 5, 7–9.
Iwashina, T. (2003). Flavonoid function and activity to plants and other organisms. Biol. Sci. Space 17, 24–44. doi: 10.2187/bss.17.24
Jimoh, M. O., Afolayan, A. J., and Lewu, F. B. (2019a). Germination response of Amaranthus caudatus L. to soil types and environmental conditions. Thaiszia-J. Bot. 29, 85–100. doi: 10.33542/TJB2019-1-07
Jimoh, M. O., Afolayan, A. J., and Lewu, F. B. (2019b). Therapeutic uses of Amaranthus caudatus L. Trop. Biomed. 36, 1038–1053.
Jones, C. G., and Hartley, S. E. (1999). A protein competition model of phenolic allocation. Oikos 86, 27–44. doi: 10.2307/3546567
Julkunen-Tiitto, R. (1985). Phenolic constituents in the leaves of northern willows: methods for the analysis of certain phenolics. J. Agric. Food Chem. 33, 213–217. doi: 10.1021/jf00062a013
Kaiser, H. F. (1960). The application of electronic computers to factor analysis. Educ. Psychol. Meas. 20, 141–151. doi: 10.1177/001316446002000116
Kariñho-Betancourt, E., Hernández-Soto, P., Rendón-Anaya, M., Calderón-Cortés, N., and Oyama, K. (2019). Differential expression of genes associated with phenolic compounds in galls of Quercus castanea induced by Amphibolips michoacaensis. J. Plant Interact. 14, 177–186. doi: 10.1080/17429145.2019.1603404
Kasera, P. K., Soni, A., Daiya, K., and Sagar, A. (2018). Evaluation of phytochemical variations in medicinal plants of Rajasthan. Int. J. Phytocos. Nat. Ingred. 5:12. doi: 10.15171/ijpni.2018.12
Kasote, D. M., Katyare, S. S., Hegde, M. V., and Bae, H. (2015). Significance of antioxidant potential of plants and its relevance to therapeutic applications. Int. J. Biol. Sci. 11, 982–991. doi: 10.7150/ijbs.12096
Kaur, D., Kaur, A., and Arora, S. (2016). Delineation of attenuation of oxidative stress and mutagenic stress by Murraya exotica L. leaves. Springerplus 5:1037. doi: 10.1186/s40064-016-2709-0
Keithley, R. B., Mark Wightman, R., and Heien, M. L. (2009). Multivariate concentration determination using principal component regression with residual analysis. TrAC Trends Anal. Chem. 28, 1127–1136. doi: 10.1016/j.trac.2009.07.002
Kesari, A. N., Gupta, R. K., and Watal, G. (2005). Hypoglycemic effects of Murraya koenigii on normal and alloxan-diabetic rabbits. J. Ethnopharmacol. 97, 247–251. doi: 10.1016/j.jep.2004.11.006
Kumar, P., Singh, S., and Ahmad, M. I. (2020). Synergistic effect of Cinnamomum zeylanicum and Murraya koenigii formulation for antiobesity and hypolipidemic activity on wistar albino rats. Adv. Tradit. Med. 21, 553–563. doi: 10.1007/s13596-020-00460-8
Lagnika, L., Amoussa, A. M. O., Adjileye, R. A. A., Laleye, A., and Sanni, A. (2016). Antimicrobial, antioxidant, toxicity and phytochemical assessment of extracts from Acmella uliginosa, a leafy-vegetable consumed in Bénin, West Africa. BMC Complement. Altern. Med. 16:34. doi: 10.1186/s12906-016-1014-3
Lestari, M. W. (2021). Yield determinants of Some sweet potato (Ipomoea batatas L) clones using principal component analysis. IOP Conf. Ser. Earth Environ. Sci. 828:012006. doi: 10.1088/1755-1315/828/1/012006
Liang, Z., Zhi, H., Fang, Z., and Zhang, P. (2021). Genetic engineering of yeast, filamentous fungi and bacteria for terpene production and applications in food industry. Food Res. Int. 147:110487. doi: 10.1016/j.foodres.2021.110487
Mahadevan, A., and Sridhar, R. (1986). Methods in Physiological Plant Pathology. Madras, India: Sivakami Publications.
Mahajan, M., Kuiry, R., and Pal, P. K. (2020). Understanding the consequence of environmental stress for accumulation of secondary metabolites in medicinal and aromatic plants. J. Appl. Res. Med. Aromat. Plants 18:100255. doi: 10.1016/j.jarmap.2020.100255
Mahipal, P., and Pawar, R. S. (2017). Nephroprotective effect of Murraya koenigii on cyclophosphamide induced nephrotoxicity in rats. Asian Pac. J. Trop. Med. 10, 808–812. doi: 10.1016/j.apjtm.2017.08.005
Mahmood, T., Anwar, F., Abbas, M., and Saari, N. (2012). Effect of maturity on Phenolics (phenolic acids and flavonoids) profile of strawberry cultivars and mulberry species from Pakistan. Int. J. Mol. Sci. 13, 4591–4607. doi: 10.3390/ijms13044591
Metlen, K. L., Aschehoug, E. T., and Callaway, R. M. (2009). Plant behavioural ecology: dynamic plasticity in secondary metabolites. Plant Cell Environ. 32, 641–653. doi: 10.1111/j.1365-3040.2008.01910.x
Mijts, B. N., and Schmidt-Dannert, C. (2003). Engineering of secondary metabolite pathways. Curr. Opin. Biotechnol. 14, 597–602. doi: 10.1016/j.copbio.2003.09.009
Murugan, K., Anandaraj, K., and Al-Sohaibani, S. (2013). Antiaflatoxigenic food additive potential of Murraya koenigii: an in vitro and molecular interaction study. Food Res. Int. 52, 8–16. doi: 10.1016/j.foodres.2013.02.001
Neudörffer, A., Desvergne, J.-P., Bonnefont-Rousselot, D., Legrand, A., Fleury, M.-B., and Largeron, M. (2006). Protective effects of 4-Hydroxycinnamic ethyl Ester derivatives and related Dehydrodimers against oxidation of LDL: radical scavengers or metal Chelators? J. Agric. Food Chem. 54, 1898–1905. doi: 10.1021/jf052923p
Nishihara, M., and Nakatsuka, T. (2011). Genetic engineering of flavonoid pigments to modify flower color in floricultural plants. Biotechnol. Lett. 33, 433–441. doi: 10.1007/s10529-010-0461-z
Nowicka, B., and Kruk, J. (2010). Occurrence, biosynthesis and function of isoprenoid quinones. Biochimica et Biophysica Acta (BBA) Bioenergetics 1797, 1587–1605. doi: 10.1016/j.bbabio.2010.06.007
Oszmiański, J., Lachowicz, S., Gorzelany, J., and Matłok, N. (2018). The effect of different maturity stages on phytochemical composition and antioxidant capacity of cranberry cultivars. Eur. Food Res. Technol. 244, 705–719. doi: 10.1007/s00217-017-2994-z
Pansera, M. R., Santos, A. C. A., Paese, K., Wasum, R., Rossato, M., Rota, L. D., et al. (2003). Análise de taninos totais em plantas aromáticas e medicinais cultivadas no Nordeste do Rio Grande do Sul. Rev. Bras. Farm. 13, 17–22. doi: 10.1590/S0102-695X2003000100002
Pietta, P.-G. (2000). Flavonoids as antioxidants. J. Nat. Prod. 63, 1035–1042. doi: 10.1021/np9904509
Pinhatti, A. V., de Matos Nunes, J., Maurmann, N., Rosa, L. M. G., von Poser, G. L., and Rech, S. B. (2010). Phenolic compounds accumulation in Hypericum ternum propagated in vitro and during plant development acclimatization. Acta Physiol. Plant. 32, 675–681. doi: 10.1007/s11738-009-0446-5
Pretti, I. R., Luz, A. C., Jamal, C. M., and Batitucci, M. C. P. (2018). Variation of biochemical and antioxidant activity with respect to the phenological stage of Tithonia diversifolia Hemsl. (Asteraceae) populations. Ind. Crop Prod. 121, 241–249. doi: 10.1016/j.indcrop.2018.04.080
Quideau, S., Deffieux, D., Douat-Casassus, C., and Pouységu, L. (2011). Plant polyphenols: chemical properties, biological activities, and synthesis. Angew. Chem. Int. Ed. 50, 586–621. doi: 10.1002/anie.201000044
Rajendran, M. P., Pallaiyan, B. B., and Selvaraj, N. (2014). Chemical composition, antibacterial and antioxidant profile of essential oil from Murraya koenigii (L.) leaves. Avicenna J. Phytomed. 4, 200–214. doi: 10.22038/AJP.2014.2564
Rana, V. S., Juyal, J. P., Rashmi, S., and Amparo Blazquez, M. (2004). Chemical constituents of the volatile oil of Murraya koenigii leaves. Int. J. Aromather. 14, 23–25. doi: 10.1016/j.ijat.2003.12.007
Rautela, R., Das, G. K., Khan, F. A., Prasad, S., Kumar, A., Prasad, J. K., et al. (2018). Antibacterial, anti-inflammatory and antioxidant effects of Aegle marmelos and Murraya koenigii in dairy cows with endometritis. Livest. Sci. 214, 142–148. doi: 10.1016/j.livsci.2018.05.015
Sandamali, J. A. N., Hewawasam, R. P., Jayatilaka, K. A. P. W., and Mudduwa, L. K. B. (2020). Cardioprotective potential of Murraya koenigii (L.) Spreng. Leaf extract against doxorubicin-induced Cardiotoxicity in rats. Evid. Based Complement. Alternat. Med. 2020:e6023737. doi: 10.1155/2020/6023737
Satyavarapu, E. M., Sinha, P. K., and Mandal, C. (2020). Preclinical development of Mahanine-enriched fraction from Indian spice Murraya koenigii for the Management of Cancer: efficacy, temperature/pH stability, pharmacokinetics, acute and chronic toxicity (14-180 days) studies. Bio. Med. Res. Int. 2020:e4638132. doi: 10.1155/2020/4638132
Shahidi, F., and Ambigaipalan, P. (2015). Phenolics and polyphenolics in foods, beverages and spices: antioxidant activity and health effects – a review. J. Funct. Foods 18, 820–897. doi: 10.1016/j.jff.2015.06.018
Shao, Y., Xu, F., Sun, X., Bao, J., and Beta, T. (2014). Phenolic acids, anthocyanins, and antioxidant capacity in rice (Oryza sativa L.) grains at four stages of development after flowering. Food Chem. 143, 90–96. doi: 10.1016/j.foodchem.2013.07.042
Sharma, O. P., and Bhat, T. K. (2009). DPPH antioxidant assay revisited. Food Chem. 113, 1202–1205. doi: 10.1016/j.foodchem.2008.08.008
Sindhu, R., and Arora, S. (2012). Evaluation of phenolic contents and antioxidant potential of Murraya Koenigii (L) spreng roots. J. App. Pharm. Sci. 2, 120–122. doi: 10.7324/JAPS.2012.21121
Sivakumar, C., and Meera, I. (2013). Antioxidant and biological activities of three morphotypes of Murraya koenigii L. from Uttarakhand. J. Food Process. Technol. 4, 1–7. doi: 10.4172/2157-7110.1000246
Suthar, P., Kumar, S., Kumar, V., Vaidya, D., Sharma, A., and Sharma, A. (2022). Murraya koenigii (L.) Spreng: speculative ethnobotanical perspectives of ubiquitous herb with versatile nutra/functional properties. S. Afr. J. Bot. 145, 111–134. doi: 10.1016/j.sajb.2021.11.025
Swain, T., and Hillis, W. E. (1959). The phenolic constituents of Prunus domestica. I.—the quantitative analysis of phenolic constituents. J. Sci. Food Agric. 10, 63–68. doi: 10.1002/jsfa.2740100110
Tanko, H., Carrier, D. J., Duan, L., and Clausen, E. (2005). Pre-and post-harvest processing of medicinal plants. Plant Genet. Resourc. 3, 304–313. doi: 10.1079/PGR200569
Tobias, S., and Carlson, J. E. (1969). Brief report: Bartlett’s test of Sphericity and chance findings in factor analysis. Multivar. Behav. Res. 4, 375–377. doi: 10.1207/s15327906mbr0403_8
Tomar, M., Bhardwaj, R., Kumar, M., Pal Singh, S., Krishnan, V., Kansal, R., et al. (2021). Nutritional composition patterns and application of multivariate analysis to evaluate indigenous pearl millet Pennisetum glaucum L. R. Br. Germplasm. J. Food Compos. Anal. 103:104086. doi: 10.1016/j.jfca.2021.104086
Uraku, A. J., and Nwankwo, V. O. (2015). Phytochemical and nutritional composition analysis of Murraya koenigii Linn leaves. J. Pharmaceut. Res. Int. 6, 174–180. doi: 10.9734/BJPR/2015/15595
Verma, R. S., Chauhan, A., Padalia, R. C., Jat, S. K., Thul, S., and Sundaresan, V. (2013). Phytochemical diversity of Murraya koenigii (L.) Spreng. From Western Himalaya. Chem. Biodivers. 10, 628–641. doi: 10.1002/cbdv.201200054
Verma, R. S., Padalia, R. C., Arya, V., and Chauhan, A. (2012). Aroma profiles of the curry leaf, Murraya koenigii (L.) Spreng. Chemotypes: variability in North India during the year. Ind. Crop Prod. 36, 343–348. doi: 10.1016/j.indcrop.2011.10.022
Verma, N., and Shukla, S. (2015). Impact of various factors responsible for fluctuation in plant secondary metabolites. J. Appl. Res. Med. Aromat. Plants 2, 105–113. doi: 10.1016/j.jarmap.2015.09.002
Vlaisavljević, S., Kaurinović, B., Popović, M., and Vasiljević, S. (2017). Profile of phenolic compounds in Trifolium pratense L. extracts at different growth stages and their biological activities. Int. J. Food Prop. 20, 3090–3101. doi: 10.1080/10942912.2016.1273235
Yadav, S., Vats, V., Dhunnoo, Y., and Grover, J. K. (2002). Hypoglycemic and antihyperglycemic activity of Murraya koenigii leaves in diabetic rats. J. Ethnopharmacol. 82, 111–116. doi: 10.1016/S0378-8741(02)00167-8
Yadav, V., Wang, Z., Wei, C., Amo, A., Ahmed, B., Yang, X., et al. (2020). Phenylpropanoid pathway engineering: An emerging approach towards plant defense. Pathogens 9:312. doi: 10.3390/pathogens9040312
Yadegarinia, D., Gachkar, L., Rezaei, M. B., Taghizadeh, M., Astaneh, S. A., and Rasooli, I. (2006). Biochemical activities of Iranian Mentha piperita L. and Myrtus communis L. essential oils. Phytochemistry 67, 1249–1255. doi: 10.1016/j.phytochem.2006.04.025
Yang, L., Wen, K.-S., Ruan, X., Zhao, Y.-X., Wei, F., and Wang, Q. (2018). Response of plant secondary metabolites to environmental factors. Molecules 23:762. doi: 10.3390/molecules23040762
Zielinski, A. A. F., Alberti, A., and Nogueira, A. (2018). “Chapter 39- characterizing fruit juices and fermented fruit beverages using Chemometrics tools,” in Fruit Juices. eds. G. Rajauria and B. K. Tiwari (San Diego: Academic Press), 823–833.
Keywords: Murraya koenigii, developmental stages, multivariate data analysis, HCA, PCA, heat map, phenols, antioxidant
Citation: Verma R, Singh N, Tomar M, Bhardwaj R, Deb D and Rana A (2022) Deciphering the growth stage specific bioactive diversity patterns in Murraya koenigii (L.) Spreng. using multivariate data analysis. Front. Plant Sci. 13:963150. doi: 10.3389/fpls.2022.963150
Edited by:
Sezai Ercisli, Atatürk University, TurkeyReviewed by:
Emelda Rohani, Universiti Kebangsaan Malaysia, MalaysiaGülçe Ilhan, Atatürk University, Turkey
Copyright © 2022 Verma, Singh, Tomar, Bhardwaj, Deb and Rana. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Maharishi Tomar, maharishi89@gmail.com