Mapping eQTL by leveraging multiple tissues and DNA methylation

BMC Bioinformatics

Methodology Article

Kouros Owzar¹  Andrew S. Allen²  Chaitanya R. Acharya² 

[1] Department of Biostatistics and Bioinformatics, Duke University, 2424 Erwin Road, Suite 1104, 27710, Durham, NC, USA;Program in Computational Biology and Bioinformatics, Duke University, 2424 Erwin Road, Suite 1104, 27710, Durham, NC, USA;Department of Biostatistics and Bioinformatics, Duke University, 2424 Erwin Road, Suite 1104, 27710, Durham, NC, USA;

关键词: eQTL;    Multiple tissues;    Tissue-specificity;    DNA methylation;    CpG islands;    Gene expression;    SNP;    Score test;    Monte Carlo simulations;    Brain;   

DOI  :  10.1186/s12859-017-1856-9

 received in 2016-10-05, accepted in 2017-10-02,  发布年份 2017

来源: Springer

【 摘 要 】

BackgroundDNA methylation is an important tissue-specific epigenetic event that influences transcriptional regulation of gene expression. Differentially methylated CpG sites may act as mediators between genetic variation and gene expression, and this relationship can be exploited while mapping multi-tissue expression quantitative trait loci (eQTL). Current multi-tissue eQTL mapping techniques are limited to only exploiting gene expression patterns across multiple tissues either in a joint tissue or tissue-by-tissue frameworks. We present a new statistical approach that enables us to model the effect of germ-line variation on tissue-specific gene expression in the presence of effects due to DNA methylation.ResultsOur method efficiently models genetic and epigenetic variation to identify genomic regions of interest containing combinations of mRNA transcripts, CpG sites, and SNPs by jointly testing for genotypic effect and higher order interaction effects between genotype, methylation and tissues. We demonstrate using Monte Carlo simulations that our approach, in the presence of both genetic and DNA methylation effects, gives an improved performance (in terms of statistical power) to detect eQTLs over the current eQTL mapping approaches. When applied to an array-based dataset from 150 neuropathologically normal adult human brains, our method identifies eQTLs that were undetected using standard tissue-by-tissue or joint tissue eQTL mapping techniques. As an example, our method identifies eQTLs by leveraging methylated CpG sites in a LIM homeobox member gene (LHX9), which may have a role in the neural development.ConclusionsOur score test-based approach does not need parameter estimation under the alternative hypothesis. As a result, our model parameters are estimated only once for each mRNA - CpG pair. Our model specifically studies the effects of non-coding regions of DNA (in this case, CpG sites) on mapping eQTLs. However, we can easily model micro-RNAs instead of CpG sites to study the effects of post-transcriptional events in mapping eQTL. Our model’s flexible framework also allows us to investigate other genomic events such as alternative gene splicing by extending our model to include gene isoform-specific data.

【 授权许可】

CC BY   
© The Author(s) 2017

【 图 表 】

12864_2017_4004_Article_IEq4.gif

12864_2017_3670_Article_IEq8.gif

12864_2017_4186_Article_IEq18.gif

12864_2017_4186_Article_IEq16.gif

12864_2017_3670_Article_IEq7.gif

12864_2016_3477_Article_IEq9.gif

12864_2017_3790_Article_IEq1.gif

12864_2017_4186_Article_IEq13.gif

12864_2017_4186_Article_IEq12.gif

12864_2017_4186_Article_IEq11.gif

12864_2015_2300_Article_IEq4.gif

12864_2017_3938_Article_IEq8.gif

12864_2016_3317_Article_IEq1.gif

12864_2017_3938_Article_IEq6.gif

12864_2015_2055_Article_IEq84.gif

12864_2017_3487_Article_IEq6.gif

12864_2017_3500_Article_IEq22.gif

12864_2017_3496_Article_IEq2.gif

12864_2017_3487_Article_IEq4.gif

12864_2017_3938_Article_IEq1.gif

12864_2017_3500_Article_IEq20.gif

12864_2017_3500_Article_IEq19.gif

12864_2016_3440_Article_IEq34.gif

12864_2017_4258_Article_IEq3.gif

12864_2016_2682_Article_IEq45.gif