Integrative analysis of lung molecular signatures reveals key drivers of idiopathic pulmonary fibrosis

High-resolution transcriptomics analysis of CXCL13+ EPSTI1+ CDK1+ cells with a specific focus on lung adenocarcinoma

Background

Programmed cell death ligand 1 (PD-L1) blocking therapy has transformed the treatment of lung adenocarcinoma (LUAD), which has significantly changed the landscape of immunotherapy. We aimed to explore specific cell subpopulations to understand tumor progression and identify markers of response to PD-L1 blocking therapy.

Methods

Bulk, fluorescence-activated cell sorting (FACS), and single-cell RNA (scRNA) sequencing were used to profile CXCL13, EPSTI1, and CDK1. The gene set variation analysis (GSVA) R package was utilized for score calculation, and prognostic analyses included receiver operating characteristic (ROC) curves, Cox proportional hazard models, and meta-analysis. Additionally, we analyzed tumor microenvironment (TME), genomics, compound perturbations, and clinical indicators. The high-dimensional analysis captured the intrinsic characteristics of the subpopulation. Furthermore, subpopulation differential genes were used for enrichment analysis of transcription factors and compounds.

Results

Literature and website analyses supported the essential role of CXCL13, CDK1, and EPSTI1 in immunotherapy. This led us to focus specifically on LUAD by representing a pan-cancer profile of immune-sensitive genes. Logically, the high-characteristic population may consist of samples positive for CXCL13, EPSTI1, and CDK1. The three-gene signature was a favorable indicator of immunotherapy response in the Stand Up to Cancer-Mark Foundation (SU2C-MARK) LUAD cohort but showed a poor prognosis before treatment in the Lung Cancer Explorer (LCE) database. Further mechanistic exploration revealed specific mutations associated with the three-gene signature in SU2C-MARK LUAD, such as STK11. In The Cancer Genome Atlas (TCGA)-LUAD cohort, the high-scoring group exhibited a higher tumor mutational burden (TMB) and global methylation but a lower fraction genome altered (FGA) and estimated tumor purity. Moreover, dasatinib demonstrated sensitivity in the high-scoring group. The co-localization of the CXCL13, EPSTI1, and CDK1 subpopulation was validated through spatial transcriptome and immunohistochemical databases. Assessment of the subpopulation depicted high-resolution intercellular communication. Maintenance of specific pathways, such as TNF, CD74, and CD44, contributed to immunotherapy sensitivity. Finally, the subpopulation-enriched targets and drugs were confirmed through ConnectivityMap (CMAP) analysis and multi-omics, respectively.

Conclusions

In this study, positive samples for CXCL13, EPSTI1, and CDK1 exhibited poor prognostic significance in treatment-naïve LUAD cases but demonstrated benefits from PD-L1 blockade and dasatinib therapies.

Tree shrews as a new animal model for systemic sclerosis research

Objective

Systemic sclerosis (SSc) is a chronic systemic disease characterized by immune dysregulation and fibrosis for which there is no effective treatment. Animal models are crucial for advancing SSc research. Tree shrews are genetically, anatomically, and immunologically closer to humans than rodents. Thus, the tree shrew model provides a unique opportunity for translational research in SSc.

Methods

In this study, a SSc tree shrew model was constructed by subcutaneous injection of different doses of bleomycin (BLM) for 21 days. We assessed the degree of inflammation and fibrosis in the skin and internal organs, and antibodies in serum. Furthermore, RNA sequencing and a series of bioinformatics analyses were performed to analyze the transcriptome changes, hub genes and immune infiltration in the skin tissues of BLM induced SSc tree shrew models. Multiple sequence alignment was utilized to analyze the conservation of selected target genes across multiple species.

Results

Subcutaneous injection of BLM successfully induced a SSc model in tree shrew. This model exhibited inflammation and fibrosis in skin and lung, and some developed esophageal fibrosis and secrum autoantibodies including antinuclear antibodies and anti-scleroderma-70 antibody. Using RNA sequencing, we compiled skin transcriptome profiles in SSc tree shrew models. 90 differentially expressed genes (DEGs) were identified, which were mainly enriched in the PPAR signaling pathway, tyrosine metabolic pathway, p53 signaling pathway, ECM receptor interaction and glutathione metabolism, all of which are closely associated with SSc. Immune infiltration analysis identified 20 different types of immune cells infiltrating the skin of the BLM-induced SSc tree shrew models and correlations between those immune cells. By constructing a protein-protein interaction (PPI) network, we identified 10 hub genes that were significantly highly expressed in the skin of the SSc models compared to controls. Furthermore, these genes were confirmed to be highly conserved in tree shrews, humans and mice.

Conclusion

This study for the first time comfirmed that tree shrew model of SSc can be used as a novel and promising experimental animal model to study the pathogenesis and translational research in SSc.

Zinc-Dependent Histone Deacetylases in Lung Endothelial Pathobiology

A monolayer of endothelial cells (ECs) lines the lumen of blood vessels and, as such, provides a semi-selective barrier between the blood and the interstitial space. Compromise of the lung EC barrier due to inflammatory or toxic events may result in pulmonary edema, which is a cardinal feature of acute lung injury (ALI) and its more severe form, acute respiratory distress syndrome (ARDS). The EC functions are controlled, at least in part, via epigenetic mechanisms mediated by histone deacetylases (HDACs). Zinc-dependent HDACs represent the largest group of HDACs and are activated by Zn2+. Members of this HDAC group are involved in epigenetic regulation primarily by modifying the structure of chromatin upon removal of acetyl groups from histones. In addition, they can deacetylate many non-histone histone proteins, including those located in extranuclear compartments. Recently, the therapeutic potential of inhibiting zinc-dependent HDACs for EC barrier preservation has gained momentum. However, the role of specific HDAC subtypes in EC barrier regulation remains largely unknown. This review aims to provide an update on the role of zinc-dependent HDACs in endothelial dysfunction and its related diseases. We will broadly focus on biological contributions, signaling pathways and transcriptional roles of HDACs in endothelial pathobiology associated mainly with lung diseases, and we will discuss the potential of their inhibitors for lung injury prevention.

Development of a Novel Biomarker for the Progression of Idiopathic Pulmonary Fibrosis

The progression of idiopathic pulmonary fibrosis (IPF) is diverse and unpredictable. We identified and validated a new biomarker for IPF progression. To identify a candidate gene to predict progression, we assessed differentially expressed genes in patients with advanced IPF compared with early IPF and controls in three lung sample cohorts. Candidate gene expression was confirmed using immunohistochemistry and Western blotting of lung tissue samples from an independent IPF clinical cohort. Biomarker potential was assessed using an enzyme-linked immunosorbent assay of serum samples from the retrospective validation cohort. We verified that the final candidate gene reflected the progression of IPF in a prospective validation cohort. In the RNA-seq comparative analysis of lung tissues, CD276, COL7A1, CTSB, GLI2, PIK3R2, PRAF2, IGF2BP3, and NUPR1 were up-regulated, and ADAMTS8 was down-regulated in the samples of advanced IPF. Only CTSB showed significant differences in expression based on Western blotting (n = 12; p < 0.001) and immunohistochemistry between the three groups of the independent IPF cohort. In the retrospective validation cohort (n = 78), serum CTSB levels were higher in the progressive group (n = 25) than in the control (n = 29, mean 7.37 ng/mL vs. 2.70 ng/mL, p < 0.001) and nonprogressive groups (n = 24, mean 7.37 ng/mL vs. 2.56 ng/mL, p < 0.001). In the prospective validation cohort (n = 129), serum CTSB levels were higher in the progressive group than in the nonprogressive group (mean 8.30 ng/mL vs. 3.00 ng/mL, p < 0.001). After adjusting for baseline FVC, we found that CTSB was independently associated with IPF progression (adjusted OR = 2.61, p < 0.001). Serum CTSB levels significantly predicted IPF progression (AUC = 0.944, p < 0.001). Serum CTSB level significantly distinguished the progression of IPF from the non-progression of IPF or healthy control.

Pulmonary fibrosis: from pathogenesis to clinical decision-making

Pulmonary fibrosis (PF) encompasses a spectrum of chronic lung diseases that progressively impact the interstitium, resulting in compromised gas exchange, breathlessness, diminished quality of life (QoL), and ultimately respiratory failure and mortality. Various diseases can cause PF, with their underlying causes primarily affecting the lung interstitium, leading to their referral as interstitial lung diseases (ILDs). The current understanding is that PF arises from abnormal wound healing processes triggered by various factors specific to each disease, leading to excessive inflammation and fibrosis. While significant progress has been made in understanding the molecular mechanisms of PF, its pathogenesis remains elusive. This review provides an in-depth exploration of the latest insights into PF pathophysiology, diagnosis, treatment, and future perspectives.

Integrating inflammatory biomarker analysis and artificial-intelligence-enabled image-based profiling to identify drug targets for intestinal fibrosis

Intestinal fibrosis, often caused by inflammatory bowel disease, can lead to intestinal stenosis and obstruction, but there are no approved treatments. Drug discovery has been hindered by the lack of screenable cellular phenotypes. To address this, we used a scalable image-based morphology assay called Cell Painting, augmented with machine learning algorithms, to identify small molecules that could reverse the activated fibrotic phenotype of intestinal myofibroblasts. We then conducted a high-throughput small molecule chemogenomics screen of approximately 5,000 compounds with known targets or mechanisms, which have achieved clinical stage or approval by the FDA. By integrating morphological analyses and AI using pathologically relevant cells and disease-relevant stimuli, we identified several compounds and target classes that are potentially able to treat intestinal fibrosis. This phenotypic screening platform offers significant improvements over conventional methods for identifying a wide range of drug targets.

Specific epigenetic regulators serve as potential therapeutic targets in idiopathic pulmonary fibrosis

Idiopathic pulmonary fibrosis (IPF), a disorder observed mostly in older human beings, is characterised by chronic and progressive lung scarring leading to an irreversible decline in lung function. This health condition has a dismal prognosis and the currently available drugs only delay but fail to reverse the progression of lung damage. Consequently, it becomes imperative to discover improved therapeutic compounds and their cellular targets to cure IPF. In this regard, a number of recent studies have targeted the epigenetic regulation by histone deacetylases (HDACs) to develop and categorise antifibrotic drugs for lungs. Therefore, this review focuses on how aberrant expression or activity of Classes I, II and III HDACs alter TGF-β signalling to promote events such as epithelial-mesenchymal transition, differentiation of activated fibroblasts into myofibroblasts, and excess deposition of the extracellular matrix to propel lung fibrosis. Further, this study describes how certain chemical compounds or dietary changes modulate dysregulated HDACs to attenuate five faulty TGF-β-dependent profibrotic processes, both in animal models and cell lines replicating IPF, thereby identifying promising means to treat this lung disorder.

Anti-Inflammatory and Anti-Fibrotic Effect of Immortalized Mesenchymal-Stem-Cell-Derived Conditioned Medium on Human Lung Myofibroblasts and Epithelial Cells

Idiopathic pulmonary fibrosis (IPF) is caused by progressive lung tissue impairment due to extended chronic fibrosis, and it has no known effective treatment. The use of conditioned media (CM) from an immortalized human adipose mesenchymal stem cell line could be a promising therapeutic strategy, as it can reduce both fibrotic and inflammatory responses. We aimed to investigate the anti-inflammatory and anti-fibrotic effect of CM on human pulmonary subepithelial myofibroblasts (hPSM) and on A549 pulmonary epithelial cells, treated with pro-inflammatory or pro-fibrotic mediators. CM inhibited the proinflammatory cytokine-induced mRNA and protein production of various chemokines in both hPSMs and A549 cells. It also downregulated the mRNA expression of IL-1α, but upregulated IL-1β and IL-6 mRNA production in both cell types. CM downregulated the pro-fibrotic-induced mRNA expression of collagen Type III and the migration rate of hPSMs, but upregulated fibronectin mRNA production and the total protein collagen secretion. CM's direct effect on the chemotaxis and cell recruitment of immune-associated cells, and its indirect effect on fibrosis through the significant decrease in the migration capacity of hPSMs, makes it a plausible candidate for further development towards a therapeutic treatment for IPF.

Identification and Validation of Aging-Related Genes in Idiopathic Pulmonary Fibrosis

Aging plays a significant role in the occurrence and development of idiopathic pulmonary fibrosis (IPF). In this study, we aimed to identify and verify potential aging-associated genes involved in IPF using bioinformatic analysis. The mRNA expression profile dataset GSE150910 available in the Gene Expression Omnibus (GEO) database and R software were used to identify the differentially expressed aging-related genes involved in IPF. Hub gene expression was validated by other GEO datasets. Gene ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were performed on differentially expressed aging-related genes. Subsequently, aging-related genes were further screened using three techniques (least absolute shrinkage and selection operator (LASSO) regression, support vector machine, and random forest), and the receiver operating characteristic curves were plotted based on screening results. Finally, real-time quantitative polymerase chain reaction (qRT-PCR) was performed to verify the RNA expression of the six differentially expressed aging-related genes using the blood samples of patients with IPF and healthy individuals. Sixteen differentially expressed aging-related genes were detected, of which the expression of 12 were upregulated and four were downregulated. GO and KEGG enrichment analyses indicated the presence of several enriched terms related to senescence and apoptotic mitochondrial changes. Further screening by LASSO regression, support vector machine, and random forest identified six genes (IGF1, RET, IGFBP2, CDKN2A, JUN, and TFAP2A) that could serve as potential diagnostic biomarkers for IPF. Furthermore, qRT-PCR analysis indicated that among the above-mentioned six aging-related genes, only the expression levels of IGF1, RET, and IGFBP2 in patients with IPF and healthy individuals were consistent with the results of bioinformatic analysis. In conclusion, bioinformatics analysis identified 16 potential aging-related genes associated with IPF, and clinical sample validation suggested that among these, IGF1, RET, and IGFBP2 might play a role in the incidence and prognosis of IPF. Our findings may help understand the pathogenesis of IPF.