Search

Understanding natural selection and other forms of non-neutrality is a major focus for the use of machine learning in population genetics. Existing methods rely on computationally intensive simulated training data. Unlike efficient neutral coalescent simulations for demographic inference, realistic simulations of selection typically require slow forward simulations. Because there are many possible modes of selection, a high dimensional parameter space must be explored, with no guarantee that the simulated models are close to the real processes. Finally, it is difficult to interpret trained neural networks, leading to a lack of understanding about what features contribute to classification. Here we develop a new approach to detect selection and other local evolutionary processes that requires relatively few selection simulations during training. We build upon a generative adversarial network trained to simulate realistic neutral data. This consists of a generator (fitted demographic model), and a discriminator (convolutional neural network) that predicts whether a genomic region is real or fake. As the generator can only generate data under neutral demographic processes, regions of real data that the discriminator recognizes as having a high probability of being "real" do not fit the neutral demographic model and are therefore candidates for targets of selection. To incentivize identification of a specific mode of selection, we fine-tune the discriminator with a small number of custom non-neutral simulations. We show that this approach has high power to detect various forms of selection in simulations, and that it finds regions under positive selection identified by state-of-the-art population genetic methods in three human populations. Finally, we show how to interpret the trained networks by clustering hidden units of the discriminator based on their correlation patterns with known summary statistics., Competing Interests: Conflicts of interest The author(s) declare no conflicts of interest., (© The Author(s) 2024. Published by Oxford University Press on behalf of The Genetics Society of America.)

Ancient DNA can directly reveal the contribution of natural selection to human genomic variation. However, while the analysis of ancient DNA has been successful at identifying genomic signals of selection, inferring the phenotypic consequences of that selection has been more difficult. Most trait-associated variants are noncoding, so we expect that a large proportion of the phenotypic effects of selection will also act through noncoding variation. Since we cannot measure gene expression directly in ancient individuals, we used an approach (Joint-Tissue Imputation [JTI]) developed to predict gene expression from genotype data. We tested for changes in the predicted expression of 17,384 protein coding genes over a time transect of 4,500 years using 91 present-day and 616 ancient individuals from Britain. We identified 28 genes at seven genomic loci with significant (false discovery rate [FDR] < 0.05) changes in predicted expression levels in this time period. We compared the results from our transcriptome-wide scan to a genome-wide scan based on estimating per-single nucleotide polymorphism (SNP) selection coefficients from time series data. At five previously identified loci, our approach allowed us to highlight small numbers of genes with evidence for significant shifts in expression from peaks that in some cases span tens of genes. At two novel loci (SLC44A5 and NUP85), we identify selection on gene expression not captured by scans based on genomic signatures of selection. Finally, we show how classical selection statistics (iHS and SDS) can be combined with JTI models to incorporate functional information into scans that use present-day data alone. These results demonstrate the potential of this type of information to explore both the causes and consequences of natural selection., (© The Author(s) 2024. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

The rules and structure of human culture impact health as much as genetics or environment. To study these relationships, we combine ancient DNA (n = 230), skeletal metrics (n = 391), palaeopathology (n = 606) and dietary stable isotopes (n = 873) to analyse stature variation in Early Neolithic Europeans from North Central, South Central, Balkan and Mediterranean regions. In North Central Europe, stable isotopes and linear enamel hypoplasias indicate high environmental stress across sexes, but female stature is low, despite polygenic scores identical to males, and suggests that cultural factors preferentially supported male recovery from stress. In Mediterranean populations, sexual dimorphism is reduced, indicating male vulnerability to stress and no strong cultural preference for males. Our analysis indicates that biological effects of sex-specific inequities can be linked to cultural influences at least as early as 7,000 yr ago, and culture, more than environment or genetics, drove height disparities in Early Neolithic Europe., (© 2023. The Author(s), under exclusive licence to Springer Nature Limited.)

To elucidate the population history of the Caucasus, we conducted a survey of genetic diversity in Samegrelo (Mingrelia), western Georgia. We collected DNA samples and genealogical information from 485 individuals residing in 30 different locations, the vast majority of whom being Mingrelian speaking. From these DNA samples, we generated mitochondrial DNA (mtDNA) control region sequences for all 485 participants (female and male), Y-short tandem repeat haplotypes for the 372 male participants, and analyzed all samples at nearly 590,000 autosomal single nucleotide polymorphisms (SNPs) plus around 33,000 on the sex chromosomes, with 27,000 SNP removed for missingness, using the GenoChip 2.0+ microarray. The resulting data were compared with those from populations from Anatolia, the Caucasus, the Near East, and Europe. Overall, Mingrelians exhibited considerable mtDNA haplogroup diversity, having high frequencies of common West Eurasian haplogroups (H, HV, I, J, K, N1, R1, R2, T, U, and W. X2) and low frequencies of East Eurasian haplogroups (A, C, D, F, and G). From a Y-chromosome standpoint, Mingrelians possessed a variety of haplogroups, including E1b1b, G2a, I2, J1, J2, L, Q, R1a, and R1b. Analysis of autosomal SNP data further revealed that Mingrelians are genetically homogeneous and cluster with other modern-day South Caucasus populations. When compared with ancient DNA samples from Bronze Age archaeological contexts in the broader region, these data indicate that the Mingrelian gene pool began taking its current form at least by this period, probably in conjunction with the formation of a distinct linguistic community., (© The Author(s) 2023. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

A new study aims to identify how genetic and physiological adaptations to altitude affect pregnancy, childbirth and neonatal health in one of the most extreme environments on Earth, the Tibetan Plateau., Competing Interests: Declaration of interests The author declares no competing interests., (Copyright © 2023 Elsevier Inc. All rights reserved.)

Ancient DNA can directly reveal the contribution of natural selection to human genomic variation. However, while the analysis of ancient DNA has been successful at identifying genomic signals of selection, inferring the phenotypic consequences of that selection has been more difficult. Most trait-associated variants are non-coding, so we expect that a large proportion of the phenotypic effects of selection will also act through non-coding variation. Since we cannot measure gene expression directly in ancient individuals, we used an approach ( Joint-Tissue Imputation ; JTI) developed to predict gene expression from genotype data. We tested for changes in the predicted expression of 17,384 protein coding genes over a time transect of 4500 years using 91 present-day and 616 ancient individuals from Britain. We identified 28 genes at seven genomic loci with significant (FDR < 0.05) changes in predicted expression levels in this time period. We compared the results from our transcriptome-wide scan to a genome-wide scan based on estimating per-SNP selection coefficients from time series data. At five previously identified loci, our approach allowed us to highlight small numbers of genes with evidence for significant shifts in expression from peaks that in some cases span tens of genes. At two novel loci ( SLC44A5 and NUP85) , we identify selection on gene expression not captured by scans based on genomic signatures of selection. Finally we show how classical selection statistics (iHS and SDS) can be combined with JTI models to incorporate functional information into scans that use present-day data alone. These results demonstrate the potential of this type of information to explore both the causes and consequences of natural selection., Competing Interests: Competing Interests The authors report no conflicts of interest.

Understanding natural selection in humans and other species is a major focus for the use of machine learning in population genetics. Existing methods rely on computationally intensive simulated training data. Unlike efficient neutral coalescent simulations for demographic inference, realistic simulations of selection typically requires slow forward simulations. Because there are many possible modes of selection, a high dimensional parameter space must be explored, with no guarantee that the simulated models are close to the real processes. Mismatches between simulated training data and real test data can lead to incorrect inference. Finally, it is difficult to interpret trained neural networks, leading to a lack of understanding about what features contribute to classification. Here we develop a new approach to detect selection that requires relatively few selection simulations during training. We use a Generative Adversarial Network (GAN) trained to simulate realistic neutral data. The resulting GAN consists of a generator (fitted demographic model) and a discriminator (convolutional neural network). For a genomic region, the discriminator predicts whether it is "real" or "fake" in the sense that it could have been simulated by the generator. As the "real" training data includes regions that experienced selection and the generator cannot produce such regions, regions with a high probability of being real are likely to have experienced selection. To further incentivize this behavior, we "fine-tune" the discriminator with a small number of selection simulations. We show that this approach has high power to detect selection in simulations, and that it finds regions under selection identified by state-of-the art population genetic methods in three human populations. Finally, we show how to interpret the trained networks by clustering hidden units of the discriminator based on their correlation patterns with known summary statistics. In summary, our approach is a novel, efficient, and powerful way to use machine learning to detect natural selection.

Germline mutation is the mechanism by which genetic variation in a population is created. Inferences derived from mutation rate models are fundamental to many population genetics methods. Previous models have demonstrated that nucleotides flanking polymorphic sites-the local sequence context-explain variation in the probability that a site is polymorphic. However, limitations to these models exist as the size of the local sequence context window expands. These include a lack of robustness to data sparsity at typical sample sizes, lack of regularization to generate parsimonious models and lack of quantified uncertainty in estimated rates to facilitate comparison between models. To address these limitations, we developed Baymer, a regularized Bayesian hierarchical tree model that captures the heterogeneous effect of sequence contexts on polymorphism probabilities. Baymer implements an adaptive Metropolis-within-Gibbs Markov Chain Monte Carlo sampling scheme to estimate the posterior distributions of sequence-context based probabilities that a site is polymorphic. We show that Baymer accurately infers polymorphism probabilities and well-calibrated posterior distributions, robustly handles data sparsity, appropriately regularizes to return parsimonious models, and scales computationally at least up to 9-mer context windows. We demonstrate application of Baymer in three ways-first, identifying differences in polymorphism probabilities between continental populations in the 1000 Genomes Phase 3 dataset, second, in a sparse data setting to examine the use of polymorphism models as a proxy for de novo mutation probabilities as a function of variant age, sequence context window size, and demographic history, and third, comparing model concordance between different great ape species. We find a shared context-dependent mutation rate architecture underlying our models, enabling a transfer-learning inspired strategy for modeling germline mutations. In summary, Baymer is an accurate polymorphism probability estimation algorithm that automatically adapts to data sparsity at different sequence context levels, thereby making efficient use of the available data., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2023 Adams et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Most variants identified in human genome-wide association studies and scans for selection are noncoding. Interpretation of their effects and the way in which they contribute to phenotypic variation and adaptation in human populations is therefore limited by our understanding of gene regulation and the difficulty of confidently linking noncoding variants to genes. To overcome this, we developed a gene-wise test for population-specific selection based on combinations of regulatory variants. Specifically, we use the QX statistic to test for polygenic selection on cis-regulatory variants based on whether the variance across populations in the predicted expression of a particular gene is higher than expected under neutrality. We then applied this approach to human data, testing for selection on 17,388 protein-coding genes in 26 populations from the Thousand Genomes Project. We identified 45 genes with significant evidence (FDR<0.1) for selection, including FADS1, KHK, SULT1A2, ITGAM, and several genes in the HLA region. We further confirm that these signals correspond to plausible population-level differences in predicted expression. While the small number of significant genes (0.2%) is consistent with most cis-regulatory variation evolving under genetic drift or stabilizing selection, it remains possible that there are effects not captured in this study. Our gene-level QX score is independent of standard genomic tests for selection, and may therefore be useful in combination with traditional selection scans to specifically identify selection on regulatory variation. Overall, our results demonstrate the utility of combining population-level genomic data with functional data to understand the evolution of gene expression., (© The Author(s) 2023. Published by Oxford University Press on behalf of The Genetics Society of America.)

Mitochondrial DNA copy number (mtCN) is often treated as a proxy for mitochondrial (dys-) function and disease risk. Pathological changes in mtCN are common symptoms of rare mitochondrial disorders, but reported associations between mtCN and common diseases vary across studies. To understand the biology of mtCN, we carried out genome- and phenome-wide association studies of mtCN in 30,666 individuals from the Penn Medicine BioBank (PMBB)-a diverse cohort of largely African and European ancestry. We estimated mtCN in peripheral blood using exome sequence data, taking cell composition into account. We replicated known genetic associations of mtCN in the PMBB and found that their effects are highly correlated between individuals of European and African ancestry. However, the heritability of mtCN was much higher among individuals of largely African ancestry ( h 2 = 0.3 ) compared with European ancestry individuals ( h 2 = 0.1 ) . Admixture mapping suggests that there are undiscovered variants underlying mtCN that are differentiated in frequency between individuals with African and European ancestry. We show that mtCN is associated with many health-related phenotypes. We discovered robust associations between mtDNA copy number and diseases of metabolically active tissues, such as cardiovascular disease and liver damage, that were consistent across African and European ancestry individuals. Other associations, such as epilepsy and prostate cancer, were only discovered in either individuals with European or African ancestry but not both. We show that mtCN-phenotype associations can be sensitive to blood cell composition and environmental modifiers, explaining why such associations are inconsistent across studies. Thus, mtCN-phenotype associations must be interpreted with care., (© 2023 The Author(s).)

Identifying genetic determinants of reproductive success may highlight mechanisms underlying fertility and identify alleles under present-day selection. Using data in 785,604 individuals of European ancestry, we identified 43 genomic loci associated with either number of children ever born (NEB) or childlessness. These loci span diverse aspects of reproductive biology, including puberty timing, age at first birth, sex hormone regulation, endometriosis and age at menopause. Missense variants in ARHGAP27 were associated with higher NEB but shorter reproductive lifespan, suggesting a trade-off at this locus between reproductive ageing and intensity. Other genes implicated by coding variants include PIK3IP1, ZFP82 and LRP4, and our results suggest a new role for the melanocortin 1 receptor (MC1R) in reproductive biology. As NEB is one component of evolutionary fitness, our identified associations indicate loci under present-day natural selection. Integration with data from historical selection scans highlighted an allele in the FADS1/2 gene locus that has been under selection for thousands of years and remains so today. Collectively, our findings demonstrate that a broad range of biological mechanisms contribute to reproductive success., (© 2023. The Author(s), under exclusive licence to Springer Nature Limited.)

Ancient DNA has revealed multiple episodes of admixture in human prehistory during geographic expansions associated with cultural innovations. One important example is the expansion of Neolithic agricultural groups out of the Near East into Europe and their consequent admixture with Mesolithic hunter-gatherers. 1 , 2 , 3 , 4 Ancient genomes from this period provide an opportunity to study the role of admixture in providing new genetic variation for selection to act upon, and also to identify genomic regions that resisted hunter-gatherer introgression and may thus have contributed to agricultural adaptations. We used genome-wide DNA from 677 individuals spanning Mesolithic and Neolithic Europe to infer ancestry deviations in the genomes of admixed individuals and to test for natural selection after admixture by testing for deviations from a genome-wide null distribution. We find that the region around the pigmentation-associated gene SLC24A5 shows the greatest overrepresentation of Neolithic local ancestry in the genome (|Z| = 3.46). In contrast, we find the greatest overrepresentation of Mesolithic ancestry across the major histocompatibility complex (MHC; |Z| = 4.21), a major immunity locus, which also shows allele frequency deviations indicative of selection following admixture (p = 1 × 10 -56 ). This could reflect negative frequency-dependent selection on MHC alleles common in Neolithic populations or that Mesolithic alleles were positively selected for and facilitated adaptation in Neolithic populations to pathogens or other environmental factors. Our study extends previous results that highlight immune function and pigmentation as targets of adaptation in more recent populations to selection processes in the Stone Age., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2023 The Author(s). Published by Elsevier Inc. All rights reserved.)

Klunk et al. analyzed ancient DNA data from individuals in London and Denmark before, during and after the Black Death [1], and argued that allele frequency changes at immune genes were too large to be produced by random genetic drift and thus must reflect natural selection. They also identified four specific variants that they claimed show evidence of selection including at ERAP2 , for which they estimate a selection coefficient of 0.39-several times larger than any selection coefficient on a common human variant reported to date. Here we show that these claims are unsupported for four reasons. First, the signal of enrichment of large allele frequency changes in immune genes comparing people in London before and after the Black Death disappears after an appropriate randomization test is carried out: the P value increases by ten orders of magnitude and is no longer significant. Second, a technical error in the estimation of allele frequencies means that none of the four originally reported loci actually pass the filtering thresholds. Third, the filtering thresholds do not adequately correct for multiple testing. Finally, in the case of the ERAP2 variant rs2549794, which Klunk et al. show experimentally may be associated with a host interaction with Y. pestis , we find no evidence of significant frequency change either in the data that Klunk et al. report, or in published data spanning 2,000 years. While it remains plausible that immune genes were subject to natural selection during the Black Death, the magnitude of this selection and which specific genes may have been affected remains unknown.

Loss or absence of hearing is common at both extremes of human lifespan, in the forms of congenital deafness and age-related hearing loss. While these are often studied separately, there is increasing evidence that their genetic basis is at least partially overlapping. In particular, both common and rare variants in genes associated with monogenic forms of hearing loss also contribute to the more polygenic basis of age-related hearing loss. Here, we directly test this model in the Penn Medicine BioBank-a healthcare system cohort of around 40,000 individuals with linked genetic and electronic health record data. We show that increased burden of predicted deleterious variants in Mendelian hearing loss genes is associated with increased risk and severity of adult-onset hearing loss. As a specific example, we identify one gene-TCOF1, responsible for a syndromic form of congenital hearing loss-in which deleterious variants are also associated with adult-onset hearing loss. We also identify four additional novel candidate genes (COL5A1, HMMR, RAPGEF3, and NNT) in which rare variant burden may be associated with hearing loss. Our results confirm that rare variants in Mendelian hearing loss genes contribute to polygenic risk of hearing loss, and emphasize the utility of healthcare system cohorts to study common complex traits and diseases., Competing Interests: I have read the journal’s policy and the authors of this manuscript have the following competing interests: Regeneron Genetics Center is part of Regeneron Pharmaceuticals. The other authors have declared that no competing interests exist., (Copyright: © 2023 Hui et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

We developed a novel method for efficiently estimating time-varying selection coefficients from genome-wide ancient DNA data. In simulations, our method accurately recovers selective trajectories and is robust to misspecification of population size. We applied it to a large data set of ancient and present-day human genomes from Britain and identified seven loci with genome-wide significant evidence of selection in the past 4500 yr. Almost all of them can be related to increased vitamin D or calcium levels, suggesting strong selective pressure on these or related phenotypes. However, the strength of selection on individual loci varied substantially over time, suggesting that cultural or environmental factors moderated the genetic response. Of 28 complex anthropometric and metabolic traits, skin pigmentation was the only one with significant evidence of polygenic selection, further underscoring the importance of phenotypes related to vitamin D. Our approach illustrates the power of ancient DNA to characterize selection in human populations and illuminates the recent evolutionary history of Britain., (© 2022 Mathieson and Terhorst; Published by Cold Spring Harbor Laboratory Press.)

Objective: The aim of this retrospective study was to review the clinical data and outcomes of patients that suffered ectopic cilium (EC)., Animals Studied: One hundred and twelve dogs from multiple private practices in France, with a clinical diagnosis of EC were included in the study., Results: The mean age of affected dogs was 2.3 years. There were 64 females and 48 males. The most represented breeds were the Shi Tzu, the French Bulldog, the English Bulldog and the Chihuahua. Eleven dogs were affected bilaterally. The upper eyelid was implicated in 93.5% of the cases, with the median portion being the most affected. No statistical difference was observed between the right and the left eye. EC were associated with distichiasis in 50% of the cases. Pigmentation of the conjunctiva at the point of exit of the EC was present in 58% of the cases. EC were short in 75% and long in 25% of the cases. Corneal complications were statistically associated with short EC. The corneal lesions associated with EC were keratitis (94%), corneal granuloma (0.8%), corneal fibrosis (2.7%), corneal degeneration (0.8%), superficial corneal ulcer (68.7%), deep corneal ulcer (8%) and perforating corneal ulcer (0.8%). The surgeries which consisted of the removal of the hair follicle was successful in 88.4% of the cases., Conclusion: EC is a rare condition which can be treated successfully by the removal of the hair follicles. It must be suspected in cases of corneal lesions unresponsive to medical treatment., (© 2021 American College of Veterinary Ophthalmologists.)