Identification of human exons overexpressed in tumors through the use of genome and expressed sequence data

Meta-analytical biomarker search of EST expression data reveals three differentially expressed candidates

Background

Researches have been conducted for the identification of differentially expressed genes (DEGs) by generating and mining of cDNA expressed sequence tags (ESTs) for more than a decade. Although the availability of public databases make possible the comprehensive mining of DEGs among the ESTs from multiple tissue types, existing studies usually employed statistics suitable only for two categories. Multi-class test has been developed to enable the finding of tissue specific genes, but subsequent search for cancer genes involves separate two-category test only on the ESTs of the tissue of interest. This constricts the amount of data used. On the other hand, simple pooling of cancer and normal genes from multiple tissue types runs the risk of Simpson's paradox. Here we presented a different approach which searched for multi-cancer DEG candidates by analyzing all pertinent ESTs in all categories and narrowing down the cancer biomarker candidates via integrative analysis with microarray data and selection of secretory and membrane protein genes as well as incorporation of network analysis. Finally, the differential expression patterns of three selected cancer biomarker candidates were confirmed by real-time qPCR analysis.

Results

Seven hundred and twenty three primary DEG candidates (p-value < 0.05 and lower bound of confidence interval of odds ratio ≧ 1.65) were selected from a curated EST database with the application of Cochran-Mantel-Haenszel statistic (CMH). GeneGO analysis results indicated this set as neoplasm enriched. Cross-examination with microarray data further narrowed the list down to 235 genes, among which 96 had membrane or secretory annotations. After examined the candidates in protein interaction network, public tissue expression databases, and literatures, we selected three genes for further evaluation by real-time qPCR with eight major normal and cancer tissues. The higher-than-normal tissue expression of COL3A1, DLG3, and RNF43 in some of the cancer tissues is in agreement with our in silico predictions.

Conclusions

Searching digitized transcriptome using CMH enabled us to identify multi-cancer differentially expressed gene candidates. Our methodology demonstrated simultaneously analysis for cancer biomarkers of multiple tissue types with the EST data. With the revived interest in digitizing the transcriptomes by NGS, cancer biomarkers could be more precisely detected from the ESTs. The three candidates identified in this study, COL3A1, DLG3, and RNF43, are valuable targets for further evaluation with a larger sample size of normal and cancer tissue or serum samples.

SAGE analysis highlights the putative role of underexpression of ribosomal proteins in GH-secreting pituitary adenomas

BACKGROUND

Although the molecular pathogenesis of pituitary adenomas has been assessed by several different techniques, it still remains partially unclear. Ribosomal proteins (RPs) have been recently related to human tumorigenesis, but they have not yet been evaluated in pituitary tumorigenesis.

OBJECTIVE

The aim of this study was to introduce serial analysis of gene expression (SAGE), a high-throughput method, in pituitary research in order to compare differential gene expression.

METHODS

Two SAGE cDNA libraries were constructed, one using a pool of mRNA obtained from five GH-secreting pituitary tumors and another from three normal pituitaries. Genes differentially expressed between the libraries were further validated by real-time PCR in 22 GH-secreting pituitary tumors and in 15 normal pituitaries.

RESULTS

Computer-generated genomic analysis tools identified 13,722 and 14,993 exclusive genes in normal and adenoma libraries respectively. Both shared 6497 genes, 2188 were underexpressed and 4309 overexpressed in tumoral library. In adenoma library, 33 genes encoding RPs were underexpressed. Among these, RPSA, RPS3, RPS14, and RPS29 were validated by real-time PCR.

CONCLUSION

We report the first SAGE library from normal pituitary tissue and GH-secreting pituitary tumor, which provide quantitative assessment of cellular transcriptome. We also validated some downregulated genes encoding RPs. Altogether, the present data suggest that the underexpression of the studied RP genes possibly collaborates directly or indirectly with other genes to modify cell cycle arrest, DNA repair, and apoptosis, leading to an environment that might have a putative role in the tumorigenesis, introducing new perspectives for further studies on molecular genesis of somatotrophinomas.

SPLOOCE: a new portal for the analysis of human splicing variants

Understanding alternative splicing is crucial to elucidate the mechanisms behind several biological phenomena, including diseases. The huge amount of expressed sequences available nowadays represents an opportunity and a challenge to catalog and display alternative splicing events (ASEs). Although several groups have faced this challenge with relative success, we still lack a computational tool that uses a simple and straightforward method to retrieve, name and present ASEs. Here we present SPLOOCE, a portal for the analysis of human splicing variants. SPLOOCE uses a method based on regular expressions for retrieval of ASEs. We propose a simple syntax that is able to capture the complexity of ASEs.

Mammary gland selective excision of c-jun identifies its role in mRNA splicing

The c-jun gene regulates cellular proliferation and apoptosis via direct regulation of cellular gene expression. Alternative splicing of pre-mRNA increases the diversity of protein functions, and alternate splicing events occur in tumors. Here, by targeting the excision of the endogenous c-jun gene within the mouse mammary epithelium, we have identified its selective role as an inhibitor of RNA splicing. Microarray-based assessment of gene expression, on laser capture microdissected c-jun(-/-) mammary epithelium, showed that endogenous c-jun regulates the expression of approximately 50 genes governing RNA splicing. In addition, genome-wide splicing arrays showed that endogenous c-jun regulated the alternate exon of approximately 147 genes, and 18% of these were either alternatively spliced in human tumors or involved in apoptosis. Endogenous c-jun also was shown to reduce splicing activity, which required the c-jun dimerization domain. Together, our findings suggest that c-jun directly attenuates RNA splicing efficiency, which may be of broad biologic importance as alternative splicing plays an important role in both cancer development and therapy resistance.

Global enhancement of nuclear localization-dependent nuclear transport in transformed cells

Fundamental to eukaryotic cell function, nucleocytoplasmic transport can be regulated at many levels, including through modulation of the importin/exportin (Imp/Exp) nuclear transport machinery itself. Although Imps/Exps are overexpressed in a number of transformed cell lines and patient tumor tissues, the efficiency of nucleocytoplasmic transport in transformed cell types compared with nontransformed cells has not been investigated. Here we use quantitative live cell imaging of 3 isogenic nontransformed/transformed cell pairs to show that nuclear accumulation of nuclear localization signal (NLS)-containing proteins, but not their NLS-mutated derivatives, is increased up to 7-fold in MCF10CA1h human epithelial breast carcinoma cells and in simian virus 40 (SV40)-transformed fibroblasts of human and monkey origin, compared with their nontransformed counterparts. The basis for this appears to be a significantly faster rate of nuclear import in transformed cell types, as revealed by analysis using fluorescence recovery after photobleaching for the human MCF10A/MCF10CA1h cell pair. Nuclear accumulation of NLS/nuclear export signal-containing (shuttling) proteins was also enhanced in transformed cell types, experiments using the nuclear export inhibitor leptomycin B demonstrating that efficient Exp-1-mediated nuclear export was not impaired in transformed compared with nontransformed cells. Enhanced nuclear import and export efficiencies were found to correlate with 2- to 4-fold higher expression of specific Imps/Exps in transformed cells, as indicated by quantitative Western blot analysis, with ectopic expression of Imps able to enhance NLS nuclear accumulation levels up to 5-fold in nontransformed MCF10A cells. The findings indicate that transformed cells possess altered nuclear transport properties, most likely due to the overexpression of Imps/Exps. The findings have important implications for the development of tumor-specific drug nanocarriers in anticancer therapy.

Alternative splicing and genetic diversity: silencers are more frequently modified by SNVs associated with alternative exon/intron borders

With the availability of a large amount of genomic data it is expected that the influence of single nucleotide variations (SNVs) in many biological phenomena will be elucidated. Here, we approached the problem of how SNVs affect alternative splicing. First, we observed that SNVs and exonic splicing regulators (ESRs) independently show a biased distribution in alternative exons. More importantly, SNVs map more frequently in ESRs located in alternative exons than in ESRs located in constitutive exons. By looking at SNVs associated with alternative exon/intron borders (by their common presence in the same cDNA molecule), we observed that a specific type of ESR, the exonic splicing silencers (ESSs), are more frequently modified by SNVs. Our results establish a clear association between genetic diversity and alternative splicing involving ESSs.

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.

Digital candidate gene approach (DigiCGA) for identification of cancer genes

The candidate gene approach is one of the most commonly used methods for identifying genes underlying disease traits. Advances in genomics have greatly contributed to the development of this approach in the past decade. More recently, with the explosion of genomic resources accessible via the public Web, digital candidate gene approach (DigiCGA) has emerged as a new development in this field. DigiCGA, an approach still in its infancy, has already achieved some primary success in cancer gene discovery. However, a detailed discussion concerning the applications of DigiCGA in cancer gene identification has not been addressed. This chapter will focus on discussing DigiCGA in a generalized sense and its applications to the identification of cancer genes, including the cancer gene resources, application status, platform and tools, challenges, and prospects.

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.

Global analysis of aberrant pre-mRNA splicing in glioblastoma using exon expression arrays

Background

Tumor-predominant splice isoforms were identified during comparative in silico sequence analysis of EST clones, suggesting that global aberrant alternative pre-mRNA splicing may be an epigenetic phenomenon in cancer. We used an exon expression array to perform an objective, genome-wide survey of glioma-specific splicing in 24 GBM and 12 nontumor brain samples. Validation studies were performed using RT-PCR on glioma cell lines, patient tumor and nontumor brain samples.

Results

In total, we confirmed 14 genes with glioma-specific splicing; seven were novel events identified by the exon expression array (A2BP1, BCAS1, CACNA1G, CLTA, KCNC2, SNCB, and TPD52L2). Our data indicate that large changes (> 5-fold) in alternative splicing are infrequent in gliomagenesis (< 3% of interrogated RefSeq entries). The lack of splicing changes may derive from the small number of splicing factors observed to be aberrantly expressed.

Conclusion

While we observed some tumor-specific alternative splicing, the number of genes showing exclusive tumor-specific isoforms was on the order of tens, rather than the hundreds suggested previously by in silico mining. Given the important role of alternative splicing in neural differentiation, there may be selective pressure to maintain a majority of splicing events in order to retain glial-like characteristics of the tumor cells.

Sense-antisense pairs in mammals: functional and evolutionary considerations

Analysis of a catalog of S-AS pairs in the human and mouse genomes revealed several putative roles for natural antisense transcripts and showed that some are artifacts of cDNA library construction.

Background

A significant number of genes in mammalian genomes are being found to have natural antisense transcripts (NATs). These sense-antisense (S-AS) pairs are believed to be involved in several cellular phenomena.

Results

Here, we generated a catalog of S-AS pairs occurring in the human and mouse genomes by analyzing different sources of expressed sequences available in the public domain plus 122 massively parallel signature sequencing (MPSS) libraries from a variety of human and mouse tissues. Using this dataset of almost 20,000 S-AS pairs in both genomes we investigated, in a computational and experimental way, several putative roles that have been assigned to NATs, including gene expression regulation. Furthermore, these global analyses allowed us to better dissect and propose new roles for NATs. Surprisingly, we found that a significant fraction of NATs are artifacts produced by genomic priming during cDNA library construction.

Conclusion

We propose an evolutionary and functional model in which alternative polyadenylation and retroposition account for the origin of a significant number of functional S-AS pairs in mammalian genomes.

Alternative splicing: a bioinformatics perspective

The degree of diversity at the transcriptome and proteome levels generated by alternative splicing is astonishing. In this review, we discuss several issues related to alternative splicing with a special emphasis on identification strategies based on bioinformatics.

Alternative splicing and differential gene expression in colon cancer detected by a whole genome exon array

Background

Alternative splicing is a mechanism for increasing protein diversity by excluding or including exons during post-transcriptional processing. Alternatively spliced proteins are particularly relevant in oncology since they may contribute to the etiology of cancer, provide selective drug targets, or serve as a marker set for cancer diagnosis. While conventional identification of splice variants generally targets individual genes, we present here a new exon-centric array (GeneChip Human Exon 1.0 ST) that allows genome-wide identification of differential splice variation, and concurrently provides a flexible and inclusive analysis of gene expression.

Results

We analyzed 20 paired tumor-normal colon cancer samples using a microarray designed to detect over one million putative exons that can be virtually assembled into potential gene-level transcripts according to various levels of prior supporting evidence. Analysis of high confidence (empirically supported) transcripts identified 160 differentially expressed genes, with 42 genes occupying a network impacting cell proliferation and another twenty nine genes with unknown functions. A more speculative analysis, including transcripts based solely on computational prediction, produced another 160 differentially expressed genes, three-fourths of which have no previous annotation. We also present a comparison of gene signal estimations from the Exon 1.0 ST and the U133 Plus 2.0 arrays.

Novel splicing events were predicted by experimental algorithms that compare the relative contribution of each exon to the cognate transcript intensity in each tissue. The resulting candidate splice variants were validated with RT-PCR. We found nine genes that were differentially spliced between colon tumors and normal colon tissues, several of which have not been previously implicated in cancer. Top scoring candidates from our analysis were also found to substantially overlap with EST-based bioinformatic predictions of alternative splicing in cancer.

Conclusion

Differential expression of high confidence transcripts correlated extremely well with known cancer genes and pathways, suggesting that the more speculative transcripts, largely based solely on computational prediction and mostly with no previous annotation, might be novel targets in colon cancer. Five of the identified splicing events affect mediators of cytoskeletal organization (ACTN1, VCL, CALD1, CTTN, TPM1), two affect extracellular matrix proteins (FN1, COL6A3) and another participates in integrin signaling (SLC3A2). Altogether they form a pattern of colon-cancer specific alterations that may particularly impact cell motility.

Data analysis of alternative splicing microarrays

The importance of alternative splicing in drug and biomarker discovery is best understood through several example genes. For most genes, the identification, detection and particularly quantification of isoforms in different tissues and conditions remain to be carried out. As a result, the focus in drug and biomarker development is increasingly on high-throughput studies of alternative splicing. Initial strategies for the parallel analysis of alternative splicing by microarrays have been recently published. The design specificities and goals of alternative splicing microarrays, in terms of identification and quantification of multiple mRNAs from one gene, are promoting the development of novel methods of analysis.

Microarray Analysis of Alternative Splicing

Alternative splicing, defined as the generation of multiple RNA transcript species from a common mRNA precursor, is one of the mechanisms for the diversification and expansion of cellular proteins from a smaller set of genes. Current estimates indicate that at least 60% of genes in the human genome exhibit alternative splicing. Over the past decade, alternative splicing has increasingly been recognized as a major regulatory process with a critical role in normal development. Furthermore, the importance of alternative splicing in disease development and treatment is starting to be appreciated. Therefore, an increasing number of high-throughput genomics and proteomics studies are being performed in order to delineate (a) the changes in alternative splicing under various conditions; (b) the properties and functions of protein isoforms; and (c) the splicing and alternative splicing regulation process. Strategies for the parallel analysis of alternative splice forms by microarray experiments have been conceived, and examples have been published. In addition to the differences in microarray probe design, the analysis of microarrays with probes for exons, exon/exon junctions as well as specific splice forms is significantly different from the standard experiment. Several methods are being developed in order to address the particular needs of alternative splicing microarrays. Many reviews have already dealt with alternative splicing. However, high-throughput analysis methods that are becoming increasingly popular have not received much attention. Here, we will provide an overview of the tools and analysis methods that were developed specifically for alternative splicing microarrays described in terms of specific experiments.

Unbalanced alternative splicing and its significance in cancer

Alternative pre-mRNA splicing leads to distinct products of gene expression in development and disease. Antagonistic splice variants of genes involved in differentiation, apoptosis, invasion and metastasis often exist in a delicate equilibrium that is found to be perturbed in tumours. In several recent examples, splice variants that are overexpressed in cancer are expressed as hyper-oncogenic proteins, which often correlate with poor prognosis, thus suggesting improved diagnosis and follow up treatment. Global gene expression technologies are just beginning to decipher the interplay between alternatively spliced isoforms and protein-splicing factors that will lead to identification of the mutations in these trans-acting factors responsible for pathogenic alternative splicing in cancer.