Identification of QTLs controlling aroma volatiles using a ‘Fortune’ x ‘Murcott’ (Citrus reticulata) population

Background Flavor is an important attribute of mandarin (Citrus reticulata Blanco), but flavor improvement via conventional breeding is very challenging largely due to the complexity of the flavor components and traits. Many aroma associated volatiles of citrus fruit have been identified, which are directly related to flavor, but knowledge of genetic linkages and relevant genes for these volatiles, along with applicable markers potentially for expeditious and economical marker-assisted selection (MAS), is very limited. The objective of this project was to identify single nucleotide polymorphism (SNP) markers associated with these volatile traits. Result Aroma volatiles were investigated in two mandarin parents (‘Fortune’ and ‘Murcott’) and their 116 F1 progeny using gas chromatography mass spectrometry in 2012 and 2013. A total of 148 volatiles were detected, including one acid, 12 alcohols, 20 aldehydes, 14 esters, one furan, three aromatic hydrocarbons, 16 ketones, one phenol, 27 sesquiterpenes, 15 monoterpenes, and 38 unknowns. A total of 206 quantitative trait loci (QTLs) were identified for 94 volatile compounds using genotyping data generated from a 1536-SNP Illumina GoldenGate assay. In detail, 25 of the QTLs were consistent over more than two harvest times. Forty-one QTLs were identified for 17 aroma active compounds that included 18 sesquiterpenes and were mapped onto four genomic regions. Fifty QTLs were for 14 monoterpenes and mapped onto five genomic regions. Candidate genes for some QTLs were also identified. A QTL interval for monoterpenes and sesquiterpenes on linkage group 2 contained four genes: geranyl diphosphate synthase 1, terpene synthase 3, terpene synthase 4, and terpene synthase 14. Conclusions Some fruit aroma QTLs were identified and the candidate genes in the terpenoid biosynthetic pathway were found within the QTL intervals. These QTLs could lead to an efficient and feasible MAS approach to mandarin flavor improvement. Electronic supplementary material The online version of this article (doi:10.1186/s12864-017-4043-5) contains supplementary material, which is available to authorized users.


Background
Citrus is one of the most economically important fruit commodities in the world, including in the US. Citrus includes diverse fruit types preferably produced in adapted regions and climates. For example, sweet oranges (Citrus sinensis L. Osb.), mandarins (C. reticulata Blanco, commonly referred to as tangerines in the US), and grapefruit (C. paradisi Macf.) are three main types of citrus fruit produced in Florida, US. Next to orange, mandarin production worldwide is the second largest of all citrus fruit. Whereas sweet orange is primarily for juice processing, diverse mandarins are primarily for fresh consumption. Citrus fruit and juice are a wellknown source of flavonoids, vitamin C, carotenoids, and volatile compounds which yield aromas and tastes specific to citrus fruits [1].
Citrus aroma is the result of a mixture of acids, alcohols, aldehydes, esters, ketones, monoterpenes, sesquiterpenes, and other volatiles [1][2][3]. The volatiles of sweet orange and grapefruit have been studied thoroughly because of their commercial importance in juice, and more than 300 volatiles have been reported in fresh orange juice [2]. For mandarins, several studies have been completed in peel oil and essence [4][5][6][7], and fresh fruits [8][9][10][11][12][13][14][15][16]. The aroma of mandarin juice is derived from a complex collection of several volatile compounds, rather than from any specific character-impact compounds [17]. In a review paper, a total of 37 consensus aroma volatiles were found in common across eight studies analyzing mandarins juice by gas chromatography mass spectrometry (GCMS), including seven alcohols, nine aldehydes, one ketone, 15 terpenes/hydrocarbons and five esters [3]. In addition, Miyazaki et al. [13,14] identified a total of 146 volatiles in 25 mandarin hybrids and reported several new aroma active compounds by gas chromatography-olfactometry (GCO).
Flavor improvement is highly prioritized in recent mandarin breeding programs. The molecular markers linked to mandarin fruit aroma could facilitate genetic improvement and the development and release of new mandarin cultivars with improved fruit flavor. Compared to crop plants and vegetables, genetic mapping of fruit aroma in citrus, especially in mandarin, has not been performed. We recently constructed mandarin genetic linkage maps and identified quantitative trait loci (QTLs) for fruit quality characteristics [18]. In the present study, we investigated aroma volatiles in a mandarin F 1 population using GCMS, and then we identified QTL regions controlling mandarin aroma volatiles. The long term objective is to understand characteristics of "good" mandarin fruit and to find fruit quality molecular markers for marker-assisted selection (MAS).

Plant materials
Mandarins (Citrus reticulata Blanco) 'Fortune' , 'Murcott' and their 116 F 1 hybrids, one tree each, were planted in 1991 for breeding purposes. Fruits were obtained from the University of Florida Citrus Research and Education Center groves in Lake Alfred, FL. Fruit of 'Fortune' , 'Murcott' and F 1 hybrids were harvested twice (a month apart) during their optimal commercial maturity in 2012 and 2013, respectively. The four harvests were presented as H1-12, H2-12, H1-13 and H2-13 in the main text subsequently.

Sample preparation
Fruits were hand-harvested in the grove and transferred to the lab within 1-2 h for phenotypic evaluation. The phenotypic evaluation was performed as described [18]: 15 sound fruits from each tree were used to measure flavedo and juice color, fruit diameter (FD, in mm) and fruit weight (FW, in g) individually. Flavedo and juice color were measured using a Minolta CR-331 colorimeter. Color measurement score was described as Hunter color space value: L (black to white), a (green to red), b (blue to yellow), and a over b ratio. Then the 15 fruits were grouped into three pools with five fruits per pool, serving as three biological replicates. For the two harvests in 2012, the flavedo layer was first collected for another experiment and then the fruits were manually peeled and juiced (pressed) using a Hamilton Beach 932 manual commercial citrus juicer (Hamilton Beach, Southern Pines, NC, USA). For the two harvests in 2013, mandarin fruits were manually peeled and then juiced directly using the same procedure and juicer. Juice percentage (JP, in ml/100 g) was calculated, then color, soluble solids content (SSC, in g/100 mL), titratable acidity (TA, in g/100 mL, presented as percentage of citric acid) and volatiles were measured.

Volatile extraction and GCMS analysis
A mixture of 3 mL of juice, 3 mL of saturated sodium chloride solution (359 g/L), and 6 μL of 1000 ppm 3hexanone as internal standard were placed in a 20 ml glass vial and capped with a magnetic crimp cap with a silicone/PTFE septum (Gerstel, Linthicum, MD, USA). The vials were stored at −20°C until analyzed. Frozen vials were thawed under tap water and loaded into the autosampler (Model MPS2, Gerstel Inc., Linthicum, MD) equipped with a cooled tray holder (Laird Tech, Sweden) controlled by a Peltier Thermostat (CTC Analytics AG, Switzerland). Samples were held 2-20 h at 4°C until analyzed. For analysis, juice samples were incubated for 30 min at 40°C. A 2.0 cm solid phase microextraction (SPME) fiber (50/30 μm DVB/Carboxen/PDMS; Supelco, Bellefonte, PA) was then exposed to the headspace for 30 min at 40°C. After exposure, the SPME fiber was inserted into the injector of a GCMS (Model 6890, Agilent, Santa Clara, CA) to desorb the extract for 15 min at 250°C. The GCMS equipment and settings were: DB-5 columns (60 m length, 0.25 mm i.d., 1.00 μm film thickness; J&W Scientific, Folsom, CA), coupled with a 5973 N MS detector (Agilent Technologies, Santa Clara, CA). The column oven was programmed to increase at 4°C min −1 from the initial 40°C to 230°C, then ramped at 100°C min −1 to 260°C and held for 11.70 min for a total run time of 60 min. Helium was used as carrier gas at a flow rate of 1.5 mL min −1 . Inlet, ionizing source and transfer line were kept at 250, 230, and 280°C, respectively. Mass units were monitored from 30 to 250 m/z and ionized at 70 eV. Data were collected using the ChemStation G1701 AA data system (Hewlett-Packard, Palo Alto, CA). A mixture of C-5 to C-18 n-alkanes was run at the beginning of each day to calculate retention indices (RIs) [19].

Volatile compound identification
Volatile compounds were identified by comparing their mass spectra with the Wiley HPCH 2205.L (Hoboken, NJ, USA), NIST05.L (NIST/EPA/NIH Mass Spectral Library; National Institute of Standards and Technology, Gaithersburg, MA, USA), and ADAMS.L [20] mass spectral database, as well as comparing RIs with published RIs on the same column. The amount of each aroma volatile was expressed as relative content (aroma volatile peak area over internal standard peak area).

SNP genotyping, statistical analysis and QTL mapping
Mandarins were genotyped using an Illumina GoldenGate 1536-SNP array platform at the University of Florida's Interdisciplinary Center for Biotechnology Research. The 1536 evenly distributed SNPs were initially mined from sweet orange bacterial artificial chromosome (BAC) end sequences [21], and have been applied to construct mandarin genetic linkage maps and identify fruit quality related QTLs using the same 'Fortune' x 'Murcott' population [18], and were then used for the QTL identification of mandarin juice volatiles in the published genetic linkage maps.
The descriptive statistical analysis of volatiles was performed using JMP Pro 11. T-test was applied to compare the means of each volatile between parents. In addition, Pearson correlation coefficients of each volatile for F 1 hybrids were calculated. The analysis of variance (ANOVA) test was performed to test for differences of volatiles among the citrus selections. Due to nonnormality, all volatile data were Log 10 -transformed and used in the subsequent statistical and QTL analysis. Principal component analysis (PCA) based on Pearson's correlations was applied to differentiate tree individuals based on their volatile content using XLSTAT software (Addinsoft, Paris, France). Finally, the corresponding heat map representation of pairwise correlations was performed with JMP Pro 11. The QTL identification was performed on the published mandarin genetic linkage maps generated from the same population as described previously [18]. The non-transformed data were analyzed first by the Kruskal-Wallis rank-sum test using MapQTL 6.0, a significance level of p = 0.001 was used as threshold. Then Log 10 -transformed data were used in the subsequent QTL analysis using JMP Genomics 7 on the parental maps separately. QTLs were detected using composite interval mapping (CIM) [22,23] with the expectation maximization (EM) mapping algorithm [24]. Logarithm (base) of odds (LOD) thresholds were calculated using a 1000-permutation test for each volatile on each map [25]. A QTL test step of 2 cM was used for CIM. QTLs were reported when their LOD values exceeded the LOD threshold, and they were literally named using volatile code followed by the number of the linkage group (LG) in which the QTL was located.

Candidate gene screening
With colocation of SNPs on both genetic maps and the Clementine genome, the detected QTLs were compared to the map position of aroma volatile-related genes annotated in the citrus Clementine pathways from the CitrusCyc Pathways Database (http://pathways.citrusgenomedb.org, part of the Citrus Genome Database) to identify potential candidate genes. All the annotated genes related to the production of corresponding volatiles were considered in the candidate screening process, for example, valencene synthase (TPS1), geranyl diphosphate synthase 1 (GPS1), lipoxygenases (LOX), hydroperoxide lyase (HPL), hexenal isomerases (HI) and alcohol dehydrogenases (ADH). Candidate genes within a QTL/ QTL cluster were searched only in the Clementine reference genome region delimited by the confidence interval of the QTL using the most proximal SNPs that were both on the Clementine genome sequence (http:// www.phytozome.net/clementine.php) and the genetic linkage maps. Whether a gene was chosen as a candidate was based on the predicted annotation of the gene from the CitrusCyc Pathways Database, and was further confirmed by BLAST of the gene sequence from the Clementine genome against the Arabidopsis Information Resource (https://www.arabidopsis.org/).
The average number of volatiles detected in F 1 hybrids was 63.7, and ranged from 28 to 111 (Table 1). The average relative content of volatiles detected in F 1 hybrids was 14.4, and ranged from 0.24 to 133.24. Three volatiles were present in all F 1 hybrids, including linalool, myrcene and limonene; and another 18 volatiles were detected in 90-99% of all F 1 hybrids (Additional file 1). On the other hand, 19 compounds were detected in less than 10% of all hybrids. Among the 148 detected     volatiles, 64 volatile compounds were detected specifically in F 1 hybrids but were not present in either of the parents. The total content of volatiles and number of volatiles in F 1 progeny showed continual variation, although some of their distributions were generally skewed and/or significantly deviated from normality (p = 0.05) based on the Shapiro-Wilk test [26] (Fig. 1, Additional file 3). For instance, the total content displayed a distribution significantly deviated from normality in H1-12, so did the number of volatiles in H1-12, H1-13 and H2-13. Transgressive segregations were found in both directions where decanal, perillaldehyde, ethyl acetate, octyl acetate, decyl acetate, pmenth-1-en-9-ol acetate, α-pinene, myrcene, limonene and terpinolene exhibited the most significant transgressions (Table 1, Additional file 2). A total of 20 volatile compounds were reported with odor activities in the F 1 individual FoMu-081 [14]. In the present study (Table 1, Additional file 2), (E)-2-decenal described as minty/piney, 1-octen-3-one described as mushroom and β-ionone described as floral, were found in 'Murcott' , but not in 'Fortune' in 2013. Four aroma active compounds were detected only in F 1 hybrids in 2013, including (E)-carveol described as minty, dodecanal described as green, ethyl 2-methylbutanoate described as fruity and nootkatone described as spicy/ woody. Two aroma active compounds were present in all F 1 hybrids in 2013, including linalool described as floral and myrcene described as metallic.

Correlation analysis of volatiles in
The pair-wise correlation coefficient for each volatile against every other compound was calculated in order to identify co-regulated compounds. Pearson correlation was performed for F 1 progeny juice samples in H1-13. Of the 11,476 pairs analyzed, 5885 resulted in significant correlations (p < 0.05) (data not shown). Of these 5885 pairs, most of them (5724) exhibited positive correlation coefficients. The strongest positive correlations were found between α-muurolene with δ-cadinene (r = 0.999) and (E)-cadina-1,4-diene (r = 0.998). Strong positive (r > 0.5) correlations were also found between (E)-2octenal with (E)-2-nonenal (r = 0.997) and heptanal (r = 0.996), between δ-cadinene and (E)-cadina-1,4-diene (r = 0.996), between nonyl acetate with octyl acetate (r = 0.995) and decyl acetate (r = 0.994). The strongest The 148 volatile compounds were categorized into 12 groups based on their chemical structure including acid, alcohol, aldehyde, ester, ether, furan, aromatic hydrocarbon, ketone, phenol, sesquiterpene, monoterpene and unknown [13]. The pair-wise correlation coefficient for each chemical group against every other group was calculated for juice samples from hybrids (Additional file 4). Strong positive correlations were observed for most of the major groups over four harvests, except for the ketone group which showed poor or no correlations. The ketone group was correlated with alcohols (r = 0.49 and 0.56) in 2013, but showed no strong correlations with any other groups.
The pair-wise correlation coefficient for each chemical group against every fruit characteristic was also calculated for hybrids (Additional file 5), and few strong correlations (|r| > 0.5) were observed. In 2013, total volatile content positively correlated with the traits TA and juice color space value L, but negatively correlated with fruit diameter, fruit weight, juice percentage, SSC/ TA ratio, flavedo and juice color space value a and a/b ratio. Significant (p < 0.05) positive correlations were found between TA with alcohols, aldehydes, esters and monoterpenes over two harvests in 2013. Juice color space value a was negatively correlated with monoterpenes over four harvests, and juice color space value b displayed relatively strong correlations with alcohols, aldehydes, esters and monoterpene in H1-13. The pairwise correlation coefficient for color characteristic against every carotenoid-derived volatile including 6methyl-5-hepten-2-one, β-cyclocitral, isophorone, neryl acetate, geranyl acetone and β-ionone, was calculated for hybrids (Additional file 6). Only a few stable significant correlations were found. Flavedo and juice color space value a and a/b ratio displayed stable significant negative correlations with β-cyclocitral and neryl acetate, and a positive correlation with geranyl acetone over two or more harvests.
Cluster validation using correlation analysis HCA (hierarchical cluster analysis) based on pairwise correlation coefficients were run individually for four harvests, and four distinct groups of volatiles were formed in H1-13 (Additional file 7). Cluster A included 26 volatiles, of which the majority were aldehydes, ketones and unknown compounds. Cluster B contained two sub-clusters B1 and B2 with 24 and 27 volatiles respectively. The B1 sub-cluster consisted of mainly aldehydes, esters, monoterpenes and unknown compounds, and the B2 comprised mainly alcohols, sesquiterpenes and monoterpenes. Cluster C included two sub-clusters C1 and C2 with 18 and 12 volatile compounds respectively. The C1 sub-cluster contained mainly sesquiterpenes, and the C2 comprised mainly sesquiterpenes and unknown compounds. Cluster D included two subclusters D1 and D2 with 20 and 23 volatiles respectively. The D1 consisted of mainly sesquiterpenes and unknown compounds, and the D2 contained mainly aldehydes, esters and ketones. The numbers can be found in the second column (volatile code) preceded by Pk in Additional files 1 and 2 matched a 3:1 ratio. All QTLs for volatile compounds were identified using Log 10 transformed data. A total of 244 significant associations were found between molecular markers and 94 volatiles ( Fig. 3 and Additional file 8).
QTLs for each volatile located in the same chromosomal regions from different harvest seasons were considered to be the same, hence the 244 associations then were summarized into 206 QTLs, with 25 (12.14%) being repeatable ( QTLs were grouped into four and three main clusters in the FOR and MUR maps, respectively. The first cluster on the FOR map was located on FOR1.1, where QTLs were identified for 47 volatiles including six alcohols, six aldehydes, five esters, one phenol, 11 sesquiterpenes, 10 monoterpenes and eight unknown compounds. Out of the 47 volatiles, nine volatiles were aroma active compounds, including linalool, heptanal, (E)-2-heptenal, decanal, dodecanal, myrcene, β-phellandrene, γ-terpinene and terpinolene. The second and third clusters were both located on FOR3.3. The second cluster was made up of nine QTLs identified for one alcohol, one ketone and seven sesquiterpenes. The R 2 was 59.82% for valencene in H1-13. Within the confidence intervals of the cluster, no candidate genes were found related to terpenoid regulation, even though the gene encoding valencene synthase (TPS1) was located on the Clementine genome scaffold 3. The third cluster of QTLs was identified for one acid, one alcohol, one ester, one ketone, two sesquiterpenes, eight monoterpenes and one unknown compound. The fourth cluster was located on FOR6.3, included QTLs for one alcohol, one aldehyde, one ester, one sesquiterpene and eight monoterpenes. Five of the eight monoterpenes were aroma active compounds including linalool, myrcene, terpinolene, β-phellandrene and γ-terpinene.
The first cluster on the MUR map was located on MUR1.2, and contained QTLs for 24 volatiles including one acid, two alcohols, three esters, three ketones, six sesquiterpenes and nine monoterpenes. The second cluster was located on MUR2, consisted of QTLs for 22 volatiles including one alcohol, one aldehyde, five esters, five sesquiterpenes, nine monoterpenes and one unknown compound. Four of them were aroma active compounds including myrcene, β-phellandrene, γterpinene and dodecanal. Within the confidence intervals of the QTL cluster on MUR2, 572 annotated genes were listed. After screening, four of these genes, GPS1, TPS3, TPS4 and TPS14, encoding geranyl diphosphate synthase 1, terpene synthase 3, terpene synthase 4 and terpene synthase 14 respectively, were probably related to terpenoid regulation based on their biological function as studied in model plant species such as Arabidopsis (Additional file 9). The third cluster was located on MUR4.1, where QTLs were identified for two esters, one aromatic hydrocarbon and one unknown compound.

Mandarin juice aroma volatile profiles
During the sample preparation, any scraping of the albedo or squeezing of the flavedo may contaminate juice samples with peel components and change the aroma volatile profiles due to high concentration of volatile compounds in mandarin peel oil. In this study, mandarin fruits were manually peeled and juiced using a juicer. For the two harvests in 2012, the flavedo layer was first collected for another experiment and then the fruits were peeled. In order to minimize the error and exhibit the truest mandarin juice volatile composition, only aroma volatiles presented in all four harvests were reported here, and the QTLs only identified in 2012 were discarded. In this study, 148 volatiles were detected by GCMS in juice samples prepared from manually peeled fruits of 'Fortune' , 'Murcott' and their 116 F 1 hybrids. A total of 55 and 48 volatiles were detected for 'Fortune' and 'Murcott' respectively in 2013. Moshonas & Shaw [27] detected and quantified 42 juice volatile compounds from unpeeled fruits of 15 mandarin selections, and 34 volatiles from 'Murcott'. In juice samples prepared from unpeeled fruits, 82 and 80 volatiles were reported from 'Murcott' and one hybrid of 'Fortune' x 'Murcott' respectively [13]. In the present study, the number of volatiles detected in F 1 hybrids ranged from 28 to 111 volatiles with average number of 63.7 in H1-13. Kerbiriou et al. [10] studied aroma volatiles in juice samples extracted from unpeeled fruits of 56 mandarin   The percentage of the total phenotypic variation explained by the QTL hybrids and detected 225 volatile compounds, ranging from 21 to 86 volatiles, and the subsequent study of 25 genetically related mandarin materials identified 146 juice aroma volatiles from unpeeled fruits, and the range for each sample was 52 to 118, with an average number of 77 per sample [13]. Limonene, ethanol, α-pinene and cymene presented in all juice samples from 56 mandarin hybrids [10], while Miyazaki et al. [13] found 15 consensus aroma volatiles from 25 mandarin hybrids. Tietel et al. [3] identified 37 consensus aroma volatiles in mandarin juice present in at least four previous studies and considered nine volatile compounds as core aroma volatiles including linalool, α-terpineol, terpinen-4-ol, nonanal, decanal, carvone, limonene, α-pinene and myrcene. Here, limonene, myrcene and linalool were present in all F 1 hybrids of 'Fortune' x 'Murcott' , and another 18 volatile compounds were detected in 90-99% of all F 1 hybrids, such as αpinene, ethanol, hexanal, α-phellandrene and terpinolene (Additional file 1). Conversely, 19 volatiles were detected in less than 10% of all F 1 hybrids and may be considered as hybrid specific volatile compounds, such as (E)-carveol and nootkatone presented in only nine hybrids.
'Fortune' was reported to be derived from a 'Clementine' x 'Orlando' tangelo [28], while 'Murcott' was reported as a mandarin sweet orange hybrid [29]. The different parentages produced significant differences in aroma volatile composition and content, and sesquiterpene was the most important contributor. A total of 13 sesquiterpenes were identified in 'Fortune' , however, only one sesquiterpene was found in 'Murcott' in 2013. The detection of fewer sesquiterpenes in 'Murcott' relative to other mandarin selections was reported, where four sesquiterpenes were identified in 'Murcott' , compared to eight, 19, 23 and 20 in 'Fallglo' , 'Temple' , 'Sanguinelli' and 'Ortanique' mandarins [13,16]. Compared to 'Fortune' , 'Murcott' produced a similar number of volatiles, but with higher concentrations. This agrees with the previous report that 'Murcott' juice samples contained a higher level of many water-soluble volatiles compared to other mandarins, especially alcohols with one to five carbons [27].

Mandarin juice aroma volatile correlations and segregation
Volatile compounds from the same biochemical pathway were normally clustered together, which indicated a coregulation of these chemical compounds. In agreement with a study in strawberry [30], most of the detected volatile compounds exhibited positive correlations. The correlations are consistent with the existence of clustering of QTLs for volatile compounds from the same chemical family and few QTL intervals controlling the content of most of the detected volatile compounds. Clustering of QTLs for volatiles with the same structures was observed in apple [31], tomato [32], almond [33], strawberry [30], peach [34] and Arabidopsis [35]. The observed QTL clusters may correspond to a tight linkage of different loci or pleiotropic effects of a single locus on several volatiles. A single locus could affect an entire pathway by either encoding transcription factors to regulate genes coordinately or encoding enzymes catalyzing limiting steps in single pathways.
Consistent correlations were found between volatile content and fruit color in the 'Fortune' x 'Murcott' hybrids. Juice color space value a negatively correlated with monoterpene content over four harvests, while juice color space value b displayed relatively strong positive correlations with alcohol, aldehyde, ester, aromatic hydrocarbon and monoterpene in H1-13. Strong associations between color and aroma volatiles were observed in watermelon and tomato [36], and the relationship is likely to be a consequence of the degradation of carotenoids into volatile compounds. The norisoprene and monoterpene volatile compounds were influenced by carotenoid pigmentation, and the volatile geranial was apparently derived from the carotenoid lycopene [37]. Carotenoid contents in tomato fruit were associated directly with the emissions of carotenoid derived volatiles [38]. Higher production of apocarotenoid volatiles and lower concentration of carotenoids were found in 'Temple' than 'Murcott' fruits [16]. In our study, three of the six carotenoid derived volatile compounds, βcyclocitral, neryl acetate, and geranyl acetone, consistently correlated with flavedo and juice color space value a and a/b ratio.
Some volatiles displayed qualitative variations and matched single locus inheritance. For example, valencene was present in only one parental line, 'Fortune' , and in about 50% of the progeny, which matched the 1:1 ratio in 2013. In addition, a large-effect QTL was identified for valencene, with R 2 being up to 59.82%. In strawberry, mesifurane and γ-decalactone were mapped as single Mendelian traits [30]. Spiller et al. [39] mapped nerol and neryl acetate as single Mendelian traits and found that geranyl acetate was controlled by two independent loci in diploid roses. It was observed in grape [40] and kiwifruit [41] that concentration of different terpenes was controlled by major genes. Transgression, which results in progeny from new combinations of multiple genes with more positive or more negative effects on a quantitative trait than were present in either parent, is often observed in the progeny derived from interspecific mating [42]. The wide variation in each volatile compound among the mandarin hybrids could be exploited to generate recombinants with markedly improved mandarin aroma flavor. Most of the detected compounds in mandarin F 1 hybrids displayed transgressive segregations, and such strong transgressive segregation for volatile compounds was also observed in apple [31], peach [34], strawberry [30,43] and tomato [32].
Terpenoids are the largest class of secondary metabolites in plants and play important roles in both plant and human health. In plants, isopentenyl diphosphate (IPP) and dimethylallyl diphosphate (DMAPP), the two C5 building blocks for all terpenoids, are derived from two independent pathways, the methylerythritol phosphate (MEP) pathway localized in the plastids [44] and the mevalonate (MVA) pathway localized in the cytosol [45]. The MEP pathway provides IPP and DMAPP for production of hemiterpene, monoterpene, and diterpene, whereas the MVA pathway provides IPP and DMAPP for sesquiterpene biosynthesis. In our study, 18 out of 27 sesquiterpenes detected in this study were grouped into cluster C based on correlation data, and the 41 QTLs identified for 18 sesquitepenes were mapped to four genomic regions. The overlapping of QTLs and strong correlations among sesquiterpenes suggested that in this mapping population the concentrations of several sesquiterpenes were controlled by a few distinct loci with pleiotropic effects, in accordance with the fact that all sesquiterpenes originate from similar biosynthetic pathways. Genes either affect the early stages of the MVA pathway or the early steps of terpene biosynthesis to have a strong effect on the numerous sesquiterpenes, or else they affect later steps in the terpene synthetic pathway to influence relatively few terpenes. Here, 11 out of 15 monoterpenes detected were grouped into cluster B based on their correlation coefficients. The 50 QTLs identified for 14 monoterpenes mapped to five genomic regions. The overlapping of QTLs and strong correlations among monoterpenes indicated pleiotropic effects of a few distinct loci. The monoterpenes and sesquiterpenes are believed to originate from quite different pathways, however, cross talk between the MEP and MVA pathways was reported in Arabidopsis [46]. In snapdragon flowers, only the plastid-localized MEP pathway was active in producing terpene compounds, providing IPP for both synthesis of monoterpenes in plastid and sesquiterpenes in cytosol [47]. In our study, total content of sesquiterpenes was positively correlated with total content of monoterpenes consistently over four harvests. The overlap of QTL clusters on FOR1.1, FOR3.4, MUR1.2 and MUR2 for both monoterpenes and sesquiterpenes may indicate tight linkage of distinct loci from MEP and MVA pathways, or the presence of major regulators in those loci, or a few genes affecting the early steps in one of the two pathways further influencing concentration of both monoterpenes and sesquiterpenes.

Candidate genes for some QTLs
In this study, 12.14% of the QTLs for volatile compounds were stable over two or more harvests, which was less than that of aroma studies in peach [34] and strawberry [30], where 67.44 and 50% of QTLs were reproducible over two or more seasons, respectively. All genes related to the production of corresponding volatiles were considered and added to the candidate screening process, and four candidates (GPS1, TPS3, TPS4 and TPS14) were eventually identified in the corresponding map regions. The terpenoids are the largest and most diverse class of plant secondary metabolites with many aroma volatiles, including hemiterpenes (C5), monoterpenes (C10), sesquiterpenes (C15), homoterpenes (C11 and C16) and diterpenes (C20) [48]. In our mandarin aroma study, we detected 27 sesquiterpenes and 15 monoterpenes, some of them presented in more than 90% of all F 1 hybrids. These could be considered as consensus volatile compounds, such as α-pinene, myrcene, limonene, α-phellandrene, α-terpinene, p-cymene, βphellandrene, terpinolene and p-mentha-2,4(8)-diene. A total of 41 and 50 QTLs were identified for 18 sesquitepenes and 14 monoterpenes respectively. A QTL interval on MUR2 was identified for one alcohol, one aldehyde, five esters, five sesquiterpenes, nine monoterpenes and one unknown compound. The QTL interval corresponded to a single genomic region that contains genes GPS1, TPS3, TPS4 and TPS14, encoding geranyl diphosphate synthase 1 (GPS1), terpene synthase 3 (TPS3), terpene synthase 4 (TPS4) and terpene synthase 14 (TPS14) respectively. GPS is responsible for generating the C10 terpene precursor geranyl pyrophosphate (GPP) by joining one IPP and DMAPP unit in the MEP pathway. Monoterpenes and sesquiterpenes are formed from terpene precursors GPP and farnesyl diphosphate (FPP) respectively catalyzed by the terpene synthases [49]. A total of 49 putative members of the terpene synthase family were reported in Citrus and grouped into six subfamilies [50]. GPP was converted to monoterpene (E)-βocimene catalyzed by a TPS3-encoded recombinant enzyme in Arabidopsis [51]. The allelic variation of terpene synthase genes, TPS2 and TPS3, produced the different level of emissions of (E)-β-ocimene and (E,E)-α-farnesene between two Arabidopsis accessions 'Wassilewskija' and 'Columbia' [52]. In Arabidopsis, caterpillar herbivory (Pieris rapae) induced the expression of the β-ocimene synthase TPS3 with increased myrcene emission [53]. TPS3 was responsible for the formation of camphene and tricyclene in tomato [54]. The monoterpene synthase TPS4 catalyzed the formation of mostly β-phellandrene from GPP in tomato [55]. The gene TPS14 expressed in Arabidopsis flowers and encoded enzymes producing linalool [56]. TPS14 catalyzed the production of several bisabolene isomers and nerolidol from FPP in vitro [54]. GPS activity was reported in tissues producing abundant quantities of monoterpenes in plant species sage (Salvia officinalis) [57], Pelargonium roseum [58], Vitis vinifera L. cv Muscat de Frontignan [59], and Abies grandis [60]. GPS existed as a heterodimer with small and large subunits in Antirrhinum majus flowers, GPS small subunit mRNA expression level and protein were correlated with monoterpene biosynthesis [61].

Conclusions
The aroma volatiles were investigated in mandarin parents 'Fortune' , 'Murcott' and their derived F 1 individuals using GCMS in 2012 and 2013 harvest seasons. A total of 148 volatiles were identified by the mass spectrometry, and among them, 28, 18 and 66 unique volatile compounds were detected in 'Fortune' , 'Murcott' and the F 1 hybrids respectively in 2013. Sesquiterpenes were the most important contributors to the difference between 'Fortune' and 'Murcott' , as 13 sesquiterpenes were identified in 'Fortune' while only one sesquiterpene was found in 'Murcott' in 2013. A total of 206 QTLs were identified for 94 volatile compounds including 17 aroma active compounds, among them, 25 (12.14%) were consistent over two or more harvests. On MUR2, a QTL interval controlling multiple monoterpenes and sesquiterpenes corresponded to a genomic region that contains genes encoding geranyl diphosphate synthase 1 (GPS1), terpene synthase 3 (TPS3), terpene synthase 4 (TPS4) and terpene synthase 14 (TPS14). All 206 QTLs were validated in 13 citrus selections in 2014 harvest season, and validations on broader breeding germplasm and segregating populations are in progress, with some of them showing the potential for MAS in citrus breeding programs (data not shown). These SNPs could be used to screen mandarin individuals for some specific volatile compounds or even aroma active compounds after further validation. However, little is known regarding the contribution of individual volatile compounds, or about interactions between volatile compounds and the juice matrix, and other volatile compounds to the whole mandarin sensory profile. Further research is needed to understand which volatile compounds constitute the 'desirable' and 'undesirable' aroma for most of the consumers. The volatile compound composition should be correlated with the mandarin sensory profile generated by trained taste panels and consumer surveys, in order to understand the possible role of each individual volatile compound in the matrix for the mandarin sensory experience. The genetic linkage map should be combined with the mandarin sensory profile from trained panel and consumer surveys, in order to identify any QTLs linked to different mandarin sensory descriptors, and 'desirable' or 'undesirable' aroma volatiles for consumers. Once the QTLs are identified for the main mandarin sensory contributors, and for the 'desirable' or 'undesirable' aroma volatiles to most consumers, the molecular marker-assisted approach for mandarin fruit quality improvement may be utilized.