Rui Wu1,2
Simcha Lev-Yadun
Hang Sun
Bo Song- 1Key Laboratory of Resource Biology and Biotechnology in Western China, Ministry of Education, College of Life Sciences, Northwest University, Xi’an, China
- 2Key Laboratory for Plant Diversity and Biogeography of East Asia, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, China
- 3Department of Biology and Environment, Faculty of Natural Sciences, University of Haifa at Oranim, Kiryat Tiv’on, Israel
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; Abdala-Roberts et al., 2016a; Descombes et al., 2017). 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 (Descombes et al., 2017), positive (e.g., Moreira et al., 2018), or no relationship between elevation and defense (Alonso-Amelot et al., 2007). 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 (Barton, 2016). 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; Abdala-Roberts et al., 2016a). 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 (Cork and Krockenberger, 1991). 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; Defossez et al., 2018; 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 (Abdala-Roberts et al., 2016a). For example, concentration of phenolics in leaves of Quercus robur increased independently with decreasing temperature toward higher elevations (Abdala-Roberts et al., 2016b). 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; Li et al., 2018), one of their most important well-known roles is to serve as a defense against herbivores, especially insects (Dimock and Kennedy, 1983; Wagner, 1991; Hanley et al., 2007). Insect herbivores can be killed directly by ingested toxic exudates produced by glandular trichomes (Dimock and Kennedy, 1983; 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., Dimock and Kennedy, 1983; Kennedy and Sorenson, 1985; Ambrosio et al., 2008; Lev-Yadun, 2014; 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., Horgan et al., 2009; 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 (Hrdina and Romportl, 2017). 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. The location of the Hengduan Mountains in southwestern China.
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 (Gao et al., 2003), 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 (Bhattarai et al., 2004). This is a commonly used method (e.g., Vetaas and Grytnes, 2002; Bhattarai et al., 2004; Zhang et al., 2009; Bhatta et al., 2018), 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 (Galmán et al., 2018). 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 (Horgan et al., 2009), 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 (Cardinale et al., 2006; Guimarães et al., 2014). We collected data on the elevational range for herbivorous insects in this area from Insects of the Hengduan Mountains Region (Chen, 1992) 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 (Hawkins and Porter, 2003; Zhang et al., 2016). Thus, we download data of mean annual temperature (MAT; °C) from the WorldClim version2 database (Fick and Hijmans, 2017). 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 (Kirkham, 2005), 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 Chen et al. (2017) 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 (Chen et al., 2017; Gao et al., 2020). 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 (Haas et al., 2018). 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 (Hu et al., 2017).
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 (Ho and Ane, 2014). 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 (Ives and Garland, 2010). For these analyses, we generated a phylogenetic tree of study species using “V.PhyloMaker” package (Jin and Qian, 2019), 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 (Chen et al., 2017; Gao et al., 2020).
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. The relationships between glandular trichomes of plant species and elevation for every 100 m belt in the Hengduan Mountains region. The graph was visualized according to the proportion of plants with glandular trichomes, but the analysis was based on the binary data of presence versus absence of glandular trichomes. The fitted line (solid line) and estimated 95% confidence interval (dashed line) display the predicted probability of plants with glandular trichomes as fit by logistic regression models. α and β value was calculated after elevation being scaled.
Table 1. 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. Relationship between elevation and herbivorous insect richness (A), mean annual temperature (B), and plant-available water (C) for every 100 m belt in the Hengduan Mountains region. Each point represents the value of herbivorous insect richness or the environmental variable in the corresponding elevation belt. The fitted lines (solid line) and estimated 95% confidence intervals (dashed line) were calculated using regression models that included elevation as both linear and quadratic terms. α and β value was calculated after elevation being scaled.
Figure 4. The relationship between glandular trichome of plant species and herbivorous insect richness (A), mean annual temperature (B), and plant-available water (C) for every 100 m belt in the Hengduan Mountains region. The graphs were visualized according to the proportion of plants with glandular trichomes, but the analysis was based on the binary data of presence versus absence of plants with glandular trichomes. The fitted lines (solid line) and estimated 95% confidence intervals (dashed line) display the predicted probability of plants with glandular trichomes as fit by logistic regression models. α and β value was calculated after predictor variables being scaled.
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. 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.
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; Ives and Garland, 2010). 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 (Becerra, 2007). 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 (Cardinale et al., 2006; Guimarães et al., 2014). 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; Baur et al., 2014; 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., Abdala-Roberts et al., 2016a; 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 (Endara and Coley, 2011; De Long et al., 2015; 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; Herms and Mattson, 1992), because the cost of losing tissues in these habitats is relatively higher (e.g., Coley et al., 1985). 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., Lauter and Munns, 1986; Elle et al., 1999). 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; Janzen, 1974; Coley et al., 1985).
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 (Körner, 2003). 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 (Körner, 2003; 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., Horgan et al., 2009). 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 (Futuyma and Agrawal, 2009; 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 (Abdala-Roberts et al., 2016a).
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.
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.
Acknowledgments
We thank Susan Lev-Yadun and the reviewers for their constructive comments on the manuscript.
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
Abdala-Roberts, L., Moreira, X., Rasmann, S., Parra-Tabla, V., and Mooney, K. A. (2016a). Test of biotic and abiotic correlates of latitudinal variation in defenses in the perennial herb Ruellia nudiflora. J. Ecol. 104, 580–590. doi: 10.1111/1365-2745.12512
Abdala-Roberts, L., Rasmann, S., Berny-Mier Y Terán, J. C., Covelo, F., Glauser, G., and Moreira, X. (2016b). Biotic and abiotic factors associated with altitudinal variation in plant traits and herbivory in a dominant oak species. Am. J. Bot. 103, 2070–2078. doi: 10.3732/ajb.1600310
Alonso-Amelot, M. E., Oliveros-Bastidas, A., and Calcagno-Pisarelli, M. P. (2007). Phenolics and condensed tannins of high altitude Pteridium arachnoideum in relation to sunlight exposure, elevation, and rain regime. Biochem. Syst. Ecol. 35, 1–10. doi: 10.1016/j.bse.2006.04.013
Ambrosio, S., Oki, Y., Heleno, V., Chaves, J., Barboni, P., Lichston, J., et al. (2008). Constituents of glandular trichomes of Tithonia diversifolia: relationships to herbivory and antifeedant activity. Phytochemistry 69, 2052–2060. doi: 10.1016/j.phytochem.2008.03.019
Barton, K. E. (2016). Tougher and thornier: general patterns in the induction of physical defense traits. Funct. Ecol. 30, 181–187. doi: 10.1111/1365-2435.12495
Baur, B., Meier, T., Baur, A., and Schmera, D. (2014). Terrestrial gastropod diversity in an alpine region: disentangling effects of elevation, area, geometric constraints, habitat type and land-use intensity. Ecography 37, 390–401. doi: 10.1111/j.1600-0587.2013.00312.x
Becerra, J. (2007). The impact of herbivore-plant coevolution on plant community structure. Proc. Natl. Acad. Sci. U.S.A. 104, 7483–7488. doi: 10.1073/pnas.0608253104
Bhatta, K., Grytnes, J. A., and Vetaas, O. (2018). Scale sensitivity of the relationship between alpha and gamma diversity along an alpine elevation gradient in central Nepal. J. Biogeography 45, 804–814. doi: 10.1111/jbi.13188
Bhattarai, K., Vetaas, O., and Grytnes, J. A. (2004). Relationship between plant species richness and biomass in an arid sub-alpine grassland of the central Himalayas. Nepal. Folia Geobot. 39, 57–71. doi: 10.1007/BF02803264
Cardinale, B., Srivastava, D., Duffy, J., Wright, J., Downing, A., Sankaran, M., et al. (2006). Effects of biodiversity on the functioning of trophic groups and ecosystems. Nature 443, 989–992. doi: 10.1038/nature05202
Chen, S. C., Cornwell, W., Zhang, H. X., and Moles, A. (2017). Plants show more flesh in the tropics: Variation in fruit type along latitudinal and climatic gradients. Ecography 40, 531–538. doi: 10.1111/ecog.02010
Coley, P., Bryant, J., and Chapin, F. S. (1985). Resource availability and plant antiherbivore defense. Science 230, 895–899. doi: 10.1126/science.230.4728.895
Cork, S., and Krockenberger, A. (1991). Methods and pitfalls of extracting condensed tannins and other phenolics from plants: insights from investigations on Eucalyptus leaves. J. Chem. Ecol. 17, 123–134. doi: 10.1007/BF00994426
De Long, J., Sundqvist, M., Gundale, M., Giesler, R., and Wardle, D. (2015). Effects of elevation and nitrogen and phosphorus fertilization on plant defense compounds in subarctic tundra heath vegetation. Funct. Ecol. 30, 314–325. doi: 10.1111/1365-2435.12493
Defossez, E., Pellissier, L., and Rasmann, S. (2018). The unfolding of plant growth form-defense syndromes along elevation gradients. Ecol. Lett. 21, 609–618. doi: 10.1111/ele.12926
Descombes, P., Marchon, J., Pradervand, J. N., Bilat, J., Guisan, A., Rasmann, S., et al. (2017). Community-level plant palatability increases with elevation as insect herbivore abundance declines. J. Ecol. 105, 142–151. doi: 10.1111/1365-2745.12664
Dimock, M., and Kennedy, G. (1983). The role of glandular trichomes in the resistance of Lycopersicon hirsutum f. Glabratum to Heliothis zea. Entomol. Expe. Appl. 33, 263–268. doi: 10.1111/j.1570-7458.1983.tb03267.x
Elle, E., van Dam, N. M., and Hare, J. D. (1999). Cost of glandular trichomes, a “Resistance” character in Datura wrightii Regel (Solanaceae). Evolution 53, 22–35. doi: 10.2307/2640917
Endara, M. J., and Coley, P. (2011). The resource availability hypothesis revisited: a meta-analysis. Funct. Ecol. 25, 389–398. doi: 10.1111/j.1365-2435.2010.01803.x
Fick, S., and Hijmans, R. (2017). WorldClim 2: new 1-km spatial resolution climate surfaces for global land areas. Int. J. Climatol. 37, 4302–4315. doi: 10.1002/joc.5086
Futuyma, D. J., and Agrawal, A. A. (2009). Evolutionary history and species interactions. Proc. Natl. Acad. Sci. U.S.A. 106, 18043–18044. doi: 10.1073/pnas.0910334106
Galmán, A., Abdala-Roberts, L., Zhang, S., Berny, J., Rasmann, S., and Moreira, X. (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. 106, 413–421. doi: 10.1111/1365-2745.12866
Gao, B. Y., Fang, W. P., Kong, X. J., Xu, J. M., Guan, Z. T., Yang, J. L., et al. (2003). Flora of Sichuan. Chengdu: Sichuan people’s Publishing House.
Gao, Y. Q., Zheng, J. X., Lin, X. Q., and Du, F. (2020). Distribution patterns of clonal plants in the subnival belt of the Hengduan Mountains SW China. Plant Div. 42, 386–392. doi: 10.1016/j.pld.2020.06.006
Guimarães, C., Viana, J., and Cornelissen, T. (2014). A meta-analysis of the effects of fragmentation on herbivorous insects. Environ. Entomol. 43, 537–545. doi: 10.1603/EN13190
Haas, S. E., Reeves, M. K., Pinkney, A. E., and Johnson, P. T. J. (2018). Continental-extent patterns in amphibian malformations linked to parasites, chemical contaminants, and their interactions. Glob. Change Biol. 24, e275–e288. doi: 10.1111/gcb.13908
Hanley, M., Lamont, B., Fairbanks, M., and Rafferty, C. (2007). Plant structural traits and their role in anti-herbivore defense. Perspect. Plant Ecol. Evol. Syst. 8, 157–178. doi: 10.1016/j.ppees.2007.01.001
Hawkins, B. A., and Porter, E. E. (2003). Relative influences of current and historical factors on mammal and bird diversity patterns in deglaciated North America. Glob. Ecol. Biogeography 12, 475–481. doi: 10.1046/j.1466-822X.2003.00060.x
Herms, D., and Mattson, W. (1992). The dilemma of plants: to grow or defend. Q. Rev. Biol. 67, 283–335. doi: 10.1086/417659
Ho, L. S. T., and Ane, C. (2014). A linear-time algorithm for gaussian and non-gaussian trait evolution models. Syst. Biol. 63, 397–408. doi: 10.1093/sysbio/syu005
Horgan, F. G., Quiring, D. T., Lagnaoui, A., and Pelletier, Y. (2009). Effects of altitude of origin on trichome-mediated anti-herbivore resistance in wild andean potatoes. Flora 204, 49–62. doi: 10.1016/j.flora.2008.01.008
Hrdina, A., and Romportl, D. (2017). Evaluating global biodiversity hotspots – very rich and even more endangered. J. Land. Ecol. 10, 108–115. doi: 10.1515/jlecol-2017-0013
Hu, Y., Jin, K., Huang, Z., Ding, Z. F., Liang, J. C., Pan, X. 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. Biogeography 44, 2764–2777. doi: 10.1111/jbi.13102
Ives, A. R., and Garland, T. (2010). Phylogenetic logistic regression for binary dependent variables. Syst. Biol. 59, 9–26. doi: 10.1093/sysbio/syp074
Janzen, D. (1974). Tropical blackwater rivers, animals, and mast fruiting by the Dipterocarpaceae. Biotropica 6, 69–103. doi: 10.2307/2989823
Jin, Y., and Qian, H. (2019). V.PhyloMaker: an R package that can generate very large phylogenies for vascular plants. Ecography 42, 1353–1359. doi: 10.1111/ecog.04434
Kennedy, G., and Sorenson, C. (1985). Role of glandular trichomes in the resistance of Lycopersicon hirsutum f. glabratum to Colorado Potato Beetle (coleoptera: chrysomelidae). J. Econ. Entomol. 78, 547–551. doi: 10.1093/jee/78.3.547
Körner, C. (2003). Alpine plant life: functional plant ecology of high mountain ecosystems. J. Ecol. 88, 727–731. doi: 10.1007/978-3-642-18970-8
Lauter, D. J., and Munns, D. N. (1986). Water loss via the glandular trichomes of Chickpea (Cicer arietinum L.). J. Exp. Botany 37, 640–649. doi: 10.1093/jxb/37.5.640
Lev-Yadun, S. (2014). Do plants use visual carrion-based aposematism to deter herbivores? Biochemist 36, 36–39. doi: 10.1042/BIO03605036
Li, S., Tosens, T., Harley, P., Jiang, Y. F., Kanagendran, A., Grosberg, M., 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. 41, 1263–1277. doi: 10.1111/pce.13128
Li, W. H., and Zhang, Y. G. (2010). Vertical Climate and Its Effect on Forest Distribution in the Hengduan Mountains. Beijing: China Meteorological Press.
Liew, T. S., Schilthuizen, M., and Lakim, M. (2010). The determinants of land snail diversity along a tropical elevational gradient: insularity, geometry and niches. J. Biogeography 37, 1071–1078. doi: 10.1111/j.1365-2699.2009.02243.x
Lihavainen, J., Hassinen, V., Keski-Saari, S., Sober, A., Oksanen, E., and Keinänen, M. (2017). Low vapor pressure deficit reduces glandular trichome density and modifies the chemical composition of cuticular waxes in silver birch leaves. Tree Physiol. 37, 1–16. doi: 10.1093/treephys/tpx045
Liu, Y., Jing, S. X., Luo, S. H., and Li, S. H. (2019). Non-volatile natural products in plant glandular trichomes: Chemistry, biological activities and biosynthesis. Nat. Prod. Rep. 36, 626–665. doi: 10.1039/C8NP00077H
McFadden, D. (1974). Conditional Logit Analysis of Qualitative Choice Behavior. Frontiers in Econometrics. Massachusetts: Academic Press, 105–142.
Metcalfe, D., Asner, G., Martin, R., Silva, E. J., Huaraca, H. W., Amézquita, F., et al. (2014). Herbivory makes major contributions to ecosystem carbon and nutrient cycling in tropical forests. Ecol. Lett. 17, 324–332. doi: 10.1111/ele.12233
Moles, A., Laffan, S., Keighery, M., Dalrymple, R., Tindall, M., and Chen, S. C. (2020). A hairy situation: plant species in warm sunny places are more likely to have pubescent leaves. J. Biogeography 47, 1934–1944. doi: 10.1111/jbi.13870
Moles, A., Wallis, I., Foley, W., Warton, D., Stegen, J., Bisigato, A., 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. 191, 777–788. doi: 10.1111/j.1469-8137.2011.03732.x
Moreira, X., Petry, W., Mooney, K., Rasmann, S., and Abdala-Roberts, L. (2018). Elevational gradients in plant defenses and insect herbivory: recent advances in the field and prospects for future research. Ecography 41, 1485–1496. doi: 10.1111/ecog.03184
Pellissier, L., Fiedler, K., Ndribe, C., Dubuis, A., Pradervand, J. N., Guisan, A., et al. (2012). Shifts in species richness, herbivore specialization, and plant resistance along elevation gradients. Ecol. Evol. 2, 1818–1825. doi: 10.1002/ece3.296
Pellissier, L., Roger, A., Bilat, J., and Rasmann, S. (2014). High elevation Plantago lanceolata plants are less resistant to herbivory than their low elevation conspecifics: is it just temperature? Ecography 37, 950–959. doi: 10.1111/ecog.00833
Pincheira-Donoso, D., Jara, M., ReAney, A., García-Roa, R., Saldarriaga, M., and Hodgson, D. (2017). Hypoxia and hypothermia as rival agents of selection driving the evolution of viviparity in lizards. Glob. Ecol. Biogeography 26, 1238–1246. doi: 10.1111/geb.12626
Pratt, J., Keefover-Ring, K., Liu, L., and Mooney, K. (2014). Genetically based latitudinal variation in Artemisia californica secondary chemistry. Oikos 123, 953–963. doi: 10.1111/oik.01156
Puterka, G., Farone, W., Palmer, T., and Barrington, A. (2003). Structure-function relationships affecting the insecticidal and miticidal activity of sugar esters. J. Econ. Entomol. 96, 636–644. doi: 10.1603/0022-0493-96.3.636
R Development Core Team. (2020). R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing.
Raffa, K., Powell, E., and Townsend, P. (2013). Temperature-driven range expansion of an irruptive insect heightened by weakly coevolved plant defenses. Proc. Natl. Acad. Sci. U.S.A. 110, 2193–2198. doi: 10.1073/pnas.1216666110
Rasmann, S., Pellissier, L., Defossez, E., Jactel, H., and Kunstler, G. (2014). Climate-driven change in plant–insect interactions along elevation gradients. Funct. Ecol. 28, 46–54. doi: 10.1111/1365-2435.12135
Reynolds, B., and Crossley, D. (1997). Spatial variation in herbivory by forest canopy arthropods along an elevation gradient. Environ. Entomol. 26, 1232–1239. doi: 10.1093/ee/26.6.1232
Ronel, M., and Lev-Yadun, S. (2012). The spiny, thorny and prickly plants in the flora of Israel. Bot. J. Linnean Soc. 168, 344–352. doi: 10.1111/j.1095-8339.2011.01211.x
Schemske, D., Mittelbach, G., Cornell, H., Sobel, J., and Roy, K. (2009). Is there a latitudinal gradient in the importance of biotic interactions? Ann. Rev. Ecol. Syst. 40, 245–269. doi: 10.1146/annurev.ecolsys.39.110707.173430
Schmera, D., and Baur, B. (2014). Gastropod communities in alpine grasslands are characterized by high beta diversity. Commun. Ecol. 15, 246–255. doi: 10.1556/ComEc.15.2014.2.13
Schoonhoven, L. M., van Loon, J., and Dicke, M. (2005). Insect-Plant Biology. Oxford: Oxford University Press.
Silva, D., Hetherington, A., and Mansfield, T. A. (1996). Where does all the calcium go? Evidence of an important regulatory role for trichomes in two calcicoles. Plant Cell Environ. 19, 880–886. doi: 10.1111/j.1365-3040.1996.tb00424.x
Simmons, A., Gurr, G., McGrath, D., Martin, P., and Nicol, H. (2004). Entrapment of Helicoverpa armigera (Hübner) (Lepidoptera: Noctuidae) on glandular trichomes of Lycopersicon species. Aust. J. Entomol. 43, 196–200. doi: 10.1111/j.1440-6055.2004.00414.x
Smith, S., and Brown, J. (2018). Constructing a broadly inclusive seed plant phylogeny. Am. J. Bot. 105, 302–314. doi: 10.1002/ajb2.1019
Song, B., Gao, Y. Q., Stöcklin, J., Song, M. S., Sun, L., and Sun, H. (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. 2, 276–286. doi: 10.1093/botlinnean/boaa005
Tattini, M., Gravano, E., Pinelli, P., Mulinacci, N., and Romani, A. (2000). Flavonoids accumulate in leaves and glandular trichomes of Phillyrea latifolia exposed to excess solar radiation. New Phytol. 148, 69–77. doi: 10.1046/j.1469-8137.2000.00743.x
Thompson, K., Askew, A., Grime, J., Dunnett, N., and Willis, A. J. (2005). Biodiversity ecosystem function and plant traits in mature and immature plant communities. Funct. Ecol. 19, 355–358. doi: 10.1111/j.0269-8463.2005.00936.x
Tindall, M. L., Thomson, F. J., Laffan, S. W., and Moles, A. T. (2016). Is there a latitudinal gradient in the proportion of species with spinescence? J. Plant Ecol. 10, 294–300. doi: 10.1093/jpe/rtw031
Tomas, J., Gil, L., Llorens-Molina, J. A., Cardona, C., García, M. T., and Llorens, L. (2019). Biogenic volatiles of rupicolous plants act as direct defenses against molluscs: the case of the endangered Clinopodium rouyanum. Flora 258:151428. doi: 10.1016/j.flora.2019.151428
Tozin, L. R., Marques, M., and Rodrigues, T. (2017). Herbivory by leaf-cutter ants changes the glandular trichomes density and the volatile components in an aromatic plant model. AoB Plants 9:lx057. doi: 10.1093/aobpla/plx057
Vetaas, O., and Grytnes, J. A. (2002). Distribution of vascular plant species richness and endemic richness along the Himalayan elevation gradient in Nepal. Glob. Ecol. Biogeography 11, 291–301. doi: 10.1046/j.1466-822X.2002.00297.x
Wagner, G. J. (1991). Secreting glandular trichomes: more than just hairs. Plant Physiol. 96, 675–679. doi: 10.1104/pp.96.3.675
Westerbergh, A., and Nyberg, A. B. (1995). Selective grazing of hairless Silene dioica plants by land gastropods. Oikos 73, 289–298. doi: 10.2307/3545953
Xu, Q., Lev-Yadun, S., Sun, L., Chen, Z., Song, B., and Sun, H. (2020). Spinescent patterns in the flora of Jiaozi Snow Mountain. Southwestern China. Plant Div. 2, 83–91. doi: 10.1016/j.pld.2019.12.002
Zhang, D. C., Zhang, Y. H., Boufford, D. E., and Sun, H. (2009). Elevational patterns of species richness and endemism for some important taxa in the Hengduan Mountains, southwestern China. Biodivers. Conserv. 18, 699–716. doi: 10.1007/s10531-008-9534-x
Zhang, S., Zhang, Y., and Ma, K. (2016). Latitudinal variation in herbivory: hemispheric asymmetries and the role of climatic drivers. J. Ecol. 104, 1089–1095. doi: 10.1111/1365-2745.12588
Zhang, Y., Schaap, M. G., and Zha, Y. (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. 54, 9774–9790. doi: 10.1029/2018WR023539
Zhang, Y. L., Li, B. Y., and Zheng, D. (2002). A discussion on the boundary and area of the tibetan plateau in China. Geogr. Res. 21, 1–10.
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.
Edited by:
Tiina Tosens, Estonian University of Life Sciences, EstoniaReviewed by:
Béla Tóthmérész, University of Debrecen, HungaryAmber Churchill, Western Sydney University, Australia
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.cn; Bo Song, songbo@mail.kib.ac.cn