Primary nasal influenza infection rewires tissue-scale memory response dynamics.

The nasal mucosa is often the initial site of respiratory viral infection, replication, and transmission. Understanding how infection shapes tissue-scale primary and memory responses is critical for designing mucosal therapeutics and vaccines. We generated a single-cell RNA-sequencing atlas of the murine nasal mucosa, sampling three regions during primary influenza infection and rechallenge. Compositional analysis revealed restricted infection to the respiratory mucosa with stepwise changes in immune and epithelial cell subsets and states. We identified and characterized a rare subset of Krt13+ nasal immune-interacting floor epithelial (KNIIFE) cells, which concurrently increased with tissue-resident memory T (TRM)-like cells. Proportionality analysis, cell-cell communication inference, and microscopy underscored the CXCL16-CXCR6 axis between KNIIFE and TRM cells. Secondary influenza challenge induced accelerated and coordinated myeloid and lymphoid responses without epithelial proliferation. Together, this atlas serves as a reference for viral infection in the upper respiratory tract and highlights the efficacy of local coordinated memory responses.

DNA methylation signatures in airway cells from adult children of asthmatic mothers reflect subtypes of severe asthma.

Maternal asthma (MA) is among the most consistent risk factors for asthma in children. Possible mechanisms for this observation are epigenetic modifications in utero that have lasting effects on developmental programs in children of mothers with asthma. To test this hypothesis, we performed differential DNA methylation analyses of 398,186 individual CpG sites in primary bronchial epithelial cells (BECs) from 42 nonasthma controls and 88 asthma cases, including 56 without MA (NMA) and 32 with MA. We used weighted gene coexpression network analysis (WGCNA) of 69 and 554 differentially methylated CpGs (DMCs) that were specific to NMA and MA cases, respectively, compared with controls. WGCNA grouped 66 NMA-DMCs and 203 MA-DMCs into two and five comethylation modules, respectively. The eigenvector of one MA-associated module (turquoise) was uniquely correlated with 85 genes expressed in BECs and enriched for 36 pathways, 16 of which discriminated between NMA and MA using machine learning. Genes in all 16 pathways were decreased in MA compared with NMA cases (P = 7.1 × 10−3), a finding that replicated in nasal epithelial cells from an independent cohort (P = 0.02). Functional interpretation of these pathways suggested impaired T cell signaling and responses to viral and bacterial pathogens. The MA-associated turquoise module eigenvector was additionally correlated with clinical features of severe asthma and reflective of type 2 (T2)-low asthma (i.e., low total serum immunoglobulin E, fractional exhaled nitric oxide, and eosinophilia). Overall, these data suggest that MA alters diverse epigenetically mediated pathways that lead to distinct subtypes of severe asthma in adults, including hard-to-treat T2-low asthma.

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.

A-to-I editing of miR-200b-3p in airway cells is associated with moderate-to-severe asthma.

BACKGROUND: Asthma is a chronic lung disease characterised by persistent airway inflammation. Altered microRNA (miRNA)-mediated gene silencing in bronchial epithelial cells (BECs) has been reported in asthma, yet adenosine deaminase acting on RNA (ADAR)-mediated miRNA editing in asthma remains unexplored. METHODS: We first identified adenosine to inosine (A-to-I) edited sites in miRNAs in BECs from 142 adult asthma cases and controls. A-to-I edited sites were tested for associations with asthma severity and clinical measures of asthma. Paired RNA sequencing data were used to perform pathway enrichments and test for associations with bioinformatically predicted target genes of the unedited and edited miRNAs. RESULTS: Of 19 A-to-I edited sites detected in these miRNAs, one site at position 5 of miR-200b-3p was edited less frequently in cases compared with controls (pcorrected=0.013), and especially compared with cases with moderate (pcorrected=0.029) and severe (pcorrected=3.9×10-4), but not mild (pcorrected=0.38), asthma. Bioinformatic prediction revealed 232 target genes of the edited miR-200b-3p, which were enriched for both interleukin-4 and interferon-γ signalling pathways, and included the SOCS1 (suppressor of cytokine signalling 1) gene. SOCS1 was more highly expressed in moderate (pcorrected=0.017) and severe (pcorrected=5.4×10-3) asthma cases compared with controls. Moreover, both miR-200b-3p editing and SOCS1 were associated with bronchoalveolar lavage eosinophil levels. CONCLUSIONS: Reduced A-to-I editing of position 5 of miR-200b-3p in lower airway cells from moderate-to-severe asthmatic subjects may lead to overexpression of SOCS1 and impaired cytokine signalling. We propose ADAR-mediated editing as an epigenetic mechanism contributing to features of moderate-to-severe asthma in adulthood.

Epigenetic landscape links upper airway microbiota in infancy with allergic rhinitis at 6 years of age.

BACKGROUND: The upper airways present a barrier to inhaled allergens and microbes, which alter immune responses and subsequent risk for diseases, such as allergic rhinitis (AR). OBJECTIVE: We tested the hypothesis that early-life microbial exposures leave a lasting signature in DNA methylation that ultimately influences the development of AR in children. METHODS: We studied upper airway microbiota at 1 week, 1 month, and 3 months of life, and measured DNA methylation and gene expression profiles in upper airway mucosal cells and assessed AR at age 6 years in children in the Copenhagen Prospective Studies on Asthma in Childhood birth cohort. RESULTS: We identified 956 AR-associated differentially methylated CpGs in upper airway mucosal cells at age 6 years, 792 of which formed 3 modules of correlated differentially methylated CpGs. The eigenvector of 1 module was correlated with the expression of genes enriched for lysosome and bacterial invasion of epithelial cell pathways. Early-life microbial diversity was lower at 1 week (richness P = .0079) in children with AR at age 6 years, and reduced diversity at 1 week was also correlated with the same module's eigenvector (ρ = -0.25; P = 3.3 × 10-5). We show that the effect of microbiota richness at 1 week on risk for AR at age 6 years was mediated in part by the epigenetic signature of this module. CONCLUSIONS: Our results suggest that upper airway microbial composition in infancy contributes to the development of AR during childhood, and this trajectory is mediated, at least in part, through altered DNA methylation patterns in upper airway mucosal cells.

Ultra-fast local-haplotype variant calling using paired-end DNA-sequencing data reveals somatic mosaicism in tumor and normal blood samples.

Somatic mosaicism refers to the existence of somatic mutations in a fraction of somatic cells in a single biological sample. Its importance has mainly been discussed in theory although experimental work has started to emerge linking somatic mosaicism to disease diagnosis. Through novel statistical modeling of paired-end DNA-sequencing data using blood-derived DNA from healthy donors as well as DNA from tumor samples, we present an ultra-fast computational pipeline, LocHap that searches for multiple single nucleotide variants (SNVs) that are scaffolded by the same reads. We refer to scaffolded SNVs as local haplotypes (LH). When an LH exhibits more than two genotypes, we call it a local haplotype variant (LHV). The presence of LHVs is considered evidence of somatic mosaicism because a genetically homogeneous cell population will not harbor LHVs. Applying LocHap to whole-genome and whole-exome sequence data in DNA from normal blood and tumor samples, we find wide-spread LHVs across the genome. Importantly, we find more LHVs in tumor samples than in normal samples, and more in older adults than in younger ones. We confirm the existence of LHVs and somatic mosaicism by validation studies in normal blood samples. LocHap is publicly available at http://www.compgenome.org/lochap.

Genome-Wide Methylation Study Identifies an IL-13-induced Epigenetic Signature in Asthmatic Airways.

RATIONALE: Epigenetic changes to airway cells have been proposed as important modulators of the effects of environmental exposures on airway diseases, yet no study to date has shown epigenetic responses to exposures in the airway that correlate with disease state. The type 2 cytokine IL-13 is a key mediator of allergic airway diseases, such as asthma, and is up-regulated in response to many asthma-promoting exposures. OBJECTIVES: To directly study the epigenetic response of airway epithelial cells (AECs) to IL-13 and test whether IL-13-induced epigenetic changes differ between individuals with and without asthma. METHODS: Genome-wide DNA methylation and gene expression patterns were studied in 58 IL-13-treated and untreated primary AEC cultures and validated in freshly isolated cells of subjects with and without asthma using the Illumina Human Methylation 450K and HumanHT-12 BeadChips. IL-13-mediated comethylation modules were identified and correlated with clinical phenotypes using weighted gene coexpression network analysis. MEASUREMENTS AND MAIN RESULTS: IL-13 altered global DNA methylation patterns in cultured AECs and were significantly enriched near genes associated with asthma. Importantly, a significant proportion of this IL-13 epigenetic signature was validated in freshly isolated AECs from subjects with asthma and clustered into two distinct modules, with module 1 correlated with asthma severity and lung function and module 2 with eosinophilia. CONCLUSIONS: These results suggest that a single exposure of IL-13 may selectively induce long-lasting DNA methylation changes in asthmatic airways that alter specific AEC pathways and contribute to asthma phenotypes.

Primary nasal viral infection rewires the tissue-scale memory response.

The nasal mucosa is frequently the initial site of respiratory viral infection, replication, and transmission. Recent work has started to clarify the independent responses of epithelial, myeloid, and lymphoid cells to viral infection in the nasal mucosa, but their spatiotemporal coordination and relative contributions remain unclear. Furthermore, understanding whether and how primary infection shapes tissue-scale memory responses to secondary challenge is critical for the rational design of nasal-targeting therapeutics and vaccines. Here, we generated a single-cell RNA-sequencing (scRNA-seq) atlas of the murine nasal mucosa sampling three distinct regions before and during primary and secondary influenza infection. Primary infection was largely restricted to respiratory mucosa and induced stepwise changes in cell type, subset, and state composition over time. Type I Interferon (IFN)-responsive neutrophils appeared 2 days post infection (dpi) and preceded transient IFN-responsive/cycling epithelial cell responses 5 dpi, which coincided with broader antiviral monocyte and NK cell accumulation. By 8 dpi, monocyte-derived macrophages (MDMs) expressing Cxcl9 and Cxcl16 arose alongside effector cytotoxic CD8 and Ifng-expressing CD4 T cells. Following viral clearance (14 dpi), rare, previously undescribed Krt13+ nasal immune-interacting floor epithelial (KNIIFE) cells expressing multiple genes with immune communication potential increased concurrently with tissue-resident memory T (TRM)-like cells and early IgG+/IgA+ plasmablasts. Proportionality analysis coupled with cell-cell communication inference, alongside validation by in situ microscopy, underscored the CXCL16-CXCR6 signaling axis between MDMs and effector CD8 T cells 8dpi and KNIIFE cells and TRM cells 14 dpi. Secondary influenza challenge with a homologous or heterologous strain administered 60 dpi induced an accelerated and coordinated myeloid and lymphoid response without epithelial proliferation, illustrating how tissue-scale memory to natural infection engages both myeloid and lymphoid cells to reduce epithelial regenerative burden. Together, this atlas serves as a reference for viral infection in the upper respiratory tract and highlights the efficacy of local coordinated memory responses upon rechallenge.

Sex-specific differences in peripheral blood leukocyte transcriptional response to LPS are enriched for HLA region and X chromosome genes.

Sex-specific differences in prevalence are well documented for many common, complex diseases, especially for immune-mediated diseases, yet the precise mechanisms through which factors associated with biological sex exert their effects throughout life are not well understood. We interrogated sex-specific transcriptional responses of peripheral blood leukocytes (PBLs) to innate immune stimulation by lipopolysaccharide (LPS) in 46 male and 66 female members of the Hutterite community, who practice a communal lifestyle. We identified 1217 autosomal and 54 X-linked genes with sex-specific responses to LPS, as well as 71 autosomal and one X-linked sex-specific expression quantitative trait loci (eQTLs). Despite a similar proportion of the 15 HLA genes responding to LPS compared to all expressed autosomal genes, there was a significant over-representation of genes with sex by treatment interactions among HLA genes. We also observed an enrichment of sex-specific differentially expressed genes in response to LPS for X-linked genes compared to autosomal genes, suggesting that HLA and X-linked genes may disproportionately contribute to sex disparities in risk for immune-mediated diseases.

Longitudinal data reveal strong genetic and weak non-genetic components of ethnicity-dependent blood DNA methylation levels.

Epigenetic architecture is influenced by genetic and environmental factors, but little is known about their relative contributions or longitudinal dynamics. Here, we studied DNA methylation (DNAm) at over 750,000 CpG sites in mononuclear blood cells collected at birth and age 7 from 196 children of primarily self-reported Black and Hispanic ethnicities to study race-associated DNAm patterns. We developed a novel Bayesian method for high-dimensional longitudinal data and showed that race-associated DNAm patterns at birth and age 7 are nearly identical. Additionally, we estimated that up to 51% of all self-reported race-associated CpGs had race-dependent DNAm levels that were mediated through local genotype and, quite surprisingly, found that genetic factors explained an overwhelming majority of the variation in DNAm levels at other, previously identified, environmentally-associated CpGs. These results indicate that race-associated blood DNAm patterns in particular, and blood DNAm levels in general, are primarily driven by genetic factors, and are not as sensitive to environmental exposures as previously suggested, at least during the first 7 years of life.

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.

Expression quantitative trait locus fine mapping of the 17q12-21 asthma locus in African American children: a genetic association and gene expression study.

BACKGROUND: African ancestry is associated with a higher prevalence and greater severity of asthma than European ancestries, yet genetic studies of the most common locus associated with childhood-onset asthma, 17q12-21, in African Americans have been inconclusive. The aim of this study was to leverage both the phenotyping of the Children's Respiratory and Environmental Workgroup (CREW) birth cohort consortium, and the reduced linkage disequilibrium in African Americans, to fine map the 17q12-21 locus. METHODS: We first did a genetic association study and meta-analysis using 17q12-21 tag single-nucleotide polymorphisms (SNPs) for childhood-onset asthma in 1613 European American and 870 African American children from the CREW consortium. Nine tag SNPs were selected based on linkage disequilibrium patterns at 17q12-21 and their association with asthma, considering the effect allele under an additive model (0, 1, or 2 effect alleles). Results were meta-analysed with publicly available summary data from the EVE consortium (on 4303 European American and 3034 African American individuals) for seven of the nine SNPs of interest. Subsequently, we tested for expression quantitative trait loci (eQTLs) among the SNPs associated with childhood-onset asthma and the expression of 17q12-21 genes in resting peripheral blood mononuclear cells (PBMCs) from 85 African American CREW children and in upper airway epithelial cells from 246 African American CREW children; and in lower airway epithelial cells from 44 European American and 72 African American adults from a case-control study of asthma genetic risk in Chicago (IL, USA). FINDINGS: 17q12-21 SNPs were broadly associated with asthma in European Americans. Only two SNPs (rs2305480 in gasdermin-B [GSDMB] and rs8076131 in ORMDL sphingolipid biosynthesis regulator 3 [ORMDL3]) were associated with asthma in African Americans, at a Bonferroni-corrected threshold of p<0·0055 (for rs2305480_G, odds ratio [OR] 1·36 [95% CI 1·12-1·65], p=0·0014; and for rs8076131_A, OR 1·37 [1·13-1·67], p=0·0010). In upper airway epithelial cells from African American children, genotype at rs2305480 was the most significant eQTL for GSDMB (eQTL effect size [ß] 1·35 [95% CI 1·25-1·46], p<0·0001), and to a lesser extent showed an eQTL effect for post-GPI attachment to proteins phospholipase 3 (ß 1·15 [1·08-1·22], p<0·0001). No SNPs were eQTLs for ORMDL3. By contrast, in PBMCs, the five core SNPs were associated only with expression of GSDMB and ORMDL3. Genotype at rs12936231 (in zona pellucida binding protein 2) showed the strongest associations across both genes (for GSDMB, eQTLß 1·24 [1·15-1·32], p<0·0001; and for ORMDL3 (ß 1·19 [1·12-1·24], p<0·0001). The eQTL effects of rs2305480 on GSDMB expression were replicated in lower airway cells from African American adults (ß 1·29 [1·15-1·44], p<0·0001). INTERPRETATION: Our study suggests that SNPs regulating GSDMB expression in airway epithelial cells have a major role in childhood-onset asthma, whereas SNPs regulating the expression levels of 17q12-21 genes in resting blood cells are not central to asthma risk. Our genetic and gene expression data in African Americans and European Americans indicated GSDMB to be the leading candidate gene at this important asthma locus. FUNDING: National Institutes of Health, Office of the Director.