Expressed Sequence Tag Contigs Research Articles

Transcriptome analysis of extant cotton progenitors revealed tetraploidization and identified genome-specific single nucleotide polymorphism in diploid and allotetraploid cotton.

BackgroundThe most widely cultivated cotton (Gossypium hirsutum L., AD-genome) is derived from tetraploidization between A- and D-genome species. G. arboreum L. (A-genome) and G. raimondii Ulbr. (D-genome) are two of closely-related extant progenitors. Gene expression studies in allotetraploid cotton are complicated by the homoeologous loci of A- and D-genome origins. To develop genomic resources for gene expression and cotton breeding, we sequenced and assembled expressed sequence tags (ESTs) derived from G. arboreum and G. raimondii.ResultsRoche/454 FLX sequencing technology was employed to sequence normalized cDNA libraries prepared from leaves, roots, bolls, ovules, and fibers in G. arboreum and G. raimondii, respectively. Sequencing reads from two independent libraries in each species were combined to assemble high-quality EST contigs. The combined sequencing reads included 1,699,776 from A-genome and 1,464,815 from D-genome, which were clustered into 89,588 contigs in the A-genome and 65,542 contigs in the D-genome. These contigs represented ~80% of EST collections in Cotton Gene Index 11 (CGI11, March 2011). Compared to the D-genome transcript database, 27,537 and 10,452 contigs were unique transcripts in A and D genomes, respectively. Further analysis using self-blastn reduced the unigene contig number by 52% in A-genome and 57% in D-genome, suggesting that 50% or more of contigs are paralogs or isoforms within each species. The majority of EST contigs (73–81%) were conserved between A- and D-genomes, whereas 27% and 19% contigs were specific to A- and D-genomes, respectively. Using these ESTs, we generated a total of 75,754 genome-specific single nucleotide polymorphism (SNP) (gSNPs or GNPs) or homoeologous-specific SNPs (hSNPs) of 10,885 contigs or genes between A and D genomes, indicating a possibility of separating allelic expression for those genes in allotetraploid cotton.ConclusionsExpressed genes are highly redundant within each diploid progenitor and between A and D progenitor species, suggesting that diploid progenitors in cotton are likely ancient tetraploids. This large set of A- and D-genome ESTs and GNPs will be valuable resources for genome annotation, gene expression, and crop improvement in allotetraploid cotton.Electronic supplementary materialThe online version of this article (doi:10.1186/1756-0500-7-493) contains supplementary material, which is available to authorized users.

Development and integration of EST–SSR markers into an established linkage map in switchgrass

Switchgrass (Panicum virgatum L.) is a model cellulosic biofuel crop in the United States. Simple sequence repeat (SSR) markers are valuable resources for genetic mapping and molecular breeding. A large number of expressed sequence tags (ESTs) of switchgrass are recently available in our sequencing project. The objectives of this study were to develop new SSR markers from the switchgrass EST sequences and to integrate them into an existing linkage map. More than 750 unique primer pairs (PPs) were designed from 243,600 EST contigs and tested for PCR amplifications, resulting in 538 PPs effectively producing amplicons of expected sizes. Of the effective PPs, 481 amplifying informative bands in NL94 were screened for polymorphisms in a panel consisting of NL94 and its seven first-generation selfed (S1) progeny. This led to the selection of 117 polymorphic EST–SSRs to genotype a mapping population encompassing 139 S1 individuals of NL94. Of 83 markers demonstrating clearly scorable alleles in the mapping population, 79 were integrated into a published linkage map, with three linked to accessory loci and one unlinked. The newly identified EST–SSR loci were distributed in 17 of 18 linkage groups with 27 (32.5 %) exhibiting distorted segregations. The integration of EST–SSRs aided in reducing the average marker interval (cM) to 3.7 from 4.2, and reduced the number of gaps (each >15 cM) to 10 from 23. Developing new EST–SSRs and constructing a higher density linkage map will facilitate quantitative trait locus mapping and provide a firm footing for marker-assisted breeding in switchgrass.

Identification of multiple genes and their expression profiles in four strains of Oreochromis spp. in response to Streptococcus iniae

Two subtractive complementary DNA libraries were constructed from Nile tilapia Oreochromis niloticus vaccinated with formalin-killed Streptococcus iniae cells, and a further two constructed from O. niloticus infected with S. iniae. Of the 68 distinct expressed sequence tag (EST) contigs and singletons, 45 and 13 EST shared high similarities with genes of known and unknown functions, respectively. Ten EST contigs and singletons had no significant similarity to any sequences. Five putative immune-relevant genes, β2m, α-ha, mmp9, pgrn and cxcr4, were selected for quantitative reverse-transcription polymerase chain reaction in four strains of Oreochromis spp.: genetically improved farmed tilapia (O. niloticus), Oreochromis aureus, O. niloticus and O. niloticus×O. aureus, with different disease resistance following infection with S. iniae. pgrn was up-regulated more significantly in disease-resistant strains than in the susceptible. α-ha was markedly down-regulated, and no significant differences in the expression level of β2m were detected. A negative correlation was observed between the expression of mmp9 and that of cxcr4. The results provide insight into the molecular response of O. niloticus to S. iniae infection.

In silico identification and validation of microsatellite markers from onion EST sequences

SummaryMicrosatellites or simple sequence repeats (SSRs) play an important role in plant genetics and breeding. The development of microsatellites is a time-consuming and expensive process. In silico mining of microsatellites from expressed sequence tags (ESTs), which are available in electronic molecular databases, is a cheaper alternative. In the present study, we used a computational approach for mining SSRs from 20,159 ESTs in Allium cepa. These onion ESTs represented a total length of 13.2 Mb and were downloaded from the dbEST database of the National Center for Biotechnology Information (NCBI) and subjected to various pre-processing steps. The pre-processed ESTs were clustered, resulting in non-redundant unigenes. These unigenes were analysed for their SSR content and distribution. In all, 1,464 SSRs consisting of di-, tri-, tetra-, penta- and hexa-nucleotide repeats were mined from the non-redundant ESTs (contigs and singletons). Tri-nucleotide SSRs were the most abundant, followed by tetra-, di-, hexa- and penta-nucleotide SSRs. Among the tri-coding repeats, leucine and serine codons were more abundant. The SSR-containing sequences were annotated and grouped into their respective functional categories. The predominant functional group among the annotated unigenes was “metabolism”, followed by “transcription factors” and “transporter proteins”. Primer pairs could be designed for 1,092 SSR-containing sequences. Of these, 51 primer pairs were validated in the laboratory. A database has been developed to store the unigenes, primer pairs, putative annotations, and BLAST results. After validation, the EST-derived microsatellite (SSR) markers can be used in studies related to marker-assisted selection, detection of polymorphism, DNA fingerprinting, and diversity analysis in onion.

Transcriptome analysis of enriched Golovinomyces orontii haustoria by deep 454 pyrosequencing

Powdery mildews are phytopathogenic ascomycetes that have an obligate biotrophic lifestyle and establish intimate relationships with their plant hosts. A crucial aspect of this plant–fungus interaction is the formation of specialized fungal infection structures termed haustoria. Although located within the cell boundaries of plant epidermal cells, haustoria remain separated from the plant cytoplasm by a host plasma membrane derivative, the extrahaustorial membrane. Haustoria are thought to represent pivotal sites of nutrient uptake and effector protein delivery. We enriched haustorial complexes from Arabidopsis thaliana plants infected with the powdery mildew fungus Golovinomyces orontii and performed in-depth transcriptome analysis by 454-based pyrosequencing of haustorial cDNAs. We assembled 7077 expressed sequence tag (EST) contigs with greater than 5-fold average coverage and analyzed these with regard to the respective predicted protein functions. We found that transcripts coding for gene products with roles in protein turnover, detoxification of reactive oxygen species and fungal pathogenesis are abundant in the haustorial EST contigs, while surprisingly transcripts encoding presumptive nutrient transporters were not highly represented in the haustorial cDNA library. A substantial proportion (∼38%) of transcripts coding for predicted secreted proteins comprises effector candidates. Our data provide valuable insights into the transcriptome of the key infection structure of a model obligate biotrophic phytopathogen.

Analysis of Porphyra Membrane Transporters Demonstrates Gene Transfer among Photosynthetic Eukaryotes and Numerous Sodium-Coupled Transport Systems

Membrane transporters play a central role in many cellular processes that rely on the movement of ions and organic molecules between the environment and the cell, and between cellular compartments. Transporters have been well characterized in plants and green algae, but little is known about transporters or their evolutionary histories in the red algae. Here we examined 482 expressed sequence tag contigs that encode putative membrane transporters in the economically important red seaweed Porphyra (Bangiophyceae, Rhodophyta). These contigs are part of a comprehensive transcriptome dataset from Porphyra umbilicalis and Porphyra purpurea. Using phylogenomics, we identified 30 trees that support the expected monophyly of red and green algae/plants (i.e. the Plantae hypothesis) and 19 expressed sequence tag contigs that show evidence of endosymbiotic/horizontal gene transfer involving stramenopiles. The majority (77%) of analyzed contigs encode transporters with unresolved phylogenies, demonstrating the difficulty in resolving the evolutionary history of genes. We observed molecular features of many sodium-coupled transport systems in marine algae, and the potential for coregulation of Porphyra transporter genes that are associated with fatty acid biosynthesis and intracellular lipid trafficking. Although both the tissue-specific and subcellular locations of the encoded proteins require further investigation, our study provides red algal gene candidates associated with transport functions and novel insights into the biology and evolution of these transporters.

Comparative transcriptome analysis of Pima and Acala cotton during boll development using 454 pyrosequencing technology

Information on gene expression during cotton fiber development for high fiber quality genotypes including Pima (Gossypium barbadense L.) and Acala cotton (G. hirsutum L.) is lacking, impeding the discovery of novel genes and DNA markers for fiber quality traits. In the present study, massively parallel sequencing technology was used to generate a total of 394,301 expressed sequence tag (EST) sequence reads, totaling over 42 million bp, from immature ovaries at 10 days post-anthesis of Pima Phy 76 and Acala 1517-99. For Pima Phy 76, 1,900 EST contigs were tentatively assembled from 21.8% (50,876) ESTs which resulted in 1,541 tentative consensus sequences (TCs). 1,005 and 1,064 EST contigs were homologous to TCs from the A2 and D5 genomes, respectively, while 738 of the EST contigs were shared by both genomes. For Acala 1517-99, 2,151 EST contigs were tentatively assembled from 24.5% (39,368) ESTs and resulted in 1,787 TCs. 1,220 and 1,206 EST contigs were homologous to known transcripts from A2 and D5 genomes and only 452 of the EST contigs were shared by both genomes. 691 EST contigs are also homologous (E-value ≤ 1 × 10−10) between Pima Phy 76 and Acala 1517-99, of which 38 contigs pairs were highly homologous (E-value ≤ 1 × 10−60) with sequence variations, predominantly single nucleotide polymorphisms (SNPs). A Blast search of a sample of contigs to Upland cotton ESTs in GenBank indicated that more than 75% EST contigs may contain SNPs. This report represents the first research in cotton using 454 pyrosequencing technology to enrich ESTs for high fiber quality cotton.

Dozens of Toxin-Related Genes Are Expressed in a Nontoxic Strain of the Dinoflagellate Heterocapsa circularisquama

The dinoflagellate Heterocapsa circularisquama is lethal to a variety of marine organisms, in particular, commercially important farmed bivalves. Unlike most dinoflagellate toxins, which are polyketides, the only described toxin from H. circularisquama (H2-a) is a porphyrin derivative that functions in light. It is unknown whether H2-a is produced specifically for its lytic properties. We searched for toxin-related genes in the transcriptome of a nontoxic strain of H. circularisquama, and surprisingly found the richest set of toxin-related genes yet described in dinoflagellates. There are 87 distinct expressed sequence tag contigs that encode polyketide synthases and nonribosomal peptide synthases, as well as 8 contigs that are involved in porphyrin biosynthesis. Phylogenomic analysis shows that many toxin-related genes are widely distributed among dinoflagellates. Our data likely indicate a variety of unknown metabolic functions for the toxin-related genes in H. circularisquama because they were identified in a nontoxic strain raised in unialgal culture.

Evaluation of 15 caespitose bamboo EST‐SSR markers for cross‐species/genera transferability and ability to identify interspecies hybrids

With 1 figure and 1 tableAbstractPublic sequence databases provide a convenient and inexpensive resource for the rapid development of microsatellite markers. We analysed 3406 publically available expressed sequence tags (ESTs) from caespitose bamboo species (Bambusa edulis and B. oldhamii) and found 245 non‐redundant simple sequence repeats (SSRs) in 205 EST contigs that were used to develop 15 EST‐SSR markers. The transferability of the markers was 59.6% among 14 additional caespitose bamboo species (4 within the same genus and 10 in other genera) and Phyllostachys pubescens. The successfully transferred markers showed 51.4% polymorphism. Markers BOM01 and BOM02 transferred successfully to most of the caespitose bamboo species showed rich polymorphism, and are therefore potentially valuable as species‐specific alleles for the identification of caespitose bamboo interspecies hybrids.

EST contig-based SSR linkage maps for Malus × domestica cv Royal Gala and an apple scab resistant accession of M. sieversii, the progenitor species of domestic apple

Malus sieversii is a progenitor species of domestic apple M. × domestica. Using population “GMAL 4595” of 188 individuals derived from a cross of Royal Gala × PI 613988 (apple scab resistant, M. sieversii), 287 SSR (simple sequence repeats) loci were mapped. Of these SSRs, 80 are published anchors and 207 are newly developed EST (expressed sequence tag) contig-based SSRs, representing 1,630 Malus EST accessions in GenBank. Putative gene functions of these EST contigs are diverse, including regulating plant growth, development and response to environmental stresses. Among the 80 published SSRs, 18 are PI 613988 specific, 38 are common and 24 are Royal Gala specific. Out of the 207 newly developed EST contig-based SSRs, 79 are PI 613988 specific, 45 are common and 83 are Royal Gala specific. These results led to the construction of a M. sieversii map (1,387.0 cM) of 180 SSR markers and a Royal Gala map (1,283.4 cM) of 190 SSR markers. Mapping of scab resistance was independently conducted in two subsets of population “GMAL 4595” that were inoculated with Ventura inaequalis races (1) and (2), respectively. In combination with the two major resistance reactions Chl (chlorotic lesions) and SN (stellate necrosis) to each race, four subsets of resistance data, i.e., Chl/race (1), SN/race (1), Chl/race (2) and SN/race (2), were constituted and analyzed, leading to four resistance loci mapped to the linkage group 2 of PI 613988; SNR1 (stellate necrosis resistance to race (1)) and SNR2 are tightly linked in a region of known scab resistance genes, and ChlR1 (Chlorotic lesion resistance to race (1)) and ChlR2 are also linked tightly but in a region without known scab resistance genes. The utility of the two linkage maps, the new EST contig-based markers and M. sieversii as sources of apple scab resistance are discussed.

Red and problematic green phylogenetic signals among thousands of nuclear genes from the photosynthetic and apicomplexa-related Chromera velia.

The photosynthetic and basal apicomplexan Chromera velia was recently described, expanding the membership of this otherwise nonphotosynthetic group of parasite protists. Apicomplexans are alveolates with secondary plastids of red algal origin, but the evolutionary history of their nuclear genes is still actively discussed. Using deep sequencing of expressed genes, we investigated the phylogenetic affinities of a stringent filtered set of 3,151 expressed sequence tag-contigs by generating clusters with eukaryotic homologs and constructing phylogenetic trees and networks. The phylogenetic positioning of this alveolate alga was determined and sets of phyla-specific proteins extracted. Phylogenetic trees provided conflicting signals, with 444 trees grouping C. velia with the apicomplexans but 354 trees grouping C. velia with the alveolate oyster pathogen Perkinsus marinus, the latter signal being reinforced from the analysis of shared genes and overall sequence similarity. Among the 513 C. velia nuclear genes that reflect a photosynthetic ancestry and for which nuclear homologs were available both from red and green lineages, 263 indicated a red photosynthetic ancestry, whereas 250 indicated a green photosynthetic ancestry. The same 1:1 signal ratio was found among the putative 255 nuclear-encoded plastid proteins identified. This finding of red and green signals for the alveolate mirrors the result observed in the heterokont lineage and supports a common but not necessarily single origin for the plastid in heterokonts and alveolates. The inference of green endosymbiosis preceding red plastid acquisition in these lineages leads to worryingly complicated evolutionary scenarios, prompting the search for other explanations for the green phylogenetic signal and the amount of hosts involved.

The Application and Performance of Single Nucleotide Polymorphism Markers for Population Genetic Analyses of Lepidoptera

Microsatellite markers are difficult to apply within lepidopteran studies due to the lack of locus-specific PCR amplification and the high proportion of “null” alleles, such that erroneous estimations of population genetic parameters often result. Herein single nucleotide polymorphism (SNP) markers are developed from Ostrinia nubilalis (Lepidoptera: Crambidae) using next generation expressed sequence tag (EST) data. A total of 2742 SNPs were predicted within a reference assembly of 7414 EST contigs, and a subset of 763 were incorporated into 24 multiplex PCR reactions. To validate this pipeline, 5 European and North American sample sites were genotyped at 178 SNP loci, which indicated 84 (47.2%) were in Hardy–Weinberg equilibrium. Locus-by-locus FST, analysis of molecular variance, and STRUCTURE analyses indicate significant genetic differentiation may exist between European and North American O. nubilalis. The observed genetic diversity was significantly lower among European sites, which may result from genetic drift, natural selection, a genetic bottleneck, or ascertainment bias due to North American origin of EST sequence data. SNPs are an abundant source of mutation data for molecular genetic marker development in non-model species, with shared ancestral SNPs showing application within closely related species. These markers offer advantages over microsatellite markers for genetic and genomic analyses of Lepidoptera, but the source of mutation data may affect the estimation of population parameters and likely need to be considered in the interpretation of empirical data.

In silico single nucleotide polymorphism discovery and application to marker-assisted selection in soybean

The general approach to discovering single nucleotide polymorphisms (SNPs) requires locus-specific PCR amplification. To enhance the efficiency of SNP discovery in soybean, we used in silico analysis prior to re-sequencing as it is both rapid and inexpensive. In silico analysis was performed to detect putative SNPs in expressed sequence tag (EST) contigs assembled using publicly available ESTs from 18 different soybean genotypes. SNP validation by direct sequencing of six soybean cultivars and a wild soybean genotype was performed with PCR primers designed from EST contigs aligned with at least 5 out of 18 soybean genotypes. The efficiency of SNP discovery among the confirmation genotypes was 81.2%. Furthermore, the efficiency of SNP discovery between Pureunkong and Jinpumkong 2 genotypes was 47.4%, a great improvement on our previous finding based on direct sequencing (22.3%). Using SNPs between Pureunkong and Jinpumkong 2 in EST contigs, which were linked to target traits, we were able to genotype 90 recombinant inbred lines by high-resolution melting (HRM) analysis. These SNPs were mapped onto the expected locations near quantitative trait loci for water-logging tolerance and seed pectin concentration. Thus, our protocol for HRM analysis can be applied successfully not only to genetic diversity studies, but also to marker-assisted selection (MAS). Our study suggests that a combination of in silico analysis and HRM can reduce the cost and labor involved in developing SNP markers and genotyping SNPs. The markers developed in this study can also easily be applied to MAS if the markers are associated with the target traits.

Rapid evolution and selection inferred from the transcriptomes of sympatric crater lake cichlid fishes

Crater lakes provide a natural laboratory to study speciation of cichlid fishes by ecological divergence. Up to now, there has been a dearth of transcriptomic and genomic information that would aid in understanding the molecular basis of the phenotypic differentiation between young species. We used next-generation sequencing (Roche 454 massively parallel pyrosequencing) to characterize the diversity of expressed sequence tags between ecologically divergent, endemic and sympatric species of cichlid fishes from crater lake Apoyo, Nicaragua: benthic Amphilophus astorquii and limnetic Amphilophus zaliosus. We obtained 24 174 A. astorquii and 21 382 A. zaliosus high-quality expressed sequence tag contigs, of which 13 106 pairs are orthologous between species. Based on the ratio of nonsynonymous to synonymous substitutions, we identified six sequences exhibiting signals of strong diversifying selection (K(a)/K(s) > 1). These included genes involved in biosynthesis, metabolic processes and development. This transcriptome sequence variation may be reflective of natural selection acting on the genomes of these young, sympatric sister species. Based on Ks ratios and p-distances between 3'-untranslated regions (UTRs) calibrated to previously published species divergence times, we estimated a neutral transcriptome-wide substitutional mutation rate of approximately 1.25 x 10(-6) per site per year. We conclude that next-generation sequencing technologies allow us to infer natural selection acting to diversify the genomes of young species, such as crater lake cichlids, with much greater scope than previously possible.

In Silico and Quantitative Analyses of MADS-Box Genes in Coffea arabica

MADS-box genes comprise a family of transcription factors that act as key regulators in many cellular development processes of several organisms. Members of this family have highly conserved regions and play important roles as transcription factors, activating target genes. The present work aimed at analyzing the MADS-box gene family present in a database generated by the Brazilian Coffee Genome Project (CAFEST) as well as of observing their respective sites of expression within these data. Through the use of bioinformatics tools, it was possible to identify and classify 26 expressed sequence tag contigs, 13 of them were expressed exclusively in vegetative tissues, 11 in reproductive tissues, and two in both. Later, quantitative analysis by quantitative reverse transcriptase PCR was carried out for three of them belonging to the groups of genes APETALA3 (B genes), AGAMOUS (C genes), and SEPALLATAS (E genes), which could be compared with the expression profile of in silico analysis and of its putative orthologous genes. Therefore, this study aims at a better understanding of the development processes performed by this family of genes in Coffea arabica, mainly in reproductive organs, as well as to compare functions of MADS-box orthologous studied in other species.

Generation and analysis of transcriptomic resources for a model system on the rise: the sea anemone Aiptasia pallida and its dinoflagellate endosymbiont

BackgroundThe most diverse marine ecosystems, coral reefs, depend upon a functional symbiosis between cnidarian hosts and unicellular dinoflagellate algae. The molecular mechanisms underlying the establishment, maintenance, and breakdown of the symbiotic partnership are, however, not well understood. Efforts to dissect these questions have been slow, as corals are notoriously difficult to work with. In order to expedite this field of research, we generated and analyzed a collection of expressed sequence tags (ESTs) from the sea anemone Aiptasia pallida and its dinoflagellate symbiont (Symbiodinium sp.), a system that is gaining popularity as a model to study cellular, molecular, and genomic questions related to cnidarian-dinoflagellate symbioses.ResultsA set of 4,925 unique sequences (UniSeqs) comprising 1,427 clusters of 2 or more ESTs (contigs) and 3,498 unclustered ESTs (singletons) was generated by analyzing 10,285 high-quality ESTs from a mixed host/symbiont cDNA library. Using a BLAST-based approach to predict which unique sequences derived from the host versus symbiont genomes, we found that the contribution of the symbiont genome to the transcriptome was surprisingly small (1.6–6.4%). This may reflect low levels of gene expression in the symbionts, low coverage of alveolate genes in the sequence databases, a small number of symbiont cells relative to the total cellular content of the anemones, or failure to adequately lyse symbiont cells. Furthermore, we were able to identify groups of genes that are known or likely to play a role in cnidarian-dinoflagellate symbioses, including oxidative stress pathways that emerged as a prominent biological feature of this transcriptome. All ESTs and UniSeqs along with annotation results and other tools have been made accessible through the implementation of a publicly accessible database named AiptasiaBase.ConclusionWe have established the first large-scale transcriptomic resource for Aiptasia pallida and its dinoflagellate symbiont. These data provide researchers with tools to study questions related to cnidarian-dinoflagellate symbioses on a molecular, cellular, and genomic level. This groundwork represents a crucial step towards the establishment of a tractable model system that can be utilized to better understand cnidarian-dinoflagellate symbioses. With the advent of next-generation sequencing methods, the transcriptomic inventory of A. pallida and its symbiont, and thus the extent of AiptasiaBase, should expand dramatically in the near future.

GarlicESTdb: an online database and mining tool for garlic EST sequences

BackgroundAllium sativum., commonly known as garlic, is a species in the onion genus (Allium), which is a large and diverse one containing over 1,250 species. Its close relatives include chives, onion, leek and shallot. Garlic has been used throughout recorded history for culinary, medicinal use and health benefits. Currently, the interest in garlic is highly increasing due to nutritional and pharmaceutical value including high blood pressure and cholesterol, atherosclerosis and cancer. For all that, there are no comprehensive databases available for Expressed Sequence Tags(EST) of garlic for gene discovery and future efforts of genome annotation. That is why we developed a new garlic database and applications to enable comprehensive analysis of garlic gene expression.DescriptionGarlicESTdb is an integrated database and mining tool for large-scale garlic (Allium sativum) EST sequencing. A total of 21,595 ESTs collected from an in-house cDNA library were used to construct the database. The analysis pipeline is an automated system written in JAVA and consists of the following components: automatic preprocessing of EST reads, assembly of raw sequences, annotation of the assembled sequences, storage of the analyzed information into MySQL databases, and graphic display of all processed data. A web application was implemented with the latest J2EE (Java 2 Platform Enterprise Edition) software technology (JSP/EJB/JavaServlet) for browsing and querying the database, for creation of dynamic web pages on the client side, and for mapping annotated enzymes to KEGG pathways, the AJAX framework was also used partially. The online resources, such as putative annotation, single nucleotide polymorphisms (SNP) and tandem repeat data sets, can be searched by text, explored on the website, searched using BLAST, and downloaded. To archive more significant BLAST results, a curation system was introduced with which biologists can easily edit best-hit annotation information for others to view. The GarlicESTdb web application is freely available at .ConclusionGarlicESTdb is the first incorporated online information database of EST sequences isolated from garlic that can be freely accessed and downloaded. It has many useful features for interactive mining of EST contigs and datasets from each library, including curation of annotated information, expression profiling, information retrieval, and summary of statistics of functional annotation. Consequently, the development of GarlicESTdb will provide a crucial contribution to biologists for data-mining and more efficient experimental studies.

NemaPath: online exploration of KEGG-based metabolic pathways for nematodes

BackgroundNematode.net is a web-accessible resource for investigating gene sequences from parasitic and free-living nematode genomes. Beyond the well-characterized model nematode C. elegans, over 500,000 expressed sequence tags (ESTs) and nearly 600,000 genome survey sequences (GSSs) have been generated from 36 nematode species as part of the Parasitic Nematode Genomics Program undertaken by the Genome Center at Washington University School of Medicine. However, these sequencing data are not present in most publicly available protein databases, which only include sequences in Swiss-Prot. Swiss-Prot, in turn, relies on GenBank/Embl/DDJP for predicted proteins from complete genomes or full-length proteins.DescriptionHere we present the NemaPath pathway server, a web-based pathway-level visualization tool for navigating putative metabolic pathways for over 30 nematode species, including 27 parasites. The NemaPath approach consists of two parts: 1) a backend tool to align and evaluate nematode genomic sequences (curated EST contigs) against the annotated Kyoto Encyclopedia of Genes and Genomes (KEGG) protein database; 2) a web viewing application that displays annotated KEGG pathway maps based on desired confidence levels of primary sequence similarity as defined by a user. NemaPath also provides cross-referenced access to nematode genome information provided by other tools available on Nematode.net, including: detailed NemaGene EST cluster information; putative translations; GBrowse EST cluster views; links from nematode data to external databases for corresponding synonymous C. elegans counterparts, subject matches in KEGG's gene database, and also KEGG Ontology (KO) identification.ConclusionThe NemaPath server hosts metabolic pathway mappings for 30 nematode species and is available on the World Wide Web at . The nematode source sequences used for the metabolic pathway mappings are available via FTP , as provided by the Genome Center at Washington University School of Medicine.

Chasing Migration Genes: A Brain Expressed Sequence Tag Resource for Summer and Migratory Monarch Butterflies (Danaus plexippus)

North American monarch butterflies (Danaus plexippus) undergo a spectacular fall migration. In contrast to summer butterflies, migrants are juvenile hormone (JH) deficient, which leads to reproductive diapause and increased longevity. Migrants also utilize time-compensated sun compass orientation to help them navigate to their overwintering grounds. Here, we describe a brain expressed sequence tag (EST) resource to identify genes involved in migratory behaviors. A brain EST library was constructed from summer and migrating butterflies. Of 9,484 unique sequences, 6068 had positive hits with the non-redundant protein database; the EST database likely represents ∼52% of the gene-encoding potential of the monarch genome. The brain transcriptome was cataloged using Gene Ontology and compared to Drosophila. Monarch genes were well represented, including those implicated in behavior. Three genes involved in increased JH activity (allatotropin, juvenile hormone acid methyltransfersase, and takeout) were upregulated in summer butterflies, compared to migrants. The locomotion-relevant turtle gene was marginally upregulated in migrants, while the foraging and single-minded genes were not differentially regulated. Many of the genes important for the monarch circadian clock mechanism (involved in sun compass orientation) were in the EST resource, including the newly identified cryptochrome 2. The EST database also revealed a novel Na+/K+ ATPase allele predicted to be more resistant to the toxic effects of milkweed than that reported previously. Potential genetic markers were identified from 3,486 EST contigs and included 1599 double-hit single nucleotide polymorphisms (SNPs) and 98 microsatellite polymorphisms. These data provide a template of the brain transcriptome for the monarch butterfly. Our “snap-shot” analysis of the differential regulation of candidate genes between summer and migratory butterflies suggests that unbiased, comprehensive transcriptional profiling will inform the molecular basis of migration. The identified SNPs and microsatellite polymorphisms can be used as genetic markers to address questions of population and subspecies structure.

EGENES: transcriptome-based plant database of genes with metabolic pathway information and expressed sequence tag indices in KEGG.

EGENES is a knowledge-based database for efficient analysis of plant expressed sequence tags (ESTs) that was recently added to the KEGG suite of databases. It links plant genomic information with higher order functional information in a single database. It also provides gene indices for each genome. The genomic information in EGENES is a collection of EST contigs constructed from assembly of ESTs. Due to the extremely large genomes of plant species, the bulk collection of data such as ESTs is a quick way to capture a complete repertoire of genes expressed in an organism. Using ESTs for reconstructing metabolic pathways is a new expansion in KEGG and provides researchers with a new resource for species in which only EST sequences are available. Functional annotation in EGENES is a process of linking a set of genes/transcripts in each genome with a network of interacting molecules in the cell. EGENES is a multispecies, integrated resource consisting of genomic, chemical, and network information containing a complete set of building blocks (genes and molecules) and wiring diagrams (biological pathways) to represent cellular functions. Using EGENES, genome-based pathway annotation and EST-based annotation can now be compared and mutually validated. The ultimate goals of EGENES will be to: bring new plant species into KEGG by clustering and annotating ESTs; abstract knowledge and principles from large-scale plant EST data; and improve computational prediction of systems of higher complexity. EGENES will be updated at least once a year. EGENES is publicly available and is accessible by the following link or by KEGG's navigation system (http://www.genome.jp/kegg-bin/create_kegg_menu?category=plants_egenes).