Abstract
Background Genome-wide association studies (GWAS) have identified hundreds of loci underlying adult-onset asthma (AOA) and childhood-onset asthma (COA). However, the causal variants, regulatory elements, and effector genes at these loci are largely unknown.
Methods We performed heritability enrichment analysis to determine relevant cell types for AOA and COA, respectively. Next, we fine-mapped putative causal variants at AOA and COA loci. To improve the resolution of fine-mapping, we integrated ATAC-seq data in blood and lung cell types to annotate variants in candidate cis-regulatory elements (CREs).
We then computationally prioritized candidate CREs underlying asthma risk, experimentally assessed their enhancer activity by massively parallel reporter assay (MPRA) in bronchial epithelial cells (BECs) and further validated a subset by luciferase assays. Combining chromatin interaction data and expression quantitative trait loci, we nominated genes targeted by candidate CREs and prioritized effector genes for AOA and COA.
Results Heritability enrichment analysis suggested a shared role of immune cells in the development of both AOA and COA while highlighting the distinct contribution of lung structural cells in COA. Functional fine-mapping uncovered 21 and 67 credible sets for AOA and COA, respectively, with only 16% shared between the two. Notably, one-third of the loci contained multiple credible sets. Our CRE prioritization strategy nominated 62 and 169 candidate CREs for AOA and COA, respectively. Over 60% of these candidate CREs showed open chromatin in multiple cell lineages, suggesting their potential pleiotropic effects in different cell types. Furthermore, COA candidate CREs were enriched for enhancers experimentally validated by MPRA in BECs. The prioritized effector genes included many genes involved in immune and inflammatory responses. Notably, multiple genes, including TNFSF4, a drug target undergoing clinical trials, were supported by two independent GWAS signals, indicating widespread allelic heterogeneity. Four out of six selected candidate CREs demonstrated allele-specific regulatory properties in luciferase assays in BECs.
Conclusions We present a comprehensive characterization of causal variants, regulatory elements, and effector genes underlying AOA and COA genetics. Our results supported a distinct genetic basis between AOA and COA and highlighted regulatory complexity at many GWAS loci marked by both extensive pleiotropy and allelic heterogeneity.
Competing Interest Statement
The authors have declared no competing interest.
Funding Statement
This study was supported by U19 AI62310 (C.O., M.A.N., A.I.S.), R01 MH110531, R01 MH116281, R01 HG010773, R01 HL163523 (X.H.), and UG3/UH3 OD023282 and UM1 AI160040 (C.O.). I.M.S. was supported by T32 HL007605. N.S. was supported by K08 HL153955.
Author Declarations
I confirm all relevant ethical guidelines have been followed, and any necessary IRB and/or ethics committee approvals have been obtained.
Yes
The details of the IRB/oversight body that provided approval or exemption for the research described are given below:
The Institutional Review Board (IRB) at the University of Chicago waived ethical approval for this study.
I confirm that all necessary patient/participant consent has been obtained and the appropriate institutional forms have been archived, and that any patient/participant/sample identifiers included were not known to anyone (e.g., hospital staff, patients or participants themselves) outside the research group so cannot be used to identify individuals.
Yes
I understand that all clinical trials and any other prospective interventional studies must be registered with an ICMJE-approved registry, such as ClinicalTrials.gov. I confirm that any such study reported in the manuscript has been registered and the trial registration ID is provided (note: if posting a prospective study registered retrospectively, please provide a statement in the trial ID field explaining why the study was not registered in advance).
Yes
I have followed all appropriate research reporting guidelines, such as any relevant EQUATOR Network research reporting checklist(s) and other pertinent material, if applicable.
Yes
Abbreviations
- GWAS
- genome-wide association study
- AOA
- adult-onset asthma
- COA
- childhood-onset asthma
- CRE
- cis-regulatory element
- MPRA
- massively parallel reporter assay
- BEC
- bronchial epithelial cell
- eQTL
- expression quantitative trait locus
- LD
- linkage disequilibrium
- TWAS
- transcriptome-wide association study
- sQTL
- splicing quantitative trait locus
- caQTL
- chromatin accessibility quantitative trait locus
- meQTL
- DNA methylation quantitative trait locus
- UKB
- UK Biobank
- SNP
- single nucleotide polymorphism
- S-LDSC
- stratified LD score regression
- OCR
- open chromatin region
- snATAC-seq
- single-nucleus ATAC-seq
- SuSiE
- sum of single effects
- PIP
- posterior inclusion probability
- ePIP
- element PIP
- PCHi-C
- promoter capture
- Hi-C ABC
- activity-by-contact
- GO
- Gene Ontology
- CPM
- counts per million
- TSS
- transcription start site
- TLR
- toll-like receptors
- OX40L
- OX40 ligand
- RV
- rhinovirus
The Chan Zuckerberg Initiative, Cold Spring Harbor Laboratory, the Sergey Brin Family Foundation, California Institute of Technology, Centre National de la Recherche Scientifique, Fred Hutchinson Cancer Center, Imperial College London, Massachusetts Institute of Technology, Stanford University, University of Washington, and Vrije Universiteit Amsterdam.