ABSTRACT
Here we report a stop-mutation in the BOD1 (Biorientation Defective 1) gene, which co-segregates with intellectual disability in a large consanguineous family, where individuals that are homozygous for the mutation have no detectable BOD1 mRNA or protein. The BOD1 protein is required for proper chromosome segregation, regulating phosphorylation of PLK1 substrates by modulating Protein Phosphatase 2A (PP2A) activity during mitosis. We report that fibroblast cell lines derived from homozygous BOD1 mutation carriers show aberrant localisation of the cell cycle kinase PLK1 and its phosphatase PP2A at mitotic kinetochores. However, in contrast to the mitotic arrest observed in BOD1-siRNA treated HeLa cells, patient-derived cells progressed through mitosis with no apparent segregation defects but at an accelerated rate compared to controls. The relatively normal cell cycle progression observed in cultured cells is in line with the absence of gross structural brain abnormalities in the affected individuals. Moreover, we found that in normal adult brain tissues BOD1 expression is maintained at considerable levels, in contrast to PLK1 expression, and provide evidence for synaptic localization of Bod1 in murine neurons. These observations suggest that BOD1 plays a cell cycle-independent role in the nervous system. To address this possibility, we established two Drosophila models, where neuron-specific knockdown of BOD1 caused pronounced learning deficits and significant abnormalities in synapse morphology. Together our results reveal novel postmitotic functions of BOD1 as well as pathogenic mechanisms that strongly support a causative role of BOD1 deficiency in the aetiology of intellectual disability. Moreover, by demonstrating its requirement for cognitive function in humans and Drosophila we provide evidence for a conserved role of BOD1 in the development and maintenance of cognitive features.
Subject(s)
Cell Cycle Proteins/genetics , Cognition , Protein Phosphatase 2/genetics , Synapses/genetics , Animals , Chromosome Segregation/genetics , Drosophila/genetics , Drosophila/physiology , Fibroblasts/metabolism , Gene Expression Regulation , Gene Knockdown Techniques , HeLa Cells , Humans , Learning , Mice , Mitosis/genetics , Neurons/metabolism , Protein Serine-Threonine Kinases/genetics , Proto-Oncogene Proteins/genetics , Synapses/pathology , Polo-Like Kinase 1ABSTRACT
PMM2-CDG, formerly known as congenital disorder of glycosylation-Ia (CDG-Ia), is caused by mutations in the gene encoding phosphomannomutase 2 (PMM2). This disease is the most frequent form of inherited CDG-diseases affecting protein N-glycosylation in human. PMM2-CDG is a multisystemic disease with severe psychomotor and mental retardation. In order to study the pathophysiology of PMM2-CDG in a human cell culture model, we generated induced pluripotent stem cells (iPSCs) from fibroblasts of a PMM2-CDG-patient (PMM2-iPSCs). Expression of pluripotency factors andin vitrodifferentiation into cell types of the three germ layers was unaffected in the analyzed clone PMM2-iPSC-C3 compared with nondiseased human pluripotent stem cells (hPSCs), revealing no broader influence of the PMM2 mutation on pluripotency in cell culture. Analysis of gene expression by deep-sequencing did not show obvious differences in the transcriptome between PMM2-iPSC-C3 and nondiseased hPSCs. By multiplexed capillary gel electrophoresis coupled to laser induced fluorescence detection (xCGE-LIF) we could show that PMM2-iPSC-C3 exhibit the common hPSC N-glycosylation pattern with high-mannose-type N-glycans as the predominant species. However, phosphomannomutase activity of PMM2-iPSC-C3 was 27% compared with control hPSCs and lectin staining revealed an overall reduced protein glycosylation. In addition, quantitative assessment of N-glycosylation by xCGE-LIF showed an up to 40% reduction of high-mannose-type N-glycans in PMM2-iPSC-C3, which was in concordance to the observed reduction of the Glc3Man9GlcNAc2 lipid-linked oligosaccharide compared with control hPSCs. Thus we could model the PMM2-CDG disease phenotype of hypoglycosylation with patient derived iPSCsin vitro Knock-down ofPMM2by shRNA in PMM2-iPSC-C3 led to a residual activity of 5% and to a further reduction of the level of N-glycosylation. Taken together we have developed human stem cell-based cell culture models with stepwise reduced levels of N-glycosylation now enabling to study the role of N-glycosylation during early human development.
Subject(s)
Congenital Disorders of Glycosylation/pathology , Glycomics/methods , Induced Pluripotent Stem Cells/metabolism , Models, Biological , Phosphotransferases (Phosphomutases)/deficiency , Cells, Cultured , Congenital Disorders of Glycosylation/metabolism , Gene Expression Profiling/methods , Glycosylation , High-Throughput Nucleotide Sequencing/methods , Humans , Induced Pluripotent Stem Cells/pathology , Phosphotransferases (Phosphomutases)/metabolism , Polysaccharides/metabolismABSTRACT
Common diseases are often complex because they are genetically heterogeneous, with many different genetic defects giving rise to clinically indistinguishable phenotypes. This has been amply documented for early-onset cognitive impairment, or intellectual disability, one of the most complex disorders known and a very important health care problem worldwide. More than 90 different gene defects have been identified for X-chromosome-linked intellectual disability alone, but research into the more frequent autosomal forms of intellectual disability is still in its infancy. To expedite the molecular elucidation of autosomal-recessive intellectual disability, we have now performed homozygosity mapping, exon enrichment and next-generation sequencing in 136 consanguineous families with autosomal-recessive intellectual disability from Iran and elsewhere. This study, the largest published so far, has revealed additional mutations in 23 genes previously implicated in intellectual disability or related neurological disorders, as well as single, probably disease-causing variants in 50 novel candidate genes. Proteins encoded by several of these genes interact directly with products of known intellectual disability genes, and many are involved in fundamental cellular processes such as transcription and translation, cell-cycle control, energy metabolism and fatty-acid synthesis, which seem to be pivotal for normal brain development and function.
Subject(s)
Cognition Disorders/genetics , Genes, Recessive/genetics , High-Throughput Nucleotide Sequencing , Intellectual Disability/genetics , Brain/metabolism , Brain/physiology , Cell Cycle , Consanguinity , DNA Mutational Analysis , Exons/genetics , Gene Regulatory Networks , Genes, Essential/genetics , Homozygote , Humans , Metabolic Networks and Pathways , Mutation/genetics , Organ Specificity , Synapses/metabolismABSTRACT
The envisioned clinical and industrial use of human pluripotent stem cells and their derivatives has given major momentum to the establishment of suspension culture protocols that enable the mass production of cells. Understanding molecular changes accompanying the transfer from adherent to suspension culture is of utmost importance because this information can have a direct effect on the development of optimized culture conditions. In this study we assessed the gene expression of human embryonic stem cells and induced pluripotent stem cells grown in surface-adherent culture (two-dimensional) versus free-floating suspension culture spheroids (three-dimensional). We combined a quantitative proteomic approach based on stable isotope labeling by amino acids in cell culture with deep-sequencing-based transcriptomics. Cells in three-dimensional culture showed reduced expression of proteins forming structural components of cell-cell and cell-extracellular matrix junctions. However, fully unexpected, we found up-regulation of secreted inhibitors of the canonical Wnt signaling pathway and, concomitantly, a reduction in the level of active ß-catenin and in the expression of Wnt target genes. In Western blot analyses the cysteine protease calpain was shown to cleave E-cadherin and ß-catenin under three-dimensional culture conditions. Our data allowed the development of a model in which calpain cleavage of E-cadherin induces the disintegration of focal cell contacts and generates a 100-kDa E-cadherin fragment required for the formation of three-dimensional cell-cell contacts in spheroids. The parallel release of ß-catenin and its potential activation by calpain cleavage are counterbalanced by the overexpression of soluble Wnt pathway inhibitors. According to this model, calpain has a key function in the interplay between E-cadherin and ß-catenin-mediated intercellular adhesion and the canonical Wnt signaling pathway. Supporting this model, we show that pharmacological modulation of calpain activity prevents spheroid formation and causes disassembly of preexisting spheroids into single cells, thereby providing novel strategies for improving suspension culture conditions for human pluripotent stem cells in the future.
Subject(s)
Cadherins/metabolism , Calpain/metabolism , Cell Culture Techniques/methods , Pluripotent Stem Cells/metabolism , Wnt Signaling Pathway , beta Catenin/metabolism , Cadherins/genetics , Calpain/antagonists & inhibitors , Cell Adhesion/drug effects , Cells, Cultured , Cysteine Proteinase Inhibitors/pharmacology , Gene Expression Regulation , Glycoproteins/pharmacology , High-Throughput Nucleotide Sequencing/methods , Humans , Isotope Labeling , Oligopeptides/pharmacology , Proteomics , Sequence Analysis, RNA/methods , Wnt Signaling Pathway/drug effects , beta Catenin/geneticsABSTRACT
With a prevalence between 1 and 3%, hereditary forms of intellectual disability (ID) are among the most important problems in health care. Particularly, autosomal-recessive forms of the disorder have a very heterogeneous molecular basis, and genes with an increased number of disease-causing mutations are not common. Here, we report on three different mutations (two nonsense mutations, c.679C>T [p.Gln227(∗)] and c.1114C>T [p.Gln372(∗)], as well as one splicing mutation, g.6622224A>C [p.Ile179Argfs(∗)192]) that cause a loss of the tRNA-methyltransferase-encoding NSUN2 main transcript in homozygotes. We identified the mutations by sequencing exons and exon-intron boundaries within the genomic region where the linkage intervals of three independent consanguineous families of Iranian and Kurdish origin overlapped with the previously described MRT5 locus. In order to gain further evidence concerning the effect of a loss of NSUN2 on memory and learning, we constructed a Drosophila model by deleting the NSUN2 ortholog, CG6133, and investigated the mutants by using molecular and behavioral approaches. When the Drosophila melanogaster NSUN2 ortholog was deleted, severe short-term-memory (STM) deficits were observed; STM could be rescued by re-expression of the wild-type protein in the nervous system. The humans homozygous for NSUN2 mutations showed an overlapping phenotype consisting of moderate to severe ID and facial dysmorphism (which includes a long face, characteristic eyebrows, a long nose, and a small chin), suggesting that mutations in this gene might even induce a syndromic form of ID. Moreover, our observations from the Drosophila model point toward an evolutionarily conserved role of RNA methylation in normal cognitive development.
Subject(s)
Codon, Nonsense , Genes, Recessive , Intellectual Disability/genetics , Methyltransferases/genetics , Adolescent , Adult , Animals , Child , Cloning, Molecular , Consanguinity , Drosophila/genetics , Exons , Female , Genetic Linkage , Genotype , Homozygote , Humans , Intellectual Disability/physiopathology , Male , Methyltransferases/metabolism , Middle Aged , Pedigree , Phenotype , Young AdultABSTRACT
BACKGROUND: Segmental duplications (SDs) are not evenly distributed along chromosomes. The reasons for this biased susceptibility to SD insertion are poorly understood. Accumulation of SDs is associated with increased genomic instability, which can lead to structural variants and genomic disorders such as the Williams-Beuren syndrome. Despite these adverse effects, SDs have become fixed in the human genome. Focusing on chromosome 7, which is particularly rich in interstitial SDs, we have investigated the distribution of SDs in the context of evolution and the three dimensional organisation of the chromosome in order to gain insights into the mutual relationship of SDs and chromatin topology. RESULTS: Intrachromosomal SDs preferentially accumulate in those segments of chromosome 7 that are homologous to marmoset chromosome 2. Although this formerly compact segment has been re-distributed to three different sites during primate evolution, we can show by means of public data on long distance chromatin interactions that these three intervals, and consequently the paralogous SDs mapping to them, have retained their spatial proximity in the nucleus. Focusing on SD clusters implicated in the aetiology of the Williams-Beuren syndrome locus we demonstrate by cross-species comparison that these SDs have inserted at the borders of a topological domain and that they flank regions with distinct DNA conformation. CONCLUSIONS: Our study suggests a link of nuclear architecture and the propagation of SDs across chromosome 7, either by promoting regional SD insertion or by contributing to the establishment of higher order chromatin organisation themselves. The latter could compensate for the high risk of structural rearrangements and thus may have contributed to their evolutionary fixation in the human genome.
Subject(s)
Chromatin/genetics , Chromosomes, Human, Pair 7 , Segmental Duplications, Genomic , Acetylation , Chromatin/metabolism , Chromosomes, Human, Pair 2 , Epistasis, Genetic , Evolution, Molecular , Genetic Loci , Genomics , Histones/metabolism , Humans , Transcription, Genetic , Williams Syndrome/geneticsABSTRACT
OBJECTIVES: To better understand the molecular pathogenesis of T-cell large granular lymphocyte leukemia (T-LGL), we decided to search for those genetic alterations in T-LGL patients and MOTN-1 cell line (established from T-LGL patient) that have an impact on gene expression and as a result can influence cell biology. METHODS: Multicolor fluorescence in situ hybridization (mFISH) analysis of the MOTN-1 cell line was performed as well as paired-end next-generation sequencing (NGS; Illumina HiSeq2000) of this cell line and one T-LGL patient. In addition, chosen 6q region was characterized in three T-LGL patients using high-resolution comparative genomic hybridization (FT-CGH) and LM-PCR. Gene expression was studied by RNA sequencing (RNAseq; SOLID5500). RESULTS: Rearrangements were detected within 1p and 2q in MOTN-1 affecting expression of FGR, ZEB2, and CASP8, and within 6q in MOTN-1 and one T-LGL patient affecting MAP3K5 and IFNGR1. Nineteen genes, among them FOXN3, RIN3, AKT1, PPP2R5C, were overexpressed as a result of an amplification in 14q in one T-LGL patient. Two novel fusion transcripts were identified: CASP8-ERBB4 in MOTN-1 and SBF1-PKHD1L1 in T-LGL patient. CONCLUSIONS: This study showed that submicroscopic genomic rearrangements change gene expression in T-LGL. Several genes involved in rearrangements were previously linked to cancer and survival pattern that characterizes T-LGL cells.
Subject(s)
Gene Expression Regulation, Leukemic , Gene Rearrangement, T-Lymphocyte , Leukemia, Large Granular Lymphocytic/genetics , Neoplasm Proteins/genetics , Cell Line, Tumor , Comparative Genomic Hybridization , Humans , In Situ Hybridization, Fluorescence , Leukemia, Large Granular Lymphocytic/pathologyABSTRACT
Mental retardation (MR) has a worldwide prevalence of around 2% and is a frequent cause of severe disability. Significant excess of MR in the progeny of consanguineous matings as well as functional considerations suggest that autosomal recessive forms of MR (ARMR) must be relatively common. To shed more light on the causes of autosomal recessive MR (ARMR), we have set out in 2003 to perform systematic clinical studies and autozygosity mapping in large consanguineous Iranian families with non-syndromic ARMR (NS-ARMR). As previously reported (Najmabadi et al. in Hum Genet 121:43-48, 2007), this led us to the identification of 12 novel ARMR loci, 8 of which had a significant LOD score (OMIM: MRT5-12). In the meantime, we and others have found causative gene defects in two of these intervals. Moreover, as reported here, tripling the size of our cohort has enabled us to identify 27 additional unrelated families with NS-ARMR and single-linkage intervals; 14 of these define novel loci for non-syndromic ARMR. Altogether, 13 out of 39 single linkage intervals observed in our cohort were found to cluster at 6 different loci on chromosomes, i.e., 1p34, 4q27, 5p15, 9q34, 11p11-q13 and 19q13, respectively. Five of these clusters consist of two significantly overlapping linkage intervals, and on chr 1p34, three single linkage intervals coincide, including the previously described MRT12 locus. The probability for this distribution to be due to chance is only 1.14 × 10(-5), as shown by Monte Carlo simulation. Thus, in contrast to our previous conclusions, these novel data indicate that common molecular causes of NS-ARMR do exist, and in the Iranian population, the most frequent ones may well account for several percent of the patients. These findings will be instrumental in the identification of the underlying genes.
Subject(s)
Intellectual Disability/genetics , Mutation , Chromosome Disorders , Family , Genes, Recessive , Iran , Monte Carlo MethodABSTRACT
OBJECTIVE: We are often asked to evaluate the postoperative elbow using MRI. In this article, we present both the normal appearance of the reconstructed ulnar collateral ligament (UCL) and a spectrum of UCL graft abnormalities at MR arthrography. CONCLUSION: The UCL is the primary medial stabilizer of the elbow during flexion and is susceptible to injury, particularly in overhead throwing athletes. In the setting of UCL injury, ligament reconstruction is the principal surgical procedure available to restore stability and function and to relieve pain with activity. We have found MR arthrography useful in assessing the reconstructed UCL.
Subject(s)
Collateral Ligaments/surgery , Elbow Joint/surgery , Magnetic Resonance Imaging/methods , Plastic Surgery Procedures , Adult , Arthroscopy , Collateral Ligaments/injuries , Contrast Media , Female , Gadolinium DTPA , Humans , Male , Treatment Outcome , Elbow InjuriesABSTRACT
Alzheimer's disease (AD) is a devastating neurodegenerative disorder characterized by early intraneuronal amyloid-ß (Aß) accumulation, extracellular deposition of Aß peptides, and intracellular hyperphosphorylated tau aggregates. These lesions cause dendritic and synaptic alterations and induce an inflammatory response in the diseased brain. Although the neuropathological characteristics of AD have been known for decades, the molecular mechanisms causing the disease are still under investigation. Studying gene expression changes in postmortem AD brain tissue can yield new insights into the molecular disease mechanisms. To that end, one can employ transgenic AD mouse models and the next-generation sequencing technology. In this study, a whole-brain transcriptome analysis was carried out using the well-characterized APP/PS1KI mouse model for AD. These mice display a robust phenotype reflected by working memory deficits at 6 months of age, a significant neuron loss in a variety of brain areas including the CA1 region of the hippocampus and a severe amyloid pathology. Based on deep sequencing, differentially expressed genes (DEGs) between 6-month-old WT or PS1KI and APP/PS1KI were identified and verified by qRT-PCR. Compared to WT mice, 250 DEGs were found in APP/PS1KI mice, while 186 DEGs could be found compared to PS1KI control mice. Most of the DEGs were upregulated in APP/PS1KI mice and belong to either inflammation-associated pathways or lysosomal activation, which is likely due to the robust intraneuronal accumulation of Aß in this mouse model. Our comprehensive brain transcriptome study further highlights APP/PS1KI mice as a valuable model for AD, covering molecular inflammatory and immune responses.
Subject(s)
Alzheimer Disease/genetics , Brain/metabolism , Transcriptome , Alzheimer Disease/pathology , Amyloid beta-Protein Precursor/genetics , Amyloid beta-Protein Precursor/metabolism , Animals , Brain/pathology , Disease Models, Animal , Gene Expression Profiling , Male , Mice , Mice, Transgenic , Nerve Degeneration/pathology , Presenilin-1/genetics , Presenilin-1/metabolismABSTRACT
Ionizing radiation is known to induce genomic lesions, such as DNA double strand breaks, whose repair can lead to mutations that can modulate cellular and organismal fate. Soon after radiation exposure, cells induce transcriptional changes and alterations of cell cycle programs to respond to the received DNA damage. Radiation-induced mutations occur through misrepair in a stochastic manner and increase the risk of developing cancers years after the incident, especially after high dose radiation exposures. Here, the authors analyzed the transcriptomic response of primary human gingival fibroblasts exposed to increasing doses of acute high dose-rate x rays. In the dataset obtained after 0.5 and 5 Gy x-ray exposures and two different repair intervals (0.5 h and 16 h), the authors discovered several radiation-induced fusion transcripts in conjunction with dose-dependent gene expression changes involving a total of 3,383 genes. Principal component analysis of repeated experiments revealed that the duration of the post-exposure repair intervals had a stronger impact than irradiation dose. Subsequent overrepresentation analyses showed a number of KEGG gene sets and WikiPathways, including pathways known to relate to radioresistance in fibroblasts (Wnt, integrin signaling). Moreover, a significant radiation-induced modulation of microRNA targets was detected. The data sets on IR-induced transcriptomic alterations in primary gingival fibroblasts will facilitate genomic comparisons in various genotoxic exposure scenarios.
Subject(s)
Fibroblasts/physiology , Fibroblasts/radiation effects , Gingiva/physiology , Gingiva/radiation effects , Transcriptome/physiology , Transcriptome/radiation effects , Cells, Cultured , Dose-Response Relationship, Drug , Fibroblasts/cytology , Gene Expression Regulation/physiology , Gene Expression Regulation/radiation effects , Humans , Radiation Dosage , Radiation Tolerance/physiology , Radiation Tolerance/radiation effects , X-RaysABSTRACT
Normal and painful stimuli are detected by specialized subgroups of peripheral sensory neurons. The understanding of the functional differences of each neuronal subgroup would be strongly enhanced by knowledge of the respective subgroup transcriptome. The separation of the subgroup of interest, however, has proven challenging as they can hardly be enriched. Instead of enriching, we now rapidly eliminated the subgroup of neurons expressing the heat-gated cation channel TRPV1 from dissociated rat sensory ganglia. Elimination was accomplished by brief treatment with TRPV1 agonists followed by the removal of compromised TRPV1(+) neurons using density centrifugation. By differential microarray and sequencing (RNA-Seq) based expression profiling we compared the transcriptome of all cells within sensory ganglia versus the same cells lacking TRPV1 expressing neurons, which revealed 240 differentially expressed genes (adj. p<0.05, fold-change>1.5). Corroborating the specificity of the approach, many of these genes have been reported to be involved in noxious heat or pain sensitization. Beyond the expected enrichment of ion channels, we found the TRPV1 transcriptome to be enriched for GPCRs and other signaling proteins involved in adenosine, calcium, and phosphatidylinositol signaling. Quantitative population analysis using a recent High Content Screening (HCS) microscopy approach identified substantial heterogeneity of expressed target proteins even within TRPV1-positive neurons. Signaling components defined distinct further subgroups within the population of TRPV1-positive neurons. Analysis of one such signaling system showed that the pain sensitizing prostaglandin PGD2 activates DP1 receptors expressed predominantly on TRPV1(+) neurons. In contrast, we found the PGD2 producing prostaglandin D synthase to be expressed exclusively in myelinated large-diameter neurons lacking TRPV1, which suggests a novel paracrine neuron-neuron communication. Thus, subgroup analysis based on the elimination rather than enrichment of the subgroup of interest revealed proteins that define subclasses of TRPV1-positive neurons and suggests a novel paracrine circuit.
Subject(s)
Ganglia, Spinal/metabolism , Nociceptors/metabolism , Paracrine Communication/genetics , Prostaglandin D2/metabolism , TRPV Cation Channels/metabolism , Animals , Base Sequence , Capsaicin/pharmacology , Cells, Cultured , Gene Expression Profiling , Male , Pain/genetics , Phosphatidylinositols/metabolism , Potassium Channels/metabolism , Prostaglandin D2/biosynthesis , Rats , Rats, Sprague-Dawley , Receptors, G-Protein-Coupled/metabolism , Sequence Analysis, DNA , TRPV Cation Channels/agonists , TRPV Cation Channels/geneticsABSTRACT
One of the central research questions on the etiology of Alzheimer's disease (AD) is the elucidation of the molecular signatures triggered by the amyloid cascade of pathological events. Next-generation sequencing allows the identification of genes involved in disease processes in an unbiased manner. We have combined this technique with the analysis of two AD mouse models: (1) The 5XFAD model develops early plaque formation, intraneuronal Aß aggregation, neuron loss, and behavioral deficits. (2) The Tg4-42 model expresses N-truncated Aß4-42 and develops neuron loss and behavioral deficits albeit without plaque formation. Our results show that learning and memory deficits in the Morris water maze and fear conditioning tasks in Tg4-42 mice at 12 months of age are similar to the deficits in 5XFAD animals. This suggested that comparative gene expression analysis between the models would allow the dissection of plaque-related and -unrelated disease relevant factors. Using deep sequencing differentially expressed genes (DEGs) were identified and subsequently verified by quantitative PCR. Nineteen DEGs were identified in pre-symptomatic young 5XFAD mice, and none in young Tg4-42 mice. In the aged cohort, 131 DEGs were found in 5XFAD and 56 DEGs in Tg4-42 mice. Many of the DEGs specific to the 5XFAD model belong to neuroinflammatory processes typically associated with plaques. Interestingly, 36 DEGs were identified in both mouse models indicating common disease pathways associated with behavioral deficits and neuron loss.
ABSTRACT
We compare the results of three different assembler programs, Celera, Phrap and Mira2, for the same set of about a hundred thousand Sanger reads derived from an unknown bacterial genome. In difference to previous assembly comparisons we do not focus on speed of computation and numbers of assembled contigs but on how the different sequence assemblies agree by content. Threefold consistently assembled genome regions are identified in order to estimate a lower bound of erroneously identified single nucleotide polymorphisms (SNP) caused by nothing but the process of mathematical sequence assembly. We identified 509 sequence triplets common to all three de-novo assemblies spanning only 34% (3.3 Mb) of the bacterial genome with 175 of these regions (~1.5 Mb) including erroneous SNPs and insertion/deletions. Within these triplets this on average leads to one error per 7,155 base pairs. Replacing the assembler Mira2 by the most recent version Mira3, the letter number even drops to 5,923. Our results therefore suggest that a considerably high number of erroneous SNPs may be present in current sequence data and mathematicians should urgently take up research on numerical stability of sequence assembly algorithms. Furthermore, even the latest versions of currently used assemblers produce erroneous SNPs that depend on the order reads are used as input. Such errors will severely hamper molecular diagnostics as well as relating genome variation and disease. This issue needs to be addressed urgently as the field is moving fast into clinical applications.