Clarifying the Pulmonary Arterial Hypertension Molecular Landscape Using Functional Genetics.

Pulmonary arterial hypertension (PAH) is a complex cardiopulmonary disease that is associated with numerous pathogenetic molecular mechanisms and results in mixed hypertrophic, plexigenic, and fibrotic vascular remodeling of distal pulmonary arterioles. Enhanced clinician awareness and early implementation of multiple PAH-specific therapies have improved the 3-year survival rate to 84% from 52% in the prior era (1). Nonetheless, PAH remains highly morbid, including impaired health-related quality of life that is akin to that of chronic obstructive pulmonary disease, particularly regarding physical inactivity and mental health burden (2). Despite widely heterogenous pathobiology, approved PAH medical therapies (still) target only nitric oxide, prostacyclin, or endothelin receptor biology. Furthermore, treatment responsiveness to PAH pharmacotherapies is highly variable even under tightly controlled circumstances customary among randomized clinical trials, leaving no doubt that as-yet undiscovered therapeutic targets exist by which to subgroup patients and modify their clinical course.

Precision-based methods for diagnosing and prognosticating PAH have focused largely on single genetic variants. In 2001, Newman and colleagues leveraged the wider availability of gene sequencing to complete an observational cohort study spanning 20 years and reported that a thymine-to-guanine transversion at position 354 in exon 3 of the BMPR2 gene was present in 18 families with PAH (3). This finding gave rise to the era of hereditary PAH and, ultimately, the description of 17 disease-causing variants (4) and important advances using genetics for PAH diagnosis, prognosis, and family screening (5). However, <30% of patients have single variants in causative genes, and posttranscriptional mechanisms in numerous cell types have been reported in PAH (4). Together, these findings suggest that, akin to other complex disorders, it is unlikely a single sentinel genetic event underlies the entire PAH phenotypic spectrum.

In 1995, findings from the first bona fide microarray technology were published by Schena and colleagues using a high-speed robotic printing of complementary DNAs on glass (6). Transcriptomic platforms have expanded greatly since then in both sophistication and availability. Greater reliance on multiplex big data platforms, however, has not necessarily been coupled with definitive progress in understanding the mechanistic basis of disease (7). Indeed, data on differentially expressed genes from array probes have been published widely in PAH, although these outputs do not in and of themselves inform the pathobiological function of specific transcripts, and numerous examples showing an uncoupling between transcript quantity and disease relevance exist.

These shortcomings in PAH science establish the following major objective for our field in the modern era: integrating genetic context with transcriptomic data to identify causative pathways underlying the clinical phenotype. In this issue of the Journal, Rhodes and colleagues (pp. 586–594) (8) use a comprehensive, clever, and sound approach to address this challenge head-on. The investigators studied the transcriptomic profile of 359 patients with PAH who were randomized to one of the following three data analysis groups: RNA discovery, RNA validation, and analytical model validation. Each of the three groups was compared with transcriptomic data from 24 distinct age- and sex-matched volunteer control subjects. They identified 507 transcripts that were differentially expressed relative to control subjects in both the discovery and validation cohorts. A LASSO (Least Absolute Shrinkage and Selection Operator) analysis, which is a statistical model reported originally in the geophysics literature to optimize linear regression fitting for variable selection (9), was then used to identify a combination of 25 RNAs that optimally discriminated patients with PAH from control subjects. This approach successfully stratified patients into low- and high-risk groups using survival as an endpoint. Additional outcome analyses yielded internally consistent findings; the RNA signature also associated with World Health Organization functional class, 6-minute-walk distance, and biochemical evidence of heart failure.

Pathway analyses affirmed that many differentially expressed RNAs share annotated function with established PAH pathobiological mechanisms, including HIF-1 signaling, DNA repair, and zinc finger–containing transcription factors (10). However, Rhodes and colleagues recognized that despite this impressive synthesis of data, additional steps were needed to decipher a molecular cue with causative bearing on PAH, and to accomplish this end, they turned to Mendelian randomization (MR). This approach focuses on the effect of genotypic variance on variance in mRNA quantity. The resulting quantitative trait locus (eQTL) map is one basis of functional genetics, which aims to filter out signals in genetic variance that may be associated with a phenotype but are less likely to be pathogenic (thus, more likely associative) (11).

The authors accessed two publicly available eQTL databases and their own previously published PAH genome-wide association study (12) to perform a two-sample MR analysis. From 293 eQTLs available for the 507 RNAs, SMAD5 was one of two genes to reach significance, and investigators focused on a specific SNP (rs4146187). They observed that in PAH, the C/C genotype was associated with decreased SMAD5 mRNA quantity and was present in ∼50% of patients with PAH, whereas the A/A genotype was linked to increased transcript quantity and greater risk reduction for developing PAH. By focusing their method on functional analyses (e.g., eQTL), the results provide a measure of specificity and boost confidence that modifying SMAD5, in this case, indeed modulates the clinical phenotype.

Identifying the relevance of SMAD5 to PAH is an important step forward, but clarifying the mechanistic implications of this finding nonetheless requires additional experimental data. As the authors assert, analyzing transcript quantity does not account for protein posttranslational modifications that are important in PAH (13) and also known to regulate SMAD5 bioactivity (14). Future avenues of research should consider transcriptomic biomarkers to predict PAH pharmacotherapy selection, escalation, or discontinuation. This, in turn, has further implications for PAH clinical trial design and patient enrollment.

Overall, Rhodes and colleagues transform the scientific landscape in PAH by expanding the continuum of biological data used to inform clinical risk. Through a multilayered and comprehensive approach culminating in MR methodology that emphasizes functional genetics via eQTL analysis, transcriptomic array data narrow toward causative molecular pathways. This work, therefore, advances knowledge on the genomic–transcriptomic axis in PAH while identifying SMAD5 and its transcript per se as novel potential therapeutic targets. Further evidence to support these data and repurpose this approach to clarify other PAH subgroups, including differences across the temporal spectrum of the disease, are just a sampling of exciting future opportunity suggested by this important work.