Rapid Annotation Using Subsystem Technology Research Articles

Genome Organization of Escherichia Phage YD-2008.s: A New Entry to Siphoviridae Family.

Malaysia is one of the countries that are loaded with mega biodiversity which includes microbial communities. Phages constitute the major component in the microbial communities and yet the numbers of discovered phages are just a minute fraction of its population in the biosphere. Taking into account of a huge numbers of waiting to be discovered phages, a new bacteriophage designated as Escherichia phage YD-2008.s was successfully isolated using Escherichia coli ATCC 11303 as the host. Phage YD-2008.s poses icosahedral head measured at 57nm in diameter with a long non-contractile flexible tail measured at 107nm; proving the phage as one of the members of Siphoviridae family under the order of Caudovirales. Genomic sequence analyses revealed phage YD-2008.s genome as linear dsDNA of 44,613 base pairs with 54.6% G+C content. Sixty-two open reading frames (ORFs) were identified on phage YD-2008.s full genome, using bioinformatics annotation software; Rapid Annotation using Subsystem Technology (RAST). Among the ORFs, twenty-eight of them code for functional proteins. Thirty two are classified as hypothetical proteins and there are two unidentified proteins. Even though majority of the coded putative proteins have high amino acids similarities to phages from the genus Hk578likevirus of the Siphoviridae family, yet phage YD-2008.s stands with its’ own distinctiveness. Therefore, this is another new finding to Siphoviridae family as well as to the growing list of viruses in International Committee on Taxonomy of Viruses (ICTV) database.

Phylogeny and putative virulence gene analysis of <i>Bartonella bovis</i>

Bartonella bovis is a small Gram-negative bacterium recognized as an etiological agent for bacteremia and endocarditis in cattle. As few reports are available on the taxonomic position of B. bovis and its mechanism of virulence, this study aims to resolve the phylogeny of B. bovis and investigate putative virulence genes based on whole genome sequence analysis. Genome-wide comparisons based on single nucleotide polymorphisms (SNP) and orthologous genes were performed in this study for phylogenetic inference of 27 Bartonella species. Rapid Annotation using Subsystem Technology (RAST) analysis was used for annotation of putative virulence genes. The phylogenetic tree generated from the genome-wide comparison of orthologous genes exhibited a topology almost similar to that of the tree generated from SNP-based comparison, indicating a high concordance in the nucleotide and amino acid sequences of Bartonella spp. The analyses show consistent grouping of B. bovis in a cluster related to ruminant-associated species, including Bartonella australis, Bartonella melophagi and Bartonella schoenbuchensis. RAST analysis revealed genes encoding flagellar components, in corroboration with the observation of flagella-like structure of BbUM strain under negative straining. Genes associated with virulence, disease and defence, prophages, membrane transport, iron acquisition, motility and chemotaxis are annotated in B. bovis genome. The flagellin (flaA) gene of B. bovis is closely related to Bartonella bacilliformis and Bartonella clarridgeiae but distinct from other Gram-negative bacteria. The absence of type IV secretion systems, the bona fide pathogenicity factors of bartonellae, in B. bovis suggests that it may have a different mechanism of pathogenicity.

Epidemiological studies based on multi-locus sequence typing genotype of methicillin susceptible Staphylococcus aureus isolated from camel's milk.

One hundred milk samples were collected from camel’s milk for the isolation of Staphylococcus aureus. Thirty-one isolates were S. aureus, 45 were other forms of staphylococci and 24 represented other bacteria. Five isolates from S. aureus were methicillin resistant S. aureus (MRSA) and 26 samples were methicillin susceptible S. aureus (MSSA). The whole genome sequence of S. aureus was annotated and visualised by rapid annotation using subsystem technology (RAST) which is a fully-automated service for annotating complete or nearly complete bacterial genomes. Four isolates from MSSA strains were subjected to multi-locus sequence typing (MLST). Three multi-locus sequences types or sequence types (MLST/ST) were found, namely ST15, ST1153 and ST130. The phylogenetic analysis of the concatenated sequences of the seven genes forming the MLST profile of S. aureus classification revealed a high degree of similarity and close relationship between the ST15 and ST1153 while the third ST (ST130) was located in a different cluster.

Complete genetic analysis of a Salmonella enterica serovar Indiana isolate accompanying four plasmids carrying mcr-1, ESBL and other resistance genes in China

One mcr-1-carrying Salmonella enterica serovar Indiana strain D90, was identified from 1320 Salmonella enterica isolates from poultry slaughterhouse in 2012 in China. The objective of this study was to verify the transferability of the mcr-1 gene and also completely characterize the sequence of the strain at the whole-genome level. Broth matting assays were carried out to detect the transferability and whole-genome sequencing (WGS) of S. enterica serovar Indiana D90 was performed using the PacBio RS II system. Open reading frames were assigned using Rapid Annotation using Subsystem Technology (RAST) and analysed by BLASTn and BLASTp. Salmonella Pathogenisity Islands (SPIs) were annotated by SPIFinder platform. The complete genome sequence of S. enterica serovar Indiana D90 contained a circular 4,779,514-bp chromosome and four plasmids. Genome analysis and sequencing revealed that 24 multi-drug resistance (MDR) genes were located on plasmids. The largest plasmid pD90-1, was found to be of an IncHI2/HI2A/Q1/N type that encoded a blaCTX-M-65 gene along with 20 additional antimicrobial resistance genes. A 60.5-kbp IncI2 plasmid pD90-2 contained a nikA-nikB-mcr-1 genetic structure, that can be successfully transferred to E. coli and S. enterica serovar Typhimurium at low transfer rates. Interestingly, comparative sequence analysis revealed the plasmids pD90-1 and pD90-2 showed considerable nucleotide similarity to pHNSHP45-2 and pHNSHP45, respectively. Moreover, the genome and the plasmid pD90-2 also showed high similarity to one carbapenem resistant S. enterica serovar Indiana strain, C629 and its plasmid pC629, respectively. This is the first report of the complete nucleotide sequence of one mcr-1-carrying MDR S. enterica serovar Indiana strain.

Next generation sequencing reveals the antibiotic resistant variants in the genome of Pseudomonas aeruginosa.

Rapid progress in next generation sequencing and allied computational tools have aided in identification of single nucleotide variants in genomes of several organisms. In the present study, we have investigated single nucleotide polymorphism (SNP) in ten multi-antibiotic resistant Pseudomonas aeruginosa clinical isolates. All the draft genomes were submitted to Rapid Annotations using Subsystems Technology (RAST) web server and the predicted protein sequences were used for comparison. Non-synonymous single nucleotide polymorphism (nsSNP) found in the clinical isolates compared to the reference genome (PAO1), and the comparison of nsSNPs between antibiotic resistant and susceptible clinical isolates revealed insights into the genome variation. These nsSNPs identified in the multi-drug resistant clinical isolates were found to be altering a single amino acid in several antibiotic resistant genes. We found mutations in genes encoding efflux pump systems, cell wall, DNA replication and genes involved in repair mechanism. In addition, nucleotide deletions in the genome and mutations leading to generation of stop codons were also observed in the antibiotic resistant clinical isolates. Next generation sequencing is a powerful tool to compare the whole genomes and analyse the single base pair variations found within the antibiotic resistant genes. We identified specific mutations within antibiotic resistant genes compared to the susceptible strain of the same bacterial species and these findings may provide insights to understand the role of single nucleotide variants in antibiotic resistance.

Draft Genome Sequence of the Deep-Sea Bacterium Moritella sp. JT01 and Identification of Biotechnologically Relevant Genes.

Deep-sea bacteria can produce various biotechnologically relevant enzymes due to their adaptations to high pressures and low temperatures. To identify such enzymes, we have sequenced the genome of the polycaprolactone-degrading bacterium Moritella sp. JT01, isolated from sediment samples from Japan Trench (6957m depth), using a Illumina HiSeq2000 sequencer (12.1 million paired-end reads) and CLC Genomics Workbench (version 6.5.1) for the assembly, resulting in a 4.83-Mb genome (42 scaffolds). The genome was annotated using Rapid Annotation using Subsystem Technology (RAST), Protein Homology/analogY Recognition Engine V 2.0 (PHYRE2), and BLAST2Go, revealing 4439 protein coding sequences and 101 RNAs. Gene products with industrial relevance, such as lipases (three) and esterases (four), were identified and are related to bacterium's ability to degrade polycaprolactone. The annotation revealed proteins related to deep-sea survival, such as cold-shock proteins (six) and desaturases (three). The presence of secondary metabolite biosynthetic gene clusters suggests that this bacterium could produce nonribosomal peptides, polyunsaturated fatty acids, and bacteriocins. To demonstrate the potential of this genome, a lipase was cloned an introduced into Escherichia coli. The lipase was purified and characterized, showing activity over a wide temperature range (over 50% at 20-60°C) and pH range (over 80% at pH6.3 to 9). This enzyme has tolerance to the surfactant action of sodium dodecyl sulfate and shows 30% increased activity when subjected to a working pressure of 200MPa. The genomic characterization of Moritella sp. JT01 reveals traits associated with survival in the deep-sea and their potential uses in biotechnology, as exemplified by the characterized lipase.

Genome Insight and Comparative Pathogenomic Analysis of Nesterenkonia jeotgali Strain CD08_7 Isolated from Duodenal Mucosa of Celiac Disease Patient.

Species of the genus Nesterenkonia have been isolated from different ecological niches, especially from saline habitats and reported as weak human pathogens causing asymptomatic bacteraemia. Here, for the first time we are reporting the genome sequence and pathogenomic analysis of a strain designated as CD08_7 isolated from the duodenal mucosa of a celiac disease patient, identified as Nesterenkonia jeotgali. To date, only five strains of the genus Nesterenkonia (N. massiliensis strain NP1T, Nesterenkonia sp. strain JCM 19054, Nesterenkonia sp. strain F and Nesterenkonia sp. strain AN1) have been whole genome sequenced and annotated. In the present study we have mapped and compared the virulence profile of N. jeotgali strain CD08_7 along with other reference genomes which showed some characteristic features that could contribute to pathogenicity. The RAST (Rapid Annotation using Subsystem Technology) based genome mining revealed more genes responsible for pathogenicity in strain CD08_7 when compared with the other four sequenced strains. The studied categories were resistance to antibiotic and toxic compounds, invasion and intracellular resistance, membrane transport, stress response, osmotic stress, oxidative stress, phages and prophages and iron acquisition. A total of 1431 protein-encoding genes were identified in the genome of strain CD08_7 among which 163 were predicted to contribute for pathogenicity. Out of 163 genes only 59 were common to other genome, which shows the higher levels of genetic richness in strain CD08_7 that may contribute to its functional versatility. This study provides a comprehensive analysis on genome of N. jeotgali strain CD08_7 and possibly indicates its importance as a clinical pathogen.

Agarolytic bacterium Persicobacter sp. CCB-QB2 exhibited a diauxic growth involving galactose utilization pathway.

The agarolytic bacterium Persicobacter sp. CCB‐QB2 was isolated from seaweed (genus Ulva) collected from a coastal area of Malaysia. Here, we report a high‐quality draft genome sequence for QB2. The Rapid Annotation using Subsystem Technology (RAST) annotation server identified four β‐agarases (PdAgaA, PdAgaB, PdAgaC, and PdAgaD) as well as galK, galE, and phosphoglucomutase, which are related to the Leloir pathway. Interestingly, QB2 exhibited a diauxic growth in the presence of two kinds of nutrients, such as tryptone and agar. In cells grown with agar, the profiles of agarase activity and growth rate were very similar. galK, galE, and phosphoglucomutase genes were highly expressed in the second growth phase of diauxic growth, indicating that QB2 cells use galactose hydrolyzed from agar by its agarases and exhibit nutrient prioritization. This is the first report describing diauxic growth for agarolytic bacteria. QB2 is a potential novel model organism for studying diauxic growth in environmental bacteria.

Elucidating genomic gaps using phenotypic profiles

Advances in genomic sequencing provide the ability to model the metabolism of organisms from their genome annotation. The bioinformatics tools developed to deduce gene function through homology-based methods are dependent on public databases; thus, novel discoveries are not readily extrapolated from current analysis tools with a homology dependence. Multi-phenotype Assay Plates (MAPs) provide a high-throughput method to profile bacterial phenotypes by growing bacteria in various growth conditions, simultaneously. More robust and accurate computational models can be constructed by coupling MAPs with current genomic annotation methods.PMAnalyzeris an online tool that analyzes bacterial growth curves from the MAP system which are then used to optimize metabolic models duringin silicogrowth simulations. UsingCitrobacter sedlakiias a prototype, the Rapid Annotation using Subsystem Technology (RAST) tool produced a model consisting of 1,367 enzymatic reactions. After the optimization, 44 reactions were added to, or modified within, the model. The model correctly predicted the outcome on 93% of growth experiments.

A molecular phylogenetic framework for Bacillus subtilis using genome sequences and its application to Bacillus subtilis subspecies stecoris strain D7XPN1, an isolate from a commercial food-waste degrading bioreactor.

A thermophilic, heterotrophic and facultatively anaerobic bacterium designated strain D7XPN1 was isolated from Baku BakuKing™, a commercial food-waste degrading bioreactor (composter). The strain grew optimally at 45 °C (growth range between 24 and 50 °C) and pH 7 (growth pH range between pH 5 and 9) in Luria Broth supplemented with 0.3 % glucose. Strain D7XPN1 tolerated up to 7 % NaCl and showed amylolytic and xylanolytic activities. 16S rRNA gene analysis placed strain D7XPN1 in the cluster represented by Bacillus subtilis and the genome analysis of the 4.1 Mb genome sequence determined using RAST (Rapid Annotation using Subsystem Technology) indicated a total of 5116 genomic features were present of which 2320 features could be grouped into several subsystem categories. Of these, 615 features were related to carbohydrate metabolism which included a range of enzymes with potential in the biodegradation of food wastes, a property consistent with the ecological habitat of the isolate. ANIb (Average Nucleotide Identity based on BLAST) analysis with 49 Bacillus subtilis genomes indicated that it was distantly related to the three currently taxonomically validated B. subtilis subspecies namely B. subtilis subsp. subtilis (95.6 %), B. subtilis subsp. spizizenii (93 %) and B. subtilis subsp. inaquosorum (92 %) and based on our current knowledge warranted that it be included as a separate cluster together with strain JS which it was closely related (98.69 %). The close relationship of strains D7XPN1 and JS is also supported from our results from electronic DNA–DNA Hybridization (e-DDH) studies. Furthermore, our additional in-depth phylogenomic analyses using three different datasets unequivocally supported the creation of a fourth B. subtilis subspecies to include strains D7XPN1 and JS for which we propose strain D7XPN1T (=KCTC 33554T, JCM 30051T) as the type strain, and designate it as B. subtilis subsp. stecoris.

PATtyFams: Protein Families for the Microbial Genomes in the PATRIC Database.

The ability to build accurate protein families is a fundamental operation in bioinformatics that influences comparative analyses, genome annotation, and metabolic modeling. For several years we have been maintaining protein families for all microbial genomes in the PATRIC database (Pathosystems Resource Integration Center, patricbrc.org) in order to drive many of the comparative analysis tools that are available through the PATRIC website. However, due to the burgeoning number of genomes, traditional approaches for generating protein families are becoming prohibitive. In this report, we describe a new approach for generating protein families, which we call PATtyFams. This method uses the k-mer-based function assignments available through RAST (Rapid Annotation using Subsystem Technology) to rapidly guide family formation, and then differentiates the function-based groups into families using a Markov Cluster algorithm (MCL). This new approach for generating protein families is rapid, scalable and has properties that are consistent with alignment-based methods.

RASTtk: a modular and extensible implementation of the RAST algorithm for building custom annotation pipelines and annotating batches of genomes.

The RAST (Rapid Annotation using Subsystem Technology) annotation engine was built in 2008 to annotate bacterial and archaeal genomes. It works by offering a standard software pipeline for identifying genomic features (i.e., protein-encoding genes and RNA) and annotating their functions. Recently, in order to make RAST a more useful research tool and to keep pace with advancements in bioinformatics, it has become desirable to build a version of RAST that is both customizable and extensible. In this paper, we describe the RAST tool kit (RASTtk), a modular version of RAST that enables researchers to build custom annotation pipelines. RASTtk offers a choice of software for identifying and annotating genomic features as well as the ability to add custom features to an annotation job. RASTtk also accommodates the batch submission of genomes and the ability to customize annotation protocols for batch submissions. This is the first major software restructuring of RAST since its inception.

Elucidating genomic gaps using phenotypic profiles

Advances in genomic sequencing provide the ability to model the metabolism of organisms from their genome annotation. The bioinformatics tools developed to deduce gene function through homology-based methods are dependent on public databases; thus, novel discoveries are not readily extrapolated from current analysis tools with a homology dependence. Multi-phenotype Assay Plates (MAPs) provide a high-throughput method to profile bacterial phenotypes by growing bacteria in various growth conditions, simultaneously. More robust and accurate computational models can be constructed by coupling MAPs with current genomic annotation methods.PMAnalyzeris an online tool that analyzes bacterial growth curves from the MAP system which are then used to optimize metabolic models duringin silicogrowth simulations. UsingCitrobacter sedlakiias a prototype, the Rapid Annotation using Subsystem Technology (RAST) tool produced a model consisting of 1,367 enzymatic reactions. After the optimization, 44 reactions were added to, or modified within, the model. The model correctly predicted the outcome on 93% of growth experiments.

Rapid construction of metabolic models for a family of Cyanobacteria using a multiple source annotation workflow

BackgroundCyanobacteria are photoautotrophic prokaryotes that exhibit robust growth under diverse environmental conditions with minimal nutritional requirements. They can use solar energy to convert CO2 and other reduced carbon sources into biofuels and chemical products. The genus Cyanothece includes unicellular nitrogen-fixing cyanobacteria that have been shown to offer high levels of hydrogen production and nitrogen fixation. The reconstruction of quality genome-scale metabolic models for organisms with limited annotation resources remains a challenging task.ResultsHere we reconstruct and subsequently analyze and compare the metabolism of five Cyanothece strains, namely Cyanothece sp. PCC 7424, 7425, 7822, 8801 and 8802, as the genome-scale metabolic reconstructions iCyc792, iCyn731, iCyj826, iCyp752, and iCyh755 respectively. We compare these phylogenetically related Cyanothece strains to assess their bio-production potential. A systematic workflow is introduced for integrating and prioritizing annotation information from the Universal Protein Resource (Uniprot), NCBI Protein Clusters, and the Rapid Annotations using Subsystems Technology (RAST) method. The genome-scale metabolic models include fully traced photosynthesis reactions and respiratory chains, as well as balanced reactions and GPR associations. Metabolic differences between the organisms are highlighted such as the non-fermentative pathway for alcohol production found in only Cyanothece 7424, 8801, and 8802.ConclusionsOur development workflow provides a path for constructing models using information from curated models of related organisms and reviewed gene annotations. This effort lays the foundation for the expedient construction of curated metabolic models for organisms that, while not being the target of comprehensive research, have a sequenced genome and are related to an organism with a curated metabolic model. Organism-specific models, such as the five presented in this paper, can be used to identify optimal genetic manipulations for targeted metabolite overproduction as well as to investigate the biology of diverse organisms.

Taxonomic and Functional Microbial Signatures of the Endemic Marine Sponge Arenosclera brasiliensis

The endemic marine sponge Arenosclera brasiliensis (Porifera, Demospongiae, Haplosclerida) is a known source of secondary metabolites such as arenosclerins A-C. In the present study, we established the composition of the A. brasiliensis microbiome and the metabolic pathways associated with this community. We used 454 shotgun pyrosequencing to generate approximately 640,000 high-quality sponge-derived sequences (∼150 Mb). Clustering analysis including sponge, seawater and twenty-three other metagenomes derived from marine animal microbiomes shows that A. brasiliensis contains a specific microbiome. Fourteen bacterial phyla (including Proteobacteria, Cyanobacteria, Actinobacteria, Bacteroidetes, Firmicutes and Cloroflexi) were consistently found in the A. brasiliensis metagenomes. The A. brasiliensis microbiome is enriched for Betaproteobacteria (e.g., Burkholderia) and Gammaproteobacteria (e.g., Pseudomonas and Alteromonas) compared with the surrounding planktonic microbial communities. Functional analysis based on Rapid Annotation using Subsystem Technology (RAST) indicated that the A. brasiliensis microbiome is enriched for sequences associated with membrane transport and one-carbon metabolism. In addition, there was an overrepresentation of sequences associated with aerobic and anaerobic metabolism as well as the synthesis and degradation of secondary metabolites. This study represents the first analysis of sponge-associated microbial communities via shotgun pyrosequencing, a strategy commonly applied in similar analyses in other marine invertebrate hosts, such as corals and algae. We demonstrate that A. brasiliensis has a unique microbiome that is distinct from that of the surrounding planktonic microbes and from other marine organisms, indicating a species-specific microbiome.

Genome Sequences of the Ethanol-Tolerant Lactobacillus vini Strains LMG 23202 T and JP7.8.9

We report on the genome sequences of Lactobacillus vini type strain LMG 23202(T) (DSM 20605) (isolated from fermenting grape musts in Spain) and the industrial strain L. vini JP7.8.9 (isolated from a bioethanol plant in northeast Brazil). All contigs were assembled using gsAssembler, and genes were predicted and annotated using Rapid Annotation using Subsystem Technology (RAST). The identified genome sequence of LMG 23202(T) had 2.201.333 bp, 37.6% G+C, and 1,833 genes, whereas the identified genome sequence of JP7.8.9 had 2.301.037 bp, 37.8% G+C, and 1,739 genes. The gene repertoire of the species L. vini offers promising opportunities for biotechnological applications.

Genome Sequence of the Bacterioplanktonic, Mixotrophic Vibrio campbellii Strain PEL22A, Isolated in the Abrolhos Bank

Vibrio campbellii PEL22A was isolated from open ocean water in the Abrolhos Bank. The genome of PEL22A consists of 6,788,038 bp (the GC content is 45%). The number of coding sequences (CDS) is 6,359, as determined according to the Rapid Annotation using Subsystem Technology (RAST) server. The number of ribosomal genes is 80, of which 68 are tRNAs and 12 are rRNAs. V. campbellii PEL22A contains genes related to virulence and fitness, including a complete proteorhodopsin cluster, complete type II and III secretion systems, incomplete type I, IV, and VI secretion systems, a hemolysin, and CTXΦ.

Analysis of the antimicrobial susceptibility of the ionizing radiation-resistant bacterium Deinococcus radiodurans: implications for bioremediation of radioactive waste

Deinococcus radiodurans R1 is one of the most resistant bacteria to chronic ionizing radiation ever identified. Therefore, Deinococcus radiodurans shows promise in the bioremediation of mixed radioactive wastes. Nuclear wastes have been demonstrated to host bacteria containing genes for antimicrobial production (Bagwell et al. PloS ONE 3:e3878, 2008; Bai et al. Chem Biol 13:387–397, 2006; Phillips et al. Int J Syst Evol Microbiol 52:933–938, 2002), which may represent a constraint on the bioremediation potential of Deinococcus radiodurans. In this context, the aim of this work is to investigate the antimicrobial susceptibility of D. radiodurans through in silico and in vitro approaches. In silico, we initiated our analyses by retrieving genes that are predicted to confer to D. radiodurans resistance to antimicrobials based on the Antibiotic Resistance Genes Database (ARDB), the Rapid Annotations using Subsystems Technology (RAST), and the automatic hierarchical classification of proteins (ProtoNet) servers. Among 19 retrieved sequences of D. radiodurans, six genes were functionally re-annotated as antimicrobial resistance genes, and only two genes were found to possess the radiation/desiccation response motif (RDRM). In vitro, we found that D. radiodurans is sensitive to 40 tested antimicrobials, and it was resistant to one macrolide (spiramycin) and to three quinolones (ofloxacin, ciprofloxacin, and nalidixic acid). The resistance of D. radiodurans to these quinolones was discussed based on different possibilities. Moreover, in silico analyses indicated that the investigated 19 genes of D. radiodurans from ARDB, RAST, and ProtoNet encode translatable messenger ribonucleic acids (mRNAs) related to resistance to antimicrobials. These data together with the antibiograms of D. radiodurans represent a suggestive evidence that the majority of these genes encode non-functional proteins or belong to inefficient biochemical pathways. The susceptibility of D. radiodurans to these antimicrobials limits its potential for bioremediation of nuclear waste and further genetic engineering is needed depending on the site to be depolluted.