Pulmonary arterial hypertension (PAH) is a rare disorder with a poor prognosis. Deleterious variation within components of the transforming growth factor-β pathway, particularly the bone morphogenetic protein type 2 receptor (BMPR2), underlies most heritable forms of PAH. To identify the missing heritability we perform whole-genome sequencing in 1038 PAH index cases and 6385 PAH-negative control subjects. Case-control analyses reveal significant overrepresentation of rare variants in ATP13A3, AQP1 and SOX17, and provide independent validation of a critical role for GDF2 in PAH. We demonstrate familial segregation of mutations in SOX17 and AQP1 with PAH. Mutations in GDF2, encoding a BMPR2 ligand, lead to reduced secretion from transfected cells. In addition, we identify pathogenic mutations in the majority of previously reported PAH genes, and provide evidence for further putative genes. Taken together these findings contribute new insights into the molecular basis of PAH and indicate unexplored pathways for therapeutic intervention.

Background: Inflammation and dysregulated immunity are important in the development of pulmonary arterial hypertension (PAH). Compelling preclinical data supports the therapeutic blockade of interleukin-6 (IL-6) signalling.Methods: We conducted a phase 2 open-label study of intravenous tocilizumab (8 mg·kg⁻¹) over 6 months in patients with group 1 PAH. Co-primary end-points were safety, defined by incidence and severity of adverse events, and change in pulmonary vascular resistance. Separately, a mendelian randomisation study was undertaken on 11 744 individuals with European ancestry including 2085 patients with idiopathic/heritable disease for the IL-6 receptor (IL6R) variant (rs7529229), known to associate with circulating IL-6R levels.Results: We recruited 29 patients (male/female 10/19; mean±sdage 54.9±11.4 years). Of these, 19 had heritable/idiopathic PAH and 10 had connective tissue disease-associated PAH. Six were withdrawn prior to drug administration; 23 patients received at least one dose of tocilizumab. Tocilizumab was discontinued in four patients owing to serious adverse events. There were no deaths. Despite evidence of target engagement in plasma IL-6 and C-reactive protein levels, both intention-to-treat and modified intention-to-treat analyses demonstrated no change in pulmonary vascular resistance. Inflammatory markers did not predict treatment response. Mendelian randomisation did not support an effect of the leadIL6Rvariant on risk of PAH (OR 0.99, p=0.88).Conclusion: Adverse events were consistent with the known safety profile of tocilizumab. Tocilizumab did not show any consistent treatment effect.

Rationale: NT-proBNP (N-terminal pro–brain natriuretic peptide), a biomarker of cardiac origin, is used to risk stratify patients with pulmonary arterial hypertension (PAH). Its limitations include poor sensitivity to early vascular pathology. Other biomarkers of vascular or systemic origin may also be useful in the management of PAH. Objectives: Identify prognostic proteins in PAH that complement NT-proBNP and clinical risk scores. Methods: An aptamer-based assay (SomaScan version 4) targeting 4,152 proteins was used to measure plasma proteins in patients with idiopathic, heritable, or drug-induced PAH from the UK National Cohort of PAH (n = 357) and the French EFORT (Evaluation of Prognostic Factors and Therapeutic Targets in PAH) study (n = 79). Prognostic proteins were identified in discovery–replication analyses of UK samples. Proteins independent of 6-minute-walk distance and NT-proBNP entered least absolute shrinkage and selection operator modeling, and the best combination in a single score was evaluated against clinical targets in EFORT. Measurements and Main Results: Thirty-one proteins robustly informed prognosis independent of NT-proBNP and 6-minute-walk distance in the UK cohort. A weighted combination score of six proteins was validated at baseline (5-yr mortality; area under the curve [AUC], 0.73; 95% confidence interval [CI], 0.63–0.85) and follow-up in EFORT (AUC, 0.84; 95% CI, 0.75–0.94; P = 9.96 × 10⁻⁶). The protein score risk stratified patients independent of established clinical targets and risk equations. The addition of the six-protein model score to NT-proBNP improved prediction of 5-year outcomes from AUC 0.762 (0.702–0.821) to 0.818 (0.767–0.869) by receiver operating characteristic analysis (P = 0.00426 for difference in AUC) in the UK replication and French samples combined. Conclusions: The plasma proteome informs prognosis beyond established factors in PAH and may provide a more sensitive measure of therapeutic response.

Idiopathic pulmonary arterial hypertension (IPAH) is a rare but fatal disease diagnosed by right heart catheterisation and the exclusion of other forms of pulmonary arterial hypertension, producing a heterogeneous population with varied treatment response. Here we show unsupervised machine learning identification of three major patient subgroups that account for 92% of the cohort, each with unique whole blood transcriptomic and clinical feature signatures. These subgroups are associated with poor, moderate, and good prognosis. The poor prognosis subgroup is associated with upregulation of the ALAS2 and downregulation of several immunoglobulin genes, while the good prognosis subgroup is defined by upregulation of the bone morphogenetic protein signalling regulator NOG, and the C/C variant of HLA-DPA1/DPB1 (independently associated with survival). These findings independently validated provide evidence for the existence of 3 major subgroups (endophenotypes) within the IPAH classification, could improve risk stratification and provide molecular insights into the pathogenesis of IPAH.

Rationale: Idiopathic and heritable pulmonary arterial hypertension (PAH) are rare but comprise a genetically heterogeneous patient group. RNA sequencing linked to the underlying genetic architecture can be used to better understand the underlying pathology by identifying key signaling pathways and stratify patients more robustly according to clinical risk. Objectives: To use a three-stage design of RNA discovery, RNA validation and model construction, and model validation to define a set of PAH-associated RNAs and a single summarizing RNA model score. To define genes most likely to be involved in disease development, we performed Mendelian randomization (MR) analysis. Methods: RNA sequencing was performed on whole-blood samples from 359 patients with idiopathic, heritable, and drug-induced PAH and 72 age- and sex-matched healthy volunteers. The score was evaluated against disease severity markers including survival analysis using all-cause mortality from diagnosis. MR used known expression quantitative trait loci and summary statistics from a PAH genome-wide association study. Measurements and Main Results: We identified 507 genes with differential RNA expression in patients with PAH compared with control subjects. A model of 25 RNAs distinguished PAH with 87% accuracy (area under the curve 95% confidence interval: 0.791–0.945) in model validation. The RNA model score was associated with disease severity and long-term survival (P = 4.66 × 10−6) in PAH. MR detected an association between SMAD5 levels and PAH disease susceptibility (odds ratio, 0.317; 95% confidence interval, 0.129–0.776; P = 0.012). Conclusions: A whole-blood RNA signature of PAH, which includes RNAs relevant to disease pathogenesis, associates with disease severity and identifies patients with poor clinical outcomes. Genetic variants associated with lower SMAD5 expression may increase susceptibility to PAH.