Gene correlation network analysis to identify regulatory factors in idiopathic pulmonary fibrosis

Long-term exposure to diesel exhaust particles induces concordant changes in DNA methylation and transcriptome in human adenocarcinoma alveolar basal epithelial cells

BACKGROUND

Diesel exhaust particles (DEP), which contain hazardous compounds, are emitted during the combustion of diesel. As approximately one-third of the vehicles worldwide use diesel, there are growing concerns about the risks posed by DEP to human health. Long-term exposure to DEP is associated with airway hyperresponsiveness, pulmonary fibrosis, and inflammation; however, the molecular mechanisms behind the effects of DEP on the respiratory tract are poorly understood. Such mechanisms can be addressed by examining transcriptional and DNA methylation changes. Although several studies have focused on the effects of short-term DEP exposure on gene expression, research on the transcriptional effects and genome-wide DNA methylation changes caused by long-term DEP exposure is lacking. Hence, in this study, we investigated transcriptional and DNA methylation changes in human adenocarcinoma alveolar basal epithelial A549 cells caused by prolonged exposure to DEP and determined whether these changes are concordant.

RESULTS

DNA methylation analysis using the Illumina Infinium MethylationEPIC BeadChips showed that the methylation levels of DEP-affected CpG sites in A549 cells changed in a dose-dependent manner; the extent of change increased with increasing dose reaching the statistical significance only in samples exposed to 30 µg/ml DEP. Four-week exposure to 30 µg/ml of DEP significantly induced DNA hypomethylation at 24,464 CpG sites, which were significantly enriched for DNase hypersensitive sites, genomic regions marked by H3K4me1 and H3K27ac, and several transcription factor binding sites. In contrast, 9,436 CpG sites with increased DNA methylation levels were significantly overrepresented in genomic regions marked by H3K27me3 as well as H3K4me1 and H3K27ac. In parallel, gene expression profiling by RNA sequencing demonstrated that long-term exposure to DEP altered the expression levels of 2,410 genes, enriching 16 gene sets including Xenobiotic metabolism, Inflammatory response, and Senescence. In silico analysis revealed that the expression levels of 854 genes correlated with the methylation levels of the DEP-affected cis-CpG sites.

CONCLUSIONS

To our knowledge, this is the first report of genome-wide transcriptional and DNA methylation changes and their associations in A549 cells following long-term exposure to DEP.

Identifying a survival-associated cell type based on multi-level transcriptome analysis in idiopathic pulmonary fibrosis

BACKGROUND

Idiopathic pulmonary fibrosis (IPF) is a progressive disease with a five-year survival rate of less than 40%. There is significant variability in survival time among IPF patients, but the underlying mechanisms for this are not clear yet.

METHODS AND RESULTS

We collected single-cell RNA sequence data of 13,223 epithelial cells taken from 32 IPF patients and bulk RNA sequence data from 456 IPF patients in GEO. Based on unsupervised clustering analysis at the single-cell level and deconvolution algorithm at bulk RNA sequence data, we discovered a special alveolar type 2 cell subtype characterized by high expression of CCL20 (referred to as ATII-CCL20), and found that IPF patients with a higher proportion of ATII-CCL20 had worse prognoses. Furthermore, we uncovered the upregulation of immune cell infiltration and metabolic functions in IPF patients with a higher proportion of ATII-CCL20. Finally, the comprehensive decision tree and nomogram were constructed to optimize the risk stratification of IPF patients and provide a reference for accurate prognosis evaluation.

CONCLUSIONS

Our study by integrating single-cell and bulk RNA sequence data from IPF patients identified a special subtype of ATII cells, ATII-CCL20, which was found to be a risk cell subtype associated with poor prognosis in IPF patients. More importantly, the ATII-CCL20 cell subtype was linked with metabolic functions and immune infiltration.

Identification of core genes affecting IMF deposition in bovine

Intramuscular fat (IMF) content is an important economic factor in beef production. However, knowledge on the key factors controlling bovine IMF is limited. In this study, using weighted gene co-expression network analysis (WGCNA), nine modules were identified and the number of transcripts in these modules ranged from 36 to 3191. Two modules were found to be significantly associated with fat deposition and three genes (TCAP, MYH7, and TNNC1) were further identified by Protein-protein interaction (PPI), which may be the hub genes regulating bovine IMF deposition. In addition, considering the genetic variation, the PCK1 gene was found by functional enrichment analysis of overlapping genes, which was previously reported to be involved in IMF deposition. We noted that the core promoter region of buffalo PCK1 binds to transcription factors involved in lipid metabolism while cattle PCK1 binds transcription factors involved in muscle development. The results suggest that PCK1 participated in IMF deposition of buffalo and cattle in different ways. In summary, gene expression networks and new candidate genes associated with IMF deposition identified in this study. This would lay the foundation for further research into the molecular regulatory mechanisms underlying bovine IMF deposition.

Iron accumulation drives fibrosis, senescence and the senescence-associated secretory phenotype

Fibrogenesis is part of a normal protective response to tissue injury that can become irreversible and progressive, leading to fatal diseases. Senescent cells are a main driver of fibrotic diseases through their secretome, known as senescence-associated secretory phenotype (SASP). Here, we report that cellular senescence, and multiple types of fibrotic diseases in mice and humans are characterized by the accumulation of iron. We show that vascular and hemolytic injuries are efficient in triggering iron accumulation, which in turn can cause senescence and promote fibrosis. Notably, we find that senescent cells persistently accumulate iron, even when the surge of extracellular iron has subdued. Indeed, under normal conditions of extracellular iron, cells exposed to different types of senescence-inducing insults accumulate abundant ferritin-bound iron, mostly within lysosomes, and present high levels of labile iron, which fuels the generation of reactive oxygen species and the SASP. Finally, we demonstrate that detection of iron by magnetic resonance imaging might allow non-invasive assessment of fibrotic burden in the kidneys of mice and in patients with renal fibrosis. Our findings suggest that iron accumulation plays a central role in senescence and fibrosis, even when the initiating events may be independent of iron, and identify iron metabolism as a potential therapeutic target for senescence-associated diseases.

Wild-type S100A3 and S100A13 restore calcium homeostasis and mitigate mitochondrial dysregulation in pulmonary fibrosis patient-derived cells

Patients with digenic S100A3 and S100A13 mutations exhibited an atypical and progressive interstitial pulmonary fibrosis, with impaired intracellular calcium homeostasis and mitochondrial dysfunction. Here we provide direct evidence of a causative effect of the mutation on receptor mediated calcium signaling and calcium store responses in control cells transfected with mutant S100A3 and mutant S100A13. We demonstrate that the mutations lead to increased mitochondrial mass and hyperpolarization, both of which were reversed by transfecting patient-derived cells with the wild type S100A3 and S100A13, or extracellular treatment with the recombinant proteins. In addition, we demonstrate increased secretion of inflammatory mediators in patient-derived cells and in control cells transfected with the mutant-encoding constructs. These findings indicate that treatment of patients' cells with recombinant S100A3 and S100A13 proteins is sufficient to normalize most of cellular responses, and may therefore suggest the use of these recombinant proteins in the treatment of this devastating disease.

Interplay of chronic obstructive pulmonary disease and colorectal cancer development: unravelling the mediating role of fatty acids through a comprehensive multi-omics analysis

BACKGROUND

Chronic obstructive pulmonary disease (COPD) patients often exhibit gastrointestinal symptoms, A potential association between COPD and Colorectal Cancer (CRC) has been indicated, warranting further examination.

METHODS

In this study, we collected COPD and CRC data from the National Health and Nutrition Examination Survey, genome-wide association studies, and RNA sequence for a comprehensive analysis. We used weighted logistic regression to explore the association between COPD and CRC incidence risk. Mendelian randomization analysis was performed to assess the causal relationship between COPD and CRC, and cross-phenotype meta-analysis was conducted to pinpoint crucial loci. Multivariable mendelian randomization was used to uncover mediating factors connecting the two diseases. Our results were validated using both NHANES and GEO databases.

RESULTS

In our analysis of the NHANES dataset, we identified COPD as a significant contributing factor to CRC development. MR analysis revealed that COPD increased the risk of CRC onset and progression (OR: 1.16, 95% CI 1.01-1.36). Cross-phenotype meta-analysis identified four critical genes associated with both CRC and COPD. Multivariable Mendelian randomization suggested body fat percentage, omega-3, omega-6, and the omega-3 to omega-6 ratio as potential mediating factors for both diseases, a finding consistent with the NHANES dataset. Further, the interrelation between fatty acid-related modules in COPD and CRC was demonstrated via weighted gene co-expression network analysis and Kyoto Encyclopedia of Genes and Genomes enrichment results using RNA expression data.

CONCLUSIONS

This study provides novel insights into the interplay between COPD and CRC, highlighting the potential impact of COPD on the development of CRC. The identification of shared genes and mediating factors related to fatty acid metabolism deepens our understanding of the underlying mechanisms connecting these two diseases.

Construction of a TFs-miRNA-mRNA network related to idiopathic pulmonary fibrosis

Background

The transcription factors (TFs)-microRNA (miRNA)-messenger RNA (mRNA) network plays an important role in a variety of diseases. However, the relationship between the TFs-miRNA-mRNA network and idiopathic pulmonary fibrosis (IPF) remains unclear.

Methods

The GSE110147 and GSE53845 datasets from the Gene Expression Omnibus (GEO) database were used to process differentially expressed genes (DEGs) analysis, gene set enrichment analysis (GSEA), weighted gene co-expression network analysis (WGCNA), as well as Gene ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses. The GSE13316 dataset was used to perform differentially expressed miRNAs (DEMs) analysis and TFs prediction. Finally, a TFs-miRNA-mRNA network related to IPF was constructed, and its function was evaluated by Gene Ontology (GO) and KEGG analyses. Also, 19 TFs in the network were verified by quantitative real time polymerase chain reaction (qRT-PCR).

Results

Through our analysis, 53 DEMs and 2,630 DEGs were screened. The GSEA results suggested these genes were mainly related to protein digestion and absorption. The WGCNA results showed that these DEGs were divided into eight modules, and the GO and KEGG analyses results of blue module genes showed that these 86 blue module genes were mainly enriched in cilium assembly and cilium organization. Moreover, a TFs-miRNA-mRNA network comprising 25 TFs, 11 miRNAs, and 60 mRNAs was constructed. Ultimately, the functional enrichment analysis showed that the TFs-miRNA-mRNA network was mainly related to the cell cycle and the phosphatidylinositol 3 kinase-protein kinase B (PI3K-Akt) signaling pathway. Furthermore, experimental verification of the TFs showed that ARNTL, TRIM28, EZH2, BCOR, and ASXL1 were sufficiently up-regulated in the transforming growth factor (TGF)-β1 treatment groups, while BCL6, BHLHE40, FOXA1, and EGR1 were significantly down-regulated.

Conclusions

The novel TFs-miRNA-mRNA network that we constructed could provide new insights into the underlying molecular mechanisms of IPF. ARNTL, TRIM28, EZH2, BCOR, ASXL1, BCL6, BHLHE40, FOXA1, and EGR1 may play important roles in IPF and become effective biomarkers for diagnosis and treatment.

Cellular and Molecular Mechanisms in Idiopathic Pulmonary Fibrosis

The respiratory system is a well-organized multicellular organ, and disruption of cellular homeostasis or abnormal tissue repair caused by genetic deficiency and exposure to risk factors lead to life-threatening pulmonary disease including idiopathic pulmonary fibrosis (IPF). Although there is no clear etiology as the name reflected, its pathological progress is closely related to uncoordinated cellular and molecular signals. Here, we review the advances in our understanding of the role of lung tissue cells in IPF pathology including epithelial cells, mesenchymal stem cells, fibroblasts, immune cells, and endothelial cells. These advances summarize the role of various cell components and signaling pathways in the pathogenesis of idiopathic pulmonary fibrosis, which is helpful to further study the pathological mechanism of the disease, provide new opportunities for disease prevention and treatment, and is expected to improve the survival rate and quality of life of patients.

Identification of key pathways and genes in the progression of silicosis based on WGCNA

Silicosis, induced by inhaling silica particles in workplaces, is one of the most common occupational diseases. The prognosis of silicosis and its consequent fibrosis is extremely poor due to limited treatment modalities and lack of understanding of the disease mechanisms. In this study, a Wistar rat model for silicosis fibrosis was established by intratracheal instillation of silica (0, 50, 100 and 200 mg/mL, 1 mL) with the evidence of Hematoxylin and Eosin (HE) and Masson staining and the expressions of inflammatory and fibrotic proteins of rats' lung tissues. RNA of lung tissues of rats exposed to 200 mg/mL silica particles and normal saline for 14 d and 28 d was extracted and sequenced to detect differentially expressed genes (DEGs) and to identify silicosis fibrosis-associated modules and hub genes by Weighted gene co-expression network analysis (WGCNA). Predictions of gene functions and signaling pathways were conducted using Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. In this study, it has been demonstrated the promising role of the Hippo signaling pathway in silicosis fibrosis, which will be conducive to elucidating the specific mechanism of pulmonary fibrosis induced by silica and to determining molecular initiating event (MIE) and adverse outcome pathway (AOP) of silicosis fibrosis.

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

Idiopathic pulmonary fibrosis (IPF) is a severe fibrotic lung disease characterized by irreversible scarring of the lung parenchyma leading to dyspnea, progressive decline in lung function, and respiratory failure. We analyzed lung transcriptomic data from independent IPF cohorts using weighted gene co-expression network analysis (WGCNA) to identify gene modules based on their preservation status in these cohorts. The consensus gene modules were characterized by leveraging existing clinical and molecular data such as lung function, biological processes, pathways, and lung cell types. From a total of 32 consensus gene modules identified, two modules were found to be significantly correlated with the disease, lung function, and preserved in other IPF datasets. The upregulated gene module was enriched for extracellular matrix, collagen metabolic process, and BMP signaling while the downregulated module consisted of genes associated with tube morphogenesis, blood vessel development, and cell migration. Using a combination of connectivity-based and trait-based significance measures, we identified and prioritized 103 "hub" genes (including 25 secretory candidate biomarkers) by their similarity to known IPF genetic markers. Our validation studies demonstrate the dysregulated expression of CRABP2, a retinol-binding protein, in multiple lung cells of IPF, and its correlation with the decline in lung function.

eNAMPT Is a Novel Damage-associated Molecular Pattern Protein That Contributes to the Severity of Radiation-induced Lung Fibrosis

The paucity of therapeutic strategies to reduce the severity of radiation-induced lung fibrosis (RILF), a life-threatening complication of intended or accidental ionizing radiation exposure, is a serious unmet need. We evaluated the contribution of eNAMPT (extracellular nicotinamide phosphoribosyltransferase), a damage-associated molecular pattern (DAMP) protein and TLR4 (Toll-like receptor 4) ligand, to the severity of whole-thorax lung irradiation (WTLI)-induced RILF. Wild-type (WT) and Nampt+/- heterozygous C57BL6 mice and nonhuman primates (NHPs, Macaca mulatta) were exposed to a single WTLI dose (9.8 or 10.7 Gy for NHPs, 20 Gy for mice). WT mice received IgG1 (control) or an eNAMPT-neutralizing polyclonal or monoclonal antibody (mAb) intraperitoneally 4 hours after WTLI and weekly thereafter. At 8-12 weeks after WTLI, NAMPT expression was assessed by immunohistochemistry, biochemistry, and plasma biomarker studies. RILF severity was determined by BAL protein/cells, hematoxylin and eosin, and trichrome blue staining and soluble collagen assays. RNA sequencing and bioinformatic analyses identified differentially expressed lung tissue genes/pathways. NAMPT lung tissue expression was increased in both WTLI-exposed WT mice and NHPs. Nampt+/- mice and eNAMPT polyclonal antibody/mAb-treated mice exhibited significantly attenuated WTLI-mediated lung fibrosis with reduced: 1) NAMPT and trichrome blue staining; 2) dysregulated lung tissue expression of smooth muscle actin, p-SMAD2/p-SMAD1/5/9, TGF-β, TSP1 (thrombospondin-1), NOX4, IL-1β, and NRF2; 3) plasma eNAMPT and IL-1β concentrations; and 4) soluble collagen. Multiple WTLI-induced dysregulated differentially expressed lung tissue genes/pathways with known tissue fibrosis involvement were each rectified in mice receiving eNAMPT mAbs.The eNAMPT/TLR4 inflammatory network is essentially involved in radiation pathobiology, with eNAMPT neutralization an effective therapeutic strategy to reduce RILF severity.

Stemness Subtypes and Scoring System Predict Prognosis and Efficacy of Immunotherapy in Soft Tissue Sarcoma

Tumor stemness has been reported to play important roles in cancers. However, a comprehensive analysis of tumor stemness remains to be performed to investigate the specific mechanisms and practical values of stemness in soft tissue sarcomas (STS). Here, we applied machine learning to muti-omic data of patients from TCGA-SARC and GSE21050 cohorts to reveal important roles of stemness in STS. We demonstrated limited roles of existing mRNAsi in clinical application. Therefore, based on stemness-related signatures (SRSs), we identified three stemness subtypes with distinct stemness, immune, and metabolic characteristics using consensus clustering. The low-stemness subtype had better prognosis, activated innate and adaptive immunity (e.g., infiltrating B, DC, Th1, CD8+ T, activated NK, gamma delta T cells, and M1 macrophages), more enrichment of metabolic pathways, more sites with higher methylation level, higher gene mutations, CNA burdens, and immunogenicity indicators. Furthermore, the 16 SRS-based stemness prognostic index (SPi) was developed, and we found that low-SPi patients with low stemness had better prognosis and other characteristics similar to those in the low-stemness subtype. Besides, low-stemness subtype and low-SPi patients could benefit from immunotherapy. The predictive value of SPi in immunotherapy was more accurate after the addition of MSI into SPi. MSI^lowSPi^low patients might be more sensitive to immunotherapy. In conclusion, we highlighted mechanisms and practical values of the stemness in STS. We also recommended the combination of MSI and SPi which is a promising tool to predict prognosis and achieve precise treatments of immunotherapy in STS.

The Prognostic Value of Integrated Analysis of Inflammation and Hypoxia-Related Genes in Idiopathic Pulmonary Fibrosis

Currently, the aetiology and pathogenesis of idiopathic pulmonary fibrosis (IPF) are still largely unclear. Moreover, patients with IPF exhibit a considerable difference in clinical presentation, treatment, and prognosis. Optimal biomarkers or models for IPF prognosis are lacking. Therefore, this study quantified the levels of various hallmarks using a single-sample gene set enrichment analysis algorithm. The hazard ration was calculated using Univariate Cox regression analysis based on the transcriptomic profile of bronchoalveolar lavage cells and clinical survival information. Afterwards, weighted Gene Co-expression Network Analysis was performed to construct a network between gene expression, inflammation response, and hypoxia. Subsequently, univariate Cox, random forest, and multivariate Cox regressions were applied to develop a robust inflammation and hypoxia-related gene signature for predicting clinical outcomes in patients with IPF. Furthermore, a nomogram was constructed to calculate risk assessment. The inflammation response and hypoxia were identified as latent risk factors for patients with IPF. Five genes, including HS3ST1, WFDC2, SPP1, TFPI, and CDC42EP2, were identified that formed the inflammation-hypoxia-related gene signature. Kaplan-Meier plotter showed that the patients with high-risk scores had a worse prognosis than those with low-risk scores in training and validation cohorts. The time-dependent concordance index and the receiver operating characteristic analysis revealed that the risk model could accurately predict the clinical outcome of patients with IPF. Therefore, this study contributes to elucidating the role of inflammation and hypoxia in IPF, which can aid in assessing individual prognosis and personalised treatment decisions.

Mucus Hypersecretion and Ciliary Impairment in Conducting Airway Contribute to Alveolar Mucus Plugging in Idiopathic Pulmonary Fibrosis

Idiopathic pulmonary fibrosis (IPF) is a chronic lung disease attributed to the complex interplay of genetic and environmental risks. The muco-ciliary clearance (MCC) system plays a critical role in maintaining the conduit for air to and from the alveoli, but it remains poorly understood whether the MCC abnormalities in conducting airway are involved in IPF pathogenesis. In this study, we obtained the surgically resected bronchi and peripheral lung tissues from 31 IPF patients and 39 control subjects, and we sought to explore the morphologic characteristics of MCC in conducting airway by using immunostaining and scanning and transmission electron microscopy. In the submucosal regions of the bronchi, we found that the areas of mucus glands (MUC5B⁺) were significantly larger in IPF patients as compared with control subjects (p < 0.05). In the surface epithelium of three airway regions (bronchi, proximal bronchioles, and distal bronchioles), increased MUC5B and MUC5AC expression of secretory cells, decreased number of ciliated cells, and increased ciliary length were observed in IPF patients than control subjects (all p < 0.05). In addition, the mRNA expression levels of MUC5B were up-regulated in both the bronchi and peripheral lung of IPF patients than those of control subjects (p < 0.05), accompanied with 93.55% IPF subjects who had obvious MUC5B⁺ mucus plugs in alveolar regions. No MUC5B rs35705950 single-nucleotide polymorphism allele was detected in both IPF patients and control subjects. Our study shows that mucus hypersecretion and ciliary impairment in conducting airway are major causes of mucus plugs in alveolar regions and may be closely related to the alveolar injuries in IPF patients.

Matrix Metalloproteinase 7 Expression and Apical Epithelial Defects in Atp8b1 Mutant Mouse Model of Pulmonary Fibrosis

Abnormalities in airway epithelia and lung parenchyma are found in Atp8b1 mutant mice, which develop pulmonary fibrosis after hyperoxic insult. Microarray and ingenuity pathway analysis (IPA) show numerous transcripts involved in ciliogenesis are downregulated in 14-month (14 M) -old Atp8b1 mouse lung compared with wild-type C57BL/6. Lung epithelium of Atp8b1 mice demonstrate apical abnormalities of ciliated and club cells in the bronchial epithelium on transmission electron microscopy (TEM). Matrix metalloproteinase 7 (MMP7) regulates of ciliogenesis and is a biomarker for idiopathic pulmonary fibrosis (IPF) in humans. Mmp7 transcript and protein expression are significantly upregulated in 14 M Atp8b1 mutant mouse lung. MMP7 expression is also increased in bronchoalveolar lavage fluid (BAL). Immunohistochemistry is localized MMP7 to bronchial epithelial cells in the Atp8b1 mutant. In conclusion, MMP7 is upregulated in the aged Atp8b1 mouse model, which displays abnormal ciliated cell and club cell morphology. This mouse model can facilitate the exploration of the role of MMP7 in epithelial integrity and ciliogenesis in IPF. The Atp8b1 mutant mouse is proposed as a model for IPF.

Urine-derived exosomes from individuals with IPF carry pro-fibrotic cargo

BACKGROUND

MicroRNAs (miRNA) and other components contained in extracellular vesicles may reflect the presence of a disease. Lung tissue, sputum, and sera of individuals with idiopathic pulmonary fibrosis (IPF) show alterations in miRNA expression. We designed this study to test whether urine and/or tissue derived exosomal miRNAs from individuals with IPF carry cargo that can promote fibrosis.

METHODS

Exosomes were isolated from urine (U-IPFexo), lung tissue myofibroblasts (MF-IPFexo), serum from individuals with IPF (n=16) and age/sex-matched controls without lung disease (n=10). We analyzed microRNA expression of isolated exosomes and their in vivo bio-distribution. We investigated the effect on ex vivo skin wound healing and in in vivo mouse lung models.

RESULTS

U-IPFexo or MF-IPFexo expressed miR-let-7d, miR-29a-5p, miR-181b-3p and miR-199a-3p consistent with previous reports of miRNA expression obtained from lung tissue/sera from patients with IPF. In vivo bio-distribution experiments detected bioluminescent exosomes in the lung of normal C57Bl6 mice within 5 min after intravenous infusion, followed by distribution to other organs irrespective of exosome source. Exosomes labeled with gold nanoparticles and imaged by transmission electron microscopy were visualized in alveolar epithelial type I and type II cells. Treatment of human and mouse lung punches obtained from control, non-fibrotic lungs with either U-IPFexo or MF-IPFexo produced a fibrotic phenotype. A fibrotic phenotype was also induced in a human ex vivo skin model and in in vivo lung models.

CONCLUSIONS

Our results provide evidence of a systemic feature of IPF whereby exosomes contain pro-fibrotic miRNAs when obtained from a fibrotic source and interfere with response to tissue injury as measured in skin and lung models.

FUNDING

This work was supported in part by Lester and Sue Smith Foundation and The Samrick Family Foundation and NIH grants R21 AG060338 (SE and MKG), U01 DK119085 (IP, RS, MTC).

Lung Microenvironments and Disease Progression in Fibrotic Hypersensitivity Pneumonitis

Rationale: Fibrotic hypersensitivity pneumonitis (fHP) is an interstitial lung disease caused by sensitization to an inhaled allergen. Objectives: To identify the molecular determinants associated with progression of fibrosis. Methods: Nine fHP explant lungs and six unused donor lungs (as controls) were systematically sampled (4 samples/lung). According to microcomputed tomography measures, fHP cores were clustered into mild, moderate, and severe fibrosis groups. Gene expression profiles were assessed using weighted gene co-expression network analysis, xCell, gene ontology, and structure enrichment analysis. Gene expression of the prevailing molecular traits was also compared with idiopathic pulmonary fibrosis (IPF). The explant lung findings were evaluated in separate clinical fHP cohorts using tissue, BAL samples, and computed tomography scans. Measurements and Main Results: We found six molecular traits that associated with differential lung involvement. In fHP, extracellular matrix and antigen presentation/sensitization transcriptomic signatures characterized lung zones with only mild structural and histological changes, whereas signatures involved in honeycombing and B cells dominated the transcriptome in the most severely affected lung zones. With increasing disease severity, endothelial function was progressively lost, and progressive disruption in normal cellular homeostatic processes emerged. All six were also found in IPF, with largely similar associations with disease microenvironments. The molecular traits correlated with in vivo disease behavior in a separate clinical fHP cohort. Conclusions: We identified six molecular traits that characterize the morphological progression of fHP and associate with in vivo clinical behavior. Comparing IPF with fHP, the transcriptome landscape was determined considerably by local disease extent rather than by diagnosis alone.

A novel prognostic signature for idiopathic pulmonary fibrosis based on five-immune-related genes

Background

Idiopathic pulmonary fibrosis (IPF) is a highly fatal lung disease of unknown etiology with a median survival after diagnosis of only 2–3 years. Its poor prognosis is due to the limited therapy options available as well as the lack of effective prognostic indicators. This study aimed to construct a novel prognostic signature for IPF to assist in the personalized management of IPF patients during treatment.

Methods

Differentially-expressed genes (DEGs) in IPF patients versus healthy individuals were analyzed using the “limma” package of R software. Immune-related genes (IRGs) were obtained from the ImmPort database. Univariate Cox regression analysis was adopted to screen significantly prognostic IRGs for IPF patients. Multiple Cox regression analysis was used to identify optimal prognostic IRGs and construct a prognostic signature.

Results

Compared with healthy individuals, there were a total of 52 prognosis-related DEGs in the bronchoalveolar lavage (BAL) samples of IPF patients, of which 37 genes were identified as IRGs. Of these, five genes (CXCL14, SLC40A1, RNASE3, CCR3, and RORA) were significantly associated with overall survival (OS) in IPF patients, and were utilized for establishment of the prognostic signature. IPF patients were divided into high- and low-risk groups based on the prognostic signature. Marked differences in the OS probability were observed between high- and low-risk IPF patients. The area under curves (AUCs) of the receiver operating characteristic (ROC) curve for the prognostic signature in the training and validation cohorts were 0.858 and 0.837, respectively. The expression levels between RNASE3 and SLC40A1 (P<0.01, r=0.394), between RORA and CXCL14 (P<0.01, r=−0.355), between CCR3 and CXCL14 (P<0.01, r=0.258), as well as between RNASE3 and CCR3 (P<0.01, r=0.293) were significantly correlated.

Conclusions

We developed a validated and reproducible IRG-based prognostic signature that should be helpful in the personalized management of patients with IPF, providing new insights into the relationship between the immune system and IPF.

Comprehensive Analysis of Regulatory Factors and Immune-Associated Patterns to Decipher Common and BRCA1/2 Mutation-Type-Specific Critical Regulation in Breast Cancer

Background: BRCA1/2 mutations are closely related to high lifetime risk of breast cancer (BC). The objective of this study was to identify the genes, regulators, and immune-associated patterns underlying disease pathology in BC with BRCA1/2 somatic mutations and their associations with clinical traits. Methods: RNA sequencing data and clinical information from The Cancer Genome Atlas (TCGA; N = 36 BRCA1-mutant BC; N = 49 BRCA2-mutant BC; and N = 117 BRCA1/2-wild-type BC samples) were used for discovery, which included consensus network analysis, function enrichment, and analysis of hub genes; other TCGA data (N = 117 triple-negative BC) and two Gene Expression Omnibus database expression profiles were used as validation cohorts. Results: Consensus network analysis helped to identify specific co-expressed modules that showed positive correlations with tumor stage, number of positive lymph nodes, and margin status in BRCA1/2-mutant BC but lacking correlations in BRCA1/2-wild-type BC. Functional enrichment suggested potential mechanisms in BRCA1/2 carriers that could regulate the cell cycle, immune response, cellular metabolic processes, and cell migration, via enriched pathways including p53 and JAK-STAT signaling. Consensus network analysis identified the specific and common carcinogenic mechanisms involving BRCA mutations. Regulators cross-linking these modules include E2F or IRF transcription factor family, associated with cell cycle or immune response regulation module, respectively. Eight hub genes, including ISG15, BUB1, and TTK, were upregulated in several BRCA1/2-mutant BC datasets and showed prognostic value in BC. Furthermore, their genetic expression was related to higher levels of immune infiltration in BRCA1/2-mutant BC, which manifested as recruitment of T helper cells (Th1 cells), follicular helper T cells, and regulatory T cells, and T cell exhaustion. Moreover, important indicators for evaluation of BC immunotherapy, tumor mutational burden and neoantigen load also positively correlated with expression of some hub genes. Conclusion: We constructed a BRCA1/2 mutation-type-specific co-expressed gene network with related transcription factors and immune-associated patterns that could regulate and influence tumor metastasis and immune microenvironment, providing novel insights into the pathological process of this disease and the corresponding BRCA mutations.

Fibromine is a multi-omics database and mining tool for target discovery in pulmonary fibrosis

Idiopathic pulmonary fibrosis is a lethal lung fibroproliferative disease with limited therapeutic options. Differential expression profiling of affected sites has been instrumental for involved pathogenetic mechanisms dissection and therapeutic targets discovery. However, there have been limited efforts to comparatively analyse/mine the numerous related publicly available datasets, to fully exploit their potential on the validation/creation of novel research hypotheses. In this context and towards that goal, we present Fibromine, an integrated database and exploration environment comprising of consistently re-analysed, manually curated transcriptomic and proteomic pulmonary fibrosis datasets covering a wide range of experimental designs in both patients and animal models. Fibromine can be accessed via an R Shiny application (http://www.fibromine.com/Fibromine) which offers dynamic data exploration and real-time integration functionalities. Moreover, we introduce a novel benchmarking system based on transcriptomic datasets underlying characteristics, resulting to dataset accreditation aiming to aid the user on dataset selection. Cell specificity of gene expression can be visualised and/or explored in several scRNA-seq datasets, in an effort to link legacy data with this cutting-edge methodology and paving the way to their integration. Several use case examples are presented, that, importantly, can be reproduced on-the-fly by a non-specialist user, the primary target and potential user of this endeavour.

Identification of the Molecular Subgroups in Idiopathic Pulmonary Fibrosis by Gene Expression Profiles

Background

Idiopathic Pulmonary Fibrosis (IPF) is one of the most common idiopathic interstitial pneumonia, which can occur all over the world. The median survival time of patients is about 3-5 years, and the mortality is relatively high.

Objective

To reveal the potential molecular characteristics of IPF and deepen the understanding of the molecular mechanism of IPF. In order to provide some guidance for the clinical treatment, new drug development, and prognosis judgment of IPF. Although the preliminary conclusion of this study has certain guiding significance for the treatment of IPF and so on, it needs more accurate analytical approaches and large sample clinical trials to verify.

Methods

220 patients with IPF were divided into different subgroups according to the gene expression profiles, which were obtained from the Gene Expression Omnibus (GEO) database. In addition, these subgroups present different expression forms and clinical features. Therefore, weighted gene coexpression analysis (WGCNA) was used to seek the differences between subtypes. And six subgroup-specific WGCNA modules were identified.

Results

Combined with the characteristics of WGCNA and KEGG enrichment modules, the autophagic pathway was only upregulated in subgroup I and enriched significantly. The differentiation pathways of Th1 and Th2 cells were only upregulated and enriched in subgroup II. At the same time, combined with clinical information, IPF patients in subgroup II were older and more serious, which may be closely related to the differentiation of Th1 and Th2 cells. In contrast, the neuroactive ligand-receptor interaction pathway and Ca⁺ signaling pathway were significantly upregulated and enriched in subgroup III. Although there was no significant difference in prognosis between subgroup I and subgroup III, their intrinsic biological characteristics were very different. These results suggest that the subtypes may represent risk factors of age and intrinsic biological characteristics and may also partly reflect the severity of the disease.

Conclusion

In conclusion, current studies have improved our understanding of IPF-related molecular mechanisms. At the same time, because the results show that patients from different subgroups may have their own unique gene expression patterns, it reminds us that patients in each subgroup should receive more personalized treatment.

Associations between mobility patterns and COVID-19 deaths during the pandemic: A network structure and rank propagation modelling approach

BACKGROUND

From February 2020, both urban and rural Ireland witnessed the rapid proliferation of the COVID-19 disease throughout its counties. During this period, the national COVID-19 responses included stay-at-home directives issued by the state, subject to varying levels of enforcement.

METHODS

In this paper, we present a new method to assess and rank the causes of Ireland COVID-19 deaths as it relates to mobility activities within each county provided by Google while taking into consideration the epidemiological confirmed positive cases reported per county. We used a network structure and rank propagation modelling approach using Personalised PageRank to reveal the importance of each mobility category linked to cases and deaths. Then a novel feature-selection method using relative prominent factors finds important features related to each county's death. Finally, we clustered the counties based on features selected with the network results using a customised network clustering algorithm for the research problem.

FINDINGS

Our analysis reveals that the most important mobility trend categories that exhibit the strongest association to COVID-19 cases and deaths include retail and recreation and workplaces. This is the first time a network structure and rank propagation modelling approach has been used to link COVID-19 data to mobility patterns. The infection determinants landscape illustrated by the network results aligns soundly with county socio-economic and demographic features. The novel feature selection and clustering method presented clusters useful to policymakers, managers of the health sector, politicians and even sociologists. Finally, each county has a different impact on the national total.

Comprehensive characterization of pathological stage-related genes of papillary thyroid cancer along with survival prediction

It is crucial to understand the differences across papillary thyroid cancer (PTC) stages, so as to provide a basis for individualized treatments. Here, comprehensive function characterization of PTC stage‐related genes was performed and a new prognostic signature was developed for advanced patients. Two gene modules were confirmed to be closely associated with PTC stages and further six hub genes were identified that yield excellent diagnostic efficiency between tumour and normal tissues. Genetic alteration analysis indicates that they are much conservative since mutations in the DNA of them rarely occur, but changes of DNA methylation on these six genes show that 12 DNA methylation sites are significantly associated with their corresponding genes' expression. Validation data set testing also suggests that these six stage‐related hub genes would be probably potential biomarkers for marking four stages. Subsequently, a 21‐mRNA‐based prognostic risk model was constructed for PTC stage III/IV patients and it could effectively predict the survival of patients with strong prognostic ability. Functional analysis shows that differential expression genes between high‐ and low‐risk patients would promote the progress of PTC to some extent. Moreover, tumour microenvironment (TME) of high‐risk patients may be more conducive to tumour growth by ESTIMATE analysis.

Pathogenesis of fibrosis in interstitial lung disease

PURPOSE OF REVIEW

Pulmonary fibrosis is a chronic and progressive lung disease involving unclear pathological mechanisms. The present review presents and discusses the major and recent advances in our knowledge of the pathogenesis of lung fibrosis.

RECENT FINDINGS

The past months have deepened our understanding on the cellular actors of fibrosis with a better characterization of the abnormal lung epithelial cells observed during lung fibrosis. Better insight has been gained into fibroblast biology and the role of immune cells during fibrosis. Mechanistically, senescence appears as a key driver of the fibrotic process. Extracellular vesicles have been discovered as participating in the impaired cellular cross-talk during fibrosis and deeper understanding has been made on developmental signaling in lung fibrosis.

SUMMARY

This review emphasizes the contribution of different cell types and mechanisms during pulmonary fibrosis, highlights new insights for identification of potential therapeutic strategies, and underlines where future research is needed to answer remaining open questions.

The function of non-coding RNAs in idiopathic pulmonary fibrosis

Non-coding ribonucleic acids (ncRNAs) are a diverse group of RNA molecules that are mostly not translated into proteins after transcription, including long non-coding RNAs (lncRNAs) with longer than 200 nucleotides non-coding transcripts and microRNAs (miRNAs) which are only 18–22 nucleotides. As families of evolutionarily conserved ncRNAs, lncRNAs activate and repress genes via a variety of mechanisms at both transcriptional and translational levels, whereas miRNAs regulate protein-coding gene expression mainly through mRNA silencing. ncRNAs are widely involved in biological functions, such as proliferation, differentiation, migration, angiogenesis, and apoptosis. Idiopathic pulmonary fibrosis (IPF) is a progressive lung disease with a poor prognosis. The etiology of IPF is still unclear. Increasing evidence shows the close correlations between the development of IPF and aberrant expressions of ncRNAs than thought previously. In this study, we provide an overview of ncRNAs participated in pathobiology of IPF, seeking the early diagnosis biomarker and aiming for potential therapeutic applications for IPF.

Construction and validation of a bronchoalveolar lavage cell-associated gene signature for prognosis prediction in idiopathic pulmonary fibrosis

BACKGROUND

Idiopathic pulmonary fibrosis (IPF) is a chronic and progressive interstitial lung disease. It is urgent to identify biomarkers to precisely predict mortality.

METHODS

Gene expression data of bronchoalveolar lavage (BAL) cells and clinical information were downloaded from the Gene Expression Omnibus database. We identified key modules associated with prognosis using weighted gene co-expression network analysis (WGCNA). Then we screened genes with the least absolute shrinkage and selection operator Cox regression. Finally, we constructed a prognostic gene signature using multivariate Cox regression. The risk model was evaluated using the time-dependent receiver operating characteristic (ROC) curve and the concordance index. Additionally, the risk model was validated using an external independent dataset.

RESULTS

Two key modules, strongly associated with inflammation and immune response, were identified by WGCNA. Four genes, including TLR2, CCR2, HTRA1, and SFN, were screened to construct the prognostic model. The patients with a high-risk score had a significantly worse prognosis than patients with a low-risk score. Time-dependent ROC analysis showed that the risk model had a moderate predictive performance for overall survival in the training and external validation datasets.

CONCLUSIONS

Our study provides new insights into the prognostic value of BAL cells in IPF and it may be helpful to assist clinicians in making treatment decisions for the personalized management of IPF.

Anti-fibrotic effects of different sources of MSC in bleomycin-induced lung fibrosis in C57BL6 male mice

BACKGROUND AND OBJECTIVE

IPF is a fatal and debilitating lung disorder increasing in incidence worldwide. To date, two approved treatments only slow disease progression, have multiple side effects and do not provide a cure. MSC have promising therapeutic potential as a cell-based therapy for many lung disorders based on the anti-fibrotic properties of the MSC.

METHODS

Critical questions remain surrounding the optimal source, timing and efficacy of cell-based therapies. The present study examines the most effective sources of MSC. Human MSC were derived from adipose, WJ, chorionic membrane (CSC) and chorionic villi (CVC). MSC were injected into the ageing mouse model of BLM-induced lung fibrosis.

RESULTS

All sources decreased Aschroft and hydroxyproline levels when injected into BLM-treated mice at day 10 with the exception of CSC cells that did not change hydroxyproline levels. There were also decreases in mRNA expression of αv -integrin and TNFα in all sources except CSC. Only ASC- and WJ-derived cells reduced AKT and MMP-2 activation, while Cav-1 was increased by ASC treatment as previously reported. BLM-induced miR dysregulation of miR-29 and miR-199 was restored only by ASC treatment.

CONCLUSION

Our data suggest that sources of MSC may differ in the pathway(s) involved in repair.

Integrated analysis of the molecular mechanisms in idiopathic pulmonary fibrosis

Rationale: Idiopathic pulmonary fibrosis (IPF) is one of the most aggressive forms of idiopathic interstitial pneumonia. Some miRNAs may be associated with IPF and may affect the occurrence and development of IPF in various pathways. Many miRNAs and genes that may be involved in the development of IPF have been discovered using chip and high throughput technologies. Methods: We analyzed one miRNA and four mRNA databases. We identified hub genes and pathways related to IPF using GO, KEGG enrichment analysis, gene set variation analysis (GSVA), PPI network construction, and hub gene analysis. A comprehensive analysis of differentially expressed miRNAs (DEMs), predicted miRNA target genes, and differentially expressed genes (DEGs) led to the creation of a miRNA-mRNA regulatory network in IPF. Results: We found 203 DEGs and 165 DEMs that were associated with IPF. The findings of enrichment analyses showed that these DEGs were mainly involved in antimicrobial humoral response, antimicrobial humoral immune response mediated by antimicrobial peptide, extracellular matrix organization, cell killing, and organ or tissue specific immune response. The VEGFA, CDH5, and WNT3A genes overlapped between hub genes and the miRNA-mRNA regulatory network. The miRNAs including miR-199b-5p, miR-140-5p, miR-199a-5p, miR-125A-5p, and miR-107 that we predicted would regulate the VEGFA, CDH5, and WNT3A genes, which were also associated with IPF or other fibrosis-related diseases. GSVA indicated that metabolic processes of UTP and IMP, immune response, regulation of Th2 cell cytokine production, and positive regulation of NK cell-mediated immunity are associated with the pathogenesis and treatment of IPF. These pathways also interact with VEGFA, CDH5, and WNT3A. Conclusion: These findings provide a new research direction for the diagnosis and treatment of IPF.

Profibrotic epithelial TGF-β1 signaling involves NOX4-mitochondria cross talk and redox-mediated activation of the tyrosine kinase FYN

Idiopathic pulmonary fibrosis (IPF) is characterized by a disturbed redox balance and increased production of reactive oxygen species (ROS), which is believed to contribute to epithelial injury and fibrotic lung scarring. The main pulmonary sources of ROS include mitochondria and NADPH oxidases (NOXs), of which the NOX4 isoform has been implicated in IPF. Non-receptor SRC tyrosine kinases (SFK) are important for cellular homeostasis and are often dysregulated in lung diseases. SFK activation by the profibrotic transforming growth factor-β (TGF-β) is thought to contribute to pulmonary fibrosis, but the relevant SFK isoform and its relationship to NOX4 and/or mitochondrial ROS in the context of profibrotic TGF-β signaling is not known. Here, we demonstrate that TGF-β1 can rapidly activate the SRC kinase FYN in human bronchial epithelial cells, which subsequently induces mitochondrial ROS (mtROS) production, genetic damage shown by the DNA damage marker γH2AX, and increased expression of profibrotic genes. Moreover, TGF-β1-induced activation of FYN involves initial activation of NOX4 and direct cysteine oxidation of FYN, and both FYN and mtROS contribute to TGF-β-induced induction of NOX4. NOX4 expression in lung tissues of IPF patients is positively correlated with disease severity, although FYN expression is down-regulated in IPF and does not correlate with disease severity. Collectively, our findings highlight a critical role for FYN in TGF-β1-induced mtROS production, DNA damage response, and induction of profibrotic genes in bronchial epithelial cells, and suggest that altered expression and activation of NOX4 and FYN may contribute to the pathogenesis of pulmonary fibrosis.

Vitamin D Metabolism Is Dysregulated in Asthma and Chronic Obstructive Pulmonary Disease

Rationale: Vitamin D deficiency is common in patients with asthma and chronic obstructive pulmonary disease (COPD). Low 25-hydroxyvitamin D (25[OH]D) levels may represent a cause or a consequence of these conditions.Objectives: To determine whether vitamin D metabolism is altered in asthma or COPD.Methods: We conducted a longitudinal study in 186 adults to determine whether the 25(OH)D response to six oral doses of 3 mg vitamin D₃, administered over 1 year, differed between those with asthma or COPD versus control subjects. Serum concentrations of vitamin D₃, 25(OH)D₃, and 1α,25-dihydroxyvitamin D₃ (1α,25[OH]₂D₃) were determined presupplementation and postsupplementation in 93 adults with asthma, COPD, or neither condition, and metabolite-to-parent compound molar ratios were compared between groups to estimate hydroxylase activity. Additionally, we analyzed 14 datasets to compare expression of 1α,25(OH)₂D₃-inducible gene expression signatures in clinical samples taken from adults with asthma or COPD versus control subjects.Measurements and Main Results: The mean postsupplementation 25(OH)D increase in participants with asthma (20.9 nmol/L) and COPD (21.5 nmol/L) was lower than in control subjects (39.8 nmol/L; P = 0.001). Compared with control subjects, patients with asthma and COPD had lower molar ratios of 25(OH)D₃-to-vitamin D₃ and higher molar ratios of 1α,25(OH)₂D₃-to-25(OH)D₃ both presupplementation and postsupplementation (P ≤ 0.005). Intergroup differences in 1α,25(OH)₂D₃-inducible gene expression signatures were modest and variable if statistically significant.Conclusions: Attenuation of the 25(OH)D response to vitamin D supplementation in asthma and COPD associated with reduced molar ratios of 25(OH)D₃-to-vitamin D₃ and increased molar ratios of 1α,25(OH)₂D₃-to-25(OH)D₃ in serum, suggesting that vitamin D metabolism is dysregulated in these conditions.

Pathology of Idiopathic Pulmonary Fibrosis Assessed by a Combination of Microcomputed Tomography, Histology, and Immunohistochemistry

Idiopathic pulmonary fibrosis (IPF) is a fibrotic disease with the histology of usual interstitial pneumonia (UIP). Although the pathologist's visual inspection is central in histologic assessments, three-dimensional microcomputed tomography (microCT) assessment may complement the pathologist's scoring. We examined associations between the histopathologic features of UIP and IPF in explanted lungs and quantitative microCT measurements, including alveolar surface density, total lung volume taken up by tissue (%), and terminal bronchiolar number. Sixty frozen samples from 10 air-inflated explanted lungs with severe IPF and 36 samples from 6 donor control lungs were scanned with microCT and processed for histologic analysis. An experienced pathologist scored three major UIP criteria (patchy fibrosis, honeycomb, and fibroblastic foci), five additional pathologic changes, and immunohistochemical staining for CD68-, CD4-, CD8-, and CD79a-positive cells, graded on a 0 to 3+ scale. The alveolar surface density and terminal bronchiolar number decreased and the tissue percentage increased in lungs with IPF compared with controls. In lungs with IPF, lower alveolar surface density and higher tissue percentage were correlated with greater scores of patchy fibrosis, fibroblastic foci, honeycomb, CD79a-positive cells, and lymphoid follicles. A decreased number of terminal bronchioles was correlated with honeycomb score but not with the other scores. The three-dimensional microCT measurements reflect the pathological UIP and IPF criteria and suggest that the reduction in the terminal bronchioles may be associated with honeycomb cyst formation.

Identification of key modules and hub genes associated with lung function in idiopathic pulmonary fibrosis

Background

Idiopathic pulmonary fibrosis (IPF) is a chronic and progressive interstitial lung disease, characterized by a decline in lung function. To date, the pathophysiologic mechanisms associated with lung dysfunction remain unclear, and no effective therapy has been identified to improve lung function.

Methods

In the present study, we used weighted gene co-expression network analysis (WGCNA) to identify key modules and hub genes associated with lung function in IPF. Three datasets, containing clinical information, were downloaded from Gene Expression Omnibus. WGCNA was performed on the GSE32537 dataset. Differentially expressed gene s (DEGs) between IPF patients and healthy controls were also identified to filter hub genes. The relationship between hub genes and lung function was then validated using the GSE47460 and GSE24206 datasets.

Results

The red module, containing 267 genes, was positively correlated with the St. George’s Respiratory Questionnaire score (r = 0.37, p < 0.001) and negatively correlated with the percent predicted forced vital capacity (FVC% predicted) (r =  − 0.46, p < 0.001) and the percent predicted diffusion capacity of the lung for carbon monoxide (Dlco% predicted) (r =  − 0.42, p < 0.001). Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analysis suggested that the genes in the red module were primarily involved in inflammation and immune pathways. Based on Module Membership and Gene Significance, 32 candidate hub genes were selected in the red module to construct a protein-protein interaction network . Based on the identified DEGs and the degree of connectivity in the network, we identified three hub genes, including interleukin 6 (IL6), suppressor of cytokine signaling-3 (SOCS3), and serpin family E member 1 (SERPINE1). In the GSE47460 dataset, Spearman correlation coefficients between Dlco% predicted and expression levels of IL6, SERPINE1, SOCS3 were –0.32, –0.41, and –0.46, respectively. Spearman correlation coefficients between FVC% predicted and expression levels of IL6, SERPINE1, SOCS3 were –0.29, –0.33, and –0.27, respectively. In the GSE24206 dataset, all three hub genes were upregulated in patients with advanced IPF.

Conclusion

We identified three hub genes that negatively correlated with the lung function of IPF patients. Our results provide insights into the pathogenesis underlying the progressive disruption of lung function, and the identified hub genes may serve as biomarkers and potential therapeutictargets for the treatment of IPF patients.

Longitudinal blood transcriptomic analysis to identify molecular regulatory patterns of bovine respiratory disease in beef cattle

Bovine respiratory disease (BRD) is the most common disease in beef cattle and leads to considerable economic losses in both beef and dairy cattle. It is important to uncover the molecular mechanisms underlying BRD and to identify biomarkers for early identification of BRD cattle in order to address its impact on production and welfare. In this study, a longitudinal transcriptomic analysis was conducted using blood samples collected from 24 beef cattle at three production stages in the feedlot: 1) arrival (Entry group); 2) when identified as sick (diagnosed as BRD) and separated for treatment (Pulled); 3) prior to marketing (Close-out, representing healthy animals). Expressed genes were significantly different in the same animal among Entry, Pulled and Close-out stages (false discovery rate (FDR) < 0.01 & |Fold Change| > 2). Beef steers at both Entry and Pulled stages presented obvious difference in GO terms (FDR < 0.05) and affected biological functions (FDR < 0.05 & |Z-score| > 2) when compared with animals at Close-out. However, no significant functional difference was observed between Entry and Pulled animals. The interferon signaling pathway showed the most significant difference between animals at Entry/Pulled and Close-out stages (P < .001 & |Z-score| > 2), suggesting the animals initiated antiviral responses at an early stage of infection. Six key genes including IFI6, IFIT3, ISG15, MX1, and OAS2 were identified as biomarkers to predict and recognize sick cattle at Entry. A gene module with 169 co-expressed genes obtained from WGCNA analysis was most positively correlated (R = 0.59, P = 6E-08) with sickness, which was regulated by 11 transcription factors. Our findings provide an initial understanding of the BRD infection process in the field and suggests a subset of novel marker genes for identifying BRD in cattle at an early stage of infection.

Toll interacting protein protects bronchial epithelial cells from bleomycin-induced apoptosis

Idiopathic pulmonary fibrosis (IPF) is characterized by altered epithelial cell phenotypes, which are associated with myofibroblast accumulation in the lung. Atypical alveolar epithelial cells in IPF express molecular markers of airway epithelium. Polymorphisms within and around Toll interacting protein (TOLLIP) are associated with the susceptibility to IPF and mortality. However, the functional role of TOLLIP in IPF is unknown. Using lung tissues from IPF and control subjects, we showed that expression of TOLLIP gene in the lung parenchyma is globally lower in IPF compared to controls. Lung cells expressing significant levels of TOLLIP include macrophages, alveolar type II, and basal cells. TOLLIP protein expression is lower in the parenchyma of IPF lungs but is expressed in the atypical epithelial cells of the distal fibrotic regions. Using overexpression and silencing approaches, we demonstrate that TOLLIP protects cells from bleomycin-induced apoptosis using primary bronchial epithelial cells and BEAS-2B cells. The protective effects are mediated by reducing mitochondrial reactive oxygen species (ROS) levels and upregulating autophagy. Therefore, global downregulation of the TOLLIP gene in IPF lungs may predispose injured lung epithelial cells to apoptosis and to the development of IPF.

Identifying the Transcriptional Regulatory Network Associated With Extrathyroidal Extension in Papillary Thyroid Carcinoma by Comprehensive Bioinformatics Analysis

Extrathyroidal extension (ETE) affects papillary thyroid cancer (PTC) prognosis. The objective of this study was to identify biomarkers for ETE and explore the mechanisms controlling its development in PTC. We performed a comprehensive bioinformatics analysis using several datasets. Differential expression analysis and weighted gene co-expression network analysis (WGCNA) on 58 paired PTC samples from The Cancer Genome Atlas (TCGA) were used to detect ETE-related mRNA and long noncoding (lnc) RNA modules and construct an lncRNA/mRNA network. An independent TCGA dataset containing 438 samples was utilized to validate and characterize the WGCNA results. Functional annotation was used to identify the biological functions and related pathways of ETE modules. Two independent RNA sequencing datasets were combined to crossvalidate relationships between lncRNAs and mRNAs by Pearson correlation analysis. Transcription factors (TFs) for affected genes were predicted using the binding motif data from Ensembl Biomart to construct a TF/lncRNA/mRNA network. Other two independent datasets were used to crossvalidate TF-mRNA associations. Finally, receiver operating characteristic, survival analyses, and Cox proportional hazard regression model were performed to explore the significance of hub genes in ETE diagnosis and PTC prognosis. Three mRNA modules and two lncRNA modules were significantly associated with ETE. Enrichment analysis showed extracellular matrix changes was closely related to the development of ETE. A TF/lncRNA/mRNA regulatory network was constructed containing 33 validated hub genes, 64 lncRNAs, and 64 TFs, all differentially expressed between ETE and non-ETE samples. Unc-5 family C-terminal like [area under the curve (AUC): 0.711], sushi repeat containing protein X-linked 2 (AUC: 0.706), lysyl oxidase (AUC: 0.704), collagen type I alpha 1 chain (AUC: 0.704), and collagen type X alpha 1 chain (AUC: 0.704) were the most highly significant hub genes for ETE diagnosis. The Cox proportional hazard regression model constructed with hub genes showed significant survival differences between low- and high-risk groups (p = 0.00025) and performed good prediction for PTC prognosis(AUC = 0.794; C-index = 0.895). The identification of 33 biomarkers and TF/lncRNA/mRNA regulatory network would provide new insights into the molecular mechanisms of ETE besides the prognosis model may have important clinical implications in the improvement of PTC risk stratification, therapeutic decision-making, and prognosis prediction.

Recapitulating idiopathic pulmonary fibrosis related alveolar epithelial dysfunction in a human iPSC-derived air-liquid interface model

Idiopathic pulmonary fibrosis (IPF) is a fatal disease of unknown cause that is characterized by progressive fibrotic lung remodeling. An abnormal emergence of airway epithelial-like cells within the alveolar compartments of the lung, herein termed bronchiolization, is often observed in IPF. However, the origin of this dysfunctional distal lung epithelium remains unknown due to a lack of suitable human model systems. In this study, we established a human induced pluripotent stem cell (iPSC)-derived air-liquid interface (ALI) model of alveolar epithelial type II (ATII)-like cell differentiation that allows us to investigate alveolar epithelial progenitor cell differentiation in vitro. We treated this system with an IPF-relevant cocktail (IPF-RC) to mimic the pro-fibrotic cytokine milieu present in IPF lungs. Stimulation with IPF-RC during differentiation increases secretion of IPF biomarkers and RNA sequencing (RNA-seq) of these cultures reveals significant overlap with human IPF patient data. IPF-RC treatment further impairs ATII differentiation by driving a shift toward an airway epithelial-like expression signature, providing evidence that a pro-fibrotic cytokine environment can influence the proximo-distal differentiation pattern of human lung epithelial cells. In conclusion, we show for the first time, the establishment of a human model system that recapitulates aspects of IPF-associated bronchiolization of the lung epithelium in vitro.

Regulation of the Extracellular Matrix by Ciliary Machinery

The primary cilium is an organelle involved in cellular signalling. Mutations affecting proteins involved in cilia assembly or function result in diseases known as ciliopathies, which cause a wide variety of phenotypes across multiple tissues. These mutations disrupt various cellular processes, including regulation of the extracellular matrix. The matrix is important for maintaining tissue homeostasis through influencing cell behaviour and providing structural support; therefore, the matrix changes observed in ciliopathies have been implicated in the pathogenesis of these diseases. Whilst many studies have associated the cilium with processes that regulate the matrix, exactly how these matrix changes arise is not well characterised. This review aims to bring together the direct and indirect evidence for ciliary regulation of matrix, in order to summarise the possible mechanisms by which the ciliary machinery could regulate the composition, secretion, remodelling and organisation of the matrix.

Hemagglutinin Cleavability, Acid Stability, and Temperature Dependence Optimize Influenza B Virus for Replication in Human Airways

Influenza A virus (IAV) and influenza B virus (IBV) cause yearly epidemics with significant morbidity and mortality. When zoonotic IAVs enter the human population, the viral hemagglutinin (HA) requires adaptation to achieve sustained virus transmission. In contrast, IBV has been circulating in humans, its only host, for a long period of time. Whether this entailed adaptation of IBV HA to the human airways is unknown. To address this question, we compared two seasonal IAVs (A/H1N1 and A/H3N2) and two IBVs (B/Victoria and B/Yamagata lineages) with regard to host-dependent activity of HA as the mediator of membrane fusion during viral entry. We first investigated proteolytic activation of HA by covering all type II transmembrane serine protease (TTSP) and kallikrein enzymes, many of which proved to be present in human respiratory epithelium. The IBV HA0 precursor is cleaved by a broader panel of TTSPs and activated with much higher efficiency than IAV HA0. Accordingly, knockdown of a single protease, TMPRSS2, abrogated spread of IAV but not IBV in human respiratory epithelial cells. Second, the HA fusion pH values proved similar for IBV and human-adapted IAVs (with one exception being the HA of 1918 IAV). Third, IBV HA exhibited higher expression at 33°C, a temperature required for membrane fusion by B/Victoria HA. This indicates pronounced adaptation of IBV HA to the mildly acidic pH and cooler temperature of human upper airways. These distinct and intrinsic features of IBV HA are compatible with extensive host adaptation during prolonged circulation of this respiratory virus in the human population.IMPORTANCE Influenza epidemics are caused by influenza A and influenza B viruses (IAV and IBV, respectively). IBV causes substantial disease; however, it is far less studied than IAV. While IAV originates from animal reservoirs, IBV circulates in humans only. Virus spread requires that the viral hemagglutinin (HA) is active and sufficiently stable in human airways. We resolve here how these mechanisms differ between IBV and IAV. Whereas human IAVs rely on one particular protease for HA activation, this is not the case for IBV. Superior activation of IBV by several proteases should enhance shedding of infectious particles. IBV HA exhibits acid stability and a preference for 33°C, indicating pronounced adaptation to the human upper airways, where the pH is mildly acidic and a cooler temperature exists. These adaptive features are rationalized by the long existence of IBV in humans and may have broader relevance for understanding the biology and evolution of respiratory viruses.

Transcriptional regulatory model of fibrosis progression in the human lung

To develop a systems biology model of fibrosis progression within the human lung we performed RNA sequencing and microRNA analysis on 95 samples obtained from 10 idiopathic pulmonary fibrosis (IPF) and 6 control lungs. Extent of fibrosis in each sample was assessed by microCT-measured alveolar surface density (ASD) and confirmed by histology. Regulatory gene expression networks were identified using linear mixed-effect models and dynamic regulatory events miner (DREM). Differential gene expression analysis identified a core set of genes increased or decreased before fibrosis was histologically evident that continued to change with advanced fibrosis. DREM generated a systems biology model (www.sb.cs.cmu.edu/IPFReg) that identified progressively divergent gene expression tracks with microRNAs and transcription factors that specifically regulate mild or advanced fibrosis. We confirmed model predictions by demonstrating that expression of POU2AF1, previously unassociated with lung fibrosis but proposed by the model as regulator, is increased in B lymphocytes in IPF lungs and that POU2AF1-knockout mice were protected from bleomycin-induced lung fibrosis. Our results reveal distinct regulation of gene expression changes in IPF tissue that remained structurally normal compared with moderate or advanced fibrosis and suggest distinct regulatory mechanisms for each stage.

Cognitive impairment and transcriptomic profile in hippocampus of young mice after multiple neonatal exposures to sevoflurane

Children with repeated inhalational anesthesia may develop cognitive disorders. This study aimed to investigate the transcriptome-wide response of hippocampus in young mice that had been exposed to multiple sevoflurane in the neonatal period. Mice received 3% sevoflurane for 2 h on postnatal day (PND) 6, 8, and 10, followed by arterial blood gas test on PND 10, behavioral experiments on PND 31–36, and RNA sequencing (RNA-seq) of hippocampus on PND 37. Functional annotation and protein-protein interaction analyses of differentially expressed genes (DEGs) and quantitative reverse transcription polymerase chain reaction (qPCR) were performed. Neonatal sevoflurane exposures induced cognitive and social behavior disorders in young mice. RNA-seq identified a total of 314 DEGs. Several enriched biological processes (ion channels, brain development, learning, and memory) and signaling pathways (oxytocin signaling pathway and glutamatergic, cholinergic, and GABAergic synapses) were highlighted. As hub-proteins, Pten was involved in nervous system development, synapse assembly, learning, memory, and behaviors, Nos3 and Pik3cd in oxytocin signaling pathway, and Cdk16 in exocytosis and phosphorylation. Some top DEGs were validated by qPCR. This study revealed a transcriptome-wide profile in mice hippocampus after multiple neonatal exposures to sevoflurane, promoting better understanding of underlying mechanisms and investigation of preventive strategies.