Characterization of type II and III restriction-modification systems from Bacillus cereus strains ATCC 10987 and ATCC 14579

Deciphering microbial gene function using natural language processing

Revealing the function of uncharacterized genes is a fundamental challenge in an era of ever-increasing volumes of sequencing data. Here, we present a concept for tackling this challenge using deep learning methodologies adopted from natural language processing (NLP). We repurpose NLP algorithms to model "gene semantics" based on a biological corpus of more than 360 million microbial genes within their genomic context. We use the language models to predict functional categories for 56,617 genes and find that out of 1369 genes associated with recently discovered defense systems, 98% are inferred correctly. We then systematically evaluate the "discovery potential" of different functional categories, pinpointing those with the most genes yet to be characterized. Finally, we demonstrate our method's ability to discover systems associated with microbial interaction and defense. Our results highlight that combining microbial genomics and language models is a promising avenue for revealing gene functions in microbes.

A GCDGC-specific DNA (cytosine-5) methyltransferase that methylates the GCWGC sequence on both strands and the GCSGC sequence on one strand

5-Methylcytosine is one of the major epigenetic marks of DNA in living organisms. Some bacterial species possess DNA methyltransferases that modify cytosines on both strands to produce fully-methylated sites or on either strand to produce hemi-methylated sites. In this study, we characterized a DNA methyltransferase that produces two sequences with different methylation patterns: one methylated on both strands and another on one strand. M.BatI is the orphan DNA methyltransferase of Bacillus anthracis coded in one of the prophages on the chromosome. Analysis of M.BatI modified DNA by bisulfite sequencing revealed that the enzyme methylates the first cytosine in sequences of 5ʹ-GCAGC-3ʹ, 5ʹ-GCTGC-3ʹ, and 5ʹ-GCGGC-3ʹ, but not of 5ʹ-GCCGC-3ʹ. This resulted in the production of fully-methylated 5ʹ-GCWGC-3ʹ and hemi-methylated 5ʹ-GCSGC-3ʹ. M.BatI also showed toxicity when expressed in E. coli, which was caused by a mechanism other than DNA modification activity. Homologs of M.BatI were found in other Bacillus species on different prophage like regions, suggesting the spread of the gene by several different phages. The discovery of the DNA methyltransferase with unique modification target specificity suggested unrevealed diversity of target sequences of bacterial cytosine DNA methyltransferase.

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.

An investigation of Burkholderia cepacia complex methylomes via SMRT sequencing and mutant analysis

Bacterial genomes can be methylated at particular motifs by methyltransferases (M). This DNA modification allows restriction endonucleases (R) to discriminate between self and foreign DNA. While the accepted primary function of such restriction modification (RM) systems is to degrade incoming foreign DNA, other roles of RM systems and lone R or M components have been found in genome protection, stability and the regulation of various phenotypes. The Burkholderia cepacia complex (Bcc) is a group of closely related opportunistic pathogens with biotechnological potential. Here, we constructed and analysed mutants lacking various RM components in the clinical Bcc isolate Burkholderia cenocepacia H111 and used SMRT sequencing of single mutants to assign the B. cenocepacia H111 Ms to their cognate motifs. DNA methylation is shown to affect biofilm formation, cell shape, motility, siderophore production and membrane vesicle production. Moreover, DNA methylation had a large effect on the maintenance of the Bcc virulence megaplasmid pC3. Our data also suggest that the gp51 M-encoding gene, which is essential in H111 and is located within a prophage, is required for maintaining the bacteriophage in a lysogenic state, thereby ensuring a constant, low level of phage production within the bacterial population.ImportanceWhile genome sequence determines an organism's proteins, methylation of the nucleotides themselves can confer additional properties. In bacteria, Ms modify specific nucleotide motifs to allow discrimination of 'self' from 'non-self' DNA, e.g. from bacteriophages. Restriction enzymes detect 'non-self' methylation patterns and cut foreign DNA. Furthermore, methylation of promoter regions can influence gene expression and hence affect various phenotypes. In this study, we determined the methylated motifs of four strains from the Burkholderia cepacia complex of opportunistic pathogens. We deleted all genes encoding the restriction and modification components in one of these strains, Burkholderia cenocepacia H111. It is shown that DNA methylation affects various phenotypic traits, the most noteworthy being lysogenicity of a bacteriophage and maintenance of a virulence megaplasmid.

In vitro Type II Restriction of Bacteriophage DNA With Modified Pyrimidines

To counteract host-encoded restriction systems, bacteriophages (phages) incorporate modified bases in their genomes. For example, phages carry in their genomes modified pyrimidines such as 5-hydroxymethyl-cytosine (5hmC) in T4gt deficient in α- and β-glycosyltransferases, glucosylated-5-hydroxymethylcytosine (5gmC) in T4, 5-methylcytosine (5mC) in Xp12, and 5-hydroxymethyldeoxyuridine (5hmdU) in SP8. In this work we sequenced phage Xp12 and SP8 genomes and examined Type II restriction of T4gt, T4, Xp12, and SP8 phage DNAs. T4gt, T4, and Xp12 genomes showed resistance to 81.9% (186 out of 227 enzymes tested), 94.3% (214 out of 227 enzymes tested), and 89.9% (196 out of 218 enzymes tested), respectively, commercially available Type II restriction endonucleases (REases). The SP8 genome, however, was resistant to only ∼8.3% of these enzymes (17 out of 204 enzymes tested). SP8 DNA could be further modified by adenine DNA methyltransferases (MTases) such as M.Dam and M.EcoGII as well as a number of cytosine DNA MTases, such as CpG methylase. The 5hmdU base in SP8 DNA was phosphorylated by treatment with a 5hmdU DNA kinase to achieve ∼20% phosphorylated 5hmdU, resulting resistance or partially resistant to more Type II restriction. This work provides a convenient reference for molecular biologists working with modified pyrimidines and using REases. The genomic sequences of phage Xp12 and SP8 lay the foundation for further studies on genetic pathways for 5mC and 5hmdU DNA base modifications and for comparative phage genomics.

Type II Restriction of Bacteriophage DNA With 5hmdU-Derived Base Modifications

To counteract bacterial defense systems, bacteriophages (phages) make extensive base modifications (substitutions) to block endonuclease restriction. Here we evaluated Type II restriction of three thymidine (T or 5-methyldeoxyuridine, 5mdU) modified phage genomes: Pseudomonas phage M6 with 5-(2-aminoethyl)deoxyuridine (5-NedU), Salmonella phage ViI (Vi1) with 5-(2-aminoethoxy)methyldeoxyuridine (5-NeOmdU) and Delftia phage phi W-14 (a.k.a. ΦW-14) with α-putrescinylthymidine (putT). Among >200 commercially available restriction endonucleases (REases) tested, phage M6, ViI, and phi W-14 genomic DNAs (gDNA) show resistance against 48.4, 71.0, and 68.8% of Type II restrictions, respectively. Inspection of the resistant sites indicates the presence of conserved dinucleotide TG or TC (TS, S=C, or G), implicating the specificity of TS sequence as the target that is converted to modified base in the genomes. We also tested a number of DNA methyltransferases (MTases) on these phage DNAs and found some MTases can fully or partially modify the DNA to confer more resistance to cleavage by REases. Phage M6 restriction fragments can be efficiently ligated by T4 DNA ligase. Phi W-14 restriction fragments show apparent reduced rate in E. coli exonuclease III degradation. This work extends previous studies that hypermodified T derived from 5hmdU provides additional resistance to host-encoded restrictions, in parallel to modified cytosines, guanine, and adenine in phage genomes. The results reported here provide a general guidance to use REases to map and clone phage DNA with hypermodified thymidine.

A novel T4SS-mediated DNA transfer used by pXO16, a conjugative plasmid from Bacillus thuringiensis serovar israelensis

The entomopathogenic Bacillus thuringiensis serovar israelensis displays peculiar conjugative transfer capabilities, accounted for by the large conjugative plasmid pXO16 (350 kb). The efficient and fast conjugative transfers are accompanied by a macroscopic aggregation of bacterial partners. Moreover, pXO16 has proven capable of effective mobilization and the retro-transfer of both mobilizable and 'non-mobilizable' plasmids. In this work, the aggregation phenomenon is shown to promote pXO16 transfer while not being mandatory for transfer. Transfer of pXO16 to B. thuringiensis recipient strains that do not display aggregation is observed as well, hence enlarging the previously defined host range. The use of variant calling analysis of transconjugants allowed for observation of up to 791 kb chromosomal regions mobilization. Previous analysis of pXO16 did not reveal any Type IV Secretion System (T4SS) homologs, which suggested the presence of an unusual conjugative system. A FtsK/SpOIIIE ATPase gene proved here to be necessary for conjugative transfer. Additionally, the analysis of natural restriction-modification systems in both conjugative partners gave credit to a ssDNA transfer mechanism. A 'transfer israelensis plasmid' (tip) region containing this ATPase gene was shown to code for other potential T4SS proteins, illustrating a conjugative system distantly related to the other known Gram-positive T4SSs.

A Novel Tool for Microbial Genome Editing Using the Restriction-Modification System

Scarless genetic manipulation of genomes is an essential tool for biological research. The restriction-modification (R-M) system is a defense system in bacteria that protects against invading genomes on the basis of its ability to distinguish foreign DNA from self DNA. Here, we designed an R-M system-mediated genome editing (RMGE) technique for scarless genetic manipulation in different microorganisms. For bacteria with Type IV REase, an RMGE technique using the inducible DNA methyltransferase gene, bceSIIM (RMGE-bceSIIM), as the counter-selection cassette was developed to edit the genome of Escherichia coli. For bacteria without Type IV REase, an RMGE technique based on a restriction endonuclease (RMGE-mcrA) was established in Bacillus subtilis. These techniques were successfully used for gene deletion and replacement with nearly 100% counter-selection efficiencies, which were higher and more stable compared to conventional methods. Furthermore, precise point mutation without limiting sites was achieved in E. coli using RMGE-bceSIIM to introduce a single base mutation of A128C into the rpsL gene. In addition, the RMGE-mcrA technique was applied to delete the CAN1 gene in Saccharomyces cerevisiae DAY414 with 100% counter-selection efficiency. The effectiveness of the RMGE technique in E. coli, B. subtilis, and S. cerevisiae suggests the potential universal usefulness of this technique for microbial genome manipulation.

CgII cleaves DNA using a mechanism distinct from other ATP-dependent restriction endonucleases

The restriction endonuclease CglI from Corynebacterium glutamicum recognizes an asymmetric 5′-GCCGC-3′ site and cleaves the DNA 7 and 6/7 nucleotides downstream on the top and bottom DNA strands, respectively, in an NTP-hydrolysis dependent reaction. CglI is composed of two different proteins: an endonuclease (R.CglI) and a DEAD-family helicase-like ATPase (H.CglI). These subunits form a heterotetrameric complex with R₂H₂ stoichiometry. However, the R₂H_2·CglI complex has only one nuclease active site sufficient to cut one DNA strand suggesting that two complexes are required to introduce a double strand break. Here, we report studies to evaluate the DNA cleavage mechanism of CglI. Using one- and two-site circular DNA substrates we show that CglI does not require two sites on the same DNA for optimal catalytic activity. However, one-site linear DNA is a poor substrate, supporting a mechanism where CglI complexes must communicate along the one-dimensional DNA contour before cleavage is activated. Based on experimental data, we propose that adenosine triphosphate (ATP) hydrolysis by CglI produces translocation on DNA preferentially in a downstream direction from the target, although upstream translocation is also possible. Our results are consistent with a mechanism of CglI action that is distinct from that of other ATP-dependent restriction-modification enzymes.

Polyvalent Proteins, a Pervasive Theme in the Intergenomic Biological Conflicts of Bacteriophages and Conjugative Elements

Intense biological conflicts between prokaryotic genomes and their genomic parasites have resulted in an arms race in terms of the molecular “weaponry” deployed on both sides. Using a recursive computational approach, we uncovered a remarkable class of multidomain proteins with 2 to 15 domains in the same polypeptide deployed by viruses and plasmids in such conflicts. Domain architectures and genomic contexts indicate that they are part of a widespread conflict strategy involving proteins injected into the host cell along with parasite DNA during the earliest phase of infection. Their unique feature is the combination of domains with highly disparate biochemical activities in the same polypeptide; accordingly, we term them polyvalent proteins. Of the 131 domains in polyvalent proteins, a large fraction are enzymatic domains predicted to modify proteins, target nucleic acids, alter nucleotide signaling/metabolism, and attack peptidoglycan or cytoskeletal components. They further contain nucleic acid-binding domains, virion structural domains, and 40 novel uncharacterized domains. Analysis of their architectural network reveals both pervasive common themes and specialized strategies for conjugative elements and plasmids or (pro)phages. The themes include likely processing of multidomain polypeptides by zincin-like metallopeptidases and mechanisms to counter restriction or CRISPR/Cas systems and jump-start transcription or replication. DNA-binding domains acquired by eukaryotes from such systems have been reused in XPC/RAD4-dependent DNA repair and mitochondrial genome replication in kinetoplastids. Characterization of the novel domains discovered here, such as RNases and peptidases, are likely to aid in the development of new reagents and elucidation of the spread of antibiotic resistance.

IMPORTANCE This is the first report of the widespread presence of large proteins, termed polyvalent proteins, predicted to be transmitted by genomic parasites such as conjugative elements, plasmids, and phages during the initial phase of infection along with their DNA. They are typified by the presence of multiple domains with disparate activities combined in the same protein. While some of these domains are predicted to assist the invasive element in replication, transcription, or protection of their DNA, several are likely to target various host defense systems or modify the host to favor the parasite's life cycle. Notably, DNA-binding domains from these systems have been transferred to eukaryotes, where they have been incorporated into DNA repair and mitochondrial genome replication systems.

Biosynthesis and Function of Modified Bases in Bacteria and Their Viruses

Naturally occurring modification of the canonical A, G, C, and T bases can be found in the DNA of cellular organisms and viruses from all domains of life. Bacterial viruses (bacteriophages) are a particularly rich but still underexploited source of such modified variant nucleotides. The modifications conserve the coding and base-pairing functions of DNA, but add regulatory and protective functions. In prokaryotes, modified bases appear primarily to be part of an arms race between bacteriophages (and other genomic parasites) and their hosts, although, as in eukaryotes, some modifications have been adapted to convey epigenetic information. The first half of this review catalogs the identification and diversity of DNA modifications found in bacteria and bacteriophages. What is known about the biogenesis, context, and function of these modifications are also described. The second part of the review places these DNA modifications in the context of the arms race between bacteria and bacteriophages. It focuses particularly on the defense and counter-defense strategies that turn on direct recognition of the presence of a modified base. Where modification has been shown to affect other DNA transactions, such as expression and chromosome segregation, that is summarized, with reference to recent reviews.

Expression and purification of a single-chain Type IV restriction enzyme Eco94GmrSD and determination of its substrate preference

The first reported Type IV restriction endonuclease (REase) GmrSD consists of GmrS and GmrD subunits. In most bacteria, however, the gmrS and gmrD genes are fused together to encode a single-chain protein. The fused coding sequence for ECSTEC94C_1402 from E. coli strain STEC_94C was expressed in T7 Express. The protein designated as Eco94GmrSD displays modification-dependent ATP-stimulated REase activity on T4 DNA with glucosyl-5-hydroxymethyl-cytosines (glc-5hmC) and T4gt DNA with 5-hydroxymethyl-cytosines (5hmC). A C-terminal 6xHis-tagged protein was purified by two-column chromatography. The enzyme is active in Mg²⁺ and Mn²⁺ buffer. It prefers to cleave large glc-5hmC- or 5hmC-modified DNA. In phage restriction assays, Eco94GmrSD weakly restricted T4 and T4gt, whereas T4 IPI*-deficient phage (Δip1) were restricted more than 10⁶-fold, consistent with IPI* protection of E. coli DH10B from lethal expression of the closely homologous E. coli CT596 GmrSD. Eco94GmrSD is proposed to belong to the His-Asn-His (HNH)-nuclease family by the identification of a putative C-terminal REase catalytic site D507-H508-N522. Supporting this, GmrSD variants D507A, H508A, and N522A displayed no endonuclease activity. The presence of a large number of fused GmrSD homologs suggests that GmrSD is an effective phage exclusion protein that provides a mechanism to thwart T-even phage infection.

Biochemical functional predictions for protein structures of unknown or uncertain function

With the exponential growth in the determination of protein sequences and structures via genome sequencing and structural genomics efforts, there is a growing need for reliable computational methods to determine the biochemical function of these proteins. This paper reviews the efforts to address the challenge of annotating the function at the molecular level of uncharacterized proteins. While sequence- and three-dimensional-structure-based methods for protein function prediction have been reviewed previously, the recent trends in local structure-based methods have received less attention. These local structure-based methods are the primary focus of this review. Computational methods have been developed to predict the residues important for catalysis and the local spatial arrangements of these residues can be used to identify protein function. In addition, the combination of different types of methods can help obtain more information and better predictions of function for proteins of unknown function. Global initiatives, including the Enzyme Function Initiative (EFI), COMputational BRidges to EXperiments (COMBREX), and the Critical Assessment of Function Annotation (CAFA), are evaluating and testing the different approaches to predicting the function of proteins of unknown function. These initiatives and global collaborations will increase the capability and reliability of methods to predict biochemical function computationally and will add substantial value to the current volume of structural genomics data by reducing the number of absent or inaccurate functional annotations.

Endonuclease specificity and sequence dependence of type IIS restriction enzymes

Restriction enzymes that recognize specific sequences but cleave unknown sequence outside the recognition site are extensively utilized tools in molecular biology. Despite this, systematic functional categorization of cleavage performance has largely been lacking. We established a simple and automatable model system to assay cleavage distance variation (termed slippage) and the sequence dependence thereof. We coupled this to massively parallel sequencing in order to provide sensitive and accurate measurement. With this system 14 enzymes were assayed (AcuI, BbvI, BpmI, BpuEI, BseRI, BsgI, Eco57I, Eco57MI, EcoP15I, FauI, FokI, GsuI, MmeI and SmuI). We report significant variation of slippage ranging from 1–54%, variations in sequence context dependence, as well as variation between isoschizomers. We believe this largely overlooked property of enzymes with shifted cleavage would benefit from further large scale classification and engineering efforts seeking to improve performance. The gained insights of in-vitro performance may also aid the in-vivo understanding of these enzymes.

DNA cleavage by CgII and NgoAVII requires interaction between N- and R-proteins and extensive nucleotide hydrolysis

The stress-sensitive restriction-modification (RM) system CglI from Corynebacterium glutamicum and the homologous NgoAVII RM system from Neisseria gonorrhoeae FA1090 are composed of three genes: a DNA methyltransferase (M.CglI and M.NgoAVII), a putative restriction endonuclease (R.CglI and R.NgoAVII, or R-proteins) and a predicted DEAD-family helicase/ATPase (N.CglI and N.NgoAVII or N-proteins). Here we report a biochemical characterization of the R- and N-proteins. Size-exclusion chromatography and SAXS experiments reveal that the isolated R.CglI, R.NgoAVII and N.CglI proteins form homodimers, while N.NgoAVII is a monomer in solution. Moreover, the R.CglI and N.CglI proteins assemble in a complex with R2N2 stoichiometry. Next, we show that N-proteins have ATPase activity that is dependent on double-stranded DNA and is stimulated by the R-proteins. Functional ATPase activity and extensive ATP hydrolysis (∼170 ATP/s/monomer) are required for site-specific DNA cleavage by R-proteins. We show that ATP-dependent DNA cleavage by R-proteins occurs at fixed positions (6-7 nucleotides) downstream of the asymmetric recognition sequence 5'-GCCGC-3'. Despite similarities to both Type I and II restriction endonucleases, the CglI and NgoAVII enzymes may employ a unique catalytic mechanism for DNA cleavage.

Unlocking the mystery of the hard-to-sequence phage genome: PaP1 methylome and bacterial immunity

Background

Whole-genome sequencing is an important method to understand the genetic information, gene function, biological characteristics and survival mechanisms of organisms. Sequencing large genomes is very simple at present. However, we encountered a hard-to-sequence genome of Pseudomonas aeruginosa phage PaP1. Shotgun sequencing method failed to complete the sequence of this genome.

Results

After persevering for 10 years and going over three generations of sequencing techniques, we successfully completed the sequence of the PaP1 genome with a length of 91,715 bp. Single-molecule real-time sequencing results revealed that this genome contains 51 N-6-methyladenines and 152 N-4-methylcytosines. Three significant modified sequence motifs were predicted, but not all of the sites found in the genome were methylated in these motifs. Further investigations revealed a novel immune mechanism of bacteria, in which host bacteria can recognise and repel modified bases containing inserts in a large scale. This mechanism could be accounted for the failure of the shotgun method in PaP1 genome sequencing. This problem was resolved using the nfi^- mutant of Escherichia coli DH5α as a host bacterium to construct a shotgun library.

Conclusions

This work provided insights into the hard-to-sequence phage PaP1 genome and discovered a new mechanism of bacterial immunity. The methylome of phage PaP1 is responsible for the failure of shotgun sequencing and for bacterial immunity mediated by enzyme Endo V activity; this methylome also provides a valuable resource for future studies on PaP1 genome replication and modification, as well as on gene regulation and host interaction.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-803) contains supplementary material, which is available to authorized users.

Comprehensive analysis of the HEPN superfamily: identification of novel roles in intra-genomic conflicts, defense, pathogenesis and RNA processing

Background

The major role of enzymatic toxins that target nucleic acids in biological conflicts at all levels has become increasingly apparent thanks in large part to the advances of comparative genomics. Typically, toxins evolve rapidly hampering the identification of these proteins by sequence analysis. Here we analyze an unexpectedly widespread superfamily of toxin domains most of which possess RNase activity.

Results

The HEPN superfamily is comprised of all α-helical domains that were first identified as being associated with DNA polymerase β-type nucleotidyltransferases in prokaryotes and animal Sacsin proteins. Using sensitive sequence and structure comparison methods, we vastly extend the HEPN superfamily by identifying numerous novel families and by detecting diverged HEPN domains in several known protein families. The new HEPN families include the RNase LS and LsoA catalytic domains, KEN domains (e.g. RNaseL and Ire1) and the RNase domains of RloC and PrrC. The majority of HEPN domains contain conserved motifs that constitute a metal-independent endoRNase active site. Some HEPN domains lacking this motif probably function as non-catalytic RNA-binding domains, such as in the case of the mannitol repressor MtlR. Our analysis shows that HEPN domains function as toxins that are shared by numerous systems implicated in intra-genomic, inter-genomic and intra-organismal conflicts across the three domains of cellular life. In prokaryotes HEPN domains are essential components of numerous toxin-antitoxin (TA) and abortive infection (Abi) systems and in addition are tightly associated with many restriction-modification (R-M) and CRISPR-Cas systems, and occasionally with other defense systems such as Pgl and Ter. We present evidence of multiple modes of action of HEPN domains in these systems, which include direct attack on viral RNAs (e.g. LsoA and RNase LS) in conjunction with other RNase domains (e.g. a novel RNase H fold domain, NamA), suicidal or dormancy-inducing attack on self RNAs (RM systems and possibly CRISPR-Cas systems), and suicidal attack coupled with direct interaction with phage components (Abi systems). These findings are compatible with the hypothesis on coupling of pathogen-targeting (immunity) and self-directed (programmed cell death and dormancy induction) responses in the evolution of robust antiviral strategies. We propose that altruistic cell suicide mediated by HEPN domains and other functionally similar RNases was essential for the evolution of kin and group selection and cell cooperation. HEPN domains were repeatedly acquired by eukaryotes and incorporated into several core functions such as endonucleolytic processing of the 5.8S-25S/28S rRNA precursor (Las1), a novel ER membrane-associated RNA degradation system (C6orf70), sensing of unprocessed transcripts at the nuclear periphery (Swt1). Multiple lines of evidence suggest that, similar to prokaryotes, HEPN proteins were recruited to antiviral, antitransposon, apoptotic systems or RNA-level response to unfolded proteins (Sacsin and KEN domains) in several groups of eukaryotes.

Conclusions

Extensive sequence and structure comparisons reveal unexpectedly broad presence of the HEPN domain in an enormous variety of defense and stress response systems across the tree of life. In addition, HEPN domains have been recruited to perform essential functions, in particular in eukaryotic rRNA processing. These findings are expected to stimulate experiments that could shed light on diverse cellular processes across the three domains of life.

Reviewers

This article was reviewed by Martijn Huynen, Igor Zhulin and Nick Grishin

Natural zinc ribbon HNH endonucleases and engineered zinc finger nicking endonuclease

Many bacteriophage and prophage genomes encode an HNH endonuclease (HNHE) next to their cohesive end site and terminase genes. The HNH catalytic domain contains the conserved catalytic residues His-Asn-His and a zinc-binding site [CxxC]₂. An additional zinc ribbon (ZR) domain with one to two zinc-binding sites ([CxxxxC], [CxxxxH], [CxxxC], [HxxxH], [CxxC] or [CxxH]) is frequently found at the N-terminus or C-terminus of the HNHE or a ZR domain protein (ZRP) located adjacent to the HNHE. We expressed and purified 10 such HNHEs and characterized their cleavage sites. These HNHEs are site-specific and strand-specific nicking endonucleases (NEase or nickase) with 3- to 7-bp specificities. A minimal HNH nicking domain of 76 amino acid residues was identified from Bacillus phage γ HNHE and subsequently fused to a zinc finger protein to generate a chimeric NEase with a new specificity (12–13 bp). The identification of a large pool of previously unknown natural NEases and engineered NEases provides more ‘tools’ for DNA manipulation and molecular diagnostics. The small modular HNH nicking domain can be used to generate rare NEases applicable to targeted genome editing. In addition, the engineered ZF nickase is useful for evaluation of off-target sites in vitro before performing cell-based gene modification.

The methylomes of six bacteria

Six bacterial genomes, Geobacter metallireducens GS-15, Chromohalobacter salexigens, Vibrio breoganii 1C-10, Bacillus cereus ATCC 10987, Campylobacter jejuni subsp. jejuni 81-176 and C. jejuni NCTC 11168, all of which had previously been sequenced using other platforms were re-sequenced using single-molecule, real-time (SMRT) sequencing specifically to analyze their methylomes. In every case a number of new N(6)-methyladenine ((m6)A) and N(4)-methylcytosine ((m4)C) methylation patterns were discovered and the DNA methyltransferases (MTases) responsible for those methylation patterns were assigned. In 15 cases, it was possible to match MTase genes with MTase recognition sequences without further sub-cloning. Two Type I restriction systems required sub-cloning to differentiate their recognition sequences, while four MTase genes that were not expressed in the native organism were sub-cloned to test for viability and recognition sequences. Two of these proved active. No attempt was made to detect 5-methylcytosine ((m5)C) recognition motifs from the SMRT® sequencing data because this modification produces weaker signals using current methods. However, all predicted (m6)A and (m4)C MTases were detected unambiguously. This study shows that the addition of SMRT sequencing to traditional sequencing approaches gives a wealth of useful functional information about a genome showing not only which MTase genes are active but also revealing their recognition sequences.

A mimicking-of-DNA-methylation-patterns pipeline for overcoming the restriction barrier of bacteria

Genetic transformation of bacteria harboring multiple Restriction-Modification (R-M) systems is often difficult using conventional methods. Here, we describe a mimicking-of-DNA-methylation-patterns (MoDMP) pipeline to address this problem in three difficult-to-transform bacterial strains. Twenty-four putative DNA methyltransferases (MTases) from these difficult-to-transform strains were cloned and expressed in an Escherichia coli strain lacking all of the known R-M systems and orphan MTases. Thirteen of these MTases exhibited DNA modification activity in Southwestern dot blot or Liquid Chromatography–Mass Spectrometry (LC–MS) assays. The active MTase genes were assembled into three operons using the Saccharomyces cerevisiae DNA assembler and were co-expressed in the E. coli strain lacking known R-M systems and orphan MTases. Thereafter, results from the dot blot and restriction enzyme digestion assays indicated that the DNA methylation patterns of the difficult-to-transform strains are mimicked in these E. coli hosts. The transformation of the Gram-positive Bacillus amyloliquefaciens TA208 and B. cereus ATCC 10987 strains with the shuttle plasmids prepared from MoDMP hosts showed increased efficiencies (up to four orders of magnitude) compared to those using the plasmids prepared from the E. coli strain lacking known R-M systems and orphan MTases or its parental strain. Additionally, the gene coding for uracil phosphoribosyltransferase (upp) was directly inactivated using non-replicative plasmids prepared from the MoDMP host in B. amyloliquefaciens TA208. Moreover, the Gram-negative chemoautotrophic Nitrobacter hamburgensis strain X14 was transformed and expressed Green Fluorescent Protein (GFP). Finally, the sequence specificities of active MTases were identified by restriction enzyme digestion, making the MoDMP system potentially useful for other strains. The effectiveness of the MoDMP pipeline in different bacterial groups suggests a universal potential. This pipeline could facilitate the functional genomics of the strains that are difficult to transform.

Approximately 95% of the genome-sequenced bacteria harbor Restriction-Modification (R-M) systems. R-M systems usually occur in pairs, i.e., DNA methyltransferases (MTases) and restriction endonucleases (REases). REases can degrade invading DNA to protect the cell from infection by phages. This protecting machinery has also become the barrier for experimental genetic manipulation, because the newly introduced DNA would be degraded by the REases of the transformed bacteria. In this study we have developed a pipeline to protect DNA by methylation from cleavage by host REases. Multiple DNA MTases were cloned from three difficult-to-transform bacterial strains and co-expressed in an E. coli strain lacking all of the known endogenous R-M systems and orphan MTases. Thus, the DNA methylation patterns of these strains have become similar to that of the difficult-to-transform strains. Ultimately, the DNA prepared from these E. coli strains can overcome the R-M barrier of the bacterial strains that are difficult to transform and achieve genetic manipulation. The effectiveness of this pipeline in different bacterial groups suggests a universal potential. This pipeline could facilitate functional genomics of bacterial strains that are difficult to transform.