Selection affecting enzyme polymorphisms in enclosed Drosophila populations maintained in a natural environment

Whole-Genome Sequencing of Native Sheep Provides Insights into Rapid Adaptations to Extreme Environments

Global climate change has a significant effect on extreme environments and a profound influence on species survival. However, little is known of the genome-wide pattern of livestock adaptations to extreme environments over a short time frame following domestication. Sheep (Ovis aries) have become well adapted to a diverse range of agroecological zones, including certain extreme environments (e.g., plateaus and deserts), during their post-domestication (approximately 8–9 kya) migration and differentiation. Here, we generated whole-genome sequences from 77 native sheep, with an average effective sequencing depth of ∼5× for 75 samples and ∼42× for 2 samples. Comparative genomic analyses among sheep in contrasting environments, that is, plateau (>4,000 m above sea level) versus lowland (<100 m), high-altitude region (>1500 m) versus low-altitude region (<1300 m), desert (<10 mm average annual precipitation) versus highly humid region (>600 mm), and arid zone (<400 mm) versus humid zone (>400 mm), detected a novel set of candidate genes as well as pathways and GO categories that are putatively associated with hypoxia responses at high altitudes and water reabsorption in arid environments. In addition, candidate genes and GO terms functionally related to energy metabolism and body size variations were identified. This study offers novel insights into rapid genomic adaptations to extreme environments in sheep and other animals, and provides a valuable resource for future research on livestock breeding in response to climate change.

Adaptations to climate-mediated selective pressures in humans

Humans inhabit a remarkably diverse range of environments, and adaptation through natural selection has likely played a central role in the capacity to survive and thrive in extreme climates. Unlike numerous studies that used only population genetic data to search for evidence of selection, here we scan the human genome for selection signals by identifying the SNPs with the strongest correlations between allele frequencies and climate across 61 worldwide populations. We find a striking enrichment of genic and nonsynonymous SNPs relative to non-genic SNPs among those that are strongly correlated with these climate variables. Among the most extreme signals, several overlap with those from GWAS, including SNPs associated with pigmentation and autoimmune diseases. Further, we find an enrichment of strong signals in gene sets related to UV radiation, infection and immunity, and cancer. Our results imply that adaptations to climate shaped the spatial distribution of variation in humans.

Classical studies that examined the global distributions of human physiological traits such as pigmentation, basal metabolic rate, and body shape and size suggested that natural selection related to climate has been important during recent human evolutionary history. We scanned the human genome using data for about 650,000 variants in 61 worldwide populations to look for correlations between allele frequencies and 9 climate variables and found evidence for adaptations to climate at the genome-wide level. In addition, we detected compelling signals for individual SNPs involved in pigmentation and immune response, as well as for pathways related to UV radiation, infection and immunity, and cancer. A particularly appealing aspect of this approach is that we identify a set of candidate advantageous SNPs associated with specific biological hypotheses, which will be useful for follow-up testing. We developed an online resource to browse the results of our data analyses, allowing researchers to quickly assess evidence for selection in a particular genomic region and to compare it across several studies.

Adaptations to climate in candidate genes for common metabolic disorders

Evolutionary pressures due to variation in climate play an important role in shaping phenotypic variation among and within species and have been shown to influence variation in phenotypes such as body shape and size among humans. Genes involved in energy metabolism are likely to be central to heat and cold tolerance. To test the hypothesis that climate shaped variation in metabolism genes in humans, we used a bioinformatics approach based on network theory to select 82 candidate genes for common metabolic disorders. We genotyped 873 tag SNPs in these genes in 54 worldwide populations (including the 52 in the Human Genome Diversity Project panel) and found correlations with climate variables using rank correlation analysis and a newly developed method termed Bayesian geographic analysis. In addition, we genotyped 210 carefully matched control SNPs to provide an empirical null distribution for spatial patterns of allele frequency due to population history alone. For nearly all climate variables, we found an excess of genic SNPs in the tail of the distributions of the test statistics compared to the control SNPs, implying that metabolic genes as a group show signals of spatially varying selection. Among our strongest signals were several SNPs (e.g., LEPR R109K, FABP2 A54T) that had previously been associated with phenotypes directly related to cold tolerance. Since variation in climate may be correlated with other aspects of environmental variation, it is possible that some of the signals that we detected reflect selective pressures other than climate. Nevertheless, our results are consistent with the idea that climate has been an important selective pressure acting on candidate genes for common metabolic disorders.

Extensive amino acid polymorphism at the pgm locus is consistent with adaptive protein evolution in Drosophila melanogaster

PGM plays a central role in the glycolytic pathway at the branch point leading to glycogen metabolism and is highly polymorphic in allozyme studies of many species. We have characterized the nucleotide diversity across the Pgm gene in Drosophila melanogaster and D. simulans to investigate the role that protein polymorphism plays at this crucial metabolic branch point shared with several other enzymes. Although D. melanogaster and D. simulans share common allozyme mobility alleles, we find these allozymes are the result of many different amino acid changes at the nucleotide level. In addition, specific allozyme classes within species contain several amino acid changes, which may explain the absence of latitudinal clines for PGM allozyme alleles, the lack of association of PGM allozymes with the cosmopolitan In(3L)P inversion, and the failure to detect differences between PGM allozymes in functional studies. We find a significant excess of amino acid polymorphisms within D. melanogaster when compared to the complete absence of fixed replacements with D. simulans. There is also strong linkage disequilibrium across the 2354 bp of the Pgm locus, which may be explained by a specific amino acid haplotype that is high in frequency yet contains an excess of singleton polymorphisms. Like G6pd, Pgm shows strong evidence for a branch point enzyme that exhibits adaptive protein evolution.

Causes and consequences of esterase 6 enzyme activity variation in pre-adult Drosophila melanogaster

We report heritable threefold differences in both larval and pupal esterase 6 activity among 17 isoallelic lines of D. melanogaster extracted from a natural population. The activity differences in the two stages are only weakly correlated with each other or with previously determined values for esterase 6 activity in adults of these lines. The pre-adult activity variation is also unrelated to polymorphisms among the lines for six esterase 6 allozymes and six restriction sites in a region encompassing the esterase 6 coding DNA and the first kbp of 5' flanking DNA. However, two insertions, of 8.0 and 6.8 kbp, located about 1.4 kbp 5' of the esterase 6 coding region are associated with low activity in larvae and, to a lesser extent, in pupae, albeit not in adults. Restriction mapping reveals similarity between the 8.0 kbp insert and the 7.4 kbp retrotransposon 17.6. The differences in larval activity among lines are positively correlated with fitness as assessed from assays of pre-adult viability and development time but no significant associations between pupal esterase 6 activity and these measures are detected. Some effects of esterase 6 allozyme differences are also found for viability and development time but these effects could be explained by linkage disequilibrium between the 8.0 kbp insert and the EST6-9 allozyme.

Evolutionary genetics of Drosophila esterases

Over 30 carboxylester hydrolases have been identified in D. melanogaster. Most are classified as acetyl, carboxyl or cholinesterases. Sequence similarities among most of the carboxyl and all the cholinesterases so far characterised from D. melanogaster and other eukaryotes justify recognition of a carboxyl/cholinesterase multigene family. This family shows minimal sequence similarities with other esterases but crystallographic data for a few non-drosophilid enzymes show that the family shares a distinctive overall structure with some other carboxyl and aryl esterases, so they are all put in one superfamily of/beta hydrolases. Fifteen esterase genes have been mapped in D. melanogaster and twelve are clustered at two chromosomal sites. The constitution of each cluster varies across Drosophila species but two carboxyl esterases in one cluster are sufficiently conserved that their homologues can be identified among enzymes conferring insecticide resistance in other Diptera. Sequence differences between two other esterases, the EST6 carboxyl esterase and acetylcholinesterase, have been interpreted against the consensus super-secondary structure for the carboxyl/cholinesterase multigene family; their sequence differences are widely dispersed across the structure and include substantial divergence in substrate binding sites and the active site gorge. This also applies when EST6 is compared across species where differences in its expression indicate a difference in function. However, comparisons within and among species where EST6 expression is conserved show that many aspects of the predicted super-secondary structure are tightly conserved. Two notable exceptions are a pair of polymorphisms in the substrate binding site of the enzyme in D. melanogaster. These polymorphisms are associated with differences in substrate interactions in vitro and demographic data indicate that the alternative forms are not selectively equivalent in vivo.

The alcohol dehydrogenase polymorphism in natural populations of Drosophila melanogaster: restriction map variation in the region of the Adh locus in populations from two hemispheres

Restriction endonuclease variation in the 12 kb region surrounding the Adh locus was measured in seven Australian and six Chinese populations of Drosophila melanogaster. There is a higher level of nucleotide-substitution variation in the Australian populations than in the Chinese, which is possibly a reflection of their origins. None of the restriction site polymorphisms, nor any of the insertions, showed a significant association with latitude. A 0.2 kb deletion varied with latitude in the Chinese populations. In accordance with previous studies, a majority of the insertions were located in a region 1.5-3.5 kb 3' from the Adh coding region, and a majority of the deletions were at a site 3 kb 5' to the Adh coding region. Two of the insertions shared homologies with known mobile elements. Overall, the data suggest that restriction endonuclease variation in the Adh region is not related to the cline in Adhs frequencies.

Genetic differentiation in populations of Drosophila melanogaster from the Peoples' Republic of China: comparison with patterns on other continents

Geographic variation at eight allozyme loci in eight natural populations of Drosophila melanogaster from the People's Republic of China (P.R.C.) was collated with data from Japanese populations and compared with previous reports for other continents. G6pdF, Est-6(1.00) and AdhS were significantly correlated with latitude. Whilst the variation in Est6(1.00) was opposite to that previously reported, AdhS showed latitudinal clines consistent with data from the northern and southern hemispheres. A thermostable variant, AdhFChD, was found at high frequency in the southern P.R.C. populations and it is suggested that the mutation occurred in this region and was then dispersed to other continents.