Transcriptomic profiling of TK2 deficient human skeletal muscle suggests a role for the p53 signalling pathway and identifies growth and differentiation factor-15 as a potential novel biomarker for mitochondrial myopathies

Background Mutations in the gene encoding thymidine kinase 2 (TK2) result in the myopathic form of mitochondrial DNA depletion syndrome which is a mitochondrial encephalomyopathy presenting in children. In order to unveil some of the mechanisms involved in this pathology and to identify potential biomarkers and therapeutic targets we have investigated the gene expression profile of human skeletal muscle deficient for TK2 using cDNA microarrays. Results We have analysed the whole transcriptome of skeletal muscle from patients with TK2 mutations and compared it to normal muscle and to muscle from patients with other mitochondrial myopathies. We have identified a set of over 700 genes which are differentially expressed in TK2 deficient muscle. Bioinformatics analysis reveals important changes in muscle metabolism, in particular, in glucose and glycogen utilisation, and activation of the starvation response which affects aminoacid and lipid metabolism. We have identified those transcriptional regulators which are likely to be responsible for the observed changes in gene expression. Conclusion Our data point towards the tumor suppressor p53 as the regulator at the centre of a network of genes which are responsible for a coordinated response to TK2 mutations which involves inflammation, activation of muscle cell death by apoptosis and induction of growth and differentiation factor 15 (GDF-15) in muscle and serum. We propose that GDF-15 may represent a potential novel biomarker for mitochondrial dysfunction although further studies are required

A homozygous MED11 C-terminal variant causes a lethal neurodegenerative disease

Purpose: The mediator (MED) multisubunit-complex modulates the activity of the transcriptional machinery, and genetic defects in different MED subunits (17, 20, 27) have been implicated in neurologic diseases. In this study, we identified a recurrent homozygous variant in MED11 (c.325C>T; p.Arg109Ter) in 7 affected individuals from 5 unrelated families. Methods: To investigate the genetic cause of the disease, exome or genome sequencing were performed in 5 unrelated families identified via different research networks and Matchmaker Exchange. Deep clinical and brain imaging evaluations were performed by clinical pediatric neurologists and neuroradiologists. The functional effect of the candidate variant on both MED11 RNA and protein was assessed using reverse transcriptase polymerase chain reaction and western blotting using fibroblast cell lines derived from 1 affected individual and controls and through computational approaches. Knockouts in zebrafish were generated using clustered regularly interspaced short palindromic repeats/Cas9. Results: The disease was characterized by microcephaly, profound neurodevelopmental impairment, exaggerated startle response, myoclonic seizures, progressive widespread neurodegeneration, and premature death. Functional studies on patient-derived fibroblasts did not show a loss of protein function but rather disruption of the C-terminal of MED11, likely impairing binding to other MED subunits. A zebrafish knockout model recapitulates key clinical phenotypes. Conclusion: Loss of the C-terminal of MED subunit 11 may affect its binding efficiency to other MED subunits, thus implicating the MED-complex stability in brain development and neurodegeneration

Bi-allelic genetic variants in the translational GTPases GTPBP1 and GTPBP2 cause a distinct identical neurodevelopmental syndrome

The homologous genes GTPBP1 and GTPBP2 encode GTP-binding proteins 1 and 2, which are involved in ribosomal homeostasis. Pathogenic variants in GTPBP2 were recently shown to be an ultra-rare cause of neurodegenerative or neurodevelopmental disorders (NDDs). Until now, no human phenotype has been linked to GTPBP1. Here, we describe individuals carrying bi-allelic GTPBP1 variants that display an identical phenotype with GTPBP2 and characterize the overall spectrum of GTP-binding protein (1/2)-related disorders. In this study, 20 individuals from 16 families with distinct NDDs and syndromic facial features were investigated by whole-exome (WES) or whole-genome (WGS) sequencing. To assess the functional impact of the identified genetic variants, semi-quantitative PCR, western blot, and ribosome profiling assays were performed in fibroblasts from affected individuals. We also investigated the effect of reducing expression of CG2017, an ortholog of human GTPBP1/2, in the fruit fly Drosophila melanogaster. Individuals with bi-allelic GTPBP1 or GTPBP2 variants presented with microcephaly, profound neurodevelopmental impairment, pathognomonic craniofacial features, and ectodermal defects. Abnormal vision and/or hearing, progressive spasticity, choreoathetoid movements, refractory epilepsy, and brain atrophy were part of the core phenotype of this syndrome. Cell line studies identified a loss-of-function (LoF) impact of the disease-associated variants but no significant abnormalities on ribosome profiling. Reduced expression of CG2017 isoforms was associated with locomotor impairment in Drosophila. In conclusion, bi-allelic GTPBP1 and GTPBP2 LoF variants cause an identical, distinct neurodevelopmental syndrome. Mutant CG2017 knockout flies display motor impairment, highlighting the conserved role for GTP-binding proteins in CNS development across species

Twist exome capture allows for lower average sequence coverage in clinical exome sequencing

Background Exome and genome sequencing are the predominant techniques in the diagnosis and research of genetic disorders. Sufficient, uniform and reproducible/consistent sequence coverage is a main determinant for the sensitivity to detect single-nucleotide (SNVs) and copy number variants (CNVs). Here we compared the ability to obtain comprehensive exome coverage for recent exome capture kits and genome sequencing techniques. Results We compared three different widely used enrichment kits (Agilent SureSelect Human All Exon V5, Agilent SureSelect Human All Exon V7 and Twist Bioscience) as well as short-read and long-read WGS. We show that the Twist exome capture significantly improves complete coverage and coverage uniformity across coding regions compared to other exome capture kits. Twist performance is comparable to that of both short- and long-read whole genome sequencing. Additionally, we show that even at a reduced average coverage of 70× there is only minimal loss in sensitivity for SNV and CNV detection. Conclusion We conclude that exome sequencing with Twist represents a significant improvement and could be performed at lower sequence coverage compared to other exome capture techniques

A Solve-RD ClinVar-based reanalysis of 1522 index cases from ERN-ITHACA reveals common pitfalls and misinterpretations in exome sequencing

Purpose Within the Solve-RD project (https://solve-rd.eu/), the European Reference Network for Intellectual disability, TeleHealth, Autism and Congenital Anomalies aimed to investigate whether a reanalysis of exomes from unsolved cases based on ClinVar annotations could establish additional diagnoses. We present the results of the “ClinVar low-hanging fruit” reanalysis, reasons for the failure of previous analyses, and lessons learned. Methods Data from the first 3576 exomes (1522 probands and 2054 relatives) collected from European Reference Network for Intellectual disability, TeleHealth, Autism and Congenital Anomalies was reanalyzed by the Solve-RD consortium by evaluating for the presence of single-nucleotide variant, and small insertions and deletions already reported as (likely) pathogenic in ClinVar. Variants were filtered according to frequency, genotype, and mode of inheritance and reinterpreted. Results We identified causal variants in 59 cases (3.9%), 50 of them also raised by other approaches and 9 leading to new diagnoses, highlighting interpretation challenges: variants in genes not known to be involved in human disease at the time of the first analysis, misleading genotypes, or variants undetected by local pipelines (variants in off-target regions, low quality filters, low allelic balance, or high frequency). Conclusion The “ClinVar low-hanging fruit” analysis represents an effective, fast, and easy approach to recover causal variants from exome sequencing data, herewith contributing to the reduction of the diagnostic deadlock

Congenital myasthenic syndrome caused by novel COL13A1 mutations.

Collagen XIII is a non-fibrillar transmembrane collagen which has been long recognized for its critical role in synaptic maturation of the neuromuscular junction. More recently, biallelic COL13A1 loss-of-function mutations were identified in three patients with congenital myasthenic syndrome (CMS), a rare inherited condition with defective neuromuscular transmission, causing abnormal fatigability and fluctuating muscle weakness and often successfully treated with acetylcholinesterase inhibitors. Here we report six additional CMS patients from three unrelated families with previously unreported homozygous COL13A1 loss-of-function mutations (p.Tyr216*, p.Glu543fs and p.Thr629fs). The phenotype of our cases was similar to the previously reported patients including respiratory distress and severe dysphagia at birth that often resolved or improved in the first days or weeks of life. All individuals had prominent eyelid ptosis with only minor ophthalmoparesis as well as generalized muscle weakness, predominantly affecting facial, bulbar, respiratory and axial muscles. Response to acetylcholinesterase inhibitor treatment was generally negative while salbutamol proved beneficial. Our data further support the causality of COL13A1 variants for CMS and suggest that this type of CMS might be clinically homogenous and requires alternative pharmacological therapy

Survival among children with “Lethal” congenital contracture syndrome 11 caused by novel mutations in the gliomedin gene (GLDN).

Wiley Periodicals, Inc. Biallelic GLDN mutations have recently been identified among infants with lethal congenital contracture syndrome 11 (LCCS11). GLDN encodes gliomedin, a protein required for the formation of the nodes of Ranvier and development of the human peripheral nervous system. We report six infants and children from four unrelated families with biallelic GLDN mutations, four of whom survived beyond the neonatal period into infancy, childhood, and late adolescence with intensive care and chronic respiratory and nutritional support. Our findings expand the genotypic and phenotypic spectrum of LCCS11 and demonstrate that the condition may not necessarily be lethal in the neonatal period