Genome-Wide Analysis of Cold Adaptation in Indigenous Siberian Populations

Following the dispersal out of Africa, where hominins evolved in warm environments for millions of years, our species has colonised different climate zones of the world, including high latitudes and cold environments. The extent to which human habitation in (sub-)Arctic regions has been enabled by cultural buffering, short-term acclimatization and genetic adaptations is not clearly understood. Present day indigenous populations of Siberia show a number of phenotypic features, such as increased basal metabolic rate, low serum lipid levels and increased blood pressure that have been attributed to adaptation to the extreme cold climate. In this study we introduce a dataset of 200 individuals from ten indigenous Siberian populations that were genotyped for 730,525 SNPs across the genome to identify genes and non-coding regions that have undergone unusually rapid allele frequency and long-range haplotype homozygosity change in the recent past. At least three distinct population clusters could be identified among the Siberians, each of which showed a number of unique signals of selection. A region on chromosome 11 (chr11:66–69 Mb) contained the largest amount of clustering of significant signals and also the strongest signals in all the different selection tests performed. We present a list of candidate cold adaption genes that showed significant signals of positive selection with our strongest signals associated with genes involved in energy regulation and metabolism (CPT1A, LRP5, THADA) and vascular smooth muscle contraction (PRKG1). By employing a new method that paints phased chromosome chunks by their ancestry we distinguish local Siberian-specific long-range haplotype signals from those introduced by admixture.     

Experimental Evolution of an RNA Virus in Wild Birds: Evidence for Host-Dependent Impacts on Population Structure and Competitive Fitness

Within hosts, RNA viruses form populations that are genetically and phenotypically complex. Heterogeneity in RNA virus genomes arises due to error-prone replication and is reduced by stochastic and selective mechanisms that are incompletely understood. Defining how natural selection shapes RNA virus populations is critical because it can inform treatment paradigms and enhance control efforts. We allowed West Nile virus (WNV) to replicate in wild-caught American crows, house sparrows and American robins to assess how natural selection shapes RNA virus populations in ecologically relevant hosts that differ in susceptibility to virus-induced mortality. After five sequential passages in each bird species, we examined the phenotype and population diversity of WNV through fitness competition assays and next generation sequencing. We demonstrate that fitness gains occur in a species-specific manner, with the greatest replicative fitness gains in robin-passaged WNV and the least in WNV passaged in crows. Sequencing data revealed that intrahost WNV populations were strongly influenced by purifying selection and the overall complexity of the viral populations was similar among passaged hosts. However, the selective pressures that control WNV populations seem to be bird species-dependent. Specifically, crow-passaged WNV populations contained the most unique mutations (~1.7× more than sparrows, ~3.4× more than robins) and defective genomes (~1.4× greater than sparrows, ~2.7× greater than robins), but the lowest average mutation frequency (about equal to sparrows, ~2.6× lower than robins). Therefore, our data suggest that WNV replication in the most disease-susceptible bird species is positively associated with virus mutational tolerance, likely via complementation, and negatively associated with the strength of selection. These differences in genetic composition most likely have distinct phenotypic consequences for the virus populations. Taken together, these results reveal important insights into how different hosts may contribute to the emergence of RNA viruses.     

Hypoxia Adaptations in the Grey Wolf (Canis lupus chanco) from Qinghai-Tibet Plateau

The Tibetan grey wolf (Canis lupus chanco) occupies habitats on the Qinghai-Tibet Plateau, a high altitude (>3000 m) environment where low oxygen tension exerts unique selection pressure on individuals to adapt to hypoxic conditions. To identify genes involved in hypoxia adaptation, we generated complete genome sequences of nine Chinese wolves from high and low altitude populations at an average coverage of 25× coverage. We found that, beginning about 55,000 years ago, the highland Tibetan grey wolf suffered a more substantial population decline than lowland wolves. Positively selected hypoxia-related genes in highland wolves are enriched in the HIF signaling pathway (P = 1.57E-6), ATP binding (P = 5.62E-5), and response to an oxygen-containing compound (P≤5.30E-4). Of these positively selected hypoxia-related genes, three genes (EPAS1, ANGPT1, and RYR2) had at least one specific fixed non-synonymous SNP in highland wolves based on the nine genome data. Our re-sequencing studies on a large panel of individuals showed a frequency difference greater than 58% between highland and lowland wolves for these specific fixed non-synonymous SNPs and a high degree of LD surrounding the three genes, which imply strong selection. Past studies have shown that EPAS1 and ANGPT1 are important in the response to hypoxic stress, and RYR2 is involved in heart function. These three genes also exhibited significant signals of natural selection in high altitude human populations, which suggest similar evolutionary constraints on natural selection in wolves and humans of the Qinghai-Tibet Plateau.     

Low Frequency Variants,Collapsed Based on Biological Knowledge,Uncover Complexity of Population Stratification in 1000 Genomes Project Data

Analyses investigating low frequency variants have the potential for explaining additional genetic heritability of many complex human traits. However, the natural frequencies of rare variation between human populations strongly confound genetic analyses. We have applied a novel collapsing method to identify biological features with low frequency variant burden differences in thirteen populations sequenced by the 1000 Genomes Project. Our flexible collapsing tool utilizes expert biological knowledge from multiple publicly available database sources to direct feature selection. Variants were collapsed according to genetically driven features, such as evolutionary conserved regions, regulatory regions genes, and pathways. We have conducted an extensive comparison of low frequency variant burden differences (MAF<0.03) between populations from 1000 Genomes Project Phase I data. We found that on average 26.87% of gene bins, 35.47% of intergenic bins, 42.85% of pathway bins, 14.86% of ORegAnno regulatory bins, and 5.97% of evolutionary conserved regions show statistically significant differences in low frequency variant burden across populations from the 1000 Genomes Project. The proportion of bins with significant differences in low frequency burden depends on the ancestral similarity of the two populations compared and types of features tested. Even closely related populations had notable differences in low frequency burden, but fewer differences than populations from different continents. Furthermore, conserved or functionally relevant regions had fewer significant differences in low frequency burden than regions under less evolutionary constraint. This degree of low frequency variant differentiation across diverse populations and feature elements highlights the critical importance of considering population stratification in the new era of DNA sequencing and low frequency variant genomic analyses.     

A New Catalog of Structural Variants in 1,301 A. thaliana Lines from Africa,Eurasia, and North America Reveals a Signature of Balancing Selection at Defense Response Genes

Genomic variation in the model plant Arabidopsis thaliana has been extensively used to understand evolutionary processes in natural populations, mainly focusing on single-nucleotide polymorphisms. Conversely, structural variation has been largely ignored in spite of its potential to dramatically affect phenotype. Here, we identify 155,440 indels and structural variants ranging in size from 1 bp to 10 kb, including presence/absence variants (PAVs), inversions, and tandem duplications in 1,301 A. thaliana natural accessions from Morocco, Madeira, Europe, Asia, and North America. We show evidence for strong purifying selection on PAVs in genes, in particular for housekeeping genes and homeobox genes, and we find that PAVs are concentrated in defense-related genes (R-genes, secondary metabolites) and F-box genes. This implies the presence of a “core” genome underlying basic cellular processes and a “flexible” genome that includes genes that may be important in spatially or temporally varying selection. Further, we find an excess of intermediate frequency PAVs in defense response genes in nearly all populations studied, consistent with a history of balancing selection on this class of genes. Finally, we find that PAVs in genes involved in the cold requirement for flowering (vernalization) and drought response are strongly associated with temperature at the sites of origin.     

Rapid evolution of cold tolerance in stickleback

Climate change is predicted to lead to increased average temperatures and greater intensity and frequency of high and low temperature extremes, but the evolutionary consequences for biological communities are not well understood. Studies of adaptive evolution of temperature tolerance have typically involved correlative analyses of natural populations or artificial selection experiments in the laboratory. Field experiments are required to provide estimates of the timing and strength of natural selection, enhance understanding of the genetics of adaptation and yield insights into the mechanisms driving evolutionary change. Here, we report the experimental evolution of cold tolerance in natural populations of threespine stickleback fish (Gasterosteus aculeatus). We show that freshwater sticklebacks are able to tolerate lower minimum temperatures than marine sticklebacks and that this difference is heritable. We transplanted marine sticklebacks to freshwater ponds and measured the rate of evolution after three generations in this environment. Cold tolerance evolved at a rate of 0.63 haldanes to a value 2.5°C lower than that of the ancestral population, matching values found in wild freshwater populations. Our results suggest that cold tolerance is under strong selection and that marine sticklebacks carry sufficient genetic variation to adapt to changes in temperature over remarkably short time scales.     

Evidence for natural selection on leukocyte immunoglobulin-like receptors for HLA class I in Northeast Asians

Human leukocyte antigen (HLA) plays a critical role in innate and adaptive immunity and is a well-known example of genes under natural selection. However, the genetic aspect of receptors recognizing HLA molecules has not yet been fully elucidated. Leukocyte immunoglobulin (Ig)-like receptors (LILRs) are a family of HLA class I-recognizing receptors comprising activating and inhibitory forms. We previously reported that the allele frequency of the 6.7 kb LILRA3 deletion is extremely high (71%) in the Japanese population, and we identified premature termination codon (PTC)-containing alleles. In this study, we observed a wide distribution of the high deletion frequency in Northeast Asians (84% in Korean Chinese, 79% in Man Chinese, 56% in Mongolian, and 76% in Buryat populations). Genotyping of the four HapMap populations revealed that LILRA3 alleles were in strong linkage disequilibrium with LILRB2 alleles in Northeast Asians. In addition, PTC-containing LILRA3 alleles were detected in Northeast Asians but not in non-Northeast Asians. Furthermore, flow-cytometric analysis revealed that the LILRB2 allele frequent in Northeast Asians was significantly associated with low levels of expression. F(ST) and extended-haplotype-homozygosity analysis for the HapMap populations provided evidence of positive selection acting on the LILRA3 and LILRB2 loci. Taken together, our results suggest that both the nonfunctional LILRA3 alleles and the low-expressing LILRB2 alleles identified in our study have increased in Northeast Asians because of natural selection. Our findings, therefore, lead us to speculate that not only HLA class I ligands but also their receptors might be sensitive to the local environment.     

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.     

Response of Two Heat Shock Genes to Selection for Knockdown Heat Resistance in Drosophila Melanogaster

To identify genes involved in stress resistance and heat hardening, replicate lines of Drosophila melanogaster were selected for increased resistance to knockdown by a 39° heat stress. Two selective regimes were used, one with and one without prior hardening. Mean knockdown times were increased from ~5 min to >20 min after 18 generations. Initial realized heritabilities were as high as 10% for lines selected without hardening, and crosses between lines indicated simple additive gene effects for the selected phenotypes. To survey allelic variation and correlated selection responses in two candidate stress genes, hsr-omega and hsp68, we applied denaturing gradient gel electrophoresis to amplified DNA sequences from small regions of these genes. After eight generations of selection, allele frequencies at both loci showed correlated responses for selection following hardening, but not without hardening. The hardening process itself was associated with a hsp68 frequency change in the opposite direction to that associated with selection that followed hardening. These stress loci are closely linked on chromosome III, and the hardening selection established a disequilibrium, suggesting an epistatic effect on resistance. The data indicate that molecular variation in both hsr-omega and hsp68 contribute to natural heritable variation for hardened heat resistance.     

Genomic Signatures of Selective Pressures and Introgression from Archaic Hominins at Human Innate Immunity Genes

Human genes governing innate immunity provide a valuable tool for the study of the selective pressure imposed by microorganisms on host genomes. A comprehensive, genome-wide study of how selective constraints and adaptations have driven the evolution of innate immunity genes is missing. Using full-genome sequence variation from the 1000 Genomes Project, we first show that innate immunity genes have globally evolved under stronger purifying selection than the remainder of protein-coding genes. We identify a gene set under the strongest selective constraints, mutations in which are likely to predispose individuals to life-threatening disease, as illustrated by STAT1 and TRAF3. We then evaluate the occurrence of local adaptation and detect 57 high-scoring signals of positive selection at innate immunity genes, variation in which has been associated with susceptibility to common infectious or autoimmune diseases. Furthermore, we show that most adaptations targeting coding variation have occurred in the last 6,000–13,000 years, the period at which populations shifted from hunting and gathering to farming. Finally, we show that innate immunity genes present higher Neandertal introgression than the remainder of the coding genome. Notably, among the genes presenting the highest Neandertal ancestry, we find the TLR6-TLR1-TLR10 cluster, which also contains functional adaptive variation in Europeans. This study identifies highly constrained genes that fulfill essential, non-redundant functions in host survival and reveals others that are more permissive to change—containing variation acquired from archaic hominins or adaptive variants in specific populations—improving our understanding of the relative biological importance of innate immunity pathways in natural conditions.     

Evaluating the Effects of Chlortetracycline on the Proliferation of Antibiotic-Resistant Bacteria in a Simulated River Water Ecosystem

Antibiotics and antibiotic metabolites have been found in the environment, but the biological activities of these compounds are uncertain, especially given the low levels that are typically detected in the environment. The objective of this study was to estimate the selection potential of chlortetracycline (CTC) on the antibiotic resistance of aerobic bacterial populations in a simulated river water ecosystem. Six replicates of a 10-day experiment using river water in continuous flow chemostat systems were conducted. Each replicate used three chemostats, one serving as a control to which no antibiotic was added and the other two receiving low and high doses of CTC (8 μg/liter and 800 μg/liter, respectively). The addition of CTC to the chemostats did not impact the overall level of cultivable aerobic bacteria (P = 0.51). The high-CTC chemostat had significantly higher tetracycline-resistant bacterial colony counts than both the low-CTC and the control chemostats (P < 0.035). The differences in resistance between the low-CTC and control chemostats were highly nonsignificant (P = 0.779). In general a greater diversity of tet resistance genes was detected in the high-CTC chemostat and with a greater frequency than in the low-CTC and control chemostats. Low levels of CTC in this in vitro experiment did not select for increased levels of tetracycline resistance among cultivable aerobic bacteria. This finding should not be equated with the absence of environmental risk, however. Low concentrations of antibiotics in the environment may select for resistant bacterial populations once they are concentrated in sediments or other locations.     

Repeatability of adaptation in experimental populations of different sizes

The degree to which evolutionary trajectories and outcomes are repeatable across independent populations depends on the relative contribution of selection, chance and history. Population size has been shown theoretically and empirically to affect the amount of variation that arises among independent populations adapting to the same environment. Here, we measure the contribution of selection, chance and history in different-sized experimental populations of the unicellular alga Chlamydomonas reinhardtii adapting to a high salt environment to determine which component of evolution is affected by population size. We find that adaptation to salt is repeatable at the fitness level in medium (N_e = 5 × 10⁴) and large (N_e = 4 × 10⁵) populations because of the large contribution of selection. Adaptation is not repeatable in small (N_e = 5 × 10³) populations because of large constraints from history. The threshold between stochastic and deterministic evolution in this case is therefore between effective population sizes of 10³ and 10⁴. Our results indicate that diversity across populations is more likely to be maintained if they are small. Experimental outcomes in large populations are likely to be robust and can inform our predictions about outcomes in similar situations.     

Allele frequency divergence reveals ubiquitous influence of positive selection in Drosophila

Resolving the role of natural selection is a basic objective of evolutionary biology. It is generally difficult to detect the influence of selection because ubiquitous non-selective stochastic change in allele frequencies (genetic drift) degrades evidence of selection. As a result, selection scans typically only identify genomic regions that have undergone episodes of intense selection. Yet it seems likely such episodes are the exception; the norm is more likely to involve subtle, concurrent selective changes at a large number of loci. We develop a new theoretical approach that uncovers a previously undocumented genome-wide signature of selection in the collective divergence of allele frequencies over time. Applying our approach to temporally resolved allele frequency measurements from laboratory and wild Drosophila populations, we quantify the selective contribution to allele frequency divergence and find that selection has substantial effects on much of the genome. We further quantify the magnitude of the total selection coefficient (a measure of the combined effects of direct and linked selection) at a typical polymorphic locus, and find this to be large (of order 1%) even though most mutations are not directly under selection. We find that selective allele frequency divergence is substantially elevated at intermediate allele frequencies, which we argue is most parsimoniously explained by positive—not negative—selection. Thus, in these populations most mutations are far from evolving neutrally in the short term (tens of generations), including mutations with neutral fitness effects, and the result cannot be explained simply as an ongoing purging of deleterious mutations.     

Polymorphic Variants of SCN1A and EPHX1 Influence Plasma Carbamazepine Concentration,Metabolism and Pharmacoresistance in a Population of Kosovar Albanian Epileptic Patients

Aim

The present study aimed to evaluate the effects of gene variants in key genes influencing pharmacokinetic and pharmacodynamic of carbamazepine (CBZ) on the response in patients with epilepsy.

Materials & Methods

Five SNPs in two candidate genes influencing CBZ transport and metabolism, namely ABCB1 or EPHX1, and CBZ response SCN1A (sodium channel) were genotyped in 145 epileptic patients treated with CBZ as monotherapy and 100 age and sex matched healthy controls. Plasma concentrations of CBZ, carbamazepine-10,11-epoxide (CBZE) and carbamazepine-10,11-trans dihydrodiol (CBZD) were determined by HPLC-UV-DAD and adjusted for CBZ dosage/kg of body weight.

Results

The presence of the SCN1A IVS5-91G>A variant allele is associated with increased epilepsy susceptibility. Furthermore, carriers of the SCN1A IVS5-91G>A variant or of EPHX1 c.337T>C variant presented significantly lower levels of plasma CBZ compared to carriers of the common alleles (0.71±0.28 vs 1.11±0.69 μg/mL per mg/Kg for SCN1A IVS5-91 AA vs GG and 0.76±0.16 vs 0.94±0.49 μg/mL per mg/Kg for EPHX1 c.337 CC vs TT; P<0.05 for both). Carriers of the EPHX1 c.416A>G showed a reduced microsomal epoxide hydrolase activity as reflected by a significantly decreased ratio of CBZD to CBZ (0.13±0.08 to 0.26±0.17, p<0.05) also of CBZD to CBZE (1.74±1.06 to 3.08±2.90; P<0.05) and CDR_CBZD (0.13±0.08 vs 0.24±0.19 μg/mL per mg/Kg; P<0.05). ABCB1 3455C>T SNP and SCN1A 3148A>G variants were not associated with significant changes in CBZ pharmacokinetic. Patients resistant to CBZ treatment showed increased dosage of CBZ (657±285 vs 489±231 mg/day; P<0.001) but also increased plasma levels of CBZ (9.84±4.37 vs 7.41±3.43 μg/mL; P<0.001) compared to patients responsive to CBZ treatment. CBZ resistance was not related to any of the SNPs investigated.

Conclusions

The SCN1A IVS5-91G>A SNP is associated with susceptibility to epilepsy. SNPs in EPHX1 gene are influencing CBZ metabolism and disposition. CBZ plasma levels are not an indicator of resistance to the therapy.     

The hemochromatosis 845 G-->A and 187 C-->G mutations: prevalence in non-Caucasian populations.

Hemochromatosis, the inherited disorder of iron metabolism, leads, if untreated, to progressive iron overload and premature death. The hemochromatosis gene, HFE, recently has been identified, and characterization of this gene has shown that it contains two mutations that result in amino acid substitutions-cDNA nucleotides 845 G-->A (C282Y) and 187 C-->G (H63D). Although hemochromatosis is common in Caucasians, affecting >=1/300 individuals of northern European origin, it has not been recognized in other populations. The present study used PCR and restriction-enzyme digestion to analyze the frequency of the 845 G-->A and 187 C-->G mutations in HLA-typed samples from non-Caucasian populations, comprising Australian Aboriginal, Chinese, and Pacific Islanders. Results showed that the 845 G-->A mutation was present in these populations (allele frequency 0.32%), and, furthermore, it was always seen in conjunction with HLA haplotypes common in Caucasians, suggesting that 845 G-->A may have been introduced into these populations by Caucasian admixture. 187 C-->G was present at an allele frequency of 2.68% in the two populations analyzed (Australian Aboriginal and Chinese). In the Australian Aboriginal samples, 187 C-->G was found to be associated with HLA haplotypes common in Caucasians, suggesting that it was introduced by recent admixture. In the Chinese samples analyzed, 187 C-->G was present in association with a wide variety of HLA haplotypes, showing this mutation to be widespread and likely to predate the more genetically restricted 845 G-->A mutation.     

The Role of Selection in Shaping Diversity of Natural M. tuberculosis Populations

Mycobacterium tuberculosis (M.tb), the cause of tuberculosis (TB), is estimated to infect a new host every second. While analyses of genetic data from natural populations of M.tb have emphasized the role of genetic drift in shaping patterns of diversity, the influence of natural selection on this successful pathogen is less well understood. We investigated the effects of natural selection on patterns of diversity in 63 globally extant genomes of M.tb and related pathogenic mycobacteria. We found evidence of strong purifying selection, with an estimated genome-wide selection coefficient equal to −9.5×10⁻⁴ (95% CI −1.1×10⁻³ to −6.8×10⁻⁴); this is several orders of magnitude higher than recent estimates for eukaryotic and prokaryotic organisms. We also identified different patterns of variation across categories of gene function. Genes involved in transport and metabolism of inorganic ions exhibited very low levels of non-synonymous polymorphism, equivalent to categories under strong purifying selection (essential and translation-associated genes). The highest levels of non-synonymous variation were seen in a group of transporter genes, likely due to either diversifying selection or local selective sweeps. In addition to selection, we identified other important influences on M.tb genetic diversity, such as a 25-fold expansion of global M.tb populations coincident with explosive growth in human populations (estimated timing 1684 C.E., 95% CI 1620–1713 C.E.). These results emphasize the parallel demographic histories of this obligate pathogen and its human host, and suggest that the dominant effect of selection on M.tb is removal of novel variants, with exceptions in an interesting group of genes involved in transportation and defense. We speculate that the hostile environment within a host imposes strict demands on M.tb physiology, and thus a substantial fitness cost for most new mutations. In this respect, obligate bacterial pathogens may differ from other host-associated microbes such as symbionts.     

Selection of very small differences in bacterial evolution.

As the Science of Biology is constantly changing due to new discoveries and advanced techniques it is essential that a systematic study of the environmental causes of natural selection on microorganisms be conducted. Very small phenotypic differences among individuals within bacterial populations arise as a result of spontaneous genetic variation, but the evolutionary importance of these small changes is frequently considered to be non-significant. Recent in vitro experiments indicate that efficient selection of these very small differences may take place in environmental compartments where a particular intensity of the selective agent is exerted. Model studies based on competition between bacterial populations only differing in one or two amino acid changes of a detoxifying antibiotic enzyme (e.g. beta-lactamase) have shown that at a narrow range of antibiotic concentrations the variant population is strongly selected over the original type, despite the extremely low phenotypic differences in antibiotic susceptibility. These selective concentrations are expected to occur in precise environmental compartments (selective compartments). Due to the high frequency of structured habitats in natural environments, the intensity of selective agents is commonly exerted along certain gradients. Each one of the points forming these gradients (or intersection among gradients) may have a particular selective ability for a specific genetic variant. Considering the environment as a composition of an extremely high number of specific selective compartments may help to understand the existence of high levels of genetic variability in natural bacterial populations. This may be one of the clues towards the unraveling of bacterial evolution.