Mapping eQTLs with RNA-seq reveals novel susceptibility genes, non-coding RNAs and alternative-splicing events in systemic lupus erythematosus
Open Access
- 28 February 2017
- journal article
- research article
- Published by Oxford University Press (OUP) in Human Molecular Genetics
- Vol. 26 (5) , 1003-1017
Abstract
Studies attempting to functionally interpret complex-disease susceptibility loci by GWAS and eQTL integration have predominantly employed microarrays to quantify gene-expression. RNA-Seq has the potential to discover a more comprehensive set of eQTLs and illuminate the underlying molecular consequence. We examine the functional outcome of 39 variants associated with Systemic Lupus Erythematosus (SLE) through the integration of GWAS and eQTL data from the TwinsUK microarray and RNA-Seq cohort in lymphoblastoid cell lines. We use conditional analysis and a Bayesian colocalisation method to provide evidence of a shared causal-variant, then compare the ability of each quantification type to detect disease relevant eQTLs and eGenes. We discovered the greatest frequency of candidate-causal eQTLs using exon-level RNA-Seq, and identified novel SLE susceptibility genes (e.g. NADSYN1 and TCF7) that were concealed using microarrays, including four non-coding RNAs. Many of these eQTLs were found to influence the expression of several genes, supporting the notion that risk haplotypes may harbour multiple functional effects. Novel SLE associated splicing events were identified in the T-reg restricted transcription factor, IKZF2, and other candidate genes (e.g. WDFY4) through asQTL mapping using the Geuvadis cohort. We have significantly increased our understanding of the genetic control of gene-expression in SLE by maximising the leverage of RNA-Seq and performing integrative GWAS-eQTL analysis against gene, exon, and splice-junction quantifications. We conclude that to better understand the true functional consequence of regulatory variants, quantification by RNA-Seq should be performed at the exon-level as a minimum, and run in parallel with gene and splice-junction level quantification.Keywords
All Related Versions
This publication has 86 references indexed in Scilit:
- An integrated encyclopedia of DNA elements in the human genomeNature, 2012
- eQTL Mapping Using RNA-seq DataStatistics in Biosciences, 2012
- GCTA: A Tool for Genome-wide Complex Trait AnalysisThe American Journal of Human Genetics, 2011
- Common genetic determinants of vitamin D insufficiency: a genome-wide association studyThe Lancet, 2010
- A role for the transcription factor Helios in human CD4+CD25+ regulatory T cellsMolecular Immunology, 2010
- Transcriptome genetics using second generation sequencing in a Caucasian populationNature, 2010
- Understanding mechanisms underlying human gene expression variation with RNA sequencingNature, 2010
- Genetic variation at the IRF7/PHRF1 locus is associated with autoantibody profile and serum interferon‐α activity in lupus patientsArthritis & Rheumatism, 2010
- Genetic variants regulating ORMDL3 expression contribute to the risk of childhood asthmaNature, 2007
- Common genetic variants account for differences in gene expression among ethnic groupsNature Genetics, 2007