Machine Learning Algorithms for Classification of MALDI-TOF MS Spectra from Phylogenetically Closely Related Species Brucella melitensis, Brucella abortus and Brucella suis.

(1) Background: MALDI-TOF mass spectrometry (MS) is the gold standard for microbial fingerprinting, however, for phylogenetically closely related species, the resolution power drops down to the genus level. In this study, we analyzed MALDI-TOF spectra from 44 strains of B. melitensis, B. suis and B. abortus to identify the optimal classification method within popular supervised and unsupervised machine learning (ML) algorithms. (2) Methods: A consensus feature selection strategy was applied to pinpoint from among the 500 MS features those that yielded the best ML model and that may play a role in species differentiation. Unsupervised k-means and hierarchical agglomerative clustering were evaluated using the silhouette coefficient, while the supervised classifiers Random Forest, Support Vector Machine, Neural Network, and Multinomial Logistic Regression were explored in a fine-tuning manner using nested k-fold cross validation (CV) with a feature reduction step between the two CV loops. (3) Results: Sixteen differentially expressed peaks were identified and used to feed ML classifiers. Unsupervised and optimized supervised models displayed excellent predictive performances with 100% accuracy. The suitability of the consensus feature selection strategy for learning system accuracy was shown. (4) Conclusion: A meaningful ML approach is here introduced, to enhance Brucella spp. classification using MALDI-TOF MS data.

In-Depth Analysis of Bacillus anthracis 16S rRNA Genes and Transcripts Reveals Intra- and Intergenomic Diversity and Facilitates Anthrax Detection.

Analysis of 16S rRNA (rRNA) genes provides a central means of taxonomic classification of bacterial species. Based on presumed sequence identity among species of the Bacillus cereus sensu lato group, the 16S rRNA genes of B. anthracis have been considered unsuitable for diagnosis of the anthrax pathogen. With the recent identification of a single nucleotide polymorphism in some 16S rRNA gene copies, specific identification of B. anthracis becomes feasible. Here, we designed and evaluated a set of in situ, in vitro, and in silico assays to assess the unknown 16S state of B. anthracis from different perspectives. Using a combination of digital PCR, fluorescence in situ hybridization, long-read genome sequencing, and bioinformatics, we were able to detect and quantify a unique 16S rRNA gene allele of B. anthracis (16S-BA-allele). This allele was found in all available B. anthracis genomes and may facilitate differentiation of the pathogen from any close relative. Bioinformatics analysis of 959 B. anthracis SRA data sets inferred that abundances and genomic arrangements of the 16S-BA-allele and the entire rRNA operon copy numbers differ considerably between strains. Expression ratios of 16S-BA-alleles were proportional to the respective genomic allele copy numbers. The findings and experimental tools presented here provide detailed insights into the intra- and intergenomic diversity of 16S rRNA genes and may pave the way for improved identification of B. anthracis and other pathogens with diverse rRNA operons. IMPORTANCE For severe infectious diseases, precise pathogen detection is crucial for antibiotic therapy and patient survival. Identification of Bacillus anthracis, the causative agent of the zoonosis anthrax, can be challenging when querying specific nucleotide sequences such as in small subunit rRNA (16S rRNA) genes, which are commonly used for typing of bacteria. This study analyzed on a broad genomic scale a cryptic and hitherto underappreciated allelic variant of the bacterium's 16S rRNA genes and their transcripts using a set of in situ, in vitro, and in silico assays and found significant intra- and intergenomic heterogeneity in the distribution of the allele and overall rRNA operon copy numbers. This allelic variation was uniquely species specific, which enabled sensitive pathogen detection on both DNA and transcript levels. The methodology used here is likely also applicable to other pathogens that are otherwise difficult to discriminate from their less harmful relatives.

CoxBase: an Online Platform for Epidemiological Surveillance, Visualization, Analysis, and Typing of Coxiella burnetii Genomic Sequences.

Q (query) fever is an infectious zoonotic disease caused by the Gram-negative bacterium Coxiella burnetii. Although the disease has been studied for decades, it still represents a threat due to sporadic outbreaks across farms in Europe. The absence of a central platform for Coxiella typing data management is an important epidemiological gap that is relevant in the case of an outbreak. To fill this gap, we have designed and implemented an online, open-source, web-based platform called CoxBase (https://coxbase.q-gaps.de). This platform includes a database that holds genotyping information on more than 400 Coxiella isolates alongside metadata that annotate them. We have also implemented features for in silico genotyping of completely or minimally assembled Coxiella sequences using five different typing methods, querying of existing isolates, visualization of isolate geodata via aggregation on a world map, and submission of new isolates. We tested our in silico typing method on 50 Coxiella genomes downloaded from the RefSeq database, and we successfully genotyped all genomes except for cases where the sequence quality was poor. We identified new spacer sequences using our implementation of the multispacer sequence typing (MST) in silico typing method and established adaA gene phenotypes for all 50 genomes as well as their plasmid types. IMPORTANCE Q fever is a zoonotic disease that is a source of active epidemiological concern due to its persistent threat to public health. In this project, we have identified areas in the field of Coxiella research, especially regarding public health and genomic analysis, where there is an inadequacy of resources to monitor, organize, and analyze genomic data from C. burnetii. Subsequently, we have created an open, web-based platform that contains epidemiological information, genome typing functions comprising all the available Coxiella typing methods, and tools for isolate data discovery and visualization that could help address the above-mentioned challenges. This is the first platform to combine all disparate genotyping systems for Coxiella burnetii as well as metadata assets with tools for genomic comparison and analyses. This platform is a valuable resource for laboratory researchers as well as research epidemiologists interested in investigating the relatedness or dissimilarity among C. burnetii strains.

Ultrasensitive Detection of Bacillus anthracis by Real-Time PCR Targeting a Polymorphism in Multi-Copy 16S rRNA Genes and Their Transcripts.

The anthrax pathogen Bacillus anthracis poses a significant threat to human health. Identification of B. anthracis is challenging because of the bacterium's close genetic relationship to other Bacillus cereus group species. Thus, molecular detection is founded on species-specific PCR targeting single-copy genes. Here, we validated a previously recognized multi-copy target, a species-specific single nucleotide polymorphism (SNP) present in 2-5 copies in every B. anthracis genome analyzed. For this, a hydrolysis probe-based real-time PCR assay was developed and rigorously tested. The assay was specific as only B. anthracis DNA yielded positive results, was linear over 9 log10 units, and was sensitive with a limit of detection (LoD) of 2.9 copies/reaction. Though not exhibiting a lower LoD than established single-copy PCR targets (dhp61 or PL3), the higher copy number of the B. anthracis-specific 16S rRNA gene alleles afforded ≤2 unit lower threshold (Ct) values. To push the detection limit even further, the assay was adapted for reverse transcription PCR on 16S rRNA transcripts. This RT-PCR assay was also linear over 9 log10 units and was sensitive with an LoD of 6.3 copies/reaction. In a dilution series of experiments, the 16S RT-PCR assay achieved a thousand-fold higher sensitivity than the DNA-targeting assays. For molecular diagnostics, we recommend a real-time RT-PCR assay variant in which both DNA and RNA serve as templates (thus, no requirement for DNase treatment). This can at least provide results equaling the DNA-based implementation if no RNA is present but is superior even at the lowest residual rRNA concentrations.

FlexTaxD: flexible modification of taxonomy databases for improved sequence classification.

SUMMARY: The Flexible Taxonomy Database framework provides a method for modification and merging official and custom taxonomic databases to create improved databases. Using such databases will increase accuracy and precision of existing methods to classify sequence reads. AVAILABILITY AND IMPLEMENTATION: Source code is freely available at https://github.com/FOI-Bioinformatics/flextaxd and installable through Bioconda.

Catching SARS-CoV-2 by Sequence Hybridization: a Comparative Analysis.

Controlling and monitoring the still ongoing severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic regarding geographical distribution, evolution, and emergence of new mutations of the SARS-CoV-2 virus is only possible due to continuous next-generation sequencing (NGS) and sharing sequence data worldwide. Efficient sequencing strategies enable the retrieval of increasing numbers of high-quality, full-length genomes and are, hence, indispensable. Two opposed enrichment methods, tiling multiplex PCR and sequence hybridization by bait capture, have been established for SARS-CoV-2 sequencing and are both frequently used, depending on the quality of the patient sample and the question at hand. Here, we focused on the evaluation of the sequence hybridization method by studying five commercially available sequence capture bait panels with regard to sensitivity and capture efficiency. We discovered the SARS-CoV-2-specific panel of Twist Bioscience to be the most efficient panel, followed by two respiratory panels from Twist Bioscience and Illumina, respectively. Our results provide on the one hand a decision basis for the sequencing community including a computation for using the full capacity of the flow cell and on the other hand potential improvements for the manufacturers. IMPORTANCE Sequencing the genomes of the circulating SARS-CoV-2 strains is the only way to monitor the viral spread and evolution of the virus. Two different approaches, namely, tiling multiplex PCR and sequence hybridization by bait capture, are commonly used to fulfill this task. This study describes for the first time a combined approach of droplet digital PCR (ddPCR) and NGS to evaluate five commercially available sequence capture panels targeting SARS-CoV-2. In doing so, we were able to determine the most sensitive and efficient capture panel, distinguish the mode of action of the various bait panels, and compute the number of read pairs needed to recover a high-quality full-length genome. By calculating the minimum number of read pairs needed, we are providing optimized flow cell loading conditions for all sequencing laboratories worldwide that are striving for maximizing sequencing output and simultaneously minimizing time, costs, and sequencing resources.

First Phylogenetic Analysis of Malian SARS-CoV-2 Sequences Provides Molecular Insights into the Genomic Diversity of the Sahel Region.

We are currently facing a pandemic of COVID-19, caused by a spillover from an animal-originating coronavirus to humans occurring in the Wuhan region of China in December 2019. From China, the virus has spread to 188 countries and regions worldwide, reaching the Sahel region on March 2, 2020. Since whole genome sequencing (WGS) data is very crucial to understand the spreading dynamics of the ongoing pandemic, but only limited sequencing data is available from the Sahel region to date, we have focused our efforts on generating the first Malian sequencing data available. Screening 217 Malian patient samples for the presence of SARS-CoV-2 resulted in 38 positive isolates, from which 21 whole genome sequences were generated. Our analysis shows that both the early A (19B) and the later observed B (20A/C) clade are present in Mali, indicating multiple and independent introductions of SARS-CoV-2 to the Sahel region.

Whole genome sequencing and phylogenetic classification of Tunisian SARS-CoV-2 strains from patients of the Military Hospital in Tunis.

In the present work, two complete genome sequences of SARS-CoV-2 were obtained from nasal swab samples of Tunisian SARS-CoV-2 PCR-positive patients using nanopore sequencing. The virus genomes of two of the patients examined, a Tunisian soldier returning from a mission in Morocco and a member of another Tunisian family, showed significant differences in analyses of the total genome and single nucleotide polymorphisms (SNPs). Phylogenetic relationships with known SARS-CoV-2 genomes in the African region, some European and Middle Eastern countries and initial epidemiological conclusions indicate that the introduction of SARS-CoV-2 into Tunisia from two independent sources was travel-related.

A headache with surprising outcome: first case of brucellosis caused by Brucella suis biovar 1 in Germany.

In July 2018, brucellosis was diagnosed in a German patient without a travel history to regions endemic for Brucella. Microbiological analysis, including whole-genome sequencing, revealed Brucella suis biovar 1 as the etiologic agent. Core-genome-based multilocus sequence-typing analysis placed the isolate in close proximity to strains originating from Argentina. Notably, despite a strong IgM response, the patient did not develop Brucella-specific IgG antibodies during infection. Here, we describe the clinical course of infection, the extensive epidemiological investigations, and discuss possible routes of transmission.

Author Correction: Genome expansion and lineage-specific genetic innovations in the forest pathogenic fungi Armillaria.

In the version of this Article originally published, it was incorrectly stated that "16,687 protein-coding genes were inferred for the most recent common ancestor (MRCA) of Armillaria"; the value was incorrect and it should have read "15,787". This has now been corrected.

Genome expansion and lineage-specific genetic innovations in the forest pathogenic fungi Armillaria.

Armillaria species are both devastating forest pathogens and some of the largest terrestrial organisms on Earth. They forage for hosts and achieve immense colony sizes via rhizomorphs, root-like multicellular structures of clonal dispersal. Here, we sequenced and analysed the genomes of four Armillaria species and performed RNA sequencing and quantitative proteomic analysis on the invasive and reproductive developmental stages of A. ostoyae. Comparison with 22 related fungi revealed a significant genome expansion in Armillaria, affecting several pathogenicity-related genes, lignocellulose-degrading enzymes and lineage-specific genes expressed during rhizomorph development. Rhizomorphs express an evolutionarily young transcriptome that shares features with the transcriptomes of both fruiting bodies and vegetative mycelia. Several genes show concomitant upregulation in rhizomorphs and fruiting bodies and share cis-regulatory signatures in their promoters, providing genetic and regulatory insights into complex multicellularity in fungi. Our results suggest that the evolution of the unique dispersal and pathogenicity mechanisms of Armillaria might have drawn upon ancestral genetic toolkits for wood-decay, morphogenesis and complex multicellularity.

Complete Circular Genome Sequence and Temperature Independent Adaptation to Anaerobiosis of Listeria weihenstephanensis DSM 24698.

The aim of this study was to analyze the adaptation of the environmental Listeria weihenstephanensis DSM 24698 to anaerobiosis. The complete circular genome sequence of this species is reported and the adaptation of L. weihenstephanensis DSM 24698 to oxygen availability was investigated by global transcriptional analyses via RNAseq at 18 and 34°C. A list of operons was created based on the transcriptional data. Forty-two genes were upregulated anaerobically and 62 genes were downregulated anaerobically. The oxygen dependent gene expression of selected genes was further validated via qPCR. Many of the differentially regulated genes encode metabolic enzymes indicating broad metabolic adaptations with respect to oxygen availability. Genes showing the strongest oxygen-dependent adaption encoded nitrate (narGHJI) and nitrite (nirBD) reductases. Together with the observation that nitrate supported anaerobic growth, these data indicate that L. weihenstephanensis DSM 24698 performs anaerobic nitrate respiration. The wide overlap between the oxygen-dependent transcriptional regulation at 18 and 34°C suggest that temperature does not play a key role in the oxygen-dependent transcriptional regulation of L. weihenstephanensis DSM 24698.

Genome Sequence of Bacillus safensis Strain Ingolstadt Isolated from the Pectoralis Pouch of a Patient with Defibrillator-Related Surgery.

We report the draft genome sequence of clindamycin-resistant Bacillus safensis strain Ingolstadt isolated from a patient with bacterial colonization after heart surgery. The draft genome comprises 3.75 Mbp and harbors 3,793 predicted protein-encoding genes and a small plasmid.

Whole genome sequencing of Brucella melitensis isolated from 57 patients in Germany reveals high diversity in strains from Middle East.

Brucellosis, a worldwide common bacterial zoonotic disease, has become quite rare in Northern and Western Europe. However, since 2014 a significant increase of imported infections caused by Brucella (B.) melitensis has been noticed in Germany. Patients predominantly originated from Middle East including Turkey and Syria. These circumstances afforded an opportunity to gain insights into the population structure of Brucella strains. Brucella-isolates from 57 patients were recovered between January 2014 and June 2016 with culture confirmed brucellosis by the National Consultant Laboratory for Brucella. Their whole genome sequences were generated using the Illumina MiSeq platform. A whole genome-based SNP typing assay was developed in order to resolve geographically attributed genetic clusters. Results were compared to MLVA typing results, the current gold-standard of Brucella typing. In addition, sequences were examined for possible genetic variation within target regions of molecular diagnostic assays. Phylogenetic analyses revealed spatial clustering and distinguished strains from different patients in either case, whereas multiple isolates from a single patient or technical replicates showed identical SNP and MLVA profiles. By including WGS data from the NCBI database, five major genotypes were identified. Notably, strains originating from Turkey showed a high diversity and grouped into seven subclusters of genotype II. MLVA analysis congruently clustered all isolates and predominantly matched the East Mediterranean genetic clade. This study confirms whole-genome based SNP-analysis as a powerful tool for accurate typing of B. melitensis. Furthermore it allows special allocation and therefore provides useful information on the geographic origin for trace-back analysis. However, the lack of reliable metadata in public databases often prevents a resolution below geographic regions or country levels and corresponding precise trace-back analysis. Once this obstacle is resolved, WGS-derived bacterial typing adds an important method to complement epidemiological surveys during outbreak investigations. This is the first report of a detailed genetic investigation of an extensive collection of B. melitensis strains isolated from human cases in Germany.

MinION as part of a biomedical rapidly deployable laboratory.

Fast turnaround times are of utmost importance for biomedical reconnaissance, particularly regarding dangerous pathogens. Recent advances in sequencing technology and its devices allow sequencing within a short time frame outside stationary laboratories close to the epicenter of the outbreak. In our study, we evaluated the portable sequencing device MinION as part of a rapidly deployable laboratory specialized in identification of highly pathogenic agents. We tested the device in the course of a NATO live agent exercise in a deployable field laboratory in hot climate conditions. The samples were obtained from bio-terroristic scenarios that formed part of the exercise and contained unknown bacterial agents. To simulate conditions of a resource-limited remote deployment site, we operated the sequencer without internet access. Using a metagenomic approach, we were able to identify the causative agent in the analyzed samples. Furthermore, depending on the obtained data, we were able to perform molecular typing down to strain level. In our study we challenged the device and discuss advances as well as remaining limitations for sequencing biological samples outside of stationary laboratories. Nevertheless, massive parallel sequencing as a non-selective methodology yields important information and is able to support outbreak investigation - even in the field.

A complete toolset for the study of Ustilago bromivora and Brachypodium sp. as a fungal-temperate grass pathosystem.

Due to their economic relevance, the study of plant pathogen interactions is of importance. However, elucidating these interactions and their underlying molecular mechanisms remains challenging since both host and pathogen need to be fully genetically accessible organisms. Here we present milestones in the establishment of a new biotrophic model pathosystem: Ustilago bromivora and Brachypodium sp. We provide a complete toolset, including an annotated fungal genome and methods for genetic manipulation of the fungus and its host plant. This toolset will enable researchers to easily study biotrophic interactions at the molecular level on both the pathogen and the host side. Moreover, our research on the fungal life cycle revealed a mating type bias phenomenon. U. bromivora harbors a haplo-lethal allele that is linked to one mating type region. As a result, the identified mating type bias strongly promotes inbreeding, which we consider to be a potential speciation driver.

The Intervening Sequence of Coxiella burnetii: Characterization and Evolution.

The intervening sequence (IVS) of Coxiella burnetii, the agent of Q fever, is a 428-nt selfish genetic element located in helix 45 of the precursor 23S rRNA. The IVS element, in turn, contains an ORF that encodes a hypothetical ribosomal S23 protein (S23p). Although S23p can be synthesized in vitro in the presence of an engineered E. coli promoter and ribosome binding site, results suggest that the protein is not synthesized in vivo. In spite of a high degree of IVS conservation among different strains of C. burnetii, the region immediately upstream of the S23p start codon is prone to change, and the S23p-encoding ORF is evidently undergoing reductive evolution. We determined that IVS excision from 23S rRNA was mediated by RNase III, and IVS RNA was rapidly degraded, thereafter. Levels of the resulting 23S rRNA fragments that flank the IVS, F1 (~1.2 kb) and F2 (~1.7 kb), were quantified over C. burnetii's logarithmic growth phase (1-5 d). Results showed that 23S F1 quantities were consistently higher than those of F2 and 16S rRNA. The disparity between levels of the two 23S rRNA fragments following excision of IVS is an interesting phenomenon of unknown significance. Based upon phylogenetic analyses, IVS was acquired through horizontal transfer after C. burnetii's divergence from an ancestral bacterium and has been subsequently maintained by vertical transfer. The widespread occurrence, maintenance and conservation of the IVS in C. burnetii imply that it plays an adaptive role or has a neutral effect on fitness.

Unexpected genomic relationships between Bacillus anthracis strains from Bangladesh and Central Europe.

The zoonosis anthrax caused by the bacterium Bacillus anthracis has a broad geographical distribution. Active enzootic areas are typically located away from central and northern Europe where cases of the disease occur only sporadically and in limited numbers. In contrast, a few out of the 64 districts of Bangladesh are hyper-endemic for anthrax and there the disease causes major losses in live-stock. In this study we genotyped eight strains of B. anthracis collected from the districts of Sirajganj and Tangail in 2013. All these strains belonged to canSNP group A.Br.001/002 Sterne differing only in a few of 31 tandem-repeat (MLVA)-markers. Whole genome sequences were obtained from five of these strains and compared with genomic information of B. anthracis strains originating from various geographical locations. Characteristic signatures were detected defining two "Bangladesh" clusters potentially useful for rapid molecular epidemiology. From this data high-resolution PCR assays were developed and subsequently tested on additional isolates from Bangladesh and Central Europe. Remarkably, this comparative genomic analysis focusing on SNP-discovery revealed a close genetic relationship between these strains from Bangladesh and historic strains collected between 1991 and 2008 in The Netherlands and Germany, respectively. Possible explanations for these phylogenetic relationships are discussed.

Coxiella-like endosymbiont in argasid ticks (Ornithodoros muesebecki) from a Socotra Cormorant colony in Umm Al Quwain, United Arab Emirates.

Coxiella burnetii is a pathogen causing Q fever in domestic animals and humans. Seabirds have been implicated as possible reservoirs of this bacterium in the Arabian Gulf and in the Western Indian Ocean. Recently, Coxiella species closely related to C. burnetii was detected from ticks collected from oil rigs used as roosting areas by Socotra Cormorants (Phalacrocorax nigrogularis) in the western Arabian Gulf. We collected ticks from the largest breeding colony of Socotra Cormorants in the United Arab Emirates on the eastern extreme of the species' breeding range to determine the prevalence of C. burnetii and evaluate its role as a wild reservoir. All ticks were identified as Ornithodoros muesebecki and genomic DNA was extracted from larval and nymph/adult tick pools. Multiplex PCR tests were performed targeting three C. burnetii specific genes. C. burnetii was not detected although a Coxiella-like endosymbiont was identified that was closely related to Coxiella symbionts from Ornithodoros capensis ticks. Because domestic and wild ungulates are the primary source of C. burnetii, we suggest that the presence of free-ranging, native and non-native ungulates in some off-shore islands in the Arabian Gulf could disseminate C. burnetii to seabirds. More comprehensive studies on seabird colonies are needed to better understand the diversity and prevalence of Coxiella symbionts and to establish if C. burnetii is endemic on some of these islands.

eggNOG 4.5: a hierarchical orthology framework with improved functional annotations for eukaryotic, prokaryotic and viral sequences.

eggNOG is a public resource that provides Orthologous Groups (OGs) of proteins at different taxonomic levels, each with integrated and summarized functional annotations. Developments since the latest public release include changes to the algorithm for creating OGs across taxonomic levels, making nested groups hierarchically consistent. This allows for a better propagation of functional terms across nested OGs and led to the novel annotation of 95 890 previously uncharacterized OGs, increasing overall annotation coverage from 67% to 72%. The functional annotations of OGs have been expanded to also provide Gene Ontology terms, KEGG pathways and SMART/Pfam domains for each group. Moreover, eggNOG now provides pairwise orthology relationships within OGs based on analysis of phylogenetic trees. We have also incorporated a framework for quickly mapping novel sequences to OGs based on precomputed HMM profiles. Finally, eggNOG version 4.5 incorporates a novel data set spanning 2605 viral OGs, covering 5228 proteins from 352 viral proteomes. All data are accessible for bulk downloading, as a web-service, and through a completely redesigned web interface. The new access points provide faster searches and a number of new browsing and visualization capabilities, facilitating the needs of both experts and less experienced users. eggNOG v4.5 is available at http://eggnog.embl.de.