Paolo Vatta

Comparative genomics of the bacterial genus Listeria: Genome evolution is characterized by limited gene acquisition and limited gene loss

BackgroundThe bacterial genus Listeria contains pathogenic and non-pathogenic species, including the pathogens L. monocytogenes and L. ivanovii, both of which carry homologous virulence gene clusters such as the prfA cluster and clusters of internalin genes. Initial evidence for multiple deletions of the prfA cluster during the evolution of Listeria indicates that this genus provides an interesting model for studying the evolution of virulence and also presents practical challenges with regard to definition of pathogenic strains.ResultsTo better understand genome evolution and evolution of virulence characteristics in Listeria, we used a next generation sequencing approach to generate draft genomes for seven strains representing Listeria species or clades for which genome sequences were not available. Comparative analyses of these draft genomes and six publicly available genomes, which together represent the main Listeria species, showed evidence for (i) a pangenome with 2,032 core and 2,918 accessory genes identified to date, (ii) a critical role of gene loss events in transition of Listeria species from facultative pathogen to saprotroph, even though a consistent pattern of gene loss seemed to be absent, and a number of isolates representing non-pathogenic species still carried some virulence associated genes, and (iii) divergence of modern pathogenic and non-pathogenic Listeria species and strains, most likely circa 47 million years ago, from a pathogenic common ancestor that contained key virulence genes.ConclusionsGenome evolution in Listeria involved limited gene loss and acquisition as supported by (i) a relatively high coverage of the predicted pan-genome by the observed pan-genome, (ii) conserved genome size (between 2.8 and 3.2 Mb), and (iii) a highly syntenic genome. Limited gene loss in Listeria did include loss of virulence associated genes, likely associated with multiple transitions to a saprotrophic lifestyle. The genus Listeria thus provides an example of a group of bacteria that appears to evolve through a loss of virulence rather than acquisition of virulence characteristics. While Listeria includes a number of species-like clades, many of these putative species include clades or strains with atypical virulence associated characteristics. This information will allow for the development of genetic and genomic criteria for pathogenic strains, including development of assays that specifically detect pathogenic Listeria strains.

Escherichia coli serotype O55:H7 diversity supports parallel acquisition of bacteriophage at Shiga toxin phage insertion sites during evolution of the O157:H7 lineage

Enteropathogenic Escherichia coli (EPEC) continues to be a leading cause of mortality and morbidity in children around the world. Two EPEC genomes have been fully sequenced: those of EPEC O127:H6 strain E2348/69 (United Kingdom, 1969) and EPEC O55:H7 strain CB9615 (Germany, 2003). The O55:H7 serotype is a recent precursor to the virulent enterohemorrhagic E. coli O157:H7. To explore the diversity of O55:H7 and better understand the clonal evolution of O157:H7, we fully sequenced EPEC O55:H7 strain RM12579 (California, 1974), which was collected 1 year before the first U.S. isolate of O157:H7 was identified in California. Phage-related sequences accounted for nearly all differences between the two O55:H7 strains. Additionally, O55:H7 and O157:H7 strains were tested for the presence and insertion sites of Shiga toxin gene (stx)-containing bacteriophages. Analysis of non-phage-associated genes supported core elements of previous O157:H7 stepwise evolutionary models, whereas phage composition and insertion analyses suggested a key refinement. Specifically, the placement and presence of lambda-like bacteriophages (including those containing stx) should not be considered stable evolutionary markers or be required in placing O55:H7 and O157:H7 strains within the stepwise evolutionary models. Additionally, we suggest that a 10.9-kb region (block 172) previously believed unique to O55:H7 strains can be used to identify early O157:H7 strains. Finally, we defined two subsets of O55:H7 strains that share an as-yet-unobserved or extinct common ancestor with O157:H7 strains. Exploration of O55:H7 diversity improved our understanding of the evolution of E. coli O157:H7 and suggested a key revision to accommodate existing and future configurations of stx-containing bacteriophages into current models.

RNA-Seq Mapping and Detection of Gene Fusions with a Suffix Array Algorithm

High-throughput RNA sequencing enables quantification of transcripts (both known and novel), exon/exon junctions and fusions of exons from different genes. Discovery of gene fusions–particularly those expressed with low abundance– is a challenge with short- and medium-length sequencing reads. To address this challenge, we implemented an RNA-Seq mapping pipeline within the LifeScope software. We introduced new features including filter and junction mapping, annotation-aided pairing rescue and accurate mapping quality values. We combined this pipeline with a Suffix Array Spliced Read (SASR) aligner to detect chimeric transcripts. Performing paired-end RNA-Seq of the breast cancer cell line MCF-7 using the SOLiD system, we called 40 gene fusions among over 120,000 splicing junctions. We validated 36 of these 40 fusions with TaqMan assays, of which 25 were expressed in MCF-7 but not the Human Brain Reference. An intra-chromosomal gene fusion involving the estrogen receptor alpha gene ESR1, and another involving the RPS6KB1 (Ribosomal protein S6 kinase beta-1) were recurrently expressed in a number of breast tumor cell lines and a clinical tumor sample.