Transposition effect of adenine (Dam) methylation on activity of O end mutants of IS50

Differential adenine methylation analysis reveals increased variability in 6mA in the absence of methyl-directed mismatch repair

IMPORTANCE

Methylation greatly influences the bacterial genome by guiding DNA repair and regulating pathogenic and stress-response phenotypes. But, the rate of epigenetic changes and their consequences on molecular phenotypes are underexplored. Through a detailed characterization of genome-wide adenine methylation in a commonly used laboratory strain of Escherichia coli, we reveal that mismatch repair deficient populations experience an increase in epimutations resulting in a genome-wide reduction of 6mA methylation in a manner consistent with genetic drift. Our findings highlight how methylation patterns evolve and the constraints on epigenetic evolution due to post-replicative DNA repair, contributing to a deeper understanding of bacterial genome evolution and how epimutations may introduce semi-permanent variation that can influence adaptation.

Bacterial DNA methyltransferase: A key to the epigenetic world with lessons learned from proteobacteria

Epigenetics modulates expression levels of various important genes in both prokaryotes and eukaryotes. These epigenetic traits are heritable without any change in genetic DNA sequences. DNA methylation is a universal mechanism of epigenetic regulation in all kingdoms of life. In bacteria, DNA methylation is the main form of epigenetic regulation and plays important roles in affecting clinically relevant phenotypes, such as virulence, host colonization, sporulation, biofilm formation et al. In this review, we survey bacterial epigenomic studies and focus on the recent developments in the structure, function, and mechanism of several highly conserved bacterial DNA methylases. These methyltransferases are relatively common in bacteria and participate in the regulation of gene expression and chromosomal DNA replication and repair control. Recent advances in sequencing techniques capable of detecting methylation signals have enabled the characterization of genome-wide epigenetic regulation. With their involvement in critical cellular processes, these highly conserved DNA methyltransferases may emerge as promising targets for developing novel epigenetic inhibitors for biomedical applications.

Clostridioides difficile specific DNA adenine methyltransferase CamA squeezes and flips adenine out of DNA helix

Clostridioides difficile infections are an urgent medical problem. The newly discovered C. difficile adenine methyltransferase A (CamA) is specified by all C. difficile genomes sequenced to date (>300), but is rare among other bacteria. CamA is an orphan methyltransferase, unassociated with a restriction endonuclease. CamA-mediated methylation at CAAAAA is required for normal sporulation, biofilm formation, and intestinal colonization by C. difficile. We characterized CamA kinetic parameters, and determined its structure bound to DNA containing the recognition sequence. CamA contains an N-terminal domain for catalyzing methyl transfer, and a C-terminal DNA recognition domain. Major and minor groove DNA contacts in the recognition site involve base-specific hydrogen bonds, van der Waals contacts and the Watson-Crick pairing of a rearranged A:T base pair. These provide sufficient sequence discrimination to ensure high specificity. Finally, the surprisingly weak binding of the methyl donor S-adenosyl-L-methionine (SAM) might provide avenues for inhibiting CamA activity using SAM analogs.

The bacterial epigenome

In all domains of life, genomes contain epigenetic information superimposed over the nucleotide sequence. Epigenetic signals control DNA-protein interactions and can cause phenotypic change in the absence of mutation. A nearly universal mechanism of epigenetic signalling is DNA methylation. In bacteria, DNA methylation has roles in genome defence, chromosome replication and segregation, nucleoid organization, cell cycle control, DNA repair and regulation of transcription. In many bacterial species, DNA methylation controls reversible switching (phase variation) of gene expression, a phenomenon that generates phenotypic cell variants. The formation of epigenetic lineages enables the adaptation of bacterial populations to harsh or changing environments and modulates the interaction of pathogens with their eukaryotic hosts.

Predicting genome terminus sequences of Bacillus cereus-group bacteriophage using next generation sequencing data

Background

Most tailed bacteriophages (phages) feature linear dsDNA genomes. Characterizing novel phages requires an understanding of complete genome sequences, including the definition of genome physical ends.

Result

We sequenced 48 Bacillus cereus phage isolates and analyzed Next-generation sequencing (NGS) data to resolve the genome configuration of these novel phages. Most assembled contigs featured reads that mapped to both contig ends and formed circularized contigs. Independent assemblies of 31 nearly identical I48-like Bacillus phage isolates allowed us to observe that the assembly programs tended to produce random cleavage on circularized contigs. However, currently available assemblers were not capable of reporting the underlying phage genome configuration from sequence data. To identify the genome configuration of sequenced phage in silico, a terminus prediction method was developed by means of ‘neighboring coverage ratios’ and ‘read edge frequencies’ from read alignment files. Termini were confirmed by primer walking and supported by phylogenetic inference of large DNA terminase protein sequences.

Conclusions

The Terminus package using phage NGS data along with the contig circularity could efficiently identify the proximal positions of phage genome terminus. Complete phage genome sequences allow a proposed characterization of the potential packaging mechanisms and more precise genome annotation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3744-0) contains supplementary material, which is available to authorized users.

Characterization of the transposase encoded by IS256, the prototype of a major family of bacterial insertion sequence elements

IS256 is the founding member of the IS256 family of insertion sequence (IS) elements. These elements encode a poorly characterized transposase, which features a conserved DDE catalytic motif and produces circular IS intermediates. Here, we characterized the IS256 transposase as a DNA-binding protein and obtained insight into the subdomain organization and functional properties of this prototype enzyme of IS256 family transposases. Recombinant forms of the transposase were shown to bind specifically to inverted repeats present in the IS256 noncoding regions. A DNA-binding domain was identified in the N-terminal part of the transposase, and a mutagenesis study targeting conserved amino acid residues in this region revealed a putative helix-turn-helix structure as a key element involved in DNA binding. Furthermore, we obtained evidence to suggest that the terminal nucleotides of IS256 are critically involved in IS circularization. Although small deletions at both ends reduced the formation of IS circles, changes at the left-hand IS256 terminus proved to be significantly more detrimental to circle production. Taken together, the data lead us to suggest that the IS256 transposase-mediated circularization reaction preferentially starts with a sequence-specific first-strand cleavage at the left-hand IS terminus.

Bias between the left and right inverted repeats during IS911 targeted insertion

IS911 is a bacterial insertion sequence composed of two consecutive overlapping open reading frames (ORFs [orfA and orfB]) encoding the transposase (OrfAB) as well as a regulatory protein (OrfA). These ORFs are bordered by terminal left and right inverted repeats (IRL and IRR, respectively) with several differences in nucleotide sequence. IS911 transposition is asymmetric: each end is cleaved on one strand to generate a free 3'-OH, which is then used as the nucleophile in attacking the opposite insertion sequence (IS) end to generate a free IS circle. This will be inserted into a new target site. We show here that the ends exhibit functional differences which, in vivo, may favor the use of one compared to the other during transposition. Electromobility shift assays showed that a truncated form of the transposase [OrfAB(1-149)] exhibits higher affinity for IRR than for IRL. While there was no detectable difference in IR activities during the early steps of transposition, IRR was more efficient during the final insertion steps. We show here that the differential activities between the two IRs correlate with the different affinities of OrfAB(1-149) for the IRs during assembly of the nucleoprotein complexes leading to transposition. We conclude that the two inverted repeats are not equivalent during IS911 transposition and that this asymmetry may intervene to determine the ordered assembly of the different protein-DNA complexes involved in the reaction.

Epigenetic gene regulation in the bacterial world

Like many eukaryotes, bacteria make widespread use of postreplicative DNA methylation for the epigenetic control of DNA-protein interactions. Unlike eukaryotes, however, bacteria use DNA adenine methylation (rather than DNA cytosine methylation) as an epigenetic signal. DNA adenine methylation plays roles in the virulence of diverse pathogens of humans and livestock animals, including pathogenic Escherichia coli, Salmonella, Vibrio, Yersinia, Haemophilus, and Brucella. In Alphaproteobacteria, methylation of adenine at GANTC sites by the CcrM methylase regulates the cell cycle and couples gene transcription to DNA replication. In Gammaproteobacteria, adenine methylation at GATC sites by the Dam methylase provides signals for DNA replication, chromosome segregation, mismatch repair, packaging of bacteriophage genomes, transposase activity, and regulation of gene expression. Transcriptional repression by Dam methylation appears to be more common than transcriptional activation. Certain promoters are active only during the hemimethylation interval that follows DNA replication; repression is restored when the newly synthesized DNA strand is methylated. In the E. coli genome, however, methylation of specific GATC sites can be blocked by cognate DNA binding proteins. Blockage of GATC methylation beyond cell division permits transmission of DNA methylation patterns to daughter cells and can give rise to distinct epigenetic states, each propagated by a positive feedback loop. Switching between alternative DNA methylation patterns can split clonal bacterial populations into epigenetic lineages in a manner reminiscent of eukaryotic cell differentiation. Inheritance of self-propagating DNA methylation patterns governs phase variation in the E. coli pap operon, the agn43 gene, and other loci encoding virulence-related cell surface functions.

Molecular genetic analysis of transposase-end DNA sequence recognition: cooperativity of three adjacent base-pairs in specific interaction with a mutant Tn5 transposase

Transposition of Tn5 and IS50 requires the specific binding of transposase (Tnp) to the end inverted repeats, the outside end (OE) and the inside end (IE). OE and IE have 12 identical base-pairs and seven non-identical base-pairs. Previously we described the isolation of a Tnp mutant, EK54, that shows an altered preference for OE versus IE compared to wild-type (wt) Tnp. EK54 enhances OE recognition and decreases IE recognition both in DNA binding and in overall transposition. Here we report that base-pairs 10, 11 and 12 of the OE are critical for the specific recognition by EK54 Tnp. These three adjacent base-pairs act cooperatively; all three must be present in order for EK54 Tnp to interact very favorably with the end DNA. The existence of only one or two of these three base-pairs decreases binding of EK54 Tnp. The combined use of EK54 Tnp and a new OE/IE mosaic end sequence containing the OE base-pairs 10, 11 and 12 gives rise to an extraordinarily high transposition frequency. Just as the Tnp is a multifunctional protein, the nucleotides in the 19 bp Tn5 ends also affect other functions besides Tnp binding. Furthermore, the fact that we were able to isolate end sequence variants that transpose at a higher frequency than the natural ends (OE and IE) with wt Tnp reveals yet another way in which the wt transposition frequency is depressed, i.e. by keeping the end sequences suboptimal.

Tn5 transposase mutants that alter DNA binding specificity

Tn5 transposase (Tnp) binds to Tn5 and IS50 end inverted repeats, the outside end (OE) and the inside end (IE), to initiate transposition. We report the isolation of four Tnp mutants (YH41, TP47, EK54 and EV54) that increase the OE-mediated transposition frequency and enhance the binding affinity of Tnp for OE DNA. In addition, two of the Tnp mutants (TP47 and EK54) appear to be change-of-specificity mutants, since they alter the recognition of OE versus IE relative to the wild-type Tnp. EK54 enhances OE recognition but decreases IE recognition. TP47 enhances both OE and IE recognition but with a much greater enhancement for IE than for OE. This change-of-specificity effect of TP47 is observed only when TP47 Tnp is synthesized in cis to the DNA that contains the ends. We propose that Lys54 makes a favorable interaction with an OE-specific nucleotide pair(s), while Pro47 may cause a more favorable interaction with an IE-specific nucleotide pair(s) than it does with the corresponding OE-specific nucleotide pair(s). A model to explain the preference of TP47 Tnp for the IE in cis but not in trans is proposed.

Tn5 insertion specificity is not influenced by IS50 end sequences in target DNA

The bacterial transposon Tn5 inserts into dozens of sites in a gene, some of which are used preferentially (hotspots). Features of certain sites and precedents provided by several other transposons had suggested that sequences in target DNA corresponding to the ends of Tn5 or of its component IS50 elements might facilitate transposition to these sites. We tested this possibility using derivatives of plasmid pBR322 carrying IS50 I or O end sequences. Tn5 inserted frequently into an IS50 I end at the major hotspot in pBR322, but not into either an I end or an O end 230 bp away from this hotspot. Adenine (dam) methylation at GATC sequences in the I end segment interferes with its use as the end of a transposon, but a dam- mutation did not affect Tn5 insertion relative to an I end sequence in target DNA. These results support models in which the ability of Tn5 to find its preferred sites depends on several features of DNA sequence and conformation, and in which target selection is distinct from recognition of the element ends during transposition.

Ability of DNA and spermidine to affect the activity of restriction endonucleases from several bacterial species

Previous work has described the novel ability to modulate in vitro the activity of restriction endonuclease NaeI from Nocardia aerocoligenes by using cleavable DNA and spermidine [Conrad & Topal (1989) Proc. Natl. Acad. Sci. U.S.A. 86, 9707-9711]. In this paper we report the results of a study of 49 type II restriction enzymes from a variety of bacterial species. On the basis of the rates of cleavage observed, we found that in addition to expected cleavable sites a number of enzymes had slow and resistant cognate recognition sites. Resistant sites were identified for BspMI, NaeI, and NarI; slow sites were identified for HpaII, NaeI, and SacII. Cleavage of these sites was found to be significantly enhanced by the addition of cleavable DNA or spermidine. We demonstrate that for BspMI, as for NaeI, activator DNAs increased Vmax without altering Km, whereas for HpaII, NarI, and SacII activator DNAs decreased Km without changing Vmax. Comparison among the Kms for NaeI cleavage of several different substrates demonstrated that distant DNA sequences can affect DNA recognition by the activated enzyme. Our observations extend DNA activation of the Nocardia NaeI endonuclease to restriction endonucleases from Nocardia argentinensis (NarI), Bacillus species M (BspMI), Haemophilus parainfluenza (HpaII), and Streptomyces achromogenes (SacII). In addition, activation has now been found to affect slow as well as resistant recognition sites.

A sequence at the inside end of IS50 down regulates transposition

Deletion analysis has shown that the segment at the IS50 inside (I) end that is needed for efficient transposition is approximately 19 bp long. Dam methylation at two 5' GATC sequences within this segment decreases I-end transposition activity. A third 5' GATC sequence is present at bp 21-24 of the I end. The comparisons presented here show that extension of the I end from 19 to 24 bp decreases its transposition activity in dam cells 5- to 50-fold, depending on the overall transposon structure.

Intramolecular transposition by a synthetic IS50 (Tn5) derivative

We report the formation of deletions and inversions by intramolecular transposition of Tn5-derived mobile elements. The synthetic transposons used contained the IS50 O and I end segments and the transposase gene, a contraselectable gene encoding sucrose sensitivity (sacB), antibiotic resistance genes, and a plasmid replication origin. Both deletions and inversions were associated with loss of a 300-bp segment that is designated the vector because it is outside of the transposon. Deletions were severalfold more frequent than inversions, perhaps reflecting constraints on DNA twisting or abortive transposition. Restriction and DNA sequence analyses showed that both types of rearrangements extended from one transposon end to many different sites in target DNA. In the case of inversions, transposition generated 9-bp direct repeats of target sequences.

Saturation mutagenesis of the inside end of insertion sequence IS50

A 19-bp segment at the inside (I) end of IS50 (Tn5) is needed for efficient transposition. The importance of each position was assayed by making at least one base substitution at each position by either chemical-or oligodeoxyribonucleotide-directed mutagenesis. Mutant I ends were paired with a wild-type (wt) segment from the outside (O) end of IS50 and the transposase (tnp) gene was placed either between the ends or 1200 bp from the O end. The frequency of transposition of the resultant elements to bacteriophage lambda was measured. At least one substitution at each of the 19 I-end positions decreased transposition activity to less than 25% of wt, and most substitutions (25 of 28) decreased it to less than 5% of wt from one or both donor plasmids. These results show that each position in the I end is important during transposition.