Complete genome assembly and methylome dissection of Methanococcus aeolicus PL15/H p.

Although restriction-modification systems are found in both Eubacterial and Archaeal kingdoms, comparatively less is known about patterns of DNA methylation and genome defense systems in archaea. Here we report the complete closed genome sequence and methylome analysis of Methanococcus aeolicus PL15/H p , a strain of the CO2-reducing methanogenic archaeon and a commercial source for MaeI, MaeII, and MaeIII restriction endonucleases. The M. aeolicus PL15/H p genome consists of a 1.68 megabase circular chromosome predicted to contain 1,615 protein coding genes and 38 tRNAs. A combination of methylome sequencing, homology-based genome annotation, and recombinant gene expression identified five restriction-modification systems encoded by this organism, including the methyltransferase and site-specific endonuclease of MaeIII. The MaeIII restriction endonuclease was recombinantly expressed, purified and shown to have site-specific DNA cleavage activity in vitro.

High-throughput sequencing of EcoWI restriction fragments maps the genome-wide landscape of phosphorothioate modification at base resolution.

Phosphorothioation (PT), in which a non-bridging oxygen is replaced by a sulfur, is one of the rare modifications discovered in bacteria and archaea that occurs on the sugar-phosphate backbone as opposed to the nucleobase moiety of DNA. While PT modification is widespread in the prokaryotic kingdom, how PT modifications are distributed in the genomes and their exact roles in the cell remain to be defined. In this study, we developed a simple and convenient technique called EcoWI-seq based on a modification-dependent restriction endonuclease to identify genomic positions of PT modifications. EcoWI-seq shows similar performance than other PT modification detection techniques and additionally, is easily scalable while requiring little starting material. As a proof of principle, we applied EcoWI-seq to map the PT modifications at base resolution in the genomes of both the Salmonella enterica cerro 87 and E. coli expressing the dnd+ gene cluster. Specifically, we address whether the partial establishment of modified PT positions is a stochastic or deterministic process. EcoWI-seq reveals a systematic usage of the same subset of target sites in clones for which the PT modification has been independently established.

Reassembling a cannon in the DNA defense arsenal: Genetics of StySA, a BREX phage exclusion system in Salmonella lab strains.

Understanding mechanisms that shape horizontal exchange in prokaryotes is a key problem in biology. A major limit on DNA entry is imposed by restriction-modification (RM) processes that depend on the pattern of DNA modification at host-specified sites. In classical RM, endonucleolytic DNA cleavage follows detection of unprotected sites on entering DNA. Recent investigation has uncovered BREX (BacteRiophage EXclusion) systems. These RM-like activities employ host protection by DNA modification, but immediate replication arrest occurs without evident of nuclease action on unmodified phage DNA. Here we show that the historical stySA RM locus of Salmonella enterica sv Typhimurium is a variant BREX system. A laboratory strain disabled for both the restriction and methylation activity of StySA nevertheless has wild type sequence in pglX, the modification gene homolog. Instead, flanking genes pglZ and brxC each carry multiple mutations (µ) in their C-terminal domains. We further investigate this system in situ, replacing the mutated pglZµ and brxCµ genes with the WT counterpart. PglZ-WT supports methylation in the presence of either BrxCµ or BrxC-WT but not in the presence of a deletion/insertion allele, ΔbrxC::cat. Restriction requires both BrxC-WT and PglZ-WT, implicating the BrxC C-terminus specifically in restriction activity. These results suggests that while BrxC, PglZ and PglX are principal components of the BREX modification activity, BrxL is required for restriction only. Furthermore, we show that a partial disruption of brxL disrupts transcription globally.

Comparative Genome Analysis of the Genus Thiothrix Involving Three Novel Species, Thiothrix subterranea sp. nov. Ku-5, Thiothrix litoralis sp. nov. AS and "Candidatus Thiothrix anitrata" sp. nov. A52, Revealed the Conservation of the Pathways of Dissimilatory Sulfur Metabolism and Variations in the Genetic Inventory for Nitrogen Metabolism and Autotrophic Carbon Fixation.

Two strains of filamentous, colorless sulfur bacteria were isolated from bacterial fouling in the outflow of hydrogen sulfide-containing waters from a coal mine (Thiothrix sp. Ku-5) and on the seashore of the White Sea (Thiothrix sp. AS). Metagenome-assembled genome (MAG) A52 was obtained from a sulfidic spring in the Volgograd region, Russia. Phylogenetic analysis based on the 16S rRNA gene sequences showed that all genomes represented the genus Thiothrix. Based on their average nucleotide identity and digital DNA-DNA hybridization data these new isolates and the MAG represent three species within the genus Thiothrix with the proposed names Thiothrix subterranea sp. nov. Ku-5T, Thiothrix litoralis sp. nov. AST, and "Candidatus Thiothrix anitrata" sp. nov. A52. The complete genome sequences of Thiothrix fructosivorans QT and Thiothrix unzii A1T were determined. Complete genomes of seven Thiothrix isolates, as well as two MAGs, were used for pangenome analysis. The Thiothrix core genome consisted of 1,355 genes, including ones for the glycolysis, the tricarboxylic acid cycle, the aerobic respiratory chain, and the Calvin cycle of carbon fixation. Genes for dissimilatory oxidation of reduced sulfur compounds, namely the branched SOX system (SoxAXBYZ), direct (soeABC) and indirect (aprAB, sat) pathways of sulfite oxidation, sulfur oxidation complex Dsr (dsrABEFHCEMKLJONR), sulfide oxidation systems SQR (sqrA, sqrF), and FCSD (fccAB) were found in the core genome. Genomes differ in the set of genes for dissimilatory reduction of nitrogen compounds, nitrogen fixation, and the presence of various types of RuBisCO.

Genome archaeology of two laboratory Salmonella enterica enterica sv Typhimurium.

The Salmonella research community has used strains and bacteriophages over decades, exchanging useful new isolates among laboratories for the study of cell surface antigens, metabolic pathways and restriction-modification (RM) studies. Here we present the sequences of two laboratory Salmonella strains (STK005, an isolate of LB5000; and its descendant ER3625). In the ancestry of LB5000, segments of ∼15 and ∼42 kb were introduced from Salmonella enterica sv Abony 803 into S. enterica sv Typhimurium LT2, forming strain SD14; this strain is thus a hybrid of S. enterica isolates. Strains in the SD14 lineage were used to define flagellar antigens from the 1950s to the 1970s, and to define three RM systems from the 1960s to the 1980s. LB5000 was also used as a host in phage typing systems used by epidemiologists. In the age of cheaper and easier sequencing, this resource will provide access to the sequence that underlies the extensive literature.

Rapid identification of methylase specificity (RIMS-seq) jointly identifies methylated motifs and generates shotgun sequencing of bacterial genomes.

DNA methylation is widespread amongst eukaryotes and prokaryotes to modulate gene expression and confer viral resistance. 5-Methylcytosine (m5C) methylation has been described in genomes of a large fraction of bacterial species as part of restriction-modification systems, each composed of a methyltransferase and cognate restriction enzyme. Methylases are site-specific and target sequences vary across organisms. High-throughput methods, such as bisulfite-sequencing can identify m5C at base resolution but require specialized library preparations and single molecule, real-time (SMRT) sequencing usually misses m5C. Here, we present a new method called RIMS-seq (rapid identification of methylase specificity) to simultaneously sequence bacterial genomes and determine m5C methylase specificities using a simple experimental protocol that closely resembles the DNA-seq protocol for Illumina. Importantly, the resulting sequencing quality is identical to DNA-seq, enabling RIMS-seq to substitute standard sequencing of bacterial genomes. Applied to bacteria and synthetic mixed communities, RIMS-seq reveals new methylase specificities, supporting routine study of m5C methylation while sequencing new genomes.

Complete Genome and Methylome Analysis of the Box-Shaped Halophilic Archaeon Haloarcula sinaiiensis ATCC 33800.

The genome of halophilic archaeon Haloarcula sinaiiensis ATCC 33800 was sequenced and assembled and comprises seven replicons. Four m6A and one m4C modified motifs and their responsible methyltransferase genes have been identified in the genome by single-molecule real-time (SMRT) sequencing and bioinformatic analysis.

Characterization of BisI Homologs.

BisI is a sequence-specific and 5-methylcytosine (m5C)-dependent restriction endonuclease (REase), that cleaves the modified DNA sequence Gm5CNGC (G indicates that the cytosine opposite to G is modified). We expressed and purified a number of BisI homologs from sequenced bacterial genomes and used Illumina sequencing to determine the Pam7902I (Esp638I-like) cleavage sites in phage Xp12 DNA. One BisI homolog KpnW2I is EcoBLMcrX-like, cleaving GCNGC/RCNGY/RCNRC sites with m5C. We also cloned and expressed three BisI homologs from metagenome sequences derived from thermophilic sources. One enzyme EsaTMI is active at 37 to 65°C. EsaHLI cleaves GCNGC sites with three to four m5C and is active up to 50°C. In addition, we determined the number and position of m5C in BisI sites for efficient cleavage. BisI cleavage efficiency of GCNGC site is as following: Gm5CNGC (two internal m5C) > Gm5CNGC (one internal m5C) > GCNGm5C (one external m5C) > > GCNGC (unmodified). Three or four m5C in GCNGC site also supports BisI cleavage although partial inhibition was observed on duplex oligos with four m5C. BisI can be used to partially cleave a desired GCNGC site targeted with a complementary oligonucleotide (hemi-methylated). The m5C-dependent BisI variants will be useful for epigenetic research.

Structural and functional diversity among Type III restriction-modification systems that confer host DNA protection via methylation of the N4 atom of cytosine.

We report a new subgroup of Type III Restriction-Modification systems that use m4C methylation for host protection. Recognition specificities for six such systems, each recognizing a novel motif, have been determined using single molecule real-time DNA sequencing. In contrast to all previously characterized Type III systems which modify adenine to m6A, protective methylation of the host genome in these new systems is achieved by the N4-methylation of a cytosine base in one strand of an asymmetric 4 to 6 base pair recognition motif. Type III systems are heterotrimeric enzyme complexes containing a single copy of an ATP-dependent restriction endonuclease-helicase (Res) and a dimeric DNA methyltransferase (Mod). The Type III Mods are beta-class amino-methyltransferases, examples of which form either N6-methyl adenine or N4-methyl cytosine in Type II RM systems. The Type III m4C Mod and Res proteins are diverged, suggesting ancient origin or that m4C modification has arisen from m6A MTases multiple times in diverged lineages. Two of the systems, from thermophilic organisms, required expression of both Mod and Res to efficiently methylate an E. coli host, unlike previous findings that Mod alone is proficient at modification, suggesting that the division of labor between protective methylation and restriction activities is atypical in these systems. Two of the characterized systems, and many homologous putative systems, appear to include a third protein; a conserved putative helicase/ATPase subunit of unknown function and located 5' of the mod gene. The function of this additional ATPase is not yet known, but close homologs co-localize with the typical Mod and Res genes in hundreds of putative Type III systems. Our findings demonstrate a rich diversity within Type III RM systems.

Production and сharacterization of the exopolysaccharide from strain Paenibacillus polymyxa 2020.

Paenibacillus spp. exopolysaccharides (EPSs) have become a growing interest recently as a source of biomaterials. In this study, we characterized Paenibacillus polymyxa 2020 strain, which produces a large quantity of EPS (up to 68 g/L),and was isolated from wasp honeycombs. Here we report its complete genome sequence and full methylome analysis detected by Pacific Biosciences SMRT sequencing. Moreover, bioinformatic analysis identified a putative levan synthetic operon. SacC and sacB genes have been cloned and their products identified as glycoside hydrolase and levansucrase respectively. The Fourier transform infrared (FT-IR) and nuclear magnetic resonance (NMR) spectra demonstrated that the EPS is a linear ß-(2→6)-linked fructan (levan). The structure and properties of levan polymer produced from sucrose and molasses were analyzed by FT-IR, NMR, scanning electron microscopy (SEM), high performance size exclusion chromatography (HPSEC), thermogravimetric analysis (TGA), cytotoxicity tests and showed low toxicity and high biocompatibility. Thus, P. polymyxa 2020 could be an exceptional cost-effective source for the industrial production of levan-type EPSs and to obtain functional biomaterials based on it for a broad range of applications, including bioengineering.

Expression of Human ACE2 N-terminal Domain, Part of the Receptor for SARS-CoV-2, in Fusion With Maltose-Binding Protein, E. coli Ribonuclease I and Human RNase A.

The SARS-CoV-2 viral genome contains a positive-strand single-stranded RNA of ∼30 kb. Human ACE2 protein is the receptor for SARS-CoV-2 virus attachment and infection. We propose to use ribonucleases (RNases) as antiviral agents to destroy the viral genome in vitro. In the virions, the RNA is protected by viral capsid proteins, membrane proteins, and nucleocapsid proteins. To utilize RNases as antiviral strategy, we set out to construct RNase fusion with human ACE2 receptor N-terminal domain (ACE2NTD). We expressed six proteins in E. coli cells: (1) MBP-ACE2NTD, (2) ACE2NTD-GFP, (3) RNase I (6×His), (4) RNase III (6×His), (5) RNase I-ACE2NTD (6×His), and (6) human RNase A-ACE2NTD (6×His). We evaluated fusion expression in different E. coli strains, partially purified MBP-ACE2NTD protein from the soluble fraction of bacterial cell lysate, and refolded MBP-ACE2NTD protein from inclusion body. The engineered RNase I-ACE2NTD (6×His) and hRNase A-ACE2NTD (6×His) fusions are active in cleaving SARS-CoV-2 RNA fragment in vitro. The recombinant RNase I (6×His) and RNase III (6×His) are active in cleaving RNA and dsRNA in test tube. This study provides a proof-of-concept for construction of fusion protein between human receptor and nuclease that may be used to degrade viral nucleic acids.

Genome-wide identification of 5-methylcytosine sites in bacterial genomes by high-throughput sequencing of MspJI restriction fragments.

Single-molecule Real-Time (SMRT) sequencing can easily identify sites of N6-methyladenine and N4-methylcytosine within DNA sequences, but similar identification of 5-methylcytosine sites is not as straightforward. In prokaryotic DNA, methylation typically occurs within specific sequence contexts, or motifs, that are a property of the methyltransferases that "write" these epigenetic marks. We present here a straightforward, cost-effective alternative to both SMRT and bisulfite sequencing for the determination of prokaryotic 5-methylcytosine methylation motifs. The method, called MFRE-Seq, relies on excision and isolation of fully methylated fragments of predictable size using MspJI-Family Restriction Enzymes (MFREs), which depend on the presence of 5-methylcytosine for cleavage. We demonstrate that MFRE-Seq is compatible with both Illumina and Ion Torrent sequencing platforms and requires only a digestion step and simple column purification of size-selected digest fragments prior to standard library preparation procedures. We applied MFRE-Seq to numerous bacterial and archaeal genomic DNA preparations and successfully confirmed known motifs and identified novel ones. This method should be a useful complement to existing methodologies for studying prokaryotic methylomes and characterizing the contributing methyltransferases.

Genome-wide methylome analysis of two strains belonging to the hypervirulent Neisseria meningitidis serogroup W ST-11 clonal complex.

A rising incidence of meningococcal serogroup W disease has been evident in many countries worldwide. Serogroup W isolates belonging to the sequence type (ST)-11 clonal complex have been associated with atypical symptoms and increased case fatality rates. The continued expansion of this clonal complex in the later part of the 2010s has been largely due to a shift from the so-called original UK strain to the 2013 strain. Here we used single-molecule real-time (SMRT) sequencing to determine the methylomes of the two major serogroup W strains belonging to ST-11 clonal complex. Five methylated motifs were identified in this study, and three of the motifs, namely 5'-GATC-3', 5'-GAAGG-3', 5'-GCGCGC-3', were found in all 13 isolates investigated. The results showed no strain-specific motifs or difference in active restriction modification systems between the two strains. Two phase variable methylases were identified and the enrichment or depletion of the methylation motifs generated by these methylases varied between the two strains. Results from this work give further insight into the low diversity of methylomes in highly related strains and encourage further research to decipher the role of regions with under- or overrepresented methylation motifs.

Complete Annotated Genome Sequence of the Salmonella enterica Serovar Typhimurium LT7 Strain STK003, Historically Used in Gene Transfer Studies.

The genome of Salmonella enterica serovar Typhimurium LT7 comprises a chromosome and two plasmids. One plasmid is very close to pSLT of Salmonella Typhimurium LT2; the second harbors a shufflon region. Prophage content is distinct: LT7 lacks Fels-1, while Gifsy-1 and Fels-2 show island-like divergence and likely programmed inversion, respectively.

Plasmid replication-associated single-strand-specific methyltransferases.

Analysis of genomic DNA from pathogenic strains of Burkholderia cenocepacia J2315 and Escherichia coli O104:H4 revealed the presence of two unusual MTase genes. Both are plasmid-borne ORFs, carried by pBCA072 for B. cenocepacia J2315 and pESBL for E. coli O104:H4. Pacific Biosciences SMRT sequencing was used to investigate DNA methyltransferases M.BceJIII and M.EcoGIX, using artificial constructs. Mating properties of engineered pESBL derivatives were also investigated. Both MTases yield promiscuous m6A modification of single strands, in the context SAY (where S = C or G and Y = C or T). Strikingly, this methylation is asymmetric in vivo, detected almost exclusively on one DNA strand, and is incomplete: typically, around 40% of susceptible motifs are modified. Genetic and biochemical studies suggest that enzyme action depends on replication mode: DNA Polymerase I (PolI)-dependent ColE1 and p15A origins support asymmetric modification, while the PolI-independent pSC101 origin does not. An MTase-PolI complex may enable discrimination of PolI-dependent and independent plasmid origins. M.EcoGIX helps to establish pESBL in new hosts by blocking the action of restriction enzymes, in an orientation-dependent fashion. Expression and action appear to occur on the entering single strand in the recipient, early in conjugal transfer, until lagging-strand replication creates the double-stranded form.

Protein Domain Guided Screen for Sequence Specific and Phosphorothioate-Dependent Restriction Endonucleases.

Modification dependent restriction endonucleases (MDREs) restrict modified DNA, typically with limited sequence specificity (∼2-4 bp). Here, we focus on MDREs that have an SRA and/or SBD (sulfur binding domain) fused to an HNH endonuclease domain, cleaving cytosine modified or phosphorothioated (PT) DNA. We independently characterized the SBD-SRA-HNH endonuclease ScoMcrA, which preferentially cleaves 5hmC modified DNA. We report five SBD-HNH endonucleases, all recognizing GpsAAC/GpsTTC sequence and cleaving outside with a single nucleotide 3' stagger: EcoWI (N7/N6), Ksp11411I (N5/N4), Bsp305I (N6/N4-5), Mae9806I [N(8-10)/N(8-9)], and Sau43800I [N(8-9)/N(7-8)]. EcoWI and Bsp305I are more specific for PT modified DNA in Mg2+ buffer, and promiscuous with Mn2+. Ksp11411I is more PT specific with Ni2+. EcoWI and Ksp11411I cleave fully- and hemi-PT modified oligos, while Bsp305I cleaves only fully modified ones. EcoWI forms a dimer in solution and cleaves more efficiently in the presence of two modified sites. In addition, we demonstrate that EcoWI PT-dependent activity has biological function: EcoWI expressing cells restrict dnd+ GpsAAC modified plasmid strongly, and GpsGCC DNA weakly. This work establishes a framework for biotechnology applications of PT-dependent restriction endonucleases (PTDRs).

Analysis of a phase-variable restriction modification system of the human gut symbiont Bacteroides fragilis.

The genomes of gut Bacteroidales contain numerous invertible regions, many of which contain promoters that dictate phase-variable synthesis of surface molecules such as polysaccharides, fimbriae, and outer surface proteins. Here, we characterize a different type of phase-variable system of Bacteroides fragilis, a Type I restriction modification system (R-M). We show that reversible DNA inversions within this R-M locus leads to the generation of eight specificity proteins with distinct recognition sites. In vitro grown bacteria have a different proportion of specificity gene combinations at the expression locus than bacteria isolated from the mammalian gut. By creating mutants, each able to produce only one specificity protein from this region, we identified the R-M recognition sites of four of these S-proteins using SMRT sequencing. Transcriptome analysis revealed that the locked specificity mutants, whether grown in vitro or isolated from the mammalian gut, have distinct transcriptional profiles, likely creating different phenotypes, one of which was confirmed. Genomic analyses of diverse strains of Bacteroidetes from both host-associated and environmental sources reveal the ubiquity of phase-variable R-M systems in this phylum.

Complete Genome Sequences and Methylome Analysis of Two Environmental Spirochaetes.

Here, we report the finished closed genomes of two environmental bacteria, Oceanispirochaeta crateria K2 and Thiospirochaeta perfilievii P (formally known as Spirochaeta perfilievii P). In addition, we provide methylation data and the associated enzymes predicted and confirmed to be responsible for each modified motif.

Complete genome and methylome analysis of Neisseria meningitidis associated with increased serogroup Y disease.

Invasive meningococcal disease (IMD) due to serogroup Y Neisseria meningitidis emerged in Europe during the 2000s. Draft genomes of serogroup Y isolates in Sweden revealed that although the population structure of these isolates was similar to other serogroup Y isolates internationally, a distinct strain (YI) and more specifically a sublineage (1) of this strain was responsible for the increase of serogroup Y IMD in Sweden. We performed single molecule real-time (SMRT) sequencing on eight serogroup Y isolates from different sublineages to unravel the genetic and epigenetic factors delineating them, in order to understand the serogroup Y emergence. Extensive comparisons between the serogroup Y sublineages of all coding sequences, complex genomic regions, intergenic regions, and methylation motifs revealed small point mutations in genes mainly encoding hypothetical and metabolic proteins, and non-synonymous variants in genes involved in adhesion, iron acquisition, and endotoxin production. The methylation motif CACNNNNNTAC was only found in isolates of sublineage 2. Only seven genes were putatively differentially expressed, and another two genes encoding hypothetical proteins were only present in sublineage 2. These data suggest that the serogroup Y IMD increase in Sweden was most probably due to small changes in genes important for colonization and transmission.

Structure of HhaI endonuclease with cognate DNA at an atomic resolution of 1.0 Å.

HhaI, a Type II restriction endonuclease, recognizes the symmetric sequence 5'-GCG↓C-3' in duplex DNA and cleaves ('↓') to produce fragments with 2-base, 3'-overhangs. We determined the structure of HhaI in complex with cognate DNA at an ultra-high atomic resolution of 1.0 Å. Most restriction enzymes act as dimers with two catalytic sites, and cleave the two strands of duplex DNA simultaneously, in a single binding event. HhaI, in contrast, acts as a monomer with only one catalytic site, and cleaves the DNA strands sequentially, one after the other. HhaI comprises three domains, each consisting of a mixed five-stranded ß sheet with a defined function. The first domain contains the catalytic-site; the second contains residues for sequence recognition; and the third contributes to non-specific DNA binding. The active-site belongs to the 'PD-D/EXK' superfamily of nucleases and contains the motif SD-X11-EAK. The first two domains are similar in structure to two other monomeric restriction enzymes, HinP1I (G↓CGC) and MspI (C↓CGG), which produce fragments with 5'-overhangs. The third domain, present only in HhaI, shifts the positions of the recognition residues relative to the catalytic site enabling this enzyme to cleave the recognition sequence at a different position. The structure of M.HhaI, the biological methyltransferase partner of HhaI, was determined earlier. Together, these two structures represent the first natural pair of restriction-modification enzymes to be characterized in atomic detail.