Multi-omics in nasal epithelium reveals three axes of dysregulation for asthma risk in the African Diaspora populations

Background: Asthma, a complex chronic lung disease affecting the airways, has striking disparities across ancestral groups, but the molecular underpinning of these differences is poorly understood and minimally studied. A major goal of the Consortium on Asthma among African-ancestry Populations in the Americas (CAAPA) is to understand multi-omics signatures of asthma risk in the nasal epithelium focusing on populations of African ancestry. Methods: RNASeq data were generated from nasal epithelium in subjects recruited from up to 7 sites (Baltimore, Washington DC, Chicago, Denver, Salvador Brazil, Barbados, and Nigeria). Current asthma cases (N=253) were compared to never-asthma controls (N=283) to identify differentially expressed genes (DEGs; q <0.05). Network analyses were performed with Ingenuity Pathway Analysis (IPA; DEGs with q<0.05) and weighted gene co-expression network analysis (WGCNA; DEGs with q<0.15). All models were fully adjusted for ancestry, sampling site, and appropriate latent factors. Findings: CAAPA represents diversity across the African Diaspora with a wide range of continental African ancestry (9%-100%). We identified 389 DEGs; the top DEG, FN1, was downregulated in asthma cases (q=3.26x10-9) and encodes fibronectin which plays a role in wound healing. Others in the top 10 DEGs have high relevance for asthma: SNTG2 (q = 1.12x10-4) is the target of multiple miRNAs related to asthma; PPP1R9A expression (q=7.60x10-5) was previously determined to be influenced by IL-13 in mouse lung; and SPTBN1 (q=1.12x10-4) plays a key role in mediating TGFβ signaling. IPA revealed networks with upstream regulators relevant for immune response (IL4; p=7.25x10-10 and TGFβ1; p=5.47x10-8) and drug response (dexamethasone; p=4.31x10-10 and fluticasone propionate; p=9.42x10-8). Among asthma cases, genes regulated by dexamethasone and fluticasone propionate were not associated with inhaled corticosteroid medication use. The top three WGCNA modules implicate networks related to immune response (CEACAM5; p=9.62x10-16 and CPA3; p=2.39x10-14) and wound healing (FN1; p=7.63x10-9). Multi-omic analysis identified FKBP5 as a key contributor to asthma risk, whereby the association between nasal epithelium gene expression is mediated through methylation and is associated with increased use of inhaled corticosteroids. FKBP5 is a co-chaperone of glucocorticoid receptor signaling and known to be involved in drug response in asthma. Interpretation: Our analyses reveal genes and networks in asthma that are differentially expressed in nasal epithelium of asthma cases of African ancestry in CAAPA. Importantly, this work reveals molecular dysregulation on three axes – increased Th2 inflammation, decreased capacity for wound healing, and impaired drug response – that may play a critical role in asthma within the African Diaspora. RNA sequencing data was generated from nasal epithelium samples in newly recruited subjects representing 7 locations across the African Diaspora: Baltimore, Washington DC, Chicago and Denver in the United States; Salvador, Brazil; Barbados in the Caribbean; and Nigeria in West Africa. We examined transcriptomic differences between asthma cases and never-asthma controls to identify differentially expressed individual genes and networks of co-expressed genes.