The rate and character of spontaneous mutation in Thermus thermophilus

Genomic attributes of thermophilic and hyperthermophilic bacteria and archaea

Thermophiles and hyperthermophiles are immensely useful in understanding the evolution of life, besides their utility in environmental and industrial biotechnology. Advancements in sequencing technologies have revolutionized the field of microbial genomics. The massive generation of data enhances the sequencing coverage multi-fold and allows to analyse the entire genomic features of microbes efficiently and accurately. The mandate of a pure isolate can also be bypassed where whole metagenome-assembled genomes and single cell-based sequencing have fulfilled the majority of the criteria to decode various attributes of microbial genomes. A boom has, therefore, been seen in analysing the extremophilic bacteria and archaea using sequence-based approaches. Due to extensive sequence analysis, it becomes easier to understand the gene flow and their evolution among the members of bacteria and archaea. For instance, sequencing unveiled that Thermotoga maritima shares around 24% of genes of archaeal origin. Comparative and functional genomics provide an analytical view to understanding the microbial diversity of thermophilic bacteria and archaea, their interactions with other microbes, their adaptations, gene flow, and evolution over time. In this review, the genomic features of thermophilic bacteria and archaea are dealt with comprehensively.

Unexpectedly high mutation rate of a deep-sea hyperthermophilic anaerobic archaeon

Deep-sea hydrothermal vents resemble the early Earth, and thus the dominant Thermococcaceae inhabitants, which occupy an evolutionarily basal position of the archaeal tree and take an obligate anaerobic hyperthermophilic free-living lifestyle, are likely excellent models to study the evolution of early life. Here, we determined that unbiased mutation rate of a representative species, Thermococcus eurythermalis, exceeded that of all known free-living prokaryotes by 1-2 orders of magnitude, and thus rejected the long-standing hypothesis that low mutation rates were selectively favored in hyperthermophiles. We further sequenced multiple and diverse isolates of this species and calculated that T. eurythermalis has a lower effective population size than other free-living prokaryotes by 1-2 orders of magnitude. These data collectively indicate that the high mutation rate of this species is not selectively favored but instead driven by random genetic drift. The availability of these unusual data also helps explore mechanisms underlying microbial genome size evolution. We showed that genome size is negatively correlated with mutation rate and positively correlated with effective population size across 30 bacterial and archaeal lineages, suggesting that increased mutation rate and random genetic drift are likely two important mechanisms driving microbial genome reduction. Future determinations of the unbiased mutation rate of more representative lineages with highly reduced genomes such as Prochlorococcus and Pelagibacterales that dominate marine microbial communities are essential to test these hypotheses.

Estimation of the Genome-Wide Mutation Rate and Spectrum in the Archaeal Species Haloferax volcanii

Organisms adapted to life in extreme habitats (extremophiles) can further our understanding of the mechanisms of genetic stability, particularly replication and repair. Despite the harsh environmental conditions they endure, these extremophiles represent a great deal of the Earth's biodiversity. Here, for the first time in a member of the archaeal domain, we report a genome-wide assay of spontaneous mutations in the halophilic species Haloferax volcanii using a direct and unbiased method: mutation accumulation experiments combined with deep whole-genome sequencing. H. volcanii is a key model organism not only for the study of halophilicity, but also for archaeal biology in general. Our methods measure the genome-wide rate, spectrum, and spatial distribution of spontaneous mutations. The estimated base substitution rate of 3.15 × 10^-10 per site per generation, or 0.0012 per genome per generation, is similar to the value found in mesophilic prokaryotes (optimal growth at ∼20-45°). This study contributes to a comprehensive phylogenetic view of how evolutionary forces and molecular mechanisms shape the rate and molecular spectrum of mutations across the tree of life.

Development of a supF-based mutation-detection system in the extreme thermophile Thermus thermophilus HB27

Thermus thermophilus (T. thermophilus) HB27 is an extreme thermophile that grows optimally at 65-72 °C. Heat-induced DNA lesions are expected to occur at a higher frequency in the genome of T. thermophilus than in those of mesophiles; however, the mechanisms underlying the maintenance of genome integrity at high temperatures remain poorly understood. The study of mutation spectra has become a powerful approach to understanding the molecular mechanisms responsible for DNA repair and mutagenesis in mesophilic species. Therefore, we developed a supF-based system to detect a broad spectrum of mutations in T. thermophilus. This system was validated by measuring spontaneous mutations in the wild type and a udgA, B double mutant deficient in uracil-DNA glycosylase (UDG) activity. We found that the mutation frequency of the udgA, B strain was 4.7-fold higher than that of the wild type and G:C→A:T transitions dominated, which was the most reasonable for the mutator phenotype associated with the loss of UDG function in T. thermophilus. These results show that this system allowed for the rapid analysis of mutations in T. thermophilus, and may be useful for studying the molecular mechanisms responsible for DNA repair and mutagenesis in this extreme thermophile.

Asexual Amoebae Escape Muller's Ratchet through Polyploidy

While some amoebae reproduce sexually, many amoebae (e.g., Acanthamoeba, Naegleria) reproduce asexually and therefore, according to popular doctrine, are likely to have been genetically disadvantaged as a consequence. In the absence of sex, mutations are proposed to accumulate by a mechanism known as Muller's ratchet. I hypothesise that amoebae can escape the ravages of accumulated mutation by virtue of their being polyploid. The polyploid state reduces spontaneous mutation accumulation by gene conversion, the freshly mutated copy being corrected by the presence of the many other wild-type copies. In this manner these amoebae reap the benefits of an asexual reproductive existence: principally, that it is rapid and convenient. Evidence for this mechanism comes from polyploid plants, bacteria, and archaea.

The apt/6-Methylpurine Counterselection System and Its Applications in Genetic Studies of the Hyperthermophilic Archaeon Sulfolobus islandicus

Sulfolobus islandicus serves as a model for studying archaeal biology as well as linking novel biology to evolutionary ecology using functional population genomics. In the present study, we developed a new counterselectable genetic marker in S. islandicus to expand the genetic toolbox for this species. We show that resistance to the purine analog 6-methylpurine (6-MP) in S. islandicus M.16.4 is due to the inactivation of a putative adenine phosphoribosyltransferase encoded by M164_0158 (apt). The application of the apt gene as a novel counterselectable marker was first illustrated by constructing an unmarked α-amylase deletion mutant. Furthermore, the 6-MP counterselection feature was employed in a forward (loss-of-function) mutation assay to reveal the profile of spontaneous mutations in S. islandicus M.16.4 at the apt locus. Moreover, the general conservation of apt genes in the crenarchaea suggests that the same strategy can be broadly applied to other crenarchaeal model organisms. These results demonstrate that the apt locus represents a new tool for genetic manipulation and sequence analysis of the hyperthermophilic crenarchaeon S. islandicus

IMPORTANCE

Currently, the pyrEF/5-fluoroorotic acid (5-FOA) counterselection system remains the sole counterselection marker in crenarchaeal genetics. Since most Sulfolobus mutants constructed by the research community were derived from genetic hosts lacking the pyrEF genes, the pyrEF/5-FOA system is no longer available for use in forward mutation assays. Demonstration of the apt/6-MP counterselection system for the Sulfolobus model renders it possible to again study the mutation profiles in mutants that have already been constructed by the use of strains with a pyrEF-deficient background. Furthermore, additional counterselectable markers will allow us to conduct more sophisticated genetic studies, i.e., investigate mechanisms of chromosomal DNA transfer and quantify recombination frequencies among S. islandicus strains.

Adaptive engineering of a hyperthermophilic archaeon on CO and discovering the underlying mechanism by multi-omics analysis

The hyperthermophilic archaeon Thermococcus onnurineus NA1 can grow and produce H₂ on carbon monoxide (CO) and its H₂ production rates have been improved through metabolic engineering. In this study, we applied adaptive evolution to enhance H₂ productivity. After over 150 serial transfers onto CO medium, cell density, CO consumption rate and H₂ production rate increased. The underlying mechanism for those physiological changes could be explained by using multi-omics approaches including genomic, transcriptomic and epigenomic analyses. A putative transcriptional regulator was newly identified to regulate the expression levels of genes related to CO oxidation. Transcriptome analysis revealed significant changes in the transcript levels of genes belonging to the categories of transcription, translation and energy metabolism. Our study presents the first genome-scale methylation pattern of hyperthermophilic archaea. Adaptive evolution led to highly enhanced H₂ productivity at high CO flow rates using synthesis gas produced from coal gasification.

A mutant ('lab strain') of the hyperthermophilic archaeon Pyrococcus furiosus, lacking flagella, has unusual growth physiology

A mutant ('lab strain') of the hyperthermophilic archaeon Pyrococcus furiosus DSM3638 exhibited an extended exponential phase and atypical cell aggregation behavior. Genomic DNA from the mutant culture was sequenced and compared to wild-type (WT) DSM3638, revealing 145 genes with one or more insertions, deletions, or substitutions (12 silent, 33 amino acid substitutions, and 100 frame shifts). Approximately, half of the mutated genes were transposases or hypothetical proteins. The WT transcriptome revealed numerous changes in amino acid and pyrimidine biosynthesis pathways coincidental with growth phase transitions, unlike the mutant whose transcriptome reflected the observed prolonged exponential phase. Targeted gene deletions, based on frame-shifted ORFs in the mutant genome, in a genetically tractable strain of P. furiosus (COM1) could not generate the extended exponential phase behavior observed for the mutant. For example, a putative radical SAM family protein (PF2064) was the most highly up-regulated ORF (>25-fold) in the WT between exponential and stationary phase, although this ORF was unresponsive in the mutant; deletion of this gene in P. furiosus COM1 resulted in no apparent phenotype. On the other hand, frame-shifting mutations in the mutant genome negatively impacted transcription of a flagellar biosynthesis operon (PF0329-PF0338).Consequently, cells in the mutant culture lacked flagella and, unlike the WT, showed minimal evidence of exopolysaccharide-based cell aggregation in post-exponential phase. Electron microscopy of PF0331-PF0337 deletions in P. furiosus COM1 showed that absence of flagella impacted normal cell aggregation behavior and, furthermore, indicated that flagella play a key role, beyond motility, in the growth physiology of P. furiosus.

Growth temperature and genome size in bacteria are negatively correlated, suggesting genomic streamlining during thermal adaptation

Prokaryotic genomes are small and compact. Either this feature is caused by neutral evolution or by natural selection favoring small genomes—genome streamlining. Three separate prior lines of evidence argue against streamlining for most prokaryotes. We find that the same three lines of evidence argue for streamlining in the genomes of thermophile bacteria. Specifically, with increasing habitat temperature and decreasing genome size, the proportion of genomic DNA in intergenic regions decreases. Furthermore, with increasing habitat temperature, generation time decreases. Genome-wide selective constraints do not decrease as in the reduced genomes of host-associated species. Reduced habitat variability is not a likely explanation for the smaller genomes of thermophiles. Genome size may be an indirect target of selection due to its association with cell volume. We use metabolic modeling to demonstrate that known changes in cell structure and physiology at high temperature can provide a selective advantage to reduce cell volume at high temperatures.

Confounders of mutation-rate estimators: selection and phenotypic lag in Thermus thermophilus

In a recent description of the rate and character of spontaneous mutation in the hyperthermophilic bacterium Thermus thermophilus, the mutation rate was observed to be substantially lower than seen in several mesophiles. Subsequently, a report appeared indicating that this bacterium maintains an average of about 4.5 genomes per cell. This number of genomes might result in a segregation lag for the expression of a recessive mutation and might therefore lead to an underestimate of the rate of mutation. Here we describe some kinds of problems that may arise when estimating mutation rates and outline ways to adjust the rates accordingly. The emphasis is mainly on differential rates of growth of mutants versus their parents and on various kinds of phenotypic lag. We then apply these methods to the T. thermophilus data and conclude that there is as yet no reliable impact on a previously described rate.

Viral diversity threshold for adaptive immunity in prokaryotes

Bacteria and archaea face continual onslaughts of rapidly diversifying viruses and plasmids. Many prokaryotes maintain adaptive immune systems known as clustered regularly interspaced short palindromic repeats (CRISPR) and CRISPR-associated genes (Cas). CRISPR-Cas systems are genomic sensors that serially acquire viral and plasmid DNA fragments (spacers) that are utilized to target and cleave matching viral and plasmid DNA in subsequent genomic invasions, offering critical immunological memory. Only 50% of sequenced bacteria possess CRISPR-Cas immunity, in contrast to over 90% of sequenced archaea. To probe why half of bacteria lack CRISPR-Cas immunity, we combined comparative genomics and mathematical modeling. Analysis of hundreds of diverse prokaryotic genomes shows that CRISPR-Cas systems are substantially more prevalent in thermophiles than in mesophiles. With sequenced bacteria disproportionately mesophilic and sequenced archaea mostly thermophilic, the presence of CRISPR-Cas appears to depend more on environmental temperature than on bacterial-archaeal taxonomy. Mutation rates are typically severalfold higher in mesophilic prokaryotes than in thermophilic prokaryotes. To quantitatively test whether accelerated viral mutation leads microbes to lose CRISPR-Cas systems, we developed a stochastic model of virus-CRISPR coevolution. The model competes CRISPR-Cas-positive (CRISPR-Cas+) prokaryotes against CRISPR-Cas-negative (CRISPR-Cas−) prokaryotes, continually weighing the antiviral benefits conferred by CRISPR-Cas immunity against its fitness costs. Tracking this cost-benefit analysis across parameter space reveals viral mutation rate thresholds beyond which CRISPR-Cas cannot provide sufficient immunity and is purged from host populations. These results offer a simple, testable viral diversity hypothesis to explain why mesophilic bacteria disproportionately lack CRISPR-Cas immunity. More generally, fundamental limits on the adaptability of biological sensors (Lamarckian evolution) are predicted.

A remarkable recent discovery in microbiology is that bacteria and archaea possess systems conferring immunological memory and adaptive immunity. Clustered regularly interspaced short palindromic repeats (CRISPR) and CRISPR-associated genes (CRISPR-Cas) are genomic sensors that allow prokaryotes to acquire DNA fragments from invading viruses and plasmids. Providing immunological memory, these stored fragments destroy matching DNA in future viral and plasmid invasions. CRISPR-Cas systems also provide adaptive immunity, keeping up with mutating viruses and plasmids by continually acquiring new DNA fragments. Surprisingly, less than 50% of mesophilic bacteria, in contrast to almost 90% of thermophilic bacteria and Archaea, maintain CRISPR-Cas immunity. Using mathematical modeling, we probe this dichotomy, showing how increased viral mutation rates can explain the reduced prevalence of CRISPR-Cas systems in mesophiles. Rapidly mutating viruses outrun CRISPR-Cas immune systems, likely decreasing their prevalence in bacterial populations. Thus, viral adaptability may select against, rather than for, immune adaptability in prokaryotes.

Proteome analyses of hydrogen-producing hyperthermophilic archaeon Thermococcus onnurineus NA1 in different one-carbon substrate culture conditions

Thermococcus onnurineus NA1, a sulfur-reducing hyperthermophilic archaeon, is capable of H(2)-producing growth, considered to be hydrogenogenic carboxydotrophy. Utilization of formate as a sole energy source has been well studied in T. onnurineus NA1. However, whether formate can be used as its carbon source remains unknown. To obtain a global view of the metabolic characteristics of H(2)-producing growth, a quantitative proteome analysis of T. onnurineus NA1 grown on formate, CO, and starch was performed by combining one-dimensional SDS-PAGE with nano UPLC-MS(E). A total of 587 proteins corresponding to 29.7% of the encoding genes were identified, and the major metabolic pathways (especially energy metabolism) were characterized at the protein level. Expression of glycolytic enzymes was common but more highly induced in starch-grown cells. In contrast, enzymes involved in key steps of the gluconeogenesis and pentose phosphate pathways were strongly up-regulated in formate-grown cells, suggesting that formate could be utilized as a carbon source by T. onnurineus NA1. In accordance with the genomic analysis, comprehensive proteomic analysis also revealed a number of hydrogenase clusters apparently associated with formate metabolism. On the other hand, CODH and CO-induced hydrogenases belonging to the Hyg4-II cluster, as well as sulfhydrogenase-I and Mbx, were prominently expressed during CO culture. Our data suggest that CO can be utilized as a sole energy source for H(2) production via an electron transport mechanism and that CO(2) produced from catabolism or CO oxidation by CODH and CO-induced hydrogenases may subsequently be assimilated into the organic carbon. Overall, proteomic comparison of formate- and CO-grown cells with starch-grown cells revealed that a single carbon compound, such as formate and CO, can be utilized as an efficient substrate to provide cellular carbon and/or energy by T. onnurineus NA1.

Measurements of spontaneous rates of mutations in the recent past and the near future

The rate of spontaneous mutation in natural populations is a fundamental parameter for many evolutionary phenomena. Because the rate of mutation is generally low, most of what is currently known about mutation has been obtained through indirect, complex and imprecise methodological approaches. However, in the past few years genome-wide sequencing of closely related individuals has made it possible to estimate the rates of mutation directly at the level of the DNA, avoiding most of the problems associated with using indirect methods. Here, we review the methods used in the past with an emphasis on next generation sequencing, which may soon make the accurate measurement of spontaneous mutation rates a matter of routine.

MutS and MutL are dispensable for maintenance of the genomic mutation rate in the halophilic archaeon Halobacterium salinarum NRC-1

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.

Avoiding dangerous missense: thermophiles display especially low mutation rates

Rates of spontaneous mutation have been estimated under optimal growth conditions for a variety of DNA-based microbes, including viruses, bacteria, and eukaryotes. When expressed as genomic mutation rates, most of the values were in the vicinity of 0.003–0.004 with a range of less than two-fold. Because the genome sizes varied by roughly 10⁴-fold, the mutation rates per average base pair varied inversely by a similar factor. Even though the commonality of the observed genomic rates remains unexplained, it implies that mutation rates in unstressed microbes reach values that can be finely tuned by evolution. An insight originating in the 1920s and maturing in the 1960s proposed that the genomic mutation rate would reflect a balance between the deleterious effect of the average mutation and the cost of further reducing the mutation rate. If this view is correct, then increasing the deleterious impact of the average mutation should be countered by reducing the genomic mutation rate. It is a common observation that many neutral or nearly neutral mutations become strongly deleterious at higher temperatures, in which case they are called temperature-sensitive mutations. Recently, the kinds and rates of spontaneous mutations were described for two microbial thermophiles, a bacterium and an archaeon. Using an updated method to extrapolate from mutation-reporter genes to whole genomes reveals that the rate of base substitutions is substantially lower in these two thermophiles than in mesophiles. This result provides the first experimental support for the concept of an evolved balance between the total genomic impact of mutations and the cost of further reducing the basal mutation rate.

Spontaneous mutations are key drivers of evolution and disease. In microbes, most mutations are deleterious, some are neutral (without significant impact), and a few are advantageous. Because deleterious mutations reduce fitness, there should be constant selection for antimutator mutations that reduce rates of spontaneous mutation. However, such reductions are necessarily achieved at some cost. Therefore, a mutation rate should converge evolutionarily on a value that reflects this trade-off. For DNA microbes, the observed genomic mutation rate is remarkably (and mysteriously) invariant, in the neighborhood of 0.003–0.004, with a range of less than two-fold despite huge variation per average base pair in organisms with a wide diversity of life histories. Would an environmental condition that increased the average deleterious impact of a mutation be balanced by additional investments in antimutator mutations? It is widely observed that many mutations with mild impacts become strongly deleterious at higher temperatures, so mutation rates were measured in two thermophiles, a bacterium and an archaeon. Remarkably, both displayed average mutation rates reduced by about five-fold from the characteristic mesophilic value, most of the decrease reflecting a 10-fold reduction in the rate of base substitutions.