Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Phylogenetic Structure of Tree Species across Different Life Stages from Seedlings to Canopy Trees in a Subtropical Evergreen Broad-Leaved Forest

  • Yi Jin,

    Affiliation College of Life Sciences, Zhejiang University, Hangzhou, Zhejiang, China

  • Hong Qian,

    Affiliation Research and Collections Center, Illinois State Museum, Springfield, IL, United States of America

  • Mingjian Yu

    fishmj202@hotmail.com

    Affiliation College of Life Sciences, Zhejiang University, Hangzhou, Zhejiang, China

Abstract

Investigating patterns of phylogenetic structure across different life stages of tree species in forests is crucial to understanding forest community assembly, and investigating forest gap influence on the phylogenetic structure of forest regeneration is necessary for understanding forest community assembly. Here, we examine the phylogenetic structure of tree species across life stages from seedlings to canopy trees, as well as forest gap influence on the phylogenetic structure of forest regeneration in a forest of the subtropical region in China. We investigate changes in phylogenetic relatedness (measured as NRI) of tree species from seedlings, saplings, treelets to canopy trees; we compare the phylogenetic turnover (measured as βNRI) between canopy trees and seedlings in forest understory with that between canopy trees and seedlings in forest gaps. We found that phylogenetic relatedness generally increases from seedlings through saplings and treelets up to canopy trees, and that phylogenetic relatedness does not differ between seedlings in forest understory and those in forest gaps, but phylogenetic turnover between canopy trees and seedlings in forest understory is lower than that between canopy trees and seedlings in forest gaps. We conclude that tree species tend to be more closely related from seedling to canopy layers, and that forest gaps alter the seedling phylogenetic turnover of the studied forest. It is likely that the increasing trend of phylogenetic clustering as tree stem size increases observed in this subtropical forest is primarily driven by abiotic filtering processes, which select a set of closely related evergreen broad-leaved tree species whose regeneration has adapted to the closed canopy environments of the subtropical forest developed under the regional monsoon climate.

Introduction

Understanding how the phylogenetic structure of a forest changes as individual plant development proceeds within the forest and what role forest gaps play in driving phylogenetic community structure of forest regeneration is crucial to understanding the assembly processes that shape forest structure and species composition [13]. Within a forest, the formation of local assemblage might be affected by multiple assembly processes, including abiotic filtering, biotic interactions and neutral processes [1, 46]. Among the species of a regional species pool that could disperse into a local assemblage, abiotic filtering tends to select ecologically similar species [78]; biotic interactions might exclude both ecologically similar species [910] and ecologically dissimilar species [1113]; and neutral processes disregard ecological differences among species [5, 14]. Depending on the degree of phylogenetic conservatism of critical ecological traits, species ecological similarities can be approximated by species phylogenetic relatedness [1516]. Specifically, a clustered phylogenetic structure could be observed when abiotic filtering is a dominant force in community assembly, causing species to be more closely related than expected by chance [1, 15], although competitive exclusion driven by fitness difference may also lead to clustered phylogenetic structure [11]; an overdispersed phylogenetic structure (i.e., co-occurring species are more distantly related than expected by chance) could be caused by niche partitioning driven by biotic interactions [11, 15], facilitation or abiotic filtering on convergently evolved clades [1]; a random phylogenetic structure can result from neutral processes [5], lack of phylogenetic conservatism, or combined effects of abiotic filtering and biotic interactions [17].

Different size classes of tree stems within a primary (or old-growth) forest stand may be considered as a putative proxy for forest temporal dynamics [3, 1819]; thus investigating changes in phylogenetic community structure across tree size classes might shed light on the community assembly processes during individual plant development [1820]. Specifically, an increasingly overdispersed phylogenetic structure as tree stem size increases might be due to differential mortality caused by competition (e.g., [2022]); an increasingly clustered phylogenetic structure might be due to differential mortality driven by abiotic filtering (e.g., [1516, 23]), or could be due to differential mortality caused by fitness difference [11]. However, nearly all relevant previous studies were conducted in tropical regions, where abiotic filtering may not play an important role in governing community assembly [24]. Furthermore, few studies have included seedlings (but see [19, 25]). Because tree seedlings represent the initial developmental stage of a forest, including tree seedlings in a study investigating changes in phylogenetic structure through individual plant development would increase our understanding of community assembly and dynamics.

The subtropical zone, which is a transitional zone between tropical and temperate zones, is narrow in all regions of the world except in China where it covers a large area (~2.5 million km2) due to the monsoon climate primarily caused by the uplift of the Qinghai-Tibetan Plateau [26]. This region is characterized by highly complex topography [26], and topographic complexity has been detected to be positively related to species niche differentiation [27]. We expect abiotic filtering in this region to be stronger than the tropics in general. Biotic interactions in the tropics are expected to be stronger than the extratropics [28]. Therefore, the relative strength of abiotic filtering might be stronger and biotic interactions might be weaker in subtropical regions than the tropics. As a result, the increasingly overdispersed phylogenetic structure as tree stem size increases, which has been frequently observed in forests of the tropics (e.g., [18, 20]), might switch to increasingly clustered phylogenetic structure as tree stem size increases in forests from the subtropical region (expectation 1). A previous study on subtropical forests has indeed shown that phylogenetic relatedness of tree species increases from saplings to canopy trees [29], although that study did not include seedlings.

Previous studies have shown that environment is a major determinant of phylogenetic structure of communities [20, 3032]. Forest gaps are commonly formed in primary (or old-growth) forests due to canopy tree fall. Forest gap environments are characterized by higher resource (e.g., light) levels and lack of canopy protection from harsh weather conditions as compared with forest understory [33], which resembles to some degree the early successional forest habitats [34]. Theoretically, early successional habitats are expected to favor pioneer species that are adapted to a high resource environment, while species interactions are expected to be relatively weak [3538]. As a result, one would expect to find higher phylogenetic relatedness among species within communities in forest gaps than in forest understory (expectation 2). Further, previous studies have shown that phylogenetic structure of communities is scale-dependent [18]. Environmental conditions in forest gaps are positively related to forest gap size [39], thus large forest gaps would resemble to a larger degree the early successional forest habitats than small forest gaps [34], one would expect to find increasing phylogenetic relatedness among species within communities in forest gap as the size of forest gap increases (expectation 3). In addition, it is expected that canopy trees in primary (or old-growth) forests mainly are primary species that are adapted to regenerate in forest understory [40]. Therefore, if forest understory favors primary species regeneration while forest gap favors pioneer species regeneration, seedlings in forest gap might be less similar to canopy trees, compared with seedlings in forest understory. We would expect to find higher phylogenetic turnover between canopy trees and forest gap seedlings than between canopy trees and forest understory seedlings (expectation 4).

In this study, corresponding to the four expectations noted above, we intend to test the following four predictions for a forest in subtropical China: (1) Phylogenetic relatedness of tree species increases (rather than decreases as observed in [18, 20] in tropical forests) with increasing size classes of tree stems in the subtropical forest of this study. (2) Phylogenetic relatedness of seedlings in a forest gap is higher than that under the canopy of the forest surrounding the forest gap. (3) The size of a forest gap is positively related to the phylogenetic relatedness of seedlings in the forest gap. (4) Phylogenetic turnover (i.e., phylogenetic beta diversity) between seedlings under the canopy of a forest and trees in the canopy of the forest is smaller than phylogenetic turnover between seedlings in a forest gap within the forest and trees in the canopy of the forest.

Materials and Methods

The study area and sampling plots

Our study site was located at the Gutianshan National Nature Reserve (29°10'19"–29°17'41" N, 118°03'50"–118°11'12" E) in Zhejiang Province, China. The reserve is ~81 km2 in area, and was set up to protect the typical old-growth evergreen broad-leaved forest (EBLF) in this humid subtropical region. Local annual mean temperature is 15.3°C, with January minimum temperature being −6.8°C and July maximum temperature being 38.1°C; annual mean precipitation recorded from 1958 to 1986 was 1964 mm [41]. In 2005, a 24-ha (400 m × 600 m) forest dynamics plot (FDP; 29°15'6"–29°15'21" N, 118°7'1"–118°7'24" E; [42]) of evergreen broad-leaved forest was set up in the core area of the reserve following the protocol of the Center for Tropical Forest Science [43]. The forest is dominated by such evergreen broad-leaved tree species as Castanopsis eyrei (maximum diameter at breast height of 87.4 cm, relative abundance 8.81%) and Schima superba (83.0 cm, 6.04%) and is a typical EBLF. All individuals of tree species with diameter at breast height (DBH) ≥ 1 cm were identified, measured, and georeferenced during the summer of 2005. Detailed descriptions for the reserve and FDP can be found in Zhu et al. [42]. The Administration Bureau of Gutianshan National Nature Reserve issued the permission of carrying out this study.

The FDP experienced a severe ice-storm damage in 2008 [44], leading to numerous canopy openings (i.e., forest gaps). During the summer of 2010, we surveyed the structure and distribution of forest gaps caused by the ice storm within the FDP. We determined forest gaps following Brokaw’s [45] definition of forest gap (i.e., a vertical hole extending from the canopy down to an average height of about two meters above ground). The area of a forest gap was estimated using Runkle’s [46] approximation method. Because the vast majority (84.3%) of the FDP was located in two topographic habitat types (i.e., low valleys and low ridges), we confined our sampling within these two topographic habitats. We divided the identified gaps into four size classes: 25–100 m2, 100–200 m2, 200–300 m2, and 300–500 m2 [4748]. Within each size class of each topographic habitat type, we randomly selected four gaps, which resulted in a total of 32 forest gaps.

For each forest gap, we set up a sampling plot which included all area within the radius of 30 m from the center of the forest gap. This 30-m distance threshold was chosen to avoid density-dependent effects within a distance of 20 m [49] on one hand and not to be too large so that including different micro-habitat types within a sampling plot was avoided or minimized on the other hand. We obtained the data of tree species composition and abundance for all individuals with DBH ≥ 1 cm from the database resulting from the 2005 survey (see above). Within each sampling plot, we set up numerous sampling quadrats of 1 × 1 m for investigating seedlings; 905 sampling quadrats were set up for the 32 sampling plots. Within each sampling plot, sampling quadrats were separated by 2.5−5 m in distance and were evenly distributed across both forest understory and forest gap. The number of sampling quadrats varies among sampling plots (28.28 ± 6.91 SD). We tagged and identified every alive seedling (< 1 cm DBH) of tree species within the quadrats in August of 2011, and reinvestigated the seedlings in August of 2013. In this study, we included those seedlings that were alive at the times of both 2011 and 2013 surveys to assure that all seedlings used in this study were established at least for two years. We divided all recorded individuals of trees in a sampling plot into four stem size classes: seedlings (DBH < 1 cm), saplings (1 cm ≤ DBH < 5 cm), treelets (5 cm ≤ DBH < 10 cm); trees (DBH ≥ 10 cm; for ease of discussion, hereafter, we called them “canopy” trees although some of them may not enter forest canopies). For each sampling plot, sampling quadrats under the forest canopy of the plot were combined to generate a list of species with the number of seedlings per species documented; similarly, sampling quadrats in the forest gap were also combined to generate a list of species with the number of seedlings per species documented; seedlings under the forest canopy were compared with those in the forest gap.

One hundred and twenty one (85.2%) of the 142 angiosperm tree species in the entirety of the 24-ha Gutianshan FDP were found in the 32 sampling plots, although the sampling plots covered only less than 37.7% of the total area of the FDP. The relative abundance of tree species in major clades (e.g., Fagales, Laurales, and Ericales) in the 32 sampling plots was consistent with that in the entire FDP not only for the combination of all tree stems of sapling, treelets, and canopy trees but also for each of the three groups of tree stems (τ = 1, P < 0.05 in all cases, Kendall’s τ test for correlations; see S1 Fig).

Species pool and phylogenetic tree

The 24-ha Gutianshan FDP comprises 144 tree species, 142 of which are angiosperms. These species were described as either trees or small trees in Flora of China (http://www.efloras.org/) and Zheng [50] and are larger than 10 cm DBH at their maturity. We followed many previous studies (e.g., [3, 51]) to exclude non-angiosperm species because non-angiosperm species are few in each sampling plot (only two gymnosperm species in the Gutianshan FDP) but their much longer branches in a phylogeny, compared to those of angiosperms, would contribute to unusually high phylogenetic measures [3] and thus may obscure phylogenetic patterns of angiosperms. The 142 angiosperm tree species were included in the species pool of all phylogenetic analyses in this study. These species belong to 89 genera in 40 families and 21 orders [52].

The phylogeny used in this study was pruned from the phylogeny constructed for the Gutianshan FDP based on three common barcoding genes (rbcLa, matK and trnH-psbA; [53]). This phylogeny was well resolved at the species level and included 159 species. Eleven tree species found in our study were absent from the phylogeny; accordingly, we added them to the phylogeny by placing them as polytomies of their most closely related species in the phylogeny (S2 Fig).

Data analyses

Net relatedness index (NRI) is a commonly used index to quantify phylogenetic relatedness (e.g., [3, 51, 54]). Accordingly, we used NRI in our study. It is defined as [1]: where MPD is the mean pairwise phylogenetic distance between all species within an assemblage, MPDsample is the observed MPD, MPDrandom is the expected MPD of randomized assemblages, and sdMPDrandom is the standard deviation of the MPD for the randomized assemblages. Randomized (null) assemblages were generated by randomly drawing species from the species pool to replace the species in the sample. A positive NRI value indicates that MPD is lower than expected by chance (i.e., species more closely related than expected by chance) and that phylogenetic clustering of species occurs; a negative NRI value results when the observed MPD is greater than expected by chance (i.e., species more distantly related than expected by chance) and thus indicates phylogenetic evenness or overdispersion. For each NRI value, 999 randomized assemblages were generated. Species were weighted by their abundance (i.e., the number of individuals) within each assemblage. For each sampling plot, we calculated NRI separately for seedlings in the forest gap, seedlings under the forest canopy, saplings, treelets, and canopy trees. In addition to calculating NRI for each assemblage of seedlings with DBH < 1 cm, we also calculated NRI for each assemblage of seedlings with height < 50 cm, which is a stricter definition for seedlings (see [19]). NRI was calculated using the “comstruct” function in the Phylocom 4.2 software [55] and null model 2 in Phylocom 4.2 (i.e., randomly drawing species from the phylogeny pool).

To determine the degree to which phylogenetic turnover between seedlings in the two habitats (i.e., forest understory vs. forest gap) and between those in either habitat and canopy trees, we utilized the βNRI index, which is related to NRI, in analyzing phylogenetic turnover between assemblages [55] with species abundance accounted for. The index represents the negative standardized effect size (−1 × SES) of inter-community mean phylogenetic distance (βMPD). The formula of this index is [55]: where βMPD is the mean pairwise phylogenetic distance of all species between two assemblages, βMPDsample is the observed βMPD, βMPDrandom is the expected βMPD between two randomized assemblages, and sdβMPDrandom is the standard deviation of the βMPD for the two randomized assemblages. For each βNRI value, 999 pairs of randomized assemblages were generated to estimate βMPDrandom and sdβMPDrandom. Randomized assemblages were generated by randomly sampling species from the species pool to replace the species in the samples (null method 2 in Phylocom 4.2, http://phylodiversity.net/phylocom; [55]). A higher value of βNRI indicates a higher phylogenetic similarity and thus a lower phylogenetic turnover between assemblages compared. βNRI was calculated using the “comdist” function in Phylocom 4.2 [55].

We used Wilocoxon signed-rank test to determine statistic significance for pairwise comparisons [18]. We used Cochran-Armitage approach to test for directional change in the phylogenetic structure with increasing tree stem size, performed by the “independence_test” function of the “coin” package in R software version 3.1.0 [5657]. We examined whether differences in gap size influenced patterns of phylogenetic structure of seedlings in forest gaps by correlating forest gap size with values of phylogenetic indices for seedlings in forest gap using Kendall’s τ as correlation coefficient [58].

Results

NRI values for seedlings in forest understory and in forest gaps were not significantly different from those expected by chance (Fig 1). NRI values for treelets and canopy trees were significantly higher than expected by chance (Fig 1). NRI values for canopy trees were significantly greater than those for seedlings and saplings (Fig 1). There was a clear trend for NRI to increase with increasing stem size of trees (P < 0.001, Cochran-Armitage test for trend for either seedlings in forest gaps, saplings, treelets and trees or seedlings in forest understory, saplings, treelets and trees; Fig 1), indicating an increasing trend in phylogenetic clustering from seedling stratum to forest canopy stratum, regardless of whether seedlings are defined as individuals with DBH < 1 cm or as individuals of < 50 cm in height (compare Fig 1 with S3 Fig). Forest gap size showed no significant influence on patterns of phylogenetic structure of seedlings (P > 0.1, Kendall’s τ coefficient for correlations between forest gap sizes and the two phylogenetic measures of seedlings in forest gaps; see S1 Table).

thumbnail
Fig 1. NRI of different size classes of tree stems.

Box represents the interquartile range and the median (the thick line); the whiskers extend to the most extreme data point within 1.5 times interquartile range from the box. Values above and below the dotted line indicate phylogenetic clustering and overdispersion, respectively. G stands for seedlings in forest gaps, NG for seedlings in forest understory, S1 for saplings (1 cm ≤ DBH < 5 cm); S2 for treelets (5 cm ≤ DBH < 10 cm), S3 for canopy trees (DBH ≥ 10 cm). Different lower case letters above the boxes indicate significant differences (P < 0.05) in multiple comparisons (P values were adjusted using the “Bonferroni” method). Asterisks above the boxes indicate that the dispersions significantly differed from zero (P < 0.05).

https://doi.org/10.1371/journal.pone.0131162.g001

The mean value of βNRI between seedlings and canopy trees was significantly greater than zero regardless of whether seedlings under forest canopies or in forest gaps were considered (Fig 2). Values of βNRI between seedlings in forest gap and canopy trees were significantly lower than those between seedlings in forest understory and canopy trees (P < 0.05, Wilocoxon signed-rank test; Fig 2). In other words, the analysis using βNRI indicates that phylogenetic turnover was higher between seedlings in forest gap and canopy trees than between seedlings in forest understory and canopy trees.

thumbnail
Fig 2. Comparison of βNRI values for paired assemblages of seedlings in forest gaps versus canopy trees (G vs. S3) with those of seedlings in forest understory versus canopy trees (NG vs. S3).

Different lower case letters above the boxes indicate significant differences (P < 0.05).

https://doi.org/10.1371/journal.pone.0131162.g002

Discussion

Since an early attempt of performing a phylogenetic analysis of community structure [15], the use of phylogenetic approaches in investigating patterns and processes of community assembly has blossomed [59]. Investigating patterns of phylogenetic structure among different stem size classes of tree species within forests can gain insights into mechanisms of forest community assembly. Several studies (e.g., [1820]) have investigated phylogenetic patterns of size classes of tree species, but nearly all of the studies are restricted to wet tropical regions where climatic conditions are generally benign and abiotic filtering is generally a minor factor driving community assembly of forests, compared to forests outside of tropical regions [24, 60]. Furthermore, seedlings of tree species are important to the formation and development of a forest, but to our knowledge only two of the studies that investigated patterns of phylogenetic structure across a series of different size classes of tree stems have included seedlings [19, 25]. To our knowledge, our study is the first to investigate patterns of phylogenetic structure across different size classes from seedlings to canopy trees in old-growth forests of subtropical regions. In addition, we have investigated the phylogenetic structure of seedlings both in forest understory and in forest gaps.

In the subtropical forest of our study, we found increasingly clustered phylogenetic structure as stem size increases. Significant phylogenetic clustering patterns were observed even in the treelet size class, a size class prior to the canopy tree size class (Fig 1). Our finding is consistent with Webb et al.’s [19] finding for a tropical forest, where they found increasingly clustered phylogenetic structures with increasing plant developmental stages from seedlings to adult trees. Our finding is also in agreement with that of Feng et al. [29] for angiosperm trees in undisturbed forests in the same subtropical region as our study, although their study did not include seedlings and differences in NRI among the three stem size classes in their study are statistically not significant.

Although 44% of the 142 angiosperm tree species included in our study site (i.e., the 24-ha Gutianshan FDP) are deciduous species, nearly all canopy trees are evergreen broad-leaved species. The relative abundance of trees with DBH ≥ 10 cm in our study site was 78.1% for evergreen broad-leaved trees, 15.4% for deciduous broad-leaved trees, and 6.5% for coniferous trees. These characterize a typical feature of the primary (or zonal) EBLFs in this region. Many of the canopy tree species are congeners, confamilials, or taxa in closely related families. For example, of the 80 evergreen broad-leaved tree species found in our study site, ten species are in Lauraceae and nine species are in Fagaceae, eight of which belong to Castanopsis, Cyclobalanopsis and Lithocarpus which are typical dominant canopy trees in Chinese subtropical forests, and Castanopsis eyrei is the most dominant canopy tree species in our study site. In addition, nine species in Theaceae found in our study site are evergreen. Evergreen broad-leaved tree species of these families and other typical subtropical-distributed families such as Aquifoliaceae (with ten evergreen species in our study site) are the dominant species forming the canopies of forests in this subtropical region [26]. Our study has shown clear shifts in the relative abundance of tree species in some major clades. For example, species in the order Fagales (especially in the family Fagaceae) are the most common taxa of forest canopies in our study site; the relative abundance of the species in Fagales clearly increases along the series from saplings to canopy trees (S1 Fig). Conversely, the relative abundance of the species in the order Ericales decreases along the same series (S1 Fig).

Previous studies have shown that phylogenetic structures of early successional forest communities are often more clustered than primary or old-growth forest communities in tropical regions [3, 21, 54, 61] as well as in the subtropical region [29] where our study was conducted. In the old-growth forest, because forest gaps are disturbed habitats and microenvironmental conditions (e.g., light, air temperature, soil moisture) in forest gaps may resemble to some degree the initial stage of forest succession [34], one would expect to find a higher degree of phylogenetic relatedness for seedlings in forest gaps than seedlings in forest understory [22, 62]. However, we did not find support for this prediction in the forest of our study. Further, because resource levels (e.g., light) in forest gaps are positively related to gap size [39], large gaps might resemble the early successional forest environments to a larger degree than small gaps and thus gap size might be positively related to seedling phylogenetic relatedness. However, we did not find support for this prediction either. Du et al. [63] compared the seed arrival patterns in forest gaps with forest understory in the studied forest, and found that both forest gaps and forest understory were subjected to dispersal limitation. We suspect that the lack of influence of forest gap as compared with forest understory and the lack of influence of forest gap size on seedling phylogenetic relatedness in this forest might be partly due to the prevalent limitation in seed dispersal in both forest gaps and forest understory of this forest [63]. Furthermore, some of the seedlings in forest gaps might be the advance regeneration that recruited before the formation of forest gaps, which would to some extent mask the effects of forest gap environments on seedling recruitment. Nonetheless, we found that phylogenetic turnover between canopy trees and seedlings in forest understory was slightly lower than that between canopy trees and seedlings in forest gaps (Fig 2), which suggests that forest gap favored the establishment of species that are phylogenetically more dissimilar with canopy trees than forest understory.

In sum, our study shows that phylogenetic relatedness generally increases from seedlings through saplings and treelets up to canopy trees in the humid subtropical forest that we investigated. Our study also shows that phylogenetic relatedness does not differ between seedlings in forest understory and those in forest gaps, but phylogenetic turnover between canopy trees and seedlings in forest understory is lower than that between canopy trees and seedlings in forest gaps. Evergreen broad-leaved forests of subtropical regions are an important type of forests not only in China [64] but also in the world [65]. Because our study is, to our knowledge, the first to investigate phylogenetic structure for humid subtropical evergreen broad-leaved forests that includes all vegetation strata from seedlings to canopy trees, we believe our study would help to gain a fuller understanding of mechanisms of community assembly of forests in extratropical regions.

Supporting Information

S1 Fig. Relative abundance of five major orders in three stem size classes for the entire 24-ha Gutianshan FDP (A) and the 32 sampling plots (B).

https://doi.org/10.1371/journal.pone.0131162.s001

(DOC)

S2 Fig. The phylogenetic tree of the 142 angiosperm tree species found in the 24-ha Gutianshan FDP.

The following species in the phylogeny are additional to Liu et al.’s (2013) phylogeny: Castanea henryi, Cerasus discoidea, Euonymus centidens, Fraxinus chinensis, Houpoea officinalis, Ilex formosana, Ilex litseifolia, Maclura tricuspidata, Mallotus apelta, Meliosma rigida, and Trema cannabina var. dielsiana.

https://doi.org/10.1371/journal.pone.0131162.s002

(DOC)

S3 Fig. NRI of different size classes of tree stems.

G stands for seedlings of < 50 cm tall in forest gaps, NG for seedlings of < 50 cm tall in forest understory. See Fig 1 for detailed interpretations.

https://doi.org/10.1371/journal.pone.0131162.s003

(DOC)

S1 Table. Correlations between forest gap sizes and phylogenetic measures of seedlings in forest gaps.

https://doi.org/10.1371/journal.pone.0131162.s004

(DOC)

S2 Table. Data of NRI and βNRI.

ψ, seedling is defined as individual with DBH < 1 cm; φ, seedling is defined as individual with height < 50 cm. See Fig 1 and Fig 2 for key to abbreviations.

https://doi.org/10.1371/journal.pone.0131162.s005

(XLS)

S3 Table. Abundances of the five major clades in the entire FDP and the 32 sampling plots.

See Fig 1 for key to abbreviations.

https://doi.org/10.1371/journal.pone.0131162.s006

(XLS)

Acknowledgments

We are grateful to anonymous reviewers for their constructive comments and Lauren E. Barry for her assistance in designing the research protocol and conducting the field work of seedling survey. We thank Bingyang Ding, Teng Fang and Shengwen Chen for their assistance in species identification, and Tufa He, Zaobing Wang, Tubao He, Youbao Liao and other local technicians, as well as Ji Zhou, Yafang Ma and other students for their contributions to the field work. We also thank Jinlong Zhang and Xiaojuan Liu for offering the barcoding phylogeny data. We acknowledge Keping Ma, Xiangcheng Mi, Jianhua Chen, Haibao Ren, Fangliang He, I-Fang Sun and other contributors for establishment and census of the 24-ha Gutianshan forest dynamics plot, and the Administration Bureau of Gutianshan National Nature Reserve for support on this study.

Author Contributions

Conceived and designed the experiments: MY YJ HQ. Performed the experiments: YJ MY. Analyzed the data: YJ HQ. Contributed reagents/materials/analysis tools: YJ HQ MY. Wrote the paper: YJ HQ MY.

References

  1. 1. Webb CO, Ackerly DD, McPeek MA, Donoghue MJ. Phylogenies and community ecology. Annual review of ecology and systematics 2002;33: 475–505.
  2. 2. Dinnage R. Disturbance alters the phylogenetic composition and structure of plant communities in an old field system. PloS ONE 2009;49: e7071. pmid:19763265
  3. 3. Letcher SG. Phylogenetic structure of angiosperm communities during tropical forest succession. Proceedings of the Royal Society B: Biological Sciences 2009;277: 97–104. pmid:19801375
  4. 4. Weiher E, Keddy P, editors. Ecological assembly rules: perspectives, advances, retreats. Cambridge: Cambridge University Press; 2001.
  5. 5. Hubbell SP. The unified neutral theory of biodiversity and biogeography. Princeton: Princeton University Press; 2001
  6. 6. Cavender-Bares J, Kozak KH, Fine PV, Kembel SW. The merging of community ecology and phylogenetic biology. Ecology Letters 2009;12: 693–715. pmid:19473217
  7. 7. Keddy PA. Assembly and response rules: two goals for predictive community ecology. Journal of Vegetation Science 1992;32: 157–164.
  8. 8. Weiher E, Clarke GP, Keddy PA. Community assembly rules, morphological dispersion, and the coexistence of plant species. Oikos 1998;81: 309–322.
  9. 9. Hardin G. The competitive exclusion principle. Science 1960;131: 1292–1297. pmid:14399717
  10. 10. Tilman D. Competition and biodiversity in spatially structured habitats. Ecology 1994;75: 2–16.
  11. 11. Mayfield MM, Levine JM. Opposing effects of competitive exclusion on the phylogenetic structure of communities. Ecology Letters 2010;13: 1085–1093. pmid:20576030
  12. 12. HilleRisLambers J, Adler PB, Harpole WS, Levine JM, Mayfield MM. Rethinking community assembly through the lens of coexistence theory. Annual Review of Ecology, Evolution, and Systematics 2012;43: 227–248.
  13. 13. Goberna M, Navarro-Cano JA, Valiente-Banuet A, García C, Verdú M. Abiotic stress tolerance and competition-related traits underlie phylogenetic clustering in soil bacterial communities. Ecology Letters 2014;17: 1191–1201. pmid:25130277
  14. 14. Bell G. Neutral macroecology. Science 2001;293: 2413–2418. pmid:11577228
  15. 15. Webb CO. Exploring the phylogenetic structure of ecological communities: an example for rain forest trees. The American Naturalist 2000;156: 145–155. pmid:10856198
  16. 16. Baraloto C, Hardy OJ, Paine C, Dexter KG, Cruaud C, Dunning LT, et al. Using functional traits and phylogenetic trees to examine the assembly of tropical tree communities. Journal of Ecology 2012;100: 690–701.
  17. 17. Soliveres S, Torices R, Maestre FT. Environmental conditions and biotic interactions acting together promote phylogenetic randomness in semi-arid plant communities: new methods help to avoid misleading conclusions. Journal of Vegetation Science 2012;23: 822–836. pmid:25983536
  18. 18. Swenson NG, Enquist BJ, Thompson J, Zimmerman JK. The influence of spatial and size scale on phylogenetic relatedness in tropical forest communities. Ecology 2007;88: 1770–1780. pmid:17645023
  19. 19. Webb CO, Gilbert GS, Donoghue MJ. Phylodiversity-dependent seedling mortality, size structure, and disease in a Bornean rain forest. Ecology 2006;87: S123–S131. pmid:16922308
  20. 20. Yang J, Zhang G, Ci X, Swenson NG, Cao M, Sha L, et al. Functional and phylogenetic assembly in a Chinese tropical tree community across size classes, spatial scales and habitats. Functional Ecology 2014;28: 520–529.
  21. 21. Ding Y, Zang R, Letcher SG, Liu S, He F. Disturbance regime changes the trait distribution, phylogenetic structure and community assembly of tropical rain forests. Oikos 2012;121: 1263–1270.
  22. 22. Norden N, Letcher SG, Boukili V, Swenson NG, Chazdon R. Demographic drivers of successional changes in phylogenetic structure across life-history stages in plant communities. Ecology 2012;93: S70–S82. pmid:22202194
  23. 23. Parmentier I, Réjou-Méchain M, Chave J, Vleminckx J, Thomas DW, Kenfack D, et al. Prevalence of phylogenetic clustering at multiple scales in an African rain forest tree community. Journal of Ecology 2014;102: 1008–1016.
  24. 24. Myers JA, Chase JM, Jiménez I, Jørgensen PM, Araujo-Murakami A, Paniagua-Zambrana N, et al. Beta-diversity in temperate and tropical forests reflects dissimilar mechanisms of community assembly. Ecology Letters 2013;16: 151–157. pmid:23113954
  25. 25. Zhu Y, Comita LS, Hubbell SP, Ma K. Conspecific and phylogenetic density-dependent survival differs across life stages in a tropical forest. Journal of Ecology. 2015 Apr 20. https://doi.org/10.1111/1365-2745.12414
  26. 26. Wu ZY. Vegetation of China. Beijing: Science Press; 1980.
  27. 27. Brown C, Burslem D, Illian JB, Bao L, Brockelman W, Cao M, et al. Multispecies coexistence of trees in tropical forests: spatial signals of topographic niche differentiation increase with environmental heterogeneity. Proceedings of the Royal Society B: Biological Sciences 2013;280: 20130502. pmid:23782876
  28. 28. Schemske DW, Mittelbach GG, Cornell HV, Sobel JM, Roy K. Is there a latitudinal gradient in the importance of biotic interactions? Annual review of ecology and systematics 2009;40: 245–269.
  29. 29. Feng G, Svenning JC, Mi X, Jia Q, Rao M, Ren H, et al. Anthropogenic disturbance shapes phylogenetic and functional tree community structure in a subtropical forest. Forest Ecology and Management 2014;313: 188–198.
  30. 30. Algar AC, Kerr JT, Currie DJ. Quantifying the importance of regional and local filters for community trait structure in tropical and temperate zones. Ecology 2011;92: 903–914. pmid:21661553
  31. 31. Qian H, Zhang Y, Zhang J, Wang X. Latitudinal gradients in phylogenetic relatedness of angiosperm trees in North America. Global Ecology and Biogeography 2013;22: 1183–1191.
  32. 32. Kembel SW, Hubbell SP. The phylogenetic structure of a neotropical forest tree community. Ecology 2006;87: S86–S99. pmid:16922305
  33. 33. Denslow JS. Tropical rainforest gaps and tree species diversity. Annual review of ecology and systematics 1987;18: 431–451.
  34. 34. Prescott CE. The influence of the forest canopy on nutrient cycling. Tree physiology 2002;22: 1193–1200. pmid:12414379
  35. 35. Chazdon RL. Chance and determinism in tropical forest succession. In: Carson WP, Schnitzer SA, editors. Tropical forest community ecology. Oxford: Blackwell Publishing Ltd; 2008. pp. 384–408.
  36. 36. Brokaw NV. Gap-phase regeneration in a tropical forest. Ecology 1985;66: 682–687.
  37. 37. Schnitzer SA, Carson WP. Treefall gaps and the maintenance of species diversity in a tropical forest. Ecology 2001;82: 913–919.
  38. 38. Fukami T. Community assembly dynamics in space. In: Verhoef HA, Morin PJ, editors. Community Ecology: Processes, Models, and Applications. New York: Oxford University Press; 2010. pp. 45–54.
  39. 39. Denslow JS, Ellison AM, Sanford RE. Treefall gap size effects on above- and below-ground processes in a tropical wet forest. Journal of Ecology 1998;86: 597–609.
  40. 40. Yamamoto SI. Forest gap dynamics and tree regeneration. Journal of Forest Research 2000;5: 223–229.
  41. 41. Yu M, Hu Z, Yu J, Ding B, Fang T. Forest vegetation types in Gutianshan natural reserve in Zhejiang. Journal of Zhejiang University (Agriculture and Life Sciences) 2000; 27: 375–380.
  42. 42. Zhu Y, Zhao G, Zhang L, Shen G, Mi X, Ren H, et al. Community composition and structure of Gutianshan forest dynamic plot in a mid-subtropical evergreen broad-leaved forest, East China. Chinese Journal of Plant Ecology 2008;32: 262–273.
  43. 43. Condit R. Tropical forest census plots. Berlin: Springer-Verlag; 1998.
  44. 44. Man X, Mi X, Ma K. Effects of an ice strom on community structure of an evergreen broad-leaved forest in Gutianshan National Natural Reserve, Zhejiang Province. Biodiversity Science 2011;19: 197–205.
  45. 45. Brokaw NV. The definition of treefall gap and its effect on measures of forest dynamics. Biotropica 1982;14: 158–160.
  46. 46. Runkle JR. Patterns of disturbance in some old-growth mesic forests of eastern North America. Ecology 1982;63: 1533–1546.
  47. 47. de Römer AH, Kneeshaw DD, Bergeron Y. Small gap dynamics in the southern boreal forest of eastern Canada: Do canopy gaps influence stand development?. Journal of Vegetation Science 2007;18: 815–826.
  48. 48. Hubbell SP, Foster RB, O'Brien ST, Harms KE, Condit R, Wechsler B, et al. Light-gap disturbances, recruitment limitation, and tree diversity in a neotropical forest. Science 1999;283: 554–557. pmid:9915706
  49. 49. Hubbell SP, Ahumada JA, Condit R, Foster RB. Local neighborhood effects on long-term survival of individual trees in a Neotropical forest. Ecological Research 2001;16: 859–875.
  50. 50. Zheng C. Keys of seed plants in Zhejiang Province. Hangzhou: Zhejiang Science and Technology Press; 2005.
  51. 51. Qian H, Swenson NG, Zhang J. Phylogenetic beta diversity of angiosperms in North America. Global Ecology and Biogeography 2013;22: 1152–1161.
  52. 52. Angiosperm Phylogeny Group. An update of the Angiosperm Phylogeny Group classification for the orders and families of flowering plants: APG III. Botanical Journal of the Linnean Society 2009;161: 105–121.
  53. 53. Liu X, Swenson NG, Zhang J, Ma K. The environment and space, not phylogeny, determine trait dispersion in a subtropical forest. Functional Ecology 2013;27: 264–272.
  54. 54. Whitfeld T, Kress W, Erickson D, Weiblen G. Change in community phylogenetic structure during tropical forest succession: evidence from New Guinea. Ecography 2012;35: 821–830.
  55. 55. Webb CO, Ackerly DD, Kembel SW. Phylocom: software for the analysis of phylogenetic community structure and trait evolution. Bioinformatics 2008;24: 2098–2100. pmid:18678590
  56. 56. Hothorn T, Hornik K, Van De Wiel MA, Zeileis A. Package ‘coin’: conditional inference procedures in a permutation test framework. Available: http://cran.R-project.org/. 2014; version 1.0–24.
  57. 57. R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. Available: http://www.R-project.org/. 2014
  58. 58. Iida Y, Kohyama TS, Swenson NG, Su SH, Chen CT, Chiang JM, et al. Linking functional traits and demographic rates in a subtropical tree community: the importance of size dependency. Journal of Ecology 2014;102: 641–650.
  59. 59. Qian H, Jiang L. Phylogenetic community ecology: integrating community ecology and evolutionary biology. Journal of Plant Ecology 2014;7: 97–100.
  60. 60. Dobzhansky T. Evolution in the tropics. American Scientist 1950;38: 209–211.
  61. 61. Letcher SG, Chazdon RL, Andrade A, Bongers F, van Breugel M, Finegan B, et al. Phylogenetic community structure during succession: evidence from three Neotropical forest sites. Perspectives in Plant Ecology, Evolution and Systematics 2012;14: 79–87.
  62. 62. Helmus MR, Keller WB, Paterson MJ, Yan ND, Cannon CH, Rusak JA. Communities contain closely related species during ecosystem disturbance. Ecology Letters 2010;13: 162–174. pmid:20015255
  63. 63. Du Y, Mi X, Ma K. Comparison of seed rain and seed limitation between community understory and gaps in a subtropical evergreen forest. Acta Oecologica 2012;44: 11–19. pmid:22623591
  64. 64. Song YC. On the global position of the evergreen broad-leaved forests of China. In: Box EO, Peet RK, Masuzawa T, Yamada I, Fujiwara K, Maycock PF, editors. Vegetation Science in Forestry. Dordrecht: Kluwer Academic Publishers; 1995. pp. 69–84.
  65. 65. Archibold OW. Ecology of world vegetation. London: Chapman & Hall Ltd; 1995.