{'fraction_non_alphanumeric': 0.037772353851649, 'fraction_numerical': 0.017263959995237527, 'mean_word_length': 5.060064935064935, 'pattern_counts': {'":': 0, '<': 5, '<?xml version=': 0, '>': 5, 'https://': 1, 'lorem ipsum': 0, 'www.': 2, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Regional Variation of Chemical Characteristics in Young Marselan (Vitis vinifera L.) Red Wines from Five Regions of China

The environmental conditions of wine regions determine the flavor characteristics of wine. The characterization of the chemical composition and sensory profiles of young Marselan wines from five wine-producing regions in China was investigated by gas chromatography-mass spectrometry (MS), high-performance liquid chromatography–triple-quadrupole MS/MS and descriptive analysis. Young Marselan wines can be successful discriminated based on concentrations of volatile compounds, but not phenolic compounds, by orthogonal partial least squares discriminant analysis according to regions. Compared to Jiaodong Peninsula (JDP) and Bohai Bay (BHB) regions, there were relatively lower average concentrations of varietal volatiles (mainly including β-citronellol, geraniol, and (E)-β-damasenone) and several fermentation aroma compounds (including isoamyl acetate, octanoic acid, decanoic acid, ethyl decanoate, etc.) but higher levels of acetic acid in Xinjiang (XJ), Loess Plateau (LP), and Huaizhuo Basin (HZB) regions, which were related to their characteristic environmental conditions. Marselan wines from HZB, LP, and XJ regions were characterized by lower L values and higher a and Cab values. Marselan wines from XJ were discriminated from the wines from other regions due to their higher concentrations of several flavonols. Sensory analysis indicated that Marselan wines from HZB region were characterized by relatively low intensities of floral and fruity aromas compared to other regions.


Introduction
Marselan (Vitis vinifera L.) is a particularly successful cross of Cabernet Sauvignon and Grenache Noir obtained by Institut national de la recherche agronomique (INRA) in 1961 ( Figure S1). This variety shows good resistance to powdery mildew, mites, and especially to botrytis bunch rot [1]. Marselan was first introduced into China in 2001 and was planted in the Huaizhuo Basin region, Hebei Province. Due to its strong disease resistance and small berries with high quality, Marselan has been popular and planted almost wherever vines are grown for wine production in China, except for cool wine-producing regions. Followed by France, China has the second-largest planting area of Marselan [2], and the planting area is still growing.
During the past decade, several researchers have focused on the phenolic compounds in Marselan grapes, such as the phenolic compounds in skins [3] and changes in the different degrees of polymerization in seeds [4]. In addition, effects of some viticultural practices on the chemical compositions in Marselan grape berries were studied, such as grape cluster shading [5] and flower cluster tip removal [6], and influences of several enological parameters on Marselan wines were investigated, such as application of non-Saccharomyces [7] and A total of 38 commercial Marselan red wines from five important wine-producing regions, namely Jiaodong Peninsula (JP), Bohai Bay (BHB), Huaizhuo Basin (HZB), Loess Plateau (LP), and Xinjiang (XJ), were used in this study (Table S1). The wines ranged in vintage from 2012 to 2016, and their accurate subregions are presented in Table S1 and plotted on a China map ( Figure S2). All were produced as premium wines, which go for a period of bottle aging. The basic wine compositions (including alcohol, reducing sugar, total acidity, pH, and glycerol) of the Marselan wines were measured by a WineScan (FT 120) rapid-scanning infrared Fourier-transform spectrometer with FOSS WineScan software version 2.2.1 (Foss Electric, Hillerød, Denmark). Wines samples were centrifuged at 4000× g for 8 min and the analyzed following the manufacturer's instructions. Colorimetric measurement, analysis of volatile compounds and phenolic compounds, and descriptive analysis were carried out on these wines and completed at the end of 2017.

Meteorological Data
The meteorological data (except for the photosynthetically active radiation (PAR)) of these five wine regions, from 2000 to 2017 vintages, were obtained from the National Meteorological Information center [15]. The PAR data of these wine regions, from 2000 to 2014 vintages, were obtained from an open access dataset, named "A dataset of reconstructed photosynthetically active radiation (PAR) in China (1961-2014)" [16]. In this study, the sunshine duration, PAR, precipitation, average daily temperature, and diurnal temperature range in the growing season (from May to October) for each wine region were calculated.

Colorimetric Measurement
Colorimetric measurement was carried out on a UV-visible spectrophotometer (Shimadzu UV-2450, Shimadzu Co., Kyoto, Japan), recording the wine absorbance spectra (380-700 nm). Ultra-pure water was used as blank. The values of the CIELab parameters, including chromaticity (Cab), lightness (L), hue (Hab), red-greenness (a), and yellow-blueness (b), were calculated according to a previous reported method [17]. Prior to analyses, all wine samples were filtered through 0.45 µm cellulose filters (Membrana Co., Wuppertal, Germany) and placed in a 2 mm path length quartz cuvette. All the experiments were performed in triplicate.

HS-SPME-GC-MS Analysis of Aroma Compounds
The extraction and identification of the volatile compounds were carried out according to the analytical method previously reported in our laboratory, with slight modifications [18]. The analysis of the volatile compounds in wines was performed on an Agilent 6890 chromatograph (GC) coupled to an Agilent 5975C mass spectrometer (MS, Agilent Technologies, Santa Clara, CA, USA). A CTC CombiPAL autosampler (CTC Analytics AG, Zwingen, Switzerland) with a 50/30 µm divinylbenzene/carboxen/polydimethylsiloxane (DVB/CAR/PDMS) fiber (2 cm, Supelco, Bellefonte, PA, USA) was used for the extraction of the volatile compounds. A wine sample of 5 mL was transferred into a 20 mL vial with 10 µL of 4-methyl-2-pentanol (internal standard, 1.0086 g/L) and 1 g NaCl. Tightly capped with a PTFE-silicon septum, the samples were equilibrated at 40 • C for 30 min before extraction by an SPME fiber at 40 • C for 30 min with stirring at 500 rpm. After extraction, the fiber was thermally desorbed by insertion into the injection port of the GC for 8 min. The injection was performed in splitless mode (0.75 min) with a GC inlet temperature of 250 • C.
The separation of the volatile compounds was carried out on a HP-INNOWAX capillary column (60 m × 0.25 mm × 0.25 µm, J&W Scientific, Folsom, CA, USA). The flow rate of the carrier gas (helium) was 1 mL/min. The oven temperature was held at 50 • C for 1 min after injection, then programmed to 220 • C at a rate of 3 • C/min, and kept at 220 • C for 5 min. The temperatures of the mass selective detector transfer line, ion source, and quadrupole were 250 • C, 230 • C, and 150 • C, respectively. The positive ion electron impact spectra at 70 eV were recorded with a mass range of 30-350 m/z for the full scan mode.
The quantification process was carried out according to our previous study [18]. A synthetic model wine solution was prepared in distilled water containing 14% v/v ethanol, 5 g/L tartaric acid, and the pH was adjusted to 3.3 with a 5 M NaOH solution. The calibration curves of each aroma standard were obtained according to fifteen dilution levels in succession with the model wine solution.

HPLC-QqQ-MS/MS Analysis of Phenolic Compounds
The analyses of the phenolic compounds were performed on an Agilent 1200 series high-performance liquid chromatographer (HPLC) equipped with an Agilent 6410B triplequadrupole (QqQ) mass spectrometer (MS) (Agilent Technologies, CA, USA), based on our previous studies [19,20]. The chromatographic separation was carried out using a Poroshell 120 EC-C18 column (150 mm × 2.1 mm, 2.7 µm; Agilent Technologies, Santa Clara, CA, USA). The mobile phases were: A = water solution containing 0.1% (v/v) formic acid; B = acetonitrile/methanol solution (50:50, v/v) containing 0.1% (v/v) formic acid. The gradient elution was: (1) from 10% to 46% B in 28 min, and (2) from 46% to 10% B in 1 min. The post time was 5 min. A wine sample of 1 mL was filtered using 0.45 µm inorganic membranes (polyether sulphone). The injection volume was 5 µL and the flow rate was 0.4 mL/min. The column was thermostatically controlled at 55 • C. An electrospray ionization source was used with 4 kV voltage in the negative mode and the positive mode for the non-anthocyanin phenolics and the anthocyanin compounds, respectively. The temperatures of the ion source and gas (N 2 ) were 150 • C and 350 • C, respectively. The gas had a flow rate of 12 L/min, while the nebulizer pressure was 35 psi. The [M + H]+ and [M + H]− ions were selected as the precursors for the anthocyanin and non-anthocyanin compounds, respectively. The multiple reaction monitoring mode was used for both identification and quantification. The identification was achieved by comparing the retention times and qualitative transition ions with the phenolic compounds in an in-home phenolic compound library [19,20]. The quantification process was carried out according to reliable methods established in our previous studies [19,20]. All anthocyanins were quantified on the basis of the calibration curve of malvidin-3-O-glucoside and non-anthocyanins were quantified according to the calibration curves of their own reference compounds.

Descriptive Analysis
Descriptive analysis was conducted in a sensory laboratory equipped with 20 individual booths at controlled room temperature (20 • C). Approximately 30 mL of each of the wines was prepared in International Standards Organization (ISO) wine tasting glasses (ISO 3591:1977). The glasses were labeled with three-digit random numbers and served in a randomized order.
A panel consisting of 16 assessors (8 males and 8 females) was convened for sensory evaluation of Marselan wines from different wine-producing regions. All panelists were student or staff or faculty at the Center for Viticulture and Enology (CFVE), China Agricultural University. All panelists attended four one-hour training sessions. In the first session, all panelists were asked to evaluate all Marselan wines used in this study and generate descriptive terms, including 21 aroma and taste attributes. In the second session, panelists were asked to discuss the terms and seven attributes (hue, color intensity, floral, overall fruity, herbaceous, acidity, and astringency) were selected which could adequately characterize the sensory difference among targeted wines. In the third session, all reference standards representative of attributes were prepared to train the panelists. In the fourth session, all panelists were asked to recognize the attributes and rate the intensities of reference standards. Panel performance was assessed after training, including discrimination ability, repeatability, and reproducibility, according to the method in our previous study [21].
For formal session, panelists were asked to rate the intensity of the attributes of wine samples from 1 (weak) to 10 (strong) using a 10-point scale according to a previous study [22].

Statistical Analysis
A one-way analysis of variance (ANOVA) was conducted to determine any significant differences between wines of different regions by using the aov and Duncan.test functions in the agricolae package of the R software environment (version 3.0.3, http://www.rproject.org/, accessed on 10 October 2021). A significance level α = 0.05 was used for the statistical evaluation. Orthogonal partial least squares discriminant analysis (OPLS-DA) was performed using soft independent modeling of class analogy (SIMCA, version 14.1 from Umetrics), to classify wine samples into different classes based on wine-producing region and identify the chemical variables most responsible for the differentiations. All the chemical variables were normalized before the multivariate statistical analysis. Boxplots for the environmental factors were prepared with the ggplot2 package in the R environment.

Meteorological Data
The climatic conditions of wine regions in China vary considerably, leading to significant variation in the quality and styles of wine. According to the Köppen-Geiger climate classification system, which is a widely used vegetation-based empirical climate classification system [23], the JDP region has a monsoon-influenced humid subtropical climate, and BHB has a monsoon-influenced hot-summer humid continental climate (Table S2). Meanwhile, the HZB, LP, and XJ regions have either a cold semi-arid climate or cold desert climate. Further investigation of the climatic conditions of wine-producing regions was achieved by analyzing the meteorological data. The climate results showed that the XJ region is characterized by the longest sunshine duration (on average, 9.7 h), highest PAR (on average, 32.98 mol/(m 2 .d) in Yanqi County), and diurnal temperature range (on average, 13.07 • C in Manasi County and 15.10 • C in Yanqi County), with the lowest rainfall (on average, 115.88 mm in Manasi County and 61.83 mm in Yanqi County) ( Figure 1). In contrast, the JDP and BHB regions displayed significantly shorter sunshine durations, lower PARs, and diurnal temperature ranges, with significantly higher rainfall. These regions have either a warm or hot climate. The Linfen and Qingdao subregions, of LP and JDP, respectively, have a relatively higher average daily temperature than the other regions due to their relatively lower latitudes (Figures 1 and S2). In addition, a comprehensive investigation of the climatic characteristics of wine-producing regions in China was carried out according to a newly established zoning index system named the FRD-DI-T system [24]. It found that the HZB, LP, and XJ regions (DI > 1.6) showed higher levels of the dryness index (DI) when compared with those of BHB and JDP (DI ≤ 1.6), even though HZB is near to the BHB and JDP regions.

Basic Wine Compositions
The ANOVA results of the basic wine compositions, including ethanol, reducing sugar, total acidity, pH, and glycerol, are presented in Table 1. The Marselan wines from XJ, LP, and HZB demonstrated a higher mean concentration of ethanol than the wines from JDP and BHB. The high diurnal temperature ranges in XJ, LP, and HZB, as well as the application of the late harvest technique in HZB, could explain these phenomena. Reducing sugars showed higher concentrations in wines from XJ, followed by wines from LP. Similar to alcohol results, glycerol in wines from XJ, LP, and HZB displayed relatively a higher mean concentration compared with wines from JDP and BHB. No differences existed between the pH levels. Excluding JDP, similar pH results were found in the total acidity of the other four regions. This could be because the vinification technique of adjusting acidity with tartaric acid is widely applied in the wine production in western regions (e.g., XJ, LP, and Ningxia).

Aroma Compounds
The OPLS-DA was applied to identify the most characteristic marker compounds for wine region discrimination. As shown in Figure 2A,B, a clear separation was obtained by a reliable OPLS-DA model (R 2 X = 0.708, R 2 Y = 0.799, and Q 2 = 0.545), except for the similarity between the XJ and LP Marselan wines. Interestingly, the OPLS-DA models showed a clear separation between two classes (Class 1 includes JDP and BHB; Class 2 includes HZB, LP, and XJ) based on the first component. It is well known that the discrimination of similar-style wines made from the same variety in different regions is typically ascribed to differences in the environmental conditions. The volatile aroma compounds that showed differences (p < 0.05) across the regions are presented in Table 2. As expected, the majority of the potential maker compounds for discrimination displayed differences in concentration. In order to interpret the volatile aroma compound differences between varieties, the potential marker compounds with differences are summarized and discussed according to their origins (e.g., grapes and The volatile aroma compounds that showed differences (p < 0.05) across the regions are presented in Table 2. As expected, the majority of the potential maker compounds for discrimination displayed differences in concentration. In order to interpret the volatile aroma compound differences between varieties, the potential marker compounds with differences are summarized and discussed according to their origins (e.g., grapes and fermentation) and chemical categories (e.g., terpenes, norisoprenoids, ethyl esters, and acetates).
For terpenes, there were higher average concentrations of β-citronellol and geraniol in JDP and BHB regions compared to those in the HZB, LP, and XJ regions, except for geraniol in HZB, which showed no difference between JDP and BHB ( Table 2). These grape-derived aroma compounds could have originated directly from the grape berry or from the release of their glycoside precursors during the fermentation and storage of wine [39]. Furthermore, factors such as cultivar, climatic condition, grape maturity, and winemaking procedure could play a vital role in modulating the volatile terpene concentrations in grapes and wines, finally determining the expression of aroma characteristics [40]. Research focusing on the comparison of the aroma compounds in Cabernet Sauvignon and Merlot wines from four regions in China showed that both varieties from HBCL (a subregion of BHB, plotted in the map shown in Figure S2) contained higher concentrations of citronellol than those in HBSC (a synonym of the HZB region, Figure S2) [41]. Additionally, according to the results of our previous studies aimed at investigating the differences in the volatile compounds in Muscat Blanc à Petits Grains grapes between two regions, we found that mature grape berries from CL contained higher concentrations of free and glycosidically bound terpenes than berries from the GT region, which is considered to have a cool continental climate due to its high altitude and is also characterized by high sunshine duration, PAR, and diurnal temperature ranges, with low rainfall [42]. Previous studies concluded that the grape berries exposed to sunlight had increased content of both free and glycosidically bound monoterpenes (e.g., linalool, β-citronellol, and geraniol) in Gewürztraminer [43], Riesling [39], and Pinot Noir [44]. However, in this study, conflicting results were found in the comparison of Marselan wines among regions with different climates. The HZB, LP, and XJ regions have relatively high sunshine duration and PAR compared to the JDP and BHB regions. As we know, the environmental factors that affect the accumulation of the aroma compounds in grapes or wines were obtained from many viticultural trials, such as leaf removal, shading, and canopy management, which influence the microclimate conditions of the fruit zone. Until now, the difference of the terpenoid metabolism in grapes under different climatic conditions, especially monsoon-influenced climates in China, has remained largely unstudied.

Norisoprenoids
Regarding norisoprenoids, (E)-β-damasenone had the highest OAVs compared to other volatiles, ranging from 189.4 (in XJ) to 306.2 (in JDP), which could be attributed to its relatively low odor threshold (0.05 µg/L in 10% w/w water/ethanol solution [22]) ( Table 2). This result was similar to that of a previous study [10], which indicated this compound could be considered the key odor-active compound in Marselan wines. (E)-β-damasenone in JDP and BHB regions presented higher concentrations compared to those in the HZB, LP, and XJ regions ( Table 2). Marais et al. (1992) reported that Riesling wines from warm or hot climates (e.g., South Africa) displayed higher concentrations of norisoprenoids (e.g., TDN and vitispirane) compared to those from cool climates (e.g., Germany and Northern Italy) [45]. Previous research in our lab was carried out on Cabernet Sauvignon grapes to investigate the differentiation in the biosynthesis of norisoprenoids between the CL and GT regions by application of the metabolomic and transcriptomics techniques [46]. However, higher concentrations of (E)-β-damascenone, as well as total norisoprenoids in mature grapes, were shown in GT (cool climate) than in CL (warm climate) [46]. According to the meteorological data, the XJ region shows similarities to the GT region in terms of the high sunshine duration, PAR, and diurnal temperature range; however, the XJ region has a significantly warmer climate during the growing season than the GT region. It seemed that a complex mechanism was involved in the effect of the environmental factors on the metabolism of norisoprenoids in grapes. Additionally, numerous researchers have investigated the influences of microclimatic conditions on the concentrations of volatile aroma compounds and their precursors. In an investigation on the effects of vine microclimate on norisoprenoid concentrations, Cabernet Sauvignon grapes without leaf removal, as well as their corresponding wines, were found to have the highest concentration of (E)-β-damascenone [47]. However, the concentrations of (E)-β-damascenone in grapes were linearly and positively correlated with increasing sunlight exposure when leaves were removed. Similar results were obtained by a leaf removal treatment performed on Pinot Noir [44]. Meanwhile, Feng et al. (2015) only found higher concentrations of (E)-βdamascenone in grapes with 100% leaf removal. Water deficiency also favors the production of bound (E)-β-damascenone in grape berries, such as in Cabernet Sauvignon [48] and Merlot [49]. It seemed that the findings from the viticultural trials related to the modulation of vineyard microclimate could not be directly applied to interpret the differences in the aroma characteristics of wines from different regions with various climatic conditions. Therefore, intensive studies should be carried out to investigate the metabolic mechanism of aroma compounds in grapes in response to either a different climate style or changes to vine microclimate by viticultural practice.

Fermentation Aroma Compounds
Acetic acid, one of the key fermentation aroma compounds, showed higher concentrations in XJ wines compared to those from other regions, excluding wines from LP. As we know, the production of ethanol from reduced sugars during alcoholic fermentation occurs along with the production of acetic acid, as well as glycerol. The higher initial content of sugar must increase the production of acetic acid as well as glycerol by yeast metabolism [50,51], even if the initial sugar ranges from 224-268 g/L, which is much lower than that in the musts for ice wine or botrytized wine production. Yeast increases its intracellular ambulation of glycerol to counterbalance the osmotic pressure induced by high sugar content. In this study, the strongly positive correlation between glycerol and ethanol (r = 0.79, n = 39, p < 0.01), as well as between acetic acid and ethanol (r = 0.69, p < 0.01), reveals that the differences in the musts from different regions may be involved in regulating the yeast metabolism during alcoholic fermentation ( Figure 3). Furthermore, higher initial must sugar concentrations (e.g., ice wine and botrytized wine musts) contribute to the hyperosmotic stress environment, which could upregulate the expression of ALD3 and GPD1 in yeast during alcoholic fermentation and lead to relatively high concentrations of acetic acid in final wines [52]. Combining the results mentioned above, relatively lower concentrations of several fermentation aroma compounds (e.g., higher alcohol acetates, ethyl esters of fatty acids, higher alcohols, fatty acids, and other esters) in wines with higher alcohol concentrations reveal that the production of volatile aroma compounds in yeast could be regulated by the initial sugar concentrations in musts. Therefore, the difference of the initial sugars in musts between cool and warm regions may play a vital role in influencing the sensory differences in final wines by modulating the yeast metabolism. ALD3 and GPD1 in yeast during alcoholic fermentation and lead to relatively high concentrations of acetic acid in final wines [52]. Combining the results mentioned above, relatively lower concentrations of several fermentation aroma compounds (e.g., higher alcohol acetates, ethyl esters of fatty acids, higher alcohols, fatty acids, and other esters) in wines with higher alcohol concentrations reveal that the production of volatile aroma compounds in yeast could be regulated by the initial sugar concentrations in musts. Therefore, the difference of the initial sugars in musts between cool and warm regions may play a vital role in influencing the sensory differences in final wines by modulating the yeast metabolism. For other fermentation aroma compounds, the Marselan wines from JDP and BHB contained relatively higher concentrations of isoamyl acetate, hexyl acetate, phenethyl acetate, octanoic acid, decanoic acid, ethyl decanoate, ethyl dodecanoate, isoamyl octanoate, 1-decanol, and 1-dodecanol (Table 2). Among these differentiated aroma compounds, the OAVs of isoamyl acetate (ranging from 21.1 in HZB to 55.0 in JDP), octanoic acid (ranging from 3.6 in HZB to 4.1 in JDP), decanoic acid (ranging from 2.3 in LP to 2.5 in JDP), and ethyl decanoate (ranging from 3.5 in LP to 7.4 in JDP) were greater than 1. These results indicated that Marselan wines from JDP and BHB regions may be characterized by a higher-intensity fresh fruity aroma. These volatile compounds are formed by the yeast metabolism from various precursors (mainly dominated by sugars, amino ac- For other fermentation aroma compounds, the Marselan wines from JDP and BHB contained relatively higher concentrations of isoamyl acetate, hexyl acetate, phenethyl acetate, octanoic acid, decanoic acid, ethyl decanoate, ethyl dodecanoate, isoamyl octanoate, 1-decanol, and 1-dodecanol (Table 2). Among these differentiated aroma compounds, the OAVs of isoamyl acetate (ranging from 21.1 in HZB to 55.0 in JDP), octanoic acid (ranging from 3.6 in HZB to 4.1 in JDP), decanoic acid (ranging from 2.3 in LP to 2.5 in JDP), and ethyl decanoate (ranging from 3.5 in LP to 7.4 in JDP) were greater than 1. These results indicated that Marselan wines from JDP and BHB regions may be characterized by a higher-intensity fresh fruity aroma. These volatile compounds are formed by the yeast metabolism from various precursors (mainly dominated by sugars, amino acids, and fatty acids). Similar to acetic acid and glycerol, the differences of the fermentation aroma compounds in wines among regions could also be affected by the differences of initial sugar in musts through yeast alcoholic fermentation. Furthermore, many studies comparing the naturally sweet wines made from dehydrated grapes with the corresponding base wines made from fresh grapes have concluded that the former contained higher concentrations of fatty acids, ethyl esters of fatty acids, and higher alcohol acetates [21,53,54]. Certainly, the formation mechanism of the volatile aroma compounds during alcoholic fermentation is complex. The initial concentrations of amino acids, fatty acids, and phenolic compounds in must could also affect the yeast metabolism during alcoholic fermentation [55]. However, how the matrix of grape must affects the volatile compounds of finished wine needs to be further studied.

Chromatic Parameters and Phenolic Compounds
According to colorimetric measurements, Marselan wines from HZB, LP, and XJ regions displayed lower L values and higher values of a and C ab compared with the wines from JDP and BHB (Table 3). These differences may contribute to the relatively deeper color intensity and higher red-hue intensity in wines from HZB, LP, and XJ regions. No significant differences in b and H ab values existed among different wines. The OPLS-DA model (R 2 X = 0.591, R 2 Y = 0.476, Q 2 = 0.183) could not discriminate well between the Marselan wines from different regions based on the quantitative data of the phenolic compounds ( Figure 2C,D). However, the Marselan wines from XJ were separated from the wines from other regions due to their higher concentrations of several flavonols (e.g., myricetin-3-O-glucoside, quercetin-3-O-galactoside, isorhamnetin-3-O-glucoside, syringetin-3-O-glucoside) (p < 0.05) ( Table 4). This suggests that aromatic and phenolic compounds in wines from different regions were influenced by climatic conditions in a different manner. The average total flavonol concentrations in XJ Marselan wines were approximately 2.98-4.31-fold higher than the wines from other regions. This is probably because the XJ region has the longest sunshine duration and highest PAR (Figure 1), which induces a more active synthesis of flavonols to resist the UV radiation [56]. Our previous study also found that more flavonols were accumulated in grapes from the GT region before the berry harvest than in those from the CL region [57]. Lower concentrations of hydroxybenzoic acids (e.g., gallic acid) in JDP and BHB wines were observed (Table 4). Our previous studies reported that Cabernet Sauvignon wines from western regions can be distinguished from those from eastern regions in China based on the quantitative results of the phenolic compounds [19,58]. Li et al. (2017) reported that the five flavanols (including catechin, epicatechin, procyanidin B1, procyanidin B2, and procyanidin C1) and gallic acid were identified as the key phenolic compounds for regional differentiation and contained higher concentrations in western regions [19]. In this study, these five flavanols had VIP scores ≥ 1 (data not shown) and displayed the highest concentration in the Marselan wines from the XJ region than from the other regions, except for epicatechin (Table 4). However, little is known that allows us to interpret the differences in aromatic and phenolic profiles of wines from different regions with various climatic conditions.

Sensory Profiles
In this study, sensory analysis of young Marselan wines from five wine-producing regions was studied. Panelists evaluated the significant differences among wines in the attributes of hue, color intensity, floral, fruity, herbaceous, acidity, and astringency, and the mean intensities of those attributes and standard deviations are summarized in Table 5. The results showed that young Marselan wines from LP and XJ regions had higher color intensity than wines from other three regions, with a decreasing order of LP > XJ > BHB > HZB > JP. This result was similar to that of the L value described above, except that the wines from HZB had relatively low L values, but with light color intensity. Floral and fruity attributes were scored lowest in wines from HZB region compared to wines from other regions in this study, which was different from the results of terpenes and esters described above. Wines from XJ region had lowest acidity compared to wines from other regions, although XJ Marselan wines had relatively high levels of total acidity. This could be explained by the higher concentrations of reducing sugar and glycerol in XJ Marselan wines, which can mask the expression of organic acids. JP wines had the highest acidity, but showed no significant difference from other three regions. There were no significant differences in hue, herbaceous, and astringency attributes among five region Marselan wines in this study, which could be due to the young Marselan wines used in this study.

Conclusions
Using OPLS-DA on the volatile aroma compound data, successful differentiation of young Marselan wines according to geographical origin was achieved. Higher concentrations of β-citronellol, geraniol, (E)-β-damascenone, and several fermentation aroma compounds were observed in young Marselan wines from JDP and BHB regions with similar climatic conditions in grape growing season (short sunshine duration, low PARs and diurnal temperature ranges, and high rainfall). Different initial sugar concentrations in musts from different regions could influence the aroma compounds in final wines by involvement in the modulation of yeast metabolism during alcoholic fermentation. These results provide a clue for further research that aims to interpret the differences in flavor compounds between different regions. However, it is difficult to differentiate the young Marselan wines between different regions based on the phenolic compound data. This research was limited by the number of Marselan wines due to the current development situation of this variety in China. For further studies, a larger dataset of quantitative results of Marselan wines obtained by high-resolution MS would be helpful to identify the aroma and flavor characteristics from different geographical origins.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/10 .3390/foods11060787/s1, Figure S1: Photographs of Marselan grape at veraison stage (A) and mature stage (B), and young wine in a glass (C), Figure S2: Regional distribution of wine production in China in this study, Table S1: Number of Marselan wine samples with various vintages from different wine regions in China, Table S2: Climatic type of wine-producing regions in China in this study according to the Köppen-Geiger climate classification system.  Data Availability Statement: All data generated or analyzed during this study are included in this article.