Production of full-length cDNA sequences by sequencing and analysis of expressed sequence tags from Schistosoma mansoni

The number of sequences generated by genome projects has increased exponentially, but gene characterization has not followed at the same rate. Sequencing and analysis of full-length cDNAs is an important step in gene characterization that has been used nowadays by several research groups. In this work, we have selected Schistosoma mansoni clones for full-length sequencing, using an algorithm that investigates the presence of the initial methionine in the parasite sequence based on the positions of alignment start between two sequences. BLAST searches to produce such alignments have been performed using parasite expressed sequence tags produced by Minas Gerais Genome Network against sequences from the database Eukaryotic Cluster of Orthologous Groups (KOG). This procedure has allowed the selection of clones representing 398 proteins which have not been deposited as S. mansoni complete CDS in any public database. Dedicated sequencing of 96 of such clones with reads from both 5' and 3' ends has been performed. These reads have been assembled using PHRAP, resulting in the production of 33 full-length sequences that represent novel S. mansoni proteins. These results shall contribute to construct a more complete view of the biology of this important parasite.

Identification and complete sequencing of novel human transcripts through the use of mouse orthologs and testis cDNA sequences

The correct identification of all human genes, and their derived transcripts, has not yet been achieved, and it remains one of the major aims of the worldwide genomics community. Computational programs suggest the existence of 30,000 to 40,000 human genes. However, definitive gene identification can only be achieved by experimental approaches. We used two distinct methodologies, one based on the alignment of mouse orthologous sequences to the human genome, and another based on the construction of a high-quality human testis cDNA library, in an attempt to identify new human transcripts within the human genome sequence. We generated 47 complete human transcript sequences, comprising 27 unannotated and 20 annotated sequences. Eight of these transcripts are variants of previously known genes. These transcripts were characterized according to size, number of exons, and chromosomal localization, and a search for protein domains was undertaken based on their putative open reading frames. In silico expression analysis suggests that some of these transcripts are expressed at low levels and in a restricted set of tissues.

Alternative splicing enriched cDNA libraries identify breast cancer-associated transcripts

Background: Alternative splicing (AS) is a central mechanism in the generation of genomic complexity and is a major contributor to transcriptome and proteome diversity. Alterations of the splicing process can lead to deregulation of crucial cellular processes and have been associated with a large spectrum of human diseases. Cancer-associated transcripts are potential molecular markers and may contribute to the development of more accurate diagnostic and prognostic methods and also serve as therapeutic targets. Alternative splicing-enriched cDNA libraries have been used to explore the variability generated by alternative splicing. In this study, by combining the use of trapping heteroduplexes and RNA amplification, we developed a powerful approach that enables transcriptome-wide exploration of the AS repertoire for identifying AS variants associated with breast tumor cells modulated by ERBB2 (HER-2/neu) oncogene expression. Results: The human breast cell line (C5.2) and a pool of 5 ERBB2 over-expressing breast tumor samples were used independently for the construction of two AS-enriched libraries. In total, 2,048 partial cDNA sequences were obtained, revealing 214 alternative splicing sequence-enriched tags (ASSETs). A subset with 79 multiple exon ASSETs was compared to public databases and reported 138 different AS events. A high success rate of RT-PCR validation (94.5%) was obtained, and 2 novel AS events were identified. The influence of ERBB2-mediated expression on AS regulation was evaluated by capillary electrophoresis and probe-ligation approaches in two mammary cell lines (Hb4a and C5.2) expressing different levels of ERBB2. The relative expression balance between AS variants from 3 genes was differentially modulated by ERBB2 in this model system. Conclusions: In this study, we presented a method for exploring AS from any RNA source in a transcriptome-wide format, which can be directly easily adapted to next generation sequencers. We identified AS transcripts that were differently modulated by ERBB2-mediated expression and that can be tested as molecular markers for breast cancer. Such a methodology will be useful for completely deciphering the cancer cell transcriptome diversity resulting from AS and for finding more precise molecular markers.

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.

Partial cDNA sequence analysis of myosin Va from rainbow trout (Oncorhynchus mykiss) and its relationship to myosin v isoforms from other vertebrates.

Partial cDNA sequences of myosin V from rainbow trout Oncorhynchus mykiss were analyzed and showed high similarity to MVa from other vertebrates. Phylogenetic analysis has shown that events resulting in the formation of paralogous copies of myosin Va, Vb, and Vc occurred before the divergence of vertebrates into different classes. Expression analysis of myosin Va, Vb, and Vc in different O. mykiss tissues revealed MVa exclusively expressed in hypophysis and brain whereas Vb and Vc were expressed in practically all tissues analyzed. The nucleotide sequence for myosin V was explored in a fish species for the first time and these results represent an important start in understanding the organization, evolution, and expression of myosins in early vertebrates. The data presented here represent contributions to the knowledge of rainbow trout genome. A better understanding of this economically important species could assist in development of improved strains of this fish for aquaculture.

HUNT: launch of a full-length cDNA database from the Helix Research Institute

The Helix Research Institute (HRI) in Japan is releasing 4356 HUman Novel Transcripts and related information in the newly established HUNT database. The institute is a joint research project principally funded by the Japanese Ministry of International Trade and Industry, and the clones were sequenced in the governmental New Energy and Industrial Technology Development Organization (NEDO) Human cDNA Sequencing Project. The HUNT database contains an extensive amount of annotation from advanced analysis and represents an essential bioinformatics contribution towards understanding of the gene function. The HRI human cDNA clones were obtained from full-length enriched cDNA libraries constructed with the oligo-capping method and have resulted in novel full-length cDNA sequences. A large fraction has little similarity to any proteins of known function and to obtain clues about possible function we have developed original analysis procedures. Any putative function deduced here can be validated or refuted by complementary analysis results. The user can also extract information from specific categories like PROSITE patterns, PFAM domains, PSORT localization, transmembrane helices and clones with GENIUS structure assignments. The HUNT database can be accessed at http://www.hri.co.jp/HUNT.

GenEST, a powerful bidirectional link between cDNA sequence data and gene expression profiles generated by cDNA-AFLP

The release of vast quantities of DNA sequence data by large-scale genome and expressed sequence tag (EST) projects underlines the necessity for the development of efficient and inexpensive ways to link sequence databases with temporal and spatial expression profiles. Here we demonstrate the power of linking cDNA sequence data (including EST sequences) with transcript profiles revealed by cDNA-AFLP, a highly reproducible differential display method based on restriction enzyme digests and selective amplification under high stringency conditions. We have developed a computer program (GenEST) that predicts the sizes of virtual transcript-derived fragments (TDFs) of in silico-digested cDNA sequences retrieved from databases. The vast majority of the resulting virtual TDFs could be traced back among the thousands of TDFs displayed on cDNA-AFLP gels. Sequencing of the corresponding bands excised from cDNA-AFLP gels revealed no inconsistencies. As a consequence, cDNA sequence databases can be screened very efficiently to identify genes with relevant expression profiles. The other way round, it is possible to switch from cDNA-AFLP gels to sequences in the databases. Using the restriction enzyme recognition sites, the primer extensions and the estimated TDF size as identifiers, the DNA sequence(s) corresponding to a TDF with an interesting expression pattern can be identified. In this paper we show examples in both directions by analyzing the plant parasitic nematode Globodera rostochiensis. Various novel pathogenicity factors were identified by combining ESTs from the infective stage juveniles with expression profiles of ∼4000 genes in five developmental stages produced by cDNA-AFLP.

Development and evaluation of an automated annotation pipeline and cDNA annotation system

Manual curation has long been held to be the gold standard for functional annotation of DNA sequence. Our experience with the annotation of more than 20,000 full-length cDNA sequences revealed problems with this approach, including inaccurate and inconsistent assignment of gene names, as well as many good assignments that were difficult to reproduce using only computational methods. For the FANTOM2 annotation of more than 60,000 cDNA clones, we developed a number of methods and tools to circumvent some of these problems, including an automated annotation pipeline that provides high-quality preliminary annotation for each sequence by introducing an uninformative filter that eliminates uninformative annotations, controlled vocabularies to accurately reflect both the functional assignments and the evidence supporting them, and a highly refined, Web-based manual annotation tool that allows users to view a wide array of sequence analyses and to assign gene names and putative functions using a consistent nomenclature. The ultimate utility of our approach is reflected in the low rate of reassignment of automated assignments by manual curation. Based on these results, we propose a new standard for large-scale annotation, in which the initial automated annotations are manually investigated and then computational methods are iteratively modified and improved based on the results of manual curation.

Structure and distribution of N-glycans on the S-7-allele stylar self-incompatibility ribonuclease of Nicotiana alata

S-RNases are the stylar products of the self-incompatibility (S)-locus in solanaceous plants (including Nicotiana alata), and as such, are involved in the prevention of self-pollination. All cDNA sequences of S-RNase products of functional S-alleles contain potential N-glycosylation sites, with one site being conserved in all cases, suggesting that N-glycosylation is important in self-incompatibility. In this study, we report on the structure and localization of the N-glycans on the S-7-allele RNase of N, alata, A total of nine N-glycans, belonging to the high-mannose- and xylosylated hybrid-classes, were identified and characterized by a combination of electrospray-ionization mass-spectrometry (ESI-MS), H-1-NMR spectroscopy, and methylation analyses. The glycosylation pattern of individual glycosylation sites was determined by ESI-MS of the glycans released from isolated chymotryptic glycopeptides, All three N-glycosylation sites showed microheterogeneity and each had a unique complement of N-glycans, The N-glycosylation pattern of the S-7-RNase is significantly different to those of the S-1- and S-2-RNases.

Transient expression of a novel serine protease in the ectoderm of the ascidian Herdmania momus during development

We have studied gene expression during ascidian embryonic development using the technique of differential display and isolated partial cDNA sequences of 12 genes. Developmental regulation of these genes has been confirmed by northern hybridization analysis. Further cDNA cloning and sequence analysis of an mRNA that is present during gastrulation, neurulation and tailbud formation reveals that it encodes a novel serine protease containing a single kringle motif and catalytic domain. The spatial expression of this gene, designated Hmserp1, is restricted to precursor cells of the epidermis. The structure and expression of Hmsery1 is discussed in relation to possible functions during development.

Novel gene containing multiple epidermal growth factor-like motifs transiently expressed in the papillae of the ascidian tadpole larvae

We have investigated molecular mechanisms of the embryonic development of an ascidian, a primitive chordate which shares features of both invertebrates and vertebrates, with a view to identifying genes involved in development and metamorphosis, We isolated 12 partial cDNA sequences which were expressed in a stage-specific manner using differential display, We report here the isolation of a full-length cDNA sequence for one of these genes which was specifically expressed during the tailbud and larval stages of ascidian development, This cDNA, 1213 bp in length, is predicted to encode a protein of 337 amino acids containing four epidermal growth factor (EGF)-like repeats and three novel cysteine-rich repeats, Characterization of its spatial expression pattern by in situ hybridisation in late tailbud and larval embryos demonstrated strong expression localised throughout the papillae and anteriormost trunk and weaker expression in the epidermis of the remainder of the embryo, As recent evidence indicates that the signal for metamorphosis originates in the anterior trunk region, these results suggest that this gene may have a role in signalling the initiation of metamorphosis. (C) 1997 Wiley-Liss, Inc.

Schistosoma mansoni encodes SMT3B and SMT3C molecules responsible for post-translational modification of cellular proteins

The sumoylation pathway is a post-translational modification of nuclear proteins widespread among several organisms. SMT3C is the main protein involved in this process and it is covalently conjugated to a diverse assortment of nuclear protein targets. To date, 3 SUMO paralogues (SMT3C, A/B) have been characterized in mammals and plants. In this work we characterized two SUMO related genes, named SMT3B and SMT3C throughout Schistosoma mansoni life cycle. The SmSMTB/C encodes for proteins sharing significant amino acid homology with SMT3. Phylogenetical analyses revealed that both SmSMT3B/C are distinct proteins. Additionally, SmSMT3B and C are expressed in cercariae, adult worms, eggs and schistosomula however SinSMT3C gene showed an expression level 7 to 9 fold higher than SmSMT3B in eggs, schistosomula and adult worms. The comparison between the SmSMT3C genomic and cDNA sequences established that the encoding sequence is interrupted by 3 introns of 70, 37 and 36 bp. Western Blot has shown SMT3 conjugates are present in nuclear and total protein fractions of adults and cercariae. Therefore our results suggest a functional sumoylation pathway, and the presence of two paralogues also suggests the specificity of substrates for SMT3 in S. mansoni. (c) 2008 Elsevier Ireland Ltd. All rights reserved.

Genes that code for T cell signaling proteins establish transcriptional regulatory networks during thymus ontogeny

Gene expression profiling by cDNA microarrays during murine thymus ontogeny has contributed to dissecting the large-scale molecular genetics of T cell maturation. Gene profiling, although useful for characterizing the thymus developmental phases and identifying the differentially expressed genes, does not permit the determination of possible interactions between genes. In order to reconstruct genetic interactions, on RNA level, within thymocyte differentiation, a pair of microarrays containing a total of 1,576 cDNA sequences derived from the IMAGE MTB library was applied on samples of developing thymuses (14-17 days of gestation). The data were analyzed using the GeneNetwork program. Genes that were previously identified as differentially expressed during thymus ontogeny showed their relationships with several other genes. The present method provided the detection of gene nodes coding for proteins implicated in the calcium signaling pathway, such as Prrg2 and Stxbp3, and in protein transport toward the cell membrane, such as Gosr2. The results demonstrate the feasibility of reconstructing networks based on cDNA microarray gene expression determinations, contributing to a clearer understanding of the complex interactions between genes involved in thymus/thymocyte development.

Eucalyptus ESTs involved in mechanisms against plant pathogens and environmental stresses

Eucalypt plantation has high economical importance in Brazil; however, it has been attacked by various pathogens under different environmental stress conditions. Disease resistance and survival under unfavorable environmental conditions have revealed that the eucalypt has developed highly efficient defense systems. Here we show the results of the Eucalyptus ESTs Genome Project (FORESTs). Using the expressed sequence tags (ESTs) obtained by the Project, contigs of similar sequences from each cDNA library induced and not induced by stress agents were formed, and cDNA sequences similar to other already known molecules, such as plant-signaling molecules, phytoalexins, lignin biosynthesis pathways, PR-proteins and putative genes corresponding to enzymes involved in the detoxification of reactive oxygen species, were identified. We also present general considerations about the mechanisms of Eucalyptus defense against biotic and abiotic stresses. These data are of extreme importance for future eucalypt breeding programs aimed at developing plants with enhanced resistance against pathogens and environmental stresses.

A simple RT-PCR-based strategy for screening connexin identity

Vertebrate gap junctions are aggregates of transmembrane channels which are composed of connexin (Cx) proteins encoded by at least fourteen distinct genes in mammals. Since the same Cx type can be expressed in different tissues and more than one Cx type can be expressed by the same cell, the thorough identification of which connexin is in which cell type and how connexin expression changes after experimental manipulation has become quite laborious. Here we describe an efficient, rapid and simple method by which connexin type(s) can be identified in mammalian tissue and cultured cells using endonuclease cleavage of RT-PCR products generated from "multi primers" (sense primer, degenerate oligonucleotide corresponding to a region of the first extracellular domain; antisense primer, degenerate oligonucleotide complementary to the second extracellular domain) that amplify the cytoplasmic loop regions of all known connexins except Cx36. In addition, we provide sequence information on RT-PCR primers used in our laboratory to screen individual connexins and predictions of extension of the "multi primer" method to several human connexins.