Junling Pang1†
Ya Luo1†
Dong Wei2†
Zhujie Cao1
Xianmei Qi1Meiyue Song3
Ying Liu1Zhaoguo Li4
Jin Zhang5
Baicun Li6*
Jingyu Chen2,5*
Jing Wang1*Chen Wang6- 1State Key Laboratory of Medical Molecular Biology, Department of Pathophysiology, Peking Union Medical College, Institute of Basic Medical Sciences, Chinese Academy of Medical Sciences, Beijing, China
- 2Transplant Center, Wuxi People’s Hospital Affiliated to Nanjing Medical University, Wuxi, China
- 3Department of Pulmonary and Critical Care Medicine/Others, Center of Respiratory Medicine, China-Japan Friendship Hospital, Beijing, China
- 4Department of Respiratory, The Second Affiliated Hospital of Harbin Medical University, Harbin, China
- 5Department of Thoracic Surgery and Lung Transplantation, China-Japan Friendship Hospital, Beijing, China
- 6State Key Laboratory of Medical Molecular Biology, Department of Physiology, Peking Union Medical College, Institute of Basic Medical Sciences, Chinese Academy of Medical Sciences, Beijing, China
Silicosis is a fatal occupational lung disease which currently has no effective clinical cure. Recent studies examining the underlying mechanism of silicosis have primarily examined experimental models, which may not perfectly reflect the nature of human silicosis progression. A comprehensive profiling of the molecular changes in human silicosis lungs is urgently needed. Here, we conducted RNA sequencing (RNA-seq) on the lung tissues of 10 silicosis patients and 7 non-diseased donors. A total of 2,605 differentially expressed genes (DEGs) and critical pathway changes were identified in human silicosis lungs. Further, the DEGs in silicosis were compared with those in idiopathic pulmonary fibrosis (IPF) and chronic obstructive pulmonary diseases (COPD), to extend current knowledge about the disease mechanisms and develop potential drugs. This analysis revealed both common and specific regulations in silicosis, along with several critical genes (e.g., MUC5AC and FGF10), which are potential drug targets for silicosis treatment. Drugs including Plerixafor and Retinoic acid were predicted as potential candidates in treating silicosis. Overall, this study provides the first transcriptomic fingerprint of human silicosis lungs. The comparative transcriptome analyses comprehensively characterize pathological regulations resulting from silicosis, and provide valuable cues for silicosis treatment.
Introduction
Silicosis is a progressive fibrotic lung disease, which arises from the accumulation of inhaled silica particles during occupational activities (Hoy and Chambers, 2020). Increasing workers are threatening by crystalline silica due to emerging various modern industries (Barnes et al., 2019). Many countries have recently experienced outbreaks of silicosis, even in some developed countries (Rose et al., 2019). Based on the Global Burden of Disease Study, the prevalence of silicosis is about 162.4 thousand from 1990 to 2017, with 23.7 thousand new cases in the year 2017 (GBD 2017 Disease and Injury Incidence and Prevalence Collaborators, 2018). Unfortunately, current treatment for silicosis is still limited, and no effective treatments could be used to reverse its progression. For advanced silicosis, lung transplantation is one of the therapeutic options to prolong survival time (Singer et al., 2012). However, considering its several limitations, such as a short median survival rate, many contraindications, and low availability of lung donors, lung transplant is not a routine recommendation for silicosis patients (Association, 2018). Although many studies have tried to identify the pathogenesis and potential therapeutic targets of silicosis through rodent models, these results have not been translated successfully into clinical applications (Lopes-Pacheco et al., 2016; Pollard, 2016).
RNA sequencing (RNA-seq) has long been used to extend our understanding of molecular functions across a wide range of diseases (Stark et al., 2019). Several previous studies have employed RNA-seq to reveal the comprehensive molecular alterations in silicosis, using either lung tissues from animal models (Chen J. et al., 2018; Sai et al., 2019) or human cell lines (Chan et al., 2017, 2018). However, because of the genetic differences between animals and humans and the fact that the development of silicosis involves diverse cell types, a comprehensive profile of the molecular changes in the lungs of human silicosis patients is still urgently needed.
Silicosis and idiopathic pulmonary fibrosis (IPF) are both progressively fibrosing interstitial lung diseases (Wollin et al., 2019). Additionally, as a chronic inflammatory lung disease, silicosis shares some common pathogenic mechanisms with chronic obstructive pulmonary diseases (COPD) (Beamer and Shepherd, 2013). Further, epidemiological studies have shown that patients with silicosis are prone to accompany with IPF and COPD (Fazen et al., 2020; Walters, 2020). Identifying the shared and specific pathways that are modulated by the three diseases will contribute to a better understanding of the mechanisms underpinning silicosis progression. Moreover, because there are no effective drugs for treating silicosis, characterizing the mechanisms that silicosis shares with IPF or COPD may point the way to a novel use of available drugs for treating silicosis.
In this project, we performed RNA-seq analyses using lung tissues from 10 silicosis patients and 7 non-diseased donors. Through bioinformatic analysis, we aimed to reveal the key genes and pathways that participate in the progression of silicosis. In addition, we conducted transcriptomic comparisons to identify the common and unique regulations associated with silicosis, IPF (Sivakumar et al., 2019) and COPD (Morrow et al., 2017) to expand current understanding of the molecular changes in silicosis.
Materials and Methods
A detailed description can be found in the online Supplementary Methods. The scripts used in the analysis can be found in the github repository1.
Human Samples
All human lung samples were obtained from Wuxi People’s Hospital Affiliated to Nanjing Medical University. Lung tissues of silicosis patients (n = 10) were obtained when they received lung transplantation. Each case (transplanted lung) was sampled close to the edge of the lung, avoiding lymph nodes, large vessels, and atmospheric tubes. Lung tissues from non-diseased donors (n = 7) were the resected parts during lung transplantation due to size incompatibility. Though the position was not completely fixed, lymph nodes, large vessels, and atmospheric tubes were avoided during sampling. The demographic and clinical characteristics of the subjects were recorded in Supplementary Tables 1, 2. This study was approved by the Institutional Review Boards of Nanjing Medical University, and Institute of Basic Medical Sciences, Chinese Academy of Medical Sciences. All participants gave written informed consent.
RNA Sequencing
The RNA extraction, mRNA library construction, and RNA-seq were all performed by BGI-Shenzhen, China. Briefly, the total RNA was extracted from each human lung tissues using Trizol (Invitrogen, Carlsbad, CA, United States) following the manual instruction. To purify mRNA, oligo(dT)-attached magnetic beads were used. After library construction, 100-bp paired-end reads were sequenced using BGIseq500 platform (BGI-Shenzhen, China). The raw sequencing data was filtered by removing sequencing adapters, and low-quality reads. The cleaned reads stored in FASTQ format were used in the following bioinformatics analysis.
Bioinformatics Analysis for the RNA-Seq Data
The human reference genome (GRCh38) and the gene annotation file in ‘‘gtf’’ format were both downloaded from Ensembl release 84 version2. The cleaned RNA-seq reads of each sample were aligned to the human reference genome using HISAT2 (v2.1.0) software with default parameters (Kim et al., 2015; Pertea et al., 2016). Only the reads that were uniquely mapped to the genome were retained for further analysis (Supplementary Table 3). R packages GenomicFeatures (v1.38.2) and GenomicAlignments (v1.22.1) were used for calculating the read counts of genes (Lawrence et al., 2013). DEGs between silicosis patients and non-diseased donors were got using DESeq2 (v1.26.0) package (Love et al., 2014). Since there was no significant difference between silicosis patients and non-diseased donors regarding to age, sex, body mass index, and smoking status (Supplementary Table 1), we did not consider these factors in the comparison between the silicosis group and the non-diseased group. Known batch effects (see “SampleInfo.xls” in above github repository) were considered when we performed DESeq2 to identify DEGs. The cutoff to select DEGs in the diseased group was set: adjusted p < 0.1, and absolute fold change > 1.5 compared to non-diseased group. The R package ggplot2 (v3.3.3) was used for the generation of volcano plot for the DEGs (Wickham, 2016). Another R package pheatmap (v1.0.12)3 was used for the generation of heatmap for the DEGs.
Differentially Expressed Genes in IPF or COPD
The DEGs (adjusted p < 0.1; absolute fold change > 1.5) in the lungs of IPF patients (n = 36) compared to non-diseased controls (n = 19) were accessed from the supplementary materials of published paper (Sivakumar et al., 2019).
The raw data of the expression profiling for lung tissues of severe COPD and controls were downloaded from NCBI GEO datasets (GSE76925) (Morrow et al., 2017). The Microarray data was background corrected, log2 transformed and quantile normalized. An empirical Bayes shrinkage method in limma package (v3.42.2) (Ritchie et al., 2015) was used to obtain the log fold change (logFC) of genes, and the p-values. The cutoff to select the significantly altered genes was also set to: adjusted p < 0.1; absolute fold change > 1.5.
Functional Enrichment Analysis
The R package clusterProfiler (v3.14.3) was used for the Gene Ontology (GO) enrichment analysis of genes [19]. Pathway analysis was performed using MetaCore online4.
Correspondence-At-Top (CAT) Curves
We used the R package “matchBox” (v1.28.0) (Irizarry et al., 2005; Ross et al., 2011) to generate the CAT curves. To compare the DEGs among diseases or datasets, we ranked genes by decreasing log fold change. Then we used “equalRank” method to compute CAT curves. The proportion of expected top ranked features was set to 0.3 when computing probability intervals. To compare the pathways between datasets, we ordered the pathways with increasing FDRs, meaning the most significant ones are at the top of the list.
Prediction of Potential Pharmaceuticals for Silicosis
(1) To performed connectivity mapping using CMAP5, the symbols of DEGs were first transferred to probe IDs based on the information of Affymetrix Human Genome U133A 2.0 Array. Then the top 500 up- and down- probe sets were uploaded as signatures in CMAP. (2) The list of DEGs in silicosis were used as input in the ‘‘Drug Lookup’’ function of MetaCore (see footnote) to search the potential drugs. (3) We also searched some key genes in DrugBank database6 to check whether any drugs may target them.
Quantitative Real-Time PCR and Immunohistochemical Staining
The detailed description of quantitative real-time PCR and immunohistochemical staining can be found in the online Supplementary Methods.
Results
Differential Gene Analysis Revealed Wide Transcriptomic Alterations in the Lungs of Silicosis Patients Compared to Non-diseased Donors
To reveal the transcriptional changes in silicosis patients, we analyzed the gene expression in the lungs of silicosis patients (n = 10) versus non-diseased donors (n = 7). A total of 2,605 significantly altered genes between the two groups were identified, in which 1,169 genes were up-regulated and 1,436 genes were down-regulated in silicosis patients (Figure 1A and Supplementary Tables 4, 5). A heatmap of the top 20 up- and down-regulated genes (based on the adjusted p-value) depicts the hierarchical clusters between the two groups (Figure 1B).
Figure 1. Transcriptional differences between the lungs of silicosis patients and non-diseased donors. (A) Volcano plot showing the significantly altered genes in silicosis lungs. Blue dots represent down-regulated genes, red dots indicate up-regulated genes, and gray dots indicate genes with no significant change. (B) Heatmap illustrating the top 20 most up- or down-regulated genes in the lungs of silicosis patients. Each row represents a gene, and each column a sample. The expression level of each gene was normalized using row z-score. Red indicates higher expression levels and blue indicates lower expression levels. NS, no significance.
Among the top 20 up-regulated genes, COL3A1, COL5A2, and COL10A1 encode collagen proteins, which are important components of the extracellular matrix (ECM) (Saito et al., 2018). In addition, we identified 2 immune response related genes, ROR2 and SLAMF6, which may play important roles in the pulmonary inflammation that accompanies silicosis (Cai et al., 2016; Wang et al., 2016). We also identified 4 genes (SULF1, FHL2, COMP, and SPP1) with profibrotic effect (Perkins et al., 2018; Posey et al., 2018; Jin et al., 2019; Morse et al., 2019). Amongst the top 20 down-regulated genes, we identified 11 genes (CFAP73, CFAP46, CFAP157, CFAP74, CFAP65, CFAP100, CCDC114, TEKT4, ARMC3, DNAAF1, and DNAH12) which are involved in encoding the cilia and flagella associated proteins that are essential for maintaining normal cilia function (Figure 1B; Bustamante-Marin and Ostrowski, 2017).
Functional Enrichment Analyses Uncovered Predominant Pathway Changes in the Lungs of Silicosis Patients
We performed GO enrichment analysis for the up-regulated and down-regulated genes separately. As shown in Figure 2A, among the top ten significantly enriched biological processes for the up-regulated genes, three (marked by ∗) are related to promoting ECM remodeling and five (marked by #) are associated with regulation of lymphocyte activation, lymphocyte differentiation, and leukocyte migration (Supplementary Table 6). In contrast, the down-regulated genes were mainly associated with processes related to cilium organization, assembly, and movement (Figure 2B and Supplementary Table 7).
Figure 2. Functional enrichments for the significantly altered genes in silicosis lungs. (A) The top 10 biological processes for up-regulated genes. Processes related to ECM remodeling are indicated with *, while those related to immune related processes are indicated with #. The size of dots reflects the gene counts for each process. Dot color reflects adjusted p-values, where the most significant result is red. (B) The top 10 biological processes for down-regulated genes. (C) The top 10 significantly enriched pathways for up-regulated genes. The x-axis indicates the value of –log10FDR for each pathway. The ratios of gene count to the total genes for each pathway are shown in parentheses. (D) The top 10 significantly enriched pathways for the down-regulated genes.
Similar results were found using pathway-based enrichment analysis. Pathways involving “ECM remodeling,” “TGFβ-induced fibroblast/myofibroblast migration and extracellular matrix production in asthmatic airways,” as well as immune-related pathways including “Role of CD8 + Tc1 cells in COPD,” “Classical complement pathway,” “Lectin induced complement pathway,” and “T cell co-signaling receptors, schema” were highly represented for the genes up-regulated in silicosis (Figure 2C and Supplementary Table 8). In contrast, the down-regulated genes mainly enriched in pathways including “Regulation of lung epithelial progenitor cell differentiation,” and “Th2-cytokines induced expression of Mucin 5AC and maintenance of mucous cell phenotype in normal airway epithelium” (Figure 2D and Supplementary Table 9).
Taken together, these findings indicate that mucociliary dysfunction, decreased lung epithelial progenitor cell differentiation, and enhanced profibrotic and immune-related pathways may play critical roles in silicosis progression. Critical genes in these pathways should be the focus of future investigations (Supplementary Table 10).
Comparative Analyses Identified Common and Specific Regulations in Silicosis Lungs Compared to IPF and COPD
To identify the shared and unique pathways among silicosis, IPF and COPD, we compared the DEGs in the lungs of silicosis with those previously reported for IPF (Sivakumar et al., 2019) and COPD (Morrow et al., 2017). As shown in Figure 3A, the significantly altered genes in the three lung diseases differed substantially from each other. Only 129 genes were significantly altered among all three lung diseases. In contrast, 1,490 genes (57.2%) were exclusively regulated in silicosis (Supplementary Table 11). The 129 genes that were differentially regulated in all three diseases were involved in pathways relevant for antigen presentation by MHC class II molecules, B cell signaling, T cell co-signaling receptors, SHH signaling, Hedgehog signaling and Wnt signaling (Supplementary Table 12). In contrast, the 1,490 genes specifically regulated in silicosis were mainly involved in antigen presentation by MHC class I molecules, CD8 + Tcl cells, mucociliary clearance function, epithelial progenitor cells differentiation, and OX40L/OX40 signaling (Figure 3B and Supplementary Table 13). A set of genes with canonical functions in these pathways are shown in Supplementary Table 14.
Figure 3. Comparison of the significantly altered genes and their pathways in silicosis, IPF, and COPD. (A) Venn diagram of the significantly altered genes in silicosis, IPF, and COPD. (B) The top 5 pathways enriched for the 1,490 genes exclusively regulated in silicosis. The enrichment results of the corresponding pathways in IPF and COPD are given for comparison. For each pathway the gene count of the differentially expressed genes compared to the total gene counts are given in parentheses. (C) CAT plot of the common DEGs shared among silicosis, IPF, and COPD. (D) CAT plot of the common significant pathways enriched for the DEGs in silicosis and IPF. (E) CAT plot of the common significant pathways enriched for the DEGs in silicosis and COPD.
We then performed a global comparison among the three diseases using Correspondence-At-Top (CAT) curves. It was shown that the DEGs in silicosis lungs more closely resemble those got from IPF lungs than those derived from COPD (Figure 3C). The CAT curves using the ordered pathway lists for the three diseases also revealed that silicosis and IPF shares more common regulatory pathways (Figure 3D) than silicosis versus COPD (Figure 3E). In the following, we will compare silicosis with IPF and COPD separately to detect potential common or different mechanisms among diseases.
Genes Regulated in Both Silicosis and IPF Suggest Common Profibrotic and Proinflammatory Mechanisms in Disease Development
There were 1,001 significantly altered genes shared between silicosis and IPF (Figure 3A). Among them, 66.9% (670 genes) were regulated in the same direction (500 up-regulated and 170 down-regulated in both diseases) (Figure 4A).
Figure 4. Functional analysis of the genes regulated in both silicosis and IPF. (A) Statistics of the gene counts regulated in both silicosis and IPF. Pink squares indicate genes regulated in the same direction (up or down) in both diseases; blue squares indicate genes regulated in different directions in the two diseases. (B) The top 10 pathways significantly enriched for the 500 genes up-regulated in both silicosis and IPF. (C) The top 10 significant biological processes enriched for the 500 genes up-regulated in both silicosis and IPF. The size of dots reflects the gene counts for each process. Dot color reflects adjusted p-values, where the most significant result is red. (D) Representative genes selected from significant pathways shared by silicosis and IPF.
To reveal the common pathological mechanisms between silicosis and IPF, we conducted functional enrichment analysis for these genes. It was shown that the 500 up-regulated genes were largely involved in the pathways and biological processes participating in pulmonary fibrosis and immune response (Figures 4B,C and Supplementary Tables 15, 16). Relevant pathways related to pulmonary fibrosis included “ECM remodeling” and “TGFβ-induced fibroblast/myofibroblast migration and extracellular matrix production in asthmatic airways.” Immune response related pathways included “Classical complement pathway,” “Lectin induced complement pathway,” “Common mechanisms of Th17 cell migration,” and “T cell co-signaling receptors” (Figure 4B). Representative genes participating in these critical pathways are presented in Figure 4D.
Unexpectedly, the 170 down-regulated genes had no significant GO enrichment results, and only two pathways were revealed through pathway-based enrichment analysis (Supplementary Table 17). As a result, we have not discussed the function of these genes further.
Genes Regulated in Both Silicosis and COPD Indicate Common Proinflammatory Mechanisms in Disease Development
There were total 243 genes that were significantly altered both in silicosis and COPD lungs (Figure 3A). About 70.0%, including 124 up-regulated and 46 down-regulated genes, were regulated in the same direction in both diseases (Figure 5A).
Figure 5. Functional analysis of the genes regulated in both silicosis and COPD. (A) Statistics of the gene counts regulated in both silicosis and COPD. Pink squares indicate genes regulated in the same direction (up or down) in both diseases; blue squares indicate genes regulated in different directions in the two diseases. (B) The top 10 significant pathways enriched for the 124 genes up-regulated in both silicosis and COPD. (C) The top 10 biological processes significantly enriched for the 124 genes up-regulated in both silicosis and COPD. The size of dots reflects the gene counts for each process. Dot color reflects adjusted p-values, where the most significant result is red. (D) Representative genes selected from significant pathways shared by silicosis and COPD.
Pathway-based analysis showed that the 124 up-regulated genes were mainly involved in immune response pathways such as “Antigen presentation by MHC class II,” “infiltrating B cells in anti-tumor immunity,” and “T follicular helper cell dysfunction in SLE” (Figure 5B and Supplementary Table 18). Consistently, GO analysis revealed the enrichment on “Lymphocyte differentiation,” “T cell activation,” and “Extracellular structure organization” (Figures 5C,D and Supplementary Table 19).
In contrast, the 46 down-regulated genes were enriched on the pathways including “Hedgehog signaling,” “Negative regulation of Wnt/β-catenin signaling in the cytoplasm,” and “Neolacto-series GSL Metabolism” (Supplementary Table 20). In addition, the two biological processes “Cilium assembly” and “Cilium organization” were also detected through GO enrichment analysis (Figure 5D and Supplementary Table 21).
A Summary of Key Pathways and Genes Promoting Silicosis, Experimental Validation, and Potential Drug Prediction
Taken together with results from previous studies, our findings provide a comprehensive view of regulatory mechanisms in silicosis (Figure 6). Except for the common immune response pathways shared with IPF and COPD (antigen presentation by MHC class II molecules, B cell signaling, and T cell co-signaling receptors), silica exposure in silicosis may specifically drive antigen presentation by MHC class I molecules, which are recognized by CD8 + T cells, leading to damage to pulmonary epithelial cells through cytotoxic pathways, followed by pulmonary fibrosis (Zhang and Zhang, 2020). The mucociliary dysfunction combined with the weakened epithelial progenitor cells differentiation may worsen the damage of respiratory epithelium (Kotton and Morrisey, 2014; Benam et al., 2018), resulting in the up-regulation of several profibrotic pathways including TGF-β signaling and Wnt/β-catenin signaling, followed by epithelial-to-mesenchymal transition (EMT) and pulmonary fibrosis (Kage and Borok, 2012; Jones et al., 2019). Further, the activation of immune cells (e.g., Th2 cells) and their related signaling pathways such as OX40L/OX40 signaling may also cause downstream profibrotic interactions (So et al., 2006). In addition, collagen and matrix deposition promoting ECM remodeling and lung fibrosis has been widely accepted as relevant to silicosis. All these pathways (1–7 in Figure 6) eventually contribute to the inflammation and fibrosis in silicosis. In contrast, the inhibited Hedgehog signaling in silicosis (8 in Figure 6) suggests a protective feedback response to chronic epithelial injury (Zhang et al., 2019) via a mesenchymal feedback mechanism, aimed at promoting regeneration (Peng et al., 2015).
Figure 6. Illustration of the proposed underlying mechanisms involved in silicosis progression based on evidence from the comparative transcriptome analysis. Silica exposure promotes silicosis progression via indicated pathways (1–8). Key pathways (1–3, and 6) and genes highlighted in red are proposed to be differentially important in silicosis compared to IPF and COPD. →: promote silicosis progression; ⊤: inhibit silicosis progression; ↑: up-regulated (pathways or genes); ↓: down-regulated (pathways or genes); ←: a cross-talk between pathways.
To validate the results got from the RNA-seq data, we performed real-time quantitative PCR and immunohistochemistry (Supplementary Tables 22, 23) for some key genes that summarized in Figure 6. As shown in Figure 7A, except one gene (FGFR2), all other genes we tested had significantly altered mRNA levels in the lungs of patients with silicosis, the same trends as revealed by RNA-seq. Further, the protein levels of MUC5AC and FGF10 detected by immunohistochemistry were significantly decreased in silicosis lungs, especially in epithelium (Figures 7B,C). These results indicated that the genes revealed by RNA-seq could be promising therapeutic targets of silicosis.
Figure 7. Validation of key genes revealing by RNA-seq through experimental approaches. (A) Relative mRNA expression levels of the genes (HLA-B, MUC5AC, FOXJ1, FGF10, FGFR2, CXCR4, CTHRC1, FHL2, MMP2, and TIMP1) in the lung tissues from silicosis patients (n = 10) or non-diseased donors (n = 7). (B) Representative images of immunohistochemical staining of MUC5AC and FGF10 in the lung tissues from silicosis patients or non-diseased donors. Red color indicates the positive staining. All scale bars are 100 μm. (C) Quantified immunoreactivity of MUC5AC and FGF10 in lung sections from silicosis patients (n = 10) or non-diseased donors (n = 7). All the quantitative results are presented as mean ± SEM. The differences were analyzed by two-tailed Student’s t-test. ns, no significance, *P < 0.05, **P < 0.01, ***P < 0.001.
Finally, to search the potential pharmaceuticals for silicosis, we performed connectivity mapping using CMAP to detect potential perturbagens that may reverse the biological states in silicosis (Supplementary Table 24). We also searched the potential drugs that target the key genes revealed in the RNA-seq data using MetaCore and DrugBank (Supplementary Table 25). Among the results, several drugs are of our interest. For example, Plerixafor inhibits the CXCR4 chemokine receptor, which is induced in silicosis; Retinoic acid may both inhibit the expression of COL2A1 and activate MUC5AC. Further drug tests on animal models are worthwhile to check their effects on silicosis.
Discussion
Our study provides the first systematic exploration of the transcriptional landscape in human silicosis lungs, with the aim of identifying relevant underlying mechanisms and potential therapeutic targets for silicosis treatment. Compared to non-diseased donors, pathways involved in pro-fibrosis and immune responses were found to be significantly enhanced in silicosis lungs, and pathways related to epithelial injury and mucociliary dysfunction were also identified. The comparison between silicosis, IPF, and COPD revealed regulated pathways that were both common across diseases and specific to silicosis.
Silicosis is mainly characterized by chronic inflammation and progressive fibrosis. A growing number of studies have revealed that immune dysfunction exerts an essential role in initiating inflammation in early-stage silicosis and promoting fibrosis in the later stage. Previous studies have shown that an increase in antigen-presenting cell (APC) activity contributes to the inflammatory process. The major disease-related genes of silicosis may be mapped near the HLA-B locus (Honda et al., 1993). A recent study found an increased expression of the activation marker HLA-DR on CD8 + T cells of silica exposed workers (Brilland et al., 2019). Similarly, our transcriptomics data showed that classic APC markers of MHC class I, especially HLA-B, were specifically increased, promoting lung inflammation. The activation and infiltration of lymphocytes are reported to have important roles in silicosis. Consistently, we also observed a large number of lymphocytes infiltrated into the lungs of the advanced silicosis patients (Supplementary Table 2). Previous research demonstrated that T-helper (Th) 1 type inflammation exists in the early stage of silicosis, which is characterized by increased levels of tumor necrosis factor alpha (TNF-α) and interferon gamma (IFN-γ) (Migliaccio et al., 2008). Subsequently, the activation of Th2 cells and related signaling pathways, such as OX40L/OX40 signaling, can promote fibrotic changes (So et al., 2006). Furthermore, the regulation of C-X-C motif chemokine receptor 4/C-X-C motif chemokine ligand 12 (CXCR4/CXCL12) axis could maintain the Th2 bias, and enhance the expression of cytokines with profibrotic effect such as IL-4 and IL-10 (Piao et al., 2012; Li et al., 2020). Additionally, research on B cells in silicosis has mainly focused on IL-10 producing regulatory B cells (B10 cells). B10 cells can reduce Th1 cells and inflammation via regulating the balance of T cells, thus alleviating silica-induced lung inflammation (Liu et al., 2016; Chen et al., 2017; Lu et al., 2017). Consistently, the transcriptomic analysis of patients with advanced silicosis revealed up-regulated expression of IFN-γ, IL-10 and CXCR4, and the OX40L/OX40 signaling pathway.
Mucins and cilia play crucial roles in protecting respiratory epithelium from injury by removing inhaled particles and pathogens from the airway (Benam et al., 2018). Previous studies found that the assembly of cilia is remarkably decreased but the secretion of mucins (encoded by MUC5AC and MUC5B) is increased in COPD, asthma, IPF, and cystic fibrosis (Benam et al., 2018; Samsuzzaman et al., 2019). In contrast, our data showed that the genes and pathways associated with both the production of mucins and the assembly of cilia are significantly down-regulated in silicosis, indicating a novel feature of mucociliary dysfunction induced by silica exposure. In support of our findings, previous research has reported that MUC5AC deficient mice have severely impaired mucociliary clearance and acute lung injury (Koeppen et al., 2013; Roy et al., 2014). Increased secretion of MUC5AC in airway epithelial cells exhibits a protective role against influenza-induced lung injury but does not cause airway obstruction (Ehre et al., 2012). Combined with our current data, improving mucociliary function via increasing the expression of MUC5AC may be an effective way to ameliorate silicosis progression.
Lung epithelial progenitor cells play a key role in effective regeneration and repair of injured lung epithelium. Failure of lung epithelium regeneration could cause pulmonary fibrosis (Kotton and Morrisey, 2014). In our results, the pathway of lung epithelial progenitor cell differentiation, including several essential genes (e.g., FGF10, FGFR2), was remarkably down-regulated in silicosis, indicating disordered respiratory epithelium regeneration in the disease progression. Fibroblast growth factor 10 (FGF10) is essential to regulate lung epithelial regeneration after injury via binding to and activating fibroblast growth factor receptor 2 (FGFR2) on epithelial progenitors (Abler et al., 2009). FGF10-FGFR2B signaling also helps to maintain the homeostasis of alveolar type 2 (AT2) stem cells through regulating the proliferation and differentiation of lung epithelial progenitor cells (Yuan et al., 2019). Moreover, FGF10-FGFR2B signaling can help to reduce elastin and fibronectin deposition (Chen S. et al., 2018). The reduction of FGF10 expression has been observed in IPF and has been identified as a cause of IPF progression (Chanda et al., 2016). Given these previous findings and our current results, we suggest that an effective treatment for silicosis may be the promotion of epithelial cell regeneration, either through activating FGF10-FGFR2B signaling (e.g., administration of exogenous FGF10) or through stem cell-based therapy, e.g., human umbilical mesenchymal stem cells (HUMSCs) transplantation, which can effectively reduce alveolar epithelial cells injury (Li et al., 2017) and attenuate silica-induced pulmonary fibrosis (Zhang et al., 2018).
Since there are currently no effective treatments for silicosis, it is of great importance to identify potential therapeutic targets and candidate drugs. In our analyses, we adopted bioinformatic methods to predict potential drugs that may attenuate silicosis. Several drugs were identified as potential candidates for further testing in silicosis (Supplementary Tables 24, 25). Plerixafor is a small molecule drug approved for patients with non-Hodgkin’s lymphoma and multiple myeloma. Previous research has shown that plerixafor inhibits CXCR4 on CD34+ cells (Uy et al., 2008). Combined with our findings that CXCR4 was increased in silicosis lungs and the CXCR4/CXCL12 axis may contribute to the enhancement of cytokines with profibrotic effects. Thus, plerixafor may target the inflammation induced by the CXCR4/CXCL12 axis. Retinoic acid has been reported to induce the expression of airway mucins (e.g., MUC5AC) and to induce mucociliary differentiation (Yoon et al., 1997; Gray et al., 2001; Choi et al., 2004). Since the epithelial injury and mucociliary dysfunction are predominant in silicosis lungs, it is also promising to investigate the effecacy of retinoic acid in the treatment of silicosis.
As the first cohort transcriptomics study of silicosis patients, we compared our results with previous studies conducted in lung tissues of animal models or silica-induced cell lines. Genes involved in immune-related pathways or processes, including leukocytes and lymphocytes activation, were strongly enriched for significantly altered genes, consistent with findings from a silicosis mice model (Chen J. et al., 2018). We also found significant up-regulation of activating transcription factor 3 (ATF3) as in macrophages exposed to silica (Chan et al., 2018). In silica-induced A549 lung epithelial cells, immune response and inflammatory pathways are induced (Chan et al., 2017), consistent with our findings. As might be expected, there were also notable differences between our findings and previous research. The overall comparison with previous report in a silicosis mice model (Chen J. et al., 2018) using CAT curves showed that very few genes (Supplementary Figure 1A) or pathways (Supplementary Figure 1B) showed consistency with our findings. Specifically, some relevant genes revealed by previous research did not show significant alterations in our results, such as the genes Ccl2 and Ccr2, which have been validated in mice silicosis models (Chen J. et al., 2018), and the DNA-binding protein inhibitor (ID) family genes, which were found to be relevant in the early responses of silicosis in epithelial cells (Chan et al., 2017). Considering differences between species and the focus of some previous research on early disease stages, our study conducted in human lungs improves our understanding of the true disease state and in particular the molecular changes underpinning silicosis.
We acknowledge limitations of our study. First, the study is limited by the small sample size. As lung transplantation is rarely conducted in silicosis patients in China, there is extremely limited availability of lung transplant tissue samples. Additionally, since the control samples were from the resected parts during lung transplantation, we were unable to decide the specific localization of sampling. The non-standardized localization of sampling from the lung may cause bias in process of comparison. To make the samples more comparable, we also sampled the silicosis lungs close to the edge of the lung. In addition, we avoided lymph nodes, large vessels, and atmospheric tubes in the sampling process. Even with these precautions, we cannot exclude the potential confounding effects caused by variability in sampling. Furthermore, the samples we used were obtained from advanced silicosis, therefore the inflammatory processes during the progression of disease may be underrepresented or underestimated. In the comparison with historical datasets (IPF and COPD), while patients in the three datasets all had advanced stage diseases and underwent lung transplantation, the different methods in sample processing and differential expression analysis may result in inaccurate conclusions, and further experimental validations are warranted. In summary, this study provides a comprehensive description of the molecular alterations in human silicosis lungs through comparative transcriptomics. The genes and pathways that characterize pathological changes in silicosis lungs may provide valuable resources and novel drug targets for silicosis.
Data Availability Statement
The RNA-seq data of this study has been deposited in the Genome Sequence Archive (GSA) (https://bigd.big.ac.cn/gsa-human/) with the Accession number: HRA000560.
Ethics Statement
The studies involving human participants were reviewed and approved by the Institutional Review Boards of Nanjing Medical University, and the Institute of Basic Medical Sciences, Chinese Academy of Medical Sciences. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
JC, JW, and CW designed the experiments. JW and CW obtained the funding. JC, DW, BL, and JZ collected the human lung samples. JP performed the bioinformatics analyses. JP and YLu prepared the original draft. JW, BL, ZC, XQ, MS, YLi, and ZL reviewed and edited the manuscript. All authors critically revised and approved the final manuscript.
Funding
This project was supported by the Chinese Academy of Medical Sciences Innovation Fund for Medical Sciences (CIFMS) (grant no. 2018-I2M-1-001), the Non-profit Central Research Institute Fund of Chinese Academy of Medical Sciences (grant nos. 2019RC330001 and 20180709), the State Key Laboratory Special Fund (grant no. 2060204), and the National Natural Science Foundation of China (grant no. 82090010).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We thank the Bioinformatics Centre of the Institute of Basic Medical Sciences, CAMS for providing supports of computing servers and MetaCore. We also thank Dr. Peiran Yang from the Institute of Basic Medical Sciences, CAMS for helping revise the manuscript.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2021.652901/full#supplementary-material
Footnotes
- ^ https://github.com/Pang-Junling/RNA-seq-analysis-for-human-silicosis-lungs
- ^ ftp://ftp.ensembl.org/pub/release-84/
- ^ https://CRAN.R-project.org/package=pheatmap
- ^ https://portal.genego.com/
- ^ https://portals.broadinstitute.org/cmap
- ^ https://go.drugbank.com/
References
Abler, L. L., Mansour, S. L., and Sun, X. (2009). Conditional gene inactivation reveals roles for Fgf10 and Fgfr2 in establishing a normal pattern of epithelial branching in the mouse lung. Dev. Dyn. 238, 1999–2013. doi: 10.1002/dvdy.22032
Association. (2018). Consensus of Chinese experts on pneumoconiosis treatment. J. Lab. Med. 35, 677–689.
Barnes, H., Goh, N. S. L., Leong, T. L., and Hoy, R. (2019). Silica-associated lung disease: an old-world exposure in modern industries. Respirology 24, 1165–1175. doi: 10.1111/resp.13695
Beamer, C. A., and Shepherd, D. M. (2013). Role of the aryl hydrocarbon receptor (AhR) in lung inflammation. Semin. Immunopathol. 35, 693–704. doi: 10.1007/s00281-013-0391-7
Benam, K. H., Vladar, E. K., Janssen, W. J., and Evans, C. M. (2018). Mucociliary defense: emerging cellular, molecular, and animal models. Ann. Am. Thoracic Soc. 15, S210–S215. doi: 10.1513/AnnalsATS.201806-439AW
Brilland, B., Beauvillain, C., Mazurkiewicz, G., Rucay, P., Roquelaure, Y., Tabiasco, J., et al. (2019). T cell dysregulation in non-silicotic silica exposed workers: a step toward immune tolerance breakdown. Front. Immunol. 10:2743. doi: 10.3389/fimmu.2019.02743
Bustamante-Marin, X. M., and Ostrowski, L. E. (2017). Cilia and mucociliary clearance. Cold Spring Harb. Perspect Biol. 9:a028241. doi: 10.1101/cshperspect.a028241
Cai, S. X., Liu, A. R., Chen, S., He, H. L., Chen, Q. H., Xu, J. Y., et al. (2016). The orphan receptor tyrosine kinase ROR2 facilitates MSCs to repair lung injury in ARDS animal model. Cell Transplant 25, 1561–1574. doi: 10.3727/096368915X689776
Chan, J. Y. W., Tsui, J. C. C., Law, P. T. W., So, W. K. W., Leung, D. Y. P., Sham, M. M. K., et al. (2017). Profiling of the silica-induced molecular events in lung epithelial cells using the RNA-Seq approach. J. Appl. Toxicol. 37, 1162–1173. doi: 10.1002/jat.3471
Chan, J. Y. W., Tsui, J. C. C., Law, P. T. W., So, W. K. W., Leung, D. Y. P., Sham, M. M. K., et al. (2018). RNA-Seq revealed ATF3-regulated inflammation induced by silica. Toxicology 393, 34–41. doi: 10.1016/j.tox.2017.11.001
Chanda, D., Kurundkar, A., Rangarajan, S., Locy, M., Bernard, K., Sharma, N. S., et al. (2016). Developmental reprogramming in mesenchymal stromal cells of human subjects with idiopathic pulmonary fibrosis. Sci. Rep. 6:37445. doi: 10.1038/srep37445
Chen, J., Yao, Y., Su, X., Shi, Y., Song, X., Xie, L., et al. (2018). Comparative RNA-Seq transcriptome analysis on silica induced pulmonary inflammation and fibrosis in mice silicosis model. J. Appl. Toxicol. 38, 773–782. doi: 10.1002/jat.3587
Chen, S., Cui, G., Peng, C., Lavin, M. F., Sun, X., Zhang, E., et al. (2018). Transplantation of adipose-derived mesenchymal stem cells attenuates pulmonary fibrosis of silicosis via anti-inflammatory and anti-apoptosis effects in rats. Stem Cell Res. Ther. 9:110.
Chen, Y., Li, C., Lu, Y., Zhuang, H., Gu, W., Liu, B., et al. (2017). IL-10-producing CD1d hi CD5 + regulatory B cells may play a critical role in modulating immune homeostasis in silicosis patients. Front. Immunol. 8:110. doi: 10.3389/fimmu.2017.00110
Choi, J. Y., Cho, K.-N., and Yoon, J.-H. (2004). All-trans retinoic acid induces mucociliary differentiation in a human cholesteatoma epithelial cell culture. Acta Otolaryngol. 124, 30–35.
Ehre, C., Worthington, E. N., Liesman, R. M., Grubb, B. R., Barbier, D., O’Neal, W. K., et al. (2012). Overexpressing mouse model demonstrates the protective role of Muc5ac in the lungs. Proc. Natl. Acad. Sci. U.S.A. 109, 16528–16533.
Fazen, L. E., Linde, B., and Redlich, C. A. (2020). Occupational lung diseases in the 21st century: the changing landscape and future challenges. Curr. Opin. Pulm Med. 26, 142–148.
GBD 2017 Disease and Injury Incidence and Prevalence Collaborators (2018). Global, regional, and national incidence, prevalence, and years lived with disability for 354 diseases and injuries for 195 countries and territories, 1990–2017: a systematic analysis for the Global Burden of Disease Study 2017. Lancet 392, 1789–1858. doi: 10.1016/S0140-6736(18)32279-7
Gray, T., Koo, J. S., and Nettesheim, P. (2001). Regulation of mucous differentiation and mucin gene expression in the tracheobronchial epithelium. Toxicology 160, 35–46.
Honda, K., Kimura, A., Dong, R. P., Tamai, H., Nagato, H., Nishimura, Y., et al. (1993). Immunogenetic analysis of silicosis in Japan. Am. J. Respir Cell Mol. Biol. 8, 106–111.
Hoy, R. F., and Chambers, D. C. (2020). Silica-related diseases in the modern world. Allergy 75, 2805–2817. doi: 10.1111/ALL.14202
Irizarry, R. A., Warren, D., Spencer, F., Kim, I. F., Biswal, S., Frank, B. C., et al. (2005). Multiple-laboratory comparison of microarray platforms. Nat. Methods 2, 345–350.
Jin, J., Togo, S., Kadoya, K., Tulafu, M., Namba, Y., Iwai, M., et al. (2019). Pirfenidone attenuates lung fibrotic fibroblast responses to transforming growth factor-beta1. Respir Res. 20:119. doi: 10.1186/s12931-019-1093-z
Jones, M. R., Dilai, S., Lingampally, A., Chao, C.-M., Danopoulos, S., Carraro, G., et al. (2019). A comprehensive analysis of fibroblast growth factor receptor 2b signaling on epithelial tip progenitor cells during early mouse lung branching morphogenesis. Front. Genet. 9:746. doi: 10.3389/fgene.2018.00746
Kage, H., and Borok, Z. (2012). EMT and interstitial lung disease: a mysterious relationship. Curr. Opin. Pulm Med. 18, 517–523.
Kim, D., Langmead, B., and Salzberg, S. L. (2015). HISAT: a fast spliced aligner with low memory requirements. Nat. Methods 12, 357–360. doi: 10.1038/nmeth.3317
Koeppen, M., McNamee, E. N., Brodsky, K. S., Aherne, C. M., Faigle, M., Downey, G. P., et al. (2013). Detrimental role of the airway mucin Muc5ac during ventilator-induced lung injury. Mucosal Immunol. 6, 762–775. doi: 10.1038/mi.2012.114
Kotton, D. N., and Morrisey, E. E. (2014). Lung regeneration: mechanisms, applications and emerging stem cell populations. Nat. Med. 20, 822–832. doi: 10.1038/nm.3642
Lawrence, M., Huber, W., Pages, H., Aboyoun, P., Carlson, M., Gentleman, R., et al. (2013). Software for computing and annotating genomic ranges. PLoS Comput. Biol. 9:e1003118. doi: 10.1371/journal.pcbi.1003118
Li, F., Xu, X., Geng, J., Wan, X., and Dai, H. (2020). The autocrine CXCR4/CXCL12 axis contributes to lung fibrosis through modulation of lung fibroblast activity. Exp. Ther. Med. 19, 1844–1854.
Li, X., Wang, Y., An, G., Liang, D., Zhu, Z., Lian, X., et al. (2017). Bone marrow mesenchymal stem cells attenuate silica-induced pulmonary fibrosis via paracrine mechanisms. Toxicol. Lett. 270, 96–107.
Liu, F., Dai, W., Li, C., Lu, X., Chen, Y., Weng, D., et al. (2016). Role of IL-10-producing regulatory B cells in modulating T-helper cell immune responses during silica-induced lung inflammation and fibrosis. Sci. Rep. 6:28911. doi: 10.1038/srep28911
Lopes-Pacheco, M., Bandeira, E., and Morales, M. M. (2016). Cell-based therapy for silicosis. Stem Cells Int. 2016, 1–9. doi: 10.1155/2016/5091838
Love, M. I., Huber, W., and Anders, S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15:550. doi: 10.1186/s13059-014-0550-8
Lu, Y., Liu, F., Li, C., Chen, Y., Weng, D., and Chen, J. (2017). IL-10-producing B cells suppress effector T cells activation and promote regulatory T cells in crystalline silica-induced inflammatory response in vitro. Mediators Inflamm. 2017:8415094. doi: 10.1155/2017/8415094
Migliaccio, C. T., Buford, M. C., Jessop, F., and Holian, A. (2008). The IL-4Ralpha pathway in macrophages and its potential role in silica-induced pulmonary fibrosis. J. Leukoc Biol. 83, 630–639.
Morrow, J. D., Zhou, X., Lao, T., Jiang, Z., DeMeo, D. L., Cho, M. H., et al. (2017). Functional interactors of three genome-wide association study genes are differentially expressed in severe chronic obstructive pulmonary disease lung tissue. Sci. Rep. 7:44232. doi: 10.1038/srep44232
Morse, C., Tabib, T., Sembrat, J., Buschur, K. L., Bittar, H. T., Valenzi, E., et al. (2019). Proliferating SPP1/MERTK-expressing macrophages in idiopathic pulmonary fibrosis. Eur. Respir J 54:1802441.
Peng, T., Frank, D. B., Kadzik, R. S., Morley, M. P., Rathi, K. S., Wang, T., et al. (2015). Hedgehog actively maintains adult lung quiescence and regulates repair and regeneration. Nature 526, 578–582.
Perkins, T. N., Peeters, P. M., Albrecht, C., Schins, R. P. F., Dentener, M. A., Mossman, B. T., et al. (2018). Crystalline silica alters Sulfatase-1 expression in rat lungs which influences hyper-proliferative and fibrogenic effects in human lung epithelial cells. Toxicol. Appl. Pharmacol. 348, 43–53.
Pertea, M., Kim, D., Pertea, G. M., Leek, J. T., and Salzberg, S. L. (2016). Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and Ballgown. Nat. Protoc. 11, 1650–1667. doi: 10.1038/nprot.2016.095
Piao, H.-L., Tao, Y., Zhu, R., Wang, S.-C., Tang, C.-L., Qiang Fu, M.-R. D., et al. (2012). The CXCL12/CXCR4 axis is involved in the maintenance of Th2 bias at the maternal/fetal interface in early human pregnancy. Cell Mol. Immunol. 9, 423–430.
Pollard, K. M. (2016). Silica, silicosis, and autoimmunity. Front. Immunol. 7:97. doi: 10.3389/fimmu.2016.00097
Posey, K. L., Coustry, F., and Hecht, J. T. (2018). Cartilage oligomeric matrix protein: COMPopathies and beyond. Matrix Biol. 7, 161–173. doi: 10.1016/j.matbio.2018.02.023
Ritchie, M. E., Phipson, B., Wu, D., Hu, Y., Law, C. W., Shi, W., et al. (2015). limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 43:e47. doi: 10.1093/nar/gkv007
Rose, C., Heinzerling, A., Patel, K., Sack, C., Wolff, J., Zell-Baran, L., et al. (2019). Severe silicosis in engineered stone fabrication workers — California, Colorado, Texas, and Washington, 2017–2019. MMWR Morb. Mortal Wkly Rep. 68, 813–818.
Ross, A. E., Marchionni, L., Vuica-Ross, M., Cheadle, C., Fan, J., Berman, D. M., et al. (2011). Gene expression pathways of high grade localized prostate cancer. Prostate 71, 1568–1577.
Roy, M. G., Livraghi-Butrico, A., Fletcher, A. A., McElwee, M. M., Evans, S. E., Boerner, R. M., et al. (2014). Muc5b is required for airway defence. Nature 505, 412–416. doi: 10.1038/nature12807
Sai, L., Yu, G., Bo, C., Zhang, Y., Du, Z., Li, C., et al. (2019). Profiling long non-coding RNA changes in silica-induced pulmonary fibrosis in rat. Toxicol. Lett. 310, 7–13. doi: 10.1016/j.toxlet.2019.04.003
Saito, A., Horie, M., and Nagase, T. (2018). TGF-beta signaling in lung health and disease. Int. J. Mol. Sci. 19:2460. doi: 10.3390/ijms19082460
Samsuzzaman, M., Uddin, M. S., Shah, M. A., and Mathew, B. (2019). Natural inhibitors on airway mucin: molecular insight into the therapeutic potential targeting MUC5AC expression and production. Life Sci. 231:116485. doi: 10.1016/j.lfs.2019.05.041
Singer, J. P., Chen, H., Phelan, T., Kukreja, J., Golden, J. A., and Blanc, P. D. (2012). Survival following lung transplantation for silicosis and other occupational lung diseases. Occup. Med. (Lond) 62, 134–137. doi: 10.1093/occmed/kqr171
Sivakumar, P., Thompson, J. R., Ammar, R., Porteous, M., McCoubrey, C., and Cantu, E. 3rd (2019). RNA sequencing of transplant-stage idiopathic pulmonary fibrosis lung reveals unique pathway regulation. ERJ Open Res. 5, 00117–2019. doi: 10.1183/23120541.00117-2019
So, T., Lee, S. W., and Croft, M. (2006). Tumor necrosis factor/tumor necrosis factor receptor family members that positively regulate immunity. Int. J. Hematol. 83, 1–11. doi: 10.1532/IJH97.05120
Stark, R., Grzelak, M., and Hadfield, J. (2019). RNA sequencing: the teenage years. Nat. Rev. Genet. 20, 631–656. doi: 10.1038/s41576-019-0150-2
Uy, G. L., Rettig, M. P., and Cashen, A. F. (2008). Plerixafor, a CXCR4 antagonist for the mobilization of hematopoietic stem cells. Expert Opin. Biol. Ther. 8, 1797–1804.
Walters, G. I. (2020). Occupational exposures and idiopathic pulmonary fibrosis. Curr. Opin. Allergy Clin. Immunol. 20, 103–111.
Wang, N., Keszei, M., Halibozek, P., Yigit, B., Engel, P., and Terhorst, C. (2016). Slamf6 negatively regulates autoimmunity. Clin. Immunol. 173, 19–26. doi: 10.1016/j.clim.2016.06.009
Wang, Y., Song, F., Zhu, J., Zhang, S., Yang, Y., Chen, T., et al. (2017). GSA: genome sequence archive ∗. Genomics Proteomics Bioinform. 15, 14–18. doi: 10.1016/j.gpb.2017.01.001
Wollin, L., Distler, J. H. W., Redente, E. F., Riches, D. W. H., Stowasser, S., Schlenker-Herceg, R., et al. (2019). Potential of nintedanib in treatment of progressive fibrosing interstitial lung diseases. Eur. Respir J 54:1900161. doi: 10.1183/13993003.00161-2019)
Yoon, J. H., Gray, T., Guzman, K., Koo, J. S., and Nettesheim, P. (1997). Regulation of the secretory phenotype of human airway epithelium by retinoic acid, triiodothyronine, and extracellular matrix. Am. J. Respir Cell Mol. Biol. 16, 724–731.
Yuan, T., Volckaert, T., Redente, E. F., Hopkins, S., Klinkhammer, K., Wasnick, R., et al. (2019). FGF10-FGFR2B signaling generates basal cells and drives alveolar epithelial regeneration by bronchial epithelial stem cells after lung injury. Stem Cell Rep. 12, 1041–1055. doi: 10.1016/j.stemcr.2019.04.003
Zhang, B., Xu, H., Zhang, Y., Yi, X., Zhang, G., Zhang, X., et al. (2019). Targeting the RAS axis alleviates silicotic fibrosis and Ang II-induced myofibroblast differentiation via inhibition of the hedgehog signaling pathway. Toxicol. Lett. 313, 30–41.
Zhang, E., Yang, Y., Chen, S., Peng, C., Lavin, M. F., Yeo, A. J., et al. (2018). Bone marrow mesenchymal stromal cells attenuate silica-induced pulmonary fibrosis potentially by attenuating Wnt/β-catenin signaling in rats. Stem Cell Res. Ther. 9:311.
Keywords: silicosis, human cohort, transcriptomics, pathological mechanisms, candidate targets
Citation: Pang J, Luo Y, Wei D, Cao Z, Qi X, Song M, Liu Y, Li Z, Zhang J, Li B, Chen J, Wang J and Wang C (2021) Comparative Transcriptome Analyses Reveal a Transcriptional Landscape of Human Silicosis Lungs and Provide Potential Strategies for Silicosis Treatment. Front. Genet. 12:652901. doi: 10.3389/fgene.2021.652901
Received: 13 January 2021; Accepted: 19 April 2021;
Published: 03 June 2021.
Edited by:
Stephen J. Bush, University of Oxford, United KingdomReviewed by:
Antje Prasse, Hannover Medical School, GermanyLuigi Marchionni, Cornell University, United States
Copyright © 2021 Pang, Luo, Wei, Cao, Qi, Song, Liu, Li, Zhang, Li, Chen, Wang and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Baicun Li, lbc19890303@126.com; Jingyu Chen, chenjy@wuxiph.com; Jing Wang, wangjing@ibms.pumc.edu.cn
†These authors have contributed equally to this work