Evolutionary patterns of volatile terpene emissions
across 202 tropical tree species

Elodie A Courtois, Kyle G. Dexter, Charles Eliot Timothy Paine, Didier

Stien, Julien Engel, Christopher Baraloto, Jérôme Chave

► To cite this version:

Elodie A Courtois, Kyle G. Dexter, Charles Eliot Timothy Paine, Didier Stien, Julien Engel, et al. Evolutionary patterns of volatile terpene emissions across 202 tropical tree species. Ecology and Evolution, 2015, 6 (9), pp.2854-2864. 10.1002/ece3.1810. hal-01296640

HAL Id: hal-01296640

https://hal.sorbonne-universite.fr/hal-01296640

Submitted on 1 Apr 2016

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Evolutionary patterns of volatile terpene emissions across 202 tropical tree species

Elodie A. Courtois<fnref n="1" />,<fnref n="2" />,<fnref n="3" />, Kyle G. Dexter<fnref n="2" />,<fnref n="4" />,<fnref n="5" />, Charles Eliot Timothy Paine<fnref n="6" />, Didier Stien<fnref n="3" />,<fnref n="7" />, Julien Engel<fnref n="8" />, Christopher Baraloto<fnref n="8" />,<fnref n="9" /> & Jérôme Chave<fnref n="2" />

**Keywords**

Chemical defense, French Guiana, herbivory, secondary metabolites, tropical forest.

**Correspondence**

Elodie A. Courtois, CNRS Guyane, USR 3456, 2, Avenue Gustave Charlery, 97300 Cayenne, France.

Tel: 00594(0)5 94 29 92 98;

Fax: 00594 (0)5 94 35 65 22;

E-mail: courtoiselodie@gmail.com

**Funding Information**

We gratefully acknowledge financial support from Agence Nationale pour la Recherche (BRIDGE, ANR-Biodiversité program) and from "Investissement d'Avenir" grants managed by Agence Nationale de la Recherche (CEBA, ref. ANR-10-LABX-25-01; TULIP: ANR-10-LABX-0041; ANAEE-France: ANR-11-INBS-0001). KGD was funded by a CNRS post-doctoral fellowship during the time this research was completed.

Received: 18 May 2015; Accepted: 4 October 2015

doi: 10.1002/ece3.1810

**Abstract**

Plant responses to natural enemies include formation of secondary metabolites acting as direct or indirect defenses. Volatile terpenes represent one of the most diverse groups of secondary metabolites. We aimed to explore evolutionary patterns of volatile terpene emission. We measured the composition of damage-induced volatile terpenes from 202 Amazonian tree species, spanning the angiosperm phylogeny. Volatile terpenes were extracted with solid-phase micro extraction and desorbed in a gas chromatography-mass spectrometry for compound identification. The chemical diversity of the terpene blend showed a strong phylogenetic signal as closely related species emitted a similar number of compounds. Closely related species also tended to have compositionally similar blends, although this relationship was weak. Meanwhile, the ability to emit a given compound showed no significant phylogenetic signal for 200 of 286 compounds, indicating a high rate of diversification in terpene synthesis and/or great variability in their expression. Three lineages (Magnoliales, Laurales, and Sapindales) showed exceptionally high rates of terpene diversification. Of the 70 compounds found in >10% of their species, 69 displayed significant correlated evolution with at least one other compound. These results provide insights into the complex evolutionary history of volatile terpenes in angiosperms, while highlighting the need for further research into this important class of compounds.

**Introduction**

The interactions of plants with pollinators, herbivores, and pathogens are mediated by the chemical compounds they emit (Hartmann 2007). Such compounds have been historically defined as secondary metabolites, that is

chemical compounds that are not used directly for growth or reproduction (Fraenkel 1959). They are characterized in plants by a large structural diversity, with over 200,000 compounds identified to date, with many more expected to be discovered (Hartmann 2007). Terpenes represent one of the largest groups of secondary metabolites, with

<footnote n="1">¹CNRS Guyane, USR 3456, 2, Avenue Gustave Charlery, 97300 Cayenne, France</footnote>

<footnote n="2">²Laboratoire Evolution et Diversité Biologique, UMR 5174 CNRS/Université Paul Sabatier, 118, route de Narbonne, 31062 Toulouse, France</footnote>

<footnote n="3">³CNRS, UMR Ecofog, Université Antilles Guyane, BP 709, 97387 Kourou Cedex, France</footnote>

<footnote n="4">⁴School of GeoSciences, University of Edinburgh, Edinburgh EH9 3FF, UK</footnote>

<footnote n="5">⁵Royal Botanic Garden Edinburgh, 20a Inverleith Row, Edinburgh EH3 5LR, UK</footnote>

<footnote n="6">⁶Biological and Environmental Sciences, University of Stirling, Stirling FK9 4LA, UK</footnote>

<footnote n="7">⁷Laboratoire de Biodiversité et Biotechnologies Microbiennes (LBBM), Observatoire Océanologique, Sorbonne Universités, UPMC Univ Paris 06, CNRS, 66650 Banyuls-sur-mer, France</footnote>

<footnote n="8">⁸INRA, UMR Ecofog, BP 701, 97387 Kourou Cedex, France</footnote>

<footnote n="9">⁹International Center for Tropical Botany, Department of Biological Sciences, Florida International University, Miami, Florida 33199</footnote>
approximately 25,000 structures reported to date (Gershzenzon and Dudareva 2007). Some terpenes (monoterpenes and sesquiterpenes) are volatile and can be stored in specialized organs in leaves, stems, and trunks and released while others are synthesized de novo after wounding (Fineschi et al. 2013; Lange 2015). The main role of these compounds is to defend against pathogens and herbivores (Unsicker et al. 2009) by direct toxicity, repulsion of herbivores, or the attraction of herbivores' enemies (Heil 2008). The blend of volatile compounds emitted by a given plant species can be highly complex, including up to 200 different compounds (Dicke and van Loon 2000).

A family of enzymes named terpene synthases (TPSs) is responsible for the synthesis of isoprene, monoterpenes, and sesquiterpenes from two isomeric 5-carbon precursor “building blocks” (Dudareva et al. 2004; Chen et al. 2011). TPS diversity appears to originate from repeated duplication and subsequent divergence of an ancestral TPS gene family important in primary metabolism (Bohlmann et al. 1998). Minor changes in the structure of TPSs can lead to major changes in product profiles, indicating that the rich terpenoid genetic diversity could constitute a means of rapid evolutionary adaptation to novel biotic interactions (Memari et al. 2013). The genomic architecture of TPS expression is complex (Chen et al. 2011) and the expression of a particular TPS usually results in multiple products (Bohlmann et al. 1998). Biosynthetic pathways for terpene production are present in the extant lineages that are most closely related to flowering plants (Chen et al. 2011), so it is reasonable to assume that they arose early in the evolution of land plants (Croteau et al. 2000; Wink 2003).

Patterns and processes of terpene diversification in plants remain subject to debate (Hanson et al. 1999; Harley et al. 1999; Fineschi et al. 2013). There may be several selective advantages to maintaining a large diversity of volatile terpenes in individuals (Gershzenzon and Dudareva 2007). First, in a world where multiple natural enemies are encountered, it could provide a more potent defense and biological activity against a spectrum of aggressors (Firn and Jones 2006). There is evidence, for example, that different monoterpenes in thyme (*Thymus*) deter different species of herbivores and inhibit different pathogens (Linhart and Thompson 1999), suggesting that a more diverse blend of compounds can protect against a broader range of enemies. Second, the blend of volatiles emitted by a plant may be a more efficient deterrent than the sum of its parts (Gershzenzon and Dudareva 2007), as several studies point out that mixtures are more active than pure compounds (Lin et al. 1987). It is therefore expected that species or lineages that invest more in this mechanism of defense should display both a higher

diversity in the emitted blend at the species level and a higher rate of diversification of compounds in the lineage.

Tropical tree species produce a large diversity of volatile terpenes, but the diversity of compounds can differ markedly among genera and families (Courtois et al. 2009). There has been considerable recent interest in exploring how sympatric species within the same assemblage partition ecological resources and functional trait space, as one of the major mechanisms of species coexistence may be niche partitioning (Cavender-Bares et al. 2009). The great diversity of moist tropical forests has motivated much of this recent research. Most recent empirical studies of niche partitioning in tropical trees have focused on leaf-level and trunk-level traits related to resource acquisition (Kraft and Ackerly 2010; Baraloto et al. 2012). However, the coexistence of diverse assemblages of tropical trees may be attributable to negative density-dependent effects mediated by natural enemies such as pathogens and herbivores (Harms et al., 2000, Wright 2002; Novotny et al. 2010), which would involve partitioning of defense niches, rather than resource-use niches. Even so, very few studies have explored the partitioning of defense niches, particularly chemical defense niches (but see Kursar et al. 2009).

In this study, we combine an original dataset of volatile terpene composition across 202 tropical tree species with a phylogenetic hypothesis for those species to explore patterns of volatile terpene evolution in tropical trees. More specifically, we test whether (1) the diversity (i.e., total number) and composition of compounds shows significant phylogenetic signal, *that is* whether closely related species tend to share a more similar blend of compounds in terms of the total number of emitted compounds and identity of compounds, (2) individual compounds display phylogenetic signal, (3) the pattern of diversification of compounds is the same in different lineages, and (4) compounds display significant correlated evolution with each other across the phylogeny.

# Materials and Methods

## Collection of plant tissue

Field sampling was conducted in nine 1-ha permanent tree plots that span coastal French Guiana (Baraloto et al. 2012). In each plot, all trees >10 cm diameter at breast height (d.b.h.) were mapped, measured, and collected for herbarium vouchers (deposited at the Cayenne herbarium) and silica gel-dried leaf DNA vouchers. All plants were subsequently identified by taxonomic specialists (Baraloto et al. 2012). For analysis of VOC composition, we selected 202 species across 172 genera,
representing nearly half of the 380 tree genera present in the Guiana Shield (de Granville 1988) and spanning the breadth of the angiosperm phylogenetic tree. A previous study on a subset of the same species revealed low to absent intraspecific variation in volatile terpene presence/absence (Courtois et al. 2009). Of the 44 species with more than three sampled individuals, only 8 displayed an ambiguous pattern in an unsupervised clustering analysis (18%; see fig. 3 of Courtois et al. 2009), with no hard conflicts. For this study, we therefore decided to include a single individual per species, which allowed us to include rare species and span a broader range of chemical diversity.

For each sampled tree, we collected ~20 mg of fresh leaf tissue and a 1-cm² piece of bark at 1 m above the ground using a leather punch. Each sample was immediately placed in a glass vial (10 mL) and sealed with a screw cap containing a Teflon-lined septum (Varian Instruments, Sunnyvale, CA, USA). Sealed vials were maintained at -4°C in a portable field freezer until arrival at the laboratory, where they were stored at -20°C until analysis.

## Phylogenetic reconstruction

We pruned the phylogenetic tree described in Baraloto et al. (2012) to species for which VOCs were sampled. Full details regarding laboratory protocols and phylogeny construction are provided by Baraloto et al. (2012), but we provide the main details here. Up to 30 mg of dry leaf tissue was ground for 2 min in a TissueLyser mixer-mill disruptor (Qiagen, Venlo, Netherlands) using tungsten beads. Lysis incubation using CTAB 1% PVP buffer was carried out at 65°C for ≥2 h for cambium tissue and 1 h for leaf tissue. Total DNA extraction was performed using a Biosprint 15 workstation (Qiagen) following the manufacturer's protocols. PCR amplification and sequencing of sections of the *rbcL* and *matK* chloroplast genes was then conducted. A total of 606 *rbcL* sequences (1320 bp), and 244 *matK* sequences (891 bp) were obtained, with at least one sequence per genus for *matK*. These sequences were deposited in GenBank/EBI (accession numbers JQ625717–JQ626579).

A dated phylogenetic hypothesis was constructed as follows. First a maximum-likelihood phylogenetic tree was generated using RAxML v7.0.1 software (Stamatakis et al. 2008), and this tree was made ultrametric using Sanderson's (1997) nonparametric rate smooth (NPRS) method. The ultrametric tree was then used as the starting tree in BEAST v1.5.3 (Drummond & Rambaut 2007). We implemented an uncorrelated relaxed molecular clock, which draws substitution rates from a log-normal distribution independently for each branch in the phylogenetic tree. In

order to calibrate the phylogenetic tree, we used fossil ages to place constraints on divergence times for minimum crown age at nine nodes spread across the phylogeny (see Appendix S2 in Baraloto et al. 2012). We ensured that the runs had converged after 30 million generations. Using the ultrametric tree, we calculated the phylogenetic distance between all pairs of species as the total branch length separating them.

## Chemical analysis

Volatiles were extracted with a solid-phase micro extraction (SPME) fiber of the type PDMS/DVB (polydimethylsiloxane/divinylbenzene) 65 µm (Supelco, Bellefonte, PA). Fibers were conditioned before the first use for 30 min at 250°C, following instructions from the manufacturer. Before extraction, the glass vials containing the tissue samples were maintained at room temperature for at least 1 h. The SPME fiber was placed into the vial with the tissue sample (bark or leaf) for 5-60 min at ambient temperature (25°C, exposure time was optimized for each species to maximize the extraction without saturating the analytical column). The fiber was inserted immediately into the 250°C inlet of a Varian 3800 Gas Chromatograph (GC) fitted with a Saturn 2000 ion-trap Mass Spectrometer (MS; Varian Instruments, Sunnyvale, CA). The GC was run with a nonpolar Varian DB-5 column (30 m × 0.25 mm ID, 0.25 µm film), and helium was the carrier gas at a constant flow of 1 mL/min. The oven temperature program of the GC started at 50°C, with 6°C/min temperature increase up to 140°C, and then with 4°C/min increments up to 160°C. This temperature was held for 1 min and increased finally to 200°C at 10°C/min. The MS was operated in electron impact (EI) mode at 70 eV, with a scan range of 30-450 m/z. After each analysis, the fibers were cleaned in the injector port for 10 min at 250°C, and each fiber was reused no more than 100 times. Control analyses (blanks) were performed every ten analyses to check for contamination of the fiber.

We used a subsample of species (*Protium* sp. Burseraceae, *Inga* sp. Mimosaceae, *Guarea* sp. Meliaceae, and *Spondias mombin* Anacardiaceae) to confirm that storage of leaf or bark tissue samples at -20°C prior to analysis did not result in drastic changes in the composition of volatiles (e.g., due to frost damage). For further details on chemical protocols, see Courtois et al. (2009).

For chemical composition, we inferred the presence/absence of chemicals from the chromatograms using a novel statistical approach implemented in the MSeasy package (Nicolè et al. 2012) in the R Statistical Environment (http://cran.r-project.org/), which is both more efficient and more accurate than visual classification procedures
when conducting a large number of analyses. This method is based on the mass spectrum of each detected molecule and the corresponding Kováts Retention Index (Kováts 1958). Molecules were identified based on the comparison of mass spectra with standards or with the NIST 98 MS library, the ADAMS library (Adams 2007), and with RI reported in the literature. Unidentified molecules were defined as morphomolecules characterized by their mass spectrum and their RI and assigned to a class of compounds (monoterpenes or sesquiterpenes) based on ion fragmentation. Overlapping peaks were checked manually. We defined the chemical composition for each individual tree as the union of all compounds found in the leaves or bark.

## Evolutionary analyses

### Measures of phylogenetic signal of the overall blend of compounds

The diversity (number of compounds of the emitted blend) can be an important component of efficiency of defense. We estimated phylogenetic signal for the total number of compounds emitted by each species using Pagel's lambda (Pagel 1999). This metric varies from 0 to 1, with 0 indicating no phylogenetic signal and 1 indicating high phylogenetic signal, or a perfect correlation of a trait with phylogeny under a Brownian motion model. We used likelihood ratio tests to determine whether the estimate of phylogenetic signal was significantly greater than zero. Moreover, we assessed whether closely related species tend to share a more similar blend of compounds using a Mantel test of the correlation between chemical dissimilarity and species pairwise phylogenetic distances. Chemical dissimilarity was measured by the pairwise Manhattan distance and the number of compounds that are present in one species but not in the other (see Legendre and Cáceres 2013).

### Measures of phylogenetic signal of individual compounds

We tested for phylogenetic signal for each compound found in more than one species using the D metric for binary traits developed by Fritz and Purvis (2010). The D statistic is equal to 1 if the observed binary trait has a phylogenetically random distribution across the tips of the phylogeny and to 0 if the observed trait is phylogenetically clustered as if it had evolved by Brownian motion (null expectation). Using 1000 permutations, a P-value is generated to test whether the observed value of D is significantly different from 1 (a random distribution) or 0 (the null expectation).

## Character reconstruction to estimate pattern of diversification in lineages

We used two approaches to reconstruct the evolutionary history of all volatile terpenes found in at least two species. For comparability with previous studies, we used maximum parsimony reconstruction, hereafter MPR (Swofford and Maddison 1987). This approach unrealistically assumes that only one change can happen on a single phylogenetic branch, no matter how long that branch is, and it underestimates uncertainty in ancestral state assignments (Bollback 2006). We therefore additionally implemented stochastic character mapping, hereafter SCM (Nielsen 2002; Huelsenbeck et al. 2003). SCM uses a probabilistic approach to reconstruct character states based on a specified model of character evolution and allows characters to change more than once along branches in a phylogenetic tree (Huelsenbeck et al. 2003). For both approaches, we determined the minimum number of gains of each compound, either across all maximally parsimonious reconstructions, in the case of MPR, or across 100 simulations, in the case of SCM. We also determined the state of compounds at the root node and considered a compound as ancestral of all taxa if it was reconstructed as unambiguously present under MPR or as present with a probability greater than 0.9 under SCM (i.e., reconstructed as present in >90% of the simulations). For two-state character reconstruction using SCM, two models can be used. In the ARD (all rates different) model, the rate of change from one state to the other can be different from the reverse change; while under ER (equal rate) model, the two rates are constrained to be equal. For each compound, we compared the likelihood of these two models with a likelihood ratio test (see below for details) and used the most likely model for simulations (Table S1).

## Correlated evolution analysis

We conducted analyses to determine the degree of correlated evolution among compounds. For the compounds that were found in more than 10% of species (i.e., over 20 species), we tested for correlated evolution using the BayesDiscrete module implemented in the BayesTraits software (Pagel and Meade 2006). Correlated evolution is tested by comparing the fit (log-likelihood) of two models: one model with independent evolution of the two traits and one with correlated evolution of the traits. In order to test which model is a better fit to the data, we computed a likelihood ratio (LR) as $LR = 2 \times (\ln_{ML2} - \ln_{ML1})$, where $\ln_{ML1}$ is the log-likelihood of the independent model and $\ln_{ML2}$ is the log-likelihood of the correlated model. The LR
follows a chi-squared distribution with four degrees of free-
dom. When the log-likelihood of the correlated evolution
model is significantly greater than the log-likelihood of the
independent model, there is support for correlated evolu-
tion of the two traits.

Evolutionary analyses were conducted in the R Statisti-
cal Environment (R Core Development Team) using the
following packages: caper (Orme 2012), ape (Paradis et al.
2004), phytools (Revell 2012), and phangorn (Schliep
2011), except for the correlated evolution analyses, which
were conducted with the BayesTraits software (Pagel and
Meade 2006).

Results

All 202 species assessed emitted at least one volatile ter-
pene compound. Yet, we observed high variation in the
diversity of emitted volatile terpenes, with species emit-
ting between 1 and 30 monoterpenes (9.2 ± 4.5 SD) and
between 0 and 59 sesquiterpenes (21.5 ± 14.3 SD). Over-
all, 376 terpenes were detected across the 202 species, 88
being monoterpenes and 288 being sesquiterpenes.
Among these, 59 monoterpenes (67%) and 114 sesquiter-
penes (40%) could be identified based on the comparison
of mass spectra with standards or with libraries and
retention indices reported in the literature. The propor-
tion of compounds found in a single species (i.e., single-
tons) was twice as high for monoterpenes (39%; 34 of 88
compounds) as for sesquiterpenes (19%; 56 of 288
compounds).

The total number of compounds found in each species exhibited high and significant phylogenetic signal ($\lambda$ = 0.65, *P* < 0.01), although this pattern was stronger for sesquiterpenes ($\lambda$ = 0.71, *P* < 0.01) than for monoterpenes ($\lambda$ = 0.14, *P* = 0.04). The composition of the overall blend of compounds showed weak but significant phylogenetic signal (Mantel *r* = 0.15, *P* < 0.001). As with the total number of compounds, the pattern with respect to sesquiterpenes was stronger (sesquiterpenes: Mantel *r* = 0.26, *P* < 0.001; monoterpenes: Mantel *r* = 0.11, *P* < 0.001, Fig. 1).

When considering each compound separately, 75 of the
232 sesquiterpenes that were found in more than one spe-
cies (32%, Table S1) showed significant phylogenetic sig-
nal, whereas only 11 of the 54 monoterpenes found in
more than one species (20%, Table S1) showed significant
phylogenetic signal.

Of the 286 compounds emitted by more than one spe-
cies, only seven compounds, including five monoterpenes
($\alpha$-pinene, limonene, $\beta$-pinene, camphene, and $p$-cymene)
and two sesquiterpenes ($\beta$-caryophyllene and $\alpha$-copaene),
were reconstructed as present at the root of the phylogeny
by both methods (SCM – stochastic character mapping

<figure><img src="image_2.png" /></figure>

<figure><img src="image_3.png" /></figure>

**Figure 1.** Representation of chemical dissimilarity for (A) monoterpenes and (B) sesquiterpenes vs. phylogenetic dissimilarity. Each point represents a pairwise comparison of two species characterized by their Manhattan chemical dissimilarity (number of compounds that are present in one species but not in the other) and their phylogenetic dissimilarity (computed using branch length, distance in My among the two species). The dotted red lines represent the regression line between the two values.

and MPR – maximum parsimony reconstruction,
Table S1). An additional four compounds (2 sesquiterpe-
nes; $\alpha$-cubebene and $\delta$-cadinene, and 2 monoterpenes;
3-carene and $\alpha$-thujene) were inferred to be present at the
root of the phylogeny using MPR (Table S1). Under both
methods, nearly all compounds (284 of 286 using SCM
and 280 using MPR) had a minimum number of gains
greater than 1 (Table S1).
The mean number of monoterpenes and sesquiterpenes per species differed among taxonomic orders (Kruskal Wallis, *P* = 0.008 and *P* < 0.001, respectively). Species with the highest number of monoterpenes and sesquiterpenes belonged to the Magnoliales, Laurales and Sapindales (Table 1, Fig. 2). These three orders also contained the highest number of unique compounds (in addition to the Malphigiales; Table 1). Overall, 102 compounds among 376 were restricted to a single order and most of those (90 of 102, 88%) are singletons, that is found in only one species.

Of the 15 monoterpenes and the 55 sesquiterpenes that are found in more than 10% of the sampled species, all but one (α-pinene) displayed significant correlated evolution with at least one other compound (Table S2 and Table S3). Moreover, 67% of these monoterpenes (10 of 15) and 45% of these sesquiterpenes (25 of 55) evolved in a correlated manner with at least five other compounds.

Discussion

The emission of volatile terpenes represents an important plant defense strategy, but little is known about its evolutionary history in Neotropical lineages. In this study, we determined the composition of the emitted volatile terpene blend for 202 tropical tree species and examined patterns of compound evolution, diversification, and correlated evolution. We found that the diversity of the emitted blend (number of compounds) in individual

species showed high phylogenetic signal and that the composition of the overall blend showed weak, albeit significant phylogenetic signal. Indeed, most of the individual compounds showed no significant phylogenetic signal. Certain lineages of plants showed significantly higher rates of compound diversification. Finally, we found a strong pattern of correlated evolution for common compounds (those found in >10% of species), with most such compounds displaying correlated evolution.

Significant phylogenetic signal of volatile blends could be explained by Brownian motion (i.e., random drift) evolution, phylogenetic conservatism, or evolutionary constraint in compound evolution (Losos 2008). Conversely, finding weaker phylogenetic signal than expected under Brownian motion suggests that diversifying selection for terpene composition may be occurring. The fact that monoterpenes displayed a weaker phylogenetic signal than sesquiterpenes might therefore suggest either a stronger diversifying selection in this group or a stronger constraining selective pressure on sesquiterpenes. Such results could also suggest, as have other studies, that the characteristic blend of volatiles emitted by wounded leaves might carry information about species identity (Courtois et al. 2009; Pearse et al. 2013) and phylogenetic relationships (Pearse and Hipp 2009). Most previous studies reporting phylogenetic conservatism in the blend of volatiles emitted by wounded leaves were based on within-genus comparisons (Pearse et al. 2013). Our study generalizes these findings across a broader taxonomic and phylogenetic range.

<table><caption>**Table 1.** Summary of the number of volatile compounds found in angiosperm orders, sorted by the mean number of terpenes per species. Note that the number of species examined differed among orders.</caption><thead><tr><th></th><th>No. of species analyzed</th><th>Mean no. of monoterpenes per species</th><th>Mean no. of sesquiterpenes per species</th><th>Mean no. of terpenes per species</th><th>No. of unique monoterpenes</th><th>No. of unique sesquiterpenes</th><th>No. of unique terpenes</th></tr></thead><tbody><tr><td>Brassicales</td><td>1</td><td>6</td><td>6</td><td>12</td><td>0</td><td>0</td><td>0</td></tr><tr><td>Lamiales</td><td>3</td><td>4.7</td><td>10.7</td><td>15.4</td><td>0</td><td>0</td><td>0</td></tr><tr><td>Santalales</td><td>5</td><td>8.2</td><td>11.4</td><td>19.6</td><td>0</td><td>0</td><td>0</td></tr><tr><td>Fabales</td><td>28</td><td>8</td><td>13.1</td><td>21.1</td><td>0</td><td>3</td><td>3</td></tr><tr><td>Rosales</td><td>8</td><td>8</td><td>14</td><td>22</td><td>0</td><td>2</td><td>2</td></tr><tr><td>Malvales</td><td>8</td><td>5.9</td><td>18</td><td>23.9</td><td>2</td><td>2</td><td>4</td></tr><tr><td>Ericales</td><td>20</td><td>8.8</td><td>16.5</td><td>25.3</td><td>3</td><td>4</td><td>7</td></tr><tr><td>Malpighiales</td><td>39</td><td>8.7</td><td>16.8</td><td>25.5</td><td>6</td><td>11</td><td>17</td></tr><tr><td>Aquifoliales</td><td>1</td><td>10</td><td>16</td><td>26</td><td>0</td><td>0</td><td>0</td></tr><tr><td>Celastrales</td><td>2</td><td>8</td><td>19</td><td>27</td><td>0</td><td>0</td><td>0</td></tr><tr><td>Gentianales</td><td>15</td><td>9.5</td><td>18.1</td><td>27.6</td><td>6</td><td>3</td><td>9</td></tr><tr><td>Myrtales</td><td>12</td><td>8.9</td><td>19</td><td>27.9</td><td>2</td><td>0</td><td>2</td></tr><tr><td>Garryales</td><td>3</td><td>11</td><td>20</td><td>31</td><td>0</td><td>2</td><td>2</td></tr><tr><td>Caryophyllales</td><td>2</td><td>7.5</td><td>24</td><td>31.5</td><td>0</td><td>0</td><td>0</td></tr><tr><td>Oxalidales</td><td>1</td><td>8</td><td>24</td><td>32</td><td>0</td><td>1</td><td>1</td></tr><tr><td>Sapindales</td><td><strong>26</strong></td><td><strong>10.7</strong></td><td><strong>33.8</strong></td><td><strong>44.5</strong></td><td><strong>4</strong></td><td><strong>16</strong></td><td><strong>20</strong></td></tr><tr><td>Laurales</td><td><strong>12</strong></td><td><strong>12.8</strong></td><td><strong>36.1</strong></td><td><strong>48.9</strong></td><td><strong>3</strong></td><td><strong>14</strong></td><td><strong>17</strong></td></tr><tr><td>Magnoliales</td><td><strong>16</strong></td><td><strong>11.7</strong></td><td><strong>40.3</strong></td><td><strong>52</strong></td><td><strong>10</strong></td><td><strong>8</strong></td><td><strong>18</strong></td></tr></tbody></table>

Bold values represent the three orders with the highest mean number of terpenes per species.
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Phylogenetic tree described in Baraloto et al. (2012) of the 202 studied species with each terminal branches colored by the number of terpenes found in the descendent node from green (low number of compounds) to red (large number of compounds). Laurales, Magnoliales, and Sapindales are highlighted.</figcaption></figure>

When considering each compound independently, we found that few (20% of monoterpenes and 32% of sesquiterpenes) displayed phylogenetic signal. The few studies that have tested for phylogenetic signal of terpene emissions in a broad set of plant species have yielded the same result (e.g., Llusià et al. 2010, 2014). In these papers, the authors used more than one individual per species and highlighted the importance of both intraspecific variation and quantification of compounds, two aspects that are not taken into account in our study, yet they found similar results. The observed low phylogenetic signal could be explained by several characteristics of terpenes. It can reflect convergent evolution of compounds, known to be frequent for this class of molecules (Bohlmann et al. 1998). For example, limonene-producing enzymes originated independently in gymnosperms and

Lamiaceae (Bohlmann et al. 1998). However, it is likely that in some cases the pathway is not lost but merely “switched off,” *that is* plants retain the ability to switch them on again (Wink 2003). For example, in *Arabidopsis thaliana*, two ecotypes differ in their emission of (*E*)-β-ocimene, and this difference is due to a mutation that inactivates one allele of the (*E*)-β-ocimene-producing synthase in one of the ecotypes (Huang et al. 2010). Moreover, it is known that VOC emissions can be regulated as elicitors released by herbivores (Mattiacci et al. 1994). Hence, the blend measured after the mechanical damage we executed in this study may not lead to the release of all the compounds that a given species can emit. Our sampling technique may therefore have led us to underestimate the actual ability of the species to express compounds. Thus, evolutionary and/or phenotypic changes in
gene regulation may underlie our finding of nonsignifi-
cant phylogenetic signal for most compounds rather than
high rates of compound evolution per se. Finally, we
acknowledge that terpene synthases usually synthesize a
spectrum of compounds and that a given compound can
be synthesized by multiple synthases (Memari et al.
2013), which could be expressed in response to various
environmental and biological cues (Penuelas and Llusia
2001). Therefore, even if the patterns we found here are
strong, it is difficult to directly relate them to evolution-
ary changes in biosynthetic pathways.

We found that the diversity of the emitted blend
showed strong phylogenetic signal and that diversification
of volatile terpenes occurred preferentially in certain plant
lineages. Three orders (Magnoliales, Laurales, and Sapin-
dales) harbored species tending to display a large number
of compounds. These same three orders also displayed
the greatest number of compounds at the order level and
the highest number of compounds that are restricted to a
single order. This suggests that the capacity to generate
novel compounds is not evenly distributed across the
angiosperm phylogeny. From an ecological standpoint, a
greater diversity of the emitted blend likely provides a
better protection against a broader range of enemies (Ger-
shenzon and Dudareva 2007). Such greater diversity in
the emitted blend at species level and a higher rate of
diversification of compounds in these lineages might
therefore indicate that their defense strategy is focused on
VOC synthesis and emissions. Other lineages may rely
upon other means of defense against natural enemies
(e.g., latex production, tough leaves).

We also found that at least two major events of diversi-
fication for volatile terpenes may have occurred in tropi-
cal trees: one corresponding to the Magnoliids
(Magnoliales and Laurales), for which the stem age is esti-
mated at 122–125 million years (Bell et al. 2010); and the
other within the Sapindales for which divergence is sug-
gested at ~ 70 million years (Magallón and Castillo 2009).
Of course, we did not include all plant orders in our
study, and similar diversification events may have
occurred in other lineages that we have not sampled (e.g.,
those in the temperate zone).

Only seven compounds are reconstructed as unambigu-
ously present at the root of the phylogeny (five monoterpe-
nes and two sesquiterpenes). Given the high rate of
changes in terpene expression, it is not surprising that we
fail to unambiguously reconstruct ancestral states of most
compounds deep in the phylogeny. Nonetheless, the fact
that, even with this limitation, the unambiguous recon-
struction of these seven compounds at the root of the
phylogeny (with two different statistical approaches) could
indicate that they are under more selective pressure to be
retained than other compounds. In fact, among the five

putatively ancestral monoterpenes, α-pinene and limonene
have been highlighted as compounds with the highest
known repellent activity against insects (Nerio et al. 2010).
They are also emitted by gymnosperms. Within the
sesquiterpenes, the two compounds that can be considered
as ancestral (β-caryophyllene and α-copaene) were already
found to be frequent in tropical tree species (Courtois et al.
2009) and they are known to be implicated in direct and
indirect defense against herbivores (Gols et al. 2008).

Most of the common compounds show evidence of
correlated evolution, which is not unexpected given that
most volatile terpenes are produced by only a few biosyn-
thetic pathways (Chen et al. 2011). Some of the detected
compounds may have no function but be present because
they covary with a compound used in defense (Chen
et al. 2011).

The emissions of terpenes have important effects on atmospheric chemistry and climate regulation (Fehsenfeld et al. 1992). Every species assessed in this study was found to emit terpenes, as found in other studies (Owen et al. 2001; Llusia et al. 2014). While we did not estimate the quantity of emitted terpenes, our study suggests that terpene emissions might be more important than previously thought in tropical trees and that forests differing in floristic composition will have important differences in volatile emissions. Our results therefore provide an important first step to parameterize models examining the effect of these compounds on atmospheric processes at broad spatial scales.

The high evolutionary variability of these compounds
and the pervasive correlated evolution suggest a complex
pattern of terpene biochemistry evolution. Nevertheless,
this study suggests two major phases of diversification in
volatile terpenes, one in the Magnoliids (Magnoliales and
Laurales) and one within the Sapindales. This might indi-
cate that volatile terpenes have been a major determinant
of the defense strategy in these lineages (Becerra et al.
2009). These lineages may therefore be considered a good
starting point for more mechanistic research into the reg-
ulation of terpene emission in tropical tree species.

Finally our study emphasizes the key role of chemistry
in mediating biotic interactions across trophic levels in the
tropics and suggests that the large diversity of compounds
across coexisting species within tropical tree communities
could be due, at least in part, to these biotic interactions.

Acknowledgments

We thank all participants of the BRIDGE project, Antoine
Stevens and the Institut Pasteur of French Guiana in
Cayenne for providing laboratory facilities. We gratefully
acknowledge financial support from Agence Nationale
pour la Recherche (BRIDGE, ANR-Biodiversité program)
and from "Investissement d'Avenir" grants managed by Agence Nationale de la Recherche (CEBA, ref. ANR-10-LABX-25-01; TULIP: ANR-10-LABX-0041; ANAEE-France: ANR-11-INBS-0001). KGD was funded by a CNRS Postdoctoral Fellowship during the time this research was completed.

## Conflict of Interest

None declared.

## References

Adams, R. P. 2007. Identification of essential oil components by gas chromatography/mass spectrometry (4). Allured publishing corporation.

Baraloto, C., O. J. Hardy, C. E. T. Paine, K. G. Dexter, C. Cruaud, L. T. Dunning, et al. 2012. Using functional traits and phylogenetic trees to examine the assembly of tropical tree communities. J. Ecol. 100:690-701.

Becerra, J. X., K. Noge, and D. L. Venable. 2009. Macroevolutionary chemical escalation in an ancient plant-herbivore arms race. Proc. Natl Acad. Sci. USA 106:18062-18066.

Bell, C. D., D. E. Soltis, and P. S. Soltis. 2010. The age and diversification of the angiosperms re-revisited. Am. J. Bot. 97:1296-1303.

Bohlmann, J., G. Meyer-Gauen, and R. Croteau. 1998. Plant terpenoid synthases: molecular biology and phylogenetic analysis. Proc. Natl Acad. Sci. 95:4126-4133.

Bollback, J. 2006. SIMMAP: stochastic character mapping of discrete traits on phylogenies. BMC Bioinformatics 7:88.

Cavender-Bares, J., K. H. Kozak, P. V. A. Fine, and S. W. Kembel. 2009. The merging of community ecology and phylogenetic biology. Ecol. Lett. 12:693-715.

Chen, F., D. Tholl, J.Bohlmann, and E. Pichersky. 2011. The family of terpene synthases in plants: a mid-size family of genes for specialized metabolism that is highly diversified throughout the kingdom. Plant J. 66:212-229.

Courtois, E. A., C. E. T. Paine, P. A. Blandinieres, D. Stien, J. M. Bessiere, E. Houel, et al. 2009. Diversity of the volatile organic compounds emitted by 55 species of tropical trees: a survey in French Guiana. J. Chem. Ecol. 35:1349-1362.

Croteau, R., T. M. Kutchan, and N. G. Lewis. 2000. Natural Products (Secondary Metabolites). Pp. 1250-1319. in Buchanan, Bob B., and Russell L. Jones, ed. Biochemistry & molecular biology of plants. Vol. 40. American Society of Plant Physiologists, Rockville.

de Granville, J. J. 1988. Phytogeographical characteristics of the Guianan forests. Taxon 37:578-594.

Dicke, M., and J. J. A. van Loon. 2000. Multitrophic effects of herbivore-induced plant volatiles in an evolutionary context. Entomol. Exp. Appl. 97:237-249.

Drummond, A. J., and A. Rambaut. 2007. BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol. Biol. 7:214.

Dudareva, N., E. Pichersky, and J. Gershenzon. 2004. Biochemistry of plant volatiles. Plant Physiol. 135:1893-1902.

Fehsenfeld, F., J. Calvert, R. Fall, P. Goldan, A. B. Guenther, C. N. Hewitt, et al. 1992. Emissions of volatile organic compounds from vegetation and the implications for atmospheric chemistry. Global Biogeochem. Cycles 6:389-430.

Fineschi, S., F. Loreto, M. Staudt, and J. Peñuelas. 2013. Diversification of volatile isoprenoid emissions from trees: evolutionary and ecological perspectives. Biology, controls and models of tree volatile organic compound emissions Springer, 1-20.

Firn, R. D., and C. G. Jones. 2006. Do we need a new hypothesis to explain plant VOC emissions? Trends Plant Sci. 11:112-113.

Fraenkel, G. S. 1959. Raison detre of secondary plant substances. Science 129:1466-1470.

Fritz, S. A., and A. Purvis. 2010. Selectivity in mammalian extinction risk and threat types: a new measure of phylogenetic signal strength in binary traits. Conserv. Biol. 24:1042-1051.

Gershenzon, J., and N. Dudareva. 2007. The function of terpene natural products in the natural world. Nat. Chem. Biol. 3:408-414.

Gols, R., L. M. A. Witjes, J. J. A. van Loon, M. A. Posthumus, M. Dicke, and J. A. Harvey. 2008. The effect of direct and indirect defenses in two wild brassicaceous plant species on a specialist herbivore and its gregarious endoparasitoid. Entomol. Exp. Appl. 128:99-108.

Hanson, D. T., S. Swanson, L. E. Graham, and T. D. Sharkey. 1999. Evolutionary significance of isoprene emission from mosses. Am. J. Bot. 86:634-639.

Harms, K. E., S. J. Wright, O. Calderón, A. Hernández, and E. A. Herre. 2000. Pervasive density-dependent recruitment enhances seedling diversity in a tropical forest. Nature 404:493-495.

Harley, P. C., R. K. Monson, and M. T. Lerdau. 1999. Ecological and evolutionary aspects of isoprene emission from plants. Oecologia 118:109-123.

Hartmann, T. 2007. From waste products to ecochemicals: fifty years research of plant secondary metabolism. Phytochemistry 68:2831-2846.

Heil, M. 2008. Indirect defence via tritrophic interactions. New Phytol. 178:41-61.

Huang, M., C. Abel, R. Sohrabi, J. Petri, I. Haupt, J. Cosimano, et al. 2010. Variation of herbivore-induced volatile terpenes among arabidopsis ecotypes depends on allelic differences and subcellular targeting of two terpene synthases, TPS02 and TPS03. Plant Physiol. 153:1293-1310.
Huelsenbeck, J. P., R. Nielsen, and J. P. Bollback. 2003. Stochastic mapping of morphological characters. *Syst. Biol.* 52:131–158.

Kováts, E. 1958. Gas-chromatographische charakterisierung organischer verbindungen. Teil 1: retentionsindices aliphatischer halogenide, alkohole, aldehyde und ketone. *Helv. Chim. Acta* 41:1915–1932.

Kraft, N. J., and D. D. Ackerly. 2010. Functional trait and phylogenetic tests of community assembly across spatial scales in an Amazonian forest. *Ecol. Monogr.* 80:401–422.

Kursar, T. A., K. G. Dexter, J. Lokvam, R. T. Pennington, J. E. Richardson, M. G. Weber, et al. 2009. The evolution of antiherbivore defenses and their contribution to species coexistence in the tropical tree genus Inga. *Proc. Natl Acad. Sci. USA* 106:18073–18078.

Lange, B. M.. 2015. The evolution of plant secretory structures and the emergence of terpenoid chemical diversity. *Annu. Rev. Plant Biol.* 66:139–159.

Legendre, P., and M. Cáceres. 2013. Beta diversity as the variance of community data: dissimilarity coefficients and partitioning. *Ecol. Lett.* 16:951–963.

Lin, S. Y., J. T. Trumble, and J. Kumamoto. 1987. Activity of volatile compounds in glandular trichomes of *Lycopersicon* species against two insect herbivores. *J. Chem. Ecol.* 13:837–850.

Linhart, Y. B., and J. Thompson. 1999. Thyme is of the essence: biochemical polymorphism and multi-species deterrence. *Evol. Ecol. Res.* 1:151–171.

Llusià, J., J. Peñuelas, J. Sardans, S. M. Owen, and Ü. Niinemets. 2010. Measurement of volatile terpene emissions in 70 dominant vascular plant species in Hawaii: aliens emit more than natives. *Glob. Ecol. Biogeogr.* 19:863–874.

Llusià, J., J. Sardans, Ü. Niinemets, S. M. Owen, and J. Peñuelas. 2014. A screening study of leaf terpene emissions of 43 rainforest species in Danum Valley Conservation Area (Borneo) and their relationships with chemical and morphological leaf traits. *Plant Biosyst. Int. Deal. Aspect. Plant Biol.* 148:307–317.

Losos, J. B. 2008. Phylogenetic niche conservatism, phylogenetic signal and the relationship between phylogenetic relatedness and ecological similarity among species. *Ecol. Lett.* 11:995–1003.

Magallón, S., and A. Castillo. 2009. Angiosperm diversification through time. *Am. J. Bot.* 96:349–365.

Mattiacci, L., M. Dicke, and M. A. Posthumus. 1994. Induction of parasitoid attracting synomone in brussels-sprouts plants by feeding of pieris-brassicae larvae – role of mechanical damage and herbivore elicitor. *J. Chem. Ecol.* 20:2229–2247.

Memari, H. R., L. Pazouki, and Ü. Niinemets 2013. The biochemistry and molecular biology of volatile messengers in trees. Pp. 47–93. in S. Fineschi, F. Loreto, M. Staudt and J. Peñuelas, eds. *Biology, controls and models of tree volatile organic compound emissions*. Springer, Netherlands.

Nerio, L. S., J. Olivero-Verbel, and E. Stashenko. 2010. Repellent activity of essential oils: a review. *Bioresour. Technol.* 101:372–378.

Nicolè, F., Y. Guitton, E. Courtois, S. Moja, L. Legendre, and M. Hossaert-McKey. 2012. MSeasy: unsupervised and untargeted GC-MS data processing. *Bioinformatics* 28:2278–2280.

Nielsen, R. 2002. Mapping mutations on phylogenies. *Syst. Biol.* 51:729–739.

Novotny, V., S. E. Miller, L. Baje, S. Balagawi, Y. Basset, L. Cizek, et al. 2010. Guild-specific patterns of species richness and host specialization in plant-herbivore food webs from a tropical forest. *J. Anim. Ecol.* 79:1193–1203.

Orme, C. D. L. 2012. *Caper: comparative analyses of phylogenetics and evolution in R. R package version 0.5.*

Owen, S. M., C. Boissard, and C. N. Hewitt. 2001. Volatile organic compounds (VOCs) emitted from 40 Mediterranean plant species: VOC speciation and extrapolation to habitat scale. *Atmos. Environ.* 35:5393–5409.

Pagel, M. 1999. Inferring the historical patterns of biological evolution. *Nature* 401:877–884.

Pagel, M., and A. Meade. 2006. Bayesian analysis of correlated evolution of discrete characters by reversible jump markov chain Monte Carlo. *Am. Nat.* 167:808–825.

Paradis, E., J. Claude, and K. Strimmer. 2004. APE: analyses of phylogenetics and evolution in R language. *Bioinformatics* 20:289–290.

Pearse, I. S., and A. L. Hipp. 2009. Phylogenetic and trait similarity to a native species predict herbivory on non-native oaks. *Proc. Natl Acad. Sci. USA* 106:18097–18102.

Pearse, I. S., W. S. Gee, and J. J. Beck. 2013. Headspace volatiles from 52 oak species advertise induction, species identity, and evolution, but not defense. *J. Chem. Ecol.* 39:90–100.

Penuelas, J., and J. Llusia. 2001. The complexity of factors driving volatile organic compound emissions by plants. *Biol. Plant.* 44:481–487.

Revell, L. J. 2012. phytools: an R package for phylogenetic comparative biology (and other things). *Methods Ecol. Evol.* 3:217–223.

Sanderson, M. J. 1997. A nonparametric approach to estimating divergence times in the absence of rate constancy. *Mol. Biol. Evol.* 14:1218–1231.

Schliep, K. P. 2011. Phangorn: phylogenetic analysis in R. *Bioinformatics* 27:592–593.

Stamatakis, A., P. Hoover, and J. Rougemont. 2008. A rapid bootstrap algorithm for the RAxML web servers. *Syst. Biol.* 57:758–771.

Swofford, D. L., and W. P. Maddison. 1987. Reconstructing ancestral character states under Wagner parsimony. *Math. Biosci.* 87:199–229.

Unsicker, S. B., G. Kunert, and J. Gershenzon. 2009. Protective perfumes: the role of vegetative volatiles in plant defense against herbivores. *Curr. Opin. Plant Biol.* 12:479–485.
Wink, M. 2003. Evolution of secondary metabolites from an
ecological and molecular phylogenetic perspective.
Phytochemistry 64:3–19.

Wright, S. J. 2002. Plant diversity in tropical forests: a review of mechanisms of species coexistence. Oecologia 130:1–14.

**Supporting Information**

Additional Supporting Information may be found in the
online version of this article:

**Table S1.** Measure of phylogenetic signal for each com-
pound using the D metric for binary traits with the *P*-

value for departure from a random distribution and
results from the SIMMAP simulation (ancestral state at
the root of the tree, minimum number of gain and mini-
mum number of losses).

**Table S2.** Coevolution of the 15 monoterpenes found in more than 20 species. The lower part of the matrix gives the LR (likelihood ratio) and the upper part of the matrix gives the *P*-value of the test.

**Table S3.** Coevolution of the 55 sesquiterpenes found in more than 20 species. The lower part of the matrix gives the LR (likelihood ratio) and the upper part of the matrix gives the *P*-value of the test.