Abstract
Understanding which factors have driven the evolutionary success of a group is a fundamental question in biology. Angiosperms are the most successful group in plants and have radiated and adapted to various habitats. Among angiosperms, legumes are a good example for such successful radiation and adaptation. We here investigated how the interplay of past climate changes, geographical expansion and habit shifts has promoted diversification of the phaseoloid legumes, one of the largest clades in the Leguminosae. Using a comprehensive genus-level phylogeny from three plastid markers, we estimate divergence times, infer habit shifts, test the phylogenetic and temporal diversification heterogeneity, and reconstruct ancestral biogeographical ranges. We found that the phaseoloid lineages underwent twice dramatic accumulation. During the Late Oligocene, at least six woody clades rapidly diverged, perhaps in response to the Late Oligocene warming and aridity, and a result of rapidly exploiting new ecological opportunities in Asia, Africa and Australia. The most speciose lineage is herbaceous and began to rapidly diversify since the Early Miocene, which was likely ascribed to arid climates, along with the expansion of seasonally dry tropical forests in Africa, Asia, and America. The phaseoloid group provides an excellent case supporting the idea that the interplay of ecological opportunities and key innovations drives the evolutionary success.
Introduction
Species diversity of different lineages inhabiting the Earth is strikingly heterogeneous. The evolutionary success is responsible for those speciose lineages, but discerning what forces have driven the evolutionary success is a major challenge in evolutionary biology (Schluter, ). Recently, ecological opportunity and/or key innovation have been considered as a cause of high diversification rates in many groups (e.g., Yoder et al., ; Duputié et al., ; Claramunt et al., ; Erkens et al., ). However, little is known about how ecological opportunity and key innovation have interplayed to promote diversification of a group.
The origin and diversification of angiosperms since the Mesozoic are the great events in the plant kingdom. Angiosperms now include about 300,000 species, have taken various ecological habitats on the Earth, and supplied most raw materials for the well-being of human beings. The family Leguminosae is especially outstanding in this term. The phaseoloid legumes represent one of the largest clades in Leguminosae and consist of 114 genera with ca. 2000 species (Table S1; Lewis et al., ). This clade contains many economically important species, such as soybean (Glycine max), common bean (Phaseolus vulgaris), cowpea (Vigna unguiculata), pigeonpea (Cajanus cajan), horse gram (Macrotyloma uniflorum), siratro (Macroptilium atropurpureum), and coral tree (Erythrina spp.) (Bruneau and Doyle, ). Molecular phylogenetic studies have contributed greatly to the delimitation of the clade (Lavin et al., ; Hu et al., ; Lee and Hymowitz, ; Doyle et al., ; Wojciechowski et al., ; Lewis et al., ; Stefanović et al., ). Now, the phaseoloid legumes contain four subtribes of Phaseoleae, viz. Phaseolinae, Glycininae, Cajaninae, and Kennediinae, and Psoraleeae and Desmodieae (Lewis et al., ; Stefanović et al., ). The aforementioned phylogenetic analyses of the clade usually sampled less than 40% of the generic diversity in the group; thus, to understand the evolutionary dynamics responsible for its current diversity we need a well-resolved phylogeny of the phaseoloid genera with more comprehensive taxon sampling.
Based on an analysis of matK sequences with 13 fossil calibration points, Lavin et al. () suggest that the phaseoloid legumes became differentiated in the Oligocene (24.2–32.1 Ma). Since the Oligocene, global climate has undergone marked changes (Zachos et al., ; Wade and Pälike, ; Pälike et al., ), which may have impacted speciation of many organisms. Egan and Crandall () assume that the recent rapid radiation of Psoraleeae of the phaseoloid legumes may be due to global climate change during the Pleistocene. The phaseoloid legumes possess both woody and herbaceous habits and are primarily distributed in tropical and temperate forests or grasslands (Table S1; Sprent, ). The habit shifts have been suggested being responsible for the diversification of some angiosperm lineages (Tiffney and Mazer, ; Verdú, ; Jabbour and Renner, ). Nevertheless, it remains puzzling and unexplored how evolution of habits fostered diversification in the phaseoloid legumes and how ecological forces have been regulating cladogenesis in different geographical areas.
In this study, we first reconstruct a genus-level phylogeny for the phaseoloid legumes using three plastid loci with a more extensive generic sampling than in any previous studies. In the improved phylogenetic framework, we then explore how the interplay of past climate changes, geographical expansion and habit shifts may have triggered diversification of the phaseoloid legumes.
Materials and methods
Taxon sampling
We sampled 85 species from 82 of the 115 genera of the phaseoloid legumes. Our worldwide taxon sampling scheme covered all tribes and subtribes of Lewis et al. () and major clades of Stefanović et al. () in the phaseoloid legumes. Our outgroups included thirteen species, representing the other three subtribes of Phaseoleae, Diocleinae (four species), Clitoriinae (two species), and Ophrestiinae (one species), which are excluded from the phaseoloid legumes; Millettieae (five species) and Abreae (one species), following the results of Wojciechowski et al. (). Voucher information and GenBank accession numbers are listed in Table S2.
Laboratory protocols
Three chloroplast markers were used in this study: rbcL, trnL-F region (trnL intron, and trnL [UAA] 3′ exon-trnF [GAA] intergenic spacer), and trnK/matK region (including most of the 3′ flanking trnK intron and the entire matK gene).
Genomic DNA was isolated from silica-gel-dried materials using a Plant Genomic DNA Kit (Beijing Biomed Co., LTD, BJ, China) or from herbarium samples following a modified CTAB procedure (Doyle and Doyle, ). Three DNA regions were amplified with the polymerase chain reaction (PCR). The primers used in this study are listed in Table S3. PCR amplifications were performed using 2 × Taq PCR MasterMix (Beijing Biomed Co., LTD) in 25-μL reactions with the following thermocycler program: 2 min at 95°C for denaturation, then 35 cycles of 30 s at 95°C, 30–60 s at 53–57°C for annealing, 2 min 30 s at 72°C for primer extension, and a 10-min incubation at 72°C following the cycles. The PCR products were purified using a GFX™ PCR DNA and Gel Band Purification Kit (Amersham Pharmacia Biotech, Piscataway, NJ, USA) and then directly sequenced. Sequencing reactions were conducted using an ABI Prism BigDye Terminator Cycle Sequencing Kit (Applied Biosystems, ABI, BJ, China). Sequences were analyzed using ABI 3730 × l DNA Analysis Systems and following the manufacturer's protocols.
Phylogeny and divergence time estimates
Sequence alignments were done using CLUSTAL X v2.0 (Larkin et al., ) and manually adjusted with BioEdit v5.0.9 (Hall, ). All alignments are available upon request from the corresponding author. We used the Bayesian relaxed clock methodology as implemented in BEAST v1.7.5 (Drummond et al., ) to generate a dated phylogeny of the phaseoloid legumes. The GTR + I + Γ model was selected as the best-fit model for each plastid region determined by ModelTest v3.7 (Posada and Crandall, ). Base frequencies were estimated. Clock rate was estimated under an uncorrelated relaxed-clock log-normal (UCLN) model. A Yule speciation model was used as a prior on the tree.
Some fruit and leaf fossils of Pueraria of the phaseoloids from the Middle Miocene of middle latitudes in Asia have been described (Wang et al., , and references therein), but we did not use them as calibration points because Pueraria is not monophyletic (Lee and Hymowitz, ; Stefanović et al., ; this study). Employing 13 fossil age constraints imposed on the matK phylogeny, Lavin et al. () provide a credible age framework for Leguminosae. Following the results of Lavin et al. (), we selected six calibration points: (1) a 45.2 Ma constraint on the root age (node 1); (2) the split between Platycyamus regnellii and the phaseoloid legumes (node 2) was set to 39.7 Ma; (3) an age of 27.8 Ma to constrain the crown group age of the phaseoloid legumes (node 3); (4) the crown group age of Desmodieae (node 4) was set to 14.2 Ma; (5) the crown group age of clade VIII (node 5, Figure S1) was set to 19.2 Ma; and (6) an age of 6.3 Ma for the crown group age of tribe Psoraleeae (node 6). A normal distribution was used for all six calibration points. The standard deviation was set to contain the lower and higher boundaries of the 95% highest posterior density values. MCMC searches were run for 100,000,000 generations, sampled every 1000 generations. Tracer v1.5 was used to monitor appropriate burn-in and the adequate effective sample sizes of the posterior distribution (>200). The maximum clade credibility tree was computed by TreeAnnotator v1.7.5 in BEAST software package (Drummond et al., ). BEAST analyses were performed in the CIPRES Web Portal 3.1 (Miller et al., ).
Habit evolution
The reconstruction of habit evolution in the phaseoloid legumes were carried out using the parsimony method with Mesquite v2.74 (Maddison and Maddison, ). The maximum clade credibility tree obtained from BEAST was used in the analysis. Two habit states were scored, herbaceous (including herbs and herbaceous climbing vines) vs. woody (including trees, woody climbers, and shrubs), based on the literature (Table S2).
Diversification analyses
To visualize the temporal variation in diversification rates, semilogarithmic lineage-through-time (LTT) plots were constructed in the R package APE v2.5-1 (Paradis et al., ). To evaluate 95% credibility interval of the empirical LTT curve, 1000 ultrametric trees randomly sampled from the converged BEAST trees were also used to calculate semilogarithmic LTT plots.
To detect rapid shifts in diversification rates at any specified time, the RC statistic was calculated with the R package GEIGER v1.3-1 (Harmon et al., ). Lineages with more or fewer descendents than expected under the constant rate model were hypothesized as a diversification rate shift. Species diversity for the phaseoloid legumes was estimated from the number of species in each genus; missing genera were assigned to corresponding clades based on previous studies (Table S1). Net diversification rates (r) for the phaseoloid legumes, nodes R1 and R2 were calculated by using BEAST chronogram under two extremes of the relative extinction rate (ε = 0 and 0.9) following the whole-clade method (Magallón and Sanderson, ). Calculations were performed using the GEIGER v1.3-1 (Harmon et al., ).
Biogeographical analyses
To reconstruct the possible ancestral ranges of the phaseoloid legumes, we conducted a Bayes-DIVA analysis (Nylander et al., ) using the software package RASP (Yu et al., ). Bayes-DIVA method can minimize the phylogenetic uncertainties by utilizing the posterior distribution of trees resulting from a BEAST analysis and generating credibility support values for alternative phylogenetic relationships (Nylander et al., ; Yu et al., ). We randomly sampled 1000 trees from the BEAST output as a “trees file” and used the maximum clade credibility (MCC) tree as a final representative tree. Biogeographical analyses were conducted on continental spatial scale at generic level, because the aim was to predicate the ancestral areas of nodes deeper down into the tree other than the ancestral areas of individual genera. Six geographic regions were coded: A, Asia; B, Africa; C, Europe; D, Australia; E, South America; F, North America (including Central America and Caribbean). Ancestral areas were reconstructed with the “maxareas” constrained to 3 because 73 of the 82 genera occur in fewer than three areas.
Results
Phylogeny and divergence times
The maximum clade credibility tree generated by BEAST analyses is well-resolved (Figure S1). Within the phaseoloid legumes, eight major clades were recognized, and the Apios is the earliest-diverging lineage (PP = 0.97). Psoraleeae and Desmodieae are strongly supported as monophyletic, both of which are imbedded within Phaseoleae.
Molecular dating shows a stem age for the phaseoloid legumes of 39.5 Ma (35.7–43.2 Ma 95% highest posterior density, HPD) (Figure 1). The earliest diverged Apios lineage (clade I) separated with the remaining phaseoloids at 28.6 Ma (HPD: 25.8–31.2 Ma). The remaining phaseoloid splits into other seven clades (clade II to VII) between 26.8 and 20.4 Ma. The most recent common ancestor (MRCA) of Pueraria phaseoloides and Pueraria lobata emerged at ca. 13.4 Ma (HPD: 9.9–17.0 Ma). The Psoraleeae crown age is estimated at 6.1 Ma (HPD: 4.8–7.7 Ma).
Figure 1
Habit evolution
Results of ancestral habit state reconstruction are shown in Figure 1. The ancestral state of growth habit in the phaseoloid legumes is woody. Within the eight early diverged clades, six are woody (clade II–VII), while clade I and clade VIII are herbaceous. The herbaceous growth habit has evolved at least ten times within the phaseoloids. Importantly, clade VIII is the largest herbaceous lineage with some derived woody species, taking up ca. Thirty percentage species of the phaseoloid legumes.
Diversification rates
The semilogarithmic lineage-through-time (LTT) plots for taxa of whole phaseoloid legumes, woody clades (clade II–VII) and herbaceous clades (clade I and VIII) are shown in Figure 2. The whole phaseoloid legumes and woody clades showed a high diversification rate at the early stages (20.4–28.6 Ma). Herbaceous lineages experienced a high diversification rate since the Early Miocene.
Figure 2

Phaseoloid divergences through time, according to habit.
The relative cladogenesis (RC) statistic indicated two significant diversification rate shifts at node R1 and R2 within the phaseoloid legumes (Figure 1). The probability of R1 and R2 that had at least maximum descendents under the null hypothesis of a birth-death process are 0.024 and 0.010 separately (Table 1). Net diversification rates of the phaseoloid legumes were estimated as 0.13 (HPD: 0.12–0.15) net speciation events per million years (sp Myr−1) under a high relative extinction rate (ε = 0.9), and 0.19 sp Myr−1 under no extinction (ε = 0). Diversificaiton rates estimated for nodes R1 and R2 are shown in Table 2.
Table 1
Relative cladogenesis (RC) test of the phaseoloid legumes. Nodes with diversification rate shift are shown on Figure 1.
Table 2
| Clades | Stem group mean age (95%HPD) | No. species | r | |||
|---|---|---|---|---|---|---|
| ε= 0 | 95% HPD | ε = 0.9 | 95% HPD | |||
| Phaseoloied legumes | 39.5 (35.8–43.3) | 2005 | 0.19 | 0.18–0.21 | 0.13 | 0.12–0.15 |
| R1 | 28.6 (25.9–31.5) | 1998 | 1.05 | 1.05–1.05 | 1.05 | 1.05–1.05 |
| R2 | 21.9 (18.9–24.9) | 679 | 0.30 | 0.26–0.35 | 0.19 | 0.17–0.22 |
Net diversification rates (r) calculated for the phaseoloid legumes, nodes R1 and R2.
Biogeographical reconstruction
The results of ancestral area reconstruction using Bayes-DIVA in RASP is shown in Figure 3. The most recent common ancestor of the phaseoloid legumes is in Asia. Two independent intercontinental dispersal events occurred in the Late Oligocene. The first dispersal is to Africa with the rise of clade V, VI, VII, and VIII. The second dispersal is to Australia giving rise to clade II including Hardenbergia,Vandasina,Kennedia, and Shuteria. In the Miocene, the ancestral range of clade VIII expanded to South and North America following multiple dispersal events.
Figure 3

Biogeographical reconstruction of the phaseoloid legumes. The pie charts show the relative probabilities of alternative ancestral distributions obtained by Bayes-DIVA optimizations over the 1000 Bayesian trees (white > red). The first four areas with highest probability are colored according to relative probability in the following order: white > red > blue > gray; and the black portion represents reconstructions with a probability <0.10.
Discussion
Our divergence time estimates (Figure 1) suggest a stem age of 39.5 Ma (HPD: 35.7–43.2 Ma) for the phaseoloid legumes and a crown age of 28.6 Ma (HPD: 25.8–31.2 Ma), which are consistent with the estimates of Lavin et al. (
The phaseoloid legumes originated in the Late Eocene, but the group became differentiated in the Late Oligocene of Asia, and a dramatic accumulation of the phaseoloid lineages immediately occurred during the Late Oligocene and Early Miocene, with two dispersal events, from Asia to Africa and Australia. The rapid divergence time post-dates a period of a drastic global cooling resulting, in part, from the development of permanent continental ice-sheets in Antarctica (Zachos et al.,
LTT plots indicate that the dramatic accumulation of the herbaceous phaseoloid lineages (clade VIII) occurred since the Early Miocene (Figure 2). The RC test found the other shift of diversification rates within the haseoloid legumes, node R2 (clade VIII), whose net diversification rate is higher than that of the whole phaseoloid legumes (Table 2). Clade VIII is herbaceous and contains about 52 genera and 679 species, greatly contributing to phaseoloid diversity (Figure 1). Our biogeographical reconstruction suggests that the MRCA of clade VIII is in Africa and subsequently multiple independent migrations from Africa to Asia, North America, and South America occurred (Figure 3). The clade VIII became diversified in the Early Miocene (19.3 Ma; HPD: 17.0–21.6 Ma), which corresponds to the time when the African plate collided with the Eurasian one (ca. 18–17 Ma; Axelrod and Raven,
Statements
Author contributions
Wei Wang and Zhiduan Chen conceived the study. Honglei Li performed the experiments. Honglei Li, Wei Wang, Li Lin, and Xinyu Zhu analyzed the data. Wei Wang, Jianhua Li, Xiangyun Zhu, and Zhiduan Chen contributed reagents/materials/analysis tools. Honglei Li and Wei Wang wrote the paper. Wei Wang, Jianhua Li, and Zhiduan Chen edited the paper.
Acknowledgments
We sincerely thank Pamela S. Soltis for carefully reading an early draft of the manuscript. This research was supported by National Basic Research Program of China (No. 2014CB954101), National Natural Science Foundation of China (Nos.31270268, 31270269, and 30990241), and Natural Science Fund for Colleges and Universities in Jiangsu Province (No. 09KJB180006).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: http://www.frontiersin.org/journal/10.3389/fpls.2013.00386/abstract
References
1
AxelrodD. I.RavenP. H. (1978). Late cretaceous and tertiary vegetation history of Africa, in Biogeography and Ecology of Southern Africa, ed WergerM. J. A. (Netherlands: Springer), 77–130. 10.1007/978-94-009-9951-0_5
2
BruneauA.DoyleJ. J. (1990). A chloroplast DNA inversion as a subtribal character in the Phaseoleae (Leguminosae). Syst. Bot. 15, 378–386. 10.2307/2419351
3
BunnA.GoetzS. J.KimballJ. S.ZhangK. (2007). Northern high-latitude ecosystems respond to climate change. EOS Trans. Am. Geophys. Union88, 333–334. 10.1029/2007EO340001
4
ClaramuntS.DerryberryE.BrumfieldR.RemsenJ. (2012). Ecological opportunity and diversification in a continental radiation of birds: climbing adaptations and cladogenesis in the Furnariidae. Am. Nat. 179, 649–666. 10.1086/664998
5
CrispM. D.CookL. G. (2009). Explosive radiation or crypticmass extinction? Interpreting signatures in molecular phylogenies. Evolution63, 2257–2265. 10.1111/j.1558-5646.2009.00728.x
6
De-NovaJ. A.MedinaR.MonteroJ. C.WeeksA.RosellJ. A.OlsonM. E.et al. (2012). Insights into the historical construction of species-rich Mesoamerican seasonally dry tropical forests: the diversification of Bursera (Burseraceae, Sapindales). New. Phytol. 193, 276–287. 10.1111/j.1469-8137.2011.03909.x
7
DoddM. E.SilvertownJ.ChaseM. W. (1999). Phylogenetic analysis of trait evolution and species diversity variation among angiosperm families. Evolution53, 732–744. 10.2307/2640713
8
DoyleJ. J.DoyleJ. L. (1987). A rapid DNA isolation procedure for small quantities of fresh leaf tissue. Phytochem. Bull. 19, 11–15.
9
DoyleJ. J.DoyleJ. L.HarbisonC. (2003). Chloroplast-expressed glutamine synthetase in Glycine and related Leguminosae: phylogeny gene duplication and ancient polyploidy. Syst. Bot. 28, 567–577.
10
DrummondA. J.SuchardM. A.XieD.RambautA. (2012). Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol. Biol. Evol. 29, 1969–1973. 10.1093/molbev/mss075
11
DuputiéA.SalickJ.McKeyD. (2011). Evolutionary biogeography of Manihot (Euphorbiaceae), a rapidly radiating Neotropical genus restricted to dry environments. J. Biogeogr. 38, 1033–1043. 10.1111/j.1365-2699.2011.02474.x
12
EganA. N.CrandallK. A. (2008). Divergence and diversification in North American Psoraleeae (Fabaceae) due to climate change. BMC Biol. 6:55. 10.1186/1741-7007-6-55
13
ErikssonO.BremerB. (1992). Pollinationsystems dispersalmodes lifeforms and diversification rates in angiosperm families. Evolution46, 258–266. 10.2307/2409820
14
ErkensR. H. J.ChatrouL. W.CouvreurT. L. P. (2012). Radiations and key innovations in an early branching angiosperm lineage (Annonaceae; Magnoliales). Bot. J. Linn. Soc. 169, 117–134. 10.1111/j.1095-8339.2012.01223.x
15
ForbesB.FauriaM. M.ZetterbergP. (2010). Russian Arctic warming and greening are closely tracked by tundra shrub willows. Glob. Change Biol. 15, 1–13.
16
GrahamA. (2010). Late Cretaceous and Cenozoic History of Latin American Vegetation and Terrestrial Environments. StLouis, MO: Missouri Botanical Garden Press.
17
GuoZ. T.RuddimanW. F.HaoQ. Z.WuH. B.QiaoY. S.ZhuR. X.et al. (2002). Onset of Asian desertification by 22 Myr ago inferred from loess deposits in China. Nature416, 159–163. 10.1038/416159a
18
HallT. A. (1999). BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp. Ser. 41, 95–98.
19
HallingerM.MantheyM.WilmkingM. (2010). Establishing a missing link: warm summers and winter snow cover promote shrub expansion into alpine tundra in Scandinavia. New. Phytol. 186, 890–899. 10.1111/j.1469-8137.2010.03223.x
20
HallingerM.WilmkingM. (2011). No change without a cause - why climate change remains the most plausible reason for shrub growth dynamics in Scandinavia. New. Phytol. 189, 902–908. 10.1111/j.1469-8137.2010.03624.x
21
HarmonL. J.WeirJ. T.BrockC. D.GlorR. E.ChallengerW. (2008). GELGER: investigating evolutionary radiations. Bioinformatics24, 129–131. 10.1093/bioinformatics/btm538
22
HuJ. M.LavinM.WojciechowskiM. F.SandersonM. J. (2000). Phylogenetic systematics of the tribe Millettieae (Leguminosae) based on chloroplast trnK/matK sequences and its implications for evolutionary patterns in Papilionoideae. Am. J. Bot. 87, 418–430. 10.2307/2656638
23
JabbourF.RennerS. S. (2012). A phylogeny of Delphinieae (Ranunculaceae). shows that Aconitum is nested within Delphinium and that Late Miocene transitions to long life cycles in the Himalayas and Southwest China coincide with bursts in diversification. Mol. Phylogenet. Evol. 62, 928–942. 10.1016/j.ympev.2011.12.005
24
JacobsB. F. (2004). Palaeobotanical studies from tropical Africa: relevance to the evolution of forest, woodland and savannah biomes. Phil. Trans. R. Soc. Lond. B Biol. Sci. 359, 1573–1583. 10.1098/rstb.2004.1533
25
LantzT. C.KokeljS. V.GergelS. E.HenryG. H. R. (2009). Relative impacts of disturbance and temperature: persistent changes in microenvironment and vegetation in retrogressive thaw slumps. Glob. Change Biol. 15, 1664–1675. 10.1111/j.1365-2486.2009.01917.x
26
LarkinM. A.BlackshieldsG.BrownN. P.ChennaR.McGettiganP. A.McWilliamH.et al. (2007). Clustal W and Clustal X v2.0. Bioinformatics23, 2947–2948. 10.1093/bioinformatics/btm404
27
LavinM.EshbaughE.HuJ. M.MathewsS.SharrockR. A. (1998). Monophyletic subgroups of the tribe Millettieae (Leguminosae) as revealed by phytochrome nucleotide sequence data. Am. J. Bot. 85, 412–433. 10.2307/2446334
28
LavinM.HerendeenP.WojciechowskiM. (2005). Evolutionary rates analysis of Leguminosae implicates a rapid diversification of lineages during the Tertiary. Syst. Biol. 54, 575–594. 10.1080/10635150590947131
29
LeeJ.HymowitzT. (2001). A molecular phylogenetic study of the subtribe Glycininae (Leguminosae) derived from the chloroplast DNA rps16 intron sequences. Am. J. Bot. 88, 2064–2073. 10.2307/3558432
30
LewisG.SchrireB.MachinderB.LockM. (2005). Legumes of the World. Kew: The Royal Botanic Gardens.
31
MaddisonW. P.MaddisonD. R. (2009). Mesquite: a modular system for evolutionary analysis v2.74, in Computer Program and Documentation Distributed by the Author. Available online at: http://mesquiteprojectorg/mesquite/mesquitehtml.
32
MagallónS.SandersonM. J. (2001). Absolute diversification rates in angiosperm clades. Evolution55, 1762–1780. 10.1111/j.0014-3820.2001.tb00826.x
33
MillerM. A.PfeifferW.SchwartzT. (2010). Creating the CIPRES Science Gateway for inference of large phylogenetic trees, in Proceedings of the Gateway Computing Environments Workshop (GCE), (New Orleans, LA), 1–8. 10.1109/GCE.2010.5676129
34
MinnichR. A. (2007). Climate, Paleoclimate, and Paleovegetation, in Terestrial vegetation of California, 3rd Edn, eds BarbourA. G.Keeler-WolfT.SchoenherrA. A. (London: University of California Press), 43–70. 10.1525/california/9780520249554.003.0002
35
MosbruggerV. (2005). Cenozoic continental climatic evolution of central Europe. Proc. Natl. Acad. Sci. U.S.A. 102, 14964–14969. 10.1073/pnas.0505267102
36
NylanderJ. A. A.OlssonU.AlströmP.SanmartínI. (2008). Accounting for phylogenetic uncertainty in biogeography: a bayesian approach to dispersal-vicariance analysis of the thrushes (Aves: Turdus). Syst. Biol. 57, 257–268. 10.1080/10635150802044003
37
PälikeH.NorrisR. D.HerrleJ. O.WilsonP. A.CoxallH. K.LearC. H.et al. (2006). The heartbeat of the Oligocene climate system. Science314, 1894–1898. 10.1126/science.1133822
38
ParadisE.ClaudeJ.StrimmerK. (2004). APE: analyses of phylogenetics and evolution in R language. Bioinformatics20, 289–290. 10.1093/bioinformatics/btg412
39
PosadaD.CrandallK. A. (1998). MODELTEST: testing the model of DNA substitution. Bioinformatics14, 817–818. 10.1093/bioinformatics/14.9.817
40
ProtheroD. R. (1994). The late Eocene-Oligocene extinctions. Ann. Rev. Earth Planet Sci. 22, 145–165. 10.1146/annurev.ea.22.050194.001045
41
RaupD. M.SepkoskiJ. J. (1986). Periodic extinction of families and genera. Science231, 833–836. 10.1126/science.11542060
42
RejmánekM.RichardsonD. M. (2013). Trees and shrubs as invasive alien species - 2013 update of the global database. Divers. Distrib. 19, 1093–1094. 10.1111/ddi.12075
43
SchluterD. (2000). The Ecology of Adaptive Radiation. Oxford: Oxford University Press.
44
ShaverG. R.GiblinA. E.NadelhofferK. J.RastetterE. B. (1997). Plant functional types and ecosystem change in arctic tundra, in Plant Functional Types, eds SmithT.WoodwardI.ShugartH. (Cambridge: Cambridge University Press), 153–173.
45
SprentJ. I. (2007). Evolving ideas of legume evolution and diversity: a taxonomic perspective on the occurrence of nodulation. New Phytol. 174, 11–25. 10.1111/j.1469-8137.2007.02015.x
46
StefanovićS.PfeilB. E.PalmerJ. D.DoyelJ. J. (2009). Relationships among phaseoloid legumes based on sequences from eight chloroplast regions. Syst. Bot. 34, 115–128. 10.1600/036364409787602221
47
StrömbergC. A. E. (2011). Evolution of grasses and grassland ecosystems. Ann. Rev. Earth Planet Sci. 39, 517–544. 10.1146/annurev-earth-040809-152402
48
SturmM.SchimelJ.MichaelsonG.WelkerJ. M.OberbauerS. F.ListonG. E.et al. (2005). Winter biological processes could help convert Arctic tundra to shrubland. Bioscience55, 17–26. 10.1641/0006-3568(2005)055[0017:WBPCHC]2.0.CO;2
49
TapeK.SturmM.RacineC. (2006). The evidence for shrub expansion in Northern Alaska and the Pan-Arctic. Glob. Change Biol. 12, 686–702. 10.1111/j.1365-2486.2006.01128.x
50
TiffneyB. H.MazerS. J. (1995). Angiosperm growth habit dispersal and diversification reconsidered. Evol. Ecol. 9, 93–117. 10.1007/BF01237700
51
TöpelM.AntonelliA.YessonC.EriksenB. (2012). Past climate change and plant evolution in western North America: a case study in Rosaceae. PLoS ONE7:e50358. 10.1371/journal.pone.0050358
52
VerdúM. (2002). Age at maturity and diversification in woody angiosperms. Evolution56, 1352–1361. 10.1111/j.0014-3820.2002.tb01449.x
53
WadeB. S.PälikeH. (2004). Oligocene climate dynamics. Paleoceanography19, PA2019. 10.1029/2004PA001042
54
WangQ.ManchesterS. R.DilcherD. L. (2010). Fruits and foliage of Pueraria (Leguminosae, Papilionoideae) from the Neogene of Eurasia and their biogeographic implications. Am. J. Bot. 97, 1982–1998. 10.3732/ajb.1000167
55
WojciechowskiM. F.LavinM.SandersonM. J. (2004). A phylogeny of legumes (Leguminosae) based on analysis of the plastid matK gene resolves many well-supported subclades within the family. Am. J. Bot. 91, 1846–1862. 10.3732/ajb.91.11.1846
56
YoderJ. B.ClanceyE.Des RochesS.EastmanJ. M.GentryL.GodsoeW.et al. (2010). Ecological opportunity and the origin of adaptive radiations. J. Evol. Biol. 23, 1581–1596. 10.1111/j.1420-9101.2010.02029.x
57
YuY.HarrisA. J.HeX. (2011). RASP (Reconstruct Ancestral State in Phylogenies) 1.1. Available online at: http://mnh.scu.edu.cn/soft/blog/RASP
58
ZachosJ.PaganiM.SloanL.ThomasE.BillupsK. (2001). Trends rhythms and aberrations in global climate 65 Ma to present. Science292, 686–693. 10.1126/science.1059412
59
ZhangP.ChenY. Q.ZhouH.LiuY. F.WangX. L.PapenfussT. J.et al. (2006). Phylogeny, evolution, and biogeography of Asiatic Salamanders (Hynobiidae). Proc. Natl. Acad. Sci. U.S.A. 103, 7360–7365. 10.1073/pnas.0602325103
Summary
Keywords
aridification, biogeography, dispersal, diversification rate, habit shift, Leguminosae, molecular dating
Citation
Li H, Wang W, Lin L, Zhu X, Li J, Zhu X and Chen Z (2013) Diversification of the phaseoloid legumes: effects of climate change, range expansion and habit shift. Front. Plant Sci. 4:386. doi: 10.3389/fpls.2013.00386
Received
26 July 2013
Accepted
11 September 2013
Published
09 October 2013
Volume
4 - 2013
Edited by
Xin Wang, Chinese Academy of Sciences, China
Reviewed by
Yang Liu, University of Connecticut, USA; Guo-Qiang Zhang, The Orchid Conservation and Research Center of Shenzhen, China
Copyright
© 2013 Li, Wang, Lin, Zhu, Li, Zhu and Chen.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Wei Wang, State Key Laboratory of Systematic and Evolutionary Botany, Institute of Botany, Chinese Academy of Sciences, 20 Nanxincun, Xiangshan, Beijing 100093, China e-mail: wangwei1127@ibcas.ac.cn
This article was submitted to Plant Evolution and Development, a section of the journal Frontiers in Plant Science.
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.