proseq: A software for preparation and evolutionary analysis of DNA sequence data sets

proseq is an integrated user‐friendly windows based program for convenient sequence editing and evolutionary analysis. It is designed to simplify preparation and analysis of DNA sequence data sets in population genetic, phylogenetic and molecular ecology studies. Sequence editor features include editing of chromatogram files, contig assembly, sequence alignment, translation and other utilities. Analysis features include calculation of genetic diversity, divergence, population subdivision and gene flow with permutation‐based significance testing and various tests of neutrality. A tool for coalescent simulations implements models with intragenic recombination, population subdivision and population growth.     

SPAms: A user-friendly software to simulate population genetics data under complex demographic models

Patterns of genetic diversity can be used to reconstruct populations demographic history. Thus, it is crucial to understand the statistical properties of data sets generated under different scenarios. SPAms is a graphical application that allows the simulation of data under a set of demographic models of increasing complexity: (i) population size change (ii) admixture and (iii) population structure (n-island or stepping stone model). SPAms generates data under the infinite site and the stepwise mutation models and computes a set of commonly used statistics. It should thus be a useful tool for both research and teaching purposes.     

Distinguishing migration from isolation: a Markov chain Monte Carlo approach

A Markov chain Monte Carlo method for estimating the relative effects of migration and isolation on genetic diversity in a pair of populations from DNA sequence data is developed and tested using simulations. The two populations are assumed to be descended from a panmictic ancestral population at some time in the past and may (or may not) after that be connected by migration. The use of a Markov chain Monte Carlo method allows the joint estimation of multiple demographic parameters in either a Bayesian or a likelihood framework. The parameters estimated include the migration rate for each population, the time since the two populations diverged from a common ancestral population, and the relative size of each of the two current populations and of the common ancestral population. The results show that even a single nonrecombining genetic locus can provide substantial power to test the hypothesis of no ongoing migration and/or to test models of symmetric migration between the two populations. The use of the method is illustrated in an application to mitochondrial DNA sequence data from a fish species: the threespine stickleback (Gasterosteus aculeatus).     

Demographic History and Population Structure of Blackfin Flounder (Glyptocephalus stelleri) in Japan Revealed by Mitochondrial Control Region Sequences

The demographic history and population genetic structure of the blackfin flounder (Glyptocephalus stelleri) along coastal regions of Japan were investigated. Genetic variation in DNA sequences was examined from the first hypervariable region of the mitochondrial DNA control region. A high level of haplotypic diversity (h = 0.99 ± 0.004) was detected, indicating a high level of intrapopulation genetic diversity. The starburst structure of the minimum spanning tree suggested a very recent origin for most haplotypes. The demographic history of G. stelleri was examined using neutrality tests and mismatch distribution analysis, which also indicated a Pleistocene population expansion at about 124,100–413,400 years ago. Hierarchical molecular variance analysis and conventional population Fst comparisons revealed no significant genetic differentiation throughout the range examined.     

Joint Estimation of Contamination,Error and Demography for Nuclear DNA from Ancient Humans

When sequencing an ancient DNA sample from a hominin fossil, DNA from present-day humans involved in excavation and extraction will be sequenced along with the endogenous material. This type of contamination is problematic for downstream analyses as it will introduce a bias towards the population of the contaminating individual(s). Quantifying the extent of contamination is a crucial step as it allows researchers to account for possible biases that may arise in downstream genetic analyses. Here, we present an MCMC algorithm to co-estimate the contamination rate, sequencing error rate and demographic parameters—including drift times and admixture rates—for an ancient nuclear genome obtained from human remains, when the putative contaminating DNA comes from present-day humans. We assume we have a large panel representing the putative contaminant population (e.g. European, East Asian or African). The method is implemented in a C++ program called ‘Demographic Inference with Contamination and Error’ (DICE). We applied it to simulations and genome data from ancient Neanderthals and modern humans. With reasonable levels of genome sequence coverage (>3X), we find we can recover accurate estimates of all these parameters, even when the contamination rate is as high as 50%.     

Population genetic structure of the point-head flounder, <Emphasis Type="Italic">Cleisthenes herzensteini</Emphasis>, in the Northwestern Pacific

Intraspecific phylogenies can provide useful insights into how populations have been shaped by historical and contemporary processes. To determine the population genetic structure and the demographic and colonization history of Cleisthenes herzensteini in the Northwestern Pacific, one hundred and twenty-one individuals were sampled from six localities along the coastal regions of Japan and the Yellow Sea of China. Mitochondrial DNA variation was analyzed using DNA sequence data from the 5′ end of control region. High levels of haplotype diversity (>0.96) were found for all populations, indicating a high level of genetic diversity. No pattern of isolation by distance was detected among the population differentiation throughout the examined range. Analyses of molecular variance (AMOVA) and the conventional population statistic Fst revealed no significant population genetic structure among populations. According to the exact test of differentiation among populations, the null hypothesis that C. herzensteini within the examined range constituted a non-differential mtDNA gene pool was accepted. The demographic history of C. herzensteini was examined using neutrality test and mismatch distribution analyses and results indicated Pleistocene population expansion (about 94–376 kya) in the species, which was consistent with the inference result of nested clade phylogeographical analysis (NCPA) showing contiguous range expansion for C. herzensteini. The lack of phylogeographical structure for the species may reflect a recent range expansion after the glacial maximum and insufficient time to attain migration-drift equilibrium.     

Domestication history in the Medicago sativa species complex: inferences from nuclear sequence polymorphism

DNA sequence polymorphism carries genealogical information and allows for testing hypotheses on selection and population history, especially through coalescent-based analysis. Understanding the evolutionary forces at work in plant domestication and subsequent selection is of critical importance for the management of genetic resources. In this study, we surveyed DNA sequence diversity at two assumed neutral nuclear loci in the wild-domesticated species complex of alfalfa (Medicago sativa L.). A high level of polymorphism was detected. The domesticated pool contains on average 31% less diversity than the wild pool, but with a high heterogeneity among loci. Coalescent simulations of the domestication process showed that this result cannot be explained by assuming a constant population size but is rather consistent with a demographic bottleneck during domestication. A very low level of divergence was detected between the wild and the domesticated forms as well as between the related subspecies of the M. sativa species complex. However, the originality of the Spanish wild populations, already observed based on mitochondrial DNA polymorphism, was confirmed. These results, together with patterns of intrapopulation polymorphism, suggest that nuclear sequence polymorphism could be a promising tool, complementary to mitochondrial DNA and phenotypic evaluations, to investigate historical demographic and evolutionary processes.     

DNA‐based approaches for evaluating historical demography in terrestrial vertebrates

Contemporary DNA sequences can provide information about the historical demography of a species. However, different molecular markers are informative under different circumstances. In particular, mitochondrial (mt)DNA is uniparentally inherited and haploid in most vertebrates and thus has a smaller effective population size than diploid, biparentally inherited nuclear (n)DNA. Here, we review the characteristics of mtDNA and nDNA in the context of historical demography. In particular, we address how their contrasting rates of evolution and sex‐biased dispersal can lead to different demographic inferences. We do so in the context of an extensive review of the vertebrate literature that describes the use of mtDNA and nDNA sequence data in demographic reconstruction. We discuss the effects of coalescence, effective population size, substitution rates, and sex‐biased dispersal on informative timeframes and expected patterns of genetic differentiation. We argue that mtDNA variationin species with male‐biased dispersal can imply deviations from neutrality that do not reflect actual population expansion or selection. By contrast, mtDNA can be more informative when coalescence has occurred within the recent past, which appears to be the case with many vertebrates. We also compare the application and interpretation of demographic and neutrality test statistics in historical demography studies. © 2014 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 112 , 367–386.     

Genetic variation and evolutionary demography of Fenneropenaeus chinensis populations, as revealed by the analysis of mitochondrial control region sequences

Genetic variation and evolutionary demography of the shrimp Fenneropenaeus chinensis were investigated using sequence data of the complete mitochondrial control region (CR). Fragments of 993 bp of the CR were sequenced for 93 individuals from five localities over most of the species' range in the Yellow Sea and the Bohai Sea. There were 84 variable sites defining 68 haplotypes. Haplotype diversity levels were very high (0.95 ± 0.03-0.99 ± 0.02) in F. chinensis populations, whereas those of nucleotide diversity were moderate to low (0.66 ± 0.36%-0.84 ± 0.46%). Analysis of molecular variance and conventional population statistics (F(ST) ) revealed no significant genetic structure throughout the range of F. chinensis. Mismatch distribution, estimates of population parameters and neutrality tests revealed that the significant fluctuations and shallow coalescence of mtDNA genealogies observed were coincident with estimated demographic parameters and neutrality tests, in implying important past-population size fluctuations or range expansion. Isolation with Migration (IM) coalescence results suggest that F. chinensis, distributed along the coasts of northern China and the Korean Peninsula (about 1000 km apart), diverged recently, the estimated time-split being 12,800 (7,400-18,600) years ago.     

Serial SimCoal: a population genetics model for data from multiple populations and points in time

SUMMARY: We present Serial SimCoal, a program that models population genetic data from multiple time points, as with ancient DNA data. An extension of SIMCOAL, it also allows simultaneous modeling of complex demographic histories, and migration between multiple populations. Further, we incorporate a statistical package to calculate relevant summary statistics, which, for the first time allows users to investigate the statistical power provided by, conduct hypothesis-testing with, and explore sample size limitations of ancient DNA data. AVAILABILITY: Source code and Windows/Mac executables at http://www.stanford.edu/group/hadlylab/ssc.html CONTACT: senka@stanford.edu.     

Contrasting selection modes at the Adh1 locus in outcrossing Miscanthus sinensis vs. inbreeding Miscanthus condensatus (Poaceae)

We estimated DNA sequence variation of the Adh1 locus in the outcrossing Miscanthus sinensis (Poaceae) and its close selfing relative, M. condensatus. Tajima's test of selection is significantly negative for both overall exons and replacement sites in M. sinensis. Among its entire sample, nucleotide diversity of nonsynonymous sites is higher than that of synonymous sites. A McDonald and Kreitman test of neutrality indicates an excess of intraspecific replacement polymorphisms, suggesting possible directional selection toward advantageous mutants. However, frequent intragenic recombination suggests both purifying and positive selection is unlikely. Recent demographic expansions coupled with relaxation of purifying selection may have resulted in elevated genetic diversity at the Adh1 locus as well as the trnL-trnF intergenic spacer of cpDNA in this outcrossing species. In contrast, low levels of genetic diversity were detected at both the Adh1 locus and the cpDNA spacer in M. condensatus, consistent with bottlenecks associated with selfing in all populations. While Tajima's D and Fu and Li's F statistics did not reveal deviation from neutrality at the Adh1 locus in M. condensatus, 12 replacements vs. 10 synonymous changes were detected. Based on pairwise comparisons of the d(N)/d(S) ratio, lineages of closely related populations of the species distributed along saline habitats appeared to be under directional selection.     

Patterns of Nucleotide Diversity at Photoperiod Related Genes in Norway Spruce [Picea abies (L.) Karst.]

The ability of plants to track seasonal changes is largely dependent on genes assigned to the photoperiod pathway, and variation in those genes is thereby important for adaptation to local day length conditions. Extensive physiological data in several temperate conifer species suggest that populations are adapted to local light conditions, but data on the genes underlying this adaptation are more limited. Here we present nucleotide diversity data from 19 genes putatively involved in photoperiodic response in Norway spruce (Picea abies). Based on similarity to model plants the genes were grouped into three categories according to their presumed position in the photoperiod pathway: photoreceptors, circadian clock genes, and downstream targets. An HKA (Hudson, Kreitman and Aquade) test showed a significant excess of diversity at photoreceptor genes, but no departure from neutrality at circadian genes and downstream targets. Departures from neutrality were also tested with Tajima''s D and Fay and Wu''s H statistics under three demographic scenarios: the standard neutral model, a population expansion model, and a more complex population split model. Only one gene, the circadian clock gene PaPRR3 with a highly positive Tajima''s D value, deviates significantly from all tested demographic scenarios. As the PaPRR3 gene harbours multiple non-synonymous variants it appears as an excellent candidate gene for control of photoperiod response in Norway spruce.     

SPLATCHE2: a spatially explicit simulation framework for complex demography, genetic admixture and recombination

SUMMARY: SPLATCHE2 is a program to simulate the demography of populations and the resulting molecular diversity for a wide range of evolutionary scenarios. The spatially explicit simulation framework can account for environmental heterogeneity and fluctuations, and it can manage multiple population sources. A coalescent-based approach is used to generate genetic markers mostly used in population genetics studies (DNA sequences, SNPs, STRs or RFLPs). Various combinations of independent, fully or partially linked genetic markers can be produced under a recombination model based on the ancestral recombination graph. Competition between two populations (or species) can also be simulated with user-defined levels of admixture between the two populations. SPLATCHE2 may be used to generate the expected genetic diversity under complex demographic scenarios and can thus serve to test null hypotheses. For model parameter estimation, SPLATCHE2 can easily be integrated into an Approximate Bayesian Computation (ABC) framework. Availability and implementation: SPLATCHE2 is a C++ program compiled for Windows and Linux platforms. It is freely available at www.splatche.com, together with its related documentation and example data. CONTACT: mathias.currat@unige.ch     

Differential population structuring and demographic history of two closely related fish species, Japanese sea bass (Lateolabrax japonicus) and spotted sea bass (Lateolabrax maculatus) in Northwestern Pacific

The Quaternary cold periods in the Northwestern Pacific are thought to have heavily influenced the amount and distribution of intraspecific genetic variation in marine fishes. To estimate the demographic history and genetic structure of Lateolabrax maculatus and L. japonicus in the Northwestern Pacific, 256 individuals were sampled from 19 localities throughout the distribution range of the two species. Mitochondrial DNA variation was analyzed using DNA sequence data from the cytochrome b gene and control region. Nucleotide diversity was much higher in L. japonicus (0.030) than in L. maculatus (0.012). The demographic history of the two species was examined using neutrality tests and mismatch distribution analyses and results indicated Pleistocene population expansion in both species. Estimates of population expansion time suggested earlier population expansion in L. japonicus than in L. maculatus. Molecular variance analyses showed differential genetic structuring for these two closely related species. The results indicated that L. japonicus is panmictic throughout its range. In contrast, populations of L. maculatus showed statistically significant levels of genetic structuring. Pattern of isolation by distance was observed in L. maculatus, suggesting that L. maculatus is in genetic equilibrium. In contrast, L. japonicus did not exhibit isolation by distance.     

Reconstructing the evolutionary history of an endangered subspecies across the changing landscape of the Great Central Valley of California

Identifying historic patterns of population genetic diversity and connectivity is a primary challenge in efforts to re‐establish the processes that have generated and maintained genetic variation across natural landscapes. The challenge of reconstructing pattern and process is even greater in highly altered landscapes where population extinctions and dramatic demographic fluctuations in remnant populations may have substantially altered, if not eliminated, historic patterns. Here, we seek to reconstruct historic patterns of diversity and connectivity in an endangered subspecies of woodrat that now occupies only 1–2 remnant locations within the highly altered landscape of the Great Central Valley of California. We examine patterns of diversity and connectivity using 14 microsatellite loci and sequence data from a mitochondrial locus and a nuclear intron. We reconstruct temporal change in habitat availability to establish several historical scenarios that could have led to contemporary patterns of diversity, and use an approximate Bayesian computation approach to test which of these scenarios is most consistent with our observed data. We find that the Central Valley populations harbour unique genetic variation coupled with a history of admixture between two well‐differentiated species of woodrats that are currently restricted to the woodlands flanking the Valley. Our simulations also show that certain commonly used analytical approaches may fail to recover a history of admixture when populations experience severe bottlenecks subsequent to hybridization. Overall our study shows the strength of combining empirical and simulation analyses to recover the history of populations occupying highly altered landscapes.     

Polymorphism Data Can Reveal the Origin of Species Abundance Statistics

What is the underlying mechanism behind the fat-tailed statistics observed for species abundance distributions? The two main hypotheses in the field are the adaptive (niche) theories, where species abundance reflects its fitness, and the neutral theory that assumes demographic stochasticity as the main factor determining community structure. Both explanations suggest quite similar species-abundance distributions, but very different histories: niche scenarios assume that a species population in the past was similar to the observed one, while neutral scenarios are characterized by strongly fluctuating populations. Since the genetic variations within a population depend on its abundance in the past, we present here a way to discriminate between the theories using the genetic diversity of noncoding DNA. A statistical test, based on the Fu-Li method, has been developed and enables such a differentiation. We have analyzed the results gathered from individual-based simulation of both types of histories and obtained clear distinction between the Fu-Li statistics of the neutral scenario and that of the niche scenario. Our results suggest that data for 10–50 species, with approximately 30 sequenced individuals for each species, may allow one to distinguish between these two theories.     

Chloroplast microsatellites reveal colonization and metapopulation dynamics in the Canary Island pine

Chloroplast microsatellites are becoming increasingly popular markers for population genetic studies in plants, but there has been little focus on their potential for demographic inference. In this work the utility of chloroplast microsatellites for the study of population expansions was explored. First, we investigated the power of mismatch distribution analysis and the F(S) test with coalescent simulations of different demographic scenarios. We then applied these methods to empirical data obtained for the Canary Island pine (Pinus canariensis). The results of the simulations showed that chloroplast microsatellites are sensitive to sudden population growth. The power of the F(S) test and accuracy of demographic parameter estimates, such as the time of expansion, were reduced proportionally to the level of homoplasy within the data. The analysis of Canary Island pine chloroplast microsatellite data indicated population expansions for almost all sample localities. Demographic expansions at the island level can be explained by the colonization of the archipelago by the pine, while population expansions of different ages in different localities within an island could be the result of local extinctions and recolonization dynamics. Comparable mitochondrial DNA sequence data from a parasite of P. canariensis, the weevil Brachyderes rugatus, supports this scenario, suggesting a key role for volcanism in the evolution of pine forest communities in the Canary Islands.     

Effect of different sampling strategies for a single geographic region in Yemen on standard genetic analyses of mitochondrial DNA sequence data

Collection of biological samples is the foundation of genetic studies ranging from estimation of genetic diversity to reconstruction of population history. Sample collections are intended to accurately represent the genetic, biological, ecological, cultural, geographic, and/or linguistic diversity of a particular region or population by providing a small, but representative, set of samples. In this study, we analyze human mitochondrial DNA variation in samples collected using four different sampling strategies to represent the same geographic region. Specifically, samples were collected from a village, a rural area, a regional clinic, and a national university in the governorate of Dhamar in Yemen. All samples were assayed for mitochondrial hypervariable region I DNA sequence variation and data were subjected to standard molecular genetic analyses. Our results suggest that analyses in which individual DNA sequences are explicitly compared or evaluated, e.g. phylogenetic and network analyses, may be more sensitive to sample collection design than analyses in which data are averaged across individuals or are analyzed more indirectly, e.g. summary statistics.     

HLA DNA sequence variation among human populations: molecular signatures of demographic and selective events

Molecular differences between HLA alleles vary up to 57 nucleotides within the peptide binding coding region of human Major Histocompatibility Complex (MHC) genes, but it is still unclear whether this variation results from a stochastic process or from selective constraints related to functional differences among HLA molecules. Although HLA alleles are generally treated as equidistant molecular units in population genetic studies, DNA sequence diversity among populations is also crucial to interpret the observed HLA polymorphism. In this study, we used a large dataset of 2,062 DNA sequences defined for the different HLA alleles to analyze nucleotide diversity of seven HLA genes in 23,500 individuals of about 200 populations spread worldwide. We first analyzed the HLA molecular structure and diversity of these populations in relation to geographic variation and we further investigated possible departures from selective neutrality through Tajima's tests and mismatch distributions. All results were compared to those obtained by classical approaches applied to HLA allele frequencies.Our study shows that the global patterns of HLA nucleotide diversity among populations are significantly correlated to geography, although in some specific cases the molecular information reveals unexpected genetic relationships. At all loci except HLA-DPB1, populations have accumulated a high proportion of very divergent alleles, suggesting an advantage of heterozygotes expressing molecularly distant HLA molecules (asymmetric overdominant selection model). However, both different intensities of selection and unequal levels of gene conversion may explain the heterogeneous mismatch distributions observed among the loci. Also, distinctive patterns of sequence divergence observed at the HLA-DPB1 locus suggest current neutrality but old selective pressures on this gene. We conclude that HLA DNA sequences advantageously complement HLA allele frequencies as a source of data used to explore the genetic history of human populations, and that their analysis allows a more thorough investigation of human MHC molecular evolution.