Genetics and precision genomics approaches to pulmonary hypertension.

Considerable progress has been made in the genomics of pulmonary arterial hypertension (PAH) since the 6th World Symposium on Pulmonary Hypertension, with the identification of rare variants in several novel genes, as well as common variants that confer a modest increase in PAH risk. Gene and variant curation by an expert panel now provides a robust framework for knowing which genes to test and how to interpret variants in clinical practice. We recommend that genetic testing be offered to specific subgroups of symptomatic patients with PAH, and to children with certain types of group 3 pulmonary hypertension (PH). Testing of asymptomatic family members and the use of genetics in reproductive decision-making require the involvement of genetics experts. Large cohorts of PAH patients with biospecimens now exist and extension to non-group 1 PH has begun. However, these cohorts are largely of European origin; greater diversity will be essential to characterise the full extent of genomic variation contributing to PH risk and treatment responses. Other types of omics data are also being incorporated. Furthermore, to advance gene- and pathway-specific care and targeted therapies, gene-specific registries will be essential to support patients and their families and to lay the foundation for genetically informed clinical trials. This will require international outreach and collaboration between patients/families, clinicians and researchers. Ultimately, harmonisation of patient-derived biospecimens, clinical and omic information, and analytic approaches will advance the field.

Integrative Multiomics in the Lung Reveals a Protective Role of Asporin in Pulmonary Arterial Hypertension.

BACKGROUND: Integrative multiomics can elucidate pulmonary arterial hypertension (PAH) pathobiology, but procuring human PAH lung samples is rare. METHODS: We leveraged transcriptomic profiling and deep phenotyping of the largest multicenter PAH lung biobank to date (96 disease and 52 control) by integration with clinicopathologic data, genome-wide association studies, Bayesian regulatory networks, single-cell transcriptomics, and pharmacotranscriptomics. RESULTS: We identified 2 potentially protective gene network modules associated with vascular cells, and we validated ASPN, coding for asporin, as a key hub gene that is upregulated as a compensatory response to counteract PAH. We found that asporin is upregulated in lungs and plasma of multiple independent PAH cohorts and correlates with reduced PAH severity. We show that asporin inhibits proliferation and transforming growth factor-ß/phosphorylated SMAD2/3 signaling in pulmonary artery smooth muscle cells from PAH lungs. We demonstrate in Sugen-hypoxia rats that ASPN knockdown exacerbated PAH and recombinant asporin attenuated PAH. CONCLUSIONS: Our integrative systems biology approach to dissect the PAH lung transcriptome uncovered asporin as a novel protective target with therapeutic potential in PAH.

Shared and Distinct Genomics of Chronic Thromboembolic Pulmonary Hypertension and Pulmonary Embolism.

Rationale: Chronic thromboembolic pulmonary hypertension involves the formation and nonresolution of thrombus, dysregulated inflammation, angiogenesis, and the development of a small-vessel vasculopathy. Objectives: We aimed to establish the genetic basis of chronic thromboembolic pulmonary hypertension to gain insight into its pathophysiological contributors. Methods: We conducted a genome-wide association study on 1,907 European cases and 10,363 European control subjects. We coanalyzed our results with existing results from genome-wide association studies on deep vein thrombosis, pulmonary embolism, and idiopathic pulmonary arterial hypertension. Measurements and Main Results: Our primary association study revealed genetic associations at the ABO, FGG, F11, MYH7B, and HLA-DRA loci. Through our coanalysis, we demonstrate further associations with chronic thromboembolic pulmonary hypertension at the F2, TSPAN15, SLC44A2, and F5 loci but find no statistically significant associations shared with idiopathic pulmonary arterial hypertension. Conclusions: Chronic thromboembolic pulmonary hypertension is a partially heritable polygenic disease, with related though distinct genetic associations with pulmonary embolism and deep vein thrombosis.

Allele-specific control of rodent and human lncRNA KMT2E-AS1 promotes hypoxic endothelial pathology in pulmonary hypertension.

Hypoxic reprogramming of vasculature relies on genetic, epigenetic, and metabolic circuitry, but the control points are unknown. In pulmonary arterial hypertension (PAH), a disease driven by hypoxia inducible factor (HIF)-dependent vascular dysfunction, HIF-2α promoted expression of neighboring genes, long noncoding RNA (lncRNA) histone lysine N-methyltransferase 2E-antisense 1 (KMT2E-AS1) and histone lysine N-methyltransferase 2E (KMT2E). KMT2E-AS1 stabilized KMT2E protein to increase epigenetic histone 3 lysine 4 trimethylation (H3K4me3), driving HIF-2α-dependent metabolic and pathogenic endothelial activity. This lncRNA axis also increased HIF-2α expression across epigenetic, transcriptional, and posttranscriptional contexts, thus promoting a positive feedback loop to further augment HIF-2α activity. We identified a genetic association between rs73184087, a single-nucleotide variant (SNV) within a KMT2E intron, and disease risk in PAH discovery and replication patient cohorts and in a global meta-analysis. This SNV displayed allele (G)-specific association with HIF-2α, engaged in long-range chromatin interactions, and induced the lncRNA-KMT2E tandem in hypoxic (G/G) cells. In vivo, KMT2E-AS1 deficiency protected against PAH in mice, as did pharmacologic inhibition of histone methylation in rats. Conversely, forced lncRNA expression promoted more severe PH. Thus, the KMT2E-AS1/KMT2E pair orchestrates across convergent multi-ome landscapes to mediate HIF-2α pathobiology and represents a key clinical target in pulmonary hypertension.

Integrative Multiomics to Dissect the Lung Transcriptional Landscape of Pulmonary Arterial Hypertension.

Pulmonary arterial hypertension (PAH) remains an incurable and often fatal disease despite currently available therapies. Multiomics systems biology analysis can shed new light on PAH pathobiology and inform translational research efforts. Using RNA sequencing on the largest PAH lung biobank to date (96 disease and 52 control), we aim to identify gene co-expression network modules associated with PAH and potential therapeutic targets. Co-expression network analysis was performed to identify modules of co-expressed genes which were then assessed for and prioritized by importance in PAH, regulatory role, and therapeutic potential via integration with clinicopathologic data, human genome-wide association studies (GWAS) of PAH, lung Bayesian regulatory networks, single-cell RNA-sequencing data, and pharmacotranscriptomic profiles. We identified a co-expression module of 266 genes, called the pink module, which may be a response to the underlying disease process to counteract disease progression in PAH. This module was associated not only with PAH severity such as increased PVR and intimal thickness, but also with compensated PAH such as lower number of hospitalizations, WHO functional class and NT-proBNP. GWAS integration demonstrated the pink module is enriched for PAH-associated genetic variation in multiple cohorts. Regulatory network analysis revealed that BMPR2 regulates the main target of FDA-approved riociguat, GUCY1A2, in the pink module. Analysis of pathway enrichment and pink hub genes (i.e. ANTXR1 and SFRP4) suggests the pink module inhibits Wnt signaling and epithelial-mesenchymal transition. Cell type deconvolution showed the pink module correlates with higher vascular cell fractions (i.e. myofibroblasts). A pharmacotranscriptomic screen discovered ubiquitin-specific peptidases (USPs) as potential therapeutic targets to mimic the pink module signature. Our multiomics integrative study uncovered a novel gene subnetwork associated with clinicopathologic severity, genetic risk, specific vascular cell types, and new therapeutic targets in PAH. Future studies are warranted to investigate the role and therapeutic potential of the pink module and targeting USPs in PAH.

Genetic counselling and testing in pulmonary arterial hypertension: a consensus statement on behalf of the International Consortium for Genetic Studies in PAH.

Pulmonary arterial hypertension (PAH) is a rare disease that can be caused by (likely) pathogenic germline genomic variants. In addition to the most prevalent disease gene, BMPR2 (bone morphogenetic protein receptor 2), several genes, some belonging to distinct functional classes, are also now known to predispose to the development of PAH. As a consequence, specialist and non-specialist clinicians and healthcare professionals are increasingly faced with a range of questions regarding the need for, approaches to and benefits/risks of genetic testing for PAH patients and/or related family members. We provide a consensus-based approach to recommendations for genetic counselling and assessment of current best practice for disease gene testing. We provide a framework and the type of information to be provided to patients and relatives through the process of genetic counselling, and describe the presently known disease causal genes to be analysed. Benefits of including molecular genetic testing within the management protocol of patients with PAH include the identification of individuals misclassified by other diagnostic approaches, the optimisation of phenotypic characterisation for aggregation of outcome data, including in clinical trials, and importantly through cascade screening, the detection of healthy causal variant carriers, to whom regular assessment should be offered.

Transferability of genetic loci and polygenic scores for cardiometabolic traits in British Pakistani and Bangladeshi individuals.

Individuals with South Asian ancestry have a higher risk of heart disease than other groups but have been largely excluded from genetic research. Using data from 22,000 British Pakistani and Bangladeshi individuals with linked electronic health records from the Genes & Health cohort, we conducted genome-wide association studies of coronary artery disease and its key risk factors. Using power-adjusted transferability ratios, we found evidence for transferability for the majority of cardiometabolic loci powered to replicate. The performance of polygenic scores was high for lipids and blood pressure, but lower for BMI and coronary artery disease. Adding a polygenic score for coronary artery disease to clinical risk factors showed significant improvement in reclassification. In Mendelian randomisation using transferable loci as instruments, our findings were consistent with results in European-ancestry individuals. Taken together, trait-specific transferability of trait loci between populations is an important consideration with implications for risk prediction and causal inference.

Integrating polygenic risk scores in the prediction of type 2 diabetes risk and subtypes in British Pakistanis and Bangladeshis: A population-based cohort study.

BACKGROUND: Type 2 diabetes (T2D) is highly prevalent in British South Asians, yet they are underrepresented in research. Genes & Health (G&H) is a large, population study of British Pakistanis and Bangladeshis (BPB) comprising genomic and routine health data. We assessed the extent to which genetic risk for T2D is shared between BPB and European populations (EUR). We then investigated whether the integration of a polygenic risk score (PRS) for T2D with an existing risk tool (QDiabetes) could improve prediction of incident disease and the characterisation of disease subtypes. METHODS AND FINDINGS: In this observational cohort study, we assessed whether common genetic loci associated with T2D in EUR individuals were replicated in 22,490 BPB individuals in G&H. We replicated fewer loci in G&H (n = 76/338, 22%) than would be expected given power if all EUR-ascertained loci were transferable (n = 101, 30%; p = 0.001). Of the 27 transferable loci that were powered to interrogate this, only 9 showed evidence of shared causal variants. We constructed a T2D PRS and combined it with a clinical risk instrument (QDiabetes) in a novel, integrated risk tool (IRT) to assess risk of incident diabetes. To assess model performance, we compared categorical net reclassification index (NRI) versus QDiabetes alone. In 13,648 patients free from T2D followed up for 10 years, NRI was 3.2% for IRT versus QDiabetes (95% confidence interval (CI): 2.0% to 4.4%). IRT performed best in reclassification of individuals aged less than 40 years deemed low risk by QDiabetes alone (NRI 5.6%, 95% CI 3.6% to 7.6%), who tended to be free from comorbidities and slim. After adjustment for QDiabetes score, PRS was independently associated with progression to T2D after gestational diabetes (hazard ratio (HR) per SD of PRS 1.23, 95% CI 1.05 to 1.42, p = 0.028). Using cluster analysis of clinical features at diabetes diagnosis, we replicated previously reported disease subgroups, including Mild Age-Related, Mild Obesity-related, and Insulin-Resistant Diabetes, and showed that PRS distribution differs between subgroups (p = 0.002). Integrating PRS in this cluster analysis revealed a Probable Severe Insulin Deficient Diabetes (pSIDD) subgroup, despite the absence of clinical measures of insulin secretion or resistance. We also observed differences in rates of progression to micro- and macrovascular complications between subgroups after adjustment for confounders. Study limitations include the absence of an external replication cohort and the potential biases arising from missing or incorrect routine health data. CONCLUSIONS: Our analysis of the transferability of T2D loci between EUR and BPB indicates the need for larger, multiancestry studies to better characterise the genetic contribution to disease and its varied aetiology. We show that a T2D PRS optimised for this high-risk BPB population has potential clinical application in BPB, improving the identification of T2D risk (especially in the young) on top of an established clinical risk algorithm and aiding identification of subgroups at diagnosis, which may help future efforts to stratify care and treatment of the disease.

Mendelian randomisation and experimental medicine approaches to interleukin-6 as a drug target in pulmonary arterial hypertension.

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-1) 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±sd age 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 lead IL6R variant 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.

Genome-Wide Association Study Identifies Risk Loci for Cluster Headache.

OBJECTIVE: This study was undertaken to identify susceptibility loci for cluster headache and obtain insights into relevant disease pathways. METHODS: We carried out a genome-wide association study, where 852 UK and 591 Swedish cluster headache cases were compared with 5,614 and 1,134 controls, respectively. Following quality control and imputation, single variant association testing was conducted using a logistic mixed model for each cohort. The 2 cohorts were subsequently combined in a merged analysis. Downstream analyses, such as gene-set enrichment, functional variant annotation, prediction and pathway analyses, were performed. RESULTS: Initial independent analysis identified 2 replicable cluster headache susceptibility loci on chromosome 2. A merged analysis identified an additional locus on chromosome 1 and confirmed a locus significant in the UK analysis on chromosome 6, which overlaps with a previously known migraine locus. The lead single nucleotide polymorphisms were rs113658130 (p = 1.92 × 10-17 , odds ratio [OR] = 1.51, 95% confidence interval [CI] = 1.37-1.66) and rs4519530 (p = 6.98 × 10-17 , OR = 1.47, 95% CI = 1.34-1.61) on chromosome 2, rs12121134 on chromosome 1 (p = 1.66 × 10-8 , OR = 1.36, 95% CI = 1.22-1.52), and rs11153082 (p = 1.85 × 10-8 , OR = 1.30, 95% CI = 1.19-1.42) on chromosome 6. Downstream analyses implicated immunological processes in the pathogenesis of cluster headache. INTERPRETATION: We identified and replicated several genome-wide significant associations supporting a genetic predisposition in cluster headache in a genome-wide association study involving 1,443 cases. Replication in larger independent cohorts combined with comprehensive phenotyping, in relation to, for example, treatment response and cluster headache subtypes, could provide unprecedented insights into genotype-phenotype correlations and the pathophysiological pathways underlying cluster headache. ANN NEUROL 2021;90:193-202.

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.

Whole Exome Sequence Analysis Provides Novel Insights into the Genetic Framework of Childhood-Onset Pulmonary Arterial Hypertension.

Pulmonary arterial hypertension (PAH) describes a rare, progressive vascular disease caused by the obstruction of pulmonary arterioles, typically resulting in right heart failure. Whilst PAH most often manifests in adulthood, paediatric disease is considered to be a distinct entity with increased morbidity and often an unexplained resistance to current therapies. Recent genetic studies have substantially increased our understanding of PAH pathogenesis, providing opportunities for molecular diagnosis and presymptomatic genetic testing in families. However, the genetic architecture of childhood-onset PAH remains relatively poorly characterised. We sought to investigate a previously unsolved paediatric cohort (n = 18) using whole exome sequencing to improve the molecular diagnosis of childhood-onset PAH. Through a targeted investigation of 26 candidate genes, we applied a rigorous variant filtering methodology to enrich for rare, likely pathogenic variants. This analysis led to the detection of novel PAH risk alleles in five genes, including the first identification of a heterozygous ATP13A3 mutation in childhood-onset disease. In addition, we provide the first independent validation of BMP10 and PDGFD as genetic risk factors for PAH. These data provide a molecular diagnosis in 28% of paediatric cases, reflecting the increased genetic burden in childhood-onset disease and highlighting the importance of next-generation sequencing approaches to diagnostic surveillance.

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.

Evaluating drug targets through human loss-of-function genetic variation.

Naturally occurring human genetic variants that are predicted to inactivate protein-coding genes provide an in vivo model of human gene inactivation that complements knockout studies in cells and model organisms. Here we report three key findings regarding the assessment of candidate drug targets using human loss-of-function variants. First, even essential genes, in which loss-of-function variants are not tolerated, can be highly successful as targets of inhibitory drugs. Second, in most genes, loss-of-function variants are sufficiently rare that genotype-based ascertainment of homozygous or compound heterozygous 'knockout' humans will await sample sizes that are approximately 1,000 times those presently available, unless recruitment focuses on consanguineous individuals. Third, automated variant annotation and filtering are powerful, but manual curation remains crucial for removing artefacts, and is a prerequisite for recall-by-genotype efforts. Our results provide a roadmap for human knockout studies and should guide the interpretation of loss-of-function variants in drug development.

Correction: A restricted spectrum of missense KMT2D variants cause a multiple malformations disorder distinct from Kabuki syndrome.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Characterising a healthy adult with a rare HAO1 knockout to support a therapeutic strategy for primary hyperoxaluria.

By sequencing autozygous human populations, we identified a healthy adult woman with lifelong complete knockout of HAO1 (expected ~1 in 30 million outbred people). HAO1 (glycolate oxidase) silencing is the mechanism of lumasiran, an investigational RNA interference therapeutic for primary hyperoxaluria type 1. Her plasma glycolate levels were 12 times, and urinary glycolate 6 times, the upper limit of normal observed in healthy reference individuals (n = 67). Plasma metabolomics and lipidomics (1871 biochemicals) revealed 18 markedly elevated biochemicals (>5 sd outliers versus n = 25 controls) suggesting additional HAO1 effects. Comparison with lumasiran preclinical and clinical trial data suggested she has <2% residual glycolate oxidase activity. Cell line p.Leu333SerfsTer4 expression showed markedly reduced HAO1 protein levels and cellular protein mis-localisation. In this woman, lifelong HAO1 knockout is safe and without clinical phenotype, de-risking a therapeutic approach and informing therapeutic mechanisms. Unlocking evidence from the diversity of human genetic variation can facilitate drug development.

A restricted spectrum of missense KMT2D variants cause a multiple malformations disorder distinct from Kabuki syndrome.

PURPOSE: To investigate if specific exon 38 or 39 KMT2D missense variants (MVs) cause a condition distinct from Kabuki syndrome type 1 (KS1). METHODS: Multiple individuals, with MVs in exons 38 or 39 of KMT2D that encode a highly conserved region of 54 amino acids flanked by Val3527 and Lys3583, were identified and phenotyped. Functional tests were performed to study their pathogenicity and understand the disease mechanism. RESULTS: The consistent clinical features of the affected individuals, from seven unrelated families, included choanal atresia, athelia or hypoplastic nipples, branchial sinus abnormalities, neck pits, lacrimal duct anomalies, hearing loss, external ear malformations, and thyroid abnormalities. None of the individuals had intellectual disability. The frequency of clinical features, objective software-based facial analysis metrics, and genome-wide peripheral blood DNA methylation patterns in these patients were significantly different from that of KS1. Circular dichroism spectroscopy indicated that these MVs perturb KMT2D secondary structure through an increased disordered to ɑ-helical transition. CONCLUSION: KMT2D MVs located in a specific region spanning exons 38 and 39 and affecting highly conserved residues cause a novel multiple malformations syndrome distinct from KS1. Unlike KMT2D haploinsufficiency in KS1, these MVs likely result in disease through a dominant negative mechanism.