Identifying digenic disease genes via machine learning in the Undiagnosed Diseases Network.

Rare diseases affect millions of people worldwide, and discovering their genetic causes is challenging. More than half of the individuals analyzed by the Undiagnosed Diseases Network (UDN) remain undiagnosed. The central hypothesis of this work is that many of these rare genetic disorders are caused by multiple variants in more than one gene. However, given the large number of variants in each individual genome, experimentally evaluating combinations of variants for potential to cause disease is currently infeasible. To address this challenge, we developed the digenic predictor (DiGePred), a random forest classifier for identifying candidate digenic disease gene pairs by features derived from biological networks, genomics, evolutionary history, and functional annotations. We trained the DiGePred classifier by using DIDA, the largest available database of known digenic-disease-causing gene pairs, and several sets of non-digenic gene pairs, including variant pairs derived from unaffected relatives of UDN individuals. DiGePred achieved high precision and recall in cross-validation and on a held-out test set (PR area under the curve > 77%), and we further demonstrate its utility by using digenic pairs from the recent literature. In contrast to other approaches, DiGePred also appropriately controls the number of false positives when applied in realistic clinical settings. Finally, to enable the rapid screening of variant gene pairs for digenic disease potential, we freely provide the predictions of DiGePred on all human gene pairs. Our work enables the discovery of genetic causes for rare non-monogenic diseases by providing a means to rapidly evaluate variant gene pairs for the potential to cause digenic disease.

Probable digenic inheritance of Diamond-Blackfan anemia.

A 26-year-old female proband with a clinical diagnosis and consistent phenotype of Diamond-Blackfan anemia (DBA, OMIM 105650) without an identified genotype was referred to the Undiagnosed Diseases Network. DBA is classically associated with monoallelic variants that have an autosomal-dominant or -recessive mode of inheritance. Intriguingly, her case was solved by a detection of a digenic interaction between non-allelic RPS19 and RPL27 variants. This was confirmed with a machine learning structural model, co-segregation analysis, and RNA sequencing. This is the first report of DBA caused by a digenic effect of two non-allelic variants demonstrated by machine learning structural model. This case suggests that atypical phenotypic presentations of DBA may be caused by digenic inheritance in some individuals. We also conclude that a machine learning structural model can be useful in detecting digenic models of possible interactions between products encoded by alleles of different genes inherited from non-affected carrier parents that can result in DBA with an unrealized 25% recurrence risk.

Undiagnosed Disease Network collaborative approach in diagnosing rare disease in a patient with a mosaic CACNA1D variant.

The Undiagnosed Disease Network (UDN) is comprised of clinical and research experts collaborating to diagnose rare disease. The UDN is funded by the National Institutes of Health and includes 12 different clinical sites (About Us, 2022). Here we highlight the success of collaborative efforts within the UDN Clinical Site at Vanderbilt University Medical Center (VUMC) in utilizing a cohort of experts in bioinformatics, structural biology, and genetics specialists in diagnosing rare disease. Our UDN team identified a de novo mosaic CACNA1D variant c.2299T>C in a 5-year-old female with a history of global developmental delay, dystonia, dyskinesis, and seizures. Using a collaborative multidisciplinary approach, our VUMC UDN team diagnosed the participant with Primary Aldosteronism, Seizures, and Neurologic abnormalities (PASNA) OMIM: 615474 due to a rare mosaic CACNA1D variant (O'Neill, 2013). Interestingly, this patient was mosaic, a phenotypic trait previously unreported in PASNA cases. This report highlights the importance of a multidisciplinary approach in diagnosing rare disease.

The Genetic Landscape of Familial Pulmonary Fibrosis.

Rationale and Objectives: Up to 20% of idiopathic interstitial lung disease is familial, referred to as familial pulmonary fibrosis (FPF). An integrated analysis of FPF genetic risk was performed by comprehensively evaluating for genetic rare variants (RVs) in a large cohort of FPF kindreds. Methods: Whole-exome sequencing and/or candidate gene sequencing from affected individuals in 569 FPF kindreds was performed, followed by cosegregation analysis in large kindreds, gene burden analysis, gene-based risk scoring, cell-type enrichment analysis, and coexpression network construction. Measurements and Main Results: It was found that 14.9-23.4% of genetic risk in kindreds could be explained by RVs in genes previously linked to FPF, predominantly telomere-related genes. New candidate genes were identified in a small number of families-including SYDE1, SERPINB8, GPR87, and NETO1-and tools were developed for evaluation and prioritization of RV-containing genes across kindreds. Several pathways were enriched for RV-containing genes in FPF, including focal adhesion and mitochondrial complex I assembly. By combining single-cell transcriptomics with prioritized candidate genes, expression of RV-containing genes was discovered to be enriched in smooth muscle cells, type II alveolar epithelial cells, and endothelial cells. Conclusions: In the most comprehensive FPF genetic study to date, the prevalence of RVs in known FPF-related genes was defined, and new candidate genes and pathways relevant to FPF were identified. However, new RV-containing genes shared across multiple kindreds were not identified, thereby suggesting that heterogeneous genetic variants involving a variety of genes and pathways mediate genetic risk in most FPF kindreds.

A medical odyssey of a 72-year-old man with Charcot-Marie-Tooth disease type 2 newly diagnosed with biallelic variants in SORD gene causing sorbitol dehydrogenase deficiency.

A 72-year-old man was referred to the Undiagnosed Diseases Network (UDN) because of gradual progressive weakness in both lower extremities for the past 45 years. He was initially diagnosed as having Charcot-Marie-Tooth disease type 2 (CMT2) without a defined molecular genetic cause. Exome sequencing (ES) failed to detect deleterious neuromuscular variants. Very recently, biallelic variants in sorbitol dehydrogenase (SORD) were discovered to be a novel cause of inherited neuropathies including CMT2 or distal hereditary motor neuropathy (dHMN) referred to as Sorbitol Dehydrogenase Deficiency with Peripheral Neuropathy (SORDD, OMIM 618912). The most common variant identified was c.757delG; p.A253Qfs*27. Through the Vanderbilt UDN clinical site, this patient was formally diagnosed with SORDD after the identification of homozygosity for the above SORD frameshift through UDN Genome Sequencing (GS). His medical odyssey was solved by GS and detection of extremely high levels of sorbitol. The diagnosis provided him the opportunity to receive potential treatment with an investigational drug in a clinical trial for SORDD. We suggest that similar studies be considered in other individuals thought to possibly have CMT2 or dHMN.

The contribution of mosaicism to genetic diseases and de novo pathogenic variants.

The contribution of mosaicism to diagnosed genetic disease and presumed de novo variants (DNV) is under investigated. We determined the contribution of mosaic genetic disease (MGD) and diagnosed parental mosaicism (PM) in parents of offspring with reported DNV (in the same variant) in the (1) Undiagnosed Diseases Network (UDN) (N = 1946) and (2) in 12,472 individuals electronic health records (EHR) who underwent genetic testing at an academic medical center. In the UDN, we found 4.51% of diagnosed probands had MGD, and 2.86% of parents of those with DNV exhibited PM. In the EHR, we found 6.03% and 2.99% and (of diagnosed probands) had MGD detected on chromosomal microarray and exome/genome sequencing, respectively. We found 2.34% (of those with a presumed pathogenic DNV) had a parent with PM for the variant. We detected mosaicism (regardless of pathogenicity) in 4.49% of genetic tests performed. We found a broad phenotypic spectrum of MGD with previously unknown phenotypic phenomena. MGD is highly heterogeneous and provides a significant contribution to genetic diseases. Further work is required to improve the diagnosis of MGD and investigate how PM contributes to DNV risk.

Commonalities across computational workflows for uncovering explanatory variants in undiagnosed cases.

PURPOSE: Genomic sequencing has become an increasingly powerful and relevant tool to be leveraged for the discovery of genetic aberrations underlying rare, Mendelian conditions. Although the computational tools incorporated into diagnostic workflows for this task are continually evolving and improving, we nevertheless sought to investigate commonalities across sequencing processing workflows to reveal consensus and standard practice tools and highlight exploratory analyses where technical and theoretical method improvements would be most impactful. METHODS: We collected details regarding the computational approaches used by a genetic testing laboratory and 11 clinical research sites in the United States participating in the Undiagnosed Diseases Network via meetings with bioinformaticians, online survey forms, and analyses of internal protocols. RESULTS: We found that tools for processing genomic sequencing data can be grouped into four distinct categories. Whereas well-established practices exist for initial variant calling and quality control steps, there is substantial divergence across sites in later stages for variant prioritization and multimodal data integration, demonstrating a diversity of approaches for solving the most mysterious undiagnosed cases. CONCLUSION: The largest differences across diagnostic workflows suggest that advances in structural variant detection, noncoding variant interpretation, and integration of additional biomedical data may be especially promising for solving chronically undiagnosed cases.

Clinical sites of the Undiagnosed Diseases Network: unique contributions to genomic medicine and science.

PURPOSE: The NIH Undiagnosed Diseases Network (UDN) evaluates participants with disorders that have defied diagnosis, applying personalized clinical and genomic evaluations and innovative research. The clinical sites of the UDN are essential to advancing the UDN mission; this study assesses their contributions relative to standard clinical practices. METHODS: We analyzed retrospective data from four UDN clinical sites, from July 2015 to September 2019, for diagnoses, new disease gene discoveries and the underlying investigative methods. RESULTS: Of 791 evaluated individuals, 231 received 240 diagnoses and 17 new disease-gene associations were recognized. Straightforward diagnoses on UDN exome and genome sequencing occurred in 35% (84/240). We considered these tractable in standard clinical practice, although genome sequencing is not yet widely available clinically. The majority (156/240, 65%) required additional UDN-driven investigations, including 90 diagnoses that occurred after prior nondiagnostic exome sequencing and 45 diagnoses (19%) that were nongenetic. The UDN-driven investigations included complementary/supplementary phenotyping, innovative analyses of genomic variants, and collaborative science for functional assays and animal modeling. CONCLUSION: Investigations driven by the clinical sites identified diagnostic and research paradigms that surpass standard diagnostic processes. The new diagnoses, disease gene discoveries, and delineation of novel disorders represent a model for genomic medicine and science.

Development and Progression of Radiologic Abnormalities in Individuals at Risk for Familial Interstitial Lung Disease.

Rationale: The preclinical natural history of progressive lung fibrosis is poorly understood.Objectives: Our goals were to identify risk factors for interstitial lung abnormalities (ILA) on high-resolution computed tomography (HRCT) scans and to determine progression toward clinical interstitial lung disease (ILD) among subjects in a longitudinal cohort of self-reported unaffected first-degree relatives of patients with familial interstitial pneumonia.Methods: Enrollment evaluation included a health history and exposure questionnaire and HRCT scans, which were categorized by visual assessment as no ILA, early/mild ILA, or extensive ILA. The study endpoint was met when ILA were extensive or when ILD was diagnosed clinically. Among subjects with adequate study time to complete 5-year follow-up HRCT, the proportion with ILD events (endpoint met or radiographic ILA progression) was calculated.Measurements and Main Results: Among 336 subjects, the mean age was 53.1 (SD, 9.9) years. Those with ILA (early/mild [n = 74] or extensive [n = 3]) were older, were more likely to be ever smokers, had shorter peripheral blood mononuclear cell telomeres, and were more likely to carry the MUC5B risk allele. Self-reported occupational or environmental exposures, including aluminum smelting, lead, birds, and mold, were independently associated with ILA. Among 129 subjects with sufficient study time, 25 (19.4%) had an ILD event by 5 years after enrollment; of these, 12 met the study endpoint and another 13 had radiologic progression of ILA. ILD events were more common among those with early/mild ILA at enrollment (63.3% vs. 6.1%; P < 0.0001).Conclusions: Rare and common environmental exposures are independent risk factors for radiologic abnormalities. In 5 years, progression of ILA occurred in most individuals with early ILA detected at enrollment.

HACER: an atlas of human active enhancers to interpret regulatory variants.

Recent studies have shown that disease-susceptibility variants frequently lie in cell-type-specific enhancer elements. To identify, interpret, and prioritize such risk variants, we must identify the enhancers active in disease-relevant cell types, their upstream transcription factor (TF) binding, and their downstream target genes. To address this need, we built HACER (http://bioinfo.vanderbilt.edu/AE/HACER/), an atlas of Human ACtive Enhancers to interpret Regulatory variants. The HACER atlas catalogues and annotates in-vivo transcribed cell-type-specific enhancers, as well as placing enhancers within transcriptional regulatory networks by integrating ENCODE TF ChIP-Seq and predicted/validated chromatin interaction data. We demonstrate the utility of HACER in (i) offering a mechanistic hypothesis to explain the association of SNP rs614367 with ER-positive breast cancer risk, (ii) exploring tumor-specific enhancers in selective MYC dysregulation and (iii) prioritizing/annotating non-coding regulatory regions targeting CCND1. HACER provides a valuable resource for studies of GWAS, non-coding variants, and enhancer-mediated regulation.

Limitations of exome sequencing in detecting rare and undiagnosed diseases.

While exome sequencing (ES) is commonly the final diagnostic step in clinical genetics, it may miss diagnoses. To clarify the limitations of ES, we investigated the diagnostic yield of genetic tests beyond ES in our Undiagnosed Diseases Network (UDN) participants. We reviewed the yield of additional genetic testing including genome sequencing (GS), copy number variant (CNV), noncoding variant (NCV), repeat expansion (RE), or methylation testing in UDN cases with nondiagnostic ES results. Overall, 36/54 (67%) of total diagnoses were based on clinical findings and coding variants found by ES and 3/54 (6%) were based on clinical findings only. The remaining 15/54 (28%) required testing beyond ES. Of these, 7/15 (47%) had NCV, 6/15 (40%) CNV, and 2/15 (13%) had a RE or a DNA methylation disorder. Thus 18/54 (33%) of diagnoses were not solved exclusively by ES. Several methods were needed to detect and/or confirm the functional effects of the variants missed by ES, and in some cases by GS. These results indicate that tests to detect elusive variants should be considered after nondiagnostic preliminary steps. Further studies are needed to determine the cost-effectiveness of tests beyond ES that provide diagnoses and insights to possible treatment.

Resequencing Study Confirms That Host Defense and Cell Senescence Gene Variants Contribute to the Risk of Idiopathic Pulmonary Fibrosis.

Rationale: Several common and rare genetic variants have been associated with idiopathic pulmonary fibrosis, a progressive fibrotic condition that is localized to the lung. Objectives: To develop an integrated understanding of the rare and common variants located in multiple loci that have been reported to contribute to the risk of disease. Methods: We performed deep targeted resequencing (3.69 Mb of DNA) in cases (n = 3,624) and control subjects (n = 4,442) across genes and regions previously associated with disease. We tested for associations between disease and 1) individual common variants via logistic regression and 2) groups of rare variants via sequence kernel association tests. Measurements and Main Results: Statistically significant common variant association signals occurred in all 10 of the regions chosen based on genome-wide association studies. The strongest risk variant is the MUC5B promoter variant rs35705950, with an odds ratio of 5.45 (95% confidence interval, 4.91-6.06) for one copy of the risk allele and 18.68 (95% confidence interval, 13.34-26.17) for two copies of the risk allele (P = 9.60 × 10-295). In addition to identifying for the first time that rare variation in FAM13A is associated with disease, we confirmed the role of rare variation in the TERT and RTEL1 gene regions in the risk of IPF, and found that the FAM13A and TERT regions have independent common and rare variant signals. Conclusions: A limited number of common and rare variants contribute to the risk of idiopathic pulmonary fibrosis in each of the resequencing regions, and these genetic variants focus on biological mechanisms of host defense and cell senescence.

Three-dimensional spatial analysis of missense variants in RTEL1 identifies pathogenic variants in patients with Familial Interstitial Pneumonia.

BACKGROUND: Next-generation sequencing of individuals with genetic diseases often detects candidate rare variants in numerous genes, but determining which are causal remains challenging. We hypothesized that the spatial distribution of missense variants in protein structures contains information about function and pathogenicity that can help prioritize variants of unknown significance (VUS) and elucidate the structural mechanisms leading to disease. RESULTS: To illustrate this approach in a clinical application, we analyzed 13 candidate missense variants in regulator of telomere elongation helicase 1 (RTEL1) identified in patients with Familial Interstitial Pneumonia (FIP). We curated pathogenic and neutral RTEL1 variants from the literature and public databases. We then used homology modeling to construct a 3D structural model of RTEL1 and mapped known variants into this structure. We next developed a pathogenicity prediction algorithm based on proximity to known disease causing and neutral variants and evaluated its performance with leave-one-out cross-validation. We further validated our predictions with segregation analyses, telomere lengths, and mutagenesis data from the homologous XPD protein. Our algorithm for classifying RTEL1 VUS based on spatial proximity to pathogenic and neutral variation accurately distinguished 7 known pathogenic from 29 neutral variants (ROC AUC = 0.85) in the N-terminal domains of RTEL1. Pathogenic proximity scores were also significantly correlated with effects on ATPase activity (Pearson r = -0.65, p = 0.0004) in XPD, a related helicase. Applying the algorithm to 13 VUS identified from sequencing of RTEL1 from patients predicted five out of six disease-segregating VUS to be pathogenic. We provide structural hypotheses regarding how these mutations may disrupt RTEL1 ATPase and helicase function. CONCLUSIONS: Spatial analysis of missense variation accurately classified candidate VUS in RTEL1 and suggests how such variants cause disease. Incorporating spatial proximity analyses into other pathogenicity prediction tools may improve accuracy for other genes and genetic diseases.

Phenotypic heterogeneity of ZMPSTE24 deficiency.

A 4-year-old girl was referred to the Undiagnosed Diseases Network with a history of short stature, thin and translucent skin, macrocephaly, small hands, and camptodactyly. She had been diagnosed with possible Hallerman-Streiff syndrome. Her evaluation showed that she was mosaic for uniparental isodisomy of chromosome 1, which harbored a pathogenic c.1077dupT variant in ZMPSTE24 which predicts p.(Leu362fsX18). ZMPSTE24 is a zinc metalloproteinase that is involved in processing farnesylated proteins and pathogenic ZMPSTE24 variants cause accumulation of abnormal farnesylated forms of prelamin A. This, in turn, causes a spectrum of disease severity which is based on enzyme activity. The current patient has an intermediate form, which is a genocopy of severe Progeria.

Clinical heterogeneity of mitochondrial NAD kinase deficiency caused by a NADK2 start loss variant.

Mitochondrial NAD kinase deficiency (NADK2D, OMIM #615787) is a rare autosomal recessive disorder of NADPH biosynthesis that can cause hyperlysinemia and dienoyl-CoA reductase deficiency (DECRD, OMIM #616034). NADK2 deficiency has been reported in only three unrelated patients. Two had severe, unremitting disease; one died at 4 months and the other at 5 years of age. The third was a 10 year old female with CNS anomalies, ataxia, and incoordination. In two cases mutations in NADK2 have been demonstrated. Here, we report the fourth known case, a 15 year old female with normal intelligence and a mild clinical and biochemical phenotype presumably without DECRD. Her clinical symptoms, which are now stable, became evident at the age of 9 with the onset of decreased visual acuity, bilateral optic atrophy, nystagmus, episodic lower extremity weakness, peripheral neuropathy, and gait abnormalities. Plasma amino acid levels were within normal limits except for mean lysine and proline levels that were 3.7 and 2.5 times the upper limits of normal. Whole exome sequencing (WES) revealed homozygosity for a g.36241900 A>G p. Met1Val start loss mutation in the primary NADK2 transcript (NM_001085411.1) encoding the 442 amino acid isoform. This presumed hypomorphic mutation has not been previously reported and is absent from the v1000GP, EVS, and ExAC databases. Our patient's normal intelligence and stable disease expands the clinical heterogeneity and the prognosis associated with NADK2 deficiency. Our findings also clarify the mechanism underlying NADK2 deficiency and suggest that this disease should be ruled out in cases of hyperlysinemia, especially those with visual loss, and neurological phenotypes.

Estrogen Metabolite 16α-Hydroxyestrone Exacerbates Bone Morphogenetic Protein Receptor Type II-Associated Pulmonary Arterial Hypertension Through MicroRNA-29-Mediated Modulation of Cellular Metabolism.

BACKGROUND: Pulmonary arterial hypertension (PAH) is a proliferative disease of the pulmonary vasculature that preferentially affects women. Estrogens such as the metabolite 16α-hydroxyestrone (16αOHE) may contribute to PAH pathogenesis, and alterations in cellular energy metabolism associate with PAH. We hypothesized that 16αOHE promotes heritable PAH (HPAH) via microRNA-29 (miR-29) family upregulation and that antagonism of miR-29 would attenuate pulmonary hypertension in transgenic mouse models of Bmpr2 mutation. METHODS AND RESULTS: MicroRNA array profiling of human lung tissue found elevation of microRNAs associated with energy metabolism, including the miR-29 family, among HPAH patients. miR-29 expression was 2-fold higher in Bmpr2 mutant mice lungs at baseline compared with controls and 4 to 8-fold higher in Bmpr2 mice exposed to 16αOHE 1.25 µg/h for 4 weeks. Blot analyses of Bmpr2 mouse lung protein showed significant reductions in peroxisome proliferator-activated receptor-γ and CD36 in those mice exposed to 16αOHE and protein derived from HPAH lungs compared with controls. Bmpr2 mice treated with anti-miR-29 (20-mg/kg injections for 6 weeks) had improvements in hemodynamic profile, histology, and markers of dysregulated energy metabolism compared with controls. Pulmonary artery smooth muscle cells derived from Bmpr2 murine lungs demonstrated mitochondrial abnormalities, which improved with anti-miR-29 transfection in vitro; endothelial-like cells derived from HPAH patient induced pluripotent stem cell lines were similar and improved with anti-miR-29 treatment. CONCLUSIONS: 16αOHE promotes the development of HPAH via upregulation of miR-29, which alters molecular and functional indexes of energy metabolism. Antagonism of miR-29 improves in vivo and in vitro features of HPAH and reveals a possible novel therapeutic target.

Bone Marrow-derived Cells Contribute to the Pathogenesis of Pulmonary Arterial Hypertension.

RATIONALE: Pulmonary arterial hypertension (PAH) is a progressive lung disease of the pulmonary microvasculature. Studies suggest that bone marrow (BM)-derived circulating cells may play an important role in its pathogenesis. OBJECTIVES: We used a genetic model of PAH, the Bmpr2 mutant mouse, to study the role of BM-derived circulating cells in its pathogenesis. METHODS: Recipient mice, either Bmpr2(R899X) mutant or controls, were lethally irradiated and transplanted with either control or Bmpr2(R899X) BM cells. Donor cells were traced in female recipient mice by Y chromosome painting. Molecular and function insights were provided by expression and cytokine arrays combined with flow cytometry, colony-forming assays, and competitive transplant assays. MEASUREMENTS AND MAIN RESULTS: We found that mutant BM cells caused PAH with remodeling and inflammation when transplanted into control mice, whereas control BM cells had a protective effect against the development of disease, when transplanted into mutant mice. Donor BM-derived cells were present in the lungs of recipient mice. Functional and molecular analysis identified mutant BM cell dysfunction suggestive of a PAH phenotype soon after activation of the transgene and long before the development of lung pathology. CONCLUSIONS: Our data show that BM cells played a key role in PAH pathogenesis and that the transplanted BM cells were able to drive the lung phenotype in a myeloablative transplant model. Furthermore, the specific cell types involved were derived from hematopoietic stem cells and exhibit dysfunction long before the development of lung pathology.

Critical Genomic Networks and Vasoreactive Variants in Idiopathic Pulmonary Arterial Hypertension.

RATIONALE: Idiopathic pulmonary arterial hypertension (IPAH) is usually without an identified genetic cause, despite clinical and molecular similarity to bone morphogenetic protein receptor type 2 mutation-associated heritable pulmonary arterial hypertension (PAH). There is phenotypic heterogeneity in IPAH, with a minority of patients showing long-term improvement with calcium channel-blocker therapy. OBJECTIVES: We sought to identify gene variants (GVs) underlying IPAH and determine whether GVs differ in vasodilator-responsive IPAH (VR-PAH) versus vasodilator-nonresponsive IPAH (VN-PAH). METHODS: We performed whole-exome sequencing (WES) on 36 patients with IPAH: 17 with VR-PAH and 19 with VN-PAH. Wnt pathway differences were explored in human lung fibroblasts. MEASUREMENTS AND MAIN RESULTS: We identified 1,369 genes with 1,580 variants unique to IPAH. We used a gene ontology approach to analyze variants and identified overrepresentation of several pathways, including cytoskeletal function and ion binding. By mapping WES data to prior genome-wide association study data, Wnt pathway genes were highlighted. Using the connectivity map to define genetic differences between VR-PAH and VN-PAH, we found enrichment in vascular smooth muscle cell contraction pathways and greater genetic variation in VR-PAH versus VN-PAH. Using human lung fibroblasts, we found increased stimulated Wnt activity in IPAH versus controls. CONCLUSIONS: A pathway-based analysis of WES data in IPAH demonstrated multiple rare GVs that converge on key biological pathways, such as cytoskeletal function and Wnt signaling pathway. Vascular smooth muscle contraction-related genes were enriched in VR-PAH, suggesting a potentially different genetic predisposition for VR-PAH. This pathway-based approach may be applied to next-generation sequencing data in other diseases to uncover the contribution of unexpected or multiple GVs to a phenotype.