Chromosome architecture constrains horizontal gene transfer in bacteria

Chromosome architecture as a determinant for biosynthetic diversity in Micromonospora

Natural products - small molecules generated by organisms to facilitate ecological interactions - are of great importance to society and are used as antibacterial, antiviral, antifungal and anticancer drugs. However, the role and evolution of these molecules and the fitness benefits they provide to their hosts in their natural habitat remain an outstanding question. In bacteria, the genes that encode the biosynthetic proteins that generate these molecules are organised into discrete loci termed biosynthetic gene clusters (BGCs). In this work, we asked the following question: How are biosynthetic gene clusters organised at the chromosomal level? We sought to answer this using publicly available high-quality assemblies of Micromonospora, an actinomycete genus with members responsible for biosynthesizing notable natural products, such as gentamicin and calicheamicin. By orienting the Micromonospora chromosome around the origin of replication, we demonstrated that Micromonospora has a conserved origin-proximal region, which becomes progressively more disordered towards the antipodes of the origin. We then demonstrated through genome mining of these organisms that the conserved origin-proximal region and the origin-distal region of Micromonospora have distinct populations of BGCs and, in this regard, parallel the organization of Streptomyces, which possesses linear chromosomes. Specifically, the origin-proximal region contains highly syntenous, conserved BGCs predicted to biosynthesize terpenes and a type III polyketide synthase. In contrast, the ori-distal region contains a highly diverse population of BGCs, with many BGCs belonging to unique gene cluster families. These data highlight that genomic plasticity in Micromonospora is locus-specific, and highlight the importance of using high-quality genome assemblies for natural product discovery and guide future natural product discovery by highlighting that biosynthetic novelty may be enriched in specific chromosomal neighbourhoods.

Characteristics of tetracycline antibiotic resistance gene enrichment and migration in soil-plant system

Tetracycline Resistance Genes (TRGs) have received widespread attention in recent years, as they are a novel environmental pollutant that can rapidly accumulate and migrate in soil plant systems through horizontal gene transfer (HGT), posing a potential threat to food safety and public health. This article systematically reviews the pollution sources, enrichment, and migration characteristics of TRGs in soil. The main sources of TRGs include livestock manure and contaminated wastewater, especially in intensive farming environments where TRGs pollution is more severe. In soil, TRGs diffuse horizontally between bacteria and migrate to plant tissues through mechanisms such as plasmid conjugation, integron mediation, and phage transduction. The migration of TRGs is not limited to the soil interior, and increasing evidence suggests that they can also enter the plant system through plant root absorption and the HGT pathway of endophytic bacteria, ultimately accumulating in plant roots, stems, leaves, fruits, and other parts. This process has a direct impact on human health, especially when TRGs are found in crops such as vegetables, which may be transmitted to the human body through the food chain. In addition, this article also deeply analyzed various factors that affect the migration of TRGs, including the residual level of tetracycline in soil, the type and concentration of microorganisms, heavy metal pollution, and the presence of new pollutants such as microplastics. These factors significantly affect the enrichment rate and migration mode of TRGs in soil. In addition, two technologies that can effectively eliminate TRGs in livestock breeding environments were introduced, providing reference for healthy agricultural production. The article concludes by summarizing the shortcomings of current research on TRGs, particularly the limited understanding of TRG migration pathways and their impact mechanisms. Future research should focus on revealing the migration mechanisms of TRGs in soil plant systems and developing effective control and governance measures to reduce the environmental transmission risks of TRGs and ensure the safety of ecosystems and human health.

Data-driven identification of major axes of functional variation in bacteria

The discovery of major axes of correlated functional variation among species and habitats has revealed the fundamental trade-offs structuring both functional and taxonomic diversity in eukaryotes such as plants. Whether such functional axes exist in the bacterial realm and whether they could explain bacterial taxonomic turnover among ecosystems remains unknown. Here, we use a data-driven approach to leverage global genomic and metagenomic datasets to reveal the existence of major axes of functional variation explaining both evolutionary differentiation within Bacteria and their ecological sorting across diverse habitats. We show that metagenomic variation among bacterial communities from various ecosystems is structured along a few axes of correlated functional pathways. Similar clusters of traits explained phylogenetic trait variation among >16,000 bacterial genomes, suggesting that functional turnover among bacterial communities from distinct habitats does not only result from the differential filtering of similar functions among communities, but also from phylogenetic correlations among these functions. Concordantly, functional pathways associated with trait clusters that were most important for defining functional turnover among bacterial communities were also those that had the highest phylogenetic signal in the bacterial genomic phylogeny. This study overall underlines the important role of evolutionary history in shaping contemporary distributions of bacteria across ecosystems.

Natural Transformation in Gram-Positive Bacteria and Its Biotechnological Relevance to Lactic Acid Bacteria

Competence refers to the specialized physiological state in which bacteria undergo transformation through the internalization of exogenous DNA in a controlled and genetically encoded process that leads to genotypic and, in many cases, phenotypic changes. Natural transformation was first described in Streptococcus pneumoniae and has since been demonstrated in numerous species, including Bacillus subtilis and Neisseria gonorrhoeae. Homologs of the genes encoding the DNA uptake machinery for natural transformation have been reported to be present in several lactic acid bacteria, including Lactobacillus spp., Streptococcus thermophilus, and Lactococcus spp. In this review, we collate current knowledge of the phenomenon of natural transformation in Gram-positive bacteria. Furthermore, we describe the mechanism of competence development and its regulation in model bacterial species. We highlight the importance and opportunities for the application of these findings in the context of bacterial starter cultures associated with food fermentations as well as current limitations in this area of research.

Genomic inversions in Escherichia coli alter gene expression and are associated with nucleoid protein binding sites

Genomic reorganization, like rearrangements and inversions, influence how genetic information is organized within bacterial genomes. Inversions in particular, facilitate genome evolution through gene gain and loss, and can alter gene expression. Previous studies investigating the impact inversions have on gene expression induced inversions targeting specific genes or examine inversions between distantly related species. This fails to encompass a genome wide perspective on naturally occurring inversions and their post adaptation impact on gene expression. Here we use bioinformatic techniques and multiple RNA-seq datasets to investigate the short- and long-range impact inversions have on genomic gene expression within <i>Escherichia coli</i>. We observed differences in gene expression between homologous inverted and non-inverted genes, even after long term exposure to adaptive selection. In 4% of inversions representing 33 genes, differential gene expression between inverted and non-inverted homologs was detected, with nearly two thirds (71%) of differentially expressed inverted genes having 9.4-85.6 fold higher gene expression. The identified inversions had more overlap than expected with nucleoid associated protein binding sites, which assist in genomic gene expression regulation. Some inversions can drastically impact gene expression even between different strains of <i>E.coli</i>, and could provide a mechanism for the diversification of genetic content through controlled expression changes.

Nucleoid-associated Rok differentially affects chromosomal transformation on Bacillus subtilis recombination-deficient cells

Rok, a Bacillus subtilis nucleoid-associated protein (NAP), negatively regulates competence development and silences xenogeneic genes. We show that rok inactivation increases rpoB482 natural intraspecies chromosomal transformation (CT) and plasmid transformation to a different extent. In ΔaddAB, ΔrecO, recF15, ΔrecU, ΔruvAB or rec⁺ cells intraspecies CT significantly increases, but the ΔrecD2 mutation reduces, and the ΔrecX, ΔradA or ΔdprA mutation further decreases CT in the Δrok context when compared to rok⁺ cells. These observations support the idea that rok inactivation, by altering the topology of the recipient DNA, differentially affects the integration of homologous DNA in rec-deficient strains, and in minor extent the competent subpopulation size. The impairment of other NAP (Hbsu or LrpC) also increased intra- and interspecies CT (nonself-DNA, ~8% nucleotide sequence divergence) in rec⁺ cells, but differentially reduced both types of CTs in certain rec-deficient strains. We describe that rok inactivation significantly stimulates intra and interspecies CT but differentially reduces them in transformation-deficient cells, perhaps by altering the nucleoid architecture. We extend the observation to other NAPs (Hbsu, LrpC).

Genetic Context Significantly Influences the Maintenance and Evolution of Degenerate Pathways

Understanding the evolution of novel physiological traits is highly relevant for expanding the characterization and manipulation of biological systems. Acquisition of new traits can be achieved through horizontal gene transfer (HGT). Here, we investigate drivers that promote or deter the maintenance of HGT-driven degeneracy, occurring when processes accomplish identical functions through nonidentical components. Subsequent evolution can optimize newly acquired functions; for example, beneficial alleles identified in an engineered Methylorubrum extorquens strain allowed it to utilize a “Foreign” formaldehyde oxidation pathway substituted for its Native pathway for methylotrophic growth. We examined the fitness consequences of interactions between these alleles when they were combined with the Native pathway or both (Dual) pathways. Unlike the Foreign pathway context where they evolved, these alleles were often neutral or deleterious when moved into these alternative genetic backgrounds. However, there were instances where combinations of multiple alleles resulted in higher fitness outcomes than individual allelic substitutions could provide. Importantly, the genetic context accompanying these allelic substitutions significantly altered the fitness landscape, shifting local fitness peaks and restricting the set of accessible evolutionary trajectories. These findings highlight how genetic context can negatively impact the probability of maintaining native and HGT-introduced functions together, making it difficult for degeneracy to evolve. However, in cases where the cost of maintaining degeneracy was mitigated by adding evolved alleles impacting the function of these pathways, we observed rare opportunities for pathway coevolution to occur. Together, our results highlight the importance of genetic context and resulting epistasis in retaining or losing HGT-acquired degenerate functions.

The Effect of Iodine-Containing Nano-Micelles, FS-1, on Antibiotic Resistance, Gene Expression and Epigenetic Modifications in the Genome of Multidrug Resistant MRSA Strain Staphylococcus aureus ATCC BAA-39

Application of supplementary drugs which increase susceptibility of pathogenic bacteria to antibiotics is a promising yet unexplored approach to overcome the global problem of multidrug-resistant infections. The discovery of a new drug, an iodine-containing nano-molecular complex FS-1, which has proven to improve susceptibility to antibiotics in various pathogens, including MRSA strain Staphylococcus aureus ATCC BAA-39TM, allowed studying this phenomenon. Chromosomal DNA and total RNA samples extracted from the FS-1 treated strain (FS) and from the negative control (NC) cultures were sequenced by PacBio SMRT and Ion Torrent technologies, respectively. PacBio DNA reads were used to assemble chromosomal DNA of the NC and FS variants of S. aureus BAA-39 and to perform profiling of epigenetically modified nucleotides. Results of transcriptional profiling, variant calling and detection of epigenetic modifications in the FS variant were compared to the NC variant. Additionally, the genetic alterations caused by the treatment of S. aureus BAA-39 with FS-1 were compared to the results of a similar experiment conducted with another model organism, E. coli ATCC BAA-196. Several commonalities in responses of these phylogenetically distant microorganisms to the treatment with FS-1 were discovered, which included metabolic transition toward anaerobiosis and oxidative/osmotic stress response. S. aureus culture appeared to be more sensitive to FS-1 due to a higher penetrability of cells by iodine bound compounds, which caused carbonyl stress associated with nucleotide damaging by FS-1, abnormal epigenetic modifications and an increased rate of mutations. It was hypothesized that the disrupted pattern of adenine methylated loci within methicillin-resistance chromosome cassettes (SCCmec) may promote excision of this antibiotic resistance determinant from chromosomes while the altered pattern of cytosine methylation was behind the adaptive gene regulation in the culture FS. The selection against the antibiotic resistance in bacterial populations caused by abnormal epigenetic modifications exemplifies possible mechanisms of antibiotic resistance reversion induced by iodine-containing compounds. These finding will facilitate development of therapeutic agents against multidrug-resistant infections.

Genome Complexity Browser: Visualization and quantification of genome variability

Comparative genomics studies may be used to acquire new knowledge regarding genome architecture, which defines the rules for combining sets of genes in the genome of living organisms. Hundreds of thousands of prokaryotic genomes have been sequenced and assembled. However, computational tools capable of simultaneously comparing large numbers of genomes are lacking. We developed the Genome Complexity Browser, a tool that allows the visualization of gene contexts, in a graph-based format, and the quantification of variability for different segments of a genome. The graph-based visualization allows the inspection of changes in gene contents and neighborhoods across hundreds of genomes, simultaneously, which may facilitate the identification of conserved and variable segments of operons or the estimation of the overall variability associated with a particular genome locus. We introduced a measure called complexity, to quantify genome variability. Intraspecies and interspecies comparisons revealed that regions with high complexity values tended to be located in areas that are conserved across different strains and species.

The comparison of genomes among different bacteria and archaea species has revealed that many species frequently exchange genes. Occasionally, such horizontal gene transfer events result in the acquisition of pathogenic properties or antibiotic resistance in the recipient organism. Previously, the probabilities of gene insertions were found to vary, with unequal distributions along a chromosome. At some loci, referred to as hotspots, changes occur with much higher frequencies compared with other regions of the chromosome. We developed a computational method and a software tool, called Genome Complexity Browser, that allows the identification of genome variability hotspots and the visualization of changes. We compared the localization of various hotspots and revealed that some demonstrate conserved localizations, even across species, whereas others are transient. Our tool allows users to visually inspect the patterns of gene changes in graph-based format, which presents the visualization in a format that is both compact and informative.

Selection for ancient periodic motifs that do not impart DNA bending

A ~10-11 bp periodicity in dinucleotides imparting DNA bending, with shorter periods found in organisms with positively-supercoiled DNA and longer periods found in organisms with negatively-supercoiled DNA, was previously suggested to assist in DNA compaction. However, when measured with more robust methods, variation in the observed periods between organisms with different growth temperatures is not consistent with that hypothesis. We demonstrate that dinucleotide periodicity does not arise solely by mutational biases but is under selection. We found variation between genomes in both the period and the suite of dinucleotides that are periodic. Whereas organisms with similar growth temperatures have highly variable periods, differences in periods increase with phylogenetic distance between organisms. In addition, while the suites of dinucleotides under selection for periodicity become more dissimilar among more distantly-related organisms, there is a core set of dinucleotides that are strongly periodic among genomes in all domains of life. Notably, this core set of periodic motifs are not involved in DNA bending. These data indicate that dinucleotide periodicity is an ancient genomic architecture which may play a role in shaping the evolution of genes and genomes.

Conserved Patterns of Symmetric Inversion in the Genome Evolution of Bordetella Respiratory Pathogens

Whooping cough (pertussis), primarily caused by Bordetella pertussis, has resurged in the United States, and circulating strains exhibit considerable chromosome structural fluidity in the form of rearrangement and deletion. The genus Bordetella includes additional pathogenic species infecting various animals, some even causing pertussis-like respiratory disease in humans; however, investigation of their genome evolution has been limited. We studied chromosome structure in complete genome sequences from 167 Bordetella species isolates, as well as 469 B. pertussis isolates, to gain a generalized understanding of rearrangement patterns among these related pathogens. Observed changes in gene order primarily resulted from large inversions and were only detected in species with genomes harboring multicopy insertion sequence (IS) elements, most notably B. holmesii and B. parapertussis While genomes of B. pertussis contain >240 copies of IS481, IS elements appear less numerous in other species and yield less chromosome structural diversity through rearrangement. These data were further used to predict all possible rearrangements between IS element copies present in Bordetella genomes, revealing that only a subset is observed among circulating strains. Therefore, while it appears that rearrangement occurs less frequently in other species than in B. pertussis, these clinically relevant respiratory pathogens likely experience similar mutation of gene order. The resulting chromosome structural fluidity presents both challenges and opportunity for the study of Bordetella respiratory pathogens.IMPORTANCE Bordetella pertussis is the primary agent of whooping cough (pertussis). The Bordetella genus includes additional pathogens of animals and humans, including some that cause pertussis-like respiratory illness. The chromosome of B. pertussis has previously been shown to exhibit considerable structural rearrangement, but insufficient data have prevented comparable investigation in related species. In this study, we analyze chromosome structure variation in several Bordetella species to gain a generalized understanding of rearrangement patterns in this genus. Just as in B. pertussis, we observed inversions in other species that likely result from common mutational processes. We used these data to further predict additional, unobserved inversions, suggesting that specific genome structures may be preferred in each species.