- Open Access
There is detectable variation in the lipidomic profile between stable and progressive patients with idiopathic pulmonary fibrosis (IPF)
Respiratory Research volume 22, Article number: 105 (2021)
Idiopathic pulmonary fibrosis (IPF) is a chronic interstitial lung disease characterized by fibrosis and progressive loss of lung function. The pathophysiological pathways involved in IPF are not well understood. Abnormal lipid metabolism has been described in various other chronic lung diseases including asthma and chronic obstructive pulmonary disease (COPD). However, its potential role in IPF pathogenesis remains unclear.
In this study, we used ultra-performance liquid chromatography-quadrupole time-of-flight mass spectrometry (UPLC-QTOF-MS) to characterize lipid changes in plasma derived from IPF patients with stable and progressive disease. We further applied a data-independent acquisition (DIA) technique called SONAR, to improve the specificity of lipid identification.
Statistical modelling showed variable discrimination between the stable and progressive subjects, revealing differences in the detection of triglycerides (TG) and phosphatidylcholines (PC) between progressors and stable IPF groups, which was further confirmed by mass spectrometry imaging (MSI) in IPF tissue.
This is the first study to characterise lipid metabolism between stable and progressive IPF, with results suggesting disparities in the circulating lipidome with disease progression.
Idiopathic pulmonary fibrosis (IPF) is a chronic, progressive lung disease characterized by alveolar epithelial cell activation and damage, resulting in proliferation of activated fibroblasts and extracellular matrix deposition leading to irreversible destruction of gas exchange units and lung remodelling . Previous studies suggested that IPF is related to abnormalities in a number of biological processes including glycolysis, fatty acid oxidation and vascular remodelling [1,2,3]. Lipidomics and metabolomics are a developing field of systems biology research that studies lipids as key intermediates of cellular mechanisms and their roles have been explored in other respiratory diseases [4, 5].
Ultra-performance liquid chromatography-quadrupole time-of-flight mass spectrometry (UPLC-QTOF-MS) is a powerful tool for qualitative characterization of chemical components, providing acquisition of MS spectra with relatively high resolution, sensitivity and mass accuracy, which has become an advanced tool for the qualitative and quantitative analysis of multiple components. Lipidomic studies of complex biological matrices (plasma and tissue extracts) are routinely performed using UPLC-QTOF-MS, resulting in excellent metabolite detection performance at high speeds and sensitivities . A key challenge in the analysis of complex lipid extracts from crude plasma samples is the co-elution of isomeric lipid species as a result of small variations in their fatty acid arrangements and hydrocarbon backbones . To help overcome this, we have applied SONAR, a rapidly scanning method in tandem with UPLC-QTOF-MS acquisition to provide additional information and quantification of metabolites in complex samples [8, 9].
There have been several studies examining circulatory molecules that characterise progressive IPF, focusing mostly on proteomics and genomics [10,11,12]. Two studies have characterised lipidome differences between IPF and healthy controls, reporting several lipids that have the ability to differentiate IPF from controls [13, 14]. We therefore hypothesize that there are differences in circulating lipids between stable and progressive IPF. To address this, we carried out UPLC-QTOF-MS implementing SONAR, to assess plasma samples from stable and progressive IPF patients.
The study cohort consisted of a total of 58 plasma samples (30 stable and 28 progressors) from the Australian IPF Registry with a clinical diagnosis of IPF (Table 1). Baseline forced vital capacity (FVC) and diffusing capacity for carbon monoxide (DLco) were assessed ± 6 months from the time of blood collection, and the longitudinal FVC and DLco trajectories were determined ± 6–12 months from the baseline lung function using a linear regression model. A decline in FVC ≥ 10% and/or DLco ≥ 15% within 6–12 months of baseline was used to define progressive IPF. All work was approved by the Sydney Local Health District Human Research Ethics Committee (Reference number: HREC 11/RPAH 439), Royal Perth Hospital Human Research Ethics Committee (Reference number: REG 15-204) and the Murdoch University Human Research Ethics Committee (Approval number: 2017/254).
Plasma lipid extraction was carried out as previously described , and each sample was transferred into two different plates corresponding to positive and negative modes of acquisition. This allows ions to be differentiated by its charge (protinate for positively charged and deprotinated for negatively charged molecules). This is useful where the stability of the molecule being characterised is affected by its charge and which mode it runs through, which may also provide complementary structural information through different fragmentation processes . Pooled samples were generated and used as the study reference or for quality control (QC) checks.
Ultra-performance liquid chromatography (UPLC) is an analytical technique used for the separation, identification and further quantitation of compounds in a mixture using pressurized organic and/or polar solvent systems. This was used for the separation, identification and further quantitation of compounds using an ACQUITY I-class system (Waters Corporation, USA) equipped with a Waters CSH C18 column (2.1 × 100 mm, 1.8 µm). Chromatographic separation was achieved with a gradient of 40 to 99% mobile phase over 18 min. Solvent flow rate and column temperature was maintained at 0.4 mL/min and 55 °C, respectively. The lock-mass compound, leucine enkaphlin (200 pg/μL) was prepared in acetonitrile/H2O (50:50, v/v), and was delivered at 10 µL/min to the reference sprayer source of the mass spectrometer.
Mass spectrometry acquisition
MS-based lipid analysis was performed using a Xevo G2-XS QTOF mass spectrometer (Waters Corporation, UK) in positive electrospray ionisation (ESI) mode. Note, positive ionization investigates positive ions in low pH, and negative ionization investigates negative ions at high pH. The source temperature and capillary voltage was set to 120 °C and 2.0 kV, respectively. The time-of-flight (TOF) mass analyser of the mass spectrometer was calibrated using a mass to size ratio (m/z) 50 to 1200.
Peak picked features were statistically analysed using EZinfo (MKS Data Analytics Solutions, Sweden) and the significant features of interest were imported back into Progenesis QI. Lipid identification was achieved by matching experimental fragments against the theoretical fragmentation product ion spectra from LIPID MAPS structure database (Lipidomics Gateway, UK). Additional manual annotation was also conducted using the Lipid Reporter toolkit and reported according to the shorthand nomenclature defined by Liebisch et al. .
The coefficient of variation (CV) was calculated across quality control (QC) samples for each feature and those with a CV > 30% were removed. The high CVs at these levels demonstrate the expected loss of precision when quantifying samples at the extreme ends of the assay's range. Two statistical modelling tools were used to discriminate diseased experimental groups, including the principal component analysis (PCA) and orthogonal projection to latent structures-discriminant analysis (OPLS-DA), based on their contribution to the variation and correlation between the two groups (stable versus progressors). These modeling tools provide insights into separations between experimental groups based on the collected high-dimensional spectral measurements. PCA was initially performed to deconstruct the stable and progressive patient dataset; however, the unsupervised multivariate approach did not reveal significant differences between both groups (Fig. 2a). The heterogeneity between these two diseased groups with similar underlying mechanisms were more likely to cluster together in a principal component space. Supervised multivariate data analysis OPLS-DA was then used to generate a regression model to disentangle group-predictive and group-unrelated variation in the measured data. The OPLS-DA model clearly distinguished the stable and progressor groups from which an S-plot was then generated (Fig. 2b, c). The S-Plot is a statistical tool for visualizing both the covariance and correlation between the endogenous features and the modelled group designation (stable versus progressor). Consolidation of the S-Plots from the OPLS-DA models were useful for the identification of biochemically atypical features with statistical significance between the groups, based on their contributions to the model and their reliabilities. Furthermore, variable importance in projection (VIP) scores of each feature were calculated from the OPLS-DA model, summarising the contributions each lipid makes to the model. Features with VIP scores ≥ 1 were considered as significant .
In order to validate whether changes in the circulation is present in the lung, mass spectrometry imaging (MSI) was carried out to spatially resolve the distribution of the biomolecules resolved by SONAR. The lipids were resolved in a total of 20 fresh frozen tissue sections (20 µm; 10 from IPF patients and 10 healthy controls) using matrix-assisted laser desorption ionisation (MALDI)-QTOF MSI and data was interrogated by High Definition Imaging (HDI) software (Waters Corporation, U.K.) to generate ion intensity maps. The heat maps of a specific ion generated corresponded to the relative abundance of ions present over the entire imaged surface. All data were acquired in positive mode operating over a mass range of m/z 50 to 1200 and were performed using the Water Synapt G2S mass spectrometer equipped with an orthogonal MALDI ion source and an Nd:YAG laser (Waters Corporation, Manchester, U.K.). Putative identification of lipids was achieved by accurate mass measurement and matched against lipid databases including LIPID MAPS and LipidBlast.
Demographic characteristics of this Australian cohort (Table 1) demonstrated for the stable cohort predominantly males (n = 19, 63%), mean age (68 ± 8 year), FVC 82 ± 20% predicted and DLco 53 ± 18% predicted; and for the progressive cohort predominantly males (n = 17, 61%), mean age (70 ± 9 year), FVC 73 ± 15% predicted and DLco 37 ± 14% predicted.
With SONAR, approximately 5000 features were resolved in each analytical run. The total ion chromatograms (TIC) also displayed a typical lipid spectral pattern consistent with previous studies . All lipid classes were resolved by SONAR (Fig. 1).
For statistical analysis, PCA was initially performed to decompress the stable and progressive patient dataset; however, unsupervised multivariate analysis did not reveal significant differences between the groups (Fig. 2a). Supervised multivariate data analysis OPLS-DA was then used to generate a regression model to disentangle group-predictive and group-unrelated variation in the measured data . The OPLS-DA model clearly distinguished the stable and progressor groups from which an S-plot was then generated (Fig. 2b, c) .
From the features exported using EZinfo, eight lipids from positive SONAR mode were putatively identified. The identified lipids included six glycerolipids (triglycerides (TG)) and two glycerophoslipids (phosphatidylcholine (PC)) (Table 2). Using SONAR, three TGs (54:5, 54:6, 53:7) and two PCs (40:6, 36:3) were resolved higher in the progressor subgroup (between 1.9–3.3-fold change). All other lipids observed relatively unchanged levels between the two subgroups (between 0.9–1.11-fold differences).
Exploration of these lipids in IPF tissue revealed an abundance of PC and TG, corresponding to circulatory findings. More specifically, the adduct species [M+K]+ showed improved abundance, while [M+H]+ and [M+NH4]+ showed poor ion spectra (Fig. 3). In this study, these enhanced signal intensities was advantageous as it allowed for the improved discrimination of lipid analytes from the background ions. This is backed up by previous research indicating that [M+K]+ ions were the only adduct type that were well-resolved for the two lipids [22, 23].
Lipids play an important role in lung pathology and physiology. Although the composition and involvement of the lipidome in various diseases is still poorly understood, abnormal lipid metabolism has been reported in few lung diseases including asthma and COPD [4, 5]. Lipids comprise diverse classes of molecules which are critically involved in cellular structure, signalling and energy storage . However, their potential role in IPF pathogenesis remains unclear. This is the first study to evaluate the differences in the lipidome of stable versus progressive IPF patients.
Another novel aspect of this study was the tandem use of SONAR in assessing differences in plasma samples between stable and progressive IPF. It is well documented that conventional quadrupole time-of-flight (QTOF)-MS methods are successful at resolving the lipid classes chromatographically, however, compound identification is still challenging due to the inaccurate assignment of lipid precursors to their corresponding product ions. In contrast, the precursor and product ions generated by SONAR contributed to the specificity of the method, increasing the probability of successful lipid library matching [8, 25]. Furthermore, the MALDI-MSI method allowed for the visualisation of these lipids within IPF tissue.
With SONAR, a number of TG and PC were identified and found to be expressed at higher levels in the IPF progressor group. Specifically, the levels of these TGs (53.7, 54:5 and 54:6) and PCs (36:3 and 40:6) appeared to be elevated in samples of progressors compared to stable patients. This trend is consistent with the findings of Yan et al.  and Kulkarni et al.  and although their observations were based on comparisons between IPF and healthy controls, higher levels in progressive IPF could relate to more active disease process. Specifically, using a bleomycin (BLM) mouse model of pulmonary fibrosis, Kulkarni et al.  observed unchanged levels of specific TGs (48:1, 52:1 and 54:2) consistent with our results, but reported a marked increase of approximately twofold in the relative content of PCs. This suggests that these species may play a role in the pathogenesis of BLM-induced pulmonary fibrosis.
TGs along with diglycerides (DGs), are the most abundant lipids found in circulating plasma. Serum total TG has been reported as a biomarker of fatty acid metabolic disturbance . TGs are stored in lipid droplet structures, formed through budding of the endoplasmic reticulum (ER), and are reported to induce the expression of endogenous ER stress markers (e.g. p-JNK, GRP 78) [27, 28]. Interestingly, ER stress is evident in alveolar epithelium of both human and mice with pulmonary fibrosis . However, neither the mechanisms causing this stress, nor its contribution to fibrosis is well understood. Further investigation of the exact roles of TGs in IPF patients would be beneficial.
Unlike TGs, the functional role of PCs has been described in lung surfactant phospholipid metabolism. Pulmonary surfactant forms the lining of the epithelial air-contact surface of the lungs and are essential for the prevention of alveolar collapse during expiration. PC is the major phospholipid comprising around 80% of surfactant lipids and alterations in their composition can cause reduced elasticity, leading to an overall decrease in lung compliance . Changes in the component phospholipids, including PCs have been described in bronchoalveolar fluids of animal models of rapidly developing pulmonary fibrosis . The up-regulated PC in this study might be partially associated with greater epithelial injury in the progressive group relative to stable.
To understand the global dynamics of metabolite differences in IPF, we performed a metabolic network analysis. We found four networks enriched in the progressive IPF samples (Linoleic acid metabolism, alpha-Linolenic acid metabolism, Arachidonic acid metabolism, Glycerophospholipid metabolism), which have been previously described in fibrotic lung [32, 33]. Based on these findings, it has been proposed that there is a fine balance between lipid metabolism and wound healing mechanisms leading to fibrosis onset and development in IPF .
Although a limitation in this study is its small cohort size and the lack of a healthy cohort, abnormalities in plasma lipids in IPF versus healthy controls have been previously described , 14] forming the basis of our lipidomic exploration. The biological significance pertaining to the fold change differences will need to be further explored in larger cohorts and with functional studies.
In conclusion, this project successfully profiled plasma samples obtained from two groups of IPF patients using the SONAR acquisition approach enhancing the specificity of unbiased lipid profiling derived from UPLC-QTOF-MS. In particular, we have identified changes related to disease progression in lipid signatures such as TGs and PCs in IPF plasma samples as whole, which have previously been associated with abnormalities in lipid metabolism through mitochondrial-beta oxidation pathways.
Availibility of data and materials
All data generated or analysed during this study are included in this published article (and its additional information files).
Chronic obstructive pulmonary disease
Coefficient of variation
Diffusing capacity for carbon monoxide
Forced vital capacity
High definition imaging
Idiopathic pulmonary fibrosis
Mass to size ratio
Matrix-assisted laser desorption ionisation
Mass spectrometry imaging
Orthogonal projections to latent structures discriminant analysis
Principal component analysis
Total ion chromatograms
Ultra-performance liquid chromatography-quadrupole time-of-flight mass spectrometry
Variable importance in projection
Kang YP, Lee SB, Lee J-M, et al. Metabolic profiling regarding pathogenesis of idiopathic pulmonary fibrosis. J Proteome Res. 2016;15(5):1717–24.
Hamanaka RB, Nigdelioglu R, Meliton AY, et al. Inhibition of phosphoglycerate dehydrogenase attenuates bleomycin-induced pulmonary fibrosis. Am J Respir Cell Mol Biol. 2018;58(5):585–93.
Zhao YD, Yin L, Archer S, et al. Metabolic heterogeneity of idiopathic pulmonary fibrosis: a metabolomic study. BMJ Open Respir Res. 2017;4(1):e000183.
Berry KAZ, Murphy RC, Kosmider B, et al. Lipidomic characterization and localization of phospholipids in the human lung. J Lipid Res. 2017;58(5):926–33.
Telenga ED, Hoffmann RF, t’Kindt R, et al. Untargeted lipidomic analysis in chronic obstructive pulmonary disease Uncovering sphingolipids. Am J Respir Crit Care Med. 2014;190(2):155–64.
Ling J, Yu Y, Long J, et al. Tentative identification of 20 (S)-protopanaxadiol metabolites in human plasma and urine using ultra-performance liquid chromatography coupled with triple quadrupole time-of-flight mass spectrometry. J Ginseng Res. 2019;43(4):539–49.
Kyle JE, Zhang X, Weitz KK, et al. Uncovering biologically significant lipid isomers with liquid chromatography, ion mobility spectrometry and mass spectrometry. Analyst. 2016;141(5):1649–59.
Gethings LA, Richardson K, Wildgoose J, et al. Lipid profiling of complex biological mixtures by liquid chromatography/mass spectrometry using a novel scanning quadrupole data-independent acquisition strategy. Rapid Commun Mass Spectrom. 2017;31(19):1599–606.
Moseley MA, Hughes CJ, Juvvadi PR, et al. Scanning quadrupole data-independent acquisition, part A: qualitative and quantitative characterization. J Proteome Res. 2018;17(2):770–9.
Nance T, Smith KS, Anaya V, et al. Transcriptome analysis reveals differential splicing events in IPF lung tissue. PLoS ONE. 2014;9(3):e92111-e.
Herazo-Maya JD, Sun J, Molyneaux PL, et al. Validation of a 52-gene risk profile for outcome prediction in patients with idiopathic pulmonary fibrosis: an international, multicentre, cohort study. Lancet Respir Med. 2017;5(11):857–68.
Boon K, Bailey NW, Yang J, et al. Molecular phenotypes distinguish patients with relatively stable from progressive idiopathic pulmonary fibrosis (IPF). PLoS ONE. 2009;4(4):e5134-e.
Yan F, Wen Z, Wang R, et al. Identification of the lipid biomarkers from plasma in idiopathic pulmonary fibrosis by Lipidomics. BMC Pulm Med. 2017;17(1):174.
Kulkarni YM, Dutta S, Iyer AKV, et al. A lipidomics approach to identifying key lipid species involved in VEGF-inhibitor mediated attenuation of bleomycin-induced pulmonary fibrosis. Proteomics Clin Appl. 2018;12(3):e1700086.
Wong MWK, Braidy N, Pickford R, et al. Comparison of single phase and biphasic extraction protocols for lipidomic studies using human plasma. Front Neurol. 2019;10:879.
Barwick V, Langley J, Mallet T, et al. Best practice guide for generating mass spectra. Teddington: LGC Limited; 2006. ISBN 978-0-948926-24-2
Liebisch G, Vizcaíno JA, Köfeler H, et al. Shorthand notation for lipid structures derived from mass spectrometry. J Lipid Res. 2013;54(6):1523–30.
Bujak R, Daghir-Wojtkowiak E, Kaliszan R, et al. PLS-based and regularization-based methods for the selection of relevant variables in non-targeted metabolomics data. Front Mol Biosci. 2016;3:35.
Cajka T, Fiehn O. Comprehensive analysis of lipids in biological systems by liquid chromatography-mass spectrometry. Trends Anal Chem. 2014;61:192–206.
Worley B, Powers R. PCA as a practical indicator of OPLS-DA model reliability. Curr Metab. 2016;4(2):97–103.
Dong S, Zhang R, Liang Y, et al. Changes of myocardial lipidomics profiling in a rat model of diabetic cardiomyopathy using UPLC/Q-TOF/MS analysis. Diabetol Metab Syndr. 2017;9(1):56.
Griffiths RL, Bunch J. A survey of useful salt additives in matrix-assisted laser desorption/ionization mass spectrometry and tandem mass spectrometry of lipids: introducing nitrates for improved analysis. Rapid Commun Mass Spectrom. 2012;26(13):1557–66.
Sugiura Y, Setou M. Selective imaging of positively charged polar and nonpolar lipids by optimizing matrix solution composition. Rapid Commun Mass Spectrom. 2009;23(20):3269–78.
Zehethofer N, Bermbach S, Hagner S, et al. Lipid analysis of airway epithelial cells for studying respiratory diseases. Chromatographia. 2015;78(5–6):403–13.
Satomi Y, Hirayama M, Kobayashi H. One-step lipid extraction for plasma lipidomics analysis by liquid chromatography mass spectrometry. J Chromatogr B. 2017;1063:93–100.
Lin X, Lu L, Liu L, et al. Blood lipids profile and lung cancer risk in a meta-analysis of prospective cohort studies. J Clin Lipidol. 2017;11(4):1073–81.
Santos CR, Schulze A. Lipid metabolism in cancer. Fed Eur Biochem Soc J. 2012;279(15):2610–23.
Kim D-S, Jeong S-K, Kim H-R, et al. Effects of triglyceride on ER stress and insulin resistance. Biochem Biophys Res Commun. 2007;363(1):140–5.
Romero F, Hong X, Shah D, et al. Lipid synthesis is required to resolve endoplasmic reticulum stress and limit fibrotic responses in the lung. Am J Respir Cell Mol Biol. 2018;59(2):225–36.
Agassandian M, Mallampalli RK. Surfactant phospholipid metabolism. Biochim Biophys Acta Mol Cell Biol Lipids. 2013;1831(3):612–25.
Low RB. Bronchoalveolar lavage lipids in idiopathic pulmonary fibrosis. Chest. 1989;95(1):3–5.
Bargagli E, Refini RM, d’Alessandro M, et al. Metabolic dysregulation in idiopathic pulmonary fibrosis. Int J Mol Sci. 2020;21(16):5663.
Nojima Y, Takeda Y, Maeda Y, et al. Metabolomic analysis of fibrotic mice combined with public RNA-Seq human lung data reveal potential diagnostic biomarker candidates for lung fibrosis. FEBS Open Bio. 2020;10(11):2427–36.
We would like to thank the co-ordinators of the Australian IPF registry (AIPFR) in each state for the collection of blood samples. Lung Foundation Australia has established the Australian IPF Registry with the generous support of unrestricted educational grant from Foundation Partners Roche Products Pty. Limited and Boehringer Ingelheim. The writers thank all participants and physicians who contribute to the Registry together with the Manager Sacha Macansh, Coordinators Amy Cashmore, Jessica Bucciarelli, Alysha Riley Karen Symons and Data Manager Faye Janice Lim.
The study is supported by NHMRC Grants (APP 1147776 and APP 1066128) and the Centre for Research Excellence in Pulmonary Fibrosis (APP 1099575).
Ethics approval and consent to participate
All work was approved by the Sydney Local Health District Human Research Ethics Committee (Reference number: HREC 11/RPAH 439), Royal Perth Hospital Human Research Ethics Committee (Reference number: REG 15-204) and the Murdoch University Human Research Ethics Committee (Approval number: 2017/254).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Nambiar, S., Clynick, B., How, B.S. et al. There is detectable variation in the lipidomic profile between stable and progressive patients with idiopathic pulmonary fibrosis (IPF). Respir Res 22, 105 (2021). https://doi.org/10.1186/s12931-021-01682-3