Genetic variants associated mRNA stability in lung

Expression quantitative trait loci (eQTLs) analyses have been widely used to identify genetic variants associated with gene expression levels to understand what molecular mechanisms underlie genetic traits. The resultant eQTLs might affect the expression of associated genes through transcriptional or post-transcriptional regulation. In this study, we attempt to distinguish these two types of regulation by identifying genetic variants associated with mRNA stability of genes (stQTLs). Specifically, we computationally inferred mRNA stability of genes based on RNA-seq data and performed association analysis to identify stQTLs. Using the Genotype-Tissue Expression (GTEx) lung RNA-Seq data, we identified a total of 142,801 stQTLs for 3,942 genes and 186,132 eQTLs for 4,751 genes from 15,122,700 genetic variants for 13,476 genes, respectively. Interesting, our results indicated that stQTLs were enriched in the CDS and 3UTR regions, while eQTLs are enriched in the CDS, 3UTR, 5UTR, and upstream regions. We also found that stQTLs are more likely than eQTLs to overlap with RNA binding protein (RBP) and microRNA (miRNA) binding sites. Our analyses demonstrate that simultaneous identification of stQTLs and eQTLs can provide more mechanistic insight on the association between genetic variants and gene expression levels. Author SummaryIn the past decade, many studies have identified genetic variants associated with gene expression level (eQTLs) in different phenotypes, including tissues and diseases. Gene expression is the result of cooperation between transcriptional regulation, such as transcriptional activity, and post-transcriptional regulation, such as mRNA stability. Here, we present a computational framework that take advantage of recently developed methods to estimate mRNA stability from RNA-Seq, which is widely used to estimate gene expression, and then to identify genetic variants associated with mRNA stability (stQTLs) in lung tissue. Compared to eQTLs, we found that genetic variants that affects mRNA stability are more significantly located in the CDS and 3UTR regions, which are known to interact with RNA-binding proteins (RBPs) or microRNAs to regulate stability. In addition, stQTLs are significantly more likely to overlap the binding sites of RBPs. We show that the six RBPs that most significantly bind to stQTLs are all known to regulate mRNA stability. This pipeline of simultaneously identifying eQTLs and stQTLs using only RNA-Seq data can provide higher resolution than traditional eQTLs study to better understand the molecular mechanisms of genetic variants on the regulation of gene expression.