Discovery of novel splice forms and functional analysis of cancer-specific alternative splicing in human expressed sequences

We report here a genome-wide analysis of alternative splicing in 2 million human expressed sequence tags (ESTs), to identify splice forms that are up-regulated in tumors relative to normal tissues. We found strong evidence (P < 0.01) of cancer-specific splice variants in 316 human genes. In total, 78% of the cancer-specific splice forms we detected are confirmed by human-curated mRNA sequences, indicating that our results are not due to random mis-splicing in tumors; 73% of the genes showed the same cancer-specific splicing changes in tissue-matched tumor versus normal datasets, indicating that the vast majority of these changes are associated with tumorigenesis, not tissue specificity. We have confirmed our EST results in an independent set of experimental data provided by human-curated mRNAs (P-value 10^–5.7). Moreover, the majority of the genes we detected have functions associated with cancer (P-value 0.0007), suggesting that their altered splicing may play a functional role in cancer. Analysis of the types of cancer-specific splicing shifts suggests that many of these shifts act by disrupting a tumor suppressor function. Sur prisingly, our data show that for a large number (190 in this study) of cancer-associated genes cloned originally from tumors, there exists a previously uncharacterized splice form of the gene that appears to be predominant in normal tissue.     

Novel splice variants of cyclin E with altered substrate specificity

Cyclin E, a G₁ cyclin, is overexpressed and present in low molecular weight (LMW) isoforms in breast cancer cells and tumor tissues. In this study we have examined the possibility that the shortened mRNA splice variants could give rise to tumor-specific cyclin E LMW proteins. We used the Splice Capture method to identify, enumerate and isolate known spliced mRNAs and to look for previously undetected mRNA forms of cyclin E that might be translated into the LMW proteins. We show that a new splice variant of cyclin E found in tumor cells isolated by the Splice Capture strategy, named Δ48, activates CDK2 more robustly than full-length cyclin E when assayed from transiently transfected cells with the natural substrate GST-Rb. We also found the Splice Capture method to be superior to the conventional RNase protection assay in analyzing the cyclin E mRNA present in normal and tumor cells. Splice Capture enumerated the relative abundance of known forms of cyclin E mRNA and easily discovered new splice variants in both normal and tumor cells. We conclude that the abundance of cyclin E splice variants in cells may represent a novel form of regulation of cyclin E, and if translated they show altered substrate specificity compared to the full length form of cyclin E.     

Isolation and differential tissue distribution of two human cDNAs encoding PDE1 splice variants.

A cDNA selection technique has been used to isolate full-length human cDNAs of the phosphodiesterase 1 (PDE1) calcium calmodulin (CaM)-regulated phosphodiesterase gene family. We isolated cDNAs representing multiple splice variants of PDE1A, 1B and 1C from a variety of tissues. Included among these were two novel splice variants for PDE1A and 1B. The first, PDE1A5, encodes a 519-residue protein, which is different from PDE1A1 by the insertion of 14 residues, a divergent carboxy terminus and also differs from PDE1A3 through a divergent amino terminus. Our second novel splice variant represents the first occurrence of a splice variant of the PDE1B gene. PDE1B2 encodes a 516-residue protein and diverges from PDE1B1 by the replacement of the first 38 residues by an alternative 18, which is predicted to be functionally significant. Using the splice variant sequence differences to perform comparative Northern analysis, we have demonstrated that each variant has a differential tissue distribution.     

Genome-wide detection of tissue-specific alternative splicing in the human transcriptome

We have developed an automated method for discovering tissue-specific regulation of alternative splicing through a genome-wide analysis of expressed sequence tags (ESTs). Using this approach, we have identified 667 tissue-specific alternative splice forms of human genes. We validated our muscle-specific and brain-specific splice forms for known genes. A high fraction (8/10) were reported to have a matching tissue specificity by independent studies in the published literature. The number of tissue-specific alternative splice forms is highest in brain, while eye-retina, muscle, skin, testis and lymph have the greatest enrichment of tissue-specific splicing. Overall, 10-30% of human alternatively spliced genes in our data show evidence of tissue-specific splice forms. Seventy-eight percent of our tissue-specific alternative splices appear to be novel discoveries. We present bioinformatics analysis of several tissue-specific splice forms, including automated protein isoform sequence and domain prediction, showing how our data can provide valuable insights into gene function in different tissues. For example, we have discovered a novel kidney-specific alternative splice form of the WNK1 gene, which appears to specifically disrupt its N-terminal kinase domain and may play a role in PHAII hypertension. Our database greatly expands knowledge of tissue-specific alternative splicing and provides a comprehensive dataset for investigating its functional roles and regulation in different human tissues.     

Alternative splicing in colon, bladder, and prostate cancer identified by exon array analysis

Alternative splicing enhances proteome diversity and modulates cancer-associated proteins. To identify tissue- and tumor-specific alternative splicing, we used the GeneChip Human Exon 1.0 ST Array to measure whole-genome exon expression in 102 normal and cancer tissue samples of different stages from colon, urinary bladder, and prostate. We identified 2069 candidate alternative splicing events between normal tissue samples from colon, bladder, and prostate and selected 15 splicing events for RT-PCR validation, 10 of which were successfully validated by RT-PCR and sequencing. Furthermore 23, 19, and 18 candidate tumor-specific splicing alterations in colon, bladder, and prostate, respectively, were selected for RT-PCR validation on an independent set of 81 normal and tumor tissue samples. In total, seven genes with tumor-specific splice variants were identified (ACTN1, CALD1, COL6A3, LRRFIP2, PIK4CB, TPM1, and VCL). The validated tumor-specific splicing alterations were highly consistent, enabling clear separation of normal and cancer samples and in some cases even of different tumor stages. A subset of the tumor-specific splicing alterations (ACTN1, CALD1, and VCL) was found in all three organs and may represent general cancer-related splicing events. In silico protein predictions suggest that the identified cancer-specific splice variants encode proteins with potentially altered functions, indicating that they may be involved in pathogenesis and hence represent novel therapeutic targets. In conclusion, we identified and validated alternative splicing between normal tissue samples from colon, bladder, and prostate in addition to cancer-specific splicing events in colon, bladder, and prostate cancer that may have diagnostic and prognostic implications.     

Characterization of a novel human breast cancer associated gene (BCA3) encoding an alternatively spliced proline-rich protein

As part of an integrated study of breast cancer gene expression, partial cDNAs were cloned from normal and tumor breast cells by subtractive-hybridization and differential display cloning. The DNA sequence for one of these breast cancer associated genes was used to construct the larger 1319 bp BCA3 cDNA sequence using ESTs without assigned names or functions. High-level BCA3 mRNA expression was found in breast and prostate tumor cell lines whereas normal breast and prostate tissues have low-level expression. Further analysis revealed possible functional domains and alternative splicing of BCA3 that we confirmed by RT-PCR analysis. Immunohistochemistry revealed that the protein is expressed in breast tumor cells in vivo, and not in surrounding stromal tissue.     

Role of splice variants in the metastatic progression of prostate cancer

AS (alternative splicing) and its role in disease, especially cancer, has come to forefront in research over the last few years. Alterations in the ratio of splice variants have been widely observed in cancer. Splice variants of cancer-associated genes have functions that can alter cellular phenotype, ultimately altering metastatic potential. As metastases are the cause of approximately 90% of all human cancer deaths, it is crucial to understand how AS is dysregulated in metastatic disease. We highlight some recent studies into the relationship between altered AS of key genes and the initiation of prostate cancer metastasis.     

Biomarkers of human cutaneous squamous cell carcinoma from tissues and cell lines identified by DNA microarrays and qRT-PCR

Squamous cell carcinoma (SCC) is the second most common form of skin cancer in Caucasians. Here we report on the identification of biomarkers of human cutaneous SCC cell lines in vitro and tissue samples in vivo using DermArray and PharmArray DNA microarrays, consisting of ca. 7400 unique human cDNAs. Differentially expressed genes were identified in two facial skin SCC cell lines (SCC 12 and SCC 13) compared to normal keratinocytes, and three cutaneous SCC tissue samples compared to normal skin. Quantitative validations of up- and down-regulated biomarkers were performed by qRT-PCR on 23 biomarker genes for the cell lines and 20 biomarker genes for the tumor tissues. In addition, three oral SCC cell lines were also included in the qRT-PCR validations for comparison, and the biomarker profiles were highly similar between the cutaneous and the oral SCC cell lines for all 23 biomarkers examined. The expression profiles for a variety of non-cutaneous SCC types, such as head-and-neck, oral, and lung, have been previously published. This report is the first to describe biomarkers for cutaneous SCC in two contexts, in vitro and in vivo. Although there was minimal overlap between the two different contexts using DNA microarrays, five genes were found common to both the cell lines and tissues, namely fibronectin 1, annexin A5, glyceraldehyde 3-phosphate dehydrogenase, zinc-finger protein 254, and huntingtin-associated protein interacting protein. Some of our previously published biomarkers of normal keratinocytes were down-regulated in SCC, suggestive of the dedifferentiated status of the transformed cells. While recent reports have identified some of the same genes as SCC biomarkers, for instance in head-and-neck cancer, thereby validating our approach, we have identified some novel biomarkers for cutaneous disease. These biomarker lists may be useful in molecular diagnostics of non-melanoma skin cancer, and a subset of the biomarkers might serve as suitable targets for drug discovery efforts of therapies for SCC.     

Use of SMART-generated cDNA for gene expression studies in multiple human tumors

We demonstrate here that SMART PCR-amplified cDNAs arrayed on a nylon membrane are suitable for high-throughput tissue expression profiling when starting biological materials are limited. We show that SMART cDNA accurately reflects gene expression patterns found in total RNA by comparing the expression level of several target genes in SMART PCR-amplified cDNAs and their corresponding total RNAs. We also arrayed cDNAs from 68 matched tumor and normal samples on a nylon membrane to determine whether SMART PCR-amplified cDNA could be used for detecting differentially expressed genes in these tissues. These arrays containing normalized tumor and normal cDNAs were hybridized with probes for glutathione peroxidase and gelsolin. The hybridization results revealed cancer-related and patient-specific gene expression differences between tumor and normal tissues for these genes. These studies show that SMART PCR-amplified cDNAs maintain the complexity of the original mRNA population and are thus suitable for high-throughput studies to compare the relative abundance of target genes and to detect differentially expressed genes in a wide variety of tissues simultaneously.     

Identification and characterization of a novel cancer/testis antigen gene CAGE

We applied serological analysis of cDNA expression library technique to identify cancer-associated genes. We screened cDNA expression libraries of human testis and gastric cancer cell lines with sera of patients with gastric cancers. We identified a gene whose expression is testis-specific among normal tissues. We cloned and characterized this novel gene. It contains D-E-A-D box domain and encodes a putative protein of 630 amino acids with possible helicase activity. It showed wide expression in various cancer tissues and cancer cell lines. The corresponding gene was named cancer-associated gene (CAGE). PCR of human x hamster Radiation Hybrids showed localization of CAGE on the human chromosome Xp22. Transient transfection of CAGE showed predominantly nuclear localization. Both Western blot and plaque assay indicated seroreactivity of CAGE protein. We found that demethylation played a role in the activation of CAGE in some cancer cell lines that do not express it. Cell synchronization experiments showed that the expression of CAGE was related with cell cycle. This suggests that CAGE might play a role in cellular proliferation. Because CAGE is expressed in a variety of cancers but not in normal tissues except testis, this gene can be a target of antitumor immunotherapy.     

Analysis of canonical and non-canonical splice sites in mammalian genomes

A set of 43 337 splice junction pairs was extracted from mammalian GenBank annotated genes. Expressed sequence tag (EST) sequences support 22 489 of them. Of these, 98.71% contain canonical dinucleotides GT and AG for donor and acceptor sites, respectively; 0.56% hold non-canonical GC-AG splice site pairs; and the remaining 0.73% occurs in a lot of small groups (with a maximum size of 0.05%). Studying these groups we observe that many of them contain splicing dinucleotides shifted from the annotated splice junction by one position. After close examination of such cases we present a new classification consisting of only eight observed types of splice site pairs (out of 256 a priori possible combinations). EST alignments allow us to verify the exonic part of the splice sites, but many non-canonical cases may be due to intron sequencing errors. This idea is given substantial support when we compare the sequences of human genes having non-canonical splice sites deposited in GenBank by high throughput genome sequencing projects (HTG). A high proportion (156 out of 171) of the human non-canonical and EST-supported splice site sequences had a clear match in the human HTG. They can be classified after corrections as: 79 GC-AG pairs (of which one was an error that corrected to GC-AG), 61 errors that were corrected to GT-AG canonical pairs, six AT-AC pairs (of which two were errors that corrected to AT-AC), one case was produced from non-existent intron, seven cases were found in HTG that were deposited to GenBank and finally there were only two cases left of supported non-canonical splice sites. If we assume that approximately the same situation is true for the whole set of annotated mammalian non-canonical splice sites, then the 99.24% of splice site pairs should be GT-AG, 0.69% GC-AG, 0.05% AT-AC and finally only 0.02% could consist of other types of non-canonical splice sites. We analyze several characteristics of EST-verified splice sites and build weight matrices for the major groups, which can be incorporated into gene prediction programs. We also present a set of EST-verified canonical splice sites larger by two orders of magnitude than the current one (22 199 entries versus approximately 600) and finally, a set of 290 EST-supported non-canonical splice sites. Both sets should be significant for future investigations of the splicing mechanism.     

Characterisation of tumour-associated antigens in colon cancer

In order to search for clinically relevant cancer-associated genes and to define further the spectrum of immunogenic proteins, we applied SEREX (serological identification of antigens by recombinant expression cloning) to analyse genes expressed in colon adenocarcinoma. Eight different serum-reactive cDNA clones were isolated by immunoscreening from a colon cancer-derived cDNA expression library. mRNA expression studies showed that 2 of them, RHAMM and AD034, have a differential tissue distribution, and that 3 genes, NAP1L1, RHAMM and AD034, are overexpressed in tumours in comparison with the adjacent non-cancerous tissues. 5' RLM-RACE analysis of AD034, a sequence with a tyrosine kinase motif, revealed a frameshifting insertion of 32 bp, most likely generated by use of cryptic splice site in tumour-derived cDNA. Analysis of full-length RHAMM cDNA sequence revealed the presence of two splice variants, which are known to have a different sub-cellular localisation; expression of these splice variants is altered in colon cancer tissues. Serological responses to three antigens (C21ORF2, EPRS and NAP1L1) were found mainly in cancer patients' sera.