纳特医学。作者手稿;于 PMC 2018 年 6 月 4 日上市。
以最终编辑形式发布为:
甲状腺激素通过改善上皮线粒体功能抑制小鼠肺纤维化
于国英, 1, * Argyris Tzouvelekis , 1, 2, *王蓉, 1 Jose D. Herazo-Maya , 1 Gabriel H. Ibarra , 1 Anup Srivastava , 1 Joao Pedro Werneck de Castro , 3, 4 Giuseppe DeIuliis , 1 Farida Ahangari , 1 Tony Woolard , 1 Nachelle Aurelien , 1 Rafael Arrojo e Drigo , 5 Ye Gan , 1 Morven Graham , 6 Xinran Liu , 6 Robert J. Homer , 7, 8 Thomas S. Scanlan , 9 Praveen Mannam , 1 Patty J. Lee , 1 Erica L. Herzog , 1 Antonio C. Bianco , 3和Naftali Kaminski 1
Guoying Yu
1Section of Pulmonary, Critical Care and Sleep Medicine, Department of Internal Medicine, Yale School of Medicine, New Haven, CT
Argyris Tzouvelekis
1Section of Pulmonary, Critical Care and Sleep Medicine, Department of Internal Medicine, Yale School of Medicine, New Haven, CT
2Division of Immunology, Biomedical Sciences Research Center “Alexander Fleming”, Athens, Greece
Rong Wang
1Section of Pulmonary, Critical Care and Sleep Medicine, Department of Internal Medicine, Yale School of Medicine, New Haven, CT
Jose D. Herazo-Maya
1Section of Pulmonary, Critical Care and Sleep Medicine, Department of Internal Medicine, Yale School of Medicine, New Haven, CT
Gabriel H. Ibarra
1Section of Pulmonary, Critical Care and Sleep Medicine, Department of Internal Medicine, Yale School of Medicine, New Haven, CT
Anup Srivastava
1Section of Pulmonary, Critical Care and Sleep Medicine, Department of Internal Medicine, Yale School of Medicine, New Haven, CT
Joao Pedro Werneck de Castro
3Division of Endocrinology/Metabolism, Rush University Medical Center, Chicago IL
4Biophysics Institute, Federal University of Rio de Janeiro, RJ, Brazil
Giuseppe DeIuliis
1Section of Pulmonary, Critical Care and Sleep Medicine, Department of Internal Medicine, Yale School of Medicine, New Haven, CT
Farida Ahangari
1Section of Pulmonary, Critical Care and Sleep Medicine, Department of Internal Medicine, Yale School of Medicine, New Haven, CT
Tony Woolard
1Section of Pulmonary, Critical Care and Sleep Medicine, Department of Internal Medicine, Yale School of Medicine, New Haven, CT
Nachelle Aurelien
1Section of Pulmonary, Critical Care and Sleep Medicine, Department of Internal Medicine, Yale School of Medicine, New Haven, CT
Rafael Arrojo e Drigo
5The Salk Institute for Biological Studies, Molecular and Cell biology laboratory, La Jolla, CA
Ye Gan
1Section of Pulmonary, Critical Care and Sleep Medicine, Department of Internal Medicine, Yale School of Medicine, New Haven, CT
Morven Graham
6CCMI Electron Microscopy Core Facility, Yale University School of Medicine, New Haven, CT
Xinran Liu
6CCMI Electron Microscopy Core Facility, Yale University School of Medicine, New Haven, CT
Robert J. Homer
7Department of Pathology, Yale University School of Medicine, New Haven, C
8Pathology and Laboratory Medicine Service, VA CT HealthCare System, West Haven, CT
Thomas S. Scanlan
9Department of Physiology and Pharmacology, Oregon Health and Science University, Portland, Oregon, USA
Praveen Mannam
1Section of Pulmonary, Critical Care and Sleep Medicine, Department of Internal Medicine, Yale School of Medicine, New Haven, CT
Patty J. Lee
1Section of Pulmonary, Critical Care and Sleep Medicine, Department of Internal Medicine, Yale School of Medicine, New Haven, CT
Erica L. Herzog
1Section of Pulmonary, Critical Care and Sleep Medicine, Department of Internal Medicine, Yale School of Medicine, New Haven, CT
Antonio C. Bianco
3Division of Endocrinology/Metabolism, Rush University Medical Center, Chicago IL
Naftali Kaminski
1Section of Pulmonary, Critical Care and Sleep Medicine, Department of Internal Medicine, Yale School of Medicine, New Haven, CT
出版商对本文的最终编辑版本可在Nat Med上找到
相关数据
- 补充材料
- NIHMS915927-supplement-1.pdf (2.2M)GUID: 07C30019-E1BA-4430-827F-C940C5BF6141
抽象的
甲状腺激素(TH)对于应激反应期间维持细胞稳态至关重要,但其在肺纤维化中的作用尚不清楚。在这里,我们发现,与对照个体相比,特发性肺纤维化患者肺部的碘甲状腺原氨酸脱碘酶 2 (DIO2)(一种激活 TH 的酶)的活性和表达较高,并且与疾病严重程度相关。我们还发现Dio2 基因敲除小鼠表现出增强的博来霉素诱导的肺纤维化。在两种小鼠肺纤维化模型(气管内博莱霉素和诱导型 TGF-β1)中,雾化 TH 递送可提高存活率并解决纤维化问题。 Sobetirome 是一种 TH 模拟物,也能减弱博来霉素诱导的肺纤维化。在博来霉素诱导的损伤后给予TH,在体内和体外都促进了肺泡上皮细胞中线粒体的生物合成,改善了线粒体生物能并减弱了线粒体调节的细胞凋亡。 TH 并没有减弱Ppargc1a或Pink1敲除小鼠的纤维化,表明对这些途径的依赖。我们的结论是,TH 的抗纤维化特性与肺泡上皮细胞的保护和线粒体功能的恢复有关,因此可能是肺纤维化的有效治疗方法。
关键词:肺纤维化,线粒体,甲状腺激素,Sobiterome,肺泡上皮细胞
肺纤维化(PF)描述了一种情况,其中正常的肺解剖结构被主动重塑、细胞外基质沉积以及成纤维细胞和肺泡上皮细胞表型的巨大变化所取代。这种情况可以是特发性的,如特发性肺纤维化 (IPF),也可以继发于遗传性疾病、自身免疫性疾病或暴露于环境毒素、化学战、药物、外来抗原或辐射1 。 IPF 是最常见的特发性 PF,在病理学上常常无法与其他形式区分,尤其是在疾病的后期。 IPF 影响着美国约 120,000 人,其发病率和死亡率都在稳步上升。 IPF不进行移植的中位生存期约为3年,使其成为预后最严重的非癌症相关肺部疾病2 。
IPF 是上皮细胞损伤和激活的多个周期的结果,其引发肌纤维母细胞病灶的形成、细胞外基质的积累和异常的伤口修复3 。多种途径,包括持续转化生长因子-β1 (TGF-β1) 激活、异常基质沉积、上皮损伤和凋亡、成纤维细胞激活和肌成纤维细胞转分化、发育途径异常重演、内质网应激、端粒酶突变、氧化损伤、金属蛋白酶激活和信号传导4以及 mRNA 5和 microRNA 6表达的广泛变化与 IPF 有关。最近有人提出,其中许多特征可以用衰老的作用来解释,实际上衰老的许多标志包括基因组不稳定、端粒磨损、表观遗传改变、蛋白质稳态丧失、营养感应失调、线粒体功能障碍、细胞衰老、干细胞耗竭和细胞间通讯改变可被视为纤维化肺的特征7 , 8 。
在本研究的背景下,线粒体功能障碍和 PF 代谢异常的新证据非常引人注目8 。代谢重编程是肌成纤维细胞分化的一个特征9 , 10 ,IPF 肺部的肺泡上皮细胞表现出大量受损的线粒体11并且在 IPF 受试者的血浆和支气管肺泡灌洗液中发现游离线粒体 DNA 水平增加9 。与健康对照相比,IPF 受试者肺部中 PTEN 诱导的推定激酶 1 (PINK1)(一种线粒体自噬调节因子)的含量较低,并且Pink1敲除小鼠表现出线粒体变形,并且比野生型小鼠更容易受到 PF 的影响11 。考虑到 IPF 中 PINK1 下调的原始观察结果源自我们生成的公共数据集( http://www.lung-genomics.org - GEO 登录号 GSE47460) 12 ,我们寻找表明基因表达的其他变化肺生物能学或代谢异常。我们发现 DIO2(一种碘甲状腺原氨酸脱碘酶,可将激素原甲状腺素 (T4) 转化为活性分子 3,5,3'-三碘甲状腺原氨酸 (T3))是 IPF 患者肺部增加最显着的基因之一。
本研究的目的是跟进这一发现,并研究甲状腺激素(TH)信号在实验和人类 PF 发病机制中的作用及其对 IPF 线粒体功能的潜在影响。简而言之,我们发现 IPF 患者肺部 DIO2 的上调很可能代表着促进 T4 向 T3 局部转化的努力,以改善纤维化肺部应激环境中肺泡上皮细胞的代谢状态。在两个实验模型中, Dio2敲除小鼠表现出增强的博来霉素诱导的肺纤维化,而雾化输送 TH 则减弱了纤维化。 TH 对小鼠肺纤维化的有益作用与肺泡 II 型上皮细胞 (AECII) 线粒体功能和形态的正常化有关。 T3介导的博来霉素诱导的线粒体异常的恢复与体内和体外肺上皮细胞中线粒体调节的凋亡的抑制有关,并且依赖于线粒体生物合成和线粒体自噬调节因子的完整表达,包括过氧化物酶体增殖物激活受体γ共激活剂α (Ppargc1a) 和 Pink1。鉴于我们发现 THR 选择性拟甲状腺化合物在博来霉素肺纤维化模型中具有抗纤维化特性,这些作用可能是通过 T3 与甲状腺激素受体 (THR) 的相互作用来介导的。
RESULTS
DIO2 is higher in IPF and its absence augments PF in mice
Reanalysis of the Lung Genomic Research Consortium (LGRC) dataset revealed that DIO2 was among the significantly highest expressed genes in lungs of patients with IPF compared to controls (Fig. 1a and Supplementary Table 1). DIO2 gene expression was negatively correlated with disease severity as reflected by diffusion capacity for carbon monoxide (DLCO % predicted) (Fig. 1b). DIO2 mRNA (Fig. 1c) and protein (Fig. 1d) were significantly higher (8- and 3-fold, respectively) in the lungs of patients with IPF compared to normal histology controls (P < 0.001). DIO2 enzymatic activity, assessed by in vitro T4 to T3 conversion, was also higher (3-fold, P < 0.05) in lung extracts from IPF patients compared to controls (Fig. 1e). DIO2 was mainly expressed in the alveolar epithelium surrounding the fibrotic interstitium in lungs of patients with IPF, and weakly expressed in the normal alveolar epithelium (Fig. 1f). In line with human data, Dio2 expression (Fig. 1g) and activity (Fig. 1h) were also significantly induced in the bleomycin model of lung fibrosis in mice.
To determine the role of DIO2 increases in human lung fibrosis we administered bleomycin to Dio2 knockout (Dio2−/−) mice that are systemically euthyroid but exhibit localized hypothyroidism in Dio2-expressing tissues13,14. Both Dio2−/− mice and wild-type littermates developed pulmonary fibrosis 14 days after intratracheal bleomycin with significant increases in hydroxyproline levels (1.7- and 1.6-fold, respectively) (Fig. 1i) compared to saline-treated animals. However, after bleomycin Dio2−/− mice exhibited significantly higher hydroxyproline levels (1.3-fold) (Fig. 1i) and increased histological evidence of fibrosis in the lung (Fig. 1j) compared to wild-type littermates, but no increases in collagen mRNA (data not shown).
TH blunts lung fibrosis in two murine models
Systemic administration of T4 (100μg/kg) at days 10, 12, 14, 16, after bleomycin administration significantly blunted fibrosis but caused significant increases in serum T3 levels (Supplementary Fig. 1). To avoid the side effects of systemic administration of T4 and maximize therapeutic efficacy we focused on the pulmonary delivery of the physiologically active hormone T3 by aerosol. Aerosolized T3 (40 μg/Kg) therapy, given every other day starting at the established fibrosis stage (days 10–20 after bleomycin), caused significant decreases in hydroxyproline levels (2.25-fold, P < 0.05) (Fig. 2a), and histologic evidence of fibrosis (Fig. 2b) without affecting serum T3 levels (Fig. 2b). The effects were comparable to those observed with oral administration of pirfenidone (100 mg/kg) or nintedanib (60 mg/kg), the two drugs currently approved by the US FDA for the treatment of human IPF15 (Fig. 2a,b).
Aerosolized T3 therapy reversed the significant reduction in static lung compliance induced by bleomycin, reflecting a functional improvement in respiratory mechanics and reduction in fibrosis (Fig. 2d). Similar functional benefits were observed with pirfenidone and nintedanib treatment (data not shown). To observe the effects of aerosolized T3 on survival we used a higher dose of bleomycin (3.0 U/kg). In contrast to pirfenidone and nintedanib, aerosolized T3 caused a substantial survival benefit with significant reduction in mortality from 90% in the untreated group to 20% in the T3 treated arm (Fig. 2e). This survival benefit suggested that T3 potentially limited the extent of acute lung injury caused by high dose bleomycin in addition to its anti-fibrotic effects. To assess the effects of aerosolized T3 in another animal model of experimental lung fibrosis, we used the triple transgenic mice in which a CC10 promoter drives the expression of active TGFβ1 in response to doxycycline treatment16. Aerosolized T3 therapy, given every other day, at days 10–20 after doxycycline treatment resulted in significantly lower pulmonary fibrosis as assessed by lower hydroxyproline levels (1.15-fold, P < 0.05) and Masson Trichrome staining compared to mice treated with doxycycline and vehicle (Fig. 2f,g).
TH reverses bleomycin-induced mitochondrial abnormalities
As previously observed11 AECIIs in bleomycin-treated mice exhibited swollen mitochondria, with severely damaged electron-lucent cristae and disrupted membranes and significantly reduced numbers of normal looking mitochondria, while aerosolized T3 reversed these changes (Fig. 3a,b). Aerosolized T3 also reversed many of the functional abnormalities induced by bleomycin in AECIIs. Specifically, primary AECIIs isolated on day 21 from the lungs of bleomycin-treated animals (AECBs) exhibited significantly lower mitochondrial membrane potential than AECIIs from mice treated with aerosolized T3 after bleomycin (AECBT3s) or AECIIs isolated from uninjured lungs (AECTRLs) (1.4- and 1.5-fold, respectively) (Fig. 3c). AECBs exhibited significantly decreased oxygen consumption rate (OCR) compared to AECTRLs (3-fold) that was restored to baseline levels in AECBT3s (Fig. 3d). In vitro, bleomycin treatment caused significant declines in the mitochondrial bioenergetic profile of primary human small airway epithelial cells (SAECs), primary mouse lung AECIIs and mouse lung epithelial cells (MLE12), but T3, added 4 hours after bleomycin, reversed this decline. SAECs, AECIIs and MLE12 that were treated with T3 after bleomycin exhibited significantly higher mitochondrial membrane potential (MMP, 1.2-, 1.3-, and 1.3-fold, respectively) (Fig 3e–g) and oxygen consumption rate (OCR, 2.1-, 1.5- and 6.9-fold, respectively) (Fig. 3h–j) compared to cells treated with vehicle after bleomycin. A detailed analysis of the mitochondrial bioenergetic profiles in lung epithelial cells, both in vivo (isolated AECIIs) and in vitro (SAECs, AECIIs, MLE12) showed similar differences between T3-treated groups and controls (Supplementary Fig 2).
To assess the effects of T3 treatment on mitogenesis we measured the Cytochrome c oxidase subunit IV (COX-IV) and Succinate Dehydrogenase Complex Flavoprotein subunit A (SDHA) ratio as previously described17. Epithelial cells treated with bleomycin exhibited a significantly lower COX-IV/SDHA ratio compared to vehicle treatment. Bleomycin treated cells that were also treated with T3 exhibited ratios similar to controls. This effect was observed in the three cell types (Fig. 3k–m, 1.4-, 1.2-, 1.4-fold, respectively). Epithelial cells treated with bleomycin exhibited significantly lower protein expression of PPARGC1A, a transcriptional co-activator of mitochondrial metabolism and biogenesis, compared to vehicle treated controls, but in cells treated with T3 following bleomycin PPARGC1A expression was similar to controls (Fig. 3n). T3 treatment after bleomycin was associated with significant higher levels of PINK1. Changes in autophagic flux, as indicated by altered autophagy protein light chain 3 (LC3B) and p62 protein levels were consistently seen in SAECs but not in AECIIs or MLE12 (Fig. 3n).
TH suppresses mitochondria-regulated apoptosis in AECIIs
After showing that TH opposes bleomycin-induced mitochondrial function abnormalities in primary mouse epithelial cells in vivo and in both human and mouse epithelial cells in vitro, we aimed to determine the effects of TH treatment on mitochondria-regulated apoptotic pathways. Immunoblot analyses revealed that exposure of AECIIs to bleomycin led to higher protein expression of the pro-apoptotic mediator BCL-2 associated X protein (BAX) and lower expression of the anti-apoptotic mediator BCL-xL (Fig. 4a). Importantly, bleomycin promoted the intracellular movement of BAX from the cytosol to mitochondria in SAECs, as indicated by co-localization of BAX with mito-tracker in the immunofluorescence analysis (Fig. 4b). BAX upregulation and mitochondrial translocation were associated with enhanced apoptotic activity as revealed by increased number of TUNEL-positive cells (Fig. 4c,d). This observation was verified in vivo by immunohistochemistry analysis showing increased TUNEL staining in the alveolar epithelium in mouse lungs challenged with bleomycin (Fig. 4e,f). T3 treatment reversed the BAX/BCL-xL ratio (Fig. 4a), attenuated BAX translocation to mitochondria (Fig. 4b) and consequently diminished bleomycin-induced apoptotic injury in AECs, both in vitro and in vivo, as indicated by significant reductions in TUNEL-positive cells (Fig. 4c–e).
Anti-fibrotic effects of TH are dependent on PPARGC1A and PINK1
Expression of PPARGC1A was lower following bleomycin in wild-type, Dio2−/− mice as well as in human IPF lungs (Supplementary Fig. 3). Mice treated with aerosolized T3 after bleomycin mice exhibited significantly higher (1.5-fold) expression of PPARGC1A mRNA in their lungs to mice treated with vehicle after bleomycin (Fig. 5a). To determine whether the anti-fibrotic effects of T3 were mediated through PPARGC1A we used Ppargc1a knockout mice (Ppargc1a−/−). After bleomycin treatment Ppargc1a−/− mice demonstrated significantly higher hydroxyproline levels (1.35-fold) compared to wild-type littermates but similar changes in Col1a1, Col3a1 mRNA levels (Fig. 5b–d). Unlike wild-type littermates, aerosolized T3 did not blunt fibrosis in Ppargc1a−/− mice and hydroxyproline content, Col1a1, Col3a1 mRNA levels and Masson Trichrome staining were unchanged (Fig. 5a–e). Consistent with previous reports18, in vitro treatment of alveolar epithelial cells (A549) induced higher expression and a shift in the cellular localization of PPARGC1A from nuclear, to both nuclear and mitochondrial (Fig. 4f,g). Dronedarone, an inhibitor of Thyroid-Hormone Receptors (THR)19 inhibited the T3-mediated induction and cellular redistribution of PPARGC1A (Fig. 5f,g). THR inhibition by dronedarone also prevented the anti-apoptotic effects of T3, as indicated by lower protein levels of BCL-xL and increased TUNEL-positive cells in the immunofluorescence analysis (Fig. 5g,h). Gold immunostaining for THRA1 and THRB in uninjured mouse lungs localized both receptors to nuclei and mitochondria of murine AECIIs (Fig. 5i). To further support our hypothesis that T3 anti-fibrotic properties are mediated through thyroid hormone receptors we investigated the therapeutic effects of sobetirome, a well-characterized thyroid hormone receptor agonist20,21. Oral administration of sobetirome started 10 days after bleomycin was associated with enhanced resolution of pulmonary fibrosis, as assessed by significantly decreased hydroxyproline content (P < 0.05) and Masson Trichrome staining (Fig 5j,k).
PINK1 is a known positive regulator of PPARGC1A-dependent mitochondrial metabolism22 and its absence renders mice susceptible to lung fibrosis due to inability to process dysfunctional mitochondria through mitophagy11. After demonstrating that the anti-fibrotic effects of TH were dependent on PPARGC1A we assessed the role of PINK1. In wild type mice aerosolized T3-treatment partially restored bleomycin-induced PINK1 downregulation (1.8-fold) (Fig. 6a). PINK1-knockout mice (Pink1−/−) exhibited higher collagen deposition measured by hydroxyproline levels (1.28-fold) (Fig. 6b), Col1a1, (1.5-fold) (Fig. 6c,d) mRNA levels and Masson Trichrome staining (Fig. 6e) compared to wild-types, indicating augmented fibrosis as previously reported11. Aerosolized T3 administration did not blunt fibrosis in Pink1−/− mice (Fig. 5b–e).
DISCUSSION
Our results here suggest that mitochondrial dysfunction, increasingly recognized as a hallmark of epithelial cell injury in pulmonary fibrosis11, can be attenuated using localized thyroid hormone delivery, leading to resolution of fibrosis in mouse models. Following the observation that both the expression and activity of DIO2 were substantially increased in lung from patients with IPF and correlated with disease severity, we studied the role of DIO2 and TH signaling in experimental pulmonary fibrosis in mice. Genetic deletion of DIO2 enhanced bleomycin-induced lung fibrosis. TH delivery to the lung blunted established pulmonary fibrosis in two experimental models as did oral therapy with sobiterome, a well characterized thyromimetic drug. TH restored bleomycin-induced mitochondrial derangements in epithelial cells in vivo and in vitro through improved mitochondrial function and homeostasis. These beneficial effects were associated with suppression of mitochondria-regulated death pathways, and dependent on intact PPARGC1A and PINK1 pathways.
TH regulates fundamental biological functions including bioenergetics and post-stress injury repair of every vertebrate tissue. Tissue-specific TH metabolism and availability are tightly regulated by DIO1-3. Among them DIO2 is mainly responsible for TH activation by converting T4 to active T323. Increased DIO2 expression and activity may either reflect TH deficiency or increased metabolic requirements23, whereas decreases reflect excess of TH in the body23. These changes are required to protect vital tissues and organs from unexpected fluctuations in TH levels. In the context of lung disease, DIO2 immunoreactivity and enzymatic activity were rapidly induced by ventilation-induced lung injury and DIO2 knockout mice exhibit increased susceptibility to acute lung injury24. TH therapy significantly attenuated experimental lung injury24 and produced similar therapeutic results in animal models of cardiac25 and renal injury26. Clinically, hypothyroidism has been associated with unfavorable prognosis in critically ill patients27, with myocardial infarction28, diabetic nephropathy29 and most recently IPF30. However, evidence for impaired TH signaling in the human IPF lung or experimental evidence in support of its therapeutic utility in IPF were not available. Our study provides novel observations that together establish the case for aerosolized TH therapy as an antifibrotic therapy. Firstly, we provide evidence that DIO2 expression and activity are increased in the lung of patients with IPF, correlate with disease severity and localize mainly to AECs, the highly metabolically active cell type thought to be at the center of the vicious cycle of injury and repair characteristic of IPF. Secondly, we establish a role for TH in experimental lung fibrosis; TH therapy effectively blunts experimental lung fibrosis in two mouse models in a manner that is replicated by the thyroid hormone mimetic sobiterome. Thirdly, the beneficial effect is observed when TH therapy is initiated during the established fibrosis phase, suggesting enhanced resolution of fibrosis. Fourthly, the antifibrotic effects are obtained with aerosolized T3 without an increase in T3 serum levels, suggesting that aerosolized delivery may be effective without the toxicity observed previously with systemic delivery31. Although some systemic deposition of TH following aerosolized delivery could not be entirely excluded as mice may clean their fur following mist exposure; yet, this appears to be non-significant, as assessed by minimal changes in T3 serum levels. Lastly, aerosolized T3 showed comparable or superior effects on lung fibrosis or survival compared to the two US FDA-approved therapeutic compounds for the treatment of IPF, pirfenidone and nintedanib. Taken together these observations establish a novel role for TH as a potential therapeutic agent in pulmonary fibrosis.
Mechanistically we established that TH exerted its antifibrotic effects in bleomycin induced pulmonary fibrosis by restoring mitochondrial function through induction of mitochondrial biogenesis and mitophagy. T3-mediated beneficial effects on mitochondrial homeostasis attenuated bleomycin-induced mitochondria-regulated cell death in alveolar epithelial cells. Recently, there has been significant interest in mitochondrial function in chronic lung diseases9,11,32. Bueno and colleagues discovered that PINK1, a regulator of clearance of dysfunctional mitochondria, was decreased in subjects with IPF11. They convincingly demonstrated abundant mitochondrial abnormalities in AECIIs derived from lungs from patients with IPF, as well as in mice treated with bleomycin11. Based on their findings we hypothesized that the effects of TH depletion or augmentation on severity of bleomycin induced pulmonary fibrosis were associated with changes in mitochondrial function. Indeed, when we isolated primary mouse AECIIs from lungs treated with bleomycin we discovered that they exhibited significant derangements in mitochondrial function including diminished ATP production, ineffective energy metabolism and mitochondrial membrane depolarization, an index of mitochondrial integrity and proton gradient across the mitochondrial inner membrane, and increased proton leak. Aerosolized TH treatment significantly improved all of the parameters of mitochondrial function, as well as mitochondrial morphology in AECIIs. TH did induce increases in proton leak but those correlated with increases in coupling efficiency and metabolic rates as has previously been observed33. This phenomenon may represent restoration of the basal mitochondrial adaptive mechanism to link higher proton conductance to higher metabolic rates resulting in protection against excessive reactive oxygen species production34.
In the context of critical illness, TH exerts cytoprotective and regenerative properties through effects on target genes leading to maintenance of cellular homeostasis, energy expenditure and modulation of stress responses to injury35. We propose a similar but more detailed model in pulmonary fibrosis; specifically, that repeated epithelial injury and increased metabolic demand lead to mitochondrial damage, increased reactive oxygen species (ROS) generation and epithelial cell apoptosis that eventually lead to downstream fibroblast activation and initiation of the fibrotic vicious cycle. TH opposes these processes by mediating mitochondrial biogenesis and effective mitophagy through induction of major transcription coactivators, PPARGC1A and PINK1. The role of PPARGC1A and PINK1 in the beneficial effect of TH is of particular interest. PPARGC1A controls mitochondrial remodeling and mitochondrial biogenesis36. PINK1 is downstream of PPARCG1A and is thought to have an important role in regulating PPARGC1A-dependent mitochondrial biogenesis through selective degradation of damaged mitochondria by mitophagy22. Both PPARGC1A and PINK1 were decreased in the lungs of animals treated with bleomycin, and TH treatment restored their levels. Aerosolized TH delivery failed to rescue the bleomycin-induced fibrotic phenotype in mice with either PPARGC1A or PINK1 deficiency suggesting that TH-mediated antifibrotic effects required intact PPARGC1A and PINK1. The relevance of this finding to human disease is enhanced because the expression of both molecules is decreased in IPF lungs, however it is unclear whether directly overexpressing either of the molecules will have any antifibrotic effect. In addition to in vivo observations, we demonstrated in vitro that T3 caused THR-dependent induction and redistribution of PPARGC1A in alveolar epithelial cells consistent with previous reports in other cells37. THR inhibition with dronedarone resulted in enhanced apoptotic activity and diminished TH-mediated PPARGC1A induction and translocation, whereas treatment with TH restored mitochondrial bioenergetics and reduced apoptotic activity in lung epithelial cells.
Epithelial cell apoptosis and impaired autophagy are increasingly considered hallmarks of pulmonary fibrosis16,38 and they are believed to reciprocally regulate each other39. Caspase-mediated apoptosis has been proposed to inhibit autophagy through cleavage of several autophagy-related proteins including LC3B40. On the other hand, mitophagy reduces the propensity of cells to undergo apoptosis and functions as an early cytoprotective response favoring adaptation to stress and injury by removing damaged mitochondria. PINK1 acts as a key cytoprotective kinase that activates anti-apoptotic members of the BCL-2 family41, which are known as positive regulators of mitochondrial bioenergetics by enhancing coupling efficiency and ATP production42. In line with these studies, we have demonstrated that alveolar epithelial cells exposed to bleomycin exhibited enhanced BAX mitochondrial translocation and diminished BCL-xL expression, mitochondrial membrane depolarization, and downregulation of PINK1 and PPARGC1A leading to impairment of mitophagy and mitochondrial biogenesis, as assessed by changes in the COX-IV/SDHA ratio. These in vitro observations were also supported by in vivo findings showing impaired mitochondrial homeostasis and enhanced apoptosis in AECIIs derived from mice challenged with bleomycin. Impressively, we observed that T3 mediated a sequence of cytoprotective anti-apoptotic events including: attenuation of mitochondrial membrane depolarization, enhancement of mitochondrial biogenesis, induction of effective mitophagy, and finally suppression of the mitochondria-regulated death pathway. These observations are in line with previous studies demonstrating that T3 protected ovarian cells from chemotherapy-induced apoptosis through downregulation of BAX and upregulation of BCL-xL levels43. Similarly, TH has been demonstrated to induce mitochondrial biogenesis and activity through positive regulation of autophagy in skeletal muscle cells44. Taken together, we propose that by restoring mitochondrial function, TH therapy allows functional recovery of alveolar epithelial cells including normalization of mitochondrial bioenergetics and attenuation of mitochondria-regulated apoptosis, which in turn stops the activation of fibroblasts by ROS and damage-associated molecular patterns (DAMPs), and thus allows resolution of fibrosis (Fig. 6f).
In this study we did not address several mechanistically relevant questions of interest. For instance, we did not establish the cause for induction of DIO2 in lung fibrosis. Demonstrating a link between mitochondrial dysfunction and upregulation of DIO2 would strengthen the evidence that the mechanism we propose is indeed also relevant in humans. We also did not study in detail the interaction between the antifibrotic effect of TH and the TGF-β signaling pathway, although we demonstrated that aerosolized T3 therapy was effective in the inducible TGF-β1 model of lung fibrosis. A recent study demonstrated antifibrotic TH effects in models of liver and skin fibrosis45. The study proposed that T3 acted as a modulator of the TGF-β signaling machinery45. We did not test this mechanism directly, and it could still be relevant in our models. Similarly, while our study focused on the beneficial effects of TH on epithelial cells, we did not directly evaluate its role in other cells. Metabolic aberrations, senescence and impaired autophagy have been reported in fibroblasts and macrophages in IPF38,46–50, and it is possible that improved mitochondrial function may also affect these cells. However, we believe that the impressive changes in mitochondrial morphology and function and the dependence of the antifibrotic effects of T3 on intact PPARGC1A and PINK1 pathways in vivo, and the impressive cytoprotective effects in vitro, support the notion that the antifibrotic effects of TH on mitochondrial homeostasis and function in epithelial cells were key mechanisms in our models, but does not rule out the role of additional mechanisms or cell types.
We believe that the observations in subjects with IPF, as well as the detailed and convincing mechanistic experiments we performed in vitro and in animal models of disease do establish a novel and significantly important observation; that TH augmentation, by localized T3 or systemic sobiterome therapy, enhances resolution of fibrosis in mouse models of pulmonary fibrosis, probably through an epithelial cytoprotective effect. Taken together, our results strongly suggest that strategies utilizing lung-targeted TH therapy should be further studied in IPF.
Online Methods
Microarray assays
For gene expression analysis we used the publically available (gene expression omnibus (GEO) GSE47460 and http://www.lung-genomics.org/) LGRC gene expression dataset, generated by our study group12,51. 123 subjects with IPF and 96 control subjects were included. Supplementary Table 1 summarizes the clinical characteristics of the subjects.
IPF Tissues
Lung tissue samples were obtained through the University of Pittsburgh Health Sciences Tissue Bank and Yale university Pathology Tissue service, which were surgical remnants of biopsies or lungs explanted from subjects with IPF who underwent pulmonary transplantation.
Human subject approval
All studies have been approved by the Institutional Review Board at the University of Pittsburgh and Yale School of Medicine.
RNA extraction
Total RNA was extracted from 30 – 50 mg of frozen lung tissue in 700 uL of Qiazol (Lysis buffer, Qiagen, Valencia, CA) according to the manufacturer’s instructions. The purity of the RNA was verified using NanoDrop at 260 nm and the quality of the RNA was assessed using the Agilent 2100 Bioanalyzer (Agilent, Technologies, Santa Clara, CA).
Quantitative Real-Time PCR
Gene expression was determined by TaqMan® (Life Technologies, Thermo Scientific Inc. Rockford IL, USA) according to manufacturer’s instruction. β-glucuronidase (Gusb) was employed as an internal standard control and the specific primers and probes were all obtained from Life Technologies (Thermo Scientific Inc. Rockford IL, USA). Each reaction was performed in triplicate. Relative gene expression was normalized to a value of 1.0 for the unstimulated control group. Control reactions without cDNA or RNA were done as a negative control. Fold change was calculated taking the mean of the controls as the baseline.
Protein extraction and western blot analysis
Proteins from tissues or cultured cells were extracted, transferred onto PVD-membranes, hybridized overnight with appropriate primary antibodies (PPARGC1A - sc-13067 and ab-54481, DIO2-sc-98716, BCL-xL- ab32124, PINK1-ab23707, LC3B-ab51520, p62/SQSTM1-ab56416, BAX-ab32503 and ab-53154, β-actin-sc-47778) obtained from Santa Cruz Biotechnology, CA (sc), or Abcam, Cambridge, UK (ab) and visualized using the Gel Doc XR+ System (Biorad, Lab. Inc. Life Science, USA), according to manufacturer’s protocol.
DIO2 activity
The DIO2 assay was established previously based on the release of radioiodide from the 125I- labelled substrate (100000 cpm-5′-125I of T4-Sigma, MO) as previously described52. Briefly, the mixture containing 25 ml tissue homogenate in 0·1mM-potassium phosphate buffer (pH 7·0) and 1mM-EDTA, and 650 ml of substrate for a final concentration of 1 nM-T4 (Sigma St Louis, MO), 20mM-cofactor diothiothreitol and 1mM-propylthiouracil (pH 7.0) was incubated at 37°C for 1 h. The reactions were stopped by the addition of 50 ml ice-cold 5% bovine serum albumin followed by 350 ml 10% ice-cold TCA, and mixtures were centrifuged at 4000 g for 20 min. The supernatant was further purified by cation exchange chromatography using 1·6ml Dowex 50 W-X2 (100–200 mesh; Sigma). The iodide was then eluted twice with 1ml 10% glacial acetic acid and counted in a g-counter. Enzymatic activity was expressed in fmol125I-released/h per mg protein. Enzymatic activity was expressed in fmol125I-released/h per mg protein.
Animal studies approval
All animal studies were conducted in accordance with the NIH guidelines for humane treatment of animals and were approved by the Institutional Animal Care and Use Committee (IACUC), Yale University (2014-11592).
Experimental models of lung fibrosis
All animal models were performed twice and data was jointly analyzed. All animals treated were included in the analysis. Interventions were not blinded, but analysis of animal samples was. a) Bleomycin-model:C57Bl/6, 9–12 weeks-old, female mice were purchased (Taconic Biosciences, Hudson, NY). Knockout mice (Dio2−/−), Ppargc1a−/− and Pink1−/− of C57/BL6 background were obtained from Jackson Lab (Bar-Harbor, ME). As previously published11,53 genetic deletion of either PINK1 or PPARGC1A does not produce any endogenous lung phenotype and is not associated with derangements of lung physiology. Ppargc1a−/− and Pink1−/− mice are not embryonically lethal, they produce litters of appropriate Mendelian size, are fertile and healthy and present with no gross morphology abnormalities. Mice were anesthetized by placing them in a chamber having paper towels soaked with 40% isoflurane solution diluted with 1,2-propanediol. Mice were randomly assigned to either intratracheal 1.5 U/kg of bleomycin (Hospira, IL) or equivalent volume (50 μL) of 0.9% normal saline was administered intratracheally as previously described54,55. To test therapeutic efficacy of TH in bleomycin-induced established fibrosis we used intraperitoneal T4 (T2376-Sigma Aldrich, 100μg/kg), or aerosolized T3 (T2877-Sigma Aldrich, 40 μg/kg). T4 was administered systemically at days 10, 12, 14, 16 and mice were sacrificed on day 19, following bleomycin-challenge. Aerosolized T3 was administered every other day at days 10–20 and mice were sacrificed on day 21. Dose regimens were based on previously published protocols24,56. Equivalent volumes of normal saline 0.9% were used as controls in all approaches. For survival analysis we used a double dose of bleomycin (3.0 U/kg) and administered aerosolized T3 (40 μg/kg) or normal saline 0.9% at days at days 10–20 (established fibrosis) following challenge and survival data was collected at day 21, respectively. Mice were randomly assigned to receive either pirfenidone (100 mg/kg) or nintedanib (60 mg/kg) or vehicle (0.9% saline) via oral gavage, as previously described57, at days 10–20 (on a daily basis) following bleomycin administration and mice were sacrificed at day 21. Aerosolized delivery of T3 was performed following a standardized protocol. Briefly: T3 was diluted to a final concentration of 40 μg/kg in 6 ml of PBS and suspension was aerosolized using a conventional aerosol nebulizer (Omron) throughout a chamber that allowed simultaneous exposure of 8 mice for 30 minutes until mist stopped forming in the nebulizer chamber. b) TGF-β1-induced lung fibrosis: Inducible lung targeted TGF-β1-overexpressing triple transgenic mice (CC10-rtTA-tTS-TGF-β1) generated as previously described16 were used. Briefly: A triple transgenic system based on the tetracycline-controlled transcriptional suppressor (tTS) and the reverse tetracycline transactivator (rtTA) was constructed. In this system, the CC10 promoter constitutively drives the expression of rtTA and tTS in a lung-specific fashion. In the absence of dox, tTS binds to and actively suppresses the expression of the tet-O–regulated TGF-β1 transgene. In the presence of dox, tTS is released allowing the activating, dox binding rtTA to bind to the tet-O and activate transgene expression. As expected induced TGF-β1 – overexpression upon addition of doxycycline water on a daily basis (days 0–20) caused airway and parenchymal fibrotic response as assessed by increased lung collagen deposition indicated by hydroxyproline levels and Masson Trichrome staining. Aerosolized T3 (40 μg/kg) or normal saline 0.9% was administered every other day on days 10–20 following addition of doxycycline and mice where sacrificed on day 21. c) Sobiterome therapeutic protocol: Sobiterome was provided by the laboratory of Dr. Scanlan as previously described21,58. C57Bl/6, 9–12 weeks-old, female mice were randomly assigned to be challenged with 1.5 U/kg of bleomycin or equivalent volume of normal saline 0.9% at day 0. Mice were then randomly assigned to treatment with 5 mg/kg of sobiterome diluted in 50 μL of normal saline 0.9% or vehicle (equivalent volume of normal saline 0.9%) administered by oral gavage at days 10, 12, 14, 16 and 18 following bleomycin administration and mice were sacrificed on day 21.
In-vivo respiratory mechanics (Flexi-Vent system-SCIREQ)
To evaluate in-vivo respiratory mechanics as an assessment tool for therapeutic efficacy we applied the Flexi-Vent computer-controlled piston-ventilator system according to manufacturer’s protocol. Briefly: Mice were anesthetized with ketamine and then mice tracheotomized and ventilated using the Flevi-Vent system with an average breathing frequency of 150 breaths/min. During the maximal vital capacity (MVC) perturbation, the lungs were inflated to a standard pressure of +30 cmH2O and then deflated to determine MVC. Maximal pressure–volume-loops (PV-loops) between TLC (+30 cmH2O) and FRC (2.8 cmH20) were finally generated to obtain static compliance (Cst) of the respiratory system. In every animal, each manoeuver was repeated until three acceptable measurements (coefficient of determination ≥0.95) were recorded. The average of three acceptable measurements was calculated.
Histology and immunohistochemistry
Tissue sections (4 μm) were stained with Masson Trichrome (collagen/connective tissue), two slices per animal, and two animals per group, as previously described. Immune staining for DIO2 (ab77481, Abcam, Cambridge, UK) or PPARGC1A (ab54481, Abcam, Cambridge, UK) was performed after paraffin removal, hydration, and blocking, following the recommendation of the manufacturer and described by us59. Sections were incubated overnight at 4°C with the primary antibody (diluted 1:100 in PBS) and during 1 hour at room temperature with the secondary antibodies (Sigma, USA). The sections were counterstained with hematoxylin. The primary antibody was replaced by non-immune serum for negative controls.
Hydroxyproline assay
Lung hydroxyproline was analyzed with hydroxyproline colorimetric assay kit from Biovision (Milpitas, CA) following manufacturer’s instruction, as previously described60. Data are expressed as μg of hydroxyproline/right lung.
Measurement of T3 serum levels
Serum from mouse peripheral blood was collected by BD vacutainers, and serum total total triiodothyronine T3 levels were measured T3 ELISA kits (ABIN2685558, Antibodies-Online, Atlanta, GA), following the manufacturer’s instructions.
Isolation of primary alveolar epithelial cells (AECIIs)
Primary AECIIs were isolated from mice treated with 0.9% saline, or bleomycin (1.5U/kg), or bleomycin + aerosolized T3 (40μg/kg), using a modified protocol of a previously published method61. AECIIs were identified as population with CD45/CD32/CD16/CD11b-ve,SSchi cells. Using this protocol, a purity exceeding 80% and viability of 90% was obtained. Isolated cells were used for functional assays.
Cell culture experiments, immunofluorescence staining and assays of apoptosis
Human small airway epithelial cells (SAECs) were purchased from Lonza, Allendale, NJ (LOT: 000470903). Primary mouse AECIIs were isolated from the lungs of unchallenged mice. Mouse lung epithelial cell line (MLE12) was purchased from ATCC (CRL-2110). Cells were seeded into 96-well plates (40.000/well), exposed to bleomycin (15 mU/ml) or PBS for 4 hours and then treated with T3 (15 ng/ml) or vehicle control for 8 hours. Immunofluorescence staining was performed using Mito-Tracker (Thermo Fisher Scientific, Waltham, MA), a cationic dye that stains active mitochondria and PPARGC1A (ab54481, Abcam, Cambridge, UK), according to manufacturer’s instructions. Detection of apoptotic cells was performed with TUNEL-assay using the in situ Cell Death Detection Kit, Fluorescein (Roche, Indianapolis, IA, USA, Catalog No 11684795910). DAPI staining was used to determine the number of nuclei and to assess gross cell morphology. The human lung adenocarcinoma cell line A549 (ATCC® CCL-185™, ATCC, Manassas, VA), negatively tested for mycoplasma, was cultured in DMEM supplemented with 10% fetal bovine serum. Cells were pre-incubated with either dronedarone 10 μM diluted in 0.04% DMSO or vehicle control (0.04% DMSO) for 24 hours and then treated with T3 (15 ng/ml) or vehicle (normal saline 0.9%) for 24 hrs.
Extracellular Flux Technology
The oxygen consumption rates (OCR) of primary mouse AECIIs, SAECs, and MLE-12 cells were measured by using a Seahorse XF96 Extracellular Flux Analyzer (Seahorse Bioscience, Billerica, MA, USA), as previously described62. All assays were performed using a seeding density of 60,000 cells/well in 200μl of DMEM in a XF96 cell culture microplate (Seahorse Bioscience). After the cells were switched to unbuffered DMEM supplemented with 2 mM sodium pyruvate and 20mM carnosine 1h prior to the beginning of the assay and maintained at 37 °C. OCR was measured after sequentially adding to each well 25μl of oligomycin (an ATP-synthase inhibitor), FCCP (a protonophore) and rotenone (inhibitors of complex I and III), to reach working concentrations of 1μg/ml, 1μM and 0.5μM respectively. OCR is reported in picomoles/minute/60,000 cells.
Determination of ATP production
Total ATP was determined using an ATP Fluorometric Assay Kit (BioVision, San Francisco, CA) as previously described62. Briefly, harvested cultured cells were lysed with a lysis buffer followed by centrifugation at 10.000g for 2 min, at 4°C. The level of ATP was determined by mixing 20 μl of the supernatant with 100 μl of luciferase reagent which catalyzed the light production from ATP and luciferin. Reactions were performed in 6-well plates at a seeding density of 106 cells/well. Luminance was measured by a monochromator microplate reader. Total ATP levels were expressed as nmol/mg protein.
Mitochondrial membrane potential (MMP)
The changes in relative MMP (ΔΨm) were measured using the JC-10 (5,5′,6,6′-tetrachloro-1,1′,3,3′-tetraethyl-benzamidazolocarbocyanine iodide) molecular probe, as previously described62.
Transmission Electron microscopy (TEM)
For TEM, mice were perfused with 4% Paraformaldehyde (PFA) and lung samples were dissected out in 2% PFA. These were fixed in 2.5% gluteraldehyde/2% PFA in 0.1M sodium cacodylate buffer pH7.4 for 1 hour, followed by rinsing in PBS buffer then post fixed in 1% osmium tetroxide. They were en-bloc stained in 2% aqueous uranyl acetate then rinsed and dehydrated in an ethanol series followed by resin infiltration Embed 812 (Electron Microscopy Sciences) and baked overnight at 60 C. Hardened blocks were cut using a Leica UltraCut UC7. 60nm sections were collected on formvar/carbon coated grids and contrast stained using 2% uranyl acetate and lead citrate. Grids were viewed FEI TencaiBiotwin TEM at 80Kv. Images were taken using Morada CCD and iTEM (Olympus) software. At least 10 cells from low- and high-magnification images (×11,500, ×26,000) were used to count the number of mitochondria per AECII (identified by the presence of lamellar bodies).
Mitobiogenesis assay
Human SAECs, primary (AECIIs) isolated from the mouse lungs and MLE-12 were seeded in 96-well microplates containing 100 μL of culture medium. Cells were first exposed to bleomycin (10mU/ml) diluted in PBS or vehicle control (PBS) for 4 hours and then were treated with T3 (15ng/ml) or vehicle control for 8 hours. The levels of two mitochondrial proteins were measured simultaneously in each well by using a colorimetric ELISA-kit, according to manufacturer’s protocol (MitoBiogenesis™ In-Cell ELISA Kit- Abcam-ab110217)17. The two proteins are each subunits of a different oxidative phosphorylation enzyme complex, one protein being subunit I of Complex IV (Cytochrome c oxidase subunit-COX-IV), which is mitochondrial-DNA-encoded, and Succinate Dehydrogenase Complex Flavoprotein subunit A-SDHA) which is a 70kDa subunit of Complex II and nuclear-DNA-encoded.
Immunogold staining protocol
The grids containing tissue samples were placed section side down on drops of 1% hydrogen peroxide for 5 minutes, rinsed and blocked for nonspecific binding on 3% bovine serum albumin in Tris buffered saline containing 1% Triton-X for 30 minutes. Grids were incubated with a primary antibody anti-rabbit thyroid hormone receptor alpha-1 and beta (PA1-211A and PA1-213A, respectively Thermo Fisher, Pierce) 1:500 overnight, rinsed in TBS then incubated with the secondary antibody 10nm protein A gold (University of Utrecht, Cell Microscopy Core, Netherlands) for 30 minutes. The grids were well rinsed in PBS, fixed using 1% gluteraldehyde for 5mins, rinsed again, dried and heavy metal contrasted using 2% aqueous uranyl acetate and lead citrate.
Statistical analysis
Data were statistically analyzed using MedCalc version 14. D’Agostino-Pearson test was used to test normal distribution. The results were analyzed by Mann-Whitney U test for comparisons between two groups when sample data were not normally distributed, by unpaired Student’s t-test for comparisons between two groups with normal distribution or by one-way ANOVA with Student-Newman-Keuls post-hoc test for pairwise comparisons, between three or more groups. Efficacy experiments were designed for 10 animals in control and treated groups, to allow for 82% power to detect a difference of 20% between the two groups at the 0.05 statistical significance level, but the actual size of groups differed because of mortality or availability of knockout mice. Data are presented as mean ± SEM, and were considered statistically significant at P < 0.05. Statistical analysis of microarray experiments was performed using BRB array tools as previously described by us63. Multiple comparisons were addressed using the FDR method.
Acknowledgments
The work was in part supported by US National Institute of Health (NIH) grants R01HL095397, R01HL127349 (N.K.), R01HL1109233, R01HL125250 (E.L.H.), RO1DK65055 (A.C.B.), RO1DK52798 (T.S.S.), Flight Attendant Medical Research Institute (P.J.L.), American Lung Association Award RT 350419, Marie Sklodowska/Curie, ERS/EU-RESPIRE 2-8860-2015 grant (A.T.). American Lung Association Award RG-415350 (A.S.), Pulmonary Fibrosis Foundation (PPF) Albert Rose Established Investigator Award 415245 (G.Y.)
A.T., G.Y., and N.K., are inventors on a pending patent on use of thyroid hormone as an antifibrotic agent entitled: “Novel Methods of Treating or Preventing Fibrotic Lung Diseases” - OCR 6368 - 047162-7029P1 (00219). N.K., consulted Biogen Idec, Boehringer Ingelheim, Numedii, MMI, and Pliant and has an ongoing collaboration with MiRagen, all outside the submitted work. E.L.H., consulted for Boehringer Ingelheim and held grant funding from Sanofi and Promedior, all outside the submitted work.
Footnotes
Authors contributions
G.Y., and A.T., performed most of the experiments with the assistance of R.W., J.D.H.-M., J.P.W.d.C., R.A.e.D., F.A., G.D., NA., T.W., J.D., Y.G., R.D., G.H.I., and A.C.B.,. A.S., and P.M., performed the mitochondrial function experiments. M.G., and X.L., performed the electron microscopy experiments and analysis. G.Y., A.T., and N.K performed the histologic analyses with R.J.H., T.S.S., G.H.I., assisted on the GC-1 animal experiments. G.Y., A.T., and N.K., conceived the project, designed the experiments and analyzed and interpreted the results. A.T., G.Y., and N.K., drafted and revised the manuscript. E.L.H., P.J.L., and A.C.B., interpreted results and revised the manuscript.Competing financial interests
All other authors declare no additional competing interests.
Data Availability
All high throughput data mentioned in the paper are publicly available (gene expression omnibus (GEO) GSE47460 and http://www.lung-genomics.org/). All other data that support the findings of this study are available from the corresponding author upon request.