Likelihood‐based inference of population history from low‐coverage de novo genome assemblies

Short‐read sequencing technologies have in principle made it feasible to draw detailed inferences about the recent history of any organism. In practice, however, this remains challenging due to the difficulty of genome assembly in most organisms and the lack of statistical methods powerful enough to discriminate between recent, nonequilibrium histories. We address both the assembly and inference challenges. We develop a bioinformatic pipeline for generating outgroup‐rooted alignments of orthologous sequence blocks from de novo low‐coverage short‐read data for a small number of genomes, and show how such sequence blocks can be used to fit explicit models of population divergence and admixture in a likelihood framework. To illustrate our approach, we reconstruct the Pleistocene history of an oak‐feeding insect (the oak gallwasp Biorhiza pallida), which, in common with many other taxa, was restricted during Pleistocene ice ages to a longitudinal series of southern refugia spanning the Western Palaearctic. Our analysis of sequence blocks sampled from a single genome from each of three major glacial refugia reveals support for an unexpected history dominated by recent admixture. Despite the fact that 80% of the genome is affected by admixture during the last glacial cycle, we are able to infer the deeper divergence history of these populations. These inferences are robust to variation in block length, mutation model and the sampling location of individual genomes within refugia. This combination of de novo assembly and numerical likelihood calculation provides a powerful framework for estimating recent population history that can be applied to any organism without the need for prior genetic resources.     

A likelihood-based comparison of population histories in a parasitoid guild

Little is known about the stability of trophic relationships in complex natural communities over evolutionary timescales. Here, we use sequence data from 18 nuclear loci to reconstruct and compare the intraspecific histories of major Pleistocene refugial populations in the Middle East, the Balkans and Iberia in a guild of four Chalcid parasitoids (Cecidostiba fungosa, Cecidostiba semifascia, Hobbya stenonota and Mesopolobus amaenus) all attacking Cynipid oak galls. We develop a likelihood method to numerically estimate models of divergence between three populations from multilocus data. We investigate the power of this framework on simulated data, and—using triplet alignments of intronic loci—quantify the support for all possible divergence relationships between refugial populations in the four parasitoids. Although an East to West order of population divergence has highest support in all but one species, we cannot rule out alternative population tree topologies. Comparing the estimated times of population splits between species, we find that one species, M. amaenus, has a significantly older history than the rest of the guild and must have arrived in central Europe at least one glacial cycle prior to other guild members. This suggests that although all four species may share a common origin in the East, they expanded westwards into Europe at different times.     

Genealogical exclusivity in geographically proximate populations of Hypochilus thorelli Marx (Araneae,Hypochilidae) on the Cumberland Plateau of North America

The issue of sampling sufficiency is too infrequently explored in phylogeographical analysis, despite both theoretical work and analytical methods that stress the importance of sampling effort. Regarding the evolutionary pattern of reciprocal monophyly, both the probability of recovering this pattern and the possible inferences derived from this pattern, are highly contingent upon the density and geographical scale of sampling. Here, we present an empirical example that relates directly to this issue. We analyse genetic structure in the southern Appalachian spider Hypochilus thorelli, using an average sample of 5 mitochondrial DNA (mtDNA) sequences per location for 19 locations. All sampled sites are reciprocally monophyletic for mtDNA variation, even when separated by geographical distances as small as 5 km. For populations separated by greater geographical distances of 20-50 km, mtDNA sequences are not only exclusive, but are also highly divergent (uncorrected p-distances exceeding 5%). Although these extreme genealogical patterns are most seemingly consistent with a complete isolation model, both a coalescent method and nested cladistic analysis suggest that other restricted, but nonzero, gene flow models may also apply. Hypochilus thorelli appears to have maintained morphological cohesion despite this limited female-based gene flow, suggesting a pattern of stasis similar to that observed at higher taxonomic levels in Hypochilus.     

Refugial isolation and secondary contact in the dyeing poison frog Dendrobates tinctorius

Recent palaeoclimactic research suggests that fluctuating environmental conditions throughout the Pleistocene of Amazonia occurred with previously unrecognized frequency. This has resulted in a theoretical shift from glacially influenced fluctuations to those driven by precessional rhythms. This theoretical revolution has a profound impact on expectations of biotic diversity within biogeographical regions that have long been based on the idea of large-scale landscape fragmentation associated with increased aridity and glacial cycles. Generally speaking, this shifts phylogeographical expectations from that of (i) large areas of sympatry of closely related (but not sister) species whose origins lie in separate refugia, and current distributions are the results of cyclic connectivity of those two refugia (refuge hypothesis), to that of (ii) fine scale genetic structure, often associated with elevation, and divergence well below expected speciation levels [disturbance-vicariance (DV) hypothesis]. To date there have been few tests of the expectations of the DV hypothesis based on empirical studies of Neotropical floral and faunal communities. Herein we examine phylogeographical structure of Dendrobates tinctorius, an amphibian species endemic to the uplands of the eastern Guiana Shield, based on sampling of 114 individuals from 24 localities. Phylogenetic, nested clade, and dispersal-vicariance (DIVA) analyses of cytochrome b sequence data reveal the presence of two mitochondrial lineages that are associated with previously identified western and eastern uplands of this area. The geographical distribution of mitochondrial haplotypes and the results of DIVA and coalescent analyses suggest that there has been extensive secondary contact between these lineages indicating a complex history of connectivity between these western and eastern highlands, supporting the predictions of the DV hypothesis.     

Molecular genetic variation and population structure in morphologically differentiated cave and surface populations of the freshwater amphipod Gammarus minus

Gammarus minus is an important component of surface spring and cave ecosystems throughout Appalachia, and is a useful indicator of the hydrology and gene flow in freshwater communities. Gammarus minus populations occupying large cave passages (> 2 km) are usually troglomorphic, having reduced eyes, fewer ommatidia, larger body size, longer antennae, and reduced pigmentation relative to surface populations. We surveyed five cave stream and 10 surface spring populations for DNA sequence variation in the cytochrome c oxidase I (COI) and internal transcribed spacer 1 (ITS-1) genes with an aim towards characterizing phylogeographical structure and comparing the nature of genetic variation in cave vs. surface populations. Although standing variation at both loci was rather low within populations, a significant degree of divergence and spatial structuring of populations was observed. Levels of genetic variation within cave and spring populations differed substantially, with caves harbouring significantly less variation at the COI locus than surface springs. Codon usage bias was significantly lower in caves, indicating that cave streams harbour smaller and/or more recently colonized populations. Overall these data indicate limited gene flow among populations and suggest that the cave populations sampled in this study are prone to bottlenecks, possibly due to larger temperature fluctuations and more frequent incidence of drought conditions associated with these particular cave habitats.     

植物亲缘地理学的研究现状与发展趋势

新生代以来剧烈的地质和气候变化,特别是第四纪冰期和间冰期的反复交错对温带植物类群的地理分布格局、种群结构和物种分化造成了深刻影响。近20年来,随着分子标记技术的发展和分析方法的成熟,植物亲缘地理学利用基因谱系关系和在种群中的时空分布来追溯种群的进化历史,在探讨植物类群对地质气候变迁的响应以及种群分化的历史成因等方面开展了大量工作,取得了一些重要进展。在取得这些成绩的同时,作者也注意到了被研究者所忽略的一些问题。试图对这些问题进行初步的总结,着重于分析方法和几个常用软件的使用,并对今后的发展趋势进行了展望。     

Statistical phylogeography

While studies of phylogeography and speciation in the past have largely focused on the documentation or detection of significant patterns of population genetic structure, the emerging field of statistical phylogeography aims to infer the history and processes underlying that structure, and to provide objective, rather than ad hoc explanations. Methods for parameter estimation are now commonly used to make inferences about demographic past. Although these approaches are well developed statistically, they typically pay little attention to geographical history. In contrast, methods that seek to reconstruct phylogeographic history are able to consider many alternative geographical scenarios, but are primarily nonstatistical, making inferences about particular biological processes without explicit reference to stochastically derived expectations. We advocate the merging of these two traditions so that statistical phylogeographic methods can provide an accurate representation of the past, consider a diverse array of processes, and yet yield a statistical estimate of that history. We discuss various conceptual issues associated with statistical phylogeographic inferences, considering especially the stochasticity of population genetic processes and assessing the confidence of phylogeographic conclusions. To this end, we present some empirical examples that utilize a statistical phylogeographic approach, and then by contrasting results from a coalescent-based approach to those from Templeton's nested cladistic analysis (NCA), we illustrate the importance of assessing error. Because NCA does not assess error in its inferences about historical processes or contemporary gene flow, we performed a small-scale study using simulated data to examine how our conclusions might be affected by such unconsidered errors. NCA did not identify the processes used to simulate the data, confusing among deterministic processes and the stochastic sorting of gene lineages. There is as yet insufficient justification of NCA's ability to accurately infer or distinguish among alternative processes. We close with a discussion of some unresolved problems of current statistical phylogeographic methods to propose areas in need of future development.