dbAPIS: a database of anti-prokaryotic immune system genes.

Anti-prokaryotic immune system (APIS) proteins, typically encoded by phages, prophages, and plasmids, inhibit prokaryotic immune systems (e.g. restriction modification, toxin-antitoxin, CRISPR-Cas). A growing number of APIS genes have been characterized and dispersed in the literature. Here we developed dbAPIS (https://bcb.unl.edu/dbAPIS), as the first literature curated data repository for experimentally verified APIS genes and their associated protein families. The key features of dbAPIS include: (i) experimentally verified APIS genes with their protein sequences, functional annotation, PDB or AlphaFold predicted structures, genomic context, sequence and structural homologs from different microbiome/virome databases; (ii) classification of APIS proteins into sequence-based families and construction of hidden Markov models (HMMs); (iii) user-friendly web interface for data browsing by the inhibited immune system types or by the hosts, and functions for searching and batch downloading of pre-computed data; (iv) Inclusion of all types of APIS proteins (except for anti-CRISPRs) that inhibit a variety of prokaryotic defense systems (e.g. RM, TA, CBASS, Thoeris, Gabija). The current release of dbAPIS contains 41 verified APIS proteins and ∼4400 sequence homologs of 92 families and 38 clans. dbAPIS will facilitate the discovery of novel anti-defense genes and genomic islands in phages, by providing a user-friendly data repository and a web resource for an easy homology search against known APIS proteins.

Restriction-Modification systems interplay causes avoidance of GATC site in prokaryotic genomes.

Palindromes are frequently underrepresented in prokaryotic genomes. Palindromic 5[Formula: see text]-GATC-3[Formula: see text] site is a recognition site of different Restriction-Modification (R-M) systems, as well as solitary methyltransferase Dam. Classical GATC-specific R-M systems methylate GATC and cleave unmethylated GATC. On the contrary, methyl-directed Type II restriction endonucleases cleave methylated GATC. Methylation of GATC by Dam methyltransferase is involved in the regulation of different cellular processes. The diversity of functions of GATC-recognizing proteins makes GATC sequence a good model for studying the reasons of palindrome avoidance in prokaryotic genomes. In this work, the influence of R-M systems and solitary proteins on the GATC site avoidance is described by a mathematical model. GATC avoidance is strongly associated with the presence of alternate (methyl-directed or classical Type II R-M system) genes in different strains of the same species, as we have shown for Streptococcus pneumoniae, Neisseria meningitidis, Eubacterium rectale, and Moraxella catarrhalis. We hypothesize that GATC avoidance can result from a DNA exchange between strains with different methylation status of GATC site within the process of natural transformation. If this hypothesis is correct, the GATC avoidance is a sign of a DNA exchange between bacteria with different methylation status in a mixed population.

The Epigenomic Landscape of Prokaryotes.

DNA methylation acts in concert with restriction enzymes to protect the integrity of prokaryotic genomes. Studies in a limited number of organisms suggest that methylation also contributes to prokaryotic genome regulation, but the prevalence and properties of such non-restriction-associated methylation systems remain poorly understood. Here, we used single molecule, real-time sequencing to map DNA modifications including m6A, m4C, and m5C across the genomes of 230 diverse bacterial and archaeal species. We observed DNA methylation in nearly all (93%) organisms examined, and identified a total of 834 distinct reproducibly methylated motifs. This data enabled annotation of the DNA binding specificities of 620 DNA Methyltransferases (MTases), doubling known specificities for previously hard to study Type I, IIG and III MTases, and revealing their extraordinary diversity. Strikingly, 48% of organisms harbor active Type II MTases with no apparent cognate restriction enzyme. These active 'orphan' MTases are present in diverse bacterial and archaeal phyla and show motif specificities and methylation patterns consistent with functions in gene regulation and DNA replication. Our results reveal the pervasive presence of DNA methylation throughout the prokaryotic kingdoms, as well as the diversity of sequence specificities and potential functions of DNA methylation systems.

To be or not to be: regulation of restriction-modification systems and other toxin-antitoxin systems.

One of the simplest classes of genes involved in programmed death is that containing the toxin-antitoxin (TA) systems of prokaryotes. These systems are composed of an intracellular toxin and an antitoxin that neutralizes its effect. These systems, now classified into five types, were initially discovered because some of them allow the stable maintenance of mobile genetic elements in a microbial population through postsegregational killing or the death of cells that have lost these systems. Here, we demonstrate parallels between some TA systems and restriction-modification systems (RM systems). RM systems are composed of a restriction enzyme (toxin) and a modification enzyme (antitoxin) and limit the genetic flux between lineages with different epigenetic identities, as defined by sequence-specific DNA methylation. The similarities between these systems include their postsegregational killing and their effects on global gene expression. Both require the finely regulated expression of a toxin and antitoxin. The antitoxin (modification enzyme) or linked protein may act as a transcriptional regulator. A regulatory antisense RNA recently identified in an RM system can be compared with those RNAs in TA systems. This review is intended to generalize the concept of TA systems in studies of stress responses, programmed death, genetic conflict and epigenetics.

DNA cleavage and methylation specificity of the single polypeptide restriction-modification enzyme LlaGI.

LlaGI is a single polypeptide restriction-modification enzyme encoded on the naturally-occurring plasmid pEW104 isolated from Lactococcus lactis ssp. cremoris W10. Bioinformatics analysis suggests that the enzyme contains domains characteristic of an mrr endonuclease, a superfamily 2 DNA helicase and a gamma-family adenine methyltransferase. LlaGI was expressed and purified from a recombinant clone and its properties characterised. An asymmetric recognition sequence was identified, 5'-CTnGAyG-3' (where n is A, G, C or T and y is C or T). Methylation of the recognition site occurred on only one strand (the non-degenerate dA residue of 5'-CrTCnAG-3' being methylated at the N6 position). Double strand DNA breaks at distant, random sites were only observed when two head-to-head oriented, unmethylated copies of the site were present; single sites or pairs in tail-to-tail or head-to-tail repeat only supported a DNA nicking activity. dsDNA nuclease activity was dependent upon the presence of ATP or dATP. Our results are consistent with a directional long-range communication mechanism that is necessitated by the partial site methylation. In the accompanying manuscript [Smith et al. (2009) The single polypeptide restriction-modification enzyme LlaGI is a self-contained molecular motor that translocates DNA loops], we demonstrate that this communication is via 1-dimensional DNA loop translocation. On the basis of this data and that in the third accompanying manuscript [Smith et al. (2009) An Mrr-family nuclease motif in the single polypeptide restriction-modification enzyme LlaGI], we propose that LlaGI is the prototype of a new sub-classification of Restriction-Modification enzymes, named Type I SP (for Single Polypeptide).

A rapid and efficient method for cloning genes of type II restriction-modification systems by use of a killer plasmid.

We present a method for cloning restriction-modification (R-M) systems that is based on the use of a lethal plasmid (pKILLER). The plasmid carries a functional gene for a restriction endonuclease having the same DNA specificity as the R-M system of interest. The first step is the standard preparation of a representative, plasmid-borne genomic library. Then this library is transformed with the killer plasmid. The only surviving bacteria are those which carry the gene specifying a protective DNA methyltransferase. Conceptually, this in vivo selection approach resembles earlier methods in which a plasmid library was selected in vitro by digestion with a suitable restriction endonuclease, but it is much more efficient than those methods. The new method was successfully used to clone two R-M systems, BstZ1II from Bacillus stearothermophilus 14P and Csp231I from Citrobacter sp. strain RFL231, both isospecific to the prototype HindIII R-M system.

Crystal structure of a putative type I restriction-modification S subunit from Mycoplasma genitalium.

The crystal structure of the eubacteria Mycoplasma genitalium ORF MG438 polypeptide, determined by multiple anomalous dispersion and refined at 2.3 A resolution, reveals the organization of S subunits from the Type I restriction and modification system. The structure consists of two globular domains, with about 150 residues each, separated by a pair of 40 residue long antiparallel alpha-helices. The globular domains correspond to the variable target recognition domains (TRDs), as previously defined for S subunits on sequence analysis, while the two helices correspond to the central (CR1) and C-terminal (CR2) conserved regions, respectively. The structure of the MG438 subunit presents an overall cyclic topology with an intramolecular 2-fold axis that superimposes the N and the C-half parts, each half containing a globular domain and a conserved helix. TRDs are found to be structurally related with the small domain of the Type II N6-adenine DNA MTase TaqI. These relationships together with the structural peculiarities of MG438, in particular the presence of the intramolecular quasi-symmetry, allow the proposal of a model for S subunits recognition of their DNA targets in agreement with previous experimental results. In the crystal, two subunits of MG438 related by a crystallographic 2-fold axis present a large contact area mainly involving the symmetric interactions of a cluster of exposed hydrophobic residues. Comparison with the recently reported structure of an S subunit from the archaea Methanococcus jannaschii highlights the structural features preserved despite a sequence identity below 20%, but also reveals important differences in the globular domains and in their disposition with respect to the conserved regions.

The biology of restriction and anti-restriction.

The phenomena of prokaryotic restriction and modification, as well as anti-restriction, were first discovered five decades ago but have yielded only gradually to rigorous analysis. Work presented at the 5th New England Biolabs Meeting on Restriction-Modification (available on REBASE, http://www.rebase.com) and several recently published genetic, biochemical and biophysical analyses indicate that these fields continue to contribute significantly to basic science. Recently, there have been several studies that have shed light on the still developing field of restriction-modification and on the newly re-emerging field of anti-restriction.

Structure of the Q237W mutant of HhaI DNA methyltransferase: an insight into protein-protein interactions.

We have determined the structure of a mutant (Q237W) of HhaI DNA methyltransferase, complexed with the methyl-donor product AdoHcy. The Q237W mutant proteins were crystallized in the monoclinic space group C2 with two molecules in the crystallographic asymmetric unit. Protein-protein interface calculations in the crystal lattices suggest that the dimer interface has the specific characteristics for homodimer protein-protein interactions, while the two active sites are spatially independent on the outer surface of the dimer. The solution behavior suggests the formation of HhaI dimers as well. The same HhaI dimer interface is also observed in the previously characterized binary (M.HhaI-AdoMet) and ternary (M.HhaI-DNA-AdoHcy) complex structures, crystallized in different space groups. The dimer is characterized either by a non-crystallographic two-fold symmetry or a crystallographic symmetry. The dimer interface involves three segments: the amino-terminal residues 2-8, the carboxy-terminal residues 313-327, and the linker (amino acids 179-184) between the two functional domains--the catalytic methylation domain and the DNA target recognition domain. Both the amino- and carboxy-terminal segments are part of the methylation domain. We also examined protein-protein interactions of other structurally characterized DNA MTases, which are often found as a 2-fold related 'dimer' with the largest dimer interface area for the group-beta MTases. A possible evolutionary link between the Type I and Type II restriction-modification systems is discussed.

Restriction and modification systems of ruminal bacteria.

A high frequency of type II restriction endonuclease activities was detected in Selenomonas ruminantium but not in other rumen bacteria tested. Eight different restriction endonucleases were characterized in 17 strains coming from genetically homogeneous local population. Chromosomal DNA isolated from S. ruminantium strains was found to be refractory to cleavage by various restriction enzymes, implying the presence of methylase activities additional to those required for protection against the cellular endonucleases. The presence of Dam methylation was detected in S. ruminantium strains as well as in several other species belonging to the Sporomusa subbranch of low G + C Gram-positive bacteria (Megasphaera elsdenii, Mitsuokella multiacidus).

In vivo restriction by LlaI is encoded by three genes, arranged in an operon with llaIM, on the conjugative Lactococcus plasmid pTR2030.

The LlaI restriction and modification (R/M) system is encoded on pTR2030, a 46.2-kb conjugative plasmid from Lactococcus lactis. The llaI methylase gene, sequenced previously, encodes a functional type IIS methylase and is located approximately 5 kb upstream from the abiA gene, encoding abortive phage resistance. In this study, the sequence of the region between llaIM and abiA was determined and revealed four consecutive open reading frames (ORFs). Northern (RNA) analysis showed that the four ORFs were part of a 7-kb operon with llaIM and the downstream abiA gene on a separate transcriptional unit. The deduced protein sequence of ORF2 revealed a P-loop consensus motif for ATP/GTP-binding sites and a three-part consensus motif for GTP-binding proteins. Data bank searches with the deduced protein sequences for all four ORFs revealed no homology except for ORF2 with MerB, in three regions that coincided with the GTP-binding motifs in both proteins. To phenotypically analyze the llaI operon, a 9.0-kb fragment was cloned into a high-copy-number lactococcal shuttle vector, pTRKH2. The resulting construct, pTRK370, exhibited a significantly higher level of in vivo restriction and modification in L. lactis NCK203 than the low-copy-number parental plasmid, pTR2030. A combination of deletion constructions and frameshift mutations indicated that the first three ORFs were involved in LlaI restriction, and they were therefore designated llaI.1, llaI.2, and llaI.3. Mutating llaI.1 completely abolished restriction, while disrupting llaI.2 or llaI.3 allowed an inefficient restriction of phage DNA to occur, manifested primarily by a variable plaque phenotype. ORF4 had no discernible effect on in vivo restriction. A frameshift mutation in llaIM proved lethal to L. lactis NCK203, implying that the restriction component was active without the modification subunit. These results suggested that the LlaI R/M system is unlike any other R/M system studied to date and has diverged from the type IIS class of restriction enzymes by acquiring some characteristics reminiscent of type I enzymes.