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.

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.

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.

The Vps10p domain receptor SorCS2 is crucial for the development and function of the nervous system and essential for brain-derived neurotrophic factor (BDNF)-induced changes in neuronal morphology and plasticity. SorCS2 regulates the subcellular trafficking of the BDNF signaling receptor TrkB as well as selected neurotransmitter receptors in a manner that is dependent on the SorCS2 intracellular domain (ICD). However, the cellular machinery and adaptor protein (AP) interactions that regulate receptor trafficking via the SorCS2 ICD are unknown. We here identify four splice variants of human SorCS2 differing in the insertion of an acidic cluster motif and/or a serine residue within the ICD. We show that each variant undergoes posttranslational proteolytic processing into a one- or two-chain receptor, giving rise to eight protein isoforms, the expression of which differs between neuronal and nonneuronal tissues and is affected by cellular stressors. We found that the only variants without the serine were able to rescue BDNF-induced branching of SorCS2 knockout hippocampal neurons, while variants without the acidic cluster showed increased interactions with clathrin-associated APs AP-1, AP-2, and AP-3. Using yeast two-hybrid screens, we further discovered that all variants bound dynein light chain Tctex-type 3; however, only variants with an acidic cluster motif bound kinesin light chain 1. Accordingly, splice variants showed markedly different trafficking properties and localized to different subcellular compartments. Taken together, our findings demonstrate the existence of eight functional SorCS2 isoforms with differential capacity for interactions with cytosolic ligands dynein light chain Tctex-type 3 and kinesin light chain 1, which potentially allows cell-type specific SorCS2 trafficking and BDNF signaling.

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.

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

Rare genetic variants cause pulmonary arterial hypertension, but the contribution of common genetic variation to disease risk and natural history is poorly characterised. We tested for genome-wide association for pulmonary arterial hypertension in large international cohorts and assessed the contribution of associated regions to outcomes.

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.

Attention deficit hyperactivity disorder (ADHD) is the most frequently diagnosed neurodevelopmental disorder worldwide. Affected individuals present with hyperactivity, inattention, and cognitive deficits and display a characteristic paradoxical response to drugs affecting the dopaminergic system. However, the underlying pathophysiology of ADHD and how this relates to dopaminergic transmission remains to be fully understood. Sorcs2-/- mice uniquely recapitulate symptoms reminiscent of ADHD in humans. Here, we show that lack of SorCS2 in mice results in lower sucrose intake, indicating general reward deficits. Using in-vivo recordings, we further find that dopaminergic transmission in the ventral tegmental area (VTA) is shifted towards a more regular firing pattern with marked reductions in the relative occurrence of irregular firing in Sorcs2-/- mice. This was paralleled by abnormal acute behavioral responses to dopamine receptor agonists, suggesting fundamental differences in dopaminergic circuits and indicating a perturbation in the balance between the activities of the postsynaptic dopamine receptor DRD1 and the presynaptic inhibitory autoreceptor DRD2. Interestingly, the hyperactivity and drug response of Sorcs2-/- mice were markedly affected by novelty. Taken together, our findings show how loss of a candidate ADHD-risk gene has marked effects on dopaminergic circuit function and the behavioral response to the environment.

SORCS2 is one of five proteins that constitute the Vps10p-domain receptor family. Members of this family play important roles in cellular processes linked to neuronal survival, differentiation and function. Genetic and functional studies implicate SORCS2 in cognitive function, as well as in neurodegenerative and psychiatric disorders. DNA damage and DNA repair deficits are linked to ageing and neurodegeneration, and transient neuronal DNA double-strand breaks (DSBs) also occur as a result of neuronal activity. Here, we report a novel role for SORCS2 in DSB formation. We show that SorCS2 loss is associated with elevated DSB levels in the mouse dentate gyrus and that knocking out SORCS2 in a human neuronal cell line increased Topoisomerase IIβ-dependent DSB formation and reduced neuronal viability. Neuronal stimulation had no impact on levels of DNA breaks in vitro, suggesting that the observed differences may not be the result of aberrant neuronal activity in these cells. Our findings are consistent with studies linking the VPS10 receptors and DNA damage to neurodegenerative conditions.

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.

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

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.

PCSK9 induces lysosomal degradation of the low-density lipoprotein (LDL) receptor (LDLR) in the liver, hereby preventing removal of LDL cholesterol from the circulation. Accordingly, PCSK9 inhibitory antibodies and siRNA potently reduce LDL cholesterol to unprecedented low levels and are approved for treatment of hypercholesterolemia. In addition, PCSK9 inactivation alters the levels of several other circulating lipid classes and species. Brain function is critically influenced by cholesterol and lipid composition. However, it remains unclear how the brain is affected long-term by the reduction in circulating lipids as achieved with potent lipid lowering therapeutics such as PCSK9 inhibitors. Furthermore, it is unknown if locally expressed PCSK9 affects neuronal circuits through regulation of receptor levels. We have studied the effect of lifelong low peripheral cholesterol levels on brain lipid composition and behavior in adult PCSK9 KO mice. In addition, we studied the effect of PCSK9 on neurons in culture and in vivo in the developing cerebral cortex. We found that PCSK9 reduced LDLR and neurite complexity in cultured neurons, but neither PCSK9 KO nor overexpression affected cortical development in vivo. Interestingly, PCSK9 deficiency resulted in changes of several lipid classes in the adult cortex and cerebellum. Despite the observed changes, PCSK9 KO mice had unchanged behavior compared to WT controls. In conclusion, our findings demonstrate that altered PCSK9 levels do not compromise brain development or function in mice, and are in line with clinical trials showing that PCSK9 inhibitors have no adverse effects on cognitive function.

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.

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.