Analysis of Microbial Composition of Edible Insect Products Available for Human Consumption within the United States Using Traditional Microbiological Methods and Whole Genome Sequencing

Edible insects offer a promising protein source for humans, but their food safety risks have not been previously investigated within the United States. Therefore, the aim of this study was to investigate the microbial content of processed edible insect products. A total of eight different types of edible insect products, including diving beetles, silkworms, grasshoppers, Jamaican crickets, mealworms, mole crickets, whole roasted crickets, and 100% pure cricket powder, were purchased from a large online retailer for the analysis. All the products were purchased in August 2022 and examined between August 2022 and November 2022. Traditional microbiological methods were employed to determine microbial counts for each product type using three replicates (total number of samples = 24). This included assessing aerobic bacterial spore, lactic acid bacteria, Enterobacteriaceae, total viable counts, and the presence of Salmonella. Additionally, whole genome sequencing was employed to further characterize selected colonies (n = 96). Microbial counts data were statistically analyzed using one-way ANOVA, while sequence data were taxonomically classified using Sepia.Bacilluscereusgroup isolates underwent additional characterization with Btyper3. Product type significantly influenced total viable counts, bacterial spore counts, and lactic acid bacteria counts (P = 0.00391, P = 0.0065, and P < 0.001, respectively), with counts ranging from < 1.70 to 6.01 Log10 CFU/g, <1.70 to 5.25 Log10 CFU/g, and < 1.70 to 4.86 Log10 CFU/g, respectively. Enterobacteriaceae were only detected in mole crickets (<2.30 Log10 CFU/g) and house cricket powder (<2.15 Log10 CFU/g). All samples were negative for Salmonella. Whole genome sequencing revealed the presence of 12 different bacterial genera among the analyzed isolates, with a majority belonging to the Bacillus genus. Some of the isolates of Bacillus cereus group were identified as biovar Emeticus. Overall, although edible insects offer a promising food alternative, the presence of Bacillus cereus group in some products could raise concerns regarding food safety.

Characterizing the gut microbiome of broilers raised under conventional and no antibiotics ever practices

Meat from broilers raised without the use of antibiotics is becoming increasingly popular among consumers. Consequently, interest in the microbial profiling of chickens produced under nonconventional practices is growing, however, research on this topic is lacking. The current study was designed to characterize the dynamics of gut microbial populations of broilers raised under conventional and no antibiotics ever (NAE) practices. Four commercial farms (2 conventional and 2 NAE) were included in this study. On each farm, cecal (n = 224) and ileal (n = 224) contents were collected from birds at different stages during the grow out of a single flock and following transportation to the processing facility. Cecal microbiota was dominated by the genera Escherichia and Enterococcus upon hatching in both conventional and NAE flocks, shifting with time toward predominantly Faecalibacterium and Bacteroides. The composition of cecal microbial communities of NAE broilers was different than that of conventional chickens (P ≤ 0.05). Conventional broilers harbored a rich, but less diverse cecal microbiota than NAE, while the ileal microbiota was primarily populated with genera previously named Lactobacillus, which exhibited a higher abundance in NAE broilers (P ≤ 0.05). In both production systems, the microbiota followed a similar temporal succession that was more evident in the ceca. Transportation to the processing plant impacted the microbial composition of the ileum (P ≤ 0.05), characterized by an increase in the relative abundance of Psychrobacter. Finally, differential abundance analysis showed a positive correlation between Campylobacter and Enorma within the cecum microbiota, and a negative correlation with Salmonella.

Phylogeny of the Bacillus altitudinis Complex and Characterization of a Newly Isolated Strain with Antilisterial Activity

ABSTRACT

Bacillus strain UTK D1-0055 was isolated from a laboratory environment and appeared to have antilisterial activity. The genome was sequenced, the strain was identified as Bacillus altitudinis, and a high-quality complete annotated genome was produced. The taxonomy was evaluated for this and related Bacillus species (B. aerophilus, B. pumilus, B. safensis, B. stratosphericus, and B. xiamenensis) because the taxonomy is unclear and contains errors in public databases such as NCBI. The included strains grouped into seven clusters based on average nucleotide identity. Strains designated as B. aerophilus, B. altitudinis, and B. stratosphericus grouped together in the cluster containing the B. altitudinis type strain, suggesting that these three species should be considered a single species, B. altitudinis. The antimicrobial activity of UTK D1-0055 was evaluated against a panel of 15 Listeria strains (nine Listeria monocytogenes serotypes, Listeria innocua, and Listeria marthii), other foodborne pathogens (six Salmonella enterica serotypes and Escherichia coli), and three representative fungi (Saccharomyces cerevisiae, Botrytis cinerea, and Hyperdermium pulvinatum). Antibacterial activity was observed against all Listeria strains, but no antibacterial effects were found against the other tested bacterial and fungal strains. Biosynthetic gene clusters were identified in silico that may be related to the observed antibacterial activity, and these clusters included genes that putatively encode bacteriocins and nonribosomally synthesized peptides. The B. altitudinis strain identified in this investigation had a broad range of antilisterial activity, suggesting that it and other related strains may be useful for biocontrol in the food industry.

HIGHLIGHTS

SneakerNet: A modular quality assurance and quality check workflow for primary genomic and metagenomic read data

Laboratories that run Whole Genome Sequencing (WGS) produce a tremendous amount of data, up to 10 gigabytes for some common instruments. There is a need to standardize the quality assurance and quality control process (QA/QC). Therefore we have created SneakerNet to automate the QA/QC for WGS.

Quantitative dynamics of Salmonella and E. coli in feces of feedlot cattle treated with ceftiofur and chlortetracycline

Antibiotic use in beef cattle is a risk factor for the expansion of antimicrobial-resistant Salmonella populations. However, actual changes in the quantity of Salmonella in cattle feces following antibiotic use have not been investigated. Previously, we observed an overall reduction in Salmonella prevalence in cattle feces associated with both ceftiofur crystalline-free acid (CCFA) and chlortetracycline (CTC) use; however, during the same time frame the prevalence of multidrug-resistant Salmonella increased. The purpose of this analysis was to quantify the dynamics of Salmonella using colony counting (via a spiral-plating method) and hydrolysis probe-based qPCR (TaqMan® qPCR). Additionally, we quantified antibiotic-resistant Salmonella by plating to agar containing antibiotics at Clinical & Laboratory Standards Institute breakpoint concentrations. Cattle were randomly assigned to 4 treatment groups across 16 pens in 2 replicates consisting of 88 cattle each. Fecal samples from Days 0, 4, 8, 14, 20, and 26 were subjected to quantification assays. Duplicate qPCR assays targeting the Salmonella invA gene were performed on total community DNA for 1,040 samples. Diluted fecal samples were spiral plated on plain Brilliant Green Agar (BGA) and BGA with ceftriaxone (4 μg/ml) or tetracycline (16 μg/ml). For comparison purposes, indicator non-type-specific (NTS) E. coli were also quantified by direct spiral plating. Quantity of NTS E. coli and Salmonella significantly decreased immediately following CCFA treatment. CTC treatment further decreased the quantity of Salmonella but not NTS E. coli. Effects of antibiotics on the imputed log10 quantity of Salmonella were analyzed via a multi-level mixed linear regression model. The invA gene copies decreased with CCFA treatment by approximately 2 log10 gene copies/g feces and remained low following additional CTC treatment. The quantities of tetracycline or ceftriaxone-resistant Salmonella were approximately 4 log10 CFU/g feces; however, most of the samples were under the quantification limit. The results of this study demonstrate that antibiotic use decreases the overall quantity of Salmonella in cattle feces in the short term; however, the overall quantities of antimicrobial-resistant NTS E. coli and Salmonella tend to remain at a constant level throughout.

Characterization of Antimicrobial Resistant (AMR) Salmonella Enterica Isolates Associated With Cattle at Harvest in Mexico

Despite being the target of control efforts for many decades, Salmonella enterica continues to be linked with a large amount of foodborne illnesses and outbreaks worldwide. Over the years, Salmonella isolated from meat products have exhibited a high level of antibiotic resistance. In this study, a total of 351 Salmonella isolates, recovered from cattle fecal samples (n = 31), hides (n = 105), and beef carcasses (n = 215) from 3 abattoirs in Mexico were analyzed for antimicrobial susceptibility. Resistance to at least one antimicrobial drug was found in 205 (58.4%) isolates and 20 different resistance phenotypes were observed among this Salmonella isolates set. Resistance to tetracycline (40.2%) and nalidixic acid (21.1%) was most commonly observed. Additionally, the most common multidrug-resistant (MDR) phenotypes shared resistance to chloramphenicol, streptomycin, tetracycline, and trimethopin/sulfamethoxazole (11.3%), resistance to ampicillin, tetracycline, and trimethopin/sulfamethoxazole (3.4%), and resistance to ampicillin, streptomycin, and tetracycline (2.5%). When it came to antimicrobial resistance phenotypes in each abattoir, we determined there was no statistical difference in the frequency of resistant vs. susceptible Salmonella isolates among the three abattoirs (P > 0.05). These data indicate that Salmonella isolates recovered from beef cattle in Mexico are commonly resistant to antimicrobials and often multiple antimicrobials. In Mexico, antimicrobial resistance, and in particular, multidrug-resistance, maybe of particular concern due to the much higher prevalence of Salmonella in retail beef. This may lead to the spread of resistance and to the reduction of antibiotic efficacy for the control of animal and human infections. Promoting control measures and inspection standards on imported animals and food products should be applied to avoid the spread of antibiotic resistance in various populations and among countries.

Global Microbial Identifier

Human and animal populations are increasingly confronted with emerging and re-emerging infections and often such infections are exchanged between these populations, e.g. through food. A more effective and uniform approach to the prevention of these microbial threats is essential. The technological advances in the next generation sequencing field and decreasing costs of these tests provide novel opportunities in understanding the dynamics of infection—even in real time—through the analysis of microbial genome diversity. The projected significant increase in whole (microbial) genome sequencing (WGS) will likely also enable a much better understanding of the pathogenesis of the infection and the molecular basis of the host response to infection. But the full potential of these advances will only transpire if the data in this area become transferable and thereby comparable, preferably in open-source systems. There is therefore an obvious need to develop a global system of whole microbial genome databases to aggregate, share, mine and use microbiological genomic data, to address global public health and clinical challenges, and most importantly to identify and diagnose infectious diseases. The global microbial identifier (GMI) initiative, aims to build a database of whole microbial genome sequencing data linked to relevant metadata, which can be used to identify microorganisms, their communities and the diseases they cause. It would be a platform for storing whole genome sequencing (WGS) data of microorganisms, for the identification of relevant genes and for the comparison of genomes to detect outbreaks and emerging pathogens. To harness the full potential of WGS, a shared global database of genomes linked to relevant metadata and the necessary software tools needs to be generated, hence the global microbial identifier (GMI) initiative. This tool will ideally be used in amongst others in the diagnosis of infectious diseases in humans and animals, in the identification of microorganisms in food and environment, and to track and trace microbial agents in all arenas globally. This will require standardization and extensive investments in computational analytical tools. In addition, the wider introduction of WGS in clinical diagnostics can accelerate developments in health care in many poor countries. This overview describes the growing network of stakeholders behind GMI, the contours of the database, and the IT structures needed to serve the GMI user community. It discusses what essentially can be done by a global GMI tool and how the GMI organization could help achieve these goals.

Rapid whole-genome sequencing for surveillance of Salmonella enterica serovar enteritidis

For Salmonella enterica serovar Enteritidis, 85% of isolates can be classified into 5 pulsed-field gel electrophoresis (PFGE) types. However, PFGE has limited discriminatory power for outbreak detection. Although whole-genome sequencing has been found to improve discrimination of outbreak clusters, whether this procedure can be used in real-time in a public health laboratory is not known. Therefore, we conducted a retrospective and prospective analysis. The retrospective study investigated isolates from 1 confirmed outbreak. Additional cases could be attributed to the outbreak strain on the basis of whole-genome data. The prospective study included 58 isolates obtained in 2012, including isolates from 1 epidemiologically defined outbreak. Whole-genome sequencing identified additional isolates that could be attributed to the outbreak, but which differed from the outbreak-associated PFGE type. Additional putative outbreak clusters were detected in the retrospective and prospective analyses. This study demonstrates the practicality of implementing this approach for outbreak surveillance in a state public health laboratory.

Carbapenem-Resistant Bacteria Recovered from Faeces of Dairy Cattle in the High Plains Region of the USA

Objective

A study was conducted to recover carbapenem-resistant bacteria from the faeces of dairy cattle and identify the underlying genetic mechanisms associated with reduced phenotypic susceptibility to carbapenems.

Methods

One hundred and fifty-nine faecal samples from dairy cattle were screened for carbapenem-resistant bacteria. Phenotypic screening was conducted on two media containing ertapenem. The isolates from the screening step were characterised via disk diffusion, Modified Hodge, and Carba NP assays. Carbapenem-resistant bacteria and carbapenemase-producing isolates were subjected to Gram staining and biochemical testing to include Gram-negative bacilli. Whole genome sequencing was performed on bacteria that exhibited either a carbapenemase-producing phenotype or were not susceptible to ertapenem and were presumptively Enterobacteriaceae.

Results

Of 323 isolates collected from the screening media, 28 were selected for WGS; 21 of which were based on a carbapenemase-producing phenotype and 7 were presumptively Enterobacteriaceae and not susceptible to ertapenem. Based on analysis of WGS data, isolates included: 3 Escherichia coli harbouring bla_CMY-2 and truncated ompF genes; 8 Aeromonas harbouring bla_cphA-like genes; 1 Acinetobacter baumannii harbouring a novel bla_OXA gene (bla_OXA-497); and 6 Pseudomonas with conserved domains of various carbapenemase-producing genes.

Conclusions

Carbapenem resistant bacteria appear to be rare in cattle. Nonetheless, carbapenem-resistant bacteria were detected across various genera and were found to harbour a variety of mechanisms conferring reduced susceptibility. The development and dissemination of carbapenem-resistant bacteria in livestock would have grave implications for therapeutic treatment options in human medicine; thus, continued monitoring of carbapenem susceptibility among enteric bacteria of livestock is warranted.

Targeted Amplicon Sequencing for Single-Nucleotide-Polymorphism Genotyping of Attaching and Effacing Escherichia coli O26:H11 Cattle Strains via a High-Throughput Library Preparation Technique

Enterohemorrhagic Escherichia coli (EHEC) O26:H11, a serotype within Shiga toxin-producing E. coli (STEC) that causes severe human disease, has been considered to have evolved from attaching and effacing E. coli (AEEC) O26:H11 through the acquisition of a Shiga toxin-encoding gene. Targeted amplicon sequencing using next-generation sequencing technology of 48 phylogenetically informative single-nucleotide polymorphisms (SNPs) and three SNPs differentiating Shiga toxin-positive (stx-positive) strains from Shiga toxin-negative (stx-negative) strains were used to infer the phylogenetic relationships of 178 E. coli O26:H11 strains (6 stx-positive strains and 172 stx-negative AEEC strains) from cattle feces to 7 publically available genomes of human clinical strains. The AEEC cattle strains displayed synonymous SNP genotypes with stx2-positive sequence type 29 (ST29) human O26:H11 strains, while stx1 ST21 human and cattle strains clustered separately, demonstrating the close phylogenetic relatedness of these Shiga toxin-negative AEEC cattle strains and human clinical strains. With the exception of seven stx-negative strains, five of which contained espK, three stx-related SNPs differentiated the STEC strains from non-STEC strains, supporting the hypothesis that these AEEC cattle strains could serve as a potential reservoir for new or existing pathogenic human strains. Our results support the idea that targeted amplicon sequencing for SNP genotyping expedites strain identification and genetic characterization of E. coli O26:H11, which is important for food safety and public health.

Genomic Epidemiology: Whole-Genome-Sequencing-Powered Surveillance and Outbreak Investigation of Foodborne Bacterial Pathogens

As we are approaching the twentieth anniversary of PulseNet, a network of public health and regulatory laboratories that has changed the landscape of foodborne illness surveillance through molecular subtyping, public health microbiology is undergoing another transformation brought about by so-called next-generation sequencing (NGS) technologies that have made whole-genome sequencing (WGS) of foodborne bacterial pathogens a realistic and superior alternative to traditional subtyping methods. Routine, real-time, and widespread application of WGS in food safety and public health is on the horizon. Technological, operational, and policy challenges are still present and being addressed by an international and multidisciplinary community of researchers, public health practitioners, and other stakeholders.

Transcriptomic Analysis of the Adaptation of Listeria monocytogenes to Growth on Vacuum-Packed Cold Smoked Salmon

The foodborne pathogen Listeria monocytogenes is able to survive and grow in ready-to-eat foods, in which it is likely to experience a number of environmental stresses due to refrigerated storage and the physicochemical properties of the food. Little is known about the specific molecular mechanisms underlying survival and growth of L. monocytogenes under different complex conditions on/in specific food matrices. Transcriptome sequencing (RNA-seq) was used to understand the transcriptional landscape of L. monocytogenes strain H7858 grown on cold smoked salmon (CSS; water phase salt, 4.65%; pH 6.1) relative to that in modified brain heart infusion broth (MBHIB; water phase salt, 4.65%; pH 6.1) at 7°C. Significant differential transcription of 149 genes was observed (false-discovery rate [FDR], <0.05; fold change, ≥2.5), and 88 and 61 genes were up- and downregulated, respectively, in H7858 grown on CSS relative to the genes in H7858 grown in MBHIB. In spite of these differences in transcriptomes under these two conditions, growth parameters for L. monocytogenes were not significantly different between CSS and MBHIB, indicating that the transcriptomic differences reflect how L. monocytogenes is able to facilitate growth under these different conditions. Differential expression analysis and Gene Ontology enrichment analysis indicated that genes encoding proteins involved in cobalamin biosynthesis as well as ethanolamine and 1,2-propanediol utilization have significantly higher transcript levels in H7858 grown on CSS than in that grown in MBHIB. Our data identify specific transcriptional profiles of L. monocytogenes growing on vacuum-packaged CSS, which may provide targets for the development of novel and improved strategies to control L. monocytogenes growth on this ready-to-eat food.

Characterization of the cytolethal distending toxin (typhoid toxin) in non-typhoidal Salmonella serovars

BACKGROUND

For many putative Salmonella enterica subsp. enterica virulence genes, functional characterization across serovars has been limited. Cytolethal distending toxin B (CdtB) is an incompletely characterized virulence factor that is found not only in Salmonella enterica subsp. enterica serovar Typhi (Salmonella Typhi) and dozens of Gram negative bacterial pathogens, but also in non-typhoidal Salmonella (NTS) serovars.

METHODS

A comparative genomics approach was performed to characterize sequence conservation of the typhoid toxin (TT), encoded in the CdtB-islet, between Salmonella Typhi and NTS serovars. The cytotoxic activity of representative Salmonella enterica subsp. enterica serovars Javiana, Montevideo and Schwarzengrund strains and their respective isogenic cdtB mutants was determined in human intestinal epithelial Henle-407 cells by assessment of cell cycle progression of infected cells using fluorescence-activated cell sorting (FACS). Two-way analysis of variance (ANOVA) was used to determine whether cdtB deletion had a significant (p < 0.05) effect on the percentage of Henle-407 cells at each stage of the cell cycle.

RESULTS

Here we show that a CdtB-islet encoding the cytolethal distending toxin B (CdtB), pertussis-like toxin A (PltA), and pertussis-like toxin B (PltB) is present in a dozen NTS serovars and that these proteins have a high level of sequence conservation and each form monophyletic clades with corresponding Salmonella Typhi genes. Human epithelial Henle-407 cells infected with three representative CdtB-encoding NTS serovars displayed G2/M phase cell cycle arrest that was absent in cells infected with corresponding isogenic cdtB null mutants (p < 0.0001 for the factor ∆cdtB deletion).

CONCLUSION

Our results show that CdtB encoded by NTS serovars has a genomic organization, amino acid sequence conservation and biological activity similar to the TT, and thus, may contribute to disease pathogenesis.

Genomic epidemiology of Salmonella enterica serotype Enteritidis based on population structure of prevalent lineages

Major lineages emerged during the 17th–18th centuries and diversified during the 1920s and 1950s.

Salmonella enterica serotype Enteritidis is one of the most commonly reported causes of human salmonellosis. Its low genetic diversity, measured by fingerprinting methods, has made subtyping a challenge. We used whole-genome sequencing to characterize 125 S. enterica Enteritidis and 3 S. enterica serotype Nitra strains. Single-nucleotide polymorphisms were filtered to identify 4,887 reliable loci that distinguished all isolates from each other. Our whole-genome single-nucleotide polymorphism typing approach was robust for S. enterica Enteritidis subtyping with combined data for different strains from 2 different sequencing platforms. Five major genetic lineages were recognized, which revealed possible patterns of geographic and epidemiologic distribution. Analyses on the population dynamics and evolutionary history estimated that major lineages emerged during the 17th–18th centuries and diversified during the 1920s and 1950s.

Listeria booriae sp. nov. and Listeria newyorkensis sp. nov., from food processing environments in the USA

Sampling of seafood and dairy processing facilities in the north-eastern USA produced 18 isolates of Listeria spp. that could not be identified at the species-level using traditional phenotypic and genotypic identification methods. Results of phenotypic and genotypic analyses suggested that the isolates represent two novel species with an average nucleotide blast identity of less than 92% with previously described species of the genus Listeria. Phylogenetic analyses based on whole genome sequences, 16S rRNA gene and sigB gene sequences confirmed that the isolates represented by type strain FSL M6-0635(T) and FSL A5-0209 cluster phylogenetically with Listeria cornellensis. Phylogenetic analyses also showed that the isolates represented by type strain FSL A5-0281(T) cluster phylogenetically with Listeria riparia. The name Listeria booriae sp. nov. is proposed for the species represented by type strain FSL A5-0281(T) ( =DSM 28860(T) =LMG 28311(T)), and the name Listeria newyorkensis sp. nov. is proposed for the species represented by type strain FSL M6-0635(T) ( =DSM 28861(T) =LMG 28310(T)). Phenotypic and genotypic analyses suggest that neither species is pathogenic.

Genomic characterization of Salmonella Cerro ST367, an emerging Salmonella subtype in cattle in the United States

Background

Within the last decade, Salmonella enterica subsp. enterica serovar Cerro (S. Cerro) has become one of the most common serovars isolated from cattle and dairy farm environments in the northeastern US. The fact that this serovar is commonly isolated from subclinically infected cattle and is rarely associated with human disease, despite its frequent isolation from cattle, has led to the hypothesis that this emerging serovar may be characterized by reduced virulence. We applied comparative and population genomic approaches to (i) characterize the evolution of this recently emerged serovar and to (ii) gain a better understanding of genomic features that could explain some of the unique epidemiological features associated with this serovar.

Results

In addition to generating a de novo draft genome for one Salmonella Cerro strain, we also generated whole genome sequence data for 26 additional S. Cerro isolates, including 16 from cattle operations in New York (NY) state, 2 from human clinical cases from NY in 2008, and 8 from diverse animal sources (7 from Washington state and 1 from Florida). All isolates sequenced in this study represent sequence type ST367. Population genomic analysis showed that isolates from the NY cattle operations form a well-supported clade within S. Cerro ST367 (designated here “NY bovine clade”), distinct from isolates from Washington state, Florida and the human clinical cases. A molecular clock analysis indicates that the most recent common ancestor of the NY bovine clade dates back to 1998, supporting the recent emergence of this clone.

Comparative genomic analyses revealed several relevant genomic features of S. Cerro ST367, that may be responsible for reduced virulence of S. Cerro, including an insertion creating a premature stop codon in sopA. In addition, patterns of gene deletion in S. Cerro ST367 further support adaptation of this clone to a unique ecological or host related niche.

Conclusions

Our results indicate that the increase in prevalence of S. Cerro ST367 is caused by a highly clonal subpopulation and that S. Cerro ST367 is characterized by unique genomic deletions that may indicate adaptation to specific ecological niches and possibly reduced virulence in some hosts.

Electronic supplementary material

The online version of this article (doi: 10.1186/1471-2164-15-427) contains supplementary material, which is available to authorized users.

Listeria floridensis sp. nov., Listeria aquatica sp. nov., Listeria cornellensis sp. nov., Listeria riparia sp. nov. and Listeria grandensis sp. nov., from agricultural and natural environments

Sampling of agricultural and natural environments in two US states (Colorado and Florida) yielded 18 Listeria-like isolates that could not be assigned to previously described species using traditional methods. Using whole-genome sequencing and traditional phenotypic methods, we identified five novel species, each with a genome-wide average BLAST nucleotide identity (ANIb) of less than 85% to currently described species. Phylogenetic analysis based on 16S rRNA gene sequences and amino acid sequences of 31 conserved loci showed the existence of four well-supported clades within the genus Listeria; (i) a clade representing Listeria monocytogenes, L. marthii, L. innocua, L. welshimeri, L. seeligeri and L. ivanovii, which we refer to as Listeria sensu stricto, (ii) a clade consisting of Listeria fleischmannii and two newly described species, Listeria aquatica sp. nov. (type strain FSL S10-1188(T) = DSM 26686(T) = LMG 28120(T) = BEI NR-42633(T)) and Listeria floridensis sp. nov. (type strain FSL S10-1187(T) = DSM 26687(T) = LMG 28121(T) = BEI NR-42632(T)), (iii) a clade consisting of Listeria rocourtiae, L. weihenstephanensis and three novel species, Listeria cornellensis sp. nov. (type strain TTU A1-0210(T) = FSL F6-0969(T) = DSM 26689(T) = LMG 28123(T) = BEI NR-42630(T)), Listeria grandensis sp. nov. (type strain TTU A1-0212(T) = FSL F6-0971(T) = DSM 26688(T) = LMG 28122(T) = BEI NR-42631(T)) and Listeria riparia sp. nov. (type strain FSL S10-1204(T) = DSM 26685(T) = LMG 28119(T) = BEI NR- 42634(T)) and (iv) a clade containing Listeria grayi. Genomic and phenotypic data suggest that the novel species are non-pathogenic.

Genomic comparison of sporeforming bacilli isolated from milk

BACKGROUND

Sporeformers in the order Bacillales are important contributors to spoilage of pasteurized milk. While only a few Bacillus and Viridibacillus strains can grow in milk at 6°C, the majority of Paenibacillus isolated from pasteurized fluid milk can grow under these conditions. To gain a better understanding of genomic features of these important spoilage organisms and to identify candidate genomic features that may facilitate cold growth in milk, we performed a comparative genomic analysis of selected dairy associated sporeformers representing isolates that can and cannot grow in milk at 6°C.

RESULTS

The genomes for seven Paenibacillus spp., two Bacillus spp., and one Viridibacillus sp. isolates were sequenced. Across the genomes sequenced, we identified numerous genes encoding antimicrobial resistance mechanisms, bacteriocins, and pathways for synthesis of non-ribosomal peptide antibiotics. Phylogenetic analysis placed genomes representing Bacillus, Paenibacillus and Viridibacillus into three distinct well supported clades and further classified the Paenibacillus strains characterized here into three distinct clades, including (i) clade I, which contains one strain able to grow at 6°C in skim milk broth and one strain not able to grow under these conditions, (ii) clade II, which contains three strains able to grow at 6°C in skim milk broth, and (iii) clade III, which contains two strains unable to grow under these conditions. While all Paenibacillus genomes were found to include multiple copies of genes encoding β-galactosidases, clade II strains showed significantly higher numbers of genes encoding these enzymes as compared to clade III strains. Genome comparison of strains able to grow at 6°C and strains unable to grow at this temperature identified numerous genes encoding features that might facilitate the growth of Paenibacillus in milk at 6°C, including peptidases with cold-adapted features (flexibility and disorder regions in the protein structure) and cold-adaptation related proteins (DEAD-box helicases, chaperone DnaJ).

CONCLUSIONS

Through a comparative genomics approach we identified a number of genomic features that may relate to the ability of selected Paenibacillus strains to cause spoilage of refrigerated fluid milk. With additional experimental evidence, these data will facilitate identification of targets to detect and control Gram positive spore formers in fluid milk.

Genome sequencing identifies Listeria fleischmannii subsp. coloradonensis subsp. nov., isolated from a ranch

Twenty Listeria -like isolates were obtained from environmental samples collected on a cattle ranch in northern Colorado; all of these isolates were found to share an identical partial sigB sequence, suggesting close relatedness. The isolates were similar to members of the genus Listeria in that they were Gram-stain-positive, short rods, oxidase-negative and catalase-positive; the isolates were similar to Listeria fleischmannii because they were non-motile at 25 °C. 16S rRNA gene sequencing for representative isolates and whole genome sequencing for one isolate was performed. The genome of the type strain of Listeria fleischmannii (strain LU2006-1T) was also sequenced. The draft genomes were very similar in size and the average MUMmer nucleotide identity across 91 % of the genomes was 95.16 %. Genome sequence data were used to design primers for a six-gene multi-locus sequence analysis (MLSA) scheme. Phylogenies based on (i) the near-complete 16S rRNA gene, (ii) 31 core genes and (iii) six housekeeping genes illustrated the close relationship of these Listeria -like isolates to Listeria fleischmannii LU2006-1T. Sufficient genetic divergence of the Listeria -like isolates from the type strain of Listeria fleischmannii and differing phenotypic characteristics warrant these isolates to be classified as members of a distinct infraspecific taxon, for which the name Listeria fleischmannii subsp. coloradonensis subsp. nov. is proposed. The type strain is TTU M1-001T ( = BAA-2414T = DSM 25391T). The isolates of Listeria fleischmannii subsp. coloradonensis subsp. nov. differ from the nominate subspecies by the inability to utilize melezitose, turanose and sucrose, and the ability to utilize inositol. The results also demonstrate the utility of whole genome sequencing to facilitate identification of novel taxa within a well-described genus. The genomes of both subspecies of Listeria fleischmannii contained putative enhancin genes; the Listeria fleischmannii subsp. coloradonensis subsp. nov. genome also encoded a putative mosquitocidal toxin. The presence of these genes suggests possible adaptation to an insect host, and further studies are needed to probe niche adaptation of Listeria fleischmannii .

Genomic characterization provides new insight into Salmonella phage diversity

Background

Salmonella is a widely distributed foodborne pathogen that causes tens of millions of salmonellosis cases globally every year. While the genomic diversity of Salmonella is increasingly well studied, our knowledge of Salmonella phage genomic diversity is still rather limited, despite the contributions of both lysogenic and lytic phages to Salmonella virulence, diversity and ecology (e.g., through horizontal gene transfer and Salmonella lysis). To gain a better understanding of phage diversity in a specific ecological niche, we sequenced 22 Salmonella phages isolated from a number of dairy farms from New York State (United States) and analyzed them using a comparative genomics approach.

Results

Classification of the 22 phages according to the presence/absence of orthologous genes allowed for classification into 8 well supported clusters. In addition to two phage clusters that represent novel virulent Salmonella phages, we also identified four phage clusters that each contained previously characterized phages from multiple continents. Our analyses also identified two clusters of phages that carry putative virulence (e.g., adhesins) and antimicrobial resistance (tellurite and bicyclomycin) genes as well as virulent and temperate transducing phages. Insights into phage evolution from our analyses include (i) identification of DNA metabolism genes that may facilitate nucleotide synthesis in phages with a G+C % distinct from Salmonella, and (ii) evidence of Salmonella phage tailspike and fiber diversity due to both single nucleotide polymorphisms and major re-arrangements, which may affect the host specificity of Salmonella phages.

Conclusions

Genomics-based characterization of 22 Salmonella phages isolated from dairy farms allowed for identification of a number of novel Salmonella phages. While the comparative genomics analyses of these phages provide a number of new insights in the evolution and diversity of Salmonella phages, they only represent a first glimpse into the diversity of Salmonella phages that is likely to be discovered when phages from different environments are characterized.

Evolutionary dynamics of the accessory genome of Listeria monocytogenes

Listeria monocytogenes, a foodborne bacterial pathogen, is comprised of four phylogenetic lineages that vary with regard to their serotypes and distribution among sources. In order to characterize lineage-specific genomic diversity within L. monocytogenes, we sequenced the genomes of eight strains from several lineages and serotypes, and characterized the accessory genome, which was hypothesized to contribute to phenotypic differences across lineages. The eight L. monocytogenes genomes sequenced range in size from 2.85-3.14 Mb, encode 2,822-3,187 genes, and include the first publicly available sequenced representatives of serotypes 1/2c, 3a and 4c. Mapping of the distribution of accessory genes revealed two distinct regions of the L. monocytogenes chromosome: an accessory-rich region in the first 65° adjacent to the origin of replication and a more stable region in the remaining 295°. This pattern of genome organization is distinct from that of related bacteria Staphylococcus aureus and Bacillus cereus. The accessory genome of all lineages is enriched for cell surface-related genes and phosphotransferase systems, and transcriptional regulators, highlighting the selective pressures faced by contemporary strains from their hosts, other microbes, and their environment. Phylogenetic analysis of O-antigen genes and gene clusters predicts that serotype 4 was ancestral in L. monocytogenes and serotype 1/2 associated gene clusters were putatively introduced through horizontal gene transfer in the ancestral population of L. monocytogenes lineage I and II.

Silage collected from dairy farms harbors an abundance of listeriaphages with considerable host range and genome size diversity

Since the food-borne pathogen Listeria monocytogenes is common in dairy farm environments, it is likely that phages infecting this bacterium ("listeriaphages") are abundant on dairy farms. To better understand the ecology and diversity of listeriaphages on dairy farms and to develop a diverse phage collection for further studies, silage samples collected on two dairy farms were screened for L. monocytogenes and listeriaphages. While only 4.5% of silage samples tested positive for L. monocytogenes, 47.8% of samples were positive for listeriaphages, containing up to >1.5 × 10(4) PFU/g. Host range characterization of the 114 phage isolates obtained, with a reference set of 13 L. monocytogenes strains representing the nine major serotypes and four lineages, revealed considerable host range diversity; phage isolates were classified into nine lysis groups. While one serotype 3c strain was not lysed by any phage isolates, serotype 4 strains were highly susceptible to phages and were lysed by 63.2 to 88.6% of phages tested. Overall, 12.3% of phage isolates showed a narrow host range (lysing 1 to 5 strains), while 28.9% of phages represented broad host range (lysing ≥11 strains). Genome sizes of the phage isolates were estimated to range from approximately 26 to 140 kb. The extensive host range and genomic diversity of phages observed here suggest an important role of phages in the ecology of L. monocytogenes on dairy farms. In addition, the phage collection developed here has the potential to facilitate further development of phage-based biocontrol strategies (e.g., in silage) and other phage-based tools.

Ribosomal RNA gene diversity, effective population size, and evolutionary longevity in asexual glomeromycota

Arbuscular mycorrhizal fungi (phylum Glomeromycota) are among the oldest and most successful symbionts of land plants. With no evidence of sexual reproduction, their evolutionary success is inconsistent with the prediction that asexual taxa are vulnerable to extinction due to accumulation of deleterious mutations. To explore why Glomeromycota defy this prediction, we studied ribosomal RNA (rRNA) gene evolution in the Claroideoglomus lineage and estimated effective population size, N(e) , in C. etunicatum. We found that rRNA genes of these fungi exhibit unusual and complex patterns of molecular evolution. In C. etunicatum, these patterns can be collectively explained by an unexpectedly large N(e) combined with imperfect genome-wide and population-level rRNA gene repeat homogenization. The mutations accumulated in rRNA gene sequences indicate that natural selection is effective at purging deleterious mutations in the Claroideoglomus lineage, which is also consistent with the large N(e) of C. etunicatum. We propose that in the near absence of recombination, asexual reproduction involving massively multinucleate spores typical for Glomeromycota is responsible for the improved efficacy of selection relative to drift. We postulate that large effective population sizes contribute to the evolutionary longevity of Glomeromycota.

Genome sequencing reveals diversification of virulence factor content and possible host adaptation in distinct subpopulations of Salmonella enterica

Background

Divergence of bacterial populations into distinct subpopulations is often the result of ecological isolation. While some studies have suggested the existence of Salmonella enterica subsp. enterica subclades, evidence for these subdivisions has been ambiguous. Here we used a comparative genomics approach to define the population structure of Salmonella enterica subsp. enterica, and identify clade-specific genes that may be the result of ecological specialization.

Results

Multi-locus sequence analysis (MLSA) and single nucleotide polymorphisms (SNPs) data for 16 newly sequenced and 30 publicly available genomes showed an unambiguous subdivision of S. enterica subsp. enterica into at least two subpopulations, which we refer to as clade A and clade B. Clade B strains contain several clade-specific genes or operons, including a β-glucuronidase operon, a S-fimbrial operon, and cell surface related genes, which strongly suggests niche specialization of this subpopulation. An additional set of 123 isolates was assigned to clades A and B by using qPCR assays targeting subpopulation-specific SNPs and genes of interest. Among 98 serovars examined, approximately 20% belonged to clade B. All clade B isolates contained two pathogenicity related genomic islands, SPI-18 and a cytolethal distending toxin islet; a combination of these two islands was previously thought to be exclusive to serovars Typhi and Paratyphi A. Presence of β-glucuronidase in clade B isolates specifically suggests an adaptation of this clade to the vertebrate gastrointestinal environment.

Conclusions

S. enterica subsp. enterica consists of at least two subpopulations that differ specifically in genes involved in host and tissue tropism, utilization of host specific carbon and nitrogen sources and are therefore likely to differ in ecology and transmission characteristics.

Listeriosis Outbreak in Dairy Cattle Caused by an Unusual Listeria Monocytogenes Serotype 4b Strain

A listeriosis outbreak, in dairy cattle, with a high case mortality and acute death after onset of symptoms was investigated using gross pathology and bacteriologic approaches, including molecular characterization of a clinical Listeria monocytogenes isolate. In a herd of 315 animals, 9 animals showed clinical symptoms consistent with listeriosis, including 3 animals that died within 2–4 days after acute onset of clinical signs, 4 animals that were euthanized, and 2 that survived. Initial EcoRI ribotyping and serotyping indicated that this outbreak was caused by an unusual L. monocytogenes serotype 4b strain, which was classified into lineage III. Further characterization of this isolate by DNA sequencing–based subtyping methods indicated that the strain responsible for this outbreak represented a unique genotype as supported by its classification into a new sigB allelic type, which has not been identified previously among >290 isolates, and by compelling phylogenetic evidence. While lineage III isolates are generally rare, they seem to be more common among L. monocytogenes isolates from animals with clinical signs of listeriosis. This is the first report of a particularly severe clinical course of disease associated with infection by a lineage III strain. The high prevalence of Listeria spp., including L. monocytogenes, in the farm environments may favor emergence and evolution of novel, and possibly more virulent, L. monocytogenes strains. Continued monitoring of animal listeriosis cases and outbreaks may not only improve animal health but also aid in the early discovery of newly emerging L. monocytogenes strains.

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

Background

The 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.

Results

To 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.

Conclusions

Genome 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.

Salt stress phenotypes in Listeria monocytogenes vary by genetic lineage and temperature

Listeria monocytogenes can survive and grow under wide-ranging environmental stress conditions encountered both in foods and in the host. The ability of certain L. monocytogenes subtypes to thrive under stress conditions present in specific niches was hypothesized to reflect genetic characteristics and phenotypic capabilities conserved among strains within a subtype. To quantify variations in salt stress phenotypes among 40 strains selected to represent the diversity of the three major L. monocytogenes genetic lineages and to determine if salt stress phenotypes were associated with genetic relatedness, we measured growth under salt stress at both 7°C and 37°C. At 7°C, in brain-heart infusion with 6% NaCl, average growth rates among the lineages were similar. A comparison of doubling times after exposure to salt stress at 7°C or 37°C indicated that growth at 7°C provided crossprotection to subsequent salt stress for strains in lineages I and II. At 37°C, in brain-heart infusion with 6% NaCl, lineage I and III strains grew significantly faster (p<0.0001) than lineage II strains. Under salt stress at 37°C, differences in growth parameters were significantly (p<0.005) associated with genetic relatedness of the strains. Compatible solute uptake is part of the L. monocytogenes salt stress response, but growth differences between the lineages were not related to differences in transcript levels of osmolyte transporter-encoding genes betL, gbuA, oppA, and opuCA. The combination of phylogenetic and phenotypic data suggests that L. monocytogenes lineage I and III strains, which are most commonly associated with human and animal disease, may be better adapted to osmotic stress at 37°C, conditions that are present in the host gastrointestinal tract.

A population genetics-based and phylogenetic approach to understanding the evolution of virulence in the genus Listeria

The genus Listeria includes (i) the opportunistic pathogens L. monocytogenes and L. ivanovii, (ii) the saprotrophs L. innocua, L. marthii, and L. welshimeri, and (iii) L. seeligeri, an apparent saprotroph that nevertheless typically contains the prfA virulence gene cluster. A novel 10-loci multilocus sequence typing scheme was developed and used to characterize 67 isolates representing six Listeria spp. (excluding L. grayi) in order to (i) provide an improved understanding of the phylogeny and evolution of the genus Listeria and (ii) use Listeria as a model to study the evolution of pathogenicity in opportunistic environmental pathogens. Phylogenetic analyses identified six well-supported Listeria species that group into two main subdivisions, with each subdivision containing strains with and without the prfA virulence gene cluster. Stochastic character mapping and phylogenetic analysis of hly, a gene in the prfA cluster, suggest that the common ancestor of the genus Listeria contained the prfA virulence gene cluster and that this cluster was lost at least five times during the evolution of Listeria, yielding multiple distinct saprotrophic clades. L. welshimeri, which appears to represent the most ancient clade that arose from an ancestor with a prfA cluster deletion, shows a considerably lower average sequence divergence than other Listeria species, suggesting a population bottleneck and a putatively different ecology than other saprotrophic Listeria species. Overall, our data suggest that, for some pathogens, loss of virulence genes may represent a selective advantage, possibly by facilitating adaptation to a specific ecological niche.

Listeria monocytogenes lineages: Genomics, evolution, ecology, and phenotypic characteristics

Listeria monocytogenes consists of at least 4 evolutionary lineages (I, II, III, and IV) with different but overlapping ecological niches. Most L. monocytogenes isolates seem to belong to lineages I and II, which harbor the serotypes more commonly associated with human clinical cases, including serotype 1/2a (lineage II) and serotypes 1/2b and 4b (lineage I). Lineage II strains are common in foods, seem to be widespread in the natural and farm environments, and are also commonly isolated from animal listeriosis cases and sporadic human clinical cases. Most human listeriosis outbreaks are associated with lineage I isolates though. In addition, a number of studies indicate that, in many countries, lineage I strains are overrepresented among human isolates, as compared to lineage II strains. Lineage III and IV strains on the other hand are rare and predominantly isolated from animal sources. The apparent differences in the distribution of strains representing the L. monocytogenes lineages has lead to a number of studies aimed at identifying phenotypic differences among the different lineages. Interestingly, lineage II isolates seem to carry more plasmids than lineage I isolates and these plasmids often confer resistance to toxic metals and possibly other compounds that may be found in the environment. Moreover, lineage II isolates seem to be more resistant to bacteriocins than lineage I isolates, which probably confers an advantage in environments where bacteriocin-producing organisms are abundant. A large number of lineage II isolates and strains have been shown to be virulence-attenuated due to premature stop codon mutations in inlA and mutations in prfA. A subset of lineage I isolates carry a listeriolysin S hemolysin, which is not present in isolates belonging to lineages II, III, or IV. While lineage II isolates also show higher recombination rates than lineage I isolates, possibly facilitating adaptation of lineage II strains to diverse environments, lineage I isolates are clonal and show a low prevalence of plasmids and IS elements, suggesting that lineage I isolates may have mechanisms that limit the acquisition of foreign DNA by horizontal gene transfer. Diversifying selection has also been shown to have played an important role during evolution of the L. monocytogenes lineages and during divergence of L. monocytogenes from the non-pathogenic species L. innocua. Overall evidence thus suggests that the 4 L. monocytogenes lineages identified so far represent distinct ecologic, genetic, and phenotypic characteristics, which appear to affect their ability to be transmitted through foods and to cause human disease. Further insights into the ecology, evolution, and characteristics of these lineages will thus not only provide an improved understanding of the evolution of this foodborne pathogen, but may also facilitate improved control of foodborne listeriosis.

Multilocus sequence typing of outbreak-associated Listeria monocytogenes isolates to identify epidemic clones

OBJECTIVE

Listeria monocytogenes is a foodborne pathogen found in a wide variety of environments. Subtype characterization of L. monocytogenes isolates from listeriosis outbreaks that have occurred over the last three decades has suggested that a number of outbreaks were caused by a small number of L. monocytogenes epidemic clones (ECs). In this study we compared the prevalence, ecology, and phylogenetic position of outbreak-associated isolates and non-outbreak-associated isolates to probe the evolutionary and ecological characteristics of outbreak-associated L. monocytogenes subtypes, including those representing previously described ECs.

METHODS

Multilocus sequence typing data for isolates from 15 listeriosis outbreaks in Europe and North America were generated and compared, using a phylogenetic framework, with 180 isolates representing a local sampling of diverse sources, including human sporadic cases.

RESULTS

Isolates from 15 listeriosis outbreaks represented eight sequence types (STs). STs corresponding to previously designated ECI (ST1 and ST93) and ECIa (ST29) represented isolates from eight outbreaks. ST17 (corresponding to ECII) was involved in two outbreaks in the United States (1998 and 2002). No other STs were involved in multiple outbreaks. While ST1 was the most common ST among sporadic human cases and non-human listeriosis-related isolates, ST29 was rare among non-human listeriosis-related isolates and was significantly overrepresented among isolates from human listeriosis outbreaks and sporadic cases as compared to isolates from other sources in our local sampling.

CONCLUSIONS

STs associated with outbreaks (and representing previously designated ECs) appear to differ in their ecology. While association of ECI with multiple human listeriosis outbreaks appears to reflect strain abundance across environments, ECIa seems to represent an L. monocytogenes EC that appears to be overrepresented among outbreaks and sporadic cases and thus may have increased transmission potential.

Listeria marthii sp. nov., isolated from the natural environment, Finger Lakes National Forest

Four isolates (FSL S4-120(T), FSL S4-696, FSL S4-710, and FSL S4-965) of Gram-positive, motile, facultatively anaerobic, non-spore-forming bacilli that were phenotypically similar to species of the genus Listeria were isolated from soil, standing water and flowing water samples obtained from the natural environment in the Finger Lakes National Forest, New York, USA. The four isolates were closely related to one another and were determined to be the same species by whole genome DNA-DNA hybridization studies (>82 % relatedness at 55 degrees C and >76 % relatedness at 70 degrees C with 0.0-0.5 % divergence). 16S rRNA gene sequence analysis confirmed their close phylogenetic relatedness to Listeria monocytogenes and Listeria innocua and more distant relatedness to Listeria welshimeri, L. seeligeri, L. ivanovii and L. grayi. Phylogenetic analysis of partial sequences for sigB, gap, and prs showed that these isolates form a well-supported sistergroup to L. monocytogenes. The four isolates were sufficiently different from L. monocytogenes and L. innocua by DNA-DNA hybridization to warrant their designation as a new species of the genus Listeria. The four isolates yielded positive reactions in the AccuProbe test that is purported to be specific for L. monocytogenes, did not ferment L-rhamnose, were non-haemolytic on blood agar media, and did not contain a homologue of the L. monocytogenes virulence gene island. On the basis of their phenotypic characteristics and their genotypic distinctiveness from L. monocytogenes and L. innocua, the four isolates should be classified as a new species within the genus Listeria, for which the name Listeria marthii sp. nov. is proposed. The type strain of L. marthii is FSL S4-120(T) (=ATCC BAA-1595(T) =BEIR NR 9579(T) =CCUG 56148(T)). L. marthii has not been associated with human or animal disease at this time.

Lineage specific recombination rates and microevolution in Listeria monocytogenes

Background

The bacterium Listeria monocytogenes is a saprotroph as well as an opportunistic human foodborne pathogen, which has previously been shown to consist of at least two widespread lineages (termed lineages I and II) and an uncommon lineage (lineage III). While some L. monocytogenes strains show evidence for considerable diversification by homologous recombination, our understanding of the contribution of recombination to L. monocytogenes evolution is still limited. We therefore used STRUCTURE and ClonalFrame, two programs that model the effect of recombination, to make inferences about the population structure and different aspects of the recombination process in L. monocytogenes. Analyses were performed using sequences for seven loci (including the house-keeping genes gap, prs, purM and ribC, the stress response gene sigB, and the virulence genes actA and inlA) for 195 L. monocytogenes isolates.

Results

Sequence analyses with ClonalFrame and the Sawyer's test showed that recombination is more prevalent in lineage II than lineage I and is most frequent in two house-keeping genes (ribC and purM) and the two virulence genes (actA and inlA). The relative occurrence of recombination versus point mutation is about six times higher in lineage II than in lineage I, which causes a higher genetic variability in lineage II. Unlike lineage I, lineage II represents a genetically heterogeneous population with a relatively high proportion (30% average) of genetic material imported from external sources. Phylograms, constructed with correcting for recombination, as well as Tajima's D data suggest that both lineages I and II have suffered a population bottleneck.

Conclusion

Our study shows that evolutionary lineages within a single bacterial species can differ considerably in the relative contributions of recombination to genetic diversification. Accounting for recombination in phylogenetic studies is critical, and new evolutionary models that account for the possibility of changes in the rate of recombination would be required. While previous studies suggested that only L. monocytogenes lineage I has experienced a recent bottleneck, our analyses clearly show that lineage II experienced a bottleneck at about the same time, which was subsequently obscured by abundant homologous recombination after the lineage II bottleneck. While lineage I and lineage II should be considered separate species from an evolutionary viewpoint, maintaining single species name may be warranted since both lineages cause the same type of human disease.

Phylogeographic patterns in Leccinum sect. Scabra and the status of the arctic-alpine species L. rotundifoliae

We investigated inter- and intraspecific phylogenetic relationships in the ectomycorrhizal fungal genus Leccinum section Scabra. Species of this section are exclusively associated with Betula and occur throughout the Northern Hemisphere. We compared the phylogenetic relationships of arctic, alpine, boreal and temperate accessions of section Scabra based on DNA sequences of the single-copy nuclear gene Gapdh and the multiple-copy nuclear region 5.8S-ITS2. Exclusively arctic lineages were not detected in species that occur both in arctic-alpine or boreal regions, except in L. rotundifoliae that was restricted to cold climates. L. scabrum and L. holopus showed an intercontinental phylogeographic pattern, and L. variicolor showed a pattern unrelated to geographical distribution. Molecular clock estimates indicated that L. rotundifoliae is as old as other species in section Scabra. Individual gene trees suggest that interspecific hybridisation occurred several times in the evolution of section Scabra.

Morphological and molecular evidence supporting an arbutoid mycorrhizal relationship in the Costa Rican páramo

This study examines evidence for a particular arbutoid mycorrhizal interaction in páramo, a high-altitude neotropical ecosystem important in hydrological regulation but poorly known in terms of its fungal communities. Comarostaphylis arbutoides Lindley (Ericaceae) often forms dense thickets in Central American páramo habitats. Based on phylogenetic classification, it has been suggested that C. arbutoides forms arbutoid mycorrhizae with diverse Basidiomycetes and Ascomycetes; however, this assumption has not previously been confirmed. Based on field data, we hypothesized an arbutoid mycorrhizal association between C. arbutoides and the recently described bolete Leccinum monticola Halling & G.M. Mueller; in this study, we applied a rigorous approach using anatomical and molecular data to examine evidence for such an association. We examined root samples collected beneath L. monticola basidiomes for mycorrhizal structures, and we also compared rDNA internal transcribed spacer (ITS) sequences between mycorrhizal root tips and leaf or basidiome material of the suspected symbionts. Root cross sections showed a thin hyphal sheath and intracellular hyphal coils typical of arbutoid mycorrhizae. DNA sequence comparisons confirmed the identity of C. arbutoides and L. monticola as the mycorrhizal symbionts. In addition, this paper provides additional evidence for the widespread presence of minisatellite-like inserts in the ITS1 spacer in Leccinum species (including a characterization of the insert in L. monticola) and reports the use of an angiosperm-specific ITS primer pair useful for amplifying plant DNA from mycorrhizal roots without co-amplifying fungal DNA.

An ITS phylogeny of Leccinum and an analysis of the evolution of minisatellite-like sequences within ITS1

Phylogenetic relationships of the European species of Leccinum (Boletales, Boletaceae) were investigated by maximum parsimony, Bayesian and likelihood analyses of nrITS1-5.8S-ITS2 and 28S sequences. The separate gene trees inferred were largely concordant, and their combined analysis indicates that several traditional sectional and species-level taxonomic schemes are artificial. In Leccinum, the nrITS region ranges in size from 694 to 1480 bp. This extreme length heterogeneity is localized to a part of the ITS1 spacer that contains a minisatellite characterized by the repeated presence of CTATTGAAAAG and CTAATAGAAAG core sequences and mutational derivatives thereof. The number of core sequences present in the minisatellite varied from 12 to 36. Intra-individual sequence variation of the minisatellite was always smaller than between different species, indicating that concerted evolution proceeds rapidly enough to retain phylogenetic signal at the infraspecific level. In contrast, the evolutionary pattern exhibited by the major ITS1 repeat types found was homoplastic when mapped onto the species lineages inferred from the combined 5.8S-ITS2 sequences. The minisatellite therefore appears not to be useful for phylogeny reconstruction at or above the species level.