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.

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.

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.

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.

HoPaCI-DB: host-Pseudomonas and Coxiella interaction database.

Bacterial infectious diseases are the result of multifactorial processes affected by the interplay between virulence factors and host targets. The host-Pseudomonas and Coxiella interaction database (HoPaCI-DB) is a publicly available manually curated integrative database (http://mips.helmholtz-muenchen.de/HoPaCI/) of host-pathogen interaction data from Pseudomonas aeruginosa and Coxiella burnetii. The resource provides structured information on 3585 experimentally validated interactions between molecules, bioprocesses and cellular structures extracted from the scientific literature. Systematic annotation and interactive graphical representation of disease networks make HoPaCI-DB a versatile knowledge base for biologists and network biology approaches.

Molecular analysis of Coxiella burnetii in Germany reveals evolution of unique clonal clusters.

The causative agent of Q fever, Coxiella burnetii, is a query agent occurring naturally all over the world. We studied 104 German Coxiella burnetii strains/DNA samples obtained between 1969 and 2011 using a 14 microsatellite marker Multiple-locus variable-number of tandem repeat (VNTR) analysis (MLVA) technique. We were able to divide our collection into 32 different genotypes clustered into four major groups (A-D). Two of these (A and C) formed predominant clonal complexes that covered 97% of all studied samples. Group C consisted exclusively of cattle-associated isolates/DNA specimens, while group A comprised all other affected species including all sheep-derived strains/DNA samples. Within this second cluster, two major genotypes (A1, A2) were identified. Genotype A2 occurred in strains isolated from ewes in northern and central Germany, whereas genotype A1 was found in most areas of Germany. MLVA analysis of C. burnetii strains from neighbouring countries revealed a close relationship to German strains. We thus hypothesize that there is a western and central European cluster of C. burnetii. We identified predominant genotypes related to relevant host species and geographic regions which is in line with findings of the Dutch Q fever outbreak (2007-2010). Furthermore three of our analyzed German strains are closely related to the Dutch outbreak clone. These findings support the theory of predominant genotypes in the context of regional outbreaks. Our results show that a combination of 8 MLVA markers provides the highest discriminatory power for attributing C. burnetii isolates to genotypes. For future epidemiological studies we propose the use of three MLVA markers for easy and rapid classification of C. burnetii into 4 main clusters.

Sequence-structure relationships in yeast mRNAs.

It is generally accepted that functionally important RNA structure is more conserved than sequence due to compensatory mutations that may alter the sequence without disrupting the structure. For small RNA molecules sequence-structure relationships are relatively well understood. However, structural bioinformatics of mRNAs is still in its infancy due to a virtual absence of experimental data. This report presents the first quantitative assessment of sequence-structure divergence in the coding regions of mRNA molecules based on recently published transcriptome-wide experimental determination of their base paring patterns. Structural resemblance in paralogous mRNA pairs quickly drops as sequence identity decreases from 100% to 85-90%. Structures of mRNAs sharing sequence identity below roughly 85% are essentially uncorrelated. This outcome is in dramatic contrast to small functional non-coding RNAs where sequence and structure divergence are correlated at very low levels of sequence similarity. The fact that very similar mRNA sequences can have vastly different secondary structures may imply that the particular global shape of base paired elements in coding regions does not play a major role in modulating gene expression and translation efficiency. Apparently, the need to maintain stable three-dimensional structures of encoded proteins places a much higher evolutionary pressure on mRNA sequences than on their RNA structures.

Control of Coxiella burnetii shedding in a dairy goat herd by annual offspring vaccination.

A Coxiella burnetii vaccination program, targeting only doelings, was introduced on a German goat farm to curb bacterial shedding. In 2018, adults were vaccinated with a C. burnetii Phase I vaccine at three-weeks apart following pathogen diagnosis, with a booster administered six months later due to sustained high shedding. From 2018 to 2021, doelings received two vaccine doses without any further boosters. To assess the program's efficacy, vaginal swabs from up to 40 animals per age group were collected during kidding seasons from 2019 to 2022. Bulk tank milk (BTM) samples were gathered monthly from January 2018 to October 2022 to monitor herd-level shedding. Real-time PCR analysis determined genome equivalents in all three sample types. Serum samples were taken before the initial immunization and during the post-kidding season from up to 40 goats per age group annually from 2018 to 2022. Phase-specific ELISAs determined IgG Phase I and Phase II antibodies. Additionally, two serum samples per age group from 2022 were analyzed using a neutralization assay. A few goats continued shedding small quantities during subsequent kidding seasons. Although positive BTM samples decreased, they displayed an undulating trend. Most age groups exhibited robust IgG Phase I responses and lower IgG Phase II levels post immunization. Mean IgG levels remained elevated until the study ended compared to pre-vaccination levels in most age groups. Additionally, neutralizing antibodies were present regardless of IgG response. Overall, double vaccination induced lasting antibody levels, but did not entirely prevent C. burnetii shedding. The resilience of the observed humoral immune activity requires further investigation.

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.

PEDANT covers all complete RefSeq genomes.

The PEDANT genome database provides exhaustive annotation of nearly 3000 publicly available eukaryotic, eubacterial, archaeal and viral genomes with more than 4.5 million proteins by a broad set of bioinformatics algorithms. In particular, all completely sequenced genomes from the NCBI's Reference Sequence collection (RefSeq) are covered. The PEDANT processing pipeline has been sped up by an order of magnitude through the utilization of precalculated similarity information stored in the similarity matrix of proteins (SIMAP) database, making it possible to process newly sequenced genomes immediately as they become available. PEDANT is freely accessible to academic users at http://pedant.gsf.de. For programmatic access Web Services are available at http://pedant.gsf.de/webservices.jsp.

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.

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.

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.

Beyond the 'best' match: machine learning annotation of protein sequences by integration of different sources of information.

MOTIVATION: Accurate automatic assignment of protein functions remains a challenge for genome annotation. We have developed and compared the automatic annotation of four bacterial genomes employing a 5-fold cross-validation procedure and several machine learning methods. RESULTS: The analyzed genomes were manually annotated with FunCat categories in MIPS providing a gold standard. Features describing a pair of sequences rather than each sequence alone were used. The descriptors were derived from sequence alignment scores, InterPro domains, synteny information, sequence length and calculated protein properties. Following training we scored all pairs from the validation sets, selected a pair with the highest predicted score and annotated the target protein with functional categories of the prototype protein. The data integration using machine-learning methods provided significantly higher annotation accuracy compared to the use of individual descriptors alone. The neural network approach showed the best performance. The descriptors derived from the InterPro domains and sequence similarity provided the highest contribution to the method performance. The predicted annotation scores allow differentiation of reliable versus non-reliable annotations. The developed approach was applied to annotate the protein sequences from 180 complete bacterial genomes. AVAILABILITY: The FUNcat Annotation Tool (FUNAT) is available on-line as Web Services at http://mips.gsf.de/proj/funat.

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.

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.

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.