Recombination mapping using Boolean logic and high-density SNP genotyping for exome sequence filtering

Clinical, genetic, and structural characterization of a novel TUBB4B tubulinopathy

Microtubules are cytoskeletal polymers of ⍺/β-tubulin heterodimers essential for a wide range of cellular processes. Pathogenic variations in microtubule-encoding genes (e.g., TUBB4B, which encodes the β-4B tubulin isotype) are responsible for a wide spectrum of cerebral malformations, collectively referred to as "tubulinopathies." The phenotypic manifestation of TUBB4B-associated tubulinopathy is Leber congenital amaurosis with early-onset deafness (LCAEOD), an autosomal dominant syndrome characterized by photoreceptor and cochlear cell loss; all known patients have pathogenic variations in amino acid R391. We present the clinical and molecular genetics findings of a 16-year-old female with a de novo missense variant in exon 1 of TUBB4B, c.32 A > G (p.Gln11Arg; Q11R). In addition to hearing loss and hyperopia without retinal abnormalities, our proband presented with two phenotypes of unknown genetic etiology, i.e., renal tubular Fanconi Syndrome (FS) and hypophosphatemic rickets (HR). The Q11R variant expands the genetic basis of early sensory hearing loss; its consequences with respect to microtubule structure are described. A mechanistic explanation for the FS and rickets, involving microtubule-mediated translocation of transporter proteins to and from the apical membrane of renal proximal tubular cells, is proposed.

Compound heterozygous KCTD7 variants in progressive myoclonus epilepsy

KCTD7 is a member of the potassium channel tetramerization domain-containing protein family and has been associated with progressive myoclonic epilepsy (PME), characterized by myoclonus, epilepsy, and neurological deterioration. Here we report four affected individuals from two unrelated families in which we identified KCTD7 compound heterozygous single nucleotide variants through exome sequencing. RNAseq was used to detect a non-annotated splicing junction created by a synonymous variant in the second family. Whole-cell patch-clamp analysis of neuroblastoma cells overexpressing the patients' variant alleles demonstrated aberrant potassium regulation. While all four patients experienced many of the common clinical features of PME, they also showed variable phenotypes not previously reported, including dysautonomia, brain pathology findings including a significantly reduced thalamus, and the lack of myoclonic seizures. To gain further insight into the pathogenesis of the disorder, zinc finger nucleases were used to generate kctd7 knockout zebrafish. Kctd7 homozygous mutants showed global dysregulation of gene expression and increased transcription of c-fos, which has previously been correlated with seizure activity in animal models. Together these findings expand the known phenotypic spectrum of KCTD7-associated PME, report a new animal model for future studies, and contribute valuable insights into the disease.

The undiagnosed diseases program: Approach to diagnosis

Undiagnosed and rare conditions are collectively common and affect millions of people worldwide. The NIH Undiagnosed Diseases Program (UDP) strives to achieve both a comprehensive diagnosis and a better understanding of the mechanisms of disease for many of these individuals. Through the careful review of records, a well-orchestrated inpatient evaluation, genomic sequencing and testing, and with the use of emerging strategies such as matchmaking programs, the UDP succeeds nearly 30 percent of the time for these highly selective cases. Although the UDP process is built on a unique set of resources, case examples demonstrate steps genetic professionals can take, in both clinical and research settings, to arrive at a diagnosis for their most challenging cases.

Cysteinyl-tRNA Synthetase Mutations Cause a Multi-System, Recessive Disease That Includes Microcephaly, Developmental Delay, and Brittle Hair and Nails

Aminoacyl-tRNA synthetases (ARSs) are essential enzymes responsible for charging tRNA molecules with cognate amino acids. Consistent with the essential function and ubiquitous expression of ARSs, mutations in 32 of the 37 ARS-encoding loci cause severe, early-onset recessive phenotypes. Previous genetic and functional data suggest a loss-of-function mechanism; however, our understanding of the allelic and locus heterogeneity of ARS-related disease is incomplete. Cysteinyl-tRNA synthetase (CARS) encodes the enzyme that charges tRNA^Cys with cysteine in the cytoplasm. To date, CARS variants have not been implicated in any human disease phenotype. Here, we report on four subjects from three families with complex syndromes that include microcephaly, developmental delay, and brittle hair and nails. Each affected person carries bi-allelic CARS variants: one individual is compound heterozygous for c.1138C>T (p.Gln380^∗) and c.1022G>A (p.Arg341His), two related individuals are compound heterozygous for c.1076C>T (p.Ser359Leu) and c.1199T>A (p.Leu400Gln), and one individual is homozygous for c.2061dup (p.Ser688Glnfs^∗2). Measurement of protein abundance, yeast complementation assays, and assessments of tRNA charging indicate that each CARS variant causes a loss-of-function effect. Compared to subjects with previously reported ARS-related diseases, individuals with bi-allelic CARS variants are unique in presenting with a brittle-hair-and-nail phenotype, which most likely reflects the high cysteine content in human keratins. In sum, our efforts implicate CARS variants in human inherited disease, expand the locus and clinical heterogeneity of ARS-related clinical phenotypes, and further support impaired tRNA charging as the primary mechanism of recessive ARS-related disease.

Novel mutations in CLN6 cause late-infantile neuronal ceroid lipofuscinosis without visual impairment in two unrelated patients

CLN6 is a transmembrane protein located in the endoplasmic reticulum that is involved in lysosomal acidification. Mutations in CLN6 cause late-infantile neuronal ceroid lipofuscinosis (LINCL), and teenage and adult onset NCL without visual impairment. Here we describe two pediatric patients with LINCL from unrelated families who were evaluated at the National Institutes of Health. Both children exhibited typical phenotypes associated with LINCL except that they lacked the expected visual impairment. Whole exome sequencing identified novel biallelic mutations in CLN6, i.e., c.218-220dupGGT (p.Trp73dup) and c.296A > G (p.Lys99Arg) in Proband 1 and homozygous c.723G > T (p.Met241Ile) in Proband 2. Expression analysis in dermal fibroblasts showed a small increase in CLN6 protein levels. Electron micrographs of these fibroblasts demonstrated large numbers of small membrane-bound vesicles, in addition to lipofuscin deposits. LysoTracker™ Red intensity was increased in fibroblasts from both patients. This study supports a role for CLN6 in lysosomal homeostasis, and highlights the importance of considering CLN6 mutations in the diagnosis of Batten Disease even in patients with normal vision.

Novel variants in SPTAN1 without epilepsy: An expansion of the phenotype

We describe two unrelated children with de novo variants in the non-erythrocytic alpha-II-spectrin (SPTAN1) gene who have hypoplastic brain structures, intellectual disability, and both fine and gross motor impairments. Using agnostic exome sequencing, we identified a nonsense variant creating a premature stop codon in exon 21 of SPTAN1, and in a second patient we identified an intronic substitution in SPTAN1 prior to exon 50 creating a new donor acceptor site. Neither of these variants has been described previously. Although some of these patients' features are consistent with the known SPTAN1 encephalopathy phenotype, these two children do not have epilepsy, in contrast to reports about nearly every other patient with heterozygous SPTAN1 variants and in all patients with a variant near the C-terminal coding region. Moreover, both children have abnormal thyroid function, which has not been previously reported in association with SPTAN1 variant. We present a detailed discussion of the clinical manifestations of these two unique SPTAN1 variants and provide evidence that both variants result in reduced mRNA expression despite different locations within the gene and clinical phenotypes. These findings expand the motor, cognitive, and behavioral spectrum of the SPTAN1-associated phenotype and invite speculation about underlying pathophysiologies.

Compound heterozygosity for loss-of-function GARS variants results in a multisystem developmental syndrome that includes severe growth retardation

Aminoacyl-tRNA synthetases (ARSs) are ubiquitously expressed enzymes that ligate amino acids onto tRNA molecules. Genes encoding ARSs have been implicated in myriad dominant and recessive disease phenotypes. Glycyl-tRNA synthetase (GARS) is a bifunctional ARS that charges tRNA^Gly in the cytoplasm and mitochondria. GARS variants have been associated with dominant Charcot-Marie-Tooth disease but have not been convincingly implicated in recessive phenotypes. Here, we describe a patient from the NIH Undiagnosed Diseases Program with a multisystem, developmental phenotype. Whole-exome sequence analysis revealed that the patient is compound heterozygous for one frameshift (p.Glu83Ilefs*6) and one missense (p.Arg310Gln) GARS variant. Using in vitro and in vivo functional studies, we show that both GARS variants cause a loss-of-function effect: the frameshift variant results in depleted protein levels and the missense variant reduces GARS tRNA charging activity. In support of GARS variant pathogenicity, our patient shows striking phenotypic overlap with other patients having ARS-related recessive diseases, including features associated with variants in both cytoplasmic and mitochondrial ARSs; this observation is consistent with the essential function of GARS in both cellular locations. In summary, our clinical, genetic, and functional analyses expand the phenotypic spectrum associated with GARS variants.

Explorations to improve the completeness of exome sequencing

BACKGROUND

Exome sequencing has advanced to clinical practice and proven useful for obtaining molecular diagnoses in rare diseases. In approximately 75 % of cases, however, a clinical exome study does not produce a definitive molecular diagnosis. These residual cases comprise a new diagnostic challenge for the genetics community. The Undiagnosed Diseases Program of the National Institutes of Health routinely utilizes exome sequencing for refractory clinical cases. Our preliminary data suggest that disease-causing variants may be missed by current standard-of-care clinical exome analysis. Such false negatives reflect limitations in experimental design, technical performance, and data analysis.

RESULTS

We present examples from our datasets to quantify the analytical performance associated with current practices, and explore strategies to improve the completeness of data analysis. In particular, we focus on patient ascertainment, exome capture, inclusion of intronic variants, and evaluation of medium-sized structural variants.

CONCLUSIONS

The strategies we present may recover previously-missed, disease causing variants in second-pass exome analysis. Understanding the limitations of the current clinical exome search space provides a rational basis to improve methods for disease variant detection using genome-scale sequencing techniques.

Phenotypic evolution of UNC80 loss of function

Failure to thrive arises as a complication of a heterogeneous group of disorders. We describe two female siblings with spastic paraplegia and global developmental delay but also, atypically for the HSPs, poor weight gain classified as failure to thrive. After extensive clinical and biochemical investigations failed to identify the etiology, we used exome sequencing to identify biallelic UNC80 mutations (NM_032504.1:c.[3983-3_3994delinsA];[2431C>T]. The paternally inherited NM_032504.1:c.3983-3_3994delinsA is predicted to encode p.Ser1328Argfs*19 and the maternally inherited NM_032504.1:c.2431C>T is predicted to encode p.Arg811*. No UNC80 mRNA was detectable in patient cultured skin fibroblasts, suggesting UNC80 loss of function by nonsense mediated mRNA decay. Further supporting the UNC80 mutations as causative of these siblings' disorder, biallelic mutations in UNC80 have recently been described among individuals with an overlapping phenotype. This report expands the disease spectrum associated with UNC80 mutations. © 2016 Wiley Periodicals, Inc.

Recurrent Mutations in the Basic Domain of TWIST2 Cause Ablepharon Macrostomia and Barber-Say Syndromes

Ablepharon macrostomia syndrome (AMS) and Barber-Say syndrome (BSS) are rare congenital ectodermal dysplasias characterized by similar clinical features. To establish the genetic basis of AMS and BSS, we performed extensive clinical phenotyping, whole exome and candidate gene sequencing, and functional validations. We identified a recurrent de novo mutation in TWIST2 in seven independent AMS-affected families, as well as another recurrent de novo mutation affecting the same amino acid in ten independent BSS-affected families. Moreover, a genotype-phenotype correlation was observed, because the two syndromes differed based solely upon the nature of the substituting amino acid: a lysine at TWIST2 residue 75 resulted in AMS, whereas a glutamine or alanine yielded BSS. TWIST2 encodes a basic helix-loop-helix transcription factor that regulates the development of mesenchymal tissues. All identified mutations fell in the basic domain of TWIST2 and altered the DNA-binding pattern of Flag-TWIST2 in HeLa cells. Comparison of wild-type and mutant TWIST2 expressed in zebrafish identified abnormal developmental phenotypes and widespread transcriptome changes. Our results suggest that autosomal-dominant TWIST2 mutations cause AMS or BSS by inducing protean effects on the transcription factor's DNA binding.

MED23-associated intellectual disability in a non-consanguineous family

Intellectual disability (ID) is a heterogeneous condition arising from a variety of environmental and genetic factors. Among these causes are defects in transcriptional regulators. Herein, we report on two brothers in a nonconsanguineous family with novel compound heterozygous, disease-segregating mutations (NM_015979.3: [3656A > G];[4006C > T], NP_057063.2: [H1219R];[R1336X]) in MED23. This gene encodes a subunit of the Mediator complex that modulates the expression of RNA polymerase II-dependent genes. These brothers, who had profound ID, spasticity, congenital heart disease, brain abnormalities, and atypical electroencephalography, represent the first case of MED23-associated ID in a non-consanguineous family. They also expand upon the clinical features previously reported for mutations in this gene.

Impaired osteoblast and osteoclast function characterize the osteoporosis of Snyder - Robinson syndrome

Background

Snyder-Robinson Syndrome (SRS) is an X-linked intellectual disability disorder also characterized by osteoporosis, scoliosis, and dysmorphic facial features. It is caused by mutations in SMS, a ubiquitously expressed gene encoding the polyamine biosynthetic enzyme spermine synthase. We hypothesized that the tissue specificity of SRS arises from differential sensitivity to spermidine toxicity or spermine deficiency.

Methods

We performed detailed clinical, endocrine, histopathologic, and morphometric studies on two affected brothers with a spermine synthase loss of function mutation (NM_004595.4:c.443A > G, p.Gln148Arg). We also measured spermine and spermidine levels in cultured human bone marrow stromal cells (hBMSCs) and fibroblasts using the Biochrom 30 polyamine protocol and assessed the osteogenic potential of hBMSCs.

Results

In addition to the known tissue-specific features of SRS, the propositi manifested retinal pigmentary changes, recurrent episodes of hyper- and hypoglycemia, nephrocalcinosis, renal cysts, and frequent respiratory infections. Bone histopathology and morphometry identified a profound depletion of osteoblasts and osteoclasts, absence of a trabecular meshwork, a low bone volume and a thin cortex. Comparison of cultured fibroblasts from affected and unaffected individuals showed relatively small changes in polyamine content, whereas comparison of cultured osteoblasts identified marked differences in spermidine and spermine content. Osteogenic differentiation of the SRS-derived hBMSCs identified a severe deficiency of calcium phosphate mineralization.

Conclusions

Our findings support the hypothesis that cell specific alterations in polyamine metabolism contribute to the tissue specificity of SRS features, and that the low bone density arises from a failure of mineralization.

Electronic supplementary material

The online version of this article (doi:10.1186/s13023-015-0235-8) contains supplementary material, which is available to authorized users.

Three rare diseases in one Sib pair: RAI1, PCK1, GRIN2B mutations associated with Smith-Magenis Syndrome, cytosolic PEPCK deficiency and NMDA receptor glutamate insensitivity

The National Institutes of Health Undiagnosed Diseases Program evaluates patients for whom no diagnosis has been discovered despite a comprehensive diagnostic workup. Failure to diagnose a condition may arise from the mutation of genes previously unassociated with disease. However, we hypothesized that this could also co-occur with multiple genetic disorders. Demonstrating a complex syndrome caused by multiple disorders, we report two siblings manifesting both similar and disparate signs and symptoms. They shared a history of episodes of hypoglycemia and lactic acidosis, but had differing exam findings and developmental courses. Clinical acumen and exome sequencing combined with biochemical and functional studies identified three genetic conditions. One sibling had Smith-Magenis Syndrome and a nonsense mutation in the RAI1 gene. The second sibling had a de novo mutation in GRIN2B, which resulted in markedly reduced glutamate potency of the encoded receptor. Both siblings had a protein-destabilizing homozygous mutation in PCK1, which encodes the cytosolic isoform of phosphoenolpyruvate carboxykinase (PEPCK-C). In summary, we present the first clinically-characterized mutation of PCK1 and demonstrate that complex medical disorders can represent the co-occurrence of multiple diseases.

An apparent homozygous deletion in maltase-glucoamylase, a lesson in the evolution of SNP arrays

Single nucleotide polymorphism (SNP) arrays possess clinical potential due to their high throughput capacity, sensitivity and versatility. We used such an array to perform a genome-wide SNP analysis of a patient with a multi-system undiagnosed disease involving peripheral neuropathies and food intolerances. The patient had a homozygous deletion within the gene encoding maltase-glucoamylase (MGAM), an intestinal starch digestion enzyme, predicting absence of enzyme activity and potential starch indigestion. We then performed validation testing using a functional MGAM analysis that involved starch ingestion followed by measuring blood glucose and insulin levels as well as hydrogen breath levels. Gastrointestinal tissue was also obtained via endoscopy and immunohistochemical staining for intestinal MGAM was performed. Our results strongly suggest the presence and functioning of MGAM which disproved deficiency predictions based on SNP array analysis findings, classifying the deletion as a functional polymorphism. This study highlights a current clinical limitation of SNP arrays, i.e., distinguishing deleterious genomic alterations from misleading functional polymorphisms. We conclude that novel findings from SNP arrays should be clinically validated and published.

Analysis of DNA sequence variants detected by high-throughput sequencing

The Undiagnosed Diseases Program at the National Institutes of Health uses high-throughput sequencing (HTS) to diagnose rare and novel diseases. HTS techniques generate large numbers of DNA sequence variants, which must be analyzed and filtered to find candidates for disease causation. Despite the publication of an increasing number of successful exome-based projects, there has been little formal discussion of the analytic steps applied to HTS variant lists. We present the results of our experience with over 30 families for whom HTS sequencing was used in an attempt to find clinical diagnoses. For each family, exome sequence was augmented with high-density SNP-array data. We present a discussion of the theory and practical application of each analytic step and provide example data to illustrate our approach. The article is designed to provide an analytic roadmap for variant analysis, thereby enabling a wide range of researchers and clinical genetics practitioners to perform direct analysis of HTS data for their patients and projects.