Abstract
Glandular trichomes are well known to participate in plant chemical and physical defenses against herbivores, especially herbivorous insects. However, little is known about large-scale geographical patterns in glandular trichome occurrence. Herbivory pressure is thought to be higher at low elevations because of warmer and more stable climates. We therefore predicted a higher proportion of species with glandular trichomes at low elevations than at higher elevations. We compiled glandular trichome data (presence/absence) for 6,262 angiosperm species from the Hengduan Mountains (a global biodiversity hotspot in southwest China). We tested the elevational gradient (800–5,000 m a.s.l.) in the occurrence of plant species with glandular trichomes, and its correlations with biotic (occurrence of herbivorous insects) and abiotic factors, potentially shaping the elevational gradient in the occurrence of glandular trichomes. We found a significantly positive relationship between elevation and the occurrence of glandular trichomes, with the proportion of species having glandular trichomes increasing from 11.89% at 800 m a.s.l. to 17.92% at above 4,700 m. This cross-species relationship remained significant after accounting for phylogenetic relationships between species. Herbivorous insect richness peaked at mid-elevations and its association with the incidence of glandular trichomes was weak. Mean annual temperature was the most important factor associated negatively with glandular trichomes. Our results do not support the hypothesis that plant defenses decrease with increasing elevation. In contrast, a higher proportion of plant species with glandular trichome toward higher elevations is observed. Our results also highlight the importance of considering the simultaneous influences of biotic and abiotic factors in testing geographical variation in multifunctional plant defenses.
Introduction
Geographical variation in species interactions (e.g., herbivory and predation) is widespread and is thought to generate concomitant patterns of species abundance, diversity, and functional traits (Thompson et al., 2005; ; ). Thus, studying the geographical patterns of biotic interactions is especially useful for our understanding of spatial variation in community composition within different environments (Schoonhoven et al., 2005). It has been widely believed that herbivory, one of the most common and important biotic (plant-animal) interactions, is more intense under warmer and more stable climates, found at lower elevations (Reynolds and Crossley, 1997; Metcalfe et al., 2014). In response to higher herbivory pressure, plants at lower elevations are predicted to invest more in their defenses (elevational herbivory defense hypothesis; Pellissier et al., 2014; Rasmann et al., 2014). Many studies have tested this hypothesis and have found a negative (), positive (e.g., Moreira et al., 2018), or no relationship between elevation and defense (). Thus, despite the widespread interest in such a pattern, there is a disagreement about the direction of a general elevational gradient in plant defense.
There are several possible and reasonable reasons why the available evidence for elevational patterns in plant defense is mixed and unconvincing. First, because different defense-related traits may be effective against different herbivores, for example, spines mainly function in deterring large mammalian herbivores, while trichomes mostly deter insect herbivores (). Different herbivore guilds with different life-histories usually respond in different ways to biotic or abiotic conditions and thus pooling various defense traits may gain contrasting results of elevational gradients in defense (Rasmann et al., 2014). Second, findings from several studies suggested that the tested elevational span was insufficient, and that studying only a relatively small number of species may also result in completely opposite results (e.g., Moles et al., 2011; ). Third, if the data used in analyses originated from sampling with very different (and often inappropriate) methods, geographic pattern of a specific defense may be diluted (Moles et al., 2011), for example, variations in the extraction methodologies are suggested to influence the quantification of condensed tannins in leaves (). Fourth, plants with different growth forms might evolve different growth- and lifestyle-related defense strategies, ultimately influencing plant defense investment (Ronel and Lev-Yadun, 2012; ; Xu et al., 2020). Unfortunately, variation in defenses among plant growth forms has frequently been overlooked or comparison of defenses among plant growth forms were limited by small datasets. Finally, it is well recognized that various defense traits respond not only to herbivore pressure, but could also be associated with abiotic conditions (Pratt et al., 2014) and resource availability (). For example, concentration of phenolics in leaves of Quercus robur increased independently with decreasing temperature toward higher elevations (). Thus, studies not considering the elevational variation in abiotic conditions may obscure the elevational gradients in plant defenses.
Indeed, overcoming the deficiencies mentioned above was very difficult in the past. For example, sampling including large number of plant species and spanning sufficient geographic breadth would be an enormously time- and effort-consuming enterprise. Now, these problems can be solved by compiling data from the literature and from data-bases. In this study, we used a consistent methodology to compare investment of a single but common and easily identified defense trait (i.e., glandular trichomes, a trait detailed described in the Flora of China) across a large number of plant species along a wide elevational range by collecting published literatures.
Glandular trichomes, an important and common trichome type, often produce, store, and secrete various secondary metabolites such as flavonoids, monoterpenes, or sesquiterpene lactones, and occur widely on many plant organs (Wagner, 1991). Although glandular trichomes may reflect an adaptation to various abiotic factors such as high-level O3 exposure, high levels of metal and non-metal ions in the soil, high vapor pressure deficit (Silva et al., 1996; Lihavainen et al., 2017; ), one of their most important well-known roles is to serve as a defense against herbivores, especially insects (; Wagner, 1991; ). Insect herbivores can be killed directly by ingested toxic exudates produced by glandular trichomes (; Liu et al., 2019). Some glandular trichomes can entrap insects arriving on plant surface by sticky exudates, and the entrapped insects usually die as a result of starvation or, in the case of small insect herbivores, of suffocation (Puterka et al., 2003; Simmons et al., 2004). Thus, glandular trichomes have attracted much attention from evolutionary ecologists who focus on how plants deter insect herbivores (e.g., ; ; ; ; Liu et al., 2019). For example, both glandular trichome density and primary volatile organic compounds produced by glandular trichomes on leaf surfaces of Ocimum gratissimum plants increased significantly after being attacked by leaf-cutter ants (Tozin et al., 2017). Although we have numerous studies quantifying the geographic variation in glandular trichomes for single species (e.g., ; Zhao et al., 2019), few studies have been conducted to test the roles of glandular trichomes in a large-scale and geographically diverse territory. These biogeographic patterns can advance our understanding of plant-herbivore interactions.
In this study, we tested the elevational variation in occurrence of glandular trichome in a large number of plant species and tested the association of the occurrence of glandular trichomes with both abiotic and biotic factors by compiling glandular trichome data (presence/absence) for 6,262 angiosperm species from the Hengduan Mountain region in southwestern China, ranging from 800 to 5,000 m a.s.l. Specifically, we asked: (1) is there an overall elevational gradient in the occurrence of plants with glandular trichomes, and does this gradient vary across growth forms? and (2) is there an association of abiotic and biotic factors with any such elevational variation in glandular trichomes?
Materials and Methods
Study Area
This study covers the geographical region of the Hengduan Mountains in southwestern China, one of the 36 most important global biodiversity hotspots (). The study area includes northwestern Yunnan, western Sichuan, southeastern Tibet, southeastern Qinghai, and southern Gansu, stretching between 24°39′ to 34°20′ N, and 89°27′ to 104°36′ E, covering about 500,000 km2 (Figure 1). The average elevation drops from 4,000 to 5,000 m in the Western Sichuan and Eastern Tibet to 2,000 m in the Northwestern Yunnan. The climate in the Hengduan Mountains region is mainly influenced by the southeast monsoon and the plateau monsoon, with clear, dry winters, and warm, humid summers (Li and Zhang, 2010). It harbors diverse vegetation types, including dry valley scrub, broad-leaved forests, coniferous forests, alpine shrublands, meadows, and scree (Zhang et al., 2002).
FIGURE 1
Data Collection
We compiled a list of all angiosperm species present in the study area, from the Vascular Flora of Hengduan Mountains (Wang, 1993), the Flora of Yunnan (Wu, 1977–2006), the Flora of Sichuan (), the Flora of Qinghai (Liu, 1996), and the Flora of Tibet (Wu, 1983-1987). In order to determine the elevational range of each species included in our dataset, we downloaded all specimen information of the studied species collected from Hengduan Mountains region from the National Plant Specimen Resource Center1. Overall, our dataset included 147,150 specimens with elevation distribution ranging from 500 to 6,300 m. The elevation gradient between 500 and 6,300 m was divided into 100 m elevation intervals. Each species was considered as present at each 100 m belt between its upper and lower elevation limits (). This is a commonly used method (e.g., Vetaas and Grytnes, 2002; ; Zhang et al., 2009; ), but it may result in an under-estimation of species richness at the extreme elevations or an over-estimation at mid-elevations (Vetaas and Grytnes, 2002). However, there is no reason to expect that any under- or over-estimation of species’ richness is likely to be biased in a way that would influence our estimation of the elevational gradient in the proportion of species with glandular trichomes. In order to determine whether elevational patterns of glandular trichomes was influenced by division of elevational intervals, the elevational gradients were also divided into 300 m and 500 m elevational intervals (also for herbivorous insect richness, mean annual temperature, and calculated plant-available water). In total, our dataset included 6,262 angiosperm species from 1,094 genera belonging to 175 plant families. We classified the species into one out of two growth form categories: herbaceous and woody plants (). Of the 6,262 species in our study, 4,010 were herbaceous (64.04%) and 2,252 were woody (35.96%). The proportion of herbaceous plants species increases with elevation, while the proportion of woody plants species decreases with elevation (Supplementary Figure 1).
Given that species with glandular trichomes must have trichomes, we began by scoring each species as having or not having any type of trichomes on any aboveground part of the plant over the entire life history according to the species’ descriptions in Flora of China (Wu et al., 1994–2012). Species with descriptions including the following terms were defined as having trichomes: pubescent, vill, bristle, puberul, toment, hair, pilose, sericeous, hirsute, setose, ciliate, strigose, velutinous, trichomes, lanate, hispid, wool, setulose, setae, fluff, and tomentellous. Species lacking any of these terms were sorted as not having trichomes. Then, we determined whether the trichomes of these species are glandular or not by checking the species descriptions in the Flora of China (Wu et al., 1994–2012) and specimens or communicating with experts. Although species may have intraspecific variation in density or phenotype of glandular trichomes within its range or between different development stages (), it is impossible to collect these data at our study scale. Thus, following the method of Moles et al. (2020), if a species was scored as having glandular trichomes, it was considered as having glandular trichomes at all elevations between its upper and lower limits and at various developmental stages.
It is impossible to gather species-specific data on insect herbivory at the huge scale of our analysis; thus, we used herbivorous insect richness as a proxy for herbivore pressure, due to their high correlation (; ). We collected data on the elevational range for herbivorous insects in this area from Insects of the Hengduan Mountains Region () supplemented by personal communications with experts. A total of 2,979 insect species with known elevational distribution were included in the dataset, including beetles, moths, cicadas, and other groups of herbivorous insects. Although some insects may change their feeding habits across developmental stages, for example, larvae of some butterflies may be phytophagous, whereas their adults may be important pollinators, we considered them as herbivorous as long as they are phytophagous at any stages. Elevational ranges for each insect species were also interpolated using the same methods described above and divided into 100 m, 300 m, and 500 m elevational belts.
Mean annual temperature and mean annual precipitation are consistently considered to be correlated with species richness (; Zhang et al., 2016). Thus, we download data of mean annual temperature (MAT; °C) from the WorldClim version2 database (). Considering that the available water for plants (PAW; cm3/cm3), defined as the soil water between field capacity and permanent wilting point, is superior to precipitation in reflecting the water actually available for plants (), we collected data on PAW from the global maps of soil hydraulic parameters (Zhang et al., 2018). Environmental layers at a 30-arcsecond resolution (∼1 km2) were projected into the same Albers Equal area coordinate system used for the study area. The longitude and latitude of each grid was obtained from the variable layer, and the corresponding elevation of each grid was obtained from the elevation layer using the longitude and latitude. The elevation and variable value of each grid were calculated, and the mean value of the data within each elevational belt of 100 m, 300 m, and 500 m was calculated to represent the variable value of this elevation belt.
Data Analyses
In all analyses, the sampling unit was the each divided elevational belt. Because there was no adequate data about herbivorous insects for the elevational belts below 800 m and above 5,000 m, in combination with low plant species numbers within these elevational belts, only the elevational belts between 800 and 5,000 m were included in our analyses.
We used simple logistic regression to test the association between glandular trichome and elevation using “glm” function following and Pincheira-Donoso et al. (2017). For logistic regression analyses, glandular trichome data were modeled as a binary response variable (0 = absence of glandular trichome, 1 = presence of glandular trichome). That is, in all logistic regression models, the response variable was the number of species with glandular trichomes relative to the number of species without glandular trichomes in each elevational belt (; ). Using the “cbind” function to combine the number of species with and without glandular trichomes in the response metric, our models weighted data points according to the number of species contained in elevational belt (). Subsequently, to test whether the elevational gradient in incidence of species with glandular trichomes is related to growth form, multiple logistic regression was conducted where glandular trichome occurrence was predicted by elevation, growth form (herbaceous or woody plants) and their interaction.
Linear polynomial regression models were performed to test the elevational patterns of herbivorous insect richness, MAT and PAW as a function of elevation along the gradient. The best-fit models were selected by comparing models that included elevation as linear term or models that included elevation as both linear and quadratic terms based on the Akaike Information Criterion (AIC) value using “AIC” function (Supplementary Table 1). We performed simple logistic regression of incidence of glandular trichomes against each of the three factors (insect richness, MAT, and PAW) in order to examine the potential mechanisms of individual factors in explaining elevational gradients of glandular trichomes. Since no collinearity was found between these factors (Supplementary Table 2 and Supplementary Figure 2), we conducted multiple logistic regression analyses in order to explore the multivariate explanations for the elevational pattern of glandular trichomes, where the incidence of glandular trichomes was predicted by all biotic and abiotic variables ().
In order to determine the relationship between glandular trichomes and elevation, herbivorous insect richness, mean annual temperature, and plant-available water in the context of phylogeny, phylogenetic logistic regression analyses were conducted with “phylolm” package (). Phylogenetic logistic regression can be used to estimate parameter alpha and calculate likelihood of trait data. The parameter alpha, which reflects the overall rate of transition between binary states, which when close to zero corresponds to stronger phylogenetic effects (). For these analyses, we generated a phylogenetic tree of study species using “V.PhyloMaker” package (), which was implemented with a mega-tree derived primarily from Smith and Brown’s (2018) phylogeny. In addition, the elevation and the corresponding herbivorous insect richness, mean annual temperature, and plant-available water for each occurrence were extracted and the species’ mean values of all these variables were used (; ).
Prior to the analysis, predictor variables were normalized using “scale” function to facilitate comparison of regression coefficients. For generalized linear models, the goodness-of-fit of each model was assessed by calculating McFadden’s R2 equivalent to R2 in ordinary least squares models (McFadden, 1974). All statistical analyses were conducted in the R (version 3.6.3) (R Development Core Team, 2020).
Results
The overall proportion of species with glandular trichomes in the Hengduan Mountains region was 14.1% (886 out of 6,262). The probability of finding species with glandular trichomes significantly increased toward higher elevations (Figure 2), ranging from 11.89% at 800 m a.s.l. to 17.92% at above 4,700 m. Multiple regression including elevation, growth form and their interaction indicated a significant effect of growth form on glandular trichomes (overall proportion of species with glandular trichomes: 15.6% for herbaceous species versus 11.6% for woody species). However, there was no significant interaction between growth form and elevation (Table 1), indicating that the relationship between glandular trichomes and elevation was not significantly different between species with different growth forms (Supplementary Figure 3).
FIGURE 2
TABLE 1
| Variable | z | p-value |
| Elevation | 9.52 | <0.01 |
| Growth form | –4.60 | <0.01 |
| Elevation × growth form | 1.236 | 0.22 |
Logistic regression of elevation, growth form (herbaceous, woody), and their interaction against glandular trichome in the Hengduan Mountain region.
Elevational gradient in herbivorous insects’ species richness was unimodal, with the highest richness occurring at mid-elevations (between 2,000 and 3,000 m) (Figure 3A). Thus, insect herbivores did not show higher species richness at lower elevations. A significant association was found between mean annual temperature and elevation, with the mean annual temperature decreasing with increasing elevation (Figure 3B). Similar to herbivorous insect richness, plant-available water also showed a unimodal association with elevation, i.e., plants at mid-elevations having the highest water availability (Figure 3C). All factors, including biotic and abiotic variables, showed significant association with the incidence of plants with glandular trichomes when each predictor was considered alone (Figure 4). The strongest correlation was with annual mean temperature (R2 = 0.40), followed by the plant-available water (R2 = 0.31). Herbivorous insect richness showed a weak correction with the probability of finding plant species with glandular trichomes (R2 = 0.08). We also conducted analyses at 300 m or 500 m elevational intervals, and found that the main findings did not change (Supplementary Figures 4, 5).
FIGURE 3
FIGURE 4
When herbivorous insect richness, mean annual temperature, and plant-available water were considered together, the mean annual temperature was negatively associated, but plant-available water and herbivorous insect richness (marginally) were positively associated with the incidence of plant species with glandular trichomes (Table 2). The multiple logistic regression model explained 43.25% of the variation in the occurrence of plant species with glandular trichomes.
TABLE 2
| Predictors | β | z | p-value |
| HIR | 0.03 | 1.98 | 0.048 |
| MAT | –0.17 | –6.71 | <0.001 |
| PAW | 0.04 | 2.27 | 0.02 |
Multiple regression models testing for the effects of biotic (HIR, herbivorous insect richness) and abiotic (MAT, mean annual temperature; PAW, plant-available water) on incidence of glandular trichomes in the Hengduan Mountains region.
β values were calculated after predictors being scaled.
The incidence of species with glandular trichomes varied among families (Supplementary Figure 6), and the presence/absence of glandular trichomes was more phylogenetically conserved than being random (alpha = 0.02; ). Phylogenetic logistic regression analyses showed that the relationships between glandular trichomes and elevation and plant-available water remained significantly positive (P < 0.001; Supplementary Figure 7), and the relationships between glandular trichomes and herbivorous insect richness and mean annual temperature remained significantly negative (P < 0.001; Supplementary Figure 7). These results indicated that our main findings are unlikely to change when the evolutionary relationships among plant species were considered.
Discussion
The question of how interactions between plants and herbivores varied along environmental gradients has interested ecologists for decades because of its importance for understanding mechanisms driving spatial variation in biodiversity and functional traits (). A long-held view in ecology is that plants from low elevations are better defended than are plants from high elevations (Schemske et al., 2009). Contrary to that prediction, our results showed that in the Hengduan Mountains region, plant species at higher elevations are more likely to have glandular trichomes even after accounting for the type of growth form and the evolutionary relationships between plant species, which is consistent with some previous studies (reviewed by Moreira et al., 2018). The congruence of these independent studies focusing on different defense traits strongly suggests that this widely accepted view needs to be reconsidered.
Defense investments are considered to reflect the historical herbivore pressure experienced by the plants, which is predicted to increase with decreasing elevations, due to the more benign and stable climate (Pellissier et al., 2012; Raffa et al., 2013). This served as the foundation for the view that low-elevation plants are thought to be better defended than are plants from higher elevations (Rasmann et al., 2014). Although we do not have direct data on herbivory pressures from insect herbivores experienced by the plants in our study area, many studies have shown that herbivorous insect richness could be used as a measure of herbivory pressure, since these two things are highly correlated (; ). In accordance with many other studies of herbivore pressure along elevational gradients (reviewed by Moreira et al., 2018), in the Hengduan Mountains region, herbivorous insect richness was not the highest at low elevations, but rather peaked at mid elevations, which means that there is no higher selection pressure for higher defense against herbivorous insects at lower elevations. Thus, it is not surprising that plants at lower elevations did not have the highest probability of having glandular trichomes.
Interestingly, although herbivorous insect richness varied along elevational gradients (i.e., peak at mid-elevation), elevational variation in herbivory pressure did not result in concomitant gradients in the proportion of plants with glandular trichomes. Furthermore, our univariate regression showed a negative, albeit weak (R2 = 0.08), relationship between herbivorous insect richness and incidence of species with glandular trichomes along the elevational gradients. In addition to herbivorous insects, glandular trichomes have also been found to function in protecting against mollusks (e.g., Westerbergh and Nyberg, 1995; Tomas et al., 2019). Unfortunately, we did not have herbivory or distribution data for mollusks in our study areas. Studies conducted at other sites have found that gastropod diversity generally decreased with increasing elevation (e.g., Liew et al., 2010; ; Schmera and Baur, 2014). If such pattern is true for the Hengduan Mountains region, higher proportion of species with glandular trichomes at higher elevations is unlikely to relate to herbivory pressure from mollusks in our study area. Our results indicated that considering only herbivory pressure by herbivores is insufficient to explain geographic variation in plant defenses in general (e.g., ; Tindall et al., 2016; Moles et al., 2020), or the occurrence of glandular trichomes in particular.
Recent work has emphasized the importance of considering abiotic factors in testing for geographic variations in interactions between plants and herbivores (; ; Zhang et al., 2016). In this study, we analyzed the effects of multiple abiotic factors associated with elevational gradients in plant traits and herbivorous insect richness on the incidence of glandular trichomes, and found that abiotic factors significantly better explained the variation (e.g., 40% and 31% for mean annual temperature and plant-available water, respectively) than herbivorous insect richness (8%). The proportion of plant species with glandular trichomes increases with decreasing temperature and with increasing water availability. Furthermore, after accounting for herbivorous insect richness in the multiple regression, mean annual temperature, and plant-available water remained significantly correlated with the occurrence of glandular trichomes. Plants growing in stressful habitats generally grow more slowly, and may adopt higher levels of resistance to herbivores (growth-defense trade-off hypothesis; ), because the cost of losing tissues in these habitats is relatively higher (e.g., ). This can explain the negative correlation between incidence of plant species with glandular trichomes and mean annual temperature. Production of glandular trichomes is costly, for example, in some plant species, production of glandular trichomes has been documented to be limited by soil water content (e.g., ; ). This can partly explain the positive correlation of incidence of glandular trichomes with plant-available water found in our study. However, compared with plant-available water, mean annual temperature explained 30% more of the variation in the incidence of species with glandular trichomes. Thus, investing more resources in defense is favored under more stressful conditions (i.e., lower temperature found at higher elevations in this study; ; ).
It is important to note that many plant traits have dual roles in protection against herbivores and environmental stressors (Moles et al., 2011). For example, in addition to protecting from herbivory, it was found that glandular trichomes play an important role in protecting plants against intense UV radiation by accumulating flavonoids (Tattini et al., 2000). Unfortunately, we cannot get data on UV radiation with sufficient high resolution along the elevational gradients in our study area. However, UV radiation has long been confirmed to increase with increasing elevation due to smaller UV absorbing air mass (). In response to higher UV radiation, plants at high elevations should have a higher level of defenses against UV radiation than plants at low elevations (; Song et al., 2020), which may also partly explain the positive correlation between the incidence of species with glandular trichomes and elevation observed in this study. Thus, in order to improve our understanding of the underlying mechanisms associated with elevational variation in the occurrence of glandular trichomes, further studies should be conducted using manipulative experiments along elevational gradients, after controlling for other relevant factors.
Some intra-specific studies have found that glandular trichomes (i.e., density and composition) may shift within a species’ range (e.g., ). Unfortunately, our study cannot account for this because we do not have corresponding data. However, as a parallel line of research, inter-specific comparisons across broad spatial, and multi-species scale are necessary for uncovering the macro-evolutionary patterns in plant defenses, which is especially important for understanding the elevational patterns of plant-herbivore interactions and the resultant elevational diversity gradients (; Schemske et al., 2009). Thus, although our analyses were based on simplified binary data (presence/absence of glandular trichomes), we have been able to provide important insights into macro-evolutionary patterns in investments of glandular trichomes by plants along elevational gradients in the Hengduan Mountains region. In order to better understand the underlying mechanisms driving macro-ecological patterns of glandular trichomes in this region, further studies assessing geographic patterns of plant intra-specific variation should be conducted.
Conclusion
In one of the largest empirical study of a single defense trait along elevational gradients to date, with more than 6,000 plant species and spanning 4,200 m of elevation, we found that plants from higher elevations had higher levels of defense against herbivores via glandular trichomes than did low-elevation plants, a finding that is contrary to the prevailing view on this topic (e.g., Rasmann et al., 2014). In addition to contributing important data to the current debate on the generality of the elevational herbivory defense hypothesis, our results highlight the importance of considering the simultaneous effects of abiotic drivers in understanding geographic variation in plant defenses (or any other traits), that have usually been overlooked ().
Statements
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author contributions
RW collected the data. RW and LS performed the analyses. BS and HS conceived the original research plan. BS, RW, and SL-Y wrote the article. All authors read and approved the final manuscript.
Funding
This work was supported by the Second Tibetan Plateau Scientific Expedition and Research (STEP) program (2019QZKK0502), the National Natural Science Foundation of China (31770249 and 32071669), the Young Academic and Technical Leader Raising Foundation of Yunnan Province (2017HB062), the Youth Innovation Promotion Association CAS (2017437), and the CAS “Light of West China” Program to BS.
Acknowledgments
We thank Susan Lev-Yadun and the reviewers for their constructive comments on the manuscript.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2021.632464/full#supplementary-material
Footnotes
References
1
Abdala-RobertsL.MoreiraX.RasmannS.Parra-TablaV.MooneyK. A. (2016a). Test of biotic and abiotic correlates of latitudinal variation in defenses in the perennial herb Ruellia nudiflora.J. Ecol.104580–590. 10.1111/1365-2745.12512
2
Abdala-RobertsL.RasmannS.Berny-Mier Y TeránJ. C.CoveloF.GlauserG.MoreiraX. (2016b). Biotic and abiotic factors associated with altitudinal variation in plant traits and herbivory in a dominant oak species.Am. J. Bot.1032070–2078. 10.3732/ajb.1600310
3
Alonso-AmelotM. E.Oliveros-BastidasA.Calcagno-PisarelliM. P. (2007). Phenolics and condensed tannins of high altitude Pteridium arachnoideum in relation to sunlight exposure, elevation, and rain regime.Biochem. Syst. Ecol.351–10. 10.1016/j.bse.2006.04.013
4
AmbrosioS.OkiY.HelenoV.ChavesJ.BarboniP.LichstonJ.et al (2008). Constituents of glandular trichomes of Tithonia diversifolia: relationships to herbivory and antifeedant activity.Phytochemistry692052–2060. 10.1016/j.phytochem.2008.03.019
5
BartonK. E. (2016). Tougher and thornier: general patterns in the induction of physical defense traits.Funct. Ecol.30181–187. 10.1111/1365-2435.12495
6
BaurB.MeierT.BaurA.SchmeraD. (2014). Terrestrial gastropod diversity in an alpine region: disentangling effects of elevation, area, geometric constraints, habitat type and land-use intensity.Ecography37390–401. 10.1111/j.1600-0587.2013.00312.x
7
BecerraJ. (2007). The impact of herbivore-plant coevolution on plant community structure.Proc. Natl. Acad. Sci. U.S.A.1047483–7488. 10.1073/pnas.0608253104
8
BhattaK.GrytnesJ. A.VetaasO. (2018). Scale sensitivity of the relationship between alpha and gamma diversity along an alpine elevation gradient in central Nepal.J. Biogeography45804–814. 10.1111/jbi.13188
9
BhattaraiK.VetaasO.GrytnesJ. A. (2004). Relationship between plant species richness and biomass in an arid sub-alpine grassland of the central Himalayas. Nepal.Folia Geobot.3957–71. 10.1007/BF02803264
10
CardinaleB.SrivastavaD.DuffyJ.WrightJ.DowningA.SankaranM.et al (2006). Effects of biodiversity on the functioning of trophic groups and ecosystems.Nature443989–992. 10.1038/nature05202
11
ChenS. C.CornwellW.ZhangH. X.MolesA. (2017). Plants show more flesh in the tropics: Variation in fruit type along latitudinal and climatic gradients.Ecography40531–538. 10.1111/ecog.02010
12
ChenS. X. (1992). Insects of the Hengduan Mountains Region. Beijing: Science Press.
13
ColeyP.BryantJ.ChapinF. S. (1985). Resource availability and plant antiherbivore defense.Science230895–899. 10.1126/science.230.4728.895
14
CorkS.KrockenbergerA. (1991). Methods and pitfalls of extracting condensed tannins and other phenolics from plants: insights from investigations on Eucalyptus leaves.J. Chem. Ecol.17123–134. 10.1007/BF00994426
15
De LongJ.SundqvistM.GundaleM.GieslerR.WardleD. (2015). Effects of elevation and nitrogen and phosphorus fertilization on plant defense compounds in subarctic tundra heath vegetation.Funct. Ecol.30314–325. 10.1111/1365-2435.12493
16
DefossezE.PellissierL.RasmannS. (2018). The unfolding of plant growth form-defense syndromes along elevation gradients.Ecol. Lett.21609–618. 10.1111/ele.12926
17
DescombesP.MarchonJ.PradervandJ. N.BilatJ.GuisanA.RasmannS.et al (2017). Community-level plant palatability increases with elevation as insect herbivore abundance declines.J. Ecol.105142–151. 10.1111/1365-2745.12664
18
DimockM.KennedyG. (1983). The role of glandular trichomes in the resistance of Lycopersicon hirsutum f. Glabratum to Heliothis zea.Entomol. Expe. Appl.33263–268. 10.1111/j.1570-7458.1983.tb03267.x
19
ElleE.van DamN. M.HareJ. D. (1999). Cost of glandular trichomes, a “Resistance” character in Datura wrightii Regel (Solanaceae).Evolution5322–35. 10.2307/2640917
20
EndaraM. J.ColeyP. (2011). The resource availability hypothesis revisited: a meta-analysis.Funct. Ecol.25389–398. 10.1111/j.1365-2435.2010.01803.x
21
FickS.HijmansR. (2017). WorldClim 2: new 1-km spatial resolution climate surfaces for global land areas.Int. J. Climatol.374302–4315. 10.1002/joc.5086
22
FutuymaD. J.AgrawalA. A. (2009). Evolutionary history and species interactions.Proc. Natl. Acad. Sci. U.S.A.10618043–18044. 10.1073/pnas.0910334106
23
GalmánA.Abdala-RobertsL.ZhangS.BernyJ.RasmannS.MoreiraX. (2018). A global analysis of elevational gradients in leaf herbivory and its underlying drivers: effects of plant growth form leaf habit and climatic correlates.J. Ecol.106413–421. 10.1111/1365-2745.12866
24
GaoB. Y.FangW. P.KongX. J.XuJ. M.GuanZ. T.YangJ. L.et al (2003). Flora of Sichuan. Chengdu: Sichuan people’s Publishing House.
25
GaoY. Q.ZhengJ. X.LinX. Q.DuF. (2020). Distribution patterns of clonal plants in the subnival belt of the Hengduan Mountains SW China.Plant Div.42386–392. 10.1016/j.pld.2020.06.006
26
GuimarãesC.VianaJ.CornelissenT. (2014). A meta-analysis of the effects of fragmentation on herbivorous insects.Environ. Entomol.43537–545. 10.1603/EN13190
27
HaasS. E.ReevesM. K.PinkneyA. E.JohnsonP. T. J. (2018). Continental-extent patterns in amphibian malformations linked to parasites, chemical contaminants, and their interactions.Glob. Change Biol.24e275–e288. 10.1111/gcb.13908
28
HanleyM.LamontB.FairbanksM.RaffertyC. (2007). Plant structural traits and their role in anti-herbivore defense.Perspect. Plant Ecol. Evol. Syst.8157–178. 10.1016/j.ppees.2007.01.001
29
HawkinsB. A.PorterE. E. (2003). Relative influences of current and historical factors on mammal and bird diversity patterns in deglaciated North America.Glob. Ecol. Biogeography12475–481. 10.1046/j.1466-822X.2003.00060.x
30
HermsD.MattsonW. (1992). The dilemma of plants: to grow or defend.Q. Rev. Biol.67283–335. 10.1086/417659
31
HoL. S. T.AneC. (2014). A linear-time algorithm for gaussian and non-gaussian trait evolution models.Syst. Biol.63397–408. 10.1093/sysbio/syu005
32
HorganF. G.QuiringD. T.LagnaouiA.PelletierY. (2009). Effects of altitude of origin on trichome-mediated anti-herbivore resistance in wild andean potatoes.Flora20449–62. 10.1016/j.flora.2008.01.008
33
HrdinaA.RomportlD. (2017). Evaluating global biodiversity hotspots – very rich and even more endangered.J. Land. Ecol.10108–115. 10.1515/jlecol-2017-0013
34
HuY.JinK.HuangZ.DingZ. F.LiangJ. C.PanX. Y.et al (2017). Elevational patterns of non-volant small mammal species richness in gyirong valley central Himalaya: evaluating multiple spatial and environmental drivers.J. Biogeography442764–2777. 10.1111/jbi.13102
35
IvesA. R.GarlandT. (2010). Phylogenetic logistic regression for binary dependent variables.Syst. Biol.599–26. 10.1093/sysbio/syp074
36
JanzenD. (1974). Tropical blackwater rivers, animals, and mast fruiting by the Dipterocarpaceae.Biotropica669–103. 10.2307/2989823
37
JinY.QianH. (2019). V.PhyloMaker: an R package that can generate very large phylogenies for vascular plants.Ecography421353–1359. 10.1111/ecog.04434
38
KennedyG.SorensonC. (1985). Role of glandular trichomes in the resistance of Lycopersicon hirsutum f. glabratum to Colorado Potato Beetle (coleoptera: chrysomelidae).J. Econ. Entomol.78547–551. 10.1093/jee/78.3.547
39
KirkhamM. B. (2005). Principles of Soil and Plant Water Relations, 2nd Edn. Amsterdam: Elsevier.
40
KörnerC. (2003). Alpine plant life: functional plant ecology of high mountain ecosystems.J. Ecol.88727–731. 10.1007/978-3-642-18970-8
41
LauterD. J.MunnsD. N. (1986). Water loss via the glandular trichomes of Chickpea (Cicer arietinum L.).J. Exp. Botany37640–649. 10.1093/jxb/37.5.640
42
Lev-YadunS. (2014). Do plants use visual carrion-based aposematism to deter herbivores?Biochemist3636–39. 10.1042/BIO03605036
43
LiS.TosensT.HarleyP.JiangY. F.KanagendranA.GrosbergM.et al (2018). Glandular trichomes as a barrier against atmospheric oxidative stress: relationships with ozone uptake leaf damage and emission of LOX products across a diverse set of species.Plant Cell Environ.411263–1277. 10.1111/pce.13128
44
LiW. H.ZhangY. G. (2010). Vertical Climate and Its Effect on Forest Distribution in the Hengduan Mountains. Beijing: China Meteorological Press.
45
LiewT. S.SchilthuizenM.LakimM. (2010). The determinants of land snail diversity along a tropical elevational gradient: insularity, geometry and niches.J. Biogeography371071–1078. 10.1111/j.1365-2699.2009.02243.x
46
LihavainenJ.HassinenV.Keski-SaariS.SoberA.OksanenE.KeinänenM. (2017). Low vapor pressure deficit reduces glandular trichome density and modifies the chemical composition of cuticular waxes in silver birch leaves.Tree Physiol.371–16. 10.1093/treephys/tpx045
47
LiuS. W.et al (1996). Flora of Qinghai. Xi’ning: Qinghai People’s Publishing House.
48
LiuY.JingS. X.LuoS. H.LiS. H. (2019). Non-volatile natural products in plant glandular trichomes: Chemistry, biological activities and biosynthesis.Nat. Prod. Rep.36626–665. 10.1039/C8NP00077H
49
McFaddenD. (1974). Conditional Logit Analysis of Qualitative Choice Behavior. Frontiers in Econometrics.Massachusetts: Academic Press, 105–142.
50
MetcalfeD.AsnerG.MartinR.SilvaE. J.HuaracaH. W.AmézquitaF.et al (2014). Herbivory makes major contributions to ecosystem carbon and nutrient cycling in tropical forests.Ecol. Lett.17324–332. 10.1111/ele.12233
51
MolesA.LaffanS.KeigheryM.DalrympleR.TindallM.ChenS. C. (2020). A hairy situation: plant species in warm sunny places are more likely to have pubescent leaves.J. Biogeography471934–1944. 10.1111/jbi.13870
52
MolesA.WallisI.FoleyW.WartonD.StegenJ.BisigatoA.et al (2011). Putting plant resistance traits on the map: a test of the idea that plants are better defended at lower latitudes.New Phytol.191777–788. 10.1111/j.1469-8137.2011.03732.x
53
MoreiraX.PetryW.MooneyK.RasmannS.Abdala-RobertsL. (2018). Elevational gradients in plant defenses and insect herbivory: recent advances in the field and prospects for future research.Ecography411485–1496. 10.1111/ecog.03184
54
PellissierL.FiedlerK.NdribeC.DubuisA.PradervandJ. N.GuisanA.et al (2012). Shifts in species richness, herbivore specialization, and plant resistance along elevation gradients.Ecol. Evol.21818–1825. 10.1002/ece3.296
55
PellissierL.RogerA.BilatJ.RasmannS. (2014). High elevation Plantago lanceolata plants are less resistant to herbivory than their low elevation conspecifics: is it just temperature?Ecography37950–959. 10.1111/ecog.00833
56
Pincheira-DonosoD.JaraM.ReAneyA.García-RoaR.SaldarriagaM.HodgsonD. (2017). Hypoxia and hypothermia as rival agents of selection driving the evolution of viviparity in lizards.Glob. Ecol. Biogeography261238–1246. 10.1111/geb.12626
57
PrattJ.Keefover-RingK.LiuL.MooneyK. (2014). Genetically based latitudinal variation in Artemisia californica secondary chemistry.Oikos123953–963. 10.1111/oik.01156
58
PuterkaG.FaroneW.PalmerT.BarringtonA. (2003). Structure-function relationships affecting the insecticidal and miticidal activity of sugar esters.J. Econ. Entomol.96636–644. 10.1603/0022-0493-96.3.636
59
R Development Core Team. (2020). R: A Language and Environment for Statistical Computing.Vienna: R Foundation for Statistical Computing.
60
RaffaK.PowellE.TownsendP. (2013). Temperature-driven range expansion of an irruptive insect heightened by weakly coevolved plant defenses.Proc. Natl. Acad. Sci. U.S.A.1102193–2198. 10.1073/pnas.1216666110
61
RasmannS.PellissierL.DefossezE.JactelH.KunstlerG. (2014). Climate-driven change in plant–insect interactions along elevation gradients.Funct. Ecol.2846–54. 10.1111/1365-2435.12135
62
ReynoldsB.CrossleyD. (1997). Spatial variation in herbivory by forest canopy arthropods along an elevation gradient.Environ. Entomol.261232–1239. 10.1093/ee/26.6.1232
63
RonelM.Lev-YadunS. (2012). The spiny, thorny and prickly plants in the flora of Israel.Bot. J. Linnean Soc.168344–352. 10.1111/j.1095-8339.2011.01211.x
64
SchemskeD.MittelbachG.CornellH.SobelJ.RoyK. (2009). Is there a latitudinal gradient in the importance of biotic interactions?Ann. Rev. Ecol. Syst.40245–269. 10.1146/annurev.ecolsys.39.110707.173430
65
SchmeraD.BaurB. (2014). Gastropod communities in alpine grasslands are characterized by high beta diversity.Commun. Ecol.15246–255. 10.1556/ComEc.15.2014.2.13
66
SchoonhovenL. M.van LoonJ.DickeM. (2005). Insect-Plant Biology.Oxford: Oxford University Press.
67
SilvaD.HetheringtonA.MansfieldT. A. (1996). Where does all the calcium go? Evidence of an important regulatory role for trichomes in two calcicoles.Plant Cell Environ.19880–886. 10.1111/j.1365-3040.1996.tb00424.x
68
SimmonsA.GurrG.McGrathD.MartinP.NicolH. (2004). Entrapment of Helicoverpa armigera (Hübner) (Lepidoptera: Noctuidae) on glandular trichomes of Lycopersicon species.Aust. J. Entomol.43196–200. 10.1111/j.1440-6055.2004.00414.x
69
SmithS.BrownJ. (2018). Constructing a broadly inclusive seed plant phylogeny.Am. J. Bot.105302–314. 10.1002/ajb2.1019
70
SongB.GaoY. Q.StöcklinJ.SongM. S.SunL.SunH. (2020). Ultraviolet screening increases with elevation in translucent bracts of Rheum nobile (Polygonaceae), an alpine ‘glasshouse’ plant from the high Himalayas.Bot. J. Linnean Soc.2276–286. 10.1093/botlinnean/boaa005
71
TattiniM.GravanoE.PinelliP.MulinacciN.RomaniA. (2000). Flavonoids accumulate in leaves and glandular trichomes of Phillyrea latifolia exposed to excess solar radiation.New Phytol.14869–77. 10.1046/j.1469-8137.2000.00743.x
72
ThompsonK.AskewA.GrimeJ.DunnettN.WillisA. J. (2005). Biodiversity ecosystem function and plant traits in mature and immature plant communities.Funct. Ecol.19355–358. 10.1111/j.0269-8463.2005.00936.x
73
TindallM. L.ThomsonF. J.LaffanS. W.MolesA. T. (2016). Is there a latitudinal gradient in the proportion of species with spinescence?J. Plant Ecol.10294–300. 10.1093/jpe/rtw031
74
TomasJ.GilL.Llorens-MolinaJ. A.CardonaC.GarcíaM. T.LlorensL. (2019). Biogenic volatiles of rupicolous plants act as direct defenses against molluscs: the case of the endangered Clinopodium rouyanum.Flora258:151428. 10.1016/j.flora.2019.151428
75
TozinL. R.MarquesM.RodriguesT. (2017). Herbivory by leaf-cutter ants changes the glandular trichomes density and the volatile components in an aromatic plant model.AoB Plants9:lx057. 10.1093/aobpla/plx057
76
VetaasO.GrytnesJ. A. (2002). Distribution of vascular plant species richness and endemic richness along the Himalayan elevation gradient in Nepal.Glob. Ecol. Biogeography11291–301. 10.1046/j.1466-822X.2002.00297.x
77
WagnerG. J. (1991). Secreting glandular trichomes: more than just hairs.Plant Physiol.96675–679. 10.1104/pp.96.3.675
78
WangW. C. (1993). Vascular Flora of Hengduan Mountains. Xining: Science Press.
79
WesterberghA.NybergA. B. (1995). Selective grazing of hairless Silene dioica plants by land gastropods.Oikos73289–298. 10.2307/3545953
80
Wu1977–2006Flora of Yunnan.Beijing: Science Publishing House.
81
WuZ. Y. (1983–1987). Flora of Tibet. Beijing: Science Publishing House.
82
WuZ. Y.PengH.LiD. Z. (1994–2012). Flora of China. Beijing: Science Publishing House.
83
XuQ.Lev-YadunS.SunL.ChenZ.SongB.SunH. (2020). Spinescent patterns in the flora of Jiaozi Snow Mountain. Southwestern China.Plant Div.283–91. 10.1016/j.pld.2019.12.002
84
ZhangD. C.ZhangY. H.BouffordD. E.SunH. (2009). Elevational patterns of species richness and endemism for some important taxa in the Hengduan Mountains, southwestern China.Biodivers. Conserv.18699–716. 10.1007/s10531-008-9534-x
85
ZhangS.ZhangY.MaK. (2016). Latitudinal variation in herbivory: hemispheric asymmetries and the role of climatic drivers.J. Ecol.1041089–1095. 10.1111/1365-2745.12588
86
ZhangY.SchaapM. G.ZhaY. (2018). A high-resolution global map of soil hydraulic properties produced by a hierarchical parameterization of a physically-based water retention model.Water Res. Res.549774–9790. 10.1029/2018WR023539
87
ZhangY. L.LiB. Y.ZhengD. (2002). A discussion on the boundary and area of the tibetan plateau in China.Geogr. Res.211–10.
88
ZhaoY.XuF. L.LiuJ.GuanF. C.QuanH.MengF. J. (2019). The adaptation strategies of Herpetospermum pedunculosum (Ser.) Baill at altitude gradient of the tibetan plateau by physiological and metabolomic methods.BMC Genomics20:451. 10.1186/s12864-019-5778-y
Summary
Keywords
biotic interaction, elevation, glandular trichome, growth form, herbivorous insect, water availability, temperature
Citation
Wu R, Lev-Yadun S, Sun L, Sun H and Song B (2021) Higher Elevations Tend to Have Higher Proportion of Plant Species With Glandular Trichomes. Front. Plant Sci. 12:632464. doi: 10.3389/fpls.2021.632464
Received
23 November 2020
Accepted
02 March 2021
Published
12 April 2021
Volume
12 - 2021
Edited by
Tiina Tosens, Estonian University of Life Sciences, Estonia
Reviewed by
Béla Tóthmérész, University of Debrecen, Hungary; Amber Churchill, Western Sydney University, Australia
Updates
Copyright
© 2021 Wu, Lev-Yadun, Sun, Sun and Song.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Hang Sun, hsun@mail.kib.ac.cnBo Song, songbo@mail.kib.ac.cn
This article was submitted to Functional Plant Ecology, a section of the journal Frontiers in Plant Science
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.