Bayesian inference of species trees from multilocus data

Evolutionary patterns of diversification in the Andean hummingbird genus Adelomyia

The patterns of genetic diversity and morphological variation are of central importance in understanding the evolutionary process that drive diversification. We use molecular, morphological, and ecological data to explore the influence of geography and ecology in promoting speciation in the widespread Andean hummingbird genus Adelomyia. Six monophyletic clades were recovered which show distributional limits at well-defined geographic barriers. Percentage sequence divergence ranged between 5.8% and 8.2% between phylogroups separated by large (>4000 km) and small (<50 km) distances respectively, suggesting that geographic isolation may be influential at very different scales. We show that morphological traits in independent phylogroups are more related to environmental heterogeneity than to geographic barriers. We provide a molecular reconstruction of relationships within Adelomyia and recommend its use in future comparative studies of historical biogeography and diversification in the Andes.

Diversification across an altitudinal gradient in the Tiny Greenbul (Phyllastrephus debilis) from the Eastern Arc Mountains of Africa

BACKGROUND

The Eastern Arc Mountains of Africa have become one of the focal systems with which to explore the patterns and mechanisms of diversification among montane species and populations. One unresolved question is the extent to which populations inhabiting montane forest interact with those of adjacent lowland forest abutting the coast of eastern Africa. The Tiny Greenbul (Phyllastephus debilis) represents the only described bird species within the Eastern Arc/coastal forest mosaic, which is polytypic across an altitudinal gradient: the subspecies albigula (green head) is distributed in the montane Usambara and Nguru Mountains whereas the subspecies rabai (grey head) is found in Tanzanian lowland and foothill forest. Using a combination of morphological and genetic data, we aim to establish if the pattern of morphological differentiation in the Tiny Greenbul (Phyllastrephus debilis) is the result of disruptive selection along an altitudinal gradient or a consequence of secondary contact following population expansion of two differentiated lineages.

RESULTS

We found significant biometric differences between the lowland (rabai) and montane (albigula) populations in Tanzania. The differences in shape are coupled with discrete differences in the coloration of the underparts. Using multi-locus data gathered from 124 individuals, we show that lowland and montane birds form two distinct genetic lineages. The divergence between the two forms occurred between 2.4 and 3.1 Myrs ago.Our coalescent analyses suggest that limited gene flow, mostly from the subspecies rabai to albigula, is taking place at three mid-altitude localities, where lowland and montane rainforest directly abut. The extent of this introgression appears to be limited and is likely a consequence of the recent expansion of rabai further inland.

CONCLUSION

The clear altitudinal segregation in morphology found within the Tiny Greenbul is the result of secondary contact of two highly differentiated lineages rather than disruptive selection in plumage pattern across an altitudinal gradient. Based on our results, we recommend albigula be elevated to species rank.

Molecular evidence for a single evolutionary origin of domesticated rice

Asian rice, Oryza sativa, is one of world's oldest and most important crop species. Rice is believed to have been domesticated ∼9,000 y ago, although debate on its origin remains contentious. A single-origin model suggests that two main subspecies of Asian rice, indica and japonica, were domesticated from the wild rice O. rufipogon. In contrast, the multiple independent domestication model proposes that these two major rice types were domesticated separately and in different parts of the species range of wild rice. This latter view has gained much support from the observation of strong genetic differentiation between indica and japonica as well as several phylogenetic studies of rice domestication. We reexamine the evolutionary history of domesticated rice by resequencing 630 gene fragments on chromosomes 8, 10, and 12 from a diverse set of wild and domesticated rice accessions. Using patterns of SNPs, we identify 20 putative selective sweeps on these chromosomes in cultivated rice. Demographic modeling based on these SNP data and a diffusion-based approach provide the strongest support for a single domestication origin of rice. Bayesian phylogenetic analyses implementing the multispecies coalescent and using previously published phylogenetic sequence datasets also point to a single origin of Asian domesticated rice. Finally, we date the origin of domestication at ∼8,200-13,500 y ago, depending on the molecular clock estimate that is used, which is consistent with known archaeological data that suggests rice was first cultivated at around this time in the Yangtze Valley of China.

Unprecedented within-species chromosome number cline in the Wood White butterfly Leptidea sinapis and its significance for karyotype evolution and speciation

Background

Species generally have a fixed number of chromosomes in the cell nuclei while between-species differences are common and often pronounced. These differences could have evolved through multiple speciation events, each involving the fixation of a single chromosomal rearrangement. Alternatively, marked changes in the karyotype may be the consequence of within-species accumulation of multiple chromosomal fissions/fusions, resulting in highly polymorphic systems with the subsequent extinction of intermediate karyomorphs. Although this mechanism of chromosome number evolution is possible in theory, it has not been well documented.

Results

We present the discovery of exceptional intraspecific variability in the karyotype of the widespread Eurasian butterfly Leptidea sinapis. We show that within this species the diploid chromosome number gradually decreases from 2n = 106 in Spain to 2n = 56 in eastern Kazakhstan, resulting in a 6000 km-wide cline that originated recently (8,500 to 31,000 years ago). Remarkably, intrapopulational chromosome number polymorphism exists, the chromosome number range overlaps between some populations separated by hundreds of kilometers, and chromosomal heterozygotes are abundant. We demonstrate that this karyotypic variability is intraspecific because in L. sinapis a broad geographical distribution is coupled with a homogenous morphological and genetic structure.

Conclusions

The discovered system represents the first clearly documented case of explosive chromosome number evolution through intraspecific and intrapopulation accumulation of multiple chromosomal changes. Leptidea sinapis may be used as a model system for studying speciation by means of chromosomally-based suppressed recombination mechanisms, as well as clinal speciation, a process that is theoretically possible but difficult to document. The discovered cline seems to represent a narrow time-window of the very first steps of species formation linked to multiple chromosomal changes that have occurred explosively. This case offers a rare opportunity to study this process before drift, dispersal, selection, extinction and speciation erase the traces of microevolutionary events and just leave the final picture of a pronounced interspecific chromosomal difference.

Estimating phylogenetic relationships despite discordant gene trees across loci: the species tree of a diverse species group of feather mites (Acari: Proctophyllodidae)

With the increased availability of multilocus sequence data, the lack of concordance of gene trees estimated for independent loci has focused attention on both the biological processes producing the discord and the methodologies used to estimate phylogenetic relationships. What has emerged is a suite of new analytical tools for phylogenetic inference--species tree approaches. In contrast to traditional phylogenetic methods that are stymied by the idiosyncrasies of gene trees, approaches for estimating species trees explicitly take into account the cause of discord among loci and, in the process, provides a direct estimate of phylogenetic history (i.e. the history of species divergence, not divergence of specific loci). We illustrate the utility of species tree estimates with an analysis of a diverse group of feather mites, the pinnatus species group (genus Proctophyllodes). Discord among four sequenced nuclear loci is consistent with theoretical expectations, given the short time separating speciation events (as evident by short internodes relative to terminal branch lengths in the trees). Nevertheless, many of the relationships are well resolved in a Bayesian estimate of the species tree; the analysis also highlights ambiguous aspects of the phylogeny that require additional loci. The broad utility of species tree approaches is discussed, and specifically, their application to groups with high speciation rates--a history of diversification with particular prevalence in host/parasite systems where species interactions can drive rapid diversification.

A single Mid-Pleistocene long-distance dispersal by a bird can explain the extreme bipolar disjunction in crowberries (Empetrum)

The proposed age of the striking biogeographic disjunction between the Arctic and southernmost South America varies from more than 65 million to a few thousand years, but no estimates based on explicit models and molecular data are available. Here we address the origin of bipolarity in crowberries (Empetrum), which are heath-forming dwarf shrubs with animal-dispersed fruits. We apply a fossil-calibrated relaxed molecular clock to model sequence evolution in two nuclear low-copy and two plastid DNA regions from 41 individual plants (420 clones for the nuclear regions) representing the entire geographic distribution of crowberries. The plastid region matK and four fossil calibration points were used to infer the ages of the crowberry stem and crown groups. All analyses resolved three major crowberry clades (A-C). Clade A contained sequences from the eastern Canadian pink-fruited crowberry (E. eamesii) as sister to clades B and C, which both contained sequences from the black-fruited northern hemisphere crowberry (E. nigrum). Clade B also contained a subclade with all sequences from the red-fruited southern hemisphere crowberry, which is often referred to as a distinct species, E. rubrum. Its closest relatives were consistently identified as black-fruited plants from northwestern North America. The median time to the most recent common ancestor for northern and southern hemisphere crowberries was estimated to 0.56-0.93 Ma, and 0.26-0.59 Ma for the southern plants only. We conclude that a single dispersal by a bird from northwestern North America to southernmost South America, taking place in the Mid-Pleistocene, is sufficient to explain the disjunction in crowberries.

Genetic population structure of the alpine species Rhododendron pseudochrysanthum sensu lato (Ericaceae) inferred from chloroplast and nuclear DNA

Background

A complex of incipient species with different degrees of morphological or ecological differentiation provides an ideal model for studying species divergence. We examined the phylogeography and the evolutionary history of the Rhododendron pseudochrysanthum s. l.

Results

Systematic inconsistency was detected between gene genealogies of the cpDNA and nrDNA. Rooted at R. hyperythrum and R. formosana, both trees lacked reciprocal monophyly for all members of the complex. For R. pseudochrysanthum s.l., the spatial distribution of the cpDNA had a noteworthy pattern showing high genetic differentiation (F_ST = 0.56-0.72) between populations in the Yushan Mountain Range and populations of the other mountain ranges.

Conclusion

Both incomplete lineage sorting and interspecific hybridization/introgression may have contributed to the lack of monophyly among R. hyperythrum, R. formosana and R. pseudochrysanthum s.l. Independent colonizations, plus low capabilities of seed dispersal in current environments, may have resulted in the genetic differentiation between populations of different mountain ranges. At the population level, the populations of Central, and Sheishan Mountains may have undergone postglacial demographic expansion, while populations of the Yushan Mountain Range are likely to have remained stable ever since the colonization. In contrast, the single population of the Alishan Mountain Range with a fixed cpDNA haplotype may have experienced bottleneck/founder's events.

Lineage-specific late pleistocene expansion of an endemic subtropical gossamer-wing damselfly, Euphaea formosa, in Taiwan

BACKGROUND

Pleistocene glacial oscillations have significantly affected the historical population dynamics of temperate taxa. However, the general effects of recent climatic changes on the evolutionary history and genetic structure of extant subtropical species remain poorly understood. In the present study, phylogeographic and historical demographic analyses based on mitochondrial and nuclear DNA sequences were used. The aim was to investigate whether Pleistocene climatic cycles, paleo-drainages or mountain vicariance of Taiwan shaped the evolutionary diversification of a subtropical gossamer-wing damselfly, Euphaea formosa.

RESULTS

E. formosa populations originated in the middle Pleistocene period (0.3 Mya) and consisted of two evolutionarily independent lineages. It is likely that they derived from the Pleistocene paleo-drainages of northern and southern Minjiang, or alternatively by divergence within Taiwan. The ancestral North-central lineage colonized northwestern Taiwan first and maintained a slowly growing population throughout much of the early to middle Pleistocene period. The ancestral widespread lineage reached central-southern Taiwan and experienced a spatial and demographic expansion into eastern Taiwan. This expansion began approximately 30,000 years ago in the Holocene interglacial period. The ancestral southern expansion into eastern Taiwan indicates that the central mountain range (CMR) formed a barrier to east-west expansion. However, E. formosa populations in the three major biogeographic regions (East, South, and North-Central) exhibit no significant genetic partitions, suggesting that river drainages and mountains did not form strong geographical barriers against gene flow among extant populations.

CONCLUSIONS

The present study implies that the antiquity of E. formosa's colonization is associated with its high dispersal ability and larval tolerance to the late Pleistocene dry grasslands. The effect of late Pleistocene climatic changes on the subtropical damselfly's historical demography is lineage-specific, depending predominantly on its colonization history and geography. It is proposed that the Riss and Würm glaciations in the late Pleistocene period had a greater impact on the evolutionary diversification of subtropical insular species than the last glacial maximum (LGM).

A next-generation sequencing method for overcoming the multiple gene copy problem in polyploid phylogenetics, applied to Poa grasses

BACKGROUND

Polyploidy is important from a phylogenetic perspective because of its immense past impact on evolution and its potential future impact on diversification, survival and adaptation, especially in plants. Molecular population genetics studies of polyploid organisms have been difficult because of problems in sequencing multiple-copy nuclear genes using Sanger sequencing. This paper describes a method for sequencing a barcoded mixture of targeted gene regions using next-generation sequencing methods to overcome these problems.

RESULTS

Using 64 3-bp barcodes, we successfully sequenced three chloroplast and two nuclear gene regions (each of which contained two gene copies with up to two alleles per individual) in a total of 60 individuals across 11 species of Australian Poa grasses. This method had high replicability, a low sequencing error rate (after appropriate quality control) and a low rate of missing data. Eighty-eight percent of the 320 gene/individual combinations produced sequence reads, and >80% of individuals produced sufficient reads to detect all four possible nuclear alleles of the homeologous nuclear loci with 95% probability.We applied this method to a group of sympatric Australian alpine Poa species, which we discovered to share an allopolyploid ancestor with a group of American Poa species. All markers revealed extensive allele sharing among the Australian species and so we recommend that the current taxonomy be re-examined. We also detected hypermutation in the trnH-psbA marker, suggesting it should not be used as a land plant barcode region. Some markers indicated differentiation between Tasmanian and mainland samples. Significant positive spatial genetic structure was detected at <100 km with chloroplast but not nuclear markers, which may be a result of restricted seed flow and long-distance pollen flow in this wind-pollinated group.

CONCLUSIONS

Our results demonstrate that 454 sequencing of barcoded amplicon mixtures can be used to reliably sample all alleles of homeologous loci in polyploid species and successfully investigate phylogenetic relationships among species, as well as to investigate phylogeographic hypotheses. This next-generation sequencing method is more affordable than and at least as reliable as bacterial cloning. It could be applied to any experiment involving sequencing of amplicon mixtures.

Estimating species trees using approximate Bayesian computation

Development of methods for estimating species trees from multilocus data is a current challenge in evolutionary biology. We propose a method for estimating the species tree topology and branch lengths using approximate Bayesian computation (ABC). The method takes as data a sample of observed rooted gene tree topologies, and then iterates through the following sequence of steps: First, a randomly selected species tree is used to compute the distribution of rooted gene tree topologies. This distribution is then compared to the observed gene topology frequencies, and if the fit between the observed and the predicted distributions is close enough, the proposed species tree is retained. Repeating this many times leads to a collection of retained species trees that are then used to form the estimate of the overall species tree. We test the performance of the method, which we call ST-ABC, using both simulated and empirical data. The simulation study examines both symmetric and asymmetric species trees over a range of branch lengths and sample sizes. The results from the simulation study show that the model performs very well, giving accurate estimates for both the topology and the branch lengths across the conditions studied, and that a sample size of 25 loci appears to be adequate for the method. Further, we apply the method to two empirical cases: a 4-taxon data set for primates and a 7-taxon data set for yeast. In both cases, we find that estimates obtained with ST-ABC agree with previous studies. The method provides efficient estimation of the species tree, and does not require sequence data, but rather the observed distribution of rooted gene topologies without branch lengths. Therefore, this method is a useful alternative to other currently available methods for species tree estimation.

Joint inference of microsatellite mutation models, population history and genealogies using transdimensional Markov Chain Monte Carlo

We provide a framework for Bayesian coalescent inference from microsatellite data that enables inference of population history parameters averaged over microsatellite mutation models. To achieve this we first implemented a rich family of microsatellite mutation models and related components in the software package BEAST. BEAST is a powerful tool that performs Bayesian MCMC analysis on molecular data to make coalescent and evolutionary inferences. Our implementation permits the application of existing nonparametric methods to microsatellite data. The implemented microsatellite models are based on the replication slippage mechanism and focus on three properties of microsatellite mutation: length dependency of mutation rate, mutational bias toward expansion or contraction, and number of repeat units changed in a single mutation event. We develop a new model that facilitates microsatellite model averaging and Bayesian model selection by transdimensional MCMC. With Bayesian model averaging, the posterior distributions of population history parameters are integrated across a set of microsatellite models and thus account for model uncertainty. Simulated data are used to evaluate our method in terms of accuracy and precision of θ estimation and also identification of the true mutation model. Finally we apply our method to a red colobus monkey data set as an example.

Coalescent histories on phylogenetic networks and detection of hybridization despite incomplete lineage sorting

Analyses of the increasingly available genomic data continue to reveal the extent of hybridization and its role in the evolutionary diversification of various groups of species. We show, through extensive coalescent-based simulations of multilocus data sets on phylogenetic networks, how divergence times before and after hybridization events can result in incomplete lineage sorting with gene tree incongruence signatures identical to those exhibited by hybridization. Evolutionary analysis of such data under the assumption of a species tree model can miss all hybridization events, whereas analysis under the assumption of a species network model would grossly overestimate hybridization events. These issues necessitate a paradigm shift in evolutionary analysis under these scenarios, from a model that assumes a priori a single source of gene tree incongruence to one that integrates multiple sources in a unifying framework. We propose a framework of coalescence within the branches of a phylogenetic network and show how this framework can be used to detect hybridization despite incomplete lineage sorting. We apply the model to simulated data and show that the signature of hybridization can be revealed as long as the interval between the divergence times of the species involved in hybridization is not too small. We reanalyze a data set of 106 loci from 7 in-group Saccharomyces species for which a species tree with no hybridization has been reported in the literature. Our analysis supports the hypothesis that hybridization occurred during the evolution of this group, explaining a large amount of the incongruence in the data. Our findings show that an integrative approach to gene tree incongruence and its reconciliation is needed. Our framework will help in systematically analyzing genomic data for the occurrence of hybridization and elucidating its evolutionary role.

Demographic signatures accompanying the evolution of selfing in Leavenworthia alabamica

The evolution of selfing from outcrossing is a common transition, yet little is known about the mutations and selective factors that promote this shift. In the mustard family, single-locus self-incompatibility (SI) enforces outcrossing. In this study, we test whether mutations causing self-compatibility (SC) are linked to the self-incompatibility locus (S-locus) in Leavenworthia alabamica, a species where two selfing races (a2 and a4) co-occur with outcrossing populations. We also infer the ecological circumstances associated with origins of selfing using molecular sequence data. Genealogical reconstruction of the Lal2 locus, the putative ortholog of the SRK locus, showed that both selfing races are fixed for one of two different S-linked Lal2 sequences, whereas outcrossing populations harbor many S-alleles. Hybrid crosses demonstrated that S-linked mutations cause SC in each selfing race. These results strongly suggest two origins of selfing in this species, a result supported by population admixture analysis of 16 microsatellite loci and by a population tree built from eight nuclear loci. One selfing race (a4) shows signs of a severe population bottleneck, suggesting that reproductive assurance might have caused the evolution of selfing in this case. In contrast, the population size of race a2 cannot be distinguished from that of outcrossing populations after correcting for differences in selfing rates. Coalescent-based analyses suggest a relatively old origin of selfing in the a4 race (∼150 ka ago), whereas selfing evolved recently in the a2 race (∼12-48 ka ago). These results imply that S-locus mutations have triggered two recent shifts to selfing in L. alabamica, but that these transitions are not always associated with a severe population bottleneck, suggesting that factors other than reproductive assurance may play a role in its evolution.

Testing spatiotemporal hypothesis of bacterial evolution using methicillin-resistant Staphylococcus aureus ST239 genome-wide data within a bayesian framework

Staphylococcus aureus is a common cause of infections that has undergone rapid global spread over recent decades. Formal phylogeographic methods have not yet been applied to the molecular epidemiology of bacterial pathogens because the limited genetic diversity of data sets based on individual genes usually results in poor phylogenetic resolution. Here, we investigated a whole-genome single nucleotide polymorphism (SNP) data set of health care-associated Methicillin-resistant S. aureus sequence type 239 (HA-MRSA ST239) strains, which we analyzed using Markov spatial models that incorporate geographical sampling distributions. The reconstructed timescale indicated a temporal origin of this strain shortly after the introduction of Methicillin, followed by global pandemic spread. The estimate of the temporal origin was robust to the molecular clock, coalescent prior, full/intergenic/synonymous SNP inclusion, and correction for excluded invariant site patterns. Finally, phylogeographic analyses statistically supported the role of human movement in the global dissemination of HA-MRSA ST239, although it was unable to conclusively resolve the location of the root. This study demonstrates that bacterial genomes can indeed contain sufficient evolutionary information to elucidate the temporal and spatial dynamics of transmission. Future applications of this approach to other bacterial strains may provide valuable epidemiological insights that may justify the cost of genome-wide typing.

Three roads diverged? Routes to phylogeographic inference

Phylogeographic methods facilitate inference of the geographical history of genetic lineages. Recent examples explore human migration and the origins of viral pandemics. There is longstanding disagreement over the use and validity of certain phylogeographic inference methodologies. In this paper, we highlight three distinct frameworks for phylogeographic inference to give a taste of this disagreement. Each of the three approaches presents a different viewpoint on phylogeography, most fundamentally on how we view the relationship between the inferred history of a sample and the history of the population the sample is embedded in. Satisfactory resolution of this relationship between history of the tree and history of the population remains a challenge for all but the most trivial models of phylogeographic processes. Intriguingly, we believe that some recent methods that entirely avoid inference about the history of the population will eventually help to reach a resolution.

A young clade repeating an old pattern: diversity in Nothonotus darters (Teleostei: Percidae) endemic to the Cumberland River

Hypotheses of diversification in eastern North American freshwater fishes have focused primarily on allopatric distributions of species between disjunct highland areas and major river systems. However, these hypotheses do not fully explain the rich diversity of species within highland regions and river systems. Relatively old diversification events at small geographic scales have been observed in the Barcheek Darter subclade that occurs in the Cumberland River drainage (CRD) in Kentucky and Tennessee, United States of America, but it is unknown if this pattern is consistent in other darter subclades. We explored phylogeographic diversity in two species of Nothonotus darters, N. microlepidus and N. sanguifluus, endemic to the CRD to compare phylogenetic patterns between Barcheek Darters and species of Nothonotus. We collected sequence data for a mitochondrial gene (cytb) and three nuclear genes (MLL, S7 and RAG1) from 19 N. microlepidus and 35 N. sanguifluus specimens. Gene trees were estimated using maximum likelihood and Bayesian methods, and a 'species tree' was inferred using a Bayesian method. These trees indicate that species diversity in Nothonotus is underestimated. Five distinct lineages were evident, despite retained ancestral polymorphism and unsampled extirpated populations. Comparison of chronograms for Barcheek Darters and Nothonotus revealed that microendemism resulting from species diversification at small geographic scales in the CRD is a consistent pattern in both old and young darter subclades. Our analyses reveal that geographic isolating mechanisms that result in similar phylogeographic patterns in the CRD are persistent through long expanses of evolutionary time.

BUCKy: gene tree/species tree reconciliation with Bayesian concordance analysis

MOTIVATION

BUCKy is a C++ program that implements Bayesian concordance analysis. The method uses a non-parametric clustering of genes with compatible trees, and reconstructs the primary concordance tree from clades supported by the largest proportions of genes. A population tree with branch lengths in coalescent units is estimated from quartet concordance factors.

AVAILABILITY

BUCKy is open source and distributed under the GNU general public license at www.stat.wisc.edu/∼ane/bucky/.

Recent trends in molecular phylogenetic analysis: where to next?

The acquisition of large multilocus sequence data is providing researchers with an unprecedented amount of information to resolve difficult phylogenetic problems. With these large quantities of data comes the increasing challenge regarding the best methods of analysis. We review the current trends in molecular phylogenetic analysis, focusing specifically on the topics of multiple sequence alignment and methods of tree reconstruction. We suggest that traditional methods are inadequate for these highly heterogeneous data sets and that researchers employ newer more sophisticated search algorithms in their analyses. If we are to best extract the information present in these data sets, a sound understanding of basic phylogenetic principles combined with modern methodological techniques are necessary.

Identifying the rooted species tree from the distribution of unrooted gene trees under the coalescent

Gene trees are evolutionary trees representing the ancestry of genes sampled from multiple populations. Species trees represent populations of individuals-each with many genes-splitting into new populations or species. The coalescent process, which models ancestry of gene copies within populations, is often used to model the probability distribution of gene trees given a fixed species tree. This multispecies coalescent model provides a framework for phylogeneticists to infer species trees from gene trees using maximum likelihood or Bayesian approaches. Because the coalescent models a branching process over time, all trees are typically assumed to be rooted in this setting. Often, however, gene trees inferred by traditional phylogenetic methods are unrooted. We investigate probabilities of unrooted gene trees under the multispecies coalescent model. We show that when there are four species with one gene sampled per species, the distribution of unrooted gene tree topologies identifies the unrooted species tree topology and some, but not all, information in the species tree edges (branch lengths). The location of the root on the species tree is not identifiable in this situation. However, for 5 or more species with one gene sampled per species, we show that the distribution of unrooted gene tree topologies identifies the rooted species tree topology and all its internal branch lengths. The length of any pendant branch leading to a leaf of the species tree is also identifiable for any species from which more than one gene is sampled.

Phylogenetic and coalescent analysis of three loci suggest that the Water Rail is divisible into two species, Rallus aquaticus and R. indicus

BACKGROUND

Water Rails (Rallus aquaticus) inhabit fragmented freshwater wetlands across their Palearctic distribution. Disjunct populations are now thought to be morphologically similar over their vast geographic range, though four subspecies had been recognized previously. The fossil record suggests that Water Rails (R. aquaticus) were already spread across the Palearctic by the Pleistocene approximately 2 million years ago, and the oldest fossil remains thought to be closely related to the common ancestor of water rails date from the Pliocene.

RESULTS

To investigate population structure in Water Rails at the genetic level we sequenced three independent loci: 686 base pairs (bp) of the mitochondrial DNA COI barcode; 618 bp of the intron ADH5; and 746 bp of the exon PTPN12. Phylogeographic analysis revealed that Water Rails breeding in eastern Asia (R. a. indicus, also known as the Brown-cheeked Rail) are strongly differentiated from the Water Rails in Western and Middle Asia and Europe (R. a. aquaticus and R. a. korejewi). The Kimura 3-parameter plus Gamma COI genetic distance between these two geographic groups was > 3%, and they differed by 18 diagnostic substitutions commensurate with differences between recently diverged sister species of birds. In spite of the low number of variable sites, the two nuclear loci supported this split. We estimated the split of the Brown-cheeked Rail and the Water Rail to have occurred approximately 534,000 years ago (95% CI 275,000-990,000 years ago). Fragmentation of the widespread ancestral population and eventual speciation of water rails is likely attributable to vicariance by a barrier formed by glacial cycles, continuous uplift of the Tibetan Plateau and increased sedimentation in deserts in southern Asia that originated in the Miocene.

CONCLUSIONS

Water Rails from East Asia were genetically differentiated from the ones breeding in Europe and Western to Middle Asia. Most of the genetic signal was from mitochondrial COI, and was corroborated by polymorphic sites in the two nuclear loci we employed. The split between these two lineages was estimated to occur in the Middle Pleistocene, when populations were isolated in disjunct wetlands with little or no gene flow. Independent evidence from differences in morphology and vocalizations in concert with genetic differentiation and a long history of isolation support recognition of the Brown-cheeked Rail breeding in East Asia as a separate species, R. indicus. The use of several independent loci is invaluable in inferring species trees from gene trees and in recognizing species limits.

Repeated genomic transfers from echovirus 30 to echovirus 6 lineages indicate co-divergence between co-circulating populations of the two human enterovirus serotypes

Human echovirus types 6 (E-6) and 30 (E-30) cause seasonal epidemics of aseptic meningitis. These two enteroviruses are frequently observed in co-circulation, an epidemiological pattern that is prerequisite for the occurrence of dual infections, which can lead to recombination between co-infecting virus strains. Viral sequences were determined at loci 1D (VP1 capsid protein) and 3CD (non structural proteins) in 49 E-6 strains recovered in a single geographical region in France from 1999 to 2007, during the epidemiological survey of enterovirus infections. They were compared with previously recorded sequences of E-30 strains to investigate their evolutionary histories and possible recombination patterns. Phylogenetic analyses identified two distinct E-6 populations and different subpopulations. Assuming a relaxed molecular clock model and a Bayesian skyline demographic model in coalescent analyses with the BEAST program, the substitution rate in E-6 was estimated at 8.597×10(-3) and 6.252×10(-3) substitution/site/year for loci 1D and 3CD respectively. Consistent estimates of divergence times (t(MRCA)) were obtained for loci 1D and 3CD indicating that two distinct E-6 populations originated in 1997 and 1999. Incongruent phylogenetic patterns inferred for the two loci were indicative of recombination events between the two populations. Phylogenies including the E-30 3CD sequences showed close genetic relationships between E-6 and discrete E-30 subpopulations. Recombination breakpoints were located with statistical significance in E-6 and E-30 genomes. Estimates of t(MRCA) of phylogenetic recombinant clades indicated directional genetic transfers from E-30 to E-6 populations and their co-divergence over the time period studied.

Bayesian species delimitation in West African forest geckos (Hemidactylus fasciatus)

Genealogical data are an important source of evidence for delimiting species, yet few statistical methods are available for calculating the probabilities associated with different species delimitations. Bayesian species delimitation uses reversible-jump Markov chain Monte Carlo (rjMCMC) in conjunction with a user-specified guide tree to estimate the posterior distribution for species delimitation models containing different numbers of species. We apply Bayesian species delimitation to investigate the speciation history of forest geckos (Hemidactylus fasciatus) from tropical West Africa using five nuclear loci (and mtDNA) for 51 specimens representing 10 populations. We find that species diversity in H. fasciatus is currently underestimated, and describe three new species to reflect the most conservative estimate for the number of species in this complex. We examine the impact of the guide tree, and the prior distributions on ancestral population sizes () and root age (τ(0)), on the posterior probabilities for species delimitation. Mis-specification of the guide tree or the prior distribution for can result in strong support for models containing more species. We describe a new statistic for summarizing the posterior distribution of species delimitation models, called speciation probabilities, which summarize the posterior support for each speciation event on the starting guide tree.

Bayesian species delimitation using multilocus sequence data

In the absence of recent admixture between species, bipartitions of individuals in gene trees that are shared across loci can potentially be used to infer the presence of two or more species. This approach to species delimitation via molecular sequence data has been constrained by the fact that genealogies for individual loci are often poorly resolved and that ancestral lineage sorting, hybridization, and other population genetic processes can lead to discordant gene trees. Here we use a Bayesian modeling approach to generate the posterior probabilities of species assignments taking account of uncertainties due to unknown gene trees and the ancestral coalescent process. For tractability, we rely on a user-specified guide tree to avoid integrating over all possible species delimitations. The statistical performance of the method is examined using simulations, and the method is illustrated by analyzing sequence data from rotifers, fence lizards, and human populations.