ORIGINAL RESEARCH article
Spatially Structured Environmental Variation Plays a Prominent Role on the Biodiversity of Freshwater Macrophytes Across China
- 1Freshwater Aquaculture Collaborative Innovation Center of Hubei Province, Hubei Provincial Engineering Laboratory for Pond Aquaculture, College of Fisheries, Huazhong Agricultural University, Wuhan, China
- 2Arctic Research Center, Hokkaido University, Sapporo, Japan
- 3Global Station for Arctic Research, Global Institution for Collaborative Research and Education, Hokkaido University, Sapporo, Japan
- 4Division of Environmental Science Development, Graduate School of Environmental Science, Hokkaido University, Sapporo, Japan
- 5Laboratoire Evolution et Diversité Biologique (EDB), Université de Toulouse, CNRS, ENFA, UPS, Toulouse, France
- 6Donghu Experimental Station of Lake Ecosystems, State Key Laboratory of Freshwater Ecology and Biotechnology of China, Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan, China
Different non-mutually exclusive mechanisms interactively shape large-scale diversity patterns. However, our understanding of multi-faceted diversity and their determinants in aquatic ecosystems is far from complete compared to terrestrial ones. Here, we use variation partitioning based on redundancy analysis to analyze the relative contribution of environmental and spatial variables to the patterns of phylogenetic, taxonomic, and functional diversity in macrophyte assemblages across 214 Chinese watersheds. We found extremely high spatial congruence among most aspects of biodiversity, with some important exceptions. We then used variation partitioning to estimate the proportions of variation in macrophyte biodiversity explained by environmental and spatial variables. All diversity facets were optimally explained by spatially structured environmental variables, not the pure environment effect, implying that macrophyte are taxonomically, phylogenetically, and functionally clustered in space, which might be the result of the interaction of environmental and/or evolutionary drives. We demonstrate that macrophytes might face extensive dispersal limitations across watersheds such as topography and habitat fragmentation and availability.
Macrophytes support many structural and functional aspects of freshwater ecosystems and their ecological properties have been extensively studied over (Westlake, 1975; Jeppesen et al., 1998; Wiegleb et al., 2015). Highly productive, macrophyte communities have a key role in carbon and nutrient fluxes, serving as sinks for organic material and sources of nutrients to the water (Jeppesen et al., 1998; Chambers et al., 2008; Hansson et al., 2010). They provide structurally complex habitat for a large diversity of organisms such as macroinvertebrate, fish and birds (Westlake, 1975; Jeppesen et al., 1998; Hansson et al., 2010). Consequently, a significant change in the structure and composition of macrophyte communities can have important knock on effects on freshwater ecosystems with important management conservation implications (Hellsten and Riihimäki, 1996; Hansson et al., 2010; Wiegleb et al., 2016; Zhang et al., 2019). The development of local macrophyte assemblages strongly depends on a variety of abiotic and biotic factors, e.g., nutrient concentrations, flow velocity, light condition, and trophic interaction (Westlake, 1975; Dodkins et al., 2005; Zhang et al., 2019). However, consensus regarding the generality of large-scale processes driving spatial variation in biodiversity of macrophytes remains elusive (Chappuis et al., 2014; Wiegleb et al., 2015).
Most macrophyte species are regarded as cosmopolitan, their broad distribution traditionally explained by common life history traits such as long-distance hydrochory, anemochory, and zoochory seed dispersal (Santamaria, 2002; Chambers et al., 2008). Such strong dispersal capacity has facilitated the evolution of their plastic response and ecological tolerances to local environmental change. For example, many macrophytes are very resilient because of their fast asexual reproduction abilities such as clonal growth and abundant propagules (Chambers et al., 2008). Consequently, under natural conditions, regional-scale taxonomic richness is generally high and relatively stable across space (Alahuhta et al., 2018). However, high environmental heterogeneity may promote species selection (environmental filtering) and persistence (niche diversification), increasing the spatial turnover of species at landscape scales (Xu et al., 2015; Alahuhta et al., 2018). Therefore, regional spatial diversity patterns in aquatic taxa may be governed by the interplay between their dispersal capacity and the spatiotemporal heterogeneity (Shmida and Wilson, 1985; Shurin, 2007; Gianuca et al., 2017; Cai et al., in press).
Three aspects of species diversity patterns (i.e., taxonomic, phylogenetic, and functional diversity) are the main focus of macroecological studies. Taxonomic diversity, the most straightforward and commonly used measurement of biodiversity in broad-scale studies (Tisseuil et al., 2013; Heino et al., 2015; Cai et al., 2018), treats all species as functionally equivalent and phylogenetically independent. However, phylogenetic diversity and functional trait variation can exert a much stronger control on biodiversity effects on ecosystem functions, such as production or nutrient cycling, than taxonomic diversity (Cadotte et al., 2011; Flynn et al., 2011). Phylogenetic diversity, incorporating evolutionary relationships between species, provides also promising way of interpreting the role of biogeographic history in community ecology (Webb et al., 2002). Therefore, studies considering all these complementary facets can provide a more complete understanding of the mechanistic links between ecological processes and evolutionary history in shaping biodiversity patterns and the provision of ecosystem services (Devictor et al., 2010; Tucker and Cadotte, 2013; Pardo et al., 2017). Here, we assess the relative contribution of spatial and environmental variables to multifaceted biodiversity patterns in macrophyte assemblages across 214 tributary drainage basins (hereafter watersheds) covering the whole China. While the existence of spatial congruence among patterns in diversity facets is an open debate, our previous work in the Yangtze River has suggested high spatial congruence of macrophyte taxonomic and functional diversity at the catchment level (Zhang et al., 2018). Whereas these patterns hold at larger, regional scales across catchments remains to be tested. On the other hand, environmental variables are themselves spatially structured (Legendre and Legendre, 2012), making it possible for spatially structured environmental variation to drive the observed variability in spatial patterns within and among diversity facets.
Materials and Methods
Data Acquisition and Key Definitions
Data on freshwater macrophyte species compiled from published sources (see below) was grouped at the watershed-scale (214 watersheds) across China (National Remote Sensing Center of China, as delimited by the National Council of China under the National Water Resources Strategic Plan1. These watersheds represent subdivisions of main river basins based on different ecohydrological criteria (e.g., river order, landscape, climate), and provide a spatial basis for the development, utilization, conservation, and management of hydrological resources in China (Cai et al., 2018). Pooling data sets into meaningful, large spatial working units such as watersheds (Pool et al., 2014) or ecoregions (Veech and Crist, 2007) is a practical, compromise solution frequently used for the analysis of macroscale diversity patterns with spatially sparse data sets (Cai et al., in press).
Although different definitions of “macrophyte” are available in the literature, we follow Cook (1974) by considering any aquatic plant that is visible to the naked eye including all higher aquatic plants, vascular cryptogams, bryophytes, and groups of algae that can be seen to be dominated by a single species. Based on this definition, we made a detailed literature review of macrophyte species in China from published (1960–2010) records related to lakes, rivers, and seasonal agricultural ponds. Documented sources included research articles and monographs together with the Scientific Database of China Plant Species2, the Database of Invasive Alien species in China3, Chinese Species Information System4, and gray research reports. This exhaustive literature review provided information for a total of 992 aquatic plant species. We then prepared a data matrix covering taxonomic information and functional traits of the species. To guarantee consistency across the data set we used five quality-control rules: (1) non-macrophyte species were filtered according to the Cook’s definition for macrophytes (Cook, 1974) and the records in Flora of China, (2) scientific names were standardized and synonyms were removed on the basis of the Chinese Virtual Herbarium4, (3) varieties were treated as a single species, (4) the distribution traits of the species were corrected according to the Flora of China, and (5) non-freshwater species were excluded. The application of these rules resulted in a total of 469 species from 214 watersheds retained for analysis, including 93 submerged species, 40 floating-leaved species, 25 free-floating species, and 311 emergent species.
Measurements of Multiple Facets of Biodiversity
We defined taxonomic richness (TRic) as the number of species recorded in each watershed over the study period. Given the lack of sufficient and consistent phylogeny information about all included macrophytes in our data set, phylogeny diversity was assessed using taxonomic hierarchies as a proxy for phylogenetic relationships (Schweiger et al., 2008; Heino et al., 2015). Taxonomic distinctness (TDis) measures the mean taxonomic (i.e., phylogenetic) distances between species (Clarke and Warwick, 1998), which was calculated by giving equal branch lengths and six supra-species taxonomic levels (i.e., genus, family, order, subclass, class, and phylum). Hence, watershed with low TDis value indicates low phylogenetic diversity, and vice versa. TDis was calculated using the functions “taxondive” and “taxa2dist” in the R package “vegan” (Oksanen et al., 2016).
To calculate functional distance between macrophyte species, we gathered information on four categorical functional traits from the Flora of China, namely life form (i.e., free-floating, emergent, floating-leaved, submerged), life cycle (i.e., annual, perennial), morphology (i.e., turion, stem, rosette, leafy), sexual propagation (monoecism, dioecy), and species’ mean adult weight (Zhang et al., 2018). Morphology was defined qualitatively based on the plant description. Leafy plants typically have more lamina, often the parts concentrating the majority of photosynthesis; stem plants are those with stem and easy to propagate due to broken branches; rosette plants have a shortened stem axis and relatively large projection area that facilitates light competition; turion plants produce winter/overwintering buds as dormant storage organs in response to unfavorable ecological conditions (Cook, 1974; Adamec, 2018). Although we acknowledge that quantitative morphological traits such as shoot height, stem diameter, specific leaf area, or leaf dry mass content, available from some local studies (Fu et al., 2014, 2018), would provide a more precise assessment, these were not available given the nature of our data set and the large scale of our study area.
Four multidimensional functional diversity indices (i.e., functional richness, functional evenness, functional divergence, and functional dispersion) were computed for each assemblage. Functional richness (FRic) describes the convex hull volume filled by a community in the multidimensional functional trait space and used as a measure of the functional richness. Functional evenness (FEve) describes the evenness of the distribution of species in a community over the functional trait space by using the minimum spanning tree. Feve quantifies the regularity with which the functional space is filled by species. Functional divergence (FDiv) describes how species distribute within the volume of functional trait space. For presence/absence data, FDiv is the highest when all the species are on the convex hull and at equal distance to its center of gravity (i.e., if the center of gravity of the convex hull is also a center of symmetry of the functional space). Finally, functional dispersion (FDis) describes the mean distance of each individual species to the centroid of all other species in the assemblage (Anderson et al., 2006), which was calculated using the function “dbFD” in the R package “FD” (Laliberté et al., 2014).
We include several major environmental factors (Supplementary Table S1) used in previous similar macroecological studies (Whittaker et al., 2001; Tisseuil et al., 2013; Cai et al., 2018), namely mean annual precipitation (MAP), mean annual temperature (MAT), solar radiation (SOLAR), and total annual runoff (RUN), as surrogates for the energy input in each watershed, together with total surface area of each watershed (AREA), and spatial variation of altitude, MAT, MAP, SOLAR, and the Shannon diversity index based on the proportions of land cover classes (forest, grass, farm, urban, water, and desert) (i.e., ALTVAR, MATVAR, MAPVAR, SOLARVAR, and LANDVAR) within each watershed representing important factors shaping biodiversity through increasing habitat diversity and availability (Whittaker et al., 2001; Tisseuil et al., 2013). ALTVAR is used as a proxy for topographic heterogeneity, calculated as the range between the maximum and minimum altitude for each watershed. The selected factors should cover the environmental drivers of macrogeographic distributions of species at the watershed scale.
All environmental data were extracted from open-access databases such as Data Sharing Infrastructure of Earth System Science, National Science and Technology Infrastructure Center5 and the Data Center of Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences. All parameters were calculated as the average of all cell values with centroids falling within each watershed at a spatial resolution of 0.5° × 0.5° over the study period (1960–2010). Prior to the analyses, variables were logarithm or square root-transformed to improve normalization when necessary and standardized to have a mean of 0 and a variance of 1. We also computed the values of variance inflation factor (VIF) for each predictor variable before analyses to assess collinearity (Cai et al., 2018).
Spatial structure in natural communities can be simultaneously generated by spatial autocorrelation in species assemblages and forcing (explanatory) variables, such as environmental and biotic controls or life history events (Borcard et al., 2004). In order to account for the contribution of the spatial structure of watershed configuration to observed variability in diversity facets, we used multiscale principal coordinates of neighbor matrices (PCNM; Borcard and Legendre, 2002) applied to a geographic distance matrix computed from the watershed centroid coordinates. In essence, distance between watersheds was first represented as a Euclidean distance matrix calculated from the watershed centroid coordinates. Principal coordinates analysis (PCoA) was then conducted on this matrix after truncation using a threshold distance equivalent to the minimum spanning tree of the distance matrix, defining what is considered to be “large” distances.
The resulting eigenvectors with positive eigenvalues (distance-based Moran’s Eigenvector Maps; dbMEMs), provide a spectral decomposition of any possible spatial relationships between the watersheds (Borcard and Legendre, 2002). That is, each eigenvector captures the dominant spatial structures, i.e., low-order eigenvectors (dbMEM1-3 in Supplementary Figure S1), associated with large eigenvalues, represent country-scale groupings whereas high-order eigenvectors (dbMEM69-71 in Supplementary Figure S1), associated with small eigenvalues, represent more regional-scale groupings. Following a constrained analysis, those positive eigenvectors found to significantly explain variation in macrophyte diversity were selected as independent explanatory variables for further analysis. dbMEMs were calculated using the function “eigenmap” with default values in the package “codep” in R (Dray et al., 2012).
Spatial congruence of watersheds with high diversity was assessed as the amount of concordance between the top 10% of watersheds (i.e., highest diversity values) for all six diversity indices (Pool et al., 2014). That is, the 21 watersheds with the highest FRic values (i.e., the top 10%) were compared to the 21 watersheds with the highest FEve values to calculate how many watersheds were shared within both groups. Subsequently, spatial congruence of pairwise diversity was assessed at the successive 10% intervals. A randomization procedure was performed (999 iterations) for all six diversity indices to determine whether watershed congruence was greater than the random expectations.
We used variation partitioning based on redundancy analysis (RDA) models, with significance assessed using 999 Monte Carlo permutations, to reveal the partial effects of environmental variables and spatial structure on each diversity index (Peres-Neto et al., 2006). This procedure decomposes the total variation in the response dataset into a pure spatial component (S|E), a pure environmental component (E|S), a component of the spatial structured environmental variation (E∩S) and the unexplained variation. Only significant predictor variables were used for variation partitioning as identified from multiple linear regression models by using forward selection procedure and two stopping rules: either exceeding the critical p-value (P = 0.05) or the adjusted R2 value of the reduced model against the global model based on 999 random permutations (Blanchet et al., 2008). We ran the variation partitioning analyses by using function “varpart” in the R package “vegan” (Oksanen et al., 2016). We then computed the Moran’s I correlograms to evaluate the degree of spatial autocorrelation of the diversity indices and the residuals from the linear models (Diniz-Filho and Bini, 2005) by using function “correlog” in the R package “pgirmess” (Giraudoux, 2016).
Given the nature of our data set, sampling bias might be expected among the 214 studied watersheds. For example, watersheds that are easy to access or largely populated may be expected to be more intensively sampled. However, Pearson correlation analysis between the number of literature sources (NOL, ranging from 4 to 23; Supplementary Figures S2a,b) compiled for each watershed, used as proxy for the total sampling effort, and species richness showed no significant correlation (Supplementary Figure S2c). Furthermore, results from a sensitivity analysis conducted by repeating the analysis only on watersheds with NOL >8 (25th quantile; n = 160) remain largely invariant to those obtained from the complete data set (Supplementary Figure S2d). Therefore, we report results for all watersheds.
Higher values of taxonomic richness and most functional indices concentrated in watersheds from central-southern China (Figures 1A–F). In contrast, higher values of macrophyte taxonomic distinctness and functional evenness appeared in watersheds from western China (Figures 1B–D). Among all six diversity indices, taxonomic, and functional richness presented the highest range of variation across watersheds (Figures 1A–C), with mean values representing, respectively, 62 and 70% of the total pool of 469 documented species and functional richness across China. At the other extreme, taxonomic distinctness of individual watersheds was very homogeneous across China (Figure 1B), but accounted on average for 96% of the total distinctness in the study region.
Figure 1. Spatial patterns of taxonomic richness (A), taxonomic distinctness (B), functional richness (C), functional evenness (D), functional divergence (E), and functional dispersion (F) of freshwater macrophyte assemblages across Chinese watersheds.
Spatial congruence for the top 20% of each type of diversity indices presented two distinct groups with high (>40%) and low (<20%) pairwise congruence (Figure 2). The highest and lowest congruence occurred, respectively, for FDiv and FDis (58.6%) and FRic and FEve (0%). At this level, all high-congruence pairs were significantly more congruent than random expectations (999 iterations; p < 0.001), and showed a significant strong positive correlation (Figure 3). On the contrary, incongruent pairs showed mainly weak or strong negative (e.g., FEve – TRic and FEve – FRic) correlations (Figure 3).
Figure 2. Pairwise congruence between taxonomic richness (TRic), taxonomic distinctness (TDis), functional richness (FRic), functional evenness (FEve), functional divergence (FDiv), and functional dispersion (FDis) of freshwater macrophyte assemblages across Chinese watersheds. Congruence was assessed by comparing the spatial concordance for each pair of biodiversity facets across watersheds grouped by percentiles (10% intervals).
Figure 3. Paired scatterplots and Pearson correlations between diversity indices of freshwater macrophyte assemblages across Chinese watersheds (Abbreviations as in Figure 2). Red and blue represents positive and negative correlations, respectively.
Based on the selected linear regression model interpreting the variation of the diversity indices (Table 1), the percentage of variation explained varied from 40% for FEve to 81% for FDis and FDiv (Figure 4). Within each diversity index, the proportion of variation explained by shared fractions (spatially structured environmental gradients) was significantly higher compared to the proportion by unique fractions (pure effects), although the pure spatial component was responsible for relatively high amounts of variation in TRic, FRic, FEve, and FDiv (Figure 4). Geographic variations of the six diversity indices were strongly spatially autocorrelated along with a steady decreasing of Moran’s I coefficient across distances (Supplementary Figure S3). Most of the residuals of the models showed weak spatial patterns, with the exception of the significantly positive autocorrelation at short distances, indicating that our models captured well the major ecological factors underlying geographic gradients of the diversity facets.
Table 1. Variables retained, adjusted R2, and significance values from the forward-selected multiple regression models examining the effect of environmental and spatial factors on six aspects of macrophyte biodiversity.
Figure 4. Variation partitioning for diversity indices of freshwater macrophyte assemblages across Chinese watersheds (Abbreviations as in Figure 2).
We found a relatively high spatial congruence between many, but not all, biodiversity facets (i.e., TRic and FRic, FRic and FDis, TDis and FEve, FDiv and FDis, and FRic and FDiv). Given the strong effects of spatially structured environmental factors on shaping the biodiversity patterns (Figure 1), this high spatial congruence can presumably be explained by topography-related dispersal limitation affecting specific functional groups and, consequently, species. Our results corroborate previous evidence on strong correlations between different components of diversity. For instance, Heino et al. (2008) found that species richness was highly correlated with functional richness in stream macroinvertebrate assemblages. Likewise, Strecker et al. (2011) and Pool et al. (2014) reported high spatial congruence among three aspects of species diversity patterns of freshwater fish assemblages. Meynard et al. (2011) found evidence that hypotheses generated for local and regional taxonomic diversity were equally applicable to both phylogenetic diversity and functional diversity. These findings suggest the possibility of using a single diversity measurement as a surrogate for other facets to optimize conservation planning. Given resources are often limited, pinpointing conservation priorities and simultaneously protecting multiple diversity facets is highly desirable (Devictor et al., 2010).
On the other hand, the extremely low spatial congruence found between FEve and both TRic and FRic (Figure 2), and the low spatial congruence between TDis and both TRic and FRic may be related to the definition of the measurements, whereby increases in taxonomic diversity and functional diversity can only cause small changes in phylogenetic diversity (Schweiger et al., 2008; Pool et al., 2014; Cai et al., 2018). The scatter of the negative correlation indicates that taxonomic and functional richness increase with decreasing phylogenetic and functional diversity (Figure 3). This suggests a stronger effect of species identity relative to that of taxonomic and functional richness. For instance, watersheds from central-southern China show highest TRic and FRic but relatively low FEve and TDis, whereas watersheds from western China show relatively high FEve and TDis but low TRic and FRic. Macrophytes from watersheds in central-southern China were phylogenetically and biologically more closely related, whereas those in western China were more distantly related to each other. An increasing number of studies reveal spatial mismatches among the three aspects of species diversity patterns (Heino et al., 2005; de Carvalho and Tejerina-Garro, 2015), suggesting that species occurring locally may originate from regional species pools with distinct biogeographic and evolutionary processes, respectively (Losos, 2008; Devictor et al., 2010). In other cases, species in the same region might respond differently to environmental variables affecting spatial TDis patterns and result in mismatches among different facets of diversity (Tucker and Cadotte, 2013). Thus, it is not surprising that great longitudinal gradients in China from west to east are associated with distinct environmental filtering conditions and dispersal limitations affecting functional traits and phylogenetic diversity of macrophytes.
Our results showed that both environmental and spatial factors influence the different facets of macrophyte biodiversity. In particular, spatially structured environmental gradients, rather than pure environmental effects, shaped the different facets of macrophyte biodiversity in watersheds across China (Figure 4). Pure spatial factors had a significant role in shaping several facets of the macrophyte biodiversity patterns suggesting that dispersal limitations exert a strong effect on macrophyte assemblage structure across the different diversity facets.
We found significant spatial autocorrelations among all six diversity metrics, whereas the selected spatially structured environmental variables optimally explained the spatial structure of all diversity facets. This results highlights the role of spatially structured environmental gradients, over and above the effect of environmental factors per se, as a major driver of biodiversity, which feeds into the debate about the effects of environmental heterogeneity and dispersal limitations on species distributions (Field et al., 2009). Our results also highlight the dominant role of climatic gradients in driving spatially structured patterns of all facets of diversity across watersheds. Climatic gradients at large spatial scales can influence biodiversity patterns through multiple mechanisms related to the physiology, energetic demand and dispersal limitations of species (Fu et al., 2018; Li et al., 2018; Zhang et al., 2019). Spatial structure in the diversity facets of TRic, FRic, FEve, and FDiv was significantly explained by broad-scale dbMEMs (Supplementary Figure S1). Macrophytes seem to generally confront dispersal limitations although they are often recognized as good dispersers (Santamaria, 2002; Chambers et al., 2008). Our findings suggest that, given the presence of mountain ranges, habitat variability, and other obstacles across the studied watersheds, macrophytes assemblages across China for most diversity aspects are strongly structured by dispersal limitation. Such a pattern of species distributions is consistent with a number of previous studies for other aquatic plant assemblages at regional scales (Capers et al., 2010; Mikulyuk et al., 2011; Alahuhta and Heino, 2013; Viana et al., 2014).
Our study on variation-partitioning analysis demonstrates that macrophyte diversity patterns in watersheds across China are not always congruent and mainly driven by spatially structured environmental determinism. This finding implies that macrophyte are taxonomically, phylogenetically, and functionally clustered in space, which might be the result of environmental and/or evolutionary forces.
JX and MZ conceived the ideas and compiled the data. MZ, GS, and JX analyzed the data. MZ, JGM, GS, HZ, and JX wrote the manuscript. All authors contributed to the final manuscript.
This research was supported by the National Key R&D Program of China (2018YFD0900904), the National Natural Science Foundations of China (Grant Nos. 31872687 and 31370473), and the Water Pollution Control and Management Project of China (Grant No. 2018ZX07208005). JGM was supported by the “Tenure-Track System Promotion Program” of the Japanese Ministry of Education, Culture, Sports Science and Technology (MEXT). GS was funded by the China Scholarship Council. HZ was supported by the Hundred-Talent Program of the Chinese Academy of Sciences. There are no actual or potential competing financial interests. The funders didn’t involve in the study design, data collection, and analysis, decision to publish, or preparation of the manuscript.
Conflict of Interest Statement
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.
We are grateful to Yingying Chen and Yi Yang at the Huazhong Aquaculture University and Xiaohu Huang at Hubei University for help in data compiling and sampling.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2019.00161/full#supplementary-material
- ^ http://www.nrscc.gov.cn/nrscc/
- ^ http://db.kib.ac.cn/eflora/Default.aspx
- ^ http://www.chinaias.cn/wjPart/index.aspx
- ^ http://db.kib.ac.cn/
- ^ http://zgswz.lifescience.com.cn/
Alahuhta, J., Lindholm, M., Bove, C. P., Chappuis, E., Clayton, J., De Winton, M., et al. (2018). Global patterns in the metacommunity structuring of lake macrophytes: regional variations and driving factors. Oecologia 188, 1167–1182. doi: 10.1007/s00442-018-4294-0
Borcard, D., and Legendre, P. (2002). All-scale spatial analysis of ecological data by means of principal coordinates of neighbour matrices. Ecol. Model. 153, 51–68. doi: 10.1016/S0304-3800(01)00501-4
Cadotte, M. W., Carscadden, K. A., and Mirotchnick, N. (2011). Beyond species: functional diversity and the maintenance of ecological processes and services. J. Appl. Ecol. 48, 1079–1087. doi: 10.1111/j.1365-2664.2011.02048.x
Cai, Y., Xu, J., Zhang, M., Wang, J., and Heino, J. (in press). Different roles for geography, energy and environment in determining three facets of freshwater molluscan beta diversity across China. Sci. Total Environ. doi: 10.1016/j.scitotenv.2018.1012.1373
Cai, Y., Zhang, M., Xu, J., and Heino, J. (2018). Geographical gradients in the biodiversity of Chinese freshwater molluscs: implications for conservation. Divers. Distribut. 24, 485–496. doi: 10.1111/ddi.12695
Capers, R. S., Selsky, R., and Bugbee, G. J. (2010). The relative importance of local conditions and regional processes in structuring aquatic plant communities. Freshw. Biol. 55, 952–966. doi: 10.1111/j.1365-2427.2009.02328.x
Chappuis, E., Gacia, E., and Ballesteros, E. (2014). Environmental factors explaining the distribution and diversity of vascular aquatic macrophytes in a highly heterogeneous Mediterranean region. Aquat. Bot. 113, 72–82. doi: 10.1016/j.aquabot.2013.11.007
de Carvalho, R. A., and Tejerina-Garro, F. L. (2015). Relationships between taxonomic and functional components of diversity: implications for conservation of tropical freshwater fishes. Freshw. Biol. 60, 1854–1862. doi: 10.1111/fwb.12616
Devictor, V., Mouillot, D., Meynard, C., Jiguet, F., Thuiller, W., and Mouquet, N. (2010). Spatial mismatch and congruence between taxonomic, phylogenetic and functional diversity: the need for integrative conservation strategies in a changing world. Ecol. Lett. 13, 1030–1040. doi: 10.1111/j.1461-0248.2010.01493.x
Diniz-Filho, J. A. F., and Bini, L. M. (2005). Modelling geographical patterns in species richness using eigenvector-based spatial filters. Glob. Ecol. Biogeogr. 14, 177–185. doi: 10.1111/j.1466-822X.2005.00147.x
Dodkins, I., Rippey, B., and Hale, P. (2005). An application of canonical correspondence analysis for developing ecological quality assessment metrics for river macrophytes. Freshw. Biol. 50, 891–904. doi: 10.1111/j.1365-2427.2005.01360.x
Dray, S., Pélissier, R., Couteron, P., Fortin, M. J., Legendre, P., Peres-Neto, P. R., et al. (2012). Community ecology in the age of multivariate multiscale spatial analysis. Ecol. Monogr. 82, 257–275. doi: 10.1890/11-1183.1
Field, R., Hawkins, B. A., Cornell, H. V., Currie, D. J., Diniz-Filho, J. A. F., Guégan, J.-F., et al. (2009). Spatial species-richness gradients across scales: a meta-analysis. J. Biogeogr. 36, 132–147. doi: 10.1111/ele.12277
Flynn, D. F. B., Mirotchnick, N., Jain, M., Palmer, M. I., and Naeem, S. (2011). Functional and phylogenetic diversity as predictors of biodiversity–ecosystem-function relationships. Ecology 92, 1573–1581. doi: 10.1890/10-1245.1
Fu, H., Yuan, G., Luo, Q., Dai, T., Xu, J., Cao, T., et al. (2018). Functional traits mediated cascading effects of water depth and light availability on temporal stability of a macrophyte species. Ecol. Indicat. 89, 168–174. doi: 10.1016/j.ecolind.2018.02.010
Fu, H., Zhong, J., Yuan, G., Ni, L., Xie, P., and Cao, T. (2014). Functional traits composition predict macrophytes community productivity along a water depth gradient in a freshwater lake. Ecol. Evol. 4, 1516–1523. doi: 10.1002/ece3.1022
Gianuca, A. T., Declerck, S. A. J., Lemmens, P., and De Meester, L. (2017). Effects of dispersal and environmental heterogeneity on the replacement and nestedness components of β-diversity. Ecology 98, 525–533. doi: 10.1002/ecy.1666
Giraudoux, P. (2016). Package ‘Pgirmess’: Data Analysis in Ecology. R package Version 1.6.5. Available at: https://cran.r-project.org/web/packages/pgirmess/index.html
Hansson, L.-A., Nicolle, A., Bronmark, C., Hargeby, A., Lindstrom, A., and Andersson, G. (2010). Waterfowl, macrophytes, and the clear water state of shallow lakes. Hydrobiologia 646, 101–109. doi: 10.1007/s004420050121
Heino, J., Alahuhta, J., and Fattorini, S. (2015). Phylogenetic diversity of regional beetle faunas at high latitudes: patterns, drivers and chance along ecological gradients. Biodivers. Conserv. 24, 2751–2767. doi: 10.1007/s10531-015-0963-z
Heino, J., Mykrä, H., and Kotanen, J. (2008). Weak relationships between landscape characteristics and multiple facets of stream macroinvertebrate biodiversity in a boreal drainage basin. Landscape Ecol. 23, 417–426. doi: 10.1007/s10980-008-9199-6
Heino, J., Soininen, J., Lappalainen, J., and Virtanen, R. (2005). The relationship between species richness and taxonomic distinctness in freshwater organisms. Limnol. Oceanogr. 30, 978–986. doi: 10.4319/lo.2005.50.3.0978
Laliberté, E., Legendre, P., and Shipley, B. (2014). FD: Measuring Functional Diversity From Multiple Traits, and Other Tools for Functional ecology. Available at: http://www.elaliberte.info/publications
Li, C., Wang, T., Zhang, M., and Xu, J. (2018). Maternal Environment Effect of Warming and Eutrophication on the Emergence of Curled Pondweed, Potamogeton crispus L. Water 10:1285. doi: 10.3390/w10091285
Losos, J. B. (2008). Phylogenetic niche conservatism, phylogenetic signal and the relationship between phylogenetic relatedness and ecological similarity among species. Ecol. Lett. 11, 995–1003. doi: 10.1111/j.1461-0248.2008.01229.x
Meynard, C. N., Devictor, V., Mouillot, D., Thuiller, W., Jiguet, F., and Mouquet, N. (2011). Beyond taxonomic diversity patterns: how do α, β, and γ components of bird functional and phylogenetic diversity respond to environmental gradients across France? Glob. Ecol. Biogeogr. 20, 893–903. doi: 10.1111/j.1466-8238.2010.00647.x
Mikulyuk, A., Sharma, S., Van Egeren, S. J., and Hauxwell, J. (2011). The relative role of environmental, spatial, and land-use patterns in explaining aquatic macrophyte community composition. Can. J. Fish. Aquat. Sci. 68, 1778–1789. doi: 10.1139/f2011-095
Oksanen, J., Blanchet, F. G., Kindt, R., Legendre, P., Minchin, P. R., O’hara, R., et al. (2016). vegan: Community Ecology Package. R package version 2.4-1. Available at: https://cran.r-project.org/web/packages/vegan/index.html
Pardo, I., Roquet, C., Lavergne, S., Olesen, J. M., Gómez, D., and García, M. B. (2017). Spatial congruence between taxonomic, phylogenetic and functional hotspots: true pattern or methodological artefact? Divers. Distribut. 23, 209–220. doi: 10.1111/ddi.12511
Peres-Neto, P., Legendre, P., Dray, S., and Borcard, D. (2006). Variation partitioning of species data matrices: estimation and comparison of fractions. Ecology 87, 2614–2625. doi: 10.1890/0012-9658(2006)87[2614:VPOSDM]2.0.CO;2
Pool, T. K., Grenouillet, G., and Villéger, S. (2014). Species contribute differently to the taxonomic, functional, and phylogenetic alpha and beta diversity of freshwater fish communities. Divers. Distribut. 20, 1235–1244. doi: 10.1111/ddi.12231
Santamaria, L. (2002). Why are most aquatic plants widely distributed? Dispersal, clonal growth and small-scale heterogeneity in a stressful environment. Acta Oecol. 23, 137–154. doi: 10.1016/S1146-609X(02)01146-3
Strecker, A., Casselman, J., Fortin, M.-J., Jackson, D., Ridgway, M., Abrams, P., et al. (2011). A multi-scale comparison of trait linkages to environmental and spatial variables in fish communities across a large freshwater lake. Oecologia 166, 819–831. doi: 10.1007/s00442-011-1924-1
Tisseuil, C., Cornu, J. F., Beauchard, O., Brosse, S., Darwall, W., Holland, R., et al. (2013). Global diversity patterns and cross-taxa convergence in freshwater systems. J. Anim. Ecol. 82, 365–376. doi: 10.1111/1365-2656.12018
Tucker, C. M., and Cadotte, M. W. (2013). Unifying measures of biodiversity: understanding when richness and phylogenetic diversity should be congruent. Divers. Distribut. 19, 845–854. doi: 10.1111/ddi.12087
Viana, D. S., Santamarĺa, L., Schwenk, K., Manca, M., Hobk, A., Mjelde, M., et al. (2014). Environment and biogeography drive aquatic plant and cladoceran species richness across Europe. Freshw. Biol. 59, 2096–2106. doi: 10.1111/fwb.12410
Veech, J. A., and Crist, T. O. (2007). Habitat and climate heterogeneity maintain beta-diversity of birds among landscapes within ecoregions. Glob. Ecol. Biogeogr. 16, 650–656. doi: 10.1111/j.1466-8238.2007.00315.x
Whittaker, R. J., Willis, K. J., and Field, R. (2001). Scale and species richness: towards a general, hierarchical theory of species diversity. J. Biogeogr. 28, 453–470. doi: 10.1046/j.1365-2699.2001.00563.x
Wiegleb, G., Gebler, D., De Weyer, K. V., and Birk, S. (2016). Comparative test of ecological assessment methods of lowland streams based on long-term monitoring data of macrophytes. Sci. Total Environ. 541, 1269–1281. doi: 10.1016/j.scitotenv.2015.10.005
Wiegleb, G., Herr, W., Zander, B., Bröring, U., Brux, H., and Van De Weyer, K. (2015). Natural variation of macrophyte vegetation of lowland streams at the regional level. Limnologica 51, 53–62. doi: 10.1016/j.limno.2014.12.005
Xu, J., Su, G., Xiong, Y., Akasaka, M., García Molinos, J., Matsuzaki, S.-I., et al. (2015). Complimentary analysis of metacommunity nestedness and diversity partitioning highlights the need for a holistic conservation strategy for highland lake fish assemblages. Glob. Ecol. Conserv. 3, 288–296. doi: 10.1016/j.gecco.2014.12.004
Zhang, M., García Molinos, J., Zhang, X., and Xu, J. (2018). Functional and taxonomic differentiation of macrophyte assemblages across the Yangtze River floodplain under human impacts. Front. Plant Sci. 9:387. doi: 10.3389/fpls.2018.00387
Zhang, P., Grutter, B. M. C., Van Leeuwen, C. H. A., Xu, J., Petruzzella, A., Van Den Berg, R. F., et al. (2019). Effects of rising temperature on the growth, stoichiometry, and palatability of aquatic plants. Front. Plant Sci. 9:1947. doi: 10.3389/fpls.2018.01947
Keywords: freshwater macrophyte, functional diversity, spatial congruence, species richness, taxonomic distinctness, phylogenetic diversity
Citation: Zhang M, García Molinos J, Su G, Zhang H and Xu J (2019) Spatially Structured Environmental Variation Plays a Prominent Role on the Biodiversity of Freshwater Macrophytes Across China. Front. Plant Sci. 10:161. doi: 10.3389/fpls.2019.00161
Received: 12 July 2018; Accepted: 29 January 2019;
Published: 22 February 2019.
Edited by:Sidinei Magela Thomaz, Universidade Estadual de Maringá, Brazil
Reviewed by:Yu An, Northeast Institute of Geography and Agroecology (CAS), China
Gustavo Gonzaga Henry-Silva, Federal Rural University of the Semi-Arid Region, Brazil
Teresa Ferreira, Universidade de Lisboa, Portugal
Copyright © 2019 Zhang, García Molinos, Su, Zhang and Xu. 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: Jun Xu, firstname.lastname@example.org