Recent Selection Research Articles

Decoding triancestral origins, archaic introgression, and natural selection in the Japanese population by whole-genome sequencing.

We generated Japanese Encyclopedia of Whole-Genome/Exome Sequencing Library (JEWEL), a high-depth whole-genome sequencing dataset comprising 3256 individuals from across Japan. Analysis of JEWEL revealed genetic characteristics of the Japanese population that were not discernible using microarray data. First, rare variant-based analysis revealed an unprecedented fine-scale genetic structure. Together with population genetics analysis, the present-day Japanese can be decomposed into three ancestral components. Second, we identified unreported loss-of-function (LoF) variants and observed that for specific genes, LoF variants appeared to be restricted to a more limited set of transcripts than would be expected by chance, with PTPRD as a notable example. Third, we identified 44 archaic segments linked to complex traits, including a Denisovan-derived segment at NKX6-1 associated with type 2 diabetes. Most of these segments are specific to East Asians. Fourth, we identified candidate genetic loci under recent natural selection. Overall, our work provided insights into genetic characteristics of the Japanese population.

Predicting Functional Consequences of Recent Natural Selection in Britain.

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.

Sex-biased gene regulation varies across human populations as a result of adaptive evolution.

Studies of human sexual dimorphism and gender disparities in health focus on ostensibly universal molecular sex differences, such as sex chromosomes and circulating hormone levels, while ignoring the extraordinary diversity in biology, behavior, and culture acquired by different human populations over their unique evolutionary histories. Using RNA-Seq data and whole genome sequences from 1000G and HGDP, we investigate variation in sex-biased gene expression across 11 human populations and test whether population-level variation in sex-biased expression may have resulted from adaptive evolution in regions containing sex-specific regulatory variants. We find that sex-biased gene expression in humans is highly variable, mostly population-specific, and demonstrates between population reversals. Expression quantitative trait locus mapping reveals sex-specific regulatory regions with evidence of recent positive natural selection, suggesting that variation in sex-biased expression may have evolved as an adaptive response to ancestral environments experienced by human populations. These results indicate that sex-biased gene expression is more flexible than previously thought and is not generally shared among human populations. Instead, molecular phenotypes associated with sex depend on complex interactions between population-specific molecular evolution and physiological responses to contemporary socioecologies.

Genomes from historical Drosophila melanogaster specimens illuminate adaptive and demographic changes across more than 200 years of evolution.

The ability to perform genomic sequencing on long-dead organisms is opening new frontiers in evolutionary research. These opportunities are especially notable in the case of museum collections, from which countless documented specimens may now be suitable for genomic analysis-if data of sufficient quality can be obtained. Here, we report 25 newly sequenced genomes from museum specimens of the model organism Drosophila melanogaster, including the oldest extant specimens of this species. By comparing historical samples ranging from the early 1800s to 1933 against modern-day genomes, we document evolution across thousands of generations, including time periods that encompass the species' initial occupation of northern Europe and an era of rapidly increasing human activity. We also find that the Lund, Sweden population underwent local genetic differentiation during the early 1800s to 1933 interval (potentially due to drift in a small population) but then became more similar to other European populations thereafter (potentially due to increased migration). Within each century-scale time period, our temporal sampling allows us to document compelling candidates for recent natural selection. In some cases, we gain insights regarding previously implicated selection candidates, such as ChKov1, for which our inferred timing of selection favors the hypothesis of antiviral resistance over insecticide resistance. Other candidates are novel, such as the circadian-related gene Ahcy, which yields a selection signal that rivals that of the DDT resistance gene Cyp6g1. These insights deepen our understanding of recent evolution in a model system, and highlight the potential of future museomic studies.

Recent natural selection conferred protection against schizophrenia by non-antagonistic pleiotropy

Schizophrenia is a debilitating psychiatric disorder associated with a reduced fertility and decreased life expectancy, yet common predisposing variation substantially contributes to the onset of the disorder, which poses an evolutionary paradox. Previous research has suggested balanced selection, a mechanism by which schizophrenia risk alleles could also provide advantages under certain environments, as a reliable explanation. However, recent studies have shown strong evidence against a positive selection of predisposing loci. Furthermore, evolutionary pressures on schizophrenia risk alleles could have changed throughout human history as new environments emerged. Here in this study, we used 1000 Genomes Project data to explore the relationship between schizophrenia predisposing loci and recent natural selection (RNS) signatures after the human diaspora out of Africa around 100,000 years ago on a genome-wide scale. We found evidence for significant enrichment of RNS markers in derived alleles arisen during human evolution conferring protection to schizophrenia. Moreover, both partitioned heritability and gene set enrichment analyses of mapped genes from schizophrenia predisposing loci subject to RNS revealed a lower involvement in brain and neuronal related functions compared to those not subject to RNS. Taken together, our results suggest non-antagonistic pleiotropy as a likely mechanism behind RNS that could explain the persistence of schizophrenia common predisposing variation in human populations due to its association to other non-psychiatric phenotypes.

GWAS on birth year infant mortality rates provides evidence of recent natural selection

Following more than a century of phenotypic measurement of natural selection processes, much recent work explores relationships between molecular genetic measurements and realized fitness in the next generation. We take an innovative approach to the study of contemporary selective pressure by examining which genetic variants are “sustained” in populations as mortality exposure increases. Specifically, we deploy a so-called “regional GWAS” (genome-wide association study) that links the infant mortality rate (IMR) by place and year in the United Kingdom with common genetic variants among birth cohorts in the UK Biobank. These cohorts (born between 1936 and 1970) saw a decline in IMR from above 65 to under 20 deaths per 1,000 live births, with substantial subnational variations and spikes alongside wartime exposures. Our results show several genome-wide significant loci, including LCT and TLR10/1/6, related to area-level cohort IMR exposure during gestation and infancy. Genetic correlations are found across multiple domains, including fertility, cognition, health behaviors, and health outcomes, suggesting an important role for cohort selection in modern populations.

Patterns of recent natural selection on genetic loci associated with sexually differentiated human body size and shape phenotypes.

Levels of sex differences for human body size and shape phenotypes are hypothesized to have adaptively reduced following the agricultural transition as part of an evolutionary response to relatively more equal divisions of labor and new technology adoption. In this study, we tested this hypothesis by studying genetic variants associated with five sexually differentiated human phenotypes: height, body mass, hip circumference, body fat percentage, and waist circumference. We first analyzed genome-wide association (GWAS) results for UK Biobank individuals (~194,000 females and ~167,000 males) to identify a total of 114,199 single nucleotide polymorphisms (SNPs) significantly associated with at least one of the studied phenotypes in females, males, or both sexes (P<5x10-8). From these loci we then identified 3,016 SNPs (2.6%) with significant differences in the strength of association between the female- and male-specific GWAS results at a low false-discovery rate (FDR<0.001). Genes with known roles in sexual differentiation are significantly enriched for co-localization with one or more of these SNPs versus SNPs associated with the phenotypes generally but not with sex differences (2.73-fold enrichment; permutation test; P = 0.0041). We also confirmed that the identified variants are disproportionately associated with greater phenotype effect sizes in the sex with the stronger association value. We then used the singleton density score statistic, which quantifies recent (within the last ~3,000 years; post-agriculture adoption in Britain) changes in the frequencies of alleles underlying polygenic traits, to identify a signature of recent positive selection on alleles associated with greater body fat percentage in females (permutation test; P = 0.0038; FDR = 0.0380), directionally opposite to that predicted by the sex differences reduction hypothesis. Otherwise, we found no evidence of positive selection for sex difference-associated alleles for any other trait. Overall, our results challenge the longstanding hypothesis that sex differences adaptively decreased following subsistence transitions from hunting and gathering to agriculture.

Molecular epidemiology and population genomics of Plasmodium knowlesi.

Molecular epidemiology has been central to uncovering P. knowlesi as an important cause of human malaria in Southeast Asia, and to understanding the complex nature of this zoonosis. Species-specific parasite detection and characterization of sequences were vital to show that P. knowlesi was distinct from the human parasite species that had been presumed to cause all malaria. With established sensitive and specific molecular detection tools, surveys subsequently indicated the distribution of P. knowlesi infections in humans, wild primate reservoir host species, and mosquito vector species. The importance of studying P. knowlesi genetic polymorphism was indicated initially by analysing a few nuclear gene loci as well as the mitochondrial genome, and subsequently by multi-locus microsatellite analyses and whole-genome sequencing. Different human infections generally have unrelated P. knowlesi genotypes, acquired from the diverse local parasite reservoirs in macaques. However, individual human infections are usually less genetically complex than those of wild macaques which experience more frequent superinfection with different P. knowlesi genotypes. Multi-locus analyses have revealed deep population subdivisions within P. knowlesi, which are structured both geographically and in relation to different macaque reservoir host species. Simplified genotypic discrimination assays now enable efficient large-scale surveillance of the sympatric P. knowlesi subpopulations within Malaysian Borneo. The whole-genome sequence analyses have also identified loci under recent positive natural selection in the P. knowlesi genome, with evidence that different loci are affected in different populations. These provide a foundation to understand recent adaptation of the zoonotic parasite populations, and to track and interpret future changes as they emerge.

Demographic history and selection at HLA loci in Native Americans.

The American continent was the last to be occupied by modern humans, and native populations bear the marks of recent expansions, bottlenecks, natural selection, and population substructure. Here we investigate how this demographic history has shaped genetic variation at the strongly selected HLA loci. In order to disentangle the relative contributions of selection and demography process, we assembled a dataset with genome-wide microsatellites and HLA-A, -B, -C, and -DRB1 typing data for a set of 424 Native American individuals. We find that demographic history explains a sizeable fraction of HLA variation, both within and among populations. A striking feature of HLA variation in the Americas is the existence of alleles which are present in the continent but either absent or very rare elsewhere in the world. We show that this feature is consistent with demographic history (i.e., the combination of changes in population size associated with bottlenecks and subsequent population expansions). However, signatures of selection at HLA loci are still visible, with significant evidence selection at deeper timescales for most loci and populations, as well as population differentiation at HLA loci exceeding that seen at neutral markers.

Probabilistic Estimation of Identity by Descent Segment Endpoints and Detection of Recent Selection

SummaryMost methods for fast detection of identity by descent (IBD) segments report identity by state segments without any quantification of the uncertainty in the endpoints and lengths of the IBD segments. We present a method for determining the posterior probability distribution of IBD segment endpoints. Our approach accounts for genotype errors, recent mutations, and gene conversions which disrupt DNA sequence identity within IBD segments, and it can be applied to large cohorts with whole-genome sequence or SNP array data. We find that our method’s estimates of uncertainty are well calibrated for homogeneous samples. We quantify endpoint uncertainty for 77.7 billion IBD segments from 408,883 individuals of white British ancestry in the UK Biobank, and we use these IBD segments to find regions showing evidence of recent natural selection. We show that many spurious selection signals are eliminated by the use of unbiased estimates of IBD segment endpoints and a pedigree-based genetic map. Eleven of the twelve regions with the greatest evidence for recent selection in our scan have been identified as selected in previous analyses using different approaches. Our computationally efficient method for quantifying IBD segment endpoint uncertainty is implemented in the open source ibd-ends software package.

Assortative mating at loci under recent natural selection in humans

Genetic correlation between mates at specific loci can greatly alter the evolutionary trajectory of a species. Genetic assortative mating has been documented in humans, but its existence beyond population stratification (shared ancestry) has been a matter of controversy. Here, we develop a method to measure assortative mating across the genome at 1,044,854 single-nucleotide polymorphisms (SNPs), controlling for population stratification and cohort-specific cryptic relatedness. Using data on 1683 human couples from two data sources, we find evidence for both assortative and disassortative mating at specific, discernible loci throughout the entire genome. Then, using the composite of multiple signals (CMS) score, we also show that the group of SNPs exhibiting the most assortativity has been under stronger recent positive selection. Simulations using realistic inputs confirm that assortative mating might indeed affect changes in allele frequency over time. These results suggest that genetic assortative mating may be speeding up evolution in humans.

Fine-Scale Inference of Ancestry Segments Without Prior Knowledge of Admixing Groups.

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

Supplemental Material for Salter-Townshend and Myers, 2019

Supplementary Material: Fine-scale Inference of Ancestry Segments without Prior Knowledge of Admixing Groups. We present an algorithm for inferring ancestry segments and characterizing admixture events, which involve an arbitrary number of genetically differentiated groups coming together. This allows inference of the demographic history of the species, properties of admixing groups, identification of signatures of natural selection, and may aid disease gene mapping. The algorithm employs nested hidden Markov models to obtain local ancestry estimation along the genome for each admixed individual. In a range of simulations, the accuracy of these estimates equals or exceeds leading existing methods. Moreover, and unlike these approaches, we do not require any prior knowledge of the relationship between sub-groups of donor reference haplotypes and the unseen mixing ancestral populations. Our approach infers these in terms of conditional copying probabilities. In application to the Human Genome Diversity Project we corroborate many previously inferred admixture events (e.g. an ancient admixture event in the Kalash). We further identify novel events such as complex 4-way admixture in San-Khomani individuals, and show that Eastern European populations possess 1 − 3% ancestry from a group resembling modern-day central Asians. We also identify evidence of recent natural selection favouring sub-Saharan ancestry at the HLA region, across North African individuals. We make available an R and C++ software library, which we term MOSAIC (which stands for MOSAIC Organises Segments of Ancestry In Chromosomes).

McSwan: A joint site frequency spectrum method to detect and date selective sweeps across multiple population genomes.

Inferring the mode and tempo of natural selection helps further our understanding of adaptation to past environmental changes. Here, we introduce McSwan, a method to detect and date past and recent natural selection events in the case of a hard sweep. The method is based on the comparison of site frequency spectra obtained under various demographic models that include selection. McSwan demonstrated high power (high sensitivity and specificity) in capturing hard selective sweep events without requiring haplotype phasing. It performed slightly better than SweeD when the recent effective population size was low and the genomic region was small. We then applied our method to a European (CEU) and an African (LWK) human re-sequencing data set. Most hard sweeps were detected in the CEU population (96%). Moreover, hard sweeps in the African population were estimated to have occurred further back in time (mode: 43,625years BP) compared to those of Europeans (mode: 24,850years BP). Most of the estimated ages of hard sweeps in Europeans were associated with the Last Glacial Maximum and were enriched in immunity-associated genes.

Analysis of type 2 diabetes and obesity genetic variants in Mexican Pima Indians: Marked allelic differentiation among Amerindians at HLA.

Prevalence of diabetes and obesity in Mexican Pima Indians is low, while prevalence in US Pima Indians is high. Although lifestyle likely accounts for much of the difference, the role of genetic factors is not well explored. To examine this, we genotyped 359 single nucleotide polymorphisms, including established type 2 diabetes and obesity variants from genome-wide association studies (GWAS) and 96 random markers, in 342 Mexican Pimas. A multimarker risk score of obesity variants was associated with body mass index (BMI; β=0.81kg/m2 per SD, P=0.0066). The mean value of the score was lower in Mexican Pimas than in US Pimas (P=4.3 × 10-11 ), and differences in allele frequencies at established loci could account for approximately 7% of the population difference in BMI; however, the difference in risk scores was consistent with evolutionary neutrality given genetic distance. To identify loci potentially under recent natural selection, allele frequencies at 283 variants were compared between US and Mexican Pimas, accounting for genetic distance. The largest differences were seen at HLA markers (e.g., rs9271720, difference=0.75, P=8.7 × 10-9 ); genetic distances at HLA were greater than at random markers (P=1.6 × 10-46 ). Analyses of GWAS data in 937 US Pimas also showed sharing of alleles identical by descent at HLA that exceeds its genomic expectation (P=7.0 × 10-10 ). These results suggest that, in addition to the widely recognized balancing selection at HLA, recent directional selection may also occur, resulting in marked allelic differentiation between closely related populations.

Searching for signals of recent natural selection in genes of the innate immune response - ancient DNA study

The last decade has seen sharp progress in the field of human evolutionary genetics and a great amount of genetic evidence of natural selection has been provided so far. Since host-pathogen co-evolution is difficult to trace due to the polygenic nature of human susceptibility to microbial diseases, of particular interest is any signal of natural selection in response to the strong selective pressure exerted by pathogens. Analysis of ancient DNA allows for the direct insight into changes of a gene pool content over time and enables monitoring allele frequency fluctuations. Among pathogenic agents, mycobacteria are proved to have remained in an intimate, long-lasting relation with humans, reflected by the current high level of host resistance. Therefore, we aimed to investigate the prevalence of several polymorphisms within innate immune response genes related to susceptibility to mycobacterial diseases (in SLC11A1, MBL2, TLR2, P2RX7, IL10, TNFA) in time series data from North and East Poland (1st–18th century AD, n = 207). The comparison of allele frequencies over time revealed a predominant role of genetic drift in shaping past gene pool of small, probably isolated groups, which was explained by the high level of population differentiation and limited gene flow. However, the trajectory of frequency fluctuations of two SNPs suggested the possibility of their non-neutral evolution and the results of applied forward simulations further strengthened the hypothesis of natural selection acting on those loci. However, we observed an unusual excess of homozygosity in the profile of several SNPs, which pinpoints to the necessity of further research on temporally and spatially diverse samples to support our inference on non-stochastic evolution, ideally employing pathway-based approaches. Nevertheless, our study confirms that time series data could help to decipher very recent human adaptation to life-threatening pathogens and assisting demographic events.

Recent natural selection causes adaptive evolution of an avian polygenic trait.

We used extensive data from a long-term study of great tits (Parus major) in the United Kingdom and Netherlands to better understand how genetic signatures of selection translate into variation in fitness and phenotypes. We found that genomic regions under differential selection contained candidate genes for bill morphology and used genetic architecture analyses to confirm that these genes, especially the collagen gene COL4A5, explained variation in bill length. COL4A5 variation was associated with reproductive success, which, combined with spatiotemporal patterns of bill length, suggested ongoing selection for longer bills in the United Kingdom. Last, bill length and COL4A5 variation were associated with usage of feeders, suggesting that longer bills may have evolved in the United Kingdom as a response to supplementary feeding.

Recent natural selection in a wild bird

Evolution Many studies have found evidence of rapid evolution in response to environmental change. In most cases, there has been some suggestion of which traits might be most responsive ahead of time. Bosse et al. turn this approach on its head by using genomic regions with a signature of selection to identify traits that are changing. In great tits ( Parus major ) in the United Kingdom, genomic regions showing selection invariably contained genes associated with bill growth. Indeed, U.K. birds not only have longer bills, but these longer bills are associated with increased fitness. These changes likely reflect an increase in domestic garden bird feeders over the past several decades. Science , this issue p. [365][1] [1]: /lookup/doi/10.1126/science.aal3298

The highly pleiotropic gene SLC39A8 as an opportunity to gain insight into the molecular pathogenesis of schizophrenia.

There is a long way from the initial discovery of a genome-wide significant signal to mechanistic understanding of the association. Identification of the gene and causal polymorphism usually requires an extensive additional effort. The schizophrenia genome-wide significant locus at 4q24 may be a rare exception to this pattern. As discussed in this review, the association at this locus is most probably driven by a functional missense variant at the metal cations transporter SLC39A8. The variant, rs13107325, is almost exclusive of European populations and is one of the most pleiotropic variants of the genome, being associated at genome-wide significant level with several additional traits, such as body mass index, Crohn's disease, blood pressure related-traits, and serum levels of manganese, N-terminal pro-B-type natriuretic peptide and HDL-cholesterol. SLC39A8 seems to be subject to recent natural selection in Europeans. It is almost ubiquitously expressed and its physiological role is beginning to be elucidated, mainly in relation to immunity. This manuscript presents arguments in favor of the rs13107325 variant as the functional variant responsible for the association of this locus with schizophrenia, reviews the genetic associations with this gene, the evidences of natural selection on the gene, and the known aspects about its structure and physiological functions. Finally, some hypotheses about putative mechanisms for its association with schizophrenia are presented based on this knowledge, including impaired immunity/inflammation, interference with glutamatergic neurotransmission, homeostasis of essential metals in brain, such as iron, zinc or manganese, or neurotoxicity by heavy metals, such as cadmium or lead.

Pooled ecotype sequencing reveals candidate genetic mechanisms for adaptive differentiation and reproductive isolation.

The early stages of speciation are often characterized by the formation of partially reproductively isolated ecotypes, which evolve as a by-product of divergent selective forces that are endemic to different habitats. Identifying the genomic regions, genes and ultimately functional polymorphisms that are involved in the processes of ecotype formation is inherently challenging, as there are likely to be many different loci involved in the process. To localize candidate regions of the genome contributing to ecotype formation, we conducted whole-genome pooled sequencing (pool-seq) with 47 coastal perennial and 50 inland annual populations of the yellow monkeyflower, Mimulus guttatus. Coastal perennial and inland annual ecotypes of M.guttatus have previously been shown to be ecologically reproductively isolated and highly locally adapted to their respective habitats. Our pool-seq results found allelic differentiation between the ecotypes for two chromosomal inversions, suggesting that frequencies of inversion heterokaryotypes are strongly differentiated between the ecotypes. Further, there were elevated levels of nonsynonymous change across chromosomal inversions. Across the genome, we identified multiple strong candidate genes potentially driving the morphological, life history and salt tolerance differences between the ecotypes. Several candidate genes coincide with previously identified quantitative trait locus regions and also show a signature of recent natural selection. Overall, the results of our study add to growing support for a major role of chromosomal inversions in adaptation and speciation and provide new insights into the genetic mechanisms underlying classic plant ecotype adaptations to wet and dry habitats.