Comparative genomics and transcriptomics of lineages I, II, and III strains of Listeria monocytogenes

Genetic distance in the whole-genome perspective on Listeria monocytogenes strains F2-382 and NIHS-28 that show similar subtyping results

Background

Genome subtyping approaches could provide useful epidemiological information regarding food pathogens. However, the full genomic diversity of strains that show similar subtyping results has not yet been completely explored. Most subtyping methods are based on the differences of only a portion of the genome. We investigated two draft genome sequences of Listeria monocytogenes strain F2-382 and NIHS-28, which have been identified as closely related strains by subtyping (identical multi-virulence-locus sequence typing and multiple-locus variable number tandem repeat analysis sequence types and very similar pulsed-field gel electrophoresis patterns), despite their different sources.

Results

Two closely related strains were compared by genome structure analysis, recombination analysis, and single nucleotide polymorphism (SNP) analysis. Both genome structure analysis and recombination analysis showed that these two strains are more closely related than other strains, from a whole-genome perspective. However, the analysis of SNPs indicated that the two strains differ at the single nucleotide level.

Conclusion

We show the relationship between the results of genome subtyping and whole-genome sequencing. It appears that the relationships among strains indicated by genome subtyping methods are in accord with the relationships indicated by whole-genome analysis. However, our results also indicate that the genetic distance between the closely related strains is greater than that between clonal strains. Our results demonstrate that subtyping methods using a part of the genome are reliable in assessing the genetic distance of the strains. Furthermore, the genetic differences in the same subtype strains may provide useful information to distinguish the bacterial strains.

Electronic supplementary material

The online version of this article (doi:10.1186/s12866-014-0309-0) contains supplementary material, which is available to authorized users.

Detection of very long antisense transcripts by whole transcriptome RNA-Seq analysis of Listeria monocytogenes by semiconductor sequencing technology

The Gram-positive bacterium Listeria monocytogenes is the causative agent of listeriosis, a severe food-borne infection characterised by abortion, septicaemia, or meningoencephalitis. L. monocytogenes causes outbreaks of febrile gastroenteritis and accounts for community-acquired bacterial meningitis in humans. Listeriosis has one of the highest mortality rates (up to 30%) of all food-borne infections. This human pathogenic bacterium is an important model organism for biomedical research to investigate cell-mediated immunity. L. monocytogenes is also one of the best characterised bacterial systems for the molecular analysis of intracellular parasitism. Recently several transcriptomic studies have also made the ubiquitous distributed bacterium as a model to understand mechanisms of gene regulation from the environment to the infected host on the level of mRNA and non-coding RNAs (ncRNAs). We have used semiconductor sequencing technology for RNA-seq to investigate the repertoire of listerial ncRNAs under extra- and intracellular growth conditions. Furthermore, we applied a new bioinformatic analysis pipeline for detection, comparative genomics and structural conservation to identify ncRNAs. With this work, in total, 741 ncRNA locations of potential ncRNA candidates are now known for L. monocytogenes, of which 611 ncRNA candidates were identified by RNA-seq. 441 transcribed ncRNAs have never been described before. Among these, we identified novel long non-coding antisense RNAs with a length of up to 5,400 nt e.g. opposite to genes coding for internalins, methylases or a high-affinity potassium uptake system, namely the kdpABC operon, which were confirmed by qRT-PCR analysis. RNA-seq, comparative genomics and structural conservation of L. monocytogenes ncRNAs illustrate that this human pathogen uses a large number and repertoire of ncRNA including novel long antisense RNAs, which could be important for intracellular survival within the infected eukaryotic host.

A novel glucose 6-phosphate isomerase from Listeria monocytogenes

D-Arabinose 5-phosphate isomerases (APIs) catalyze the interconversion of D-ribulose 5-phosphate and D-arabinose 5-phosphate (A5P). A5P is an intermediate in the biosynthesis of 3-deoxy-D-manno-octulosonate (Kdo), an essential component of lipopolysaccharide, the lipopolysaccharide found in the outer membrane of Gram-negative bacteria. The genome of the Gram-positive pathogen Listeria monocytogenes contains a gene encoding a putative sugar isomerase domain API, Q723E8, with significant similarity to c3406, the only one of four APIs from Escherichia coli CFT073 that lacks a cystathionine-β-synthase domain. However, L. monocytogenes lacks genes encoding any of the other enzymes of the Kdo biosynthesis pathway. Realizing that the discovery of an API in a Gram-positive bacterium could provide insight into an alternate physiological role of A5P in the cell, we prepared and purified recombinant Q723E8. We found that Q723E8 does not possess API activity, but instead is a novel GPI (D-glucose 6-phosphate isomerase). However, the GPI activity of Q723E8 is weak compared with previously described GPIS. L. monocytogenes contains an ortholog of the well-studied two-domain bacterial GPI, so this maybe redundant. Based on this evidence glucose utilization is likely not the primary physiological role of Q723E8.

Bacterial and cellular RNAs at work during Listeria infection

ABSTRACT 

Listeria monocytogenes is an intracellular pathogen that can enter and invade host cells. In the course of its infection, RNA-mediated regulatory mechanisms provide a fast and versatile response for both the bacterium and the host. They regulate a variety of processes, such as environment sensing and virulence in pathogenic bacteria, as well as development, cellular differentiation, metabolism and immune responses in eukaryotic cells. The aim of this article is to summarize first the RNA-mediated regulatory mechanisms that play a role in the Listeria lifestyle and in its virulence, and then the host miRNA responses to Listeria infection. Finally, we discuss the potential cross-talk between bacterial RNAs and host RNA regulatory mechanisms as new mechanisms of bacterial virulence.

Genome comparison of Listeria monocytogenes serotype 4a strain HCC23 with selected lineage I and lineage II L. monocytogenes strains and other Listeria strains

More than 98% of reported human listeriosis cases are caused by specific serotypes within genetic lineages I and II. The genome sequence of Listeria monocytogenes lineage III strain HCC23 (serotype 4a) enables whole genomic comparisons across all three L. monocytogenes lineages. Protein cluster analysis indicated that strain HCC23 has the most unique protein pairs with nonpathogenic species Listeria innocua. Orthology analysis of the genome sequences of representative strains from the three L. monocytogenes genetic lineages and L. innocua (CLIP11262) identified 319 proteins unique to nonpathogenic strains HCC23 and CLIP11262 and 58 proteins unique to pathogenic strains F2365 and EGD-e. BLAST comparison of these proteins with all the sequenced L. monocytogenes and L. innocua revealed 126 proteins unique to serotype 4a and/or L. innocua; 14 proteins were only found in pathogenic serotypes. Some of the 58 proteins unique to pathogenic strains F2365 and EGD-e were previously published and are already known to contribute to listerial virulence.

Genome sequencing of Listeria monocytogenes "Quargel" listeriosis outbreak strains reveals two different strains with distinct in vitro virulence potential

A large listeriosis outbreak occurred in Austria, Germany and the Czech Republic in 2009 and 2010. The outbreak was traced back to a traditional Austrian curd cheese called “Quargel” which was contaminated with two distinct serovar 1/2a Listeria monocytogenes strains (QOC1 and QOC2). In this study we sequenced and analysed the genomes of both outbreak strains in order to investigate the extent of genetic diversity between the two strains belonging to MLST sequence types 398 (QOC2) and 403 (QOC1). Both genomes are highly similar, but also display distinct properties: The QOC1 genome is approximately 74 kbp larger than the QOC2 genome. In addition, the strains harbour 93 (QOC1) and 45 (QOC2) genes encoding strain-specific proteins. A 21 kbp region showing highest similarity to plasmid pLMIV encoding three putative internalins is integrated in the QOC1 genome. In contrast to QOC1, strain QOC2 harbours a vip homologue, which encodes a LPXTG surface protein involved in cell invasion. In accordance, in vitro virulence assays revealed distinct differences in invasion efficiency and intracellular proliferation within different cell types. The higher virulence potential of QOC1 in non-phagocytic cells may be explained by the presence of additional internalins in the pLMIV-like region, whereas the higher invasion capability of QOC2 into phagocytic cells may be due to the presence of a vip homologue. In addition, both strains show differences in stress-related gene content. Strain QOC1 encodes a so-called stress survival islet 1, whereas strain QOC2 harbours a homologue of the uncharacterized LMOf2365_0481 gene. Consistently, QOC1 shows higher resistance to acidic, alkaline and gastric stress. In conclusion, our results show that strain QOC1 and QOC2 are distinct and did not recently evolve from a common ancestor.

The Listeria monocytogenes LPXTG surface protein Lmo1413 is an invasin with capacity to bind mucin

Many Gram-positive bacterial pathogens use surface proteins covalently anchored to the peptidoglycan to cause disease. Bacteria of the genus Listeria have the largest number of surface proteins of this family. Every Listeria genome sequenced to date contains more than forty genes encoding surface proteins bearing anchoring-domains with an LPXTG motif that is recognized for covalent linkage to the peptidoglycan. About one-third of these proteins are present exclusively in pathogenic Listeria species, with some of them acting as adhesins or invasins that promote bacterial entry into eukaryotic cells. Here, we investigated two LPXTG surface proteins of the pathogen L. monocytogenes, Lmo1413 and Lmo2085, of unknown function and absent in non-pathogenic Listeria species. Lack of these two proteins does not affect bacterial adhesion or invasion of host cells using in vitro infection models. However, expression of Lmo1413 promotes entry of the non-invasive species L. innocua into non-phagocytic host cells, an effect not observed with Lmo2085. Moreover, overproduction of Lmo1413, but not Lmo2085, increases the invasion rate in non-phagocytic eukaryotic cells of an L. monocytogenes mutant deficient in the acting-binding protein ActA. Unexpectedly, production of full-length Lmo1413 and InlA exhibited opposite trends in a high percentage of L. monocytogenes isolates obtained from different sources. The idea of Lmo1413 playing a role as a new auxiliary invasin was also sustained by assays revealing that purified Lmo1413 binds to mucin via its MucBP domains. Taken together, these data indicate that Lmo1413, which we rename LmiA, for Listeria-mucin-binding invasin-A, may promote interaction of bacteria with adhesive host protective components and, in this manner, facilitate bacterial entry.

Ultra deep sequencing of Listeria monocytogenes sRNA transcriptome revealed new antisense RNAs

Listeria monocytogenes, a gram-positive pathogen, and causative agent of listeriosis, has become a widely used model organism for intracellular infections. Recent studies have identified small non-coding RNAs (sRNAs) as important factors for regulating gene expression and pathogenicity of L. monocytogenes. Increased speed and reduced costs of high throughput sequencing (HTS) techniques have made RNA sequencing (RNA-Seq) the state-of-the-art method to study bacterial transcriptomes. We created a large transcriptome dataset of L. monocytogenes containing a total of 21 million reads, using the SOLiD sequencing technology. The dataset contained cDNA sequences generated from L. monocytogenes RNA collected under intracellular and extracellular condition and additionally was size fractioned into three different size ranges from <40 nt, 40–150 nt and >150 nt. We report here, the identification of nine new sRNAs candidates of L. monocytogenes and a reevaluation of known sRNAs of L. monocytogenes EGD-e. Automatic comparison to known sRNAs revealed a high recovery rate of 55%, which was increased to 90% by manual revision of the data. Moreover, thorough classification of known sRNAs shed further light on their possible biological functions. Interestingly among the newly identified sRNA candidates are antisense RNAs (asRNAs) associated to the housekeeping genes purA, fumC and pgi and potentially their regulation, emphasizing the significance of sRNAs for metabolic adaptation in L. monocytogenes.

Comparison of the major virulence-related genes of Listeria monocytogenes in internalin A truncated strain 36-25-1 and a clinical wild-type strain

Background

Internalin A (InlA) facilitates the invasion of Listeria monocytogenes into a host cell. Some strains of Listeria monocytogenes express truncated forms of InlA, which reduces invasiveness. However, few virulence-related genes other than inlA have been analyzed in InlA-truncated strains. In the present study, we sequenced the draft genome of strain 36-25-1, an InlA-truncated strain, with pyrosequencing and compared 36 major virulence-related genes in this strain and a clinical wild-type strain.

Results

Strain 36-25-1 possessed all of the virulence-related genes analyzed. Of the analyzed genes, only 4 genes (dltA, gtcA, iap, and inlA) differed when the nucleotide sequences of strain 36-25-1 and the clinical wild-type strain were compared. Analysis of the deduced amino acid sequences found no mutations that significantly influenced virulence in genes other than inlA.

Conclusions

The virulence-associated genes in strain 36-25-1 differ little from those of the clinical wild-type strain, indicating that a slight mutation in the nucleotide sequence determines the virulence of the InlA-truncated strain. In addition, the results suggest that, aside from InlA-mediated cell invasiveness, there is almost no difference between the virulence of strain 36-25-1 and that of the clinical wild-type strain.

A PNPase dependent CRISPR System in Listeria

The human bacterial pathogen Listeria monocytogenes is emerging as a model organism to study RNA-mediated regulation in pathogenic bacteria. A class of non-coding RNAs called CRISPRs (clustered regularly interspaced short palindromic repeats) has been described to confer bacterial resistance against invading bacteriophages and conjugative plasmids. CRISPR function relies on the activity of CRISPR associated (cas) genes that encode a large family of proteins with nuclease or helicase activities and DNA and RNA binding domains. Here, we characterized a CRISPR element (RliB) that is expressed and processed in the L. monocytogenes strain EGD-e, which is completely devoid of cas genes. Structural probing revealed that RliB has an unexpected secondary structure comprising basepair interactions between the repeats and the adjacent spacers in place of canonical hairpins formed by the palindromic repeats. Moreover, in contrast to other CRISPR-Cas systems identified in Listeria, RliB-CRISPR is ubiquitously present among Listeria genomes at the same genomic locus and is never associated with the cas genes. We showed that RliB-CRISPR is a substrate for the endogenously encoded polynucleotide phosphorylase (PNPase) enzyme. The spacers of the different Listeria RliB-CRISPRs share many sequences with temperate and virulent phages. Furthermore, we show that a cas-less RliB-CRISPR lowers the acquisition frequency of a plasmid carrying the matching protospacer, provided that trans encoded cas genes of a second CRISPR-Cas system are present in the genome. Importantly, we show that PNPase is required for RliB-CRISPR mediated DNA interference. Altogether, our data reveal a yet undescribed CRISPR system whose both processing and activity depend on PNPase, highlighting a new and unexpected function for PNPase in “CRISPRology”.

CRISPR-Cas systems confer to bacteria and archaea an adaptive immunity that protects them against invading bacteriophages and plasmids. In this study, we characterize a CRISPR (RliB-CRISPR) that is present in all L. monocytogenes strains at the same genomic locus but is never associated with a cas operon. It is an unusual CRISPR that, as we demonstrate, has a secondary structure consisting of basepair interactions between the repeat sequence and the adjacent spacer. We show that the RliB-CRISPR is processed by the endogenously encoded polynucleotide phosphorylase enzyme (PNPase). In addition, we show that the RliB-CRISPR system requires PNPase and presence of trans encoded cas genes of a second CRISPR-Cas system, to mediate DNA interference directed against a plasmid carrying a matching protospacer. Altogether, our data reveal a novel type of CRISPR system in bacteria that requires endogenously encoded PNPase enzyme for its processing and interference activity.

Virulence characterization of Listeria monocytogenes

Listeria monocytogenes, a foodborne intracellular pathogen, is responsible for about 300 deaths every year in the USA. It has the ability to escape host defense mechanisms and causes listeriosis more frequently in immunocompromised individuals. Virulence mechanisms in L. monocytogenes are highly regulated and tightly controlled. A number of virulence factors that play important roles in pathogenesis of listeriosis have been identified and characterized. This review highlights the power of comparative genomics and functional genomics in identifying genes and proteins involved in the infection process. These genes and proteins are potentially useful as biomarkers for detecting virulent L. monocytogenes. This review also focuses on developments in the in vivo and in vitro models used in characterization of listerial virulence.

Molecular approaches to the identification of pathogenic and nonpathogenic listeriae

The genus Listeria consists of a closely related group of Gram-positive bacteria that commonly occur in the environment and demonstrate varied pathogenic potential. Of the 10 species identified to date, L. monocytogenes is a facultative intracellular pathogen of both humans and animals, L. ivanovii mainly infects ungulates (eg., sheep and cattle), while other species (L. innocua, L. seeligeri, L. welshimeri, L. grayi, L. marthii, L. rocourtiae, L. fleischmannii and L. weihenstephanensis) are essentially saprophytes. Within the species of L. monocytogenes, several serovars (e.g., 4b, 1/2a, 1/2b and 1/2c) are highly pathogenic and account for a majority of clinical isolations. Due to their close morphological, biological, biochemical and genetic similarities, laboratory identification of pathogenic and nonpathogenic Listeria organisms is technically challenging. With the development and application of various molecular approaches, accurate and rapid discrimination of pathogenic and nonpathogenic Listeria organisms, as well as pathogenic and nonpathogenic L. monocytogenes strains, has become possible.

Listeria phages: Genomes, evolution, and application

Listeria is an important foodborne pathogen and the causative agent of Listeriosis, a potentially fatal infection. Several hundred Listeria bacteriophages have been described over the past decades, but only few have actually been characterized in some detail, and genome sequences are available for less than twenty of them. We here present an overview of what is currently known about Listeria phage genomics, their role in host evolution and pathogenicity, and their various applications in biotechnology and diagnostics.

Fluorescence amplified fragment length polymorphism compared to pulsed field gel electrophoresis for Listeria monocytogenes subtyping

Background

Listeriosis is a severe infection which mainly affects pregnant women, neonates and immuno-compromised adults. ANSES’s Laboratory for Food safety has been the European Union Reference Laboratory (EURL) for L. monocytogenes in the food chain since 2006. Pulsed Field Gel Electrophoresis (PFGE) is routinely used in the EURL for the surveillance of L. monocytogenes isolated from foods, animals and the environment. One of the main EURL activities is to evaluate alternative molecular subtyping methods to PFGE, and integrate their use within the National Reference Laboratories (NRL) network. Since 2008, the United Kingdom (UK)-NRL for L. monocytogenes at the Health Protection Agency (HPA), London, has used fluorescent Amplified Fragment Length Polymorphism (fAFLP) for the routine surveillance of L. monocytogenes isolated from human clinical cases, food and food processing environments in the UK. This study compares fAFLP with PFGE for subtyping L. monocytogenes.

Results

A panel of 109 L. monocytogenes isolates from either human cases of listeriosis, foods, food processing environments and animals were used for the comparative evaluation. Among these, 2 strains were tested from duplicate culture by both methods. The panel also included field isolates, isolates associated with outbreaks or sporadic cases and reference strains. The two strains tested in duplicate displayed the same fAFLP and PFGE types. Strains known to be epidemiologically associated with one another were found to have unique PFGE and fAFLP types. FAFLP and PFGE divided the strains into 76 and 82 distinct profiles, or types, respectively. The discriminatory index calculated was 0.993 and 0.996 for fAFLP and PFGE, respectively.

Conclusions

The discriminatory ability of fAFLP was similar to that of PFGE for the subtyping of L. monocytogenes isolates. As a less labour intensive technique fAFLP may be a better method to use than PFGE in investigating outbreaks of human listeriosis and tracking the source of contamination in food processing facilities in real time.

Reassessment of the Listeria monocytogenes pan-genome reveals dynamic integration hotspots and mobile genetic elements as major components of the accessory genome

Background

Listeria monocytogenes is an important food-borne pathogen and model organism for host-pathogen interaction, thus representing an invaluable target considering research on the forces governing the evolution of such microbes. The diversity of this species has not been exhaustively explored yet, as previous efforts have focused on analyses of serotypes primarily implicated in human listeriosis. We conducted complete genome sequencing of 11 strains employing 454 GS FLX technology, thereby achieving full coverage of all serotypes including the first complete strains of serotypes 1/2b, 3c, 3b, 4c, 4d, and 4e. These were comparatively analyzed in conjunction with publicly available data and assessed for pathogenicity in the Galleria mellonella insect model.

Results

The species pan-genome of L. monocytogenes is highly stable but open, suggesting an ability to adapt to new niches by generating or including new genetic information. The majority of gene-scale differences represented by the accessory genome resulted from nine hyper variable hotspots, a similar number of different prophages, three transposons (Tn916, Tn554, IS3-like), and two mobilizable islands. Only a subset of strains showed CRISPR/Cas bacteriophage resistance systems of different subtypes, suggesting a supplementary function in maintenance of chromosomal stability. Multiple phylogenetic branches of the genus Listeria imply long common histories of strains of each lineage as revealed by a SNP-based core genome tree highlighting the impact of small mutations for the evolution of species L. monocytogenes. Frequent loss or truncation of genes described to be vital for virulence or pathogenicity was confirmed as a recurring pattern, especially for strains belonging to lineages III and II. New candidate genes implicated in virulence function were predicted based on functional domains and phylogenetic distribution. A comparative analysis of small regulatory RNA candidates supports observations of a differential distribution of trans-encoded RNA, hinting at a diverse range of adaptations and regulatory impact.

Conclusions

This study determined commonly occurring hyper variable hotspots and mobile elements as primary effectors of quantitative gene-scale evolution of species L. monocytogenes, while gene decay and SNPs seem to represent major factors influencing long-term evolution. The discovery of common and disparately distributed genes considering lineages, serogroups, serotypes and strains of species L. monocytogenes will assist in diagnostic, phylogenetic and functional research, supported by the comparative genomic GECO-LisDB analysis server (http://bioinfo.mikrobio.med.uni-giessen.de/geco2lisdb).