RNA Sequencing as a Diagnostic Tool Ketkar, Shamika; Burrage, Lindsay C; Lee, Brendan
JAMA : the journal of the American Medical Association,
01/2023, Volume:
329, Issue:
1
Journal Article
BACKGROUNDTranscriptome sequencing (RNA-seq) improves diagnostic rates in individuals with suspected Mendelian conditions to varying degrees, primarily by directing the prioritization of candidate ...DNA variants identified on exome or genome sequencing (ES/GS). Here we implemented an RNA-seq-guided method to diagnose individuals across a wide range of ages and clinical phenotypes.METHODSOne hundred fifteen undiagnosed adult and pediatric patients with diverse phenotypes and 67 family members (182 total individuals) underwent RNA-seq from whole blood and skin fibroblasts at the Baylor College of Medicine (BCM) Undiagnosed Diseases Network clinical site from 2014 to 2020. We implemented a workflow to detect outliers in gene expression and splicing for cases that remained undiagnosed despite standard genomic and transcriptomic analysis.RESULTSThe transcriptome-directed approach resulted in a diagnostic rate of 12% across the entire cohort, or 17% after excluding cases solved on ES/GS alone. Newly diagnosed conditions included Koolen-de Vries syndrome (KANSL1), Renpenning syndrome (PQBP1), TBCK-associated encephalopathy, NSD2- and CLTC-related intellectual disability, and others, all with negative conventional genomic testing, including ES and chromosomal microarray (CMA). Skin fibroblasts exhibited higher and more consistent expression of clinically relevant genes than whole blood. In solved cases with RNA-seq from both tissues, the causative defect was missed in blood in half the cases but none from fibroblasts.CONCLUSIONSFor our cohort of undiagnosed individuals with suspected Mendelian conditions, transcriptome-directed genomic analysis facilitated diagnoses, primarily through the identification of variants missed on ES and CMA.TRIAL REGISTRATIONNot applicable.FUNDINGNIH Common Fund, BCM Intellectual and Developmental Disabilities Research Center, Eunice Kennedy Shriver National Institute of Child Health & Human Development.
Branched-chain amino acid (BCAA) metabolism plays a central role in the pathophysiology of both rare inborn errors of metabolism and the more common multifactorial diseases. Although deficiency of ...the branched-chain ketoacid dehydrogenase (BCKDC) and associated elevations in the BCAAs and their ketoacids have been recognized as the cause of maple syrup urine disease (MSUD) for decades, treatment options for this disorder have been limited to dietary interventions. In recent years, the discovery of improved leucine tolerance after liver transplantation has resulted in a new therapeutic strategy for this disorder. Likewise, targeting the regulation of the BCKDC activity may be an alternative potential treatment strategy for MSUD. The regulation of the BCKDC by the branched-chain ketoacid dehydrogenase kinase has also been implicated in a new inborn error of metabolism characterized by autism, intellectual disability and seizures. Finally, there is a growing body of literature implicating BCAA metabolism in more common disorders such as the metabolic syndrome, cancer and hepatic disease. This review surveys the knowledge acquired on the topic over the past 50 years and focuses on recent developments in the field of BCAA metabolism.
While congenital malformations and genetic diseases are a leading cause of early infant death, to our knowledge, the contribution of single-gene disorders in this group is undetermined.
To determine ...the diagnostic yield and use of clinical exome sequencing in critically ill infants.
Clinical exome sequencing was performed for 278 unrelated infants within the first 100 days of life who were admitted to Texas Children's Hospital in Houston, Texas, during a 5-year period between December 2011 and January 2017. Exome sequencing types included proband exome, trio exome, and critical trio exome, a rapid genomic assay for seriously ill infants.
Indications for testing, diagnostic yield of clinical exome sequencing, turnaround time, molecular findings, patient age at diagnosis, and effect on medical management among a group of critically ill infants who were suspected to have genetic disorders.
The mean (SEM) age for infants participating in the study was 28.5 (1.7) days; of these, the mean (SEM) age was 29.0 (2.2) days for infants undergoing proband exome sequencing, 31.5 (3.9) days for trio exome, and 22.7 (3.9) days for critical trio exome. Clinical indications for exome sequencing included a range of medical concerns. Overall, a molecular diagnosis was achieved in 102 infants (36.7%) by clinical exome sequencing, with relatively low yield for cardiovascular abnormalities. The diagnosis affected medical management for 53 infants (52.0%) and had a substantial effect on informed redirection of care, initiation of new subspecialist care, medication/dietary modifications, and furthering life-saving procedures in select patients. Critical trio exome sequencing revealed a molecular diagnosis in 32 of 63 infants (50.8%) at a mean (SEM) of 33.1 (5.6) days of life with a mean (SEM) turnaround time of 13.0 (0.4) days. Clinical care was altered by the diagnosis in 23 of 32 patients (71.9%). The diagnostic yield, patient age at diagnosis, and medical effect in the group that underwent critical trio exome sequencing were significantly different compared with the group who underwent regular exome testing. For deceased infants (n = 81), genetic disorders were molecularly diagnosed in 39 (48.1%) by exome sequencing, with implications for recurrence risk counseling.
Exome sequencing is a powerful tool for the diagnostic evaluation of critically ill infants with suspected monogenic disorders in the neonatal and pediatric intensive care units and its use has a notable effect on clinical decision making.
Untargeted metabolomics is a global molecular profiling technology that can be used to screen for inborn errors of metabolism (IEMs). Metabolite perturbations are evaluated based on current knowledge ...of specific metabolic pathway deficiencies, a manual diagnostic process that is qualitative, has limited scalability, and is not equipped to learn from accumulating clinical data. Our purpose was to improve upon manual diagnosis of IEMs in the clinic by developing novel computational methods for analyzing untargeted metabolomics data. We employed CTD, an automated computational diagnostic method that "connects the dots" between metabolite perturbations observed in individual metabolomics profiling data and modules identified in disease-specific metabolite co-perturbation networks learned from prior profiling data. We also extended CTD to calculate distances between any two individuals (CTDncd) and between an individual and a disease state (CTDdm), to provide additional network-quantified predictors for use in diagnosis. We show that across 539 plasma samples, CTD-based network-quantified measures can reproduce accurate diagnosis of 16 different IEMs, including adenylosuccinase deficiency, argininemia, argininosuccinic aciduria, aromatic L-amino acid decarboxylase deficiency, cerebral creatine deficiency syndrome type 2, citrullinemia, cobalamin biosynthesis defect, GABA-transaminase deficiency, glutaric acidemia type 1, maple syrup urine disease, methylmalonic aciduria, ornithine transcarbamylase deficiency, phenylketonuria, propionic acidemia, rhizomelic chondrodysplasia punctata, and the Zellweger spectrum disorders. Our approach can be used to supplement information from biochemical pathways and has the potential to significantly enhance the interpretation of variants of uncertain significance uncovered by exome sequencing. CTD, CTDdm, and CTDncd can serve as an essential toolset for biological interpretation of untargeted metabolomics data that overcomes limitations associated with manual diagnosis to assist diagnosticians in clinical decision-making. By automating and quantifying the interpretation of perturbation patterns, CTD can improve the speed and confidence by which clinical laboratory directors make diagnostic and treatment decisions, while automatically improving performance with new case data.
Previous studies have shown that nitric oxide (NO) supplements may prevent bone loss and fractures in preclinical models of estrogen deficiency. However, the mechanisms by which NO modulates bone ...anabolism remain largely unclear. Argininosuccinate lyase (ASL) is the only mammalian enzyme capable of synthesizing arginine, the sole precursor for nitric oxide synthase-dependent (NOS-dependent) NO synthesis. Moreover, ASL is also required for channeling extracellular arginine to NOS for NO production. ASL deficiency (ASLD) is thus a model to study cell-autonomous, NOS-dependent NO deficiency. Here, we report that loss of ASL led to decreased NO production and impairment of osteoblast differentiation. Mechanistically, the bone phenotype was at least in part driven by the loss of NO-mediated activation of the glycolysis pathway in osteoblasts that led to decreased osteoblast differentiation and function. Heterozygous deletion of caveolin 1, a negative regulator of NO synthesis, restored NO production, osteoblast differentiation, glycolysis, and bone mass in a hypomorphic mouse model of ASLD. The translational significance of these preclinical studies was further reiterated by studies conducted in induced pluripotent stem cells from an individual with ASLD. Taken together, our findings suggest that ASLD is a unique genetic model for studying NO-dependent osteoblast function and that the NO/glycolysis pathway may be a new target to modulate bone anabolism.
Objectives To test whether follow-up testing for very long-chain acyl-CoA dehydrogenase (VLCAD) deficiency uncovers a diagnosis in patients with elevations of C14:1 and C14:2 plasma acylcarnitines ...after a controlled fasting study performed for clinically suspected hypoglycemia and to compare the acylcarnitine profiles from fasted patients without VLCAD deficiency vs patients with known VLCAD deficiency to determine whether metabolite testing distinguishes these groups. Study design We performed a retrospective chart review and identified 17 patients with elevated C14:1 and C14:2 plasma acylcarnitine levels after a controlled fast and with testing for VLCAD deficiency ( ACADVL sequencing or fibroblast fatty acid oxidation studies). The follow-up testing in all patients was inconsistent with a diagnosis of VLCAD deficiency. We compared the plasma acylcarnitine profiles from these fasted patients vs patients with VLCAD deficiency. Results C14:1/C12:1 was significantly lower ( P < .001) in fasted patients vs patients with VLCAD deficiency. Metabolomics analysis performed in 2 fasted patients and 1 patient with VLCAD deficiency demonstrated evidence for up-regulated lipolysis and β-oxidation in the fasted state. Conclusions Elevations of plasma C14:1 and C14:2 acylcarnitines appear to be a physiologic result of lipolysis that occurs with fasting. Both metabolomics analysis and/or C14:1/C12:1 may distinguish C14:1 elevations from physiologic fasting-induced lipolysis vs VLCAD deficiency.
The genetic architecture of complex traits underlying physiology and disease in most organisms remains elusive. We still know little about the number of genes that underlie these traits, the ...magnitude of their effects, or the extent to which they interact. Chromosome substitution strains (CSSs) enable statistically powerful studies based on testing engineered inbred strains that have single, unique, and nonoverlapping genetic differences, thereby providing measures of phenotypic effects that are attributable to individual chromosomes. Here, we report a study of phenotypic effects and gene interactions for 90 blood, bone, and metabolic traits in a mouse CSS panel and 54 traits in a rat CSS panel. Two key observations emerge about the genetic architecture of these traits. First, the traits tend to be highly polygenic: across the genome, many individual chromosome substitutions each had significant phenotypic effects and, within each of the chromosomes studied, multiple distinct loci were found. Second, strong epistasis was found among the individual chromosomes. Specifically, individual chromosome substitutions often conferred surprisingly large effects (often a substantial fraction of the entire phenotypic difference between the parental strains), with the result that the sum of these individual effects often dramatically exceeded the difference between the parental strains. We suggest that strong, pervasive epistasis may reflect the presence of several phenotypically-buffered physiological states. These results have implications for identification of complex trait genes, developmental and physiological studies of phenotypic variation, and opportunities to engineer phenotypic outcomes in complex biological systems.
Kinesin motor proteins transport intracellular cargo, including mRNA, proteins, and organelles. Pathogenic variants in kinesin-related genes have been implicated in neurodevelopmental disorders and ...skeletal dysplasias. We identified de novo, heterozygous variants in KIF5B, encoding a kinesin-1 subunit, in four individuals with osteogenesis imperfecta. The variants cluster within the highly conserved kinesin motor domain and are predicted to interfere with nucleotide binding, although the mechanistic consequences on cell signaling and function are unknown.
To understand the in vivo genetic mechanism of KIF5B variants, we modeled the p.Thr87Ile variant that was found in two patients in the C. elegans ortholog, unc-116, at the corresponding position (Thr90Ile) by CRISPR/Cas9 editing and performed functional analysis. Next, we studied the cellular and molecular consequences of the recurrent p.Thr87Ile variant by microscopy, RNA and protein analysis in NIH3T3 cells, primary human fibroblasts and bone biopsy.
C. elegans heterozygous for the unc-116 Thr90Ile variant displayed abnormal body length and motility phenotypes that were suppressed by additional copies of the wild type allele, consistent with a dominant negative mechanism. Time-lapse imaging of GFP-tagged mitochondria showed defective mitochondria transport in unc-116 Thr90Ile neurons providing strong evidence for disrupted kinesin motor function. Microscopy studies in human cells showed dilated endoplasmic reticulum, multiple intracellular vacuoles, and abnormal distribution of the Golgi complex, supporting an intracellular trafficking defect. RNA sequencing, proteomic analysis, and bone immunohistochemistry demonstrated down regulation of the mTOR signaling pathway that was partially rescued with leucine supplementation in patient cells.
We report dominant negative variants in the KIF5B kinesin motor domain in individuals with osteogenesis imperfecta. This study expands the spectrum of kinesin-related disorders and identifies dysregulated signaling targets for KIF5B in skeletal development.