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

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.

RNA polymerase collisions and their role in transcription

RNA polymerases are the central enzymes of gene expression and function frequently in either a head-on or co-directional manner on the busy DNA track. Whether and how these collisions between RNA polymerases contribute to transcriptional regulation is mysterious. Increasing evidence from biochemical and single-molecule studies suggests that RNA polymerase collisions function as an important regulator to fine-tune transcription, rather than creating deleterious "traffic jams". This review summarizes the recent progress on elucidating the consequences of RNA polymerase collisions during transcription and highlights the significance of cooperation and coordination between RNA polymerases.

Mutation bias and adaptation in bacteria

Genetic mutation, which provides the raw material for evolutionary adaptation, is largely a stochastic force. However, there is ample evidence showing that mutations can also exhibit strong biases, with some mutation types and certain genomic positions mutating more often than others. It is becoming increasingly clear that mutational bias can play a role in determining adaptive outcomes in bacteria in both the laboratory and the clinic. As such, understanding the causes and consequences of mutation bias can help microbiologists to anticipate and predict adaptive outcomes. In this review, we provide an overview of the mechanisms and features of the bacterial genome that cause mutational biases to occur. We then describe the environmental triggers that drive these mechanisms to be more potent and outline the adaptive scenarios where mutation bias can synergize with natural selection to define evolutionary outcomes. We conclude by describing how understanding mutagenic genomic features can help microbiologists predict areas sensitive to mutational bias, and finish by outlining future work that will help us achieve more accurate evolutionary forecasts.

Mutational signatures and increased retrotransposon insertions in xeroderma pigmentosum variant skin tumors

Xeroderma pigmentosum variant (XP-V) is an autosomal recessive disease with an increased risk of developing cutaneous neoplasms in sunlight-exposed regions. These cells are deficient in the translesion synthesis (TLS) DNA polymerase eta, responsible for bypassing different types of DNA lesions. From the exome sequencing of 11 skin tumors of a genetic XP-V patients' cluster, classical mutational signatures related to sunlight exposure, such as C>T transitions targeted to pyrimidine dimers, were identified. However, basal cell carcinomas also showed distinct C>A mutation spectra reflecting a mutational signature possibly related to sunlight-induced oxidative stress. Moreover, four samples carry different mutational signatures, with C>A mutations associated with tobacco chewing or smoking usage. Thus, XP-V patients should be warned of the risk of these habits. Surprisingly, higher levels of retrotransposon somatic insertions were also detected when the tumors were compared with non-XP skin tumors, revealing other possible causes for XP-V tumors and novel functions for the TLS polymerase eta in suppressing retrotransposition. Finally, the expected high mutation burden found in most of these tumors renders these XP patients good candidates for checkpoint blockade immunotherapy.

Selective Pressure by Rifampicin Modulates Mutation Rates and Evolutionary Trajectories of Mycobacterial Genomes

Resistance to the frontline antibiotic rifampicin constitutes a challenge to the treatment and control of tuberculosis. Here, we analyzed the mutational landscape of Mycobacterium smegmatis during long-term evolution with increasing concentrations of rifampicin, using a mutation accumulation assay combined with whole-genome sequencing. Antibiotic treatment enhanced the acquisition of mutations, doubling the genome-wide mutation rate of the wild-type cells. While antibiotic exposure led to extinction of almost all wild-type lines, the hypermutable phenotype of the ΔnucS mutant strain (noncanonical mismatch repair deficient) provided an efficient response to the antibiotic, leading to high rates of survival. This adaptative advantage resulted in the emergence of higher levels of rifampicin resistance, an accelerated acquisition of drug resistance mutations in rpoB (β RNA polymerase), and a wider diversity of evolutionary pathways that led to drug resistance. Finally, this approach revealed a subset of adaptive genes under positive selection with rifampicin that could be associated with the development of antibiotic resistance. IMPORTANCE Rifampicin is the most important first-line antibiotic against mycobacterial infections, including tuberculosis, one of the top causes of death worldwide. Acquisition of rifampicin resistance constitutes a major global public health problem that makes the control of the disease challenging. Here, we performed an experimental evolution assay under antibiotic selection to analyze the response and adaptation of mycobacteria, leading to the acquisition of rifampicin resistance. This approach explored the total number of mutations that arose in the mycobacterial genomes under long-term rifampicin exposure, using whole-genome sequencing. Our results revealed the effect of rifampicin at a genomic level, identifying different mechanisms and multiple pathways leading to rifampicin resistance in mycobacteria. Moreover, this study detected that an increase in the rate of mutations led to enhanced levels of drug resistance and survival. In summary, all of these results could be useful to understand and prevent the emergence of drug-resistant isolates in mycobacterial infections.

RNase H genes cause distinct impacts on RNA:DNA hybrid formation and mutagenesis genome wide

RNA:DNA hybrids compromise replication fork progression and genome integrity in all cells. The overall impacts of naturally occurring RNA:DNA hybrids on genome integrity, and the relative contributions of ribonucleases H to mitigating the negative effects of hybrids, remain unknown. Here, we investigate the contributions of RNases HII (RnhB) and HIII (RnhC) to hybrid removal, DNA replication, and mutagenesis genome wide. Deletion of either rnhB or rnhC triggers RNA:DNA hybrid accumulation but with distinct patterns of mutagenesis and hybrid accumulation. Across all cells, hybrids accumulate strongly in noncoding RNAs and 5'-UTRs of coding sequences. For ΔrnhB, hybrids accumulate preferentially in untranslated regions and early in coding sequences. We show that hybrid accumulation is particularly sensitive to gene expression in ΔrnhC cells. DNA replication in ΔrnhC cells is disrupted, leading to transversions and structural variation. Our results resolve the outstanding question of how hybrids in native genomic contexts cause mutagenesis and shape genome organization.

Interplay between chromosomal architecture and termination of DNA replication in bacteria

Faithful transmission of the genome from one generation to the next is key to life in all cellular organisms. In the majority of bacteria, the genome is comprised of a single circular chromosome that is normally replicated from a single origin, though additional genetic information may be encoded within much smaller extrachromosomal elements called plasmids. By contrast, the genome of a eukaryote is distributed across multiple linear chromosomes, each of which is replicated from multiple origins. The genomes of archaeal species are circular, but are predominantly replicated from multiple origins. In all three cases, replication is bidirectional and terminates when converging replication fork complexes merge and 'fuse' as replication of the chromosomal DNA is completed. While the mechanics of replication initiation are quite well understood, exactly what happens during termination is far from clear, although studies in bacterial and eukaryotic models over recent years have started to provide some insight. Bacterial models with a circular chromosome and a single bidirectional origin offer the distinct advantage that there is normally just one fusion event between two replication fork complexes as synthesis terminates. Moreover, whereas termination of replication appears to happen in many bacteria wherever forks happen to meet, termination in some bacterial species, including the well-studied bacteria Escherichia coli and Bacillus subtilis, is more restrictive and confined to a 'replication fork trap' region, making termination even more tractable. This region is defined by multiple genomic terminator (ter) sites, which, if bound by specific terminator proteins, form unidirectional fork barriers. In this review we discuss a range of experimental results highlighting how the fork fusion process can trigger significant pathologies that interfere with the successful conclusion of DNA replication, how these pathologies might be resolved in bacteria without a fork trap system and how the acquisition of a fork trap might have provided an alternative and cleaner solution, thus explaining why in bacterial species that have acquired a fork trap system, this system is remarkably well maintained. Finally, we consider how eukaryotic cells can cope with a much-increased number of termination events.

RNase H genes cause distinct impacts on RNA:DNA hybrid formation and mutagenesis genome-wide

RNA:DNA hybrids such as R-loops affect genome integrity and DNA replication fork progression. The overall impacts of naturally occurring RNA:DNA hybrids on genome integrity, and the relative contributions of ribonucleases H to mitigating the negative effects of hybrids, remain unknown. Here, we investigate the contributions of RNases HII (RnhB) and HIII (RnhC) to hybrid removal, DNA replication, and mutagenesis genome-wide. Deletion of either rnhB or rnhC triggers RNA:DNA hybrid accumulation, but with distinct patterns of mutagenesis and hybrid accumulation. Across all cells, hybrids accumulate most strongly in non-coding RNAs and 5'-UTRs of coding sequences. For Δ rnhB , hybrids accumulate preferentially in untranslated regions and early in coding sequences. Hybrid accumulation is particularly sensitive to gene expression in Δ rnhC ; in cells lacking RnhC, DNA replication is disrupted leading to transversions and structural variation. Our results resolve the outstanding question of how hybrids in native genomic contexts interact with replication to cause mutagenesis and shape genome organization.

Dynamic state of plasmid genomic architectures resulting from XerC/D-mediated site-specific recombination in Acinetobacter baumannii Rep_3 superfamily resistance plasmids carrying blaOXA-58 - and TnaphA6-resistance modules

The acquisition of bla _OXA genes encoding different carbapenem-hydrolyzing class-D β-lactamases (CHDL) represents a main determinant of carbapenem resistance in the nosocomial pathogen Acinetobacter baumannii. The bla_OXA-58 gene, in particular, is generally embedded in similar resistance modules (RM) carried by plasmids unique to the Acinetobacter genus lacking self-transferability. The ample variations in the immediate genomic contexts in which bla_OXA-58 -containing RMs are inserted among these plasmids, and the almost invariable presence at their borders of non-identical 28-bp sequences potentially recognized by the host XerC and XerD tyrosine recombinases (pXerC/D-like sites), suggested an involvement of these sites in the lateral mobilization of the gene structures they encircle. However, whether and how these pXerC/D sites participate in this process is only beginning to be understood. Here, we used a series of experimental approaches to analyze the contribution of pXerC/D-mediated site-specific recombination to the generation of structural diversity between resistance plasmids carrying pXerC/D-bounded bla _OXA-58- and TnaphA6-containing RM harbored by two phylogenetically- and epidemiologically-closely related A. baumannii strains of our collection, Ab242 and Ab825, during adaptation to the hospital environment. Our analysis disclosed the existence of different bona fide pairs of recombinationally-active pXerC/D sites in these plasmids, some mediating reversible intramolecular inversions and others reversible plasmid fusions/resolutions. All of the identified recombinationally-active pairs shared identical GGTGTA sequences at the cr spacer separating the XerC- and XerD-binding regions. The fusion of two Ab825 plasmids mediated by a pair of recombinationally-active pXerC/D sites displaying sequence differences at the cr spacer could be inferred on the basis of sequence comparison analysis, but no evidence of reversibility could be obtained in this case. The reversible plasmid genome rearrangements mediated by recombinationally-active pairs of pXerC/D sites reported here probably represents an ancient mechanism of generating structural diversity in the Acinetobacter plasmid pool. This recursive process could facilitate a rapid adaptation of an eventual bacterial host to changing environments, and has certainly contributed to the evolution of Acinetobacter plasmids and the capture and dissemination of bla _OXA-58 genes among Acinetobacter and non-Acinetobacter populations co-residing in the hospital niche.

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.

OUP accepted manuscript

Conflicts between transcription and replication machinery are a potent source of replication stress and genome instability; however, no technique currently exists to identify endogenous genomic locations prone to transcription–replication interactions. Here, we report a novel method to identify genomic loci prone to transcription–replication interactions termed transcription–replication immunoprecipitation on nascent DNA sequencing, TRIPn-Seq. TRIPn-Seq employs the sequential immunoprecipitation of RNA polymerase 2 phosphorylated at serine 5 (RNAP2s5) followed by enrichment of nascent DNA previously labeled with bromodeoxyuridine. Using TRIPn-Seq, we mapped 1009 unique transcription–replication interactions (TRIs) in mouse primary B cells characterized by a bimodal pattern of RNAP2s5, bidirectional transcription, an enrichment of RNA:DNA hybrids, and a high probability of forming G-quadruplexes. TRIs are highly enriched at transcription start sites and map to early replicating regions. TRIs exhibit enhanced Replication Protein A association and TRI-associated genes exhibit higher replication fork termination than control transcription start sites, two marks of replication stress. TRIs colocalize with double-strand DNA breaks, are enriched for deletions, and accumulate mutations in tumors. We propose that replication stress at TRIs induces mutations potentially contributing to age-related disease, as well as tumor formation and development.

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.

Methyltransferase-directed orthogonal tagging and sequencing of miRNAs and bacterial small RNAs

Background

Targeted installation of designer chemical moieties on biopolymers provides an orthogonal means for their visualisation, manipulation and sequence analysis. Although high-throughput RNA sequencing is a widely used method for transcriptome analysis, certain steps, such as 3′ adapter ligation in strand-specific RNA sequencing, remain challenging due to structure- and sequence-related biases introduced by RNA ligases, leading to misrepresentation of particular RNA species. Here, we remedy this limitation by adapting two RNA 2′-O-methyltransferases from the Hen1 family for orthogonal chemo-enzymatic click tethering of a 3′ sequencing adapter that supports cDNA production by reverse transcription of the tagged RNA.

Results

We showed that the ssRNA-specific DmHen1 and dsRNA-specific AtHEN1 can be used to efficiently append an oligonucleotide adapter to the 3′ end of target RNA for sequencing library preparation. Using this new chemo-enzymatic approach, we identified miRNAs and prokaryotic small non-coding sRNAs in probiotic Lactobacillus casei BL23. We found that compared to a reference conventional RNA library preparation, methyltransferase-Directed Orthogonal Tagging and RNA sequencing, mDOT-seq, avoids misdetection of unspecific highly-structured RNA species, thus providing better accuracy in identifying the groups of transcripts analysed. Our results suggest that mDOT-seq has the potential to advance analysis of eukaryotic and prokaryotic ssRNAs.

Conclusions

Our findings provide a valuable resource for studies of the RNA-centred regulatory networks in Lactobacilli and pave the way to developing novel transcriptome and epitranscriptome profiling approaches in vitro and inside living cells. As RNA methyltransferases share the structure of the AdoMet-binding domain and several specific cofactor binding features, the basic principles of our approach could be easily translated to other AdoMet-dependent enzymes for the development of modification-specific RNA-seq techniques.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01053-w.

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.