Genome-wide Selection Scans Research Articles

Signals of selection in outlier loci in a widely dispersing species across an environmental mosaic

Local adaptation reflects a balance between natural selection and gene flow and is classically thought to require the retention of locally adapted alleles. However, organisms with high dispersal potential across a spatially or temporally heterogeneous landscape pose an interesting challenge to this view requiring local selection every generation or when environmental conditions change to generate adaptation in adults. Here, we test for geographical and sequence-based signals of selection in five putatively adaptive and two putatively neutral genes identified in a previous genome scan of the highly dispersing purple sea urchin, Strongylocentrotus purpuratus. Comparing six populations spanning the species' wide latitudinal range from Canada to Baja California, Mexico, we find positive tests for selection in the putative adaptive genes and not in the putative neutral genes. Specifically, we find an excess of low-frequency and nonsynonymous polymorphisms in two transcription factors and a transporter protein, and an excess of common amino acid polymorphisms in the two transcription factors, suggestive of spatially balancing selection. We test for a genetic correlation with temperature, a dominant environmental variable in this coastal ecosystem. We find mild clines and a stronger association of genetic variation with temperature than latitude in four of the five putative adaptive loci and a signal of local adaptation in the Southern California Bight. Overall, patterns of genetic variation match predictions based on spatially or temporally balancing selection in a heterogeneous landscape and illustrate the value of geographical and coalescent tests on candidate loci identified in a genome-wide scan for selection.

How to Make a Dolphin: Molecular Signature of Positive Selection in Cetacean Genome.

Cetaceans are unique in being the only mammals completely adapted to an aquatic environment. This adaptation has required complex changes and sometimes a complete restructuring of physiology, behavior and morphology. Identifying genes that have been subjected to selection pressure during cetacean evolution would greatly enhance our knowledge of the ways in which genetic variation in this mammalian order has been shaped by natural selection. Here, we performed a genome-wide scan for positive selection in the dolphin lineage. We employed models of codon substitution that account for variation of selective pressure over branches on the tree and across sites in a sequence. We analyzed 7,859 nuclear-coding ortholog genes and using a series of likelihood ratio tests (LRTs), we identified 376 genes (4.8%) with molecular signatures of positive selection in the dolphin lineage. We used the cow as the sister group and compared estimates of selection in the cetacean genome to this using the same methods. This allowed us to define which genes have been exclusively under positive selection in the dolphin lineage. The enrichment analysis found that the identified positively selected genes are significantly over-represented for three exclusive functional categories only in the dolphin lineage: segment specification, mesoderm development and system development. Of particular interest for cetacean adaptation to an aquatic life are the following GeneOntology targets under positive selection: genes related to kidney, heart, lung, eye, ear and nervous system development.

Multiple instances of ancient balancing selection shared between humans and chimpanzees.

Instances in which natural selection maintains genetic variation in a population over millions of years are thought to be extremely rare. We conducted a genome-wide scan for long-lived balancing selection by looking for combinations of SNPs shared between humans and chimpanzees. In addition to the major histocompatibility complex, we identified 125 regions in which the same haplotypes are segregating in the two species, all but two of which are noncoding. In six cases, there is evidence for an ancestral polymorphism that persisted to the present in humans and chimpanzees. Regions with shared haplotypes are significantly enriched for membrane glycoproteins, and a similar trend is seen among shared coding polymorphisms. These findings indicate that ancient balancing selection has shaped human variation and point to genes involved in host-pathogen interactions as common targets.

Improving Genome-Wide Scans of Positive Selection by Using Protein Isoforms of Similar Length

Large-scale evolutionary studies often require the automated construction of alignments of a large number of homologous gene families. The majority of eukaryotic genes can produce different transcripts due to alternative splicing or transcription initiation, and many such transcripts encode different protein isoforms. As analyses tend to be gene centered, one single-protein isoform per gene is selected for the alignment, with the de facto approach being to use the longest protein isoform per gene (Longest), presumably to avoid including partial sequences and to maximize sequence information. Here, we show that this approach is problematic because it increases the number of indels in the alignments due to the inclusion of nonhomologous regions, such as those derived from species-specific exons, increasing the number of misaligned positions. With the aim of ameliorating this problem, we have developed a novel heuristic, Protein ALignment Optimizer (PALO), which, for each gene family, selects the combination of protein isoforms that are most similar in length. We examine several evolutionary parameters inferred from alignments in which the only difference is the method used to select the protein isoform combination: Longest, PALO, the combination that results in the highest sequence conservation, and a randomly selected combination. We observe that Longest tends to overestimate both nonsynonymous and synonymous substitution rates when compared with PALO, which is most likely due to an excess of misaligned positions. The estimation of the fraction of genes that have experienced positive selection by maximum likelihood is very sensitive to the method of isoform selection employed, both when alignments are constructed with MAFFT and with Prank+F. Longest performs better than a random combination but still estimates up to 3 times more positively selected genes than the combination showing the highest conservation, indicating the presence of many false positives. We show that PALO can eliminate the majority of such false positives and thus that it is a more appropriate approach for large-scale analyses than Longest. A web server has been set up to facilitate the use of PALO given a user-defined set of gene families; it is available at http://evolutionarygenomics.imim.es/palo.

Studies of Natural Selection in the Era of Next-generation Sequencing

A deeper understanding of positive selection is of fundamental importance. Evolutionary pressures acting on the genome as a whole and the specific role of environmental pressures on divergence shape the natural selection. However, despite the intense interest in genome-wide scans, a coherent study of recent human evolutionary history has yet to emerge due to the limitations of data, models, and analytical tools. Fast and cheaper Next Generation Sequencing (NGS) technologies will generate unprecedentedly massive and highly-dimensional genomic variation data. NGS will offer unprecedented opportunities for population genetic and natural selection studies, but also raise great challenges. In the conventional population genetics, most researches have primarily focused on natural selection acting on a single locus. Little attention has been paid to determining how the natural selection acts on multiple interacted genes in response to environmental perturbation. In addition, the current paradigm for analysis of natural selection on gene-expression is to use a single value of summarizing statistic to represent gene expression level and overlook all information on expression difference in exons, genomic position and alleles. To address these critical limitations, we develop a unified framework and statistical methods for genome-wide scans for natural selection and investigation of natural selection on pathways. To explore observed expression variation in exons or in genomic position across the genes, we extend one dimensional diffusion process to multidimensional diffusion process and a single variate stochastic differential equation to multivariate stochastic differential equations. We then use the extended multidimensional diffusion processes to model the evolution of gene-expression acted by natural selection. We hope that the present new development of natural selection analysis for NGS data will open a new avenue for natural selection analysis with NGS data.

LARGE-SCALE CANDIDATE GENE SCAN REVEALS THE ROLE OF CHEMORECEPTOR GENES IN HOST PLANT SPECIALIZATION AND SPECIATION IN THE PEA APHID

Understanding the drivers of speciation is critical to interpreting patterns of biodiversity. The identification of the genetic changes underlying adaptation and reproductive isolation is necessary to link barriers to gene flow to the causal origins of divergence. Here, we present a novel approach to the genetics of speciation, which should complement the commonly used approaches of quantitative trait locus mapping and genome-wide scans for selection. We present a large-scale candidate gene approach by means of sequence capture, applied to identifying the genetic changes underlying reproductive isolation in the pea aphid, a model system for the study of ecological speciation. Targeted resequencing enabled us to scale up the candidate gene approach, specifically testing for the role of chemosensory gene families in host plant specialization. Screening for the signature of divergence under selection at 172 candidate and noncandidate loci, we revealed a handful of loci that show high levels of differentiation among host races, which almost all correspond to odorant and gustatory receptor genes. This study offers the first indication that some chemoreceptor genes, often tightly linked together in the genome, could play a key role in local adaptation and reproductive isolation in the pea aphid and potentially other phytophagous insects. Our approach opens a new route toward the functional genomics of ecological speciation.

Rehh: an R package to detect footprints of selection in genome-wide SNP data from haplotype structure

With the development of next-generation sequencing and genotyping approaches, large single nucleotide polymorphism haplotype datasets are becoming available in a growing number of both model and non-model species. Identifying genomic regions with unexpectedly high local haplotype homozygosity relatively to neutral expectation represents a powerful strategy to ascertain candidate genes responding to natural or artificial selection. To facilitate genome-wide scans of selection based on the analysis of long-range haplotypes, we developed the R package rehh. It provides a versatile tool to detect the footprints of recent or ongoing selection with several graphical functions that help visual interpretation of the results. Stable version is available from CRAN: http://cran.r-project.org/. Development version is available from the R-forge repository: http://r-forge.r-project.org/projects/rehh. Both versions can be installed directly from R. Function documentation and example data files are provided within the package and a tutorial is available as Supplementary Material. rehh is distributed under the GNU General Public Licence (GPL ≥ 2).

Loss and Gain of Function in SERPINB11: An Example of a Gene under Selection on Standing Variation, with Implications for Host-Pathogen Interactions

Serine protease inhibitors (SERPINs) are crucial in the regulation of diverse biological processes including inflammation and immune response. SERPINB11, located in the 18q21 gene cluster, is a polymorphic gene/pseudogene coding for a non-inhibitory SERPIN. In a genome-wide scan for recent selection, SERPINB11 was identified as a potential candidate gene for adaptive evolution in Yoruba. The present study sought a better understanding of the evolutionary history of SERPINB11, with special focus on evaluating its selective signature. Through the resequencing of coding and noncoding regions of SERPINB11 in 20 Yorubans and analyzing primate orthologous sequences, we identified a full-length SERPINB11 variant encoding a non-inhibitory SERPIN as the putative candidate of selection – probably driven to higher frequencies by an adaptive response using preexisting variation. In addition, we detected contrasting evolutionary features of SERPINB11 in primates: While primate phylogeny as a whole is under purifying selection, the human lineage shows evidence of positive selection in a few codons, all associated with the active SERPINB11. Comparative modeling studies suggest that positively selected codons reduce SERPINB11's ability to undergo the conformational changes typical of inhibitory SERPINs – suggesting that it is evolving towards a new non-inhibitory function in humans. Significant correlations between SERPINB11 variants and the environmental variables, pastoralism and pathogen richness, have led us to propose a selective advantage through host-pathogen interactions, possibly linked to an adaptive response combating the emergence of infectious diseases in recent human evolution. This work represents the first description of a resurrected gene in humans, and may well exemplify selection on standing variation triggered by drastic ecological shifts.

Genetic adaptation to high altitude in the Ethiopian highlands

BackgroundGenomic analysis of high-altitude populations residing in the Andes and Tibet has revealed several candidate loci for involvement in high-altitude adaptation, a subset of which have also been shown to be associated with hemoglobin levels, including EPAS1, EGLN1, and PPARA, which play a role in the HIF-1 pathway. Here, we have extended this work to high- and low-altitude populations living in Ethiopia, for which we have measured hemoglobin levels. We genotyped the Illumina 1M SNP array and employed several genome-wide scans for selection and targeted association with hemoglobin levels to identify genes that play a role in adaptation to high altitude.ResultsWe have identified a set of candidate genes for positive selection in our high-altitude population sample, demonstrated significantly different hemoglobin levels between high- and low-altitude Ethiopians and have identified a subset of candidate genes for selection, several of which also show suggestive associations with hemoglobin levels.ConclusionsWe highlight several candidate genes for involvement in high-altitude adaptation in Ethiopia, including CBARA1, VAV3, ARNT2 and THRB. Although most of these genes have not been identified in previous studies of high-altitude Tibetan or Andean population samples, two of these genes (THRB and ARNT2) play a role in the HIF-1 pathway, a pathway implicated in previous work reported in Tibetan and Andean studies. These combined results suggest that adaptation to high altitude arose independently due to convergent evolution in high-altitude Amhara populations in Ethiopia.

Systematic underestimation of the age of selected alleles

A common interpretation of genome-wide selection scans is that the dispersal of anatomically modern humans out of Africa and into diverse environments led to a number of genetic adaptations. If so, patterns of polymorphism from non-African individuals should show the signature of adaptations dating to 40,000–100,000 Kya, coinciding with the main exodus from Africa. However, scans of polymorphism data from a few populations have yielded conflicting results about the chronology of local, population-specific adaptations. In particular, a number of papers report very recent ages for selected alleles in humans, which postdate the development of agriculture 10 Kya, and suggest that adaptive differences among human populations are much more recent. I present an analysis of simulations suggesting a downward bias in methods commonly used to estimate the age of selected alleles. These findings indicate that an estimate of a time to the most recent common ancestor (tMRCA) obtained using standard methods (used as a proxy for the age of an allele) of less than 10 Kya is consistent with an allele that actually became selected before the onset of agriculture and potentially as early as 50 Kya. These findings suggest that the genomic scans for selection may be consistent with selective pressures tied to the Out of Africa expansion of modern human populations.

Signatures of Environmental Genetic Adaptation Pinpoint Pathogens as the Main Selective Pressure through Human Evolution

Previous genome-wide scans of positive natural selection in humans have identified a number of non-neutrally evolving genes that play important roles in skin pigmentation, metabolism, or immune function. Recent studies have also shown that a genome-wide pattern of local adaptation can be detected by identifying correlations between patterns of allele frequencies and environmental variables. Despite these observations, the degree to which natural selection is primarily driven by adaptation to local environments, and the role of pathogens or other ecological factors as selective agents, is still under debate. To address this issue, we correlated the spatial allele frequency distribution of a large sample of SNPs from 55 distinct human populations to a set of environmental factors that describe local geographical features such as climate, diet regimes, and pathogen loads. In concordance with previous studies, we detected a significant enrichment of genic SNPs, and particularly non-synonymous SNPs associated with local adaptation. Furthermore, we show that the diversity of the local pathogenic environment is the predominant driver of local adaptation, and that climate, at least as measured here, only plays a relatively minor role. While background demography by far makes the strongest contribution in explaining the genetic variance among populations, we detected about 100 genes which show an unexpectedly strong correlation between allele frequencies and pathogenic environment, after correcting for demography. Conversely, for diet regimes and climatic conditions, no genes show a similar correlation between the environmental factor and allele frequencies. This result is validated using low-coverage sequencing data for multiple populations. Among the loci targeted by pathogen-driven selection, we found an enrichment of genes associated to autoimmune diseases, such as celiac disease, type 1 diabetes, and multiples sclerosis, which lends credence to the hypothesis that some susceptibility alleles for autoimmune diseases may be maintained in human population due to past selective processes.

Recent advances in understanding the role of nutrition in human genome evolution.

Dietary transitions in human history have been suggested to play important roles in the evolution of mankind. Genetic variations caused by adaptation to diet during human evolution could have important health consequences in current society. The advance of sequencing technologies and the rapid accumulation of genome information provide an unprecedented opportunity to comprehensively characterize genetic variations in human populations and unravel the genetic basis of human evolution. Series of selection detection methods, based on various theoretical models and exploiting different aspects of selection signatures, have been developed. Their applications at the species and population levels have respectively led to the identification of human specific selection events that distinguish human from nonhuman primates and local adaptation events that contribute to human diversity. Scrutiny of candidate genes has revealed paradigms of adaptations to specific nutritional components and genome-wide selection scans have verified the prevalence of diet-related selection events and provided many more candidates awaiting further investigation. Understanding the role of diet in human evolution is fundamental for the development of evidence-based, genome-informed nutritional practices in the era of personal genomics.

Genomewide Scan Reveals Amplification of mdr1 as a Common Denominator of Resistance to Mefloquine, Lumefantrine, and Artemisinin in Plasmodium chabaudi Malaria Parasites

Multidrug-resistant Plasmodium falciparum malaria parasites pose a threat to effective drug control, even to artemisinin-based combination therapies (ACTs). Here we used linkage group selection and Solexa whole-genome resequencing to investigate the genetic basis of resistance to component drugs of ACTs. Using the rodent malaria parasite P. chabaudi, we analyzed the uncloned progeny of a genetic backcross between the mefloquine-, lumefantrine-, and artemisinin-resistant mutant AS-15MF and a genetically distinct sensitive clone, AJ, following drug treatment. Genomewide scans of selection showed that parasites surviving each drug treatment bore a duplication of a segment of chromosome 12 (translocated to chromosome 04) present in AS-15MF. Whole-genome resequencing identified the size of the duplicated segment and its position on chromosome 4. The duplicated fragment extends for ∼393 kbp and contains over 100 genes, including mdr1, encoding the multidrug resistance P-glycoprotein homologue 1. We therefore show that resistance to chemically distinct components of ACTs is mediated by the same genetic mutation, highlighting a possible limitation of these therapies.

Detecting Directional Selection in the Presence of Recent Admixture in African-Americans

We investigate the performance of tests of neutrality in admixed populations using plausible demographic models for African-American history as well as resequencing data from African and African-American populations. The analysis of both simulated and human resequencing data suggests that recent admixture does not result in an excess of false-positive results for neutrality tests based on the frequency spectrum after accounting for the population growth in the parental African population. Furthermore, when simulating positive selection, Tajima's D, Fu and Li's D, and haplotype homozygosity have lower power to detect population-specific selection using individuals sampled from the admixed population than from the nonadmixed population. Fay and Wu's H test, however, has more power to detect selection using individuals from the admixed population than from the nonadmixed population, especially when the selective sweep ended long ago. Our results have implications for interpreting recent genome-wide scans for positive selection in human populations.

A Grand Challenge in Evolutionary and Population Genetics: New Paradigms for Exploring the Past and Charting the Future in the Post-Genomic era

A Grand Challenge in Evolutionary and Population Genetics: New Paradigms for Exploring the Past and Charting the Future in the Post-Genomic era

MADS-box out of the black box

The compelling elegance of using genome-wide scans to detect the signature of selection is difficult to resist, but is countered by the low demonstrated efficacy of pinpointing the actual genes and traits that are the targets of selection in nonmodel species. While the difficulty of going from a suggestive signature to a functional nucleotide polymorphism should not prevent researchers from using genome scans, it does lessen their long-term utility within and across study systems. In a new study published in this issue of Molecular Ecology (Mariac etal. 2011), researchers have gone a long way towards increasing the relevance of genome-wide scans for selection via two approaches: (i) they tailored the markers used in the scan to target a family of developmental genes that were good candidates for controlling a trait of interest and (ii) they used an independent mapping population to confirm the association of the gene with polymorphism in the trait of interest. All of this was completed in the nonmodel system of pearl millet (Pennisetum glaucum) and may provide a road map for other researchers hoping to pin down solid candidate genes for selected traits in natural or cultivated systems. Outside of these broad methodological innovations, the paper specifically focuses on a trait (flowering time) that varies across an environmental gradient (rainfall). This environmental gradient potentially serves as a model for environmental change over time, and allele frequencies at the gene can therefore be used to track how populations of pearl millet will adapt to future climate shifts at the genetic level.

Scanning for Genomic Regions Subject to Selective Sweeps Using SNP-MaP Strategy

Population genomic approaches, which take advantages of high-throughput genotyping, are powerful yet costly methods to scan for selective sweeps. DNA-pooling strategies have been widely used for association studies because it is a cost-effective alternative to large-scale individual genotyping. Here, we performed an SNP-MaP (single nucleotide polymorphism microarrays and pooling) analysis using samples from Eurasia to evaluate the efficiency of pooling strategy in genome-wide scans for selection. By conducting simulations of allelotype data, we first demonstrated that the boxplot with average heterozygosity (HET) is a promising method to detect strong selective sweeps with a moderate level of pooling error. Based on this, we used a sliding window analysis of HET to detect the large contiguous regions (LCRs) putatively under selective sweeps from Eurasia datasets. This survey identified 63 LCRs in a European population. These signals were further supported by the integrated haplotype score (iHS) test using HapMap II data. We also confirmed the European-specific signatures of positive selection from several previously identified genes (KEL, TRPV5, TRPV6, EPHB6). In summary, our results not only revealed the high credibility of SNP-MaP strategy in scanning for selective sweeps, but also provided an insight into the population differentiation.

A Novel PHD2 Mutation Associated with Tibetan Genetic Adaptation to High Altitude Hypoxia.

AbstractAbstract 2602The ability of the body to maintain oxygen hemostasis is essential for its survival. The body is required to adapt to the physiologic constraints imposed by a hypoxic environment in order to ensure optimal physiological functions. The study of indigenous high-altitude residents provides an opportunity to identify genes that may have play a role in hypoxia adaptation, and recent studies indicate that Tibetans have genetically adapted to the extremely hypoxic environment at high altitude. Interestingly, Tibetan adaption includes protection from polycythemia. We recently reported genome-wide scans of selection (Affymetrix Human SNP chips v6.0) for Tibetans (Simonson, Science July, 2010). At least 10 chromosomal regions were identified as likely targets for high-altitude adaptation, and two of these regions, EGLN1 (PHD2) and PPARA, were associated with protection against high-altitude polycythemia which is unique to this indigenous highland population. PHD2 is a negative regulator of the alpha subunits of HIF-1 and HIF-2 transcription factors regulating many essential functions including erythropoiesis, iron metabolism, development, in addition to other functions. A gain of function mutation of the alpha subunit of HIF-2 is associated with congenital polycythemia. We sequenced the exons of PPARA and PHD2 by direct sequencing and found two missense mutations in exon 1 of the PHD2 gene: the first is located within codon 4, and changes Aspartate to Glutamate (Asp4Glu) and the second is an unvalidated SNP (rs12097901) at codon 127 that changes Cysteine to Serine (Cy127Ser). The samples with high Hb show a heterozygous change, while the subject with normal Hb exhibited a homozygous change in both codons of the PHD2 gene.The association of these two missense PHD2 mutations to hemoglobin level is now being evaluated by studies of a larger population of Tibetan natives and the functional effects on levels and activities of HIF-1 and HIF-2 genes evaluated by transfection of these mutated PHD2 cDNA constructs in reporter cells.Table 1List of samples and controls screened for mutation in PHD2.SampleRaceSexHbPHD2 MutationAsp4GluCys127SerT03TibetanF194HeteroHeteroT43TibetanM194HeteroHeteroT48TibetanM101HeteroHomoT53TibetanM159HomoHomoControlCaucasianM?WTWTControlCaucasianF?WTWT Disclosure:No relevant conflicts of interest to declare.

The Impact of Natural Selection on an ABCC11 SNP Determining Earwax Type

A nonsynonymous single nucleotide polymorphism (SNP), rs17822931-G/A (538G>A; Gly180Arg), in the ABCC11 gene determines human earwax type (i.e., wet or dry) and is one of most differentiated nonsynonymous SNPs between East Asian and African populations. A recent genome-wide scan for positive selection revealed that a genomic region spanning ABCC11, LONP2, and SIAH1 genes has been subjected to a selective sweep in East Asians. Considering the potential functional significance as well as the population differentiation of SNPs located in that region, rs17822931 is the most plausible candidate polymorphism to have undergone geographically restricted positive selection. In this study, we estimated the selection intensity or selection coefficient of rs17822931-A in East Asians by analyzing two microsatellite loci flanking rs17822931 in the African (HapMap-YRI) and East Asian (HapMap-JPT and HapMap-CHB) populations. Assuming a recessive selection model, a coalescent-based simulation approach suggested that the selection coefficient of rs17822931-A had been approximately 0.01 in the East Asian population, and a simulation experiment using a pseudo-sampling variable revealed that the mutation of rs17822931-A occurred 2006 generations (95% credible interval, 1,023-3,901 generations) ago. In addition, we show that absolute latitude is significantly associated with the allele frequency of rs17822931-A in Asian, Native American, and European populations, implying that the selective advantage of rs17822931-A is related to an adaptation to a cold climate. Our results provide a striking example of how local adaptation has played a significant role in the diversification of human traits.

Plasmodium falciparum genome-wide scans for positive selection, recombination hot spots and resistance to antimalarial drugs.

Antimalarial drugs impose strong pressure on Plasmodium falciparum parasites and leave signatures of selection in the parasite genome 1,2. Search for signals of selection may lead to genes encoding drug or immune targets 3. The lack of high-throughput genotyping methods, inadequate knowledge of parasite population history, and time-consuming adaptations of parasites to in vitro culture have hampered genome-wide association studies (GWAS) of parasite traits. Here we report genotyping of DNA from 189 culture-adapted P. falciparum parasites using a custom-built array with thousands of single nucleotide polymorphisms (SNPs). Population structure, variation in recombination rate, and loci under recent positive selection were detected. Parasite half maximum inhibitory concentrations (IC50) to seven antimalarial drugs were obtained and used in GWAS to identify genes associated with drug responses. The SNP array and genome-wide parameters provide valuable tools and information for new advances in P. falciparum genetics.