Search

ABSTRACTIn 2021, a white-handed gibbon (Hylobates lar) succumbed to illness shortly after transfer from one zoo to another in Germany, due to Francisella tularensis subsp. holarctica infection. To determine the source of infection, whole genome sequencing of the gibbon-derived isolate was performed and wild pest rodents (and captive squirrels) from both zoos were screened for F. tularensis. The F. tularensis whole genome sequence obtained from the gibbon was closely related to previous subclade B.281 sequences obtained from hares from Baden-Wuerttemberg, the same region where the gibbon was first housed. However, F. tularensis DNA was detected in one Norway rat from the receiving zoo. Therefore, neither zoo can be excluded as the source of infection.

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

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

BackgroundBovine mastitis is a common disease of dairy cattle causing major economic losses due to reduced yield and poor quality of milk worldwide. The current investigation aimed to gain insight into the genetic diversity, antimicrobial resistance profiles and virulence associated factors of Staphylococcus (S.) aureus isolated from clinical bovine mastitis in dairy farms in Thuringia, Germany.MethodsForty Staphylococcus aureus isolates collected from clinical bovine mastitis cases from 17 Thuringian dairy farms were phenotyped and genetically characterized using whole genome sequencing.ResultsOut of 40 S. aureus, 30 (75%) were confirmed as methicillin resistant isolates. The isolates showed elevated antimicrobial resistance against penicillin, tetracycline and oxacillin, i.e., 77.5, 77.5, and 75%, respectively. Lower resistance rates were found against moxifloxacin, ciprofloxacin, gentamicin and trimethoprim/sulfamethoxazole, i.e., 35, 35, 30, and 22.5%, respectively. While resistance against clindamycin and erythromycin was rarely found (5 and 2.5%, respectively). All isolates were susceptible to linezolid, teicoplanin, vancomycin, tigecycline, fosfomycin, fusidic acid and rifampicin. These isolates were further allocated into five different sequence types: ST398 (n = 31), ST1074 (n = 4), ST504 (n = 3), ST582 (CC15) (n = 1) and ST479 (n = 1). These isolates were also assigned to seven clusters with up to 100 SNP which has facilitated geographical mapping and epidemiological distribution in Thuringia. Strains belonging to ST398 were classified into clusters 1, 2, 3, 4 and 7. The isolates of ST504 were of cluster 5, those of ST1074 were belonging to cluster 6. Resistance genes blaZ, blaI and blaR associated with penicillin resistance were found in 32 (80%) strains, all except one were belonging to ST398. Methicillin resistance associated mecA was identified in 30 (96.8%) isolates of ST398. All tetracycline and erythromycin resistant isolates were of ST398, and all harbored both tetM and ermA. About 90.3% of tetracycline resistant isolates assigned to ST398 were also carrying tetK gene. The point mutations parC_S80F, gyrA_S84L and parC_S80Y in gyrA and parC associated with quinolone resistance were found in all phenotypically resistant isolates to ciprofloxacin and moxifloxacin (n = 14). Sixty-eight virulence genes were identified among isolates. Both lukD/E and lukM/F-PV-P83 were identified in 22.5% of isolates, all were non-ST398.ConclusionIn this study, ST398 had the highest potential to cause disease and had a massive prevalence in bovine mastitis cases. Five different sequence types and seven clusters were identified in the federal state of Thuringia. The circulation of some clusters in the same region over several years shows the persistence of cluster-associated infection despite the intensive veterinary care. On the other hand, some regions had different clusters at the same year or in different consecutive years. Different sequence types and associated different clusters of S. aureus were geographically widely distributed among dairy farms in Thuringia. The findings of this study show that various clusters have the potential to spread over a large geographical scale. The detection of LA-MRSA on dairy farms, which is known for cabapility to widely spread among different groups of animals, humans and their environment urges for the implementation of national wide strategic programs. The identification of CA-MRSA among the isolates such as ST398 poses a significant risk for the transmission of such strains between animals and humans on dairy farms.

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.

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.

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 (blaOXA), 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 (blaOXA 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.

Bacillus cereus is a ubiquitous spore-forming bacterium causing food spoilage. In this study, seventeen B. cereus were isolated from dairy processing equipment in Algeria and characterized. The pathogenic potential genes encoding hemolysin, nonhaemolytic enterotoxin, cytotoxin K and emetic toxin were investigated by PCR. Phylogeny and genetic markers for antimicrobial resistance and virulence were analyzed using whole-genome sequencing data using the Illumina MiSeq® technology. The identity of the isolated strains was confirmed as B. cereus sensu lato and the sequencing of the panC gene clustered 64.7% of them in phylogenetic group III and 35.3% in group IV. The results of PFGE showed that 8 (47.1%), 4 (23.5%), 3 (17.6%) and 2 (11.8%) were differentiated into four groups A, B, C and D, respectively. Each of the two panC phylogenetic groups contains two different types of isolates corresponding to PFGE group (A and C) and (B and D) for phylogenetic groups III and IV, respectively. None of the 14 isolates carried the emetic toxin (ces) gene, whereas other toxin genes were variably detected. Genotyping assigned all isolates to three different sequence types (ST 2226, ST1018 and ST1431) while one clonal complex (CC-142) was identified in three isolates. Hierarchical clustering with a threshold of 500 core genome single nucleotide polymorphism (cgSNPs) differences grouped the strains into three clusters while two strains were less than 20 SNPs different. Three B. cereus (25%) contained three enterotoxic HBL complex-encoding genes hblA/C/D. Enterotoxic NHE complex-encoding genes (nheA/B/C) and the ctyK1 gene were detected in all isolates. The overall results highlight the risk due to toxins produced by the B. cereus group in the food safety and dairy industries in Algeria. This study emphasized the potential of whole-genome sequencing for genotyping and predicting the virulence-associated genes and antimicrobial resistance. The present study will help to better assess the health and spoilage risk associated with B. cereus in dairy processing plants and to incorporate adequate preventive measures.

Campylobacter fetus subsp. venerealis (Cfv) causes bovine genital campylobacteriosis (BGC), a World Organization for Animal Health (WOAH)-listed trade-relevant disease characterized by severe reproductive losses, such as infertility, early embryonic death and abortion in cattle. BGC has significant economic implications that have prompted several countries to adopt stringent eradication and surveillance measures to contain the disease. In Germany, there has been a low incidence of BGC cases over the past 28 years. This study aimed to illustrate the genomic diversity of German Cfv strains isolated from different federal states in Germany. This study analyzed 63 Cfv strains, collected between 1985 and 2015, by whole-genome sequencing and compared them with genome data of 91 international Cfv isolates. The phylogenetic analysis showed that the Cfv population is genetically conserved and has geographic clusters. In Germany, one phylogenetic lineage comprising all strains was identified. This German lineage was part of a subclade that probably emerged in the nineteenth century and diversified over time. The results of this study point to a non-recurrent cross-border introduction of Cfv in Germany. The BGC control interventions in Germany can be considered successful as no outbreaks were reported since 2015.

Abstract Clostridium perfringens causes a plethora of devastating infections, with toxin production being the underlying mechanism of pathogenicity in various hosts. Genomic analyses of 206 public-available C. perfringens strains´ sequence data identified a substantial degree of genomic variability in respect to episome content, chromosome size and mobile elements. However, the position and order of the local collinear blocks on the chromosome showed a considerable degree of preservation. The strains were divided into five stable phylogroups (I–V). Phylogroup I contained human food poisoning strains with chromosomal enterotoxin (cpe) and a Darmbrand strain characterized by a high frequency of mobile elements, a relatively small genome size and a marked loss of chromosomal genes, including loss of genes encoding virulence traits. These features might correspond to the adaptation of these strains to a particular habitat, causing human foodborne illnesses. This contrasts strains that belong to phylogroup II where the genome size points to the acquisition of genetic material. Most strains of phylogroup II have been isolated from enteric lesions in horses and dogs. Phylogroups III, IV and V are heterogeneous groups containing a variety of different strains, with phylogroup III being the most abundant (65.5%). In conclusion, C. perfringens displays five stable phylogroups reflecting different disease involvements, prompting further studies on the evolution of this highly important pathogen.

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.

ABSTRACT Salmonella enterica subsp. enterica serovar Dublin is a cattle-adapted serovar that causes enteritis and systemic diseases in animals. In Germany, S. Dublin is not detected or is very rarely detected in some federal states but is endemic in certain regions. Information on detailed genetic characteristics of S. Dublin is not available. An understanding of the paths and spreading of S. Dublin within and between regions and over time is essential to establish effective control strategies. Whole-genome sequencing (WGS) and bioinformatic analysis were used to explore the genetic traits of S. Dublin and to determine their epidemiological context. Seventy-four S. Dublin strains collected in 2005 to 2018 from 10 federal states were studied. The phylogeny was analyzed using core-genome single-nucleotide polymorphisms (cgSNPs) and core-genome multilocus sequence typing. Genomic clusters at 100 cgSNPs, 40 cgSNPs, and 15 cgSNPs were selected for molecular epidemiology. WGS-based genoserotyping confirmed serotyping. Important specific virulence determinants were detected in all strains, but multidrug resistance in German S. Dublin organisms is uncommon. Use of different thresholds for cgSNP analysis enabled a broad view and also a detailed view of the occurrence of S. Dublin in Germany. Genomic clusters could be allocated nationwide, to a limited number of federal states, or to special regions only. Results indicate both persistence and spread of S. Dublin within and between federal states in short and longer time periods. However, to detect possible routes of infection or persistence of S. Dublin indicated by genomic analysis, information on the management of the cattle farms and contacts with corresponding farms are essential. IMPORTANCE Salmonella enterica subsp. enterica serovar Dublin is a bovine host-adapted serovar that causes up to 50% of all registered outbreaks of salmonellosis in cattle in Germany. S. Dublin is not detected or is only rarely detected in some federal states but has been endemic in certain regions of the country for a long time. Information on genetic traits of the causative strains is essential to determine routes of infection. WGS and bioinformatic analysis should be used to explore the genetic characteristics of S. Dublin. Combining the genomic features of S. Dublin strains with information on the management of the cattle farms concerned should enable the detection of possible routes of infection or persistence of S. Dublin. This approach is regarded as a prerequisite to developing effective intervention strategies.

ABSTRACT Clostridium perfringens is a spore-forming anaerobic pathogen responsible for a variety of histotoxic and intestinal infections in humans and animals. High-resolution genotyping aiming to identify bacteria at strain level has become increasingly important in modern microbiology to understand pathogen transmission pathways and to tackle infection sources. This study aimed at establishing a publicly available genome-wide multilocus sequence-typing (MLST) scheme for C. perfringens. A total of 1,431 highly conserved core genes (1.34 megabases; 50% of the reference genome genes) were indexed for a core genome-based MLST (cgMLST) scheme for C. perfringens. The scheme was applied to 282 ecologically and geographically diverse genomes, showing that the genotyping results of cgMLST were highly congruent with the core genome-based single-nucleotide-polymorphism typing in terms of resolution and tree topology. In addition, the cgMLST provided a greater discrimination than classical MLST methods for C. perfringens. The usability of the scheme for outbreak analysis was confirmed by reinvestigating published outbreaks of C. perfringens-associated infections in the United States and the United Kingdom. In summary, a publicly available scheme and an allele nomenclature database for genomic typing of C. perfringens have been established and can be used for broad-based and standardized epidemiological studies. IMPORTANCE Global epidemiological surveillance of bacterial pathogens is enhanced by the availability of standard tools and sharing of typing data. The use of whole-genome sequencing has opened the possibility for high-resolution characterization of bacterial strains down to the clonal and subclonal levels. Core genome multilocus sequence typing is a robust system that uses highly conserved core genes for deep genotyping. The method has been successfully and widely used to describe the epidemiology of various bacterial species. Nevertheless, a cgMLST typing scheme for Clostridium perfringens is currently not publicly available. In this study, we (i) developed a cgMLST typing scheme for C. perfringens, (ii) evaluated the performance of the scheme on different sets of C. perfringens genomes from different hosts and geographic regions as well as from different outbreak situations, and, finally, (iii) made this scheme publicly available supported by an allele nomenclature database for global and standard genomic typing.

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.

Burkholderia (B.) mallei is a host-adapted equine pathogen that causes glanders, a re-emerging zoonotic disease, which is endemic in Pakistan and other developing countries and seriously impacts the global equine movement. Due to globalization, the geographical restriction of diseases vanishes and the lack of awareness of and experience with eradicated diseases in industrialized countries also promotes the re-introduction of infections in these regions. Owing to the high equine population, the Pakistani province Punjab is a potential hotspot where several glanders outbreaks have been seen over last two decades. For determining the genomic diversity of B. mallei in this and other equine-populated prefectures, the genomes of 19 B. mallei strains isolated between 1999 and 2020 in different locations were sequenced and their genotypes were determined. Particularly, for genetically highly homogenous pathogens like B. mallei genotyping techniques require a high discriminatory power for enabling differentiation on the strain level. Thus, core-genome single nucleotide polymorphism (cgSNP) analysis was applied for distinguishing the highly similar strains. Furthermore, a whole-genome sequence-based core genome multi locus sequence typing (cgMLST) scheme, specific to B. mallei, was developed and additionally applied to the data. It was found that B. mallei genotypes in Pakistan persisted over time and space and genotype clusters preferred connection with a time point rather than the place of isolation, probably due to frequent equine movement, which promotes the spread of glanders. The cgMLST approach proved to work in accord with SNP typing and may help to investigate future glanders outbreaks.

Despite extensive monitoring programs and preventative measures, Salmonella spp. continue to cause tens of thousands human infections per year, as well as many regional and international food-borne outbreaks, that are of great importance for public health and cause significant socio-economic costs. In Germany, salmonellosis is the second most common cause of bacterial diarrhea in humans and is associated with high hospitalization rates. Whole-genome sequencing (WGS) combined with data analysis is a high throughput technology with an unprecedented discriminatory power, which is particularly well suited for targeted pathogen monitoring, rapid cluster detection and assignment of possible infection sources. However, an effective implementation of WGS methods for large-scale microbial pathogen detection and surveillance has been hampered by the lack of standardized methods, uniform quality criteria and strategies for data sharing, all of which are essential for a successful interpretation of sequencing data from different sources. To overcome these challenges, the national GenoSalmSurv project aims to establish a working model for an integrated genome-based surveillance system of Salmonella spp. in Germany, based on a decentralized data analysis. Backbone of the model is the harmonization of laboratory procedures and sequencing protocols, the implementation of open-source bioinformatics tools for data analysis at each institution and the establishment of routine practices for cross-sectoral data sharing for a uniform result interpretation. With this model, we present a working solution for cross-sector interpretation of sequencing data from different sources (such as human, veterinarian, food, feed and environmental) and outline how a decentralized data analysis can contribute to a uniform cluster detection and facilitate outbreak investigations.

Aliarcobacter butzleri is an emerging foodborne and zoonotic pathogen that is usually transmitted via contaminated food or water. A. butzleri is not only the most prevalent Aliarcobacter species, it is also closely related to thermophilic Campylobacter, which have shown increasing resistance in recent years. Therefore, it is important to assess its resistance and virulence profiles. In this study, 45 Aliarcobacter butzleri strains from water poultry farms in Thuringia, Germany, were subjected to an antimicrobial susceptibility test using the gradient strip diffusion method and whole-genome sequencing. In the phylogenetic analysis, the genomes of the German strains showed high genetic diversity. Thirty-three isolates formed 11 subgroups containing two to six strains. The antimicrobial susceptibility testing showed that 32 strains were resistant to erythromycin, 26 to doxycycline, and 20 to tetracycline, respectively. Only two strains were resistant to ciprofloxacin, while 39 strains were resistant to streptomycin. The in silico prediction of the antimicrobial resistance profiles identified a large repertoire of potential resistance mechanisms. A strong correlation between a gyrA point mutation (Thr-85-Ile) and ciprofloxacin resistance was found in 11 strains. A partial correlation was observed between the presence of the bla3 gene and ampicillin resistance. In silico virulence profiling revealed a broad spectrum of putative virulence factors, including a complete lipid A cluster in all studied genomes.

Subspecies of the species Campylobacter fetus are associated with specific host niches including mammals and reptiles. Campylobacter fetus subsp. fetus is a zoonotic pathogen infecting humans. Infections can vary from an acute intestinal illness to severe systemic infections, with sheep and cattle as major reservoirs. In contrast, Campylobacter fetus subsp. venerealis causes bovine genital campylobacteriosis, which leads to abortion in cattle and a high economic burden for the farmers. Therefore, high-quality molecular subtyping is indispensable for interventional epidemiology. We used whole-genome sequencing (WGS) data of 283 Campylobacter fetus strains from 18 countries and compared several methods for Campylobacter fetus subtyping, including WGS, multilocus sequence typing, PCR assays, and the presence of the insertion element ISCfe1. We identified a highly clonal clade (designated as clade 1) that harbors the insertion sequence ISCfe1. The presence of this insertion sequence is an essential diagnostic tool for the identification of the subspecies Campylobacter fetus subsp. venerealis, serving as a target for several PCR assays. However, we have found a high sequence variability for the ISCfe1 besides the presence of ISCfe1-paralogues in certain other genomes (n = 7) which may cause incorrect diagnostic results. Clade 1 seems to be the cattle-specific clade of this species. We propose that only this clade might be designated as Campylobacter fetus subsp. venerealis as it harbors the ISCfe1 marker sequence, which is a major target for molecular methods currently used for Campylobacter fetus subspecies identification. Fostering this proposal, we defined eleven stable nucleotide markers specific for this clade. Additionally, we developed a bioinformatics toolbox for the fast identification of this clade based on WGS data. In conclusion, our results demonstrate that WGS can be used for Campylobacter fetus subtyping overcoming limitations of current PCR and MLST protocols.

Summary: High-throughput RNA sequencing (RNA-seq) is routinely applied to study diverse biological processes; however, when performed separately on interacting organisms, systemic noise intrinsic to RNA extraction, library preparation, and sequencing hampers the identification of cross-species interaction nodes. Here, we develop triple RNA-seq to simultaneously detect transcriptomes of monocyte-derived dendritic cells (moDCs) infected with the frequently co-occurring pulmonary pathogens Aspergillus fumigatus and human cytomegalovirus (CMV). Comparing expression patterns after co-infection with those after single infections, our data reveal synergistic effects and mutual interferences between host responses to the two pathogens. For example, CMV attenuates the fungus-mediated activation of pro-inflammatory cytokines through NF-κB (nuclear factor κB) and NFAT (nuclear factor of activated T cells) cascades, while A. fumigatus impairs viral clearance by counteracting viral nucleic acid-induced activation of type I interferon signaling. Together, the analytical power of triple RNA-seq proposes molecular hubs in the differential moDC response to fungal/viral single infection or co-infection that contribute to our understanding of the etiology and, potentially, clearance of post-transplant infections.

ABSTRACT Only four species, Candida albicans, C. glabrata, C. parapsilosis, and C. tropicalis, together account for about 90% of all Candida bloodstream infections and are among the most common causes of invasive fungal infections of humans. However, virulence potential varies among these species, and the phylogenetic tree reveals that their pathogenicity may have emerged several times independently during evolution. We therefore tested these four species in a human whole-blood infection model to determine, via comprehensive dual-species RNA-sequencing analyses, which fungal infection strategies are conserved and which are recent evolutionary developments. The ex vivo infection progressed from initial immune cell interactions to nearly complete killing of all fungal cells. During the course of infection, we characterized important parameters of pathogen-host interactions, such as fungal survival, types of interacting immune cells, and cytokine release. On the transcriptional level, we obtained a predominantly uniform and species-independent human response governed by a strong upregulation of proinflammatory processes, which was downregulated at later time points after most of the fungal cells were killed. In stark contrast, we observed that the different fungal species pursued predominantly individual strategies and showed significantly different global transcriptome patterns. Among other findings, our functional analyses revealed that the fungal species relied on different metabolic pathways and virulence factors to survive the host-imposed stress. These data show that adaptation of Candida species as a response to the host is not a phylogenetic trait, but rather has likely evolved independently as a prerequisite to cause human infections. IMPORTANCE To ensure their survival, pathogens have to adapt immediately to new environments in their hosts, for example, during the transition from the gut to the bloodstream. Here, we investigated the basis of this adaptation in a group of fungal species which are among the most common causes of hospital-acquired infections, the Candida species. On the basis of a human whole-blood infection model, we studied which genes and processes are active over the course of an infection in both the host and four different Candida pathogens. Remarkably, we found that, while the human host response during the early phase of infection is predominantly uniform, the pathogens pursue largely individual strategies and each one regulates genes involved in largely disparate processes in the blood. Our results reveal that C. albicans, C. glabrata, C. parapsilosis, and C. tropicalis all have developed individual strategies for survival in the host. This indicates that their pathogenicity in humans has evolved several times independently and that genes which are central for survival in the host for one species may be irrelevant in another.

The emergence of Klebsiella pneumoniae (K. pneumoniae) in German healthcare is worrying. It is not well-investigated in the veterinary world and food chains. In the current study, antibiotic susceptibility profiles of 24 K. pneumoniae strains isolated from powdered milk samples produced in Germany were investigated by a microdilution test. Next-generation sequencing (NGS) was applied to identify genomic determinants for antimicrobial resistance (AMR), virulence-associated genes and plasmids replicons. All isolates were susceptible to the majority (14/18) of tested antibiotics. Resistance to colistin, fosfomycin, chloramphenicol and piperacillin was found. The ambler class A ß-lactamase, blaSHV variants were identified in all isolates, of which blaSHV-187 was most prevalent and found in 50% of isolates. Single-nucleotide-variants of oqxA and oqxB conferring resistance to phenicol/quinolone were found in all isolates, and the oqxB17 was the most prevalent found in 46% of isolates. 67% of isolates harbored fosA genes; however, only one was fosfomycin-resistant. Two isolates harbored genes conferring resistance to colistin, despite being susceptible. The majority of identified virulome genes were iron uptake siderophores. Two enterobactins (entB, fepC), six adherence-related genes belonging to E. coli common pilus (ECP) and one secretion system (ompA gene) were found in all isolates. In contrast, yersiniabactin was found in two isolates. One ST23 strain was susceptible to all tested antibiotics, and harbored determinants discriminatory for hypervirulent strains, e.g., aerobactin, salmochelin, yersiniabactin, enterobactin and regulator of mucoid phenotype A genes that are highly associated with hypervirulent K. pneumoniae. The IncF plasmid family was found in all strains, while almost half of the isolates harbored Col440I-type plasmids and nine isolates harbored various Inc-type plasmids. The presence of K. pneumoniae carrying different resistomes and major virulent specific virulomes in powdered milk samples is alarming. This could threaten public health, particularly of neonates and infants consuming dried milk.

Anthrax is a recurrent zoonosis in the Ukraine with outbreaks occurring repeatedly in certain areas. For determining whether several Bacillus anthracis genotypes are circulating in this region, four strains from various sources isolated from different regions of the Ukraine were investigated. By combining long- and short-read next-generation sequencing techniques, highly accurate genomes were reconstructed, enabling detailed in silico genotyping. Thus, the strains could be assigned to the Tsiankovskii subgroup of the “TransEurAsia” clade, which is commonly found in this region. Their high genetic similarity suggests that the four strains are members of the endemic population whose progenitor was once introduced in the Ukraine and bordering regions. This study provides information on B. anthracis strains from a region where there is little knowledge of the local population, thereby adding to the picture of global B. anthracis genotype distribution. We also emphasize the importance of surveillance and prevention methods regarding anthrax outbreaks, as other studies predicted a higher number of cases in the future due to global warming.

Delayed natural killer (NK) cell reconstitution after allogeneic stem cell transplantation (alloSCT) is associated with a higher risk of developing invasive aspergillosis. The interaction of NK cells with the human pathogen Aspergillus (A.) fumigatus is mediated by the fungal recognition receptor CD56, which is relocated to the fungal interface after contact. Blocking of CD56 signaling inhibits the fungal mediated chemokine secretion of MIP-1α, MIP-1β, and RANTES and reduces cell activation, indicating a functional role of CD56 in fungal recognition. We collected peripheral blood from recipients of an allograft at defined time points after alloSCT (day 60, 90, 120, 180). NK cells were isolated, directly challenged with live A. fumigatus germ tubes, and cell function was analyzed and compared to healthy age and gender-matched individuals. After alloSCT, NK cells displayed a higher percentage of CD56brightCD16dim cells throughout the time of blood collection. However, CD56 binding and relocalization to the fungal contact side were decreased. We were able to correlate this deficiency to the administration of corticosteroid therapy that further negatively influenced the secretion of MIP-1α, MIP-1β, and RANTES. As a consequence, the treatment of healthy NK cells ex vivo with corticosteroids abrogated chemokine secretion measured by multiplex immunoassay. Furthermore, we analyzed NK cells regarding their actin cytoskeleton by Structured Illumination Microscopy (SIM) and flow cytometry and demonstrate an actin dysfunction of NK cells shown by reduced F-actin content after fungal co-cultivation early after alloSCT. This dysfunction remains until 180 days post-alloSCT, concluding that further actin-dependent cellular processes may be negatively influenced after alloSCT. To investigate the molecular pathomechansism, we compared CD56 receptor mobility on the plasma membrane of healthy and alloSCT primary NK cells by single-molecule tracking. The results were very robust and reproducible between tested conditions which point to a different molecular mechanism and emphasize the importance of proper CD56 mobility.

During the last decade, Salmonella enterica subspecies enterica serovar Infantis (S. Infantis) has become more prevalent across Europe with an increased capability to persist in broiler farms. In this study, we aimed to identify potential genetic causes for the increased emergence and longer persistence of S. Infantis in German poultry farms by high-throughput-sequencing. Broiler derived S. Infantis strains from two decades, the 1990s (n = 12) and the 2010s (n = 18), were examined phenotypically and genotypically to detect potential differences responsible for increased prevalence and persistence. S. Infantis organisms were characterized by serotyping and determining antimicrobial susceptibility using the microdilution method. Genotypic characteristics were analyzed by whole genome sequencing (WGS) to detect antimicrobial resistance and virulence genes as well as plasmids. To detect possible clonal relatedness within S. Infantis organisms, 17 accessible genomes from previous studies about emergent S. Infantis were downloaded and analyzed using complete genome sequence of SI119944 from Israel as reference. In contrast to the broiler derived antibiotic-sensitive S. Infantis strains from the 1990s, the majority of strains from the 2010s (15 out of 18) revealed a multidrug-resistance (MDR) phenotype that encodes for at least three antimicrobials families: aminoglycosides [ant(3“)-Ia], sulfonamides (sul1), and tetracyclines [tet(A)]. Moreover, these MDR strains carry a virulence gene pattern missing in strains from the 1990s. It includes genes encoding for fimbriae clusters, the yersiniabactin siderophore, mercury and disinfectants resistance and toxin/antitoxin complexes. In depth genomic analysis confirmed that the 15 MDR strains from the 2010s carry a pESI-like megaplasmid with resistance and virulence gene patterns detected in the emerged S. Infantis strain SI119944 from Israel and clones inside and outside Europe. Genotyping analysis revealed two sequence types (STs) among the resistant strains from the 2010s, ST2283 (n = 13) and ST32 (n = 2). The sensitive strains from the 1990s, belong to sequence type ST32 (n = 10) and ST1032 (n = 2). Therefore, this study confirms the emergence of a MDR S. Infantis pESI-like clone of ST2283 in German broiler farms with presumably high tendency of dissemination. Further studies on the epidemiology and control of S. Infantis in broilers are needed to prevent the transfer from poultry into the human food chain.

Brucellosis, caused by the bacteria of the genus Brucella, is one of the most neglected common zoonotic diseases globally with a public health significance and a high economic loss among the livestock industry worldwide. Since little is known about the distribution of B. abortus in Egypt, a total of 46 B. abortus isolates recovered between 2012–2020, plus one animal isolate from 2006, were analyzed by examining the whole core genome single nucleotide polymorphism (cgSNP) in comparison to the in silico multilocus variable number of tandem repeat analysis (MLVA). Both cgSNP analysis and MLVA revealed three clusters and one isolate only was distantly related to the others. One cluster identified a rather widely distributed outbreak strain which is repeatedly occurring for at least 16 years with marginal deviations in cgSNP analysis. The other cluster of isolates represents a rather newly introduced outbreak strain. A separate cluster comprised RB51 vaccine related strains, isolated from aborted material. The comparison with MLVA data sets from public databases reveals one near relative from Argentina to the oldest outbreak strain and a related strain from Spain to a newly introduced outbreak strain in Egypt. The distantly related isolate matches with a strain from Portugal in the MLVA profile. Based on cgSNP analysis the oldest outbreak strain clusters with strains from the UK. Compared to the in silico analysis of MLVA, cgSNP analysis using WGS data provides a much higher resolution of genotypes and, when correlated to the associated epidemiological metadata, cgSNP analysis allows the differentiation of outbreaks by defining different outbreak strains. In this respect, MLVA data are error-prone and can lead to incorrect interpretations of outbreak events.

Brucellosis is a highly contagious zoonosis that occurs worldwide. Whole-genome sequencing (WGS) has become a widely accepted molecular typing method for outbreak tracing and genomic epidemiology of brucellosis. Twenty-nine Brucella spp. (eight B. abortus biovar 1 and 21 B. melitensis biovar 3) were isolated from lymph nodes, milk, and fetal abomasal contents of infected cattle, buffaloes, sheep, and goats originating from nine districts in Egypt. The isolates were identified by microbiological methods and matrix-assisted laser desorption ionization time-of-flight mass spectrometry (MALDI-TOF MS). Differentiation and genotyping were confirmed using multiplex PCR. Illumina MiSeq® was used to sequence the 29 Brucella isolates. Using MLST typing, ST11 and ST1 were identified among B. melitensis and B. abortus, respectively. Brucella abortus and B. melitensis isolates were divided into two main clusters (clusters 1 and 2) containing two and nine distinct genotypes by core-genome SNP analysis, respectively. The genotypes were irregularly distributed over time and space in the study area. Both Egyptian B. abortus and B. melitensis isolates proved to be genomically unique upon comparison with publicly available sequencing from strains of neighboring Mediterranean, African, and Asian countries. The antimicrobial resistance mechanism caused by mutations in rpoB, gyrA, and gyrB genes associated with rifampicin and ciprofloxacin resistance were identified. To the best of our knowledge, this is the first study investigating the epidemiology of Brucella isolates from livestock belonging to different localities in Egypt based on whole genome analysis.

Acinetobacter baumannii (A. baumannii) is a major cause of severe nosocomial infections worldwide. The emergence of infections associated with A. baumannii poses a significant health risk in Germany. A. baumannii is part of the ACB complex and is difficult to distinguish from other species phenotypically, necessitating its reliable identification. The current study analyzed 89 A. baumannii strains from human and non-human origins by matrix-assisted laser desorption/ionization (MALDI–TOF) and PCR detection of intrinsic blaOXA-51-like carbapenemase, blaOXA-23-like, blaOXA-24-like, blaOXA-58-like, and ISAba 1 genes. Whole-genome sequencing (WGS) was applied for species confirmation and strain type determination. Combining the molecular detection of the intrinsic blaOXA-51-like carbapenemase gene together with MALDI–TOF with a score value of >2.300 proved to be a suitable tool for A. baumannii identification. WGS data for all of the sequenced strains confirmed the identity of all A. baumannii strains. The Pasteur scheme successfully assigned 79.7% of the strains into distinct STs, while the Oxford scheme succeeded in allocating only 42.7% of isolates. Multilocus sequence typing (MLST) analysis based on the Pasteur scheme identified 16 STs. ST/241 was the most prevalent in samples from non-human origin, whereas ST/2 was predominant in human samples. Furthermore, eight isolates of non-human origin were allocated to seven new STs (ST/1410, ST/1414, ST/1416, ST/1417, ST/1418, ST/1419, and ST/1421). Ten isolates from non-human origin could not be typed since new alleles were observed in the loci Pas_cpn60, Pas_rpoB, and Pas_gltA. MLST analysis based on the Pasteur scheme was more appropriate than the Oxford scheme for the current group of A. baumannii.

Carbapenem-resistant Acinetobacter baumannii (A. baumannii, CRAb) is an emerging global threat for healthcare systems, particularly in Southeast Asia. Next-generation sequencing (NGS) technology was employed to map genes associated with antimicrobial resistance (AMR) and to identify multilocus sequence types (MLST). Eleven strains isolated from humans in Vietnam were sequenced, and their AMR genes and MLST were compared to published genomes of strains originating from Southeast Asia, i.e., Thailand (n = 49), Myanmar (n = 38), Malaysia (n = 11), Singapore (n = 4) and Taiwan (n = 1). Ten out of eleven Vietnamese strains were CRAb and were susceptible only to colistin. All strains harbored ant(3”)-IIa, armA, aph(6)-Id and aph(3”) genes conferring resistance to aminoglycosides, and blaOXA-51 variants and blaADC-25 conferring resistance to ß-lactams. More than half of the strains harbored genes that confer resistance to tetracyclines, sulfonamides and macrolides. The strains showed high diversity, where six were assigned to sequence type (ST)/2, and two were allocated to two new STs (ST/1411-1412). MLST analyses of 108 strains from Southeast Asia identified 19 sequence types (ST), and ST/2 was the most prevalent found in 62 strains. A broad range of AMR genes was identified mediating resistance to ß-lactams, including cephalosporins and carbapenems (e.g., blaOXA-51-like, blaOXA-23, blaADC-25, blaADC-73, blaTEM-1, blaNDM-1), aminoglycosides (e.g., ant(3”)-IIa, aph(3”)-Ib, aph(6)-Id, armA and aph(3’)-Ia), phenicoles (e.g., catB8), tetracyclines (e.g., tet.B and tet.39), sulfonamides (e.g., sul.1 and sul.2), macrolides and lincosamide (e.g., mph.E, msr.E and abaF). MLST and core genome MLST (cgMLST) showed an extreme diversity among the strains. Several strains isolated from different countries clustered together by cgMLST; however, different clusters shared the same ST. Developing an action plan on AMR, increasing awareness and prohibiting the selling of antibiotics without prescription must be mandatory for this region. Such efforts are critical for enforcing targeted policies on the rational use of carbapenem compounds and controlling AMR dissemination and emergence in general.

Alternative splicing (AS) is an important regulatory mechanism in eukaryotes but only little is known about its impact in fungi. Human fungal pathogens are of high clinical interest causing recurrent or life-threatening infections. AS can be well-investigated genome-wide and quantitatively with the powerful technology of RNA-Seq. Here, we systematically studied AS in human fungal pathogens based on RNA-Seq data. To do so, we investigated its effect in seven fungi during conditions simulating ex vivo infection processes and during in vitro stress. Genes undergoing AS are species-specific and act independently from differentially expressed genes pointing to an independent mechanism to change abundance and functionality. Candida species stand out with a low number of introns with higher and more varying lengths and more alternative splice sites. Moreover, we identified a functional difference between response to host and other stress conditions: During stress, AS affects more genes and is involved in diverse regulatory functions. In contrast, during response-to-host conditions, genes undergoing AS have membrane functionalities and might be involved in the interaction with the host. We assume that AS plays a crucial regulatory role in pathogenic fungi and is important in both response to host and stress conditions.

Schizophyllum commune is a filamentous basidiomycete causing white-rot in many wood species with the help of a broad range of enzymes including multicopper oxidases such as laccases and laccase-like oxidases. Since these enzymes exhibit a broad substrate range, their ability to oxidatively degrade black slate was investigated. Both haploid monokaryotic, and mated dikaryotic strains were able to grow on black slate rich in organic carbon as sole carbon source. On defined media, only the monokaryon showed growth promotion by addition of slate. At the same time, metals were released from the slate and, after reaching a threshold concentration, inhibited further growth of the fungus. The proteome during decomposition of the black slate showed induction of proteins potentially involved in rock degradation and stress resistance, and the gene for laccase-like oxidase mco2 was up-regulated. Specifically in the dikaryon, the laccase gene lcc1 was induced, while lcc2 as well as mco1, mco3, and mco4 expression levels remained similar. Spectrophotometric analysis revealed that both life forms were able to degrade the rock and produce smaller particles.

Francisella tularensis is the causative agent of the zoonotic disease tularemia. In Germany, most human infections are caused by contact with infected hares. The aim of this study was to characterize Francisella tularensis subsp. holarctica strains isolated from hares in Germany and to develop bioinformatics tools to analyze their genetic relatedness. In total, 257 German isolates—obtained mainly from hares (n = 233), other vertebrate animals, and ticks, but also from humans (n = 3)—were analyzed within this study. Publically available sequence data from 49 isolates were used to put our isolates into an epidemiological context and to compare isolates from natural foci and humans. Whole-genome sequences were analyzed using core-genome Multi-Locus-Sequence-Typing, canonical Single Nucleotide Polymorphism (SNP) typing and whole-genome SNP typing. An overall conformity of genotype clustering between the typing methods was found, albeit with a lower resolution for canonical single SNP typing. The subclade distribution, both on local and national levels, among strains from humans and hares was similar, suggesting circulation of the same genotypes both in animals and humans. Whilst close to identical isolates of the same subclade were found distributed over large areas, small geographical foci often harbored members of different subclades. In conclusion, although genomic high-resolution typing was shown to be robust, reproducible and allowed the identification of highly closely related strains, genetic profiling alone is not always conclusive for epidemiological linkage of F. tularensis strains.

Anthrax is a serious infectious disease caused by the gram-positive and spore-forming bacterium Bacillus anthracis. In Italy, anthrax is an endemic disease with sporadic cases each year and few outbreaks, especially in Southern Italy. However, new foci have been discovered in zones without previous history of anthrax. During summer 2016, an outbreak of anthrax caused the death of four goats in the Abruzzo region, where the disease had not been reported before. In order to investigate the outbreak, we sequenced one strain and compared it to 19 Italian B. anthracis genomes. Furthermore, we downloaded 71 whole genome sequences representing the global distribution of canonical SNP lineages and used them to verify the phylogenetic positioning. To this end, we analyzed and compared the genome sequences using canonical SNPs and the whole genome SNP-based analysis. Our results demonstrate that the outbreak strain belonged to the Trans-Eurasian (TEA) group A.Br.011/009, which is the predominant clade in Central-Southern Italy. In conclusion, the high genomic relatedness of the Italian TEA strains suggests their evolution from a common ancestor, while the spread is supposedly driven by trade as well as human and transhumance activities. Here, we demonstrated the capabilities of whole genome sequencing (WGS), which can be used as a tool for outbreak analyses and surveillance activities.

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.

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

ABSTRACT Iron is an essential micronutrient for both pathogens and their hosts, which restrict iron availability during infections in an effort to prevent microbial growth. Successful human pathogens like the yeast Candida glabrata have thus developed effective iron acquisition strategies. Their regulation has been investigated well for some pathogenic fungi and in the model organism Saccharomyces cerevisiae, which employs an evolutionarily derived system. Here, we show that C. glabrata uses a regulation network largely consisting of components of the S. cerevisiae regulon but also of elements of other pathogenic fungi. Specifically, similarly to baker’s yeast, Aft1 is the main positive regulator under iron starvation conditions, while Cth2 degrades mRNAs encoding iron-requiring enzymes. However, unlike the case with S. cerevisiae, a Sef1 ortholog is required for full growth under iron limitation conditions, making C. glabrata an evolutionary intermediate to SEF1-dependent fungal pathogens. Therefore, C. glabrata has evolved an iron homeostasis system which seems to be unique within the pathogenic fungi. IMPORTANCE The fungus Candida glabrata represents an evolutionarily close relative of the well-studied and benign baker’s yeast and model organism Saccharomyces cerevisiae. On the other hand, C. glabrata is an important opportunistic human pathogen causing both superficial and systemic infections. The ability to acquire trace metals, in particular, iron, and to tightly regulate this process during infection is considered an important virulence attribute of a variety of pathogens. Importantly, S. cerevisiae uses a highly derivative regulatory system distinct from those of other fungi. Until now, the regulatory mechanism of iron homeostasis in C. glabrata has been mostly unknown. Our study revealed a hybrid iron regulation network that is unique to C. glabrata and is placed at an evolutionary midpoint between those of S. cerevisiae and related fungal pathogens. We thereby show that, in the host, even a successful human pathogen can rely largely on a strategy normally found in nonpathogenic fungi from a terrestrial environment.

Severe influenza associated with strong symptoms and lung inflammation can be caused by intra-host evolution of quasispecies with aspartic acid or glycine in hemagglutinin position 222 (HA-222D/G; H1 numbering). To gain insights into the dynamics of host response to this coevolution and to identify key mechanisms contributing to copathogenesis, the lung transcriptional response of BALB/c mice infected with an A(H1N1)pdm09 isolate consisting HA-222D/G quasispecies was analyzed from day 1 to 12 post infection (p.i). At day 2 p.i. 968 differentially expressed genes (DEGs) were detected. The DEG number declined to 359 at day 4 and reached 1001 at day 7 p.i. prior to recovery. Interestingly, a biphasic expression profile was shown for the majority of these genes. Cytokine assays confirmed these results on protein level exemplarily for two key inflammatory cytokines, interferon gamma and interleukin 6. Using a reverse engineering strategy, a regulatory network was inferred to hypothetically explain the biphasic pattern for selected DEGs. Known regulatory interactions were extracted by Pathway Studio 9.0 and integrated during network inference. The hypothetic gene regulatory network revealed a positive feedback loop of Ifng, Stat1, and Tlr3 gene signaling that was triggered by the HA-G222 variant and correlated with a clinical symptom score indicating disease severity.

Mitogen activated protein kinases (MAPKs) are highly conserved in eukaryotic organisms. In pathogenic fungi, their activities were assigned to different physiological functions including drug adaptation and resistance. Aspergillus fumigatus is a human pathogenic fungus, which causes life-threatening invasive infections. Therapeutic options against invasive mycoses are still limited. One of the clinically used drugs is caspofungin, which specifically targets the fungal cell wall biosynthesis. A systems biology approach, based on comprehensive transcriptome data sets and mathematical modeling, was employed to infer a regulatory network and identify key interactions during adaptation to caspofungin stress in A. fumigatus. Mathematical modeling and experimental validations confirmed an intimate cross talk occurring between the cell wall-integrity and the high osmolarity-glycerol signaling pathways. Specifically, increased concentrations of caspofungin promoted activation of these signalings. Moreover, caspofungin affected the intracellular transport, which caused an additional osmotic stress that is independent of glucan inhibition. High concentrations of caspofungin reduced this osmotic stress, and thus decreased its toxic activity. Our results demonstrated that MAPK signaling pathways play a key role during caspofungin adaptation and are contributing to the paradoxical effect exerted by this drug.

The white-rot fungus Schizophyllum commune (Agaricomycetes) was used to study the cell biology of microtubular trafficking during mating interactions, when the two partners exchange nuclei, which are transported along microtubule tracks. For this transport activity, the motor protein dynein is required. In S. commune, the dynein heavy chain is encoded in two parts by two separate genes, dhc1 and dhc2. The N-terminal protein Dhc1 supplies the dimerization domain, while Dhc2 encodes the motor machinery and the microtubule binding domain. This split motor protein is unique to Basidiomycota, where three different sequence patterns suggest independent split events during evolution. To investigate the function of the dynein heavy chain, the gene dhc1 and the motor domain in dhc2 were deleted. Both resulting mutants were viable, but revealed phenotypes in hyphal growth morphology and mating behavior as well as in sexual development. Viability of strain Δdhc2 is due to the higher expression of kinesin-2 and kinesin-14, which was proven via RNA sequencing.

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.

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 (

Nitrogen is one of the key nutrients for microbial growth. During infection, pathogenic fungi like C. albicans need to acquire nitrogen from a broad range of different and changing sources inside the host. Detecting the available nitrogen sources and adjusting the expression of genes for their uptake and degradation is therefore crucial for survival and growth as well as for establishing an infection. Here, we analyzed the transcriptional response of C. albicans to nitrogen starvation and feeding with the infection-relevant nitrogen sources arginine and bovine serum albumin (BSA), representing amino acids and proteins, respectively. The response to nitrogen starvation was marked by an immediate repression of protein synthesis and an up-regulation of general amino acid permeases, as well as an up-regulation of autophagal processes in its later stages. Feeding with arginine led to a fast reduction in expression of general permeases for amino acids and to resumption of protein synthesis. The response to BSA feeding was generally slower, and was additionally characterized by an up-regulation of oligopeptide transporter genes. From time-series data, we inferred network interaction models for genes relevant in nitrogen detection and uptake. Each individual network was found to be largely specific for the experimental condition (starvation or feeding with arginine or BSA). In addition, we detected several novel connections between regulator and effector genes, with putative roles in nitrogen uptake. We conclude that C. albicans adopts a particular nitrogen response network, defined by sets of specific gene-gene connections for each environmental condition. All together, they form a grid of possible gene regulatory networks, increasing the transcriptional flexibility of C. albicans.

Candida albicans and Candida dubliniensis are pathogenic fungi that are highly related but differ in virulence and in some phenotypic traits. During in vitro growth on certain nutrient-poor media, C. albicans and C. dubliniensis are the only yeast species which are able to produce chlamydospores, large thick-walled cells of unknown function. Interestingly, only C. dubliniensis forms pseudohyphae with abundant chlamydospores when grown on Staib medium, while C. albicans grows exclusively as a budding yeast. In order to further our understanding of chlamydospore development and assembly, we compared the global transcriptional profile of both species during growth in liquid Staib medium by RNA sequencing. We also included a C. albicans mutant in our study which lacks the morphogenetic transcriptional repressor Nrg1. This strain, which is characterized by its constitutive pseudohyphal growth, specifically produces masses of chlamydospores in Staib medium, similar to C. dubliniensis. This comparative approach identified a set of putatively chlamydospore-related genes. Two of the homologous C. albicans and C. dubliniensis genes (CSP1 and CSP2) which were most strongly upregulated during chlamydospore development were analysed in more detail. By use of the green fluorescent protein as a reporter, the encoded putative cell wall related proteins were found to exclusively localize to C. albicans and C. dubliniensis chlamydospores. Our findings uncover the first chlamydospore specific markers in Candida species and provide novel insights in the complex morphogenetic development of these important fungal pathogens.