Background: The genome of Bacillus anthracis, the etiological agent of anthrax, is highly monomorphic which makes differentiation between strains difficult. A Multiple Locus Variable-number tandem repeats (VNTR) Analysis (MLVA) assay based on 20 markers was previously described. It has considerable discrimination power, reproducibility, and low cost, especially since the markers proposed can be typed by agarose-gel electrophoresis. However in an emergency situation, faster genotyping and access to representative databases is necessary.

Anthrax is a zoonotic disease that occurs naturally in wild and domestic animals but has been used by both state-sponsored programs and terrorists as a biological weapon. A Soviet industrial production facility in Sverdlovsk, USSR, proved deficient in 1979 when a plume of spores was accidentally released and resulted in one of the largest known human anthrax outbreaks. In order to understand this outbreak and others, we generated a Bacillus anthracis population genetic database based upon whole-genome analysis to identify all single-nucleotide polymorphisms (SNPs) across a reference genome. Phylogenetic analysis has defined three major clades (A, B, and C), B and C being relatively rare compared to A. The A clade has numerous subclades, including a major polytomy named the trans-Eurasian (TEA) group. The TEA radiation is a dominant evolutionary feature of B. anthracis, with many contemporary populations having resulted from a large spatial dispersal of spores from a single source. Two autopsy specimens from the Sverdlovsk outbreak were deep sequenced to produce draft B. anthracis genomes. This allowed the phylogenetic placement of the Sverdlovsk strain into a clade with two Asian live vaccine strains, including the Russian Tsiankovskii strain. The genome was examined for evidence of drug resistance manipulation or other genetic engineering, but none was found. The Soviet Sverdlovsk strain genome is consistent with a wild-type strain from Russia that had no evidence of genetic manipulation during its industrial production. This work provides insights into the world’s largest biological weapons program and provides an extensive B. anthracis phylogenetic reference.

1314

New SARS-CoV-2 mutations are constantly emerging, raising concerns of increased transmissibility, virulence or escape from host immune response. We describe a nested RT-PCR assay (~1500 bps) to detect multiple nucleotide changes resulting in key spike protein mutations distinctive of the major known circulating SARS-CoV-2 variants, including the three Variants of Concern (VOCs) 20I/501Y.V1 (United Kingdom), 20H/501Y.V2 (South Africa), and 20 J/501Y.V3 (Brazil), as well as the 20E.EU1 variant (Spain), the CAL.20C recently identified in California, and the mink-associated variant (GR, lineage B.1.1.298). Prior to application to field samples, the discriminatory potential of this PCR assay was explored using GISAID and Nextclade. To extend variant detection to challenging matrices such as sewage, where the amplification of long fragments is problematic, two short nested RT-PCR assays (~300 bps) were also designed, targeting portions of the region spanned by the long nested assay. The three newly-designed assays were then tested on field samples, including 31 clinical samples (7 fully-sequenced swab samples, and 24 uncharacterized ones) and 34 urban wastewater samples, some of which collected in areas where circulation of VOCs had been reported. The long assay successfully amplified 29 of the 31 swabs (93%), allowing the correct identification of variants 20I/501Y.V1 and 20E.EU1 present in the panel of previously characterized samples. The Spanish variant was detected in 14/24 of the uncharacterized samples as well. The sequences obtained using the short assays were consistent with those obtained with the long assay. Mutations characteristic of VOCs (UK and Brazilian variant) and of other variant (Spanish) were detected in sewage samples. To our knowledge, this is the first evidence of the presence of sequences harboring key mutations of 20I/501Y.V1 and 20 J/501Y.V3 in urban wastewaters, highlighting the potential contribution of wastewater surveillance to explore SARS-CoV-2 diversity. The developed nested RT-PCR assays can be used as an initial rapid screening test to select clinical samples containing mutations of interest. This can speed up diagnosis and optimize resources since it allows full genome sequencing to be done only on clinically relevant specimens. The assays can be also employed for a rapid and cost-effective detection of VOCs or other variants in sewage for the purposes of wastewater-based epidemiology. The approach proposed here can be used to better understand SARS-CoV-2 variant diversity, geographic distribution and impact worldwide.

We used multiple-locus variable-number tandem repeat analysis (MLVA) to type 64 Bacillus anthracis isolates from outbreaks that have occurred during the past 40 years in Italy. MLVA of the 64 isolates revealed 10 unique genotypes; 9 of these genotypes and the majority of isolates (63/64) belonged to the previously described genetic cluster A1.a. Within the A1.a isolates, two previously described genotypes (G1 and G3), which differ by a single mutation in the pX01 locus, account for the majority of isolates in the country (53/63). The low diversity of B. anthracis genotypes in Italy suggests a single, dominant historical introduction, followed by limited localized differentiation.