European Union Reference Laboratories support the National food, feed and veterinary Reference Laboratories with rolling out whole genome sequencing in Europe.

The Inter European Union Reference Laboratories (EURLs) Working Group on Next Generation Sequencing (NGS) involves eight EURLs for microbiological food and feed hazards and has been working since 2017 to promote the adoption of NGS by the National Reference Laboratories (NRLs) in the European Union. This work illustrates the results of the first 5 years of activity. By working together, the EURLs involved have released guidance documents for assisting NRLs in all the steps of NGS, helping the transition from classical molecular methods towards whole genome sequencing while ensuring harmonization, with the final aim of improving preparedness in the use of NGS to characterize microbial hazards and trace the sources of infection.

The efficiency of Nextera XT tagmentation depends on G and C bases in the binding motif leading to uneven coverage in bacterial species with low and neutral GC-content.

Whole-genome sequencing (WGS) is becoming the new standard for bacterial high-resolution typing and the performance of laboratories is being evaluated in interlaboratory comparisons. The use of the Illumina Nextera XT library preparation kit has been found to be associated with poorer performance due to a GC-content-dependent coverage bias. The bias is especially strong when sequencing low GC-content species. Here, we have made an in-depth analysis of the Nextera XT coverage bias problem using data from a proficiency test of the low GC-content species Campylobacter jejuni. We have compared Nextera XT with Nextera Flex/DNA Prep and examined the consequences on downstream WGS analysis when using different quantities of raw data. We have also analyzed how the coverage bias relates to differential usage of tagmentation cleavage sites. We found that the tagmentation site was characterized by a symmetrical motif with a central AT-rich region surrounded by Gs and Cs. The Gs and Cs appeared to be the main determinant for cleavage efficiency and the genomic regions that were associated with low coverage only contained low-efficiency cleavage sites. This explains why low GC-content genomes and regions are more subjected to coverage bias. We furthermore extended our analysis to other datasets representing other bacterial species. We visualized how the coverage bias was large in low GC-content species such as C. jejuni, C. coli, Staphylococcus aureus, and Listeria monocytogenes, whereas species with neutral GC-content such as Salmonella enterica and Escherichia coli were only affected in certain regions. Species with high GC-content such as Mycobacterium tuberculosis and Pseudomonas aeruginosa were hardly affected at all. The coverage bias associated with Nextera XT was not found when Nextera Flex/DNA Prep had been used.

The Human Adenovirus Type 2 Transcriptome: An Amazing Complexity of Alternatively Spliced mRNAs.

We have used the Nanopore long-read sequencing platform to demonstrate how amazingly complex the human adenovirus type 2 (Ad2) transcriptome is with a flexible splicing machinery producing a range of novel mRNAs both from the early and late transcription units. In total we report more than 900 alternatively spliced mRNAs produced from the Ad2 transcriptome whereof more than 850 are novel mRNAs. A surprising finding was that more than 50% of all E1A transcripts extended upstream of the previously defined transcriptional start site. The novel start sites mapped close to the inverted terminal repeat (ITR) and within the E1A enhancer region. We speculate that novel promoters or enhancer driven transcription, so-called eRNA transcription, is responsible for producing these novel mRNAs. Their existence was verified by a peptide in the Ad2 proteome that was unique for the E1A ITR mRNA. Although we show a high complexity of alternative splicing from most early and late regions, the E3 region was by far the most complex when expressed at late times of infection. More than 400 alternatively spliced mRNAs were observed in this region alone. These mRNAs included extended L4 mRNAs containing E3 and L5 sequences and readthrough mRNAs combining E3 and L5 sequences. Our findings demonstrate that the virus has a remarkable capacity to produce novel exon combinations, which will offer the virus an evolutionary advantage to change the gene expression repertoire and protein production in an evolving environment.IMPORTANCE Work in the adenovirus system led to the groundbreaking discovery of RNA splicing and alternative RNA splicing in 1977. These mechanisms are essential in mammalian evolution by increasing the coding capacity of a genome. Here, we have used a long-read sequencing technology to characterize the complexity of human adenovirus pre-mRNA splicing in detail. It is mindboggling that the viral genome, which only houses around 36,000 bp, not being much larger than a single cellular gene, generates more than 900 alternatively spliced mRNAs. Recently, adenoviruses have been used as the backbone in several promising SARS-CoV-2 vaccines. Further improvement of adenovirus-based vaccines demands that the virus can be tamed into an innocent carrier of foreign genes. This requires a full understanding of the components that govern adenovirus replication and gene expression.

A shotgun metagenomic investigation of the microbiota of udder cleft dermatitis in comparison to healthy skin in dairy cows.

Udder cleft dermatitis (UCD) is a skin condition affecting the fore udder attachment of dairy cows. UCD may be defined as mild (eczematous skin changes) or severe (open wounds, large skin changes). Our aims were to compare the microbiota of mild and severe UCD lesions with the microbiota of healthy skin from the fore udder attachment of control cows, and to investigate whether mastitis-causing pathogens are present in UCD lesions. Samples were obtained from cows in six dairy herds. In total, 36 UCD samples categorized as mild (n = 17) or severe (n = 19) and 13 control samples were sequenced using a shotgun metagenomic approach and the reads were taxonomically classified based on their k-mer content. The Wilcoxon rank sum test was used to compare the abundance of different taxa between different sample types, as well as to compare the bacterial diversity between samples. A high proportion of bacteria was seen in all samples. Control samples had a higher proportion of archaeal reads, whereas most samples had low proportions of fungi, protozoa and viruses. The bacterial microbiota differed between controls and mild and severe UCD samples in both composition and diversity. Subgroups of UCD samples were visible, characterized by increased proportion of one or a few bacterial genera or species, e.g. Corynebacterium, Staphylococcus, Brevibacterium luteolum, Trueperella pyogenes and Fusobacterium necrophorum. Bifidobacterium spp. were more common in controls compared to UCD samples. The bacterial diversity was higher in controls compared to UCD samples. Bacteria commonly associated with mastitis were uncommon. In conclusion, a dysbiosis of the microbiota of mild and severe UCD samples was seen, characterized by decreased diversity and an increased proportion of certain bacteria. There was no evidence of a specific pathogen causing UCD or that UCD lesions are important reservoirs for mastitis-causing bacteria.

The Most Frequently Used Sequencing Technologies and Assembly Methods in Different Time Segments of the Bacterial Surveillance and RefSeq Genome Databases.

Whole genome sequencing has become a powerful tool in modern microbiology. Especially bacterial genomes are sequenced in high numbers. Whole genome sequencing is not only used in research projects, but also in surveillance projects and outbreak investigations. Many whole genome analysis workflows begins with the production of a genome assembly. To accomplish this, a number of different sequencing technologies and assembly methods are available. Here, a summarization is provided over the most frequently used sequence technology and genome assembly approaches reported for the bacterial RefSeq genomes and for the bacterial genomes submitted as belonging to a surveillance project. The data is presented both in total and broken up on a per year basis. Information associated with over 400,000 publically available genomes dated April 2020 and prior were used. The information summarized include (i) the most frequently used sequencing technologies, (ii) the most common combinations of sequencing technologies, (iii) the most reported sequencing depth, and (iv) the most frequently used assembly software solutions. In all, this mini review provides an overview of the currently most common workflows for producing bacterial whole genome sequence assemblies.

Mesenchymal transition and increased therapy resistance of glioblastoma cells is related to astrocyte reactivity.

Grade IV astrocytoma/glioblastoma multiforme (GBM) is essentially incurable, partly due to its heterogenous nature, demonstrated even within the glioma-initiating cell (GIC) population. Increased therapy resistance of GICs is coupled to transition into a mesenchymal (MES) cell state. The GBM MES molecular signature displays a pronounced inflammatory character and its expression vary within and between tumors. Herein, we investigate how MES transition of GBM cells relates to inflammatory responses of normal astroglia. In response to CNS insults astrocytes enter a reactive cell state and participate in directing neuroinflammation and subsequent healing processes. We found that the MES signature show strong resemblance to gene programs induced in reactive astrocytes. Likewise, astrocyte reactivity gene signatures were enriched in therapy-resistant MES-like GIC clones. Variable expression of astrocyte reactivity related genes also largely defined intratumoral GBM cell heterogeneity at the single-cell level and strongly correlated with our previously defined therapy-resistance signature (based on linked molecular and functional characterization of GIC clones). In line with this, therapy-resistant MES-like GIC secreted immunoregulatory and tissue repair related proteins characteristic of astrocyte reactivity. Moreover, sensitive GIC clones could be made reactive through long-term exposure to the proinflammatory cytokine interleukin 1 beta (IL1ß). IL1ß induced a slow MES transition, increased therapy resistance, and a shift in DNA methylation profile towards that of resistant clones, which confirmed a slow reprogramming process. In summary, GICs enter through MES transition a reactive-astrocyte-like cell state, connected to therapy resistance. Thus, from a biological point of view, MES GICs would preferably be called 'reactive GICs'. The ability of GBM cells to mimic astroglial reactivity contextualizes the immunomodulatory and microenvironment reshaping abilities of GBM cells that generate a tumor-promoting milieu. This insight will be important to guide the development of future sensitizing therapies targeting treatment-resistant relapse-driving cell populations as well as enhancing the efficiency of immunotherapies in GBM. © 2019 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Clonal Variation in Drug and Radiation Response among Glioma-Initiating Cells Is Linked to Proneural-Mesenchymal Transition.

Intratumoral heterogeneity is a hallmark of glioblastoma multiforme and thought to negatively affect treatment efficacy. Here, we establish libraries of glioma-initiating cell (GIC) clones from patient samples and find extensive molecular and phenotypic variability among clones, including a range of responses to radiation and drugs. This widespread variability was observed as a continuum of multitherapy resistance phenotypes linked to a proneural-mesenchymal shift in the transcriptome. Multitherapy resistance was associated with a semi-stable cell state that was characterized by an altered DNA methylation pattern at promoter regions of mesenchymal master regulators and enhancers. The gradient of cell states within the GIC compartment constitutes a distinct form of heterogeneity. Our findings may open an avenue toward the development of new therapeutic rationales designed to reverse resistant cell states.

Genetic analysis of a Treponema phagedenis locus encoding antigenic lipoproteins with potential for antigenic variation.

Digital dermatitis (DD) is a painful and debilitating claw disease in cattle. Spirochetes of the genus Treponema are found in high numbers in the lesions and are likely to be involved in the pathogenesis. The occurrence of Treponema phagedenis in DD lesions, especially near the interface of healthy and diseased tissue, suggests that this species contributes to the development and/or progression of the lesions. In this study we characterized a genetic locus in T. phagedenis that contains coding regions for three antigenic proteins, PrrA, VpsA, and VpsB. Comparative analysis of homologous loci from fifteen strains suggests that prrA may be transposed into or out of this locus. Alterations in the copy number of TA repeats within the putative promoter region may regulate VpsA/B expression. The vpsA and prrA genes occur in allelic variants in different T. phagedenis isolates and may provide one explanation for the antigenic variation observed in T. phagedenis DD isolates.

Plasmidome interchange between Clostridium botulinum, Clostridium novyi and Clostridium haemolyticum converts strains of independent lineages into distinctly different pathogens.

Clostridium botulinum (group III), Clostridium novyi and Clostridium haemolyticum are well-known pathogens causing animal botulism, gas gangrene/black disease, and bacillary hemoglobinuria, respectively. A close genetic relationship exists between the species, which has resulted in the collective term C. novyi sensu lato. The pathogenic traits in these species, e.g., the botulinum neurotoxin and the novyi alpha toxin, are mainly linked to a large plasmidome consisting of plasmids and circular prophages. The plasmidome of C. novyi sensu lato has so far been poorly characterized. In this study we explored the genomic relationship of a wide range of strains of C. novyi sensu lato with a special focus on the dynamics of the plasmidome. Twenty-four genomes were sequenced from strains selected to represent as much as possible the genetic diversity in C. novyi sensu lato. Sixty-one plasmids were identified in these genomes and 28 of them were completed. The genomic comparisons revealed four separate lineages, which did not strictly correlate with the species designations. The plasmids were categorized into 13 different plasmid groups on the basis of their similarity and conservation of plasmid replication or partitioning genes. The plasmid groups, lineages and species were to a large extent entwined because plasmids and toxin genes had moved across the lineage boundaries. This dynamic process appears to be primarily driven by phages. We here present a comprehensive characterization of the complex species group C. novyi sensu lato, explaining the intermixed genetic properties. This study also provides examples how the reorganization of the botulinum toxin and the novyi alpha toxin genes within the plasmidome has affected the pathogenesis of the strains.

Small RNA sequence analysis of adenovirus VA RNA-derived miRNAs reveals an unexpected serotype-specific difference in structure and abundance.

Human adenoviruses (HAds) encode for one or two highly abundant virus-associated RNAs, designated VA RNAI and VA RNAII, which fold into stable hairpin structures resembling miRNA precursors. Here we show that the terminal stem of the VA RNAs originating from Ad4, Ad5, Ad11 and Ad37, all undergo Dicer dependent processing into virus-specific miRNAs (so-called mivaRNAs). We further show that the mivaRNA duplex is subjected to a highly asymmetric RISC loading with the 3'-strand from all VA RNAs being the favored strand, except for the Ad37 VA RNAII, where the 5'-mivaRNAII strand was preferentially assembled into RISC. Although the mivaRNA seed sequences are not fully conserved between the HAds a bioinformatics prediction approach suggests that a large fraction of the VA RNAII-, but not the VA RNAI-derived mivaRNAs still are able to target the same cellular genes. Using small RNA deep sequencing we demonstrate that the Dicer processing event in the terminal stem of the VA RNAs is not unique and generates 3'-mivaRNAs with a slight variation of the position of the 5' terminal nucleotide in the RISC loaded guide strand. Also, we show that all analyzed VA RNAs, except Ad37 VA RNAI and Ad5 VA RNAII, utilize an alternative upstream A start site in addition to the classical +1 G start site. Further, the 5'-mivaRNAs with an A start appears to be preferentially incorporated into RISC. Although the majority of mivaRNA research has been done using Ad5 as the model system our analysis demonstrates that the mivaRNAs expressed in Ad11- and Ad37-infected cells are the most abundant mivaRNAs associated with Ago2-containing RISC. Collectively, our results show an unexpected variability in Dicer processing of the VA RNAs and a serotype-specific loading of mivaRNAs into Ago2-based RISC.

Microevolution during an Anthrax outbreak leading to clonal heterogeneity and penicillin resistance.

Anthrax is a bacterial disease primarily affecting grazing animals but it can also cause severe disease in humans. We have used genomic epidemiology to study microevolution of the bacterium in a confined outbreak in cattle which involved emergence of an antibiotic-resistant phenotype. At the time of death, the animals contained a heterogeneous population of Single Nucleotide Variants (SNVs), some being clonal but most being subclonal. We found that independent isolates from the same carcass had similar levels of SNV differences as isolates from different animals. Furthermore the relative levels of subclonal populations were different in different locations in the same carcass. The heterogeneity appeared to be derived in part from heterogeneity in the infectious dose. The resistance phenotype was linked to clonal mutations in an anti-sigma factor gene and in one case was preceded by an acquisition of a hypermutator phenotype. In another animal, small subclonal populations were observed with counteracting mutations that had turned off the resistance genes. In summary, this study shows the importance of accounting for both acquired and inherited heterogeneity when doing high-resolution infection tracing and when estimating the risks associated with penicillin treatment.

In silico and in vitro evaluation of PCR-based assays for the detection of Bacillus anthracis chromosomal signature sequences.

Bacillus anthracis, the causative agent of anthrax, is a zoonotic pathogen that is relatively common throughout the world and may cause life threatening diseases in animals and humans. There are many PCR-based assays in use for the detection of B. anthracis. While most of the developed assays rely on unique markers present on virulence plasmids pXO1 and pXO2, relatively few assays incorporate chromosomal DNA markers due to the close relatedness of B. anthracis to the B. cereus group strains. For the detection of chromosomal DNA, different genes have been used, such as BA813, rpoB, gyrA, plcR, S-layer, and prophage-lambda. Following a review of the literature, an in silico analysis of all signature sequences reported for identification of B. anthracis was conducted. Published primer and probe sequences were compared for specificity against 134 available Bacillus spp. genomes. Although many of the chromosomal targets evaluated are claimed to be specific to B. anthracis, cross-reactions with closely related B. cereus and B. thuringiensis strains were often observed. Of the 35 investigated PCR assays, only 4 were 100% specific for the B. anthracis chromosome. An interlaboratory ring trial among five European laboratories was then performed to evaluate six assays, including the WHO recommended procedures, using a collection of 90 Bacillus strains. Three assays performed adequately, yielding no false positive or negative results. All three assays target chromosomal markers located within the lambdaBa03 prophage region (PL3, BA5345, and BA5357). Detection limit was further assessed for one of these highly specific assays.

Management of animal botulism outbreaks: from clinical suspicion to practical countermeasures to prevent or minimize outbreaks.

Botulism is a severe neuroparalytic disease that affects humans, all warm-blooded animals, and some fishes. The disease is caused by exposure to toxins produced by Clostridium botulinum and other botulinum toxin-producing clostridia. Botulism in animals represents a severe environmental and economic concern because of its high mortality rate. Moreover, meat or other products from affected animals entering the food chain may result in a public health problem. To this end, early diagnosis is crucial to define and apply appropriate veterinary public health measures. Clinical diagnosis is based on clinical findings eliminating other causes of neuromuscular disorders and on the absence of internal lesions observed during postmortem examination. Since clinical signs alone are often insufficient to make a definitive diagnosis, laboratory confirmation is required. Botulinum antitoxin administration and supportive therapies are used to treat sick animals. Once the diagnosis has been made, euthanasia is frequently advisable. Vaccine administration is subject to health authorities' permission, and it is restricted to a small number of animal species. Several measures can be adopted to prevent or minimize outbreaks. In this article we outline all phases of management of animal botulism outbreaks occurring in wet wild birds, poultry, cattle, horses, and fur farm animals.

Multiplex real-time PCR for detecting and typing Clostridium botulinum group III organisms and their mosaic variants.

Botulism is a neuroparalytic disease that can occur in all warm-blooded animals, birds, and fishes. The disease in animals is mainly caused by toxins produced by Clostridium botulinum strains belonging to group III, although outbreaks due to toxins produced by group I and II organisms have been recognized. Group III strains are capable of producing botulinum toxins of type C, D, and C/D and D/C mosaic variants. Definitive diagnosis of animal botulism is made by combining clinical findings with laboratory investigations. Detection of toxins in clinical specimens and feed is the gold standard for laboratory diagnosis. Since toxins may be degraded by organisms contained in the gastrointestinal tract or may be present at levels below the detection limit, the recovery of C. botulinum from sick animal specimens is consistent for laboratory confirmation. In this article we report the development and in-house validation of a new multiplex real-time PCR for detecting and typing the neurotoxin genes found in C. botulinum group III organisms. Validation procedures have been carried out according to ISO 16140, using strains and samples recovered from cases of animal botulism in Italy and France.

Genome-wide relatedness of Treponema pedis, from gingiva and necrotic skin lesions of pigs, with the human oral pathogen Treponema denticola.

Treponema pedis and T. denticola are two genetically related species with different origins of isolation. Treponema denticola is part of the human oral microbiota and is associated with periodontitis while T. pedis has been isolated from skin lesions in animals, e.g., digital dermatitis in cattle and necrotic ulcers in pigs. Although multiple Treponema phylotypes may exist in ulcerative lesions in pigs, T. pedis appears to be a predominant spirochete in these lesions. Treponema pedis can also be present in pig gingiva. In this study, we determined the complete genome sequence of T. pedis strain T A4, isolated from a porcine necrotic ear lesion, and compared its genome with that of T. denticola. Most genes in T. pedis were homologous to those in T. denticola and the two species were similar in general genomic features such as size, G+C content, and number of genes. In addition, many homologues of specific virulence-related genes in T. denticola were found in T. pedis. Comparing a selected pair of strains will usually not give a complete picture of the relatedness between two species. We therefore complemented the analysis with draft genomes from six T. pedis isolates, originating from gingiva and necrotic ulcers in pigs, and from twelve T. denticola strains. Each strain carried a considerable amount of accessory genetic material, of which a large part was strain specific. There was also extensive sequence variability in putative virulence-related genes between strains belonging to the same species. Signs of lateral gene-transfer events from bacteria known to colonize oral environments were found. This suggests that the oral cavity is an important habitat for T. pedis. In summary, we found extensive genomic similarities between T. pedis and T. denticola but also large variability within each species.

The adenovirus VA RNA-derived miRNAs are not essential for lytic virus growth in tissue culture cells.

At late times during a lytic infection human adenovirus type 5 produces ∼10(8) copies per cell of virus-associated RNA I (VA RNAI). This short highly structured RNA polymerase III transcript has previously been shown to be essential for lytic virus growth. A fraction of VA RNAI is processed by Dicer into small RNAs, so-called mivaRNAIs, which are efficiently incorporated into the RNA-induced silencing complex. Here, we constructed recombinant adenoviruses with mutations in the seed sequence of both the 5'- and the 3'-strand of the mivaRNAI duplex. The results showed that late viral protein synthesis, as well as new virus progeny formation, was essentially unaffected by the seed sequence mutations under lytic replicative conditions in HeLa or HEK293 cells. Collectively, our results suggest that either strand of the mivaRNAI duplex does not have target mRNA interactions that are critical for the establishment of virus growth under lytic conditions. Further, by depletion of protein kinase R (PKR) in HEK293 cells, we show that the suppressive effect of VA RNAI on the interferon-induced PKR pathway is most critical for late gene expression.

The genetic integrity of bacterial species: the core genome and the accessory genome, two different stories.

Strains within a bacterial species typically have a set of conserved core genes and a variable set of accessory genes. The accessory genes often appear to move laterally between strains, thereby forming new trait combinations. Sometimes, genetic material also moves laterally between species, thereby resulting in diffuse borders between them. The growing number of genome sequences offers new possibilities to study these processes. Ten species for which abundant genomic data exists were here selected for analysis of the species border integrity. The average core genome similarities and relative core genome sizes (RCGSs) were determined for strain pairs within the species and for strain pairs crossing the species border. The variability within the species as well as the border integrity varies for different bacterial species. Some have very distinct borders while others are more or less indefinable. From the growing amount of genomic data, it becomes even clearer that the concept of bacterial species is, in many cases, far from absolute.

Gegenees: fragmented alignment of multiple genomes for determining phylogenomic distances and genetic signatures unique for specified target groups.

The rapid development of Next Generation Sequencing technologies leads to the accumulation of huge amounts of sequencing data. The scientific community faces an enormous challenge in how to deal with this explosion. Here we present a software tool, 'Gegenees', that uses a fragmented alignment approach to facilitate the comparative analysis of hundreds of microbial genomes. The genomes are fragmented and compared, all against all, by a multithreaded BLAST control engine. Ready-made alignments can be complemented with new genomes without recalculating the existing data points. Gegenees gives a phylogenomic overview of the genomes and the alignment can then be mined for genomic regions with conservation patterns matching a defined target group and absent from a background group. The genomic regions are given biomarker scores forming a uniqueness signature that can be viewed and explored, graphically and in tabular form. A primer/probe alignment tool is also included for specificity verification of currently used or new primers. We exemplify the use of Gegenees on the Bacillus cereus group, on Foot and Mouth Disease Viruses, and on strains from the 2011 Escherichia coli O104:H4 outbreak. Gegenees contributes towards an increased capacity of fast and efficient data mining as more and more genomes become sequenced.

Complete genome sequence of Brachyspira intermedia reveals unique genomic features in Brachyspira species and phage-mediated horizontal gene transfer.

BACKGROUND: Brachyspira spp. colonize the intestines of some mammalian and avian species and show different degrees of enteropathogenicity. Brachyspira intermedia can cause production losses in chickens and strain PWS/AT now becomes the fourth genome to be completed in the genus Brachyspira. RESULTS: 15 classes of unique and shared genes were analyzed in B. intermedia, B. murdochii, B. hyodysenteriae and B. pilosicoli. The largest number of unique genes was found in B. intermedia and B. murdochii. This indicates the presence of larger pan-genomes. In general, hypothetical protein annotations are overrepresented among the unique genes. A 3.2 kb plasmid was found in B. intermedia strain PWS/AT. The plasmid was also present in the B. murdochii strain but not in nine other Brachyspira isolates. Within the Brachyspira genomes, genes had been translocated and also frequently switched between leading and lagging strands, a process that can be followed by different AT-skews in the third positions of synonymous codons. We also found evidence that bacteriophages were being remodeled and genes incorporated into them. CONCLUSIONS: The accessory gene pool shapes species-specific traits. It is also influenced by reductive genome evolution and horizontal gene transfer. Gene-transfer events can cross both species and genus boundaries and bacteriophages appear to play an important role in this process. A mechanism for horizontal gene transfer appears to be gene translocations leading to remodeling of bacteriophages in combination with broad tropism.

Clostridium botulinum group III: a group with dual identity shaped by plasmids, phages and mobile elements.

BACKGROUND: Clostridium botulinum strains can be divided into four physiological groups that are sufficiently diverged to be considered as separate species. Here we present the first complete genome of a C. botulinum strain from physiological group III, causing animal botulism. We also compare the sequence to three new draft genomes from the same physiological group. RESULTS: The 2.77 Mb chromosome was highly conserved between the isolates and also closely related to that of C. novyi. However, the sequence was very different from the human C. botulinum group genomes. Replication-directed translocations were rare and conservation of synteny was high. The largest difference between C. botulinum group III isolates occurred within their surprisingly large plasmidomes and in the pattern of mobile elements insertions. Five plasmids, constituting 13.5% of the total genetic material, were present in the completed genome. Interestingly, the set of plasmids differed compared to other isolates. The largest plasmid, the botulinum-neurotoxin carrying prophage, was conserved at a level similar to that of the chromosome while the medium-sized plasmids seemed to be undergoing faster genetic drift. These plasmids also contained more mobile elements than other replicons. Several toxins and resistance genes were identified, many of which were located on the plasmids. CONCLUSIONS: The completion of the genome of C. botulinum group III has revealed it to be a genome with dual identity. It belongs to the pathogenic species C. botulinum, but as a genotypic species it should also include C. novyi and C. haemolyticum. The genotypic species share a conserved chromosomal core that can be transformed into various pathogenic variants by modulation of the highly plastic plasmidome.