Molecular Transducers of Physical Activity Consortium (MoTrPAC): Mapping the Dynamic Responses to Exercise.

Exercise provides a robust physiological stimulus that evokes cross-talk among multiple tissues that when repeated regularly (i.e., training) improves physiological capacity, benefits numerous organ systems, and decreases the risk for premature mortality. However, a gap remains in identifying the detailed molecular signals induced by exercise that benefits health and prevents disease. The Molecular Transducers of Physical Activity Consortium (MoTrPAC) was established to address this gap and generate a molecular map of exercise. Preclinical and clinical studies will examine the systemic effects of endurance and resistance exercise across a range of ages and fitness levels by molecular probing of multiple tissues before and after acute and chronic exercise. From this multi-omic and bioinformatic analysis, a molecular map of exercise will be established. Altogether, MoTrPAC will provide a public database that is expected to enhance our understanding of the health benefits of exercise and to provide insight into how physical activity mitigates disease.

Non-Abelian topological order and anyons on a trapped-ion processor.

Non-Abelian topological order is a coveted state of matter with remarkable properties, including quasiparticles that can remember the sequence in which they are exchanged1-4. These anyonic excitations are promising building blocks of fault-tolerant quantum computers5,6. However, despite extensive efforts, non-Abelian topological order and its excitations have remained elusive, unlike the simpler quasiparticles or defects in Abelian topological order. Here we present the realization of non-Abelian topological order in the wavefunction prepared in a quantum processor and demonstrate control of its anyons. Using an adaptive circuit on Quantinuum's H2 trapped-ion quantum processor, we create the ground-state wavefunction of D4 topological order on a kagome lattice of 27 qubits, with fidelity per site exceeding 98.4 per cent. By creating and moving anyons along Borromean rings in spacetime, anyon interferometry detects an intrinsically non-Abelian braiding process. Furthermore, tunnelling non-Abelions around a torus creates all 22 ground states, as well as an excited state with a single anyon-a peculiar feature of non-Abelian topological order. This work illustrates the counterintuitive nature of non-Abelions and enables their study in quantum devices.

Interaction molecular QTL mapping discovers cellular and environmental modifiers of genetic regulatory effects.

Bulk-tissue molecular quantitative trait loci (QTLs) have been the starting point for interpreting disease-associated variants, and context-specific QTLs show particular relevance for disease. Here, we present the results of mapping interaction QTLs (iQTLs) for cell type, age, and other phenotypic variables in multi-omic, longitudinal data from the blood of individuals of diverse ancestries. By modeling the interaction between genotype and estimated cell-type proportions, we demonstrate that cell-type iQTLs could be considered as proxies for cell-type-specific QTL effects, particularly for the most abundant cell type in the tissue. The interpretation of age iQTLs, however, warrants caution because the moderation effect of age on the genotype and molecular phenotype association could be mediated by changes in cell-type composition. Finally, we show that cell-type iQTLs contribute to cell-type-specific enrichment of diseases that, in combination with additional functional data, could guide future functional studies. Overall, this study highlights the use of iQTLs to gain insights into the context specificity of regulatory effects.

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.

Rare EIF4A2 variants are associated with a neurodevelopmental disorder characterized by intellectual disability, hypotonia, and epilepsy.

Eukaryotic initiation factor-4A2 (EIF4A2) is an ATP-dependent RNA helicase and a member of the DEAD-box protein family that recognizes the 5' cap structure of mRNAs, allows mRNA to bind to the ribosome, and plays an important role in microRNA-regulated gene repression. Here, we report on 15 individuals from 14 families presenting with global developmental delay, intellectual disability, hypotonia, epilepsy, and structural brain anomalies, all of whom have extremely rare de novo mono-allelic or inherited bi-allelic variants in EIF4A2. Neurodegeneration was predominantly reported in individuals with bi-allelic variants. Molecular modeling predicts these variants would perturb structural interactions in key protein domains. To determine the pathogenicity of the EIF4A2 variants in vivo, we examined the mono-allelic variants in Drosophila melanogaster (fruit fly) and identified variant-specific behavioral and developmental defects. The fruit fly homolog of EIF4A2 is eIF4A, a negative regulator of decapentaplegic (dpp) signaling that regulates embryo patterning, eye and wing morphogenesis, and stem cell identity determination. Our loss-of-function (LOF) rescue assay demonstrated a pupal lethality phenotype induced by loss of eIF4A, which was fully rescued with human EIF4A2 wild-type (WT) cDNA expression. In comparison, the EIF4A2 variant cDNAs failed or incompletely rescued the lethality. Overall, our findings reveal that EIF4A2 variants cause a genetic neurodevelopmental syndrome with both LOF and gain of function as underlying mechanisms.

Whole-genome analysis of plasma fibrinogen reveals population-differentiated genetic regulators with putative liver roles.

Genetic studies have identified numerous regions associated with plasma fibrinogen levels in Europeans, yet missing heritability and limited inclusion of non-Europeans necessitates further studies with improved power and sensitivity. Compared with array-based genotyping, whole genome sequencing (WGS) data provides better coverage of the genome and better representation of non-European variants. To better understand the genetic landscape regulating plasma fibrinogen levels, we meta-analyzed WGS data from the NHLBI's Trans-Omics for Precision Medicine (TOPMed) program (n=32,572), with array-based genotype data from the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium (n=131,340) imputed to the TOPMed or Haplotype Reference Consortium panel. We identified 18 loci that have not been identified in prior genetic studies of fibrinogen. Of these, four are driven by common variants of small effect with reported MAF at least 10 percentage points higher in African populations. Three signals (SERPINA1, ZFP36L2, and TLR10) contain predicted deleterious missense variants. Two loci, SOCS3 and HPN, each harbor two conditionally distinct, non-coding variants. The gene region encoding the fibrinogen protein chain subunits (FGG;FGB;FGA), contains 7 distinct signals, including one novel signal driven by rs28577061, a variant common in African ancestry populations but extremely rare in Europeans (MAFAFR=0.180; MAFEUR=0.008). Through phenome-wide association studies in the VA Million Veteran Program, we found associations between fibrinogen polygenic risk scores and thrombotic and inflammatory disease phenotypes, including an association with gout. Our findings demonstrate the utility of WGS to augment genetic discovery in diverse populations and offer new insights for putative mechanisms of fibrinogen regulation.

Hepatic NADH reductive stress underlies common variation in metabolic traits.

The cellular NADH/NAD+ ratio is fundamental to biochemistry, but the extent to which it reflects versus drives metabolic physiology in vivo is poorly understood. Here we report the in vivo application of Lactobacillus brevis (Lb)NOX1, a bacterial water-forming NADH oxidase, to assess the metabolic consequences of directly lowering the hepatic cytosolic NADH/NAD+ ratio in mice. By combining this genetic tool with metabolomics, we identify circulating α-hydroxybutyrate levels as a robust marker of an elevated hepatic cytosolic NADH/NAD+ ratio, also known as reductive stress. In humans, elevations in circulating α-hydroxybutyrate levels have previously been associated with impaired glucose tolerance2, insulin resistance3 and mitochondrial disease4, and are associated with a common genetic variant in GCKR5, which has previously been associated with many seemingly disparate metabolic traits. Using LbNOX, we demonstrate that NADH reductive stress mediates the effects of GCKR variation on many metabolic traits, including circulating triglyceride levels, glucose tolerance and FGF21 levels. Our work identifies an elevated hepatic NADH/NAD+ ratio as a latent metabolic parameter that is shaped by human genetic variation and contributes causally to key metabolic traits and diseases. Moreover, it underscores the utility of genetic tools such as LbNOX to empower studies of 'causal metabolism'.

Canonical correlation analysis for multi-omics: Application to cross-cohort analysis.

Integrative approaches that simultaneously model multi-omics data have gained increasing popularity because they provide holistic system biology views of multiple or all components in a biological system of interest. Canonical correlation analysis (CCA) is a correlation-based integrative method designed to extract latent features shared between multiple assays by finding the linear combinations of features-referred to as canonical variables (CVs)-within each assay that achieve maximal across-assay correlation. Although widely acknowledged as a powerful approach for multi-omics data, CCA has not been systematically applied to multi-omics data in large cohort studies, which has only recently become available. Here, we adapted sparse multiple CCA (SMCCA), a widely-used derivative of CCA, to proteomics and methylomics data from the Multi-Ethnic Study of Atherosclerosis (MESA) and Jackson Heart Study (JHS). To tackle challenges encountered when applying SMCCA to MESA and JHS, our adaptations include the incorporation of the Gram-Schmidt (GS) algorithm with SMCCA to improve orthogonality among CVs, and the development of Sparse Supervised Multiple CCA (SSMCCA) to allow supervised integration analysis for more than two assays. Effective application of SMCCA to the two real datasets reveals important findings. Applying our SMCCA-GS to MESA and JHS, we identified strong associations between blood cell counts and protein abundance, suggesting that adjustment of blood cell composition should be considered in protein-based association studies. Importantly, CVs obtained from two independent cohorts also demonstrate transferability across the cohorts. For example, proteomic CVs learned from JHS, when transferred to MESA, explain similar amounts of blood cell count phenotypic variance in MESA, explaining 39.0% ~ 50.0% variation in JHS and 38.9% ~ 49.1% in MESA. Similar transferability was observed for other omics-CV-trait pairs. This suggests that biologically meaningful and cohort-agnostic variation is captured by CVs. We anticipate that applying our SMCCA-GS and SSMCCA on various cohorts would help identify cohort-agnostic biologically meaningful relationships between multi-omics data and phenotypic traits.

Whole-exome sequencing study identifies four novel gene loci associated with diabetic kidney disease.

Diabetic kidney disease (DKD) is recognized as an important public health challenge. However, its genomic mechanisms are poorly understood. To identify rare variants for DKD, we conducted a whole-exome sequencing (WES) study leveraging large cohorts well-phenotyped for chronic kidney disease and diabetes. Our two-stage WES study included 4372 European and African ancestry participants from the Chronic Renal Insufficiency Cohort and Atherosclerosis Risk in Communities studies (stage 1) and 11 487 multi-ancestry Trans-Omics for Precision Medicine participants (stage 2). Generalized linear mixed models, which accounted for genetic relatedness and adjusted for age, sex and ancestry, were used to test associations between single variants and DKD. Gene-based aggregate rare variant analyses were conducted using an optimized sequence kernel association test implemented within our mixed model framework. We identified four novel exome-wide significant DKD-related loci through initiating diabetes. In single-variant analyses, participants carrying a rare, in-frame insertion in the DIS3L2 gene (rs141560952) exhibited a 193-fold increased odds [95% confidence interval (CI): 33.6, 1105] of DKD compared with noncarriers (P = 3.59 × 10-9). Likewise, each copy of a low-frequency KRT6B splice-site variant (rs425827) conferred a 5.31-fold higher odds (95% CI: 3.06, 9.21) of DKD (P = 2.72 × 10-9). Aggregate gene-based analyses further identified ERAP2 (P = 4.03 × 10-8) and NPEPPS (P = 1.51 × 10-7), which are both expressed in the kidney and implicated in renin-angiotensin-aldosterone system modulated immune response. In the largest WES study of DKD, we identified novel rare variant loci attaining exome-wide significance. These findings provide new insights into the molecular mechanisms underlying DKD.

Dephospho-CoA kinase, a nuclear-encoded apicoplast protein, remains active and essential after Plasmodium falciparum apicoplast disruption.

Malaria parasites contain an essential organelle called the apicoplast that houses metabolic pathways for fatty acid, heme, isoprenoid, and iron-sulfur cluster synthesis. Surprisingly, malaria parasites can survive without the apicoplast as long as the isoprenoid precursor isopentenyl pyrophosphate (IPP) is supplemented in the growth medium, making it appear that isoprenoid synthesis is the only essential function of the organelle in blood-stage parasites. In the work described here, we localized an enzyme responsible for coenzyme A synthesis, DPCK, to the apicoplast, but we were unable to delete DPCK, even in the presence of IPP. However, once the endogenous DPCK was complemented with the E. coli DPCK (EcDPCK), we were successful in deleting it. We were then able to show that DPCK activity is required for parasite survival through knockdown of the complemented EcDPCK. Additionally, we showed that DPCK enzyme activity remains functional and essential within the vesicles present after apicoplast disruption. These results demonstrate that while the apicoplast of blood-stage P. falciparum parasites can be disrupted, the resulting vesicles remain biochemically active and are capable of fulfilling essential functions.

Guanylate Kinase 1 Deficiency: A Novel and Potentially Treatable Mitochondrial DNA Depletion/Deletions Disease.

OBJECTIVE: Mitochondrial DNA (mtDNA) depletion/deletions syndrome (MDDS) comprises a group of diseases caused by primary autosomal defects of mtDNA maintenance. Our objective was to study the etiology of MDDS in 4 patients who lack pathogenic variants in known genetic causes. METHODS: Whole exome sequencing of the probands was performed to identify pathogenic variants. We validated the mitochondrial defect by analyzing mtDNA, mitochondrial dNTP pools, respiratory chain activities, and GUK1 activity. To confirm pathogenicity of GUK1 deficiency, we expressed 2 GUK1 isoforms in patient cells. RESULTS: We identified biallelic GUK1 pathogenic variants in all 4 probands who presented with ptosis, ophthalmoparesis, and myopathic proximal limb weakness, as well as variable hepatopathy and altered T-lymphocyte profiles. Muscle biopsies from all probands showed mtDNA depletion, deletions, or both, as well as reduced activities of mitochondrial respiratory chain enzymes. GUK1 encodes guanylate kinase, originally identified as a cytosolic enzyme. Long and short isoforms of GUK1 exist. We observed that the long isoform is intramitochondrial and the short is cytosolic. In probands' fibroblasts, we noted decreased GUK1 activity causing unbalanced mitochondrial dNTP pools and mtDNA depletion in both replicating and quiescent fibroblasts indicating that GUK1 deficiency impairs de novo and salvage nucleotide pathways. Proband fibroblasts treated with deoxyguanosine and/or forodesine, a purine phosphatase inhibitor, ameliorated mtDNA depletion, indicating potential pharmacological therapies. INTERPRETATION: Primary GUK1 deficiency is a new and potentially treatable cause of MDDS. The cytosolic isoform of GUK1 may contribute to the T-lymphocyte abnormality, which has not been observed in other MDDS disorders. ANN NEUROL 2024.

A genetically encoded tool to increase cellular NADH/NAD+ ratio in living cells.

Impaired redox metabolism is a key contributor to the etiology of many diseases, including primary mitochondrial disorders, cancer, neurodegeneration and aging. However, mechanistic studies of redox imbalance remain challenging due to limited strategies that can perturb redox metabolism in various cellular or organismal backgrounds. Most studies involving impaired redox metabolism have focused on oxidative stress; consequently, less is known about the settings where there is an overabundance of NADH reducing equivalents, termed reductive stress. Here we introduce a soluble transhydrogenase from Escherichia coli (EcSTH) as a novel genetically encoded tool to promote reductive stress in living cells. When expressed in mammalian cells, EcSTH, and a mitochondrially targeted version (mitoEcSTH), robustly elevated the NADH/NAD+ ratio in a compartment-specific manner. Using this tool, we determined that metabolic and transcriptomic signatures of the NADH reductive stress are cellular background specific. Collectively, our novel genetically encoded tool represents an orthogonal strategy to promote reductive stress.

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.

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.

Proteomic Risk Score of Increased Respiratory Susceptibility: A Multi-Cohort Study.

RATIONALE: Accelerated decline in lung function is associated with incident COPD, hospitalizations and death. However, identifying this trajectory with longitudinal spirometry measurements is challenging in clinical practice. OBJECTIVE: To determine whether a proteomic risk score trained on accelerated decline in lung function can assess risk of future respiratory disease and mortality. METHODS: In CARDIA, a population-based cohort starting in young adulthood, longitudinal measurements of FEV1 percent predicted (up to six timepoints over 30 years) were used to identify accelerated and normal decline trajectories. Protein aptamers associated with an accelerated decline trajectory were identified with multivariable logistic regression followed by LASSO regression. The proteomic respiratory susceptibility score was derived based on these circulating proteins and applied to the UK Biobank and COPDGene studies to examine associations with future respiratory morbidity and mortality. MEASUREMENTS AND RESULTS: Higher susceptibility score was independently associated with all-cause mortality (UKBB: HR 1.56, 95%CI 1.50-1.61; COPDGene: HR 1.75, 95%CI 1.63-1.88), respiratory mortality (UKBB: HR 2.39, 95% CI 2.16-2.64; COPDGene: HR 1.83, 95%CI 1.33-2.51), incident COPD (UKBB: HR 1.84, 95%CI 1.71-1.98), incident respiratory exacerbation (COPDGene: OR 1.11, 95%CI 1.03-1.20), and incident exacerbation requiring hospitalization (COPDGene: OR 1.18, 95%CI 1.08-1.28). CONCLUSIONS: A proteomic signature of increased respiratory susceptibility identifies people at risk of respiratory death, incident COPD, and respiratory exacerbations. This susceptibility score is comprised of proteins with well-known and novel associations with lung health and holds promise for the early detection of lung disease without requiring years of spirometry measurements.

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 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 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. Measurements and Main Results: Twenty-two proteins mediated the chronological age-survival relationship after adjustment for false discovery in the idiopathic pulmonary fibrosis discovery cohort (n = 874), with 19 remaining significant mediators of this relationship in the ILD validation cohort (n = 983) and one mediating this relationship in the chronic obstructive pulmonary disease cohort. Latent transforming growth factor-ß 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.

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.

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.

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.

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.