Ancestral Effective Population Size Research Articles

In silico phylogenomics using complete genomes: a case study on the evolution of hominoids.

The increasing availability of complete genome data is facilitating the acquisition of phylogenomic data sets, but the process of obtaining orthologous sequences from other genomes and assembling multiple sequence alignments remains piecemeal and arduous. We designed software that performs these tasks and outputs anonymous loci (AL) or anchored enrichment/ultraconserved element loci (AE/UCE) data sets in ready-to-analyze formats. We demonstrate our program by applying it to the hominoids. Starting with human, chimpanzee, gorilla, and orangutan genomes, our software generated an exhaustive data set of 292 ALs (∼1 kb each) in ∼3 h. Not only did analyses of our AL data set validate the program by yielding a portrait of hominoid evolution in agreement with previous studies, but the accuracy and precision of our estimated ancestral effective population sizes and speciation times represent improvements. We also used our program with a published set of 512 vertebrate-wide AE “probe” sequences to generate data sets consisting of 171 and 242 independent loci (∼1 kb each) in 11 and 13 min, respectively. The former data set consisted of flanking sequences 500 bp from adjacent AEs, while the latter contained sequences bordering AEs. Although our AE data sets produced the expected hominoid species tree, coalescent-based estimates of ancestral population sizes and speciation times based on these data were considerably lower than estimates from our AL data set and previous studies. Accordingly, we suggest that loci subjected to direct or indirect selection may not be appropriate for coalescent-based methods. Complete in silico approaches, combined with the burgeoning genome databases, will accelerate the pace of phylogenomics.

Estimation of linkage disequilibrium levels and haplotype block structure in Chinese Simmental and Wagyu beef cattle using high-density genotypes

With the advance of high-throughput genotyping technologies, it is possible to estimate the pattern of linkage disequilibrium (LD) and haplotype block structure at genome wide level in farm animals. Understanding LD characteristics and haplotype blocks can provide useful information to guide genome wide association study (GWAS) and genomic selection (GS). In this study, we investigated the LD pattern and haplotype block structure in two beef cattle populations (Chinese Simmental and Wagyu) using the BovineHD BeadChip. We found the average LD measured by r2 in Chinese Wagyu was higher than that in Simmental. We observed LD decayed more slowly on the X chromosome than autosomes for both populations. The persistence of LD phase varied from 0.96 at a distance of <2.5kb to 0.51 at a distance from 400kb to 500kb. The estimated ancestral effective population sizes (Ne) in Chinese Simmental and Wagyu 5 years ago were approximately 74 and 23, respectively. Moreover, we detected 50,475 (991.1Mb, 69.3%) haplotype blocks shared by two populations, 18,440 (238.9Mb, 16.7%) unique blocks in Chinese Simmental and 15,417 (199.8Mb, 14%) in Wagyu. Our study revealed that the r2 dropped below 0.2 at distances of 34kb and 40kb in Chinese Simmental and Wagyu, which indicated that the implementation of GS for both populations require at least 77,941 and 66,250 markers, respectively. The pattern of LD mirrored the slightly different selection histories for Chinese Simmental and Wagyu, and our result suggested the shared haplotype blocks may offer valuable insights for the implementation of genome wide association and genomic selection studies based on multi-population.

Distinguishing noise from signal in patterns of genomic divergence in a highly polymorphic avian radiation.

Recently diverged taxa provide the opportunity to search for the genetic basis of the phenotypes that distinguish them. Genomic scans aim to identify loci that are diverged with respect to an otherwise weakly differentiated genetic background. These loci are candidates for being past targets of selection because they behave differently from the rest of the genome that has either not yet differentiated or that may cross species barriers through introgressive hybridization. Here we use a reduced-representation genomic approach to explore divergence among six species of southern capuchino seedeaters, a group of recently radiated sympatric passerine birds in the genus Sporophila. For the first time in these taxa, we discovered a small proportion of markers that appeared differentiated among species. However, when assessing the significance of these signatures of divergence, we found that similar patterns can also be recovered from random grouping of individuals representing different species. A detailed demographic inference indicates that genetic differences among Sporophila species could be the consequence of neutral processes, which include a very large ancestral effective population size that accentuates the effects of incomplete lineage sorting. As these neutral phenomena can generate genomic scan patterns that mimic those of markers involved in speciation and phenotypic differentiation, they highlight the need for caution when ascertaining and interpreting differentiated markers between species, especially when large numbers of markers are surveyed. Our study provides new insights into the demography of the southern capuchino radiation and proposes controls to distinguish signal from noise in similar genomic scans.

Similarity thresholds used in DNA sequence assembly from short reads can reduce the comparability of population histories across species.

Comparing inferences among datasets generated using short read sequencing may provide insight into the concerted impacts of divergence, gene flow and selection across organisms, but comparisons are complicated by biases introduced during dataset assembly. Sequence similarity thresholds allow the de novo assembly of short reads into clusters of alleles representing different loci, but the resulting datasets are sensitive to both the similarity threshold used and to the variation naturally present in the organism under study. Thresholds that require high sequence similarity among reads for assembly (stringent thresholds) as well as highly variable species may result in datasets in which divergent alleles are lost or divided into separate loci (‘over-splitting’), whereas liberal thresholds increase the risk of paralogous loci being combined into a single locus (‘under-splitting’). Comparisons among datasets or species are therefore potentially biased if different similarity thresholds are applied or if the species differ in levels of within-lineage genetic variation. We examine the impact of a range of similarity thresholds on assembly of empirical short read datasets from populations of four different non-model bird lineages (species or species pairs) with different levels of genetic divergence. We find that, in all species, stringent similarity thresholds result in fewer alleles per locus than more liberal thresholds, which appears to be the result of high levels of over-splitting. The frequency of putative under-splitting, conversely, is low at all thresholds. Inferred genetic distances between individuals, gene tree depths, and estimates of the ancestral mutation-scaled effective population size (θ) differ depending upon the similarity threshold applied. Relative differences in inferences across species differ even when the same threshold is applied, but may be dramatically different when datasets assembled under different thresholds are compared. These differences not only complicate comparisons across species, but also preclude the application of standard mutation rates for parameter calibration. We suggest some best practices for assembling short read data to maximize comparability, such as using more liberal thresholds and examining the impact of different thresholds on each dataset.

Human paternal and maternal demographic histories: insights from high-resolution Y chromosome and mtDNA sequences.

BackgroundComparisons of maternally-inherited mitochondrial DNA (mtDNA) and paternally-inherited non-recombining Y chromosome (NRY) variation have provided important insights into the impact of sex-biased processes (such as migration, residence pattern, and so on) on human genetic variation. However, such comparisons have been limited by the different molecular methods typically used to assay mtDNA and NRY variation (for example, sequencing hypervariable segments of the control region for mtDNA vs. genotyping SNPs and/or STR loci for the NRY). Here, we report a simple capture array method to enrich Illumina sequencing libraries for approximately 500 kb of NRY sequence, which we use to generate NRY sequences from 623 males from 51 populations in the CEPH Human Genome Diversity Panel (HGDP). We also obtained complete mtDNA genome sequences from the same individuals, allowing us to compare maternal and paternal histories free of any ascertainment bias.ResultsWe identified 2,228 SNPs in the NRY sequences and 2,163 SNPs in the mtDNA sequences. Our results confirm the controversial assertion that genetic differences between human populations on a global scale are bigger for the NRY than for mtDNA, although the differences are not as large as previously suggested. More importantly, we find substantial regional variation in patterns of mtDNA versus NRY variation. Model-based simulations indicate very small ancestral effective population sizes (<100) for the out-of-Africa migration as well as for many human populations. We also find that the ratio of female effective population size to male effective population size (Nf/Nm) has been greater than one throughout the history of modern humans, and has recently increased due to faster growth in Nf than Nm.ConclusionsThe NRY and mtDNA sequences provide new insights into the paternal and maternal histories of human populations, and the methods we introduce here should be widely applicable for further such studies.

Sex-biased dispersal produces high error rates in mitochondrial distance-based and tree-based species delimitation

Abstract Species delimitation using mitochondrial sequences aims to identify species as morphological expertise and biodiversity both decline. Species delimitation in animals relies completely or in part on 2 criteria: genetic distance and reciprocal monophyly. Using coalescent simulations of populations experiencing continuous and interrupted gene flow we show that these commonly applied criteria incur both high false-positive and high false-negative error rates in species delimitation when dispersal is sex biased, as it is in most mammals. The combination of distance- and topology-based criteria will incur false-positive error rates well above 5% when ancestral effective population sizes are large and when population structure has been in place for 10,000 years or more. This effect persists even with a subtle bias in dispersal of 1 : 4 female : males dispersing in each generation. High false-positive error rates in genetic distances can be overcome by correcting for within-population sequence divergence...

Estimation of the ancestral effective population sizes of African great apes under different selection regimes.

Reliable estimates of ancestral effective population sizes are necessary to unveil the population-level phenomena that shaped the phylogeny and molecular evolution of the African great apes. Although several methods have previously been applied to infer ancestral effective population sizes, an analysis of the influence of the selective regime on the estimates of ancestral demography has not been thoroughly conducted. In this study, three independent data sets under different selective regimes were used were composed to tackle this issue. The results showed that selection had a significant impact on the estimates of ancestral effective population sizes of the African great apes. The inference of the ancestral demography of African great apes was affected by the selection regime. The effects, however, were not homogeneous along the ancestral populations of great apes. The effective population size of the ancestor of humans and chimpanzees was more impacted by the selection regime when compared to the same parameter in the ancestor of humans, chimpanzees and gorillas. Because the selection regime influenced the estimates of ancestral effective population size, it is reasonable to assume that a portion of the discrepancy found in previous studies that inferred the ancestral effective population size may be attributable to the differential action of selection on the genes sampled.

Detecting past changes of effective population size.

Understanding and predicting population abundance is a major challenge confronting scientists. Several genetic models have been developed using microsatellite markers to estimate the present and ancestral effective population sizes. However, to get an overview on the evolution of population requires that past fluctuation of population size be traceable. To address the question, we developed a new model estimating the past changes of effective population size from microsatellite by resolving coalescence theory and using approximate likelihoods in a Monte Carlo Markov Chain approach. The efficiency of the model and its sensitivity to gene flow and to assumptions on the mutational process were checked using simulated data and analysis. The model was found especially useful to provide evidence of transient changes of population size in the past. The times at which some past demographic events cannot be detected because they are too ancient and the risk that gene flow may suggest the false detection of a bottleneck are discussed considering the distribution of coalescence times. The method was applied on real data sets from several Atlantic salmon populations. The method called VarEff (Variation of Effective size) was implemented in the R package VarEff and is made available at https://qgsp.jouy.inra.fr and at http://cran.r-project.org/web/packages/VarEff.

The limiting distribution of the effective population size of the ancestor of humans and chimpanzees

The effective population size is a fundamental parameter for the understanding of microevolutionary process. Indeed, the consideration of population-level phenomena within phylogenies provides insight into the influence of the past evolutionary demography on the genetic diversity of living species. Although the effective population size of the last common ancestor of humans and chimpanzees has been extensively investigated by molecular evolutionists, variance in the estimates of this parameter among studies is large. However, with the availability of genome sequences, the estimation of evolutionary parameters may be conducted with minimum stochastic error, and the limiting distribution of the estimates may be obtained. This statistical property was utilized in the present study and coupled with analytical derivations from the coalescent theory to examine the limiting distribution of the ancestral effective population size of Homo–Pan. The mean ancestral effective population size of Homo–Pan was inferred at approximately 47,500, and the results showed that the uncertainty of the estimates was large, even under the limiting distribution. Further reductions of the estimates are feasible only if additional calibration information from the fossil record is provided and if a probabilistic model of ancestral generation time is envisioned.

Radically different phylogeographies and patterns of genetic variation in two European brown frogs, genus Rana

We reconstruct range–wide phylogeographies of two widespread and largely co-occurring Western Palearctic frogs, Rana temporaria and R. dalmatina. Based on tissue or saliva samples of over 1000 individuals, we compare a variety of genetic marker systems, including mitochondrial DNA, single-copy protein-coding nuclear genes, microsatellite loci, and single nucleotide polymorphisms (SNPs) of transcriptomes of both species. The two focal species differ radically in their phylogeographic structure, with R. temporaria being strongly variable among and within populations, and R. dalmatina homogeneous across Europe with a single strongly differentiated population in southern Italy. These differences were observed across the various markers studied, including microsatellites and SNP density, but especially in protein-coding nuclear genes where R. dalmatina had extremely low heterozygosity values across its range, including potential refugial areas. On the contrary, R. temporaria had comparably high range-wide values, including many areas of probable postglacial colonization. A phylogeny of R. temporaria based on various concatenated mtDNA genes revealed that two haplotype clades endemic to Iberia form a paraphyletic group at the base of the cladogram, and all other haplotypes form a monophyletic group, in agreement with an Iberian origin of the species. Demographic analysis suggests that R. temporaria and R. dalmatina have genealogies of roughly the same time to coalescence (TMRCA ∼3.5mya for both species), but R. temporaria might have been characterized by larger ancestral and current effective population sizes than R. dalmatina. The high genetic variation in R. temporaria can therefore be explained by its early range expansion out of Iberia, with subsequent cycles of differentiation in cryptic glacial refugial areas followed by admixture, while the range expansion of R. dalmatina into central Europe is a probably more recent event.

Multilocus estimation of divergence times and ancestral effective population sizes of Oryza species and implications for the rapid diversification of the genus

· Despite substantial investigations into Oryza phylogeny and evolution, reliable estimates of the divergence times and ancestral effective population sizes of major lineages in Oryza are challenging. · We sampled sequences of 106 single-copy nuclear genes from all six diploid genomes of Oryza to investigate the divergence times through extensive relaxed molecular clock analyses and estimated the ancestral effective population sizes using maximum likelihood and Bayesian methods. · We estimated that Oryza originated in the middle Miocene (c. 13-15 million years ago; Ma) and obtained an explicit time frame for two rapid diversifications in this genus. The first diversification involving the extant F-/G-genomes and possibly the extinct H-/J-/K-genomes occurred in the middle Miocene immediately after (within < 1 Myr) the origin of Oryza. The second giving rise to the A-/B-/C-genomes happened c. 5-6 Ma. We found that ancestral effective population sizes were much larger than those of extant species in Oryza. · We suggest that the climate fluctuations during the period from the middle Miocene to Pliocene may have contributed to the two rapid diversifications of Oryza species. Such information helps better understand the evolutionary history of Oryza and provides further insights into the pattern and mechanism of diversification in plants in general.

The Pipid Root

The estimation of phylogenetic relationships is an essential component of understanding evolution. Accurate phylogenetic estimation is difficult, however, when internodes are short and old, when genealogical discordance is common due to large ancestral effective population sizes or ancestral population structure, and when homoplasy is prevalent. Inference of divergence times is also hampered by unknown and uneven rates of evolution, the incomplete fossil record, uncertainty in relationships between fossil and extant lineages, and uncertainty in the age of fossils. Ideally, these challenges can be overcome by developing large "phylogenomic" data sets and by analyzing them with methods that accommodate features of the evolutionary process, such as genealogical discordance, recurrent substitution, recombination, ancestral population structure, gene flow after speciation among sampled and unsampled taxa, and variation in evolutionary rates. In some phylogenetic problems, it is possible to use information that is independent of fossils, such as the geological record, to identify putative triggers for diversification whose associated estimated divergence times can then be compared a posteriori with estimated relationships and ages of fossils. The history of diversification of pipid frog genera Pipa, Hymenochirus, Silurana, and Xenopus, for instance, is characterized by many of these evolutionary and analytical challenges. These frogs diversified dozens of millions of years ago, they have a relatively rich fossil record, their distributions span continental plates with a well characterized geological record of ancient connectivity, and there is considerable disagreement across studies in estimated evolutionary relationships. We used high throughput sequencing and public databases to generate a large phylogenomic data set with which we estimated evolutionary relationships using multilocus coalescence methods. We collected sequence data from Pipa, Hymenochirus, Silurana, and Xenopus and the outgroup taxon Rhinophrynus dorsalis from coding sequence of 113 autosomal regions, averaging ∼300 bp in length (range: 102-1695 bp) and also a portion of the mitochondrial genome. Analysis of these data using multiple approaches recovers strong support for the ((Xenopus, Silurana)(Pipa, Hymenochirus)) topology, and geologically calibrated divergence time estimates that are consistent with estimated ages and phylogenetic affinities of many fossils. These results provide new insights into the biogeography and chronology of pipid diversification during the breakup of Gondwanaland and illustrate how phylogenomic data may be necessary to tackle tough problems in molecular systematics. [Coalescence; gene tree; high-throughout sequencing; lineage sorting; pipid; species tree; Xenopus.].

The importance of environmental heterogeneity and spatial distance in generating phylogeographic structure in edaphic specialist and generalist tree species of Protium (Burseraceae) across the Amazon Basin

AbstractAim Edaphic heterogeneity may be an important driver of population differentiation in the Amazon but remains to be investigated in trees. We compared the phylogeographic structure across the geographic distribution of two Protium (Burseraceae) species with different degrees of edaphic specialization: Protium alvarezianum, an edaphic specialist of white‐sand habitat islands; and Protium subserratum, an edaphic generalist found in white sand as well as in more widespread soil types. We predicted that in the edaphic specialist, geographic distance would structure populations more strongly than in the edaphic generalist, and that soil type would not structure populations in the edaphic generalist unless habitat acts as a barrier promoting population differentiation.Location Tropical rain forests of the Peruvian and Brazilian Amazon, Guyana and French Guiana.Methods We sequenced 1209–1211 bp of non‐coding nuclear ribosomal DNA (internal transcribed spacer and external transcribed spacer) and a neutral low‐copy nuclear gene (phytochrome C) from P. subserratum (n = 65, 10 populations) and P. alvarezianum (n = 19, three populations). We conducted a Bayesian phylogenetic analysis, constructed maximum parsimony haplotype networks and assessed population differentiation among groups (soil type or geographic locality) using analysis of molecular variance and spatial analysis of molecular variance.Results The edaphic specialist exhibited considerable genetic differentiation among geographically distant populations. The edaphic generalist showed significant genetic differentiation between the Guianan and Amazon Basin populations. Within Peru, soil type and not geographic distance explained most of the variation among populations. Non‐white‐sand populations in Peru exhibited lower haplotype/nucleotide diversity than white‐sand populations, were each other’s close relatives, and formed an unresolved clade derived from within the white‐sand populations.Main conclusions Geographic distance is a stronger driver of population differentiation in the edaphic specialist than in the generalist. However, this difference did not appear to be related to edaphic generalism per se as adjacent populations from both soil types in the edaphic generalist did not share many haplotypes. Populations of the edaphic generalist in white‐sand habitats exhibited high haplotype diversity and shared haplotypes with distant white‐sand habitat islands, indicating that they have either efficient long‐distance dispersal and/or larger ancestral effective population sizes and thus retain ancestral polymorphisms. These results highlight the importance of edaphic heterogeneity in promoting population differentiation in tropical trees.

Nuclear and plastid haplotypes suggest rapid diploid and polyploid speciation in the N Hemisphere Achillea millefolium complex (Asteraceae)

BackgroundSpecies complexes or aggregates consist of a set of closely related species often of different ploidy levels, whose relationships are difficult to reconstruct. The N Hemisphere Achillea millefolium aggregate exhibits complex morphological and genetic variation and a broad ecological amplitude. To understand its evolutionary history, we study sequence variation at two nuclear genes and three plastid loci across the natural distribution of this species complex and compare the patterns of such variations to the species tree inferred earlier from AFLP data.ResultsAmong the diploid species of A. millefolium agg., gene trees of the two nuclear loci, ncpGS and SBP, and the combined plastid fragments are incongruent with each other and with the AFLP tree likely due to incomplete lineage sorting or secondary introgression. In spite of the large distributional range, no isolation by distance is found. Furthermore, there is evidence for intragenic recombination in the ncpGS gene. An analysis using a probabilistic model for population demographic history indicates large ancestral effective population sizes and short intervals between speciation events. Such a scenario explains the incongruence of the gene trees and species tree we observe. The relationships are particularly complex in the polyploid members of A. millefolium agg.ConclusionsThe present study indicates that the diploid members of A. millefolium agg. share a large part of their molecular genetic variation. The findings of little lineage sorting and lack of isolation by distance is likely due to short intervals between speciation events and close proximity of ancestral populations. While previous AFLP data provide species trees congruent with earlier morphological classification and phylogeographic considerations, the present sequence data are not suited to recover the relationships of diploid species in A. millefolium agg. For the polyploid taxa many hybrid links and introgression from the diploids are suggested.

Comparative multi-locus phylogeography confirms multiple vicariance events in co-distributed rainforest frogs

Though Pleistocene refugia are frequently cited as drivers of species diversification, comparisons of molecular divergence among sister species typically indicate a continuum of divergence times from the Late Miocene, rather than a clear pulse of speciation events at the Last Glacial Maximum. Community-scale inference methods that explicitly test for multiple vicariance events, and account for differences in ancestral effective population size and gene flow, are well suited for detecting heterogeneity of species' responses to past climate fluctuations. We apply this approach to multi-locus sequence data from five co-distributed frog species endemic to the Wet Tropics rainforests of northeast Australia. Our results demonstrate at least two episodes of vicariance owing to climate-driven forest contractions: one in the Early Pleistocene and the other considerably older. Understanding how repeated cycles of rainforest contraction and expansion differentially affected lineage divergence among co-distributed species provides a framework for identifying evolutionary processes that underlie population divergence and speciation.

Genetic diversity, population genetic structure and demographic history of Przewalski’s gazelle (Procapra przewalskii): implications for conservation

The Przewalski’s gazelle (Procapra przewalskii) is one of the most endangered antelope species in the world. It is endemic to China and is a flagship species in the eastern part of the Qinghai–Tibet plateau. To establish effective conservation measures on this species, genetic information such as genetic structure is needed. However, there has not been a comprehensive genetic assessment on this gazelle using nuclear DNA markers yet. Here, we employed 13 microsatellite loci to investigate genetic diversity, population genetic structure and demographic history of Przewalski’s gazelle using noninvasive samples of 169 wild gazelles collected from nine populations. A total of 76 alleles were detected from the entire samples, mean allele number was 5.85, and overall HO and HE were 0.525 and 0.552, respectively. Structure and GENELAND analyses found six genetic groups in the nine populations. Between the inferred genetic groups, significant genetic differentiation and low migration rates were detected. Demographic analyses indicated that Przewalski’s gazelle experienced genetic bottleneck and severe population decline, with the ancestral effective population size reducing to less than one percent. Based on the results of this study, we provide several conservation recommendations for Przewalski’s gazelle, such as six management units, periodic monitoring and special conservation consideration on the Qiejitan population.

Population growth of Mexican free-tailed bats (Tadarida brasiliensis mexicana) predates human agricultural activity

BackgroundHuman activities, such as agriculture, hunting, and habitat modification, exert a significant effect on native species. Although many species have suffered population declines, increased population fragmentation, or even extinction in connection with these human impacts, others seem to have benefitted from human modification of their habitat. Here we examine whether population growth in an insectivorous bat (Tadarida brasiliensis mexicana) can be attributed to the widespread expansion of agriculture in North America following European settlement. Colonies of T. b. mexicana are extremely large (~106 individuals) and, in the modern era, major agricultural insect pests form an important component of their food resource. It is thus hypothesized that the growth of these insectivorous bat populations was coupled to the expansion of agricultural land use in North America over the last few centuries.ResultsWe sequenced one haploid and one autosomal locus to determine the rate and time of onset of population growth in T. b. mexicana. Using an approximate Maximum Likelihood method, we have determined that T. b. mexicana populations began to grow ~220 kya from a relatively small ancestral effective population size before reaching the large effective population size observed today.ConclusionsOur analyses reject the hypothesis that T. b. mexicana populations grew in connection with the expansion of human agriculture in North America, and instead suggest that this growth commenced long before the arrival of humans. As T. brasiliensis is a subtropical species, we hypothesize that the observed signals of population growth may instead reflect range expansions of ancestral bat populations from southern glacial refugia during the tail end of the Pleistocene.

Estimating Divergence Time and Ancestral Effective Population Size of Bornean and Sumatran Orangutan Subspecies Using a Coalescent Hidden Markov Model

Due to genetic variation in the ancestor of two populations or two species, the divergence time for DNA sequences from two populations is variable along the genome. Within genomic segments all bases will share the same divergence—because they share a most recent common ancestor—when no recombination event has occurred to split them apart. The size of these segments of constant divergence depends on the recombination rate, but also on the speciation time, the effective population size of the ancestral population, as well as demographic effects and selection. Thus, inference of these parameters may be possible if we can decode the divergence times along a genomic alignment. Here, we present a new hidden Markov model that infers the changing divergence (coalescence) times along the genome alignment using a coalescent framework, in order to estimate the speciation time, the recombination rate, and the ancestral effective population size. The model is efficient enough to allow inference on whole-genome data sets. We first investigate the power and consistency of the model with coalescent simulations and then apply it to the whole-genome sequences of the two orangutan sub-species, Bornean (P. p. pygmaeus) and Sumatran (P. p. abelii) orangutans from the Orangutan Genome Project. We estimate the speciation time between the two sub-species to be thousand years ago and the effective population size of the ancestral orangutan species to be , consistent with recent results based on smaller data sets. We also report a negative correlation between chromosome size and ancestral effective population size, which we interpret as a signature of recombination increasing the efficacy of selection.

Deep Divergences of Human Gene Trees and Models of Human Origins

Two competing hypotheses are at the forefront of the debate on modern human origins. In the first scenario, known as the recent Out-of-Africa hypothesis, modern humans arose in Africa about 100,000-200,000 years ago and spread throughout the world by replacing the local archaic human populations. By contrast, the second hypothesis posits substantial gene flow between archaic and emerging modern humans. In the last two decades, the young time estimates--between 100,000 and 200,000 years--of the most recent common ancestors for the mitochondrion and the Y chromosome provided evidence in favor of a recent African origin of modern humans. However, the presence of very old lineages for autosomal and X-linked genes has often been claimed to be incompatible with a simple, single origin of modern humans. Through the analysis of a public DNA sequence database, we find, similar to previous estimates, that the common ancestors of autosomal and X-linked genes are indeed very old, living, on average, respectively, 1,500,000 and 1,000,000 years ago. However, contrary to previous conclusions, we find that these deep gene genealogies are consistent with the Out-of-Africa scenario provided that the ancestral effective population size was approximately 14,000 individuals. We show that an ancient bottleneck in the Middle Pleistocene, possibly arising from an ancestral structured population, can reconcile the contradictory findings from the mitochondrion on the one hand, with the autosomes and the X chromosome on the other hand.

Diversification at high latitudes: speciation of buntings in the genusPlectrophenaxinferred from mitochondrial and nuclear markers

High-latitude diversification is a process characterized by speciation and extinction due to climatically driven vicariance and dispersal events. McKay's buntings (Plectrophenax hyperboreus) are high-latitude island endemic songbirds, and their global range is restricted to Beringia. Snow buntings (P. nivalis), their closest relatives, are distributed throughout the Holarctic, breeding in available habitat surrounding the island range of McKay's buntings. We sequenced 1123 base pairs of mitochondrial DNA for 40 individuals of each species and analysed a total of 913 AFLPs for 57 individuals. Both marker types suggested weak but significant genetic differentiation. Analysis of sequence data indicated divergence occurring when the current breeding range of McKay's buntings was a hill on the Beringian steppe (approximately 18,400 to approximately 73,700 years before present), suggesting that snow buntings were restricted to lower latitudes by ice sheets. Ancestral effective population size estimates indicate a founder event in McKay's buntings followed by an expansion and then a reduction in effective size. Rising sea levels and asymmetric hybridization from McKay's buntings into the postglacially-colonizing population of snow buntings could account for this reduction. Reproductive isolation is likely maintained through differential arrival dates on breeding grounds and the high breeding density of McKay's buntings. This recent, high-latitude divergence best fits a model of founder event speciation driven by vicariance and oscillations in habitat due to climate change.