FitSNPs: highly differentially expressed genes are more likely to have variants associated with disease

Background

Candidate single nucleotide polymorphisms (SNPs) from genome-wide association studies (GWASs) were often selected for validation based on their functional annotation, which was inadequate and biased. We propose to use the more than 200,000 microarray studies in the Gene Expression Omnibus to systematically prioritize candidate SNPs from GWASs.

Results

We analyzed all human microarray studies from the Gene Expression Omnibus, and calculated the observed frequency of differential expression, which we called differential expression ratio, for every human gene. Analysis conducted in a comprehensive list of curated disease genes revealed a positive association between differential expression ratio values and the likelihood of harboring disease-associated variants. By considering highly differentially expressed genes, we were able to rediscover disease genes with 79% specificity and 37% sensitivity. We successfully distinguished true disease genes from false positives in multiple GWASs for multiple diseases. We then derived a list of functionally interpolating SNPs (fitSNPs) to analyze the top seven loci of Wellcome Trust Case Control Consortium type 1 diabetes mellitus GWASs, rediscovered all type 1 diabetes mellitus genes, and predicted a novel gene (KIAA1109) for an unexplained locus 4q27. We suggest that fitSNPs would work equally well for both Mendelian and complex diseases (being more effective for cancer) and proposed candidate genes to sequence for their association with 597 syndromes with unknown molecular basis.

Conclusions

Our study demonstrates that highly differentially expressed genes are more likely to harbor disease-associated DNA variants. FitSNPs can serve as an effective tool to systematically prioritize candidate SNPs from GWASs.     

Balancing Selection on a Regulatory Region Exhibiting Ancient Variation That Predates Human–Neandertal Divergence

Ancient population structure shaping contemporary genetic variation has been recently appreciated and has important implications regarding our understanding of the structure of modern human genomes. We identified a ∼36-kb DNA segment in the human genome that displays an ancient substructure. The variation at this locus exists primarily as two highly divergent haplogroups. One of these haplogroups (the NE1 haplogroup) aligns with the Neandertal haplotype and contains a 4.6-kb deletion polymorphism in perfect linkage disequilibrium with 12 single nucleotide polymorphisms (SNPs) across diverse populations. The other haplogroup, which does not contain the 4.6-kb deletion, aligns with the chimpanzee haplotype and is likely ancestral. Africans have higher overall pairwise differences with the Neandertal haplotype than Eurasians do for this NE1 locus (p<10⁻¹⁵). Moreover, the nucleotide diversity at this locus is higher in Eurasians than in Africans. These results mimic signatures of recent Neandertal admixture contributing to this locus. However, an in-depth assessment of the variation in this region across multiple populations reveals that African NE1 haplotypes, albeit rare, harbor more sequence variation than NE1 haplotypes found in Europeans, indicating an ancient African origin of this haplogroup and refuting recent Neandertal admixture. Population genetic analyses of the SNPs within each of these haplogroups, along with genome-wide comparisons revealed significant F_ST (p = 0.00003) and positive Tajima''s D (p = 0.00285) statistics, pointing to non-neutral evolution of this locus. The NE1 locus harbors no protein-coding genes, but contains transcribed sequences as well as sequences with putative regulatory function based on bioinformatic predictions and in vitro experiments. We postulate that the variation observed at this locus predates Human–Neandertal divergence and is evolving under balancing selection, especially among European populations.     

Detection of selection signatures within candidate regions underlying trypanotolerance in outbred cattle populations

Breeding indigenous African taurine cattle tolerant to trypanosomosis is a straightforward approach to control costs generated by this disease. A recent study identified quantitative trait loci (QTL) underlying trypanotolerance traits in experimental crosses between tolerant N'Dama and susceptible Boran zebu cattle. As trypanotolerance is thought to result from local adaptation of indigenous cattle breeds, we propose an alternative and complementary approach to study the genetic architecture of this trait, based on the identification of selection signatures within QTL or candidate genes. A panel of 92 microsatellite markers was genotyped on 509 cattle belonging to four West African trypanotolerant taurine breeds and 10 trypanosusceptible European or African cattle breeds. Some of these markers were located within previously identified QTL regions or candidate genes, while others were chosen in regions assumed to be neutral. A detailed analysis of the genetic structure of these different breeds was carried out to confirm a priori grouping of populations based on previous data. Tests based on the comparison of the observed heterozygosities and variances in microsatellite allelic size among trypanotolerant and trypanosusceptible breeds led to the identification of two significantly less variable microsatellite markers. BM4440, one of these two outlier loci, is located within the confidence interval of a previously described QTL underlying a trypanotolerance-related trait.
Detection of selection signatures appears to be a straightforward approach for unravelling the molecular determinism of trypanosomosis pathogenesis. We expect that a whole genome approach will help confirm these results and achieve a higher resolving power.     

Cold adaptation in the Asian tiger mosquito's native range precedes its invasion success in temperate regions

Adaptation to environmental conditions within the native range of exotic species can condition the invasion success of these species outside their range. The striking success of the Asian tiger mosquito, Aedes albopictus, to invade temperate regions has been attributed to the winter survival of diapause eggs in cold environments. In this study, we evaluate genetic polymorphisms (SNPs) and wing morphometric variation among three biogeographical regions of the native range of A. albopictus. Reconstructed demographic histories of populations show an initial expansion in Southeast Asia and suggest that marine regression during late Pleistocene and climate warming after the last glacial period favored expansion of populations in southern and northern regions, respectively. Searching for genomic signatures of selection, we identified significantly differentiated SNPs among which several are located in or within 20 kb distance from candidate genes for cold adaptation. These genes involve cellular and metabolic processes and several of them have been shown to be differentially expressed under diapausing conditions. The three biogeographical regions also differ for wing size and shape, and wing size increases with latitude supporting Bergmann's rule. Adaptive genetic and morphometric variation observed along the climatic gradient of A. albopictus native range suggests that colonization of northern latitudes promoted adaptation to cold environments prior to its worldwide invasion.     

Genomic and phenotypic differentiation of Arabidopsis thaliana along altitudinal gradients in the North Italian Alps

Altitudinal gradients in mountain regions are short‐range clines of different environmental parameters such as temperature or radiation. We investigated genomic and phenotypic signatures of adaptation to such gradients in five Arabidopsis thaliana populations from the North Italian Alps that originated from 580 to 2350 m altitude by resequencing pools of 19–29 individuals from each population. The sample includes two pairs of low‐ and high‐altitude populations from two different valleys. High‐altitude populations showed a lower nucleotide diversity and negative Tajima's D values and were more closely related to each other than to low‐altitude populations from the same valley. Despite their close geographic proximity, demographic analysis revealed that low‐ and high‐altitude populations split between 260 000 and 15 000 years before present. Single nucleotide polymorphisms whose allele frequencies were highly differentiated between low‐ and high‐altitude populations identified genomic regions of up to 50 kb length where patterns of genetic diversity are consistent with signatures of local selective sweeps. These regions harbour multiple genes involved in stress response. Variation among populations in two putative adaptive phenotypic traits, frost tolerance and response to light/UV stress was not correlated with altitude. Taken together, the spatial distribution of genetic diversity reflects a potentially adaptive differentiation between low‐ and high‐altitude populations, whereas the phenotypic differentiation in the two traits investigated does not. It may resemble an interaction between adaptation to the local microhabitat and demographic history influenced by historical glaciation cycles, recent seed dispersal and genetic drift in local populations.     

Resequencing candidate genes implicates rare variants in asthma susceptibility

Common variation in over 100 genes has been implicated in the risk of developing asthma, but the contribution of rare variants to asthma susceptibility remains largely unexplored. We selected nine genes that showed the strongest signatures of weak purifying selection from among 53 candidate asthma-associated genes, and we sequenced the coding exons and flanking noncoding regions in 450 asthmatic cases and 515 nonasthmatic controls. We observed an overall excess of p values <0.05 (p = 0.02), and rare variants in four genes (AGT, DPP10, IKBKAP, and IL12RB1) contributed to asthma susceptibility among African Americans. Rare variants in IL12RB1 were also associated with asthma susceptibility among European Americans, despite the fact that the majority of rare variants in IL12RB1 were specific to either one of the populations. The combined evidence of association with rare noncoding variants in IL12RB1 remained significant (p = 3.7 × 10(-4)) after correcting for multiple testing. Overall, the contribution of rare variants to asthma susceptibility was predominantly due to noncoding variants in sequences flanking the exons, although nonsynonymous rare variants in DPP10 and in IL12RB1 were associated with asthma in African Americans and European Americans, respectively. This study provides evidence that rare variants contribute to asthma susceptibility. Additional studies are required for testing whether prioritizing genes for resequencing on the basis of signatures of purifying selection is an efficient means of identifying novel rare variants that contribute to complex disease.     

Larger genetic differences within africans than between Africans and Eurasians

The worldwide pattern of single nucleotide polymorphism (SNP) variation is of great interest to human geneticists, population geneticists, and evolutionists, but remains incompletely understood. We studied the pattern in noncoding regions, because they are less affected by natural selection than are coding regions. Thus, it can reflect better the history of human evolution and can serve as a baseline for understanding the maintenance of SNPs in human populations. We sequenced 50 noncoding DNA segments each approximately 500 bp long in 10 Africans, 10 Europeans, and 10 Asians. An analysis of the data suggests that the sampling scheme is adequate for our purpose. The average nucleotide diversity (pi) for the 50 segments is only 0.061% +/- 0.010% among Asians and 0.064% +/- 0.011% among Europeans but almost twice as high (0.115% +/- 0.016%) among Africans. The African diversity estimate is even higher than that between Africans and Eurasians (0.096% +/- 0.012%). From available data for noncoding autosomal regions (total length = 47,038 bp) and X-linked regions (47,421 bp), we estimated the pi-values for autosomal regions to be 0.105, 0.070, 0.069, and 0.097% for Africans, Asians, Europeans, and between Africans and Eurasians, and the corresponding values for X-linked regions to be 0.088, 0.042, 0.053, and 0.082%. Thus, Africans differ from one another slightly more than from Eurasians, and the genetic diversity in Eurasians is largely a subset of that in Africans, supporting the out of Africa model of human evolution. Clearly, one must specify the geographic origins of the individuals sampled when studying pi or SNP density.     

Lineage‐specific adaptation to climate involves flowering time in North American Arabidopsis lyrata

Adaptation to local climatic conditions is commonly found within species, but whether it involves the same intraspecific genomic variants is unknown. We studied this question in North American Arabidopsis lyrata, whose current distribution is shaped by post‐glacial range expansion from two refugia, resulting in two distinct genetic clusters covering comparable climatic gradients. Using pooled whole‐genome sequence data of 41 outcrossing populations, we identified loci associated with three niche‐determining climatic variables in the two clusters and compared these outliers. Little evidence was found for parallelism in climate adaptation for single nucleotide polymorphisms (SNPs) and for genes with an accumulation of outlier SNPs. Significantly increased selection coefficients supported them as candidates of climate adaptation. However, the fraction of gene ontology (GO) terms shared between clusters was higher compared to outlier SNPs and outlier genes, suggesting that selection acts on similar pathways but not necessarily the same genes. Enriched GO terms involved responses to abiotic and biotic stress, circadian rhythm and development, with flower development and reproduction being among the most frequently detected. In line with GO enrichment, regulators of flowering time were detected as outlier genes. Our results suggest that while adaptation to environmental gradients on the genomic level are lineage‐specific in A. lyrata, similar biological processes seem to be involved. Differential loss of standing genetic variation, probably driven by genetic drift, can in part account for the lack of parallel evolution on the genomic level.     

Large Genomic Region Free of GWAS-Based Common Variants Contains Fertility-Related Genes

DNA variants, such as single nucleotide polymorphisms (SNPs) and copy number variants (CNVs), are unevenly distributed across the human genome. Currently, dbSNP contains more than 6 million human SNPs, and whole-genome genotyping arrays can assay more than 4 million of them simultaneously. In our study, we first questioned whether published genome-wide association studies (GWASs) assays cover all regions well in the genome. Using dbSNP build 135 data, we identified 50 genomic regions longer than 100 Kb that do not contain any common SNPs, i.e., those with minor allele frequency (MAF)≥1%. Secondly, because conserved regions are generally of functional importance, we tested genes in those large genomic regions without common SNPs. We found 97 genes and were enriched for reproduction function. In addition, we further filtered out regions with CNVs listed in the Database of Genomic Variants (DGV), segmental duplications from Human Genome Project and common variants identified by personal genome sequencing (UCSC). No region survived after those filtering. Our analysis suggests that, while there may not be many large genomic regions free of common variants, there are still some “holes” in the current human genomic map for common SNPs. Because GWAS only focused on common SNPs, interpretation of GWAS results should take this limitation into account. Particularly, two recent GWAS of fertility may be incomplete due to the map deficit. Additional SNP discovery efforts should pay close attention to these regions.     

Linkage-disequilibrium-based binning affects the interpretation of GWASs

Genome-wide association studies (GWASs) are critically dependent on detailed knowledge of the pattern of linkage disequilibrium (LD) in the human genome. GWASs generate lists of variants, usually SNPs, ranked according to the significance of their association to a trait. Downstream analyses generally focus on the gene or genes that are physically closest to these SNPs and ignore their LD profile with other SNPs. We have developed a flexible R package (LDsnpR) that efficiently assigns SNPs to genes on the basis of both their physical position and their pairwise LD with other SNPs. We used the positional-binning and LD-based-binning approaches to investigate whether including these "LD-based" SNPs would affect the interpretation of three published GWASs on bipolar affective disorder (BP) and of the imputed versions of two of these GWASs. We show how including LD can be important for interpreting and comparing GWASs. In the published, unimputed GWASs, LD-based binning effectively "recovered" 6.1%-8.3% of Ensembl-defined genes. It altered the ranks of the genes and resulted in nonnegligible differences between the lists of the top 2,000 genes emerging from the two binning approaches. It also improved the overall gene-based concordance between independent BP studies. In the imputed datasets, although the increases in coverage (>0.4%) and rank changes were more modest, even greater concordance between the studies was observed, attesting to the potential of LD-based binning on imputed data as well. Thus, ignoring LD can result in the misinterpretation of the GWAS findings and have an impact on subsequent genetic and functional studies.     

Whole‐genome resequencing uncovers molecular signatures of natural and sexual selection in wild bighorn sheep

The identification of genes influencing fitness is central to our understanding of the genetic basis of adaptation and how it shapes phenotypic variation in wild populations. Here, we used whole‐genome resequencing of wild Rocky Mountain bighorn sheep (Ovis canadensis) to >50‐fold coverage to identify 2.8 million single nucleotide polymorphisms (SNPs) and genomic regions bearing signatures of directional selection (i.e. selective sweeps). A comparison of SNP diversity between the X chromosome and the autosomes indicated that bighorn males had a dramatically reduced long‐term effective population size compared to females. This probably reflects a long history of intense sexual selection mediated by male–male competition for mates. Selective sweep scans based on heterozygosity and nucleotide diversity revealed evidence for a selective sweep shared across multiple populations at RXFP2, a gene that strongly affects horn size in domestic ungulates. The massive horns carried by bighorn rams appear to have evolved in part via strong positive selection at RXFP2. We identified evidence for selection within individual populations at genes affecting early body growth and cellular response to hypoxia; however, these must be interpreted more cautiously as genetic drift is strong within local populations and may have caused false positives. These results represent a rare example of strong genomic signatures of selection identified at genes with known function in wild populations of a nonmodel species. Our results also showcase the value of reference genome assemblies from agricultural or model species for studies of the genomic basis of adaptation in closely related wild taxa.     

Copy number variants and common disorders: filling the gaps and exploring complexity in genome-wide association studies

Genome-wide association scans (GWASs) using single nucleotide polymorphisms (SNPs) have been completed successfully for several common disorders and have detected over 30 new associations. Considering the large sample sizes and genome-wide SNP coverage of the scans, one might have expected many of the common variants underpinning the genetic component of various disorders to have been identified by now. However, these studies have not evaluated the contribution of other forms of genetic variation, such as structural variation, mainly in the form of copy number variants (CNVs). Known CNVs account for over 15% of the assembled human genome sequence. Since CNVs are not easily tagged by SNPs, might have a wide range of copy number variability, and often fall in genomic regions not well covered by whole-genome arrays or not genotyped by the HapMap project, current GWASs have largely missed the contribution of CNVs to complex disorders. In fact, some CNVs have already been reported to show association with several complex disorders using candidate gene/region approaches, underpinning the importance of regions not investigated in current GWASs. This reveals the need for new generation arrays (some already in the market) and the use of tailored approaches to explore the full dimension of genome variability beyond the single nucleotide scale.     

A global barley panel revealing genomic signatures of breeding in modern Australian cultivars

The future of plant cultivar improvement lies in the evaluation of genetic resources from currently available germplasm. Today’s gene pool of crop genetic diversity has been shaped during domestication and more recently by breeding. Recent efforts in plant breeding have been aimed at developing new and improved varieties from poorly adapted crops to suit local environments. However, the impact of these breeding efforts is poorly understood. Here, we assess the contributions of both historical and recent breeding efforts to local adaptation and crop improvement in a global barley panel by analysing the distribution of genetic variants with respect to geographic region or historical breeding category. By tracing the impact that breeding had on the genetic diversity of Hordeum vulgare (barley) released in Australia, where the history of barley production is relatively young, we identify 69 candidate regions within 922 genes that were under selection pressure. We also show that modern Australian barley varieties exhibit 12% higher genetic diversity than historical cultivars. Finally, field-trialling and phenotyping for agriculturally relevant traits across a diverse range of Australian environments suggests that genomic regions under strong breeding selection and their candidate genes are closely associated with key agronomic traits. In conclusion, our combined data set and germplasm collection provide a rich source of genetic diversity that can be applied to understanding and improving environmental adaptation and enhanced yields.     

Adaptive evolution of the Populus tremula photoperiod pathway

Perennial plants monitor seasonal changes through changes in environmental conditions such as the quantity and quality of light and genes in the photoperiodic pathway are known to be involved in controlling these processes. Here, we examine 25 of genes from the photoperiod pathway in Populus tremula (Salicaceae) for signatures of adaptive evolution. Overall, levels of synonymous polymorphism in the 25 genes are lower than at control loci selected randomly from the genome. This appears primarily to be caused by lower levels of synonymous polymorphism in genes associated with the circadian clock. Natural selection appears to play an important role in shaping protein evolution at several of the genes in the photoperiod pathways, which is highlighted by the fact that approximately 40% of the genes from the photoperiod pathway have estimates of selection on nonsynonymous polymorphisms that are significantly different from zero. A surprising observation we make is that circadian clock-associated genes appear to be over-represented among the genes showing elevated rates of protein evolution; seven genes are evolving under positive selection and all but one of these genes are involved in the circadian clock of Populus.     

Variation in the human TAS1R taste receptor genes

We have performed a comprehensive evaluation of single-nucleotide polymorphisms (SNPs) and haplotypes in the human TAS1R gene family, which encodes receptors for sweet and umami tastes. Complete DNA sequences of TAS1R1-, TAS1R2-, and TAS1R3-coding regions, obtained from 88 individuals of African, Asian, European, and Native American origin, revealed substantial coding and noncoding diversity: polymorphisms are common in these genes, and polymorphic sites and SNP frequencies vary widely in human populations. The genes TAS1R1 and TAS1R3, which encode proteins that act as a dimer to form the umami (glutamate) taste receptor, showed less variation than the TAS1R2 gene, which acts as a dimer with TAS1R3 to form the sweet taste receptor. The TAS1R3 gene, which encodes a subunit common to both the sweet and umami receptors, was the most conserved. Evolutionary genetic analysis indicates that these variants have come to their current frequencies under natural selection during population growth and support the view that the coding sequence variants affect receptor function. We propose that human populations likely vary little with respect to umami perception, which is controlled by one major form of the receptor that is optimized for detecting glutamate but may vary much more with respect to sweet perception.