Complex nature of SNP genotype effects on gene expression in primary human leucocytes

Meta-analyses identify 13 loci associated with age at menopause and highlight DNA repair and immune pathways

To newly identify loci for age at natural menopause, we carried out a meta-analysis of 22 genome-wide association studies (GWAS) in 38,968 women of European descent, with replication in up to 14,435 women. In addition to four known loci, we identified 13 loci newly associated with age at natural menopause (at P < 5 × 10(-8)). Candidate genes located at these newly associated loci include genes implicated in DNA repair (EXO1, HELQ, UIMC1, FAM175A, FANCI, TLK1, POLG and PRIM1) and immune function (IL11, NLRP11 and PRRC2A (also known as BAT2)). Gene-set enrichment pathway analyses using the full GWAS data set identified exoDNase, NF-κB signaling and mitochondrial dysfunction as biological processes related to timing of menopause.

Unraveling the regulatory mechanisms underlying tissue-dependent genetic variation of gene expression

It is known that genetic variants can affect gene expression, but it is not yet completely clear through what mechanisms genetic variation mediate this expression. We therefore compared the cis-effect of single nucleotide polymorphisms (SNPs) on gene expression between blood samples from 1,240 human subjects and four primary non-blood tissues (liver, subcutaneous, and visceral adipose tissue and skeletal muscle) from 85 subjects. We characterized four different mechanisms for 2,072 probes that show tissue-dependent genetic regulation between blood and non-blood tissues: on average 33.2% only showed cis-regulation in non-blood tissues; 14.5% of the eQTL probes were regulated by different, independent SNPs depending on the tissue of investigation. 47.9% showed a different effect size although they were regulated by the same SNPs. Surprisingly, we observed that 4.4% were regulated by the same SNP but with opposite allelic direction. We show here that SNPs that are located in transcriptional regulatory elements are enriched for tissue-dependent regulation, including SNPs at 3' and 5' untranslated regions (P = 1.84×10(-5) and 4.7×10(-4), respectively) and SNPs that are synonymous-coding (P = 9.9×10(-4)). SNPs that are associated with complex traits more often exert a tissue-dependent effect on gene expression (P = 2.6×10(-10)). Our study yields new insights into the genetic basis of tissue-dependent expression and suggests that complex trait associated genetic variants have even more complex regulatory effects than previously anticipated.

Feature Selection with Conjunctions of Decision Stumps and Learning from Microarray Data

One of the objectives of designing feature selection learning algorithms is to obtain classifiers that depend on a small number of attributes and have verifiable future performance guarantees. There are few, if any, approaches that successfully address the two goals simultaneously. To the best of our knowledge, such algorithms that give theoretical bounds on the future performance have not been proposed so far in the context of the classification of gene expression data. In this work, we investigate the premise of learning a conjunction (or disjunction) of decision stumps in Occam's Razor, Sample Compression, and PAC-Bayes learning settings for identifying a small subset of attributes that can be used to perform reliable classification tasks. We apply the proposed approaches for gene identification from DNA microarray data and compare our results to those of the well-known successful approaches proposed for the task. We show that our algorithm not only finds hypotheses with a much smaller number of genes while giving competitive classification accuracy but also having tight risk guarantees on future performance, unlike other approaches. The proposed approaches are general and extensible in terms of both designing novel algorithms and application to other domains.

Comparing methods for mapping cis acting polymorphisms using allelic expression ratios

Genome wide association studies frequently reveal associations between disease susceptibility and polymorphisms outside coding regions. Such associations cannot always be explained by linkage disequilibrium with changes affecting the transcription products. This has stimulated the interest in characterising sequence variation influencing gene expression levels, in particular in changes acting in cis. Differences in transcription between the two alleles at an autosomal locus can be used to test the association between candidate polymorphisms and the modulation of gene expression in cis. This type of approach requires at least one transcribed polymorphism and one candidate polymorphism. In the past five years, different methods have been proposed to analyse such data. Here we use simulations and real data sets to compare the power of some of these methods. The results show that when it is not possible to determine the phase between the transcribed and potentially cis acting allele there is some advantage in using methods that estimate phased genotype and effect on expression simultaneously. However when the phase can be determined, simple regression models seem preferable because of their simplicity and flexibility. The simulations and the analysis of experimental data suggest that in the majority of situations, methods that assume a lognormal distribution of the allelic expression ratios are both robust to deviations from this assumption and more powerful than alternatives that do not make these assumptions.

Haplotype in the IBD5 region is associated with refractory Crohn's disease in Slovenian patients and modulates expression of the SLC22A5 gene

BACKGROUND

The IBD5 locus (OMIM ID 606348) on chromosome 5 was suggested to be one of the most important genetic factors involved in the pathogenesis of inflammatory bowel diseases (IBDs). However the main contributor from this region is still unknown.

METHODS

We investigated the possible association of the IBD5 locus with IBD in Slovenian patients and correlation between disease-associated single nucleotide polymorphisms (SNPs) and quantitative gene expression (eQTL) of candidate genes from the IBD5 locus in peripheral blood lymphocytes and colon tissue biopsies from IBD patients. We genotyped SNPs from the IBD5 locus in 312 healthy controls and 632 IBD patients.

RESULTS

We found statistically significant association of polymorphisms rs1050152 in gene SLC22A4 (p = 0.005, OR = 2.177, 95% CI = 1.270-3.526) and rs2631372 in gene SLC22A5 (p = 0.001, OR = 0.473, 95% CI = 0.307-0.731) and TC haplotype of both polymorphisms (p = 0.006, OR = 1,541, 95% CI = 1.130-2.100) with refractory Crohn's disease (CD) in Slovenian patients who do not respond to standard therapy, including patients who develop fistulas. We found decreased expression of SLC22A4 and SLC22A5 genes in peripheral blood lymphocytes from IBD patients compared to control group and decreased expression of SLC22A5 gene in inflamed tissue biopsies compared to noninflamed colon (p = 0.009). We found lower expression of SLC22A5 gene in IBD patients with disease-susceptible genotypes for both disease-associated SNPs.

CONCLUSIONS

Our data suggest that SNPs and haplotype in the IBD5 SLC22A4/SLC22A5 region contribute to the development of particularly refractory Crohn's disease in the Slovenian population, and expression studies in blood lymphocytes and colon tissue biopsies and eQTL analysis suggest that SLC22A5 is the main gene in the IBD5 region contributing to the IBD pathogenesis.

Trans-eQTLs reveal that independent genetic variants associated with a complex phenotype converge on intermediate genes, with a major role for the HLA

For many complex traits, genetic variants have been found associated. However, it is still mostly unclear through which downstream mechanism these variants cause these phenotypes. Knowledge of these intermediate steps is crucial to understand pathogenesis, while also providing leads for potential pharmacological intervention. Here we relied upon natural human genetic variation to identify effects of these variants on trans-gene expression (expression quantitative trait locus mapping, eQTL) in whole peripheral blood from 1,469 unrelated individuals. We looked at 1,167 published trait- or disease-associated SNPs and observed trans-eQTL effects on 113 different genes, of which we replicated 46 in monocytes of 1,490 different individuals and 18 in a smaller dataset that comprised subcutaneous adipose, visceral adipose, liver tissue, and muscle tissue. HLA single-nucleotide polymorphisms (SNPs) were 10-fold enriched for trans-eQTLs: 48% of the trans-acting SNPs map within the HLA, including ulcerative colitis susceptibility variants that affect plausible candidate genes AOAH and TRBV18 in trans. We identified 18 pairs of unlinked SNPs associated with the same phenotype and affecting expression of the same trans-gene (21 times more than expected, P<10⁻¹⁶). This was particularly pronounced for mean platelet volume (MPV): Two independent SNPs significantly affect the well-known blood coagulation genes GP9 and F13A1 but also C19orf33, SAMD14, VCL, and GNG11. Several of these SNPs have a substantially higher effect on the downstream trans-genes than on the eventual phenotypes, supporting the concept that the effects of these SNPs on expression seems to be much less multifactorial. Therefore, these trans-eQTLs could well represent some of the intermediate genes that connect genetic variants with their eventual complex phenotypic outcomes.

Many genetic variants have been found associated with diseases. However, for many of these genetic variants, it remains unclear how they exert their effect on the eventual phenotype. We investigated genetic variants that are known to be associated with diseases and complex phenotypes and assessed whether these variants were also associated with gene expression levels in a set of 1,469 unrelated whole blood samples. For several diseases, such as type 1 diabetes and ulcerative colitis, we observed that genetic variants affect the expression of genes, not implicated before. For complex traits, such as mean platelet volume and mean corpuscular volume, we observed that independent genetic variants on different chromosomes influence the expression of exactly the same genes. For mean platelet volume, these genes include well-known blood coagulation genes but also genes with still unknown functions. These results indicate that, by systematically correlating genetic variation with gene expression levels, it is possible to identify downstream genes, which provide important avenues for further research.

Insights from GWAS into the quantitative genetics of transcription in humans

Human gene expression profiles have emerged as an effective model system for the dissection of quantitative genetic traits. Peripheral blood and transformed lymphoblasts are particularly attractive for their ready availability and repeatability, respectively, and the advent of relatively inexpensive genotyping and microarray analysis technologies has facilitated genome-wide association for transcript abundance in numerous settings. Thousands of genes have been shown to harbour regulatory polymorphisms that have large local effects on transcription, explaining 20% or more of the variance in many cases, but the focus on such results obscures the reality that the vast majority of the genetic component of transcriptional variance remains to be ascertained. This mini-review surveys the inferences derived from genome-wide association studies (GWAS) for gene expression to date, and discusses some of the issues we face in finding the remainder of the heritability and understanding how environmental and genetic regulatory factors orchestrate the highly structured architecture of transcriptional variation.

MixupMapper: correcting sample mix-ups in genome-wide datasets increases power to detect small genetic effects

MOTIVATION

Sample mix-ups can arise during sample collection, handling, genotyping or data management. It is unclear how often sample mix-ups occur in genome-wide studies, as there currently are no post hoc methods that can identify these mix-ups in unrelated samples. We have therefore developed an algorithm (MixupMapper) that can both detect and correct sample mix-ups in genome-wide studies that study gene expression levels.

RESULTS

We applied MixupMapper to five publicly available human genetical genomics datasets. On average, 3% of all analyzed samples had been assigned incorrect expression phenotypes: in one of the datasets 23% of the samples had incorrect expression phenotypes. The consequences of sample mix-ups are substantial: when we corrected these sample mix-ups, we identified on average 15% more significant cis-expression quantitative trait loci (cis-eQTLs). In one dataset, we identified three times as many significant cis-eQTLs after correction. Furthermore, we show through simulations that sample mix-ups can lead to an underestimation of the explained heritability of complex traits in genome-wide association datasets.

AVAILABILITY AND IMPLEMENTATION

MixupMapper is freely available at http://www.genenetwork.nl/mixupmapper/

Peripheral blood gene expression profiles in metabolic syndrome, coronary artery disease and type 2 diabetes

To determine if individuals with metabolic disorders possess unique gene expression profiles, we compared transcript levels in peripheral blood from patients with coronary artery disease, type 2 diabetes and their precursor state, metabolic syndrome to those of control subjects and subjects with rheumatoid arthritis. The gene expression profile of each metabolic state was distinguishable from controls and correlated with other metabolic states more than with rheumatoid arthritis. Of note, subjects in the metabolic cohorts over-expressed gene sets that participate in the innate immune response. Genes involved in activation of the pro-inflammatory transcription factor, NF-κB, were over-expressed in coronary artery disease while genes differentially expressed in type 2 diabetes play key roles in T cell activation and signaling. RT-PCR validation confirmed microarray results. Furthermore, several genes differentially expressed in human metabolic disorders have been previously shown to participate in inflammatory responses in murine models of obesity and Type 2 diabetes. Taken together, these data demonstrate that peripheral blood from individuals with metabolic disorders display overlapping and non-overlapping patterns of gene expression indicative of unique, underlying immune processes.

Gene expression and genetic variation data implicate PCLO in bipolar disorder

BACKGROUND

Genetic variation may contribute to differential gene expression in the brain of individuals with psychiatric disorders. To test this hypothesis, we identified genes that were differentially expressed in individuals with bipolar disorder, along with nearby single nucleotide polymorphisms (SNPs) that were associated with expression of the same genes. We then tested these SNPs for association with bipolar disorder in large case-control samples.

METHODS

We used the Stanley Genomics Database to extract gene expression and SNP microarray data from individuals with bipolar disorder (n = 40) and unaffected controls (n = 43). We identified 367 genes that were differentially expressed in the prefrontal cortex of cases vs. controls (fold change > 1.3 and FDR q-value < .05) and 45 nearby SNPs that were associated with expression of those same genes (FDR q-value < .05). We tested these SNPs for association with bipolar disorder in a meta-analysis of genome-wide association studies (GWAS) including 4,936 cases and 6,654 healthy controls.

RESULTS

We identified 45 SNPs that were associated with expression of differentially expressed genes, including HBS1L (15 SNPs), HLA-DPB1 (15 SNPs), AMFR (8 SNPs), PCLO (2 SNPs) and WDR41 (2 SNPs). Of these, one SNP (rs13438494), in an intron of the piccolo (PCLO) gene, was significantly associated with bipolar disorder (FDR adjusted p < .05) in the meta-analysis of GWAS.

CONCLUSIONS

These results support the previous findings implicating PCLO in mood disorders and demonstrate the utility of combining gene expression and genetic variation data to improve our understanding of the genetic contribution to bipolar disorder.

Specific T cell frequency and cytokine expression profile do not correlate with protection against tuberculosis after bacillus Calmette-Guérin vaccination of newborns

RATIONALE

Immunogenicity of new tuberculosis (TB) vaccines is commonly assessed by measuring the frequency and cytokine expression profile of T cells.

OBJECTIVES

We tested whether this outcome correlates with protection against childhood TB disease after newborn vaccination with bacillus Calmette-Guérin (BCG).

METHODS

Whole blood from 10-week-old infants, routinely vaccinated with BCG at birth, was incubated with BCG for 12 hours, followed by cryopreservation for intracellular cytokine analysis. Infants were followed for 2 years to identify those who developed culture-positive TB-these infants were regarded as not protected against TB. Infants who did not develop TB disease despite exposure to TB in the household, and another group of randomly selected infants who were never evaluated for TB, were also identified-these groups were regarded as protected against TB. Cells from these groups were thawed, and CD4, CD8, and γδ T cell-specific expression of IFN-γ, TNF-α, IL-2, and IL-17 measured by flow cytometry.

MEASUREMENTS AND MAIN RESULTS

A total of 5,662 infants were enrolled; 29 unprotected and two groups of 55 protected infants were identified. There was no difference in frequencies of BCG-specific CD4, CD8, and γδ T cells between the three groups of infants. Although BCG induced complex patterns of intracellular cytokine expression, there were no differences between protected and unprotected infants.

CONCLUSIONS

The frequency and cytokine profile of mycobacteria-specific T cells did not correlate with protection against TB. Critical components of immunity against Mycobacterium tuberculosis, such as CD4 T cell IFN-γ production, may not necessarily translate into immune correlates of protection against TB disease.

Multiple common variants for celiac disease influencing immune gene expression

We performed a second-generation genome-wide association study of 4,533 individuals with celiac disease (cases) and 10,750 control subjects. We genotyped 113 selected SNPs with P(GWAS) < 10(-4) and 18 SNPs from 14 known loci in a further 4,918 cases and 5,684 controls. Variants from 13 new regions reached genome-wide significance (P(combined) < 5 x 10(-8)); most contain genes with immune functions (BACH2, CCR4, CD80, CIITA-SOCS1-CLEC16A, ICOSLG and ZMIZ1), with ETS1, RUNX3, THEMIS and TNFRSF14 having key roles in thymic T-cell selection. There was evidence to suggest associations for a further 13 regions. In an expression quantitative trait meta-analysis of 1,469 whole blood samples, 20 of 38 (52.6%) tested loci had celiac risk variants correlated (P < 0.0028, FDR 5%) with cis gene expression.

XGAP: a uniform and extensible data model and software platform for genotype and phenotype experiments

XGAP, a software platform for the integration and analysis of genotype and phenotype data.

We present an extensible software model for the genotype and phenotype community, XGAP. Readers can download a standard XGAP (http://www.xgap.org) or auto-generate a custom version using MOLGENIS with programming interfaces to R-software and web-services or user interfaces for biologists. XGAP has simple load formats for any type of genotype, epigenotype, transcript, protein, metabolite or other phenotype data. Current functionality includes tools ranging from eQTL analysis in mouse to genome-wide association studies in humans.

Genome-wide analysis of allelic expression imbalance in human primary cells by high-throughput transcriptome resequencing

Many disease-associated variants identified by genome-wide association (GWA) studies are expected to regulate gene expression. Allele-specific expression (ASE) quantifies transcription from both haplotypes using individuals heterozygous at tested SNPs. We performed deep human transcriptome-wide resequencing (RNA-seq) for ASE analysis and expression quantitative trait locus discovery. We resequenced double poly(A)-selected RNA from primary CD4(+) T cells (n = 4 individuals, both activated and untreated conditions) and developed tools for paired-end RNA-seq alignment and ASE analysis. We generated an average of 20 million uniquely mapping 45 base reads per sample. We obtained sufficient read depth to test 1371 unique transcripts for ASE. Multiple biases inflate the false discovery rate which we estimate to be approximately 50% for random SNPs. However, after controlling for these biases and considering the subset of SNPs that pass HapMap QC, 4.6% of heterozygous SNP-sample pairs show evidence of imbalance (P < 0.001). We validated four findings by both bacterial cloning and Sanger sequencing assays. We also found convincing evidence for allelic imbalance at multiple reporter exonic SNPs in CD6 for two samples heterozygous at the multiple sclerosis-associated variant rs17824933, linking GWA findings with variation in gene expression. Finally, we show in CD4(+) T cells from a further individual that high-throughput sequencing of genomic DNA and RNA-seq following enrichment for targeted gene sequences by sequence capture methods offers an unbiased means to increase the read depth for transcripts of interest, and therefore a method to investigate the regulatory role of many disease-associated genetic variants.

Modeling expression quantitative trait loci in data combining ethnic populations

Background

Combining data from different ethnic populations in a study can increase efficacy of methods designed to identify expression quantitative trait loci (eQTL) compared to analyzing each population independently. In such studies, however, the genetic diversity of minor allele frequencies among populations has rarely been taken into account. Due to the fact that allele frequency diversity and population-level expression differences are present in populations, a consensus regarding the optimal statistical approach for analysis of eQTL in data combining different populations remains inconclusive.

Results

In this report, we explored the applicability of a constrained two-way model to identify eQTL for combined ethnic data that might contain genetic diversity among ethnic populations. In addition, gene expression differences resulted from ethnic allele frequency diversity between populations were directly estimated and analyzed by the constrained two-way model. Through simulation, we investigated effects of genetic diversity on eQTL identification by examining gene expression data pooled from normal quantile transformation of each population. Using the constrained two-way model to reanalyze data from Caucasians and Asian individuals available from HapMap, a large number of eQTL were identified with similar genetic effects on the gene expression levels in these two populations. Furthermore, 19 single nucleotide polymorphisms with inter-population differences with respect to both genotype frequency and gene expression levels directed by genotypes were identified and reflected a clear distinction between Caucasians and Asian individuals.

Conclusions

This study illustrates the influence of minor allele frequencies on common eQTL identification using either separate or combined population data. Our findings are important for future eQTL studies in which different datasets are combined to increase the power of eQTL identification.

eQTL analysis in humans

Improving human health is a major aim of medical research, but it requires that variation between individuals be taken into account since each person carries a different combination of gene variants and is exposed to different environmental conditions, which can cause differences in susceptibility to diseases. With the advent of molecular markers in the 1980s, it became possible to genotype individuals (i.e., to detect the presence or absence of local DNA sequence variants at each of hundreds of genome positions). This DNA sequence variation could then be related to disease susceptibility by using pedigree data. Such linkage analyses proved to be difficult for more complex diseases. Recently, with the decreasing costs of genotyping, analyses of large natural populations of unrelated individuals became possible and resulted in the association of many genes (and genetic variants in these genes) with complex diseases. Unfortunately, for a considerable proportion of these genes and their proteins, it is not yet clear what their downstream effects are. Studying the expression of these genes and proteins can help to uncover the effects of these variants on the expression of these and other genes, proteins, metabolites, and phenotypes. In this chapter, we focus on the high-throughput and genome-wide measurement of gene expression in a natural population of unrelated humans, and on the subsequent association of variation in expression to "expression quantitative trait loci" (eQTLs) on DNA using oligonucleotide arrays with hundreds of thousands of single-nucleotide polymorphism (SNP) markers that capture most of the human genetic variation well. This strategy has been successfully applied to several diseases such as celiac disease (Hunt et al. 2008, Nat Genet 40, 395-402) and asthma (Moffatt et al. 2007, Nature 448, 470-473): associated genetic variants have been identified that affect levels of gene expression in cis or in trans, providing insight into the biological pathways affected by these diseases.

Geographical genomics of human leukocyte gene expression variation in southern Morocco

Studies of the genetics of gene expression can identify expression SNPs (eSNPs) that explain variation in transcript abundance. Here we address the robustness of eSNP associations to environmental geography and population structure in a comparison of 194 Arab and Amazigh individuals from a city and two villages in southern Morocco. Gene expression differed between pairs of locations for up to a third of all transcripts, with notable enrichment of transcripts involved in ribosomal biosynthesis and oxidative phosphorylation. Robust associations were observed in the leukocyte samples: cis eSNPs (P < 10(-08)) were identified for 346 genes, and trans eSNPs (P < 10(-11)) for 10 genes. All of these associations were consistent both across the three sample locations and after controlling for ancestry and relatedness. No evidence of large-effect trans-acting mediators of the pervasive environmental influence was found; instead, genetic and environmental factors acted in a largely additive manner.