Dynamics of genome rearrangement in bacterial populations

Linear dicentric chromosomes in bacterial natural isolates reveal common constraints for replicon fusion

Multipartite bacterial genome organization can confer advantages, including coordinated gene regulation and faster genome replication, but is challenging to maintain. Agrobacterium tumefaciens lineages often contain a circular chromosome (Ch1), a linear chromosome (Ch2), and multiple plasmids. We previously observed that in some stocks of the C58 lab model, Ch1 and Ch2 were fused into a linear dicentric chromosome. Here we analyzed Agrobacterium natural isolates from the French Collection for Plant-Associated Bacteria and identified two strains distinct from C58 with fused chromosomes. Chromosome conformation capture identified integration junctions that were different from the C58 fusion strain. Genome-wide DNA replication profiling showed that both replication origins remained active. Transposon sequencing revealed that partitioning systems of both chromosome centromeres were essential. Importantly, the site-specific recombinase XerCD is required for the survival of the strains containing the fusion chromosome. Our findings show that replicon fusion occurs in natural environments and that balanced replication arm sizes and proper resolution systems enable the survival of such strains.

IMPORTANCE

Most bacterial genomes are monopartite with a single, circular chromosome. However, some species, like Agrobacterium tumefaciens, carry multiple chromosomes. Emergence of multipartite genomes is often related to adaptation to specific niches, including pathogenesis or symbiosis. Multipartite genomes confer certain advantages; however, maintaining this complex structure can present significant challenges. We previously reported a laboratory-propagated lineage of A. tumefaciens strain C58 in which the circular and linear chromosomes fused to form a single dicentric chromosome. Here we discovered two geographically separated environmental isolates of A. tumefaciens containing fused chromosomes with integration junctions different from the C58 fusion chromosome, revealing the constraints and diversification of this process. We found that balanced replication arm sizes and the repurposing of multimer resolution systems enable the survival and stable maintenance of dicentric chromosomes. These findings reveal how multipartite genomes function across different bacterial species and the role of genomic plasticity in bacterial genetic diversification.

Establishment of a RecET-Assisted CRISPR-Cas12a System for Large Deoxyribonucleic Acid-Fragment Manipulation in Zymomonas mobilis

The lack of effective and efficient genome-editing toolkits for large DNA-fragment manipulation impedes the development of robust cell factories to meet the needs of fast-growing biomanufacturing. Zymomonas mobilis is an important nonmodel polyploid industrial strain with excellent industrial characteristics. Although diverse CRISPR-Cas systems have been established in Z. mobilis for gene deletion, replacement, and ssDNA recombination, it is still challenging to achieve large DNA-fragment manipulation due to its low recombination and repair efficiencies for double-strand DNA breaks. In this study, a RecET-assisted CRISPR-Cas12a genome editing system was developed using a chromosome-borne cas12a and recET-encoded recombinase, as well as an all-in-one editing plasmid. Different promoters were used for recET and cas12a to determine optimal expression. The combination of PB-cas12a_Pt-recET had the highest efficiency of 97.92 ± 2.95% for 9-kb DNA-fragment deletion, which also had efficiencies about 100%, 80%, and 5%, respectively, for the deletion of 9-16, 20-25, and 30 kb DNA fragments. The RecET-assisted CRISPR-Cas12a was further applied for deletions of different large gene clusters and had the potential for efficient pathway knock-in. This study highlights the importance of the Cas12a nuclease expression levels and the combination of the RecET system in improving the double-strand DNA repair capability for large DNA-fragment manipulation in Z. mobilis. The RecET-assisted CRISPR-Cas12a system established in this study provides a versatile and powerful tool for large DNA-fragment manipulation in Z. mobilis, which is beneficial for functional genomic research, strain improvement, as well as the development of synthetic microbial chassis.

Laboratory evolution of the bacterial genome structure through insertion sequence activation

The genome structure fundamentally shapes bacterial physiology, ecology, and evolution. Though insertion sequences (IS) are known drivers of drastic evolutionary changes in the genome structure, the process is typically slow and challenging to observe in the laboratory. Here, we developed a system to accelerate IS-mediated genome structure evolution by introducing multiple copies of a high-activity IS in Escherichia coli. We evolved the bacteria under relaxed neutral conditions, simulating those leading to IS expansion in host-restricted endosymbionts and pathogens. Strains accumulated a median of 24.5 IS insertions and underwent over 5% genome size changes within ten weeks, comparable to decades-long evolution in wild-type strains. The detected interplay of frequent small deletions and rare large duplications updates the view of genome reduction under relaxed selection from a simple consequence of the deletion bias to a nuanced picture including transient expansions. The high IS activity resulted in structural variants of IS and the emergence of composite transposons, illuminating potential evolutionary pathways for ISs and composite transposons. The extensive genome rearrangements we observed establish a baseline for assessing the fitness effects of IS insertions, genome size changes, and rearrangements, advancing our understanding of how mobile elements shape bacterial genomes.

Bacterial genome structural variation: prevalence, mechanisms, and consequences

A vast number of bacterial genome sequences are publicly available. However, the majority were generated using short-read sequencing, producing fragmented assemblies. Long-read sequencing can generate closed assemblies, and they reveal that bacterial genome structure, the order and orientation of genes on the chromosome, is highly variable for many species. Growing evidence suggests that genome structure is a determinant of genome-wide gene expression levels and thus phenotype. We review this developing picture of genome structure variation among bacteria, the challenges for the study of this phenomenon, and its impact on adaptation and evolution, including virulence and infection.

Most bacterial gene families are biased toward specific chromosomal positions

The arrangement of genes along bacterial chromosomes influences their expression through growth rate-dependent gene copy number changes during DNA replication. Although translation- and transcription-related genes often cluster near the origin of replication, the extent of positional biases across gene families remains unclear. We hypothesized that natural selection broadly favors specific chromosomal positions to optimize growth rate-dependent expression. Analyzing 910 bacterial species and proteomics data from Escherichia coli and Bacillus subtilis, we found that about two-thirds of bacterial gene families are positionally biased. Natural selection drives genes mainly toward the origin or terminus of replication, with the strongest selection in fast-growing species. Our findings reveal chromosomal positioning as a fundamental mechanism for coordinating gene expression with growth rate, highlighting evolutionary constraints on bacterial genome architecture.

Linear dicentric chromosomes in bacterial natural isolates reveal common constraints for replicon fusion

Multipartite bacterial genome organization can confer advantages including coordinated gene regulation and faster genome replication but is challenging to maintain. Agrobacterium tumefaciens lineages often contain a circular chromosome (Ch1), a linear chromosome (Ch2), and multiple plasmids. We previously observed that in some stocks of the lab model strain C58, Ch1 and Ch2 were fused into a linear dicentric chromosome. Here we analyzed Agrobacterium natural isolates from the French Collection for Plant-Associated Bacteria (CFBP) and identified two strains with fused chromosomes. Chromosome conformation capture identified integration junctions that were different from the C58 fusion strain. Genome-wide DNA replication profiling showed both replication origins remain active. Transposon sequencing revealed that partitioning systems of both chromosome centromeres are essential. Importantly, the site-specific recombinases XerCD are required for the survival of the strains containing the fusion chromosome. Our findings show that replicon fusion occurs in natural environments and that balanced replication arm sizes and proper resolution systems enable the survival of such strains.

Importance

Most bacterial genomes are monopartite with a single, circular chromosome. But some species, like Agrobacterium tumefaciens, carry multiple chromosomes. Emergence of multipartite genomes is often related to adaptation to specific niches including pathogenesis or symbiosis. Multipartite genomes confer certain advantages, however, maintaining this complex structure can present significant challenges. We previously reported a laboratory-propagated lineage of A. tumefaciens strain C58 in which the circular and linear chromosomes fused to form a single dicentric chromosome. Here we discovered two environmental isolates of A. tumefaciens containing fused chromosomes derived from a different route, revealing the constraints and diversification of this process. We found that balanced replication arm sizes and the repurposing of multimer resolution systems enable the survival and stable maintenance of dicentric chromosomes. These findings help us better understand how multipartite genomes function across different bacterial species and the role of genomic plasticity in bacterial genetic diversification.

Forward-in-time simulation of chromosomal rearrangements: The invisible backbone that sustains long-term adaptation

While chromosomal rearrangements are ubiquitous in all domains of life, very little is known about their evolutionary significance, mostly because, apart from a few specifically studied and well-documented mechanisms (interaction with recombination, gene duplication, etc.), very few models take them into account. As a consequence, we lack a general theory to account for their direct and indirect contributions to evolution. Here, we propose Aevol, a forward-in-time simulation platform specifically dedicated to unravelling the evolutionary significance of chromosomal rearrangements (CR) compared to local mutations (LM). Using the platform, we evolve populations of organisms in four conditions characterized by an increasing diversity of mutational operators-from substitutions alone to a mix of substitutions, InDels and CR-but with a constant global mutational rate. Despite being almost invisible in the phylogeny owing to the scarcity of their fixation in the lineages, we show that CR make a decisive contribution to the evolutionary dynamics by comparing the outcome in these four conditions. As expected, chromosomal rearrangements allow fast expansion of the gene repertoire through gene duplication, but they also reduce the effect of diminishing-returns epistasis, hence sustaining adaptation on the long-run. At last, we show that chromosomal rearrangements tightly regulate the size of the genome through indirect selection for reproductive robustness. Overall, these results confirm the need to improve our theoretical understanding of the contribution of chromosomal rearrangements to evolution and show that dedicated platforms like Aevol can efficiently contribute to this agenda.

Insertion sequence elements and unique symmetrical genomic regions mediate chromosomal inversions in Streptococcus pyogenes

Chromosomal inversions are a phenomenon in many bacterial species, often across the axis of replication. Inversions have been shown to alter gene expression, changing persistence of colonisation and infection following environmental stresses. In Streptococcus pyogenes, inversions have been reported. However, frequency and molecular markers of inversions have not been systematically examined. Here, 249 complete S.pyogenes genomes were analysed using a pangenomic core gene synteny framework to identify sequences associated with inversions. 47% of genomes (118/249) contained at least one inversion, from 23 unique inversion locations. Chromosomal locations enabling inversions were usually associated with mobile elements (insertion sequences n = 9 and prophages n = 7). Two insertion sequences, IS1548 and IS1239, accounted for >80% of insertion sequences and were the only insertion sequences associated with inversions. The most observed inversion location (n = 104 genomes, 88% of genomes with an inversion) occurs between two conserved regions encoding rRNAs, tRNAs and sigma factor genes. The regions are symmetrically placed around the origin of replication forming a unique chromosomal structure in S. pyogenes, relative to other streptococci. Cataloging of the chromosomal location and frequency of inversions can direct dissection of phenotypic changes following chromosomal inversions. The framework used here can be transferred to other bacterial species to characterise chromosomal inversions.

GWarrange: a pre- and post- genome-wide association studies pipeline for detecting phenotype-associated genome rearrangement events

The use of k-mers to capture genetic variation in bacterial genome-wide association studies (bGWAS) has demonstrated its effectiveness in overcoming the plasticity of bacterial genomes by providing a comprehensive array of genetic variants in a genome set that is not confined to a single reference genome. However, little attempt has been made to interpret k-mers in the context of genome rearrangements, partly due to challenges in the exhaustive and high-throughput identification of genome structure and individual rearrangement events. Here, we present GWarrange, a pre- and post-bGWAS processing methodology that leverages the unique properties of k-mers to facilitate bGWAS for genome rearrangements. Repeat sequences are common instigators of genome rearrangements through intragenomic homologous recombination, and they are commonly found at rearrangement boundaries. Using whole-genome sequences, repeat sequences are replaced by short placeholder sequences, allowing the regions flanking repeats to be incorporated into relatively short k-mers. Then, locations of flanking regions in significant k-mers are mapped back to complete genome sequences to visualise genome rearrangements. Four case studies based on two bacterial species (Bordetella pertussis and Enterococcus faecium) and a simulated genome set are presented to demonstrate the ability to identify phenotype-associated rearrangements. GWarrange is available at https://github.com/DorothyTamYiLing/GWarrange.

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.

Pathogenomes and virulence profiles of representative big six non-O157 serogroup Shiga toxin-producing Escherichia coli

Shiga toxin (Stx)-producing Escherichia coli (STEC) of non-O157:H7 serotypes are responsible for global and widespread human food-borne disease. Among these serogroups, O26, O45, O103, O111, O121, and O145 account for the majority of clinical infections and are colloquially referred to as the "Big Six." The "Big Six" strain panel we sequenced and analyzed in this study are reference type cultures comprised of six strains representing each of the non-O157 STEC serogroups curated and distributed by the American Type Culture Collection (ATCC) as a resource to the research community under panel number ATCC MP-9. The application of long- and short-read hybrid sequencing yielded closed chromosomes and a total of 14 plasmids of diverse functions. Through high-resolution comparative phylogenomics, we cataloged the shared and strain-specific virulence and resistance gene content and established the close relationship of serogroup O26 and O103 strains featuring flagellar H-type 11. Virulence phenotyping revealed statistically significant differences in the Stx-production capabilities that we found to be correlated to the strain's individual stx-status. Among the carried Stx1a, Stx2a, and Stx2d phages, the Stx2a phage is by far the most responsive upon RecA-mediated phage mobilization, and in consequence, stx2a + isolates produced the highest-level of toxin in this panel. The availability of high-quality closed genomes for this "Big Six" reference set, including carried plasmids, along with the recorded genomic virulence profiles and Stx-production phenotypes will provide a valuable foundation to further explore the plasticity in evolutionary trajectories in these emerging non-O157 STEC lineages, which are major culprits of human food-borne disease.

Genome mining of Escherichia coli WG5D from drinking water source: unraveling antibiotic resistance genes, virulence factors, and pathogenicity

BACKGROUND

Escherichia coli, a ubiquitous inhabitant of the gut microbiota, has been recognized as an indicator of fecal contamination and a potential reservoir for antibiotic resistance genes. Its prevalence in drinking water sources raises concerns about the potential dissemination of antibiotic resistance within aquatic ecosystems and the subsequent impact on public health. The ability of E. coli to acquire and transfer resistance genes, coupled with the constant exposure to low levels of antibiotics in the environment, underscores the need for comprehensive surveillance and rigorous antimicrobial stewardship strategies to safeguard the quality and safety of drinking water supplies, ultimately mitigating the escalation of antibiotic resistance and its implications for human well-being.

METHODS

WG5D strain, isolated from a drinking water distribution source in North-West Province, South Africa, underwent genomic analysis following isolation on nutrient agar, anaerobic cultivation, and DNA extraction. Paired-end Illumina sequencing with a Nextera XT Library Preparation kit was performed. The assembly, annotation, and subsequent genomic analyses, including phylogenetic analysis using TYGS, pairwise comparisons, and determination of genes related to antimicrobial resistance and virulence, were carried out following standard protocols and tools, ensuring comprehensive insights into the strain's genomic features.

RESULTS

This study explores the notable characteristics of E. coli strain WG5D. This strain stands out because it possesses multiple antibiotic resistance genes, encompassing tetracycline, cephalosporin, vancomycin, and aminoglycoside resistances. Additionally, virulence-associated genes indicate potential heightened pathogenicity, complemented by the identification of mobile genetic elements that underscore its adaptability. The intriguing possibility of bacteriophage involvement and factors contributing to pathogenicity further enriches our understanding. We identified E. coli WG5D as a potential human pathogen associated with a drinking water source in South Africa. The analysis provided several antibiotic resistance-associated genes/mutations and mobile genetic elements. It further identified WG5D as a potential human pathogen. The occurrence of E. coli WG5D raised the awareness of the potential pathogens and the carrying of antibiotic resistance in drinking water.

CONCLUSIONS

The findings of this study have highlighted the advantages of the genomic approach in identifying the bacterial species and antibiotic resistance genes of E. coli and its potential as a human pathogen.

The evolutionary mechanism and function analysis of two subgroups of histamine-producing and non-histamine-producing Tetragenococcus halophilus

Tetragenococcus halophilus is a halophilic bacterium that existed in the fermentation of soy sauce and miso for flavor production and probiotic benefits. However, it is composed of two subgroups, histamine-producing and non-histamine-producing, with the former causing histamine accumulation and offering risks to food safety. Exploring the evolutionary mechanisms and physiological function of histamine-biosynthesis is of significance for understanding the formative mechanism of T. halophilus's strain-specificity and is helpful for microbial control. Using systematic genomic analysis, we found that plasmid acquisition and loss is the evolutionary form resulting in the two subgroups of T. halophilus. Two plasmids, plasmid α with 30 kb and plasmid β with 4 kb existed in histamine-producing T. halophilus. We investigated the whole genetic information and proposed their genetic function in both two plasmids. The acquisition of histamine-producing plasmid enhanced the acid tolerance of histamine-producing T. halophilus but did not affect salt tolerance. More interestingly, we found that the existence of plasmid will promote the co-culture growth of T. halophilus. This study deepens our understanding of the formative mechanism of microbial species diversity, and provides our knowledge of the physiological function of histamine-producing plasmids.

Structural Evolution of Primate Glutamate Dehydrogenase 2 as Revealed by In Silico Predictions and Experimentally Determined Structures

Glutamate dehydrogenase (GDH) interconverts glutamate to a-ketoglutarate and ammonia, interconnecting amino acid and carbohydrate metabolism. In humans, two functional GDH genes, GLUD1 and GLUD2, encode for hGDH1 and hGDH2, respectively. GLUD2 evolved from retrotransposition of the GLUD1 gene in the common ancestor of modern apes. These two isoenzymes are involved in the pathophysiology of human metabolic, neoplastic, and neurodegenerative disorders. The 3D structures of hGDH1 and hGDH2 have been experimentally determined; however, no information is available about the path of GDH2 structure changes during primate evolution. Here, we compare the structures predicted by the AlphaFold Colab method for the GDH2 enzyme of modern apes and their extinct primate ancestors. Also, we analyze the individual effect of amino acid substitutions emerging during primate evolution. Our most important finding is that the predicted structure of GDH2 in the common ancestor of apes was the steppingstone for the structural evolution of primate GDH2s. Two changes with a strong functional impact occurring at the first evolutionary step, Arg443Ser and Gly456Ala, had a destabilizing and stabilizing effect, respectively, making this step the most important one. Subsequently, GDH2 underwent additional modifications that fine-tuned its enzymatic properties to adapt to the functional needs of modern-day primate tissues.

New algorithms for structure informed genome rearrangement

We define two new computational problems in the domain of perfect genome rearrangements, and propose three algorithms to solve them. The rearrangement scenarios modeled by the problems consider Reversal and Block Interchange operations, and a PQ-tree is utilized to guide the allowed operations and to compute their weights. In the first problem, [Formula: see text] ([Formula: see text]), we define the basic structure-informed rearrangement measure. Here, we assume that the gene order members of the gene cluster from which the PQ-tree is constructed are permutations. The PQ-tree representing the gene cluster is ordered such that the series of gene IDs spelled by its leaves is equivalent to that of the reference gene order. Then, a structure-informed genome rearrangement distance is computed between the ordered PQ-tree and the target gene order. The second problem, [Formula: see text] ([Formula: see text]), generalizes [Formula: see text], where the gene order members are not necessarily permutations and the structure informed rearrangement measure is extended to also consider up to [Formula: see text] and [Formula: see text] gene insertion and deletion operations, respectively, when modelling the PQ-tree informed divergence process from the reference gene order to the target gene order. The first algorithm solves [Formula: see text] in [Formula: see text] time and [Formula: see text] space, where [Formula: see text] is the maximum number of children of a node, n is the length of the string and the number of leaves in the tree, and [Formula: see text] and [Formula: see text] are the number of P-nodes and Q-nodes in the tree, respectively. If one of the penalties of [Formula: see text] is 0, then the algorithm runs in [Formula: see text] time and [Formula: see text] space. The second algorithm solves [Formula: see text] in [Formula: see text] time and [Formula: see text] space, where [Formula: see text] is the maximum number of children of a node, n is the length of the string, m is the number of leaves in the tree, [Formula: see text] and [Formula: see text] are the number of P-nodes and Q-nodes in the tree, respectively, and allowing up to [Formula: see text] deletions from the tree and up to [Formula: see text] deletions from the string. The third algorithm is intended to reduce the space complexity of the second algorithm. It solves a variant of the problem (where one of the penalties of [Formula: see text] is 0) in [Formula: see text] time and [Formula: see text] space. The algorithm is implemented as a software tool, denoted MEM-Rearrange, and applied to the comparative and evolutionary analysis of 59 chromosomal gene clusters extracted from a dataset of 1487 prokaryotic genomes.

Examination of Large Chromosomal Inversions in the Genome of Erwinia amylovora Strains Reveals Worldwide Distribution and North America-Specific Types

Erwinia amylovora is a relatively homogeneous species with low genetic diversity at the nucleotide level. However, phenotypic differences and genomic structural variations among E. amylovora strains have been documented. In this study, we identified 10 large chromosomal inversion (LCI) types in the Spiraeoideae-infecting (SI) E. amylovora strains by combining whole genome sequencing and PCR-based molecular markers. It was found that LCIs were mainly caused by homologous recombination events among seven rRNA operons (rrns) in SI E. amylovora strains. Although ribotyping results identified inter- and intra-variations in the internal transcribed spacer (ITS1 and ITS2) regions among rrns, LCIs tend to occur between rrns transcribed in the opposite directions and with the same tRNA content (tRNA-Glu or tRNA-Ile/Ala) in ITS1. Based on the LCI types, physical/estimated replichore imbalance (PRI/ERI) was examined and calculated. Among the 117 SI strains evaluated, the LCI types of Ea1189, CFBP1430, and Ea273 were the most common, with ERI values at 1.31, 7.87, and 4.47°, respectively. These three LCI types had worldwide distribution, whereas the remaining seven LCI types were restricted to North America (or certain regions of the United States). Our results indicated ongoing chromosomal recombination events in the SI E. amylovora population and showed that LCI events are mostly symmetrical, keeping the ERI less than 15°. These findings provide initial evidence about the prevalence of certain LCI types in E. amylovora strains, how LCI occurs, and its potential evolutionary advantage and history, which might help track the movement of the pathogen.

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.

Rearrangement Events on Circular Genomes

Early literature on genome rearrangement modelling views the problem of computing evolutionary distances as an inherently combinatorial one. In particular, attention is given to estimating distances using the minimum number of events required to transform one genome into another. In hindsight, this approach is analogous to early methods for inferring phylogenetic trees from DNA sequences such as maximum parsimony-both are motivated by the principle that the true distance minimises evolutionary change, and both are effective if this principle is a true reflection of reality. Recent literature considers genome rearrangement under statistical models, continuing this parallel with DNA-based methods, with the goal of using model-based methods (for example maximum likelihood techniques) to compute distance estimates that incorporate the large number of rearrangement paths that can transform one genome into another. Crucially, this approach requires one to decide upon a set of feasible rearrangement events and, in this paper, we focus on characterising well-motivated models for signed, uni-chromosomal circular genomes, where the number of regions remains fixed. Since rearrangements are often mathematically described using permutations, we isolate the sets of permutations representing rearrangements that are biologically reasonable in this context, for example inversions and transpositions. We provide precise mathematical expressions for these rearrangements, and then describe them in terms of the set of cuts made in the genome when they are applied. We directly compare cuts to breakpoints, and use this concept to count the distinct rearrangement actions which apply a given number of cuts. Finally, we provide some examples of rearrangement models, and include a discussion of some questions that arise when defining plausible models.

Persistent Salmonella enterica serovar Typhi sub-populations within host interrogated by whole genome sequencing and metagenomics

Salmonella enterica serovar Typhi (S. Typhi) causes typhoid fever and, in some cases, chronic carriage after resolution of acute disease. This study examined sequential isolates of S. Typhi from a single host with persistent asymptomatic infection. These isolates, along with another S. Typhi isolate recovered from a household contact with typhoid fever, were subjected to whole genome sequencing and analysis. In addition, direct sequencing of the bile fluid from the host with persistent infection was also performed. Comparative analysis of isolates revealed three sub-populations of S. Typhi with distinct genetic patterns. Metagenomic sequencing recognised only two of the three sub-populations within the bile fluid. The detection and investigation of insertion sequences IS10R and associated deletions complemented analysis of single nucleotide polymorphisms. These findings improve our understanding of within-host dynamics of S. Typhi in cases of persistent infection and inform epidemiological investigations of transmission events associated with chronic carriers.

An algebraic model for inversion and deletion in bacterial genome rearrangement

Inversions, also sometimes called reversals, are a major contributor to variation among bacterial genomes, with studies suggesting that those involving small numbers of regions are more likely than larger inversions. Deletions may arise in bacterial genomes through the same biological mechanism as inversions, and hence a model that incorporates both is desirable. However, while inversion distances between genomes have been well studied, there has yet to be a model which accounts for the combination of both deletions and inversions. To account for both of these operations, we introduce an algebraic model that utilises partial permutations. This leads to an algorithm for calculating the minimum distance to the most recent common ancestor of two bacterial genomes evolving by inversions (of adjacent regions) and deletions. The algebraic model makes the existing short inversion models more complete and realistic by including deletions, and also introduces new algebraic tools into evolutionary distance problems.

Genomic Diversity among Actinomyces naeslundii Strains and Closely Related Species

The aim of this study was to investigate and clarify the ambiguous taxonomy of Actinomyces naeslundii and its closely related species using state-of-the-art high-throughput sequencing techniques, and, furthermore, to determine whether sub-clusters identified within Actinomyces oris and Actinomyces naeslundii in a previous study by multi locus sequence typing (MLST) using concatenation of seven housekeeping genes should either be classified as subspecies or distinct species. The strains in this study were broadly classified under Actinomyces naeslundii group as A. naeslundii genospecies I and genospecies II. Based on MLST data analysis, these were further classified as A. oris and A. naeslundii. The whole genome sequencing of selected strains of A. oris (n = 17) and A. naeslundii (n = 19) was carried out using Illumina Genome Analyzer IIxe and Roche 454 allowing paired-end and single-reads sequencing, respectively. The sequences obtained were aligned using CLC Genomic workbench version 5.1 and annotated using RAST (Rapid Annotation using Subsystem Technology) release version 59 accessible online. Additionally, genomes of seven publicly available strains of Actinomyces (k20, MG1, c505, OT175, OT171, OT170, and A. johnsonii) were also included. Comparative genomic analysis (CGA) using Mauve, Progressive Mauve, gene-by-gene, Core, and Pan Genome, and finally Digital DNA-DNA homology (DDH) analysis was carried out. DDH values were obtained using in silico genome-genome comparison. Evolutionary analysis using ClonalFrame was also undertaken. The mutation and recombination events were compared using chi-square test among A. oris and A. naeslundii isolates (analysis methods are not included in the study). CGA results were consistent with previous traditional classification using MLST. It was found that strains of Actinomyces k20, MG1, c505, and OT175 clustered in A. oris group of isolates, while OT171, OT170, and A. johnsonii appeared as separate branches. Similar clustering to MLST was observed for other isolates. The mutation and recombination events were significantly higher in A. oris than A. naeslundii, highlighting the diversity of A. oris strains in the oral cavity. These findings suggest that A. oris forms six distinct groups, whereas A. naeslundii forms three. The correct designation of isolates will help in the identification of clinical Actinomyces isolates found in dental plaque. Easily accessible online genomic sequence data will also accelerate the investigation of the biochemical characterisation and pathogenesis of this important group of micro-organisms.

Impact of Salmonella genome rearrangement on gene expression

In addition to nucleotide variation, many bacteria also undergo changes at a much larger scale via rearrangement of their genome structure (GS) around long repeat sequences. These rearrangements result in genome fragments shifting position and/or orientation in the genome without necessarily affecting the underlying nucleotide sequence. To date, scalable techniques have not been applied to GS identification, so it remains unclear how extensive this variation is and the extent of its impact upon gene expression. However, the emergence of multiplexed, long-read sequencing overcomes the scale problem, as reads of several thousand bases are routinely produced that can span long repeat sequences to identify the flanking chromosomal DNA, allowing GS identification. Genome rearrangements were generated in Salmonella enterica serovar Typhi through long-term culture at ambient temperature. Colonies with rearrangements were identified via long-range PCR and subjected to long-read nanopore sequencing to confirm genome variation. Four rearrangements were investigated for differential gene expression using transcriptomics. All isolates with changes in genome arrangement relative to the parent strain were accompanied by changes in gene expression. Rearrangements with similar fragment movements demonstrated similar changes in gene expression. The most extreme rearrangement caused a large imbalance between the origin and terminus of replication and was associated with differential gene expression as a factor of distance moved toward or away from the origin of replication. Genome structure variation may provide a mechanism through which bacteria can quickly adapt to new environments and warrants routine assessment alongside traditional nucleotide-level measures of variation.

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.

Genome-guided prediction of acid resistance mechanisms in acidophilic methanotrophs of phylogenetically deep-rooted Verrucomicrobia isolated from geothermal environments

Verrucomicrobia are a group of microorganisms that have been proposed to be deeply rooted in the Tree of Life. Some are methanotrophs that oxidize the potent greenhouse gas methane and are thus important in decreasing atmospheric concentrations of the gas, potentially ameliorating climate change. They are widespread in various environments including soil and fresh or marine waters. Recently, a clade of extremely acidophilic Verrucomicrobia, flourishing at pH < 3, were described from high-temperature geothermal ecosystems. This novel group could be of interest for studies about the emergence of life on Earth and to astrobiologists as homologs for possible extraterrestrial life. In this paper, we describe predicted mechanisms for survival of this clade at low pH and suggest its possible evolutionary trajectory from an inferred neutrophilic ancestor. Extreme acidophiles are defined as organisms that thrive in extremely low pH environments (≤ pH 3). Many are polyextremophiles facing high temperatures and high salt as well as low pH. They are important to study for both providing fundamental insights into biological mechanisms of survival and evolution in such extreme environments and for understanding their roles in biotechnological applications such as industrial mineral recovery (bioleaching) and mitigation of acid mine drainage. They are also, potentially, a rich source of novel genes and pathways for the genetic engineering of microbial strains. Acidophiles of the Verrucomicrobia phylum are unique as they are the only known aerobic methanotrophs that can grow optimally under acidic (pH 2–3) and moderately thermophilic conditions (50–60°C). Three moderately thermophilic genera, namely Methylacidiphilum, Methylacidimicrobium, and Ca. Methylacidithermus, have been described in geothermal environments. Most of the investigations of these organisms have focused on their methane oxidizing capabilities (methanotrophy) and use of lanthanides as a protein cofactor, with no extensive study that sheds light on the mechanisms that they use to flourish at extremely low pH. In this paper, we extend the phylogenetic description of this group of acidophiles using whole genome information and we identify several mechanisms, potentially involved in acid resistance, including “first line of defense” mechanisms that impede the entry of protons into the cell. These include the presence of membrane-associated hopanoids, multiple copies of the outer membrane protein (Slp), and inner membrane potassium channels (kup, kdp) that generate a reversed membrane potential repelling the intrusion of protons. Acidophilic Verrucomicrobia also display a wide array of proteins potentially involved in the “second line of defense” where protons that evaded the first line of defense and entered the cell are expelled or neutralized, such as the glutamate decarboxylation (gadAB) and phosphate-uptake systems. An exclusive N-type ATPase F₀-F₁ was identified only in acidophiles of Verrucomicrobia and is predicted to be a specific adaptation in these organisms. Phylogenetic analyses suggest that many predicted mechanisms are evolutionarily conserved and most likely entered the acidophilic lineage of Verrucomicrobia by vertical descent from a common ancestor. However, it is likely that some defense mechanisms such as gadA and kup entered the acidophilic Verrucomicrobia lineage by horizontal gene transfer.

Analysis of the Taxonomy, Synteny, and Virulence Factors for Soft Rot Pathogen Pectobacterium aroidearum in Amorphophallus konjac Using Comparative Genomics

Bacterial soft rot is a devastating disease for a wide range of crops, vegetables, and ornamental plants including konjac (Amorphophallus konjac). However, the pangenome and genomic plasticity of the konjac soft rot pathogens is little explored. In this study, we reported the complete genome sequences of 11 bacterial isolates that can cause typical soft rot symptoms in konjac by in vitro and in vivo pathogenicity tests. Based on in silico DNA-DNA hybridization, average nucleotide identity and phylogenomic analysis, all 11 isolates were determined to be Pectobacterium aroidearum. In addition, synteny analysis of these genomes revealed considerable chromosomal inversions, one of which is triggered by homologous recombination of ribose operon. Pangenome analysis and COG enrichment analysis showed that the pangenome of P. aroidearum is open and that accessory genes are enriched in replication, recombination, and repair. Variations in type IV secretion system and type VI secretion system were found, while plant cell wall degrading enzymes were conserved. Furthermore, sequence analyses also provided evidence for the presence of a type V secretion system in Pectobacterium. These findings advance our understanding of the pathogenicity determinants, genomic plasticity, and evolution of P. aroidearum.

The Microbiology of Non-aeruginosa Pseudomonas Isolated From Adults With Cystic Fibrosis: Criteria to Help Determine the Clinical Significance of Non-aeruginosa Pseudomonas in CF Lung Pathology

Introduction: There is a paucity of reports on non-aeruginosa Pseudomonas (NAPs) in cystic fibrosis, hence this study wished 1). to examine the diversity/frequency of NAPs in an adult CF population, 2) to compare/contrast the microbiology and genomics of NAPs to P. aeruginosa and 3) to propose clinical and laboratory criteria to help determine their clinical significance in CF lung pathology.

Materials and Methods: Microbiological data was examined from 100 adult patients with cystic fibrosis from birth to present (31/12/2021), equating to 2455 patient years. 16S rDNA phylogenetic relatedness of NAPs was determined, as well as bioinformatical comparison of whole genomes of P. aeruginosa against P. fluorescens.

Results: Ten species were isolated from this patient cohort during this time period, with three species, i.e., P. fluorescens, P. putida and P. stutzeri, accounting for the majority (87.5%) of non-aeruginosa reports. This is the first report of the isolation of P. fragi, P. nitroreducens, P. oryzihabitans and P. veronii in patients with cystic fibrosis. The mean time to first detection of any non-aeruginosa species was 183 months (15.25 years) [median = 229 months (19.1 years)], with a range from 11 months to 338 months (28.2 years). Several of the NAPs were closely related to P. aeruginosa.

Discussion: NAPs were isolated infrequently and were transient colonisers of the CF airways, in those patients with CF in which they were isolated. A set of ten clinical and laboratory criteria are proposed to provide key indicators, as to the clinical importance of the non-aeruginosa species isolated.

Genomic structural plasticity of rodent-associated Bartonella in nature

Rodent‐associated Bartonella species have shown a remarkable genetic diversity and pathogenic potential. To further explore the extent of the natural intraspecific genomic variation and its potential role as an evolutionary driver, we focused on a single genetically diverse Bartonella species, Bartonella krasnovii, which circulates among gerbils and their associated fleas. Twenty genomes from 16 different B. krasnovii genotypes were fully characterized through a genome sequencing assay (using short and long read sequencing), pulse field gel electrophoresis (PFGE), and PCR validation. Genomic analyses were performed in comparison to the B. krasnovii strain OE 1–1. While, single nucleotide polymorphism represented only a 0.3% of the genome variation, structural diversity was identified in these genomes, with an average of 51 ± 24 structural variation (SV) events per genome. Interestingly, a large proportion of the SVs (>40%) was associated with prophages. Further analyses revealed that most of the SVs, and prophage insertions were found at the chromosome replication termination site (ter), suggesting this site as a plastic zone of the B. krasnovii chromosome. Accordingly, six genomes were found to be unbalanced, and essential genes near the ter showed a shift between the leading and lagging strands, revealing the SV effect on these genomes. In summary, our findings demonstrate the extensive genomic diversity harbored by wild B. krasnovii strains and suggests that its diversification is initially promoted by structural changes, probably driven by phages. These events may constantly feed the system with novel genotypes that ultimately lead to inter‐ and intraspecies competition and adaptation.

Dynamic interactions between mega symbiosis ICEs and bacterial chromosomes maintain genome architecture

Acquisition of mobile genetic elements can confer novel traits to bacteria. Some integrative and conjugative elements confer upon members of Bradyrhizobium the capacity to fix nitrogen in symbiosis with legumes. These so-called symbiosis integrative conjugative elements (symICEs) can be extremely large and vary as monopartite and polypartite configurations within chromosomes of related strains. These features are predicted to impose fitness costs and have defied explanation. Here, we show that chromosome architecture is largely conserved despite diversity in genome composition, variations in locations of attachment sites recognized by integrases of symICEs, and differences in large-scale chromosomal changes that occur upon integration. Conversely, many simulated nonnative chromosome–symICE combinations are predicted to result in lethal deletions or disruptions to architecture. Findings suggest that there is compatibility between chromosomes and symICEs. We hypothesize that the size and structural flexibility of symICEs are important for generating combinations that maintain chromosome architecture across a genus of nitrogen-fixing bacteria with diverse and dynamic genomes.

Directed yeast genome evolution by controlled introduction of trans-chromosomic structural variations

Naturally occurring structural variations (SVs) are a considerable source of genomic variation that can reshape the 3D architecture of chromosomes. Controllable methods aimed at introducing the complex SVs and their related molecular mechanisms have remained farfetched. In this study, an SV-prone yeast strain was developed using Synthetic Chromosome Rearrangement and Modification by LoxP-mediated Evolution (SCRaMbLE) technology with two synthetic chromosomes, namely synV and synX. The biosynthesis of astaxanthin is used as a readout and a proof of concept for the application of SVs in industries. Our findings showed that complex SVs, including a pericentric inversion and a trans-chromosome translocation between synV and synX, resulted in two neo-chromosomes and a 2.7-fold yield of astaxanthin. Also, genetic targets were mapped, which resulted in a higher astaxanthin yield, thus demonstrating the SVs' ability to reorganize genetic information along the chromosomes. The rational design of trans-chromosome translocation and pericentric inversion enabled precise induction of these phenomena. Collectively, this study provides an effective tool to not only accelerate the directed genome evolution but also to reveal the mechanistic insight of complex SVs for altering phenotypes.

Positive selection during niche adaptation results in large-scale and irreversible rearrangement of chromosomal gene order in bacteria

Analysis of bacterial genomes shows that, whereas diverse species share many genes in common, their linear order on the chromosome is often not conserved. Whereas rearrangements in gene order could occur by genetic drift, an alternative hypothesis is rearrangement driven by positive selection during niche adaptation (SNAP). Here, we provide the first experimental support for the SNAP hypothesis. We evolved Salmonella to adapt to growth on malate as the sole carbon source and followed the evolutionary trajectories. The initial adaptation to growth in the new environment involved the duplication of 1.66 Mb, corresponding to one-third of the Salmonella chromosome. This duplication is selected to increase the copy number of a single gene, dctA, involved in the uptake of malate. Continuing selection led to the rapid loss or mutation of duplicate genes from either copy of the duplicated region. After 2000 generations, only 31% of the originally duplicated genes remained intact and the gene order within the Salmonella chromosome has been significantly and irreversibly altered. These results experientially validate predictions made by the SNAP hypothesis and show that SNAP can be a strong driving force for rearrangements in chromosomal gene order.

A symmetry-inclusive algebraic approach to genome rearrangement

Of the many modern approaches to calculating evolutionary distance via models of genome rearrangement, most are tied to a particular set of genomic modeling assumptions and to a restricted class of allowed rearrangements. The "position paradigm", in which genomes are represented as permutations signifying the position (and orientation) of each region, enables a refined model-based approach, where one can select biologically plausible rearrangements and assign to them relative probabilities/costs. Here, one must further incorporate any underlying structural symmetry of the genomes into the calculations and ensure that this symmetry is reflected in the model. In our recently-introduced framework of genome algebras, each genome corresponds to an element that simultaneously incorporates all of its inherent physical symmetries. The representation theory of these algebras then provides a natural model of evolution via rearrangement as a Markov chain. Whilst the implementation of this framework to calculate distances for genomes with "practical" numbers of regions is currently computationally infeasible, we consider it to be a significant theoretical advance: one can incorporate different genomic modeling assumptions, calculate various genomic distances, and compare the results under different rearrangement models. The aim of this paper is to demonstrate some of these features.

Genome structural variation in Escherichia coli O157:H7

The human zoonotic pathogen Escherichia coli O157:H7 is defined by its extensive prophage repertoire including those that encode Shiga toxin, the factor responsible for inducing life-threatening pathology in humans. As well as introducing genes that can contribute to the virulence of a strain, prophage can enable the generation of large-chromosomal rearrangements (LCRs) by homologous recombination. This work examines the types and frequencies of LCRs across the major lineages of the O157:H7 serotype. We demonstrate that LCRs are a major source of genomic variation across all lineages of E. coli O157:H7 and by using both optical mapping and Oxford Nanopore long-read sequencing prove that LCRs are generated in laboratory cultures started from a single colony and that these variants can be recovered from colonized cattle. LCRs are biased towards the terminus region of the genome and are bounded by specific prophages that share large regions of sequence homology associated with the recombinational activity. RNA transcriptional profiling and phenotyping of specific structural variants indicated that important virulence phenotypes such as Shiga-toxin production, type-3 secretion and motility can be affected by LCRs. In summary, E. coli O157:H7 has acquired multiple prophage regions over time that act to continually produce structural variants of the genome. These findings raise important questions about the significance of this prophage-mediated genome contingency to enhance adaptability between environments.

PaReBrick: PArallel REarrangements and BReaks identification toolkit

MOTIVATION

High plasticity of bacterial genomes is provided by numerous mechanisms including horizontal gene transfer and recombination via numerous flanking repeats. Genome rearrangements such as inversions, deletions, insertions and duplications may independently occur in different strains, providing parallel adaptation or phenotypic diversity. Specifically, such rearrangements might be responsible for virulence, antibiotic resistance and antigenic variation. However, identification of such events requires laborious manual inspection and verification of phyletic pattern consistency.

RESULTS

Here, we define the term 'parallel rearrangements' as events that occur independently in phylogenetically distant bacterial strains and present a formalization of the problem of parallel rearrangements calling. We implement an algorithmic solution for the identification of parallel rearrangements in bacterial populations as a tool PaReBrick. The tool takes a collection of strains represented as a sequence of oriented synteny blocks and a phylogenetic tree as input data. It identifies rearrangements, tests them for consistency with a tree, and sorts the events by their parallelism score. The tool provides diagrams of the neighbors for each block of interest, allowing the detection of horizontally transferred blocks or their extra copies and the inversions in which copied blocks are involved. We demonstrated PaReBrick's efficiency and accuracy and showed its potential to detect genome rearrangements responsible for pathogenicity and adaptation in bacterial genomes.

AVAILABILITY AND IMPLEMENTATION

PaReBrick is written in Python and is available on GitHub: https://github.com/ctlab/parallel-rearrangements.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Whole-genome microsynteny-based phylogeny of angiosperms

Plant genomes vary greatly in size, organization, and architecture. Such structural differences may be highly relevant for inference of genome evolution dynamics and phylogeny. Indeed, microsynteny-the conservation of local gene content and order-is recognized as a valuable source of phylogenetic information, but its use for the inference of large phylogenies has been limited. Here, by combining synteny network analysis, matrix representation, and maximum likelihood phylogenetic inference, we provide a way to reconstruct phylogenies based on microsynteny information. Both simulations and use of empirical data sets show our method to be accurate, consistent, and widely applicable. As an example, we focus on the analysis of a large-scale whole-genome data set for angiosperms, including more than 120 available high-quality genomes, representing more than 50 different plant families and 30 orders. Our 'microsynteny-based' tree is largely congruent with phylogenies proposed based on more traditional sequence alignment-based methods and current phylogenetic classifications but differs for some long-contested and controversial relationships. For instance, our synteny-based tree finds Vitales as early diverging eudicots, Saxifragales within superasterids, and magnoliids as sister to monocots. We discuss how synteny-based phylogenetic inference can complement traditional methods and could provide additional insights into some long-standing controversial phylogenetic relationships.

High Rates of Genome Rearrangements and Pathogenicity of Shigella spp

Shigella are pathogens originating within the Escherichia lineage but frequently classified as a separate genus. Shigella genomes contain numerous insertion sequences (ISs) that lead to pseudogenisation of affected genes and an increase of non-homologous recombination. Here, we study 414 genomes of E. coli and Shigella strains to assess the contribution of genomic rearrangements to Shigella evolution. We found that Shigella experienced exceptionally high rates of intragenomic rearrangements and had a decreased rate of homologous recombination compared to pathogenic and non-pathogenic E. coli. The high rearrangement rate resulted in independent disruption of syntenic regions and parallel rearrangements in different Shigella lineages. Specifically, we identified two types of chromosomally encoded E3 ubiquitin-protein ligases acquired independently by all Shigella strains that also showed a high level of sequence conservation in the promoter and further in the 5′-intergenic region. In the only available enteroinvasive E. coli (EIEC) strain, which is a pathogenic E. coli with a phenotype intermediate between Shigella and non-pathogenic E. coli, we found a rate of genome rearrangements comparable to those in other E. coli and no functional copies of the two Shigella-specific E3 ubiquitin ligases. These data indicate that the accumulation of ISs influenced many aspects of genome evolution and played an important role in the evolution of intracellular pathogens. Our research demonstrates the power of comparative genomics-based on synteny block composition and an important role of non-coding regions in the evolution of genomic islands.

Comparative in silico genome analysis of Clostridium perfringens unravels stable phylogroups with different genome characteristics and pathogenic potential

Clostridium perfringens causes a plethora of devastating infections, with toxin production being the underlying mechanism of pathogenicity in various hosts. Genomic analyses of 206 public-available C. perfringens strains´ sequence data identified a substantial degree of genomic variability in respect to episome content, chromosome size and mobile elements. However, the position and order of the local collinear blocks on the chromosome showed a considerable degree of preservation. The strains were divided into five stable phylogroups (I–V). Phylogroup I contained human food poisoning strains with chromosomal enterotoxin (cpe) and a Darmbrand strain characterized by a high frequency of mobile elements, a relatively small genome size and a marked loss of chromosomal genes, including loss of genes encoding virulence traits. These features might correspond to the adaptation of these strains to a particular habitat, causing human foodborne illnesses. This contrasts strains that belong to phylogroup II where the genome size points to the acquisition of genetic material. Most strains of phylogroup II have been isolated from enteric lesions in horses and dogs. Phylogroups III, IV and V are heterogeneous groups containing a variety of different strains, with phylogroup III being the most abundant (65.5%). In conclusion, C. perfringens displays five stable phylogroups reflecting different disease involvements, prompting further studies on the evolution of this highly important pathogen.

Comparative Analysis of Genome of Ehrlichia sp. HF, a Model Bacterium to Study Fatal Human Ehrlichiosis

BACKGROUND

The genus Ehrlichia consists of tick-borne obligatory intracellular bacteria that can cause deadly diseases of medical and agricultural importance. Ehrlichia sp. HF, isolated from Ixodes ovatus ticks in Japan [also referred to as I. ovatus Ehrlichia (IOE) agent], causes acute fatal infection in laboratory mice that resembles acute fatal human monocytic ehrlichiosis caused by Ehrlichia chaffeensis. As there is no small laboratory animal model to study fatal human ehrlichiosis, Ehrlichia sp. HF provides a needed disease model. However, the inability to culture Ehrlichia sp. HF and the lack of genomic information have been a barrier to advance this animal model. In addition, Ehrlichia sp. HF has several designations in the literature as it lacks a taxonomically recognized name.

RESULTS

We stably cultured Ehrlichia sp. HF in canine histiocytic leukemia DH82 cells from the HF strain-infected mice, and determined its complete genome sequence. Ehrlichia sp. HF has a single double-stranded circular chromosome of 1,148,904 bp, which encodes 866 proteins with a similar metabolic potential as E. chaffeensis. Ehrlichia sp. HF encodes homologs of all virulence factors identified in E. chaffeensis, including 23 paralogs of P28/OMP-1 family outer membrane proteins, type IV secretion system apparatus and effector proteins, two-component systems, ankyrin-repeat proteins, and tandem repeat proteins. Ehrlichia sp. HF is a novel species in the genus Ehrlichia, as demonstrated through whole genome comparisons with six representative Ehrlichia species, subspecies, and strains, using average nucleotide identity, digital DNA-DNA hybridization, and core genome alignment sequence identity.

CONCLUSIONS

The genome of Ehrlichia sp. HF encodes all known virulence factors found in E. chaffeensis, substantiating it as a model Ehrlichia species to study fatal human ehrlichiosis. Comparisons between Ehrlichia sp. HF and E. chaffeensis will enable identification of in vivo virulence factors that are related to host specificity, disease severity, and host inflammatory responses. We propose to name Ehrlichia sp. HF as Ehrlichia japonica sp. nov. (type strain HF), to denote the geographic region where this bacterium was initially isolated.

A Path-Deformation Framework for Determining Weighted Genome Rearrangement Distance

Measuring the distance between two bacterial genomes under the inversion process is usually done by assuming all inversions to occur with equal probability. Recently, an approach to calculating inversion distance using group theory was introduced, and is effective for the model in which only very short inversions occur. In this paper, we show how to use the group-theoretic framework to establish minimal distance for any weighting on the set of inversions, generalizing previous approaches. To do this we use the theory of rewriting systems for groups, and exploit the Knuth–Bendix algorithm, the first time this theory has been introduced into genome rearrangement problems. The central idea of the approach is to use existing group theoretic methods to find an initial path between two genomes in genome space (for instance using only short inversions), and then to deform this path to optimality using a confluent system of rewriting rules generated by the Knuth–Bendix algorithm.

Network Rewiring: Physiological Consequences of Reciprocally Exchanging the Physical Locations and Growth-Phase-Dependent Expression Patterns of the Salmonella fis and dps Genes

We assessed the impact on Salmonella physiology of reciprocally translocating the genes encoding the Fis and Dps nucleoid-associated proteins (NAPs) and of inverting their growth-phase production patterns such that Fis was produced in stationary phase (like Dps) and Dps was produced in exponential phase (like Fis). Changes to peak binding of Fis were detected by ChIP-seq on the chromosome, as were widespread impacts on the transcriptome, especially when Fis production mimicked Dps production. Virulence gene expression and the expression of a virulence phenotype were altered. Overall, these radical changes to NAP gene expression were well tolerated, revealing the robust and well-buffered nature of global gene regulation networks in the bacterium.

The Fis nucleoid-associated protein controls the expression of a large and diverse regulon of genes in Gram-negative bacteria. Fis production is normally maximal in bacteria during the early exponential phase of batch culture growth, becoming almost undetectable by the onset of stationary phase. We tested the effect on the Fis regulatory network in Salmonella of moving the complete fis gene from its usual location near the origin of chromosomal replication to the position normally occupied by the dps gene in the right macrodomain of the chromosome, and vice versa, creating the gene exchange (GX) strain. In a parallel experiment, we tested the effect of rewiring the Fis regulatory network by placing the fis open reading frame under the control of the stationary-phase-activated dps promoter at the dps genetic location within the right macrodomain, and vice versa, creating the open reading frame exchange (OX) strain. Chromatin immunoprecipitation sequencing (ChIP-seq) was used to measure global Fis protein binding levels and to determine gene expression patterns. Strain GX showed few changes compared with the wild type, although we did detect increased Fis binding at Ter, accompanied by reduced binding at Ori. Strain OX displayed a more pronounced version of this distorted Fis protein-binding pattern together with numerous alterations in the expression of genes in the Fis regulon. OX, but not GX, had a reduced ability to infect cultured mammalian cells. These findings illustrate the inherent robustness of the Fis regulatory network with respect to the effects of rewiring based on gene repositioning alone and emphasize the importance of fis expression signals in phenotypic determination.

Mycoavidus sp. Strain B2-EB: Comparative Genomics Reveals Minimal Genomic Features Required by a Cultivable Burkholderiaceae-Related Endofungal Bacterium

Obligate bacterial endosymbionts are critical to the existence of many eukaryotes. Such endobacteria are usually characterized by reduced genomes and metabolic dependence on the host, which may cause difficulty in isolating them in pure cultures. Family Burkholderiaceae-related endofungal bacteria affiliated with the Mycoavidus-Glomeribacter clade can be associated with the fungal subphyla Mortierellomycotina and Glomeromycotina. In this study, a cultivable endosymbiotic bacterium, Mycoavidus sp. strain B2-EB, present in the fungal host Mortierella parvispora was obtained successfully. The B2-EB genome (1.88 Mb) represents the smallest genome among the endofungal bacterium Mycoavidus cysteinexigens (2.64-2.80 Mb) of Mortierella elongata and the uncultured endosymbiont "Candidatus Glomeribacter gigasporarum" (1.37 to 2.36 Mb) of arbuscular mycorrhizal fungi. Despite a reduction in genome size, strain B2-EB displays a high genome completeness, suggesting a nondegenerative reduction in the B2-EB genome. Compared with a large proportion of transposable elements (TEs) in other known Mycoavidus genomes (7.2 to 11.5% of the total genome length), TEs accounted for only 2.4% of the B2-EB genome. This pattern, together with a high proportion of single-copy genes in the B2-EB genome, suggests that the B2-EB genome reached a state of relative evolutionary stability. These results represent the most streamlined structure among the cultivable endofungal bacteria and suggest the minimal genome features required by both an endofungal lifestyle and artificial culture. This study allows us to understand the genome evolution of Burkholderiaceae-related endosymbionts and to elucidate microbiological interactions.IMPORTANCE This study attempted the isolation of a novel endobacterium, Mycoavidus sp. B2-EB (JCM 33615), harbored in the fungal host Mortierella parvispora E1425 (JCM 39028). We report the complete genome sequence of this strain, which possesses a reduced genome size with relatively high genome completeness and a streamlined genome structure. The information indicates the minimal genomic features required by both the endofungal lifestyle and artificial cultivation, which furthers our understanding of genome reduction in fungal endosymbionts and extends the culture resources for biotechnological development on engineering synthetic microbiomes.

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.

The SNAP hypothesis: Chromosomal rearrangements could emerge from positive Selection during Niche Adaptation

The relative linear order of most genes on bacterial chromosomes is not conserved over evolutionary timescales. One explanation is that selection is weak, allowing recombination to randomize gene order by genetic drift. However, most chromosomal rearrangements are deleterious to fitness. In contrast, we propose the hypothesis that rearrangements in gene order are more likely the result of selection during niche adaptation (SNAP). Partial chromosomal duplications occur very frequently by recombination between direct repeat sequences. Duplicated regions may contain tens to hundreds of genes and segregate quickly unless maintained by selection. Bacteria exposed to non-lethal selections (for example, a requirement to grow on a poor nutrient) can adapt by maintaining a duplication that includes a gene that improves relative fitness. Further improvements in fitness result from the loss or inactivation of non-selected genes within each copy of the duplication. When genes that are essential in single copy are lost from different copies of the duplication, segregation is prevented even if the original selection is lifted. Functional gene loss continues until a new genetic equilibrium is reached. The outcome is a rearranged gene order. Mathematical modelling shows that this process of positive selection to adapt to a new niche can rapidly drive rearrangements in gene order to fixation. Signature features (duplication formation and divergence) of the SNAP model were identified in natural isolates from multiple species showing that the initial two steps in the SNAP process can occur with a remarkably high frequency. Further bioinformatic and experimental analyses are required to test if and to which extend the SNAP process acts on bacterial genomes.

All life on earth has evolved from a universal common ancestor with a specific order of genes on the chromosome. This order is not maintained in modern species and the standard hypothesis is that changes reflect a lack of strong selection on gene order. Here, we propose an alternative hypothesis, SNAP. The occupation of a novel environment by bacteria is generally a trade-off situation. For example, while the bacteria may not be adapted to grow well under the new conditions, they may benefit by not having to share available resources with other microorganisms. Bacterial populations frequently acquire duplications of chromosomal segments containing genes that can help them adapt to a new environment. Other genes that are also duplicated are not required in two copies so that over time a superfluous copy can be lost. Eventually, the process of duplication and gene loss can lead to the rearrangement of the gene order in the chromosomal segment. The major benefit of this model over the standard hypothesis is that the process is driven by positive selection and can reach fixation rapidly.

The number of k-mer matches between two DNA sequences as a function of k and applications to estimate phylogenetic distances

We study the number Nk of length-k word matches between pairs of evolutionarily related DNA sequences, as a function of k. We show that the Jukes-Cantor distance between two genome sequences-i.e. the number of substitutions per site that occurred since they evolved from their last common ancestor-can be estimated from the slope of a function F that depends on Nk and that is affine-linear within a certain range of k. Integers kmin and kmax can be calculated depending on the length of the input sequences, such that the slope of F in the relevant range can be estimated from the values F(kmin) and F(kmax). This approach can be generalized to so-called Spaced-word Matches (SpaM), where mismatches are allowed at positions specified by a user-defined binary pattern. Based on these theoretical results, we implemented a prototype software program for alignment-free sequence comparison called Slope-SpaM. Test runs on real and simulated sequence data show that Slope-SpaM can accurately estimate phylogenetic distances for distances up to around 0.5 substitutions per position. The statistical stability of our results is improved if spaced words are used instead of contiguous words. Unlike previous alignment-free methods that are based on the number of (spaced) word matches, Slope-SpaM produces accurate results, even if sequences share only local homologies.

The Progress of Multi-Omics Technologies: Determining Function in Lactic Acid Bacteria Using a Systems Level Approach

Lactic Acid Bacteria (LAB) have long been recognized as having a significant impact ranging from commercial to health domains. A vast amount of research has been carried out on these microbes, deciphering many of the pathways and components responsible for these desirable effects. However, a large proportion of this functional information has been derived from a reductionist approach working with pure culture strains. This provides limited insight into understanding the impact of LAB within intricate systems such as the gut microbiome or multi strain starter cultures. Whole genome sequencing of strains and shotgun metagenomics of entire systems are powerful techniques that are currently widely used to decipher function in microbes, but they also have their limitations. An available genome or metagenome can provide an image of what a strain or microbiome, respectively, is potentially capable of and the functions that they may carry out. A top-down, multi-omics approach has the power to resolve the functional potential of an ecosystem into an image of what is being expressed, translated and produced. With this image, it is possible to see the real functions that members of a system are performing and allow more accurate and impactful predictions of the effects of these microorganisms. This review will discuss how technological advances have the potential to increase the yield of information from genomics, transcriptomics, proteomics and metabolomics. The potential for integrated omics to resolve the role of LAB in complex systems will also be assessed. Finally, the current software approaches for managing these omics data sets will be discussed.

Comparative Genomics of Streptococcus thermophilus Support Important Traits Concerning the Evolution, Biology and Technological Properties of the Species

Streptococcus thermophilus is a major starter for the dairy industry with great economic importance. In this study we analyzed 23 fully sequenced genomes of S. thermophilus to highlight novel aspects of the evolution, biology and technological properties of this species. Pan/core genome analysis revealed that the species has an important number of conserved genes and that the pan genome is probably going to be closed soon. According to whole genome phylogeny and average nucleotide identity (ANI) analysis, most S. thermophilus strains were grouped in two major clusters (i.e., clusters A and B). More specifically, cluster A includes strains with chromosomes above 1.83 Mbp, while cluster B includes chromosomes below this threshold. This observation suggests that strains belonging to the two clusters may be differentiated by gene gain or gene loss events. Furthermore, certain strains of cluster A could be further subdivided in subgroups, i.e., subgroup I (ASCC 1275, DGCC 7710, KLDS SM, MN-BM-A02, and ND07), II (MN-BM-A01 and MN-ZLW-002), III (LMD-9 and SMQ-301), and IV (APC151 and ND03). In cluster B certain strains formed one distinct subgroup, i.e., subgroup I (CNRZ1066, CS8, EPS, and S9). Clusters and subgroups observed for S. thermophilus indicate the existence of lineages within the species, an observation which was further supported to a variable degree by the distribution and/or the architecture of several genomic traits. These would include exopolysaccharide (EPS) gene clusters, Clustered Regularly Interspaced Short Palindromic Repeats (CRISPRs)-CRISPR associated (Cas) systems, as well as restriction-modification (R-M) systems and genomic islands (GIs). Of note, the histidine biosynthetic cluster was found present in all cluster A strains (plus strain NCTC12958^T) but was absent from all strains in cluster B. Other loci related to lactose/galactose catabolism and urea metabolism, aminopeptidases, the majority of amino acid and peptide transporters, as well as amino acid biosynthetic pathways were found to be conserved in all strains suggesting their central role for the species. Our study highlights the necessity of sequencing and analyzing more S. thermophilus complete genomes to further elucidate important aspects of strain diversity within this starter culture that may be related to its application in the dairy industry.

Phylogeography of the second plague pandemic revealed through analysis of historical Yersinia pestis genomes

The second plague pandemic, caused by Yersinia pestis, devastated Europe and the nearby regions between the 14th and 18th centuries AD. Here we analyse human remains from ten European archaeological sites spanning this period and reconstruct 34 ancient Y. pestis genomes. Our data support an initial entry of the bacterium through eastern Europe, the absence of genetic diversity during the Black Death, and low within-outbreak diversity thereafter. Analysis of post-Black Death genomes shows the diversification of a Y. pestis lineage into multiple genetically distinct clades that may have given rise to more than one disease reservoir in, or close to, Europe. In addition, we show the loss of a genomic region that includes virulence-related genes in strains associated with late stages of the pandemic. The deletion was also identified in genomes connected with the first plague pandemic (541-750 AD), suggesting a comparable evolutionary trajectory of Y. pestis during both events.

Evolutionary mechanism leading to the multi-cagA genotype in Helicobacter pylori

Infection with CagA+ Helicobacter pylori strains is linked to an increased risk for gastric diseases, including gastric cancer. Recent evidence indicates that dynamic expansion and contraction of cagA copy number may serve as a novel mechanism to enhance disease development. Herein, comparative genomic analysis divided hpEurope into two groups: hpEurope/type-A and type-B. Only hpEurope/type-B displayed the multi-cagA genotype. Further analysis showed that cagPAI appears to have been independently introduced into two different H. pylori types, termed pre-type-A and pre-type-B, which consequently evolved to cagPAI type-A and type-B, respectively; importantly, all multi-cagA genotype strains displayed cagPAI type-B. Two direct cagA-flanking repeats of a genetic element termed CHA-ud were essential for the multi-cagA genotype in strain PMSS1 (hpEurope/type-B and cagPAI type-B). Furthermore, introduction of this genetic element into strain G27 (hpEurope/type-A and cagPAI type-A) was sufficient to generate the multi-cagA genotype. The critical steps in the evolution of the multi-cagA genotype involved creation of CHA-ud at cagA upstream in cagPAI type-B strains followed by its duplication to cagA downstream. En masse, elucidation of the mechanism by which H. pylori evolved to carry multiple copies of cagA helps to provide a better understanding of how this ancient pathogen interacts with its host.