Standard Population Genetics Research Articles

Stable population structure in Europe since the Iron Age, despite high mobility.

Ancient DNA research in the past decade has revealed that European population structure changed dramatically in the prehistoric period (14,000-3000 years before present, YBP), reflecting the widespread introduction of Neolithic farmer and Bronze Age Steppe ancestries. However, little is known about how population structure changed from the historical period onward (3000 YBP - present). To address this, we collected whole genomes from 204 individuals from Europe and the Mediterranean, many of which are the first historical period genomes from their region (e.g. Armenia and France). We found that most regions show remarkable inter-individual heterogeneity. At least 7% of historical individuals carry ancestry uncommon in the region where they were sampled, some indicating cross-Mediterranean contacts. Despite this high level of mobility, overall population structure across western Eurasia is relatively stable through the historical period up to the present, mirroring geography. We show that, under standard population genetics models with local panmixia, the observed level of dispersal would lead to a collapse of population structure. Persistent population structure thus suggests a lower effective migration rate than indicated by the observed dispersal. We hypothesize that this phenomenon can be explained by extensive transient dispersal arising from drastically improved transportation networks and the Roman Empire's mobilization of people for trade, labor, and military. This work highlights the utility of ancient DNA in elucidating finer scale human population dynamics in recent history.

Can Unlikely Neanderthal Chloride Channel CLC-2 Gene Variants Provide Insights in Modern Human Infertility?

Neanderthals, although well adapted to local environments, were rapidly replaced by anatomically modern humans (AMH) for unknown reasons. Genetic information on Neanderthals is limited restricting applicability of standard population genetics. Here, we apply a novel combination of restricted genetic analyses on preselected physiological key players (ion channels), electrophysiological analyses of gene variants of unclear significance expressed in Xenopus laevis oocytes using two electrode voltage clamp and transfer of results to AMH genetics. Using genetic screening in infertile men identified a loss of CLC-2 associated with sperm deficiency. Increased genetic variation caused functionally impaired Neanderthals CLC-2 channels. Increased genetic variation could reflect an adaptation to different local salt supplies at the cost of reduced sperm density. Interestingly and consistent with this hypothesis, lack of CLC-2 protein in a patient associates with high blood K+ concentration and azoospermia.

High-throughput sequencing on preservative ethanol is effective at jointly examining infraspecific and taxonomic diversity, although bioinformatics pipelines do not perform equally.

High‐throughput sequencing of amplicons (HTSA) has been proposed as an effective approach to evaluate taxonomic and genetic diversity at the same time. However, there are still uncertainties as to how the results produced by different bioinformatics treatments impact the conclusions drawn on biodiversity and population genetics indices.We evaluated the ability of six bioinformatics pipelines to recover taxonomic and genetic diversity from HTSA data obtained from controlled assemblages. To that end, 20 assemblages were produced using 354 colonies of Botrylloides spp., sampled in the wild in ten marinas around Brittany (France). We used DNA extracted from preservative ethanol (ebDNA) after various time of storage (3, 6, and 12 months), and from a bulk of preserved specimens (bulkDNA). DNA was amplified with primers designed for targeting this ascidian genus. Results obtained from HTSA data were compared with Sanger sequencing on individual zooids (i.e., individual barcoding).Species identification and relative abundance determined with HTSA data from either ebDNA or bulkDNA were similar to those obtained with traditional individual barcoding. However, after 12 months of storage, the correlation between HTSA and individual‐based data was lower than after shorter durations. The six bioinformatics pipelines were able to depict accurately the genetic diversity using standard population genetics indices (HS and FST), despite producing false positives and missing rare haplotypes. However, they did not perform equally and dada2 was the only pipeline able to retrieve all expected haplotypes.This study showed that ebDNA is a nondestructive alternative for both species identification and haplotype recovery, providing storage does not last more than 6 months before DNA extraction. Choosing the bioinformatics pipeline is a matter of compromise, aiming to retrieve all true haplotypes while avoiding false positives. We here recommend to process HTSA data using dada2, including a chimera‐removal step. Even if the possibility to use multiplexed primer sets deserves further investigation to expand the taxonomic coverage in future similar studies, we showed that primers targeting a particular genus allowed to reliably analyze this genus within a complex community.

Population data for 23 autosomal STR loci in White British population.

We report allele frequencies and standard population genetics parameters for 23 forensic autosomal STR loci tested among White British population. In addition, we use generated allele frequency data to compare the extent of genetic variation with 19 other European populations.

A population genetics perspective on the evolutionary histories of three clonal, endemic, and dominant grass species of the Qinghai-Tibet Plateau: Orinus (Poaceae).

We performed analyses of amplified fragment length polymorphism (AFLP) in order to characterize the evolutionary history of Orinus according to its population genetic structure, as well as to investigate putative hybrid origins of O. intermedius and to provide additional insights into relationships among species. The genus Orinus comprises three clonal grasses that are dominant species within xeric alpine grasslands of the Qinghai–Tibet Plateau (QTP). Here, we used eight selectively obtained primer pairs of EcoRI/MseI to perform amplifications in 231 individuals of Orinus representing 48 populations and all three species. We compared our resulting data to genetic models of hybridization using a Bayesian algorithm within NewHybrids software. We determined that genetic variation in Orinus was 56.65% within populations while the among‐species component was 30.04% using standard population genetics statistics. Nevertheless, we detected that species of Orinus were clustered into three highly distinct genetic groups corresponding to classic species identities. Our results suggest that there is some introgression among species. Thus, we tested explicit models of hybridization using a Bayesian approach within NewHybrids software. However, O. intermedius likely derives from a common ancestor with O. kokonoricus and is probably not the result of hybrid speciation between O. kokonoricus and O. thoroldii. We suspect that recent isolation of species of Orinus in allopatry via vicariance may explain the patterns in diversity that we observed, and this is corroborated by a Mantel test that showed significant positive correlation between geographic and genetic distance (r = 0.05, p < 0.05). Recent isolation may explain why Orinus differs from many other clonal species by exhibiting the highest diversity within populations rather than among them.

Spatial genetics of brown hares (Lepus europaeus Pallas, 1778) from Turkey: Different gene pool architecture on either side of the Bosphorus?

We determined allelic variation at eleven microsatellite loci to study the effect of the Bosphorus, the Sea of Marmara, and the Dardanelles on the population genetic architecture of brown hares, Lepus europaeus, from Anatolia and Turkish Thrace in the southeastern most Balkans. The latter region was connected with Anatolia during periods of the late Pleistocene and the early Holocene (until ca. 8000 y BP) enabling natural gene flow. Following a recent mitochondrial (mt) DNA-based phylogeographic model that suggested unidirectional ancient gene flow from Anatolia to the SE Balkans, we expected lower genetic diversity in hares from Turkish Thrace than Anatolia. The Turkish Thrace population might have been established only relatively recently as natural northwestward expansion of a long existing Anatolian population. Alternatively, it might have received late Pleistocene or Holocene gene flow by hares from eastern Europe from north of the Black Sea, even after the formation of the marine barrier. This might have produced a notable genetic distinction between the hares from either side of the Bosphorus, the Sea of Marmara, and the Dardanelles. We found significantly lower allelic richness in hares from Turkish Thrace than Anatolia, very little genetic admixture in the Turkish Thrace population, and only one allele at low frequency in Turkish Thrace that was not observed in Anatolian hares. On the contrary, many alleles of Anatolian hares were not recovered in hares from Turkish Thrace. Nevertheless, genetic divergence (Fst, Cavalli-Sforza & Edwards distances, AMOVA) between the hares from either side of the current marine barrier was low (Fst values <10%). All findings accord to our initial hypothesis that the hares from Turkish Thrace represent merely a somewhat genetically depleted version of the Anatolian hares. Seemingly, no strong drift has occurred in the Turkish Thrace population since the break-up of the Holocene landbridge, and there is no indication of gene flow from other regions in (eastern) Europe into the Turkish Thrace population. Spatial genetic approaches and standard population genetics (BAPS, GENELAND, STRUCTURE, F-statistics, Cavalli-Sforza & Edwards distances) indicated significant spatial genetic differentiation of hares across Anatolia, but also on a very low level. Overall, our results are concordant with the latest mtDNA based model of ancient gene flow in hares only from Anatolia to the southeastern Balkans across the Pleistocene and early Holocene landbridge; this implies that most or all of the Balkans have been colonized by hares exclusively originating eventually from Anatolia.

An ancient and a recent colonization of islands by an Australian sap‐feeding insect

AbstractAimTo assess the genetic structure, biogeography, and the potential for speciation, of a highly host‐specific insect pest with mainland and island populations.LocationEast coast of Australia, Lord Howe Island (LHI), and New Zealand.MethodsWe focussed on Mycopsylla fici, a plant sap‐feeding insect host‐specific to the fig tree Ficus macrophylla. We genotyped 152 insects from across the natural and extended host plant range at 14 microsatellite loci and analysed the data using standard population genetics statistics, Discriminant Analysis of Principal Components, genetic autocorrelation, and two Bayesian clustering approaches.ResultsGenetic analyses revealed that the northeastern Australian mainland population (Brisbane) is the centre of genetic diversity. Northeastern and southeastern (Sydney) mainland populations are genetically differentiated and interconnected in a stepping‐stone pattern. The LHI population is the most distinct genetically and Bayesian estimates indicated that the most recent colonization occurred c. 2,000‐17,500 years ago from a northeastern mainland origin. In contrast, the New Zealand population is little differentiated from the Sydney population and probably diverged by colonization within the past 200 years.Main conclusionsThe strong differentiation in nuclear microsatellites mirrors previous evidence for divergence of the LHI population from both mtDNA and endosymbiont DNA. The LHI population may be undergoing speciation from the mainland populations, with an oceanic barrier to gene flow. In contrast, the geographically isolated population in Auckland represents a far more recent colonization reflecting the contemporary naturalization of the plant host in New Zealand.

Comparative phylogeography of Aedes mosquitoes and the role of past climatic change for evolution within Africa.

The study of demographic processes involved in species diversification and evolution ultimately provides explanations for the complex distribution of biodiversity on earth, indicates regions important for the maintenance and generation of biodiversity, and identifies biological units important for conservation or medical consequence. African and forest biota have both received relatively little attention with regard to understanding their diversification, although one possible mechanism is that this has been driven by historical climate change. To investigate this, we implemented a standard population genetics approach along with Approximate Bayesian Computation, using sequence data from two exon‐primed intron‐crossing (EPIC) nuclear loci and mitochondrial cytochrome oxidase subunit I, to investigate the evolutionary history of five medically important and inherently forest dependent mosquito species of the genus Aedes. By testing different demographic hypotheses, we show that Aedes bromeliae and Aedes lilii fit the same model of lineage diversification, admixture, expansion, and recent population structure previously inferred for Aedes aegypti. In addition, analyses of population structure show that Aedes africanus has undergone lineage diversification and expansion while Aedes hansfordi has been impacted by population expansion within Uganda. This congruence in evolutionary history is likely to relate to historical climate‐driven habitat change within Africa during the late Pleistocene and Holocene epoch. We find differences in the population structure of mosquitoes from Tanzania and Uganda compared to Benin and Uganda which could relate to differences in the historical connectivity of forests across the continent. Our findings emphasize the importance of recent climate change in the evolution of African forest biota.

Dinucleotide (CA)n tandem repeats on the human X-chromosome and the history of the Mediterranean populations

Background: Tandem repeats (STRs) are genomic markers of particular interest in forensic and population genetics. Most of the population data currently available correspond to the variation of STRs of forensic panels, which barely include dinucleotide tandem repeats.Aims: The aim of the study is to test the usefulness of a battery of dinucleotide STRs on the X chromosome for population and forensic studies.Subjects and methods: A total of 672 individuals from 12 Mediterranean populations and two external references were analysed for 15 X-STR following the instructions of the commercial company and using control DNA from the CEPH centre whose sequences are published in GenBank. Genotypic results were analysed using standard population genetics methods including estimates of linkage disequilibrium, population structure and gene flow. Common forensic efficiency parameters were calculated.Results: The analysed X-STRs show high values of genetic diversity, comparable to other STRs of more common use. No significant associations between markers were found. A slight population structure was detected between the two shores of the Mediterranean. The X-STRs studied here present a similar degree of variability to that of other X-STRs used in forensics.Conclusion: Tandem-repeated dinucleotides are a good tool for evidencing population differences here. Forensic parameters indicate that the dinucleotide X-STRs are suitable for forensic use.

Population Genetics with Fluctuating Population Sizes

Standard neutral population genetics theory with a strictly fixed population size has important limitations. An alternative model that allows independently fluctuating population sizes and reproduces the standard neutral evolution is reviewed. We then study a situation such that the competing species are neutral at the equilibrium population size but population size fluctuations nevertheless favor fixation of one species over the other. In this case, a separation of timescales emerges naturally and allows adiabatic elimination of a fast population size variable to deduce the fluctuations-induced selection dynamics near the equilibrium population size. The results highlight the incompleteness of the standard population genetics with a strictly fixed population size.

Latitudinal structured populations of the Mexican wild squash Cucurbita argyrosperma subsp. sororia revealed by microsatellite markers

Crop wild relatives represent an important agronomic resource for crop improvement and biodiversity conservation. The wild squash Cucurbita argyrosperma subsp. sororia. (Cucurbitaceae) has been considered the wild ancestor of cultivated forms of C. argyrosperma. In order to characterise the geographic patterns of genetic variation in this wild cucurbit and to identify priority areas for conservation, we analysed the genetic diversity and structure of natural populations along the Mexican Pacific coast. By using 14 polymorphic microsatellites, we genotyped 378 individuals sampled from 61 locations. Standard population genetics analyses and group testing were conducted on the genotypes with the aid of principal coordinate analysis and Bayesian analysis. Overall, we found an average of 12.3 alleles per locus and an expected heterozygosity of 0.756. We found greater genetic diversity in southern populations. The fixation index was 0.113, suggesting a mixed mating system. The Mantel test revealed a minor distance effect on genetic differentiation between individuals (r = 0.321). Finally, we found three main groups of populations arranged in a mostly latitudinal pattern, from Sinaloa (north-west) to Oaxaca–Guerrero (south-east). The greater genetic diversity and heterogeneity among southern populations (Guerrero–Oaxaca), suggests that this region is an important centre of diversity of this wild squash with important implications for conservation.

Weakly Deleterious Mutations and Low Rates of Recombination Limit the Impact of Natural Selection on Bacterial Genomes.

ABSTRACTFree-living bacteria are usually thought to have large effective population sizes, and so tiny selective differences can drive their evolution. However, because recombination is infrequent, “background selection” against slightly deleterious alleles should reduce the effective population size (Ne) by orders of magnitude. For example, for a well-mixed population with 1012 individuals and a typical level of homologous recombination (r/m = 3, i.e., nucleotide changes due to recombination [r] occur at 3 times the mutation rate [m]), we predict that Ne is <107. An argument for high Ne values for bacteria has been the high genetic diversity within many bacterial “species,” but this diversity may be due to population structure: diversity across subpopulations can be far higher than diversity within a subpopulation, which makes it difficult to estimate Ne correctly. Given an estimate of Ne, standard population genetics models imply that selection should be sufficient to drive evolution if Ne × s is >1, where s is the selection coefficient. We found that this remains approximately correct if background selection is occurring or when population structure is present. Overall, we predict that even for free-living bacteria with enormous populations, natural selection is only a significant force if s is above 10−7 or so.

Epigenetics and Evolution: Transposons and the Stochastic Epigenetic Modification Model

AbstractIn addition to genetic variation, epigenetic variation and transposons can greatly affect the evolutionary fitnesses landscape and gene expression. Previously we proposed a mathematical treatment of a general epigenetic variation model that we called Stochastic Epigenetic Modification (SEM) model. In this study we follow up with a special case, the Transposon Silencing Model (TSM), with, once again, emphasis on quantitative treatment. We have investigated the evolutionary effects of epigenetic changes due to transposon (T) insertions; in particular, we have considered a typical gene locus A and postulated that (i) the expression level of gene A depends on the epigenetic state (active or inactive) of a cis-located transposon element T, (ii) stochastic variability in the epigenetic silencing of T occurs only in a short window of opportunity during development, (iii) the epigenetic state is then stable during further development, and (iv) the epigenetic memory is fully reset at each generation. We develop the model using two complementary approaches: a standard analytical population genetics framework (diffusion equations) and Monte-Carlo simulations. Both approaches led to similar estimates for the probability of fixation and time of fixation of locus TA with initial frequency P in a randomly mating diploid population of effective size Ne . We have ascertained the effect that ρ, the probability of transposon modification during the developmental window, has on the population (species). One of our principal conclusions is that as ρ increases, the pattern of fixation of the combined TA locus goes from “neutral” to “dominant” to “over-dominant”. We observe that, under realistic values of ρ, epigenetic modifications can provide an efficient mechanism for more rapid fixation of transposons and cis-located gene alleles. The results obtained suggest that epigenetic silencing, even if strictly transient (being reset at each generation), can still have significant macro-evolutionary effects. Importantly, this conclusion also holds for the static fitness landscape. To the best of our knowledge, no previous analytical modeling has treated stochastic epigenetic changes during a window of opportunity.

Phenotypic plasticity, the Baldwin effect, and the speeding up of evolution: the computational roots of an illusion.

An increasing number of dissident voices claim that the standard neo-Darwinian view of genes as 'leaders' and phenotypes as 'followers' during the process of adaptive evolution should be turned on its head. This idea is older than the rediscovery of Mendel's laws of inheritance, with the turn-of-the-twentieth-century notion eventually labeled as the 'Baldwin effect' as one of the many ways in which the standard neo-Darwinian view can be turned around. A condition for this effect is that environmentally induced variation such as phenotypic plasticity or learning is crucial for the initial establishment of a trait. This gives the additional time for natural selection to act on genetic variation and the adaptive trait can be eventually encoded in the genotype. An influential paper published in the late 1980s claimed the Baldwin effect to happen in computer simulations, and avowed that it was crucial to solve a difficult adaptive task. This generated much excitement among scholars in various disciplines that regard neo-Darwinian accounts to explain the evolutionary emergence of high-order phenotypic traits such as consciousness or language almost hopeless. Here, we use analytical and computational approaches to show that a standard population genetics treatment can easily crack what the scientific community has granted as an unsolvable adaptive problem without learning. Evolutionary psychologists and linguists have invoked the (claimed) Baldwin effect to make wild assertions that should not be taken seriously. What the Baldwin effect needs are plausible case-histories.

Riverine habitat specificity constrains dispersion in a Neotropical fish (Characidae) along Southern Brazilian drainages

Freshwater fishes often display a marked phylogeographic structure strongly associated with historical and ecological changes in the aquatic environment. Different ecological conditions in the same river drainage may act as permeable barriers to dispersion and gene flow. Previous studies recognized two discrete spatial components for the ichthyofauna in the freshwater coastal drainages of southern Brazil: the lowland fish fauna in the lagoons and the fish fauna of the rivers flowing in the valleys. In order to test if the coastal lagoons may limit the dispersion of a riverine species, we describe the phylogeographic structure among populations of Cyanocharax itaimbe, a species endemic to this region. We analysed 55 specimens characterized for two mitochondrial and one nuclear genes. Sequences were analysed using gene trees and species tree approaches, together with standard population genetics methods. Molecular analyses indicated three evolutionary groups which diverged from each other between an estimated 1,600,000 and 450,000 years before the present. However, two currently isolated river systems share the same evolutionary clade, whereas a single drainage contains two different lineages. Our results indicate strong genetic structure among populations along with generally conserved morphology. The strong genetic structure among populations living in the same drainage system may be explained by ecological differences between lagoons and rivers (or palaeochannels) that act as barriers to dispersion.

Purifying Selection, Drift, and Reversible Mutation with Arbitrarily High Mutation Rates

Some species exhibit very high levels of DNA sequence variability; there is also evidence for the existence of heritable epigenetic variants that experience state changes at a much higher rate than sequence variants. In both cases, the resulting high diversity levels within a population (hyperdiversity) mean that standard population genetics methods are not trustworthy. We analyze a population genetics model that incorporates purifying selection, reversible mutations, and genetic drift, assuming a stationary population size. We derive analytical results for both population parameters and sample statistics and discuss their implications for studies of natural genetic and epigenetic variation. In particular, we find that (1) many more intermediate-frequency variants are expected than under standard models, even with moderately strong purifying selection, and (2) rates of evolution under purifying selection may be close to, or even exceed, neutral rates. These findings are related to empirical studies of sequence and epigenetic variation.

Parentage and sibship inference from markers in polyploids

Many plants and some animal species are polyploids. Nondisomically inherited markers (e.g. microsatellites) in such species cannot be analysed directly by standard population genetics methods developed for diploid species. One solution is to transform the polyploid codominant genotypes to pseudodiploid-dominant genotypes, which can then be analysed by standard methods for various purposes such as spatial genetic structure, individual relatedness and relationship. Although this data transformation approach has been used repeatedly in the literature, no systematic study has been conducted to investigate how efficient it is, how much marker information is lost and thus how much analysis accuracy is reduced. More specifically, it is unknown whether or not the transformed data can be used to infer parentage and sibship jointly, and how different sampling schemes (number and polymorphism of markers, number of individuals) and ploidy level affect the inference accuracy. This study analyses both simulated and empirical data to examine the effects of polyploid levels, actual pedigree structures and marker number and polymorphism on the accuracy of joint parentage and sibship assignments in polyploid species. We show that sibship, parentage and selfing rates in polyploids can be inferred accurately from a typical set of microsatellite loci. We also show that inferences can be substantially improved by allowing for a small genotyping error rate to accommodate the distortion in assumed Mendelian inheritance of the converted markers when large sibship groups are involved. The results are discussed in the context of polyploid data analysis in molecular ecology.

The genetic diversity of Sardinian myrtle (Myrtus communis L.) populations

Background: The myrtle ( Myrtus communis ) is a common shrub widespread in the Mediterranean Basin. Its fruit and leaves exhibit antioxidant, antibacterial and antifungal properties, and are used for their content of essential oils and for their medicinal properties, but most commonly as an ingredient in locally made liquor. The uncontrolled exploitation of natural stands has reduced both the species' geographical coverage and the size of individual populations. The selection of genotypes for controlled cultivation requires a characterization of the genetic diversity present both within and between populations. Results: Genotypic variation was evaluated using ISSR profiling and genetic diversity characterized using standard population genetics approaches. Two major clusters were identified: one capturing all the candidate cultivars selected from various Sardinian localities, and the other wild individuals collected from Asinara, Corsica and Surigheddu. A moderate level of gene flow between the Sardinian and Corsican populations was identified. Discriminant analysis of principal components revealed a level of separation among the wild populations, confirming the population structure identified by the clustering methods. Conclusions: The wild accessions were well differentiated from the candidate cultivars. The level of genetic variability was high. The genetic data were compatible with the notion that myrtle has a mixed pollination system, including both out-pollination by insects and self-pollination. The candidate cultivars are suggested to represent an appropriate basis for directed breeding.

Using neutral cline decay to estimate contemporary dispersal: a generic tool and its application to a major crop pathogen

Dispersal is a key parameter of adaptation, invasion and persistence. Yet standard population genetics inference methods hardly distinguish it from drift and many species cannot be studied by direct mark-recapture methods. Here, we introduce a method using rates of change in cline shapes for neutral markers to estimate contemporary dispersal. We apply it to the devastating banana pest Mycosphaerella fijiensis, a wind-dispersed fungus for which a secondary contact zone had previously been detected using landscape genetics tools. By tracking the spatio-temporal frequency change of 15 microsatellite markers, we find that σ, the standard deviation of parent–offspring dispersal distances, is 1.2 km/generation1/2. The analysis is further shown robust to a large range of dispersal kernels. We conclude that combining landscape genetics approaches to detect breaks in allelic frequencies with analyses of changes in neutral genetic clines offers a powerful way to obtain ecologically relevant estimates of dispersal in many species.

Bayesian computation via empirical likelihood.

Approximate Bayesian computation has become an essential tool for the analysis of complex stochastic models when the likelihood function is numerically unavailable. However, the well-established statistical method of empirical likelihood provides another route to such settings that bypasses simulations from the model and the choices of the approximate Bayesian computation parameters (summary statistics, distance, tolerance), while being convergent in the number of observations. Furthermore, bypassing model simulations may lead to significant time savings in complex models, for instance those found in population genetics. The Bayesian computation with empirical likelihood algorithm we develop in this paper also provides an evaluation of its own performance through an associated effective sample size. The method is illustrated using several examples, including estimation of standard distributions, time series, and population genetics models.