Consensus Gene Co-Expression Network Analysis Identifies Novel Genes Associated with Severity of Fibrotic Lung Disease

SEMA3B inhibits TGFβ-induced extracellular matrix protein production and its reduced levels are associated with a decline in lung function in IPF

Idiopathic pulmonary fibrosis (IPF) is marked by the activation of fibroblasts, leading to excessive production and deposition of extracellular matrix (ECM) within the lung parenchyma. Despite the pivotal role of ECM overexpression in IPF, potential negative regulators of ECM production in fibroblasts have yet to be identified. Semaphorin class 3B (SEMA3B), a secreted protein highly expressed in lung tissues, has established roles in axonal guidance and tumor suppression. However, the role of SEMA3B in ECM production by fibroblasts in the pathogenesis of IPF remains unexplored. Here, we show the downregulation of SEMA3B and its cognate binding receptor, neuropilin 1 (NRP1), in IPF lungs compared with healthy controls. Notably, the reduced expression of SEMA3B and NRP1 is associated with a decline in lung function in IPF. The downregulation of SEMA3B and NRP1 transcripts was validated in the lung tissues of patients with IPF, and two alternative mouse models of pulmonary fibrosis. In addition, we show that transforming growth factor-β (TGFβ) functions as a negative regulator of SEMA3B and NRP1 expression in lung fibroblasts. Furthermore, we demonstrate the antifibrotic effects of SEMA3B against TGFβ-induced ECM production in IPF lung fibroblasts. Overall, our findings uncovered a novel role of SEMA3B in the pathogenesis of pulmonary fibrosis and provided novel insights into modulating the SEMA3B-NRP1 axis to attenuate pulmonary fibrosis.NEW & NOTEWORTHY The excessive production and secretion of collagens and other extracellular matrix proteins by fibroblasts lead to the scarring of the lung in severe fibrotic lung diseases. This study unveils an antifibrotic role for semaphorin class 3B (SEMA3B) in the pathogenesis of idiopathic pulmonary fibrosis. SEMA3B functions as an inhibitor of transforming growth factor-β-driven fibroblast activation and reduced levels of SEMA3B and its receptor, neuropilin 1, are associated with decreased lung function in idiopathic pulmonary fibrosis.

Genomic analysis of severe COVID-19 considering or not asthma comorbidity: GWAS insights from the BQC19 cohort

BACKGROUND

The severity of COVID-19 is influenced by various factors including the presence of respiratory diseases. Studies have indicated a potential relationship between asthma and COVID-19 severity.

OBJECTIVE

This study aimed to conduct a genome-wide association study (GWAS) to identify genetic and clinical variants associated with the severity of COVID-19, both among patients with and without asthma.

METHODS

We analyzed data from 2131 samples sourced from the Biobanque québécoise de la COVID-19 (BQC19), with 1499 samples from patients who tested positive for COVID-19. Among these, 1110 exhibited mild-to-moderate symptoms, 389 had severe symptoms, and 58 had asthma. We conducted a comparative analysis of clinical data from individuals in these three groups and GWAS using a logistic regression model. Phenotypic data analysis resulted in the refined covariates integrated into logistic models for genetic studies.

RESULTS

Considering a significance threshold of 1 × 10-6, seven genetic variants were associated with severe COVID-19. These variants were located proximal to five genes: sodium voltage-gated channel alpha subunit 1 (SCN10A), desmoplakin (DSP), RP1 axonemal microtubule associated (RP1), IGF like family member 1 (IGFL1), and docking protein 5 (DOK5). The GWAS comparing individuals with severe COVID-19 with asthma to those without asthma revealed four genetic variants in transmembrane protein with EGF like and two follistatin like domains 2 (TMEFF2) and huntingtin interacting protein-1 (HIP1) genes.

CONCLUSION

This study provides significant insights into the genetic profiles of patients with severe forms of the disease, whether accompanied by asthma or not. These findings enhance our comprehension of the genetic factors that affect COVID-19 severity.

KEY MESSAGES

Seven genetic variants were associated with the severe form of COVID-19; Four genetic variants were associated with the severe form of COVID-19 in individuals with comorbid asthma; These findings help define the genetic component of the severe form of COVID-19 in relation to asthma as a comorbidity.

ER exit in physiology and disease

The biosynthetic secretory pathway is comprised of multiple steps, modifications and interactions that form a highly precise pathway of protein trafficking and secretion, that is essential for eukaryotic life. The general outline of this pathway is understood, however the specific mechanisms are still unclear. In the last 15 years there have been vast advancements in technology that enable us to advance our understanding of this complex and subtle pathway. Therefore, based on the strong foundation of work performed over the last 40 years, we can now build another level of understanding, using the new technologies available. The biosynthetic secretory pathway is a high precision process, that involves a number of tightly regulated steps: Protein folding and quality control, cargo selection for Endoplasmic Reticulum (ER) exit, Golgi trafficking, sorting and secretion. When deregulated it causes severe diseases that here we categorise into three main groups of aberrant secretion: decreased, excess and altered secretion. Each of these categories disrupts organ homeostasis differently, effecting extracellular matrix composition, changing signalling events, or damaging the secretory cells due to aberrant intracellular accumulation of secretory proteins. Diseases of aberrant secretion are very common, but despite this, there are few effective therapies. Here we describe ER exit sites (ERES) as key hubs for regulation of the secretory pathway, protein quality control and an integratory hub for signalling within the cell. This review also describes the challenges that will be faced in developing effective therapies, due to the specificity required of potential drug candidates and the crucial need to respect the fine equilibrium of the pathway. The development of novel tools is moving forward, and we can also use these tools to build our understanding of the acute regulation of ERES and protein trafficking. Here we review ERES regulation in context as a therapeutic strategy.

Unagi: Deep Generative Model for Deciphering Cellular Dynamics and In-Silico Drug Discovery in Complex Diseases

Human diseases are characterized by intricate cellular dynamics. Single-cell sequencing provides critical insights, yet a persistent gap remains in computational tools for detailed disease progression analysis and targeted in-silico drug interventions. Here, we introduce UNAGI, a deep generative neural network tailored to analyze time-series single-cell transcriptomic data. This tool captures the complex cellular dynamics underlying disease progression, enhancing drug perturbation modeling and discovery. When applied to a dataset from patients with Idiopathic Pulmonary Fibrosis (IPF), UNAGI learns disease-informed cell embeddings that sharpen our understanding of disease progression, leading to the identification of potential therapeutic drug candidates. Validation via proteomics reveals the accuracy of UNAGI's cellular dynamics analyses, and the use of the Fibrotic Cocktail treated human Precision-cut Lung Slices confirms UNAGI's predictions that Nifedipine, an antihypertensive drug, may have antifibrotic effects on human tissues. UNAGI's versatility extends to other diseases, including a COVID dataset, demonstrating adaptability and confirming its broader applicability in decoding complex cellular dynamics beyond IPF, amplifying its utility in the quest for therapeutic solutions across diverse pathological landscapes.

Wnt signaling regulates ion channel expression to promote smooth muscle and cartilage formation in developing mouse trachea

Ion channels play critical roles in the physiology and function of the nervous system and contractile tissue; however, their role in noncontractile tissue and embryonic development has yet to be understood. Tracheobronchomalacia (TBM) and complete tracheal rings (CTR) are disorders affecting the muscle and cartilage of the trachea and bronchi, whose etiology remains poorly understood. We demonstrated that trachealis muscle organization and polarity are disrupted after epithelial ablation of Wntless (Wls), a cargo receptor critical for the Wnt signaling pathway, in developing trachea. The phenotype resembles the anomalous trachealis muscle observed after deletion of ion channel encoding genes in developing mouse trachea. We sought to investigate whether and how the deletion of Wls affects ion channels during tracheal development. We hypothesize that Wnt signaling influences the expression of ion channels to promote trachealis muscle cell assembly and patterning. Deleting Wls in developing trachea causes differential regulation of genes mediating actin binding, cytoskeleton organization, and potassium ion channel activity. Wnt signaling regulates the expression of Kcnj13, Kcnd3, Kcnj8, and Abcc9 as demonstrated by in vitro studies and in vivo analysis in Wnt5a and β-catenin-deficient tracheas. Pharmacological inhibition of potassium ion channels and Wnt signaling impaired contractility of developing trachealis smooth muscle and formation of cartilaginous mesenchymal condensation. Thus, in mice, epithelial-induced Wnt/β-catenin signaling mediates trachealis muscle and cartilage development via modulation of ion channel expression, promoting trachealis muscle architecture, contractility, and cartilaginous extracellular matrix. In turn, ion channel activity may influence tracheal morphogenesis underlying TBM and CTR.NEW & NOTEWORTHY Ion channels play critical roles in the physiology and function of the nervous system and contractile tissue; however, their role in noncontractile tissue and embryonic development has yet to be understood. In this study, we focused on the role of ion channels in the differentiation and patterning of the large airways of the developing respiratory tract. We identify a mechanism by which Wnt-beta-catenin signaling controls levels of ion channel-encoding genes to promote tracheal differentiation.

Wnt signaling regulates ion channel expression to promote smooth muscle and cartilage formation in developing mouse trachea

Ion channels play critical roles in the physiology and function of the nervous system and contractile tissue; however, their role in non-contractile tissue and embryonic development has yet to be understood. Tracheobronchomalacia (TBM) and complete tracheal rings (CTR) are disorders affecting the muscle and cartilage of the trachea and bronchi, whose etiology remains poorly understood. We demonstrated that trachealis muscle organization and polarity are disrupted after epithelial ablation of Wls, a cargo receptor critical for the Wnt signaling pathway, in developing trachea. The phenotype resembles the anomalous trachealis muscle observed after deletion of ion channel encoding genes in developing mouse trachea. We sought to investigate whether and how the deletion of Wls affects ion channels during tracheal development. We hypothesize that Wnt signaling influences the expression of ion channels to promote trachealis muscle cell assembly and patterning. Deleting Wls in developing trachea causes differential regulation of genes mediating actin binding, cytoskeleton organization, and potassium ion channel activity. Wnt signaling regulated expression of Kcnj13, Kcnd3, Kcnj8, and Abcc9 as demonstrated by in vitro studies and in vivo analysis in Wnt5a and β-catenin deficient tracheas. Pharmacological inhibition of potassium ion channels and Wnt signaling impaired contractility of developing trachealis smooth muscle and formation of cartilaginous mesenchymal condensation. Thus, in mice, epithelial-induced Wnt/β-catenin signaling mediates trachealis muscle and cartilage development via modulation of ion channel expression, promoting trachealis muscle architecture, contractility, and cartilaginous extracellular matrix. In turn, ion channel activity may influence tracheal morphogenesis underlying TBM and CTR.

A fibroblast-dependent TGFβ1/sFRP2 noncanonical Wnt signaling axis underlies epithelial metaplasia in idiopathic pulmonary fibrosis

Reciprocal interactions between alveolar fibroblasts and epithelial cells are crucial for lung homeostasis, injury repair, and fibrogenesis, but underlying mechanisms remain unclear. To investigate this, we administered the fibroblast-selective TGFβ1 signaling inhibitor, epigallocatechin gallate (EGCG), to Interstitial Lung Disease (ILD) patients undergoing diagnostic lung biopsy and conducted single-cell RNA sequencing on spare tissue. Unexposed biopsy samples showed higher fibroblast TGFβ1 signaling compared to non-disease donor or end-stage ILD tissues. In vivo, EGCG significantly downregulated TGFβ1 signaling and several pro-inflammatory and stress pathways in biopsy samples. Notably, EGCG reduced fibroblast secreted Frizzle-like Receptor Protein 2 (sFRP2), an unrecognized TGFβ1 fibroblast target gene induced near type II alveolar epithelial cells (AEC2s). In human AEC2-fibroblast coculture organoids, sFRP2 was essential for AEC2 trans-differentiation to basal cells. Precision cut lung slices (PCLS) from normal donors demonstrated that TGFβ1 promoted KRT17 expression and AEC2 morphological change, while sFRP2 was necessary for KRT5 expression in AEC2-derived basaloid cells. Wnt-receptor Frizzled 5 (Fzd5) expression and downstream calcineurin-related signaling in AEC2s were required for sFRP2-induced KRT5 expression. These findings highlight stage-specific TGFβ1 signaling in ILD, the therapeutic potential of EGCG in reducing IPF-related transcriptional changes, and identify the TGFβ1-non-canonical Wnt pathway crosstalk via sFRP2 as a novel mechanism for dysfunctional epithelial signaling in Idiopathic Pulmonary Fibrosis/ILD.

Wilms Tumor 1-Driven Fibroblast Activation and Subpleural Thickening in Idiopathic Pulmonary Fibrosis

Idiopathic pulmonary fibrosis (IPF) is a progressive fibrotic lung disease that is often fatal due to the formation of irreversible scar tissue in the distal areas of the lung. Although the pathological and radiological features of IPF lungs are well defined, the lack of insight into the fibrogenic role of fibroblasts that accumulate in distinct anatomical regions of the lungs is a critical knowledge gap. Fibrotic lesions have been shown to originate in the subpleural areas and extend into the lung parenchyma through processes of dysregulated fibroproliferation, migration, fibroblast-to-myofibroblast transformation, and extracellular matrix production. Identifying the molecular targets underlying subpleural thickening at the early and late stages of fibrosis could facilitate the development of new therapies to attenuate fibroblast activation and improve the survival of patients with IPF. Here, we discuss the key cellular and molecular events that contribute to (myo)fibroblast activation and subpleural thickening in IPF. In particular, we highlight the transcriptional programs involved in mesothelial to mesenchymal transformation and fibroblast dysfunction that can be targeted to alter the course of the progressive expansion of fibrotic lesions in the distal areas of IPF lungs.

Neonatal hyperoxia induces activated pulmonary cellular states and sex-dependent transcriptomic changes in a model of experimental bronchopulmonary dysplasia

Hyperoxia disrupts lung development in mice and causes bronchopulmonary dysplasia (BPD) in neonates. To investigate sex-dependent molecular and cellular programming involved in hyperoxia, we surveyed the mouse lung using single cell RNA sequencing (scRNA-seq), and validated our findings in human neonatal lung cells in vitro. Hyperoxia-induced inflammation in alveolar type (AT) 2 cells gave rise to damage-associated transient progenitors (DATPs). It also induced a new subpopulation of AT1 cells with reduced expression of growth factors normally secreted by AT1 cells, but increased mitochondrial gene expression. Female alveolar epithelial cells had less EMT and pulmonary fibrosis signaling in hyperoxia. In the endothelium, expansion of Car4+ EC (Cap2) was seen in hyperoxia along with an emergent subpopulation of Cap2 with repressed VEGF signaling. This regenerative response was increased in females exposed to hyperoxia. Mesenchymal cells had inflammatory signatures in hyperoxia, with a new distal interstitial fibroblast subcluster characterized by repressed lipid biosynthesis and a transcriptomic signature resembling myofibroblasts. Hyperoxia-induced gene expression signatures in human neonatal fibroblasts and alveolar epithelial cells in vitro resembled mouse scRNA-seq data. These findings suggest that neonatal exposure to hyperoxia programs distinct sex-specific stem cell progenitor and cellular reparative responses that underpin lung remodeling in BPD.

A potential three-gene-based diagnostic signature for idiopathic pulmonary fibrosis

Background: Idiopathic pulmonary fibrosis (IPF) is a life-threatening disease whose etiology remains unknown. This study aims to explore diagnostic biomarkers and pathways involved in IPF using bioinformatics analysis. Methods: IPF-related gene expression datasets were retrieved and downloaded from the NCBI Gene Expression Omnibus database. Differentially expressed genes (DEGs) were screened, and weighted correlation network analysis (WGCNA) was performed to identify key module and genes. Functional enrichment analysis was performed on genes in the clinically significant module. Then least absolute shrinkage and selection operator (LASSO) logistic regression and support vector machine-recursive feature elimination (SVM-RFE) algorithms were run to screen candidate biomarkers. The expression and diagnostic value of the biomarkers in IPF were further validated in external test datasets (GSE110147). Results: 292 samples and 1,163 DEGs were screened to construct WGCNA. In WGCNA, the blue module was identified as the key module, and 59 genes in this module correlated highly with IPF. Functional enrichment analysis of blue module genes revealed the importance of extracellular matrix-associated pathways in IPF. IL13RA2, CDH3, and COMP were identified as diagnostic markers of IPF via LASSO and SVM-RFE. These genes showed good diagnostic value for IPF and were significantly upregulated in IPF. Conclusion: This study indicates that IL13RA2, CDH3, and COMP could serve as diagnostic signature for IPF and might offer new insights in the underlying diagnosis of IPF.

Gastroesophageal Reflux Disease in Idiopathic Pulmonary Fibrosis: Viewer or Actor? To Treat or Not to Treat?

Background: Idiopathic pulmonary fibrosis (IPF) is a rare and severe disease with a median survival of ∼3 years. Several risk factors have been identified, such as age, genetic predisposition, tobacco exposure, and gastro-oesophageal reflux disease (GERD). Prevalence of GERD in IPF is high and may affect 87% of patients, of whom only half (47%) report symptoms. Objective: The aim of this study is to review current evidence regarding the correlation between GERD and IPF and to evaluate the current studies regarding treatments for GERD-IPF. Methods: A review to identify research papers documenting an association between GERD and IPF was performed. Results: We identified several studies that have confirmed the association between GERD and IPF, with an increased acid exposure, risk of gastric aspiration and bile acids levels in these patients. Few studies focused their attention on GERD treatment, showing how antiacid therapy was not able to change IPF evolution. Conclusions: This review investigating the correlation between GERD and IPF has confirmed the hypothesized association. However, further large prospective studies are needed to corroborate and elucidate these findings with a focus on preventative and treatment strategies.

Spatially resolved deconvolution of the fibrotic niche in lung fibrosis

A defining pathological feature of human lung fibrosis is localized tissue heterogeneity, which challenges the interpretation of transcriptomic studies that typically lose spatial information. Here we investigate spatial gene expression in diagnostic tissue using digital profiling technology. We identify distinct, region-specific gene expression signatures as well as shared gene signatures. By integration with single-cell data, we spatially map the cellular composition within and distant from the fibrotic niche, demonstrating discrete changes in homeostatic and pathologic cell populations even in morphologically preserved lung, while through ligand-receptor analysis, we investigate cellular cross-talk within the fibrotic niche. We confirm findings through bioinformatic, tissue, and in vitro analyses, identifying that loss of NFKB inhibitor zeta in alveolar epithelial cells dysregulates the TGFβ/IL-6 signaling axis, which may impair homeostatic responses to environmental stress. Thus, spatially resolved deconvolution advances understanding of cell composition and microenvironment in human lung fibrogenesis.