The replication-related organization of bacterial genomes

Differentially used codons among essential genes in bacteria identified by machine learning-based analysis

Codon usage bias (CUB), the uneven usage of synonymous codons encoding the same amino acid, differs among genes within and across bacteria genomes. CUB is known to be influenced by gene expression and accordingly, CUB differs between the high-expression and low-expression genes in several bacteria. In this article, we have extended codon usage study considering gene essentiality as a feature. Using machine learning (ML) based approaches, we have analysed Relative Synonymous Codon Usage (RSCU) values between essential and non-essential genes in Escherichia coli and thirty-four other bacterial genomes whose gene essentiality features were available in public databases. We observed significant differences in codon usage patterns between essential and non-essential genes for majority of the bacterial genomes and accordingly, ML based classifiers achieved high area under curve (AUC) scores, with a minimum score of 70.0 across twenty-eight organisms. Further, importance of the codons towards classifying genes found to differ among the codons in each genome. Arg codon CGT and Gly codon GGT were observed to be the most preferred codons among essential genes in Escherichia coli. Interestingly, some of the codons like CGT, ATA, GGT and GGG observed to be contributing consistently towards classifying essential genes across thirty-five bacteria genomes studied. In other hand, codons TGY and CAY encoding amino acids Cys and His respectively were among the least contributing codons towards classification among all these bacteria. This study demonstrates the gene essentiality based differences in synonymous codon usage in bacteria genomes and presents a common codon usage pattern across bacteria.

Structural Differences between the Genomes of Deinococcus radiodurans Strains from Different Laboratories

The bacterium Deinococcus radiodurans is known to efficiently and accurately reassemble its genome after hundreds of DNA double-strand breaks (DSBs). Only at very large amounts of radiation-induced DSBs is this accuracy affected in the wild-type D. radiodurans, causing rearrangements in its genome structure. However, changes in its genome structure may also be possible during the propagation and storage of cell cultures. We investigate this possibility by listing structural differences between three completely sequenced genomes of D. radiodurans strains with a recent common ancestor-the type strain stored and sequenced in two different laboratories (of the ATCC 13939 lineage) and the first sequenced strain historically used as the reference (ATCC BAA-816). We detected a number of structural differences and found the most likely mechanisms behind them: (i) transposition/copy number change in mobile interspersed repeats-insertion sequences and small non-coding repeats, (ii) variable number of monomers within tandem repeats, (iii) deletions between long direct DNA repeats, and (iv) deletions between short (4-10 bp) direct DNA repeats. The most surprising finding was the deletions between short repeats because it indicates the utilization of a less accurate DSB repair mechanism in conditions in which a more accurate one should be both available and preferred. The detected structural differences, as well as SNPs and short indels, while being important footprints of deinococcal DNA metabolism and repair, are also a valuable resource for researchers using these D. radiodurans strains.

On the evolution of chromosomal regions with high gene strand bias in bacteria

On circular bacterial chromosomes, the majority of genes are coded on the leading strand. This gene strand bias (GSB) ranges from up to 85% in some Bacillota to a little more than 50% in other phyla. The factors determining the extent of the strand bias remain to be found. Here, we report that species in the phylum Gemmatimonadota share a unique chromosome architecture, distinct from neighboring phyla: in a conserved 600-kb region around the terminus of replication, almost all genes were located on the leading strands, while on the remaining part of the chromosome, the strand preference was more balanced. The high strand bias (HSB) region harbors the rRNA clusters, core, and highly expressed genes. Selective pressure for reduction of collisions with DNA replication to minimize detrimental mutations can explain the conservation of essential genes in this region. Repetitive and mobile elements are underrepresented, suggesting reduced recombination frequency by structural isolation from other parts of the chromosome. We propose that the HSB region forms a distinct chromosomal domain. Gemmatimonadota chromosomes evolved mainly by expansion through horizontal gene transfer and duplications outside of the ancient high strand bias region. In support of our hypothesis, we could further identify two Spiroplasma strains on a similar evolutionary path.IMPORTANCEOn bacterial chromosomes, a preferred location of genes on the leading strand has evolved to reduce conflicts between replication and transcription. Despite a vast body of research, the question why bacteria show large differences in their gene strand bias is still not solved. The discovery of "hybrid" chromosomes in different phyla, including Gemmatimonadota, in which a conserved high strand bias is found exclusively in a region at ter, points toward a role of nucleoid structure, additional to replication, in the evolution of strand preferences. A fine-grained structural analysis of the ever-increasing number of available bacterial genomes could help to better understand the forces that shape the sequential and spatial organization of the cell's information content.

Origin, evolution, and maintenance of gene-strand bias in bacteria

Gene-strand bias is a characteristic feature of bacterial genome organization wherein genes are preferentially encoded on the leading strand of replication, promoting co-orientation of replication and transcription. This co-orientation bias has evolved to protect gene essentiality, expression, and genomic stability from the harmful effects of head-on replication-transcription collisions. However, the origin, variation, and maintenance of gene-strand bias remain elusive. Here, we reveal that the frequency of inversions that alter gene orientation exhibits large variation across bacterial populations and negatively correlates with gene-strand bias. The density, distance, and distribution of inverted repeats show a similar negative relationship with gene-strand bias explaining the heterogeneity in inversions. Importantly, these observations are broadly evident across the entire bacterial kingdom uncovering inversions and inverted repeats as primary factors underlying the variation in gene-strand bias and its maintenance. The distinct catalytic subunits of replicative DNA polymerase have co-evolved with gene-strand bias, suggesting a close link between replication and the origin of gene-strand bias. Congruently, inversion frequencies and inverted repeats vary among bacteria with different DNA polymerases. In summary, we propose that the nature of replication determines the fitness cost of replication-transcription collisions, establishing a selection gradient on gene-strand bias by fine-tuning DNA sequence repeats and, thereby, gene inversions.

Cost-effective in-house COVID-19 reverse transcription-polymerase chain reaction testing with yeast-derived Taq polymerase

BACKGROUND

Despite the decline of the COVID-19 pandemic, there continues to be a persistent requirement for reliable testing methods that can be adapted to future outbreaks and areas with limited resources. While the standard approach of using reverse transcription-polymerase chain reaction (RT-PCR) with Taq polymerase is effective, it faces challenges such as limited access to high-quality enzymes and the presence of bacterial DNA contamination in commercial kits, which can impact the accuracy of test results.

METHODS

This study investigates the production of recombinant Taq polymerase in yeast cells and assesses its crude lysate in a multiplex RT-PCR assay for detecting the SARS-CoV-2 RNA-dependent RNA polymerase (RdRP) and N genes, with human Ribonuclease P serving as an internal control.

RESULTS

The unpurified yeast Taq polymerase demonstrates sensitivity comparable to commercially purified bacterial Taq polymerase and unpurified bacterial counterparts in detecting the RdRP and N genes. It exhibits the highest specificity, with 100% accuracy, for the N gene. The specificity for the RdRP gene closely aligns with that of commercially purified bacterial Taq polymerase and unpurified bacterial Taq polymerase.

CONCLUSIONS

The use of unpurified recombinant yeast Taq polymerase shows promise as a cost-effective approach for conducting in-house COVID-19 RT-PCR testing. By eliminating the need for chromatography purification steps, the production of RT-PCR kits can be streamlined, potentially improving accessibility and scalability, especially in resource-limited settings and future pandemics.

Genome evolution of Buchnera aphidicola (Gammaproteobacteria): Insights into strand compositional asymmetry, codon usage bias, and phylogenetic implications

Taxa of Buchnera aphidicola (hereafter "Buchnera") are mutualistic intracellular symbionts of aphids, known for their remarkable biological traits such as genome reduction, strand compositional asymmetry, and symbiont-host coevolution. With the growing availability of genomic data, we performed a comprehensive analysis of 103 genomes of Buchnera strains from 12 host subfamilies, focusing on the genomic characterizations, codon usage patterns, and phylogenetic implications. Our findings revealed consistent features among all genomes, including small genome sizes, low GC contents, and gene losses. We also identified strong strand compositional asymmetries in all strains at the genome level. Further investigation suggested that mutation pressure may have played a crucial role in shaping codon usage of Buchnera. Moreover, the genomic asymmetries were reflected in asymmetric codon usage preferences within chromosomal genes. Notably, the levels of these asymmetries were varied among strains and were significantly influenced by the degrees of genome shrinkages. Lastly, our phylogenetic analyses presented an alternative topology of Aphididae, based on the Buchnera symbionts, providing robust confirmation of the paraphylies of Eriosomatinae, and Macrosiphini. Our objectives are to further understand the strand compositional asymmetry and codon usage bias of Buchnera taxa, and provide new perspectives for phylogenetic studies of Aphididae.

Comparative genomics analysis of Streptococcus iniae isolated from Trachinotus ovatus: novel insight into antimicrobial resistance and virulence differentiation

BACKGROUND

Streptococcus iniae is an important fish pathogen that cause significant economic losses to the global aquaculture industry every year. Although there have some reports on the genotype of S.iniae and its relationship with virulence, no genome-scale comparative analysis has been performed so far. In our previous work, we characterized 17 isolates of S.iniae from Trachinotus ovatus and divided them into two genotypes using RAPD and rep-PCR methods. Among them, BH15-2 was classified as designated genotype A (in RAPD) and genotype 1 (in rep-PCR), while BH16-24 was classified as genotype B and genotype 2. Herein, we compared the differences in growth, drug resistance, virulence, and genome between BH15-2 and BH16-24.

RESULTS

The results showed that the growth ability of BH16-24 was significantly faster than that of BH15-2 at the exponential stage. Antimicrobial tests revealed that BH15-2 was susceptible to most of the tested antibiotics except neomycin and gentamycin. In contrast, BH16-24 was resistant to 7 antibiotics including penicillin, sulfasomizole, compound sulfamethoxazole tablets, polymyxin B, spectinomycin, rifampin and ceftazidime. Intraperitoneal challenge of T.ovatus, showed that the LD50 value of BH15-2 was 4.0 × 102 CFU/g, while that of BH16-24 was 1.2 × 105 CFU/g. The genome of S.iniae BH15-2 was 2,175,659 bp with a GC content of 36.80%. Meanwhile, the genome of BH16-24 was 2,153,918 bp with a GC content of 36.83%. Comparative genome analysis indicated that compared with BH15-2, BH16-24 genome had a large-scale genomic inversion fragment, at the location from 502,513 bp to 1,788,813 bp, resulting in many of virulence and resistance genes differentially expression. In addition, there was a 46 kb length, intact phage sequence in BH15-2 genome, which was absent in BH16-24.

CONCLUSION

Comparative genomic studies of BH15-2 and BH16-24 showed that the main difference is a 1.28 Mbp inversion fragment. The inversion fragment may lead to abnormal expression of drug resistant and virulence genes, which is believed to be the main reason for the multiple resistance and weakened virulence of BH16-24. Our study revealed the potential mechanisms in underlying the differences of multidrug resistance and virulence among different genotypes of S.iniae.

Potential of Selected African Medicinal Plants as Alternative Therapeutics against Multi-Drug-Resistant Bacteria

Antimicrobial resistance is considered a "One-Health" problem, impacting humans, animals, and the environment. The problem of the rapid development and spread of bacteria resistant to multiple antibiotics is a rising global health threat affecting both rich and poor nations. Low- and middle-income countries are at highest risk, in part due to the lack of innovative research on the surveillance and discovery of novel therapeutic options. Fast and effective drug discovery is crucial towards combatting antimicrobial resistance and reducing the burden of infectious diseases. African medicinal plants have been used for millennia in folk medicine to cure many diseases and ailments. Over 10% of the Southern African vegetation is applied in traditional medicine, with over 15 species being partially or fully commercialized. These include the genera Euclea, Ficus, Aloe, Lippia. And Artemisia, amongst many others. Bioactive compounds from indigenous medicinal plants, alone or in combination with existing antimicrobials, offer promising solutions towards overcoming multi-drug resistance. Secondary metabolites have different mechanisms and modes of action against bacteria, such as the inhibition and disruption of cell wall synthesis; inhibition of DNA replication and ATP synthesis; inhibition of quorum sensing; inhibition of AHL or oligopeptide signal generation, broadcasting, and reception; inhibition of the formation of biofilm; disruption of pathogenicity activities; and generation of reactive oxygen species. The aim of this review is to highlight some promising traditional medicinal plants found in Africa and provide insights into their secondary metabolites as alternative options in antibiotic therapy against multi-drug-resistant bacteria. Additionally, synergism between plant secondary metabolites and antibiotics has been discussed.

Expression of Microcystis Biosynthetic Gene Clusters in Natural Populations Suggests Temporally Dynamic Synthesis of Novel and Known Secondary Metabolites in Western Lake Erie

Microcystis spp. produce diverse secondary metabolites within freshwater cyanobacterial harmful algal blooms (cyanoHABs) around the world. In addition to the biosynthetic gene clusters (BGCs) encoding known compounds, Microcystis genomes harbor numerous BGCs of unknown function, indicating a poorly understood chemical repertoire. While recent studies show that Microcystis produces several metabolites in the lab and field, little work has focused on analyzing the abundance and expression of its broader suite of BGCs during cyanoHAB events. Here, we use metagenomic and metatranscriptomic approaches to track the relative abundance of Microcystis BGCs and their transcripts throughout the 2014 western Lake Erie cyanoHAB. The results indicate the presence of several transcriptionally active BGCs that are predicted to synthesize both known and novel secondary metabolites. The abundance and expression of these BGCs shifted throughout the bloom, with transcript abundance levels correlating with temperature, nitrate, and phosphorus concentrations and the abundance of co-occurring predatory and competitive eukaryotic microorganisms, suggesting the importance of both abiotic and biotic controls in regulating expression. This work highlights the need for understanding the chemical ecology and potential risks to human and environmental health posed by secondary metabolites that are produced but often unmonitored. It also indicates the prospects for identifying pharmaceutical-like molecules from cyanoHAB-derived BGCs. IMPORTANCE Microcystis spp. dominate cyanobacterial harmful algal blooms (cyanoHABs) worldwide and pose significant threats to water quality through the production of secondary metabolites, many of which are toxic. While the toxicity and biochemistry of microcystins and several other compounds have been studied, the broader suite of secondary metabolites produced by Microcystis remains poorly understood, leaving gaps in our understanding of their impacts on human and ecosystem health. We used community DNA and RNA sequences to track the diversity of genes encoding synthesis of secondary metabolites in natural Microcystis populations and assess patterns of transcription in western Lake Erie cyanoHABs. Our results reveal the presence of both known gene clusters that encode toxic secondary metabolites as well as novel ones that may encode cryptic compounds. This research highlights the need for targeted studies of the secondary metabolite diversity in western Lake Erie, a vital freshwater source to the United States and Canada.

Chromosomal Position of Ribosomal Protein Genes Affects Long-Term Evolution of Vibrio cholerae

It is unclear how gene order within the chromosome influences genome evolution. Bacteria cluster transcription and translation genes close to the replication origin (oriC). In Vibrio cholerae, relocation of s10-spc-α locus (S10), the major locus of ribosomal protein genes, to ectopic genomic positions shows that its relative distance to the oriC correlates to a reduction in growth rate, fitness, and infectivity. To test the long-term impact of this trait, we evolved 12 populations of V. cholerae strains bearing S10 at an oriC-proximal or an oriC-distal location for 1,000 generations. During the first 250 generations, positive selection was the main force driving mutation. After 1,000 generations, we observed more nonadaptative mutations and hypermutator genotypes. Populations fixed inactivating mutations at many genes linked to virulence: flagellum, chemotaxis, biofilm, and quorum sensing. Throughout the experiment, all populations increased their growth rates. However, those bearing S10 close to oriC remained the fittest, indicating that suppressor mutations cannot compensate for the genomic position of the main ribosomal protein locus. Selection and sequencing of the fastest-growing clones allowed us to characterize mutations inactivating, among other sites, flagellum master regulators. Reintroduction of these mutations into the wild-type context led to a ≈10% growth improvement. In conclusion, the genomic location of ribosomal protein genes conditions the evolutionary trajectory of V. cholerae. While genomic content is highly plastic in prokaryotes, gene order is an underestimated factor that conditions cellular physiology and evolution. A lack of suppression enables artificial gene relocation as a tool for genetic circuit reprogramming. IMPORTANCE The bacterial chromosome harbors several entangled processes such as replication, transcription, DNA repair, and segregation. Replication begins bidirectionally at the replication origin (oriC) until the terminal region (ter) organizing the genome along the ori-ter axis gene order along this axis could link genome structure to cell physiology. Fast-growing bacteria cluster translation genes near oriC. In Vibrio cholerae, moving them away was feasible but at the cost of losing fitness and infectivity. Here, we evolved strains harboring ribosomal genes close or far from oriC. Growth rate differences persisted after 1,000 generations. No mutation was able to compensate for the growth defect, showing that ribosomal gene location conditions their evolutionary trajectory. Despite the high plasticity of bacterial genomes, evolution has sculpted gene order to optimize the ecological strategy of the microorganism. We observed growth rate improvement throughout the evolution experiment that occurred at expense of energetically costly processes such the flagellum biosynthesis and virulence-related functions. From the biotechnological point of view, manipulation of gene order enables altering bacterial growth with no escape events.

Form and function of archaeal genomes

A key maxim in modernist architecture is that 'form follows function'. While modernist buildings are hopefully the product of intelligent design, the architectures of chromosomes have been sculpted by the forces of evolution over many thousands of generations. In the following, I will describe recent advances in our understanding of chromosome architecture in the archaeal domain of life. Although much remains to be learned about the mechanistic details of archaeal chromosome organization, some general principles have emerged. At the 10-100 kb level, archaeal chromosomes have a conserved local organization reminiscent of bacterial genomes. In contrast, lineage-specific innovations appear to have imposed distinct large-scale architectural features. The ultimate functions of genomes are to store and to express genetic information. Gene expression profiles have been shown to influence chromosome architecture, thus their form follows function. However, local changes to chromosome conformation can also influence gene expression and therefore, in these instances, function follows form.

Patterns of abundance, chromosomal localization, and domain organization among c-di-GMP-metabolizing genes revealed by comparative genomics of five alphaproteobacterial orders

BACKGROUND

Bis-(3'-5')-cyclic dimeric guanosine monophosphate (c-di-GMP) is a bacterial second messenger that affects diverse processes in different bacteria, including the cell cycle, motility, and biofilm formation. Its cellular levels are controlled by the opposing activities of two types of enzymes, with synthesis by diguanylate cyclases containing a GGDEF domain and degradation by phosphodiesterases containing either an HD-GYP or an EAL domain. These enzymes are ubiquitous in bacteria with up to 50 encoded in some genomes, the specific functions of which are mostly unknown.

RESULTS

We used comparative analyses to identify genomic patterns among genes encoding proteins with GGDEF, EAL, and HD-GYP domains in five orders of the class Alphaproteobacteria. GGDEF-containing sequences and GGDEF-EAL hybrids were the most abundant and had the highest diversity of co-occurring auxiliary domains while EAL and HD-GYP containing sequences were less abundant and less diverse with respect to auxiliary domains. There were striking patterns in the chromosomal localizations of the genes found in two of the orders. The Rhodobacterales' EAL-encoding genes and Rhizobiales' GGDEF-EAL-encoding genes showed opposing patterns of distribution compared to the GGDEF-encoding genes. In the Rhodobacterales, the GGDEF-encoding genes showed a tri-modal distribution with peaks mid-way between the origin (ori) and terminus (ter) of replication and at ter while the EAL-encoding genes peaked near ori. The patterns were more complex in the Rhizobiales, but the GGDEF-encoding genes were biased for localization near ter.

CONCLUSIONS

The observed patterns in the chromosomal localizations of these genes suggest a coupling of synthesis and hydrolysis of c-di-GMP with the cell cycle. Moreover, the higher proportions and diversities of auxiliary domains associated with GGDEF domains and GGDEF-EAL hybrids compared to EAL or HD-GYP domains could indicate that more stimuli affect synthesis compared to hydrolysis of c-di-GMP.

Shared properties of gene transfer agent and core genes revealed by comparative genomics of Alphaproteobacteria

Gene transfer agents (GTAs) are phage-like particles that transfer pieces of cellular genomic DNA to other cells. Homologues of the Rhodobacter capsulatus GTA (RcGTA) structural genes are widely distributed in the Alphaproteobacteria and particularly well conserved in the order Rhodobacterales. Possible reasons for their widespread conservation are still being discussed. It has been suggested that these alphaproteobacterial elements originate from a prophage that was present in an ancestral bacterium and subsequently evolved into a GTA that is now widely maintained in extant descendant lineages. Here, we analysed genomic properties that might relate to the conservation of these alphaproteobacterial GTAs. This revealed that the chromosomal locations of the GTA gene clusters are biased. They primarily occur on the leading strand of DNA replication, at large distances from long repetitive elements, and thus are in regions of lower plasticity, and in areas of extreme GC skew, which also accumulate core genes. These extreme GC skew regions arise from the preferential use of codons with an excess of G over C, a distinct phenomenon from the elevated GC content that has previously been found to be associated with GTA genes. The observed properties, along with their high level of conservation, show that GTA genes share multiple features with core genes in the examined lineages of the Alphaproteobacteria.

The Impact of RNA-DNA Hybrids on Genome Integrity in Bacteria

During the essential processes of DNA replication and transcription, RNA-DNA hybrid intermediates are formed that pose significant risks to genome integrity when left unresolved. To manage RNA-DNA hybrids, all cells rely on RNase H family enzymes that specifically cleave the RNA portion of the many different types of hybrids that form in vivo. Recent experimental advances have provided new insight into how RNA-DNA hybrids form and the consequences to genome integrity that ensue when persistent hybrids remain unresolved. Here we review the types of RNA-DNA hybrids, including R-loops, RNA primers, and ribonucleotide misincorporations, that form during DNA replication and transcription and discuss how each type of hybrid can contribute to genome instability in bacteria. Further, we discuss how bacterial RNase HI, HII, and HIII and bacterial FEN enzymes contribute to genome maintenance through the resolution of hybrids.

RecBCD enzyme and Chi recombination hotspots as determinants of self vs. non-self: Myths and mechanisms

Bacteria face a challenge when DNA enters their cells by transformation, mating, or phage infection. Should they treat this DNA as an invasive foreigner and destroy it, or consider it one of their own and potentially benefit from incorporating new genes or alleles to gain useful functions? It is frequently stated that the short nucleotide sequence Chi (5' GCTGGTGG 3'), a hotspot of homologous genetic recombination recognized by Escherichia coli's RecBCD helicase-nuclease, allows E. coli to distinguish its DNA (self) from any other DNA (non-self) and to destroy non-self DNA, and that Chi is "over-represented" in the E. coli genome. We show here that these latter statements (dogmas) are not supported by available evidence. We note Chi's wide-spread occurrence and activity in distantly related bacterial species and phages. We illustrate multiple, highly non-random features of the genomes of E. coli and coliphage P1 that account for Chi's high frequency and genomic position, leading us to propose that P1 selects for Chi's enhancement of recombination, whereas E. coli selects for the preferred codons in Chi. We discuss other, previously described mechanisms for self vs. non-self determination involving RecBCD and for RecBCD's destruction of DNA that cannot recombine, whether foreign or domestic, with or without Chi.

Replication-Dependent Organization Constrains Positioning of Long DNA Repeats in Bacterial Genomes

Bacterial genome organization is primarily driven by chromosomal replication from a single origin of replication. However, chromosomal rearrangements, which can disrupt such organization, are inevitable in nature. Long DNA repeats are major players mediating rearrangements, large and small, via homologous recombination. Since changes to genome organization affect bacterial fitness-and more so in fast-growing than slow-growing bacteria-and are under selection, it is reasonable to expect that genomic positioning of long DNA repeats is also under selection. To test this, we identified identical DNA repeats of at least 100 base pairs across ∼6,000 bacterial genomes and compared their distribution in fast- and slow-growing bacteria. We found that long identical DNA repeats are distributed in a non-random manner across bacterial genomes. Their distribution differs in the overall number, orientation, and proximity to the origin of replication, between fast- and slow-growing bacteria. We show that their positioning-which might arise from a combination of the processes that produce repeats and selection on rearrangements that recombination between repeat elements might cause-permits less disruption to the replication-dependent genome organization of bacteria compared with random suggesting it as a major constraint to positioning of long DNA repeats.

Metagenomic and Metatranscriptomic Insights into Population Diversity of Microcystis Blooms: Spatial and Temporal Dynamics of mcy Genotypes, Including a Partial Operon That Can Be Abundant and Expressed

Cyanobacterial harmful algal blooms (cyanoHABs) degrade freshwater ecosystems globally. Microcystis aeruginosa often dominates cyanoHABs and produces microcystin (MC), a class of hepatotoxins that poses threats to human and animal health. Microcystin toxicity is influenced by distinct structural elements across a diversity of related molecules encoded by variant mcy operons. However, the composition and distribution of mcy operon variants in natural blooms remain poorly understood. Here, we characterized the variant composition of mcy genes in western Lake Erie Microcystis blooms from 2014 and 2018. Sampling was conducted across several spatial and temporal scales, including different bloom phases within 2014, extensive spatial coverage on the same day (2018), and frequent, autonomous sampling over a 2-week period (2018). Mapping of metagenomic and metatranscriptomic sequences to reference sequences revealed three Microcystis mcy genotypes: complete (all genes present [mcyA-J]), partial (truncated mcyA, complete mcyBC, and missing mcyD-J), and absent (no mcy genes). We also detected two different variants of mcyB that may influence the production of microcystin congeners. The relative abundance of these genotypes was correlated with pH and nitrate concentrations. Metatranscriptomic analysis revealed that partial operons were, at times, the most abundant genotype and expressed in situ, suggesting the potential biosynthesis of truncated products. Quantification of genetic divergence between genotypes suggests that the observed strains are the result of preexisting heterogeneity rather than de novo mutation during the sampling period. Overall, our results show that natural Microcystis populations contain several cooccurring mcy genotypes that dynamically shift in abundance spatiotemporally via strain succession and likely influence the observed diversity of the produced congeners. IMPORTANCE Cyanobacteria are responsible for producing microcystins (MCs), a class of potent and structurally diverse toxins, in freshwater systems around the world. While microcystins have been studied for over 50 years, the diversity of their chemical forms and how this variation is encoded at the genetic level remain poorly understood, especially within natural populations of cyanobacterial harmful algal blooms (cyanoHABs). Here, we leverage community DNA and RNA sequences to track shifts in mcy genes responsible for producing microcystin, uncovering the relative abundance, expression, and variation of these genes. We studied this phenomenon in western Lake Erie, which suffers annually from cyanoHAB events, with impacts on drinking water, recreation, tourism, and commercial fishing.

Identification of Group II Intron RmInt1 Binding Sites in a Bacterial Genome

RmInt1 is a group II intron encoding a reverse transcriptase protein (IEP) lacking the C-terminal endonuclease domain. RmInt1 is an efficient mobile retroelement that predominantly reverse splices into the transient single-stranded DNA at the template for lagging strand DNA synthesis during host replication, a process facilitated by the interaction of the RmInt1 IEP with DnaN at the replication fork. It has been suggested that group II intron ribonucleoprotein particles bind DNA nonspecifically, and then scan for their correct target site. In this study, we investigated RmInt1 binding sites throughout the Sinorhizobium meliloti genome, by chromatin-immunoprecipitation coupled with next-generation sequencing. We found that RmInt1 binding sites cluster around the bidirectional replication origin of each of the three replicons comprising the S. meliloti genome. Our results provide new evidence linking group II intron mobility to host DNA replication.

Multiple classes and isoforms of the RNA polymerase recycling motor protein HelD

Efficient control of transcription is essential in all organisms. In bacteria, where DNA replication and transcription occur simultaneously, the replication machinery is at risk of colliding with highly abundant transcription complexes. This can be exacerbated by the fact that transcription complexes pause frequently. When pauses are long-lasting, the stalled complexes must be removed to prevent collisions with either another transcription complex or the replication machinery. HelD is a protein that represents a new class of ATP-dependent motor proteins distantly related to helicases. It was first identified in the model Gram-positive bacterium Bacillus subtilis and is involved in removing and recycling stalled transcription complexes. To date, two classes of HelD have been identified: one in the low G+C and the other in the high G+C Gram-positive bacteria. In this work, we have undertaken the first comprehensive investigation of the phylogenetic diversity of HelD proteins. We show that genes in certain bacterial classes have been inherited by horizontal gene transfer, many organisms contain multiple expressed isoforms of HelD, some of which are associated with antibiotic resistance, and that there is a third class of HelD protein found in Gram-negative bacteria. In summary, HelD proteins represent an important new class of transcription factors associated with genome maintenance and antibiotic resistance that are conserved across the Eubacterial kingdom.

DNA Replication-Transcription Conflicts Do Not Significantly Contribute to Spontaneous Mutations Due to Replication Errors in Escherichia coli

Encounters between DNA replication and transcription can cause genomic disruption, particularly when the two meet head-on. Whether these conflicts produce point mutations is debated. This paper presents detailed analyses of a large collection of mutations generated during mutation accumulation experiments with mismatch repair (MMR)-defective Escherichia coli. With MMR absent, mutations are primarily due to DNA replication errors. Overall, there were no differences in the frequencies of base pair substitutions or small indels (i.e., insertion and deletions of ≤4 bp) in the coding sequences or promoters of genes oriented codirectionally versus head-on to replication. Among a subset of highly expressed genes, there was a 2- to 3-fold bias for indels in genes oriented head-on to replication, but this difference was almost entirely due to the asymmetrical genomic locations of tRNA genes containing mononucleotide runs, which are hot spots for indels. No additional orientation bias in mutation frequencies occurred when MMR⁻ strains were also defective for transcription-coupled repair (TCR). However, in contrast to other reports, loss of TCR slightly increased the overall mutation rate, meaning that TCR is antimutagenic. There was no orientation bias in mutation frequencies among the stress response genes that are regulated by RpoS or induced by DNA damage. Thus, biases in the locations of mutational targets can account for most, if not all, apparent biases in mutation frequencies between genes oriented head-on versus codirectional to replication. In addition, the data revealed a strong correlation of the frequency of base pair substitutions with gene length but no correlation with gene expression levels.

GRINS: Genetic elements that recode assembly-line polyketide synthases and accelerate their diversification

Assembly-line polyketide synthases (PKSs) are large and complex enzymatic machineries with a multimodular architecture, typically encoded in bacterial genomes by biosynthetic gene clusters. Their modularity has led to an astounding diversity of biosynthesized molecules, many with medical relevance. Thus, understanding the mechanisms that drive PKS evolution is fundamental for both functional prediction of natural PKSs as well as for the engineering of novel PKSs. Here, we describe a repetitive genetic element in assembly-line PKS genes which appears to play a role in accelerating the diversification of closely related biosynthetic clusters. We named this element GRINS: genetic repeats of intense nucleotide skews. GRINS appear to recode PKS protein regions with a biased nucleotide composition and to promote gene conversion. GRINS are present in a large number of assembly-line PKS gene clusters and are particularly widespread in the actinobacterial genus Streptomyces While the molecular mechanisms associated with GRINS appearance, dissemination, and maintenance are unknown, the presence of GRINS in a broad range of bacterial phyla and gene families indicates that these genetic elements could play a fundamental role in protein evolution.

Characterization of gross genome rearrangements in Deinococcus radiodurans recA mutants

Genome stability in radioresistant bacterium Deinococcus radiodurans depends on RecA, the main bacterial recombinase. Without RecA, gross genome rearrangements occur during repair of DNA double-strand breaks. Long repeated (insertion) sequences have been identified as hot spots for ectopic recombination leading to genome rearrangements, and single-strand annealing (SSA) postulated to be the most likely mechanism involved in this process. Here, we have sequenced five isolates of D. radiodurans recA mutant carrying gross genome rearrangements to precisely characterize the rearrangements and to elucidate the underlying repair mechanism. The detected rearrangements consisted of large deletions in chromosome II in all the sequenced recA isolates. The mechanism behind these deletions clearly differs from the classical SSA; it utilized short (4-11 bp) repeats as opposed to insertion sequences or other long repeats. Moreover, it worked over larger linear DNA distances from those previously tested. Our data are most compatible with alternative end-joining, a recombination mechanism that operates in eukaryotes, but is also found in Escherichia coli. Additionally, despite the recA isolates being preselected for different rearrangement patterns, all identified deletions were found to overlap in a 35 kb genomic region. We weigh the evidence for mechanistic vs. adaptive reasons for this phenomenon.

Connection Between Chromosomal Location and Function of CtrA Phosphorelay Genes in Alphaproteobacteria

Most bacterial chromosomes are circular, with replication starting at one origin (ori) and proceeding on both replichores toward the terminus (ter). Several studies have shown that the location of genes relative to ori and ter can have profound effects on regulatory networks and physiological processes. The CtrA phosphorelay is a gene regulatory system conserved in most alphaproteobacteria. It was first discovered in Caulobacter crescentus where it controls replication and division into a stalked and a motile cell in coordination with other factors. The locations of the ctrA gene and targets of this response regulator on the chromosome affect their expression through replication-induced DNA hemi-methylation and specific positioning along a CtrA activity gradient in the dividing cell, respectively. Here we asked to what extent the location of CtrA regulatory network genes might be conserved in the alphaproteobacteria. We determined the locations of the CtrA phosphorelay and associated genes in closed genomes with unambiguously identifiable ori from members of five alphaproteobacterial orders. The location of the phosphorelay genes was the least conserved in the Rhodospirillales followed by the Sphingomonadales. In the Rhizobiales a trend toward certain chromosomal positions could be observed. Compared to the other orders, the CtrA phosphorelay genes were conserved closer to ori in the Caulobacterales. In contrast, the genes were highly conserved closer to ter in the Rhodobacterales. Our data suggest selection pressure results in differential positioning of CtrA phosphorelay and associated genes in alphaproteobacteria, particularly in the orders Rhodobacterales, Caulobacterales and Rhizobiales that is worth deeper investigation.

Understanding the codon usage patterns of mitochondrial CO genes among Amphibians

A universal phenomenon of using synonymous codons unequally in coding sequences known as codon usage bias (CUB) is observed in all forms of life. Mutation and natural selection drive CUB in many species but the relative role of evolutionary forces varies across species, genes and genomes. We studied the CUB in mitochondrial (mt) CO genes from three orders of Amphibia using bioinformatics approach as no work was reported yet. We observed that CUB of mt CO genes of Amphibians was weak across different orders. Order Caudata had higher CUB followed by Gymnophiona and Anura for all genes and CUB also varied across genes. Nucleotide composition analysis showed that CO genes were AT-rich. The AT content in Caudata was higher than that in Gymnophiona while Anura showed the least content. Multiple investigations namely nucleotide composition, correspondence analysis, parity plot analysis showed that the interplay of mutation pressure and natural selection caused CUB in these genes. Neutrality plot suggested the involvement of natural selection was more than the mutation pressure. The contribution of natural selection was higher in Anura than Gymnophiona and the lowest in Caudata. The codons CGA, TGA, AAA were found to be highly favoured by nature across all genes and orders.

Evolutionary Trajectory of the Replication Mode of Bacterial Replicons

Chromosome replication is an essential process for cell division. The mode of chromosome replication has important impacts on the structure of the chromosome and replication speed.

As typical bacterial replicons, circular chromosomes replicate bidirectionally and circular plasmids replicate either bidirectionally or unidirectionally. Whereas the finding of chromids (plasmid-derived chromosomes) in multiple bacterial lineages provides circumstantial evidence that chromosomes likely evolved from plasmids, all experimentally assayed chromids were shown to use bidirectional replication. Here, we employed a model system, the marine bacterial genus Pseudoalteromonas, members of which consistently carry a chromosome and a chromid. We provide experimental and bioinformatic evidence that while chromids in a few strains replicate bidirectionally, most replicate unidirectionally. This is the first experimental demonstration of the unidirectional replication mode in bacterial chromids. Phylogenomic and comparative genomic analyses showed that the bidirectional replication evolved only once from a unidirectional ancestor and that this transition was associated with insertions of exogenous DNA and relocation of the replication terminus region (ter2) from near the origin site (ori2) to a position roughly opposite it. This process enables a plasmid-derived chromosome to increase its size and expand the bacterium’s metabolic versatility while keeping its replication synchronized with that of the main chromosome. A major implication of our study is that the uni- and bidirectionally replicating chromids may represent two stages on the evolutionary trajectory from unidirectionally replicating plasmids to bidirectionally replicating chromosomes in bacteria. Further bioinformatic analyses predicted unidirectionally replicating chromids in several unrelated bacterial phyla, suggesting that evolution from unidirectionally to bidirectionally replicating replicons occurred multiple times in bacteria.

Molecular basis for RNA polymerase-dependent transcription complex recycling by the helicase-like motor protein HelD

In bacteria, transcription complexes stalled on DNA represent a major source of roadblocks for the DNA replication machinery that must be removed in order to prevent damaging collisions. Gram-positive bacteria contain a transcription factor HelD that is able to remove and recycle stalled complexes, but it was not known how it performed this function. Here, using single particle cryo-electron microscopy, we have determined the structures of Bacillus subtilis RNA polymerase (RNAP) elongation and HelD complexes, enabling analysis of the conformational changes that occur in RNAP driven by HelD interaction. HelD has a 2-armed structure which penetrates deep into the primary and secondary channels of RNA polymerase. One arm removes nucleic acids from the active site, and the other induces a large conformational change in the primary channel leading to removal and recycling of the stalled polymerase, representing a novel mechanism for recycling transcription complexes in bacteria.

Long inverted repeats around the chromosome replication terminus in the model strain Bacillus thuringiensis serovar israelensis BGSC 4Q7

Bacillus thuringiensis serovar israelensis is the most widely used natural biopesticide against mosquito larvae worldwide. Its lineage has been actively studied and a plasmid-free strain, B. thuringiensis serovar israelensis BGSC 4Q7 (4Q7), has been produced. Previous sequencing of the genome of this strain has revealed the persistent presence of a 235 kb extrachromosomal element, pBtic235, which has been shown to be an inducible prophage, although three putative chromosomal prophages have been lost. Moreover, a 492 kb region, potentially including the standard replication terminus, has also been deleted in the 4Q7 strain, indicating an absence of essential genes in this area. We reanalysed the genome coverage distribution of reads for the previously sequenced variant strain, and sequenced two independently maintained samples of the 4Q7 strain. A 553 kb area, close to the 492 kb deletion, was found to be duplicated. This duplication presumably restored the equal sizes of the replichores, and a balanced functioning of replication termination. An analysis of genome assembly graphs revealed a transient association of the host chromosome with the pBtic235 element. This association may play a functional role in the replication of the bacterial chromosome, and the termination of this process in particular. The genome-restructuring events detected may modify the genetic status of cytotoxic or haemolytic toxins, potentially influencing strain virulence. Twelve of the single-nucleotide variants identified in 4Q7 were probably due to the procedure used for strain construction or were present in the precursor of this strain. No sequence variants were found in pBtic235, but the distribution of the corresponding 4Q7 reads indicates a significant difference from counterparts in natural B. thuringiensis serovar israelensis strains, suggesting a duplication or over-replication in 4Q7. Thus, the 4Q7 strain is not a pure plasmid-less offshoot, but a highly genetically modified derivative of its natural ancestor. In addition to potentially influencing virulence, genome-restructuring events can modify the replication termination machinery. These findings have potential implications for the conclusions of virulence studies on 4Q7 as a model, but they also raise interesting fundamental questions about the functioning of the Bacillus genome.

Comparative Analysis of the Circular and Highly Asymmetrical Marseilleviridae Genomes

Marseilleviridae members are large dsDNA viruses with icosahedral particles 250 nm in diameter infecting Acanthamoeba. Their 340 to 390 kb genomes encode 450 to 550 protein-coding genes. Since the discovery of marseillevirus (the prototype of the family) in 2009, several strains were isolated from various locations, among which 13 are now fully sequenced. This allows the organization of their genomes to be deciphered through comparative genomics. Here, we first experimentally demonstrate that the Marseilleviridae genomes are circular. We then acknowledge a strong bias in sequence conservation, revealing two distinct genomic regions. One gathers most Marseilleviridae paralogs and has undergone genomic rearrangements, while the other, enriched in core genes, exhibits the opposite pattern. Most of the genes whose protein products compose the viral particles are located in the conserved region. They are also strongly biased toward a late gene expression pattern. We finally discuss the potential advantages of Marseilleviridae having a circular genome, and the possible link between the biased distribution of their genes and the transcription as well as DNA replication mechanisms that remain to be characterized.

Vibrionaceae core, shell and cloud genes are non-randomly distributed on Chr 1: An hypothesis that links the genomic location of genes with their intracellular placement

Background

The genome of Vibrionaceae bacteria, which consists of two circular chromosomes, is replicated in a highly ordered fashion. In fast-growing bacteria, multifork replication results in higher gene copy numbers and increased expression of genes located close to the origin of replication of Chr 1 (ori1). This is believed to be a growth optimization strategy to satisfy the high demand of essential growth factors during fast growth. The relationship between ori1-proximate growth-related genes and gene expression during fast growth has been investigated by many researchers. However, it remains unclear which other gene categories that are present close to ori1 and if expression of all ori1-proximate genes is increased during fast growth, or if expression is selectively elevated for certain gene categories.

Results

We calculated the pangenome of all complete genomes from the Vibrionaceae family and mapped the four pangene categories, core, softcore, shell and cloud, to their chromosomal positions. This revealed that core and softcore genes were found heavily biased towards ori1, while shell genes were overrepresented at the opposite part of Chr 1 (i.e., close to ter1). RNA-seq of Aliivibrio salmonicida and Vibrio natriegens showed global gene expression patterns that consistently correlated with chromosomal distance to ori1. Despite a biased gene distribution pattern, all pangene categories contributed to a skewed expression pattern at fast-growing conditions, whereas at slow-growing conditions, softcore, shell and cloud genes were responsible for elevated expression.

Conclusion

The pangene categories were non-randomly organized on Chr 1, with an overrepresentation of core and softcore genes around ori1, and overrepresentation of shell and cloud genes around ter1. Furthermore, we mapped our gene distribution data on to the intracellular positioning of chromatin described for V. cholerae, and found that core/softcore and shell/cloud genes appear enriched at two spatially separated intracellular regions. Based on these observations, we hypothesize that there is a link between the genomic location of genes and their cellular placement.

The roles of replication-transcription conflict in mutagenesis and evolution of genome organization

Replication-transcription conflicts promote mutagenesis and give rise to evolutionary signatures, with fundamental importance to genome stability ranging from bacteria to metastatic cancer cells. This review focuses on the interplay between replication-transcription conflicts and the evolution of gene directionality. In most bacteria, the majority of genes are encoded on the leading strand of replication such that their transcription is co-directional with the direction of DNA replication fork movement. This gene strand bias arises primarily due to negative selection against deleterious consequences of head-on replication-transcription conflict. However, many genes remain head-on. Can head-on orientation provide some benefit? We combine insights from both mechanistic and evolutionary studies, review published work, and analyze gene expression data to evaluate an emerging model that head-on genes are temporal targets for adaptive mutagenesis during stress. We highlight the alternative explanation that genes in the head-on orientation may simply be the result of genomic inversions and relaxed selection acting on nonessential genes. We seek to clarify how the mechanisms of replication-transcription conflict, in concert with other mutagenic mechanisms, balanced by natural selection, have shaped bacterial genome evolution.

Evolutionary advantage of anti-parallel strand orientation of duplex DNA

DNA in all living systems shares common properties that are remarkably well suited to its function, suggesting refinement by evolution. However, DNA also shares some counter-intuitive properties which confer no obvious benefit, such as strand directionality and anti-parallel strand orientation, which together result in the complicated lagging strand replication. The evolutionary dynamics that led to these properties of DNA remain unknown but their universality suggests that they confer as yet unknown selective advantage to DNA. In this article, we identify an evolutionary advantage of anti-parallel strand orientation of duplex DNA, within a given set of plausible premises. The advantage stems from the increased rate of replication, achieved by dividing the DNA into predictable, independently and simultaneously replicating segments, as opposed to sequentially replicating the entire DNA, thereby parallelizing the replication process. We show that anti-parallel strand orientation is essential for such a replicative organization of DNA, given our premises, the most important of which is the assumption of the presence of sequence-dependent asymmetric cooperativity in DNA.

Spatial Patterns of Gene Expression in Bacterial Genomes

Gene expression in bacteria is a remarkably controlled and intricate process impacted by many factors. One such factor is the genomic position of a gene within a bacterial genome. Genes located near the origin of replication generally have a higher expression level, increased dosage, and are often more conserved than genes located farther from the origin of replication. The majority of the studies involved with these findings have only noted this phenomenon in a single gene or cluster of genes that was re-located to pre-determined positions within a bacterial genome. In this work, we look at the overall expression levels from eleven bacterial data sets from Escherichia coli, Bacillus subtilis, Streptomyces, and Sinorhizobium meliloti. We have confirmed that gene expression tends to decrease when moving away from the origin of replication in majority of the replicons analysed in this study. This study sheds light on the impact of genomic location on molecular trends such as gene expression and highlights the importance of accounting for spatial trends in bacterial molecular analysis.

Codon usage trend in genes associated with obesity

Obesity is not only a social menace but also an economic burden as it reduces productivity and increases health care cost. We used bioinformatic tools to analyze the CUB of obesity associated genes and compared with housekeeping genes (control) to explore the similarities and differences between two data sets as no work was reported yet. The mean effective number of codons (ENC) in genes associated with obesity and housekeeping gene was 50.45 and 52.03 respectively, indicating low CUB. The relative synonymous codon usage (RSCU) suggested that codons namely CTG and GTG were over-represented in both obesity and housekeeping genes while under-represented codons were TCG, TTA, CTA, CCG, CAA, CGT, ATA, ACG, GTA and GCG in obesity genes and TCG, TTA, CCG, ATA, ACG, GTA, and GCG in housekeeping genes. t test analysis suggested that 11 codons namely TTA (Leu), TTG (Leu), CCG (Pro), CAC (His), CAA (Gln), CAG (Gln), CGT (Arg), AGA (Arg), ATA (Ile), ATT (Ile) and GCG (Ala) were significantly differed (p < 0.05 or p < 0.01) between obesity and housekeeping genes. Highly significant correlation was observed between GC12 and GC3 in obesity and housekeeping genes i.e. r = 0.580** and r = 0.498** (p < 0.01) respectively indicating the effect of directional mutation pressure present in all codon positions.

Operon-based approach for the inference of rRNA and tRNA evolutionary histories in bacteria

Background

In bacterial genomes, rRNA and tRNA genes are often organized into operons, i.e. segments of closely located genes that share a single promoter and are transcribed as a single unit. Analyzing how these genes and operons evolve can help us understand what are the most common evolutionary events affecting them and give us a better picture of ancestral codon usage and protein synthesis.

Results

We introduce BOPAL, a new approach for the inference of evolutionary histories of rRNA and tRNA genes in bacteria, which is based on the identification of orthologous operons. Since operons can move around in the genome but are rarely transformed (e.g. rarely broken into different parts), this approach allows for a better inference of orthologous genes in genomes that have been affected by many rearrangements, which in turn helps with the inference of more realistic evolutionary scenarios and ancestors.

Conclusions

From our comparisons of BOPAL with other gene order alignment programs using simulated data, we have found that BOPAL infers evolutionary events and ancestral gene orders more accurately than other methods based on alignments. An analysis of 12 Bacillus genomes also showed that BOPAL performs just as well as other programs at building ancestral histories in a minimal amount of events.

Genome replication dynamics of a bacteriophage and its satellite reveal strategies for parasitism and viral restriction

Phage-inducible chromosomal island-like elements (PLEs) are bacteriophage satellites found in Vibrio cholerae. PLEs parasitize the lytic phage ICP1, excising from the bacterial chromosome, replicating, and mobilizing to new host cells following cell lysis. PLEs protect their host cell populations by completely restricting the production of ICP1 progeny. Previously, it was found that ICP1 replication was reduced during PLE(+) infection. Despite robust replication of the PLE genome, relatively few transducing units are produced. We investigated if PLE DNA replication itself is antagonistic to ICP1 replication. Here we identify key constituents of PLE replication and assess their role in interference of ICP1. PLE encodes a RepA_N initiation factor that is sufficient to drive replication from the PLE origin of replication during ICP1 infection. In contrast to previously characterized bacteriophage satellites, expression of the PLE initiation factor was not sufficient for PLE replication in the absence of phage. Replication of PLE was necessary for interference of ICP1 DNA replication, but replication of a minimalized PLE replicon was not sufficient for ICP1 DNA replication interference. Despite restoration of ICP1 DNA replication, non-replicating PLE remained broadly inhibitory against ICP1. These results suggest that PLE DNA replication is one of multiple mechanisms contributing to ICP1 restriction.

Significance and roles of synonymous codon usage in the evolutionary process of Proteus

Proteus spp. bacteria frequently serve as opportunistic pathogens that can infect many animals and show positive survival and existence in various natural environments. The evolutionary pattern of Proteus spp. is an unknown topic, which benefits understanding the different evolutionary dynamics for excellent bacterial adaptation to various environments. Here, the eight whole genomes of different Proteus species were analyzed for the interplay between nucleotide usage and synonymous codon usage. Although the orthologous average nucleotide identity and average nucleotide identity display the genetic diversity of these Proteus species at the genome level, the principal component analysis further shows that these species sustain the specific genetic niche at the aspect of synonymous codon usage patterns. Interestingly, although these Proteus species have A/T rich genes with underrepresented G (guanine) or C (cytosine) at the third codon positions and overrepresented A or T at these positions, some synonymous codons with A or T end are obviously suppressed in usage. The overall codon usage pattern reflected by the effective number of codons (ENC) has a significantly positive correlation with GC3 content (GC content at the third codon position), and ENC has a significantly negative correlation with the adaptation index for these species. These results suggest that the mutation pressure caused by nucleotide composition constraint serves as a dominant evolutionary dynamic driving evolutionary trend of Proteus spp., along with other selections related to natural selection, replication and fine-tune translation, and so on. Taken together, the analyses help to understand the evolutionary interplay between nucleotide and codon usage at the gene level of Proteus.

Laboratory Evolution Experiments Help Identify a Predominant Region of Constitutive Stable DNA Replication Initiation

The bacterium Escherichia coli can initiate replication in the absence of the replication initiator protein DnaA and/or the canonical origin of replication oriC in a ΔrnhA background. This phenomenon, which can be primed by R-loops, is called constitutive stable DNA replication (cSDR). Whether DNA replication during cSDR initiates in a stochastic manner through the length of the chromosome or at specific sites and how E. coli can find adaptations to loss of fitness caused by cSDR remain inadequately answered. We use laboratory evolution experiments of ΔrnhA-ΔdnaA strains followed by deep sequencing to show that DNA replication preferentially initiates within a broad region located ∼0.4 to 0.7 Mb clockwise of oriC. This region includes many bisulfite-sensitive sites, which have been previously defined as R-loop-forming regions, and includes a site containing sequence motifs that favor R-loop formation. Initiation from this region would result in head-on replication-transcription conflicts at rRNA loci. Inversions of these rRNA loci, which can partly resolve these conflicts, help the bacterium suppress the fitness defects of cSDR. These inversions partially restore the gene expression changes brought about by cSDR. The inversion, however, increases the possibility of conflicts at essential mRNA genes, which would utilize only a minuscule fraction of RNA polymerase molecules, most of which transcribe rRNA genes. Whether subsequent adaptive strategies would attempt to resolve these conflicts remains an open question.IMPORTANCE The bacterium E. coli can replicate its DNA even in the absence of the molecules that are required for canonical replication initiation. This often requires the formation of RNA-DNA hybrid structures and is referred to as constitutive stable DNA replication (cSDR). Where on the chromosome does cSDR initiate? We answer this question using laboratory evolution experiments and genomics and show that selection favors cSDR initiation predominantly at a region ∼0.6 Mb clockwise of oriC. Initiation from this site will result in more head-on collisions of DNA polymerase with RNA polymerase operating on rRNA loci. The bacterium adapts to this problem by inverting a region of the genome including several rRNA loci such that head-on collisions between the two polymerases are minimized. Understanding such evolutionary strategies in the context of cSDR can provide insights into the potential causes of resistance to antibiotics that target initiation of DNA replication.

Mutation and Recombination Rates Vary Across Bacterial Chromosome

Bacteria evolve as a result of mutations and acquisition of foreign DNA by recombination processes. A growing body of evidence suggests that mutation and recombination rates are not constant across the bacterial chromosome. Bacterial chromosomal DNA is organized into a compact nucleoid structure which is established by binding of the nucleoid-associated proteins (NAPs) and other proteins. This review gives an overview of recent findings indicating that the mutagenic and recombination processes in bacteria vary at different chromosomal positions. Involvement of NAPs and other possible mechanisms in these regional differences are discussed. Variations in mutation and recombination rates across the bacterial chromosome may have implications in the evolution of bacteria.

Symmetric activity of DNA polymerases at and recruitment of exonuclease ExoR and of PolA to the Bacillus subtilis replication forks

DNA replication forks are intrinsically asymmetric and may arrest during the cell cycle upon encountering modifications in the DNA. We have studied real time dynamics of three DNA polymerases and an exonuclease at a single molecule level in the bacterium Bacillus subtilis. PolC and DnaE work in a symmetric manner and show similar dwell times. After addition of DNA damage, their static fractions and dwell times decreased, in agreement with increased re-establishment of replication forks. Only a minor fraction of replication forks showed a loss of active polymerases, indicating relatively robust activity during DNA repair. Conversely, PolA, homolog of polymerase I and exonuclease ExoR were rarely present at forks during unperturbed replication but were recruited to replications forks after induction of DNA damage. Protein dynamics of PolA or ExoR were altered in the absence of each other during exponential growth and during DNA repair, indicating overlapping functions. Purified ExoR displayed exonuclease activity and preferentially bound to DNA having 5′ overhangs in vitro. Our analyses support the idea that two replicative DNA polymerases work together at the lagging strand whilst only PolC acts at the leading strand, and that PolA and ExoR perform inducible functions at replication forks during DNA repair.

Biological Impact of a Large-Scale Genomic Inversion That Grossly Disrupts the Relative Positions of the Origin and Terminus Loci of the Streptococcus pyogenes Chromosome

A large-scale genomic inversion encompassing 0.79 Mb of the 1.816-Mb-long Streptococcus pyogenes serotype M49 strain NZ131 chromosome spontaneously occurs in a minor subpopulation of cells, and in this report genetic selection was used to obtain a stable lineage with this chromosomal rearrangement. This inversion, which drastically displaces the ori site relative to the terminus, changes the relative length of the replication arms so that one replichore is approximately 0.41 Mb while the other is about 1.40 Mb in length. Genomic reversion to the original chromosome constellation is not observed in PCR-monitored analyses after 180 generations of growth in rich medium. Compared to the parental strain, the inversion surprisingly demonstrates a nearly identical growth pattern in the first phase of the exponential phase, but differences do occur when resources in the medium become limited. When cultured separately in rich medium during prolonged stationary phase or in an experimental acute infection animal model (Galleria mellonella), the parental strain and the invertant have equivalent survival rates. However, when they are coincubated together, both in vitro and in vivo, the survival of the invertant declines relative to the level for the parental strain. The accompanying aspect of the study suggests that inversions taking place near oriC always happen to secure the linkage of oriC to DNA sequences responsible for chromosome partition. The biological relevance of large-scale inversions is also discussed.IMPORTANCE Based on our previous work, we created to our knowledge the largest asymmetric inversion, covering 43.5% of the S. pyogenes genome. In spite of a drastic replacement of origin of replication and the unbalanced size of replichores (1.4 Mb versus 0.41 Mb), the invertant, when not challenged with its progenitor, showed impressive vitality for growth in vitro and in pathogenesis assays. The mutant supports the existing idea that slightly deleterious mutations can provide the setting for secondary adaptive changes. Furthermore, comparative analysis of the mutant with previously published data strongly indicates that even large genomic rearrangements survive provided that the integrity of the oriC and the chromosome partition cluster is preserved.

Clustered Core- and Pan-Genome Content on Rhodobacteraceae Chromosomes

In Bacteria, chromosome replication starts at a single origin of replication and proceeds on both replichores. Due to its asymmetric nature, replication influences chromosome structure and gene organization, mutation rate, and expression. To date, little is known about the distribution of highly conserved genes over the bacterial chromosome. Here, we used a set of 101 fully sequenced Rhodobacteraceae representatives to analyze the relationship between conservation of genes within this family and their distance from the origin of replication. Twenty-two of the analyzed species had core genes clustered significantly closer to the origin of replication with representatives of the genus Celeribacter being the most apparent example. Interestingly, there were also eight species with the opposite organization. In particular, Rhodobaca barguzinensis and Loktanella vestfoldensis showed a significant increase of core genes with distance from the origin of replication. The uneven distribution of low-conserved regions is in particular pronounced for genomes in which the halves of one replichore differ in their conserved gene content. Phage integration and horizontal gene transfer partially explain the scattered nature of Rhodobacteraceae genomes. Our findings lay the foundation for a better understanding of bacterial genome evolution and the role of replication therein.

Design of a microaerobically inducible replicon for high-yield plasmid DNA production

A pUC-derived replicon inducible by oxygen limitation was designed and tested in fed-batch cultures of Escherichia coli. It included the addition of a second inducible copy of rnaII, the positive replication control element. The rnaII gene was expressed from P_trc and cloned into pUC18 to test the hypothesis that the ratio of the positive control molecule RNAII to the negative control element, RNAI, was the determinant of plasmid copy number per chromosome (PCN). The construct was evaluated in several E. coli strains. Evaluations of the RNAII/RNAI ratio, PCN and plasmid yield normalized to biomass (Y_pDNA/X ) were performed and the initial hypothesis was probed. Furthermore, in high cell-density cultures in shake flasks, an outstanding amount of 126 mg/L of plasmid was produced. The microaerobically inducible plasmid was obtained by cloning the rnaII gene under the control of the oxygen-responsive Vitreoscilla stercoraria hemoglobin promoter. For this plasmid, but not for pUC18, the RNAII/RNAI ratio, PCN and Y_pDNA/X efficiently increased after the shift to the microaerobic regime in fed-batch cultures in a 1 L bioreactor. The Y_pDNA/X of the inducible plasmid reached 12 mg/g at the end of the fed-batch but the original pUC18 only reached ca. 6 mg/g. The proposed plasmid is a valuable alternative for the operation and scale-up of plasmid DNA production processes in which mass transfer limitations will not represent an issue.

The mechanism of DNA unwinding by the eukaryotic replicative helicase

Accurate DNA replication is tightly regulated in eukaryotes to ensure genome stability during cell division and is performed by the multi-protein replisome. At the core an AAA+ hetero-hexameric complex, Mcm2-7, together with GINS and Cdc45 form the active replicative helicase Cdc45/Mcm2-7/GINS (CMG). It is not clear how this replicative ring helicase translocates on, and unwinds, DNA. We measure real-time dynamics of purified recombinant Drosophila melanogaster CMG unwinding DNA with single-molecule magnetic tweezers. Our data demonstrates that CMG exhibits a biased random walk, not the expected unidirectional motion. Through building a kinetic model we find CMG may enter up to three paused states rather than unwinding, and should these be prevented, in vivo fork rates would be recovered in vitro. We propose a mechanism in which CMG couples ATP hydrolysis to unwinding by acting as a lazy Brownian ratchet, thus providing quantitative understanding of the central process in eukaryotic DNA replication.

Pseudomonas putida rDNA is a favored site for the expression of biosynthetic genes

Since high-value bacterial secondary metabolites, including antibiotics, are often naturally produced in only low amounts, their efficient biosynthesis typically requires the transfer of entire metabolic pathways into suitable bacterial hosts like Pseudomonas putida. Stable maintenance and sufficient expression of heterologous pathway-encoding genes in host microbes, however, still remain key challenges. In this study, the 21 kb prodigiosin gene cluster from Serratia marcescens was used as a reporter to identify genomic sites in P. putida KT2440 especially suitable for maintenance and expression of pathway genes. After generation of a strain library by random Tn5 transposon-based chromosomal integration of the cluster, 50 strains exhibited strong prodigiosin production. Remarkably, chromosomal integration sites were exclusively identified in the seven rRNA-encoding rrn operons of P. putida. We could further demonstrate that prodigiosin production was mainly dependent on (i) the individual rrn operon where the gene cluster was inserted as well as (ii) the distance between the rrn promoter and the inserted prodigiosin biosynthetic genes. In addition, the recombinant strains showed high stability upon subculturing for many generations. Consequently, our findings demonstrate the general applicability of rDNA loci as chromosomal integration sites for gene cluster expression and recombinant pathway implementation in P. putida KT2440.

Quantitative analysis and assessment of base composition asymmetry and gene orientation bias in bacterial genomes

Base composition asymmetry and gene orientation bias are two common genomic structures in bacterial genomes. Here, correlation coefficients between nucleotide disparities and coding sequence (CDS) skew have been calculated, which provides insights into their relationship from an individual genome perspective. Consequently, we find GC and RY disparities correlate significantly with CDS skew, since around 60% of the bacterial genomes under study have correlation coefficients > 0.9. Then, we present a model for quantitative assessment of nucleotide disparity and CDS skew in which a numerical index R² is used for evaluation. We find that skew curves with higher R² perform better on the prediction of replication origins in bacteria.

Carbapenemases on the move: it's good to be on ICEs

Background

The evolution and spread of antibiotic resistance is often mediated by mobile genetic elements. Integrative and conjugative elements (ICEs) are the most abundant conjugative elements among prokaryotes. However, the contribution of ICEs to horizontal gene transfer of antibiotic resistance has been largely unexplored.

Results

Here we report that ICEs belonging to mating-pair formation (MPF) classes G and T are highly prevalent among the opportunistic pathogen Pseudomonas aeruginosa, contributing to the spread of carbapenemase-encoding genes (CEGs). Most CEGs of the MPF_G class were encoded within class I integrons, which co-harbour genes conferring resistance to other antibiotics. The majority of the integrons were located within Tn3-like and composite transposons. Conserved attachment site could be predicted for the MPF_G class ICEs. MPF_T class ICEs carried the CEGs within composite transposons which were not associated with integrons.

Conclusions

The data presented here provides a global snapshot of the different CEG-harbouring ICEs and sheds light on the underappreciated contribution of these elements to the evolution and dissemination of antibiotic resistance on P. aeruginosa.

Electronic supplementary material

The online version of this article (10.1186/s13100-018-0141-4) contains supplementary material, which is available to authorized users.

How Acts of Infidelity Promote DNA Break Repair: Collision and Collusion Between DNA Repair and Transcription

Transcription is a fundamental cellular process and the first step in gene regulation. Although RNA polymerase (RNAP) is highly processive, in growing cells the progression of transcription can be hindered by obstacles on the DNA template, such as damaged DNA. The authors recent findings highlight a trade-off between transcription fidelity and DNA break repair. While a lot of work has focused on the interaction between transcription and nucleotide excision repair, less is known about how transcription influences the repair of DNA breaks. The authors suggest that when the cell experiences stress from DNA breaks, the control of RNAP processivity affects the balance between preserving transcription integrity and DNA repair. Here, how the conflict between transcription and DNA double-strand break (DSB) repair threatens the integrity of both RNA and DNA are discussed. In reviewing this field, the authors speculate on cellular paradigms where this equilibrium is well sustained, and instances where the maintenance of transcription fidelity is favored over genome stability.

Chromosome architecture constrains horizontal gene transfer in bacteria

Despite significant frequencies of lateral gene transfer between species, higher taxonomic groups of bacteria show ecological and phenotypic cohesion. This suggests that barriers prevent panmictic dissemination of genes via lateral gene transfer. We have proposed that most bacterial genomes have a functional architecture imposed by Architecture IMparting Sequences (AIMS). AIMS are defined as 8 base pair sequences preferentially abundant on leading strands, whose abundance and strand-bias are positively correlated with proximity to the replication terminus. We determined that inversions whose endpoints lie within a single chromosome arm, which would reverse the polarity of AIMS in the inverted region, are both shorter and less frequent near the replication terminus. This distribution is consistent with the increased selection on AIMS function in this region, thus constraining DNA rearrangement. To test the hypothesis that AIMS also constrain DNA transfer between genomes, AIMS were identified in genomes while ignoring atypical, potentially laterally-transferred genes. The strand-bias of AIMS within recently acquired genes was negatively correlated with the distance of those genes from their genome’s replication terminus. This suggests that selection for AIMS function prevents the acquisition of genes whose AIMS are not found predominantly in the permissive orientation. This constraint has led to the loss of at least 18% of genes acquired by transfer in the terminus-proximal region. We used completely sequenced genomes to produce a predictive road map of paths of expected horizontal gene transfer between species based on AIMS compatibility between donor and recipient genomes. These results support a model whereby organisms retain introgressed genes only if the benefits conferred by their encoded functions outweigh the detriments incurred by the presence of foreign DNA lacking genome-wide architectural information.

The potential success of horizontal gene transfer events is historically equated to the benefits conferred by encoded products. Here we show that gene transfer events are observed less frequently if the introduced genes disrupt important patterns of genomic information, suggesting that this disruption would confer an unacceptable cost. As a result, gene transfer events are less likely to be successful if the potential donor genomes have incompatible genome architecture. Because more distantly-related genes are less compatible, chromosome architecture serves as a mechanism to bias gene transfer events to those involving closer relatives, thereby providing a mechanism for the genotypic and phenotypic cohesion of higher taxonomic groups.