Conflicts targeting epigenetic systems and their resolution by cell death: novel concepts for methyl-specific and other restriction systems

Dependence of post-segregational killing mediated by Type II restriction-modification systems on the lifetime of restriction endonuclease effective activity

Plasmid-borne Type II restriction-modification (RM) systems mediate post-segregational killing (PSK). PSK is thought to be caused by the dilution of restriction and modification enzymes during cell division, resulting in accumulation of unmethylated DNA recognition sites and their cleavage by restriction endonucleases. PSK is the likely reason for stabilization of plasmids carrying RM systems in the absence of selection for plasmid maintenance. In this study, we developed a CRISPR interference-based method to eliminate RM-carrying plasmids and study PSK-related phenomena with minimal perturbation to the Escherichia coli host. Plasmids carrying the EcoRV, Eco29kI, and EcoRI RM systems were highly stable, and their loss resulted in SOS response and PSK. In contrast, plasmids carrying the Esp1396I system were poorly stabilized; their loss led to a temporary cessation of growth, followed by full recovery. We demonstrate that this unusual behavior is due to a limited lifetime of the Esp1396I restriction endonuclease activity, which, upon Esp1396I plasmid loss, disappears approximately after two cycles of cell division, i.e., before unmethylated sites appear in significant numbers. Our results indicate that whenever PSK induced by a loss of RM systems, and, possibly, other toxin-antitoxin systems, is considered, the lifetimes of individual system components and the growth rate of host cells shall be taken in account. Mathematical modeling shows, that unlike the situation with classical toxin-antitoxin systems, RM system-mediated PSK is possible when the lifetimes of restriction endonuclease and methyltransferase activities are similar, as long as the toxic restriction endonuclease activity persists for more than two chromosome replication cycles.IMPORTANCEIt is widely accepted that many Type II restriction-modification (RM) systems mediate post-segregational killing (PSK) if plasmids that encode them are lost. In this study, we harnessed an inducible CRISPR-Cas system to remove RM plasmids from Escherichia coli cells to study PSK while minimally perturbing cell physiology. We demonstrate that PSK depends on restriction endonuclease activity lifetime and is not observed when it is less than two replication cycles. We present a mathematical model that explains experimental data and shows that unlike the case of toxin-antitoxin-mediated PSK, the loss of an RM system induced PSK even when the RM enzymes have identical lifetimes.

Functionally comparable but evolutionarily distinct nucleotide-targeting effectors help identify conserved paradigms across diverse immune systems

While nucleic acid-targeting effectors are known to be central to biological conflicts and anti-selfish element immunity, recent findings have revealed immune effectors that target their building blocks and the cellular energy currency-free nucleotides. Through comparative genomics and sequence-structure analysis, we identified several distinct effector domains, which we named Calcineurin-CE, HD-CE, and PRTase-CE. These domains, along with specific versions of the ParB and MazG domains, are widely present in diverse prokaryotic immune systems and are predicted to degrade nucleotides by targeting phosphate or glycosidic linkages. Our findings unveil multiple potential immune systems associated with at least 17 different functional themes featuring these effectors. Some of these systems sense modified DNA/nucleotides from phages or operate downstream of novel enzymes generating signaling nucleotides. We also uncovered a class of systems utilizing HSP90- and HSP70-related modules as analogs of STAND and GTPase domains that are coupled to these nucleotide-targeting- or proteolysis-induced complex-forming effectors. While widespread in bacteria, only a limited subset of nucleotide-targeting effectors was integrated into eukaryotic immune systems, suggesting barriers to interoperability across subcellular contexts. This work establishes nucleotide-degrading effectors as an emerging immune paradigm and traces their origins back to homologous domains in housekeeping systems.

A host of armor: Prokaryotic immune strategies against mobile genetic elements

Prokaryotic adaptation is strongly influenced by the horizontal acquisition of beneficial traits via mobile genetic elements (MGEs), such as viruses/bacteriophages and plasmids. However, MGEs can also impose a fitness cost due to their often parasitic nature and differing evolutionary trajectories. In response, prokaryotes have evolved diverse immune mechanisms against MGEs. Recently, our understanding of the abundance and diversity of prokaryotic immune systems has greatly expanded. These defense systems can degrade the invading genetic material, inhibit genome replication, or trigger abortive infection, leading to population protection. In this review, we highlight these strategies, focusing on the most recent discoveries. The study of prokaryotic defenses not only sheds light on microbial evolution but also uncovers novel enzymatic activities with promising biotechnological applications.

Antimicrobial resistance and mechanisms of epigenetic regulation

The rampant use of antibiotics in animal husbandry, farming and clinical disease treatment has led to a significant issue with pathogen resistance worldwide over the past decades. The classical mechanisms of resistance typically investigate antimicrobial resistance resulting from natural resistance, mutation, gene transfer and other processes. However, the emergence and development of bacterial resistance cannot be fully explained from a genetic and biochemical standpoint. Evolution necessitates phenotypic variation, selection, and inheritance. There are indications that epigenetic modifications also play a role in antimicrobial resistance. This review will specifically focus on the effects of DNA modification, histone modification, rRNA methylation and the regulation of non-coding RNAs expression on antimicrobial resistance. In particular, we highlight critical work that how DNA methyltransferases and non-coding RNAs act as transcriptional regulators that allow bacteria to rapidly adapt to environmental changes and control their gene expressions to resist antibiotic stress. Additionally, it will delve into how Nucleolar-associated proteins in bacteria perform histone functions akin to eukaryotes. Epigenetics, a non-classical regulatory mechanism of bacterial resistance, may offer new avenues for antibiotic target selection and the development of novel antibiotics.

OUP accepted manuscript

Epigenetic DNA methylation plays an important role in bacteria by influencing gene expression and allowing discrimination between self-DNA and intruders such as phages and plasmids. Restriction–modification (RM) systems use a methyltransferase (MTase) to modify a specific sequence motif, thus protecting host DNA from cleavage by a cognate restriction endonuclease (REase) while leaving invading DNA vulnerable. Other REases occur solitarily and cleave methylated DNA. REases and RM systems are frequently mobile, influencing horizontal gene transfer by altering the compatibility of the host for foreign DNA uptake. However, whether mobile defence systems affect pre-existing host defences remains obscure. Here, we reveal an epigenetic conflict between an RM system (PcaRCI) and a methylation-dependent REase (PcaRCII) in the plant pathogen Pectobacterium carotovorum RC5297. The PcaRCI RM system provides potent protection against unmethylated plasmids and phages, but its methylation motif is targeted by the methylation-dependent PcaRCII. This potentially lethal co-existence is enabled through epigenetic silencing of the PcaRCII-encoding gene via promoter methylation by the PcaRCI MTase. Comparative genome analyses suggest that the PcaRCII-encoding gene was already present and was silenced upon establishment of the PcaRCI system. These findings provide a striking example for selfishness of RM systems and intracellular competition between different defences.

DNA cytosine methyltransferase enhances viability during prolonged stationary phase in Escherichia coli

In Escherichia coli, DNA cytosine methyltransferase (Dcm) methylates the second cytosine in the sequence 5'CCWGG3' generating 5-methylcytosine. Dcm is not associated with a cognate restriction enzyme, suggesting Dcm impacts facets of bacterial physiology outside of restriction-modification systems. Other than gene expression changes, there are few phenotypes that have been identified in strains with natural or engineered Dcm loss, and thus Dcm function has remained an enigma. Herein, we demonstrate that Dcm does not impact bacterial growth under optimal and selected stress conditions. However, Dcm does impact viability in long-term stationary phase competition experiments. Dcm+ cells outcompete cells lacking dcm under different conditions. Dcm knockout cells have more RpoS-dependent HPII catalase activity than wild-type cells. Thus, the impact of Dcm on stationary phase may involve changes in RpoS activity. Overall, our data reveal a new role for Dcm during long-term stationary phase.

Comparative Methylome Analysis of Campylobacter jejuni Strain YH002 Reveals a Putative Novel Motif and Diverse Epigenetic Regulations of Virulence Genes

Campylobacter jejuni is a major cause of foodborne gastroenteritis worldwide inflicting palpable socioeconomic costs. The ability of this pathogen to successfully infect its hosts is determined not only by the presence of specific virulence genes but also by the pathogen’s capacity to appropriately regulate those virulence genes. Therefore, DNA methylation can play a critical role in both aspects of this process because it serves as both a means to protect the integrity of the cellular DNA from invasion and as a mechanism to control transcriptional regulation within the cell. In the present study we report the comparative methylome data of C. jejuni YH002, a multidrug resistant strain isolated from retail beef liver. Investigation into the methylome identified a putative novel motif (CGCGA) of a type II restriction-modification (RM) system. Comparison of methylomes of the strain to well-studied C. jejuni strains highlighted non-uniform methylation patterns among the strains though the existence of the typical type I and type IV RM systems were also observed. Additional investigations into the existence of DNA methylation sites within gene promoters, which may ultimately result in altered levels of transcription, revealed several virulence genes putatively regulated using this mode of action. Of those identified, a flagella gene (flhB), a RNA polymerase sigma factor (rpoN), a capsular polysaccharide export protein (kpsD), and a multidrug efflux pump were highly notable.

The N-terminal domain of Staphylothermus marinus McrB shares structural homology with PUA-like RNA binding proteins

McrBC is a conserved modification-dependent restriction system that in Escherichia coli specifically targets foreign DNA containing methylated cytosines. Crystallographic data show that the N-terminal domain of Escherichia coli McrB binds substrates via a base flipping mechanism. This region is poorly conserved among the plethora of McrB homologs, suggesting that other species may use alternative binding strategies and/or recognize different targets. Here we present the crystal structure of the N-terminal domain from Stayphlothermus marinus McrB (Sm3-180) at 1.92 Å, which adopts a PUA-like EVE fold that is closely related to the YTH and ASCH RNA binding domains. Unlike most PUA-like domains, Sm3-180 binds DNA and can associate with different modified substrates. We find the canonical 'aromatic cage' binding pocket that confers specificity for methylated bases in other EVE/YTH domains is degenerate and occluded in Sm3-180, which may contribute to its promiscuity in target recognition. Further structural comparison between different PUA-like domains identifies motifs and conformational variations that correlate with the preference for binding either DNA or RNA. Together these data have important implications for PUA-like domain specificity and suggest a broader biological versatility for the McrBC family than previously described.

Highly regulated, diversifying NTP-dependent biological conflict systems with implications for the emergence of multicellularity

Social cellular aggregation or multicellular organization pose increased risk of transmission of infections through the system upon infection of a single cell. The generality of the evolutionary responses to this outside of Metazoa remains unclear. We report the discovery of several thematically unified, remarkable biological conflict systems preponderantly present in multicellular prokaryotes. These combine thresholding mechanisms utilizing NTPase chaperones (the MoxR-vWA couple), GTPases and proteolytic cascades with hypervariable effectors, which vary either by using a reverse transcriptase-dependent diversity-generating system or through a system of acquisition of diverse protein modules, typically in inactive form, from various cellular subsystems. Conciliant lines of evidence indicate their deployment against invasive entities, like viruses, to limit their spread in multicellular/social contexts via physical containment, dominant-negative interactions or apoptosis. These findings argue for both a similar operational 'grammar' and shared protein domains in the sensing and limiting of infections during the multiple emergences of multicellularity.

The structure of the Thermococcus gammatolerans McrB N-terminal domain reveals a new mode of substrate recognition and specificity among McrB homologs

McrBC is a two-component, modification-dependent restriction system that cleaves foreign DNA-containing methylated cytosines. Previous crystallographic studies have shown that Escherichia coli McrB uses a base-flipping mechanism to recognize these modified substrates with high affinity. The side chains stabilizing both the flipped base and the distorted duplex are poorly conserved among McrB homologs, suggesting that other mechanisms may exist for binding modified DNA. Here we present the structures of the Thermococcus gammatolerans McrB DNA-binding domain (TgΔ185) both alone and in complex with a methylated DNA substrate at 1.68 and 2.27 Å resolution, respectively. The structures reveal that TgΔ185 consists of a YT521-B homology (YTH) domain, which is commonly found in eukaryotic proteins that bind methylated RNA and is structurally unrelated to the E. coli McrB DNA-binding domain. Structural superposition and co-crystallization further show that TgΔ185 shares a conserved aromatic cage with other YTH domains, which forms the binding pocket for a flipped-out base. Mutational analysis of this aromatic cage supports its role in conferring specificity for the methylated adenines, whereas an extended basic surface present in TgΔ185 facilitates its preferential binding to duplex DNA rather than RNA. Together, these findings establish a new binding mode and specificity among McrB homologs and expand the biological roles of YTH domains.

Methylation Warfare: Interaction of Pneumococcal Bacteriophages with Their Host

With antimicrobial drug resistance becoming an increasing burden on human health, much attention has been focused on the potential use of bacteriophages and their enzymes as therapeutics. However, the investigations into the physiology of the complex interactions of bacteriophages with their hosts have attracted far less attention, in comparison. This work describes the molecular characterization of the infectious cycle of a bacteriophage in the important human pathogen Streptococcus pneumoniae and explores the intricate relationship between phase-variable host defense mechanisms and the virus. This is the first report showing how a phase-variable type I restriction-modification system is involved in bacteriophage restriction while it also provides an additional level of infection control through abortive infection.

Virus-host interactions are regulated by complex coevolutionary dynamics. In Streptococcus pneumoniae, phase-variable type I restriction-modification (R-M) systems are part of the core genome. We hypothesized that the ability of the R-M systems to switch between six target DNA specificities also has a key role in preventing the spread of bacteriophages. Using the streptococcal temperate bacteriophage SpSL1, we show that the variants of both the SpnIII and SpnIV R-M systems are able to restrict invading bacteriophage with an efficiency approximately proportional to the number of target sites in the bacteriophage genome. In addition to restriction of lytic replication, SpnIII also led to abortive infection in the majority of host cells. During lytic infection, transcriptional analysis found evidence of phage-host interaction through the strong upregulation of the nrdR nucleotide biosynthesis regulon. During lysogeny, the phage had less of an effect on host gene regulation. This research demonstrates a novel combined bacteriophage restriction and abortive infection mechanism, highlighting the importance that the phase-variable type I R-M systems have in the multifunctional defense against bacteriophage infection in the respiratory pathogen S. pneumoniae.

IMPORTANCE With antimicrobial drug resistance becoming an increasing burden on human health, much attention has been focused on the potential use of bacteriophages and their enzymes as therapeutics. However, the investigations into the physiology of the complex interactions of bacteriophages with their hosts have attracted far less attention, in comparison. This work describes the molecular characterization of the infectious cycle of a bacteriophage in the important human pathogen Streptococcus pneumoniae and explores the intricate relationship between phase-variable host defense mechanisms and the virus. This is the first report showing how a phase-variable type I restriction-modification system is involved in bacteriophage restriction while it also provides an additional level of infection control through abortive infection.

Unexpected Evolution of Lesion-Recognition Modules in Eukaryotic NER and Kinetoplast DNA Dynamics Proteins from Bacterial Mobile Elements

The provenance of several components of major uniquely eukaryotic molecular machines are increasingly being traced back to prokaryotic biological conflict systems. Here, we demonstrate that the N-terminal single-stranded DNA-binding domain from the anti-restriction protein ArdC, deployed by bacterial mobile elements against their host, was independently acquired twice by eukaryotes, giving rise to the DNA-binding domains of XPC/Rad4 and the Tc-38-like proteins in the stem kinetoplastid. In both instances, the ArdC-N domain tandemly duplicated forming an extensive DNA-binding interface. In XPC/Rad4, the ArdC-N domains (BHDs) also fused to the inactive transglutaminase domain of a peptide-N-glycanase ultimately derived from an archaeal conflict system. Alongside, we delineate several parallel acquisitions from conjugative elements/bacteriophages that gave rise to key components of the kinetoplast DNA (kDNA) replication apparatus. These findings resolve two outstanding questions in eukaryote biology: (1) the origin of the unique DNA lesion-recognition component of NER and (2) origin of the unusual, plasmid-like features of kDNA.

The SCO1731 methyltransferase modulates actinorhodin production and morphological differentiation of Streptomyces coelicolor A3(2)

Streptomyces coelicolor is a Gram-positive microorganism often used as a model of physiological and morphological differentiation in streptomycetes, prolific producers of secondary metabolites with important biological activities. In the present study, we analysed Streptomyces coelicolor growth and differentiation in the presence of the hypo-methylating agent 5′-aza-2′-deoxycytidine (5-aza-dC) in order to investigate whether cytosine methylation has a role in differentiation. We found that cytosine demethylation caused a delay in spore germination, aerial mycelium development, sporulation, as well as a massive impairment of actinorhodin production. Thus, we searched for putative DNA methyltransferase genes in the genome and constructed a mutant of the SCO1731 gene. The analysis of the SCO1731::Tn5062 mutant strain demonstrated that inactivation of SCO1731 leads to a strong decrease of cytosine methylation and almost to the same phenotype obtained after 5-aza-dC treatment. Altogether, our data demonstrate that cytosine methylation influences morphological differentiation and actinorhodin production in S. coelicolor and expand our knowledge on this model bacterial system.

Genome content, metabolic pathways and biotechnological potential of the psychrophilic Arctic bacterium Psychrobacter sp. DAB_AL43B, a source and a host of novel Psychrobacter-specific vectors

Psychrobacter sp. DAB_AL43B, isolated from ornithogenic soil collected on the Arctic island of Spitsbergen, is a newly sequenced psychrophilic strain susceptible to conjugation and electrotransformation. Its genome consists of a circular chromosome (3.3 Mb) and four plasmids (4.4-6.4kb). In silico genome mining and microarray-based phenotypic analysis were performed to describe the metabolic potential of this strain and identify possible biotechnological applications. Metabolic reconstruction indicated that DAB_AL43B prefers low-molecular-weight carboxylates and amino acids as carbon and energy sources. Genetic determinants of heavy-metal resistance, anthracene degradation and possible aerobic denitrification were also identified. Comparative analyses revealed a relatively close relationship between DAB_AL43B and other sequenced Psychrobacter species. In addition, the plasmids of this strain were used as the basis for the construction of Escherichia coli-Psychrobacter spp. shuttle vectors. Taken together, the results of this work suggest that DAB_AL43B is a promising candidate as a new model strain for studies on Psychrobacter spp.

The Evolution of Antiviral Defense Systems

Self-replicating genetic material presumably provided the architecture necessary for generating the last universal ancestor of all nucleic-acid-based life. As biological complexity increased in the billions of years that followed, the same genetic material also morphed into a wide spectrum of viruses and other parasitic genetic elements. The resulting struggle for existence drove the evolution of host defenses, giving rise to a perpetual arms race. This Perspective summarizes the antiviral mechanisms evident across the tree of life, discussing each in their evolutionary context to postulate how the coevolution of host and pathogen shaped the cellular antiviral defenses we know today.

Identification of DNA Methyltransferase Genes in Human Pathogenic Bacteria by Comparative Genomics

DNA methylation plays an important role in gene expression and virulence in some pathogenic bacteria. In this report, we describe DNA methyltransferases (MTases) present in human pathogenic bacteria and compared them with related species, which are not pathogenic or less pathogenic, based in comparative genomics. We performed a search in the KEGG database of the KEGG database orthology groups associated with adenine and cytosine DNA MTase activities (EC: 2.1.1.37, EC: 2.1.1.113 and EC: 2.1.1.72) in 37 human pathogenic species and 18 non/less pathogenic relatives and performed comparisons of the number of these MTases sequences according to their genome size, the DNA MTase type and with their non-less pathogenic relatives. We observed that Helicobacter pylori and Neisseria spp. presented the highest number of MTases while ten different species did not present a predicted DNA MTase. We also detected a significant increase of adenine MTases over cytosine MTases (2.19 vs. 1.06, respectively, p < 0.001). Adenine MTases were the only MTases associated with restriction modification systems and DNA MTases associated with type I restriction modification systems were more numerous than those associated with type III restriction modification systems (0.84 vs. 0.17, p < 0.001); additionally, there was no correlation with the genome size and the total number of DNA MTases, indicating that the number of DNA MTases is related to the particular evolution and lifestyle of specific species, regulating the expression of virulence genes in some pathogenic bacteria.

Streptacidiphilus toruniensis sp. nov., isolated from a pine forest soil

Two acidophilic actinobacteria, isolates NA14 and NF37^T, were the subject of a polyphasic taxonomic study. Chemotaxonomic and morphological properties of the isolates were characteristic of the genus Streptacidiphilus. The isolates were shown to have identical 16S rRNA gene sequences and to be closely related to Streptacidiphilus neutrinimicus DSM 41755^T (>99.9 %). However, DNA:DNA relatedness between isolate NF37^T and the type strain of S. neutrinimicus was found to be low at 11.1 (±3.5) %. A broad range of phenotypic features were shown to distinguish the isolates from their close phylogenetic neighbours. These data shown that the isolates form a novel species of Streptacidiphilus for which the name Streptacidiphilus toruniensis sp. nov. is proposed. The type strain is NF37^T (= DSM 102291^T = NCIMB 15025^T).

Transposons to toxins: the provenance, architecture and diversification of a widespread class of eukaryotic effectors

Enzymatic effectors targeting nucleic acids, proteins and other cellular components are the mainstay of conflicts across life forms. Using comparative genomics we identify a large class of eukaryotic proteins, which include effectors from oomycetes, fungi and other parasites. The majority of these proteins have a characteristic domain architecture with one of several N-terminal 'Header' domains, which are predicted to play a role in trafficking of these effectors, including a novel version of the Ubiquitin fold. The Headers are followed by one or more diverse C-terminal domains, such as restriction endonuclease (REase), protein kinase, HNH endonuclease, LK-nuclease (a RNase) and multiple distinct peptidase domains, which are predicted to carry their toxicity determinants. The most common types of these proteins appear to have originated from prokaryotic transposases (e.g. TN7 and Mu) and combine a CDC6/ORC1-STAND clade NTPase domain with a C-terminal REase domain. Other than the so-called Crinkler effectors of oomycetes and fungi, these effectors are encoded by other eukaryotic parasites such as trypanosomatids (the RHS proteins) and the rhizarian Plasmodiophora, and symbionts like Capsaspora Remarkably, we also find these proteins in free-living eukaryotes, including several viridiplantae, fungi, amoebozoans and animals. These versions might either still be transposons or function in other poorly understood eukaryote-specific inter-organismal and inter-genomic conflicts. These include the Medea1 selfish element of Tribolium that spreads via post-zygotic killing. We present a unified mechanism for the recombination-dependent diversification and action of this widespread class of molecular weaponry deployed across diverse conflicts ranging from parasitic to free-living forms.

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.

DNA Methylation Assessed by SMRT Sequencing Is Linked to Mutations in Neisseria meningitidis Isolates

The Gram-negative bacterium Neisseria meningitidis features extensive genetic variability. To present, proposed virulence genotypes are also detected in isolates from asymptomatic carriers, indicating more complex mechanisms underlying variable colonization modes of N. meningitidis. We applied the Single Molecule, Real-Time (SMRT) sequencing method from Pacific Biosciences to assess the genome-wide DNA modification profiles of two genetically related N. meningitidis strains, both of serogroup A. The resulting DNA methylomes revealed clear divergences, represented by the detection of shared and of strain-specific DNA methylation target motifs. The positional distribution of these methylated target sites within the genomic sequences displayed clear biases, which suggest a functional role of DNA methylation related to the regulation of genes. DNA methylation in N. meningitidis has a likely underestimated potential for variability, as evidenced by a careful analysis of the ORF status of a panel of confirmed and predicted DNA methyltransferase genes in an extended collection of N. meningitidis strains of serogroup A. Based on high coverage short sequence reads, we find phase variability as a major contributor to the variability in DNA methylation. Taking into account the phase variable loci, the inferred functional status of DNA methyltransferase genes matched the observed methylation profiles. Towards an elucidation of presently incompletely characterized functional consequences of DNA methylation in N. meningitidis, we reveal a prominent colocalization of methylated bases with Single Nucleotide Polymorphisms (SNPs) detected within our genomic sequence collection. As a novel observation we report increased mutability also at 6mA methylated nucleotides, complementing mutational hotspots previously described at 5mC methylated nucleotides. These findings suggest a more diverse role of DNA methylation and Restriction-Modification (RM) systems in the evolution of prokaryotic genomes.

Classification of prokaryotic genetic replicators: between selfishness and altruism

Prokaryotes harbor a variety of genetic replicators, including plasmids, viruses, and chromosomes, each having different effects on the phenotype of the hosting cell. Here, we propose a classification for replicators of bacteria and archaea on the basis of their horizontal-transfer potential and the type of relationships (mutualistic, symbiotic, commensal, or parasitic) that they have with the host cell vehicle. Horizontal movement of replicators can be either active or passive, reflecting whether or not the replicator encodes the means to mediate its own transfer from one cell to another. Some replicators also have an infectious extracellular state, thus separating viruses from other mobile elements. From the perspective of the cell vehicle, the different types of replicators form a continuum from genuinely mutualistic to completely parasitic replicators. This classification provides a general framework for dissecting prokaryotic systems into evolutionarily meaningful components.

Restriction-modification system with methyl-inhibited base excision and abasic-site cleavage activities

The restriction-modification systems use epigenetic modification to distinguish between self and nonself DNA. A modification enzyme transfers a methyl group to a base in a specific DNA sequence while its cognate restriction enzyme introduces breaks in DNA lacking this methyl group. So far, all the restriction enzymes hydrolyze phosphodiester bonds linking the monomer units of DNA. We recently reported that a restriction enzyme (R.PabI) of the PabI superfamily with half-pipe fold has DNA glycosylase activity that excises an adenine base in the recognition sequence (5′-GTAC). We now found a second activity in this enzyme: at the resulting apurinic/apyrimidinic (AP) (abasic) site (5′-GT#C, # = AP), its AP lyase activity generates an atypical strand break. Although the lyase activity is weak and lacks sequence specificity, its covalent DNA–R.PabI reaction intermediates can be trapped by NaBH₄ reduction. The base excision is not coupled with the strand breakage and yet causes restriction because the restriction enzyme action can impair transformation ability of unmethylated DNA even in the absence of strand breaks in vitro. The base excision of R.PabI is inhibited by methylation of the target adenine base. These findings expand our understanding of genetic and epigenetic processes linking those in prokaryotes and eukaryotes.

A comparative analysis of methylome profiles of Campylobacter jejuni sheep abortion isolate and gastroenteric strains using PacBio data

Campylobacter jejuni is a leading cause of human gastrointestinal disease and small ruminant abortions in the United States. The recent emergence of a highly virulent, tetracycline-resistant C. jejuni subsp. jejuni sheep abortion clone (clone SA) in the United States, and that strain's association with human disease, has resulted in a heightened awareness of the zoonotic potential of this organism. Pacific Biosciences' Single Molecule, Real-Time sequencing technology was used to explore the variation in the genome-wide methylation patterns of the abortifacient clone SA (IA3902) and phenotypically distinct gastrointestinal-specific C. jejuni strains (NCTC 11168 and 81-176). Several notable differences were discovered that distinguished the methylome of IA3902 from that of 11168 and 81-176: identification of motifs novel to IA3902, genome-specific hypo- and hypermethylated regions, strain level variability in genes methylated, and differences in the types of methylation motifs present in each strain. These observations suggest a possible role of methylation in the contrasting disease presentations of these three C. jejuni strains. In addition, the methylation profiles between IA3902 and a luxS mutant were explored to determine if variations in methylation patterns could be identified that might explain the role of LuxS-dependent methyl recycling in IA3902 abortifacient potential.

PsasM2I, a type II restriction-modification system in Pseudomonas savastanoi pv. savastanoi: differential distribution of carrier strains in the environment and the evolutionary history of homologous RM systems in the Pseudomonas syringae complex

A type II restriction-modification system was found in a native plasmid of Pseudomonas savastanoi pv. savastanoi MLLI2. Functional analysis of the methyltransferase showed that the enzyme acts by protecting the DNA sequence CTGCAG from cleavage. Restriction endonuclease expression in recombinant Escherichia coli cells resulted in mutations in the REase sequence or transposition of insertion sequence 1A in the coding sequence, preventing lethal gene expression. Population screening detected homologous RM systems in other P. savastanoi strains and in the Pseudomonas syringae complex. An epidemiological survey carried out by sampling olive and oleander knots in two Italian regions showed an uneven diffusion of carrier strains, whose presence could be related to a selective advantage in maintaining the RM system in particular environments or subpopulations. Moreover, carrier strains can coexist in the same orchards, plants, and knot tissues with non-carriers, revealing unexpected genetic variability on a very small spatial scale. Phylogenetic analysis of the RM system and housekeeping gene sequences in the P. syringae complex demonstrated the ancient acquisition of the RM systems. However, the evolutionary history of the gene complex also showed the involvement of horizontal gene transfer between related strains and recombination events.

Dysregulation of genome-wide gene expression and DNA methylation in abnormal cloned piglets

Background

Epigenetic modifications (especially altered DNA methylation) resulting in altered gene expression may be one reason for development failure or abnormalities in cloned animals, but the underlying mechanism of the abnormal phenotype in cloned piglets remains unknown. Some cloned piglets in our study showed abnormal phenotypes such as large tongue (longer and thicker), weak muscles, and exomphalos. Here we conducted DNA methylation (DNAm) immunoprecipitation and high throughput sequencing (MeDIP-seq) and RNA sequencing (RNA-seq) of muscle tissues of cloned piglets to investigate the relationship of abnormal DNAm with gene dysregulation and the unusual phenotypes in cloned piglets.

Results

Analysis of the methylomes revealed that abnormal cloned piglets suffered more hypomethylation than hypermethylation compared to the normal cloned piglets, although the DNAm level in the CpG Island was higher in the abnormal cloned piglets. Some repetitive elements, such as SINE/tRNA-Glu Satellite/centr also showed differences. We detected 1,711 differentially expressed genes (DEGs) between the two groups, of which 243 genes also changed methylation level in the abnormal cloned piglets. The altered DNA methylation mainly affected the low and silently expressed genes. There were differences in both pathways and genes, such as the MAPK signalling pathway, the hypertrophic cardiomyopathy pathway, and the imprinted gene PLAGL1; all of which may play important roles in development of the abnormal phenotype.

Conclusions

The abnormal cloned piglets showed substantial changes both in the DNAm and the gene expression. Our data may provide new insights into understanding the molecular mechanisms of the reprogramming of genetic information in cloned animals.

Electronic supplementary material

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

Bacterial genome instability

Bacterial genomes are remarkably stable from one generation to the next but are plastic on an evolutionary time scale, substantially shaped by horizontal gene transfer, genome rearrangement, and the activities of mobile DNA elements. This implies the existence of a delicate balance between the maintenance of genome stability and the tolerance of genome instability. In this review, we describe the specialized genetic elements and the endogenous processes that contribute to genome instability. We then discuss the consequences of genome instability at the physiological level, where cells have harnessed instability to mediate phase and antigenic variation, and at the evolutionary level, where horizontal gene transfer has played an important role. Indeed, this ability to share DNA sequences has played a major part in the evolution of life on Earth. The evolutionary plasticity of bacterial genomes, coupled with the vast numbers of bacteria on the planet, substantially limits our ability to control disease.

Cellular localization and dynamics of the Mrr type IV restriction endonuclease of Escherichia coli

In this study, we examined the intracellular whereabouts of Mrr, a cryptic type IV restriction endonuclease of Escherichia coli K12, in response to different conditions. In absence of stimuli triggering its activity, Mrr was found to be strongly associated with the nucleoid as a number of discrete foci, suggesting the presence of Mrr hotspots on the chromosome. Previously established elicitors of Mrr activity, such as exposure to high (hydrostatic) pressure (HP) or expression of the HhaII methyltransferase, both caused nucleoid condensation and an unexpected coalescence of Mrr foci. However, although the resulting Mrr/nucleoid complex was stable when triggered with HhaII, it tended to be only short-lived when elicited with HP. Moreover, HP-mediated activation of Mrr typically led to cellular blebbing, suggesting a link between chromosome and cellular integrity. Interestingly, Mrr variants could be isolated that were specifically compromised in either HhaII- or HP-dependent activation, underscoring a mechanistic difference in the way both triggers activate Mrr. In general, our results reveal that Mrr can take part in complex spatial distributions on the nucleoid and can be engaged in distinct modes of activity.

Restriction-Modification Systems as a Barrier for Genetic Manipulation of Staphylococcus aureus

Genetic manipulation is a powerful approach to study fundamental aspects of bacterial physiology, metabolism, and pathogenesis. Most Staphylococcus aureus strains are remarkably difficult to genetically manipulate as they possess strong host defense mechanisms that protect bacteria from cellular invasion by foreign DNA. In S. aureus these bacterial "immunity" mechanisms against invading genomes are mainly associated with restriction-modification systems. To date, prokaryotic restriction-modification systems are classified into four different types (Type I-IV), all of which have been found in the sequenced S. aureus genomes. This chapter describes the roles, classification, mechanisms of action of different types of restriction-modification systems and the recent advances in the biology of restriction and modification in S. aureus.

Studying epigenetic DNA modifications in undergraduate laboratories using complementary bioinformatic and molecular approaches

Epigenetic inheritance is the inheritance of genetic information that is not based on DNA sequence alone. One type of epigenetic information that has come to the forefront in the last few years is modified DNA bases. The most common modified DNA base in nature is 5-methylcytosine. Herein, we describe a laboratory experiment that combines bioinformatic and molecular approaches to study the presence and abundance of 5-methylcytosine in different organisms. Students were originally provided with the protein sequence of the Xenopus laevis DNMT1 cytosine-5 DNA methyltransferase and used BLASTP searches to detect the presence of protein orthologs in the genomes of several organisms including Homo sapiens, Mus musculus, Plasmodium falciparum, Drosophila melanogaster, Saccharomyces cerevisiae, Arabidopsis thaliana, and Caenorhabditis elegans. Students generated hypotheses regarding the presence and abundance of 5-methylcytosine in these organisms based on their bioinformatics data, and directly tested their predictions on a subset of DNAs using restriction enzyme isoschizomer assays. A southern blotting assay to answer the same question is also presented. In addition to exposure to the field of epigenetics, the strengths of the laboratory are students are able to make predictions using bioinformatic tools and quickly test them in the laboratory. In addition, students are exposed to two potential misinterpretations of bioinformatic search data. The laboratory is easily modified to incorporate outside research interests in epigenetics.

The other face of restriction: modification-dependent enzymes

The 1952 observation of host-induced non-hereditary variation in bacteriophages by Salvador Luria and Mary Human led to the discovery in the 1960s of modifying enzymes that glucosylate hydroxymethylcytosine in T-even phages and of genes encoding corresponding host activities that restrict non-glucosylated phage DNA: rglA and rglB (restricts glucoseless phage). In the 1980’s, appreciation of the biological scope of these activities was dramatically expanded with the demonstration that plant and animal DNA was also sensitive to restriction in cloning experiments. The rgl genes were renamed mcrA and mcrBC (modified cytosine restriction). The new class of modification-dependent restriction enzymes was named Type IV, as distinct from the familiar modification-blocked Types I–III. A third Escherichia coli enzyme, mrr (modified DNA rejection and restriction) recognizes both methylcytosine and methyladenine. In recent years, the universe of modification-dependent enzymes has expanded greatly. Technical advances allow use of Type IV enzymes to study epigenetic mechanisms in mammals and plants. Type IV enzymes recognize modified DNA with low sequence selectivity and have emerged many times independently during evolution. Here, we review biochemical and structural data on these proteins, the resurgent interest in Type IV enzymes as tools for epigenetic research and the evolutionary pressures on these 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.

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

Comparative genomics of defense systems in archaea and bacteria

Our knowledge of prokaryotic defense systems has vastly expanded as the result of comparative genomic analysis, followed by experimental validation. This expansion is both quantitative, including the discovery of diverse new examples of known types of defense systems, such as restriction-modification or toxin-antitoxin systems, and qualitative, including the discovery of fundamentally new defense mechanisms, such as the CRISPR-Cas immunity system. Large-scale statistical analysis reveals that the distribution of different defense systems in bacterial and archaeal taxa is non-uniform, with four groups of organisms distinguishable with respect to the overall abundance and the balance between specific types of defense systems. The genes encoding defense system components in bacterial and archaea typically cluster in defense islands. In addition to genes encoding known defense systems, these islands contain numerous uncharacterized genes, which are candidates for new types of defense systems. The tight association of the genes encoding immunity systems and dormancy- or cell death-inducing defense systems in prokaryotic genomes suggests that these two major types of defense are functionally coupled, providing for effective protection at the population level.

Diverse functions of restriction-modification systems in addition to cellular defense

Restriction-modification (R-M) systems are ubiquitous and are often considered primitive immune systems in bacteria. Their diversity and prevalence across the prokaryotic kingdom are an indication of their success as a defense mechanism against invading genomes. However, their cellular defense function does not adequately explain the basis for their immaculate specificity in sequence recognition and nonuniform distribution, ranging from none to too many, in diverse species. The present review deals with new developments which provide insights into the roles of these enzymes in other aspects of cellular function. In this review, emphasis is placed on novel hypotheses and various findings that have not yet been dealt with in a critical review. Emerging studies indicate their role in various cellular processes other than host defense, virulence, and even controlling the rate of evolution of the organism. We also discuss how R-M systems could have successfully evolved and be involved in additional cellular portfolios, thereby increasing the relative fitness of their hosts in the population.

An overview of the domestication and impact of the Salmonella mobilome

Salmonella spp. are accountable for a large fraction of the global infectious disease burden, with most of their infections being food- or water-borne. The phenotypic features and adaptive potential of Salmonella spp. appear to be driven to a large extent by mobile or laterally acquired genetic elements. A better understanding of the conduct and diversification of these important pathogens consequently requires a more profound insight into the different mechanisms by which these pivotal elements establish themselves in the cell and affect its behavior. This review, therefore, provides an overview of the physiological impact and domestication of the Salmonella mobilome.

Extensive Inter-Domain Lateral Gene Transfer in the Evolution of the Human Commensal Methanosphaera stadtmanae

Methanosphaera stadtmanae is a commensal methanogenic archaeon found in the human gut. As most of its niche-neighbors are bacteria, it is expected that lateral gene transfer (LGT) from bacteria might have contributed to the evolutionary history of this organism. We performed a phylogenomic survey of putative LGT events in M. stadtmanae, using a phylogenetic pipeline. Our analysis indicates that a substantial fraction of the proteins of M. stadtmanae are inferred to have been involved in inter-domain LGT. Laterally acquired genes have had a large contribution to surface functions, by providing novel glycosyltransferase functions. In addition, several ABC transporters seem to be of bacterial origin, including the molybdate transporter. Thus, bacterial genes contributed to the adaptation of M. stadtmanae to a host-dependent lifestyle by allowing a larger variation in surface structures and increasing transport efficiency in the gut niche which is diverse and competitive.

Movement of DNA sequence recognition domains between non-orthologous proteins

Comparisons of proteins show that they evolve through the movement of domains. However, in many cases, the underlying mechanisms remain unclear. Here, we observed the movements of DNA recognition domains between non-orthologous proteins within a prokaryote genome. Restriction-modification (RM) systems, consisting of a sequence-specific DNA methyltransferase and a restriction enzyme, contribute to maintenance/evolution of genomes/epigenomes. RM systems limit horizontal gene transfer but are themselves mobile. We compared Type III RM systems in Helicobacter pylori genomes and found that target recognition domain (TRD) sequences are mobile, moving between different orthologous groups that occupy unique chromosomal locations. Sequence comparisons suggested that a likely underlying mechanism is movement through homologous recombination of similar DNA sequences that encode amino acid sequence motifs that are conserved among Type III DNA methyltransferases. Consistent with this movement, incongruence was observed between the phylogenetic trees of TRD regions and other regions in proteins. Horizontal acquisition of diverse TRD sequences was suggested by detection of homologs in other Helicobacter species and distantly related bacterial species. One of these RM systems in H. pylori was inactivated by insertion of another RM system that likely transferred from an oral bacterium. TRD movement represents a novel route for diversification of DNA-interacting proteins.

Gene flow and biological conflict systems in the origin and evolution of eukaryotes

The endosymbiotic origin of eukaryotes brought together two disparate genomes in the cell. Additionally, eukaryotic natural history has included other endosymbiotic events, phagotrophic consumption of organisms, and intimate interactions with viruses and endoparasites. These phenomena facilitated large-scale lateral gene transfer and biological conflicts. We synthesize information from nearly two decades of genomics to illustrate how the interplay between lateral gene transfer and biological conflicts has impacted the emergence of new adaptations in eukaryotes. Using apicomplexans as example, we illustrate how lateral transfer from animals has contributed to unique parasite-host interfaces comprised of adhesion- and O-linked glycosylation-related domains. Adaptations, emerging due to intense selection for diversity in the molecular participants in organismal and genomic conflicts, being dispersed by lateral transfer, were subsequently exapted for eukaryote-specific innovations. We illustrate this using examples relating to eukaryotic chromatin, RNAi and RNA-processing systems, signaling pathways, apoptosis and immunity. We highlight the major contributions from catalytic domains of bacterial toxin systems to the origin of signaling enzymes (e.g., ADP-ribosylation and small molecule messenger synthesis), mutagenic enzymes for immune receptor diversification and RNA-processing. Similarly, we discuss contributions of bacterial antibiotic/siderophore synthesis systems and intra-genomic and intra-cellular selfish elements (e.g., restriction-modification, mobile elements and lysogenic phages) in the emergence of chromatin remodeling/modifying enzymes and RNA-based regulation. We develop the concept that biological conflict systems served as evolutionary “nurseries” for innovations in the protein world, which were delivered to eukaryotes via lateral gene flow to spur key evolutionary innovations all the way from nucleogenesis to lineage-specific adaptations.

Polymorphic toxin systems: Comprehensive characterization of trafficking modes, processing, mechanisms of action, immunity and ecology using comparative genomics

BACKGROUND

Proteinaceous toxins are observed across all levels of inter-organismal and intra-genomic conflicts. These include recently discovered prokaryotic polymorphic toxin systems implicated in intra-specific conflicts. They are characterized by a remarkable diversity of C-terminal toxin domains generated by recombination with standalone toxin-coding cassettes. Prior analysis revealed a striking diversity of nuclease and deaminase domains among the toxin modules. We systematically investigated polymorphic toxin systems using comparative genomics, sequence and structure analysis.

RESULTS

Polymorphic toxin systems are distributed across all major bacterial lineages and are delivered by at least eight distinct secretory systems. In addition to type-II, these include type-V, VI, VII (ESX), and the poorly characterized "Photorhabdus virulence cassettes (PVC)", PrsW-dependent and MuF phage-capsid-like systems. We present evidence that trafficking of these toxins is often accompanied by autoproteolytic processing catalyzed by HINT, ZU5, PrsW, caspase-like, papain-like, and a novel metallopeptidase associated with the PVC system. We identified over 150 distinct toxin domains in these systems. These span an extraordinary catalytic spectrum to include 23 distinct clades of peptidases, numerous previously unrecognized versions of nucleases and deaminases, ADP-ribosyltransferases, ADP ribosyl cyclases, RelA/SpoT-like nucleotidyltransferases, glycosyltranferases and other enzymes predicted to modify lipids and carbohydrates, and a pore-forming toxin domain. Several of these toxin domains are shared with host-directed effectors of pathogenic bacteria. Over 90 families of immunity proteins might neutralize anywhere between a single to at least 27 distinct types of toxin domains. In some organisms multiple tandem immunity genes or immunity protein domains are organized into polyimmunity loci or polyimmunity proteins. Gene-neighborhood-analysis of polymorphic toxin systems predicts the presence of novel trafficking-related components, and also the organizational logic that allows toxin diversification through recombination. Domain architecture and protein-length analysis revealed that these toxins might be deployed as secreted factors, through directed injection, or via inter-cellular contact facilitated by filamentous structures formed by RHS/YD, filamentous hemagglutinin and other repeats. Phyletic pattern and life-style analysis indicate that polymorphic toxins and polyimmunity loci participate in cooperative behavior and facultative 'cheating' in several ecosystems such as the human oral cavity and soil. Multiple domains from these systems have also been repeatedly transferred to eukaryotes and their viruses, such as the nucleo-cytoplasmic large DNA viruses.

CONCLUSIONS

Along with a comprehensive inventory of toxins and immunity proteins, we present several testable predictions regarding active sites and catalytic mechanisms of toxins, their processing and trafficking and their role in intra-specific and inter-specific interactions between bacteria. These systems provide insights regarding the emergence of key systems at different points in eukaryotic evolution, such as ADP ribosylation, interaction of myosin VI with cargo proteins, mediation of apoptosis, hyphal heteroincompatibility, hedgehog signaling, arthropod toxins, cell-cell interaction molecules like teneurins and different signaling messengers.

The recognition domain of the methyl-specific endonuclease McrBC flips out 5-methylcytosine

DNA cytosine methylation is a widespread epigenetic mark. Biological effects of DNA methylation are mediated by the proteins that preferentially bind to 5-methylcytosine (5mC) in different sequence contexts. Until now two different structural mechanisms have been established for 5mC recognition in eukaryotes; however, it is still unknown how discrimination of the 5mC modification is achieved in prokaryotes. Here we report the crystal structure of the N-terminal DNA-binding domain (McrB-N) of the methyl-specific endonuclease McrBC from Escherichia coli. The McrB-N protein shows a novel DNA-binding fold adapted for 5mC-recognition. In the McrB-N structure in complex with methylated DNA, the 5mC base is flipped out from the DNA duplex and positioned within a binding pocket. Base flipping elegantly explains why McrBC system restricts only T4-even phages impaired in glycosylation [Luria, S.E. and Human, M.L. (1952) A nonhereditary, host-induced variation of bacterial viruses. J. Bacteriol., 64, 557–569]: flipped out 5-hydroxymethylcytosine is accommodated in the binding pocket but there is no room for the glycosylated base. The mechanism for 5mC recognition employed by McrB-N is highly reminiscent of that for eukaryotic SRA domains, despite the differences in their protein folds.

Conservation of Dcm-mediated cytosine DNA methylation in Escherichia coli

In Escherichia coli, cytosine DNA methylation is catalyzed by the DNA cytosine methyltransferase (Dcm) protein and occurs at the second cytosine in the sequence 5'CCWGG3'. Although the presence of cytosine DNA methylation was reported over 35 years ago, the biological role of 5-methylcytosine in E. coli remains unclear. To gain insight into the role of cytosine DNA methylation in E. coli, we (1) screened the 72 strains of the ECOR collection and 90 recently isolated environmental samples for the presence of the full-length dcm gene using the polymerase chain reaction; (2) examined the same strains for the presence of 5-methylcytosine at 5'CCWGG3' sites using a restriction enzyme isoschizomer digestion assay; and (3) quantified the levels of 5-methyl-2'-deoxycytidine in selected strains using liquid chromatography tandem mass spectrometry. Dcm-mediated cytosine DNA methylation is conserved in all 162 strains examined, and the level of 5-methylcytosine ranges from 0.86% to 1.30% of the cytosines. We also demonstrate that Dcm reduces the expression of ribosomal protein genes during stationary phase, and this may explain the highly conserved nature of this DNA modification pathway.

Functional characterization of the type II PamI restriction-modification system derived from plasmid pAMI7 of Paracoccus aminophilus JCM 7686

Plasmid pAMI7 of the methylotrophic bacterium Paracoccus aminophilus JCM 7686 (Alphaproteobacteria) encodes a functional type II restriction-modification (R-M) system designated PamI. Homologous systems were identified in the genomes of distinct taxonomic groups of Bacteria and Archaea, which provides evidence that horizontal gene transfer has contributed to the wide dissemination of R-M modules - even between domains. Analysis of the cleavage specificity of the R.PamI endonuclease revealed that this protein is an isoschizomer of restriction enzyme NcoI. Interestingly, bioinformatic analyses suggest that R.PamI and NcoI are accompanied by methyltransferases of different methylation specificities (C5-methylcytosine and N4-methylcytosine methyltransferases, respectively), which possibly exemplifies recombinational shuffling of genes coding for individual components of R-M systems. The PamI system can stabilize plasmid pAMI7 in a bacterial population, most probably at the postsegregational level. Therefore, it functions in an analogous manner to plasmid-encoded toxin-antitoxin (TA) systems. Since the TA system of pAMI7 is nonfunctional, it is highly probable that this lack is compensated by the stabilizing activity of PamI. This indicates the crucial role of the analyzed R-M system in the stable maintenance of pAMI7, which is, to our knowledge, the first report of 'symbiosis' between a R-M system and a plasmid in the Alphaproteobacteria.

Defense islands in bacterial and archaeal genomes and prediction of novel defense systems

The arms race between cellular life forms and viruses is a major driving force of evolution. A substantial fraction of bacterial and archaeal genomes is dedicated to antivirus defense. We analyzed the distribution of defense genes and typical mobilome components (such as viral and transposon genes) in bacterial and archaeal genomes and demonstrated statistically significant clustering of antivirus defense systems and mobile genes and elements in genomic islands. The defense islands are enriched in putative operons and contain numerous overrepresented gene families. A detailed sequence analysis of the proteins encoded by genes in these families shows that many of them are diverged variants of known defense system components, whereas others show features, such as characteristic operonic organization, that are suggestive of novel defense systems. Thus, genomic islands provide abundant material for the experimental study of bacterial and archaeal antivirus defense. Except for the CRISPR-Cas systems, different classes of defense systems, in particular toxin-antitoxin and restriction-modification systems, show nonrandom clustering in defense islands. It remains unclear to what extent these associations reflect functional cooperation between different defense systems and to what extent the islands are genomic "sinks" that accumulate diverse nonessential genes, particularly those acquired via horizontal gene transfer. The characteristics of defense islands resemble those of mobilome islands. Defense and mobilome genes are nonrandomly associated in islands, suggesting nonadaptive evolution of the islands via a preferential attachment-like mechanism underpinned by the addictive properties of defense systems such as toxins-antitoxins and an important role of horizontal mobility in the evolution of these islands.

Evolutionary genome engineering using a restriction-modification system

Modification of complex microbial cellular processes is often necessary to obtain organisms with particularly favorable characteristics, but such experiments can take many generations to achieve. In the present article, we accelerated the experimental evolution of Escherichia coli populations under selection for improved growth using one of the restriction-modification systems, which have shaped bacterial genomes. This resulted in faster evolutionary changes in both the genome and bacterial growth. Transcriptome/genome analysis at various stages enabled prompt identification of sequential genome rearrangements and dynamic gene-expression changes associated with growth improvement. The changes were related to cell-to-cell communication, the cell death program, as well as mass production and energy consumption. These observed changes imply that improvements in microorganism population growth can be achieved by inactivating the cellular mechanisms regulating fraction of active cells in a population. Some of the mutations were shown to have additive effects on growth. These results open the way for the application of evolutionary genome engineering to generate organisms with desirable properties.