In-depth analysis of Bacillus anthracis 16S rRNA genes and transcripts reveals intra- and intergenomic diversity and facilitates anthrax detection

Analysis of 16S ribosomal RNA (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 yet 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 genome 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.