Robust Identification of Regulatory Variants (eQTLs) Using a Differential Expression Framework Developed for RNA-sequencing

Overview

Journal J Anim Sci Biotechnol

Publisher Biomed Central

Date 2023 May 4

PMID 37143150

Authors

Mackenzie A Marrella

Fernando H Biase

Affiliations

Soon will be listed here.

Abstract

Background: A gap currently exists between genetic variants and the underlying cell and tissue biology of a trait, and expression quantitative trait loci (eQTL) studies provide important information to help close that gap. However, two concerns that arise with eQTL analyses using RNA-sequencing data are normalization of data across samples and the data not following a normal distribution. Multiple pipelines have been suggested to address this. For instance, the most recent analysis of the human and farm Genotype-Tissue Expression (GTEx) project proposes using trimmed means of M-values (TMM) to normalize the data followed by an inverse normal transformation.

Results: In this study, we reasoned that eQTL analysis could be carried out using the same framework used for differential gene expression (DGE), which uses a negative binomial model, a statistical test feasible for count data. Using the GTEx framework, we identified 35 significant eQTLs (P < 5 × 10) following the ANOVA model and 39 significant eQTLs (P < 5 × 10) following the additive model. Using a differential gene expression framework, we identified 930 and six significant eQTLs (P < 5 × 10) following an analytical framework equivalent to the ANOVA and additive model, respectively. When we compared the two approaches, there was no overlap of significant eQTLs between the two frameworks. Because we defined specific contrasts, we identified trans eQTLs that more closely resembled what we expect from genetic variants showing complete dominance between alleles. Yet, these were not identified by the GTEx framework.

Conclusions: Our results show that transforming RNA-sequencing data to fit a normal distribution prior to eQTL analysis is not required when the DGE framework is employed. Our proposed approach detected biologically relevant variants that otherwise would not have been identified due to data transformation to fit a normal distribution.

References

Elston R, Satagopan J, Sun S . Genetic terminology. Methods Mol Biol. 2012; 850:1-9. PMC: 4450815. DOI: 10.1007/978-1-61779-555-8_1. View

Loguercio S, Overall R, Michaelson J, Wiltshire T, Pletcher M, Miller B . Integrative analysis of low- and high-resolution eQTL. PLoS One. 2010; 5(11):e13920. PMC: 2978079. DOI: 10.1371/journal.pone.0013920. View

Price A, Patterson N, Plenge R, Weinblatt M, Shadick N, Reich D . Principal components analysis corrects for stratification in genome-wide association studies. Nat Genet. 2006; 38(8):904-9. DOI: 10.1038/ng1847. View

Dendrou C, Plagnol V, Fung E, Yang J, Downes K, Cooper J . Cell-specific protein phenotypes for the autoimmune locus IL2RA using a genotype-selectable human bioresource. Nat Genet. 2009; 41(9):1011-5. PMC: 2749506. DOI: 10.1038/ng.434. View

Fang L, Cai W, Liu S, Canela-Xandri O, Gao Y, Jiang J . Comprehensive analyses of 723 transcriptomes enhance genetic and biological interpretations for complex traits in cattle. Genome Res. 2020; 30(5):790-801. PMC: 7263193. DOI: 10.1101/gr.250704.119. View

Wood A, Esko T, Yang J, Vedantam S, Pers T, Gustafsson S . Defining the role of common variation in the genomic and biological architecture of adult human height. Nat Genet. 2014; 46(11):1173-86. PMC: 4250049. DOI: 10.1038/ng.3097. View

Lam S, Zeidan J, Miglior F, Suarez-Vega A, Gomez-Redondo I, Fonseca P . Development and comparison of RNA-sequencing pipelines for more accurate SNP identification: practical example of functional SNP detection associated with feed efficiency in Nellore beef cattle. BMC Genomics. 2020; 21(1):703. PMC: 7545862. DOI: 10.1186/s12864-020-07107-7. View

Goddard M, Kemper K, MacLeod I, Chamberlain A, Hayes B . Genetics of complex traits: prediction of phenotype, identification of causal polymorphisms and genetic architecture. Proc Biol Sci. 2016; 283(1835). PMC: 4971198. DOI: 10.1098/rspb.2016.0569. View

Gjuvsland A, Plahte E, Adnoy T, Omholt S . Allele interaction--single locus genetics meets regulatory biology. PLoS One. 2010; 5(2):e9379. PMC: 2826424. DOI: 10.1371/journal.pone.0009379. View

10.

Albert F, Bloom J, Siegel J, Day L, Kruglyak L . Genetics of -regulatory variation in gene expression. Elife. 2018; 7. PMC: 6072440. DOI: 10.7554/eLife.35471. View

11.

Kreitmaier P, Katsoula G, Zeggini E . Insights from multi-omics integration in complex disease primary tissues. Trends Genet. 2022; 39(1):46-58. DOI: 10.1016/j.tig.2022.08.005. View

12.

McCarthy D, Chen Y, Smyth G . Differential expression analysis of multifactor RNA-Seq experiments with respect to biological variation. Nucleic Acids Res. 2012; 40(10):4288-97. PMC: 3378882. DOI: 10.1093/nar/gks042. View

13.

Yang J, Wang D, Yang Y, Yang W, Jin W, Niu X . A systematic comparison of normalization methods for eQTL analysis. Brief Bioinform. 2021; 22(6). DOI: 10.1093/bib/bbab193. View

14.

Wagner G, Kin K, Lynch V . Measurement of mRNA abundance using RNA-seq data: RPKM measure is inconsistent among samples. Theory Biosci. 2012; 131(4):281-5. DOI: 10.1007/s12064-012-0162-3. View

15.

Hunt S, Moore B, Amode R, Armean I, Lemos D, Mushtaq A . Annotating and prioritizing genomic variants using the Ensembl Variant Effect Predictor-A tutorial. Hum Mutat. 2021; 43(8):986-997. PMC: 7613081. DOI: 10.1002/humu.24298. View

16.

McLaren W, Gil L, Hunt S, Riat H, Ritchie G, Thormann A . The Ensembl Variant Effect Predictor. Genome Biol. 2016; 17(1):122. PMC: 4893825. DOI: 10.1186/s13059-016-0974-4. View

17.

Peer I, Yelensky R, Altshuler D, Daly M . Estimation of the multiple testing burden for genomewide association studies of nearly all common variants. Genet Epidemiol. 2008; 32(4):381-5. DOI: 10.1002/gepi.20303. View

18.

Strunz T, Grassmann F, Gayan J, Nahkuri S, Souza-Costa D, Maugeais C . A mega-analysis of expression quantitative trait loci (eQTL) provides insight into the regulatory architecture of gene expression variation in liver. Sci Rep. 2018; 8(1):5865. PMC: 5897392. DOI: 10.1038/s41598-018-24219-z. View

19.

Horita N, Kaneko T . Genetic model selection for a case-control study and a meta-analysis. Meta Gene. 2015; 5:1-8. PMC: 4443430. DOI: 10.1016/j.mgene.2015.04.003. View

20.

Shabalin A . Matrix eQTL: ultra fast eQTL analysis via large matrix operations. Bioinformatics. 2012; 28(10):1353-8. PMC: 3348564. DOI: 10.1093/bioinformatics/bts163. View