The contribution of 700,000 ORF sequence tags to the definition of the human transcriptome

In vitro and in silico validation of CA3 and FHL1 downregulation in oral cancer

Background

Aberrant methylation is a frequent event in oral cancer.

Methods

In order to better characterize these alterations, a search for genes downregulated by aberrant methylation in oral squamous cell carcinoma (OSCC) was conducted through the mining of ORESTES dataset. Findings were further validated in OSCC cell lines and patients’ samples and confirmed using TCGA data. Differentially expressed genes were identified in ORESTES libraries and validated in vitro using RT-PCR in HNSCC cell-lines and OSCC tumor samples. Further confirmation of these results was performed using mRNA expression and methylation data from The Cancer Genome Atlas (TCGA) data.

Results

From the set of genes selected for validation, CA3 and FHL1 were downregulated in 60% (12/20) and 75% (15/20) of OSCC samples, respectively, and in HNSCC cell lines. The treatment of cell lines JHU-13 and FaDu with the demethylating agent 5'-aza-dC was efficient in restoring CA3 and FHL1 expression. TCGA expression and methylation data on OSCC confirms the downregulation of these genes in OSCC samples and also suggests that expression of CA3 and FHL1 is probably regulated by methylation. The downregulation of CA3 and FHL1 observed in silico was validated in HNSCC cell lines and OSCC samples, showing the feasibility of integrating different datasets to select differentially expressed genes in silico.

Conclusions

These results showed that the downregulation of CA3 and FHL1 data observed in the ORESTES libraries was validated in HNSCC cell lines and OSCC samples and in a large cohort of samples from the TCGA database. Moreover, it suggests that expression of CA3 and FHL1 could probably be regulated by methylation having an important role the oral carcinogenesis.

Electronic supplementary material

The online version of this article (10.1186/s12885-018-4077-3) contains supplementary material, which is available to authorized users.

Analysis of a novel protein in human colorectal adenocarcinoma

Colorectal adenocarcinoma (CRC) is the third most common type of cancer worldwide with a low 5‑year survival rate. The present study aimed to investigate the structure and function of a novel protein identified from human colorectal adenocarcinoma (CRC). A differentially expressed sequence tag (GenBank accession number, ES274081) was collected from GenBank. Bioinformatics tools were employed to obtain the sequence of the full‑length cDNA in order to localize the open reading frame and to predict the protein sequence. Mass spectro-metry was used to analyze the structure of this novel protein and western blot analysis was used to confirm the expression of this protein in human CRC tissue samples. The full‑length cDNA was composed of 4,283‑bp nucleotides and the sequence information was obtained (GenBank accession number, NM_001013649). The corresponding protein molecule contained 165 amino acids, with a monoisotopic molecular weight of 18.6033 kDa and an isoelectric point of 8.43, determined by mass spectrometry. The protein structure and its function in adenocarcinoma were further explored. In the present study, a novel protein, which may be involved in nuclear signal transduction, was identified using bioinformatics, mass spectrometry and western blot analysis.

Identification of upregulated genes in oral squamous cell carcinomas

BACKGROUND

Oral cancer is the most common subset of head and neck squamous cell carcinomas (HNSCC). These tumors often have an aggressive clinical outcome hallmarked by a propensity for local invasion and regional nodal metastasis. Upregulated genes could be useful as markers for diagnosis, prognosis, and as new drug targets for these tumors.

METHODS

To identify upregulated genes in oral squamous cell carcinomas (OSSCs), we examined the ORESTES public database and used a quantitative reverse transcription-polymerase chain reaction (qRT-PCR) approach to determine the expression level of selected genes in tumor samples.

RESULTS AND CONCLUSIONS

The ORESTES data mining analysis indicated 40 upregulated genes in HNSCC. Nine of these candidate genes were selected for further qRT-PCR validation and 3 of them (ALDOA, AHSA1, and POLQ) were frequently found upregulated in OSCC samples, which may indicate an association of these genes with the carcinogenesis process in this tumor site and they can constitute potential new targets for therapy.

Expression of human protein S100A7 (psoriasin), preparation of antibody and application to human larynx squamous cell carcinoma

BACKGROUND

Up-regulation of S100A7 (Psoriasin), a small calcium-binding protein, is associated with the development of several types of carcinomas, but its function and possibility to serve as a diagnostic or prognostic marker have not been fully defined. In order to prepare antibodies to the protein for immunohistochemical studies we produced the recombinant S100A7 protein in E. coli. mRNA extracted from human tracheal tumor tissue which was amplified by RT-PCR to provide the region coding for the S100A7 gene. The amplified fragment was cloned in the vector pCR2.1-TOPO and sub-cloned in the expression vector pAE. The protein rS100A7 (His-tag) was expressed in E. coli BL21::DE3, purified by affinity chromatography on an Ni-NTA column, recovered in the 2.0 to 3.5 mg/mL range in culture medium, and used to produce a rabbit polyclonal antibody anti-rS100A7 protein. The profile of this polyclonal antibody was evaluated in a tissue microarray.

RESULTS

The rS100A7 (His-tag) protein was homogeneous by SDS-PAGE and mass spectrometry and was used to produce an anti-recombinant S100A7 (His-tag) rabbit serum (polyclonal antibody anti-rS100A7). The molecular weight of rS100A7 (His-tag) protein determined by linear MALDI-TOF-MS was 12,655.91 Da. The theoretical mass calculated for the nonapeptide attached to the amino terminus is 12,653.26 Da (delta 2.65 Da). Immunostaining with the polyclonal anti-rS100A7 protein generated showed reactivity with little or no background staining in head and neck squamous cell carcinoma cells, detecting S100A7 both in nucleus and cytoplasm. Lower levels of S100A7 were detected in non-neoplastic tissue.

CONCLUSIONS

The polyclonal anti-rS100A7 antibody generated here yielded a good signal-to-noise contrast and should be useful for immunohistochemical detection of S100A7 protein. Its potential use for other epithelial lesions besides human larynx squamous cell carcinoma and non-neoplastic larynx should be explored in future.

Long noncoding intronic RNAs are differentially expressed in primary and metastatic pancreatic cancer

Background

Pancreatic ductal adenocarcinoma (PDAC) is known by its aggressiveness and lack of effective therapeutic options. Thus, improvement in current knowledge of molecular changes associated with pancreatic cancer is urgently needed to explore novel venues of diagnostics and treatment of this dismal disease. While there is mounting evidence that long noncoding RNAs (lncRNAs) transcribed from intronic and intergenic regions of the human genome may play different roles in the regulation of gene expression in normal and cancer cells, their expression pattern and biological relevance in pancreatic cancer is currently unknown. In the present work we investigated the relative abundance of a collection of lncRNAs in patients' pancreatic tissue samples aiming at identifying gene expression profiles correlated to pancreatic cancer and metastasis.

Methods

Custom 3,355-element spotted cDNA microarray interrogating protein-coding genes and putative lncRNA were used to obtain expression profiles from 38 clinical samples of tumor and non-tumor pancreatic tissues. Bioinformatics analyses were performed to characterize structure and conservation of lncRNAs expressed in pancreatic tissues, as well as to identify expression signatures correlated to tissue histology. Strand-specific reverse transcription followed by PCR and qRT-PCR were employed to determine strandedness of lncRNAs and to validate microarray results, respectively.

Results

We show that subsets of intronic/intergenic lncRNAs are expressed across tumor and non-tumor pancreatic tissue samples. Enrichment of promoter-associated chromatin marks and over-representation of conserved DNA elements and stable secondary structure predictions suggest that these transcripts are generated from independent transcriptional units and that at least a fraction is under evolutionary selection, and thus potentially functional.

Statistically significant expression signatures comprising protein-coding mRNAs and lncRNAs that correlate to PDAC or to pancreatic cancer metastasis were identified. Interestingly, loci harboring intronic lncRNAs differentially expressed in PDAC metastases were enriched in genes associated to the MAPK pathway. Orientation-specific RT-PCR documented that intronic transcripts are expressed in sense, antisense or both orientations relative to protein-coding mRNAs. Differential expression of a subset of intronic lncRNAs (PPP3CB, MAP3K14 and DAPK1 loci) in metastatic samples was confirmed by Real-Time PCR.

Conclusion

Our findings reveal sets of intronic lncRNAs expressed in pancreatic tissues whose abundance is correlated to PDAC or metastasis, thus pointing to the potential relevance of this class of transcripts in biological processes related to malignant transformation and metastasis in pancreatic cancer.

Temporal blastemal cell gene expression analysis in the kidney reveals new Wnt and related signaling pathway genes to be essential for Wilms' tumor onset

Wilms' tumors (WTs) originate from metanephric blastema cells that are unable to complete differentiation, resulting in triphasic tumors composed of epithelial, stromal and blastemal cells, with the latter harboring molecular characteristics similar to those of the earliest kidney development stages. Precise regulation of Wnt and related signaling pathways has been shown to be crucial for correct kidney differentiation. In this study, the gene expression profile of Wnt and related pathways was assessed in laser-microdissected blastemal cells in WTs and differentiated kidneys, in human and in four temporal kidney differentiation stages (i.e. E15.5, E17.5, P1.5 and P7.5) in mice, using an orthologous cDNA microarray platform. A signaling pathway-based gene signature was shared between cells of WT and of earliest kidney differentiation stages, revealing genes involved in the interruption of blastemal cell differentiation in WT. Reverse transcription-quantitative PCR showed high robustness of the microarray data demonstrating 75 and 56% agreement in the initial and independent sample sets, respectively. The protein expression of CRABP2, IGF2, GRK7, TESK1, HDGF, WNT5B, FZD2 and TIMP3 was characterized in WTs and in a panel of human fetal kidneys displaying remarkable aspects of differentiation, which was recapitulated in the tumor. Taken together, this study reveals new genes candidate for triggering WT onset and for therapeutic treatment targets.

Gene network analyses point to the importance of human tissue kallikreins in melanoma progression

Background

A wide variety of high-throughput microarray platforms have been used to identify molecular targets associated with biological and clinical tumor phenotypes by comparing samples representing distinct pathological states.

Methods

The gene expression profiles of human cutaneous melanomas were determined by cDNA microarray analysis. Next, a robust analysis to determine functional classifications and make predictions based on data-oriented hypotheses was performed. Relevant networks that may be implicated in melanoma progression were also considered.

Results

In this study we aimed to analyze coordinated gene expression changes to find molecular pathways involved in melanoma progression. To achieve this goal, ontologically-linked modules with coordinated expression changes in melanoma samples were identified. With this approach, we detected several gene networks related to different modules that were induced or repressed during melanoma progression. Among them we observed high coordinated expression levels of genes involved in a) cell communication (KRT4, VWF and COMP); b) epidermal development (KLK7, LAMA3 and EVPL); and c) functionally related to kallikreins (EVPL, KLK6, KLK7, KLK8, SERPINB13, SERPING1 and SLPI). Our data also indicated that hKLK7 protein expression was significantly associated with good prognosis and survival.

Conclusions

Our findings, derived from a different type of analysis of microarray data, highlight the importance of analyzing coordinated gene expression to find molecular pathways involved in melanoma progression.

Functional microarray analysis suggests repressed cell-cell signaling and cell survival-related modules inhibit progression of head and neck squamous cell carcinoma

Background

Cancer shows a great diversity in its clinical behavior which cannot be easily predicted using the currently available clinical or pathological markers. The identification of pathways associated with lymph node metastasis (N+) and recurrent head and neck squamous cell carcinoma (HNSCC) may increase our understanding of the complex biology of this disease.

Methods

Tumor samples were obtained from untreated HNSCC patients undergoing surgery. Patients were classified according to pathologic lymph node status (positive or negative) or tumor recurrence (recurrent or non-recurrent tumor) after treatment (surgery with neck dissection followed by radiotherapy). Using microarray gene expression, we screened tumor samples according to modules comprised by genes in the same pathway or functional category.

Results

The most frequent alterations were the repression of modules in negative lymph node (N0) and in non-recurrent tumors rather than induction of modules in N+ or in recurrent tumors. N0 tumors showed repression of modules that contain cell survival genes and in non-recurrent tumors cell-cell signaling and extracellular region modules were repressed.

Conclusions

The repression of modules that contain cell survival genes in N0 tumors reinforces the important role that apoptosis plays in the regulation of metastasis. In addition, because tumor samples used here were not microdissected, tumor gene expression data are represented together with the stroma, which may reveal signaling between the microenvironment and tumor cells. For instance, in non-recurrent tumors, extracellular region module was repressed, indicating that the stroma and tumor cells may have fewer interactions, which disable metastasis development. Finally, the genes highlighted in our analysis can be implicated in more than one pathway or characteristic, suggesting that therapeutic approaches to prevent tumor progression should target more than one gene or pathway, specially apoptosis and interactions between tumor cells and the stroma.

Upstream - news in genomics

Since our last issue, several important genomes have been completely or ‘almost completely’ sequenced. The debate over the number of human genes has flared up once more, with one computational and one experimental study into the annotation of the human genome. The mouse genome project has a clone fingerprint map to aid their sequencing effort. The SAGE technique has been applied to Drosophila and the US National Science Foundation announced increased spending on plant genome research.

Transcriptome analysis of Taenia solium cysticerci using Open Reading Frame ESTs (ORESTES)

BACKGROUND

Human infection by the pork tapeworm Taenia solium affects more than 50 million people worldwide, particularly in underdeveloped and developing countries. Cysticercosis which arises from larval encystation can be life threatening and difficult to treat. Here, we investigate for the first time the transcriptome of the clinically relevant cysticerci larval form.

RESULTS

Using Expressed Sequence Tags (ESTs) produced by the ORESTES method, a total of 1,520 high quality ESTs were generated from 20 ORESTES cDNA mini-libraries and its analysis revealed fragments of genes with promising applications including 51 ESTs matching antigens previously described in other species, as well as 113 sequences representing proteins with potential extracellular localization, with obvious applications for immune-diagnosis or vaccine development.

CONCLUSION

The set of sequences described here will contribute to deciphering the expression profile of this important parasite and will be informative for the genome assembly and annotation, as well as for studies of intra- and inter-specific sequence variability. Genes of interest for developing new diagnostic and therapeutic tools are described and discussed.

The DING family of proteins: ubiquitous in eukaryotes, but where are the genes?

PstS and DING proteins are members of a superfamily of secreted, high-affinity phosphate-binding proteins. Whereas microbial PstS have a well-defined role in phosphate ABC transporters, the physiological function of DING proteins, named after their DINGGG N termini, still needs to be determined. PstS and DING proteins co-exist in some Pseudomonas strains, to which they confer a highly adhesive and virulent phenotype. More than 30 DING proteins have now been purified, mostly from eukaryotes. They are often associated with infections or with dysregulation of cell proliferation. Consequently, eukaryotic DING proteins could also be involved in cell-cell communication or adherence. The ubiquitous presence in eukaryotes of proteins structurally and functionally related to bacterial virulence factors is intriguing, as is the absence of eukaryotic genes encoding DING proteins in databases. DING proteins in eukaryotes could originate from unidentified commensal or symbiotic bacteria and could contribute to essential functions. Alternatively, DING proteins could be encoded by eukaryotic genes sharing special features that prevent their cloning. Both hypotheses are discussed.

A global meta-analysis of microarray expression data to predict unknown gene functions and estimate the literature-data divide

MOTIVATION

Approximately 9334 (37%) of human genes have no publications documenting their function and, for those that are published, the number of publications per gene is highly skewed. Furthermore, for reasons not clear, the entry of new gene names into the literature has slowed in recent years. If we are to better understand human/mammalian biology and complete the catalog of human gene function, it is important to finish predicting putative functions for these genes based upon existing experimental evidence.

RESULTS

A global meta-analysis (GMA) of all publicly available GEO two-channel human microarray datasets (3551 experiments total) was conducted to identify genes with recurrent, reproducible patterns of co-regulation across different conditions. Patterns of co-expression were divided into parallel (i.e. genes are up and down-regulated together) and anti-parallel. Several ranking methods to predict a gene's function based on its top 20 co-expressed gene pairs were compared. In the best method, 34% of predicted Gene Ontology (GO) categories matched exactly with the known GO categories for approximately 5000 genes analyzed versus only 3% for random gene sets. Only 2.4% of co-expressed gene pairs were found as co-occurring gene pairs in MEDLINE.

CONCLUSIONS

Via a GO enrichment analysis, genes co-expressed in parallel with the query gene were frequently associated with the same GO categories, whereas anti-parallel genes were not. Combining parallel and anti-parallel genes for analysis resulted in fewer significant GO categories, suggesting they are best analyzed separately. Expression databases contain much unexpected genetic knowledge that has not yet been reported in the literature. A total of 1642 Human genes with unknown function were differentially expressed in at least 30 experiments.

AVAILABILITY

Data matrix available upon request.

Transcriptome analysis of the salivary glands of the female tick Amblyomma americanum (Acari: Ixodidae)

Ticks infest a wide range of hosts while bypassing their immune, inflammatory and haemostatic responses during their extended feeding, which may last for more than two weeks. Here, we present a transcriptome analysis of 3868 expressed sequence tags (ESTs) from three cDNA libraries generated from the salivary glands of adult female Ambyomma americanum ticks at different stages of feeding. We applied a normalization step for one library, significantly decreasing the abundance of mitochondrial sequences amongst the 2292 sequences from the normalized library. Our ESTs include homologues that may modulate haemostatic, immune and inflammatory responses of the hosts. Other ESTs probably represent important components of the highly efficient secretory pathways for salivary proteins and concomitantly transmitted pathogens.

No-match ORESTES explored as tumor markers

Sequencing technologies and new bioinformatics tools have led to the complete sequencing of various genomes. However, information regarding the human transcriptome and its annotation is yet to be completed. The Human Cancer Genome Project, using ORESTES (open reading frame EST sequences) methodology, contributed to this objective by generating data from about 1.2 million expressed sequence tags. Approximately 30% of these sequences did not align to ESTs in the public databases and were considered no-match ORESTES. On the basis that a set of these ESTs could represent new transcripts, we constructed a cDNA microarray. This platform was used to hybridize against 12 different normal or tumor tissues. We identified 3421 transcribed regions not associated with annotated transcripts, representing 83.3% of the platform. The total number of differentially expressed sequences was 1007. Also, 28% of analyzed sequences could represent noncoding RNAs. Our data reinforces the knowledge of the human genome being pervasively transcribed, and point out molecular marker candidates for different cancers. To reinforce our data, we confirmed, by real-time PCR, the differential expression of three out of eight potentially tumor markers in prostate tissues. Lists of 1007 differentially expressed sequences, and the 291 potentially noncoding tumor markers were provided.

Transcriptome-guided characterization of genomic rearrangements in a breast cancer cell line

We have identified new genomic alterations in the breast cancer cell line HCC1954, using high-throughput transcriptome sequencing. With 120 Mb of cDNA sequences, we were able to identify genomic rearrangement events leading to fusions or truncations of genes including MRE11 and NSD1, genes already implicated in oncogenesis, and 7 rearrangements involving other additional genes. This approach demonstrates that high-throughput transcriptome sequencing is an effective strategy for the characterization of genomic rearrangements in cancers.

Exploiting ESTs in human health

Expressed Sequence Tags (ESTs) are fragments of cDNA clones. They correspond to the most abundant type of cDNA information available in the public databases. ESTs have been used for expression profiling, gene identification, characterization of differentially expressed genes, and identification of transcript variants among other utilities. In this review I will discuss the major features of the collection of ESTs available in the public domain giving a special emphasis on how this dataset has been used in studies about human diseases.

Identification of protein-coding and intronic noncoding RNAs down-regulated in clear cell renal carcinoma

The clear cell subtype of renal cell carcinoma (RCC) is the most lethal and prevalent cancer of the urinary system. To investigate the molecular changes associated with malignant transformation in clear cell RCC, the gene expression profiles of matched samples of tumor and adjacent non-neoplastic tissue were obtained from six patients. A custom-built cDNA microarray platform was used, comprising 2292 probes that map to exons of genes and 822 probes for noncoding RNAs mapping to intronic regions. Intronic transcription was detected in all normal and neoplastic renal tissues. A subset of 55 transcripts was significantly down-regulated in clear cell RCC relative to the matched nontumor tissue as determined by a combination of two statistical tests and leave-one-out patient cross-validation. Among the down-regulated transcripts, 49 mapped to untranslated or coding exons and 6 were intronic relative to known exons of protein-coding genes. Lower levels of expression of SIN3B, TRIP3, SYNJ2BP and NDE1 (P < 0.02), and of intronic transcripts derived from SND1 and ACTN4 loci (P < 0.05), were confirmed in clear cell RCC by Real-time RT-PCR. A subset of 25 transcripts was deregulated in additional six nonclear cell RCC samples, pointing to common transcriptional alterations in RCC irrespective of the histological subtype or differentiation state of the tumor. Our results indicate a novel set of tumor suppressor gene candidates, including noncoding intronic RNAs, which may play a significant role in malignant transformations of normal renal cells.

Complete genome of Phenylobacterium zucineum--a novel facultative intracellular bacterium isolated from human erythroleukemia cell line K562

Background

Phenylobacterium zucineum is a recently identified facultative intracellular species isolated from the human leukemia cell line K562. Unlike the known intracellular pathogens, P. zucineum maintains a stable association with its host cell without affecting the growth and morphology of the latter.

Results

Here, we report the whole genome sequence of the type strain HLK1^T. The genome consists of a circular chromosome (3,996,255 bp) and a circular plasmid (382,976 bp). It encodes 3,861 putative proteins, 42 tRNAs, and a 16S-23S-5S rRNA operon. Comparative genomic analysis revealed that it is phylogenetically closest to Caulobacter crescentus, a model species for cell cycle research. Notably, P. zucineum has a gene that is strikingly similar, both structurally and functionally, to the cell cycle master regulator CtrA of C. crescentus, and most of the genes directly regulated by CtrA in the latter have orthologs in the former.

Conclusion

This work presents the first complete bacterial genome in the genus Phenylobacterium. Comparative genomic analysis indicated that the CtrA regulon is well conserved between C. crescentus and P. zucineum.

Comparative ORESTES-sampling of transcriptomes of immune-challenged Biomphalaria glabrata snails

The snail Biomphalaria glabrata (Gastropoda, Mollusca) is an important intermediate host for the human parasite Schistosoma mansoni (Digenea, Trematoda). Anti-pathogen responses of B. glabrata were studied towards a better understanding of snail immunity and host-parasite compatibility. Open reading frame ESTs (ORESTES) were sampled from different transcriptomes of M line strain B. glabrata, 12h post-challenge with Escherichia coli (Gram-negative), Micrococcus luteus (Gram-positive) bacteria or compatible S. mansoni, and controls. The resulting 3123 ORESTES represented 2129 unique sequences (373 clusters, 1756 singletons). Of these, 175 (8.1%) were putative defense factors, including lectins, antimicrobial peptides and components of various immune-effector systems. Comparison of biological processes (GO-terms) within different transcriptomes indicated that B. glabrata increased oxygen transport and metal binding in reaction to all challenges. Comprehensive comparisons of transcriptomes revealed that responses of B. glabrata against bacteria were similar to each other and differed from the ineffective response to S. mansoni. Furthermore, the response to S. mansoni infection was less comprehensive than that to bacteria. Many novel (unknown) sequences were recovered in association with particular challenges. B. glabrata possesses multi-faceted, potent immune defenses. This agrees with the notion that S. mansoni is capable of immune-evasion and prevents effective host defense responses in order to survive in B. glabrata. Future analysis of the numerous unknown sequences recovered from challenged snails may reveal novel immune factors and provide increased understanding of immunity of B. glabrata in relation to parasite-host compatibility.

Porcine transcriptome analysis based on 97 non-normalized cDNA libraries and assembly of 1,021,891 expressed sequence tags

A resource consisting of one million porcine ESTs is described, providing an essential resource for annotation, comparative genomics, assembly of the pig genome sequence, and further porcine transcription studies.

Background

Knowledge of the structure of gene expression is essential for mammalian transcriptomics research. We analyzed a collection of more than one million porcine expressed sequence tags (ESTs), of which two-thirds were generated in the Sino-Danish Pig Genome Project and one-third are from public databases. The Sino-Danish ESTs were generated from one normalized and 97 non-normalized cDNA libraries representing 35 different tissues and three developmental stages.

Results

Using the Distiller package, the ESTs were assembled to roughly 48,000 contigs and 73,000 singletons, of which approximately 25% have a high confidence match to UniProt. Approximately 6,000 new porcine gene clusters were identified. Expression analysis based on the non-normalized libraries resulted in the following findings. The distribution of cluster sizes is scaling invariant. Brain and testes are among the tissues with the greatest number of different expressed genes, whereas tissues with more specialized function, such as developing liver, have fewer expressed genes. There are at least 65 high confidence housekeeping gene candidates and 876 cDNA library-specific gene candidates. We identified differential expression of genes between different tissues, in particular brain/spinal cord, and found patterns of correlation between genes that share expression in pairs of libraries. Finally, there was remarkable agreement in expression between specialized tissues according to Gene Ontology categories.

Conclusion

This EST collection, the largest to date in pig, represents an essential resource for annotation, comparative genomics, assembly of the pig genome sequence, and further porcine transcription studies.

Identification of unannotated exons of low abundance transcripts in Drosophila melanogaster and cloning of a new serine protease gene upregulated upon injury

BACKGROUND

The sequencing of the D.melanogaster genome revealed an unexpected small number of genes (~ 14,000) indicating that mechanisms acting on generation of transcript diversity must have played a major role in the evolution of complex metazoans. Among the most extensively used mechanisms that accounts for this diversity is alternative splicing. It is estimated that over 40% of Drosophila protein-coding genes contain one or more alternative exons. A recent transcription map of the Drosophila embryogenesis indicates that 30% of the transcribed regions are unannotated, and that 1/3 of this is estimated as missed or alternative exons of previously characterized protein-coding genes. Therefore, the identification of the variety of expressed transcripts depends on experimental data for its final validation and is continuously being performed using different approaches. We applied the Open Reading Frame Expressed Sequence Tags (ORESTES) methodology, which is capable of generating cDNA data from the central portion of rare transcripts, in order to investigate the presence of hitherto unnanotated regions of Drosophila transcriptome.

RESULTS

Bioinformatic analysis of 1,303 Drosophila ORESTES clusters identified 68 sequences derived from unannotated regions in the current Drosophila genome version (4.3). Of these, a set of 38 was analysed by polyA+ northern blot hybridization, validating 17 (50%) new exons of low abundance transcripts. For one of these ESTs, we obtained the cDNA encompassing the complete coding sequence of a new serine protease, named SP212. The SP212 gene is part of a serine protease gene cluster located in the chromosome region 88A12-B1. This cluster includes the predicted genes CG9631, CG9649 and CG31326, which were previously identified as up-regulated after immune challenges in genomic-scale microarray analysis. In agreement with the proposal that this locus is co-regulated in response to microorganisms infection, we show here that SP212 is also up-regulated upon injury.

CONCLUSION

Using the ORESTES methodology we identified 17 novel exons from low abundance Drosophila transcripts, and through a PCR approach the complete CDS of one of these transcripts was defined. Our results show that the computational identification and manual inspection are not sufficient to annotate a genome in the absence of experimentally derived data.

Real-time definition of non-randomness in the distribution of genomic events

Features such as mutations or structural characteristics can be non-randomly or non-uniformly distributed within a genome. So far, computer simulations were required for statistical inferences on the distribution of sequence motifs. Here, we show that these analyses are possible using an analytical, mathematical approach. For the assessment of non-randomness, our calculations only require information including genome size, number of (sampled) sequence motifs and distance parameters. We have developed computer programs evaluating our analytical formulas for the real-time determination of expected values and p-values. This approach permits a flexible cluster definition that can be applied to most effectively identify non-random or non-uniform sequence motif distribution. As an example, we show the effectivity and reliability of our mathematical approach in clinical retroviral vector integration site distribution.

Hippocampal gene expression analysis using the ORESTES methodology shows that homer 1a mRNA is upregulated in the acute period of the pilocarpine epilepsy model

In the study of temporal lobe epilepsy (TLE) the characterization of genes expressed in the hippocampus is of central importance for understanding their roles in epileptogenic mechanisms. Although several large-scale studies on TLE gene expression have been reported, precise assignment of individual genes associated with this syndrome is still debatable. Here we investigated differentially expressed genes by comparison of mRNAs from normal and epileptic rat hippocampus in the pilocarpine model of epilepsy. For this we used a powerful EST sequencing methodology, ORESTES (Open Reading frame Expressed Sequence Tags), which generates sequence datasets enriched for mRNAs open reading frames (ORFs) rather than simple 5' and 3' ends of mRNAs. Analysis of our sequences shows that ORESTES readily enables the identification of epilepsy associated ORFs. PFAM analysis of protein motifs present in our ORESTES epilepsy database revealed diverse important protein family domains, such as cytoskeletal, cell signaling and protein kinase domains, which could be involved in processes underlying epileptogenesis. More importantly, we show that the expression of homer 1a, known to be coupled to mGluR and NMDA synaptic transmission, is associated with pilocarpine induced status epilepticus (SE). The combined use of the pilocarpine model of epilepsy with the ORESTES technique can significantly contribute to the identification of specific genes and proteins related to TLE. This is the first study applying a large-scale method for rapid shotgun sequencing directed to ORFs in epilepsy research.

Biomphalaria glabrata transcriptome: identification of cell-signalling, transcriptional control and immune-related genes from open reading frame expressed sequence tags (ORESTES)

Biomphalaria glabrata is the major intermediate snail host for Schistosoma mansoni, one of the important schistosomes infecting man. Much remains to be discovered concerning specific molecules mediating the defence events in these intermediate hosts, triggered by invading schistosomes. An expressed sequence tag (EST) gene discovery strategy known as ORESTES has been employed to identify transcripts that might be involved in snail-schistosome interactions in order to examine gene expression patterns in infected B. glabrata. Over 3930 ESTs were sequenced from cDNA libraries made from both schistosome-exposed and unexposed snails using different tissue types, producing a database of 1843 non-redundant clones. The non-redundant set has been assessed for gene ontology and KEGG pathway assignments. This approach has revealed a number of signalling, antioxidant and immune-related gene homologues that, based on current understanding of molluscan and other comparative systems, might play an important role in the molluscan defence response towards infection.

Pan-genome isolation of low abundance transcripts using SAGE tag

The SAGE (serial analysis of gene expression) method is sensitive at detecting the lower abundance transcripts. More than a third of human SAGE tags identified are novel representing the low abundance unknown transcripts. Using the GLGI method (generation of longer 3' EST from SAGE tag for gene identification), we converted 1009 low-copy, human X chromosome-specific SAGE tags into 10210 3' ESTs. We identified 3418 unique 3' ESTs, 46% of which are novel and originated from the lower abundance transcripts. However, nearly all 3' ESTs were mapped to various regions across the genome but not X chromosome. Detailed analysis indicates that those 3' ESTs were isolated by SAGE tag mis-priming to the non-parent transcripts. Replacing SAGE tags with non-transcribed genomic DNA tags resulted in poor amplification, indicating that the sequence similarity between different transcripts contributed to the amplification. Our study shows the prevalence of novel low abundance transcripts that can be isolated efficiently through SAGE tags mis-priming.

Auto-antibodies in prostate cancer: humoral immune response to antigenic determinants coded by the differentially expressed transcripts FLJ23438 and VAMP3

BACKGROUND

Here we evaluate auto-antibody response against two potential antigenic determinants of genes highly expressed in low Gleason Score prostate cancer (PC) tumor samples, namely FLJ23438 and VAMP3.

METHODS

RT-PCR assays were used to analyze mRNA expression profiles of FLJ23438 and VAMP3 transcripts. The auto-antibody response against FLJ23438 and VAMP3 recombinant proteins was tested by immunoblot assays using PC, benign prostate hyperplasia (BPH), healthy donors (HD), and other human cancers plasma samples.

RESULTS

Our data showed that 37% (10/27) and 7.4% (2/27) of PC plasma samples presented auto-antibodies against FLJ23438 and VAMP3, respectively. Only 8.3% (1/12) of BPH plasma samples were reactive for both auto-antibodies, while none (0/12) of HD plasma samples tested were reactive.

CONCLUSIONS

The prevalence of 37% of positive PC plasma samples for anti-FLJ23438 antibodies suggests that humoral immune response against this antigenic determinant could be a potential serum marker for this cancer.

Identification of candidates for tumor-specific alternative splicing in the thyroid

Alternative splicing is the differential processing of exon junctions to produce a new transcript variant from one gene. Some aberrant splicing, however, has been shown to be cancer specific. Identification of these specific splice variations will provide important insight into the molecular mechanism of normal cellular physiology as well as the disease processes. To gain knowledge about whether alternative splicing is linked to thyroid tumorigenesis, we used our prediction database to select targets for analysis. Fifteen putatively new alternative splicing isoforms were selected on the basis of their expression in thyroid libraries and/or their origin in genes previously associated with carcinogenesis. Using a set of 66 normal, benign, and malignant thyroid tissue samples, new splicing events were confirmed by RT-PCR for 13 of 15 genes (a validation rate of 87%). In addition, new alternative splicing isoforms not predicted by the system and not previously described in public databases were identified. Five genes (PTPN18, ABI3BP, PFDN5, SULF2, and ST5) presented new and/or additional unpredicted isoforms differentially expressed between malignant and benign or normal thyroid tissues, confirmed by sequencing. PTPN18, ABI3BP, and PFDN5 revealed a statistically significant differential splicing profile. In addition, real-time PCR analysis revealed that expression of an alternative PFDN5 variant was higher in malignant lesions than in benign lesions or normal tissues.

Class distinction between follicular adenomas and follicular carcinomas of the thyroid gland on the basis of their signature expression

BACKGROUND

Nodules of the thyroid gland are observed frequently in patients who undergo ultrasound studies. The majority of these nodules are benign, corresponding to goiters or adenomas, and only a small fraction corresponds to carcinomas. Among thyroid tumors, the diagnosis of follicular adenocarcinomas by preoperative fine-needle aspiration biopsy is a major challenge, because it requires inspection of the entire capsule to differentiate it from adenoma. Consequently, large numbers of patients undergo unnecessary thyroidectomy.

METHODS

Using data from gene expression analysis, the authors applied Fisher linear discriminant analysis and searched for expression signatures of individual samples of adenomas and follicular carcinomas that could be used as molecular classifiers for the precise classification of malignant and nonmalignant lesions.

RESULTS

Fourteen trios of genes were described that fulfilled the criteria for the correct classification of 100% of samples. The robustness of these trios was verified by using leave-1-out cross-validation and bootstrap analyses. The results demonstrated that, by combining trios, better classifiers could be generated that correctly classified >92% of samples.

CONCLUSIONS

The strategy of classifiers based on individual signatures was a useful strategy for distinguishing between samples with very similar expression profiles.

Evolutionary change - patterns and processes

The present review considered: (a) the factors that conditioned the early transition from non-life to life; (b) genome structure and complexity in prokaryotes, eukaryotes, and organelles; (c) comparative human chromosome genomics; and (d) the Brazilian contribution to some of these studies. Understanding the dialectical conflict between freedom and organization is fundamental to give meaning to the patterns and processes of organic evolution.

Gene expression arrays in cancer research: methods and applications

During the last 5 years, the number of papers describing data obtained by microarray technology increased exponentially with about 3000 papers in 2003. Undoubtedly, cancer is by far the disease that received most of the attention as far as the amount of data generated. As array technology is rather new and highly dependent on bioinformatics, mathematics and statistics, a clear understanding of the knowledge and information derived from array-based experiments is not widely appreciated. We shall review herein some of the issues related to the construction of DNA arrays, quantities and heterogeneity of probes and targets, the consequences of the physical characteristics of the probes, data extraction and data analysis as well as the applications of array technology. Our goal is to bring to the general audience, some of the basics of array technology and its possible application in oncology. By discussing some of the basic aspects of the methodology, we hope to stimulate criticism concerning the conclusions proposed by authors, especially in the light of the very low degree of reproducibility already proven when commercially available platforms were compared . Regardless of its pitfalls, it is unquestionable that array technology will have a great impact in the management of cancer and its applications will range from the discovery of new drug targets, new molecular tools for diagnosis and prognosis as well as for a tailored treatment that will take into account the molecular determinants of a given tumor. Hence, we shall also highlight some of the already available and promising applications of array technology on the day-to-day practice of oncology.

Libraries enriched for alternatively spliced exons reveal splicing patterns in melanocytes and melanomas

It is becoming increasingly clear that alternative splicing enables the complex development and homeostasis of higher organisms. To gain a better understanding of how splicing contributes to regulatory pathways, we have developed an alternative splicing library approach for the identification of alternatively spliced exons and their flanking regions by alternative splicing sequence enriched tags sequencing. Here, we have applied our approach to mouse melan-c melanocyte and B16-F10Y melanoma cell lines, in which 5,401 genes were found to be alternatively spliced. These genes include those encoding important regulatory factors such as cyclin D2, Ilk, MAPK12, MAPK14, RAB4, melastatin 1 and previously unidentified splicing events for 436 genes. Real-time PCR further identified cell line-specific exons for Tmc6, Abi1, Sorbs1, Ndel1 and Snx16. Thus, the ASL approach proved effective in identifying splicing events, which suggest that alternative splicing is important in melanoma development.

Medical genetics services in the city of Sao Paulo, Brazil

The city of Sao Paulo is located in the center of a metropolitan area with nearly 18 million inhabitants and 300,000 births/year. The currently existing medical genetics services are unable to meet the demand, due to their insufficient physical and personnel infrastructure. Institutions and experts in medical genetics could give short training and refresher courses to health professionals to enable them to work in the public health network. The city has a reasonably well developed health care network, represented by the Single Health System (Sistema Unico de Saude - SUS) and by the Family Health Program (Programa de Saude da Familia - PSF). The financial resources for such actions originate in the budget of the managing agencies of such systems. The limitations of genetic services provided to the population of the city could be overcome in a short period of time by developing programs within the public health care network. The city has institutions, professionals and financial resources to make this project feasible. To that end, the competent authorities of the Sao Paulo State and City Secretariats of Health should take managerial responsibility for the genetic services in the city.

Genetic Services in Latin America

This special issue of Community Genetics reviews some of the most important developments in medical genetics in key countries of Latin America. Contributions to this issue were prepared for a special consultation of the World Health Organization held in Porto Alegre, Brazil, on June 19, 2003. Latin America is a region of medium- to low-income countries characterized by socioeconomic problems, with large segments of the population living in poverty and extreme disparities in the distribution of wealth. A rise in chronic diseases typical of the processes of industrialization and urbanization coexists with the persistence of nutritional and infectious diseases characteristic of poverty and underdevelopment. Over the last 2 decades of the 20th century, birth defects and genetic disorders have increased their share of morbidity and mortality, and tertiary-care-based genetic services have developed in urban areas. Although privatization of health care is eroding the public sector, the public institutions continue to be the main providers of genetic services for the bulk of the population and the leaders in research. The development of clinical genetics in the region is concentrated in tertiary-care centers in large cities, although a recent trend began extending genetic services to the community.

The use of Open Reading frame ESTs (ORESTES) for analysis of the honey bee transcriptome

Background

The ongoing efforts to sequence the honey bee genome require additional initiatives to define its transcriptome. Towards this end, we employed the Open Reading frame ESTs (ORESTES) strategy to generate profiles for the life cycle of Apis mellifera workers.

Results

Of the 5,021 ORESTES, 35.2% matched with previously deposited Apis ESTs. The analysis of the remaining sequences defined a set of putative orthologs whose majority had their best-match hits with Anopheles and Drosophila genes. CAP3 assembly of the Apis ORESTES with the already existing 15,500 Apis ESTs generated 3,408 contigs. BLASTX comparison of these contigs with protein sets of organisms representing distinct phylogenetic clades revealed a total of 1,629 contigs that Apis mellifera shares with different taxa. Most (41%) represent genes that are in common to all taxa, another 21% are shared between metazoans (Bilateria), and 16% are shared only within the Insecta clade. A set of 23 putative genes presented a best match with human genes, many of which encode factors related to cell signaling/signal transduction. 1,779 contigs (52%) did not match any known sequence. Applying a correction factor deduced from a parallel analysis performed with Drosophila melanogaster ORESTES, we estimate that approximately half of these no-match ESTs contigs (22%) should represent Apis-specific genes.

Conclusions

The versatile and cost-efficient ORESTES approach produced minilibraries for honey bee life cycle stages. Such information on central gene regions contributes to genome annotation and also lends itself to cross-transcriptome comparisons to reveal evolutionary trends in insect genomes.

A comprehensive transcript index of the human genome generated using microarrays and computational approaches

BACKGROUND

Computational and microarray-based experimental approaches were used to generate a comprehensive transcript index for the human genome. Oligonucleotide probes designed from approximately 50,000 known and predicted transcript sequences from the human genome were used to survey transcription from a diverse set of 60 tissues and cell lines using ink-jet microarrays. Further, expression activity over at least six conditions was more generally assessed using genomic tiling arrays consisting of probes tiled through a repeat-masked version of the genomic sequence making up chromosomes 20 and 22.

RESULTS

The combination of microarray data with extensive genome annotations resulted in a set of 28,456 experimentally supported transcripts. This set of high-confidence transcripts represents the first experimentally driven annotation of the human genome. In addition, the results from genomic tiling suggest that a large amount of transcription exists outside of annotated regions of the genome and serves as an example of how this activity could be measured on a genome-wide scale.

CONCLUSIONS

These data represent one of the most comprehensive assessments of transcriptional activity in the human genome and provide an atlas of human gene expression over a unique set of gene predictions. Before the annotation of the human genome is considered complete, however, the previously unannotated transcriptional activity throughout the genome must be fully characterized.

Antisense intronic non-coding RNA levels correlate to the degree of tumor differentiation in prostate cancer

A large fraction of transcripts are expressed antisense to introns of known genes in the human genome. Here we show the construction and use of a cDNA microarray platform enriched in intronic transcripts to assess their biological relevance in pathological conditions. To validate the approach, prostate cancer was used as a model, and 27 patient tumor samples with Gleason scores ranging from 5 to 10 were analyzed. We find that a considerably higher fraction (6.6%, [23/346]) of intronic transcripts are significantly correlated (P< or =0.001) to the degree of prostate tumor differentiation (Gleason score) when compared to transcripts from unannotated genomic regions (1%, [6/539]) or from exons of known genes (2%, [27/1369]). Among the top twelve transcripts most correlated to tumor differentiation, six are antisense intronic messages as shown by orientation-specific RT-PCR or Northern blot analysis with strand-specific riboprobe. Orientation-specific real-time RT-PCR with six tumor samples, confirmed the correlation (P=0.024) between the low/high degrees of tumor differentiation and antisense intronic RASSF1 transcript levels. The need to use intron arrays to reveal the transcriptome profile of antisense intronic RNA in cancer has clearly emerged.

Characterization of open reading frame-expressed sequence tags generated from Bos indicus and B. taurus mammary gland cDNA libraries

Sequence-based gene expression data are used to interpret results from functional genomic and proteomics studies. Although more than 300 000 bovine-expressed sequence tags (ESTs) are available in public databases, a more thorough and directed sampling of the expressed genome is needed to identify new transcripts and improve assembly and annotation of existing transcript sequences. Accordingly, we examined the utility of constructing cDNA libraries synthesized by arbitrarily primed RT-PCR of mRNA from tissues not well represented in the publicly available bovine EST database. A total of 33 cDNA libraries were constructed from healthy and infected mammary gland tissues of Brazilian Gir and Holstein cattle. This series of libraries was used to generate 6481 open reading frame-expressed sequence tags (ORESTES) that assembled into 1798 unique sequence elements of which, 1157 did not significantly match sequence assemblies available in the Bos taurus gene index. However, a total of 264 of these 1157 sequence elements aligned with mouse and human expressed sequences demonstrating that ORESTES is an effective resource for discovery of novel expressed sequences in cattle. Furthermore, comparison of the alignment position of bovine ORESTES-derived sequence elements to human gene reference sequences suggested that the priming events for cDNA synthesis more often occurred at the central portion of a transcript, which may have contributed to the relatively high rate of novel sequence discovery.

A transcript finishing initiative for closing gaps in the human transcriptome

We report the results of a transcript finishing initiative, undertaken for the purpose of identifying and characterizing novel human transcripts, in which RT-PCR was used to bridge gaps between paired EST clusters, mapped against the genomic sequence. Each pair of EST clusters selected for experimental validation was designated a transcript finishing unit (TFU). A total of 489 TFUs were selected for validation, and an overall efficiency of 43.1% was achieved. We generated a total of 59,975 bp of transcribed sequences organized into 432 exons, contributing to the definition of the structure of 211 human transcripts. The structure of several transcripts reported here was confirmed during the course of this project, through the generation of their corresponding full-length cDNA sequences. Nevertheless, for 21% of the validated TFUs, a full-length cDNA sequence is not yet available in public databases, and the structure of 69.2% of these TFUs was not correctly predicted by computer programs. The TF strategy provides a significant contribution to the definition of the complete catalog of human genes and transcripts, because it appears to be particularly useful for identification of low abundance transcripts expressed in a restricted set of tissues as well as for the delineation of gene boundaries and alternatively spliced isoforms.

Using the miraEST assembler for reliable and automated mRNA transcript assembly and SNP detection in sequenced ESTs

We present an EST sequence assembler that specializes in reconstruction of pristine mRNA transcripts, while at the same time detecting and classifying single nucleotide polymorphisms (SNPs) occuring in different variations thereof. The assembler uses iterative multipass strategies centered on high-confidence regions within sequences and has a fallback strategy for using low-confidence regions when needed. It features special functions to assemble high numbers of highly similar sequences without prior masking, an automatic editor that edits and analyzes alignments by inspecting the underlying traces, and detection and classification of sequence properties like SNPs with a high specificity and a sensitivity down to one mutation per sequence. In addition, it includes possibilities to use incorrectly preprocessed sequences, routines to make use of additional sequencing information such as base-error probabilities, template insert sizes, strain information, etc., and functions to detect and resolve possible misassemblies. The assembler is routinely used for such various tasks as mutation detection in different cell types, similarity analysis of transcripts between organisms, and pristine assembly of sequences from various sources for oligo design in clinical microarray experiments.

Gopinath G, Graudens E, Hahn Y, Han M, Han ZG, Hanada K, Hanaoka H, Harada E, Hashimoto K, Hinz U, Hirai M, Hishiki T, Hopkinson I, Imbeaud S, Inoko H, Kanapin A, Kaneko Y, Kasukawa T, Kelso J, Kersey P, Kikuno R, Kimura K, Korn B, Kuryshev V, Makalowska I, Makino T, Mano S, Mariage-Samson R, Mashima J, Matsuda H, Mewes HW, Minoshima S, Nagai K, Nagasaki H, Nagata N, Nigam R, Ogasawara O, Ohara O, Ohtsubo M, Okada N, Okido T, Oota S, Ota M, Ota T, Otsuki T, Piatier-Tonneau D, Poustka A, Ren SX, Saitou N, Sakai K, Sakamoto S, Sakate R, Schupp I, Servant F, Sherry S, Shiba R, Shimizu N, Shimoyama M, Simpson AJ, Soares B, Steward C, Suwa M, Suzuki M, Takahashi A, Tamiya G, Tanaka H, Taylor T, Terwilliger JD, Unneberg P, Veeramachaneni V, Watanabe S, Wilming L, Yasuda N, Yoo HS, Stodolsky M, Makalowski W, Go M, Nakai K, Takagi T, Kanehisa M, Sakaki Y, Quackenbush J, Okazaki Y, Hayashizaki Y, Hide W, Chakraborty R, Nishikawa K, Sugawara H, Tateno Y, Chen Z, Oishi M, Tonellato P, Apweiler R, Okubo K, Wagner L, Wiemann S, Strausberg RL, Isogai T, Auffray C, Nomura N, Gojobori T, Sugano S. Integrative annotation of 21,037 human genes validated by full-length cDNA clones

The human genome sequence defines our inherent biological potential; the realization of the biology encoded therein requires knowledge of the function of each gene. Currently, our knowledge in this area is still limited. Several lines of investigation have been used to elucidate the structure and function of the genes in the human genome. Even so, gene prediction remains a difficult task, as the varieties of transcripts of a gene may vary to a great extent. We thus performed an exhaustive integrative characterization of 41,118 full-length cDNAs that capture the gene transcripts as complete functional cassettes, providing an unequivocal report of structural and functional diversity at the gene level. Our international collaboration has validated 21,037 human gene candidates by analysis of high-quality full-length cDNA clones through curation using unified criteria. This led to the identification of 5,155 new gene candidates. It also manifested the most reliable way to control the quality of the cDNA clones. We have developed a human gene database, called the H-Invitational Database (H-InvDB; http://www.h-invitational.jp/). It provides the following: integrative annotation of human genes, description of gene structures, details of novel alternative splicing isoforms, non-protein-coding RNAs, functional domains, subcellular localizations, metabolic pathways, predictions of protein three-dimensional structure, mapping of known single nucleotide polymorphisms (SNPs), identification of polymorphic microsatellite repeats within human genes, and comparative results with mouse full-length cDNAs. The H-InvDB analysis has shown that up to 4% of the human genome sequence (National Center for Biotechnology Information build 34 assembly) may contain misassembled or missing regions. We found that 6.5% of the human gene candidates (1,377 loci) did not have a good protein-coding open reading frame, of which 296 loci are strong candidates for non-protein-coding RNA genes. In addition, among 72,027 uniquely mapped SNPs and insertions/deletions localized within human genes, 13,215 nonsynonymous SNPs, 315 nonsense SNPs, and 452 indels occurred in coding regions. Together with 25 polymorphic microsatellite repeats present in coding regions, they may alter protein structure, causing phenotypic effects or resulting in disease. The H-InvDB platform represents a substantial contribution to resources needed for the exploration of human biology and pathology.

RASL11A, member of a novel small monomeric GTPase gene family, is down-regulated in prostate tumors

We performed a genome-wide search for novel loci encoding for Ras-related proteins based on the genome mapping coordinates of the cancer-derived EST dataset at GenBank. Partial sequences from two novel human genes were identified and subsequently used for full length transcript cloning. RASL11A and ARL9 belong to two novel subfamilies coding for small GTPases that we found to be highly conserved among eukaryotes. The Arl9/Arl10 subfamily displays a conserved interswitch toggle that places it evolutionarily closer to the Arf family. Rasl11 proteins are more closely related to the Ras branch of GTPases. All orthologues newly identified here exhibit an Asn residue in place of the highly conserved Thr35 of the G domain, suggesting that the universal switch mechanism of small GTPases may be structurally different in this subfamily. We determined by Northern blot that RASL11A is transcribed in several human tissues and that it is down-regulated in prostate tumors as measured by quantitative real-time PCR. These results highlight a previously uncharacterized subfamily of Ras-related genes that may have a tumor suppressor role in prostate cancer.

Marcadores biológicos de tumores tiroidianos

Um marcador biológico ideal deve ser específico e sensível para identificar o tipo tumoral e caracterizar o estágio da progressão neoplásica. Os tumores de tiróide originam-se de dois tipos celulares: 1) carcinoma medular originário de células parafoliculares; e 2) as neoplasias de células epiteliais foliculares, que incluem bócio, adenomas, carcinomas diferenciados (carcinoma papilífero e carcinoma folicular) e carcinoma indiferenciado (carcinoma anaplásico). O comportamento biológico distinto faz com que cada tipo tumoral necessite de uma conduta terapêutica específica. O conhecimento acumulado ao longo destes anos, utilizando métodos de biologia molecular e, mais recentemente, a genômica, identificou mutações específicas de câncer de tiróide e, atualmente, entendemos muito das alterações que ocorrem na expressão de fatores de crescimento, seus receptores e proteínas sinalizadoras intracelular nas neoplasias tiroidianas. Contudo, apesar desses, até o momento não dispomos de um marcador eficiente que auxilie no diagnóstico e prognóstico e, conseqüentemente, para indicação de uma terapêutica mais adequada. Nesta revisão, discutiremos os principais aspectos relacionados à tumorigênese tiroidiana, avaliando o potencial destes fatores como marcador em neoplasia folicular de tiróide.

trome, trEST and trGEN: databases of predicted protein sequences

We previously introduced two new protein databases (trEST and trGEN) of hypothetical protein sequences predicted from EST and HTG sequences, respectively. Here, we present the updates made on these two databases plus a new database (trome), which uses alignments of EST data to HTG or full genomes to generate virtual transcripts and coding sequences. This new database is of higher quality and since it contains the information in a much denser format it is of much smaller size. These new databases are in a Swiss-Prot-like format and are updated on a weekly basis (trEST and trGEN) or every 3 months (trome). They can be downloaded by anonymous ftp from ftp://ftp.isrec.isb-sib.ch/pub/databases.

Genome-wide survey of human alternative pre-mRNA splicing with exon junction microarrays

Alternative pre-messenger RNA (pre-mRNA) splicing plays important roles in development, physiology, and disease, and more than half of human genes are alternatively spliced. To understand the biological roles and regulation of alternative splicing across different tissues and stages of development, systematic methods are needed. Here, we demonstrate the use of microarrays to monitor splicing at every exon-exon junction in more than 10,000 multi-exon human genes in 52 tissues and cell lines. These genome-wide data provide experimental evidence and tissue distributions for thousands of known and novel alternative splicing events. Adding to previous studies, the results indicate that at least 74% of human multi-exon genes are alternatively spliced.

High frequency of alternative splicing of human genes participating in the HIV-1 life cycle: a model using TSG101, betaTrCP, PPIA, INI1, NAF1, and PML

Alternative splicing may generate splice forms with different biologic roles or missing protein domains implicated in the interaction with HIV-1. To address this issue, 6 human genes were investigated-tumor suppressor gene 101 (TSG101), beta-transducin repeats containing protein (betaTrCP), peptidyl-proly cis-trans isomerase, cyclophilin A (PPIA), integrase interactor 1 protein (INI1), Nef-associated factor 1 (NAF1), and promyelacytic leukemia (PML)-involved in the viral life cycle and HIV-1 pathogenesis. All 6 genes presented alternative splicing, and a combined bioinformatic and reverse transcription polymerase chain reaction (RT-PCR) analysis identified 27 new variants for a total of 53 splice forms (an average of 9 variants per gene). The predicted frequency of the various splice forms based on expressed sequence tags (EST) analysis corresponded to the semiquantitative findings on RT-PCR analysis for the cell culture systems and for native CD4 cells investigated. Interindividual variation in the frequencies of various splice forms in CD4 T cells from blood donors was observed for INI1. Cell type-specific variation of splice pattern was observed for NAF1. Eight splice forms lacked or modified motifs implicated in the interaction with HIV-1, underscoring the potential interest of assessing alternative splicing when investigating viral cell biology and pathogenesis.