Integrated longitudinal multiomics study identifies immune programs associated with acute COVID-19 severity and mortality

BACKGROUNDPatients hospitalized for COVID-19 exhibit diverse clinical outcomes, with outcomes for some individuals diverging over time even though their initial disease severity appears similar to that of other patients. A systematic evaluation of molecular and cellular profiles over the full disease course can link immune programs and their coordination with progression heterogeneity.METHODSWe performed deep immunophenotyping and conducted longitudinal multiomics modeling, integrating 10 assays for 1,152 Immunophenotyping Assessment in a COVID-19 Cohort (IMPACC) study participants and identifying several immune cascades that were significant drivers of differential clinical outcomes.RESULTSIncreasing disease severity was driven by a temporal pattern that began with the early upregulation of immunosuppressive metabolites and then elevated levels of inflammatory cytokines, signatures of coagulation, formation of neutrophil extracellular traps, and T cell functional dysregulation. A second immune cascade, predictive of 28-day mortality among critically ill patients, was characterized by reduced total plasma Igs and B cells and dysregulated IFN responsiveness. We demonstrated that the balance disruption between IFN-stimulated genes and IFN inhibitors is a crucial biomarker of COVID-19 mortality, potentially contributing to failure of viral clearance in patients with fatal illness.CONCLUSIONOur longitudinal multiomics profiling study revealed temporal coordination across diverse omics that potentially explain the disease progression, providing insights that can inform the targeted development of therapies for patients hospitalized with COVID-19, especially those who are critically ill.TRIAL REGISTRATIONClinicalTrials.gov NCT04378777.FUNDINGNIH (5R01AI135803-03, 5U19AI118608-04, 5U19AI128910-04, 4U19AI090023-11, 4U19AI118610-06, R01AI145835-01A1S1, 5U19AI062629-17, 5U19AI057229-17, 5U19AI125357-05, 5U19AI128913-03, 3U19AI077439-13, 5U54AI142766-03, 5R01AI104870-07, 3U19AI089992-09, 3U19AI128913-03, and 5T32DA018926-18); NIAID, NIH (3U19AI1289130, U19AI128913-04S1, and R01AI122220); and National Science Foundation (DMS2310836).

Host-microbe multiomic profiling reveals age-dependent immune dysregulation associated with COVID-19 immunopathology

Age is a major risk factor for severe coronavirus disease 2019 (COVID-19), yet the mechanisms behind this relationship have remained incompletely understood. To address this, we evaluated the impact of aging on host immune response in the blood and the upper airway, as well as the nasal microbiome in a prospective, multicenter cohort of 1031 vaccine-naïve patients hospitalized for COVID-19 between 18 and 96 years old. We performed mass cytometry, serum protein profiling, anti-severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) antibody assays, and blood and nasal transcriptomics. We found that older age correlated with increased SARS-CoV-2 viral abundance upon hospital admission, delayed viral clearance, and increased type I interferon gene expression in both the blood and upper airway. We also observed age-dependent up-regulation of innate immune signaling pathways and down-regulation of adaptive immune signaling pathways. Older adults had lower naïve T and B cell populations and higher monocyte populations. Over time, older adults demonstrated a sustained induction of pro-inflammatory genes and serum chemokines compared with younger individuals, suggesting an age-dependent impairment in inflammation resolution. Transcriptional and protein biomarkers of disease severity differed with age, with the oldest adults exhibiting greater expression of pro-inflammatory genes and proteins in severe disease. Together, our study finds that aging is associated with impaired viral clearance, dysregulated immune signaling, and persistent and potentially pathologic activation of pro-inflammatory genes and proteins.

Analytical challenges in omics research on asthma and allergy: A National Institute of Allergy and Infectious Diseases workshop

Studies of asthma and allergy are generating increasing volumes of omics data for analysis and interpretation. The National Institute of Allergy and Infectious Diseases (NIAID) assembled a workshop comprising investigators studying asthma and allergic diseases using omics approaches, omics investigators from outside the field, and NIAID medical and scientific officers to discuss the following areas in asthma and allergy research: genomics, epigenomics, transcriptomics, microbiomics, metabolomics, proteomics, lipidomics, integrative omics, systems biology, and causal inference. Current states of the art, present challenges, novel and emerging strategies, and priorities for progress were presented and discussed for each area. This workshop report summarizes the major points and conclusions from this NIAID workshop. As a group, the investigators underscored the imperatives for rigorous analytic frameworks, integration of different omics data types, cross-disciplinary interaction, strategies for overcoming current limitations, and the overarching goal to improve scientific understanding and care of asthma and allergic diseases.

Activated sputum eosinophils associated with exacerbations in children on mepolizumab

BACKGROUND

MUPPITS-2 was a randomized, placebo-controlled clinical trial that demonstrated mepolizumab (anti-IL-5) reduced exacerbations and blood and airway eosinophils in urban children with severe eosinophilic asthma. Despite this reduction in eosinophilia, exacerbation risk persisted in certain patients treated with mepolizumab. This raises the possibility that subpopulations of airway eosinophils exist that contribute to breakthrough exacerbations.

OBJECTIVE

We aimed to determine the effect of mepolizumab on airway eosinophils in childhood asthma.

METHODS

Sputum samples were obtained from 53 MUPPITS-2 participants. Airway eosinophils were characterized using mass cytometry and grouped into subpopulations using unsupervised clustering analyses of 38 surface and intracellular markers. Differences in frequency and immunophenotype of sputum eosinophil subpopulations were assessed based on treatment arm and frequency of exacerbations.

RESULTS

Median sputum eosinophils were significantly lower among participants treated with mepolizumab compared with placebo (58% lower, 0.35% difference [95% CI 0.01, 0.74], P = .04). Clustering analysis identified 3 subpopulations of sputum eosinophils with varied expression of CD62L. CD62Lint and CD62Lhi eosinophils exhibited significantly elevated activation marker and eosinophil peroxidase expression, respectively. In mepolizumab-treated participants, CD62Lint and CD62Lhi eosinophils were more abundant in participants who experienced exacerbations than in those who did not (100% higher for CD62Lint, 0.04% difference [95% CI 0.0, 0.13], P = .04; 93% higher for CD62Lhi, 0.21% difference [95% CI 0.0, 0.77], P = .04).

CONCLUSIONS

Children with eosinophilic asthma treated with mepolizumab had significantly lower sputum eosinophils. However, CD62Lint and CD62Lhi eosinophils were significantly elevated in children on mepolizumab who had exacerbations, suggesting that eosinophil subpopulations exist that contribute to exacerbations despite anti-IL-5 treatment.

Host-Microbe Multiomic Profiling Reveals Age-Dependent COVID-19 Immunopathology

Age is a major risk factor for severe coronavirus disease-2019 (COVID-19), yet the mechanisms responsible for this relationship have remained incompletely understood. To address this, we evaluated the impact of aging on host and viral dynamics in a prospective, multicenter cohort of 1,031 patients hospitalized for COVID-19, ranging from 18 to 96 years of age. We performed blood transcriptomics and nasal metatranscriptomics, and measured peripheral blood immune cell populations, inflammatory protein expression, anti-SARS-CoV-2 antibodies, and anti-interferon (IFN) autoantibodies. We found that older age correlated with an increased SARS-CoV-2 viral load at the time of admission, and with delayed viral clearance over 28 days. This contributed to an age-dependent increase in type I IFN gene expression in both the respiratory tract and blood. We also observed age-dependent transcriptional increases in peripheral blood IFN-γ, neutrophil degranulation, and Toll like receptor (TLR) signaling pathways, and decreases in T cell receptor (TCR) and B cell receptor signaling pathways. Over time, older adults exhibited a remarkably sustained induction of proinflammatory genes (e.g., CXCL6) and serum chemokines (e.g., CXCL9) compared to younger individuals, highlighting a striking age-dependent impairment in inflammation resolution. Augmented inflammatory signaling also involved the upper airway, where aging was associated with upregulation of TLR, IL17, type I IFN and IL1 pathways, and downregulation TCR and PD-1 signaling pathways. Metatranscriptomics revealed that the oldest adults exhibited disproportionate reactivation of herpes simplex virus and cytomegalovirus in the upper airway following hospitalization. Mass cytometry demonstrated that aging correlated with reduced naïve T and B cell populations, and increased monocytes and exhausted natural killer cells. Transcriptional and protein biomarkers of disease severity markedly differed with age, with the oldest adults exhibiting greater expression of TLR and inflammasome signaling genes, as well as proinflammatory proteins (e.g., IL6, CXCL8), in severe COVID-19 compared to mild/moderate disease. Anti-IFN autoantibody prevalence correlated with both age and disease severity. Taken together, this work profiles both host and microbe in the blood and airway to provide fresh insights into aging-related immune changes in a large cohort of vaccine-naïve COVID-19 patients. We observed age-dependent immune dysregulation at the transcriptional, protein and cellular levels, manifesting in an imbalance of inflammatory responses over the course of hospitalization, and suggesting potential new therapeutic targets.

IgM N-glycosylation correlates with COVID-19 severity and rate of complement deposition

The glycosylation of IgG plays a critical role during human severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, activating immune cells and inducing cytokine production. However, the role of IgM N-glycosylation has not been studied during human acute viral infection. The analysis of IgM N-glycosylation from healthy controls and hospitalized coronavirus disease 2019 (COVID-19) patients reveals increased high-mannose and sialylation that correlates with COVID-19 severity. These trends are confirmed within SARS-CoV-2-specific immunoglobulin N-glycan profiles. Moreover, the degree of total IgM mannosylation and sialylation correlate significantly with markers of disease severity. We link the changes of IgM N-glycosylation with the expression of Golgi glycosyltransferases. Lastly, we observe antigen-specific IgM antibody-dependent complement deposition is elevated in severe COVID-19 patients and modulated by exoglycosidase digestion. Taken together, this work links the IgM N-glycosylation with COVID-19 severity and highlights the need to understand IgM glycosylation and downstream immune function during human disease.

Features of acute COVID-19 associated with post-acute sequelae of SARS-CoV-2 phenotypes: results from the IMPACC study

Post-acute sequelae of SARS-CoV-2 (PASC) is a significant public health concern. We describe Patient Reported Outcomes (PROs) on 590 participants prospectively assessed from hospital admission for COVID-19 through one year after discharge. Modeling identified 4 PRO clusters based on reported deficits (minimal, physical, mental/cognitive, and multidomain), supporting heterogenous clinical presentations in PASC, with sub-phenotypes associated with female sex and distinctive comorbidities. During the acute phase of disease, a higher respiratory SARS-CoV-2 viral burden and lower Receptor Binding Domain and Spike antibody titers were associated with both the physical predominant and the multidomain deficit clusters. A lower frequency of circulating B lymphocytes by mass cytometry (CyTOF) was observed in the multidomain deficit cluster. Circulating fibroblast growth factor 21 (FGF21) was significantly elevated in the mental/cognitive predominant and the multidomain clusters. Future efforts to link PASC to acute anti-viral host responses may help to better target treatment and prevention of PASC.

Molecular pathways underlying lung-brain axis signaling in asthma: Relevance for psychopathology and neuroinflammation

BACKGROUND

Accumulating evidence indicates that asthma has systemic effects and affects brain function. Although airway inflammation is proposed to initiate afferent communications with the brain, the signaling pathways have not been established.

OBJECTIVE

We sought to identify the cellular and molecular pathways involved in afferent lung-brain communication during airway inflammation in asthma.

METHODS

In 23 adults with mild asthma, segmental bronchial provocation with allergen (SBP-Ag) was used to provoke airway inflammation and retrieve bronchoalveolar lavage fluid for targeted protein analysis and RNA sequencing to determine gene expression profiles. Neural responses to emotional cues in nodes of the salience network were assessed with functional magnetic resonance imaging at baseline and 48 hours after SBP-Ag.

RESULTS

Cell deconvolution and gene coexpression network analysis identified 11 cell-associated gene modules that changed in response to SBP-Ag. SBP-Ag increased bronchoalveolar lavage eosinophils and expression of an eosinophil-associated module enriched for genes related to TH17-type inflammation (eg, IL17A), as well as cell proliferation in lung and brain (eg, NOTCH1, VEGFA, and LIF). Increased expression of genes in this module, as well as several TH17-type inflammation-related proteins, was associated with an increase from baseline in salience network reactivity.

CONCLUSIONS

Our results identify a specific inflammatory pathway linking asthma-related airway inflammation and emotion-related neural function. Systemically, TH17-type inflammation has been implicated in both depression and neuroinflammation, with impacts on long-term brain health. Thus, our data emphasize that inflammation in the lung in asthma may have profound effects outside of the lung that may be targetable with novel therapeutic approaches.

Integrated longitudinal multi-omics study identifies immune programs associated with COVID-19 severity and mortality in 1152 hospitalized participants

Hospitalized COVID-19 patients exhibit diverse clinical outcomes, with some individuals diverging over time even though their initial disease severity appears similar. A systematic evaluation of molecular and cellular profiles over the full disease course can link immune programs and their coordination with progression heterogeneity. In this study, we carried out deep immunophenotyping and conducted longitudinal multi-omics modeling integrating ten distinct assays on a total of 1,152 IMPACC participants and identified several immune cascades that were significant drivers of differential clinical outcomes. Increasing disease severity was driven by a temporal pattern that began with the early upregulation of immunosuppressive metabolites and then elevated levels of inflammatory cytokines, signatures of coagulation, NETosis, and T-cell functional dysregulation. A second immune cascade, predictive of 28-day mortality among critically ill patients, was characterized by reduced total plasma immunoglobulins and B cells, as well as dysregulated IFN responsiveness. We demonstrated that the balance disruption between IFN-stimulated genes and IFN inhibitors is a crucial biomarker of COVID-19 mortality, potentially contributing to the failure of viral clearance in patients with fatal illness. Our longitudinal multi-omics profiling study revealed novel temporal coordination across diverse omics that potentially explain disease progression, providing insights that inform the targeted development of therapies for hospitalized COVID-19 patients, especially those critically ill.

Nasal and blood transcriptomic pathways underpinning the clinical response to grass pollen immunotherapy

BACKGROUND

Allergen immunotherapy (AIT) is a well-established disease-modifying therapy for allergic rhinitis, yet the fundamental mechanisms underlying its clinical effect remain inadequately understood. Gauging Response in Allergic Rhinitis to Sublingual and Subcutaneous Immunotherapy was a randomized, double-blind, placebo-controlled trial of individuals allergic to timothy grass who received 2 years of placebo (n = 30), subcutaneous immunotherapy (SCIT) (n = 27), or sublingual immunotherapy (SLIT) (n = 27) and were then followed for 1 additional year.

OBJECTIVE

We used yearly biospecimens from the Gauging Response in Allergic Rhinitis to Sublingual and Subcutaneous Immunotherapy study to identify molecular mechanisms of response.

METHODS

We used longitudinal transcriptomic profiling of nasal brush and PBMC samples after allergen provocation to uncover airway and systemic expression pathways mediating responsiveness to AIT.

TRIAL REGISTRATION

ClinicalTrials.gov Identifier: NCT01335139, EudraCT Number: 2010-023536-16.

RESULTS

SCIT and SLIT demonstrated similar changes in gene module expression over time. In nasal samples, alterations included downregulation of pathways of mucus hypersecretion, leukocyte migration/activation, and endoplasmic reticulum stress (log2 fold changes -0.133 to -0.640, false discovery rates [FDRs] <0.05). We observed upregulation of modules related to epithelial development, junction formation, and lipid metabolism (log2 fold changes 0.104 to 0.393, FDRs <0.05). In PBMCs, modules related to cellular stress response and type 2 cytokine signaling were reduced by immunotherapy (log2 fold changes -0.611 to -0.828, FDRs <0.05). Expression of these modules was also significantly associated with both Total Nasal Symptom Score and peak nasal inspiratory flow, indicating important links between treatment, module expression, and allergen response.

CONCLUSIONS

Our results identify specific molecular responses of the nasal airway impacting barrier function, leukocyte migration activation, and mucus secretion that are affected by both SCIT and SLIT, offering potential targets to guide novel strategies for AIT.

Type 2 inflammation reduces SARS-CoV-2 replication in the airway epithelium in allergic asthma through functional alteration of ciliated epithelial cells

BACKGROUND

Despite well-known susceptibilities to other respiratory viral infections, individuals with allergic asthma have shown reduced susceptibility to severe coronavirus disease 2019 (COVID-19).

OBJECTIVE

We sought to identify mechanisms whereby type 2 inflammation in the airway protects against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) by using bronchial airway epithelial cells (AECs) from aeroallergen-sensitized children with asthma and healthy nonsensitized children.

METHODS

We measured SARS-CoV-2 replication and ACE2 protein and performed bulk and single-cell RNA sequencing of ex vivo infected AEC samples with SARS-CoV-2 infection and with or without IL-13 treatment.

RESULTS

We observed that viral replication was lower in AECs from children with allergic asthma than those from in healthy nonsensitized children and that IL-13 treatment reduced viral replication only in children with allergic asthma and not in healthy children. Lower viral transcript levels were associated with a downregulation of functional pathways of the ciliated epithelium related to differentiation as well as cilia and axoneme production and function, rather than lower ACE2 expression or increases in goblet cells or mucus secretion pathways. Moreover, single-cell RNA sequencing identified specific subsets of relatively undifferentiated ciliated epithelium (which are common in allergic asthma and highly responsive to IL-13) that directly accounted for impaired viral replication.

CONCLUSION

Our results identify a novel mechanism of innate protection against SARS-CoV-2 in allergic asthma that provides important molecular and clinical insights during the ongoing COVID-19 pandemic.

Airway transcriptome networks identify susceptibility to frequent asthma exacerbations in children

BACKGROUND

Frequent asthma exacerbators, defined as those experiencing more than 1 hospitalization in a year for an asthma exacerbation, represent an important subgroup of individuals with asthma. However, this group remains poorly defined and understudied in children.

OBJECTIVE

Our aim was to determine the molecular mechanisms underlying asthma pathogenesis and exacerbation frequency.

METHODS

We performed RNA sequencing of upper airway cells from both frequent and nonfrequent exacerbators enrolled in the Ohio Pediatric Asthma Repository.

RESULTS

Through molecular network analysis, we found that nonfrequent exacerbators display an increase in modules enriched for immune system processes, including type 2 inflammation and response to infection. In contrast, frequent exacerbators showed expression of modules enriched for nervous system processes, such as synaptic formation and axonal outgrowth.

CONCLUSION

These data suggest that the upper airway of frequent exacerbators undergoes peripheral nervous system remodeling, representing a novel mechanism underlying pediatric asthma exacerbation.

Multi-omic longitudinal study reveals immune correlates of clinical course among hospitalized COVID-19 patients

The IMPACC cohort, composed of >1,000 hospitalized COVID-19 participants, contains five illness trajectory groups (TGs) during acute infection (first 28 days), ranging from milder (TG1-3) to more severe disease course (TG4) and death (TG5). Here, we report deep immunophenotyping, profiling of >15,000 longitudinal blood and nasal samples from 540 participants of the IMPACC cohort, using 14 distinct assays. These unbiased analyses identify cellular and molecular signatures present within 72 h of hospital admission that distinguish moderate from severe and fatal COVID-19 disease. Importantly, cellular and molecular states also distinguish participants with more severe disease that recover or stabilize within 28 days from those that progress to fatal outcomes (TG4 vs. TG5). Furthermore, our longitudinal design reveals that these biologic states display distinct temporal patterns associated with clinical outcomes. Characterizing host immune responses in relation to heterogeneity in disease course may inform clinical prognosis and opportunities for intervention.

Distinct Epithelial-Innate Immune Cell Transcriptional Circuits Underlie Airway Hyperresponsiveness in Asthma

Rationale: Indirect airway hyperresponsiveness (AHR) is a highly specific feature of asthma, but the underlying mechanisms responsible for driving indirect AHR remain incompletely understood. Objectives: To identify differences in gene expression in epithelial brushings obtained from individuals with asthma who were characterized for indirect AHR in the form of exercise-induced bronchoconstriction (EIB). Methods: RNA-sequencing analysis was performed on epithelial brushings obtained from individuals with asthma with EIB (n = 11) and without EIB (n = 9). Differentially expressed genes (DEGs) between the groups were correlated with measures of airway physiology, sputum inflammatory markers, and airway wall immunopathology. On the basis of these relationships, we examined the effects of primary airway epithelial cells (AECs) and specific epithelial cell-derived cytokines on both mast cells (MCs) and eosinophils (EOS). Measurements and Main Results: We identified 120 DEGs in individuals with and without EIB. Network analyses suggested critical roles for IL-33-, IL-18-, and IFN-γ-related signaling among these DEGs. IL1RL1 expression was positively correlated with the density of MCs in the epithelial compartment, and IL1RL1, IL18R1, and IFNG were positively correlated with the density of intraepithelial EOS. Subsequent ex vivo modeling demonstrated that AECs promote sustained type 2 (T2) inflammation in MCs and enhance IL-33-induced T2 gene expression. Furthermore, EOS increase the expression of IFNG and IL13 in response to both IL-18 and IL-33 as well as exposure to AECs. Conclusions: Circuits involving epithelial interactions with MCs and EOS are closely associated with indirect AHR. Ex vivo modeling indicates that epithelial-dependent regulation of these innate cells may be critical in indirect AHR and modulating T2 and non-T2 inflammation in asthma.

Rhinovirus infection of the airway epithelium enhances mast cell immune responses via epithelial-derived interferons

BACKGROUND

Mast cells (MCs) within the airway epithelium in asthma are closely related to airway dysfunction, but cross talk between airway epithelial cells (AECs) and MCs in asthma remains incompletely understood. Human rhinovirus (RV) infections are key triggers for asthma progression, and AECs from individuals with asthma may have dysregulated antiviral responses.

OBJECTIVE

We utilized primary AECs in an ex vivo coculture model system to examine cross talk between AECs and MCs after epithelial rhinovirus infection.

METHODS

Primary AECs were obtained from 11 children with asthma and 10 healthy children, differentiated at air-liquid interface, and cultured in the presence of laboratory of allergic diseases 2 (LAD2) MCs. AECs were infected with rhinovirus serogroup A 16 (RV16) for 48 hours. RNA isolated from both AECs and MCs underwent RNA sequencing. Direct effects of epithelial-derived interferons on LAD2 MCs were examined by real-time quantitative PCR.

RESULTS

MCs increased expression of proinflammatory and antiviral genes in AECs. AECs demonstrated a robust antiviral response after RV16 infection that resulted in significant changes in MC gene expression, including upregulation of genes involved in antiviral responses, leukocyte activation, and type 2 inflammation. Subsequent ex vivo modeling demonstrated that IFN-β induces MC type 2 gene expression. The effects of AEC donor phenotype were small relative to the effects of viral infection and the presence of MCs.

CONCLUSIONS

There is significant cross talk between AECs and MCs, which are present in the epithelium in asthma. Epithelial-derived interferons not only play a role in viral suppression but also further alter MC immune responses including specific type 2 genes.

A functional genomics pipeline to identify high-value asthma and allergy CpGs in the human methylome

BACKGROUND

DNA methylation of cytosines at cytosine-phosphate-guanine (CpG) dinucleotides (CpGs) is a widespread epigenetic mark, but genome-wide variation has been relatively unexplored due to the limited representation of variable CpGs on commercial high-throughput arrays.

OBJECTIVES

To explore this hidden portion of the epigenome, this study combined whole-genome bisulfite sequencing with in silico evidence of gene regulatory regions to design a custom array of high-value CpGs. This study focused on airway epithelial cells from children with and without allergic asthma because these cells mediate the effects of inhaled microbes, pollution, and allergens on asthma and allergic disease risk.

METHODS

This study identified differentially methylated regions from whole-genome bisulfite sequencing in nasal epithelial cell DNA from a total of 39 children with and without allergic asthma of both European and African ancestries. This study selected CpGs from differentially methylated regions, previous allergy or asthma epigenome-wide association studies (EWAS), or genome-wide association study loci, and overlapped them with functional annotations for inclusion on a custom Asthma&Allergy array. This study used both the custom and EPIC arrays to perform EWAS of allergic sensitization (AS) in nasal epithelial cell DNA from children in the URECA (Urban Environment and Childhood Asthma) birth cohort and using the custom array in the INSPIRE [Infant Susceptibility to Pulmonary Infections and Asthma Following RSV Exposure] birth cohort. Each CpG on the arrays was assigned to its nearest gene and its promotor capture Hi-C interacting gene and performed expression quantitative trait methylation (eQTM) studies for both sets of genes.

RESULTS

Custom array CpGs were enriched for intermediate methylation levels compared to EPIC CpGs. Intermediate methylation CpGs were further enriched among those associated with AS and for eQTMs on both arrays.

CONCLUSIONS

This study revealed signature features of high-value CpGs and evidence for epigenetic regulation of genes at AS EWAS loci that are robust to race/ethnicity, ascertainment, age, and geography.

Kimma: flexible linear mixed effects modeling with kinship covariance for RNA-seq data

MOTIVATION

The identification of differentially expressed genes (DEGs) from transcriptomic datasets is a major avenue of research across diverse disciplines. However, current bioinformatic tools do not support covariance matrices in DEG modeling. Here, we introduce kimma (Kinship In Mixed Model Analysis), an open-source R package for flexible linear mixed effects modeling including covariates, weights, random effects, covariance matrices, and fit metrics.

RESULTS

In simulated datasets, kimma detects DEGs with similar specificity, sensitivity, and computational time as limma unpaired and dream paired models. Unlike other software, kimma supports covariance matrices as well as fit metrics like Akaike information criterion (AIC). Utilizing genetic kinship covariance, kimma revealed that kinship impacts model fit and DEG detection in a related cohort. Thus, kimma equals or outcompetes current DEG pipelines in sensitivity, computational time, and model complexity.

AVAILABILITY AND IMPLEMENTATION

Kimma is freely available on GitHub https://github.com/BIGslu/kimma with an instructional vignette at https://bigslu.github.io/kimma_vignette/kimma_vignette.html.

Associations between outdoor air pollutants and non-viral asthma exacerbations and airway inflammatory responses in children and adolescents living in urban areas in the USA: a retrospective secondary analysis

BACKGROUND

Asthma prevalence and severity have markedly increased with urbanisation, and children in low-income urban centres have among the greatest asthma morbidity. Outdoor air pollution has been associated with adverse respiratory effects in children with asthma. However, the mechanisms by which air pollution exposure exacerbates asthma, and how these mechanisms compare with exacerbations induced by respiratory viruses, are poorly understood. We aimed to investigate the associations between regional air pollutant concentrations, respiratory illnesses, lung function, and upper airway transcriptional signatures in children with asthma, with particular focus on asthma exacerbations occurring in the absence of respiratory virus.

METHODS

We performed a retrospective analysis of data from the MUPPITS1 cohort and validated our findings in the ICATA cohort. The MUPPITS1 cohort recruited 208 children aged 6-17 years living in urban areas across nine US cities with exacerbation-prone asthma between Oct 7, 2015, and Oct 18, 2016, and monitored them during reported respiratory illnesses. The last MUPPITS1 study visit occurred on Jan 6, 2017. The ICATA cohort recruited 419 participants aged 6-20 years with persistent allergic asthma living in urban sites across eight US cities between Oct 23, 2006, and March 25, 2008, and the last study visit occurred on Dec 30, 2009. We included participants from the MUPPITS1 cohort who reported a respiratory illness at some point during the follow-up and participants from the ICATA cohort who had nasal samples collected during respiratory illness or at a scheduled visit. We used air quality index values and air pollutant concentrations for PM2·5, PM10, O3, NO2, SO2, CO, and Pb from the US Environmental Protection Agency spanning the years of both cohorts, and matched values and concentrations to each illness for each participant. We investigated the associations between regional air pollutant concentrations and respiratory illnesses and asthma exacerbations, pulmonary function, and upper airway transcriptional signatures by use of a combination of generalised additive models, case crossover analyses, and generalised linear mixed-effects models.

FINDINGS

Of the 208 participants from the MUPPITS1 cohort and 419 participants from the ICATA cohort, 168 participants in the MUPPITS1 cohort (98 male participants and 70 female participants) and 189 participants in the ICATA cohort (115 male participants and 74 female participants) were included in our analysis. We identified that increased air quality index values, driven predominantly by increased PM2·5 and O3 concentrations, were significantly associated with asthma exacerbations and decreases in pulmonary function that occurred in the absence of a provoking viral infection. Moreover, individual pollutants were significantly associated with altered gene expression in coordinated inflammatory pathways, including PM2·5 with increased epithelial induction of tissue kallikreins, mucus hypersecretion, and barrier functions and O3 with increased type-2 inflammation.

INTERPRETATION

Our findings suggest that air pollution is an important independent risk factor for asthma exacerbations in children living in urban areas and is potentially linked to exacerbations through specific inflammatory pathways in the airway. Further investigation of these potential mechanistic pathways could inform asthma prevention and management approaches.

FUNDING

National Institutes of Health, National Institute of Allergy and Infectious Diseases.

Effects of combination treatment with tezepelumab and allergen immunotherapy on nasal responses to allergen: A randomized controlled trial

BACKGROUND

Thymic stromal lymphopoietin (TSLP) has been shown to play a central role in the initiation and persistence of allergic responses.

OBJECTIVE

We evaluated whether tezepelumab, a human monoclonal anti-TSLP antibody, improved the efficacy of subcutaneous allergen immunotherapy (SCIT) and promoted the development of tolerance in patients with allergic rhinitis.

METHODS

We conducted a double-blind parallel design trial in patients with cat allergy. A total of 121 patients were randomized to receive either intravenous tezepelumab plus subcutaneous cat SCIT, cat SCIT alone, tezepelumab alone, or placebo for 52 weeks, followed by 52 weeks of observation. Nasal allergen challenge (NAC), skin testing, and blood and nasal samples were obtained throughout the study.

RESULTS

At week 52, the NAC-induced total nasal symptom scores (TNSS) (calculated as area under the curve [AUC0-1h] and as peak score [Peak0-1h] during the first hour after NAC) were significantly reduced in patients receiving tezepelumab/SCIT compared to SCIT alone. At week 104, one year after stopping treatment, the primary end point TNSS AUC0-1h was not significantly different in the tezepelumab/SCIT group compared to SCIT alone, while TNSS Peak0-1h was significantly lower in those receiving combination treatment versus SCIT. Transcriptomic analysis of nasal epithelial samples demonstrated that treatment with the combination of SCIT/tezepelumab, but neither monotherapy, caused persistent downregulation of a gene network related to type 2 inflammation that was associated with improvement in NAC responses.

CONCLUSIONS

Inhibition of TSLP augments the efficacy of SCIT during therapy and may promote tolerance after a 1-year course of treatment. (ClinicalTrials.gov NCT02237196).

Multi-omic association study identifies DNA methylation-mediated genotype and smoking exposure effects on lung function in children living in urban settings

Impaired lung function in early life is associated with the subsequent development of chronic respiratory disease. Most genetic associations with lung function have been identified in adults of European descent and therefore may not represent those most relevant to pediatric populations and populations of different ancestries. In this study, we performed genome-wide association analyses of lung function in a multiethnic cohort of children (n = 1,035) living in low-income urban neighborhoods. We identified one novel locus at the TDRD9 gene in chromosome 14q32.33 associated with percent predicted forced expiratory volume in one second (FEV1) (p = 2.4x10-9; βz = -0.31, 95% CI = -0.41- -0.21). Mendelian randomization and mediation analyses revealed that this genetic effect on FEV1 was partially mediated by DNA methylation levels at this locus in airway epithelial cells, which were also associated with environmental tobacco smoke exposure (p = 0.015). Promoter-enhancer interactions in airway epithelial cells revealed chromatin interaction loops between FEV1-associated variants in TDRD9 and the promoter region of the PPP1R13B gene, a stimulator of p53-mediated apoptosis. Expression of PPP1R13B in airway epithelial cells was significantly associated the FEV1 risk alleles (p = 1.3x10-5; β = 0.12, 95% CI = 0.06-0.17). These combined results highlight a potential novel mechanism for reduced lung function in urban youth resulting from both genetics and smoking exposure.

African-specific alleles modify risk for asthma at the 17q12-q21 locus in African Americans

BACKGROUND

Asthma is the most common chronic disease in children, occurring at higher frequencies and with more severe disease in children with African ancestry.

METHODS

We tested for association with haplotypes at the most replicated and significant childhood-onset asthma locus at 17q12-q21 and asthma in European American and African American children. Following this, we used whole-genome sequencing data from 1060 African American and 100 European American individuals to identify novel variants on a high-risk African American-specific haplotype. We characterized these variants in silico using gene expression and ATAC-seq data from airway epithelial cells, functional annotations from ENCODE, and promoter capture (pc)Hi-C maps in airway epithelial cells. Candidate causal variants were then assessed for correlation with asthma-associated phenotypes in African American children and adults.

RESULTS

Our studies revealed nine novel African-specific common variants, enriched on a high-risk asthma haplotype, which regulated the expression of GSDMA in airway epithelial cells and were associated with features of severe asthma. Using ENCODE annotations, ATAC-seq, and pcHi-C, we narrowed the associations to two candidate causal variants that are associated with features of T2 low severe asthma.

CONCLUSIONS

Previously unknown genetic variation at the 17q12-21 childhood-onset asthma locus contributes to asthma severity in individuals with African ancestries. We suggest that many other population-specific variants that have not been discovered in GWAS contribute to the genetic risk for asthma and other common diseases.

Phenotypes of disease severity in a cohort of hospitalized COVID-19 patients: Results from the IMPACC study

BACKGROUND

Better understanding of the association between characteristics of patients hospitalized with coronavirus disease 2019 (COVID-19) and outcome is needed to further improve upon patient management.

METHODS

Immunophenotyping Assessment in a COVID-19 Cohort (IMPACC) is a prospective, observational study of 1164 patients from 20 hospitals across the United States. Disease severity was assessed using a 7-point ordinal scale based on degree of respiratory illness. Patients were prospectively surveyed for 1 year after discharge for post-acute sequalae of COVID-19 (PASC) through quarterly surveys. Demographics, comorbidities, radiographic findings, clinical laboratory values, SARS-CoV-2 PCR and serology were captured over a 28-day period. Multivariable logistic regression was performed.

FINDINGS

The median age was 59 years (interquartile range [IQR] 20); 711 (61%) were men; overall mortality was 14%, and 228 (20%) required invasive mechanical ventilation. Unsupervised clustering of ordinal score over time revealed distinct disease course trajectories. Risk factors associated with prolonged hospitalization or death by day 28 included age ≥ 65 years (odds ratio [OR], 2.01; 95% CI 1.28-3.17), Hispanic ethnicity (OR, 1.71; 95% CI 1.13-2.57), elevated baseline creatinine (OR 2.80; 95% CI 1.63- 4.80) or troponin (OR 1.89; 95% 1.03-3.47), baseline lymphopenia (OR 2.19; 95% CI 1.61-2.97), presence of infiltrate by chest imaging (OR 3.16; 95% CI 1.96-5.10), and high SARS-CoV2 viral load (OR 1.53; 95% CI 1.17-2.00). Fatal cases had the lowest ratio of SARS-CoV-2 antibody to viral load levels compared to other trajectories over time (p=0.001). 589 survivors (51%) completed at least one survey at follow-up with 305 (52%) having at least one symptom consistent with PASC, most commonly dyspnea (56% among symptomatic patients). Female sex was the only associated risk factor for PASC.

INTERPRETATION

Integration of PCR cycle threshold, and antibody values with demographics, comorbidities, and laboratory/radiographic findings identified risk factors for 28-day outcome severity, though only female sex was associated with PASC. Longitudinal clinical phenotyping offers important insights, and provides a framework for immunophenotyping for acute and long COVID-19.

FUNDING

NIH.

Genome-wide study of early and severe childhood asthma identifies interaction between CDHR3 and GSDMB

BACKGROUND

Asthma with severe exacerbation is one of the most common causes of hospitalization among young children. Exacerbations are typically triggered by respiratory infections, but the host factors causing recurrent infections and exacerbations in some children are poorly understood. As a result, current treatment options and preventive measures are inadequate.

OBJECTIVE

We sought to identify genetic interaction associated with the development of childhood asthma.

METHODS

We performed an exhaustive search for pairwise interaction between genetic single nucleotide polymorphisms using 1204 cases of a specific phenotype of early childhood asthma with severe exacerbations in patients aged 2 to 6 years combined with 5328 nonasthmatic controls. Replication was attempted in 3 independent populations, and potential underlying immune mechanisms were investigated in the COPSAC2010 and COPSAC2000 birth cohorts.

RESULTS

We found evidence of interaction, including replication in independent populations, between the known childhood asthma loci CDHR3 and GSDMB. The effect of CDHR3 was dependent on the GSDMB genotype, and this interaction was more pronounced for severe and early onset of disease. Blood immune analyses suggested a mechanism related to increased IL-17A production after viral stimulation.

CONCLUSIONS

We found evidence of interaction between CDHR3 and GSDMB in development of early childhood asthma, possibly related to increased IL-17A response to viral infections. This study demonstrates the importance of focusing on specific disease subtypes for understanding the genetic mechanisms of asthma.

Mepolizumab for urban children with exacerbation-prone eosinophilic asthma in the USA (MUPPITS-2): a randomised, double-blind, placebo-controlled, parallel-group trial

BACKGROUND

Black and Hispanic children living in urban environments in the USA have an excess burden of morbidity and mortality from asthma. Therapies directed at the eosinophilic phenotype reduce asthma exacerbations in adults, but few data are available in children and diverse populations. Furthermore, the molecular mechanisms that underlie exacerbations either being prevented by, or persisting despite, immune-based therapies are not well understood. We aimed to determine whether mepolizumab, added to guidelines-based care, reduced the number of asthma exacerbations during a 52-week period compared with guidelines-based care alone.

METHODS

This is a randomised, double-blind, placebo-controlled, parallel-group trial done at nine urban medical centres in the USA. Children and adolescents aged 6-17 years, who lived in socioeconomically disadvantaged neighbourhoods and had exacerbation-prone asthma (defined as ≥two exacerbations in the previous year) and blood eosinophils of at least 150 cells per μL were randomly assigned 1:1 to mepolizumab (6-11 years: 40 mg; 12-17 years: 100 mg) or placebo injections once every 4 weeks, plus guideline-based care, for 52 weeks. Randomisation was done using a validated automated system. Participants, investigators, and the research staff who collected outcome measures remained masked to group assignments. The primary outcome was the number of asthma exacerbations that were treated with systemic corticosteroids during 52 weeks in the intention-to-treat population. The mechanisms of treatment response were assessed by study investigators using nasal transcriptomic modular analysis. Safety was assessed in the intention-to-treat population. This trial is registered with ClinicalTrials.gov, NCT03292588.

FINDINGS

Between Nov 1, 2017, and Mar 12, 2020, we recruited 585 children and adolescents. We screened 390 individuals, of whom 335 met the inclusion criteria and were enrolled. 290 met the randomisation criteria, were randomly assigned to mepolizumab (n=146) or placebo (n=144), and were included in the intention-to-treat analysis. 248 completed the study. The mean number of asthma exacerbations within the 52-week study period was 0·96 (95% CI 0·78-1·17) with mepolizumab and 1·30 (1·08-1·57) with placebo (rate ratio 0·73; 0·56-0·96; p=0·027). Treatment-emergent adverse events occurred in 42 (29%) of 146 participants in the mepolizumab group versus 16 (11%) of 144 participants in the placebo group. No deaths were attributed to mepolizumab.

INTERPRETATION

Phenotype-directed therapy with mepolizumab in urban children with exacerbation-prone eosinophilic asthma reduced the number of exacerbations.

FUNDING

US National Institute of Allergy and Infectious Diseases and GlaxoSmithKline.

Seasonal airway microbiome and transcriptome interactions promote childhood asthma exacerbations

BACKGROUND

Seasonal variation in respiratory illnesses and exacerbations in pediatric populations with asthma is well described, though whether upper airway microbes play season-specific roles in these events is unknown.

OBJECTIVE

We hypothesized that nasal microbiota composition is seasonally dynamic and that discrete microbe-host interactions modify risk of asthma exacerbation in a season-specific manner.

METHODS

Repeated nasal samples from children with exacerbation-prone asthma collected during periods of respiratory health (baseline; n = 181 samples) or first captured respiratory illness (n = 97) across all seasons, underwent bacterial (16S ribosomal RNA gene) and fungal (internal transcribed spacer region 2) biomarker sequencing. Virus detection was performed by multiplex PCR. Paired nasal transcriptome data were examined for seasonal dynamics and integrative analyses.

RESULTS

Upper airway bacterial and fungal microbiota and rhinovirus detection exhibited significant seasonal dynamics. In seasonally adjusted analysis, variation in both baseline and respiratory illness microbiota related to subsequent exacerbation. Specifically, in the fall, when respiratory illness and exacerbation events were most frequent, several Moraxella and Haemophilus members were enriched both in virus-positive respiratory illnesses and those that progressed to exacerbations. The abundance of 2 discrete bacterial networks, characteristically comprising either Streptococcus or Staphylococcus, exhibited opposing interactions with an exacerbation-associated SMAD3 nasal epithelial transcriptional module to significantly increase the odds of subsequent exacerbation (odds ratio = 14.7, 95% confidence interval = 1.50-144, P = .02; odds ratio = 39.17, 95% confidence interval = 2.44-626, P = .008, respectively).

CONCLUSIONS

Upper airway microbiomes covary with season and with seasonal trends in respiratory illnesses and asthma exacerbations. Seasonally adjusted analyses reveal specific bacteria-host interactions that significantly increase risk of asthma exacerbation in these children.

Fine-mapping studies distinguish genetic risks for childhood- and adult-onset asthma in the HLA region

BACKGROUND

Genome-wide association studies of asthma have revealed robust associations with variation across the human leukocyte antigen (HLA) complex with independent associations in the HLA class I and class II regions for both childhood-onset asthma (COA) and adult-onset asthma (AOA). However, the specific variants and genes contributing to risk are unknown.

METHODS

We used Bayesian approaches to perform genetic fine-mapping for COA and AOA (n=9432 and 21,556, respectively; n=318,167 shared controls) in White British individuals from the UK Biobank and to perform expression quantitative trait locus (eQTL) fine-mapping in immune (lymphoblastoid cell lines, n=398; peripheral blood mononuclear cells, n=132) and airway (nasal epithelial cells, n=188) cells from ethnically diverse individuals. We also examined putatively causal protein coding variation from protein crystal structures and conducted replication studies in independent multi-ethnic cohorts from the UK Biobank (COA n=1686; AOA n=3666; controls n=56,063).

RESULTS

Genetic fine-mapping revealed both shared and distinct causal variation between COA and AOA in the class I region but only distinct causal variation in the class II region. Both gene expression levels and amino acid variation contributed to risk. Our results from eQTL fine-mapping and amino acid visualization suggested that the HLA-DQA1*03:01 allele and variation associated with expression of the nonclassical HLA-DQA2 and HLA-DQB2 genes accounted entirely for the most significant association with AOA in GWAS. Our studies also suggested a potentially prominent role for HLA-C protein coding variation in the class I region in COA. We replicated putatively causal variant associations in a multi-ethnic cohort.

CONCLUSIONS

We highlight roles for both gene expression and protein coding variation in asthma risk and identified putatively causal variation and genes in the HLA region. A convergence of genomic, transcriptional, and protein coding evidence implicates the HLA-DQA2 and HLA-DQB2 genes and HLA-DQA1*03:01 allele in AOA.

Location of eosinophils in the airway wall is critical for specific features of airway hyperresponsiveness and T2 inflammation in asthma

Eosinophils are implicated as effector cells in asthma but the functional implications of the precise location of eosinophils in the airway wall is poorly understood. We aimed to quantify eosinophils in the different compartments of the airway wall and associate these findings with clinical features of asthma and markers of airway inflammation.In this cross-sectional study, we utilised design-based stereology to accurately partition the numerical density of eosinophils in both the epithelial compartment and the subepithelial space (airway wall area below the basal lamina including the submucosa) in individuals with and without asthma and related these findings to airway hyperresponsiveness (AHR) and features of airway inflammation.Intraepithelial eosinophils were linked to the presence of asthma and endogenous AHR, the type of AHR that is most specific for asthma. In contrast, both intraepithelial and subepithelial eosinophils were associated with type-2 (T2) inflammation, with the strongest association between IL5 expression and intraepithelial eosinophils. Eosinophil infiltration of the airway wall was linked to a specific mast cell phenotype that has been described in asthma. We found that IL-33 and IL-5 additively increased cysteinyl leukotriene (CysLT) production by eosinophils and that the CysLT LTC₄ along with IL-33 increased IL13 expression in mast cells and altered their protease profile.We conclude that intraepithelial eosinophils are associated with endogenous AHR and T2 inflammation and may interact with intraepithelial mast cells via CysLTs to regulate airway inflammation.

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.

Heterogeneity of magnitude, allergen immunodominance, and cytokine polarization of cockroach allergen-specific T cell responses in allergic sensitized children

BACKGROUND

Characterization of allergic responses to cockroach (CR), a common aeroallergen associated with asthma, has focused mainly on IgE reactivity, but little is known about T cell responses, particularly in children. We conducted a functional evaluation of CR allergen-specific T cell reactivity in a cohort of CR allergic children with asthma.

METHODS

Peripheral blood mononuclear cells (PBMCs) were obtained from 71 children, with mild-to-moderate asthma who were enrolled in a CR immunotherapy (IT) clinical trial, prior to treatment initiation. PBMC were stimulated with peptide pools derived from 11 CR allergens, and CD4+ T cell responses assessed by intracellular cytokine staining.

RESULTS

Highly heterogeneous responses in T cell reactivity were observed among participants, both in terms of the magnitude of cytokine response and allergen immunodominance. Reactivity against Bla g 9 and Bla g 5 was most frequent. The phenotype of the T cell response was dominated by IL-4 production and a Th2 polarized profile in 54.9% of participants, but IFNγ production and Th1 polarization was observed in 25.3% of the participants. The numbers of regulatory CD4+ T cells were also highly variable and the magnitude of effector responses and Th2 polarization were positively correlated with serum IgE levels specific to a clinical CR extract.

CONCLUSIONS

Our results demonstrate that in children with mild-to-moderate asthma, CR-specific T cell responses display a wide range of magnitude, allergen dominance, and polarization. These results will enable examination of whether any of the variables measured are affected by IT and/or are predictive of clinical outcomes.

Management Strategies to Reduce Exacerbations in non-T2 Asthma

There have been considerable advances in our understanding of asthmatic airway inflammation, resulting in a paradigm shift of classifying individuals on the basis of either the presence or the absence of type 2 (T2) inflammatory markers. Several novel monoclonal antibody therapies targeting T2 cytokines have demonstrated significant clinical effects including reductions in acute exacerbations and improvements in asthma-related quality of life and lung function for individuals with T2-high asthma. However, there have been fewer advancements in developing therapies for those without evidence of T2 airway inflammation (so-called non-T2 asthma). Here, we review the heterogeneity of molecular mechanisms responsible for initiation and regulation of non-T2 inflammation and discuss both current and potential future therapeutic options for individuals with non-T2 asthma.

Development of a fixed module repertoire for the analysis and interpretation of blood transcriptome data

As the capacity for generating large-scale molecular profiling data continues to grow, the ability to extract meaningful biological knowledge from it remains a limitation. Here, we describe the development of a new fixed repertoire of transcriptional modules, BloodGen3, that is designed to serve as a stable reusable framework for the analysis and interpretation of blood transcriptome data. The construction of this repertoire is based on co-clustering patterns observed across sixteen immunological and physiological states encompassing 985 blood transcriptome profiles. Interpretation is supported by customized resources, including module-level analysis workflows, fingerprint grid plot visualizations, interactive web applications and an extensive annotation framework comprising functional profiling reports and reference transcriptional profiles. Taken together, this well-characterized and well-supported transcriptional module repertoire can be employed for the interpretation and benchmarking of blood transcriptome profiles within and across patient cohorts. Blood transcriptome fingerprints for the 16 reference cohorts can be accessed interactively via:  https://drinchai.shinyapps.io/BloodGen3Module/ .

Inducible expression quantitative trait locus analysis of the MUC5AC gene in asthma in urban populations of children

BACKGROUND

Mucus plugging can worsen asthma control, lead to reduced lung function and fatal exacerbations. MUC5AC is the secretory mucin implicated in mucus plugging, and MUC5AC gene expression has been associated with development of airway obstruction and asthma exacerbations in urban children with asthma. However, the genetic determinants of MUC5AC expression are not established.

OBJECTIVES

This study sought to assess single-nucleotide polymorphisms (SNPs) that influence MUC5AC expression and relate to pulmonary functions in childhood asthma.

METHODS

This study used RNA-sequencing data from upper airway samples and performed cis-expression quantitative trait loci (eQTL) and allele-specific expression analyses in 2 cohorts of predominantly Black and Hispanic urban children, a high asthma-risk birth cohort, and an exacerbation-prone asthma cohort. Inducible MUC5AC eQTLs were further investigated during incipient asthma exacerbations. Significant eQTLs SNPs were tested for associations with lung function measurements and their functional consequences were investigated in DNA regulatory databases.

RESULTS

Two independent groups of SNPs in the MUC5AC gene that were significantly associated with MUC5AC expression were identified. Moreover, these SNPs showed stronger eQTL associations with MUC5AC expression during asthma exacerbations, which is consistent with inducible expression. SNPs in 1 group also showed significant association with decreased pulmonary functions. These SNPs included multiple EGR1 transcription factor binding sites, suggesting a mechanism of effect.

CONCLUSIONS

These findings demonstrate the applicability of organ-specific RNA-sequencing data to determine genetic factors contributing to a key disease pathway. Specifically, they suggest important genetic variations that may underlie propensity to mucus plugging in asthma and could be important in targeted asthma phenotyping and disease management strategies.

A consequentialist argument for considering age in triage decisions during the coronavirus pandemic

Most ethics guidelines for distributing scarce medical resources during the coronavirus pandemic seek to save the most lives and the most life-years. A patient's prognosis is determined using a SOFA or MSOFA score to measure likelihood of survival to discharge, as well as a consideration of relevant comorbidities and their effects on likelihood of survival up to one or five years. Although some guidelines use age as a tiebreaker when two patients' prognoses are identical, others refuse to consider age for fear of discriminating against the elderly. In this paper, I argue that age is directly relevant for maximizing health benefits, so current ethics guidelines are wrongly excluding or deemphasizing life-stage in their triage algorithms. Research on COVID-19 has shown that age is a risk factor in adverse outcomes, independent of comorbidities. And limiting a consideration of life-years to only one or five years past discharge does not maximize health benefits. Therefore, based on their own stated values, triage algorithms for coronavirus patients ought to include life-stage as a primary consideration, along with the SOFA score and comorbidities, rather than excluding it or using it merely as a tiebreaker. This is not discriminatory because patients ought to have equal opportunity to experience life-stages. The equitable enforcement of that right justifies unequal treatment based on age in cases when there is a scarcity of life-saving resources. A consideration of life-stage would thus allow healthcare workers to responsibly steward public resources in order to maximize lives and life-years saved.

Endotype of allergic asthma with airway obstruction in urban children

BACKGROUND

Black and Hispanic children growing up in disadvantaged urban neighborhoods have the highest rates of asthma and related morbidity in the United States.

OBJECTIVES

This study sought to identify specific respiratory phenotypes of health and disease in this population, associations with early life exposures, and molecular patterns of gene expression in nasal epithelial cells that underlie clinical disease.

METHODS

The study population consisted of 442 high-risk urban children who had repeated assessments of wheezing, allergen-specific IgE, and lung function through 10 years of age. Phenotypes were identified by developing temporal trajectories for these data, and then compared to early life exposures and patterns of nasal epithelial gene expression at 11 years of age.

RESULTS

Of the 6 identified respiratory phenotypes, a high wheeze, high atopy, low lung function group had the greatest respiratory morbidity. In early life, this group had low exposure to common allergens and high exposure to ergosterol in house dust. While all high-atopy groups were associated with increased expression of a type-2 inflammation gene module in nasal epithelial samples, an epithelium IL-13 response module tracked closely with impaired lung function, and a MUC5AC hypersecretion module was uniquely upregulated in the high wheeze, high atopy, low lung function group. In contrast, a medium wheeze, low atopy group showed altered expression of modules of epithelial integrity, epithelial injury, and antioxidant pathways.

CONCLUSIONS

In the first decade of life, high-risk urban children develop distinct phenotypes of respiratory health versus disease that link early life environmental exposures to childhood allergic sensitization and asthma. Moreover, unique patterns of airway gene expression demonstrate how specific molecular pathways underlie distinct respiratory phenotypes, including allergic and nonallergic asthma.

BloodGen3Module: Blood transcriptional module repertoire analysis and visualization using R

Motivation

We previously described the construction and characterization of fixed reusable blood transcriptional module repertoires. More recently we released a third iteration (‘BloodGen3’ module repertoire) that comprises 382 functionally annotated modules and encompasses 14 168 transcripts. Custom bioinformatic tools are needed to support downstream analysis, visualization and interpretation relying on such fixed module repertoires.

Results

We have developed and describe here an R package, BloodGen3Module. The functions of our package permit group comparison analyses to be performed at the module-level, and to display the results as annotated fingerprint grid plots. A parallel workflow for computing module repertoire changes for individual samples rather than groups of samples is also available; these results are displayed as fingerprint heatmaps. An illustrative case is used to demonstrate the steps involved in generating blood transcriptome repertoire fingerprints of septic patients. Taken together, this resource could facilitate the analysis and interpretation of changes in blood transcript abundance observed across a wide range of pathological and physiological states.

Availability and implementation

The BloodGen3Module package and documentation are freely available from Github: https://github.com/Drinchai/BloodGen3Module.

Supplementary information

Supplementary data are available at Bioinformatics online.

Definition of erythroid cell-positive blood transcriptome phenotypes associated with severe respiratory syncytial virus infection

Biomarkers to assess the risk of developing severe respiratory syncytial virus (RSV) infection are needed. We conducted a meta-analysis of 490 unique profiles from six public RSV blood transcriptome datasets. A repertoire of 382 well-characterized transcriptional modules was used to define dominant host responses to RSV infection. The consolidated RSV cohort was stratified according to four traits: "interferon response" (IFN), "neutrophil-driven inflammation" (Infl), "cell cycle" (CC), and "erythrocytes" (Ery). We identified eight prevalent blood transcriptome phenotypes, of which three Ery+ phenotypes comprised higher proportions of patients requiring intensive care. This finding confirms on a larger scale data from one of our earlier reports describing an association between an erythrocyte signature and RSV disease severity. Further contextual interpretation made it possible to associate this signature with immunosuppressive states (late stage cancer, pharmacological immunosuppression), and with a population of fetal glycophorin A+ erythroid precursors. Furthermore, we posit that this erythrocyte cell signature may be linked to a population of immunosuppressive erythroid cells previously described in the literature, and that overabundance of this cell population in RSV patients may underlie progression to severe disease. These findings outline potential priority areas for biomarker development and investigations into the immune biology of RSV infection. The approach that we developed and employed here should also permit to delineate prevalent blood transcriptome phenotypes in other settings.

A modular framework for the development of targeted Covid-19 blood transcript profiling panels

BACKGROUND

Covid-19 morbidity and mortality are associated with a dysregulated immune response. Tools are needed to enhance existing immune profiling capabilities in affected patients. Here we aimed to develop an approach to support the design of targeted blood transcriptome panels for profiling the immune response to SARS-CoV-2 infection.

METHODS

We designed a pool of candidates based on a pre-existing and well-characterized repertoire of blood transcriptional modules. Available Covid-19 blood transcriptome data was also used to guide this process. Further selection steps relied on expert curation. Additionally, we developed several custom web applications to support the evaluation of candidates.

RESULTS

As a proof of principle, we designed three targeted blood transcript panels, each with a different translational connotation: immunological relevance, therapeutic development relevance and SARS biology relevance.

CONCLUSION

Altogether the work presented here may contribute to the future expansion of immune profiling capabilities via targeted profiling of blood transcript abundance in Covid-19 patients.

Evolving concepts in how viruses impact asthma: A Work Group Report of the Microbes in Allergy Committee of the American Academy of Allergy, Asthma & Immunology

Over the past decade, there have been substantial advances in our understanding about how viral infections regulate asthma. Important lessons have been learned from birth cohort studies examining viral infections and subsequent asthma and from understanding the relationships between host genetics and viral infections, the contributions of respiratory viral infections to patterns of immune development, the impact of environmental exposure on the severity of viral infections, and how the viral genome influences host immune responses to viral infections. Further, there has been major progress in our knowledge about how bacteria regulate host immune responses in asthma pathogenesis. In this article, we also examine the dynamics of bacterial colonization of the respiratory tract during viral upper respiratory tract infection, in addition to the relationship of the gut and respiratory microbiomes with respiratory viral infections. Finally, we focus on potential interventions that could decrease virus-induced wheezing and asthma. There are emerging therapeutic options to decrease the severity of wheezing exacerbations caused by respiratory viral infections. Primary prevention is a major goal, and a strategy toward this end is considered.

Nasal allergen challenge and environmental exposure chamber challenge: A randomized trial comparing clinical and biological responses to cat allergen

BACKGROUND

The direct-instillation nasal allergen challenge (NAC) and the environmental exposure chamber (EEC) are 2 methods of conducting controlled allergen provocations. The clinical and biological comparability of these methods has not been thoroughly investigated.

OBJECTIVE

We sought to compare clinical and immunologic responses to cat allergen in NAC versus EEC.

METHODS

Twenty-four participants were randomized to receive either NAC followed by a 2-day challenge in an EEC or a 2-day challenge in an EEC followed by NAC. Challenges were separated by 28-day washout periods. We measured total nasal symptom scores, peak nasal inspiratory flow, nasal (0-8 hours) and serum cytokines, serum antibodies, peripheral blood antigen-specific T lymphocytes, and gene expression in nasal scrapings. The primary outcome was the total nasal symptom score area under the curve for the first 3 hours after allergen exposure in NAC or after initiation of exposure in EEC.

RESULTS

Both challenges increased IL-5 and IL-13 in nasal fluids and serum and resulted in altered nasal cell expression of gene modules related to mucosal biology and transcriptional regulation. Changes in gene modules, more so than cytokine measurements, showed significant associations with total nasal symptom score and peak nasal inspiratory flow. Overall, EEC exposure generated larger responses and more early terminations compared with NAC. Although the 2 challenges did not correlate in symptom magnitude or temporality, striking correlations were observed in cytokine levels.

CONCLUSIONS

Although clinical outcomes of NAC and EEC were temporally different and nonequivalent in magnitude, immunologic responses were similar. Selection of a particular allergen challenge method should depend on considerations of study objectives and cost.

A 71-year-old man with recurrent pulmonary mycobacterial avium complex infections and lymphopenia

Mycobacterium avium complex (MAC) infections, generally viewed as opportunistic infections, often trigger an evaluation for an underlying immunodeficiency disorder. However, MAC infections can occur in patients who presumably are immunocompetent, particularly in those with an underlying structural lung disease. T-cell immunity plays a critical role in controlling MAC infection. We presented a case of lymphopenia, which complicated the clinical course of a pulmonary MAC infection in a patient who was negative for human immunodeficiency virus.

Airway epithelium-shifted mast cell infiltration regulates asthmatic inflammation via IL-33 signaling

Asthma is a heterogeneous syndrome that has been subdivided into physiologic phenotypes and molecular endotypes. The most specific phenotypic manifestation of asthma is indirect airway hyperresponsiveness (AHR), and a prominent molecular endotype is the presence of type 2 inflammation. The underlying basis for type 2 inflammation and its relationship to AHR are incompletely understood. We assessed the expression of type 2 cytokines in the airways of subjects with and without asthma who were extensively characterized for AHR. Using quantitative morphometry of the airway wall, we identified a shift in mast cells from the submucosa to the airway epithelium specifically associated with both type 2 inflammation and indirect AHR. Using ex vivo modeling of primary airway epithelial cells in organotypic coculture with mast cells, we show that epithelial-derived IL-33 uniquely induced type 2 cytokines in mast cells, which regulated the expression of epithelial IL33 in a feed-forward loop. This feed-forward loop was accentuated in epithelial cells derived from subjects with asthma. These results demonstrate that type 2 inflammation and indirect AHR in asthma are related to a shift in mast cell infiltration to the airway epithelium, and that mast cells cooperate with epithelial cells through IL-33 signaling to regulate type 2 inflammation.

Transcriptome networks identify mechanisms of viral and nonviral asthma exacerbations in children

Respiratory infections are common precursors to asthma exacerbations in children, but molecular immune responses that determine whether and how an infection causes an exacerbation are poorly understood. By using systems-scale network analysis, we identify repertoires of cellular transcriptional pathways that lead to and underlie distinct patterns of asthma exacerbation. Specifically, in both virus-associated and nonviral exacerbations, we demonstrate a set of core exacerbation modules, among which epithelial-associated SMAD3 signaling is upregulated and lymphocyte response pathways are downregulated early in exacerbation, followed by later upregulation of effector pathways including epidermal growth factor receptor signaling, extracellular matrix production, mucus hypersecretion, and eosinophil activation. We show an additional set of multiple inflammatory cell pathways involved in virus-associated exacerbations, in contrast to squamous cell pathways associated with nonviral exacerbations. Our work introduces an in vivo molecular platform to investigate, in a clinical setting, both the mechanisms of disease pathogenesis and therapeutic targets to modify exacerbations.