REBASE--a database for DNA restriction and modification: enzymes, genes and genomes

Development of Leptolyngbya sp. BL0902 into a model organism for synthetic biological research in filamentous cyanobacteria

Cyanobacteria have great potential in CO2-based bio-manufacturing and synthetic biological studies. The filamentous cyanobacterium, Leptolyngbya sp. strain BL0902, is comparable to Arthrospira (Spirulina) platensis in commercial-scale cultivation while proving to be more genetically tractable. Here, we report the analyses of the whole genome sequence, gene inactivation/overexpression in the chromosome and deletion of non-essential chromosomal regions in this strain. The genetic manipulations were performed via homologous double recombination using either an antibiotic resistance marker or the CRISPR/Cpf1 editing system for positive selection. A desD-overexpressing strain produced γ-linolenic acid in an open raceway photobioreactor with the productivity of 0.36 g·m-2·d-1. Deletion mutants of predicted patX and hetR, two genes with opposite effects on cell differentiation in heterocyst-forming species, were used to demonstrate an analysis of the relationship between regulatory genes in the non-heterocystous species. Furthermore, a 50.8-kb chromosomal region was successfully deleted in BL0902 with the Cpf1 system. These results supported that BL0902 can be developed into a stable photosynthetic cell factory for synthesizing high value-added products, or used as a model strain for investigating the functions of genes that are unique to filamentous cyanobacteria, and could be systematically modified into a genome-streamlined chassis for synthetic biological purposes.

Generation of tools for expression and purification of the phage-encoded Type I restriction enzyme inhibitor, Ocr

DNA manipulation is an essential tool in molecular microbiology research that is dependent on the ability of bacteria to take up and preserve foreign DNA by horizontal gene transfer. This process can be significantly impaired by the activity of bacterial restriction modification systems; bacterial operons comprising paired enzymatic activities that protectively methylate host DNA, while cleaving incoming unmodified foreign DNA. Ocr is a phage-encoded protein that inhibits Type I restriction modification systems, the addition of which significantly improves bacterial transformation efficiency. We recently established an improved and highly efficient transformation protocol for the important human pathogen group A Streptococcus using commercially available recombinant Ocr protein, manufacture of which has since been discontinued. In order to ensure the continued availability of Ocr protein within the research community, we have generated tools and methods for in-house Ocr production and validated the activity of the purified recombinant protein.

The binding affinity-dependent inhibition of cell growth and viability by DNA sulfur-binding domains

Increasing evidence suggests that DNA phosphorothioate (PT) modification serves several purposes in the bacterial host, and some restriction enzymes specifically target PT-DNA. PT-dependent restriction enzymes (PDREs) bind PT-DNA through their DNA sulfur binding domain (SBD) with dissociation constants (KD) of 5 nM~1 μM. Here, we report that SprMcrA, a PDRE, failed to dissociate from PT-DNA after cleavage due to high binding affinity, resulting in low DNA cleavage efficiency. Expression of SBDs in Escherichia coli cells with PT modification induced a drastic loss of cell viability at 25°C when both DNA strands of a PT site were bound, with one SBD on each DNA strand. However, at this temperature, SBD binding to only one PT DNA strand elicited a severe growth lag rather than lethality. This cell growth inhibition phenotype was alleviated by raising the growth temperature. An in vitro assay mimicking DNA replication and RNA transcription demonstrated that the bound SBD hindered the synthesis of new DNA and RNA when using PT-DNA as the template. Our findings suggest that DNA modification-targeting proteins might regulate cellular processes involved in DNA metabolism in addition to being components of restriction-modification systems and epigenetic readers.

Mini-review of DNA Methylation Detection Techniques and Their Potential Applications in Disease Diagnosis, Prognosis, and Treatment

DNA methylation is the dominant epigenetic mechanism for regulating gene expression in mammals, playing crucial roles in development, differentiation, and tissue homeostasis. Aberrations in DNA methylation are closely associated with the potential onset of various diseases. Consequently, numerous DNA methylation detection techniques have been successively developed. These methods not only facilitate the exploration of disease mechanisms but also hold significant promise for the development of diagnostic and prognostic strategies. In this Perspective, we present a comprehensive overview of commonly employed DNA methylation detection techniques as well as biosensing based on their underlying analytical techniques. For its medical applications, we begin by examining the pathogenesis of different diseases and then proceed to discuss how relevant technologies are applied in the context of these specific medical conditions. Additionally, we briefly discuss the current limitations of these techniques and highlight future challenges in advancing methylation detection and analysis methodologies.

DNA Sequence Control of Enzyme Filamentation and Activation of the SgrAI Endonuclease

Enzyme polymerization (also known as filamentation) has emerged as a new layer of enzyme regulation. SgrAI is a sequence-dependent DNA endonuclease that forms polymeric filaments with enhanced DNA cleavage activity as well as altered DNA sequence specificity. To better understand this unusual regulatory mechanism, full global kinetic modeling of the reaction pathway, including the enzyme filamentation steps, has been undertaken. Prior work with the primary DNA recognition sequence cleaved by SgrAI has shown how the kinetic rate constants of each reaction step are tuned to maximize activation and DNA cleavage while minimizing the extent of DNA cleavage to the host genome. In the current work, we expand on our prior study by now including DNA cleavage of a secondary recognition sequence, to understand how the sequence of the bound DNA modulates filamentation and activation of SgrAI. The work shows that an allosteric equilibrium between low and high activity states is modulated by the sequence of bound DNA, with primary sequences more prone to activation and filament formation, while SgrAI bound to secondary recognition sequences favor the low (and nonfilamenting) state by up to 40-fold. In addition, the degree of methylation of secondary sequences in the host organism, Streptomyces griseus, is now reported for the first time and shows that as predicted, these sequences are left unprotected from the SgrAI endonuclease making sequence specificity critical in this unusual filament-forming enzyme.

An Expanded Genetic Toolbox to Accelerate the Creation of Acholeplasma laidlawii Driven by Synthetic Genomes

We have developed genetic tools for the atypical bacterium Acholeplasma laidlawii. A. laidlawii is a member of the class Mollicutes, which lacks cell walls, has small genomes, and has limited metabolic capabilities, requiring many metabolites from their hosts. Several of these traits have facilitated the development of genome transplantation for some Mollicutes, consequently enabling the generation of synthetic cells. Here, we propose the development of genome transplantation for A. laidlawii. We first investigated a donor-recipient relationship between two strains, PG-8A and PG-8195, through whole-genome sequencing. We then created multihost shuttle plasmids and used them to optimize an electroporation protocol. We also evolved a superior strain for DNA uptake via electroporation. We created a PG-8A donor strain with a Tn5 transposon carrying a tetracycline resistance gene. These tools will enhance Acholeplasma research and accelerate the effort toward creating A. laidlawii strains with synthetic genomes.

Antimicrobial resistance level and conjugation permissiveness shape plasmid distribution in clinical enterobacteria

Conjugative plasmids play a key role in the dissemination of antimicrobial resistance (AMR) genes across bacterial pathogens. AMR plasmids are widespread in clinical settings, but their distribution is not random, and certain associations between plasmids and bacterial clones are particularly successful. For example, the globally spread carbapenem resistance plasmid pOXA-48 can use a wide range of enterobacterial species as hosts, but it is usually associated with a small number of specific Klebsiella pneumoniae clones. These successful associations represent an important threat for hospitalized patients. However, knowledge remains limited about the factors determining AMR plasmid distribution in clinically relevant bacteria. Here, we combined in vitro and in vivo experimental approaches to analyze pOXA-48-associated AMR levels and conjugation dynamics in a collection of wild-type enterobacterial strains isolated from hospitalized patients. Our results revealed significant variability in these traits across different bacterial hosts, with Klebsiella spp. strains showing higher pOXA-48-mediated AMR and conjugation frequencies than Escherichia coli strains. Using experimentally determined parameters, we developed a simple mathematical model to interrogate the contribution of AMR levels and conjugation permissiveness to plasmid distribution in bacterial communities. The simulations revealed that a small subset of clones, combining high AMR levels and conjugation permissiveness, play a critical role in stabilizing the plasmid in different polyclonal microbial communities. These results help to explain the preferential association of plasmid pOXA-48 with K. pneumoniae clones in clinical settings. More generally, our study reveals that species- and strain-specific variability in plasmid-associated phenotypes shape AMR evolution in clinically relevant bacterial communities.

Genomic and transcriptomic insights into complex virus-prokaryote interactions in marine biofilms

Marine biofilms are complex communities of microorganisms that play a crucial ecological role in oceans. Although prokaryotes are the dominant members of these biofilms, little is known about their interactions with viruses. By analysing publicly available and newly sequenced metagenomic data, we identified 2446 virus-prokaryote connections in 84 marine biofilms. Most of these connections were between the bacteriophages in the Uroviricota phylum and the bacteria of Proteobacteria, Cyanobacteria and Bacteroidota. The network of virus-host pairs is complex; a single virus can infect multiple prokaryotic populations or a single prokaryote is susceptible to several viral populations. Analysis of genomes of paired prokaryotes and viruses revealed the presence of 425 putative auxiliary metabolic genes (AMGs), 239 viral genes related to restriction-modification (RM) systems and 38,538 prokaryotic anti-viral defence-related genes involved in 15 defence systems. Transcriptomic evidence from newly established biofilms revealed the expression of viral genes, including AMGs and RM, and prokaryotic defence systems, indicating the active interplay between viruses and prokaryotes. A comparison between biofilms and seawater showed that biofilm prokaryotes have more abundant defence genes than seawater prokaryotes, and the defence gene composition differs between biofilms and the surrounding seawater. Overall, our study unveiled active viruses in natural biofilms and their complex interplay with prokaryotes, which may result in the blooming of defence strategists in biofilms. The detachment of bloomed defence strategists may reduce the infectivity of viruses in seawater and result in the emergence of a novel role of marine biofilms.

Comparison of Yersinia enterocolitica DNA Methylation at Ambient and Host Temperatures

Pathogenic bacteria recognize environmental cues to vary gene expression for host adaptation. Moving from ambient to host temperature, Yersinia enterocolitica responds by immediately repressing flagella synthesis and inducing the virulence plasmid (pYV)-encoded type III secretion system. In contrast, shifting from host to ambient temperature requires 2.5 generations to restore motility, suggesting a link to the cell cycle. We hypothesized that differential DNA methylation contributes to temperature-regulated gene expression. We tested this hypothesis by comparing single-molecule real-time (SMRT) sequencing of Y. enterocolitica DNA from cells growing exponentially at 22 °C and 37 °C. The inter-pulse duration ratio rather than the traditional QV scoring was the kinetic metric to compare DNA from cells grown at each temperature. All 565 YenI restriction sites were fully methylated at both temperatures. Among the 27,118 DNA adenine methylase (Dam) sites, 42 had differential methylation patterns, while 17 remained unmethylated regardless of the temperature. A subset of the differentially methylated Dam sites localized to promoter regions of predicted regulatory genes including LysR-type and PadR-like transcriptional regulators and a cyclic-di-GMP phosphodiesterase. The unmethylated Dam sites localized with a bias to the replication terminus, suggesting they were protected from Dam methylase. No cytosine methylation was detected at Dcm sites.

The DNA cytosine methylome revealed two methylation motifs in the upstream regions of genes related to morphological and physiological differentiation in Streptomyces coelicolor A(3)2 M145

DNA methylation is an epigenetic modification detected in both prokaryotic and eukaryotic genomic DNAs. In bacteria, the importance of 5-methylcytosine (m5C) in gene expression has been less investigated than in eukaryotic systems. Through dot-blot analysis employing m5C antibodies against chromosomal DNA, we have previously demonstrated that m5C influences the differentiation of Streptomyces coelicolor A(3)2 M145 in solid sporulating and liquid non-sporulating complex media. Here, we mapped the methylated cytosines of the M145 strain growing in the defined Maltose Glutamate (MG) liquid medium. Sequencing of the M145 genome after bisulfite treatment (BS-sequencing) evidenced 3360 methylated cytosines and the two methylation motifs, GGCmCGG and GCCmCG, in the upstream regions of 321 genes. Besides, the role of cytosine methylation was investigated using the hypo-methylating agent 5'-aza-2'-deoxycytidine (5-aza-dC) in S. coelicolor cultures, demonstrating that m5C affects both growth and antibiotic biosynthesis. Finally, quantitative reverse-transcription polymerase-chain-reaction (RT-qPCR) analysis of genes containing the methylation motifs in the upstream regions showed that 5-aza-dC treatment influenced their transcriptional levels and those of the regulatory genes for two antibiotics. To the best of our knowledge, this is the first study that reports the cytosine methylome of S. coelicolor M145, supporting the crucial role ascribed to cytosine methylation in controlling bacterial gene expression.

Applications of bioinformatics in epigenetics

Epigenetic modifications such as DNA methylation, post-translational chromatin modifications and non-coding RNA-mediated mechanisms are responsible for epigenetic inheritance. Change in gene expression due to these epigenetic modifications are responsible for new traits in different organisms leading to various diseases including cancer, diabetic kidney disease (DKD), diabetic nephropathy (DN) and renal fibrosis. Bioinformatics is an effective approach for epigenomic profiling. These epigenomic data can be analyzed by a large number of bioinformatics tools and software. Many databases are available online, which comprises huge amount of information regarding these modifications. Recent methodologies include many sequencing and analytical techniques to extrapolate different types of epigenetic data. This data can be used to design drugs against diseases linked to epigenetic modifications. This chapter briefly highlights different epigenetics databases (MethDB, REBASE, Pubmeth, MethPrimerDB, Histone Database, ChromDB, MeInfoText database, EpimiR, Methylome DB, and dbHiMo), and tools (compEpiTools, CpGProD, MethBlAST, EpiExplorer, and BiQ analyzer), which are being utilized to retrieve the data and mechanistically analysis of epigenetics modifications.

Characterization and genomic analysis of a novel Synechococcus phage S-H9-2 belonging to Bristolvirus genus isolated from the Yellow Sea

Cyanophages are known to influence the population dynamics and community structure of cyanobacteria and thus play an important role in biogeochemical cycles in aquatic ecosystems. In this study, a novel Synechococcus phage S-H9-2 infecting Synechococcus sp. WH 8102 was isolated from the coastal water of the Yellow Sea. Synechococcus phage S-H9-2 contains a 187,320 bp genome of double-stranded DNA with a G + C content of 40.3%, 202 potential open reading frames (ORFs), and 15 tRNAs. Phylogenetic analysis and nucleotide-based intergenomic similarity suggest that Synechococcus phage S-H9-2 belongs to the Bristolvirus genus under the family Kyanoviridae. Homologs of the S-H9-2 open reading frame can be found in a variety of marine environments, as shown by the results of mapping the genome sequence of S-H9-2 to the Global Ocean Viromes 2.0 dataset. The presence of auxiliary metabolic genes (AMGs) related to photosynthesis, carbon metabolism, and phosphorus assimilation, as well as phylogenetic relationships based on complete genome sequences, reflect the mechanism of phage-host interaction and host-specific strategies for adaptation to environmental conditions. This study enriches the current genomic database of cyanophage and contributed to our understanding of the virus-host interactions and their adaption to the environment.

Extension of bacterial rDNA sequencing for simultaneous methylation detection and its application in microflora analysis

Although polymerase chain reaction (PCR) amplification and sequencing of the bacterial 16S rDNA region has numerous scientific applications, it does not provide DNA methylation information. Herein, we propose a simple extension for bisulfite sequencing to investigate 5-methylcytosine residues in the bacterial 16S rDNA region from clinical isolates or flora. Multiple displacement amplification without DNA denaturation was used to preferentially pre-amplify single-stranded bacterial DNA after bisulfite conversion. Following the pre-amplification, the 16S rDNA region was analyzed using nested bisulfite PCR and sequencing, enabling the simultaneous identification of DNA methylation status and sequence data. We used this approach (termed sm16S rDNA PCR/sequencing) to identify novel methylation sites and a methyltransferase (M. MmnI) in Morganella morganii and different methylation motifs among Enterococcus faecalis strains from small volumes of clinical specimens. Further, our analysis suggested that M. MmnI may be correlated to erythromycin resistance. Thus, sm16S rDNA PCR/sequencing is a useful extension method for analyzing the DNA methylation of 16S rDNA regions in a microflora, providing additional information not provided by conventional PCR. Given the relationship between DNA methylation status and drug resistance in bacteria, we believe this technique can be effectively applied in clinical sample testing.

Growth condition-dependent differences in methylation imply transiently differentiated DNA methylation states in Escherichia coli

DNA methylation in bacteria frequently serves as a simple immune system, allowing recognition of DNA from foreign sources, such as phages or selfish genetic elements. However, DNA methylation also affects other cell phenotypes in a heritable manner (i.e. epigenetically). While there are several examples of methylation affecting transcription in an epigenetic manner in highly localized contexts, it is not well-established how frequently methylation serves a more general epigenetic function over larger genomic scales. To address this question, here we use Oxford Nanopore sequencing to profile DNA modification marks in three natural isolates of Escherichia coli. We first identify the DNA sequence motifs targeted by the methyltransferases in each strain. We then quantify the frequency of methylation at each of these motifs across the entire genome in different growth conditions. We find that motifs in specific regions of the genome consistently exhibit high or low levels of methylation. Furthermore, we show that there are replicable and consistent differences in methylated regions across different growth conditions. This suggests that during growth, E. coli transiently differentiate into distinct methylation states that depend on the growth state, raising the possibility that measuring DNA methylation alone can be used to infer bacterial growth states without additional information such as transcriptome or proteome data. These results show the utility of using Oxford Nanopore sequencing as an economic means to infer DNA methylation status. They also provide new insights into the dynamics of methylation during bacterial growth and provide evidence of differentiated cell states, a transient analog to what is observed in the differentiation of cell types in multicellular organisms.

DNA Methyltransferase Regulates Nitric Oxide Homeostasis and Virulence in a Chronically Adapted Pseudomonas aeruginosa Strain

Opportunistic pathogens such as Pseudomonas aeruginosa adapt their genomes rapidly during chronic infections. Understanding their epigenetic regulation may provide biomarkers for diagnosis and reveal novel regulatory mechanisms. We performed single-molecule real-time sequencing (SMRT-seq) to characterize the methylome of a chronically adapted P. aeruginosa clinical strain, TBCF10839. Two N⁶-methyladenine (6mA) methylation recognition motifs (RCCANNNNNNNTGAR and TRGANNNNNNTGC [modification sites are in bold]) were identified and predicted as new type I methylation sites using REBASE analysis. We confirmed that the motif TRGANNNNNNTGC was methylated by the methyltransferase (MTase) M.PaeTBCFII, according to methylation sensitivity assays in vivo and vitro. Transcriptomic analysis showed that a ΔpaeTBCFIIM knockout mutant significantly downregulated nitric oxide reductase (NOR) regulation and expression of coding genes such as nosR and norB, which contain methylated motifs in their promoters or coding regions. The ΔpaeTBCFIIM strain exhibited reduced intercellular survival capacity in NO-producing RAW264.7 macrophages and attenuated virulence in a Galleria mellonella infection model; the complemented strain recovered these defective phenotypes. Further phylogenetic analysis demonstrated that homologs of M.PaeTBCFII occur frequently in P. aeruginosa as well as other bacterial species. Our work therefore provided new insights into the relationship between DNA methylation, NO detoxification, and bacterial virulence, laying a foundation for further exploring the molecular mechanism of DNA methyltransferase in regulating the pathogenicity of P. aeruginosa. IMPORTANCE Pseudomonas aeruginosa is an opportunistic pathogen which causes acute and chronic infections that are difficult to treat. Comparative genomic analysis has showed broad genome diversity among P. aeruginosa clinical strains and revealed their different regulatory traits compared to the laboratory strains. While current investigation of the epigenetics of P. aeruginosa is still lacking, understanding epigenetic regulation may provide biomarkers for diagnosis and facilitate development of novel therapies. Denitrification capability is critical for microbial versatility in response to different environmental stress conditions, including the bacterial infection process, where nitric oxide (NO) can be generated by phagocytic cells. The denitrification regulation mechanisms have been studied intensively at genetic and biochemical levels. However, there is very little evidence about the epigenetic regulation of bacterial denitrification mechanism. P. aeruginosa TBCF10839 is a chronically host-adapted strain isolated from a cystic fibrosis (CF) patient with special antiphagocytosis characteristics. Here, we investigated the regulatory effect of an orphan DNA MTase, M.PaeTBCFII, in P. aeruginosa TBCF10839. We demonstrated that the DNA MTase regulates the transcription of denitrification genes represented by NOR and affects antiphagocytic ability in bacteria. In silico analysis suggested that DNA methylation modification may enhance gene expression by affecting the binding of transacting factors such as DNR and RpoN. Our findings not only deepen the understanding of the role of DNA MTase in transcriptional regulation in P. aeruginosa but also provide a theoretical foundation for the in-depth study of the molecular mechanism of the epigenetic regulation on denitrification, virulence, and host-pathogen interaction.

Mining bacterial NGS data vastly expands the complete genomes of temperate phages

Temperate phages (active prophages induced from bacteria) help control pathogenicity, modulate community structure, and maintain gut homeostasis. Complete phage genome sequences are indispensable for understanding phage biology. Traditional plaque techniques are inapplicable to temperate phages due to their lysogenicity, curbing their identification and characterization. Existing bioinformatics tools for prophage prediction usually fail to detect accurate and complete temperate phage genomes. This study proposes a novel computational temperate phage detection method (TemPhD) mining both the integrated active prophages and their spontaneously induced forms (temperate phages) from next-generation sequencing raw data. Applying the method to the available dataset resulted in 192 326 complete temperate phage genomes with different host species, expanding the existing number of complete temperate phage genomes by more than 100-fold. The wet-lab experiments demonstrated that TemPhD can accurately determine the complete genome sequences of the temperate phages, with exact flanking sites, outperforming other state-of-the-art prophage prediction methods. Our analysis indicates that temperate phages are likely to function in the microbial evolution by (i) cross-infecting different bacterial host species; (ii) transferring antibiotic resistance and virulence genes and (iii) interacting with hosts through restriction-modification and CRISPR/anti-CRISPR systems. This work provides a comprehensively complete temperate phage genome database and relevant information, which can serve as a valuable resource for phage research.

Plant DNA Methylation Responds to Nutrient Stress

Nutrient stress as abiotic stress has become one of the important factors restricting crop yield and quality. DNA methylation is an essential epigenetic modification that can effectively regulate genome stability. Exploring DNA methylation responses to nutrient stress could lay the foundation for improving plant tolerance to nutrient stress. This article summarizes the plant DNA methylation patterns, the effects of nutrient stress, such as nitrogen, phosphorus, iron, zinc and sulfur stress, on plant DNA methylation and research techniques for plant DNA methylation, etc. Our discussion provides insight for further research on epigenetics response to nutrient stress in the future.

Intracellular Salmonella Paratyphi A is motile and differs in the expression of flagella-chemotaxis, SPI-1 and carbon utilization pathways in comparison to intracellular S. Typhimurium

Although Salmonella Typhimurium (STM) and Salmonella Paratyphi A (SPA) belong to the same phylogenetic species, share large portions of their genome and express many common virulence factors, they differ vastly in their host specificity, the immune response they elicit, and the clinical manifestations they cause. In this work, we compared their intracellular transcriptomic architecture and cellular phenotypes during human epithelial cell infection. While transcription induction of many metal transport systems, purines, biotin, PhoPQ and SPI-2 regulons was similar in both intracellular SPA and STM, we identified 234 differentially expressed genes that showed distinct expression patterns in intracellular SPA vs. STM. Surprisingly, clear expression differences were found in SPI-1, motility and chemotaxis, and carbon (mainly citrate, galactonate and ethanolamine) utilization pathways, indicating that these pathways are regulated differently during their intracellular phase. Concurring, on the cellular level, we show that while the majority of STM are non-motile and reside within Salmonella-Containing Vacuoles (SCV), a significant proportion of intracellular SPA cells are motile and compartmentalized in the cytosol. Moreover, we found that the elevated expression of SPI-1 and motility genes by intracellular SPA results in increased invasiveness of SPA, following exit from host cells. These findings demonstrate unexpected flagellum-dependent intracellular motility of a typhoidal Salmonella serovar and intriguing differences in intracellular localization between typhoidal and non-typhoidal salmonellae. We propose that these differences facilitate new cycles of host cell infection by SPA and may contribute to the ability of SPA to disseminate beyond the intestinal lamina propria of the human host during enteric fever.

Salmonella enterica is a ubiquitous, facultative intracellular animal and human pathogen. Although non-typhoidal Salmonella (NTS) and typhoidal Salmonella serovars belong to the same phylogenetic species and share many virulence factors, the disease they cause in humans is very different. While the underlying mechanisms for these differences are not fully understood, one possible reason expected to contribute to their different pathogenicity is a distinct expression pattern of genes involved in host-pathogen interactions. Here, we compared the global gene expression and intracellular phenotypes, during human epithelial cell infection of S. Paratyphi A (SPA) and S. Typhimurium (STM), as prototypical serovars of typhoidal and NTS, respectively. Interestingly, we identified different expression patterns in key virulence and metabolic pathways, cytosolic motility and increased reinvasion of SPA, following exit from infected cells. We hypothesize that these differences contribute to the invasive and systemic disease developed following SPA infection in humans.

A high-quality reference genome for the fish pathogen Streptococcus iniae

Fish mortality caused by Streptococcus iniae is a major economic problem in aquaculture in warm and temperate regions globally. There is also risk of zoonotic infection by S. iniae through handling of contaminated fish. In this study, we present the complete genome sequence of S. iniae strain QMA0248, isolated from farmed barramundi in South Australia. The 2.12 Mb genome of S. iniae QMA0248 carries a 32 kb prophage, a 12 kb genomic island and 92 discrete insertion sequence (IS) elements. These include nine novel IS types that belong mostly to the IS3 family. Comparative and phylogenetic analysis between S. iniae QMA0248 and publicly available complete S. iniae genomes revealed discrepancies that are probably due to misassembly in the genomes of isolates ISET0901 and ISNO. Long-range PCR confirmed five rRNA loci in the PacBio assembly of QMA0248, and, unlike S. iniae 89353, no tandemly repeated rRNA loci in the consensus genome. However, we found sequence read evidence that the tandem rRNA repeat existed within a subpopulation of the original QMA0248 culture. Subsequent nanopore sequencing revealed that the tandem rRNA repeat was the most prevalent genotype, suggesting that there is selective pressure to maintain fewer rRNA copies under uncertain laboratory conditions. Our study not only highlights assembly problems in existing genomes, but provides a high-quality reference genome for S. iniae QMA0248, including manually curated mobile genetic elements, that will assist future S. iniae comparative genomic and evolutionary studies.

Genomic Features and Pervasive Negative Selection in Rhodanobacter Strains Isolated from Nitrate and Heavy Metal Contaminated Aquifer

Rhodanobacter species dominate in the Oak Ridge Reservation (ORR) subsurface environments contaminated with acids, nitrate, metal radionuclides, and other heavy metals. To uncover the genomic features underlying adaptations to these mixed-waste environments and to guide genetic tool development, we sequenced the whole genomes of eight Rhodanobacter strains isolated from the ORR site. The genome sizes ranged from 3.9 to 4.2 Mb harboring 3,695 to 4,035 protein-coding genes and GC contents approximately 67%. Seven strains were classified as R. denitrificans and one strain, FW510-R12, as R. thiooxydans based on full length 16S rRNA sequences. According to gene annotation, the top two Cluster of Orthologous Groups (COGs) with high pan-genome expansion rates (Pan/Core gene ratio) were “replication, recombination and repair” and “defense mechanisms.” The denitrifying genes had high DNA homologies except the predicted protein structure variances in NosZ. In contrast, heavy metal resistance genes were diverse with between 7 to 34% of them were located in genomic islands, and these results suggested origins from horizontal gene transfer. Analysis of the methylation patterns in four strains revealed the unique 5mC methylation motifs. Most orthologs (78%) had ratios of nonsynonymous to synonymous substitutions (dN/dS) less than one when compared to the type strain 2APBS1, suggesting the prevalence of negative selection. Overall, the results provide evidence for the important roles of horizontal gene transfer and negative selection in genomic adaptation at the contaminated field site. The complex restriction-modification system genes and the unique methylation motifs in Rhodanobacter strains suggest the potential recalcitrance to genetic manipulation.

IMPORTANCE Despite the dominance of Rhodanobacter species in the subsurface of the contaminated Oak Ridge Reservation (ORR) site, very little is known about the mechanisms underlying their adaptions to the various stressors present at ORR. Recently, multiple Rhodanobacter strains have been isolated from the ORR groundwater samples from several wells with varying geochemical properties. Using Illumina, PacBio, and Oxford Nanopore sequencing platforms, we obtained the whole genome sequences of eight Rhodanobacter strains. Comparison of the whole genomes demonstrated the genetic diversity, and analysis of the long nanopore reads revealed the heterogeneity of methylation patterns in strains isolated from the same well. Although all strains contained a complete set of denitrifying genes, the predicted tertiary structures of NosZ differed. The sequence comparison results demonstrate the important roles of horizontal gene transfer and negative selection in adaptation. In addition, these strains may be recalcitrant to genetic manipulation due to the complex restriction-modification systems and methylations.

Detecting DNA Methylations in the Hyperthermoacidophilic Crenarchaeon Sulfolobus acidocaldarius Using SMRT Sequencing

DNA methylations are one of the most well-known epigenetic modifications along with histone modifications and noncoding RNAs. They are found at specific sites along the DNA in all domains of life, with 5-mC and 6-mA/4-mC being well-characterized in eukaryotes and bacteria respectively, and they have not only been described as contributing to the structure of the double helix itself but also as regulators of DNA-based processes such as replication, transcription, and recombination. Different methods have been developed to accurately identify and/or map methylated motifs to decipher the involvement of DNA methylations in regulatory networks that affect the cellular state.Although DNA methylations have been detected along archaeal genomes, their involvement as regulators of DNA-based processes remains the least known. To highlight the importance of DNA methylations in the control of key cellular mechanisms and their dynamics in archaea cells, we have used single-molecule real-time (SMRT) sequencing. This sequencing technology allows the identification and direct mapping of the methylated motifs along the genome of an organism. In this chapter, we present a step-by-step protocol for detecting DNA methylations in the hyperthermophilic crenarchaeon Sulfolobus acidocaldarius using SMRT sequencing. This protocol can easily be adapted to other prokaryotes.

Characterization of the Type I Restriction Modification System Broadly Conserved among Group A Streptococci

Although prokaryotic DNA methylation investigations have long focused on immunity against exogenous DNA, it has been recently recognized that DNA methylation impacts gene expression and phase variation in Streptococcus pneumoniae and Streptococcus suis. A comprehensive analysis of DNA methylation is lacking for beta-hemolytic streptococci, and thus we sought to examine DNA methylation in the major human pathogen group A Streptococcus (GAS). Using a database of 224 GAS genomes encompassing 80 emm types, we found that nearly all GAS strains encode a type I restriction modification (RM) system that lacks the hsdS′ alleles responsible for impacting gene expression in S. pneumoniae and S. suis. The GAS type I system is located on the core chromosome, while sporadically present type II orphan methyltransferases were identified on prophages. By combining single-molecule real-time (SMRT) analyses of 10 distinct emm types along with phylogenomics of 224 strains, we were able to assign 13 methylation patterns to the GAS population. Inactivation of the type I RM system, occurring either naturally through phage insertion or through laboratory-induced gene deletion, abrogated DNA methylation detectable via either SMRT or MinION sequencing. Contrary to a previous report, inactivation of the type I system did not impact transcript levels of the gene (mga) encoding the key multigene activator protein (Mga) or Mga-regulated genes. Inactivation of the type I system significantly increased plasmid transformation rates. These data delineate the breadth of the core chromosomal type I RM system in the GAS population and clarify its role in immunity rather than impacting Mga regulon expression.

IMPORTANCE The advent of whole-genome approaches capable of detecting DNA methylation has markedly expanded appreciation of the diverse roles of epigenetic modification in prokaryotic physiology. For example, recent studies have suggested that DNA methylation impacts gene expression in some streptococci. The data described herein are from the first systematic analysis of DNA methylation in a beta-hemolytic streptococcus and one of the few analyses to comprehensively characterize DNA methylation across hundreds of strains of the same bacterial species. We clarify that DNA methylation in group A Streptococcus (GAS) is primarily due to a type I restriction modification (RM) system present in the core genome and does not impact mga-regulated virulence gene expression, but does impact immunity against exogenous DNA. The identification of the DNA motifs recognized by each type I RM system may assist with optimizing methods for GAS genetic manipulation and help us understand how bacterial pathogens acquire exogenous DNA elements.

Development of an in vivo methylation system for transformation of Ruminiclostridium cellulolyticum

AIMS

Ruminiclostridium cellulolyticum, an anaerobic cellulolytic bacterium producing an efficient cellulolytic extracellular complex named cellulosome, is a promising host for biofuel production from lignocellulose. This study aims to develop a rapid transformation method for R. cellulolyticum avoiding its restriction system.

METHODS AND RESULTS

The CceI restriction system is a major barrier to introduction of foreign DNA into R. cellulolyticum cells. To improve the transformation efficiency of R. cellulolyticum, the gene encoding CceI methyltransferase (M.CceI) of R. cellulolyticum H10 was functionally expressed in Escherichia coli, resulting in an in vivo methylation system for transformation of R. cellulolyticum. The electrotransformation experiments of R. cellulolyticum H10 with the E. coli-Clostridium shuttle plasmid pMTC6 showed that the transformation efficiency reached up to 2.6 × 10³ ±0.23 × 10³  CFU per μg plasmid DNA. The results demonstrated that the system is able to confer the M.CceI-specific DNA methylation pattern to its resident plasmid, which makes the plasmid resistant to the CceI restriction and efficiently transferred into R. cellulolyticum.

CONCLUSIONS

In this study, we generated an in vivo methylation system of R. cellulolyticum, allowing interspecies DNA transfer and improving transformation efficiency.

SIGNIFICANCE AND IMPACT OF THE STUDY

This research result will greatly facilitate the metabolic engineering of R. cellulolyticum for biofuel production directly from cellulose.

Evolutionary Ecology and Interplay of Prokaryotic Innate and Adaptive Immune Systems

Like many organisms, bacteria and archaea have both innate and adaptive immune systems to defend against infection by viruses and other parasites. Innate immunity most commonly relies on the endonuclease-mediated cleavage of any incoming DNA that lacks a specific epigenetic modification, through a system known as restriction-modification. CRISPR-Cas-mediated adaptive immunity relies on the insertion of short DNA sequences from parasite genomes into CRISPR arrays on the host genome to provide sequence-specific protection. The discovery of each of these systems has revolutionised our ability to carry out genetic manipulations, and, as a consequence, the enzymes involved have been characterised in exquisite detail. In comparison, much less is known about the importance of these two arms of the defence for the ecology and evolution of prokaryotes and their parasites. Here, we review our current ecological and evolutionary understanding of these systems in isolation, and discuss the need to study how innate and adaptive immune responses are integrated when they coexist in the same cell.

Realizing the role of N-acyl-homoserine lactone-mediated quorum sensing in nitrification and denitrification: A review

Nitrification and denitrification are crucial processes in the nitrogen cycle, a vital microbially driven biogeochemical cycle. N-acyl-homoserine lactone (AHL)-mediated quorum sensing (QS) is widespread in bacteria and plays a key role in their physiological status. Recently, there has been an increase in research into how the AHL-mediated QS system is involved in nitrification and denitrification. Consequentially, the AHL-mediated QS system has been considered a promising regulatory approach in nitrogen metabolism processes, with high potential for real-world applications. In this review, the universal presence of QS in nitrifiers and denitrifiers is summarized. Many microorganisms taking part in nitrification and denitrification harbor QS genes, and they may produce AHLs with different chain lengths. The phenotypes and processes affected by QS in real-world applications are also reviewed. In wastewater bioreactors, QS could affect nitrogen metabolism efficiency, granule aggregation, and biofilm formation. Furthermore, methods commonly used to identify the existence and functions of QS, including physiological tests, genetic manipulation and omics analyses are discussed.

Epigenomic Landscape of Lyme Disease Spirochetes Reveals Novel Motifs

Borrelia burgdorferi, the etiological agent of Lyme disease, persists in nature through an enzootic cycle consisting of a vertebrate host and an Ixodes tick vector. The sequence motifs modified by two well-characterized restriction/modification (R/M) loci of B. burgdorferi type strain B31 were recently described, but the methylation profiles of other Lyme disease Borrelia bacteria have not been characterized. Here, the methylomes of B. burgdorferi type strain B31 and 7 clonal derivatives, along with B. burgdorferi N40, B. burgdorferi 297, B. burgdorferi CA-11, B. afzelii PKo, B. afzelii BO23, and B. garinii PBr, were defined through PacBio single-molecule real-time (SMRT) sequencing. This analysis revealed 9 novel sequence motifs methylated by the plasmid-encoded restriction/modification enzymes of these Borrelia strains. Furthermore, while a previous analysis of B. burgdorferi B31 revealed an epigenetic impact of methylation on the global transcriptome, the current data contradict those findings; our analyses of wild-type B. burgdorferi B31 revealed no consistent differences in gene expression among isogenic derivatives lacking one or more restriction/modification enzymes. IMPORTANCE The principal causative agent of Lyme disease in humans in the United States is Borrelia burgdorferi, while B. burgdorferi, B. afzelii, and B. garinii, collectively members of the Borrelia burgdorferi sensu lato species complex, cause Lyme disease in Europe and Asia. Two plasmid-encoded restriction/modification systems have been shown to limit the genetic transformation of B. burgdorferi type strain B31 with foreign DNA, but little is known about the restriction/modification systems of other Lyme disease Borrelia bacteria. This paper describes the methylation motifs present on genomic DNAs of multiple B. burgdorferi, B. afzelii, and B. garinii strains. Contrary to a previous report, we did not find evidence for an epigenetic impact on gene expression by methylation. Knowledge of the motifs recognized and methylated by the restriction/modification enzymes of Lyme disease Borrelia will facilitate molecular genetic investigations of these important human pathogens. Additionally, the similar motifs methylated by orthologous restriction/modification systems of Lyme disease Borrelia bacteria and the presence of these motifs within recombinogenic loci suggest a biological role for these ubiquitous restriction/modification systems in horizontal gene transfer.

The Influence of 5'R and 5'S cdA and cdG on the Activity of BsmAI and SspI Restriction Enzymes

Restriction endonucleases (REs) are intra-bacterial scissors that are considered tools in the fight against foreign genetic material. SspI and BsmAI, examined in this study, cleave dsDNA at their site of recognition or within a short distance of it. Both enzymes are representatives of type II REs, which have played an extremely important role in research on the genetics of organisms and molecular biology. Therefore, the study of agents affecting their activity has become highly important. Ionizing radiation may damage basic cellular mechanisms by inducing lesions in the genome, with 5',8-cyclo-2'-deoxypurines (cdPus) as a model example. Since cdPus may become components of clustered DNA lesions (CDLs), which are unfavorable for DNA repair pathways, their impact on other cellular mechanisms is worthy of attention. This study investigated the influence of cdPus on the elements of the bacterial restriction-modification system. In this study, it was shown that cdPus present in DNA affect the activity of REs. SspI was blocked by any cdPu lesion present at the enzyme's recognition site. When lesions were placed near the recognition sequence, the SspI was inhibited up to 46%. Moreover, (5'S)-5',8-cyclo-2'-deoxyadenosine (ScdA) present in the oligonucleotide sequence lowered BsmAI activity more than (5'R)-5',8-cyclo-2'-deoxyadenosine (RcdA). Interestingly, in the case of 5',8-cyclo-2'-deoxyguanosine (cdG), both 5'S and 5'R diastereomers inhibited BsmAI activity (up to 55% more than cdA). The inhibition was weaker when cdG was present at the recognition site rather than the cleavage site.

Nutrient Loading and Viral Memory Drive Accumulation of Restriction Modification Systems in Bloom-Forming Cyanobacteria

The mechanisms driving cyanobacterial harmful algal blooms (HABs) like those caused by Microcystis aeruginosa remain elusive, but improved defense against viral predation has been implicated for success in eutrophic environments. Our genus-level analyses of 139,023 genomes revealed that HAB-forming cyanobacteria carry vastly more restriction modification systems per genome (RMPG) than nearly all other prokaryotic genera, suggesting that viral defense is a cornerstone of their ecological success. In contrast, picocyanobacteria that numerically dominate nutrient-poor systems have the fewest RMPG within the phylum Cyanobacteria. We used classic resource competition models to explore the hypothesis that nutrient enrichments drive ecological selection for high RMPG due to increased host-phage contact rate. These classic models, agnostic to the mechanism of defense, explain how nutrient loading can select for increased RMPG but, importantly, fail to explain the extreme accumulation of these defense systems. However, extreme accumulation of RMPG can be achieved in a novel “memory” model that accounts for a unique activity of restriction modification systems: the accidental methylation of viral DNA by the methyltransferase. The methylated virus “remembers” the RM defenses of its former host and can evade these defenses if they are present in the next host. This viral memory leads to continual RM system devaluation; RMs accumulate extensively because the benefit of each addition is diminished. Our modeling leads to the hypothesis that nutrient loading and virion methylation drive the extreme accumulation of RMPG in HAB-forming cyanobacteria. Finally, our models suggest that hosts with different RMPG values can coexist when hosts have unique sets of RM systems.

Persistent Bacterial Coinfection of a COVID-19 Patient Caused by a Genetically Adapted Pseudomonas aeruginosa Chronic Colonizer

Pseudomonas aeruginosa is a biofilm-forming opportunistic pathogen which causes chronic infections in immunocompromised patients and leads to high mortality rate. It is identified as a common coinfecting pathogen in COVID-19 patients causing exacerbation of illness. In our hospital, P. aeruginosa is one of the top coinfecting bacteria identified among COVID-19 patients. We collected a strong biofilm-forming P. aeruginosa strain displaying small colony variant morphology from a severe COVID-19 patient. Genomic and transcriptomic sequencing analyses were performed with phenotypic validation to investigate its adaptation in SARS-CoV-2 infected environment. Genomic characterization predicted specific genomic islands highly associated with virulence, transcriptional regulation, and DNA restriction-modification systems. Epigenetic analysis revealed a specific N₆-methyl adenine (m₆A) methylating pattern including methylation of alginate, flagellar and quorum sensing associated genes. Differential gene expression analysis indicated that this isolate formed excessive biofilm by reducing flagellar formation (7.4 to 1,624.1 folds) and overproducing extracellular matrix components including CdrA (4.4 folds), alginate (5.2 to 29.1 folds) and Pel (4.8–5.5 folds). In summary, we demonstrated that P. aeuginosa clinical isolates with novel epigenetic markers could form excessive biofilm, which might enhance its antibiotic resistance and in vivo colonization in COVID-19 patients.

Novel Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV2) and Other Coronaviruses: A Genome-wide Comparative Annotation and Analysis

Novel strain of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV2) causes mild to severe respiratory illness. The early symptoms may be fever, dry cough, sour throat, and difficulty in breathing which may lead to death in severe cases. Compared to previous outbreaks like SARS-CoV and Middle East Respiratory Syndrome (MERS), SARS-CoV2 disease (COVID-19) outbreak has been much distressing due to its high rate of infection but low infection fatality rate (IFR) with 1.4% around the world. World Health Organization (WHO) has declared (COVID-19) a pandemic on March 11, 2020. In the month of January 2020, the whole genome of SARS-CoV2 was sequenced which made work easy for researchers to develop diagnostic kits and to carry out drug repurposing to effectively alleviate the pandemic situation in the world. Now, it is important to understand why this virus has high rate of infectivity or is there any factor involved at the genome level which actually facilitates this virus infection globally? In this study, we have extensively analyzed the whole genomes of different coronaviruses infecting humans and animals in different geographical locations around the world. The main aim of the study is to identify the similarity and the mutational adaptation of the coronaviruses from different host and geographical locations to the SARS-CoV2 and provide a better strategy to understand the mutational rate for specific target-based drug designing. This study is focused to every annotation in a comparative manner which includes SNPs, repeat analysis with the different categorization of the short-sequence repeats and long-sequence repeats, different UTR’s, transcriptional factors, and the predicted matured peptides with the specific length and positions on the genomes. The extensive analysis on SNPs revealed that Wuhan SARS-CoV2 and Indian SARS-CoV2 are having only eight SNPs. Collectively, phylogenetic analysis, repeat analysis, and the polymorphism revealed the genomic conserveness within the SARS-CoV2 and few other coronaviruses with very less mutational chances and the huge distance and mutations from the few other species.

Liquid Chromatography-Mass Spectrometry Analysis of Cytosine Modifications

Liquid chromatography-tandem mass spectrometry (LC-MS/MS) is a widely used technique in the global analysis of epigenetic DNA modifications. The high-resolution chromatographic separation along with sensitive MS detection permits the identification and quantification of deoxyribonucleosides with precision and reliability. Although there have been tremendous advances in LC and MS instrumentation in recent years, sample preparation has not experienced a similar rate of development and is often a bottleneck to chemical analysis. Here we present a protocol for identification and quantification of cytosine modifications that combines a robust and efficient method to generate single nucleosides from genomic DNA samples followed by direct LC-MS/MS analysis.

Multiplex genetic engineering improves endogenous expression of mesophilic α-amylase gene in a wild strain Bacillus amyloliquefaciens 205

A wild strain Bacillus amyloliquefaciens 205 was screened for its high activity of α-amylase. A mesophilic α-amylase encoding gene amyE-205 was revealed and analyzed by genome sequencing. In order to facilitate plasmid transformation to strain 205, an interspecific plasmid transformation method was improved with 5-13 times higher in transformants than that of electronic transformation. A series of CRISPR genome editing tools have been successfully constructed for gene knockout, transcript repression and activation in 205 genome. At this basis, sporulation related genes spo0A and spoIIAC were knockout and suppressed with CRISPR/Cas9 and CRISPR/dCas9 respectively. The double knockout strain 205spo^- was eliminated sporulation with 22.8% increasing of α-amylase activity. The optimal binding site G8 for dCas9-ω has been confirmed in the transcript activation. When amyE-205 was over-expressed with high copy plasmid pUC980-2, its whole upstream sequences containing G8 were also cloned. Whereafter, dCas9-ω was used to activate amyE-205 expression both at genome and plasmid. The final engineered strain 205PG8spo^- achieved 784.3% promotion on α-amylase activity than the starting strain 205. The novel genetic tool box containing an efficient interspecific transformation method and functional CRISPR systems, superadded the multiplex regulation strategies used in strain modification would be also applicative in many Bacillus species.

Methylase-assisted subcloning for high throughput BioBrick assembly

The BioBrick standard makes possible iterated pairwise assembly of cloned parts without any depletion of unique restriction sites. Every part that conforms to the standard is compatible with every other part, thereby fostering a worldwide user community. The assembly methods, however, are labor intensive or inefficient compared to some newer ones so the standard may be falling out of favor. An easier way to assemble BioBricks is described herein. Plasmids encoding BioBrick parts are purified from Escherichia coli cells that express a foreign site-specific DNA methyltransferase, so that each is subsequently protected in vitro from the activity of a particular restriction endonuclease. Each plasmid is double-digested and all resulting restriction fragments are ligated together without gel purification. The ligation products are subsequently double-digested with another pair of restriction endonucleases so only the desired insert-recipient vector construct retains the capacity to transform E. coli. This 4R/2M BioBrick assembly protocol is more efficient and accurate than established workflows including 3A assembly. It is also much easier than gel purification to miniaturize, automate and perform more assembly reactions in parallel. As such, it should streamline DNA assembly for the existing community of BioBrick users, and possibly encourage others to join.

Experimental approaches to tracking mobile genetic elements in microbial communities

Horizontal gene transfer is an important mechanism of microbial evolution and is often driven by the movement of mobile genetic elements between cells. Due to the fact that microbes live within communities, various mechanisms of horizontal gene transfer and types of mobile elements can co-occur. However, the ways in which horizontal gene transfer impacts and is impacted by communities containing diverse mobile elements has been challenging to address. Thus, the field would benefit from incorporating community-level information and novel approaches alongside existing methods. Emerging technologies for tracking mobile elements and assigning them to host organisms provide promise for understanding the web of potential DNA transfers in diverse microbial communities more comprehensively. Compared to existing experimental approaches, chromosome conformation capture and methylome analyses have the potential to simultaneously study various types of mobile elements and their associated hosts. We also briefly discuss how fermented food microbiomes, given their experimental tractability and moderate species complexity, make ideal models to which to apply the techniques discussed herein and how they can be used to address outstanding questions in the field of horizontal gene transfer in microbial communities.

Differential DNA methylomes of clinical MDR, XDR and XXDR Mycobacterium tuberculosis isolates revealed by using single-molecule real-time sequencing

Post-replicative DNA methylation is essential for diverse biological processes in both eukaryotes and prokaryotes. Mycobacterium tuberculosis (M. tuberculosis), the causative agent of tuberculosis, remains one of the most formidable threats worldwide. Although DNA methylation of M. tuberculosis has been documented, little information is available for clinical drug-resistant M. tuberculosis. Single-molecule real-time (SMRT) sequencing was used to profile the core methylome of three clinical isolates, namely multidrug-resistant (MDR), extensively drug-resistant (XDR) and extremely drug-resistant (XXDR) strains. 3812, 6808 and 6041 DNA methylated sites were identified in MDR-MTB, XDR-MTB and XXDR-MTB genome, respectively. There are two types of methylated motifs, namely N⁶-methyladenine (m6A) and N⁴-methylcytosine (m4C). A novel widespread 6 mA methylation motif 5'-CACGCAG-3' was found in XDR-MTB and XXDR-MTB. The methylated genes are involved in multiple cellular processes, especially metabolic enzymes engaged in glucose metabolism, fatty acid and TCA cycle. Many methylated genes are involved in mycobacterial virulence, antibiotic resistance and tolerance. This provided a comprehensive list of methylated genes in drug-resistant clinical isolates and the basis for further functional elucidation.

Avoidance of Self during CRISPR Immunization

The battle between microbes and their viruses is ancient and ongoing. Clustered regularly interspaced short palindromic repeat (CRISPR) immunity, the first and, to date, only form of adaptive immunity found in prokaryotes, represents a flexible mechanism to recall past infections while also adapting to a changing pathogenic environment. Critical to the role of CRISPR as an adaptive immune mechanism is its capacity for self versus non-self recognition when acquiring novel immune memories. Yet, CRISPR systems vary widely in both how and to what degree they can distinguish foreign from self-derived genetic material. We document known and hypothesized mechanisms that bias the acquisition of immune memory towards non-self targets. We demonstrate that diversity is the rule, with many widespread but no universal mechanisms for self versus non-self recognition.

Linking high GC content to the repair of double strand breaks in prokaryotic genomes

Genomic GC content varies widely among microbes for reasons unknown. While mutation bias partially explains this variation, prokaryotes near-universally have a higher GC content than predicted solely by this bias. Debate surrounds the relative importance of the remaining explanations of selection versus biased gene conversion favoring GC alleles. Some environments (e.g. soils) are associated with a high genomic GC content of their inhabitants, which implies that either high GC content is a selective adaptation to particular habitats, or that certain habitats favor increased rates of gene conversion. Here, we report a novel association between the presence of the non-homologous end joining DNA double-strand break repair pathway and GC content; this observation suggests that DNA damage may be a fundamental driver of GC content, leading in part to the many environmental patterns observed to-date. We discuss potential mechanisms accounting for the observed association, and provide preliminary evidence that sites experiencing higher rates of double-strand breaks are under selection for increased GC content relative to the genomic background.

The overall nucleotide composition of an organism’s genome varies greatly between species. Previous work has identified certain environmental factors (e.g., oxygen availability) associated with the relative number of GC bases as opposed to AT bases in the genomes of species. Many of these environments that are associated with high GC content are also associated with relatively high rates of DNA damage. We show that organisms possessing the non-homologous end-joining DNA repair pathway, which is one mechanism to repair DNA double-strand breaks, have an elevated GC content relative to expectation. We also show that certain sites on the genome that are particularly susceptible to double strand breaks have an elevated GC content. This leads us to suggest that an important underlying driver of variability in nucleotide composition across environments is the rate of DNA damage (specifically double-strand breaks) to which an organism living in each environment is exposed.

Functional screening of a soil metagenome for DNA endonucleases by acquired resistance to bacteriophage infection

Endonucleases play a crucial role as reagents in laboratory research and diagnostics. Here, metagenomics was used to functionally screen a fosmid library for endonucleases. A fosmid library was constructed using metagenomic DNA isolated from soil sampled from the unique environment of the Kogelberg Nature Reserve in the Western Cape of South Africa. The principle of acquired immunity against phage infection was used to develop a plate-based screening technique for the isolation of restriction endonucleases from the library. Using next-generation sequencing and bioinformatics tools, sequence data were generated and analysed, revealing 113 novel open reading frames (ORFs) encoding putative endonuclease genes and ORFs of unknown identity and function. One endonuclease designated Endo52 was selected from the putative endonuclease ORFs and was recombinantly produced in Escherichia coli Rosetta™ (DE3) pLysS. Endo52 was purified by immobilised metal affinity chromatography and yielded 0.437 g per litre of cultivation volume. Its enzyme activity was monitored by cleaving lambda DNA and pUC19 plasmid as substrates, and it demonstrated non-specific endonuclease activity. In addition to endonuclease-like genes, the screen identified several unknown genes. These could present new phage resistance mechanisms and are an opportunity for future investigations.

A Stable Genetic Transformation System and Implications of the Type IV Restriction System in the Nitrogen-Fixing Plant Endosymbiont Frankia alni ACN14a

Genus Frankia is comprised primarily of nitrogen-fixing actinobacteria that form root nodule symbioses with a group of hosts known as the actinorhizal plants. These plants are evolutionarily closely related to the legumes that are nodulated by the rhizobia. Both host groups utilize homologs of nodulation genes for root-nodule symbiosis, derived from common plant ancestors. The corresponding endosymbionts, Frankia and the rhizobia, however, are distantly related groups of bacteria, leading to questions about their symbiotic mechanisms and evolutionary history. To date, a stable system of electrotransformation has been lacking in Frankia despite numerous attempts by research groups worldwide. We have identified type IV methyl-directed restriction systems, highly-expressed in a range of actinobacteria, as a likely barrier to Frankia transformation. Here we report the successful electrotransformation of the model strain F. alni ACN14a with an unmethylated, broad host-range replicating plasmid, expressing chloramphenicol-resistance for selection and GFP as a marker of gene expression. This system circumvented the type IV restriction barrier and allowed the stable maintenance of the plasmid. During nitrogen limitation, Frankia differentiates into two cell types: the vegetative hyphae and nitrogen-fixing vesicles. When the expression of egfp under the control of the nif gene cluster promoter was localized using fluorescence imaging, the expression of nitrogen fixation in nitrogen-limited culture was localized in Frankia vesicles but not in hyphae. The ability to separate gene expression patterns between Frankia hyphae and vesicles will enable deeper comparisons of molecular signaling and metabolic exchange between Frankia-actinorhizal and rhizobia-legume symbioses to be made, and may broaden potential applications in agriculture. Further downstream applications are possible, including gene knock-outs and complementation, to open up a range of experiments in Frankia and its symbioses. Additionally, in the transcriptome of F. alni ACN14a, type IV restriction enzymes were highly expressed in nitrogen-replete culture but their expression strongly decreased during symbiosis. The down-regulation of type IV restriction enzymes in symbiosis suggests that horizontal gene transfer may occur more frequently inside the nodule, with possible new implications for the evolution of Frankia.

Genome-wide systematic identification of methyltransferase recognition and modification patterns

Genome-wide analysis of DNA methylation patterns using single molecule real-time DNA sequencing has boosted the number of publicly available methylomes. However, there is a lack of tools coupling methylation patterns and the corresponding methyltransferase genes. Here we demonstrate a high-throughput method for coupling methyltransferases with their respective motifs, using automated cloning and analysing the methyltransferases in vectors carrying a strain-specific cassette containing all potential target sites. To validate the method, we analyse the genomes of the thermophile Moorella thermoacetica and the mesophile Acetobacterium woodii, two acetogenic bacteria having substantially modified genomes with 12 methylation motifs and a total of 23 methyltransferase genes. Using our method, we characterize the 23 methyltransferases, assign motifs to the respective enzymes and verify activity for 11 of the 12 motifs.

Single molecule real-time DNA sequencing allows genome-wide identification of DNA methylation patterns. Here, Jensen et al. present a high-throughput method that allows rapid coupling of DNA methylation patterns with their corresponding methyltransferase genes in bacteria.

SMRT sequencing reveals differential patterns of methylation in two O111:H- STEC isolates from a hemolytic uremic syndrome outbreak in Australia

In 1995 a severe haemolytic-uremic syndrome (HUS) outbreak in Adelaide occurred. A recent genomic analysis of Shiga toxigenic Escherichia coli (STEC) O111:H- strains 95JB1 and 95NR1 from this outbreak found that the more virulent isolate, 95NR1, harboured two additional copies of the Shiga toxin 2 (Stx2) genes encoded within prophage regions. The structure of the Stx2-converting prophages could not be fully resolved using short-read sequence data alone and it was not clear if there were other genomic differences between 95JB1 and 95NR1. In this study we have used Pacific Biosciences (PacBio) single molecule real-time (SMRT) sequencing to characterise the genome and methylome of 95JB1 and 95NR1. We completely resolved the structure of all prophages including two, tandemly inserted, Stx2-converting prophages in 95NR1 that were absent from 95JB1. Furthermore we defined all insertion sequences and found an additional IS1203 element in the chromosome of 95JB1. Our analysis of the methylome of 95NR1 and 95JB1 identified hemi-methylation of a novel motif (5′-CTGC^m6AG-3′) in more than 4000 sites in the 95NR1 genome. These sites were entirely unmethylated in the 95JB1 genome, and included at least 177 potential promoter regions that could contribute to regulatory differences between the strains. IS1203 mediated deactivation of a novel type IIG methyltransferase in 95JB1 is the likely cause of the observed differential patterns of methylation between 95NR1 and 95JB1. This study demonstrates the capability of PacBio SMRT sequencing to resolve complex prophage regions and reveal the genetic and epigenetic heterogeneity within a clonal population of bacteria.

Visualization and prediction of CRISPR incidence in microbial trait-space to identify drivers of antiviral immune strategy

Bacteria and archaea are locked in a near-constant battle with their viral pathogens. Despite previous mechanistic characterization of numerous prokaryotic defense strategies, the underlying ecological drivers of different strategies remain largely unknown and predicting which species will take which strategies remains a challenge. Here, we focus on the CRISPR immune strategy and develop a phylogenetically-corrected machine learning approach to build a predictive model of CRISPR incidence using data on over 100 traits across over 2600 species. We discover a strong but hitherto-unknown negative interaction between CRISPR and aerobicity, which we hypothesize may result from interference between CRISPR-associated proteins and non-homologous end-joining DNA repair due to oxidative stress. Our predictive model also quantitatively confirms previous observations of an association between CRISPR and temperature. Finally, we contrast the environmental associations of different CRISPR system types (I, II, III) and restriction modification systems, all of which act as intracellular immune systems.

DNA phosphorothioate modification-a new multi-functional epigenetic system in bacteria

Synthetic phosphorothioate (PT) internucleotide linkages, in which a nonbridging oxygen is replaced by a sulphur atom, share similar physical and chemical properties with phosphodiesters but confer enhanced nuclease tolerance on DNA/RNA, making PTs a valuable biochemical and pharmacological tool. Interestingly, PT modification was recently found to occur naturally in bacteria in a sequence-selective and RP configuration-specific manner. This oxygen-sulphur swap is catalysed by the gene products of dndABCDE, which constitute a defence barrier with DndFGH in some bacterial strains that can distinguish and attack non-PT-modified foreign DNA, resembling DNA methylation-based restriction-modification (R-M) systems. Despite their similar defensive mechanisms, PT- and methylation-based R-M systems have evolved to target different consensus contexts in the host cell because when they share the same recognition sequences, the protective function of each can be impeded. The redox and nucleophilic properties of PT sulphur render PT modification a versatile player in the maintenance of cellular redox homeostasis, epigenetic regulation and environmental fitness. The widespread presence of dnd systems is considered a consequence of extensive horizontal gene transfer, whereas the lability of PT during oxidative stress and the susceptibility of PT to PT-dependent endonucleases provide possible explanations for the ubiquitous but sporadic distribution of PT modification in the bacterial world.

A new type of DNA phosphorothioation-based antiviral system in archaea

Archaea and Bacteria have evolved different defence strategies that target virtually all steps of the viral life cycle. The diversified virion morphotypes and genome contents of archaeal viruses result in a highly complex array of archaea-virus interactions. However, our understanding of archaeal antiviral activities lags far behind our knowledges of those in bacteria. Here we report a new archaeal defence system that involves DndCDEA-specific DNA phosphorothioate (PT) modification and the PbeABCD-mediated halt of virus propagation via inhibition of DNA replication. In contrast to the breakage of invasive DNA by DndFGH in bacteria, DndCDEA-PbeABCD does not degrade or cleave viral DNA. The PbeABCD-mediated PT defence system is widespread and exhibits extensive interdomain and intradomain gene transfer events. Our results suggest that DndCDEA-PbeABCD is a new type of PT-based virus resistance system, expanding the known arsenal of defence systems as well as our understanding of host-virus interactions.

Adaptive dif Modules in Permafrost Strains of Acinetobacter lwoffii and Their Distribution and Abundance Among Present Day Acinetobacter Strains

The dif/Xer system of site-specific recombination allows resolution of chromosomal dimers during bacterial DNA replication. Recently, it was also shown to be involved in horizontal transfer of a few known Xer-dependent mobile elements. Here, we show that plasmids of various Acinetobacter species, including clinically important strains, often contain multiple pdif sites that are mainly located within their accessory regions. Chromosomes of Acinetobacter strains may also contain additional dif sites, and their similarity with plasmid pdif sites is higher than with the main chromosomal site dif1. We further identify putative mobile genetic elements containing pdif sites on both flanks of adaptive genes and analyze their distribution in Acinetobacter species. In total, we describe seven mobile elements containing genes with various adaptive functions from permafrost strains of A. lwoffii group. All of them are also spread in modern plasmids of different Acinetobacter species including A. baumannii. We could not detect pdif sites and corresponding mobile elements in closely related bacterial genera, including Psychrobacter and Moraxella. Thus, the widespread distribution of dif modules is a characteristic feature of Acinetobacter species and may contribute to their high adaptability both in the environment and in the clinic.

The Need for Speed: Run-On Oligomer Filament Formation Provides Maximum Speed with Maximum Sequestration of Activity

Here, we investigate an unusual antiviral mechanism developed in the bacterium Streptomyces griseus SgrAI is a type II restriction endonuclease that forms run-on oligomer filaments when activated and possesses both accelerated DNA cleavage activity and expanded DNA sequence specificity. Mutations disrupting the run-on oligomer filament eliminate the robust antiphage activity of wild-type SgrAI, and the observation that even relatively modest disruptions completely abolish this anti-viral activity shows that the greater speed imparted by the run-on oligomer filament mechanism is critical to its biological function. Simulations of DNA cleavage by SgrAI uncover the origins of the kinetic advantage of this newly described mechanism of enzyme regulation over more conventional mechanisms, as well as the origin of the sequestering effect responsible for the protection of the host genome against damaging DNA cleavage activity of activated SgrAI.IMPORTANCE This work is motivated by an interest in understanding the characteristics and advantages of a relatively newly discovered enzyme mechanism involving filament formation. SgrAI is an enzyme responsible for protecting against viral infections in its host bacterium and was one of the first such enzymes shown to utilize such a mechanism. In this work, filament formation by SgrAI is disrupted, and the effects on the speed of the purified enzyme as well as its function in cells are measured. It was found that even small disruptions, which weaken but do not destroy filament formation, eliminate the ability of SgrAI to protect cells from viral infection, its normal biological function. Simulations of enzyme activity were also performed and show how filament formation can greatly speed up an enzyme's activation compared to that of other known mechanisms, as well as to better localize its action to molecules of interest, such as invading phage DNA.