Figures
Abstract
Background
The genome of the halophilic archaeon Halobacterium salinarum NRC-1 encodes for homologs of MutS and MutL, which are key proteins of a DNA mismatch repair pathway conserved in Bacteria and Eukarya. Mismatch repair is essential for retaining the fidelity of genetic information and defects in this pathway result in the deleterious accumulation of mutations and in hereditary diseases in humans.
Methodology/Principal Findings
We calculated the spontaneous genomic mutation rate of H. salinarum NRC-1 using fluctuation tests targeting genes of the uracil monophosphate biosynthesis pathway. We found that H. salinarum NRC-1 has a low incidence of mutation suggesting the presence of active mechanisms to control spontaneous mutations during replication. The spectrum of mutational changes found in H. salinarum NRC-1, and in other archaea, appears to be unique to this domain of life and might be a consequence of their adaption to extreme environmental conditions. In-frame targeted gene deletions of H. salinarum NRC-1 mismatch repair genes and phenotypic characterization of the mutants demonstrated that the mutS and mutL genes are not required for maintenance of the observed mutation rate.
Conclusions/Significance
We established that H. salinarum NRC-1 mutS and mutL genes are redundant to an alternative system that limits spontaneous mutation in this organism. This finding leads to the puzzling question of what mechanism is responsible for maintenance of the low genomic mutation rates observed in the Archaea, which for the most part do not have MutS and MutL homologs.
Citation: Busch CR, DiRuggiero J (2010) MutS and MutL Are Dispensable for Maintenance of the Genomic Mutation Rate in the Halophilic Archaeon Halobacterium salinarum NRC-1. PLoS ONE 5(2): e9045. https://doi.org/10.1371/journal.pone.0009045
Editor: Rodolfo Aramayo, Texas A&M University, United States of America
Received: October 16, 2009; Accepted: January 5, 2010; Published: February 4, 2010
Copyright: © 2010 Busch, DiRuggiero. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This work was supported by grants from NASA (EXB04-0000-0055), NSF (MCB-0425825), and AFOSR ((FA95500710158) to J.D.R. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
DNA mismatch repair (MMR) is the major pathway for the repair of DNA replication errors such as nucleotide mismatches, insertions, and deletions [1]. Defects in the MMR pathway lead to genomic instability that can cause a 10 to 1000-fold increase in spontaneous mutability, meiotic defects in eukaryotes, and tolerance to DNA alkylating agents [1], [2], [3]. In humans, inactivation of the MMR pathway leads to a predisposition to hereditary nonpolyposis colon cancer and other types of tumors [2], [3]. The MMR pathway also plays an important role in preventing recombination events between divergent sequences [1], [2].
The key proteins of the MMR pathway, MutS and MutL, are highly conserved between Bacteria and Eukarya. The pathway has been characterized in bacterial and eukaryal systems and comprises three basic steps: (1) MutS/L recognition of mismatch, (2) excision of the mismatched base and surrounding DNA, and (3) repair synthesis [1], [2], [4], [5]. The three-dimensional structure of MutS has been resolved for the Escherichia coli and Thermus aquaticus proteins [6], [7]. MutS is a 95kDa protein that functions as a dimer in vivo [3], [8]. MutS has ATPase activity with Walker A/B sequence motifs and a highly conserved Phe-X-Glu motif responsible for binding DNA [9]. MutL is a 68kDa protein that exists as a dimer in solution and is a member of the Bergerat-fold ATPase/kinase family [3], [10]. Eukaryotes have multiple homologs of the MutS and MutL proteins that form heterodimers suggesting a more complex system than in bacteria with multiple interactions [1], [2]. In contrast to E. coli and several other gram-negative bacteria, eukaryotes and most bacteria do not have a methylation-directed MMR system for strand discrimination or a MutH homolog. Studies suggest a nick-directed mechanism using Okazaki fragments produced during replication of the lagging strand or a strand discrimination mechanism directed by the proliferating cell nuclear antigen (PCNA), thus coupling replication and MMR [2], [11].
MMR has not been investigated in the Archaea but studies of the genomic mutation rate in the thermophilic acidophile, Sulfolobus acidocaldarius, and the halophile, Haloferax volcanii, revealed rates of spontaneous mutation very similar to rates previously reported for DNA-based microorganisms with 3.4×10−3 spontaneous mutations per genome, per replication, suggesting that DNA mismatches resulting from DNA replication errors are actively corrected in those organisms [12]. Surprisingly, only 11 out of the 54 archaeal genomes sequenced so far encode for homologs of the conserved MutS1 protein subfamily found in Bacteria and Eukarya [13]. Archaea with MutS1 homologs include halophiles and methanogens, all members of the domain Euryarchaeota. These archaeal MutS1 proteins share identical domain structure with their bacterial counterparts and are likely the result of a lateral gene transfer event [13]. Also detected in the Archaea are MutS2-like proteins [13]. In the hyperthermophilic archaeon Pyrococcus furiosus, MutS2 has been shown to have ATPase and DNA binding activity but no specific DNA mismatch binding activity [14]. Proteins from the MutS2 subfamily, which are not thought to be involved in MMR, are different in structure and sequence to the MutS1 subfamily proteins except for the MutSAc domain essential for dimerization, ATPase and DNA binding activities [13].
Halobacterium salinarum NRC-1 is an extremely halophilic archaeon growing optimally in 4M NaCl [15], [16], [17]. The high osmotic pressure from its environment is counterbalanced by a 4M intracellular concentration of KCl [18]. Previous studies revealed the exceptional resistance of H. salinarum NRC-1 to desiccation, UV and ionizing radiation, which was attributed to efficient DNA repair and detoxification systems and to its adaptation to hypersaline environments, characterized by high levels of solar radiation and periodic desiccation [19], [20], [21]. The genome of H. salinarum NRC-1 has been sequenced [22] and encodes for proteins of conserved DNA repair pathways that include damage reversal, base excision repair (BER), nucleotide excision repair (NER), homologous recombination, and the bacterial-like MMR proteins MutS and MutL [22]. We demonstrated that the eukaryal-like homologous recombination protein, Mre11, is essential for the repair of DNA double strand breaks in H. salinarum NRC-1, whereas Rad50 is dispensable, representing a shift from the eukaryotic model of recombinational repair [23]. Crowley et al. [24] showed that the bacterial NER homologs UvrA/B/C encoded in the genome of H. salinarum NRC-1 were essential for the survival of the organism to UV irradiation. These studies demonstrate the mosaic nature of the DNA repair pathways in H. salinarum NRC-1, and in the Archaea in general, and raise questions about the nature of DNA MMR in this organism. Through computational analysis we found that H. salinarum NRC-1 has three bacterial-like mutS genes, a bacterial-like mutL gene, 4 bacterial-like recJ exonuclease genes, 1 eukaryotic-like rad2 5′-3′ exonuclease gene, and a bacterial-like uvrD helicase gene, all potentially involved in MMR. Two of the MutS proteins in H. salinarum NRC-1, MutS1 and MutS2, are homologous to the MutS1 protein subfamily, and have been renamed MutS1A and MutS1B in this study, while the third MutS protein, MutS3, is homologous to proteins found in the MutS2 subfamily [13]. Whole-genome transcriptomic studies conducted on H. salinarum cells exposed to UV and gamma radiation, and to oxidative stress revealed no significant changes in mRNA level for mutS1A, mutS1B, and mutL when compared to untreated cells [20], [21, unpublished].
Here, we used a genetic approach to determine the spontaneous genomic mutation rate in H. salinarum NRC-1 and to determine the cellular role of the bacterial-like MMR proteins MutS and MutL encoded in its genome. Our analysis, using fluctuation tests targeting genes of the uridine monophosphate (UMP) biosynthesis pathway, revealed a genomic mutation rate similar to that of other DNA-based microorganisms and a markedly different spectrum of mutational changes. The phenotypic analysis of deletion mutants for the mutL, mutS1A, mutS1B, and uvrD genes and a mutS1A/mutS1B double mutant showed little difference between the mutant and background strains indicating that the MutS and MutL protein homologs found in H. salinarum NRC-1 are not essential for maintaining the low incidence of spontaneous mutations observed in this organism.
Results
Genomic Mutation Rate
We calculated the spontaneous genomic mutation rate of H. salinarum NRC-1 to determine the replication fidelity in this mesophilic archaeon. We performed six independent fluctuation tests [12], targeting forward mutations in genes of the UMP biosynthetic pathway producing 5-fluoroorotic acid (5-FOA) resistant mutants. The mutation rate was calculated using the equation μ = ln(m/Nav) [25] and was found to be 3.73×10−7+/−1.44×10−7 mutations per replication. Sequencing of purified 5-FOA-resistant mutants revealed that only 13 out of 55 sequenced mutants had a mutation in the pyrF gene (orotate decarboxylase), only 42 out of 69 had a mutation in the pyrE2 gene, and none out of 61 had mutation in the pyrE1 gene (orotate phosphoribosyl transferases). We therefore adjusted the gene mutation rates for the pyrF and pyrE2 genes by factors of 0.24 (13/55) and 0.60 (42/69), respectively, and did not use the pyrE1 gene in our calculation. The resulting spontaneous mutation rates at the gene level were 8.95×10−8 and 2.24×10−7 mutations per gene per replication, for the pyrF and pyrE2 genes, respectively. To correct this rate for the fraction of undetected mutations producing no phenotypic effect, we adjusted the total number of base pair substitutions (BPS) using published information on BPS detection efficiency (approximately 0.2) [4]. The resulting rate estimate per gene was calculated as follows:This rate was then converted into a genomic rate by dividing by gene size (pyrF = 803bp, and pyrE2 = 527bp) and multiplying by genome size (2,571,010bp), resulting in an average genomic mutation rate, corrected for undetected mutations, for H. salinarum NRC-1 of 1.67×10−3±1.4×10−3 mutation per replication (0.62×10−3 for pyrF and 2.7×10−3 for pyrE2).
Mutational Spectrum
One hundred forty-nine 5-FOA-resistant mutants were recovered from two fluctuation tests and of those 55 were sequenced using primers for pyrF, 61 using primers for pyrE1, and 65 using primers for pyrE2. Mutations were only found in the coding regions of the pyrF and pyrE2 genes and a number of mutants had mutations in more than one gene (Table 1). No mutation was found within a 100-bp region upstream of either gene. A disproportionate number of deletions were found in the pyrE2 gene when compared to the pyrF gene (Figure 1, Table 2). All of the pyrE2 mutations occurred at a single hotspot in the gene at positions 382–397, with a 7-nucleotide (nt) deletion (GTCGACG) found in 42 mutants; two 1-nt deletions and two single BPS also contributed to the changes observed in this gene (Table 2, Figure 2A). A 9-nt sequence found at the pyrE2 gene mutational hotspot was the direct repeat of a sequence located immediately upstream (Figure 2A). Mutations were distributed throughout the pyrF gene with a concentration of insertions at position 354–365 (Table 2, Figure 2B). None of the mutations were the result of a transposon element insertion. Indels out numbered BPS and approximately 80% of the BPS resulted in non-synonymous amino acid changes. BPS found in the pyrF and pyrE2 genes were mostly transversions (80%) (Table 2).
Insertions, deletions and base pair substitutions (BPS) in the pyrF and pyrE2 genes were obtained by sequencing 5-FOA-resistant uracil auxotrophs of H. salinarum NRC-1.
Mutations were identified by sequencing 5-FOA-resistant mutants. One-letter code for amino acid is under the gene nucleotide sequence; start and stop codons are boxed with solid lines; putative TATA box is boxed with a dotted line; ▿ indicates insertion of the base(s) specified in parenthesis next to the symbol; ▵ indicates deletion of bases directly located below the symbol; BPS changes are indicated above the sequence in bold; highlighted in green and yellow in pyrE2 are the two 9-nt direct repeats. (A) pyrE2 gene and (B) pyrF gene.
MutS and MutL Are Not Essential for the Low Incidence of Mutation Observed in H. salinarum NRC-1
To determine whether the bacterial-like MMR proteins encoded in the genome of H. salinarum NRC-1 were essential in maintaining the low genomic mutation rate we observed, we carried out targeted gene deletions of the mutS1A, mutS1B, mutL, and uvrD genes, with a double deletion of the mutS1A and mutS1B genes, using the background strains AK07 (Δura3) and CB08 (Δura3Δzim) [26], [27]. Genotypes of mutant strains CB071 (ΔmutS1A), CB072 (ΔmutS1B), CB073 (ΔmutS1AΔmutS1B), CB074 (ΔmutL), and CB081 (ΔuvrD) were confirmed by Southern blot hybridization after initial screening by PCR (Figure 3). Phenotypic characterization of the mutant strains revealed no growth defects at 37, 42, and 45°C when compared to the background strain AK07 (Δura3). While tolerance to alkylating agents is a hallmark of bacterial and eukaryal MMR systems, we found no increased tolerance to alkylation when the mutants were exposed to 600mg/L N-Methyl-N′-Nitro-N-Nitrosoguanidine (MNNG) (Figure 4). Mismatch repair deletion mutants were also constructed in another background strain, CB08 (Δura3Δzim), to confirm the validity of our gene deletions. All mismatch deletion mutants in the Δura3Δzim background showed the same phenotypes as the mutants constructed in the Δura3 background, which together with PCR and Southern blot analyses established that we constructed true gene deletions.
Probes for Southern blot analysis were designed to hybridize to regions 500 nt downstream of the target genes coding region. PCR analysis primers were located 500 nt upstream of the start codon of the targeted gene and 1000 nt downstream of the stop codon. (A) Southern hybridizations. (B) PCR analysis: positive lanes template was wildtype H. salinarum NRC-1 DNA, negative lanes had no template DNA.
H. salinarum NRC-1 background strain Δura3 (AK07) and mutant strains ΔmutL (CB074), ΔmutS1A (CB071), ΔmutS1B (CB072), ΔmutS1AΔmutS1B (CB073), and ΔuvrD (CB081) were exposed to 50, 100, and 600mg/L of MNNG. Survival was calculated as the average ratio (N/No) of surviving CFU from treated cultures (N) and untreated (No) cultures. Data are the average of a least three independent experiments, with standard errors shown.
Significant increases in mutation rate following deletion of mismatch repair genes has been demonstrated in both Bacteria and Eukarya [1], [2], [3]. Here we designed a 5-FOA mutation frequency assay to test whether our H. salinarum NRC-1 MMR deficient strains showed increased mutation rates. This assay, targeting forward mutations in a plasmid copy of the pyrF (ura3) gene, allowed us to compare the mutation frequencies for each of the mutant strains tested with that of the background strain. The use of a plasmid assay was born out of necessity because our MMR deletions were constructed in a pyrF deletion background; the limited genetic markers available for H. salinarum NRC-1 constrained our ability to restore the pyrF gene into the genome of H. salinarum NRC-1. Using fluctuation tests, we obtained comparable mutation frequencies (ratio of 5-FOA-resistant colonies to the average number of colonies plated) for the background and deletion mutant strains (Table 3). Sequencing of 20 mutants for strain AK07 (Δura3), 22 mutants for strain CB071 (ΔmutS1A), 19 mutants for strain CB073 (ΔmutS1AΔmutS1B), and 33 mutants for strain CB074 (ΔmutL) showed that approximately 70% of the mutants had changes in the plasmid copy of the pyrF gene and that the distribution of BPS, insertions and deletions in the pyrF gene of the MMR mutants was similar to that found for the control strain AK07 (Table 3). We found that all the changes in the plasmid copies of the pyrF gene were BPS, with the exception of the insertion of a T in one of the CB071 (ΔmutS1A) mutants, and that most BPS were C/G or G/C transversions (55 to 75%).
Discussion
The MMR pathway is essential for maintaining genome stability by correcting errors introduced by DNA polymerases during DNA replication [1], [2], [3]. We determined the genomic mutation rate of the halophilic mesophile, H. salinarum NRC-1, one of the few archaea to encode for a bacterial-like version of the conserved MutS and MutL proteins. We found that H. salinarum NRC-1 mutation rate was similar to that previously calculated for DNA-based microorganisms, with 1.67×10−3 mutations per genome per replication, suggesting the presence of high fidelity replication machinery in this organism [4], [12], [28], [29], [30], [31], [32]. The DNA-based organisms that formed the basis for this genomic mutation rate comparison are all mesophiles, belong to the Bacteria and the Eukarya domains, and include several bacteriophages. When we compared the mutation rate of H. salinarum NRC-1 with that of organisms adapted to high temperature, we found that it was almost twice that of the acidophilic archaeon, S. acidocaldarius, and more than an order of magnitude higher than that of the bacterium Thermus thermophilus, with 1.8×10−3 and 9.7×10−4 mutations per genome per replication, respectively [4], [33]. While the low genomic mutation rate found in thermophiles might be an adaption to extreme temperature conditions [33], the mutation rate of the mesophilic halophile, H. volcanii was also found to be extremely low with 4.5×10−4 mutations per genome per replication [5]. It is possible that this extremely low genomic mutation rate is the result of phenotypic lag or the effects of H. volcanii polyploidy [5]. However, H. salinarum NRC-1, similarly to H. volcanii, has approximately 20 copies of its chromosome per cell [5], [34] suggesting that high gene redundancy does not necessarily result in mutations not being efficiently detected if they occur at low copy number.
H. volcanii has an unusual mutational spectrum with a prevalence of in-frame indels flanked by direct repeats [5]. In H. salinarum NRC-1, a 9-nt direct repeat in the pyrE2 gene resulted in a mutational hotspot that strongly biased the mutational spectrum of this gene toward deletions. Short direct repeats, tandem repeats and monotonic runs are known to promote strand misalignments during DNA replication and might explain the increased deletion frequency we observed in the H. salinarum NRC-1 pyrE2 gene [35]. An alternative explanation might be that this 7-nt deletion and the 4 other mutations found in the pyrE2 gene – shared by 42 of the mutants - were a pre-existing set of mutations present in the initial cultures used for the fluctuation tests, which were a subset of the 6 fluctuation tests performed to calculate the genomic mutation rate of H. salinarum NRC-1. With the exception of the pyrE2 gene hotspot most of the indels found in H. salinarum NRC-1 were ±1 frameshifts, which were also found to be predominant in mesophiles and in the two thermophilic organisms investigated so far [4], [33]. A low overall number of BPS was found in H. salinarum NRC-1, in the two archaea, S. acidocaldarius and H. volcanii, and in the thermophilic bacterium, T. thermophilus [4], [5], [33]. This is a departure from data reported for mesophilic bacteria and eukaryotes where BPS constituted the overwhelming majority, 60–80%, of mutations [4], [12]. The low occurrence of BPS in Archaea might be related to the high GC content (>60%) found in both H. salinarum NRC-1 and H. volcanii, which results in long runs of G and C nucleotides [22]. During replication these GC runs may cause polymerase slippage leading to the formation of indels. However, S. acidocaldarius only has a 36% GC content arguing against this idea [36]. An alternative explanation for the low BPS found in these archaea might be a consequence of their adaptations to extreme environmental conditions – high temperature and high salt – resulting in intrinsic properties of their replication machinery [37], [38], [39]. This argument is strengthened by the recent finding of a low mutation rate and a low incidence of BPS in the thermophilic bacterium, T. thermophilus [33].
Phenotypic analyses of MMR deletion mutants in H. salinarum NRC-1 did not show an increase in alkylation tolerance to MNNG. This was in contrast to previous studies showing that increased alkylation tolerance is a hallmark of MMR systems in most bacteria and eukarya [2], [40]. However, in Saccharomyces cerevisiae no tolerance to MNNG was demonstrated in MMR deficient yeast strains unless the MGT1 methyltransferase, responsible for correcting O6-methyl guanine damage, was also absent [41]. The H. salinarum NRC-1 genome does not encode for a yeast MGT1 methyltransferase homolog but this does not rule out the repair of O6-methyl guanine lesions by another glycosylase yet to be identified.
An increased mutation rate following deletion of mismatch repair genes has been demonstrated both in Bacteria and Eukarya. In D. radiodurans and E. coli, the mutation rate was calculated in cells deficient in MutS, MutL, or UvrD proteins and was found to increase 7 to 1,000-fold compared to wildtype cells [42], [43]. In E. coli lacking MutH, MutS, and MutL proteins, forward mutation studies in the lacI gene showed a 200-fold increase in mutation rate [32], [44], [45]. In S. cerevisiae, both forward and reverse mutation rate studies have shown that the MutS and MutL homologs are required for base correction [1], [2]. In contrast to those studies, the mutation frequencies obtained from our MMR deletion mutants showed that the bacterial-like MMR pathway encoded in the genome of H. salinarum NRC-1 is not essential for maintaining the low genomic mutation rate we observed in this organism. We found that the mutational spectrum of H. salinarum NRC-1 MMR mutants is identical to that of the control strain, indicating that if functional, the bacterial-like MMR system does not have a major role in this organism. This finding is a departure from the reported role of the bacterial-like NER pathway in H. salinarum NRC-1, where the uvrA, uvrB or uvrC genes were required for the repair of photoproducts induced by UV light [24].
No homologs of the conserved mismatch repair genes, mutS and mutL, have been found in the genome of S. acidocaldarius despite its low genomic mutation rate, leading to the puzzling question of what is responsible for the maintenance of the low genomic mutation rate observed in this organism and in the Archaea as a whole [4], [5], [46]. One explanation might come from the adaptive mechanisms these organisms have evolved to thrive in extreme environments. Studies calculating replication fidelity showed that the polymerase from P. furiosus has 10-fold higher replication fidelity than the polymerase from T. aquaticus, a thermophilic bacterium, and 5-fold higher replication fidelity than E. coli DNA polymerase III holoenzyme [37], [38], [39], [47], [48]. The replicative polymerases in the Archaea are members of the B-family and are more similar to eukaryotic polymerases than to bacterial ones [49]. The higher fidelity of archaeal polymerases, resulting from differences in sequence and structure, could play a role in maintaining the genomic integrity of the Archaea. An alternative hypothesis might be that MMR is carried out by DNA repair enzymes specifically recruited to this function rather than by a canonical MMR pathway [50]. For example, deamination of cytosine to uracil in a GC base pair is a major mutagenic event that generates a G.C to A.T mutations [51]. Polymerases in the Archaea possess the unique ability to stall when a uracil residue is encountered [52], the uracil is then removed by an uracil-DNA glycosylase [52], [53]. Furthermore, direct interaction between uracil-DNA glycosylase and a PCNA homolog from Pyrobaculum aerophilum has been documented [54], [55], [56], suggesting the possibility of recruitment of DNA glycosylases, or other DNA repair proteins, to damage sites by the replication machinery.
Further studies, in particular to elucidate if there is any MMR activity in the protein complement of H. salinarum NRC-1, are necessary to differentiate if the low mutation rate observed in this organism is the result of intrinsic properties of its DNA polymerase or if DNA repair activities are recruited to the replication site to remove mismatches. The presence of bacterial-like MMR genes in several archaea is thought to be the result of lateral gene transfer from bacteria [13], [20] and, although we showed here that the mutS and mutL genes are not essential for the maintenance of a low mutation rate in H. salinarum NRC-1, it does not mean that this pathway is not functional, existing possibly as an alternate or minor pathway. Until the mechanism(s) maintaining genomic stability in the archaea are elucidated, this hypothesis remains difficult to test.
Materials and Methods
Organism and Growth Conditions
Halobacterium salinarum NRC-1 (ATCC 700922) was grown in standard GN101 medium (250g/L NaCl, 20g/L MgSO4, 2g/L KCl, 3g/L sodium citrate, 10g/L Oxoid brand bacteriological peptone), pH 7.2, with the addition of 1 mL/L trace elements solution (31.5mg/L FeSO4·7H2O, 4.4mg/L ZnSO4·7H2O, 3.3mg/L MnSO4·H2O, 0.1mg/L CuSO4·5H2O) at 42°C shaking in a Gyromax 737 shaker (Amerex Instruments; LaFayette, CA) at 220rpm. When specified, the GN101 media was supplemented with 50mg/L uracil and 350mg/L 5-fluoroorotic acid (5-FOA) (Sigma; St. Louis, MO), final concentrations. Basal salts solution (BSS), of the same composition as GN101 without the peptone, was used for culture dilutions.
Targeted Gene Deletion
In-frame gene deletions, ΔmutL (CB074), ΔmutS1A (CB071), ΔmutS1B (CB072), ΔmutS1AΔmutS1B (CB073), and ΔuvrD (CB081), were constructed as described previously [26], [27]. In short, uracil dropout medium was used for selection for uracil prototrophy following transformation of strain AK07 (Δura3) [23] and strain CB08 (Δura3Δzim) (DiRuggiero et al. unpublished) with plasmid pNBK07 (from M.P. Krebs, Illinois State University, Normal, Illinois) bearing the knockout gene constructs and the ura3 marker for uracil biosynthesis. The zim gene was originally thought to have a function similar to the dam methylase from E. coli. Further experiments showed that it is a putative CTAG methylase from an uncharacterized restriction and modification system. 5-FOA was subsequently used to select intra-molecular recombinants that lost the plasmid. Recombinants were screened by PCR and confirmed by Southern blot analysis. The GN101 medium was supplemented with 50mg/L uracil for all Δura3 strains.
Mutation Rate Assay and Mutant Isolation
Fluctuation tests were performed as previously described [4], [57] with modifications. Tubes of GN101 medium supplemented with 50mg/L uracil and inoculated with one colony of H. salinarum NRC-1 each were incubated at 42°C, with shaking, for two days. The resulting cultures were diluted to 1×102 cells/mL, 150µL aliquots were dispensed into 60 wells of a 96 well flat-bottomed plate (VWR; West Chester, PA), and incubated for 3 days at 42°C without shaking until cell density reached approximately 1×105 cells/mL. The entire contents from each well were spread on GN101 plates with 50mg/L uracil and 350mg/L 5-FOA and colonies were counted 11 days after incubation at 42°C. The spontaneous mutation rate for H. salinarum NRC-1 was calculated using the relationship μ = ln(m/Nav) [4], [5]. The mutation rate, μ, is equal to the natural log of m, number of mutational events per culture, divided by Nav, average number of cells per culture. The m value was calculated using the MSS Maximum-Likelihood method as previously described [25]. Colonies from the 5-FOA plates were randomly selected, clonally purified by restreaking twice onto GN101 plates with 50mg/L uracil and 350mg/L 5-FOA and stored at −80°C for sequence analysis.
Sequence Analysis
The UMP biosynthesis genes, pyrE1, pyrE2, and pyrF, were amplified by colony PCR using clonally purified 5-FOA resistant mutants as described above and the following primers: pyrE1-F (5′CCTCGTCCTGGAGAACAAAG3′), pyrE1-R (5′ATCGAAGGCCATGTCCCACCGT3′), pyrE2-F (5′GGTTCATACCGACCACACG3′), pyrE2-R (5′TCGGCGACACCTTCGGGCTG3′), pyrF-F (5′GCGCGCCTCGTGGTGTTCGT3′), and pyrF-R (5′AGCGTCGTCTGTGACACCCA3′). Primers were located 200 bp upstream and 100 downstream of the coding region for these genes. PCR conditions were 2 minutes at 94°C, 30 cycles of 35 seconds at 94°C, 40 seconds at 56°C, and 60 seconds at 72°C, and a final 5 minutes at 72°C using FastTaq DNA polymerase (Roche; Indianapolis, IN). PCR products were purified using ExoSap-IT (GE Healthcare; Piscataway, NJ) and sequencing was performed using an Applied Biosystems 3730xl DNA Analyzer with BigDye Terminator reagents (Applied Biosystems; Foster City, CA).
Mutation Frequency Assay by the 5-FOA Selection Method
The ura3 (pyrF) gene and its native promoter were PCR amplified from H. salinarum NRC-1 genomic DNA, cloned into plasmid pNBPA, and the resulting construct was transformed into strain AK07 (Δura3) and into the deletion mutant strains CB071 (ΔmutS1A), CB073 (ΔmutS1AΔmutS1B), and CB074 (ΔmutL). Fluctuation tests were performed as described above and cultures were plated in toto on GN101 medium supplemented with 50mg/L uracil, 50µM mevinolin, and 350mg/L 5-FOA. Colonies were screened by PCR to ensure presence of the ura3 gene on plasmid pNBPA and the ura3 gene was sequenced using primers: ura3-Acc65I (5′GCGGGTACCGTCGGCTGGCGGGCACGCGGT3′) and ura3-SpeI (5′GCGACTAGTCTACCGGTGGCGGTTCAGGCG3′). Mutation frequency was calculated as the ratio of 5-FOA-resistant colonies to the average number of colonies plated.
Survival Assays
Temperature growth defect assay.
Single colonies of wildtype H. salinarum NRC-1 and deletion mutant strains were grown to mid-log phase in GN101 medium supplemented with 50mg/L uracil. Cultures were diluted to 104, 103, and 102 cells/mL and 5µL of each spotted in triplicate on GN101 medium supplemented with 50mg/L uracil. Plates were observed for growth at 37°C, 42°C, and 45°C after 7 days of incubation.
N-Methyl-N′-Nitro-N-Nitrosoguanidine (MNNG) survival assay.
cultures were grown to OD600 0.6 in GN101 medium supplemented with 50mg/L uracil and diluted to OD600 0.4. Cultures were divided into 5mL aliquots, treated with the addition of 0 or 600mg/L MNNG, in triplicate, and incubated in the dark at 42°C with shaking for 1 hour. Cells were washed, diluted in BSS, and plated on GN101 medium supplemented with 50mg/mL uracil. Plates were incubated at 42°C in the dark for 7–10 days and colonies counted. Survival was calculated as N/No where N is the number of viable cells after MNNG treatment and No is the number of viable cells without treatment. MNNG (TimTech; New Zealand) stock solution was made at 50g/L in DMSO.
Author Contributions
Conceived and designed the experiments: CRB JD. Performed the experiments: CRB. Analyzed the data: CRB JD. Wrote the paper: JD.
References
- 1. Harfe BD, Jinks-Robertson S (2000) DNA mismatch repair and genetic instability. Annu Rev Genet 34: 359–399.
- 2. Schofield MJ, Hsieh P (2003) DNA mismatch repair: molecular mechanisms and biological function. Annu Rev Microbiol 57: 579–608.
- 3. Iyer RR, Pluciennik A, Burdett V, Modrich PL (2006) DNA mismatch repair: functions and mechanisms. Chem Rev 106: 302–323.
- 4. Grogan DW, Carver GT, Drake JW (2001) Genetic fidelity under harsh conditions: Analysis of spontaneous mutation in the thermoacidophilic archaeon Sulfolobus acidocaldarius. Proc Natl Acad Sci U S A 98: 7928–7933.
- 5. Mackwan RR, Carver GT, Drake JW, Grogan DW (2007) An unusual pattern of spontaneous mutations recovered in the halophilic archaeon Haloferax volcanii. Genetics 176: 697–702.
- 6. Bjornson KP, Blackwell LJ, Sage H, Baitinger C, Allen D, et al. (2003) Assembly and molecular activities of the MutS tetramer. J Biol Chem 278: 34667–34673.
- 7. Obmolova G, Ban C, Hsieh P, Yang W (2000) Crystal structures of mismatch repair protein MutS and its complex with a substrate DNA. Nature 407: 703–710.
- 8. Mendillo ML, Putnam CD, Polodner RD (2007) Escherichia coli MutS tetramerization domain structure reveals that stable dimers but not tetramers are essential for DNA mismatch repair in vivo. J Biol Chem 282: 16345–16354.
- 9. Acharya S, Foster PL, Brooks P, Fishel R (2003) The coordinated functions of the E. coli MutS and MutL proteins in mismatch repair. Mol Cell 12: 233–246.
- 10. Bergerat A, Demassy B, Gadelle D, Varoutas PC, Nicolas A, et al. (1997) An atypical topoisomerase II from archaea with implications for meiotic recombination. Nature 386: 414–417.
- 11. Umar A, Buermeyer AB, Simon JA, Thomas DC, Clark AB, et al. (1996) Requirement for PCNA in DNA mismatch repair at a step preceding DNA resynthesis. Cell 87: 65–73.
- 12. Drake J (1991) A constant rate of spontaneous mutation in DNA-based microbes. Proc Natl Acad Sci U S A 88: 7160–7164.
- 13. Lin Z, Nei M, Ma H (2007) The origins and early evolution of DNA mismatch repair genes multiple horizontal gene transfers and co-evolution. Nucl Acids Res 35: 7591–7603.
- 14. Vijayvargia R, Biswas I (2002) MutS2 family protein from Pyrococcus furiosus. Curr Microbiol 44: 224–228.
- 15. McCready S, Muller JA, Boubriak I, Berquist BR, Ng WL, et al. (2005) UV irradiation induces homologous recombination genes in the model archaeon, Halobacterium sp. NRC-1. Saline Systems 1: 3.
- 16. Martin EL, Reinhardt RL, Baum LL, Becker MR, Shaffer JJ, et al. (2000) The effects of ultraviolet radiation on the moderate halophile Halomonas elongata and the extreme halophile Halobacterium salinarum. Can J Microbiol 46: 180–187.
- 17. Potts M (1994) Desiccation tolerance of prokaryotes. Microbiol Rev 58: 755–805.
- 18. Engel MB, Catchpole HR (2005) A microprobe analysis of inorganic elements in Halobacterium salinarum. Cell Biol International 29: 616–622.
- 19. Kottemann M, Kish A, Iloanusi C, Bjork S, DiRuggiero J (2005) Physiological responses of the halophilic archaeon Halobacterium sp. strain NRC1 to desiccation and gamma irradiation. Extremophiles 9: 219–227.
- 20. Baliga NS, Bjork SJ, Bonneau R, Pan M, Iloanusi C, et al. (2004) Systems level insights into the stress response to UV radiation in the halophilic archaeon Halobacterium NRC-1. Genome Res 14: 1025–1035.
- 21. Whitehead K, Kish A, Pan M, Kaur A, Reiss DJ, et al. (2006) An integrated systems approach for understanding cellular responses to gamma radiation. Mol Syst Biol 2: 47.
- 22. Ng WV, Kennedy SP, Mahairas GG, Berquist B, Pan M, et al. (2000) Genome sequence of Halobacterium species NRC-1. Proc Natl Acad Sci U S A 97: 12176–12181.
- 23. Kish A, DiRuggiero J (2008) Rad50 is not essential for the Mre11-dependant repair of DNA double strand breaks in Halobacterium sp. str. NRC-1. J Bacteriol 190: 5210–5216.
- 24. Crowley DJ, Boubriak I, Berquist BR, Clark M, Richard E, et al. (2006) The uvrA, uvrB and uvrC genes are required for repair of ultraviolet light induced DNA photoproducts in Halobacterium sp. NRC-1. Saline Systems 2: 11.
- 25. Rosche WA, Foster PL (2000) Determining mutation rates in bacterial populations. Methods 20: 4–17.
- 26. Peck RF, Dassarma S, Krebs MP (2000) Homologous gene knockout in the archaeon Halobacterium salinarum with ura3 as a counterselectable marker. Mol Microbiol 35: 667–676.
- 27. Wang G, Kennedy SP, Fasiludeen S, Rensing C, DasSarma S (2004) Arsenic resistance in Halobacterium sp. strain NRC-1 examined by using an improved gene knockout system. J Bacteriol 186: 3187–3194.
- 28. Schaaper RM, Danforth BN, Glickman BW (1986) Mechanisms of spontaneous mutagenesis: an analysis of the spectrum of spontaneous mutation in the Escherichia coli lacI gene. J Mol Biol 189: 273–284.
- 29. Drake JW, Charlesworth B, Charlesworth D, Crow JF (1998) Rates of spontaneous mutation. Genetics 148: 1667–1686.
- 30. Halliday JA, Glickman BW (1991) Mechanisms of spontaneous mutation in DNA repair-proficient Escherichia coli. Mutat Res 250: 55–71.
- 31. Farabaugh PJ, Fchmeissner U, Hofer M, Miller JH (1978) Genetic studies of the lac repressor. VII. On the molecular nature of spontaneous hotspots in the lacI gene of Escherichia coli. J Mol Biol 126: 847–857.
- 32. Schaaper RM, Dunn RL (1991) Spontaneous mutation in the Escherichia coli lacI gene. Genetics 129: 317–326.
- 33. Mackwan RR, Carver GT, Kissling GE, Drake JW, Grogan DW (2008) The rate and character of spontaneous mutation in Thermus thermophilus. Genetics 180: 17–25.
- 34. Breuert S, Allers T, Spohn G, Soppa J (2006) Regulated polyploidy in halophilic Archaea. PLoS ONE 1: e92.
- 35. Ripley LS (1990) Frameshift mutation: determinants of specificity. Annu Rev Genet 23: 189–213.
- 36. Chen L, Brugger K, Skovgaard M, Redder P, She Q, et al. (2005) The genome of Sulfolobus acidocaldarius, a model organism of the Crenarchaeota. J Bacteriol 187: 4992–4999.
- 37. Cline J, Braman JC, Hogrefe HH (1996) PCR fidelity of Pfu DNA polymerase and other thermostable DNA polymerases. Nucleic Acids Res 24: 3546–3551.
- 38. Lundberg KS, Shoemaker DD, Adams MWW, Short JM, Sorge JA, et al. (1991) High-fidelity amplification using a thermostable DNA polymerase isolated from Pyrococcus furiosus. Gene 108: 1–6.
- 39. Mattila P, Korpela J, Tenkanen T, Pitkanen K (1991) Fidelity of DNA synthesis by the Thermococcus litoralis DNA polymerase - an extremely heat stable enzyme with proofreading activity. Nucleic Acids Res 19: 4967–4973.
- 40. Bawa S, Xiao W (1997) A mutation in the MSH5 gene results in alkylation tolerance. Cancer Res 57: 2715–2720.
- 41. Xiao W, Rathgeber L, Fontanie T, Bawa S (1995) DNA mismatch repair mutants do not increase N-methyl-N′-nitro-N-nitrosoguanidine tolerance in O6-methylguanine DNA methyltransferase-deficient yeast cells. Carcinogenesis 16: 1933–1939.
- 42. Mennecier S, Coste G, Servant P, Bailone A, Sommer S (2004) Mismatch repair ensures fidelity of replication and recombination in the radioresistant organism Deinococcus radiodurans. Mol Genet Genomics 272: 460–469. Epub.
- 43. LeClerc JE, Li B, Payne WL, Cebula TA (1996) High mutation frequencies among Escherichia coli and Salmonella pathogens. Science 274: 1208–1211.
- 44. Schaaper RM, Dunn RL (1987) Spectra of spontaneous mutations in Escherichia coli strains defective in mismatch correction: the nature of in vivo DNA replication errors. Proc Natl Acad Sci U S A 84: 6220–6224.
- 45. Schaaper RM (1993) Base selection, proofreading, and mismatch repair during DNA replication in Escherichia coli. J Biol Chem 268: 23762–23765.
- 46. Berkner S, Lipps G (2008) Mutation and reversion frequencies of different Sulfolobus species and strains. Extremophiles 12: 263–270.
- 47. Bloom LB, Chen X, Fygenson DK, Turner J, O'Donnell M, et al. (1997) Fidelity of Escherichia coli DNA polymerase III holoenzyme. The effects of beta, gamma complex processivity proteins and epsilon proofreading exonuclease on nucleotide misincorporation efficiencies. J Biol Chem 272: 27919–27930.
- 48. Kunkel TA (1992) DNA replication fidelity. J Biol Chem 267: 18251–18254.
- 49. Barry ER, Bell SD (2006) DNA replication in the Archaea. Microbiol Mol Biol Rev 70: 876–887.
- 50. Grogan DW (2000) The question of DNA repair in hyperthermophilic archaea. Trends Microbiol 8: 180–185.
- 51. Zharkov DO (2008) Base excision DNA repair. Cell Mol Life Sci 65: 1544–1565.
- 52. Lasken RS, Schuster DM, Rashtchian A (1996) Archaebacterial DNA polymerases tightly bind uracil-containing DNA. J Biol Chem 271: 17692–17696.
- 53. Yang H, Fitz-Gibbon S, Marcotte EM, Tai JH, Hyman EC, et al. (2000) Characterization of a thermostable DNA Glycolase specific for U/G and T/G mismatches from the hyperthermophilic archeon Pyrobaculum aerophilum. J Bacteriol 182: 1272–1279.
- 54. Yang H, Chiang J-H, Fitz-Gibbon S, Lebel M, Sartori AA, et al. (2002) Direct interaction between uracil-DNA glycosylase and a PCNA homolog in the crenarchaeon Pyrobaculum aerophilum. J Biol Chem M201820200.
- 55. Connolly BA, Fogg MJ, Shuttleworth G, Wilson BT (2003) Uracil recognition by archaeal family B DNA polymerases. Biochem Soc Trans 31: 699–702.
- 56. Fogg MJ, Pearl LH, Connolly BA (2002) Structural basis for uracil recognition by archaeal family B DNA polymerases. Nat Struct Biol 9: 922–927.
- 57. Jacobs KL, Grogan DW (1997) Rates of spontaneous mutation in an archaeon from geothermal environments. J Bacteriol 179: 3298–3303.