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.     

Recombination and population structure in Salmonella enterica

Salmonella enterica is a bacterial pathogen that causes enteric fever and gastroenteritis in humans and animals. Although its population structure was long described as clonal, based on high linkage disequilibrium between loci typed by enzyme electrophoresis, recent examination of gene sequences has revealed that recombination plays an important evolutionary role. We sequenced around 10% of the core genome of 114 isolates of enterica using a resequencing microarray. Application of two different analysis methods (Structure and ClonalFrame) to our genomic data allowed us to define five clear lineages within S. enterica subspecies enterica, one of which is five times older than the other four and two thirds of the age of the whole subspecies. We show that some of these lineages display more evidence of recombination than others. We also demonstrate that some level of sexual isolation exists between the lineages, so that recombination has occurred predominantly between members of the same lineage. This pattern of recombination is compatible with expectations from the previously described ecological structuring of the enterica population as well as mechanistic barriers to recombination observed in laboratory experiments. In spite of their relatively low level of genetic differentiation, these lineages might therefore represent incipient species.     

Long-term prevalence of antimicrobial-resistant Salmonella enterica subspecies enterica Serovar infantis in the broiler chicken industry in Japan

Eleven broiler isolates of Salmonella Infantis obtained between 1989 and 1998 were examined for antimicrobial susceptibility and pulse field gel electrophoresis (PFGE) profiles. Seven strains of S. Infantis isolated after 1993 harbored similar antimicrobial susceptibilities to the recent isolates between 2001 and 2003. In comparison of PFGE profile with 22 isolates obtained from 22 apparently healthy broiler chickens between 2001 and 2003, the predominant cluster included the seven strains isolated after 1993. We could not clarify the reasons why the serovar has been prevalent in the broiler industry for a long time, but current antimicrobial usage is not always linked to its prevalence.     

Comparative genomics of 28 Salmonella enterica isolates: evidence for CRISPR-mediated adaptive sublineage evolution

Despite extensive surveillance, food-borne Salmonella enterica infections continue to be a significant burden on public health systems worldwide. As the S. enterica species comprises sublineages that differ greatly in antigenic representation, virulence, and antimicrobial resistance phenotypes, a better understanding of the species' evolution is critical for the prediction and prevention of future outbreaks. The roles that virulence and resistance phenotype acquisition, exchange, and loss play in the evolution of S. enterica sublineages, which to a certain extent are represented by serotypes, remains mostly uncharacterized. Here, we compare 17 newly sequenced and phenotypically characterized nontyphoidal S. enterica strains to 11 previously sequenced S. enterica genomes to carry out the most comprehensive comparative analysis of this species so far. These phenotypic and genotypic data comparisons in the phylogenetic species context suggest that the evolution of known S. enterica sublineages is mediated mostly by two mechanisms, (i) the loss of coding sequences with known metabolic functions, which leads to functional reduction, and (ii) the acquisition of horizontally transferred phage and plasmid DNA, which provides virulence and resistance functions and leads to increasing specialization. Matches between S. enterica clustered regularly interspaced short palindromic repeats (CRISPR), part of a defense mechanism against invading plasmid and phage DNA, and plasmid and prophage regions suggest that CRISPR-mediated immunity could control short-term phenotype changes and mediate long-term sublineage evolution. CRISPR analysis could therefore be critical in assessing the evolutionary potential of S. enterica sublineages and aid in the prediction and prevention of future S. enterica outbreaks.     

Comparative genomics reveals genes significantly associated with woody hosts in the plant pathogen Pseudomonas syringae

The diversification of lineages within Pseudomonas syringae has involved a number of adaptive shifts from herbaceous hosts onto various species of tree, resulting in the emergence of highly destructive diseases such as bacterial canker of kiwi and bleeding canker of horse chestnut. This diversification has involved a high level of gene gain and loss, and these processes are likely to play major roles in the adaptation of individual lineages onto their host plants. In order to better understand the evolution of P. syringae onto woody plants, we have generated de novo genome sequences for 26 strains from the P. syringae species complex that are pathogenic on a range of woody species, and have looked for statistically significant associations between gene presence and host type (i.e. woody or herbaceous) across a phylogeny of 64 strains. We have found evidence for a common set of genes associated with strains that are able to colonize woody plants, suggesting that divergent lineages have acquired similarities in genome composition that may form the genetic basis of their adaptation to woody hosts. We also describe in detail the gain, loss and rearrangement of specific loci that may be functionally important in facilitating this adaptive shift. Overall, our analyses allow for a greater understanding of how gene gain and loss may contribute to adaptation in P. syringae.     

Macrophage apoptosis associated with Salmonella enterica serovar Typhi plasmid

The present study was undertaken to investigate the relationship between plasmid isolated from S. enterica serovar Typhi (pR(ST98)) and macrophage apoptosis. pR(ST98) was transferred into an attenuated S. enterica serovar Typhimurium strain RIA to create a transconjugant pRsT98/RIA. Standard S. enterica serovar Typhimurium virulence strain SR-11 was used as a positive control, and RIA as a negative one. Murine macrophage-like cell line (J774A.1) was used as an infectious cell model in vitro. In order to determine the inhibition and bactericidal effect of amikacin (AMK) to extracellular bacteria and the best optimization co-culture ratio between Salmonella and J774A.1, the minimal inhibitory concentration (MIC) and minimal bactericidal concentration (MBC) of AMK to strains SR-11, pR(ST98)/RIA and RIA and multiplicity of infection (MOI) were detected first, and then J774A.1 was infected by the above three serovar Typhimurium strains. Apoptosis of J774A.1 was examined with electron microscopy and flow cytometry after annexin-V/propidium iodide labeling at 0, 1, 3, 6, 12 and 24 h. Mitochondrial membrane potential was detected by JC-1 staining method. It was demonstrated that MIC of AMK to the three strains was 10 microg/ml, MBC was 80 microg/ml, and optimal MOI was 100:1. pR(ST98)/RIA resulted in a higher apoptosis of J774A.1 than RIA, apoptotic features such as chromatin margination could be observed after 3 h, and death of J774A.1 cells was associated with the loss of mitochondrial membrane potential. These results indicated that pR(ST98) could enhance the virulence of its host bacteria, evidenced by increased macrophage apoptosis.     

Comparative population genomics identified genomic regions and candidate genes associated with fruit domestication traits in peach

Crop evolution is a long‐term process involving selection by natural evolutionary forces and anthropogenic influences; however, the genetic mechanisms underlying the domestication and improvement of fruit crops have not been well studied to date. Here, we performed a population structure analysis in peach (Prunus persica) based on the genome‐wide resequencing of 418 accessions and confirmed the presence of an obvious domestication event during evolution. We identified 132 and 106 selective sweeps associated with domestication and improvement, respectively. Analysis of their tissue‐specific expression patterns indicated that the up‐regulation of selection genes during domestication occurred mostly in fruit and seeds as opposed to other organs. However, during the improvement stage, more up‐regulated selection genes were identified in leaves and seeds than in the other organs. Genome‐wide association studies (GWAS) using 4.24 million single nucleotide polymorphisms (SNPs) revealed 171 loci associated with 26 fruit domestication traits. Among these loci, three candidate genes were highly associated with fruit weight and the sorbitol and catechin content in fruit. We demonstrated that as the allele frequency of the SNPs associated with high polyphenol composition decreased during peach evolution, alleles associated with high sugar content increased significantly. This indicates that there is genetic potential for the breeding of more nutritious fruit with enhanced bioactive polyphenols without disturbing a harmonious sugar and acid balance by crossing with wild species. This study also describes the development of the genomic resources necessary for evolutionary research in peach and provides the large‐scale characterization of key agronomic traits in this crop species.     

Comparative genomics of Helicobacter pylori strains of China associated with different clinical outcome

In this study, a whole-genome CombiMatrix Custom oligonucleotide tiling microarray with 90,000 probes covering six sequenced Helicobacter pylori (H. pylori) genomes was designed. This microarray was used to compare the genomic profiles of eight unsequenced strains isolated from patients with different gastroduodenal diseases in Heilongjiang province of China. Since significant genomic variation was found among these strains, an additional 76 H. pylori strains associated with different clinical outcomes were isolated from various provinces of China. These strains were tested by polymerase chain reaction to demonstrate this distinction. We identified several highly variable regions in strains associated with gastritis, gastric ulceration, and gastric cancer. These regions are associated with genes involved in the bacterial type I, type II, and type III R-M systems. They were also associated with the virB gene, which lies on the well-studied cag pathogenic island. While previous studies have reported on the diverse genetic characterization of this pathogenic island, in this study, we find that it is conserved in all strains tested by microarray. Moreover, a number of genes involved in the type IV secretion system, which is related to horizontal DNA transfer between H. pylori strains, were identified in the comparative analysis of the strain-specific genes. These findings may provide insight into new biomarkers for the prediction of gastric diseases.     

Evolution and population structure of Salmonella enterica serovar Newport

Salmonellosis caused by Salmonella enterica serovar Newport is a major global public health concern, particularly because S. Newport isolates that are resistant to multiple drugs (MDR), including third-generation cephalosporins (MDR-AmpC phenotype), have been commonly isolated from food animals. We analyzed 384 S. Newport isolates from various sources by a multilocus sequence typing (MLST) scheme to study the evolution and population structure of the serovar. These were compared to the population structure of S. enterica serovars Enteritidis, Kentucky, Paratyphi B, and Typhimurium. Our S. Newport collection fell into three lineages, Newport-I, Newport-II, and Newport-III, each of which contained multiple sequence types (STs). Newport-I has only a few STs, unlike Newport-II or Newport-III, and has possibly emerged recently. Newport-I is more prevalent among humans in Europe than in North America, whereas Newport-II is preferentially associated with animals. Two STs of Newport-II encompassed all MDR-AmpC isolates, suggesting recent global spread after the acquisition of the bla(CMY-2) gene. In contrast, most Newport-III isolates were from humans in North America and were pansusceptible to antibiotics. Newport was intermediate in population structure to the other serovars, which varied from a single monophyletic lineage in S. Enteritidis or S. Typhimurium to four discrete lineages within S. Paratyphi B. Both mutation and homologous recombination are responsible for diversification within each of these lineages, but the relative frequencies differed with the lineage. We conclude that serovars of S. enterica provide a variety of different population structures.     

Comparative genomics reveals unusually long motifs in mammalian genomes

MOTIVATION: The recent discovery of the first small modulatory RNA (smRNA) presents the challenge of finding other molecules of similar length and conservation level. Unlike short interfering RNA (siRNA) and micro-RNA (miRNA), effective computational and experimental screening methods are not currently known for this species of RNA molecule, and the discovery of the one known example was partly fortuitous because it happened to be complementary to a well-studied DNA binding motif (the Neuron Restrictive Silencer Element). RESULTS: The existing comparative genomics approaches (e.g., phylogenetic footprinting) rely on alignments of orthologous regions across multiple genomes. This approach, while extremely valuable, is not suitable for finding motifs with highly diverged "non-alignable" flanking regions. Here we show that several unusually long and well conserved motifs can be discovered de novo through a comparative genomics approach that does not require an alignment of orthologous upstream regions. These motifs, including Neuron Restrictive Silencer Element, were missed in recent comparative genomics studies that rely on phylogenetic footprinting. While the functions of these motifs remain unknown, we argue that some may represent biologically important sites. AVAILABILITY: Our comparative genomics software, a web-accessible database of our results and a compilation of experimentally validated binding sites for NRSE can be found at http://www.cse.ucsd.edu/groups/bioinformatics.     

Genomic population structure associated with repeated escape of Salmonella enterica ATCC14028s from the laboratory into nature

Salmonella enterica serovar Typhimurium strain ATCC14028s is commercially available from multiple national type culture collections, and has been widely used since 1960 for quality control of growth media and experiments on fitness (“laboratory evolution”). ATCC14028s has been implicated in multiple cross-contaminations in the laboratory, and has also caused multiple laboratory infections and one known attempt at bioterrorism. According to hierarchical clustering of 3002 core gene sequences, ATCC14028s belongs to HierCC cluster HC20_373 in which most internal branch lengths are only one to three SNPs long. Many natural Typhimurium isolates from humans, domesticated animals and the environment also belong to HC20_373, and their core genomes are almost indistinguishable from those of laboratory strains. These natural isolates have infected humans in Ireland and Taiwan for decades, and are common in the British Isles as well as the Americas. The isolation history of some of the natural isolates confirms the conclusion that they do not represent recent contamination by the laboratory strain, and 10% carry plasmids or bacteriophages which have been acquired in nature by HGT from unrelated bacteria. We propose that ATCC14028s has repeatedly escaped from the laboratory environment into nature via laboratory accidents or infections, but the escaped micro-lineages have only a limited life span. As a result, there is a genetic gap separating HC20_373 from its closest natural relatives due to a divergence between them in the late 19^th century followed by repeated extinction events of escaped HC20_373.     

Comparative genomic analysis of two Burkholderia glumae strains from different geographic origins reveals a high degree of plasticity in genome structure associated with genomic islands

Burkholderia glumae is the major causal agent of bacterial panicle blight of rice, a growing disease problem in global rice production. To better understand its genome-scale characteristics, the genome of the highly virulent B. glumae strain 336gr-1 isolated from Louisiana, USA was sequenced using the Illumina Genome Analyser II system. De novo assembled 336gr-1 contigs were aligned and compared with the previously sequenced genome of B. glumae strain BGR1, which was isolated from an infected rice plant in South Korea. Comparative analysis of the whole genomes of B. glumae 336gr-1 and B. glumae BGR1 revealed numerous unique genomic regions present only in one of the two strains. These unique regions contained accessory genes including mobile elements and phage-related genes, and some of the unique regions in B. glumae BGR1 corresponded to predicted genomic islands. In contrast, little variation was observed in known and potential virulence genes between the two genomes. The considerable amount of plasticity largely based on accessory genes and genome islands observed from the comparison of the genomes of these two strains of B. glumae may explain the versatility of this bacterial species in various environmental conditions and geographic locations.     

Comparative genomics reveals new functional insights in uncultured MAST species

Heterotrophic lineages of stramenopiles exhibit enormous diversity in morphology, lifestyle, and habitat. Among them, the marine stramenopiles (MASTs) represent numerous independent lineages that are only known from environmental sequences retrieved from marine samples. The core energy metabolism characterizing these unicellular eukaryotes is poorly understood. Here, we used single-cell genomics to retrieve, annotate, and compare the genomes of 15 MAST species, obtained by coassembling sequences from 140 individual cells sampled from the marine surface plankton. Functional annotations from their gene repertoires are compatible with all of them being phagocytotic. The unique presence of rhodopsin genes in MAST species, together with their widespread expression in oceanic waters, supports the idea that MASTs may be capable of using sunlight to thrive in the photic ocean. Additional subsets of genes used in phagocytosis, such as proton pumps for vacuole acidification and peptidases for prey digestion, did not reveal particular trends in MAST genomes as compared with nonphagocytotic stramenopiles, except a larger presence and diversity of V-PPase genes. Our analysis reflects the complexity of phagocytosis machinery in microbial eukaryotes, which contrasts with the well-defined set of genes for photosynthesis. These new genomic data provide the essential framework to study ecophysiology of uncultured species and to gain better understanding of the function of rhodopsins and related carotenoids in stramenopiles.Subject terms: Genomics, Microbiology     

Comparative genomics reveals insights into cyanobacterial evolution and habitat adaptation

Cyanobacteria are photosynthetic prokaryotes that inhabit diverse aquatic and terrestrial environments. However, the evolutionary mechanisms involved in the cyanobacterial habitat adaptation remain poorly understood. Here, based on phylogenetic and comparative genomic analyses of 650 cyanobacterial genomes, we investigated the genetic basis of cyanobacterial habitat adaptation (marine, freshwater, and terrestrial). We show: (1) the expansion of gene families is a common strategy whereby terrestrial cyanobacteria cope with fluctuating environments, whereas the genomes of many marine strains have undergone contraction to adapt to nutrient-poor conditions. (2) Hundreds of genes are strongly associated with specific habitats. Genes that are differentially abundant in genomes of marine, freshwater, and terrestrial cyanobacteria were found to be involved in light sensing and absorption, chemotaxis, nutrient transporters, responses to osmotic stress, etc., indicating the importance of these genes in the survival and adaptation of organisms in specific habitats. (3) A substantial fraction of genes that facilitate the adaptation of Cyanobacteria to specific habitats are contributed by horizontal gene transfer, and such genetic exchanges are more frequent in terrestrial cyanobacteria. Collectively, our results further our understandings of the adaptations of Cyanobacteria to different environments, highlighting the importance of ecological constraints imposed by the environment in shaping the evolution of Cyanobacteria.Subject terms: Phylogenetics, Microbial ecology     

Comparative genomics reveals birth and death of fragile regions in mammalian evolution

Background  

An important question in genome evolution is whether there exist fragile regions (rearrangement hotspots) where chromosomal rearrangements are happening over and over again. Although nearly all recent studies supported the existence of fragile regions in mammalian genomes, the most comprehensive phylogenomic study of mammals raised some doubts about their existence.     

Extracellular polysaccharides associated with thin aggregative fimbriae of Salmonella enterica serovar enteritidis

Lipopolysaccharide (LPS) O polysaccharide was identified as the principle factor impeding intercellular formation of intact thin aggregative fimbriae (Tafi) in Salmonella enterica serovar Enteritidis. The extracellular nucleation-precipitation assembly pathway for these organelles was investigated by quantifying fimbrial formation between deltaagfA (AgfA recipient) and deltaagfB (AgfA donor) cells harboring mutations in LPS (galE::Tn10) and/or cellulose (deltabcsA) synthesis. Intercellular complementation could be detected between deltaagfA and deltaagfB strains only when both possessed the galE mutation. LPS O polysaccharide appears to be an impenetrable barrier to AgfA assembly between cells but not within individual cells. The presence of cellulose did not restrict Tafi formation between cells. Transmission electron microscopy of w+ S. enterica serovar Enteritidis 3b cells revealed diffuse Tafi networks without discernible fine structure. In the absence of cellulose, however, individual Tafi fibers were clearly visible, appeared to be occasionally branched, and showed the generally distinctive appearance described for Escherichia coli K-12 curli. A third extracellular matrix component closely associated with cellulose and Tafi was detected on Western blots by using immune serum raised to whole, purified Tafi aggregates. Cellulose was required to tightly link this material to cells. Antigenically similar material was also detected in S. enterica serovar Typhimurium and one diarrheagenic E. coli isolate. Preliminary analysis indicated that this material represented an anionic, extracellular polysaccharide that was distinct from colanic acid. Therefore, Tafi in their native state appear to exist as a complex with cellulose and at least one other component.     

ChIP‐seq and transcriptome analysis of the OmpR regulon of Salmonella enterica serovars Typhi and Typhimurium reveals accessory genes implicated in host colonization

OmpR is a multifunctional DNA binding regulator with orthologues in many enteric bacteria that exhibits classical regulator activity as well as nucleoid‐associated protein‐like characteristics. In the enteric pathogen Salmonella enterica, using chromatin immunoprecipitation of OmpR:FLAG and nucleotide sequencing, 43 putative OmpR binding sites were identified in S. enterica serovar Typhi, 22 of which were associated with OmpR‐regulated genes. Mutation of a sequence motif (TGTWACAW) that was associated with the putative OmpR binding sites abrogated binding of OmpR:6×His to the tviA upstream region. A core set of 31 orthologous genes were found to exhibit OmpR‐dependent expression in both S. Typhi and S. Typhimurium. S. Typhimurium‐encoded orthologues of two divergently transcribed OmpR‐regulated operons (SL1068–71 and SL1066–67) had a putative OmpR binding site in the inter‐operon region in S. Typhi, and were characterized using in vitro and in vivo assays. These operons are widely distributed within S. enterica but absent from the closely related Escherichia coli. SL1066 and SL1067 were required for growth on N‐acetylmuramic acid as a sole carbon source. SL1068–71 exhibited sequence similarity to sialic acid uptake systems and contributed to colonization of the ileum and caecum in the streptomycin‐pretreated mouse model of colitis.