The DNMT3A ADD domain is required for efficient de novo DNA methylation and maternal imprinting in mouse oocytes.

Establishment of a proper DNA methylation landscape in mammalian oocytes is important for maternal imprinting and embryonic development. De novo DNA methylation in oocytes is mediated by the DNA methyltransferase DNMT3A, which has an ATRX-DNMT3-DNMT3L (ADD) domain that interacts with histone H3 tail unmethylated at lysine-4 (H3K4me0). The domain normally blocks the methyltransferase domain via intramolecular interaction and binding to histone H3K4me0 releases the autoinhibition. However, H3K4me0 is widespread in chromatin and the role of the ADD-histone interaction has not been studied in vivo. We herein show that amino-acid substitutions in the ADD domain of mouse DNMT3A cause dwarfism. Oocytes derived from homozygous females show mosaic loss of CG methylation and almost complete loss of non-CG methylation. Embryos derived from such oocytes die in mid-to-late gestation, with stochastic and often all-or-none-type CG-methylation loss at imprinting control regions and misexpression of the linked genes. The stochastic loss is a two-step process, with loss occurring in cleavage-stage embryos and regaining occurring after implantation. These results highlight an important role for the ADD domain in efficient, and likely processive, de novo CG methylation and pose a model for stochastic inheritance of epigenetic perturbations in germ cells to the next generation.

The DNMT3A PWWP domain is essential for the normal DNA methylation landscape in mouse somatic cells and oocytes.

DNA methylation at CG sites is important for gene regulation and embryonic development. In mouse oocytes, de novo CG methylation requires preceding transcription-coupled histone mark H3K36me3 and is mediated by a DNA methyltransferase DNMT3A. DNMT3A has a PWWP domain, which recognizes H3K36me2/3, and heterozygous mutations in this domain, including D329A substitution, cause aberrant CG hypermethylation of regions marked by H3K27me3 in somatic cells, leading to a dwarfism phenotype. We herein demonstrate that D329A homozygous mice show greater CG hypermethylation and severer dwarfism. In oocytes, D329A substitution did not affect CG methylation of H3K36me2/3-marked regions, including maternally methylated imprinting control regions; rather, it caused aberrant hypermethylation in regions lacking H3K36me2/3, including H3K27me3-marked regions. Thus, the role of the PWWP domain in CG methylation seems similar in somatic cells and oocytes; however, there were cell-type-specific differences in affected regions. The major satellite repeat was also hypermethylated in mutant oocytes. Contrary to the CA hypomethylation in somatic cells, the mutation caused hypermethylation at CH sites, including CA sites. Surprisingly, oocytes expressing only the mutated protein could support embryonic and postnatal development. Our study reveals that the DNMT3A PWWP domain is important for suppressing aberrant CG hypermethylation in both somatic cells and oocytes but that D329A mutation has little impact on the developmental potential of oocytes.

CMIC: predicting DNA methylation inheritance of CpG islands with embedding vectors of variable-length k-mers.

BACKGROUND: Epigenetic modifications established in mammalian gametes are largely reprogrammed during early development, however, are partly inherited by the embryo to support its development. In this study, we examine CpG island (CGI) sequences to predict whether a mouse blastocyst CGI inherits oocyte-derived DNA methylation from the maternal genome. Recurrent neural networks (RNNs), including that based on gated recurrent units (GRUs), have recently been employed for variable-length inputs in classification and regression analyses. One advantage of this strategy is the ability of RNNs to automatically learn latent features embedded in inputs by learning their model parameters. However, the available CGI dataset applied for the prediction of oocyte-derived DNA methylation inheritance are not large enough to train the neural networks. RESULTS: We propose a GRU-based model called CMIC (CGI Methylation Inheritance Classifier) to augment CGI sequence by converting it into variable-length k-mers, where the length k is randomly selected from the range [Formula: see text] to [Formula: see text], N times, which were then used as neural network input. N was set to 1000 in the default setting. In addition, we proposed a new embedding vector generator for k-mers called splitDNA2vec. The randomness of this procedure was higher than the previous work, dna2vec. CONCLUSIONS: We found that CMIC can predict the inheritance of oocyte-derived DNA methylation at CGIs in the maternal genome of blastocysts with a high F-measure (0.93). We also show that the F-measure can be improved by increasing the parameter N, that is, the number of sequences of variable-length k-mers derived from a single CGI sequence. This implies the effectiveness of augmenting input data by converting a DNA sequence to N sequences of variable-length k-mers. This approach can be applied to different DNA sequence classification and regression analyses, particularly those involving a small amount of data.

Aging of spermatogonial stem cells by Jnk-mediated glycolysis activation.

Because spermatogonial stem cells (SSCs) are immortal by serial transplantation, SSC aging in intact testes is considered to be caused by a deteriorated microenvironment. Here, we report a cell-intrinsic mode of SSC aging by glycolysis activation. Using cultured SSCs, we found that aged SSCs proliferated more actively than young SSCs and showed enhanced glycolytic activity. Moreover, they remained euploid and exhibited stable androgenetic imprinting patterns with robust SSC activity despite having shortened telomeres. Aged SSCs showed increased Wnt7b expression, which was associated with decreased Polycomb complex 2 activity. Our results suggest that aberrant Wnt7b expression activated c-jun N-terminal kinase (JNK), which down-regulated mitochondria numbers by suppressing Ppargc1a Down-regulation of Ppargc1a probably decreased reactive oxygen species and enhanced glycolysis. Analyses of the Klotho-deficient aging mouse model and 2-y-old aged rats confirmed JNK hyperactivation and increased glycolysis. Therefore, not only microenvironment but also intrinsic activation of JNK-mediated glycolysis contributes to SSC aging.

Non-transmissible MV Vector with Segmented RNA Genome Establishes Different Types of iPSCs from Hematopoietic Cells.

Recent advances in gene therapy technologies have enabled the treatment of congenital disorders and cancers and facilitated the development of innovative methods, including induced pluripotent stem cell (iPSC) production and genome editing. We recently developed a novel non-transmissible and non-integrating measles virus (MV) vector capable of transferring multiple genes simultaneously into a wide range of cells through the CD46 and CD150 receptors. The MV vector expresses four genes for iPSC generation and the GFP gene for a period of time sufficient to establish iPSCs from human fibroblasts as well as peripheral blood T cells. The transgenes were expressed differentially depending on their gene order in the vector. Human hematopoietic stem/progenitor cells were directly and efficiently reprogrammed to naive-like cells that could proliferate and differentiate into primed iPSCs by the same method used to establish primed iPSCs from other cell types. The novel MV vector has several advantages for establishing iPSCs and potential future applications in gene therapy.

Role of UHRF1 in de novo DNA methylation in oocytes and maintenance methylation in preimplantation embryos.

The methylation of cytosine at CG sites in the mammalian genome is dynamically reprogrammed during gametogenesis and preimplantation development. It was previously shown that oocyte-derived DNMT1 (a maintenance methyltransferase) is essential for maintaining and propagating CG methylation at imprinting control regions in preimplantation embryos. In mammalian somatic cells, hemimethylated-CG-binding protein UHRF1 plays a critical role in maintaining CG methylation by recruiting DNMT1 to hemimethylated CG sites. However, the role of UHRF1 in oogenesis and preimplantation development is unknown. In the present study, we show that UHRF1 is mainly, but not exclusively, localized in the cytoplasm of oocytes and preimplantation embryos. However, smaller amounts of UHRF1 existed in the nucleus, consistent with the expected role in DNA methylation. We then generated oocyte-specific Uhrf1 knockout (KO) mice and found that, although oogenesis was itself unaffected, a large proportion of the embryos derived from the KO oocytes died before reaching the blastocyst stage (a maternal effect). Whole genome bisulfite sequencing revealed that blastocysts derived from KO oocytes have a greatly reduced level of CG methylation, suggesting that maternal UHRF1 is essential for maintaining CG methylation, particularly at the imprinting control regions, in preimplantation embryos. Surprisingly, UHRF1 was also found to contribute to de novo CG and non-CG methylation during oocyte growth: in Uhrf1 KO oocytes, transcriptionally-inactive regions gained less methylation, while actively transcribed regions, including the imprinting control regions, were unaffected or only slightly affected. We also found that de novo methylation was defective during the late stage of oocyte growth. To the best of our knowledge, this is the first study to demonstrate the role of UHRF1 in de novo DNA methylation in vivo. Our study reveals multiple functions of UHRF1 during the global epigenetic reprogramming of oocytes and early embryos.

An atomic model of Zfp57 recognition of CpG methylation within a specific DNA sequence.

Zinc finger transcription factor Zfp57 recognizes the methylated CpG within the TGCCGC element. We determined the structure of the DNA-binding domain of Zfp57, consisting of two adjacent zinc fingers, in complex with fully methylated DNA at 1.0 Å resolution. The first zinc finger contacts the 5' half (TGC), and the second recognizes the 3' half (CGC) of the recognition sequence. Zfp57 recognizes the two 5-methylcytosines (5mCs) asymmetrically: One involves hydrophobic interactions with Arg178, which also interacts with the neighboring 3' guanine and forms a 5mC-Arg-G interaction, while the other involves a layer of ordered water molecules. Two point mutations in patients with transient neonatal diabetes abolish DNA-binding activity. Zfp57 has reduced binding affinity for unmodified DNA and the oxidative products of 5mC.

Allele-Specific Methylome and Transcriptome Analysis Reveals Widespread Imprinting in the Human Placenta.

DNA methylation is globally reprogrammed after fertilization, and as a result, the parental genomes have similar DNA-methylation profiles after implantation except at the germline differentially methylated regions (gDMRs). We and others have previously shown that human blastocysts might contain thousands of transient maternally methylated gDMRs (transient mDMRs), whose maternal methylation is lost in embryonic tissues after implantation. In this study, we performed genome-wide allelic DNA methylation analyses of purified trophoblast cells from human placentas and, surprisingly, found that more than one-quarter of the transient-in-embryo mDMRs maintained their maternally biased DNA methylation. RNA-sequencing-based allelic expression analyses revealed that some of the placenta-specific mDMRs were associated with expression of imprinted genes (e.g., TIGAR, SLC4A7, PROSER2-AS1, and KLHDC10), and three imprinted gene clusters were identified. This approach also identified some X-linked gDMRs. Comparisons of the data with those from other mammals revealed that genomic imprinting in the placenta is highly variable. These findings highlight the incomplete erasure of germline DNA methylation in the human placenta; understanding this erasure is important for understanding normal placental development and the pathogenesis of developmental disorders with imprinting effects.

Roles of MIWI, MILI and PLD6 in small RNA regulation in mouse growing oocytes.

The mouse PIWI-interacting RNA (piRNA) pathway produces a class of 26-30-nucleotide (nt) small RNAs and is essential for spermatogenesis and retrotransposon repression. In oocytes, however, its regulation and function are poorly understood. In the present study, we investigated the consequences of loss of piRNA-pathway components in growing oocytes. When MILI (or PIWIL2), a PIWI family member, was depleted by gene knockout, almost all piRNAs disappeared. This severe loss of piRNA was accompanied by an increase in transcripts derived from specific retrotransposons, especially IAPs. MIWI (or PIWIL1) depletion had a smaller effect. In oocytes lacking PLD6 (or ZUCCHINI or MITOPLD), a mitochondrial nuclease/phospholipase involved in piRNA biogenesis in male germ cells, the piRNA level was decreased to 50% compared to wild-type, a phenotype much milder than that in males. Since PLD6 is essential for the creation of the 5΄ ends of primary piRNAs in males, the presence of mature piRNA in PLD6-depleted oocytes suggests the presence of compensating enzymes. Furthermore, we identified novel 21-23-nt small RNAs, termed spiRNAs, possessing a 10-nt complementarity with piRNAs, which were produced dependent on MILI and independent of DICER. Our study revealed the differences in the biogenesis and function of the piRNA pathway between sexes.

Software updates in the Illumina HiSeq platform affect whole-genome bisulfite sequencing.

BACKGROUND: Methylation of cytosine in genomic DNA is a well-characterized epigenetic modification involved in many cellular processes and diseases. Whole-genome bisulfite sequencing (WGBS), such as MethylC-seq and post-bisulfite adaptor tagging sequencing (PBAT-seq), uses the power of high-throughput DNA sequencers and provides genome-wide DNA methylation profiles at single-base resolution. However, the accuracy and consistency of WGBS outputs in relation to the operating conditions of high-throughput sequencers have not been explored. RESULTS: We have used the Illumina HiSeq platform for our PBAT-based WGBS, and found that different versions of HiSeq Control Software (HCS) and Real-Time Analysis (RTA) installed on the system provided different global CpG methylation levels (approximately 5% overall difference) for the same libraries. This problem was reproduced multiple times with different WGBS libraries and likely to be associated with the low sequence diversity of bisulfite-converted DNA. We found that HCS was the major determinant in the observed differences. To determine which version of HCS is most suitable for WGBS, we used substrates with predetermined CpG methylation levels, and found that HCS v2.0.5 is the best among the examined versions. HCS v2.0.12 showed the poorest performance and provided artificially lower CpG methylation levels when 5-methylcytosine is read as guanine (first read of PBAT-seq and second read of MethylC-seq). In addition, paired-end sequencing of low diversity libraries using HCS v2.2.38 or the latest HCS v2.2.58 was greatly affected by cluster densities. CONCLUSIONS: Software updates in the Illumina HiSeq platform can affect the outputs from low-diversity sequencing libraries such as WGBS libraries. More recent versions are not necessarily the better, and HCS v2.0.5 is currently the best for WGBS among the examined HCS versions. Thus, together with other experimental conditions, special care has to be taken on this point when CpG methylation levels are to be compared between different samples by WGBS.

Bifidobacteria can protect from enteropathogenic infection through production of acetate.

The human gut is colonized with a wide variety of microorganisms, including species, such as those belonging to the bacterial genus Bifidobacterium, that have beneficial effects on human physiology and pathology. Among the most distinctive benefits of bifidobacteria are modulation of host defence responses and protection against infectious diseases. Nevertheless, the molecular mechanisms underlying these effects have barely been elucidated. To investigate these mechanisms, we used mice associated with certain bifidobacterial strains and a simplified model of lethal infection with enterohaemorrhagic Escherichia coli O157:H7, together with an integrated 'omics' approach. Here we show that genes encoding an ATP-binding-cassette-type carbohydrate transporter present in certain bifidobacteria contribute to protecting mice against death induced by E. coli O157:H7. We found that this effect can be attributed, at least in part, to increased production of acetate and that translocation of the E. coli O157:H7 Shiga toxin from the gut lumen to the blood was inhibited. We propose that acetate produced by protective bifidobacteria improves intestinal defence mediated by epithelial cells and thereby protects the host against lethal infection.

Mouse oocyte methylomes at base resolution reveal genome-wide accumulation of non-CpG methylation and role of DNA methyltransferases.

DNA methylation is an epigenetic modification that plays a crucial role in normal mammalian development, retrotransposon silencing, and cellular reprogramming. Although methylation mainly occurs on the cytosine in a CG site, non-CG methylation is prevalent in pluripotent stem cells, brain, and oocytes. We previously identified non-CG methylation in several CG-rich regions in mouse germinal vesicle oocytes (GVOs), but the overall distribution of non-CG methylation and the enzymes responsible for this modification are unknown. Using amplification-free whole-genome bisulfite sequencing, which can be used with minute amounts of DNA, we constructed the base-resolution methylome maps of GVOs, non-growing oocytes (NGOs), and mutant GVOs lacking the DNA methyltransferase Dnmt1, Dnmt3a, Dnmt3b, or Dnmt3L. We found that nearly two-thirds of all methylcytosines occur in a non-CG context in GVOs. The distribution of non-CG methylation closely resembled that of CG methylation throughout the genome and showed clear enrichment in gene bodies. Compared to NGOs, GVOs were over four times more methylated at non-CG sites, indicating that non-CG methylation accumulates during oocyte growth. Lack of Dnmt3a or Dnmt3L resulted in a global reduction in both CG and non-CG methylation, showing that non-CG methylation depends on the Dnmt3a-Dnmt3L complex. Dnmt3b was dispensable. Of note, lack of Dnmt1 resulted in a slight decrease in CG methylation, suggesting that this maintenance enzyme plays a role in non-dividing oocytes. Dnmt1 may act on CG sites that remain hemimethylated in the de novo methylation process. Our results provide a basis for understanding the mechanisms and significance of non-CG methylation in mammalian oocytes.

Neuraminidase Amino Acid Sequences of Influenza A/H3N2 and B Viruses Isolated from Influenza Patients in the 2014/15 Japanese Influenza Season.

Background: Neuraminidase (NA) is a surface protein essential for influenza virus replication. NA inhibitors are commonly used for the treatment of influenza patients in Japan. Several mutations that reduce the effect of NA inhibitors have been reported. We sequenced the whole NA segment of isolated virus from influenza patients and investigated the relation between the NA amino acid sequence and the 50％ inhibitory concentration (IC_50) of four NA inhibitors. Materials and Methods: Forty A/H3N2 and 19 B influenza virus isolated from patients in the 2014/15 influenza season were analyzed. The IC_50 was determined by a neuraminidase inhibition assay using a fluorescent substrate. Viral RNA was amplified by RT-PCR and the genome was sequenced using a next generation sequencer. The deduced amino acid sequences were analyzed. Results: There was no AA change in the NA catalytic site of the A/H3N2 and B viruses isolated in the 2014-15 influenza season. There was no significant relation between the NA amino acids and the IC_50 of the four NA inhibitors for A/H3N2 or B viruses. Conclusion: The catalytic site of NA was highly conserved for these A/H3N2 and B viruses. No emergence of NA amino acid mutations related to the sensitivity of the four currently used NA inhibitors was observed.

DNA methylation and gene expression dynamics during spermatogonial stem cell differentiation in the early postnatal mouse testis.

BACKGROUND: In the male germline, neonatal prospermatogonia give rise to spermatogonia, which include stem cell population (undifferentiated spermatogonia) that supports continuous spermatogenesis in adults. Although the levels of DNA methyltransferases change dynamically in the neonatal and early postnatal male germ cells, detailed genome-wide DNA methylation profiles of these cells during the stem cell formation and differentiation have not been reported. RESULTS: To understand the regulation of spermatogonial stem cell formation and differentiation, we examined the DNA methylation and gene expression dynamics of male mouse germ cells at the critical stages: neonatal prospermatogonia, and early postntal (day 7) undifferentiated and differentiating spermatogonia. We found large partially methylated domains similar to those found in cancer cells and placenta in all these germ cells, and high levels of non-CG methylation and 5-hydroxymethylcytosines in neonatal prospermatogonia. Although the global CG methylation levels were stable in early postnatal male germ cells, and despite the reported scarcity of differential methylation in the adult spermatogonial stem cells, we identified many regions showing stage-specific differential methylation in and around genes important for stem cell function and spermatogenesis. These regions contained binding sites for specific transcription factors including the SOX family members. CONCLUSIONS: Our findings show a distinctive and dynamic regulation of DNA methylation during spermatogonial stem cell formation and differentiation in the neonatal and early postnatal testes. Furthermore, we revealed a unique accumulation and distribution of non-CG methylation and 5hmC marks in neonatal prospermatogonia. These findings contrast with the reported scarcity of differential methylation in adult spermatogonial stem cell differentiation and represent a unique phase of male germ cell development.

[Analysis of the Neuraminidase Amino Acid Sequences of Influenza A/H1N1pdm09, A/H3N2, and B Viruses Isolated from Influenza Patients in the 2013/14 Japanese Influenza Season].

BACKGROUND: Neuraminidase (NA) is an essential surface protein for influenza virus replication. NA inhibitors are commonly used for the treatment of influenza patients in Japan. Several mutations that reduce the effect of NA inhibitors have been reported. We sequenced the whole NA segment of isolated virus from influenza patients and investigated the relation between the NA amino acid sequence and the 50% inhibitory concentration (IC50) of four NA inhibitors. MATERIALS AND METHODS: A total of 20 viruses that showed high or low IC50 of NA inhibitors were selected from A/H1N1pdm09, A/H3N2, and B isolates from the viruses isolated from patients in the 2013-14 influenza season. Viral RNA was extracted and RT-PCR was done. The amplified genome was sequenced using a next generation sequencer", and the deduced amino acid sequences were analyzed. RESULTS: Two A/H1N1pdm09 viruses that showed very high IC50 for oseltamivir (150 nM and 130 nM) contained the H275Y mutation. Otherwise, no significant relation was found between the NA amino acids and the IC50 of the four NA inhibitors. There was no significant relation between the NA amino acids and the IC50 of the four NA inhibitors for A/H3N2 viruses. The B viruses that showed a high IC50 for oseltamivir and laninamivir shared some amino acids. The B viruses that showed a high IC50 of zanamivir and peramivir also shared some amino acids. They were different from the shared amino acids found for oseltamivir and laninamivir. CONCLUSION: The previously reported H275Y mutation that causes oseltamivir resistance was found in the two A/H1N1pdm09 viruses that showed a very high IC50 for oseltamivir. No additional NA amino acid sequences related to the IC50 of the four NA inhibitors was found. The meaning of the shared amino acids among B viruses that showed a high IC50 would be an interesting target for further investigation.

[Analysis of influenza A/H3N2 neuraminidase genes obtained from influenza patients in the 2011/12 and 2012/13 seasons in Japan].

BACKGROUND: Influenza virus has neuraminidase (NA), a surface protein with enzymatic activity that is essential for virus replication. Mutation may affect the effectiveness of NA inhibitors that are used for the treatment of influenza patients. In this study, we determined the NA gene sequences from the clinical isolates of influenza patients to examine the chronological genetic changes and the relation to drug susceptibility. METHODS: For 96 A/H3N2 virus isolates the 50% inhibitory concentration (IC50) (48 each from the 2011-12 and 12-13 influenza seasons) was measured. RT-PCR was done with extracted viral RNA, followed by nucleotide sequencing. RESULTS: One putative amino acid mutation, D151N, was found in an NA activity-related cite in five of ninety-six tested isolate. The mutation did not affect the IC50 value. The mutations identified at amino acid positions 387 and 400 were statistically correlated with an increased IC50 value, although the change was less than ten times, suggesting no significant difference in the clinical effectiveness. A small number .of isolates showed mutation in the T and/or B cell epitope region of NA. CONCLUSION: No mutation that affected the IC50 value or effectiveness of NAIs was detected. Antigenic mutations of NA, which influence the selection of epidemic strains, were not determined. Continuous observation will be necessary to further clarify the genetic features of NA.

Spatiotemporal DNA methylation dynamics shape megabase-scale methylome landscapes.

DNA methylation is an essential epigenetic mechanism that regulates cellular reprogramming and development. Studies using whole-genome bisulfite sequencing have revealed distinct DNA methylome landscapes in human and mouse cells and tissues. However, the factors responsible for the differences in megabase-scale methylome patterns between cell types remain poorly understood. By analyzing publicly available 258 human and 301 mouse whole-genome bisulfite sequencing datasets, we reveal that genomic regions rich in guanine and cytosine, when located near the nuclear center, are highly susceptible to both global DNA demethylation and methylation events during embryonic and germline reprogramming. Furthermore, we found that regions that generate partially methylated domains during global DNA methylation are more likely to resist global DNA demethylation, contain high levels of adenine and thymine, and are adjacent to the nuclear lamina. The spatial properties of genomic regions, influenced by their guanine-cytosine content, are likely to affect the accessibility of molecules involved in DNA (de)methylation. These properties shape megabase-scale DNA methylation patterns and change as cells differentiate, leading to the emergence of different megabase-scale methylome patterns across cell types.

Isolation and genomic analysis of a type IV pili-independent Thermus thermophilus phage, φMN1 from a Japanese hot spring.

A Thermus thermophilus lytic phage was isolated from a Japanese hot spring using a type IV pili-deficient strain as an indicator host, and designated as φMN1. Electron microscopic (EM) examination revealed that φMN1 had an icosahedral head and a contractile tail, suggesting that φMN1 belonged to Myoviridae. An EM analysis focused on φMN1 adsorption to the Thermus host cell showed that the receptor molecules for the phage were uniformly distributed on the outer surface of the cells. The circular double-stranded DNA of φMN1 was 76,659 base pairs in length, and the guanine and cytosine content was 61.8%. It was predicted to contain 99 open reading frames, and its putative distal tail fiber protein, which is essential for non-piliated host cell surface receptor recognition, was dissimilar in terms of sequence and length with its counterpart in the type IV pili-dependent φYS40. A phage proteomic tree revealed that φMN1 and φYS40 are in the same cluster, but many genes had low sequence similarities and some seemed to be derived from both mesophilic and thermophilic organisms. The gene organization suggested that φMN1 evolved from a non-Thermus phage through large-scale recombination events of the genes determining the host specificity, followed by gradual evolution by recombination of both the thermophilic and mesophilic DNAs assimilated by the host Thermus cells. This newly isolated phage will provide evolutionary insights into thermophilic phages.

Complete genome sequence of Bacillus cereus NC7401, which produces high levels of the emetic toxin cereulide.

We report the complete and annotated genome sequence of Bacillus cereus NC7401, a representative of the strain group that causes emetic-type food poisoning. The emetic toxin, cereulide, is produced by a nonribosomal protein synthesis (NRPS) system that is encoded by a gene cluster on a large resident plasmid, pNCcld.

Comparative genomics reveal the mechanism of the parallel evolution of O157 and non-O157 enterohemorrhagic Escherichia coli.

Among the various pathogenic Escherichia coli strains, enterohemorrhagic E. coli (EHEC) is the most devastating. Although serotype O157:H7 strains are the most prevalent, strains of different serotypes also possess similar pathogenic potential. Here, we present the results of a genomic comparison between EHECs of serotype O157, O26, O111, and O103, as well as 21 other, fully sequenced E. coli/Shigella strains. All EHECs have much larger genomes (5.5-5.9 Mb) than the other strains and contain surprisingly large numbers of prophages and integrative elements (IEs). The gene contents of the 4 EHECs do not follow the phylogenetic relationships of the strains, and they share virulence genes for Shiga toxins and many other factors. We found many lambdoid phages, IEs, and virulence plasmids that carry the same or similar virulence genes but have distinct evolutionary histories, indicating that independent acquisition of these mobile genetic elements has driven the evolution of each EHEC. Particularly interesting is the evolution of the type III secretion system (T3SS). We found that the T3SS of EHECs is composed of genes that were introduced by 3 different types of genetic elements: an IE referred to as the locus of enterocyte effacement, which encodes a central part of the T3SS; SpLE3-like IEs; and lambdoid phages carrying numerous T3SS effector genes and other T3SS-related genes. Our data demonstrate how E. coli strains of different phylogenies can independently evolve into EHECs, providing unique insights into the mechanisms underlying the parallel evolution of complex virulence systems in bacteria.