Comparison of Illumina and Oxford Nanopore Technology for genome analysis of Francisella tularensis, Bacillus anthracis, and Brucella suis.

BACKGROUND: Bacterial epidemiology needs to understand the spread and dissemination of strains in a One Health context. This is important for highly pathogenic bacteria such as Bacillus anthracis, Brucella species, and Francisella tularensis. Whole genome sequencing (WGS) has paved the way for genetic marker detection and high-resolution genotyping. While such tasks are established for Illumina short-read sequencing, Oxford Nanopore Technology (ONT) long-read sequencing has yet to be evaluated for such highly pathogenic bacteria with little genomic variations between strains. In this study, three independent sequencing runs were performed using Illumina, ONT flow cell version 9.4.1, and 10.4 for six strains of each of Ba. anthracis, Br. suis and F. tularensis. Data from ONT sequencing alone, Illumina sequencing alone and two hybrid assembly approaches were compared. RESULTS: As previously shown, ONT produces ultra-long reads, while Illumina produces short reads with higher sequencing accuracy. Flow cell version 10.4 improved sequencing accuracy over version 9.4.1. The correct (sub-)species were inferred from all tested technologies, individually. Moreover, the sets of genetic markers for virulence, were almost identical for the respective species. The long reads of ONT allowed to assemble not only chromosomes of all species to near closure, but also virulence plasmids of Ba. anthracis. Assemblies based on nanopore data alone, Illumina data alone, and both hybrid assemblies correctly detected canonical (sub-)clades for Ba. anthracis and F. tularensis as well as multilocus sequence types for Br. suis. For F. tularensis, high-resolution genotyping using core-genome MLST (cgMLST) and core-genome Single-Nucleotide-Polymorphism (cgSNP) typing produced highly comparable results between data from Illumina and both ONT flow cell versions. For Ba. anthracis, only data from flow cell version 10.4 produced similar results to Illumina for both high-resolution typing methods. However, for Br. suis, high-resolution genotyping yielded larger differences comparing Illumina data to data from both ONT flow cell versions. CONCLUSIONS: In summary, combining data from ONT and Illumina for high-resolution genotyping might be feasible for F. tularensis and Ba. anthracis, but not yet for Br. suis. The ongoing improvement of nanopore technology and subsequent data analysis may facilitate high-resolution genotyping for all bacteria with highly stable genomes in future.

Whole genome sequencing of Salmonella enterica serovars isolated from humans, animals, and the environment in Lagos, Nigeria.

BACKGROUND: Salmonella infections remain an important public health issue worldwide. Some serovars of non-typhoidal Salmonella (NTS) have been associated with bloodstream infections and gastroenteritis, especially in children in Sub-Saharan Africa with circulating S. enterica serovars with drug resistance and virulence genes. This study identified and verified the clonal relationship of Nigerian NTS strains isolated from humans, animals, and the environment. METHODS: In total, 2,522 samples were collected from patients, animals (cattle and poultry), and environmental sources between December 2017 and May 2019. The samples were subjected to a standard microbiological investigation. All the isolates were identified using Microbact 24E, and MALDI-TOF MS. The isolates were serotyped using the Kauffmann-White scheme. Antibiotic susceptibility testing was conducted using the disc diffusion method and the Vitek 2 compact system. Virulence and antimicrobial resistance genes, sequence type, and cluster analysis were investigated using WGS data. RESULTS: Forty-eight (48) NTS isolates (1.9%) were obtained. The prevalence of NTS from clinical sources was 0.9%, while 4% was recorded for animal sources. The serovars identified were S. Cotham (n = 17), S. Give (n = 16), S. Mokola (n = 6), S. Abony (n = 4), S. Typhimurium (n = 4), and S. Senftenberg (n = 1). All 48 Salmonella isolates carried intrinsic and acquired resistant genes such as aac.6…Iaa, mdf(A), qnrB, qnrB19 genes and golT, golS, pcoA, and silP, mediated by plasmid Col440I_1, incFIB.B and incFII. Between 100 and 118 virulence gene markers distributed across several Salmonella pathogenicity islands (SPIs), clusters, prophages, and plasmid operons were found in each isolate. WGS revealed that strains of each Salmonella serovar could be assigned to a single 7-gene MLST cluster, and strains within the clusters were identical strains and closely related as defined by the 0 and 10 cgSNPs and likely shared a common ancestor. The dominant sequence types were S. Give ST516 and S. Cotham ST617. CONCLUSION: We found identical Salmonella sequence types in human, animal, and environmental samples in the same locality, which demonstrates the great potential of the applied tools to trace back outbreak strains. Strategies to control and prevent the spread of NTS in the context of one's health are essential to prevent possible outbreaks.

[Assessment of available and currently applied typing methods including genome-based methods for zoonotic pathogens with a focus on Salmonella enterica]. / Bestandsaufnahme der verfügbaren und aktuell eingesetzten Typisierungsmethoden einschließlich genombasierter Verfahren von Zoonoseerregern am Beispiel von Salmonella enterica.

BACKGROUND: In recent years, whole genome sequencing (WGS) in combination with bioinformatic analyses has become state of the art in evaluating the pathogenicity/resistance potential and relatedness of bacteria. WGS analysis thus represents a central tool in the investigation of the resistance and virulence potential of pathogens, as well as their dissemination via outbreak clusters and transmission chains within the framework of molecular epidemiology. In order to gain an overview of the available genotypic and phenotypic methods used for pathogen typing of Salmonella and Shiga toxin-producing and enterohemorrhagic Escherichia coli (STEC/EHEC) in Germany at state and federal level, along with the availability of WGS-based typing and corresponding analytical methods, a survey of laboratories was conducted. METHODS: An electronic survey of laboratories working for public health protection and consumer health protection was conducted from February to June 2020. RESULTS AND CONCLUSION: The results of the survey showed that many of the participating laboratories provide a wide range of phenotypic and molecular methods. Molecular typing is most commonly used for species identification of Salmonella. In many cases, WGS-based methods have already been established at federal and state institutions or are in the process of being established. The Illumina sequencing technology is the most widely used technology. The survey confirms the importance of molecular biology and whole genome typing technologies for laboratories in the diagnosis of bacterial zoonotic pathogens.

The fungivorous amoeba Protostelium aurantium targets redox homeostasis and cell wall integrity during intracellular killing of Candida parapsilosis.

Predatory interactions among microbes are major evolutionary driving forces for biodiversity. The fungivorous amoeba Protostelium aurantium has a wide fungal food spectrum including foremost pathogenic members of the genus Candida. Here we show that upon phagocytic ingestion by the amoeba, Candida parapsilosis is confronted with an oxidative burst and undergoes lysis within minutes of processing in acidified phagolysosomes. On the fungal side, a functional genomic approach identified copper and redox homeostasis as primary targets of amoeba predation, with the highly expressed copper exporter gene CRP1 and the peroxiredoxin gene PRX1 contributing to survival when encountered with P. aurantium. The fungicidal activity was largely retained in intracellular vesicles of the amoebae. Following their isolation, the content of these vesicles induced immediate killing and lysis of C. parapsilosis in vitro. Proteomic analysis identified 56 vesicular proteins from P. aurantium. Although completely unknown proteins were dominant, many of them could be categorised as hydrolytic enzymes targeting the fungal cell wall, indicating that fungal cell wall structures are under selection pressure by predatory phagocytes in natural environments. TAKE AWAY: The amoeba Protostelium aurantium feeds on fungi, such as Candida parapsilosis. Ingested yeast cells are exposed to reactive oxygen species. A copper exporter and a peroxiredoxin contribute to fungal defence. Yeast cells undergo intracellular lysis. Lysis occurs via a cocktail of hydrolytic enzymes from intracellular vesicles.

Sequence Analysis of Novel Staphylococcus aureus Lineages from Wild and Captive Macaques.

Staphylococcus aureus is a widespread and common opportunistic bacterium that can colonise or infect humans as well as a wide range of animals. There are a few studies of both methicillin-susceptible S. aureus (MSSA) and methicillin-resistant S. aureus (MRSA) isolated from monkeys, apes, and lemurs, indicating a presence of a number of poorly or unknown lineages of the pathogen. In order to obtain insight into staphylococcal diversity, we sequenced strains from wild and captive individuals of three macaque species (Macaca mulatta, M. assamensis, and M. sylvanus) using Nanopore and Illumina technologies. These strains were previously identified by microarray as poorly or unknown strains. Isolates of novel lineages ST4168, ST7687, ST7688, ST7689, ST7690, ST7691, ST7692, ST7693, ST7694, ST7695, ST7745, ST7746, ST7747, ST7748, ST7749, ST7750, ST7751, ST7752, ST7753, and ST7754 were sequenced and characterised for the first time. In addition, isolates belonging to ST2990, a lineage also observed in humans, and ST3268, a MRSA strain already known from macaques, were also included into the study. Mobile genetic elements, genomic islands, and carriage of prophages were analysed. There was no evidence for novel host-specific virulence factors. However, a conspicuously high rate of carriage of a pathogenicity island harbouring edinB and etD2/etE as well as a higher number of repeat units within the gene sasG (encoding an adhesion factor) than in human isolates were observed. None of the strains harboured the genes encoding Panton-Valentine leukocidin. In conclusion, wildlife including macaques may harbour an unappreciated diversity of S. aureus lineages that may be of clinical relevance for humans, livestock, or for wildlife conservation, given the declining state of many wildlife populations.

A Whole-Genome-Based Gene-by-Gene Typing System for Standardized High-Resolution Strain Typing of Bacillus anthracis.

Whole-genome sequencing (WGS) has been established for bacterial subtyping and is regularly used to study pathogen transmission, to investigate outbreaks, and to perform routine surveillance. Core-genome multilocus sequence typing (cgMLST) is a bacterial subtyping method that uses WGS data to provide a high-resolution strain characterization. This study aimed at developing a novel cgMLST scheme for Bacillus anthracis, a notorious pathogen that causes anthrax in livestock and humans worldwide. The scheme comprises 3,803 genes that were conserved in 57 B. anthracis genomes spanning the whole phylogeny. The scheme has been evaluated and applied to 584 genomes from 50 countries. On average, 99.5% of the cgMLST targets were detected. The cgMLST results confirmed the classical canonical single-nucleotide-polymorphism (SNP) grouping of B. anthracis into major clades and subclades. Genetic distances calculated based on cgMLST were comparable to distances from whole-genome-based SNP analysis with similar phylogenetic topology and comparable discriminatory power. Additionally, the application of the cgMLST scheme to anthrax outbreaks from Germany and Italy led to a definition of a cutoff threshold of five allele differences to trace epidemiologically linked strains for cluster typing and transmission analysis. Finally, the association of two clusters of B. anthracis with human cases of injectional anthrax in four European countries was confirmed using cgMLST. In summary, this study presents a novel cgMLST scheme that provides high-resolution strain genotyping for B. anthracis. This scheme can be used in parallel with SNP typing methods to facilitate rapid and harmonized interlaboratory comparisons, essential for global surveillance and outbreak analysis. The scheme is publicly available for application by users, including those with little bioinformatics knowledge.

Comparative Genomics of Serial Candida glabrata Isolates and the Rapid Acquisition of Echinocandin Resistance during Therapy.

The opportunistic pathogen Candida glabrata shows a concerning increase in drug resistance. Here, we present the analysis of two serial bloodstream isolates, obtained 12 days apart. Both isolates show pan-azole resistance and echinocandin resistance was acquired during the sampling interval. Genome sequencing identified nine nonsynonymous SNVs between the strains, including a S663P substitution in FKS2 and previously undescribed SNVs in MDE1 and FPR1, offering insight into how C. glabrata acquires drug resistance and adapts to a human host.

Correction: Gene Expansion Shapes Genome Architecture in the Human Pathogen Lichtheimia corymbifera: An Evolutionary Genomics Analysis in the Ancient Terrestrial Mucorales (Mucoromycotina).

[This corrects the article DOI: 10.1371/journal.pgen.1004496.].

Smelling the difference: Transcriptome, proteome and volatilome changes after mating.

Mushrooms, such as Schizophyllum commune, have a specific odor. Whether this is linked to mating, prerequisite for mushroom formation, or also found in monokaryotic, unmated strains, was investigated with a comprehensive study on the transcriptome and proteome of this model organism. Mating interactions were investigated using a complete, cytosolic proteome map for unmated, monokaryotic, as well as for mated, dikaryotic mycelia. The regulations of the proteome were compared to transcriptional changes upon mating and to changes in smell by volatilome studies. We could show a good overlap between proteome and transcriptome data, but extensive posttranslational regulation was identified for more than 80% of transcripts. This suggests down-stream regulation upon interaction of mating partners and formation of the dikaryon that is competent to form fruiting bodies. The volatilome was shown to respond to mating by a broader spectrum of volatiles and increased emission of the mushroom smell molecules 3-octanone and 1-octen-3-ol, as well as ethanol and ß-bisabolol in the dikaryon. Putatively involved biosynthetic proteins like alcohol dehydrogenases, Ppo-like oxygenases, or sesquiterpene synthases showed correlating transcriptional regulation depending on either mono- or dikaryotic stages.

LASSIM-A network inference toolbox for genome-wide mechanistic modeling.

Recent technological advancements have made time-resolved, quantitative, multi-omics data available for many model systems, which could be integrated for systems pharmacokinetic use. Here, we present large-scale simulation modeling (LASSIM), which is a novel mathematical tool for performing large-scale inference using mechanistically defined ordinary differential equations (ODE) for gene regulatory networks (GRNs). LASSIM integrates structural knowledge about regulatory interactions and non-linear equations with multiple steady state and dynamic response expression datasets. The rationale behind LASSIM is that biological GRNs can be simplified using a limited subset of core genes that are assumed to regulate all other gene transcription events in the network. The LASSIM method is implemented as a general-purpose toolbox using the PyGMO Python package to make the most of multicore computers and high performance clusters, and is available at https://gitlab.com/Gustafsson-lab/lassim. As a method, LASSIM works in two steps, where it first infers a non-linear ODE system of the pre-specified core gene expression. Second, LASSIM in parallel optimizes the parameters that model the regulation of peripheral genes by core system genes. We showed the usefulness of this method by applying LASSIM to infer a large-scale non-linear model of naïve Th2 cell differentiation, made possible by integrating Th2 specific bindings, time-series together with six public and six novel siRNA-mediated knock-down experiments. ChIP-seq showed significant overlap for all tested transcription factors. Next, we performed novel time-series measurements of total T-cells during differentiation towards Th2 and verified that our LASSIM model could monitor those data significantly better than comparable models that used the same Th2 bindings. In summary, the LASSIM toolbox opens the door to a new type of model-based data analysis that combines the strengths of reliable mechanistic models with truly systems-level data. We demonstrate the power of this approach by inferring a mechanistically motivated, genome-wide model of the Th2 transcription regulatory system, which plays an important role in several immune related diseases.

Candida albicans infection leads to barrier breakdown and a MAPK/NF-κB mediated stress response in the intestinal epithelial cell line C2BBe1.

Intestinal epithelial cells (IEC) form a tight barrier to the gut lumen. Paracellular permeability of the intestinal barrier is regulated by tight junction proteins and can be modulated by microorganisms and other stimuli. The polymorphic fungus Candida albicans, a frequent commensal of the human mucosa, has the capacity of traversing this barrier and establishing systemic disease within the host. Infection of polarized C2BBe1 IEC with wild-type C. albicans led to a transient increase of transepithelial electric resistance (TEER) before subsequent barrier disruption, accompanied by a strong decline of junctional protein levels and substantial, but considerably delayed cytotoxicity. Time-resolved microarray-based transcriptome analysis of C. albicans challenged IEC revealed a prominent role of NF-κB and MAPK signalling pathways in the response to infection. Hence, we inferred a gene regulatory network based on differentially expressed NF-κB and MAPK pathway components and their predicted transcriptional targets. The network model predicted activation of GDF15 by NF-κB was experimentally validated. Furthermore, inhibition of NF-κB activation in C. albicans infected C2BBe1 cells led to enhanced cytotoxicity in the epithelial cells. Taken together our study identifies NF-κB activation as an important protective signalling pathway in the response of epithelial cells to C. albicans.

Defining the transcriptomic landscape of Candida glabrata by RNA-Seq.

Candida glabrata is the second most common pathogenic Candida species and has emerged as a leading cause of nosocomial fungal infections. Its reduced susceptibility to antifungal drugs and its close relationship to Saccharomyces cerevisiae make it an interesting research focus. Although its genome sequence was published in 2004, little is known about its transcriptional dynamics. Here, we provide a detailed RNA-Seq-based analysis of the transcriptomic landscape of C. glabrata in nutrient-rich media, as well as under nitrosative stress and during pH shift. Using RNA-Seq data together with state-of-the-art gene prediction tools, we refined the annotation of the C. glabrata genome and predicted 49 novel protein-coding genes. Of these novel genes, 14 have homologs in S. cerevisiae and six are shared with other Candida species. We experimentally validated four novel protein-coding genes of which two are differentially regulated during pH shift and interaction with human neutrophils, indicating a potential role in host-pathogen interaction. Furthermore, we identified 58 novel non-protein-coding genes, 38 new introns and condition-specific alternative splicing. Finally, our data suggest different patterns of adaptation to pH shift and nitrosative stress in C. glabrata, Candida albicans and S. cerevisiae and thus further underline a distinct evolution of virulence in yeast.

Gene expansion shapes genome architecture in the human pathogen Lichtheimia corymbifera: an evolutionary genomics analysis in the ancient terrestrial mucorales (Mucoromycotina).

Lichtheimia species are the second most important cause of mucormycosis in Europe. To provide broader insights into the molecular basis of the pathogenicity-associated traits of the basal Mucorales, we report the full genome sequence of L. corymbifera and compared it to the genome of Rhizopus oryzae, the most common cause of mucormycosis worldwide. The genome assembly encompasses 33.6 MB and 12,379 protein-coding genes. This study reveals four major differences of the L. corymbifera genome to R. oryzae: (i) the presence of an highly elevated number of gene duplications which are unlike R. oryzae not due to whole genome duplication (WGD), (ii) despite the relatively high incidence of introns, alternative splicing (AS) is not frequently observed for the generation of paralogs and in response to stress, (iii) the content of repetitive elements is strikingly low (<5%), (iv) L. corymbifera is typically haploid. Novel virulence factors were identified which may be involved in the regulation of the adaptation to iron-limitation, e.g. LCor01340.1 encoding a putative siderophore transporter and LCor00410.1 involved in the siderophore metabolism. Genes encoding the transcription factors LCor08192.1 and LCor01236.1, which are similar to GATA type regulators and to calcineurin regulated CRZ1, respectively, indicating an involvement of the calcineurin pathway in the adaption to iron limitation. Genes encoding MADS-box transcription factors are elevated up to 11 copies compared to the 1-4 copies usually found in other fungi. More findings are: (i) lower content of tRNAs, but unique codons in L. corymbifera, (ii) Over 25% of the proteins are apparently specific for L. corymbifera. (iii) L. corymbifera contains only 2/3 of the proteases (known to be essential virulence factors) in comparison to R. oryzae. On the other hand, the number of secreted proteases, however, is roughly twice as high as in R. oryzae.

Microevolution of Candida albicans in macrophages restores filamentation in a nonfilamentous mutant.

Following antifungal treatment, Candida albicans, and other human pathogenic fungi can undergo microevolution, which leads to the emergence of drug resistance. However, the capacity for microevolutionary adaptation of fungi goes beyond the development of resistance against antifungals. Here we used an experimental microevolution approach to show that one of the central pathogenicity mechanisms of C. albicans, the yeast-to-hyphae transition, can be subject to experimental evolution. The C. albicans cph1Δ/efg1Δ mutant is nonfilamentous, as central signaling pathways linking environmental cues to hyphal formation are disrupted. We subjected this mutant to constant selection pressure in the hostile environment of the macrophage phagosome. In a comparatively short time-frame, the mutant evolved the ability to escape macrophages by filamentation. In addition, the evolved mutant exhibited hyper-virulence in a murine infection model and an altered cell wall composition compared to the cph1Δ/efg1Δ strain. Moreover, the transcriptional regulation of hyphae-associated, and other pathogenicity-related genes became re-responsive to environmental cues in the evolved strain. We went on to identify the causative missense mutation via whole genome- and transcriptome-sequencing: a single nucleotide exchange took place within SSN3 that encodes a component of the Cdk8 module of the Mediator complex, which links transcription factors with the general transcription machinery. This mutation was responsible for the reconnection of the hyphal growth program with environmental signals in the evolved strain and was sufficient to bypass Efg1/Cph1-dependent filamentation. These data demonstrate that even central transcriptional networks can be remodeled very quickly under appropriate selection pressure.

FungiFun2: a comprehensive online resource for systematic analysis of gene lists from fungal species.

SUMMARY: Systematically extracting biological meaning from omics data is a major challenge in systems biology. Enrichment analysis is often used to identify characteristic patterns in candidate lists. FungiFun is a user-friendly Web tool for functional enrichment analysis of fungal genes and proteins. The novel tool FungiFun2 uses a completely revised data management system and thus allows enrichment analysis for 298 currently available fungal strains published in standard databases. FungiFun2 offers a modern Web interface and creates interactive tables, charts and figures, which users can directly manipulate to their needs. AVAILABILITY AND IMPLEMENTATION: FungiFun2, examples and tutorials are publicly available at https://elbe.hki-jena.de/fungifun/. CONTACT: steffen.priebe@hki-jena.de or joerg.linde@hki-jena.de.

The novel globin protein fungoglobin is involved in low oxygen adaptation of Aspergillus fumigatus.

The human pathogenic fungus Aspergillus fumigatus normally lives as a soil saprophyte. Its environment includes poorly oxygenated substrates that also occur during tissue invasive growth of the fungus in the human host. Up to now, few cellular factors have been identified that allow the fungus to efficiently adapt its energy metabolism to hypoxia. Here, we cultivated A. fumigatus in an O2 -controlled fermenter and analysed its responses to O2 limitation on a minute timescale. Transcriptome sequencing revealed several genes displaying a rapid and highly dynamic regulation. One of these genes was analysed in detail and found to encode fungoglobin, a previously uncharacterized member of the sensor globin protein family widely conserved in filamentous fungi. Besides low O2 , iron limitation also induced transcription, but regulation was not entirely dependent on the two major transcription factors involved in adaptation to iron starvation and hypoxia, HapX and SrbA respectively. The protein was identified as a functional haemoglobin, as binding of this cofactor was detected for the recombinant protein. Gene deletion in A. fumigatus confirmed that haem-binding fungoglobins are important for growth in microaerobic environments with O2 levels far lower than in hypoxic human tissue.

Epidemiological Analysis on the Occurrence of Salmonella enterica Subspecies enterica Serovar Dublin in the German Federal State Schleswig-Holstein Using Whole-Genome Sequencing.

The cattle-adapted serovar Salmonella Dublin (S. Dublin) causes enteritis and systemic diseases in animals. In the German federal state Schleswig-Holstein, S. Dublin is the most important serovar in cattle indicating an endemic character of the infection. To gain information on dissemination and routes of infection, whole-genome sequencing (WGS) was used to explore the genetic traits of 78 S. Dublin strains collected over a period of six years. The phylogeny was analysed using core-genome single nucleotide polymorphisms (cgSNPs). Genomic clusters at 100, 15 and 1 cgSNPs were selected for molecular analysis. Important specific virulence determinants were detected in all strains but multidrug resistance in S. Dublin organisms was not found. Using 15 cgSNPs epidemiological links between herds were identified, clusters at 1 cgSNPs provided clear evidence on both persistence of S. Dublin at single farms in consecutive years and transmission of the organisms between herds in different distances. A possible risk factor for the repeated occurrence of S. Dublin in certain districts of Schleswig-Holstein might be the spreading of manure on pastures and grassland. Effective control of S. Dublin requires farm-specific analysis of the management supplemented by WGS of outbreak causing S. Dublin strains to clearly identify routes of infection.

Oxford nanopore technologies-a valuable tool to generate whole-genome sequencing data for in silico serotyping and the detection of genetic markers in Salmonella.

Bacteria of the genus Salmonella pose a major risk to livestock, the food economy, and public health. Salmonella infections are one of the leading causes of food poisoning. The identification of serovars of Salmonella achieved by their diverse surface antigens is essential to gain information on their epidemiological context. Traditionally, slide agglutination has been used for serotyping. In recent years, whole-genome sequencing (WGS) followed by in silico serotyping has been established as an alternative method for serotyping and the detection of genetic markers for Salmonella. Until now, WGS data generated with Illumina sequencing are used to validate in silico serotyping methods. Oxford Nanopore Technologies (ONT) opens the possibility to sequence ultra-long reads and has frequently been used for bacterial sequencing. In this study, ONT sequencing data of 28 Salmonella strains of different serovars with epidemiological relevance in humans, food, and animals were taken to investigate the performance of the in silico serotyping tools SISTR and SeqSero2 compared to traditional slide agglutination tests. Moreover, the detection of genetic markers for resistance against antimicrobial agents, virulence, and plasmids was studied by comparing WGS data based on ONT with WGS data based on Illumina. Based on the ONT data from flow cell version R9.4.1, in silico serotyping achieved an accuracy of 96.4 and 92% for the tools SISTR and SeqSero2, respectively. Highly similar sets of genetic markers comparing both sequencing technologies were identified. Taking the ongoing improvement of basecalling and flow cells into account, ONT data can be used for Salmonella in silico serotyping and genetic marker detection.

Clonal relation between Salmonella enterica subspecies enterica serovar Dublin strains of bovine and food origin in Germany.

Salmonella enterica subspecies enterica serovar Dublin (S. Dublin) is a host-adapted serovar causing enteritis and/or systemic diseases in cattle. As the serovar is not host-restricted, it may cause infections in other animals, including humans with severe illness and higher mortality rates than other non-typhoidal serovars. As human infections are mainly caused by contaminated milk, milk products and beef, information on the genetic relationship of S. Dublin strains from cattle and food should be evaluated. Whole-genome sequencing (WGS) of 144 S. Dublin strains from cattle and 30 strains from food origin was performed. Multilocus sequence typing (MLST) revealed mostly sequence type ST-10 from both, cattle and food isolates. In total, 14 of 30 strains from food origin were clonally related to at least one strain from cattle, as detected by core-genome single nucleotide polymorphisms typing as well as core-genome MLST. The remaining 16 foodborne strains fit into the genome structure of S. Dublin in Germany without outliers. WGS proved to be a powerful tool not only to gain information on the epidemiology of Salmonella strains but also to detect clonal relations between organisms isolated from different stages of production. This study has shown a high genetic correlation between S. Dublin strains from cattle and food and, therefore, the potential to cause human infections. S. Dublin strains of both origins share an almost identical set of virulence factors, emphasizing their potential to cause severe clinical manifestations in animals, but also in humans and thus the need for effective control of S. Dublin in a farm-to-fork strategy.

Genomic insight into Campylobacter jejuni isolated from commercial turkey flocks in Germany using whole-genome sequencing analysis.

Campylobacter (C.) jejuni is a zoonotic bacterium of public health significance. The present investigation was designed to assess the epidemiology and genetic heterogeneity of C. jejuni recovered from commercial turkey farms in Germany using whole-genome sequencing. The Illumina MiSeq® technology was used to sequence 66 C. jejuni isolates obtained between 2010 and 2011 from commercial meat turkey flocks located in ten German federal states. Phenotypic antimicrobial resistance was determined. Phylogeny, resistome, plasmidome and virulome profiles were analyzed using whole-genome sequencing data. Genetic resistance markers were identified with bioinformatics tools (AMRFinder, ResFinder, NCBI and ABRicate) and compared with the phenotypic antimicrobial resistance. The isolates were assigned to 28 different sequence types and 11 clonal complexes. The average pairwise single nucleotide-polymorphisms distance of 14,585 SNPs (range: 0-26,540 SNPs) revealed a high genetic distinction between the isolates. Thirteen virulence-associated genes were identified in C. jejuni isolates. Most of the isolates harbored the genes flaA (83.3%) and flaB (78.8%). The wlaN gene associated with the Guillain-Barré syndrome was detected in nine (13.6%) isolates. The genes for resistance to ampicillin (bla OXA), tetracycline [tet(O)], neomycin [aph(3')-IIIa], streptomycin (aadE) and streptothricin (sat4) were detected in isolated C. jejuni using WGS. A gene cluster comprising the genes sat4, aph(3')-IIIa and aadE was present in six isolates. The single point mutation T86I in the housekeeping gene gyrA conferring resistance to quinolones was retrieved in 93.6% of phenotypically fluoroquinolone-resistant isolates. Five phenotypically erythromycin-susceptible isolates carried the mutation A103V in the gene for the ribosomal protein L22 inferring macrolide resistance. An assortment of 13 ß-lactam resistance genes (bla OXA variants) was detected in 58 C. jejuni isolates. Out of 66 sequenced isolates, 28 (42.4%) carried plasmid-borne contigs. Six isolates harbored a pTet-like plasmid-borne contig which carries the tet(O) gene. This study emphasized the potential of whole-genome sequencing to ameliorate the routine surveillance of C. jejuni. Whole-genome sequencing can predict antimicrobial resistance with a high degree of accuracy. However, resistance gene databases need curation and updates to revoke inaccuracy when using WGS-based analysis pipelines for AMR detection.