Historical population size change of bowhead whales inferred from DNA sequence polymorphism data

Nucleotide sequence data from the mitochondrial control region were used from a phylogenetic context to investigate the long-term history of a population of bowhead whales (Balaena mysticetus). In addition, the coalescence time of these sequences was used to estimate the age of the inferred patterns of population size change. The results indicate that mitochondrial genetic polymorphism was not affected by a recent bottleneck that occurred near the turn of the 20th century, thereby preserving the signature of historical population size change in the mitochondrial genome. Further analysis showed that this population underwent an expansion initiated in the Middle to Late Pleistocene. As such, early Holocene changes in Arctic sea ice distribution appear to have had little influence on patterns of genetic variability in this population.     

Comparative phylogeography of two African carnivorans presumably introduced into Europe: disentangling natural versus human‐mediated dispersal across the Strait of Gibraltar

Aim Natural processes of colonization and human‐mediated introductions have shaped current patterns of biodiversity in the Mediterranean Basin. We use a comparative phylogeographic approach to investigate the genetic structure of Herpestes ichneumon and Genetta genetta (Carnivora) across the Strait of Gibraltar, and test for their supposedly contemporaneous introduction into Iberia. Location Mediterranean Basin and Africa. Methods We sequenced two mitochondrial fragments (cytochrome b and control region) of 91 (H. ichneumon) and 185 (G. genetta) individuals, including the sole archaeological record of G. genetta in Iberia, dating from the Muslim occupation. We used phylogenetic and tokogenetic methods, summary statistics, neutrality tests, geographic–genetic pairwise comparisons and coalescent estimates to explore the history of the two species in the Mediterranean Basin. Results In North Africa, an autochthonous (Clade I) and a western African mtDNA clade, coalescing in the Middle to Late Pleistocene, co‐occurred in both species. Only Clade I was present in Europe. In H. ichneumon, the European pool showed deep coalescence (median = 335 kyr) and high genetic differentiation and diversity compared with its North African counterpart, suggesting long‐term stability of female effective population size. In sharp contrast, G. genetta in Europe exhibited lower genetic diversity, weak differentiation with North Africa and recent demographic expansion; however, Andalusia and Catalonia (Spain) showed distinctly higher genetic diversity, and the archaeological specimen had the predominant European haplotype. Main conclusions The co‐occurrence of autochthonous and sub‐Saharan lineages in North Africa (1) supports a new, emerging biogeographic scenario in North Africa, and (2) suggests a connection through the Sahara, possibly from the Middle Pleistocene onwards. Our results refute the idea that H. ichneumon was introduced into Europe contemporaneously with G. genetta. Instead, they support a scenario of sweepstake dispersal during Late Pleistocene sea‐level fluctuations, followed by long‐term in situ evolution throughout the last glaciation cycles. Genetta genetta appears to have undergone a recent spread from at least two independent introduction ‘hotspots’ in Catalonia and Andalusia, possibly following antique trade routes and/or Muslim invasions. Despite their contrasting histories, the European gene pools of both species represent unusual cases leading to the preservation of autochthonous, North African lineages.     

Genetic admixture supports an ancient hybrid origin of the endangered Hawaiian duck

Speciation is regarded primarily as a bifurcation from an ancestral species into two distinct taxonomic units, but gene flow can create complex signals of phylogenetic relationships, especially among different loci. We evaluated several hypotheses that could account for phylogenetic discord between mitochondrial DNA (mtDNA) and nuclear DNA (nuDNA) within Hawaiian duck (Anas wyvilliana), including stochastic lineage sorting, mtDNA capture and widespread genomic introgression. Our results best support the hypothesis that the contemporary Hawaiian duck is descended from an ancient hybridization event between the mallard (Anas platyrhynchos) and Laysan duck (Anas laysanensis). Whereas mtDNA clearly shows a sister relationship between Hawaiian duck and mallard, nuDNA is consistent with a genetic mosaic with nearly equal contributions from Laysan duck and mallard. In addition, coalescent analyses suggest that gene flow from either mallard or Laysan duck, depending on the predefined tree topology, is necessary to explain contemporary genetic diversity in Hawaiian ducks, and these estimates are more consistent with ancient, rather than contemporary, hybridization. Time since divergence estimates suggest that the genetic admixture event occurred around the Pleistocene–Holocene boundary, which is further supported by circumstantial evidence from the Hawaiian subfossil record. Although the extent of reproductive isolation from either putative parental taxon is not currently known, these species are phenotypically, genetically and ecologically different, and they meet primary criteria used in avian taxonomy for species designation. Thus, the available data are consistent with an admixed origin and support the hypothesis that the Hawaiian duck may represent a young hybrid species.     

Species Delimitation--a Geneious plugin for the exploration of species boundaries

Species Delimitation is a plugin to the Geneious software to support the exploration of species boundaries in a gene tree. The user assigns taxa to putative species and the plugin computes statistics relating to the probability of the observed monophyly or exclusivity having occurred by chance in a coalescent process. It also assesses the within and between species genetic distances to infer the probability with which members of a putative species might be identified successfully with tree-based methods.     

Factors Influencing Ascertainment Bias of Microsatellite Allele Sizes: Impact on Estimates of Mutation Rates

Microsatellite loci play an important role as markers for identification, disease gene mapping, and evolutionary studies. Mutation rate, which is of fundamental importance, can be obtained from interspecies comparisons, which, however, are subject to ascertainment bias. This bias arises, for example, when a locus is selected on the basis of its large allele size in one species (cognate species 1), in which it is first discovered. This bias is reflected in average allele length in any noncognate species 2 being smaller than that in species 1. This phenomenon was observed in various pairs of species, including comparisons of allele sizes in human and chimpanzee. Various mechanisms were proposed to explain observed differences in mean allele lengths between two species. Here, we examine the framework of a single-step asymmetric and unrestricted stepwise mutation model with genetic drift. Analysis is based on coalescent theory. Analytical results are confirmed by simulations using the simuPOP software. The mechanism of ascertainment bias in this model is a tighter correlation of allele sizes within a cognate species 1 than of allele sizes in two different species 1 and 2. We present computations of the expected average allele size difference, given the mutation rate, population sizes of species 1 and 2, time of separation of species 1 and 2, and the age of the allele. We show that when the past demographic histories of the cognate and noncognate taxa are different, the rate and directionality of mutations affect the allele sizes in the two taxa differently from the simple effect of ascertainment bias. This effect may exaggerate or reverse the effect of difference in mutation rates. We reanalyze literature data, which indicate that despite the bias, the microsatellite mutation rate estimate in the ancestral population is consistently greater than that in either human or chimpanzee and the mutation rate estimate in human exceeds or equals that in chimpanzee with the rate of allele length expansion in human being greater than that in chimpanzee. We also demonstrate that population bottlenecks and expansions in the recent human history have little impact on our conclusions.     

Asymptotic Distributions of Coalescence Times and Ancestral Lineage Numbers for Populations with Temporally Varying Size

The distributions of coalescence times and ancestral lineage numbers play an essential role in coalescent modeling and ancestral inference. Both exact distributions of coalescence times and ancestral lineage numbers are expressed as the sum of alternating series, and the terms in the series become numerically intractable for large samples. More computationally attractive are their asymptotic distributions, which were derived in Griffiths (1984) for populations with constant size. In this article, we derive the asymptotic distributions of coalescence times and ancestral lineage numbers for populations with temporally varying size. For a sample of size n, denote by T_m the mth coalescent time, when m + 1 lineages coalesce into m lineages, and A_n(t) the number of ancestral lineages at time t back from the current generation. Similar to the results in Griffiths (1984), the number of ancestral lineages, A_n(t), and the coalescence times, T_m, are asymptotically normal, with the mean and variance of these distributions depending on the population size function, N(t). At the very early stage of the coalescent, when t → 0, the number of coalesced lineages n − A_n(t) follows a Poisson distribution, and as m → n, n(n ? 1)T_m/2N(0) follows a gamma distribution. We demonstrate the accuracy of the asymptotic approximations by comparing to both exact distributions and coalescent simulations. Several applications of the theoretical results are also shown: deriving statistics related to the properties of gene genealogies, such as the time to the most recent common ancestor (TMRCA) and the total branch length (TBL) of the genealogy, and deriving the allele frequency spectrum for large genealogies. With the advent of genomic-level sequencing data for large samples, the asymptotic distributions are expected to have wide applications in theoretical and methodological development for population genetic inference.     

Phylogeography of the Trans‐Volcanic bunchgrass lizard (Sceloporus bicanthalis) across the highlands of south‐eastern Mexico

We quantify the population divergence processes that shaped population genetic structure in the Trans‐Volcanic bunchgrass lizard (Sceloporus bicanthalis) across the highlands of south‐eastern Mexico. Multilocus genetic data from nine nuclear loci and mitochondrial (mt)DNA were used to estimate the population divergence history for 47 samples of S. bicanthalis. Bayesian clustering methods partitioned S. bicanthalis into three populations: (1) a southern population in Oaxaca and southern Puebla; (2) a population in western Puebla; and (3) a northern population with a broad distribution across Hidalgo, Puebla, and Veracruz. The multilocus nuclear data and mtDNA both supported a Late Pleistocene increase in effective population size, and the nuclear data revealed low levels of unidirectional gene flow from the widespread northern population into the southern and western populations. Populations of S. bicanthalis experienced different demographic histories during the Pleistocene, and phylogeographical patterns were similar to those observed in many co‐distributed highland taxa. Although we recommend continuing to recognize S. bicanthalis as a single species, future research on the evolution of viviparity could gain novel insights by contrasting physiological and genomic patterns among the different populations located across the highlands of south‐eastern Mexico. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2013, 110 , 852–865.     

Cooperation is less likely to evolve in a finite population with a highly skewed distribution of family size

In the context of the finitely repeated Prisoner's Dilemma with the possibility of cooperating or defecting each time, the strategy tit-for-tat (TFT) consists in cooperating the first time and copying the strategy previously used by the opponent the next times. Assuming random pairwise interactions in a finite population of always defecting individuals, TFT can be favoured by selection to go to fixation following its introduction as a mutant strategy. We deduce the condition for this to be the case under weak selection in the framework of a general reproduction scheme in discrete time. In fact, we show when and why the one-third rule for the evolution of cooperation holds, and how it extends to a more general rule. The condition turns out to be more stringent when the numbers of descendants left by the individuals from one time-step to the next may substantially differ. This suggests that the evolution of cooperation is made more difficult in populations with a highly skewed distribution of family size. This is illustrated by two examples.     

No resistance to male-killing Wolbachia after thousands of years of infection

Maternally transmitted male-killing endosymbionts can exert strong and relentless selection pressure on their hosts to evolve resistance to these infections. Surveys of current infection prevalence and mtDNA diversity indicate that Drosophila innubila is and has been infected with male-killing Wolbachia at moderate frequencies for extended evolutionary periods. Here, we use coalescent simulations to infer the minimum age of the Wolbachia infection in this species, and estimate that the infection is at least 15,000 and perhaps over 700,000 years old. We also surveyed this species for genetic variation for resistance to the male-killing effects of infection. Our surveys revealed no evidence for any resistance polymorphism, such that all flies are completely susceptible to male killing. Given the general assumption that Drosophila can be selected for anything, the lack of resistance, despite thousands of years of strong selection, is an apparent evolutionary conundrum. We hypothesize that resistance requires a mutation of major effect that acts early in development, and that the adverse pleiotropic consequences of such mutations in both infected and uninfected individuals may exceed the possible benefit to infected flies.