Widespread interspecies homologous recombination reveals reticulate evolution within the genus Streptomyces

Current Methods for Recombination Detection in Bacteria

The role of genetic exchanges, i.e., homologous recombination (HR) and horizontal gene transfer (HGT), in bacteria cannot be overestimated for it is a pivotal mechanism leading to their evolution and adaptation, thus, tracking the signs of recombination and HGT events is importance both for fundamental and applied science. To date, dozens of bioinformatics tools for revealing recombination signals are available, however, their pros and cons as well as the spectra of solvable tasks have not yet been systematically reviewed. Moreover, there are two major groups of software. One aims to infer evidence of HR, while the other only deals with horizontal gene transfer (HGT). However, despite seemingly different goals, all the methods use similar algorithmic approaches, and the processes are interconnected in terms of genomic evolution influencing each other. In this review, we propose a classification of novel instruments for both HR and HGT detection based on the genomic consequences of recombination. In this context, we summarize available methodologies paying particular attention to the type of traceable events for which a certain program has been designed.

Streptomyces apricus sp. nov., isolated from soil

A novel Streptomyces strain, SUN51T, was isolated from soils sampled in Wisconsin, USA, as part of a Streptomyces biogeography survey. Genome sequencing revealed that this strain had less than 90 % average nucleotide identity (ANI) to type species of Streptomyces : SUN51T was most closely related to Streptomyces dioscori A217T (99.5 % 16S rRNA gene identity, 89.4 % ANI). Genome size was estimated at 8.81 Mb, and the genome DNA G+C content was 72 mol%. The strain possessed the cellular fatty acids anteiso-C15 : 0, iso-C16 : 0, 16 : 1 ω7c, anteiso-C17 : 0, iso-C14 : 0 and C16 : 0. The predominant menaquinones were MK-9 H4, MK-9 H6 and MK-9 H8. Strain SUN51T contained the polar lipids phosphatidic acid, phosphatidyl ethanolamine, phosphatidyl glycerol and diphosphatidyl glycerol. The cell wall contained ll-diaminopimelic acid. The strain could grow on a broad range of carbon sources and tolerate temperatures of up to 40 °C. The results of the polyphasic study confirmed that this isolate represents a novel species of the genus Streptomyces , for which the name Streptomyces apricus sp. nov. is proposed. The type strain of this species is SUN51T (=NRRL B-65543T=JCM 33736T).

Distribution and Evolutionary History of Sialic Acid Catabolism in the Phylum Actinobacteria

Sialic acids are present in humans and other metazoans, playing essential roles in physiological and pathological processes. Commensal and pathogenic bacteria have evolved the capacity to utilize sialic acids as nutrient and energy sources. However, in some actinobacteria, sialic acid catabolism (SAC) is associated with free-living populations. To unravel the distribution and evolutionary history of SAC in the phylum Actinobacteria, we analyzed the presence and diversity of the putative SAC gene cluster (nan) in 7,180 high-quality, nonredundant actinobacterial genomes that covered 1,969 species. The results showed that ∼13% of actinobacterial species had the potential to utilize sialic acids, with 45 species capable of anhydro-SAC, all except two of them through the canonical pathway. These species belonged to 20 orders and 81 genera, with ∼36% of them from four genera, Actinomyces, Bifidobacterium, Corynebacterium, and Streptomyces. Moreover, ∼40% of the nan-positive species are free living. Phylogenetic analysis of the key nan genes, nanA, nanK, and nanE, revealed a strong signal of horizontal gene transfer (HGT), accompanied with vertical inheritance and gene loss. This evolutionary pattern led to high diversity and differential distribution of nan among actinobacterial taxa and might cause the cluster to spread to some free-living species while losing in some host-associated species. The evolution of SAC in actinobacteria probably represents the evolution of certain kinds of noncore bacterial functions for environmental adaptation and lifestyle switch, in which HGT plays a dominant role. IMPORTANCE Sialic acids play essential roles in the physiology of humans and other metazoan animals, and microbial sialic acid catabolism (SAC) is one of the processes critical for pathogenesis. To date, microbial SAC is studied mainly in commensals and pathogens, while its distribution in free-living microbes and evolutionary pathway remain largely unexplored. Here, by examining all actinobacterial genomes available, we demonstrate that putative SAC is present in a small proportion of actinobacterial species, of which, however, ∼40% are free-living species. We also reveal remarkable difference in the distribution of SAC among actinobacterial taxa and high diversity of the putative SAC gene clusters. HGT plays a significant role in the evolution of SAC, accompanied with vertical inheritance and gene loss. Our results provide a comprehensive and systematic picture of the distribution and evolutionary history of SAC in actinobacteria, expanding the current knowledge on bacterial adaptation and diversification.

Comparative Genomic Insights Into the Taxonomic Classification, Diversity, and Secondary Metabolic Potentials of Kitasatospora, a Genus Closely Related to Streptomyces

While the genus Streptomyces (family Streptomycetaceae) has been studied as a model for bacterial secondary metabolism and genetics, its close relatives have been less studied. The genus Kitasatospora is the second largest genus in the family Streptomycetaceae. However, its taxonomic position within the family remains under debate and the secondary metabolic potential remains largely unclear. Here, we performed systematic comparative genomic and phylogenomic analyses of Kitasatospora. Firstly, the three genera within the family Streptomycetaceae (Kitasatospora, Streptomyces, and Streptacidiphilus) showed common genomic features, including high G + C contents, high secondary metabolic potentials, and high recombination frequencies. Secondly, phylogenomic and comparative genomic analyses revealed phylogenetic distinctions and genome content differences among these three genera, supporting Kitasatospora as a separate genus within the family. Lastly, the pan-genome analysis revealed extensive genetic diversity within the genus Kitasatospora, while functional annotation and genome content comparison suggested genomic differentiation among lineages. This study provided new insights into genomic characteristics of the genus Kitasatospora, and also uncovered its previously underestimated and complex secondary metabolism.

Proposal of Carbonactinosporaceae fam. nov. within the class Actinomycetia. Reclassification of Streptomyces thermoautotrophicus as Carbonactinospora thermoautotrophica gen. nov., comb. nov

Streptomyces thermoautotrophicus UBT1^T has been suggested to merit generic status due to its phylogenetic placement and distinctive phenotypes among Actinomycetia. To evaluate whether 'S. thermoautotrophicus' represents a higher taxonomic rank, 'S. thermoautotrophicus' strains UBT1^T and H1 were compared to Actinomycetia using 16S rRNA gene sequences and comparative genome analyses. The UBT1^T and H1 genomes each contain at least two different 16S rRNA sequences, which are closely related to those of Acidothermus cellulolyticus (order Acidothermales). In multigene-based phylogenomic trees, UBT1^T and H1 typically formed a sister group to the Streptosporangiales-Acidothermales clade. The Average Amino Acid Identity, Percentage of Conserved Proteins, and whole-genome Average Nucleotide Identity (Alignment Fraction) values were ≤58.5%, ≤48%, ≤75.5% (0.3) between 'S. thermoautotrophicus' and Streptosporangiales members, all below the respective thresholds for delineating genera. The values for genomics comparisons between strains UBT1^T and H1 with Acidothermales, as well as members of the genus Streptomyces, were even lower. A review of the 'S. thermoautotrophicus' proteomic profiles and KEGG orthology demonstrated that UBT1^T and H1 present pronounced differences, both tested and predicted, in phenotypic and chemotaxonomic characteristics compared to its sister clades and Streptomyces. The distinct phylogenetic position and the combination of genotypic and phenotypic characteristics justify the proposal of Carbonactinospora gen. nov., with the type species Carbonactinospora thermoautotrophica comb. nov. (type strain UBT1^T, = DSM 100163^T = KCTC 49540^T) belonging to Carbonactinosporaceae fam. nov. within Actinomycetia.

The biogeography of Streptomyces in New Zealand enabled by high-throughput sequencing of genus-specific rpoB amplicons

We evaluated Streptomyces biogeography in soils along a 1200 km latitudinal transect across New Zealand (NZ). Streptomyces diversity was examined using high-throughput sequencing of rpoB amplicons generated with a Streptomyces specific primer set. We detected 1287 Streptomyces rpoB operational taxonomic units (OTUs) with 159 ± 92 (average ± SD) rpoB OTUs per site. Only 12% (n = 149) of these OTUs matched rpoB sequences from cultured specimens (99% nucleotide identity cutoff). Streptomyces phylogenetic diversity (Faith's PD) was correlated with soil pH, mean annual temperature and plant community richness (Spearman's r: 0.77, 0.64 and -0.79, respectively; P < 0.05), but not with latitude. In addition, soil pH and plant community richness both explained significant variation in Streptomyces beta diversity. Streptomyces communities exhibited both high dissimilarity and strong dominance of one or a few species at each site. Taken together, these results suggest that dispersal limitation due to competitive interactions limits the colonization success of spores that relocate to new sites. Cultivated Streptomyces isolates represent a major source of clinically useful antibiotics, but only a small fraction of extant diversity within the genus have been identified and most species of Streptomyces have yet to be described.

Telomeric and sub-telomeric regions undergo rapid turnover within a Streptomyces population

Genome dynamics was investigated within natural populations of the soil bacterium Streptomyces. The exploration of a set of closely related strains isolated from micro-habitats of a forest soil exhibited a strong diversity of the terminal structures of the linear chromosome, i.e. terminal inverted repeats (TIRs). Large insertions, deletions and translocations could be observed along with evidence of transfer events between strains. In addition, the telomere and its cognate terminal protein complexes required for terminal replication and chromosome maintenance, were shown to be variable within the population probably reflecting telomere exchanges between the chromosome and other linear replicons (i.e., plasmids). Considering the close genetic relatedness of the strains, these data suggest that the terminal regions are prone to a high turnover due to a high recombination associated with extensive horizontal gene transfer.

Evolution and Diversity of Assembly-Line Polyketide Synthases

Assembly-line polyketide synthases (PKSs) are among the most complex protein machineries known in nature, responsible for the biosynthesis of numerous compounds used in the clinic. Their present-day diversity is the result of an evolutionary path that has involved the emergence of a multimodular architecture and further diversification of assembly-line PKSs. In this review, we provide an overview of previous studies that investigated PKS evolution and propose a model that challenges the currently prevailing view that gene duplication has played a major role in the emergence of multimodularity. We also analyze the ensemble of orphan PKS clusters sequenced so far to evaluate how large the entire diversity of assembly-line PKS clusters and their chemical products could be. Finally, we examine the existing techniques to access the natural PKS diversity in natural and heterologous hosts and describe approaches to further expand this diversity through engineering.

Population Genomics Insights into Adaptive Evolution and Ecological Differentiation in Streptomycetes

Deciphering the genomic variation that represents microevolutionary processes toward species divergence is key to understanding microbial speciation, which has long been under debate. Streptomycetes are filamentous bacteria that are ubiquitous in nature and the richest source of antibiotics; however, their speciation processes remain unknown. To tackle this issue, we performed a comprehensive population genomics analysis on Streptomyces albidoflavus residing in different habitats and with a worldwide distribution and identified and characterized the foundational changes within the species. We detected three well-defined phylogenomic clades, of which clades I and III mainly contained free-living (soil/marine) and insect-associated strains, respectively, and clade II had a mixed origin. By performing genome-wide association studies (GWAS), we identified a number of genetic variants associated with free-living or entomic (denoting or relating to insects) habitats in both the accessory and core genomes. These variants contributed collectively to the population structure and had annotated or confirmed functions that likely facilitate differential adaptation of the species. In addition, we detected higher levels of homologous recombination within each clade and in the free-living group than within the whole species and in the entomic group. A subset of the insect-associated strains (clade III) showed a relatively independent evolutionary trajectory with more symbiosis-favorable genes but little genetic interchange with the other lineages. Our results demonstrate that ecological adaptation promotes genetic differentiation in S. albidoflavus, suggesting a model of ecological speciation with gene flow in streptomycetes.IMPORTANCE Species are the fundamental units of ecology and evolution, and speciation leads to the astounding diversity of life on Earth. Studying speciation is thus of great significance to understand, protect, and exploit biodiversity, but it is a challenge in the microbial world. In this study, using population genomics, we placed Streptomyces albidoflavus strains in a spectrum of speciation and showed that the genetic differences between phylogenomic clusters evolved mainly by environmental selection and gene-specific sweeps. These findings highlight the role of ecology in structuring recombining bacterial species, making a step toward a deeper understanding of microbial speciation. Our results also raise concerns of an underrated microbial diversity at the intraspecies level, which can be utilized for mining of ecologically relevant natural products.

Within-Species Genomic Variation and Variable Patterns of Recombination in the Tetracycline Producer Streptomyces rimosus

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.

Impact of Homologous Recombination on the Evolution of Prokaryotic Core Genomes

Microbial populations exchange genetic material through a process called homologous recombination. Although this process has been studied in particular organisms, we lack an understanding of its differential impact over the genome and across microbes with different life-styles. We used a common analytical framework to assess this process in a representative set of microorganisms. Our results uncovered important trends. First, microbes with different lifestyles are differentially impacted, with endosymbionts and obligate pathogens being those less prone to undergo this process. Second, certain genetic elements such as restriction-modification systems seem to be associated with higher rates of recombination. Most importantly, recombined genomes show the footprints of natural selection in which recombined regions preferentially contain genes that can be related to specific ecological adaptations. Taken together, our results clarify the relative contributions of factors modulating homologous recombination and show evidence for a clear a role of this process in shaping microbial genomes and driving ecological adaptations.

Homologous recombination (HR) enables the exchange of genetic material between and within species. Recent studies suggest that this process plays a major role in the microevolution of microbial genomes, contributing to core genome homogenization and to the maintenance of cohesive population structures. However, we still have a very poor understanding of the possible adaptive roles of intraspecific HR and of the factors that determine its differential impact across clades and lifestyles. Here we used a unified methodological framework to assess HR in 338 complete genomes from 54 phylogenetically diverse and representative prokaryotic species, encompassing different lifestyles and a broad phylogenetic distribution. Our results indicate that lifestyle and presence of restriction-modification (RM) machineries are among the main factors shaping HR patterns, with symbionts and intracellular pathogens having the lowest HR levels. Similarly, the size of exchanged genomic fragments correlated with the presence of RM and competence machineries. Finally, genes exchanged by HR showed functional enrichments which could be related to adaptations to different environments and ecological strategies. Taken together, our results clarify the factors underlying HR impact and suggest important adaptive roles of genes exchanged through this mechanism. Our results also revealed that the extent of genetic exchange correlated with lifestyle and some genomic features. Moreover, the genes in exchanged regions were enriched for functions that reflected specific adaptations, supporting identification of HR as one of the main evolutionary mechanisms shaping prokaryotic core genomes.

Streptomyces carminius sp. nov., a novel actinomycete isolated from Sophora alopecuroides in Xinjiang, China

A novel actinobacterial strain, designated TRM SA0054^T, was isolated from the roots of Sophora alopecuroides grown in Alar, Xinjiang, north-west China, and characterised by a polyphasic taxonomic approach. Phylogenetic analysis showed that strain TRM SA0054^T has 16S rRNA gene sequence similarity of 98.22% with Streptomyces barkulensis RC 1831^T. Whole cell hydrolysates of strain TRM SA0054^T were found to contain LL-diaminopimelic acid as the diagnostic diamino acid and ribose, glucose and xylose as the major whole cell sugars. The major fatty acids were identified as iso-C_16:0, iso-C_16:1 G, anteiso-C_17:0, anteiso-C_15:0, C_16:0, anteiso-C_17:1 ω9c. The main menaquinones were determined to be MK-8 (H₄), MK-9 (H₆) and MK-9 (H₈). The polar lipids were identified as diphosphatidylglycerol, phosphatidylmethylethanolamine, phosphatidylethanolamine, phosphatidylcholine, phosphatidylinositol, phosphatidylinositol mannoside and an unidentified lipid. The G + C content of the genomic DNA was determined to be 73.04%. Strain TRM SA0054^T has a relatively low DNA-DNA relatedness value with Streptomyces barkulensis RC 1831^T as determined by calculating the average nucleotide identity value (ANI = 84.1%). Based on the phenotypic, chemotaxonomic and phylogenetic data, it is concluded that strain TRM SA0054^T should be designated as a novel species of the genus Streptomyces, for which the name Streptomyces carminius sp. nov. is proposed, with type strain TRM SA0054^T (= CCTCC AA 2016041^T = KCTC39903^T).

Dibutyl phthalate alters the metabolic pathways of microbes in black soils

Dibutyl phthalate (DBP) is well known as a high-priority pollutant. This study explored the impacts of DBP on the metabolic pathways of microbes in black soils in the short term (20 days). The results showed that the microbial communities were changed in black soils with DBP. In nitrogen cycling, the abundances of the genes were elevated by DBP. DBP contamination facilitated 3'-phosphoadenosine-5'-phosphosulfate (PAPS) formation, and the gene flux of sulfate metabolism was increased. The total abundances of ABC transporters and the gene abundances of the monosaccharide-transporting ATPases MalK and MsmK were increased by DBP. The total abundance of two-component system (TCS) genes and the gene abundances of malate dehydrogenase, histidine kinase and citryl-CoA lyase were increased after DBP contamination. The total abundance of phosphotransferase system (PTS) genes and the gene abundances of phosphotransferase, Crr and BglF were raised by DBP. The increased gene abundances of ABC transporters, TCS and PTS could be the reasons for the acceleration of nitrogen, carbon and sulfate metabolism. The degrading-genes of DBP were increased markedly in soil exposed to DBP. In summary, DBP contamination altered the microbial community and enhanced the gene abundances of the carbon, nitrogen and sulfur metabolism in black soils in the short term.

Multilocus sequence analysis of homologous recombination and diversity in Arthrobacter sensu lato named species and glacier-inhabiting strains

Members of the bacterial genus Arthrobacter sensu lato are Gram-positive actinomycetes distributed worldwide and found in numerous environments including soil, water, glacier ice, and sewage. Homologous recombination is an important driving force in bacterial evolution, but its impact on Arthrobacter sensu lato evolution is poorly understood. We evaluated homologous recombination among 41 Arthrobacter sensu lato named species, using multilocus sequence analysis (MLSA). A high level of recombination was found, associated with strong diversification and a reticulate evolutionary pattern of Arthrobacter sensu lato. We also collected a total of 31 cold-adapted Arthrobacter sensu lato strains from two cold glaciers located in northwest China and two temperate glaciers in southwest China, and evaluated their diversity and population structure by MLSA. The glacier strains displayed high diversity, but rates of recombination among the four glacier groups were quite low, indicating that barriers to homologous recombination formed in the past among the populations on different glaciers. Our findings indicate that historical glaciation events shaped the contemporary distributions, taxonomic relationships, and phylogeographic patterns of Arthrobacter sensu lato species on glaciers.

Phylogenomic Insight into Salinispora (Bacteria, Actinobacteria) Species Designations

Bacteria represent the most genetically diverse kingdom of life. While great progress has been made in describing this diversity, it remains difficult to identify the phylogenetic and ecological characteristics that delineate groups of bacteria that possess species-like properties. One major challenge associated with species delineations is that not all shared genes have the same evolutionary history, and thus the choice of loci can have a major impact on phylogenetic reconstruction. Sequencing the genomes of large numbers of closely related strains provides new opportunities to distinguish ancestral from acquired alleles and assess the effects of recombination on phylogenetic inference. Here we analyzed the genomes of 119 strains of the marine actinomycete genus Salinispora, which is currently comprised of three named species that share 99% 16S rRNA gene sequence identity. While 63% of the core genome showed evidence of recombination, this had no effect on species-level phylogenomic resolution. Recombination did however blur intra-species relationships and biogeographic resolution. The genome-wide average nucleotide identity provided a new perspective on Salinispora diversity, revealing as many as seven new species. Patterns of orthologous group distributions reveal a genetic basis to delineation the candidate taxa and insight into the levels of genetic cohesion associated with bacterial species.

Evolution and Distribution of C7-Cyclitol Synthases in Prokaryotes and Eukaryotes

2-Epi-5-epi-valiolone synthase (EEVS), a C₇-sugar phosphate cyclase (SPC) homologous to 3-dehydroquinate synthase (DHQS), was discovered during studies of the biosynthesis of the C₇N-aminocyclitol family of natural products. EEVS was originally thought to be present only in certain actinomycetes, but analyses of genome sequences showed that it is broadly distributed in both prokaryotes and eukaryotes, including vertebrates. Another SPC, desmethyl-4-deoxygadusol synthase (DDGS), was later discovered as being involved in the biosynthesis of mycosporine-like amino acid sunscreen compounds. Current database annotations are quite unreliable, with many EEVSs reported as DHQS, and most DDGSs reported as EEVS, DHQS, or simply hypothetical proteins. Here, we identify sequence features useful for distinguishing these enzymes, report a crystal structure of a representative DDGS showing the high similarity of the EEVS and DDGS enzymes, identify notable active site differences, and demonstrate the importance of two of these active site residues for catalysis by point mutations. Further, we functionally characterized two representatives of a distinct clade equidistant from known EEVS and known DDGS groups and show them to be authentic EEVSs. Moreover, we document and discuss the distribution of genes that encode EEVS and DDGS in various prokaryotes and eukaryotes, including pathogenic bacteria, plant symbionts, nitrogen-fixing bacteria, myxobacteria, cyanobacteria, fungi, stramenopiles, and animals, suggesting their broad potential biological roles in nature.

Recent advances in understanding Streptomyces

About 2,500 papers dated 2014–2016 were recovered by searching the PubMed database for Streptomyces, which are the richest known source of antibiotics. This review integrates around 100 of these papers in sections dealing with evolution, ecology, pathogenicity, growth and development, stress responses and secondary metabolism, gene expression, and technical advances. Genomic approaches have greatly accelerated progress. For example, it has been definitively shown that interspecies recombination of conserved genes has occurred during evolution, in addition to exchanges of some of the tens of thousands of non-conserved accessory genes. The closeness of the association of Streptomyces with plants, fungi, and insects has become clear and is reflected in the importance of regulators of cellulose and chitin utilisation in overall Streptomyces biology. Interestingly, endogenous cellulose-like glycans are also proving important in hyphal growth and in the clumping that affects industrial fermentations. Nucleotide secondary messengers, including cyclic di-GMP, have been shown to provide key input into developmental processes such as germination and reproductive growth, while late morphological changes during sporulation involve control by phosphorylation. The discovery that nitric oxide is produced endogenously puts a new face on speculative models in which regulatory Wbl proteins (peculiar to actinobacteria) respond to nitric oxide produced in stressful physiological transitions. Some dramatic insights have come from a new model system for Streptomyces developmental biology, Streptomyces venezuelae, including molecular evidence of very close interplay in each of two pairs of regulatory proteins. An extra dimension has been added to the many complexities of the regulation of secondary metabolism by findings of regulatory crosstalk within and between pathways, and even between species, mediated by end products. Among many outcomes from the application of chromosome immunoprecipitation sequencing (ChIP-seq) analysis and other methods based on “next-generation sequencing” has been the finding that 21% of Streptomyces mRNA species lack leader sequences and conventional ribosome binding sites. Further technical advances now emerging should lead to continued acceleration of knowledge, and more effective exploitation, of these astonishing and critically important organisms.