Next Article in Journal
Genotype and Seasonal Variation Affect Yield and Oil Quality of Safflower (Carthamus tinctorius L.) under Mediterranean Conditions
Previous Article in Journal
Ecosystem Services Provided by Cover Crops and Biofumigation in Sunflower Cultivation
Previous Article in Special Issue
The Perfect Match: Adjusting High Tree Density to Rootstock Vigor for Improving Cropping and Land Use Efficiency of Sweet Orange
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Flavonoids and Limonoids Profiles Variation in Leaves from Mandarin Cultivars and Its Relationship with Alternate Bearing

by
Cecilia Rodríguez Ceraolo
1,
Valeria Vázquez
1,
Ignacio Migues
1,2,
María Verónica Cesio
1,
Fernando Rivas
3 and
Horacio Heinzen
1,*
1
Grupo de Análisis Metabolómicos (GAMe), Laboratorio de Farmacognosia y Productos Naturales, Departamento de Química Orgánica, Facultad de Química, Universidad de la República, Montevideo 11800, Uruguay
2
Graduate Program in Chemistry, Facultad de Química, Universidad de la República, Montevideo 11800, Uruguay
3
Programa Nacional de Investigación en Produccción Citrícola, Instituto Nacional de Investigación Agropecuaria (INIA), Salto 50000, Uruguay
*
Author to whom correspondence should be addressed.
Submission received: 31 October 2021 / Revised: 22 November 2021 / Accepted: 2 December 2021 / Published: 4 January 2022

Abstract

:
Alternate bearing in citrus trees has been extensively studied as a key feature for citrus growers. Although the genetic and the biochemical process occurring during alternate bearing has been studied extensively, there is a lack of information identifying the presence of metabolic indicators during “on” and “off” years. In citrus plants, leaves play a central role in the metabolic pathway triggering the flowering induction process. To investigate the changes during this transition, a liquid chromatography coupled to tandem mass spectrometry (LC-MS/MS) analysis of the leaf profiles of 20 compounds (17 polyphenols, two limonoids, and one furanocoumarin), in bearing and non-bearing branches arising from four different mandarin genotypes, was performed. The same metabolites were found in all the genotypes at both stages: both limonoids and 11 polyphenols. Using these compounds, the chemotaxonomic differentiation between cultivars was assessed. The levels of flavanones and limonoids showed differences in both bearing stages and the transition from vegetative to flowering could be shown by the activation of the polyphenol biosynthetic pathway, from precursors like naringenin to metabolic end-points such as narirutin and polymethoxyflavones. Narirutin levels showed significant differences between both stages, suggesting it as a possible marker of the physiological status of the branch.

1. Introduction

Citrus is one of the main fruit crops worldwide with more than 143 million tons produced in 2019 [1]. Alternate bearing in Citrus is one of the major limiting factors affecting production profits and, thus, tackling the maximal potential benefits for the global citrus industry is essential. Therefore, the genetic, biochemical, and metabolic determinations of alternate bearing have been widely studied, aiming to gain knowledge about the regulatory process in the transition from bearing to non-bearing trees. However, the mechanism by which the developing crop influences blooming and fruit yield the forthcoming year has not been fully understood yet. One hypothesis holds that return bloom and yield are related to tree carbohydrate status and, thus, lack of carbohydrate in the “on” year directly or indirectly reduces flowering the following year [2]. Other authors claim that developing fruit produces an inhibitor that directly or indirectly reduces flowering in the spring following the “on” crop [2,3]. Although extensive research has been conducted to unravel the hormonal and genetic regulation of flowering induction in Citrus [4,5], there is still a lack of information describing metabolites’ changes (in terms of quality and quantity) occurring during “on” and “off” years as putative markers depicting this transition.
Flowering represents a transition phase from vegetative to reproductive growth and, thus, it is one of the most important events in the plant life cycle. During vegetative to reproductive transition, citrus undergoes metabolic regulations that have been extensively reported [6].
Flowering is regulated by endogenous and exogenous factors. The first one has been mainly attributed to the genetic and hormonal modulation of the LFY gen and Gibberellins metabolic pathway [7,8,9], identifying the fruit-mediated inhibition of flowering as the main cause of alternate bearing in citrus [10]. Similarly, carbon metabolism has also been reported to be involved in flower induction mediated by fruit presence in citrus [2]. Exogenous regulation of flowering has been mainly attributed to environmental stress factors including light intensity and wavelength, temperature, mineral and carbon nutrition, drought, low oxygen, crowding, root removal, and mechanical stimulation, as well as biotic stress factors such as infection by pests and pathogens [11,12,13].
A possible approach to find the metabolic changes involved in the process of alternate bearing is the study of the metabolome of the leaves in both stages. In this way, polyphenols and limonoids fulfill definite biological and ecological properties in plants [14]. As an example, when the carbon/nutrient ratio is high, photosynthesis is primarily diverted to the production of carbon-based metabolites such as phenolics [15,16] and, thus, high levels of phenolics are associated with in-excess accumulation of leaf carbohydrates. In this way, phenyl-propanoid biosynthesis serves as an ‘energy overflow’ mechanism by diverting photosynthate and cellular reducing power into stable product pools [17], preventing further oxidative damage during photoinhibition caused by high carbohydrates accumulation in leaves [18].
Polyphenols are a remarkably diverse group of secondary products with a vast array of biological functions. They protect plants from different biotic and abiotic stress and act as unique UV filters, signal molecules, allopathic compounds, phytoalexins, detoxifying agents and antimicrobial defensive molecules [19]. They are classified in to two major classes: flavonoids and nonflavonoids. The latter include structurally simple molecules such as phenolic acids (hydroxybenzoic acids and hydroxycinnamic acids) and stilbenes, along with highly complex molecules such as stilbene oligomers, tannins, and lignans [20]. The former is the most studied class of polyphenols, sharing the same general C6-C3-C6 structure. The most abundant have a flavan nucleus, consisting of two aromatic rings (A and B) linked by a pyran ring (C). Differences in the location of the B-ring to C-ring linkage make it possible to distinguish between flavonoids (2-phenylbenzopyrans), isoflavonoids (3-phenylbenzopyrans), and neoflavonoids (4-phenylbenzopyrans). The first one is by far the most abundant group and may be further divided into 3-hydroxyflavonoids (flavonols, flavanols, anthocyanidins, and dihydroflavonols), and flavonoids without substituent at C3 (flavanones and flavones). The flavonoids may be among the most important aromatic compounds. In this regard, they are the most ancient and widespread of the polyphenols, synthesized even in mosses and ferns, and have a wide range of potent physiological activities [21].
Citrus limonoids (CLs) are a crowded family of polycyclic secondary metabolites found mostly in seeds, fruits, and peel tissues of citrus fruits, including lemons, limes, oranges (both sour and sweet), pummels, grapefruits, bergamots, and mandarins [22]. CLs are highly oxygenated and modified triterpenes classified as tetranotriterpenoids, found as both free aglycones and their corresponding β-D-glycosides. Up to date, about 55 limonoid aglycones and 18 limonoid glycosides were identified in Citrus. Limonin is the most abundant aglycone in most citrus species followed by Nomilin, and limonin glycoside is the most representative limonoid glycoside [23]. CLs possess a wide array of biological activities [23], and have been acknowledged as high anti-insect activity compounds, both by their toxicity and antifeedant properties [22,24,25,26,27,28]. They probably protect leaflets from insects during bud development, justifying their abundancy in young leaves and fruit [29]. The first studies of CLs originate on the attempt to find bitter compounds in citrus fruit [30,31,32], but few reports on the occurrence of CLs in citrus tree leaves are found.
Metabolomic research strategies can be faced by two distinct approaches, untargeted and targeted, each with their own inherent advantages and disadvantages. Untargeted metabolomics is the comprehensive analysis of all measurable analytes in a sample. Although it offers the opportunity for novel target compounds discovery, the importance of rare molecules can be underestimated [33]. In contrast, the targeted approach is the precise and quantitative measurement of the concentrations of a defined group of chemically characterized metabolites. Hypothesis-driven targeted experiments can go together with the risk of losing the real natural products that are responsible for an observed biological property. For targeted analysis, sample preparation is optimized for obtaining a defined group of natural products.
Seeking the identification of possible representative phytomarkers of the transition from “on” to “off” bearing, the profiling of polyphenols and limonoids in leaves from bearing and non-bearing branches was studied using a targeted metabolomic approach.

2. Materials and Methods

2.1. Plant Material

The study was carried out at the National Agricultural Research Institute (INIA) of Uruguay, in Salto Province (31°16′ S; 57°53′ W, altitude: 47 m), on adult trees from cv. ‘Ellendale’ (Citrus reticulata Blanco), cv. ‘Satsuma’ (Citrus unshiu Marc.), cv. ‘Willowleaf’ (Citrus deliciosa Ten.), and cv. ‘Page’ (Citrus reticulata Blanco × Citrus paradisi Macf.) grafted onto Poncirus trifoliata (L.) Raf. rootstock, under accurate horticulture cares (irrigation and pest management). From each variety, five bearing and five non-bearing trees were sampled. Adult leaves between 6 and 12 months old (5 leaves/branch; 10 branches sampled) were collected from bearing (B) and non-bearing (NB) branches. In bearing branches, leaves were collected immediately below the fruits (5–20 cm from the peduncle) during full fruit color stage, whereas in non-bearing branches, leaves were selected to have the same characteristics and position as in bearing branches. All trees were located in the same orchard and growing conditions. Therefore, the differences in the composition between the varieties could not be ascribed to changes in the tree environment, but due to the genetic background and physiological state of each tree. In all the branches for sampling, leaves were taken randomly all around the outer part of the tree canopy, exposed to full sunlight, and immediately stored at −20 °C until extraction.

2.1.1. Reagents and Equipment

All mobile phases were prepared using High Performance Liquid Chromatography-grade solvents, supplied by Merck KGaA (Darmstadt, Germany), as well as analytical grade solvents used for extraction. Ultrapure water (Millipore Milli-Q® Ultrapure Water Solutions Type 1 by Merckmillipore, Molsheim, France) was used for the extraction and HPLC analysis. All standards were selected based on previous studies of flavonoid content in mandarins [34] and purchased from PhytoLab GmbH & Co.KG (Vestenbergsgreuth, Germany).
Polyphenols and limonoids were analyzed in an Agilent 1200 LC system (Agilent Technologies, Palo Alto, CA, USA) coupled to a 4000 QTRAP® LC-MS/MS System from AB SCIEX™ (Framingham, MA, USA). LC-Separation was performed on a Gold Aq-Hypersil (5 µm 150 × 4.6 mm) column from Thermo Scientific™ (USA). The injection volume was 5 μL. MS/MS detection was performed in the multiple reaction monitoring (MRM) mode using an electrospray ionization (ESI) interface in both positive and negative ion mode.

2.1.2. Standard Solutions Preparation

Individual standard solutions of the target compounds were prepared and stored at −20 °C. Each standard was weighed to approximately 20 mg and dissolved in 10.00 mL of HPLC-grade methanol, to obtain a 2000 mg/L stock solution. A working solution mix of all targeted compounds (20 mg/L) was prepared by appropriate dilution of the stock solutions in acetonitrile and used for the calibration curve construction. The selected concentrations for the curve were 10, 50, 100, 200, 400, 600, 800, 1000, and 1500 μg/L.

2.1.3. Sample Preparation

The extraction procedure for mandarin leaves was performed using an Ultrasound-assisted extraction (UAE) adapted from Safdar M.N. et al. [35]. A total of 50 g of freeze-dried leaves was pool together and ground using a cut mill for each variety. Then, 1 g of the powdered leaf samples was extracted with methanol 80% in a sample/solvent ratio 1:10 for 15 min in an Ultrasound bath (Transsonic 460; Elma, Wetzikon, Switzerland) at 35 kHz frequency. This procedure was repeated twice using fresh solvent. The extracts were combined and subjected to bulk filtration, solvent vacuum evaporation, and lyophilization, followed by redissolution in 80% methanol HPLC grade and filtration through 0.22 μm Polyvinylidene fluoride (PVDF) filters.

2.1.4. Chromatographic and MS/MS Conditions

The mobile phase was ultrapure water with 0.1% formic acid (A) and acetonitrile (B), with the subsequent gradient: 0–23 min 5% to 50% B, 23–26 min 50% to 100% B, 26–31 100% B, 31–33 min 100% to 5% B, 33–38 5% B, at 600 μL/min. The ionization voltage was 5000 V, the nebulizer and the curtain gases were nitrogen at 50 psi each. The solvent evaporation in the source was assisted by a drying gas (heated nitrogen at 500 °C per 50 psi).
A multiple reaction monitoring (MRM) mode method was used for quantification. The optimal MRM transitions, collision energies (CE), and declustering potentials (DP) for each investigated compound were determined infusing the standard solutions directly into the instrument at a constant flow, using a syringe.

2.1.5. Data Treatment

The statistical treatment was performed using Metaboanalyst 5.0 [36] and R-Studio software package (version 1.1.463). Each of the studied compounds were considered as a single variable and were included in all the treatments. Data were normalized using pareto scaling. ANOVA, Wilcoxon/Mann–Whitney, and multiple comparisons test using a Bonferroni correction were performed. Principal component analysis (PCA) and the supervised test partial least squares discriminant analysis (PLS-DA) were used for the classification of the varieties. The latter was validated to avoid overfitting using cross-validation and permutation test. Heatmap based on Student’s t-test was also performed.

3. Results and Discussion

3.1. Optimization of MS/MS Parameters of the Target Compounds

The acquisition method was ESI positive for limonoids and polymethoxyflavones, and negative for flavanones, flavone, and flavonols (Table 1). The chromatogram of the optimized compounds is shown in Supplementary Materials.

3.2. Identification and Quantification of Phytochemicals in Leaves

A total of 11 polyphenols were identified and quantified among the four genotypes, together with two limonoids. The identified phytochemicals were the flavanones didymin, eriocitrin, naringenin, hesperidin, and narirutin; two polimethoxyflavones, tangeretin and nobiletin; a flavonol, quercitrin, and the flavone luteolin, as well as the hydroxycinnamic acids chlorogenic and caffeic acid and the limonoids limonin and nomilin (Figure 1). The chromatogram of the identified compounds is shown in Supplementary Materials.

3.3. Chemotaxonomic Characterization of Bearing Branches

Leaves from bearing branches of the four studied genotypes contained high concentrations of hesperidin (2824.43 μg/g *), followed by nobiletin (583.71 μg/g *), tangeretin (326.60 μg/g *), narirutin (63.22 μg/g *), and eriocitrin (37.23 μg/g *). Meanwhile, the most abundant limonoid was limonin (0.69 μg/g *). The phytochemical concentration of each variety is displayed in Table 2. * Average of the four genotypes.
A PCA was performed to check the grouping of the samples and detect similarities and differences between the four varieties. Figure 2a shows the 2D plot of the PCA, where PC1 and PC2 explain 97.7% of the total data variance. It shows a good classification of all genotypes, with ‘Page’ and ‘Satsuma’ being the ones that have the most similarities between each other.
After performing a PLS-DA, the top three scores of the Variables Importance in the Projection (VIP) are nobiletin for the differentiation of END-PG, tangeretin for END-PG, and didymin for PG-WL (Figure 2b).

3.4. Chemotaxonomic Characterization of Non-Bearing Branches

The same polyphenols and limonoids found in bearing branches were found in the non-bearing branches’ leaves but with different concentrations (Table 3). Narirutin was the principal flavonoid with the highest concentration for non-bearing ‘Page’ and ‘Willowleaf’ (42198.36 μg/g in NB-PG and 10049.66 μg/g in NB-WL, respectively). Hesperidin is the most abundant flavonoid for the non-bearing ‘Satsuma’ and ‘Ellendale’ (2591.54 μg/g in NB-SAT and 2433.21 μg/g in NB-END, respectively), as shown also in bearing branches.
Non-bearing branches of ‘Page’ were characterized by the highest luteolin, narirutin, and quercitrin levels, a situation which was also observed in ‘Willowleaf’. In ‘Satsuma’, eriocitrin and chlorogenic acid were characteristic of this genotype, whereas in ‘Ellendale’, the polymethoxyflavones tangeretin and nobiletin were the metabolites that characterize the leaves of non-bearing branches.
In the PCA plot, it can be seen that the four genotypes are well classified from each other. PC1 and PC2 accounted for the 97.3% of the observed variance (Figure 3a). ‘Willowleaf’ and ‘Page’ showed more dispersion compared to the other two varieties due to the variability in the levels of some flavonoids.
The concentration of both polymethoxyflavones were different, being highest to lowest: NB-END, NB-WL, NB-SAT, and NB-PG. The flavanone hesperidin was highest for the duad WL-PG; the same behavior was observed for flavanones naringin and naringenin, flavonols, and flavones. Didymin had high concentration values in NB-END and NB-PG, but not in NB-WL and NB-SAT. Hydroxycinnamic acids are present in the highest concentration in NB-SAT variety (Figure 3a, Table 2 and Table 3).
The most important compounds according to VIP scores include luteolin, hesperidin, and naringenin. (Figure 3b). ‘Willowleaf’ has the highest concentration of luteolin, and ‘Page’ of hesperidin and naringenin.

3.5. Metobolites’ Profile Comparison between Leaves from Bearing and Non-Bearing Branches

The chemotaxonomic characterization of the four genotypes in both bearing stages (“on” and “off”) allowed the identification of characteristic compounds that could differentiate between bearing stages. In order to determine if there is a significant difference between the means of these two groups, a Wilcoxon/Mann–Whitney test for each compound between bearing and non-bearing leaves was performed for each genotype. Figure 4 shows that nobiletin is not useful in differentiating between bearing and non-bearing stages, for any genotype. Hesperidin, luteolin, naringenin, narirutin, and quercitrin are the compounds that could differentiate between stages; however, it depends on the genotype considered. For example, narirutin is useful for differentiating all the genotypes except for ‘Satsuma’. Using this test, a set of compounds could not be selected to classify both stages for all the varieties, due to their high variability. Instead, a pairwise study of the “on”–“off” physiological situations of each cultivar showed significant differences between them. These facts could be ascribed to the different genetic origins of each genotype. Interestingly, there were significant differences in eriocitrin levels between both physiological stages for three of the genotypes excepting ‘Satsuma’.
In order to identify the compounds responsible for the differentiation of stages in all the genotypes, other statistical approaches were performed and, thus, a Heatmap based on a Student’s t-test was developed and shown in Figure 5a. Narirutin, hesperidin, luteolin, naringenin, and quercitrin could be used to differentiate both stages in ‘Page’ and ‘Willowleaf’, as the non-bearing samples have the highest concentration of these compounds. Interestingly, caffeic acid is a common marker for the bearing branches of the other varieties. In Figure 5b, the VIP scores of the PLS-DA performed are shown. This analysis allowed the identification of the compounds responsible for the classification of the “on” an “off” bearing stages. Narirutin and luteolin are the variables with the highest VIP scores, and could be identified as candidate markers of alternate bearing.
‘Willowleaf’ and ‘Page’ showed clear differences in the concentration level of some metabolites between “on” and “off” phases, that could easily distinguish both bearing stages. Quercitrin presented two orders of magnitude higher concentration in “off” than in “on” phase. The ratio hesperidin/narirutin could be indicative of the bearing stage of ‘Willowleaf’ and ‘Page’, when it is <1 the tree is in the non-bearing phase, whereas when the quotient is >10, the tree is at the bearing stage.
Some differences in the metabolites’ profile detected between “on” and “off” mandarin branches are due to the position of the biosynthetic pathways of the prevailing metabolites found in them. For instance, nomilin is the biosynthetic precursor of Limonin [37,38]. Limonoids have not been studied as markers of alternate bearing in mandarin trees. In the case of the limonoid biosynthetic pathway the sequence Nomilin–Limonin reaches higher concentration levels during the non-bearing phase, indicating the possibility that the tree is preparing to protect itself against pest and insects before the flowering stage. Naringenin, which is the first polycyclic flavonoid skeleton, is found in trace amounts in leaves arising from non-bearing branches of cvs. ‘Satsuma’ and ‘Ellendale’, but in significant levels in ‘Page’ and “Willowleaf’ cultivars. In bearing branches naringenin was found in very low quantities in all four varieties. In the former case, the biosynthetic pathway to flavanones is favored by 7-O-glycosylation with hesperoside, a disaccharide [1→6]glucose-rhamnose producing narirutin, which is the prevailing flavanone in non-bearing leaves of ‘Page’ and ‘Willowleaf’ [39,40].
From naringenin, hesperidin and narirutin are formed. Narirutin seems to be a marker of non-bearing trees, as it showed significant differences between bearing and non-bearing samples, particularly in ‘Page’ and ‘Willowleaf’. The biosynthesis of both hesperosides (Figure 6) goes through the 7-O-glycosidation of naringenin, which is then specifically glycosidated in the position 6-O of the sugar moiety by a rhamonsyl residue. In the non-bearing phase, when the tree is recovering after the energy sink caused by fruiting, the glycosidation reaction yielded higher levels of the naringenin glycosides. When the fruit is developing, the metabolism is fully expressed and active, and all the naringenin is driven through the oxydation-methylation pathway to hesperidin, the most ubiquitous citrus flavonoid.
Interestingly Mou et al. [41], in an attempt to decipher the flavonoid biosynthetic pathway in C. reticulata, did not describe the pathway leading to narirutin, but they described the naringenin–naringin route. The mandarin species studied produce the tasteless narirutin instead of the bitter naringin. The difference between both compounds is the union between the rhamnose and the glucose moiety, in narirutin the link is (1→6) (rutinoside), whereas in naringin the link is (1→2) (neohesperoside).
Another route downstream of the flavonoid biosynthesis in mandarin trees leads to the polymethoxyflavones (PMF), suggested as markers of bearing stages [42,43]. Luteolin, an acknowledged precursor of PMF is also distinctive polyphenol of the non-bearing phase.
The general trend discussed above, allowed the assessment of qualitative and quantitative differences between the four genotypes, depending on the physiological status of the branch. Flavonoid biosynthesis in citrus is regulated by a combination of environmental, genetic, and hormonal factors. As the former is the same for all the individuals in this study, the two latter variables are the ones responsible for the observed differences. It has been accepted that in “off” branches, buds are directed to morph to leaves instead of flowers, enhancing primary and secondary metabolism, which confirms the general trend identified by Shalom et al. [2], that flavonoid biosynthesis is enhanced during the “off” phase in ‘Murcott’ mandarins. They found that the expression of all the enzymes of the flavonoid pathway were significatively higher in “off” branches, as well as the expression of mRNA from FT, AP1, and LFY related genes, which was correlated with the up-regulation effect of these genes in the “off” phase. The obtained results of our study in four different mandarin cultivars are in accordance with the expression of the genes described by Shalom [2]. Moreover, high levels of abscisic acid, which is high in buds, leaves, and fruits during autumn in ‘Murcott’ (“on” phase) and other alternating citrus varieties, have been shown to up regulate the flavonoid biosynthesis in another non-climateric fruits, such as blueberries and tomatoes [44,45].
Furthermore, Narirutin is identified as a possible marker of the “on/off” phase as it presents statistically different levels in all varieties during both stages. The only exception is ‘Satsuma’, a variety that does not present the typical alternate bearing behavior as exhibited by the rest of the studied varieties.

4. Conclusions

A specific analytical protocol using LC-MS-(QqQ)-MS, applying different acquisition modes, permitted the unambiguous determination of 13 secondary metabolites in leaves of mandarin trees. A chemosystematic study of the profile of these metabolites in leaves from non-bearing (“off”) and bearing (“on”) mandarin branches belonging to four different genotypes was performed. The different genotypes could be chemotaxonomically distinguished from each other at both developmental stages using multivariate data analysis and the representative metabolites were identified.
This study allowed the differentiation between leaves from non-bearing and bearing branches. Leaves from both physiological stages accumulated limonin, but nomilin showed an increased leaf concentration in non-bearing branches compared to bearing ones. “On” and “off” ‘Willowleaf’ and ‘Page’ mandarin trees can be distinguished due to the high narirutin content in leaves of non-bearing branches, and an easy relationship has been proposed to define each phase. One of the most important observations found is the increase of Narirutin content in non-bearing leaves for genotypes that show a typical alternate bearing behavior which could lead to the use of this compound as a biomarker of the stage of the branch. The presence of flavonoids as markers of the different physiological stages of mandarin trees also hints to possible new ways of crop management. Future studies should be performed following the developed methodology to control the levels of these compounds, specially Narirutin, during the ontogeny of mandarin branches along the plant growth.
The results reported in the present research can contribute towards a useful technology for growers and scientists to determine the physiological stage of a tree depicting the transition from “off” to “on” phase in Citrus.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/agronomy12010121/s1, Table S1: Target compounds selected for the analysis, Figure S1: Chromatogram of the mix of standards at 800 g/L injected into the LC-MS/MS in negative (Top) and positive 46 (Bottom) mode, and Figure S2: Chromatograms of ‘Willowleaf’ genotype at non-bearing (“off”) stage (A) and at a bearing (“on”) stage (B) in negative acquisition mode.

Author Contributions

Conceptualization, H.H. and I.M.; methodology, C.R.C. and V.V.; validation, C.R.C. and I.M.; formal analysis, C.R.C., I.M. and V.V.; investigation, C.R.C., I.M. and H.H.; writing—original draft preparation, C.R.C. and I.M.; writing—review and editing, H.H. and F.R.; supervision, H.H. and M.V.C.; project administration, H.H. and F.R.; funding acquisition, F.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by INIA through the L4 funding program, CSIC grant number 1822, the Agencia Nacional de Investigación e Innovación grant number POS_NAC_2017_1_140316, and the Programa de Desarrollo de las Ciencias Básicas.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

We would like to thank Natalia Besil for the assistance during the optimization of the MS analyses.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. FAOSTAT. Available online: http://www.fao.org/faostat/en/#data/QCL/visualize (accessed on 5 October 2021).
  2. Shalom, L.; Samuels, S.; Zur, N.; Shlizerman, L.; Zemach, H.; Weissberg, M.; Ophir, R.; Blumwald, E.; Sadka, A. Alternate Bearing in Citrus: Changes in the Expression of Flowering Control Genes and in Global Gene Expression in ON- versus OFF-Crop Trees. PLoS ONE 2012, 7, e46930. [Google Scholar] [CrossRef] [Green Version]
  3. Agust, M.; Almela, V.; Pons, J. Effects of girdling on alternate bearing in citrus. J. Hortic. Sci. 1992, 67, 203–210. [Google Scholar] [CrossRef]
  4. Talon, M.; Tadeo, F.R.; Ben-Cheikh, W.; Gomez-Cadenas, A.; Mehouachi, J.; Pérez-Botella, J.; Primo-Millo, E. Hormonal regulation of fruit set and abscission in citrus: Classical concepts and new Evidence. Acta Hortic. 1998, 463, 209–217. [Google Scholar] [CrossRef]
  5. Bower, J.P.; Cutting, J.G.M.; Lovatt, C.J.; Blanke, M.M. Interaction of plant growth regulator and carbohydrate in flowering and fruit set. Acta Hortic. 1990, 425–434. [Google Scholar] [CrossRef]
  6. Davenport, T.L. Citrus Flowering. In Horticultural Reviews; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2011; pp. 349–408. [Google Scholar]
  7. Peña, L.; Martín-Trillo, M.; Juárez, J.; Pina, J.A.; Navarro, L.; Martínez-Zapater, J.M. Constitutive expression of Arabidopsis LEAFY or APETALA1 genes in citrus reduces their generation time. Nat. Biotechnol. 2001, 19, 263–267. [Google Scholar] [CrossRef]
  8. Pillitteri, L.J.; Lovatt, C.J.; Walling, L.L. Isolation and characterization of LEAFY and APETALA1 homologues from Citrus sinensis L. Osbeck “Washington”. J. Am. Soc. Hortic. Sci. 2004, 129, 846–856. [Google Scholar] [CrossRef] [Green Version]
  9. Araki, T. Transition from vegetative to reproductive phase. Curr. Opin. Plant Biol. 2001, 4, 63–68. [Google Scholar] [CrossRef]
  10. Moss, G.I. Effect of fruit on flowering in relation to biennial bearing in sweet orange (Citrus sinensis). J. Hortic. Sci. 1971, 46, 177–184. [Google Scholar] [CrossRef]
  11. Wada, K.C.; Takeno, K. Stress-induced flowering. Plant Signal. Behav. 2010, 5, 944–947. [Google Scholar] [CrossRef] [PubMed]
  12. Takeno, K. Stress-Induced Flowering. Abiotic Stress Responses Plants 2012, 331–345. [Google Scholar] [CrossRef]
  13. Kazan, K.; Lyons, R. The link between flowering time and stress tolerance. J. Exp. Bot. 2016, 67, 47–60. [Google Scholar] [CrossRef] [Green Version]
  14. Zagoskina, N.V.; Kazantseva, V.V.; Fesenko, A.N.; Shirokova, A.V. Accumulation of Phenolic Compounds at the Initial Steps of Ontogenesis of Fagopyrum esculentum Plants That Differ in Their Ploidy Levels. Biol. Bull. 2018, 45, 171–178. [Google Scholar] [CrossRef]
  15. Coley, P.D.; Bryant, J.P.; Chapin, F.S. Resource Availability and Plant Antiherbivore Defense. Science 1985, 230, 895–899. [Google Scholar] [CrossRef] [Green Version]
  16. Waterman, P.; Mole, S. Analysis of Phenolic Plant Metabolites; Blackwell Scientific Publications: Oxford, UK, 1994. [Google Scholar]
  17. Grace, S.G.; Logan, B.A. Energy Dissipation and Radical Scavenging by the Plant Phenylpropanoid Pathway. Philos. Trans. R. Soc. Lond. B Biol. Sci. 2000, 355, 1499–1510. [Google Scholar] [CrossRef] [Green Version]
  18. Rivas, F.; Fornes, F.; Agustí, M. Girdling induces oxidative damage and triggers enzymatic and non-enzymatic antioxidative defences in Citrus leaves. Environ. Exp. Bot. 2008, 64, 256–263. [Google Scholar] [CrossRef]
  19. Panche, A.N.; Diwan, A.D.; Chandra, S.R. Flavonoids: An overview. J. Nutr. Sci. 2016, 5, 1–15. [Google Scholar] [CrossRef] [Green Version]
  20. Cheynier, V. Polyphenols in foods are more complex than often thought. Am. J. Clin. Nutr. 2005, 81, 223S–229S. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. McClure, J.W. Physiology and Functions of Flavonoids. In The Flavonoids; Harbone, J.B., Mabry, T.J., Mabry, H., Eds.; Springer: Boston, MA, USA, 1975; pp. 970–1055. [Google Scholar]
  22. Gualdani, R.; Cavalluzzi, M.; Lentini, G.; Habtemariam, S. The Chemistry and Pharmacology of Citrus Limonoids. Molecules 2016, 21, 1530. [Google Scholar] [CrossRef]
  23. Shi, Y.S.; Zhang, Y.; Li, H.T.; Wu, C.H.; El-Seedi, H.R.; Ye, W.K.; Wang, Z.W.; Li, C.B.; Zhang, X.F.; Kai, G.Y. Limonoids from Citrus: Chemistry, anti-tumor potential, and other bioactivities. J. Funct. Foods 2020, 75, 104213. [Google Scholar] [CrossRef]
  24. Roy, A.; Saraf, S. Limonoids: Overview of Significant Bioactive Triterpenes Distributed in Plants Kingdom. Biol. Pharm. Bull. 2006, 29, 191–201. [Google Scholar] [CrossRef] [Green Version]
  25. Hafeez, F.; Akram, W.; Shaalan, E.A.-S. Mosquito larvicidal activity of citrus limonoids against Aedes albopictus. Parasitol. Res. 2011, 109, 221–229. [Google Scholar] [CrossRef] [PubMed]
  26. Warikoo, R.; Ray, A.; Sandhu, J.K.; Samal, R.; Wahab, N.; Kumar, S. Larvicidal and irritant activities of hexane leaf extracts of Citrus sinensis against dengue vector Aedes aegypti L. Asian Pac. J. Trop. Biomed. 2012, 2, 152–155. [Google Scholar] [CrossRef] [Green Version]
  27. Khan, M.K.; Dangles, O. A comprehensive review on flavanones, the major citrus polyphenols. J. Food Compos. Anal. 2014, 33, 85–104. [Google Scholar] [CrossRef]
  28. Yu, X.; Shi, D.; Zhi, X.; Li, Q.; Yao, X.; Xu, H. Synthesis and quantitative structure-activity relationship (QSAR) study of C7-oxime ester derivatives of obacunone as insecticidal agents. RSC Adv. 2015, 5, 31700–31707. [Google Scholar] [CrossRef]
  29. Hasegawa, S.; Berhow, M.A.; Manners, G.D. Citrus limonoid research: An overview. ACS Symp. Ser. 2000, 758, 1–8. [Google Scholar] [CrossRef] [Green Version]
  30. Celano, R.; Campone, L.; Pagano, I.; Carabetta, S.; Di Sanzo, R.; Rastrelli, L.; Piccinelli, A.L.; Russo, M. Characterisation of nutraceutical compounds from different parts of particular species of Citrus sinensis ‘Ovale Calabrese’ by UHPLC-UV-ESI-HRMS. Nat. Prod. Res. 2019, 33, 244–251. [Google Scholar] [CrossRef]
  31. Russo, M.; Arigò, A.; Calabrò, M.L.; Farnetti, S.; Mondello, L.; Dugo, P. Bergamot (Citrus bergamia Risso) as a source of nutraceuticals: Limonoids and flavonoids. J. Funct. Foods 2016, 20, 10–19. [Google Scholar] [CrossRef]
  32. Wang, F.; Yu, X.; Liu, X.; Shen, W.; Zhu, S.; Zhao, X. Temporal and spatial variations on accumulation of nomilin and limonin in the pummelos. Plant Physiol. Biochem. 2016, 106, 23–29. [Google Scholar] [CrossRef] [PubMed]
  33. Roberts, L.D.; Souza, A.L.; Gerszten, R.E.; Clish, C.B. Targeted metabolomics. Curr. Protoc. Mol. Biol. 2012, 1, 30.2.1–30.2.24. [Google Scholar] [CrossRef]
  34. Xi, W.; Zhang, Y.; Sun, Y.; Shen, Y.; Ye, X.; Zhou, Z. Phenolic composition of Chinese wild mandarin (Citrus reticulata Balnco.) pulps and their antioxidant properties. Ind. Crops Prod. 2014, 52, 466–474. [Google Scholar] [CrossRef]
  35. Safdar, M.N.; Kausar, T.; Jabbar, S.; Mumtaz, A.; Ahad, K.; Saddozai, A.A. Extraction and quantification of polyphenols from kinnow (Citrus reticulate L.) peel using ultrasound and maceration techniques. J. Food Drug Anal. 2017, 25, 488–500. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Pang, Z.; Chong, J.; Zhou, G.; de Lima Morais, D.A.; Chang, L.; Barrette, M.; Gauthier, C.; Jacques, P.-É.; Li, S.; Xia, J. MetaboAnalyst 5.0: Narrowing the gap between raw spectra and functional insights. Nucleic Acids Res. 2021, 49, 388–396. [Google Scholar] [CrossRef] [PubMed]
  37. Ou, P.; Hasegawa, S.; Herman, Z.; Fong, C.H. Limonoid biosynthesis in the stem of Citrus limon. Phytochemistry 1988, 27, 115–118. [Google Scholar] [CrossRef]
  38. Gerolino, E.; Chierrito, T.; Santana-Filho, A.; Rorigues, S.; Gonçalves, R.; Oliveira, A. Evaluation of limonoid production in suspension cell culture of Citrus sinensis. Rev. Bras. Farmacogn. 2015, 56. [Google Scholar] [CrossRef] [Green Version]
  39. Berim, A.; Gang, D.R. Methoxylated flavones: Occurrence, importance, biosynthesis. Phytochem. Rev. 2016, 15, 363–390. [Google Scholar] [CrossRef]
  40. Abad-García, B.; Garmón-Lobato, S.; Berrueta, L.A.; Gallo, B.; Vicente, F. On line characterization of 58 phenolic compounds in Citrus fruit juices from Spanish cultivars by high-performance liquid chromatography with photodiode-array detection coupled to electrospray ionization triple quadrupole mass spectrometry. Talanta 2012, 99, 213–224. [Google Scholar] [CrossRef]
  41. Mou, J.; Zhang, Z.; Qiu, H.; Lu, Y.; Zhu, X.; Fan, Z.; Zhang, Q.; Ye, J.; Fernie, A.R.; Cheng, Y.; et al. Multiomics-based dissection of citrus flavonoid metabolism using a Citrus reticulata × Poncirus trifoliata population. Hortic. Res. 2021, 8, 56. [Google Scholar] [CrossRef] [PubMed]
  42. Zhao, C.; Wang, F.; Lian, Y.; Xiao, H.; Zheng, J. Biosynthesis of citrus flavonoids and their health effects. Crit. Rev. Food Sci. Nutr. 2020, 60, 566–583. [Google Scholar] [CrossRef]
  43. Frydman, A.; Weisshaus, O.; Bar-Peled, M.; Huhman, D.V.; Sumner, L.W.; Marin, F.R.; Lewinsohn, E.; Fluhr, R.; Gressel, J.; Eyal, Y. Citrus fruit bitter flavors: Isolation and functional characterization of the gene Cm1,2RhaT encoding a 1,2 rhamnosyltransferase, a key enzyme in the biosynthesis of the bitter flavonoids of citrus. Plant J. 2004, 40, 88–100. [Google Scholar] [CrossRef]
  44. Chung, S.W.; Yu, D.J.; Oh, H.D.; Ahn, J.H.; Huh, J.H.; Lee, H.J. Transcriptional regulation of abscisic acid biosynthesis and signal transduction, and anthocyanin biosynthesis in ‘Bluecrop’ highbush blueberry fruit during ripening. PLoS ONE 2019, 14, e0220015. [Google Scholar] [CrossRef] [Green Version]
  45. Monselise, S.P.; Goldschmidt, E.E. Alternate bearing in citrus and ways of control. In Proceedings of the International Society of Citriculture, Tokyo, Japan, 9–12 November 1981; pp. 239–242. [Google Scholar]
Figure 1. Structure and classification of polyphenols quantified in commercial varieties of mandarin leaves. Polymethoxyflavones, tangeretin (a) and nobiletin (b); hydroxycinnamic acids, caffeic acid (c) and chlorogenic acid (d); flavanones, naringenin (e), didymin (f), narirutin (g), eriocitrin (i), and hesperidin (k); flavone, luteolin (h); flavonol, quercitrin (j); and limonoids, limonin (m) and nomilin (n).
Figure 1. Structure and classification of polyphenols quantified in commercial varieties of mandarin leaves. Polymethoxyflavones, tangeretin (a) and nobiletin (b); hydroxycinnamic acids, caffeic acid (c) and chlorogenic acid (d); flavanones, naringenin (e), didymin (f), narirutin (g), eriocitrin (i), and hesperidin (k); flavone, luteolin (h); flavonol, quercitrin (j); and limonoids, limonin (m) and nomilin (n).
Agronomy 12 00121 g001
Figure 2. (a) PCA plot of the four varieties of mandarin leaf at a bearing stage ‘Ellendale’ (B-END), ‘Willowleaf’ (B-WL), ‘Page’ (B-PG), and ‘Satsuma’ (B-SAT); (b) VIP scores of the PLS-DA performed over the four varieties of bearing mandarin leaves.
Figure 2. (a) PCA plot of the four varieties of mandarin leaf at a bearing stage ‘Ellendale’ (B-END), ‘Willowleaf’ (B-WL), ‘Page’ (B-PG), and ‘Satsuma’ (B-SAT); (b) VIP scores of the PLS-DA performed over the four varieties of bearing mandarin leaves.
Agronomy 12 00121 g002
Figure 3. (a) PCA plot of the four varieties of mandarin leaf in their non-bearing stage: ‘Ellendale’ (NB-END), ‘Willowleaf’ (NB-WL), ‘Page’ (NB-PG), and ‘Satsuma’ (NB-SAT); (b) PLS-DA VIP scores.
Figure 3. (a) PCA plot of the four varieties of mandarin leaf in their non-bearing stage: ‘Ellendale’ (NB-END), ‘Willowleaf’ (NB-WL), ‘Page’ (NB-PG), and ‘Satsuma’ (NB-SAT); (b) PLS-DA VIP scores.
Agronomy 12 00121 g003
Figure 4. Visual representation of Wilcoxon test for determination of significant differences between bearing and non-bearing stages, for each compound and each variety. References: ‘Ellendale’ (END), ‘Page’ (PG), ‘Satsuma’ (SAT), ‘Willowleaf’ (WL), bearing stage (B), non-bearing stage (NB), Agronomy 12 00121 i001 means does not differ significantly using Wilcoxon test, p ≤ 0.05, and Agronomy 12 00121 i002 means differs significantly using Wilcoxon test, p ≤ 0.05.
Figure 4. Visual representation of Wilcoxon test for determination of significant differences between bearing and non-bearing stages, for each compound and each variety. References: ‘Ellendale’ (END), ‘Page’ (PG), ‘Satsuma’ (SAT), ‘Willowleaf’ (WL), bearing stage (B), non-bearing stage (NB), Agronomy 12 00121 i001 means does not differ significantly using Wilcoxon test, p ≤ 0.05, and Agronomy 12 00121 i002 means differs significantly using Wilcoxon test, p ≤ 0.05.
Agronomy 12 00121 g004
Figure 5. (a) Heatmap of bearing and non-bearing leaf samples of the four varieties of mandarin: ‘Ellendale’ (END), ‘Willowleaf’ (WL), ‘Page’ (PG), and ‘Satsuma’ (SAT); (b) VIP scores of the PLS-DA performed to identify the metabolites that differentiate between both bearing stages.
Figure 5. (a) Heatmap of bearing and non-bearing leaf samples of the four varieties of mandarin: ‘Ellendale’ (END), ‘Willowleaf’ (WL), ‘Page’ (PG), and ‘Satsuma’ (SAT); (b) VIP scores of the PLS-DA performed to identify the metabolites that differentiate between both bearing stages.
Agronomy 12 00121 g005
Figure 6. Biosynthetic pathway from naringenin to hesperidin through eriodictyol and hesperitin. The aglycone is then glycosidated with glucose and rhamnose. In a similar sequence, intact naringenin is glycosidated two times; first, to form prunin, which is later glycosidated by the 6-O-rhamnosyl transferase to form narirutin.
Figure 6. Biosynthetic pathway from naringenin to hesperidin through eriodictyol and hesperitin. The aglycone is then glycosidated with glucose and rhamnose. In a similar sequence, intact naringenin is glycosidated two times; first, to form prunin, which is later glycosidated by the 6-O-rhamnosyl transferase to form narirutin.
Agronomy 12 00121 g006
Table 1. Results of optimization parameters for each evaluated compound.
Table 1. Results of optimization parameters for each evaluated compound.
Analyte[M-H][M-H] ⟶ m/zDWDPEPCECXPRt (min)
Protocatechuic acid153108.8150−40−10−22−159.4
90.7150−40−10−42−1
Chlorogenic acid353.7190.8150−50−10−24−3111.1
84.7150−50−10−68−9
Caffeic acid179133.3150−45−10−40.7−21.112.7
88.7150−45−10−48−14.1
Eriocitrin 595.2286.6150−80−10−38−914.5
134.6150−80−10−90−19
Narirutin579.2150.5150−95−10−75−20.315.9
270.8150−95−10−39.9−8.4
Naringin 579.2270.8150−105−10−54−716.3
150.7150−122.03−10−75.03−21
Quercitrin 447.1300.5150−90−10−36−1116.3
270.8150−90−10−66−41
Hesperidin609.2300.8150−85−10−46−4916.5
285.8150−85−10−70−41
163.5150−85−10−88−25
Neohesperidin609.2300.7150−115−10−56−9.3216.9
285.7150−115−10−72.89−42.3
150.5150−115−10−100−21
Didymin 593.2284.8150−100−10−48−4519.8
164150−100−10−92−27
Poncirin 593.3284.8150−120−10−54−920.0
163.4150−120−10−92−23
Luteolin 285150.6150−85−10−41.6−20.920.5
132.9150−85−10−54−17
Quercetin 301150.8150−75−10−36−2120.7
178.7150−75−10−28−27
120.8150−75−10−44−16.3
Kaempferol 284.992.9150−80−10−51.4−14.123.1
116.88150−80−10−67−18.8
Naringenin271150.7150−60−10−28.87−23.1423.2
118.7150−60−10−40.9−16.7
106.6150−60−10−37.8−16.9
Analyte[M-H]+[M-H]+⟶ m/zDWDPEPCECXPRt (min)
Bergapten217.220215064.11030.910.925.8
174.315064.11041.46.1
90.115064.11057.93.1
Limonin 471.2425.3150711028.62226.7
161.1150711037.712.7
Nobiletin 403.1373.11506110411426.7
388.115061103116
Nomilin 515.2469.315044.91020.1811.427.9
161.115044.91035.111.2
Tangeretin372.973431507110392628.4
297.115071104710
211.1150711060.877.3
M, molecular weight; m/z, selected ions monitored; DW, dwell time in msec; DP, declustering potential; EP, entrance potential; CE, collision energy; CXP, collision cell exit potential; and Rt, retention time.
Table 2. Identified phytochemicals in bearing leaf of ‘Ellendale’ (B-END), ‘Willowleaf’ (B-WL), ‘Page’ (B-PG), and ‘Satsuma’ (B-SAT) and their concentration expressed in μg/g of dried leaf.
Table 2. Identified phytochemicals in bearing leaf of ‘Ellendale’ (B-END), ‘Willowleaf’ (B-WL), ‘Page’ (B-PG), and ‘Satsuma’ (B-SAT) and their concentration expressed in μg/g of dried leaf.
CompoundB-SATB-ENDB-WLB-PG
(μg/g) *
Chlorogenic Acid3.88 b ± 0.403.53 b ± 0.376.70 a ± 0.785.04 b ± 1.39
Caffeic Acid3.87 b ± 0.214.34 ab ± 0.235.38 a ± 1.055.18 a ± 0.51
Didymin10.57 b ± 1.1120.19 a ± 3.273.30 c ± 0.4321.17 a ± 3.76
Hesperidin3008.58 a ± 151.892512.39 b ± 182.612624.31 b ± 114.093152.43 a ± 298.13
Eriocitrin47.93 a ± 6.4213.90 c ± 1.1929.85 b ± 2.4457.24 a ± 8.02
Luteolin1.45 a ± 0.051.21 a ± 0.081.81 a ± 0.282.31 a ± 1.26
Naringenin0.38 b ± 0.020.57 a ± 0.010.34 b ± 0.010.63 a ± 0.15
Narirutin87.03 a ± 12.6850.08 b ± 7.5817.58 c ± 7.5898.20 a ± 14.35
Quercitrin2.29 b ± 0.072.43 b ± 0.132.98 ab ± 0.553.69 a ± 0.61
Nobiletin258.17 c ± 16.481121.36 a ± 34.94747.15 b ± 55.16208.17 c± 13.63
Tangeretin64.96 c ± 6.15731.28 a ± 30.15475.84 b ± 56.2434.34 c ± 5.01
Limonin0.45 c ± 0.040.88 a ± 0.080.72 b ± 0.120.76 ab ± 0.08
Nomilin0.42 b ± 0.030.62 a ± 0.080.57 ab ± 0.130.52 ab ± 0.04
* In each row, means (5 replicates) ± SD followed by the same letter do not differ significantly using significance t-test with Bonferroni correction, p ≤ 0.05.
Table 3. Identified phytochemicals in non-bearing leaves of ‘Ellendale’ (NB-END), ‘Willowleaf’ (NB-WL), ‘Page’ (NB-PG), and ‘Satsuma’ (NB-SAT) and their concentration expressed in μg/g of dried leaf.
Table 3. Identified phytochemicals in non-bearing leaves of ‘Ellendale’ (NB-END), ‘Willowleaf’ (NB-WL), ‘Page’ (NB-PG), and ‘Satsuma’ (NB-SAT) and their concentration expressed in μg/g of dried leaf.
CompoundNB-SATNB-ENDNB-WLNB-PG
(μg/g) *
Chlorogenic Acid126.02 a ± 6.4011.31 b ± 0.626.75 b ± 1.416.35 b ± 1.50
Caffeic Acid2.32 a ± 0.122.00 b ± 0.090.12 c ± 0.040.04 c ± 0.02
Didymin6.44 bc ± 0.4717.40 a ± 2.975.03 c ± 2.5614.07 ab ± 8.77
Hesperidin2591.54 b ± 55.632433.21 b ± 179.794224.72 a ± 450.984270.21 a ± 400.09
Eriocitrin132.39 a ± 15.0117.09 d ± 1.8239.62 c ± 5.3469.98 b ± 17.91
Luteolin3.16 c ± 0.602.13 c ± 0.152659.07 a ± 267.402023.38 b ± 121.10
Naringenin0.32 b ± 0.020.68 b ± 0.12509.61 a ± 100.14574.03 a ± 101.15
Narirutin97.45 c ± 8.0862.85 c ± 5.0410049.66 b ± 1005.6442198.36 a ± 3700.22
Quercitrin2.80 b ± 0.212.44 b ± 0.16350.39 a ± 122.31362.45 a ± 50.83
Nobiletin290.48 e ± 24.071257.23 a ± 87.181028.88 b ± 151.80228.71 c ± 12.86
Tangeretin99.89 c ± 10.05517.97 a ± 78.41347.23 b ± 82.5933.27 c ± 5.98
Limonin3.40 a ± 0.440.77 c ± 0.051.30 b ± 0.091.28 b ± 0.03
Nomilin0.6 ab ± 0.010.54 b ± 0.050.19 a ± 0.140.66 ab ± 0.03
* In each row, means (5 replicates) ± SD followed by the same letter do not differ significantly using significance t-test with Bonferroni correction, p ≤ 0.05.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Rodríguez Ceraolo, C.; Vázquez, V.; Migues, I.; Cesio, M.V.; Rivas, F.; Heinzen, H. Flavonoids and Limonoids Profiles Variation in Leaves from Mandarin Cultivars and Its Relationship with Alternate Bearing. Agronomy 2022, 12, 121. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy12010121

AMA Style

Rodríguez Ceraolo C, Vázquez V, Migues I, Cesio MV, Rivas F, Heinzen H. Flavonoids and Limonoids Profiles Variation in Leaves from Mandarin Cultivars and Its Relationship with Alternate Bearing. Agronomy. 2022; 12(1):121. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy12010121

Chicago/Turabian Style

Rodríguez Ceraolo, Cecilia, Valeria Vázquez, Ignacio Migues, María Verónica Cesio, Fernando Rivas, and Horacio Heinzen. 2022. "Flavonoids and Limonoids Profiles Variation in Leaves from Mandarin Cultivars and Its Relationship with Alternate Bearing" Agronomy 12, no. 1: 121. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy12010121

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop