Batch Blast Extractor: an automated blastx parser application

MotivationBLAST programs are very efficient in finding similarities for sequences. However for large datasets such as ESTs, manual extraction of the information from the batch BLAST output is needed. This can be time consuming, insufficient, and inaccurate. Therefore implementation of a parser application would be extremely useful in extracting information from BLAST outputs. ResultsWe have developed a java application, Batch Blast Extractor, with a user friendly graphical interface to extract information from BLAST output. The application generates a tab delimited text file that can be easily imported into any statistical package such as Excel or SPSS for further analysis. For each BLAST hit, the program obtains and saves the essential features from the BLAST output file that would allow further analysis. The program was written in Java and therefore is OS independent. It works on both Windows and Linux OS with java 1.4 and higher. It is freely available from: http://mcbc.usm.edu/BatchBlastExtractor

RiboaptDB: A Comprehensive Database of Ribozymes and Aptamers

BACKGROUND: Catalytic RNA molecules are called ribozymes. The aptamers are DNA or RNA molecules that have been selected from vast populations of random sequences, through a combinatorial approach known as SELEX. The selected oligo-nucleotide sequences (~200 bp in length) have the ability to recognize a broad range of specific ligands by forming binding pockets. These novel aptamer sequences can bind to nucleic acids, proteins or small organic and inorganic chemical compounds and have many potential uses in medicine and technology. RESULTS: The comprehensive sequence information on aptamers and ribozymes that have been generated by in vitro selection methods are included in this RiboaptDB database. Such types of unnatural data generated by in vitro methods are not available in the public 'natural' sequence databases such as GenBank and EMBL. The amount of sequence data generated by in vitro selection experiments has been accumulating exponentially. There are 370 artificial ribozyme sequences and 3842 aptamer sequences in the total 4212 sequences from 423 citations in this RiboaptDB. We included general search feature, and individual feature wise search, user submission form for new data through online and also local BLAST search. CONCLUSION: This database, besides serving as a storehouse of sequences that may have diagnostic or therapeutic utility in medicine, provides valuable information for computational and theoretical biologists. The RiboaptDB is extremely useful for garnering information about in vitro selection experiments as a whole and for better understanding the distribution of functional nucleic acids in sequence space. The database is updated regularly and is publicly available at

Improving the Modeling of Sea Surface Currents in the Persian Gulf and the Oman Sea Using Data Assimilation of Satellite Altimetry and Hydrographic Observations

Sea surface currents are often modeled using numerical models without adequately addressing the issue of model calibration at the regional scale. The aim of this study is to calibrate the MIKE 21 numerical ocean model for the Persian Gulf and the Oman Sea to improve the sea surface currents obtained from the model. The calibration was performed through data assimilation of the model with altimetry and hydrographic observations using variational data assimilation, where the weights of the objective functions were defined based on the type of observations and optimized using metaheuristic optimization methods. According to the results, the calibration of the model generally led the model results closer to the observations. This was reflected in an improvement of about 0.09 m/s in the obtained sea surface currents. It also allowed for more accurate evaluations of model parameters, such as Smagorinsky and Manning coefficients. Moreover, the root mean square error values between the satellite altimetry observations at control stations and the assimilated model varied between 0.058 and 0.085 m. We further showed that the kinetic energy produced by sea surface currents could be used for generating electricity in the Oman Sea and near Jask harbor

Generation, Analysis and Functional Annotation of Expressed Sequence Tags From the Sheepshead Minnow (\u3ci\u3eCyprinodon variegatus\u3c/i\u3e)

Background: Sheepshead minnow (Cyprinodon variegatus) are small fish capable of withstanding exposure to very low levels of dissolved oxygen, as well as extreme temperatures and salinities. It is an important model in understanding the impacts and biological response to hypoxia and co-occurring compounding stressors such as polycyclic aromatic hydrocarbons, endocrine disrupting chemicals, metals and herbicides. Here, we initiated a project to sequence and analyze over 10,000 ESTs generated from the Sheepshead minnow (Cyprinodon variegatus) as a resource for investigating stressor responses. Results: We sequenced 10,858 EST clones using a normalized cDNA library made from larval, embryonic and adult suppression subtractive hybridization-PCR (SSH) libraries. Post- sequencing processing led to 8,099 high quality sequences. Clustering analysis of these ESTs indentified 4,223 unique sequences containing 1,053 contigs and 3,170 singletons. BLASTX searches produced 1,394 significant (E-value \u3c 10(-5)) hits and further Gene Ontology (GO) analysis annotated 388 of these genes. All the EST sequences were deposited by Expressed Sequence Tags database (dbEST) in GenBank (GenBank: GE329585 to GE337683). Gene discovery and annotations are presented and discussed. This set of ESTs represents a significant proportion of the Sheepshead minnow (Cyprinodon variegatus) transcriptome, and provides a material basis for the development of microarrays useful for further gene expression studies in association with stressors such as hypoxia, cadmium, chromium and pyrene

Transcriptomic analysis of RDX and TNT interactive sublethal effects in the earthworm Eisenia fetida

<p>Abstract</p> <p>Background</p> <p>Explosive compounds such as TNT and RDX are recalcitrant contaminants often found co-existing in the environment. In order to understand the joint effects of TNT and RDX on earthworms, an important ecological and bioindicator species at the molecular level, we sampled worms (<it>Eisenia fetida</it>) exposed singly or jointly to TNT (50 mg/kg soil) and RDX (30 mg/kg soil) for 28 days and profiled gene expression in an interwoven loop designed microarray experiment using a 4k-cDNA array. Lethality, growth and reproductive endpoints were measured.</p> <p>Results</p> <p>Sublethal doses of TNT and RDX had no significant effects on the survival and growth of earthworms, but significantly reduced cocoon and juvenile counts. The mixture exhibited more pronounced reproductive toxicity than each single compound, suggesting an additive interaction between the two compounds. In comparison with the controls, we identified 321 differentially expressed transcripts in TNT treated worms, 32 in RDX treated worms, and only 6 in mixture treated worms. Of the 329 unique differentially expressed transcripts, 294 were affected only by TNT, 24 were common to both TNT and RDX treatments, and 3 were common to all treatments. The reduced effects on gene expression in the mixture exposure suggest that RDX might interact in an antagonistic manner with TNT at the gene expression level. The disagreement between gene expression and reproduction results may be attributed to sampling time, absence of known reproduction-related genes, and lack of functional information for many differentially expressed transcripts. A gene potentially related to reproduction (echinonectin) was significantly depressed in TNT or RDX exposed worms and may be linked to reduced fecundity.</p> <p>Conclusions</p> <p>Sublethal doses of TNT and RDX affected many biological pathways from innate immune response to oogenesis, leading to reduced reproduction without affecting survival and growth. A complex interaction between mixtures of RDX and TNT was observed at the gene expression level that requires further study of the dynamics of gene expression and reproductive activities in <it>E. fetida</it>. These efforts will be essential to gain an understanding of the additive reproductive toxicity between RDX and TNT.</p

ILOOP – a web application for two-channel microarray interwoven loop design

Microarray technology is widely applied to address complex scientific questions. However, there remain fundamental issues on how to design experiments to ensure that the resulting data enables robust statistical analysis. Interwoven loop design has several advantages over other designs. However it suffers in the complexity of design. We have implemented an online web application which allows users to find optimal loop designs for two-color microarray experiments. Given a number of conditions (such as treatments or time points) and replicates, the application will find the best possible design of the experiment and output experimental parameters. It is freely available from

Cloning, analysis and functional annotation of expressed sequence tags from the Earthworm Eisenia fetida

<p>Abstract</p> <p>Background</p> <p><it>Eisenia fetida</it>, commonly known as red wiggler or compost worm, belongs to the Lumbricidae family of the Annelida phylum. Little is known about its genome sequence although it has been extensively used as a test organism in terrestrial ecotoxicology. In order to understand its gene expression response to environmental contaminants, we cloned 4032 cDNAs or expressed sequence tags (ESTs) from two <it>E. fetida </it>libraries enriched with genes responsive to ten ordnance related compounds using suppressive subtractive hybridization-PCR.</p> <p>Results</p> <p>A total of 3144 good quality ESTs (GenBank dbEST accession number <ext-link ext-link-type="gen" ext-link-id="EH669363">EH669363</ext-link>–<ext-link ext-link-type="gen" ext-link-id="EH672369">EH672369</ext-link> and <ext-link ext-link-type="gen" ext-link-id="EL515444">EL515444</ext-link>–<ext-link ext-link-type="gen" ext-link-id="EL515580">EL515580</ext-link>) were obtained from the raw clone sequences after cleaning. Clustering analysis yielded 2231 unique sequences including 448 contigs (from 1361 ESTs) and 1783 singletons. Comparative genomic analysis showed that 743 or 33% of the unique sequences shared high similarity with existing genes in the GenBank nr database. Provisional function annotation assigned 830 Gene Ontology terms to 517 unique sequences based on their homology with the annotated genomes of four model organisms <it>Drosophila melanogaster</it>, <it>Mus musculus</it>, <it>Saccharomyces cerevisiae</it>, and <it>Caenorhabditis elegans</it>. Seven percent of the unique sequences were further mapped to 99 Kyoto Encyclopedia of Genes and Genomes pathways based on their matching Enzyme Commission numbers. All the information is stored and retrievable at a highly performed, web-based and user-friendly relational database called EST model database or ESTMD version 2.</p> <p>Conclusion</p> <p>The ESTMD containing the sequence and annotation information of 4032 <it>E. fetida </it>ESTs is publicly accessible at <url>http://mcbc.usm.edu/estmd/</url>.</p

Cloning, Analysis and Functional Annotation of Expressed Sequence Tags from the Earthworm \u3ci\u3eEisenia fetida\u3c/i\u3e

Background Eisenia fetida, commonly known as red wiggler or compost worm, belongs to the Lumbricidae family of the Annelida phylum. Little is known about its genome sequence although it has been extensively used as a test organism in terrestrial ecotoxicology. In order to understand its gene expression response to environmental contaminants, we cloned 4032 cDNAs or expressed sequence tags (ESTs) from two E. fetida libraries enriched with genes responsive to ten ordnance related compounds using suppressive subtractive hybridization-PCR. Results A total of 3144 good quality ESTs (GenBank dbEST accession number EH669363–EH672369 and EL515444–EL515580) were obtained from the raw clone sequences after cleaning. Clustering analysis yielded 2231 unique sequences including 448 contigs (from 1361 ESTs) and 1783 singletons. Comparative genomic analysis showed that 743 or 33% of the unique sequences shared high similarity with existing genes in the GenBank nr database. Provisional function annotation assigned 830 Gene Ontology terms to 517 unique sequences based on their homology with the annotated genomes of four model organisms Drosophila melanogaster, Mus musculus, Saccharomyces cerevisiae, and Caenorhabditis elegans. Seven percent of the unique sequences were further mapped to 99 Kyoto Encyclopedia of Genes and Genomes pathways based on their matching Enzyme Commission numbers. All the information is stored and retrievable at a highly performed, web-based and user-friendly relational database called EST model database or ESTMD version 2. Conclusion The ESTMD containing the sequence and annotation information of 4032 E. fetida ESTs is publicly accessible at http://mcbc.usm.edu/estmd