Intraspecies comparison of Streptomyces pratensis genomes reveals high levels of recombination and gene conservation between strains of disparate geographic origin

Background

Streptomyces are widespread bacteria that contribute to the terrestrial carbon cycle and produce the majority of clinically useful antibiotics. While interspecific genomic diversity has been investigated among Streptomyces, information is lacking on intraspecific genomic diversity. Streptomyces pratensis has high rates of homologous recombination but the impact of such gene exchange on genome evolution and the evolution of natural product gene clusters remains uncharacterized.

Results

We report draft genome sequences of four S. pratensis strains and compare to the complete genome of Streptomyces flavogriseus IAF-45-CD (=ATCC 33331), a strain recently reclassified to S. pratensis. Despite disparate geographic origins, the genomes are highly similar with 85.9% of genes present in the core genome and conservation of all natural product gene clusters. Natural products include a novel combination of carbapenem and beta-lactamase inhibitor gene clusters. While high intraspecies recombination rates abolish the phylogenetic signal across the genome, intraspecies recombination is suppressed in two genomic regions. The first region is centered on an insertion/deletion polymorphism and the second on a hybrid NRPS-PKS gene. Finally, two gene families accounted for over 25% of the divergent genes in the core genome. The first includes homologs of bldB (required for spore development and antibiotic production) while the second includes homologs of an uncharacterized protein with a helix-turn-helix motif (hpb). Genes from these families co-occur with fifteen pairs spread across the genome. These genes have evidence for co-evolution of co-localized pairs, supporting previous assertions that these genes may function akin to a toxin-antitoxin system.

Conclusions

S. pratensis genomes are highly similar with exceptional levels of recombination which erase phylogenetic signal among strains of the species. This species has a large core genome and variable terminal regions that are smaller than those found in interspecies comparisons. There is no geographic differentiation between these strains, but there is evidence for local linkage disequilibrium affecting two genomic regions. We have also shown further observational evidence that the DUF397-HTH (bldB and hpb) are a novel toxin-antitoxin pair.     

Genome mining of Streptomyces ambofaciens

Since the discovery of the streptomycin produced by Streptomyces griseus in the middle of the last century, members of this bacterial genus have been largely exploited for the production of secondary metabolites with wide uses in medicine and in agriculture. They have even been recognized as one of the most prolific producers of natural products among microorganisms. With the onset of the genomic era, it became evident that these microorganisms still represent a major source for the discovery of novel secondary metabolites. This was highlighted with the complete genome sequencing of Streptomyces coelicolor A3(2) which revealed an unexpected potential of this organism to synthesize natural products undetected until then by classical screening methods. Since then, analysis of sequenced genomes from numerous Streptomyces species has shown that a single species can carry more than 30 secondary metabolite gene clusters, reinforcing the idea that the biosynthetic potential of this bacterial genus is far from being fully exploited. This review highlights our knowledge on the potential of Streptomyces ambofaciens ATCC 23877 to synthesize natural products. This industrial strain was known for decades to only produce the drug spiramycin and another antibacterial compound, congocidine. Mining of its genome allowed the identification of 23 clusters potentially involved in the production of other secondary metabolites. Studies of some of these clusters resulted in the characterization of novel compounds and of previously known compounds but never characterized in this Streptomyces species. In addition, genome mining revealed that secondary metabolite gene clusters of phylogenetically closely related Streptomyces are mainly species-specific.     

Strain-Level Diversity of Secondary Metabolism in Streptomyces albus

Streptomyces spp. are robust producers of medicinally-, industrially- and agriculturally-important small molecules. Increased resistance to antibacterial agents and the lack of new antibiotics in the pipeline have led to a renaissance in natural product discovery. This endeavor has benefited from inexpensive high quality DNA sequencing technology, which has generated more than 140 genome sequences for taxonomic type strains and environmental Streptomyces spp. isolates. Many of the sequenced streptomycetes belong to the same species. For instance, Streptomyces albus has been isolated from diverse environmental niches and seven strains have been sequenced, consequently this species has been sequenced more than any other streptomycete, allowing valuable analyses of strain-level diversity in secondary metabolism. Bioinformatics analyses identified a total of 48 unique biosynthetic gene clusters harboured by Streptomyces albus strains. Eighteen of these gene clusters specify the core secondary metabolome of the species. Fourteen of the gene clusters are contained by one or more strain and are considered auxiliary, while 16 of the gene clusters encode the production of putative strain-specific secondary metabolites. Analysis of Streptomyces albus strains suggests that each strain of a Streptomyces species likely harbours at least one strain-specific biosynthetic gene cluster. Importantly, this implies that deep sequencing of a species will not exhaust gene cluster diversity and will continue to yield novelty.     

Pleiotropic Control of Secondary Metabolism and Morphological Development by KsbC,a Butyrolactone Autoregulator Receptor Homologue in Kitasatospora setae

The γ-butyrolactone autoregulator signaling cascades have been shown to control secondary metabolism and/or morphological development among many Streptomyces species. However, the conservation and variation of the regulatory systems among actinomycetes remain to be clarified. The genome sequence of Kitasatospora setae, which also belongs to the family Streptomycetaceae containing the genus Streptomyces, has revealed the presence of three homologues of the autoregulator receptor: KsbA, which has previously been confirmed to be involved only in secondary metabolism; KsbB; and KsbC. We describe here the characterization of ksbC, whose regulatory cluster closely resembles the Streptomyces virginiae barA locus responsible for the autoregulator signaling cascade. Deletion of the gene ksbC resulted in lowered production of bafilomycin and a defect of aerial mycelium formation, together with the early and enhanced production of a novel β-carboline alkaloid named kitasetaline. A putative kitasetaline biosynthetic gene cluster was identified, and its expression in a heterologous host led to the production of kitasetaline together with JBIR-133, the production of which is also detected in the ksbC disruptant, and JBIR-134 as novel β-carboline alkaloids, indicating that these genes were biosynthetic genes for β-carboline alkaloid and thus are the first such genes to be discovered in bacteria.     

Genome based analysis of type-I polyketide synthase and nonribosomal peptide synthetase gene clusters in seven strains of five representative Nocardia species

Background

Actinobacteria of the genus Nocardia usually live in soil or water and play saprophytic roles, but they also opportunistically infect the respiratory system, skin, and other organs of humans and animals. Primarily because of the clinical importance of the strains, some Nocardia genomes have been sequenced, and genome sequences have accumulated. Genome sizes of Nocardia strains are similar to those of Streptomyces strains, the producers of most antibiotics. In the present work, we compared secondary metabolite biosynthesis gene clusters of type-I polyketide synthase (PKS-I) and nonribosomal peptide synthetase (NRPS) among genomes of representative Nocardia species/strains based on domain organization and amino acid sequence homology.

Results

Draft genome sequences of Nocardia asteroides NBRC 15531^T, Nocardia otitidiscaviarum IFM 11049, Nocardia brasiliensis NBRC 14402^T, and N. brasiliensis IFM 10847 were read and compared with published complete genome sequences of Nocardia farcinica IFM 10152, Nocardia cyriacigeorgica GUH-2, and N. brasiliensis HUJEG-1. Genome sizes are as follows: N. farcinica, 6.0 Mb; N. cyriacigeorgica, 6.2 Mb; N. asteroides, 7.0 Mb; N. otitidiscaviarum, 7.8 Mb; and N. brasiliensis, 8.9 - 9.4 Mb. Predicted numbers of PKS-I, NRPS, and PKS-I/NRPS hybrid clusters ranged between 4–11, 7–13, and 1–6, respectively, depending on strains, and tended to increase with increasing genome size. Domain and module structures of representative or unique clusters are discussed in the text.

Conclusion

We conclude the following: 1) genomes of Nocardia strains carry as many PKS-I and NRPS gene clusters as those of Streptomyces strains, 2) the number of PKS-I and NRPS gene clusters in Nocardia strains varies substantially depending on species, and N. brasiliensis strains carry the largest numbers of clusters among the species studied, 3) the seven Nocardia strains studied in the present work have seven common PKS-I and/or NRPS clusters, some of whose products are yet to be studied, and 4) different N. brasiliensis strains have some different gene clusters of PKS-I/NRPS, although the rest of the clusters are common within the N. brasiliensis strains. Genome sequencing suggested that Nocardia strains are highly promising resources in the search of novel secondary metabolites.

Electronic supplementary material

The online version of this article (doi: 10.1186/1471-2164-15-323) contains supplementary material, which is available to authorized users.     

Genome Sequence of Kitasatospora setae NBRC 14216T: An Evolutionary Snapshot of the Family Streptomycetaceae

Kitasatospora setae NBRC 14216^T (=KM-6054^T) is known to produce setamycin (bafilomycin B1) possessing antitrichomonal activity. The genus Kitasatospora is morphologically similar to the genus Streptomyces, although they are distinguishable from each other on the basis of cell wall composition and the 16S rDNA sequence. We have determined the complete genome sequence of K. setae NBRC 14216^T as the first Streptomycetaceae genome other than Streptomyces. The genome is a single linear chromosome of 8 783 278 bp with terminal inverted repeats of 127 148 bp, predicted to encode 7569 protein-coding genes, 9 rRNA operons, 1 tmRNA and 74 tRNA genes. Although these features resemble those of Streptomyces, genome-wide comparison of orthologous genes between K. setae and Streptomyces revealed smaller extent of synteny. Multilocus phylogenetic analysis based on amino acid sequences unequivocally placed K. setae outside the Streptomyces genus. Although many of the genes related to morphological differentiation identified in Streptomyces were highly conserved in K. setae, there were some differences such as the apparent absence of the AmfS (SapB) class of surfactant protein and differences in the copy number and variation of paralogous components involved in cell wall synthesis.     

Identification of MicroRNAs in Helicoverpa armigera and Spodoptera litura Based on Deep Sequencing and Homology Analysis

The current identification of microRNAs (miRNAs) in insects is largely dependent on genome sequences. However, the lack of available genome sequences inhibits the identification of miRNAs in various insect species. In this study, we used a miRNA database of the silkworm Bombyx mori as a reference to identify miRNAs in Helicoverpa armigera and Spodoptera litura using deep sequencing and homology analysis. Because all three species belong to the Lepidoptera, the experiment produced reliable results. Our study identified 97 and 91 conserved miRNAs in H. armigera and S. litura, respectively. Using the genome of B. mori and BAC sequences of H. armigera as references, 1 novel miRNA and 8 novel miRNA candidates were identified in H. armigera, and 4 novel miRNA candidates were identified in S. litura. An evolutionary analysis revealed that most of the identified miRNAs were insect-specific, and more than 20 miRNAs were Lepidoptera-specific. The investigation of the expression patterns of miR-2a, miR-34, miR-2796-3p and miR-11 revealed their potential roles in insect development. miRNA target prediction revealed that conserved miRNA target sites exist in various genes in the 3 species. Conserved miRNA target sites for the Hsp90 gene among the 3 species were validated in the mammalian 293T cell line using a dual-luciferase reporter assay. Our study provides a new approach with which to identify miRNAs in insects lacking genome information and contributes to the functional analysis of insect miRNAs.     

Complete Mitochondrial Genomes of Carcinoscorpius rotundicauda and Tachypleus tridentatus (Xiphosura,Arthropoda) and Implications for Chelicerate Phylogenetic Studies

Horseshoe crabs (order Xiphosura) are often referred to as an ancient order of marine chelicerates and have been considered as keystone taxa for the understanding of chelicerate evolution. However, the mitochondrial genome of this order is only available from a single species, Limulus polyphemus. In the present study, we analyzed the complete mitochondrial genomes from two Asian horseshoe crabs, Carcinoscorpius rotundicauda and Tachypleus tridentatus to offer novel data for the evolutionary relationship within Xiphosura and their position in the chelicerate phylogeny. The mitochondrial genomes of C. rotundicauda (15,033 bp) and T. tridentatus (15,006 bp) encode 13 protein-coding genes, two ribosomal RNA (rRNA) genes, and 22 transfer RNA (tRNA) genes. Overall sequences and genome structure of two Asian species were highly similar to that of Limulus polyphemus, though clear differences among three were found in the stem-loop structure of the putative control region. In the phylogenetic analysis with complete mitochondrial genomes of 43 chelicerate species, C. rotundicauda and T. tridentatus were recovered as a monophyly, while L. polyphemus solely formed an independent clade. Xiphosuran species were placed at the basal root of the tree, and major other chelicerate taxa were clustered in a single monophyly, clearly confirming that horseshoe crabs composed an ancestral taxon among chelicerates. By contrast, the phylogenetic tree without the information of Asian horseshoe crabs did not support monophyletic clustering of other chelicerates. In conclusion, our analyses may provide more robust and reliable perspective on the study of evolutionary history for chelicerates than earlier analyses with a single Atlantic species.     

Dissecting the fungal biology of Bipolaris papendorfii: from phylogenetic to comparative genomic analysis

Bipolaris papendorfii has been reported as a fungal plant pathogen that rarely causes opportunistic infection in humans. Secondary metabolites isolated from this fungus possess medicinal and anticancer properties. However, its genetic fundamental and basic biology are largely unknown. In this study, we report the first draft genome sequence of B. papendorfii UM 226 isolated from the skin scraping of a patient. The assembled 33.4 Mb genome encodes 11,015 putative coding DNA sequences, of which, 2.49% are predicted transposable elements. Multilocus phylogenetic and phylogenomic analyses showed B. papendorfii UM 226 clustering with Curvularia species, apart from other plant pathogenic Bipolaris species. Its genomic features suggest that it is a heterothallic fungus with a putative unique gene encoding the LysM-containing protein which might be involved in fungal virulence on host plants, as well as a wide array of enzymes involved in carbohydrate metabolism, degradation of polysaccharides and lignin in the plant cell wall, secondary metabolite biosynthesis (including dimethylallyl tryptophan synthase, non-ribosomal peptide synthetase, polyketide synthase), the terpenoid pathway and the caffeine metabolism. This first genomic characterization of B. papendorfii provides the basis for further studies on its biology, pathogenicity and medicinal potential.     

A novel spliceosome-mediated trans-splicing can change our view on genome complexity of the divergent eukaryote Giardia intestinalis

Although spliceosomal introns are an abundant landmark in eukaryotic genomes, the nuclear genome of the divergent eukaryote Giardia intestinalis, the causative agent of giardiasis, has been considered as “intron-poor” with only five canonical (cis-spliced) introns. However, three research groups (including ours) have independently reported a novel class of spliceosomal introns in the G. intestinalis genome. Three protein-coding genes are split into pieces in the G. intestinalis genome, and each of the partial coding regions was independently transcribed into polyadenylated premature mRNAs (pre-mRNAs). The two pre-mRNAs directly interact with each other by an intermolecular-stem structure formed between their non-coding portions, and are then processed into mature mRNAs by spliceosome-mediated trans-splicing. Here, we summarize the recently published works on split introns (“splintrons”) in the G. intestinalis genome, and then provide our speculation on the functional property of the Giardia spliceosomes based on the putative ratio of splintrons to canonical introns. Finally, we discuss a scenario for the transition from typical GT-AG boundaries to non-typical AT-AC boundaries in a particular splintron of Giardia.     

Systematic investigation of Amphioxus (Branchiostoma floridae) microRNAs

Ascaris lumbricoides and Ascaris suum are parasitic nematodes living in the small intestine of humans and pigs, and can cause the disease ascariasis. For long, there has been controversy as to whether the two ascaridoid taxa represent the same species due to their significant resemblances in morphology. However, the complete mitochondrial (mt) genome data have been lacking for A. lumbricoides in spite of human and animal health significance and socio-economic impact globally of these parasites. In the present study, we sequenced the complete mt genomes of A. lumbricoides and A. suum (China isolate), which was 14,303 bp and 14,311 bp in size, respectively. The identity of the mt genomes was 98.1% between A. lumbricoides and A. suum (China isolate), and 98.5% between A. suum (China isolate) and A. suum (USA isolate). Both genomes are circular, and consist of 36 genes, including 12 genes for proteins, 2 genes for rRNA and 22 genes for tRNA, which are consistent with that of all other species of ascaridoid studied to date. All genes are transcribed in the same direction and have a nucleotide composition high in A and T (71.7% for A. lumbricoides and 71.8% for A. suum). The AT bias had a significant effect on both the codon usage pattern and amino acid composition of proteins. Phylogenetic analyses of A. lumbricoides and A. suum using concatenated amino acid sequences of 12 protein-coding genes, with three different computational algorithms (Bayesian analysis, maximum likelihood and maximum parsimony) all clustered in a clade with high statistical support, indicating that A. lumbricoides and A. suum was very closely related. These mt genome data and the results provide some additional genetic evidence that A. lumbricoides and A. suum may represent the same species. The mt genome data presented in this study are also useful novel markers for studying the molecular epidemiology and population genetics of Ascaris.     

Cellulose-degrading bacteria associated with the invasive woodwasp Sirex noctilio

Sirex noctilio is an invasive wood-feeding wasp that threatens the world''s commercial and natural pine forests. Successful tree colonization by this insect is contingent on the decline of host defenses and the ability to utilize the woody substrate as a source of energy. We explored its potential association with bacterial symbionts that may assist in nutrient acquisition via plant biomass deconstruction using growth assays, culture-dependent and -independent analysis of bacterial frequency of association and whole-genome analysis. We identified Streptomyces and γ-Proteobacteria that were each associated with 94% and 88% of wasps, respectively. Streptomyces isolates grew on all three cellulose substrates tested and across a range of pH 5.6 to 9. On the basis of whole-genome sequencing, three Streptomyces isolates have some of the highest proportions of genes predicted to encode for carbohydrate-active enzymes (CAZyme) of sequenced Actinobacteria. γ-Proteobacteria isolates grew on a cellulose derivative and a structurally diverse substrate, ammonia fiber explosion-treated corn stover, but not on microcrystalline cellulose. Analysis of the genome of a Pantoea isolate detected genes putatively encoding for CAZymes, the majority predicted to be active on hemicellulose and more simple sugars. We propose that a consortium of microorganisms, including the described bacteria and the fungal symbiont Amylostereum areolatum, has complementary functions for degrading woody substrates and that such degradation may assist in nutrient acquisition by S. noctilio, thus contributing to its ability to be established in forested habitats worldwide.     

Editor's choice: The complete genome sequencing of Prevotella intermedia strain OMA14 and a subsequent fine-scale,intra-species genomic comparison reveal an unusual amplification of conjugative and mobile transposons and identify a novel Prevotella-lineage-specific repeat

Prevotella intermedia is a pathogenic bacterium involved in periodontal diseases. Here, we present the complete genome sequence of a clinical strain, OMA14, of this bacterium along with the results of comparative genome analysis with strain 17 of the same species whose genome has also been sequenced, but not fully analysed yet. The genomes of both strains consist of two circular chromosomes: the larger chromosomes are similar in size and exhibit a high overall linearity of gene organizations, whereas the smaller chromosomes show a significant size variation and have undergone remarkable genome rearrangements. Unique features of the Pre. intermedia genomes are the presence of a remarkable number of essential genes on the second chromosomes and the abundance of conjugative and mobilizable transposons (CTns and MTns). The CTns/MTns are particularly abundant in the second chromosomes, involved in its extensive genome rearrangement, and have introduced a number of strain-specific genes into each strain. We also found a novel 188-bp repeat sequence that has been highly amplified in Pre. intermedia and are specifically distributed among the Pre. intermedia-related species. These findings expand our understanding of the genetic features of Pre. intermedia and the roles of CTns and MTns in the evolution of bacteria.     

Genome mining of the Streptomyces avermitilis genome and development of genome-minimized hosts for heterologous expression of biosynthetic gene clusters

To date, several actinomycete genomes have been completed and annotated. Among them, Streptomyces microorganisms are of major pharmaceutical interest because they are a rich source of numerous secondary metabolites. S. avermitilis is an industrial microorganism used for the production of an anthelmintic agent, avermectin, which is a commercially important antiparasitic agent in human and veterinary medicine, and agricultural pesticides. Genome analysis of S. avermitilis provides significant information for not only industrial applications but also understanding the features of this genus. On genome mining of S. avermitilis, the microorganism has been found to harbor at least 38 secondary metabolic gene clusters and 46 insertion sequence (IS)-like sequences on the genome, which have not been searched so far. A significant use of the genome data of Streptomyces microorganisms is the construction of a versatile host for heterologous expression of exogenous biosynthetic gene clusters by genetic engineering. Since S. avermitilis is used as an industrial microorganism, the microorganism is already optimized for the efficient supply of primary metabolic precursors and biochemical energy to support multistep biosynthesis. The feasibility of large-deletion mutants of S. avermitilis has been confirmed by heterologous expression of more than 20 exogenous biosynthetic gene clusters.     

Extensive Chromosome Homoeology among Brassiceae Species Were Revealed by Comparative Genetic Mapping with High-Density EST-Based SNP Markers in Radish (Raphanus sativus L.)

A linkage map of expressed sequence tag (EST)-based markers in radish (Raphanus sativus L.) was constructed using a low-cost and high-efficiency single-nucleotide polymorphism (SNP) genotyping method named multiplex polymerase chain reaction–mixed probe dot-blot analysis developed in this study. Seven hundred and forty-six SNP markers derived from EST sequences of R. sativus were assigned to nine linkage groups with a total length of 806.7 cM. By BLASTN, 726 markers were found to have homologous genes in Arabidopsis thaliana, and 72 syntenic regions, which have great potential for utilizing genomic information of the model species A. thaliana in basic and applied genetics of R. sativus, were identified. By construction and analysis of the genome structures of R. sativus based on the 24 genomic blocks within the Brassicaceae ancestral karyotype, 23 of the 24 genomic blocks were detected in the genome of R. sativus, and half of them were found to be triplicated. Comparison of the genome structure of R. sativus with those of the A, B, and C genomes of Brassica species and that of Sinapis alba L. revealed extensive chromosome homoeology among Brassiceae species, which would facilitate transfer of the genomic information from one Brassiceae species to another.     

High quality draft genome sequence of Streptomyces sp. strain AW19M42 isolated from a sea squirt in Northern Norway

Here we report the 8 Mb high quality draft genome of Streptomyces sp. strain AW19M42, together with specific properties of the organism and the generation, annotation and analysis of its genome sequence. The genome encodes 7,727 putative open reading frames, of which 6,400 could be assigned with COG categories. Also, 62 tRNA genes and 8 rRNA operons were identified. The genome harbors several gene clusters involved in the production of secondary metabolites. Functional screening of the isolate was positive for several enzymatic activities, and some candidate genes coding for those activities are listed in this report. We find that this isolate shows biotechnological potential and is an interesting target for bioprospecting.     

大连海域沉积物中可培养海洋链霉菌活性及其多样性

以大肠杆菌、金黄色葡萄球菌和尖孢镰刀枯萎病菌作为测试靶目标,采用9种分离培养基从大连海域13个不同采样点的海洋沉积物样品中分离到165株海洋链霉菌。从165株海洋放线菌中筛选到对金黄色葡萄球菌具有抑制活性的菌株85株,占总菌株数的51.5%;对大肠杆菌具有抑制活性的菌株27株,占总菌株数的16.4%;对尖孢镰刀枯萎病菌具有抑制活性的菌株仅有6株,占总菌株数的3.6%。因此,海洋链霉菌的活性更多地表现为对细菌的抗性,尤其对革兰氏阳性细菌具有更高的抑制活性。对其中具有抑制活性或形态独特的菌株进行了16S r DNA序列分析,并构建系统发育树,显示活性海洋链霉菌具有丰富的种类多样性和广谱抗菌活性。同种海洋链霉菌与土壤链霉菌活性比较结果也表明,海洋链霉菌多表现抗革兰氏阳性细菌活性。     

Evolutionary Relationships among Actinophages and a Putative Adaptation for Growth in Streptomyces spp.

The genome sequences of eight Streptomyces phages are presented, four of which were isolated for this study. Phages R4, TG1, ϕHau3, and SV1 were isolated previously and have been exploited as tools for understanding and genetically manipulating Streptomyces spp. We also extracted five apparently intact prophages from recent Streptomyces spp. genome projects and, together with six phage genomes in the database, we analyzed all 19 Streptomyces phage genomes with a view to understanding their relationships to each other and to other actinophages, particularly the mycobacteriophages. Fifteen of the Streptomyces phages group into four clusters of related genomes. Although the R4-like phages do not share nucleotide sequence similarity with other phages, they clearly have common ancestry with cluster A mycobacteriophages, sharing many protein homologues, common gene syntenies, and similar repressor-stoperator regulatory systems. The R4-like phage ϕHau3 and the prophage StrepC.1 (from Streptomyces sp. strain C) appear to have hijacked a unique adaptation of the streptomycetes, i.e., use of the rare UUA codon, to control translation of the essential phage protein, the terminase. The Streptomyces venezuelae generalized transducing phage SV1 was used to predict the presence of other generalized transducing phages for different Streptomyces species.