plants Article Metabolome-Wide, Phylogenetically Controlled Comparison Indicates Higher Phenolic Diversity in Tropical Tree Species Guille Peguero 1,2,3, *,† , Albert Gargallo-Garriga 1,2,4,† , Joan Maspons 1,2,† , Karel Klem 4 , Otmar Urban 4 , Jordi Sardans 1,2 and Josep Peñuelas 1,2 1 CSIC, Global Ecology Unit CREAF-CSIC-UAB, E-08913 Bellaterra, Spain; albert.gargallo@gmail.com (A.G.-G.); j.maspons@creaf.uab.cat (J.M.); j.sardans@creaf.uab.cat (J.S.); josep.penuelas@uab.cat (J.P.) 2 Ecological and Forestry Applications Research Centre, Autonomous University of Barcelona, E-08913 Cerdanyola del Vallès, Spain 3 Departament de Biologia Animal, Biologia Vegetal i Ecologia, Universitat Autònoma de Barcelona, E-08193 Bellaterra, Spain 4 Global Change Research Institute, Czech Academy of Sciences, CZ-60300 Brno, Czech Republic; klem.k@czechglobe.cz (K.K.); urban.o@czechglobe.cz (O.U.) * Correspondence: guille.peguero@gmail.com † These authors contributed equally. Abstract: Tropical plants are expected to have a higher variety of defensive traits, such as a more diverse array of secondary metabolic compounds in response to greater pressures of antagonistic interactions, than their temperate counterparts. We test this hypothesis using advanced metabolomics   linked to a novel stoichiometric compound classification to analyze the complete foliar metabolomes of four tropical and four temperate tree species, which were selected so that each subset contained the Citation: Peguero, G.; same amount of phylogenetic diversity and evenness. We then built Bayesian phylogenetic multilevel Gargallo-Garriga, A.; Maspons, J.; models to test for tropical–temperate differences in metabolite diversity for the entire metabolome Klem, K.; Urban, O.; Sardans, J.; Peñuelas, J. Metabolome-Wide, and for four major families of secondary compounds. We found strong evidence supporting that Phylogenetically Controlled the leaves of tropical tree species have a higher phenolic diversity. The functionally closer group Comparison Indicates Higher of polyphenolics also showed moderate evidence of higher diversity in tropical species, but there Phenolic Diversity in Tropical Tree were no differences either for the entire metabolome or for the other major families of compounds Species. Plants 2021, 10, 554. https:// analyzed. This supports the interpretation that this tropical–temperate contrast must be related to doi.org/10.3390/plants10030554 the functional role of phenolics and polyphenolics. Academic Editors: Andrea Ghirardo Keywords: antagonistic interactions; bayesian phylogenetic models; latitudinal biodiversity gradient; and Alberto Minassi metabolomics; phenolics; plant defense Received: 12 February 2021 Accepted: 14 March 2021 Published: 16 March 2021 1. Introduction Publisher’s Note: MDPI stays neutral Attempting to explain the latitudinal gradient of plant diversity has been a major with regard to jurisdictional claims in focus of interest since the age of the pioneer naturalists [1]. Building upon these early published maps and institutional affil- observations, the modern synthesizers laid out the hypothesis that the strength of bi- iations. otic interactions, particularly antagonistic relationships with competitors, herbivores and pathogens, increase toward the tropics [2]. Accordingly, these stronger biotic interactions in the tropics must have had consequences in diversity generation accelerating diversification through coevolution and arms races [2–4]. However, they also must have consequences Copyright: © 2021 by the authors. in diversity maintenance increasing the species-carrying capacity of tropical ecosystems Licensee MDPI, Basel, Switzerland. through multiple mechanisms that facilitate species coexistence, such as a better niche This article is an open access article partitioning or increased niche packing and stronger negative density-dependent interac- distributed under the terms and tions [3–7]. However, all these arguments are sustained by a simple though sometimes conditions of the Creative Commons contentious prediction, i.e., tropical species must be better defended against antagonistic Attribution (CC BY) license (https:// agents because these are more diverse and exert a greater pressure [2]. In the case of creativecommons.org/licenses/by/ tropical plants, the wide array of secondary metabolic compounds is perhaps the most 4.0/). paradigmatic example of traits against herbivores and pathogens [8,9]. Plants 2021, 10, 554. https://doi.org/10.3390/plants10030554 https://www.mdpi.com/journal/plants Plants 2021, 10, 554 2 of 9 Many aspects of this overarching hypothesis have been well supported so far [10–17], but recent global studies and meta-analyses have questioned the evidence for some of its main predictions, for instance claiming that tropical plants are neither better defended against nor more damaged by herbivores [18,19]. Some suggestions have been posed to reconcile these conflicting views [20]. Accounting for the entire complexity of plant metabolomes, instead of focusing only on specific chemical compounds, could provide a broader perspective when conducting tropical–temperate comparisons [21]. Additionally, carefully considering the phylogenetic relationships between the species included in the analysis by properly applying comparative methods should allow to assess plant defenses within an evolutionary context [22]. Secondary metabolism in plants produces thousands of compounds in nearly 20 broad chemical families [23]. This extensive phytochemical diversity has been both the main driver of interest and the main hindrance to the comprehensive study of plant defense. The advances in mass spectrometry and nuclear magnetic resonance allow for unprece- dented characterizations of complete metabolomes, which are renewing interest in the links between taxonomic and phytochemical diversity [21,24–26]. In addition to these breakthroughs in analytical chemistry, the development of Bayesian comparative methods facilitates the modeling of the phylogenetic structure of the data and the inclusion of intraspecific variability in trait analyses, thereby seamlessly incorporating the evolutionary relationships among species within the inferential process [27]. Here, we test the hypothesis that tropical tree species present a higher diversity of secondary metabolites in their leaves using state-of-the-art liquid chromatography and mass spectrometry (LC–MS) linked to a novel, stoichiometrically based classification of compounds into the major chemical families [28]. We selected eight tree species from distinct tropical and temperate distributions while controlling for the phylogenetic diversity and evenness within each subset of species, and by means of Bayesian phylogenetic multilevel models we tested the prediction that the richness and diversity of metabolites should be greater for species at lower latitudes, particularly for the families of chemical compounds associated with protection against antagonistic biotic interactions, such as phenolics and polyphenolics. Additionally, we analyzed the similarity of the complete foliar metabolomes and of their phenolic profiles to check for differences according to region of origin and to assess variations within and between species in chemical composition. 2. Material and Methods 2.1. Species Selection and Sample Processing To exclude any evolutionary effect on the metabolomic differences between tropical and temperate species (e.g., contrasting species ages or patterns of diversification leading to a recent escalation or divergence of defenses [29]), we selected four tropical and four temperate tree species so that each subset encompassed the same amount of phylogenetic diversity and evenness, i.e., the total branch length and mean distances between the species in each phylogenetic subtree were equivalent (Table 1 and Figure S1). We collected leaves (2–10 per species from 1 to 10 individuals, n = 53 samples; see complete sample list in Table 2) directly from the canopies of mature trees (including sunlit and shaded leaves) in French Guiana and Spain, and immediately froze them in liquid nitrogen (N) directly in the field to stop all biochemical activity. The environmental conditions during the samplings were within the range of typical conditions for each one of the sites. The leaves were then lyophilized and ground with a ball mill in the laboratory. The metabolites of the homogenized samples were extracted using a 1:1 methanol:water solution, and the LC–MS raw chemical spectra were processed and compared using XCMS [30]. Peaks were identified in each sample and then matched across samples to calculate the deviations of retention time and to compare relative ion intensities. MS1 spectra were filtered, establish- ing a noise threshold at 2 × 104 and a minimum peak height at 2 × 105 . Chromatograms were aligned using an algorithm with a tolerance of 5 ppm of the mass-to-charge ratio (m/z) and a retention time of 0.2 min. Normalized areas under the peaks were used for Plants 2021, 10, 554 3 of 9 quantification and were log-transformed before statistical analysis (see Supplementary Methods in Supplementary S1 for further details). Table 1. Phylogenetic relationships of the tree species. Region Species Richness Phylodiversity Mean Pairwise Distance Mean Neighbor Taxon Distance Tropical 4 377 200 146 Temperate 4 392 208 157 Phylogenetic metrics (in millions of years) are derived from the plant phylogeny published by [31] and updated by [32]. Phylodiversity refers to the sum of the lengths of all branches on the phylogeny that span the members of the species set. Temperate species are Sambucus nigra (Adoxaceae), Salix alba (Salicaceae), Quercus ilex and Fagus sylvatica (Fagaceae). Tropical species are Goupia glabra (Goupiaceae), Protium decandrum (Burseraceae), Symphonia globulifera and Moronobea coccinea (Clusiaceae). See Figure S1 in Supplementary S1. Table 2. List of samples analyzed. Species Names and Family Sampled Leaves Individuals Fagus sylvatica L. (Fagaceae) 9 9 Goupia glabra Aubl. (Goupiaceae) 2 1 Moronobea coccinea Aubl. (Clusiaceae) 4 1 Protium decandrum Aubl. (Burseraceae) 4 1 Quercus ilex L. (Fagaceae) 10 10 Salix alba L. (Salicaceae) 10 4 Sambucus nigra L. (Adoxaceae) 10 10 Symphonia globulifera L.f. (Clusiaceae) 4 1 Temperate species were sampled from different sites in Catalonia (Spain), while tropical species were sampled from the Paracou and Nouragues research stations in French Guiana (France). We obtained a total of 2461 metabolic compounds, which were further classified based on the stoichiometry of N and hydrogen:carbon and oxygen:carbon ratios [28,33] into four broad families of secondary metabolites: highly unsaturated polyphenols (92 compounds), polycyclic aromatics (42 compounds), aliphatics (864 compounds) and phenolics (80 com- pounds). It must be noted that this classifying procedure is useful when metabolites have a clear stoichiometric signature that allows grouping them into broad groups. Important compounds known to be related with plant defense such as terpenoids are mostly included within aliphatics (e.g., those without oxygen atoms such as isoprene, and most mono-, sesqui-, di- and triterpenoids), while some other complex terpene compounds having aromaticity may be included in other categories, such as unsaturated polyphenols. For each sample, we calculated its metabolite richness (i.e., its number of different compounds) and its chemical diversity (CD) as: n CD = − ∑ pi × ln pi i =0 where pi is the peak area of a compound relative to the total peak area, and the sum includes all compounds detected (after [25]). 2.2. Data Analysis To test for differences in chemical richness and diversity for the entire metabolome and for each family of compounds, we built Bayesian multilevel models using the R package brms [34]. The richness and chemical diversity of the entire metabolome and for each family of metabolites calculated from each individual sample were modeled with region of origin (temperate vs. tropical) as a two-level, fixed-effects categorical factor, and the phylogenetic covariance matrix and the species identity as random factors, thus accounting for evolutionary relationships and the intraspecific variance [27,34,35]. We used the phylogenetic tree from Zanne et al. [31] to introduce the evolutionary relationships between the species in our models. We used a prior with a normal distribution with mean 0 and variance 10 times the standard deviation of the response variable to sample each Plants 2021, 10, 554 4 of 9 coefficient, and a prior with a half Student-t distribution with 3 degrees of freedom for each of the random factors. We ran the models with 4 chains of 2000 iterations, burning-in the first 1000. All models converged with a potential scale-reduction statistic Rˆ < 1.01. Hypothesis testing of the fixed-effect parameters was carried out by checking 95% credible intervals and computing their associated evidence ratios. An evidence ratio of X can be interpreted as the alternative hypothesis (e.g., that tropical species have a greater diversity of phenolic compounds) being X times more probable than the null hypothesis. This was calculated using the hypothesis function from the R package brms [34]. Finally, we analyzed the similarity in chemical composition by means of a partial least squares discriminant analysis (PLS-DA) with the mixOmics R package [36], and tested for differences between region of origin and species by means of permutational multivariate analysis of variance (PERMANOVA) using Bray–Curtis distance matrices with the R package vegan [37]. 3. Results The metabolomic profiles of the species indicated that tropical trees tended to have more metabolites than did temperate trees (2280 ± 148 vs. 2111 ± 147, respectively; hereafter mean ± standard deviation of the posterior distribution estimate). This trend toward a greater richness of metabolites was modest but consistent when considering each family of biochemical compounds separately (794 ± 56 vs. 730 ± 53 for aliphatics, 40 ± 2 vs. 38 ± 2 for polycyclic aromatics, 84 ± 7 vs. 75 ± 7 for polyphenolics and 73 ± 6 vs. 67 ± 6 for phenolics for tropical vs. temperate species, respectively). This tropical effect on the number of metabolites was nevertheless weak or inconclusive, mainly due to the particularly rich metabolic profiles of the two temperate Fagaceae species (see the marginal-density plots on the X-axes in Figure 1). In contrast to the number of metabolites, we found strong evidence that the diversity of phenolic compounds was higher in all four tropical species (the diversity increase in tropical species was 0.95 ± 0.44; (0.29–1.61) 95% CI, with an associated evidence ratio of 74.5; see the marginal-density plots on the Y-axis in Figure 1d). These compounds included procyanidins, flavone glycosides, chlorogenic acids and many other molecules with defensive properties against herbivores and pathogens [38]. Likewise, large unsatu- rated metabolites with many phenyl groups such as tannic acid, ellagitannins and other hydrolyzable tannins also tended to be more diverse in tropical species, although this difference was lower than with phenolic compounds (the diversity increase in tropical species for polyphenolics was 0.47 ± 0.49; (–0.29 to 1.23) 95% CI, with a moderate evidence ratio of 7.2; Figure 1c). In contrast to phenolics and polyphenolics, there were no differences in phytochemical diversity between tropical and temperate tree species either for the entire metabolome or for aliphatic and aromatic compounds (Figure 1a,b). The analysis of the complete foliar metabolomes by means of a PLS-DA showed that the tropical and temperate species clustered together, thus indicating a distinctive chemical composition of their leaves (Figure 2a). The clearly greater spread of the samples of tropical species suggests that for the whole foliar metabolome, both the interspecific and intraspecific variation seemed greater between and within the tropical species than in the temperate ones. An equivalent ordination using only the subset of phenolic compounds also showed that the tropical and temperate tree species had distinctive phenolic profiles, but interestingly, the tropical species apparently had less (or at least similar) intraspecific variation of their phenolic profiles than the temperate species (Figure 2b). The region of origin of the tree species was significant (p < 0.001) in both PERMANOVA tests, explaining 11% and 13% of the variation in the composition of the complete and of the phenolic com- ponent of the foliar metabolome, respectively. The differences between species, irrespective of their region of origin, were also significant (p < 0.001) and accounted for 49% and 53% of the variation in the composition of the complete and of the phenolic component of the foliar metabolome, respectively. Plants 2021, 10, 554 5 of 9 Plants 2021, 10, x FOR PEER REVIEW 5 of 10 (a) Aliphatics (b) Aromatics (c) Polyphenolics (d) phenolics Figure Figure 1. Relationshipsbetween 1. Relationships betweenthe the number number of of metabolites metabolites(metabolomic (metabolomic richness) with richness) metabolite with diversity metabolite (chemical- diversity (chemical- diversity index) for the main families of chemical compounds in the tropical and temperate tree species studied. Phenolic diversity index) for the main families of chemical compounds in the tropical and temperate tree species studied. Phenolic and polyphenolic compounds were more diverse in tropical species, with an evidence ratio of 74.5 and 7.2, respectively, andbased polyphenolic on Bayesiancompounds multilevel were more diverse models in tropical species, with an differences evidence ratio in theofmarginal-density 74.5 and 7.2, respectively, Plants 2021, 10, phylogenetic x FOR PEER REVIEW (note the tropical–temperate plots 6 of 10 based onY-axis of the Bayesian multilevel in panels c and phylogenetic d, and see textmodels (note for further the tropical–temperate differences in the marginal-density plots of details). the Y-axis in panels c and d, and see text for further details). The analysis of the complete foliar metabolomes by means of a PLS-DA showed that (a) tropical and temperate species clustered together, thus indicating a distinctive chemi- the cal composition of their leaves (Figure 2a). The clearly greater spread of the samples of tropical species suggests that for the whole foliar metabolome, both the interspecific and intraspecific variation seemed greater between and within the tropical species than in the temperate ones. An equivalent ordination using only the subset of phenolic compounds also showed that the tropical and temperate tree species had distinctive phenolic profiles, but interestingly, the tropical species apparently had less (or at least similar) intraspecific variation of their phenolic profiles than the temperate species (Figure 2b). The region of origin of the tree species was significant (p < 0.001) in both PERMANOVA tests, explaining 11% and 13% of the variation in the composition of the complete and of the phenolic com- ponent of the foliar metabolome, respectively. The differences between species, irrespec- tive of their region of origin, were also significant (p < 0.001) and accounted for 49% and 53% of the variation in the composition of the complete and of the phenolic component of the foliar metabolome, respectively. (b) Figure 2. Cont. Plants 2021, 10, 554 6 of 9 (b) FigureFigure 2. Similarity 2. Similarity in chemical in chemical composition composition for thefoliar for the complete complete foliar metabolome metabolome (a) and (b) (a) and for phenolics forbetween phenolicsthe tropical (b)and temperate between thetree speciesand tropical studied. The ordination temperate of individual tree species studied.samples resulted from The ordination of the partial least individual squares samples discriminant analyses (PLS-DAs) with species and region as group classifiers. The statistical significance and explained resulted from the partial least squares discriminant analyses (PLS-DAs) with species and region as variance of the insets resulted from a permutational multivariate analysis of variance using Bray–Curtis distance matrices. group classifiers. The statistical significance and explained variance of the insets resulted from a permutational multivariate analysis of variance using Bray–Curtis distance matrices. 4. Discussion This report provides strong evidence supporting that the diversity of phenolic com- 4. Discussion pounds, which are a key class of constitutive and induced plant defenses against patho- This report provides strong evidence gens and herbivores [9], was supporting higher for thethat the diversity tropical tree speciesofstudied, phenolic com- of irrespective pounds, which are a key class of constitutive and induced plant defenses against pathogens and herbivores [9], was higher for the tropical tree species studied, irrespective of their evo- lutionary history. The functionally closely related polyphenolics also provided moderate evidence supporting the same pattern, but importantly, this “tropical effect” was not de- tected for the entire metabolomic profile or for any other family of biochemical compounds assessed in this study. The fact that this tropical–temperate difference in phytochemical diversity is mostly due to the phenolic part of the foliar metabolome, and to a lower extent by the polyphenolic compounds, indicates that this metabolic contrast must be linked to the functional role that these compounds play. Since it is well established that the primary function of these compounds is protection against abiotic and biotic pressures, it seems reasonable to presume that tropical tree species must have been selected to keep a more diverse array of phenolic diversity than their temperate counterparts, likely due to stronger abiotic and biotic pressures [2]. In fact, numerous studies suggest that tropical trees must cope with stronger biotic interactions and a more diverse array of either generalist or specialized enemies [11,14,17,39,40]. This higher diversity in phenolic compounds in the leaves of tropical trees was neither the result of having more metabolites, which actually is hardly the case in the species studied here, nor the outcome of a divergent evolutionary history between the tropical and temperate subsets of species. The control for the phylogenetic relationships among the species included in the analysis (see Table 1) allowed us to exclude the emergence of the observed pattern only from differences in the timing and mode of diversification between the two latitudinal regions [7]. Additionally, the higher diversity of phenolic compounds in the leaves of the tropical species was coupled with a distinctive composition. Richards et al. [41] assessed the metabolomic profile, herbivore community and herbivory pressure of 22 Piper species from Costa Rica and found that phytochemical diversity re- duced herbivore damage, while it could also promote herbivore specialization. Similarly, Salazar et al. [38] reported that phenolic compounds were the most important chemical Plants 2021, 10, 554 7 of 9 defenses in 31 sympatric species within the monophyletic Protieae clade in Peru, and that their high diversity was due to cumulative diffuse interactions with a broad range of generalist herbivores. The only other latitudinal comparative study with such a broad metabolomics approach reported greater chemical divergence among tropical species, sug- gesting a stronger biotic selective force in the tropics and therefore more rapidly evolving defenses [42]. Our results are therefore in line with these previous studies, and also with the only assessment with a geographically and methodologically comparable approach [42]. The limited scope of this study, however, precludes a proper test of the mechanism behind the observed greater phytochemical diversity in phenolic and polyphenolic compounds. The evolution of plant defenses can escalate (i.e., the accumulation of chemical traits during the macroevolution of a plant lineage) and diverge (i.e., the increase in chemical dissimilarity between closely related plant species), both by generalist and specialized biotic antagonists, respectively [22,29]. Previous research has shown that that interspecific variation is greater than intraspecific variation among tropical tree species, even suggesting the existence of species-specific “metabolomic niches” [42–44]. This study, however, was not conceived to assess intraspecific variation in phytochemical diversity between tropical and temperate tree species. Thus, the limitations of our dataset in terms of number of species, samples per species, and particularly in the number of individuals per species, severely hinders the ability to draw firm conclusions on the potential role of intraspecific and interspecific variation on leaf chemical composition and the evolution of plant de- fenses. Even so, our finding that for the whole foliar metabolome, both the interspecific and intraspecific variation seems greater in the tropical species than in the temperate ones, agrees with what we would expect by their generally greater chemical richness. Moreover, when analyzing only the phenolic profiles, we found that tropical tree species apparently showed less (or similar) intraspecific variation than temperate ones, and this likely could be due to a lower capacity to exhibit inducible variation in defenses [45]. These suggestions should be properly investigated in further research efforts with more balanced sampling designs. It is also important to acknowledge that endophytic fungi may also synthesize defense compounds that cannot be easily distinguished from those produced by plants. Hence, part of the greater phytochemical diversity observed in the tropical tree species analyzed could also result from a greater contribution by mutualistic fungi. Neverthe- less, future research with a similar metabolomic approach applied to a phylogenetically broader dataset and including more sources of variation of intraspecific chemical com- position [43–45] should allow assessing phytochemical diversity across multiple plant lineages with contrasting patterns of diversification. This could help to disentangle the mechanisms driving the evolution of plant chemical defenses and the role that they may have in plant-species coexistence. Supplementary Materials: The following are available online at https://www.mdpi.com/2223 -7747/10/3/554/s1, Supplementary Materials ‘S1’: Extended Methods Figure S1: Phylogenetic relationships of the species. References [46,47] are cited in the supplementary materials. Author Contributions: G.P., A.G.-G., J.M. and J.P. designed the study; A.G.-G., K.K. and O.U. compiled the metabolomic data sets; G.P., A.G.-G. and J.M. analyzed the data; G.P. led the writing; and A.G.-G., J.M., K.K., O.U., J.S. and J.P. contributed ideas and provided substantial inputs. All authors have read and agreed to the published version of the manuscript. Funding: This research was supported by the European Research Council Synergy grant ERC-2013- SyG 610028-IMBALANCE-P, the Spanish Government grant PID2019-110521GB-I00 and the Catalan Government grant SGR 2017-1005. Institutional Review Board Statement: Not applicable. Informed Consent Statement: Not applicable. Data Availability Statement: Raw data supporting the results of this study are available at Figshare repository https://doi.org/10.6084/m9.figshare.12240581.v1. Last accessed date 15 March 2021. Conflicts of Interest: The authors no conflict of interest. Plants 2021, 10, 554 8 of 9 References 1. Brown, J.H. Why are there so many species in the tropics? J. Biogeogr. 2014, 41, 8–22. [CrossRef] [PubMed] 2. Schemske, D.W.; Mittelbach, G.G.; Cornell, H.V.; Sobel, J.M.; Roy, K. Is there a latitudinal gradient in the importance of biotic interactions? Annu. Rev. Ecol. Evol. Syst. 2009, 40, 245–269. [CrossRef] 3. Dobzhansky, T. Evolution in the tropics. Am. Sci. 1950, 38, 209–221. 4. Ehrlich, P.R.; Raven, P.H. Butterflies and Plants: A Study in Coevolution. Evolution 1964, 18, 586–608. [CrossRef] 5. Janzen, D.H. Herbivores and the number of tree species in tropical forests. Am. Nat. 1970, 104, 501–528. [CrossRef] 6. MacArthur, R.H. Geographical Ecology: Patterns in The Distribution of Species; Princeton University Press: Princeton, NJ, USA, 1972. 7. Mittelbach, G.G.; Schemske, D.W.; Cornell, H.V.; Allen, A.P.; Brown, J.M.; Bush, M.B.; Harrison, S.P.; Hurlbert, A.H.; Knowlton, N.; Lessios, H.A.; et al. Evolution and the latitudinal diversity gradient: Speciation, extinction and biogeography. Ecol. Lett. 2007, 10, 315–331. [CrossRef] 8. Coley, P.D.; Barone, J.A. Herbivory and Plant Defenses in Tropical Forests. Annu. Rev. Ecol. Syst. 1996, 27, 305–335. [CrossRef] 9. Levin, D.A. The Chemical Defenses of Plants to Pathogens and Herbivores. Annu. Rev. Ecol. Syst. 1976, 7, 121–159. [CrossRef] 10. Bagchi, R.; Gallery, R.E.; Gripenberg, S.; Gurr, S.J.; Narayan, L.; Addis, C.E.; Freckleton, R.P.; Lewis, O.T. Pathogens and insect herbivores drive rainforest plant diversity and composition. Nature 2014, 506, 85–88. [CrossRef] [PubMed] 11. Salazar, D.; Marquis, R.J. Herbivore pressure increases toward the equator. Proc. Natl. Acad. Sci. USA 2012, 109, 12616–12620. [CrossRef] [PubMed] 12. Kursar, T.A.; Dexter, K.G.; Lokvam, J.; Pennington, R.T.; Richardson, J.E.; Weber, M.G.; Murakami, E.T.; Drake, C.; McGregor, R.; Coley, P.D. The evolution of antiherbivore defenses and their contribution to species coexistence in the tropical tree genus Inga. Proc. Natl. Acad. Sci. USA 2009, 106, 18073–18078. [CrossRef] [PubMed] 13. Becerra, J.X. On the factors that promote the diversity of herbivorous insects and plants in tropical forests. Proc. Natl. Acad. Sci. USA 2015, 112, 6098–6103. [CrossRef] 14. Lim, J.Y.; Fine, P.V.A.; Mittelbach, G.G. Assessing the latitudinal gradient in herbivory. Glob. Ecol. Biogeogr. 2015, 24, 1106–1112. [CrossRef] 15. Roslin, T.; Hardwick, B.; Novotny, V.; Petry, W.K.; Andrew, N.R.; Asmus, A.; Barrio, I.C.; Basset, Y.; Boesing, A.L.; Bonebrake, T.C.; et al. Higher predation risk for insect prey at low latitudes and elevations. Science 2017, 356, 742–744. [CrossRef] 16. Fine, P.V.A.; Mesones, I.; Coley, P.D. Herbivores Promote Habitat Specialization by Trees in Amazonian Forests. Science 2004, 305, 663–665. [CrossRef] 17. Mangan, S.A.; Schnitzer, S.A.; Herre, E.A.; Mac, K.K.M.L.; Valencia, M.C.; Sanchez, E.I.; Bever, J.D. Negative plant-soil feedback predicts tree-species relative abundance in a tropical forest. Nature 2010, 466, 752–755. [CrossRef] 18. Moles, A.T.; Bonser, S.P.; Poore, A.G.B.; Wallis, I.R.; Foley, W.J. Assessing the evidence for latitudinal gradients in plant defence and herbivory. Funct. Ecol. 2011, 25, 380–388. [CrossRef] 19. Moles, A.T.; Wallis, I.R.; Foley, W.J.; Warton, D.I.; Stegen, J.C.; Bisigato, A.J.; Pizarro, L.C.; Clark, C.J.; Cohen, P.S.; Cornwell, W.K.; et al. Putting plant resistance traits on the map: A test of the idea that plants are better defended at lower latitudes. New Phytol. 2011, 191, 777–788. [CrossRef] 20. Anstett, D.N.; Nunes, K.A.; Baskett, C.; Kotanen, P.M. Sources of Controversy Surrounding Latitudinal Patterns in Herbivory and Defense. Trends Ecol. Evol. 2016, 31, 789–802. [CrossRef] 21. Sedio, B.E. Recent breakthroughs in metabolomics promise to reveal the cryptic chemical traits that mediate plant community composition, character evolution and lineage diversification. New Phytol. 2017, 214, 952–958. [CrossRef] [PubMed] 22. Agrawal, A.A. Macroevolution of plant defense strategies. Trends Ecol. Evol. 2006, 22, 103–109. [CrossRef] [PubMed] 23. Wink, M. Evolution of secondary metabolites from an ecological and molecular phylogenetic perspective. Phytochemistry 2003, 64, 3–19. [CrossRef] 24. Hall, R.D. Plant metabolomics: From holistic hope, to hype, to hot topic. New Phytol. 2005, 169, 468–543. [CrossRef] 25. Hilker, M. New Synthesis: Parallels Between Biodiversity and Chemodiversity. J. Chem. Ecol. 2014, 40, 225–226. [CrossRef] 26. Dyer, L.A.; Philbin, C.S.; Ochsenrider, K.M.; Richards, L.A.; Massad, T.J.; Smilanich, A.M.; Forister, M.L.; Parchman, T.L.; Galland, L.M.; Hurtado, P.J.; et al. Modern approaches to study plant–insect interactions in chemical ecology. Nat. Rev. Chem. 2018, 2, 50–64. [CrossRef] 27. Hadfield, J.D.; Nakagawa, S. General quantitative genetic methods for comparative biology: Phylogenies, taxonomies and multi-trait models for continuous and categorical characters. J. Evol. Biol. 2010, 23, 494–508. [CrossRef] 28. Rivas-Ubach, A.; Liu, Y.; Bianchi, T.S.; Tolić, N.; Jansson, C.; Paša-Tolić, L. Moving beyond the van Krevelen Diagram: A New Stoichiometric Approach for Compound Classification in Organisms. Anal. Chem. 2018, 90, 6152–6160. [CrossRef] [PubMed] 29. Volf, M.; Segar, S.T.; Miller, S.E.; Isua, B.; Sisol, M.; Aubona, G.; Šimek, P.; Moos, M.; Laitila, J.; Kim, J.; et al. Community structure of insect herbivores is driven by conservatism, escalation and divergence of defensive traits in Ficus. Ecol. Lett. 2018, 21, 83–92. [CrossRef] 30. Smith, C.A.; Want, E.J.; O’Maille, G.; Abagyan, R.; Siuzdak, G. XCMS: Processing mass spectrometry data for metabolite profiling using nonlinear peak alignment, matching, and identification. Anal. Chem. 2006, 78, 779–787. [CrossRef] 31. Zanne, A.E.; Tank, D.C.; Cornwell, W.K.; Eastman, J.M.; Smith, S.A.; FitzJohn, R.G.; McGlinn, D.J.; O’Meara, B.C.; Moles, A.T.; Reich, P.B.; et al. Three keys to the radiation of angiosperms into freezing environments. Nature 2014, 506, 89–92. [CrossRef] Plants 2021, 10, 554 9 of 9 32. Kellerman, A.M.; Dittmar, T.; Kothawala, D.N.; Tranvik, L.J. Chemodiversity of dissolved organic matter in lakes driven by climate and hydrology. Nat. Commun. 2014, 5, 3804. [CrossRef] 33. Bürkner, P.C. Brms: An R package for Bayesian multilevel models using Stan. J. Stat. Softw. 2017, 80. [CrossRef] 34. Gelman, A.; Hill, J. Data Analysis Using Regression and Multilevel/Hierarchical Models; Cambridge University Press: Cambridge, MA, USA, 2006. 35. Qian, H.; Jin, Y. An updated megaphylogeny of plants, a tool for generating plant phylogenies and an analysis of phylogenetic community structure. J. Plant Ecol. 2016, 9, 233–239. [CrossRef] 36. Rohart, F.; Gautier, B.; Singh, A.; Cao, K.A.L. MixOmics: An R package for ‘omics feature selection and multiple data integration. PLoS Comput. Biol. 2017, 13, e1005752. [CrossRef] [PubMed] 37. Oksanen, J.; Blanchet, F.G.; Kindt, R.; Friendly, M.; Legendre, P.; McGlinn, D.; Minchin, P.R.; O’Hara, R.B.; Simpson, G.L.; Solymos, P.; et al. Vegan: Community Ecology Package. R Package Version 2.4-3. Vienna: R Foundation for Statistical Computing. 2016. Available online: https://www.researchgate.net/publication/323265822_vegan_Community_Ecology_Package_R_package_ version_24-3_2017_accessed_2016_Jan_1 (accessed on 13 March 2021). 38. Salazar, D.; Lokvam, J.; Mesones, I.; Pilco, M.V.; Zuñiga, J.M.A.; de Valpine, P.; Fine, P.V.A. Origin and maintenance of chemical diversity in a species-rich tropical tree lineage. Nat. Ecol. Evol. 2018, 2, 983–990. [CrossRef] 39. Forister, M.L.; Novotny, V.; Panorska, A.K.; Baje, L.; Basset, Y.; Butterill, P.T.; Cizek, L.; Coley, P.D.; Dem, F.; Diniz, I.R.; et al. The global distribution of diet breadth in insect herbivores. Proc. Natl. Acad. Sci. USA 2015, 112, 442–447. [CrossRef] [PubMed] 40. Peguero, G.; Bonal, R.; Sol, D.; Muñoz, A.; Sork, V.L.; Espelta, J.M. Tropical insect diversity: Evidence of greater host specialization in seed-feeding weevils. Ecology 2017, 98, 2180–2190. [CrossRef] [PubMed] 41. Richards, L.A.; Dyer, L.A.; Forister, M.L.; Smilanich, A.M.; Dodson, C.D.; Leonard, M.D.; Jeffrey, C.S. Phytochemical diversity drives plant–insect community diversity. Proc. Natl. Acad. Sci. USA 2015, 112, 10973–10978. [CrossRef] 42. Sedio, B.E.; Parker, J.D.; McMahon, S.M.; Wright, S.J. Comparative foliar metabolomics of a tropical and a temperate forest community. Ecology 2018, 99, 2647–2653. [CrossRef] 43. Sedio, B.E.; Echeverri, J.C.R.; Boya, C.A.; Wright, S.J. Sources of variation in foliar secondary chemistry in a tropical forest tree community. Ecology 2017, 98, 616–623. [CrossRef] 44. Garriga, A.G.; Sardans, J.; Granda, V.; Llusià, J.; Peguero, G.; Asensio, D.; Ogaya, R.; Urbina, I.; van Langenhove, L.; Verryckt, L.T.; et al. Different “metabolomic niches” of the highly diverse tree species of the French Guiana rainforests. Sci. Rep. 2020, 10, 6937. [CrossRef] [PubMed] 45. Sedio, B.E.; Archibold, A.D.; Echeverri, J.C.R.; Debyser, C.; Cristopher, A.B.P.; Wright, S.J. A comparison of inducible, ontogenetic, and interspecific sources of variation in the foliar metabolome in tropical trees. PeerJ 2019, 7, e7536. [CrossRef] [PubMed] 46. Rivas-Ubach, A.; Pérez-Trujillo, M.; Sardans, J.; Gargallo-Garriga, A.; Parella, T.; Peñuelas, J. Ecometabolomics: Optimized NMR-based method. Methods Ecol. Evol. 2013, 4, 464–473. [CrossRef] 47. Koch, B.P.; Dittmar, T.; Witt, M.; Kattner, G. Fundamentals of molecular formula assignment to ultrahigh resolution mass data of natural organic matter. Anal. Chem. 2007, 79, 1758–1763. [CrossRef]