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

Information Gradient among Nucleotide Sequences of Essential RNAs from an Evolutionary Perspective

We hypothesize that the first ancestral "protocell" molecular structures, i.e., the first RNAs and peptides that gradually transformed into real cells once the Earth had cooled sufficiently for organic molecules to appear, have left traces in the RNAs and the genes in present cells. We propose a circular RNA that could have been one of these ancestral structures whose vestigial pentameric subsequences would mark the evolution from this key moment when the protocells began to join with living organisms. In particular, we propose that, in present RNAs (ribosomal or messenger), which play an important role in the metabolism of current cells, we look for traces of the proposed primitive structure in the form of pentamers (or longer fragments) that belong to their nucleotide sequence. The result obtained can be summarized in the existence of a gradient of occurrence of such pentamers, with a high frequency for the most vital functions (protein synthesis, nucleic synthesis, cell respiration, etc.). This gradient is also visible between organisms, from the oldest (Archaea) to the most recent (Eukaryotes) in the evolution of species.

Regulatory sequence-based discovery of anti-defense genes in archaeal viruses

In silico identification of viral anti-CRISPR proteins (Acrs) has relied largely on the guilt-by-association method using known Acrs or anti-CRISPR associated proteins (Acas) as the bait. However, the low number and limited spread of the characterized archaeal Acrs and Aca hinders our ability to identify Acrs using guilt-by-association. Here, based on the observation that the few characterized archaeal Acrs and Aca are transcribed immediately post viral infection, we hypothesize that these genes, and many other unidentified anti-defense genes (ADG), are under the control of conserved regulatory sequences including a strong promoter, which can be used to predict anti-defense genes in archaeal viruses. Using this consensus sequence based method, we identify 354 potential ADGs in 57 archaeal viruses and 6 metagenome-assembled genomes. Experimental validation identified a CRISPR subtype I-A inhibitor and the first virally encoded inhibitor of an archaeal toxin-antitoxin based immune system. We also identify regulatory proteins potentially akin to Acas that can facilitate further identification of ADGs combined with the guilt-by-association approach. These results demonstrate the potential of regulatory sequence analysis for extensive identification of ADGs in viruses of archaea and bacteria.

Restriction modification systems in archaea: A panoramic outlook

Restriction modification (RM) systems are one of the ubiquitous yet primitive defense responses employed by bacteria and archaea with the primary role of safeguarding themselves against invading bacteriophages. Protection of the host occurs by the cleavage of the invading foreign DNA via restriction endonucleases with concomitant methylation of host DNA with the aid of a methyltransferase counterpart. RM systems have been extensively studied in bacteria, however, in the case of archaea there are limited reports of RM enzymes that are investigated to date owing to their inhospitable growth demands. This review aims to broaden the knowledge about what is known about the diversity of RM systems in archaea and encapsulate the current knowledge on restriction and modification enzymes characterized in archaea so far and the role of RM systems in the milieu of archaeal biology.

Interplay of intracellular and trans-cellular DNA methylation in natural archaeal consortia

DNA methylation serves a variety of functions across all life domains. In this study, we investigated archaeal methylomics within a tripartite xylanolytic halophilic consortium. This consortium includes Haloferax lucertense SVX82, Halorhabdus sp. SVX81, and an ectosymbiotic Candidatus Nanohalococcus occultus SVXNc, a nano-sized archaeon from the DPANN superphylum. We utilized PacBio SMRT and Illumina cDNA sequencing to analyse samples from consortia of different compositions for methylomics and transcriptomics. Endogenous cTAG methylation, typical of Haloferax, was accompanied in this strain by methylation at four other motifs, including GDGcHC methylation, which is specific to the ectosymbiont. Our analysis of the distribution of methylated and unmethylated motifs suggests that autochthonous cTAG methylation may influence gene regulation. The frequency of GRAGAaG methylation increased in highly expressed genes, while CcTTG and GTCGaGG methylation could be linked to restriction-modification (RM) activity. Generally, the RM activity might have been reduced during the evolution of this archaeon to balance the protection of cells from intruders, the reduction of DNA damage due to self-restriction in stressful environments, and the benefits of DNA exchange under extreme conditions. Our methylomics, transcriptomics and complementary electron cryotomography (cryo-ET) data suggest that the nanohaloarchaeon exports its methyltransferase to methylate the Haloferax genome, unveiling a new aspect of the interaction between the symbiont and its host.

The extremophile Picrophilus torridus carries a DNA adenine methylase M.PtoI that is part of a Type I restriction-modification system

DNA methylation events mediated by orphan methyltransferases modulate various cellular processes like replication, repair and transcription. Bacteria and archaea also harbor DNA methyltransferases that are part of restriction-modification systems, which serve to protect the host genome from being cleaved by the cognate restriction enzyme. While DNA methylation has been exhaustively investigated in bacteria it remains poorly understood in archaea. Picrophilus torridus is a euryarchaeon that can thrive under conditions of extremely low pH (0.7), and thus far no reports have been published regarding DNA methylation in this extremophile. This study reports the first experimentation examining DNA methylation in P. torridus. We find the genome to carry methylated adenine (m6A) but not methylated cytosine (m5C) residues. The m6A modification is absent at GATC sites, indicating the absence of an active Dam methylase even though the dam gene has been annotated in the genome sequence. Two other methylases have also been annotated in the P. torridus genome sequence. One of these is a part of a Type I restriction-modification system. Considering that all Type I modification methylases characterized to date target adenine residues, the modification methylase of this Type I system has been examined. The genes encoding the S subunit (that is responsible for DNA recognition) and M subunit (that is responsible for DNA methylation) have been cloned and the recombinant protein purified from E.coli, and regions involved in M-S interactions have been identified. The M.PtoI enzyme harbors all the motifs that typify Type I modification methylases, and displays robust adenine methylation in in vitro assays under a variety of conditions. Interestingly, magnesium is essential for enzyme activity. The enzyme displays substrate inhibition at higher concentrations of AdoMet. Mutational analyses reveal that Motif I plays a role in AdoMet binding, and Motif IV is critical for methylation activity. The data presented here lays the foundation for further research in the area of DNA methylation and restriction-modification research in this most unusual microorganism.

DNA Methyltransferases: From Evolution to Clinical Applications

DNA methylation is an epigenetic mark that living beings have used in different environments. The MTases family catalyzes DNA methylation. This process is conserved from archaea to eukaryotes, from fertilization to every stage of development, and from the early stages of cancer to metastasis. The family of DNMTs has been classified into DNMT1, DNMT2, and DNMT3. Each DNMT has been duplicated or deleted, having consequences on DNMT structure and cellular function, resulting in a conserved evolutionary reaction of DNA methylation. DNMTs are conserved in the five kingdoms of life: bacteria, protists, fungi, plants, and animals. The importance of DNMTs in whether methylate or not has a historical adaptation that in mammals has been discovered in complex regulatory mechanisms to develop another padlock to genomic insurance stability. The regulatory mechanisms that control DNMTs expression are involved in a diversity of cell phenotypes and are associated with pathologies transcription deregulation. This work focused on DNA methyltransferases, their biology, functions, and new inhibitory mechanisms reported. We also discuss different approaches to inhibit DNMTs, the use of non-coding RNAs and nucleoside chemical compounds in recent studies, and their importance in biological, clinical, and industry research.

The Hyperthermophilic Restriction-Modification Systems of Thermococcus kodakarensis Protect Genome Integrity

Thermococcus kodakarensis (T. kodakarensis), a hyperthermophilic, genetically accessible model archaeon, encodes two putative restriction modification (R-M) defense systems, TkoI and TkoII. TkoI is encoded by TK1460 while TkoII is encoded by TK1158. Bioinformative analysis suggests both R-M enzymes are large, fused methyltransferase (MTase)-endonuclease polypeptides that contain both restriction endonuclease (REase) activity to degrade foreign invading DNA and MTase activity to methylate host genomic DNA at specific recognition sites. In this work, we demonsrate T. kodakarensis strains deleted for either or both R-M enzymes grow more slowly but display significantly increased competency compared to strains with intact R-M systems, suggesting that both TkoI and TkoII assist in maintenance of genomic integrity in vivo and likely protect against viral- or plasmid-based DNA transfers. Pacific Biosciences single molecule real-time (SMRT) sequencing of T. kodakarensis strains containing both, one or neither R-M systems permitted assignment of the recognition sites for TkoI and TkoII and demonstrated that both R-M enzymes are TypeIIL; TkoI and TkoII methylate the N⁶ position of adenine on one strand of the recognition sequences GTGAAG and TTCAAG, respectively. Further in vitro biochemical characterization of the REase activities reveal TkoI and TkoII cleave the DNA backbone GTGAAG(N)₂₀/(N)₁₈ and TTCAAG(N)₁₀/(N)₈, respectively, away from the recognition sequences, while in vitro characterization of the MTase activities reveal transfer of tritiated S-adenosyl methionine by TkoI and TkoII to their respective recognition sites. Together these results demonstrate TkoI and TkoII restriction systems are important for protecting T. kodakarensis genome integrity from invading foreign DNA.

Insights into gene expression changes under conditions that facilitate horizontal gene transfer (mating) of a model archaeon

Horizontal gene transfer is a means by which bacteria, archaea, and eukaryotes are able to trade DNA within and between species. While there are a variety of mechanisms through which this genetic exchange can take place, one means prevalent in the archaeon Haloferax volcanii involves the transient formation of cytoplasmic bridges between cells and is referred to as mating. This process can result in the exchange of very large fragments of DNA between the participating cells. Genes governing the process of mating, including triggers to initiate mating, mechanisms of cell fusion, and DNA exchange, have yet to be characterized. We used a transcriptomic approach to gain a more detailed knowledge of how mating might transpire. By examining the differential expression of genes expressed in cells harvested from mating conditions on a filter over time and comparing them to those expressed in a shaking culture, we were able to identify genes and pathways potentially associated with mating. These analyses provide new insights into both the mechanisms and barriers of mating in Hfx. volcanii.

4-Methylcytosine DNA modification is critical for global epigenetic regulation and virulence in the human pathogen Leptospira interrogans

In bacteria, DNA methylation can be facilitated by 'orphan' DNA methyltransferases lacking cognate restriction endonucleases, but whether and how these enzymes control key cellular processes are poorly understood. The effects of a specific modification, 4-methylcytosine (4mC), are even less clear, as this epigenetic marker is unique to bacteria and archaea, whereas the bulk of epigenetic research is currently performed on eukaryotes. Here, we characterize a 4mC methyltransferase from the understudied pathogen Leptospira spp. Inactivating this enzyme resulted in complete abrogation of CTAG motif methylation, leading to genome-wide dysregulation of gene expression. Mutants exhibited growth defects, decreased adhesion to host cells, higher susceptibility to LPS-targeting antibiotics, and, importantly, were no longer virulent in an acute infection model. Further investigation resulted in the discovery of at least one gene, that of an ECF sigma factor, whose transcription was altered in the methylase mutant and, subsequently, by mutation of the CTAG motifs in the promoter of the gene. The genes that comprise the regulon of this sigma factor were, accordingly, dysregulated in the methylase mutant and in a strain overexpressing the sigma factor. Our results highlight the importance of 4mC in Leptospira physiology, and suggest the same of other understudied species.

Haloferax volcanii-a model archaeon for studying DNA replication and repair

The tree of life shows the relationship between all organisms based on their common ancestry. Until 1977, it comprised two major branches: prokaryotes and eukaryotes. Work by Carl Woese and other microbiologists led to the recategorization of prokaryotes and the proposal of three primary domains: Eukarya, Bacteria and Archaea. Microbiological, genetic and biochemical techniques were then needed to study the third domain of life. Haloferax volcanii, a halophilic species belonging to the phylum Euryarchaeota, has provided many useful tools to study Archaea, including easy culturing methods, genetic manipulation and phenotypic screening. This review will focus on DNA replication and DNA repair pathways in H. volcanii, how this work has advanced our knowledge of archaeal cellular biology, and how it may deepen our understanding of bacterial and eukaryotic processes.

DNA Methylation Epigenetically Regulates Gene Expression in Burkholderia cenocepacia and Controls Biofilm Formation, Cell Aggregation, and Motility

CF patients diagnosed with Burkholderia cenocepacia infections often experience rapid deterioration of lung function, known as cepacia syndrome. B. cenocepacia has a large multireplicon genome, and much remains to be learned about regulation of gene expression in this organism. From studies in other (model) organisms, it is known that epigenetic changes through DNA methylation play an important role in this regulation. The identification of B. cenocepacia genes of which the expression is regulated by DNA methylation and identification of the regulatory systems involved in this methylation are likely to advance the biological understanding of B. cenocepacia cell adaptation via epigenetic regulation. In time, this might lead to novel approaches to tackle B. cenocepacia infections in CF patients.

Respiratory tract infections by the opportunistic pathogen Burkholderia cenocepacia often lead to severe lung damage in cystic fibrosis (CF) patients. New insights in how to tackle these infections might emerge from the field of epigenetics, as DNA methylation is an important regulator of gene expression. The present study focused on two DNA methyltransferases (MTases) in B. cenocepacia strains J2315 and K56-2 and their role in regulating gene expression. In silico predicted DNA MTase genes BCAL3494 and BCAM0992 were deleted in both strains, and the phenotypes of the resulting deletion mutants were studied: deletion mutant ΔBCAL3494 showed changes in biofilm structure and cell aggregation, while ΔBCAM0992 was less motile. B. cenocepacia wild-type cultures treated with sinefungin, a known DNA MTase inhibitor, exhibited the same phenotype as DNA MTase deletion mutants. Single-molecule real-time sequencing was used to characterize the methylome of B. cenocepacia, including methylation at the origin of replication, and motifs CACAG and GTWWAC were identified as targets of BCAL3494 and BCAM0992, respectively. All genes with methylated motifs in their putative promoter region were identified, and qPCR experiments showed an upregulation of several genes, including biofilm- and motility-related genes, in MTase deletion mutants with unmethylated motifs, explaining the observed phenotypes in these mutants. In summary, our data confirm that DNA methylation plays an important role in regulating the expression of B. cenocepacia genes involved in biofilm formation, cell aggregation, and motility.

IMPORTANCE CF patients diagnosed with Burkholderia cenocepacia infections often experience rapid deterioration of lung function, known as cepacia syndrome. B. cenocepacia has a large multireplicon genome, and much remains to be learned about regulation of gene expression in this organism. From studies in other (model) organisms, it is known that epigenetic changes through DNA methylation play an important role in this regulation. The identification of B. cenocepacia genes of which the expression is regulated by DNA methylation and identification of the regulatory systems involved in this methylation are likely to advance the biological understanding of B. cenocepacia cell adaptation via epigenetic regulation. In time, this might lead to novel approaches to tackle B. cenocepacia infections in CF patients.

The complex phylogenetic relationships of a 4mC/6mA DNA methyltransferase in prokaryotes

DNA methyltransferases are proteins that modify DNA via attachment of methyl groups to nucleobases and are ubiquitous across the bacterial, archaeal, and eukaryotic domains of life. Here, we investigated the complex evolutionary history of the large and consequential 4mC/6mA DNA methyltransferase protein family using phylogenetic reconstruction of amino acid sequences. We present a well-supported phylogeny of this family based on systematic sampling of taxa across superphyla of bacteria and archaea. We compared the phylogeny to a current representation of the species tree of life and found that the 4mC/6mA methyltransferase family has a strikingly complex evolutionary history that likely began sometime after the last universal common ancestor of life diverged into the bacterial and archaeal lineages and probably involved many horizontal gene transfers within and between domains. Despite the complexity of its evolutionary history, we inferred that only one significant shift in molecular evolutionary rate characterizes the diversification of this protein family.

Comparative Genomics of Two New HF1-like Haloviruses

Few genomes of the HF1-group of viruses are currently available, and further examples would enhance the understanding of their evolution, improve their gene annotation, and assist in understanding gene function and regulation. Two novel HF1-group haloviruses, Serpecor1 and Hardycor2, were recovered from widely separated hypersaline lakes in Australia. Both are myoviruses with linear dsDNA genomes and infect the haloarchaeon Halorubrum coriense. Both genomes possess long, terminal direct repeat (TDR) sequences (320 bp for Serpecor1 and 306 bp for Hardycor2). The Serpecor1 genome is 74,196 bp in length, 57.0% G+C, and has 126 annotated coding sequences (CDS). Hardycor2 has a genome of 77,342 bp, 55.6% G+C, and 125 annotated CDS. They show high nucleotide sequence similarity to each other (78%) and with HF1 (>75%), and carry similar intergenic repeat (IR) sequences to those originally described in HF1 and HF2. Hardycor2 carries a DNA methyltransferase gene in the same genomic neighborhood as the methyltransferase genes of HF1, HF2 and HRTV-5, but is in the opposite orientation, and the inferred proteins are only distantly related. Comparative genomics allowed us to identify the candidate genes mediating cell attachment. The genomes of Serpecor1 and Hardycor2 encode numerous small proteins carrying one or more CxxC motifs, a signature feature of zinc-finger domain proteins that are known to participate in diverse biomolecular interactions.

Whole-genome comparison between the type strain of Halobacterium salinarum (DSM 3754T ) and the laboratory strains R1 and NRC-1

Halobacterium salinarum is an extremely halophilic archaeon that is widely distributed in hypersaline environments and was originally isolated as a spoilage organism of salted fish and hides. The type strain 91-R6 (DSM 3754T ) has seldom been studied and its genome sequence has only recently been determined by our group. The exact relationship between the type strain and two widely used model strains, NRC-1 and R1, has not been described before. The genome of Hbt. salinarum strain 91-R6 consists of a chromosome (2.17 Mb) and two large plasmids (148 and 102 kb, with 39,230 bp being duplicated). Cytosine residues are methylated (m4 C) within CTAG motifs. The genomes of type and laboratory strains are closely related, their chromosomes sharing average nucleotide identity (ANIb) values of 98% and in silico DNA-DNA hybridization (DDH) values of 95%. The chromosomes are completely colinear, do not show genome rearrangement, and matching segments show <1% sequence difference. Among the strain-specific sequences are three large chromosomal replacement regions (>10 kb). The well-studied AT-rich island (61 kb) of the laboratory strains is replaced by a distinct AT-rich sequence (47 kb) in 91-R6. Another large replacement (91-R6: 78 kb, R1: 44 kb) codes for distinct homologs of proteins involved in motility and N-glycosylation. Most (107 kb) of plasmid pHSAL1 (91-R6) is very closely related to part of plasmid pHS3 (R1) and codes for essential genes (e.g. arginine-tRNA ligase and the pyrimidine biosynthesis enzyme aspartate carbamoyltransferase). Part of pHS3 (42.5 kb total) is closely related to the largest strain-specific sequence (164 kb) in the type strain chromosome. Genome sequencing unraveled the close relationship between the Hbt. salinarum type strain and two well-studied laboratory strains at the DNA and protein levels. Although an independent isolate, the type strain shows a remarkably low evolutionary difference to the laboratory strains.

The Role of Archaeal Chromatin in Transcription

Genomic organization impacts accessibility and movement of information processing systems along DNA. DNA-bound proteins dynamically dictate gene expression and provide regulatory potential to tune transcription rates to match ever-changing environmental conditions. Archaeal genomes are typically small, circular, gene dense, and organized either by histone proteins that are homologous to their eukaryotic counterparts, or small basic proteins that function analogously to bacterial nucleoid proteins. We review here how archaeal genomes are organized and how such organization impacts archaeal gene expression, focusing on conserved DNA-binding proteins within the clade and the factors that are known to impact transcription initiation and elongation within protein-bound genomes.

Methylomes of Two Extremely Halophilic Archaea Species, Haloarcula marismortui and Haloferax mediterranei

The genomes of two extremely halophilic Archaea species, Haloarcula marismortui and Haloferax mediterranei, were sequenced using single-molecule real-time sequencing. The ∼4-Mbp genomes are GC rich with multiple large plasmids and two 4-methyl-cytosine patterns. Methyl transferases were incorporated into the Restriction Enzymes Database (REBASE), and gene annotation was incorporated into the Haloarchaeal Genomes Database (HaloWeb).

The genomes of two extremely halophilic Archaea species, Haloarcula marismortui and Haloferax mediterranei, were sequenced using single-molecule real-time sequencing. The ∼4-Mbp genomes are GC rich with multiple large plasmids and two 4-methyl-cytosine patterns. Methyl transferases were incorporated into the Restriction Enzymes Database (REBASE), and gene annotation was incorporated into the Haloarchaeal Genomes Database (HaloWeb).

The Patchy Distribution of Restriction⁻Modification System Genes and the Conservation of Orphan Methyltransferases in Halobacteria

Restriction⁻modification (RM) systems in bacteria are implicated in multiple biological roles ranging from defense against parasitic genetic elements, to selfish addiction cassettes, and barriers to gene transfer and lineage homogenization. In bacteria, DNA-methylation without cognate restriction also plays important roles in DNA replication, mismatch repair, protein expression, and in biasing DNA uptake. Little is known about archaeal RM systems and DNA methylation. To elucidate further understanding for the role of RM systems and DNA methylation in Archaea, we undertook a survey of the presence of RM system genes and related genes, including orphan DNA methylases, in the halophilic archaeal class Halobacteria. Our results reveal that some orphan DNA methyltransferase genes were highly conserved among lineages indicating an important functional constraint, whereas RM systems demonstrated patchy patterns of presence and absence. This irregular distribution is due to frequent horizontal gene transfer and gene loss, a finding suggesting that the evolution and life cycle of RM systems may be best described as that of a selfish genetic element. A putative target motif (CTAG) of one of the orphan methylases was underrepresented in all of the analyzed genomes, whereas another motif (GATC) was overrepresented in most of the haloarchaeal genomes, particularly in those that encoded the cognate orphan methylase.