EST-PAC HPC - a web portal for high-throughput EST annotation and protein sequence prediction

Expressed Sequence Tags (ESTs) are short DNA sequences generated by sequencing the transcribed cDNAs coming from a gene expression. They can provide significant functional, structural and evolutionary information and thus are a primary resource for gene discovery. EST annotation basically refers to the analysis of unknown ESTs that can be performed by database similarity search for possible identities and database search for functional prediction of translation products. Such kind of annotation typically consists of a series of repetitive tasks which should be automated, and be customizable and amenable to using distributed computing resources. Furthermore, processing of EST data should be done efficiently using a high performance computing platform. In this paper, we describe an EST annotator, EST-PACHPC, which has been developed for harnessing HPC resources potentially from Grid and Cloud systems for high throughput EST annotations. The performance analysis of EST-PACHPC has shown that it provides substantial performance gain in EST annotation.

Eucalyptus ESTs involved in mechanisms against plant pathogens and environmental stresses

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Identification of human chromosome 22 transcribed sequences with ORF expressed sequence tags

Transcribed sequences in the human genome can be identified with confidence only by alignment with sequences derived from cDNAs synthesized from naturally occurring mRNAs. We constructed a set of 250,000 cDNAs that represent partial expressed gene sequences and that are biased toward the central coding regions of the resulting transcripts. They are termed ORF expressed sequence tags (ORESTES). The 250,000 ORESTEs were assembled into 81,429 contigs. of these, 1,181 (1.45%) were found to match sequences in chromosome 22 with at least one ORESTES contig for 162 (65.6%) of the 247 known genes, for 67 (44.6%) of the 150 related genes, and for 45 of the 148 (30.4%) EST-predicted genes on this chromosome. Using a set of stringent criteria to validate our sequences, we identified a further 219 previously unannotated transcribed sequences on chromosome 22. of these, 171 were in fact also defined by EST or full length cDNA sequences available in GenBank but not utilized in the initial annotation of the first human chromosome sequence. Thus despite representing less than 15% of all expressed human sequences in the public databases at the time of the present analysis, ORESTEs sequences defined 48 transcribed sequences on chromosome 22 not defined by other sequences. All of the transcribed sequences defined by ORESTEs coincided with DNA regions predicted as encoding exons by GENSCAN.

Expressed sequence tags from Atta laevigata and identification of candidate genes for the control of pest leaf-cutting ants.

Leafcutters are the highest evolved within Neotropical ants in the tribe Attini and model systems for studying caste formation, labor division and symbiosis with microorganisms. Some species of leafcutters are agricultural pests controlled by chemicals which affect other animals and accumulate in the environment. Aiming to provide genetic basis for the study of leafcutters and for the development of more specific and environmentally friendly methods for the control of pest leafcutters, we generated expressed sequence tag data from Atta laevigata, one of the pest ants with broad geographic distribution in South America. Results: The analysis of the expressed sequence tags allowed us to characterize 2,006 unique sequences in Atta laevigata. Sixteen of these genes had a high number of transcripts and are likely positively selected for high level of gene expression, being responsible for three basic biological functions: energy conservation through redox reactions in mitochondria; cytoskeleton and muscle structuring; regulation of gene expression and metabolism. Based on leafcutters lifestyle and reports of genes involved in key processes of other social insects, we identified 146 sequences potential targets for controlling pest leafcutters. The targets are responsible for antixenobiosis, development and longevity, immunity, resistance to pathogens, pheromone function, cell signaling, behavior, polysaccharide metabolism and arginine kynase activity. Conclusion: The generation and analysis of expressed sequence tags from Atta laevigata have provided important genetic basis for future studies on the biology of leaf-cutting ants and may contribute to the development of a more specific and environmentally friendly method for the control of agricultural pest leafcutters.

Eucalyptus ESTs associated with resistance to herbicide inhibitors of aromatic and branched-chain amino acid synthesis

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Gene2EST: a BLAST2 server for searching expressed sequence tag (EST) databases with eukaryotic gene-sized queries

Expressed sequence tags (ESTs) are randomly sequenced cDNA clones. Currently, nearly 3 million human and 2 million mouse ESTs provide valuable resources that enable researchers to investigate the products of gene expression. The EST databases have proven to be useful tools for detecting homologous genes, for exon mapping, revealing differential splicing, etc. With the increasing availability of large amounts of poorly characterised eukaryotic (notably human) genomic sequence, ESTs have now become a vital tool for gene identification, sometimes yielding the only unambiguous evidence for the existence of a gene expression product. However, BLAST-based Web servers available to the general user have not kept pace with these developments and do not provide appropriate tools for querying EST databases with large highly spliced genes, often spanning 50 000–100 000 bases or more. Here we describe Gene2EST (http://woody.embl-heidelberg.de/gene2est/), a server that brings together a set of tools enabling efficient retrieval of ESTs matching large DNA queries and their subsequent analysis. RepeatMasker is used to mask dispersed repetitive sequences (such as Alu elements) in the query, BLAST2 for searching EST databases and Artemis for graphical display of the findings. Gene2EST combines these components into a Web resource targeted at the researcher who wishes to study one or a few genes to a high level of detail.

Expressed Sequence Tags from a Root-Hair-Enriched Medicago truncatula cDNA Library1

The root hair is a specialized cell type involved in water and nutrient uptake in plants. In legumes the root hair is also the primary site of recognition and infection by symbiotic nitrogen-fixing Rhizobium bacteria. We have studied the root hairs of Medicago truncatula, which is emerging as an increasingly important model legume for studies of symbiotic nodulation. However, only 27 genes from M. truncatula were represented in GenBank/EMBL as of October, 1997. We report here the construction of a root-hair-enriched cDNA library and single-pass sequencing of randomly selected clones. Expressed sequence tags (899 total, 603 of which have homology to known genes) were generated and made available on the Internet. We believe that the database and the associated DNA materials will provide a useful resource to the community of scientists studying the biology of roots, root tips, root hairs, and nodulation.

Investigation of molecular mechanisms regulating biomineralization of pearl oyster Pinctada maxima

Biomineralization is a process encompassing all mineral containing tissues produced within an organism. The most dynamic example of this process is the formation of the mollusk shell, comprising a variety of crystal phases and microstructures. The organic component incorporated within the shell is said to dictate this remarkable architecture. Subsequently, for the past decade considerable research have been undertaken to identify and characterize the protein components involved in biomineralization. Despite these efforts the general understanding of the process remains ambiguous. This study employs a novel molecular approach to further the elucidation of the shell biomineralization. A microarray platform has been custom generated (PmaxArray 1.0) from the pearl oyster Pinctada maxima. PmaxArray 1.0 consists of 4992 expressed sequence tags (ESTs) originating from the mantle, an organ involved in shell formation. This microarray has been used as the primary tool for three separate investigations in an effort to associate transcriptional gene expression from P. maxima to the process of shell biomineralization. The first investigation analyzes the spatial expression of ESTs throughout the mantle organ. The mantle was dissected into five discrete regions and each analyzed for gene expression with PmaxArray 1.0. Over 2000 ESTs were differentially expressed among the tissue sections, identifying five major expression regions. Three of these regions have been proposed to have shell formation functions belonging to nacre, prismatic calcite and periostracum. The spatial gene expression map was confirmed by in situ hybridization, localizing a subset of ESTs from each expression region to the same mantle area. Comparative sequence analysis of ESTs expressed in the proposed shell formation regions with the BLAST tool, revealed a number of the transcripts were novel while others showed significant sequence similarities to previously characterized shell formation genes. The second investigation correlates temporal EST expression during P. maxima larval ontogeny with transitions in shell mineralization during the same period. A timeline documenting the morphologicat microstructural and mineralogical shell characteristics of P. maxima throughout larval ontogeny has been established. Three different shell types were noted based on the physical characters and termed, prodissoconch I, prodissoconch 11 and dissoconch. PmaxArray 1.0 analyzed ESTs expression of animals throughout the larval development of P. maxima, noting up-regulation of 359 ESTs in association with the shell transitions from prodissoconch 1 to prodissoconch 11 to dissoconch. Comparative sequence analysis of these ESTs indicates a number of the transcripts are novel as well as showing significant sequence similarities between ESTs and known shell matrix associated genes and proteins. These ESTs are discussed in relation to the shell characters associated with their temporal expression. The third investigation uses PmaxArray 1.0 to analyze gene expression in the mantle tissue of P. maxima specimens exposed to sub-lethal concentrations of a shell-deforming toxin, tributyltin (TBT). The shell specific effects of TBT are used in this investigation to interpret differential expression of ESTs with respect to shell formation functions. A lethal and sublethal TBT concentration range was established for P. maxima, noting a concentration of 50 ng L- 1 TBT as sub-lethal over a 21 day period. Mantle tissue from P. maxima animals treated with 50 ng L- 1 TBT was assessed for differential EST expression with untreated control animals. A total of 102 ESTs were identified as differentially expressed in association with TBT exposure, comparative sequence identities included an up-regulation of immunity and detoxification related genes and down-regulation of several shell matrix genes. A number of transcripts encoding novel peptides were additionally identified. The potential actions of these genes are discussed with reference to TBT toxicity and shell biomineralization. This thesis has used a microarray platform to analyze gene expression in spatial, temporal and toxicity investigations, revealing the involvement of numerous gene transcripts in specific shell formation functions. Investigation of thousands of transcripts simultaneously has provided a holistic interpretation of the organic components regulating shell biomineralization.

Spatial analysis of biomineralization associated gene expression from the mantle organ of the pearl oyster Pinctada maxima

Background: Biomineralization is a process encompassing all mineral containing tissues produced within an organism. One of the most dynamic examples of this process is the formation of the mollusk shell, comprising a variety of crystal phases and microstructures. The organic component incorporated within the shell is said to dictate this architecture. However general understanding of how this process is achieved remains ambiguous. The mantle is a conserved organ involved in shell formation throughout molluscs. Specifically the mantle is thought to be responsible for secreting the protein component of the shell. This study employs molecular approaches to determine the spatial expression of genes within the mantle tissue to further the elucidation of the shell biomineralization. Results: A microarray platform was custom generated (PmaxArray 1.0) from the pearl oyster Pinctada maxima. PmaxArray 1.0 consists of 4992 expressed sequence tags (ESTs) originating from mantle tissue. This microarray was used to analyze the spatial expression of ESTs throughout the mantle organ. The mantle was dissected into five discrete regions and analyzed for differential gene expression with PmaxArray 1.0. Over 2000 ESTs were determined to be differentially expressed among the tissue sections, identifying five major expression regions. In situ hybridization validated and further localized the expression for a subset of these ESTs. Comparative sequence similarity analysis of these ESTs revealed a number of the transcripts were novel while others showed significant sequence similarities to previously characterized shell related genes.

Effect of 5′UTR introns on gene expression in Arabidopsis thaliana

Background The majority of introns in gene transcripts are found within the coding sequences (CDSs). A small but significant fraction of introns are also found to reside within the untranslated regions (5′UTRs and 3′UTRs) of expressed sequences. Alignment of the whole genome and expressed sequence tags (ESTs) of the model plant Arabidopsis thaliana has identified introns residing in both coding and non-coding regions of the genome. Results A bioinformatic analysis revealed some interesting observations: (1) the density of introns in 5′UTRs is similar to that in CDSs but much higher than that in 3′UTRs; (2) the 5′UTR introns are preferentially located close to the initiating ATG codon; (3) introns in the 5′UTRs are, on average, longer than introns in the CDSs and 3′UTRs; and (4) 5′UTR introns have a different nucleotide composition to that of CDs and 3′UTR introns. Furthermore, we show that the 5′UTR intron of the A. thaliana EFIα-A3 gene affects the gene expression and the size of the 5′UTR intron influences the level of gene expression. Conclusion Introns within the 5′UTR show specific features that distinguish them from introns that reside within the coding sequence and the 3′UTR. In the EFIα-A3 gene, the presence of a long intron in the 5′UTR is sufficient to enhance gene expression in plants in a size dependent manner.

Functionally associated molecular genetic marker map construction in perennial ryegrass (Lolium perenne L.)

A molecular marker-based map of perennial ryegrass (Lolium perenne L.) has been constructed through the use of polymorphisms associated with expressed sequence tags (ESTs). A pair-cross between genotypes from a North African ecotype and the cultivar Aurora was used to generate a two-way pseudo-testcross population. A selection of 157 cDNAs assigned to eight different functional categories associated with agronomically important biological processes was used to detect polymorphic EST–RFLP loci in the F1(NA6 × AU6) population. A comprehensive set of EST–SSR markers was developed from the analysis of 14,767 unigenes, with 310 primer pairs showing efficient amplification and detecting 113 polymorphic loci. Two parental genetic maps were produced: the NA6 genetic map contains 88 EST–RFLP and 71 EST–SSR loci with a total map length of 963 cM, while the AU6 genetic map contains 67 EST–RFLP and 58 EST–SSR loci with a total map length of 757 cM. Bridging loci permitted the alignment of homologous chromosomes between the parental maps, and a sub-set of genomic DNA-derived SSRs was used to relate linkage groups to the perennial ryegrass reference map. Regions of segregation distortion were identified, in some instances in common with other perennial ryegrass maps. The EST-derived marker-based map provides the basis for in silico comparative genetic mapping, as well as the evaluation of co-location between QTLs and functionally associated genetic loci.

Quantitative and molecular genetics of juvenile wood traits in radiata and slash/Caribbean pines.

The Juvenile Wood Initiative (JWI) project has been running successfully since July 2003 under a Research Agreement with FWPA and Letters of Association with the consortium partners STBA (Southern Tree Breeding Association), ArborGen and FPQ (Forestry Plantations Queensland). Over the last five and half years, JWI scientists in CSIRO, FPQ, and STBA have completed all 12 major milestones and 28 component milestones according to the project schedule. We have made benchmark progress in understanding the genetic control of wood formation and interrelationships among wood traits. The project has made 15 primary scientific findings and several results have been adopted by industry as summarized below. This progress was detailed in 10 technical reports to funding organizations and industry clients. Team scientists produced 16 scientific manuscripts (8 published, 1 in press, 2 submitted, and several others in the process of submission) and 15 conference papers or presentations. Primary Scientific Findings. The 15 major scientific findings related to wood science, inheritance and the genetic basis of juvenile wood traits are: 1. An optimal method to predict stiffness of standing trees in slash/Caribbean pine is to combine gravimetric basic density from 12 mm increment cores with a standing tree prediction of MoE using a time of flight acoustic tool. This was the most accurate and cheapest way to rank trees for breeding selection for slash/Caribbean hybrid pine. This method was also recommended for radiata pine. 2. Wood density breeding values were predicted for the first time in the STBA breeding population using a large sample of 7,078 trees (increment cores) and it was estimated that selection of the best 250 trees for deployment will produce wood density gains of 12.4%. 3. Large genetic variation for a suite of wood quality traits including density, MFA, spiral grain, shrinkage, acoustic and non-acoustic stiffness (MoE) for clear wood and standing trees were observed. Genetic gains of between 8 and 49% were predicted for these wood quality traits with selection intensity between 1 to 10% for radiata pine. 4. Site had a major effect on juvenile-mature wood transition age and the effect of selective breeding for a shorter juvenile wood formation phase was only moderate (about 10% genetic gain with 10% selection intensity, equivalent to about 2 years reduction of juvenile wood). 5. The study found no usable site by genotype interactions for the wood quality traits of density, MFA and MoE for both radiata and slash/Caribbean pines, suggesting that assessment of wood properties on one or two sites will provide reliable estimates of the genetic worth of individuals for use in future breeding. 6. There were significant and sizable genotype by environment interactions between the mainland and Tasmanian regions and within Tasmania for DBH and branch size. 7. Strong genetic correlations between rings for density, MFA and MoE for both radiata and slash/Caribbean pines were observed. This suggests that selection for improved wood properties in the innermost rings would also result in improvement of wood properties in the subsequent rings, as well as improved average performance of the entire core. 8. Strong genetic correlations between pure species and hybrid performance for each of the wood quality traits were observed in the hybrid pines. Parental performance can be used to identify the hybrid families which are most likely to have superior juvenile wood properties of the slash/Caribbean F1 hybrid in southeast Queensland. 9. Large unfavourable genetic correlations between growth and wood quality traits were a prominent feature in radiata pine, indicating that overcoming this unfavourable genetic correlation will be a major technical issue in progressing radiata pine breeding. 10. The project created the first radiata pine 18 k cDNA microarray and generated 5,952 radiata pine xylogenesis expressed sequence tags (ESTs) which assembled into 3,304 unigenes. 11. A total of 348 genes were identified as preferentially expressed genes in earlywood or latewood while a total of 168 genes were identified as preferentially expressed genes in either juvenile or mature wood. 12. Juvenile earlywood has a distinct transcriptome relative to other stages of wood development. 13. Discovered rapid decay of linkage disequilibrium (LD) in radiata pine with LD decaying to approximately 50% within 1,700 base pairs (within a typical gene). A total of 913 SNPS from sequencing 177,380 base pairs were identified for association genetic studies. 14. 149 SNPs from 44 genes and 255 SNPs from a further 51 genes (total 95 genes) were selected for association analysis with 62 wood traits, and 30 SNPs were shortlisted for their significant association with variation of wood quality traits (density, MFA and MoE) with individual significant SNPs accounting for between 1.9 and 9.7% of the total genetic variation in traits. 15. Index selection using breeding objectives was the most profitable selection method for radiata pine, but in the long term it may not be the most effective in dealing with negative genetic correlations between wood volume and quality traits. A combination of economic and biological approaches may be needed to deal with the strong adverse correlation.

Gene expression evidence for off-target effects caused by RNA interference-mediated gene silencing of Ubiquitin-63E in the cattle tick Rhipicephalus microplus.

Knowledge of cattle tick (Rhipicephalus (Boophilus) microplus; Acari: Ixodidae) molecular and cellular pathways has been hampered by the lack of an annotated genome. In addition, most of the tick expressed sequence tags (ESTs) available to date consist of similar to 50% unassigned sequences without predicted functions. The most common approach to address this has been the application of RNA interference (RNAi) methods to investigate genes and their pathways. This approach has been widely adopted in tick research despite minimal knowledge of the tick RNAi pathway and double-stranded RNA (dsRNA) uptake mechanisms. A strong knockdown phenotype of adult female ticks had previously been observed using a 594 bp dsRNA targeting the cattle tick homologue for the Drosophila Ubiquitin-63E gene leading to nil or deformed eggs. A NimbleGen cattle tick custom microarray based on the BmiGI.V2 database of R. microplus ESTs was used to evaluate the expression of mRNAs harvested from ticks treated with the tick Ubiquitin-63E 594 bp dsRNA compared with controls. A total of 144 ESTs including TC6372 (Ubiquitin-63E) were down-regulated with 136 ESTs up-regulated following treatment. The results obtained substantiated the knockdown phenotype with ESTs identified as being associated with ubiquitin proteolysis as well as oogenesis, embryogenesis, fatty acid synthesis and stress responses. A bioinformatics analysis was undertaken to predict off-target effects (OTE) resulting from the in silico dicing of the 594 bp Ubiquitin-63E dsRNA which identified 10 down-regulated ESTs (including TC6372) within the list of differentially expressed probes on the microarrays. Subsequent knockdown experiments utilising 196 and 109 bp dsRNAs, and a cocktail of short hairpin RNAs (shRNA) targeting Ubiquitin-63E, demonstrated similar phenotypes for the dsRNAs but nil effect following shRNA treatment. Quantitative reverse transcriptase PCR analysis confirmed differential expression of TC6372 and selected ESTs. Our study demonstrated the minimisation of predicted OTEs in the shorter dsRNA treatments (similar to 100-200 bp) and the usefulness of microarrays to study knockdown phenotypes.

Genome-wide identification and characterization of cytochrome P450 monooxygenase genes in the ciliate Tetrahymena thermophila

Background: Cytochrome P450 monooxygenases play key roles in the metabolism of a wide variety of substrates and they are closely associated with endocellular physiological processes or detoxification metabolism under environmental exposure. To date, however, none has been systematically characterized in the phylum Ciliophora. T. thermophila possess many advantages as a eukaryotic model organism and it exhibits rapid and sensitive responses to xenobiotics, making it an ideal model system to study the evolutionary and functional diversity of the P450 monooxygenase gene family. Results: A total of 44 putative functional cytochrome P450 genes were identified and could be classified into 13 families and 21 sub-families according to standard nomenclature. The characteristics of both the conserved intron-exon organization and scaffold localization of tandem repeats within each P450 family clade suggested that the enlargement of T. thermophila P450 families probably resulted from recent separate small duplication events. Gene expression patterns of all T. thermophila P450s during three important cell physiological stages (vegetative growth, starvation and conjugation) were analyzed based on EST and microarray data, and three main categories of expression patterns were postulated. Evolutionary analysis including codon usage preference, sit-especific selection and gene-expression evolution patterns were investigated and the results indicated remarkable divergences among the T. thermophila P450 genes. Conclusion: The characterization, expression and evolutionary analysis of T. thermophila P450 monooxygenase genes in the current study provides useful information for understanding the characteristics and diversities of the P450 genes in the Ciliophora, and provides the baseline for functional analyses of individual P450 isoforms in this model ciliate species.

Identification of differentially expressed genes from the cross-subfamily cloned embryos derived from zebrafish nuclei and rare minnow enucleated eggs

Cross-species nuclear transfer (NT) has been used to retain the genetic viability of a species near extinction. However, unlike intra-species NT, most embryos produced by cross-species NT were unable to develop to later stages due to incompatible nucleocytoplasmic interactions between the donor nuclei and the recipient cytoplasm from different species. To study the early nucleocytoplasmic interaction in cross-species NT, two laboratory fish species (zebrafish and rare minnow) from different subfamilies were used to generate cross-subfamily NT embryos in the present study. Suppression subtractive hybridization (SSH) was performed to screen out differentially expressed genes from the forward and reverse subtractive cDNA libraries. After dot blot and real-time PCR analysis, 80 of 500 randomly selective sequences were proven to be differentially expressed in the cloned embryos. Among them, 45 sequences shared high homology with 28 zebrafish known genes, and 35 sequences were corresponding to 22 novel expressed sequence tags (ESTs). Based on functional clustering and literature mining analysis, up-and down-regulated genes in the cross-subfamily cloned embryos were mostly relevant to transcription and translation initiation, cell cycle regulation, protein binding, etc. To our knowledge, this is the first report on the determination of genes involved in the early development of cross-species NT embryos of fish. (C) 2007 Elsevier Inc. All rights reserved.