Streptococcus pyogenes can cause a wide variety of acute infections throughout the body of its human host. An underlying transcriptional regulatory network (TRN) is responsible for altering the physiological state of the bacterium to adapt to each unique host environment. Consequently, an in-depth understanding of the comprehensive dynamics of the S. pyogenes TRN could inform new therapeutic strategies. Here, we compiled 116 existing high-quality RNA sequencing data sets of invasive S. pyogenes serotype M1 and estimated the TRN structure in a top-down fashion by performing independent component analysis (ICA). The algorithm computed 42 independently modulated sets of genes (iModulons). Four iModulons contained the nga-ifs-slo virulence-related operon, which allowed us to identify carbon sources that control its expression. In particular, dextrin utilization upregulated the nga-ifs-slo operon by activation of two-component regulatory system CovRS-related iModulons, altering bacterial hemolytic activity compared to glucose or maltose utilization. Finally, we show that the iModulon-based TRN structure can be used to simplify the interpretation of noisy bacterial transcriptome data at the infection site. IMPORTANCE S. pyogenes is a pre-eminent human bacterial pathogen that causes a wide variety of acute infections throughout the body of its host. Understanding the comprehensive dynamics of its TRN could inform new therapeutic strategies. Since at least 43 S. pyogenes transcriptional regulators are known, it is often difficult to interpret transcriptomic data from regulon annotations. This study shows the novel ICA-based framework to elucidate the underlying regulatory structure of S. pyogenes allows us to interpret the transcriptome profile using data-driven regulons (iModulons). Additionally, the observations of the iModulon architecture lead us to identify the multiple regulatory inputs governing the expression of a virulence-related operon. The iModulons identified in this study serve as a powerful guidepost to further our understanding of S. pyogenes TRN structure and dynamics.

Despite progress understanding microbial communities involved in terrestrial vertebrate decomposition, little is known about the microbial decomposition of aquatic vertebrates from a functional and environmental context. Here, we analyzed temporal changes in the "necrobiome" of rainbow darters, which are common North American fish that are sensitive indicators of water quality. By combining 16S rRNA gene and shotgun metagenomic sequence data from four time points, we studied the progression of decomposers from both taxonomic and functional perspectives. The 16S rRNA gene profiles revealed strong community succession, with early decomposition stages associated with Aeromonas and Clostridium taxa and later stages dominated by members of the Rikenellaceae (i.e., Alistipes/Acetobacteroides genera). These results were reproducible and independent of environmental perturbation, given that exposure to wastewater treatment plant effluent did not substantially influence the necrobiome composition of fish or the associated water sample microbiota. Metagenomic analysis revealed significant changes throughout decomposition in degradation pathways for amino acids, carbohydrates/glycans, and other compounds, in addition to putrefaction pathways for production of putrescine, cadaverine, and indole. Binning of contigs confirmed a predominance of Aeromonas genome assemblies, including those from novel strains related to the pathogen Aeromonas veronii These bins of Aeromonas genes also encoded known hemolysin toxins (e.g., aerolysin) that were particularly abundant early in the process, potentially contributing to host cell lysis during decomposition. Overall, our results demonstrate that wild-caught fish have a reproducible decomposer succession and that the fish necrobiome serves as a potential source of putative pathogens and toxigenic bacteria.IMPORTANCE The microbial decomposition of animal tissues is an important ecological process that impacts nutrient cycling in natural environments. We studied the microbial decomposition of a common North American fish (rainbow darters) over four time points, combining 16S rRNA gene and shotgun metagenomic sequence data to obtain both taxonomic and functional perspectives. Our data revealed a strong community succession that was reproduced across different fish and environments. Decomposition time point was the main driver of community composition and functional potential; fish environmental origin (upstream or downstream of a wastewater treatment plant) had a secondary effect. We also identified strains related to the putative pathogen Aeromonas veronii as dominant members of the decomposition community. These bacteria peaked early in decomposition and coincided with the metagenomic abundance of hemolytic toxin genes. Our work reveals a strong decomposer succession in wild-caught fish, providing functional and taxonomic insights into the vertebrate necrobiome.

The Bacillus cereus group (Bacillus cereus sensu lato) has a diverse ecology, including various species that are vertebrate or invertebrate pathogens. Few isolates from the B. cereus group have however been demonstrated to benefit plant growth. Therefore, it is crucial to explore how bacterial development and pathogenesis evolve during plant colonization. Herein, we investigated Bacillus thuringiensis (Cry-) adaptation to the colonization of Arabidopsis thaliana roots and monitored changes in cellular differentiation in experimentally evolved isolates. Isolates from two populations displayed improved iterative ecesis on roots and increased virulence against insect larvae. Molecular dissection and recreation of a causative mutation revealed the importance of a nonsense mutation in the rho transcription terminator gene. Transcriptome analysis revealed how Rho impacts various B. thuringiensis genes involved in carbohydrate metabolism and virulence. Our work suggests that evolved multicellular aggregates have a fitness advantage over single cells when colonizing plants, creating a trade-off between swimming and multicellularity in evolved lineages, in addition to unrelated alterations in pathogenicity. IMPORTANCE Biologicals-based plant protection relies on the use of safe microbial strains. During application of biologicals to the rhizosphere, microbes adapt to the niche, including genetic mutations shaping the physiology of the cells. Here, the experimental evolution of Bacillus thuringiensis lacking the insecticide crystal toxins was examined on the plant root to reveal how adaptation shapes the differentiation of this bacterium. Interestingly, evolution of certain lineages led to increased hemolysis and insect larva pathogenesis in B. thuringiensis driven by transcriptional rewiring. Further, our detailed study reveals how inactivation of the transcription termination protein Rho promotes aggregation on the plant root in addition to altered differentiation and pathogenesis in B. thuringiensis.

Conventional bacterial genome annotation provides information about coding sequences but ignores untranslated regions and operons. However, untranslated regions contain important regulatory elements as well as targets for many regulatory factors, such as small RNAs. Operon maps are also essential for functional gene analysis. In the last decade, considerable progress has been made in the study of bacterial transcriptomes through transcriptome sequencing (RNA-seq). Given the compact nature of bacterial genomes, many challenges still cannot be resolved through short reads generated using classical RNA-seq because of fragmentation and loss of the full-length information. Direct RNA sequencing is a technology that sequences the native RNA directly without information loss or bias. Here, we employed direct RNA sequencing to annotate the Vibrio parahaemolyticus transcriptome with its full features, including transcription start sites (TSSs), transcription termination sites, and operon maps. A total of 4,103 TSSs were identified. In comparison to short-read sequencing, full-length information provided a deeper view of TSS classification, showing that most internal and antisense TSSs were actually a result of gene overlap. Sequencing the transcriptome of V. parahaemolyticus grown with bile allowed us to study the landscape of pathogenicity island Vp-PAI. Some genes in this region were reannotated, providing more accurate annotation to increase precision in their characterization. Quantitative detection of operons in V. parahaemolyticus showed high complexity in some operons, shedding light on a greater extent of regulation within the same operon. Our study using direct RNA sequencing provides a quantitative and high-resolution landscape of the V. parahaemolyticus transcriptome. IMPORTANCE Vibrio parahaemolyticus is a halophilic bacterium found in the marine environment. Outbreaks of gastroenteritis resulting from seafood poisoning by these pathogens have risen over the past 2 decades. Upon ingestion by humans-often through the consumption of raw or undercooked seafood-V. parahaemolyticus senses the host environment and expresses numerous genes, the products of which synergize to synthesize and secrete toxins that can cause acute gastroenteritis. To understand the regulation of such adaptive response, mRNA transcripts must be mapped accurately. However, due to the limitations of common sequencing methods, not all features of bacterial transcriptomes are always reported. We applied direct RNA sequencing to analyze the V. parahaemolyticus transcriptome. Mapping the full features of the transcriptome is anticipated to enhance our understanding of gene regulation in this bacterium and provides a data set for future work. Additionally, this study reveals a deeper view of a complicated transcriptome landscape, demonstrating the importance of applying such methods to other bacterial models.

Elimination or alteration of select members of the gut microbiota is key to therapeutic efficacy. However, the complexity of these microbial inhabitants makes it challenging to precisely target bacteria. Here, we deliver exogenous genes to specific bacteria by genomic integration of temperate phage for long-lasting modification. As a real-world therapeutic test, we engineered λ phage to transcriptionally repress Shiga toxin by using genetic hybrids between λ and other lambdoid phages to overcome resistance encoded by the virulence-expressing prophage. We show that a single dose of engineered phage propagates throughout the bacterial community and reduces Shiga toxin production in an enteric mouse model of infection without markedly affecting bacterial concentrations. Our work reveals a new framework for transferring functions to bacteria within their native environment.IMPORTANCE With the increasing frequency of antibiotic resistance, it is critical to explore new therapeutic strategies for treating bacterial infections. Here, we use a temperate phage, i.e., one that integrates itself into the bacterial genome, to neutralize the expression of a virulence factor by modifying bacterial function at the genetic level. We show that Shiga toxin production can be significantly reduced in vitro and in the mammalian gut. Alternative to traditional applications of phage therapy that rely on killing bacteria, our genetics-based antivirulence approach introduces a new framework for treating bacterial infections.

Gram-negative bacteria frequently possess type VI secretion systems (T6SSs), protein complexes that are able to inject toxic proteins into nearby cells. Many aspects of T6SS structure and function have been characterized for model species, but less is known about the evolutionary processes that shape T6SS and effector (toxin) diversity in host-associated microbial communities. The bee gut microbiota is a simple community that has codiversified with bees for >80 million years. This study investigated how complements of T6SSs and effectors within the bee microbiota changed as bacteria and their hosts diversified into isolated species. We used protein homology to survey 198 isolate genomes of 9 Gram-negative species for genes encoding T6SS structural components; Rhs toxins, which are common T6SS effectors; and VgrG proteins, which are structural components associated with specific toxins. T6SS loci were present in 5 species clusters found only in bees, namely Apibacter spp., Gilliamella spp., Frischella perrara, "Candidatus Schmidhempelia bombi," and Snodgrassella alvi The distribution of T6SS loci suggests that at least 3 were present in the microbiota of the common ancestor of social bees and that loss of these genes in some bacterial lineages was linked to both host and bacterial speciation. Isolates differed enormously in repertoires of Rhs and VgrG proteins. We found that bacterial species employ different mechanisms for toxin acquisition and diversification and that species and strains sometimes lose the T6SS entirely, likely causing shifts in competitive dynamics within these communities.IMPORTANCE Antagonistic interactions between bacteria affect diversity and dynamics of host-associated communities, including gut communities that are linked to host health. In many bacterial communities, including human and honey bee gut microbiotas, antagonism is mediated by type VI secretion systems (T6SSs) that deliver lethal toxins to competing strains. In this study, we explored how T6SSs and associated toxins have evolved in the simple, host-specific gut microbiota of honey bees and bumble bees. Using comparative genomics, we explored the conservation, recombination, horizontal transfer, and loss of T6SSs and effectors during 80 million years of evolution of this bee-associated community. We find that that patterns of T6SS loss and retention are linked to differences in biology across host species, while trends in effector diversification are mostly specific to bacterial lineages.

Type VI secretion systems (T6SSs) play a major role in interbacterial competition and in bacterial interactions with eukaryotic cells. The distribution of T6SSs and the effectors they secrete vary between strains of the same bacterial species. Therefore, a pan-genome investigation is required to better understand the T6SS potential of a bacterial species of interest. Here, we performed a comprehensive, systematic analysis of T6SS gene clusters and auxiliary modules found in the pan-genome of Vibrio parahaemolyticus, an emerging pathogen widespread in marine environments. We identified 4 different T6SS gene clusters within genomes of this species; two systems appear to be ancient and widespread, whereas the other 2 systems are rare and appear to have been more recently acquired via horizontal gene transfer. In addition, we identified diverse T6SS auxiliary modules containing putative effectors with either known or predicted toxin domains. Many auxiliary modules are possibly horizontally shared between V. parahaemolyticus genomes, since they are flanked by DNA mobility genes. We further investigated a DUF4225-containing protein encoded on an Hcp auxiliary module, and we showed that it is an antibacterial T6SS effector that exerts its toxicity in the bacterial periplasm, leading to cell lysis. Computational analyses of DUF4225 revealed a widespread toxin domain associated with various toxin delivery systems. Taken together, our findings reveal a diverse repertoire of T6SSs and auxiliary modules in the V. parahaemolyticus pan-genome, as well as novel T6SS effectors and toxin domains that can play a major role in the interactions of this species with other cells. IMPORTANCE Gram-negative bacteria employ toxin delivery systems to mediate their interactions with neighboring cells. Vibrio parahaemolyticus, an emerging pathogen of humans and marine animals, was shown to deploy antibacterial toxins into competing bacteria via the type VI secretion system (T6SS). Here, we analyzed 1,727 V. parahaemolyticus genomes and revealed the pan-genome T6SS repertoire of this species, including the T6SS gene clusters, horizontally shared auxiliary modules, and toxins. We also identified a role for a previously uncharacterized domain, DUF4225, as a widespread antibacterial toxin associated with diverse toxin delivery systems.

Clostridioides difficile is a Gram-positive, anaerobic, spore-forming bacterium responsible for antibiotic-associated pseudomembranous colitis. Clostridioides difficile infection (CDI) symptoms can range from diarrhea to life-threatening colon damage. Toxins produced by C. difficile (TcdA and TcdB) cause intestinal epithelial injury and lead to severe gut barrier dysfunction, stem cell damage, and impaired regeneration of the gut epithelium. Current treatment options for intestinal repair are limited. In this study, we demonstrate that treatment with the microbial metabolite urolithin A (UroA) attenuates CDI-induced adverse effects on the colon epithelium in a preclinical model of CDI-induced colitis. Moreover, our analysis suggests that UroA treatment protects against C. difficile-induced inflammation, disruption of gut barrier integrity, and intestinal tight junction proteins in the colon of CDI mice. Importantly, UroA treatment significantly reduced the expression and release of toxins from C. difficile without inducing bacterial cell death. These results indicate the direct regulatory effects of UroA on bacterial gene regulation. Overall, our findings reveal a novel aspect of UroA activity, as it appears to act at both the bacterial and host levels to protect against CDI-induced colitis pathogenesis. This research sheds light on a promising avenue for the development of novel treatments for C. difficile infection.IMPORTANCETherapy for Clostridioides difficile infections includes the use of antibiotics, immunosuppressors, and fecal microbiota transplantation. However, these treatments have several drawbacks, including the loss of colonization resistance, the promotion of autoimmune disorders, and the potential for unknown pathogens in donor samples. To date, the potential benefits of microbial metabolites in CDI-induced colitis have not been fully investigated. Here, we report for the first time that the microbial metabolite urolithin A has the potential to block toxin production from C. difficile and enhance gut barrier function to mitigate CDI-induced colitis.

Protein toxins secreted by prokaryotes have been found to affect the pathogenicity of pathogens or directly mediate antagonistic interactions between prokaryotes. PAAR proteins are important carriers of toxic effectors and are located at the forefront of either the type VI secretion system (T6SS) or the extracellular contractile injection system (eCIS). This study systematically investigated PAAR homologues and related toxic effectors. We found that PAAR homologues were divided into 8 types and 16 subtypes and distributed in 23.1% of bacterial genomes and 7.8% of archaeal genomes. PAAR proteins of all types fold into a highly similar conical structure, even from relatively diverse underlying sequences. PAAR homologues associated with different secretion systems display a mixed phylogenetic relationship, indicating that PAAR proteins from such a subtype can be assembled on either a T6SS or an eCIS. More than 1,300 PAAR-related toxic effector genes were identified; one PAAR subtype can be associated with toxins of over 40 families, and toxins from one family can be associated with more than 10 PAAR subtypes. A large-scale comparison of Earth Microbiome Project data and prokaryotic genomes revealed that prokaryotes encoding PAAR genes are widely present in diverse environments worldwide, and taxa encoding multiple PAAR gene copies exhibit a wider distribution in environments than other taxa. Overall, our studies highlighted that PAAR proteins are versatile clips loaded with antimicrobial toxin bullets for secretion weapons (T6SS and eCIS), greatly enriching the weapon arsenal of prokaryotes, which, often together with VgrG, help prokaryotes fight for survival advantages in crowded environments. IMPORTANCE Infectious diseases caused by microbial pathogens are severe threats to human health and economic development. To respond to these threats, it is necessary to understand how microorganisms survive in and adapt to complex environments. Microorganic toxins, which are widely distributed in nature, are the key weapons in life domain interactions. PAAR proteins are important carriers of prokaryotic toxic effectors. We reveal the versatility of PAAR proteins between secretory systems and the massive diversity of toxic effectors carried by PAAR proteins, which helps prokaryotes enrich their arsenal and expand their ability to attack their neighbors. A large number of PAAR homologues and related toxic effectors enhance the survival competitiveness of prokaryotic populations. In conclusion, our work provides an example for large-scale analysis of the global distribution and ecological functions of prokaryotic functional genes.

The innate immune system is the body's first line of defense against pathogens and its protection against infectious diseases. On the surface of host myeloid cells, Toll-like receptor 4 (TLR4) senses lipopolysaccharide (LPS), the major outer membrane component of Gram-negative bacteria. Intracellularly, LPS is recognized by caspase 11 through the noncanonical inflammasome to induce pyroptosis-an inflammatory form of lytic cell death. While TLR4-mediated signaling perturbations result in secretion of cytokines and chemokines that help clear infection and facilitate adaptive immunity, caspase 11-mediated pyroptosis leads to the release of damage-associated molecular patterns and inflammatory mediators. Although the core signaling events and many associated proteins in the TLR4 signaling pathway are known, the complex signaling events and protein networks within the noncanonical inflammasome pathway remain obscure. Moreover, there is mounting evidence for pathogen-specific innate immune tuning. We characterized the major LPS structures from two different pathogens, modeled their binding to the surface receptors, systematically examined macrophage inflammatory responses to these LPS molecules, and surveyed the temporal differences in global protein secretion resulting from TLR4 and caspase 11 activation in macrophages using mass spectrometry (MS)-based quantitative proteomics. This integrated strategy, spanning functional activity assays, top-down structural elucidation of endotoxins, and secretome analysis of stimulated macrophages, allowed us to identify crucial differences in TLR4- and caspase 11-mediated protein secretion in response to two Gram-negative bacterial endotoxins. IMPORTANCE Macrophages and monocytes are innate immune cells playing an important role in orchestrating the initial innate immune response to bacterial infection and the tissue damage. This response is facilitated by specific receptors on the cell surface and intracellularly. One of the bacterial molecules recognized is a Gram-negative bacteria cell wall component, lipopolysaccharide (LPS). The structure of LPS differs between different species. We have characterized the innate immune responses to the LPS molecules from two bacteria, Escherichia coli and Bordetella pertussis, administered either extracellularly or intracellularly, whose structures we first determined. We observed marked differences in the temporal dynamics and amounts of proteins secreted by the innate immune cells stimulated by any of these molecules and routes. This suggests that there is specificity in the first line of response to different Gram-negative bacteria that can be explored to tailor specific therapeutic interventions.

Yersinia enterocolitica is a zoonotic pathogen and an important cause of bacterial gastrointestinal infections in humans. Large-scale population genomic analyses revealed genetic and phenotypic diversity of this bacterial species, but little is known about the differences in the transcriptome organization, small RNA (sRNA) repertoire, and transcriptional output. Here, we present the first comparative high-resolution transcriptome analysis of Y. enterocolitica strains representing highly pathogenic phylogroup 2 (serotype O:8) and moderately pathogenic phylogroup 3 (serotype O:3) grown under four infection-relevant conditions. Our transcriptome sequencing (RNA-seq) approach revealed 1,299 and 1,076 transcriptional start sites and identified strain-specific sRNAs that could contribute to differential regulation among the phylogroups. Comparative transcriptomics further uncovered major gene expression differences, in particular, in the temperature-responsive regulon. Multiple virulence-relevant genes are differentially regulated between the two strains, supporting an ecological separation of phylogroups with certain niche-adapted properties. Strong upregulation of the ystA enterotoxin gene in combination with constitutive high expression of cell invasion factor InvA further showed that the toxicity of recent outbreak O:3 strains has increased. Overall, our report provides new insights into the specific transcriptome organization of phylogroups 2 and 3 and reveals gene expression differences contributing to the substantial phenotypic differences that exist between the lineages. IMPORTANCE Yersinia enterocolitica is a major diarrheal pathogen and is associated with a large range of gut-associated diseases. Members of this species have evolved into different phylogroups with genotypic variations. We performed the first characterization of the Y. enterocolitica transcriptional landscape and tracked the consequences of the genomic variations between two different pathogenic phylogroups by comparing their RNA repertoire, promoter usage, and expression profiles under four different virulence-relevant conditions. Our analysis revealed major differences in the transcriptional outputs of the closely related strains, pointing to an ecological separation in which one is more adapted to an environmental lifestyle and the other to a mostly mammal-associated lifestyle. Moreover, a variety of pathoadaptive alterations, including alterations in acid resistance genes, colonization factors, and toxins, were identified which affect virulence and host specificity. This illustrates that comparative transcriptomics is an excellent approach to discover differences in the functional output from closely related genomes affecting niche adaptation and virulence, which cannot be directly inferred from DNA sequences.

The complexity of microbial communities hinders our understanding of how microbial diversity and microbe-microbe interactions impact community functions. Here, using six independent communities originating from the refuse dumps of leaf-cutter ants and enriched using the plant polymer cellulose as the sole source of carbon, we examine how changes in bacterial diversity and interactions impact plant biomass decomposition. Over up to 60 serial transfers (∼8 months) using Whatman cellulose filter paper, cellulolytic ability increased and then stabilized in four enrichment lines and was variable in two lines. Bacterial community characterization using 16S rRNA gene amplicon sequencing showed community succession differed between the highly cellulolytic enrichment lines and those that had slower and more variable cellulose degradation rates. Metagenomic and metatranscriptomic analyses revealed that Cellvibrio and/or Cellulomonas dominated each enrichment line and produced the majority of cellulase enzymes, while diverse taxa were retained within these communities over the duration of transfers. Interestingly, the less cellulolytic communities had a higher diversity of organisms competing for the cellulose breakdown product cellobiose, suggesting that cheating slowed cellulose degradation. In addition, we found competitive exclusion as an important factor shaping all of the communities, with a negative correlation of Cellvibrio and Cellulomonas abundance within individual enrichment lines and the expression of genes associated with the production of secondary metabolites, toxins, and other antagonistic compounds. Our results provide insights into how microbial diversity and competition affect the stability and function of cellulose-degrading communities. IMPORTANCE Microbial communities are a key driver of the carbon cycle through the breakdown of complex polysaccharides in diverse environments including soil, marine systems, and the mammalian gut. However, due to the complexity of these communities, the species-species interactions that impact community structure and ultimately shape the rate of decomposition are difficult to define. Here, we performed serial enrichment on cellulose using communities inoculated from leaf-cutter ant refuse dumps, a cellulose-rich environment. By concurrently tracking cellulolytic ability and community composition and through metagenomic and metatranscriptomic sequencing, we analyzed the ecological dynamics of the enrichment lines. Our data suggest that antagonism is prevalent in these communities and that competition for soluble sugars may slow degradation and lead to community instability. Together, these results help reveal the relationships between competition and polysaccharide decomposition, with implications in diverse areas ranging from microbial community ecology to cellulosic biofuels production.

AtxA, the master virulence regulator of Bacillus anthracis, regulates the expression of three toxins and genes for capsule formation that are required for the pathogenicity of B. anthracis. Recent transcriptome analyses showed that AtxA affects a large number of genes on the chromosome and plasmids, suggesting a role as a global regulator. However, information on genes directly regulated by AtxA is scarce. In this work, we conducted genome-wide analyses and cataloged the binding sites of AtxA in vivo and transcription start sites on the B. anthracis genome. By integrating these results, we detected eight genes as direct regulons of AtxA. These consisted of five protein-coding genes, including two of the three toxin genes, and three genes encoding the small RNAs XrrA and XrrB and a newly discovered 95-nucleotide small RNA, XrrC. Transcriptomes from single-knockout mutants of these small RNAs revealed changes in the transcription levels of genes related to the aerobic electron transport chain, heme biosynthesis, and amino acid metabolism, suggesting their function for the control of cell physiology. These results reveal the first layer of the gene regulatory network for the pathogenicity of B. anthracis and provide a data set for the further study of the genomics and genetics of B. anthracis. IMPORTANCE Bacillus anthracis is the Gram-positive bacterial species that causes anthrax. Anthrax is still prevalent in countries mainly in Asia and Africa, where it causes economic damage and remains a public health issue. The mechanism of pathogenicity is mainly explained by the three toxin proteins expressed from the pXO1 plasmid and by proteins involved in capsule formation expressed from the pXO2 plasmid. AtxA is a protein expressed from the pXO1 plasmid that is known to upregulate genes involved in toxin production and capsule formation and is thus considered the master virulence regulator of B. anthracis. Therefore, understanding the detailed mechanism of gene regulation is important for the control of anthrax. The significance of this work lies in the identification of genes that are directly regulated by AtxA via genome-wide analyses. The results reveal the first layer of the gene regulatory network for the pathogenicity of B. anthracis and provide useful resources for a further understanding of B. anthracis.

For mammals, oxalate enters the body through the diet or is endogenously produced by the liver; it is removed by microbial oxalate metabolism in the gut and/or excretion in feces or urine. Deficiencies in any one of the these pathways can lead to complications, such as calcium oxalate urinary stones. While considerable research has been conducted on individual oxalate-degrading bacterial isolates, interactions between oxalate and the gut microbiota as a whole are unknown. We examined the reduction in oxalate excretion in a rat model following oral administration of fecal microbes from a mammalian herbivore adapted to a high oxalate diet or to fecal transplants consisting of two different formulations of mixed oxalate-degrading isolates. While all transplants elicited a significant reduction in oxalate excretion initially, the greatest effect was seen with fecal microbial transplants, which persisted even in the absence of dietary oxalate. The reduction in oxalate excretion in animals given fecal transplants corresponded with the establishment of diverse bacteria, including known oxalate-degrading bacteria and a cohesive network of bacteria centered on oxalate-degrading specialists from the Oxalobacteraceae family. Results suggested that the administration of a complete community of bacteria facilitates a cohesive balance in terms of microbial interactions. Our work offers important insights into the development of targeted bacteriotherapies intended to reduce urinary oxalate excretion in patients at risk for recurrent calcium oxalate stones as well as bacteriotherapies targeting other toxins for elimination. IMPORTANCE Oxalate is a central component in 80% of kidney stones. While mammals do not possess the enzymes to degrade oxalate, many gastrointestinal bacteria are efficient oxalate degraders. We examined the role of cohesive microbial networks for oxalate metabolism, using Sprague-Dawley rats as a model host. While the transplantation of oxalate-degrading bacteria alone to the Sprague-Dawley hosts did increase oxalate metabolism, fecal transplants from a wild mammalian herbivore, Neotoma albigula, had a significantly greater effect. Furthermore, the boost for oxalate metabolism persisted only in animals that received fecal transplants. Animals receiving fecal transplants had a more diverse and cohesive network of bacteria associated with the Oxalobacteraceae, a family known to consist of specialist oxalate-degrading bacteria, than did animals that received oxalate-degrading bacteria alone. Our results indicate that fecal transplants are more effective at transferring specific functions than are microbial specialists alone, which has broad implications for the development of bacteriotherapies.

Non-aureus staphylococci (NAS) are the most frequently isolated pathogens from intramammary infection (IMI) in dairy cattle. Virulence factors (VFs) and mechanisms by which NAS cause IMI are not fully known. Herein, we analyzed the distribution of 191 VFs in 441 genomes of 25 NAS species, after classifying VFs into functional categories: adherence (n = 28), exoenzymes (n = 21), immune evasion (n = 20), iron metabolism (n = 29), and toxins (n = 93). In addition to establishing VF gene profiles, associations of VF genes between and among functional categories were computed, revealing distinctive patterns of association among VFs for various NAS species. Associations were also computed for low, medium, and high somatic cell count (SCC) and clinical mastitis (CM) isolates, demonstrating distinctive patterns of associations for low SCC and CM isolates, but no differences between high SCC and CM isolates. To determine whether VF distributions had any association with SCC or CM, various clustering approaches, including complete linkages, Ward clustering, and t-distributed stochastic neighbor embedding, were applied. However, no clustering of isolates representing low SCC, medium SCC, or high SCC or CM was identified. Regression analysis to test for associations with individual VF functional categories demonstrated that each additional toxin and host immune evasion gene increased the odds of having high SCC or CM, although an overall increase in the number of VFs was not associated with increased SCC or occurrence of CM. In conclusion, we established comprehensive VF gene profiling, determined VF gene distributions and associations, calculated pathogenic potentials of all NAS species, and detected no clear link between VF genes and mastitis. IMPORTANCE Non-aureus staphylococci (NAS) are the most frequently isolated pathogens from milk in dairy cattle worldwide. The virulence factors (VFs) and mechanisms by which these bacteria cause udder infection are not fully known. We determined the distribution and associations of 191 VFs in 25 NAS species and investigated the relationship between VFs and disease. Although the overall number of VFs was not associated with disease severity, increasing numbers of toxin and host immune evasion genes specifically were associated with more severe disease outcomes. These findings suggest that the development of disease and the interactions of VFs with the host are complex and determined by the interplay of genes rather than just the presence of virulence genes. Together, our results provide foundational genetic knowledge to other researchers to design and conduct further experiments, focusing on understanding the synergy between VFs and roles of individual NAS species in IMI and characterizing species-specific effects on udder health.

Staphylococcus aureus causes persistent clinical and subclinical bovine intramammary infections (IMI) worldwide. However, there is a lack of comprehensive information regarding genetic diversity, the presence of antimicrobial resistance (AMR), and virulence genes for S. aureus in bovine milk in Canada. Here, we performed whole-genome sequencing (WGS) of 119 Canadian bovine milk S. aureus isolates and determined they belonged to 8 sequence types (ST151, ST352, ST351, ST2187, ST2270, ST126, ST133, and ST8), 5 clonal complexes (CC151, CC97, CC126, CC133, and CC8), and 18 distinct Spa types. Pan-, core, and accessory genomes were composed of 6,340, 1,279, and 2,431 genes, respectively. Based on phenotypic screening for AMR, resistance was common against beta-lactams (19% of isolates) and sulfonamides (7% of isolates), whereas resistance against pirlimycin, tetracycline, ceftiofur, and erythromycin and to the combination of penicillin and novobiocin was uncommon (3, 3, 3, 2, and 2% of all isolates, respectively). We also determined distributions of 191 virulence factors (VFs) in 119 S. aureus isolates after classifying them into 5 functional categories (adherence [n = 28], exoenzymes [n = 21], immune evasion [n = 20], iron metabolism [n = 29], and toxins [n = 93]). Additionally, we calculated the pathogenic potential of distinct CCs and STs and determined that CC151 (ST151 and ST351) had the highest pathogenic potential (calculated by subtracting core-VFs from total VFs), followed by CC97 (ST352 and ST2187) and CC126 (ST126 and ST2270), potentially linked to their higher prevalence in bovine IMI worldwide. However, there was no statistically significant link between the presence of VF genes and mastitis.IMPORTANCE Staphylococcus aureus is a major cause of bovine intramammary infections, leading to significant economic losses to dairy industry in Canada and worldwide. There is a lack of knowledge regarding genetic diversity, the presence of antimicrobial resistance (AMR), and virulence genes for S. aureus isolated from bovine milk in Canada. Based on whole-genome sequencing and genomic analysis, we have determined the phylogeny and diversity of S. aureus in bovine milk and concluded that it had a large accessory genome, limited distribution of AMR genes, variable VF gene profiles and sequence types (ST), and clonal complex (CC)-specific pathogenic potentials. Comprehensive information on the population structure, as well as the virulence and resistance characteristics of S. aureus from bovine milk, will allow for source attribution, risk assessment, and improved therapeutic approaches in cattle.