The Presence of the DNA Repair Genes mutM, mutY, mutL, and mutS is Related to Proteome Size in Bacterial Genomes

Metagenome mining and functional analysis reveal oxidized guanine DNA repair at the Lost City Hydrothermal Field

The GO DNA repair system protects against GC → TA mutations by finding and removing oxidized guanine. The system is mechanistically well understood but its origins are unknown. We searched metagenomes and abundantly found the genes encoding GO DNA repair at the Lost City Hydrothermal Field (LCHF). We recombinantly expressed the final enzyme in the system to show MutY homologs function to suppress mutations. Microbes at the LCHF thrive without sunlight, fueled by the products of geochemical transformations of seafloor rocks, under conditions believed to resemble a young Earth. High levels of the reductant H2 and low levels of O2 in this environment raise the question, why are resident microbes equipped to repair damage caused by oxidative stress? MutY genes could be assigned to metagenome-assembled genomes (MAGs), and thereby associate GO DNA repair with metabolic pathways that generate reactive oxygen, nitrogen and sulfur species. Our results indicate that cell-based life was under evolutionary pressure to cope with oxidized guanine well before O2 levels rose following the great oxidation event.

Host specific adaptations of Ligilactobacillus aviarius to poultry

The genus Ligilactobacillus encompasses species adapted to vertebrate hosts and fermented food. Their genomes encode adaptations to the host lifestyle. Reports of gut microbiota from chicken and turkey gastrointestinal tract have shown a high persistence of Ligilactobacillus aviarius along the digestive system compared to other species found in the same host. However, its adaptations to poultry as a host has not yet been described. In this work, the pan-genome of Ligilactobacillus aviarius was explored to describe the functional adaptability to the gastrointestinal environment. The core genome is composed of 1179 gene clusters that are present at least in one copy that codifies to structural, ribosomal and biogenesis proteins. The rest of the identified regions were classified into three different functional clusters of orthologous groups (clusters) that codify carbohydrate metabolism, envelope biogenesis, viral defence mechanisms, and mobilome inclusions. The pan-genome of Ligilactobacillus aviarius is a closed pan-genome, frequently found in poultry and highly prevalent across chicken faecal samples. The genome of L. aviarius codifies different clusters of glycoside hydrolases and glycosyltransferases that mediate interactions with the host cells. Accessory features, such as antiviral mechanisms and prophage inclusions, variate amongst strains from different GIT sections. This information provides hints about the interaction of this species with viral particles and other bacterial species. This work highlights functional adaptability traits present in L. aviarius that make it a dominant key member of the poultry gut microbiota and enlightens the convergent ecological relation of this species to the poultry gut environment.

Epigenetic status of FBXW7 gene and its role in Ovarian cancer pathogenesis

BACKGROUND

Chromatin immunoprecipitation (ChIP) analysis revealed that the FBXW7 gene and the long non-coding RNA (LINC01588) are potential candidates in epithelial ovarian cancer (EOC) pathogenesis. However, their exact role in EOC is not yet known. Thus, the present study sheds light on the impact of the mutations/ methylation status of the FBXW7 gene.

MATERIALS AND METHODS

We used public databases to assess the correlation between mutations/ methylation status and the FBXW7 expression. Furthermore, we performed Pearson's correlation analysis between the FBXW7 gene and LINC01588. We performed gene panel exome sequencing and Methylation-specific PCR (MSP) in HOSE 6-3, MCAS, OVSAHO, and eight EOC patients' samples to validate the bioinformatics results.

RESULTS

The FBXW7 gene was less expressed in EOC, particularly in stages III and IV, compared to healthy tissues. Furthermore, bioinformatics analysis, gene panel exome sequencing, and MSP revealed that the FBXW7 gene is neither mutated nor methylated in EOC cell lines and tissues, suggesting alternative mechanisms for FBXW7 gene regulation. Interestingly, Pearson's correlation analysis showed an inverse, significant correlation between the FBXW7 gene and LINC01588  expression, suggesting a potential regulatory role of LINC01588.

CONCLUSION

Neither mutations nor methylation is the causative mechanism for the FBXW7 downregulation in EOC, suggesting alternative means involving the lncRNA LINC01588.

Evidence of gene nucleotide composition favoring replication and growth in a fastidious plant pathogen

Nucleotide composition (GC content) varies across bacteria species, genome regions, and specific genes. In Xylella fastidiosa, a vector-borne fastidious plant pathogen infecting multiple crops, GC content ranges between ∼51-52%; however, these values were gathered using limited genomic data. We evaluated GC content variations across X. fastidiosa subspecies fastidiosa (N = 194), subsp. pauca (N = 107), and subsp. multiplex (N = 39). Genomes were classified based on plant host and geographic origin; individual genes within each genome were classified based on gene function, strand, length, ortholog group, Core vs. Accessory, and Recombinant vs. Non-recombinant. GC content was calculated for each gene within each evaluated genome. The effects of genome and gene level variables were evaluated with a mixed effect ANOVA, and the marginal-GC content was calculated for each gene. Also, the correlation between gene-specific GC content vs. natural selection (dN/dS) and recombination/mutation (r/m) was estimated. Our analyses show that intra-genomic changes in nucleotide composition in X. fastidiosa are small and influenced by multiple variables. Higher AT-richness is observed in genes involved in replication and translation, and genes in the leading strand. In addition, we observed a negative correlation between high-AT and dN/dS in subsp. pauca. The relationship between recombination and GC content varied between core and accessory genes. We hypothesize that distinct evolutionary forces and energetic constraints both drive and limit these small variations in nucleotide composition.

Proteome size reduction in Apicomplexans is linked with loss of DNA repair and host redundant pathways

Apicomplexans are alveolate parasites which include Plasmodium falciparum, the main cause of malaria, one of the world's biggest killers from infectious disease. Apicomplexans are characterized by a reduction in proteome size, which appears to result from metabolic and functional simplification, commensurate with their parasitic lifestyle. However, other factors may also help to explain gene loss such as population bottlenecks experienced during transmission, and the effect of reducing the overall genomic information content. The latter constitutes an 'informational constraint', which is proposed to exert a selective pressure to evolve and maintain genes involved in informational fidelity and error correction, proportional to the quantity of information in the genome (which approximates to proteome size). The dynamics of gene loss was examined in 41 Apicomplexan genomes using orthogroup analysis. We show that loss of genes involved in amino acid metabolism and steroid biosynthesis can be explained by metabolic redundancy with the host. We also show that there is a marked tendency to lose DNA repair genes as proteome size is reduced. This may be explained by a reduction in size of the informational constraint and can help to explain elevated mutation rates in pathogens with reduced genome size. Multiple Sequentially Markovian Coalescent (MSMC) analysis indicates a recent bottleneck, consistent with predictions generated using allele-based population genetics approaches, implying that relaxed selection pressure due to reduced population size might have contributed to gene loss. However, the non-randomness of pathways that are lost challenges this scenario. Lastly, we identify unique orthogroups in malaria-causing Plasmodium species that infect humans, with a high proportion of membrane associated proteins. Thus, orthogroup analysis appears useful for identifying novel candidate pathogenic factors in parasites, when there is a wide sample of genomes available.

Synergism between the Black Queen effect and the proteomic constraint on genome size reduction in the photosynthetic picoeukaryotes

The photosynthetic picoeukaryotes (PPEs) comprise a rare example of free-living eukaryotes that have undergone genome reduction. Here, we examine a duality in the process; the proposed driver of genome reduction (the Black Queen hypothesis, BQH), and the resultant impact of genome information loss (the Proteomic Constraint hypothesis, PCH). The BQH predicts that some metabolites may be shared in the open ocean, thus driving loss of redundant metabolic pathways in individual genomes. In contrast, the PCH predicts that as the information content of a genome is reduced, the total mutation load is also reduced, leading to loss of DNA repair genes due to the resulting reduction in selective constraint. Consistent with the BQH, we observe that biosynthetic pathways involved with soluble metabolites such as amino acids and carotenoids are preferentially lost from the PPEs, in contrast to biosynthetic pathways involved with insoluble metabolites, such as lipids, which are retained. Consistent with the PCH, a correlation between proteome size and the number of DNA repair genes, and numerous other informational categories, is observed. While elevated mutation rates resulting from the loss of DNA repair genes have been linked to reduced effective population sizes in intracellular bacteria, this remains to be established. This study shows that in microbial species with large population sizes, an underlying factor in modulating their DNA repair capacity appears to be information content.

Temperature responses of mutation rate and mutational spectrum in an Escherichia coli strain and the correlation with metabolic rate

Background

Temperature is a major determinant of spontaneous mutation, but the precise mode, and the underlying mechanisms, of the temperature influences remain less clear. Here we used a mutation accumulation approach combined with whole-genome sequencing to investigate the temperature dependence of spontaneous mutation in an Escherichia coli strain. Experiments were performed under aerobic conditions at 25, 28 and 37 °C, three temperatures that were non-stressful for the bacterium but caused significantly different bacterial growth rates.

Results

Mutation rate did not differ between 25 and 28 °C, but was higher at 37 °C. Detailed analyses of the molecular spectrum of mutations were performed; and a particularly interesting finding is that higher temperature led to a bias of mutation to coding, relative to noncoding, DNA. Furthermore, the temperature response of mutation rate was extremely similar to that of metabolic rate, consistent with an idea that metabolic rate predicts mutation rate.

Conclusions

Temperature affects mutation rate and the types of mutation supply, both being crucial for the opportunity of natural selection. Our results help understand how temperature drives evolutionary speed of organisms and thus the global patterns of biodiversity. This study also lend support to the metabolic theory of ecology for linking metabolic rate and molecular evolution rate.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1252-8) contains supplementary material, which is available to authorized users.

Plasmodium parasites of birds have the most AT-rich genes of eukaryotes

The genomic architecture of organisms, including nucleotide composition, can be highly variable, even among closely-related species. To better understand the causes leading to structural variation in genomes, information on distinct and diverse genomic features is needed. Malaria parasites are known for encompassing a wide range of genomic GC-content and it has long been thought that Plasmodium falciparum, the virulent malaria parasite of humans, has the most AT-biased eukaryotic genome. Here, I perform comparative genomic analyses of the most AT-rich eukaryotes sequenced to date, and show that the avian malaria parasites Plasmodium gallinaceum, P. ashfordi, and P. relictum have the most extreme coding sequences in terms of AT-bias. Their mean GC-content is 21.21, 21.22 and 21.60 %, respectively, which is considerably lower than the transcriptome of P. falciparum (23.79 %) and other eukaryotes. This information enables a better understanding of genome evolution and raises the question of how certain organisms are able to prosper despite severe compositional constraints.

MicroRNAs as biomarkers associated with bladder cancer

Bladder cancer is the fifth most common cancer with significant morbidity and mortality. Recently, numerous studies demonstrated that microRNAs are emerging as diagnostic biomarkers for bladder cancer. Specific miRNA profiles have been identified for several samples from patients with bladder cancer. MicroRNAs are noncoding RNA molecules of approximately 23 nucleotides that play important roles in multiple steps during the progression of bladder cancer. Here, we review the expression profiles of miRNAs and their biological functions, regulation, and clinical implications in bladder cancer. Either downregulation or upregulation of miRNAs occurs in bladder cancer through epigenetic changes or defects of the biogenesis apparatus. Deregulation of miRNAs is involved in cell cycle arrest, apoptosis, proliferation, metastasis, drug resistance, and other functions in bladder cancer. A number of miRNAs, have been associated with tumor type, stage, or patient survival, and miRNAs might be developed as diagnostic or prognostic markers. A better understanding of the roles of miRNAs in bladder cancer will shed light on the molecular mechanisms of bladder cancer.

The Rate and Spectrum of Spontaneous Mutations in Mycobacterium smegmatis, a Bacterium Naturally Devoid of the Postreplicative Mismatch Repair Pathway

Mycobacterium smegmatis is a bacterium that is naturally devoid of known postreplicative DNA mismatch repair (MMR) homologs, mutS and mutL, providing an opportunity to investigate how the mutation rate and spectrum has evolved in the absence of a highly conserved primary repair pathway. Mutation accumulation experiments of M. smegmatis yielded a base-substitution mutation rate of 5.27 × 10⁻¹⁰ per site per generation, or 0.0036 per genome per generation, which is surprisingly similar to the mutation rate in MMR-functional unicellular organisms. Transitions were found more frequently than transversions, with the A:T→G:C transition rate significantly higher than the G:C→A:T transition rate, opposite to what is observed in most studied bacteria. We also found that the transition-mutation rate of M. smegmatis is significantly lower than that of other naturally MMR-devoid or MMR-knockout organisms. Two possible candidates that could be responsible for maintaining high DNA fidelity in this MMR-deficient organism are the ancestral-like DNA polymerase DnaE1, which contains a highly efficient DNA proofreading histidinol phosphatase (PHP) domain, and/or the existence of a uracil-DNA glycosylase B (UdgB) homolog that might protect the GC-rich M. smegmatis genome against DNA damage arising from oxidation or deamination. Our results suggest that M. smegmatis has a noncanonical Dam (DNA adenine methylase) methylation system, with target motifs differing from those previously reported. The mutation features of M. smegmatis provide further evidence that genomes harbor alternative routes for improving replication fidelity, even in the absence of major repair pathways.

The properties of spontaneous mutations in the opportunistic pathogen Pseudomonas aeruginosa

Background

Natural genetic variation ultimately arises from the process of mutation. Knowledge of how the raw material for evolution is produced is necessary for a full understanding of several fundamental evolutionary concepts. We performed a mutation accumulation experiment with wild-type and mismatch-repair deficient, mutator lines of the pathogenic bacterium Pseudomonas aeruginosa, and used whole-genome sequencing to reveal the genome-wide rate, spectrum, distribution, leading/lagging bias, and context-dependency of spontaneous mutations.

Results

Wild-type base-pair mutation and indel rates were ~10⁻¹⁰ and ~10⁻¹¹ per nucleotide per generation, respectively, and deficiencies in the mismatch-repair system caused rates to increase by over two orders of magnitude. A universal bias towards AT was observed in wild-type lines, but was reversed in mutator lines to a bias towards GC. Biases for which types of mutations occurred during replication of the leading versus lagging strand were detected reciprocally in both replichores. The distribution of mutations along the chromosome was non-random, with peaks near the terminus of replication and at positions intermediate to the replication origin and terminus. A similar distribution bias was observed along the chromosome in natural populations of P. aeruginosa. Site-specific mutation rates were higher when the focal nucleotide was immediately flanked by C:G pairings.

Conclusions

Whole-genome sequencing of mutation accumulation lines allowed the comprehensive identification of mutations and revealed what factors of molecular and genomic architecture affect the mutational process. Our study provides a more complete view of how several mechanisms of mutation, mutation repair, and bias act simultaneously to produce the raw material for evolution.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2244-3) contains supplementary material, which is available to authorized users.

Mutation--The Engine of Evolution: Studying Mutation and Its Role in the Evolution of Bacteria

Mutation is the engine of evolution in that it generates the genetic variation on which the evolutionary process depends. To understand the evolutionary process we must therefore characterize the rates and patterns of mutation. Starting with the seminal Luria and Delbruck fluctuation experiments in 1943, studies utilizing a variety of approaches have revealed much about mutation rates and patterns and about how these may vary between different bacterial strains and species along the chromosome and between different growth conditions. This work provides a critical overview of the results and conclusions drawn from these studies, of the debate surrounding some of these conclusions, and of the challenges faced when studying mutation and its role in bacterial evolution.

DNA Repair Is Associated with Information Content in Bacteria, Archaea, and DNA Viruses

The concept of a "proteomic constraint" proposes that DNA repair capacity is positively correlated with the information content of a genome, which can be approximated to the size of the proteome (P). This in turn implies that DNA repair genes are more likely to be present in genomes with larger values of P. This stands in contrast to the common assumption that informational genes have a core function and so are evenly distributed across organisms. We examined the presence/absence of 18 DNA repair genes in bacterial genomes. A positive relationship between gene presence and P was observed for 17 genes in the total dataset, and 16 genes when only nonintracellular bacteria were examined. A marked reduction of DNA repair genes was observed in intracellular bacteria, consistent with their reduced value of P. We also examined archaeal and DNA virus genomes, and show that the presence of DNA repair genes is likewise related to a larger value of P. In addition, the products of the bacterial genes mutY, vsr, and ndk, involved in the correction of GC/AT mutations, are strongly associated with reduced genome GC content. We therefore propose that a reduction in information content leads to a loss of DNA repair genes and indirectly to a reduction in genome GC content in bacteria by exposure to the underlying AT mutation bias. The reduction in P may also indirectly lead to the increase in substitution rates observed in intracellular bacteria via loss of DNA repair genes.

Genetic code evolution reveals the neutral emergence of mutational robustness, and information as an evolutionary constraint

The standard genetic code (SGC) is central to molecular biology and its origin and evolution is a fundamental problem in evolutionary biology, the elucidation of which promises to reveal much about the origins of life. In addition, we propose that study of its origin can also reveal some fundamental and generalizable insights into mechanisms of molecular evolution, utilizing concepts from complexity theory. The first is that beneficial traits may arise by non-adaptive processes, via a process of "neutral emergence". The structure of the SGC is optimized for the property of error minimization, which reduces the deleterious impact of point mutations. Via simulation, it can be shown that genetic codes with error minimization superior to the SGC can emerge in a neutral fashion simply by a process of genetic code expansion via tRNA and aminoacyl-tRNA synthetase duplication, whereby similar amino acids are added to codons related to that of the parent amino acid. This process of neutral emergence has implications beyond that of the genetic code, as it suggests that not all beneficial traits have arisen by the direct action of natural selection; we term these "pseudaptations", and discuss a range of potential examples. Secondly, consideration of genetic code deviations (codon reassignments) reveals that these are mostly associated with a reduction in proteome size. This code malleability implies the existence of a proteomic constraint on the genetic code, proportional to the size of the proteome (P), and that its reduction in size leads to an "unfreezing" of the codon - amino acid mapping that defines the genetic code, consistent with Crick's Frozen Accident theory. The concept of a proteomic constraint may be extended to propose a general informational constraint on genetic fidelity, which may be used to explain variously, differences in mutation rates in genomes with differing proteome sizes, differences in DNA repair capacity and genome GC content between organisms, a selective pressure in the evolution of sexual reproduction, and differences in translational fidelity. Lastly, the utility of the concept of an informational constraint to other diverse fields of research is explored.

Genome evolution in an ancient bacteria-ant symbiosis: parallel gene loss among Blochmannia spanning the origin of the ant tribe Camponotini

Stable associations between bacterial endosymbionts and insect hosts provide opportunities to explore genome evolution in the context of established mutualisms and assess the roles of selection and genetic drift across host lineages and habitats. Blochmannia, obligate endosymbionts of ants of the tribe Camponotini, have coevolved with their ant hosts for ∼40 MY. To investigate early events in Blochmannia genome evolution across this ant host tribe, we sequenced Blochmannia from two divergent host lineages, Colobopsis obliquus and Polyrhachis turneri, and compared them with four published genomes from Blochmannia of Camponotus sensu stricto. Reconstructed gene content of the last common ancestor (LCA) of these six Blochmannia genomes is reduced (690 protein coding genes), consistent with rapid gene loss soon after establishment of the symbiosis. Differential gene loss among Blochmannia lineages has affected cellular functions and metabolic pathways, including DNA replication and repair, vitamin biosynthesis and membrane proteins. Blochmannia of P. turneri (i.e., B. turneri) encodes an intact DnaA chromosomal replication initiation protein, demonstrating that loss of dnaA was not essential for establishment of the symbiosis. Based on gene content, B. obliquus and B. turneri are unable to provision hosts with riboflavin. Of the six sequenced Blochmannia, B. obliquus is the earliest diverging lineage (i.e., the sister group of other Blochmannia sampled) and encodes the fewest protein-coding genes and the most pseudogenes. We identified 55 genes involved in parallel gene loss, including glutamine synthetase, which may participate in nitrogen recycling. Pathways for biosynthesis of coenzyme A, terpenoids and riboflavin were lost in multiple lineages, suggesting relaxed selection on the pathway after inactivation of one component. Analysis of Illumina read datasets did not detect evidence of plasmids encoding missing functions, nor the presence of coresident symbionts other than Wolbachia. Although gene order is strictly conserved in four Blochmannia of Camponotus sensu stricto, comparisons with deeply divergent lineages revealed inversions in eight genomic regions, indicating ongoing recombination despite ancestral loss of recA. In sum, the addition of two Blochmannia genomes of divergent host lineages enables reconstruction of early events in evolution of this symbiosis and suggests that Blochmannia lineages may experience distinct, host-associated selective pressures. Understanding how evolutionary forces shape genome reduction in this system may help to clarify forces driving gene loss in other bacteria, including intracellular pathogens.

An evolutionary analysis of genome expansion and pathogenicity in Escherichia coli

BACKGROUND

There are several studies describing loss of genes through reductive evolution in microbes, but how selective forces are associated with genome expansion due to horizontal gene transfer (HGT) has not received similar attention. The aim of this study was therefore to examine how selective pressures influence genome expansion in 53 fully sequenced and assembled Escherichia coli strains. We also explored potential connections between genome expansion and the attainment of virulence factors. This was performed using estimations of several genomic parameters such as AT content, genomic drift (measured using relative entropy), genome size and estimated HGT size, which were subsequently compared to analogous parameters computed from the core genome consisting of 1729 genes common to the 53 E. coli strains. Moreover, we analyzed how selective pressures (quantified using relative entropy and dN/dS), acting on the E. coli core genome, influenced lineage and phylogroup formation.

RESULTS

Hierarchical clustering of dS and dN estimations from the E. coli core genome resulted in phylogenetic trees with topologies in agreement with known E. coli taxonomy and phylogroups. High values of dS, compared to dN, indicate that the E. coli core genome has been subjected to substantial purifying selection over time; significantly more than the non-core part of the genome (p<0.001). This is further supported by a linear association between strain-wise dS and dN values (β = 26.94 ± 0.44, R2~0.98, p<0.001). The non-core part of the genome was also significantly more AT-rich (p<0.001) than the core genome and E. coli genome size correlated with estimated HGT size (p<0.001). In addition, genome size (p<0.001), AT content (p<0.001) as well as estimated HGT size (p<0.005) were all associated with the presence of virulence factors, suggesting that pathogenicity traits in E. coli are largely attained through HGT. No associations were found between selective pressures operating on the E. coli core genome, as estimated using relative entropy, and genome size (p~0.98).

CONCLUSIONS

On a larger time frame, genome expansion in E. coli, which is significantly associated with the acquisition of virulence factors, appears to be independent of selective forces operating on the core genome.

MutY-glycosylase: an overview on mutagenesis and activities beyond the GO system

MutY is a glycosylase known for its role in DNA base excision repair (BER). It is critically important in the prevention of DNA mutations derived from 7,8-dihydro-8-oxoguanine (8-oxoG), which are the major lesions resulting from guanine oxidation. MutY has been described as a DNA repair enzyme in the GO system responsible for removing adenine residues misincorporated in 8-oxoG:A mispairs, avoiding G:C to T:A mutations. Further studies have shown that this enzyme binds to other mispairs, interacts with several enzymes, avoids different transversions/transitions in DNA, and is involved in different repair pathways. Additional activities have been reported for MutY, such as the repair of replication errors in newly synthesized DNA strands through its glycosylase activity. Moreover, MutY is a highly conserved enzyme present in several prokaryotic and eukaryotic organisms. MutY defects are associated with a hereditary colorectal cancer syndrome termed MUTYH-associated polyposis (MAP). Here, we have reviewed the roles of MutY in the repair of mispaired bases in DNA as well as its activities beyond the GO system.

Coregulated genes link sulfide:quinone oxidoreductase and arsenic metabolism in Synechocystis sp. strain PCC6803

Although the biogeochemistry of the two environmentally hazardous compounds arsenic and sulfide has been extensively investigated, the biological interference of these two toxic but potentially energy-rich compounds has only been hypothesized and indirectly proven. Here we provide direct evidence for the first time that in the photosynthetic model organism Synechocystis sp. strain PCC6803 the two metabolic pathways are linked by coregulated genes that are involved in arsenic transport, sulfide oxidation, and probably in sulfide-based alternative photosynthesis. Although Synechocystis sp. strain PCC6803 is an obligate photoautotrophic cyanobacterium that grows via oxygenic photosynthesis, we discovered that specific genes are activated in the presence of sulfide or arsenite to exploit the energy potentials of these chemicals. These genes form an operon that we termed suoRSCT, located on a transposable element of type IS4 on the plasmid pSYSM of the cyanobacterium. suoS (sll5036) encodes a light-dependent, type I sulfide:quinone oxidoreductase. The suoR (sll5035) gene downstream of suoS encodes a regulatory protein that belongs to the ArsR-type repressors that are normally involved in arsenic resistance. We found that this repressor has dual specificity, resulting in 200-fold induction of the operon upon either arsenite or sulfide exposure. The suoT gene encodes a transmembrane protein similar to chromate transporters but in fact functioning as an arsenite importer at permissive concentrations. We propose that the proteins encoded by the suoRSCT operon might have played an important role under anaerobic, reducing conditions on primordial Earth and that the operon was acquired by the cyanobacterium via horizontal gene transfer.

The quest for a unified view of bacterial land colonization

Exploring molecular mechanisms underlying bacterial water-to-land transition represents a critical start toward a better understanding of the functioning and stability of the terrestrial ecosystems. Here, we perform comprehensive analyses based on a large variety of bacteria by integrating taxonomic, phylogenetic and metagenomic data, in the quest for a unified view that elucidates genomic, evolutionary and ecological dynamics of the marine progenitors in adapting to nonaquatic environments. We hypothesize that bacterial land colonization is dominated by a single-gene sweep, that is, the emergence of dnaE2 derived from an early duplication event of the primordial dnaE, followed by a series of niche-specific genomic adaptations, including GC content increase, intensive horizontal gene transfer and constant genome expansion. In addition, early bacterial radiation may be stimulated by an explosion of land-borne hosts (for example, plants and animals) after initial land colonization events.

Amino acid usage is asymmetrically biased in AT- and GC-rich microbial genomes

INTRODUCTION

Genomic base composition ranges from less than 25% AT to more than 85% AT in prokaryotes. Since only a small fraction of prokaryotic genomes is not protein coding even a minor change in genomic base composition will induce profound protein changes. We examined how amino acid and codon frequencies were distributed in over 2000 microbial genomes and how these distributions were affected by base compositional changes. In addition, we wanted to know how genome-wide amino acid usage was biased in the different genomes and how changes to base composition and mutations affected this bias. To carry this out, we used a Generalized Additive Mixed-effects Model (GAMM) to explore non-linear associations and strong data dependences in closely related microbes; principal component analysis (PCA) was used to examine genomic amino acid- and codon frequencies, while the concept of relative entropy was used to analyze genomic mutation rates.

RESULTS

We found that genomic amino acid frequencies carried a stronger phylogenetic signal than codon frequencies, but that this signal was weak compared to that of genomic %AT. Further, in contrast to codon usage bias (CUB), amino acid usage bias (AAUB) was differently distributed in AT- and GC-rich genomes in the sense that AT-rich genomes did not prefer specific amino acids over others to the same extent as GC-rich genomes. AAUB was also associated with relative entropy; genomes with low AAUB contained more random mutations as a consequence of relaxed purifying selection than genomes with higher AAUB.

CONCLUSION

Genomic base composition has a substantial effect on both amino acid- and codon frequencies in bacterial genomes. While phylogeny influenced amino acid usage more in GC-rich genomes, AT-content was driving amino acid usage in AT-rich genomes. We found the GAMM model to be an excellent tool to analyze the genomic data used in this study.

Rate and molecular spectrum of spontaneous mutations in the bacterium Escherichia coli as determined by whole-genome sequencing

Knowledge of the rate and nature of spontaneous mutation is fundamental to understanding evolutionary and molecular processes. In this report, we analyze spontaneous mutations accumulated over thousands of generations by wild-type Escherichia coli and a derivative defective in mismatch repair (MMR), the primary pathway for correcting replication errors. The major conclusions are (i) the mutation rate of a wild-type E. coli strain is ~1 × 10(-3) per genome per generation; (ii) mutations in the wild-type strain have the expected mutational bias for G:C > A:T mutations, but the bias changes to A:T > G:C mutations in the absence of MMR; (iii) during replication, A:T > G:C transitions preferentially occur with A templating the lagging strand and T templating the leading strand, whereas G:C > A:T transitions preferentially occur with C templating the lagging strand and G templating the leading strand; (iv) there is a strong bias for transition mutations to occur at 5'ApC3'/3'TpG5' sites (where bases 5'A and 3'T are mutated) and, to a lesser extent, at 5'GpC3'/3'CpG5' sites (where bases 5'G and 3'C are mutated); (v) although the rate of small (≤4 nt) insertions and deletions is high at repeat sequences, these events occur at only 1/10th the genomic rate of base-pair substitutions. MMR activity is genetically regulated, and bacteria isolated from nature often lack MMR capacity, suggesting that modulation of MMR can be adaptive. Thus, comparing results from the wild-type and MMR-defective strains may lead to a deeper understanding of factors that determine mutation rates and spectra, how these factors may differ among organisms, and how they may be shaped by environmental conditions.