<?xml version="1.0" encoding="utf-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" article-type="research-article" dtd-version="2.3" xml:lang="EN">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Microbiol.</journal-id>
<journal-title>Frontiers in Microbiology</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Microbiol.</abbrev-journal-title>
<issn pub-type="epub">1664-302X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fmicb.2023.1126750</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Microbiology</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>The extremophile <italic>Picrophilus torridus</italic> carries a DNA adenine methylase M.PtoI that is part of a Type I restriction-modification system</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author"><name><surname>Gulati</surname>
<given-names>Pallavi</given-names></name><xref rid="aff1" ref-type="aff"><sup>1</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/2144254/overview"/>
</contrib>
<contrib contrib-type="author"><name><surname>Singh</surname>
<given-names>Ashish</given-names></name><xref rid="aff1" ref-type="aff"><sup>1</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/2144264/overview"/>
</contrib>
<contrib contrib-type="author"><name><surname>Goel</surname>
<given-names>Manisha</given-names></name><xref rid="aff2" ref-type="aff"><sup>2</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/31106/overview"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name><surname>Saha</surname>
<given-names>Swati</given-names></name><xref rid="aff1" ref-type="aff"><sup>1</sup></xref>
<xref rid="c001" ref-type="corresp"><sup>&#x002A;</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/2144037/overview"/>
</contrib>
</contrib-group>
<aff id="aff1"><sup>1</sup><institution>Department of Microbiology, University of Delhi South Campus</institution>, <addr-line>New Delhi</addr-line>, <country>India</country></aff>
<aff id="aff2"><sup>2</sup><institution>Department of Biophysics, University of Delhi South Campus</institution>, <addr-line>New Delhi</addr-line>, <country>India</country></aff>
<author-notes>
<fn id="fn0001" fn-type="edited-by">
<p>Edited by: Arthur Charles-Orszag, University of California, San Francisco, United States</p>
</fn>
<fn id="fn0002" fn-type="edited-by">
<p>Reviewed by: Zhenfeng Zhang, Institute of Microbiology (CAS), China; Haruyuki Atomi, Kyoto University, Japan</p>
</fn>
<corresp id="c001">&#x002A;Correspondence: Swati Saha, <email>ss5gp@yahoo.co.in</email>; <email>ssaha@south.du.ac.in</email></corresp>
<fn id="fn0003" fn-type="other">
<p>This article was submitted to Biology of Archaea, a section of the journal Frontiers in Microbiology</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>15</day>
<month>03</month>
<year>2023</year>
</pub-date>
<pub-date pub-type="collection">
<year>2023</year>
</pub-date>
<volume>14</volume>
<elocation-id>1126750</elocation-id>
<history>
<date date-type="received">
<day>18</day>
<month>12</month>
<year>2022</year>
</date>
<date date-type="accepted">
<day>27</day>
<month>02</month>
<year>2023</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x00A9; 2023 Gulati, Singh, Goel and Saha.</copyright-statement>
<copyright-year>2023</copyright-year>
<copyright-holder>Gulati, Singh, Goel and Saha</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<p>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. <italic>Picrophilus torridus</italic> 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 <italic>P. torridus</italic>. 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 <italic>dam</italic> gene has been annotated in the genome sequence. Two other methylases have also been annotated in the <italic>P. torridus</italic> 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 <italic>E.coli</italic>, 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 <italic>in vitro</italic> 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.</p>
</abstract>
<kwd-group>
<kwd>DNA adenine methylase</kwd>
<kwd>Type I restriction-modification system</kwd>
<kwd>modification methylase</kwd>
<kwd>
<italic>Picrophilus torridus</italic>
</kwd>
<kwd>acidophile</kwd>
<kwd>extremophile</kwd>
</kwd-group>
<counts>
<fig-count count="8"/>
<table-count count="0"/>
<equation-count count="0"/>
<ref-count count="57"/>
<page-count count="16"/>
<word-count count="10740"/>
</counts>
</article-meta>
</front>
<body>
<sec id="sec1" sec-type="intro">
<label>1.</label>
<title>Introduction</title>
<p>DNA methylation typically occurs at adenine and cytosine residues. While methylation at adenine occurs exocyclically at the sixth position (N<sup>6</sup>-methyladenine: m6A), methylation at cytosine occurs either exocyclically at the fourth position (N<sup>4</sup>-methylcytosine: m4C), or endocyclically at the fifth position (5-methylcytosine: m5C). The primary DNA methylation mark in eukaryotes is m5C, while in bacteria the m6A methylation mark is the most widely prevalent. The marking of DNA by methylation has two distinctive roles: one, it modulates DNA interactions with proteins, thereby having an impact on DNA replication, DNA mismatch repair and gene expression, and second, it protects the bacterial and archaeal genome from self-destruction by their own restriction systems (<xref ref-type="bibr" rid="ref44">Pingoud et al., 2014</xref>; <xref ref-type="bibr" rid="ref4">Anton and Roberts, 2021</xref>).</p>
<p>The enzymes that mediate DNA methylation in bacteria and archaea use S-adenosyl methionine (AdoMet) as the methyl-group donor and are either orphan methyltransferases or a component of restriction-modification systems. Orphan methyltransferases (which do not have cognate restriction endonucleases) are associated with methylation events that regulate DNA-related transactions like replication and transcription. The Dam methylase which targets adenine methylation in the sequence 5&#x2032;-GATC-3&#x2032; is the best known of these (<xref ref-type="bibr" rid="ref1">Adhikari and Curtis, 2016</xref>). Methylases that are a component of restriction-modification (R-M) systems are designed to protect the genome from being attacked by cognate restriction endonucleases harbored by the organism. R-M systems are linked to the cell&#x2019;s defense against incoming DNA, with the restriction endonuclease (RE) cleaving the foreign DNA that is not methylated at its specific recognition sequence. The host genome remains resistant to cleavage due to the recognition sites on the genome being methylated by the cognate modification methylase (MTase). Such methyltransferases can be a component of three types of R-M systems: Types I-III (<xref ref-type="bibr" rid="ref32">Loenen et al., 2014b</xref>).</p>
<p>The Types I-III R-M systems differ in their subunit structure and organization, recognition sequence and cleavage site, and cofactor requirements. In the widely found Type II R-M systems, the modification methylase and the restriction endonuclease functions reside in independent enzymes that recognize the same palindromic sequence, with the endonuclease cleaving DNA within the (non-methylated) recognition sequence or just beyond. Type I and Type III R-M systems are complex multisubunit enzymes, with the modification methylase and restriction endonuclease functions residing in the same enzyme albeit different subunits (<xref ref-type="bibr" rid="ref48">Rao et al., 2014</xref>; <xref ref-type="bibr" rid="ref31">Loenen et al., 2014a</xref>). Type III R-M systems comprise two subunits: Mod and Res. The Mod is responsible for the recognition of the short asymmetric recognition sequence on DNA as well as its methylation, while the Res mediates DNA cleavage 25&#x2013;27&#x2009;bp downstream of the site. The Type I R-M systems comprise three subunits: S, M, and R, typically existing as pentamers of 1S, 2M, and 2R subunits. They recognize a specific bipartite DNA sequence <italic>via</italic> their S (specificity) subunit, modify the DNA at the recognition sequence <italic>via</italic> their M (modification) subunits, and cleave DNA through their R (restriction) subunits hundreds of base pairs away. While the restriction endonuclease exists in active form only as part of the multisubunit 1S2M2R enzyme, the modification methylase can also exist in active form independent of the R subunits, as a complex with the S subunit which it depends on for DNA sequence recognition (M2S1).</p>
<p>While DNA methylation and its effects have been extensively investigated in bacteria and eukaryotes, it remains a poorly understood area in archaea. The primary methylation marks in archaea are m6A and m4C, with m5C remaining undetected so far in these organisms. Single molecule real time (SMRT) sequencing has led to the analysis of the DNA methylomes of several bacteria and archaea. In a landmark effort, <xref ref-type="bibr" rid="ref11">Blow et al. (2016)</xref> analyzed the DNA methylomes of 230 bacterial and archaeal species. Their study identified over 800 methylated motifs, emphasizing the extent of diversity among these groups of organisms (<xref ref-type="bibr" rid="ref11">Blow et al., 2016</xref>). Other studies include the analysis of the DNA methylome of the hyperthermophilic archaeon <italic>Thermococcus onnurineus</italic> where over 2,000 sites of m6A were detected in addition to several sites of m4C (<xref ref-type="bibr" rid="ref30">Lee et al., 2016</xref>), and the identification of DNA methylation at over 99% of the motifs 5&#x2032;-GATC-3&#x2032; and 5&#x2032;-AGCT-3&#x2032; in the genome of the marine thaumarchaeote <italic>Candidatus Nitrosomarinus catalina</italic> SPOT01 at the adenine and cytosine, respectively (<xref ref-type="bibr" rid="ref2">Ahlgren et al., 2017</xref>). The DNA methylome of the hyperthermoacidophile <italic>Sulfolobus acidocaldarius</italic> was found to carry the m6A mark in the 5&#x2032;-GATC-3&#x2032; motif occurring in two different genome contexts: 5&#x2032;-AGATCC-3&#x2032; and 5&#x2032;-GGATCT/C-3&#x2032; (<xref ref-type="bibr" rid="ref14">Couturier and Lindas, 2018</xref>), and the m4C mark was also detected at the 5&#x2032;-GGCC-3&#x2032; motif that was previously identified as the target site of the SuaI/M.SuaI restriction modification system (<xref ref-type="bibr" rid="ref46">Prangishvili et al., 1985</xref>; <xref ref-type="bibr" rid="ref20">Grogan, 2003</xref>).</p>
<p>Several R-M systems have been annotated in archaeal genomes but very few have been characterized. Among the earliest to be purified and studied were the Type II restriction endonuclease ThaI from <italic>Thermoplasma acidophilum</italic> (<xref ref-type="bibr" rid="ref36">McConnell et al., 1978</xref>), three Type II restriction endonucleases from <italic>Methanococcus aeolicus</italic> (<xref ref-type="bibr" rid="ref51">Schmid et al., 1984</xref>), a Type II R-M system from <italic>Methanobacterium wolfei</italic> (<xref ref-type="bibr" rid="ref33">Lunnen et al., 1989</xref>), three Type II R-M systems from <italic>Methanobacterium thermoformicicum</italic> strains (<xref ref-type="bibr" rid="ref39">Nolling and de Vos, 1992a</xref>,<xref ref-type="bibr" rid="ref40">b</xref>), the PspG1 Type II R-M system from a <italic>Pyrococcus</italic> species (<xref ref-type="bibr" rid="ref38">Morgan et al., 1998</xref>; <xref ref-type="bibr" rid="ref43">Pingoud et al., 2003</xref>) and the SuaI/M.SuaI Type II R-M system in <italic>Sulfolobus acidocaldarius</italic> (<xref ref-type="bibr" rid="ref46">Prangishvili et al., 1985</xref>; <xref ref-type="bibr" rid="ref20">Grogan, 2003</xref>). Other Type II R-M systems that have been studied in archaea are PabI/M.PabI in <italic>Pyrococcus abyssi</italic> (<xref ref-type="bibr" rid="ref21">Ishikawa et al., 2005</xref>; <xref ref-type="bibr" rid="ref55">Watanabe et al., 2006</xref>) a type II restriction endonuclease SuiI from <italic>Sulfolobus islandicus</italic> (<xref ref-type="bibr" rid="ref53">Suzuki and Kurosawa, 2016</xref>) and two Type IIG R-M systems TkoI and TkoII in <italic>Thermococcus kodakarensis</italic> (<xref ref-type="bibr" rid="ref57">Zatopek et al., 2021</xref>). The only Type I R-M system studied in archaea to date is in <italic>Haloferax volcanii</italic> (<xref ref-type="bibr" rid="ref41">Ouellette et al., 2018</xref>).</p>
<p>The present study was initiated with the aim of examining DNA methylation in the thermoacidophilic archaeon <italic>Picrophilus torridus</italic>. Originally isolated from a solfataric field in northern Japan, the <italic>Picrophilus</italic> species are the most acidophilic organisms isolated to date, growing optimally at a pH of 0.7 at 55&#x2013;60&#x00B0;C, with an intracellular pH of ~4.6 (<xref ref-type="bibr" rid="ref50">Schleper et al., 1995</xref>). When the 1.55&#x2009;Mb genome of <italic>Picrophilus torridus</italic> was sequenced (<xref ref-type="bibr" rid="ref17">Futterer et al., 2004</xref>) it was found to be AT-rich (64%) even though 92% of the genome was coding sequence. Our study found that the genome carried m6A methylation, but Dam-mediated methylation was not detectable. The investigation was taken forward with the examination of the modification methylase of the Type I R-M system annotated in its genome sequence, as all Type I R-M systems characterized to date demonstrate m6A methylation activity. The recombinant protein (M.PtoI) was found to mediate DNA adenine methylation under a wide range of conditions <italic>in vitro,</italic> and residues critical to methylation activity were identified.</p>
</sec>
<sec id="sec2" sec-type="materials|methods">
<label>2.</label>
<title>Materials and methods</title>
<sec id="sec3">
<label>2.1.</label>
<title><italic>Picrophilus torridus</italic> cultures and genomic DNA isolation</title>
<p><italic>Picrophilus torridus</italic> (strain DSM 9790) was cultured as described earlier (<xref ref-type="bibr" rid="ref5">Arora et al., 2014</xref>). Briefly, cells were grown at 55&#x00B0;C in liquid culture with aeration, in medium of pH 1.0 comprising yeast extract (0.2%, Difco, United States), glucose (1%, Sigma, United States), potassium dihydrogen phosphate (0.3%), magnesium sulphate (0.05%), calcium chloride (0.025%) and ammonium sulphate (0.02%). For growth analysis, cultures were initiated with 1% (v/v) inoculum in media of different pH, from a starter culture at pH 1.0 that was at OD<sub>600</sub>&#x2009;~&#x2009;1.5. Growth was monitored by measuring OD<sub>600</sub> every 24&#x2009;h over a period of 14&#x2009;days. Three experiments were set up in parallel and the average values of the three are plotted graphically, with error bars representing standard deviation. Genomic DNA was isolated from exponentially growing <italic>Picrophilus torridus</italic> cultures as described earlier (<xref ref-type="bibr" rid="ref5">Arora et al., 2014</xref>). Briefly, harvested cells were resuspended in TEN buffer [20&#x2009;mM Tris-Cl (pH 8.0), 1&#x2009;mM EDTA, 100&#x2009;mM NaCl], lysed by the addition of sodium lauroylsarcosine (1.6%) and Triton X-100 (0.12%), the lysate extracted with phenol, and genomic DNA precipitated using ethanol.</p>
</sec>
<sec id="sec4">
<label>2.2.</label>
<title>Cloning of M.PtoI genes</title>
<p>The ~1.3&#x2009;kb gene encoding the S subunit was amplified off <italic>Picrophilus torridus</italic> genomic DNA using primers PtS-F (5&#x2032;-TAAGATCTCATGAATAAAAAGGATTATAAT-3&#x2032;) and PtS-R (5&#x2032;-TACTCGAGATCACCACTTATTTTCAC-3&#x2032;), with the proof-reading enzyme Phusion DNA polymerase (NEB, United States). The amplicon was cloned into pUC19, followed by subcloning into the BglII-XhoI sites of the pETDuet-1 vector (Novagen, United States), thus obtaining the pET-Duet/PtS clone. The ~1.7&#x2009;kb gene encoding the M subunit was similarly amplified using the primers PtM-F (5&#x2032;-TAGAATTCGATGGCAGATTTAAACTGG-3&#x2032;) and PtM-R (5&#x2032;-TAGCGGCCGCTTC ATAGTAACCAAG-3&#x2032;), the amplicon cloned into pUC19, and then subcloned into the EcoRI-NotI sites of pET-Duet and pET-Duet/PtS, thus obtaining plasmids pET-Duet/PtM and pET-Duet/PtSM, respectively.</p>
</sec>
<sec id="sec5">
<label>2.3.</label>
<title>Creation of PtM mutants</title>
<p>All mutageneses were carried out through PCR using Phusion DNA polymerase. For creating PtM mutants the plasmid pUC/PtM was used as template. The C-terminal deletion mutant was created by PCR using the PtM-F primer in combination with PtM<sub>&#x0394;537-576</sub>-R (5&#x2032;-TAGCGGCCGCTTTATTGTCTATATATAAAG-3&#x2032;). The ~1.6&#x2009;kb amplicon thus obtained was cloned into pJET vector followed by subcloning into the EcoRI-NotI sites of pETDuet/PtS, yielding the plasmid pETDuet/PtSM<sub>&#x0394;537-576</sub>. The NPPW and xxGxxG mutants were created by overlap PCR. To mutate the NPPW motif, the N-terminal part of the gene was amplified using primers PtM-F and PtM-N360A-W363A-R (5&#x2032;-CTGGTTC<underline>GC</underline>TGGAGG<underline>CGCC</underline>G CAAC-3&#x2032;) while the C-terminal part of the gene was amplified using primers PtM-N360A-W363A-F (5&#x2032;-GTTGC<underline>GGCG</underline>CCTCCA<underline>GC</underline>GAACCAG-3&#x2032;) and PtM-R. The full length amplicon obtained by amplification with PtM-F and PtM-R primers was cloned into pJET vector followed by subcloning into the EcoRI-NotI sites of pETDuet/PtS, yielding the plasmid pETDuet/PtSM<sub>N360A-W363A</sub>. To mutate the xxGxxG motif, the N-terminal part of the gene was amplified using primers PtM-F and PtM-G284S-R (5&#x2032;-CCTGCAGTT<underline>GA</underline>ACATGCCGG-3&#x2032;) while the C-terminal part of the gene was amplified using primers PtM-G284S-F (5&#x2032;-CCGGCATGT<underline>TC</underline>AACTGCAGG-3&#x2032;) and PtM-R. The full length amplicon obtained using primers PtM-F and PtM-R was cloned into pJET vector followed by subcloning into the EcoRI-NotI sites of pETDuet/PtS, yielding the plasmid pETDuet/PtSM<sub>G284S</sub>.</p>
</sec>
<sec id="sec6">
<label>2.4.</label>
<title>Creation of PtS mutant</title>
<p>M.PtoI/S<sub>211-222(A-G)</sub> was created using Overlap PCR. The N-terminal fragment of the mutant amplicon was produced using primers PtS-F and PtS<sub>211-222(A-G)</sub>-R (5&#x2032;-GGCCCCGGCCC CAGCTCCACCTGCCCCTGCTCCTGCCAACTGTATTTCCTT-3&#x2032;) and the C-terminal fragment of the mutant amplicon was produced using primers PtS<sub>211-222(A-G)</sub>-F (5&#x2032;-GCAGGAGCAGGGGCAGGTGGAGCTGGGGCCGGGGCCTTCTGAAGGCATC-3&#x2032;) and PtS-R. The full-length mutated gene was amplified using PtS-F and PtS-R and ligated into the BglII-XhoI sites of pET-Duet vector to create plasmid pET-Duet/PtS<sub>211-222(A-G)</sub>. This was followed by cloning the PtM gene into the EcoRI-NotI sites of pETDuet/PtS<sub>211-222(A-G)</sub> to create plasmid pETDuet/PtS<sub>211-222(A-G)</sub>M.</p>
</sec>
<sec id="sec7">
<label>2.5.</label>
<title>Expression and purification of recombinant M.PtoI</title>
<p>The M.PtoI enzymes were expressed in <italic>E.coli</italic> BL21 Codon Plus cells (Stratagene, United States) that had been transformed with their pET-Duet clones. Expression of the wild type or mutant proteins was induced in cells growing exponentially (OD<sub>600</sub>&#x2009;=&#x2009;0.6&#x2013;0.8) at 37&#x00B0;C using 1&#x2009;mM IPTG. Cells were allowed to grow for a further 3&#x2009;h at 37&#x00B0;C before harvesting by centrifugation and resuspension in lysis buffer (100&#x2009;mM Tris.Cl pH8, 150&#x2009;mM NaCl, 10% glycerol) containing protease inhibitors and lysozyme. This was followed by incubation of the cell suspension on ice for 20&#x2009;min and lysis by sonication. After clarification of the lysates by high speed centrifugation they were subjected to cobalt-affinity chromatography using TALON metal affinity resin (Clonetech). The lysates were loaded on the column resin, the column extensively washed (with 100&#x2009;mM Tris.Cl (pH 8), 1&#x2009;M NaCl, 10% glycerol), and the bound proteins eluted using 100&#x2009;mM Tris.Cl (pH 8), 250&#x2009;mM imidazole, 10% glycerol. The M subunit of M.PtoI was similarly purified from bacteria transformed with pET-Duet/PtM.</p>
</sec>
<sec id="sec8">
<label>2.6.</label>
<title><italic>Picrophilus torridus</italic> genomic DNA methylation analysis</title>
<p>Dot blot assay was used to analyze genomic DNA that was isolated from <italic>P. torridus</italic> cultures as described above. For this, 1&#x2009;&#x03BC;g genomic DNA was digested with EcoRI and then denatured with 0.3&#x2009;N sodium hydroxide for 30&#x2009;min. This was followed by neutralization with 1&#x2009;M ammonium acetate before spotting on Hybond N<sup>+</sup> nylon membrane (GE Healthcare, United States) using the Bio-Dot apparatus (Bio-Rad Laboratories, United States). The membrane was baked at 80&#x00B0;C for 2&#x2009;h before probing with anti-m6A antibody (Sigma Aldrich, Cat. No. ABE572) or anti-m5C antibody (EpiGentek, Cat.No. A-1014-010). This was done by blocking the baked dot blots with 10% skim milk (Difco Laboratories) for an hour at room temperature, washing the blots twice with 1X PBS-T, and then incubating them overnight at 4&#x00B0;C with the anti-m6A or anti-m5C antibodies (1:1000 dil in 1X PBS). Following three washes with 1X PBS-T the blots were incubated with HRP-labelled secondary antibody (1:10000 dil, Jackson Laboratory, United States) for an hour and developed using a chemiluminescence method.</p>
</sec>
<sec id="sec9">
<label>2.7.</label>
<title>DNA methylation assay</title>
<p>DNA methylation assays were performed using a dot blot-based assay with lambda DNA (as the substrate) that had been isolated from phage particles hosted by an <italic>E.coli dam<sup>&#x2212;</sup>dcm<sup>&#x2212;</sup></italic> strain (DNA purchased from Thermo, United States; Cat no: SD0021). Reactions (10&#x2009;&#x03BC;l) were typically carried out in a Tris-acetate based buffer [30&#x2009;mM Tris-acetate (pH 5.0), 10&#x2009;mM potassium acetate, 10&#x2009;mM magnesium acetate, 100&#x2009;&#x03BC;g/ml BSA] containing 100&#x2009;&#x03BC;M&#x2009;S-adenosyl methionine (AdoMet; purchased from NEB, United States), 12&#x2009;mM 2-mercaptoethanol, 20&#x2009;mM sodium chloride, 1&#x2009;&#x03BC;g DNA and 1.4&#x2009;&#x03BC;M&#x2009;M.PtoI, at 55&#x00B0;C for 1&#x2009;h. Reactions were stopped with 100&#x2009;mM EDTA. Reactions were analyzed by dot blot assay as described above. The blots were quantified using ImageJ analysis, and the ratios of m6A with reference to <italic>E.coli</italic> genomic DNA spotted on the same blot was determined. Each experiment was done thrice. The graphs depict average values of three experiments with error bars representing standard deviation.</p>
</sec>
<sec id="sec10">
<label>2.8.</label>
<title>CD spectroscopy analysis</title>
<p>CD spectroscopy was carried out as earlier (<xref ref-type="bibr" rid="ref5">Arora et al., 2014</xref>). Briefly, proteins (100&#x2009;&#x03BC;g in 10&#x2009;mM potassium phosphate buffer (pH 5.8)) were analyzed with a Jasco J-815 spectropolarimeter using cuvette of 0.1&#x2009;cm path length. CD spectra were recorded over 250&#x2013;190&#x2009;nm at room temperature, in 1&#x2009;nm steps at a scan speed of 200&#x2009;nm/min. The presented spectra are the average of 20 scans. The data is depicted as the mean residue ellipticity (MRE).</p>
</sec>
<sec id="sec11">
<label>2.9.</label>
<title>Homology modelling and structural analysis</title>
<p>The structures of the M and S subunits of M.PtoI were predicted using the Phyre2 (<italic>p</italic>rotein <italic>h</italic>omology/analog<italic>y r</italic>ecognition <italic>e</italic>ngine version 2.0) online server (<xref ref-type="bibr" rid="ref24">Kelley et al., 2015</xref>).<xref rid="fn0004" ref-type="fn">
<sup>1</sup></xref> The reliability of the predicted models were evaluated using the QMean and LG scores obtained (<xref ref-type="bibr" rid="ref8">Benkert et al., 2009</xref>).<xref rid="fn0005" ref-type="fn">
<sup>2</sup></xref> In both cases the models that were considered demonstrated high coverage against their respective templates. All structures were visualized and illustrations created using PyMOL<xref rid="fn0006" ref-type="fn">
<sup>3</sup></xref> (<xref ref-type="bibr" rid="ref37">Mooers, 2016</xref>).</p>
</sec>
</sec>
<sec id="sec12" sec-type="results">
<label>3.</label>
<title>Results</title>
<sec id="sec13">
<label>3.1.</label>
<title><italic>Picrophilus torridus</italic> genome carries methylated adenine residues</title>
<p>The study was initiated with examining the growth pattern of the organism in media of different pH. The data obtained revealed that while <italic>P. torridus</italic> grew comparably at pH 0.7 and 1.0, growth was considerably slower at pH 2.0 and no growth was discernible at pH 3&#x2013;5 (<xref rid="fig1" ref-type="fig">Figure 1A</xref>). To determine if DNA methylation exists in this microorganism under the extreme conditions in which it thrives, genomic DNA was isolated from <italic>P. torridus</italic> cells grown in media whose pH ranged over 0.7 to 2 and the DNA analyzed by dot blot assay as described above. The results obtained clearly demonstrated that the <italic>Picrophilus torridus</italic> genome carries methylated adenine residues (m6A) regardless of the pH at which it is growing, although displaying lower levels of methylation at pH 0.7 (<xref rid="fig1" ref-type="fig">Figure 1B</xref>). No m5C methylation was detectable (<xref rid="fig1" ref-type="fig">Figure 1C</xref>), and we were unable to assess m4C methylation due to lack of quality antibodies.</p>
<fig position="float" id="fig1">
<label>Figure 1</label>
<caption>
<p><bold>(A)</bold> Analysis of growth of <italic>P. torridus</italic> at different pH. Cultures were initiated with 1% (v/v) inoculum from a starter culture at pH 1.0 that was at OD<sub>600</sub>&#x2009;~&#x2009;1.5. Growth was monitored every 24&#x2009;h. Graph represents average values of three experiments and error bars represent standard deviation. Left panel: data plotted against linear y-axis. Right panel: data plotted against log y-axis. <bold>(B)</bold> Dot blot analysis of <italic>P. torridus</italic> genomic DNA for m6A mark. 1&#x2009;&#x03BC;g of genomic DNA was used in each case. Upper left panel: blot probed with anti-m6A antibody (1:1000 dil). Loading control: agarose gel analysis of input genomic DNA digested with EcoRI. Upper right panel: loading scheme of dot blot. Lower panel: quantification of <italic>P. torridus</italic> genomic DNA signal as a ratio of <italic>E.coli</italic> genomic DNA signal, done using ImageJ software. Bars represent average of three experiments. Error bars represent standard deviation. Student&#x2019;s two-tailed <italic>t-</italic>test was applied, &#x002A;&#x002A;&#x002A;<italic>p</italic>&#x2009;&#x003C;&#x2009;0.0005. <bold>(C)</bold> Dot blot analysis of <italic>P. torridus</italic> genomic DNA for m5C mark. 1&#x2009;&#x03BC;g of genomic DNA was used in each case. Upper panel: blot probed with anti-m5C antibody (1:1000 dil). Lower panel: loading scheme of dot blot. <bold>(D)</bold> Analysis of 5&#x2032;-GA<sup>m6</sup>TC in <italic>P. torridus</italic> genomic DNA. In each case, 1&#x2009;&#x03BC;g of genomic DNA was treated with restriction enzyme. M: molecular weight marker. <bold>(E)</bold> Line diagram of <italic>P. torridus</italic> Type I RM system: PTO0076: gene for R subunit, PTO0077: gene for S subunit, PTO0078: gene for M subunit.</p>
</caption>
<graphic xlink:href="fmicb-14-1126750-g001.tif"/>
</fig>
<p>A previous study by <xref ref-type="bibr" rid="ref27">Koike et al. (2005)</xref> had established the presence of the <italic>dam</italic> gene in several archaea including <italic>P. torridus</italic> by computational analysis, but Dam-mediated adenine methylation was not experimentally verified in <italic>P. torridus</italic>. To check the methylation status of the <italic>P. torridus</italic> genome at Dam target sites (5&#x2032;-GATC-3&#x2032;) we exploited the property of the DpnI restriction enzyme that cleaves methylated but not unmethylated GATC sites. Genomic DNA isolated from <italic>Picrophilus</italic> cells grown in medium of pH 1.0 was thus subjected to DpnI digestion, and it was observed that while <italic>E.coli</italic> genomic DNA was cleaved by the enzyme the <italic>P. torridus</italic> genomic DNA was not (<xref rid="fig1" ref-type="fig">Figure 1D</xref>), nor was <italic>S.pombe</italic> genomic DNA (which lacks m6A methylation as evident from <xref rid="fig1" ref-type="fig">Figure 1B</xref>). These data indicate that while the <italic>P. torridus</italic> genome carries the m6A mark the organism does not harbor Dam-mediated adenine methylation.</p>
<p>An examination of the <italic>P. torridus</italic> genome sequence (KEGG database)<xref rid="fn0007" ref-type="fn">
<sup>4</sup></xref> revealed that in addition to the <italic>dam</italic> gene, the euryarchaeon has genes encoding a Type II R-M system and a Type I R-M system. All Type I modification methylases characterized to date are m6A methylases, and thus we investigated the modification methylase of the <italic>P. torridus</italic> Type I R-M system. The three genes of the <italic>P. torridus</italic> Type I R-M system encoding for R, M and S subunits lie in a single operon on the lower strand of the circular genome (<xref rid="fig1" ref-type="fig">Figure 1E</xref>). The investigation of the modification methylase (henceforth referred to as M.PtoI) commenced with cloning the genes encoding M and S subunits.</p>
</sec>
<sec id="sec14">
<label>3.2.</label>
<title>M.PtoI carries all the motifs that typify Type I modification methylases</title>
<p>The M.PtoI genes (designated herewith as PtM and PtS for M and S subunit genes respectively) were cloned by amplification using primers designed against the published PTO0078 and PTO0077 sequences (encoding M and S subunits respectively), as described in Materials and Methods. The pUC clones obtained were sequenced for verification. The sequences of the M and S subunits of M.PtoI were compared with those of typical Type I modification methylases EcoKI and EcoR124I as well as the identified archaeal Hvo MTase (HVO_2270&#x2013;71), using Clustal Omega and BLASTp analyses (<xref ref-type="bibr" rid="ref52">Sievers et al., 2011</xref>).<xref rid="fn0008" ref-type="fn">
<sup>5</sup></xref> While the M subunit shared ~25&#x2013;28% sequence identity over a coverage of 55&#x2013;85% with the M subunits of the other enzymes (<xref rid="fig2" ref-type="fig">Figure 2</xref>), the S subunit shared between 23&#x2013;25% identity over a coverage of 60&#x2013;93% with their corresponding S subunits (<xref ref-type="supplementary-material" rid="SM1">Supplementary Figure S1</xref>). A comparison of the sequence of M subunit with the sequences of modification methylases of typical Type II systems using the same tools revealed that although all of them harbor the characteristic motifs discussed below (<xref ref-type="supplementary-material" rid="SM1">Supplementary Figure S2</xref>), overall there was no significant sequence identity (data not shown), reflecting the differences in the structure and subunit composition of Type I and Type II modification methylases.</p>
<fig position="float" id="fig2">
<label>Figure 2</label>
<caption>
<p>Comparison of derived amino acid sequence of M subunit of M.PtoI with the sequences of M subunits of other Type I modification methylases: Clustal Omega analysis viewed using Jalview multiple alignment editor (<xref ref-type="bibr" rid="ref52">Sievers et al., 2011</xref>). The conserved motifs are demarcated in black rectangular boxes. Colors are indicative of the physicochemical properties of the amino acids. Pink, aliphatic/hydrophobic; orange/ochre, aromatic; purple, glycine/proline; dark blue, basic; green, hydrophilic; red, acidic; yellow, cysteine.</p>
</caption>
<graphic xlink:href="fmicb-14-1126750-g002.tif"/>
</fig>
<p>m6A and m4C methylases are typified by the presence of nine conserved motifs, identified more than 25&#x2009;years ago by analyzing the sequences of over 40 DNA MTases (<xref ref-type="bibr" rid="ref34">Malone et al., 1995</xref>). By examining these motifs in the context of published structures of MTases it is now known that Motifs I-III and X are involved in the binding of the methyl group donor S-adenosyl methionine (AdoMet). It has been suggested that Motif I and motif X form a pocket that lodges the methionine group of AdoMet while motifs II and III are involved in interactions with the ribose and adenine moieties of the same. On the basis of structural as well as biochemical analyses it is also understood that Motifs IV-VIII are part of the active site pocket and are involved in catalysis of the methyl-transfer reaction, while the Target Recognition Domain (TRD) is responsible for recognition of the specific DNA (<xref ref-type="bibr" rid="ref34">Malone et al., 1995</xref>; <xref ref-type="bibr" rid="ref9">Bheemanaik et al., 2006</xref>). In Type I enzymes where DNA recognition is not a part of the functions of the M subunit the TRD lies in the S subunit. The recognition sequence of Type I enzymes is bipartite, and accordingly, the S subunits harbor two TRD domains. The MTases have been classified into six categories: &#x03B1;, &#x03B2;, &#x03B3;, &#x03B4;, &#x03B5;, and &#x03B6; based on the order of these signature motifs in their primary sequence (<xref ref-type="bibr" rid="ref13">Bujnicki, 2002</xref>). The majority of the m6A and m4C methylases fall in the &#x03B1;, &#x03B2;, and &#x03B3; categories. All the Type I modification methylases as well as several Type II methylases are &#x03B3; MTases, whose order of the motifs is X-I-II-III-IV-V-VI-VII-VIII, and M.PtoI conforms to this as seen in <xref rid="fig2" ref-type="fig">Figure 2</xref>. No clearcut motif II and motif III were identifiable in M.PtoI.</p>
</sec>
<sec id="sec15">
<label>3.3.</label>
<title>The structures of the M and S subunits of M.PtoI are broadly conserved with those of other Type I modification methylases</title>
<p>While the structures of the individual subunits of Type I R-M systems have been solved in case of several enzymes, the holoenzyme (M2S1) structures of only two Type I modification methylases are available to date. The EcoKI modification methylase structures (pdb IDs: 2Y7C and 2Y7H) were proposed based on single particle electron microscopy constructions at 18&#x2009;&#x00C5; resolution (<xref ref-type="bibr" rid="ref25">Kennaway et al., 2009</xref>). The more recent EcoR124I modification methylase and restriction enzyme structures were proposed based on single particle cryo-EM at 4.54&#x2009;&#x00C5; resolution (<xref ref-type="bibr" rid="ref18">Gao et al., 2020</xref>). The basic architecture of the M and S subunits of the two MTases and the interfaces of interaction among the different subunits are largely conserved between the two enzymes. The structures of the M and S subunits of M.PtoI were predicted by modelling using Phyre2. The modelled structures were viewed using PyMOL, which allowed us to superimpose the obtained structures on their respective templates to analyze structural alignments.</p>
<p>The most reliable model obtained for the M subunit used the EcoR124I enzyme as template (pdb ID: 7BTQ; <xref ref-type="bibr" rid="ref18">Gao et al., 2020</xref>). The predicted structure of the M subunit displayed 27% identity with the structure of the M subunits in the 7BTQ cryo-EM structure (4.54&#x2009;&#x00C5;), over a coverage of ~83% (residues 89&#x2013;153 and 164&#x2013;569). As seen in <xref rid="fig3" ref-type="fig">Figure 3</xref>, seven motifs that were aligned in the primary sequence (<xref rid="fig2" ref-type="fig">Figure 2</xref>) were also aligned structurally. The M subunits of Type I RM systems are typified by a C-terminal helix that is involved in interactions with the S subunit, and the M subunit of M.PtoI too harbored this conserved helical C-terminal domain (<xref rid="fig3" ref-type="fig">Figure 3</xref>).</p>
<fig position="float" id="fig3">
<label>Figure 3</label>
<caption>
<p>Structural analysis of M subunit of M.PtoI: Upper panel: Alignment of sequence of M subunit of M.PtoI with M subunit of EcoR124I, using Phyre2. The boxes demarcate the positions of the motifs conserved between the two subunits. CTD: C-terminal domain. Lower left panel: Ribbon representation derived from cryo-EM structure (4.54&#x2009;&#x00C5;) of EcoR124I [pdb ID: 7BTQ reported by <xref ref-type="bibr" rid="ref18">Gao et al. (2020)</xref>]. Lower right panel: Ribbon representation of 3D model of M subunit of M.PtoI, modelled against the M subunit of EcoR124I (7BTQ) as template using Phyre2. CTD: C-terminal domain. Roman numerals indicate motif positions. Arabic numerals depict amino acid residues that were excluded from the models.</p>
</caption>
<graphic xlink:href="fmicb-14-1126750-g003.tif"/>
</fig>
<p>The most reliable structure obtained for the S subunit used the S subunit of the <italic>Methanocaldococcus jannaschii</italic> Type I R-M system (pdb ID: 1YF2; <xref ref-type="bibr" rid="ref26">Kim et al., 2005</xref>) as template. The S subunit of M.PtoI showed 30% identity with the 1YF2 X-ray crystal structure of the <italic>M. jannaschii</italic> S subunit (2.4&#x2009;&#x00C5;), over ~97% coverage (amino acids 6&#x2013;106 and 121&#x2013;437). The S subunit of Type I systems harbors two TRDs (target recognition domains) that are oriented inverted with respect to each other. Each TRD carries a globular domain and an alpha-helical dimerization domain (<xref ref-type="bibr" rid="ref31">Loenen et al., 2014a</xref>). The globular domains (TRDI and TRDII in <xref rid="fig4" ref-type="fig">Figure 4</xref>) are responsible for DNA binding, with the N-terminal one binding to the 5&#x2032; segment (specific sequence) of the bipartite recognition site and the C-terminal one binding to the 3&#x2032; segment (specific sequence) of the same. The alpha-helical domains of the two TRDs (named central conserved region or CCR, and distal conserved region or DCR; see <xref rid="fig4" ref-type="fig">Figure 4</xref>) associate with each other such that the two globular domains are held apart at a fixed distance, thus fixing the length of the non-specific spacer sequence in the bipartite recognition site. The architecture of the S subunit of M.PtoI was predicted to be similar, with two TRDs each harboring a globular domain and an alpha helical domain (<xref rid="fig4" ref-type="fig">Figure 4</xref>).</p>
<fig position="float" id="fig4">
<label>Figure 4</label>
<caption>
<p>Structural analysis of S subunit of M.PtoI: Upper panel: Alignment of sequence of S subunit of M.PtoI with S subunit of <italic>Methanocaldococcus jannaschii</italic> Type I R-M system, using Phyre2. The boxes demarcate the positions of the domains conserved between the two subunits. Lower left panel: Ribbon representation derived from X-ray crystallography structure (2.4&#x2009;&#x00C5;) of <italic>M. jannaschii</italic> Type I RM system&#x2019;s S subunit [pdb ID: 1YF2, reported by <xref ref-type="bibr" rid="ref26">Kim et al. (2005)</xref>]. Lower right panel: Ribbon representation of 3D model of S subunit of M.PtoI, modelled against the 1YF2 S subunit structure using Phyre2. Arabic numerals depict amino acid residues that were excluded from the models. TRDI and TRDII: Target Recognition Domain I and Target Recognition Domain II representing the globular domains of the N-terminal and C-terminal TRDs, respectively. CCR: Central Conserved Region, DCR: Distal Conserved Region.</p>
</caption>
<graphic xlink:href="fmicb-14-1126750-g004.tif"/>
</fig>
<p>The S and M subunits of Type I enzymes characteristically interact with each other (either in the M2S1 or in the R2M2S1 complex) through a four-helix bundle that comprises the C-terminal helical domains of the two M subunits (depicted as CTD in <xref rid="fig3" ref-type="fig">Figure 3</xref>) and the two alpha-helical domains of the S subunit (depicted as CCR and DCR in <xref rid="fig4" ref-type="fig">Figure 4</xref>). As these regions are conserved in the M and S subunits of M.PtoI, the structure of the M.PtoI holoenzyme (M2S1) was predicted using template-based docking, <italic>via</italic> superimposition of the predicted 3D structures of M and S of M.PtoI (<xref rid="fig3" ref-type="fig">Figures 3</xref>, <xref rid="fig4" ref-type="fig">4</xref>) on the M.EcoR124I holoenzyme structure (pdb:7BTQ). The M.PtoI docked structure obtained predicted the two M subunits of the M2S1 holoenzyme to interact with each other <italic>via</italic> their N-terminal domains, and the M and S subunits to interact with each other <italic>via</italic> the four-helix bundle (<xref rid="fig5" ref-type="fig">Figure 5</xref>).</p>
<fig position="float" id="fig5">
<label>Figure 5</label>
<caption>
<p>Predicted structure of M.PtoI: Predicted 3D structures of M and S subunits of M.PtoI were superimposed on the 7BTQ M.EcoR124I holoenzyme structure. M.PtoI structure was thus obtained using template-based docking. M: M subunit, CTD: C-terminal domain, TRDI and TRDII: Target Recognition Domain I and II respectively, CCR: Central Conserved Region, DCR: Distal Conserved Region.</p>
</caption>
<graphic xlink:href="fmicb-14-1126750-g005.tif"/>
</fig>
<p>The role of the C-terminal helical domain of M in mediating interactions with the S subunit was experimentally validated by creating a C-terminal deletion mutant of M: M.PtoI/M<sub>&#x0394;537-576</sub>, and analyzing it for loss of interactions with S. For this, the recombinant M.PtoI proteins were purified from <italic>E.coli</italic>. Recombinant M.PtoI (full length) was expressed in <italic>E.coli</italic> using the pET-Duet expression vector which allows the simultaneous expression of two proteins from two separate T7<italic>lac</italic> promoters, with one of the two proteins being expressed in fusion with the His tag at the N-terminus. The M and S subunits were thus simultaneously expressed in BL21 Codon Plus <italic>E. coli</italic> cells from the pETDuet/PtSM clone (as described in Methods) such that M was His-tagged. The cell lysates were subjected to cobalt affinity-based chromatography. SDS-PAGE analysis revealed the eluate to carry both subunits in stoichiometric amounts, leading us to conclude that the M.PtoI holoenzyme (M2S1) assembles in the host bacterium, allowing the S subunit to co-purify with the His-tagged M (<xref rid="fig6" ref-type="fig">Figure 6A</xref>). In creating the C-terminal deletion of M, amino acids 537&#x2013;576 were excluded (mutagenesis detailed in Methods) and the M<sub>&#x0394;537-576</sub> co-expressed with S in BL21 Codon Plus cells from the pETDuet/PtSM<sub>&#x0394;537-576</sub> clone. When these cell lysates were put through TALON beads the eluate carried only the truncated M subunit, clearly demonstrating loss of interaction with S in M.PtoI/M<sub>&#x0394;537-576</sub> (<xref rid="fig6" ref-type="fig">Figure 6B</xref>, left panel, fourth lane). The deletion did not cause any gross structural changes in the M subunit based on a comparison of the CD spectroscopy profiles of purified M (<xref rid="fig6" ref-type="fig">Figure 6B</xref>, right panel, fourth lane) and M.PtoI/M<sub>&#x0394;537-576</sub> (<xref rid="fig6" ref-type="fig">Figure 6C</xref>). Taken together, these results indicate that the C-terminal helical domain of M is functionally conserved in M.PtoI.</p>
<fig position="float" id="fig6">
<label>Figure 6</label>
<caption>
<p><bold>(A,B,D)</bold> Analysis of purified recombinant M.PtoI proteins. Coomassie-staining of SDS-PAGE. Mr.: molecular weight marker, preload: bacterial whole cell lysates after induction of expression of recombinant proteins, flowthrough: unbound protein fraction after incubation with cobalt affinity beads, eluate: protein fraction after elution with 250&#x2009;mM imidazole. Arrowheads indicate subunits in eluate fractions. <bold>(A)</bold> Wild type. Eluate: M and S subunits (molecular weights 63&#x2009;kDa and 48&#x2009;kDa respectively). <bold>(B)</bold> Left panel: M.PtoI/M<sub>&#x0394;537-576</sub>. Eluate: M<sub>&#x0394;537-576</sub> subunit only (molecular weights ~58&#x2009;kDa). Right panel: M subunit of M.PtoI. Eluate: M subunit (molecular weight 63&#x2009;kDa). <bold>(C)</bold> CD spectra of M.PtoI M subunit and M.PtoI/M<sub>&#x0394;537-576</sub> depicted as a measure of mean residue ellipticity (MRE). <bold>(D)</bold> M.PtoI/S<sub>211-222(A-G).</sub> Eluate: M and S subunits (molecular weights 63&#x2009;kDa and 48&#x2009;kDa respectively). <bold>(E)</bold> CD spectra of M.PtoI (left panel) and M.PtoI/S<sub>211-222(A-G)</sub> (right panel) depicted as a measure of mean residue ellipticity (MRE).</p>
</caption>
<graphic xlink:href="fmicb-14-1126750-g006.tif"/>
</fig>
<p>The role of the central helical region of the S subunit in mediating M-S interactions was examined by creating a mutant S subunit where 12 amino acids of the central conserved region (CCR) were replaced by a 12 amino acid stretch of alternate glycine and alanine residues to create M.PtoI/S<sub>211-222(A-G)</sub> (described in &#x201C;Methods&#x201D;). When the recombinant M.PtoI/S<sub>211-222(A-G)</sub> protein was purified from <italic>E.coli</italic> using cobalt affinity chromatography it was found that there was a substantial increase in the ratio of M:S in the eluate fractions (<xref rid="fig6" ref-type="fig">Figure 6D</xref>, fourth lane), implying that the replaced region played a significant role in mediating M-S interactions. The difference in CD spectroscopy profiles of the wild type and mutant M.PtoI proteins most likely reflects the differences in M:S stoichiometry in the two proteins (<xref rid="fig6" ref-type="fig">Figure 6E</xref>).</p>
<p>M-S interactions in wild type and mutant M.PtoI proteins were also checked using purified proteins. For this, M and S subunits were separately expressed in <italic>E.coli</italic> cells from respective pET/Duet clones carrying either M or S genes. As all M subunits were tagged with His at their N-terminus, they were immobilized on Co-affinity resin. Whole cell lysates harboring recombinant S subunit (wild type or mutant) were added to the reaction, and after washing off the unbound fractions, the bound fractions (M subunits plus any interacting S subunits) were eluted using imidazole and analyzed. As seen in <xref ref-type="supplementary-material" rid="SM1">Supplementary Figure S3</xref>, while S subunit interacts with wild type M subunit (lane 5), it fails to interact with mutant M<sub>&#x0394;537-576</sub> subunit (lane 6). The mutant S<sub>211-222(A-G)</sub> subunit interacts with M subunit to a significantly lesser extent than wild type S (compare lane 7 with lane 5). These data confirm the role of the C-terminus of M subunit and residues 211&#x2013;222 of S subunit in maintaining M-S interactions.</p>
</sec>
<sec id="sec16">
<label>3.4.</label>
<title>M.PtoI methylates DNA at adenine residues <italic>in vitro</italic></title>
<p>The ability of the purified M.PtoI to methylate DNA at adenine residues was analyzed using unmethylated lambda DNA isolated from phage that had been grown in a <italic>dam<sup>&#x2212;</sup>dcm<sup>&#x2212;</sup> E.coli</italic> strain. The reactions were analyzed in dot blot assays using anti-m6A antibody, as described in Methods. Activity assays were initiated at 55&#x00B0;C for 1&#x2009;h in Tris-acetate buffers whose pH ranged over 3 to 7, in the presence of sodium chloride ranging over 0&#x2013;100&#x2009;mM. As seen in <xref rid="fig7" ref-type="fig">Figure 7A</xref> M.PtoI was able to methylate lambda DNA (producing m6A) at pH ranging from 4 to 7, with optimum activity being detected at pH 5, which is near the organism&#x2019;s intracellular pH of 4.6. Higher salt concentrations were detrimental to activity, with optimal activity seen at 0&#x2013;20&#x2009;mM NaCl. Thus, all further assays were carried out in Tris-acetate buffer (pH 5) containing 20&#x2009;mM NaCl. The optimal concentration of magnesium ions that supported the DNA methylation reaction was determined over 0&#x2013;50&#x2009;mM magnesium acetate. It was found that while magnesium was necessary for the methylation reaction, concentrations of magnesium over 5&#x2013;25&#x2009;mM supported the reaction more or less equivalently, with a higher concentration of 50&#x2009;mM being inhibitory (<xref rid="fig7" ref-type="fig">Figure 7B</xref>). The methylation reaction displayed an initial lag before proceeding almost linearly for approximately the first 45&#x2009;min and plateaued beyond 60&#x2009;min (<xref rid="fig7" ref-type="fig">Figure 7C</xref>). The ability of M.PtoI to methylate DNA was concentration-dependent, increasing almost linearly between 200 and 600&#x2009;nM before plateauing (<xref rid="fig7" ref-type="fig">Figure 7D</xref>). High concentrations of enzyme (2,000&#x2009;nM and beyond) proved to be detrimental, perhaps due to protein aggregation occurring at such high concentrations. In carrying out the assays at different temperatures ranging over 37&#x2013;70&#x00B0;C it was found that M.PtoI was active at temperatures between 37&#x2013;65&#x00B0;C with maximal activity at 55 and 60&#x00B0;C, in keeping with the conditions in which it thrives (<xref rid="fig7" ref-type="fig">Figure 7E</xref>). AdoMet supported the methylation reaction over concentrations varying from 10&#x2013;500&#x2009;&#x03BC;M, with 50&#x2013;100&#x2009;&#x03BC;M being optimal (<xref rid="fig7" ref-type="fig">Figure 7F</xref>). Interestingly, the enzyme displayed substrate inhibition at higher AdoMet concentration of 1&#x2009;mM. Furthermore, the enzyme was able to methylate DNA, although to a much lesser extent, even in the absence of exogenously added AdoMet, suggesting that endogenous AdoMet remains bound to the M subunit through the purification process. This is not unusual for such enzymes, with similar findings in multiple enzymes such as in BpmI, BseMII, and EcoR124II (<xref ref-type="bibr" rid="ref15">Dreier et al., 1996</xref>; <xref ref-type="bibr" rid="ref22">Jurenaite-Urbanaviciene et al., 2001</xref>; <xref ref-type="bibr" rid="ref7">Bath et al., 2002</xref>). Unsurprisingly, the addition of ATP had no effect on methylation activity (data not shown). Taken together, the results presented in <xref rid="fig7" ref-type="fig">Figure 7</xref> establish M.PtoI as an active m6A MTase.</p>
<fig position="float" id="fig7">
<label>Figure 7</label>
<caption>
<p><bold>(A&#x2013;F)</bold> Dot blot analyses of methylation assays performed with M.PtoI. Assays were performed with 1&#x2009;&#x03BC;g lambda DNA substrate, in Tris-acetate based buffer (pH 5) carrying 10 mM magnesium acetate, 100 &#x03BC;M S-adenosyl methionine (AdoMet), 20&#x2009;mM sodium chloride, and 1.4 &#x03BC;M M.PtoI, at 55&#x00B0;C, unless otherwise specified. 1 &#x03BC;g <italic>E.coli</italic> genomic DNA was used as control on every dot blot. Upper panels in each figure: loading scheme of dot blot. Middle panel: Dot blot of assay. Lower panel: quantification of signals from methylation reactions as a ratio of <italic>E.coli</italic> genomic DNA signal, done using ImageJ software. Data represent average of three experiments. Error bars represent standard deviation. Student&#x2019;s two-tailed <italic>t-</italic>test was applied, &#x002A;<italic>p</italic>&#x2009;&#x003C;&#x2009;0.05, &#x002A;&#x002A;<italic>p</italic>&#x2009;&#x003C;&#x2009;0.005, &#x002A;&#x002A;&#x002A;<italic>p</italic>&#x2009;&#x003C;&#x2009;0.0005, ns-not significant. <bold>(A)</bold> Reaction as a function of pH and salt concentration: Activity was assayed over pH 3&#x2013;7 in presence of 0&#x2013;100 mM sodium chloride, using 1.4 &#x03BC;M M.PtoI at 55&#x00B0;C for 1&#x2009;h. <bold>(B)</bold> Reaction as a function of magnesium ion concentration: Activity was assayed over 0&#x2013;50 mM magnesium acetate, using 1.4 &#x03BC;M M.PtoI at 55&#x00B0;C for 1&#x2009;h. Numbers in the boxes above the blot represent magnesium acetate concentrations (mM). Reaction without protein had 10&#x2009;mM magnesium acetate. <bold>(C)</bold> Progress curve of the reaction: Activity was monitored over 2&#x2013;120&#x2009;min using 1.4 &#x03BC;M M.PtoI, at 55&#x00B0;C. Numbers in the boxes above the blot represent incubation time (min). Reaction without protein was incubated for 120 min. <bold>(D)</bold> Reaction as a function of enzyme concentration: Activity was assayed with 1&#x2013;2,500 nM M.PtoI, at 55&#x00B0;C for 45 min. Numbers in the boxes above the blot represent enzyme concentrations (nM). <bold>(E)</bold> Reaction as a function of temperature: Activity was assayed over 37&#x2013;70&#x00B0; C, using 600 nM M.PtoI for 45&#x2009;min. Numbers in the boxes above the blot represent reaction temperature. Reaction without protein was incubated at 55&#x00B0;C. <bold>(F)</bold> Reaction as a function of AdoMet concentration: Activity was assayed over 0&#x2013;1,000 &#x03BC;M S-adenosyl methionine (AdoMet), using 400 nM M.PtoI at 55&#x00B0;C for 1 h. Numbers in the boxes above the blot represent AdoMet concentrations (&#x03BC;M). Reaction without protein had 100 &#x03BC;M AdoMet.</p>
</caption>
<graphic xlink:href="fmicb-14-1126750-g007.tif"/>
</fig>
</sec>
<sec id="sec17">
<label>3.5.</label>
<title>The NPPW and xxGxxG motifs in M.PtoI are critical for methylation activity</title>
<p>Type I modification methylases, like most adenine MTases, are typified by two primary domains they harbor in their M subunits: the AdoMet binding domain (motifs X, I-III in <xref rid="fig2" ref-type="fig">Figures 2</xref>, <xref rid="fig3" ref-type="fig">3</xref>), and the catalytic domain that is responsible for the transfer of the methyl group from the donor to the nucleotide which is to be methylated (motifs IV-VIII in <xref rid="fig2" ref-type="fig">Figures 2</xref>, <xref rid="fig3" ref-type="fig">3</xref>). Two motifs that are of singular importance are the AdoMet-binding motif I (FxGxxG) and the catalytic motif IV (NPPY/F/W). The catalytic motif IV lies in a pocket on the surface of the M subunit and at the time of methylation the specific adenine residue flips out of the duplex DNA substrate into the pocket where it stacks with the aromatic residue as well as forms hydrogen bonds with the first amino acid of the motif (which could by D, S, or N), stabilizing it for receiving the methyl group which is directly transferred to the target amino group on the adenine (<xref ref-type="bibr" rid="ref9">Bheemanaik et al., 2006</xref>). Based on the sequence alignment with other Type I MTases (<xref rid="fig2" ref-type="fig">Figure 2</xref>), the FxGxxG motif in M is located at residues 282&#x2013;287 and the NPPW motif lies at residues 360&#x2013;363 (<xref rid="fig8" ref-type="fig">Figure 8A</xref>). The importance of these two motifs for the methyltransferase activity of M.PtoI was assessed by developing mutant enzymes and analyzing their activity. Accordingly, the glycine residue at position 284 was mutated to serine, and the asparagine and tryptophan residues at positions 360 and 363 were mutated to alanine (described in Methods). The mutated M subunits were co-expressed with S subunit in <italic>E.coli</italic> cells from pETDuet/PtSM<sub>G284S</sub> and pETDuet/PtSM<sub>N360A-W363A</sub>, respectively. The mutant enzymes M.PtoI/M<sub>G284S</sub> and M.PtoI/M<sub>N360A-W363A</sub> were purified (<xref rid="fig8" ref-type="fig">Figure 8B</xref>) and assessed for any gross structural changes using CD spectroscopy; no gross structural changes were detectable (<xref rid="fig8" ref-type="fig">Figure 8C</xref>).</p>
<fig position="float" id="fig8">
<label>Figure 8</label>
<caption>
<p><bold>(A)</bold> Ribbon representation of motif I and motif IV derived from the predicted structure of M subunit of M.PtoI. <bold>(B)</bold> Coomassie stained SDS-PAGE analysis of purified recombinant M.PtoI/M<sub>N360A-W363A</sub> and M.PtoI/M<sub>G284S.</sub> Mr.: molecular weight marker. Arrowheads depict M and S subunits (63&#x2009;kDa and 48&#x2009;kDa respectively). <bold>(C)</bold> CD spectra of M.PtoI/M<sub>N360A-W363A</sub> and M.PtoI/M<sub>G284S</sub> depicted as a measure of mean residue ellipticity (MRE). <bold>(D,E)</bold> Dot blot analyses of methylation assays with mutant M.PtoI proteins. Assays were performed as described in the legend to <xref rid="fig7" ref-type="fig">Figure 7</xref> unless otherwise specified. Left panel: loading scheme of dot blot. Middle panel: Dot blot of assay. Right panel: quantification of signals from methylation reactions as a ratio of <italic>E.coli</italic> genomic DNA signal, done using ImageJ software. Data represent average of three experiments. Error bars represent standard deviation. Student&#x2019;s two-tailed <italic>t-</italic>test was applied, &#x002A;<italic>p</italic>&#x2009;&#x003C;&#x2009;0.05, &#x002A;&#x002A;<italic>p</italic>&#x2009;&#x003C;&#x2009;0.005, &#x002A;&#x002A;&#x002A;<italic>p</italic>&#x2009;&#x003C;&#x2009;0.0005, ns-not significant. <bold>(D)</bold> Reactions with M.PtoI, M.PtoI/M<sub>N360A-W363A</sub> and M.PtoI/M<sub>G284S</sub>. Activity assayed as a function of sodium chloride concentration (0&#x2013;50&#x2009;mM) at pH 5. Numbers in boxes represent NaCl concentrations (mM). Reaction without protein had 20&#x2009;mM NaCl. <bold>(E)</bold> Reactions with M.PtoI and M.PtoI/M<sub>G284S</sub>. Activity assayed as a function of AdoMet concentration (0&#x2013;1,000&#x2009;&#x03BC;M). Numbers in boxes represent AdoMet concentrations (&#x03BC;M). Reaction without protein had 100&#x2009;&#x03BC;M AdoMet.</p>
</caption>
<graphic xlink:href="fmicb-14-1126750-g008.tif"/>
</fig>
<p>The two mutant enzymes were assessed for their ability to methylate DNA in Tris-acetate buffer at pH 5 over sodium chloride concentrations ranging from 0&#x2013;50&#x2009;mM. M.PtoI/M<sub>N360A-W363A</sub> was found to be virtually inactive under all tested conditions, while M.PtoI/M<sub>G284S</sub> was less active than wild type enzyme at all tested conditions (<xref rid="fig8" ref-type="fig">Figure 8D</xref>). The activity of M.PtoI/M<sub>G284S</sub> was compared with that of the wild type M.PtoI at different concentrations of AdoMet. It was observed that while the wild type enzyme supported DNA methylation more or less equivalently over 50&#x2013;500&#x2009;&#x03BC;M AdoMet, the mutant enzyme while displaying an overall lower activity showed maximal support of DNA methylation at 100&#x2013;200&#x2009;&#x03BC;M AdoMet (<xref rid="fig8" ref-type="fig">Figure 8E</xref>). These data establish that the NPPW and xxGxxG motifs in M.PtoI are critical to M.PtoI activity.</p>
</sec>
</sec>
<sec id="sec18" sec-type="discussions">
<label>4.</label>
<title>Discussion</title>
<p>DNA methylation is widely prevalent in bacteria and archaea, and known to modulate various cellular processes as well as exist as part of the organisms&#x2019; defense mechanisms against foreign DNA. <italic>Picrophilus torridus</italic> grows in conditions of extreme pH, thriving at pH 0.7 to 2 (<xref rid="fig1" ref-type="fig">Figure 1A</xref>), and to date no other microorganism growing under these extreme conditions has been isolated. We initiated the study with examining whether the <italic>Picrophilus</italic> genome is methylated, and found it to carry the m6A modification mark but not the m5C mark (<xref rid="fig1" ref-type="fig">Figures 1B</xref>,<xref rid="fig1" ref-type="fig">C</xref>). Interestingly, the extent of adenine methylation was significantly lower when the organism was grown at pH 0.7 as compared to when grown at pH 1 or 2. The commonest mediator of adenine methylation in prokaryotes is the Dam methylase. The <italic>E. coli</italic> Dam methylase regulates the timing of replication origin firing <italic>via</italic> methylation of the GATC sites in the <italic>OriC</italic> region. Dam-mediated GATC methylation also modulates methyl-directed mismatch repair (<xref ref-type="bibr" rid="ref1">Adhikari and Curtis, 2016</xref>), and has been identified as a regulator of gene expression in <italic>E. coli</italic>, <italic>Yersinia</italic>, and <italic>Salmonella typhimurium</italic> (<xref ref-type="bibr" rid="ref16">Falker et al., 2007</xref>; <xref ref-type="bibr" rid="ref12">Brunet et al., 2020</xref>; <xref ref-type="bibr" rid="ref23">Keceli Oguz et al., 2022</xref>). Widely prevalent across bacteria, some archaea also harbor the Dam. <italic>Halobacterium saccharovorum</italic> and <italic>Methanobacterium</italic> strain Ivanov have been reported to be Dam<sup>+</sup> as early as 1984 (<xref ref-type="bibr" rid="ref6">Barbeyron et al., 1984</xref>), and the <italic>in vitro</italic> methylation activity of the Dam methylase of <italic>Pyrococcus horikoshii</italic> has also been examined (<xref ref-type="bibr" rid="ref35">Maynard-Smith et al., 2011</xref>). <xref ref-type="bibr" rid="ref27">Koike et al. (2005)</xref> identified the Dam methylase in multiple archaea species through annotation of genome sequences. Their experimental analysis found <italic>Thermoplasma volcanium, Thermoplasma acidophilum</italic>, and <italic>Pyrocococcus</italic> species OT3 genomes to carry methylated 5&#x2032;-GATC sequences, while <italic>Sulfolobus solfataricus</italic> and <italic>Sulfolobus shibatae</italic> did not show evidence of 5&#x2032;-GATC methylation. A subsequent study by <xref ref-type="bibr" rid="ref14">Couturier and Lindas (2018)</xref> identified the methylation of 5&#x2032;-GATC sites in <italic>Sulfolobus acidocaldarius</italic>. While the Dam methylase was identified by sequence annotation in <italic>P. torridus</italic> in the study by <xref ref-type="bibr" rid="ref27">Koike et al. (2005)</xref>, 5&#x2032;-GATC methylation was not experimentally verified in that study. Our results indicated lack of 5&#x2032;-GATC methylation in the <italic>P. torridus</italic> genome (<xref rid="fig1" ref-type="fig">Figure 1D</xref>), suggesting the absence of an active Dam in this organism and implicating the possible role of an R-M system in mediating adenine methylation in the organism.</p>
<p>The genes encoding the components of the single Type I R-M system identified through sequence annotation of the <italic>P. torridus</italic> genome lie in a cluster, on the lower strand of the genome between positions 85,223 and 91,251 (<xref rid="fig1" ref-type="fig">Figure 1E</xref>), and the modification methylase component of this system, named M.PtoI, was investigated for its ability to methylate adenine residues. The structure of M.PtoI and the interacting interfaces of the M and S subunits were predicted to be conserved with those of other Type I R-M systems (<xref rid="fig3" ref-type="fig">Figures 3</xref>&#x2013;<xref rid="fig5" ref-type="fig">5</xref>), and the predictions were experimentally verified in part by the creation and analysis of suitable M and S mutant proteins (<xref rid="fig6" ref-type="fig">Figure 6</xref>), which confirmed that the M subunits interacted with the central conserved region (CCR) of the S subunit through their C-terminal helical domains (CTDs). The ability of M.PtoI to methylate DNA was assessed using the recombinant protein in <italic>in vitro</italic> assays. Activity was found to be optimal at pH and temperature conditions that reflected the organism&#x2019;s lifestyle and physiological state (<xref rid="fig7" ref-type="fig">Figures 7A</xref>,<xref rid="fig7" ref-type="fig">E</xref>). The M.PabI Type II modification methylase of <italic>Pyrococcus abyssi</italic> has been found to methylate DNA at 95&#x00B0;C as much as at 65&#x00B0;C although optimally at 85&#x00B0;C, and demonstrated optimal activity at pH 5.8 to 6.7, similarly reflecting the growth and intracellular conditions of the hyperthermophilic organism (<xref ref-type="bibr" rid="ref55">Watanabe et al., 2006</xref>). Magnesium ions were found to be essential for M.PtoI activity (<xref rid="fig7" ref-type="fig">Figure 7B</xref>). Magnesium ions (but not calcium or manganese ions) have been found to be essential for M.EcoP15I methyltransferase activity as well. CD spectroscopy analysis of M.EcoP15I titrated against increasing concentrations of magnesium revealed a secondary structure alteration in the enzyme in response to magnesium ions, suggesting that magnesium ion-induced conformational changes preceded the catalysis reaction. Mutation of the magnesium binding motif resulted in enzyme inactivity (<xref ref-type="bibr" rid="ref10">Bist and Rao, 2003</xref>). Earlier studies with EcoBI as well as MmeI methyltransferases indicated that the presence of magnesium ions, though not essential to the methylation reaction, stimulated it (<xref ref-type="bibr" rid="ref29">Lautenberger and Linn, 1972</xref>; <xref ref-type="bibr" rid="ref54">Tucholski et al., 1995</xref>). M.PtoI was found to be active over AdoMet concentrations 25&#x2013;500&#x2009;&#x03BC;M, exhibiting substrate inhibition at 1&#x2009;mM AdoMet. The purified protein most likely harbored endogenous AdoMet as it exhibited methylation activity even when AdoMet was not added to the reaction (<xref rid="fig7" ref-type="fig">Figure 7F</xref>), a feature of a few other methyltransferases as well (<xref ref-type="bibr" rid="ref15">Dreier et al., 1996</xref>; <xref ref-type="bibr" rid="ref22">Jurenaite-Urbanaviciene et al., 2001</xref>; <xref ref-type="bibr" rid="ref7">Bath et al., 2002</xref>).</p>
<p>M.PtoI harbored all the motifs typifying adenine methyltransferases (<xref rid="fig2" ref-type="fig">Figure 2</xref>). The functional roles of these motifs have been uncovered by a combination of structural and biochemical studies of other DNA methyltransferases. The methylation reaction is mediated by a base-flipping mechanism, wherein the target adenine residue flips out of the double helix into the catalytic pocket (<xref ref-type="bibr" rid="ref19">Goedecke et al., 2001</xref>). Motif IV (NPPW) is a primary constituent of the catalytic pocket, and the interaction of the flipped target adenine with the asparagine residue is crucial for catalysis. The data from structural and biochemical studies lead us to believe that hydrogen bonding between the N<sup>6</sup>-amino group of the adenine (which serves as the donor) and the asparagine side chain of the NPPW motif possibly polarizes the amino group, promoting the direct transfer of the methyl group from AdoMet to the adenine N<sup>6</sup> position (<xref ref-type="bibr" rid="ref45">Pogolotti et al., 1988</xref>; <xref ref-type="bibr" rid="ref28">Labahn et al., 1994</xref>; <xref ref-type="bibr" rid="ref19">Goedecke et al., 2001</xref>; <xref ref-type="bibr" rid="ref9">Bheemanaik et al., 2006</xref>). The tryptophan residue at the last position of the motif can be replaced with any other aromatic residue without any significant impact on catalytic activity, but replacement with other residues has a negative impact on activity (<xref ref-type="bibr" rid="ref56">Willcock et al., 1994</xref>; <xref ref-type="bibr" rid="ref3">Ahmad et al., 1995</xref>; <xref ref-type="bibr" rid="ref47">Pues et al., 1999</xref>). In keeping with these findings, the ability of the M.PtoI/M<sub>N360A-W363A</sub> motif IV mutant to methylate DNA was severely compromised (<xref rid="fig8" ref-type="fig">Figure 8D</xref>). Structural studies of AdoMet-bound enzyme as well as in-depth mutational analyses coupled to biochemical assays have together revealed that motif I plays a major role in AdoMet binding, though it is not directly involved in catalysis <italic>per se</italic> (<xref ref-type="bibr" rid="ref28">Labahn et al., 1994</xref>; <xref ref-type="bibr" rid="ref56">Willcock et al., 1994</xref>; <xref ref-type="bibr" rid="ref3">Ahmad et al., 1995</xref>; <xref ref-type="bibr" rid="ref49">Roth et al., 1998</xref>; <xref ref-type="bibr" rid="ref19">Goedecke et al., 2001</xref>). It is possible that the M.PtoI/M<sub>G284A</sub> mutant displays weaker binding of AdoMet, as the mutant (unlike the wild type enzyme) does not appear to be harboring endogenous AdoMet when purified from <italic>E.coli</italic> since no catalysis is detected at 0&#x2009;&#x03BC;M AdoMet (<xref rid="fig8" ref-type="fig">Figure 8E</xref>). Collectively, the data suggest that the functions of motif I and motif IV may be conserved in M.PtoI.</p>
<p>To date, only one Type I R-M system has been investigated in archaea: that of <italic>Haloferax volcanii.</italic> While several MTases have been annotated in the <italic>Haloferax</italic> genome sequence, a combination of single molecule real time sequencing and gene knockout analyses have demonstrated that only two motifs are methylated in the genome, one of which is methylated at a cytosine residue <italic>via</italic> a Type II enzyme and the second of which (5&#x2032;-GCAm6BN6VTGC-3&#x2032;) is methylated by a Type I adenine modification methylase (<xref ref-type="bibr" rid="ref42">Ouellette et al., 2015</xref>, <xref ref-type="bibr" rid="ref41">2018</xref>). The Type I modification methylase harbors all the motifs that typify these enzymes. Interestingly and unusually, adenine methylation was detected only in the first half of the bipartite recognition sequence, suggesting that either the second half of the sequence does not exhibit methylation on the lower strand, or methylation does not occur equivalently on the two strands and thus remains undetected during sequencing.</p>
<p>Taken together, the data presented here establish that even organisms growing in such extreme conditions (as <italic>Picrophilus torridus</italic> does) exhibit DNA methylation. We find evidence for m6A but not m5C methylation in <italic>P. torridus,</italic> and the absence of GATC methylation indicates that the organism lacks an active Dam. Our results suggest that the Type I R-M system annotated in the <italic>P. torridus</italic> genome sequence is active and mediates m6A methylation. We were unable to analyze the <italic>in vivo</italic> role of M.PtoI as <italic>Picrophilus</italic> species are not amenable to genetic manipulations. Attempts to carry out complementation analyses in <italic>E.coli dam</italic><sup>&#x2212;</sup> strain (a kind gift from the Lab Collection of Manjula Reddy at CSIR-CCMB, India) did not succeed as the M and S subunits did not get expressed (from the tetracycline-driven promoter system in pASKIBA43PLUS) in this strain, which is deficient in <italic>argU</italic>, <italic>ileY</italic>, and <italic>leuW</italic> tRNA genes (<xref ref-type="supplementary-material" rid="SM1">Supplementary Figure S4</xref>). Future studies will be directed toward identifying the target site of M.PtoI and characterizing the PtoI restriction enzyme (M2R2S1).</p>
</sec>
<sec id="sec19" sec-type="data-availability">
<title>Data availability statement</title>
<p>The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.</p>
</sec>
<sec id="sec20">
<title>Author contributions</title>
<p>SS and MG designed the project. PG and AS performed the research and prepared the figures. PG, AS, MG, and SS analyzed the data. SS wrote the paper. All authors contributed to the article and approved the submitted version.</p>
</sec>
<sec id="conf1" sec-type="COI-statement">
<title>Conflict of interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec id="sec100" sec-type="disclaimer">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
</body>
<back>
<ack>
<p>We thank the Central Instrumentation Facility, University of Delhi South Campus for DNA sequencing and the use of CD spectrometer. We thank Manjula Reddy (CSIR-CCMB, Hyderabad, India) for providing us with an <italic>E.coli dam</italic><sup>&#x2212;</sup> strain from her Lab Collection.</p>
</ack>
<sec id="sec22" sec-type="supplementary-material">
<title>Supplementary material</title>
<p>The Supplementary material for this article can be found online at: <ext-link xlink:href="https://www.frontiersin.org/articles/10.3389/fmicb.2023.1126750/full#supplementary-material" ext-link-type="uri">https://www.frontiersin.org/articles/10.3389/fmicb.2023.1126750/full#supplementary-material</ext-link></p>
<supplementary-material xlink:href="Data_Sheet_1.pdf" id="SM1" mimetype="application/pdf" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</sec>
<ref-list>
<title>References</title>
<ref id="ref1">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Adhikari</surname> <given-names>S.</given-names></name> <name><surname>Curtis</surname> <given-names>P. D.</given-names></name></person-group> (<year>2016</year>). <article-title>DNA methyltransferases and epigenetic regulation in bacteria</article-title>. <source>FEMS Microbiol. Rev.</source> <volume>40</volume>, <fpage>575</fpage>&#x2013;<lpage>591</lpage>. doi: <pub-id pub-id-type="doi">10.1093/femsre/fuw023</pub-id></citation>
</ref>
<ref id="ref2">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ahlgren</surname> <given-names>N. A.</given-names></name> <name><surname>Chen</surname> <given-names>Y.</given-names></name> <name><surname>Needham</surname> <given-names>D. M.</given-names></name> <name><surname>Parada</surname> <given-names>A. E.</given-names></name> <name><surname>Sachdeva</surname> <given-names>R.</given-names></name> <name><surname>Trinh</surname> <given-names>V.</given-names></name> <etal/></person-group>. (<year>2017</year>). <article-title>Genome and epigenome of a novel marine Thaumarchaeota strain suggest viral infection, phosphorothioation DNA modification and multiple restriction systems</article-title>. <source>Environ. Microbiol.</source> <volume>19</volume>, <fpage>2434</fpage>&#x2013;<lpage>2452</lpage>. doi: <pub-id pub-id-type="doi">10.1111/1462-2920.13768</pub-id>, PMID: <pub-id pub-id-type="pmid">28418097</pub-id></citation>
</ref>
<ref id="ref3">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ahmad</surname> <given-names>I.</given-names></name> <name><surname>Krishnamurthy</surname> <given-names>V.</given-names></name> <name><surname>Rao</surname> <given-names>D. N.</given-names></name></person-group> (<year>1995</year>). <article-title>DNA recognition by the EcoP15I and EcoPI modification methyltransferases</article-title>. <source>Gene</source> <volume>157</volume>, <fpage>143</fpage>&#x2013;<lpage>147</lpage>. doi: <pub-id pub-id-type="doi">10.1016/0378-1119(95)00671-R</pub-id>, PMID: <pub-id pub-id-type="pmid">7607479</pub-id></citation>
</ref>
<ref id="ref4">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Anton</surname> <given-names>B. P.</given-names></name> <name><surname>Roberts</surname> <given-names>R. J.</given-names></name></person-group> (<year>2021</year>). <article-title>Beyond restriction modification: Epigenomic roles of DNA methylation in prokaryotes</article-title>. <source>Annu. Rev. Microbiol.</source> <volume>75</volume>, <fpage>129</fpage>&#x2013;<lpage>149</lpage>. doi: <pub-id pub-id-type="doi">10.1146/annurev-micro-040521-035040</pub-id>, PMID: <pub-id pub-id-type="pmid">34314594</pub-id></citation>
</ref>
<ref id="ref5">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Arora</surname> <given-names>J.</given-names></name> <name><surname>Goswami</surname> <given-names>K.</given-names></name> <name><surname>Saha</surname> <given-names>S.</given-names></name></person-group> (<year>2014</year>). <article-title>Characterization of the replication initiator Orc1/Cdc6 from the archaeon <italic>Picrophilus torridus</italic></article-title>. <source>J. Bacteriol.</source> <volume>196</volume>, <fpage>276</fpage>&#x2013;<lpage>286</lpage>. doi: <pub-id pub-id-type="doi">10.1128/JB.01020-13</pub-id>, PMID: <pub-id pub-id-type="pmid">24187082</pub-id></citation>
</ref>
<ref id="ref6">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Barbeyron</surname> <given-names>T.</given-names></name> <name><surname>Kean</surname> <given-names>K.</given-names></name> <name><surname>Forterre</surname> <given-names>P.</given-names></name></person-group> (<year>1984</year>). <article-title>DNA adenine methylation of GATC sequences appeared recently in the <italic>Escherichia coli</italic> lineage</article-title>. <source>J. Bacteriol.</source> <volume>160</volume>, <fpage>586</fpage>&#x2013;<lpage>590</lpage>. doi: <pub-id pub-id-type="doi">10.1128/jb.160.2.586-590.1984</pub-id>, PMID: <pub-id pub-id-type="pmid">6094478</pub-id></citation>
</ref>
<ref id="ref7">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bath</surname> <given-names>A. J.</given-names></name> <name><surname>Milsom</surname> <given-names>S. E.</given-names></name> <name><surname>Gormley</surname> <given-names>N. A.</given-names></name> <name><surname>Halford</surname> <given-names>S. E.</given-names></name></person-group> (<year>2002</year>). <article-title>Many type IIs restriction endonucleases interact with two recognition sites before cleaving DNA</article-title>. <source>J. Biol. Chem.</source> <volume>277</volume>, <fpage>4024</fpage>&#x2013;<lpage>4033</lpage>. doi: <pub-id pub-id-type="doi">10.1074/jbc.M108441200</pub-id>, PMID: <pub-id pub-id-type="pmid">11729187</pub-id></citation>
</ref>
<ref id="ref8">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Benkert</surname> <given-names>P.</given-names></name> <name><surname>Kunzli</surname> <given-names>M.</given-names></name> <name><surname>Schwede</surname> <given-names>T.</given-names></name></person-group> (<year>2009</year>). <article-title>QMEAN server for protein model quality estimation</article-title>. <source>Nucleic Acids Res.</source> <volume>37</volume>, <fpage>W510</fpage>&#x2013;<lpage>W514</lpage>. doi: <pub-id pub-id-type="doi">10.1093/nar/gkp322</pub-id>, PMID: <pub-id pub-id-type="pmid">19429685</pub-id></citation>
</ref>
<ref id="ref9">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bheemanaik</surname> <given-names>S.</given-names></name> <name><surname>Reddy</surname> <given-names>Y. V.</given-names></name> <name><surname>Rao</surname> <given-names>D. N.</given-names></name></person-group> (<year>2006</year>). <article-title>Structure, function and mechanism of exocyclic DNA methyltransferases</article-title>. <source>Biochem. J.</source> <volume>399</volume>, <fpage>177</fpage>&#x2013;<lpage>190</lpage>. doi: <pub-id pub-id-type="doi">10.1042/BJ20060854</pub-id>, PMID: <pub-id pub-id-type="pmid">16987108</pub-id></citation>
</ref>
<ref id="ref10">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bist</surname> <given-names>P.</given-names></name> <name><surname>Rao</surname> <given-names>D. N.</given-names></name></person-group> (<year>2003</year>). <article-title>Identification and mutational analysis of Mg2+ binding site in EcoP15I DNA methyltransferase: involvement in target base eversion</article-title>. <source>J. Biol. Chem.</source> <volume>278</volume>, <fpage>41837</fpage>&#x2013;<lpage>41848</lpage>. doi: <pub-id pub-id-type="doi">10.1074/jbc.M307053200</pub-id>, PMID: <pub-id pub-id-type="pmid">12917398</pub-id></citation>
</ref>
<ref id="ref11">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Blow</surname> <given-names>M. J.</given-names></name> <name><surname>Clark</surname> <given-names>T. A.</given-names></name> <name><surname>Daum</surname> <given-names>C. G.</given-names></name> <name><surname>Deutschbauer</surname> <given-names>A. M.</given-names></name> <name><surname>Fomenkov</surname> <given-names>A.</given-names></name> <name><surname>Fries</surname> <given-names>R.</given-names></name> <etal/></person-group>. (<year>2016</year>). <article-title>The Epigenomic landscape of prokaryotes</article-title>. <source>PLoS Genet.</source> <volume>12</volume>:<fpage>e1005854</fpage>. doi: <pub-id pub-id-type="doi">10.1371/journal.pgen.1005854</pub-id>, PMID: <pub-id pub-id-type="pmid">26870957</pub-id></citation>
</ref>
<ref id="ref12">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Brunet</surname> <given-names>Y. R.</given-names></name> <name><surname>Bernard</surname> <given-names>C. S.</given-names></name> <name><surname>Cascales</surname> <given-names>E.</given-names></name></person-group> (<year>2020</year>). <article-title>Fur-dam regulatory interplay at an internal promoter of the Enteroaggregative <italic>Escherichia coli</italic> type VI secretion sci1 gene cluster</article-title>. <source>J. Bacteriol.</source> <volume>202</volume>:<fpage>20</fpage>. doi: <pub-id pub-id-type="doi">10.1128/JB.00075-20</pub-id>, PMID: <pub-id pub-id-type="pmid">32152218</pub-id></citation>
</ref>
<ref id="ref13">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bujnicki</surname> <given-names>J. M.</given-names></name></person-group> (<year>2002</year>). <article-title>Sequence permutations in the molecular evolution of DNA methyltransferases</article-title>. <source>BMC Evol. Biol.</source> <volume>2</volume>:<fpage>3</fpage>. doi: <pub-id pub-id-type="doi">10.1186/1471-2148-2-3</pub-id>, PMID: <pub-id pub-id-type="pmid">11914127</pub-id></citation>
</ref>
<ref id="ref14">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Couturier</surname> <given-names>M.</given-names></name> <name><surname>Lindas</surname> <given-names>A. C.</given-names></name></person-group> (<year>2018</year>). <article-title>The DNA methylome of the hyperthermoacidophilic crenarchaeon <italic>Sulfolobus acidocaldarius</italic></article-title>. <source>Front. Microbiol.</source> <volume>9</volume>:<fpage>137</fpage>. doi: <pub-id pub-id-type="doi">10.3389/fmicb.2018.00137</pub-id>, PMID: <pub-id pub-id-type="pmid">29472906</pub-id></citation>
</ref>
<ref id="ref15">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Dreier</surname> <given-names>J.</given-names></name> <name><surname>Macwilliams</surname> <given-names>M. P.</given-names></name> <name><surname>Bickle</surname> <given-names>T. A.</given-names></name></person-group> (<year>1996</year>). <article-title>DNA cleavage by the type IC restriction-modification enzyme EcoR124II</article-title>. <source>J. Mol. Biol.</source> <volume>264</volume>, <fpage>722</fpage>&#x2013;<lpage>733</lpage>. doi: <pub-id pub-id-type="doi">10.1006/jmbi.1996.0672</pub-id>, PMID: <pub-id pub-id-type="pmid">8980681</pub-id></citation>
</ref>
<ref id="ref16">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Falker</surname> <given-names>S.</given-names></name> <name><surname>Schilling</surname> <given-names>J.</given-names></name> <name><surname>Schmidt</surname> <given-names>M. A.</given-names></name> <name><surname>Heusipp</surname> <given-names>G.</given-names></name></person-group> (<year>2007</year>). <article-title>Overproduction of DNA adenine methyltransferase alters motility, invasion, and the lipopolysaccharide O-antigen composition of <italic>Yersinia enterocolitica</italic></article-title>. <source>Infect. Immun.</source> <volume>75</volume>, <fpage>4990</fpage>&#x2013;<lpage>4997</lpage>. doi: <pub-id pub-id-type="doi">10.1128/IAI.00457-07</pub-id>, PMID: <pub-id pub-id-type="pmid">17682042</pub-id></citation>
</ref>
<ref id="ref17">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Futterer</surname> <given-names>O.</given-names></name> <name><surname>Angelov</surname> <given-names>A.</given-names></name> <name><surname>Liesegang</surname> <given-names>H.</given-names></name> <name><surname>Gottschalk</surname> <given-names>G.</given-names></name> <name><surname>Schleper</surname> <given-names>C.</given-names></name> <name><surname>Schepers</surname> <given-names>B.</given-names></name> <etal/></person-group>. (<year>2004</year>). <article-title>Genome sequence of <italic>Picrophilus torridus</italic> and its implications for life around pH 0</article-title>. <source>Proc. Natl. Acad. Sci. U. S. A.</source> <volume>101</volume>, <fpage>9091</fpage>&#x2013;<lpage>9096</lpage>. doi: <pub-id pub-id-type="doi">10.1073/pnas.0401356101</pub-id>, PMID: <pub-id pub-id-type="pmid">15184674</pub-id></citation>
</ref>
<ref id="ref18">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gao</surname> <given-names>Y.</given-names></name> <name><surname>Cao</surname> <given-names>D.</given-names></name> <name><surname>Zhu</surname> <given-names>J.</given-names></name> <name><surname>Feng</surname> <given-names>H.</given-names></name> <name><surname>Luo</surname> <given-names>X.</given-names></name> <name><surname>Liu</surname> <given-names>S.</given-names></name> <etal/></person-group>. (<year>2020</year>). <article-title>Structural insights into assembly, operation and inhibition of a type I restriction-modification system</article-title>. <source>Nat. Microbiol.</source> <volume>5</volume>, <fpage>1107</fpage>&#x2013;<lpage>1118</lpage>. doi: <pub-id pub-id-type="doi">10.1038/s41564-020-0731-z</pub-id>, PMID: <pub-id pub-id-type="pmid">32483229</pub-id></citation>
</ref>
<ref id="ref19">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Goedecke</surname> <given-names>K.</given-names></name> <name><surname>Pignot</surname> <given-names>M.</given-names></name> <name><surname>Goody</surname> <given-names>R. S.</given-names></name> <name><surname>Scheidig</surname> <given-names>A. J.</given-names></name> <name><surname>Weinhold</surname> <given-names>E.</given-names></name></person-group> (<year>2001</year>). <article-title>Structure of the N6-adenine DNA methyltransferase M.TaqI in complex with DNA and a cofactor analog</article-title>. <source>Nat. Struct. Biol.</source> <volume>8</volume>, <fpage>121</fpage>&#x2013;<lpage>125</lpage>. doi: <pub-id pub-id-type="doi">10.1038/84104</pub-id></citation>
</ref>
<ref id="ref20">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Grogan</surname> <given-names>D. W.</given-names></name></person-group> (<year>2003</year>). <article-title>Cytosine methylation by the SuaI restriction-modification system: implications for genetic fidelity in a hyperthermophilic archaeon</article-title>. <source>J. Bacteriol.</source> <volume>185</volume>, <fpage>4657</fpage>&#x2013;<lpage>4661</lpage>. doi: <pub-id pub-id-type="doi">10.1128/JB.185.15.4657-4661.2003</pub-id>, PMID: <pub-id pub-id-type="pmid">12867480</pub-id></citation>
</ref>
<ref id="ref21">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ishikawa</surname> <given-names>K.</given-names></name> <name><surname>Watanabe</surname> <given-names>M.</given-names></name> <name><surname>Kuroita</surname> <given-names>T.</given-names></name> <name><surname>Uchiyama</surname> <given-names>I.</given-names></name> <name><surname>Bujnicki</surname> <given-names>J. M.</given-names></name> <name><surname>Kawakami</surname> <given-names>B.</given-names></name> <etal/></person-group>. (<year>2005</year>). <article-title>Discovery of a novel restriction endonuclease by genome comparison and application of a wheat-germ-based cell-free translation assay: PabI (5&#x2032;-GTA/C) from the hyperthermophilic archaeon <italic>Pyrococcus abyssi</italic></article-title>. <source>Nucleic Acids Res.</source> <volume>33</volume>:<fpage>e112</fpage>. doi: <pub-id pub-id-type="doi">10.1093/nar/gni113</pub-id>, PMID: <pub-id pub-id-type="pmid">16040595</pub-id></citation>
</ref>
<ref id="ref22">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Jurenaite-Urbanaviciene</surname> <given-names>S.</given-names></name> <name><surname>Kazlauskiene</surname> <given-names>R.</given-names></name> <name><surname>Urbelyte</surname> <given-names>V.</given-names></name> <name><surname>Maneliene</surname> <given-names>Z.</given-names></name> <name><surname>Petrusyte</surname> <given-names>M.</given-names></name> <name><surname>Lubys</surname> <given-names>A.</given-names></name> <etal/></person-group>. (<year>2001</year>). <article-title>Characterization of BseMII, a new type IV restriction-modification system, which recognizes the pentanucleotide sequence 5&#x2032;-CTCAG(N)(10/8)</article-title>. <source>Nucleic Acids Res.</source> <volume>29</volume>, <fpage>895</fpage>&#x2013;<lpage>903</lpage>. doi: <pub-id pub-id-type="doi">10.1093/nar/29.4.895</pub-id>, PMID: <pub-id pub-id-type="pmid">11160921</pub-id></citation>
</ref>
<ref id="ref23">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Keceli Oguz</surname> <given-names>S.</given-names></name> <name><surname>Has</surname> <given-names>E. G.</given-names></name> <name><surname>Akcelik</surname> <given-names>N.</given-names></name> <name><surname>Akcelik</surname> <given-names>M.</given-names></name></person-group> (<year>2022</year>). <article-title>Phenotypic impacts and genetic regulation characteristics of the DNA adenine methylase gene (dam) in <italic>Salmonella typhimurium</italic> biofilm forms</article-title>. <source>Res. Microbiol.</source> <volume>174</volume>:<fpage>103991</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.resmic.2022.103991</pub-id></citation>
</ref>
<ref id="ref24">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kelley</surname> <given-names>L. A.</given-names></name> <name><surname>Mezulis</surname> <given-names>S.</given-names></name> <name><surname>Yates</surname> <given-names>C. M.</given-names></name> <name><surname>Wass</surname> <given-names>M. N.</given-names></name> <name><surname>Sternberg</surname> <given-names>M. J.</given-names></name></person-group> (<year>2015</year>). <article-title>The Phyre2 web portal for protein modeling, prediction and analysis</article-title>. <source>Nat. Protoc.</source> <volume>10</volume>, <fpage>845</fpage>&#x2013;<lpage>858</lpage>. doi: <pub-id pub-id-type="doi">10.1038/nprot.2015.053</pub-id>, PMID: <pub-id pub-id-type="pmid">25950237</pub-id></citation>
</ref>
<ref id="ref25">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kennaway</surname> <given-names>C. K.</given-names></name> <name><surname>Obarska-Kosinska</surname> <given-names>A.</given-names></name> <name><surname>White</surname> <given-names>J. H.</given-names></name> <name><surname>Tuszynska</surname> <given-names>I.</given-names></name> <name><surname>Cooper</surname> <given-names>L. P.</given-names></name> <name><surname>Bujnicki</surname> <given-names>J. M.</given-names></name> <etal/></person-group>. (<year>2009</year>). <article-title>The structure of M.EcoKI type I DNA methyltransferase with a DNA mimic antirestriction protein</article-title>. <source>Nucleic Acids Res.</source> <volume>37</volume>, <fpage>762</fpage>&#x2013;<lpage>770</lpage>. doi: <pub-id pub-id-type="doi">10.1093/nar/gkn988</pub-id>, PMID: <pub-id pub-id-type="pmid">19074193</pub-id></citation>
</ref>
<ref id="ref26">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kim</surname> <given-names>J. S.</given-names></name> <name><surname>Degiovanni</surname> <given-names>A.</given-names></name> <name><surname>Jancarik</surname> <given-names>J.</given-names></name> <name><surname>Adams</surname> <given-names>P. D.</given-names></name> <name><surname>Yokota</surname> <given-names>H.</given-names></name> <name><surname>Kim</surname> <given-names>R.</given-names></name> <etal/></person-group>. (<year>2005</year>). <article-title>Crystal structure of DNA sequence specificity subunit of a type I restriction-modification enzyme and its functional implications</article-title>. <source>Proc. Natl. Acad. Sci. U. S. A.</source> <volume>102</volume>, <fpage>3248</fpage>&#x2013;<lpage>3253</lpage>. doi: <pub-id pub-id-type="doi">10.1073/pnas.0409851102</pub-id>, PMID: <pub-id pub-id-type="pmid">15728358</pub-id></citation>
</ref>
<ref id="ref27">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Koike</surname> <given-names>H.</given-names></name> <name><surname>Yokoyama</surname> <given-names>K.</given-names></name> <name><surname>Kawashima</surname> <given-names>T.</given-names></name> <name><surname>Yamasaki</surname> <given-names>T.</given-names></name> <name><surname>Makino</surname> <given-names>S.</given-names></name> <name><surname>Clowney</surname> <given-names>L.</given-names></name> <etal/></person-group>. (<year>2005</year>). <article-title>GATC methylation by dam methylase in bacteria: its roles and possible transcription regulation by an FFRP</article-title>. <source>Proc. Jpn. Acad.</source> <volume>81</volume>, <fpage>278</fpage>&#x2013;<lpage>290</lpage>. doi: <pub-id pub-id-type="doi">10.2183/pjab.81.278</pub-id></citation>
</ref>
<ref id="ref28">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Labahn</surname> <given-names>J.</given-names></name> <name><surname>Granzin</surname> <given-names>J.</given-names></name> <name><surname>Schluckebier</surname> <given-names>G.</given-names></name> <name><surname>Robinson</surname> <given-names>D. P.</given-names></name> <name><surname>Jack</surname> <given-names>W. E.</given-names></name> <name><surname>Schildkraut</surname> <given-names>I.</given-names></name> <etal/></person-group>. (<year>1994</year>). <article-title>Three-dimensional structure of the adenine-specific DNA methyltransferase M.Taq I in complex with the cofactor S-adenosylmethionine</article-title>. <source>Proc. Natl. Acad. Sci. U. S. A.</source> <volume>91</volume>, <fpage>10957</fpage>&#x2013;<lpage>10961</lpage>. doi: <pub-id pub-id-type="doi">10.1073/pnas.91.23.10957</pub-id>, PMID: <pub-id pub-id-type="pmid">7971991</pub-id></citation>
</ref>
<ref id="ref29">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lautenberger</surname> <given-names>J. A.</given-names></name> <name><surname>Linn</surname> <given-names>S.</given-names></name></person-group> (<year>1972</year>). <article-title>The deoxyribonucleic acid modification and restriction enzymes of <italic>Escherichia coli</italic> B. I. Purification, subunit structure, and catalytic properties of the modification methylase</article-title>. <source>J. Biol. Chem.</source> <volume>247</volume>, <fpage>6176</fpage>&#x2013;<lpage>6182</lpage>. doi: <pub-id pub-id-type="doi">10.1016/S0021-9258(19)44779-0</pub-id></citation>
</ref>
<ref id="ref30">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lee</surname> <given-names>S. H.</given-names></name> <name><surname>Kim</surname> <given-names>M. S.</given-names></name> <name><surname>Lee</surname> <given-names>J. H.</given-names></name> <name><surname>Kim</surname> <given-names>T. W.</given-names></name> <name><surname>Bae</surname> <given-names>S. S.</given-names></name> <name><surname>Lee</surname> <given-names>S. M.</given-names></name> <etal/></person-group>. (<year>2016</year>). <article-title>Adaptive engineering of a hyperthermophilic archaeon on CO and discovering the underlying mechanism by multi-omics analysis</article-title>. <source>Sci. Rep.</source> <volume>6</volume>:<fpage>22896</fpage>. doi: <pub-id pub-id-type="doi">10.1038/srep22896</pub-id>, PMID: <pub-id pub-id-type="pmid">26975345</pub-id></citation>
</ref>
<ref id="ref31">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Loenen</surname> <given-names>W. A.</given-names></name> <name><surname>Dryden</surname> <given-names>D. T.</given-names></name> <name><surname>Raleigh</surname> <given-names>E. A.</given-names></name> <name><surname>Wilson</surname> <given-names>G. G.</given-names></name></person-group> (<year>2014a</year>). <article-title>Type I restriction enzymes and their relatives</article-title>. <source>Nucleic Acids Res.</source> <volume>42</volume>, <fpage>20</fpage>&#x2013;<lpage>44</lpage>. doi: <pub-id pub-id-type="doi">10.1093/nar/gkt847</pub-id>, PMID: <pub-id pub-id-type="pmid">24068554</pub-id></citation>
</ref>
<ref id="ref32">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Loenen</surname> <given-names>W. A.</given-names></name> <name><surname>Dryden</surname> <given-names>D. T.</given-names></name> <name><surname>Raleigh</surname> <given-names>E. A.</given-names></name> <name><surname>Wilson</surname> <given-names>G. G.</given-names></name> <name><surname>Murray</surname> <given-names>N. E.</given-names></name></person-group> (<year>2014b</year>). <article-title>Highlights of the DNA cutters: a short history of the restriction enzymes</article-title>. <source>Nucleic Acids Res.</source> <volume>42</volume>, <fpage>3</fpage>&#x2013;<lpage>19</lpage>. doi: <pub-id pub-id-type="doi">10.1093/nar/gkt990</pub-id>, PMID: <pub-id pub-id-type="pmid">24141096</pub-id></citation>
</ref>
<ref id="ref33">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lunnen</surname> <given-names>K. D.</given-names></name> <name><surname>Morgan</surname> <given-names>R. D.</given-names></name> <name><surname>Timan</surname> <given-names>C. J.</given-names></name> <name><surname>Krzycki</surname> <given-names>J. A.</given-names></name> <name><surname>Reeve</surname> <given-names>J. N.</given-names></name> <name><surname>Wilson</surname> <given-names>G. G.</given-names></name></person-group> (<year>1989</year>). <article-title>Characterization and cloning of MwoI (GCN7GC), a new type-II restriction-modification system from <italic>Methanobacterium wolfei</italic></article-title>. <source>Gene</source> <volume>77</volume>, <fpage>11</fpage>&#x2013;<lpage>19</lpage>. doi: <pub-id pub-id-type="doi">10.1016/0378-1119(89)90354-5</pub-id>, PMID: <pub-id pub-id-type="pmid">2663652</pub-id></citation>
</ref>
<ref id="ref34">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Malone</surname> <given-names>T.</given-names></name> <name><surname>Blumenthal</surname> <given-names>R. M.</given-names></name> <name><surname>Cheng</surname> <given-names>X.</given-names></name></person-group> (<year>1995</year>). <article-title>Structure-guided analysis reveals nine sequence motifs conserved among DNA amino-methyltransferases, and suggests a catalytic mechanism for these enzymes</article-title>. <source>J. Mol. Biol.</source> <volume>253</volume>, <fpage>618</fpage>&#x2013;<lpage>632</lpage>. doi: <pub-id pub-id-type="doi">10.1006/jmbi.1995.0577</pub-id>, PMID: <pub-id pub-id-type="pmid">7473738</pub-id></citation>
</ref>
<ref id="ref35">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Maynard-Smith</surname> <given-names>M. D.</given-names></name> <name><surname>Mckelvie</surname> <given-names>J. C.</given-names></name> <name><surname>Wood</surname> <given-names>R. J.</given-names></name> <name><surname>Harmer</surname> <given-names>J. E.</given-names></name> <name><surname>Ranasinghe</surname> <given-names>R. T.</given-names></name> <name><surname>Williams</surname> <given-names>C. L.</given-names></name> <etal/></person-group>. (<year>2011</year>). <article-title>Direct and continuous fluorescence-based measurements of <italic>Pyrococcus horikoshii</italic> DNA N-6 adenine methyltransferase activity</article-title>. <source>Anal. Biochem.</source> <volume>418</volume>, <fpage>204</fpage>&#x2013;<lpage>212</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.ab.2011.07.023</pub-id>, PMID: <pub-id pub-id-type="pmid">21839719</pub-id></citation>
</ref>
<ref id="ref36">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mcconnell</surname> <given-names>D. J.</given-names></name> <name><surname>Searcy</surname> <given-names>D. G.</given-names></name> <name><surname>Sutcliffe</surname> <given-names>J. G.</given-names></name></person-group> (<year>1978</year>). <article-title>A restriction enzyme Tha I from the thermophilic mycoplasma <italic>Thermoplasma acidophilum</italic></article-title>. <source>Nucleic Acids Res.</source> <volume>5</volume>, <fpage>1729</fpage>&#x2013;<lpage>1739</lpage>. doi: <pub-id pub-id-type="doi">10.1093/nar/5.6.1729</pub-id>, PMID: <pub-id pub-id-type="pmid">673834</pub-id></citation>
</ref>
<ref id="ref37">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mooers</surname> <given-names>B. H. M.</given-names></name></person-group> (<year>2016</year>). <article-title>Simplifying and enhancing the use of PyMOL with horizontal scripts</article-title>. <source>Protein Sci.</source> <volume>25</volume>, <fpage>1873</fpage>&#x2013;<lpage>1882</lpage>. doi: <pub-id pub-id-type="doi">10.1002/pro.2996</pub-id>, PMID: <pub-id pub-id-type="pmid">27488983</pub-id></citation>
</ref>
<ref id="ref38">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Morgan</surname> <given-names>R.</given-names></name> <name><surname>Xiao</surname> <given-names>J.</given-names></name> <name><surname>Xu</surname> <given-names>S.</given-names></name></person-group> (<year>1998</year>). <article-title>Characterization of an extremely thermostable restriction enzyme, PspGI, from a Pyrococcus strain and cloning of the PspGI restriction-modification system in <italic>Escherichia coli</italic></article-title>. <source>Appl. Environ. Microbiol.</source> <volume>64</volume>, <fpage>3669</fpage>&#x2013;<lpage>3673</lpage>. doi: <pub-id pub-id-type="doi">10.1128/AEM.64.10.3669-3673.1998</pub-id>, PMID: <pub-id pub-id-type="pmid">9758783</pub-id></citation>
</ref>
<ref id="ref39">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Nolling</surname> <given-names>J.</given-names></name> <name><surname>De Vos</surname> <given-names>W. M.</given-names></name></person-group> (<year>1992a</year>). <article-title>Characterization of the archaeal, plasmid-encoded type II restriction-modification system MthTI from <italic>Methanobacterium thermoformicicum</italic> THF: homology to the bacterial NgoPII system from Neisseria gonorrhoeae</article-title>. <source>J. Bacteriol.</source> <volume>174</volume>, <fpage>5719</fpage>&#x2013;<lpage>5726</lpage>. doi: <pub-id pub-id-type="doi">10.1128/jb.174.17.5719-5726.1992</pub-id>, PMID: <pub-id pub-id-type="pmid">1512204</pub-id></citation>
</ref>
<ref id="ref40">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Nolling</surname> <given-names>J.</given-names></name> <name><surname>De Vos</surname> <given-names>W. M.</given-names></name></person-group> (<year>1992b</year>). <article-title>Identification of the CTAG-recognizing restriction-modification systems MthZI and MthFI from <italic>Methanobacterium thermoformicicum</italic> and characterization of the plasmid-encoded mthZIM gene</article-title>. <source>Nucleic Acids Res.</source> <volume>20</volume>, <fpage>5047</fpage>&#x2013;<lpage>5052</lpage>. doi: <pub-id pub-id-type="doi">10.1093/nar/20.19.5047</pub-id>, PMID: <pub-id pub-id-type="pmid">1408820</pub-id></citation>
</ref>
<ref id="ref41">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ouellette</surname> <given-names>M.</given-names></name> <name><surname>Gogarten</surname> <given-names>J. P.</given-names></name> <name><surname>Lajoie</surname> <given-names>J.</given-names></name> <name><surname>Makkay</surname> <given-names>A. M.</given-names></name> <name><surname>Papke</surname> <given-names>R. T.</given-names></name></person-group> (<year>2018</year>). <article-title>Characterizing the DNA methyltransferases of <italic>Haloferax volcanii</italic> via bioinformatics, gene deletion, and SMRT sequencing</article-title>. <source>Genes (Basel)</source> <volume>9</volume>:<fpage>129</fpage>. doi: <pub-id pub-id-type="doi">10.3390/genes9030129</pub-id>, PMID: <pub-id pub-id-type="pmid">29495512</pub-id></citation>
</ref>
<ref id="ref42">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ouellette</surname> <given-names>M.</given-names></name> <name><surname>Jackson</surname> <given-names>L.</given-names></name> <name><surname>Chimileski</surname> <given-names>S.</given-names></name> <name><surname>Papke</surname> <given-names>R. T.</given-names></name></person-group> (<year>2015</year>). <article-title>Genome-wide DNA methylation analysis of <italic>Haloferax volcanii</italic> H26 and identification of DNA methyltransferase related PD-(D/E)XK nuclease family protein HVO_A0006</article-title>. <source>Front. Microbiol.</source> <volume>6</volume>:<fpage>251</fpage>. doi: <pub-id pub-id-type="doi">10.3389/fmicb.2015.00251</pub-id></citation>
</ref>
<ref id="ref43">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Pingoud</surname> <given-names>V.</given-names></name> <name><surname>Conzelmann</surname> <given-names>C.</given-names></name> <name><surname>Kinzebach</surname> <given-names>S.</given-names></name> <name><surname>Sudina</surname> <given-names>A.</given-names></name> <name><surname>Metelev</surname> <given-names>V.</given-names></name> <name><surname>Kubareva</surname> <given-names>E.</given-names></name> <etal/></person-group>. (<year>2003</year>). <article-title>PspGI, a type II restriction endonuclease from the extreme thermophile Pyrococcus sp.: structural and functional studies to investigate an evolutionary relationship with several mesophilic restriction enzymes</article-title>. <source>J. Mol. Biol.</source> <volume>329</volume>, <fpage>913</fpage>&#x2013;<lpage>929</lpage>. doi: <pub-id pub-id-type="doi">10.1016/S0022-2836(03)00523-0</pub-id>, PMID: <pub-id pub-id-type="pmid">12798682</pub-id></citation>
</ref>
<ref id="ref44">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Pingoud</surname> <given-names>A.</given-names></name> <name><surname>Wilson</surname> <given-names>G. G.</given-names></name> <name><surname>Wende</surname> <given-names>W.</given-names></name></person-group> (<year>2014</year>). <article-title>Type II restriction endonucleases&#x2014;a historical perspective and more</article-title>. <source>Nucleic Acids Res.</source> <volume>42</volume>, <fpage>7489</fpage>&#x2013;<lpage>7527</lpage>. doi: <pub-id pub-id-type="doi">10.1093/nar/gku447</pub-id>, PMID: <pub-id pub-id-type="pmid">24878924</pub-id></citation>
</ref>
<ref id="ref45">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Pogolotti</surname> <given-names>A. L.</given-names> <suffix>Jr.</suffix></name> <name><surname>Ono</surname> <given-names>A.</given-names></name> <name><surname>Subramaniam</surname> <given-names>R.</given-names></name> <name><surname>Santi</surname> <given-names>D. V.</given-names></name></person-group> (<year>1988</year>). <article-title>On the mechanism of DNA-adenine methylase</article-title>. <source>J. Biol. Chem.</source> <volume>263</volume>, <fpage>7461</fpage>&#x2013;<lpage>7464</lpage>. doi: <pub-id pub-id-type="doi">10.1016/S0021-9258(18)68520-5</pub-id></citation>
</ref>
<ref id="ref46">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Prangishvili</surname> <given-names>D. A.</given-names></name> <name><surname>Vashakidze</surname> <given-names>R. P.</given-names></name> <name><surname>Chelidze</surname> <given-names>M. G.</given-names></name> <name><surname>Gabriadze</surname> <given-names>I.</given-names></name></person-group> (<year>1985</year>). <article-title>A restriction endonuclease SuaI from the thermoacidophilic archaebacterium <italic>Sulfolobus acidocaldarius</italic></article-title>. <source>FEBS Lett.</source> <volume>192</volume>, <fpage>57</fpage>&#x2013;<lpage>60</lpage>. doi: <pub-id pub-id-type="doi">10.1016/0014-5793(85)80042-9</pub-id>, PMID: <pub-id pub-id-type="pmid">2996942</pub-id></citation>
</ref>
<ref id="ref47">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Pues</surname> <given-names>H.</given-names></name> <name><surname>Bleimling</surname> <given-names>N.</given-names></name> <name><surname>Holz</surname> <given-names>B.</given-names></name> <name><surname>Wolcke</surname> <given-names>J.</given-names></name> <name><surname>Weinhold</surname> <given-names>E.</given-names></name></person-group> (<year>1999</year>). <article-title>Functional roles of the conserved aromatic amino acid residues at position 108 (motif IV) and position 196 (motif VIII) in base flipping and catalysis by the N6-adenine DNA methyltransferase from Thermus aquaticus</article-title>. <source>Biochemistry</source> <volume>38</volume>, <fpage>1426</fpage>&#x2013;<lpage>1434</lpage>. doi: <pub-id pub-id-type="doi">10.1021/bi9818016</pub-id>, PMID: <pub-id pub-id-type="pmid">9931007</pub-id></citation>
</ref>
<ref id="ref48">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Rao</surname> <given-names>D. N.</given-names></name> <name><surname>Dryden</surname> <given-names>D. T.</given-names></name> <name><surname>Bheemanaik</surname> <given-names>S.</given-names></name></person-group> (<year>2014</year>). <article-title>Type III restriction-modification enzymes: a historical perspective</article-title>. <source>Nucleic Acids Res.</source> <volume>42</volume>, <fpage>45</fpage>&#x2013;<lpage>55</lpage>. doi: <pub-id pub-id-type="doi">10.1093/nar/gkt616</pub-id>, PMID: <pub-id pub-id-type="pmid">23863841</pub-id></citation>
</ref>
<ref id="ref49">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Roth</surname> <given-names>M.</given-names></name> <name><surname>Helm-Kruse</surname> <given-names>S.</given-names></name> <name><surname>Friedrich</surname> <given-names>T.</given-names></name> <name><surname>Jeltsch</surname> <given-names>A.</given-names></name></person-group> (<year>1998</year>). <article-title>Functional roles of conserved amino acid residues in DNA methyltransferases investigated by site-directed mutagenesis of the EcoRV adenine-N6-methyltransferase</article-title>. <source>J. Biol. Chem.</source> <volume>273</volume>, <fpage>17333</fpage>&#x2013;<lpage>17342</lpage>. doi: <pub-id pub-id-type="doi">10.1074/jbc.273.28.17333</pub-id>, PMID: <pub-id pub-id-type="pmid">9651316</pub-id></citation>
</ref>
<ref id="ref50">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Schleper</surname> <given-names>C.</given-names></name> <name><surname>Puehler</surname> <given-names>G.</given-names></name> <name><surname>Holz</surname> <given-names>I.</given-names></name> <name><surname>Gambacorta</surname> <given-names>A.</given-names></name> <name><surname>Janekovic</surname> <given-names>D.</given-names></name> <name><surname>Santarius</surname> <given-names>U.</given-names></name> <etal/></person-group>. (<year>1995</year>). <article-title>Picrophilus gen. Nov., fam. Nov.: a novel aerobic, heterotrophic, thermoacidophilic genus and family comprising archaea capable of growth around pH 0</article-title>. <source>J. Bacteriol.</source> <volume>177</volume>, <fpage>7050</fpage>&#x2013;<lpage>7059</lpage>. doi: <pub-id pub-id-type="doi">10.1128/jb.177.24.7050-7059.1995</pub-id>, PMID: <pub-id pub-id-type="pmid">8522509</pub-id></citation>
</ref>
<ref id="ref51">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Schmid</surname> <given-names>K.</given-names></name> <name><surname>Thomm</surname> <given-names>M.</given-names></name> <name><surname>Laminet</surname> <given-names>A.</given-names></name> <name><surname>Laue</surname> <given-names>F. G.</given-names></name> <name><surname>Kessler</surname> <given-names>C.</given-names></name> <name><surname>Stetter</surname> <given-names>K. O.</given-names></name> <etal/></person-group>. (<year>1984</year>). <article-title>Three new restriction endonucleases MaeI, MaeII and MaeIII from <italic>Methanococcus aeolicus</italic></article-title>. <source>Nucleic Acids Res.</source> <volume>12</volume>, <fpage>2619</fpage>&#x2013;<lpage>2628</lpage>. doi: <pub-id pub-id-type="doi">10.1093/nar/12.6.2619</pub-id>, PMID: <pub-id pub-id-type="pmid">6324124</pub-id></citation>
</ref>
<ref id="ref52">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Sievers</surname> <given-names>F.</given-names></name> <name><surname>Wilm</surname> <given-names>A.</given-names></name> <name><surname>Dineen</surname> <given-names>D.</given-names></name> <name><surname>Gibson</surname> <given-names>T. J.</given-names></name> <name><surname>Karplus</surname> <given-names>K.</given-names></name> <name><surname>Li</surname> <given-names>W.</given-names></name> <etal/></person-group>. (<year>2011</year>). <article-title>Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal omega</article-title>. <source>Mol. Syst. Biol.</source> <volume>7</volume>:<fpage>539</fpage>. doi: <pub-id pub-id-type="doi">10.1038/msb.2011.75</pub-id>, PMID: <pub-id pub-id-type="pmid">21988835</pub-id></citation>
</ref>
<ref id="ref53">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Suzuki</surname> <given-names>S.</given-names></name> <name><surname>Kurosawa</surname> <given-names>N.</given-names></name></person-group> (<year>2016</year>). <article-title>Disruption of the gene encoding restriction endonuclease SuaI and development of a host-vector system for the thermoacidophilic archaeon <italic>Sulfolobus acidocaldarius</italic></article-title>. <source>Extremophiles</source> <volume>20</volume>, <fpage>139</fpage>&#x2013;<lpage>148</lpage>. doi: <pub-id pub-id-type="doi">10.1007/s00792-016-0807-0</pub-id>, PMID: <pub-id pub-id-type="pmid">26791382</pub-id></citation>
</ref>
<ref id="ref54">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Tucholski</surname> <given-names>J.</given-names></name> <name><surname>Skowron</surname> <given-names>P. M.</given-names></name> <name><surname>Podhajska</surname> <given-names>A. J.</given-names></name></person-group> (<year>1995</year>). <article-title>MmeI, a class-IIS restriction endonuclease: purification and characterization</article-title>. <source>Gene</source> <volume>157</volume>, <fpage>87</fpage>&#x2013;<lpage>92</lpage>. doi: <pub-id pub-id-type="doi">10.1016/0378-1119(94)00787-S</pub-id>, PMID: <pub-id pub-id-type="pmid">7607532</pub-id></citation>
</ref>
<ref id="ref55">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Watanabe</surname> <given-names>M.</given-names></name> <name><surname>Yuzawa</surname> <given-names>H.</given-names></name> <name><surname>Handa</surname> <given-names>N.</given-names></name> <name><surname>Kobayashi</surname> <given-names>I.</given-names></name></person-group> (<year>2006</year>). <article-title>Hyperthermophilic DNA methyltransferase M.PabI from the archaeon <italic>Pyrococcus abyssi</italic></article-title>. <source>Appl. Environ. Microbiol.</source> <volume>72</volume>, <fpage>5367</fpage>&#x2013;<lpage>5375</lpage>. doi: <pub-id pub-id-type="doi">10.1128/AEM.00433-06</pub-id>, PMID: <pub-id pub-id-type="pmid">16885288</pub-id></citation>
</ref>
<ref id="ref56">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Willcock</surname> <given-names>D. F.</given-names></name> <name><surname>Dryden</surname> <given-names>D. T.</given-names></name> <name><surname>Murray</surname> <given-names>N. E.</given-names></name></person-group> (<year>1994</year>). <article-title>A mutational analysis of the two motifs common to adenine methyltransferases</article-title>. <source>EMBO J.</source> <volume>13</volume>, <fpage>3902</fpage>&#x2013;<lpage>3908</lpage>. doi: <pub-id pub-id-type="doi">10.1002/j.1460-2075.1994.tb06701.x</pub-id>, PMID: <pub-id pub-id-type="pmid">8070417</pub-id></citation>
</ref>
<ref id="ref57">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zatopek</surname> <given-names>K. M.</given-names></name> <name><surname>Burkhart</surname> <given-names>B. W.</given-names></name> <name><surname>Morgan</surname> <given-names>R. D.</given-names></name> <name><surname>Gehring</surname> <given-names>A. M.</given-names></name> <name><surname>Scott</surname> <given-names>K. A.</given-names></name> <name><surname>Santangelo</surname> <given-names>T. J.</given-names></name> <etal/></person-group>. (<year>2021</year>). <article-title>The Hyperthermophilic restriction-modification Systems of <italic>Thermococcus kodakarensis</italic> protect genome integrity</article-title>. <source>Front. Microbiol.</source> <volume>12</volume>:<fpage>657356</fpage>. doi: <pub-id pub-id-type="doi">10.3389/fmicb.2021.657356</pub-id>, PMID: <pub-id pub-id-type="pmid">34093470</pub-id></citation>
</ref>
</ref-list>
<fn-group>
<fn id="fn0004">
<p>
<sup>1</sup>
<ext-link xlink:href="http://www.sbg.bio.ic.ac.uk/phyre2" ext-link-type="uri">http://www.sbg.bio.ic.ac.uk/phyre2</ext-link>
</p>
</fn>
<fn id="fn0005">
<p>
<sup>2</sup>
<ext-link xlink:href="https://proq.bioinfo.se/ProQ/ProQ.html" ext-link-type="uri">https://proq.bioinfo.se/ProQ/ProQ.html</ext-link>
</p>
</fn>
<fn id="fn0006">
<p>
<sup>3</sup>
<ext-link xlink:href="http://pymol.org" ext-link-type="uri">http://pymol.org</ext-link>
</p>
</fn>
<fn id="fn0007">
<p>
<sup>4</sup>
<ext-link xlink:href="http://www.genome.jp/kegg" ext-link-type="uri">www.genome.jp/kegg</ext-link>
</p>
</fn>
<fn id="fn0008">
<p>
<sup>5</sup>
<ext-link xlink:href="http://www.blast.ncbi.nlm.nih.gov" ext-link-type="uri">www.blast.ncbi.nlm.nih.gov</ext-link>
</p>
</fn>
</fn-group>
</back>
</article>