Whole blood RNA profiles associated with pulmonary arterial hypertension and clinical outcome

Rationale: Idiopathic and hereditary 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 signalling pathways and stratify patients more robustly according to clinical risk.

Objectives: Using a three-stage design of RNA discovery, RNA validation/model construction and model validation to define a set of PAH-associated RNAs and a single summarising RNA model score. To define genes most likely to be involved in disease development, we performed Mendelian randomisation (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 eQTL and summary statistics from a PAH GWAS.

Measurements and Main Results: We identified 507 genes with differential RNA expression in PAH patients compared to controls. A model of 25 RNAs was able to distinguish PAH with 87% accuracy (AUC 95% CI: 0.791-0.945) in model validation. The RNA model score was associated with disease severity and long-term survival (p=4.66x10-6) in PAH. MR detected an association between SMAD5 levels and PAH disease susceptibility (OR:0.317, 95%CI: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.