ORIGINAL RESEARCH article

Front. Med., 20 October 2021

Sec. Respiratory Medicine

Volume 8 - 2021 | https://doi.org/10.3389/fmed.2021.744239

Identification of Genetic Signature Associated With Aging in Pulmonary Fibrosis

  • 1. Department of Respiratory and Critical Care Medicine, National Clinical Research Center of Respiratory Disease, Key Laboratory of Pulmonary Diseases of Health Ministry, Tongji Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China

  • 2. Department of Science, Western University, London, ON, Canada

  • 3. Department of Pulmonary and Critical Care Medicine, The First Affiliated Hospital of Sun Yat-sen University, Guangzhou, China

  • 4. Institute of Respiratory Diseases of Sun Yat-sen University, Guangzhou, China

Abstract

Background: Aging is a strong risk factor and an independent prognostic factor in idiopathic pulmonary fibrosis (IPF). In this study, we aimed to conduct a comprehensive analysis based on gene expression profiles for the role of aging in pulmonary fibrosis.

Method: Four datasets (GSE21411, GSE24206, GSE47460, and GSE101286) for patients with clinical IPF and one dataset for bleomycin (BLM)-induced pulmonary fibrosis (BIPF) mouse model (GSE123293) were obtained from Gene Expression Omnibus (GEO). According to different age ranges, both patients with IPF and BIPF mice were divided into young and aged groups. The differently expressed genes (DEGs) were systemically analyzed using Gene Ontology (GO) functional, Kyoto Encyclopedia of Genes and Genomes (KEGG), and hub genes analysis. Finally, we verified the role of age and core genes associated with age in vivo.

Results:Via the expression profile comparisons of aged and young patients with IPF, we identified 108 aging-associated DEGs, with 21 upregulated and 87 downregulated. The DEGs were associated with “response to glucocorticoid,” “response to corticosteroid,” and “rhythmic process” in GO biological process (BP). For KEGG analysis, the top three significantly enriched KEGG pathways of the DEGs included “IL-17 signaling pathway,” “Mineral absorption,” and “HIF-1-signaling pathway.” Through the comparisons of aged and young BIPF mice, a total number of 778 aging-associated DEGs were identified, with 453 genes increased and 325 genes decreased. For GO and KEGG analysis, the DEGs were enriched in extracellular matrix (ECM) and collagen metabolism. The common DEGs of patients with IPF and BIPF mice were enriched in the BP category, including “induction of bacterial agglutination,” “hyaluronan biosynthetic process,” and “positive regulation of heterotypic cell-cell adhesion.” We confirmed that aged BIPF mice developed more serious pulmonary fibrosis. Finally, the four aging-associated core genes (Slc2a3, Fga, Hp, and Thbs1) were verified in vivo.

Conclusion: This study provides new insights into the impact of aging on pulmonary fibrosis. We also identified four aging-associated core genes (Slc2a3, Fga, Hp, and Thbs1) related to the development of pulmonary fibrosis.

Introduction

The process of aging is featured by gradual functional impairments of tissue and organisms of the body, leading to reduced resilience to environmental damages and increased risks of occurrence of diseases and death (). Idiopathic pulmonary fibrosis (IPF) is an age-associated chronic lethal disease caused by an unknown reason (). It is commonly diagnosed in patients over 50 years old, and the incidence of this disease significantly increases in the population older than 50 years old. (). For this reason, better understandings of the pathophysiology of IPF and how the process of aging is related to the incidence of IPF are necessary for the development of new therapies for clinical treatments of this disease.

Despite the repeated proposals of the association between aging and IPF, the mechanism underlying this association is still unclear (, ). Lopez-Otin et al. enumerate nine tentative aging hallmarks: genomic instability, telomere attrition, epigenetic alterations, loss of proteostasis, deregulated nutrient sensing, mitochondrial dysfunction, cellular senescence, stem cell exhaustion, and altered intercellular communication (). Meanwhile, heterozygous mutations in four telomere-related genes are also found to associate with the development of pulmonary fibrosis (). Moreover, biomarkers of aging, including p21 and p16, are also found significantly upregulated in both epithelial cells and fibroblasts in lung tissue of patients with IPF as compared with that of healthy people (). In animal studies of IPF using the bleomycin (BLM)-induced pulmonary fibrosis (BIPF), aged BIPF mice failed to resolve fibrosis as compared with young mice (). However, despite the wide-recognized association between aging and IPF, its potential mechanism is unclear ().

In this study, to investigate the mechanisms underlying the association between pulmonary fibrosis and aging, using gene expression profiles of both patients with IPF and BIPF mice, we comprehensively analyzed the effects of aging in the development of pulmonary fibrosis and verified our findings in in vivo studies.

Method

Microarray Data Acquisition and Process

To identify the genes associated with aging in patients with IPF, four microarray datasets (GSE21411, GSE24206, GSE47460, and GSE101286) were downloaded from Gene Expression Omnibus (GEO, http://www.ncbi.nlm.nih.gov/geo) (). IPF was diagnosed according to the ATS/ERS criteria published (, ). Young patients with IPF were defined as age ≤ 55 years old, and aged patients with IPF were defined as age ≥70 years old. Gender- and age-matched healthy controls were gotten from GSE47460. In order to explore the genes associated with age in BIPF mice, we downloaded one microarray dataset (GSE123293) from GEO (). In this dataset, mice were divided into young (2–4 months old) and aged (18–20 months old) groups. Gender- and age-matched mouse controls were from GSE123293. Detailed information of GSE datasets was summarized in Table 1. Preprocessed expression matrix and probe annotation files of the four GSE datasets were obtained from the GEO repository. Probe annotations and the corresponding expression profiles from different datasets were mapped into a common gene expression list according to the method described by Shi et al. (). Probes obtained from different datasets were listed using official gene symbols. Multiple expression results of the certain gene were replaced by the median value of the expression results (, ). Log2 fold-change (FC) of all expression results was normalized using the Cross-platform Normalization (xpn) method by platforms.

Table 1

DatasetPlatformType of sampleSample (Young group)Sample (Aged group)Reference
GSE21411GPL570IPF lung tissue20()
GSE24206GPL570IPF lung tissue11()
GSE47460GPL6480IPF lung tissue1629()
GPL14550
Healthy lung tissue2323
GSE101286GPL6947IPF lung tissue03()
GSE123293GPL16570BIPF lung tissue33()
Healthy lung tissue33

Characteristics of studies included.

IPF, idiopathic pulmonary fibrosis; BIPF, bleomycin induced pulmonary fibrosis.

Differently expressed genes (DEGs) analysis, Gene Ontology (GO) functional analysis, and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis were performed.

To identify DEGs between young and aged groups, the limma R package (http://www.bioconductor.org/packages) was used to perform the negative binomial distribution method according to the absolute value of FC (>1.5) and false discovery rate (FDR) <0.05. According to the hypergeometric distribution algorithm, pathway enrichment analyses of the GO biological process (BP), molecular function (MF), and cell component (CC) were carried out by the “cluster-profler” R package. The same method was applied for the KEGG analysis. The cutoff value was p < 0.05. The Cytoscape (version 3.8.0) was used to select hub or core genes.

Animal

Specific pathogen-free (SPF) C57 BL/6N male mice (young mice, 2–4 months old; aged mice, 18–20 months old) were purchased from the HFK Bioscience company (Beijing, China). All animals were kept under the SPF environment at Tongji Medical College. All procedures for animal experiments were approved by the Animal Care and Use Committee of Tongji Hospital.

The animal model of pulmonary fibrosis was established by single intratracheal administration of BLM (2 mg/kg, MCE, USA) or normal saline (NS) as control. Mice were grouped as follows: NS (young mice intratracheal with NS), aged NS (aged mice intratracheal with NS), BLM (young mice intratracheal with BLM), and aged BLM (aged mice intratracheal with BLM). Mice were sacrificed 21 days after the establishment of the BIPF mice model. Bronchoalveolar lavage fluid (BALF) was carried out. Samples of lung tissue were collected. The left lung was paraffin embedded for hematoxylin-eosin (HE) staining and masson trichrome (MT) staining (). Lung hydroxyproline contents were measured by hydroxyproline colorimetric assay (Biovision, Milpitas, USA).

Quantitative Real-Time Polymerase Chain Reaction (RT-PCR)

Briefly, the total RNA of the lung tissue was extracted by RNAiso plus kit (TaKaRa) and reversely transcripted to cDNA. RT-PCR was carried out using SYBR Premix Ex Taq (Tarkara, Shiga, Japan). The relative mRNA expression was calculated using the 2−ΔΔCT method normalized to the level of β-actin. The primers of four core genes are listed in Supplementary Table 1.

Statistical Analysis

Data analysis was carried out by Rstudio (version 3.6.2) and GraphPad Prism 7. The microarray data were analyzed by different R packages. Results were expressed as means ± SEM. A two-way analysis of variance (ANOVA) followed by Tukey's multiple comparison post-test was used. p < 0.05 was considered as statistically significant.

Results

DEGs Associated With Age in IPF

To determine the core gene profiles related to age in patients with IPF, four microarray datasets were obtained from GEO. Information of the included studies is listed in Table 1; Supplementary Table 2. According to the age criterion, this study included 19 young patients with IPF and 33 aged patients with IPF. Gene expression profiles were compared between the two groups, generating 108 DEGs, with 21 upregulated and 87 downregulated (Figure 1A). For further understanding of the function related to these 108 DEGs, GO enrichment analysis, including BP, MF, and CC, was performed. In BP analysis, the DEGs were found to associate with “response to glucocorticoid,” “response to corticosteroid,” and “rhythmic process” (Figure 1B). In the CC category, the DEGs were related to the “integral component of postsynaptic membrane,” “intrinsic component of postsynaptic membrane,” and “integral component of synaptic membrane” (Figure 1C). Moreover, DEGs were enriched in the MF category related to “receptor ligand activity” and “growth factor activity” (Figure 1D). For KEGG pathway enrichment analysis, the top three significant KEGG pathways of the DEGs included “IL-17 signaling pathway,” “Mineral absorption,” and “HIF-1-signaling pathway” (Figure 1E). More information on IPF DEGs analysis was shown in Supplementary File 1. The same biological analysis was performed in human young healthy control (HC) and aged HC (Supplementary File 3; Supplementary Figure 1). Different from young and aged IPF, DEGs of young and aged HC were enriched in inflammation and immune.

Figure 1

DEGs Associated With Age in BIPF Mice

For further identification of the expression features associated with aging in BIPF mice, we obtained one dataset from GEO and compared the DEGs between aged and young BIPF mice. A total number of 778 DEGs were identified, with 453 genes increased and 325 genes decreased (Figure 2A). Slightly different from the results obtained from patients with IPF, the top five BP terms were enriched in “extracellular matrix organization,” “extracellular structure organization,” “cell chemotaxis,” “collagen metabolic process,” and “neutrophil chemotaxis” (Figure 2B). In the CC category, the DEGs were associated with “extracellular matrix,” “collagen-containing extracellular matrix,” and “extracellular matrix component” (Figure 2C). “Extracellular matrix component,” “glycosaminoglycan binding,” and “collagen binding” were the most important MF terms (Figure 2D). Furthermore, KEGG pathway analysis showed enrichments in “ECM-receptor interaction,” “Cytokine-cytokine receptor interaction,” and “Amoebiasis” (Figure 2E). Detailed information of BIPF DEGs analysis was displayed in Supplementary File 2. The same bioinformatical analysis was conducted in young mouse control (MC) and aged MC (Supplementary File 4; Supplementary Figure 2). Similar to HC age-associated DEGs, inflammation and immune were enriched in DEGs of young and aged MC.

Figure 2

Common DEGs Associated With Age in Pulmonary Fibrosis

To comprehensively explore the effects of aging in pulmonary fibrosis, we combined the results of patients with IPF and BIPF mice. As shown in the Venn diagram, eight common DEGs were found in both aged patients with IPF and aged BIPF mice (Figure 3A). The eight common DEGs were enriched in the BP category related to “induction of bacterial agglutination,” “hyaluronan biosynthetic process,” and “positive regulation of heterotypic cell-cell adhesion” (Figure 3B). In CC analysis, the eight common genes were associated with “extracellular space,” “extracellular exosome,” and “extracellular region” (Figure 3C). Subsequently, to confirm the aging-related core genes, we overlapped the eight common genes with the top 100 core DEGs in BIPF mice. In the results, the four aging-related core genes were identified: solute carrier family 2 member 3 (Slc2a3), fibrinogen alpha chain (Fga), haptoglobin (Hp), and thrombospondin 1 (Thbs1) (Figure 3D). The heatmap showed the four core genes expression in both patients with IPF and BIPF mice (Figures 3E,F). Detailed information on the four core genes is displayed in Table 2. There was no common gene among aged and young IPF DEGs, BIPF DEGs, HC DEGs, and MC DECs (Supplementary Figure 3; Supplementary File 5).

Figure 3

Table 2

HumanMouse
GeneIDlog2 FCP valueGeneIDlog2 FCP value
SLC2A36515−0.8010.00068Slc2a3205271.149134.74E−05
FGA2243−1.4620.00017Fga141611.421230.0001638
HP3240−1.1620.02032Hp154390.861720.000122
THBS17057−0.6110.03689Thbs121825−0.949970.038588

List of 4 core genes.

Human: aged IPF vs. young IPF; mouse: aged BLM vs. young BLM; ID, entrez ID; FC, fold change.

In vivo Validations of the Impact of Aging and Four Core Genes

Aging is an important risk factor in patients with IPF. In this study, we compared the fibrosis degree between young and aged mice. In histopathological analysis, we found that, with the increase of age, fibrosis appeared in lung tissue spontaneously (Figures 4A,B). Under the stimulation of BLM, more serious pulmonary fibrosis was developed in aged mice. This finding was consistent with the previous study (). For further confirmation of the expression changes in the four core genes (Slc2a3, Fga, Hp, and Thbs1) in NS, aged NS, BLM, and aged BLM groups, qRT-PCR verification was conducted. The results demonstrated that Slc2a3 and Fga were significantly increased with age in both NS and BLM-stimulated mice (Figures 5A,B). Meanwhile, the Hp and Thbs1 were found to significantly decreased with the increase of age (Figures 5C,D).

Figure 4

Figure 5

Discussion

Aging is a strong risk factor and an independent prognostic biomarker for progressive IPF (). However, a comprehensive analysis based on gene expression profiles for the investigation of the role played by aging in IPF was missing. In this study, using four human IPF datasets, we identified 108 altered aging-related DEGs in IPF. Similarly, through the analysis of a dataset of the BIPF mouse model, we identified 778 aging-related DEGs. Followed by that, both GO and KEGG analyses of the DEGs were conducted separately. The resulting DEGs were subsequently combined and generated eight common DEGs and four core genes. Finally, we confirmed that aging aggravated pulmonary fibrosis and validated the four aging-associated core genes in experimental pulmonary fibrosis.

In the present study, the 108 aging-associated DEGs of human IPF were found to be involved functionally in “response to glucocorticoid” and “response to corticosteroid” BP terms. Given the clinical practice at the time some of these samples were collected, it is likely some of these patients may have been receiving steroids, which could further influence this result. Although glucocorticoid showed no benefit to patients with IPF in several retrospective studies (, ), Zhang et al. implied that IPF with different gene types may benefit differently from glucocorticoid therapy (). Previous research reported that dexamethasone treatments targeting the macrophages in long significantly alleviated pulmonary fibrosis through the regulation of the immune microenvironment (). Generally, more comprehensive experimental designs might apply for exploring the role of glucocorticoids in IPF. Postsynaptic components are involved in most top five enriched CC terms, indicating an important role of the nerve system in aged patients with IPF. A longitudinal IPF cohort demonstrated that patients with IPF treated with α1 adrenoreceptor blockade displayed improvements in survival rates (). Moreover, vagotomy attenuated symptoms of pulmonary fibrosis in BIPF mice by enhancing fibrogenic factors and fibrogenic cells (). Growth factor activity was one significantly enriched MF term. Genes related to growth factor activity were downregulated in aged-related disease (, ). Nutrient metabolism imbalance was enriched in KEGG analysis. Loss of pyruvate kinase M2 resulted in cone degeneration in an age-dependent manner (). Moreover, mitochondrial energy production reduction and cellular redox hemostasis loss played important roles in hearing loss associated with aging ().

Different from DEGs of IPF, the significant GO enrichment of DEGs in the mouse model of BIPF was extracellular matrix (ECM) and collagen metabolism. Persistent ECM accumulation was considered to participate in the development of several diseases, including IPF, causing difficulties in reversing the disease progression of such disease (). The ECM of aged lungs was featured by the upregulated transcriptional level of fibronectin extracellular domain A (Fn-EDA) and matrix metalloproteinases (MMPs), including MMP-2 and MMP-9 (). Collagen contents were found higher in adult mice; however, its transcriptional level was found higher in young mice. This phenomenon was caused by dysregulations in collagen degradation (). Despite the wide recognition of the association between ECM dysregulations and the process of aging, its underlying mechanisms were still unclear.

Age-associated DEGs of IPF and BIPF exhibited differences, while age-associated DEGs of HC and MC overlapped more. Inflammation and immunomodulation were enriched in age-associated DEGs of HC and MC. Age-related chronic inflammation is a major contributor to diseases with advancing age (). Extensive evidence indicates that persistent low-grade inflammation is present in the aged population and that age-associated inflammation occurs in the lungs (, ). BALF collected from aged mice and aged humans contains increased levels of IL-6, TNF, and complement components, reflective of a disordered innate immunity (). Although collagen and matrix synthesis are important in aged pulmonary fibrosis, the disordered inflammation and immune response cannot be ignored in aged pulmonary fibrosis (Figure 2B).

The GO terms of the common eight genes were enriched in the hyaluronan (HA) biosynthetic process. Hp, one of the four core genes, produced important effects in the process of lung injuries and repairments (). Extracellular matrix Hp on type 2 alveolar epithelial cells (AEC2s) was critical for AEC2 renewal, repairments of lung injuries, and limitations of the development of fibrosis (). Thbs1 is associated with aging and promoted matrix homeostasis through the interaction with collagen, lysyl oxidase precursors, and sites of collagen cross-linking (, ). The relationship between Thbs1 and aged pulmonary fibrosis was worth studying. Slc2a3, mediating the uptake of various monosaccharides, was reported to be associated with pulmonary fibrosis. Aging altered the uptakes of glucose and the transcription profile of Slc2a3/Slc2a4 before and after spinal cord injury, which is potentially related to the inflammation level (42). Fga was cleaved by the protease thrombin to produce monomers. Subsequently, the monomers polymerized with fibrinogen beta (Fgb) and fibrinogen gamma (Fgg), forming the insoluble fibrin matrix. The level of Fga was higher in patients with IPF compared with that in healthy controls (43). Moreover, previous research found that the lack of a C-terminal fragment of Fga preceded the progress of fibrosis in patients with liver diseases (44). Most importantly, differences in the expression profiles of the four core genes were found in both patients with IPF and BIPF mice. In order to clarify that the expression changes of the four core genes were due to aging in the context of disease or aging alone, we analyzed the age-associated DEGs in HC and MC. Comprehensively, analysis of our research results, the bioinformatics analysis results, and the expression change of four core genes were due to aging in the context of pulmonary fibrosis.

Therefore, considering the crucial roles of these four core genes on the basis of this study and previous studies above, further studies may be focused on exploring the mechanisms under pulmonary fibrosis in the context of aging. There were still several limitations in the study. Firstly, we could not assess the role of other aging marks in pulmonary fibrosis, such as telomere attrition, epigenetic alterations, and genomic instability. They need more professional research methods, whereas we only explored the role of age in pulmonary fibrosis from the perspective of gene transcription. Secondly, the expression of the four core genes needs to be verified in humans. Thirdly, while our study was able to identify a number of novel pathways in aged IPF, it was limited by the available clinical data (such as the history of glucocorticoid use) and biases inherent in enrichment analyses. Fourthly, some gene information is lost in the process of data sets integration because it cannot be repeated in different datasets. Furthermore, more and larger-scale studies should be conducted.

It was widely recognized that aging was a dependent risk factor in pulmonary fibrosis. As shown by our study, aging-aggravated BLM induced pulmonary fibrosis in BIPF mice. The pathogenesis of young and aged pulmonary fibrosis was not exactly the same from the perspective of gene expression profiles. Four aging-associated core genes (Slc2a3, Fga, Hp, and Thbs1) were related to the development of pulmonary fibrosis.

Funding

This study was supported by the National Natural Science Foundation of China (Nos. 81973986, 82170049, 82070032, and 81800041), Health Research Fund of Wuhan (No. WX21Q07), Health and family planning research project of Hubei (No. WJ2019M116), and the National Major Science and Technology Project for the Control and Prevention of Major Infectious Diseases of China (2017ZX10103004).

Publisher's Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Statements

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

Ethics statement

The animal study was reviewed and approved by Animal Care and Use Committee of Tongji Hospital.

Author contributions

YL and JC: conceptualization, investigation, formal analysis, and writing original draft. SW: conceptualization, formal analysis, and writing original draft. ZT and MW: formal analysis. YF: formal analysis and writing original draft. JZ and JX: conceptualization and supervision. KT: conceptualization, investigation, and formal analysis. All authors contributed to the article and approved the submitted version.

Acknowledgments

We thank all the patients with IPF and healthy controls involved in our study.

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.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmed.2021.744239/full#supplementary-material

    Abbreviations

  • IPF

    idiopathic pulmonary fibrosis

  • BLM

    bleomycin

  • BIPF

    bleomycin (BLM)-induced pulmonary fibrosis

  • GEO

    gene expression omnibus

  • FC

    fold change

  • DEGs

    differently expressed genes

  • GO

    gene ontology

  • KEGG

    Kyoto Encyclopedia of Genes and Genomes

  • BP

    biological process

  • CC

    cell component

  • MF

    molecular function

  • SPF

    specific pathogen-free

  • NS

    normal saline

  • BALF

    bronchoalveolar lavage fluid

  • HE

    hematoxylin-eosin

  • MT

    masson trichrome

  • RT-PCR

    real-time polymerase chain reaction

  • ANOVA

    one-way analysis of variance

  • HC

    healthy control

  • MC

    mouse control

  • Slc2a3

    solute carrier family 2 member 3

  • Fga

    fibrinogen alpha chain

  • Hp

    haptoglobin

  • Thbs1

    thrombospondin 1

  • ECM

    extracellular matrix

  • Fn-EDA

    fibronectin extracellular domain A

  • MMPs

    matrix metalloproteinases

  • AEC2

    type 2 alveolar epithelial cells

  • Fgb

    fibrinogen beta

  • Fgg

    fibrinogen gamma.

References

Summary

Keywords

aging, pulmonary fibrosis, genetic signature, idiopathic pulmonary fibrosis, glucocorticoids

Citation

Lu Y, Chen J, Wang S, Tian Z, Fan Y, Wang M, Zhao J, Tang K and Xie J (2021) Identification of Genetic Signature Associated With Aging in Pulmonary Fibrosis. Front. Med. 8:744239. doi: 10.3389/fmed.2021.744239

Received

22 July 2021

Accepted

20 September 2021

Published

20 October 2021

Volume

8 - 2021

Edited by

Marta Bueno, University of Pittsburgh, United States

Reviewed by

Eleanor Valenzi, University of Pittsburgh, United States; Yuli Huang, Southern Medical University, China; Mingyi Wang, National Institutes of Health (NIH), United States

Updates

Copyright

*Correspondence: Jianping Zhao Kun Tang Jungang Xie

†These authors share first authorship

This article was submitted to Pulmonary Medicine, a section of the journal Frontiers in Medicine

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics