Identification of skewed X chromosome inactivation using exome and transcriptome sequencing in patients with suspected rare genetic disease.

BACKGROUND: X-chromosome inactivation (XCI) is an epigenetic process that occurs during early development in mammalian females by randomly silencing one of two copies of the X chromosome in each cell. The preferential inactivation of either the maternal or paternal copy of the X chromosome in a majority of cells results in a skewed or non-random pattern of X inactivation and is observed in over 25% of adult females. Identifying skewed X inactivation is of clinical significance in patients with suspected rare genetic diseases due to the possibility of biased expression of disease-causing genes present on the active X chromosome. The current clinical test for the detection of skewed XCI relies on the methylation status of the methylation-sensitive restriction enzyme (Hpall) binding site present in proximity of short tandem polymorphic repeats on the androgen receptor (AR) gene. This approach using one locus results in uninformative or inconclusive data for 10-20% of tests. Further, recent studies have shown inconsistency between methylation of the AR locus and the state of inactivation of the X chromosome. Herein, we develop a method for estimating X inactivation status, using exome and transcriptome sequencing data derived from blood in 227 female samples. We built a reference model for evaluation of XCI in 135 females from the GTEx consortium. We tested and validated the model on 11 female individuals with different types of undiagnosed rare genetic disorders who were clinically tested for X-skew using the AR gene assay and compared results to our outlier-based analysis technique. RESULTS: In comparison to the AR clinical test for identification of X inactivation, our method was concordant with the AR method in 9 samples, discordant in 1, and provided a measure of X inactivation in 1 sample with uninformative clinical results. We applied this method on an additional 81 females presenting to the clinic with phenotypes consistent with different hereditary disorders without a known genetic diagnosis. CONCLUSIONS: This study presents the use of transcriptome and exome sequencing data to provide an accurate and complete estimation of X-inactivation and skew status in a cohort of female patients with different types of suspected rare genetic disease.

Semiautomated approach focused on new genomic information results in time and effort-efficient reannotation of negative exome data.

Most rare disease patients (75-50%) undergoing genomic sequencing remain unsolved, often due to lack of information about variants identified. Data review over time can leverage novel information regarding disease-causing variants and genes, increasing this diagnostic yield. However, time and resource constraints have limited reanalysis of genetic data in clinical laboratories setting. We developed RENEW, (REannotation of NEgative WES/WGS) an automated reannotation procedure that uses relevant new information in on-line genomic databases to enable rapid review of genomic findings. We tested RENEW in an unselected cohort of 1066 undiagnosed cases with a broad spectrum of phenotypes from the Mayo Clinic Center for Individualized Medicine using new information in ClinVar, HGMD and OMIM between the date of previous analysis/testing and April of 2022. 5741 variants prioritized by RENEW were rapidly reviewed by variant interpretation specialists. Mean analysis time was approximately 20 s per variant (32 h total time). Reviewed cases were classified as: 879 (93.0%) undiagnosed, 63 (6.6%) putatively diagnosed, and 4 (0.4%) definitively diagnosed. New strategies are needed to enable efficient review of genomic findings in unsolved cases. We report on a fast and practical approach to address this need and improve overall diagnostic success in patient testing through a recurrent reannotation process.

De novo PBX1 variant in a patient with glaucoma, kidney anomalies, and developmental delay: An expansion of the CAKUTHED phenotype.

An infant was referred for evaluation of congenital glaucoma and corneal clouding. In addition, he had a pelvic kidney, hypotonia, patent ductus arteriosus, abnormal pinnae, and developmental delay. Exome sequencing identified a previously unpublished de novo single nucleotide insertion in PBX1 c.400dupG (NM_002585.3), predicted to cause a frameshift resulting in a truncated protein with loss of function (p.Ala134Glyfs*65). Identification of this loss of function variant supports the diagnosis of congenital anomalies of the kidney and urinary tract syndrome with or without hearing loss, abnormal ears, or developmental delay (CAKUTHED). Here, we propose glaucoma as an extra-renal manifestation associated with PBX1-related disease due to the relationship of PBX1 with MEIS1, MEIS2, and FOXC1 transcription factors associated with eye development.

Monoallelic Mutations to DNAJB11 Cause Atypical Autosomal-Dominant Polycystic Kidney Disease.

Autosomal-dominant polycystic kidney disease (ADPKD) is characterized by the progressive development of kidney cysts, often resulting in end-stage renal disease (ESRD). This disorder is genetically heterogeneous with ∼7% of families genetically unresolved. We performed whole-exome sequencing (WES) in two multiplex ADPKD-like pedigrees, and we analyzed a further 591 genetically unresolved, phenotypically similar families by targeted next-generation sequencing of 65 candidate genes. WES identified a DNAJB11 missense variant (p.Pro54Arg) in two family members presenting with non-enlarged polycystic kidneys and a frameshifting change (c.166_167insTT) in a second family with small renal and liver cysts. DNAJB11 is a co-factor of BiP, a key chaperone in the endoplasmic reticulum controlling folding, trafficking, and degradation of secreted and membrane proteins. Five additional multigenerational families carrying DNAJB11 mutations were identified by the targeted analysis. The clinical phenotype was consistent in the 23 affected members, with non-enlarged cystic kidneys that often evolved to kidney atrophy; 7 subjects reached ESRD from 59 to 89 years. The lack of kidney enlargement, histologically evident interstitial fibrosis in non-cystic parenchyma, and recurring episodes of gout (one family) suggested partial phenotypic overlap with autosomal-dominant tubulointerstitial diseases (ADTKD). Characterization of DNAJB11-null cells and kidney samples from affected individuals revealed a pathogenesis associated with maturation and trafficking defects involving the ADPKD protein, PC1, and ADTKD proteins, such as UMOD. DNAJB11-associated disease is a phenotypic hybrid of ADPKD and ADTKD, characterized by normal-sized cystic kidneys and progressive interstitial fibrosis resulting in late-onset ESRD.

UBA2 variants underlie a recognizable syndrome with variable aplasia cutis congenita and ectrodactyly.

PURPOSE: The human chromosome 19q13.11 deletion syndrome is associated with a variable phenotype that includes aplasia cutis congenita (ACC) and ectrodactyly as specific features. UBA2 (ubiquitin-like modifier-activating enzyme 2) lies adjacent to the minimal deletion overlap region. We aimed to define the UBA2-related phenotypic spectrum in humans and zebrafish due to sequence variants and to establish the mechanism of disease. METHODS: Exome sequencing was used to detect UBA2 sequence variants in 16 subjects in 7 unrelated families. uba2 loss of function was modeled in zebrafish. Effects of human missense variants were assessed in zebrafish rescue experiments. RESULTS: Seven human UBA2 loss-of-function and missense sequence variants were detected. UBA2-phenotypes included ACC, ectrodactyly, neurodevelopmental abnormalities, ectodermal, skeletal, craniofacial, cardiac, renal, and genital anomalies. uba2 was expressed in zebrafish eye, brain, and pectoral fins; uba2-null fish showed deficient growth, microcephaly, microphthalmia, mandibular hypoplasia, and abnormal fins. uba2-mRNAs with human missense variants failed to rescue nullizygous zebrafish phenotypes. CONCLUSION: UBA2 variants cause a recognizable syndrome with a wide phenotypic spectrum. Our data suggest that loss of UBA2 function underlies the human UBA2 monogenic disorder and highlights the importance of SUMOylation in the development of affected tissues.

Impact of integrated translational research on clinical exome sequencing.

PURPOSE: Exome sequencing often identifies pathogenic genetic variants in patients with undiagnosed diseases. Nevertheless, frequent findings of variants of uncertain significance necessitate additional efforts to establish causality before reaching a conclusive diagnosis. To provide comprehensive genomic testing to patients with undiagnosed disease, we established an Individualized Medicine Clinic, which offered clinical exome testing and included a Translational Omics Program (TOP) that provided variant curation, research activities, or research exome sequencing. METHODS: From 2012 to 2018, 1101 unselected patients with undiagnosed diseases received exome testing. Outcomes were reviewed to assess impact of the TOP and patient characteristics on diagnostic rates through descriptive and multivariate analyses. RESULTS: The overall diagnostic yield was 24.9% (274 of 1101 patients), with 174 (15.8% of 1101) diagnosed on the basis of clinical exome sequencing alone. Four hundred twenty-three patients with nondiagnostic or without access to clinical exome sequencing were evaluated by the TOP, with 100 (9% of 1101) patients receiving a diagnosis, accounting for 36.5% of the diagnostic yield. The identification of a genetic diagnosis was influenced by the age at time of testing and the disease phenotype of the patient. CONCLUSION: Integration of translational research activities into clinical practice of a tertiary medical center can significantly increase the diagnostic yield of patients with undiagnosed disease.

Recurrent ganglioneuroma in PTPN11-associated Noonan syndrome: A case report and literature review.

Noonan syndrome (NS) is an autosomal dominant condition with variable expressivity most commonly due to a germline pathogenic variant in PTPN11, which encodes the protein tyrosine phosphatase SHP-2. Gain-of-function variants in PTPN11 are known to promote oncogenic behavior in affected tissues. We report the clinical description of a young adult male presenting with relapsing ganglioneuromas, dysmorphic features, cardiac abnormalities, and multiple lentigines, strongly suspicious for NS. Solid tumor testing identified the recurrent pathogenic c.922G>A (p.Asn308Asp) in PTPN11. Proband and parental blood sampling testing confirmed c.922G>A as a de novo germline alteration. Comprehensive literature review of solid tumors specifically associated to PTPN11, indicates that this is the first documentation of ganglioneuroma and its clinical recurrence after resection in conjunction with a genetically confirmed NS diagnosis. The findings in our patient further extend the list of neuroblastic and neural crest-derived neoplasms associated with this condition.

Synergistic Genetic Interactions between Pkhd1 and Pkd1 Result in an ARPKD-Like Phenotype in Murine Models.

BACKGROUND: Autosomal recessive polycystic kidney disease (ARPKD) and autosomal dominant polycystic kidney disease (ADPKD) are genetically distinct, with ADPKD usually caused by the genes PKD1 or PKD2 (encoding polycystin-1 and polycystin-2, respectively) and ARPKD caused by PKHD1 (encoding fibrocystin/polyductin [FPC]). Primary cilia have been considered central to PKD pathogenesis due to protein localization and common cystic phenotypes in syndromic ciliopathies, but their relevance is questioned in the simple PKDs. ARPKD's mild phenotype in murine models versus in humans has hampered investigating its pathogenesis. METHODS: To study the interaction between Pkhd1 and Pkd1, including dosage effects on the phenotype, we generated digenic mouse and rat models and characterized and compared digenic, monogenic, and wild-type phenotypes. RESULTS: The genetic interaction was synergistic in both species, with digenic animals exhibiting phenotypes of rapidly progressive PKD and early lethality resembling classic ARPKD. Genetic interaction between Pkhd1 and Pkd1 depended on dosage in the digenic murine models, with no significant enhancement of the monogenic phenotype until a threshold of reduced expression at the second locus was breached. Pkhd1 loss did not alter expression, maturation, or localization of the ADPKD polycystin proteins, with no interaction detected between the ARPKD FPC protein and polycystins. RNA-seq analysis in the digenic and monogenic mouse models highlighted the ciliary compartment as a common dysregulated target, with enhanced ciliary expression and length changes in the digenic models. CONCLUSIONS: These data indicate that FPC and the polycystins work independently, with separate disease-causing thresholds; however, a combined protein threshold triggers the synergistic, cystogenic response because of enhanced dysregulation of primary cilia. These insights into pathogenesis highlight possible common therapeutic targets.

Prevalence, Penetrance, and Phenotypic Manifestation of Cardiomyopathy-Associated Genetic Variants in the General Population: Insights from a Mayo Clinic Biobank Study.

OBJECTIVE: To determine the prevalence, penetrance, and disease expression of cardiomyopathy-related genetic variants in an unselected, richly phenotyped Mayo Clinic population in the setting of preemptive sequencing, with return of incidental findings following the American College of Medical Genetics and Genomics recommendations. PATIENTS AND METHODS: We analyzed a quaternary medical center-based biobank cohort (n=983) for reportable variants in 15 cardiomyopathy genes. Prioritization of genetic variants was performed using an internally developed pipeline to identify potentially reportable variants. Prioritized variants were then manually curated. The correlation of likely pathogenic/pathogenic (LP/P) variants with clinical phenotypes and outcomes was established. Artificial intelligence-enabled electrocardiographic predictions of reduced left ventricular ejection fraction and hypertrophic cardiomyopathy were applied to genotype-positive (G+) participants. RESULTS: Of the 983 patients, 11 (1%) were G+, with 11 LP/P variants found in the MYBPC3, DSG2, MYH7, DSP, and PKP2 genes. All G+ participants underwent electrocardiography, and 10 (90%) underwent echocardiography. Most patients (10 [90%]) did not have a prior diagnosis of cardiomyopathy. Definitive disease penetrance (heart failure or cardiomyopathy) was present in 4 (36%), while 3 (27%) had possible penetrance (structural heart disease identified by echocardiography). Arrhythmias and/or cardiac conduction disease was present in 4 of 11 G+ individuals (36%). Artificial intelligence-electrocardiography was positive for hypertrophic cardiomyopathy or reduced left ventricular ejection fraction in 5 of the G+ participants (45%), of whom 4 (80%) had definitive or possible disease penetrance. CONCLUSION: Cardiomyopathy-associated LP/P variants are present in a small subset of a quaternary medical center population, and disease penetrance in G+ individuals is high in the form of cardiac structural abnormalities and heart failure.

Bi-allelic variants in CELSR3 are implicated in central nervous system and urinary tract anomalies.

CELSR3 codes for a planar cell polarity protein. We describe twelve affected individuals from eleven independent families with bi-allelic variants in CELSR3. Affected individuals presented with an overlapping phenotypic spectrum comprising central nervous system (CNS) anomalies (7/12), combined CNS anomalies and congenital anomalies of the kidneys and urinary tract (CAKUT) (3/12) and CAKUT only (2/12). Computational simulation of the 3D protein structure suggests the position of the identified variants to be implicated in penetrance and phenotype expression. CELSR3 immunolocalization in human embryonic urinary tract and transient suppression and rescue experiments of Celsr3 in fluorescent zebrafish reporter lines further support an embryonic role of CELSR3 in CNS and urinary tract formation.

Identification of AFG3L2 dominant optic atrophy following reanalysis of clinical exome sequencing.

Purpose: To highlight the importance of the utility of clinical exome sequencing, and show how it led to the diagnosis of nonsyndromic autosomal dominant optic atrophy arising from an autosomal dominant variant in AFG3L2. Observations: A healthy father and daughter of East African heritage experienced the onset of vision loss in the first decade of life due to optic atrophy. No additional neurologic or neuroimaging abnormalities were detected. Clinical exome sequencing was initially performed and provided a negative result. Reanalysis of the sequencing data revealed an autosomal dominant pathogenic variant in AFG3L2, c.1064C>T (p.Thr355Met), a gene that was recently identified to be associated with non-syndromic optic atrophy. This variant has previously been reported in a patient with optic atrophy, motor disturbances, and an abnormal brain MRI. Conclusions: As the causes of dominant optic atrophy continue to expand, accurate genetic diagnosis is aided by an iterative reanalysis process for individuals and families when initial exome and genome testing does not provide an answer.

Neurodevelopmental and Epilepsy Phenotypes in Individuals With Missense Variants in the Voltage-Sensing and Pore Domains of KCNH5.

BACKGROUND AND OBJECTIVES: KCNH5 encodes the voltage-gated potassium channel EAG2/Kv10.2. We aimed to delineate the neurodevelopmental and epilepsy phenotypic spectrum associated with de novo KCNH5 variants. METHODS: We screened 893 individuals with developmental and epileptic encephalopathies for KCNH5 variants using targeted or exome sequencing. Additional individuals with KCNH5 variants were identified through an international collaboration. Clinical history, EEG, and imaging data were analyzed; seizure types and epilepsy syndromes were classified. We included 3 previously published individuals including additional phenotypic details. RESULTS: We report a cohort of 17 patients, including 9 with a recurrent de novo missense variant p.Arg327His, 4 with a recurrent missense variant p.Arg333His, and 4 additional novel missense variants. All variants were located in or near the functionally critical voltage-sensing or pore domains, absent in the general population, and classified as pathogenic or likely pathogenic using the American College of Medical Genetics and Genomics criteria. All individuals presented with epilepsy with a median seizure onset at 6 months. They had a wide range of seizure types, including focal and generalized seizures. Cognitive outcomes ranged from normal intellect to profound impairment. Individuals with the recurrent p.Arg333His variant had a self-limited drug-responsive focal or generalized epilepsy and normal intellect, whereas the recurrent p.Arg327His variant was associated with infantile-onset DEE. Two individuals with variants in the pore domain were more severely affected, with a neonatal-onset movement disorder, early-infantile DEE, profound disability, and childhood death. DISCUSSION: We describe a cohort of 17 individuals with pathogenic or likely pathogenic missense variants in the voltage-sensing and pore domains of Kv10.2, including 14 previously unreported individuals. We present evidence for a putative emerging genotype-phenotype correlation with a spectrum of epilepsy and cognitive outcomes. Overall, we expand the role of EAG proteins in human disease and establish KCNH5 as implicated in a spectrum of neurodevelopmental disorders and epilepsy.

Childhood-onset dystonia-causing KMT2B variants result in a distinctive genomic hypermethylation profile.

BACKGROUND: Dystonia is a clinically and genetically heterogeneous movement disorder characterized by sustained or intermittent muscle contractions causing abnormal, often repetitive, movements and/or postures. Heterozygous variants in lysine methyltransferase 2B (KMT2B), encoding a histone H3 methyltransferase, have been associated with a childhood-onset, progressive and complex form of dystonia (dystonia 28, DYT28). Since 2016, more than one hundred rare KMT2B variants have been reported, including frameshift, nonsense, splice site, missense and other in-frame changes, many having an uncertain clinical impact. RESULTS: We characterize the genome-wide peripheral blood DNA methylation profiles of a cohort of 18 patients with pathogenic and unclassified KMT2B variants. We resolve the "episignature" associated with KMT2B haploinsufficiency, proving that this approach is robust in diagnosing clinically unsolved cases, properly classifying them with respect to other partially overlapping dystonic phenotypes, other rare neurodevelopmental disorders and healthy controls. Notably, defective KMT2B function in DYT28 causes a non-random DNA hypermethylation across the genome, selectively involving promoters and other regulatory regions positively controlling gene expression. CONCLUSIONS: We demonstrate a distinctive DNA hypermethylation pattern associated with DYT28, provide an epigenetic signature for this disorder enabling accurate diagnosis and reclassification of ambiguous genetic findings and suggest potential therapeutic approaches.

Utility of survival motor neuron ELISA for spinal muscular atrophy clinical and preclinical analyses.

OBJECTIVES: Genetic defects leading to the reduction of the survival motor neuron protein (SMN) are a causal factor for Spinal Muscular Atrophy (SMA). While there are a number of therapies under evaluation as potential treatments for SMA, there is a critical lack of a biomarker method for assessing efficacy of therapeutic interventions, particularly those targeting upregulation of SMN protein levels. Towards this end we have engaged in developing an immunoassay capable of accurately measuring SMN protein levels in blood, specifically in peripheral blood mononuclear cells (PBMCs), as a tool for validating SMN protein as a biomarker in SMA. METHODS: A sandwich enzyme-linked immunosorbent assay (ELISA) was developed and validated for measuring SMN protein in human PBMCs and other cell lysates. Protocols for detection and extraction of SMN from transgenic SMA mouse tissues were also developed. RESULTS: The assay sensitivity for human SMN is 50 pg/mL. Initial analysis reveals that PBMCs yield enough SMN to analyze from blood volumes of less than 1 mL, and SMA Type I patients' PBMCs show ∼90% reduction of SMN protein compared to normal adults. The ELISA can reliably quantify SMN protein in human and mouse PBMCs and muscle, as well as brain, and spinal cord from a mouse model of severe SMA. CONCLUSIONS: This SMN ELISA assay enables the reliable, quantitative and rapid measurement of SMN in healthy human and SMA patient PBMCs, muscle and fibroblasts. SMN was also detected in several tissues in a mouse model of SMA, as well as in wildtype mouse tissues. This SMN ELISA has general translational applicability to both preclinical and clinical research efforts.