The Impact of Cell Type and Context-dependent Regulatory Variants on Human Immune Traits

Overview

Journal Genome Biol

Specialties Biology
Genetics

Date 2021 Apr 30

PMID 33926512

Citations 36

Authors

Zepeng Mu

Wei Wei

Benjamin Fair

Jinlin Miao

Ping Zhu

Yang I Li

Affiliations

Soon will be listed here.

Abstract

Background: The vast majority of trait-associated variants identified using genome-wide association studies (GWAS) are noncoding, and therefore assumed to impact gene regulation. However, the majority of trait-associated loci are unexplained by regulatory quantitative trait loci (QTLs).

Results: We perform a comprehensive characterization of the putative mechanisms by which GWAS loci impact human immune traits. By harmonizing four major immune QTL studies, we identify 26,271 expression QTLs (eQTLs) and 23,121 splicing QTLs (sQTLs) spanning 18 immune cell types. Our colocalization analyses between QTLs and trait-associated loci from 72 GWAS reveals that genetic effects on RNA expression and splicing in immune cells colocalize with 40.4% of GWAS loci for immune-related traits, in many cases increasing the fraction of colocalized loci by two fold compared to previous studies. Notably, we find that the largest contributors of this increase are splicing QTLs, which colocalize on average with 14% of all GWAS loci that do not colocalize with eQTLs. By contrast, we find that cell type-specific eQTLs, and eQTLs with small effect sizes contribute very few new colocalizations. To investigate the 60% of GWAS loci that remain unexplained, we collect H3K27ac CUT&Tag data from rheumatoid arthritis and healthy controls, and find large-scale differences between immune cells from the different disease contexts, including at regions overlapping unexplained GWAS loci.

Conclusion: Altogether, our work supports RNA splicing as an important mediator of genetic effects on immune traits, and suggests that we must expand our study of regulatory processes in disease contexts to improve functional interpretation of as yet unexplained GWAS loci.

Citing Articles

Integration of functional genomics and statistical fine-mapping systematically characterizes adult-onset and childhood-onset asthma genetic associations.

Zhong X, Mitchell R, Billstrand C, Thompson E, Sakabe N, Aneas I medRxiv. 2025; .

PMID: 40034789 PMC: 11875274. DOI: 10.1101/2025.02.11.25322088.

Long-read RNA sequencing atlas of human microglia isoforms elucidates disease-associated genetic regulation of splicing.

Humphrey J, Brophy E, Kosoy R, Zeng B, Coccia E, Mattei D Nat Genet. 2025; 57(3):604-615.

PMID: 40033057 DOI: 10.1038/s41588-025-02099-0.

Implications of gene × environment interactions in post-traumatic stress disorder risk and treatment.

Seah C, Sidamon-Eristoff A, Huckins L, Brennand K J Clin Invest. 2025; 135(5).

PMID: 40026250 PMC: 11870735. DOI: 10.1172/JCI185102.

Widespread gene-environment interactions shape the immune response to SARS-CoV-2 infection in hospitalized COVID-19 patients.

Randolph H, Aguirre-Gamboa R, Brunet-Ratnasingham E, Nakanishi T, Locher V, Ketter E bioRxiv. 2024; .

PMID: 39677792 PMC: 11642875. DOI: 10.1101/2024.12.03.626676.

Pleiotropic expression quantitative trait loci are enriched in enhancers and transcription factor binding sites and impact more genes.

Gonzalez A, Paul P Comput Struct Biotechnol J. 2024; 23:4260-4270.

PMID: 39669750 PMC: 11635986. DOI: 10.1016/j.csbj.2024.11.019.

References

Liu J, van Sommeren S, Huang H, Ng S, Alberts R, Takahashi A . Association analyses identify 38 susceptibility loci for inflammatory bowel disease and highlight shared genetic risk across populations. Nat Genet. 2015; 47(9):979-986. PMC: 4881818. DOI: 10.1038/ng.3359. View

Katz Y, Wang E, Airoldi E, Burge C . Analysis and design of RNA sequencing experiments for identifying isoform regulation. Nat Methods. 2010; 7(12):1009-15. PMC: 3037023. DOI: 10.1038/nmeth.1528. View

Watanabe K, Stringer S, Frei O, Mirkov M, de Leeuw C, Polderman T . A global overview of pleiotropy and genetic architecture in complex traits. Nat Genet. 2019; 51(9):1339-1348. DOI: 10.1038/s41588-019-0481-0. View

Foley C, Staley J, Breen P, Sun B, Kirk P, Burgess S . A fast and efficient colocalization algorithm for identifying shared genetic risk factors across multiple traits. Nat Commun. 2021; 12(1):764. PMC: 7858636. DOI: 10.1038/s41467-020-20885-8. View

Liu B, Gloudemans M, Rao A, Ingelsson E, Montgomery S . Abundant associations with gene expression complicate GWAS follow-up. Nat Genet. 2019; 51(5):768-769. PMC: 6904208. DOI: 10.1038/s41588-019-0404-0. View

Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N . The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009; 25(16):2078-9. PMC: 2723002. DOI: 10.1093/bioinformatics/btp352. View

Karaky M, Fedetz M, Potenciano V, Andres-Leon E, Codina A, Barrionuevo C . SP140 regulates the expression of immune-related genes associated with multiple sclerosis and other autoimmune diseases by NF-κB inhibition. Hum Mol Genet. 2018; 27(23):4012-4023. DOI: 10.1093/hmg/ddy284. View

Delaneau O, Ongen H, Brown A, Fort A, Panousis N, Dermitzakis E . A complete tool set for molecular QTL discovery and analysis. Nat Commun. 2017; 8:15452. PMC: 5454369. DOI: 10.1038/ncomms15452. View

Chen L, Ge B, Casale F, Vasquez L, Kwan T, Garrido-Martin D . Genetic Drivers of Epigenetic and Transcriptional Variation in Human Immune Cells. Cell. 2016; 167(5):1398-1414.e24. PMC: 5119954. DOI: 10.1016/j.cell.2016.10.026. View

10.

Wang X, Goldstein D . Enhancer Domains Predict Gene Pathogenicity and Inform Gene Discovery in Complex Disease. Am J Hum Genet. 2020; 106(2):215-233. PMC: 7010980. DOI: 10.1016/j.ajhg.2020.01.012. View

11.

Schmiedel B, Singh D, Madrigal A, Valdovino-Gonzalez A, White B, Zapardiel-Gonzalo J . Impact of Genetic Polymorphisms on Human Immune Cell Gene Expression. Cell. 2018; 175(6):1701-1715.e16. PMC: 6289654. DOI: 10.1016/j.cell.2018.10.022. View

12.

Visscher P, Wray N, Zhang Q, Sklar P, McCarthy M, Brown M . 10 Years of GWAS Discovery: Biology, Function, and Translation. Am J Hum Genet. 2017; 101(1):5-22. PMC: 5501872. DOI: 10.1016/j.ajhg.2017.06.005. View

13.

Calderon D, Nguyen M, Mezger A, Kathiria A, Muller F, Nguyen V . Landscape of stimulation-responsive chromatin across diverse human immune cells. Nat Genet. 2019; 51(10):1494-1505. PMC: 6858557. DOI: 10.1038/s41588-019-0505-9. View

14.

Hu X, Kim H, Raj T, Brennan P, Trynka G, Teslovich N . Regulation of gene expression in autoimmune disease loci and the genetic basis of proliferation in CD4+ effector memory T cells. PLoS Genet. 2014; 10(6):e1004404. PMC: 4072514. DOI: 10.1371/journal.pgen.1004404. View

15.

Courtney P, Doherty M . Joint aspiration and injection and synovial fluid analysis. Best Pract Res Clin Rheumatol. 2013; 27(2):137-69. DOI: 10.1016/j.berh.2013.02.005. View

16.

Paternoster L, Standl M, Waage J, Baurecht H, Hotze M, Strachan D . Multi-ancestry genome-wide association study of 21,000 cases and 95,000 controls identifies new risk loci for atopic dermatitis. Nat Genet. 2015; 47(12):1449-1456. PMC: 4753676. DOI: 10.1038/ng.3424. View

17.

van de Geijn B, McVicker G, Gilad Y, Pritchard J . WASP: allele-specific software for robust molecular quantitative trait locus discovery. Nat Methods. 2015; 12(11):1061-3. PMC: 4626402. DOI: 10.1038/nmeth.3582. View

18.

Liao Y, Smyth G, Shi W . featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2013; 30(7):923-30. DOI: 10.1093/bioinformatics/btt656. View

19.

Kim-Hellmuth S, Bechheim M, Putz B, Mohammadi P, Nedelec Y, Giangreco N . Genetic regulatory effects modified by immune activation contribute to autoimmune disease associations. Nat Commun. 2017; 8(1):266. PMC: 5559603. DOI: 10.1038/s41467-017-00366-1. View

20.

Li Y, van de Geijn B, Raj A, Knowles D, Petti A, Golan D . RNA splicing is a primary link between genetic variation and disease. Science. 2016; 352(6285):600-4. PMC: 5182069. DOI: 10.1126/science.aad9417. View