Comparative Genomics of Interreplichore Translocations in Bacteria: A Measure of Chromosome Topology?

Replication-associated inversions are the dominant form of bacterial chromosome structural variation

The structural arrangements of bacterial chromosomes vary widely between closely related species and can result in significant phenotypic outcomes. The appearance of large-scale chromosomal inversions that are symmetric relative to markers for the origin of replication (OriC) has been previously observed; however, the overall prevalence of replication-associated structural rearrangements (RASRs) in bacteria and their causal mechanisms are currently unknown. Here, we systematically identify the locations of RASRs in species with multiple complete-sequenced genomes and investigate potential mediating biological mechanisms. We found that 247 of 313 species contained sequences with at least one large (>50 Kb) inversion in their sequence comparisons, and the aggregated inversion distances away from symmetry were normally distributed with a mean of zero. Many inversions that were offset from dnaA were found to be centered on a different marker for the OriC Instances of flanking repeats provide evidence that breaks formed during the replication process could be repaired to opposing positions. We also found a strong relationship between the later stages of replication and the range in distance variation from symmetry.

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.

Myxobacterial Genomics and Post-Genomics: A Review of Genome Biology, Genome Sequences and Related 'Omics Studies

Myxobacteria are fascinating and complex microbes. They prey upon other members of the soil microbiome by secreting antimicrobial proteins and metabolites, and will undergo multicellular development if starved. The genome sequence of the model myxobacterium Myxococcus xanthus DK1622 was published in 2006 and 15 years later, 163 myxobacterial genome sequences have now been made public. This explosion in genomic data has enabled comparative genomics analyses to be performed across the taxon, providing important insights into myxobacterial gene conservation and evolution. The availability of myxobacterial genome sequences has allowed system-wide functional genomic investigations into entire classes of genes. It has also enabled post-genomic technologies to be applied to myxobacteria, including transcriptome analyses (microarrays and RNA-seq), proteome studies (gel-based and gel-free), investigations into protein–DNA interactions (ChIP-seq) and metabolism. Here, we review myxobacterial genome sequencing, and summarise the insights into myxobacterial biology that have emerged as a result. We also outline the application of functional genomics and post-genomic approaches in myxobacterial research, highlighting important findings to emerge from seminal studies. The review also provides a comprehensive guide to the genomic datasets available in mid-2021 for myxobacteria (including 24 genomes that we have sequenced and which are described here for the first time).

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.

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.

Non-Random Inversion Landscapes in Prokaryotic Genomes Are Shaped by Heterogeneous Selection Pressures

Inversions are a major contributor to structural genome evolution in prokaryotes. Here, using a novel alignment-based method, we systematically compare 1,651 bacterial and 98 archaeal genomes to show that inversion landscapes are frequently biased toward (symmetric) inversions around the origin–terminus axis. However, symmetric inversion bias is not a universal feature of prokaryotic genome evolution but varies considerably across clades. At the extremes, inversion landscapes in Bacillus–Clostridium and Actinobacteria are dominated by symmetric inversions, while there is little or no systematic bias favoring symmetric rearrangements in archaea with a single origin of replication. Within clades, we find strong but clade-specific relationships between symmetric inversion bias and different features of adaptive genome architecture, including the distance of essential genes to the origin of replication and the preferential localization of genes on the leading strand. We suggest that heterogeneous selection pressures have converged to produce similar patterns of structural genome evolution across prokaryotes.

The Transcriptome of Streptococcus pneumoniae Induced by Local and Global Changes in Supercoiling

The bacterial chromosome is compacted in a manner optimal for DNA transactions to occur. The degree of compaction results from the level of DNA-supercoiling and the presence of nucleoid-binding proteins. DNA-supercoiling is homeostatically maintained by the opposing activities of relaxing DNA topoisomerases and negative supercoil-inducing DNA gyrase. DNA-supercoiling acts as a general cis regulator of transcription, which can be superimposed upon other types of more specific trans regulatory mechanism. Transcriptomic studies on the human pathogen Streptococcus pneumoniae, which has a relatively small genome (∼2 Mb) and few nucleoid-binding proteins, have been performed under conditions of local and global changes in supercoiling. The response to local changes induced by fluoroquinolone antibiotics, which target DNA gyrase subunit A and/or topoisomerase IV, involves an increase in oxygen radicals which reduces cell viability, while the induction of global supercoiling changes by novobiocin (a DNA gyrase subunit B inhibitor), or by seconeolitsine (a topoisomerase I inhibitor), has revealed the existence of topological domains that specifically respond to such changes. The control of DNA-supercoiling in S. pneumoniae occurs mainly via the regulation of topoisomerase gene transcription: relaxation triggers the up-regulation of gyrase and the down-regulation of topoisomerases I and IV, while hypernegative supercoiling down-regulates the expression of topoisomerase I. Relaxation affects 13% of the genome, with the majority of the genes affected located in 15 domains. Hypernegative supercoiling affects 10% of the genome, with one quarter of the genes affected located in 12 domains. However, all the above domains overlap, suggesting that the chromosome is organized into topological domains with fixed locations. Based on its response to relaxation, the pneumococcal chromosome can be said to be organized into five types of domain: up-regulated, down-regulated, position-conserved non-regulated, position-variable non-regulated, and AT-rich. The AT content is higher in the up-regulated than in the down-regulated domains. Genes within the different domains share structural and functional characteristics. It would seem that a topology-driven selection pressure has defined the chromosomal location of the metabolism, virulence and competence genes, which suggests the existence of topological rules that aim to improve bacterial fitness.

Bioinformatic Analyses of Unique (Orphan) Core Genes of the Genus Acidithiobacillus: Functional Inferences and Use As Molecular Probes for Genomic and Metagenomic/Transcriptomic Interrogation

Using phylogenomic and gene compositional analyses, five highly conserved gene families have been detected in the core genome of the phylogenetically coherent genus Acidithiobacillus of the class Acidithiobacillia. These core gene families are absent in the closest extant genus Thermithiobacillus tepidarius that subtends the Acidithiobacillus genus and roots the deepest in this class. The predicted proteins encoded by these core gene families are not detected by a BLAST search in the NCBI non-redundant database of more than 90 million proteins using a relaxed cut-off of 1.0e⁻⁵. None of the five families has a clear functional prediction. However, bioinformatic scrutiny, using pI prediction, motif/domain searches, cellular location predictions, genomic context analyses, and chromosome topology studies together with previously published transcriptomic and proteomic data, suggests that some may have functions associated with membrane remodeling during cell division perhaps in response to pH stress. Despite the high level of amino acid sequence conservation within each family, there is sufficient nucleotide variation of the respective genes to permit the use of the DNA sequences to distinguish different species of Acidithiobacillus, making them useful additions to the armamentarium of tools for phylogenetic analysis. Since the protein families are unique to the Acidithiobacillus genus, they can also be leveraged as probes to detect the genus in environmental metagenomes and metatranscriptomes, including industrial biomining operations, and acid mine drainage (AMD).