- Open Access
Accumulation of damaged mitochondria in alveolar macrophages with reduced OXPHOS related gene expression in IPF
Respiratory Research volume 20, Article number: 264 (2019)
Impaired mitochondria homeostasis and function are established hallmarks of aging and increasing evidence suggests a link with lung fibrosis. Mitochondria homeostasis may be also affected in alveolar macrophages (AMs) in idiopathic pulmonary fibrosis (IPF). In this study, we used bronchoalveolar lavage (BAL), a tool for both clinical and research purposes, and a rich source of AMs.
BAL samples were examined from 52 patients with IPF and 19 healthy individuals. Measurements of mitochondria reactive oxygen species (mtROS), mitochondria morphology and related gene expression were performed. Additionally, autophagy and mitophagy levels were analysed.
Mitochondria in AMs from IPF patients had prominent morphological defects and impaired transcription paralleled to a significant reduction of mitochondria homeostasis regulators PINK1, PARK2 and NRF1. mtROS, was significantly higher in IPF and associated with reduced expression of mitochondria-encoded oxidative phosphorylation (OXPHOS) genes. Age and decline in lung function correlated with higher mtROS levels. Augmentation of damaged, oxidised mitochondria in IPF AMs however was not coupled to increased macroautophagy and mitophagy, central processes in the maintenance of healthy mitochondria levels.
Our results suggest a perturbation of mitochondria homeostasis in alveolar macrophages in IPF.
Idiopathic Pulmonary Fibrosis (IPF) is a chronic lung disease, characterized by poor prognosis, primarily affecting older individuals . IPF pathogenesis is attributed to an accelerated apoptosis of Type-II pneumocytes coupled to an aberrant wound healing process, involving defective re-epithelization, abnormal fibroblast/myofibroblast apoptosis and excessive accumulation of extracellular matrix components . Age associated changes, including telomere attrition, defective autophagy, mitochondrial dysfunction and cell senescence, are linked to IPF pathogenesis .
Mitochondria drive ATP production through oxidative phosphorylation (OXPHOS) and they represent the major generator of reactive oxygen species (ROS). In the last decade, it has become evident that cellular metabolism is intertwined with immune responses while, mitochondria activity largely shapes immune cell functions . Mitochondrial function declines with age, resulting in morphological alterations, loss of cristae, reduced mitochondria respiratory chain (MRC) function and increased ROS production which is associated with mtDNA damage, fueling a vicious cycle of mitochondria impairment . Mitochondrial homeostasis is essential for cellular health, and is tightly regulated by biogenesis and degradation/recycling . PTEN-induced kinase 1 (PINK1) and E3 ubiquitin ligase (Parkin - encoded by PARK2 gene) represent key effector proteins of mitophagy, a selective form of macroautophagy, participating in the clearance of oxidized and damaged mitochondria .
Accumulation of damaged mitochondria is associated with IPF and the aging lung [3, 8, 9]. PINK1 deficiency and insufficient mitophagy in alveolar epithelial cells leads to increased mitochondrial damage and progression of fibrosis . Conversely, PINK1 upregulation in mice contributes to elimination of fibrosis . PINK1 decrease in fibrotic foci in IPF, and TGF-β treatment of lung fibroblasts lowers PINK1 and increases ROS levels . PARK2 knockdown induces mitochondrial damage, increases mtROS production and promotes fibrosis through enhancement of apoptotic resistant phenotype of fibroblasts and myofibroblast differentiation . MtROS is linked to mtDNA damage and evidence from a recent study in IPF showed a higher mutation rate in mitochondrial tRNA, probably affecting intra-mitochondrial translation . Importantly, a widely used antidiabetic medication, metformin, ameliorated lung fibrosis in mice through beneficial effects on mitochondria, highlighting their pathogenetic role in pulmonary fibrosis .
Lung macrophages are a fluctuating and diverse population with a central role in inflammation, host defense, and tissue homeostasis . It is well documented that macrophages within the lung orchestrate the downstream progression and maintenance of fibrosis . A substantial expansion of macrophages in distal lung tissue is observed in IPF and a profibrotic macrophage population that appears in the lung during the development of fibrosis is currently under investigation [18, 19]. ROS production is required for the alternative activation of macrophages and can be attenuated by ROS inhibitors including N-acetyl-cysteine . Interestingly, generation of mtROS from AMs has been causatively related to the development of asbestos-induced pulmonary fibrosis . Mitochondria homeostasis in IPF AMs is not well characterized and in this study we sought to evaluate mitochondrial ROS levels, morphology, and mitophagy in AMs isolated from BAL in IPF patients.
Materials and methods
Seventy one (71) patients were prospectively enrolled in the Department of Thoracic Medicine, University Hospital of Heraklion, Crete, Greece from 2014 to 2018. Fifty two  patients with IPF and nineteen  healthy controls underwent bronchoscopy and Bronchoalveolar Lavage Fluid (BAL) was obtained.
IPF group: The diagnosis of IPF was based on ATS/ERS criteria or on multidisciplinary discussion . Bronchoscopy was performed in treatment naive IPF patients as a part of the diagnostic approach.
All patients were evaluated within 1 month from bronchoscopy, with complete pulmonary function tests (PFTs). Lung volumes (forced expiratory volume in 1 sec – FEV1, forced vital capacity – FVC), using the helium-dilution technique and diffusion capacity (DLco, corrected for haemoglobin) using the single breath technique. The computerized system (Jaeger 2.12; MasterLab, Würzburg, Germany) was used and predicted values were obtained from the standardized lung function testing of the European Coal and Steel Community, Luxembourg (1993).
Control group: Control subjects were patients undergoing bronchoscopy for the investigation of haemoptysis and BAL sampling was performed at least 6 weeks following hospitalization. No signs of infection, ILD abnormalities or putative malignant lesions in the CT scans were observed. Furthermore, during bronchoscopy, no malignant lesions in the airways were observed, or abnormal mucus secretion was noted. BAL cellularity and differential BAL cell count were within the normal range. Additionally, BAL samples testing positive for bacteria growth were excluded. Since controls were healthy subjects no PFTs were performed.
Patients were classified as non-smokers, current smokers or former smokers (defined as having smoked a minimum of one cigarette a day for a minimum of 1 year, stopping at least 6 months before presentation).
Demographic data, pulmonary function tests, and BAL differential cell count are summarised in Table 1.
BAL cell isolation and determination of cellular composition
Bronchoalveolar Lavage (BAL) was obtained as previously reported . In brief, BAL was obtained from all patients at room temperature. A flexible bronchoscope was wedged into a sub-segmental bronchus of a predetermined region of interest based on radiographical findings. A BAL technique was performed by instilling a total of 180 ml of normal saline in 60-mL aliquots, each retrieved by low suction.
The BAL samples were subsequently kept on ice and were processed within 2 h of collection. Samples were filtered through sterile 70 nm cell strainers (BD) and centrifuged at 1500 rpm for 5 min at 4οC. Cell pellets were washed and re-suspended with cold PBS. Total cell count and cell viability were subsequently assessed using Trypan blue (ICN). Differential cell population count was analysed following May-Grunewald-Giemsa staining as previously described .
BAL cell isolation and AM culture
In order to study alveolar macrophages, by immune fluorescence or western blot, we enriched their population by pre-culture on glass slides or plastic dishes respectively, followed by washing of non-adherent cells such as lymphocytes, eosinophils and neutrophils, minor components of the BAL cell population. The remaining attached cell population comprises of macrophages and monocytes from the alveolar space. 0.5 million freshly isolated BAL cells were cultured in 24 well plates in DMEM (Biosera) growth media supplemented with 2% FCS (Biosera) and penicillin-streptomycin in a humidified incubator at 37 °C containing 5% CO2 for 1 h, with subsequent washes to remove non-adherent cells, allowing enrichment for alveolar macrophages/ monocytes.
BAL cell surface marker expression analysis
0.5 million freshly isolated BAL cells were resuspended in PBS, 2% FCS and 2 mM EDTA buffer and labelled with pan-leucocyte marker CD45-FITC (#IQP-124F, IQProducts), macrophage/DC marker CD11c-PE-Cy5 (#301610, Biolegend), or CD45-FITC, monocyte marker CD14-PE-Cy5 (#A07765, IOTest, Beckman Coulter). CD163-PE-Cy7 (#333614, Biolegend) and appropriate Isotype controls. Data were acquired with Beckman Coulter flow cytometer and analyzed with FlowJo 8.7 (Treestar, Ashland, OR).
Mitochondrial ROS measurement
Mitochondrial ROS was measured by MitoSOX™Red (Invitrogen) staining. MitoSOX™red is targeted to mitochondria in live cells and is readily oxidized by superoxide reactive oxygen species. Initially, the staining of oxidized mitochondria with MitoSOX following treatment with hydrogen peroxide (H2O2) was assessed. A549 cells, a non-small cell lung cancer (NSCLC)-derived human alveolar epithelial cell line, were treated with 1000 μM H2O2 and the characteristic staining for the mitochondrial network in the cytoplasm was observed which was further enhanced upon H202 treatment and mitochondria oxidation (Additional file 1: Figure S1a). Secondly, flow cytometry analysis of MitoSOX staining was tested on of THP1 human monocytes. THP1 cells showed a strong fluorescent MitoSOX signal (Additional file 1: Figure S1b) which was greatly enhanced following 1000 μM H2O2 treatment. Necrotic and apoptotic cells were expected to contain oxidized mitochondria and since H2O2 treatment is known to induce apoptosis, thus Propidium Iodide (PI) staining was analysed in parallel which allowed for the exclusion of the MitoSOX positive necrotic/apoptotic population from the analysis of the alive MitoSOX positive cells (Additional file 1: Figure S1b). A high proportion of MitoSOX positive cells upon H2O2 treatment was present despite the exclusion of the PI negative cell population (Additional file 1: Figure S1b) and the same PI positive population exclusion strategy was used subsequently in human BAL cells.
0.5 million freshly isolated BAL cells were cells were resuspended in RPMI-1640 supplemented with 2% FCS were stained with MitoSOX Red at a final concentration of 5 μM and CD45-FITC for 10 min at 37 °C. For FL-2, MitoSOX staining quantification, alveolar macrophages populations was selected according to high side scatter SSChigh and CD45+(Additional file 1: Figure S1c and d). Propidium Iodide (PI) was added to cells not previously stained with MitoSOX at a final concentration of 1 ng/ml, for 5 min immediately before flow cytometry analysis, for the detection of necrotic/apoptotic cells. A CD45-FITC stained control for each patient sample served as the unstained control sample for subsequent quantifications. PI staining allowed the exclusion of apoptotic/necrotic cells. Despite the PI exclusion gating strategy, a small percentage of Propidium Iodide (PI) positive cells was still among the MitoSOX positive analysed population and it was subtracted from the MitoSOX positive cell percentage during quantifications shown in Fig. 2b.
The percentage of MitoSOX positive cells was determined by the percent of cells showing FL-2 fluorescence higher than the unstained control, followed by subtraction of the FL-2 PI positive percentage of cells. Relative mean fluorescence intensity (MFI) was calculated by normalizing the MFI of the FL-2 channel/MitoSOX positive cells by the MFI of the FL-2 channel of the unstained cells since patient samples displayed wide ranges of autofluorescence. Data were acquired with Beckman Coulter flow cytometer and analyzed with FlowJo 8.7 (Treestar, Ashland, OR).
RNA extraction and mRNA expression
1–1.5 million cells were centrifuged and cell pellets were homogenised in TriReagent™(MBL) for total RNA, followed by storage at -80οC. Total RNA extraction and cDNA synthesis were performed as previously described . Probe and primer sequences are summarized in Additional file 2: Table S1. GAPDH levels were used as endogenous control for the normalization of mRNA expression levels in BAL samples. Gene expression analysis was performed following incorporation of relative expression values in average (duplicates) normalized by GAPDH. Relative expression values for the patient cohort were calculated by 2-ΔΔCt method, where ΔΔCt = (sample Ct GOI-sample Ct GAPDH)-(Calibrator Ct GOI- Calibrator Ct GAPDH), GOI = gene of interest, calibrator = mean of all Cts.
mtDNA/ gDNA ratio determination
Total DNA isolation from 104 to 105 BAL cells was performed using NucleoSpin Tissue (Macherey-Nagel) according to the manufacturer’s protocol. To quantify the mtDNA/gDNA ratio, qPCR was used to quantify the mitochondria DNA encoded genes ND1 and ND5 and the haemoglobin B1 gene. Primer sequences are included in Additional file 2: Table S1. Relative amounts of mitochondria DNA to HGB1 genomic sequence was calculated according to 2-ΔΔCt method, where ΔΔCt = (sample Ct GOI-sample Ct HGB1)-(Calibrator Ct GOI- Calibrator Ct HGB1), GOI = gene of interest ie ND1 or ND5, calibrator = mean of all Cts.
Western blot analysis
1–1.5 million cells were centrifuged and cell pellets were homogenised in RIPA buffer (Invitrogen) containing protease and phosphatase inhibitors, Pierce), followed by storage at -80οC. 40–60 μg Total protein lysates of BAL samples were separated in 12% SDS-PAGE, transferred to 0.45 nm nitrocellulose membrane (Biorad), followed by detection of PINK1[anti-PINK1, mouse monoclonal antibody, (Novus Biologicals)], TOMM20 [anti-TOMM20, rabbit polyclonal (Abnova)], p62 [anti-p62 mouse monoclonal antibody (MBL)], LC3 [anti-LC3 mouse monoclonal antibody (Abgent)] and b-actin [anti-b-actin mouse monoclonal antibody (Sigma)]. Appropriate anti-mouse HRP conjugated secondary antibody (Chemicon) was used and immunodetection was performed with enhanced chemiluminescence reagent Luminata™ (Millipore). Bands were visualised with the ChemiDocXRS+ system (Biorad) and densitometry analyses were performed using Image Lab™ software (Biorad).
Permeabilisation buffer (0.5%FBS, 0.2%Triton in TBS) was added for 10 min followed by blocking buffer (0.5%FBS, 0.1%Triton, 2 mg/ml BSA in TBS) for 10 min at room temperature. Primary antibody incubations were carried out for 60 min at room temperature followed by washing with TBS. Secondary antibody incubations were performed for 30 min followed by washing with TBS. ToPro or DAPI was added for nuclear staining. Finally, TrueBlack (Biotium) was added for 30 s so as to eliminate autofluorescence signal due to intracellular lipofuscin granules. Images were taken by Leica confocal microscope. Mean intensity of protein levels/cell was performed with ImageJ v2.0.0-rc-69/1.52i. Colocalization of mitochondria with autophagosomes or lysosomes was analysed by coloc2 plugin of Image J software.
Lysosome inhibition for autophagy and mitophagy flux measurements
BAL cells were plated on glass slides or plastic dishes for 1 h as described above. Following removal of non-adherent cells, AMs were treated with 0.1 mM chloroquine, for 30 min and fixed with 4% formaldehyde for 20 min followed by washing with PBS for immunofluorescence analysis or lysed in RIPA buffer for western blot analyses.
Data were analyzed using Prism 5 (Graph Pad) software and comparisons were made by t test, Mann-Whitney test, or chi-square test, as appropriate. Spearman’s correlation coefficient (r) analysis measured the association between two variables. A p value less than 0.05 was considered statistically significant (*p < 0.05, **p < 0.01, ***p < 0.001).
PINK1 and PARK2 expression is reduced in IPF BAL cells
The downregulation of both PINK1 and PARK2, two key molecules in mitochondrial homeostasis has been reported in IPF [10, 13]. Analysis of total PINK1protein levels in whole BAL cells showed a significant decrease of PINK1 in IPF (Fig. 1a and b). Fluorescence intensity quantification of immune-stained BAL AM cultures also showed lower levels of PINK1 in IPF (Fig. 1c-e). Consistently with protein levels, PINK1 mRNA in total BAL cell lysates were significantly lower in IPF (Fig. 1f). Additionally, PARK2 mRNA was also significantly lower in IPF (Fig. 1g) and positively correlated with PINK1 mRNA levels (p < 0.001, r = 0.5).
Increased levels of mtROS in IPF alveolar macrophages
The levels of mtROS were measured in fresh, BAL AMs (CD45+FShighSShighcells) by flow cytometry using MitoSOX™red superoxide indicator (Fig. 2a). The CD45+FShighSShigh population consisted of CD11c+ and/or CD14+ macrophages as determined in parallel (Additional file 1: Figure S1c). An increase in mtROS levels was demonstrated in IPF AMs, as shown by higher mean fluorescence intensity (MFI) (Fig. 2b, Additional file 1: Figure S1d) and percentage of positive cells (Fig. 2c). In parallel, analysis of cell viability, using propidium iodide staining (Fig. 2d), revealed no difference between IPF and control groups suggesting that the elevated mtROS in IPF were not due to an increase of apoptotic/necrotic cells.
Importantly, mtROS at the time of diagnosis showed a mild correlation with Composite Physiologic Index (CPI) (p = 0.04, spearman r = 0.33) (Fig. 2e), lower diffusion capacity for carbon monoxide (DLco%) (p = 0.03, spearman r = − 0.35)(Fig. 2f) and total lung capacity (TLC%) (p = 0.03, spearman r = − 0.35)(Fig. 2g). Furthermore, the levels of mtROS significantly correlated with patients’ age (r = 0.37, p = 0.01) (Fig. 2h) and not smoking habits.
CD163 expression is upregulated IPF AMs
MtROS is required for alternative activation of macrophages. We therefore examined the surface expression of scavenger receptor CD163, a frequently used marker of alternatively activated macrophages  by flow cytometry in the CD45+FShighSShighCD11c+AM population of the BAL cells. CD163 expression was significantly higher in IPF relative to controls (Fig. 3a and b). However, no correlation between CD163 and the levels of mtROS in the IPF cohort (r mtROS-MFI = -0.11, r mtROS-% = − 0.27, pNS) was observed. TGFb mRNA expression was not higher in IPF BAL cells (Fig. 3c), in contrast to Collagen1a levels (Fig. 3d) as previously described .
Mitochondria encoded OXPHOS related gene expression is lower in IPF BAL cells
Persistently elevated levels of mtROS observed during ageing and chronic lung diseases such as IPF may drive mitochondria damage and impair their function. Since alternatively activated macrophages rely on OXPHOS for ATP production, we analyzed the transcription of mt-encoded gene members of OXPHOS in particular. Transcription of MRC complex I and V members, NADH dehydrogenase subunit 1 (ND1) and ATP synthase 6 (ATP6) respectively, was significantly lower in IPF BAL cells (Fig. 4a and b). Furthermore, the mt-encoded 12 s ribosomal RNA required for the intra-mitochondrial translation was also significantly decreased in IPF (Fig. 4c). Additionally, the expression of nuclear respiratory factor 1(NRF1), that activates the expression of mitochondrial transcription factor A (TFAM) and in turn mt-OXPHOS gene expression, was also significantly lower in IPF BAL cells (Fig. 4d). Importantly, the levels of mtROS were associated with a reduction of MRC transcription of both complex I (ND1r = − 0.37, p = 0.03) and complex V (ATP6 r = − 0.4,p = 0.01).
Accumulation of dysmorphic mitochondria in IPF AMs
Mitochondria morphology in IPF AMs was evaluated using transmission electron microscopy (TEM). Ultrastructural qualitative analysis revealed that mitochondria in IPF were dysmorphic and contained disorganized cristae (Fig. 5a). Quantitative morphometry in a subset of samples showed a significant increase in the frequency of swollen mitochondria with higher median diameter size ascompared to controls (Fig. 5b and c). In parallel, in IPF, mitochondria were shorteras shown by their significantly decreased length (Fig. 5d and e).
Mitochondria mass was not altered in IPF AMs
Previous studies suggested that mitochondria mass increases in IPF lungs [10, 26]. We therefore determined the levels of mtDNA and mitochondria protein marker, Translocase Of Outer Mitochondrial Membrane 20 (TOMM20). The mitochondria encoded genes, ND1 and ND5, relative to the hemoglobin 1 (HGB1) genomic locus were analyzed by qPCR in whole BAL cell lysates. No differences between controls and IPF patients (Additional file 3: Figure S2a and b), or mtDNA correlation with mtROS were observed. Additionally, TOMM20 levels, measured either by western blot in whole BAL cell lysates or immune fluorescence in AMs, were similar between the two groups (Additional file 3: Figure S2c-g).
LC3B II and autophagy levels in IPF alveolar macrophages
In view of the accumulation of oxidized and dysmorphic mitochondria in AMs we subsequently examined commonly used markers of autophagy and mitophagy. Initially, we assessed endogenous basal expression of autophagy markers; microtubule-associated protein light chain 3B (LC3B), forms I and II, and the selective autophagy receptor p62/sequestosome1 (SQSTM1) protein in BAL cells. Western Blot analysis from IPF and control groups showed no significant differences in the cytosolic form of LC3B (LC3B I), the lipidated form present mostly on autophagosomes (LC3B II), the LC3B II/I ratio or the p62 protein (Fig. 6a-e and Additional file 4: Figure S3a-b). mRNA levels of p62 and Beclin 1 were also measured and no significant differences were observed (Additional file 4: Figure S3c-d).
Autophagy flux was measured following treatment with chloroquine, a lysosomal acidification inhibitor that enables the estimation of the turnover of proteins participating in the degradation of autophagy cargo. Similarly to the basal levels, we did not observe significant differences in the accumulation of autophagy markers LC3B II and p62 in IPF relative to controls (Fig. 6f-h). LC3B levels were also analysed in AM cultures, following chloroquine treatment. In agreement with the western blot analysis, in addition to LC3B I, LC3B II puncta were readily observed in control and IPF AMs, in untreated and following chloroquine treatment (Additional file 4: Figure S3e).
LC3B/PINK1 dependent mitophagy levels were similar in IPF and controls
As an indication of mitophagy induction and mitochondria degradation rate, PINK1 and TOMM20 turnover was measured following lysosome inhibition. PINK1 levels did not show any increase following chloroquine treatment (Fig. 7a and c) while TOMM20 accumulation varied among IPF patients and overall increased significantly (Fig. 7b and c). In AMs significant variation in TOMM20 accumulation following lysosome inhibition was observed while, overall no significant increase was observed among the two groups (Fig. 7d-e).
Additionally, the colocalization of TOMM20 with LC3B II puncta in AMs, by immunofluorescence and confocal microscopy was analyzed, as a more direct assessment of mitophagy. TOMM20 and LC3B II colocalization was not more prominent in IPF than controls following chloroquine treatment. This was confirmed by the lack of significant differences in the colocalization coefficients’ analysis that included Pearson coefficient of colocalization, Costes’s threshold corrected Pearson coefficient of colocalization and M1 and M2 Menders colocalization values(Fig. 7f).
Impaired mitochondria homeostasis and function are established hallmarks of aging and contribute to cell senescence, chronic inflammation and age-dependent decline in stem cell activity  while age associated changes have been linked to IPF pathogenesis [3, 28]. In this study, we used BAL, a tool for both clinical and research purposes [29, 30], and a rich source of AMs. We showed that mitochondria in AMs from IPF patients have prominent morphological defects and impaired transcription. We observed a significant reduction of mitochondria homeostasis regulators PINK1, PARK2 and NRF1. Additionally, a significant increase in mtROS was associated with reduced expression of mitochondria-encoded OXPHOS genes.
PINK1 plays a significant role in mitochondria homeostasis via mitophagy dependent and independent pathways . A decrease in PINK1 levels in type II AECs from IPF patients, IPF lung tissue, aging mice and mice challenged with bleomycin has been associated with the development of fibrosis [9,10,11,12, 32]. PINK1-knockout results in increased mtROS production, mitochondria fragmentation/fission and activation of autophagy, prompting the clearance of oxidized mitochondria [9, 33]. We demonstrated a significant decrease in PINK1 levels, coupled to lower PARK2 transcription, in IPF AMs. Consistently with decreased PINK1 and increased mtROS, morphological characteristics of mitochondria, including prominent association with vesicular structures, swollen, irregular shapes, shorter length, and disorganized cristae, suggestive of mitochondrial fission, were also observed.
The accumulation of oxidized mitochondria could be attributed to a dysregulation of mitochondria quality control mechanism that includes clearance of damaged mitochondria mainly via mitophagy . Autophagosome formation and LC3B activation were previously shown reduced in IPF lungs and following TGF-βtreatment , while fibroblast foci show elevated activation of mTOR, suggesting a downregulation of autophagy . However, a survival mechanism of pro-fibrotic AMs via increased LC3B-dependent mitophagy of PINK1 associated mitochondria has been previously demonstrated . In our study, the endogenous levels of LC3B II were quite prominent in several samples suggesting the presence of the lipid-associated form of LC3B in AMs from both controls and IPF. We did not observe increased levels of LC3B II in IPF as previously shown , while autophagy flux was similar in the two groups. Furthermore, we assessed mitophagy levels by the degree of colocalization of LC3B puncta with mitochondrial protein TOMM20 and no significant differences were observed. Interestingly, we observed a significant increase in lysosomal degradation of TOMM20 in IPF cells. Conceivably, the increased levels of mitochondria degradation could be attributed to increased phagocytosis of cell free mitochondria accumulated in BAL . Alternatively, our results may reflect an imbalance in mitochondria biogenesis and degradation, since a recent study has shown that the biogenesis of mitochondria is also defective in IPF . Additionally, other mitochondria clearance mechanisms have been described such as mitochondria derived vesicles and spheroids that drive mitochondria lysosomal degradation independent of the classic autophagy pathway and may play a role in IPF .
Furthermore, mtROS was significantly increased and associated with increased age and diminished lung function in IPF, suggesting a correlation with disease severity. Persistently increased mtROS, has been associated with mitochondria damage that could blunt macrophage activation and function. Indeed, increased basal levels of mtROS in Chronic Obstructive Pulmonary Disease (COPD) resulted in lower acute mtROS production upon bacterial infection leading to reduced apoptosis and bactericidal activity of AMs . Notably, the persistently increased mtROS in IPF AMs described in this study may explain previous observations showing that innate immune responses such as NLRP3 inflammasome activation and antiviral gene expression are altered in the BAL cells of IPF patients [40, 41].
Immune cell functions including macrophage plasticity are largely shaped by mitochondria activity and cellular metabolism [4, 42]. Activation of Toll-like receptors, such as TLR4, leads to a switch from OXPHOS to glycolysis in immune cells, similar to that occurring in tumors  whereas macrophages activated with IL-4 commit to OXPHOS and fatty acid oxidation (FAO) showing minimal metabolic adaptation when compared to homeostasis . Two populations of AMs exist in the injured lung, resident AMs, and recruited AMs that originate from circulating monocytes. Recruited AMs produce inflammatory cytokines in association with increased glycolytic and arginine metabolism whereas resident AMs proliferate locally and are characterized by increased tricarboxylic acid cycle and amino acid metabolism . Pro-fibrotic mouse AMs showed an increase in glycolysis and glycolytic enzyme transcription, but no dependence on FAO or OXPHOS . Importantly, reversal of macrophage phenotype by 2-deoxy-D-glucose, an inhibitor of glycolysis was shown to reduce fibrosis in the same study . To date, there are no reports of macrophage metabolic profiles from IPF patients. Interestingly, expression of MRC I and V mt-encoded OXPHOS genes was significantly reduced in IPF BAL cells. In addition, the NRF1 transcription factor that regulates mt-respiratory gene expression and mitochondria biogenesis  was significantly downregulated in IPF. Intriguingly, loss of PINK1 also results in decreased MRC complex 1 activity independently of mitophagy . Therefore, our results may suggest a decrease in mitochondria respiratory activity in IPF alveolar macrophages.
According to the previous paradigm of macrophage categorisation into classical/proinflammatory, alternative/resolving and pro-repair/fibrotic population, alternatively activated macrophages were suggested to predominate in IPF as indicated by the increased surface expression of manose receptor CD206 . In our study, the levels of CD45+CD11c+CD163+ alveolar macrophages were significantly higher in IPF BAL. Recent studies propose that a diverse population of macrophages exists in IPF lungs with both pro-inflammatory and pro-fibrotic gene expression [18, 19, 49, 50]. Further studies are required to elucidate the phenotypic and functional characteristics of the AM population in the BAL showing increased mtROS.
Some limitations of our study should be considered, mainly regarding the differences in the number of patients included in each experiment, which is based on the number of cells required for each method, and the diversity in total cell count of cells in each BAL sample. However, it needs to be stressed that patients were randomly selected for each method, and furthermore, mtROS levels were evaluated in most patients included in the study.
The impact of dysregulated mitochondrial complex expression and increased oxidation levels on the antimicrobial action of lung macrophages in health and disease has been recently described [51,52,53]. Furthermore, the association of increased mtROS with aging, a major cofactor in the development of IPF was evident in this series of patients. Our study provides novel findings about mitochondrial gene expression, morphology, and oxidation levels, in alveolar macrophages in IPF that could be indirectly linked with IPF pathogenesis and disease progression. Conclusively, we report that AMs in IPF accumulate oxidized and dysmorphic mitochondria, which correlated with clinical variables of the disease. Impaired mitochondria activity may provide a better understanding of AMs activation in IPF. Importantly, targeted antioxidant therapies may need re-evaluation in IPF with respect to macrophage functions.
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
Beclin-1coiled-coil myosin-like BCL2-interacting protein
Haemoglobin 1 gene
Idiopathic pulmonary fibrosis
Mitochondria respiratory chain
Mitochondria respiratory chain
ATP synthase 6
Mitochondrial encoded NADH dehydrogenase 1
Mitochondrial encoded NADH dehydrogenase 5
Nuclear respiratory factor 1
Selective autophagy receptor p62/sequestosome1 SQSTM1
Parkin E3 ubiquitin ligase
PTEN-induced kinase 1
Reactive oxygen species
Richeldi L, Collard HR, Jones MG. Idiopathic pulmonary fibrosis. Lancet. 2017;389(10082):1941–52.
Wuyts WA, Agostini C, Antoniou KM, Bouros D, Chambers RC, Cottin V, et al. The pathogenesis of pulmonary fibrosis: a moving target. Eur Respir J. 2013;41(5):1207–18.
Zank DC, Bueno M, Mora AL, Rojas M. Idiopathic pulmonary fibrosis: aging, mitochondrial dysfunction, and cellular bioenergetics. Front Med. 2018;5:10.
Angajala A, Lim S, Phillips JB, Kim J-H, Yates C, You Z, et al. Diverse Roles of Mitochondria in Immune Responses: Novel Insights Into Immuno-Metabolism. Front Immunol. 2018;9:1605.
Bratic A, Larsson N-G. The role of mitochondria in aging. J Clin Invest. 2013;123(3):951–7.
Palikaras K, Lionaki E, Tavernarakis N. Balancing mitochondrial biogenesis and mitophagy to maintain energy metabolism homeostasis. Cell Death Differ. 2015;22(9):1399–401.
Palikaras K, Tavernarakis N. Mitochondrial homeostasis: the interplay between mitophagy and mitochondrial biogenesis. Exp Gerontol. 2014;56:182–8.
Mora AL, Bueno M, Rojas M. Mitochondria in the spotlight of aging and idiopathic pulmonary fibrosis. J Clin Invest. 2017;127(2):405–14.
Patel AS, Song JW, Chu SG, Mizumura K, Osorio JC, Shi Y, et al. Epithelial cell mitochondrial dysfunction and PINK1 are induced by transforming growth factor-beta1 in pulmonary fibrosis. PLoS One. 2015;10(3):e0121246.
Bueno M, Lai YC, Romero Y, Brands J, St Croix CM, Kamga C, et al. PINK1 deficiency impairs mitochondrial homeostasis and promotes lung fibrosis. J Clin Invest. 2015;125(2):521–38.
Yu G, Tzouvelekis A, Wang R, Herazo-Maya JD, Ibarra GH, Srivastava A, et al. Thyroid hormone inhibits lung fibrosis in mice by improving epithelial mitochondrial function. Nat Med. 2018;24(1):39–49.
Sosulski ML, Gongora R, Danchuk S, Dong C, Luo F, Sanchez CG. Deregulation of selective autophagy during aging and pulmonary fibrosis: the role of TGFbeta1. Aging Cell. 2015;14(5):774–83.
Kobayashi K, Araya J, Minagawa S, Hara H, Saito N, Kadota T, et al. Involvement of PARK2-Mediated Mitophagy in Idiopathic Pulmonary Fibrosis Pathogenesis. J Immunol. 2016;197(2):504–16.
Daniil Z, Kotsiou OS, Grammatikopoulos A, Peletidou S, Gkika H, Malli F, et al. Detection of mitochondrial transfer RNA (mt-tRNA) gene mutations in patients with idiopathic pulmonary fibrosis and sarcoidosis. Mitochondrion. 2018;43:43–52.
Sato N, Takasaka N, Yoshida M, Tsubouchi K, Minagawa S, Araya J, et al. Metformin attenuates lung fibrosis development via NOX4 suppression. Respir Res. 2016;17(1):107.
Murray LA, Chen Q, Kramer MS, Hesson DP, Argentieri RL, Peng X, et al. TGF-beta driven lung fibrosis is macrophage dependent and blocked by serum amyloid P. Int J Biochem Cell Biol. 2011;43(1):154–62.
Zhou X, Moore BB. Location or origin? What is critical for macrophage propagation of lung fibrosis? Eur Respir J. 2018;51(3).
Aran D, Looney AP, Liu L, Wu E, Fong V, Hsu A, et al. Reference-based analysis of lung single-cell sequencing reveals a transitional profibrotic macrophage. Nat Immunol. 2019;20(2):163–72.
Misharin AV, Morales-Nebreda L, Reyfman PA, Cuda CM, Walter JM, McQuattie-Pimentel AC, et al. Monocyte-derived alveolar macrophages drive lung fibrosis and persist in the lung over the life span. J Exp Med. 2017;214(8):2387–404.
Zhang Y, Choksi S, Chen K, Pobezinskaya Y, Linnoila I, Liu ZG. ROS play a critical role in the differentiation of alternatively activated macrophages and the occurrence of tumor-associated macrophages. Cell Res. 2013;23(7):898–914.
Osborn-Heaford HL, Ryan AJ, Murthy S, Racila A-M, He C, Sieren JC, et al. Mitochondrial Rac1 GTPase import and Electron transfer from cytochrome c are required for pulmonary fibrosis. J Biol Chem. 2012;287(5):3301–12.
Raghu G, Collard HR, Egan JJ, Martinez FJ, Behr J, Brown KK, et al. An official ATS/ERS/JRS/ALAT statement: idiopathic pulmonary fibrosis: evidence-based guidelines for diagnosis and management. Am J Respir Crit Care Med. 2011;183(6):788–824.
Samara KD, Trachalaki A, Tsitoura E, Koutsopoulos AV, Lagoudaki ED, Lasithiotaki I, et al. Upregulation of citrullination pathway: from autoimmune to idiopathic lung fibrosis. Respir Res. 2017;18(1):218.
Tsitoura E, Wells AU, Karagiannis K, Lasithiotaki I, Vasarmidi E, Bibaki E, et al. MiR-185/AKT and miR-29a/collagen 1a pathways are activated in IPF BAL cells. Oncotarget. 2016;7(46):74569–81.
Kaku Y, Imaoka H, Morimatsu Y, Komohara Y, Ohnishi K, Oda H, et al. Overexpression of CD163, CD204 and CD206 on alveolar macrophages in the lungs of patients with severe chronic obstructive pulmonary disease. PLoS One. 2014;9(1):e87400.
Ryu C, Sun H, Gulati M, Herazo-Maya JD, Chen Y, Osafo-Addo A, et al. Extracellular mitochondrial DNA is generated by fibroblasts and predicts death in idiopathic pulmonary fibrosis. Am J Respir Crit Care Med. 2017;196(12):1571–81.
Sun N, Youle RJ, Finkel T. The mitochondrial basis of aging. Mol Cell. 2016;61(5):654–66.
Antoniou KM, Hansell DM, Rubens MB, Marten K, Desai SR, Siafakas NM, et al. Idiopathic pulmonary fibrosis: outcome in relation to smoking status. Am J Respir Crit Care Med. 2008;177(2):190–4.
Raghu G, Remy-Jardin M, Myers JL, Richeldi L, Ryerson CJ, Lederer DJ, et al. Diagnosis of idiopathic pulmonary fibrosis. An official ATS/ERS/JRS/ALAT clinical practice guideline. Am J Respir Crit Care Med. 2018;198(5):e44–68.
Prasse A, Binder H, Schupp JC, Kayser G, Bargagli E, Jaeger B, et al. BAL cell gene expression is indicative of outcome and airway basal cell involvement in idiopathic pulmonary fibrosis. Am J Respir Crit Care Med. 2019;199(5):622–30.
Voigt A, Berlemann LA, Winklhofer KF. The mitochondrial kinase PINK1: functions beyond mitophagy. J Neurochem. 2016;139(Suppl 1):232–9.
Bueno M, Brands J, Voltz L, Fiedler K, Mays B, St Croix C, et al. ATF3 represses PINK1 gene transcription in lung epithelial cells to control mitochondrial homeostasis. Aging Cell. 2018;17(2).
Dagda RK, Cherra SJ 3rd, Kulich SM, Tandon A, Park D, Chu CT. Loss of PINK1 function promotes mitophagy through effects on oxidative stress and mitochondrial fission. J Biol Chem. 2009;284(20):13843–55.
Hara H, Kuwano K, Araya J. Mitochondrial quality control in COPD and IPF. Cells. 2018;7(8).
Patel AS, Lin L, Geyer A, Haspel JA, An CH, Cao J, et al. Autophagy in idiopathic pulmonary fibrosis. PLoS One. 2012;7(7):e41394.
Gui YS, Wang L, Tian X, Li X, Ma A, Zhou W, et al. mTOR Overactivation and compromised autophagy in the pathogenesis of pulmonary fibrosis. PLoS One. 2015;10(9):e0138625.
Larson-Casey JL, Deshane JS, Ryan AJ, Thannickal VJ, Carter AB. Macrophage Akt1 kinase-mediated Mitophagy modulates apoptosis resistance and pulmonary fibrosis. Immunity. 2016;44(3):582–96.
McLelland GL, Soubannier V, Chen CX, McBride HM, Fon EA. Parkin and PINK1 function in a vesicular trafficking pathway regulating mitochondrial quality control. EMBO J. 2014;33(4):282–95.
Bewley MA, Preston JA, Mohasin M, Marriott HM, Budd RC, Swales J, et al. Impaired mitochondrial Microbicidal responses in chronic obstructive pulmonary disease macrophages. Am J Respir Crit Care Med. 2017;196(7):845–55.
Lasithiotaki I, Giannarakis I, Tsitoura E, Samara KD, Margaritopoulos GA, Choulaki C, et al. NLRP3 inflammasome expression in idiopathic pulmonary fibrosis and rheumatoid lung. Eur Respir J. 2016;47(3):910–8.
Lasithiotaki I, Antoniou KM, Vlahava VM, Karagiannis K, Spandidos DA, Siafakas NM, et al. Detection of herpes simplex virus type-1 in patients with fibrotic lung diseases. PLoS One. 2011;6(12):e27800.
O'Neill LA, Pearce EJ. Immunometabolism governs dendritic cell and macrophage function. J Exp Med. 2016;213(1):15–23.
Rodriguez-Prados JC, Traves PG, Cuenca J, Rico D, Aragones J, Martin-Sanz P, et al. Substrate fate in activated macrophages: a comparison between innate, classic, and alternative activation. J Immunol. 2010;185(1):605–14.
Mould KJ, Barthel L, Mohning MP, Thomas SM, McCubbrey AL, Danhorn T, et al. Cell origin dictates programming of resident versus recruited macrophages during acute lung injury. Am J Respir Cell Mol Biol. 2017;57(3):294–306.
Xie N, Cui H, Ge J, Banerjee S, Guo S, Dubey S, et al. Metabolic characterization and RNA profiling reveal glycolytic dependence of profibrotic phenotype of alveolar macrophages in lung fibrosis. Am J Physiol Lung Cell Mol Physiol. 2017;313(5):L834–l44.
Vercauteren K, Pasko RA, Gleyzer N, Marino VM, Scarpulla RC. PGC-1-related coactivator: immediate early expression and characterization of a CREB/NRF-1 binding domain associated with cytochrome c promoter occupancy and respiratory growth. Mol Cell Biol. 2006;26(20):7409–19.
Morais VA, Haddad D, Craessaerts K, De Bock PJ, Swerts J, Vilain S, et al. PINK1 loss-of-function mutations affect mitochondrial complex I activity via NdufA10 ubiquinone uncoupling. Science. 2014;344(6180):203–7.
Pechkovsky DV, Prasse A, Kollert F, Engel KM, Dentler J, Luttmann W, et al. Alternatively activated alveolar macrophages in pulmonary fibrosis-mediator production and intracellular signal transduction. Clin Immunol. 2010;137(1):89–101.
Lee J, Arisi I, Puxeddu E, Mramba LK, Amicosante M, Swaisgood CM, et al. Bronchoalveolar lavage (BAL) cells in idiopathic pulmonary fibrosis express a complex pro-inflammatory, pro-repair, angiogenic activation pattern, likely associated with macrophage iron accumulation. PLoS One. 2018;13(4):e0194803.
Allden SJ, Ogger PP, Ghai P, McErlean P, Hewitt R, Toshner R, et al. The transferrin receptor CD71 delineates functionally distinct airway macrophage subsets during idiopathic pulmonary fibrosis. Am J Respir Crit Care Med. 2019;200(2):209–19.
Plataki M, Cho SJ, Harris RM, Huang HR, Yun HS, Schiffer KT, et al. Mitochondrial dysfunction in aged macrophages and lung during primary Streptococcus pneumoniae infection is improved with Pirfenidone. Sci Rep. 2019;9(1):971.
Liu Y, Lu F, Kang L, Wang Z, Wang Y. Pirfenidone attenuates bleomycin-induced pulmonary fibrosis in mice by regulating Nrf2/Bach1 equilibrium. BMC Pulm Med. 2017;17(1):63.
Belchamber KBR, Singh R, Batista CM, Whyte MK, Dockrell DH, Kilty I, et al. Defective bacterial phagocytosis is associated with dysfunctional mitochondria in COPD macrophages. Eur Respir J. 2019;54.
The authors would like to thank Margarita Stapounzi and the confocal microscopy unit of FORTH-IMBB, Herakleion.
This work was supported by a Hellenic Thoracic Society Research Award (2015).
Ethics approval and consent to participate
The study was approved by the Ethics Committees of the University Hospital of Heraklion (IRB number: 1045 and17030) and all patients provided written informed consent.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Figure S1a. A549 cells treated with 1000 μM H2O2 were stained with 5 μm MitoSOXTMred and examined using confocal microscopy. The characteristic mitochondrial network/web staining in the cytoplasm observed, was significantly enhanced upon H202 treatment. Figure S1b. PMA treated THP1 cells were stained with 5 μm MitoSOXTMred or Propidium Iodide and analysed by flow cytometry. The non-treated PI negative cell population stained positive with MitoSOX while treatment with 1000 μM H2O2 resulted in stronger MitoSOX staining (pale pink histograms). Following H2O2 treatment the PI positive population showed the highest MitoSOX staining (dark red histogram). Figure S1c. Representative example of BAL sample stained with CD45-FITC/CD11c-PC5 or CD45-FITC/CD14-PC. The CD45 + FShighSShigh population in panel (i) comprises of CD11c positive cells (panel ii) and a lower percentage of CD14 positive cells (iii). The CD45 + FShighSShigh population of alveolar macrophages/monocyte derived phagocytes, was used for MitoSOXTMred analyses. Figure S1d. Representative examples of control and IPF BAL MitoSOX analyses.
Additional file 2: Table S1. Primer sequences for qPCR. Table S2. Clinicopathological characteristics of the subjects included in Fig. 1b. Values are expressed as means±SD or medians with range. Table S3. Clinicopathological characteristics of the subjects included in Fig. 1f and g, Fig. 3c-d and Fig. 4a-d. Values are expressed as means±SD or medians with range. Table S4. Clinicopathological characteristics of the subjects included in Fig. 2b-d Values are expressed as means±SD or medians with range. Table S5. Clinicopathological characteristics of the subjects included in Fig. 3b. Values are expressed as means±SD or medians with range.
Additional file 3: Figure S2. Mitochondrial DNA content assessed by mtDNA/gDNA ratio measuring both ND1(A) and ND5(B) relative to HGB1 in BAL cells. No statistical difference between control (n = 10) and IPF (n = 18). Representative western blot (C) and densitometry analysis (D) of mitochondria protein TOMM20 showing no statistical difference between control (n = 5) and IPF (n = 11). TOMM20 mean fluorescence intensity/cell measured from images acquired by confocal microscopy, per patient (E) and group (F). Non significant differences are highlighted by doted lines, all other comparisons showed statistically significant differences according to Kruskal Wallis test with Dunn’s test for multiple comparisons. (G) Representative image of TOMM20 staining of alveolar macrophages from IPF patient (p1). With the exception of (E) all for all other comparisons Mann-Whitney test was used, data are represented as median with interquartile ranges.
Additional file 4: Figure S3. A, B: (i) LC3b I, (ii)II and (iii) LC3bII/I densitometry analyses from western blots performed on 2 independent IPF/control cohorts, C and D: relative mRNA expression of p62 and Beclin1. E. Representative images of LC3b immunostained alveolar macrophages from IPF and control samples treated with CQ or left untreated. All graphs represent median and interquartile ranges, no significant differences observed according to Mann-Whitney test.
About this article
Cite this article
Tsitoura, E., Vasarmidi, E., Bibaki, E. et al. Accumulation of damaged mitochondria in alveolar macrophages with reduced OXPHOS related gene expression in IPF. Respir Res 20, 264 (2019) doi:10.1186/s12931-019-1196-6
- Idiopathic pulmonary fibrosis
- Alveolar macrophages