Detecting the Genetic Signature of Natural Selection in Human Populations: Models, Methods, and Data

Integration of selective sweeps across the sheep genome: understanding the relationship between production and adaptation traits

BACKGROUND

Livestock populations are under constant selective pressure for higher productivity levels for different selective purposes. This pressure results in the selection of animals with unique adaptive and production traits. The study of genomic regions associated with these unique characteristics has the potential to improve biological knowledge regarding the adaptive process and how it is connected to production levels and resilience, which is the ability of an animal to adapt to stress or an imbalance in homeostasis. Sheep is a species that has been subjected to several natural and artificial selective pressures during its history, resulting in a highly specialized species for production and adaptation to challenging environments. Here, the data from multiple studies that aim at mapping selective sweeps across the sheep genome associated with production and adaptation traits were integrated to identify confirmed selective sweeps (CSS).

RESULTS

In total, 37 studies were used to identify 518 CSS across the sheep genome, which were classified as production (147 prodCSS) and adaptation (219 adapCSS) CSS based on the frequency of each type of associated study. The genes within the CSS were associated with relevant biological processes for adaptation and production. For example, for adapCSS, the associated genes were related to the control of seasonality, circadian rhythm, and thermoregulation. On the other hand, genes associated with prodCSS were related to the control of feeding behaviour, reproduction, and cellular differentiation. In addition, genes harbouring both prodCSS and adapCSS showed an interesting association with lipid metabolism, suggesting a potential role of this process in the regulation of pleiotropic effects between these classes of traits.

CONCLUSIONS

The findings of this study contribute to a deeper understanding of the genetic link between productivity and adaptability in sheep breeds. This information may provide insights into the genetic mechanisms that underlie undesirable genetic correlations between these two groups of traits and pave the way for a better understanding of resilience as a positive ability to respond to environmental stressors, where the negative effects on production level are minimized.

Multilocus phylogeography, population genetics and niche evolution of Australian brown and black-tailed treecreepers (Aves: Climacteris)

The Carpentarian barrier across north-eastern Australia is a major biogeographic barrier and a generator of biodiversity within the Australian Monsoonal Tropics. Here we present a continent-wide analysis of mitochondrial (control region) and autosomal (14 anonymous loci) sequence and indel variation and niche modelling of brown and black-tailed treecreepers (Climacteris picumnus and Climacteris melanurus), a clade with a classic distribution on either side of the Carpentarian barrier. mtDNA control region sequences exhibited reciprocal monophyly and strong differentiation (Fst = 0.91), and revealed a signature of a recent selective sweep in C. picumnus. A variety of tests support an isolation-with-migration model of divergence, albeit with low levels of gene flow across the Carpentarian barrier and a divergence time between species of ~1.7–2.8 Mya. Palaeoecological niche models show that both range size as measured by available habitat and estimated historical population sizes of both species declined in the past ~600 kyr and that the area of interspecific range overlap was never historically large, perhaps decreasing opportunities for extensive gene flow. The relatively long divergence time and low opportunity for gene flow may have facilitated speciation more so than in other co-distributed bird taxa across the Australian Monsoonal Tropics.

Identifying signatures of natural selection in Indian populations

In this study, we present the results of a genome-wide scan for signatures of positive selection using data from four tribal groups (Kokana, Warli, Bhil, and Pawara) and two caste groups (Deshastha Brahmin and Kunbi Maratha) from West of the Maharashtra State In India, as well as two samples of South Asian ancestry from the 1KG project (Gujarati Indian from Houston, Texas and Indian Telugu from UK). We used an outlier approach based on different statistics, including PBS, xpEHH, iHS, CLR, Tajima’s D, as well as two recently developed methods: Graph-aware Retrieval of Selective Sweeps (GRoSS) and Ascertained Sequentially Markovian Coalescent (ASMC). In order to minimize the risk of false positives, we selected regions that are outliers in all the samples included in the study using more than one method. We identified putative selection signals in 107 regions encompassing 434 genes. Many of the regions overlap with only one gene. The signals observed using microarray-based data are very consistent with our analyses using high-coverage sequencing data, as well as those identified with a novel coalescence-based method (ASMC). Importantly, at least 24 of these genomic regions have been identified in previous selection scans in South Asian populations or in other population groups. Our study highlights genomic regions that may have played a role in the adaptation of anatomically modern humans to novel environmental conditions after the out of Africa migration.

Rare Alleles and Signatures of Selection on the Immunodominant Domains of Pfs230 and Pfs48/45 in Malaria Parasites From Western Kenya

Background: Malaria elimination and eradication efforts can be advanced by including transmission-blocking or reducing vaccines (TBVs) alongside existing interventions. Key transmission-blocking vaccine candidates, such as Pfs230 domain one and Pfs48/45 domain 3, should be genetically stable to avoid developing ineffective vaccines due to antigenic polymorphisms. We evaluated genetic polymorphism and temporal stability of Pfs230 domain one and Pfs48/45 domain three in Plasmodium falciparum parasites from western Kenya. Methods: Dry blood spots on filter paper were collected from febrile malaria patients reporting to community health facilities in endemic areas of Homa Bay and Kisumu Counties and an epidemic-prone area of Kisii County in 2018 and 2019. Plasmodium speciation was performed using eluted DNA and real-time PCR. Amplification of the target domains of the two Pfs genes was performed on P. falciparum positive samples. We sequenced Pfs230 domain one on 156 clinical isolates and Pfs48/45 domain three on 118 clinical isolates to infer the levels of genetic variability, signatures of selection, genetic diversity indices and perform other evolutionary analyses. Results: Pfs230 domain one had low nucleotide diversity (π = 0.15 × 10^-2) with slight variation per study site. Six polymorphic sites with nonsynonymous mutations and eight haplotypes were discovered. I539T was a novel variant, whereas G605S was nearing fixation. Pfs48/45 domain three had a low π (0.063 × 10^-2), high conservation index, and three segregating sites, resulting in nonsynonymous mutation and four haplotypes. Some loci of Pfs230 D1 were in positive or negative linkage disequilibrium, had negative or positive selection signatures, and others (1813, 1955) and (1813, 1983) had a history of recombination. Mutated loci pairs in Pfs48/45 domain three had negative linkage disequilibrium, and some had negative and positive Tajima's D values with no history of recombination events. Conclusion: The two transmission blocking vaccine candidates have low nucleotide diversity, a small number of zone-specific variants, high nucleotide conservation index, and high frequency of rare alleles. With the near fixation a polymorphic site and the proximity of mutated codons to antibody binding epitopes, it will be necessary to continue monitoring sequence modifications of these domains when designing TBVs that include Pfs230 and Pfs48/45 antigens.

Worldwide diversity, association potential, and natural selection in the superimposed taste genes, CD36 and GNAT3

CD36 and GNAT3 mediate taste responses, with CD36 acting as a lipid detector and GNAT3 acting as the α subunit of gustducin, a G protein governing sweet, savory, and bitter transduction. Strikingly, the genes encoding CD36 and GNAT3 are genomically superimposed, with CD36 completely encompassing GNAT3. To characterize genetic variation across the CD36-GNAT3 region, its implications for phenotypic diversity, and its recent evolution, we analyzed from ~2,500 worldwide subjects sequenced by the 1000 Genomes Project (1000GP). CD36-GNAT3 harbored extensive diversity including 8,688 single-nucleotide polymorphisms (SNPs), 414 indels, and other complex variants. Sliding window analyses revealed that nucleotide diversity and population differentiation across CD36-GNAT3 were consistent with genome-wide trends in the 1000GP (π = 0.10%, P = 0.64; FST = 9.0%, P = 0.57). In addition, functional predictions using SIFT and PolyPhen-2 identified 60 variants likely to alter protein function, and they were in weak linkage disequilibrium (r2 < 0.17), suggesting their effects are largely independent. However, the frequencies of predicted functional variants were low (P¯ = 0.0013), indicating their contributions to phenotypic variance on population scales are limited. Tests using Tajima's D statistic revealed that pressures from natural selection have been relaxed across most of CD36-GNAT3 during its recent history (0.39 < P < 0.67). However, CD36 exons showed signs of local adaptation consistent with prior reports (P < 0.035). Thus, CD36 and GNAT3 harbor numerous variants predicted to affect taste sensitivity, but most are rare and phenotypic variance on a population level is likely mediated by a small number of sites.

Genome-Wide Detection of Structural Variations Reveals New Regions Associated with Domestication in Small Ruminants

During domestication processes, changes in selective pressures induce multiple phenotypical, physiological, and behavioral changes in target species. The rise of next-generation sequencing has provided a chance to study the genetics bases of these changes, most of the time based on single nucleotide polymorphisms (SNPs). However, several studies have highlighted the impact of structural variations (SVs) on individual fitness, particularly in domestic species. We aimed at unraveling the role of SVs during the domestication and later improvement of small ruminants by analyzing whole-genome sequences of 40 domestic sheep and 11 of their close wild relatives (Ovis orientalis), and 40 goats and 18 of their close wild relatives (Capra aegagrus). Using a combination of detection tools, we called 45,796 SVs in Ovis and 15,047 SVs in Capra genomes, including insertions, deletions, inversions, copy number variations, and chromosomal translocations. Most of these SVs were previously unreported in small ruminants. 69 and 45 SVs in sheep and goats, respectively, were in genomic regions with neighboring SNPs highly differentiated between wilds and domestics (i.e., putatively related to domestication). Among them, 25 and 20 SVs were close to or overlapping with genes related to physiological and morpho-anatomical traits linked with productivity (e.g., size, meat or milk quality, wool color), reproduction, or immunity. Finally, several of the SVs differentiated between wilds and domestics would not have been detected by screening only the differentiation of SNPs surrounding them, highlighting the complementarity of SVs and SNPs based approaches to detect signatures of selection.

The genomics of rapid climatic adaptation and parallel evolution in North American house mice

Parallel changes in genotype and phenotype in response to similar selection pressures in different populations provide compelling evidence of adaptation. House mice (Mus musculus domesticus) have recently colonized North America and are found in a wide range of environments. Here we measure phenotypic and genotypic differentiation among house mice from five populations sampled across 21° of latitude in western North America, and we compare our results to a parallel latitudinal cline in eastern North America. First, we show that mice are genetically differentiated between transects, indicating that they have independently colonized similar environments in eastern and western North America. Next, we find genetically-based differences in body weight and nest building behavior between mice from the ends of the western transect which mirror differences seen in the eastern transect, demonstrating parallel phenotypic change. We then conduct genome-wide scans for selection and a genome-wide association study to identify targets of selection and candidate genes for body weight. We find some genomic signatures that are unique to each transect, indicating population-specific responses to selection. However, there is significant overlap between genes under selection in eastern and western house mouse transects, providing evidence of parallel genetic evolution in response to similar selection pressures across North America.

Rapid evolutionary changes in gene expression in response to climate fluctuations

There is now abundant evidence of rapid evolution in natural populations, but the genetic mechanisms of these changes remain unclear. One possible route to rapid evolution is through changes in the expression of genes that influence traits under selection. We examined contemporary evolutionary gene expression changes in plant populations responding to environmental fluctuations. We compared genome‐wide gene expression, using RNA‐seq, in two populations of Brassica rapa collected over four time points between 1997 and 2014, during which precipitation in southern California fluctuated dramatically and phenotypic and genotypic changes occurred. By combining transcriptome profiling with the resurrection approach, we directly examined evolutionary changes in gene expression over time. For both populations, we found a substantial number of differentially expressed genes between generations, indicating rapid evolution in the expression of many genes. Using existing gene annotations, we found that many changes occurred in genes involved in regulating stress responses and flowering time. These appeared related to the fluctuations in precipitation and were potentially adaptive. However, the evolutionary changes in gene expression differed across generations within and between populations, indicating largely independent evolutionary trajectories across populations and over time. Our study provides strong evidence for rapid evolution in gene expression, and indicates that changes in gene expression can be one mechanism of rapid evolutionary responses to selection episodes. This study also illustrates that combining resurrection studies with transcriptomics is a powerful approach for investigating evolutionary changes at the gene regulatory level, and will provide new insights into the genetic basis of contemporary evolution.

see also the Perspective by Emily B. Josephs.

Evaluation of genetic diversity, agronomic traits, and anthracnose resistance in the NPGS Sudan Sorghum Core collection

Background

The United States Department of Agriculture (USDA) National Plant Germplasm System (NPGS) sorghum core collection contains 3011 accessions randomly selected from 77 countries. Genomic and phenotypic characterization of this core collection is necessary to encourage and facilitate its utilization in breeding programs and to improve conservation efforts. In this study, we examined the genome sequences of 318 accessions belonging to the NPGS Sudan sorghum core set, and characterized their agronomic traits and anthracnose resistance response.

Results

We identified 183,144 single nucleotide polymorphisms (SNPs) located within or in proximity of 25,124 annotated genes using the genotyping-by-sequencing (GBS) approach. The core collection was genetically highly diverse, with an average pairwise genetic distance of 0.76 among accessions. Population structure and cluster analysis revealed five ancestral populations within the Sudan core set, with moderate to high level of genetic differentiation. In total, 171 accessions (54%) were assigned to one of these populations, which covered 96% of the total genomic variation. Genome scan based on Tajima’s D values revealed two populations under balancing selection. Phenotypic analysis showed differences in agronomic traits among the populations, suggesting that these populations belong to different ecogeographical regions. A total of 55 accessions were resistant to anthracnose; these accessions could represent multiple resistance sources. Genome-wide association study based on fixed and random model Circulating Probability (farmCPU) identified genomic regions associated with plant height, flowering time, panicle length and diameter, and anthracnose resistance response. Integrated analysis of the Sudan core set and sorghum association panel indicated that a large portion of the genetic variation in the Sudan core set might be present in breeding programs but remains unexploited within some clusters of accessions.

Conclusions

The NPGS Sudan core collection comprises genetically and phenotypically diverse germplasm with multiple anthracnose resistance sources. Population genomic analysis could be used to improve screening efforts and identify the most valuable germplasm for breeding programs. The new GBS data set generated in this study represents a novel genomic resource for plant breeders interested in mining the genetic diversity of the NPGS sorghum collection.

Multiple selective sweeps of ancient polymorphisms in and around LTα located in the MHC class III region on chromosome 6

BACKGROUND

Lymphotoxin-α (LTα), located in the Major Histocompatibility Complex (MHC) class III region on chromosome 6, encodes a cytotoxic protein that mediates a variety of antiviral responses among other biological functions. Furthermore, several genotypes at this gene have been implicated in the onset of a number of complex diseases, including myocardial infarction, autoimmunity, and various types of cancer. However, little is known about levels of nucleotide variation and linkage disequilibrium (LD) in and near LTα, which could also influence phenotypic variance. To address this gap in knowledge, we examined sequence variation across ~ 10 kilobases (kbs), encompassing LTα and the upstream region, in 2039 individuals from the 1000 Genomes Project originating from 21 global populations.

RESULTS

Here, we observed striking patterns of diversity, including an excess of intermediate-frequency alleles, the maintenance of multiple common haplotypes and a deep coalescence time for variation (dating > 1.0 million years ago), in global populations. While these results are generally consistent with a model of balancing selection, we also uncovered a signature of positive selection in the form of long-range LD on chromosomes with derived alleles primarily in Eurasian populations. To reconcile these findings, which appear to support different models of selection, we argue that selective sweeps (particularly, soft sweeps) of multiple derived alleles in and/or near LTα occurred in non-Africans after their ancestors left Africa. Furthermore, these targets of selection were predicted to alter transcription factor binding site affinity and protein stability, suggesting they play a role in gene function. Additionally, our data also showed that a subset of these functional adaptive variants are present in archaic hominin genomes.

CONCLUSIONS

Overall, this study identified candidate functional alleles in a biologically-relevant genomic region, and offers new insights into the evolutionary origins of these loci in modern human populations.

"Omics" data integration and functional analyses link Enoyl-CoA hydratase, short chain 1 to drug refractory dilated cardiomyopathy

Background

Large-scale “omics” datasets have not been leveraged and integrated with functional analyses to discover potential drivers of cardiomyopathy. This study addresses the knowledge gap.

Methods

We coupled RNA sequence (RNA-Seq) variant detection and transcriptome profiling with pathway analysis to model drug refractory dilated cardiomyopathy (drDCM) using the BaseSpace sequencing hub and Ingenuity Pathway Analysis. We used RNA-Seq case-control datasets (n = 6 cases, n = 4 controls), exome sequence familial DCM datasets (n = 3 Italians, n = 5 Italians, n = 5 Chinese), and controls from the HapMap project (n = 5 Caucasians, and n = 5 Asians) for disease modeling and putative mutation discovery. Variant replication datasets: n = 128 cases and n = 15 controls. Source of datasets: NCBI Sequence Read Archive. Statistics: Pairwise differential expression analyses to determine differentially expressed genes and t-tests to calculate p-values. We adjusted for false discovery rates and reported q-values. We used chi-square tests to assess independence among variables, the Fisher’s Exact Tests and overlap p-values for the pathways and p-scores to rank network.

Results

Data revealed that ECHS1(enoyl-CoA hydratase, short chain 1(log₂(foldchange) = 1.63329) hosts a mirtron, MIR3944 expressed in drDCM (FPKM = 5.2857) and not in controls (FPKM = 0). Has-miR3944-3p is a putative target of BAG1 (BCL2 associated athanogene 1(log₂(foldchange) = 1.31978) and has-miR3944-5p of ITGAV (integrin subunit alpha V(log₂(foldchange) = 1.46107) and RHOD (ras homolog family member D(log₂(foldchange) = 1.28851). There is an association between ECHS1:11 V/A(rs10466126) and drDCM (p = 0.02496). The interaction (p = 2.82E-07) between ECHS1:75 T/I(rs1049951) and ECHS1:rs10466126 is associated with drDCM (p < 2.2e-16). ECHS1:rs10466126 and ECHS1:rs1049951 are in linkage disequilibrium (D’ = 1). The interaction (p = 7.84E-08) between ECHS1:rs1049951 and the novel ECHS1:c.41insT variant is associated with drDCM (p < 2.2e-16). The interaction (p = 0.001096) between DBT (Dihydrolipoamide branched chain transacylase E2):384G/S(rs12021720) and ECHS1:rs10466126 is associated with drDCM (p < 2.2e-16). At the mRNA level, there is an association between ECHS1 (log₂(foldchange) = 1.63329; q = 0.013927) and DBT (log₂(foldchange) = 0.955072; q = 0.0368792) with drDCM. ECHS1 is involved in valine (−log (p = 3.39E00)), isoleucine degradation (p = 0.00457), fatty acid β-oxidation (−log(p) = 2.83E00), and drug metabolism:cytochrome P450 (z-score = 2.07985196) pathways. The mitochondria (−log(p) = 8.73E00), oxidative phosphorylation (−log(p) = 5.35E00) and TCA-cycle II (−log(p) = 2.70E00) are dysfunctional.

Conclusions

We introduce an integrative data strategy that considers the interplay between the DNA, mRNA, and associated pathways, which represents a possible diagnostic, prognostic, biomarker, and personalized treatment discovery approach in genomically heterogeneous diseases.

Electronic supplementary material

The online version of this article (10.1186/s12920-018-0439-6) contains supplementary material, which is available to authorized users.

Integration of genome wide association studies and whole genome sequencing provides novel insights into fat deposition in chicken

Excessive fat deposition is a negative factor for poultry production because it reduces feed efficiency, increases the cost of meat production and is a health concern for consumers. We genotyped 497 birds from a Brazilian F₂ Chicken Resource Population, using a high-density SNP array (600 K), to estimate the genomic heritability of fat deposition related traits and to identify genomic regions and positional candidate genes (PCGs) associated with these traits. Selection signature regions, haplotype blocks and SNP data from a previous whole genome sequencing study in the founders of this chicken F2 population were used to refine the list of PCGs and to identify potential causative SNPs. We obtained high genomic heritabilities (0.43-0.56) and identified 22 unique QTLs for abdominal fat and carcass fat content traits. These QTLs harbored 26 PCGs involved in biological processes such as fat cell differentiation, insulin and triglyceride levels, and lipid biosynthetic process. Three of these 26 PCGs were located within haplotype blocks there were associated with fat traits, five overlapped with selection signature regions, and 12 contained predicted deleterious variants. The identified QTLs, PCGs and potentially causative SNPs provide new insights into the genetic control of fat deposition and can lead to improved accuracy of selection to reduce excessive fat deposition in chickens.

Evolutionary lessons from California plant phylogeography

Phylogeography documents the spatial distribution of genetic lineages that result from demographic processes, such as population expansion, population contraction, and gene movement, shaped by climate fluctuations and the physical landscape. Because most phylogeographic studies have used neutral markers, the role of selection may have been undervalued. In this paper, we contend that plants provide a useful evolutionary lesson about the impact of selection on spatial patterns of neutral genetic variation, when the environment affects which individuals can colonize new sites, and on adaptive genetic variation, when environmental heterogeneity creates divergence at specific loci underlying local adaptation. Specifically, we discuss five characteristics found in plants that intensify the impact of selection: sessile growth form, high reproductive output, leptokurtic dispersal, isolation by environment, and the potential to evolve longevity. Collectively, these traits exacerbate the impact of environment on movement between populations and local selection pressures-both of which influence phylogeographic structure. We illustrate how these unique traits shape these processes with case studies of the California endemic oak, Quercus lobata, and the western North American lichen, Ramalina menziesii Obviously, the lessons we learn from plant traits are not unique to plants, but they highlight the need for future animal, plant, and microbe studies to incorporate its impact. Modern tools that generate genome-wide sequence data are now allowing us to decipher how evolutionary processes affect the spatial distribution of different kinds of genes and also to better model future spatial distribution of species in response to climate change.

Myosin XI is associated with fitness and adaptation to aridity in wild pearl millet

Phenotypic changes in plants can be observed along many environmental gradients and are determined by both environmental and genetic factors. The identification of alleles associated with phenotypic variations is a rapidly developing area of research. We studied the genetic basis of phenotypic variations in 11 populations of wild pearl millet (Pennisetum glaucum) on two North-South aridity gradients, one in Niger and one in Mali. Most of the 11 phenotypic traits assessed in a common garden experiment varied between the populations studied. Moreover, the size of the inflorescence, the number of flowers and aboveground dry mass co-varied positively with a decrease in rainfall. To decipher the genetic basis of these phenotypes, we used an association mapping strategy with a mixed model. We found two SNPs on the same myosin XI contig significantly associated with variations in the average number of flowers. Both the allele frequency of the two SNPs and the average number of flowers co-varied with the rainfall gradient on the two gradients. Interestingly, this gene was also a target of selection during domestication. The Myosin XI gene is thus a good candidate for fitness-related adaptation in wild populations.

A Meta-Assembly of Selection Signatures in Cattle

Since domestication, significant genetic improvement has been achieved for many traits of commercial importance in cattle, including adaptation, appearance and production. In response to such intense selection pressures, the bovine genome has undergone changes at the underlying regions of functional genetic variants, which are termed “selection signatures”. This article reviews 64 recent (2009–2015) investigations testing genomic diversity for departure from neutrality in worldwide cattle populations. In particular, we constructed a meta-assembly of 16,158 selection signatures for individual breeds and their archetype groups (European, African, Zebu and composite) from 56 genome-wide scans representing 70,743 animals of 90 pure and crossbred cattle breeds. Meta-selection-scores (MSS) were computed by combining published results at every given locus, within a sliding window span. MSS were adjusted for common samples across studies and were weighted for significance thresholds across and within studies. Published selection signatures show extensive coverage across the bovine genome, however, the meta-assembly provides a consensus profile of 263 genomic regions of which 141 were unique (113 were breed-specific) and 122 were shared across cattle archetypes. The most prominent peaks of MSS represent regions under selection across multiple populations and harboured genes of known major effects (coat color, polledness and muscle hypertrophy) and genes known to influence polygenic traits (stature, adaptation, feed efficiency, immunity, behaviour, reproduction, beef and dairy production). As the first meta-assembly of selection signatures, it offers novel insights about the hotspots of selective sweeps in the bovine genome, and this method could equally be applied to other species.

Ancient DNA reveals selection acting on genes associated with hypoxia response in pre-Columbian Peruvian Highlanders in the last 8500 years

Archaeological evidence shows that humans began living in the high altitude Andes approximately 12,000 years ago. Andean highlanders are known to have developed the most complex societies of pre-Columbian South America despite challenges to their health and reproductive success resulting from chronic exposure to hypoxia. While the physiological adaptations to this environmental stressor are well studied in contemporary Andean highlanders, the molecular evolutionary processes associated with such adaptations remain unclear. We aim to better understand how humans managed to demographically establish in this harsh environment by addressing a central question: did exposure to hypoxia drive adaptation via natural selection within Andean populations or did an existing phenotype –characterized by reduced susceptibility to hypoxic stress–enable human settlement of the Andes? We genotyped three variable loci within the NOS3 and EGLN1 genes previously associated with adaptation to high altitude in 150 ancient human DNA samples from Peruvian high altitude and coastal low altitude sites in a time frame between ~8500–560 BP. We compare the data of 109 successful samples to forward simulations of genetic drift with natural selection and find that selection, rather than drift, explains the gradual frequency changes observed in the highland populations for two of the three SNPs.

Loss and Gain of Natural Killer Cell Receptor Function in an African Hunter-Gatherer Population

Modulating natural killer cell functions in human immunity and reproduction are diverse interactions between the killer cell immunoglobulin-like receptors (KIR) of Natural Killer (NK) cells and HLA class I ligands on the surface of tissue cells. Dominant interactions are between KIR2DL1 and the C2 epitope of HLA-C and between KIR2DL2/3 and the C1 epitope of HLA-C. KhoeSan hunter-gatherers of Southern Africa represent the earliest population divergence known and are the most genetically diverse indigenous people, qualities reflected in their KIR and HLA genes. Of the ten KhoeSan KIR2DL1 alleles, KIR2DL1*022 and KIR2DL1*026 likely originated in the KhoeSan, and later were transmitted at low frequency to the neighboring Zulus through gene flow. These alleles arose by point mutation from other KhoeSan KIR2DL1 alleles that are more widespread globally. Mutation of KIR2DL1*001 gave rise to KIR2DL1*022, causing loss of C2 recognition and gain of C1 recognition. This makes KIR2DL1*022 a more avid and specific C1 receptor than any KIR2DL2/3 allotype. Mutation of KIR2DL1*012 gave rise to KIR2DL1*026, causing premature termination of translation at the end of the transmembrane domain. This makes KIR2DL1*026 a membrane-associated receptor that lacks both a cytoplasmic tail and signaling function. At higher frequencies than their parental allotypes, the combined effect of the KhoeSan-specific KIR2DL1*022 and KIR2DL1*026 is to reduce the frequency of strong inhibitory C2 receptors and increase the frequency of strong inhibitory C1 receptors. Because interaction of KIR2DL1 with C2 is associated with risk of pregnancy disorder, these functional changes are potentially advantageous. Whereas all other KhoeSan KIR2DL1 alleles are present on a wide diversity of centromeric KIR haplotypes, KIR2DL1*026 is present on a single KIR haplotype and KIR2DL1*022 is present on two very similar haplotypes. The high linkage disequilibrium across their haplotypes is consistent with a recent emergence for these KIR2DL1 alleles that have distinctive functions.

Composite Selection Signals for Complex Traits Exemplified Through Bovine Stature Using Multibreed Cohorts of European and African Bos taurus

Understanding the evolution and molecular architecture of complex traits is important in domestic animals. Due to phenotypic selection, genomic regions develop unique patterns of genetic diversity called signatures of selection, which are challenging to detect, especially for complex polygenic traits. In this study, we applied the composite selection signals (CSS) method to investigate evidence of positive selection in a complex polygenic trait by examining stature in phenotypically diverse cattle comprising 47 European and 8 African Bos taurus breeds, utilizing a panel of 38,033 SNPs genotyped on 1106 animals. CSS were computed for phenotypic contrasts between multibreed cohorts of cattle by classifying the breeds according to their documented wither height to detect the candidate regions under selection. Using the CSS method, clusters of signatures of selection were detected at 26 regions (9 in European and 17 in African cohorts) on 13 bovine autosomes. Using comparative mapping information on human height, 30 candidate genes mapped at 12 selection regions (on 8 autosomes) could be linked to bovine stature diversity. Of these 12 candidate gene regions, three contained known genes (i.e., NCAPG-LCORL, FBP2-PTCH1, and PLAG1-CHCHD7) related to bovine stature, and nine were not previously described in cattle (five in European and four in African cohorts). Overall, this study demonstrates the utility of CSS coupled with strategies of combining multibreed datasets in the identification and discovery of genomic regions underlying complex traits. Characterization of multiple signatures of selection and their underlying candidate genes will elucidate the polygenic nature of stature across cattle breeds.

Adaptation to High Altitude: Phenotypes and Genotypes

Populations residing for millennia on the high-altitude plateaus of the world started natural experiments that we can evaluate to address questions about the processes of evolution and adaptation. A 2001 assessment in this journal summarized abundant evidence that Tibetan and Andean high-altitude natives had different phenotypes, and the article made a case for the hypothesis that different genetic bases underlie traits in the two populations. Since then, knowledge of the prehistory of high-altitude populations has grown, information about East African highlanders has become available, genomic science has grown exponentially, and the genetic and molecular bases of oxygen homeostasis have been clarified. Those scientific advances have transformed the study of high-altitude populations. The present review aims to summarize recent advances in understanding with an emphasis on the genetic bases of adaptive phenotypes, particularly hemoglobin concentration among Tibetan highlanders. EGLN1 and EPAS1 encode two crucial proteins contributing to oxygen homeostasis, the oxygen sensor PHD2 and the transcription factor subunit HIF-2α, respectively; they show signals of natural selection such as marked allele frequency differentiation between Tibetans and lowland populations. EPAS1 genotypes associated in several studies with the dampened hemoglobin phenotype that is characteristic of Tibetans at high altitude but did not associate with the dampened response among Amhara from Ethiopia or the vigorous elevation of hemoglobin concentration among Andean highlanders. Future work will likely develop understanding of the integrative biology leading from genotype to phenotype to population in all highland areas.

Detecting recent positive selection with high accuracy and reliability by conditional coalescent tree

Studies of natural selection, followed by functional validation, are shedding light on understanding of genetic mechanisms underlying human evolution and adaptation. Classic methods for detecting selection, such as the integrated haplotype score (iHS) and Fay and Wu's H statistic, are useful for candidate gene searching underlying positive selection. These methods, however, have limited capability to localize causal variants in selection target regions. In this study, we developed a novel method based on conditional coalescent tree to detect recent positive selection by counting unbalanced mutations on coalescent gene genealogies. Extensive simulation studies revealed that our method is more robust than many other approaches against biases due to various demographic effects, including population bottleneck, expansion, or stratification, while not sacrificing its power. Furthermore, our method demonstrated its superiority in localizing causal variants from massive linked genetic variants. The rate of successful localization was about 20-40% higher than that of other state-of-the-art methods on simulated data sets. On empirical data, validated functional causal variants of four well-known positive selected genes were all successfully localized by our method, such as ADH1B, MCM6, APOL1, and HBB. Finally, the computational efficiency of this new method was much higher than that of iHS implementations, that is, 24-66 times faster than the REHH package, and more than 10,000 times faster than the original iHS implementation. These magnitudes make our method suitable for applying on large sequencing data sets. Software can be downloaded from https://github.com/wavefancy/scct.

Signatures of volcanism and aridity in the evolution of an insular pine (Pinus canariensis Chr. Sm. Ex DC in Buch)

Oceanic islands of volcanic origin provide useful templates for the study of evolution because they are subjected to recurrent perturbations that generate steep environmental gradients that may promote adaptation. Here we combine population genetic data from nuclear genes with the analysis of environmental variation and phenotypic data from common gardens to disentangle the confounding effects of demography and selection to identify the factors of importance for the evolution of the insular pine P. canariensis. Eight nuclear genes were partially sequenced in a survey covering the entire species range, and phenotypic traits were measured in four common gardens from contrasting environments. The explanatory power of population substrate age and environmental indices were assessed against molecular and phenotypic diversity estimates. In addition, neutral genetic variability (FST) and the genetic differentiation of phenotypic variation (QST) were compared in order to identify the evolutionary forces acting on these traits. Two key factors in the evolution of the species were identified: (1) recurrent volcanic activity has left an imprint in the genetic diversity of the nuclear genes; (2) aridity in southern slopes promotes local adaptation in the driest localities of P. canariensis, despite high levels of gene flow among populations.

Opportunity for selection in human health

Natural selection defined by differential survival and reproduction of individuals in populations is influenced by genetic, developmental, and environmental factors operating at every age and stage in human life history: generation of gametes, conception, birth, maturation, reproduction, senescence, and death. Biological systems are built upon a hierarchical organization nesting subcellular organelles, cells, tissues, and organs within individuals, individuals within families, and families within populations, and the latter among other populations. Natural selection often acts simultaneously at more than one level of biological organization and on specific traits, which we define as multilevel selection. Under this model, the individual is a fundamental unit of biological organization and also of selection, imbedded in a larger evolutionary context, just as it is a unit of medical intervention imbedded in larger biological, cultural, and environmental contexts. Here, we view human health and life span as necessary consequences of natural selection, operating at all levels and phases of biological hierarchy in human life history as well as in sociological and environmental milieu. An understanding of the spectrum of opportunities for natural selection will help us develop novel approaches to improving healthy life span through specific and global interventions that simultaneously focus on multiple levels of biological organization. Indeed, many opportunities exist to apply multilevel selection models employed in evolutionary biology and biodemography to improving human health at all hierarchical levels. Multilevel selection perspective provides a rational theoretical foundation for a synthesis of medicine and evolution that could lead to discovering effective predictive, preventive, palliative, potentially curative, and individualized approaches in medicine and in global health programs.

An ABC method for estimating the rate and distribution of effects of beneficial mutations

Determining the distribution of adaptive mutations available to natural selection is a difficult task. These are rare events and most of them are lost by chance. Some theoretical works propose that the distribution of newly arising beneficial mutations should be close to exponential. Empirical data are scarce and do not always support an exponential distribution. Analysis of the dynamics of adaptation in asexual populations of microorganisms has revealed that these can be summarized by two effective parameters, the effective mutation rate, U_e, and the effective selection coefficient of a beneficial mutation, S_e. Here, we show that these effective parameters will not always reflect the rate and mean effect of beneficial mutations, especially when the distribution of arising mutations has high variance, and the mutation rate is high. We propose a method to estimate the distribution of arising beneficial mutations, which is motivated by a common experimental setup. The method, which we call One Biallelic Marker Approximate Bayesian Computation, makes use of experimental data consisting of periodic measures of neutral marker frequencies and mean population fitness. Using simulations, we find that this method allows the discrimination of the shape of the distribution of arising mutations and that it provides reasonable estimates of their rates and mean effects in ranges of the parameter space that may be of biological relevance.

Detecting signatures of positive selection in partial sequences generated on a large scale: pitfalls, procedures and resources

Studying the actions of selection provides insight into adaptation, population divergence and gene function. Next-generation sequencing produces large amounts of partial sequences, potentially facilitating efforts to detect signatures of selection based on comparisons between synonymous (d(S)) and nonsynonymous (d(N)) substitutions, and single nucleotide polymorphism assays placed in selected genes would improve the ability to study adaptation in population surveys. However, sequences generated by these technologies are typically short. In nonmodel organisms that are a focus of evolutionary studies, the lack of a reference genome that facilitates the assembly of short sequences has limited surveys of positive selection in large numbers of genes. Here, we describe a series of steps to facilitate these surveys. We provide PERL scripts to assist data analysis, and describe the use of commonly available programs. We demonstrate these approaches in six salmon species, which have partially duplicated genomes. We recommend using multiway blast to optimize the number of alignments between partial coding sequences. Reading frames should be manually detected after alignment with sequences in Genbank using the BLASTX program. We encourage the use of a phylogenetic approach to separate orthologs from paralogs in duplicated genomes. Simple simulations on a gene known to have undergone selection in salmon species, transferrin, showed that the ability to detect selection in short sequences (<600 bp) depended on the proportion of codons under selection (1-2%) within that sequence. This relationship was less relevant in longer sequences. In this exploratory study, we detected 11 genes showing evidence of positive selection.

Adaptations to new environments in humans: the role of subtle allele frequency shifts

Humans show tremendous phenotypic diversity across geographically distributed populations, and much of this diversity undoubtedly results from genetic adaptations to different environmental pressures. The availability of genome-wide genetic variation data from densely sampled populations offers unprecedented opportunities for identifying the loci responsible for these adaptations and for elucidating the genetic architecture of human adaptive traits. Several approaches have been used to detect signals of selection in human populations, and these approaches differ in the assumptions they make about the underlying mode of selection. We contrast the results of approaches based on haplotype structure and differentiation of allele frequencies to those from a method for identifying single nucleotide polymorphisms strongly correlated with environmental variables. Although the first group of approaches tends to detect new beneficial alleles that were driven to high frequencies by selection, the environmental correlation approach has power to identify alleles that experienced small shifts in frequency owing to selection. We suggest that the first group of approaches tends to identify only variants with relatively strong phenotypic effects, whereas the environmental correlation methods can detect variants that make smaller contributions to an adaptive trait.

The power to detect artificial selection acting on single loci in recently domesticated species

Background

An increasing number of aquaculture species are subjected to artificial selection in systematic breeding programs. Rapid improvements of important commercial traits are reported, but little is known about the effects of the strong directional selection applied, on gene level variation. Large numbers of genetic markers are becoming available, making it feasible to detect and estimate these effects. Here a simulation tool was developed in order to explore the power by which single genetic loci subjected to uni-directional selection in parallel breeding populations may be detected.

Findings

Two simulation models were pursued: 1) screening for loci displaying higher genetic differentiation than expected (high-F_ST outliers), from neutral evolution between a pool of domesticated populations and a pool of wild populations; 2) screening for loci displaying lower genetic differentiation (low-F_ST outliers) between domesticated strains than expected from neutral evolution. The premise for both approaches was that the isolated domesticated strains are subjected to the same breeding goals. The power to detect outlier loci was calculated under the following parameter values: number of populations, effective population size per population, number of generations since onset of selection, initial F_ST, and the selection coefficient acting on the locus. Among the parameters investigated, selection coefficient, the number of generation since onset of selection, and number of populations, had the largest impact on power. The power to detect loci subjected to directional in breeding programmes was high when applying the between farmed and wild population approach, and low for the between farmed populations approach.

Conclusions

A simulation tool was developed for estimating the power to detect artificial selection acting directly on single loci. The simulation tool should be applicable to most species subject to domestication, as long as a reasonable high accuracy in input parameters such as effective population size, number of generations since the initiation of selection, and initial differentiation (F_ST) can be obtained. Identification of genetic loci under artificial selection would be highly valuable, since such loci could be used to monitor maintenance of genetic variation in the breeding populations and monitoring possible genetic changes in wild populations from genetic interaction between escapees and their wild counterpart.

The evolution of human genetic and phenotypic variation in Africa

Africa is the birthplace of modern humans, and is the source of the geographic expansion of ancestral populations into other regions of the world. Indigenous Africans are characterized by high levels of genetic diversity within and between populations. The pattern of genetic variation in these populations has been shaped by demographic events occurring over the last 200,000 years. The dramatic variation in climate, diet, and exposure to infectious disease across the continent has also resulted in novel genetic and phenotypic adaptations in extant Africans. This review summarizes some recent advances in our understanding of the demographic history and selective pressures that have influenced levels and patterns of diversity in African populations.

Sodium, blood pressure, and ethnicity: what have we learned?

An enormous amount of research has yielded significant knowledge about ethnic differences in sodium homeostasis and blood pressure regulation. Consistent findings such as greater sodium-sensitivity, lower potassium excretion and high higher serum sodium levels in African Americans need further exploration to define more precise physiological mechanisms. The genetic alleles associated with sodium homeostasis in relation to blood pressure have accounted for only a small proportion of the variance in blood pressure. Several allelic variants differ in frequency among ethnic groups and heat-adapted genetic variants have a high prevalence in low latitudes and hot, wet climates which lends support to the "sodium retention" hypothesis. The blood pressure disparities between African Americans and whites may, in part, be due to different allelic frequencies of genes associated with sodium homeostasis. However, with advances in genomics, environmental factors tend to be neglected in research. Better measures of environmental stress have recently been developed by anthropologists and should be included in research designs by investigators in other disciplines. Public health efforts should encourage food producers to reduce sodium content of its products, and physicians should encourage patients to reduce consumption of high sodium packaged and fast foods.

Colloquium papers: Numbering the hairs on our heads: the shared challenge and promise of phenomics

Evolution and medicine share a dependence on the genotype-phenotype map. Although genotypes exist and are inherited in a discrete space convenient for many sorts of analyses, the causation of key phenomena such as natural selection and disease takes place in a continuous phenotype space whose relationship to the genotype space is only dimly grasped. Direct study of genotypes with minimal reference to phenotypes is clearly insufficient to elucidate these phenomena. Phenomics, the comprehensive study of phenotypes, is therefore essential to understanding biology. For all of the advances in knowledge that a genomic approach to biology has brought, awareness is growing that many phenotypes are highly polygenic and susceptible to genetic interactions. Prime examples are common human diseases. Phenomic thinking is starting to take hold and yield results that reveal why it is so critical. The dimensionality of phenotypic data are often extremely high, suggesting that attempts to characterize phenotypes with a few key measurements are unlikely to be completely successful. However, once phenotypic data are obtained, causation can turn out to be unexpectedly simple. Phenotypic data can be informative about the past history of selection and unexpectedly predictive of long-term evolution. Comprehensive efforts to increase the throughput and range of phenotyping are an urgent priority.