Expressed sequence tag analysis of expression profiles of zebrafish testis and ovary

In the present study, two gonad cDNA libraries from zebrafish testes and ovaries were constructed and a total of 1025 expressed sequence tag (EST) clones were generated from the two libraries: 501 from the testis library and 524 from the ovary library. A total of 641 of the EST clones were identified to share significant sequence identity with known sequences in GenBank, representing at least 478 different zebrafish genes. In order to understand the molecular compositions of the two gonad organs, the expression profiles of the identified clones in these two gonad cDNA libraries were analyzed. Both gonad libraries have a higher portion of clones for nuclear proteins and a lower portion for proteins in translational machinery, cytoskeleton and mitochondria than our previously characterized whole-adult cDNA library. Most abundant cDNA clones in the two gonad libraries were identified and over 10% of ovary clones were found to encode egg membrane proteins (zona pellucida or ZP proteins). Furthermore, the testis library showed a more even distribution of cDNA clones with relatively fewer abundant clones that tend to contribute redundant clones in EST projects; thus, the testis library can supply more unique and novel cDNA sequences in a zebrafish EST project. Another aim of this study is to identify cDNA clones that can be used as molecular markers for the analysis of the gonad development in zebrafish. Eleven potential clones were selected to analyze their expression patterns by Northern blot hybridization. Most of them showed a specific or predominant expression in the expected testis or ovary tissue. At last, four of the clones were found, by section in situ hybridization, to be expressed specifically in the germ cells of the testis or ovary and thus they are suitable molecular markers for analyses of spermatogenesis and oogenesis.     

Gene expression profile of rabbit cartilage by expressed sequence tag analysis

Although the rabbit is commonly used as an animal model for the in vivo study of cartilage formation or regeneration, genetic approaches to the rabbit cartilage are rare. We constructed an expressed sequence tag (EST) library from rabbit cartilage tissue for the first time to establish the foundations for genetic study on rabbit cartilage. From our results, we identified 2387 unique genes among 4885 clones, corresponding to 1839 matched to characterized genes including 1618 genes with known function and 548 uncharacterized and novel genes. Gene expression profiles based on EST frequency show that type II collagen (COL2A1) and type X collagen (COL10A1) among collagen clones, proteoglycan 4 (PRG4) and decorin (DCN) among proteoglycan clones, and cartilage oligomeric matrix protein (COMP) and matrix Gla protein (MGP) among other extracellular matrix clones, are highly expressed in rabbit cartilage. In addition, gene expression analysis based on real-time PCR of these major extracellular matrix constituents showed that expression of col2a1 and col10a1 remains constant whereas the expression of prg4, dcn, and comp reveals substantial change with rabbit age. This EST library will provide a valuable resource with which to identify genes involved in the biochemical and physiological functions of rabbit cartilage, and will contribute to establishing the rabbit as an animal model for cartilage research.     

Ciona intestinalis cDNA projects: expressed sequence tag analyses and gene expression profiles during embryogenesis

Ascidians are primitive chordates. Their fertilized egg develops quickly into a tadpole-type larva, which consists of a small number but distinct types of cells, including those of epidermis, central nervous system with two sensory organs, endoderm and mesenchyme in the trunk, and notochord and muscle in the tail. This configuration of the ascidian tadpole is thought to represent the most simplified and primitive chordate body plan. In addition, the free-swimming and non-feeding larvae metamorphose into sessile and filter-feeding adults. The genome size of Ciona intestinalis is estimated to be about 160 Mb, and the number of genes approximately 15,500. The present Ciona cDNA projects focused on gene expression profiles of fertilized eggs, 32-110-cell stage embryos, tailbud embryos, larvae, and young adults. Expressed sequence tags (ESTs) of the 5'-most end and 3'-most end of more than 3000 clones were determined at each developmental stage, and the clones were categorized into independent clusters using the 3'-end sequences. Nearly 1000 clusters of them were then analyzed in detail of their sequences against a BLASTX search. This analysis demonstrates that, on average, half of the clusters showed proteins with sequence similarities to known proteins and the other half did not show sequence similarities to known proteins. Genes with sequence similarities were further categorized into three major subclasses, depending on their functions. Furthermore, the expression profiles of all of the clusters were analyzed by whole-mount in situ hybridization. This analysis highlights gene expression patterns characteristic to each developmental stage. As a result, the present study provides many new molecular markers for each of the tissues and/or organs that constitutes the Ciona tailbud embryo. This sequence information will be used for further comparative genome studies to explore molecular mechanisms involved in the formation of one of the most primitive chordate body plans. All of the data fully characterized may be viewed at the web site http://ghost.zool.kyoto-u.ac.jp.     

Comparative gene expression in the symbiotic and aposymbiotic Aiptasia pulchella by expressed sequence tag analysis

Intracellular symbiotic relationships are prevalent between cnidarians, such as corals and sea anemones, and the photosynthetic dinoflagellate symbionts. However, there is little understanding about how the genes express when the symbiotic relationship is set up. To characterize genes involved in this association, the endosymbiosis between sea anemone, Aiptasia pulchella, and dinoflagellate zooxanthellae, Symbiodinium spp., was employed as a model. Two complementary DNA (cDNA) libraries were constructed from RNA isolated from symbiotic and aposymbiotic A. pulchella. Using single-pass sequencing of cDNA clones, a total of 870 expressed sequence tags (ESTs) clones were generated from the two libraries: 474 from symbiotic animal and 396 from aposymbiotic animal. The initial ESTs consisted of 143 clusters and 231 singletons. A BLASTX search revealed that 147 unique genes had similarities with protein sequences available from databases; 120 of these clones were categorized according to their putative function. However, many ESTs could not assign functionally. The putative roles of some of the identified genes relative to endosymbiosis were discussed. This is the first report of the use of EST analysis to examine the gene expression in symbiotic and aposymbiotic states of the cnidarians. The systematic analysis of EST from this study provides a useful database for future investigations of the molecular mechanisms involved in algal-cnidarian symbiosis.     

Redundancy based detection of sequence polymorphisms in expressed sequence tag data using autoSNP

AutoSNP is a program to detect single nucleotide polymorphisms (SNPs) and insertion/deletion polymorphisms (indels) in expressed sequence tag (EST) data. The program uses d2cluster and cap3 to cluster and align EST sequences, and uses redundancy to differentiate between candidate SNPs and sequence errors. Candidate polymorphisms are identified as occurring in multiple reads within an alignment. For each candidate SNP, two measures of confidence are calculated, the redundancy of the polymorphism at a SNP locus and the co segregation of the candidate SNP with other SNPs in the alignment. AVAILABILITY: The program was written in PERL and is freely available to non-commercial users by request from the authors.     

Identification and validation of Asteraceae miRNAs by the expressed sequence tag analysis

MicroRNAs (miRNAs) are small non-coding RNA molecules that play a vital role in the regulation of gene expression. Despite their identification in hundreds of plant species, few miRNAs have been identified in the Asteraceae, a large family that comprises approximately one tenth of all flowering plants. In this study, we used the expressed sequence tag (EST) analysis to identify potential conserved miRNAs and their putative target genes in the Asteraceae. We applied quantitative Real-Time PCR (qRT-PCR) to confirm the expression of eight potential miRNAs in Carthamus tinctorius and Helianthus annuus. We also performed qRT-PCR analysis to investigate the differential expression pattern of five newly identified miRNAs during five different cotyledon growth stages in safflower. Using these methods, we successfully identified and characterized 151 potentially conserved miRNAs, belonging to 26 miRNA families, in 11 genus of Asteraceae. EST analysis predicted that the newly identified conserved Asteraceae miRNAs target 130 total protein-coding ESTs in sunflower and safflower, as well as 433 additional target genes in other plant species. We experimentally confirmed the existence of seven predicted miRNAs, (miR156, miR159, miR160, miR162, miR166, miR396, and miR398) in safflower and sunflower seedlings. We also observed that five out of eight miRNAs are differentially expressed during cotyledon development. Our results indicate that miRNAs may be involved in the regulation of gene expression during seed germination and the formation of the cotyledons in the Asteraceae. The findings of this study might ultimately help in the understanding of miRNA-mediated gene regulation in important crop species.     

Mining for single nucleotide polymorphisms and insertions/deletions in maize expressed sequence tag data

We have developed a computer based method to identify candidate single nucleotide polymorphisms (SNPs) and small insertions/deletions from expressed sequence tag data. Using a redundancy-based approach, valid SNPs are distinguished from erroneous sequence by their representation multiple times in an alignment of sequence reads. A second measure of validity was also calculated based on the cosegregation of the SNP pattern between multiple SNP loci in an alignment. The utility of this method was demonstrated by applying it to 102,551 maize (Zea mays) expressed sequence tag sequences. A total of 14,832 candidate polymorphisms were identified with an SNP redundancy score of two or greater. Segregation of these SNPs with haplotype indicates that candidate SNPs with high redundancy and cosegregation confidence scores are likely to represent true SNPs. This was confirmed by validation of 264 candidate SNPs from 27 loci, with a range of redundancy and cosegregation scores, in four inbred maize lines. The SNP transition/transversion ratio and insertion/deletion size frequencies correspond to those observed by direct sequencing methods of SNP discovery and suggest that the majority of predicted SNPs and insertion/deletions identified using this approach represent true genetic variation in maize.     

Integration of expressed sequence tag data flanking predicted RNA secondary structures facilitates novel non-coding RNA discovery

Many computational methods have been used to predict novel non-coding RNAs (ncRNAs), but none, to our knowledge, have explicitly investigated the impact of integrating existing cDNA-based Expressed Sequence Tag (EST) data that flank structural RNA predictions. To determine whether flanking EST data can assist in microRNA (miRNA) prediction, we identified genomic sites encoding putative miRNAs by combining functional RNA predictions with flanking ESTs data in a model consistent with miRNAs undergoing cleavage during maturation. In both human and mouse genomes, we observed that the inclusion of flanking ESTs adjacent to and not overlapping predicted miRNAs significantly improved the performance of various methods of miRNA prediction, including direct high-throughput sequencing of small RNA libraries. We analyzed the expression of hundreds of miRNAs predicted to be expressed during myogenic differentiation using a customized microarray and identified several known and predicted myogenic miRNA hairpins. Our results indicate that integrating ESTs flanking structural RNA predictions improves the quality of cleaved miRNA predictions and suggest that this strategy can be used to predict other non-coding RNAs undergoing cleavage during maturation.     

Comparison of quantitative trait loci for adaptive traits between oak and chestnut based on an expressed sequence tag consensus map

A comparative genetic and QTL mapping was performed between Quercus robur L. and Castanea sativa Mill., two major forest tree species belonging to the Fagaceae family. Oak EST-derived markers (STSs) were used to align the 12 linkage groups of the two species. Fifty-one and 45 STSs were mapped in oak and chestnut, respectively. These STSs, added to SSR markers previously mapped in both species, provided a total number of 55 orthologous molecular markers for comparative mapping within the Fagaceae family. Homeologous genomic regions identified between oak and chestnut allowed us to compare QTL positions for three important adaptive traits. Colocation of the QTL controlling the timing of bud burst was significant between the two species. However, conservation of QTL for height growth was not supported by statistical tests. No QTL for carbon isotope discrimination was conserved between the two species. Putative candidate genes for bud burst can be identified on the basis of colocations between EST-derived markers and QTL.     

Genome-scale analysis of human mRNA 5' coding sequences based on expressed sequence tag (EST) database

The "5' end mRNA artifact" issue refers to the incorrect assignment of the first AUG codon in an mRNA, due to the incomplete determination of its 5' end sequence. We performed a systematic identification of coding regions at the 5' end of all human known mRNAs, using an automated expressed sequence tag (EST)-based approach. Following parsing of more than 7 million BLAT alignments, we found 477 human loci, out of 18,665 analyzed, in which an extension of the mRNA 5' coding region was identified. Proof-of-concept confirmation was obtained by in vitro cloning and sequencing for GNB2L1, QARS and TDP2 cDNAs, and the consequences for the functional studies of these loci are discussed. We also generated a list of 20,775 human mRNAs where the presence of an in-frame stop codon upstream of the known start codon indicates completeness of the coding sequence at 5' in the current form.     

Comparison of gene representation between diploid and triploid Paragonimus westermani by expressed sequence tag analyses

Expressed sequence tag (EST) analysis of the diploid and triploid Paragonimus westermani genes was done to have a rapid and informative outlook of the gene-expression profiles of the parasites. Totals of 506 and 505 ESTs were generated from the diploid and triploid P. westermani cDNA libraries. Based on the BLASTx search results of the diploid P. westermani ESTs, 308 (60.9%) matched significantly with formerly identified genes and 198 (39.1%) showed no significant homology in the GenBank database. A similar homology pattern was shown from the triploid EST BLASTx search results with 346 (68.5%) sharing homology with previously identified genes and 159 (31.5%) showing no significant homology. The EST data from both libraries were analyzed and grouped into 9 categories. Comparison of the 2 EST pools revealed high similarities among the categories of the significantly matched genes. Single genes matched repeatedly were also observed in the 2 EST data. Some genes were found that are not yet characterized in P. westermani; these genes were matched by both the diploid and triploid ESTs. Further study of these genes may provide us with more understanding on the parasite's biology and their specific functions in the 2 strains.     

油料作物EST资源的生物信息学分析

利用生物信息学方法,收集整理GenBank数据库中截至2008年5月收录的油料作物油菜、花生、芝麻、大豆、向日葵、蓖麻、亚麻、棕榈等八种油料作物的表达序列标签(EST)序列信息,共获得1,185,911条EST序列,使用Crosmatch、RepeatMask-er、Phrap、CAP3、EMBOSS、Blast、EST-pipeline、ORF finder、Interproscan、blast2go、IdentiCS等软件,基于Linux操作系统,进行了综合及分类分析。共获得289,892条UniEST序列,通过以上对EST序列信息的基因注释信息,筛选出与油脂代谢相关的基因信息,并以此为基础构建了油料作物油脂代谢途径比较结构图。本研究为油料作物油脂代谢相关基因数据库的构建和不同油料作物油脂代谢异同的比较打下基础。     

Assembly and functional expression of murine glutathione reductase cDNA: a sequence missing in expressed sequence tag libraries

Glutathione reductase (GR) is a chemotherapeutic target. Murine GRcDNA, which contains 85% GC in the 38 codons following the start codon, was assembled from the PCR-amplified exon 1 and a downstream cDNA prior to expression in Escherichia coli as a His(6)-tagged protein. Recombinant GR, an FAD-containing homodimer, corresponds in its enzymic and spectral properties to GR isolated from murine Ehrlich ascites tumor cells. Another cDNA, representing GR with a mitochondrial targeting sequence, yielded two distinct enzymically active expression products.     

An expressed sequence tag database of T-cell-enriched activated chicken splenocytes: sequence analysis of 5251 clones

The cDNA and gene sequences of many mammalian cytokines and their receptors are known. However, corresponding information on avian cytokines is limited due to the lack of cross-species activity at the functional level or strong homology at the molecular level. To improve the efficiency of identifying cytokines and novel chicken genes, a directionally cloned cDNA library from T-cell-enriched activated chicken splenocytes was constructed, and the partial sequence of 5251 clones was obtained. Sequence clustering indicates that 2357 (42%) of the clones are present as a single copy, and 2961 are distinct clones, demonstrating the high level of complexity of this library. Comparisons of the sequence data with known DNA sequences in GenBank indicate that approximately 25% of the clones match known chicken genes, 39% have similarity to known genes in other species, and 11% had no match to any sequence in the database. Several previously uncharacterized chicken cytokines and their receptors were present in our library. This collection provides a useful database for cataloging genes expressed in T cells and a valuable resource for future investigations of gene expression in avian immunology. A chicken EST Web site (http://udgenome. ags.udel. edu/chickest/chick.htm) has been created to provide access to the data, and a set of unique sequences has been deposited with GenBank (Accession Nos. AI979741-AI982511). Our new Web site (http://www. chickest.udel.edu) will be active as of March 3, 2000, and will also provide keyword-searching capabilities for BLASTX and BLASTN hits of all our clones.     

Use of a large-scale Triticeae expressed sequence tag resource to reveal gene expression profiles in hexaploid wheat (Triticum aestivum L.).

The US Wheat Genome Project, funded by the National Science Foundation, developed the first large public Triticeae expressed sequence tag (EST) resource. Altogether, 116,272 ESTs were produced, comprising 100,674 5' ESTs and 15 598 3' ESTs. These ESTs were derived from 42 cDNA libraries, which were created from hexaploid bread wheat (Triticum aestivum L.) and its close relatives, including diploid wheat (T. monococcum L. and Aegilops speltoides L.), tetraploid wheat (T. turgidum L.), and rye (Secale cereale L.), using tissues collected from various stages of plant growth and development and under diverse regimes of abiotic and biotic stress treatments. ESTs were assembled into 18,876 contigs and 23,034 singletons, or 41,910 wheat unigenes. Over 90% of the contigs contained fewer than 10 EST members, implying that the ESTs represented a diverse selection of genes and that genes expressed at low and moderate to high levels were well sampled. Statistical methods were used to study the correlation of gene expression patterns, based on the ESTs clustered in the 1536 contigs that contained at least 10 5' EST members and thus representing the most abundant genes expressed in wheat. Analysis further identified genes in wheat that were significantly upregulated (p < 0.05) in tissues under various abiotic stresses when compared with control tissues. Though the function annotation cannot be assigned for many of these genes, it is likely that they play a role associated with the stress response. This study predicted the possible functionality for 4% of total wheat unigenes, which leaves the remaining 96% with their functional roles and expression patterns largely unknown. Nonetheless, the EST data generated in this project provide a diverse and rich source for gene discovery in wheat.     

Analysis of expressed sequence tag loci on wheat chromosome group 4

A total of 1918 loci, detected by the hybridization of 938 expressed sequence tag unigenes (ESTs) from 26 Triticeae cDNA libraries, were mapped to wheat (Triticum aestivum L.) homoeologous group 4 chromosomes using a set of deletion, ditelosomic, and nulli-tetrasomic lines. The 1918 EST loci were not distributed uniformly among the three group 4 chromosomes; 41, 28, and 31% mapped to chromosomes 4A, 4B, and 4D, respectively. This pattern is in contrast to the cumulative results of EST mapping in all homoeologous groups, as reported elsewhere, that found the highest proportion of loci mapped to the B genome. Sixty-five percent of these 1918 loci mapped to the long arms of homoeologous group 4 chromosomes, while 35% mapped to the short arms. The distal regions of chromosome arms showed higher numbers of loci than the proximal regions, with the exception of 4DL. This study confirmed the complex structure of chromosome 4A that contains two reciprocal translocations and two inversions, previously identified. An additional inversion in the centromeric region of 4A was revealed. A consensus map for homoeologous group 4 was developed from 119 ESTs unique to group 4. Forty-nine percent of these ESTs were found to be homoeologous to sequences on rice chromosome 3, 12% had matches with sequences on other rice chromosomes, and 39% had no matches with rice sequences at all. Limited homology (only 26 of the 119 consensus ESTs) was found between wheat ESTs on homoeologous group 4 and the Arabidopsis genome. Forty-two percent of the homoeologous group 4 ESTs could be classified into functional categories on the basis of blastX searches against all protein databases.