PatientMatcher: A customizable Python-based open-source tool for matching undiagnosed rare disease patients via the Matchmaker Exchange network.

The amount of data available from genomic medicine has revolutionized the approach to identify the determinants underlying many rare diseases. The task of confirming a genotype-phenotype causality for a patient affected with a rare genetic disease is often challenging. In this context, the establishment of the Matchmaker Exchange (MME) network has assumed a pivotal role in bridging heterogeneous patient information stored on different medical and research servers. MME has made it possible to solve rare disease cases by "matching" the genotypic and phenotypic characteristics of a patient of interest with patient data available at other clinical facilities participating in the network. Here, we present PatientMatcher (https://github.com/Clinical-Genomics/patientMatcher), an open-source Python and MongoDB-based software solution developed by Clinical Genomics facility at the Science for Life Laboratory in Stockholm. PatientMatcher is designed as a standalone MME server, but can easily communicate via REST API with external applications managing genetic analyses and patient data. The MME node is being implemented in clinical routine in collaboration with the Genomic Medicine Center Karolinska at the Karolinska University Hospital. PatientMatcher is written to implement the MME API and provides several customizable settings, including a custom-fit similarity score algorithm and adjustable matching results notifications.

Severe congenital lactic acidosis and hypertrophic cardiomyopathy caused by an intronic variant in NDUFB7.

Mutations in structural subunits and assembly factors of complex I of the oxidative phosphorylation system constitute the most common cause of mitochondrial respiratory chain defects. Such mutations can present a wide range of clinical manifestations, varying from mild deficiencies to severe, lethal disorders. We describe a patient presenting intrauterine growth restriction and anemia, which displayed postpartum hypertrophic cardiomyopathy, lactic acidosis, encephalopathy, and a severe complex I defect with fatal outcome. Whole genome sequencing revealed an intronic biallelic mutation in the NDUFB7 gene (c.113-10C>G) and splicing pattern alterations in NDUFB7 messenger RNA were confirmed by RNA Sequencing. The detected variant resulted in a significant reduction of the NDUFB7 protein and reduced complex I activity. Complementation studies with expression of wild-type NDUFB7 in patient fibroblasts normalized complex I function. Here we report a case with a primary complex I defect due to a homozygous mutation in an intron region of the NDUFB7 gene.

Clinical Presentation, Genetic Etiology, and Coenzyme Q10 Levels in 55 Children with Combined Enzyme Deficiencies of the Mitochondrial Respiratory Chain.

OBJECTIVES: To evaluate the clinical symptoms and biochemical findings and establish the genetic etiology in a cohort of pediatric patients with combined deficiencies of the mitochondrial respiratory chain complexes. STUDY DESIGN: Clinical and biochemical data were collected from 55 children. All patients were subjected to sequence analysis of the entire mitochondrial genome, except when the causative mutations had been identified based on the clinical picture. Whole exome sequencing/whole genome sequencing (WES/WGS) was performed in 32 patients. RESULTS: Onset of disease was generally early in life (median age, 6 weeks). The most common symptoms were muscle weakness, hypotonia, and developmental delay/intellectual disability. Nonneurologic symptoms were frequent. Disease causing mutations were found in 20 different nuclear genes, and 7 patients had mutations in mitochondrial DNA. Causative variants were found in 18 of the 32 patients subjected to WES/WGS. Interestingly, many patients had low levels of coenzyme Q10 in muscle, irrespective of genetic cause. CONCLUSIONS: Children with combined enzyme defects display a diversity of clinical symptoms with varying age of presentation. We established the genetic diagnosis in 35 of the 55 patients (64%). The high diagnostic yield was achieved by the introduction of massive parallel sequencing, which also revealed novel genes and enabled elucidation of new disease mechanisms.

High diagnostic yield in skeletal ciliopathies using massively parallel genome sequencing, structural variant screening and RNA analyses.

Skeletal ciliopathies are a heterogenous group of disorders with overlapping clinical and radiographic features including bone dysplasia and internal abnormalities. To date, pathogenic variants in at least 30 genes, coding for different structural cilia proteins, are reported to cause skeletal ciliopathies. Here, we summarize genetic and phenotypic features of 34 affected individuals from 29 families with skeletal ciliopathies. Molecular diagnostic testing was performed using massively parallel sequencing (MPS) in combination with copy number variant (CNV) analyses and in silico filtering for variants in known skeletal ciliopathy genes. We identified biallelic disease-causing variants in seven genes: DYNC2H1, KIAA0753, WDR19, C2CD3, TTC21B, EVC, and EVC2. Four variants located in non-canonical splice sites of DYNC2H1, EVC, and KIAA0753 led to aberrant splicing that was shown by sequencing of cDNA. Furthermore, CNV analyses showed an intragenic deletion of DYNC2H1 in one individual and a 6.7 Mb de novo deletion on chromosome 1q24q25 in another. In five unsolved cases, MPS was performed in family setting. In one proband we identified a de novo variant in PRKACA and in another we found a homozygous intragenic deletion of IFT74, removing the first coding exon and leading to expression of a shorter message predicted to result in loss of 40 amino acids at the N-terminus. These findings establish IFT74 as a new skeletal ciliopathy gene. In conclusion, combined single nucleotide variant, CNV and cDNA analyses lead to a high yield of genetic diagnoses (90%) in a cohort of patients with skeletal ciliopathies.

Loqusdb: added value of an observations database of local genomic variation.

BACKGROUND: Exome and genome sequencing is becoming the method of choice for rare disease diagnostics. One of the key challenges remaining is distinguishing the disease causing variants from the benign background variation. After analysis and annotation of the sequencing data there are typically thousands of candidate variants requiring further investigation. One of the most effective and least biased ways to reduce this number is to assess the rarity of a variant in any population. Currently, there are a number of reliable sources of information for major population frequencies when considering single nucleotide variants (SNVs) and small insertion and deletions (INDELs), with gnomAD as the most prominent public resource available. However, local variation or frequencies in sub-populations may be underrepresented in these public resources. In contrast, for structural variation (SV), the background frequency in the general population is more or less unknown mostly due to challenges in calling SVs in a consistent way. Keeping track of local variation is one way to overcome these problems and significantly reduce the number of potential disease causing variants retained for manual inspection, both for SNVs and SVs. RESULTS: Here, we present loqusdb, a tool to solve the challenge of keeping track of any type of variant observations from genome sequencing data. Loqusdb was designed to handle a large flow of samples and unlike other solutions, samples can be added continuously to the database without rebuilding it, facilitating improvements and additions. We assessed the added value of a local observations database using 98 samples annotated with information from a background of 888 unrelated individuals. CONCLUSIONS: We show both how powerful SV analysis can be when filtering for population frequencies and how the number of apparently rare SNVs/INDELs can be reduced by adding local population information even after annotating the data with other large frequency databases, such as gnomAD. In conclusion, we show that a local frequency database is an attractive, and a necessary addition to the publicly available databases that facilitate the analysis of exome and genome data in a clinical setting.

Epilepsy syndromes, etiologies, and the use of next-generation sequencing in epilepsy presenting in the first 2 years of life: A population-based study.

OBJECTIVE: Population-based data on epilepsy syndromes and etiologies in early onset epilepsy are scarce. The use of next-generation sequencing (NGS) has hitherto not been reported in this context. The aim of this study is to describe children with epilepsy onset before 2 years of age, and to explore to what degree whole exome and whole genome sequencing (WES/WGS) can help reveal a molecular genetic diagnosis. METHODS: Children presenting with a first unprovoked epileptic seizure before age 2 years and registered in the Stockholm Incidence Registry of Epilepsy (SIRE) between September 1, 2001 and December 31, 2006, were retrieved and their medical records up to age 7 years reviewed. Children who met the epilepsy criteria were included in the study cohort. WES/WGS was offered in cases of suspected genetic etiology regardless of whether a structural or metabolic diagnosis had been established. RESULTS: One hundred sixteen children were included, of which 88 had seizure onset during the first year of life and 28 during the second, corresponding to incidences of 139 and 42/100 000 person-years, respectively. An epilepsy syndrome could be diagnosed in 54% of cases, corresponding to a birth prevalence of 1/1100. Structural etiology was revealed in 34% of cases, a genetic cause in 20%, and altogether etiology was known in 65% of children. The highest diagnostic yield was seen in magnetic resonance imaging (MRI) with 65% revealing an etiology. WES/WGS was performed in 26/116 cases (22%), with a diagnostic yield of 58%. SIGNIFICANCE: Epilepsy syndromes can be diagnosed and etiologies revealed in a majority of early onset cases. NGS can identify a molecular diagnosis in a substantial number of children, and should be included in the work-up, especially in cases of epileptic encephalopathy, cerebral malformation, or metabolic disease without molecular diagnosis. A genetic diagnosis is essential to genetic counselling, prenatal diagnostics, and precision therapy.

Absence of the Autophagy Adaptor SQSTM1/p62 Causes Childhood-Onset Neurodegeneration with Ataxia, Dystonia, and Gaze Palsy.

SQSTM1 (sequestosome 1; also known as p62) encodes a multidomain scaffolding protein involved in various key cellular processes, including the removal of damaged mitochondria by its function as a selective autophagy receptor. Heterozygous variants in SQSTM1 have been associated with Paget disease of the bone and might contribute to neurodegeneration in amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). Using exome sequencing, we identified three different biallelic loss-of-function variants in SQSTM1 in nine affected individuals from four families with a childhood- or adolescence-onset neurodegenerative disorder characterized by gait abnormalities, ataxia, dysarthria, dystonia, vertical gaze palsy, and cognitive decline. We confirmed absence of the SQSTM1/p62 protein in affected individuals' fibroblasts and found evidence of a defect in the early response to mitochondrial depolarization and autophagosome formation. Our findings expand the SQSTM1-associated phenotypic spectrum and lend further support to the concept of disturbed selective autophagy pathways in neurodegenerative diseases.

Intra-mitochondrial Methylation Deficiency Due to Mutations in SLC25A26.

S-adenosylmethionine (SAM) is the predominant methyl group donor and has a large spectrum of target substrates. As such, it is essential for nearly all biological methylation reactions. SAM is synthesized by methionine adenosyltransferase from methionine and ATP in the cytoplasm and subsequently distributed throughout the different cellular compartments, including mitochondria, where methylation is mostly required for nucleic-acid modifications and respiratory-chain function. We report a syndrome in three families affected by reduced intra-mitochondrial methylation caused by recessive mutations in the gene encoding the only known mitochondrial SAM transporter, SLC25A26. Clinical findings ranged from neonatal mortality resulting from respiratory insufficiency and hydrops to childhood acute episodes of cardiopulmonary failure and slowly progressive muscle weakness. We show that SLC25A26 mutations cause various mitochondrial defects, including those affecting RNA stability, protein modification, mitochondrial translation, and the biosynthesis of CoQ10 and lipoic acid.

Rescue of primary ubiquinone deficiency due to a novel COQ7 defect using 2,4-dihydroxybensoic acid.

BACKGROUND: Coenzyme Q is an essential mitochondrial electron carrier, redox cofactor and a potent antioxidant in the majority of cellular membranes. Coenzyme Q deficiency has been associated with a range of metabolic diseases, as well as with some drug treatments and ageing. METHODS: We used whole exome sequencing (WES) to investigate patients with inherited metabolic diseases and applied a novel ultra-pressure liquid chromatography-mass spectrometry approach to measure coenzyme Q in patient samples. RESULTS: We identified a homozygous missense mutation in the COQ7 gene in a patient with complex mitochondrial deficiency, resulting in severely reduced coenzyme Q levels We demonstrate that the coenzyme Q analogue 2,4-dihydroxybensoic acid (2,4DHB) was able to specifically bypass the COQ7 deficiency, increase cellular coenzyme Q levels and rescue the biochemical defect in patient fibroblasts. CONCLUSION: We report the first patient with primary coenzyme Q deficiency due to a homozygous COQ7 mutation and a potentially beneficial treatment using 2,4DHB.

Rapid pulsed whole genome sequencing for comprehensive acute diagnostics of inborn errors of metabolism.

BACKGROUND: Massively parallel DNA sequencing (MPS) has the potential to revolutionize diagnostics, in particular for monogenic disorders. Inborn errors of metabolism (IEM) constitute a large group of monogenic disorders with highly variable clinical presentation, often with acute, nonspecific initial symptoms. In many cases irreversible damage can be reduced by initiation of specific treatment, provided that a correct molecular diagnosis can be rapidly obtained. MPS thus has the potential to significantly improve both diagnostics and outcome for affected patients in this highly specialized area of medicine. RESULTS: We have developed a conceptually novel approach for acute MPS, by analysing pulsed whole genome sequence data in real time, using automated analysis combined with data reduction and parallelization. We applied this novel methodology to an in-house developed customized work flow enabling clinical-grade analysis of all IEM with a known genetic basis, represented by a database containing 474 disease genes which is continuously updated. As proof-of-concept, two patients were retrospectively analysed in whom diagnostics had previously been performed by conventional methods. The correct disease-causing mutations were identified and presented to the clinical team after 15 and 18 hours from start of sequencing, respectively. With this information available, correct treatment would have been possible significantly sooner, likely improving outcome. CONCLUSIONS: We have adapted MPS to fit into the dynamic, multidisciplinary work-flow of acute metabolic medicine. As the extent of irreversible damage in patients with IEM often correlates with timing and accuracy of management in early, critical disease stages, our novel methodology is predicted to improve patient outcome. All procedures have been designed such that they can be implemented in any technical setting and to any genetic disease area. The strategy conforms to international guidelines for clinical MPS, as only validated disease genes are investigated and as clinical specialists take responsibility for translation of results. As follow-up in patients without any known IEM, filters can be lifted and the full genome investigated, after genetic counselling and informed consent.

Individualised human phenotype ontology gene panels improve clinical whole exome and genome sequencing analytical efficacy in a cohort of developmental and epileptic encephalopathies.

BACKGROUND: The majority of genetic epilepsies remain unsolved in terms of specific genotype. Phenotype-based genomic analyses have shown potential to strengthen genomic analysis in various ways, including improving analytical efficacy. METHODS: We have tested a standardised phenotyping method termed 'Phenomodels' for integrating deep-phenotyping information with our in-house developed clinical whole exome/genome sequencing analytical pipeline. Phenomodels includes a user-friendly epilepsy phenotyping template and an objective measure for selecting which template terms to include in individualised Human Phenotype Ontology (HPO) gene panels. In a pilot study of 38 previously solved cases of developmental and epileptic encephalopathies, we compared the sensitivity and specificity of the individualised HPO gene panels with the clinical epilepsy gene panel. RESULTS: The Phenomodels template showed high sensitivity for capturing relevant phenotypic information, where 37/38 individuals' HPO gene panels included the causative gene. The HPO gene panels also had far fewer variants to assess than the epilepsy gene panel. CONCLUSION: We have demonstrated a viable approach for incorporating standardised phenotype information into clinical genomic analyses, which may enable more efficient analysis.

Comparison of total and cytoplasmic mRNA reveals global regulation by nuclear retention and miRNAs.

BACKGROUND: The majority of published gene-expression studies have used RNA isolated from whole cells, overlooking the potential impact of including nuclear transcriptome in the analyses. In this study, mRNA fractions from the cytoplasm and from whole cells (total RNA) were prepared from three human cell lines and sequenced using massive parallel sequencing. RESULTS: For all three cell lines, of about 15000 detected genes approximately 400 to 1400 genes were detected in different amounts in the cytoplasmic and total RNA fractions. Transcripts detected at higher levels in the total RNA fraction had longer coding sequences and higher number of miRNA target sites. Transcripts detected at higher levels in the cytoplasmic fraction were shorter or contained shorter untranslated regions. Nuclear retention of transcripts and mRNA degradation via miRNA pathway might contribute to this differential detection of genes. The consequence of the differential detection was further investigated by comparison to proteomics data. Interestingly, the expression profiles of cytoplasmic and total RNA correlated equally well with protein abundance levels indicating regulation at a higher level. CONCLUSIONS: We conclude that expression levels derived from the total RNA fraction be regarded as an appropriate estimate of the amount of mRNAs present in a given cell population, independent of the coding sequence length or UTRs.

Classification of DNA sequences using Bloom filters.

MOTIVATION: New generation sequencing technologies producing increasingly complex datasets demand new efficient and specialized sequence analysis algorithms. Often, it is only the 'novel' sequences in a complex dataset that are of interest and the superfluous sequences need to be removed. RESULTS: A novel algorithm, fast and accurate classification of sequences (FACSs), is introduced that can accurately and rapidly classify sequences as belonging or not belonging to a reference sequence. FACS was first optimized and validated using a synthetic metagenome dataset. An experimental metagenome dataset was then used to show that FACS achieves comparable accuracy as BLAT and SSAHA2 but is at least 21 times faster in classifying sequences. AVAILABILITY: Source code for FACS, Bloom filters and MetaSim dataset used is available at http://facs.biotech.kth.se. The Bloom::Faster 1.6 Perl module can be downloaded from CPAN at http://search.cpan.org/ approximately palvaro/Bloom-Faster-1.6/ CONTACTS: henrik.stranneheim@biotech.kth.se; joakiml@biotech.kth.se SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Novel Mutation m.10372A>G in MT-ND3 Causing Sensorimotor Axonal Polyneuropathy.

OBJECTIVE: To investigate the pathogenicity of a novel MT-ND3 mutation identified in a patient with adult-onset sensorimotor axonal polyneuropathy and report the clinical, morphologic, and biochemical findings. METHODS: Clinical assessments and morphologic and biochemical investigations of skeletal muscle and cultured myoblasts from the patient were performed. Whole-genome sequencing (WGS) of DNA from skeletal muscle and Sanger sequencing of mitochondrial DNA (mtDNA) from both skeletal muscle and cultured myoblasts were performed. Heteroplasmic levels of mutated mtDNA in different tissues were quantified by last-cycle hot PCR. RESULTS: Muscle showed ragged red fibers, paracrystalline inclusions, a significant reduction in complex I (CI) respiratory chain (RC) activity, and decreased adenosine triphosphate (ATP) production for all substrates used by CI. Sanger sequencing of DNA from skeletal muscle detected a unique previously unreported heteroplasmic mutation in mtDNA encoded MT-ND3, coding for a subunit in CI. WGS confirmed the mtDNA mutation but did not detect any other mutation explaining the disease. Cultured myoblasts, however, did not carry the mutation, and RC activity measurements in myoblasts were normal. CONCLUSIONS: We report a case with adult-onset sensorimotor axonal polyneuropathy caused by a novel mtDNA mutation in MT-ND3. Loss of heteroplasmy in blood, cultured fibroblasts and myoblasts from the patient, and normal measurement of RC activity of the myoblasts support pathogenicity of the mutation. These findings highlight the importance of mitochondrial investigations in patients presenting with seemingly idiopathic polyneuropathy, especially if muscle also is affected.

Integration of whole genome sequencing into a healthcare setting: high diagnostic rates across multiple clinical entities in 3219 rare disease patients.

BACKGROUND: We report the findings from 4437 individuals (3219 patients and 1218 relatives) who have been analyzed by whole genome sequencing (WGS) at the Genomic Medicine Center Karolinska-Rare Diseases (GMCK-RD) since mid-2015. GMCK-RD represents a long-term collaborative initiative between Karolinska University Hospital and Science for Life Laboratory to establish advanced, genomics-based diagnostics in the Stockholm healthcare setting. METHODS: Our analysis covers detection and interpretation of SNVs, INDELs, uniparental disomy, CNVs, balanced structural variants, and short tandem repeat expansions. Visualization of results for clinical interpretation is carried out in Scout-a custom-developed decision support system. Results from both singleton (84%) and trio/family (16%) analyses are reported. Variant interpretation is done by 15 expert teams at the hospital involving staff from three clinics. For patients with complex phenotypes, data is shared between the teams. RESULTS: Overall, 40% of the patients received a molecular diagnosis ranging from 19 to 54% for specific disease groups. There was heterogeneity regarding causative genes (n = 754) with some of the most common ones being COL2A1 (n = 12; skeletal dysplasia), SCN1A (n = 8; epilepsy), and TNFRSF13B (n = 4; inborn errors of immunity). Some causative variants were recurrent, including previously known founder mutations, some novel mutations, and recurrent de novo mutations. Overall, GMCK-RD has resulted in a large number of patients receiving specific molecular diagnoses. Furthermore, negative cases have been included in research studies that have resulted in the discovery of 17 published, novel disease-causing genes. To facilitate the discovery of new disease genes, GMCK-RD has joined international data sharing initiatives, including ClinVar, UDNI, Beacon, and MatchMaker Exchange. CONCLUSIONS: Clinical WGS at GMCK-RD has provided molecular diagnoses to over 1200 individuals with a broad range of rare diseases. Consolidation and spread of this clinical-academic partnership will enable large-scale national collaboration.

SLC12A2 mutations cause NKCC1 deficiency with encephalopathy and impaired secretory epithelia.

OBJECTIVE: To describe the phenotype in 2 sisters with a rare constellation of neurologic symptoms and secretory impairments and to identify the etiology by the use of whole-genome sequencing (WGS). METHODS: After an extensive workup failed to reveal the cause of disease, in a girl with a previously not reported phenotype, WGS of the proband, her diseased older sister, an older healthy brother, and their parents was performed, and potentially pathogenic variants were analyzed. RESULTS: The proband and her older sister both presented with neonatal Staphylococcus aureus parotitis, apneas, disappearance of the Moro reflex, and hypotonia. The proband survived. Her brain MRI showed white matter and basal ganglia abnormalities, and CSF damage biomarkers were increased. At age 8 years, she exhibits a constellation of symptoms including severe neurodevelopmental disorder, hearing impairment, gastrointestinal problems, and a striking lack of tear fluid, saliva, and sweat. Her respiratory mucosa is dry with potentially life-threatening mucus plugging. Through WGS, 2 loss-of-function variants in SLC12A2 were identified that follow an autosomal recessive inheritance pattern. CONCLUSIONS: Taken together with a single previously reported case and the close resemblance to the phenotypes of corresponding mouse models, our study firmly establishes biallelic variants in SLC12A2 as causing human disease and adds data regarding the neurologic phenotype.

A comparison between protein profiles of B cell subpopulations and mantle cell lymphoma cells.

BACKGROUND: B-cell lymphomas are thought to reflect different stages of B-cell maturation. Based on cytogenetics and molecular markers, mantle cell lymphoma (MCL) is presumed to derive predominantly from naïve, pre-germinal centre (pre-GC) B lymphocytes. The aim of this study was to develop a method to investigate the similarity between MCL cells and different B-cell compartments on a protein expression level. METHODS: Subpopulations of B cells representing the germinal centre (GC), the pre-GC mantle zone and the post-GC marginal zone were isolated from tonsils using automated magnetic cell sorting (AutoMACS) of cells based on their expression of CD27 and IgD. Protein profiling of the B cell subsets, of cell lines representing different lymphomas and of primary MCL samples was performed using top-down proteomics profiling by surface-enhanced laser detection/ionization time-of-flight mass spectrometry (SELDI-TOF-MS). RESULTS: Quantitative MS data of significant protein peaks (p-value < 0.05) separating the three B-cell subpopulations were generated. Together, hierarchical clustering and principal component analysis (PCA) showed that the primary MCL samples clustered together with the pre- and post-GC subpopulations. Both primary MCL cells and MCL cell lines were clearly separated from the B cells representing the GC compartment. CONCLUSION: AutoMACS sorting generates sufficient purity to enable a comparison between protein profiles of B cell subpopulations and malignant B lymphocytes applying SELDI-TOF-MS. Further validation with an increased number of patient samples and identification of differentially expressed proteins would enable a search for possible treatment targets that are expressed during the early development of MCL.

SQSTM1/p62-Directed Metabolic Reprogramming Is Essential for Normal Neurodifferentiation.

Neurodegenerative disorders are an increasingly common and irreversible burden on society, often affecting the aging population, but their etiology and disease mechanisms are poorly understood. Studying monogenic neurodegenerative diseases with known genetic cause provides an opportunity to understand cellular mechanisms also affected in more complex disorders. We recently reported that loss-of-function mutations in the autophagy adaptor protein SQSTM1/p62 lead to a slowly progressive neurodegenerative disease presenting in childhood. To further elucidate the neuronal involvement, we studied the cellular consequences of loss of p62 in a neuroepithelial stem cell (NESC) model and differentiated neurons derived from reprogrammed p62 patient cells or by CRISPR/Cas9-directed gene editing in NESCs. Transcriptomic and proteomic analyses suggest that p62 is essential for neuronal differentiation by controlling the metabolic shift from aerobic glycolysis to oxidative phosphorylation required for neuronal maturation. This shift is blocked by the failure to sufficiently downregulate lactate dehydrogenase expression due to the loss of p62, possibly through impaired Hif-1α downregulation and increased sensitivity to oxidative stress. The findings imply an important role for p62 in neuronal energy metabolism and particularly in the regulation of the shift between glycolysis and oxidative phosphorylation required for normal neurodifferentiation.

From cytogenetics to cytogenomics: whole-genome sequencing as a first-line test comprehensively captures the diverse spectrum of disease-causing genetic variation underlying intellectual disability.

BACKGROUND: Since different types of genetic variants, from single nucleotide variants (SNVs) to large chromosomal rearrangements, underlie intellectual disability, we evaluated the use of whole-genome sequencing (WGS) rather than chromosomal microarray analysis (CMA) as a first-line genetic diagnostic test. METHODS: We analyzed three cohorts with short-read WGS: (i) a retrospective cohort with validated copy number variants (CNVs) (cohort 1, n = 68), (ii) individuals referred for monogenic multi-gene panels (cohort 2, n = 156), and (iii) 100 prospective, consecutive cases referred to our center for CMA (cohort 3). Bioinformatic tools developed include FindSV, SVDB, Rhocall, Rhoviz, and vcf2cytosure. RESULTS: First, we validated our structural variant (SV)-calling pipeline on cohort 1, consisting of three trisomies and 79 deletions and duplications with a median size of 850 kb (min 500 bp, max 155 Mb). All variants were detected. Second, we utilized the same pipeline in cohort 2 and analyzed with monogenic WGS panels, increasing the diagnostic yield to 8%. Next, cohort 3 was analyzed by both CMA and WGS. The WGS data was processed for large (> 10 kb) SVs genome-wide and for exonic SVs and SNVs in a panel of 887 genes linked to intellectual disability as well as genes matched to patient-specific Human Phenotype Ontology (HPO) phenotypes. This yielded a total of 25 pathogenic variants (SNVs or SVs), of which 12 were detected by CMA as well. We also applied short tandem repeat (STR) expansion detection and discovered one pathologic expansion in ATXN7. Finally, a case of Prader-Willi syndrome with uniparental disomy (UPD) was validated in the WGS data. Important positional information was obtained in all cohorts. Remarkably, 7% of the analyzed cases harbored complex structural variants, as exemplified by a ring chromosome and two duplications found to be an insertional translocation and part of a cryptic unbalanced translocation, respectively. CONCLUSION: The overall diagnostic rate of 27% was more than doubled compared to clinical microarray (12%). Using WGS, we detected a wide range of SVs with high accuracy. Since the WGS data also allowed for analysis of SNVs, UPD, and STRs, it represents a powerful comprehensive genetic test in a clinical diagnostic laboratory setting.

Respiratory chain complex III deficiency due to mutated BCS1L: a novel phenotype with encephalomyopathy, partially phenocopied in a Bcs1l mutant mouse model.

BACKGROUND: Mitochondrial diseases due to defective respiratory chain complex III (CIII) are relatively uncommon. The assembly of the eleven-subunit CIII is completed by the insertion of the Rieske iron-sulfur protein, a process for which BCS1L protein is indispensable. Mutations in the BCS1L gene constitute the most common diagnosed cause of CIII deficiency, and the phenotypic spectrum arising from mutations in this gene is wide. RESULTS: A case of CIII deficiency was investigated in depth to assess respiratory chain function and assembly, and brain, skeletal muscle and liver histology. Exome sequencing was performed to search for the causative mutation(s). The patient's platelets and muscle mitochondria showed respiration defects and defective assembly of CIII was detected in fibroblast mitochondria. The patient was compound heterozygous for two novel mutations in BCS1L, c.306A > T and c.399delA. In the cerebral cortex a specific pattern of astrogliosis and widespread loss of microglia was observed. Further analysis showed loss of Kupffer cells in the liver. These changes were not found in infants suffering from GRACILE syndrome, the most severe BCS1L-related disorder causing early postnatal mortality, but were partially corroborated in a knock-in mouse model of BCS1L deficiency. CONCLUSIONS: We describe two novel compound heterozygous mutations in BCS1L causing CIII deficiency. The pathogenicity of one of the mutations was unexpected and points to the importance of combining next generation sequencing with a biochemical approach when investigating these patients. We further show novel manifestations in brain, skeletal muscle and liver, including abnormality in specialized resident macrophages (microglia and Kupffer cells). These novel phenotypes forward our understanding of CIII deficiencies caused by BCS1L mutations.