Chromosomal rearrangements in enteric bacteria

Whole genome sequence analysis of Salmonella Typhi in Papua New Guinea reveals an established population of genotype 2.1.7 sensitive to antimicrobials

BACKGROUND

Typhoid fever, a systemic infection caused by Salmonella enterica serovar Typhi, remains a considerable public health threat in impoverished regions within many low- and middle-income settings. However, we still lack a detailed understanding of the emergence, population structure, molecular mechanisms of antimicrobial resistance (AMR), and transmission dynamics of S. Typhi across many settings, particularly throughout the Asia-Pacific islands. Here we present a comprehensive whole genome sequence (WGS) based overview of S. Typhi populations circulating in Papua New Guinea (PNG) over 30 years.

PRINCIPLE FINDINGS

Bioinformatic analysis of 86 S. Typhi isolates collected between 1980-2010 demonstrated that the population structure of PNG is dominated by a single genotype (2.1.7) that appears to have emerged in the Indonesian archipelago in the mid-twentieth century with minimal evidence of inter-country transmission. Genotypic and phenotypic data demonstrated that the PNG S. Typhi population appears to be susceptible to former first line drugs for treating typhoid fever (chloramphenicol, ampicillin and co-trimoxazole), as well as fluoroquinolones, third generation cephalosporins, and macrolides. PNG genotype 2.1.7 was genetically conserved, with very few deletions, and no evidence of plasmid or prophage acquisition. Genetic variation among this population was attributed to either single point mutations, or homologous recombination adjacent to repetitive ribosomal RNA operons.

SIGNIFICANCE

Antimicrobials remain an effective option for the treatment of typhoid fever in PNG, along with other intervention strategies including improvements to water, sanitation and hygiene (WaSH) related infrastructure and potentially the introduction of Vi-conjugate vaccines. However, continued genomic surveillance is warranted to monitor for the emergence of AMR within local populations, or the introduction of AMR associated genotypes of S. Typhi in this setting.

Chromosomal Location Determines the Rate of Intrachromosomal Homologous Recombination in Salmonella

Homologous recombination is an important mechanism directly involved in the repair, organization, and evolution of prokaryotic and eukaryotic chromosomes. We developed a system, based on two genetic cassettes, that allows the measurement of recombinational repair rates between different locations on the chromosome. Using this system, we analyzed 81 different positional combinations throughout the chromosome to answer the question of how the position and orientation of sequences affect intrachromosomal homologous recombination. Our results show that recombination was possible between any two locations tested in this study and that recombinational repair rates varied by just above an order of magnitude. The observed differences in rate do not correlate with distance between the recombination cassettes or with distance from the origin of replication but could be explained if each location contributes individually to the recombination event. The relative levels of accessibility for recombination vary 5-fold between the various cassette locations, and we found that the nucleoid structure of the chromosome may be the major factor influencing the recombinational accessibility of each chromosomal site. Furthermore, we found that the orientation of the recombination cassettes had a significant impact on recombination. Recombinational repair rates for the cassettes inserted as direct repeats are, on average, 2.2-fold higher than those for the same sets inserted as inverted repeats. These results suggest that the bacterial chromosome is not homogenous with regard to homologous recombination, with regions that are more or less accessible, and that the orientation of genes affects recombination rates. IMPORTANCE Bacterial chromosomes frequently carry multiple copies of genes at separate chromosomal locations. In Salmonella, these include the 7 rrn operons and the duplicate tuf genes. Genes within these families coevolve by homologous recombination, but it is not obvious whether their rates of recombination reflect general rates of intrachromosomal recombination or are an evolved property particularly associated with these conserved genes and locations. Using a novel experimental system, we show that recombination is possible between all tested pairs of locations at rates that vary by just above 1 order of magnitude. Differences in rate do not correlate with distance between the sites or distance to the origin of replication but may be explained if each location contributes individually to the recombination event. Our results suggest the existence of bacterial chromosomal domains that are differentially available for recombination and that gene orientation affects recombination rates.

socru: typing of genome-level order and orientation around ribosomal operons in bacteria

Rearrangements of large genome fragments occur in bacteria between repeat sequences and can impact on growth and gene expression. Homologous recombination resulting in inversion between indirect repeats and excision/translocation between direct repeats enables these structural changes. One form of rearrangement occurs around ribosomal operons, found in multiple copies across many bacteria, but identification of these rearrangements by sequencing requires reads of several thousand bases to span the ribosomal operons. With long-read sequencing aiding the routine generation of complete bacterial assemblies, we have developed socru, a typing method for the order and orientation of genome fragments between ribosomal operons. It allows for a single identifier to convey the order and orientation of genome-level structure and we have successfully applied this typing to 433 of the most common bacterial species. In a focused analysis, we observed the presence of multiple structural genotypes in nine bacterial pathogens, underscoring the importance of routinely assessing this form of variation alongside traditional single-nucleotide polymorphism (SNP) typing.

Synchronous termination of replication of the two chromosomes is an evolutionary selected feature in Vibrionaceae

Vibrio cholerae, the causative agent of the cholera disease, is commonly used as a model organism for the study of bacteria with multipartite genomes. Its two chromosomes of different sizes initiate their DNA replication at distinct time points in the cell cycle and terminate in synchrony. In this study, the time-delayed start of Chr2 was verified in a synchronized cell population. This replication pattern suggests two possible regulation mechanisms for other Vibrio species with different sized secondary chromosomes: Either all Chr2 start DNA replication with a fixed delay after Chr1 initiation, or the timepoint at which Chr2 initiates varies such that termination of chromosomal replication occurs in synchrony. We investigated these two models and revealed that the two chromosomes of various Vibrionaceae species terminate in synchrony while Chr2-initiation timing relative to Chr1 is variable. Moreover, the sequence and function of the Chr2-triggering crtS site recently discovered in V. cholerae were found to be conserved, explaining the observed timing mechanism. Our results suggest that it is beneficial for bacterial cells with multiple chromosomes to synchronize their replication termination, potentially to optimize chromosome related processes as dimer resolution or segregation.

Spatial and Temporal Control of Evolution through Replication-Transcription Conflicts

Evolution could potentially be accelerated if an organism could selectively increase the mutation rate of specific genes that are actively under positive selection. Recently, a mechanism that cells can use to target rapid evolution to specific genes was discovered. This mechanism is driven by gene orientation-dependent encounters between DNA replication and transcription machineries. These encounters increase mutagenesis in lagging-strand genes, where replication-transcription conflicts are severe. Due to the orientation and transcription-dependent nature of this process, conflict-driven mutagenesis can be used by cells to spatially (gene-specifically) and temporally (only upon transcription induction) regulate the rate of gene evolution. Here, I summarize recent findings on this topic, and discuss the implications of increasing mutagenesis rates and accelerating evolution through active mechanisms.

Genetic characterization of atypical Citrobacter freundii

The ability of a bacterial population to survive in different niches, as well as in stressful and rapidly changing environmental conditions, depends greatly on its genetic content. To survive such fluctuating conditions, bacteria have evolved different mechanisms to modulate phenotypic variations and related strategies to produce high levels of genetic diversity. Laboratories working in microbiological diagnosis have shown that Citrobacter freundii is very versatile in its colony morphology, as well as in its biochemical, antigenic and pathogenic behaviours. This phenotypic versatility has made C. freundii difficult to identify and it is frequently confused with both Salmonella enterica and Escherichia coli. In order to determine the genomic events and to explain the mechanisms involved in this plasticity, six C. freundii isolates were selected from a phenotypic variation study. An I-CeuI genomic cleavage map was created and eight housekeeping genes, including 16S rRNA, were sequenced. In general, the results showed a range of both phenotypes and genotypes among the isolates with some revealing a greater similarity to C. freundii and some to S. enterica, while others were identified as phenotypic and genotypic intermediary states between the two species. The occurrence of these events in natural populations may have important implications for genomic diversification in bacterial evolution, especially when considering bacterial species boundaries. In addition, such events may have a profound impact on medical science in terms of treatment, course and outcomes of infectious diseases, evading the immune response, and understanding host-pathogen interactions.

Homologous Recombination-Experimental Systems, Analysis, and Significance

Homologous recombination is the most complex of all recombination events that shape genomes and produce material for evolution. Homologous recombination events are exchanges between DNA molecules in the lengthy regions of shared identity, catalyzed by a group of dedicated enzymes. There is a variety of experimental systems in Escherichia coli and Salmonella to detect homologous recombination events of several different kinds. Genetic analysis of homologous recombination reveals three separate phases of this process: pre-synapsis (the early phase), synapsis (homologous strand exchange), and post-synapsis (the late phase). In E. coli, there are at least two independent pathway of the early phase and at least two independent pathways of the late phase. All this complexity is incongruent with the originally ascribed role of homologous recombination as accelerator of genome evolution: there is simply not enough duplication and repetition in enterobacterial genomes for homologous recombination to have a detectable evolutionary role and therefore not enough selection to maintain such a complexity. At the same time, the mechanisms of homologous recombination are uniquely suited for repair of complex DNA lesions called chromosomal lesions. In fact, the two major classes of chromosomal lesions are recognized and processed by the two individual pathways at the early phase of homologous recombination. It follows, therefore, that homologous recombination events are occasional reflections of the continual recombinational repair, made possible in cases of natural or artificial genome redundancy.

Citrobacter rodentium is an unstable pathogen showing evidence of significant genomic flux

Citrobacter rodentium is a natural mouse pathogen that causes attaching and effacing (A/E) lesions. It shares a common virulence strategy with the clinically significant human A/E pathogens enteropathogenic E. coli (EPEC) and enterohaemorrhagic E. coli (EHEC) and is widely used to model this route of pathogenesis. We previously reported the complete genome sequence of C. rodentium ICC168, where we found that the genome displayed many characteristics of a newly evolved pathogen. In this study, through PFGE, sequencing of isolates showing variation, whole genome transcriptome analysis and examination of the mobile genetic elements, we found that, consistent with our previous hypothesis, the genome of C. rodentium is unstable as a result of repeat-mediated, large-scale genome recombination and because of active transposition of mobile genetic elements such as the prophages. We sequenced an additional C. rodentium strain, EX-33, to reveal that the reference strain ICC168 is representative of the species and that most of the inactivating mutations were common to both isolates and likely to have occurred early on in the evolution of this pathogen. We draw parallels with the evolution of other bacterial pathogens and conclude that C. rodentium is a recently evolved pathogen that may have emerged alongside the development of inbred mice as a model for human disease.

Selection for Chromosome Architecture in Bacteria

Bacterial chromosomes are immense polymers whose faithful replication and segregation are crucial to cell survival. The ability of proteins such as FtsK to move unidirectionally toward the replication terminus, and direct DNA translocation into the appropriate daughter cell during cell division, requires that bacterial genomes maintain an architecture for the orderly replication and segregation of chromosomes. We suggest that proteins that locate the replication terminus exploit strand-biased sequences that are overrepresented on one DNA strand, and that selection increases with decreased distance to the replication terminus. We report a generalized method for detecting these architecture imparting sequences (AIMS) and have identified AIMS in nearly all bacterial genomes. Their increased abundance on leading strands and decreased abundance on lagging strands toward replication termini are not the result of changes in mutational bias; rather, they reflect a gradient of long-term positive selection for AIMS. The maintenance of the pattern of AIMS across the genomes of related bacteria independent of their positions within individual genes suggests a well-conserved role in genome biology. The stable gradient of AIMS abundance from replication origin to terminus suggests that the replicore acts as a target of selection, where selection for chromosome architecture results in the maintenance of gene order and in the lack of high-frequency DNA inversion within replicores.

Genomic diversity of Erwinia carotovora subsp. carotovora and its correlation with virulence

We used genetic and biochemical methods to examine the genomic diversity of the enterobacterial plant pathogen Erwinia carotovora subsp. carotovora. The results obtained with each method showed that E. carotovora subsp. carotovora strains isolated from one ecological niche, potato plants, are surprisingly diverse compared to related pathogens. A comparison of 23 partial mdh sequences revealed a maximum pairwise difference of 10.49% and an average pairwise difference of 2.13%, values which are much greater than the maximum variation (1.81%) and average variation (0.75%) previously reported for Escherichia coli. Pulsed-field gel electrophoresis analysis of I-CeuI-digested genomic DNA revealed seven rrn operons in all E. carotovora subsp. carotovora strains examined except strain WPP17, which had only six copies. We identified 26 I-CeuI restriction fragment length polymorphism patterns and observed significant polymorphism in fragment sizes ranging from 100 to 450 kb for all strains. We detected large plasmids in two strains, including the model strain E. carotovora subsp. carotovora 71. The two least virulent strains had an unusual chromosomal structure, suggesting that a particular pulsotype is correlated with virulence. To compare chromosomal organization of multiple enterobacterial genomes, several genes were mapped onto I-CeuI fragments. We identified portions of the genome that appear to be conserved across enterobacteria and portions that have undergone genome rearrangements. We found that the least virulent strain, WPP17, failed to oxidize cellobiose and was missing several hrp and hrc genes. The unexpected variability among isolates obtained from clonal hosts in one region and in one season suggests that factors other than the host plant, potato, drive the evolution of this common environmental bacterium and key plant pathogen.

Comparative genomics of Salmonella enterica serovar Typhi strains Ty2 and CT18

We present the 4.8-Mb complete genome sequence of Salmonella enterica serovar Typhi strain Ty2, a human-specific pathogen causing typhoid fever. A comparison with the genome sequence of recently isolated S. enterica serovar Typhi strain CT18 showed that 29 of the 4,646 predicted genes in Ty2 are unique to this strain, while 84 genes are unique to CT18. Both genomes contain more than 200 pseudogenes; 9 of these genes in CT18 are intact in Ty2, while 11 intact CT18 genes are pseudogenes in Ty2. A half-genome interreplichore inversion in Ty2 relative to CT18 was confirmed. The two strains exhibit differences in prophages, insertion sequences, and island structures. While CT18 carries two plasmids, one conferring multiple drug resistance, Ty2 has no plasmids and is sensitive to antibiotics.

Unlocking the genome of the human typhoid bacillus

Molecular studies are shedding new light on the pathogenesis of human typhoid fever, which is still a very common disease in developing countries. For example, the total genome DNA sequence has recently been determined for a multiple-drug-resistant Salmonella typhi, the serotype that is the cause of typhoid fever. The genome sequence showed many distinguishing features, including clusters of S typhi specific genes and a large number--over 200--of pseudogenes. This information, together with other molecular studies, has provided vital clues in several important areas of typhoid biology. We have new insights into the mechanisms underpinning the human host specificity of S typhi, and have exploitable new routes to improved diagnostics and a better understanding of the epidemiology of the disease.

Physical and genetic map of Enterococcus faecium ATCC19434 and demonstration of intra- and interspecific genomic diversity in enterococci

A physical map of the Enterococcus faecium ATCC19434 chromosome was constructed by NotI, I-CeuI and Sse8387I. The chromosome was a circular DNA of 2600 kb in size, and contained six rRNA operons (rrn). The locations and orientations of the six rrn operons and 24 different determinants were mapped. Genomes of three additional E. faecium strains were also analyzed by I-CeuI digestion, and the genome sizes were found to vary from 2550 to 2995 kb. We further investigated the genome sizes and number of rrn operons in four E. faecalis, one E. avium, and one E. durans strains. The genome sizes were larger than E. faecium: 3000-3250 kb in E. faecalis, 3445 kb in E. avium, and 3070 kb in E. durans. E. avium and E. durans contained six rrn operons as in E. faecium, but all the E. faecalis strains possessed four rrn operons.

Complete genome sequence of a multiple drug resistant Salmonella enterica serovar Typhi CT18

Salmonella enterica serovar Typhi (S. typhi) is the aetiological agent of typhoid fever, a serious invasive bacterial disease of humans with an annual global burden of approximately 16 million cases, leading to 600,000 fatalities. Many S. enterica serovars actively invade the mucosal surface of the intestine but are normally contained in healthy individuals by the local immune defence mechanisms. However, S. typhi has evolved the ability to spread to the deeper tissues of humans, including liver, spleen and bone marrow. Here we have sequenced the 4,809,037-base pair (bp) genome of a S. typhi (CT18) that is resistant to multiple drugs, revealing the presence of hundreds of insertions and deletions compared with the Escherichia coli genome, ranging in size from single genes to large islands. Notably, the genome sequence identifies over two hundred pseudogenes, several corresponding to genes that are known to contribute to virulence in Salmonella typhimurium. This genetic degradation may contribute to the human-restricted host range for S. typhi. CT18 harbours a 218,150-bp multiple-drug-resistance incH1 plasmid (pHCM1), and a 106,516-bp cryptic plasmid (pHCM2), which shows recent common ancestry with a virulence plasmid of Yersinia pestis.

Genotypic analyses of Escherichia coli O157:H7 and O157 nonmotile isolates recovered from beef cattle and carcasses at processing plants in the Midwestern states of the United States

Escherichia coli O157:H7 and O157 nonmotile isolates (E. coli O157) previously were recovered from feces, hides, and carcasses at four large Midwestern beef processing plants (R. O. Elder, J. E. Keen, G. R. Siragusa, G. A. Barkocy-Gallagher, M. Koohmaraie, and W. W. Laegreid, Proc. Natl. Acad. Sci. USA 97:2999-3003, 2000). The study implied relationships between cattle infection and carcass contamination within single-source lots as well as between preevisceration and postprocessing carcass contamination, based on prevalence. These relationships now have been verified based on identification of isolates by genomic fingerprinting. E. coli O157 isolates from all positive samples were analyzed by pulsed-field gel electrophoresis of genomic DNA after digestion with XbaI. Seventy-seven individual subtypes (fingerprint patterns) grouping into 47 types were discerned among 343 isolates. Comparison of the fingerprint patterns revealed three clusters of isolates, two of which were closely related to each other. Remarkably, isolates carrying both Shiga toxin genes and nonmotile isolates largely fell into specific clusters. Within lots analyzed, 68.2% of the postharvest (carcass) isolates matched preharvest (animal) isolates. For individual carcasses, 65.3 and 66.7% of the isolates recovered postevisceration and in the cooler, respectively, matched those recovered preevisceration. Multiple isolates were analyzed from some carcass samples and were found to include strains with different genotypes. This study suggests that most E. coli O157 carcass contamination originates from animals within the same lot and not from cross-contamination between lots. In addition, the data demonstrate that most carcass contamination occurs very early during processing.

Evolution of prokaryotic gene order: genome rearrangements in closely related species

Conservation of gene order in prokaryotes has become important in predicting protein function because, over the evolutionary timescale, genomes are shuffled so that local gene-order conservation reflects the functional constraints within the protein. Here, we compare closely related genomes to identify the rate with which gene order is disrupted and to infer the genes involved in the genome rearrangement.

Tracing of Salmonella spp. in two pork slaughter and cutting plants using serotyping and macrorestriction genotyping

AIMS

The origin of Salmonella contamination of pork products is not well established. In order to further this knowledge, the transmission of Salmonella spp. from live pigs to pork cuts was investigated in two pork slaughter and cutting plants.

METHODS AND RESULTS

Salmonella spp. were isolated from both pork (pigs, carcasses, cuts) and the environment before and during slaughterhouse activities. Eight serotypes were identified. XbaI and SpeI macrorestriction distinguished 20 genotypes of Salmonella Typhimurium and 16 genotypes of Salmonella Derby. A major cluster of Salmonella Typhimurium genotypes was common to both plants and all pig-related genotypes, while a predominant pig-related Salmonella Derby genotype was common to both plants.

CONCLUSION

None of the Salmonella strains persisted for long periods in the pork-processing environments.

SIGNIFICANCE AND IMPACT OF THE STUDY

This work shows that contaminated live pigs, because of bacterial spread due to the process and ineffective cleaning procedures, are involved in Salmonella contamination.

Evaluating genome dynamics: the constraints on rearrangements within bacterial genomes

Inversions and translocations distinguish the genomes of closely related bacterial species, but most of these rearrangements preserve the relationship between the rearranged fragments and the axis of chromosome replication. Within species, such rearrangements are found less frequently, except in the case of clinical isolates of human pathogens, where rearrangements are very frequent.

Genome plasticity in Enterobacteriaceae

The comparative analysis of multiple representatives of the genomes of particular species are leading us away from a view of bacterial genomes as static, monolithic structures towards the view that they are relatively variable, fluid structures. This plasticity is mainly the result of the rearrangement of genes within the genome and the acquisition of novel genes by horizontal transfer systems, e. g. plasmids, bacteriophages, transposons or gene cassettes. These mechanisms often act in concert thus generating a complex genetic structure. Genomic variations are not a phenomenon at the DNA level alone, they influence the phenotype of a bacterium as well and can render a formerly harmless organism into a hazardous pathogen. This review deals not only with the mechanisms of genome rearrangements and the horizontal transfer of genes in Enterobacteriaceae but also points out that mobile genetic elements themselves are subjected to variation.

Rates and patterns of chromosome evolution in enteric bacteria

Although several types of large-scale alterations potentially affect the structure and organization of bacterial genomes, recent analyses of physical maps and complete genomic sequences reveal that chromosome heterogeneity in enteric bacteria has resulted from the acquisition and deletion of large segments of DNA. These acquired sequences can provide novel functions immediately upon their introduction and play a significant role in the diversification of bacterial species.