Next Article in Journal
Defining Acute Coronary Syndrome through Metabolomics
Previous Article in Journal
Fungal Depsides—Naturally Inspiring Molecules: Biosynthesis, Structural Characterization, and Biological Activities
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Oak Species Quercus robur L. and Quercus petraea Liebl. Identification Based on UHPLC-HRMS/MS Molecular Networks

1
Institut de Chimie Organique et Analytique, Université d’Orléans-CNRS, UMR 7311 BP 6759, CEDEX 2, 45067 Orléans, France
2
Centre de Biophysique Moléculaire, CNRS-Université d’Orléans, UPR 4301, CEDEX 2, 45071 Orléans, France
*
Author to whom correspondence should be addressed.
Metabolites 2021, 11(10), 684; https://doi.org/10.3390/metabo11100684
Submission received: 10 September 2021 / Revised: 30 September 2021 / Accepted: 4 October 2021 / Published: 6 October 2021
(This article belongs to the Section Plant Metabolism)

Abstract

:
Two species of oak are dominant in French forests: pedunculate oak (Quercus robur L.) and sessile oak (Quercus petraea Liebl.). Their differentiation is not straightforward but is essential to better understand their respective molecular content in order to better valorize them. Thus, to improve oak species identification, an untargeted UHPLC-HRMS/MS method associated with a two-step data treatment was developed to analyze a wide range of specialized metabolites enabling the comparison of both species of oak extracts. Pooled extracts from sessile and pedunculate oaks, composed of extracts from several trees of pure species from various origins, were compared using first the Venn diagram, as a quick way to get an initial idea of how close the extracts are, and then using a molecular network to visualize, on the one hand, the ions shared between the two species and, on the other hand, the compounds specific to one species. The molecular network showed that the two species shared common clusters mainly representative of tannins derivatives and that each species has specific molecules with similar fragmentation patterns, associated in specific clusters. This methodology was then applied to compare these two pooled extracts to unknown individuals in order to determine the species. The Venn diagram allowed for the quick presumption of the species of the individual and then the species could be assigned more precisely with the molecular network, at the level of specific clusters. This method, developed for the first time, has several interests. First, it makes it possible to discriminate the species and to correctly assign the species of unknown samples. Moreover, it gave an overview of the metabolite composition of each sample to better target oak tree utilization and valorization.

Graphical Abstract

1. Introduction

Oak wood is very widespread and used for many applications. Two species are predominant in French forests, sessile oak (Quercus petraea Liebl.) and pedunculate oak (Quercus robur L.). Some morphological differences, in particular acorns or even leaves, allow foresters to differentiate them in the forest [1]. However, these morphological criteria are not always easy to visualize, especially for coopers, who select the logs they purchase on the roadside. Cooperage turned to a macroscopic anatomical criterion of the wood, the width of tree annual growth or “grain”, which is easily identifiable, and could be correlated with the composition of the wood. Nevertheless, wines and spirits aged in oak barrels with the same grains still present a significant sensory heterogeneity.
Oak barrels are not just containers used for transport but have become an essential tool for aging wines and spirits. The oak used for wine aging modifies the sensory characteristics of wine by releasing various wood compounds, such as aromatic compounds [2,3,4,5], tannins [6,7,8], bitter [9,10], and sweet [11,12,13] tasting compounds into the beverage and by the slow oxygenation process occurring during aging [14,15,16]. The differentiation of species was initially based on tannin content or on the analysis of volatile compounds. Statistically, pedunculate oak is richer in tannin [17], while sessile oak is richer in aromatic compounds [18]. These methods have proven to be unreliable, however, due to the variability of oak composition of individual trees belonging to the same species [17,19,20]. This variability highlights the difficulty of differentiating botanical species based solely on a few chemical parameters.
There are now other reliable methods including genetic analysis [21,22,23], which requires the use of databases and fresh tissue or targeted chemical analyses based on a calculated ratio of quercotriterpenosid (QTT)/glycosylated bartogenic acid (GluBA) content that defines a triterpenoid index used to determine the oak species [24,25]. Other non-targeted analysis methods have been developed to identify metabolite families and highlight discriminating compounds of oak species [26,27]. Still others try to determine the probability of a species’ distribution by non-parametric classification approaches [28].
In recent years, molecular networks emerged. This approach consists of organizing and visualizing tandem mass spectrometry data through spectral similarities. Compounds presenting a similar fragmentation pathway are grouped in clusters that highlight the structural relationship between compounds belonging to the same molecular family [29,30]. Crude extracts of various origins containing many metabolites can thus be more quickly characterized by comparing experimental spectral data to each other and to databases. Moreover, molecular networks help to point out compounds of interest in one or several extracts in comparison in order to focus the analysis on these compounds specifically without spending time on known compounds.
The aim of the present study was to develop a new method able to compare sessile and pedunculate oak samples and to determine the species of unknown samples. Firstly, HHPLCMS/MS analyses of oak wood extracts from sessile and pedunculate oaks were compared using the Venn diagram combined with molecular networking to highlight the differences in profiles between the two species. Secondly, these species profiles were then assessed to determine the species of unknown individual samples due to their distribution between extracts of sessile and pedunculate oaks in the Venn diagram and in the network.
This methodology allowed for an easier comparison of samples from different species of wood, highlighting the differences in molecular content and a rapid determination of the species for unknown samples supported by the detection of numerous relevant compounds belonging to specific families. Moreover, it affords a better understanding of the molecular composition of the samples that would improve the valorization of oak wood by foresters who cannot determine the species in forests reliably and better target the relevant applications for this content. In addition, it would allow the determination and optimization of the sensory quality of casks, and thus helps cooperage to make a reasoned choice when selecting the type of wood according to the desired sensory characteristics.

2. Results and Discussion

2.1. Comparison of Sessile and Pedunculate Oak Chemical Composition

Sessile and pedunculate chemical compositions were first compared using a Venn diagram, presented in Figure 1.
For both samples, more than 1700 ions were detected in negative ionization mode. The Venn diagram shows that the majority of the ions (1644) are common to sessile and pedunculate oaks, corresponding to molecules related to the genus Quercus [31]. Some ions are specific to each of the species, namely, 129 for sessile oaks and 130 for pedunculate oaks.
In order to correlate this first comparison of the dataset with the mass spectra and to describe the ions that compose it, a molecular network was created using MS/MS analysis of pedunculate and sessile oak samples. The aim was to determine the differences in molecular content between these two species without focusing on a particular molecular family.
First, MS/MS spectra were collected and processed on MetaboScape. Then, a quantitative molecular network was generated with Global Natural Product Social Molecular Networking (GNPS) and was visualized using Cytoscape. Molecules generating similar MS/MS spectra are clustered due to similarities in their fragmentation patterns. The resulting molecular networks, presented in Figure 2, visualize chemical relationships of the compounds contained in pedunculate (in red) and sessile (in green) extracts and their relative abundance in the two samples by the proportion of red and green in each node.
Standards were not available for each cluster, so the clusters were annotated thanks to online databases (PubChem, Lotus, SciFinder, GNPS), the interpretation of the spectra, and the literature. For each cluster, two attempts of annotations are reported in Table 1. Molecules with the same fragmentation pattern and therefore belonging to the same molecular family are associated in the same cluster. As a result, the nodes close to the annotated molecule can be more easily associated with the same type of compounds.
The molecular network shows 11 major clusters. Those numbered from 1 to 7 are composed of ions present in equal amounts in the two oak species, corresponding to the shared molecules highlighted in the Venn diagram. The average retention times of these compounds vary between 3 and 7 min, corresponding to the elution of the most polar compounds in the extracts. Molecules present in the same cluster show common fragment ions that support their association. For each cluster, representative MS/MS spectra are available in supplementary data, Figure S1.
For Cluster 1, the main similar fragments observed were [C27H19O18] and [C14H5O8] corresponding to castalin or vescalin and to ellagic acid, respectively. This cluster includes compounds from the ellagitannin family.
Cluster 2 is composed of ellagic acid derivatives with [C14H5O8] as the common fragment of ellagic acids.
Cluster 3 contains multi-charged ions that have a larger m/z. Fragments were identified as vescalagin or castalagin [C41H25O26], castalin or vescalin [C27H19O18], and ions corresponding to grandinin or roburin E [C46H33O30]. The molecules in this cluster are assimilated to complex tannins of high molecular weight derived from the ellagitannin family.
Cluster 4 groups molecules derived from gallic acid. Common losses of 170 mass units corresponding to the loss of gallic acid or even losses of 152 mass units corresponding to galloyl units were observed.
Cluster 5 shows common losses of 44 mass units corresponding to a decarboxylation, highlighting carboxylic acid functions. Well-known fragments of vescalagin or castalagin and castalin or vescalin were also observed for these compounds of Cluster 5. Losses of 16 mass units (oxygen atom) were observed for deoxy-carboxyvescalagin fragments compared to vescavaloninic acid ones. Cluster 5 groups vescalagin or castalagin acid derivatives.
Cluster 6 is made of molecules built with a unitary brick of ellagic acid [C14H5O8] and hexahydroxy diphenoyl glucose (HHDP-glucose) [C27H14O9]. This cluster is again a combination of hydrolyzable tannins from the ellagitannin family.
Cluster 7 groups together compounds with two common ions 183.0655 m/z and 313.0574 m/z. Depending on the molecules, these fragments are not the same moiety. For 3,4,5-trimethoxyphenyl-(6′-O-galloyl)-O-β-glucopyranoside, the first example in Table 1, 183.0655 m/z corresponds to a fragment of trimethoxyphenol, while for 3-methoxy-4-hydroxyphenol 1-O-β-d-(6′-O-galloyl)glucopyranosid, the second example, it is due to a methylgallate moiety.
The common fragment 313.0574 m/z corresponds to the glucogallin moiety (sugar linked to a gallic acid) for the first compound, while it corresponds to the sugar moiety linked to the methoxyhydroxyphenol [Glucogallin−H2O−CH2] for the second.
These fragments demonstrate a grouping of compounds formed from a glucogallin and a phenol derivative from the family of phenol glucosides.
All these clusters are common to both species. While pedunculate oaks are often considered to be richer in tannin than sessile oaks, the results obtained show that this cannot be applied to all individuals, since with a pooled sample composed of trees from different origins the mean tannin composition of sessile and pedunculate oaks is quite similar. The high variability of individual tannin composition was already described [17,19,20] and our results confirm that this criterion cannot be used alone to differentiate the two oak species.
Two clusters clearly appear to be specific to each oak species: Clusters 8 and 9 present compounds specific to sessile oaks while Clusters 10 and 11 present compounds specific to pedunculate oaks.
Cluster 8 contains ions specific to sessile oak, but unlike the other clusters, which are almost exclusively made of compounds specific to a species, this cluster is nevertheless shared to some extent. Among the specific molecules, a lactone precursor C22H32O12 well described in sessile oaks was observed. It was also demonstrated in a previous study that this compound is a molecular marker of sessile oaks [26,38].
Cluster 9, also composed of molecules belonging to the triterpenoid family, is specific of sessile oaks with glycosylated terpenoids such as quercotriterpenoids [11,13]. Currently, around ten quercotriterpenoids have been identified. This cluster suggests the presence of other compounds of lower intensity, which may have related structures.
Clusters 10 and 11. The ions contained in these clusters are triterpenoid compounds: Cluster 10 contains bartogenic acid and oleanane derivatives and Cluster 11 contains triterpenes of larger size in the form of a dimer with the unitary brick of roburgenic or bartogenic acid.
Molecular networks provide an overview of hydro-ethanolic extracts of sessile and pedunculate oaks by rapidly highlighting common compounds and species-specific ones. These results are consistent with previous studies that reported ubiquitous tannin derivatives in both species; a higher content in GluBA and terpene derivatives in pedunculate oaks, whereas sessile oaks contain more QTT and lactone precursors [24,25,26].

2.2. Species Assignment of an Unknown Sample

In order to assess whether the molecular network could be used to determine the species of an unknown oak sample and thus be a tool to better identify a species on the basis of its molecular composition, extracts from individual oak trees were introduced into the molecular network and compared to the pooled samples of both species. The methodology was implemented to determine the species of 10 samples collected in various forests. The results presented illustrate the distribution of one pedunculate and one sessile oak individual, in the Venn diagram and in the network. Similar results were obtained for the other samples.

2.2.1. Addition of Pedunculate Oak Individuals to the Network

First, a Venn diagram was built in order to compare the MS/MS data of the pedunculate oak individual with MS/MS data of pooled sessile and pedunculate oak samples (Figure 3).
Again, most of the ions are common for the three samples (1486). The number of ions shared only with the pedunculate oak pool (109) is greater than the number of ions shared with the sessile oak pool (68), in favor of a suitable assignment. Some ions (67) are specific to the individual, which can be explained by the high intra-individual variability even within the same species.
A quantitative molecular network was again generated with addition of the individual to see how it was distributed in the network (Figure 4).
On the one hand, the individual (in blue) is fairly distributed within the sessile pool and the pedunculate pool for the common Clusters 1 to 7 described in Table 1. On the other hand, the individual shares very few ions with triterpenoids derived from quercotriterpenosids, which are specific of sessile oaks in Cluster 9. It is also not distributed in Cluster 8 regarding the lactone precursors (m/z 487.1821 and 639.1908), which are also markers of the sessile oak pool. However, it is well distributed with the other ions of Cluster 8, which are shared between the two pools of species, highlighting the lack of selectivity of these compounds in the extracts compared to the lactone precursors. As for the two pedunculate-specific Clusters 10 and 11, the distribution is clearly visible across all nodes. Different relative abundances of ions are observed. Some molecules are fairly well shared between the pedunculate pool sample and the individual sample, whereas others are more abundant in either sample due to the specific chemical composition of the individual. Considering both the sharing of the pedunculate-specific molecules and the non-sharing sessile ones, this individual was assigned to the pedunculate species, which was confirmed by the genetic analyses of the individuals.

2.2.2. Addition of Sessile Oak Individuals to the Network

In order to illustrate the method with a sessile oak individual, the same approach was implemented. The Venn diagram is shown in Figure 5.
It again highlights a majority of common ions (1288). The number of ions shared only with the sessile oak pool (91) is greater than the number of ions shared with the pedunculate oak pool (45).
The distribution of the individual in the network is shown in Figure 6.
The individual (in blue) is well distributed within the sessile pool and the pedunculate pool in the common clusters (1 to 7) described in Table 1. The individual is not present in Clusters 10 and 11 that are specific to pedunculate oak (only 2 nodes out of the 32 in the cluster). On the contrary, it is distributed in Cluster 9 composed of the triterpenoids derived from quercotriterpenosid of the sessile oak pool and also in the few specific compounds (m/z 487.1821 and 639.1908) in Cluster 8. Thus, this sample was assigned to the sessile species which was confirmed by the genetic analysis.
These two individuals illustrated the possibility to well determine both species of sessile and pedunculated oak. First, the Venn diagram points out twice more common ions between the individual sample and pooled extract of the suited species. Then, the network consolidates species assignation confirming that a large number of species-specific molecules are well detected in the sample. The molecular network makes it possible to validate this distribution on a large number of markers and not only on a few major molecules in contrast to methods using targeted compounds.

3. Materials and Methods

3.1. Chemicals

Ethanol used for oak wood extraction was HPLC grade and was purchased from VWR (Fontenay-sous-Bois, France). Ultra-pure water was produced with a Purelab Flex system from Veolia (Wissous, France). Acetonitrile, water, and formic acid used for UHPLC-HRMS analysis were of Optima LC-MS grade from Fisher Scientific (Illkirch-Graffenstaden, France).

3.2. Oak Wood Extract Preparation

3.2.1. Oak Wood Sampling

In order to constitute a representative oak extract which takes into account the intra-species variability which may be due in particular to the geographical location of the oaks, 20 samples were taken in 7 forests in the Centre-Val de Loire region in France. Thirty-two logs were also sampled at the cooperage (oak stave mill) on available wood in the timber yard with as many different woods as possible (growth ring widths/diameters). A wooden disk was cut from each freshly felled tree and the heartwood was collected.

3.2.2. Genetic Analyses for Species Assignation

To distinguish oak species, genetic analyses were conducted by the CGAF-ONF laboratory (UMR BioForA, Orléans, France). For this study, total genomic DNA was extracted on the same trees on the same forest stands and cambial tissues from the same cooperage samples as previously, using the NucleoSpin® 96 Plant II Kit (Macherey-Nagel, Hoerdt, France). Genotype data were obtained on the capillary of an ABI 3500 automatic sequencer (Applied Biosystems, Foster City, CA, USA) using 18 nuclear microsatellite markers developed by Guichoux et al. and analyzed with GeneMapper™ software v4.1 (Applied Biosystems) [45].
The determination of the species was carried out using an assignment method from Pritchard et al. using structured software. The method consists in comparing the multi-locus genotypes with the reference data for each species. If the percentage of assignment is above 80%, the sample is assigned to the species. In total, the sample set was composed of 26 pedunculate oaks and 26 sessile oaks: 21 were used to make pooled extracts and 5 were used as individuals (cf. supplementary data).

3.2.3. Oak Wood Extraction

Two pooled extracts were prepared: one pooled extract of pedunculate oak containing the 21 pedunculate samples and another of sessile oaks was prepared with the same number of sessile samples.
The powdered mixture (2 g) of 21 oak wood samples was extracted using ultrasound-assisted extraction with 30 mL of water/ethanol mixture (85:15 v/v) during 1 h. The extract was then centrifuged at 5000× g during 5 min. The supernatant was recovered and analyzed.
Per species, 5 samples were used to demonstrate the distribution of unknowns introduced into the network. The 5 individuals of sessile oaks came from Aboncourt, Russy, Orléans, and Abbayes forests, while the 5 individuals of pedunculate oaks came from Sully la Chapelle, Chateauroux, Lisledon, and Abbayes forests.

3.3. UHPLC-HRMS/MS Analysis

Chromatographic analyses were performed using an Ultimate 3000 RSLC system equipped with an autosampler, a binary pump, a thermostated column compartment, and a DAD detector (Dionex, Germering, Germany). The column was a Luna® Omega C18 (150 × 2.1 mm; 1.6 µm) (Phenomenex, Le Pecq, France). The mobile phase was composed of water (A) and acetonitrile (B), both acidified with 0.1% of formic acid. Elution was performed at a flow rate of 500 μL min−1 and with the following binary gradient program: starting with 3% of solvent B during 0.2 min, 3–45% B from 0.2 to 12 min, 45–90% B from 12 to 14 min, 90–3% B from 15 to 15.5 min. Then the column was re-equilibrated with 3% of solvent B during 3 min. The column temperature was set at 40 °C and 5 µL of the oak extracts were injected.
The MS and MS/MS experiments were carried out on a maXis UHR-Q-TOF mass spectrometer (Bruker, Bremen, Germany) with an electrospray ion source (ESI), working in negative ionization mode that allows for better detection and isolation of the molecular ion compared to the positive ionization mode. The pressure of nebulizing gas was set to 2 bar, and the flow rate and temperature of dry gas were set at 9.0 L min−1 and 200 °C. The capillary voltage was set at 4 kV. Mass spectra were summed during 400 ms in the m/z range 50–2250. All the MS data were processed using DataAnalysis software version 4.4 (Bruker). Molecular formulae were generated using the SmartFormula algorithm with an elemental composition of C, H, O to an infinite number and N ≤ 4 with a mass accuracy ≤3 ppm and were submitted to the SciFinder, PubChem, Lotus, and GNPS databases in order to propose compound structures.
MS/MS experiments were conducted using Data Dependent Acquisition mode with 3 precursor ions in the m/z range 150–1600; ions were excluded after 6 s. The analysis included 2 segments: in the first one, from 0 to 8 min, collision energies (CE) of 40 eV for monocharged ions and 30 eV for doubly charged ions were applied. In the second one, from 8 min to the end of the analysis, collision energies of 75 eV for monocharged ions and 40 eV for doubly charged ions were applied. CE values were adapted to the kinds of compounds: before 8 min, the molecules belong to the tannin family, which is mostly composed of C-O bonds, weaker than the C-C bonds of the terpenes, which elutes after 8 min. CE values for doubly charged ions are lower than for single-charged ones due to the proximity of both charges. MS/MS spectra were summed during 400 ms, so the total cycle time for MS and MS/MS was 1.6 s.

3.4. Venn Diagrams

A Venn diagram shows the logical relation between sets; it illustrates sample relationships by comparing common ions between them. InteractiVenn software (http://www.interactivenn.net/index2.html, accessed for the last figures on 30 September 2021) was used to build the Venn diagrams [46].
After analyses of the pooled extracts, sets were built using a bucket associating for each detected compound its retention time and the m/z of the most abundant ion. One set is represented by a circle, the overlap between two circles or sets is the ions they have in common, whereas the others are ions that are specific to a set.

3.5. Molecular Network Design

Molecular networks are visual representations of the chemical space present in tandem mass spectrometry (MS/MS) experiments by comparing mass spectra in pairs to map an extract. Therefore, considering that ions with closed structures will give similar neutral losses or fragment ions, each cluster associates molecules with a similar fragmentation pathway that possess high structural similarities and are likely to belong to the same chemical family.
The bucket table was built using MetaboScape software version 4.0 (Bruker). The T-ReX 3D algorithm detected precursor ions with an intensity threshold of 10,000 .a.u. and associated retention time, m/z, and area for each analyte.
The mascot generic format (mgf) file and a quantification table were exported to the GNPS platform (http://gnps.ucsd.edu, accessed for the last treatments on 7 and 8 September 2021) using feature-based molecular networking (FBMN) in order to build the molecular network [29,47]. Optimized parameters used for molecular network design were: mass tolerance 0.02 Da for parent and fragment mass, min pair cos 0.75, network TopK 10, maximum connected component size 100, minimum matched fragment ions 6, minimum cluster size 2, yes run MSCluster.
Quantitative molecular networks were visualized using Cytoscape software version 3.8.1 (https://cytoscape.org, accessed for the last figures on 7 and 8 September 2021). Within the network, one node corresponds to one MS/MS spectrum. Nodes are represented by pie charts. The proportion of each color in the nodes correlates with the relative abundance of the ion in each sample. The more abundant an ion is in a sample (species), the more dominant the color of this sample is in the pie chart.

4. Conclusions

The developed method based on an untargeted UHPLC-MS/MS analysis of several families of specialized metabolites of oak species (phenolic and terpenic compounds) was able to discriminate sessile and pedunculate chemical compositions. Both species showed specific molecules with similar fragmentation patterns associated in specific clusters in the molecular network. The comparison of individual unknown samples with pooled samples composed of extracts from several trees of pure species from various origins enabled the sample species to be assigned using the proposed two-step method: firstly, with the Venn diagram as a quick way to obtain an initial idea of how close an individual is to a species by comparing the number of similar ions detected in the sample vs. the pooled extracts that allow for the presumption of the sample species. Then with the molecular network, it is possible to visualize more precisely, at the level of specific clusters, the ions that are shared or not shared with one of the species. Due to the intra-species variability observed on individual samples, the confidence level of the species identification based on the detection of numerous compounds with related structures belonging to the same molecular family is improved, compared to the ones based on only few markers. Compared to genetic analysis, this method also presents advantages such as easier sample collection and conservation, higher convenience, and lower price. This method could now be applied on any new oak samples (one individual or a set of samples) to help in species identification.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/metabo11100684/s1, Table S1: Composition of sessile pooled extract, Table S2: Composition of sessile pooled extract, Table S3: Listing of individual samples, Figure S1: Representative MS/MS spectra of each cluster.

Author Contributions

Conceptualization, G.B., L.F., and E.D.; methodology, G.B., L.F., and E.D.; software, G.B., C.C., and L.F.; validation, G.B.; formal analysis, G.B.; investigation, G.B.; resources, E.D.; data curation, G.B.; writing—original draft preparation, G.B.; writing—review and editing, G.B., C.C., L.F., and E.D.; visualization, G.B.; supervision, E.D.; project administration, E.D.; funding acquisition, E.D. All authors have read and agreed to the published version of the manuscript.

Funding

The authors gratefully acknowledge major financial support from the region Centre-Val de Loire through the “Chêne et vin” program which made this study possible, and also the projects Valbiocosm (FEDER-FSE 2014-2020-EX003202), CHemBio (FEDER-FSE 2014-2020-EX003677), Techsab (FEDER-FSE 2014-2020-EX011313), the RTR Motivhealth (2019-00131403) and the Labex programs SYNORG (ANR-11-LABX-0029), and IRON (ANR-11-LABX-0018-01) for their financial support of ICOA, UMR 7311, University of Orléans, CNRS.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to [industrial partnerships of the program].

Acknowledgments

We are thankful to Thomas Giordanengo from the Tonnellerie Radoux for the many discussions that helped to carry out this study and for the cooperage sampling of oak woods. We are also grateful for Brigitte Musch and Bénédicte Le Guerroué from the CGAF-ONF laboratory UMR BioForA (Ardon, France) for the genetic assignment of species. We thank Claire Quinones and Xavier Mandret from the ONF Centre-Ouest-Nouvelle-Aquitaine territorial direction (Boigny-sur-Bionne, France) for the forest sampling of oak woods.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kremer, A.; Dupouey, J.L.; Deans, J.D.; Cottrell, J.; Csaikl, U.M.; Finkeldey, R.; Espinel, S.; Jensen, J.S.; Kleinschmit, J.; Van Dam, B.; et al. Leaf morphological differentiation between Quercus robur and Quercus petraea is stable across western European mixed oak stands. Ann. For. Sci. 2002, 59, 777–787. [Google Scholar] [CrossRef]
  2. Boidron, J.-N.; Chatonnet, P.; Pons, M. Influence du bois sur certaines substances odorantes des vins. OENO One 1988, 22, 275–294. [Google Scholar] [CrossRef]
  3. Jarauta, I.; Cacho, J.; Ferreira, V. Concurrent Phenomena Contributing to the Formation of the Aroma of Wine during Aging in Oak Wood: An Analytical Study. J. Agric. Food Chem. 2005, 53, 4166–4177. [Google Scholar] [CrossRef] [PubMed]
  4. Spillman, P.J.; Sefton, M.A.; Gawel, R. The effect of oak wood source, location of seasoning and coopering on the composition of volatile compounds in oak-matured wines. Aust. J. Grape Wine Res. 2008, 10, 216–226. [Google Scholar] [CrossRef]
  5. Zamora, F. Barrel aging; Types of wood. In Red Wine Technology; Academic Press: Cambridge, MA, USA, 2019; pp. 125–147. [Google Scholar] [CrossRef]
  6. Zhang, B.; Cai, J.; Duan, C.-Q.; Reeves, M.J.; He, F. A Review of Polyphenolics in Oak Woods. Int. J. Mol. Sci. 2015, 16, 6978–7014. [Google Scholar] [CrossRef]
  7. Glabasnia, A.; Hofmann, T. Sensory-Directed Identification of Taste-Active Ellagitannins in American (Quercus alba L.) and European Oak Wood (Quercus robur L.) and Quantitative Analysis in Bourbon Whiskey and Oak-Matured Red Wines. J. Agric. Food Chem. 2006, 54, 3380–3390. [Google Scholar] [CrossRef] [PubMed]
  8. Chira, K.; Anguellu, L.; Da Costa, G.; Richard, T.; Pedrot, E.; Jourdes, M.; Teissedre, P.-L. New C-Glycosidic Ellagitannins Formed upon Oak Wood Toasting; Identification and Sensory Evaluation. Foods 2020, 9, 1477. [Google Scholar] [CrossRef]
  9. Winstel, D.; Marchal, A. Lignans in Spirits: Chemical Diversity, Quantification, and Sensory Impact of (±)-Lyoniresinol. Molecules 2018, 24, 117. [Google Scholar] [CrossRef] [Green Version]
  10. Marchal, A.; Cretin, B.N.; Sindt, L.; Waffo-Téguo, P.; Dubourdieu, D. Contribution of oak lignans to wine taste: Chemical identification, sensory characterization and quantification. Tetrahedron 2015, 71, 3148–3156. [Google Scholar] [CrossRef]
  11. Marchal, A.; Génin, E.; Waffo-Téguo, P.; Bibès, A.; Da Costa, G.; Mérillon, J.-M.; Dubourdieu, D. Development of an analytical methodology using Fourier transform mass spectrometry to discover new structural analogs of wine natural sweeteners. Anal. Chim. Acta 2015, 853, 425–434. [Google Scholar] [CrossRef]
  12. Marchal, A.; Waffo-Téguo, P.; Génin, E.; Mérillon, J.-M.; Dubourdieu, D. Identification of New Natural Sweet Compounds in Wine Using Centrifugal Partition Chromatography–Gustatometry and Fourier Transform Mass Spectrometry. Anal. Chem. 2011, 83, 9629–9637. [Google Scholar] [CrossRef] [PubMed]
  13. Gammacurta, M.; Waffo-Teguo, P.; Winstel, D.; Cretin, B.N.; Sindt, L.; Dubourdieu, D.; Marchal, A. Triterpenoids from Quercus petraea: Identification in Wines and Spirits and Sensory Assessment. J. Nat. Prod. 2019, 82, 265–275. [Google Scholar] [CrossRef]
  14. Vivas, N. Modélisation et calcul du bilan des apports d’oxygène au cours de l’élevage des vins rouges. IV-Elevage des vins rouges en conditions d’oxydations ménagées controlées. Progres. Agric. Vitic. 1999, 116, 305–311. [Google Scholar]
  15. Ribéreau-Gayon, P.; Dubourdieu, D.; Glories, Y.; Maujean, A. Handbook of Enology, Volume 2: The Chemistry of Wine Stabilization and Treatments; John Wiley & Sons: Hoboken, NJ, USA, 2021; ISBN 978-1-119-58776-7. [Google Scholar]
  16. Del Alamo-Sanza, M.; Nevares, I. Oak wine barrel as an active vessel: A critical review of past and current knowledge. Crit. Rev. Food Sci. Nutr. 2017, 58, 2711–2726. [Google Scholar] [CrossRef]
  17. Prida, A.; Boulet, J.-C.; Ducousso, A.; Nepveu, G.; Puech, J.-L. Effect of species and ecological conditions on ellagitannin content in oak wood from an even-aged and mixed stand of Quercus robur L. and Quercus petraea Liebl. Ann. For. Sci. 2006, 63, 415–424. [Google Scholar] [CrossRef] [Green Version]
  18. Prida, A.; Ducousso, A.; Petit, R.J.; Nepveu, G.; Puech, J.-L. Variation in wood volatile compounds in a mixed oak stand: Strong species and spatial differentiation in whisky-lactone content. Ann. For. Sci. 2007, 64, 313–320. [Google Scholar] [CrossRef]
  19. Doussot, F.; De Jéso, B.; Quideau, S.; Pardon, P. Extractives Content in Cooperage Oak Wood during Natural Seasoning and Toasting; Influence of Tree Species, Geographic Location, and Single-Tree Effects. J. Agric. Food Chem. 2002, 50, 5955–5961. [Google Scholar] [CrossRef]
  20. Mosedale, J.R.; Savill, P.S. Variation of heartwood phenolics and oak lactones between the species and phenological types of Quercus petraea and Q. robur. Forestry 1996, 69, 47–56. [Google Scholar] [CrossRef]
  21. Kremer, A.; Petit, R.-J.; Ducousso, A. Biologie Évolutive et Diversité Génétique des Chênes Sessile et Pédonculé; INRAE: Cestas, France, 2002; p. 111. [Google Scholar] [CrossRef] [Green Version]
  22. Muir, G.; Fleming, C.C.; Schltterer, C. Species status of hybridizing oaks. Nature 2000, 405, 1016. [Google Scholar] [CrossRef] [PubMed]
  23. Guichoux, E.; Petit, R. Méthode de traçabilité géographique des bois de chêne; Déclaration d’invention (n DIRV-13-00566); INPI: Bordeaux, France, 2014. [Google Scholar]
  24. Marchal, A.; Prida, A.; Dubourdieu, D. New Approach for Differentiating Sessile and Pedunculate Oak: Development of a LC-HRMS Method to Quantitate Triterpenoids in Wood. J. Agric. Food Chem. 2016, 64, 618–626. [Google Scholar] [CrossRef] [PubMed]
  25. Marchal, A.; Dubourdieu, D. Method for Identifying the Oak Species of an Oak Wood Sample. U.S. Patent 20160011166A1, 14 January 2016. [Google Scholar]
  26. Buche, G.; Colas, C.; Fougère, L.; Giordanengo, T.; Destandau, E. Untargeted UHPLC-Q-TOF-HRMS based determination of discrimating compounds for oak species Quercus robur L. and Quercus petraea Liebl. identification. Phytochem. Anal. 2020. [Google Scholar] [CrossRef]
  27. Gougeon, R.D.; Lucio, M.; De Boel, A.; Frommberger, M.; Hertkorn, N.; Peyron, D.; Chassagne, D.; Feuillat, F.; Cayot, P.; Voilley, A.; et al. Expressing Forest Origins in the Chemical Composition of Cooperage Oak Woods and Corresponding Wines by Using FTICR-MS. Chem. Eur. J. 2008, 15, 600–611. [Google Scholar] [CrossRef]
  28. Vayssières, M.P.; Plant, R.E.; Allen-Diaz, B.H. Classification trees: An alternative non-parametric approach for predicting species distributions. J. Veg. Sci. 2000, 11, 679–694. [Google Scholar] [CrossRef]
  29. Wang, M.; Carver, J.J.; Phelan, V.; Sanchez, L.; Garg, N.; Peng, Y.; Nguyen, D.D.; Watrous, J.; A Kapono, C.; Luzzatto-Knaan, T.; et al. Sharing and community curation of mass spectrometry data with Global Natural Products Social Molecular Networking. Nat. Biotechnol. 2016, 34, 828–837. [Google Scholar] [CrossRef] [Green Version]
  30. Yang, J.Y.; Sanchez, L.M.; Rath, C.M.; Liu, X.; Boudreau, P.D.; Bruns, N.; Glukhov, E.; Wodtke, A.; de Felicio, R.; Fenner, A.; et al. Molecular Networking as a Dereplication Strategy. J. Nat. Prod. 2013, 76, 1686–1699. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. de Simón, M.B.F.; Sanz, M.; Cadahía, E.; Poveda, A.P.; Broto-Cartagena, M. Chemical Characterization of Oak Heartwood from Spanish Forests of Quercus pyrenaica (Wild.). Ellagitannins, Low Molecular Weight Phenolic, and Volatile Compounds. J. Agric. Food Chem. 2006, 54, 8314–8321. [Google Scholar] [CrossRef] [PubMed]
  32. Vivas, N.; Glories, Y.; Bourgeois, G.; Vitry, C. The Heartwood Ellagitannins of Different Oaks (Quercus Sp.) and Chestnut Species (Castanea Sativa Mill.). Quantity Analysis of Red Wines Aging in Barrels. J. Des Sci. Tech. Tonnelerie 1996, 25–75. [Google Scholar]
  33. Venter, P. Analysis of Complex Tannins by Multidimensional Techniques. Ph.D. Thesis, Stellenbosch University, Stellenbosch, South Africa, 2019. [Google Scholar]
  34. Cerulli, A.; Napolitano, A.; Hošek, J.; Masullo, M.; Pizza, C.; Piacente, S. Antioxidant and In Vitro Preliminary Anti-Inflammatory Activity of Castanea sativa (Italian Cultivar “Marrone di Roccadaspide” PGI) Burs, Leaves and Chestnuts Extracts and their Metabolite Profiles by LC-ESI/LTQOrbitrap/MS/MS. Antioxidants 2021, 10, 278. [Google Scholar] [CrossRef]
  35. Mämmelä, P.; Savolainen, H.; Lindroos, L.; Kangas, J.; Vartiainen, T. Analysis of Oak Tannins by Liquid Chromatography-Electrospray Ionisation Mass Spectrometry. Available online: https://reader.elsevier.com/reader/sd/pii/S0021967300006245?token=2868AF0C859030365A67BD9580C8D918DE45468890B2BBF50E64B69B3C16913BF46EA1E94311438DAD78BA73AE69031D&originRegion=eu-west-1&originCreation=20210628133609 (accessed on 28 June 2021).
  36. Cammann, J.; Denzel, K.; Schilling, G.; Gross, G.G. Biosynthesis of gallotannins: β-Glucogallin-dependent formation of 1,2,3,4,6-pentagalloylglucose by enzymatic galloylation of 1,2,3,6-tetragalloylglucose. Arch. Biochem. Biophys. 1989, 273, 58–63. [Google Scholar] [CrossRef]
  37. Saucier, C.; Jourdes, M.; Glories, Y.; Quideau, S. Extraction, Detection, and Quantification of Flavano-Ellagitannins and Ethylvescalagin in a Bordeaux Red Wine Aged in Oak Barrels. J. Agric. Food Chem. 2006, 54, 7349–7354. [Google Scholar] [CrossRef] [PubMed]
  38. Centrifugal Partition Chromatography Applied to the Isolation of Oak Wood Aroma Precursors—Elsevier Enhanced Reader. Available online: https://reader.elsevier.com/reader/sd/pii/S0308814613005153?token=E993B2DAF29A3BBE1B6F2EB9836E50326FC688D5DA80E6036E1EB04F1BC6F141B831A98064361AB4E5DE019E76585EDB&originRegion=eu-west-1&originCreation=20210707140737 (accessed on 7 July 2021).
  39. Sindt, L.; Gammacurta, M.; Waffo-Teguo, P.; Dubourdieu, D.; Marchal, A. Taste-Guided Isolation of Bitter Lignans from Quercus petraea and Their Identification in Wine. J. Nat. Prod. 2016, 79, 2432–2438. [Google Scholar] [CrossRef] [PubMed]
  40. Tanaka, T.; Kouno, I. Whisky Lactone Precursors from the Wood of Platycarya strobilacea. J. Nat. Prod. 1996, 59, 997–999. [Google Scholar] [CrossRef]
  41. Masson, E.; Baumes, R.; Le Guerneve, C.; Puech, J.-L. Identification of a Precursor of β-Methyl-γ-octalactone in the Wood of Sessile Oak (Quercus petraea (Matt.) Liebl.). J. Agric. Food Chem. 2000, 48, 4306–4309. [Google Scholar] [CrossRef] [PubMed]
  42. Hayasaka, Y.; Wilkinson, K.L.; Elsey, G.M.; Raunkjær, M.; Sefton, M.A. Identification of Natural Oak Lactone Precursors in Extracts of American and French Oak Woods by Liquid Chromatography—Tandem Mass Spectrometry. J. Agric. Food Chem. 2007, 55, 9195–9201. [Google Scholar] [CrossRef]
  43. Arramon, G.; Saucier, C.; Colombani, D.; Glories, Y. Identification of triterpene saponins in Quercus robur L. and Q. petraea Liebl. heartwood by LC-ESI/MS and NMR. Phytochem. Anal. 2002, 13, 305–310. [Google Scholar] [CrossRef]
  44. Pérez, A.J.; Pecio, Ł.; Kowalczyk, M.; Kontek, R.; Gajek, G.; Stopinsek, L.; Mirt, I.; Oleszek, W.; Stochmal, A. Triterpenoid Components from Oak Heartwood (Quercus Robur) and Their Potential Health Benefits. J. Agric. Food Chem. 2017, 65, 23. [Google Scholar] [CrossRef]
  45. Guichoux, E.; Lagache, L.; Wagner, S.; Chaumeil, P.; Léger, P.; Lepais, O.; Lepoittevin, C.; Malausa, T.; Revardel, E.; Salin, F.; et al. Current trends in microsatellite genotyping. Mol. Ecol. Resour. 2011, 11, 591–611. [Google Scholar] [CrossRef]
  46. Heberle, H.; Meirelles, G.V.; Da Silva, F.R.; Telles, G.P.; Minghim, R. InteractiVenn: A web-based tool for the analysis of sets through Venn diagrams. BMC Bioinform. 2015, 16, 169. [Google Scholar] [CrossRef]
  47. Feature-Based Molecular Networking in the GNPS Analysis Environment—Nature Methods. Available online: https://www.nature.com/articles/s41592-020-0933-6 (accessed on 5 July 2021).
Figure 1. Venn diagram of sessile pooled extracts (21 samples) and pedunculate pooled extracts (21 samples). Numbers above species represent total precursor ions detected in each pooled extract. The circle overlap represents common ions, and the non-overlapping part represents specific ions.
Figure 1. Venn diagram of sessile pooled extracts (21 samples) and pedunculate pooled extracts (21 samples). Numbers above species represent total precursor ions detected in each pooled extract. The circle overlap represents common ions, and the non-overlapping part represents specific ions.
Metabolites 11 00684 g001
Figure 2. Mapping of pooled extracts (21 samples) of sessile (green) and (21 samples) pedunculate (red) oaks with cos 0.75 and matched fragment 6 parameters. (a) Clusters 1–7 representing common clusters between sessile and pedunculated oaks, (b) Clusters 8 and 9 containing specific ions to sessile oaks, (c) Clusters 10 and 11 containing specific ions to pedunculated oaks.
Figure 2. Mapping of pooled extracts (21 samples) of sessile (green) and (21 samples) pedunculate (red) oaks with cos 0.75 and matched fragment 6 parameters. (a) Clusters 1–7 representing common clusters between sessile and pedunculated oaks, (b) Clusters 8 and 9 containing specific ions to sessile oaks, (c) Clusters 10 and 11 containing specific ions to pedunculated oaks.
Metabolites 11 00684 g002aMetabolites 11 00684 g002b
Figure 3. Venn diagram of an unknown individual sample (1 sample) attributed pedunculated compared to sessile pooled extracts (21 samples) and pedunculate pooled extracts (21 samples). Numbers above species represent total precursor ions detected in each extract. The circle overlap represents common ions, and the non-overlapping part represents specific ions.
Figure 3. Venn diagram of an unknown individual sample (1 sample) attributed pedunculated compared to sessile pooled extracts (21 samples) and pedunculate pooled extracts (21 samples). Numbers above species represent total precursor ions detected in each extract. The circle overlap represents common ions, and the non-overlapping part represents specific ions.
Metabolites 11 00684 g003
Figure 4. Molecular network representing the distribution of an individual (1 sample) of pedunculate oak (blue) compared to the pool (21 samples) of sessile (green) and (21 samples) pedunculate (red) oaks with cos 0.75 and matched fragment 6 parameters. (a) Clusters 1–7 representing common clusters shared between sessile, pedunculated and unknown individual, (b) Clusters 8 and 9 containing specific ions to sessile oaks and (c) Clusters 10 and 11 containing specific ions to pedunculated oaks shared with unknown individual.
Figure 4. Molecular network representing the distribution of an individual (1 sample) of pedunculate oak (blue) compared to the pool (21 samples) of sessile (green) and (21 samples) pedunculate (red) oaks with cos 0.75 and matched fragment 6 parameters. (a) Clusters 1–7 representing common clusters shared between sessile, pedunculated and unknown individual, (b) Clusters 8 and 9 containing specific ions to sessile oaks and (c) Clusters 10 and 11 containing specific ions to pedunculated oaks shared with unknown individual.
Metabolites 11 00684 g004
Figure 5. Venn diagram of an unknown individual sample (1 sample) attributed sessile compared to sessile pooled extracts (21 samples) and pedunculate pooled extracts (21 samples). Numbers above species represent total precursor ions detected in each extract. The circle overlap represents common ions, and the non-overlapping part represents specific ions.
Figure 5. Venn diagram of an unknown individual sample (1 sample) attributed sessile compared to sessile pooled extracts (21 samples) and pedunculate pooled extracts (21 samples). Numbers above species represent total precursor ions detected in each extract. The circle overlap represents common ions, and the non-overlapping part represents specific ions.
Metabolites 11 00684 g005
Figure 6. Molecular network representing the distribution of an individual (1 sample) of sessile oak (blue) compared to the pool (21 samples) of sessile (green) and (21 samples) pedunculate (red) oaks with cos 0.75 and matched fragment 6 parameters. (a) Clusters 1-7 representing common clusters shared between sessile, pedunculated and unknown individual, (b) Clusters 8 and 9 containing specific ions to sessile oaks shared with unknown individual and (c) Clusters 10 and 11 containing specific ions to pedunculated oaks.
Figure 6. Molecular network representing the distribution of an individual (1 sample) of sessile oak (blue) compared to the pool (21 samples) of sessile (green) and (21 samples) pedunculate (red) oaks with cos 0.75 and matched fragment 6 parameters. (a) Clusters 1-7 representing common clusters shared between sessile, pedunculated and unknown individual, (b) Clusters 8 and 9 containing specific ions to sessile oaks shared with unknown individual and (c) Clusters 10 and 11 containing specific ions to pedunculated oaks.
Metabolites 11 00684 g006
Table 1. Description of the global network that compares pedunculate oak and sessile oak.
Table 1. Description of the global network that compares pedunculate oak and sessile oak.
Cluster NumberRT (min)Measured m/z
[M−H] or [M−2H]2−
Formula [M]Error (ppm)MS/MS Fragments
[M−H]
Proposed Annotation for MoleculesMolecular Family of the Cluster
14.10933.06251−C41H26O26−0.1631.0582 C27H19O18
300.9998 C14H5O8
Vescalagin or Castalagin [32]Ellagitannins
3.16924.05712−C82H50O51−1.4631.0.581C27H19O18
300.9993 C14H5O8
Roburin D or A [32]
27.15433.04011−C19H14O12−1.5300.9910 C14H5O8Ellagic acid pentoside [33]Ellagic acid derivatives
7.40447.05561−C20H16O120.8300.9913 C14H5O8Methyl ellagic acid pentoside [34]
33.46999.08112−C111H172O310.11065.1018 C46H33O30
975.0744 C43H27O27
933.0650 C41H25O26
631.0575 C27H19O18
Grandinin or Roburin E derivativeEllagitannin derivatives
3.331227.41063−xx1065.1093 C46H33O30
933.0631 C41H25O26
915.0539 C41H23O25
783.0690 C34H23O22
631.0574 C27H19O18
Glycosylated grandinin
46.26635.08731−C27H24O18−0.5483.0678 C20H19O14
465.0678 C20H17O13
313.0568 C13H13O9
169.0148 C7H5O5
Trigalloyl glucose [35]Gallic acid derivatives
7.16939.11581−C34H28O22−0.3769.0892 C34H25O21
617.0753 C34H17O12
447.0573 C20H15O12
169.0151 C7H5O5
Penta-galloylglucose [36]
53.571101.06671−C48H30O31−0.91057.0797 C47H29O29
933.0605 C41H25O26
631.0574 C27H19O18
425.0141 C20H9O11
Vescavaloninic acid
[37]
Vescalagin or
castalagin acid
derivatives
3.71961.06411−C42H26O27−0.6917.0683 C41H25O25
873.0778 C40H25O23
615.0631 C27H19O17
491.062 C25H15O11
Deoxy-carboxy
vescalagin
[7]
61.89481.06111−C20H18O14 0.1300.9951 C14H5O8
275.0161 C13H7O7
HHDP glucose
[34]
Ellagitannins HHDP glucose
derivatives
3.87783.06621−C34H24O22−0.6481.0640 C20H17O14
300.9994 C14H5O8
Pedunculagin
[36]
77.30497.12891−C22H26O13−0.9313.0574 C13H13O9
183.0655 C9H11O4
169.0136 C7H5O5
3,4,5-trimethoxyphenyl-(6′-O-galloyl)-O-β-glucopyranoside [38] Phenol
glucosides
5.58453.10271−C20H22O120.6313.0569 C13H13O9
327.0706 C14H15O9
183.0293 C8H7O5
169.0152 C7H5O5
3-methoxy-4-hydroxyphenol 1-O-β-d-(6′-O-galloyl)glucopyranoside
[39]
89.01487.18211−C22H32O120.1211.0266 C9H7O6
168.0074 C7H4O5
124.0233 C6H4O3
3-Methyl-4-[[6-O-(3,4,5-trihydroxybenzoyl)-β-d-glucopyranosyl]oxy]
octanoic acid
[40,41,42]
Lactone
precursors
10.40639.19081−C29H36O16−0.2271.0454 C11H11O8
211.0244 C9H7O6
169.0145 C7H5O5
125.0249 C6H5O3
Not identified
910.60817.39911−C43H62O15−0.6655.3483 C37H51O10
611.3593 C36H51O8
169.0139 C7H5O5
Quercotriterpenosid [11,12]Triterpenoids
derived from
quercotriterpenosids
9.95979.45181−C49H72O20−0.6817.4011 C43H61O15
755.4010 C42H59O12
Quercotriterpenosid derivative
[11,12]
1010.90679.36801−C36H56O12−0.5517.3207 C30H45O7
455.3171 C29H43O4
437.3075 C29H41O3
Glucosylated
bartogenic acid
[43,44]
Triterpenoids
derived from
bartogenic acid
14.30669.32721−C37H50O11−2.8517.3173 C30H45O7
455.3171 C29H43O4
437.3060 C29H41O3
Galloyl bartogenic acid
[43,44]
1113.501373.72481−C72H110O25−0.1695.3649 C36H55O13
647.3443 C35H51O11
485.2909 C29H41O6
471.3113 C29H43O5
Roburosid B or C
[43,44]
Triterpenoids
13.801403.73411−C72H110O240.2679.3705 C36H55O12
517.3173 C30H45O7
Roburosid A
[43,44]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Buche, G.; Colas, C.; Fougère, L.; Destandau, E. Oak Species Quercus robur L. and Quercus petraea Liebl. Identification Based on UHPLC-HRMS/MS Molecular Networks. Metabolites 2021, 11, 684. https://doi.org/10.3390/metabo11100684

AMA Style

Buche G, Colas C, Fougère L, Destandau E. Oak Species Quercus robur L. and Quercus petraea Liebl. Identification Based on UHPLC-HRMS/MS Molecular Networks. Metabolites. 2021; 11(10):684. https://doi.org/10.3390/metabo11100684

Chicago/Turabian Style

Buche, Gaëlle, Cyril Colas, Laëtitia Fougère, and Emilie Destandau. 2021. "Oak Species Quercus robur L. and Quercus petraea Liebl. Identification Based on UHPLC-HRMS/MS Molecular Networks" Metabolites 11, no. 10: 684. https://doi.org/10.3390/metabo11100684

APA Style

Buche, G., Colas, C., Fougère, L., & Destandau, E. (2021). Oak Species Quercus robur L. and Quercus petraea Liebl. Identification Based on UHPLC-HRMS/MS Molecular Networks. Metabolites, 11(10), 684. https://doi.org/10.3390/metabo11100684

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

Article Metrics

Back to TopTop