Abstract

Distribution patterns of alien species in nonnative ranges might be driven by a combination of various mechanisms, including phylogenetic history, competition for resource, environmental filtering, and so on. Both phylogenetic and functional limitations might work synergistically to determine the distribution of alien species. In this report, by utilizing the information of provincial distribution, functional attributes, and phylogeny for 95 alien species of China, the corresponding phylogenetic and functional community structures are evaluated. The results show that introduction pathway, life form, and flowering time of alien plants of China processed significant phylogenetic clustering patterns, while both the origin of distribution and reproduction mode of alien species showed phylogenetic overdispersion patterns, as revealed by NRI/NTI indices. The phylogenetic signal tests using Pagel’s statistic and Blomberg et al.’s statistic further verified the previous patterns, even though there are some inconsistencies. Through partial Mantel test, it is found that compositional patterns of alien plant community were mainly affected by phylogenetic limitation but not functional limitation. Conclusively, phylogeny plays a more important role in structuring provincial distribution of alien plants in China.

1. Introduction

Alien species may possess novel functional traits for successful colonization of new habitats, in addition to random events or similar traits to native species [13]. However, some recent studies have pointed out that phylogeny is also playing some role on influencing alien diversity and distribution patterns [47]. Functional attributes and phylogenetic history reflected different facets of biological diversity, and the growing trend in ecological studies nowadays is to combine both to address the questions related to the diversity [810], distribution [11], and ecosystem functioning [12, 13].

The application of combining both functional and phylogenetic filtering simultaneously to understand the distributional pattern of alien species in nonnative ranges has been concerned in recent studies [14]. It would be of great value to explore their interplay using alien species as the studied model. Alien species have expanded their distributional ranges in recent decades due to globalization and modernization [15], the processes of which required the rapid ecological adaptation [16, 17]. The rapid range expansion indicated the functional uniqueness or commonness of traits when compared alien and native plants [1, 3].

At another perspective, in addition to novel functional traits, species’ evolutionary history might be critical to promote the invasiveness, as evidenced by some previous studies [7, 1719]. Phylogenetic signals have been widely found in the interaction between native and alien species [20, 21]. Phylogenetical similarity for alien species with respect to native species were predicted to be more successful to invade novel areas [4, 7]. However, there are still some controversial debates on this issue, based on the predictions of Darwin’s naturalization and enemy release hypotheses [22, 23].

The success of invasiveness of species in non-native areas might be due to a combination of various physiological (e.g., life form, flowering time, and reproduction mode), morphological (invasive mode), and physical (origin of distribution and colonizing habitats) factors, which could shape the distributional patterns and phylogenetic community structure of species [7]. Of course, these factors could also constitute the functional diversity patterns of species [24]. At the meantime, the taxonomic classification and trait information of alien plants in China have been well documented [15, 2527]. Environmental correlates of alien plant at provincial levels of China have been revealed too [28]. However, up to now, there are no comparative studies working on the phylogenetic community structure analysis of alien plants of China as far as I know. In light with the previous information and to fill the just-mentioned knowledge gap, the central objective of the present study is to explore functional and phylogenetic community structures of the provincial distributions of alien plants in China.

The following principal hypotheses would be tested accordingly. First, alien species’ provincial distribution in China was predicted to have similar phylogenetic and functional clustering (or overdispersion) patterns because the convergence of functional traits might be related to adaptive evolution and phylogenetic conservatism [12, 29]. Further, I predicted that the imprint of phylogenetic history is more important (i.e., clustering patterns should be more remarkable in phylogenetic structure of species in comparison to functional structure). The assumption underpinning this prediction is that alien species’ traits and distribution were synergistically determined by the interaction of evolutionary history, environmental filtering, and ecological adaptation because of rapid range expansion [5, 7].

Second, I predicted that alien plants showed phylogenetic clustering patterns within each of the functional groups for a given trait, because closely related species in their traits was expected to be closely related in their phylogenetic history [7, 18] for facilitating invasive success [7, 30]. I considered five functional attributes for testing the existence of phylogenetic structure: origin of distribution, introduction pathway, flowering time, reproduction mode, and life form. I utilized and compared a suite of statistical methods for this test, including NRI/NTI indices [31], Pagel’s stati1stic [32], and Blomberg et al.'s statistic [33].

2. Materials and Methods

2.1. Dataset

A dataset, including both the trait information and distribution records at 28 provinces, is collected for 95 alien plant species found in China (there should be more than 100 alien plants presented in China, but only the previous 95 species have detailed information as far as I can get access) [15, 25, 26, 28, 34]. The following functional attributes were considered: life form, reproduction mode, flowering time, origin of distribution and introduction pathway. Rigorously speaking, origin of distribution, and introduction pathway are not functional traits, but because they are related to suitable habitat and environmental tolerance of species [35, 36], both are still treated as functional attributes (not say “traits”).

When performing analyses, origin of distribution is set to have five categories: America, Europe, Asia, Africa, and Mediterranean. Flowering time is set to have five categories: spring, summer, autumn, winter, and year-around. Life form is set to have five categories: annual grass, perennial grass, liane, shrub, and arbor. Introduction pathway is set to have four categories: human, animal, human-animal, and abiotic factors. Finally, reproduction mode is set to have two categories: seed and seed-vegetative.

2.2. Construction of the Phylogenetic Tree

The phylogenetic tree for the alien species was constructed by inputting the checklist of alien plants found in China into the plant phylogeny database PHYLOMATIC [37]. The Nexus files obtained from PHYLOMATIC was imported into the community phylogenetic software PHYLOCOM (version 4.1, available at http://phylodiversity.net/phylocom/) [31]. The Angiosperm Phylogeny Group III classification [38] was taken as the background supertree [39, 40]. Branch lengths of the resultant phylogeny was calibrated using the PHYLOCOM command “bladj” with known molecular and fossil dates of plants [39, 41, 42]. The reconstructed phylogenetic relationship for 95 alien plants was presented in Figure 1.

2.3. Construction of the Functional Dendrogram

I followed the methods introduced by previous studies [13, 24, 25, 43] to calculate the functional dendrogram using clustering analysis with alternative algorithms (including Ward, average, and complete methods). The distance matrix was constructed using the Gower distance [44]. The constructed dendrograms [44] were analyzed in subsequent analyses.

2.4. Phylogenetic and Functional Structures of Alien Plants in China

I calculated two indices, the standardized net relatedness index (NRI) and the nearest taxon index (NTI) for the distribution of 95 alien plants in China across different provinces to measure their phylogenetic and functional relatedness. These two indices [18, 31] were not only initially applied to phylogenetic community structure analysis but also applied to functional community analysis in my study: where represents the mean pairwise phylogenetic distance in which it finds the average distance to all other taxa in the sample for each taxon and is to calculate the nearest phylogenetic neighbor in the sample for each taxon [31]. The and represent the mean and mean from randomly generated assemblages. Negative values of both metrics indicate overdispersion, while positive values of both metrics indicate clustering [45]. These two standardized indices should follow the unit normal distribution with mean = 0 and variance = 1. As such, if any of the indices is larger than 1.96 (or less than −1.96), the clustering (or overdispersion) pattern identified is regarded to be statistically significant.

2.5. Phylogenetic and Functional Beta Diversity of the Distributional Patterns of Alien Plants in China

I calculated a suite of distance indices to measure phylogenetic and functional beta diversity of alien species in different provinces of China. The following presence-absence-based metrics were calculated [10, 42, 46]. Please note that all these metrics were applied to both functional dendrogram and phylogenetic cladogram comparatively. (1) Phylogenetic Sorensen index: where is the total length of the branches shared between communities and , and and are the total branch lengths found in communities and , respectively. (2)UniFrac index: where is the number of branches in the phylogeny, is the length of branch , and and are the numbers of species descended from branch in communities and . and are the total numbers of species in communities and , respectively. (3) (mean nearest phylogenetic neighbor distance between two communities): where is the distance from the nearest phylogenetic neighbor in community to species in community . is the number of species in community . The definition of and are analogous.(4) (mean pairwise phylogenetic distance between two communities): where is the mean pairwise phylogenetic distance between species in community to all the species in community . The definition of is analogous.

It is worthy to mention here that the last two indices actually were very similar to and indices, but being applied to pairwise communities (i.e., and indices). In contrast, and indices are alpha indices and applied to a single community.

After generating these phylogenetic/functional distance matrices, I employed partial Mantel test [47, 48] to examine the independent influence of phylogenetic and functional constraints on structuring species distributional patterns of alien plants. The compositional dissimilarity (or distance) matrix constructed from species-provincial distribution matrix was used as the response matrix.

3. Results

3.1. Phylogenetic and Functional Structures of the Provincial Distribution of Alien Plants

As showed in Table 1, phylogenetic clustering patterns were prevailing for the distribution across different provinces (except species found in Henan Province, indicated by negative values in both metrics). In contrast, weak functional clustering patterns were observed, and two provinces showed overdispersion patterns (Shaanxi and Shanxi). Half part of other provinces has clustering patterns, while the others showed random (no clustering or overdispersion) patterns since one of the NRI/NTI indices is negative, while the other is positive. As such, the first hypothesis was supported, showing that phylogenetic history determined provincial distribution of alien plants in China but not functional attributes.

3.2. Alien Plant Phylogenetic Signals for Functional Attributes

As revealed by the NRI/NTI indices (Table 2), the following functional attributes showed remarkable phylogenetic clustering patterns, including introduction pathway (except human-mediated species, indicated by negative values in both metrics), the life form (except liane species, indicated by negative values in both metrics), and flowering time (except summer flowering species, indicated by negative values in both metrics) for the alien species. In contrast, it was found that both the origin of distribution (except America-originated species, indicated by positive values in both metrics) and reproduction mode of alien species had phylogenetic overdispersion patterns in principle (Table 2).

The phylogenetic signal tests using Pagel’s statistic and Blomberg et al.'s statistic basically verified the findings for independent functional attributes (Table 3), but there existed some inconsistencies as well. For Pagel’s statistic, life form was found to have marginally significant phylogenetic signal, but NRI/NTI indices and Blomberg et al.'s statistic supported the existence of significant phylogenetic signal. In contrast, both Pagel’s and Blomberg et al.'s statistics identified that the origin of distribution processed significant phylogenetic clustering patterns, while NRI/NTI indices supported phylogenetic overdispersion patterns. Except the previous cases, all methods consistently identified phylogenetic clustering patterns in other functional attributes (Tables 2 and 3).

3.3. Relative Importance of Phylogenetic and Functional Constraints Structuring Alien Distribution

The partial Mantel test only had a significant result for the correlation between species’ Sorensen’s distributional dissimilarity (calculated from species-site matrix) and phylogenetic distance matrices from index when functional dendrogram was created using Ward method (Table 4). Such a significant correlation between phylogenetic distance (from index) and distributional dissimilarity matrices was found again in another situation when functional clustering of species was recalculated using the average method. Across all the partial Mantel tests when different clustering algorithms were introduced for calculating functional distance matrix, there were no significant correlations between distributional distance and functional distance matrices (Table 4).

4. Discussion

Alien plants, including their distribution and ecological/economic implications in China, have been broadly studied previously [15, 25, 26, 34]. Compared to these pervious studies, the present study made a first attempt to reveal the association of distributional patterns and functional attributes of the species by considering both phylogenetic and functional relatedness of species. Three methods for testing phylogenetic clusering signals consistently identified that the life form, flowering time, introduction pathway, and provincial distribution had strong phylogenetic clustering signals, while reproduction mode did not present any phylogenetic clustering signals. Origin of range of species did not present strong phylogenetic clustering structure, because NRI/NTI indices failed to identify possible phylogenetic clustering structure for it. As such, not all the functional attributes presented phylogenetic signals, contradictory to the second prediction.

4.1. Phylogenetic versus Functional Constraints on Alien Plant Community Structure

Through partial Mantel tests, it was ascertained that phylogenetic constraint played a more important role shaping alien species distribution patterns in China. The conclusion was still valid even when I considered performing different clustering strategies (including Ward, average, and complete methods) to obtain different functional dendrograms and recalculate functional diversity and distances matrix. The role of functional constraint could not be observed throughout, but the significant signal of phylogenetic constraint emerged in two situations (Table 4). This result further supported the first hypothesis, indicating that phylogeny of species has played a vital role on structuring alien plants’ distribution.

4.2. On the Comparison of Different Methods for Measuring Phylogenetic Signals

NRI/NTI indices were more feasible for evaluating phylogenetic signal patterns for categorical or discrete data since their calculations could allow checking the clustering patterns of a single category of a trait. However, the indices are not applicable to continuous traits.

Pagel’s and Blomberg et al.'s statistics are more general and could be applied to both discrete and continuous data. They provided an overall evaluation of phylogenetic signals of species distribution across communities. However, these tests did not separate each category of a trait and investigate individual clustering/overdispersion patterns. The potential issue of doing so as NRI/NTI indices was that the dependence of communities could not be taken into account.

There were some incongruent aspects among three methods (Tables 2 and 3). Therefore, it was recommended to utilize different methods for comparison so as to ensure the accurate detection of phylogenetic signals in the trait data.

Acknowledgments

The author is grateful to Dr. Xiangcheng Mi (Institute of Botany, Chinese Academy of Sciences) for his helpful comments and polishing on an earlier version of the paper. Associate editors and anonymous reviewers are appreciated too for their constructive comments. The work was supported by the University of British Columbia.