Systematic Chemical Analysis Approach Reveals Superior Antioxidant Capacity via the Synergistic Effect of Flavonoid Compounds in Red Vegetative Tissues

The flavonoid system comprises an abundance of compounds with multiple functions; however, their potential synergism in antioxidant function remains unclear. We established an approach using ever-red (RL) and ever-green leaves (GL) of crabapple cultivars during their development to determine interrelationships among flavonoid compounds. RL scored significantly better than GL in terms of the type, composition, and diversity of flavonoids than GL. Principal component analysis predicted flavonoids in RL to have positive interaction effects, and the total antioxidant capacity was significantly higher than the sum of antioxidant capacities of the individual compounds. This synergy was verified by the high antioxidant capacity in rat serum after feeding on red leaves. Our findings suggest that the synergistic effect is a result of the high transcription levels regulated by McMYBs in RL. In summary, individual flavonoids cooperate in a flavonoid system, thus producing a synergistic antioxidant effect, and the approach used herein can provide insights into the roles of flavonoids and other compounds in future studies.

Introduction

Flavonoids are the most abundant class of polyphenol compounds in the human diet and are widely found in grains, vegetables, fruits, and their processed products (wine, juice) and health foods (Bondonno et al., 2012; Czank et al., 2013). The basic flavonoid structure consists of A and B benzene rings, an oxygen-containing pyran C ring and a polyhydroxyl structure, which enables a single compound to have strong and diverse biological functions, including antioxidant, anti-aging, and anticarcinogenic capacities (Tsuji et al., 2013). For example, catechin from green tea (Chyu et al., 2004), phloridzin from Malus plants (Nair et al., 2014), and anthocyanins from berries (Bagchi et al., 2004) have been shown to reduce atherosclerotic lesions, to have anticancer effects, and to have anti-oxidant functions, respectively. Furthermore, various flavonoid monomers have been actively developed as drugs used in human disease control. For example, the quercetin monomer can regulate hyperthyroidism (Panda and Kar, 2007).

Thousands of flavonoid phytochemicals have been identified throughout the plant kingdom (Crozier et al., 2009). The main sub-classes of flavonoids are flavones, flavanones, isoflavones, flavonols and their glycosides, flavanonols, flavan-3-ols, and anthocyanidins, while the minor components include flavan-3,4-diols, biflavones, dihydrochalcones, and flavans (Pietta, 2000; Cushnie and Lamb, 2005; Crozier et al., 2009). Flavonoids cooperate to form a biochemical reaction system. When considering individual flavonoid compounds as drugs for disease treatment and stress release, it is important to note that although these flavonoid compounds are present in most dietary sources and are beneficial for human health, how exactly the flavonoid system functions and the effects of intake of these compounds in animals remain poorly understood.

The flavonoid system has several properties. First, there is a top-down metabolic flux from flavones, via flavonols, to anthocyanidins in the phenylpropanoid biosynthesis pathway (Shirley, 1996). Second, this system contains several important metabolic branches, from dihydrochalcones to phloridzin and its glycoside, dihydroflavonols to flavonols (e.g., quercetin, kaempferol, rutin, and their ligands), leucoanthocyanidins to flavan-3-ols (e.g., catechinic acid, epicatechin), and anthocyanidins to anthocyanins (Falcone Ferreyra et al., 2012; Henry-Kirk et al., 2012). Third, there are complex interrelations among compounds, including positive and negative feedback regulation in metabolic processes (Yin et al., 2012). Fourth, this system changes constantly in response to environmental and developmental cues. Fifth, most compounds have powerful biological functions, such as antioxidant, anti-aging, and anti-cancer capacities (Czank et al., 2013; Sáez-Ayala et al., 2015; Srinivasan et al., 2016). Moreover, different plants, organs, and tissues have different compositions and contents. For example, flavonoid compounds are significantly more abundant in red (fruit peels, leaves, berries) than in green tissues (Revilla and Ryan, 2000; Tsao et al., 2003; Tian et al., 2015). These contents are associated with the transcriptional regulation of structural genes in the flavonoid biosynthesis pathway. Flavonoids are synthesized from three molecules of malonyl CoA and one molecule of p-coumaroyl CoA via the phenylpropanoid pathway. Anthocyanin accumulation in ever-red leaf crabapple (Malus spp. “Royalty”) is regulated by the transcription factor McMYB10, which activates McF3′H and other structural genes involved in anthocyanin biosynthesis (Tian et al., 2015). In white Chilean strawberry (Fragaria chiloensis ssp. chiloensis f. chiloensis) fruit, down regulation of FcMYB1 expression results in the up-regulation of ANS and strong repression of anthocyanidin reductase (ANR) and leucoanthocyanidin reductase (LAR) transcript accumulation (Salvatierra et al., 2013). In developing seedlings of Arabidopsis thaliana, AtMYB11, AtMYB12, and AtMYB111 control flavonol accumulation through transcriptional regulation of their target genes (Stracke et al., 2007; Pandey et al., 2014).

Free radicals, as strong oxidizers, can cause degenerative diseases and stress, such as cancer, aging, and apoptosis. The antioxidant' properties of flavonoids have been well-described (Candeias et al., 1993; Roleira et al., 2015). The mechanisms underlying antioxidant activities of flavonoids include (1) suppression and scavenging of reactive oxygen species by inhibiting enzymes or chelating trace elements involved in free-radical production; (2) inhibition of lipid peroxidation by decreasing malondialdehyde (MDA) levels; and (3) up-regulation or protection of antioxidant defenses by stimulating the activities of antioxidant enzymes, such as superoxide dismutase (SOD), catalase (CAT), peroxidase (POD), ascorbate peroxidase (APX), and glutathione peroxidase (GPX) (Symons and Gutteridge, 1998). Procyanidins have been found to decrease MDA levels, maintain a high reduced/oxidized glutathione ratio, and increase CAT/SOD and GPX/SOD ratios as well as glutathione reductase and glutathione transferase activities (Roig et al., 2002). The actual antioxidant capacity, depends on the concentration of flavonoid compounds and the positive interaction/synergistic effects among the antioxidants in a mixture (Hernandez et al., 2009). However, synergistic interactions of flavonoid system compounds with respect to antioxidant capacity remain unclear.

To explore such synergistic effects, we focused on the Malus crabapple, which belongs to the family Rosaceae, Malus Mill. Malus crabapple originated in Europe and have been cultivated worldwide for 240 years, and the genome of crabapple is highly similar to that of Malus domestica (Zhang, 2007; Velasco et al., 2010). In this study, for comparison, we selected a Malus cultivar, M. spp. “Royalty,” with typical ever-red leaves, and M. spp. “Flame,” with ever-green traits and leaves similar to those of M. domestica. In our previous AFLP-based genetic study, we found that both cultivars are closely related (Zhang, 2007). The highly abundant flavonoids in crabapple leaves have antitumor activity and have been suggested to represent an excellent source of antioxidants for use as food additives (Qin et al., 2015). Our preliminary study revealed that the diversity and contents of phenolic compounds in “Royalty” leaves are greater than those in M. domestica leaves, and Malus crabapple represents a useful model system as one of the most economically important ornamental apple germplasm resources owing to the high flavonoid levels in leaves, flowers, and fruits (Tian et al., 2011). Here, we developed a systematic approach to evaluate the potential synergistic antioxidant effect of flavonoids in crabapple, including: (1) identification of individual compound classes; (2) quantification of the temporal dynamics of all compounds at various leaf developmental stages; (3) estimation of diversity properties of the flavonoid system; (4) clustering of leaves based on diversity indices and selection of optimal leaves for following measurements and/or applications; (5) analysis of the effects, contributions, and interrelations of individual compounds to identify predominant compounds using principal component analysis (PCA); (6) exploration of the molecular mechanism of flavonoid system establishment; (7) detection of the antioxidant capacity of flavonoid compounds in on-line and off-line systems to compare differences in antioxidant between total antioxidant capacity (T-AOC) and the sum of capacities of individual compounds; (8) verification of the antioxidant capacity in animal experiments.

Results

Types and Compositions of Flavonoid Compounds

To identify the flavonoid compounds in crabapple red leaves (RL) and green leaves (GL), we subjected leaf extracts to high-performance liquid chromatography coupled with tandem mass spectrometry (HPLC-MS) and compared the findings to the standards and reported references. In total, 19 compounds were identified (Text S1), including 16 flavonoid compounds, one phenolic acid and two terpenoids (Table 1 and Supplementary Figure S1). RL contained 17 of the 19 compounds, and GL contained 16 of the 19 compounds, with 14 compounds in common. P2, P4, and P7 were detected only in RL, while P6 and P19 were detected only in GL.

TABLE 1

Table 1. Main classes of compounds identified by LC-UV-ESI-MS² in crabapple leaves.

The concentrations of total flavonoid compounds were significantly higher in RL than in GL at the S1–S7 and S10 leaf nodes (representing leaf developmental stages, Figure 1B). The average contents of flavonoid compounds P1, P2, P4, P7, P9, P11, P13, P15, P16, and P18 were significantly higher in RL than in GL (Figure 1C) during all leaf developmental stages (Figure 1C). The dynamics of the changes in the individual flavonoid compounds differed significantly between RL and GL (Supplementary Figure S2). The profiles of these compounds may reflect differences in the distribution and composition of flavonoid metabolites in the leaves of the two cultivars. Luteolin-5-O-rutinoside, a flavone produced from the direct precursor naringenin, by flavone synthase (FNS), was detected only in GL.

FIGURE 1

Figure 1. Analysis of phenotypes and compound compositions detected by HPLC in Malus ever-red (RL) and ever-green leaves (GL). (A) Phenotypes of leaves at 10 developmental stages of RL and GL. Scale bar represents 1.0 cm. (B) HPLC analysis of total detected compounds at the 10 developmental stages in RL and GL. The data (mean ± SD) are the sum for each detected compound. S1–S10, the leaves at the 1st to 10th nodes on the annual shoot from the tip toward the base, represent the 10 leaf developmental stages. (C) HPLC analysis of each compound detected in the leaves at the 10 developmental stages of RL and GL. The data show the average values (mean ± SD) of the compounds in S1–S10. P1–P19 is the same as in Table 1.

Properties of Flavonoid Systems

To determine the diversity of the flavonoid compounds in the leaves, we analyzed the diversity (H), dominance (D), and evenness (E) indices. H differed significantly between RL and GL, with high diversity during early leaf developmental stages. D was significantly higher for RL than for GL, except for S4. E also differed significantly between RL and GL in the S1, S2, S5, S6, S8, S9 developmental stages (Table 2). The three indices for S1, S2, S6, and S8 leaves showed significant differences between RL and GL. These results indicated that the diversity of the complex flavonoid systems in RL and GL throughout leaf development differs.

TABLE 2

Table 2. The diversity properties of flavonoid compounds in crabapple leaves.

Based on the above indices, we conducted CA to identify groups with similar temporal diversity of the flavonoid composition. When the similarity coefficient r ≈ 1 (Chauhan and Vaish, 2012), these temporal patterns for flavonoid composition in RL and GL could be divided into four groups (Figures 2A,B). In RL, the groups were S1–S6, S7, S10, and S8–S9; in GL, the groups were S1–S3, S4–S9, S7, and S10. These results suggest that the leaves at different nodes in each group have obvious differences in the spatiotemporal pattern of flavonoid diversity. S1–S6 in RL and S4–S9 in GL were the groups encompassing the most nodes, suggesting the changes in flavonoids diversity followed a similar spatiotemporal pattern among these leaf groups. The S4, S5, and S6 leaves were comparable in both RL and GL, which may inform leaf material selection in future studies.

FIGURE 2

Figure 2. Clustering analysis of the diversity characteristics of flavonoid compounds in RL (A) and GL (B) and principal component analysis of various types of compounds in the flavonoid systems of RL (C) and GL (D); r represents the similarity coefficient. P1-P19 in (C,D) are the same as in Table 1.

To determine the contribution rates (CRs) [usually used to analyze the effects of various factors on system functions; (Baerenfaller et al., 2012)] of the individual flavonoids to the function of the flavonoid system and the differences between the two cultivars, we performed PCA based on the combined data sets for each flavonoid compound plus the 10 developmental periods. For RL (Figure 2C), the first and second principal components (PCs) accounted for 87.65% of the total variance, where PC₁ = 59.59% and PC₂ = 28.06%. In PC₁, the CRs of P1, P4, P10, P15, and P18 were the greatest, and in PC₂, the CRs of P11, P12, P14, P16, and P17 were the greatest. For GL (Figure 2D), the first and second PCs accounted for 73.72% of the total variance, where PC₁ = 41.95% and PC₂ = 31.77%. In PC₁, the CRs of P1, P3, P6, P9, P10, and P14 were the greatest, and in PC₂, the CRs of P13, P16, and P17 were the greatest. These results indicate that the numbers of flavonoid compounds in PC₁ and PC₂ differ between RL and GL and that different compounds have different contribution to the respective flavonoid systems.

Flavonoid Biosynthesis Related Gene Profiling

To explore whether the differences in the levels of flavonoid compounds can be explained by different biosynthesis mechanisms in RL and GL, we selected the leaves at S6 (based on the above diversity index and CA results) to compare the expression levels of the main flavonoids biosynthetic genes using real-time PCR. The expression of McCHS was two-fold higher in RL than in GL (Figure 3), resulting in the accumulation of dihydrochalcone compounds, such as phloridzin dihydrates (Supplementary Figure S2). The expression levels of McCHI, McF3H, McF3′H, and McFLS were higher in RL than in GL (Figure 3). McDFR, McANS, and McUFGT showed higher expression in RL, leading to strong accumulation of cyanidin-3-O-galactoside (Supplementary Figure S2). McMYB10, encoding an R2R3 transcription factor (TF) involved in activating anthocyanin biosynthesis (Tian et al., 2015), showed significantly higher expression in RL than in GL. However, McMYB16, which is homologous to FcMYB1, a repressor of anthocyanin biosynthesis in strawberry fruit (Salvatierra et al., 2013), showed lower expression in RL than in GL. McMYB4, homologous to AtMYB12, which is involved in flavonol biosynthesis in Arabidopsis (Stracke et al., 2007), displayed a dramatically higher expression in RL than in GL.

FIGURE 3

Figure 3. Analysis of the relative expression levels of the main structural genes and the transcription factor MYBs related to the flavonoid biosynthesis pathway in RL and GL. CHS, chalcone synthase; CHI, chalcone isomerase; F3H, flavanone 3β-hydroxylase; F3′H, flavonoid 3′-monooxygenase; DFR, dihydroflavonol 4-reductase; ANS, anthocyanidin synthase; UFGT, flavonoid 3-O-glycosyltransferase; FLS, flavonol synthase. McMYB4, McMYB10, and McMYB16 are R2R3 transcription factors in crabapple. The data shown are the mean ± SD. ** and * indicate significance at P < 0.01 and P < 0.05 by t-test, respectively.

In Vitro Antioxidant Capacity

To compare the antioxidant capacity of RL and GL, we measured SOD and CAT activities, MDA content and T-AOC. SOD activity was significantly higher and MDA content significantly lower in RL than in GL, while CAT activity showed no significant difference between RL and GL. T-AOC was significantly higher in RL than in GL (Figures 4A–D).

FIGURE 4

Figure 4. Antioxidant capacity of the detected flavonoid compounds in the extracts of RL and GL determined by HPLC-ABTS in vitro. (A,B) Superoxide dismutase (SOD) and catalase (CAT) activities. (C,D) Malondialdehyde (MDA) and total antioxidant capacity (T-AOC) levels. (E,F) Antioxidant activity (negative peaks) of each detected flavonoid compound (positive peaks) in RL and GL, respectively. (G,H) Total antioxidant activities (negative peaks) of the detected flavonoid compounds (positive peaks) in RL and GL, respectively. For each sample, the extract solution was diluted five-fold in RL and three-fold in GL. (I) Individual flavonoid compounds related to the Vc antioxidant coefficient in RL and GL. (J) Comparison of T-AOC (HPLC-ABTS) levels and sum of the antioxidant activities of the individual compounds (SAIC) (negative peak area) in RL and GL. The UV profile was monitored at 350 nm on the chromatogram, while the scavenging activity profile was determined using HPLC-ABTS and monitored at 517 nm (negative peaks). The values shown are the mean ± SD (N = 10/group). ** and * indicate significance at P < 0.01 and P < 0.05 by t-test, respectively.

To further explore the differences in antioxidant capacity for each flavonoid compound, we improved the HPLC with radical scavenging methods (HPLC-ABTS) system for crabapple leaf extracts based on screening results from the on-line HPLC-ABTS system (Figures 4E–H). The antioxidant activities of P3 (phenolic acid), P5, P9, P11, P13, P16, and P18 were significantly higher in RL than GL. P2 and P19 were detected only in RL and GL, respectively, and both had obvious antioxidant capacities (Supplementary Table S1). In addition, the antioxidant activities relative to the ascorbic acid (Vc) of P3, P5, P9, P11, P13, P15, P16, and P18 were higher in RL than in GL, while P12 activity was lower in RL than in GL. The differences in the activities of P3, P5, P9, P12, and P15 were statistically significant (Figure 4I). T-AOC (HPC-ABTS) was significantly higher than the sum of the antioxidant activities of the individual compounds (SAIC) in both RL and GL, T-AOC (HPC-ABTS) was higher than the sum of the antioxidant activities of the individual compounds (SAIC) in both RL and GL. T-AOC (HPC-ABTS) was three-fold higher than SAIC in RL and two-fold higher than T-AOC (HPC-ABTS) in GL. T-AOC (HPLC-ABTS) was also nearly three-fold greater than SAIC in GL (Figure 4J). These results showed that the high antioxidant capacity of RL was associated with the high concentrations of flavonoids, and from the high antioxidant activities of individual flavonoids and flavonoid synergistic effects.

Antioxidant Activity in the Blood of Rats Fed Leaf Extracts

With the objective of introducing and producing a new nutrient, we conducted an animal model experiment to demonstrate that the crude extract of RL has powerful antioxidant capacity. The results indicated that the SOD and CAT activities and T-AOC remained significantly higher and the MDA content remained significantly lower in the serum of rat fed a diet with RL extract (FRL group) than in the serum of the rat fed a diet with GL extract (FGL group) or standard diet (CK group) (Figure 5). We also measured blood glucose (GLU), cholesterol (TC), and triglyceride (TG) levels in rat sera. Serum TC and TG levels were lower in the FRL group than in the FGL group (Figure 5). GLU, TC, and TG levels were negatively correlated with T-AOC (Figures 5E–G). These results indicated that consumption of RL extract might improve antioxidant capacity, decreasing pathological indices in rat blood.

FIGURE 5

Figure 5. Antioxidant capacity in rat serum. (A–G) CK, the normal control group of rats receiving the basic diet and distilled water; FRL, the treatment group of rats receiving the basic diet and red leaf extract solution (50 mg/kg/day, dissolved in distilled water); FGL, the treatment group of rats receiving the basic diet and green leaf extract solution (50 mg/kg/day, dissolved in distilled water). GLU, blood glucose; TC, cholesterol; TG, triglyceride. (H–J) Regression analysis of T-AOC to GLU (H), TC (I), and TG (J) in rat serum. The values are expressed as the mean ± SD (N = 10/group). Lower-case letters indicate significance at P < 0.05 by Duncan's new multiple range test.

Discussion

Flavonoid compounds play important roles in various aspects of plant growth and development and have positive effects on human quality of life (Crozier et al., 2009). Their metabolism and function in plants are affected by not only environmental, but also inherent factors, such as transcriptional regulation and the interactions of the compounds with each other and their surroundings (Winkel-Shirley, 2001). Flavonoids are involved in defense strategies and signaling, particularly, in interactions between plants and their environment. Each class of flavonoid has been shown to have biological activity, epidemiological studies have revealed that flavonoids have antioxidant activities toward cancer, cardiovascular, and neurodegenerative diseases, or an anti-aging effect that is useful for application in cosmetic products (Grotewold, 2006; Boudet, 2007). As a chemical reaction system, the compositions and structures of flavonoid compounds vary greatly, and the interrelation of compounds is complex and diverse due to the metabolic flux and branches catalyzed by the enzyme system (Jorgensen et al., 2005). Further, the flavonoid compositions and contents of different plants vary depending on growth and developmental stages. This complex regulation includes top-down metabolic flux, bottom-up feedback regulation, and interactions among compounds, thus, different of flavonoids coordinate with each other in metabolism to form a flavonoids system. The flavonoid system is established through the interactions and interdependence of multiple compounds possessing special biological functions as well as the external plant environment. Understanding the biological function of the flavonoid system from a metabolic perspective and, importantly, explaining the effects of their intake are research goals that have been attracting increasing attention. Here, we attempted to establish a systematic, multistep approach to address these questions.

PCA implied that the synergistic effect of the flavonoid system played an important role in the antioxidant capacity of crabapple leaves. The percentage of PC₁, the number of compounds distributed in the first quadrant, and the correlation among the first five compounds with high contribution rates to PC₁ were all higher in RL than in GL. This result supported the conclusion that a significant cooperative relationship was involved in the flavonoid system in RL, whereas this relationship was inferior in GL. Furthermore, analysis of the diversity, evenness and dominance indices, combined with CA, revealed distinguishing characteristics of the flavonoid system at the leaf developmental stages of RL and GL. However, the leaves at certain stages, such as S6, stood out among the leaves from all other10 stages, displaying optimal flavonoid system characteristics in the two cultivars. This finding may provide guidance regarding the selection of optimal leaf materials for further study and application.

When tracing the origin of the synergistic antioxidant effect of the flavonoid system, the flavonoid biosynthesis pathway should be considered. A top-down flavonoid metabolic flux with positive or feedback regulation on the process exists (Yin et al., 2012), where flavonoid accumulation is controlled by structural genes (such as CHS, F3H, and DFR) and transcription factors (such as R2R3-MYBs) (Nakabayashi et al., 2014). McMYB10, an R2R3-MYB TF, is an activator of anthocyanin biosynthesis via regulating the genes F3'H, DFR, and ANS in red tissues (Tian et al., 2015), while McMYB16 may be a repressor of red formation, as downregulation of the homologous FcMYB1 in strawberry fruits induced the up-regulation of ANS and a strong repression of ANR and LAR transcript accumulation, resulting in increased anthocyanins and undetectable levels of flavan-3-ols (Salvatierra et al., 2013). McMYB4, homologous to AtMYB111 of Arabidopsis thaliana, may participate in flavonol biosynthesis in crabapple. AtMYB111 was reported to control flavonol accumulation through interaction with target genes, such as FLS and UGT91A1 and UGT84A1 (Pandey et al., 2014). In addition, slight UFGT expression was detected in GL, which differed from a report in white grapes, whose skin showed loss of UFGT expression (Boss et al., 1996) but was similar to the results of UFGT expression in white apple flesh and skin as well as red and green leaves (Espley et al., 2007). These results indicated that UFGT may be involved in the biosynthesis of other flavonoid glycosides, such as quercetin 3-O-β-D-glucuronide and flavone 3-O-glucoside (Xiao et al., 2014). Interestingly, luteolin-5-O-rutinoside was detected only in GL in this study. The reason might be the high FNS activity in GL and absence of FNS activity in RL, similar to the finding that red petals of Gerbera Hybrida show no FNS II expression while white petals do (Martens and Mithofer, 2005). The molecular mechanism by which UFGT regulates flavonoid glycoside biosynthesis and FNS regulates flavones (e.g., luteolin-5-O-rutinoside in GL) should be explored further.

The on-line HPLC-ABTS system provided fast, reliable and high-efficiency detection of the antioxidant capacity of individual compounds, as demonstrated by Shi et al. (2009). In the present study, we identified 12 negative peaks for 19 types of flavonoid compounds and demonstrated that T-AOC was higher than that indicated by the measurements of each compound (Supplementary Table S1), suggesting a synergistic effect among the compounds. However, using online HPLC-ABTS, we were unable to determine the antioxidant capacity of some compounds because no negative peaks appeared in P1, P4, P6, P7, P8, P10, and P18. The lack of negative peaks may have been due to the low content or low antioxidant capacity of these compounds. Alternatively, the corresponding negative peaks may not have been detected using HPLC-ABTS if the online reaction time was insufficient to allow reaction between ABTS+ and the antioxidants in the test system, such as the phloridzin (P18). Therefore, we transformed the phloridzin values determined by ultraviolet spectrophotometry to the values corresponding to the online analysis (Supplementary Table S1). The results indicated that this online system has potential for analyzing total and individual antioxidant capacities of functional compounds, but requires further improvement.

GLU, TC, and TG were lower, while SOD, CAT, and T-AOC activities were significantly higher, and MDA content was lower in the sera of treated rats than in those of non-treated control rats. This finding indicated that the enhancement of the antioxidant capacity helped reduce GLU, TC, and TG in rats fed RL extracts owing to some extent to the synergistic protective effect of total flavonoids (Lewis et al., 2013). Based on previous reports on the health benefits of flavonoids, these compounds should reduce fat accumulation and facilitate GLU movement as they participate in scavenging reactive oxygen species, inhibiting lipid peroxidation (which generates MDA), and promoting antioxidant defenses through the stimulation of antioxidative enzymes (Sugiura et al., 2012).

In view of the synergistic, it has been previously reported that, propolis induces a high bactericidal activity when consumed as a whole, while isolated compounds, shows low activity due to the lack of a stable and reproducible chemical composition. (Scheller et al., 1999; Bachiega et al., 2012). In addition, when foods are consumed together, the T-AOC may be modified via synergistic, additive, or antagonistic interactions among the constituents (Wang et al., 2011). Strong support for this conclusion was provided by our findings, which revealed that the T-AOC was significantly higher than the SAIC, with the synergistic effects being stronger in RL than in GL. Meanwhile, the antioxidant capacities of most of the individual compounds in RL were superior to those in GL. These findings could be explained by the synergy being a result of both interactive effects among and individual effects of the flavonoids compounds.

The systemic approach we established for comparing differences in the flavonoid system between cultivars, which describes the combined principles for deep analysis of results, is a novel approach compared to systems theory (Bertalanffy, 1968). Based on our findings, we recommend the application of this approach for other compounds in plant secondary metabolism and even in the areas of nutrient quality composition analysis in food and quantitative trait analysis in plant breeding.

Materials and Methods

Plant Materials

Leaves were collected from 4-year-old Malus cv. “Royalty” and “Flame” trees; grafted on Malus hupehensis and planted in the Crabapple Germplasm Resource Garden of BUA, Beijing, China. On June 10, 2013, when the annual shoots stopped growing, we randomly selected 10–20 annual shoots with similar growth vigor from the south side of each plant and collected 10 leaves from the tip to the base nodes of each shoot. The leaves representing the 10 developmental stages, were labeled S1 to S10. The weight at S1 was 0.1–0.2 g and that at S10 was 0.8–1.0 g, and the width of the collected leaves at the same stage was similar for each cultivar. The experiment was performed using a randomized block design of two cultivars with three biological replicates, and five trees were used in each plot. After collection, 30–40 leaf discs (1 cm diameter, punched in the middle of each side of the main vein) were collected randomly from leaves at each stage; per plot, frozen in liquid nitrogen, and stored at −80°C for compound analysis and RNA extraction. The remaining non-perforated leaves were stored at −20°C for enzyme analysis and animal experiments.

Animals and Experimental Design

Weaning male Sprague-Dawley rats (7-week-old, 180–200 g) were purchased from Beijing Xinglong Experimental Animal Farms and kept five per cage under standard housing conditions (12/12 h light/dark cycle, lights on at 7:00 a.m., humidity 50–60%, and at 22 ± 2°C), with ad libitum access to food and water throughout the study period. One week after acclimatization to a basic diet (containing crude protein ≥ 20%, crude ash ≤ 8%, crude fiber ≤ 5%, 1.8% ≥ calcium ≥ 1%, 1.2% ≥ phosphorus ≥ 0.6%, lysine acuity ≥ 1.35%, salt = 1.35%; purchased from Beijing Keao Xieli Feed Company Limited), the rats were randomly divided into three groups of six rats per group: Group 1 (CK), the normal control group received basic diet and distilled water; Group 2 (FRL), received basic diet and RL extract (50 mg/kg/day, dissolved in distilled water); and Group 3 (FGL), received basic diet and GL extract (50 mg/kg/day, dissolved in distilled water).

At the end of the 5-week experimental period, after 12 h of starvation, all rats were anesthetized with intraperitoneal chloral hydrate (300–350 mg/kg) at the same time of the day to avoid circadian metabolic fluctuations. Blood samples were obtained from the carotid artery and centrifuged immediately at 800 × g for 15 min at −4°C to separate the serum samples (Jia et al., 2014). All study procedures performed on the animals were approved by the Animal Care and Use Ethics Committees of China Agricultural University, and all rats received humane care in accordance with the guidelines.

Preparation of Leaf Extracts for the Rat Diet and Extraction of Flavonoids

To prepare the diet, RL and GL at S6 were extracted. Samples were extracted using an 80% methanol solution and a material-to-solvent ratio of 1:10 under sonication at a frequency of 12,000 Hz at 45°C. Extraction was performed three times in total for 2.0, 1.5, and 1.0 h. Next, the solutions were filtered and concentrated at 40°C with a rotary evaporator under reduced pressure before being frozen at low temperature in a ModulyoD-230 freeze dryer (Thermo Fisher, New York, USA) to obtain a powder. Approximately 50 g of leaf extract powder per sample type (RL and GL) was added to the rats' diet.

For flavonoid extraction, 1.0 g of frozen RL and GL was powdered in liquid nitrogen using a mortar and pestle. Then, 3 ml of methanol/water/formic acid (80:19:1, v/v/v) was added to the sample, followed by 50 min of sonication at 12,000 Hz at 45°C to extract the flavonoid compounds. The compounds were filtered through a 0.22-μm membrane (Shanghai ANPEL, Shanghai, China) before HPLC-DAD and HPLC-MS analyses. We used mixed leaf samples from S1 to S10 for the identification of flavonoid types and for the detection of temporal changes in the leaf flavonoid composition during leaf development. We used leaf samples at S6 for further experimental analyses after identifying it as the optimal stage.

HPLC-ESI (±)-MS² Analysis and Identification of Flavonoids

For HPLC-ESI-MS analyses of the flavonoids, we followed the method of Yi et al. (2016) with modification. Experimental standards of cyanidin-3-O-galactoside, quercetin-3-O-glucoside, astilbin, phloridzin, rutin, quercetin-3-O-arabinoside, quercetin, and luteolin were purchased from Sigma-Aldrich (Steinheim, Germany); experimental standards of 4-O-coumaroyl quinic acid, procyanidin B2, (-)-epicatechin, and taxifolin-3-O-glucoside were purchased from Sigma-Aldrich (Poole, U.K.); data obtained for experimental standards and information regarding the identified compounds were combined to identify the flavonoids.

HPLC-DAD Analysis

An Agilent 1100 series HPLC system (Agilent Technologies, Wilmington, DE, USA) was used in our experiments. For chromatographic separation, eluent A was 0.5% aqueous formic acid, eluent B was 100% acetonitrile, and the elution gradient was as follows 90% A at 0 min, 90% A to 81% A from 0 to 40 min, 81% A to 60% A from 40 to 50 min, 90% A at 50.01 min, and 90% A at 60 min. The flow rate was 1.0 ml/min, and injection volume was 20 μl. The ZORBAX Eclipse XDB-C18 (5 μm, 150 × 4.6 mm, Agilent, USA) analytical column temperature was 27°C for all analyses. Chromatograms were acquired at 520 nm for anthocyanins and at 350 and 280 nm for all other flavonoids. DAD data were recorded from 200 to 600 nm.

Quantitative Analysis of Individual Flavonoid Compounds

We used quantitative analysis of multi-components by single-marker (QAMS) to calculate the contents of phytochemical products in the leaves (Zhou et al., 2008; Hou et al., 2011). We analyzed eight classes of flavonoid compounds as follows. For semi-quantification of the anthocyanins, a standard curve was drawn from 0 to 200 μg/g with a standard compound of cyanidin-3-O-galactoside at 20 μg/g. The equation of the calibration curve describing the results was Y₁ (equivalent content μg/g) = (0.003*A₁+0.3)/m (m ≤ 1.0 g, R² = 0.9999), where A₁ is the peak area at 520 nm. For the semi-quantification of flavonol and phenolic acid, a standard curve was drawn from 0 to 3,000 μg/g with a standard compound of quercetin at 300 μg/g. The equation for the calibration curve describing the results was Y₂ (equivalent content μg/g) = (0.02*A₂-10.6)/m (m ≤ 1.0 g, R² = 0.9996), where A₂ is the peak area at 350 nm. For semi-quantification of the dihydrochalcone, proanthocyanidin, dihydroflavonol, flavones, and flavan-3-ols, a standard curve was drawn from 0 to 2,000 μg/g with a standard compound of phloridzin at 200 μg/g. The equation of the calibration curve describing the results was Y₃ (equivalent content μg/g) = (0.0005*A₃-0.05)/m (m ≤ 1.0 g, R² = 0.9998), where A₃ is the peak area at 280 nm.

Antioxidant Activity and Other Physical Characteristics of Leaf Samples at S6 in Vivo and in Vitro

On-Line HPLC-ABTS Assay

HPLC coupled with the ABTS assay was performed using the method (Arnao et al., 2001), with some modifications. ABTS [2, 2′-azino-di-(3-ethyl-benzothiazoline-6-sulphonic acid)] was purchased from Sigma-Aldrich (St Louis, MO, USA). A stock solution containing 140 mM potassium persulphate and 7 mM ABTS was prepared and kept at room temperature in darkness for 12–16 h to stabilize the radicals. The radical reagent was prepared by diluting the stock solution with absolute ethyl alcohol to an absorbance of 0.70 ± 0.02 at 734 or 414 nm. The S6 leaf sample was extracted in the same way as in the HPLC analysis [1.0 g frozen leaves (RL and GL) were extracted with 3 ml of methanol/water/formic acid (80:19:1, v/v/v)]. The HPLC-ABTS system used consists of a Waters e2695 HPLC separation system, post-column reaction system, Waters e2998 photodiode array detector (PAD), and Waters 2489 UV/Visible detector. We used the same ZORBAX Eclipse XDB-C18 (5 μm, 150 × 4.6 mm, Agilent, USA) analytical column in an Agilent 1100 HPLC system to separate the compounds from the samples. The HPLC separation conditions were the same as above to maintain the same retention time. The HPLC eluates from the column arrived at a T-junction, where the ABTS reagent was added, at 0.4 ml/min by a Waters Reagent Pump (Waters Corporation, Massachusetts, USA). After the eluates were mixed with ABTS reagent in a reaction coil (15 m, 0.25 mm i.d. PEEK tubing), a Waters temperature control module (Waters Corporation, Massachusetts, USA) maintained the temperature at 40°C. The negative peaks were measured by DAD at 517 nm. We replaced the Milli-Q water with ABTS⁺ as the control. The data were analyzed using the Waters ChemStation Software.

Analysis of the T-AOC Values of Mixed Compound Samples by an On-Line HPLC-ABTS System

In the on-line HPLC-ABT system, we used a Waters X-Terra MS C18 guard column (to hold the samples without separation of the compounds) instead of the ZORBAX Eclipse XDB-C18 (5 μm, 150 × 4.6 mm, Agilent, USA) analytical column to analyze the T-AOC (HPLC-ABTS) of the mixed compound samples. For a given sample, we diluted the extract solution five-fold for RL and three-fold for GL because the antioxidant activities of the undiluted extract solution from 1.0 g of leaf samples appeared as a flat negative peak using the UV/Visible chromatographic detector.

Detection of Antioxidant Capacity Relative to that of Ascorbic acid (Vc) of Phloridzin and Other Compounds

Phloridzin is prevalent in Malus plants and has been reported to have antioxidant activity both in vitro and in vivo (Jugde et al., 2008; Li et al., 2014). However, the on-line HPLC-ABTS system has certain shortcomings or limitations. On one hand, the on-line post-column reaction time was set shortly (most < 1 min) that might cause false-negative results for antioxidants with slow kinetics of reactions (Lu et al., 2017). For example, phloridzin as an excellent antioxidant, may stay in the reaction coil for longer than the limited time in the system, preventing its detection, its negative peak could not be detected in on-line ABTS analysis. On the other hand, the post-column reaction system consisted of a post-column reaction model (peak reaction coil) in which the antioxidants reacted with the free radicals, and the reaction might be not complete in a small coil. Therefore, we used a Rapid ABTS method (Total Antioxidant Capacity Assay Kit, Nanjing Jianchen Bioengineering Institute, Nanjing, China) to detect the antioxidant capacities of phloridzin and ascorbic acid (Vc) by A414 ELISA in an off-line system. After a 10 mM Trolox standard solution was diluted to 0.15, 0.3, 0.6, 0.9, 1.2, and 1.5 mM, and a linear equation for the absorbance/concentration ratio was established, the absorbance values of 1.0 mM phloridzin and of Vc were each determined relative to the Trolox-equivalent antioxidant capacity. We detected the concentrations of the 1.0 mM Vc sample negative peak area using the HPLC-ABTS system. According to the method described above, the values for other compounds relative to the Vc antioxidant values were calculated by the formula: relative to Vc = (each compound negative peak area/each compound positive peak area)/(Vc negative peak area/ Vc positive peak area).

Determination of SOD and, CAT Activities, MDA, and T-AOC

All enzyme extractions were conducted at 4°C. Each sample (0.1 g fresh weight, at S6) was thoroughly homogenized in 1.0 ml of 50 mM phosphate buffer (pH 7.8) containing 0.1 mM ethylene diaminetetraacetic acid (EDTA) and 0.05 g of quartz for homogenization. The homogenates for SOD, CAT and MDA were centrifuged at 8,000 × g for 10 min; and the homogenate for T-AOC was centrifuged at 10,000 × g for 10 min. The crude supernatant was used to determine the SOD and CAT activities, total antioxidant capacity (T-AOC) was measured using the ferric reducing activity of plasma (FRAP) assay, and MDA concentrations in each tissue, using commercial kits (Product Codes: FY1, FY2, EY2-1, FG5, respectively; Suzhou Comin Biotechnology Co., Ltd, Jiangsu, China). For measuring these four parameters in animal serum samples (centrifuged immediately at 800 × g for 15 min at 4°C from blood samples), the above kits were also used.

Determination of Glucose (GLU), TC, and TG Levels in Vivo

The serum level of GLU was measured according to a standard procedure using a Vet Test 8008 Automatic Biochemical Analyzer (IDEXX, USA). TC and TG were measured using a Hitachi 7020 Automatic Analyzer (Hitachi, Tokyo, Japan).

RNA Extraction and Real-Time PCR Analysis

Total RNA was extracted from leaves at S6 and used and real-time PCR analysis following by Yi et al. (2016). The primers used are listed in Supplementary Table S2. PCRs were carried out for 40 cycles of 95°C for 15 s and 59°C for 20 s and were run with three biological replicates and two technical replicates. Gene expression levels were calculated using the 2^(−ΔΔCt) method. Malus 18S ribosomal RNA gene (GenBank accession number: DQ341382) was used as a reference gene because of its consistent transcript level among “Royalty” and “Flame” leaves.

Statistical Analysis

For diversity analysis of flavonoid systems in both leaf types, we created matrices of the values of the measured content of each type of flavonoid compound at different developmental stages for RL and GL as well as the calculated Shannon-Wiener diversity index, the evenness index and Simpson's index according to the methods described by Magurran (2003) and Morris et al. (2014). The Shannon-Wiener diversity index (H) was applied to measure the diversity of the flavonoid systems: H' = $-{\displaystyle {\sum }_{\text{i}}^{\text{s}}{\text{P}}_{\text{i}}\ln {\text{P}}_{\text{i}}},\text{i}=1,2,3,\dots \text{s},{\text{p}}_{\text{i}}={\text{n}}_{\text{i}}/\text{N},$ where n_i is the total amount of the ith type of flavonoid compound, N is the total amount of all types of flavonoid compounds, and s is the number of types of flavonoid compounds in the sample. The Simpson index (D) was used to calculate the dominance index for each type of flavonoid compound: $D=-{\displaystyle {\sum }_{\text{i=1}}^{\text{s}}\frac{{\text{n}}_{\text{i}}\left({\text{n}}_{\text{i}}-\text{1}\right)}{\text{N }\left(\text{N}-\text{1}\right)}}$ where n_i is the amount of the ith type of flavonoid compound, and N is the total amount of all types of flavonoid compounds. The evenness (E) of the types of flavonoid systems was measured using the formula E = H'/ln S, where H' is the Shannon-Wiener diversity index, and S is the amount of all types of flavonoid compounds.

Clustering analysis (CA): CA is the most similar objects are grouped first, as the similarity decreases, all subgroups are fused into a single cluster (Simao et al., 2010). We employed system CA using the diversity index (H), dominance (D), and evenness index (E) as the indices using the SPSS 17.0 software (SPSS Inc., Chicago, IL, USA). We obtained clustering results for diversity in the flavonoid system of the leaves at 10 developmental stages using the Euclidean shortest distance method.

Principal component analysis (PCA): PCA convert a set of correlated variables into a set of values of linearly uncorrelated variables (Patras et al., 2011). The distributions of the compounds in RL and GL were then drawn in PCA diagrams using SPSS 17.0. These diagrams revealed the percentage variations of PC₁ and PC₂ in flavonoid systems, the contribution rates of individual compounds to PC₁ and PC₂, and the interrelations of the individual flavonoid compounds.

The data are presented as the mean ± standard deviation (SD). Analysis of variance (ANOVA, SPSS 17.0 software, SPSS Inc., Chicago, IL, USA) of all values was used to assess differences in the means among samples (P < 0.05). Duncan's multiple analysis and Student's t-test were used to identify significant differences among groups (P < 0.05, P < 0.01). Graphs were prepared in Origin Pro 8.0 SR4 (Origin Lab, Northampton, MA, USA) and Microsoft Office PowerPoint 2007.

Author Contributions

YY: designed the study; XQ and YL: performed the research; SF: provided new methods and materials; XQ, YL, and ZP: wrote the manuscript; and all authors approved the revisions.

Conflict of Interest Statement

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.

Acknowledgments

This study was supported by the National Natural Science Foundation of China (31372057), Beijing Postdoctoral Research Foundation (2017-22-006), the Beijing Collaborative Innovation Center for Eco-environmental Improvement with Forestry and Fruit Trees (No. CEFF-PXM2017_014207_000043), the Agricultural Science and Technology Project of Beijing Municipal Commission of Rural Affairs (No. 20150121), the Opening Project of Beijing Key Laboratory of New Technology in Agricultural Application (No. kf2017021 and No. KF2016024), the Beijing Key Project of Science and Technology Plan (D161100000716003), the Scientific Research Improvement Project of Beijing University of Agriculture (GZL2014003).

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fchem.2018.00009/full#supplementary-material

References

Arnao, M. B., Cano, A., and Acosta, M. (2001). The hydrophilic and lipophilic contribution to total antioxidant activity. Food Chem. 73, 239–244. doi: 10.1016/S0308-8146(00)00324-1

CrossRef Full Text | Google Scholar

Bachiega, T. F., Orsatti, C. L., Pagliarone, A. C., and Sforcin, J. M. (2012). The effects of propolis and its isolated compounds on cytokine production by murine macrophages. Phytother. Res. 26, 1308–1313. doi: 10.1002/ptr.3731

PubMed Abstract | CrossRef Full Text | Google Scholar

Baerenfaller, K., Massonnet, C., Walsh, S., Baginsky, S., Buhlmann, P., Hennig, L., et al. (2012). Systems-based analysis of Arabidopsis leaf growth reveals adaptation to water deficit. Mol. Syst. Biol. 8, 606. doi: 10.1038/msb.2012.39

PubMed Abstract | CrossRef Full Text | Google Scholar

Bagchi, D., Sen, C. K., Bagchi, M., and Atalay, M. (2004). Anti-angiogenic, antioxidant, and anti-carcinogenic properties of a novel anthocyanin-rich berry extract formula. Biochemistry 69, 75–80. doi: 10.1023/B:BIRY.0000016355.19999.93

PubMed Abstract | CrossRef Full Text | Google Scholar

Bertalanffy, L. V. (1968). General System Theory [Online]. Available online at: http://www.panarchy.org/vonbertalanffy/systems.1968.html

Bondonno, C. P., Yang, X., Croft, K. D., Considine, M. J., Ward, N. C., Rich, L., et al. (2012). Flavonoid-rich apples and nitrate-rich spinach augment nitric oxide status and improve endothelial function in healthy men and women: a randomized controlled trial. Free Radic. Biol. Med. 52, 95–102. doi: 10.1016/j.freeradbiomed.2011.09.028

PubMed Abstract | CrossRef Full Text | Google Scholar

Boss, P. K., Davies, C., and Robinson, S. P. (1996). Expression of anthocyanin biosynthesis pathway genes in red and white grapes. Plant Mol. Biol. 32, 565–569. doi: 10.1007/BF00019111

PubMed Abstract | CrossRef Full Text | Google Scholar

Boudet, A. M. (2007). Evolution and current status of research in phenolic compounds. Phytochemistry 68, 2722–2735. doi: 10.1016/j.phytochem.2007.06.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Candeias, L. P., Everett, S. A., and Wardman, P. (1993). Free radical intermediates in the oxidation of flavone-8-acetic acid: possible involvement in its antitumour activity. Free Radic. Biol. Med. 15, 385–394. doi: 10.1016/0891-5849(93)90038-V

PubMed Abstract | CrossRef Full Text | Google Scholar

Chauhan, A., and Vaish, R. (2012). Magnetic material selection using multiple attribute decision making approach. Mater. Des. 36, 1–5. doi: 10.1016/j.matdes.2011.11.021

CrossRef Full Text | Google Scholar

Chyu, K. Y., Babbidge, S. M., Zhao, X., Dandillaya, R., Rietveld, A. G., Yano, J., et al. (2004). Differential effects of green tea-derived catechin on developing versus established atherosclerosis in apolipoprotein E-null mice. Circulation 109, 2448–2453. doi: 10.1161/01.CIR.0000128034.70732.C2

PubMed Abstract | CrossRef Full Text | Google Scholar

Crozier, A., Jaganath, I. B., and Clifford, M. N. (2009). Dietary phenolics: chemistry, bioavailability and effects on health. Nat. Prod. Rep. 26, 1001–1043. doi: 10.1039/b802662a

PubMed Abstract | CrossRef Full Text | Google Scholar

Cushnie, T. P., and Lamb, A. J. (2005). Antimicrobial activity of flavonoids. Int. J. Antimicrob. Agents 26, 343–356. doi: 10.1016/j.ijantimicag.2005.09.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Czank, C., Cassidy, A., Zhang, Q., Morrison, D. J., Preston, T., Kroon, P. A., et al. (2013). Human metabolism and elimination of the anthocyanin, cyanidin-3-glucoside: a 13C-tracer study. Am. J. Clin. Nutr. 97, 995–1003. doi: 10.3945/ajcn.112.049247

CrossRef Full Text | Google Scholar

Espley, R. V., Hellens, R. P., Putterill, J., Stevenson, D. E., Kutty-Amma, S., and Allan, A. C. (2007). Red colouration in apple fruit is due to the activity of the MYB transcription factor, MdMYB10. Plant J. 49, 414–427. doi: 10.1111/j.1365-313X.2006.02964.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Falcone Ferreyra, M. L., Rius, S. P., and Casati, P. (2012). Flavonoids: biosynthesis, biological functions, and biotechnological applications. Front. Plant Sci. 3:222. doi: 10.3389/fpls.2012.00222

PubMed Abstract | CrossRef Full Text

Grotewold, E. (2006). The Science of Flavonoids. New York, NY: Springer.

Henry-Kirk, R. A., McGhie, T. K., Andre, C. M., Hellens, R. P., and Allan, A. C. (2012). Transcriptional analysis of apple fruit proanthocyanidin biosynthesis. J. Exp. Bot. 63, 5437–5450. doi: 10.1093/jxb/ers193

PubMed Abstract | CrossRef Full Text | Google Scholar

Hernandez, I., Alegre, L., Van Breusegem, F., and Munne-Bosch, S. (2009). How relevant are flavonoids as antioxidants in plants? Trends Plant Sci. 14, 125–132. doi: 10.1016/j.tplants.2008.12.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Hou, J. J., Wu, W. Y., Da, J., Yao, S., Long, H.-L., Yang, Z., et al. (2011). Ruggedness and robustness of conversion factors in method of simultaneous determination of multi-components with single reference standard. J. Chromatogr. A 1218, 5618–5627. doi: 10.1016/j.chroma.2011.06.058

PubMed Abstract | CrossRef Full Text | Google Scholar

Jia, H., Aw, W., Saito, K., Hanate, M., Hasebe, Y., and Kato, H. (2014). Eggshell membrane ameliorates hepatic fibrogenesis in human C3A cells and rats through changes in PPARγ-Endothelin 1 signaling. Sci. Rep. 4:7473. doi: 10.1038/srep07473

PubMed Abstract | CrossRef Full Text | Google Scholar

Jorgensen, K., Rasmussen, A. V., Morant, M., Nielsen, A. H., Bjarnholt, N., Zagrobelny, M., et al. (2005). Metabolon formation and metabolic channeling in the biosynthesis of plant natural products. Curr. Opin. Plant Biol. 8, 280–291. doi: 10.1016/j.pbi.2005.03.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Jugde, H., Nguy, D., Moller, I., Cooney, J. M., and Atkinson, R. G. (2008). Isolation and characterization of a novel glycosyltransferase that converts phloretin to phlorizin, a potent antioxidant in apple. FEBS J. 275, 3804–3814. doi: 10.1111/j.1742-4658.2008.06526.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Lewis, B. J., Herrlinger, K. A., Craig, T. A., Mehring-Franklin, C. E., Defreitas, Z., and Hinojosa-Laborde, C. (2013). Antihypertensive effect of passion fruit peel extract and its major bioactive components following acute supplementation in spontaneously hypertensive rats. J. Nutr. Biochem. 24, 1359–1366. doi: 10.1016/j.jnutbio.2012.11.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, D., Peng, C., Xie, X., Mao, Y., Li, M., Cao, Z., et al. (2014). Antidiabetic effect of flavonoids from Malus toringoides (Rehd.) Hughes leaves in diabetic mice and rats. J. Ethnopharmacol. 153, 561–567. doi: 10.1016/j.jep.2014.02.026

PubMed Abstract | CrossRef Full Text | Google Scholar

Lommen, A., Godejohann, M., Venema, D. P., Hollman, P. C., and Spraul, M. (2000). Application of directly coupled HPLC-NMR-MS to the identification and confirmation of quercetin glycosides and phloretin glycosides in apple peel. Anal. Chem. 72, 1793–1797. doi: 10.1021/ac9912303

PubMed Abstract | CrossRef Full Text | Google Scholar

Lu, Y., Wu, N., Fang, Y., Shaheen, N., and Wei, Y. (2017). An automatic on-line 2,2-diphenyl-1-picrylhydrazyl-high performance liquid chromatography method for high-throughput screening of antioxidants from natural products. J. Chromatogr. A 1521, 100–109. doi: 10.1016/j.chroma.2017.09.030

PubMed Abstract | CrossRef Full Text | Google Scholar

Magurran, A. E. (2003). Measuring Biological Diversity. Malden, MA: Wiley-Blackwell.

Google Scholar

Martens, S., and Mithofer, A. (2005). Flavones and flavone synthases. Phytochemistry 66, 2399–2407. doi: 10.1016/j.phytochem.2005.07.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Morris, E. K., Caruso, T., Buscot, F., Fischer, M., Hancock, C., Maier, T. S., et al. (2014). Choosing and using diversity indices: insights for ecological applications from the German Biodiversity Exploratories. Ecol. Evol. 4, 3514–3524. doi: 10.1002/ece3.1155

PubMed Abstract | CrossRef Full Text | Google Scholar

Nair, S., Ziaullah, Z., and Rupasinghe, H. P. V. (2014). Phloridzin fatty acid esters induce apoptosis and alters gene expression in human liver cancer cells (261.2). FASEB J. 28(1Suppl.), 261.2. doi: 10.1096/fasebj.28.1_supplement.261.2

CrossRef Full Text

Nakabayashi, R., Yonekura-Sakakibara, K., Urano, K., Suzuki, M., Yamada, Y., Nishizawa, T., et al. (2014). Enhancement of oxidative and drought tolerance in Arabidopsis by overaccumulation of antioxidant flavonoids. Plant J. 77, 367–379. doi: 10.1111/tpj.12388

PubMed Abstract | CrossRef Full Text | Google Scholar

Panda, S., and Kar, A. (2007). Annona squamosa seed extract in the regulation of hyperthyroidism and lipid-peroxidation in mice: possible involvement of quercetin. Phytomedicine 14, 799–805. doi: 10.1016/j.phymed.2006.12.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Pandey, A., Misra, P., Bhambhani, S., Bhatia, C., and Trivedi, P. K. (2014). Expression of Arabidopsis MYB transcription factor, AtMYB111, in tobacco requires light to modulate flavonol content. Sci. Rep. 4:5018. doi: 10.1038/srep05018

PubMed Abstract | CrossRef Full Text | Google Scholar

Patras, A., Brunton, N. P., Downey, G., Rawson, A., Warriner, K., and Gernigon, G. (2011). Application of principal component and hierarchical cluster analysis to classify fruits and vegetables commonly consumed in Ireland based on in vitro antioxidant activity. J. Food Compos. Anal. 24, 250–256. doi: 10.1016/j.jfca.2010.09.012

CrossRef Full Text | Google Scholar

Pietta, P. G. (2000). Flavonoids as antioxidants. J. Nat. Prod. 63, 1035–1042. doi: 10.1021/np9904509

PubMed Abstract | CrossRef Full Text | Google Scholar

Qin, X., Xing, Y. F., Zhou, Z., and Yao, Y. (2015). Dihydrochalcone compounds isolated from crabapple leaves showed anticancer effects on human cancer cell lines. Molecules 20, 21193–21203. doi: 10.3390/molecules201219754

PubMed Abstract | CrossRef Full Text | Google Scholar

Revilla, E., and Ryan, J. M. (2000). Analysis of several phenolic compounds with potential antioxidant properties in grape extracts and wines by high-performance liquid chromatography-photodiode array detection without sample preparation. J. Chromatogr. A 881, 461–469. doi: 10.1016/S0021-9673(00)00269-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Roig, R., Cascon, E., Arola, L., Blade, C., and Salvado, M. J. (2002). Procyanidins protect Fao cells against hydrogen peroxide-induced oxidative stress. Biochim. Biophys. Acta 1572, 25–30. doi: 10.1016/S0304-4165(02)00273-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Roleira, F. M. F., Tavares-da-Silva, E. J., Varela, C. L., Costa, S. C., Silva, T., Garrido, J., et al. (2015). Plant derived and dietary phenolic antioxidants: anticancer properties. Food Chem. 183, 235–258. doi: 10.1016/j.foodchem.2015.03.039

PubMed Abstract | CrossRef Full Text | Google Scholar

Sáez-Ayala, M., Sánchez-Del-Campo, L., Montenegro, M. F., Chazarra, S., Tárraga, A., and Cabezas-Herrera, J. (2015). Comparison of a pair of synthetic tea-catechin-derived epimers: synthesis, antifolate activity, and tyrosinase-mediated activation in melanoma. ChemMedChem 6, 440–449. doi: 10.1002/cmdc.201000482

PubMed Abstract | CrossRef Full Text

Salvatierra, A., Pimentel, P., Moya-Leon, M. A., and Herrera, R. (2013). Increased accumulation of anthocyanins in Fragaria chiloensis fruits by transient suppression of FcMYB1 gene. Phytochemistry 90, 25–36. doi: 10.1016/j.phytochem.2013.02.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Sanchez-Rabaneda, F., Jauregui, O., Lamuela-Raventos, R. M., Viladomat, F., Bastida, J., and Codina, C. (2004). Qualitative analysis of phenolic compounds in apple pomace using liquid chromatography coupled to mass spectrometry in tandem mode. Rapid Commun. Mass Spectrom. 18, 553–563. doi: 10.1002/rcm.1370

PubMed Abstract | CrossRef Full Text | Google Scholar

Scheller, S., Dworniczak, S., Waldemarklimmek, K., Rajca, M., Tomczyk, A., and Shani, J. (1999). Synergism between ethanolic extract of propolis (EEP) and anti-tuberculosis drugs on growth of mycobacteria. Z. Naturforschung C 54, 549–553. doi: 10.1515/znc-1999-7-814

PubMed Abstract | CrossRef Full Text | Google Scholar

Shi, S.-Y., Zhang, Y.-P., Jiang, X.-Y., Chen, X.-Q., Huang, K.-L., and Zhou, H.-H. (2009). Coupling HPLC to on-line, post-column (bio)chemical assays for high-resolution screening of bioactive compounds from complex mixtures. TrAC Trends Analyt. Chem. 28, 865–877. doi: 10.1016/j.trac.2009.03.009

CrossRef Full Text | Google Scholar

Shirley, B. W. (1996). Flavonoid biosynthesis: ‘new’ functions for an ‘old’ pathway. Trends Plant Sci. 1, 377–382.

Google Scholar

Simao, M. C. M., Flory, S. L., and Rudgers, J. A. (2010). Experimental plant invasion reduces arthropod abundance and richness across multiple trophic levels. Oikos 119, 1553–1562. doi: 10.1111/j.1600-0706.2010.18382.x

CrossRef Full Text | Google Scholar

Srinivasan, A., Thangavel, C., Liu, Y., Shoyele, S., Den, R. B., Selvakumar, P., et al. (2016). Quercetin regulates β-catenin signaling and reduces the migration of triple negative breast cancer. Mol. Carcinog. 55:743. doi: 10.1002/mc.22318

PubMed Abstract | CrossRef Full Text | Google Scholar

Stracke, R., Ishihara, H., Huep, G., Barsch, A., Mehrtens, F., Niehaus, K., et al. (2007). Differential regulation of closely related R2R3-MYB transcription factors controls flavonol accumulation in different parts of the Arabidopsis thaliana seedling. Plant J. 50, 660–677. doi: 10.1111/j.1365-313X.2007.03078.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Sugiura, C., Nishimatsu, S., Moriyama, T., Ozasa, S., Kawada, T., and Sayama, K. (2012). Catechins and caffeine inhibit fat accumulation in mice through the improvement of hepatic lipid metabolism. J. Obes. 2012:520510. doi: 10.1155/2012/520510

PubMed Abstract | CrossRef Full Text | Google Scholar

Symons, M. C. R., and Gutteridge, J. M. C. (1998). Free Radicals and Iron: Chemistry, Biology and Medicine. Oxford: Oxford University Press.

Google Scholar

Tian, J., Peng, Z., Zhang, J., Song, T., Wan, H., Zhang, M., et al. (2015). McMYB10 regulates coloration via activating McF3'H and later structural genes in ever-red leaf crabapple. Plant Biotechnol. J. 13, 948–961. doi: 10.1111/pbi.12331

PubMed Abstract | CrossRef Full Text | Google Scholar

Tian, J., Shen, H., Zhang, J., Song, T., and Yao, Y. (2011). Characteristics of chalcone synthase promoters from different leaf-color malus crabapple cultivars. Sci. Hortic. 129, 449–458. doi: 10.1016/j.scienta.2011.04.015

CrossRef Full Text | Google Scholar

Tsao, R. Y. R., Young, J. C., and Zhu, H. (2003). Polyphenolic profiles in eight apple cultivars using high-performance liquid chromatography (HPLC). J. Agric. Food Chem. 51, 6347–6353. doi: 10.1021/jf0346298

PubMed Abstract | CrossRef Full Text | Google Scholar

Tsuji, P. A., Stephenson, K. K., Wade, K. L., Liu, H., and Fahey, J. W. (2013). Structure-activity analysis of flavonoids: direct and indirect antioxidant, and antiinflammatory potencies and toxicities. Nutr. Cancer 65, 1014–1025. doi: 10.1080/01635581.2013.809127

PubMed Abstract | CrossRef Full Text | Google Scholar

Velasco, R., Zharkikh, A., Affourtit, J., Dhingra, A., Cestaro, A., Kalyanaraman, A., et al. (2010). The genome of the domesticated apple (Malus × domestica Borkh.). Nat. Genet. 42, 833–839. doi: 10.1038/ng.654

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, S., Meckling, K. A., Marcone, M. F., Kakuda, Y., and Tsao, R. (2011). Synergistic, additive, and antagonistic effects of food mixtures on total antioxidant capacities. J. Agric. Food Chem. 59, 960–968. doi: 10.1021/jf1040977

PubMed Abstract | CrossRef Full Text | Google Scholar

Winkel-Shirley, B. (2001). It takes a garden. how work on diverse plant species has contributed to an understanding of flavonoid metabolism. Plant Physiol. 127, 1399–1404. doi: 10.1104/pp.010675

PubMed Abstract | CrossRef Full Text | Google Scholar

Winterhalter, P., Güldner, A., Jakob, U., and Schreier, P. (1994). 3-Hydroxy-5, 6-epoxy-β-ionyl-β-D-glucopyranoside from Malus domestica leaves. Nat. Prod. Lett. 4, 57–60. doi: 10.1080/10575639408043892

CrossRef Full Text | Google Scholar

Xiao, J., Muzashvili, T. S., and Georgiev, M. I. (2014). Advances in the biotechnological glycosylation of valuable flavonoids. Biotechnol. Adv. 32, 1145–1156. doi: 10.1016/j.biotechadv.2014.04.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Yin, R., Messner, B., Faus-Kessler, T., Hoffmann, T., Schwab, W., Hajirezaei, M. R., et al. (2012). Feedback inhibition of the general phenylpropanoid and flavonol biosynthetic pathways upon a compromised flavonol-3-O-glycosylation. J. Exp. Bot. 63, 2465–2478. doi: 10.1093/jxb/err416

PubMed Abstract | CrossRef Full Text | Google Scholar

Yi, H., Chen, Y., Liu, J., Zhang, J., Guo, W., Xiao, W., et al. (2016). Extraction and separation of active ingredients in Schisandra chinensis (Turcz.) Baill and the study of their antifungal effects. PLoS ONE 11:e0154731. doi: 10.1371/journal.pone.0154731

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, N. (2007). Phylogenetic relationship between ornamental and wild species of Malus in china. Acta Hortic. Sin. 34, 1227–1234. doi: 10.16420/j.issn.0513-353x.2007.05.031

CrossRef Full Text | Google Scholar

Zhou, J.-L., Li, P., Li, H.-J., Jiang, Y., Ren, M.-T., and Liu, Y. (2008). Development and validation of a liquid chromatography/electrospray ionization time-of-flight mass spectrometry method for relative and absolute quantification of steroidal alkaloids in Fritillaria species. J. Chromatogr. A 1177, 126–137. doi: 10.1016/j.chroma.2007.11.030

PubMed Abstract | CrossRef Full Text | Google Scholar