Abstract
Extensive studies have been conducted on the analysis of genome function, especially on the expression quantitative trait loci (eQTL). These studies offered promising results for characterization of the functional sequencing variation and understanding of the basic processes of gene regulation. Parent of origin effect (POE) is an important epigenetic phenomenon describing that the expression of certain genes depends on their allelic parent-of-origin and it is known to play important roles in human complex diseases. However, traditional eQTL mapping approaches do not allow for the detection of imprinting, or they focus on modeling the additive genetic effect thereby ignoring the estimation of the dominance genetic effect. In this study, we proposed a statistical framework to test the additive and dominance genetic effects of the candidate eQTLs along with detection of the POE with a functional model and an orthogonal model for RNA-seq data. We demonstrated the desirable power and preserved Type I errors of the methods in most scenarios, especially the orthogonal model with un-biased estimation of the genetic effects and over-dispersion of the RNA-seq data. The application to a HapMap project trio dataset validated existing imprinting genes and discovered two novel imprinting genes with potential dominance genetic effect and RB1 and IGF1R genes. This study provides new insights into the next generation statistical modeling of eQTL mapping for better understanding of the genetic architecture underlying the mechanisms of gene expression regulation.
Similar content being viewed by others
References
Álvarez-Castro JM (2014) Dissecting genetic effects with imprinting. Front Ecol Evol. https://doi.org/10.3389/fevo.2014.00051
Alvarez-Castro JM, Carlborg O (2007) A unified model for functional and statistical epistasis and its application in quantitative trait Loci analysis. Genetics 176:1151–1167. https://doi.org/10.1534/genetics.106.067348
Benjamini Y, Hochberg Y (1995) controlling the false discovery rate—a practical and powerful approach to multiple testing. J R Stat Soc B 57:289–300
Browning BL, Browning SR (2009) A unified approach to genotype imputation and haplotype-phase inference for large data sets of trios and unrelated individuals. Am J Hum Genet 84:210–223. https://doi.org/10.1016/j.ajhg.2009.01.005
Chinnam M, Goodrich DW (2011) Rb1, development, and cancer. Curr Top Dev Biol 94:129–169. https://doi.org/10.1016/B978-0-12-380916-2.00005-X
Cookson W, Liang L, Abecasis G, Moffatt M, Lathrop M (2009) Mapping complex disease traits with global gene expression. Nat Rev Genet 10:184–194. https://doi.org/10.1038/nrg2537
Delaneau O, Marchini J, Genomes Project C, Genomes Project C (2014) Integrating sequence and array data to create an improved 1000 Genomes Project haplotype reference panel. Nat Commun 5:3934. https://doi.org/10.1038/ncomms4934
Duenk P, Calus MPL, Wientjes YCJ, Bijma P (2017) Benefits of dominance over additive models for the estimation of average effects in the presence of dominance. G3 (Bethesda) 7:3405–3414. https://doi.org/10.1534/g3.117.300113
Ellis SE, Gupta S, Ashar FN, Bader JS, West AB, Arking DE (2013) RNA-Seq optimization with eQTL gold standards. BMC Genom 14:892. https://doi.org/10.1186/1471-2164-14-892
Ertl J, Legarra A, Vitezica ZG, Varona L, Edel C, Emmerling R, Gotz KU (2014) Genomic analysis of dominance effects on milk production and conformation traits in Fleckvieh cattle. Genet Sel Evol 46:40. https://doi.org/10.1186/1297-9686-46-40
Genomes Project C et al (2012) An integrated map of genetic variation from 1092 human genomes. Nature 491:56–65. https://doi.org/10.1038/nature11632
Genomes Project C et al (2015) A global reference for human genetic variation. Nature 526:68–74. https://doi.org/10.1038/nature15393
Green PJ (1984) Iteratively reweighted least squares for maximum likelihood estimation, and some robust and resistant alternatives. J R Stat Soc Ser B (Methodol) 46:149–192
Howie B, Fuchsberger C, Stephens M, Marchini J, Abecasis GR (2012) Fast and accurate genotype imputation in genome-wide association studies through pre-phasing. Nat Genet 44:955. https://doi.org/10.1038/ng.2354
Huang W, Mackay TF (2016) The genetic architecture of quantitative traits cannot be inferred from variance component analysis. PLoS Genet 12:e1006421. https://doi.org/10.1371/journal.pgen.1006421
Jadhav B et al (2019) RNA-Seq in 296 phased trios provides a high-resolution map of genomic imprinting. BMC Biol 17:50. https://doi.org/10.1186/s12915-019-0674-0
Kanber D et al (2009) The Human Retinoblastoma Gene Is Imprinted. Plos Genet 5:e1000790. https://doi.org/10.1371/journal.pgen.1000790
Kong A et al (2009) Parental origin of sequence variants associated with complex diseases. Nature 462:868–U859. https://doi.org/10.1038/nature08625
Ma J et al (2012) Natural and orthogonal interaction framework for modeling gene-environment interactions with application to lung cancer. Hum Hered 73:185–194. https://doi.org/10.1159/000339906
Palowitch J, Shabalin A, Zhou YH, Nobel AB, Wright FA (2018) Estimation of cis-eQTL effect sizes using a log of linear model. Biometrics 74:616–625. https://doi.org/10.1111/biom.12810
Perry JRB et al (2014) Parent-of-origin-specific allelic associations among 106 genomic loci for age at menarche. Nature 514:92. https://doi.org/10.1038/nature13545
Rao CR (2005) Score test: historical review and recent developments. In: Balakrishnan N, Nagaraja HN, Kannan N (eds) Advances in ranking and selection, multiple comparisons, and reliability. Statistics for Industry and Technology, Birkhäuser, Boston, pp 3–20. https://doi.org/10.1007/0-8176-4422-9_1
Su G, Christensen OF, Ostersen T, Henryon M, Lund MS (2012) Estimating additive and non-additive genetic variances and predicting genetic merits using genome-wide dense single nucleotide polymorphism markers. PLoS ONE 7:e45293. https://doi.org/10.1371/journal.pone.0045293
Wallace C, Smyth DJ, Maisuria-Armer M, Walker NM, Todd JA, Clayton DG (2010) The imprinted DLK1-MEG3 gene region on chromosome 14q32.2 alters susceptibility to type 1 diabetes. Nat Genet 42:68–U85. https://doi.org/10.1038/ng.493
Wolf JB, Cheverud JM (2009) A framework for detecting and characterizing genetic background-dependent imprinting effects. Mamm Genome 20:681–698. https://doi.org/10.1007/s00335-009-9209-2
Xiao F, Ma J, Amos CI (2013) A unified framework integrating parent-of-origin effects for association study. PLoS ONE 8:e72208. https://doi.org/10.1371/journal.pone.0072208
Xiao F, Ma J, Cai G, Fang S, Lee JE, Wei Q, Amos CI (2014) Natural and orthogonal model for estimating gene-gene interactions applied to cutaneous melanoma. Hum Genet 133:559–574. https://doi.org/10.1007/s00439-013-1392-2
Xiang T, Christensen OF, Vitezica ZG, Legarra A (2018) Genomic model with correlation between additive and dominance effects. Genetics 209:711–723. https://doi.org/10.1534/genetics.118.301015
Zhabotynsky V, Inoue K, Magnuson T, Mauro Calabrese J, Sun W (2019) A statistical method for joint estimation of cis-eQTLs and parent-of-origin effects under family trio design. Biometrics. https://doi.org/10.1111/biom.13026
Acknowledgements
The preparation of this manuscript was supported by the internal funding from the University of South Carolina for Dr. Feifei Xiao. We are extremely grateful for the generous support of Dr. Wei Sun and Dr. Zhabotynsky, who kindly provided the phased and imputed genotypes and pre-processed RNA-seq data in the application study. We also acknowledge Dr. Guoshuai Cai for his valuable inputs in the RNA-seq data statistical modeling and analysis.
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Conflict of interest
On behalf of all authors, the corresponding author states that there is no conflict of interest.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Below is the link to the electronic supplementary material.
Rights and permissions
About this article
Cite this article
Deng, S., Hardin, J., Amos, C.I. et al. Joint modeling of eQTLs and parent-of-origin effects using an orthogonal framework with RNA-seq data. Hum Genet 139, 1107–1117 (2020). https://doi.org/10.1007/s00439-020-02162-2
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00439-020-02162-2