Whole genome sequencing based analysis of inflammation biomarkers in the Trans-Omics for Precision Medicine (TOPMed) consortium.

Inflammation biomarkers can provide valuable insight into the role of inflammatory processes in many diseases and conditions. Sequencing based analyses of such biomarkers can also serve as an exemplar of the genetic architecture of quantitative traits. To evaluate the biological insight, which can be provided by a multi-ancestry, whole-genome based association study, we performed a comprehensive analysis of 21 inflammation biomarkers from up to 38 465 individuals with whole-genome sequencing from the Trans-Omics for Precision Medicine (TOPMed) program (with varying sample size by trait, where the minimum sample size was n = 737 for MMP-1). We identified 22 distinct single-variant associations across 6 traits-E-selectin, intercellular adhesion molecule 1, interleukin-6, lipoprotein-associated phospholipase A2 activity and mass, and P-selectin-that remained significant after conditioning on previously identified associations for these inflammatory biomarkers. We further expanded upon known biomarker associations by pairing the single-variant analysis with a rare variant set-based analysis that further identified 19 significant rare variant set-based associations with 5 traits. These signals were distinct from both significant single variant association signals within TOPMed and genetic signals observed in prior studies, demonstrating the complementary value of performing both single and rare variant analyses when analyzing quantitative traits. We also confirm several previously reported signals from semi-quantitative proteomics platforms. Many of these signals demonstrate the extensive allelic heterogeneity and ancestry-differentiated variant-trait associations common for inflammation biomarkers, a characteristic we hypothesize will be increasingly observed with well-powered, large-scale analyses of complex traits.

Blood-based Transcriptomic and Proteomic Biomarkers of Emphysema.

Rationale: Emphysema is a chronic obstructive pulmonary disease phenotype with important prognostic implications. Identifying blood-based biomarkers of emphysema will facilitate early diagnosis and development of targeted therapies. Objectives: To discover blood omics biomarkers for chest computed tomography-quantified emphysema and develop predictive biomarker panels. Methods: Emphysema blood biomarker discovery was performed using differential gene expression, alternative splicing, and protein association analyses in a training sample of 2,370 COPDGene participants with available blood RNA sequencing, plasma proteomics, and clinical data. Internal validation was conducted in a COPDGene testing sample (n = 1,016), and external validation was done in the ECLIPSE study (n = 526). Because low body mass index (BMI) and emphysema often co-occur, we performed a mediation analysis to quantify the effect of BMI on gene and protein associations with emphysema. Elastic net models with bootstrapping were also developed in the training sample sequentially using clinical, blood cell proportions, RNA-sequencing, and proteomic biomarkers to predict quantitative emphysema. Model accuracy was assessed by the area under the receiver operating characteristic curves for subjects stratified into tertiles of emphysema severity. Measurements and Main Results: Totals of 3,829 genes, 942 isoforms, 260 exons, and 714 proteins were significantly associated with emphysema (false discovery rate, 5%) and yielded 11 biological pathways. Seventy-four percent of these genes and 62% of these proteins showed mediation by BMI. Our prediction models demonstrated reasonable predictive performance in both COPDGene and ECLIPSE. The highest-performing model used clinical, blood cell, and protein data (area under the receiver operating characteristic curve in COPDGene testing, 0.90; 95% confidence interval, 0.85-0.90). Conclusions: Blood transcriptome and proteome-wide analyses revealed key biological pathways of emphysema and enhanced the prediction of emphysema.

Biological Age, Chronological Age and Survival in Pulmonary Fibrosis: A Causal Mediation Analysis.

RATIONALE: Accelerated biological aging has been implicated in the development of interstitial lung disease (ILD) and other diseases of aging but remains poorly understood. OBJECTIVES: To identify plasma proteins that mediate the relationship between chronological age and survival association in patients with ILD. METHODS: Causal mediation analysis was performed to identify plasma proteins that mediated the chronological age-survival relationship in an idiopathic pulmonary fibrosis (IPF) discovery cohort. Proteins mediating this relationship after adjustment for false discovery were advanced for testing in an independent ILD validation cohort and explored in a chronic obstructive pulmonary disease (COPD) cohort. A proteomic-based measure of biological age was constructed and survival analysis performed assessing the impact of biological age and peripheral blood telomere length on the chronological age-survival relationship. RESULTS: Twenty-two proteins mediated the chronological age-survival relationship after adjustment for false discovery in the IPF discovery cohort (n=874), with nineteen remaining significant mediators of this relationship in the ILD validation cohort (n=983) and one mediating this relationship in the COPD cohort. Latent transforming growth factor beta binding protein 2 and ectodysplasin A2 receptor showed the strongest mediation across cohorts. A proteomic measure of biological age completely attenuated the chronological age-survival association and better discriminated survival than chronological age. Results were robust to adjustment for peripheral blood telomere length, which did not mediate the chronological age-survival relationship. CONCLUSIONS: Molecular measures of aging completely mediate the relationship between chronological age and survival, suggesting that chronological age has no direct effect on ILD survival.

Proteomic Biomarkers of Quantitative Interstitial Abnormalities in COPDGene and CARDIA Lung Study.

Rationale: Quantitative interstitial abnormalities (QIAs) are early measures of lung injury automatically detected on chest computed tomography scans. QIAs are associated with impaired respiratory health and share features with advanced lung diseases, but their biological underpinnings are not well understood. Objectives: To identify novel protein biomarkers of QIAs using high-throughput plasma proteomic panels within two multicenter cohorts. Methods: We measured the plasma proteomics of 4,383 participants in an older, ever-smoker cohort (COPDGene [Genetic Epidemiology of Chronic Obstructive Pulmonary Disease]) and 2,925 participants in a younger population cohort (CARDIA [Coronary Artery Disease Risk in Young Adults]) using the SomaLogic SomaScan assays. We measured QIAs using a local density histogram method. We assessed the associations between proteomic biomarker concentrations and QIAs using multivariable linear regression models adjusted for age, sex, body mass index, smoking status, and study center (Benjamini-Hochberg false discovery rate-corrected P ⩽ 0.05). Measurements and Main Results: In total, 852 proteins were significantly associated with QIAs in COPDGene and 185 in CARDIA. Of the 144 proteins that overlapped between COPDGene and CARDIA, all but one shared directionalities and magnitudes. These proteins were enriched for 49 Gene Ontology pathways, including biological processes in inflammatory response, cell adhesion, immune response, ERK1/2 regulation, and signaling; cellular components in extracellular regions; and molecular functions including calcium ion and heparin binding. Conclusions: We identified the proteomic biomarkers of QIAs in an older, smoking population with a higher prevalence of pulmonary disease and in a younger, healthier community cohort. These proteomics features may be markers of early precursors of advanced lung diseases.

The Effect of Chronic Altitude Exposure on COPD Outcomes in the SPIROMICS Cohort.

RATIONALE: Individuals with COPD have airflow obstruction and maldistribution of ventilation. For those living at high altitude, any gas exchange abnormality is compounded by reduced partial pressures of inspired oxygen. OBJECTIVES: Does residence at higher-altitude exposure affect COPD outcomes, including lung function, imaging characteristics, symptoms, health status, functional exercise capacity, exacerbations, or mortality? METHODS: From the SPIROMICS cohort, we identified individuals with COPD living below 1,000 ft (305 m) elevation (n= 1,367) versus above 4,000 ft (1,219 m) elevation (n= 288). Multivariable regression models were used to evaluate associations of exposure to high altitude with COPD-related outcomes. MEASUREMENTS AND MAIN RESULTS: Living at higher altitude was associated with reduced functional exercise capacity as defined by 6MWD (-32.3 m, (-55.7 to -28.6)). There were no differences in patient-reported outcomes as defined by symptoms (CAT, mMRC), or health status (SGRQ). Higher altitude was not associated with a different rate of FEV1 decline. Higher altitude was associated with lower odds of severe exacerbations (IRR 0.65, (0.46 to 0.90)). There were no differences in small airway disease, air trapping, or emphysema. In longitudinal analyses, higher altitude was associated with increased mortality (HR 1.25, (1.0 to 1.55)); however, this association was no longer significant when accounting for air pollution. CONCLUSIONS: Chronic altitude exposure is associated with reduced functional exercise capacity in individuals with COPD, but this did not translate into differences in symptoms or health status. Additionally, chronic high-altitude exposure did not affect progression of disease as defined by longitudinal changes in spirometry.

Preacinar Arterial Dilation Mediates Outcomes of Quantitative Interstitial Abnormalities in COPDGene.

RATIONALE: Quantitative interstitial abnormalities (QIA) are a computed tomography (CT) measure of early parenchymal lung disease associated with worse clinical outcomes including exercise capacity and symptoms. The presence of pulmonary vasculopathy in QIA and its role in the QIA-outcome relationship is unknown. OBJECTIVES: To quantify radiographic pulmonary vasculopathy in quantitative interstitial abnormalities (QIA) and determine if this vasculopathy mediates the QIA-outcome relationship. METHODS: Ever-smokers with QIA, outcome, and pulmonary vascular mediator data were identified from the COPDGene cohort. CT-based vascular mediators were: right ventricle-to-left ventricle ratio (RV/LV), pulmonary artery-to-aorta ratio (PA/Ao), and pre-acinar intraparenchymal arterial dilation (PA volume 5-20mm2 in cross-sectional area, normalized to total arterial volume). Outcomes were: six-minute walk distance (6MWD) and modified Medical Council Research Council (mMRC) Dyspnea score ≥2. Adjusted causal mediation analyses were used to determine if the pulmonary vasculature mediated the QIA effect on outcomes. Associations of pre-acinar arterial dilation with select plasma biomarkers of pulmonary vascular dysfunction were examined. MAIN RESULTS: Among 8,200 participants, QIA burden correlated positively with vascular damage measures including pre-acinar arterial dilation. Pre-acinar arterial dilation mediated 79.6% of the detrimental impact of QIA on 6MWD (56.2-100%, p<0.001). PA/Ao was a weak mediator and RV/LV was a suppressor. Similar results were observed in the QIA-mMRC relationship. Pre-acinar arterial dilation correlated with increased pulmonary vascular dysfunction biomarker levels including angiopoietin-2 and NT-proBNP. CONCLUSIONS: Parenchymal quantitative interstitial abnormalities (QIA) deleteriously impact outcomes primarily through pulmonary vasculopathy. Pre-acinar arterial dilation may be a novel marker of pulmonary vasculopathy in QIA.

Association of Ground Glass Opacities with Systemic Inflammation and Progression of Emphysema.

RATIONAL: Ground glass opacities (GGO) in the absence of interstitial lung disease are understudied. OBJECTIVE: To assess the association of GGO with white blood cells (WBCs) and progression of quantified chest CT emphysema. METHODS: We analyzed data of participants in the Subpopulations and Intermediate Outcome Measures In COPD Study (SPIROMICS). Chest radiologists and pulmonologists labeled regions of the lung as GGO and adaptive multiple feature method (AMFM) trained the computer to assign those labels to image voxels and quantify the volume of the lung with GGO (%GGOAMFM). We used multivariable linear regression, zero-inflated negative binomial, and proportional hazards regression models to assess the association of %GGOAMFM with WBC, changes in %emphysema, and clinical outcomes. MEASUREMENTS AND MAIN RESULTS: Among 2,714 participants, 1,680 had COPD and 1,034 had normal spirometry. Among COPD participants, based on the multivariable analysis, current smoking and chronic productive cough was associated with higher %GGOAMFM. Higher %GGOAMFM was cross-sectionally associated with higher WBCs and neutrophils levels. Higher %GGOAMFM per interquartile range at visit 1 (baseline) was associated with an increase in emphysema at one-year follow visit by 11.7% (Relative increase; 95%CI 7.5-16.1%;P<0.001). We found no association between %GGOAMFM and one-year FEV1 decline but %GGOAMFM was associated with exacerbations and all-cause mortality during a median follow-up time of 1,544 days (Interquartile Interval=1,118-2,059). Among normal spirometry participants, we found similar results except that %GGOAMFM was associated with progression to COPD at one-year follow-up. CONCLUSIONS: Our findings suggest that GGOAMFM is associated with increased systemic inflammation and emphysema progression.

Loss of Airway Phylogenetic Diversity Is Associated with Clinical and Pathobiological Markers of Disease Development in Chronic Obstructive Pulmonary Disease.

Rationale: The airway microbiome has the potential to shape chronic obstructive pulmonary disease (COPD) pathogenesis, but its relationship to outcomes in milder disease is unestablished. Objectives: To identify sputum microbiome characteristics associated with markers of COPD in participants of the Subpopulations and Intermediate Outcome Measures of COPD Study (SPIROMICS). Methods: Sputum DNA from 877 participants was analyzed using 16S ribosomal RNA gene sequencing. Relationships between baseline airway microbiota composition and clinical, radiographic, and mucoinflammatory markers, including longitudinal lung function trajectory, were examined. Measurements and Main Results: Participant data represented predominantly milder disease (Global Initiative for Chronic Obstructive Lung Disease stage 0-2 obstruction in 732 of 877 participants). Phylogenetic diversity (i.e., range of different species within a sample) correlated positively with baseline lung function, decreased with higher Global Initiative for Chronic Obstructive Lung Disease stage, and correlated negatively with symptom burden, radiographic markers of airway disease, and total mucin concentrations (P < 0.001). In covariate-adjusted regression models, organisms robustly associated with better lung function included Alloprevotella, Oribacterium, and Veillonella species. Conversely, lower lung function, greater symptoms, and radiographic measures of small airway disease were associated with enrichment in members of Streptococcus, Actinobacillus, Actinomyces, and other genera. Baseline sputum microbiota features were also associated with lung function trajectory during SPIROMICS follow-up (stable/improved, decline, or rapid decline groups). The stable/improved group (slope of FEV1 regression ⩾66th percentile) had greater bacterial diversity at baseline associated with enrichment in Prevotella, Leptotrichia, and Neisseria species. In contrast, the rapid decline group (FEV1 slope ⩽33rd percentile) had significantly lower baseline diversity associated with enrichment in Streptococcus species. Conclusions: In SPIROMICS, baseline airway microbiota features demonstrate divergent associations with better or worse COPD-related outcomes.

The Regulation of Fatty Acid Synthase by Exosomal miR-143-5p and miR-342-5p in Idiopathic Pulmonary Fibrosis.

Idiopathic pulmonary fibrosis (IPF) is a chronic and progressive disease caused by an aberrant repair of injured alveolar epithelial cells. The maintenance of the alveolar epithelium and its regeneration after the damage is fueled by alveolar type II (ATII) cells. Injured cells release exosomes containing microRNAs (miRNAs), which can alter the recipient cells' function. Lung tissue, ATII cells, fibroblasts, plasma, and exosomes were obtained from naive patients with IPF, patients with IPF taking pirfenidone or nintedanib, and control organ donors. miRNA expression was analyzed to study their impact on exosome-mediated effects in IPF. High miR-143-5p and miR-342-5p levels were detected in ATII cells, lung tissue, plasma, and exosomes in naive patients with IPF. Decreased FASN (fatty acid synthase) and ACSL-4 (acyl-CoA-synthetase long-chain family member 4) expression was found in ATII cells. miR-143-5p and miR-342-5p overexpression or ATII cell treatment with IPF-derived exosomes containing these miRNAs lowered FASN and ACSL-4 levels. Also, this contributed to ATII cell injury and senescence. However, exosomes isolated from patients with IPF taking nintedanib or pirfenidone increased FASN expression in ATII cells compared with naive patients with IPF. Furthermore, fibroblast treatment with exosomes obtained from naive patients with IPF increased SMAD3, CTGF, COL3A1, and TGFß1 expression. Our results suggest that IPF-derived exosomes containing miR-143-5p and miR-342-5p inhibited the de novo fatty acid synthesis pathway in ATII cells. They also induced the profibrotic response in fibroblasts. Pirfenidone and nintedanib improved ATII cell function and inhibited fibrogenesis. This study highlights the importance of exosomes in IPF pathophysiology.

NetSHy: network summarization via a hybrid approach leveraging topological properties.

MOTIVATION: Biological networks can provide a system-level understanding of underlying processes. In many contexts, networks have a high degree of modularity, i.e. they consist of subsets of nodes, often known as subnetworks or modules, which are highly interconnected and may perform separate functions. In order to perform subsequent analyses to investigate the association between the identified module and a variable of interest, a module summarization, that best explains the module's information and reduces dimensionality is often needed. Conventional approaches for obtaining network representation typically rely only on the profiles of the nodes within the network while disregarding the inherent network topological information. RESULTS: In this article, we propose NetSHy, a hybrid approach which is capable of reducing the dimension of a network while incorporating topological properties to aid the interpretation of the downstream analyses. In particular, NetSHy applies principal component analysis (PCA) on a combination of the node profiles and the well-known Laplacian matrix derived directly from the network similarity matrix to extract a summarization at a subject level. Simulation scenarios based on random and empirical networks at varying network sizes and sparsity levels show that NetSHy outperforms the conventional PCA approach applied directly on node profiles, in terms of recovering the true correlation with a phenotype of interest and maintaining a higher amount of explained variation in the data when networks are relatively sparse. The robustness of NetSHy is also demonstrated by a more consistent correlation with the observed phenotype as the sample size decreases. Lastly, a genome-wide association study is performed as an application of a downstream analysis, where NetSHy summarization scores on the biological networks identify more significant single nucleotide polymorphisms than the conventional network representation. AVAILABILITY AND IMPLEMENTATION: R code implementation of NetSHy is available at https://github.com/thaovu1/NetSHy. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Causes of and Clinical Features Associated with Death in Tobacco Cigarette Users by Lung Function Impairment.

Rationale: Cigarette smoking contributes to the risk of death through different mechanisms. Objectives: To determine how causes of and clinical features associated with death vary in tobacco cigarette users by lung function impairment. Methods: We stratified current and former tobacco cigarette users enrolled in Genetic Epidemiology of Chronic Obstructive Pulmonary Disease (COPDGene) into normal spirometry, PRISm (Preserved Ratio Impaired Spirometry), Global Initiative for Chronic Obstructive Lung Disease (GOLD) 1-2 COPD, and GOLD 3-4 COPD. Deaths were identified via longitudinal follow-up and Social Security Death Index search. Causes of death were adjudicated after a review of death certificates, medical records, and next-of-kin interviews. We tested associations between baseline clinical variables and all-cause mortality using multivariable Cox proportional hazards models. Measurements and Main Results: Over a 10.1-year median follow-up, 2,200 deaths occurred among 10,132 participants (age 59.5 ± 9.0 yr; 46.6% women). Death from cardiovascular disease was most frequent in PRISm (31% of deaths). Lung cancer deaths were most frequent in GOLD 1-2 (18% of deaths vs. 9-11% in other groups). Respiratory deaths outpaced competing causes of death in GOLD 3-4, particularly when BODE index ⩾7. St. George's Respiratory Questionnaire score ⩾25 was associated with higher mortality in all groups: Hazard ratio (HR), 1.48 (1.20-1.84) normal spirometry; HR, 1.40 (1.05-1.87) PRISm; HR, 1.80 (1.49-2.17) GOLD 1-2; HR, 1.65 (1.26-2.17) GOLD 3-4. History of respiratory exacerbations was associated with higher mortality in GOLD 1-2 and GOLD 3-4, quantitative emphysema in GOLD 1-2, and airway wall thickness in PRISm and GOLD 3-4. Conclusions: Leading causes of death vary by lung function impairment in tobacco cigarette users. Worse respiratory-related quality of life is associated with all-cause mortality regardless of lung function.

Alpha-1 antitrypsin inhibits fractalkine-mediated monocyte-lung endothelial cell interactions.

Chronic obstructive pulmonary disease (COPD) is characterized by nonresolving inflammation fueled by breach in the endothelial barrier and leukocyte recruitment into the airspaces. Among the ligand-receptor axes that control leukocyte recruitment, the full-length fractalkine ligand (CX3CL1)-receptor (CX3CR1) ensures homeostatic endothelial-leukocyte interactions. Cigarette smoke (CS) exposure and respiratory pathogens increase expression of endothelial sheddases, such as a-disintegrin-and-metalloproteinase-domain 17 (ADAM17, TACE), inhibited by the anti-protease α-1 antitrypsin (AAT). In the systemic endothelium, TACE cleaves CX3CL1 to release soluble CX3CL1 (sCX3CL1). During CS exposure, it is not known whether AAT inhibits sCX3CL1 shedding and CX3CR1+ leukocyte transendothelial migration across lung microvasculature. We investigated the mechanism of sCX3CL1 shedding, its role in endothelial-monocyte interactions, and AAT effect on these interactions during acute inflammation. We used two, CS and lipopolysaccharide (LPS) models of acute inflammation in transgenic Cx3cr1gfp/gfp mice and primary human endothelial cells and monocytes to study sCX3CL1-mediated CX3CR1+ monocyte adhesion and migration. We measured sCX3CL1 levels in plasma and bronchoalveolar lavage (BALF) of individuals with COPD. Both sCX3CL1 shedding and CX3CR1+ monocytes transendothelial migration were triggered by LPS and CS exposure in mice, and were significantly attenuated by AAT. The inhibition of monocyte-endothelial adhesion and migration by AAT was TACE-dependent. Compared with healthy controls, sCX3CL1 levels were increased in plasma and BALF of individuals with COPD, and were associated with clinical parameters of emphysema. Our results indicate that inhibition of sCX3CL1 as well as AAT augmentation may be effective approaches to decrease excessive monocyte lung recruitment during acute and chronic inflammatory states.NEW & NOTEWORTHY Our novel findings that AAT and other inhibitors of TACE, the sheddase that controls full-length fractalkine (CX3CL1) endothelial expression, may provide fine-tuning of the CX3CL1-CX3CR1 axis specifically involved in endothelial-monocyte cross talk and leukocyte recruitment to the alveolar space, suggests that AAT and inhibitors of sCX3CL1 signaling may be harnessed to reduce lung inflammation.

Electronic cigarette menthol flavoring is associated with increased inhaled micro and sub-micron particles and worse lung function in combustion cigarette smokers.

Flavored electronic cigarettes (ECs) present a serious health challenge globally. Currently, it is unknown whether the addition of highly popular menthol flavoring to e-liquid is associated with changes in the number of aerosolized particles generated or altered lung function. Here, we first performed preclinical studies using our novel robotic platform Human Vaping Mimetic Real-Time Particle Analyzer (HUMITIPAA). HUMITIPAA generates fresh aerosols for any desired EC in a very controlled and user-definable manner and utilizes an optical sensing system to quantitate and analyze sub-micron and microparticles from every puff over the course of vaping session in real-time while emulating clinically relevant breathing mechanics and vaping topography. We discovered that addition of menthol flavoring to freshly prepared e-liquid base propylene glycol-vegetable glycerin leads to enhanced particle counts in all tested size fractions, similar to the effect of adding vitamin E acetate to e-liquid we previously reported. Similarly, we found that menthol vs. non-menthol (tobacco) flavored pods from commercially available ECs leads to generation of significantly higher quantities of 1-10 µm particles upon inhalation. We then retrospectively analyzed data from the COPDGene study and identified an association between the use of menthol flavored ECs and reduced FEV1% predicted and FEV1/FVC independent of age, gender, race, pack-years of smoking, and use of nicotine or cannabis-containing vaping products. Our results reveal an association between enhanced inhaled particle due to menthol addition to ECs and worse lung function indices. Detailed causal relation remains to be demonstrated in future large-scale prospective clinical studies. Importantly, here we demonstrate utility of the HUMITIPAA as a predictive enabling technology to identify inhalation toxicological potential of emerging ECs as the chemical formulation of e-liquid gets modified.

Plasma metabolomics and quantitative interstitial abnormalities in ever-smokers.

BACKGROUND: Quantitative interstitial abnormalities (QIA) are an automated computed tomography (CT) finding of early parenchymal lung disease, associated with worse lung function, reduced exercise capacity, increased respiratory symptoms, and death. The metabolomic perturbations associated with QIA are not well known. We sought to identify plasma metabolites associated with QIA in smokers. We also sought to identify shared and differentiating metabolomics features between QIA and emphysema, another smoking-related advanced radiographic abnormality. METHODS: In 928 former and current smokers in the Genetic Epidemiology of COPD cohort, we measured QIA and emphysema using an automated local density histogram method and generated metabolite profiles from plasma samples using liquid chromatography-mass spectrometry (Metabolon). We assessed the associations between metabolite levels and QIA using multivariable linear regression models adjusted for age, sex, body mass index, smoking status, pack-years, and inhaled corticosteroid use, at a Benjamini-Hochberg False Discovery Rate p-value of ≤ 0.05. Using multinomial regression models adjusted for these covariates, we assessed the associations between metabolite levels and the following CT phenotypes: QIA-predominant, emphysema-predominant, combined-predominant, and neither- predominant. Pathway enrichment analyses were performed using MetaboAnalyst. RESULTS: We found 85 metabolites significantly associated with QIA, with overrepresentation of the nicotinate and nicotinamide, histidine, starch and sucrose, pyrimidine, phosphatidylcholine, lysophospholipid, and sphingomyelin pathways. These included metabolites involved in inflammation and immune response, extracellular matrix remodeling, surfactant, and muscle cachexia. There were 75 metabolites significantly different between QIA-predominant and emphysema-predominant phenotypes, with overrepresentation of the phosphatidylethanolamine, nicotinate and nicotinamide, aminoacyl-tRNA, arginine, proline, alanine, aspartate, and glutamate pathways. CONCLUSIONS: Metabolomic correlates may lend insight to the biologic perturbations and pathways that underlie clinically meaningful quantitative CT measurements like QIA in smokers.

Levels of Urinary Mercapturic Acids of Acrolein, Methacrolein, Crotonaldehyde, and Methyl Vinyl Ketone in Relationship to Chronic Obstructive Pulmonary Disease in Cigarette Smokers of the Subpopulations and Intermediate Outcome Measures in COPD Study (SPIROMICS).

Cigarette smoking is an established cause of chronic obstructive pulmonary disease (COPD). Numerous studies implicate acrolein, which occurs in relatively high concentrations in cigarette smoke and reacts readily with proteins, as one causative factor for COPD in smokers. Far less is known about the possible roles in COPD of the related α,ß-unsaturated carbonyl compounds of cigarette smoke crotonaldehyde, methacrolein, and methyl vinyl ketone. In the study reported here, we analyzed mercapturic acids of these α,ß-unsaturated compounds in the urine of 413 confirmed cigarette smokers in the Subpopulations and Intermediate Outcome Measures in COPD Study (SPIROMICS)─202 with COPD and 211 without COPD. The mercapturic acids analyzed were 3-hydroxypropyl mercapturic acid (3-HPMA) from acrolein, 3-hydroxy-1-methylpropyl mercapturic acid (HMPMA-1) from crotonaldehyde, 3-hydroxy-2-methylpropyl mercapturic acid (HMPMA-2) from methacrolein, and 3-hydroxy-3-methylpropyl mercapturic acid (HMPMA-3) from methyl vinyl ketone. In models adjusting for age, sex, race, pack years of tobacco use, and BMI, all four mercapturic acids were increased in individuals with COPD but not significantly. Stratified by the GOLD status, there were increased levels of the metabolites associated with GOLD 3-4 compared to that with GOLD 0, with the methacrolein metabolite HMPMA-2 reaching statistical significance (adjusted odds ratio 1.23 [95% CI: 1.00-1.53]). These results highlight the possible role of methacrolein, which has previously received little attention in this regard, as a causative factor in COPD in cigarette smokers.

The Proteomic Profile of Interstitial Lung Abnormalities.

Rationale: Knowledge on biomarkers of interstitial lung disease is incomplete. Interstitial lung abnormalities (ILAs) are radiologic changes that may present in its early stages. Objectives: To uncover blood proteins associated with ILAs using large-scale proteomics methods. Methods: Data from two prospective cohort studies, the AGES-Reykjavik (Age, Gene/Environment Susceptibility-Reykjavik) study (N = 5,259) for biomarker discovery and the COPDGene (Genetic Epidemiology of COPD) study (N = 4,899) for replication, were used. Blood proteins were measured using DNA aptamers, targeting more than 4,700 protein analytes. The association of proteins with ILAs and ILA progression was assessed with regression modeling, as were associations with genetic risk factors. Adaptive Least Absolute Shrinkage and Selection Operator models were applied to bootstrap data samples to discover sets of proteins predictive of ILAs and their progression. Measurements and Main Results: Of 287 associations, SFTPB (surfactant protein B) (odds ratio [OR], 3.71 [95% confidence interval (CI), 3.20-4.30]; P = 4.28 × 10-67), SCGB3A1 (Secretoglobin family 3A member 1) (OR, 2.43 [95% CI, 2.13-2.77]; P = 8.01 × 10-40), and WFDC2 (WAP four-disulfide core domain protein 2) (OR, 2.42 [95% CI, 2.11-2.78]; P = 4.01 × 10-36) were most significantly associated with ILA in AGES-Reykjavik and were replicated in COPDGene. In AGES-Reykjavik, concentrations of SFTPB were associated with the rs35705950 MUC5B (mucin 5B) promoter polymorphism, and SFTPB and WFDC2 had the strongest associations with ILA progression. Multivariate models of ILAs in AGES-Reykjavik, ILAs in COPDGene, and ILA progression in AGES-Reykjavik had validated areas under the receiver operating characteristic curve of 0.880, 0.826, and 0.824, respectively. Conclusions: Novel, replicated associations of ILA, its progression, and genetic risk factors with numerous blood proteins are demonstrated as well as machine-learning-based models with favorable predictive potential. Several proteins are revealed as potential markers of early fibrotic lung disease.

Lung Microbiota and Metabolites Collectively Associate with Clinical Outcomes in Milder Stage Chronic Obstructive Pulmonary Disease.

Rationale: Chronic obstructive pulmonary disease (COPD) is variable in its development. Lung microbiota and metabolites collectively may impact COPD pathophysiology, but relationships to clinical outcomes in milder disease are unclear. Objectives: Identify components of the lung microbiome and metabolome collectively associated with clinical markers in milder stage COPD. Methods: We analyzed paired microbiome and metabolomic data previously characterized from bronchoalveolar lavage fluid in 137 participants in the SPIROMICS (Subpopulations and Intermediate Outcome Measures in COPD Study), or (GOLD [Global Initiative for Chronic Obstructive Lung Disease Stage 0-2). Datasets used included 1) bacterial 16S rRNA gene sequencing; 2) untargeted metabolomics of the hydrophobic fraction, largely comprising lipids; and 3) targeted metabolomics for a panel of hydrophilic compounds previously implicated in mucoinflammation. We applied an integrative approach to select features and model 14 individual clinical variables representative of known associations with COPD trajectory (lung function, symptoms, and exacerbations). Measurements and Main Results: The majority of clinical measures associated with the lung microbiome and metabolome collectively in overall models (classification accuracies, >50%, P < 0.05 vs. chance). Lower lung function, COPD diagnosis, and greater symptoms associated positively with Streptococcus, Neisseria, and Veillonella, together with compounds from several classes (glycosphingolipids, glycerophospholipids, polyamines and xanthine, an adenosine metabolite). In contrast, several Prevotella members, together with adenosine, 5'-methylthioadenosine, sialic acid, tyrosine, and glutathione, associated with better lung function, absence of COPD, or less symptoms. Significant correlations were observed between specific metabolites and bacteria (Padj < 0.05). Conclusions: Components of the lung microbiome and metabolome in combination relate to outcome measures in milder COPD, highlighting their potential collaborative roles in disease pathogenesis.

Longitudinal Follow-Up of Participants With Tobacco Exposure and Preserved Spirometry.

Importance: People who smoked cigarettes may experience respiratory symptoms without spirometric airflow obstruction. These individuals are typically excluded from chronic obstructive pulmonary disease (COPD) trials and lack evidence-based therapies. Objective: To define the natural history of persons with tobacco exposure and preserved spirometry (TEPS) and symptoms (symptomatic TEPS). Design, Setting, and Participants: SPIROMICS II was an extension of SPIROMICS I, a multicenter study of persons aged 40 to 80 years who smoked cigarettes (>20 pack-years) with or without COPD and controls without tobacco exposure or airflow obstruction. Participants were enrolled in SPIROMICS I and II from November 10, 2010, through July 31, 2015, and followed up through July 31, 2021. Exposures: Participants in SPIROMICS I underwent spirometry, 6-minute walk distance testing, assessment of respiratory symptoms, and computed tomography of the chest at yearly visits for 3 to 4 years. Participants in SPIROMICS II had 1 additional in-person visit 5 to 7 years after enrollment in SPIROMICS I. Respiratory symptoms were assessed with the COPD Assessment Test (range, 0 to 40; higher scores indicate more severe symptoms). Participants with symptomatic TEPS had normal spirometry (postbronchodilator ratio of forced expiratory volume in the first second [FEV1] to forced vital capacity >0.70) and COPD Assessment Test scores of 10 or greater. Participants with asymptomatic TEPS had normal spirometry and COPD Assessment Test scores of less than 10. Patient-reported respiratory symptoms and exacerbations were assessed every 4 months via phone calls. Main Outcomes and Measures: The primary outcome was assessment for accelerated decline in lung function (FEV1) in participants with symptomatic TEPS vs asymptomatic TEPS. Secondary outcomes included development of COPD defined by spirometry, respiratory symptoms, rates of respiratory exacerbations, and progression of computed tomographic-defined airway wall thickening or emphysema. Results: Of 1397 study participants, 226 had symptomatic TEPS (mean age, 60.1 [SD, 9.8] years; 134 were women [59%]) and 269 had asymptomatic TEPS (mean age, 63.1 [SD, 9.1] years; 134 were women [50%]). At a median follow-up of 5.76 years, the decline in FEV1 was -31.3 mL/y for participants with symptomatic TEPS vs -38.8 mL/y for those with asymptomatic TEPS (between-group difference, -7.5 mL/y [95% CI, -16.6 to 1.6 mL/y]). The cumulative incidence of COPD was 33.0% among participants with symptomatic TEPS vs 31.6% among those with asymptomatic TEPS (hazard ratio, 1.05 [95% CI, 0.76 to 1.46]). Participants with symptomatic TEPS had significantly more respiratory exacerbations than those with asymptomatic TEPS (0.23 vs 0.08 exacerbations per person-year, respectively; rate ratio, 2.38 [95% CI, 1.71 to 3.31], P < .001). Conclusions and Relevance: Participants with symptomatic TEPS did not have accelerated rates of decline in FEV1 or increased incidence of COPD vs those with asymptomatic TEPS, but participants with symptomatic TEPS did experience significantly more respiratory exacerbations over a median follow-up of 5.8 years.

Efficient Variant Set Mixed Model Association Tests for Continuous and Binary Traits in Large-Scale Whole-Genome Sequencing Studies.

With advances in whole-genome sequencing (WGS) technology, more advanced statistical methods for testing genetic association with rare variants are being developed. Methods in which variants are grouped for analysis are also known as variant-set, gene-based, and aggregate unit tests. The burden test and sequence kernel association test (SKAT) are two widely used variant-set tests, which were originally developed for samples of unrelated individuals and later have been extended to family data with known pedigree structures. However, computationally efficient and powerful variant-set tests are needed to make analyses tractable in large-scale WGS studies with complex study samples. In this paper, we propose the variant-set mixed model association tests (SMMAT) for continuous and binary traits using the generalized linear mixed model framework. These tests can be applied to large-scale WGS studies involving samples with population structure and relatedness, such as in the National Heart, Lung, and Blood Institute's Trans-Omics for Precision Medicine (TOPMed) program. SMMATs share the same null model for different variant sets, and a virtue of this null model, which includes covariates only, is that it needs to be fit only once for all tests in each genome-wide analysis. Simulation studies show that all the proposed SMMATs correctly control type I error rates for both continuous and binary traits in the presence of population structure and relatedness. We also illustrate our tests in a real data example of analysis of plasma fibrinogen levels in the TOPMed program (n = 23,763), using the Analysis Commons, a cloud-based computing platform.

Plasma sRAGE levels strongly associate with centrilobular emphysema assessed by HRCT scans.

BACKGROUND: There is a strong need for biomarkers to better characterize individuals with COPD and to take into account the heterogeneity of COPD. The blood protein sRAGE has been put forward as promising biomarker for COPD in general and emphysema in particular. Here, we measured plasma sRAGE levels using quantitative LC-MS and assessed whether the plasma sRAGE levels associate with (changes in) lung function, radiological emphysema parameters, and radiological subtypes of emphysema. METHODS: Three hundred and twenty-four COPD patients (mean FEV1: 63%predicted) and 185 healthy controls from the COPDGene study were selected. Plasma sRAGE was measured by immunoprecipitation in 96-well plate methodology to enrich sRAGE, followed by targeted quantitative liquid chromatography-mass spectrometry. Spirometry and HRCT scans (inspiration and expiration) with a 5-year follow-up were used; both subjected to high quality control standards. RESULTS: Lower sRAGE values significantly associated with the presence of COPD, the severity of airflow obstruction, the severity of emphysema on HRCT, the heterogeneous distribution of emphysema, centrilobular emphysema, and 5-year progression of emphysema. However, sRAGE values did not associate with airway wall thickness or paraseptal emphysema. CONCLUSIONS: Rather than being a general COPD biomarker, sRAGE is especially a promising biomarker for centrilobular emphysema. Follow-up studies should elucidate whether sRAGE can be used as a biomarker for other COPD phenotypes as well.