Streptomyces bacteria are known for their prolific production of secondary metabolites, many of which have been widely used in human medicine, agriculture and animal health. To guide the effective prioritization of specific biosynthetic gene clusters (BGCs) for drug development and targeting the most prolific producer strains, knowledge about phylogenetic relationships of Streptomyces species, genomewide diversity and distribution patterns of BGcs is critical. We used genomic and phylogenetic methods to elucidate the diversity of major classes of BGCs in 1,110 publicly available Streptomyces genomes. Genome mining of Streptomyces reveals high diversity of BGcs and variable distribution patterns in the Streptomyces phylogeny, even among very closely related strains. The most common BGCs are non-ribosomal peptide synthetases, type 1 polyketide synthases, terpenes, and lantipeptides. We also found that numerous Streptomyces species harbor BGcs known to encode antitumor compounds. We observed that strains that are considered the same species can vary tremendously in the BGcs they carry, suggesting that strain-level genome sequencing can uncover high levels of BGC diversity and potentially useful derivatives of any one compound. These findings suggest that a strain-level strategy for exploring secondary metabolites for clinical use provides an alternative or complementary approach to discovering novel pharmaceutical compounds from microbes.

Streptomyces rimosus is best known as the primary source of the tetracycline class of antibiotics, most notably oxytetracycline, which have been widely used against many gram-positive and gram-negative pathogens and protozoan parasites. However, despite the medical and agricultural importance of S. rimosus , little is known of its evolutionary history and genome dynamics. In this study, we aim to elucidate the pan-genome characteristics and phylogenetic relationships of 32 S. rimosus genomes. The S. rimosus pan-genome contains more than 22,000 orthologous gene clusters, and approximately 8.8% of these genes constitutes the core genome. A large part of the accessory genome is composed of 9,646 strain-specific genes. S. rimosus exhibits an open pan-genome (decay parameter α = 0.83) and high gene diversity between strains (genomic fluidity φ = 0.12). We also observed strain-level variation in the distribution and abundance of biosynthetic gene clusters (BGCs) and that each individual S. rimosus genome has a unique repertoire of BGCs. Lastly, we observed variation in recombination, with some strains donating or receiving DNA more often than others, strains that tend to frequently recombine with specific partners, genes that often experience recombination more than others, and variable sizes of recombined DNA sequences. We conclude that the high levels of inter-strain genomic variation in S. rimosus is partly explained by differences in recombination among strains. These results have important implications on current efforts for natural drug discovery, the ecological role of strain-level variation in microbial populations, and addressing the fundamental question of why microbes have pan-genomes.

Salmonella is responsible for many nontyphoidal foodborne infections and enteric (typhoid) fever in humans. Of the two Salmonella species, Salmonella enterica is highly diverse and includes 10 known subspecies and approximately 2,600 serotypes. Understanding the evolutionary processes that generate the tremendous diversity in Salmonella is important in reducing and controlling the incidence of disease outbreaks and the emergence of virulent strains. In this study, we aim to elucidate the impact of homologous recombination in the diversification of S. enterica subspecies. Using a data set of previously published 926 Salmonella genomes representing the 10 S. enterica subspecies and Salmonella bongori, we calculated a genus-wide pan-genome composed of 84,041 genes and the S. enterica pan-genome of 81,371 genes. The size of the accessory genomes varies between 12,429 genes in S. enterica subsp. arizonae (subsp. IIIa) to 33,257 genes in S. enterica subsp. enterica (subsp. I). A total of 12,136 genes in the Salmonella pan-genome show evidence of recombination, representing 14.44% of the pan-genome. We identified genomic hot spots of recombination that include genes associated with flagellin and the synthesis of methionine and thiamine pyrophosphate, which are known to influence host adaptation and virulence. Last, we uncovered within-species heterogeneity in rates of recombination and preferential genetic exchange between certain donor and recipient strains. Frequent but biased recombination within a bacterial species may suggest that lineages vary in their response to environmental selection pressure. Certain lineages, such as the more uncommon non-enterica subspecies (non-S. enterica subsp. enterica), may also act as a major reservoir of genetic diversity for the wider population. IMPORTANCE S. enterica is a major foodborne pathogen, which can be transmitted via several distinct routes from animals and environmental sources to human hosts. Multiple subspecies and serotypes of S. enterica exhibit considerable differences in virulence, host specificity, and colonization. This study provides detailed insights into the dynamics of recombination and its contributions to S. enterica subspecies evolution. Widespread recombination within the species means that new adaptations arising in one lineage can be rapidly transferred to another lineage. We therefore predict that recombination has been an important factor in the emergence of several major disease-causing strains from diverse genomic backgrounds and their ability to adapt to disparate environments.

Campylobacter jejuni is one of the leading causes of bacterial gastroenteritis worldwide. In the United States, New Hampshire was one of the 18 states that reported cases in the 2016 to 2018 multistate outbreak of multidrug-resistant C. jejuni. Here, we aimed to elucidate the baseline diversity of the wider New Hampshire C. jejuni population during the outbreak. We used genome sequences of 52 clinical isolates sampled in New Hampshire in 2017, including 1 of the 2 isolates from the outbreak. Results revealed a remarkably diverse population composed of at least 28 sequence types, which are mostly represented by 1 or a few strains. A comparison of our isolates with 249 clinical C. jejuni from other states showed frequent phylogenetic intermingling, suggesting a lack of geographical structure and minimal local diversification within the state. Multiple independent acquisitions of resistance genes from 5 classes of antibiotics characterize the population, with 47/52 (90.4%) of the genomes carrying at least 1 horizontally acquired resistance gene. Frequently recombining genes include those associated with heptose biosynthesis, colonization, and stress resistance. We conclude that the diversity of clinical C. jejuni in New Hampshire in 2017 was driven mainly by the coexistence of phylogenetically diverse antibiotic-resistant lineages, widespread geographical mixing, and frequent recombination. This study provides an important baseline census of the standing pangenomic variation and drug resistance to aid the development of a statewide database for epidemiological studies and clinical decision making. Continued genomic surveillance will be necessary to accurately assess how the population of C. jejuni changes over the long term.

One free-ranging Gray fox (Urocyon cinereoargenteus) underwent autopsy following neurologic disease, with findings including morbilliviral inclusions and associated lesions in numerous tissues, adenoviral intranuclear inclusions in bronchial epithelial cells, and septic pleuropneumonia, hepatitis, splenitis, and meningoencephalitis. Molecular diagnostics on fresh lung identified a strain within a distinct clade of canine distemper that is currently unique to wildlife in New England, as well as the emerging multi-host viral pathogen skunk adenovirus-1. Bacterial culture of fresh liver resulted in a pure growth of Listeria monocytogenes, with whole genome sequencing indicating that the isolate had a vast array of antimicrobial resistance and virulence-associated genes. One year later, a second fox was euthanized for inappropriate behavior in a residential area, and diagnostic workup revealed canine distemper and septic L. monocytogenes, with the former closely related to the distemper virus found in the previous fox and the latter divergent from the L. monocytogenes from the previous fox.