On the historical biogeography of global Galliformes：ancestral range and diversification patterns

Youhua Chen

On the historical biogeography of global Galliformes：ancestral range and diversification patterns

Youhua Chen

Background:In this study,the ancestral distributional ranges and the tempo of diversification patterns of global Galliformes were investigated.

Game birds,Macro-evolution,Rate shifting,Zoogeography,Phylogenetics

Background

Galliformes are ground-living birds,an order consisting of around 280 species worldwide(Chen,2014).This group of birds is found to have around 23 species listed as endangered and 6 as critically endangered in the IUCN red list(The World Conservation Union 2010). Phylogenetic relationships among species have been widely studied in recent years(Crowe et al.,2006; Frank-Hoeflich et al.,2007;Hugall and Stuart-Fox, 2012).Although current phylogenetic information is still limited and does not cover the systematics and affinities of all 280 Galliformes,a recent investigation has established a phylogeny of Galliformes for up to 197 species(Hugall and Stuart-Fox,2012),accounting for over 70%of global Galliformes.

Macro-evolutionary patterns of this large-size and attractive-appearance bird assemblage might be initiated by utilizing the available phylogenetic affinity information of the 197 species(Hugall and Stuart-Fox,2012). The proposal to study macro-evolutionary patterns of Galliformes is to reveal in a more comprehensive way the extinction mechanism of this species assemblage from a long-term evolutionary perspective,for the purpose of better conserving them.In the present study, several macro-evolutionary attributes relevant to the phylogeny and diversity of Galliformes are considered.

In first instance,clade age has been thought to link up with species richness.The relationship between clade age and species richness is one of recent interest in macroevolutionary studies(Rabosky et al.,2012;Rabosky andAdams,2012).Clade age has been thought to relate to species richness because of the fact that older clades could have more time to diversify(McPeek and Brown,2007; Etienne et al.,2012),implying that the older the age of the clade,the higher its richness,resulting in a positive clade age-species richness relationship.However,whether such a relationship is universal is still controversial.Several previous studies have shown that clade age could predict species richness(McPeek and Brown,2007;Etienne et al., 2012),while others argued that there is no clear relationship between clade age and species diversity(Rabosky et al.,2012).In the present study,I wanted to test whether there is a relationship between the evolutionary age of the ancestors of Galliformes and the number of their externally living descendants.

Secondly,a shifting pattern in the rate of diversification has been broadly observed in many taxa,by showing patterns in which this rate is initially high but decreases over time(Rabosky and Lovette,2008a,2008b).Such a declining-diversification model is predicted by adaptive evolution(Schluter,2000),because openings of new vacant niches are limited.Filling vacant niches would have been accomplished at an early evolutionary time,leading to a decline in the rate of diversification of species (Rabosky and Lovette,2008b).

Finally,the size of their range might be related to the rate of diversification of species(Chen,2013),since range size of a species is a trait jointly affected by species dispersal,colonization and reproduction.The importance of geographic isolation in shaping speciation has been debated for a long time(Fitzpatrick and Turelli, 2006).A phylogenetic comparative method,referred to as“age-range correlation”,is introduced to quantify the relative importance of sympatric and allopatric speciation on structuring contemporary species diversity patterns(Fitzpatrick and Turelli,2006;Bolnick and Fitzpatrick,2007).

Considerable progress has been made in the development of robust statistical methods to infer ancestral ranges by incorporating a variety of biological processes.For example,a dispersal-extinction-cladogenesis model(DEC)has been proposed to estimate explicitly and infer historical changing patterns of ancestral ranges of species when projected on the phylogeny under a maximum likelihood framework (Ree et al.,2005;Ree and Smith,2008).In addition to maximum likelihood-based methods for reconstructing ancestral distributional ranges,traditional methods are derived from the parsimony principle.Dispersal-vicariance analysis (DIVA)(Ronquist,1997)and some of its extensions(for example,statistical DIVA(S-DIVA),Yu et al.,2010)are built on parsimony algorithms and still widely cited in current literature of phylogeographic studies(Ali et al.,2012).Recently,a Markov Chain Monte Carlo(MCMC)method has been proposed by Yu et al.(2013).

Methods

Historical biogeographical analyses

The ancestral ranges of the Galliformes group were estimated by utilizing the well-established phylogenetic tree for 197 Galliformes(Hugall and Stuart-Fox,2012).The following terrestrial regions are considered in biogeographical analyses:East Asia(A),South Asia(B),Southeast Asia (C),West Asia(D),North America(E),South America(F), Africa(G),Europe(H)and Oceania(I).Distribution of each species over these terrestrial regions was obtained from the Avibase database(http://avibase.bsc-eoc.org/ avibase.jsp?lang=EN).

In comparing ancestral ranges of Galliformes species, three analytical methods are used,i.e.,dispersal-vicariance analysis(DIVA),Bayesian binary MCMC analysis(BBM) and dispersal-extinction-cladogenesis analysis(DEC).All three methods are carried out by using the software RASP (Yu et al.,2010,2011).

Temporal analyses of rates of diversification

Different diversification rate-shifting models have been tried to fit the Galliformes phylogeny,as proposed in previous studies(Rabosky,2006b;Rabosky and Lovette, 2008a,2008b).Specifically,four methods from an R package“laser”(Rabosky,2006a)are implemented for comparative purposes,consisting of a constant-speciation and constant-extinction model(CONSTANT),a decreasingspeciation and constant-extinction model(SPVAR),a constant-speciation and increasing-extinction model (EXVAR)and a decreasing-speciation and increasingextinction model(BOTHVAR)(Chen,2013).These four models have been used to test the temporal shifting patterns in rates of diversification of different taxa(Rabosky and Lovette,2008a,2008b).

Each of the models require four parameters for estimation(Rabosky and Lovette,2008b),which can be obtained by maximizing the following likelihood equation (Rabosky,2006b):

wheretis the vector of observed branch times from the phylogeny,tnthe branch time for the lineage,whilen,λ(t)and μ(t)are time-dependent speciation and extinction rates,respectively.The time-dependent rate of diversification is defined asr(t)=λ(t)-μ(t)andNis the number of external tips in the tree.

In the CONSTANT model,λ(t)and μ(t)are assumed to be constant over the entire phylogenetic tree(i.e., λ(t)=λ0,μ(t)=μ0),where λ0and μ0are the constants to be estimated.In the SPVAR model,μ(t)is assumed to be constant over the entire tree(i.e.,μ(t)=μ0),while λ(t)is assumed to decrease continuously from the root to the tips of the tree,defined as follows:λ(t)=λ0exp(-kt).As seen in the SPVAR model,the additional parameterkis required to model the declining trend of rate of speciation over the tree.For this model,the rate of diversification is predicted to decline over the evolutionary time as(r(t)=λ0exp(-kt)-μ0).In the EXVAR model,the rate of speciation is assumed to be constant over the tree while the rate of extinction is assumed to decline over the tree as follows:u(t)=u0(1-exp(-zt)).As well,this model has an additional parameter,i.e.,z,to be estimated.Finally,in the BOTHVAR model,both λ(t)and μ(t) are assumed to change over time as follows:λ(t)=λ0exp (-kt)andu(t)=u0(1-exp(-zt))(Rabosky,2006b;Rabosky and Lovette,2008b).

Clade age and clade richness relationships

The evolutionary age for each clade is calculated as the phylogenetic distance between the root and the internal node leading to the focused clade.The corresponding clade species richness is defined as the number of external tips(living species)for that specific clade.

To reveal the possible relationship between clade age and clade species richness and/or phylogenetic diversity, I performed both a non-phylogenetic ordinary leastsquares regression analysis(OLS)and a phylogenetic general least-squares regression analysis(PGLS).Given that various clades are not independent from each other, it is necessary to remove the impacts of phylogenetic inertia by performing PGLS,i.e.,a method to introduce a phylogenetic variance-covariance matrix in the fitting formula,a matrix missing in the OLS method.For the OLS method,the vector of coefficients is fitted using the following identity:

while for the PGLS method,the vector of coefficients is estimated from the following equation:

The superscriptTdenotes the transpose of a matrix, while-1 denotes the inverse of a matrix.X is a matrix with columns indicating the explanatory variables,whileyis a column vector storing the values for the response variable and W is the phylogenetic variance-covariance matrix.The calculation of W is only related to the branch lengths of the phylogenetic tree(Revell,2010).

Results

Historical dispersal,vicariance and extinction events

Both the Bayesian and maximum likelihood methods (BBI and Lagrange)identified SE Asia as the most likely origin of the distribution of the most common ancestor for all Galliformes species.The S-DIVA method failed to run because of unknown errors(out of memory when using RASP software).

Because BBI and Lagrange share some levels of similarity,ancestral ranges estimated by Lagrange are the focus in the subsequent biogoegraphic discussion(Figure 1).As seen,a number of dispersal and vicariance events have occurred in the distribution of Galliformes species.

At root node 393(Figure 1),the most likely ancestral range is SE Asia and N America with a marginal probability of 29%(light orange color,symbol:CE).One vicariance event is identified,in which the Galliformes lineages in N America and SE Asia are separated in subsequent evolutionary times.

For the lineages distributed in N America,in later times at node 378(Figure 1),two dispersal events occurred.Starting from N America,one ancestral lineage dispersed to S America(node 233,Figure 1,light blue color,symbol:EF)while another dispersed to Africa (node 377,Figure 1,purple color,symbol:EG).Again some time later(node 366),some lineages dispersed to E Asia(blue color,symbol:A).

For the lineage distributed in SE Asia at the root (Figure 1),it continued to inhabit that region and later underwent local radiation up until our contemporary era.During some evolutionary time points at nodes 391 and 389,some ancestral lineages of Galliformes dispersed to Oceania(green color,symbol:I),leading to the contemporary distribution ofAlectura lathamiandLeipoa ocellatain Australia,Megapodius pritchardiiin Tonga andMegapodius layardiin Vanuatu.

Diversification rate-shifting patterns during evolutionary history

No significant diversification rate-shifting pattern is evident in Galliformes phylogeny.As shown by the comparison of different diversification models,the constant-rate model received the lowest AIC value(-621.62,Table 1)and thus became the best model.

Clade age-richness relationship

There is a significant and positive relationship between clade age and Galliformes species richness(Figure 2), regardless of whether the situation or whether phylogenetic inertia is controlled or not.For the OLS,the fitted equation ofclade richness=1.24×clade age-8.29 (p＜0.05),while for the PGLS,the best fitted equation of clade richness=3.111×clade age-101.48(p＜0.05). This supports the prediction that older clades processhigher species diversity since they have more time to diversify.

Figure 1 Ancestral range reconstruction of Galliformes with color legends using Lagrange maximum likelihood method.Codes for terrestrial regions:E Asia(A),S Asia(B),SE Asia(C),W Asia(D),N America (E),S America(F),Africa(G),Europe(H)and Oceania(I).

Discussion

Interestingly,no significant diversification-shifting trend has been observed for the phylogeny of the 197 Galliformes species(Table 1).Several previous studies working on other avian taxa suggest a temporally diversification ratedeclining pattern,for example that of the North American Wood-warblers(Rabosky and Lovette,2008a).However, other studies also show that rates of diversification could have shown a temporally increasing trend up to our contemporary era for some specific species assemblages,for example Tiger Beetles(Barraclough and Vogler,2002)and angiosperm plants(Magallon and Castillo,2008).It is still controversial whether the rate of species diversification has a density-dependent trend throughout evolutionary times.From my observation of the Galliformes phylogeny, I conclude that this avian assemblage has a relatively constant rate of diversification over time,contradicting any time-dependent shifting trends of diversification.

There is a strong correlation between clade age and Galliformes richness as shown in Figure 2.As such,the currentstudy confirms thatclade age predicts Galliformes diversity throughout evolutionary times, but not rates of diversification.This conclusion is consistent with those of several previous studies (McPeek and Brown,2007;Etienne et al.,2012).However, other studies have argued that there is no clear positive relationship between clade age and species richness(Rabosky et al.,2012;Rabosky,2009)due to variation in rates of diversification among clades(Rabosky et al.,2007).Because I found that rates of diversification tend to be constant over evolutionary times for this Galliformes assemblage(Table 1), species richness of this avian group is principally driven by clade age(Figure 2).

As evidenced by the results of the Lagrange maximum likelihood analysis,it is found that SE Asia and N America are two disjunctive ancestral distributional origins for the earliest ancestor of Galliformes species.Dispersal frequency is very high for ancestral lineages of Galliformes.At some point in time,the lineage distributed in SE Asia then dispersed over Oceania while another lineage from N America dispersed to S America,Africa,E Asia and Europe.Active dispersals of Galliformes ancestors over the various continents might be an important driver of species diversity,because new vacant niches were available in these new terrestrial regions.Therefore,instead of rates of diversification,ecological opportunity might have played,implicitly,a role in species richness of Galliformes (Rabosky,2009;Mahler et al.,2010;Setiadi et al.,2011).

Table 1 A comparison of different time-dependent diversification models for fitting the Galliformes phylogeny

Some limitations apply to this study.First,the phylogeny used contains only 197 Galliformes species, which might be not sufficient to unravel the true macro-evolutionary pattern of the assemblage because the remaining 83 species have not been included in this analysis.It has been suggested that estimating rates of diversification is very sensitive to the complete status of the tree,because incomplete taxonomic sampling could generate artificial diversification rate-shifting patterns(Cusimano and Renner,2010;Brock et al.,2011). Second,estimation of the ancestral range is widely applied in plant biogeographical studies.Whether it is legitimate to infer historically ancestral ranges for bird taxa using planttailored statistical methods requires further elaboration. However,there is a growing trend in inferring ancestral life-history states of avian groups(Winger et al.,2012).To a certain extent therefore,it should be rational to estimate the ancestral distribution of Galliformes using the analytical DEC,S-DIVA and BBM methods.Third,utilization of country-level distributional information might not be sufficient to quantify range-clade relationships,given that country-level distributional records do not accurately reflect the true distributional ranges of Galliformes species. For example,some species might be present in a small area of a large country,leading to the over-representation of the distribution of the species.

In implementing further studies,the correlation of life-history traits and rates of diversification would be of interest,since it is reported that the ability of migration of avian species might slow down the rate of speciation of taxa(Ikeda et al.,2012;Claramunt et al.,2012).Analyses of the correlation between functional traits and rates of diversification would offer some insights intothe relationship between dispersal capability and speciation patterns for Galliformes(Claramunt et al.,2012).

Figure 2 Clade age-richness relationships.The black fitted line indicated the result from ordinary least-squares regression analysis(OLS), while the red line is from the phylogenetical general least-squares regression analysis(PGLS).Both lines have significant slopes(p＜0.05).For OLS, the fitted equation is clade richness=1.24×clade age-8.29;while for PGLS,the fitted equation is clade richness=3.111×clade age-101.48.

Conclusions

The constant diversification rate for global Galliforme species implied that there were no diversification rateshifting trends for Galliformes species.The present study may contribute to the understanding of the ecology and diversity patterns of Galliformes from the perspective ofhistoricalbiogeography,although some limitations existed.

Competing interests

The author declares that he has no competing interests.

Acknowledgements

This work is supported by the China Scholarship Council(CSC).I like to thank two anonymous reviewers for their insightful comments to improve the quality of the present work.

Received：31 July 2014 Accepted：6 August 2014

Ali S,Yan Y,Pfosser M,Wetschnig W(2012)Inferences of biogeographical histories within subfamily Hyacinthoideae using S-DIVA and Bayesian binary MCMC analysis implemented in RASP(Reconstruct Ancestral State in Phylogenies).Ann Bot 109:95-107

Barraclough T,Vogler A(2002)Recent diversification rates in North American Tiger beetles estimated from a dated mtDNA phylogenetic tree.Mol Biol Evol 19:1706-1716

Bolnick D,Fitzpatrick B(2007)Sympatric speciation:models and empirical evidence.Annu Rev Ecol Syst 38:459-487

Brock C,Harmon L,Alfaro M(2011)Testing for temporal variation in diversification rates when sampling is incomplete and nonrandom.Syst Biol 60:410-419

Chen Y(2014)Conservation priority of global Galliformes species based on phylogenetic diversity.Integr Zool 9:340-348.

Chen Y(2013)A phylogenetic subclade analysis of range sizes of endemic woody see plant species of China:trait conservatism,diversification rates and evolutionary models.J Syst Evol 51:590-600

Claramunt S,Derryberry E,Remsen J,Brumfield R(2012)High dispersal ability inhabits speciation in continental raidation of passerine birds.Proc R Soc B Biol Sci 279:1567-1574

Crowe T,Bowie R,Bloomer P,Mandiwana T,Hedderson T,Randi E,Pereira S, Wakeling J(2006)Phylogenetics,biogeography and classification of,and character evolution in gamebirds(Aves:Galliformes):effects of character exclusion,data partitioning and missing data.Cladistics 22:495-532

Cusimano N,Renner S(2010)Slowdowns in diverisfication rates from real phylogenies may not be real.Syst Biol 59:458-464

Etienne R,de Visser S,Janzen T,Olsen J,Olff H,Rosindell J(2012)Can clade age alone explain the relationship between body size and diversity?Interface Focus 2:170-179

Fitzpatrick B,Turelli M(2006)The geography of mammalian speciation:mixed signals from phylogenies and range maps.Evolution(N Y)60:601-615

Frank-Hoeflich K,Silveira L,Estudillo-Lopez J,Garcia-Koch A,Ongay-Larios L, Pinero D(2007)Increased taxon and character sampling reveals novel intergeneric relationships in the Cracidae(Aves:Galliformes).J Zool Syst Evol Res 45:242-254

Hugall A,Stuart-Fox D(2012)Accelerated speciation in colour-polymorphic birds. Nature 485:631-634

Ikeda H,Nishikawa M,Sota T(2012)Loss of flight promotes beetle diversification. Nat Commun 3:648

Magallon S,Castillo A(2008)Angiosperm diversification through time.Am J Bot 96:349-365

Mahler D,Revell L,Glor R,Losos J(2010)Ecological opportunity and the rate of morphological evolution in the diversification of greater antillean anoles. Evolution(N Y)64:2731-2745

McPeek M,Brown J(2007)Clade age and not diversification rate explains species richenss among animal taxa.Am Nat 169:E97-E106

Rabosky D(2006a)LASER:a maximum likelihood toolkit for detecting temporal shifts in diversification rates from molecular phylogenies.Evol Bioinform Online 2:247-250

Rabosky D(2006b)Likelihood methods for detecting temporal shifts in diversification rates.Evolution(N Y)60:1152-1164

Rabosky D,Adams D(2012)Rates of morphological evolution are correlated with species richness in salamanders.Evolution(N Y)66:1807-1818

Rabosky D,Donnellan S,Talaba A,Lovette I(2007)Exceptional among-lineage variation in diversification rates during the radiation of Australia’s most diverse vertebrate clade.Proc R Soc B Biol Sci 274:2915-2923

Rabosky D,Lovette I(2008a)Density dependent diversification in North American wood-warblers.Proc R Soc B Biol Sci 275:2363-2371

Rabosky D,Lovette I(2008b)Explosive evolutionary radiations:decreasing speciation or increasing extinction through time?Evolution(N Y)62:1866-1875

Rabosky D,Slater G,Alfaro M(2012)Clade age and species richness are decoupled across the eukaryotic tree of life.PLoS Biol 10:e1001381

Rabosky DL(2009)Ecological limits and diversification rate:alternative paradigms to explain the variation in species richness among clades and regions. Ecol Lett 12:35-43

Ree R,Moore B,Webb C,Donoghue M(2005)A likelihood framework for inferring the evolution of geogrpahic range on phylogenetic trees.Evolution (N Y)59:2299-2311

Ree R,Smith S(2008)Maximum likelihood inference of geographic range evolution by dispersal,local extinction and cladogenesis.Syst Biol 57:4-14

Revell L(2010)Phylogenetic signal and linear regression on species data. Methods Ecol Evol 1:319-329

Ronquist F(1997)Dispersal-vicariance analysis:a new approach to the quantification of historical biogeography.Syst Biol 46:195-203

Schluter D(2000)The Ecology of Adaptive Radiation.Oxford University Press, New York

Setiadi M,McGuire J,Brown R,Zubairi M,Iskandar D,Andayani N,Supriatna J, Evans B(2011)Adaptive radiation and ecological opportunity in Sulawesi and Philippine fanged frog(Limnonectes)communities.Am Nat 178:221-240

The World Conservation Union(2010)IUCN Red List of Threatened Species. [http://www.iucnredlist.org/]

Winger B,Lovette I,Winkler D(2012)Ancestry and evolution of seasonal migration in the Parulidae.Proc R Soc B Biol Sci 279:610-618

Yu Y,Harris A,He X(2010)S-DIVA(statistical dispersal-vicariance analysis):a tool for inferring biogeographic histories.Mol Phylogenet Evol 56:848-850

Yu Y,Harris A,He X-J(2013)RASP(Reconstruct Ancestral State in Phylogenetics) 2.1.beta.[http://mnh.scu.edu.cn/soft/blog/RASP]

Yu Y,Harris A,He X-J(2011)RASP(Reconstruct Ancestral State in Phylogenies) 2.0 beta.[http://mnh.scu.edu.cn/soft/blog/RASP]

Cite this article as:Chen:On the historical biogeography of global

Galliformes：ancestral range and diversification patterns.Avian Research 2014 5:3.

10.1186/s40657-014-0003-9

Correspondence:haydi@126.com

Department of Renewable Resources,University of Alberta,Edmonton T6G 2H1,Canada

©2014 Chen;licensee BioMed Central Ltd.This is an Open Access article distributed under the terms of the Creative

Commons Attribution License(http://creativecommons.org/licenses/by/4.0),which permits unrestricted use,distribution,and reproduction in any medium,provided the original work is properly credited.The Creative Commons Public Domain Dedication waiver(http://creativecommons.org/publicdomain/zero/1.0/)applies to the data made available in this article, unless otherwise stated.

Methods:Different diversification models characterizing possible tempo patterns were fitted and compared onto the phylogenetic tree for the 197 Galliforme species,consisting of a constant-speciation and constant-extinction model(CONSTANT),a decreasing-speciation and constant-extinction model(SPVAR),a constant-speciation and increasing-extinction model(EXVAR)and a decreasing-speciation and increasing-extinction model(BOTHVAR). Ancestral range reconstruction was conducted using the dispersal-extinction-cladogenesis model.

Results:A constant-diversification-rate(CONSTANT)model best quantified the historical speciation patterns of this avian assemblage through model selection.Clade age and species richness are significantly and positively correlated.The most recent common ancestor for Galliformes species was originally found in the disjunctive regions between Southeast Asia and North America.High-frequency dispersal events were identified across the whole evolutionary time.

Conclusions:The constant diversification rate for global Galliforme species implied that there were no diversification rate-shifting trends for Galliformes species.The present study may contribute to the understanding of the ecology and diversity patterns of Galliformes from the perspective of historical biogeography,although some limitations existed.