Spatial Patterns and Drivers of Angiosperm Sexual Systems in China Differ Between Woody and Herbaceous Species

Plant sexual systems play an important role in the evolution of angiosperm diversity. However, large-scale patterns in the frequencies of sexual systems (i.e. dioecy, monoecy, and hermaphroditism) and their drivers for species with different growth forms remain poorly known. Here, using a newly compiled database on the sexual systems and distributions of 19780 angiosperm species in China, we map the large-scale geographical patterns in frequencies of the sexual systems of woody and herbaceous species separately. We use these data to test the following two hypotheses: (1) the prevalence of sexual systems differs between woody and herbaceous assemblies because woody plants have taller canopies and are found in warm and humid climates; (2) the relative contributions of different drivers (specifically climate, evolutionary age, and mature plant height) to these patterns differ between woody and herbaceous species. We show that geographical patterns in proportions of different sexual systems (especially dioecy) differ between woody and herbaceous species. Geographical variations in sexual systems of woody species were influenced by climate, evolutionary age and plant height. In contrast, these have only weakly significant effects on the patterns of sexual systems of herbaceous species. We suggest that differences between species with woody and herbaceous growth forms in terms of biogeographic patterns of sexual systems, and their drivers, may reflect their differences in physiological and ecological adaptions, as well as the coevolution of sexual system with vegetative traits in response to environmental changes.

Introduction

The sexual systems of plant species play a significant role in the evolution of angiosperm diversity (Charlesworth and Wright, 2001; Nazareno et al., 2013; Käfer et al., 2014; Sabath et al., 2016; Rosche et al., 2018), variations in reproductive strategies (Bawa, 1980), and population and community dynamics in response to climate changes (Queenborough et al., 2009; Wang et al., 2018). Understanding the spatial distribution of plant sexual systems at large scales is vital for understanding the functional biogeography of reproductive traits. Previous studies have shown that the geographical distributions of angiosperm sexual systems may ultimately reflect coevolution with vegetative traits that are tied closely to sexual systems (i.e. growth form, Vamosi et al., 2003; Moeller et al., 2017; Wang et al., 2020). Nonetheless, the determinants of geographical patterns of sexual systems remain controversial.

The evolution of plant sexual systems has been widely hypothesized to be closely linked to the evolution of plant growth forms (e.g. woody vs. herbaceous growth forms, Bawa et al., 1985; Renner and Ricklefs, 1995; Vamosi et al., 2003; Vamosi and Vamosi, 2004; Mitchell and Diggle, 2005; Soza et al., 2012). Generally, dioecy is significantly more frequent in woody flowering plant species with long lifespan than in other growth forms (Renner and Ricklefs, 1995; Sakai et al., 1995; Vamosi et al., 2003; Ward et al., 2005). In contrast, hermaphroditism is more common in herbaceous species with short lifespan (Senarath, 2008; Moeller et al., 2017). This is likely because obligately outcrossing dioecious species need a long lifespan to find mates, reproduce and complete their life cycle (Morgan et al., 1997; Aarssen, 2000). In contrast, hermaphrodites, being able to self-pollinate, are likely to accumulate more genetic load and suffer from an increased turnover rate (Klekowski and Godfrey, 1989; Chen and Li, 2008). Moreover, growth form may affect relationships between traits and the environment (Fitzjohn et al., 2014). For example, based on more than 88,400 species and six plant traits, Šímová et al. (2018) have found significant differences in trait-climate correlations between woody and herbaceous species. However, the biogeographical distributions of sexual systems of different growth forms across large-scale environmental gradients remain poorly understood.

Plant sexual systems and sexual reproduction are sensitive to climate variations (Hedhly et al., 2009; Eckert et al., 2010). Recent studies indicate that reductions in precipitation can change the allocation of resources to male and female functions within individual plants (Hultine et al., 2016), which further influences the expression and composition of plant sexual system in local floras. For example, dioecy is more common in humid areas, especially in rainforest trees (Lloyd, 1980; Sakai et al., 1995; Friedman and Barrett, 2009; Vary et al., 2011), although aridity has been shown to have contributed to the evolutionary transition from hermaphroditism to dioecy (Ashman, 2006). Climate warming can also affect the composition of sexual systems in local floras by causing a mismatch between the timing of flowering and pollinator abundance/presence, which may lead to declines of plants populations (Etterson and Mazer, 2016). Moreover, different growth forms have contrasting climate preferences. Warm and humid climates are more suitable for woody plants (Wang et al., 2011), while herbaceous plants often have broader climate adaptations and hence have widespread distributions (Curtis and Bradley, 2016). Therefore, the influences of climate on sexual system compositions may differ between woody and herbaceous species. However, the role of climate in driving spatial variations in sexual system compositions of woody and herbaceous species remains poorly understood (Gamble et al., 2018).

Evolutionary history might be expected to influence the distribution of sexual systems (Obbard et al., 2006; Barrett, 2013; Käfer et al., 2014). Self-pollinated populations likely accumulate harmful mutants, which may increase their extinction rate (Etterson and Mazer, 2016). Therefore, self-pollinated species may be more common in temperate environments with more short-lived species and younger flora than in tropical environments. In contrast, dioecy (obligate outcrossing) can reduce the expression of recessive deleterious mutations, which tends to reduce their extinction rate (Obbard et al., 2006). Therefore, dioecy may be more frequent in regions with older floras. These findings suggest that the composition of sexual systems within communities is likely associated with the evolutionary ages of species in a region. Additionally, woody and herbaceous species differ significantly in their evolutionary rates (Laroche et al., 1997; Kay et al., 2006). Woody species have longer reproductive cycles and lifespans, and tend to accumulate genetic changes more slowly than herbaceous species. Differences in evolutionary rates may lead to differences in the speed of climate-niche evolution between these two growth forms. For example, using more than 5,000 plant species, Smith and Beaulieu (2009) found that woody plants adapted to new climates at a rate two to ten times slower than herbs over the course of their evolution. Therefore, exploring the differences between the effects of evolutionary time on frequencies of sexual systems of woody and herbaceous species would improve our understanding of the evolutionary mechanisms underlying large-scale patterns in sexual system composition.

Here, using data on sexual systems and spatial distributions of 19,780 angiosperm species across China, we compare the geographical patterns in plant sexual system compositions between woody and herbaceous species, and explore the ecological and evolutionary determinants of these patterns. Specifically, we test the following hypotheses: (1) the geographical patterns in sexual system composition differ between woody and herbaceous species; (2) the relative contributions of different drivers (specifically climate, evolutionary age, and mature plant height) to these patterns differ between woody and herbaceous species. Dioecious species are more common in humid areas and in floras with older and more woody species, whereas hermaphroditic species are more common in temperate arid areas of northwestern China and in floras with younger and more herbaceous species.

Materials and Methods

Sexual System

We compiled a database on the sexual systems of angiosperms in China using published sources: Flora of China (Wu et al., 1994-2013), Flora Republicae Popularis Sinicae (126 issues of 80 volumes), Seeds of Woody Plants in China (Chen and Huang, 2001), and efloras (http://efloras.org/). The Tree of Sex (Ashman et al., 2014), TYR Plant Trait Database (www.try-db.org; Kattge et al., 2011), Botanical Information and Ecology Network (BIEN, Enquist et al., 2016, http://bien.nceas.ucsb.edu/bien/biendata/bien-3/) and journal publications (Sabath et al., 2016; Goldberg et al., 2017) were used to check and supplement the sexual system information in the database. For species with information of sexual systems from multiple sources, conflicting records were checked and removed. In total, we compiled data of sexual systems for 19,780 species from 2,506 genera and 262 families in China (Table 1A).

TABLE 1

Table 1 The correlations between the geographical patterns in proportions of sexual systems among different growth forms.

We classified all species into three categories according to their sexual systems following Cardoso et al. (2018): dioecy (separate male and female individuals), monoecy (having pistils and stamens on separate flowers of the same plant), and hermaphrodites (having both male and female sex organs on the same flower). The category of dioecy includes the androdioecious, gynodioecious, and polygamodioecious species, while monoecy includes the monoecious, andromonoecious, and gynomonoecious species. Considering that a few species vary intraspecifically in their sexual systems (e.g., Dorken et al., 2017) in response to local abiotic or biotic conditions (e.g., climate variables or pollinator densities; Barrett and Harder, 2017), we excluded such species from the final dataset used in this study.

Species Distribution Data

Data on the distributions of plant species were compiled from all national, provincial and local floras, including Flora of China (both Chinese and English versions), Higher Plants of China, The Atlas of Woody Plants in China (Fang et al., 2011; see Wang et al., 2009 for more information), etc. From these resources, we compiled only species distribution records at county-level or finer scales. Species distributions were further supplemented with data from National Specimen Information Infrastructure (NSII, http://www.nsii.org.cn/), and screened. NSII included over fifteen million specimen records that are published recently. In total, our dataset consists the distributions of 19,780 Chinese plant species with data for sexual systems.

Species distributions were rasterized to a grid with a spatial resolution of 100 × 100 km to eliminate the potential bias of area on subsequent analyses. We removed the grid cells with less than 50% of their area along the country borders and coastal areas. In all, 949 grid cells were included in the analyses reported below.

Growth Form

Data on the growth forms of all angiosperm species in China were compiled from the Flora of China (Wu et al., 1994-2013) and Flora Reipublicae Popularis Bulgaricae (Kuzmanov and Kožuharov, 1968). All species were categorized into two growth forms: “woody” and “herbaceous”. Woody species included trees, shrubs and woody lianas, whereas “herbaceous” species included those recorded as herbs, herbaceous lianas, and subshrubs. Species with information on both sexual systems and geographical distributions included 10,200 woody species, 9,479 herbaceous species, and 101 undefined.

Climate Data

Climate can influence flower morphology through its influences on reproductive performance (e.g., pollen production) during the flowering period (Zhang, 2006), and hence could influence sexual systems of plants. Moreover, the reproductive (gametophytic) phase in flowering plants is often highly sensitive to temperature or water stresses. Thus we selected four variables to evaluate the influences of climate, namely, mean temperature of coldest quarter (MTCQ), mean precipitation of warmest quarter (MPWQ), actual evapotranspiration (AET), and aridity index (AI). AET reflects the influences of both water availability and environmental energy, and is significantly associated with ecosystem primary productivity. AET was calculated using the Thornthwaite equation (Thornthwaite and Hare, 1955) from the average monthly temperature and precipitation that were obtained from the WorldClim database (Hijmans et al., 2005). AI reflects the level of aridity and was calculated as the ratio of mean annual precipitation to annual potential evapotranspiration. Higher AI values represent lower aridity. The average and the full range of values for each climatic variable within each grid cell were estimated with the zonal statistics tool in ArcGIS 10.0.

Genus Age

To evaluate the effect of evolutionary time (i.e., genus age) on biogeographical patterns in sexual system compositions, we used the average genus age derived from three phylogenies including Zanne et al. (2014), Lu et al. (2018), and Smith and Brown (2018). This is because (1) different dating methods may lead to different estimations of genus ages; and (2) incomplete phylogenies tend to overestimate the ages of some genera when their close relatives are missed from the phylogenies, and none of these phylogenies is complete at the genus level. For comparison, we also repeated all analyses using the genus ages of each phylogeny separately. It is noteworthy that the geographical patterns in mean genus age per grid cell estimated using the genus ages derived from each phylogeny were consistent with each other (Figure 1A). Moreover, all major findings based on these four age estimations were also highly consistent. Therefore, we included the results based on the average genus ages across the three phylogenies in the main text, and those based on the genus ages derived from each phylogeny in the supplementary materials for comparison.

FIGURE 1

Figure 1 Spatial patterns in proportions of angiosperm species with different sexual systems at 100 × 100 km scale. From top down, the three rows represent the proportions of species with dioecy, monoecy, and hermaphroditism respectively. From left to right, the three columns represent all species, woody species and herbaceous species respectively. Grid cells without data are shown in white. The patterns in the proportions of sexual systems for woody species were updated from Wang et al., 2020.

Mature Plant Height

Plant mature height and longevity of plant species are expected to be positively associated (Marbà et al., 2007; Moles and Leishman, 2008; Moles et al., 2009). Here we used mature height as a proxy for longevity to test the association between sexual system and longevity. Mature height of plant species was obtained from Flora of China (http://frps.eflora.cn/, accessed in November 2013; http://www.efloras.org/flora_page.aspx?flora_id=2, accessed in February 2014). As lower and upper limits of mature height are normally reported for most species in floras, we used the average of lower and upper limits to represent the mature height of species. Species without erect stems (e.g., woody lianas, scandent shrubs, climbers or epiphytes) were excluded from our analyses following Moles et al. (2009). For each growth-form, we averaged the mature height across species within each grid cell to evaluate the effect of plant height on the biogeographical patterns of sexual system compositions (proportions of the three sexual systems per grid cell).

Statistical Analyses

First, based on the data of sexual systems, growth forms, and distributions of all species, we calculated the proportions of species with different sexual systems within each grid cell for all species together and for each growth form separately. We used Pearson correlation coefficients to evaluate the similarity between the geographical patterns in the proportions of different sexual systems of woody and herbaceous species (e.g. woody dioecy vs. herbaceous dioecy, woody monoecy vs. herbaceous monoecy, woody hermaphroditism vs. herbaceous hermaphroditism) and used Dutilleul’s t test to evaluate the significance of these correlations due the influences of spatial autocorrelation on significance tests (using R package of “SpatialPack”; Dutilleul et al., 1993).

We used generalized linear models (GLMs) with quasi-Poisson residuals to evaluate the explanatory power of each predictor on the proportions of sexual systems per grid cell and this was conducted separately for each of the three species groups (i.e. overall, woody, and herbaceous). The climate variable (i.e., MTCQ, MPWQ, AET, and AI), genus age, and plant height were used as predictors and the proportions of each sexual systems per grid cell were used as response variables. The explanatory power of each variable was estimated as the adjusted R²_adj (%) of the GLM model. Modified t tests were used to correct for the effect of spatial autocorrelation on p values (Clifford et al., 1989).

To evaluate the influences of growth forms on the relationships between the proportion of sexual systems and different predictors, we first built spatial linear models (SLM) for the combined dataset of sexual system for both growth forms together using spatial simultaneous autoregressive error (SAR) models. SLMs could account for the effects of residual spatial autocorrelation on the significance tests of regression slopes (Kissling and Carl, 2008). The climate variables (i.e., AET, AI, MTCQ, and MPWQ), genus age, plant height, growth form (i.e., woody vs. herbaceous), and interaction terms between growth form and other predictors were included as the explanatory variables. All the predictors were standardized before fitting the models. Similarly, we also conducted SLMs for each growth form separately. Moran’s I was estimated for the residuals of SLMs, and indicated no spatial autocorrelations.

All analyses were performed in R 3.4.3 (Core Team R, 2017).

Results

Patterns in the Proportions of Sexual Systems for Different Growth Forms

The spatial patterns in the proportions of hermaphroditic and monoecious species per grid cell were consistent between woody and herbaceous growth forms (r = 0.304 and p < 0.05 for proportion of hermaphroditic species; 0.382 and p < 0.05 for proportion of monoecious species; Table 1). For all species and both growth forms, the proportion of hermaphroditic species per grid cell was higher in northwest China than in other regions, while the proportion of monoecious species was higher in eastern China (Figure 1). In contrast, spatial patterns in the proportion of dioecious species were not significantly correlated between woody and herbaceous growth forms (P > 0.05; Table 1). The woody dioecious species were the most frequent in northeast China and medium in southern China, while herbaceous dioecious species were the most frequent in southwest China (Figure 1). Growth form was the strongest predictor for the proportions of all three sexual systems, and the effects of growth form on the proportions of sexual systems were higher than the other predictors by one order of magnitude (Tables 2, 3, and A2; Figure 2).

TABLE 2

Table 2 The slopes of growth forms (woody vs. herbaceous), climate, genus age, plant height, and their interaction on proportions of sexual systems (i.e. dioecy, monoecy, and hermaphroditism) evaluated using spatial linear models (SLMs) with simultaneous autoregressive errors (SAR).

TABLE 3

Table 3 The slopes and significance of different predictors on proportions of sexual systems evaluated using spatial linear models (SLMs) with simultaneous autoregressive errors (SAR).

FIGURE 2

Figure 2 Proportions of species with different sexual systems among different growth forms in China. D, Dioecy; M, monoecy; H, hermaphroditism.

Effects of Climate on Distributions of Proportions of Sexual Systems

The proportions of all three sexual systems were relatively weakly related to climate for herbaceous species than for both all species and woody species (Table 3). The growth form–climate interaction terms were significant in almost all cases (Table 2). This indicates that the relationship between proportions of sexual systems and particular climate variables differed significantly between woody and herbaceous assemblages.

The effects of climate on proportions of sexual systems were significantly greater for woody species than for herbaceous species. For woody species, the proportion of hermaphroditism decreased with AET and MPWQ but increased with AI and MTCQ, while the proportion of dioecious showed almost opposite trends (Table 3). For herbaceous species, all three sexual systems were weakly correlated with climate (R² ranged from 0.217 to 0.474 in Table 3).

Effects of Genus Age on Geographical Distributions of Sexual Systems

Genus age was consistently a significant predictor of sexual system frequencies for woody species, but had only weak effect on herbaceous species. For woody species, the proportion of hermaphroditic was negatively correlated with genus age, while the proportion of dioecious species and monoecious species were positively correlated with genus age (all P < 0.05, Table 3). The explanatory power of genus age was much stronger on the proportions of dioecious and hermaphroditic species than on monoecious species (Table 3).

Effects of Mature Plant Height on the Geographical Frequency of Sexual Systems

Mature plant height per grid cell was a significant predictor for geographical patterns in the proportions of almost all sexual systems examined here (Tables 2 and 3), but its effect was much stronger for woody species than for herbaceous species (Tables 3 and A2). For both growth forms, the proportion of hermaphroditic species decreased. For woody species, the proportion of dioecious species increased but the monoecious decreased with plant height per grid cell, whereas herbaceous species showed opposite trend (Table 3).

Discussion

Because different sexual systems are associated with variation of plant growth form, there are large-differences in the distributions of plant sexual systems between contrasting woody and herbaceous species. Compared with Wang et al. (2020), we advanced our understanding of the spatial distribution of plant sexual systems and the determinants underlying these patterns, especially between growth forms. Growth form significantly differed the relationships between geographic patterns in the proportions of all three sexual systems with different biotic and abiotic factors. Climate, evolutionary age, and mature plant height were all significantly correlates of the distributions of sexual systems of woody species than herbaceous ones. These results are partially consistent with the existing evidence that species diversity of woody plants is more closely related with family age and climate than that of herbaceous plants (Oberle et al., 2009; Hawkins et al., 2011).

Effect of Plant Height on Proportions of Sexual Systems

Plant height was a significant predictor of the proportions of the three sexual systems in almost all cases, which partly confirmed previous hypotheses about the coevolution between sexual systems and growth forms (Renner and Ricklefs, 1995; Vamosi et al., 2003). Both plant height and sexual systems could reflect the adaptation and evolutionary trend of plant life-history strategies to a certain extent, and they also tend to coevolve (Renner, 2014; Wang et al., 2020). Owing to limited mating opportunities, selection has tended to favor dioecy in long-lived woody species which are tall (Renner and Ricklefs, 1995; Obbard et al., 2006). Taller species with height dimorphism and greater dispersal investment may also benefit dioecious species (e.g. Salicaceae) by improving the efficiency in pollen and seed dispersal, which further promote the maintenance and continuation of species populations (Pickup and Barrett, 2012; Thomson et al., 2018). In contrast, hermaphroditism tends to be associated with short-lifespan species with low height because hermaphroditic species can self-pollinate and thus tend to accumulate much more deleterious mutations (Etterson and Mazer, 2016), and which may have higher mortality and extinction risks (Obbard et al., 2006; Cheptou, 2019).

The explanatory power of plant height on the frequency of sexual systems was stronger for woody than for herbaceous species. This may be likely because woody and herbaceous plant species have different adaptive strategies in resource allocation between reproductive and vegetative tissues (Abrahamson and Gadgil, 1973; Bonser and Aarssen, 2003). Woody species may have more resources for the investment to sexual systems and reproduction (Porth and El-Kassaby, 2014), while herbaceous species (e.g. the annual Mercurialis annua) tend to exhibit short life-cycles, fast rates of population but limited resource accumulation (Lanfear et al., 2013).

The Effects of Climate on the Proportions of Sexual Systems

The proportion of hermaphroditic species for all species combined was high in areas with high AI or MTCQ, consistent with Wang et al. (2020). The proportion of dioecious species among all species combined (e.g. Anaphalis) was higher in southern China than in other regions. This finding supports the hypothesis that dioecy is more common in tropical and subtropical floras where climate is warm and humid (Freeman et al., 1976; Bawa, 1980; Bawa et al., 1985; Sakai et al., 1995; Vary et al., 2011). In contrast, the proportion of woody dioecious species (e.g. Rhamnus and Acer) is higher in Northeast China with relatively higher AET and MPWQ compared to other regions in China. These results reveal that the differences in requirements for hydrothermal conditions between woody and herbaceous growth forms may influence their geographical patterns in proportions of sexual systems, consistent with our hypothesis. Generally, woody species often have deep roots, and herbaceous species have relatively shallow roots. Therefore, herbaceous species have lower efficiencies in soil water usage and are more sensitive to reduced water supply compared with woody species (Franco, 2002). Correspondingly, the interaction terms between climate and growth form in the spatial linear models explaining the frequencies of different sexual systems were significant in most cases. The differences in the effects of climate on proportions of sexual systems between woody and herbaceous growth forms may be partly due to the coevolution of sexual systems with growth forms (Renner and Ricklefs, 1995; Vamosi et al., 2003; Case and Barrett, 2004). Moreover, our results further corroborate previous findings that climate may have influenced the spatial pattern of sexual systems via its effect on the composition of plant growth forms (Wang et al., 2020).

We also found that the sexual system composition of woody species was more closely associated with contemporary climate than that of herbaceous species, which is consistent with previous studies at broad scales (Oberle et al., 2009; Wang et al., 2009). Compared with herbaceous species that can be protected from frost by being annual/biennial or by the production of underground buds, woody plant species typically with longer life cycles and large sizes expose their stems and buds (with reproductive organs) in climate. Thus climate has been found to have significant influences on the spatial pattern of woody species (Currie and Paquin, 1987) compared with herbaceous species. These results suggest that the differences in the effects of climate on geographical patterns in proportions of sexual systems for woody and herbaceous species may also be associated with the differences in ecophysiological strategies related to resource acquisition and utilization between growth forms.

Effect of Evolutionary Time on Proportions of Sexual Systems

For both woody and all species, the proportion of dioecious species was positively correlated with mean genus age per grid cell, which support our hypothesis that dioecy tends to be maintained in regions with older floras. Similarly, the proportion of woody monoecious species was also positively correlated with mean genus age. These results partly support and expand the hypothesis that dioecy and monoecy are often found to be significantly associated with each other across phylogenies of woody angiosperms (Renner and Ricklefs, 1995). In contrast, the proportions of sexual systems for herbaceous species were weakly correlated with mean genus age (e.g. Apiaceae, Asteraceae, Cyperaceae), which might be due to the higher population turnover rate and faster adaptation to climate of herbaceous species than those of woody species (e.g. Acer and Salicaceae, Obbard et al., 2006; Smith and Beaulieu, 2009). Moreover, the effect of genus age on the proportions of woody sexual systems was much stronger than on those of herbaceous species. The differences in the effect of evolutionary age on proportions of sexual systems between woody and herbaceous growth forms might reflect their differences in evolutionary rates (i.e. speciation and extinction rates). Molecular studies have revealed that compared with herbaceous species, woody species such as trees tend to have longer generation times and lower evolutionary rate (Smith and Beaulieu, 2009), and thus may have lower chances to transit between sexual systems in response to environmental changes.

Conclusions

Here we mapped the geographical distributions of the frequencies of sexual systems and compared the effects of abiotic and biotic factors on sexual system frequencies for both woody and herbaceous flowering plant species in China. The results suggest that woody and herbaceous species differed significantly in their patterns in the frequencies of sexual systems, and also in the determinants underlying these patterns. The proportions of sexual systems of woody species were more strongly influenced by climate, evolutionary age, and mature height, whereas those of herbaceous species were less influenced by these factors. These findings shed light on previous hypotheses about the association between sexual systems and growth forms, and further demonstrate the necessity to differentiate woody and herbaceous growth forms when investigating the evolution and ecology of sexual systems. It is noteworthy that sexual systems may also vary intraspecifically across populations of a single species (e.g. Mercurialis annua, Dorken et al., 2017). Exploring the intraspecific sexual system variations for species with different growth forms along environmental gradients would further improve our understanding on the drivers of sexual system variations. Our study also suggests that compared with both dioecious and hermaphroditic species, more attention should be paid to monoecious species from both evolutionary and ecological perspectives. This study improves our understanding of the mechanisms underlying geographical patterns of reproductive trait diversity across different growth forms and their responses to climate.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation, to any qualified researcher.

Author Contributions

ZW and YW conceived the idea and design of the study. YW, TL, AL, YLi, and YLiu organized the database. YW and TL performed the statistical analysis. YW wrote the first draft of the manuscript. ZW has comprehensively revised the article. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by the National Key Research Development Program of China (#2017YFA0605101; #2018YFA0606104), the Strategic Priority Research Program of Chinese Academy of Sciences (#XDB31000000), and National Natural Science Foundation of China (#31901216, #31988102, #31911530102).

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.2020.01222/full#supplementary-material

References

Aarssen, L. W. (2000). Why are most selfers annuals? A new hypothesis for the fitness benefit of selfing. Oikos 89 (3), 606–612. doi: 10.1034/j.1600-0706.2000.890321.x

CrossRef Full Text | Google Scholar

Abrahamson, W. A., Gadgil, M. (1973). Growth Form and Reproductive Effort in Goldenrods (Solidago, Compositae). Am. Nat. 107 (957), 651–661. doi: 10.1086/282864

CrossRef Full Text | Google Scholar

Ashman, T. L., Bachtrog, D., Blackmon, H., Goldberg, E., Hahn, M., Kirkpatrick, M., et al. (2014). The Tree of Sex Consortium. Tree of Sex: A database of sexual systems. Sci. Data 1, 140015. doi: 10.1038/sdata.2014.15

PubMed Abstract | CrossRef Full Text | Google Scholar

Ashman, T. L. (2006). “The evolution of separate sexes: A focus on the ecological context,” in Ecology and Evolution of Flowers. Eds. Harder, L. D., Barrett, S. C. H. (Oxford: Oxford University Press).

Google Scholar

Barrett, S. C. H. (2013). The evolution of plant reproductive systems: how often are transitions irreversible? Proc. Royal Soc. B: Biol. Sci. 280 (1765), 20130913. doi: 10.1098/rspb.2013.0913

CrossRef Full Text | Google Scholar

Barrett, S. C. H., Harder, L. D. (2017). The Ecology of Mating and Its Evolutionary Consequences in Seed Plants. Annu. Rev. Ecol. Evol. Syst. 48, 135–157. doi: 10.1146/annurev-ecolsys-110316-023021

CrossRef Full Text | Google Scholar

Bawa, K. S., Bullock, S. H., Perry, D. R., Coville, R. E., Grayum, M. H. (1985). Reproduction biology of tropical lowland rain forest trees. II. Pollination systems. Am. J. Bot. 72, 346–356. doi: 10.1002/j.1537-2197.1985.tb05358.x

CrossRef Full Text | Google Scholar

Bawa, K. S. (1980). Evolution of dioecy in flowering plants. Annu. Rev. Ecol. Evol. Syst. 11, 15–39. doi: 10.1146/annurev.es.11.110180.000311

CrossRef Full Text | Google Scholar

Bonser, S. P., Aarssen, L. W. (2003). Allometry and development in herbaceous plants: functional responses of meristem allocation to light and nutrient availability. Am. J. Bot. 90 (3), 404–412. doi: 10.3732/ajb.90.3.404

PubMed Abstract | CrossRef Full Text | Google Scholar

Cardoso, J. C. F., Viana, M. L., Matias, R., Furtado, M., Caetano, A. P. S., Consolaro, H., et al. (2018). Towards a unified terminology for angiosperm reproductive systems. Acta Bot. Bras. 32 (3), 329–348. doi: 10.1590/0102-33062018abb0124

CrossRef Full Text | Google Scholar

Case, A. L., Barrett, S. C. H. (2004). Environmental stress and the evolution of dioecy: Wurmbea dioica (Colchicaceae) in Western Australia. Evol. Ecol. 18, 145–164. doi: 10.1023/B:EVEC.0000021152.34483.77

CrossRef Full Text | Google Scholar

Charlesworth, D., Wright, S., II (2001). Breeding systems and genome evolution. Curr. Opin. Genet. Dev. 11 (6), 685–690. doi: 10.1016/S0959-437X(00)00254-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Y. S., Huang, P. C. (2001). Seeds of woody plants in China (Beijing: China Woody Plant Press).

Google Scholar

Chen, X. S., Li, Q.-J. (2008). Patterns of plant sexual systems in subtropical evergreen broad-leaved forests in Ailao Mountains, SW China. J. Plant Ecol. 1 (3), 179–185. doi: 10.1093/jpe/rtn019

CrossRef Full Text | Google Scholar

Cheptou, P. O. (2019). Does the evolution of self-fertilization rescue populations or increase the risk of extinction? Ann. Bot. 123 (2), 337–345. doi: 10.1093/aob/mcy144

PubMed Abstract | CrossRef Full Text | Google Scholar

Clifford, P., Richardson, S., Hemon, D. (1989). Assessing the significance of the correlation between two spatial processes. Biometrics 45 (1), 123–144. doi: 10.2307/2532039

PubMed Abstract | CrossRef Full Text | Google Scholar

Core Team, R. (2017). R: A Language and Environment for Statistical Computing. (Vienna: R Foundation for Statistical Computing). https://www.R-project.org/.

Google Scholar

Currie, D. J., Paquin, V. (1987). Large-scale biogeographical patterns of species richness of trees. Nature 329, 326–327. doi: 10.1038/329326a0

CrossRef Full Text | Google Scholar

Curtis, C. A., Bradley, B. A. (2016). Plant Distribution Data Show Broader Climatic Limits than Expert-Based Climatic Tolerance Estimates. PloS One 11 (11), e0166407. doi: 10.1371/journal.pone.0166407

PubMed Abstract | CrossRef Full Text | Google Scholar

Dorken, M. E., Freckleton, R. P., Pannell, J. R. (2017). Small-scale and regional spatial dynamics of an annual plant with contrasting sexual systems. J. Ecol. 105 (4), 1044–1057. doi: 10.1111/1365-2745.12719

CrossRef Full Text | Google Scholar

Dutilleul, P., Clifford, P., Richardson, S., Hemon, D. (1993). Modifying the t test for assessing the correlation between two spatial processes. Biometrics 49 (1), 305–314. doi: 10.2307/2532625

CrossRef Full Text | Google Scholar

Eckert, C. G., Kalisz, S., Geber, M. A., Sargent, R., Elle, E., Cheptou, P. O., et al. (2010). Plant mating systems in a changing world. Trends Ecol. Evol. 25 (1), 35–43. doi: 10.1016/j.tree.2009.06.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Enquist, B. J., Condit, R., Peet, R. K., Schildhauer, M., Thiers, B. M. (2016). Cyber infrastructure for an integrated botanical information network to investigate the ecological impacts of global climate change on plant biodiversity. PeerJ 4, e2615v2. doi: 10.7287/peerj.preprints.2615v2

CrossRef Full Text | Google Scholar

Etterson, J. R., Mazer, S. J. (2016). How climate change affects plants' sex lives. Science 353 (6294), 32–33. doi: 10.1126/science.aag1624

CrossRef Full Text | Google Scholar

Franco, A. C (2002). “Ecophysiology of woody plants: Ecology and Natural History of a Neotropical Savanna,” in Thecerrados of Brazil. Eds. Oliveira, P. S., Marquis, R. J. (New York: Columbia University Press), 178–197. doi: 10.7312/oliv12042-009

CrossRef Full Text | Google Scholar

Fang, J. Y., Wang, Z. H., Tang, Z. Y. (Eds.) (2011). “Atlas of Woody Plants in China,” in Distribution and Climate (Berlin Heidelberg: Springer-Verlag).

Google Scholar

Fitzjohn, R. G., Pennell, M. W., Zanne, A. E., Stevens, P. F., Tank, D. C., Cornwell, W. K. (2014). How much of the world is woody? J. Ecol. 102 (5), 1266–1272. doi: 10.1111/1365-2745.12260

CrossRef Full Text | Google Scholar

Freeman, D. C., Klikoff, L. G., Harper, K. T. (1976). Dirrerential resource utilization by the sexes of dioecious plants. Science 193, 597–599. doi: 10.1126/science.193.4253.597

PubMed Abstract | CrossRef Full Text | Google Scholar

Friedman, J., Barrett, S. C. H. (2009). Wind of change: new insights on the ecology and evolution of pollination and mating in wind-pollinated plants. Ann. Bot. 103 (9), 1515–1527. doi: 10.1093/aob/mcp035

PubMed Abstract | CrossRef Full Text | Google Scholar

Gamble, D. E., Bontrager, M., Angert, A. L. (2018). Floral trait variation and links to climate in the mixed-mating annual Clarkia pulchella. Botany 96 (7), 425–435. doi: 10.1139/cjb-2017-0234

CrossRef Full Text | Google Scholar

Goldberg, E. E., Otto, S. P., Vamosi, J. C., Mayrose, I., Sabath, N., Ming, R., et al. (2017). Macroevolutionary synthesis of flowering plant sexual systems. Evolution 71 (4), 898–912. doi: 10.1111/evo.13181

PubMed Abstract | CrossRef Full Text | Google Scholar

Hawkins, B. A., Rodríguez, M. Á., Weller, A. G. (2011). Global angiosperm family richness revisited: linking ecology and evolution to climate. J. Biogeogr. 38 (7), 1253–1266. doi: 10.1111/j.1365-2699.2011.02490.x

CrossRef Full Text | Google Scholar

Hedhly, A., Hormaza, J., II, Herrero, M. (2009). Global warming and sexual plant reproduction. Trends Plant Sci. 14 (1), 30–36. doi: 10.1016/j.tplants.2008.11.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Hijmans, R. J., Cameron, S. E., Parra, J. L., Jones, P. G., Jarvis, A. (2005). Very high-resolution interpolated climate surfaces for global land areas. Int. J. Climatol. 25 (15), 1965–1978. doi: 10.1002/joc.1276

CrossRef Full Text | Google Scholar

Hultine, K. R., Grady, K. C., Wood, T. E., Shuster, S. M., Stella, J. C., Whitham, T. G. (2016). Climate change perils for dioecious plant species. Nat. Plants 2, 16109. doi: 10.1038/nplants.2016.109

PubMed Abstract | CrossRef Full Text | Google Scholar

Käfer, J., de Boer, H. J., Mousset, S., Kool, A., Dufay, M., Marais, G. A. B. (2014). Dioecy is associated with higher diversification rates in flowering plants. J. Evol. Biol. 27 (7), 1478–1490. doi: 10.1111/jeb.12385

PubMed Abstract | CrossRef Full Text | Google Scholar

Kattge, J., Díaz, S., Lavorel, S., Prentice, I. C., Leadley, P., Bönisch, G., et al. (2011). TRY – a global database of plant traits. Glob. Chang. Biol. 17 (9), 2905–2935. doi: 10.1111/j.1365-2486.2011.02451.x

CrossRef Full Text | Google Scholar

Kay, K. M., Whittall, J. B., Hodges, S. A. (2006). A survey of nuclear ribosomal internal transcribed spacer substitution rates across angiosperms: an approximate molecular clock with life history effects. BMC Evol. Biol. 6, 36. doi: 10.1186/1471-2148-6-36

PubMed Abstract | CrossRef Full Text | Google Scholar

Kissling, W. D., Carl, G. (2008). Spatial autocorrelation and the selection of simultaneous autoregressive models. Glob. Ecol. Biogeogr. 17 (1), 59–71. doi: 10.1111/j.1466-8238.2007.00334.x

CrossRef Full Text | Google Scholar

Klekowski, E. J., Godfrey, P. J. (1989). Aging and mutation in plants. Nature 340, 389–391. doi: 10.1038/340389a0

CrossRef Full Text | Google Scholar

Kuzmanov, B., Kožuharov, S. (1968). Flora Reipublicae Popularis Bulgaricae. Taxon. 17 (2), 205–210.

Google Scholar

Lanfear, R., Simon, Y. W., Ho, S. Y. W., Davies, T. J., Moles, A. T., Aarssen, L., et al. (2013). Taller plants have lower rates of molecular evolution. Nat. Commun. 4, 1879. doi: 10.1038/ncomms2836

PubMed Abstract | CrossRef Full Text | Google Scholar

Laroche, J., Li, P., Maggia, L., Bousquet, J. (1997). Molecular evolution of angiosperm mitochondrial introns and exons. PNAS 94 (11), 5722–5727. doi: 10.1073/pnas.94.11.5722

PubMed Abstract | CrossRef Full Text | Google Scholar

Lloyd, D. G. (1980). Sexual strategies in plants. I. An hypothesis of serial adjustment of maternal investment during one reproductive session. New Phytol. 86 (1), 69–79. doi: 10.1111/j.1469-8137.1980.tb00781.x

CrossRef Full Text | Google Scholar

Lu, L. M., Mao, L. F., Yang, T., Ye, J. F., Liu, B., Li, H. L., et al. (2018). Evolutionary history of the angiosperm flora of China. Nature 554, 234–238. doi: 10.1038/nature25485

PubMed Abstract | CrossRef Full Text | Google Scholar

Marbà, N., Duarte, M. C., Agustí, S. (2007). Allometric scaling of plant life history. PNAS 104 (40), 15777–15780. doi: 10.1073/pnas.0703476104

PubMed Abstract | CrossRef Full Text | Google Scholar

Mitchell, C. H., Diggle, P. K. (2005). The evolution of unisexual flowers: morphological and functional convergence results from diverse developmental transitions. Am. J. Bot. 92 (7), 1068–1076. doi: 10.3732/ajb.92.7.1068

PubMed Abstract | CrossRef Full Text | Google Scholar

Moeller, D. A., Runquist, B., D, R., Moe, A. M., Geber, M. A., Goodwillie, C., et al. (2017). Global biogeography of mating system variation in seed plants. Ecol. Lett. 20 (3), 375–384. doi: 10.1111/ele.12738

PubMed Abstract | CrossRef Full Text | Google Scholar

Moles, A. T., Leishman, M. R. (2008). “The seedling as part of a plant’s life history strategy,” in Seedling Ecology and Evolution. Eds. Leck, M. A., Parker, V. T., Simpson, R. L. (Cambridge: Cambridge University Press).

Google Scholar

Moles, A. T., Warton, D., II, Warman, L., Swenson, N. G., Laffan, S. W., Zanne, A. E., et al. (2009). Global patterns in plant height. J. Ecol. 97 (5), 923–932. doi: 10.1111/j.1365-2745.2009.01526.x

CrossRef Full Text | Google Scholar

Morgan, M. T., Schoen, D. J., Bataillon, T. M. (1997). The evolution of self-fertilization in perennials. Am. Nat. 150 (5), 618–638. doi: 10.1086/286085

PubMed Abstract | CrossRef Full Text | Google Scholar

Nazareno, A. G., Alzate–Marin, A. L., Pereira, R. A. S. (2013). Dioecy, more than monoecy, affects plant spatial genetic structure: the case study of Ficus. Ecol. Evol. 3 (10), 3495–3508.

PubMed Abstract | Google Scholar

Obbard, D. J., Harris, S. A., Pannell, J. R. (2006). Sexual systems and population genetic structure in an annual plant: testing the metapopulation model. Am. Nat. 167, 354–366. doi: 10.1086/499546

PubMed Abstract | CrossRef Full Text | Google Scholar

Oberle, B., Grace, J. B., Chase, J. M. (2009). Beneath the veil: plant growth form influences the strength of species richness–productivity relationships in forests. Global Ecol. Biogeography 18, 416–425. doi: 10.1111/j.1466-8238.2009.00457.x

CrossRef Full Text | Google Scholar

Pickup, M., Barrett, S. C. H. (2012). Reversal of height dimorphism promotes pollen and seed dispersal in a wind-pollinated dioecious plant. Biol. Lett. 8 (2), 245–248. doi: 10.1098/rsbl.2011.0950

PubMed Abstract | CrossRef Full Text | Google Scholar

Porth, I., El-Kassaby, Y. A. (2014). Assessment of the Genetic Diversity in Forest Tree Populations Using Molecular Markers. Diversity 6 (2), 283–295. doi: 10.3390/d6020283

CrossRef Full Text | Google Scholar

Queenborough, S. A., Mazer, S. J., Vamosi, S. M., Garwood, N. C., Valencia, R., Freckleton, R. P. (2009). Seed mass, abundance and breeding system among tropical forest species: Do dioecious species exhibit compensatory reproduction or abundances? J. Ecol. 97 (3), 555–566. doi: 10.1111/j.1365-2745.2009.01485.x

CrossRef Full Text | Google Scholar

Renner, S. S., Ricklefs, R. E. (1995). Dioecy and its correlates in the flowering plants. Am. J. Bot. 82 (5), 596–606. doi: 10.1002/j.1537-2197.1995.tb11504.x

CrossRef Full Text | Google Scholar

Renner, S. S. (2014). The relative and absolute frequencies of angiosperm sexual systems: dioecy, monoecy, gynodioecy, and an updated online database. Am. J. Bot. 101 (10), 1588–1596. doi: 10.3732/ajb.1400196

PubMed Abstract | CrossRef Full Text | Google Scholar

Rosche, C., Schrieber, K., Lachmuth, S., Durka, W., Hirsch, H., Wagner, V., et al. (2018). Sex ratio rather than population size affects genetic diversity in Antennaria dioica. Plant Biol. 20, 789–796. doi: 10.1111/plb.12716

CrossRef Full Text | Google Scholar

Sabath, N., Goldberg, E. E., Glick, L., Einhorn, M., Ashman, T.-L., Ming, R., et al. (2016). Dioecy does not consistently accelerate or slow lineage diversification across multiple genera of angiosperms. New Phytol. 209 (3), 1290–1300. doi: 10.1111/nph.13696

PubMed Abstract | CrossRef Full Text | Google Scholar

Sakai, A. K., Wagner, W. L., Ferguson, D. M., Herbst, D. R. (1995). Biogeographical and ecological correlates of dioecy in the Hawaiian flora. Ecology 76 (8), 2530–2543. doi: 10.2307/2265826

CrossRef Full Text | Google Scholar

Sakai, A. K., Wagner, W. L., Ferguson, D. M., Herbst, D. R. (1995). Origins of Dioecy in the Hawaiian Flora. Ecology 76 (8), 2517–2529. doi: 10.2307/2265825

CrossRef Full Text | Google Scholar

Senarath, W. T. P. S. K. (2008). Dioecy and monoecy in the flora of Sri Lanka and their evolutionary correlations to endemism, growth form, fruit type, seed number and flower size. Bangl. J. Plant Taxon. 15 (1), 13–19. doi: 10.3329/bjpt.v15i1.907

CrossRef Full Text | Google Scholar

Šímová, I., Violle, C., Svenning, J. C., Kattge, J., Engemann, K., Sandel, B., et al. (2018). Spatial patterns and climate relationships of major plant traits in the New World differ between woody and herbaceous species. J. Biogeogr. 45 (4), 895–916. doi: 10.1111/jbi.13171

CrossRef Full Text | Google Scholar

Smith, S. A., Beaulieu, J. M. (2009). Life-history influences rates of climatic niche evolution in flowering plants. P. R. Soc Lon. B. Bio. 276 (1677), 4345–4352. doi: 10.1098/rspb.2009.1176

CrossRef Full Text | Google Scholar

Smith, S. A., Brown, J. W. (2018). Constructing a broadly inclusive seed plant phylogeny. Am. J. Bot. 105, 302–314. doi: 10.1002/ajb2.1019

PubMed Abstract | CrossRef Full Text | Google Scholar

Soza, V. L., Brunet, J., Liston, A., Smith, P. S., Di Stilio, V. S. (2012). Phylogenetic insights into the correlates of dioecy in meadow-rues (Thalictrum, Ranunculaceae). Mol. Phylogenet. Evol. 63 (1), 180–192. doi: 10.1016/j.ympev.2012.01.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Thomson, F. J., Letten, A. D., Tamme, R., Edwards, W., Moles, A. T. (2018). Can dispersal investment explain why tall plant species achieve longer dispersal distances than short plant species? New Phytol. 217 (1), 407–415. doi: 10.1111/nph.14735

PubMed Abstract | CrossRef Full Text | Google Scholar

Thornthwaite, C. W., Hare, F. K. (1955). Climatic classification in forest. Unasylva. 9, 51–59.

Google Scholar

Vamosi, J. C., Vamosi, S. M. (2004). The role of diversification in causing the correlations of dioecy. Evolution 58 (4), 723–731. doi: 10.1111/j.0014-3820.2004.tb00405.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Vamosi, J. C., Otto, S. P., Barrett, S. C. (2003). Phylogenetic analysis of the ecological correlates of dioecy in angiosperms. J. Evol. Biol. 16 (5), 1006–1018. doi: 10.1046/j.1420-9101.2003.00559.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Vary, L. B., Gillen, D. L., Randrianjanahary, M., Lowry, P. P., II, Sakai, A. K., Weller, S. G. (2011). Dioecy, monoecy, and their ecological correlates in the littoral forest of Madagascar. Biotropica 43 (5), 582–590. doi: 10.1111/j.1744-7429.2010.00742.x

CrossRef Full Text | Google Scholar

Wang, Z. H., Brown, J. H., Tang, Z. Y., Fang, J. Y. (2009). Temperature dependence, spatial scale, and tree species diversity in eastern Asia and North America. PNAS 106 (32), 13388–13392. doi: 10.1073/pnas.0905030106

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Z. H., Fang, J. Y., Tang, Z. Y., Lin, X. (2011). Patterns, terminants and models of woody plant diversity in China. P. R. Soc Lon. B. Bio. 278 (1715), 2122–2132. doi: 10.1098/rspb.2010.1897

CrossRef Full Text | Google Scholar

Wang, Y. Y., Freckleton, R. P., Wang, B. J., Kuang, X., Yuan, Z. Q., Lin, F., et al. (2018). The role of breeding system in community dynamics: Growth and mortality in forests of different successional stages. Ecol. Evol. 8 (15), 7285–7296. doi: 10.1002/ece3.4190

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y. Y., Lyu, T., Shrestha, N., Lyu, L. S., Li, Y. Q., Schmid, B., et al. (2020). Drivers of large-scale geographical variation in sexual systems of woody plants. Glob. Ecol. Biogeogr. 29 (3), 546–557. doi: 10.1111/geb.13052

CrossRef Full Text | Google Scholar

Ward, M., Dick, C. W., Gribel, R., Lowe, A. J. (2005). To self, or not to self… A review of outcrossing and pollen-mediated gene flow in neotropical trees. Heredity 95 (4), 246–254. doi: 10.1038/sj.hdy.6800712

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, Z., Raven, P., Hong, D. Y. (1994–2013). Flora of China. (Beijing and St Louis: Science Press and Missouri Botanical Garden Press). Online at eFloras.org.

Google Scholar

Zanne, A., Tank, D., Cornwell, W., Eastman, J. M., Smith, S. A., FitzJohn, R. G., et al. (2014). Three keys to the radiation of angiosperms into freezing environments. Nature 506, 89–92. doi: 10.1038/nature12872

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, D. Y. (2006). “Evolutionary stable reproductive investment and sex allocation in plants,” in Ecology and Evolution of Flowers. Eds. Harder, L., Barrett, S. C. H. (Oxford: Oxford University Press), 41–60. doi: 10.1034/j.1600-0706.2001.950118.x

CrossRef Full Text | Google Scholar

Wu, Z. Y., Raven, P. H., Hong, D. Y. (Eds.) (2007). “Flora of China,” in Clusiaceae through Araliaceae, vol. 13 (Beijing: Science Press, and St. Louis: Missouri Botanical Garden Press).

Google Scholar