Barnard, John

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 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.

Our study shows an early increase in kynurenine pathway metabolism in at-risk subjects with systemic sclerosis who develop pulmonary arterial hypertension (PAH). We show that kynurenine pathway upregulation precedes clinical diagnosis and that this metabolic shift is associated with increased disease severity and shorter survival times. We also show that gene expression of TDO2, an enzyme that generates kynurenine from tryptophan, rises with PAH development.

Background: While NKX2-5 plays pivotal roles in human cardiac development and disease, no ChIP-Seq studies of NKX2-5 in human cardiac tissues currently exist, resulting in an incomplete understanding of its direct gene targets. Modern computational methods which identify binding sites using transcription factor motif and tissue-specific experimental data can help to fill this knowledge gap. Objective: To use computational methods to identify likely NKX2-5 binding sites and downstream gene targets using human heart-specific experimental data. Methods: Human cardiomyocyte DNAse hypersensitivity data (2 replicates) were downloaded from the Encyclopedia of DNA Elements (ENCODE) database. The position weight matrix (PWM) representing the transcription factor motif of NKX2-5 was obtained from JASPAR. We applied the Protein Interaction Quantification (PIQ) algorithm to detect NKX2-5 binding sites using the PWM and DNAse hypersensitivity data as inputs. RNA-Seq data from 108 human heart-specific samples (atrial appendages and left ventricles) were downloaded from the Genotype-Tissue Expression (GTEx) database. Protein-coding genes significantly expressed in the heart (RPKM <= 1 based on the GTEx RNA-Seq data) that were within 100kb of the predicted binding sites were then identified. Pathway analysis of these genes was performed using Ingenuity Pathway Analysis (IPA). Results: 1283 binding sites for NKX2-5 were discovered by PIQ. Of 12698 protein-coding, heart-expressed genes, 625 were within 100kb of these binding sites. The identified genes were highly enriched in physiologic categories like “Vasculogenesis” and “Development of Cardiovascular Tissue” (p = 2.36x10-9 and 2.31x10-8 respectively). Notable genes included: cardiac transcription factors (MEF2A, TBX20); growth factors (TGFB2, BMP2); muscle and ion channel function (ACTA2, BIN1), and; calcium signaling (CALM2, CAMK2D). Conclusion: By using computational analyses of transcription factor motif and human heart-specific experimental data, we have identified candidate downstream targets of NKX2-5. Future work will include validation studies in an external cohort and analysis of associations between these candidate genes and cardiac disease.

Abstract Traditional Hardy–Weinberg equilibrium (HWE) tests (the χ2 test and the exact test) have long been used as a metric for evaluating genotype quality, as technical artifacts leading to incorrect genotype calls often can be identified as deviations from HWE. However, in data sets composed of individuals from diverse ancestries, HWE can be violated even without genotyping error, complicating the use of HWE testing to assess genotype data quality. In this manuscript, we present the Robust Unified Test for HWE (RUTH) to test for HWE while accounting for population structure and genotype uncertainty, and to evaluate the impact of population heterogeneity and genotype uncertainty on the standard HWE tests and alternative methods using simulated and real sequence data sets. Our results demonstrate that ignoring population structure or genotype uncertainty in HWE tests can inflate false-positive rates by many orders of magnitude. Our evaluations demonstrate different tradeoffs between false positives and statistical power across the methods, with RUTH consistently among the best across all evaluations. RUTH is implemented as a practical and scalable software tool to rapidly perform HWE tests across millions of markers and hundreds of thousands of individuals while supporting standard VCF/BCF formats. RUTH is publicly available at https://www.github.com/statgen/ruth.

Key Points In the largest study of HMAs in RR-AML to date, 16% of patients achieved CR/CRi and experienced a median OS of 21 months. Low proliferative disease (peripheral blood blasts <5%) was associated with improved response and OS.