Illustration of the variation in the content of flavanone rutinosides in various citrus germplasms from genetic and enzymatic perspectives

Abstract In citrus, 1,6-rhamnosytransferase (1,6RhaT) and 1,2-rhamnosytransferase (1,2RhaT) catalyze flavanone-7-O-glucosides to form nonbitter flavanone rutinosides (FRs) and bitter flavanone neohesperidosides (FNs), respectively. As revealed in this study of fruit peels from 36 citrus accessions, FRs varied from undetectable levels in pummelo and kumquat to being the dominant flavonoids in sweet orange and loose-skin mandarins. Furthermore, a previously annotated full-length 1,6RhaT-like gene was identified as another 1,6RhaT-encoding gene by in vitro experiments. In total, 28 alleles of full-length 1,6RhaTs were isolated and classified into A, B and C types with only type A alleles encoding a functional protein. Coincidently, only the accessions that contained FRs harbored type A alleles, as was further verified in two F1 hybrid populations. Moreover, the inferior substrate conversion efficiency of 1,6RhaTs in comparison with that of 1,2RhaT in vitro might partly explain the lower proportions of FRs to total flavanone disaccharides in citrus hybrids harboring both functional rhamnosyltransferases. Our findings provide a better understanding of FR content variations among citrus and are meaningful for a mechanistic illustration of citrus flavonoid metabolism and fruit quality improvement practices.

A considerable number of studies have investigated the specific accumulation of FRs and FNs in citrus germplasms and the changes in the contents of these compounds among different developmental stages or tissues [15][16][17][18]. High levels of FRs are present in orange, mandarin and lemon; low levels of FRs are present in grapefruit and bitter orange, which contain high levels of FNs; and FRs were not detected in pummelo [19]. The available evidence shows that the accumulation of FRs varies among citrus germplasms, but the mechanisms are not well understood. Fortunately, based on our previous elucidation of the molecular mechanism underlining the variations in the FN profiles among different citrus germplasms [20], revealing the mechanisms underlying the variations in the FR profiles in different citrus germplasms is feasible.
In this study, the FRs in various citrus accessions covering the main cultivated species, wild species and their relatives were investigated. Another 1,6RhaT gene was identified, which implies that there are multiple 1,6RhaT genes in citrus. The 1,6RhaT genes and their allelic variations were then characterized in vitro. Furthermore, the production of FNs and FRs in two hybrid cross populations harboring both functional 1,2RhaT and 1,6RhaT genes was further investigated. This study will contribute to a better understanding of the accumulation of f lavanone disaccharides with different bioactive characteristics in various citrus germplasms and is thus meaningful for improving fruit quality.

Plant materials and compounds
Thirty-four accessions were obtained from the National Citrus Breeding Centre (NCBC) at Huazhong Agricultural University, Wuhan, China, and two other accessions were obtained from Haohuahong Town, Huishui County, Guizhou Province, China (Table S1). For each accession, young leaves and 15-18 fruits were randomly collected from at least three healthy trees at their respective commercial maturation stages at approximately 230 d post anthesis (DPA).
All fruits and leaves were randomly divided into three biological replications. The fruit peel (f lavedo and albedo) tissues were immediately separated, frozen in liquid nitrogen, ground into a fine powder and stored at −80 • C. All operations were performed under dim light conditions, and long-term exposure to open air was avoided to prevent the degradation of f lavonoids.
Transgenic tobacco bright yellow (BY2) calli harboring Cm1,2RhaT were obtained as described by Chen

Extraction of flavonoids and determination of FRs and FNs
The extraction and separation of the flavanone disaccharides were conducted as described previously [17]. Flavonoids were separated using a gradient elution program with a flow rate of 1 ml/min, a column temperature of 35 • C and a Waters 1525 binary system equipped with a 2996 photodiode array detector, a 717 Plus autosampler (Waters Corp., Milford, MA, USA) and a C18 column (250 mm × 4.6 mm, 5 μm, Thermo Scientific, Waltham, MA, USA). Subsequently, 0.15% formic acid diluted in water was prepared as mobile phase A, and 0.15% formic acid-acetonitrile was prepared as mobile phase B. The data were analyzed using Empower Chromatography Manager software (Waters Co., Milford, MA, USA).

RNA isolations, cDNA synthesis and real-time PCR
Total RNA was isolated using an EASYspin Plus Plant RNA Kit (Aidlab Biotechnologies Co., Ltd.). HiScript II Reverse Transcriptase (Aidlab Biotechnologies Co., Ltd.) was used for the amplification of cDNA. qRT-PCR was performed using a Prism 7900HT (Applied Biosystems) with the SYBR Green system (Hieff qPCR SYBR Green Master Mix, Yeasen Biotech Co., Ltd.). All gene expression analyses were performed with three independent biological replicates. The primers used in the analysis are listed in Supporting Table S2. The actin gene was used as the reference control.

BLAST analysis, phylogenetic analysis and gene cloning
BLAST analysis was performed against the Phytozome 12 database (https://phytozome.jgi.doe.gov/pz/portal. html#!search?show=BLAST) and Citrus sinensis Genome Annotation Project database (http://citrus.hzau.edu.cn/ cgi-bin/orange/blast) using Cs1,6RhaT (accession number: DQ119035) as the reference sequence. ClustalW (http://www.genome.jp/tools-bin/clustalw) was utilized for multiple alignments of the nucleotide sequences of the alleles in 36 accessions. Phylogenetic analysis was generated using MEGA 7.0 according to the maximum likelihood method. There were 1000 bootstrap replications to test the reliability.
Young leaves of 36 accessions and two F1 hybrid populations were collected for genomic DNA extraction according to the protocol described by Cheng et al. [21] To amplify the complete coding sequence of 1,6RhaT genes from both gDNA and cDNA templates, specific PCR primers were designed based on Cs1,6RhaT (Supporting Table S2). The primers were synthesized by Sangon Biotech Co., Ltd. (Shanghai, China). To ensure accuracy, at least eight positive clones of each accession were sequenced by TSINGKE Biological Technology Co., Ltd. (Beijing, China).

Expression and biotransformation assay in tobacco BY2 cells
Seven 1,6RhaT alleles that were found at a higher frequency (A1, A5, A10, B2 and C4) during isolation or appeared in some unique accessions (A8 and C6) were cloned into the pH7WG2D vector under the control of a CAMV-35S promoter harboring green f luorescent protein (GFP) and transferred into Agrobacterium GV3101. The resulting construct was used to stably transform BY 2 callus suspensions according to a published protocol [20,22]. The transgenic calluses were confirmed by fluorescence detection through epif luorescence stereomicroscopy at 488 nm, PCR amplification and sequencing. The positively transformed calli were subjected to subsequent substrate feeding analysis.
To carry out the biotransformation assays, a substrate feeding experiment using BY2 cells was performed according to the protocol described by Frydman et al. [13]. Brief ly, transgenic and wild-type BY 2 calli were suspended in Erlenmeyer f lasks for one week at 28 • C with shaking at 140 rpm. Subsequently, 1 g of callus was diluted in 20 mL of MS liquid media and incubated under the same conditions for 4 d. Flavonoid substrates (mg/ml) were added to the f lask and then incubated for another 2 d. The resulting BY 2 cells were harvested and used for extraction and determination of f lavonoids.

Protein expression and purification
The expression and purification of Cme1,6RhaT-like was performed as described by Peng et al. [23]. In brief, Cme1,6RhaT-like was cloned into the pMal-c2x expression vector with a maltose-binding protein tag. Recombinant proteins were expressed in BL21 (DE3) cells and induced with IPTG at 16 • C. The harvested cells were lysed in phosphate-buffered saline pH 7.4 containing 400 mM NaCl by sonication at 4 • C. The protein was transferred into a column containing dextrin beads 6FF and purified by washing buffer (20 mM Tris-HCl, 20 mM NaCl, 1 mM EDTA and 1 mM DTT, pH 7.4) and elution buffer (20 mM Tris-HCl, 1 mM EDTA, 10 mM maltose, and 1 mM DTT, pH 7.4).
Cp1,2RhaT and Cp1,6RhaT were expressed and purified following a similar protocol.

Assay enzyme activity and kinetics
The enzyme assay was performed according to Peng et al. [23]. The glycosylation reaction was performed in a total volume of 100 μl containing 5 μg of purified glycosyltransferases, 5 μg of f lavonoid substrates, 5 μg of UDPrhamnose, and 0.25 mM MgCl 2 in phosphate buffered saline (50 mM, pH 7.4). The reaction was then incubated for 30 min at 35 • C. The reaction was stopped by the addition of 300 μl of ice-cold methanol, and the resulting products were subjected to LC-MS analysis.

FR accumulation is species-dependent in various citrus germplasms
In our previous study, flavonoids were found to accumulate at higher levels in fruit peel (flavedo and albedo) than in the pulp of most citrus germplasms. Consequently, substantial variations in the contents of FRs and FNs were observed among the fruit peels of 36 accessions (Supplement Table S1). These accessions cover the main cultivated species, wild species and their relatives, including trifoliata orange (Poncirus trifoliata), kumquat (Citrus japonica), loose-skin mandarin (Citrus reticulata), sweet orange (Citrus sinensis), pummelo (Citrus maxima), grapefruit (Citrus paradisi), lemon (Citrus limon), citron (Citrus medica) and ichang papeda (Citrus ichangensis).
These citrus accessions could be divided into three groups according to their contents of FRs (Fig. 1). The accessions in Group 1, such as sweet orange, looseskin mandarin and citrons, had high levels of FRs but undetectable levels of FNs. Group 2 included accessions with both FRs and FNs, such as ichang papeda, trifoliate orange and grapefruit. However, the proportions of FR among the total flavanone disaccharides were lower than those of FNs. Group 3 was comprised of pummelo and kumquat and exhibited undetectable levels of FRs; moreover, pummelo was enriched in FNs, whereas FNs were undetectable in kumquat. These results showed that the citrus FR profiles are species-dependent. The associations among the contents and proportions of the FRs and the genetic background in various citrus accessions was then investigated.

Identification of a new 1,6RhaT gene in citrus genomes
A 1,6RhaT-like gene with high nucleotide sequence homology to Cs1,6RhaT (GenBank accession No. AY048882) was identified in the genomes of various citrus species (Fig. 2). The 1,6RhaT-like gene was annotated in the genomes of various citrus species with the exception of sweet orange and loose-skin mandarin and it was located on different loci from those of 1,6RhaT. The genomic sequences revealed that 1,6RhaT-like shared high nucleotide sequence identities (approximately 96%) with Cs1,6RhaT, and various base deletions and insertions were found in Atalantia buxifolia, C. maxima  and C. ichangensis, whereas an intact 1,6RhaT-like was only found in Citrus medica with an ORF of 1428 bp. Cme1,6RhaT-like from Eureka lemon was isolated, sequenced and subjected to an analysis of its catalytic activity. Cme1,6RhaT-like was expressed in Escherichia coli and characterized. After the application of naringenin-7-O-glucoside and UDP-rhamnose as substrates, a new peak was detected and identified as narirutin by comparing the retention time and the MS2 spectra of the corresponding product of Cs1,6RhaT (Fig. 3a)  its catalytic activity, we used other types of f lavonoids, including f lavone (kaempferol-3-O-glucoside) and flavonol (diosmetin-7-O-glucoside), as substrates ( Fig. 3b-c).
The results showed that Cme1,6RhaT-like also exhibited rhamnosylation activity toward these two types of flavonoids and converted these substrates into their corresponding f lavanone disaccharides. Thus, our results showed that Cme1,6RhaT-like is another 1,6RhaT gene, and multiple 1,6RhaT genes are found in citrus. Consequently, we defined 1,6RhaT and  1,6RhaT-like as 1,6RhaTs in the rest of the manuscript.

Allelic sequence variations of citrus 1,6RhaTs
To examine the correlation between the 1,6RhaT genotype and the rutinoside accumulation phenotype, 1,6RhaT was amplified from the DNA of 36 citrus accessions. A total of 28 different alleles were isolated with high similarity, and insertions, deletions and single nucleotide polymorphisms (SNPs) were also found. These alleles can be divided into three types, A, B and C, based on their sequences. Type A contains 18 alleles defined as A1 to A18, type B consists of B1 to B4, and C contains C1 to C6. Each accession investigated possessed two different  ranged from 1420-1433 bp in length, shared 95.3%-96.0% nucleotide identity with the reference gene, and they generated premature proteins, which potentially do not exhibit 1,6RhaT functions (Supporting Table S4). Nonsense and frameshift mutations of the type B alleles were caused by the insertion of five base pairs, whereas base deletions at positions 185 or 984 led to frameshift mutations among the type C alleles. Additionally, more than 50 SNPs were observed for each member of the type B and C alleles. These alleles could be amplified from the genomes of pummelo, kumquat and hybrids such as grapefruit and citrons. Among the type A alleles, A10 appeared most frequently and could be isolated from 18 accessions. For type B and type C alleles, the B2 allele showed the highest frequency during isolation.

Citrus 1,6RhaT type a alleles encode proteins with 1,6RhaT function
As indicated above, premature stop codons in the type B and C alleles were hypothesized to produce truncated proteins without 1,6RhaT enzymatic function, which would lead to the unsuccessful accumulation of FRs. To validate this hypothesis, the catalytic activities of representative 1,6RhaT alleles were analyzed in suspensions of transgenic BY2 tobacco calli. Wild-type BY2 and transgenic cell lines for the A10 and B2 alleles were given hesperetin or naringenin substrate for 48 h, and the resulting biotransformation products were extracted and identified by HPLC (Fig. 4b-c). New peaks were detected only in the transgenic cell lines for the A10 allele but not in the B2 allele and BY2 control. The products in the A10 allele cell lines were identified as hesperedin (hesperetin-7-O-rutinoside) and narirutin (naringenin-7-O-rutinoside) by comparison with standards. Thus, the substrate could be transformed into FRs in cells transformed with the A10 allele but this failed in cells with B2 and in the wild type, which indicated that only type A alleles have 1,6RhaT function.

FR accumulation is correlated with the presence of type a alleles in citrus
The results from the 1,6RhaT allele analysis with the FR contents among the 36 accessions suggested a significant  ZJ72  ZJ80  ZJ58  ZJ54  ZJ56  ZJ7  ZJ97  ZJ10  ZJ41  ZJ71  ZJ1  ZJ50  ZJ26  ZJ16  ZJ20  ZJ82  ZJ84  ZJ61  ZJ63  ZJ33  ZJ49  ZJ29  ZJ25  ZJ22  ZJ37  ZJ77  ZJ40  ZJ85  ZJ21  ZJ60  ZJ31  ZJ96  ZJ106  ZJ87  ZJ19  ZJ27  ZJ62  ZJ39  ZJ89  ZJ14  ZJ79  ZJ102  ZJ11  ZJ9  ZJ90  ZJ81  ZJ104  ZJ47  ZJ5  ZJ34  ZJ55  ZJ95  ZJ109  ZJ30  ZJ66  ZJ99  ZJ59  ZJ98  ZJ35  ZJ57  ZJ36  ZJ70  ZJ108  ZJ2  ZJ67  association between allelic diversity and FR phenotypes (Supporting Table S5). FRs were detectable in accessions containing at least one type A allele, such as accessions belonging to sweet orange, loose-skin mandarin, trifoliate orange, ichang papeda and citrons. The undetectable FRs in pummelo and kumquat possessed type B and/or C but not type A alleles. Moreover, the contents of FRs in the accessions with two type A alleles (21.81 mg/g DW) were significantly higher than those with one type A allele (5.39 mg/g DW, P < 0.05). We hypothesized that the accumulation of FR was correlated with the presence of type A alleles in citrus.
To confirm the hypothesis of the association between allelic diversity and FR phenotypes, populations of "HB pummelo", "Fairchild tangelo" and their F 1 hybrids were analyzed (Fig. 5a). The female parent HB pummelo contained B2 and C4 alleles with undetectable FRs in fruit peel, and the male parent "Fairchild tangelo" contained A5 and A10 with high levels of FRs in fruits (54.37 mg/g DW). Each F1 progeny obtained from the crossing carries a type B or C allele copy, and a type A allele copy, and thus it should accumulate FR (Supporting Fig. S1). HPLC analysis showed that all 37 F 1 hybrids accumulated FR in the peel, which confirms our hypothesis.

The proportion of FRs among flavanone disaccharides is affected by the genotypes of both 1,6RhaT and 1,2RhaT in the hybrids
The crossing of "HB pummelo" and "Fairchild tangelo" resulted in extremely significantly lower levels of FRs in the F1 hybrid population, which ranged from 1.05 to 5.71 mg/g DW with an average content of 2.54 mg/g DW in peel. The same result was also found in natural hybrids of trifoliate orange and grapefruit among the Group 2 accessions (Fig. 4). Therefore, a question raised in this study is what is the cause of the relatively low contents of FRs in these hybrids.
The abovementioned progenies contain functional 1,2RhaT to produce FNs, which are a subgroup of flavanone disaccharides that compete with FRs converted from the same substrates of flavanone-7-Oglucosides (Supporting Fig. S2). Further determination of the contents of FNs (e.g. naringin and neohesperidin) showed that FRs were found at lower proportions, ranging from 5.53% to 30.20%, whereas FNs were found at proportions of 69.80%-94.47% (Fig. 5a), which implied that the synthesis of FNs might have negative influences on the accumulation of FRs.
To determine the above influences, F1 progenies resulting from another crossing of "red tangerine" and "trifoliata orange" were employed. Both the female and male parents in the crossing combination contained two type A 1,6RhaT alleles (red tangerine: A4 and A10, trifoliate orange: A8 and A9) (Supplementary Fig. S3). However, 1,2RhaT exists in the male parent "trifoliate orange" and 39 progenies but does not appear in the female parent "red tangerine" and 27 other progenies. High contents of FRs and undetectable levels of FNs were detected in the peel of "red tangerine", and both FRs and FNs were detected in the peel of "trifoliate orange". HPLC analysis showed that all 66 F 1 hybrids accumulated FR in the peel, and both FRs and FNs were detectable in the abovementioned 39 progenies harboring 1,2RhaT (Fig. 5b) Table S6). These results showed that due to the involvement of FN biosynthesis, low FR contents were found in 39 of 66 F1 progenies from crossings of "red tangerine" and "trifoliata orange", whereas high contents were detected in 27 other F1 progenies without FN biosynthesis. Therefore, it seems that the lower proportion of FRs among f lavanone disaccharides was caused by the synthesis of FNs, which requires further in-depth research.

Enzymatic analyses of Cp1,6-RhaT (A10) and Cp1,2-RhaT
Furthermore, to clarify the inf luences of both 1,6RhaT and 1,2RhaT on FR production by considering substrate preference or competition, enzyme kinetic assays were performed in vitro. First, 1,6RhaT and 1,2RhaT were cloned from C. paradisi, which contains both genes. After the optimal temperature and pH for the reactions catalyzed by 1,2RhaT and 1,6RhaT were identified (Supporting Fig. S4), naringenin-7-O-glucoside (N7G) and hesperetin-7-O-glucoside (H7G) were used as substrates to determine the kinetic parameters because most citrus flavonoids are produced from these two substrates. As shown in Table 2, the Km values of 1,6RhaT for N7G and H7G were 14.7 and 12.9 μM, respectively, whereas those of 1,2RhaT were 11.7 and 25.5 μM, respectively, which indicated that both enzymes exhibited strong binding affinity for these two compounds. The conversion efficiency (Kcat) of 1,2RhaT for N7G and H7G was approximately 1.5-fold higher and approximately 1.6-fold higher than that of 1,6RhaT, respectively. Both 1,2RhaT and 1,6RhaT consumed N7G more efficiently than H7G. Analysis of the ratio of the catalytic efficiency to specificity (Kcat/Km) for N7G showed that the ratio of 1,2RhaT was approximately 1.9-fold higher than that of 1,6RhaT. In contrast, the Kcat/Km ratio of 1,2RhaT for H7G was lower due to the weaker affinity between 1,2RhaT and H7G.
These results showed that the efficiency of 1,2RhaT to form FNs was higher than that of 1,6RhaT to form FRs, which might partly explain the lower FR contents found in hybrids containing both 1,2RhaT and 1,6RhaT.

Discussion
In our study, 1,6RhaT-like genes were found in the genomes of relatively primitive species such as A. buxifolia and C. medica but were missing in model species such as loose-skin mandarin and sweet orange. Many reports have shown that plant flavanone disaccharides are related to resistance to environmental stress [24,25]. For example, hesperidin accumulates in young leaves and fruits and plays an important role in defense against animal and microbial intrusion [26]. Naringin, hesperidin, and quercetin-3-O-rutinoside together with other active compounds can prevent Papilio machaon from laying eggs on young citrus leaves and can reduce the occurrence of citrus diseases and pests [27]. Identification of another 1,6RhaT (Cme1,6RhaT-like) shows that multiple 1,6RhaT genes work together to guarantee the synthesis of FRs in resisting the harsh environment of ancient times. As citrus began to be cultivated, 1,6RhaTlike might have gradually became redundant due to the improvement of cultivation conditions and is now absent in modern germplasms.
A phylogenetic tree was constructed based on the allelic variations among the 36 accessions, which not only successfully separated type A alleles from the other two types of alleles but also reflected complex genetic backgrounds among these accessions (Fig. 6). The results show that the type B and C alleles of hybrids such as grapefruits and "citrumelo" (ZY) were derived from pummelo, as these type B and C alleles clustered together. Type A alleles of grapefruits and citrange (ZC) were clustered with A10 of sweet orange, indicating that type A alleles of these accessions were derived from sweet orange. These results are consistent with the kinship of   Supplementary Table S1. these species and accessions (grapefruit: sweet orange × pummelo; citrumelo: trifoliata orange × grapefruit; citrange: sweet orange × trifoliata orange). Interestingly, sweet orange results from a hybrid of pummelo and mandarin followed by a backcross with a mandarin [28] and therefore should have inherited the type B or C allele from the pummelo parent. Notably, we found that modern sweet orange does not contain type B or C alleles. We speculated that type B and C alleles existed in the initial hybrids of pummelo and mandarin; however, the hybrid not harboring type B and type C alleles that appeared during the backcross was selected for cultivation, which consequently led to the disappearance of type B and type C alleles in modern sweet orange.
The transcript levels of key genes are usually used to elucidate the accumulation of metabolites in various plants. However, the expression levels of 1,6RhaTs showed marked differences among the 36 citrus accessions but were not well associated with the FR contents (Supporting Fig. S5). For example, the expression of 1,6RhaT was detected in species such as pummelo and kumquat that do not accumulate FRs. This may be due to the simultaneous amplification of transcripts from nonfunctional 1,6RhaTs (type B and C alleles) in these citrus species. High expression levels of 1,6RhaT were detected in citrons that accumulate low levels of FRs, and previous research has shown that low enzyme levels of 1,6RhaT result in a low content of FRs in citron [19]; moreover, flavonoid-7-O-glucosides were not detected in citrus, indicating that the catalytic efficiency of 1,2RhaT and 1,6RhaT is very strong and that the supply of substrate is not saturated. Therefore, an insufficient substrate supply of 1,6RhaT may also be a cause of the low content of FR in citron.
The synthesis of citrus FNs and FRs involves two distinct but closely related competing subpathways during the rhamnosylation step that are catalyzed by 1,6RhaT and 1,2RhaT, respectively. Unlike 1,6RhaTs, 1,2RhaT has only one locus and it is directly related to the accumulation of FNs in citrus fruits, as demonstrated through an association analysis of 1,2RhaT with the FN accumulation phenotypes of 50 citrus accessions [20]. It is well known that FRs and FNs produce citrus fruits with nonbitter and bitter tastes, respectively, and these compounds have thus become important targets in citrus flavor-related breeding.
Our study found that the normal 1,6RhaT (type A alleles) genotype is correlated with high levels of FRs in all investigated accessions, whereas accessions with early terminations of 1,6RhaT resulting from frameshift mutations (type B and C alleles) had undetectable levels of FRs. Frameshift mutations are widespread in kiwifruit [29], Brassicaceae [30], pear [31] and apricot [32] crops and result in nonfunctional truncated proteins. However, exceptions were also found in grapefruit, trifoliata orange and its citrumelo hybrid harboring the type A allele and a lower FR content. These exceptions might be due to the inf luence of other genes related to the substrate supply, such as 7-O-glucosyltransferase, which catalyzes the formation of substrates of f lavanone 7-O-glucosides, or to genes sharing the same substrates, such as 1,2RhaT, which would account for the accumulation of FRs. The accumulation characteristics of f lavanone glycosides in the two artificial hybrid populations also indicated that other factors might affect the accumulation of FRs.
Previous studies have shown that low enzyme levels of 1,6RhaT and high enzyme levels of 1,6RhaT were detected by corresponding antibodies in grapefruit and bitter orange, and different enzyme levels in germplasms containing both 1,6RhaT and 1,2RhaT may be one of the causes of germplasms with type A alleles but with low FR content [13,19]. Moreover, we think that the catalytic efficiency of 1,2RhaT and 1,6RhaT is also a factor that affects the f lavonoid content. Our results showed that the efficiency of 1,2RhaT to form FNs was higher than that of 1,6RhaT. This result is consistent with the flavonoid content in grapefruit, which might explain the lower FR contents from another perspective.
It is worth noting that kumquat, which contains functional 1,2RhaT and nonfunctional type C 1,6RhaT alleles, harbors another f lavanone accumulation mode but no detectable FNs. This finding might be due to the loss or nonfunction of upstream genes that results in a lack of flavanone-7-O-glucoside supply [20].
Collectively, our results obtained in this study could partially explain the FR metabolic characteristics of various citrus species and can provide a better understanding of the biosynthesis and regulation of FRs, which is meaningful in theory and practice for improving citrus flavor and the functional components of f lavonoid production. Using two distinct but closely related competing subpathways to reduce the bitter taste of citrus, we hope to obtain a certain degree of freedom to modulate the biosynthesis of f lavanone disaccharides between FRs and FNs in the future.