Genome-based phylogenetic analysis of Streptomyces and its relatives

16S rRNA phylogeny and clustering is not a reliable proxy for genome-based taxonomy in Streptomyces

Streptomyces is among the most extensively studied genera of bacteria but its complex taxonomy remains contested and is suspected to contain significant species-level misclassification. Resolving the classification of Streptomyces would benefit many areas of applied microbiology that rely on an accurate ground truth for grouping of related organisms, including comparative genomics-based searches for novel antimicrobials. We survey taxonomic conflicts between 16S rRNA and whole genome-based Streptomyces classifications using 2276 publicly available Streptomyces genome assemblies and 48 981 publicly available full-length 16S rRNA Streptomyces sequences from silva, Greengenes, Ribosomal Database Project (RDP), and NCBI (National Centre for Biotechnology Information) databases. We construct a full-length 16S gene tree for 14 239 distinct Streptomyces sequences that resolves three major lineages of Streptomyces, but whose topology is not consistent with existing taxonomic assignments. We use these sequence data to delineate 16S and whole genome landscapes for Streptomyces, demonstrating that 16S and whole-genome classifications are frequently in disagreement, and that 16S zero-radius Operational Taxonomic Units (zOTUs) are often inconsistent with Average Nucleotide Identity (ANI)-based taxonomy. Our results strongly imply that 16S rRNA sequence data does not map to taxonomy sufficiently well to delineate Streptomyces species routinely. We propose that alternative marker sequences should be adopted by the community for classification and metabarcoding. Insofar as Streptomyces taxonomy has been determined or supported by 16S sequence data and may in parts be in error, we also propose that reclassification of the genus by alternative approaches may benefit the Streptomyces community.

Present Status, Limitations, and Prospects of Using Streptomyces Bacteria as a Potential Probiotic Agent in Aquaculture

Streptomyces is a Gram-positive bacterium, belonging to the family Streptomycetaceae and order Streptomycetales. Several strains from different species of Streptomyces can be used to promote the health and growth of artificially cultured fish and shellfish by producing secondary metabolites including antibiotics, anticancer agents, antiparasitic agents, antifungal agents, and enzymes (protease and amylase). Some Streptomyces strains also exhibit antagonistic and antimicrobial activity against aquaculture-based pathogens by producing inhibitory compounds such as bacteriocins, siderophores, hydrogen peroxide, and organic acids to compete for nutrients and attachment sites in the host. The administration of Streptomyces in aquaculture could also induce an immune response, disease resistance, quorum sensing/antibiofilm activity, antiviral activity, competitive exclusion, modification in gastrointestinal microflora, growth enhancement, and water quality amelioration via nitrogen fixation and degradation of organic residues from the culture system. This review provides the current status and prospects of Streptomyces as potential probiotics in aquaculture, their selection criteria, administrative methods, and mechanisms of action. The limitations of Streptomyces as probiotics in aquaculture are highlighted and the solutions to these limitations are also discussed.

Characterization of the mitochondrial genome of Chlorolobion braunii ITBB-AG6, an azolla-associated green alga isolated from sanitary sewage

Sphaeropleales have the characteristics of rapid growth, high oil content, and efficient removal rates of nitrogen and phosphorus in sewage waters, and is potentially valuable in biodiesel production and environmental remediation. In this study, we isolated a strain of Sphaeropleales, Chlorolobion braunii strain ITBB-AG6 from an azolla community in a sewage pond. Its mitochondrial genome contains 110,124 bp and harbors at least 40 genes, including 15 protein-coding genes, 20 tRNA genes, and three rRNA genes. The protein-coding genes include two for ATP synthases, seven for NAD(P)H-quinone oxidoreductases (nad), three for cytochrome c oxidase subunits (coxs), and one for cytochrome b (cob). Transfer RNA genes for 18 amino acids were identified, in which the tRNA genes for leucine and serine are doubled, but the tRNA genes for threonine and valine are not annotated. Phylogenetic analysis using the mitochondrial genomes of seven families of Sphaeropleales indicated that ITBB-AG6 is closely related to Monoraphidium neglectum, and falls in the family Selenastraceae with 100% bootstrap support. Two species in the family Neochloridaceae are separated by a species in Hydrodictyaceae, indicating a polyphyletic nature. These findings revealed the complicated phylogenetic relationships of the Sphaeropleales and the necessity of genome sequences in the taxonomy of microalgae.

A panoramic view of the genomic landscape of the genus Streptomyces

We delineate the evolutionary plasticity of the ecologically and biotechnologically important genus Streptomyces, by analysing the genomes of 213 species. Streptomycetes genomes demonstrate high levels of internal homology, whereas the genome of their last common ancestor was already complex. Importantly, we identify the species-specific fingerprint proteins that characterize each species. Even among closely related species, we observed high interspecies variability of chromosomal protein-coding genes, species-level core genes, accessory genes and fingerprints. Notably, secondary metabolite biosynthetic gene clusters (smBGCs), carbohydrate-active enzymes (CAZymes) and protein-coding genes bearing the rare TTA codon demonstrate high intraspecies and interspecies variability, which emphasizes the need for strain-specific genomic mining. Highly conserved genes, such as those specifying genus-level core proteins, tend to occur in the central region of the chromosome, whereas those encoding proteins with evolutionarily volatile species-level fingerprints, smBGCs, CAZymes and TTA-codon-bearing genes are often found towards the ends of the linear chromosome. Thus, the chromosomal arms emerge as the part of the genome that is mainly responsible for rapid adaptation at the species and strain level. Finally, we observed a moderate, but statistically significant, correlation between the total number of CAZymes and three categories of smBGCs (siderophores, e-Polylysin and type III lanthipeptides) that are related to competition among bacteria.

Plant catalase in silico characterization and phylogenetic analysis with structural modeling

Background

Catalase (EC 1.11.1.6) is a heme-containing tetrameric enzyme that plays a critical role in signaling and hydrogen peroxide metabolism. It was the first enzyme to be crystallized and isolated. Catalase is a well-known industrial enzyme used in diagnostic and analytical methods in the form of biomarkers and biosensors, as well as in the textile, paper, food, and pharmaceutical industries. In silico analysis of CAT genes and proteins has gained increased interest, emphasizing the development of biomarkers and drug designs. The present work aims to understand the catalase evolutionary relationship of plant species and analyze its physicochemical characteristics, homology, phylogenetic tree construction, secondary structure prediction, and 3D modeling of protein sequences and its validation using a variety of conventional computational methods to assist researchers in better understanding the structure of proteins.

Results

Around 65 plant catalase sequences were computationally evaluated and subjected to bioinformatics assessment for physicochemical characterization, multiple sequence alignment, phylogenetic construction, motif and domain identification, and secondary and tertiary structure prediction. The phylogenetic tree revealed six unique clusters where diversity of plant catalases was found to be the largest for Oryza sativa. The thermostability and hydrophilic nature of these proteins were primarily observed, as evidenced by a relatively high aliphatic index and negative GRAVY value. The distribution of 5 sequence motifs was uniformly distributed with a width length of 50 with the best possible amino residue sequences that resemble the plant catalase PLN02609 superfamily. Using SOPMA, the predicted secondary structure of the protein sequences revealed the predominance of the random coil. The predicted 3D CAT model from Arabidopsis thaliana was a homotetramer, thermostable protein with 59-KDa weight, and its structural validation was confirmed by PROCHECK, ERRAT, Verify3D, and Ramachandran plot. The functional relationships of our query sequence revealed the glutathione reductase as the closest interacting protein of query protein.

Conclusions

This theoretical plant catalases in silico analysis provide insight into its physiochemical characteristics and functional and structural understanding and its evolutionary behavior and exploring protein structure-function relationships when crystal structures are unavailable.

Revealing Genome-Based Biosynthetic Potential of Streptomyces sp. BR123 Isolated from Sunflower Rhizosphere with Broad Spectrum Antimicrobial Activity

Actinomycetes, most notably the genus Streptomyces, have great importance due to their role in the discovery of new natural products, especially for finding antimicrobial secondary metabolites that are useful in the medicinal science and biotechnology industries. In the current study, a genome-based evaluation of Streptomyces sp. isolate BR123 was analyzed to determine its biosynthetic potential, based on its in vitro antimicrobial activity against a broad range of microbial pathogens, including gram-positive and gram-negative bacteria and fungi. A draft genome sequence of 8.15 Mb of Streptomyces sp. isolate BR123 was attained, containing a GC content of 72.63% and 8103 protein coding genes. Many antimicrobial, antiparasitic, and anticancerous compounds were detected by the presence of multiple biosynthetic gene clusters, which was predicted by in silico analysis. A novel metabolite with a molecular mass of 1271.7773 in positive ion mode was detected through a high-performance liquid chromatography linked with mass spectrometry (HPLC-MS) analysis. In addition, another compound, meridamycin, was also identified through a HPLC-MS analysis. The current study reveals the biosynthetic potential of Streptomyces sp. isolate BR123, with respect to the synthesis of bioactive secondary metabolites through genomic and spectrometric analysis. Moreover, the comparative genome study compared the isolate BR123 with other Streptomyces strains, which may expand the knowledge concerning the mechanism involved in novel antimicrobial metabolite synthesis.

Whole genome sequencing of Streptomyces actuosus ISP-5337, Streptomyces sioyaensis B-5408, and Actinospica acidiphila B-2296 reveals secondary metabolomes with antibiotic potential

•

Whole genome sequencing of Actinomycetes reveals metabolic potential.

•

High quality genomes are necessary for mining of biosynthetic gene clusters.

•

Characterization of thiopeptides by high resolution mass spectrometry.

•

Thiopeptides are potent antibacterials against Staphylococcus aureus.

Streptomycetes are bacteria of biotechnological importance since they are avid producers of secondary metabolites, including antibiotics. Progress in genome mining has recently shown that Streptomyces species encode for many biosynthetic gene clusters which are mostly unexplored. Here, we selected three Actinomycetes species for whole genome sequencing that are known to produce potent thiopeptide antibiotics. Streptomyces actuosus biosynthesizes nosiheptide, Streptomyces sioyaensis produces siomycin, and Actinospica acidiphila is a member of the Actinomycete subfamily. Bioinformatic analyses demonstrated diverse secondary metabolomes with multiple antibiotic-encoding gene clusters. Detailed mass spectrometry analysis of metabolite extracts verified the active expression of nosiheptide and siomycin from S. actuosus and S. sioyaensis while fractionation of the bacterial extracts and subsequent challenge against Staphylococcus aureus demonstrated potent antibiotic activity of fractions containing these compounds. Whole genome sequencing of these species facilitates future bioengineering efforts for thiopeptides and characterization of relevant secondary metabolites.

Characterization of the mitochondrion genome of a Chlorella vulgaris strain isolated from rubber processing wastewater

Chlorella vulgaris ITBBA3-12 was isolated from the rubber processing wastewater and has a role in wastewater purification. Its complete mitogenome contains 88754 bp, with a G + C content of 29.7%. A total of 64 genes were annotated, including 34 protein-coding genes, 27 tRNA genes, three rRNA (rrn23, rrn16, and rrn5). Phylogenetic analysis using the mitogenomes of Trebouxiophyceae species indicated that the strain ITBBA3-12 is closely related to C. vulgaris strain UTEX259 and NJ-7, and they clustered in the Chlorella lineage.

Comparative Genomics among Closely Related Streptomyces Strains Revealed Specialized Metabolite Biosynthetic Gene Cluster Diversity

Specialized metabolites are of great interest due to their possible industrial and clinical applications. The increasing number of antimicrobial resistant infectious agents is a major health threat and therefore, the discovery of chemical diversity and new antimicrobials is crucial. Extensive genomic data from Streptomyces spp. confirm their production potential and great importance. Genome sequencing of the same species strains indicates that specialized metabolite biosynthetic gene cluster (SMBGC) diversity is not exhausted, and instead, a pool of novel specialized metabolites still exists. Here, we analyze the genome sequence data from six phylogenetically close Streptomyces strains. The results reveal that the closer strains are phylogenetically, the number of shared gene clusters is higher. Eight specialized metabolites comprise the core metabolome, although some strains have only six core gene clusters. The number of conserved gene clusters common between the isolated strains and their closest phylogenetic counterparts varies from nine to 23 SMBGCs. However, the analysis of these phylogenetic relationships is not affected by the acquisition of gene clusters, probably by horizontal gene transfer events, as each strain also harbors strain-specific SMBGCs. Between one and 15 strain-specific gene clusters were identified, of which up to six gene clusters in a single strain are unknown and have no identifiable orthologs in other species, attesting to the existing SMBGC novelty at the strain level.

Analysis of metabolic networks of Streptomyces leeuwenhoekii C34 by means of a genome scale model: Prediction of modifications that enhance the production of specialized metabolites

The first genome scale model (GSM) for Streptomyces leeuwenhoekii C34 was developed to study the biosynthesis pathways of specialized metabolites and to find metabolic engineering targets for enhancing their production. The model, iVR1007, consists of 1,722 reactions, 1,463 metabolites, and 1,007 genes, it includes the biosynthesis pathways of chaxamycins, chaxalactins, desferrioxamines, ectoine, and other specialized metabolites. iVR1007 was validated using experimental information of growth on 166 different sources of carbon, nitrogen and phosphorous, showing an 83.7% accuracy. The model was used to predict metabolic engineering targets for enhancing the biosynthesis of chaxamycins and chaxalactins. Gene knockouts, such as sle03600 (L-homoserine O-acetyltransferase), and sle39090 (trehalose-phosphate synthase), that enhance the production of the specialized metabolites by increasing the pool of precursors were identified. Using the algorithm of flux scanning based on enforced objective flux (FSEOF) implemented in python, 35 and 25 over-expression targets for increasing the production of chaxamycin A and chaxalactin A, respectively, that were not directly associated with their biosynthesis routes were identified. Nineteen over-expression targets that were common to the two specialized metabolites studied, like the over-expression of the acetyl carboxylase complex (sle47660 (accA) and any of the following genes: sle44630 (accA_1) or sle39830 (accA_2) or sle27560 (bccA) or sle59710) were identified. The predicted knockouts and over-expression targets will be used to perform metabolic engineering of S. leeuwenhoekii C34 and obtain overproducer strains.

A systematic reconstruction and constraint-based analysis of Leishmania donovani metabolic network: identification of potential antileishmanial drug targets

Visceral leishmaniasis, a lethal parasitic disease, is caused by the protozoan parasite Leishmania donovani. The absence of an effective vaccine, drug toxicity and parasite resistance necessitates the identification of novel drug targets. Reconstruction of genome-scale metabolic models and their simulation has been established as an important tool for systems-level understanding of a microorganism's metabolism. In this work, amalgamating the tools and techniques of computational systems biology with rigorous manual curation, a constraint-based metabolic model for Leishmania donovani BPK282A1 has been developed. New functional annotations for 18 formerly hypothetical or erroneously annotated genes (encountered during iterative refinement of the model) have been proposed. Further, to formulate an accurate biomass objective function, experimental determination of previously uncharacterized biomass constituents was performed. The developed model is a highly compartmentalized metabolic model, comprising 1159 reactions, 1135 metabolites and 604 genes. The model exhibited around 76% accuracy for the prediction of experimental phenotypes of gene knockout studies and drug inhibition assays. Employing in silico gene knockout studies, we identified 28 essential genes with negligible sequence identity to the human proteins. Moreover, by dissecting the functional interdependencies of metabolic pathways, 70 synthetic lethal pairs were identified. Finally, in order to delineate stage-specific metabolism, gene-expression data of the amastigote stage residing in human macrophages were integrated into the model. By comparing the flux distribution, we illustrated the stage-specific differences in metabolism and environmental conditions that are in good agreement with the experimental findings. The developed model can serve as a highly enriched knowledgebase of legacy data and an important tool for generating experimentally verifiable hypotheses.

Promiscuous Pathogenicity Islands and Phylogeny of Pathogenic Streptomyces spp

Approximately 10 Streptomyces species cause disease on underground plant structures. The most economically important of these is potato scab, and the most studied of these pathogens is Streptomyces scabiei (syn. S. scabies). The main pathogenicity determinant of scab-causing Streptomyces species is a nitrated diketopiperazine, known as thaxtomin A (ThxA). In the pathogenic species Streptomyces turgidiscabies, ThxA biosynthetic genes reside on a mobile pathogenicity island (PAI). However, the mobilization of PAIs in other Streptomyces species remains uncharacterized. Here, we investigated the mobilization of the PAI of S. scabiei 87-22. Based on whole genome sequences, we inferred the evolutionary relationships of pathogenic Streptomyces species and discovered that Streptomyces sp. strain 96-12, a novel pathogenic species isolated from potatoes in Egypt, was phylogenetically grouped with nonpathogenic species rather than with known pathogenic species. We also found that Streptomyces sp. strain 96-12 contains a PAI that is almost identical to the PAI in S. scabiei 87-22, despite significant differences in their genome sequences. This suggested direct or indirect in vivo mobilization of the PAI between S. scabiei and nonpathogenic Streptomyces species. To test whether the S. scabiei 87-22 PAI could, indeed, be mobilized, S. scabiei 87-22 deletion mutants containing antibiotic resistance markers in the PAI were mated with Streptomyces diastatochromogenes, a nonpathogenic species. The PAI of S. scabiei was site-specifically inserted into the aviX1 gene of S. diastatochromogenes and conferred pathogenicity in radish seedling assays. Our results demonstrated that S. scabiei, the earliest described Streptomyces pathogen, could be the source of a PAI responsible for the emergence of novel pathogenic species.

Developmental biology of Streptomyces from the perspective of 100 actinobacterial genome sequences

To illuminate the evolution and mechanisms of actinobacterial complexity, we evaluate the distribution and origins of known Streptomyces developmental genes and the developmental significance of actinobacteria-specific genes. As an aid, we developed the Actinoblast database of reciprocal blastp best hits between the Streptomyces coelicolor genome and more than 100 other actinobacterial genomes (http://streptomyces.org.uk/actinoblast/). We suggest that the emergence of morphological complexity was underpinned by special features of early actinobacteria, such as polar growth and the coupled participation of regulatory Wbl proteins and the redox-protecting thiol mycothiol in transducing a transient nitric oxide signal generated during physiologically stressful growth transitions. It seems that some cell growth and division proteins of early actinobacteria have acquired greater importance for sporulation of complex actinobacteria than for mycelial growth, in which septa are infrequent and not associated with complete cell separation. The acquisition of extracellular proteins with structural roles, a highly regulated extracellular protease cascade, and additional regulatory genes allowed early actinobacterial stationary phase processes to be redeployed in the emergence of aerial hyphae from mycelial mats and in the formation of spore chains. These extracellular proteins may have contributed to speciation. Simpler members of morphologically diverse clades have lost some developmental genes.

Comparative genomics of actinomycetes with a focus on natural product biosynthetic genes

Background

Actinomycetes are a diverse group of medically, industrially and ecologically important bacteria, studied as much for the diseases they cause as for the cures they hold. The genomes of actinomycetes revealed that these bacteria have a large number of natural product gene clusters, although many of these are difficult to tie to products in the laboratory. Large scale comparisons of these clusters are difficult to perform due to the presence of highly similar repeated domains in the most common biosynthetic machinery: polyketide synthases (PKSs) and nonribosomal peptide synthetases (NRPSs).

Results

We have used comparative genomics to provide an overview of the genomic features of a set of 102 closed genomes from this important group of bacteria with a focus on natural product biosynthetic genes. We have focused on well-represented genera and determine the occurrence of gene cluster families therein. Conservation of natural product gene clusters within Mycobacterium, Streptomyces and Frankia suggest crucial roles for natural products in the biology of each genus. The abundance of natural product classes is also found to vary greatly between genera, revealing underlying patterns that are not yet understood.

Conclusions

A large-scale analysis of natural product gene clusters presents a useful foundation for hypothesis formulation that is currently underutilized in the field. Such studies will be increasingly necessary to study the diversity and ecology of natural products as the number of genome sequences available continues to grow.

Phylogenetic analyses of phytoplasmas based on whole-genome comparison

It is no longer as difficult to determine genomic DNA sequences of uncultured bacteria as it once was, due to the development of DNA sequencing technology. It is likely that the number of whole-genome sequences of phytoplasmas will increase. In this chapter, two major strategies of whole-genome comparison studies, viz. gene content and orthologous protein sequence comparisons, are described. In general, horizontal gene transfer has greater influence on gene content-based phylogenetic analysis than orthologous protein sequence-based analysis. However, horizontal gene transfer has occurred rarely during the evolution of Mollicutes. Thus, the two phylogenetic topologies of the Mollicutes based on the two different strategies are similar.

PhyloPhlAn is a new method for improved phylogenetic and taxonomic placement of microbes

New microbial genomes are constantly being sequenced, and it is crucial to accurately determine their taxonomic identities and evolutionary relationships. Here we report PhyloPhlAn, a new method to assign microbial phylogeny and putative taxonomy using >400 proteins optimized from among 3,737 genomes. This method measures the sequence diversity of all clades, classifies genomes from deep-branching candidate divisions through closely related subspecies and improves consistency between phylogenetic and taxonomic groupings. PhyloPhlAn improved taxonomic accuracy for existing and newly sequenced genomes, detecting 157 erroneous labels, correcting 46 and placing or refining 130 new genomes. We provide examples of accurate classifications from subspecies (Sulfolobus spp.) to phyla, and of preliminary rooting of deep-branching candidate divisions, including consistent statistical support for Caldiserica (formerly candidate division OP5). PhyloPhlAn will thus be useful for both phylogenetic assessment and taxonomic quality control of newly sequenced genomes. The final phylogenies, conserved protein sequences and open-source implementation are available online.

Diversity of Streptomyces spp. in Eastern Himalayan region - computational RNomics approach to phylogeny

Isolation and characterization of actinomycetes from soil samples from altitudinal gradient of North-East India were investigated for computational RNomics based phylogeny. A total of 52 diverse isolates of Streptomyces from the soil samples were isolated on four different media and from these 6 isolates were selected on the basis of cultural characteristics, microscopic and biochemical studies. Sequencing of 16S rDNA of the selected isolates identified them to belong to six different species of Streptomyces. The molecular morphometric and physico-kinetic analysis of 16S rRNA sequences were performed to predict the diversity of the genus. The computational RNomics study revealed the significance of the structural RNA based phylogenetic analysis in a relatively diverse group of Streptomyces.

Lateral gene transfer as a support for the tree of life

Lateral gene transfer (LGT), the acquisition of genes from other species, is a major evolutionary force. However, its success as an adaptive process makes the reconstruction of the history of life an intricate puzzle: If no gene has remained unaffected during the course of life's evolution, how can one rely on molecular markers to reconstruct the relationships among species? Here, we take a completely different look at LGT and its impact for the reconstruction of the history of life. Rather than trying to remove the effect of LGT in phylogenies, and ignoring as a result most of the information of gene histories, we use an explicit phylogenetic model of gene transfer to reconcile gene histories with the tree of species. We studied 16 bacterial and archaeal phyla, representing a dataset of 12,000 gene families distributed in 336 genomes. Our results show that, in most phyla, LGT provides an abundant phylogenetic signal on the pattern of species diversification and that this signal is robust to the choice of gene families under study. We also find that LGT brings an abundant signal on the location of the root of species trees, which has been previously overlooked. Our results quantify the great variety of gene transfer rates among lineages of the tree of life and provide strong support for the "complexity hypothesis," which states that genes whose products participate to macromolecular protein complexes are relatively resistant to transfer.

Phylogenetic framework and molecular signatures for the main clades of the phylum Actinobacteria

The phylum Actinobacteria harbors many important human pathogens and also provides one of the richest sources of natural products, including numerous antibiotics and other compounds of biotechnological interest. Thus, a reliable phylogeny of this large phylum and the means to accurately identify its different constituent groups are of much interest. Detailed phylogenetic and comparative analyses of >150 actinobacterial genomes reported here form the basis for achieving these objectives. In phylogenetic trees based upon 35 conserved proteins, most of the main groups of Actinobacteria as well as a number of their superageneric clades are resolved. We also describe large numbers of molecular markers consisting of conserved signature indels in protein sequences and whole proteins that are specific for either all Actinobacteria or their different clades (viz., orders, families, genera, and subgenera) at various taxonomic levels. These signatures independently support the existence of different phylogenetic clades, and based upon them, it is now possible to delimit the phylum Actinobacteria (excluding Coriobacteriia) and most of its major groups in clear molecular terms. The species distribution patterns of these markers also provide important information regarding the interrelationships among different main orders of Actinobacteria. The identified molecular markers, in addition to enabling the development of a stable and reliable phylogenetic framework for this phylum, also provide novel and powerful means for the identification of different groups of Actinobacteria in diverse environments. Genetic and biochemical studies on these Actinobacteria-specific markers should lead to the discovery of novel biochemical and/or other properties that are unique to different groups of Actinobacteria.

Multifractal analysis of HIV-1 genomes

Pathogens like HIV-1, which evolve into many closely related variants displaying differential infectivity and evolutionary dynamics in a short time scale, require fast and accurate classification. Conventional whole genome sequence alignment-based methods are computationally expensive and involve complex analysis. Alignment-free methodologies are increasingly being used to effectively differentiate genomic variations between viral species. Multifractal analysis, which explores the self-similar nature of genomes, is an alignment-free methodology that has been applied to study such variations. However, whether multifractal analysis can quantify variations between closely related genomes, such as the HIV-1 subtypes, is an open question. Here we address the above by implementing the multifractal analysis on four retroviral genomes (HIV-1, HIV-2, SIVcpz, and HTLV-1), and demonstrate that individual multifractal properties can differentiate between different retrovirus types easily. However, the individual multifractal measures do not resolve within-group variations for different known subtypes of HIV-1 M group. We show here that these known subtypes can instead be classified correctly using a combination of the crucial multifractal measures. This method is simple and computationally fast in comparison to the conventional alignment-based methods for whole genome phylogenetic analysis.

Prioritizing orphan proteins for further study using phylogenomics and gene expression profiles in Streptomyces coelicolor

Background

Streptomyces coelicolor, a model organism of antibiotic producing bacteria, has one of the largest genomes of the bacterial kingdom, including 7825 predicted protein coding genes. A large number of these genes, nearly 34%, are functionally orphan (hypothetical proteins with unknown function). However, in gene expression time course data, many of these functionally orphan genes show interesting expression patterns.

Results

In this paper, we analyzed all functionally orphan genes of Streptomyces coelicolor and identified a list of "high priority" orphans by combining gene expression analysis and additional phylogenetic information (i.e. the level of evolutionary conservation of each protein).

Conclusions

The prioritized orphan genes are promising candidates to be examined experimentally in the lab for further characterization of their function.

Whole-genome comparison clarifies close phylogenetic relationships between the phyla Dictyoglomi and Thermotogae

The anaerobic thermophilic bacterial genus Dictyoglomus is characterized by the ability to produce useful enzymes such as amylase, mannanase, and xylanase. Despite the significance, the phylogenetic position of Dictyoglomus has not yet been clarified, since it exhibits ambiguous phylogenetic positions in a single gene sequence comparison-based analysis. The number of substitutions at the diverging point of Dictyoglomus is insufficient to show the relationships in a single gene comparison-based analysis. Hence, we studied its evolutionary trait based on whole-genome comparison. Both gene content and orthologous protein sequence comparisons indicated that Dictyoglomus is most closely related to the phylum Thermotogae and it forms a monophyletic group with Coprothermobacter proteolyticus (a constituent of the phylum Firmicutes) and Thermotogae. Our findings indicate that C. proteolyticus does not belong to the phylum Firmicutes and that the phylum Dictyoglomi is not closely related to either the phylum Firmicutes or Synergistetes but to the phylum Thermotogae.

Comparative genome-scale metabolic modeling of actinomycetes: the topology of essential core metabolism

Actinomycetes are highly important bacteria. On one hand, some of them cause severe human and plant diseases, on the other hand, many species are known for their ability to produce antibiotics. Here we report the results of a comparative analysis of genome-scale metabolic models of 37 species of actinomycetes. Based on in silico knockouts we generated topological and genomic maps for each organism. Combining the collection of genome-wide models, we constructed a global enzyme association network to identify both a conserved "core network" and an "essential core network" of the entire group. As has been reported for low-degree metabolites in several organisms, low-degree enzymes (in linear pathways) turn out to be generally more essential than high-degree enzymes (in metabolic hubs).

Presence and diversity of Streptomyces in Dendroctonus and sympatric bark beetle galleries across North America

Recent studies have revealed several examples of intimate associations between insects and Actinobacteria, including the Southern Pine Beetle Dendroctonus frontalis and the Spruce Beetle Dendroctonus rufipennis. Here, we surveyed Streptomyces Actinobacteria co-occurring with 10 species of Dendroctonus bark beetles across the United States, using both phylogenetic and community ecology approaches. From these 10 species, and 19 other scolytine beetles that occur in the same trees, we obtained 154 Streptomyces-like isolates and generated 16S sequences from 134 of those. Confirmed 16S sequences of Streptomyces were binned into 36 distinct strains using a threshold of 0.2% sequence divergence. The 16S rDNA phylogeny of all isolates does not correlate with the distribution of strains among beetle species, localities, or parts of the beetles or their galleries. However, we identified three Streptomyces strains occurring repeatedly on Dendroctonus beetles and in their galleries. Identity of these isolates was corroborated using a house-keeping gene sequence (efTu). These strains are not confined to a certain species of beetle, locality, or part of the beetle or their galleries. However, their role as residents in the woodboring insect niche is supported by the repeated association of their 16S and efTu from across the continent, and also having been reported in studies of other subcortical insects.