Grünig, Gabriele

Background The current classification of pulmonary hypertension (PH), based largely on expert opinion, has limitations in prognostication and guiding therapies. We hypothesize that novel PH clusters that predict survival will reveal mechanistic phenotypes associated with biomarkers of vascular health across all PH groups. Methods We first identify novel PH clinical clusters by performing unsupervised clustering analysis on the CC‐PH (Cleveland Clinic PH) registry (N=1529). We develop classification models to predict the new PH clusters and then apply them to the multicenter PVDOMICS (Pulmonary Vascular Diseases Phenomics) cohort (N=853) for validation. We compare transplantation‐free survival across the new PH clusters. We quantify metabolites of the arginine‐nitric oxide pathway and D‐dimer levels and calculate global arginine bioavailability (arginine/[ornithine+citrulline]) to assess endothelial function and activation in the new clusters and link these biomarkers to clinical outcomes. Results Clustering analysis identify 3 clear clusters in CC ‐ PH that are validated in PVDOMICS and outperform conventional classifications in predicting transplantation‐free survival. The phenotype associated with the worst survival is characterized by reduced lung diffusion capacity, decreased arginine bioavailability and nitrate levels, and elevated D‐dimer levels, consistent with loss of pulmonary microcirculation and endothelial dysfunction. Conclusions We identify new informative PH phenotypes associated with mortality and defined by biomarkers of endothelial function and activation. Loss of endothelial health and pronounced pulmonary vascular rarefication contribute more substantially to mortality across the spectrum of PH than right heart function. REGISTRATION: URL: https://www.clinicaltrials.gov ; Unique identifier: NCT02980887.

ABSTRACT Background Pulmonary arterial hypertension (PAH) is characterized by circulating metabolic alterations, but whether these reflect disease-specific metabolic programs or reorganization of normal metabolic architecture, and how they relate to right ventricular-pulmonary vascular function (RV-PV), remains unclear. We hypothesized that the PAH metabolome is organized into biologically coherent, co-regulated metabolic modules whose relationships to RV-PV function would provide insight into known and novel metabolic pathways. Methods We applied weighted gene co-expression network analysis (WGCNA) to untargeted metabolomic data from 412 PAH patients enrolled in the multicenter PVDOMICS study. Module preservation analysis was performed in 85 healthy controls, with external replication in an independent single-center pulmonary hypertension cohort of 89 patients. Results WGCNA identified 16 distinct metabolic modules organized around biologically coherent programs. A coherent fatty acid axis, spanning substrate pools, β-oxidation intermediates, and conjugated fatty acid disposal products, formed a central organizing structure, with downstream fatty acid oxidation modules strongly associated with adverse hemodynamics and worse RV-pulmonary artery (PA) coupling. Acylcholine-enriched and 5α-reduced androgen metabolite modules were associated with favorable hemodynamic indices. Module architecture was largely preserved in healthy controls, with subtle disease-associated modular reorganization, rather than emergence of novel modules, observed in PAH. Core modules were recovered in the replication cohort with conserved hub metabolites. Conclusions These findings establish a systems-level framework demonstrating that PAH involves structured intensification and reorganization of interconnected metabolic programs associated with favorable and adverse RV-PV phenotypes. This work provides new insight into the metabolic architecture underlying PAH and identifies coordinated metabolic pathways linked to pulmonary vascular and right ventricular function.

Abstract Introduction Interstitial lung disease-associated pulmonary hypertension (ILD-PH) is a clinically challenging entity that lies on a pathophysiologic spectrum between pulmonary arterial hypertension (PAH) and interstitial lung disease (ILD). We compared the metabolomes of patients with ILD-PH to those with PAH and ILD to characterize the PAH-ILD axis and determine whether ILD-PH patients can be phenotyped based on their metabolomic similarity to PAH. Methods Using untargeted metabolomics profiles from patients with PAH, ILD, and ILD-PH in the large, U.S.-based, NIH-funded, multi-center PVDOMICS cohort, we trained an anchor-based classifier on PAH and ILD metabolomes using supervised PLS-DA with a 75/25 train/test split. Primary and secondary PH or risk-group classifications in PVDOMICS were based on expert adjudication of clinical data including right heart catheterization hemodynamics, pulmonary function testing, chest imaging, and laboratory results. We then applied the trained model to each ILD-PH patient to assign a PAH similarity score. We classified patients as having a PAH-like or ILD-like metabolomic phenotype using an optimal PAH similarity score cutoff of 0.5 then compared PAH metabolomic similarity to expert-adjudicated primary and secondary WSPH classifications in PVDOMICS. We implemented FELLA, a network-based metabolomics enrichment method, to holistically characterize the biologic underpinnings of PAH, ILD-PH, and ILD metabolomes. Results We analyzed 244 patients with PAH, 54 with ILD, and 47 with ILD-PH. The anchor-based classifier demonstrated high train-test performance with a c-statistic of 0.81 and 78% agreement with expert-adjudicated diagnoses. Among the ILD-PH patients, 23 were classified as PAH-like and 24 as ILD-like (74% agreement with expert adjudication, Cohen’s kappa 0.49, Hedge’s g estimate of 0.86). FELLA enrichment revealed core glycolytic and bioenergetic pathways were enriched in both pure PAH and ILD-PH. Comparison of diffusion Z-scores between groups demonstrated subtle but biologically coherent differences in several reaction-level nodes, including lactate, NADPH oxidoreductase, and phosphofructokinase, which were enriched in ILD-PH. Conclusion Patients with ILD-PH demonstrate clinically relevant metabolomic phenotypes. Examination of the network biology underlying ILD-PH may point to potential therapeutic targets in this vulnerable subgroup with lagging clinical outcomes. This abstract is funded by: Pulmonary Fibrosis Foundation PFF Scholars Program

Abstract Rationale Pulmonary hypertension (PH) is a heterogeneous condition classified into 5 WSPH groups but lacking detailed molecular distinction. The PVDOMICS cohort includes all PH groups alongside healthy and disease comparators, is deeply clinically phenotyped and has RNA sequencing data. Pseudogenes are conserved genome elements that due to sequence-defects cannot be translated into functional proteins. Some pseudogenes never-the-less produce RNA-transcripts which can regulate the expression of other, protein coding genes. Therefore, our study was focused on identifying differentially expressed pseudogene-transcripts within the dataset from the PVDOMICS study and gaining insights into their potential function. Methods Peripheral venous whole blood was collected from PVDOMICS participants [n = 525 PH patients; n = 312 controls without PH who were comparators (n = 256) or healthy (n = 56)]. RNA was isolated, RNA sequencing was performed. After quality control, differential transcript expression between PH and control groups was analyzed using the limma-voom pipeline, rSEA with adjustment for age, sex, race, smoking, WSPH group, latent factors and false discovery rate. We focused on a) transcripts expressed at increased levels (n = 218): the first 200 DEG (differentially-expressed-genes)-transcripts based on log2 difference plus increased DEG within the first 200 most significant p-values; and b) transcripts expressed at decreased levels (n = 292): last 200 DEG based on log2 difference plus decreased DEG within the first 200 most significant p-values. Results Six pseudogene transcripts were among the increased and 35 pseudogene transcripts were among the decreased DEG in the PH group compared to controls. Some of the differentially expressed pseudogenes (DEP) belonged to gene families, like the ribosomal protein family (7 DEP), mitochondria expressed transcripts (3 DEP), zinc finger proteins (2 DEP). Associations of pseudogenes with biological function or disease were extracted via GeneCards(R) and information from GWAS analysis. The DUSP8P5 pseudogene was among the DEP and is known for its association with left heart disease and hypertension; 2 additional DEP have known association with hypertension. Further known associations were with inhibitory or inflammatory immune processes (8 DEP), bone density (6 DEP), metabolism and diabetes (6 DEP), respiratory function or disease (4 DEP), erythrocyte biology (3 DEP), cancer (5 DEP), kidney disease (1 DEP), and demographic parameters like eye or hair color (2 DEP). Ten DEP had no known association with function or disease. Conclusions We identified pseudogene RNA-transcripts associated with PH across WSPH groups. The local vs. systemic effects of these pseudogene transcripts on PH remain to be determined and regulators of pseudogene-transcript expression (e.g. environmental exposures) identified. This abstract is funded by: 1 X01 HL148611 (TopMed), PVDOMICS NHLBI, PVDOMICS Pulmonary Hypertension Association, HT9425-25-1-0533 (GG), 1R21OH012955-01 (GG)

BACKGROUND: The metabolic mechanisms underlying right ventricular (RV) dysfunction are poorly understood, particularly outside of group 1 pulmonary hypertension (PH). We aimed to identify metabolites and pathways associated with RV systolic function and explored whether associations differed by pulmonary vascular resistance, PH group 1 status, and sex. METHODS: We analyzed data from the multicenter PVDOMICS (Pulmonary Vascular Disease Phenomics) cohort. RV systolic function metrics included fractional area change (echo), global longitudinal strain (echo), and ejection fraction (cardiac magnetic resonance). We used linear regression adjusted for age, sex, body mass index, and PH group to assess associations between metabolites and RV function. Pathway enrichment analyses were used to identify pathways significantly associated with RV function. Interaction terms were assessed to determine whether metabolite associations were modified by pulmonary vascular resistance, group 1 PH status, or sex. Least absolute shrinkage and selection operator regression was used to develop metabolite-based scores for RV function, and prognostic performance was assessed. RESULTS: There were 979 participants with plasma metabolomics and RV function data. Linear regression identified 170 metabolites that were significantly associated with all 3 RV metrics. Androgenic steroid, gamma-glutamyl amino acid, polyamine, vitamin A, fatty acid, and sterol pathways are most strongly associated with RV systolic function. Two metabolites interacted with group 1 PH status, and 6 interacted with pulmonary vascular resistance. Four androgenic steroids are associated more strongly with RV systolic function in women compared with men. Metabolite-based scores were prognostically equivalent to RV systolic function metrics and less accurate than REVEAL Lite 2 scores. CONCLUSIONS: We provide a blueprint of metabolites and metabolic pathways associated with RV systolic function across the spectrum of PH. Novel links to vitamin A and glutathione metabolites were observed. We detected few metabolites that associated with RV systolic function differentially by group 1 PH status or degree of pulmonary vascular resistance elevation. Androgenic steroids may associate more strongly with RV systolic function in women compared with men.

Abstract Rationale Pulmonary hypertension (PH) is associated with significant morbidity and mortality. Ground-glass opacities (GGOs) are common in Group 1 PH, but their clinical significance is unclear. Objectives We sought to characterize the clinical features and outcomes of patients with Group 1 PH with and without GGOs in the PVDOMICS study, a prospective multicenter cohort study aimed at deep phenotyping PH. Methods Patients with incident and prevalent PH were enrolled across seven U.S. centers. We included patients with Group 1 PH and excluded those with parenchymal lung disease or without chest imaging, resulting in a cohort of 242 patients. Results GGOs were common among patients with Group 1 PH (43% prevalence), associated with female sex, younger age, prostanoid use, and longer disease duration. GGOs were more common among patients with familial pulmonary arterial hypertension and pulmonary veno-occlusive disease. GGOs were associated with established markers of disease severity, including echocardiographic (right ventricular systolic pressure and tricuspid annular plane systolic excursion), biomarkers (N-terminal pro B-type natriuretic peptide), and worse hemodynamics (higher mean pulmonary artery pressure, pulmonary vascular resistance, and pulmonary artery wedge pressure). GGOs were associated with worse transplant-free survival (hazard ratio, 2.49; 95% confidence interval = 1.43–4.32; P = 0.001) and had independent prognostic value for predicting transplant-free survival after adjusting for European Society of Cardiology and European Respiratory Society risk stratification (hazard ratio, 2.19; 95% confidence interval = 1.20–3.99; P = 0.01). Conclusions Overall, GGOs were associated with specific clinical characteristics and disease phenotypes, as well as worse hemodynamics, longer disease duration, prostanoid use, and worse survival. Future studies evaluating the pathophysiology and “omic” correlates of GGOs are warranted. Clinical trial registered with www.clinicaltrials.gov (NCT02980887).

Abstract Rationale A subset of patients with group 1 pulmonary hypertension (PH) have superimposed left heart abnormalities with unclear metabolic implications. Objectives To compare serum/transpulmonary metabolome between group 1 PH stratified by heart failure with preserved ejection fraction (HFpEF) probability. Methods Patients with group 1 PH were stratified into low (<25%) and high (⩾75%) HFpEF-ABA (age, body mass index, and atrial fibrillation) probability, with healthy control subjects and subjects with clinical HFpEF used for comparison of venous and transpulmonary metabolomics. Measurements and Main Results Group 1 PH + high HFpEF probability (n = 131) was associated with a significant increase in 207 metabolites (false discovery rate [FDR] P < 0.05 and fold change >1) (n = 193, t test) and a significant decrease in 231 metabolites (FDR P < 0.05 and fold change <1) (n = 193, t test) compared with group 1 PH + low HFpEF probability (n = 62). Group 1 PH + high HFpEF probability was associated with enhanced tryptophan metabolism with higher downstream kynurenine metabolite concentrations and lower serotonin concentrations (FDR P < 0.002 for all, n = 193, t test). Linoleate (precursor to arachidonic acid and prostaglandins) and arginine and homoarginine (precursors to nitric oxide) were all lower in group 1 PH + high HFpEF probability (FDR P < 0.03 for all, n = 193, t test). Metabolome changes in group 1 PH + high HFpEF probability overlapped with clinical HFpEF (n = 240) but were abnormal relative to control subjects (n = 85) (P < 0.0001 for all, n = 456, t test). There was no evidence of differential transpulmonary uptake/release of most metabolites, suggesting probable nonpulmonary origin (except for serotonin, interaction P = 0.04; and kynurenine, interaction P = 0.03; n = 433, mixed model). Conclusions Patients with group 1 PH + high HFpEF probability have a unique metabolome characterized by enhanced tryptophan-kynurenine pathway breakdown, deficiency of amino acids (such as glycine and serine), lower serotonin, and decreased prostaglandin and nitric oxide precursors. Despite fulfilling clinical criteria for group 1 PH, these metabolome changes were comparable with clinical HFpEF, supporting biological overlap between these two forms of PH.

Background Right ventricular (RV) maladaptation to elevated pulmonary afterload is the primary determinant of outcomes in pulmonary artery (PA) hypertension; however, the pathobiological mechanisms underlying RV decompensation remain poorly understood. Methods We performed global untargeted metabolomics on plasma from 55 patients who underwent gold‐standard RV‐PA coupling measurements using multibeat pressure volume loop assessment in a single‐center cohort and from 1027 patients with coupling surrogate measurements in a larger multicenter cohort, the PVDOMICS (Pulmonary Vascular Disease Phenomics) study. Age and sex‐adjusted linear regression was performed to identify associations between metabolites and coupling metrics. Additionally, we performed a metabolic flux analysis using gene expression data from RV tissue in an independent cohort of 32 patients. Partial least squares–discriminant analysis was used to identify metabolites and reactions characteristic of the decompensated RV. Results RV‐PA coupling was negatively associated with tricarboxylic acid (TCA) cycle intermediate levels. Specifically, plasma α‐ketoglutarate and fumarate were significantly associated with all coupling metrics in both cohorts. Metabolic flux analysis indicated that decompensated RVs exhibited aberrant TCA cycle activity, including reduced acetyl coenzyme A entry and increased lactate elimination, suggesting a shift from the TCA cycle toward glycolysis at the RV tissue level. Conclusions We identify an association between circulating TCA cycle intermediate levels and RV‐PA uncoupling in 2 independent cohorts, and dysregulated TCA cycle metabolism in decompensated PA hypertension RVs, suggesting that aberrant TCA cycle metabolism could represent a hallmark of RV maladaptation in PA hypertension. Further study of this pathway is warranted to develop novel biomarkers of RV function or RV‐targeted therapies.

Aims: Alterations of mitochondrial bioenergetics and arginine metabolism are universally present and mechanistically linked to pulmonary arterial hypertension (PAH), but there is little knowledge of arginine metabolism and mitochondrial functions across the different pulmonary hypertension (PH) groups. We hypothesize that abnormalities in mitochondrial functions are present across all PH groups and associated with clinical phenotypes. We test the hypothesis in PH patients and healthy controls from the Pulmonary Vascular Disease Phenomics Program cohort, who had comprehensive clinical phenotyping and follow-up for at least 4 years for death or transplant status. Mitochondrial transmembrane potential, superoxide production, and mass were measured by flow cytometry in fresh platelets. Metabolomics analysis was performed on plasma samples. Global arginine bioavailability was calculated as the ratio of arginine/(ornithine+citrulline). Results: Global arginine bioavailability is consistently lower than controls in all PH groups. Although the mitochondrial mass is similar across all PH groups and controls, superoxide production and transmembrane potential vary across groups. Mitochondrial superoxide is higher in group 1 PAH and lowest in group 3 compared with other groups, while transmembrane potential is lower in group 1 PAH than controls or group 3. The alterations in mitochondrial functions of group 1 PAH are associated with changes in fatty acid metabolism. Mitochondrial transmembrane potential in group 1 PAH is associated with transplant-free survival. Conclusion: While alterations in mitochondrial function are found in all PH groups, group 1 PAH has a unique mitochondrial phenotype with greater superoxide and lower transmembrane potential linked to fatty acid metabolism, and clinically to survival. Antioxid. Redox Signal. 42, 361–377.