Blood DNA methylation profiling identifies cathepsin Z dysregulation in pulmonary arterial hypertension

Pulmonary arterial hypertension (PAH) is characterised by pulmonary vascular remodelling causing premature death from right heart failure. Established DNA variants influence PAH risk, but susceptibility from epigenetic changes is unknown. We addressed this through epigenome-wide association study (EWAS), testing 865,848 CpG sites for association with PAH in 429 individuals with PAH and 1226 controls. Three loci, at Cathepsin Z (CTSZ, cg04917472), Conserved oligomeric Golgi complex 6 (COG6, cg27396197), and Zinc Finger Protein 678 (ZNF678, cg03144189), reached epigenome-wide significance (p < 10-7) and are hypermethylated in PAH, including in individuals with PAH at 1-year follow-up. Of 16 established PAH genes, only cg10976975 in BMP10 shows hypermethylation in PAH. Hypermethylation at CTSZ is associated with decreased blood cathepsin Z mRNA levels. Knockdown of CTSZ expression in human pulmonary artery endothelial cells increases caspase-3/7 activity (p < 10-4). DNA methylation profiles are altered in PAH, exemplified by the pulmonary endothelial function modifier CTSZ, encoding protease cathepsin Z.

Understanding what drives genetic study participation: Perspectives of patients, carers, and relatives

Genetic research's growing importance in understanding pulmonary arterial hypertension (PAH) and developing effective treatments prompted the RAPID-PAH study. This study sought feedback from stakeholders who participated in two genomic studies to enhance genetic study delivery and clinical integration. Stakeholders from nine UK PH centres, representing various roles, ages, genders, and mutation statuses, took part in 53 semi-structured interviews and focus groups. Transcripts were thematically coded using inductive analysis. Clustering analysis was conducted to identify patient groups that shared attitudes. In this paper, we focus on patients', carers', and relatives' perspectives. The key interview themes revealed insights into participants' attitudes toward genetic research and testing more generally, expertise and knowledge of the disease itself, motivations and barriers to participating in genetic research, awareness of and interest in consent procedures and the use of personal and genetic data, as well as the process of communicating individual genetic results. Factors influencing genetic research participation included altruistic motives, personal diagnostic experiences, and family-related hopes. Clustering analysis produced distinct clusters based on the presence of barriers and motivators for research participation; however, hardly any patients shared identical sets of attitudes, emphasising the need for personalised approaches to recruitment. Most patients reported poor engagement with study-related materials. Patients who received individual genetic results expressed satisfaction with the process, whereas those who did not were disappointed with the lack of feedback. Reflecting on patient perspectives, we offer recommendations to improve the genetic study delivery process. Enhancing genetic research integration into clinical practice requires tailored engagement, clear communication, and support from healthcare stakeholders.

The Medical Action Ontology: A tool for annotating and analyzing treatments and clinical management of human disease

BACKGROUND

Navigating the clinical literature to determine the optimal clinical management for rare diseases presents significant challenges. We introduce the Medical Action Ontology (MAxO), an ontology specifically designed to organize medical procedures, therapies, and interventions.

METHODS

MAxO incorporates logical structures that link MAxO terms to numerous other ontologies within the OBO Foundry. Term development involves a blend of manual and semi-automated processes. Additionally, we have generated annotations detailing diagnostic modalities for specific phenotypic abnormalities defined by the Human Phenotype Ontology (HPO). We introduce a web application, POET, that facilitates MAxO annotations for specific medical actions for diseases using the Mondo Disease Ontology.

FINDINGS

MAxO encompasses 1,757 terms spanning a wide range of biomedical domains, from human anatomy and investigations to the chemical and protein entities involved in biological processes. These terms annotate phenotypic features associated with specific disease (using HPO and Mondo). Presently, there are over 16,000 MAxO diagnostic annotations that target HPO terms. Through POET, we have created 413 MAxO annotations specifying treatments for 189 rare diseases.

CONCLUSIONS

MAxO offers a computational representation of treatments and other actions taken for the clinical management of patients. Its development is closely coupled to Mondo and HPO, broadening the scope of our computational modeling of diseases and phenotypic features. We invite the community to contribute disease annotations using POET (https://poet.jax.org/). MAxO is available under the open-source CC-BY 4.0 license (https://github.com/monarch-initiative/MAxO).

FUNDING

NHGRI 1U24HG011449-01A1 and NHGRI 5RM1HG010860-04.

Natural history of chronic thromboembolic pulmonary disease with no or mild pulmonary hypertension

BACKGROUND

We describe baseline characteristics, disease progression and mortality in chronic thromboembolic pulmonary disease patients as a function of mean pulmonary artery pressure (mPAP) according to new and previous definitions of pulmonary hypertension.

METHODS

All patients diagnosed with chronic thromboembolic pulmonary disease between January, 2015 and December, 2019 were dichotomized according to initial mPAP: ≤ 20 mmHg ('normal') vs 21-24 mmHg ('mildly-elevated'). Baseline features were compared between the groups, and pairwise analysis performed to determine changes in clinical endpoints at 1-year, excluding those who underwent pulmonary endarterectomy or did not attend follow-up. Mortality was assessed for the whole cohort over the entire study period.

RESULTS

One hundred thirteen patients were included; 57 had mPAP ≤ 20 mmHg and 56 had mPAP 21-24 mmHg. Normal mPAP patients had lower pulmonary vascular resistance (1.6 vs 2.5WU, p < 0.01) and right ventricular end-diastolic pressure (5.9 vs 7.8 mmHg, p < 0.01) at presentation. At 3 years, no major deterioration was seen in either group. No patients were treated with pulmonary artery vasodilators. Eight had undergone pulmonary endarterectomy. Over 37 months median follow-up, mortality was 7.0% in the normal mPAP group and 8.9% in the mildly-elevated mPAP group. Cause of death was malignancy in 62.5% of cases.

CONCLUSIONS

Chronic thromboembolic pulmonary disease patients with mild pulmonary hypertension have statistically higher right ventricular end-diastolic pressure and pulmonary vascular resistance than those with mPAP ≤ 20 mmHg. Baseline characteristics were otherwise similar. Neither group displayed disease progression on non-invasive tests up to 3 years. Mortality over 37 months follow-up is 8%, and mainly attributable to malignancy. Further prospective study is required to validate these findings.

The Medical Action Ontology: A Tool for Annotating and Analyzing Treatments and Clinical Management of Human Disease

Navigating the vast landscape of clinical literature to find optimal treatments and management strategies can be a challenging task, especially for rare diseases. To address this task, we introduce the Medical Action Ontology (MAxO), the first ontology specifically designed to organize medical procedures, therapies, and interventions in a structured way. Currently, MAxO contains 1757 medical action terms added through a combination of manual and semi-automated processes. MAxO was developed with logical structures that make it compatible with several other ontologies within the Open Biological and Biomedical Ontologies (OBO) Foundry. These cover a wide range of biomedical domains, from human anatomy and investigations to the chemical and protein entities involved in biological processes. We have created a database of over 16000 annotations that describe diagnostic modalities for specific phenotypic abnormalities as defined by the Human Phenotype Ontology (HPO). Additionally, 413 annotations are provided for medical actions for 189 rare diseases. We have developed a web application called POET (https://poet.jax.org/) for the community to use to contribute MAxO annotations. MAxO provides a computational representation of treatments and other actions taken for the clinical management of patients. The development of MAxO is closely coupled to the Mondo Disease Ontology (Mondo) and the Human Phenotype Ontology (HPO) and expands the scope of our computational modeling of diseases and phenotypic features to include diagnostics and therapeutic actions. MAxO is available under the open-source CC-BY 4.0 license (https://github.com/monarch-initiative/MAxO).

Angiopoietin 2 and hsCRP are associated with pulmonary hemodynamics and long-term mortality respectively in CTEPH-Results from a prospective discovery and validation biomarker study

INTRODUCTION

Chronic thromboembolic pulmonary hypertension (CTEPH) is an underdiagnosed disease of uncertain etiology. Altered endothelial homeostasis, defective angiogenesis and inflammation are implicated. Angiopoietin 2 (Ang2) impairs acute thrombus resolution and is associated with vasculopathy in idiopathic pulmonary arterial hypertension.

METHODS

We assessed circulating proteins associated with these processes in serum from patients with CTEPH (n = 71) before and after pulmonary endarterectomy (PEA), chronic thromboembolic pulmonary disease without pulmonary hypertension (CTEPD, n = 9) and healthy controls (n = 20) using Luminex multiplex arrays. Comparisons between groups were made using multivariable rank regression models. Ang2 and high-sensitivity C-reactive protein (hsCRP) were measured in a larger validation dataset (CTEPH = 277, CTEPD = 26). Cox proportional hazards models were used to identify markers predictive of survival.

RESULTS

In CTEPH patients, Ang2, interleukin (IL) 8, tumor necrosis factor α, and hsCRP were elevated compared to controls, while vascular endothelial growth factor (VEGF) c was lower (p < 0.05). Ang2 fell post-PEA (p < 0.05) and was associated with both pre- and post-PEA pulmonary hemodynamic variables and functional assessments (p < 0.05). In the validation dataset, Ang2 was significantly higher in CTEPH compared to CTEPD. Pre-operative hsCRP was an independent predictor of mortality.

CONCLUSIONS

We hypothesize that CTEPH patients have significant distal micro-vasculopathy and consequently high circulating Ang2. Patients with CTEPD without pulmonary hypertension have no discernible distal micro-vasculopathy and therefore have low circulating Ang2. This suggests Ang2 may be critical to CTEPH disease pathogenesis (impaired thrombus organization and disease severity).

First Genotype-Phenotype Study in TBX4 Syndrome: Gain-of-Function Mutations Causative for Lung Disease

Rationale: Despite the increased recognition of TBX4 (T-BOX transcription factor 4)-associated pulmonary arterial hypertension (PAH), genotype-phenotype associations are lacking and may provide important insights. Objectives: To compile and functionally characterize all TBX4 variants reported to date and undertake a comprehensive genotype-phenotype analysis. Methods: We assembled a multicenter cohort of 137 patients harboring monoallelic TBX4 variants and assessed the pathogenicity of missense variation (n = 42) using a novel luciferase reporter assay containing T-BOX binding motifs. We sought genotype-phenotype correlations and undertook a comparative analysis with patients with PAH with BMPR2 (Bone Morphogenetic Protein Receptor type 2) causal variants (n = 162) or no identified variants in PAH-associated genes (n = 741) genotyped via the National Institute for Health Research BioResource-Rare Diseases. Measurements and Main Results: Functional assessment of TBX4 missense variants led to the novel finding of gain-of-function effects associated with older age at diagnosis of lung disease compared with loss-of-function effects (P = 0.038). Variants located in the T-BOX and nuclear localization domains were associated with earlier presentation (P = 0.005) and increased incidence of interstitial lung disease (P = 0.003). Event-free survival (death or transplantation) was shorter in the T-BOX group (P = 0.022), although age had a significant effect in the hazard model (P = 0.0461). Carriers of TBX4 variants were diagnosed at a younger age (P < 0.001) and had worse baseline lung function (FEV1, FVC) (P = 0.009) than the BMPR2 and no identified causal variant groups. Conclusions: We demonstrated that TBX4 syndrome is not strictly the result of haploinsufficiency but can also be caused by gain of function. The pleiotropic effects of TBX4 in lung disease may be in part explained by the differential effect of pathogenic mutations located in critical protein domains.

A pilot study to examine association of BMI with functional class and 6 min walk distance in idiopathic and heritable PAH: Possible association with estrogen metabolism

The hypothesis that a relationship exists between body mass index (BMI), functional class, and 6 min walk distance (6MWD) in Group 1-pulmonary arterial hypertension (PAH) was examined. Analysis of data from the UK National Cohort Study for heritable pulmonary arterial/idiopathic PAH suggests increased BMI is a predictor of worse functional class and shorter 6MWD; increased body-weight in mice and man may be associated with increased estrogen metabolism.

Autoimmunity Is a Significant Feature of Idiopathic Pulmonary Arterial Hypertension

Rationale: Autoimmunity is believed to play a role in idiopathic pulmonary arterial hypertension (IPAH). It is not clear whether this is causative or a bystander of disease and if it carries any prognostic or treatment significance. Objectives: To study autoimmunity in IPAH using a large cross-sectional cohort. Methods: Assessment of the circulating immune cell phenotype was undertaken using flow cytometry, and the profile of serum immunoglobulins was generated using a standardized multiplex array of 19 clinically validated autoantibodies in 473 cases and 946 control subjects. Additional glutathione S-transferase fusion array and ELISA data were used to identify a serum autoantibody to BMPR2 (bone morphogenetic protein receptor type 2). Clustering analyses and clinical correlations were used to determine associations between immunogenicity and clinical outcomes. Measurements and Main Results: Flow cytometric immune profiling demonstrates that IPAH is associated with an altered humoral immune response in addition to raised IgG3. Multiplexed autoantibodies were significantly raised in IPAH, and clustering demonstrated three distinct clusters: "high autoantibody," "low autoantibody," and a small "intermediate" cluster exhibiting high concentrations of ribonucleic protein complex. The high-autoantibody cluster had worse hemodynamics but improved survival. A small subset of patients demonstrated immunoglobulin reactivity to BMPR2. Conclusions: This study establishes aberrant immune regulation and presence of autoantibodies as key features in the profile of a significant proportion of patients with IPAH and is associated with clinical outcomes.

Mining the Plasma Proteome for Insights into the Molecular Pathology of Pulmonary Arterial Hypertension

Rationale: Pulmonary arterial hypertension (PAH) is characterized by structural remodeling of pulmonary arteries and arterioles. Underlying biological processes are likely reflected in a perturbation of circulating proteins. Objectives: To quantify and analyze the plasma proteome of patients with PAH using inherited genetic variation to inform on underlying molecular drivers. Methods: An aptamer-based assay was used to measure plasma proteins in 357 patients with idiopathic or heritable PAH, 103 healthy volunteers, and 23 relatives of patients with PAH. In discovery and replication subgroups, the plasma proteomes of PAH and healthy individuals were compared, and the relationship to transplantation-free survival in PAH was determined. To examine causal relationships to PAH, protein quantitative trait loci (pQTL) that influenced protein levels in the patient population were used as instruments for Mendelian randomization (MR) analysis. Measurements and Main Results: From 4,152 annotated plasma proteins, levels of 208 differed between patients with PAH and healthy subjects, and 49 predicted long-term survival. MR based on cis-pQTL located in proximity to the encoding gene for proteins that were prognostic and distinguished PAH from health estimated an adverse effect for higher levels of netrin-4 (odds ratio [OR], 1.55; 95% confidence interval [CI], 1.16-2.08) and a protective effect for higher levels of thrombospondin-2 (OR, 0.83; 95% CI, 0.74-0.94) on PAH. Both proteins tracked the development of PAH in previously healthy relatives and changes in thrombospondin-2 associated with pulmonary arterial pressure at disease onset. Conclusions: Integrated analysis of the plasma proteome and genome implicates two secreted matrix-binding proteins, netrin-4 and thrombospondin-2, in the pathobiology of PAH.

Using the Plasma Proteome for Risk Stratifying Patients with Pulmonary Arterial Hypertension

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^-6). 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.

Severe Pulmonary Arterial Hypertension Is Characterized by Increased Neutrophil Elastase and Relative Elafin Deficiency

BACKGROUND

Preclinical evidence implicates neutrophil elastase (NE) in pulmonary arterial hypertension (PAH) pathogenesis, and the NE inhibitor elafin is under early therapeutic investigation.

RESEARCH QUESTION

Are circulating NE and elafin levels abnormal in PAH and are they associated with clinical severity?

STUDY DESIGN AND METHODS

In an observational Stanford University PAH cohort (n = 249), plasma NE and elafin levels were measured in comparison with those of healthy control participants (n = 106). NE and elafin measurements were then related to PAH clinical features and relevant ancillary biomarkers. Cox regression models were fitted with cubic spline functions to associate NE and elafin levels with survival. To validate prognostic relationships, we analyzed two United Kingdom cohorts (n = 75 and n = 357). Mixed-effects models evaluated NE and elafin changes during disease progression. Finally, we studied effects of NE-elafin balance on pulmonary artery endothelial cells (PAECs) from patients with PAH.

RESULTS

Relative to control participants, patients with PAH were found to have increased NE levels (205.1 ng/mL [interquartile range (IQR), 123.6-387.3 ng/mL] vs 97.6 ng/mL [IQR, 74.4-126.6 ng/mL]; P < .0001) and decreased elafin levels (32.0 ng/mL [IQR, 15.3-59.1 ng/mL] vs 45.5 ng/mL [IQR, 28.1-92.8 ng/mL]; P < .0001) independent of PAH subtype, illness duration, and therapies. Higher NE levels were associated with worse symptom severity, shorter 6-min walk distance, higher N-terminal pro-type brain natriuretic peptide levels, greater right ventricular dysfunction, worse hemodynamics, increased circulating neutrophil levels, elevated cytokine levels, and lower blood BMPR2 expression. In Stanford patients, NE levels of > 168.5 ng/mL portended increased mortality risk after adjustment for known clinical predictors (hazard ratio [HR], 2.52; CI, 1.36-4.65, P = .003) or prognostic cytokines (HR, 2.63; CI, 1.42-4.87; P = .001), and the NE level added incremental value to established PAH risk scores. Similar prognostic thresholds were identified in validation cohorts. Longitudinal NE changes tracked with clinical trends and outcomes. PAH PAECs exhibited increased apoptosis and attenuated angiogenesis when exposed to NE at the level observed in patients' blood. Elafin rescued PAEC homeostasis, yet the required dose exceeded levels found in patients.

INTERPRETATION

Blood levels of NE are increased while elafin levels are deficient across PAH subtypes. Higher NE levels are associated with worse clinical disease severity and outcomes, and this target-specific biomarker could facilitate therapeutic development of elafin.

Rare variant analysis of 4241 pulmonary arterial hypertension cases from an international consortium implicates FBLN2, PDGFD, and rare de novo variants in PAH

BACKGROUND

Pulmonary arterial hypertension (PAH) is a lethal vasculopathy characterized by pathogenic remodeling of pulmonary arterioles leading to increased pulmonary pressures, right ventricular hypertrophy, and heart failure. PAH can be associated with other diseases (APAH: connective tissue diseases, congenital heart disease, and others) but often the etiology is idiopathic (IPAH). Mutations in bone morphogenetic protein receptor 2 (BMPR2) are the cause of most heritable cases but the vast majority of other cases are genetically undefined.

METHODS

To identify new risk genes, we utilized an international consortium of 4241 PAH cases with exome or genome sequencing data from the National Biological Sample and Data Repository for PAH, Columbia University Irving Medical Center, and the UK NIHR BioResource - Rare Diseases Study. The strength of this combined cohort is a doubling of the number of IPAH cases compared to either national cohort alone. We identified protein-coding variants and performed rare variant association analyses in unrelated participants of European ancestry, including 1647 IPAH cases and 18,819 controls. We also analyzed de novo variants in 124 pediatric trios enriched for IPAH and APAH-CHD.

RESULTS

Seven genes with rare deleterious variants were associated with IPAH with false discovery rate smaller than 0.1: three known genes (BMPR2, GDF2, and TBX4), two recently identified candidate genes (SOX17, KDR), and two new candidate genes (fibulin 2, FBLN2; platelet-derived growth factor D, PDGFD). The new genes were identified based solely on rare deleterious missense variants, a variant type that could not be adequately assessed in either cohort alone. The candidate genes exhibit expression patterns in lung and heart similar to that of known PAH risk genes, and most variants occur in conserved protein domains. For pediatric PAH, predicted deleterious de novo variants exhibited a significant burden compared to the background mutation rate (2.45×, p = 2.5e-5). At least eight novel pediatric candidate genes carrying de novo variants have plausible roles in lung/heart development.

CONCLUSIONS

Rare variant analysis of a large international consortium identified two new candidate genes-FBLN2 and PDGFD. The new genes have known functions in vasculogenesis and remodeling. Trio analysis predicted that ~ 15% of pediatric IPAH may be explained by de novo variants.

Bayesian Inference Associates Rare KDR Variants with Specific Phenotypes in Pulmonary Arterial Hypertension

Background - Approximately 25% of patients with pulmonary arterial hypertension (PAH) have been found to harbor rare mutations in disease-causing genes. To identify missing heritability in PAH we integrated deep phenotyping with whole-genome sequencing data using Bayesian statistics. Methods - We analyzed 13,037 participants enrolled in the NIHR BioResource - Rare Diseases (NBR) study, of which 1,148 were recruited to the PAH domain. To test for genetic associations between genes and selected phenotypes of pulmonary hypertension (PH), we used the Bayesian rare-variant association method BeviMed. Results - Heterozygous, high impact, likely loss-of-function variants in the Kinase Insert Domain Receptor (KDR) gene were strongly associated with significantly reduced transfer coefficient for carbon monoxide (KCO, posterior probability (PP)=0.989) and older age at diagnosis (PP=0.912). We also provide evidence for familial segregation of a rare nonsense KDR variant with these phenotypes. On computed tomographic imaging of the lungs, a range of parenchymal abnormalities were observed in the five patients harboring these predicted deleterious variants in KDR. Four additional PAH cases with rare likely loss-of-function variants in KDR were independently identified in the US PAH Biobank cohort with similar phenotypic characteristics. Conclusions - The Bayesian inference approach allowed us to independently validate KDR, which encodes for the Vascular Endothelial Growth Factor Receptor 2 (VEGFR2), as a novel PAH candidate gene. Furthermore, this approach specifically associated high impact likely loss-of-function variants in the genetically constrained gene with distinct phenotypes. These findings provide evidence for KDR being a clinically actionable PAH gene and further support the central role of the vascular endothelium in the pathobiology of PAH.

'There and Back Again'-Forward Genetics and Reverse Phenotyping in Pulmonary Arterial Hypertension

Although the invention of right heart catheterisation in the 1950s enabled accurate clinical diagnosis of pulmonary arterial hypertension (PAH), it was not until 2000 when the landmark discovery of the causative role of bone morphogenetic protein receptor type II (BMPR2) mutations shed new light on the pathogenesis of PAH. Since then several genes have been discovered, which now account for around 25% of cases with the clinical diagnosis of idiopathic PAH. Despite the ongoing efforts, in the majority of patients the cause of the disease remains elusive, a phenomenon often referred to as "missing heritability". In this review, we discuss research approaches to uncover the genetic architecture of PAH starting with forward phenotyping, which in a research setting should focus on stable intermediate phenotypes, forward and reverse genetics, and finally reverse phenotyping. We then discuss potential sources of "missing heritability" and how functional genomics and multi-omics methods are employed to tackle this problem.

Expression Quantitative Trait Locus Mapping in Pulmonary Arterial Hypertension

Expression quantitative trait loci (eQTL) can provide a link between disease susceptibility variants discovered by genetic association studies and biology. To date, eQTL mapping studies have been primarily conducted in healthy individuals from population-based cohorts. Genetic effects have been known to be context-specific and vary with changing environmental stimuli. We conducted a transcriptome- and genome-wide eQTL mapping study in a cohort of patients with idiopathic or heritable pulmonary arterial hypertension (PAH) using RNA sequencing (RNAseq) data from whole blood. We sought confirmation from three published population-based eQTL studies, including the GTEx Project, and followed up potentially novel eQTL not observed in the general population. In total, we identified 2314 eQTL of which 90% were cis-acting and 75% were confirmed by at least one of the published studies. While we observed a higher GWAS trait colocalization rate among confirmed eQTL, colocalisation rate of novel eQTL reported for lung-related phenotypes was twice as high as that of confirmed eQTL. Functional enrichment analysis of genes with novel eQTL in PAH highlighted immune-related processes, a suspected contributor to PAH. These potentially novel eQTL specific to or active in PAH could be useful in understanding genetic risk factors for other diseases that share common mechanisms with PAH.

Plasma metabolomics exhibit response to therapy in chronic thromboembolic pulmonary hypertension

Pulmonary hypertension is a condition with limited effective treatment options. Chronic thromboembolic pulmonary hypertension (CTEPH) is a notable exception, with pulmonary endarterectomy (PEA) often proving curative. This study investigated the plasma metabolome of CTEPH patients, estimated reversibility to an effective treatment and explored the source of metabolic perturbations.We performed untargeted analysis of plasma metabolites in CTEPH patients compared to healthy controls and disease comparators. Changes in metabolic profile were evaluated in response to PEA. A subset of patients were sampled at three anatomical locations and plasma metabolite gradients calculated.We defined and validated altered plasma metabolite profiles in patients with CTEPH. 12 metabolites were confirmed by receiver operating characteristic analysis to distinguish CTEPH and both healthy (area under the curve (AUC) 0.64-0.94, all p<2×10^-5) and disease controls (AUC 0.58-0.77, all p<0.05). Many of the metabolic changes were notably similar to those observed in idiopathic pulmonary arterial hypertension (IPAH). Only five metabolites (5-methylthioadenosine, N1-methyladenosine, N1-methylinosine, 7-methylguanine, N-formylmethionine) distinguished CTEPH from chronic thromboembolic disease or IPAH. Significant corrections (15-100% of perturbation) in response to PEA were observed in some, but not all metabolites. Anatomical sampling identified 188 plasma metabolites, with significant gradients in tryptophan, sphingomyelin, methionine and Krebs cycle metabolites. In addition, metabolites associated with CTEPH and gradients showed significant associations with clinical measures of disease severity.We identified a specific metabolic profile that distinguishes CTEPH from controls and disease comparators, despite the observation that most metabolic changes were common to both CTEPH and IPAH patients. Plasma metabolite gradients implicate cardiopulmonary tissue metabolism of metabolites associated with pulmonary hypertension and metabolites that respond to PEA surgery could be a suitable noninvasive marker for evaluating future targeted therapeutic interventions.

Functional respiratory imaging identifies redistribution of pulmonary blood flow in patients with COVID-19

An increasing observation is that some patients with COVID-19 have normal lung compliance but significant hypoxaemia different from typical acute respiratory distress syndrome (ARDS). We hypothesised that changes in pulmonary blood distribution may be partially responsible and used functional respiratory imaging on CT scans to calculate pulmonary blood volume. We found that patients with COVID-19 had significantly reduced blood volume in the smaller calibre blood vessels (here defined as <5 mm² cross-sectional area) compared with matched ARDS patients and healthy controls. This suggests that using high levels of PEEP may not alone be enough to oxygenate these patients and that additional management strategies may be needed.

Whole-Blood RNA Profiles Associated with Pulmonary Arterial Hypertension and Clinical Outcome

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.

Characterization of GDF2 Mutations and Levels of BMP9 and BMP10 in Pulmonary Arterial Hypertension

Rationale: Recently, rare heterozygous mutations in GDF2 were identified in patients with pulmonary arterial hypertension (PAH). GDF2 encodes the circulating BMP (bone morphogenetic protein) type 9, which is a ligand for the BMP2 receptor.Objectives: Here we determined the functional impact of GDF2 mutations and characterized plasma BMP9 and BMP10 levels in patients with idiopathic PAH.Methods: Missense BMP9 mutant proteins were expressed in vitro and the impact on BMP9 protein processing and secretion, endothelial signaling, and functional activity was assessed. Plasma BMP9 and BMP10 levels and activity were assayed in patients with PAH with GDF2 variants and in control subjects. Levels were also measured in a larger cohort of control subjects (n = 120) and patients with idiopathic PAH (n = 260).Measurements and Main Results: We identified a novel rare variation at the GDF2 and BMP10 loci, including copy number variation. In vitro, BMP9 missense proteins demonstrated impaired cellular processing and secretion. Patients with PAH who carried these mutations exhibited reduced plasma levels of BMP9 and reduced BMP activity. Unexpectedly, plasma BMP10 levels were also markedly reduced in these individuals. Although overall BMP9 and BMP10 levels did not differ between patients with PAH and control subjects, BMP10 levels were lower in PAH females. A subset of patients with PAH had markedly reduced plasma levels of BMP9 and BMP10 in the absence of GDF2 mutations.Conclusions: Our findings demonstrate that GDF2 mutations result in BMP9 loss of function and are likely causal. These mutations lead to reduced circulating levels of both BMP9 and BMP10. These findings support therapeutic strategies to enhance BMP9 or BMP10 signaling in PAH.

Mendelian randomisation analysis of red cell distribution width in pulmonary arterial hypertension

Pulmonary arterial hypertension (PAH) is a rare disease that leads to premature death from right heart failure. It is strongly associated with elevated red cell distribution width (RDW), a correlate of several iron status biomarkers. High RDW values can signal early-stage iron deficiency or iron deficiency anaemia. This study investigated whether elevated RDW is causally associated with PAH.A two-sample Mendelian randomisation (MR) approach was applied to investigate whether genetic predisposition to higher levels of RDW increases the odds of developing PAH. Primary and secondary MR analyses were performed using all available genome-wide significant RDW variants (n=179) and five genome-wide significant RDW variants that act via systemic iron status, respectively.We confirmed the observed association between RDW and PAH (OR 1.90, 95% CI 1.80-2.01) in a multicentre case-control study (cases n=642, disease controls n=15 889). The primary MR analysis was adequately powered to detect a causal effect (odds ratio) between 1.25 and 1.52 or greater based on estimates reported in the RDW genome-wide association study or from our own data. There was no evidence for a causal association between RDW and PAH in either the primary (OR_causal 1.07, 95% CI 0.92-1.24) or the secondary (OR_causal 1.09, 95% CI 0.77-1.54) MR analysis.The results suggest that at least some of the observed association of RDW with PAH is secondary to disease progression. Results of iron therapeutic trials in PAH should be interpreted with caution, as any improvements observed may not be mechanistically linked to the development of PAH.

Traffic exposures, air pollution and outcomes in pulmonary arterial hypertension: a UK cohort study analysis

While traffic and air pollution exposure is associated with increased mortality in numerous diseases, its association with disease severity and outcomes in pulmonary arterial hypertension (PAH) remains unknown.Exposure to particulate matter with a 50% cut-off aerodynamic diameter ≤2.5 μm (PM_2.5), nitrogen dioxide (NO₂) and indirect measures of traffic-related air pollution (distance to main road and length of roads within buffer zones surrounding residential addresses) were estimated for 301 patients with idiopathic/heritable PAH recruited in the UK National Cohort Study of Idiopathic and Heritable PAH. Associations with transplant-free survival and pulmonary haemodynamic severity at baseline were assessed, adjusting for confounding variables defined a prioriHigher estimated exposure to PM_2.5 was associated with higher risk of death or lung transplant (unadjusted hazard ratio (HR) 2.68 (95% CI 1.11-6.47) per 3 μg·m^-3; p=0.028). This association remained similar when adjusted for potential confounding variables (HR 4.38 (95% CI 1.44-13.36) per 3 μg·m^-3; p=0.009). No associations were found between NO₂ exposure or other traffic pollution indicators and transplant-free survival. Conversely, indirect measures of exposure to traffic-related air pollution within the 500-1000 m buffer zones correlated with the European Society of Cardiology/European Respiratory Society risk categories as well as pulmonary haemodynamics at baseline. This association was strongest for pulmonary vascular resistance.In idiopathic/heritable PAH, indirect measures of exposure to traffic-related air pollution were associated with disease severity at baseline, whereas higher PM_2.5 exposure may independently predict shorter transplant-free survival.

Identification of rare sequence variation underlying heritable pulmonary arterial hypertension

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.