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.

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.

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.

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.

Extensive phenotyping of individuals at risk for familial interstitial pneumonia reveals clues to the pathogenesis of interstitial lung disease.

RATIONALE: Asymptomatic relatives of patients with familial interstitial pneumonia (FIP), the inherited form of idiopathic interstitial pneumonia, carry increased risk for developing interstitial lung disease. OBJECTIVES: Studying these at-risk individuals provides a unique opportunity to investigate early stages of FIP pathogenesis and develop predictive models of disease onset. METHODS: Seventy-five asymptomatic first-degree relatives of FIP patients (mean age, 50.8 yr) underwent blood sampling and high-resolution chest computed tomography (HRCT) scanning in an ongoing cohort study; 72 consented to bronchoscopy with bronchoalveolar lavage (BAL) and transbronchial biopsies. Twenty-seven healthy individuals were used as control subjects. MEASUREMENTS AND MAIN RESULTS: Eleven of 75 at-risk subjects (14%) had evidence of interstitial changes by HRCT, whereas 35.2% had abnormalities on transbronchial biopsies. No differences were noted in inflammatory cells in BAL between at-risk individuals and control subjects. At-risk subjects had increased herpesvirus DNA in cell-free BAL and evidence of herpesvirus antigen expression in alveolar epithelial cells (AECs), which correlated with expression of endoplasmic reticulum stress markers in AECs. Peripheral blood mononuclear cell and AEC telomere length were shorter in at-risk individuals than healthy control subjects. The minor allele frequency of the Muc5B rs35705950 promoter polymorphism was increased in at-risk subjects. Levels of several plasma biomarkers differed between at-risk subjects and control subjects, and correlated with abnormal HRCT scans. CONCLUSIONS: Evidence of lung parenchymal remodeling and epithelial dysfunction was identified in asymptomatic individuals at risk for FIP. Together, these findings offer new insights into the early pathogenesis of idiopathic interstitial pneumonia and provide an ongoing opportunity to characterize presymptomatic abnormalities that predict progression to clinical disease.

Rare variants in RTEL1 are associated with familial interstitial pneumonia.

RATIONALE: Up to 20% of cases of idiopathic interstitial pneumonia cluster in families, comprising the syndrome of familial interstitial pneumonia (FIP); however, the genetic basis of FIP remains uncertain in most families. OBJECTIVES: To determine if new disease-causing rare genetic variants could be identified using whole-exome sequencing of affected members from FIP families, providing additional insights into disease pathogenesis. METHODS: Affected subjects from 25 kindreds were selected from an ongoing FIP registry for whole-exome sequencing from genomic DNA. Candidate rare variants were confirmed by Sanger sequencing, and cosegregation analysis was performed in families, followed by additional sequencing of affected individuals from another 163 kindreds. MEASUREMENTS AND MAIN RESULTS: We identified a potentially damaging rare variant in the gene encoding for regulator of telomere elongation helicase 1 (RTEL1) that segregated with disease and was associated with very short telomeres in peripheral blood mononuclear cells in 1 of 25 families in our original whole-exome sequencing cohort. Evaluation of affected individuals in 163 additional kindreds revealed another eight families (4.7%) with heterozygous rare variants in RTEL1 that segregated with clinical FIP. Probands and unaffected carriers of these rare variants had short telomeres (<10% for age) in peripheral blood mononuclear cells and increased T-circle formation, suggesting impaired RTEL1 function. CONCLUSIONS: Rare loss-of-function variants in RTEL1 represent a newly defined genetic predisposition for FIP, supporting the importance of telomere-related pathways in pulmonary fibrosis.

Functional prostacyclin synthase promoter polymorphisms. Impact in pulmonary arterial hypertension.

RATIONALE: Pulmonary arterial hypertension (PAH) is a progressive disease characterized by elevated pulmonary artery pressure, vascular remodeling, and ultimately right ventricular heart failure. PAH can have a genetic component (heritable PAH), most often through mutations of bone morphogenetic protein receptor 2, and idiopathic and associated forms. Heritable PAH is not completely penetrant within families, with approximately 20% concurrence of inactivating bone morphogenetic protein receptor 2 mutations and delayed onset of PAH disease. Because one of the treatment options is using prostacyclin analogs, we hypothesized that prostacyclin synthase promoter sequence variants associated with increased mRNA expression may play a protective role in the bone morphogenetic protein receptor 2 unaffected carriers. OBJECTIVES: To characterize the range of prostacyclin synthase promoter variants and assess their transcriptional activities in PAH-relevant cell types. To determine the distribution of prostacyclin synthase promoter variants in PAH, unaffected carriers in heritable PAH families, and control populations. METHODS: Polymerase chain reaction approaches were used to genotype prostacyclin synthase promoter variants in more than 300 individuals. Prostacyclin synthase promoter haplotypes' transcriptional activities were determined with luciferase reporter assays. MEASUREMENTS AND MAIN RESULTS: We identified a comprehensive set of prostacyclin synthase promoter variants and tested their transcriptional activities in PAH-relevant cell types. We demonstrated differences of prostacyclin synthase promoter activities dependent on their haplotype. CONCLUSIONS: Prostacyclin synthase promoter sequence variants exhibit a range of transcriptional activities. We discovered a significant bias for more active prostacyclin synthase promoter variants in unaffected carriers as compared with affected patients with PAH.