Population genetics from 1966 to 2016

scalepopgen: Bioinformatic Workflow Resources Implemented in Nextflow for Comprehensive Population Genomic Analyses

Population genomic analyses such as inference of population structure and identifying signatures of selection usually involve the application of a plethora of tools. The installation of tools and their dependencies, data transformation, or series of data preprocessing in a particular order sometimes makes the analyses challenging. While the usage of container-based technologies has significantly resolved the problems associated with the installation of tools and their dependencies, population genomic analyses requiring multistep pipelines or complex data transformation can greatly be facilitated by the application of workflow management systems such as Nextflow and Snakemake. Here, we present scalepopgen, a collection of fully automated workflows that can carry out widely used population genomic analyses on the biallelic single nucleotide polymorphism data stored in either variant calling format files or the plink-generated binary files. scalepopgen is developed in Nextflow and can be run locally or on high-performance computing systems using either Conda, Singularity, or Docker. The automated workflow includes procedures such as (i) filtering of individuals and genotypes; (ii) principal component analysis, admixture with identifying optimal K-values; (iii) running TreeMix analysis with or without bootstrapping and migration edges, followed by identification of an optimal number of migration edges; (iv) implementing single-population and pair-wise population comparison-based procedures to identify genomic signatures of selection. The pipeline uses various open-source tools; additionally, several Python and R scripts are also provided to collect and visualize the results. The tool is freely available at https://github.com/Popgen48/scalepopgen.

Speciation, natural selection, and networks: three historians versus theoretical population geneticists

In 1913, the geneticist William Bateson called for a halt in studies of genetic phenomena until evolutionary fundamentals had been sufficiently addressed at the molecular level. Nevertheless, in the 1960s, the theoretical population geneticists celebrated a "modern synthesis" of the teachings of Mendel and Darwin, with an exclusive role for natural selection in speciation. This was supported, albeit with minor reservations, by historians Mark Adams and William Provine, who taught it to generations of students. In subsequent decades, doubts were raised by molecular biologists and, despite the deep influence of various mentors, Adams and Provine noted serious anomalies and began to question traditional "just-so-stories." They were joined in challenging the genetic orthodoxy by a scientist-historian, Donald Forsdyke, who suggested that a "collective variation" postulated by Darwin's young research associate, George Romanes, and a mysterious "residue" postulated by Bateson, might relate to differences in short runs of DNA bases (oligonucleotides). The dispute between a small network of historians and a large network of geneticists can be understood in the context of national politics. Contrasts are drawn between democracies, where capturing the narrative makes reversal difficult, and dictatorships, where overthrow of a supportive dictator can result in rapid reversal.

Conservation genetics of Notelaea lloydii (Oleaceae) in south-eastern Queensland, Australia

Habitat fragmentation can increase the chance of population bottlenecks and inbreeding, and may ultimately lead to reduced fitness and local extinction. Notelaea lloydii is a native olive species endemic to Australia and listed as vulnerable due to its restricted distribution. A recent molecular systematics study has revealed there might be some geographic structuring among N. lloydii populations. Therefore, we undertook a genome-wide single nucleotide polymorphism (SNP) analysis to determine levels and patterns of genetic diversity, inbreeding and gene flow within and among N. lloydii populations in south-eastern Queensland. Furthermore, as the reproductive phase of a plant's life history has a profound influence on genetic diversity, life history reproductive traits were also studied. Our SNP analysis revealed low genetic diversity, inbreeding and significant genetic structuring even among proximate populations. Results of a flower and fruit bagging experiment in two consecutive seasons revealed that N. lloydii produced many flowers but only a few fruits survived to maturity. There were no differences in bagged and un-bagged flowering and fruiting rates, and therefore, we conclude that the high fruit abortion rate was probably due to inbreeding depression and/or suboptimal conditions, rather than pollinator availability and insect attack. Overall, results of this study indicate that the populations of N. lloydii are small, inbred and genetically isolated and represent unique management units that require local conservation management due to ongoing threats associated with urbanisation.

Harnessing deep learning for population genetic inference

In population genetics, the emergence of large-scale genomic data for various species and populations has provided new opportunities to understand the evolutionary forces that drive genetic diversity using statistical inference. However, the era of population genomics presents new challenges in analysing the massive amounts of genomes and variants. Deep learning has demonstrated state-of-the-art performance for numerous applications involving large-scale data. Recently, deep learning approaches have gained popularity in population genetics; facilitated by the advent of massive genomic data sets, powerful computational hardware and complex deep learning architectures, they have been used to identify population structure, infer demographic history and investigate natural selection. Here, we introduce common deep learning architectures and provide comprehensive guidelines for implementing deep learning models for population genetic inference. We also discuss current challenges and future directions for applying deep learning in population genetics, focusing on efficiency, robustness and interpretability.

Immunogenetics of autism spectrum disorder: A systematic literature review

The aetiology of autism spectrum disorder (ASD) is complex and, partly, accounted by genetic factors. Nonetheless, the genetic underpinnings of ASD are poorly defined. The presence of immune dysregulations in autistic individuals, and their families, supports a role of the immune system and its genetic regulators. Albeit immune responses belong either to the innate or adaptive arms, the overall immune system genetics is broad, and encompasses a multitude of functionally heterogenous pathways which may have different influences on ASD. Hence, to gain insights on the immunogenetic underpinnings of ASD, we conducted a systematic literature review of previous immune genetic and transcription studies in ASD. We defined a list of immune genes relevant to ASD and explored their neuro-immune function. Our review confirms the presence of immunogenetic variability in ASD, accounted by inherited variations of innate and adaptive immune system genes and genetic expression changes in the blood and post-mortem brain of autistic individuals. Besides their immune function, the identified genes control neurodevelopment processes (neuronal and synaptic plasticity) and are highly expressed in pre/peri-natal periods. Hence, our synthesis bolsters the hypothesis that perturbation in immune genes may contribute to ASD by derailing the typical trajectory of neurodevelopment. Our review also helped identifying some of the limitations of prior immunogenetic research in ASD. Thus, alongside clarifying the neurodevelopment role of immune genes, we outline key considerations for future work into the aetiology of ASD and possible novel intervention targets.

Genetic footprint of population diversity and genetic structure of Venturia inaequalis infecting apple (Malus × domestica Borkh.)

Apple scab instigated by Venturia inaequalis impels remarkable losses to apple fruit production. In an effort to comprehend the key mechanisms of evolutionary potential defining V. inaequalis population, 132 isolates of V. inaequalis from five commercial apple orchards were collected and assayed using 14 microsatellite markers. The average diversity was observed within the individuals of populations based on the Shannon-Wieners index (I) and observed heterozygosity (Ho) was average but considerably lower than expected heterozygosity (He). The genetic differentiation based on FST values was revealed as an average measure of divergence between populations and had varying proportions of gene flow and migration among themselves. Analysis of Molecular Variance (AMOVA) revealed that variance (94%) was dispersed across individuals with a significant (6%) variation between populations from different regions. To examine host specialization within the V. inaequalis population, the assignment approach based on K-means of clustering (an unsupervised machine learning approach), revealed that the clustering method supported three clusters at (K = 3) and three major clusters were also observed in Principle Component Analysis (PCA). Additionally, Nei's genetic distance values, pairwise estimates of genetic differentiation, dendrogram using the neighbor-joining and PCoA revealed the random distribution of V. Inaequalis isolates that depicted a high proportion of genotypic diversity within populations and population genetic structure.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03672-2.

Genomic diversity landscapes in outcrossing and selfing Caenorhabditis nematodes

Caenorhabditis nematodes form an excellent model for studying how the mode of reproduction affects genetic diversity, as some species reproduce via outcrossing whereas others can self-fertilize. Currently, chromosome-level patterns of diversity and recombination are only available for self-reproducing Caenorhabditis, making the generality of genomic patterns across the genus unclear given the profound potential influence of reproductive mode. Here we present a whole-genome diversity landscape, coupled with a new genetic map, for the outcrossing nematode C. remanei. We demonstrate that the genomic distribution of recombination in C. remanei, like the model nematode C. elegans, shows high recombination rates on chromosome arms and low rates toward the central regions. Patterns of genetic variation across the genome are also similar between these species, but differ dramatically in scale, being tenfold greater for C. remanei. Historical reconstructions of variation in effective population size over the past million generations echo this difference in polymorphism. Evolutionary simulations demonstrate how selection, recombination, mutation, and selfing shape variation along the genome, and that multiple drivers can produce patterns similar to those observed in natural populations. The results illustrate how genome organization and selection play a crucial role in shaping the genomic pattern of diversity whereas demographic processes scale the level of diversity across the genome as a whole.

Inferring genetic structure when there is little: population genetics versus genomics of the threatened bat Miniopterus schreibersii across Europe

Despite their paramount importance in molecular ecology and conservation, genetic diversity and structure remain challenging to quantify with traditional genotyping methods. Next-generation sequencing holds great promises, but this has not been properly tested in highly mobile species. In this article, we compared microsatellite and RAD-sequencing (RAD-seq) analyses to investigate population structure in the declining bent-winged bat (Miniopterus schreibersii) across Europe. Both markers retrieved general patterns of weak range-wide differentiation, little sex-biased dispersal, and strong isolation by distance that associated with significant genetic structure between the three Mediterranean Peninsulas, which could have acted as glacial refugia. Microsatellites proved uninformative in individual-based analyses, but the resolution offered by genomic SNPs illuminated on regional substructures within several countries, with colonies sharing migrators of distinct ancestry without admixture. This finding is consistent with a marked philopatry and spatial partitioning between mating and rearing grounds in the species, which was suspected from marked-recaptured data. Our study advocates that genomic data are necessary to properly unveil the genetic footprints left by biogeographic processes and social organization in long-distant flyers, which are otherwise rapidly blurred by their high levels of gene flow.

Balancing selection on genomic deletion polymorphisms in humans

A key question in biology is why genomic variation persists in a population for extended periods. Recent studies have identified examples of genomic deletions that have remained polymorphic in the human lineage for hundreds of millennia, ostensibly owing to balancing selection. Nevertheless, genome-wide investigation of ancient and possibly adaptive deletions remains an imperative exercise. Here, we demonstrate an excess of polymorphisms in present-day humans that predate the modern human-Neanderthal split (ancient polymorphisms), which cannot be explained solely by selectively neutral scenarios. We analyze the adaptive mechanisms that underlie this excess in deletion polymorphisms. Using a previously published measure of balancing selection, we show that this excess of ancient deletions is largely owing to balancing selection. Based on the absence of signatures of overdominance, we conclude that it is a rare mode of balancing selection among ancient deletions. Instead, more complex scenarios involving spatially and temporally variable selective pressures are likely more common mechanisms. Our results suggest that balancing selection resulted in ancient deletions harboring disproportionately more exonic variants with GWAS (genome-wide association studies) associations. We further found that ancient deletions are significantly enriched for traits related to metabolism and immunity. As a by-product of our analysis, we show that deletions are, on average, more deleterious than single nucleotide variants. We can now argue that not only is a vast majority of common variants shared among human populations, but a considerable portion of biologically relevant variants has been segregating among our ancestors for hundreds of thousands, if not millions, of years.

Distinct signals of clinal and seasonal allele frequency change at eQTLs in Drosophila melanogaster

Populations of short-lived organisms can respond to spatial and temporal environmental heterogeneity through local adaptation. Local adaptation can be reflected on both phenotypic and genetic levels, and it has been documented in many organisms. Although complex fitness-related phenotypes have been shown to vary across latitudinal clines and seasons in similar ways in Drosophila melanogaster populations, the comparative signals of local adaptation across space and time remain poorly understood. Here, we examined patterns of allele frequency change across a latitudinal cline and between seasons at previously reported expression quantitative trait loci (eQTLs). We divided eQTLs into groups by using differential expression profiles of fly populations collected across latitudinal clines or exposed to different environmental conditions. In general, we find that eQTLs are enriched for clinally varying polymorphisms, and that these eQTLs change in frequency in concordant ways across the cline and in response to starvation and chill-coma. The enrichment of eQTLs among seasonally varying polymorphisms is more subtle, and the direction of allele frequency change at eQTLs appears to be somewhat idiosyncratic. Taken together, we suggest that clinal adaptation at eQTLs is at least partially distinct from seasonal adaptation.

Epigenetic variation: A major player in facilitating plant fitness under changing environmental conditions

Recent research in plant epigenetics has increased our understanding of how epigenetic variability can contribute to adaptive phenotypic plasticity in natural populations. Studies show that environmental changes induce epigenetic switches either independently or in complementation with the genetic variation. Although most of the induced epigenetic variability gets reset between generations and is short-lived, some variation becomes transgenerational and results in heritable phenotypic traits. The short-term epigenetic responses provide the first tier of transient plasticity required for local adaptations while transgenerational epigenetic changes contribute to stress memory and help the plants respond better to recurring or long-term stresses. These transgenerational epigenetic variations translate into an additional tier of diversity which results in stable epialleles. In recent years, studies have been conducted on epigenetic variation in natural populations related to various biological processes, ecological factors, communities, and habitats. With the advent of advanced NGS-based technologies, epigenetic studies targeting plants in diverse environments have increased manifold to enhance our understanding of epigenetic responses to environmental stimuli in facilitating plant fitness. Taking all points together in a frame, the present review is a compilation of present-day knowledge and understanding of the role of epigenetics and its fitness benefits in diverse ecological systems in natural populations.

Population genomics reveals moderate genetic differentiation between populations of endangered Forest Musk Deer located in Shaanxi and Sichuan

BACKGROUND

Many endangered species exist in small, genetically depauperate, or inbred populations, hence promoting genetic differentiation and reducing long-term population viability. Forest Musk Deer (Moschus berezovskii) has been subject to illegal hunting for hundreds of years due to the medical and commercial values of musk, resulting in a significant decline in population size. However, it is still unclear to what extent the genetic exchange and inbreeding levels are between geographically isolated populations. By using whole-genome data, we reconstructed the demographic history, evaluated genetic diversity, and characterized the population genetic structure of Forest Musk Deer from one wild population in Sichuan Province and two captive populations from two ex-situ centers in Shaanxi Province.

RESULTS

SNP calling by GATK resulted in a total of 44,008,662 SNPs. Principal component analysis (PCA), phylogenetic tree (NJ tree), ancestral component analysis (ADMIXTURE) and the ABBA-BABA test separated Sichuan and Shaanxi Forest Musk Deer as two genetic clusters, but no obvious genetic differentiation was observed between the two captive populations. The average pairwise FST value between the populations in Sichuan and Shaanxi ranged from 0.05-0.07, suggesting a low to moderate genetic differentiation. The mean heterozygous SNPs rate was 0.14% (0.11%-0.15%) for Forest Musk Deer at the genomic scale, and varied significantly among three populations (Chi-square = 1.22, p < 0.05, Kruskal-Wallis Test), with the Sichuan population having the lowest (0.11%). The nucleotide diversity of three populations varied significantly (p < 0.05, Kruskal-Wallis Test), with the Sichuan population having the lowest genetic θπ (1.69 × 10-3).

CONCLUSIONS

Genetic diversity of Forest Musk Deer was moderate at the genomic scale compared with other endangered species. Genetic differentiation between populations in Sichuan and Shaanxi may not only result from historical biogeographical factors but also be associated with contemporary human disturbances. Our findings provide scientific aid for the conservation and management of Forest Musk Deer. They can extend the proposed measures at the genomic level to apply to other musk deer species worldwide.

How methodological changes have influenced our understanding of population structure in threatened species: insights from tiger populations across India

Unprecedented advances in sequencing technology in the past decade allow a better understanding of genetic variation and its partitioning in natural populations. Such inference is critical to conservation: to understand species biology and identify isolated populations. We review empirical population genetics studies of Endangered Bengal tigers within India, where 60-70% of wild tigers live. We assess how changes in marker type and sampling strategy have impacted inferences by reviewing past studies, and presenting three novel analyses including a single-nucleotide polymorphism (SNP) panel, genome-wide SNP markers, and a whole-mitochondrial genome network. At a broad spatial scale, less than 100 SNPs revealed the same patterns of population clustering as whole genomes (with the exception of one additional population sampled only in the SNP panel). Mitochondrial DNA indicates a strong structure between the northeast and other regions. Two studies with more populations sampled revealed further substructure within Central India. Overall, the comparison of studies with varied marker types and sample sets allows more rigorous inference of population structure. Yet sampling of some populations is limited across all studies, and these should be the focus of future sampling efforts. We discuss challenges in our understanding of population structure, and how to further address relevant questions in conservation genetics. This article is part of the theme issue 'Celebrating 50 years since Lewontin's apportionment of human diversity'.

The apportionment of citations: a scientometric analysis of Lewontin 1972

'The apportionment of human diversity' (1972) is the most highly cited research article published by geneticist Richard Lewontin in his career. This study's primary result-that most genetic diversity in humans can be accounted for by within-population differences, not between-population differences-along with Lewontin's outspoken, politically charged interpretations thereof, has become foundational to the scientific and cultural discourse pertaining to human genetic variation. The article has an unusual bibliometric trajectory in that it is much more salient in the bibliographic record today compared to the first 20 years after its publication. Here, we highlight four factors that may have played a role in shaping the paper's fame: (i) citations in influential publications across several disciplines; (ii) Lewontin's own popular books and media appearances; (iii) the renaissance of population genetics research of the early 1990s; and (iv) the serendipitous collision of scientific progress, influential books and papers, and heated controversies around the year 1994. We conclude with an analysis of Twitter data to characterize the communities and conversations that continue to keep this study at the centre of discussions about race and genetics, prompting new challenges for scientists who have inherited Lewontin's legacy. This article is part of the theme issue 'Celebrating 50 years since Lewontin's apportionment of human diversity'.

Whole-Genome Genotyping Using DNA Microarrays for Population Genetics

The field of population genetics has exploded in the last two decades following the sequencing of the human genome in 2001 (Green et al. Nature 526:29-31, 2015). Tools to measure genetic variation have matured significantly throughout this advancement in knowledge (Lenoir and Giannella. J Biomed Discov Collab 1:11, 2006; Marzancola et al. Methods Mol Biol 1368:161-178, 2016). In this chapter, the focus is on the laboratory methods developed to perform genome-wide genotyping utilizing DNA microarrays, which is one of the most commonly used molecular techniques to assess global genetic variation (Heller MJ, Annu Rev Biomed Eng 4:129-153, 2002). DNA microarrays allow for the interrogation of hundreds of thousands of SNPs (single nucleotide polymorphisms) at once utilizing array-based technology in conjunction with fluorescent molecular labels in a process referred to as genotyping (Marzancola et al. Methods Mol Biol 1368:161-178, 2016). Genotype data can be utilized to associate certain phenotypes in relation with specific genetic variants within a population in a process known as genome-wide association studies or GWAS (Charlesworth and Charlesworth. Heredity (Edinb) 118(1):2-9, 2017; Casillas and Barbadilla. Genetics 205(3):1003-1035, 2017). This experimental technique is a multiple-day process involving the combination of DNA extraction, amplification, fragmentation, binding, and staining (Illumina Infinium HTS Assay Protocol Guide, 2013). Many vendors supply platforms and products to assess global genetic variation using DNA microarrays (Illumina Infinium HTS Assay Protocol Guide, 2013). In this chapter, the focus is on the methods utilized to generate high-quality genotype data with the Illumina^® Infinium Global Screening Array. Although data analysis and quality control are not the focus for this chapter, they are also briefly addressed.

Drivers of neutral and adaptive differentiation in pike (Esox lucius) populations from contrasting environments

Understanding how eco-evolutionary processes and environmental factors drive population differentiation and adaptation are key challenges in evolutionary biology of relevance for biodiversity protection. Differentiation requires at least partial reproductive separation, which may result from different modes of isolation such as geographic isolation (allopatry) or isolation by distance (IBD), resistance (IBR), and environment (IBE). Despite that multiple modes might jointly influence differentiation, studies that compare the relative contributions are scarce. Using RADseq, we analyse neutral and adaptive genetic diversity and structure in 11 pike (Esox lucius) populations from contrasting environments along a latitudinal gradient (54.9-63.6°N), to investigate the relative effects of IBD, IBE and IBR, and to assess whether the effects differ between neutral and adaptive variation, or across structural levels. Patterns of neutral and adaptive variation differed, probably reflecting that they have been differently affected by stochastic and deterministic processes. The importance of the different modes of isolation differed between neutral and adaptive diversity, yet were consistent across structural levels. Neutral variation was influenced by interactions among all three modes of isolation, with IBR (seascape features) playing a central role, wheares adaptive variation was mainly influenced by IBE (environmental conditions). Taken together, this and previous studies suggest that it is common that multiple modes of isolation interactively shape patterns of genetic variation, and that their relative contributions differ among systems. To enable identification of general patterns and understand how various factors influence the relative contributions, it is important that several modes are simultaneously investigated in additional populations, species and environmental settings.

Comparisons of historical Dutch commons inform about the long-term dynamics of social-ecological systems

Human societies and natural ecosystems are under threat by growing populations, overexploitation of natural resources and climate change. This calls for more sustainable utilization of resources based on past experiences and insights from many different disciplines. Interdisciplinary approaches to studies of historical commons have potential to identify drivers of change and keys to success in the past, and offer advice about the management and use of shared resources in contemporary and future systems. We address these issues by applying an ecological perspective to historical data on social-ecological systems. We perform comparisons and time series analyses for nine successful Dutch commons for which high-resolution data on the regulatory activities and use of shared resources is available for on average 380 years (range 236 to 568) during the period 1300 to 1972. Within commons, institutional developments were oscillating, with periods of intense regulatory activity being separated by periods of low activity, and with the dynamics of regulations being largely independent across commons. Ecological theory posits that species that occupy similar niches should show correlated responses to environmental challenges; however, commons using more similar resources did not have more parallel or similar institutional developments. One notable exception was that sanctioning was more frequent in commons that directed more regulatory activities towards non-renewable subsoil resources, whereas there was no association between sanctioning and the use of renewable resources. This might indicate that commoners were aware of potential resource depletion and attempted to influence freeriding by actively trying to solve the underlying social dilemmas. Sanctioning regulations were more frequent during the first than during the second part of a common’s life, indicating that while sanctioning might have been important for the establishment of commons it was not key to the long-term persistence of historical commons.

Richard Lewontin and the "complications of linkage"

During the 1960s and 1970s population geneticists pushed beyond models of single genes to grapple with the effect on evolution of multiple genes associated by linkage. The resulting models of multiple interacting loci suggested that blocks of genes, maybe even entire chromosomes or the genome itself, should be treated as a unit. In this context, Richard Lewontin wrote his famous 1974 book The Genetic Basis of Evolutionary Change, which concludes with an argument for considering the entire genome as the unit of selection as a result of linkage. Why did Lewontin and others devote so much intellectual energy to the "complications of linkage" in the 1960s and 1970s? We argue that this attention to linkage should be understood in the context of research on chromosomal inversions and co-adapted gene complexes that occupied mid-century evolutionary genetics. For Lewontin, the complications of linkage were an extension of this chromosomal focus expressed in the new language of models for linkage disequilibrium.

Natural selection and the distribution of chromosomal inversion lengths

Chromosomal inversions contribute substantially to genome evolution, yet the processes governing their evolutionary dynamics remain poorly understood. Theory suggests that a readily measurable property of inversions-their length-can potentially affect their evolutionary fates. Emerging data on the lengths of polymorphic and fixed inversions may therefore provide clues to the evolutionary processes promoting inversion establishment. However, formal predictions for the distribution of inversion lengths remain incomplete, making empirical patterns difficult to interpret. We model the relation between inversion length and establishment probability for four inversion types: (1) neutral, (2) underdominant, (3) directly beneficial, and (4) indirectly beneficial, with selection favouring the latter because they capture locally adapted alleles at migration-selection balance and suppress recombination between them. We also consider how deleterious mutations affect the lengths of established inversions. We show that length distributions of common polymorphic and fixed inversions systematically differ among inversion types. Small rearrangements contribute the most to genome evolution under neutral and underdominant scenarios of selection, with the lengths of neutral inversion substitutions increasing, and those of underdominant substitutions decreasing, with effective population size. Among directly beneficial inversions, small rearrangements are preferentially fixed, whereas intermediate-to-large inversions are maintained as balanced polymorphisms via associative overdominance. Finally, inversions established under the local adaptation scenario are predominantly intermediate-to-large. Such inversions remain polymorphic or approach fixation within the local populations where they are favoured. Our models clarify how inversion length distributions relate to processes of inversion establishment, providing a platform for testing how natural selection shapes the evolution of genome structure.

How well do genetic markers inform about responses to intraspecific admixture? A comparative analysis of microsatellites and RADseq

Background

Fitness consequences of intraspecific genetic admixture can vary from positive to negative depending on the genetic composition of the populations and environmental conditions. Because admixture has potential to influence the success of management and conservation efforts, genetic similarity has been suggested to be used as a proxy to predict the outcome. Studies utilizing microsatellites (a neutral marker) to investigate associations between genetic distance and admixture effects show conflicting results. Marker types that yield information on genome-wide and/or adaptive variation might be more useful for predicting responses to inter-population hybridization. In this study we utilized published data for three populations of pike (Esox lucius) to investigate associations between offspring performance (hatching success) and parental genetic similarity in experimentally purebred and admixed families, based on neutral (microsatellites), genome-wide neutral (RADseq SNPs), and adaptive (SNPs under selection) markers.

Results

Estimated similarity varied among the markers, likely reflecting differences in their inherent properties, but was consistently higher in purebred than admixed families. A significant interaction between marker type and admixture treatment reflected that neutral SNPs yielded higher estimates than adaptive SNPs for admixed families whereas no difference was found for purebred families, which indicates that neutral similarity was not reflective of adaptive similarity. When all samples were pooled, no association between similarity and performance was found for any marker. For microsatellites, similarity was positively correlated with hatching success in purebred families, whereas no association was found in admixed families; however, the direction of the effect differed between the population combinations.

Conclusions

The results strengthen the notion that, as of today, there is no proxy that can reliably predicted the outcome of admixture. This emphasizes the need of further studies to advance knowledge that can shed light on how to safeguard against negative consequences of admixture, and thereby inform management and promote conservation of biological diversity.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-021-00974-3.

Unboxing mutations: Connecting mutation types with evolutionary consequences

A key step in understanding the genetic basis of different evolutionary outcomes (e.g., adaptation) is to determine the roles played by different mutation types (e.g., SNPs, translocations and inversions). To do this we must simultaneously consider different mutation types in an evolutionary framework. Here, we propose a research framework that directly utilizes the most important characteristics of mutations, their population genetic effects, to determine their relative evolutionary significance in a given scenario. We review known population genetic effects of different mutation types and show how these may be connected to different evolutionary outcomes. We provide examples of how to implement this framework and pinpoint areas where more data, theory and synthesis are needed. Linking experimental and theoretical approaches to examine different mutation types simultaneously is a critical step towards understanding their evolutionary significance.

The Pop-Gen Pipeline Platform: A Software Platform for Population Genomic Analyses

The Pop-Gen Pipeline Platform (PPP) is a software platform for population genomic analyses. The PPP was designed as a collection of scripts that facilitate common population genomic workflows in a consistent and standardized Python environment. Functions were developed to encompass entire workflows, including input preparation, file format conversion, various population genomic analyses, and output generation. The platform has also been developed with reproducibility and extensibility of analyses in mind. The PPP is an open-source package that is available for download and use at https://ppp.readthedocs.io/en/latest/PPP_pages/install.html.

Chromosome-Scale Genome Assemblies of Aphids Reveal Extensively Rearranged Autosomes and Long-Term Conservation of the X Chromosome

Chromosome rearrangements are arguably the most dramatic type of mutations, often leading to rapid evolution and speciation. However, chromosome dynamics have only been studied at the sequence level in a small number of model systems. In insects, Diptera and Lepidoptera have conserved genome structure at the scale of whole chromosomes or chromosome arms. Whether this reflects the diversity of insect genome evolution is questionable given that many species exhibit rapid karyotype evolution. Here, we investigate chromosome evolution in aphids-an important group of hemipteran plant pests-using newly generated chromosome-scale genome assemblies of the green peach aphid (Myzus persicae) and the pea aphid (Acyrthosiphon pisum), and a previously published assembly of the corn-leaf aphid (Rhopalosiphum maidis). We find that aphid autosomes have undergone dramatic reorganization over the last 30 My, to the extent that chromosome homology cannot be determined between aphids from the tribes Macrosiphini (Myzus persicae and Acyrthosiphon pisum) and Aphidini (Rhopalosiphum maidis). In contrast, gene content of the aphid sex (X) chromosome remained unchanged despite rapid sequence evolution, low gene expression, and high transposable element load. To test whether rapid evolution of genome structure is a hallmark of Hemiptera, we compared our aphid assemblies with chromosome-scale assemblies of two blood-feeding Hemiptera (Rhodnius prolixus and Triatoma rubrofasciata). Despite being more diverged, the blood-feeding hemipterans have conserved synteny. The exceptional rate of structural evolution of aphid autosomes renders them an important emerging model system for studying the role of large-scale genome rearrangements in evolution.

Population Genomics on the Fly: Recent Advances in Drosophila

Drosophila melanogaster, a small dipteran of African origin, represents one of the best-studied model organisms. Early work in this system has uniquely shed light on the basic principles of genetics and resulted in a versatile collection of genetic tools that allow to uncover mechanistic links between genotype and phenotype. Moreover, given its worldwide distribution in diverse habitats and its moderate genome-size, Drosophila has proven very powerful for population genetics inference and was one of the first eukaryotes whose genome was fully sequenced. In this book chapter, we provide a brief historical overview of research in Drosophila and then focus on recent advances during the genomic era. After describing different types and sources of genomic data, we discuss mechanisms of neutral evolution including the demographic history of Drosophila and the effects of recombination and biased gene conversion. Then, we review recent advances in detecting genome-wide signals of selection, such as soft and hard selective sweeps. We further provide a brief introduction to background selection, selection of noncoding DNA and codon usage and focus on the role of structural variants, such as transposable elements and chromosomal inversions, during the adaptive process. Finally, we discuss how genomic data helps to dissect neutral and adaptive evolutionary mechanisms that shape genetic and phenotypic variation in natural populations along environmental gradients. In summary, this book chapter serves as a starting point to Drosophila population genomics and provides an introduction to the system and an overview to data sources, important population genetic concepts and recent advances in the field.

Sympatric divergence of the ergot fungus, Claviceps purpurea, populations infecting agricultural and nonagricultural grasses in North America

The ergot diseases of agricultural and nonagricultural grasses are caused by the infection of Claviceps spp. (Hypocreales, Ascomycota) on florets, producing dark spur-like sclerotia on spikes that are toxic to humans and animals, leading to detrimental impacts on agriculture and economy due to the downgrading of cereal grains, import-export barriers, reduced yield, and ecological concerns. At least seven phylogenetic lineages (phylogenetic species) were identified within the premolecular concept of C. purpurea s.l. (sensu lato) in agricultural areas and vicinities in Canada and the Western United States. Claviceps purpurea s.s (sensu stricto) remained as the most prevalent species with a wide host range, including cereal crops, native, invasive, and weedy grasses. The knowledge on genetic diversity and distribution of C. purpurea s.s. in North America is lacking. The objective of the present study was to shed light on genetic differentiation and evolution of the natural populations of C. purpurea s.s. Multilocus DNA sequences of samples from Canada and the Western USA were analyzed using a phylogenetic network approach, and population demographic parameters were investigated. Results showed that three distinct genetically subdivided populations exist, and the subdivision is not correlated with geographic or host differentiations. Potential intrinsic mechanisms that might play roles in leading to the cessation of gene flows among the subpopulations, that is, mating and/or vegetative incompatibility, genomic adaptation, were discussed. The neutrality of two house-keeping genes that are widely used for DNA barcoding, that is, translation elongation factor 1-α (TEF1-α) and RNA polymerase II second largest subunit (RPB2), was challenged and discussed.

Population Genomics of Filamentous Plant Pathogens-A Brief Overview of Research Questions, Approaches, and Pitfalls

With ever-decreasing sequencing costs, research on the population biology of plant pathogens is transitioning from population genetics-using dozens of genetic markers or polymorphism data of several genes-to population genomics-using several hundred to tens of thousands of markers or whole-genome sequence data. The field of population genomics is characterized by rapid theoretical and methodological advances and by numerous steps and pitfalls in its technical and analytical workflow. In this article, we aim to provide a brief overview of topics relevant to the study of population genomics of filamentous plant pathogens and direct readers to more extensive reviews for in-depth understanding. We briefly discuss different types of population genomics-inspired research questions and give insights into the sampling strategies that can be used to answer such questions. We then consider different sequencing strategies, the various options available for data processing, and some of the currently available tools for population genomic data analysis. We conclude by highlighting some of the hurdles along the population genomic workflow, providing cautionary warnings relative to assumptions and technical challenges, and presenting our own future perspectives of the field of population genomics for filamentous plant pathogens.

On the History and Applications of Congenic Strains in Cryptococcus Research

Congenic strains have been utilized in numerous model organisms to determine the genetic underpinning of various phenotypic traits. Congenic strains are usually derived after 10 backcrosses to a recipient parent, at which point they are 99.95% genetically identical to the parental strain. In recent decades, congenic pairs have provided an invaluable tool for genetics and molecular biology research in the Cryptococcus neoformans species complex. Here, we summarize the history of Cryptococcus congenic pairs and their application in Cryptococcus research on topics including the impact of the mating type locus on unisexual reproduction, virulence, tissue tropism, uniparental mitochondrial inheritance, and the genetic underpinning of other various traits. We also discuss the limitations of these approaches and other biological questions, which could be explored by employing congenic pairs.

Variation in reproductive effort, genetic diversity and mating systems across Posidonia australis seagrass meadows in Western Australia

Populations at the edges of their geographical range tend to have lower genetic diversity, smaller effective population sizes and limited connectivity relative to centre of range populations. Range edge populations are also likely to be better adapted to more extreme conditions for future survival and resilience in warming environments. However, they may also be most at risk of extinction from changing climate. We compare reproductive and genetic data of the temperate seagrass, Posidonia australis on the west coast of Australia. Measures of reproductive effort (flowering and fruit production and seed to ovule ratios) and estimates of genetic diversity and mating patterns (nuclear microsatellite DNA loci) were used to assess sexual reproduction in northern range edge (low latitude, elevated salinities, Shark Bay World Heritage Site) and centre of range (mid-latitude, oceanic salinity, Perth metropolitan waters) meadows in Western Australia. Flower and fruit production were highly variable among meadows and there was no significant relationship between seed to ovule ratio and clonal diversity. However, Shark Bay meadows were two orders of magnitude less fecund than those in Perth metropolitan waters. Shark Bay meadows were characterized by significantly lower levels of genetic diversity and a mixed mating system relative to meadows in Perth metropolitan waters, which had high genetic diversity and a completely outcrossed mating system. The combination of reproductive and genetic data showed overall lower sexual productivity in Shark Bay meadows relative to Perth metropolitan waters. The mixed mating system is likely driven by a combination of local environmental conditions and pollen limitation. These results indicate that seagrass restoration in Shark Bay may benefit from sourcing plant material from multiple reproductive meadows to increase outcrossed pollen availability and seed production for natural recruitment.

Population genetics: past, present, and future

We present selected topics of population genetics and molecular phylogeny. As several excellent review articles have been published and generally focus on European and American scientists, here, we emphasize contributions by Japanese researchers. Our review may also be seen as a belated 50-year celebration of Motoo Kimura's early seminal paper on the molecular clock, published in 1968.

Tools for Evolutionary and Genetic Analysis (TEGA): A new platform for the management of molecular and environmental data

Population genetics studies the distributions and changes in population allele frequencies in response to processes, such as mutation, natural selection, gene flow, and genetic drift. Researchers daily manage genetic, biological, and environmental data of the samples, storing them in text files or spreadsheets, which makes it difficult to maintain consistency and traceability. Here we present TEGA, a WEB-based stand-alone software developed for the easy analysis and management of population genetics data. It was designed to: 1) facilitate data management, 2) provide a way to execute the analysis procedures, and 3) supply a means to publish data, procedures, and results. TEGA is distributed under the GNU AGPL v3 license. The documentation, source code, and screenshots are available at https://github.com/darioelias/TEGA. In addition, we present Rabid Fish, the first implementation of TEGA in the Genetics Labortory of the Faculty of Humanities and Sciences at the National University of the Litoral, where research focuses on population genetics studies applied to non-model organisms.

Comparing the Performance of Microsatellites and RADseq in Population Genetic Studies: Analysis of Data for Pike (Esox lucius) and a Synthesis of Previous Studies

Population genetic studies reveal biodiversity patterns and inform about drivers of evolutionary differentiation and adaptation, including gene flow, drift and selection. This can advance our understanding and aid decision making regarding management and conservation efforts. Microsatellites have long been used in population genetic studies. Thanks to the development of newer techniques, sequencing approaches such as restriction site associated DNA sequencing (RADseq) are on their way to replace microsatellites for some applications. However, the performance of these two marker types in population genetics have rarely been systematically compared. We utilized three neutrally and adaptively differentiated populations of anadromous pike (Esox lucius) to assess the relative performance of microsatellites and RADseq with respect to resolution and conclusiveness of estimates of population differentiation and genetic structure. To this end, the same set of individuals (N = 64) were genotyped with both RADseq and microsatellite markers. To assess effects of sample size, the same subset of 10 randomly chosen individuals from each population (N = 30 in total) were also genotyped with both methods. Comparisons of estimated genetic diversity and structure showed that both markers were able to uncover genetic structuring. The full RADseq dataset provided the clearest detection of the finer scaled genetic structuring, and the other three datasets (full and subset microsatellite, and subset RADseq) provided comparable results. A search for outlier loci performed on the full SNP dataset pointed to signs of selection potentially associated with salinity and temperature, exemplifying the utility of RADseq to inform about the importance of different environmental factors. To evaluate whether performance differences between the markers are general or context specific, the results of previous studies that have investigated population structure using both marker types were synthesized. The synthesis revealed that RADseq performed as well as, or better than microsatellites in detecting genetic structuring in the included studies. The differences in the ability to detect population structure, both in the present and the previous studies, are likely explained by the higher number of loci typically utilized in RADseq compared to microsatellite analysis, as increasing the number of markers will (regardless of the marker type) increase power and allow for clearer detection and higher resolution of genetic structure.

Differences in lipidome and metabolome organization of prefrontal cortex among human populations

Human populations, despite their overwhelming similarity, contain some distinct phenotypic, genetic, epigenetic, and gene expression features. In this study, we explore population differences at yet another level of molecular phenotype: the abundance of non-polar and polar low molecular weight compounds, lipids and metabolites in the prefrontal cortical region of the brain. We assessed the abundance of 1,670 lipids and 258 metabolites in 146 Han Chinese, 97 Western European, and 60 African American individuals of varying ages, covering most of the lifespan. The statistical analysis and logistic regression models both demonstrated extensive lipid and metabolic divergence of the Han Chinese individuals from the other two populations. This divergence was age-dependent, peaking in young adults, and involved metabolites and lipids clustering in specific metabolic pathways.

Utilizing field collected insects for next generation sequencing: Effects of sampling, storage, and DNA extraction methods

DNA sequencing technologies continue to advance the biological sciences, expanding opportunities for genomic studies of non-model organisms for basic and applied questions. Despite these opportunities, many next generation sequencing protocols have been developed assuming a substantial quantity of high molecular weight DNA (>100 ng), which can be difficult to obtain for many study systems. In particular, the ability to sequence field-collected specimens that exhibit varying levels of DNA degradation remains largely unexplored. In this study we investigate the influence of five traditional insect capture and curation methods on Double-Digest Restriction Enzyme Associated DNA (ddRAD) sequencing success for three wild bee species. We sequenced a total of 105 specimens (between 7-13 specimens per species and treatment). We additionally investigated how different DNA quality metrics (including pre-sequence concentration and contamination) predicted downstream sequencing success, and also compared two DNA extraction methods. We report successful library preparation for all specimens, with all treatments and extraction methods producing enough highly reliable loci for population genetic analyses. Although results varied between species, we found that specimens collected by net sampling directly into 100% EtOH, or by passive trapping followed by 100% EtOH storage before pinning tended to produce higher quality ddRAD assemblies, likely as a result of rapid specimen desiccation. Surprisingly, we found that specimens preserved in propylene glycol during field sampling exhibited lower-quality assemblies. We provide recommendations for each treatment, extraction method, and DNA quality assessment, and further encourage researchers to consider utilizing a wider variety of specimens for genomic analyses.

Whole-genome DNA similarity and population structure of Plasmodiophora brassicae strains from Canada

Background

Clubroot is an important disease of brassica crops world-wide. The causal agent, Plasmodiophora brassicae, has been present in Canada for over a century but was first identified on canola (Brassica napus) in Alberta, Canada in 2003. Genetic resistance to clubroot in an adapted canola cultivar has been available since 2009, but resistance breakdown was detected in 2013 and new pathotypes are increasing rapidly. Information on genetic similarity among pathogen populations across Canada could be useful in estimating the genetic variation in pathogen populations, predicting the effect of subsequent selection pressure on changes in the pathogen population over time, and even in identifying the origin of the initial pathogen introduction to canola in Alberta.

Results

The genomic sequences of 43 strains (34 field collections, 9 single-spore isolates) of P. brassicae from Canada, the United States, and China clustered into five clades based on SNP similarity. The strains from Canada separated into four clades, with two containing mostly strains from the Prairies (provinces of Alberta, Saskatchewan, and Manitoba) and two that were mostly from the rest of Canada or the USA. Several strains from China formed a separate clade. More than one pathotype and host were present in all four Canadian clades. The initial pathotypes from canola on the Prairies clustered separately from the pathotypes on canola that could overcome resistance to the initial pathotypes. Similarly, at one site in central Canada where resistance had broken down, about half of the genes differed (based on SNPs) between strains before and after the breakdown.

Conclusion

Clustering based on genome-wide DNA sequencing demonstrated that the initial pathotypes on canola on the Prairies clustered separately from the new virulent pathotypes on the Prairies. Analysis indicated that these ‘new’ pathotypes were likely present in the pathogen population at very low frequency, maintained through balancing selection, and increased rapidly in response to selection from repeated exposure to host resistance.

Fitness effects but no temperature-mediated balancing selection at the polymorphic Adh gene of Drosophila melanogaster

What factors maintain genetic variation in natural populations? Opposing selection pressures on protein stability and catalytic activity are thought to maintain variation along thermal gradients in many enzymes. We examined a classic hypothesis of temperature-mediated balancing selection, the alcohol dehydrogenase enzyme of Drosophila melanogaster, in which 2 latitudinally distributed variants are thought to be maintained by an activity/stability trade-off. Using in vitro and in vivo assays and population genetic analyses, we found no evidence of the predicted biochemical or fitness trade-offs and no signature of balancing selection. Rather, one variant confers greater activity and survival in the presence of ethanol, irrespective of temperature. Variation in Adh, and possibly other enzymes, must therefore be caused by other factors correlated with temperature.

Polymorphism in the alcohol dehydrogenase (ADH) protein of Drosophila melanogaster, like genetic variation in many other enzymes, has long been hypothesized to be maintained by a selective trade-off between thermostability and enzyme activity. Two major Adh variants, named Fast and Slow, are distributed along latitudinal clines on several continents. The balancing selection trade-off hypothesis posits that Fast is favored at high latitudes because it metabolizes alcohol faster, whereas Slow is favored at low latitudes because it is more stable at high temperatures. Here we use biochemical and physiological assays of precisely engineered genetic variants to directly test this hypothesis. As predicted, the Fast protein has higher catalytic activity than Slow, and both the Fast protein and regulatory variants linked to it confer greater ethanol tolerance on transgenic animals. But we found no evidence of a temperature-mediated trade-off: The Fast protein is not less stable or active at high temperatures, and Fast alleles increase ethanol tolerance and survivorship at all temperatures tested. Further, analysis of a population genomic dataset reveals no signature of balancing selection in the Adh gene. These results provide strong evidence against balancing selection driven by a stability/activity trade-off in Adh, and they justify caution about this hypothesis for other enzymes except those for which it has been directly tested. Our findings tentatively suggest that environment-specific selection for the Fast allele, coupled with demographic history, may have produced the observed pattern of Adh variation.

Speciation, gene flow, and seasonal migration in Catharus thrushes (Aves:Turdidae)

New World thrushes in the genus Catharus are small, insectivorous or omnivorous birds that have been used to explore several important questions in avian evolution, including the evolution of seasonal migration and plumage variation. Within Catharus, members of a clade of obligate long-distance migrants (C. fuscescens, C. minimus, and C. bicknelli) have also been used in the development of heteropatric speciation theory, a divergence process in which migratory lineages (which might occur in allopatry or sympatry during portions of their annual cycle) diverge despite low levels of gene flow. However, research on Catharus relationships has thus far been restricted to the use of small genetic datasets, which provide limited resolution of both phylogenetic and demographic histories. We used a large, multi-locus dataset from loci containing ultraconserved elements (UCEs) to study the demographic histories of the migratory C. fuscescens-minimus-bicknelli clade and to resolve the phylogeny of the migratory species of Catharus. Our dataset included more than 2000 loci and over 1700 variable genotyped sites, and analyses supported our prediction of divergence with gene flow in the fully migratory clade, with significant gene flow among all three species. Our phylogeny of the genus differs from past work in its placement of C. ustulatus, and further analyses suggest historic gene flow throughout the genus, producing genetically reticulate (or network) phylogenies. This raises questions about trait origins and suggests that seasonal migration and the resulting migratory condition of heteropatry is likely to promote hybridization not only during pairwise divergence and speciation, but also among non-sisters.

Genetic Variation for Ontogenetic Shifts in Metabolism Underlies Physiological Homeostasis in Drosophila

Organismal physiology emerges from metabolic pathways and subcellular structures like the mitochondria that can vary across development and among individuals. Here, we tested whether genetic variation at one level of physiology can be buffered at higher levels of biological organization during development by the inherent capacity for homeostasis in physiological systems. We found that the fundamental scaling relationship between mass and metabolic rate, as well as the oxidative capacity per mitochondria, changed significantly across development in the fruit fly Drosophila However, mitochondrial respiration rate was maintained at similar levels across development. Furthermore, larvae clustered into two types-those that switched to aerobic, mitochondrial ATP production before the second instar, and those that relied on anaerobic, glycolytic production of ATP through the second instar. Despite genetic variation for the timing of this metabolic shift, metabolic rate in second-instar larvae was more robust to genetic variation than was the metabolic rate of other instars. We found that larvae with a mitochondrial-nuclear incompatibility that disrupts mitochondrial function had increased aerobic capacity and relied more on anaerobic ATP production throughout development relative to larvae from wild-type strains. By taking advantage of both ways of making ATP, larvae with this mitochondrial-nuclear incompatibility maintained mitochondrial respiratory capacity, but also had higher levels of whole-body reactive oxygen species and decreased mitochondrial membrane potential, potentially as a physiological defense mechanism. Thus, genetic defects in core physiology can be buffered at the organismal level via physiological plasticity, and natural populations may harbor genetic variation for distinct metabolic strategies in development that generate similar organismal outcomes.

Contrasting patterns of neutral and functional genetic diversity in stable and disturbed environments

Genetic structure among and diversity within natural populations is influenced by a combination of ecological and evolutionary processes. These processes can differently influence neutral and functional genetic diversity and also vary according to environmental settings. To investigate the roles of interacting processes as drivers of population-level genetic diversity in the wild, we compared neutral and functional structure and diversity between 20 Tetrix undulata pygmy grasshopper populations in disturbed and stable habitats. Genetic differentiation was evident among the different populations, but there was no genetic separation between stable and disturbed environments. The incidence of long-winged phenotypes was higher in disturbed habitats, indicating that these populations were recently established by flight-capable colonizers. Color morph diversity and dispersion of outlier genetic diversity, estimated using AFLP markers, were higher in disturbed than in stable environments, likely reflecting that color polymorphism and variation in other functionally important traits increase establishment success. Neutral genetic diversity estimated using AFLP markers was lower in disturbed habitats, indicating stronger eroding effects on neutral diversity of genetic drift associated with founding events in disturbed compared to stable habitats. Functional diversity and neutral diversity were negatively correlated across populations, highlighting the utility of outlier loci in genetics studies and reinforcing that estimates of genetic diversity based on neutral markers do not infer evolutionary potential and the ability of populations and species to cope with environmental change.

Ultraconserved elements (UCEs) illuminate the population genomics of a recent, high-latitude avian speciation event

Using a large, consistent set of loci shared by descent (orthologous) to study relationships among taxa would revolutionize among-lineage comparisons of divergence and speciation processes. Ultraconserved elements (UCEs), highly conserved regions of the genome, offer such genomic markers. The utility of UCEs for deep phylogenetics is clearly established and there are mature analytical frameworks available, but fewer studies apply UCEs to recent evolutionary events, creating a need for additional example datasets and analytical approaches. We used UCEs to study population genomics in snow and McKay's buntings (Plectrophenax nivalis and P. hyperboreus). Prior work suggested divergence of these sister species during the last glacial maximum (∼18-74 Kya). With a sequencing depth of ∼30× from four individuals of each species, we used a series of analysis tools to genotype both alleles, obtaining a complete dataset of 2,635 variable loci (∼3.6 single nucleotide polymorphisms/locus) and 796 invariable loci. We found no fixed allelic differences between the lineages, and few loci had large allele frequency differences. Nevertheless, individuals were 100% diagnosable to species, and the two taxa were different genetically (F ST = 0.034; P = 0.03). The demographic model best fitting the data was one of divergence with gene flow. Estimates of demographic parameters differed from published mtDNA research, with UCE data suggesting lower effective population sizes (∼92,500-240,500 individuals), a deeper divergence time (∼241,000 years), and lower gene flow (2.8-5.2 individuals per generation). Our methods provide a framework for future population studies using UCEs, and our results provide additional evidence that UCEs are useful for answering questions at shallow evolutionary depths.

Disentangling genetic structure for genetic monitoring of complex populations

Genetic monitoring estimates temporal changes in population parameters from molecular marker information. Most populations are complex in structure and change through time by expanding or contracting their geographic range, becoming fragmented or coalescing, or increasing or decreasing density. Traditional approaches to genetic monitoring rely on quantifying temporal shifts of specific population metrics-heterozygosity, numbers of alleles, effective population size-or measures of geographic differentiation such as FST. However, the accuracy and precision of the results can be heavily influenced by the type of genetic marker used and how closely they adhere to analytical assumptions. Care must be taken to ensure that inferences reflect actual population processes rather than changing molecular techniques or incorrect assumptions of an underlying model of population structure. In many species of conservation concern, true population structure is unknown, or structure might shift over time. In these cases, metrics based on inappropriate assumptions of population structure may not provide quality information regarding the monitored population. Thus, we need an inference model that decouples the complex elements that define population structure from estimation of population parameters of interest and reveals, rather than assumes, fine details of population structure. Encompassing a broad range of possible population structures would enable comparable inferences across biological systems, even in the face of range expansion or contraction, fragmentation, or changes in density. Currently, the best candidate is the spatial Λ-Fleming-Viot (SLFV) model, a spatially explicit individually based coalescent model that allows independent inference of two of the most important elements of population structure: local population density and local dispersal. We support increased use of the SLFV model for genetic monitoring by highlighting its benefits over traditional approaches. We also discuss necessary future directions for model development to support large genomic datasets informing real-world management and conservation issues.

Decimation by sea star wasting disease and rapid genetic change in a keystone species, Pisaster ochraceus

Standing genetic variation enables or restricts a population's capacity to respond to changing conditions, including the extreme disturbances expected to increase in frequency and intensity with continuing anthropogenic climate change. However, we know little about how populations might respond to extreme events with rapid genetic shifts, or how population dynamics may influence and be influenced by population genomic change. We use a range-wide epizootic, sea star wasting disease, that onset in mid-2013 and caused mass mortality in Pisaster ochraceus to explore how a keystone marine species responded to an extreme perturbation. We integrated field surveys with restriction site-associated DNA sequencing data to (i) describe the population dynamics of mortality and recovery, and (ii) compare allele frequencies in mature P. ochraceus before the disease outbreak with allele frequencies in adults and new juveniles after the outbreak, to identify whether selection may have occurred. We found P. ochraceus suffered 81% mortality in the study region between 2012 and 2015, and experienced a concurrent 74-fold increase in recruitment beginning in late 2013. Comparison of pre- and postoutbreak adults revealed significant allele frequency changes at three loci, which showed consistent changes across the large majority of locations. Allele frequency shifts in juvenile P. ochraceus (spawned from premortality adults) were consistent with those seen in adult survivors. Such parallel shifts suggest detectable signals of selection and highlight the potential for persistence of this change in subsequent generations, which may influence the resilience of this keystone species to future outbreaks.

Genetic constraints on adaptation: a theoretical primer for the genomics era

Genetic constraints are features of inheritance systems that slow or prohibit adaptation. Several population genetic mechanisms of constraint have received sustained attention within the field since they were first articulated in the early 20th century. This attention is now reflected in a rich, and still growing, theoretical literature on the genetic limits to adaptive change. In turn, empirical research on constraints has seen a rapid expansion over the last two decades in response to changing interests of evolutionary biologists, along with new technologies, expanding data sets, and creative analytical approaches that blend mathematical modeling with genomics. Indeed, one of the most notable and exciting features of recent progress in genetic constraints is the close connection between theoretical and empirical research. In this review, we discuss five major population genetic contexts of genetic constraint: genetic dominance, pleiotropy, fitness trade-offs between types of individuals of a population, sign epistasis, and genetic linkage between loci. For each, we outline historical antecedents of the theory, specific contexts where constraints manifest, and their quantitative consequences for adaptation. From each of these theoretical foundations, we discuss recent empirical approaches for identifying and characterizing genetic constraints, each grounded and motivated by this theory, and outline promising areas for future work.