Andre Villegas

Pan-genome sequence analysis using Panseq: an online tool for the rapid analysis of core and accessory genomic regions

BackgroundThe pan-genome of a bacterial species consists of a core and an accessory gene pool. The accessory genome is thought to be an important source of genetic variability in bacterial populations and is gained through lateral gene transfer, allowing subpopulations of bacteria to better adapt to specific niches. Low-cost and high-throughput sequencing platforms have created an exponential increase in genome sequence data and an opportunity to study the pan-genomes of many bacterial species. In this study, we describe a new online pan-genome sequence analysis program, Panseq.ResultsPanseq was used to identify Escherichia coli O157:H7 and E. coli K-12 genomic islands. Within a population of 60 E. coli O157:H7 strains, the existence of 65 accessory genomic regions identified by Panseq analysis was confirmed by PCR. The accessory genome and binary presence/absence data, and core genome and single nucleotide polymorphisms (SNPs) of six L. monocytogenes strains were extracted with Panseq and hierarchically clustered and visualized. The nucleotide core and binary accessory data were also used to construct maximum parsimony (MP) trees, which were compared to the MP tree generated by multi-locus sequence typing (MLST). The topology of the accessory and core trees was identical but differed from the tree produced using seven MLST loci. The Loci Selector module found the most variable and discriminatory combinations of four loci within a 100 loci set among 10 strains in 1 s, compared to the 449 s required to exhaustively search for all possible combinations; it also found the most discriminatory 20 loci from a 96 loci E. coli O157:H7 SNP dataset.ConclusionPanseq determines the core and accessory regions among a collection of genomic sequences based on user-defined parameters. It readily extracts regions unique to a genome or group of genomes, identifies SNPs within shared core genomic regions, constructs files for use in phylogeny programs based on both the presence/absence of accessory regions and SNPs within core regions and produces a graphical overview of the output. Panseq also includes a loci selector that calculates the most variable and discriminatory loci among sets of accessory loci or core gene SNPs.AvailabilityPanseq is freely available online at http://76.70.11.198/panseq. Panseq is written in Perl.

Genome sequence of adherent-invasive Escherichia coli and comparative genomic analysis with other E. coli pathotypes

BackgroundAdherent and invasive Escherichia coli (AIEC) are commonly found in ileal lesions of Crohns Disease (CD) patients, where they adhere to intestinal epithelial cells and invade into and survive in epithelial cells and macrophages, thereby gaining access to a typically restricted host niche. Colonization leads to strong inflammatory responses in the gut suggesting that AIEC could play a role in CD immunopathology. Despite extensive investigation, the genetic determinants accounting for the AIEC phenotype remain poorly defined. To address this, we present the complete genome sequence of an AIEC, revealing the genetic blueprint for this disease-associated E. coli pathotype.ResultsWe sequenced the complete genome of E. coli NRG857c (O83:H1), a clinical isolate of AIEC from the ileum of a Crohns Disease patient. Our sequence data confirmed a phylogenetic linkage between AIEC and extraintestinal pathogenic E. coli causing urinary tract infections and neonatal meningitis. The comparison of the NRG857c AIEC genome with other pathogenic and commensal E. coli allowed for the identification of unique genetic features of the AIEC pathotype, including 41 genomic islands, and unique genes that are found only in strains exhibiting the adherent and invasive phenotype.ConclusionsUp to now, the virulence-like features associated with AIEC are detectable only phenotypically. AIEC genome sequence data will facilitate the identification of genetic determinants implicated in invasion and intracellular growth, as well as enable functional genomic studies of AIEC gene expression during health and disease.

Genomic and Proteomic Characterization of the Broad-Host-Range Salmonella Phage PVP-SE1: Creation of a New Phage Genus

ABSTRACT (Bacterio)phage PVP-SE1, isolated from a German wastewater plant, presents a high potential value as a biocontrol agent and as a diagnostic tool, even compared to the well-studied typing phage Felix 01, due to its broad lytic spectrum against different Salmonella strains. Sequence analysis of its genome (145,964 bp) shows it to be terminally redundant and circularly permuted. Its G+C content, 45.6 mol%, is lower than that of its hosts (50 to 54 mol%). We found a total of 244 open reading frames (ORFs), representing 91.6% of the coding capacity of the genome. Approximately 46% of encoded proteins are unique to this phage, and 22.1% of the proteins could be functionally assigned. This myovirus encodes a large number of tRNAs (n=24), reflecting its lytic capacity and evolution through different hosts. Tandem mass spectrometric analysis using electron spray ionization revealed 25 structural proteins as part of the mature phage particle. The genome sequence was found to share homology with 140 proteins of the Escherichia coli bacteriophage rV5. Both phages are unrelated to any other known virus, which suggests that an “rV5-like virus” genus should be created within the Myoviridae to contain these two phages.

Genome sequence and analysis of a broad-host range lytic bacteriophage that infects the Bacillus cereus group

BackgroundComparatively little information is available on members of the Myoviridae infecting low G+C content, Gram-positive host bacteria of the family Firmicutes. While numerous Bacillus phages have been isolated up till now only very few Bacillus cereus phages have been characterized in detail.ResultsHere we present data on the large, virulent, broad-host-range B. cereus phage vB_BceM_Bc431v3 (Bc431v3). Bc431v3 features a 158,618 bp dsDNA genome, encompassing 239 putative open reading frames (ORFs) and, 20 tRNA genes encoding 17 different amino acids. Since pulsed-field gel electrophoresis indicated that the genome of this phage has a mass of 155-158 kb Bc431v3 DNA appears not to contain long terminal repeats that are found in the genome of Bacillus phage SPO1.ConclusionsBc431v3 displays significant sequence similarity, at the protein level, to B. cereus phage BCP78, Listeria phage A511 and Enterococcus phage ØEF24C and other morphologically related phages infecting Firmicutes such as Staphylococcus phage K and Lactobacillus phage LP65. Based on these data we suggest that Bc431v3 should be included as a member of the Spounavirinae; however, because of all the diverse taxonomical information has been addressed recently, it is difficult to determine the genus. The Bc431v3 phage contains some highly unusual genes such as gp143 encoding putative tRNAHis guanylyltransferase. In addition, it carries some genes that appear to be related to the host sporulation regulators. These are: gp098, which encodes a putative segregation protein related to FstK/SpoIIIE DNA transporters; gp105, a putative segregation protein; gp108, RNA polymerase sigma factor F/B; and, gp109 encoding RNA polymerase sigma factor G.

Genomic, Proteomic and Physiological Characterization of a T5-like Bacteriophage for Control of Shiga Toxin-Producing Escherichia coli O157:H7

Despite multiple control measures, Escherichia coli O157:H7 (STEC O157:H7) continues to be responsible for many food borne outbreaks in North America and elsewhere. Bacteriophage therapy may prove useful for controlling this pathogen in the host, their environment and food. Bacteriophage vB_EcoS_AKFV33 (AKFV33), a T5-like phage of Siphoviridae lysed common phage types of STEC O157:H7 and not non-O157 E. coli. Moreover, STEC O157:H7 isolated from the same feedlot pen from which the phage was obtained, were highly susceptible to AKFV33. Adsorption rate constant and burst size were estimated to be 9.31×10−9 ml/min and 350 PFU/infected cell, respectively. The genome of AKVF33 was 108,853 bp (38.95% G+C), containing 160 open reading frames (ORFs), 22 tRNA genes and 32 strong promoters recognized by host RNA polymerase. Of 12 ORFs without homologues to T5-like phages, 7 predicted novel proteins while others exhibited low identity (<60%) to proteins in the National Centre for Biotechnology Information database. AKVF33 also lacked the L-shaped tail fiber protein typical of T5, but was predicted to have tail fibers comprised of 2 novel proteins with low identity (37–41%) to tail fibers of E. coli phage phiEco32 of Podoviridae, a putative side tail fiber protein of a prophage from E. coli IAI39 and a conserved domain protein of E. coli MS196-1. The receptor-binding tail protein (pb5) shared an overall identify of 29–72% to that of other T5-like phages, with no region coding for more than 6 amino acids in common. Proteomic analysis identified 4 structural proteins corresponding to the capsid, major tail, tail fiber and pore-forming tail tip (pb2). The genome of AKFV33 lacked regions coding for known virulence factors, integration-related proteins or antibiotic resistance determinants. Phage AKFV33 is a unique, highly lytic STEC O157:H7-specific T5-like phage that may have considerable potential as a pre- and post-harvest biocontrol agent.

A Shigella boydii bacteriophage which resembles Salmonella phage ViI

BackgroundLytic bacteriophages have been applied successfully to control the growth of various foodborne pathogens. Sequencing of their genomes is considered as an important preliminary step to ensure their safety prior to food applications.ResultsThe lytic bacteriophage, ΦSboM-AG3, targets the important foodborne pathogen, Shigella. It is morphologically similar to phage ViI of Salmonella enterica serovar Typhi and a series of phages of Acinetobacter calcoaceticus and Rhizobium meliloti. The complete genome of ΦSboM-AG3 was determined to be 158 kb and was terminally redundant and circularly permuted. Two hundred and sixteen open reading frames (ORFs) were identified and annotated, most of which displayed homology to proteins of Salmonella phage ViI. The genome also included four genes specifying tRNAs.ConclusionsThis is the first time that a Vi-specific phage for Shigella has been described. There is no evidence for the presence of virulence and lysogeny-associated genes. In conclusion, the genome analysis of ΦSboM-AG3 indicates that this phage can be safely used for biocontrol purposes.

An analysis of initiation codon utilization in the Domain Bacteria - concerns about the quality of bacterial genome annotation.

Using custom software (Inidon) we have examined the initiation codon utilization in 620 complete bacterial genomes downloaded from the National Center for Biotechnology Information (NCBI). The mean utilization of ATG, GTG and TTG codons is 80.1, 11.6 and 7.8 %, respectively. In most cases in which similar species or strains have been analysed the utilization percentages of the three initiation codons are remarkably similar, but in certain cases the results exhibit significant differences.

Multi-laboratory evaluation of the rapid genoserotyping array (SGSA) for the identification of Salmonella serovars

Salmonella serotyping is an essential first step for identification of isolates associated with disease outbreaks. The Salmonella genoserotyping array (SGSA) is a microarray-based alternative to standard serotyping designed to rapidly identify 57 of the most commonly reported serovars through detection of the genes encoding surface O and H antigens and reporting the corresponding serovar in accordance with the existing White-Kaufmann-Le Minor serotyping scheme. In this study, we evaluated the SGSA at 4 laboratories in 3 countries by testing 1874 isolates from human and non-human sources. The SGSA correctly identified 96.7% of isolates from the target 57 serovars. For the prevalent and clinically important Salmonella serovars Enteritidis and Typhimurium, test specificity and sensitivity were greater than 98% and 99%, respectively. Due to its high-throughput nature, the SGSA is a rapid and cost-effective alternative to standard serotyping for identifying the most prevalent serovars of Salmonella.

The genome and proteome of the Kluyvera bacteriophage Kvp1 – another member of the T7-like Autographivirinae

BackgroundKluyvera, a genus within the family Enterobacteriaceae, is an infrequent cause of human infections. Bacteriophage Kvp1, the only bacteriophage isolated for one of its species, Kluyvera cryocrescens, is a member of the viral family Podoviridae.ResultsThe genome of Kvp1, the first Kluyvera cryocrescens-specific bacteriophage, was sequenced using pyrosequencing (454 technology) at the McGill University and Genome Québec Innovation Centre. The two contigs were closed using PCR and the sequence of the terminal repeats completed by primer walking off the phage DNA. The phage structural proteome was investigated by SDS-PAGE and mass spectrometry.ConclusionAt 39,472 bp, the annotated genome revealed a closer relationship to coliphage T3 than T7 with Kvp1 containing homologs to T3 early proteins S-adenosyl-L-methionine hydrolase (0.3) and protein kinase (0.7). The quantitative nature of the relationships between Kvp1 and the other members of the T7-like virus genus (T7, T3, φA1122, φYeO3-12, Berlin, K1F, VP4 and gh-1) was confirmed using CoreGenes.

In silico Characterization of DNA Motifs with Particular Reference to Promoters and Terminators

Knowledge of the regulatory elements contained within bacteriophage genomes forms the basis for understanding genomic expression and organization. The in silico prediction of promoter and terminator sequences in phage genomes is a first step towards this understanding. In this chapter, a number of programs and resources to identify regulatory elements are listed and discussed. Combining the available web-resources and literature data optimizes these predictions and can thus aid in a more directed experimental identification of these regulatory elements.