Gene network approach reveals co-expression patterns in nasal and bronchial epithelium

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.

The functional role of CST1 and CCL26 in asthma development

BACKGROUND

Asthma is the most common chronic disease in children with an increasing prevalence. Its development is caused by genetic and environmental factors and allergic sensitization is a known trigger. Dog allergens affect up to 30% of all children and dog dander-sensitized children show increased expression of cystatin-1 (CST1) and eotaxin-3 (CCL26) in nasal epithelium. The aim of our study was to investigate the functional mechanism of CST1 and CCL26 in the alveolar basal epithelial cell line A549.

METHODS

A549 cells were transfected with individual overexpression vectors for CST1 and CCL26 and RNA sequencing was performed to examine the transcriptomics. edgeR was used to identify differentially expressed genes (= DEG, |log2 FC | ≥ 2, FDR < 0.01). The protein expression levels of A549 cells overexpressing CST1 and CCL26 were analyzed using the Target 96 inflammation panel from OLINK (antibody-mediated proximity extension-based assay; OLINK Proteomics). Differentially expressed proteins were considered with a |log2 FC| ≥ 1, p < .05.

RESULTS

The overexpression of CST1 resulted in a total of 27 DEG (1 upregulated and 26 downregulated) and the overexpression of CCL26 in a total of 137 DEG (0 upregulated and 137 downregulated). The gene ontology enrichment analysis showed a significant downregulation of type I and III interferon signaling pathway genes as well as interferon-stimulated genes. At the protein level, overexpression of CST1 induced a significantly increased expression of CCL3, whereas CCL26 overexpression led to increased expression of HGF, and a decrease of CXCL11, CCL20, CCL3 and CXCL10.

CONCLUSION

Our results indicate that an overexpression of CST1 and CCL26 cause a downregulation of interferon related genes and inflammatory proteins. It might cause a higher disease susceptibility, mainly for allergic asthma, as CCL26 is an agonist for CCR-3-carrying cells, such as eosinophils and Th2 lymphocytes, mostly active in allergic asthma.

An epithelial gene signature of trans-IL-6 signaling defines a subgroup of type 2-low asthma

BACKGROUND

Asthma is stratified into type 2-high and type 2-low inflammatory phenotypes. Limited success has been achieved in developing drugs that target type 2-low inflammation. Previous studies have linked IL-6 signaling to severe asthma. IL-6 cooperates with soluble-IL-6Rα to activate cell signaling in airway epithelium.

OBJECTIVE

We sought to study the role of sIL-6Rα amplified IL-6 signaling in airway epithelium and to develop an IL-6+ sIL-6Rα gene signature that may be used to select asthma patients who potentially respond to anti-IL-6 therapy.

METHODS

Human airway epithelial cells were stimulated with combinations of IL-6, sIL-6Rα, and inhibitors, sgp130 (Olamkicept), and anti-IL-6R (Tocilizumab), to assess effects on pathway activation, epithelial barrier integrity, and gene expression. A gene signature was generated to identify IL-6 high patients using bronchial biopsies and nasal brushes.

RESULTS

Soluble-IL-6Rα amplified the activation of the IL-6 pathway, shown by the increase of STAT3 phosphorylation and stronger gene induction in airway epithelial cells compared to IL-6 alone. Olamkicept and Tocilizumab inhibited the effect of IL-6 + sIL-6Rα on gene expression. We developed an IL-6 + sIL-6Rα gene signature and observed enrichment of this signature in bronchial biopsies but not nasal brushes from asthma patients compared to healthy controls. An IL-6 + sIL-6Rα gene signature score was associated with lower levels of sputum eosinophils in asthma.

CONCLUSION

sIL-6Rα amplifies IL-6 signaling in bronchial epithelial cells. Higher local airway IL-6 + sIL-6Rα signaling is observed in asthma patients with low sputum eosinophils.

Exploring intrinsic variability between cultured nasal and bronchial epithelia in cystic fibrosis

The nasal and bronchial epithelium are unified parts of the respiratory tract that are affected in the monogenic disorder cystic fibrosis (CF). Recent studies have uncovered that nasal and bronchial tissues exhibit intrinsic variability, including differences in mucociliary cell composition and expression of unique transcriptional regulatory proteins which relate to germ layer origin. In the present study, we explored whether intrinsic differences between nasal and bronchial epithelial cells persist in cell cultures and affect epithelial cell functioning in CF. Comparison of air-liquid interface (ALI) differentiated epithelial cells from subjects with CF revealed distinct mucociliary differentiation states of nasal and bronchial cultures. Moreover, using RNA sequencing we identified cell type-specific signature transcription factors in differentiated nasal and bronchial epithelial cells, some of which were already poised for expression in basal progenitor cells as evidenced by ATAC sequencing. Analysis of differentiated nasal and bronchial epithelial 3D organoids revealed distinct capacities for fluid secretion, which was linked to differences in ciliated cell differentiation. In conclusion, we show that unique phenotypical and functional features of nasal and bronchial epithelial cells persist in cell culture models, which can be further used to investigate the effects of tissue-specific features on upper and lower respiratory disease development in CF.

A persistent neutrophil-associated immune signature characterizes post-COVID-19 pulmonary sequelae

Interstitial lung disease and associated fibrosis occur in a proportion of individuals who have recovered from severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection through unknown mechanisms. We studied individuals with severe coronavirus disease 2019 (COVID-19) after recovery from acute illness. Individuals with evidence of interstitial lung changes at 3 to 6 months after recovery had an up-regulated neutrophil-associated immune signature including increased chemokines, proteases, and markers of neutrophil extracellular traps that were detectable in the blood. Similar pathways were enriched in the upper airway with a concomitant increase in antiviral type I interferon signaling. Interaction analysis of the peripheral phosphoproteome identified enriched kinases critical for neutrophil inflammatory pathways. Evaluation of these individuals at 12 months after recovery indicated that a subset of the individuals had not yet achieved full normalization of radiological and functional changes. These data provide insight into mechanisms driving development of pulmonary sequelae during and after COVID-19 and provide a rational basis for development of targeted approaches to prevent long-term complications.

GeneNetTools: tests for Gaussian graphical models with shrinkage

MOTIVATION

Gaussian graphical models (GGMs) are network representations of random variables (as nodes) and their partial correlations (as edges). GGMs overcome the challenges of high-dimensional data analysis by using shrinkage methodologies. Therefore, they have become useful to reconstruct gene regulatory networks from gene-expression profiles. However, it is often ignored that the partial correlations are 'shrunk' and that they cannot be compared/assessed directly. Therefore, accurate (differential) network analyses need to account for the number of variables, the sample size, and also the shrinkage value, otherwise, the analysis and its biological interpretation would turn biased. To date, there are no appropriate methods to account for these factors and address these issues.

RESULTS

We derive the statistical properties of the partial correlation obtained with the Ledoit-Wolf shrinkage. Our result provides a toolbox for (differential) network analyses as (i) confidence intervals, (ii) a test for zero partial correlation (null-effects) and (iii) a test to compare partial correlations. Our novel (parametric) methods account for the number of variables, the sample size and the shrinkage values. Additionally, they are computationally fast, simple to implement and require only basic statistical knowledge. Our simulations show that the novel tests perform better than DiffNetFDR-a recently published alternative-in terms of the trade-off between true and false positives. The methods are demonstrated on synthetic data and two gene-expression datasets from Escherichia coli and Mus musculus.

AVAILABILITY AND IMPLEMENTATION

The R package with the methods and the R script with the analysis are available in https://github.com/V-Bernal/GeneNetTools.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Bone Marrow Derived Mesenchymal Stromal Cells Promote Vascularization and Ciliation in Airway Mucosa Tri-Culture Models in Vitro

Patients suffering from irresectable tracheal stenosis often face limited treatment options associated with low quality of life. To date, an optimal tracheal replacement strategy does not exist. A tissue-engineered tracheal substitute promises to overcome limitations such as implant vascularization, functional mucociliary clearance and mechanical stability. In order to advance a tracheal mucosa model recently developed by our group, we examined different supporting cell types in fibrin-based tri-culture with primary human umbilical vein endothelial cells (HUVEC) and primary human respiratory epithelial cells (HRE). Bone marrow-derived mesenchymal stromal cells (BM-MSC), adipose-derived mesenchymal stromal cells (ASC) and human nasal fibroblasts (HNF) were compared regarding their ability to promote mucociliary differentiation and vascularization in vitro. Three-dimensional co-cultures of the supporting cell types with either HRE or HUVEC were used as controls. Mucociliary differentiation and formation of vascular-like structures were analyzed by scanning electron microscopy (SEM), periodic acid Schiff’s reaction (PAS reaction), two-photon laser scanning microscopy (TPLSM) and immunohistochemistry. Cytokine levels of vascular endothelial growth factor (VEGF), epidermal growth factor (EGF), interleukin-6 (IL6), interleukin-8 (IL8), angiopoietin 1, angiopoietin 2, fibroblast growth factor basic (FGF-b) and placenta growth factor (PIGF) in media supernatant were investigated using LEGENDplex™ bead-based immunoassay. Epithelial morphology of tri-cultures with BM-MSC most closely resembled native respiratory epithelium with respect to ciliation, mucus production as well as expression and localization of epithelial cell markers pan-cytokeratin, claudin-1, α-tubulin and mucin5AC. This was followed by tri-cultures with HNF, while ASC-supported tri-cultures lacked mucociliary differentiation. For all supporting cell types, a reduced ciliation was observed in tri-cultures compared to the corresponding co-cultures. Although formation of vascular-like structures was confirmed in all cultures, vascular networks in BM-MSC-tri-cultures were found to be more branched and extended. Concentrations of pro-angiogenic and inflammatory cytokines, in particular VEGF and angiopoietin 2, revealed to be reduced in tri-cultures compared to co-cultures. With these results, our study provides an important step towards a vascularized and ciliated tissue-engineered tracheal replacement. Additionally, our tri-culture model may in the future contribute to an improved understanding of cell-cell interactions in diseases associated with impaired mucosal function.

S100A7 Co-localization and Up-regulation of Filaggrin in Human Sinonasal Epithelial Cells

OBJECTIVE

Filaggrin (FLG) is a protein expressed in the epidermis and involved in the maintenance of the epidermal barrier. However, the expression and localization of FLG in the upper airway remain controversial. The present study aimed to determine the significance of FLG and the effect of S100A7 on FLG expression in the upper respiratory mucosa.

METHODS

Human nasal epithelial cells (HNECs) were cultured and examined for FLG expression and S100A7 effects by real-time polymerase chain reaction and Western blotting. The localization and distribution of FLG were assessed using sinonasal mucosa.

RESULTS

A significant expression of FLG was detected at the mRNA and protein levels in HNECs. A moderate FLG immunoreactivity was observed in the epithelial cells, but no staining was seen in epithelial goblet cells. S100A7 increased the FLG mRNA level in HNECs in a dose-dependent manner and also up-regulated the FLG protein in a dose-dependent manner.

CONCLUSION

This study significantly contributes to a better understanding of the role of FLG in the pathogenesis of airway inflammation from the viewpoint of the epithelial barrier function. FLG-related events in response to S100A7 protein may represent novel therapeutic targets for the treatment of upper airway inflammation.

Small Airway Susceptibility to Chemical and Particle Injury

Small airways (SA) in humans are commonly defined as those conducting airways <2 mm in diameter. They are susceptible to particle- and chemical-induced injury and play a major role in the development of airway disease such as COPD and asthma. Susceptibility to injury can be attributed in part to structural features including airflow dynamics and tissue architecture, but recent evidence may indicate a more prominent role for cellular composition in directing toxicological responses. Animal studies support the hypothesis that inherent cellular differences across the tracheobronchial tree, including metabolic CYP450 expression in the distal conducting airways, can influence SA susceptibility to injury. Currently, there is insufficient information in humans to make similar conclusions, prompting further necessary work in this area. An understanding of why the SA are more susceptible to certain chemical and particle exposures than other airway regions is fundamental to our ability to identify hazardous materials, their properties, and accompanying exposure scenarios that compromise lung function. It is also important for the ability to develop appropriate models for toxicity testing. Moreover, it is central to our understanding of SA disease aetiology and how interventional strategies for treatment may be developed. In this review, we will document the structural and cellular airway regional differences that are likely to influence airway susceptibility to injury, including the role of secretory club cells. We will also describe recent advances in single-cell sequencing of human airways, which have provided unprecedented details of cell phenotype, likely to impact airway chemical and particle injury.

Role of sex hormones in lung cancer

Lung cancer represents the world's leading cause of cancer deaths. Sex differences in the incidence and mortality rates for various types of lung cancers have been identified, but the biological and endocrine mechanisms implicated in these disparities have not yet been determined. While some cancers such as lung adenocarcinoma are more commonly found among women than men, others like squamous cell carcinoma display the opposite pattern or show no sex differences. Associations of tobacco product use rates, susceptibility to carcinogens, occupational exposures, and indoor and outdoor air pollution have also been linked to differential rates of lung cancer occurrence and mortality between sexes. While roles for sex hormones in other types of cancers affecting women or men have been identified and described, little is known about the influence of sex hormones in lung cancer. One potential mechanism identified to date is the synergism between estrogen and some tobacco compounds, and oncogene mutations, in inducing the expression of metabolic enzymes, leading to enhanced formation of reactive oxygen species and DNA adducts, and subsequent lung carcinogenesis. In this review, we present the literature available regarding sex differences in cancer rates, associations of male and female sex hormones with lung cancer, the influence of exogenous hormone therapy in women, and potential mechanisms mediated by male and female sex hormone receptors in lung carcinogenesis. The influence of biological sex on lung disease has recently been established, thus new research incorporating this variable will shed light on the mechanisms behind the observed disparities in lung cancer rates, and potentially lead to the development of new therapeutics to treat this devastating disease.

The 'un-shrunk' partial correlation in Gaussian graphical models

BACKGROUND

In systems biology, it is important to reconstruct regulatory networks from quantitative molecular profiles. Gaussian graphical models (GGMs) are one of the most popular methods to this end. A GGM consists of nodes (representing the transcripts, metabolites or proteins) inter-connected by edges (reflecting their partial correlations). Learning the edges from quantitative molecular profiles is statistically challenging, as there are usually fewer samples than nodes ('high dimensional problem'). Shrinkage methods address this issue by learning a regularized GGM. However, it remains open to study how the shrinkage affects the final result and its interpretation.

RESULTS

We show that the shrinkage biases the partial correlation in a non-linear way. This bias does not only change the magnitudes of the partial correlations but also affects their order. Furthermore, it makes networks obtained from different experiments incomparable and hinders their biological interpretation. We propose a method, referred to as 'un-shrinking' the partial correlation, which corrects for this non-linear bias. Unlike traditional methods, which use a fixed shrinkage value, the new approach provides partial correlations that are closer to the actual (population) values and that are easier to interpret. This is demonstrated on two gene expression datasets from Escherichia coli and Mus musculus.

CONCLUSIONS

GGMs are popular undirected graphical models based on partial correlations. The application of GGMs to reconstruct regulatory networks is commonly performed using shrinkage to overcome the 'high-dimensional problem'. Besides it advantages, we have identified that the shrinkage introduces a non-linear bias in the partial correlations. Ignoring this type of effects caused by the shrinkage can obscure the interpretation of the network, and impede the validation of earlier reported results.

Network and co-expression analysis of airway smooth muscle cell transcriptome delineates potential gene signatures in asthma

Airway smooth muscle (ASM) is known for its role in asthma exacerbations characterized by acute bronchoconstriction and remodeling. The molecular mechanisms underlying multiple gene interactions regulating gene expression in asthma remain elusive. Herein, we explored the regulatory relationship between ASM genes to uncover the putative mechanism underlying asthma in humans. To this end, the gene expression from human ASM was measured with RNA-Seq in non-asthmatic and asthmatic groups. The gene network for the asthmatic and non-asthmatic group was constructed by prioritizing differentially expressed genes (DEGs) (121) and transcription factors (TFs) (116). Furthermore, we identified differentially connected or co-expressed genes in each group. The asthmatic group showed a loss of gene connectivity due to the rewiring of major regulators. Notably, TFs such as ZNF792, SMAD1, and SMAD7 were differentially correlated in the asthmatic ASM. Additionally, the DEGs, TFs, and differentially connected genes over-represented in the pathways involved with herpes simplex virus infection, Hippo and TGF-β signaling, adherens junctions, gap junctions, and ferroptosis. The rewiring of major regulators unveiled in this study likely modulates the expression of gene-targets as an adaptive response to asthma. These multiple gene interactions pointed out novel targets and pathways for asthma exacerbations.

A Single-Cell Atlas of the Human Healthy Airways

Rationale: The respiratory tract constitutes an elaborate line of defense that is based on a unique cellular ecosystem.Objectives: We aimed to investigate cell population distributions and transcriptional changes along the airways by using single-cell RNA profiling.Methods: We have explored the cellular heterogeneity of the human airway epithelium in 10 healthy living volunteers by single-cell RNA profiling. A total of 77,969 cells were collected at 35 distinct locations, from the nose to the 12th division of the airway tree.Measurements and Main Results: The resulting atlas is composed of a high percentage of epithelial cells (89.1%) but also immune (6.2%) and stromal (4.7%) cells with distinct cellular proportions in different regions of the airways. It reveals differential gene expression between identical cell types (suprabasal, secretory, and multiciliated cells) from the nose (MUC4, PI3, SIX3) and tracheobronchial (SCGB1A1, TFF3) airways. By contrast, cell-type-specific gene expression is stable across all tracheobronchial samples. Our atlas improves the description of ionocytes, pulmonary neuroendocrine cells, and brush cells and identifies a related population of NREP-positive cells. We also report the association of KRT13 with dividing cells that are reminiscent of previously described mouse "hillock" cells and with squamous cells expressing SCEL and SPRR1A/B.Conclusions: Robust characterization of a single-cell cohort in healthy airways establishes a valuable resource for future investigations. The precise description of the continuum existing from the nasal epithelium to successive divisions of the airways and the stable gene expression profile of these regions better defines conditions under which relevant tracheobronchial proxies of human respiratory diseases can be developed.

Determinants of SARS-CoV-2 receptor gene expression in upper and lower airways

The recent outbreak of the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), which causes coronavirus disease 2019 (COVID-19), has led to a worldwide pandemic. One week after initial symptoms develop, a subset of patients progresses to severe disease, with high mortality and limited treatment options. To design novel interventions aimed at preventing spread of the virus and reducing progression to severe disease, detailed knowledge of the cell types and regulating factors driving cellular entry is urgently needed. Here we assess the expression patterns in genes required for COVID-19 entry into cells and replication, and their regulation by genetic, epigenetic and environmental factors, throughout the respiratory tract using samples collected from the upper (nasal) and lower airways (bronchi). Matched samples from the upper and lower airways show a clear increased expression of these genes in the nose compared to the bronchi and parenchyma. Cellular deconvolution indicates a clear association of these genes with the proportion of secretory epithelial cells. Smoking status was found to increase the majority of COVID-19 related genes including ACE2 and TMPRSS2 but only in the lower airways, which was associated with a significant increase in the predicted proportion of goblet cells in bronchial samples of current smokers. Both acute and second hand smoke were found to increase ACE2 expression in the bronchus. Inhaled corticosteroids decrease ACE2 expression in the lower airways. No significant effect of genetics on ACE2 expression was observed, but a strong association of DNA- methylation with ACE2 and TMPRSS2- mRNA expression was identified in the bronchus.