Search

Raw reads of the 34 genomes of S. Bovismorbificans and 12 genomes of S. Typhimurium from human and food origin in this study were submitted to the European Nucleotide Archive under project numbers PRJEB41645 and PRJEB49565.; International audience; The French National Reference Centre for Escherichia coli, Shigella and Salmonella (FNRC-ESS) detected two human clusters of 33 cases (median age: 10 years; 17 females) infected by Salmonella enterica serotype Bovismorbificans, ST142, HC5_243255 (EnteroBase HierCC‑cgMLST scheme) in September-November 2020 and of 11 cases (median age: 11 years; seven males) infected by S. enterica serotype 4,12:i:-, ST34, HC5_198125 in October-December 2020. Epidemiological investigations conducted by Santé publique France linked these outbreaks to the consumption of dried pork sausages from the same manufacturer. S. Bovismorbificans and S. 4,12:i:- were isolated by the National Reference Laboratory from different food samples, but both strains were identified in a single food sample only by qPCR. Three recalls and withdrawals of dried pork products were issued by the French general directorate of food of the French ministry for agriculture and food in November 2020, affecting eight supermarket chains. A notification on the European Rapid Alert System for Food and Feed and a European urgent enquiry on the Epidemic Intelligence Information System for Food and Waterborne Diseases and Zoonoses (EPIS-FWD) were launched. No cases were reported outside France. Outbreaks caused by multiple serotypes of Salmonella may go undetected by protocols in standard procedures in microbiology laboratories.

With the advent of next-generation whole-genome sequencing (WGS), the need for good-quality and well-characterised

The European epidemic monophasic variant of Salmonella enterica serovar Typhimurium (S. 1,4,[5],12:i:-) characterized by the multi locus sequence type ST34 and the antimicrobial resistance ASSuT profile has become one of the most common serovars in Europe (EU) and the United States (US). In this study, we reconstructed the time-scaled phylogeny and evolution of this Salmonella in Europe. The epidemic S. 1,4,[5],12:i:- ST34 emerged in the 1980s by an acquisition of the Salmonella Genomic Island (SGI)-4 at the 3′ end of the phenylalanine phe tRNA locus conferring resistance to copper and arsenic toxicity. Subsequent integration of the Tn21 transposon into the fljAB locus gave resistance to mercury toxicity and several classes of antibiotics used in food-producing animals (ASSuT profile). The second step of the evolution occurred in the 1990s, with the integration of mTmV and mTmV-like prophages carrying the perC and/or sopE genes involved in the ability to reduce nitrates in intestinal contents and facilitate the disruption of the junctions of the host intestinal epithelial cells. Heavy metals are largely used as food supplements or pesticide for cultivation of seeds intended for animal feed so the expansion of the epidemic S. 1,4,[5],12:i:- ST34 was strongly related to the multiple-heavy metal resistance acquired by transposons, integrative and conjugative elements and facilitated by the escape until 2011 from the regulatory actions applied in the control of S. Typhimurium in Europe. The genomic plasticity of the epidemic S. 1,4,[5],12:i:- was demonstrated in our study by the analysis of the plasmidome. We were able to identify plasmids harboring genes mediating resistance to phenicols, colistin, and fluoroquinolone and also describe for the first time in six of the analyzed genomes the presence of two plasmids (pERR1744967-1 and pERR2174855-2) previously described only in strains of enterotoxigenic Escherichia coli and E. fergusonii.

The partitioning of pathogenic strains isolated in environmental or human cases to their sources is challenging. The pathogens usually colonize multiple animal hosts, including livestock, which contaminate the food-production chain and the environment (e.g. soil and water), posing an additional public-health burden and major challenges in the identification of the source. Genomic data opens up new opportunities for the development of statistical models aiming to indicate the likely source of pathogen contamination. Here, we propose a computationally fast and efficient multinomial logistic regression source-attribution classifier to predict the animal source of bacterial isolates based on ‘source-enriched’ loci extracted from the accessory-genome profiles of a pangenomic dataset. Depending on the accuracy of the model’s self-attribution step, the modeller selects the number of candidate accessory genes that best fit the model for calculating the likelihood of (source) category membership. The Accessory genes-Based Source Attribution (AB_SA) method was applied to a dataset of strains ofSalmonella entericaTyphimurium and its monophasic variant (S.enterica1,4,[5],12:i:-). The model was trained on 69 strains with known animal-source categories (i.e. poultry, ruminant and pig). The AB_SA method helped to identify 8 genes as predictors among the 2802 accessory genes. The self-attribution accuracy was 80 %. The AB_SA model was then able to classify 25 of the 29S.entericaTyphimurium andS.enterica1,4,[5],12:i:- isolates collected from the environment (considered to be of unknown source) into a specific category (i.e. animal source), with more than 85 % of probability. The AB_SA method herein described provides a user-friendly and valuable tool for performing source-attribution studies in only a few steps. AB_SA is written in R and freely available athttps://github.com/lguillier/AB_SA.

Volume 8, no. 1, e01361-18, 2019, [https://doi.org/10.1128/MRA.01361-18][1]. Page 1: The article title should read as given above. Page 1, abstract, line 2: “St. Nectaire and Morbier cheeses” should read “two raw milk cheeses.” Page 1, abstract, lines 3 and 4: “fatal outbreaks” should

The partitioning of pathogenic strains isolated in environmental or human cases to their original source is challenging. The pathogens usually colonize multiple animal hosts, including livestock, which contaminate food-producing and environment (e.g. soil and water), posing additional public health burden and major challenges in the identification of the source. Genomic data opens new opportunities for the development of statistical models aiming to infer the likely source of pathogen contamination. Here, we propose a computationally fast and efficient multinomial logistic regression (MLR) source attribution classifier to predict the animal source of bacterial isolates based on “source-enriched” loci extracted from the accessory-genome profiles of a pangenomic dataset. Depending on the accuracy of the model’s self-attribution step, the modeler selects the number of candidate accessory genes that better fit the model for calculating the likelihood of (source) category membership. The accessory genes-based source attribution (AB_SA) method was applied on a dataset of strains of Salmonella Typhimurium and its monophasic variants (S. 1,4,[5],12:i:-). The model was trained on 69 strains with known animal source categories (i.e., poultry, ruminant, and pig). The AB_SA method helped to identify eight genes as predictors among the 2,802 accessory genes. The self-attribution accuracy was 80%. The AB_SA model was then able to classify 25 over 29 S. Typhimurium and S. 1,4,[5],12:i:-isolates collected from the environment (considered as unknown source) into a specific category (i.e., animal source), with more than 85% of probability. The AB_SA method herein described provides a user-friendly and valuable tool to perform source attribution studies in few steps. AB_SA is written in R and freely available at https://github.com/lguillier/AB_SA.Author NotesAll supporting data, code, and protocols have been provided within the article and through supplementary data files.Supplementary material is available with the online version of this article.AbbreviationsAB_SA, accessory-based source attribution; MLR, multinomial logistic regression; SNPs, single nucleotide polymorphisms; GFF, general feature format; AIC, Akaike information criteria.Data SummaryThe AB_SA model is written in R, open-source and freely available Github under the GNU GPLv3 licence (https://github.com/lguillier/AB_SA).All sequencing reads used to generate the assemblies analyzed in this study have been deposited in the European Nucleotide Archive (ENA) (http://www.ebi.ac.uk/ena) under project number PRJEB16326. Genome metadata and ENA run accession ID for all the assemblies are reported in the supplementary material.Impact StatementThis article describes AB_SA (“Accessory-Based Source Attribution method”), a novel approach for source attribution based on “source enriched” accessory genomics data and unsupervised multinomial logistic regression. We demonstrate that the AB_SA method enables the animal source prediction of large-scale datasets of bacterial populations through rapid and easy identification of source predictors from the non-core genomic regions. Herein, AB_SA correctly self-attribute the animal source of a set of S. Typhimurium and S. 1,4,[5],12:i:- isolates and further classifies the 84% of strains contaminating natural environments in the pig category (with high probability ranging between ∼85 and ∼99%).

International audience; The investigation of foodborne outbreaks (FBOs) from genomic data typically relies on inspecting the relatedness of samples through a phylogenomic tree computed on either SNPs, genes, kmers, or alleles (i.e., cgMLST and wgMLST). The phylogenomic reconstruction is often time-consuming, computation-intensive and depends on hidden assumptions, pipelines implementation and their parameterization. In the context of FBO investigations, robust links between isolates are required in a timely manner to trigger appropriate management actions. Here, we propose a non-parametric statistical method to assert the relatedness of samples (i.e., outbreak cases) or whether to reject them (i.e., non-outbreak cases). With typical computation running within minutes on a desktop computer, we benchmarked the ability of three non-parametric statistical tests (i.e., Wilcoxon rank-sum, Kolmogorov-Smirnov and Kruskal-Wallis) on six different genomic features (i.e., SNPs, SNPs excluding recombination events, genes, kmers, cgMLST alleles, and wgMLST alleles) to discriminate outbreak cases (i.e., positive control: C+) from non-outbreak cases (i.e., negative control: C−). We leveraged four well-characterized and retrospectively investigated FBOs of Salmonella Typhimurium and its monophasic variant S. 1,4,[5],12:i:-from France, setting positive and negative controls in all the assays. We show that the approaches relying on pairwise SNP differences distinguished all four considered outbreaks in contrast to the other tested genomic features (i.e., genes, kmers, cgMLST alleles, and wgMLST alleles). The freely available non-parametric method written in R has been designed to be independent of both the phylogenomic reconstruction and the detection methods of genomic features (i.e., SNPs, genes, kmers, or alleles), making it widely and easily usable to anybody working on genomic data from suspected samples.

Salmonella Derby (S. Derby) is emerging in Europe as a predominant serovar in fattening turkey flocks. This serovar was recorded as being predominant in the turkey sector in 2014 in the United Kingdom (UK). Only two years later, in 2016, it was also recorded in the turkey and broiler sectors in Ireland and Spain. These S. Derby isolates were characterised as members of the multilocus sequence type (MLST) profile 71 (ST71). For the first time, we characterise by whole genome sequencing (WGS) analysis a panel of 90 S. Derby ST71 genomes to understand the routes of transmission of this emerging pathogen within the poultry/turkey food trade. Selected panel included strains isolated as early as 2010 in five leading European g countries for turkey meat production. Twenty-one of the 90 genomes were extracted from a public database&mdash, Enterobase. Five of these originated from the United States (n=3), China (n=1) and Taiwan (n=1) isolated between 1986 and 2016. A phylogenomic analysis at the core-genome level revealed the presence of three groups. The largest group contained 97.5% of the European strains and included both, turkey and human isolates that were genetically related by an average of 35 &plusmn, 15 single nucleotide polymorphism substitutions (SNPs). To illustrate the diversity, the presence of antimicrobial resistance genes and phages were characteised in 30, S. Derby ST71 genomes, including 11 belonging to this study This study revealed an emergent turkey-related S. Derby ST71 clone circulating in at least five European countries (the UK, Germany, Poland, Italy, and France) since 2010 that causes human gastroenteritis. A matter of concern is the identification of a gyrA mutation involved in resistance to quinolone, present in the Italian genomes. Interestingly, the diversity of phages seems to be related to the geographic origins. These results constitute a baseline for following the spread of this emerging pathogen and identifying appropriate monitoring and prevention measures.

Here, we report the genome sequences of two Salmonella enterica serovar Dublin strains, 03EB8736SAL and 03EB8994SAL, isolated from raw-milk cheese and milk filtrate, respectively. Analysis of the draft genomes of the two isolates reveals the presence of the viaA, viaB, and ompB loci of the Vi capsular polysaccharide antigen (Vi antigen).

Here, we report the genome sequences of two Salmonella enterica serovar Dublin strains, 03EB8736SAL and 03EB8994SAL, isolated from raw-milk cheese and milk filtrate, respectively. Analysis of the draft genomes of the two isolates reveals the presence of the viaA , viaB, and ompB loci of the Vi capsular polysaccharide antigen (Vi antigen).

On 18 January 2016, the French National Reference Centre for Salmonella reported to Santé publique France an excess of Salmonella enterica serotype Dublin (S. Dublin) infections. We investigated to identify the source of infection and implement control measures. Whole genome sequencing (WGS) and multilocus variable-number tandem repeat analysis (MLVA) were performed to identify microbiological clusters and links among cases, animal and food sources. Clusters were defined as isolates with less than 15 single nucleotide polymorphisms determined by WGS and/or with identical MLVA pattern. We compared different clusters of cases with other cases (case–case study) and controls recruited from a web-based cohort (case–control study) in terms of food consumption. We interviewed 63/83 (76%) cases; 2,914 controls completed a questionnaire. Both studies’ findings indicated that successive S. Dublin outbreaks from different sources had occurred between November 2015 and March 2016. In the case–control study, cases of distinct WGS clusters were more likely to have consumed Morbier (adjusted odds ratio (aOR): 14; 95% confidence interval (CI): 4.8–42) or Vacherin Mont d’Or (aOR: 27; 95% CI: 6.8–105), two bovine raw-milk cheeses. Based on these results, the Ministry of Agriculture launched a reinforced control plan for processing plants of raw-milk cheeses in the production region, to prevent future outbreaks.

On 18 January 2016, the French National Reference Centre for Salmonella reported to Santé publique France an excess of Salmonella enterica serotype Dublin (S. Dublin) infections. We investigated to identify the source of infection and implement control measures. Whole genome sequencing (WGS) and multilocus variable-number tandem repeat analysis (MLVA) were performed to identify microbiological clusters and links among cases, animal and food sources. Clusters were defined as isolates with less than 15 single nucleotide polymorphisms determined by WGS and/or with identical MLVA pattern. We compared different clusters of cases with other cases (case–case study) and controls recruited from a web-based cohort (case–control study) in terms of food consumption. We interviewed 63/83 (76%) cases; 2,914 controls completed a questionnaire. Both studies’ findings indicated that successive S. Dublin outbreaks from different sources had occurred between November 2015 and March 2016. In the case–control study, cases of distinct WGS clusters were more likely to have consumed Morbier (adjusted odds ratio (aOR): 14; 95% confidence interval (CI): 4.8–42) or Vacherin Mont d’Or (aOR: 27; 95% CI: 6.8–105), two bovine raw-milk cheeses. Based on these results, the Ministry of Agriculture launched a reinforced control plan for processing plants of raw-milk cheeses in the production region, to prevent future outbreaks.

In the European Union, Salmonella enterica subsp. enterica serovar Derby is the most abundant serotype isolated from pork., In the European Union, Salmonella enterica subsp. enterica serovar Derby is the most abundant serotype isolated from pork. Recent studies have shown that this serotype is polyphyletic. However, one main genomic lineage, characterized by sequence type 40 (ST40), the presence of the Salmonella pathogenicity island 23, and showing resistance to streptomycin, sulphonamides, and tetracycline (STR-SSS-TET), is pork associated. Here, we describe the complete genome sequence of a strain from this lineage isolated in France.

We report here the draft genome sequences of 2 Salmonella enterica subsp. enterica serovar Dublin strains from St. Nectaire and Morbier cheeses having multilocus variable-number tandem-repeat analysis (MLVA) profiles identified during the fatal outbreaks that occurred in France in 2012 and 2015 to 2016, respectively., We report here the draft genome sequences of 2 Salmonella enterica subsp. enterica serovar Dublin strains from St. Nectaire and Morbier cheeses having multilocus variable-number tandem-repeat analysis (MLVA) profiles identified during the fatal outbreaks that occurred in France in 2012 and 2015 to 2016, respectively. These draft genome sequences will help uncover the virulence determinants in invasive S. Dublin strains.

In France, Salmonella Derby is one of the most prevalent serotypes in pork and poultry meat. Since 2006, it has ranked among the 10 most frequent Salmonella serotypes isolated in humans. In previous publications, Salmonella Derby isolates have been characterized by pulsed field gel electrophoresis (PFGE) and antimicrobial resistance (AMR) profiles revealing the existence of different pulsotypes and AMR phenotypic groups. However, these results suffer from the low discriminatory power of these typing methods. In the present study, we built a collection of 140 strains of S. Derby collected in France from 2014 to 2015 representative of the pork and poultry food sectors. The whole collection was characterized using whole genome sequencing (WGS), providing a significant contribution to the knowledge of this underrepresented serotype, with few genomes available in public databases. The genetic diversity of the S. Derby strains was analyzed by single-nucleotide polymorphism (SNP). We also investigated AMR by both genome and phenotype, the main Salmonella pathogenicity island (SPI) and the fimH gene sequences. Our results show that this S. Derby collection is spread across four different lineages genetically distant by an average of 15k SNPs. These lineages correspond to four multilocus sequence typing (MLST) types (ST39, ST40, ST71, and ST682), which were found to be associated with specific animal hosts: pork and poultry. While the ST71 and ST682 strains are pansusceptible, ST40 isolates are characterized by the multidrug resistant profile STR-SSS-TET. Considering virulence determinants, only ST39 and ST40 present the SPI-23, which has previously been associated with pork enterocyte invasion. Furthermore, the pork ST682 isolates were found to carry mutations in the fimH sequence that could participate in the host tropism of this group. Our phylogenetic analysis demonstrates the polyphyletic nature of the Salmonella serotype Derby and provides an opportunity to identify genetic factors associated with host adaptation and markers for the monitoring of these different lineages within the corresponding animal sectors. The recognition of these four lineages is of primary importance for epidemiological surveillance throughout the food production chains and constitutes the first step toward refining monitoring and preventing dispersal of this pathogen.

International audience; In the European Union, Salmonella enterica subsp. enterica serovar Derby is the most abundant serotype isolated from pork. Recent studies have shown that this serotype is polyphyletic. However, one main genomic lineage, characterized by sequence type 40 (ST40), the presence of the Salmonella pathogenicity island 23, and showing resistance to streptomycin, sulphonamides, and tetracycline (STR-SSS-TET), is pork associated. Here, we describe the complete genome sequence of a strain from this lineage isolated in France.

Salmonella enterica subspecies enterica serovar Dublin (S. Dublin) figures among the most frequently isolated Salmonella strains in humans in France. This serovar may affect production and animal health mainly in cattle herds with corresponding high economic losses. Given that the current gold standard method, pulsed-field gel electrophoresis (PFGE), provides insufficient discrimination for epidemiological investigations, we propose a standard operating procedure in this study for multiple-locus variable number tandem repeat analysis (MLVA) of S. Dublin, suitable for inter-laboratory surveillance. An in silico analysis on the genome of S. Dublin strains CT_02021853 was performed to identify appropriate microsatellite regions. Of 21 VNTR loci screened, six were selected and 401 epidemiologically unrelated and related strains, isolated from humans, food and animals were analyzed to assess performance criteria such as typeability, discriminatory power and epidemiological concordance. The MLVA scheme developed was applied to an outbreak involving Saint-Nectaire cheese for which investigations were conducted in France in 2012, making it possible to discriminate between epidemiologically related strains and sporadic case strains, while PFGE assigned only a single profile. The six loci selected were sequenced on a large set of strains to determine the sequence of the repeated units and flanking regions, and their stability was evaluated in vivo through the analysis of the strains investigated from humans, food and the farm environment during the outbreak. The six VNTR selected were found to be stable and the discriminatory power of the MLVA method developed was calculated to be 0.954 compared with that for PFGE, which was only 0.625. Twenty-four reference strains were selected from the 401 examined strains in order to represent most of the allele diversity observed for each locus. This reference set can be used to harmonize MLVA results and allow data exchange between laboratories. This original MLVA protocol could be used easily and routinely for monitoring of serovar Dublin isolates and for conducting outbreak investigations.

Most of the bacterial typing methods used to discriminate isolates in medical or food safety microbiology are based on genetic markers used as targets in PCR or hybridization experiments. These DNA typing methods are important tools for studying prevalence and epidemiology, for conducting surveillance, investigations and control of biological hazard sources. In that perspective, it is crucial to insure that the chosen genetic markers have the greatest specificity and sensitivity. The wealth of whole-genome sequences available for many bacterial species offers the opportunity to evaluate the performance of these genetic markers. In the present study, we have developed GTEvaluator, a bioinformatics workflow which ranks genetic markers depending on their sensitivity and specificity towards groups of well-defined genomes. GTEvaluator identifies the most performant genetic markers to target individuals among a population. The individuals (i.e. a group of genomes within a collection) are defined by any kind of particular phenotypic or biological properties inside a related population (i.e. collection of genomes). The performance of the genetic markers is computed by a distance value which takes into account both sensitivity and specificity. In this study we report two examples of GTEvaluator application. In the first example Bacillus phenotypic markers were evaluated for their capacity to distinguish B. cereus from B. thuringiensis. In the second experiment, GTEvaluator measured the performance of genetic markers dedicated to the molecular serotyping of Salmonella enterica. In one in silico experiment it was possible to test 64 markers onto 134 genomes corresponding to 14 different serotypes.

So far only a few publications have explored the development of extraction methods of cyanotoxin extracted from complex matrices. With regard to cyanobacterial microcystins (MCs), the data on the contamination of the flesh of aquatic organisms is hard to compare and very limited due to the lack of validated methods. In recent years, evidence that both free and bound fractions of toxin are found in these tissues has highlighted the need to develop effective methods of quantification. Several techniques do exist, but only the Lemieux oxidation has so far been used to investigate complex tissue matrices. In this study, protocols based on the Lemieux approach were adapted for the quantitative chemical analysis of free MC-LR and MMPB derived from bound toxin in the tissues of juvenile trout gavaged with MC-LR. Afterwards, the NF V03 110 guideline was used to characterize the protocols elaborated and evaluate their effectiveness.

Aquatic ecosystems are facing increasing environmental pressures, leading to an increasing frequency of cyanobacterial Harmful Algal Blooms (cHABs) that have emerged as a worldwide concern due to their growing frequency and their potential toxicity to the fauna that threatens the functioning of ecosystems. Cyanobacterial blooms raise concerns due to the fact that several strains produce potent bioactive or toxic secondary metabolites, such as the microcystins (MCs), which are hepatotoxic to vertebrates. These strains of cyanobacteria may be potentially toxic to fish via gastrointestinal ingestion and also by direct absorption of the toxin MC from the water. The purpose of our study was to investigate toxic effects observed in fish taken from several lakes in the Ile-de-France region, where MCs-producing blooms occur. This study comprises histological studies and the measurement of MC concentrations in various organs. The histological findings are similar to those obtained following laboratory exposure of medaka fish to MCs: hepatic lesions predominate and include cell lysis and cell detachment. MC concentrations in the organs revealed that accumulation was particularly high in the digestive tract and the liver, which are known to be classical targets of MCs. In contrast concentrations were very low in the muscles. Differences in the accumulation of MC variants produced by blooms indicate that in order to more precisely evaluate the toxic potential of a specific bloom it is necessary not only to consider the concentration of toxins, but also the variants produced.

A prolonged outbreak of Salmonella enterica serotype Enteritidis occurred in northern France between December 2014 and April 2015. Epidemiological investigations following the initial notification on 30 December 2014 of five cases of salmonellosis (two confirmed S. Enteritidis) in young children residing in the Somme department revealed that all cases frequented the same food bank A. Further epidemiological, microbiological and food trace-back investigations indicated frozen beefburgers as the source of the outbreak and the suspected lot originating from Poland was recalled on 22 January 2015. On 2 March 2015 a second notification of S. Enteritidis cases in the Somme reinitiated investigations that confirmed a link with food bank A and with consumption of frozen beefburgers from the same Polish producer. In the face of a possible persistent source of contamination, all frozen beefburgers distributed by food bank A and from the same origin were blocked on 3 March 2015. Microbiological analyses confirmed contamination by S. Enteritidis of frozen beefburgers from a second lot remaining in cases’ homes. A second recall was initiated on 6 March 2015 and all frozen beefburgers from the Polish producer remain blocked after analyses identified additional contaminated lots over several months of production.

The aim of this study was to identify and characterise Bacillus cereus from a unique national collection of 564 strains associated with 140 strong-evidence food-borne outbreaks (FBOs) occurring in France during 2007 to 2014. Starchy food and vegetables were the most frequent food vehicles identified; 747 of 911 human cases occurred in institutional catering contexts. Incubation period was significantly shorter for emetic strains compared with diarrhoeal strains A sub-panel of 149 strains strictly associated to 74 FBOs and selected on Coliphage M13-PCR pattern, was studied for detection of the genes encoding cereulide, diarrhoeic toxins (Nhe, Hbl, CytK1 and CytK2) and haemolysin (HlyII), as well as panC phylogenetic classification. This clustered the strains into 12 genetic signatures (GSs) highlighting the virulence potential of each strain. GS1 (nhe genes only) and GS2 (nhe, hbl and cytK2), were the most prevalent GS and may have a large impact on human health as they were present in 28% and 31% of FBOs, respectively. Our study provides a convenient molecular scheme for characterisation of B. cereus strains responsible for FBOs in order to improve the monitoring and investigation of B. cereus-induced FBOs, assess emerging clusters and diversity of strains.

Motivation: MLVA (Multiple Loci VNTR Analysis) is a typing method used today to characterize several major pathogens such as Brucella, Mycobacterium tuberculosis or Salmonella. It takes advantage from the comparison of the size of specific genomic loci constituted by tandem repeat sequences. Unfortunately the raw size estimate is instrument dependent and consequently the results obtained cannot be compared without normalization between laboratories involved in disease monitoring, surveillance and official controls. Results: To overcome this problem we developed a workflow tool, MLVA_normalizer, conceived to normalize MLVA results. This normalization workflow tool is designed to be applied to any bacterial genera and does not depend on the MLVA protocol used. Availability: MLVA_normalizer is available under the GNU general public license (version 2, June 1991). Source code is available at: https://github.com/afelten-Anses/MLVA_normalizer.

Cyanobacteria of the genus Microcystis are known to produce secondary metabolites of large structural diversity by nonribosomal peptide synthetase (NRPS) pathways. For a number of such compounds, halogenated congeners have been reported along with nonhalogenated ones. In the present study, chlorinated cyanopeptolin- and/or aeruginosin-type peptides were detected by mass spectrometry in 17 out of 28 axenic strains of Microcystis. In these strains, a halogenase gene was identified between 2 genes coding for NRPS modules in respective gene clusters, whereas it was consistently absent when the strains produced only nonchlorinated corresponding congeners. Nucleotide sequences were obtained for 12 complete halogenase genes and 14 intermodule regions of gene clusters lacking a halogenase gene or containing only fragments of it. When a halogenase gene was found absent, a specific, identical excision pattern was observed for both synthetase gene clusters in most strains. A phylogenetic analysis including other bacterial halogenases showed that the NRPS-related halogenases of Microcystis form a monophyletic group divided into 2 subgroups, corresponding to either the cyanopeptolin or the aeruginosin peptide synthetases. The distribution of these peptide synthetase gene clusters, among the tested Microcystis strains, was found in relative agreement with their phylogeny reconstructed from 16S-23S rDNA intergenic spacer sequences, whereas the distribution of the associated halogenase genes appears to be sporadic. The presented data suggest that in cyanobacteria these prevalent halogenase genes originated from an ancient horizontal gene transfer followed by duplication in the cyanobacterial lineage. We propose an evolutionary scenario implying repeated gene losses to explain the distribution of halogenase genes in 2 NRPS gene clusters that subsequently defines the seemingly erratic production of halogenated and nonhalogenated aeruginosins and cyanopeptolins among Microcystis strains.

The proliferation of microcystins (MCs)-producing cyanobacteria (MCs) can have detrimental effects on the food chain in aquatic environments. Until recently, few studies had focused on the fate of MCs in exposed organisms, such as primary consumers of cyanobacteria. In this study, we investigate the impact of an MC-producing strain of the cyanobacterium Planktothrix agardhii on the growth and physiology of a Nassula sp. ciliate isolated from a non-toxic cyanobacterial bloom. We show that this Nassula sp. strain was able to consume and grow while feeding exclusively on an MC-producing cyanobacterium over a prolonged period of time (8 months). In short-term exposure experiments (8 days), ciliates consuming an MC-producing cyanobacterial strain displayed slower growth rate and higher levels of antioxidant enzymes than ciliates feeding on two non-MC-producing strains. Three high-performance methods (LC/MS, LC/MS-MS and ELISA) were used to quantify the free and bound MCs in the culture medium and in the cells. We show that ciliate grazing led to a marked decrease in free MCs (methanol extractable) in cells, the MCs were therefore no longer found in the surrounding culture medium. These findings suggest that MCs may have undergone redistribution (free vs bound MCs) or chemical degradation within the ciliates.

Lake Bourget (France) provides drinking water and is a place for professional and recreational fishing. Since the mid 1990s, the lake has been exhibiting blooms of the filamentous cyanobacterium Planktothrix rubescens . This species is able to produce microcystin-LR and RR, toxins that contaminate different fish tissues and, if concentrated in the liver, can induce fish mortality. However, data on fish exposure to these toxins in a natural environment are scare and comparisons of spatial distribution between P. rubescens and exploited fish are needed to determine whether these fish avoid or converge in zones affected by the cyanobacteria. From June to November 2009, diurnal data on P. rubescens and whitefish ( Coregonus lavaretus ) spatial distributions have been monitored by hydroacoustic and BBE probe sampling. For all water samples, intracellular microcystin concentration of P. rubescens was quantified by HPLC/PDA. Furthermore, an arbitrary sample of eight whitefish captured at the moment of highest P. rubescens concentration in the lake were analyzed to investigate both the presence of this cyanobacterium in the gastrointestinal tract by optical microscopy and the possible bioaccumulation of the microcystin in their tissues by liquid chromatography–tandem mass spectrometry. Results show that P. rubescens abundance was at a maximum between the end of July and the beginning of September. During this period, P. rubescens abundance was vertically stratified with a maximum around 14–22 m, depending on the sampling station. The horizontal distributions of P. rubescens were heterogeneous at the scale of the lake. Results indicate that the presence of P. rubescens , for the observed cyanobacterial abundance, does not exert a significant pressure on the distribution pattern of the whitefish. Whitefish were present in the same areas as P. rubescens maxima during daytime, and they do not avoid or seek out the zone of high P. rubescens abundance. Filaments of P. rubescens have been observed in intestinal tracts of whitefish and the presence of microcystin-LR has been detected in their intestine and liver. Consequently, because of a direct contact between these organisms, toxins can be incorporated into whitefish by ingestion of P. rubescens filaments, leading to potential adverse effects on the health of this species.

We report the identification of a sequence from the genome of Oscillatoria sp. strain PCC 6506 coding for a polyketide synthase. Using 50 axenic cyanobacteria, we found this sequence only in the genomes of Oscillatoria strains producing anatoxin-a or homoanatoxin-a, indicating its likely involvement in the biosynthesis of these toxins.

Aeruginosins are bioactive oligopeptides that are produced in high structural diversity by strains of the bloom-forming cyanobacterial genera Microcystis and Planktothrix . A hallmark of aeruginosins is the unusual Choi moiety central to the tetrapeptides, while other positions are occupied by variable moieties in individual congeners. Here we report on three aeruginosin synthetase gene clusters ( aer ) of Microcystis aeruginosa (strains PCC 7806, NIES-98, and NIES-843). The analysis and comparison the aer gene clusters provide the first insight into the molecular basis of biosynthetic and structural plasticity in aeruginosin pathways. Major parts of the aer gene clusters are highly similar in all strains, particularly the genes coding for the first three nonribosomal peptide synthetase (NRPS) modules except for the region coding for the second adenylation domain. However, the gene clusters differ largely in genes coding for tailoring enzymes such as halogenases and sulfotransferases, reflecting structural peculiarities in aeruginosin congeners produced by the individual strains. Significant deviations were further observed in the C-terminal NRPS modules, suggesting two distinct release mechanisms. The architecture of the gene clusters is in agreement with the particular aeruginosin variants that are produced by individual strains, the structures of two of which (aeruginosins 686 A and 686 B) were elucidated. The aer gene clusters of Microcystis and Planktothrix are proposed to originate from a common ancestor and to have evolved to their present-day diversity largely through horizontal gene transfer and recombination events.