- Research
- Open access
- Published:
Proteomic associations with forced expiratory volume: a Mendelian randomisation study
Respiratory Research volume 25, Article number: 44 (2024)
Abstract
Background
A decline in forced expiratory volume (FEV1) is a hallmark of respiratory diseases that are an important cause of morbidity among the elderly. While some data exist on biomarkers that are related to FEV1, we sought to do a systematic analysis of causal relations of biomarkers with FEV1.
Methods
Data from the population-based AGES-Reykjavik study were used. Serum proteomic measurements were done using 4782 DNA aptamers (SOMAmers). Data from 1479 participants with spirometric data were used to assess the association of SOMAmer measurements with FEV1 using linear regression. Bi-directional two-sample Mendelian randomisation (MR) analyses were done to assess causal relations of observationally associated SOMAmers with FEV1, using genotype and SOMAmer data from 5368 AGES-Reykjavik participants and genetic associations with FEV1 from a publicly available GWAS (n = 400,102).
Results
In observational analyses, 530 SOMAmers were associated with FEV1 after multiple testing adjustment (FDR < 0.05). The most significant were Retinoic Acid Receptor Responder 2 (RARRES2), R-Spondin 4 (RSPO4) and Alkaline Phosphatase, Placental Like 2 (ALPPL2). Of the 257 SOMAmers with genetic instruments available, eight were associated with FEV1 in MR analyses. Three were directionally consistent with the observational estimate, Thrombospondin 2 (THBS2), Endoplasmic Reticulum Oxidoreductase 1 Beta (ERO1B) and Apolipoprotein M (APOM). THBS2 was further supported by a colocalization analysis. Analyses in the reverse direction, testing whether changes in SOMAmer levels were caused by changes in FEV1, were performed but no significant associations were found after multiple testing adjustments.
Conclusions
In summary, this large scale proteogenomic analyses of FEV1 reveals circulating protein markers of FEV1, as well as several proteins with potential causality to lung function.
Background
Chronic respiratory diseases such as chronic obstructive pulmonary disease (COPD) are a leading global cause of mortality and morbidity, with their relative importance increasing in the last decades [1]. Diagnosis of COPD is based on pulmonary function testing, by a low forced expiratory volume in one second (FEV1) relative to the forced vital capacity (FVC), and a progressive decline in pulmonary function is a feature of the disease [2]. Therefore, a decline in FEV1 is a hallmark of respiratory diseases that bring a great burden of disease to individuals and societies. While it is undisputed that exposure to external harmful stimuli such as cigarette smoke and biomass fumes are substantial risk factors for pulmonary function decline, intrinsic factors such as genetics and gene-environment interactions play a significant part as well [2,3,4]. Genome wide association studies (GWAS) have found several genetic polymorphisms that are associated with COPD or lung function decline, including polymorphisms in or near genes encoding matrix metalloproteinase 12, nicotinic acetylcholine receptor, hedgehog interacting protein, glutathione S-transferase, C-terminal domain–containing protein [5,6,7,8,9] and the antiprotease alpha-1-antitrypsin [10]. In addition to genetic polymorphisms, biomarkers that predict COPD or lung function have been discovered. Among them are inflammatory markers, soluble receptor for advanced glycoprotein end products (AGER), club cell secretory protein 16 (SCGB1A1) and surfactant protein D (SFTPD) [11,12,13,14]. For several of these biomarkers, estimations of potential causality have been made. AGER, SCGB1A1 and SFTPD have been suggested to be causally associated with COPD or with lung function, while analyses have pointed against such an association for the inflammatory markers CRP and IL-6 [15,16,17,18,19]. Still, while substantial epidemiologic data exist regarding single genetic markers and biomarkers that predict lung function or diseases whose diagnosis is based on impaired lung function, a systematic large-scale analysis of observational and causal associations of protein markers with lung function has not been undertaken to our knowledge.
Proteomics have emerged as a way of exploring molecular signatures of disease, especially with the advent of methods that allow for measurement and evaluation of thousands of proteins in biological samples from participants of large cohort studies [20]. In addition to predicting disease and disease-related outcomes, integration of genetic data allows one to make assumptions regarding the causality of proteomic markers. Mendelian randomization (MR) is such a method and utilizes genetic polymorphisms as instrumental variables to assess the relationship of an exposure with an outcome. As chromosomal alleles are randomly allocated during gamete formation, this methodology allows one to avoid the effect of confounders and to infer causality from epidemiologic data [21].
The aim of the study was to systematically explore the associations of a multitude of protein markers with lung function in an elderly population, focusing on FEV1 as the main outcome. Then, the aim was to assess the potential causal relationships of these protein markers with FEV1 by use of bi-directional two-sample Mendelian randomization.
Methods
Study phenotyping
The Age/Gene Environment Susceptibility (AGES)-Reykjavik study is a population-based cohort study of 5,764 elderly Icelanders that was carried out between 2002 and 2006. The participants, aged from 66 to 96 years (mean 76 years), were all prior participants of the Reykjavik Study done decades earlier. As part of AGES-Reykjavik, the participants underwent extensive phenotyping by questionnaires, physiological measurements, imaging studies and laboratory measurements, during a three-day period. The study was approved by the Icelandic National Bioethics Committee (VSN-00-063), in accordance with the Helsinki Declaration and the Institutional Review Board of the Intramural Program of the National Institute for Aging, with informed consent obtained from all participants. Further details of the study design are previously published [22].
A random subset of the study participants underwent lung function testing in a standardised manner. The device used was a Vitalograph Gold Standard Plus (Vitalograph Ltd., UK). Each participant completed three attempts. Participants with spirometry of acceptable quality were included. Spirometry measures from participants that completed at least two attempts with a no more than 300 ml difference between the attempts and exhalation for at least 6 s were deemed acceptable [23]. Smoking history was ascertained from questionnaires while anthropometric measurements were done during the clinic visit. Protein measurements were done in serum samples from participants using a high throughput proteomics technology, the SOMAscan (SomaLogic, Boulder, CO) in which DNA aptamers (Slow-Off Rate Modified Aptamers (SOMAmers)) bind to target protein epitopes and are then quantified with the help of fluorescence after wash-out of unbound aptamers and proteins. Measurements in AGES-Reykjavik were done with a 5034 SOMAmer platform in serum from 5457 AGES-Reykjavik participants. For the analyses, 4782 SOMAmers targeting 4135 human proteins were used, excluding SOMAmers annotated to non-human proteins. Measurement data were transformed using Box-Cox transformation and extreme outliers were excluded, as previously described [24].
Statistical analyses
A flow chart of study design is shown in Fig. 1. Descriptive statistics were compiled for participants with acceptable pulmonary function tests, demographic covariate data and protein measurements available (n = 1479). Using data from these participants, the association of all measured human SOMAmers (n = 4782) with FEV1 was assessed using linear regression modelling. These analyses were adjusted for variables that are commonly used to predict FEV1 in clinical practice: age, sex, height, age squared, and height squared [25]. Adjustment for multiple testing was done using a Benjamini–Hochberg False Discovery Rate (FDR), where FDR < 0.05 was considered statistically significant. To understand how history of tobacco smoking, the most important lifestyle factor influencing lung function, affected these associations, the analyses were repeated with participants stratified by smoking history (ever-smokers versus never-smokers). In a secondary analysis, all SOMAmers associated (FDR < 0.05) with FEV1 were tested for association with FEV1/FVC ratio and a FEV1/FVC ratio under 0.7, a value used in the clinical diagnosis of COPD [2], using a linear and logistic regression, respectively. These analyses were adjusted for the same covariates as the primary analysis for FEV1. Genes encoding proteins significantly associated with FEV1 after adjustment for multiple testing (FDR < 0.05) were subjected to over-representation analysis of Gene Ontology terms [26].
For SOMAmers that were associated with FEV1 after adjustment for multiple testing (FDR < 0.05), a bi-directional two-sample MR analysis was performed to assess support for causal relationships between SOMAmers and FEV1, in either direction. When testing the causal effect of SOMAmer on FEV1 (forward direction), the MR analysis was restricted to SOMAmers with genetic instruments available, which were defined as follows. The associations of single-nucleotide polymorphisms (SNPs) with SOMAmers were calculated using data from 5,368 participants of AGES-Reykjavik, as previously published and described in detail [24]. The possible instrumental variables for each SOMAmer were defined as SNPs located within a cis-window for the gene coding for the protein measured by the SOMAmer, defined as within 500 kb up- and downstream of the gene. SNPs with a window-wide significant association (p < 0.05/number of SNPs in cis-window) with a given SOMAmer were considered as potential instruments. For each gene, SNPs were filtered based on linkage disequilibrium (LD; r2 < 0.2) or distance (> 1 mb) using Plink v1.9. As the forward MR analysis was restricted to a single region on the genome, a more inclusive LD threshold compared to the reverse MR was chosen to increase statistical power [27]. Instruments were considered valid if F > 10 for the association of instruments with SOMAmers [28]. The associations of the instrumental variables with FEV1 were obtained from a GWAS of lung function [29] in which data from two cohorts, UKBiobank and SpiroMeta Consortium, were meta-analysed. The total number of participants in that analysis was 400,102 [29]. A MR analysis in the reverse direction to assess the causal effect of FEV1 on SOMAmer levels was performed using all SOMAmers associated with FEV1 in the observational linear regression model (FDR < 0.05) as outcomes. Genetic instruments for FEV1 were selected based on the same GWAS of lung function as described above [29]. For these analyses, SNPs associated with FEV1 (p < 5 × 10–8) were defined as instruments after clumping by an LD threshold of r2 < 0.01. The associations of the instrumental variables with SOMAmer levels were obtained from AGES-Reykjavik results [24]. For the MR analyses, SNP reference alleles were harmonised between studies using the TwoSampleMR R package [30]. The MR estimates for FEV1 were obtained using the generalized weighted least squares (GWLS) method [31], which accounts for correlation between instruments, except for instances where only one instrument was available, in which Wald ratios were calculated. For each MR analysis with more than two instruments, sensitivity analyses using the weighted median and Egger estimators were done to assess the validity of instruments and limit the effect of pleiotropic associations, respectively. Results were considered to pass these sensitivity analyses when the following conditions were met. For the weighted median estimator, the weighted median estimate had to be significant and directionally consistent with the GWLS estimate and for Egger, the Egger estimate had to be directionally consistent with the GWLS estimate and the intercept not significant. Additionally, for SOMAmers associated with FEV1 in the forward MR analysis, a leave-one-out analysis was performed. For analyses in the reverse direction, estimates were obtained using the inverse variance weighted method. Adjustment for multiple testing was done using the Benjamini–Hochberg False Discovery Rate (FDR).
A colocalization analysis was performed to provide additional causal support for analytes associated with FEV1 in MR analyses. Here, AGES-Reykjavik serum protein quantitative trait loci (pQTLs) [24], summary statistics for FEV1 [29] and plasma pQTLs from a cohort study of over 35,000 Icelanders designed to associate genetics, proteins and disease, published by deCODE Genetics [32] were harmonized to account for strand orientation and differences in genome builds. All studies were lifted over to build GRCh38 for colocalization when needed. For the QTLs, only putative cis-regulatory variants were examined defined by 500kb away from the gene body. All regions with variants that had associations of P < 1 × 10–5 were fine-mapped using Sum of Single Effects (SuSiE) [33] based on 1000 genome population reference [34]. 95% credible-sets were filtered and colocalization between GWAS and QTLs was performed using fastENLOC using the posterior inclusion probabilities estimated in SuSiE [33, 35]. Regional colocalization probability (RCP) was calculated by summing the colocalization probability within each 95% credible-set and filtered colocalization results at RCP > 80% for further examination. RCP represents the probability that a given genomic region contains a single colocalized variant [36]. Visualization was done using Julia 1.7 and libraries within [37].
Results
The protein profile of pulmonary function
Descriptive statistics for participants that had pulmonary function testing data are shown in Table 1. A majority of participants were female (54%) and participants were on average 76 years old. Most participants had a history of smoking (60%) while a minority had evidence of obstruction on spirometry (37%).
Of the 4,782 SOMAmers tested, 530 were observationally associated (FDR < 0.05) with FEV1 after adjustment for multiple testing (Table 2, Additional file 5: Table S1, Fig. 2). The most significantly associated SOMAmers measured were Retinoic Acid Receptor Responder 2 (RARRES2, β = − 0.103, p = 7.76 × 10–12), R-Spondin 4 (RSPO4, β = − 0.094, p = 1.86 × 10–11), Alkaline Phosphatase, Placental Like 2 (ALPPL2, β = − 0.087, p = 1.23 × 10–10), Complement C9 (C9, β = − 0.084, p = 3.04 × 10–10) and Hematopoietic Prostaglandin D Synthase (HPGDS, β = 0.083, p = 5.12 × 10–10). Results for proteins that have been previously suggested as biomarkers of FEV1 [11, 13] are shown in Additional file 6: Table S2. Prior associations of SCGB1A1 (β = − 0.034, p = 0.02), SFTPD (β = − 0.041, p = 1.67 × 10–3), CRP (β = − 0.065, p = 2.72 × 10–7), fibrinogen (β = − 0.041, p = 2.8 × 10–3 for the stronger associated SOMAmer), IL6 (β = − 0.054, p = 1.5 × 10–4 for the stronger associated SOMAmer), eotaxin (β = − 0.051, p = 1.1 × 10–4) and TNF (β = − 0.028, p = 0.04 for the stronger associated SOMAmer) were reproduced in the AGES-Reykjavik data while associations for other suggested biomarkers of FEV1 were not.
The proteins most significantly associated with FEV1 overall also had strong associations with FEV1 among ever-smokers, such as RARRES2 (β = − 0.12, p = 5.85 × 10–9), RSPO4 (β = − 0.095, p = 9.61 × 10–7), ALPPL2 (β = − 0.097, p = 1.05 × 10–7) and Complement 9 (C9, β = − 0.096, p = 2.05 × 10–7). Trefoil Factor 2 (TFF2, β = − 0.104, p = 3.37 × 10–8), and Neurotrophic Receptor Tyrosine Kinase 3 (NTRK3, β = 0.107, p = 2.47 × 10–8) also had notably strong associations with the outcome. However, associations of these proteins were much weaker among never-smokers, with the strongest associations observed for two SOMAmers measuring SVEP1 (β = − 0.082, p = 1.64 × 10–5 for the stronger associated SOMAmer) (Additional file 6: Table S3) in this subgroup. Of the 530 proteins associated with FEV1, 224 were also associated (FDR < 0.05) with FEV1/FVC (42%) and 160 (30%) were associated with an obstructive deficit on spirometry (FEV1/FVC < 0.7; Additional file 1: Figure S1, Additional file 5: Table S1). Among the 25 proteins most strongly associated with FEV1, including RARRES2, RSPO4 and ALPPL2, 19 were associated with continuous FEV1/FVC and 16 were associated with FEV1/FVC under 0.7 (Additional file 6: Table S4). The 530 proteins associated with FEV1 in linear regression analyses were most strongly enriched for Gene Ontology terms related to compliment activation, extracellular matrix organization, peptidase regulator activity, humoral immune response and regulation of neurogenesis (Additional file 2: Figure S2, Additional file 5: Table S5).
Mendelian randomization analysis
Of the 530 SOMAmers observationally associated with FEV1, 257 (49%) had genetic instruments available (Additional file 5: Table S1) and were included in MR analyses to evaluate potentially causal associations of SOMAmers with FEV1. The instruments are shown in Additional file 5: Table S6 for SOMAmers with nominally significant (unadjusted p < 0.05) MR associations. Table 3 shows the 35 SOMAmers that were nominally associated with FEV1 (unadjusted p < 0.05) in the MR analysis and passed weighted median sensitivity testing. Eight were significantly associated with FEV1 in the MR analysis (FDR for MR estimate < 0.05; Table 3), suggesting they may have a causal effect on lung function. Of the seven associations based on more than two genetic instruments, none were solely driven by a single variant (Additional file 4: Figure S4). Three SOMAmers that were significant (FDR < 0.05) in the MR analysis, Thrombospondin 2 (THBS2, β = − 0.037, p = 9.53 × 10–5), Endoplasmic Reticulum Oxidoreductase 1 Beta (ERO1B, β = − 0.025, p = 8.05 × 10–4) and Apolipoprotein M (APOM, β = 0.053, p = 9.72 × 10–4) were directionally consistent with the observational analyses. The other five SOMAmers with significant causal estimates, R-Spondin-2 (RSPO2), TIMP Metallopeptidase Inhibitor 4 (TIMP4), interleukin 1 receptor antagonist (IL1RN), CD14 and Heparin Binding Growth Factor (HDGF) were directionally inconsistent (Fig. 3). Among all 35 nominally associated SOMAmers, the directional consistency between causal and observational estimates was low, or 34%, suggesting that the potentially causal effects of the proteins are generally not reflected in the observational estimates. However, these nominally associated SOMAmers included SERPINA1, which measures the level of alpha-1-antitrypsin, deficiency of which is a well-established causal determinant of impaired lung function via emphysema formation [38] (Table 3). Of the five previously suggested biomarkers of FEV1 listed in Additional file 6: Table S2 that were associated with FEV1 in linear regression analyses (FDR < 0.05), only three proteins had genetic instruments available, the acute phase reactants CRP and fibrinogen as well as eotaxin. None of these proteins were associated with FEV1 in the MR analysis (Additional file 6: Table S7).
The 35 proteins representing the 35 SOMAmers identified from the MR analysis (unadjusted p < 0.05; Table 3) were further examined for colocalization evidence between genetic associations for protein levels [24, 32] and FEV1 [29]. Prioritizing the 8 proteins with significant (FDR < 0.05) MR association, we found strong colocalization support (RCP > 0.9) for a single protein, THBS2 (Table 4, Fig. 4, Additional file 5: Table S8). The genetic cis-signal for THBS2 protein expression in AGES-Reykjavik (rs3253 β = − 0.18, p = 1.4 × 10–19) colocalized with a signal for FEV1, with the 3’ UTR variant rs3253 being the strongest shared variant (RCP = 0.9976). The THBS2 protein association for this variant replicated in the deCODE cohort (p = 2.1 × 10–186, β = − 0.26) (Table 4, Fig. 4). Among the remaining proteins with only nominal associations (p < 0.05) in the MR analysis, TNFSF12 had colocalization support where an upstream variant (rs4968200) associated with TNFSF12 protein levels in AGES-Reykjavik (β = 0.63, p = 5.6 × 10–133) colocalized with a signal for FEV1 (RCP = 0.9995; Additional file 3: Figure S3). The other proteins had colocalization probability less than 80%.
Finally, all 530 SOMAmers observationally associated with FEV1 were included in a MR analysis in the reverse direction, i.e., to evaluate if the changes in SOMAmer levels are downstream of changes in lung function. FEV1 was not causally associated with levels of any SOMAmers after adjustment for multiple testing. Data for the 25 SOMAmers that FEV1 was nominally associated with (p < 0.05) in the reverse-MR analysis are shown in Additional file 6: Table S9.
Discussion
We present findings from a proteomic analysis of pulmonary function with more candidate protein analytes than previously published to our knowledge [11], highlighting several proteins as strong markers of FEV1. Stratification by smoking shows that most associations are driven by ever-smokers. Mendelian randomisation was systematically applied to the candidate markers, revealing proteins whose levels may have a causal effect on lung function. Of those, probabilistic colocalization supported a role of THBS2 and TNFSF12 in affecting lung function. Reverse causation analyses failed to demonstrate that protein level changes associated with FEV1 occur downstream of the phenotype change, although this could partly be due to insufficient power.
Eight proteins were suggested to be causally implicated in lung function based on the MR analysis. However, only three (THBS2, ERO1B and APOM) had consistent direction of effect for the observational and causal estimates. Such discrepancy has been observed when comparing causal and observational estimates [39, 40] for serum proteins. Based on probabilistic quantification, a 3’ UTR variant within THBS2 was identified as a putative causal variant affecting FEV1 and this lead variant colocalized with THBS2 protein expression in AGES-Reykjavik. The effect of this variant was replicated in an independent cohort for THBS2 protein levels. Because directions of effects were consistent across the datasets and THBS2 had non-revertible causal association to FEV1, supported by colocalization, THBS2 may have biological importance in impaired lung function in the elderly, and could represent a therapeutic target for some forms of respiratory disease. THBS2 is an extracellular matrix protein that has been implicated in various cardiovascular disorders and is also a candidate biomarker for non-small cell lung cancer [41, 42]. THBS2 is involved in tissue repair and interacts with many different ligands in the extracellular matrix, among them matrix metalloproteases and elastase [43]. Although mechanism of THBS2 needs further experimental validation, it is possible that protein levels of THBS2 may influence lung function via extracellular matrix and regenerative pathways. Meanwhile, ERO1B is a disulfide oxidase in the endoplasmic reticulum that is shown to predict survival in pancreatic and pulmonary cancers [44,45,46] and APOM is an apolipoprotein that is mainly a component of high density lipoproteins and has been associated with COPD severity [47]. Genetic variants flanking the APOM gene have been associated with obstructive spirometry measurements [48]. However, it must be kept in mind that the causal association of APOM in our study is based on a single SNP (rs2736163, intronic to PRRC2A), thus complicating the interpretation of the MR results.
Despite not reaching the study threshold for statistical significance, some proteins that were nominally associated with lung function in the MR analysis are of interest. First, alpha-1-antitrypsin (SERPINA1) is the best-known protein known to cause COPD as severe deficiency of alpha-1-antitrypsin results in obstructive lung disease [38]. However, in our data, serum levels of SERPINA1 are inversely associated with FEV1 in both observational and MR analyses, contrary to what would possibly be expected, although this directionality is known from previous observational analyses of FEV1 and explained by alpha-1-antitrypsin’s role as an acute phase reactant [49]. Also, polymorphisms that cause mild or intermediate alpha-1-antitrypsin deficiency are not consistently associated with decreased lung function, suggesting that levels of the protein may only affect lung function below a threshold level [50]. Second, TNFSF12 was a protein with a colocalizing pQTL variant with the FEV1 GWAS (Table 4, Fig. 4). TNFSF12 is a member of the Tumor Necrosis Factor (TNF) superfamily, of which one key cytokine, TNF-⍺, is a well-known protein whose levels are disrupted in COPD patients [51]. Presented here are MR based causal evidence and probabilistic colocalization findings that suggest that diseases that present with impaired lung function could be impacted by TNF-⍺ associated pathobiology via TNFSF12. Notable other nominally significant or directionally inconsistent proteins in the findings are matrix metalloproteinase 8 (MMP-8), one of the proteases observationally associated with lung function and implicated in the pathogenesis of COPD [52, 53], TIMP4, an inhibitor of matrix metalloproteinases that has been shown to be upregulated in COPD patients [54] and CD14, levels of which have been shown to be elevated in lungs of smokers [55]. While only eight proteins were significantly associated with FEV1 in MR analyses after multiple testing, FEV1 was not associated with any proteins in analyses in the reverse direction. The data were therefore unable to support prior causal analyses involving inflammatory markers that have suggested reverse causality of FEV1 with inflammatory markers [18] (Additional file 6: Table S7).
The study reveals novel markers with strong observational relations to lung function such as RSPO4, a signalling molecule that is part of the Wnt signalling pathway [56], the tumor marker ALPPL2 [57], the adipokine chemerin (RARRES2), and SVEP1, a protein that is thought to play a role in inflammation in atherosclerosis [58]. In addition, the findings validate the observational associations of some of the previously suggested protein markers of FEV1, such as SFTPD, fibrinogen, IL6, eotaxin and CRP (Additional file 6: Table S2) [13], while the associations of other previously suggested biomarkers of FEV1 such as AGER were not corroborated in this study [13, 59]. Secondary analyses showed that most SOMAmers with the strongest associations with FEV1 were also associated with FEV1/FVC and/or FEV1/FVC under 0.7 in a directionally consistent manner (Additional file 1: Figure S1 and Additional file 6: Table S4). However, under half of all 530 FEV1-associated SOMAmers had statistically significant (FDR < 0.05) associations with the secondary outcomes (Additional file 5: Table S1). Also notable is the finding that among the eight SOMAmers with support for a causal effect on FEV1 from MR analyses (FDR < 0.05), three (RSPO2, TIMP4 and CD14) were observationally associated with FEV1/FVC and/or FEV1/FVC under 0.7, while the remaining five (APOM, THBS2, ERO1B, HDGF, IL1RN) were not (Additional file 5: Table S1). Collectively, these results from secondary analyses suggest that some SOMAmer associations with FEV1 may be explained by other disease mechanisms than obstructive pulmonary disease, such as lung aging. Finally, our findings show that proteins that take part in immune responses, peptidase regulation and extracellular matrix modulation are over-represented among the proteins related to FEV1.
This work is subject to a number of limitations. First, the study is based on SOMAmer technology, a relatively novel aptamer-based technology for protein measurements. While many of these SOMAmers have been validated with encouraging results, it has been pointed out that a minority of SOMAmers could be subject to cross-reactivity with related or homologous proteins [20]. Second, both the AGES-Reykjavik cohort and the UKBiobank and SpiroMeta Consortium are of European ancestry [29]. Therefore, the observed associations could not be generalizable to other populations. Additionally, the AGES-Reykjavik cohort is older than the UKBiobank and SpiroMeta Consortium which could distort comparisons between observational and genetic findings. Third, the processes which mediate the association of SOMAmers with lung function cannot be elucidated from these results. The measure of lung function used in this paper, FEV1, is disproportionally impaired in obstructive lung disease such as COPD. Still, less than half of the FEV1-associated SOMAmers were associated with spirometric obstruction, as discussed above. So, while some associations of SOMAmers with FEV1 may reflect obstructive lung disease, this is likely not the case for all associated SOMAmers. Fourth, subtle differences in genetic structure and datasets between AGES-Reykjavik, UKBiobank, and SpiroMeta cohorts are present and may be contributing to the lack of colocalization for some of the MR identified proteins. For instance, very few MR instruments overlapped with 95% credible-sets identified (Additional file 5: Table S8). In comparisons to UKBiobank, missing variants in the AGES-Reykjavik cohort may contribute to colocalization false negatives. Lastly, many causal estimates are directionally inconsistent with observational estimates. While this phenomenon is known from previous proteogenomic studies, its reasons are unclear.
Conclusions
In conclusion, this proteogenomic analysis reveals several proteins that are potentially causally related to lung function, most notably THBS2, ERO1B and APOM.
Availability of data and materials
While study materials are not publicly available due to participant privacy, further results or data generated in this study are available upon reasonable request to the authors.
References
Vos T, Lim SS, Abbafati C, Abbas KM, Abbasi M, Abbasifard M, et al. Global burden of 369 diseases and injuries in 204 countries and territories, 1990–2019: a systematic analysis for the Global Burden of Disease Study 2019. The Lancet. 2020;396(10258):1204–22.
Vestbo J, Hurd SS, Agustí AG, Jones PW, Vogelmeier C, Anzueto A, et al. Global strategy for the diagnosis, management, and prevention of chronic obstructive pulmonary disease. Am J Respir Crit Care Med. 2013;187(4):347–65.
Rennard SI, Vestbo J. COPD: the dangerous underestimate of 15%. The Lancet. 2006;367(9518):1216–9.
Gordon SB, Bruce NG, Grigg J, Hibberd PL, Kurmi OP, Lam KB, et al. Respiratory risks from household air pollution in low and middle income countries. Lancet Respir Med. 2014;2(10):823–60.
Hunninghake GM, Cho MH, Tesfaigzi Y, Soto-Quiros ME, Avila L, Lasky-Su J, et al. MMP12, lung function, and COPD in high-risk populations. N Engl J Med. 2009;361(27):2599–608.
Cho MH, Boutaoui N, Klanderman BJ, Sylvia JS, Ziniti JP, Hersh CP, et al. Variants in FAM13A are associated with chronic obstructive pulmonary disease. Nat Genet. 2010;42(3):200–2.
Pillai SG, Ge D, Zhu G, Kong X, Shianna KV, Need AC, et al. A genome-wide association study in chronic obstructive pulmonary disease (COPD): identification of two major susceptibility loci. PLoS Genet. 2009;5(3): e1000421.
Soler Artigas M, Wain LV, Repapi E, Obeidat M, Sayers I, Burton PR, et al. Effect of five genetic variants associated with lung function on the risk of chronic obstructive lung disease, and their joint effects on lung function. Am J Respir Crit Care Med. 2011;184(7):786–95.
Cho MH, McDonald M-LN, Zhou X, Mattheisen M, Castaldi PJ, Hersh CP, et al. Risk loci for chronic obstructive pulmonary disease: a genome-wide association study and meta-analysis. Lancet Respir Med. 2014;2(3):214–25.
Wood AM, Stockley RA. Alpha one antitrypsin deficiency: from gene to treatment. Respiration. 2007;74(5):481–92.
Serban KA, Pratte KA, Bowler RP. Protein biomarkers for COPD outcomes. Chest. 2021;159(6):2244–53.
Faner R, Tal-Singer R, Riley JH, Celli B, Vestbo J, MacNee W, et al. Lessons from ECLIPSE: a review of COPD biomarkers. Thorax. 2014;69(7):666.
Zemans RL, Jacobson S, Keene J, Kechris K, Miller BE, Tal-Singer R, et al. Multiple biomarkers predict disease severity, progression and mortality in COPD. Respir Res. 2017;18(1):117.
Gan WQ, Man SFP, Senthilselvan A, Sin DD. Association between chronic obstructive pulmonary disease and systemic inflammation: a systematic review and a meta-analysis. Thorax. 2004;59(7):574.
Keefe J, Yao C, Hwang S-J, Courchesne P, Lee GY, Dupuis J, et al. An integrative genomic strategy identifies sRAGE as a causal and protective biomarker of lung function. Chest. 2022;161(1):76–84.
Milne S, Li X, Hernandez Cordero AI, Yang CX, Cho MH, Beaty TH, et al. Protective effect of club cell secretory protein (CC-16) on COPD risk and progression: a Mendelian randomisation study. Thorax. 2020;75(11):934.
Obeidat M, Li X, Burgess S, Zhou G, Fishbane N, et al. Surfactant protein D is a causal risk factor for COPD: results of Mendelian randomisation. Eur Respir J. 2017;50(5):1700657.
Dahl M, Vestbo J, Zacho J, Lange P, Tybjærg-Hansen A, Nordestgaard BG. C reactive protein and chronic obstructive pulmonary disease: a Mendelian randomisation approach. Thorax. 2011;66(3):197.
van Durme YMTA, Lahousse L, Verhamme KMC, Stolk L, Eijgelsheim M, Loth DW, et al. Mendelian randomization study of interleukin-6 in chronic obstructive pulmonary disease. Respiration. 2011;82(6):530–8.
Emilsson V, Ilkov M, Lamb JR, Finkel N, Gudmundsson EF, Pitts R, et al. Co-regulatory networks of human serum proteins link genetics to disease. Science. 2018;361(6404):769–73.
Emdin CA, Khera AV, Kathiresan S. Mendelian randomization. JAMA. 2017;318(19):1925–6.
Harris TB, Launer LJ, Eiriksdottir G, Kjartansson O, Jonsson PV, Sigurdsson G, et al. Age, gene/environment susceptibility—Reykjavik study: multidisciplinary applied phenomics. Am J Epidemiol. 2007;165(9):1076–87.
Gudmundsson G, Margretardottir OB, Sigurdsson MI, Harris TB, Launer LJ, Sigurdsson S, et al. Airflow obstruction, atherosclerosis and cardiovascular risk factors in the AGES Reykjavik study. Atherosclerosis. 2016;252:122–7.
Gudjonsson A, Gudmundsdottir V, Axelsson GT, Gudmundsson EF, Jonsson BG, Launer LJ, et al. A genome-wide association study of serum proteins reveals shared loci with common diseases. Nat Commun. 2022;13(1):480.
Hankinson JL, Odencrantz JR, Fedan KB. Spirometric reference values from a sample of the general U.S. population. Am J Respir Crit Care Med. 1999;159(1):179–87.
Boyle EI, Weng S, Gollub J, Jin H, Botstein D, Cherry JM, et al. GO::TermFinder–open source software for accessing Gene Ontology information and finding significantly enriched Gene Ontology terms associated with a list of genes. Bioinformatics. 2004;20(18):3710–5.
Gkatzionis A, Burgess S, Newcombe PJ. Statistical methods for cis-Mendelian randomization with two-sample summary-level data. Genet Epidemiol. 2023;47(1):3–25.
Sanderson E, Glymour MM, Holmes MV, Kang H, Morrison J, Munafò MR, et al. Mendelian randomization. Nat Rev Methods Primers. 2022;2.
Shrine N, Guyatt AL, Erzurumluoglu AM, Jackson VE, Hobbs BD, Melbourne CA, et al. New genetic signals for lung function highlight pathways and chronic obstructive pulmonary disease associations across multiple ancestries. Nat Genet. 2019;51(3):481–93.
Hemani G, Zheng J, Elsworth B, Wade KH, Haberland V, Baird D, et al. The MR-Base platform supports systematic causal inference across the human phenome. Elife. 2018;7:e34408.
Burgess S, Dudbridge F, Thompson SG. Combining information on multiple instrumental variables in Mendelian randomization: comparison of allele score and summarized data methods. Stat Med. 2016;35(11):1880–906.
Ferkingstad E, Sulem P, Atlason BA, Sveinbjornsson G, Magnusson MI, Styrmisdottir EL, et al. Large-scale integration of the plasma proteome with genetics and disease. Nat Genet. 2021;53(12):1712–21.
Zou Y, Carbonetto P, Wang G, Stephens M. Fine-mapping from summary data with the “Sum of Single Effects” model. PLoS Genet. 2022;18(7): e1010299.
Auton A, Brooks LD, Durbin RM, Garrison EP, Kang HM, Korbel JO, et al. A global reference for human genetic variation. Nature. 2015;526(7571):68–74.
Wen X, Pique-Regi R, Luca F. Integrating molecular QTL data into genome-wide genetic association analysis: probabilistic assessment of enrichment and colocalization. PLoS Genet. 2017;13(3): e1006646.
Hukku A, Pividori M, Luca F, Pique-Regi R, Im HK, Wen X. Probabilistic colocalization of genetic variants from complex and molecular traits: promise and limitations. Am J Hum Genet. 2021;108(1):25–35.
Kim M, Vo DD, Kumagai ME, Jops CT, Gandal MJ. GeneticsMakie.jl: a versatile and scalable toolkit for visualizing locus-level genetic and genomic data. Bioinformatics. 2023;39(1).
Strnad P, McElvaney NG, Lomas DA. Alpha1-antitrypsin deficiency. N Engl J Med. 2020;382(15):1443–55.
Sun BB, Maranville JC, Peters JE, Stacey D, Staley JR, Blackshaw J, et al. Genomic atlas of the human plasma proteome. Nature. 2018;558(7708):73–9.
Gudmundsdottir V, Zaghlool SB, Emilsson V, Aspelund T, Ilkov M, Gudmundsson EF, et al. Circulating protein signatures and causal candidates for type 2 diabetes. Diabetes. 2020;69(8):1843–53.
Zhang K, Li M, Yin L, Fu G, Liu Z. Role of thrombospondin-1 and thrombospondin-2 in cardiovascular diseases (Review). Int J Mol Med. 2020;45(5):1275–93.
Jiang YM, Yu DL, Hou GX, Jiang JL, Zhou Q, Xu XF. Serum thrombospondin-2 is a candidate diagnosis biomarker for early non-small-cell lung cancer. Biosci Rep. 2019;39(7).
Calabro NE, Kristofik NJ, Kyriakides TR. Thrombospondin-2 and extracellular matrix assembly. Biochim Biophys Acta. 2014;1840(8):2396–402.
Zito E, Chin KT, Blais J, Harding HP, Ron D. ERO1-beta, a pancreas-specific disulfide oxidase, promotes insulin biogenesis and glucose homeostasis. J Cell Biol. 2010;188(6):821–32.
Asada K, Kobayashi K, Joutard S, Tubaki M, Takahashi S, Takasawa K, et al. Uncovering prognosis-related genes and pathways by multi-omics analysis in lung cancer. Biomolecules. 2020;10(4):524.
Zhu T, Gao YF, Chen YX, Wang ZB, Yin JY, Mao XY, et al. Genome-scale analysis identifies GJB2 and ERO1LB as prognosis markers in patients with pancreatic cancer. Oncotarget. 2017;8(13):21281–9.
Li H, Liu Y, Wang L, Shen T, Du W, Liu Z, et al. High apolipoprotein M serum levels correlate with chronic obstructive pulmonary disease. Lipids Health Dis. 2016;15:59.
Burkart KM, Manichaikul A, Wilk JB, Ahmed FS, Burke GL, Enright P, et al. APOM and high-density lipoprotein cholesterol are associated with lung function and per cent emphysema. Eur Respir J. 2014;43(4):1003–17.
Senn O, Russi EW, Schindler C, Imboden M, von Eckardstein A, Brändli O, et al. Circulating alpha1-antitrypsin in the general population: determinants and association with lung function. Respir Res. 2008;9(1):35.
Thun G-A, Ferrarotti I, Imboden M, Rochat T, Gerbase M, Kronenberg F, et al. SERPINA1 PiZ and PiS heterozygotes and lung function decline in the SAPALDIA cohort. PLoS ONE. 2012;7(8): e42728.
Tanni SE, Pelegrino NR, Angeleli AY, Correa C, Godoy I. Smoking status and tumor necrosis factor-alpha mediated systemic inflammation in COPD patients. J Inflamm (Lond). 2010;7:29.
Koo HK, Hong Y, Lim MN, Yim JJ, Kim WJ. Relationship between plasma matrix metalloproteinase levels, pulmonary function, bronchodilator response, and emphysema severity. Int J Chron Obstruct Pulmon Dis. 2016;11:1129–37.
Churg A, Zhou S, Wright JL. Matrix metalloproteinases in COPD. Eur Respir J. 2012;39(1):197.
Navratilova Z, Zatloukal J, Kriegova EVA, Kolek V, Petrek M. Simultaneous up-regulation of matrix metalloproteinases 1, 2, 3, 7, 8, 9 and tissue inhibitors of metalloproteinases 1, 4 in serum of patients with chronic obstructive pulmonary disease. Respirology. 2012;17(6):1006–12.
Regueiro V, Campos MA, Morey P, Sauleda J, Agustí AGN, Garmendia J, et al. Lipopolysaccharide-binding protein and CD14 are increased in the bronchoalveolar lavage fluid of smokers. Eur Respir J. 2009;33(2):273.
de Lau WB, Snel B, Clevers HC. The R-spondin protein family. Genome Biol. 2012;13(3):242.
Su Y, Zhang X, Bidlingmaier S, Behrens CR, Lee NK, Liu B. ALPPL2 is a highly specific and targetable tumor cell surface antigen. Cancer Res. 2020;80(20):4552–64.
Jung IH, Elenbaas JS, Alisio A, Santana K, Young EP, Kang CJ, et al. SVEP1 is a human coronary artery disease locus that promotes atherosclerosis. Sci Transl Med. 2021;13(586).
Merali S, Barrero CA, Bowler RP, Chen DE, Criner G, Braverman A, et al. Analysis of the plasma proteome in COPD: novel low abundance proteins reflect the severity of lung remodeling. COPD J Chron Obstruct Pulmon Dis. 2014;11(2):177–89.
Acknowledgements
The authors thank all IHA staff that contributed to data acquisition but were not directly involved in the design or writing of this study. In addition, the authors wholeheartedly thank all participants of the AGES-Reykjavik Study for their participation.
Funding
The AGES-Reykjavik study was funded by Icelandic Heart Association, National Institute on Aging contract N01-AG-12100 and HHSN271201200022C, and Althingi (the Icelandic Parliament). GTA was supported by the Eimskip University Fund from the University of Iceland. GG was supported by the University of Iceland and by Landspitali University Hospital. TJ was supported by the Icelandic Research fund (Grant no. 206692-052) and VaG by the University of Iceland postdoctoral fund. All protein measurements were funded by Novartis Biomedical Research.
Author information
Authors and Affiliations
Contributions
Conceptualization, methodology, and design of the study: GTA, TJ, VaG, GG, VE and ViG. Data acquisition: TA, LLJ, JJL, APO, VE, ViG. Analysis of the data: GTA, TJ, VaG, EAF, TA, YJW, JJL. Administration and supervision: VaG, LLJ, GG, VE, ViG. Drafting of the initial manuscript: GTA, TJ, VaG and ViG. All authors revised and/or edited the manuscript for scientific content.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
Approval for the AGES-Reykjavik study was obtained from the National Bioethics Committee in Iceland that acts as the Institutional Review Board for the Icelandic Heart Association (approval number: VSN-00-063), and by the National Institute on Aging Intramural Institutional Review Board.
Consent for publication
Not applicable.
Competing interests
JJL, YJW and LLJ are employees and stockholders of Novartis. GTA and GG report travel support from Boehringer-Ingelheim for work not related to the study. No other potential conflicts of interest relevant to this article were reported.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1:
Protein associations with FEV1 (linear regression) compared with (A) continuous FEV1/FVC (linear regression) and (B) FEV1/FVC under 0.70 (logistic regression).
Additional file 2:
Over-representation analyses of Gene Ontology (GO) terms associated with genes annotated to FEV1-associated SOMAmers.
Additional file 3:
Colocalization plot for TNFSF12 protein levels and FEV1.
Additional file 4:
The results of a leave-one-out analysis for the seven proteins that had significant (FDR < 0.05) causal estimates for FEV1 in the MR analysis and had three or more SNPs as instruments (A: THBS2; B: ILRN; C: TIMP4; D: ERO1B; E: RSPO2; F: HDGF; G: CD14). The original causal estimate is shown in red. Each remaining x- and y-axis pair represents a causal estimate and its standard error evaluated without the listed SNP.
Additional file 5:
Supplementary Tables 1, 5, 6 and 8.
Additional file 6:
Supplementary Material, including Supplementary Tables 2-4, 7 and 9, as well as Supplementary Table and Figure legends.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Axelsson, G.T., Jonmundsson, T., Woo, Y. et al. Proteomic associations with forced expiratory volume: a Mendelian randomisation study. Respir Res 25, 44 (2024). https://doi.org/10.1186/s12931-023-02587-z
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12931-023-02587-z