Search

Abstract Background Though next-generation sequencing (NGS) tests like exome sequencing (ES), genome sequencing (GS), and panels derived from exome and genome data (EGBP) are effective for rare diseases, the ideal diagnostic approach is debated. Limited research has explored reanalyzing raw ES and GS data post-negative EGBP results for diagnostics. Results: We analyzed complete ES/GS raw sequencing data from Mayo Clinic's Program for Rare and Undiagnosed Diseases (PRaUD) patients to assess whether supplementary findings could augment diagnostic yield. ES data from 80 patients (59 adults) and GS data from 20 patients (10 adults), averaging 43 years in age, were analyzed. Most patients had renal (n=44) and auto-inflammatory (n=29) phenotypes. Ninety-six cases had negative findings and in four cases additional genetic variants were found, including a variant related to a recently described disease (RRAGD-related hypomagnesemia), a variant missed due to discordant inheritance pattern (COL4A3), a variant with high allelic frequency (NPHS2) in the general population, and a variant associated with an initially untargeted phenotype (HNF1A). Conclusion: ES and GS show diagnostic yields comparable to EGBP for single-system diseases. However, EGBP's limitations in detecting new disease-associated genes underscore the necessity for periodic updates.

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

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

Abstract Background In the United States, rare disease (RD) is defined as a condition that affects fewer than 200,000 individuals. Collectively, RD affects an estimated 30 million Americans. A significant portion of RD has an underlying genetic cause; however, this may go undiagnosed. To better serve these patients, the Mayo Clinic Program for Rare and Undiagnosed Diseases (PRaUD) was created under the auspices of the Center for Individualized Medicine (CIM) aiming to integrate genomics into subspecialty practice including targeted genetic testing, research, and education. Methods Patients were identified by subspecialty healthcare providers from 11 clinical divisions/departments. Targeted multi-gene panels or custom exome/genome-based panels were utilized. To support the goals of PRaUD, a new clinical service model, the Genetic Testing and Counseling (GTAC) unit, was established to improve access and increase efficiency for genetic test facilitation. The GTAC unit includes genetic counselors, genetic counseling assistants, genetic nurses, and a medical geneticist. Patients receive abbreviated point-of-care genetic counseling and testing through a partnership with subspecialty providers. Results Implementation of PRaUD began in 2018 and GTAC unit launched in 2020 to support program expansion. Currently, 29 RD clinical indications are included in 11 specialty divisions/departments with over 142 referring providers. To date, 1152 patients have been evaluated with an overall solved or likely solved rate of 17.5% and as high as 66.7% depending on the phenotype. Noteworthy, 42.7% of the solved or likely solved patients underwent changes in medical management and outcome based on genetic test results. Conclusion Implementation of PRaUD and GTAC have enabled subspecialty practices advance expertise in RD where genetic counselors have not historically been embedded in practice. Democratizing access to genetic testing and counseling can broaden the reach of patients with RD and increase the diagnostic yield of such indications leading to better medical management as well as expanding research opportunities.

ABSTRACTSevere acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has had an unprecedented impact on human health and highlights the need for genomic epidemiology studies to increase our understanding of virus evolution and spread and to inform policy decisions. We sequenced viral genomes from over 22,000 patient samples tested at Mayo Clinic Laboratories between 2020 and 2022 and used Bayesian phylodynamics to describe county and regional spread in Minnesota. The earliest calculated introduction into Minnesota was to Hennepin County from a domestic source around 22 January 2020; 6 weeks before the first confirmed case in the state. This led to the virus spreading to Northern Minnesota and, eventually, the rest of the state. International introductions were most abundant in Hennepin (home to the Minneapolis/St. Paul International Airport) totaling 45 (out of 107) over the 2-year period. Southern Minnesota counties were most common for domestic introductions, with 19 (out of 64), potentially driven by bordering states such as Iowa and Wisconsin as well as Illinois, which is nearby. Hennepin also was, by far, the most dominant source of in-state transmissions to other Minnesota locations (n = 772) over the 2-year period. We also analyzed the diversity of the location source of SARS-CoV-2 viruses in each county and noted the timing of state-wide policies as well as trends in clinical cases. Neither the number of clinical cases nor the major policy decisions, such as the end of the lockdown period in 2020 or the end of all restrictions in 2021, appeared to have an impact on virus diversity across each individual county.IMPORTANCEWe analyzed over 22,000 severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genomes of patient samples tested at Mayo Clinic Laboratories during a 2-year period in the COVID-19 pandemic, which included Alpha, Delta, and Omicron variants of concern to examine the roles and relationships of Minnesota virus transmission. We found that Hennepin County, the most populous county, drove the transmission of SARS-CoV-2 viruses in the state after including the formation of earlier clades including 20A, 20C, and 20G, as well as variants of concern Alpha and Delta. We also found that Hennepin County was the source for most of the county-to-county introductions after an initial predicted introduction with the virus in early 2020 from an international source, while other counties acted as transmission “sinks.” In addition, major policies, such as the end of the lockdown period in 2020 or the end of all restrictions in 2021, did not appear to have an impact on virus diversity across individual counties.

Abstract Background The PI3K/AKT pathway, extensively studied in cancer, is vital for regulating cell metabolism, differentiation, and proliferation. Pathogenic variants in the PIK3R1 gene, which encodes three regulatory units of class IA PI3Ks, have been found in affected tissue of individuals with vascular lesions. These variants predominantly occur in the iSH2 domain, disrupting inhibitory contacts with the catalytic unit and leading to PI3K activation. Germline variants in this gene are also linked to an immunological condition called Activated PI3K delta syndrome type 2 (APDS2). Methods This is a case report and literature review. Clinical data were retrieved from medical records. Results A male patient presented with extensive vascular malformation covering over 90% of his body, along with complete 2–3 toe syndactyly, suggesting a vascular malformation syndrome called PROS. Low levels of IgA and IgG were detected. The patient achieved his developmental milestones and had above‐average weight, height, and head circumference. Exome sequencing of skin and blood DNA revealed a de novo variant in PIK3R1 (c.1746‐2A>G, p.?) in 9% of the patient's blood cells and 25% of cultured fibroblasts. Initially, classified as a variant of uncertain significance, this variant was later confirmed to be the cause. Conclusions This is the first intronic SNV in a canonical splice site within iSH2 described, highlighting the importance of iSH2 in the regulation of the PI3K/AKT pathway and its involvement in the development of vascular overgrowth and antibody deficiency.

Summary: Autoimmune polyendocrine syndrome type 1 (APS-1) is an autosomal recessive disease characterized by severe and childhood onset organ-specific autoimmunity caused by mutations in the autoimmune regulator (AIRE) gene. More recently, dominant-negative mutations within the PHD1, PHD2, and SAND domains have been associated with an incompletely penetrant milder phenotype with later onset familial clustering, often masquerading as organ-specific autoimmunity. Patients with immunodeficiencies or autoimmunity where genetic analyses revealed heterozygous AIRE mutations were included in the study and the dominant-negative effects of the AIRE mutations were functionally assessed in vitro. We here report additional families with phenotypes ranging from immunodeficiency, enteropathy, and vitiligo to asymptomatic carrier status. APS-1-specific autoantibodies can hint to the presence of these pathogenic AIRE variants although their absence does not rule out their presence. Our findings suggest functional studies of heterozygous AIRE variants and close follow-up of identified individuals and their families.

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.

Abstract Background VEXAS syndrome (vacuoles, E1 enzyme, X-linked, autoinflammatory, somatic syndrome) is a recently described syndrome caused by a somatic missense variant at the methionine-41 (p.(Met41)) position in the ubiquitin-like modifier activating enzyme 1 (UBA1) in Xp11.3. Germline pathogenic variants in UBA1 are associated with a distinct phenotype: a syndrome with severe neurologic features associated with loss of anterior horn cells and infantile death denominated X-Linked Spinal Muscular Atrophy 2 (SMAX2) (OMIM 301,830). Case presentation We report a male individual with the phenotype of VEXAS syndrome that was initially identified through exome sequencing (ES) as having a hemizygous germline variant in UBA1 due to high variant allele frequency (VAF). Research Sanger sequencing was able to confirm the absence of the p.(Met41Val) variant in a skin biopsy and in gastric mucosa tissue sample confirming the variant happened as a postzygotic event. Conclusions The present case exemplifies the diagnostic challenge that was imposed by the high VAF detected by ES that failed to correctly demonstrate that the variant was in a mosaic state. Sequencing of different tissues should be considered when there is conflict between the UBA1 variant status and the clinical findings.

Abstract Background Genomics enables individualized diagnosis and treatment, but large challenges remain to functionally interpret rare variants. To date, only one causative variant has been described for KCNK9 imprinting syndrome (KIS). The genotypic and phenotypic spectrum of KIS has yet to be described and the precise mechanism of disease fully understood. Methods This study discovers mechanisms underlying KCNK9 imprinting syndrome (KIS) by describing 15 novel KCNK9 alterations from 47 KIS-affected individuals. We use clinical genetics and computer-assisted facial phenotyping to describe the phenotypic spectrum of KIS. We then interrogate the functional effects of the variants in the encoded TASK3 channel using sequence-based analysis, 3D molecular mechanic and dynamic protein modeling, and in vitro electrophysiological and functional methodologies. Results We describe the broader genetic and phenotypic variability for KIS in a cohort of individuals identifying an additional mutational hotspot at p.Arg131 and demonstrating the common features of this neurodevelopmental disorder to include motor and speech delay, intellectual disability, early feeding difficulties, muscular hypotonia, behavioral abnormalities, and dysmorphic features. The computational protein modeling and in vitro electrophysiological studies discover variability of the impact of KCNK9 variants on TASK3 channel function identifying variants causing gain and others causing loss of conductance. The most consistent functional impact of KCNK9 genetic variants, however, was altered channel regulation. Conclusions This study extends our understanding of KIS mechanisms demonstrating its complex etiology including gain and loss of channel function and consistent loss of channel regulation. These data are rapidly applicable to diagnostic strategies, as KIS is not identifiable from clinical features alone and thus should be molecularly diagnosed. Furthermore, our data suggests unique therapeutic strategies may be needed to address the specific functional consequences of KCNK9 variation on channel function and regulation.

Background: Rare genetic disease studies have benefited from the era of high throughput sequencing. DNA sequencing results in genetic diagnosis of 18-40% of previously unsolved cases, while the incorporation of RNA-Seq analysis has more recently been shown to generate significant numbers of previously unattainable diagnoses. While DNA methylation remains less explored, multiple inborn diseases resulting from disorders of genomic imprinting are well characterized and a growing body of literature suggests the causative or correlative role of aberrant methylation in diverse rare inherited conditions. Complex pictures of methylation patterning are also emerging, including the association of regional, multiple specific-site or even single-site methylation, with disease. The systematic application of genomic-wide methylation-based sequencing for undiagnosed cases of rare diseases is a logical progression from current testing paradigms. Similar to the rationale previously exploited in RNA-based rare disease studies, we can assume that disease-associated or causative methylation aberrations in an individual will demonstrate significant differences from other individuals with unrelated phenotypes. Thus, aberrantly methylated sites will be outliers from a heterogeneous cohort of individuals. Methods: Based on this rationale, we present BOREALIS: Bisulfite-seq Outlie R M Ethylation At Sing Le S Ite Re Solution. BOREALIS uses a beta binomial model to identify outlier methylation at single CpG site resolution from bisulfite sequencing data. Results: Utilizing power analyses, we demonstrate that BOREALIS can identify outlier CpG methylation within a cohort of samples. Furthermore, we show that BOREALIS is tolerant to the inclusion of multiple identical outliers with sufficient cohort size and sequencing depth. Conclusions: The method demonstrates improved performance versus standard statistical testing and is suited for single or multi-site downstream analysis.

Abstract The changing landscape of genomics research and clinical practice has created a need for computational pipelines capable of efficiently orchestrating complex analysis stages while handling large volumes of data across heterogeneous computational environments. Workflow Management Systems (WfMSs) are the software components employed to fill this gap. This work provides an approach and systematic evaluation of key features of popular bioinformatics WfMSs in use today: Nextflow, CWL, and WDL and some of their executors, along with Swift/T, a workflow manager commonly used in high-scale physics applications. We employed two use cases: a variant-calling genomic pipeline and a scalability-testing framework, where both were run locally, on an HPC cluster, and in the cloud. This allowed for evaluation of those four WfMSs in terms of language expressiveness, modularity, scalability, robustness, reproducibility, interoperability, ease of development, along with adoption and usage in research labs and healthcare settings. This article is trying to answer, which WfMS should be chosen for a given bioinformatics application regardless of analysis type?. The choice of a given WfMS is a function of both its intrinsic language and engine features. Within bioinformatics, where analysts are a mix of dry and wet lab scientists, the choice is also governed by collaborations and adoption within large consortia and technical support provided by the WfMS team/community. As the community and its needs continue to evolve along with computational infrastructure, WfMSs will also evolve, especially those with permissive licenses that allow commercial use. In much the same way as the dataflow paradigm and containerization are now well understood to be very useful in bioinformatics applications, we will continue to see innovations of tools and utilities for other purposes, like big data technologies, interoperability, and provenance.

Rationale & Objective: The etiology of kidney disease remains unknown in many individuals with chronic kidney disease (CKD). We created the Mayo Clinic Nephrology Genomics Clinic to improve our ability to integrate genomic and clinical data to identify the etiology of unexplained CKD. Study Design: Retrospective study. Setting & Participants: An essential component of our program is the Nephrology Genomics Board which consists of nephrologists, geneticists, pathologists, translational omics scientists, and trainees who interpret the patient’s clinical and genetic data. Since September 2016, the Board has reviewed 163 cases (15 cystic, 100 glomerular, 6 congenital anomalies of kidney and urinary tract (CAKUT), 20 stones, 15 tubulointerstitial, and 13 other). Analytical Approach: Testing was performed with targeted panels, single gene analysis, or analysis of kidney-related genes from exome sequencing. Variant classification was obtained based on the 2015 American College of Medical Genetics and Genomics and the Association for Molecular Pathology guidelines. Results: A definitive genetic diagnosis was achieved for 50 families (30.7%). The highest diagnostic yield was obtained in individuals with tubulointerstitial diseases (53.3%), followed by congenital anomalies of the kidney and urological tract (33.3%), glomerular (31%), cysts (26.7%), stones (25%), and others (15.4%). A further 20 (12.3%) patients had variants of interest, and variant segregation, and research activities (exome, genome, or transcriptome sequencing) are ongoing for 44 (40%) unresolved families. Limitations: Possible overestimation of diagnostic rate due to inclusion of individuals with variants with evidence of pathogenicity but classified as of uncertain significance by the clinical laboratory. Conclusions: Integration of genomic and research testing and multidisciplinary evaluation in a nephrology cohort with CKD of unknown etiology or suspected monogenic disease provided a diagnosis in a third of families. These diagnoses had prognostic implications, and often changes in management were implemented.

Abstract Background Modern Next Generation- and Third Generation- Sequencing methods such as Illumina and PacBio Circular Consensus Sequencing platforms provide accurate sequencing data. Parallel developments in Deep Learning have enabled the application of Deep Neural Networks to variant calling, surpassing the accuracy of classical approaches in many settings. DeepVariant, arguably the most popular among such methods, transforms the problem of variant calling into one of image recognition where a Deep Neural Network analyzes sequencing data that is formatted as images, achieving high accuracy. In this paper, we explore an alternative approach to designing Deep Neural Networks for variant calling, where we use meticulously designed Deep Neural Network architectures and customized variant inference functions that account for the underlying nature of sequencing data instead of converting the problem to one of image recognition. Results Results from 27 whole-genome variant calling experiments spanning Illumina, PacBio and hybrid Illumina-PacBio settings suggest that our method allows vastly smaller Deep Neural Networks to outperform the Inception-v3 architecture used in DeepVariant for indel and substitution-type variant calls. For example, our method reduces the number of indel call errors by up to 18%, 55% and 65% for Illumina, PacBio and hybrid Illumina-PacBio variant calling respectively, compared to a similarly trained DeepVariant pipeline. In these cases, our models are between 7 and 14 times smaller. Conclusions We believe that the improved accuracy and problem-specific customization of our models will enable more accurate pipelines and further method development in the field. HELLO is available at https://github.com/anands-repo/hello

Abstract Background Achalasia‐addisonianism‐alacrima syndrome, frequently referred to as Allgrove syndrome or Triple A syndrome, is a multisystem disorder resulting from homozygous or compound heterozygous pathogenic variants in the gene encoding aladin (AAAS). Aladin is a member of the WD‐repeat family of proteins and is a component of the nuclear pore complex. It is thought to be involved in nuclear import and export of molecules. Here, we describe an individual with a paternally inherited truncating variant and a maternally inherited, novel missense variant in AAAS presenting with alacrima, achalasia, anejaculation, optic atrophy, muscle weakness, dysarthria, and autonomic dysfunction. Methods Whole‐exome sequencing was performed in the proband, sister, and parents. Variants were confirmed by Sanger sequencing. The localization of aladin to the nuclear pore was assessed in cells expressing the patient variant. Results Functional testing of the maternally inherited variant, p.(Arg270Pro), demonstrated decreased localization of aladin to the nuclear pore in cells expressing the variant, indicating a deleterious effect. Follow‐up testing in the proband's affected sister revealed that she also inherited the biallelic AAAS variants. Conclusions Review of the patient's clinical, pathological, and genetic findings resulted in a diagnosis of Triple A syndrome.

The chromokinesin KIF22 generates forces that contribute to mitotic chromosome congression and alignment. Mutations in the α2 helix of the motor domain of KIF22 have been identified in patients with abnormal skeletal development, and we report the identification of a patient with a novel mutation in the KIF22 tail. We demonstrate that pathogenic mutations do not result in a loss of KIF22’s functions in early mitosis. Instead, mutations disrupt chromosome segregation in anaphase, resulting in reduced proliferation, abnormal daughter cell nuclear morphology, and, in a subset of cells, cytokinesis failure. This phenotype could be explained by a failure of KIF22 to inactivate in anaphase. Consistent with this model, constitutive activation of the motor via a known site of phosphoregulation in the tail phenocopied the effects of pathogenic mutations. These results suggest that the motor domain α2 helix may be an important site for regulation of KIF22 activity at the metaphase to anaphase transition. In support of this conclusion, mimicking phosphorylation of α2 helix residue T158 also prevents inactivation of KIF22 in anaphase. These findings demonstrate the importance of both the head and tail of the motor in regulating the activity of KIF22 and offer insight into the cellular consequences of preventing KIF22 inactivation and disrupting force balance in anaphase.

Abstract T cell prolymphocytic leukemia (T-PLL) is a rare disease with aggressive clinical course. Cytogenetic analysis, whole-exome and whole-genome sequencing have identified primary structural alterations in T-PLL, including inversion, translocation and copy number variation. Recurrent somatic mutations were also identified in genes encoding chromatin regulators and those in the JAK-STAT signaling pathway. Epigenetic alterations are the hallmark of many cancers. However, genome-wide epigenomic profiles have not been reported in T-PLL, limiting the mechanistic study of its carcinogenesis. We hypothesize epigenetic mechanisms also play a key role in T-PLL pathogenesis. To systematically test this hypothesis, we generated genome-wide maps of regulatory regions using H3K4me3 and H3K27ac ChIP-seq, as well as RNA-seq data in both T-PLL patients and healthy individuals. We found that genes down-regulated in T-PLL are mainly associated with defense response, immune system or adaptive immune response, while up-regulated genes are enriched in developmental process, as well as WNT signaling pathway with crucial roles in cell fate decision. In particular, our analysis revealed a global alteration of regulatory landscape in T-PLL, with differential peaks highly enriched for binding motifs of immune related transcription factors, supporting the epigenetic regulation of oncogenes and genes involved in DNA damage response and T-cell activation. Together, our work reveals a causal role of epigenetic dysregulation in T-PLL.

Abstract Background Damaging variants in TRIO have been associated with moderate to severe neurodevelopmental disorders in humans. While recent work has delineated the positional effect of missense variation on the resulting phenotype, the clinical spectrum associated with loss-of-function variation has yet to be fully defined. Case presentation We report on two probands with novel loss-of-function variants in TRIO. Patient 1 presents with a severe neurodevelopmental disorder and macrocephaly. The TRIO variant is inherited from his affected mother. Patient 2 presents with moderate developmental delays, microcephaly, and cutis aplasia with a frameshift variant of unknown inheritance. Conclusions We describe two patients with neurodevelopmental disorder, macro/microcephaly, and cutis aplasia in one patient. Both patients have loss-of-function variants, helping to further characterize how these types of variants affect the phenotypic spectrum associated with TRIO. We also present the third reported case of autosomal dominant inheritance of a damaging variant in TRIO.

Abstract Background Genetic changes in the LIM homeobox transcription factor 1 beta (LMX1B) have been associated with focal segmental glomerulosclerosis (FSGS) without the extra-renal or ultrastructural manifestations of Nail-patella syndrome (NPS) known as Nail-patella-like renal disease (NPLRD). Fabry disease (FD) is an X-linked lysosomal disease caused by the deficiency of alpha-galactosidase A. The classic form of the disease is characterized by acroparesthesia, angiokeratomas, cornea verticillata, hypertrophic cardiomyopathy, strokes, and chronic kidney disease. Podocyte myelin bodies on ultrastructural examination of kidney tissue are very characteristic of FD; however some medications and other conditions may mimic this finding. Case presentation Here, we report on a female patient with chronic kidney disease (CKD), positive family history for kidney disease and kidney biopsy showing a FSGS lesion and presence of focal myelin figures within podocytes concerning for FD. However, genetic testing for FD was negative. After comprehensive clinical, biochemical, and genetic evaluation, including whole exome and RNA sequencing, she was ultimately diagnosed with NPLRD. Conclusions This case illustrates the difficulties of diagnosing atypical forms of rare Mendelian kidney diseases and the role of a multidisciplinary team in an individualized medicine clinic setting in combination with state-of-the-art sequencing technologies to reach a definitive diagnosis.

Abstract Mitochondria play a variety of roles in the cell, far beyond their widely recognized role in ATP generation. One such role is the regulation and sequestration of calcium, which is done with the help of the mitochondrial calcium uniporter (MCU) and its regulators, MICU1 and MICU2. Genetic variations in MICU1 and MICU2 have been reported to cause myopathy, developmental disability and neurological symptoms typical of mitochondrial disorders. The symptoms of MICU1/2 deficiency have generally been attributed to calcium regulation in the metabolic and biochemical roles of mitochondria. Here, we report a female child with heterozygous MICU1 variants and multiple congenital brain malformations on MRI. Specifically, she shows anterior perisylvian polymicrogyria, dysmorphic basal ganglia, and cerebellar dysplasia in addition to white matter abnormalities. These novel findings suggest that MICU1 is necessary for proper neurodevelopment through a variety of potential mechanisms, including calcium‐mediated regulation of the neuronal cytoskeleton, Miro1‐MCU complex‐mediated mitochondrial movement, or enhancing ATP production. This case provides new insight into the molecular pathogenesis of MCU dysfunction and may represent a novel diagnostic feature of calcium‐based mitochondrial disease.

UDP-glucuronic acid is a component of the extracellular matrix. Here, the authors report biallelic variants in the gene encoding UDP-Glucose 6-Dehydrogenase (UGDH) in individuals affected by developmental epileptic encephalopathies that impair UGDH stability, oligomerization, or enzymatic activity in vitro.

The hereditary ataxias are a heterogenous group of disorders with an increasing number of causative genes being described. Due to the clinical and genetic heterogeneity seen in these conditions, the majority of such individuals endure a diagnostic odyssey or remain undiagnosed. Defining the molecular etiology can bring insights into the responsible molecular pathways and eventually the identification of therapeutic targets. Here, we describe the identification of biallelic variants in the GEMIN5 gene among seven unrelated families with nine affected individuals presenting with spastic ataxia and cerebellar atrophy. GEMIN5, an RNA-binding protein, has been shown to regulate transcription and translation machinery. GEMIN5 is a component of small nuclear ribonucleoprotein (snRNP) complexes and helps in the assembly of the spliceosome complexes. We found that biallelic GEMIN5 variants cause structural abnormalities in the encoded protein and reduce expression of snRNP complex proteins in patient cells compared with unaffected controls. Finally, knocking out endogenous Gemin5 in mice caused early embryonic lethality, suggesting that Gemin5 expression is crucial for normal development. Our work further expands on the phenotypic spectrum associated with GEMIN5-related disease and implicates the role of GEMIN5 among patients with spastic ataxia, cerebellar atrophy, and motor predominant developmental delay.

Detecting gene fusions involving driver oncogenes is pivotal in clinical diagnosis and treatment of cancer patients. Recent developments in next-generation sequencing (NGS) technologies have enabled improved assays for bioinformatics-based gene fusions detection. In clinical applications, where a small number of fusions are clinically actionable, targeted polymerase chain reaction (PCR)-based NGS chemistries, such as the QIAseq RNAscan assay, aim to improve accuracy compared to standard RNA sequencing. Existing informatics methods for gene fusion detection in NGS-based RNA sequencing assays traditionally use a transcriptome-based spliced alignment approach or a de-novo assembly approach. Transcriptome-based spliced alignment methods face challenges with short read mapping yielding low quality alignments. De-novo assembly-based methods yield longer contigs from short reads that can be more sensitive for genomic rearrangements, but face performance and scalability challenges. Consequently, there exists a need for a method to efficiently and accurately detect fusions in targeted PCR-based NGS chemistries. We describe SeekFusion, a highly accurate and computationally efficient pipeline enabling identification of gene fusions from PCR-based NGS chemistries. Utilizing biological samples processed with the QIAseq RNAscan assay and in-silico simulated data we demonstrate that SeekFusion gene fusion detection accuracy outperforms popular existing methods such as STAR-Fusion, TOPHAT-Fusion and JAFFA-hybrid. We also present results from 4,484 patient samples tested for neurological tumors and sarcoma, encompassing details on some novel fusions identified.

Abstract Background Biallelic pathogenic variants in HTRA1 cause CARASIL. More recently, monoallelic variants have been associated with the autosomal dominant disorder CADASIL2 but not all carriers develop disease manifestations. We describe the clinicoradiologic and mutation spectrum of four new CADASIL2 individuals. Methods Medical records at Mayo Clinic between 2013 and 2020 were retrospectively reviewed to identify patients with cerebral small vessel disease related to monoallelic HTRA1 variants. Results Four patients met the study inclusion criteria for cerebral small vessel disease related to HTRA1 monoallelic variants. The mean age at onset of first clinical stroke was 51.25 years (range 41–64 years). The mean disease duration was 6.5 years (range 4–12). All individuals had recurrent strokes within the duration of follow‐up with a mean number of strokes per patient being 5.5 (range 2–12). Three individuals had leukoencephalopathy with brain stem involvement. Microhemorrhages were seen on brain MRI in three patients. HTRA1 monoallelic variants identified in our cohort were missense variants in three patients and a novel frameshift variation in one patient. Interestingly, two of these missense variants were previously reported in an autosomal recessive pattern of inheritance and here are associated with a dominant effect. Conclusions Clinicoradiologic characteristics of heterozygous HTRA1‐related CSVD may overlap with sporadic CSVD. Heterozygous HTRA1 variants can contribute to dominant or recessive disease mechanisms.

BackgroundCOVID-19 is caused by the SARS-CoV-2 virus and has strikingly heterogeneous clinical manifestations, with most individuals contracting mild disease but a substantial minority experiencing fulminant cardiopulmonary symptoms or death. The clinical covariates and the laboratory tests performed on a patient provide robust statistics to guide clinical treatment. Deep learning approaches on a data set of this nature enable patient stratification and provide methods to guide clinical treatment. ObjectiveHere, we report on the development and prospective validation of a state-of-the-art machine learning model to provide mortality prediction shortly after confirmation of SARS-CoV-2 infection in the Mayo Clinic patient population. MethodsWe retrospectively constructed one of the largest reported and most geographically diverse laboratory information system and electronic health record of COVID-19 data sets in the published literature, which included 11,807 patients residing in 41 states of the United States of America and treated at medical sites across 5 states in 3 time zones. Traditional machine learning models were evaluated independently as well as in a stacked learner approach by using AutoGluon, and various recurrent neural network architectures were considered. The traditional machine learning models were implemented using the AutoGluon-Tabular framework, whereas the recurrent neural networks utilized the TensorFlow Keras framework. We trained these models to operate solely using routine laboratory measurements and clinical covariates available within 72 hours of a patient’s first positive COVID-19 nucleic acid test result. ResultsThe GRU-D recurrent neural network achieved peak cross-validation performance with 0.938 (SE 0.004) as the area under the receiver operating characteristic (AUROC) curve. This model retained strong performance by reducing the follow-up time to 12 hours (0.916 [SE 0.005] AUROC), and the leave-one-out feature importance analysis indicated that the most independently valuable features were age, Charlson comorbidity index, minimum oxygen saturation, fibrinogen level, and serum iron level. In the prospective testing cohort, this model provided an AUROC of 0.901 and a statistically significant difference in survival (P

Abstract Background Use of the Genome Analysis Toolkit (GATK) continues to be the standard practice in genomic variant calling in both research and the clinic. Recently the toolkit has been rapidly evolving. Significant computational performance improvements have been introduced in GATK3.8 through collaboration with Intel in 2017. The first release of GATK4 in early 2018 revealed rewrites in the code base, as the stepping stone toward a Spark implementation. As the software continues to be a moving target for optimal deployment in highly productive environments, we present a detailed analysis of these improvements, to help the community stay abreast with changes in performance. Results We re-evaluated multiple options, such as threading, parallel garbage collection, I/O options and data-level parallelization. Additionally, we considered the trade-offs of using GATK3.8 and GATK4. We found optimized parameter values that reduce the time of executing the best practices variant calling procedure by 29.3% for GATK3.8 and 16.9% for GATK4. Further speedups can be accomplished by splitting data for parallel analysis, resulting in run time of only a few hours on whole human genome sequenced to the depth of 20X, for both versions of GATK. Nonetheless, GATK4 is already much more cost-effective than GATK3.8. Thanks to significant rewrites of the algorithms, the same analysis can be run largely in a single-threaded fashion, allowing users to process multiple samples on the same CPU. Conclusions In time-sensitive situations, when a patient has a critical or rapidly developing condition, it is useful to minimize the time to process a single sample. In such cases we recommend using GATK3.8 by splitting the sample into chunks and computing across multiple nodes. The resultant walltime will be nnn.4 hours at the cost of $41.60 on 4 c5.18xlarge instances of Amazon Cloud. For cost-effectiveness of routine analyses or for large population studies, it is useful to maximize the number of samples processed per unit time. Thus we recommend GATK4, running multiple samples on one node. The total walltime will be ∼34.1 hours on 40 samples, with 1.18 samples processed per hour at the cost of $2.60 per sample on c5.18xlarge instance of Amazon Cloud.

Genetic variants in ionotropic glutamate receptors have been implicated in neurodevelopmental disorders. Here, the authors report heterozygous de novo mutations in the GRIA2 gene in 28 individuals with intellectual disability and neurodevelopmental abnormalities associated with reduced Ca2+ transport and AMPAR currents.”

Key Clinical Message We report a patient with developmental delay, autism, epilepsy, macrocephaly, facial dysmorphism, gastrointestinal, and behavioral issues due to EXT2 compound heterozygous likely pathogenic variants. This case report expands the EXT2 gene mutation database and the clinical spectrum of patients with deficiencies in the heparan sulfate pathway.

Gestational trophoblastic disease (GTD) is a heterogeneous group of lesions arising from placental tissue. Epithelioid trophoblastic tumor (ETT), derived from chorionic-type trophoblast, is the rarest form of GTD with only approximately 130 cases described in the literature. Due to its morphologic mimicry of epithelioid smooth muscle tumors and carcinoma, ETT can be misdiagnosed. To date, molecular characterization of ETTs is lacking. Furthermore, ETT is difficult to treat when disease spreads beyond the uterus. Here using RNA-Seq analysis in a cohort of ETTs and other gestational trophoblastic lesions we describe the discovery of LPCAT1-TERT fusion transcripts that occur in ETTs and coincide with underlying genomic deletions. Through cell-growth assays we demonstrate that LPCAT1-TERT fusion proteins can positively modulate cell proliferation and therefore may represent future treatment targets. Furthermore, we demonstrate that TERT upregulation appears to be a characteristic of ETTs, even in the absence of LPCAT1-TERT fusions, and that it appears linked to copy number gains of chromosome 5. No evidence of TERT upregulation was identified in other trophoblastic lesions tested, including placental site trophoblastic tumors and placental site nodules, which are thought to be the benign chorionic-type trophoblast counterpart to ETT. These findings indicate that LPCAT1-TERT fusions and copy-number driven TERT activation may represent novel markers for ETT, with the potential to improve the diagnosis, treatment, and outcome for women with this rare form of GTD.

Abstract Background GNB1 encodes a subunit of a heterotrimeric G‐protein complex that transduces intracellular signaling cascades. Disruptions to the gene have previously been shown to be embryonic lethal in knockout mice and to cause complex neurodevelopmental disorders in humans. To date, the majority of variants associated with disease in humans have been missense variants in exons 5‐7. Methods Genetic sequencing was performed on two patients presenting with complex neurological phenotypes including intellectual disability, hypotonia, and in one patient seizures. Reported variants were assessed using RNA sequencing and functional BRET/BiFC assays. Results A splice variant reported in patient 1 was confirmed to cause usage of a cryptic splice site leading to a truncated protein product. Patient 2 was reported to have a truncating variant. BRET and BiFC assays of both patient variants confirmed both were deficient in inducing GPCR‐induced G protein activation due to lack of dimer formation with the Gγ subunit. Conclusion Here, we report two patients with functionally confirmed loss of function variants in GNB1 and neurodevelopmental phenotypes including intellectual disability, hypotonia, and seizures in one patient. These results suggest haploinsufficiency of GNB1 is a mechanism for neurodevelopmental disorders in humans.

Summary: Chondrogenic tumors involving the temporal bone are rare and typically arise spontaneously with unilateral presentation. Somatic IDH mutations are common in these tumors, but germline inheritance has not been documented to our knowledge. We describe familial chondrosarcoma, grade 1, of the mastoid with unilateral presentation in the mother and bilateral presentation in each of her two children. Each individual presented with headaches, facial paresis, and conductive hearing loss between the ages of 9–12. Exome sequencing of all three affected family members identified a shared germline heterozygous c.299G>A (p.Arg100Gln) missense variant in IDH1. The p.Arg100Gln variant has only rarely been observed as a somatic mutation in glial tumors, and previous in vitro experiments have shown that p.Arg100Gln produces small amounts of the oncometabolite D-2-hydroxyglutarate (D2HG). Biochemical testing in all three affected family members on urine and plasma was unable to detect increases in D2HG in these sample types. Due to insufficient tumor for methylation studies, we performed genome-wide methylation analysis of an IDH1 p.Arg100Gln mutant brain tumor from an unrelated individual to functionally evaluate this variant. These studies demonstrated a global hypermethylation phenotype consistent with other known isocitrate dehydrogenase (IDH) mutant brain tumors, suggesting that this variant has neomorphic activity despite low-level production of D2HG. The bones of the facial skeleton are formed by membranous ossification and we hypothesize that abnormal embryonic cartilage that rests within the suture lines may be involved in this tumor entity. Testing of additional individuals with similar presentations is needed to confirm this finding and clarify the associated phenotypes.

Abstract Background RNA polymerase III (Pol III)‐related disorders are autosomal recessive neurodegenerative disorders caused by variants in POLR3A or POLR3B. Recently, a novel phenotype of adult‐onset spastic ataxia was identified in individuals with the c.1909+22G>A POLR3A variant in compound heterozygosity. Methods Whole‐exome sequencing was performed in the proband and parents. Variants were confirmed by Sanger sequencing. RNA sequencing was performed to evaluate splicing implications. Results A 42‐year‐old female was evaluated for unexplained neurological findings with a slow progressive decline in gait and walking speed since adolescence. WES revealed a novel missense variant (c.3593A>C, p.Lys1198Arg) in exon 27 of POLR3A in compound heterozygosity with the c.1909+22G>A variant. Summary of previously reported clinical features from individuals with pathogenic biallelic alterations in POLR3A and adult‐onset phenotype is consistent with our findings. RNA analysis revealed c.3593A>G drives the production of four RNA transcript products each with different functional impacts. Conclusion The novel dual‐class c.3593A>C variant in POLR3A causes an amino acid substitution and complex disruption of splicing. Our report supports the need to investigate variants near splice junctions for proper interpretation. Current interpretation guidelines need to address best practices for inclusion of predicted or measured transcriptional disruption pending functional activity or reliable transcript abundance estimates.

Several recent studies have demonstrated the utility of RNA-Seq in the diagnosis of rare inherited disease. Diagnostic rates 35% higher than those previously achievable with DNA-Seq alone have been attained. These studies have primarily profiled gene expression and splicing defects, however, some have also shown that fusion transcripts are diagnostic or phenotypically relevant in patients with constitutional disorders. Fusion transcripts have traditionally been studied as oncogenic phenomena, with relevance only to cancer testing. Consequently, fusion detection algorithms were biased toward the detection of well-known oncogenic fusions, hindering their application to rare Mendelian genetic disease studies. A recent methodology published by the authors successfully tailored a traditional algorithm to the detection of pathogenic fusion events in inherited disease. A key mechanism of decreasing false positive or biologically benign events was comparison to a database of events detected in normal tissues. This approach is akin to population frequency-based filtering of genetic variants. It is predicated on the idea that pathogenic fusion transcripts are absent from normal tissue. We report on an analysis of RNA-Seq data from the genotype-tissue expression (GTEx) project in which known pathogenic fusions are computationally detected at low levels in normal tissues unassociated with the disease phenotype. Examples include archetypal cancer fusion transcripts, as well as fusions responsible for rare inherited disease. We consider potential explanations for the detectability of such transcripts and discuss the bearing such results have on the future profiling of genetic disease patients for pathogenic gene fusions.

Lysosomal diseases (LD) comprise a group of approximately 60 hereditary conditions caused by progressive accumulation of metabolites due to defects in lysosomal enzymes and degradation pathways, which lead to a wide range of clinical manifestations. The estimated combined incidence of LD is between 1 in 4000 to 1 in 13,000 live births, with recent data from pilot newborn screening studies showing even higher incidence. We aimed to determine the prevalence of the classical LD and other diseases caused by lysosome-related genes in our cohort of diagnostic odyssey patients. The Individualized Medicine Clinic at Mayo Clinic is increasingly utilizing whole exome sequencing (WES) to determine the genetic etiology of undiagnosed Mendelian disease. From September 2012 to April 2017, WES results from 350 patients with unexplained symptoms were reviewed. Disease-causing variants were identified in MYO6, CLN6, LRBA, KCTD7, and ARSB revealing a genetic diagnosis of a LD in 8 individuals from 5 families. Based on our findings, lysosome-related disorders may be collectively common, reaching up to 1.5% prevalence in a cohort of patients with undiagnosed diseases presenting to a genetics clinic.

In humans, de novo truncating variants in WASF1 (Wiskott–Aldrich syndrome protein family member 1) have been linked to presentations of moderate-to-profound intellectual disability (ID), autistic features, and epilepsy. Apart from one case series, there is limited information on the phenotypic spectrum and genetic landscape of WASF1-related neurodevelopmental disorder (NDD). In this report, we describe detailed clinical characteristics of six individuals with WASF1-related NDD. We demonstrate a broader spectrum of neurodevelopmental impairment including more mildly affected individuals. Further, we report new variant types, including a copy number variant (CNV), resulting in the partial deletion of WASF1 in monozygotic twins, and three missense variants, two of which alter the same residue, p.W161. This report adds further evidence that de novo variants in WASF1 cause an autosomal dominant NDD.

Abstract Background Multiple endocrine neoplasia type 1 (MEN1) is a hereditary cancer syndrome associated with several endocrine as well as non-endocrine tumors and is caused by mutations in the MEN1 gene. Primary hyperparathyroidism affects the majority of MEN1 individuals by age 50 years. Additionally, MEN1 mutations trigger familial isolated hyperparathyroidism. We describe a seemingly unaffected 76-year-old female who presented to our Genetics Clinic with a family history of primary hyperparathyroidism and the identification of a pathogenic MEN1 variant. Case Presentation The patient was a 76 year-old woman who appeared to be unaffected. She had a family history of a known MEN1 pathogenic variant. Molecular testing for the known MEN1 mutation c.1A > G, as well as, biochemical testing, MRI of the brain and abdomen were all performed using standard methods. Molecular testing revealed our patient possessed the MEN1 pathogenic variant previously identified in her two offspring. Physical exam revealed red facial papules with onset in her seventies, involving her cheeks, nose and upper lip. Formerly, she was diagnosed with rosacea by a dermatologist and noted no improvement with treatment. Clinically, these lesions appeared to be facial angiofibromas. Brain MRI was normal. However, an MRI of her abdomen revealed a 1.5 cm lesion at the tail of the pancreas with normal adrenal glands. Glucagon was mildly elevated and pancreatic polypeptide was nearly seven times the upper limit of the normal range. The patient underwent spleen sparing distal pancreatectomy and subsequent pathology was consistent with a well-differentiated pancreatic neuroendocrine tumor (pNET). Conclusions Age-related penetrance and variable expressivity are well documented in families with MEN1. It is thought that nearly all individuals with MEN1 manifest disease by age 40. We present a case of late-onset MEN1 in the absence of the most common feature, primary hyperparathyroidism, but with the presence of a pNET and cutaneous findings. This family expands the phenotype associated with the c.1A > G pathogenic variant and highlights the importance of providing comprehensive assessment of MEN1 mutation carriers in families that at first blush may appear to have isolated hyperparathyroidism.

Following publication of the original article [1], the author explained that Table 2 is displayed incorrectly. The correct Table 2 is given below. The original article has been corrected.

As reliable, efficient genome sequencing becomes ubiquitous, the need for similarly reliable and efficient variant calling becomes increasingly important. The Genome Analysis Toolkit (GATK), maintained by the Broad Institute, is currently the widely accepted standard for variant calling software. However, alternative solutions may provide faster variant calling without sacrificing accuracy. One such alternative is Sentieon DNASeq, a toolkit analogous to GATK but built on a highly optimized backend. We conducted an independent evaluation of the DNASeq single-sample variant calling pipeline in comparison to that of GATK. Our results support the near-identical accuracy of the two software packages, showcase optimal scalability and great speed from Sentieon, and describe computational performance considerations for the deployment of DNASeq.

BackgroundRNA sequencing has been proposed as a means of increasing diagnostic rates in studies of undiagnosed rare inherited disease. Recent studies have reported diagnostic improvements in the range of 7.5-35% by profiling splicing, gene expression quantification and allele specific expression. To-date however, no study has systematically assessed the presence of gene-fusion transcripts in cases of germline disease. Fusion transcripts are routinely identified in cancer studies and are increasingly recognized as having diagnostic, prognostic or therapeutic relevance. Isolated reports exist of fusion transcripts being detected in cases of developmental and neurological phenotypes, and thus, systematic application of fusion detection to germline conditions may further increase diagnostic rates. However, current fusion detection methods are unsuited to the investigation of germline disease due to performance biases arising from their development using tumor, cell-line or in-silico data.MethodsWe describe a tailored approach to fusion candidate identification and prioritization in a cohort of 47 undiagnosed, suspected inherited disease patients. We modify an existing fusion transcript detection algorithm by eliminating its cell line-derived filtering steps, and instead, prioritize candidates using a custom workflow that integrates genomic and transcriptomic sequence alignment, biological and technical annotations, customized categorization logic, and phenotypic prioritization.ResultsWe demonstrate that our approach to fusion transcript identification and prioritization detects genuine fusion events excluded by standard analyses and efficiently removes phenotypically unimportant candidates and false positive events, resulting in a reduced candidate list enriched for events with potential phenotypic relevance. We describe the successful genetic resolution of two previously undiagnosed disease cases through the detection of pathogenic fusion transcripts. Furthermore, we report the experimental validation of five additional cases of fusion transcripts with potential phenotypic relevance.ConclusionsThe approach we describe can be implemented to enable the detection of phenotypically relevant fusion transcripts in studies of rare inherited disease. Fusion transcript detection has the potential to increase diagnostic rates in rare inherited disease and should be included in RNA-based analytical pipelines aimed at genetic diagnosis.