Cell Specific eQTL Analysis without Sorting Cells.

The functional consequences of trait associated SNPs are often investigated using expression quantitative trait locus (eQTL) mapping. While trait-associated variants may operate in a cell-type specific manner, eQTL datasets for such cell-types may not always be available. We performed a genome-environment interaction (GxE) meta-analysis on data from 5,683 samples to infer the cell type specificity of whole blood cis-eQTLs. We demonstrate that this method is able to predict neutrophil and lymphocyte specific cis-eQTLs and replicate these predictions in independent cell-type specific datasets. Finally, we show that SNPs associated with Crohn's disease preferentially affect gene expression within neutrophils, including the archetypal NOD2 locus.

A meta-analysis of gene expression signatures of blood pressure and hypertension.

Genome-wide association studies (GWAS) have uncovered numerous genetic variants (SNPs) that are associated with blood pressure (BP). Genetic variants may lead to BP changes by acting on intermediate molecular phenotypes such as coded protein sequence or gene expression, which in turn affect BP variability. Therefore, characterizing genes whose expression is associated with BP may reveal cellular processes involved in BP regulation and uncover how transcripts mediate genetic and environmental effects on BP variability. A meta-analysis of results from six studies of global gene expression profiles of BP and hypertension in whole blood was performed in 7017 individuals who were not receiving antihypertensive drug treatment. We identified 34 genes that were differentially expressed in relation to BP (Bonferroni-corrected p<0.05). Among these genes, FOS and PTGS2 have been previously reported to be involved in BP-related processes; the others are novel. The top BP signature genes in aggregate explain 5%-9% of inter-individual variance in BP. Of note, rs3184504 in SH2B3, which was also reported in GWAS to be associated with BP, was found to be a trans regulator of the expression of 6 of the transcripts we found to be associated with BP (FOS, MYADM, PP1R15A, TAGAP, S100A10, and FGBP2). Gene set enrichment analysis suggested that the BP-related global gene expression changes include genes involved in inflammatory response and apoptosis pathways. Our study provides new insights into molecular mechanisms underlying BP regulation, and suggests novel transcriptomic markers for the treatment and prevention of hypertension.

A common 16p11.2 inversion underlies the joint susceptibility to asthma and obesity.

The prevalence of asthma and obesity is increasing worldwide, and obesity is a well-documented risk factor for asthma. The mechanisms underlying this association and parallel time trends remain largely unknown but genetic factors may be involved. Here, we report on a common ~0.45 Mb genomic inversion at 16p11.2 that can be accurately genotyped via SNP array data. We show that the inversion allele protects against the joint occurrence of asthma and obesity in five large independent studies (combined sample size of 317 cases and 543 controls drawn from a total of 5,809 samples; combined OR = 0.48, p = 5.5 × 10(-6)). Allele frequencies show remarkable worldwide population stratification, ranging from 10% in East Africa to 49% in Northern Europe, consistent with discordant and extreme genetic drifts or adaptive selections after human migration out of Africa. Inversion alleles strongly correlate with expression levels of neighboring genes, especially TUFM (p = 3.0 × 10(-40)) that encodes a mitochondrial protein regulator of energy balance and inhibitor of type 1 interferon, and other candidates for asthma (IL27) and obesity (APOB48R and SH2B1). Therefore, by affecting gene expression, the ~0.45 Mb 16p11.2 inversion provides a genetic basis for the joint susceptibility to asthma and obesity, with a population attributable risk of 39.7%. Differential mitochondrial function and basal energy balance of inversion alleles might also underlie the potential selection signature that led to their uneven distribution in world populations.

Human basonuclin 2 up-regulates a cascade set of interferon-stimulated genes with anti-cancerous properties in a lung cancer model.

BACKGROUND: Human basonuclin 2 (BNC2) acts as a tumor suppressor in multiple cancers in an as yet unidentified manner. The role and expression of the BNC2 gene in lung cancer has not yet been investigated. METHODS: BNC2 expression was studied in the A549 and BEAS-2B cell lines, as well as in lung cancer tissue. Illumina array analysis and a viability assay were used to study the effects of transient transfection of BNC2 in A549 cells. Ingenuity pathway analysis and g:Profiler were applied to identify affected pathways and networks. RT-qPCR was used to validate the array results. RESULTS: We showed the reduced mRNA expression of BNC2 in non-small cell lung cancer tissue and lung cancer cell line A549 compared to non-cancerous lung tissue and BEAS-2B cells, respectively. Further array analysis demonstrated that the transfection of BNC2 into A549 cells resulted in the increased expression of 139 genes and the down-regulation of 13 genes. Pathway analysis revealed that half of the up-regulated genes were from the interferon/signal transducer and activator of transcription signaling pathways. The differential expression of selected sets of genes, including interferon-stimulated and tumor suppressor genes of the XAF1 and OAS families, was confirmed by RT-qPCR. In addition, we showed that the over-expression of BNC2 inhibited the proliferation of A549 cells. CONCLUSION: Our data suggest that human BNC2 is an activator of a subset of IFN-regulated genes and might thereby act as a tumor suppressor.

Human disease-associated genetic variation impacts large intergenic non-coding RNA expression.

Recently it has become clear that only a small percentage (7%) of disease-associated single nucleotide polymorphisms (SNPs) are located in protein-coding regions, while the remaining 93% are located in gene regulatory regions or in intergenic regions. Thus, the understanding of how genetic variations control the expression of non-coding RNAs (in a tissue-dependent manner) has far-reaching implications. We tested the association of SNPs with expression levels (eQTLs) of large intergenic non-coding RNAs (lincRNAs), using genome-wide gene expression and genotype data from five different tissues. We identified 112 cis-regulated lincRNAs, of which 45% could be replicated in an independent dataset. We observed that 75% of the SNPs affecting lincRNA expression (lincRNA cis-eQTLs) were specific to lincRNA alone and did not affect the expression of neighboring protein-coding genes. We show that this specific genotype-lincRNA expression correlation is tissue-dependent and that many of these lincRNA cis-eQTL SNPs are also associated with complex traits and diseases.

Polymorphisms in ESR1, ESR2 and HSD17B1 genes are associated with fertility status in endometriosis.

OBJECTIVE: To investigate whether polymorphisms in genes involved in biosynthesis and signalling of sex steroids influence susceptibility to endometriosis and to infertility associated with it. MATERIALS AND METHODS: Patients with endometriosis (n = 150) and fertile controls (n = 199) were genotyped for polymorphisms in oestrogen receptor genes ESR1 (rs2234693 - T/C single nucleotide polymorphism (SNP), dinucleotide (TA)(n) repeat) and ESR2 (dinucleotide (CA)(n) repeat), progesterone receptor gene PGR (rs10895068 - G/A SNP, 306-bp Alu-insertion), 17ß-hydroxysteroid dehydrogenase type 1 gene HSD17B1 (rs605059 - A/G SNP), and aromatase gene CYP19A1 (rs10046 - C/T SNP, (TTTA)(n) tetranucleotide repeat, 3-bp TCT insertion/deletion polymorphism). RESULTS: The HSD17B1 A/G SNP A allele increased overall endometriosis risk and the risk of stage I-II disease, while ESR1 longer (TA)(n) repeats only correlated with susceptibility to stage I-II endometriosis. When considering patients' fertility status, HSD17B1 A/G SNP A allele and ESR1 longer (TA)(n) repeats were associated with endometriosis accompanied by infertility, while ESR2 shorter (CA)(n) repeats were linked with endometriosis without infertility. Other polymorphisms were distributed similarly among patients and controls. CONCLUSIONS: Genetic variants in ESR1, ESR2, and HSD17B1 genes could modify susceptibility to endometriosis and might influence the fertility status in endometriosis patients.

Extensive alterations of the whole-blood transcriptome are associated with body mass index: results of an mRNA profiling study involving two large population-based cohorts.

BACKGROUND: Obesity, defined as pathologically increased body mass index (BMI), is strongly related to an increased risk for numerous common cardiovascular and metabolic diseases. It is particularly associated with insulin resistance, hyperglycemia, and systemic oxidative stress and represents the most important risk factor for type 2 diabetes (T2D). However, the pathophysiological mechanisms underlying these associations are still not completely understood. Therefore, in order to identify potentially disease-relevant BMI-associated gene expression signatures, a transcriptome-wide association study (TWAS) on BMI was performed. METHODS: Whole-blood mRNA levels determined by array-based transcriptional profiling were correlated with BMI in two large independent population-based cohort studies (KORA F4 and SHIP-TREND) comprising a total of 1977 individuals. RESULTS: Extensive alterations of the whole-blood transcriptome were associated with BMI: More than 3500 transcripts exhibited significant positive or negative BMI-correlation. Three major whole-blood gene expression signatures associated with increased BMI were identified. The three signatures suggested: i) a ratio shift from mature erythrocytes towards reticulocytes, ii) decreased expression of several genes essentially involved in the transmission and amplification of the insulin signal, and iii) reduced expression of several key genes involved in the defence against reactive oxygen species (ROS). CONCLUSIONS: Whereas the first signature confirms published results, the other two provide possible mechanistic explanations for well-known epidemiological findings under conditions of increased BMI, namely attenuated insulin signaling and increased oxidative stress. The putatively causative BMI-dependent down-regulation of the expression of numerous genes on the mRNA level represents a novel finding. BMI-associated negative transcriptional regulation of insulin signaling and oxidative stress management provide new insights into the pathogenesis of metabolic syndrome and T2D.

Novel insights into the genetics of smoking behaviour, lung function, and chronic obstructive pulmonary disease (UK BiLEVE): a genetic association study in UK Biobank.

BACKGROUND: Understanding the genetic basis of airflow obstruction and smoking behaviour is key to determining the pathophysiology of chronic obstructive pulmonary disease (COPD). We used UK Biobank data to study the genetic causes of smoking behaviour and lung health. METHODS: We sampled individuals of European ancestry from UK Biobank, from the middle and extremes of the forced expiratory volume in 1 s (FEV1) distribution among heavy smokers (mean 35 pack-years) and never smokers. We developed a custom array for UK Biobank to provide optimum genome-wide coverage of common and low-frequency variants, dense coverage of genomic regions already implicated in lung health and disease, and to assay rare coding variants relevant to the UK population. We investigated whether there were shared genetic causes between different phenotypes defined by extremes of FEV1. We also looked for novel variants associated with extremes of FEV1 and smoking behaviour and assessed regions of the genome that had already shown evidence for a role in lung health and disease. We set genome-wide significance at p<5 × 10(-8). FINDINGS: UK Biobank participants were recruited from March 15, 2006, to July 7, 2010. Sample selection for the UK BiLEVE study started on Nov 22, 2012, and was completed on Dec 20, 2012. We selected 50,008 unique samples: 10,002 individuals with low FEV1, 10,000 with average FEV1, and 5002 with high FEV1 from each of the heavy smoker and never smoker groups. We noted a substantial sharing of genetic causes of low FEV1 between heavy smokers and never smokers (p=2.29 × 10(-16)) and between individuals with and without doctor-diagnosed asthma (p=6.06 × 10(-11)). We discovered six novel genome-wide significant signals of association with extremes of FEV1, including signals at four novel loci (KANSL1, TSEN54, TET2, and RBM19/TBX5) and independent signals at two previously reported loci (NPNT and HLA-DQB1/HLA-DQA2). These variants also showed association with COPD, including in individuals with no history of smoking. The number of copies of a 150 kb region containing the 5' end of KANSL1, a gene that is important for epigenetic gene regulation, was associated with extremes of FEV1. We also discovered five new genome-wide significant signals for smoking behaviour, including a variant in NCAM1 (chromosome 11) and a variant on chromosome 2 (between TEX41 and PABPC1P2) that has a trans effect on expression of NCAM1 in brain tissue. INTERPRETATION: By sampling from the extremes of the lung function distribution in UK Biobank, we identified novel genetic causes of lung function and smoking behaviour. These results provide new insight into the specific mechanisms underlying airflow obstruction, COPD, and tobacco addiction, and show substantial shared genetic architecture underlying airflow obstruction across individuals, irrespective of smoking behaviour and other airway disease. FUNDING: Medical Research Council.

The transcriptional landscape of age in human peripheral blood.

Disease incidences increase with age, but the molecular characteristics of ageing that lead to increased disease susceptibility remain inadequately understood. Here we perform a whole-blood gene expression meta-analysis in 14,983 individuals of European ancestry (including replication) and identify 1,497 genes that are differentially expressed with chronological age. The age-associated genes do not harbor more age-associated CpG-methylation sites than other genes, but are instead enriched for the presence of potentially functional CpG-methylation sites in enhancer and insulator regions that associate with both chronological age and gene expression levels. We further used the gene expression profiles to calculate the 'transcriptomic age' of an individual, and show that differences between transcriptomic age and chronological age are associated with biological features linked to ageing, such as blood pressure, cholesterol levels, fasting glucose, and body mass index. The transcriptomic prediction model adds biological relevance and complements existing epigenetic prediction models, and can be used by others to calculate transcriptomic age in external cohorts.

Gene expression alterations in the cerebellum and granule neurons of Cstb(-/-) mouse are associated with early synaptic changes and inflammation.

Progressive myoclonus epilepsy of Unverricht-Lundborg type (EPM1) is an autosomal recessively inherited neurodegenerative disease, manifesting with myoclonus, seizures and ataxia, caused by mutations in the cystatin B (CSTB) gene. With the aim of understanding the molecular basis of pathogenetic events in EPM1 we characterized gene expression changes in the cerebella of pre-symptomatic postnatal day 7 (P7) and symptomatic P30 cystatin B -deficient (Cstb(-/-) ) mice, a model for the disease, and in cultured Cstb(-/-) cerebellar granule cells using a pathway-based approach. Differentially expressed genes in P7 cerebella were connected to synaptic function and plasticity, and in cultured cerebellar granule cells, to cell cycle, cytoskeleton, and intracellular transport. In particular, the gene expression data pinpointed alterations in GABAergic pathway. Electrophysiological recordings from Cstb(-/-) cerebellar Purkinje cells revealed a shift of the balance towards decreased inhibition, yet the amount of inhibitory interneurons was not declined in young animals. Instead, we found diminished number of GABAergic terminals and reduced ligand binding to GABAA receptors in Cstb(-/-) cerebellum. These results suggest that alterations in GABAergic signaling could result in reduced inhibition in Cstb(-/-) cerebellum leading to the hyperexcitable phenotype of Cstb(-/-) mice. At P30, the microarray data revealed a marked upregulation of immune and defense response genes, compatible with the previously reported early glial activation that precedes neuronal degeneration. This further implies the role of early-onset neuroinflammation in the pathogenesis of EPM1.

Mapping the genetic architecture of gene regulation in whole blood.

BACKGROUND: We aimed to assess whether whole blood expression quantitative trait loci (eQTLs) with effects in cis and trans are robust and can be used to identify regulatory pathways affecting disease susceptibility. MATERIALS AND METHODS: We performed whole-genome eQTL analyses in 890 participants of the KORA F4 study and in two independent replication samples (SHIP-TREND, N = 976 and EGCUT, N = 842) using linear regression models and Bonferroni correction. RESULTS: In the KORA F4 study, 4,116 cis-eQTLs (defined as SNP-probe pairs where the SNP is located within a 500 kb window around the transcription unit) and 94 trans-eQTLs reached genome-wide significance and overall 91% (92% of cis-, 84% of trans-eQTLs) were confirmed in at least one of the two replication studies. Different study designs including distinct laboratory reagents (PAXgene™ vs. Tempus™ tubes) did not affect reproducibility (separate overall replication overlap: 78% and 82%). Immune response pathways were enriched in cis- and trans-eQTLs and significant cis-eQTLs were partly coexistent in other tissues (cross-tissue similarity 40-70%). Furthermore, four chromosomal regions displayed simultaneous impact on multiple gene expression levels in trans, and 746 eQTL-SNPs have been previously reported to have clinical relevance. We demonstrated cross-associations between eQTL-SNPs, gene expression levels in trans, and clinical phenotypes as well as a link between eQTLs and human metabolic traits via modification of gene regulation in cis. CONCLUSIONS: Our data suggest that whole blood is a robust tissue for eQTL analysis and may be used both for biomarker studies and to enhance our understanding of molecular mechanisms underlying gene-disease associations.

Systematic identification of trans eQTLs as putative drivers of known disease associations.

Identifying the downstream effects of disease-associated SNPs is challenging. To help overcome this problem, we performed expression quantitative trait locus (eQTL) meta-analysis in non-transformed peripheral blood samples from 5,311 individuals with replication in 2,775 individuals. We identified and replicated trans eQTLs for 233 SNPs (reflecting 103 independent loci) that were previously associated with complex traits at genome-wide significance. Some of these SNPs affect multiple genes in trans that are known to be altered in individuals with disease: rs4917014, previously associated with systemic lupus erythematosus (SLE), altered gene expression of C1QB and five type I interferon response genes, both hallmarks of SLE. DeepSAGE RNA sequencing showed that rs4917014 strongly alters the 3' UTR levels of IKZF1 in cis, and chromatin immunoprecipitation and sequencing analysis of the trans-regulated genes implicated IKZF1 as the causal gene. Variants associated with cholesterol metabolism and type 1 diabetes showed similar phenomena, indicating that large-scale eQTL mapping provides insight into the downstream effects of many trait-associated variants.