Search

The molecular mechanism of the Ca(2+)-mediated formation of competent cells in Escherichia coli remains unclear. In this study, transcriptome and proteomics techniques were used to screen genes in response to Ca(2+) treatment. A total of 333 differentially expressed genes (317 upregulated and 16 downregulated) and 145 differentially expressed proteins (54 upregulated and 91 downregulated) were obtained. These genes and proteins are mainly enriched in cell membrane components, transmembrane transport, and stress response-related functional terms. Fifteen genes with these functions, including yiaW, ygiZ, and osmB, are speculated to play a key role in the cellular response to Ca(2+). Three single-gene deletion strains were constructed with the Red homologous recombination method to verify its function in genetic transformation. The transformation efficiencies of yiaW, ygiZ, and osmB deletion strains for different-size plasmids were significantly increased. None of the three gene deletion strains changed in size, which is one of the main elements of microscopic morphology, but they exhibited different membrane permeabilities and transformation efficiencies. This study demonstrates that Ca(2+)-mediated competence formation in E. coli is not a simple physicochemical process and may involve the regulation of genes in response to Ca(2+). This study lays the foundation for further in-depth analyses of the molecular mechanism of Ca(2+)-mediated transformation. IMPORTANCE Using transcriptome and proteome techniques and association analysis, we identified several key genes involved in the formation of Ca(2+)-mediated E. coli DH5α competent cells. We used Red homologous recombination technology to construct three single-gene deletion strains and found that the transformation efficiencies of yiaW, ygiZ, and osmB deletion strains for different-size plasmids were significantly increased. These results proved that the genetic transformation process is not only a physicochemical process but also a reaction process involving multiple genes. These results suggest ways to improve the horizontal gene transfer mechanism of foodborne microorganisms and provide new ideas for ensuring the safety of food preservation and processing.

Xylella fastidiosa is an important bacterial pathogen of plants causing high consequence diseases in agricultural crops around the world. Although as a species X. fastidiosa can infect an extremely broad range of host plants, significant variability exists between strains and subspecies groups in virulence on specific host plant species, and other traits such as growth habits. Natural competence and horizontal gene transfer are believed to occur frequently in X. fastidiosa, and likely influences the evolution of this pathogen. However, some X. fastidiosa strains are extremely difficult or impossible to manipulate genetically using standard transformation techniques. Several restriction-modification systems are encoded in the X. fastidiosa genome, including multiple Type I R-M systems that may influence horizontal gene transfer and recombination. In this study, several conserved Type I R-M systems were compared across 129 X. fastidiosa genome assemblies representing all known subspecies and 32 sequence types. Considerable allelic variation among strains was identified among the single specificity subunit (hsdS) of each Type I R-M system, with a unique hsdS allele profile generally conserved within a monophyletic cluster of strains. Inactivating mutations were identified in Type I R-M systems of specific strains, showing heterogeneity in the complement of functional Type I R-M systems across X. fastidiosa. Genomic DNA methylation patterns were characterized in 20 X. fastidiosa strains and associated with Type I R-M system allele profiles. Overall, this study describes epigenetic modifications in X. fastidiosa associated with functional Type I R-M systems and characterizes the diversity in these systems across X. fastidiosa lineages.ImportanceEconomic impacts on agricultural production due to X. fastidiosa have been severe in the Americas, Europe, and parts of Asia. Despite a long history of research on this pathogen, certain fundamental questions regarding the biology, pathogenicity, and evolution of X. fastidiosa have still not been answered. Wide scale whole genome sequencing has begun to provide a more insight into X. fastidiosa genetic diversity and horizontal gene transfer but the mechanics of genomic recombination in natural settings and extent to which this directly influences bacterial phenotypes such as plant host range are not well understood. Genome methylation is an important factor in horizontal gene transfer and bacterial recombination that has not been comprehensively studied in X. fastidiosa. This study characterizes methylation associated with Type I restriction-modification systems across a wide range of X. fastidiosa strains and lays the groundwork for a better understanding of X. fastidiosa biology and evolution through epigenetics.

Enterococcus faecalis is an important intestinal colonizing bacteria and can cause various tissue infections, including invasive blood infection (BI). The annual incidence of E. faecalis BI has been estimated to be ~4.5 per 100,000, with a fatality rate that can reach 20%. However, whether bacterial colonization or invasive infections are tissue based has not been thoroughly studied. In this study, we analyzed 537 clinical isolates from 7 different tissues to identify the key genomic elements that facilitate the colonization and invasive infection of E. faecalis. Comparative genomic analysis revealed that the BI E. faecalis isolates had the largest genome size but the lowest GC content, fsr quorum-sensing system genes were enriched in the BI E. faecalis, and the fsr gene cluster could enhance biofilm formation and serum resistance ability. Our findings also provide deep insight into the genomic differences between different tissue isolates, and the fsr quorum-sensing systems could be a key factor promoting E. faecalis invasion into the blood. IMPORTANCE First, we conducted an advanced study on the genomic differences between colonizing and infecting E. faecalis, which provides support and evidence for early and accurate diagnoses. Second, we discovered that fsr was significantly associated with blood infections, which also provides additional information for studies exploring the invasiveness of E. faecalis. Most importantly, we found that fsr played an important role in both biofilm formation and serum resistance ability in E. faecalis.

Angucyclines are a family of structurally diverse, aromatic polyketides with some members that exhibit potent bioactivity. Angucyclines have also attracted considerable attention due to the intriguing biosynthetic origins that underlie their structural complexity and diversity. Balmoralmycin (compound 1) represents a unique group of angucyclines that contain an angular benz[α]anthracene tetracyclic system, a characteristic C-glycosidic bond-linked deoxy-sugar (d-olivose), and an unsaturated fatty acid chain. In this study, we identified a Streptomyces strain that produces balmoralmycin and seven biosynthetically related coproducts (compounds 2−8). Four of the coproducts (compounds 5−8) are novel compounds that feature a highly oxygenated or fragmented lactone ring, and three of them (compounds 3−5) exhibited cytotoxicity against the human pancreatic cancer cell line MIA PaCa-2 with IC(50) values ranging from 0.9 to 1.2 μg/mL. Genome sequencing and CRISPR/dCas9-assisted gene knockdown led to the identification of the ~43 kb balmoralmycin biosynthetic gene cluster (bal BGC). The bal BGC encodes a type II polyketide synthase (PKS) system for assembling the angucycline aglycone, six enzymes for generating the deoxysugar d-olivose, and a hybrid type II/III PKS system for synthesizing the 2,4-decadienoic acid chain. Based on the genetic and chemical information, we propose a mechanism for the biosynthesis of balmoralmycin and the shunt products. The chemical and genetic studies yielded insights into the biosynthetic origin of the structural diversity of angucyclines. IMPORTANCE Angucyclines are structurally diverse aromatic polyketides that have attracted considerable attention due to their potent bioactivity and intriguing biosynthetic origin. Balmoralmycin is a representative of a small family of angucyclines with unique structural features and an unknown biosynthetic origin. We report a newly isolated Streptomyces strain that produces balmoralmycin in a high fermentation titer as well as several structurally related shunt products. Based on the chemical and genetic information, a biosynthetic pathway that involves a type II polyketide synthase (PKS) system, cyclases/aromatases, oxidoreductases, and other ancillary enzymes was established. The elucidation of the balmoralmycin pathway enriches our understanding of how structural diversity is generated in angucyclines and opens the door for the production of balmoralmycin derivatives via pathway engineering.

The exchange of bacterial extracellular vesicles facilitates molecular exchange between cells, including the horizontal transfer of genetic material. Given the implications of such transfer events on cell physiology and adaptation, some bacterial cells have likely evolved mechanisms to regulate vesicle exchange. Past work has identified mechanisms that influence the formation of extracellular vesicles, including the production of small molecules that modulate membrane structure; however, whether these mechanisms also modulate vesicle uptake and have an overall impact on the rate of vesicle exchange is unknown. Here, we show that membrane-binding molecules produced by microbes influence both the formation and uptake of extracellular vesicles and have the overall impact of increasing the vesicle exchange rate within a bacterial coculture. In effect, production of compounds that increase vesicle exchange rates encourage gene exchange between neighboring cells. The ability of several membrane-binding compounds to increase vesicle exchange was demonstrated. Three of these compounds, nisin, colistin, and polymyxin B, are antimicrobial peptides added at sub-inhibitory concentrations. These results suggest that a potential function of exogenous compounds that bind to membranes may be the regulation of vesicle exchange between cells. IMPORTANCE The exchange of bacterial extracellular vesicles is one route of gene transfer between bacteria, although it was unclear if bacteria developed strategies to modulate the rate of gene transfer within vesicles. In eukaryotes, there are many examples of specialized molecules that have evolved to facilitate the production, loading, and uptake of vesicles. Recent work with bacteria has shown that some small molecules influence membrane curvature and induce vesicle formation. Here, we show that similar compounds facilitate vesicle uptake, thereby increasing the overall rate of vesicle exchange within bacterial populations. The addition of membrane-binding compounds, several of them antibiotics at subinhibitory concentrations, to a bacterial coculture increased the rate of horizontal gene transfer via vesicle exchange.

The advent of high-throughput sequencing and population genomics has enabled researchers to investigate selection pressure at hypervariable genomic loci encoding pathogen-associated molecular pattern (PAMP) molecules like lipopolysaccharide (LPS). Xanthomonas is a model and a major group of phytopathogenic bacteria that infect hosts in tissue-specific manner. Our in-depth population-based genomic investigation revealed the emergence of major lineages in two Xanthomonas pathogens that infect xylem of rice and sugarcane is associated with the acquisition and later large-scale replacement by distinct type of LPS cassettes. In the population of the rice xylem pathogen, Xanthomonas oryzae pv. oryzae (Xoo) and sugarcane pathogens Xanthomonas sacchari (Xsac) and Xanthomonas vasicola (Xvv), the BXO8 type of LPS cassette is replaced by a BXO1 type of cassette in Xoo and by Xvv type LPS cassette in Xsac and Xvv. These findings suggest a wave of parallel evolution at an LPS locus mediated by horizontal gene transfer (HGT) events during its adaptation and emergence. Aside from xylem pathogens, two closely related lineages of Xoo that infect parenchyma of rice and Leersia hexandra grass have acquired an LPS cassette from Xanthomonas pathogens that infect parenchyma of citrus, walnut, and strawberries, indicating yet another instance of parallel evolution mediated by HGT at an LPS locus. Our targeted and megapopulation-based genome dynamic studies revealed the acquisition and dominance of specific types of LPS cassettes in adaptation and success of a major group of phytopathogenic bacteria. IMPORTANCE Lipopolysaccharide (LPS) is a major microbe associated molecular pattern and hence a major immunomodulator. As a major and outer member component, it is expected that LPS is a frontline defense mechanism to deal with different host responses. Limited studies have indicated that LPS loci are also highly variable at strain and species level in plant-pathogenic bacteria, suggesting strong selection pressure from plants and associated niches. The advent of high-throughput genomics has led to the availability of a large set of genomic resources at taxonomic and population levels. This provides an exciting and important opportunity to carryout megascale targeted and population-based comparative genomic/association studies at important loci like those encoding LPS biosynthesis to understand their role in the evolution of the host, tissue specificity, and also predominant lineages. Such studies will also fill major gap in understanding host and tissue specificity in pathogenic bacteria. Our pioneering study uses the Xanthomonas group of phytopathogens that are known for their characteristic host and tissue specificity. The present deep phylogenomics of diverse Xanthomonas species and its members revealed lineage association and dominance of distinct types of LPS in accordance with their origin, host, tissue specificity, and evolutionary success.

In Europe, genetically distinct ecotypes of the tick-vectored bacterium Anaplasma phagocytophilum circulate among mammals in three discrete enzootic cycles. To date, potential ecological factors that contributed to the emergence of these divergent ecotypes have been poorly studied. Here, we show that the ecotype that predominantly infects roe deer (Capreolus capreolus) is evolutionarily derived. Its divergence from a host generalist ancestor occurred after the last glacial maximum as mammal populations, including roe deer, recolonized the European mainland from southern refugia. We also provide evidence that this host specialist ecotype’s effective population size (N(e)) has tracked changes in the population of its roe deer host. Specifically, both host and bacterium have undergone substantial increases in N(e) over the past 1,500 years. In contrast, we show that while it appears to have undergone a major population expansion starting ~3,500 years ago, in the past 500 years, the contemporary host generalist ecotype has experienced a substantial reduction in genetic diversity levels, possibly as a result of reduced opportunities for transmission between competent hosts. IMPORTANCE The findings of this study reveal specific events important for the evolution of host specialization in a naturally occurring, obligately intracellular bacterial pathogen. Specifically, they show that host range shifts and the emergence of host specialization may occur during periods of population growth in a generalist ancestor. Our results also demonstrate the close correlation between demographic patterns in host and pathogen for a specialist system. These findings have important relevance for understanding the evolution of host range diversity. They may inform future work on host range dynamics, and they provide insights for understanding the emergence of pathogens that have human and veterinary health implications.

Salmonella enterica serovar Typhimurium strains from passerines have caused wild bird deaths and human salmonellosis outbreaks in Europe, Oceania, and North America. Here, we performed comparative genomic analysis to explore the emergence, genetic relationship, and evolution of geographically dispersed passerine isolates. We found that passerine isolates from Europe and the United States clustered to form two lineages (EU and US passerine lineages), which were distinct from major S. Typhimurium lineages circulating in other diverse hosts (e.g., humans, cattle, pigs, chickens, and other avian hosts, such as pigeons and ducks). Further, passerine isolates from New Zealand clustered to form a sublineage (NZ passerine lineage) of the US passerine lineage. We inferred that the passerine isolates mutated at a rate of 3.2 × 10(−7) substitutions/site/year, and the US, EU, and NZ passerine lineages emerged in approximately 1952, 1970, and 1996, respectively. Isolates from the three lineages presented genetic similarity, such as lack of antimicrobial resistance genes and accumulation of the same virulence pseudogenes. In addition, genetic diversity due to microevolution existed in the three passerine lineages. Specifically, pseudogenization in the type 1 fimbrial gene fimC (deletion of G at position 87) was detected only in the US and NZ passerine isolates, while single-base deletions in type 3 secretion system effector genes (i.e., gogB, sseJ, and sseK2) cooccurred solely in the EU passerine isolates. These findings provide insights into the evolution, host adaptation, and epidemiology of S. Typhimurium in passerines. IMPORTANCE Passerine-associated S. Typhimurium strains have been linked to human salmonellosis outbreaks in recent years. Here, we investigated the phylogenetic relationship of globally distributed passerine isolates and profiled their genomic similarity and diversity. Our study reveals two passerine-associated S. Typhimurium lineages circulating in Europe, Oceania, and North America. Isolates from the two lineages presented phylogenetic and genetic signatures that were distinct from those of isolates from other hosts. The findings shed light on the host adaptation of S. Typhimurium in passerines and are important for source attribution of S. Typhimurium strains to avian hosts. Further, we found that S. Typhimurium definitive phage type 160 (DT160) from passerines, which caused decades-long human salmonellosis outbreaks in New Zealand and Australia, formed a sublineage of the US passerine lineage, suggesting that DT160 might have originated from passerines outside Oceania. Our study demonstrates the importance of whole-genome sequencing and genomic analysis of historical microbial collections to modern epidemiologic surveillance.

Xylella fastidiosa is an insect-transmitted bacterial plant pathogen found across the Americas and, more recently, worldwide. X. fastidiosa infects plants of at least 563 species belonging to 82 botanical families. While the species X. fastidiosa infects many plants, particular strains have increased plant specificity. Understanding the molecular underpinnings of plant host specificity in X. fastidiosa is vital for predicting host shifts and epidemics. While there may exist multiple genetic determinants of host range in X. fastidiosa, the drivers of the unique relationships between X. fastidiosa and its hosts should be elucidated. Our objective with this study was to predict the ancestral plant hosts of this pathogen using phylogenetic and genomic methods based on a large data set of pathogen whole-genome data from agricultural hosts. We used genomic data to construct maximum-likelihood (ML) phylogenetic trees of subsets of the core and pan-genomes. With those trees, we ran ML ancestral state reconstructions of plant host at two taxonomic scales (genus and multiorder clades). Both the core and pan-genomes were informative in terms of predicting ancestral host state, giving new insight into the history of the plant hosts of X. fastidiosa. Subsequently, gene gain and loss in the pan-genome were found to be significantly correlated with plant host through genes that had statistically significant associations with particular hosts. IMPORTANCE Xylella fastidiosa is a globally important bacterial plant pathogen with many hosts; however, the underpinnings of host specificity are not known. This paper contains important findings about the usage of phylogenetics to understand the history of host specificity in this bacterial species, as well as convergent evolution in the pan-genome. There are strong signals of historical host range that give us insights into the history of this pathogen and its various invasions. The data from this paper are relevant in making decisions for quarantine and eradication, as they show the historical trends of host switching, which can help us predict likely future host shifts. We also demonstrate that using multilocus sequence type (MLST) genes in this system, which is still a commonly used process for policymaking, does not reconstruct the same phylogenetic topology as whole-genome data.

Bacterial porin-encoding genes are often found under positive selection. Local recombination has also been identified in a few of them to facilitate bacterial rapid adaptation, although it remains unknown whether it is a common evolutionary mechanism for the porins or outer membrane proteins in Gram-negative bacteria. In this study, we investigated the beta-barrel (β-barrel) porin-encoding genes in Escherichia coli that were reported under positive Darwinian selection. Besides fhuA that was found with ingenic local recombination previously, we identified four other genes, i.e., lamB, ompA, ompC, and ompF, all showing the similar mosaic evolution patterns. Comparative analysis of the protein sequences disclosed a list of highly variable regions in each family, which are mostly located in the convex of extracellular loops and coinciding with the binding sites of bacteriophages. For each of the porin families, mosaic recombination leads to unique combinations of the variable regions with different sequence patterns, generating diverse protein groups. Structural modeling indicated a conserved global topology among the different porins, with the extracellular surface varying a lot due to individual or combinatorial variable regions. The conservation of global tertiary structure would ensure the channel activity, while the wide diversity of variable regions may represent selection to avoid the invasion of phages, antibiotics or immune surveillance factors. Our study identified multiple bacterial porin genes with mosaic evolution. We hypothesize that this could be generalized strategy for outer membrane proteins to both maintain normal life processes and evade the attack of unfavored factors rapidly. IMPORTANCE Microevolution studies can disclose more elaborate evolutionary mechanisms of genes, appearing especially important for genes with multifaceted function such as those encoding outer membrane proteins. However, in most cases, the gene is considered as a whole unit, and the evolutionary patterns are disclosed. Here, we report that multiple bacterial porin proteins follow mosaic evolution, with local ingenic recombination combined with spontaneous mutations based on positive Darwinian selection, and conservation for most structural regions. This could represent a common mechanism for bacterial outer membrane proteins. The variable regions within each porin family showed large coincidence with the binding sites of bacteriophages, antibiotics, and immune factors and therefore would represent effective targets for the development of new antibacterial agents or vaccines.

Lyme borreliosis is the most common vector-borne disease in the Northern Hemisphere, caused by spirochetes belonging to the Borrelia burgdorferi sensu lato species complex, which are transmitted by ixodid ticks. B. burgdorferi sensu lato species produce a family of proteins on the linear plasmid 54 (PFam54), some of which confer the functions of cell adhesion and inactivation of complement, the first line of host defense. However, the impact of PFam54 in promoting B. burgdorferi sensu lato pathogenesis remains unclear because of the hurdles to simultaneously knock out all PFam54 proteins in a spirochete. Here, we describe two Borrelia bavariensis strains, PBN and PNi, isolated from patients naturally lacking PFam54 but maintaining the rest of the genome with greater than 95% identity to the reference B. bavariensis strain, PBi. We found that PBN and PNi less efficiently survive in human serum than PBi. Such defects were restored by introducing two B. bavariensis PFam54 recombinant proteins, BGA66 and BGA71, confirming the role of these proteins in providing complement evasion of B. bavariensis. Further, we found that all three strains remain detectable in various murine tissues 21 days post-subcutaneous infection, supporting the nonessential role of B. bavariensis PFam54 in promoting spirochete persistence. This study identified and utilized isolates deficient in PFam54 to associate the defects with the absence of these proteins, building the foundation to further study the role of each PFam54 protein in contributing to B. burgdorferi sensu lato pathogenesis. IMPORTANCE To establish infections, Lyme borreliae utilize various means to overcome the host’s immune system. Proteins encoded by the PFam54 gene array play a role in spirochete survival in vitro and in vivo. Moreover, this gene array has been described in all currently available Lyme borreliae genomes. By investigating the first two Borrelia bavariensis isolates naturally lacking the entire PFam54 gene array, we showed that both patient isolates display an increased susceptibility to human serum, which can be rescued in the presence of two PFam54 recombinant proteins. However, both isolates remain infectious to mice after intradermal inoculation, suggesting the nonessential role of PFam54 during the long-term, but may differ slightly in the colonization of specific tissues. Furthermore, these isolates show high genomic similarity to type strain PBi (>95%) and could be used in future studies investigating the role of each PFam54 protein in Lyme borreliosis pathogenesis.

Clostridium thermocellum is a thermophilic, anaerobic bacterium that natively ferments cellulose to ethanol and is a candidate for cellulosic biofuel production. Recently, we identified a hypermutator strain of C. thermocellum with a C669Y mutation in the polC gene, which encodes a DNA polymerase III enzyme. Here, we reintroduced this mutation using recently developed CRISPR tools to demonstrate that this mutation is sufficient to recreate the hypermutator phenotype. The resulting strain shows an approximately 30-fold increase in the mutation rate. This mutation is hypothesized to function by interfering with metal ion coordination in the PHP (polymerase and histidinol phosphatase) domain, which is responsible for proofreading. The ability to selectively increase the mutation rate in C. thermocellum is a useful tool for future directed evolution experiments. IMPORTANCE Cellulosic biofuels are a promising approach to decarbonize the heavy-duty-transportation sector. A longstanding barrier to cost-effective cellulosic biofuel production is the recalcitrance of cellulose to solubilization. Native cellulose-consuming organisms, such as Clostridium thermocellum, are promising candidates for cellulosic biofuel production; however, they often need to be genetically modified to improve product formation. One approach is adaptive laboratory evolution. Our findings demonstrate a way to increase the mutation rate in this industrially relevant organism, which can reduce the time needed for adaptive evolution experiments.

Horizontal gene transfer was long thought to be marginal in Mollicutes, but the capacity of some of these wall-less bacteria to exchange large chromosomal regions has been recently documented. Mycoplasma chromosomal transfer (MCT) is an unconventional mechanism that relies on the presence of a functional integrative conjugative element (ICE) in at least one partner and involves the horizontal acquisition of small and large chromosomal fragments from any part of the donor genome, which results in progenies composed of an infinite variety of mosaic genomes. The present study focuses on Mycoplasma bovis, an important pathogen of cattle responsible for major economic losses worldwide. By combining phylogenetic tree reconstructions and detailed comparative genome analyses of 36 isolates collected in Spain (2016 to 2018), we confirmed the mosaic nature of 16 field isolates and mapped chromosomal transfers exchanged between their hypothetical ancestors. This study provides evidence that MCT can take place in the field, most likely during coinfections by multiple strains. Because mobile genetic elements (MGEs) are classical contributors of genome plasticity, the presence of phages, insertion sequences (ISs), and ICEs was also investigated. Data revealed that these elements are widespread within the M. bovis species and evidenced classical horizontal transfer of phages and ICEs in addition to MCT. These events contribute to wide-genome diversity and reorganization within this species and may have a tremendous impact on diagnostic and disease control. IMPORTANCE Mycoplasma bovis is a major pathogen of cattle that has significant detrimental effects on economics and animal welfare in cattle rearing worldwide. Understanding the evolution and the adaptative potential of pathogenic mycoplasma species in the natural host is essential to combating them. In this study, we documented the occurrence of mycoplasma chromosomal transfer, an atypical mechanism of horizontal gene transfer, in field isolates of M. bovis that provide new insights into the evolution of this pathogenic species in their natural host. Although these events are expected to occur at low frequency, their impact is accountable for genome-wide variety and reorganization within M. bovis species, which may compromise both diagnostic and disease control.

Escherichia coli is a commensal species of the lower intestine, but also a major pathogen causing intestinal and extra-intestinal infections, increasingly prevalent and resistant to antibiotics. Most studies on genomic evolution of E. coli used isolates from infections. Here instead, we whole-genome sequenced a collection of 403 commensal E. coli isolated from fecal samples of healthy adult volunteers in France (1980-2010). These isolates were distributed mainly in phylogroups A and B2 (30% each) and belonged to 152 sequence types (STs), the five most frequent being ST10 (phylogroup A) (16.3%), ST73 and ST95 (phylogroup B2) (6.3 and 5.0%, respectively), ST69 (phylogroup D) (4.2%) and ST59 (phylogroup F) (3.9%), and 224 O:H serotypes. ST and serotype diversity increased over time. The O1, O2, O6 and O25-groups used in bioconjugate O-antigen vaccine against extra-intestinal infections were found in 23% of the strains of our collection. The increase in frequency of virulence-associated genes and antibiotic resistances was driven by two evolutionary mechanisms. Evolution of virulence gene frequency was driven by both clonal expansion of STs with more virulence genes (“ST-driven”) and increases in gene frequency within STs independently of changes in ST frequencies (“gene-driven”). In contrast, the evolution of resistance was dominated by increases in frequency within STs (“gene-driven”). This study provides a unique picture of the phylogenomic evolution of E. coli in its human commensal habitat over 30 years and will have implications for the development of preventive strategies.IMPORTANCEEscherichia coli is an opportunistic pathogen with the greatest burden of antibiotic resistance, one of the main causes of bacterial infections and an increasing concern in an ageing population. Deciphering the evolutionary dynamics of virulence and antibiotic resistance in commensal E. coli is important to understand adaptation and anticipate future changes. The gut of vertebrates is the primary habitat of E. coli and probably where selection for virulence and resistance take place. Unfortunately, most whole-genome sequenced strains are isolated from pathogenic conditions. Here, we whole genome sequenced 403 E. coli commensals isolated from healthy French subjects on a 30-year period. Virulence genes increased in frequency by both clonal expansion of clones carrying them and increases in frequency within clones whereas resistance genes increased by within clone increased frequency. Prospective studies of E. coli commensals should be performed worldwide to have a broader picture of evolution and adaptation of this species.

Tailed bacteriophages constitute the bulk of the intestinal viromes of the vertebrate animals. However, the relationships between lytic and lysogenic lifestyles of the phages in these ecosystems are not always clear and may vary between the species or even between the individuals. The human intestinal (fecal) viromes are believed to be dominated by temperate phages, while in the horse feces the virulent phages are more prevalent. Almost all the isolates of horse fecal coliphages are virulent. Phage Hf4s is the first temperate equine intestinal coliphage characterized. It was isolated from the horse feces on the indigenous equine E. coli 4s strain. It is a podovirus, related to Lederbergvirus genus (including the well–characterized Salmonella phage P22). Hf4s recognizes the host O antigen as its primary receptor and possesses a functional O-antigen seroconversion cluster that renders the lysogens protected from the superinfection by the same phage and also abolishes the adsorption of some indigenous equine virulent coliphages, such as DT57C, while the other phages, such as G7C or phiKT retain the ability to infect E. coli 4s (Hf4s) lysogens.ImportanceThe relationships between virulent and temperate bacteriophages and their impact on high-density symbiotic microbial ecosystems of the animal are not always clear and may vary between the species or even between the individuals. The horse intestinal virome is dominated by the virulent phages, and Hf4s is the first temperate equine intestinal coliphage characterized. It recognizes the host O antigen as its primary receptor and possesses a functional O-antigen seroconversion cluster that renders the lysogens protected from the superinfection by some indigenous equine virulent coliphages, such as DT57C, while the other phages, such as G7C or phiKT retain the ability to infect E. coli 4s (Hf4s) lysogens. microbial viruses in the mammal intestinal ecosystems.

Bacteria rapidly adapt to their environment by integrating external stimuli through diverse signal transduction systems. Pseudomonas aeruginosa, for example, senses surface contact through the Wsp signal transduction system to trigger the production of cyclic di-GMP. Diverse mutations in wsp genes that manifest enhanced biofilm formation are frequently reported in clinical isolates of P. aeruginosa and in biofilm studies of Pseudomonas spp. and Burkholderia cenocepacia. In contrast to the convergent phenotypes associated with comparable wsp mutations, we demonstrate that the Wsp system in B. cenocepacia does not impact intracellular cyclic di-GMP levels, unlike that in Pseudomonas spp. Our current mechanistic understanding of the Wsp system is based entirely on the study of four Pseudomonas spp., and its phylogenetic distribution remains unknown. Here, we present a broad phylogenetic analysis to show that the Wsp system originated in the betaproteobacteria and then horizontally transferred to Pseudomonas spp., the sole member of the gammaproteobacteria. Alignment of 794 independent Wsp systems with reported mutations from the literature identified key amino acid residues that fall within and outside annotated functional domains. Specific residues that are highly conserved but uniquely modified in B. cenocepacia likely define mechanistic differences among Wsp systems. We also find the greatest sequence variation in the extracellular sensory domain of WspA, indicating potential adaptations to diverse external stimuli beyond surface contact sensing. This study emphasizes the need to better understand the breadth of functional diversity of the Wsp system as a major regulator of bacterial adaptation beyond B. cenocepacia and select Pseudomonas spp. IMPORTANCE The Wsp signal transduction system serves as an important model system for studying how bacteria adapt to living in densely structured communities known as biofilms. Biofilms frequently cause chronic infections and environmental fouling, and they are very difficult to eradicate. In Pseudomonas aeruginosa, the Wsp system senses contact with a surface, which in turn activates specific genes that promote biofilm formation. We demonstrate that the Wsp system in Burkholderia cenocepacia regulates biofilm formation uniquely from that in Pseudomonas species. Furthermore, a broad phylogenetic analysis reveals the presence of the Wsp system in diverse bacterial species, and sequence analyses of 794 independent systems suggest that the core signaling components function similarly but with key differences that may alter what or how they sense. This study shows that Wsp systems are highly conserved and more broadly distributed than previously thought, and their unique differences likely reflect adaptations to distinct environments.

Phaeobacter inhibens 2.10 is an effective biofilm former on marine surfaces and has the ability to outcompete other microorganisms, possibly due to the production of the plasmid-encoded secondary metabolite tropodithietic acid (TDA). P. inhibens 2.10 biofilms produce phenotypic variants with reduced competitiveness compared to the wild type. In the present study, we used longitudinal, genome-wide deep sequencing to uncover the genetic foundation that contributes to the emergent phenotypic diversity in P. inhibens 2.10 biofilm dispersants. Our results show that phenotypic variation is not due to the loss of the plasmid that carries the genes for TDA synthesis but instead show that P. inhibens 2.10 biofilm populations become rapidly enriched in single nucleotide variations in genes involved in the synthesis of TDA. While variants in genes previously linked to other phenotypes, such as lipopolysaccharide production (i.e., rfbA) and cellular persistence (i.e., metG), also appear to be selected for during biofilm dispersal, the number and consistency of variations found for genes involved in TDA production suggest that this metabolite imposes a burden on P. inhibens 2.10 cells. Our results indicate a strong selection pressure for the loss of TDA in monospecies biofilm populations and provide insight into how competition (or a lack thereof) in biofilms might shape genome evolution in bacteria. IMPORTANCE Biofilm formation and dispersal are important survival strategies for environmental bacteria. During biofilm dispersal, cells often display stable and heritable variants from the parental biofilm. Phaeobacter inhibens is an effective colonizer of marine surfaces, in which a subpopulation of its biofilm dispersal cells displays a noncompetitive phenotype. This study aimed to elucidate the genetic basis of these phenotypic changes. Despite the progress made to date in characterizing the dispersal variants in P. inhibens, little is understood about the underlying genetic changes that result in the development of the specific variants. Here, P. inhibens phenotypic variation was linked to single nucleotide polymorphisms (SNPs), in particular in genes affecting the competitive ability of P. inhibens, including genes related to the production of the antibiotic tropodithietic acid (TDA) and bacterial cell-cell communication (e.g., quorum sensing). This work is significant as it reveals how the biofilm lifestyle might shape genome evolution in a cosmopolitan bacterium.

The concept of bacterial dark matter stems from our inability to culture most microbes and represents a fundamental gap in our knowledge of microbial diversity. Here, we present the domestication of such an organism: a previously uncultured, novel species from the rare Actinomycetes genus Verrucosispora. Although initial recovery took >4 months, isolation of phenotypically distinct, domesticated generations occurred within weeks. Two isolates were subjected to phenogenomic analyses, revealing domestication correlated with enhanced growth rates in nutrient-rich media but diminished capacity to metabolize diverse amino acids. This is seemingly mediated by genomic atrophy through a mixed approach of pseudogenization and reversion of pseudogenization of amino acid metabolism genes. Conversely, later generational strains had enhanced spore germination rates, potentially through the reversion of a sporulation-associated kinase from pseudogene to true gene status. We observed that our most wild-type isolate had the greatest potential for antibacterial activity, which correlated with extensive mutational attrition of biosynthetic gene clusters in domesticated strains. Comparative analyses revealed wholesale genomic reordering in strains, with widespread single nucleotide polymorphism, indel, and pseudogene-impactful mutations observed. We hypothesize that domestication of this previously unculturable organism resulted from the shedding of genomic flexibility required for life in a dynamic marine environment, parsing out genetic redundancy to allow for a newfound cultivable amenability. IMPORTANCE The majority of environmental bacteria cannot be cultured within the laboratory. Understanding why only certain environmental isolates can be recovered is key to unlocking the abundant microbial dark matter that is widespread on our planet. In this study, we present not only the culturing but domestication of just such an organism. Although initial recovery took >4 months, we were able to isolate distinct, subpassaged offspring from the originating colony within mere weeks. A phenotypic and genotypic analysis of our generational strains revealed that adaptation to life in the lab occurred as a result of wholesale mutational changes. These permitted an enhanced ability for growth in nutrient rich media but came at the expense of reduced genomic flexibility. We suggest that without dynamic natural environmental stressors our domesticated strains effectively underwent genomic atrophy as they adapted to static conditions experienced in the laboratory.

Cell-cell adhesion between oral bacteria plays a key role in the development of polymicrobial communities such as dental plaque. Oral streptococci such as Streptococcus gordonii and Streptococcus oralis are important early colonizers of dental plaque and bind to a wide range of different oral microorganisms, forming multispecies clumps or “coaggregates.” S. gordonii actively responds to coaggregation by regulating gene expression. To further understand these responses, we assessed gene regulation in S. gordonii and S. oralis following coaggregation in 25% human saliva. Coaggregates were formed by mixing, and after 30 min, RNA was extracted for dual transcriptome sequencing (RNA-Seq) analysis. In S. oralis, 18 genes (6 upregulated and 12 downregulated) were regulated by coaggregation. Significantly downregulated genes encoded functions such as amino acid and antibiotic biosynthesis, ribosome, and central carbon metabolism. In total, 28 genes were differentially regulated in Streptococcus gordonii (25 upregulated and 3 downregulated). Many genes associated with transporters and a two-component (NisK/SpaK) regulatory system were upregulated following coaggregation. Our comparative analyses of S. gordonii-S. oralis with different previously published S. gordonii pairings (S. gordonii-Fusobacterium nucleatum and S. gordonii-Veillonella parvula) suggest that the gene regulation is specific to each pairing, and responses do not appear to be conserved. This ability to distinguish between neighboring bacteria may be important for S. gordonii to adapt appropriately during the development of complex biofilms such as dental plaque. IMPORTANCE Dental plaque is responsible for two of the most prevalent diseases in humans, dental caries and periodontitis. Controlling the formation of dental plaque and preventing the transition from oral health to disease requires a detailed understanding of microbial colonization and biofilm development. Streptococci are among the most common colonizers of dental plaque. This study identifies key genes that are regulated when oral streptococci bind to one another, as they do in the early stages of dental plaque formation. We show that specific genes are regulated in two different oral streptococci following the formation of mixed-species aggregates. The specific responses of S. gordonii to coaggregation with S. oralis are different from those to coaggregation with other oral bacteria. Targeting the key genes that are upregulated during interspecies interactions may be a powerful approach to control the development of biofilm and maintain oral health.

Salmonella enterica serovar Heidelberg is isolated from poultry-producing regions around the world. In Brazil, S. Heidelberg has been frequently detected in poultry flocks, slaughterhouses, and chicken meat. The goal of the present study was to assess the population structure, recent temporal evolution, and some important genetic characteristics of S. Heidelberg isolated from Brazilian poultry farms. Phylogenetic analysis of 68 S. Heidelberg genomes sequenced here and additional whole-genome data from NCBI demonstrated that all isolates from the Brazilian poultry production chain clustered into a monophyletic group, here called S. Heidelberg Brazilian poultry lineage (SH-BPL). Bayesian analysis defined the time of the most recent common ancestor (tMRCA) as 2004, and the overall population size (N(e)) was constant until 2008, when an ∼10-fold N(e) increase was observed until circa 2013. SH-BPL presented at least two plasmids with replicons ColpVC (n = 68; 100%), IncX1 (n = 66; 97%), IncA/C2 (n = 65; 95.5%), ColRNAI (n = 43; 63.2%), IncI1 (n = 32; 47%), ColMG828, Col156, IncHI2A, IncHI2, IncQ1, IncX4, IncY, and TrfA (each with n

Multidrug-resistant (MDR) Escherichia coli strains that carry extended-spectrum β-lactamases (ESBLs) or colistin resistance gene mcr-1 have been identified in the human gut at an increasing incidence worldwide. In this study, we isolated and characterized MDR Enterobacteriaceae from the gut microbiota of healthy Singaporeans and show that the detection rates for ESBL-producing and mcr-positive Enterobacteriaceae are 25.7% (28/109) and 7.3% (8/109), respectively. Whole-genome sequencing analysis of the 37 E. coli isolates assigned them into 25 sequence types and 6 different phylogroups, suggesting that the MDR E. coli gut colonizers are highly diverse. We then analyzed the genetic context of the resistance genes and found that composite transposons played important roles in the cotransfer of bla(CTX-M-15/55) and qnrS1, as well as the acquisition of mcr-1. Furthermore, comparative genomic analysis showed that 12 of the 37 MDR E. coli isolates showed high similarity to ESBL-producing E. coli isolates from raw meat products in local markets. By analyzing the core genome single nucleotide polymorphisms (SNPs) shared by these isolates, we identified possible clonal transmission of an MDR E. coli clone between human and raw meat, as well as a group of highly similar IncI2 (Delta) plasmids that might be responsible for the dissemination of mcr-1 in a much wider geographic region. Together, these results suggest that antibiotic resistance may be transmitted between different environmental settings by the expansion of MDR E. coli clones, as well as by the dissemination of resistance plasmids. IMPORTANCE The human gut can harbor both antibiotic-resistant and virulent Escherichia coli which may subsequently cause infections. In this study, we found that multidrug-resistant (MDR) E. coli isolates from the gut of healthy Singaporeans carry a diverse range of antibiotic resistance mechanisms and virulence factor genes and are highly diverse. By comparing their genomes with the extended-spectrum β-lactamase (ESBL)-producing E. coli isolates from raw meat products that were sampled at a similar time from local markets, we detected an MDR E. coli clone that was possibly transmitted between humans and raw meat products. Furthermore, we also found that a group of resistance plasmids might be responsible for the dissemination of colistin resistance gene mcr-1 in Singapore, Malaysia, and Europe. Our findings call for better countermeasures to block the transmission of antibiotic resistance.

In Serratia marcescens JNB5-1, prodigiosin was highly produced at 30°C, but it was noticeably repressed at ≥37°C. Our initial results demonstrated that both the production and the stability of the O-methyl transferase (PigF) and oxidoreductase (PigN) involved in the prodigiosin pathway in S. marcescens JNB5-1 sharply decreased at ≥37°C. Therefore, in this study, we improved mRNA stability and protein production using de novo polynucleotide fragments (PNFs) and the introduction of disulfide bonds, respectively, and observed their effects on prodigiosin production. Our results demonstrate that adding PNFs at the 3′ untranslated regions of pigF and pigN significantly improved the mRNA half-lives of these genes, leading to an increase in the transcript and expression levels. Subsequently, the introduction of disulfide bonds in pigF improved the thermal stability, pH stability, and copper ion resistance of PigF. Finally, shake flask fermentation showed that the prodigiosin titer with the engineered S. marcescens was increased by 61.38% from 5.36 to 8.65 g/liter compared to the JNB5-1 strain at 30°C and, significantly, the prodigiosin yield increased 2.05-fold from 0.38 to 0.78 g/liter at 37°C. In this study, we revealed that the introduction of PNFs and disulfide bonds greatly improved the expression and stability of pigF and pigN, hence efficiently enhancing prodigiosin production with S. marcescens at 30 and 37°C. IMPORTANCE This study highlights a promising strategy to improve mRNA/enzyme stability and to increase production using de novo PNF libraries and the introduction of disulfide bonds into the protein. PNFs could increase the half-life of target gene mRNA and effectively prevent its degradation. Moreover, PNFs could increase the relative intensity of target genes without affecting the expression of other genes; as a result, it could alleviate the cellular burden compared to other regulatory elements such as promoters. In addition, we obtained a PigF variant with improved activity and stability by the introduction of disulfide bonds into PigF. Collectively, we demonstrate here a novel approach for improving mRNA/enzyme stability using PNFs, which results in enhanced prodigiosin production in S. marcescens at 30°C.

Approaches for recovering and analyzing genomes belonging to novel, hitherto-unexplored bacterial lineages have provided invaluable insights into the metabolic capabilities and ecological roles of yet-uncultured taxa. The phylum Acidobacteria is one of the most prevalent and ecologically successful lineages on Earth, yet currently, multiple lineages within this phylum remain unexplored. Here, we utilize genomes recovered from Zodletone Spring, an anaerobic sulfide and sulfur-rich spring in southwestern Oklahoma, as well as from multiple disparate soil and nonsoil habitats, to examine the metabolic capabilities and ecological role of members of family UBA6911 (group 18) Acidobacteria. The analyzed genomes clustered into five distinct genera, with genera Gp18_AA60 and QHZH01 recovered from soils, genus Ga0209509 from anaerobic digestors, and genera Ga0212092 and UBA6911 from freshwater habitats. All genomes analyzed suggested that members of Acidobacteria group 18 are metabolically versatile heterotrophs capable of utilizing a wide range of proteins, amino acids, and sugars as carbon sources, possess respiratory and fermentative capacities, and display few auxotrophies. Soil-dwelling genera were characterized by larger genome sizes, higher numbers of CRISPR loci, an expanded carbohydrate active enzyme (CAZyme) machinery enabling debranching of specific sugars from polymers, possession of a C(1) (methanol and methylamine) degradation machinery, and a sole dependence on aerobic respiration. In contrast, nonsoil genomes encoded a more versatile respiratory capacity for oxygen, nitrite, sulfate, and trimethylamine N-oxide (TMAO) respiration, as well as the potential for utilizing the Wood-Ljungdahl (WL) pathway as an electron sink during heterotrophic growth. Our results not only expand our knowledge of the metabolism of a yet-uncultured bacterial lineage but also provide interesting clues on how terrestrialization and niche adaptation drive metabolic specialization within the Acidobacteria. IMPORTANCE Members of the Acidobacteria are important players in global biogeochemical cycles, especially in soils. A wide range of acidobacterial lineages remain currently unexplored. We present a detailed genomic characterization of genomes belonging to family UBA6911 (also known as group 18) within the phylum Acidobacteria. The genomes belong to different genera and were obtained from soil (genera Gp18_AA60 and QHZH01), freshwater habitats (genera Ga0212092 and UBA6911), and an anaerobic digestor (genus Ga0209509). While all members of the family shared common metabolic features, e.g., heterotrophic respiratory abilities, broad substrate utilization capacities, and few auxotrophies, distinct differences between soil and nonsoil genera were observed. Soil genera were characterized by expanded genomes, higher numbers of CRISPR loci, a larger carbohydrate active enzyme (CAZyme) repertoire enabling monomer extractions from polymer side chains, and methylotrophic (methanol and methylamine) degradation capacities. In contrast, nonsoil genera encoded more versatile respiratory capacities for utilizing nitrite, sulfate, TMAO, and the WL pathway, in addition to oxygen as electron acceptors. Our results not only broaden our understanding of the metabolic capacities within the Acidobacteria but also provide interesting clues on how terrestrialization shaped Acidobacteria evolution and niche adaptation.

Perturbation of natural microbial communities by antimicrobials, such as triclosan, can result in selection for antibiotic tolerance, which is of particular concern when pathogens are present. Members of the genus Pseudomonas are found in many natural microbial communities and frequently demonstrate increased abundance following triclosan exposure. The pathogen and well-studied model organism Pseudomonas aeruginosa exhibits high triclosan tolerance; however, it is unknown if all Pseudomonas species share this trait or if there are susceptible strains. We characterized the triclosan tolerance phenotypes of diverse Pseudomonas isolates obtained from triclosan-exposed built environments and identified both tolerant and sensitive strains. High tolerance is associated with carriage of the enoyl-acyl carrier reductase (ENR) isozyme gene fabV, compared to the lesser protective effects of efflux or presence of ENRs. Given its unique importance, we examined fabV distribution throughout Pseudomonas species using large-scale phylogenomic analyses. We find fabV presence or absence is largely invariant at the species level but demonstrates multiple gain and loss events in its evolutionary history. We further provide evidence of its presence on mobile genetic elements. Our results demonstrate the surprising variability in triclosan tolerance in Pseudomonas and confirm fabV to be a useful indicator for high triclosan tolerance in Pseudomonas. These findings provide a framework for better monitoring of Pseudomonas in triclosan-exposed environments and interpreting effects on species and gene composition. IMPORTANCE Closely related species are typically assumed to demonstrate similar phenotypes driven by underlying conserved genotypes. When monitoring for the effect of antimicrobials on the types of species that may be selected for, this assumption may prove to be incorrect, and identification of additional genetic markers may be necessary. We isolated several phylogenetically diverse members of Pseudomonas from indoor environments and tested their phenotypic tolerance toward the commonly used antimicrobial triclosan. Although Pseudomonas isolates are broadly regarded to be highly triclosan tolerant, we demonstrate the presence of both triclosan-tolerant and -susceptible strains, separated by a difference in tolerance of nearly 3 orders of magnitude. Bioinformatic and experimental investigation demonstrated that the presence of the gene fabV was associated with high tolerance. We demonstrate that fabV is not evenly distributed in all Pseudomonas species and that its presence could be a useful predictor of high triclosan tolerance suitable for antimicrobial monitoring efforts involving triclosan.

Mycobacterium avium subsp. paratuberculosis (MAP) is the causative agent of Johne's disease in ruminants, which has important health consequences for dairy cattle. The Regional Dairy Quality Management Alliance (RDQMA) project is a multistate research program involving MAP isolates taken from three intensively studied commercial dairy farms in the northeastern United States, which emphasized longitudinal data collection of both MAP isolates and animal health in three regional dairy herds for a period of about 7 years. This paper reports the results of a pan-GWAS analysis involving 318 MAP isolates and dairy cow Johne’s disease phenotypes, taken from these three farms. Based on our highly curated accessory gene count, the pan-GWAS analysis identified several MAP genes associated with bovine Johne’s disease phenotypes scored from these three farms, with some of the genes having functions suggestive of possible cause/effect relationships with these phenotypes. This paper reports a pangenomic comparative analysis between MAP and Mycobacterium tuberculosis, assessing functional Gene Ontology category enrichments between these taxa. Finally, we also provide a population genomic perspective on the effectiveness of herd isolation, involving closed dairy farms, in preventing MAP interfarm cross infection on a microgeographic scale. IMPORTANCE Mycobacterium avium subsp. paratuberculosis (MAP) is the causative agent of Johne's disease in ruminants, which has important health consequences for dairy cattle and enormous economic consequences for the dairy industry. Understanding which genes in this bacterium are correlated with key disease phenotypes can lead to functional experiments targeting these genes and ultimately lead to improved control strategies. This study represents a rare example of a prolonged longitudinal study of dairy cattle where the disease was measured and the bacteria were isolated from the same cows. The genome sequences of over 300 MAP isolates were analyzed for genes that were correlated with a wide range of Johne’s disease phenotypes. A number of genes were identified that were significantly associated with several aspects of the disease and suggestive of further experimental follow-up.

This study showed the genetic diversity and population structure of S. aureus presenting the same spa type, t899, but belonging to different STs. Our findings revealed that these isolates vary deeply in their core and accessory genomes, contrary to what is regularly inferred from studies using spa typing only., Methicillin-resistant Staphylococcus aureus (MRSA) presenting spa type t899 is commonly associated with sequence type 9 (ST9) but is also increasingly linked to ST398. This study provides genomic insight into the diversity of t899 isolates using core genome multilocus sequence typing (cgMLST), single nucleotide polymorphism (SNP)-based phylogeny, and the description of selected antimicrobial resistance and virulence markers. The SNP-based phylogenic tree showed that isolates sharing the same spa type (t899) but different STs highly diverged in their core and accessory genomes, revealing discriminant antimicrobial resistance (AMR) and virulence markers. Our results highlighted the idea that in a surveillance context where only spa typing is used, an additional multiplex PCR for the detection of the tet(M), sak, and seg genes would be valuable in helping distinguish ST9 from ST398 isolates on a routine basis. IMPORTANCE This study showed the genetic diversity and population structure of S. aureus presenting the same spa type, t899, but belonging to different STs. Our findings revealed that these isolates vary deeply in their core and accessory genomes, contrary to what is regularly inferred from studies using spa typing only. Given that identical spa types can be associated with different STs and that spa typing only is not appropriate for S. aureus isolates that have undergone major recombination events which include the passage of the spa gene (such as in t899-positive MRSA), the combination of both MLST and spa typing methods is recommended. However, spa typing alone is still largely used in surveillance studies and basic characterization. Our data suggest that additional markers, such as tet(M), sak, and seg genes, could be implemented in an easy and inexpensive manner in order to identify S. aureus lineages with a higher accuracy.

Genomic data reveal single-nucleotide polymorphisms (SNPs) that may carry information about the evolutionary history of bacteria. However, it remains unclear what inferences about selection can be made from genomic SNP data. Bacterial species are often sampled during epidemic outbreaks or within hosts during the course of chronic infections. SNPs obtained from genomic analysis of these data are not necessarily fixed. Treating them as fixed during analysis by using measures such as the ratio of nonsynonymous to synonymous evolutionary changes (dN/dS) may lead to incorrect inferences about the strength and direction of selection. In this study, we consider data from a range of whole-genome sequencing studies of bacterial pathogens and explore patterns of nonsynonymous variation to assess whether evidence of selection can be identified by investigating SNP counts alone across multiple WGS studies. We visualize these SNP data in ways that highlight their relationship to neutral baseline expectations. These neutral expectations are based on a simple model of mutation, from which we simulate SNP accumulation to investigate how SNP counts are distributed under alternative assumptions about positive and negative selection. We compare these patterns with empirical SNP data and illustrate the general difficulty of detecting positive selection from SNP data. Finally, we consider whether SNP counts observed at the between-host population level differ from those observed at the within-host level and find some evidence that suggests that dynamics across these two scales are driven by different underlying processes. IMPORTANCE Identifying selection from SNP data obtained from whole-genome sequencing studies is challenging. Some current measures used to identify and quantify selection acting on genomes rely on fixed differences; thus, these are inappropriate for SNP data where variants are not fixed. With the increase in whole-genome sequencing studies, it is important to consider SNP data in the context of evolutionary processes. How SNPs are counted and analyzed can help in understanding mutation accumulation and trajectories of strains. We developed a tool for identifying possible evidence of selection and for comparative analysis with other SNP data. We propose a model that provides a rule-of-thumb guideline and two new visualization techniques that can be used to interpret and compare SNP data. We quantify the expected proportion of nonsynonymous SNPs in coding regions under neutrality and demonstrate its use in identifying evidence of positive and negative selection from simulations and empirical data.

Bacteriophages (phages) are currently available for use by the food industry to control the foodborne pathogen Listeria monocytogenes. Although phage biocontrols are effective under specific conditions, their use can select for phage-resistant bacteria that repopulate phage-treated environments. Here, we performed short-term coevolution experiments to investigate the impact of single phages and a two-phage cocktail on the regrowth of phage-resistant L. monocytogenes and the adaptation of the phages to overcome this resistance. We used whole-genome sequencing to identify mutations in the target host that confer phage resistance and in the phages that alter host range. We found that infections with Listeria phages LP-048, LP-125, or a combination of both select for different populations of phage-resistant L. monocytogenes bacteria with different regrowth times. Phages isolated from the end of the coevolution experiments were found to have gained the ability to infect phage-resistant mutants of L. monocytogenes and L. monocytogenes strains previously found to be broadly resistant to phage infection. Phages isolated from coinfected cultures were identified as recombinants of LP-048 and LP-125. Interestingly, recombination events occurred twice independently in a locus encoding two proteins putatively involved in DNA binding. We show that short-term coevolution of phages and their hosts can be utilized to obtain mutant and recombinant phages with adapted host ranges. These laboratory-evolved phages may be useful for limiting the emergence of phage resistance and for targeting strains that show general resistance to wild-type (WT) phages. IMPORTANCEListeria monocytogenes is a life-threatening bacterial foodborne pathogen that can persist in food processing facilities for years. Phages can be used to control L. monocytogenes in food production, but phage-resistant bacterial subpopulations can regrow in phage-treated environments. Coevolution experiments were conducted on a Listeria phage-host system to provide insight into the genetic variation that emerges in both the phage and bacterial host under reciprocal selective pressure. As expected, mutations were identified in both phage and host, but additionally, recombination events were shown to have repeatedly occurred between closely related phages that coinfected L. monocytogenes. This study demonstrates that in vitro evolution of phages can be utilized to expand the host range and improve the long-term efficacy of phage-based control of L. monocytogenes. This approach may also be applied to other phage-host systems for applications in biocontrol, detection, and phage therapy.

Gardnerella vaginalis is described as a common anaerobic vaginal bacterium whose presence may correlate with vaginal dysbiotic conditions. In the current study, we performed phylogenomic analyses of 72 G. vaginalis genome sequences, revealing noteworthy genome differences underlying a polyphyletic organization of this taxon. Particularly, the genomic survey revealed that this species may actually include nine distinct genotypes (GGtype1 to GGtype9). Furthermore, the observed link between sialidase and phylogenomic grouping provided clues of a connection between virulence potential and the evolutionary history of this microbial taxon. Specifically, based on the outcomes of these in silico analyses, GGtype3, GGtype7, GGtype8, and GGtype9 appear to have virulence potential since they exhibited the sialidase gene in their genomes. Notably, the analysis of 34 publicly available vaginal metagenomic samples allowed us to trace the distribution of the nine G. vaginalis genotypes identified in this study among the human population, highlighting how differences in genetic makeup could be related to specific ecological properties. Furthermore, comparative genomic analyses provided details about the G. vaginalis pan- and core genome contents, including putative genetic elements involved in the adaptation to the ecological niche as well as many putative virulence factors. Among these putative virulence factors, particularly noteworthy genes identified were the gene encoding cholesterol-dependent cytolysin (CDC) toxin vaginolysin and genes related to microbial biofilm formation, iron uptake, adhesion to the vaginal epithelium, as well as macrolide antibiotic resistance. IMPORTANCE The identification of nine different genotypes among members of G. vaginalis allowed us to distinguish an uneven distribution of virulence-associated genetic traits within this taxon and thus suggest the potential occurrence of putative pathogen and commensal G. vaginalis strains. These findings, coupled with metagenomics microbial profiling of human vaginal microbiota, permitted us to get insights into the distribution of the genotypes among the human population, highlighting the presence of different structural communities in terms of G. vaginalis genotypes.

Within the landscape of globally distributed pathogens, populations differentiate via both adaptive and nonadaptive forces. Individual populations are likely to show unique trends of genetic diversity, host-pathogen interaction, and ecological adaptation. In plant pathogens, allopatric divergence may occur particularly rapidly within simplified agricultural monoculture landscapes. As such, the study of plant pathogen populations in monocultures can highlight the distinct evolutionary mechanisms that lead to local genetic differentiation. Xylella fastidiosa is a plant pathogen known to infect and damage multiple monocultures worldwide. One subspecies, Xylella fastidiosa subsp. fastidiosa, was first introduced to the United States ∼150 years ago, where it was found to infect and cause disease in grapevines (Pierce’s disease of grapevines, or PD). Here, we studied PD-causing subsp. fastidiosa populations, with an emphasis on those found in the United States. Our study shows that following their establishment in the United States, PD-causing strains likely split into populations on the East and West Coasts. This diversification has occurred via both changes in gene content (gene gain/loss events) and variations in nucleotide sequence (mutation and recombination). In addition, we reinforce the notion that PD-causing populations within the United States acted as the source for subsequent subsp. fastidiosa outbreaks in Europe and Asia. IMPORTANCE Compared to natural environments, the reduced diversity of monoculture agricultural landscapes can lead bacterial plant pathogens to quickly adapt to local biological and ecological conditions. Because of this, accidental introductions of microbial pathogens into naive regions represents a significant economic and environmental threat. Xylella fastidiosa is a plant pathogen with an expanding host and geographic range due to multiple intra- and intercontinental introductions. X. fastidiosa subsp. fastidiosa infects and causes disease in grapevines (Pierce’s disease of grapevines [PD]). This study focused on PD-causing X. fastidiosa populations, particularly those found in the United States but also invasions into Taiwan and Spain. The analysis shows that PD-causing X. fastidiosa has diversified via multiple cooccurring evolutionary forces acting at an intra- and interpopulation level. This analysis enables a better understanding of the mechanisms leading to the local adaptation of X. fastidiosa and how a plant pathogen diverges allopatrically after multiple and sequential introduction events.

This research offers a new route for a wider understanding of the dependency between pathogenicity and phylogeny of a natural bacterial population, specifically within Bacillus cereus sensu lato, that is widely distributed around the world and easily transferred into food products. Our study indicates differences in the phylogenetic and geographical distributions of potential enterotoxigenic B. cereus sensu lato strains. Hence, these bacilli possess a risk for human health, and rapid testing methods for their identification are greatly needed. In particular, the detection of the CytK enterotoxin should be a supporting strategy for the identification of pathogenic B. cereus sensu lato., Bacillus cereus sensu lato comprises Gram-positive spore-forming bacteria producing toxins associated with foodborne diseases. Three pore-forming enterotoxins, nonhemolytic enterotoxin (Nhe), hemolysin BL (Hbl), and cytotoxin K (CytK), are considered the primary factors in B. cereus sensu lato diarrhea. The aim of this study was to determine the potential risk of enterotoxicity among soil B. cereus sensu lato isolates representing diverse phylogroups and originated from different geographic locations with various climates (Burkina Faso, Kenya, Argentina, Kazakhstan, and Poland). While nheA- and hblA-positive isolates were present among all B. cereus sensu lato populations and distributed across all phylogenetic groups, cytK-2-positive strains predominated in geographic regions with an arid hot climate (Africa) and clustered together on a phylogenetic tree mainly within mesophilic groups III and IV. The highest in vitro cytotoxicity to Caco-2 and HeLa cells was demonstrated by the strains clustered within phylogroups II and IV. Overall, our results suggest that B. cereus sensu lato pathogenicity is a comprehensive process conditioned by many intracellular factors and diverse environmental conditions. IMPORTANCE This research offers a new route for a wider understanding of the dependency between pathogenicity and phylogeny of a natural bacterial population, specifically within Bacillus cereus sensu lato, that is widely distributed around the world and easily transferred into food products. Our study indicates differences in the phylogenetic and geographical distributions of potential enterotoxigenic B. cereus sensu lato strains. Hence, these bacilli possess a risk for human health, and rapid testing methods for their identification are greatly needed. In particular, the detection of the CytK enterotoxin should be a supporting strategy for the identification of pathogenic B. cereus sensu lato.

Cerato-platanins (CPs) form a family of fungal small secreted cysteine-rich proteins (SSCPs) and are of particular interest not only because of their surface activity but also their abundant secretion by fungi. We performed an evolutionary analysis of 283 CPs from 157 fungal genomes with the focus on the environmental opportunistic plant-beneficial and mycoparasitic fungus Trichoderma. Our results revealed a long evolutionary history of CPs in Dikarya fungi that have undergone several events of lateral gene transfer and gene duplication. Three genes were maintained in the core genome of Trichoderma, while some species have up to four CP-encoding genes. All Trichoderma CPs evolve under stabilizing natural selection pressure. The functional genomic analysis of CPs in Trichoderma guizhouense and Trichoderma harzianum revealed that only epl1 is active at all stages of development but that it plays a minor role in interactions with other fungi and bacteria. The deletion of this gene results in increased colonization of tomato roots by Trichoderma spp. Similarly, biochemical tests of EPL1 heterologously produced by Pichia pastoris support the claims described above. Based on the results obtained, we conclude that the function of CPs is probably linked to their surfactant properties and the ability to modify the hyphosphere of submerged mycelia and, thus, facilitate the nutritional versatility of fungi. The effector-like functions do not sufficiently describe the diversity and evolution of these proteins in fungi, as they are also maintained, duplicated, or laterally transferred in the genomes of nonherbivore fungi. IMPORTANCE Cerato-platanins (CPs) are surface-active small proteins abundantly secreted by filamentous fungi. Consequently, immune systems of plants and other organisms recognize CPs and activate defense mechanisms. Some CPs are toxic to plants and act as virulence factors in plant-pathogenic fungi. Our analysis, however, demonstrates that the interactions with plants do not explain the origin and evolution of CPs in the fungal kingdom. We revealed a long evolutionary history of CPs with multiple cases of gene duplication and events of interfungal lateral gene transfers. In the mycoparasitic Trichoderma spp., CPs evolve under stabilizing natural selection and hamper the colonization of roots. We propose that the ability to modify the hydrophobicity of the fungal hyphosphere is a key to unlock the evolutionary and functional paradox of these proteins.

Staphylococcus agnetis has been previously associated with subclinical or clinically mild cases of mastitis in dairy cattle and is one of several staphylococcal species that have been isolated from the bones and blood of lame broilers. We reported that S. agnetis could be obtained frequently from bacterial chondronecrosis with osteomyelitis (BCO) lesions of lame broilers (A. Al-Rubaye et al., PLoS One 10:e0143336, 2015 [https://doi.org/10.1371/journal.pone.0143336]). A particular isolate, S. agnetis 908, can induce lameness in over 50% of exposed chickens, exceeding normal BCO incidences in broiler operations. We reported the assembly and annotation of the genome of isolate 908. To better understand the relationship between dairy cattle and broiler isolates, we assembled 11 additional genomes for S. agnetis isolates, an additional chicken BCO strain, and ten isolates from cattle milk, mammary gland secretions, or udder skin from the collection at the University of Missouri. To trace phylogenetic relationships, we constructed phylogenetic trees based on multilocus sequence typing and genome-to-genome distance comparisons. Chicken isolate 908 clustered with two of the cattle isolates, along with three isolates from chickens in Denmark and an isolate of S. agnetis we isolated from a BCO lesion on a commercial broiler farm in Arkansas. We used a number of BLAST tools to compare the chicken isolates to those from cattle and identified 98 coding sequences distinguishing isolate 908 from the cattle isolates. None of the identified genes explain the differences in host or tissue tropism. These analyses are critical to understanding how staphylococci colonize and infect different hosts and potentially how they can transition to alternative niches (bone versus dermis). IMPORTANCE Staphylococcus agnetis has been recently recognized as associated with disease in dairy cattle and meat-type chickens. The infections appear to be limited in cattle and systemic in broilers. This report details the molecular relationships between cattle and chicken isolates in order to understand how this recently recognized species infects different hosts with different disease manifestations. The data show that the chicken and cattle isolates are very closely related, but the chicken isolates all cluster together, suggesting a single jump from cattle to chickens.

Pathogen introductions have led to numerous disease outbreaks in naive regions of the globe. The plant pathogen Xylella fastidiosa has been associated with various recent epidemics in Europe affecting agricultural crops, such as almond, grapevine, and olive, but also endemic species occurring in natural forest landscapes and ornamental plants. We compared whole-genome sequences of X. fastidiosa subspecies multiplex from America and strains associated with recent outbreaks in southern Europe to infer their likely origins and paths of introduction within and between the two continents. Phylogenetic analyses indicated multiple introductions of X. fastidiosa subspecies multiplex into Italy, Spain, and France, most of which emerged from a clade with limited genetic diversity with a likely origin in California, USA. The limited genetic diversity observed in X. fastidiosa subspecies multiplex strains originating from California is likely due to the clade itself being an introduction from X. fastidiosa subspecies multiplex populations in the southeastern United States, where this subspecies is most likely endemic. Despite the genetic diversity found in some areas in Europe, there was no clear evidence of recombination occurring among introduced X. fastidiosa strains in Europe. Sequence type taxonomy, based on multilocus sequence typing (MLST), was shown, at least in one case, to not lead to monophyletic clades of this pathogen; whole-genome sequence data were more informative in resolving the history of introductions than MLST data. Although additional data are necessary to carefully tease out the paths of these recent dispersal events, our results indicate that whole-genome sequence data should be considered when developing management strategies for X. fastidiosa outbreaks., This work was funded by the following awards: XF-ACTORS (Xylella fastidiosa Active Containment Through a Multidisciplinary-Oriented Research Strategy, grant 727987 from the European Union’s Horizon 2020 Framework Research Programme), E-RTA2017-00004-C02 (Desarrollo de estrategias de erradicación, contención y control de X. fastidiosa en España from “Programa Estatal de I+D+I Orientada a los Retos de la Sociedad of the Spanish Government” and FEDER), the Spanish Olive Oil Interprofessional, and the California Department of Food and Agriculture Pierce’ Disease Research Program.

Since antimicrobial usage is low in Norwegian animal husbandry, Norway is an ideal country to study antimicrobial resistance in the absence of selective pressure from antimicrobial usage. In particular, the usage of quinolones is very low, which makes it possible to investigate the spread and development of quinolone resistance in natural environments. Comparison of quinolone-resistant E. coli (QREC) isolates from livestock and wild animals in light of this low quinolone usage provides new insights into the development and dissemination of QREC in both natural and production environments. With this information, preventive measures may be taken to prevent further dissemination within Norwegian livestock and between other animals, thus maintaining the favorable situation in Norway., In Norway, the use of quinolones in livestock populations is very low, and prophylactic use is prohibited. Despite this, quinolone-resistant Escherichia coli (QREC) isolates are present at low levels in several animal species. The source of these QREC isolates is unknown. The aim of this study was to characterize and compare QREC isolates from different animal species to identify putative factors that may promote the occurrence of QREC. A total of 280 QREC isolates, from broilers, pigs, red foxes, and wild birds, were whole-genome sequenced and analyzed. Well-known chromosomal and plasmid-mediated resistance mechanisms were identified. In addition, mutations in marR, marA, and rpoB causing novel amino acid substitutions in their respective proteins were detected. Phylogenetic analyses were used to determine the relationships between the isolates. Quinolone resistance mechanism patterns appeared to follow sequence type groups. Similar QREC isolates with similar resistance mechanism patterns were detected from the samples, and further phylogenetic analysis indicated close evolutionary relationships between specific isolates from different sources. This suggests the dissemination of highly similar QREC isolates between animal species and also the persistence of QREC strains within the broiler production chain. This highlights the importance of both control measures at the top of the production chain as well as biosecurity measures to avoid the further dissemination and persistence of QREC in these environments. IMPORTANCE Since antimicrobial usage is low in Norwegian animal husbandry, Norway is an ideal country to study antimicrobial resistance in the absence of selective pressure from antimicrobial usage. In particular, the usage of quinolones is very low, which makes it possible to investigate the spread and development of quinolone resistance in natural environments. Comparison of quinolone-resistant E. coli (QREC) isolates from livestock and wild animals in light of this low quinolone usage provides new insights into the development and dissemination of QREC in both natural and production environments. With this information, preventive measures may be taken to prevent further dissemination within Norwegian livestock and between other animals, thus maintaining the favorable situation in Norway.

Aerobic anoxygenic phototrophic bacteria (AAPB) are thought to be important players in oceanic carbon and energy cycling in the euphotic zone of the ocean. The genus Citromicrobium , widely found in oligotrophic oceans, is a member of marine alphaproteobacterial AAPB. Nine Citromicrobium strains isolated from the South China Sea, the Mediterranean Sea, or the tropical South Atlantic Ocean were found to harbor identical 16S rRNA sequences. The sequencing of their genomes revealed high synteny in major regions. Nine genetic islands (GIs) involved mainly in type IV secretion systems, flagellar biosynthesis, prophage, and integrative conjugative elements, were identified by a fine-scale comparative genomics analysis. These GIs played significant roles in genomic evolution and divergence. Interestingly, the coexistence of two different photosynthetic gene clusters (PGCs) was not only found in the analyzed genomes but also confirmed, for the first time, to our knowledge, in environmental samples. The prevalence of the coexistence of two different PGCs may suggest an adaptation mechanism for Citromicrobium members to survive in the oceans. Comparison of genomic characteristics (e.g., GIs, average nucleotide identity [ANI], single-nucleotide polymorphisms [SNPs], and phylogeny) revealed that strains within a marine region shared a similar evolutionary history that was distinct from that of strains isolated from other regions (South China Sea versus Mediterranean Sea). Geographic differences are partly responsible for driving the observed genomic divergences and allow microbes to evolve through local adaptation. Three Citromicrobium strains isolated from the Mediterranean Sea diverged millions of years ago from other strains and evolved into a novel group. IMPORTANCE Aerobic anoxygenic phototrophic bacteria are a widespread functional group in the upper ocean, and their abundance could be up to 15% of the total heterotrophic bacteria. To date, a great number of studies display AAPB biogeographic distribution patterns in the ocean; however, little is understood about the geographic isolation impact on the genome divergence of marine AAPB. In this study, we compare nine Citromicrobium genomes of strains that have identical 16S rRNA sequences but different ocean origins. Our results reveal that strains isolated from the same marine region share a similar evolutionary history that is distinct from that of strains isolated from other regions. These Citromicrobium strains diverged millions of years ago. In addition, the coexistence of two different PGCs is prevalent in the analyzed genomes and in environmental samples.

Lactobacillus crispatus and Lactobacillus iners are common inhabitants of the healthy human vagina. These two species are closely related and are thought to perform similar ecological functions in the vaginal environment. Temporal data on the vaginal microbiome have shown that nontransient instances of cooccurrence are uncommon, while transitions from an L. iners -dominated community to one dominated by L. crispatus , and vice versa, occur often. This suggests that there is substantial overlap in the fundamental niches of these species. Given this apparent niche overlap, it is unclear how they have been maintained as common inhabitants of the human vagina. In this study, we characterized and compared the genomes of L. iners and L. crispatus to gain insight into possible mechanisms driving the maintenance of this species diversity. Our results highlight differences in the genomes of these two species that may facilitate the partitioning of their shared niche space. Many of the identified differences may impact the protective benefits provided to the host by these two species. IMPORTANCE The microbial communities that inhabit the human vagina play a critical role in the maintenance of vaginal health through the production of lactic acid and lowering the environmental pH. This precludes the growth of nonindigenous organisms and protects against infectious disease. The two most common types of vaginal communities are dominated by either Lactobacillus iners or Lactobacillus crispatus , while some communities alternate between the two over time. We combined ecological theory with state-of-the-art genome analyses to characterize how these two species might partition their shared niche space in the vagina. We show that the genomes of L. iners and L. crispatus differ in many respects, several of which may drive differences in their competitive abilities in the vagina. Our results provide insight into factors that drive the complicated temporal dynamics of the vaginal microbiome and demonstrate how closely related microbial species partition shared fundamental niche space.

Bacteria belonging to the Pseudovibrio genus are widespread, metabolically versatile, and able to thrive as both free-living and host-associated organisms. Although more than 50 genomes are available, a comprehensive comparative genomics study to resolve taxonomic inconsistencies is currently missing. We analyzed all available genomes and used 552 core genes to perform a robust phylogenomic reconstruction. This in-depth analysis revealed the divergence of two monophyletic basal lineages of strains isolated from polyclad flatworm hosts, namely, Pseudovibrio hongkongensis and Pseudovibrio stylochi. These strains have reduced genomes and lack sulfur-related metabolisms and major biosynthetic gene clusters, and their environmental distribution appears to be tightly associated with invertebrate hosts. We showed experimentally that the divergent strains are unable to utilize various sulfur compounds that, in contrast, can be utilized by the type strain Pseudovibrio denitrificans. Our analyses suggest that the lineage leading to these two strains has been subject to relaxed purifying selection resulting in great gene loss. Overall genome relatedness indices (OGRI) indicate substantial differences between the divergent strains and the rest of the genus. While 16S rRNA gene analyses do not support the establishment of a different genus for the divergent strains, their substantial genomic, phylogenomic, and physiological differences strongly suggest a divergent evolutionary trajectory and the need for their reclassification. Therefore, we propose the novel genus Polycladidibacter gen. nov. IMPORTANCE The genus Pseudovibrio is commonly associated with marine invertebrates, which are essential for ocean health and marine nutrient cycling. Traditionally, the phylogeny of the genus has been based on 16S rRNA gene analysis. The use of the 16S rRNA gene or any other single marker gene for robust phylogenetic placement has recently been questioned. We used a large set of marker genes from all available Pseudovibrio genomes for in-depth phylogenomic analyses. We identified divergent monophyletic basal lineages within the Pseudovibrio genus, including two strains isolated from polyclad flatworms. These strains showed reduced sulfur metabolism and biosynthesis capacities. The phylogenomic analyses revealed distinct evolutionary trajectories and ecological adaptations that differentiate the divergent strains from the other Pseudovibrio members and suggest that they fall into a novel genus. Our data show the importance of widening the use of phylogenomics for better understanding bacterial physiology, phylogeny, and evolution.

Feeding a global population of 8 billion people and climate change are the primary challenges facing agriculture today. Ruminant livestock are important food-producing animals, and maximizing their productivity requires an understanding of their digestive systems and the roles played by rumen microbes in plant polysaccharide degradation. Members of the genera Butyrivibrio and Pseudobutyrivibrio are a phylogenetically diverse group of bacteria and are commonly found in the rumen, where they are a substantial source of polysaccharide-degrading enzymes for the depolymerization of lignocellulosic material. Our findings have highlighted the immense enzymatic machinery of Butyrivibrio and Pseudobutyrivibrio species for the degradation of plant fiber, suggesting that these bacteria occupy similar niches but apply different degradation strategies in order to coexist in the competitive rumen environment., Plant polysaccharide breakdown by microbes in the rumen is fundamental to digestion in ruminant livestock. Bacterial species belonging to the rumen genera Butyrivibrio and Pseudobutyrivibrio are important degraders and utilizers of lignocellulosic plant material. These bacteria degrade polysaccharides and ferment the released monosaccharides to yield short-chain fatty acids that are used by the ruminant for growth and the production of meat, milk, and fiber products. Although rumen Butyrivibrio and Pseudobutyrivibrio species are regarded as common rumen inhabitants, their polysaccharide-degrading and carbohydrate-utilizing enzymes are not well understood. In this study, we analyzed the genomes of 40 Butyrivibrio and 6 Pseudobutyrivibrio strains isolated from the plant-adherent fraction of New Zealand dairy cows to explore the polysaccharide-degrading potential of these important rumen bacteria. Comparative genome analyses combined with phylogenetic analysis of their 16S rRNA genes and short-chain fatty acid production patterns provide insight into the genomic diversity and physiology of these bacteria and divide Butyrivibrio into 3 species clusters. Rumen Butyrivibrio bacteria were found to encode a large and diverse spectrum of degradative carbohydrate-active enzymes (CAZymes) and binding proteins. In total, 4,421 glycoside hydrolases (GHs), 1,283 carbohydrate esterases (CEs), 110 polysaccharide lyases (PLs), 3,605 glycosyltransferases (GTs), and 1,706 carbohydrate-binding protein modules (CBM) with predicted activities involved in the depolymerization and transport of the insoluble plant polysaccharides were identified. Butyrivibrio genomes had similar patterns of CAZyme families but varied greatly in the number of genes within each category in the Carbohydrate-Active Enzymes database (CAZy), suggesting some level of functional redundancy. These results suggest that rumen Butyrivibrio species occupy similar niches but apply different degradation strategies to be able to coexist in the rumen. IMPORTANCE Feeding a global population of 8 billion people and climate change are the primary challenges facing agriculture today. Ruminant livestock are important food-producing animals, and maximizing their productivity requires an understanding of their digestive systems and the roles played by rumen microbes in plant polysaccharide degradation. Members of the genera Butyrivibrio and Pseudobutyrivibrio are a phylogenetically diverse group of bacteria and are commonly found in the rumen, where they are a substantial source of polysaccharide-degrading enzymes for the depolymerization of lignocellulosic material. Our findings have highlighted the immense enzymatic machinery of Butyrivibrio and Pseudobutyrivibrio species for the degradation of plant fiber, suggesting that these bacteria occupy similar niches but apply different degradation strategies in order to coexist in the competitive rumen environment.

The chicken gastrointestinal tract harbors microorganisms that play a role in the health and disease status of the host. The cecum is the part of the gut that carries the highest microbial densities, has the longest residence time of digesta, and is a vital site for urea recycling and water regulation. Therefore, the cecum provides a rich environment for bacteria to horizontally transfer genes between one another via mobile genetic elements such as plasmids and bacteriophages. In this study, we used broiler chicken cecum as a model to investigate antibiotic resistance genes that can be transferred in vitro from cecal flora to Salmonella enterica serovar Heidelberg. We used whole-genome sequencing and resistome enrichment to decipher the interactions between S. Heidelberg, the gut microbiome, and acquired antibiotic resistance. After 48 h of incubation of ceca under microaerophilic conditions, we recovered one S. Heidelberg isolate with an acquired IncK2 plasmid (88 kb) carrying an extended-spectrum-β-lactamase gene (bla(CMY-2)). In vitro, this plasmid was transferable between Escherichia coli and S. Heidelberg strains but transfer was unsuccessful between S. Heidelberg strains. An in-depth genetic characterization of transferred plasmids suggests that they share significant homology with P1-like phages. This study contributes to our understanding of horizontal gene transfer between an important foodborne pathogen and the chicken gut microbiome. IMPORTANCE S. Heidelberg is a clinically important serovar, linked to foodborne illness and among the top 5 serovars isolated from poultry in the United States and Canada. Acquisition of new genetic material from the microbial flora in the gastrointestinal tract of food animals, including broilers, may contribute to increased fitness of pathogens like S. Heidelberg and may increase their level of antibiotic tolerance. Therefore, it is critical to gain a better understanding of the interactions that occur between important pathogens and the commensals present in the animal gut and other agroecosystems. In this report, we show that the native flora in broiler ceca were capable of transferring mobile genetic elements carrying the AmpC β-lactamase (bla(CMY-2)) gene to an important foodborne pathogen, S. Heidelberg. The potential role for bacteriophage transduction is also discussed.

Salmonella enterica serovar Dublin is a cattle-adapted S. enterica serovar causing both intestinal and systemic infection in its bovine host, and it is also a serious threat to human health. The present study aimed to determine the population structure of S. Dublin isolates obtained from Danish cattle herds and to investigate how cattle isolates relate to Danish human isolates, as well as to non-Danish human and bovine isolates. Phylogenetic analysis of 197 Danish cattle isolates from 1996 to 2016 identified three major clades corresponding to distinct geographical regions of cattle herds. Persistence of closely related isolates within the same herd and their circulation between epidemiologically linked herds for a period of more than 20 years were demonstrated. These findings suggest that a lack of internal biosecurity and, to some extent, also a lack of external biosecurity in the herds have played an important role in the long-term persistence of S. Dublin in Danish cattle herds in the period investigated. Global population analysis revealed that Danish cattle isolates clustered separately from bovine isolates from other countries, whereas human isolates were geographically spread. Resistance genes were not commonly demonstrated in Danish bovine isolates; only the isolates within one Danish clade were found to often harbor two plasmids of IncFII/IncFIB and IncN types, the latter plasmid carrying bla(TEM-1), tetA, strA, and strB antibiotic resistance genes. IMPORTANCE S. Dublin causes economic losses in cattle production, and the bacterium is a public health concern. A surveillance and control program has been in place in Denmark since 2002 with the ultimate goal to eradicate S. Dublin from Danish cattle herds; however, a small proportion of herds have remained positive for many years. In this study, we demonstrate that herds with persistent infection often were infected with the same strain for many years, indicating that internal biosecurity has to be improved to curb the infection. Further, domestic cases of S. Dublin infection in humans were found to be caused both by Danish cattle isolates and by isolates acquired abroad. This study shows the strength of whole-genome sequencing to obtain detailed information on epidemiology of S. Dublin and allows us to suggest internal biosecurity as a main way to control this bacterium in Danish cattle herds.

While microbiome studies have focused on diversity on the species or higher level, bacterial species in microbiomes are represented by different, often multiple strains. These strains could be clonally and phenotypically very different, making assessment of strain content vital to a full understanding of microbiome function. This is especially important with respect to antibiotic resistant strains, the clonal spread of which may be dependent on competition between them and susceptible strains from the same species. The pandemic, multi-drug resistant, and highly pathogenic E. coli subclone ST131-H30 (H30) is of special interest, as it has already been found persisting in the gut and bladder of healthy people. In order to rapidly assess E. coli clonal diversity, we developed a novel method based on deep sequencing of two loci used for sequence typing, along with an algorithm for analysis of resulting data. Using this method, we assessed fecal and urinary samples from healthy women carrying H30, and were able to uncover considerable diversity, including strains with frequencies at E. coli population. We also found that even in the absence of antibiotic use, H30 could complete dominate the gut and, especially, urine of healthy carriers. Our study offers a novel tool for assessing a species’ clonal diversity (clonobiome) within the microbiome, that could be useful in studying population structure and dynamics of multi-drug resistant and/or highly pathogenic strains in their natural environments.IMPORTANCEBacterial species in the microbiome are often represented by multiple genetically and phenotypically different strains, making insight into subspecies diversity critical to a full understanding of the microbiome, especially with respect to opportunistic pathogens. However, methods allowing efficient high-throughput clonal typing are not currently available. This study combines a conventional E. coli typing method with deep amplicon sequencing to allow analysis of many samples concurrently. While our method was developed for E. coli, it may be adapted for other species, allowing for microbiome researchers to assess clonal strain diversity in natural samples. Since assessment of subspecies diversity is particularly important for understanding the spread of antibiotic resistance, we applied our method to study of a pandemic multidrug-resistant E. coli clone. The results we present suggest that this clone could be highly competitive in healthy carriers, and that the mechanisms of colonization by such clones need to be studied.

It has been 30 years since the initial emergence and subsequent rapid global spread of multidrug-resistant Salmonella enterica serovar Typhimurium DT104 (MDR DT104). Nonetheless, its origin and transmission route have never been revealed. We used whole-genome sequencing (WGS) and temporally structured sequence analysis within a Bayesian framework to reconstruct temporal and spatial phylogenetic trees and estimate the rates of mutation and divergence times of 315 S . Typhimurium DT104 isolates sampled from 1969 to 2012 from 21 countries on six continents. DT104 was estimated to have emerged initially as antimicrobial susceptible in ∼1948 (95% credible interval [CI], 1934 to 1962) and later became MDR DT104 in ∼1972 (95% CI, 1972 to 1988) through horizontal transfer of the 13-kb Salmonella genomic island 1 (SGI1) MDR region into susceptible strains already containing SGI1. This was followed by multiple transmission events, initially from central Europe and later between several European countries. An independent transmission to the United States and another to Japan occurred, and from there MDR DT104 was probably transmitted to Taiwan and Canada. An independent acquisition of resistance genes took place in Thailand in ∼1975 (95% CI, 1975 to 1990). In Denmark, WGS analysis provided evidence for transmission of the organism between herds of animals. Interestingly, the demographic history of Danish MDR DT104 provided evidence for the success of the program to eradicate Salmonella from pig herds in Denmark from 1996 to 2000. The results from this study refute several hypotheses on the evolution of DT104 and suggest that WGS may be useful in monitoring emerging clones and devising strategies for prevention of Salmonella infections.

Staphylococcus aureus is a Gram-positive pathogenic bacterium that colonizes an estimated one-third of the human population and can cause a wide spectrum of disease, ranging from superficial skin infections to life-threatening sepsis. The adaptive mechanisms that contribute to the success of this pathogen remain obscure partially due to a lack of knowledge of its metabolic requirements. Systems biology approaches can be extremely useful in predicting and interpreting metabolic phenotypes; however, such approaches rely on a chemically defined minimal medium as a basis to investigate the requirements of the cell. In this study, a chemically defined minimal medium formulation, termed synthetic minimal medium (SMM), was investigated and validated to support growth of three S. aureus strains: LAC and TCH1516 (USA300 lineage), as well as D592 (USA100 lineage). The formulated SMM was used in an adaptive laboratory evolution experiment to probe the various mutational trajectories of all three strains leading to optimized growth capabilities. The evolved strains were phenotypically characterized for their growth rate and antimicrobial susceptibility. Strains were also resequenced to examine the genetic basis for observed changes in phenotype and to design follow-up metabolite supplementation assays. Our results reveal evolutionary trajectories that arose from strain-specific metabolic requirements. SMM and the evolved strains can also serve as important tools to study antibiotic resistance phenotypes of S. aureus. IMPORTANCE As researchers try to understand and combat the development of antibiotic resistance in pathogens, there is a growing need to thoroughly understand the physiology and metabolism of the microbes. Staphylococcus aureus is a threatening pathogen with increased antibiotic resistance and well-studied virulence mechanisms. However, the adaptive mechanisms used by this pathogen to survive environmental stresses remain unclear, mostly due to the lack of information about its metabolic requirements. Defining the minimal metabolic requirements for S. aureus growth is a first step toward unraveling the mechanisms by which it adapts to metabolic stresses. Here, we present the development of a chemically defined minimal medium supporting growth of three S. aureus strains, and we reveal key genetic mutations contributing to improved growth in minimal medium.

Survival and growth of the anaerobic gut fungi (AGF; Neocallimastigomycota) in the herbivorous gut necessitate the possession of multiple abilities absent in other fungal lineages. We hypothesized that horizontal gene transfer (HGT) was instrumental in forging the evolution of AGF into a phylogenetically distinct gut-dwelling fungal lineage. The patterns of HGT were evaluated in the transcriptomes of 27 AGF strains, 22 of which were isolated and sequenced in this study, and 4 AGF genomes broadly covering the breadth of AGF diversity. We identified 277 distinct incidents of HGT in AGF transcriptomes, with subsequent gene duplication resulting in an HGT frequency of 2 to 3.5% in AGF genomes. The majority of HGT events were AGF specific (91.7%) and wide (70.8%), indicating their occurrence at early stages of AGF evolution. The acquired genes allowed AGF to expand their substrate utilization range, provided new venues for electron disposal, augmented their biosynthetic capabilities, and facilitated their adaptation to anaerobiosis. The majority of donors were anaerobic fermentative bacteria prevalent in the herbivorous gut. This study strongly indicates that HGT indispensably forged the evolution of AGF as a distinct fungal phylum and provides a unique example of the role of HGT in shaping the evolution of a high-rank taxonomic eukaryotic lineage. IMPORTANCE The anaerobic gut fungi (AGF) represent a distinct basal phylum lineage (Neocallimastigomycota) commonly encountered in the rumen and alimentary tracts of herbivores. Survival and growth of anaerobic gut fungi in these anaerobic, eutrophic, and prokaryote-dominated habitats necessitates the acquisition of several traits absent in other fungal lineages. We assess here the role of horizontal gene transfer as a relatively fast mechanism for trait acquisition by the Neocallimastigomycota postsequestration in the herbivorous gut. Analysis of 27 transcriptomes that represent the broad diversity of Neocallimastigomycota identified 277 distinct HGT events, with subsequent gene duplication resulting in an HGT frequency of 2 to 3.5% in AGF genomes. These HGT events have allowed AGF to survive in the herbivorous gut by expanding their substrate utilization range, augmenting their biosynthetic pathway, providing new routes for electron disposal by expanding fermentative capacities, and facilitating their adaptation to anaerobiosis. HGT in the AGF is also shown to be mainly a cross-kingdom affair, with the majority of donors belonging to the bacteria. This study represents a unique example of the role of HGT in shaping the evolution of a high-rank taxonomic eukaryotic lineage.