Nasal epithelial gene expression identifies relevant asthma endotypes in the ATLANTIS study.

INTRODUCTION: Asthma is an inflammatory airways disease encompassing multiple phenotypes and endotypes. Several studies suggested gene expression in nasal epithelium to serve as a proxy for bronchial epithelium, being a non-invasive approach to investigate lung diseases. We hypothesised that molecular differences in upper airway epithelium reflect asthma-associated differences in the lower airways and are associated with clinical expression of asthma. METHODS: We analysed nasal epithelial gene expression data from 369 patients with asthma and 58 non-asthmatic controls from the Assessment of Small Airways Involvement in Asthma study. Unsupervised hierarchical clustering was performed on asthma-associated genes. Asthma-associated gene signatures were replicated in independent cohorts with nasal and bronchial brushes data by comparing Gene Set Variation Analysis scores between asthma patients and non-asthmatic controls. RESULTS: We identified 67 higher expressed and 59 lower expressed genes in nasal epithelium from asthma patients compared with controls (false discovery rate<0.05), including CLCA1, CST1 and POSTN, genes well known to reflect asthma in bronchial airway epithelium. Hierarchical clustering revealed several molecular asthma endotypes with distinct clinical characteristics, including an endotype with higher blood and sputum eosinophils, high fractional exhaled nitric oxide, and more severe small airway dysfunction, as reflected by lower forced expiratory flow at 50%. In an independent cohort, we demonstrated that genes higher expressed in the nasal epithelium reflect asthma-associated changes in the lower airways. CONCLUSION: Our results show that the nasal epithelial gene expression profile reflects asthma-related processes in the lower airways. We suggest that nasal epithelium may be a useful non-invasive tool to identify asthma endotypes and may advance personalised management of the disease.

SARS-CoV-2-infected human airway epithelial cell cultures uniquely lack interferon and immediate early gene responses caused by other coronaviruses.

Objectives: Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a member of a class of highly pathogenic coronaviruses. The large family of coronaviruses, however, also includes members that cause only mild symptoms, like human coronavirus-229E (HCoV-229E) or OC43 (HCoV-OC43). Unravelling how molecular (and cellular) pathophysiology differs between highly and low pathogenic coronaviruses is important for the development of therapeutic strategies. Methods: Here, we analysed the transcriptome of primary human bronchial epithelial cells (PBEC), differentiated at the air-liquid interface (ALI) after infection with SARS-CoV-2, SARS-CoV, Middle East Respiratory Syndrome (MERS)-CoV and HCoV-229E using bulk RNA sequencing. Results: ALI-PBEC were efficiently infected by all viruses, and SARS-CoV, MERS-CoV and HCoV-229E infection resulted in a largely similar transcriptional response. The response to SARS-CoV-2 infection differed markedly as it uniquely lacked the increase in expression of immediate early genes, including FOS, FOSB and NR4A1 that was observed with all other coronaviruses. This finding was further confirmed in publicly available experimental and clinical datasets. Interfering with NR4A1 signalling in Calu-3 lung epithelial cells resulted in a 100-fold reduction in extracellular RNA copies of SARS-CoV-2 and MERS-CoV, suggesting an involvement in virus replication. Furthermore, a lack in induction of interferon-related gene expression characterised the main difference between the highly pathogenic coronaviruses and low pathogenic viruses HCoV-229E and HCoV-OC43. Conclusion: Our results demonstrate a previously unknown suppression of a host response gene set by SARS-CoV-2 and confirm a difference in interferon-related gene expression between highly pathogenic and low pathogenic coronaviruses.

Lung Tissue Multi-Layer Network Analysis Uncovers the Molecular Heterogeneity of COPD.

BACKGROUND: Chronic Obstructive Pulmonary Disease (COPD) is a heterogeneous condition. We hypothesized that the unbiased integration of different COPD lung omics using a novel multi-layer approach may unravel mechanisms associated with clinical characteristics. METHODS: We profiled mRNA, miRNA and methylome in lung tissue samples from 135 former smokers with COPD. For each omic (layer) we built a patient network based on molecular similarity. The three networks were used to build a multi-layer network, and optimization of multiplex-modularity was employed to identify patient communities across the three distinct layers. Uncovered communities were related to clinical features. RESULTS: We identified five patient communities in the multi-layer network which were molecularly distinct and related to clinical characteristics, such as FEV1 and blood eosinophils. Two communities (C#3 and C#4) had both similarly low FEV1 values and emphysema, but were molecularly different: C#3, but not C#4, presented B and T cell signatures and a downregulation of secretory (SCGB1A1/SCGB3A1) and ciliated cells. A machine learning model was set up to discriminate C#3 and C#4 in our cohort, and to validate them in an independent cohort. Finally, using spatial transcriptomics we characterized the small airway differences between C#3 and C#4, identifying an upregulation of T/B cell homing chemokines, and bacterial response genes in C#3. CONCLUSIONS: A novel multi-layer network analysis is able to identify clinically relevant COPD patient communities. Patients with similarly low FEV1 and emphysema can have molecularly distinct small airways and immune response patterns, indicating that different endotypes can lead to similar clinical presentation.

Alternative Splicing Is a Major Factor Shaping Transcriptome Diversity in Mild and Severe Chronic Obstructive Pulmonary Disease.

The role of alternative splicing in chronic obstructive pulmonary disease (COPD) is still largely unknown. We aimed to investigate the differences in alternatively splicing events between patients with mild-to-moderate and severe COPD compared with non-COPD control subjects and to identify splicing factors associated with aberrant alternative splicing in COPD. For this purpose, we performed genome-wide RNA-sequencing analysis of bronchial brushings from 23 patients with mild-to-moderate COPD, 121 with severe COPD, and 23 non-COPD control subjects. We found a significant difference in the frequency of alternative splicing events in patients with mild-to-moderate and severe COPD compared with non-COPD control subjects. There were from two to eight times (depending on event type) more differential alternative splicing events in the severe than in the mild-to-moderate stage. The severe COPD samples showed less intron retention and more exon skipping. It is interesting that the transcript levels of the top 10 differentially expressed splicing factors were significantly correlated with the percentage of many alternatively spliced transcripts in severe COPD. The aberrant alternative splicing in severe COPD was predicted to increase the overall protein-coding capacity of gene products. In conclusion, we observed large and significant differences in alternative splicing between bronchial samples of patients with COPD and control subjects, with more events observed in severe than in mild-to-moderate COPD. The changes in the expression of several splicing factors correlated with prevalence of alternative splicing in severe COPD. Alternative splicing can indirectly impact gene expression by changing the relative abundance of protein-coding isoforms potentially influencing pathophysiological changes. The results provide a better understanding of COPD-related alternative splicing changes.

Faecal microbial transfer and complex carbohydrates mediate protection against COPD.

OBJECTIVE: Chronic obstructive pulmonary disease (COPD) is a major cause of global illness and death, most commonly caused by cigarette smoke. The mechanisms of pathogenesis remain poorly understood, limiting the development of effective therapies. The gastrointestinal microbiome has been implicated in chronic lung diseases via the gut-lung axis, but its role is unclear. DESIGN: Using an in vivo mouse model of cigarette smoke (CS)-induced COPD and faecal microbial transfer (FMT), we characterised the faecal microbiota using metagenomics, proteomics and metabolomics. Findings were correlated with airway and systemic inflammation, lung and gut histopathology and lung function. Complex carbohydrates were assessed in mice using a high resistant starch diet, and in 16 patients with COPD using a randomised, double-blind, placebo-controlled pilot study of inulin supplementation. RESULTS: FMT alleviated hallmark features of COPD (inflammation, alveolar destruction, impaired lung function), gastrointestinal pathology and systemic immune changes. Protective effects were additive to smoking cessation, and transfer of CS-associated microbiota after antibiotic-induced microbiome depletion was sufficient to increase lung inflammation while suppressing colonic immunity in the absence of CS exposure. Disease features correlated with the relative abundance of Muribaculaceae, Desulfovibrionaceae and Lachnospiraceae family members. Proteomics and metabolomics identified downregulation of glucose and starch metabolism in CS-associated microbiota, and supplementation of mice or human patients with complex carbohydrates improved disease outcomes. CONCLUSION: The gut microbiome contributes to COPD pathogenesis and can be targeted therapeutically.

TLR7 promotes smoke-induced experimental lung damage through the activity of mast cell tryptase.

Toll-like receptor 7 (TLR7) is known for eliciting immunity against single-stranded RNA viruses, and is increased in both human and cigarette smoke (CS)-induced, experimental chronic obstructive pulmonary disease (COPD). Here we show that the severity of CS-induced emphysema and COPD is reduced in TLR7-deficient mice, while inhalation of imiquimod, a TLR7-agonist, induces emphysema without CS exposure. This imiquimod-induced emphysema is reduced in mice deficient in mast cell protease-6, or when wild-type mice are treated with the mast cell stabilizer, cromolyn. Furthermore, therapeutic treatment with anti-TLR7 monoclonal antibody suppresses CS-induced emphysema, experimental COPD and accumulation of pulmonary mast cells in mice. Lastly, TLR7 mRNA is increased in pre-existing datasets from patients with COPD, while TLR7+ mast cells are increased in COPD lungs and associated with severity of COPD. Our results thus support roles for TLR7 in mediating emphysema and COPD through mast cell activity, and may implicate TLR7 as a potential therapeutic target.

Gene signatures in U-BIOPRED severe asthma for molecular phenotyping and precision medicine: time for clinical use.

INTRODUCTION: The use and generation of gene signatures have been established as a method to define molecular endotypes in complex diseases such as severe asthma. Bioinformatic approaches have now been applied to large omics datasets to define the various co-existing inflammatory and cellular functional pathways driving or characterizing a particular molecular endotype. AREAS COVERED: Molecular phenotypes and endotypes of Type 2 inflammatory pathways and also of non-Type 2 inflammatory pathways, such as IL-6 trans-signaling, IL-17 activation, and IL-22 activation, have been defined in the Unbiased Biomarkers for the Prediction of Respiratory Disease Outcomes dataset. There has also been the identification of the role of mast cell activation and of macrophage dysfunction in various phenotypes of severe asthma. EXPERT OPINION: Phenotyping on the basis of clinical treatable traits is not sufficient for understanding of mechanisms driving the disease in severe asthma. It is time to consider whether certain patients with severe asthma, such as those non-responsive to current therapies, including Type 2 biologics, would be better served using an approach of molecular endotyping using gene signatures for management purposes rather than the current sole reliance on blood eosinophil counts or exhaled nitric oxide measurements.

A bronchial gene signature specific for severe COPD that is retained in the nose.

Introduction: A subset of COPD patients develops advanced disease with severe airflow obstruction, hyperinflation and extensive emphysema. We propose that the pathogenesis in these patients differs from mild-moderate COPD and is reflected by bronchial gene expression. The aim of the present study was to identify a unique bronchial epithelial gene signature for severe COPD patients. Methods: We obtained RNA sequencing data from bronchial brushes from 123 ex-smokers with severe COPD, 23 with mild-moderate COPD and 23 non-COPD controls. We identified genes specific to severe COPD by comparing severe COPD to non-COPD controls, followed by removing genes that were also differentially expressed between mild-moderate COPD and non-COPD controls. Next, we performed a pathway analysis on these genes and evaluated whether this signature is retained in matched nasal brushings. Results: We identified 219 genes uniquely differentially expressed in severe COPD. Interaction network analysis identified VEGFA and FN1 as the key genes with the most interactions. Genes were involved in extracellular matrix regulation, collagen binding and the immune response. Of interest were 10 genes (VEGFA, DCN, SPARC, COL6A2, MGP, CYR61, ANXA6, LGALS1, C1QA and C1QB) directly connected to fibronectin 1 (FN1). Most of these genes were lower expressed in severe COPD and showed the same effect in nasal brushings. Conclusions: We found a unique severe COPD bronchial gene signature with key roles for VEGFA and FN1, which was retained in the upper airways. This supports the hypothesis that severe COPD, at least partly, comprises a different pathology and supports the potential for biomarker development based on nasal brushes in COPD.

IL-33 Expression Is Lower in Current Smokers at both Transcriptomic and Protein Levels.

Rationale: IL-33 is a proinflammatory cytokine thought to play a role in the pathogenesis of asthma and chronic obstructive pulmonary disease (COPD). A recent clinical trial using an anti-IL-33 antibody showed a reduction in exacerbation and improved lung function in ex-smokers but not current smokers with COPD. Objectives: This study aimed to understand the effects of smoking status on IL-33. Methods: We investigated the association of smoking status with the level of gene expression of IL-33 in the airways in eight independent transcriptomic studies of lung airways. Additionally, we performed Western blot analysis and immunohistochemistry for IL-33 in lung tissue to assess protein levels. Measurements and Main Results: Across the bulk RNA-sequencing datasets, IL-33 gene expression and its signaling pathway were significantly lower in current versus former or never-smokers and increased upon smoking cessation (P < 0.05). Single-cell sequencing showed that IL-33 is predominantly expressed in resting basal epithelial cells and decreases during the differentiation process triggered by smoke exposure. We also found a higher transitioning of this cellular subpopulation into a more differentiated cell type during chronic smoking, potentially driving the reduction of IL-33. Protein analysis demonstrated lower IL-33 levels in lung tissue from current versus former smokers with COPD and a lower proportion of IL-33-positive basal cells in current versus ex-smoking controls. Conclusions: We provide strong evidence that cigarette smoke leads to an overall reduction in IL-33 expression in transcriptomic and protein level, and this may be due to the decrease in resting basal cells. Together, these findings may explain the clinical observation that a recent antibody-based anti-IL-33 treatment is more effective in former than current smokers with COPD.

Acute cigarette smoke exposure leads to higher viral infection in human bronchial epithelial cultures by altering interferon, glycolysis and GDF15-related pathways.

BACKGROUND: Acute exacerbations of chronic inflammatory lung diseases, such as chronic obstructive pulmonary disease (COPD), are frequently associated with rhinovirus (RV) infections. Despite these associations, the pathogenesis of virus-induced exacerbations is incompletely understood. We aimed to investigate effects of cigarette smoke (CS), a primary risk factor for COPD, on RV infection in airway epithelium and identify novel mechanisms related to these effects. METHODS: Primary bronchial epithelial cells (PBEC) from COPD patients and controls were differentiated by culture at the air-liquid interface (ALI) and exposed to CS and RV-A16. Bulk RNA sequencing was performed using samples collected at 6 and 24 h post infection (hpi), and viral load, mediator and L-lactate levels were measured at 6, 24 and 48hpi. To further delineate the effect of CS on RV-A16 infection, we performed growth differentiation factor 15 (GDF15) knockdown, L-lactate and interferon pre-treatment in ALI-PBEC. We performed deconvolution analysis to predict changes in the cell composition of ALI-PBEC after the various exposures. Finally, we compared transcriptional responses of ALI-PBEC to those in nasal epithelium after human RV-A16 challenge. RESULTS: CS exposure impaired antiviral responses at 6hpi and increased viral replication at 24 and 48hpi in ALI-PBEC. At 24hpi, CS exposure enhanced expression of RV-A16-induced epithelial interferons, inflammation-related genes and CXCL8. CS exposure increased expression of oxidative stress-related genes, of GDF15, and decreased mitochondrial membrane potential. GDF15 knockdown experiments suggested involvement of this pathway in the CS-induced increase in viral replication. Expression of glycolysis-related genes and L-lactate production were increased by CS exposure, and was demonstrated to contribute to higher viral replication. No major differences were demonstrated between COPD and non-COPD-derived cultures. However, cellular deconvolution analysis predicted higher secretory cells in COPD-derived cultures at baseline. CONCLUSION: Altogether, our findings demonstrate that CS exposure leads to higher viral infection in human bronchial epithelium by altering not only interferon responses, but likely also through a switch to glycolysis, and via GDF15-related pathways.

Expression of targets of the RNA-binding protein AUF-1 in human airway epithelium indicates its role in cellular senescence and inflammation.

Introduction: The RNA-binding protein AU-rich-element factor-1 (AUF-1) participates to posttranscriptional regulation of genes involved in inflammation and cellular senescence, two pathogenic mechanisms of chronic obstructive pulmonary disease (COPD). Decreased AUF-1 expression was described in bronchiolar epithelium of COPD patients versus controls and in vitro cytokine- and cigarette smoke-challenged human airway epithelial cells, prompting the identification of epithelial AUF-1-targeted transcripts and function, and investigation on the mechanism of its loss. Results: RNA immunoprecipitation-sequencing (RIP-Seq) identified, in the human airway epithelial cell line BEAS-2B, 494 AUF-1-bound mRNAs enriched in their 3'-untranslated regions for a Guanine-Cytosine (GC)-rich binding motif. AUF-1 association with selected transcripts and with a synthetic GC-rich motif were validated by biotin pulldown. AUF-1-targets' steady-state levels were equally affected by partial or near-total AUF-1 loss induced by cytomix (TNFα/IL1ß/IFNγ/10 nM each) and siRNA, respectively, with differential transcript decay rates. Cytomix-mediated decrease in AUF-1 levels in BEAS-2B and primary human small-airways epithelium (HSAEC) was replicated by treatment with the senescence- inducer compound etoposide and associated with readouts of cell-cycle arrest, increase in lysosomal damage and senescence-associated secretory phenotype (SASP) factors, and with AUF-1 transfer in extracellular vesicles, detected by transmission electron microscopy and immunoblotting. Extensive in-silico and genome ontology analysis found, consistent with AUF-1 functions, enriched RIP-Seq-derived AUF-1-targets in COPD-related pathways involved in inflammation, senescence, gene regulation and also in the public SASP proteome atlas; AUF-1 target signature was also significantly represented in multiple transcriptomic COPD databases generated from primary HSAEC, from lung tissue and from single-cell RNA-sequencing, displaying a predominant downregulation of expression. Discussion: Loss of intracellular AUF-1 may alter posttranscriptional regulation of targets particularly relevant for protection of genomic integrity and gene regulation, thus concurring to airway epithelial inflammatory responses related to oxidative stress and accelerated aging. Exosomal-associated AUF-1 may in turn preserve bound RNA targets and sustain their function, participating to spreading of inflammation and senescence to neighbouring cells.

CRISPR single base-editing: in silico predictions to variant clonal cell lines.

Engineering single base edits using CRISPR technology including specific deaminases and single-guide RNA (sgRNA) is a rapidly evolving field. Different types of base edits can be constructed, with cytidine base editors (CBEs) facilitating transition of C-to-T variants, adenine base editors (ABEs) enabling transition of A-to-G variants, C-to-G transversion base editors (CGBEs) and recently adenine transversion editors (AYBE) that create A-to-C and A-to-T variants. The base-editing machine learning algorithm BE-Hive predicts which sgRNA and base editor combinations have the strongest likelihood of achieving desired base edits. We have used BE-Hive and TP53 mutation data from The Cancer Genome Atlas (TCGA) ovarian cancer cohort to predict which mutations can be engineered, or reverted to wild-type (WT) sequence, using CBEs, ABEs or CGBEs. We have developed and automated a ranking system to assist in selecting optimally designed sgRNA that considers the presence of a suitable protospacer adjacent motif (PAM), the frequency of predicted bystander edits, editing efficiency and target base change. We have generated single constructs containing ABE or CBE editing machinery, an sgRNA cloning backbone and an enhanced green fluorescent protein tag (EGFP), removing the need for co-transfection of multiple plasmids. We have tested our ranking system and new plasmid constructs to engineer the p53 mutants Y220C, R282W and R248Q into WT p53 cells and shown that these mutants cannot activate four p53 target genes, mimicking the behaviour of endogenous p53 mutations. This field will continue to rapidly progress, requiring new strategies such as we propose to ensure desired base-editing outcomes.

An integrated cell atlas of the lung in health and disease.

Single-cell technologies have transformed our understanding of human tissues. Yet, studies typically capture only a limited number of donors and disagree on cell type definitions. Integrating many single-cell datasets can address these limitations of individual studies and capture the variability present in the population. Here we present the integrated Human Lung Cell Atlas (HLCA), combining 49 datasets of the human respiratory system into a single atlas spanning over 2.4 million cells from 486 individuals. The HLCA presents a consensus cell type re-annotation with matching marker genes, including annotations of rare and previously undescribed cell types. Leveraging the number and diversity of individuals in the HLCA, we identify gene modules that are associated with demographic covariates such as age, sex and body mass index, as well as gene modules changing expression along the proximal-to-distal axis of the bronchial tree. Mapping new data to the HLCA enables rapid data annotation and interpretation. Using the HLCA as a reference for the study of disease, we identify shared cell states across multiple lung diseases, including SPP1+ profibrotic monocyte-derived macrophages in COVID-19, pulmonary fibrosis and lung carcinoma. Overall, the HLCA serves as an example for the development and use of large-scale, cross-dataset organ atlases within the Human Cell Atlas.

Transcriptomic Comparison of Human Peripartum and Dilated Cardiomyopathy Identifies Differences in Key Disease Pathways.

Peripartum cardiomyopathy (PPCM) is a rare form of acute onset heart failure that presents in otherwise healthy pregnant women around the time of delivery. While most of these women respond to early intervention, about 20% progress to end-stage heart failure that symptomatically resembles dilated cardiomyopathy (DCM). In this study, we examined two independent RNAseq datasets from the left ventricle of end-stage PPCM patients and compared gene expression profiles to female DCM and non-failing donors. Differential gene expression, enrichment analysis and cellular deconvolution were performed to identify key processes in disease pathology. PPCM and DCM display similar enrichment in metabolic pathways and extracellular matrix remodeling suggesting these are similar processes across end-stage systolic heart failure. Genes involved in golgi vesicles biogenesis and budding were enriched in PPCM left ventricles compared to healthy donors but were not found in DCM. Furthermore, changes in immune cell populations are evident in PPCM but to a lesser extent compared to DCM, where the latter is associated with pronounced pro-inflammatory and cytotoxic T cell activity. This study reveals several pathways that are common to end-stage heart failure but also identifies potential targets of disease that may be unique to PPCM and DCM.

Smoking increases expression of the SARS-CoV-2 spike protein-binding long ACE2 isoform in bronchial epithelium.

After more than two years the COVID-19 pandemic, that is caused by infection with the respiratory SARS-CoV-2 virus, is still ongoing. The risk to develop severe COVID-19 upon SARS-CoV-2 infection is increased in individuals with a high age, high body mass index, and who are smoking. The SARS-CoV-2 virus infects cells of the upper respiratory tract by entering these cells upon binding to the Angiotensin-converting enzyme 2 (ACE2) receptor. ACE2 is expressed in various cell types in the lung but the expression is especially high in goblet and ciliated cells. Recently, it was shown that next to its full-length isoform, ACE2 also has a short isoform. The short isoform is unable to bind SARS-CoV-2 and does not facilitate viral entry. In the current study we investigated whether active cigarette smoking increases the expression of the long or the short ACE2 isoform. We showed that in active smokers the expression of the long, active isoform, but not the short isoform of ACE2 is higher compared to never smokers. Additionally, it was shown that the expression of especially the long, active isoform of ACE2 was associated with secretory, club and goblet epithelial cells. This study increases our understanding of why current smokers are more susceptible to SARS-CoV-2 infection, in addition to the already established increased risk to develop severe COVID-19.