On the heterozygosity of an admixed population

Identification of SNPs in the second intron of IGF2BP1 and their Association with growth traits in Nanjiang Yellow goat

Insulin-like Growth Factor 2 mRNA-binding Protein 1 (IGF2BP1) is a candidate gene of significant interest for modulating economically important traits in livestock and poultry. The second intron of IGF2BP1 has been implicated in growth-related traits, though its precise mechanistic role remains elusive. Initial resequencing analyses in our laboratory indicated strong selective pressures on the IGF2BP1 genomic region, prompting the selection and identification of several single nucleotide polymorphisms (SNPs). Seven SNPs were mapped to the conserved region of the second intron, necessitating further investigation into their functional relevance and association with growth traits. In this study, 348 Nanjiang Yellow goats were analyzed, and the association analysis via the GLM program in SAS 9.4 identified five SNPs significantly correlated with growth traits. Notably, rs652062749(A > G) emerged as a critical locus influencing later-stage growth traits. Furthermore, strong linkage disequilibrium was observed among three SNPs, with the rs638185407 (T > A) variant markedly enhancing luciferase activity in H293T cells. Combination genotypes TTAACT, TTCCCC, and ATCACT were identified as superior for growth traits, offering theoretical insights for genetic co-breeding. This study underscores the potential utility of IGF2BP1 as a functional genetic marker in Nanjiang Yellow goat breeding programs.

Determining population structure from k-mer frequencies

Background

Understanding population structure within species provides information on connections among different populations and how they evolve over time. This knowledge is important for studies ranging from evolutionary biology to large-scale variant-trait association studies. Current approaches to determining population structure include model-based approaches, statistical approaches, and distance-based ancestry inference approaches.

Methods

In this work, we identify population structure from DNA sequence data using an alignment-free approach. We use the frequencies of short DNA substrings from across the genome (k-mers) with principal component analysis (PCA). K-mer frequencies can be viewed as a summary statistic of a genome and have the advantage of being easily derived from a genome by counting the number of times a k-mer occurred in a sequence. In contrast, most population structure work employing PCA uses multi-locus genotype data (SNPs, microsatellites, or haplotypes). No genetic assumptions must be met to generate k-mers, whereas current population structure approaches often depend on several genetic assumptions and can require careful selection of ancestry informative markers to identify populations. We compare our k-mer based approach to population structure estimated using SNPs with both empirical and simulated data.

Results

In this work, we show that PCA is able to determine population structure just from the frequency of k-mers found in the genome. The application of PCA and a clustering algorithm to k-mer profiles of genomes provides an easy approach to detecting the number and composition of populations (clusters) present in the dataset. Using simulations, we show that results are at least comparable to population structure estimates using SNPs. When using human genomes from populations identified by the 1000 Genomes Project, the results are better than population structure estimates using SNPs from the same samples, and comparable to those found by a model-based approach using genetic markers from larger numbers of samples.

Conclusions

This study shows that PCA, together with the clustering algorithm, is able to detect population structure from k-mer frequencies and can separate samples of admixed and non-admixed origin. Using k-mer frequencies to determine population structure has the potential to avoid some challenges of existing methods and may even improve on estimates from small samples.

Patterns of Genetic Diversity and Gene Flow Associated With an Aridity Gradient in Populations of Common Mole-rats, Cryptomys hottentotus hottentotus

Genetic adaptation is the change of a population toward a phenotype that best fits the present ecological conditions of the environment it inhabits. As environmental conditions change, allele frequencies shift, resulting in different populations of the same species possessing genetic variation and divergent phenotypes. Cooperatively breeding common mole-rats (Cryptomys hottentotus hottentotus) inhabit environments along an aridity gradient in South Africa, which provides an opportunity for local genetic adaptations to occur. Using one mitochondrial gene (cytochrome b) and 3,540 SNP loci across the whole genome, we determined the phylogenetic relationship, population structure and genetic diversity of five populations of C. h. hottentotus located along an aridity gradient. Mitochondrial data identified population-specific clades that were less distinct in the two mesic populations, potentially indicating historical or recent gene flow, or the retention of ancestral haplotypes. Arid and semi-arid populations formed a distinct cluster from the non-arid populations. Genetic diversity and gene flow were higher in arid-dwelling individuals, suggesting greater connectivity and interactions between colonies in arid regions in comparison to mesic ones. Using an Aridity Index, we determined that isolation by environment, rather than isolation by geographical distance, best explains the genetic distance between the populations. Further analyses using target loci may determine if there are differing underlying genetic adaptations among populations of C. h. hottentotus. These analyses could help unravel population differences in response to environmental factors within a subspecies of bathyergid mole-rat and determine the adaptive capacity of this small nonmigratory subterranean rodent species in response to aridification in the face of climate change.

Population genomics provides insights into the genetic diversity and adaptation of the Pieris rapae in China

The cabbage white butterfly (Pieris rapae), a major agricultural pest, has become one of the most abundant and destructive butterflies in the world. It is widely distributed in a large variety of climates and terrains of China due to its strong adaptability. To gain insight into the population genetic characteristics of P. rapae in China, we resequenced the genome of 51 individuals from 19 areas throughout China. Using population genomics approaches, a dense variant map of P. rapae was observed, indicating a high level of polymorphism that could result in adaptation to a changing environment. The feature of the genetic structure suggested considerable genetic admixture in different geographical groups. Additionally, our analyses suggest that physical barriers may have played a more important role than geographic distance in driving genetic differentiation. Population history showed the effective population size of P. rapae was greatly affected by global temperature changes, with mild periods (i.e., temperatures warmer than those during glaciation but not excessively hot) leading to an increase in population size. Furthermore, by comparing populations from south and north China, we have identified selected genes related to sensing temperature, growth, neuromodulation and immune response, which may reveal the genetic basis of adaptation to different environments. Our study is the first to illustrate the genetic signatures of P. rapae in China at the population genomic level, providing fundamental knowledge of the genetic diversity and adaptation of P. rapae.

Hybridization of Atlantic puffins in the Arctic coincides with 20th-century climate change

The Arctic is experiencing the fastest rates of global warming, leading to shifts in the distribution of its biota and increasing the potential for hybridization. However, genomic evidence of recent hybridization events in the Arctic remains unexpectedly rare. Here, we use whole-genome sequencing of contemporary and 122-year-old historical specimens to investigate the origin of an Arctic hybrid population of Atlantic puffins (Fratercula arctica) on Bjørnøya, Norway. We show that the hybridization between the High Arctic, large-bodied subspecies F. a. naumanni and the temperate, smaller-sized subspecies F. a. arctica began as recently as six generations ago due to an unexpected southward range expansion of F. a. naumanni. Moreover, we find a significant temporal loss of genetic diversity across Arctic and temperate puffin populations. Our observations provide compelling genomic evidence of the impacts of recent distributional shifts and loss of diversity in Arctic communities during the 20th century.

The Impact of Modern Admixture on Archaic Human Ancestry in Human Populations

Admixture, the genetic merging of parental populations resulting in mixed ancestry, has occurred frequently throughout the course of human history. Numerous admixture events have occurred between human populations across the world, which have shaped genetic ancestry in modern humans. For example, populations in the Americas are often mosaics of different ancestries due to recent admixture events as part of European colonization. Admixed individuals also often have introgressed DNA from Neanderthals and Denisovans that may have come from multiple ancestral populations, which may affect how archaic ancestry is distributed across an admixed genome. In this study, we analyzed admixed populations from the Americas to assess whether the proportion and location of admixed segments due to recent admixture impact an individual's archaic ancestry. We identified a positive correlation between non-African ancestry and archaic alleles, as well as a slight increase of Denisovan alleles in Indigenous American segments relative to European segments in admixed genomes. We also identify several genes as candidates for adaptive introgression, based on archaic alleles present at high frequency in admixed American populations but low frequency in East Asian populations. These results provide insights into how recent admixture events between modern humans redistributed archaic ancestry in admixed genomes.

Genomic evidence for homoploid hybrid speciation in a marine mammal apex predator

Hybridization is widespread and constitutes an important source of genetic variability and evolution. In animals, its role in generating novel and independent lineages (hybrid speciation) has been strongly debated, with only a few cases supported by genomic data. The South American fur seal (SAfs) Arctocephalus australis is a marine apex predator of Pacific and Atlantic waters, with a disjunct set of populations in Peru and Northern Chile [Peruvian fur seal (Pfs)] with controversial taxonomic status. We demonstrate, using complete genome and reduced representation sequencing, that the Pfs is a genetically distinct species with an admixed genome that originated from hybridization between the SAfs and the Galapagos fur seal (Arctocephalus galapagoensis) ~400,000 years ago. Our results strongly support the origin of Pfs by homoploid hybrid speciation over alternative introgression scenarios. This study highlights the role of hybridization in promoting species-level biodiversity in large vertebrates.

Deleterious mutation load in the admixed mice population

Deleterious mutation loads are known to correlate negatively with effective population size (Ne). Due to this reason, previous studies observed a higher proportion of harmful mutations in small populations than that in large populations. However, the mutational load in an admixed population that derived from introgression between individuals from two populations with vastly different Ne is not known. We investigated this using the whole genome data from two subspecies of the mouse (Mus musculus castaneus and Mus musculus musculus) with significantly different Ne. We used the ratio of diversities at nonsynonymous and synonymous sites (dN/dS) to measure the harmful mutation load. Our results showed that this ratio observed for the admixed population was intermediate between those of the parental populations. The dN/dS ratio of the hybrid population was significantly higher than that of M. m. castaneus but lower than that of M. m. musculus. Our analysis revealed a significant positive correlation between the proportion of M. m. musculus ancestry in admixed individuals and their dN/dS ratio. This suggests that the admixed individuals with high proportions of M. m. musculus ancestry have large dN/dS ratios. We also used the proportion of deleterious nonsynonymous SNVs as a proxy for deleterious mutation load, which also produced similar results. The observed results were in concordance with those expected by theory. We also show a shift in the distribution of fitness effects of nonsynonymous SNVs in the admixed genomes compared to the parental populations. These findings suggest that the deleterious mutation load of the admixed population is determined by the proportion of the ancestries of the subspecies. Therefore, it is important to consider the status and the level of genetic admixture of the populations whilst estimating the mutation loads.

The impact of modern admixture on archaic human ancestry in human populations

Admixture, the genetic merging of parental populations resulting in mixed ancestry, has occurred frequently throughout the course of human history. Numerous admixture events have occurred between human populations across the world, as well as introgression between humans and archaic humans, Neanderthals and Denisovans. One example are genomes from populations in the Americas, as these are often mosaics of different ancestries due to recent admixture events as part of European colonization. In this study, we analyzed admixed populations from the Americas to assess whether the proportion and location of admixed segments due to recent admixture impact an individual’s archaic ancestry. We identified a positive correlation between non-African ancestry and archaic alleles, as well as a slight enrichment of Denisovan alleles in Indigenous American segments relative to European segments in admixed genomes. We also identify several genes as candidates for adaptive introgression, based on archaic alleles present at high frequency in admixed American populations but low frequency in East Asian populations. These results provide insights into how recent admixture events between modern humans redistributed archaic ancestry in admixed genomes.

A central limit theorem concerning uncertainty in estimates of individual admixture

The concept of individual admixture (IA) assumes that the genome of individuals is composed of alleles inherited from K ancestral populations. Each copy of each allele has the same chance q_k to originate from population k, and together with the allele frequencies p in all populations at all M markers, comprises the admixture model. Here, we assume a supervised scheme, i.e. allele frequencies p are given through a reference database of size N, and q is estimated via maximum likelihood for a single sample. We study laws of large numbers and central limit theorems describing effects of finiteness of both, M and N, on the estimate of q. We recall results for the effect of finite M, and provide a central limit theorem for the effect of finite N, introduce a new way to express the uncertainty in estimates in standard barplots, give simulation results, and discuss applications in forensic genetics.

Human genetic admixture through the lens of population genomics

Over the past 50 years, geneticists have made great strides in understanding how our species' evolutionary history gave rise to current patterns of human genetic diversity classically summarized by Lewontin in his 1972 paper, 'The Apportionment of Human Diversity'. One evolutionary process that requires special attention in both population genetics and statistical genetics is admixture: gene flow between two or more previously separated source populations to form a new admixed population. The admixture process introduces ancestry-based structure into patterns of genetic variation within and between populations, which in turn influences the inference of demographic histories, identification of genetic targets of selection and prediction of complex traits. In this review, we outline some challenges for admixture population genetics, including limitations of applying methods designed for populations without recent admixture to the study of admixed populations. We highlight recent studies and methodological advances that aim to overcome such challenges, leveraging genomic signatures of admixture that occurred in the past tens of generations to gain insights into human history, natural selection and complex trait architecture. This article is part of the theme issue 'Celebrating 50 years since Lewontin's apportionment of human diversity'.

Genetic Diversity and Population Structure of Portunustrituberculatus in Released and Wild Populations Based on Microsatellite DNA Markers from the Yangtze Estuary

Portunus trituberculatus is an important economic species of crab that is artificially bred and released in the Yangtze River Estuary and its adjacent sea areas. Based on six microsatellite markers, we investigate the genetic diversity and structure of 101 P. trituberculatus specimens collected from two hatcheries in Nantong and Zhoushan that participated in stock enhancement in the year 2019. We compared these with 124 wild specimens caught from 13 localities in the estuary. Analysis of several genetic diversity parameters (NA, RS, I, HO, HE, FIS, and FST) for the 15 populations demonstrates that both released and wild populations possess relatively rich genetic diversity. Furthermore, the released groups demonstrate no less genetic variation between themselves than do the wild crabs. Most FIS values are greater than zero, which shows inbreeding is common among specimens with geographically open sites. However, insufficient sampling may have led to a wide distribution of null alleles, a Hardy–Weinberg test disequilibrium in microsatellite markers PN22 and P04, and a lack of crab genetic diversity in site 14. All populations (except locality 14) have not suffered the bottleneck effect. Four subgroups can be seen to roughly spread longitudinally along the sample area by performing pairwise comparisons of genetic distance and FST values among the populations. No obvious topological heterogeneity is discovered among the four subgroups in a phylogenetic tree. The existence of genetic exchange and differentiation among the subgroups is also verified using structure analysis. Therefore, based on this evidence, we propose that the hatchery stock enhancements performed in Nantong and Zhoushan result in no reduction in genetic diversity for wild populations in the Yangtze Estuary in 2019.

A geometric relationship of F 2 , F 3 and F 4 -statistics with principal component analysis

Principal component analysis (PCA) and F -statistics sensu Patterson are two of the most widely used population genetic tools to study human genetic variation. Here, I derive explicit connections between the two approaches and show that these two methods are closely related. F -statistics have a simple geometrical interpretation in the context of PCA, and orthogonal projections are a key concept to establish this link. I show that for any pair of populations, any population that is admixed as determined by an F 3 -statistic will lie inside a circle on a PCA plot. Furthermore, the F 4 -statistic is closely related to an angle measurement, and will be zero if the differences between pairs of populations intersect at a right angle in PCA space. I illustrate my results on two examples, one of Western Eurasian, and one of global human diversity. In both examples, I find that the first few PCs are sufficient to approximate most F -statistics, and that PCA plots are effective at predicting F -statistics. Thus, while F -statistics are commonly understood in terms of discrete populations, the geometric perspective illustrates that they can be viewed in a framework of populations that vary in a more continuous manner.

This article is part of the theme issue ‘Celebrating 50 years since Lewontin's apportionment of human diversity’.

An Overview of Strategies for Detecting Genotype-Phenotype Associations Across Ancestrally Diverse Populations

Genome-wide association studies (GWAS) have been very successful at identifying genetic variants influencing a large number of traits. Although the great majority of these studies have been performed in European-descent individuals, it has been recognised that including populations with differing ancestries enhances the potential for identifying causal SNPs due to their differing patterns of linkage disequilibrium. However, when individuals from distinct ethnicities are included in a GWAS, it is necessary to implement a number of control steps to ensure that the identified associations are real genotype-phenotype relationships. In this Review, we discuss the analyses that are required when performing multi-ethnic studies, including methods for determining ancestry at the global and local level for sample exclusion, controlling for ancestry in association testing, and post-GWAS interrogation methods such as genomic control and meta-analysis. We hope that this overview provides a primer for those researchers interested in including distinct populations in their studies.

A century of genetic homogenization in Baltic salmon-evidence from archival DNA

Intra-species genetic homogenization arising from anthropogenic impacts is a major threat to biodiversity. However, few taxa have sufficient historical material to systematically quantify long-term genetic changes. Using archival DNA collected over approximately 100 years, we assessed spatio-temporal genetic change in Atlantic salmon populations across the Baltic Sea, an area heavily impacted by hydropower exploitation and associated with large-scale mitigation stocking. Analysis was carried out by screening 82 SNPs in 1680 individuals from 13 Swedish rivers. We found an overall decrease in genetic divergence and diminished isolation by distance among populations, strongly indicating genetic homogenization over the past century. We further observed an increase in genetic diversity within populations consistent with increased gene flow. The temporal genetic change was lower in larger wild populations than in smaller wild and hatchery-reared ones, indicating that larger populations have been able to support a high number of native spawners in relation to immigrants. Our results demonstrate that stocking practices of salmon in the Baltic Sea have led to the homogenization of populations over the last century, potentially compromising their ability to adapt to environmental change. Stocking of reared fish is common worldwide, and our study is a cautionary example of the potentially long-term negative effects of such activities.