Epigenetic responses to rhinovirus exposure in airway epithelial cells are correlated with key transcriptional pathways in chronic rhinosinusitis.

BACKGROUND: Viruses may drive immune mechanisms responsible for chronic rhinosinusitis with nasal polyposis (CRSwNP), but little is known about the underlying molecular mechanisms. OBJECTIVES: To identify epigenetic and transcriptional responses to a common upper respiratory pathogen, rhinovirus (RV), that are specific to patients with CRSwNP using a primary sinonasal epithelial cell culture model. METHODS: Airway epithelial cells were collected at surgery from patients with CRSwNP (cases) and from controls without sinus disease, cultured, and then exposed to RV or vehicle for 48 h. Differential gene expression and DNA methylation (DNAm) between cases and controls in response to RV were determined using linear mixed models. Weighted gene co-expression analysis (WGCNA) was used to identify (a) co-regulated gene expression and DNAm signatures, and (b) genes, pathways, and regulatory mechanisms specific to CRSwNP. RESULTS: We identified 5585 differential transcriptional and 261 DNAm responses (FDR <0.10) to RV between CRSwNP cases and controls. These differential responses formed three co-expression/co-methylation modules that were related to CRSwNP and three that were related to RV (Bonferroni corrected p < .01). Most (95%) of the differentially methylated CpGs (DMCs) were in modules related to CRSwNP, whereas the differentially expressed genes (DEGs) were more equally distributed between the CRSwNP- and RV-related modules. Genes in the CRSwNP-related modules were enriched in known CRS and/or viral response immune pathways. CONCLUSION: RV activates specific epigenetic programs and correlated transcriptional networks in the sinonasal epithelium of individuals with CRSwNP. These novel observations suggest epigenetic signatures specific to patients with CRSwNP modulate response to viral pathogens at the mucosal environmental interface. Determining how viral response pathways are involved in epithelial inflammation in CRSwNP could lead to therapeutic targets for this burdensome airway disorder.

Multi-omics colocalization with genome-wide association studies reveals a context-specific genetic mechanism at a childhood onset asthma risk locus.

BACKGROUND: Genome-wide association studies (GWASs) have identified thousands of variants associated with asthma and other complex diseases. However, the functional effects of most of these variants are unknown. Moreover, GWASs do not provide context-specific information on cell types or environmental factors that affect specific disease risks and outcomes. To address these limitations, we used an upper airway epithelial cell (AEC) culture model to assess transcriptional and epigenetic responses to rhinovirus (RV), an asthma-promoting pathogen, and provide context-specific functional annotations to variants discovered in GWASs of asthma. METHODS: Genome-wide genetic, gene expression, and DNA methylation data in vehicle- and RV-treated upper AECs were collected from 104 individuals who had a diagnosis of airway disease (n=66) or were healthy participants (n=38). We mapped cis expression and methylation quantitative trait loci (cis-eQTLs and cis-meQTLs, respectively) in each treatment condition (RV and vehicle) in AECs from these individuals. A Bayesian test for colocalization between AEC molecular QTLs and adult onset asthma and childhood onset asthma GWAS SNPs, and a multi-ethnic GWAS of asthma, was used to assign the function to variants associated with asthma. We used Mendelian randomization to demonstrate DNA methylation effects on gene expression at asthma colocalized loci. RESULTS: Asthma and allergic disease-associated GWAS SNPs were specifically enriched among molecular QTLs in AECs, but not in GWASs from non-immune diseases, and in AEC eQTLs, but not among eQTLs from other tissues. Colocalization analyses of AEC QTLs with asthma GWAS variants revealed potential molecular mechanisms of asthma, including QTLs at the TSLP locus that were common to both the RV and vehicle treatments and to both childhood onset and adult onset asthma, as well as QTLs at the 17q12-21 asthma locus that were specific to RV exposure and childhood onset asthma, consistent with clinical and epidemiological studies of these loci. CONCLUSIONS: This study provides evidence of functional effects for asthma risk variants in AECs and insight into RV-mediated transcriptional and epigenetic response mechanisms that modulate genetic effects in the airway and risk for asthma.