Search

Inference of demographic histories of species and populations is one of the central problems in population genetics. It is usually stated as an optimization problem: find a model's parameters that maximize a certain log-likelihood. This log-likelihood is often expensive to evaluate in terms of time and hardware resources, critically more so for larger population counts. Although genetic algorithm-based solution has proven efficient for demographic inference in the past, it struggles to deal with log-likelihoods in the setting of more than three populations. Different tools are therefore needed to handle such scenarios. We introduce a new optimization pipeline for demographic inference with time consuming log-likelihood evaluations. It is based on Bayesian optimization, a prominent technique for optimizing expensive black box functions. Comparing to the existing widely used genetic algorithm solution, we demonstrate new pipeline's superiority in the limited time budget setting with four and five populations, when using the log-likelihoods provided by the moments tool., G3: Genes, Genomes, Genetics, 13 (7), ISSN:2160-1836

The establishment of European colonies across the world had important demographic consequences because it brought together diverse and distant civilizations for the first time. One clear example of this phenomenon is observed in the Canary Islands. The modern Canarian population is mainly the result of the admixture of natives of North African origin and European colonizers. However, additional migratory flows reached the islands due to the importation of enslaved Africans to cultivate sugarcane and the intense commercial contact with the American continent. In this review, we evaluate how the genetic analysis of indigenous, historical and current populations has provided a glimpse into the Canary Islands’ complex genetic composition. We show that each island subpopulation’s characterization is needed to fully disentangle the demographic history of the Canarian archipelago. Finally, we discuss what research avenues remain to be explored to improve our knowledge of the impact that the European colonization had on its native population.

The common hermit crab Pagurus bernhardus (Linnaeus, 1758) is an abundant and ecologically important benthic crustacean in the northeastern Atlantic Ocean. While this species has been intensively studied in terms of its ecology, physiology, behavior, and larval development, knowledge about its population structure and demographic history is still lacking. We examined, for the first time, the genetic variability of P. bernhardus by analyzing two mitochondrial gene fragments (CO1 and 16S) from more than 150 specimens collected from various locations from the Iberian Peninsula to Norway. Our results provide evidence for a significant genetic structure according to the sampled regions for both genetic markers. Furthermore, a comprehensive demographic history reconstruction, mainly based on neutrality tests and a Bayesian Skyline Plot (CO1), revealed a recent demographic expansion of P. bernhardus that preceded the Last Glacial Maximum. Such pattern of retrieved demographic trend could have been likely a successive process to historical contraction of the species into potential climate refugia within the surveyed geographic spectrum.

Motivation Population admixture is an important subject in population genetics. Inferring population demographic history with admixture under the so-called admixture network model from population genetic data is an established problem in genetics. Existing admixture network inference approaches work with single genetic polymorphisms. While these methods are usually very fast, they do not fully utilize the information [e.g. linkage disequilibrium (LD)] contained in population genetic data. Results In this article, we develop a new admixture network inference method called GTmix. Different from existing methods, GTmix works with local gene genealogies that can be inferred from population haplotypes. Local gene genealogies represent the evolutionary history of sampled haplotypes and contain the LD information. GTmix performs coalescent-based maximum likelihood inference of admixture networks with inferred local genealogies based on the well-known multispecies coalescent (MSC) model. GTmix utilizes various techniques to speed up the likelihood computation on the MSC model and the optimal network search. Our simulations show that GTmix can infer more accurate admixture networks with much smaller data than existing methods, even when these existing methods are given much larger data. GTmix is reasonably efficient and can analyze population genetic datasets of current interests. Availability and implementation The program GTmix is available for download at: https://github.com/yufengwudcs/GTmix. Supplementary information Supplementary data are available at Bioinformatics online.

Evolutionary processes related to climatic changes and ecological factors, such as microhabitat affinities and food specialization, can be important contributors to phylogeographic discordance between codistributed and related species. Here, we evaluate the evolutionary histories of two cactophilic and codistributed Drosophila species (Diptera: Drosophilidae) from South America, Drosophila antonietae and Drosophila meridionalis, where they use mainly Cereus hildmaniannus (Cactoideae: Cereeae) as host, using mitochondrial DNA sequences and species distribution modeling. The diversification of both species was estimated during the Pleistocene. For both species, the distribution of suitable areas through the Last Glacial period to the present showed a similar dynamic from Andes Valley through east and through the Paraná-Paraguay river basin to the Atlantic coastline. The current distribution of D. antonietae was influenced by demographic expansion and putative migration route from northwest to south and then to coast, with two genetic incipient groups with bidirectional genetic flow between them. For D. meridionalis, we suggested a migration route from south to north as well as to coast, with three genetic groups deeply structured with no evidence of demographic expansion. Our comparative results showed that the Quaternary paleoclimatic dynamic has had a similar role in both species (displacement of the high suitability areas) with similar routes but in different directions. Additionally, the Araucaria forest represents a putative biogeographic barrier for Drosophila species and also for host C. hildmaniannus. The phylogeographical differences between these species related to geographical distribution, genetic structure, and demographic history could be explained for differences to adaptation and plasticity to explore a new host.

Chromosomal evolution is widely considered to be an important driver of speciation, as karyotypic reorganization can bring about the establishment of reproductive barriers between incipient species. One textbook example for genetic mechanisms of speciation are large-scale chromosomal rearrangements such as Robertsonian (Rb) fusions, a common class of structural variants that can drastically change the recombination landscape by suppressing crossing-over and influence gene expression by altering regulatory networks. Here, we explore the population structure and demographic patterns of a well-known house mouse Rb system in the Aeolian archipelago in Southern Italy using genome-wide data. By analyzing chromosomal regions characterized by different levels of recombination, we trace the evolutionary history of a set of Rb chromosomes occurring in different geographical locations and test whether chromosomal fusions have a single shared origin or occurred multiple times. Using a combination of phylogenetic and population genetic approaches, we find support for multiple, independent origins of three focal Rb chromosomes. The elucidation of the demographic patterns of the mouse populations within the Aeolian archipelago shows that an interplay between fixation of newly formed Rb chromosomes and hybridization events has contributed to shaping their current karyotypic distribution. Overall, our results illustrate that chromosome structure is much more dynamic than anticipated and emphasize the importance of large-scale chromosomal translocations in speciation.

The prehistoric demography of human populations is an essential piece of information for illustrating our evolution. Despite its importance and the advancement of ancient DNA studies, our knowledge of human evolution is still limited, which is also the case for relatively recent population dynamics during and around the Holocene. Here, we inferred detailed demographic histories from 1 to 40 ka for 24 population samples using an improved model-flexible method with 36 million genome-wide noncoding CpG sites. Our results showed many population growth events that were likely due to the Neolithic Revolution (i.e., the shift from hunting and gathering to agriculture and settlement). Our results help to provide a clearer picture of human prehistoric demography, confirming the significant impact of agriculture on population expansion, and provide new hypotheses and directions for future research.

The relative evolutionary roles of adaptive and non-adaptive processes remain a central question in population genetics. Resolution of this debate has been difficult as an appropriate null model incorporating... The question of the relative evolutionary roles of adaptive and nonadaptive processes has been a central debate in population genetics for nearly a century. While advances have been made in the theoretical development of the underlying models, and statistical methods for estimating their parameters from large-scale genomic data, a framework for an appropriate null model remains elusive. A model incorporating evolutionary processes known to be in constant operation, genetic drift (as modulated by the demographic history of the population) and purifying selection, is lacking. Without such a null model, the role of adaptive processes in shaping within- and between-population variation may not be accurately assessed. Here, we investigate how population size changes and the strength of purifying selection affect patterns of variation at “neutral” sites near functional genomic components. We propose a novel statistical framework for jointly inferring the contribution of the relevant selective and demographic parameters. By means of extensive performance analyses, we quantify the utility of the approach, identify the most important statistics for parameter estimation, and compare the results with existing methods. Finally, we reanalyze genome-wide population-level data from a Zambian population of Drosophila melanogaster, and find that it has experienced a much slower rate of population growth than was inferred when the effects of purifying selection were neglected. Our approach represents an appropriate null model, against which the effects of positive selection can be assessed.

Natural populations of Populus euphratica in northwest China are divided into four clades exhibiting strong geographical distribution patterns. A total of 38 single nucleotide polymorphisms were associated with salinity tolerance, located in or near 82 genes., Populus euphratica is a dominant tree species in desert riparian forests and possesses extraordinary adaptation to salinity stress. Exploration of its genomic variation and molecular underpinning of salinity tolerance is important for elucidating population evolution and identifying stress-related genes. Here, we identify approximately 3.15 million single nucleotide polymorphisms using whole-genome resequencing. The natural populations of P. euphratica in northwest China are divided into four distinct clades that exhibit strong geographical distribution patterns. Pleistocene climatic fluctuations and tectonic deformation jointly shaped the extant genetic patterns. A seed germination rate-based salinity tolerance index was used to evaluate seed salinity tolerance of P. euphratica and a genome-wide association study was implemented. A total of 38 single nucleotide polymorphisms were associated with seed salinity tolerance and were located within or near 82 genes. Expression profiles showed that most of these genes were regulated under salt stress, revealing the genetic complexity of seed salinity tolerance. Furthermore, DEAD-box ATP-dependent RNA helicase 57 and one undescribed gene (CCG029559) were demonstrated to improve the seed salinity tolerance in transgenic Arabidopsis. These results provide new insights into the demographic history and genetic architecture of seed salinity tolerance in desert poplar.

Neutral genetic diversity across the genome is determined by the complex interplay of mutation, demographic history, and natural selection. While the direct action of natural selection is limited to functional loci across the genome, its impact can have effects on nearby neutral loci due to genetic linkage. These effects of selection at linked sites, referred to as genetic hitchhiking and background selection (BGS), are pervasive across natural populations. However, only recently has there been a focus on the joint consequences of demography and selection at linked sites, and some empirical studies have come to apparently contradictory conclusions as to their combined effects. To understand the relationship between demography and selection at linked sites, we conducted an extensive forward simulation study of BGS under a range of demographic models. We found that the relative levels of diversity in BGS and neutral regions vary over time and that the initial dynamics after a population size change are often in the opposite direction of the long-term expected trajectory. Our detailed observations of the temporal dynamics of neutral diversity in the context of selection at linked sites in nonequilibrium populations provide new intuition about why patterns of diversity under BGS vary through time in natural populations and help reconcile previously contradictory observations. Most notably, our results highlight that classical models of BGS are poorly suited for predicting diversity in nonequilibrium populations.

Pangolins are among the most critically endangered animals due to heavy poaching and worldwide trafficking. However, their demographic histories and the genomic consequences of their recent population declines remain unknown. We generated high-quality de novo reference genomes for critically endangered Malayan (Manis javanica, MJ) and Chinese (M. pentadactyla, MP) pangolins and re-sequencing population genomic data from 74 MJs and 23 MPs. We recovered the population identities of illegally traded pangolins and previously unrecognized genetic populations that should be protected as evolutionarily distinct conservation units. Demographic reconstruction suggested environmental changes have resulted in a population size fluctuation of pangolins. Additionally, recent population size declines due to human activities have resulted in an increase in inbreeding and genetic load. Deleterious mutations were enriched in genes related to cancer/diseases and cholesterol homeostasis, which may have increased their susceptibility to diseases and decreased their survival potential to adapt to environmental changes and high-cholesterol diets. This comprehensive study provides not only high-quality pangolin reference genomes, but also valuable information concerning the driving factors of long-term population size fluctuations and the genomic impact of recent population size declines due to human activities, which is essential for pangolin conservation management and global action planning.

The North American mountain goat (Oreamnos americanus) is an iconic alpine species that faces stressors from climate change, industrial development, and recreational activities. This species’ phylogenetic position within the Caprinae lineage has not been resolved and their phylogeographic history is dynamic and controversial. Genomic data could be used to address these questions and provide valuable insights to conservation and management initiatives. We sequenced short-read genomic libraries constructed from a DNA sample of a 2.5-year-old female mountain goat at 80X coverage. We improved the short-read assembly by generating Chicago library data and scaffolding using the HiRise approach. The final assembly was 2,506 Mbp in length with an N50 of 66.6 Mbp, which is within the length range and in the upper quartile for N50 published ungulate genome assemblies. Comparative analysis identified 84 gene families unique to the mountain goat. The species demographic history in terms of effective population size generally mirrored climatic trends over the past one hundred thousand years and showed a sharp decline during the last glacial maximum. This genome assembly will provide a reference basis for future population and comparative genomic analyses.

Kangaroo rats in the genus Dipodomys are found in a variety of habitat types in western North America, including deserts, arid and semi-arid grasslands, and scrublands. Many Dipodomys species are experiencing strong population declines due to increasing habitat fragmentation, with two species listed as federally endangered. The precarious state of many Dipodomys populations, including those occupying extreme environments, make species of this genus valuable subjects for studying the impacts of habitat degradation and fragmentation on population genomic patterns and for characterizing the genomic bases of adaptation to harsh conditions. To facilitate exploration of such questions, we assembled and annotated a reference genome for the banner-tailed kangaroo rat (D. spectabilis) using PacBio HiFi sequencing reads, providing a more contiguous genomic resource than two previously assembled Dipodomys genomes. Using the HiFi data for D. spectabilis and publicly available sequencing data for two other Dipodomys species (D. ordii and D. stephensi), we demonstrate the utility of this new assembly for studies of congeners by conducting inference of historic effective population sizes (Ne) and linking these patterns to the species’ current extinction risk statuses. The genome assembly presented here will serve as a valuable resource for population and conservation genomic studies of Dipodomys species, comparative genomic research within mammals and rodents, and investigations into genomic adaptation to extreme environments and changing landscapes.Significance statementKangaroo rats in the genus Dipodomys occur in a wide variety of habitat types, ranging from scrublands to arid deserts, and are increasingly impacted by habitat fragmentation with populations of many species in strong decline. To facilitate population and conservation genomic studies of Dipodomys species, we generated the first reference genome assembly for the extensively studied banner-tailed kangaroo rat (D. spectabilis) from long read PacBio sequencing data. The genome assembly presented here will serve as a valuable resource for studies of Dipodomys species—which have long served as ecological and physiological models for the study of osmoregulation—comparative genomic surveys of mammals and rodents, and investigations into genomic adaptation to extreme environments and changing landscapes.

The southwestern and central United States serve as an ideal region to test alternative hypotheses regarding biotic diversification. Genomic data can now be combined with sophisticated computational models to quantify the impacts of paleoclimate change, geographic features, and habitat heterogeneity on spatial patterns of genetic diversity. In this study, we combine thousands of genotyping-by-sequencing (GBS) loci with mtDNA sequences (ND1) from the Texas horned lizard (Phrynosoma cornutum) to quantify relative support for different catalysts of diversification. Phylogenetic and clustering analyses of the GBS data indicate support for at least three primary populations. The spatial distribution of populations appears concordant with habitat type, with desert populations in AZ and NM showing the largest genetic divergence from the remaining populations. The mtDNA data also support a divergent desert population, but other relationships differ and suggest mtDNA introgression. Genotype–environment association with bioclimatic variables supports divergence along precipitation gradients more than along temperature gradients. Demographic analyses support a complex history, with introgression and gene flow playing an important role during diversification. Bayesian multispecies coalescent analyses with introgression (MSci) analyses also suggest that gene flow occurred between populations. Paleo-species distribution models support two southern refugia that geographically correspond to contemporary lineages. We find that divergence times are underestimated and population sizes are overestimated when introgression occurred and is ignored in coalescent analyses, and furthermore, inference of ancient introgression events and demographic history is sensitive to inclusion of a single recently admixed sample. Our analyses cannot refute the riverine barrier or glacial refugia hypotheses. Results also suggest that populations are continuing to diverge along habitat gradients. Finally, the strong evidence of admixture, gene flow, and mtDNA introgression among populations suggests that P. cornutum should be considered a single widespread species under the General Lineage Species Concept.

Previous studies of the brown bear (Ursus arctos) on Hokkaido Island, Japan, have detected three geographically distinct subpopulations representing different mitochondrial lineages and shown that gene flow between subpopulations has occurred due to male-biased dispersal. In this study, we determined whole-genomic sequences for six Hokkaido brown bears and analyzed these data along with previously published genomic sequences of 17 brown bears from other parts of the world. We found that the Hokkaido population is genetically distinct from the other populations, keeping genetic diversity higher than the endangered populations in western Europe but lower than most populations on the continents. A reconstruction of historical demography showed no increase in population size for the Hokkaido population during the Eemian interglacial period (130,000–114,000 years ago). In a phylogenetic analysis of the autosomal data, the Hokkaido population formed a clade distinct from North American and European populations, showing that it has maintained genetic diversity independently from continental populations following geographical isolation on the island. This autosomal genetic similarity contrasts with the geographically separate mitochondrial lineages on Hokkaido and indicates the occurrence of male-driven gene flow between subpopulations.

The Patterson F- and D-statistics are commonly used measures for quantifying population relationships and for testing hypotheses about demographic history. These statistics make use of allele frequency information across populations to infer different aspects of population history, such as population structure and introgression events. Inclusion of related or inbred individuals can bias such statistics, which may often lead to the filtering of such individuals. Here, we derive statistical properties of the F- and D-statistics, including their biases due to the inclusion of related or inbred individuals, their variances, and their corresponding mean squared errors. Moreover, for those statistics that are biased, we develop unbiased estimators and evaluate the variances of these new quantities. Comparisons of the new unbiased statistics to the originals demonstrates that our newly derived statistics often have lower error across a wide population parameter space. Furthermore, we apply these unbiased estimators using several global human populations with the inclusion of related individuals to highlight their application on an empirical dataset. Finally, we implement these unbiased estimators in open-source software package funbiased for easy application by the scientific community.

A long-standing enigma concerns the geographic and ecological origins of the intensively studied vinegar fly, Drosophila melanogaster. This globally distributed human commensal is thought to originate from sub-Saharan Africa, yet until recently, it had never been reported from undisturbed wilderness environments that could reflect its precommensal niche. Here, we document the collection of 288 D. melanogaster individuals from multiple African wilderness areas in Zambia, Zimbabwe, and Namibia. The presence of D. melanogaster in these remote woodland environments is consistent with an ancestral range in southern-central Africa, as opposed to equatorial regions. After sequencing the genomes of 17 wilderness-collected flies collected from Kafue National Park in Zambia, we found reduced genetic diversity relative to town populations, elevated chromosomal inversion frequencies, and strong differences at specific genes including known insecticide targets. Combining these genomes with existing data, we probed the history of this species’ geographic expansion. Demographic estimates indicated that expansion from southern-central Africa began ∼13,000 years ago, with a Saharan crossing soon after, but expansion from the Middle East into Europe did not begin until roughly 1,800 years ago. This improved model of demographic history will provide an important resource for future evolutionary and genomic studies of this key model organism. Our findings add context to the history of D. melanogaster, while opening the door for future studies on the biological basis of adaptation to human environments.

Southeast Asian rainforests at upper hill elevations are increasingly vulnerable to degradation because most lowland forest areas have been converted to different land uses. As such, understanding the genetics of upper hill species is becoming more crucial for their future management and conservation. Shorea platyclados is an important, widespread upper hill dipterocarp in Malaysia. To elucidate the genetic structure of S. platyclados and ultimately provide guidelines for a conservation strategy for this species, we carried out a comprehensive study of the genetic diversity and demographic history of S. platyclados. Twenty-seven populations of S. platyclados across its range in Malaysia were genotyped at 15 polymorphic microsatellite loci and sequenced at seven noncoding chloroplast DNA (cpDNA) regions. A total of 303 alleles were derived from the microsatellite loci, and 29 haplotypes were identified based on 2892 bp of concatenated cpDNA sequences. The populations showed moderately high genetic diversity (mean HE = 0.680 for microsatellite gene diversity and HT = 0.650 for total haplotype diversity) and low genetic differentiation (FST = 0.060). Bayesian clustering divided the studied populations into two groups corresponding to western and eastern Malaysia. Bottleneck analysis did not detect any recent bottleneck events. Extended Bayesian skyline analyses showed a model of constant size for the past population history of this species. Based on our findings, priority areas for in situ and ex situ conservation and a minimum population size are recommended for the sustainable utilization of S. platyclados.

Short-scale local adaptation is a complex process involving selection, migration, and drift. The expected effects on the genome are well grounded in theory but examining these on an empirical level has proven difficult, as it requires information about local selection, demographic history, and recombination rate variation. Here, we use locally adapted and phenotypically differentiated Arabidopsis lyrata populations from two altitudinal gradients in Norway to test these expectations at the whole-genome level. Demography modeling indicates that populations within the gradients diverged

Background and AimsPolyploidy is an important driver of plant diversification and adaptation to novel environments. As a consequence of genome doubling, polyploids often exhibit greater colonizing ability or occupy a wider ecological niche than diploids. Although elevation has been traditionally considered as a key driver structuring ploidy variation, we do not know if environmental and phenotypic differentiation among ploidy cytotypes varies along an elevational gradient. Here, we tested for the consequences of genome duplication on genetic diversity, phenotypic variation and habitat preferences on closely related diploid and tetraploid populations that coexist along approx. 2300 m of varying elevation.MethodsWe sampled and phenotyped 45 natural diploid and tetraploid populations of Arabidopsis arenosa in one mountain range in Central Europe (Western Carpathians) and recorded abiotic and biotic variables at each collection site. We inferred genetic variation, population structure and demographic history in a sub-set of 29 populations genotyped for approx. 36 000 single nucleotide polymorphisms.Key ResultsWe found minor effects of polyploidy on colonization of alpine stands and low genetic differentiation between the two cytotypes, mirroring recent divergence of the polyploids from the local diploid lineage and repeated reticulation events among the cytotypes. This pattern was corroborated by the absence of ecological niche differentiation between the two cytotypes and overall phenotypic similarity at a given elevation.ConclusionsThe case of A. arenosa contrasts with previous studies that frequently showed clear niche differentiation between cytotypes. Our work stresses the importance of considering genetic structure and past demographic processes when interpreting the patterns of ploidy distributions, especially in species that underwent recent polyploidization events.

Salter-Townshend and Myers present an open source tool for modelling multi-way admixture events using dense haplotype data. Their Hidden Markov Model approach is scalable to thousands of samples and, unlike existing methods..., We present an algorithm for inferring ancestry segments and characterizing admixture events, which involve an arbitrary number of genetically differentiated groups coming together. This allows inference of the demographic history of the species, properties of admixing groups, identification of signatures of natural selection, and may aid disease gene mapping. The algorithm employs nested hidden Markov models to obtain local ancestry estimation along the genome for each admixed individual. In a range of simulations, the accuracy of these estimates equals or exceeds leading existing methods. Moreover, and unlike these approaches, we do not require any prior knowledge of the relationship between subgroups of donor reference haplotypes and the unseen mixing ancestral populations. Our approach infers these in terms of conditional “copying probabilities.” In application to the Human Genome Diversity Project, we corroborate many previously inferred admixture events (e.g., an ancient admixture event in the Kalash). We further identify novel events such as complex four-way admixture in San-Khomani individuals, and show that Eastern European populations possess 1−3% ancestry from a group resembling modern-day central Asians. We also identify evidence of recent natural selection favoring sub-Saharan ancestry at the human leukocyte antigen (HLA) region, across North African individuals. We make available an R and C++ software library, which we term MOSAIC (which stands for MOSAIC Organizes Segments of Ancestry In Chromosomes).

The archaic ancestry present in the human genome has captured the imagination of both scientists and the wider public in recent years. This excitement is the result of new studies pushing the envelope of what we can learn from the archaic genetic information that has survived for over 50,000 years in the human genome. Here, we review the most recent ten years of literature on the topic of archaic introgression, including the current state of knowledge on Neanderthal and Denisovan introgression, as well as introgression from other as-yet unidentified archaic populations. We focus this review on four topics: 1) a reimagining of human demographic history, including evidence for multiple admixture events between modern humans, Neanderthals, Denisovans, and other archaic populations; 2) state-of-the-art methods for detecting archaic ancestry in population-level genomic data; 3) how these novel methods can detect archaic introgression in modern African populations; and 4) the functional consequences of archaic gene variants, including how those variants were co-opted into novel function in modern human populations. The goal of this review is to provide a simple-to-access reference for the relevant methods and novel data, which has changed our understanding of the relationship between our species and its siblings. This body of literature reveals the large degree to which the genetic legacy of these extinct hominins has been integrated into the human populations of today.

The tight interaction between pathogens and their hosts results in reciprocal selective forces that impact the genetic diversity of the interacting species. The footprints of this selection differ between pathosystems because of distinct life-history traits, demographic histories, or genome architectures. Here, we studied the genome-wide patterns of genetic diversity of 22 isolates of the causative agent of the corn smut disease, Ustilago maydis, originating from five locations in Mexico, the presumed center of origin of this species. In this species, many genes encoding secreted effector proteins reside in so-called virulence clusters in the genome, an arrangement that is so far not found in other filamentous plant pathogens. Using a combination of population genomic statistical analyses, we assessed the geographical, historical and genome-wide variation of genetic diversity in this fungal pathogen.We report evidence of two partially admixed subpopulations that are only loosely associated with geographic origin. Using the multiple sequentially Markov coalescent model, we inferred the demographic history of the two pathogen subpopulations over the last 0.5 million years. We show that both populations experienced a recent strong bottleneck starting around 10,000 years ago, coinciding with the assumed time of maize domestication. While the genome average genetic diversity is low compared to other fungal pathogens, we estimated that the rate of non-synonymous adaptive substitutions is three times higher in genes located within virulence clusters compared to non-clustered genes, including non-clustered effector genes. These results highlight the role that these singular genomic regions play in the evolution of this pathogen.Significance statementThe maize pathogen Ustilago maydis is a model species to study fungal cell biology and biotrophic host-pathogen interactions. Population genetic studies of this species, however, were so far restricted to using a few molecular markers, and genome-wide comparisons involved species that diverged more than 20 million years ago. Here, we sequenced the genomes of 22 Mexican U. maydis isolates to study the recent evolutionary history of this species. We identified two co-existing populations that went through a recent bottleneck and whose divergence date overlaps with the time of maize domestication. Contrasting the patterns of genetic diversity in different categories of genes, we further showed that effector genes in virulence clusters display a high rate of adaptive mutations, highlighting the importance of these effector arrangements for the adaptation of U. maydis to its host.

The cohabitation of Drosophila melanogaster with humans is nearly ubiquitous. Though it has been well established that this fly species originated in sub-Saharan Africa, and only recently has spread globally, many details of its swift expansion remain unclear. Elucidating the demographic history of D. melanogaster provides a unique opportunity to investigate how human movement might have impacted patterns of genetic diversity in a commensal species, as well as providing neutral null models for studies aimed at identifying genomic signatures of local adaptation. Here, we use whole-genome data from five populations (Africa, North America, Europe, Central Asia, and the South Pacific) to carry out demographic inferences, with particular attention to the inclusion of migration and admixture. We demonstrate the importance of these parameters for model fitting and show that how previous estimates of divergence times are likely to be significantly underestimated as a result of not including them. Finally, we discuss how human movement along early shipping routes might have shaped the present-day population structure of D. melanogaster.

For organisms to survive and prosper in a harsh environment, particularly under rapid climate change, poses tremendous challenges. Recent studies have highlighted the continued loss of megafauna in terrestrial ecosystems and the subsequent surge of small mammals, such as rodents, bats, lagomorphs, and insectivores. However, the ecological partitioning of these animals will likely lead to large variation in their responses to environmental change. In the present study, we investigated the evolutionary history and genetic adaptations of white-bellied rats (Niviventer Marshall, 1976), which are widespread in the natural terrestrial ecosystems in Asia but also known as important zoonotic pathogen vectors and transmitters. The southeastern Qinghai-Tibet Plateau was inferred as the origin center of this genus, with parallel diversification in temperate and tropical niches. Demographic history analyses from mitochondrial and nuclear sequences of Niviventer demonstrated population size increases and range expansion for species in Southeast Asia, and habitat generalists elsewhere. Unexpectedly, population increases were seen in N. eha, which inhabits the highest elevation among Niviventer species. Genome scans of nuclear exons revealed that among the congeneric species, N. eha has the largest number of positively selected genes. Protein functions of these genes are mainly related to olfaction, taste, and tumor suppression. Extensive genetic modification presents a major strategy in response to global changes in these alpine species.

Genetic variation and phylogeographic studies have been crucial for understanding mechanisms of speciation. We analyzed genetic variation and phylogeography to reconstruct the demographic history of the Rivoli's Hummingbird (Eugenes fulgens) species complex and also evaluated their morphological differentiation. This widely distributed species inhabits the highlands of Mexico and northern Central America, with 2 subspecies separated by the Isthmus of Tehuantepec (west: E. f. fulgens, east: E. f. viridiceps). We surveyed genetic variation in 2 mitochondrial DNA markers (mtDNA, with 129 individuals) and nuclear DNA (6 microsatellites, with 85 individuals). We also inferred the demographic history, estimated divergence times, and analyzed morphological variation using 470 vouchered specimens. We modeled the current potential distribution of the species using ecological niche modeling and projected it into the past to model the effects of the Pleistocene climatic cycles. Haplotype networks, pairwise FST comparisons, AMOVA, and morphological analysis revealed differences between geographically isolated populations separated by the Isthmus of Tehuantepec (IT; corresponding to the 2 recognized subspecies: fulgens and viridiceps), and by the Motagua-Polochic-Jocotan (MPJ) system fault. Demographic scenarios revealed a contraction in distribution during the last interglacial, and expansion during the Last Glacial Maximum (LGM) with little change since the LGM. Divergence between groups separated by the Isthmus of Tehuantepec ∼59,600 yr ago occurred in the presence of gene flow, suggesting that the Isthmus of Tehuantepec is a semipermeable barrier to gene flow. STRUCTURE analyses of microsatellite data detected 3 genetically differentiated groups. Several results fit a model of recent lineage divergence, including a significant signal of genetic differentiation, demographic expansion, decreased gene flow from past to present, and northward expansion during the LGM and contraction during the interglacial periods. We conclude that the genetic differentiation of E. fulgens in the Madrean Pine-Oak Woodlands resulted from recent geographical isolation of populations separated by natural barriers (IT and MPJ).LAY SUMMARYThe present work allowed us to reconstruct the phylogeographic pattern and evolutionary history of Eugenes fulgens in the highlands of Mesoamerica.It revealed the presence of 3 main lineages: populations west of the Isthmus of Tehuantepec (E. fulgens) and 2 populations east of the isthmus (E. viridiceps), geographically isolated from each other by the Isthmus of Tehuantepec and by the Motagua-Polochic-Jocotan system fault.Our results support the hypothesis that the divergence and demographic expansion within the E. fulgens species complex are associated with the Pleistocene glacial–interglacial cycles.

Background The demographic history of any population is imprinted in the genomes of the individuals that make up the population. One of the most popular and convenient representations of genetic information is the allele frequency spectrum (AFS), the distribution of allele frequencies in populations. The joint AFS is commonly used to reconstruct the demographic history of multiple populations, and several methods based on diffusion approximation (e.g., ∂a∂i) and ordinary differential equations (e.g., moments) have been developed and applied for demographic inference. These methods provide an opportunity to simulate AFS under a variety of researcher-specified demographic models and to estimate the best model and associated parameters using likelihood-based local optimizations. However, there are no known algorithms to perform global searches of demographic models with a given AFS. Results Here, we introduce a new method that implements a global search using a genetic algorithm for the automatic and unsupervised inference of demographic history from joint AFS data. Our method is implemented in the software GADMA (Genetic Algorithm for Demographic Model Analysis, https://github.com/ctlab/GADMA). Conclusions We demonstrate the performance of GADMA by applying it to sequence data from humans and non-model organisms and show that it is able to automatically infer a demographic model close to or even better than the one that was previously obtained manually. Moreover, GADMA is able to infer multiple demographic models at different local optima close to the global one, providing a larger set of possible scenarios to further explore demographic history.

Larger populations are expected to have larger genetic diversity. However, as pointed out by Lewontin in 1974, the range of population sizes exceeds the range of genetic diversity by many orders of magnitude (a.k.a. "Lewontin's paradox," LP). The reasons for LP remain obscure. Here, This paper reports an extreme case of LP in astronomically large populations of the ubiquitous unicellular marine phytoplankton species Emiliania huxleyi (Haptophyta)-the species that accounts for 10-20% of primary productivity in the oceans and its blooms are so extensive that they are visible from space. This study demonstrates that despite the wide distribution and enormous population size, the world-wide sample of E. huxleyi strains with sequenced genomes represents a single cohesive species and contains surprisingly limited genetic diversity (π ∼ 0.006 per silent site). The patterns of polymorphism reveal even larger populations in the past, and frequent recombination (ρ ∼ 0.006) throughout the genome, ruling out demographic history and asexual reproduction as possible causes of low polymorphism in E. huxleyi. Natural selection wiping out genetic diversity at linked sites (a.k.a. "genetic draft") must be strong and frequent to account for low polymorphism in E. huxleyi. This study sheds the first light on poorly understood evolutionary genetic processes in astronomically large populations of marine microplankton.

As one of the most commonly utilized organisms in the study of local adaptation, an accurate characterization of the demographic history of Drosophila melanogaster remains as an important research question. This owes both to the inherent interest in characterizing the population history of this model organism, as well as to the well-established importance of an accurate null demographic model for increasing power and decreasing false positive rates in genomic scans for positive selection. Although considerable attention has been afforded to this issue in non-African populations, less is known about the demographic history of African populations, including from the ancestral range of the species. While qualitative predictions and hypotheses have previously been forwarded, we here present a quantitative model fitting of the population history characterizing both the ancestral Zambian population range as well as the subsequently colonized west African populations, which themselves served as the source of multiple non-African colonization events. We here report the split time of the West African population at 72 kya, a date corresponding to human migration into this region as well as a period of climatic changes in the African continent. Furthermore, we have estimated population sizes at this split time. These parameter estimates thus represent an important null model for future investigations in to African and non-African D. melanogaster populations alike.

Variation in LINE composition is one of the major determinants for the substantial size and structural differences among vertebrate genomes. In particular, the larger genomes of mammals are characterized by hundreds of thousands of copies from a single LINE clade, L1, whereas nonmammalian vertebrates possess a much greater diversity of LINEs, yet with orders of magnitude less in copy number. It has been proposed that such variation in copy number among vertebrates is due to differential effect of LINE insertions on host fitness. To investigate LINE selection, we deployed a framework of demographic modeling, coalescent simulations, and probabilistic inference against population-level whole-genome data sets for four model species: one population each of threespine stickleback, green anole, and house mouse, as well as three human populations. Specifically, we inferred a null demographic background utilizing SNP data, which was then exploited to simulate a putative null distribution of summary statistics that was compared with LINE data. Subsequently, we applied the inferred null demographic model with an additional exponential size change parameter, coupled with model selection, to test for neutrality as well as estimate the strength of either negative or positive selection. We found a robust signal for purifying selection in anole and mouse, but a lack of clear evidence for selection in stickleback and human. Overall, we demonstrated LINE insertion dynamics that are not in accordance to a mammalian versus nonmammalian dichotomy, and instead the degree of existing LINE activity together with host-specific demographic history may be the main determinants of LINE abundance.

International audience; Human-driven global environmental changes have considerably increased the risk of biological invasions, especially the spread of human parasites and their vectors. Among exotic species that have major impacts on public health, the dengue fever mosquito Aedes aegypti originating from Africa has spread worldwide during the last three centuries. Although considerable progress has been recently made in understanding the history of this invasion, the respective roles of human and abiotic factors in shaping patterns of genetic diversity remain largely unexplored. Using a genome-wide sample of genetic variants (3,530 ddRAD SNPs), we analyzed the genetic structure of Ae. aegypti populations in the Caribbean, the first introduced territories in the Americas. Fourteen populations were sampled in Guyane and in four islands of the Antilles that differ in climatic conditions, intensity of urbanization, and vector control history. The genetic diversity in the Caribbean was low (He = 0.14-0.17), as compared with a single African collection from Benin (He = 0.26) and site-frequency spectrum analysis detected an ancient bottleneck dating back ∼300 years ago, supporting a founder event during the introduction of Ae. aegypti. Evidence for a more recent bottleneck may be related to the eradication program undertaken on the American continent in the 1950s. Among 12 loci detected as FST-outliers, two were located in candidate genes for insecticide resistance (cytochrome P450 and voltage-gated sodium channel). Genome-environment association tests identified additional loci associated with human density and/or deltamethrin resistance. Our results highlight the high impact of human pressures on the demographic history and genetic variation of Ae. aegypti Caribbean populations.

Analyses of DNA from fossil specimens can generate critical insights into the demographic histories of natural populations. Previous analyses of fossil specimens from a single cave site in the Limay Valley in Patagonia have revealed striking historical differences in genetic variability between two species of Ctenomys, the colonial tuco-tuco (C. sociabilis) and the Patagonian tuco-tuco (C. haigi). As a first step towards identifying environmental or other factors contributing to this outcome, we assessed whether these differences in variability are generally characteristic of these species. We sequenced a 136-bp fragment of the cytochrome b locus for fossil specimens of Ctenomys excavated from two additional cave sites in the Limay Valley. Analyses of this expanded data set revealed that while C. sociabilis has undergone a pronounced loss of genetic diversity at all three cave sites over the past ~12 000 years, genetic diversity in C. haigi has remained relatively constant over the same temporal and geographic scales. The generality of these patterns suggests that although the factors affecting genetic diversity in C. sociabilis were widespread in the Limay Valley, impacts on the study species differed, probably due to known behavioural, ecological and demographic differences between these taxa. Fil: Tammone, Mauro Nicolás. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Centro Nacional Patagónico. Instituto de Diversidad y Evolución Austral; Argentina. Programa de Estudios Aplicados a la Conservación del Parque Nacional Nahuel Huapi; Argentina Fil: Pardiñas, Ulises Francisco J.. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Centro Nacional Patagónico. Instituto de Diversidad y Evolución Austral; Argentina Fil: Lacey, Eileen Anne. University of California at Berkeley; Estados Unidos

In some species, meiotic recombination is concentrated in small genomic regions. These “recombination hotspots” leave signatures in fine-scale patterns of linkage disequilibrium, raising the prospect that the genomic landscape of hotspots can be characterized from sequence variation. This approach has led to the inference that hotspots evolve rapidly in some species, but are conserved in others. Historic demographic events, such as population bottlenecks, are known to affect patterns of linkage disequilibrium across the genome, violating population genetic assumptions of this approach. Although such events are prevalent, demographic history is generally ignored when making inferences about the evolution of recombination hotspots. To determine the effect of demography on the detection of recombination hotspots, we use the coalescent to simulate haplotypes with a known recombination landscape. We measure the ability of popular linkage disequilibrium-based programs to detect hotspots across a range of demographic histories, including population bottlenecks, hidden population structure, population expansions, and population contractions. We find that demographic events have the potential to greatly reduce the power and increase the false positive rate of hotspot discovery. Neither the power nor the false positive rate of hotspot detection can be predicted without also knowing the demographic history of the sample. Our results suggest that ignoring demographic history likely overestimates the power to detect hotspots and therefore underestimates the degree of hotspot sharing between species. We suggest strategies for incorporating demographic history into population genetic inferences about recombination hotspots.

Identifying genetically different groups of animals, occupying specific geographical areas, is a prerequisite for conservation and management priorities. In the present study, the genetic structure of Atlantic walruses (Odobenus rosmarus rosmarus) occupying the Pechora Sea (PEC) in the western Russian Arctic, including walruses from Svalbard-Franz Josef Land (SVA-FJL) and East Greenland (EGR) regions, was investigated using 14 microsatellites (N = 159) and mtDNA sequences (N = 212). Bayesian-based clustering analysis identified two clusters: EGR and the other Northeast Atlantic areas. Pairwise F-ST analyses based on microsatellites revealed low but significant genetic differences between walruses from the PEC and SVA-FJL groups, which was supported by mtDNA analysis. F-ST was not significant for all sampling years, indicating a temporal effect or male-biased gene flow. Extended Bayesian Skyline Plots suggested a constant female subpopulation size (N-ef) for EGR and an increase for the SVA-FJL and PEC groups that commenced around 40-30 Kyr ago, indicating different demographic histories for walruses in the EGR. Further, the evolutionary phylogenetic relationship between Atlantic and Pacific walruses (O. r. divergence), based on mtDNA sequences, showed a monophyletic Atlantic clade, suggesting that Atlantic and Pacific walruses diverged similar to 949 Kyr. The principal finding suggests that PEC walruses show low, but significant genetic distinction from walruses in SVA-FJL and should be managed conservatively, as a separate, small population.

Siberia is one of the coldest environments on Earth and has great seasonal temperature variation. Long-term settlement in northern Siberia undoubtedly required biological adaptation to severe cold stress, dramatic variation in photoperiod, and limited food resources. In addition, recent archeological studies show that humans first occupied Siberia at least 45,000 years ago; yet our understanding of the demographic history of modern indigenous Siberians remains incomplete. In this study, we use whole-exome sequencing data from the Nganasans and Yakuts to infer the evolutionary history of these two indigenous Siberian populations. Recognizing the complexity of the adaptive process, we designed a model-based test to systematically search for signatures of polygenic selection. Our approach accounts for stochasticity in the demographic process and the hitchhiking effect of classic selective sweeps, as well as potential biases resulting from recombination rate and mutation rate heterogeneity. Our demographic inference shows that the Nganasans and Yakuts diverged ∼12,000-13,000 years ago from East-Asian ancestors in a process involving continuous gene flow. Our polygenic selection scan identifies seven candidate gene sets with Siberian-specific signals. Three of these gene sets are related to diet, especially to fat metabolism, consistent with the hypothesis of adaptation to a fat-rich animal diet. Additional testing rejects the effect of hitchhiking and favors a model in which selection yields small allele frequency changes at multiple unlinked genes.

Climate change during the late Quaternary period (LQP), in both temperate and tropical zones, has been a major driver in the shaping of species distributions and abundances. Our understanding of the various effects of climate change on population dynamics of marine species in temperate zones (such as the Gulf of California) is growing. However, studies on the demographic history of seabirds are rare and there is no description of how regional climate change has affected high-trophic-level marine species such as Heermann's Gull (Larus heermanni; Charadriiformes: Laridae). We investigated whether the demographic history of Heermann's Gull reflects population change consistent with past changes in climate in the Gulf of California during the LQP. We also explored whether those past changes affected the demographic history of codistributed marine organisms in a similar way as found for Heermann's Gulls. To test our hypotheses, we sequenced a fragment of the mitochondrial DNA cytochrome b gene (1,033 ...

The demographic history of most species is complex, with multiple evolutionary processes combining to shape the observed patterns of genetic diversity. To infer this history, the discipline of phylogeography has (to date) used models that simplify the historical demography of the focal organism, for example by assuming or ignoring ongoing gene flow between populations or by requiring a priori specification of divergence history. Since no single model incorporates every possible evolutionary process, researchers rely on intuition to choose the models that they use to analyze their data. Here, we describe an approximate likelihood approach that reduces this reliance on intuition. PHRAPL allows users to calculate the probability of a large number of complex demographic histories given a set of gene trees, enabling them to identify the most likely underlying model and estimate parameters for a given system. Available model parameters include coalescence time among populations or species, gene flow, and population size. We describe the method and test its performance in model selection and parameter estimation using simulated data. We also compare model probabilities estimated using our approximate likelihood method to those obtained using standard analytical likelihood. The method performs well under a wide range of scenarios, although this is sometimes contingent on sampling many loci. In most scenarios, as long as there are enough loci and if divergence among populations is sufficiently deep, PHRAPL can return the true model in nearly all simulated replicates. Parameter estimates from the method are also generally accurate in most cases. PHRAPL is a valuable new method for phylogeographic model selection and will be particularly useful as a tool to more extensively explore demographic model space than is typically done or to estimate parameters for complex models that are not readily implemented using current methods. Estimating relevant parameters using the most appropriate demographic model can help to sharpen our understanding of the evolutionary processes giving rise to phylogeographic patterns. [AIC; grid search; isolation-with-migration; migration rate; multispecies coalescent; parameter optimization; population genetics; tree topologies.].

New study systems and tools are needed to understand how divergence and speciation occur between lineages with gene flow. Migratory birds often exhibit divergence despite seasonal migration, which brings populations into contact with one another. We studied divergence between 2 subspecies of Northern Saw-whet Owl (Aegolius acadicus), in which a sedentary population on the islands of Haida Gwaii, British Columbia (A. a. brooksi), exists in the presence of the other form (A. a. acadicus) during migration but not during the breeding season. Prior research showed fixed mtDNA divergence but left open the question of nuclear gene flow. We used 2,517 ultraconserved element loci to examine the demographic history of this young taxon pair. Although we did not observe fixed single nucleotide polymorphism differences between populations among our genotyped individuals, 100% of the birds were diagnosable and δaδI analyses suggested the demographic model best fitting the data was one of split-bidirectional-migration (i.e. speciation with gene flow). We dated the split between brooksi and acadicus to ~278 Kya, and our analyses suggested gene flow between groups was skewed, with ~0.7 individuals per generation coming from acadicus into brooksi and ~4.4 going the opposite direction. Coupled with an absence of evidence of phenotypic hybrids and the birds’ natural history, these data suggest brooksi may be a young biological species arising despite historic gene flow.

Motivation Haplotypes from one or multiple related populations share a common genealogical history. If this shared genealogy can be inferred from haplotypes, it can be very useful for many population genetics problems. However, with the presence of recombination, the genealogical history of haplotypes is complex and cannot be represented by a single genealogical tree. Therefore, inference of genealogical history with recombination is much more challenging than the case of no recombination. Results In this paper, we present a new approach called RENT+ for the inference of local genealogical trees from haplotypes with the presence of recombination. RENT+ builds on a previous genealogy inference approach called RENT, which infers a set of related genealogical trees at different genomic positions. RENT+ represents a significant improvement over RENT in the sense that it is more effective in extracting information contained in the haplotype data about the underlying genealogy than RENT. The key components of RENT+ are several greatly enhanced genealogy inference rules. Through simulation, we show that RENT+ is more efficient and accurate than several existing genealogy inference methods. As an application, we apply RENT+ in the inference of population demographic history from haplotypes, which outperforms several existing methods. Availability and Implementation RENT+ is implemented in Java, and is freely available for download from: https://github.com/SajadMirzaei/RentPlus. Supplementary information Supplementary data are available at Bioinformatics online.

In this study, we evaluate phylogeographic patterns and predictions of ecological niche modelling (ENM) for Eugenia uniflora (Myrtaceae), a widely distributed taxon in the Atlantic forest domain, to understand the effect of past climatic oscillations on the demographic history of this species. An analysis of phylogeographic population structure and demography was conducted on E. uniflora from 46 localities in natural environments across the distribution range of the species based on three plastid markers. ENM was also performed to predict suitable environments and areas of dramatic decrease in future suitability for the species under distinct representative concentration pathways (RCPs). Eugenia uniflora exhibited higher haplotype and nucleotide diversity in the southern part of its distribution than in the northern part. Two divergent lineages were revealed in the phylogenetic analysis of haplotypes, with an estimated divergence at c. 4.9 Mya. The populations in the northern and central regions of the range probably experienced population growth, whereas populations in the southern region are marked by historical demographic stability. ENM results indicate that the distribution of E. uniflora was fragmented in cool periods and was broader and more connected during warm periods during Pleistocene. The results suggest distinct evolutionary histories in southern to northern populations, indicating region-specific responses to changes.

We demonstrate the advantages of using information at many unlinked loci to better calibrate estimates of the time to the most recent common ancestor (TMRCA) at a given locus. To this end, we apply a simple empirical Bayes method to estimate the TMRCA. This method is both asymptotically optimal, in the sense that the estimator converges to the true value when the number of unlinked loci for which we have information is large, and has the advantage of not making any assumptions about demographic history. The algorithm works as follows: we first split the sample at each locus into inferred left and right clades to obtain many estimates of the TMRCA, which we can average to obtain an initial estimate of the TMRCA. We then use nucleotide sequence data from other unlinked loci to form an empirical distribution that we can use to improve this initial estimate.

Population declines resulting from anthropogenic activities are of major consequence for the long-term survival of species because the resulting loss of genetic diversity can lead to extinction via the effects of inbreeding depression, fixation of deleterious mutations, and loss of adaptive potential. Otariid pinnipeds have been exploited commercially to near extinction with some species showing higher demographic resilience and recolonization potential than others. The New Zealand fur seal (NZFS) was heavily impacted by commercial sealing between the late 18th and early 19th centuries, but has recolonized its former range in southern Australia. The species has also recolonized its former range in New Zealand, yet little is known about the pattern of recolonization. Here, we first used 11 microsatellite markers (n = 383) to investigate the contemporary population structure and dispersal patterns in the NZFS (Arctocephalus forsteri). Secondly, we model postsealing recolonization with 1 additional mtDNA cytochrome b (n = 261) marker. Our data identified 3 genetic clusters: an Australian, a subantarctic, and a New Zealand one, with a weak and probably transient subdivision within the latter cluster. Demographic history scenarios supported a recolonization of the New Zealand coastline from remote west coast colonies, which is consistent with contemporary gene flow and with the species' high resilience. The present data suggest the management of distinct genetic units in the North and South of New Zealand along a genetic gradient. Assignment of individuals to their colony of origin was limited (32%) with the present data indicating the current microsatellite markers are unlikely sufficient to assign fisheries bycatch of NZFSs to colonies.

The term ‘ancient DNA’ (aDNA) is coming of age, with over 1,200 hits in the PubMed database, beginning in the early 1980s with the studies of ‘molecular paleontology’. Rooted in cloning and limited sequencing of DNA from ancient remains during the pre-PCR era, the field has made incredible progress since the introduction of PCR and next-generation sequencing. Over the last decade, aDNA analysis ushered in a new era in genomics and became the method of choice for reconstructing the history of organisms, their biogeography, and migration routes, with applications in evolutionary biology, population genetics, archaeogenetics, paleo-epidemiology, and many other areas. This change was brought by development of new strategies for coping with the challenges in studying aDNA due to damage and fragmentation, scarce samples, significant historical gaps, and limited applicability of population genetics methods. In this review, we describe the state-of-the-art achievements in aDNA studies, with particular focus on human evolution and demographic history. We present the current experimental and theoretical procedures for handling and analysing highly degraded aDNA. We also review the challenges in the rapidly growing field of ancient epigenomics. Advancement of aDNA tools and methods signifies a new era in population genetics and evolutionary medicine research.

The Andean degu, Octodontomys gliroides Gervais & d'Orbigny, 1844, has a broad distribution inhabiting pre-Andean pre-Puna and Puna environments of tropical South America. In order to understand the phylogeographic patterns of Octodontomys gliroides, we sequenced 579 bp of the mitochondrial DNA control region from 100 individuals collected from 20 populations across its entire distributional range. The phylogenetic and parsimony network, in conjunction with analysis of molecular variance (AMOVA), revealed a structured pattern of geographic differentiation of O. gliroides, with the occurrence of two well-defined evolutionary lineages: lineage A, restricted to Bolivia and Chile, and lineage B, restricted mainly to Argentina. Analysis of population structure inferred three genetic clusters along the distribution of O. gliroides that mostly agree with the four major barriers inferred by BARRIER analysis (e.g. rivers, salt flats, deserts, and mountain systems). In addition to the significant differentiation found among all levels studied, a positive correlation was identified between genetic and geographic distance, similar to as expected under the isolation-by-distance model. The most recent common ancestor of O. gliroides was estimated as c. 5.99 Mya, and the divergence between lineages A and B is estimated to have occurred by the Middle Pleistocene, about 0.69 Mya. The mismatch distributions and neutrality tests suggested a signal of population range expansion for both lineages coincident with major climatic changes that occurred during the wet–dry events of the Pleistocene in the Andean Puna region. Bayesian skyline plots (BSPs) for lineage A suggest a long history of constant population size followed by a period of slight to moderate demographic expansion at c. 0.04 Mya, whereas lineage B remained unclear after BSP analysis, probably because of the limited sample size.

Understanding the causes of pronounced losses of genetic diversity in natural populations may provide important insights into the evolutionary significance of these events. However, such analyses are typically based on post-reduction levels and patterns of variability in modern populations, which often lead to results that are biased towards more recent demographic events. In this context, population data recovered from the fossil record provide a powerful resource for studying historical processes of genetic decline. Using radiocarbon dating, stratigraphic analyses and DNA sequencing, we compared genetic variation and relative abundances of fossil specimens from the mid-Holocene to today to evaluate potential explanations for the marked historical loss of genetic diversity in the colonial tuco-tuco (Ctenomys sociabilis), a subterranean rodent that is endemic to Patagonia. We conclude that a general pattern of climatic change during the mid-Holocene, particularly changes in precipitation, led to changes in abundance of suitable habitats for this species. Loss of suitable habitat combined with the unusual demographic structure of this species may have facilitated the decline of populations of C. sociabilis by decreasing gene flow and increasing the potential for fixation of haplotypes due to genetic drift. Our analyses of temporal changes in abundance and genetic diversity in ctenomyids have implications for understanding more widespread patterns of Holocene change in the mammalian fauna of Patagonia. Fil: Tammone, Mauro Nicolás. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Centro Nacional Patagónico. Instituto de Diversidad y Evolución Austral; Argentina Fil: Pardiñas, Ulises Francisco J.. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Centro Nacional Patagónico. Instituto de Diversidad y Evolución Austral; Argentina Fil: Lacey, Eileen Anne. University of California at Berkeley; Estados Unidos

Using South African Families (SAF), a new database of settler genealogies, we provide the first comprehensive analysis of women’s cohort fertility in settler South Africa between 1700 and 1900. We find parity rates of approximately seven children per woman until a decline begins starting with women born in the 1850s. We date the start of South Africa’s fertility transition to cohorts born in the 1850s, having children from the 1870s. Both average parity and the timing of the transition are similar to other settler societies suggesting that although the sample suffers from selection it does not bias the parity estimates.

In the face of climate change, there is a growing need for research into the ability of organisms to persist at the limits of their bioclimatic envelope. American pikas ( Ochotona princeps ) have emerged as a focal mammalian species for investigating extinction risk related to climate change; however, most studies have occurred in characteristic alpine talus habitat within the range core. In the Columbia River Gorge (CRG), Oregon, American pikas inhabit low-elevation talus slopes previously considered outside the species’ bioclimatic range. We used microsatellite genotypic data to reconstruct levels of genetic variation, population connectivity, and demographic history at 11 CRG sites spanning an elevational gradient (104–1,292 m). Sampled sites separated into 2 genetic clusters largely explained by elevation, topography, and geographic proximity, with pairwise estimates of differentiation and migration rates suggesting little gene flow may be occurring. Sites were characterized by levels of allelic richness and heterozygosity substantially lower than values reported at characteristic alpine sites from the range core. Evidence of recent demographic contraction was found almost exclusively at high-elevation sites despite these areas being considered refuges from climate warming in more typical habitat. Given their unique genetic characteristics and persistence in an atypical environment, the CRG pika populations likely constitute a significant component of intraspecific biodiversity with high conservation value.