Sex-stratified and interaction multi-omics analysis of DNA methylation and gene expression in asthma

Background: Sex differences in asthma are well recognized, but sex-associated airway DNA methylation (DNAm) patterns and their relationships with gene expression remain incompletely characterized. We therefore conducted cis- and trans-expression quantitative trait methylation (eQTM) analyses using sex-stratified and CpG-by-sex interaction models. Methods: We analyzed nasal airway epithelial DNAm and transcriptomic data from 298 participants in the CAAPA discovery cohort. Asthma-associated differentially methylated CpGs (DMCs) were identified using sex-combined and sex-stratified models. Cis- and trans-eQTM analyses were then performed in participants with asthma using sex-combined, sex-stratified, and CpG-by-sex interaction models. To prioritize asthma-associated regulatory relationships, eQTMs that were also significant in non-asthmatic controls were excluded from downstream analyses. Replication in an independent bronchial epithelial cohort was defined by concordant effect direction and P < 0.05. Replicated targets underwent pathway and network analyses. Findings: Asthma-associated DMCs were predominantly hypomethylated in females and hypermethylated in males. Among trans-eQTMs, 379 associations were shared across all four models, whereas the majority were model-specific, indicating extensive sex-dependent regulatory architecture. We replicated 3, 1, 9, and 2 cis-eQTMs and 96, 84, 216, and 41 trans-eQTMs in the sex-combined, female-stratified, male-stratified, and interaction models, respectively. Sex-combined targets implicated innate immune and inflammatory signaling (TNF, IRF8); female-stratified targets implicated cytoskeletal remodeling, focal adhesion, and PI3K-Akt signaling (ITGA10); male-stratified targets implicated transcriptional regulation and post-translational modification (MOV10, HDAC7, CSF1, KAT2A); and interaction targets implicated lipid-mediator metabolism and inflammatory signaling and resolution (CYP2C9, LTB4R, CCL19). Our sex-dependent analyses showed that the majority of eQTMs resided within the trans-regulatory architecture, suggesting that sex-dependent gene expression differences are associated predominantly with trans rather than cis regulatory mechanisms. Pathway analyses revealed enrichment of innate immunity in the sex-combined model, structural remodeling in females, transcriptional control in males, and inflammatory resolution in the sex-interaction model. Interpretation: Sex-specific and sex-dependent airway methylation-expression associations were largely undetected by sex-combined analysis. Sex-aware analytical models may help resolve epigenetic mechanisms underlying sex differences in asthma biology.