Uplift-driven diversification in the Hengduan Mountains, a temperate biodiversity hotspot

Edited by Scott V. Edwards, Harvard University, Cambridge, MA, and approved February 27, 2017 (received for review September 26, 2016)
April 3, 2017
114 (17) E3444-E3451

Significance

Why do so many species occur in mountains? A popular but little-tested hypothesis is that tectonic uplift creates environmental conditions (new habitats, dispersal barriers, etc.) that increase the rate at which resident species divide and evolve to form new ones. In China’s Hengduan Mountains region, a biodiversity hotspot uplifted over the last 8 million years, this rate does in fact show a significant increase during that time, relative to the rate for adjacent older mountains, and to the rate of species immigration. The Hengduan Mountains flora is thus made up disproportionately of species that evolved within the region during its uplift, supporting the original hypothesis and helping to explain the prevalence of mountains as global biodiversity hotspots.

Abstract

A common hypothesis for the rich biodiversity found in mountains is uplift-driven diversification—that orogeny creates conditions favoring rapid in situ speciation of resident lineages. We tested this hypothesis in the context of the Qinghai–Tibetan Plateau (QTP) and adjoining mountain ranges, using the phylogenetic and geographic histories of multiple groups of plants to infer the tempo (rate) and mode (colonization versus in situ diversification) of biotic assembly through time and across regions. We focused on the Hengduan Mountains region, which in comparison with the QTP and Himalayas was uplifted more recently (since the late Miocene) and is smaller in area and richer in species. Time-calibrated phylogenetic analyses show that about 8 million y ago the rate of in situ diversification increased in the Hengduan Mountains, significantly exceeding that in the geologically older QTP and Himalayas. By contrast, in the QTP and Himalayas during the same period the rate of in situ diversification remained relatively flat, with colonization dominating lineage accumulation. The Hengduan Mountains flora was thus assembled disproportionately by recent in situ diversification, temporally congruent with independent estimates of orogeny. This study shows quantitative evidence for uplift-driven diversification in this region, and more generally, tests the hypothesis by comparing the rate and mode of biotic assembly jointly across time and space. It thus complements the more prevalent method of examining endemic radiations individually and could be used as a template to augment such studies in other biodiversity hotspots.
Central to understanding global patterns of biodiversity are considerations of biotic assembly: For a given region, when and how did resident species accumulate? Of primary interest is tempo (the rate of accumulation) and mode (the process, whether by colonization via dispersal or in situ lineage diversification). We wish to know how and why these vary in time and space.
For mountains, well-known for harboring a disproportionate fraction of terrestrial species, a common hypothesis is that of uplift-driven diversification—that orogeny creates conditions favoring in situ speciation of resident lineages (16). Among global biodiversity hotspots, the mountain ranges surrounding the Qinghai–Tibetan Plateau (QTP) are unusual and enigmatic: They harbor one of the world’s richest temperate floras, and (unlike other hotspots) they are neither tropical nor Mediterranean in climate. Moreover, despite increasing interest from biogeographers, their biotic assembly remains poorly understood (3, 4, 7). The mountains form three distinct hotspots of biodiversity that respectively lie to the west, south, and east of the QTP’s central high desert: the Central Asian mountains (Altai and Tianshan ranges), the Himalayas, and the Hengduan Mountains region (4) (Fig. 1). Of these, the richest in plant diversity is the Hengduan Mountains, with a vascular flora of about 12,000 species in an area of about 500,000 km2 (810). In this study we seek to better understand the origins of this remarkable flora through the lens of historical biogeography. In particular, was in situ diversification accelerated by uplift of the Hengduan Mountains?
Fig. 1.
Map of the Hengduan Mountains region in relation to the QTP and Himalayas.
The geological history of the QTP and its surrounding ranges is complex, and many details remain controversial, with different lines of evidence often yielding conflicting inferences. However, some points of consensus have emerged from recent syntheses (4, 1113). One is that the central plateau was uplifted first, forming a “proto-QTP” as early as 40 Mya, with subsequent outward extensions by the early Miocene (11, 14). By the late Miocene, 8 to 10 Mya, all of the mountains surrounding the QTP to the south, west, and north had reached their current elevations (12, 1517). By contrast, uplift of the Hengduan Mountains region, at the southeastern margin (Fig. 1), is generally believed to have been rapid and recent, occurring mainly between the late Miocene and late Pliocene (1823).
The Hengduan Mountains are thus younger than the rest of the QTP and relatively small in area, but conspicuously rich in species. If floristic assembly tracks orogeny, this suggests an elevated tempo since the late Miocene, compared with adjacent regions, but what of mode? Either colonization, in situ diversification, or both processes must have been accelerated. In biogeographic studies of QTP-associated clades, the hypothesis of uplift-driven in situ diversification has been clearly favored over colonization, being commonly invoked to explain phylogenetic and phylogeographic divergences (e.g., refs. 2427). However, conspicuously lacking from these studies are quantitative analyses that explicitly measure rates of diversification and colonization and compare them across regions and time (3, 4). Consequently, the idea that uplift-driven diversification has contributed disproportionately to floristic assembly in the Hengduan Mountains has yet to be rigorously tested.
In this study we used the evolutionary histories of multiple plant groups to study the floristic assembly of the Hengduan Mountains region, focusing on comparisons to adjacent regions, especially the Himalayas and other geologically older parts of the QTP (Fig. 1). We inferred regional rates of diversification and colonization through time from fossil-calibrated molecular chronograms and reconstructions of ancestral range and rates of lineage diversification, using data from 19 clades of vascular plants chosen for their potential to inform the biogeographic history of the Hengduan Mountains. Our primary aim was to discover the differences in tempo and mode of biotic assembly that help account for the remarkable diversity of the Hengduan flora.

Results

Contrasting Histories of Floristic Assembly.

Reconstructed biogeographic histories of the selected clades reveal distinct patterns of floristic assembly across regions when plotted as the cumulative number of colonization and in situ speciation events through time (Fig. 2). We focus here on contrasting the Hengduan Mountains region with a combined Himalayas–QTP region to the west and temperate/boreal East Asia to the east (Delimitation of Biogeographic Regions). Given geological and paleontological evidence that the Hengduan Mountains are younger than the Himalayas–QTP, we expected to find that its flora was assembled more recently. Instead, we found surprisingly deep phylogenetic histories of Hengduan occupancy, with initial assembly of its flora predating that of the Himalayas–QTP. In more than 95% of the pseudoreplicated joint histories the earliest colonization event occurs by 64 Ma for the Hengduan Mountains and by 59 Ma for the Himalayas–QTP; in situ speciation begins later, starting by 34 Ma and 19 Ma, respectively. By contrast, in temperate/boreal East Asia, which was commonly reconstructed as the root ancestral area for clades, in situ speciation initially occurs by 155 Ma and colonization initially occurs by 60 Ma.
Fig. 2.
Assembly of regional floras by colonization and in situ speciation events in 19 plant clades, inferred from ancestral-range reconstructions on time-calibrated molecular phylogenies. Shaded regions indicate the 5 to 95% quantile intervals for the cumulative number of events through time from 500 pseudoreplicated joint biogeographic histories designed to account for phylogenetic uncertainty (discussed in the text). Panels on the right focus on the last 20 Ma, in which differences in regional assembly are most apparent. In the Hengduan Mountains region, cumulative in situ speciation overtakes colonization about 8 Ma, whereas for the Himalayas–QTP colonization remains the dominant process. In situ speciation thus seems to have played a disproportionately large role in assembling the Hengduan Mountains flora since the late Miocene compared with the Himalayas–QTP, consistent with the theory of uplift-driven diversification in the Hengduan Mountains region.
For all regions, species accumulation by each process is approximately exponential (log-linear) following initiation until about the last 8 Ma. During this earlier “constant” phase, assembly in both the Hengduan Mountains and Himalayas–QTP is dominated by colonization, whereas in temperate/boreal East Asia the dominant process is in situ speciation. Following this phase, the assembly dynamics of the Hengduan Mountains and Himalayas diverge considerably, with relatively little change in temperate/boreal East Asia. In the Hengduan Mountains, the cumulative number of in situ speciation events overtakes that of colonization around 8 Ma, and by the present the median number is 508 versus 335. In the Himalayas–QTP in situ speciation never overtakes colonization, and by the present the median number of events is 192 versus 411. The relative contribution of in situ speciation to floristic assembly is thus 508/(508+335)=0.603 for the Hengduan Mountains, and 192/(192+411)=0.318 for the Himalayas–QTP. In other words, in situ speciation has contributed almost twice as much to the assembly of the Hengduan Mountains flora as it has to the Himalayas–QTP flora, especially since the late Miocene.

Regional Rates of in Situ Speciation and Dispersal Through Time.

The rate of in situ speciation for the Hengduan Mountains region, in events per resident lineage per million years, increased almost twofold over the past 10 Ma, whereas for the Himalayas–QTP, it remained more or less constant (Fig. 3). Dispersal rates between the Hengduan Mountains and Himalayas–QTP increase gradually over the last 4 to 5 Ma, and sharply in the last 1 Ma. By contrast, dispersal from each region to temperate/boreal East Asia increases only slightly in the same time period (Fig. 4).
Fig. 3.
Rolling estimates of in situ speciation rates through time for the Hengduan Mountains and Himalayas–QTP regions from inferred biogeographic histories of 19 plant clades. Lines indicate medians and shaded areas indicate 5 to 95% quantile intervals from 500 pseudoreplicated joint histories designed to account for phylogenetic uncertainty (discussed in the text). Regional rates begin to diverge about 8 Ma, with the Hengduan Mountains showing a striking increase in in situ speciation relative to the Himalayas–QTP.
Fig. 4.
Rolling estimates of colonization rates through time for the Hengduan Mountains, Himalayas–QTP, and temperate/boreal East Asia regions from inferred biogeographic histories of 19 plant clades. Lines indicate medians and shaded areas indicate 5 to 95% quantile intervals from 500 pseudoreplicated joint histories designed to account for phylogenetic uncertainty (discussed in the text). Dispersal between the Hengduan Mountains and Himalayas–QTP increases in the last 2 Ma relative to dispersal between either region and temperate/boreal East Asia.

Shifts in Diversification Rate.

Seven out of the 19 clades showed strong evidence (Bayes factor >15) for one or more shifts in diversification regime (Table 1). However, inspection of branch-specific evidence (the cumulative probability of occurring in the credible set of shift configurations, and the marginal odds ratio in favor of a shift) in light of ancestral-range reconstructions does not reveal any clear geographic patterns in the phylogenetic locations of regime shifts (SI Appendix, Fig. S2). That is, across clades, macroevolutionary jumps in diversification rate are not obviously associated with the colonization or occupancy of any particular region.
Table 1.
Bayes factor support for shifts in diversification in sampled clades, relative to the null hypothesis of no shifts
 No. of regime shifts
Clade1234567891011
Acer31.7671.6226.075.591.050.12     
Allium14.24182.481,048.163,836.486,668.167,038.725,401.603,015.681,382.40409.60163.84
Clematidinae + Anemoninae6.746.824.101.720.440.160.05    
Cyananthus0.960.520.180.060.020.00     
Delphineae377.435,871.716,889.714,316.571,872.00566.86137.1436.57   
Isodon8.188.605.232.230.940.190.06    
Ligularia–Cremanthodium–           
Parasenecio complex1.641.410.800.320.120.040.01    
Lilium + Nomocharis0.170.0410.01        
Meconopsis0.770.320.100.020.01      
Microsoroideae3.834.372.791.190.380.110.07    
Pinaceae40.7829.3610.152.170.340.05     
Polygoneae550.121,841.187,241.418,086.595,225.412,304.00783.06165.65   
Primulaceae2.132.181.540.770.360.090.060.02   
Rhodiola1.941.921.330.710.260.080.02    
Rhododendron17.7687.97247.00289.51215.75111.6440.479.762.30  
Rosa1.310.900.420.180.030.01     
Saussurea7.3020.2121.6813.566.262.010.35    
Saxifragaceae + Grossulariaceae1.911.520.720.260.040.00     
Thalictrum1.571.250.670.280.100.020.01    
Strong support (value > 15) is indicated in bold.
A notable exception to this general pattern is Rhododendron, in which results from BAMM (Bayesian Analysis of Macroevolutionary Mixtures) and Lagrange suggest that in a window of about 9 to 15 Ma net diversification increased independently in two clades that each originated in the Hengduan region and are currently dominated by Hengduan species (Fig. 5). In Saussurea, an increase in net diversification is inferred about 1.7 Ma along the stem of a branch containing most of the clade’s Hengduan species (SI Appendix, Fig. S2O). Similarly, in Delphineae, two separate increases are inferred in clades that are ancestrally Hengduan and together hold most of the Hengduan species; however, the times of these rate shifts (about 28 Ma and 37 Ma, respectively) predate the uplift of the Hengduan Mountains (SI Appendix, Fig. S2D). Finally, in Isodon and the Ligularia–Cremanthodium–Parasenecio complex, branch-specific measures show one increase in net diversification for each clade in the context of Hengduan ancestry about 7 and 5 Ma, respectively, but in both cases models with regime shifts are not supported by Bayes factors, and the Hengduan Mountains region is reconstructed as ancestral across most of the phylogeny, rendering the geographic context of the shift less informative.
Fig. 5.
Reconstructions of ancestral geographic range (Left) and net diversification rate (Right) on the maximum clade credibility tree, with branch lengths set to posterior means, for Rhododendron. Ancestral ranges are maximum-likelihood estimates at the start and end of each branch. Net diversification values are branch-segment means of the posterior distribution estimated by BAMM. Filled circles on the right indicate branches that appear in the 95% credible set of distinct shift configurations, with the size and label of a circle indicating the cumulative probability of the branch over all configurations in the credible set. On the left, the marginal odds ratio for a shift in diversification regime along a branch is drawn for branches where the ratio exceeds 20. Geographic regions are coded as follows: AUS, Australasia; CWA, central/western Asia; EAS, temperate-boreal East Asia; EUR, Europe; HEN, Hengduan Mountains; HIM, Himalayas–QTP; NAM, North America; and SEA, Southeast Asia. Hengduan species cluster primarily in two clades, both of which show evidence of ancestral shifts to higher diversification rate in the mid-to-late Miocene.

Discussion

Our analysis quantifies the relative contributions of in situ lineage diversification and colonization to the assembly of one of the world’s richest temperate floras, that of the Hengduan Mountains. We expected regional differences to reflect contrasting times of orogeny, in particular the prediction of the uplift-driven diversification hypothesis that in situ speciation increases with mountain-building activity. In this context, the late Miocene emerges as an important reference point, because previous studies of geology and paleontology indicate that the Hengduan Mountains achieved their current height only after this time, whereas the Himalayas and central QTP did so before (4, 13, 15, 17, 21). Here, our phylogenetic inferences, which make no prior assumptions about the timing of geological events, show that after about 8 Ma the rate of in situ diversification increased in the Hengduan Mountains (Fig. 3), yielding a remarkable inflection point at which cumulative speciation overtakes colonization (Fig. 2). This indicates that the Hengduan Mountains flora has been assembled disproportionately by recent in situ diversification that coincides temporally with independent estimates of orogeny.
Rapid Himalayan orogeny during the early to middle Miocene (15, 21, 28) would predict a similar pulse of uplift-driven diversification, but none is evident whether the Himalayas are treated as part of the QTP (Fig. 3) or separately (SI Appendix, Fig. S4). It is possible that the Himalayas were uplifted more gradually; alternatively, the signal of such a pulse could be masked by the greater uncertainty associated with estimating clade ages and ancestral ranges in deeper time, and/or subsequent extinction and turnover in the affected lineages.
In a recent review, Renner (13) argued that numerous phylogenetic studies have misinterpreted the Earth sciences literature and incorrectly attributed relatively young (Miocene and later) estimates of lineage divergence to QTP uplift, and recommended that “Biogeographical discussions should also keep in mind expected differences in the biota of the TP and those in the Himalayas, the Kunlun range to the N, the Pamir and Karakoram to the W, and Hengduan Mountains to the E” (ref. 13, p. 7). Here, our focus on distinguishing the younger and more species-rich Hengduan Mountains region from the Himalayas–QTP is in the spirit of that advice, and our results suggest that at least some of the referenced cases of diversification could be more accurately interpreted as coincident with orogeny of the Hengduan region, not the older QTP. Indeed, several of the studies cited in Renner (13) are the original sources of molecular data for the clades analyzed here and include many Hengduan species. It thus seems important to emphasize that, despite precedent in the literature for considering the Hengduan Mountains to be simply part of a greater biogeographic region that includes the Himalayas and QTP (e.g., refs. 26, 27, 29, and 30), the Hengduan region actually has a very distinct history of assembly that reflects its younger age.
Why did in situ diversification in the Hengduan Mountains increase since the late Miocene, and, more specifically, to what extent was it driven by orogeny? Further studies are needed. One might look for contrasts between the Hengduan Mountains and adjacent regions in the evolution of ecological traits and environmental tolerances (e.g., ref. 31) across multiple clades. These may reveal, for example, the extent to which speciation within a region is associated with adaptive divergence and niche-filling, as opposed to nonadaptive processes such as genetic isolation and drift arising from dispersal limitation and topographic effects (vicariance, emergence of sky islands, etc.). These analyses would require denser and more fine-grained taxonomic and geographic sampling than was possible here, in addition to the requisite trait data. However, we suspect that evidence for common processes will prove elusive, because the proximate causes of diversification are likely to be idiosyncratic and involve a variety of contingencies and associated variables (e.g., refs. 32 and 33). For many clades, factors unrelated to uplift per se, such as biotic interactions, are almost certain to have played important roles (3, 13, 34). For example, in Pedicularis (Orobanchaceae) diversification of the many Hengduan species may have been facilitated by recurrent divergence in floral traits associated with pollinator sharing and reproductive interference (e.g., ref. 35), in combination with other mechanisms promoting population isolation, including uplift and associated aspects of landscape and resource heterogeneity.
A coarser-grained view is that diversification in the Hengduan region reflects a macroevolutionary response to the rapid expansion of moist temperate montane and alpine habitats in the late Miocene. At this level the effects of orogeny and climate change are not clearly separable, because they concurrently and jointly set a stage of ecological/evolutionary opportunity. During the early to middle Miocene, the southern QTP had already achieved its current elevation (15) and was covered by coniferous and deciduous-leaved forests (36, 37). As Hengduan uplift proceeded, falling temperatures would have driven the treeline lower, increasing the extent of shrublands, meadows, and scree slopes—habitats where most of the species in our dataset occur and would have diversified. Similar changes would also have occurred in the Himalayas and QTP, but in comparison with those areas moisture in the Hengduan region would have been abundant, because by then the summer monsoon was established, and the central QTP and Central Asian interior had aridified (13). This predicts that carrying capacities—and by extension, diversification potential—of the Hengduan region’s newly formed habitats were enhanced by greater summer rainfall.
Although not a driver of diversification per se, it is important to consider the effects of glaciation and climatic oscillations during the Pleistocene, and the possibility that our results can be explained, at least in part, by regional differences in extinction caused by these factors. In the Hengduan region, the north–south orientation of major valleys suggests that extinction was relatively low, because species would have been able to persist by migrating relatively unhindered toward southern ice-free areas during glacial cycles; as a result, the region’s history of in situ diversification would be better preserved in time-calibrated phylogenies, simply because fewer branches would have been “pruned” by extinction. By contrast, in the Himalayas–QTP, the east–west orientation of valleys would have inhibited southern displacement, leading to greater extinction. In this case, the pruning effect of extinction would yield lower estimates of in situ diversification during the Pleistocene.
This pattern is probably not due solely to differences in glacial extent. Numerous studies of species distributions (e.g., refs. 38 and 39) and phylogeography (e.g., refs. 4043) support the idea that refugia occurred along the southeastern margin of the QTP as well as in areas further west and north, even at high elevation (e.g., refs. 4446). This aligns with geological evidence that ice-free areas were present across the QTP and adjacent mountains throughout the Quaternary (47). It may be that greater Pleistocene extinction in the Himalayas–QTP compared with the Hengduan region was driven as much (or more) by aridification as glaciation, because the former experienced reduced summer monsoon strength during interglacials (48).
These effects may be reflected in temporal patterns of dispersal between the Hengduan Mountains and Himalayas–QTP. The older age of the latter predicts that its flora should have been an early source of lineage dispersal to the younger Hengduan Mountains. However, we do not find any marked asymmetry in dispersal between these regions through the end of the Miocene; instead, rates in both directions are relatively low and increase only gradually until about 2 to 3 Ma, when both increase more sharply, but more so from the Hengduan Mountains to the Himalayas–QTP (Fig. 4). That is, in the Quaternary, the Hengduan Mountains seem to have acted as a biogeographic source, from which lineages buffered from extinction colonized the relatively depauperate regions to the west and north.
Other than plants, the group that has received the most biogeographic scrutiny in relation to the QTP and adjacent regions is birds, which present similar patterns (49), notably that in the species-rich Hengduan Mountains region and adjacent forests of SE Asia, in situ diversification during and subsequent to the Miocene led to repeated colonization of the Himalayas and QTP (e.g., refs. 31, 50, and 51). In the Himalayas, the lack of in situ diversification of songbirds has been attributed to the filling of niches by preadapted colonizing lineages (52).
Among studies inferring biodiversity hotspot assembly from multiple clades, quantitative analyses of in situ diversification and/or colonization are still rare, and ours estimates rates of both processes through time and across regions. This comparative framework has allowed us to characterize, at a very basic level, changes in the tempo (rate) and mode (process) of biotic assembly in relation to independently inferred timelines of regional orogeny and climate change. By contrast, studies of other hotspots have tended to focus on single regions or biomes, and subsets of these dimensions (time, space, tempo, and mode). For example, the high species diversity in the California Floristic Province is associated with lower rates of extinction more than elevated speciation or immigration, using models in which rates of lineage birth, death, and movement are dependent on geographic occupancy but not time (53). Assembly of the Cerrado biome in Brazil has been characterized by recent (late Miocene onward) in situ adaptation of lineages to fire resistance (54), with some clades, especially in the campos rupestres highlands, showing evidence of endemic radiation [e.g., Mimosa (54, 55), Calliandra (56), and Chamaecrista (57)]. In the Páramo biome of the Andes, analyses of net diversification (but not immigration) yielded the highest average rate compared with eight other hotspots and appear driven more by Pleistocene climate oscillations than orogeny (58). Clearly, for the Hengduan Mountains, much remains to be studied about how patterns of assembly through time break down at a finer scale of particular biomes and habitats within the region.
The Andes, perhaps the world’s richest mountain hotspot (2), are a compelling subject for comparison with the present results. The tropical montane/alpine habitats of the Andes are geographically isolated, whereas the temperate Hengduan region is embedded in a network of high mountains surrounding the QTP and proximate to cool high-latitude habitats. It thus stands to reason that dispersal figures prominently in Hengduan biotic assembly, compared with the Andes, in which endemic radiations are the more dominant pattern (e.g., refs. 5, 7, 34, 59, and 60).
Prevalence of dispersal is likely the main reason that we find little evidence of macroevolutionary jumps in diversification rate associated with Hengduan colonization in individual clades. In our phylogenies, clades tend not to be entirely endemic to the Hengduan Mountains or Himalayas–QTP; dispersal is common and reflects the fact that these regions are not defined by differences in biome but share broad physiographic similarities that presumably have facilitated biogeographic exchange. That is, dispersal between these regions should not necessarily require extensive ecophysiological adaptations that may be difficult to evolve (61), countering expectations of infrequent colonization events followed by endemic radiations (cf. ref. 60). For this reason our data are perhaps not well-suited for BAMM, which models diversification shifts as relatively rare events. Moreover, the link between dispersal events and shifts in diversification regime is further obscured by the inherent uncertainty associated with inferring their phylogenetic locations. The most appropriate interpretation may simply be that the signal of higher diversification in the Hengduan region is more or less diffuse across clades and emerges only when their biogeographic histories are considered jointly. This helps explain why previous analyses of single clades have not yielded the pattern inferred here.
The validity of our results hinges on the degree to which the data are representative of the focal regions (Hengduan Mountains and Himalayas–QTP). The clades we selected span a wide taxonomic range, and the species exhibit a diversity of growth forms, life histories, and habitat preferences. Taxon sampling, although incomplete, is not generally biased with respect to the focal regions (SI Appendix, Table S1) and therefore should not unduly affect inferences of the relative rates of assembly processes across regions. We also expect the phylogenetic distribution of unsampled species within each clade to be effectively random, which means that reconstructed phylogenies would tend to be missing more nodes—and thus inferred biogeographic events—closer to the present. In other words, the extent to which we have underestimated the frequencies of those events is likely greater in the Pliocene and Quaternary than in the Miocene and earlier.
Another potential source of error is in the demarcation of discrete areas as required for biogeographic analysis. Although the QTP, Himalayas, and Hengduan Mountains may have broadly distinct orogenic histories, this does not translate into clearly defined borders, and to the extent such borders existed they surely would have shifted through time, leading to ambiguity in the present. Our definitions of these areas (Fig. 1) are thus necessarily somewhat subjective and represent our best attempts to capture the essential elements of this complex system. That being said, our results are robust to the treatment of the Himalayas and QTP as a single biogeographic area; an analyis of recoded species range data treating them as separate areas yielded the same signature of uplift-driven Hengduan diversification (SI Appendix, Additional Analyses of Biotic Assembly and Figs. S3–S5).

Materials and Methods

Clade Selection and Phylogeny Reconstruction.

Our criteria were that each clade (i) included a substantial number of species that occur in the Hengduan Mountains, as well as their closest known relatives in other biogeographic regions, (ii) had sufficient molecular data available to infer a phylogeny that was broadly representative of the clade’s taxonomic diversity and geographic range, and (iii) had fossil data suitable for molecular clock calibration or secondary calibrations inferred from fossil dated phylogenies. We found 19 clades—17 of angiosperms and one each of gymnosperms (Pinaceae) and ferns (Microsorioideae)—that met these criteria. The combined taxon sample included 4,668 ingroup species, of which 930 occur in the Hengduan Mountains region, 703 in the Himalayas–QTP, and 1,045 in the remainder of temperate/boreal East Asia (Delimitation of Biogeographic Regions). About 370 species are shared between the Hengduan Mountains and Himalayas–QTP regions. Across clades, the proportion of species sampled ranged from 20 to 97% globally and 26 to 94% for the Hengduan Mountains region (SI Appendix, Table S1).
For each clade, we assembled molecular sequence alignments (Dataset S1) and fossil calibration data and used relaxed molecular clock models implemented in BEAST (Bayesian Evolutionary Analysis Sampling Trees) (62, 63) to generate a Bayesian posterior sample of time-calibrated phylogenies and the associated maximum-clade-credibility tree (SI Appendix).

Inference of Range Evolution and Lineage Diversification.

Delimitation of biogeographic regions.

We delimited 11 biogeographic regions (SI Appendix, Fig. S1) and scored the range of each sampled species as its presence or absence in these regions based on floras, online databases, and published datasets (SI Appendix and Dataset S2). Delimitations reflected our need to accommodate the global distributions of the selected clades and the granularity of available species range data, and most importantly to focus on distinguishing the Hengduan Mountains region from adjacent, geologically older parts of the QTP, especially the Himalayas (Fig. 1).
The Hengduan Mountains region, following Boufford (8), is bounded to the west by the Yarlung Tsangpo River in eastern Xizang (Tibet), to the northwest by the high plateau in Qinghai, to the north by the Tao River in southern Gansu, to the east by the Sichuan Basin, and to the south by subtropical forests and the Yunnan–Guizhou Plateau (Fig. 1). We delimited the complementary “Himalayas–QTP” region as including the plateau itself, the Himalayas to the south, the Kunlun Mountains to the north, and the Qilian Mountains to the northeast. The distinct geological histories of the Himalayas and QTP argue against their treatment as a single unit for biogeographic analysis; however, the granularity of available data on species ranges made it difficult to differentiate between them. We attempted to do so, with no effect on outcomes (SI Appendix, Additional Analyses of Biotic Assembly and Figs. S3–S5), but we are most confident with analyses based on the combined delimitation, which effectively serves the purpose of contrasting with the Hengduan region.
The other regions we delimited are temperate/boreal East Asia (the eastern boreal part of Russia and temperate regions of East Asia, including Japan and Taiwan), central/western Asia (including the Xinjiang Uyghur Autonomous Region to the north of the Kunlun Range), southeast Asia (including tropical regions of China, Malesia, and Papuasia), Australasia (Australia, New Zealand, and the southwestern Pacific), India (south of the Himalyas, including Sri Lanka), Africa, Europe, North America, and South America (south of the Panama Canal).

Ancestral range reconstruction.

We estimated ancestral geographic ranges for all ingroups using a dispersal-extinction-cladogenesis model in Lagrange (64, 65) in which dispersal was allowed based on spatial adjacency (SI Appendix, Fig. S1) and the maximum range size was set to 3. We explicitly avoided placing any temporal constraints on dispersal, so that biogeographic inferences were independent of prior beliefs about the ages of the Hengduan Mountains or Himalayas–QTP regions. For each ingroup we inferred a distribution of biogeographic histories to account for phylogenetic uncertainty. A single history was generated in three steps. First, a chronogram was randomly sampled from the clade’s Bayesian posterior distribution. Second, the maximum-likelihood set of ancestral ranges (geographic speciation scenarios) at internal nodes was estimated using Lagrange. Third, parsimonious sequences of dispersal and local extinction events were randomly interpolated along branches having different ancestral and descendant ranges. This procedure yielded a complete chronology of biogeographic events (i.e., where and when ancestral lineages moved and underwent speciation and local extinction). To account for phylogenetic uncertainty, one history was generated for each of 500 chronograms sampled from the posterior for each clade.

Regional assembly processes through time.

We focused our analysis on the dynamics of the Hengduan Mountains region in comparison with the Himalayas–QTP and temperate East Asia. Our objective was to estimate, for each region, rates of in situ diversification and colonization and their cumulative contributions to biotic assembly through time, while accounting for phylogenetic uncertainy. Taking the biogeographic histories of clades from the ancestral-range analysis, we generated 500 pseudoreplicated joint histories—sets of one history drawn randomly from each clade’s pool without replacement. From each joint history we extracted the chronology of in situ speciation, colonization, and local extinction events affecting the Hengduan Mountains, Himalayas–QTP, and temperate/boreal East Asia, and binned them into 1-My periods. This allowed region-specific estimates of assembly processes through time. We calculated rolling estimates of per capita in situ speciation rates as λ(t)=s(t)/n(t1), where s(t) is the number of in situ speciation events inferred in a region in a 1-My period t and n(t1) is the number of inferred lineages in the region in the previous period (the cumulative sum of in situ speciation and colonization minus local extinction). We also calculated rolling per capita rates of colonization between regions as dij(t)=cij(t)/ni(t1), where cij(t) is the number of inferred colonization events of area j from area i. For all estimates, confidence intervals (5 to 95% quantiles) were calculated from the pseudoreplicated joint histories. Our results thus do not condition on any particular tree topologies or branch lengths and reflect the levels of clade support, and confidence in divergence times, provided by the sequence data.

Geography-independent rates of diversification.

To complement the regional-process analysis and to better understand heterogeneity in rates of lineage diversification independently of geography, we generated Bayesian inferences of diversification rate for each clade’s maximum-clade-credibility tree, with posterior mean branch lengths, using BAMM version 2.5 and the BAMMtools R package (66). BAMM uses Markov chain Monte Carlo procedures to jointly estimate the number, parameters, and locations of distinct macroevolutionary regimes (rates of lineage birth and death, possibly time-dependent) on a phylogenetic tree. It can account statistically for incomplete taxon sampling, so we assigned sampling fractions at the finest level of taxonomic resolution possible (SI Appendix). Priors for speciation and extinction were set empirically using the setBAMMpriors function. For all clades smaller than 500 species we set the geometric prior on the expected number of regime shifts to 1, as recommended in the BAMM documentation; for larger clades we also tested a prior of 10 expected shifts. We ran each BAMM Markov chain Monte Carlo analysis for 10 million generations with a sampling frequency of 1/1,000 and assessed convergence by visually inspecting plots of the likelihood trace and calculating the effective sample size after discarding the first 10% of the run as burn-in. For each clade, we calculated Bayes factors for the distinct numbers of regime shifts sampled, marginal odds ratios (relative evidence in favor of a shift) for individual branches, and the 95% credible set of distinct shift configurations.

Acknowledgments

We thank D. Boufford and three anonymous reviewers for valuable comments on previous drafts. This work was supported by a Boyd Postdoctoral Fellowship at the Field Museum, Swiss Advanced Postdoc.Mobility Fellowship P300P3_158528, the Pioneer Hundred Talents Program of the Chinese Academy of Sciences (Y.X.), and NSF Grant DEB-1119098 (to R.H.R.).

Supporting Information

Appendix (PDF)
Dataset_S01 (XLSX)
Dataset_S02 (XLSX)

References

1
C Hoorn, V Mosbrugger, A Mulch, A Antonelli, Biodiversity from mountain building. Nat Geosci 6, 154–154 (2013).
2
CE Hughes, The tropical Andean plant diversity powerhouse. New Phytol 210, 1152–1154 (2016).
3
J Wen, JQ Zhang, ZL Nie, Y Zhong, H Sun, Evolutionary diversifications of plants on the Qinghai-Tibetan Plateau. Front Genet 5, 1–16 (2014).
4
A Favre, et al., The role of the uplift of the Qinghai-Tibetan Plateau for the evolution of Tibetan biotas. Biol Rev 90, 236–253 (2015).
5
LP Lagomarsino, FL Condamine, A Antonelli, A Mulch, CC Davis, The abiotic and biotic drivers of rapid diversification in Andean bellflowers (Campanulaceae). New Phytol 210, 1430–1442 (2016).
6
O Schwery, et al., As old as the mountains: The radiations of the Ericaceae. New Phytol 207, 355–367 (2015).
7
CE Hughes, GW Atchison, The ubiquity of alpine plant radiations: From the Andes to the Hengduan Mountains. New Phytol 207, 275–282 (2015).
8
D Boufford, Biodiversity hotspot: China’s Hengduan Mountains. Arnoldia 72, 24–35 (2014).
9
XW Li, J Li, A preliminary floristic study on the seed plants from the region of Hengduan Mountain. Yunnan-zhiwu-yanjiu 15, 217–231 (1993).
10
ZY Wu, Hengduan mountain flora and her significance. J Jpn Bot 63, 297–311 (1988).
11
C Wang, et al., Outward-growth of the Tibetan Plateau during the Cenozoic: A review. Tectonophysics 621, 1–43 (2014).
12
T Deng, L Ding, Paleoaltimetry reconstructions of the Tibetan Plateau: Progress and contradictions. Natl Sci Rev 2, 417–437 (2015).
13
SS Renner, Available data point to a 4-km-high Tibetan Plateau by 40 Ma but 100 molecular-clock papers have linked supposed recent uplift to young node ages. J Biogeogr 43, 1479–1487 (2016).
14
DB Rowley, BS Currie, Palaeo-altimetry of the late Eocene to Miocene Lunpola basin, central Tibet. Nature 439, 677–681 (2006).
15
RA Spicer, et al., Constant elevation of southern Tibet over the past 15 million years. Nature 421, 622–624 (2003).
16
X Fang, et al., Magnetostratigraphy of the late Cenozoic Laojunmiao anticline in the northern Qilian Mountains and its implications for the northern Tibetan Plateau uplift. Sci China Ser D 48, 1040–1051 (2005).
17
Y Wang, et al., Cenozoic uplift of the Tibetan Plateau: Evidence from the tectonic–sedimentary evolution of the western Qaidam Basin. Geosci Front 3, 175–187 (2012).
18
E Kirby, et al., Late Cenozoic evolution of the eastern margin of the Tibetan Plateau: Inferences from 40ar/39ar and (U-Th)/He thermochronology. Tectonics 21, 1–20 (2002).
19
MK Clark, et al., Late Cenozoic uplift of southeastern Tibet. Geology 33, 525–528 (2005).
20
E Wang, et al., Two-phase growth of high topography in eastern Tibet during the Cenozoic. Nat Geosci 5, 640–645 (2012).
21
P Wang, et al., Tectonic control of Yarlung Tsangpo Gorge revealed by a buried canyon in Southern Tibet. Science 346, 978–981 (2014).
22
K Meng, E Wang, G Wang, Uplift of the Emei Shan, western Sichuan Basin: Implication for eastward propagation of the Tibetan Plateau in Early Miocene. J Asian Earth Sci 115, 29–39 (2016).
23
BN Sun, et al., Reconstructing Neogene vegetation and climates to infer tectonic uplift in western Yunnan, China. Palaeogeogr, Palaeoclimatol, Palaeoecol 304, 328–336 (2011).
24
JQ Liu, YJ Wang, AL Wang, O Hideaki, RJ Abbott, Radiation and diversification within the Ligularia–Cremanthodium–Parasenecio complex (Asteraceae) triggered by uplift of the Qinghai-Tibetan Plateau. Mol Phylogenet Evol 38, 31–49 (2006).
25
Y Wang, A Susanna, E von Raab-Straube, R Milne, J Liu, Island-like radiation of Saussurea (Asteraceae: Cardueae) triggered by uplifts of the Qinghai–Tibetan Plateau. Biol J Linn Soc 97, 893–903 (2009).
26
JQ Zhang, SY Meng, GA Allen, J Wen, GY Rao, Rapid radiation and dispersal out of the Qinghai-Tibetan Plateau of an alpine plant lineage Rhodiola (Crassulaceae). Mol Phylogenet Evol 77, 147–158 (2014).
27
YD Gao, AJ Harris, SD Zhou, XJ He, Evolutionary events in Lilium (including Nomocharis, Liliaceae) are temporally correlated with orogenies of the Q–T plateau and the Hengduan Mountains. Mol Phylogenet Evol 68, 443–460 (2013).
28
MP Searle, Geological evolution of the Karakoram Ranges. Ital J Geosci 130, 147–159 (2011).
29
ZL Nie, et al., Molecular phylogeny of Anaphalis (Asteraceae, Gnaphalieae) with biogeographic implications in the Northern Hemisphere. J Plant Res 126, 17–32 (2013).
30
S Matuszak, AN Muellner-Riehl, H Sun, A Favre, Dispersal routes between biodiversity hotspots in Asia: The case of the mountain genus Tripterospermum (Gentianinae, Gentianaceae) and its close relatives. J Biogeogr 43, 580–590 (2016).
31
Y Liu, et al., Sino-Himalayan mountains act as cradles of diversity and immigration centres in the diversification of parrotbills (Paradoxornithidae). J Biogeogr 43, 1488–1501 (2016).
32
MJ Donoghue, MJ Sanderson, Confluence, synnovation, and depauperons in plant diversification. New Phytol 207, 260–274 (2015).
33
Y Bouchenak-Khelladi, RE Onstein, Y Xing, O Schwery, HP Linder, On the complexity of triggering evolutionary radiations. New Phytol 207, 313–326 (2015).
34
F Luebert, M Weigend, Phylogenetic insights into Andean plant diversification. Evol Popul Genet 2, 27 (2014).
35
DAR Eaton, CB Fenster, J Hereford, SQ Huang, RH Ree, Floral diversity and community structure in Pedicularis (Orobanchaceae). Ecology 93, S182–S194 (2012).
36
J Sun, et al., Palynological evidence for the latest Oligocene-early Miocene paleoelevation estimate in the Lunpola Basin, central Tibet. Palaeogeogr, Palaeoclimatol, Palaeoecol 399, 21–30 (2014).
37
H Li, S Guo, The Miocene flora Namling of Xizang. Gu sheng wu xue bao 15, 7–17 (1976).
38
U Srinivasan, K Tamma, U Ramakrishnan, Past climate and species ecology drive nested species richness patterns along an east-west axis in the Himalaya: Nestedness in Himalayan fauna. Global Ecol Biogeogr 23, 52–60 (2014).
39
J López-Pujol, FM Zhang, HQ Sun, TS Ying, S Ge, Centres of plant endemism in China: Places for survival or for speciation? J Biogeogr 38, 1267–1280 (2011).
40
YZ Cun, XQ Wang, Plant recolonization in the Himalaya from the southeastern Qinghai-Tibetan Plateau: Geographical isolation contributed to high population differentiation. Mol Phylogenet Evol 56, 972–982 (2010).
41
B Wang, JF Mao, J Gao, W Zhao, XR Wang, Colonization of the Tibetan Plateau by the homoploid hybrid pine Pinus densata. Mol Evol 20, 3796–3811 (2011).
42
F Lei, Y Qu, G Song, Species diversification and phylogeographical patterns of birds in response to the uplift of the Qinghai-Tibet Plateau and Quaternary glaciations. Curr Zool 60, 149–161 (2014).
43
L Meng, et al., Refugial isolation and range expansions drive the genetic structure of Oxyria sinensis (Polygonaceae) in the Himalaya-Hengduan Mountains. Sci Rep 5, 10396 (2015).
44
L Wang, et al., History and evolution of alpine plants endemic to the Qinghai-Tibetan Plateau: Aconitum gymnandrum (Ranunculaceae). Mol Evol 18, 709–721 (2009).
45
YS Sun, H Ikeda, YJ Wang, JQ Liu, Phylogeography of Potentilla fruticosa (Rosaceae) in the Qinghai-Tibetan Plateau revisited: A reappraisal and new insights. Plant Ecol Divers 3, 249–257 (2010).
46
L Opgenoorth, et al., Tree endurance on the Tibetan Plateau marks the world’s highest known tree line of the Last Glacial Maximum. New Phytol 185, 332–342 (2010).
47
LA Owen, JM Dortch, Nature and timing of Quaternary glaciation in the Himalayan–Tibetan orogen. Quaternary Sci Rev 88, 14–54 (2014).
48
LA Owen, MW Caffee, RC Finkel, YB Seong, Quaternary glaciation of the Himalayan–Tibetan orogen. J Quaternary Sci 23, 513–531 (2008).
49
M Päckert, J Martens, YH Sun, DT Tietze, Evolutionary history of passerine birds (Aves: Passeriformes) from the Qinghai–Tibetan plateau: From a pre-Quarternary perspective to an integrative biodiversity assessment. J Ornithol 156, 355–365 (2015).
50
DT Tietze, M Päckert, J Martens, H Lehmann, YH Sun, Complete phylogeny and historical biogeography of true rosefinches (Aves: Carpodacus). Zool J Linn Soc 169, 215–234 (2013).
51
US Johansson, et al., Build-up of the Himalayan avifauna through immigration: A biogeographical analysis of the Phylloscopus and Seicercus warblers. Evolution 61, 324–333 (2007).
52
TD Price, et al., Niche filling slows the diversification of Himalayan songbirds. Nature 509, 222–225 (2014).
53
LT Lancaster, KM Kay, Origin and diversification of the California flora: Re-examining classic hypotheses with molecular phylogenies. Evolution 67, 1041–1054 (2013).
54
MF Simon, et al., Recent assembly of the Cerrado, a neotropical plant diversity hotspot, by in situ evolution of adaptations to fire. Proc Natl Acad Sci USA 106, 20359–20364 (2009).
55
EJM Koenen, et al., Exploring the tempo of species diversification in legumes. S Afr J Bot 89, 19–30 (2013).
56
ER de Souza, et al., Phylogeny of Calliandra (Leguminosae: Mimosoideae) based on nuclear and plastid molecular markers. Taxon 62, 1200–1219 (2013).
57
JG Rando, et al., Phylogeny of Chamaecrista ser. Coriaceae (Leguminosae) unveils a lineage recently diversified in Brazilian campo rupestre vegetation. Int J Plant Sci 177, 3–17 (2015).
58
S Madriñán, AJ Cortés, JE Richardson, Páramo is the world’s fastest evolving and coolest biodiversity hotspot. Front Genet 4, 192 (2013).
59
RT Pennington, et al., Contrasting plant diversification histories within the Andean biodiversity hotspot. Proc Natl Acad Sci USA 107, 13783–13787 (2010).
60
C Hughes, R Eastwood, Island radiation on a continental scale: Exceptional rates of plant diversification after uplift of the Andes. Proc Natl Acad Sci USA 103, 10334–10339 (2006).
61
MJ Donoghue, EJ Edwards, Biome shifts and niche evolution in plants. Annu Rev Ecol Evol Systemat 45, 547–572 (2014).
62
AJ Drummond, MA Suchard, D Xie, A Rambaut, Bayesian phylogenetics with BEAUti and the BEAST 1.7. Molecular Biology And Evolution 29, 1969–1973 (2012).
63
R Bouckaert, et al., BEAST 2: A software platform for Bayesian evolutionary analysis. PLoS Comput Biol 10, e1003537 (2014).
64
RH Ree, BR Moore, CO Webb, MJ Donoghue, A likelihood framework for inferring the evolution of geographic range on phylogenetic trees. Evolution 59, 2299–2311 (2005).
65
RH Ree, SA Smith, Maximum likelihood inference of geographic range evolution by dispersal, local extinction, and cladogenesis. Syst Biol 57, 4–14 (2008).
66
DL Rabosky, Automatic detection of key innovations, rate shifts, and diversity-dependence on phylogenetic trees. PloS One 9, e89543 (2014).

Information & Authors

Information

Published in

Go to Proceedings of the National Academy of Sciences
Go to Proceedings of the National Academy of Sciences
Proceedings of the National Academy of Sciences
Vol. 114 | No. 17
April 25, 2017
PubMed: 28373546

Classifications

Submission history

Published online: April 3, 2017
Published in issue: April 25, 2017

Keywords

  1. biogeography
  2. vascular plants
  3. molecular clocks
  4. dispersal
  5. speciation

Acknowledgments

We thank D. Boufford and three anonymous reviewers for valuable comments on previous drafts. This work was supported by a Boyd Postdoctoral Fellowship at the Field Museum, Swiss Advanced Postdoc.Mobility Fellowship P300P3_158528, the Pioneer Hundred Talents Program of the Chinese Academy of Sciences (Y.X.), and NSF Grant DEB-1119098 (to R.H.R.).

Notes

This article is a PNAS Direct Submission.
See Commentary on page 4275.

Authors

Affiliations

Yaowu Xing1
Life Sciences Section, Integrative Research Center, The Field Museum, Chicago, IL 60605;
Key Laboratory of Tropical Forest Ecology, Xishuangbanna Tropical Botanical Garden, Chinese Academy of Sciences, Yunnan 666303, China;
Richard H. Ree2,1 [email protected]
Life Sciences Section, Integrative Research Center, The Field Museum, Chicago, IL 60605;
National Institute of Ecology, 1210 Geumgang-ro, Maseo-myeon, Seocheon-gun, Chungcheongnam-do, South Korea

Notes

2
To whom correspondence should be addressed. Email: [email protected].
Author contributions: Y.X. and R.H.R. designed research, performed research, analyzed data, and wrote the paper.
1
Y.X. and R.H.R. contributed equally to this work.

Competing Interests

The authors declare no conflict of interest.

Metrics & Citations

Metrics

Note: The article usage is presented with a three- to four-day delay and will update daily once available. Due to ths delay, usage data will not appear immediately following publication. Citation information is sourced from Crossref Cited-by service.


Citation statements




Altmetrics

Citations

If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click Download.

Cited by

    Loading...

    View Options

    View options

    PDF format

    Download this article as a PDF file

    DOWNLOAD PDF

    Get Access

    Login options

    Check if you have access through your login credentials or your institution to get full access on this article.

    Personal login Institutional Login

    Recommend to a librarian

    Recommend PNAS to a Librarian

    Purchase options

    Purchase this article to get full access to it.

    Single Article Purchase

    Uplift-driven diversification in the Hengduan Mountains, a temperate biodiversity hotspot
    Proceedings of the National Academy of Sciences
    • Vol. 114
    • No. 17
    • pp. 4267-E3585

    Media

    Figures

    Tables

    Other

    Share

    Share

    Share article link

    Share on social media