Enhanced 5-methylcytosine detection in single-molecule, real-time sequencing via Tet1 oxidation

Fluorimetric detection of DNA methylation by cerium oxide nanoparticles for early cancer diagnosis

In this study, a very sensitive fluorescence nano-biosensor was developed using CeO2 nanoparticles for the rapid detection of DNA methylation. The characteristics of CeO2 nanoparticles were determined by transmission electron microscopy (TEM), scanning electron microscopy (SEM), energy dispersive spectroscopy (EDS), X-ray diffraction (XRD) spectroscopy, UV-visible spectroscopy, and fluorescence spectroscopy. The CeO2 nanoparticles were reacted with a single-stranded DNA (ssDNA) probe, and then methylated and unmethylated target DNAs hybridized with an ssDNA probe, and the fluorescence emission was measured. Upon adding the target unmethylated and methylated ssDNA, the fluorescence intensity increased in the linear range of concentration from 2 × 10-13 - 10-18 M. The limit of detection (LOD) was 1.597 × 10-6 M for methylated DNA and 1.043 × 10-6 M for unmethylated DNA. The fluorescence emission intensity of methylated sequences was higher than of that unmethylated sequences. The fabricated DNA nanobiosensor showed a fluorescence emission at 420 nm with an excitation wavelength of 280 nm. The impact of CeO2 binding on methylated and unmethylated DNA was further demonstrated by agarose gel electrophoresis. Finally, the actual sample analysis suggested that the nanobiosensor could have practical applications for detecting methylation in the human plasma samples.

5-Hydroxymethylcytosine: the many faces of the sixth base of mammalian DNA

Epigenetic phenomena play a central role in cell regulatory processes and are important factors for understanding complex human disease. One of the best understood epigenetic mechanisms is DNA methylation. In the mammalian genome, cytosines (C) in CpG dinucleotides were long known to undergo methylation at the 5-position of the pyrimidine ring (mC). Later it was found that mC can be oxidized to 5-hydroxymethylcytosine (hmC) or even further to 5-formylcytosine (fC) and to 5-carboxylcytosine (caC) by the action of 2-oxoglutarate-dependent dioxygenases of the TET family. These findings unveiled a long elusive mechanism of active DNA demethylation and bolstered a wave of studies in the area of epigenetic regulation in mammals. This review is dedicated to critical assessment of recent data on biochemical and chemical aspects of the formation and conversion of hmC in DNA, analytical techniques used for detection and mapping of this nucleobase in mammalian genomes as well as epigenetic roles of hmC in DNA replication, transcription, cell differentiation and human disease.

Detecting DNA hydroxymethylation: exploring its role in genome regulation

DNA methylation is one of the most extensively studied epigenetic regulatory mechanisms, known to play crucial roles in various organisms. It has been implicated in the regulation of gene expression and chromatin changes, ranging from global alterations during cell state transitions to locus-specific modifications. 5-hydroxymethylcytosine (5hmC) is produced by a major oxidation, from 5-methylcytosine (5mC), catalyzed by the ten-eleven translocation (TET) enzymes, and is gradually being recognized for its significant role in genome regulation. With the development of state-of-the-art experimental techniques, it has become possible to detect and distinguish 5mC and 5hmC at base resolution. Various techniques have evolved, encompassing chemical and enzymatic approaches, as well as thirdgeneration sequencing techniques. These advancements have paved the way for a thorough exploration of the role of 5hmC across a diverse array of cell types, from embryonic stem cells (ESCs) to various differentiated cells. This review aims to comprehensively report on recent techniques and discuss the emerging roles of 5hmC. [BMB Reports 2024; 57(3): 135-142].

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.

Chemical Methods to Identify Epigenetic Modifications in Cytosine Bases

Chemical modifications to Cytosine bases are among the most studied epigenetic markers and their detection in the human genome plays a crucial role in gaining more insights about gene regulation, prognosis of genetic disorders and unraveling genetic inheritance patterns. The Cytosine methylated at the 5th position and oxidized derivatives thereof generated in the demethylation pathways, perform separate and unique epigenetic functions in an organism. As the presence of various Cytosine modifications is associated with diverse diseases, including cancer, there has been a strong focus on developing methods, both chemical and alternative approaches, capable of detecting these modifications at a single-base resolution across the entire genome. In this comprehensive review, we aim to consolidate the various chemical methods and understanding their chemistry that have been established to date for the detection of various Cytosine modifications.

The Helicobacter pylori methylome is acid-responsive due to regulation by the two-component system ArsRS and the type I DNA methyltransferase HsdM1 (HP0463)

In addition to its role in genome protection, DNA methylation can regulate gene expression. In this study, we characterized the impact of acidity, phase variation, and the ArsRS TCS on the expression of the Type I m6A DNA methyltransferase HsdM1 (HP0463) of Helicobacter pylori 26695 and their subsequent effects on the methylome. Transcription of hsdM1 increases at least fourfold in the absence of the sensory histidine kinase ArsS, the major acid-sensing protein of H. pylori. hsdM1 exists in the phase-variable operon hsdR1-hsdM1. Phase-locking hsdR1 (HP0464), the restriction endonuclease gene, has significant impacts on the transcription of hsdM1. To determine the impacts of methyltransferase transcription patterns on the methylome, we conducted methylome sequencing on samples cultured at pH 7 or pH 5. We found differentially methylated motifs between these growth conditions and that deletions of arsS and/or hsdM1 interfere with the epigenetic acid response. Deletion of arsS leads to altered activity of HsdM1 and multiple other methyltransferases under both pH conditions indicating that the ArsRS TCS, in addition to direct effects on regulon transcription during acid acclimation, may also indirectly impact gene expression via regulation of the methylome. We determined the target motif of HsdM1 (HP0463) to be the complementary bipartite sequence pair 5'-TCAm6AVN6TGY-3' and 3'-AGTN6GAm6ACA-5'. This complex regulation of DNA methyltransferases, and thus differential methylation patterns, may have implications for the decades-long persistent infection by H. pylori. IMPORTANCE This study expands the possibilities for complex, epigenomic regulation in Helicobacter pylori. We demonstrate that the H. pylori methylome is plastic and acid sensitive via the two-component system ArsRS and the DNA methyltransferase HsdM1. The control of a methyltransferase by ArsRS may allow for a layered response to changing acidity. Likely, an early response whereby ArsR~P affects regulon expression, including the methyltransferase hsdM1. Then, a somewhat later effect as the altered methylome, due to altered HsdM1 expression, subsequently alters the expression of other genes involved in acclimation. The intermediate methylation of certain motifs supports the hypothesis that methyltransferases play a regulatory role. Untangling this additional web of regulation could play a key role in understanding H. pylori colonization and persistence.

Genomics in the long-read sequencing era

Long-read sequencing (LRS) technologies have provided extremely powerful tools to explore genomes. While in the early years these methods suffered technical limitations, they have recently made significant progress in terms of read length, throughput, and accuracy and bioinformatics tools have strongly improved. Here, we aim to review the current status of LRS technologies, the development of novel methods, and the impact on genomics research. We will explore the most impactful recent findings made possible by these technologies focusing on high-resolution sequencing of genomes and transcriptomes and the direct detection of DNA and RNA modifications. We will also discuss how LRS methods promise a more comprehensive understanding of human genetic variation, transcriptomics, and epigenetics for the coming years.

Navigating the pitfalls of mapping DNA and RNA modifications

Chemical modifications to nucleic acids occur across the kingdoms of life and carry important regulatory information. Reliable high-resolution mapping of these modifications is the foundation of functional and mechanistic studies, and recent methodological advances based on next-generation sequencing and long-read sequencing platforms are critical to achieving this aim. However, mapping technologies may have limitations that sometimes lead to inconsistent results. Some of these limitations are technical in nature and specific to certain types of technology. Here, however, we focus on common (yet not always widely recognized) pitfalls that are shared among frequently used mapping technologies and discuss strategies to help technology developers and users mitigate their effects. Although the emphasis is primarily on DNA modifications, RNA modifications are also discussed.

Third-Generation Sequencing of Epigenetic DNA

The discovery of epigenetic bases has revolutionised the understanding of disease and development. Among the most studied epigenetic marks are cytosines covalently modified at the 5 position. In order to gain insight into their biological significance, the ability to determine their spatiotemporal distribution within the genome is essential. Techniques for sequencing on "next-generation" platforms often involve harsh chemical treatments leading to sample degradation. Third-generation sequencing promises to further revolutionise the field by providing long reads, enabling coverage of highly repetitive regions of the genome or structural variants considered unmappable by next generation sequencing technology. While the ability of third-generation platforms to directly detect epigenetic modifications is continuously improving, at present chemical or enzymatic derivatisation presents the most convenient means of enhancing reliability. This Review presents techniques available for the detection of cytosine modifications on third-generation platforms.

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

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

Methyl-SNP-seq reveals dual readouts of methylome and variome at molecule resolution while enabling target enrichment

Covalent modifications of genomic DNA are crucial for most organisms to survive. Amplicon-based high-throughput sequencing technologies erase all DNA modifications to retain only sequence information for the four canonical nucleobases, necessitating specialized technologies for ascertaining epigenetic information. To also capture base modification information, we developed Methyl-SNP-seq, a technology that takes advantage of the complementarity of the double helix to extract the methylation and original sequence information from a single DNA molecule. More specifically, Methyl-SNP-seq uses bisulfite conversion of one of the strands to identify cytosine methylation while retaining the original four-bases sequence information on the other strand. As both strands are locked together to link the dual readouts on a single paired-end read, Methyl-SNP-seq allows detecting the methylation status of any DNA even without a reference genome. Because one of the strands retains the original four nucleotide composition, Methyl-SNP-seq can also be used in conjunction with standard sequence-specific probes for targeted enrichment and amplification. We show the usefulness of this technology in a broad spectrum of applications ranging from allele-specific methylation analysis in humans to identification of methyltransferase specificity in complex bacterial communities.

Rapid identification of methylase specificity (RIMS-seq) jointly identifies methylated motifs and generates shotgun sequencing of bacterial genomes

DNA methylation is widespread amongst eukaryotes and prokaryotes to modulate gene expression and confer viral resistance. 5-Methylcytosine (m5C) methylation has been described in genomes of a large fraction of bacterial species as part of restriction-modification systems, each composed of a methyltransferase and cognate restriction enzyme. Methylases are site-specific and target sequences vary across organisms. High-throughput methods, such as bisulfite-sequencing can identify m5C at base resolution but require specialized library preparations and single molecule, real-time (SMRT) sequencing usually misses m5C. Here, we present a new method called RIMS-seq (rapid identification of methylase specificity) to simultaneously sequence bacterial genomes and determine m5C methylase specificities using a simple experimental protocol that closely resembles the DNA-seq protocol for Illumina. Importantly, the resulting sequencing quality is identical to DNA-seq, enabling RIMS-seq to substitute standard sequencing of bacterial genomes. Applied to bacteria and synthetic mixed communities, RIMS-seq reveals new methylase specificities, supporting routine study of m5C methylation while sequencing new genomes.

Quantification and mapping of DNA modifications

Apart from the four canonical nucleobases, DNA molecules carry a number of natural modifications. Substantial evidence shows that DNA modifications can regulate diverse biological processes. Dynamic and reversible modifications of DNA are critical for cell differentiation and development. Dysregulation of DNA modifications is closely related to many human diseases. The research of DNA modifications is a rapidly expanding area and has been significantly stimulated by the innovations of analytical methods. With the recent advances in methods and techniques, a series of new DNA modifications have been discovered in the genomes of prokaryotes and eukaryotes. Deciphering the biological roles of DNA modifications depends on the sensitive detection, accurate quantification, and genome-wide mapping of modifications in genomic DNA. This review provides an overview of the recent advances in analytical methods and techniques for both the quantification and genome-wide mapping of natural DNA modifications. We discuss the principles, advantages, and limitations of these developed methods. It is anticipated that new methods and techniques will resolve the current challenges in this burgeoning research field and expedite the elucidation of the functions of DNA modifications.

Comparison of long-read sequencing technologies in interrogating bacteria and fly genomes

The newest generation of DNA sequencing technology is highlighted by the ability to generate sequence reads hundreds of kilobases in length. Pacific Biosciences (PacBio) and Oxford Nanopore Technologies (ONT) have pioneered competitive long read platforms, with more recent work focused on improving sequencing throughput and per-base accuracy. We used whole-genome sequencing data produced by three PacBio protocols (Sequel II CLR, Sequel II HiFi, RS II) and two ONT protocols (Rapid Sequencing and Ligation Sequencing) to compare assemblies of the bacteria Escherichia coli and the fruit fly Drosophila ananassae. In both organisms tested, Sequel II assemblies had the highest consensus accuracy, even after accounting for differences in sequencing throughput. ONT and PacBio CLR had the longest reads sequenced compared to PacBio RS II and HiFi, and genome contiguity was highest when assembling these datasets. ONT Rapid Sequencing libraries had the fewest chimeric reads in addition to superior quantification of E. coli plasmids versus ligation-based libraries. The quality of assemblies can be enhanced by adopting hybrid approaches using Illumina libraries for bacterial genome assembly or polishing eukaryotic genome assemblies, and an ONT-Illumina hybrid approach would be more cost-effective for many users. Genome-wide DNA methylation could be detected using both technologies, however ONT libraries enabled the identification of a broader range of known E. coli methyltransferase recognition motifs in addition to undocumented D. ananassae motifs. The ideal choice of long read technology may depend on several factors including the question or hypothesis under examination. No single technology outperformed others in all metrics examined.

Genome-wide detection of cytosine methylation by single molecule real-time sequencing

Single molecule real-time (SMRT) sequencing theoretically offers the opportunity to directly assess certain base modifications of native DNA molecules without any prior chemical/enzymatic conversions and PCR amplification, using kinetic signals of a DNA polymerase. However, the kinetic signal changes caused by 5mC modification are extremely subtle. Hence, the robust genome-wide measurement of 5mC modification has not been achieved. We enhanced 5mC detection using SMRT sequencing by holistically analyzing kinetic signals of a DNA polymerase and sequence context for every base within a measurement window. We employed a convolutional neural network to train a methylation classification model, leading to genome-wide 5mC detection. The sensitivity and specificity reached 90% and 94%, with a 99% correlation of overall methylation level with bisulfite sequencing.

5-Methylcytosine (5mC) is an important type of epigenetic modification. Bisulfite sequencing (BS-seq) has limitations, such as severe DNA degradation. Using single molecule real-time sequencing, we developed a methodology to directly examine 5mC. This approach holistically examined kinetic signals of a DNA polymerase (including interpulse duration and pulse width) and sequence context for every nucleotide within a measurement window, termed the holistic kinetic (HK) model. The measurement window of each analyzed double-stranded DNA molecule comprised 21 nucleotides with a cytosine in a CpG site in the center. We used amplified DNA (unmethylated) and M.SssI-treated DNA (methylated) (M.SssI being a CpG methyltransferase) to train a convolutional neural network. The area under the curve for differentiating methylation states using such samples was up to 0.97. The sensitivity and specificity for genome-wide 5mC detection at single-base resolution reached 90% and 94%, respectively. The HK model was then tested on human–mouse hybrid fragments in which each member of the hybrid had a different methylation status. The model was also tested on human genomic DNA molecules extracted from various biological samples, such as buffy coat, placental, and tumoral tissues. The overall methylation levels deduced by the HK model were well correlated with those by BS-seq (r = 0.99; P < 0.0001) and allowed the measurement of allele-specific methylation patterns in imprinted genes. Taken together, this methodology has provided a system for simultaneous genome-wide genetic and epigenetic analyses.

Genome-wide identification of 5-methylcytosine sites in bacterial genomes by high-throughput sequencing of MspJI restriction fragments

Single-molecule Real-Time (SMRT) sequencing can easily identify sites of N6-methyladenine and N4-methylcytosine within DNA sequences, but similar identification of 5-methylcytosine sites is not as straightforward. In prokaryotic DNA, methylation typically occurs within specific sequence contexts, or motifs, that are a property of the methyltransferases that "write" these epigenetic marks. We present here a straightforward, cost-effective alternative to both SMRT and bisulfite sequencing for the determination of prokaryotic 5-methylcytosine methylation motifs. The method, called MFRE-Seq, relies on excision and isolation of fully methylated fragments of predictable size using MspJI-Family Restriction Enzymes (MFREs), which depend on the presence of 5-methylcytosine for cleavage. We demonstrate that MFRE-Seq is compatible with both Illumina and Ion Torrent sequencing platforms and requires only a digestion step and simple column purification of size-selected digest fragments prior to standard library preparation procedures. We applied MFRE-Seq to numerous bacterial and archaeal genomic DNA preparations and successfully confirmed known motifs and identified novel ones. This method should be a useful complement to existing methodologies for studying prokaryotic methylomes and characterizing the contributing methyltransferases.

Methylation analysis of Klebsiella pneumoniae from Portuguese hospitals

Klebsiella pneumoniae is an important nosocomial infectious agent with a high antimicrobial resistance (AMR) burden. The application of long read sequencing technologies is providing insights into bacterial chromosomal and putative extra-chromosomal genetic elements (PEGEs) associated with AMR, but also epigenetic DNA methylation, which is thought to play a role in cleavage of foreign DNA and expression regulation. Here, we apply the PacBio sequencing platform to eight Portuguese hospital isolates, including one carbapenemase producing isolate, to identify methylation motifs. The resulting assembled chromosomes were between 5.2 and 5.5Mbp in length, and twenty-six PEGEs were found. Four of our eight samples carry blaCTX-M-15, a dominant Extended Spectrum Beta Lactamase in Europe. We identified methylation motifs that control Restriction-Modification systems, including GATC of the DNA adenine methylase (Dam), which methylates N6-methyladenine (m6A) across all our K. pneumoniae assemblies. There was a consistent lack of methylation by Dam of the GATC motif downstream of two genes: fosA, a locus associated with low level fosfomycin resistance, and tnpB transposase on IncFIB(K) plasmids. Overall, we have constructed eight high quality reference genomes of K. pneumoniae, with insights into horizontal gene transfer and methylation m6A motifs.

Detection of cytosine methylation in Burkholderia cenocepacia by single-molecule real-time sequencing and whole-genome bisulfite sequencing

Research on prokaryotic epigenetics, the study of heritable changes in gene expression independent of sequence changes, led to the identification of DNA methylation as a versatile regulator of diverse cellular processes. Methylation of adenine bases is often linked to regulation of gene expression in bacteria, but cytosine methylation is also frequently observed. In this study, we present a complete overview of the cytosine methylome in Burkholderia cenocepacia, an opportunistic respiratory pathogen in cystic fibrosis patients. Single-molecule real-time (SMRT) sequencing was used to map all 4mC-modified cytosines, as analysis of the predicted MTases in the B. cenocepacia genome revealed the presence of a 4mC-specific phage MTase, M.BceJII, targeting GGCC sequences. Methylation motif GCGGCCGC was identified, and out of 6850 motifs detected across the genome, 2051 (29.9 %) were methylated at the fifth position. Whole-genome bisulfite sequencing (WGBS) was performed to map 5mC methylation and 1635 5mC-modified cytosines were identified in CpG motifs. A comparison of the genomic positions of the modified bases called by each method revealed no overlap, which confirmed the authenticity of the detected 4mC and 5mC methylation by SMRT sequencing and WGBS, respectively. Large inter-strain variation of the 4mC-methylated cytosines was observed when B. cenocepacia strains J2315 and K56-2 were compared, which suggests that GGCC methylation patterns in B. cenocepacia are strain-specific. It seems likely that 4mC methylation of GGCC is not involved in regulation of gene expression but rather is a remnant of bacteriophage invasion, in which methylation of the phage genome was crucial for protection against restriction-modification systems of B. cenocepacia.

A positive perspective on DNA methylation: regulatory functions of DNA methylation outside of host defense in Gram-positive bacteria

The presence of post-replicative DNA methylation is pervasive among both prokaryotic and eukaryotic organisms. In bacteria, the study of DNA methylation has largely been in the context of restriction-modification systems, where DNA methylation serves to safeguard the chromosome against restriction endonuclease cleavage intended for invading DNA. There has been a growing recognition that the methyltransferase component of restriction-modification systems can also regulate gene expression, with important contributions to virulence factor gene expression in bacterial pathogens. Outside of restriction-modification systems, DNA methylation from orphan methyltransferases, which lack cognate restriction endonucleases, has been shown to regulate important processes, including DNA replication, DNA mismatch repair, and the regulation of gene expression. The majority of research and review articles have been focused on DNA methylation in the context of Gram-negative bacteria, with emphasis toward Escherichia coli, Caulobacter crescentus, and related Proteobacteria. Here we summarize the epigenetic functions of DNA methylation outside of host defense in Gram-positive bacteria, with a focus on the regulatory effects of both phase variable methyltransferases and DNA methyltransferases from traditional restriction-modification systems.

Next-generation DNA damage sequencing

Cellular DNA is constantly chemically altered by exogenous and endogenous agents. As all processes of life depend on the transmission of the genetic information, multiple biological processes exist to ensure genome integrity. Chemically damaged DNA has been linked to cancer and aging, therefore it is of great interest to map DNA damage formation and repair to elucidate the distribution of damage on a genome-wide scale. While the low abundance and inability to enzymatically amplify DNA damage are obstacles to genome-wide sequencing, new developments in the last few years have enabled high-resolution mapping of damaged bases. Recently, a number of DNA damage sequencing library construction strategies coupled to new data analysis pipelines allowed the mapping of specific DNA damage formation and repair at high and single nucleotide resolution. Strikingly, these advancements revealed that the distribution of DNA damage is heavily influenced by chromatin states and the binding of transcription factors. In the last seven years, these novel approaches have revealed new genomic maps of DNA damage distribution in a variety of organisms as generated by diverse chemical and physical DNA insults; oxidative stress, chemotherapeutic drugs, environmental pollutants, and sun exposure. Preferred sequences for damage formation and repair have been elucidated, thus making it possible to identify persistent weak spots in the genome as locations predicted to be vulnerable for mutation. As such, sequencing DNA damage will have an immense impact on our ability to elucidate mechanisms of disease initiation, and to evaluate and predict the efficacy of chemotherapeutic drugs.

In Vivo Genome and Methylome Adaptation of cag-Negative Helicobacter pylori during Experimental Human Infection

Exceptional genetic diversity and variability are hallmarks of Helicobacter pylori, but the biological role of this plasticity remains incompletely understood. Here, we had the rare opportunity to investigate the molecular evolution during the first weeks of H. pylori infection by comparing the genomes and epigenomes of H. pylori strain BCS 100 used to challenge human volunteers in a vaccine trial with those of bacteria reisolated from the volunteers 10 weeks after the challenge. The data provide molecular insights into the process of establishment of this highly versatile pathogen in 10 different human individual hosts, showing, for example, selection for changes in host-interaction molecules as well as changes in epigenetic methylation patterns. The data provide important clues to the early adaptation of H. pylori to new host niches after transmission, which we believe is vital to understand its success as a chronic pathogen and develop more efficient treatments and vaccines.

Multiple studies have demonstrated rapid bacterial genome evolution during chronic infection with Helicobacter pylori. In contrast, little was known about genetic changes during the first stages of infection, when selective pressure is likely to be highest. Using single-molecule, real-time (SMRT) and Illumina sequencing technologies, we analyzed genome and methylome evolution during the first 10 weeks of infection by comparing the cag pathogenicity island (cagPAI)-negative H. pylori challenge strain BCS 100 with pairs of H. pylori reisolates from gastric antrum and corpus biopsy specimens of 10 human volunteers who had been infected with this strain as part of a vaccine trial. Most genetic changes detected in the reisolates affected genes with a surface-related role or a predicted function in peptide uptake. Apart from phenotypic changes of the bacterial envelope, a duplication of the catalase gene was observed in one reisolate, which resulted in higher catalase activity and improved survival under oxidative stress conditions. The methylomes also varied in some of the reisolates, mostly by activity switching of phase-variable methyltransferase (MTase) genes. The observed in vivo mutation spectrum was remarkable for a very high proportion of nonsynonymous mutations. Although the data showed substantial within-strain genome diversity in the challenge strain, most antrum and corpus reisolates from the same volunteers were highly similar to each other, indicating that the challenge infection represents a major selective bottleneck shaping the transmitted population. Our findings suggest rapid in vivo selection of H. pylori during early-phase infection providing adaptation to different individuals by common mechanisms of genetic and epigenetic alterations.

SMRT Sequencing of Paramecium Bursaria Chlorella Virus-1 Reveals Diverse Methylation Stability in Adenines Targeted by Restriction Modification Systems

Chloroviruses (family Phycodnaviridae) infect eukaryotic, freshwater, unicellular green algae. A unique feature of these viruses is an abundance of DNA methyltransferases, with isolates dedicating up to 4.5% of their protein coding potential to these genes. This diversity highlights just one of the long-standing values of the chlorovirus model system; where group-wide epigenomic characterization might begin to elucidate the function(s) of DNA methylation in large dsDNA viruses. We characterized DNA modifications in the prototype chlorovirus, PBCV-1, using single-molecule real time (SMRT) sequencing (aka PacBio). Results were compared to total available sites predicted in silico based on DNA sequence alone. SMRT-software detected N6-methyl-adenine (m6A) at GATC and CATG recognition sites, motifs previously shown to be targeted by PBCV-1 DNA methyltransferases M.CviAI and M. CviAII, respectively. At the same time, PacBio analyses indicated that 10.9% of the PBCV-1 genome had large interpulse duration ratio (ipdRatio) values, the primary metric for DNA modification identification. These events represent 20.6x more sites than can be accounted for by all available adenines in GATC and CATG motifs, suggesting base or backbone modifications other than methylation might be present. To define methylation stability, we cross-compared methylation status of each GATC and CATG sequence in three biological replicates and found ∼81% of sites were stably methylated, while ∼2% consistently lack methylation. The remaining 17% of sites were stochastically methylated. When methylation status was analyzed for both strands of each target, we show that palindromes existed in completely non-methylated states, fully-methylated states, or hemi-methylated states, though GATC sites more often lack methylation than CATG sequences. Given that both sequences are targeted by not just methyltransferases, but by restriction endonucleases that are together encoded by PBCV-1 as virus-originating restriction modification (RM) systems, there is strong selective pressure to modify all target sites. The finding that most instances of non-methylation are associated with hemi-methylation is congruent with observations that hemi-methylated palindromes are resistant to cleavage by restriction endonucleases. However, sites where hemi-methylation is conserved might represent a unique regulatory function for PBCV-1. This study serves as a baseline for future investigation into the epigenomics of chloroviruses and their giant virus relatives.

A Deep Neural Network for Identifying DNA N4-Methylcytosine Sites

Motivation: N4-methylcytosine (4mC) plays an important role in host defense and transcriptional regulation. Accurate identification of 4mc sites provides a more comprehensive understanding of its biological effects. At present, the traditional machine learning algorithms are used in the research on 4mC sites prediction, but the complexity of the algorithms is relatively high, which is not suitable for the processing of large data sets, and the accuracy of prediction needs to be improved. Therefore, it is necessary to develop a new and effective method to accurately identify 4mC sites.

Results: In this work, we found a large number of 4mC sites and non 4mC sites of Caenorhabditis elegans (C. elegans) from the latest MethSMRT website, which greatly expanded the dataset of C. elegans, and developed a hybrid deep neural network framework named 4mcDeep-CBI, aiming to identify 4mC sites. In order to obtain the high latitude information of the feature, we input the preliminary extracted features into the Convolutional Neural Network (CNN) and Bidirectional Long Short Term Memory network (BLSTM) to generate advanced features. Taking the advanced features as algorithm input, we have proposed an integrated algorithm to improve feature representation. Experimental results on large new dataset show that the proposed predictor is able to achieve generally better performance in identifying 4mC sites as compared to the state-of-art predictor. Notably, this is the first study of identifying 4mC sites using deep neural network. Moreover, our model runs much faster than the state-of-art predictor.

Recent advances in the detection of base modifications using the Nanopore sequencer

DNA and RNA modifications have important functions, including the regulation of gene expression. Existing methods based on short-read sequencing for the detection of modifications show difficulty in determining the modification patterns of single chromosomes or an entire transcript sequence. Furthermore, the kinds of modifications for which detection methods are available are very limited. The Nanopore sequencer is a single-molecule, long-read sequencer that can directly sequence RNA as well as DNA. Moreover, the Nanopore sequencer detects modifications on long DNA and RNA molecules. In this review, we mainly focus on base modification detection in the DNA and RNA of mammals using the Nanopore sequencer. We summarize current studies of modifications using the Nanopore sequencer, detection tools using statistical tests or machine learning, and applications of this technology, such as analyses of open chromatin, DNA replication, and RNA metabolism.

Latest techniques to study DNA methylation

Bisulfite sequencing is a powerful technique to detect 5-methylcytosine in DNA that has immensely contributed to our understanding of epigenetic regulation in plants and animals. Meanwhile, research on other base modifications, including 6-methyladenine and 4-methylcytosine that are frequent in prokaryotes, has been impeded by the lack of a comparable technique. Bisulfite sequencing also suffers from a number of drawbacks that are difficult to surmount, among which DNA degradation, lack of specificity, or short reads with low sequence diversity. In this review, we explore the recent refinements to bisulfite sequencing protocols that enable targeting genomic regions of interest, detecting derivatives of 5-methylcytosine, and mapping single-cell methylomes. We then present the unique advantage of long-read sequencing in detecting base modifications in native DNA and highlight the respective strengths and weaknesses of PacBio and Nanopore sequencing for this application. Although analysing epigenetic data from long-read platforms remains challenging, the ability to detect various modified bases from a universal sample preparation, in addition to the mapping and phasing advantages of the longer read lengths, provide long-read sequencing with a decisive edge over short-read bisulfite sequencing for an expanding number of applications across kingdoms.

Genetic variation regulates the activation and specificity of Restriction-Modification systems in Neisseria gonorrhoeae

Restriction-Modification systems (RMS) are one of the main mechanisms of defence against foreign DNA invasion and can have an important role in the regulation of gene expression. The obligate human pathogen Neisseria gonorrhoeae carries one of the highest loads of RMS in its genome; between 13 to 15 of the three main types. Previous work has described their organization in the reference genome FA1090 and has inferred the associated methylated motifs. Here, we studied the structure of RMS and target methylated motifs in 25 gonococcal strains sequenced with Single Molecule Real-Time (SMRT) technology, which provides data on DNA modification. The results showed a variable picture of active RMS in different strains, with phase variation switching the activity of Type III RMS, and both the activity and specificity of a Type I RMS. Interestingly, the Dam methylase was found in place of the NgoAXI endonuclease in two of the strains, despite being previously thought to be absent in the gonococcus. We also identified the real methylation target of NgoAXII as 5′-GCAGA-3′, different from that previously described. Results from this work give further insights into the diversity and dynamics of RMS and methylation patterns in N. gonorrhoeae.

Microfluidic epigenomic mapping technologies for precision medicine

Epigenomic mapping of tissue samples generates critical insights into genome-wide regulations of gene activities and expressions during normal development and disease processes. Epigenomic profiling using a low number of cells produced by patient and mouse samples presents new challenges to biotechnologists. In this review, we first discuss the rationale and premise behind profiling epigenomes for precision medicine. We then examine the existing literature on applying microfluidics to facilitate low-input and high-throughput epigenomic profiling, with emphasis on technologies enabling interfacing with next-generation sequencing. We detail assays on studies of histone modifications, DNA methylation, 3D chromatin structures and non-coding RNAs. Finally, we discuss what the future may hold in terms of method development and translational potential.

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.

Rational development of transformation in Clostridium thermocellum ATCC 27405 via complete methylome analysis and evasion of native restriction-modification systems

A major barrier to both metabolic engineering and fundamental biological studies is the lack of genetic tools in most microorganisms. One example is Clostridium thermocellum ATCC 27405^T, where genetic tools are not available to help validate decades of hypotheses. A significant barrier to DNA transformation is restriction–modification systems, which defend against foreign DNA methylated differently than the host. To determine the active restriction–modification systems in this strain, we performed complete methylome analysis via single-molecule, real-time sequencing to detect 6-methyladenine and 4-methylcytosine and the rarely used whole-genome bisulfite sequencing to detect 5-methylcytosine. Multiple active systems were identified, and corresponding DNA methyltransferases were expressed from the Escherichia coli chromosome to mimic the C. thermocellum methylome. Plasmid methylation was experimentally validated and successfully electroporated into C. thermocellum ATCC 27405. This combined approach enabled genetic modification of the C. thermocellum-type strain and acts as a blueprint for transformation of other non-model microorganisms.

Deciphering bacterial epigenomes using modern sequencing technologies

Prokaryotic DNA contains three types of methylation: N6-methyladenine, N4-methylcytosine and 5-methylcytosine. The lack of tools to analyse the frequency and distribution of methylated residues in bacterial genomes has prevented a full understanding of their functions. Now, advances in DNA sequencing technology, including single-molecule, real-time sequencing and nanopore-based sequencing, have provided new opportunities for systematic detection of all three forms of methylated DNA at a genome-wide scale and offer unprecedented opportunities for achieving a more complete understanding of bacterial epigenomes. Indeed, as the number of mapped bacterial methylomes approaches 2,000, increasing evidence supports roles for methylation in regulation of gene expression, virulence and pathogen-host interactions.

Formation and biological consequences of 5-Formylcytosine in genomic DNA

5-Formyl-2'-deoxycytidine (5fdC) is a naturally occurring nucleobase that is broadly distributed in genomic DNA. 5fdC is produced via the oxidation of 5-methylcytosine (5mdC) by ten-eleven translocation enzyme (TET) and can be further converted to 5-carboxylcytosine (5cadC) by TET. Both 5fdC and 5cadC can be restored to dC by TDG-mediated base excision repair and direct deformylation/decarboxylation. Thus, 5fdC is considered an intermediate in the TET-mediated DNA demethylation pathway. 5fdC also alters the structure and stability of genomic DNA and affects genetic expression. This review summarizes the recent research on 5fdC, detailing its formation, detection and distribution, biological functions and transformation in cells. The challenges and future prospects to further explore the function and metabolism of 5fdC are briefly discussed at the end.

Detection of DNA base modifications by deep recurrent neural network on Oxford Nanopore sequencing data

DNA base modifications, such as C5-methylcytosine (5mC) and N6-methyldeoxyadenosine (6mA), are important types of epigenetic regulations. Short-read bisulfite sequencing and long-read PacBio sequencing have inherent limitations to detect DNA modifications. Here, using raw electric signals of Oxford Nanopore long-read sequencing data, we design DeepMod, a bidirectional recurrent neural network (RNN) with long short-term memory (LSTM) to detect DNA modifications. We sequence a human genome HX1 and a Chlamydomonas reinhardtii genome using Nanopore sequencing, and then evaluate DeepMod on three types of genomes (Escherichia coli, Chlamydomonas reinhardtii and human genomes). For 5mC detection, DeepMod achieves average precision up to 0.99 for both synthetically introduced and naturally occurring modifications. For 6mA detection, DeepMod achieves ~0.9 average precision on Escherichia coli data, and have improved performance than existing methods on Chlamydomonas reinhardtii data. In conclusion, DeepMod performs well for genome-scale detection of DNA modifications and will facilitate epigenetic analysis on diverse species.

DNA modification generates unique electric signals in Oxford Nanopore sequencing data but the signals can be complicated to decipher. Here, the authors develop a deep learning framework, DeepMod, to detect DNA base modifications including 5mC and 6mA using Nanopore sequencing data

Sources of artifact in measurements of 6mA and 4mC abundance in eukaryotic genomic DNA

BACKGROUND

Directed DNA methylation on N6-adenine (6mA), N4-cytosine (4mC), and C5-cytosine (5mC) can potentially increase DNA coding capacity and regulate a variety of biological functions. These modifications are relatively abundant in bacteria, occurring in about a percent of all bases of most bacteria. Until recently, 5mC and its oxidized derivatives were thought to be the only directed DNA methylation events in metazoa. New and more sensitive detection techniques (ultra-high performance liquid chromatography coupled with mass spectrometry (UHPLC-ms/ms) and single molecule real-time sequencing (SMRTseq)) have suggested that 6mA and 4mC modifications could be present in a variety of metazoa.

RESULTS

Here, we find that both of these techniques are prone to inaccuracies, which overestimate DNA methylation concentrations in metazoan genomic DNA. Artifacts can arise from methylated bacterial DNA contamination of enzyme preparations used to digest DNA and contaminating bacterial DNA in eukaryotic DNA preparations. Moreover, DNA sonication introduces a novel modified base from 5mC that has a retention time near 4mC that can be confused with 4mC. Our analyses also suggest that SMRTseq systematically overestimates 4mC in prokaryotic and eukaryotic DNA and 6mA in DNA samples in which it is rare. Using UHPLC-ms/ms designed to minimize and subtract artifacts, we find low to undetectable levels of 4mC and 6mA in genomes of representative worms, insects, amphibians, birds, rodents and primates under normal growth conditions. We also find that mammalian cells incorporate exogenous methylated nucleosides into their genome, suggesting that a portion of 6mA modifications could derive from incorporation of nucleosides from bacteria in food or microbiota. However, gDNA samples from gnotobiotic mouse tissues found rare (0.9-3.7 ppm) 6mA modifications above background.

CONCLUSIONS

Altogether these data demonstrate that 6mA and 4mC are rarer in metazoa than previously reported, and highlight the importance of careful sample preparation and measurement, and need for more accurate sequencing techniques.

Compilation of Modern Technologies To Map Genome-Wide Cytosine Modifications in DNA

Over the past few decades, various DNA modification detection methods have been developed; many of the high-resolution methods are based on bisulfite treatment, which leads to DNA degradation, to a degree. Thus, novel bisulfite-free approaches have been developed in recent years and shown to be useful for epigenome analysis in otherwise difficult-to-handle, but important, DNA samples, such as hmC-seal and hmC-CATCH. Herein, an overview of advances in the development of epigenome sequencing methods for these important DNA modifications is provided.

NanoMod: a computational tool to detect DNA modifications using Nanopore long-read sequencing data

BACKGROUND

Recent advances in single-molecule sequencing techniques, such as Nanopore sequencing, improved read length, increased sequencing throughput, and enabled direct detection of DNA modifications through the analysis of raw signals. These DNA modifications include naturally occurring modifications such as DNA methylations, as well as modifications that are introduced by DNA damage or through synthetic modifications to one of the four standard nucleotides.

METHODS

To improve the performance of detecting DNA modifications, especially synthetically introduced modifications, we developed a novel computational tool called NanoMod. NanoMod takes raw signal data on a pair of DNA samples with and without modified bases, extracts signal intensities, performs base error correction based on a reference sequence, and then identifies bases with modifications by comparing the distribution of raw signals between two samples, while taking into account of the effects of neighboring bases on modified bases ("neighborhood effects").

RESULTS

We evaluated NanoMod on simulation data sets, based on different types of modifications and different magnitudes of neighborhood effects, and found that NanoMod outperformed other methods in identifying known modified bases. Additionally, we demonstrated superior performance of NanoMod on an E. coli data set with 5mC (5-methylcytosine) modifications.

CONCLUSIONS

In summary, NanoMod is a flexible tool to detect DNA modifications with single-base resolution from raw signals in Nanopore sequencing, and will facilitate large-scale functional genomics experiments that use modified nucleotides.

Single-molecule sequencing detection of N6-methyladenine in microbial reference materials

The DNA base modification N6-methyladenine (m6A) is involved in many pathways related to the survival of bacteria and their interactions with hosts. Nanopore sequencing offers a new, portable method to detect base modifications. Here, we show that a neural network can improve m6A detection at trained sequence contexts compared to previously published methods using deviations between measured and expected current values as each adenine travels through a pore. The model, implemented as the mCaller software package, can be extended to detect known or confirm suspected methyltransferase target motifs based on predictions of methylation at untrained contexts. We use PacBio, Oxford Nanopore, methylated DNA immunoprecipitation sequencing (MeDIP-seq), and whole-genome bisulfite sequencing data to generate and orthogonally validate methylomes for eight microbial reference species. These well-characterized microbial references can serve as controls in the development and evaluation of future methods for the identification of base modifications from single-molecule sequencing data.

Study of the whole genome, methylome and transcriptome of Cordyceps militaris

The complete genome of Cordyceps militaris was sequenced using single-molecule real-time (SMRT) sequencing technology at a coverage over 300×. The genome size was 32.57 Mb, and 14 contigs ranging from 0.35 to 4.58 Mb with an N50 of 2.86 Mb were assembled, including 4 contigs with telomeric sequences on both ends and an additional 8 contigs with telomeric sequences on either the 5′ or 3′ end. A methylome database of the genome was constructed using SMRT and m4C and m6A methylated nucleotides, and many unknown modification types were identified. The major m6A methylation motif is GA and GGAG, and the major m4C methylation motif is GC or CG/GC. In the C. militaris genome DNA, there were four types of methylated nucleotides that we confirmed using high-resolution LCMS-IT-TOF. Using PacBio Iso-Seq, a total of 31,133 complete cDNA sequences were obtained in the fruiting body. The conserved domains of the nontranscribed regions of the genome include TATA boxes, which are the initial regions of genome replication. There were 406 structural variants between the HN and CM01 strains, and there were 1,114 structural variants between the HN and ATCC strains.

Insights into the Biochemistry, Evolution, and Biotechnological Applications of the Ten-Eleven Translocation (TET) Enzymes

A tight link exists between patterns of DNA methylation at carbon 5 of cytosine and differential gene expression in mammalian tissues. Indeed, aberrant DNA methylation results in various human diseases, including neurologic and immune disorders, and contributes to the initiation and progression of various cancers. Proper DNA methylation depends on the fidelity and control of the underlying mechanisms that write, maintain, and erase these epigenetic marks. In this Perspective, we address one of the key players in active demethylation: the ten-eleven translocation enzymes or TETs. These enzymes belong to the Fe²⁺/α-ketoglutarate-dependent dioxygenase superfamily and iteratively oxidize 5-methylcytosine (5mC) in DNA to produce 5-hydroxymethylcytosine, 5-formylcytosine, and 5-carboxycytosine. The latter three bases may convey additional layers of epigenetic information in addition to being intermediates in active demethylation. Despite the intense interest in understanding the physiological roles TETs play in active demethylation and cell regulation, less has been done, in comparison, to illuminate details of the chemistry and factors involved in regulating the three-step oxidation mechanism. Herein, we focus on what is known about the biochemical features of TETs and explore questions whose answers will lead to a more detailed understanding of the in vivo modus operandi of these enzymes. We also summarize the membership and evolutionary history of the TET/JBP family and highlight the prokaryotic homologues as a reservoir of potentially diverse functionalities awaiting discovery. Finally, we spotlight sequencing methods that utilize TETs for mapping 5mC and its oxidation products in genomic DNA and comment on possible improvements in these approaches.

Advances in detection and quantification of methylcytosine and its derivatives

Methylation of the fifth carbon atom in cytosine is an epigenetic modification of deoxyribonucleic acid that plays important roles in numerous cellular processes and disease pathogenesis. Three additional states of cytosine, that is, 5-hydroxymethylcytosine, 5-formylcytosine and 5-carboxylcytosine, have been identified and associated with the diagnosis and/or prognosis of diseases. However, accurate measurement of those intermediates is a challenge since their global levels are relatively low. A number of innovative methods have been developed to detect and quantify these compounds in biological samples, such as blood, tissue and urine, etc. This review focuses on recent advancement in detection and quantification of four cytosine modifications, based on which, the development, diagnosis, and prognosis of diseases could be monitored through non-invasive procedures.

Identification of DNA Base Modifications by Means of Pacific Biosciences RS Sequencing Technology

Whole phage genomes can be sequenced readily using one or a combination of next generation sequencing (NGS) technologies. One of the most recently developed NGS platforms, the so-called Single-Molecule Real-Time (SMRT) sequencing approach provided by the PacBio RS platform, is particularly useful in providing complete (i.e., un-gapped) genome sequences, but differs from other technologies in that the platform also allows for downstream analysis to identify nucleotides that have been modified by DNA methylation. Here, we describe the methodological approach for the detection of genomic methylation motifs by means of SMRT sequencing.

Epigenetic modification of nucleic acids: from basic studies to medical applications

The epigenetic modification of nucleic acids represents one of the most significant areas of study in the field of nucleic acids because it makes gene regulation more complex and heredity more complicated, thus indicating its profound impact on aspects of heredity, growth, and diseases. The recent characterization of epigenetic modifications of DNA and RNA using chemical labelling strategies has promoted the discovery of these modifications, and the newly developed single-base or single-cell resolution mapping strategies have enabled large-scale epigenetic studies in eukaryotes. Due to these technological breakthroughs, several new epigenetic marks have been discovered that have greatly extended the scope and impact of epigenetic modifications in nucleic acids over the past few years. Because epigenetics is reversible and susceptible to environmental factors, it could potentially be a promising direction for clinical medicine research. In this review, we have comprehensively discussed how these epigenetic marks are involved in disease, including the pathogenesis, prevention, diagnosis and treatment of disease. These findings have revealed that the epigenetic modification of nucleic acids has considerable significance in various areas from methodology to clinical medicine and even in biomedical applications.

Resolving the complete genome of Kuenenia stuttgartiensis from a membrane bioreactor enrichment using Single-Molecule Real-Time sequencing

Anaerobic ammonium-oxidizing (anammox) bacteria are a group of strictly anaerobic chemolithoautotrophic microorganisms. They are capable of oxidizing ammonium to nitrogen gas using nitrite as a terminal electron acceptor, thereby facilitating the release of fixed nitrogen into the atmosphere. The anammox process is thought to exert a profound impact on the global nitrogen cycle and has been harnessed as an environment-friendly method for nitrogen removal from wastewater. In this study, we present the first closed genome sequence of an anammox bacterium, Kuenenia stuttgartiensis MBR1. It was obtained through Single-Molecule Real-Time (SMRT) sequencing of an enrichment culture constituting a mixture of at least two highly similar Kuenenia strains. The genome of the novel MBR1 strain is different from the previously reported Kuenenia KUST reference genome as it contains numerous structural variations and unique genomic regions. We find new proteins, such as a type 3b (sulf)hydrogenase and an additional copy of the hydrazine synthase gene cluster. Moreover, multiple copies of ammonium transporters and proteins regulating nitrogen uptake were identified, suggesting functional differences in metabolism. This assembly, including the genome-wide methylation profile, provides a new foundation for comparative and functional studies aiming to elucidate the biochemical and metabolic processes of these organisms.

Characterizing the DNA Methyltransferases of Haloferax volcanii via Bioinformatics, Gene Deletion, and SMRT Sequencing

DNA methyltransferases (MTases), which catalyze the methylation of adenine and cytosine bases in DNA, can occur in bacteria and archaea alongside cognate restriction endonucleases (REases) in restriction-modification (RM) systems or independently as orphan MTases. Although DNA methylation and MTases have been well-characterized in bacteria, research into archaeal MTases has been limited. A previous study examined the genomic DNA methylation patterns (methylome) of the halophilic archaeon Haloferax volcanii, a model archaeal system which can be easily manipulated in laboratory settings, via single-molecule real-time (SMRT) sequencing and deletion of a putative MTase gene (HVO_A0006). In this follow-up study, we deleted other putative MTase genes in H. volcanii and sequenced the methylomes of the resulting deletion mutants via SMRT sequencing to characterize the genes responsible for DNA methylation. The results indicate that deletion of putative RM genes HVO_0794, HVO_A0006, and HVO_A0237 in a single strain abolished methylation of the sole cytosine motif in the genome (C^m4TAG). Amino acid alignments demonstrated that HVO_0794 shares homology with characterized cytosine CTAG MTases in other organisms, indicating that this MTase is responsible for C^m4TAG methylation in H. volcanii. The CTAG motif has high density at only one of the origins of replication, and there is no relative increase in CTAG motif frequency in the genome of H. volcanii, indicating that CTAG methylation might not have effectively taken over the role of regulating DNA replication and mismatch repair in the organism as previously predicted. Deletion of the putative Type I RM operon rmeRMS (HVO_2269-2271) resulted in abolished methylation of the adenine motif in the genome (GCA^m6BN₆VTGC). Alignments of the MTase (HVO_2270) and site specificity subunit (HVO_2271) demonstrate homology with other characterized Type I MTases and site specificity subunits, indicating that the rmeRMS operon is responsible for adenine methylation in H. volcanii. Together with HVO_0794, these genes appear to be responsible for all detected methylation in H. volcanii, even though other putative MTases (HVO_C0040, HVO_A0079) share homology with characterized MTases in other organisms. We also report the construction of a multi-RM deletion mutant (ΔRM), with multiple RM genes deleted and with no methylation detected via SMRT sequencing, which we anticipate will be useful for future studies on DNA methylation in H. volcanii.