CYP8B1 downregulation mediates the metabolic effects of vertical sleeve gastrectomy in mice.

BACKGROUND AND AIMS: Although the benefits of vertical sleeve gastrectomy (VSG) surgery are well known, the molecular mechanisms by which VSG alleviates obesity and its complications remain unclear. We aim to determine the role of CYP8B1 (cytochrome P450, family 8, subfamily B, polypeptide 1) in mediating the metabolic benefits of VSG. APPROACH AND RESULTS: We found that expression of CYP8B1, a key enzyme in controlling the 12α-hydroxylated (12α-OH) bile acid (BA) to non-12α-OH BA ratio, was strongly downregulated after VSG. Using genetic mouse models of CYP8B1 overexpression, knockdown, and knockout, we demonstrated that overexpression of CYP8B1 dampened the metabolic improvements associated with VSG. In contrast, short hairpin RNA-mediated CYP8B1 knockdown improved metabolism similar to those observed after VSG. Cyp8b1 deficiency diminished the metabolic effects of VSG. Further, VSG-induced alterations to the 12α-OH/non-12α-OH BA ratio in the BA pool depended on CYP8B1 expression level. Consequently, intestinal lipid absorption was restricted, and the gut microbiota (GM) profile was altered. Fecal microbiota transplantation from wild type-VSG mice (vs. fecal microbiota transplantation from wild-type-sham mice) improved metabolism in recipient mice, while there were no differences between mice that received fecal microbiota transplantation from knockout-sham and knockout-VSG mice. CONCLUSIONS: CYP8B1 is a critical downstream target of VSG. Modulation of BA composition and gut microbiota profile by targeting CYP8B1 may provide novel insight into the development of therapies that noninvasively mimic bariatric surgery to treat obesity and its complications.

A Phase II Clinical Trial of Pembrolizumab and Enobosarm in Patients with Androgen Receptor-Positive Metastatic Triple-Negative Breast Cancer.

LESSONS LEARNED: The combination of enobosarm and pembrolizumab was well tolerated and showed a modest clinical benefit rate of 25% at 16 weeks. Future trials investigating androgen receptor-targeted therapy in combination with immune checkpoint inhibitors are warranted. BACKGROUND: Luminal androgen receptor is a distinct molecular subtype of triple-negative breast cancer (TNBC) defined by overexpression of androgen receptor (AR). AR-targeted therapy has shown modest activity in AR-positive (AR+) TNBC. Enobosarm (GTx-024) is a nonsteroidal selective androgen receptor modulator (SARM) that demonstrates preclinical and clinical activity in AR+ breast cancer. The current study was designed to explore the safety and efficacy of the combination of enobosarm and pembrolizumab in patients with AR+ metastatic TNBC (mTNBC). METHODS: This study was an open-label phase II study for AR+ (≥10%, 1+ by immunohistochemistry [IHC]) mTNBC. Eligible patients received pembrolizumab 200 mg intravenous (IV) every 3 weeks and enobosarm 18 mg oral daily. The primary objective was to evaluate the safety of enobosarm plus pembrolizumab and determine the response rate. Peripheral blood, tumor biopsies, and stool samples were collected for correlative analysis. RESULTS: The trial was stopped early because of the withdrawal of GTx-024 drug supply. Eighteen patients were enrolled, and 16 were evaluable for responses. Median age was 64 (range 36-81) years. The combination was well tolerated, with only a few grade 3 adverse events: one dry skin, one diarrhea, and one musculoskeletal ache. The responses were 1 of 16 (6%) complete response (CR), 1 of 16 (6%) partial response (PR), 2 of 16 (13%) stable disease (SD), and 12 of 16 (75%) progressive disease (PD). Response rate (RR) was 2 of 16 (13%). Clinical benefit rate (CBR) at 16 weeks was 4 of 16 (25%). Median follow-up was 24.9 months (95% confidence interval [CI], 17.5-30.9). Progression-free survival (PFS) was 2.6 months (95% CI, 1.9-3.1) and overall survival (OS) was 25.5 months (95% CI, 10.4-not reached [NR]). CONCLUSION: The combination of enobosarm and pembrolizumab was well tolerated, with a modest clinical benefit rate of 25% at 16 weeks in heavily pretreated AR+ TNBC without preselected programmed death ligand-1 (PD-L1). Future clinical trials combining AR-targeted therapy with immune checkpoint inhibitor (ICI) for AR+ TNBC warrant investigation.

Rapid Drug Susceptibility Testing of Drug-Resistant Mycobacterium tuberculosis Isolates Directly from Clinical Samples by Use of Amplicon Sequencing: a Proof-of-Concept Study.

Increasingly complex drug-resistant tuberculosis (DR-TB) is a major global health concern and one of the primary reasons why TB is now the leading infectious cause of death worldwide. Rapid characterization of a DR-TB patient's complete drug resistance profile would facilitate individualized treatment in place of empirical treatment, improve treatment outcomes, prevent amplification of resistance, and reduce the transmission of DR-TB. The use of targeted next-generation sequencing (NGS) to obtain drug resistance profiles directly from patient sputum samples has the potential to enable comprehensive evidence-based treatment plans to be implemented quickly, rather than in weeks to months, which is currently needed for phenotypic drug susceptibility testing (DST) results. In this pilot study, we evaluated the performance of amplicon sequencing of Mycobacterium tuberculosis DNA from patient sputum samples using a tabletop NGS technology and automated data analysis to provide a rapid DST solution (the Next Gen-RDST assay). One hundred sixty-six out of 176 (94.3%) sputum samples from the Republic of Moldova yielded complete Next Gen-RDST assay profiles for 7 drugs of interest. We found a high level of concordance of our Next Gen-RDST assay results with phenotypic DST (97.0%) and pyrosequencing (97.8%) results from the same clinical samples. Our Next Gen-RDST assay was also able to estimate the proportion of resistant-to-wild-type alleles down to mixtures of ≤1%, which demonstrates the ability to detect very low levels of resistant variants not detected by pyrosequencing and possibly below the threshold for phenotypic growth methods. The assay as described here could be used as a clinical or surveillance tool.

Francisella tularensis subsp. tularensis group A.I, United States.

We used whole-genome analysis and subsequent characterization of geographically diverse strains using new genetic signatures to identify distinct subgroups within Francisella tularensis subsp. tularensis group A.I: A.I.3, A.I.8, and A.I.12. These subgroups exhibit complex phylogeographic patterns within North America. The widest distribution was observed for A.I.12, which suggests an adaptive advantage.

Whole-genome analysis of Exserohilum rostratum from an outbreak of fungal meningitis and other infections.

Exserohilum rostratum was the cause of most cases of fungal meningitis and other infections associated with the injection of contaminated methylprednisolone acetate produced by the New England Compounding Center (NECC). Until this outbreak, very few human cases of Exserohilum infection had been reported, and very little was known about this dematiaceous fungus, which usually infects plants. Here, we report using whole-genome sequencing (WGS) for the detection of single nucleotide polymorphisms (SNPs) and phylogenetic analysis to investigate the molecular origin of the outbreak using 22 isolates of E. rostratum retrieved from 19 case patients with meningitis or epidural/spinal abscesses, 6 isolates from contaminated NECC vials, and 7 isolates unrelated to the outbreak. Our analysis indicates that all 28 isolates associated with the outbreak had nearly identical genomes of 33.8 Mb. A total of 8 SNPs were detected among the outbreak genomes, with no more than 2 SNPs separating any 2 of the 28 genomes. The outbreak genomes were separated from the next most closely related control strain by ∼136,000 SNPs. We also observed significant genomic variability among strains unrelated to the outbreak, which may suggest the possibility of cryptic speciation in E. rostratum.

Real-time PCR assays for genotyping of Cryptococcus gattii in North America.

BACKGROUND: Cryptococcus gattii has been the cause of an ongoing outbreak starting in 1999 on Vancouver Island, British Columbia and spreading to mainland Canada and the US Pacific Northwest. In the course of the outbreak, C. gattii has been identified outside of its previously documented climate, habitat, and host disease. Genotyping of C. gattii is essential to understand the ecological and geographical expansion of this emerging pathogen. METHODS: We developed and validated a mismatch amplification mutation assay (MAMA) real-time PCR panel for genotyping C. gattii molecular types VGI-VGIV and VGII subtypes a,b,c. Subtype assays were designed based on whole-genome sequence of 20 C. gattii strains. Publically available multilocus sequence typing (MLST) data from a study of 202 strains was used for the molecular type (VGI-VGIV) assay design. All assays were validated across DNA from 112 strains of diverse international origin and sample types, including animal, environmental and human. RESULTS: Validation revealed each assay on the panel is 100% sensitive, specific and concordant with MLST. The assay panel can detect down to 0.5 picograms of template DNA. CONCLUSIONS: The (MAMA) real-time PCR panel for C. gattii accurately typed a collection of 112 diverse strains and demonstrated high sensitivity. This is a time and cost efficient method of genotyping C. gattii best suited for application in large-scale epidemiological studies.

The Role of Geography, Diet, and Host Phylogeny on the Gut Microbiome in the Hawaiian Honeycreeper Radiation.

The animal gut microbiome can have a strong influence on the health, fitness, and behavior of its hosts. The composition of the gut microbial community can be influenced by factors such as diet, environment, and evolutionary history (phylosymbiosis). However, the relative influence of these factors is unknown in most bird species. Furthermore, phylosymbiosis studies have largely focused on clades that diverged tens of millions of years ago, and little is known about the degree of gut microbiome divergence in more recent species radiations. This study explores the drivers of microbiome variation across the unique and recent Hawaiian honeycreeper radiation (Fringillidae: Drepanidinae). Fecal samples were collected from 14 extant species spanning the main islands of the Hawaiian archipelago and were sequenced using three metabarcoding markers to characterize the gut microbiome, invertebrate diet, and plant diet of Hawaiian honeycreepers. We then used these metabarcoding data and the honeycreeper host phylogeny to evaluate their relative roles in shaping the gut microbiome. Microbiome variation across birds was highly individualized; however, source island had a small but significant effect on microbiome structure. The microbiomes did not recapitulate the host phylogenetic tree, indicating that evolutionary history does not strongly influence microbiome structure in the honeycreeper clade. These results expand our understanding of the roles of diet, geography, and phylogeny on avian microbiome structure, while also providing important ecological information about the diet and gut microbiota of wild Hawaiian honeycreepers.

Multi-faceted metagenomic analysis of spacecraft associated surfaces reveal planetary protection relevant microbial composition.

The National Aeronautics and Space Administration (NASA) has been monitoring the microbial burden of spacecraft since the 1970's Viking missions. Originally culture-based and then focused 16S sequencing techniques were used, but we have now applied whole metagenomic sequencing to a variety of cleanroom samples at the Jet Propulsion Lab (JPL), including the Spacecraft Assembly Facility (SAF) with the goals of taxonomic identification and for functional assignment. Our samples included facility pre-filters, cleanroom vacuum debris, and surface wipes. The taxonomic composition was carried out by three different analysis tools to contrast marker, k-mer, and true alignment approaches. Hierarchical clustering analysis of the data separated vacuum particles from other SAF DNA samples. Vacuum particle samples were the most diverse while DNA samples from the ISO (International Standards Organization) compliant facilities and the SAF were the least diverse; all three were dominated by Proteobacteria. Wipe samples had higher diversity and were predominated by Actinobacteria, including human commensals Cutibacterium acnes and Corynebacterium spp. Taxa identified by the three methods were not identical, supporting the use of multiple methods for metagenome characterization. Likewise, functional annotation was performed using multiple methods. Vacuum particles and SAF samples contained strong signals of the tricarboxylic acid cycle and of amino acid biosynthesis, suggesting that many of the identified microorganisms have the ability to grow in nutrient-limited environments. In total, 18 samples generated high quality metagenome assembled genomes (MAG), which were dominated by Moraxella osloensis or Malassezia restricta. One M. osloensis MAG was assembled into a single circular scaffold and gene annotated. This study includes a rigorous quantitative determination of microbial loads and a qualitative dissection of microbial composition. Assembly of multiple specimens led to greater confidence for the identification of particular species and their predicted functional roles.

Comparison of Brucella abortus population structure based on genotyping methods with different levels of resolution.

Numerous genotyping techniques based on different principles and with different costs and levels of resolution are currently available for understanding the transmission dynamics of brucellosis worldwide. We aimed to compare the population structure of the genomes of 53 Brazilian Brucella abortus isolates using eight different genotyping methods: multiple-locus variable-number tandem-repeat analysis (MLVA8, MLVA11, MLVA16), multilocus sequence typing (MLST9, MLST21), core genome MLST (cgMLST) and two techniques based on single nucleotide polymorphism (SNP) detection (parSNP and NASP) from whole genomes. The strains were isolated from six different Brazilian states between 1977 and 2008 and had previously been analyzed using MLVA8, MLVA11, and MLVA16. Their whole genomes were sequenced, assembled, and subjected to MLST9 MLST21, cgMLST, and SNP analyses. All the genotypes were compared by hierarchical grouping method based on the average distances between the correlation matrices of each technique. MLST9 and MLST21 had the lowest level of resolution, both revealing only four genotypes. MLVA8, MLVA11, and MLVA16 had progressively increasing levels of resolution as more loci were analyzed, identifying 6, 16, and 44 genotypes, respectively. cgMLST showed the highest level of resolution, identifying 45 genotypes, followed by the SNP-based methods, both of which had 44 genotypes. In the assessed population, MLVA was more discriminatory than MLST and was easier and cheaper to perform. SNP techniques and cgMLST provided the highest levels of resolution and the results from the two methods were in close agreement. In conclusion, the choice of genotyping technique can strongly affect one's ability to make meaningful epidemiological conclusions but is dependent on available resources: while the VNTR based techniques are more indicated to high prevalence scenarios, the WGS methods are the ones with the best discriminative power and therefore recommended for outbreaks investigation.

Genome sequence of Staphylococcus aureus strain 11819-97, an ST80-IV European community-acquired methicillin-resistant isolate.

The European methicillin-resistant Staphylococcus aureus (MRSA) clone ST80-IV has historically dominated community-associated infections in major parts of Europe and is a lineage strongly linked to skin and soft tissue infections. Here, we report the genome sequence of an ST80-IV representative, 11819-97, isolated from a skin infection in Denmark in 1997.

Genome sequence of "Candidatus Microthrix parvicella" Bio17-1, a long-chain-fatty-acid-accumulating filamentous actinobacterium from a biological wastewater treatment plant.

"Candidatus Microthrix" bacteria are deeply branching filamentous actinobacteria which occur at the water-air interface of biological wastewater treatment plants, where they are often responsible for foaming and bulking. Here, we report the first draft genome sequence of a strain from this genus: "Candidatus Microthrix parvicella" strain Bio17-1.

Randomized trial assessing impact of probiotic supplementation on gut microbiome and clinical outcome from targeted therapy in metastatic renal cell carcinoma.

Studies suggest a link between the gut microbiome and metastatic renal cell carcinoma (mRCC) outcomes, including evidence that mRCC patients possess a lower abundance of Bifidobacterium spp. compared to healthy adults. We sought to assess if a Bifidobacterium-containing yogurt product could modulate the gut microbiome and clinical outcome from vascular endothelial growth factor-tyrosine kinase inhibitors (VEGF-TKIs). mRCC patients initiating VEGF-TKIs, regardless of the line of therapy, were randomized to probiotic-supplemented (two 4 oz. servings of the probiotic yogurt product daily) or probiotic-restricted arms. Stool samples were collected prior to therapy and at weeks 2, 3, 4, and 12. Microbiome composition was assessed using whole-metagenome sequencing. A total of 20 patients were randomized. Bifidobacterium animalis, the active ingredient of the probiotic supplement, reached detectable levels in all patients in the probiotic-supplemented arm versus two patients in the probiotic-restricted arm. Clinical benefit rate was similar in probiotic-supplemented versus probiotic-restricted arms (70% vs. 80%, p = 0.606). Linear discriminant analysis (LDA) effect size analysis of MetaPhIAn2 abundance data predicted 25 enriched species demonstrating an LDA score >3 in either clinical benefit or no clinical benefit. In patients with clinical benefit (vs. no clinical benefit), Barnesiella intestinihominis and Akkermansia muciniphila were significantly more abundant (p = 7.4 × 10-6 and p = 5.6 × 10-3 , respectively). This is the first prospective randomized study demonstrating modulation of the gut microbiome with a probiotic in mRCC. Probiotic supplementation successfully increased the Bifidobacterium spp. levels. Analysis of longitudinal stool specimens identified an association between B. intestinihominis, A. muciniphila, and clinical benefit with therapy. Trial Registration: NCT02944617.

Analysis of Gut Microbiome Using Explainable Machine Learning Predicts Risk of Diarrhea Associated With Tyrosine Kinase Inhibitor Neratinib: A Pilot Study.

Neratinib has great efficacy in treating HER2+ breast cancer but is associated with significant gastrointestinal toxicity. The objective of this pilot study was to understand the association of gut microbiome and neratinib-induced diarrhea. Twenty-five patients (age ≥ 60) were enrolled in a phase II trial evaluating safety and tolerability of neratinib in older adults with HER2+ breast cancer (NCT02673398). Fifty stool samples were collected from 11 patients at baseline and during treatment. 16S rRNA analysis was performed and relative abundance data were generated. Shannon's diversity was calculated to examine gut microbiome dysbiosis. An explainable tree-based approach was utilized to classify patients who might experience neratinib-related diarrhea (grade ≥ 1) based on pre-treatment baseline microbial relative abundance data. The hold-out Area Under Receiver Operating Characteristic and Area Under Precision-Recall Curves of the model were 0.88 and 0.95, respectively. Model explanations showed that patients with a larger relative abundance of Ruminiclostridium 9 and Bacteroides sp. HPS0048 may have reduced risk of neratinib-related diarrhea and was confirmed by Kruskal-Wallis test (p ≤ 0.05, uncorrected). Our machine learning model identified microbiota associated with reduced risk of neratinib-induced diarrhea and the result from this pilot study will be further verified in a larger study. CLINICAL TRIAL REGISTRATION: ClinicalTrials.gov, identifier NCT02673398.

Roles of bacteriophages, plasmids and CRISPR immunity in microbial community dynamics revealed using time-series integrated meta-omics.

Viruses and plasmids (invasive mobile genetic elements (iMGEs)) have important roles in shaping microbial communities, but their dynamic interactions with CRISPR-based immunity remain unresolved. We analysed generation-resolved iMGE-host dynamics spanning one and a half years in a microbial consortium from a biological wastewater treatment plant using integrated meta-omics. We identified 31 bacterial metagenome-assembled genomes encoding complete CRISPR-Cas systems and their corresponding iMGEs. CRISPR-targeted plasmids outnumbered their bacteriophage counterparts by at least fivefold, highlighting the importance of CRISPR-mediated defence against plasmids. Linear modelling of our time-series data revealed that the variation in plasmid abundance over time explained more of the observed community dynamics than phages. Community-scale CRISPR-based plasmid-host and phage-host interaction networks revealed an increase in CRISPR-mediated interactions coinciding with a decrease in the dominant 'Candidatus Microthrix parvicella' population. Protospacers were enriched in sequences targeting genes involved in the transmission of iMGEs. Understanding the factors shaping the fitness of specific populations is necessary to devise control strategies for undesirable species and to predict or explain community-wide phenotypes.

Stool Microbiome Profiling of Patients with Metastatic Renal Cell Carcinoma Receiving Anti-PD-1 Immune Checkpoint Inhibitors.

Preclinical models and early clinical data suggest an interplay between the gut microbiome and response to immunotherapy in solid tumors including metastatic renal cell carcinoma (mRCC). We sought to characterize the stool microbiome of mRCC patients receiving a checkpoint inhibitor (CPI) and to assess treatment-related changes in microbiome composition over the course of CPI therapy. Stool was collected from 31 patients before initiation of nivolumab (77%) or nivolumab plus ipilimumab (23%) therapy, of whom 58% experienced clinical benefit. Greater microbial diversity was associated with clinical benefit from CPI therapy (p = 0.001), and multiple species were associated with clinical benefit or lack thereof. Temporal profiling of the microbiome indicated that the relative abundance of Akkermansia muciniphila increased in patients deriving clinical benefit from CPIs. This study substantiates results from previous CPI-related microbiome profiling studies in mRCC. Temporal changes in microbiome composition suggest potential utility in modulating the microbiome for more successful CPI outcomes. PATIENT SUMMARY: We compared the composition and diversity of the gut microbiome in patients receiving immunotherapy for renal cell carcinoma. We found that higher microbial diversity is associated with better treatment outcomes. Treatment response is characterized by changes in microbial species over the course of treatment.

Integration of time-series meta-omics data reveals how microbial ecosystems respond to disturbance.

The development of reliable, mixed-culture biotechnological processes hinges on understanding how microbial ecosystems respond to disturbances. Here we reveal extensive phenotypic plasticity and niche complementarity in oleaginous microbial populations from a biological wastewater treatment plant. We perform meta-omics analyses (metagenomics, metatranscriptomics, metaproteomics and metabolomics) on in situ samples over 14 months at weekly intervals. Based on 1,364 de novo metagenome-assembled genomes, we uncover four distinct fundamental niche types. Throughout the time-series, we observe a major, transient shift in community structure, coinciding with substrate availability changes. Functional omics data reveals extensive variation in gene expression and substrate usage amongst community members. Ex situ bioreactor experiments confirm that responses occur within five hours of a pulse disturbance, demonstrating rapid adaptation by specific populations. Our results show that community resistance and resilience are a function of phenotypic plasticity and niche complementarity, and set the foundation for future ecological engineering efforts.

Mapping and comparing bacterial microbiota in the sinonasal cavity of healthy, allergic rhinitis, and chronic rhinosinusitis subjects.

BACKGROUND: The role of microbiota in sinonasal inflammation can be further understood by targeted sampling of healthy and diseased subjects. We compared the microbiota of the middle meatus (MM) and inferior meatus (IM) in healthy, allergic rhinitis (AR), and chronic rhinosinusitis (CRS) subjects to characterize intrasubject, intersubject, and intergroup differences. METHODS: Subjects were recruited in the office, and characterized into healthy, AR, and CRS groups. Endoscopically-guided swab samples were obtained from the MM and IM bilaterally. Bacterial microbiota were characterized by sequencing the V3-V4 region of the 16S ribosomal RNA (rRNA) gene. RESULTS: Intersubject microbiome analyses were conducted in 65 subjects: 8 healthy, 11 AR, and 46 CRS (25 CRS with nasal polyps [CRSwNP]; 21 CRS without nasal polyps [CRSsNP]). Intrasubject analyses were conducted for 48 individuals (4 controls, 11 AR, 8 CRSwNP, and 15 CRSwNP). There was considerable intersubject microbiota variability, but intrasubject profiles were similar (p = 0.001, nonparametric t test). Intrasubject bacterial diversity was significantly reduced in MM of CRSsNP subjects compared to IM samples (p = 0.022, nonparametric t test). CRSsNP MM samples were enriched in Streptococcus, Haemophilus, and Fusobacterium spp. but exhibited loss of diversity compared to healthy, CRSwNP, and AR subject-samples (p < 0.05; nonparametric t test). CRSwNP patients were enriched in Staphylococcus, Alloiococcus, and Corynebacterium spp. CONCLUSION: This study presents the sinonasal microbiome profile in one of the larger populations of non-CRS and CRS subjects, and is the first office-based cohort in the literature. In contrast to healthy, AR, and CRSwNP subjects, CRSsNP MM samples exhibited decreased microbiome diversity and anaerobic enrichment. CRSsNP MM samples had reduced diversity compared to same-subject IM samples, a novel finding.

First draft genome sequence of a strain belonging to the Zoogloea genus and its gene expression in situ.

The Gram-negative beta-proteobacterium Zoogloea sp. LCSB751 (LMG 29444) was newly isolated from foaming activated sludge of a municipal wastewater treatment plant. Here, we describe its draft genome sequence and annotation together with a general physiological and genomic analysis, as the first sequenced representative of the Zoogloea genus. Moreover, Zoogloea sp. gene expression in its environment is described using metatranscriptomic data obtained from the same treatment plant. The presented genomic and transcriptomic information demonstrate a pronounced capacity of this genus to synthesize poly-ß-hydroxyalkanoate within wastewater.

Phylogeographic, genomic, and meropenem susceptibility analysis of Burkholderia ubonensis.

The bacterium Burkholderia ubonensis is commonly co-isolated from environmental specimens harbouring the melioidosis pathogen, Burkholderia pseudomallei. B. ubonensis has been reported in northern Australia and Thailand but not North America, suggesting similar geographic distribution to B. pseudomallei. Unlike most other Burkholderia cepacia complex (Bcc) species, B. ubonensis is considered non-pathogenic, although its virulence potential has not been tested. Antibiotic resistance in B. ubonensis, particularly towards drugs used to treat the most severe B. pseudomallei infections, has also been poorly characterised. This study examined the population biology of B. ubonensis, and includes the first reported isolates from the Caribbean. Phylogenomic analysis of 264 B. ubonensis genomes identified distinct clades that corresponded with geographic origin, similar to B. pseudomallei. A small proportion (4%) of strains lacked the 920kb chromosome III replicon, with discordance of presence/absence amongst genetically highly related strains, demonstrating that the third chromosome of B. ubonensis, like other Bcc species, probably encodes for a nonessential pC3 megaplasmid. Multilocus sequence typing using the B. pseudomallei scheme revealed that one-third of strains lack the "housekeeping" narK locus. In comparison, all strains could be genotyped using the Bcc scheme. Several strains possessed high-level meropenem resistance (≥32 µg/mL), a concern due to potential transmission of this phenotype to B. pseudomallei. In silico analysis uncovered a high degree of heterogeneity among the lipopolysaccharide O-antigen cluster loci, with at least 35 different variants identified. Finally, we show that Asian B. ubonensis isolate RF23-BP41 is avirulent in the BALB/c mouse model via a subcutaneous route of infection. Our results provide several new insights into the biology of this understudied species.

A Bacillus anthracis Genome Sequence from the Sverdlovsk 1979 Autopsy Specimens.

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. IMPORTANCE: The 1979 Russian anthrax outbreak resulted from an industrial accident at the Soviet anthrax spore production facility in the city of Sverdlovsk. Deep genomic sequencing of two autopsy specimens generated a draft genome and phylogenetic placement of the Soviet Sverdlovsk anthrax strain. While it is known that Soviet scientists had genetically manipulated Bacillus anthracis with the potential to evade vaccine prophylaxis and antibiotic therapeutics, there was no genomic evidence of this from the Sverdlovsk production strain genome. The whole-genome SNP genotype of the Sverdlovsk strain was used to precisely identify it and its close relatives in the context of an extensive global B. anthracis strain collection. This genomic identity can now be used for forensic tracking of this weapons material on a global scale and for future anthrax investigations.