- Research article
- Open Access
Verification of predicted alternatively spliced Wnt genes reveals two new splice variants (CTNNB1 and LRP5) and altered Axin-1 expression during tumour progression
BMC Genomicsvolume 7, Article number: 148 (2006)
Splicing processes might play a major role in carcinogenesis and tumour progression. The Wnt pathway is of crucial relevance for cancer progression. Therefore we focussed on the Wnt/β-catenin signalling pathway in order to validate the expression of sequences predicted as alternatively spliced by bioinformatic methods. Splice variants of its key molecules were selected, which may be critical components for the understanding of colorectal tumour progression and may have the potential to act as biological markers. For some of the Wnt pathway genes the existence of splice variants was either proposed (e.g. β-Catenin and CTNNB1) or described only in non-colon tissues (e.g. GSK3β) or hitherto not published (e.g. LRP5).
Both splice variants – normal and alternative form – of all selected Wnt pathway components were found to be expressed in cell lines as well as in samples derived from tumour, normal and healthy tissues. All splice positions corresponded totally with the bioinformatical prediction as shown by sequencing. Two hitherto not described alternative splice forms (CTNNB1 and LRP5) were detected. Although the underlying EST data used for the bioinformatic analysis suggested a tumour-specific expression neither a qualitative nor a significant quantitative difference between the expression in tumour and healthy tissues was detected. Axin-1 expression was reduced in later stages and in samples from carcinomas forming distant metastases.
We were first to describe that splice forms of crucial genes of the Wnt-pathway are expressed in human colorectal tissue. Newly described splicefoms were found for β-Catenin, LRP5, GSK3β, Axin-1 and CtBP1. However, the predicted cancer specificity suggested by the origin of the underlying ESTs was neither qualitatively nor significant quantitatively confirmed. That let us to conclude that EST sequence data can give adequate hints for the existence of alternative splicing in tumour tissues. That no difference in the expression of these splice forms between cancerous tissues and normal mucosa was found, may indicate that the existence of different splice forms is of less significance for cancer formation as suggested by the available EST data. The currently available EST source is still insufficient to clearly deduce colon cancer specificity. More EST data from colon (tumour and healthy) is required to make reliable predictions.
Alternative splicing is an important source for complexity and diversity of eukaryotic proteomes. A couple of databases hold information of published alternatively spliced forms (e.g. ASDB , Alternative Splicing Database ). A new kind of database contains splice variants detected by sequence comparison of transcripts and ESTs. These alignment based databases compare expressed sequence tags (ESTs) either with genomic sequences , with UniGene clusters [4–6] or with Ensembl transcripts . An important advantage of this approach is the possibility to detect undiscovered splice forms . In spite of some (partially solved) problems with EST data (genomic contamination or incomplete mRNA processing; see ), they may reflect different environmental or developmental conditions. Therefore, ESTs are an important source for the analysis of alternative splicing which is often interpreted as a regulatory mechanism for certain stages in the life span of organisms. In addititon, alternative splicing has been found to be associated with various diseases including cancer .
The motivation for the study presented here was (1) to validate the expression of sequences predicted as alternatively spliced by bioinformatic methods, (2) to analyse the expression of alternatively spliced genes which are relevant for pathways involved in cancer formation in colon, and (3) to try to identify splice forms which show a deregulated expression in human carcinomas.
For our experimental studies we focussed on the Wnt/β-catenin signalling pathway, which has a central role for colorectal tumourigenesis. Deregulated activation of this pathway is also observed in hepatocellular carcinomas, melanomas and uterine and ovarian carcinomas. For instance, in about 80% of all colorectal carcinomas, the tumour suppressor gene APC of the Wnt pathway carries a mutation . In general, the Wnt/β-catenin pathway plays an important role for cell proliferation, differentiation, and cell motility. Wnt proteins act by activating either the "classical" Wnt pathway or one of two non-canonical signalling cascades. Constitutive activation of the Wnt signalling pathway due to mutations in the regulator genes of the pathway (β-catenin, APC and axin) leads to an increase in the cytoplasmic level of β-catenin and subsequent transcription of Wnt target genes such as c-Myc, cyclin D1 and WISP-1. An overview of Wnt pathway components which have been found to be mutated in human tumours is presented by Thorstensen and Lothe . According to these authors mutations in GSK-3β, another relevant regulator molecule, have never been reported in human cancers. This might be due to the central role of GSK-3β in pathways other than the Wnt pathway and due to the fact that mutations in this gene may be incompatible with cell viability.
By bioinformatic analysis a list of putative splice variants of relevant molecules of the Wnt pathway genes was generated. RT-PCR was then used to examine gene expression of the putative sequences in several colorectal tumour cell lines, but also in patient samples derived from colon cancer and matching normal tissue.
Results and discussions
Splicing processes might play a major role in carcinogenesis. For example CD44, Wilms tumour gene WT1, BRCA1/2 and kallikrein family genes are well known and encode cancer type related splice variants . Due to the crucial relevance of the Wnt pathway for cancer progression, splice variants of the relevant molecules may be critical components for the understanding of colorectal tumour progression and may have the potential to act as biological markers. For the Wnt signalling pathway, 22 genes out of 25 are predicted to have alternatively spliced forms by EASED. Five genes involved in the Wnt signal transduction pathway were selected with respect to their role in the pathway. RT-PCR was used to examine gene expression of splice forms of β-catenin (CTNNB1), axin-1, glycogen synthase kinase (GSK-3β), C-terminal binding protein 1 (CtBP1) und low density lipoprotein receptor related protein 5 (LRP5) in several colorectal tumour cell lines, but also in patient samples derived from colon cancer and matching normal tissue. All those genes (cf. Table 3) are important key players and differ in their function and in their intracellular localization (Fig. 1). β-catenin (CTNNB1) is carrying two different splice sites. Thus, we analysed 6 EASED splice sites (Table 3). Sequence information is given in Fig. 2A–2F. We took great care in determining optimal RT-PCR conditions in order not to miss any less expressed splice forms. This is demonstrated for GSK-3β (Figure 3). Annealing temperature (3A), as well as cycling time (3B) and annealing/extension time (3C) were optimized in order to detect both splice forms. With the exception of CTNNB1, all splice sites are situated within the coding DNA and lead to changes in amino acid sequences. We used one primer set for the detection of up to three transcripts simultaneously (Figure 4). On the one hand, the presence of the different splice forms can easily compared using this procedure. On the other hand the size of the amplicon might impact the efficiency of the reaction and the expression of longer splice forms might be underestimated. As already known, the quantitative value of "hand-made" RT-PCR is limited and small expression differences will be neglected. Therefore, we used also quantitative RT-PCR (TaqMan) for one of the genes (Axin-1).
The PCR results for the 6 selected splice variants give strong evidence that the algorithm used to predict potential alternative splice forms is very reliable. The selected Wnt pathway was used as a model to test the splice variant prediction algorithm. In summary, we found that always either (1) both putative splice forms or (2) only the normal (longer) isoform are expressed in tumour and in normal mucosa. Although the underlying EST data used for the bioinformatic analysis suggested an exclusive expression in tumours neither a qualitative nor a significant quantitative difference between the expression in tumour and healthy tissues was detected. However, this result is in good agreement with Gupta et al.  who described that predicted tumour-specific expression of isoforms derived from ESTs usually tends not to reflect the experimentally validated expression pattern.
Accumulation of β-catenin, a key molecule of the Wnt pathway, in the nucleus and interaction with TCF/LEF transcription factors leads to transcription activation of a number of target genes such as c-Myc, cyclin D1, fra1, c-Jun, matrilysin, CD44, urokinase-type plasminogen activator receptor and others . These molecules are clearly relevant for tumour formation, because of their roles in proliferation, apoptosis and cell-cycle regulation. High concentration of cytoplasmic β-catenin is critical for growth of colorectal tumour cells . Furthermore, the results of Hlubek et al.  suggest that nuclear β-catenin activates the coordinated expression of the interacting proinvasive proteins laminin gamma 2 chain and MT1-MMP, thereby leading to a promigratory activity at the invasive front of colorectal cancers. Besides its role for tumour formation β-catenin is relevant for cell-cell adhesion in conjunction with E-cadherin . The existence of an additional alternative form of CTNNB1 was proposed by Nollet et al. . We found this second and hitherto not described alternatively spliced isoform (Fig. 2B and Fig. 4). The sequence of this splice variant is registered under GenBank accession number AB062292. This isoform contains the complete intron 15 sequence as well as the complete exon 16. The splice site is not located within the coding region, but may affect regulatory functions of the RNA-sequence.
CtBP1 might selectively repress epithelial cell adhesion and transcription of several pro-apoptotic genes. CtBP interacts with TCF-4, histone deacetylases and is able to decrease the expression of the endogenous Wnt targets . Grooteclaes et al. report that CtBP-knockout cells were hypersensitive to apoptosis . The expression of both CtBP1 splice variants was shown in colorectal tumour as well as in normal mucosa. However, the expression of the normal 63 bp isoform was lower. The alternative isoform contains an insert of 196 bp-length within exon 2. Due to a change of the start of the coding sequence the predicted protein is smaller than that of the normal isoform.
The serine/threonine kinase GSK-3β binds to and phosphorylates several proteins in the Wnt pathway. Phosphorylation of β-catenin by GSK-3β leads to ubiquitination of β-catenin and degradation in proteasomes . No mutations of the GSK-3β gene were found in three colorectal cancer cells , but further studies of cells and cancer tissues are needed to test this hypothesis. The expression of both putative GSK-3β isoforms was detected in colorectal carcinomas and normal mucosa (Fig. 2D and Fig. 5). These isoforms have been annotated in Ensembl (ENST00000264235, ENST00000316626). The proteinkinase-domain (56–353 amino acids [aa]) of the alternative isoform is shorter than that of the normal isoform (13 aa between 304–316 aa) as shown before. This GSK-3β splice variant was already described by Mukai et al. , who localised it in human brain. Our results in colon cancer show that this splice form is not exclusively expressed in brain. Moreover, the hypothesis of Thorstensen and Lothe  that GSK-3α may substitute for loss of GSK-3β in tumour tissue is in agreement with our findings.
LDL receptor related proteins are cell surface molecules recently identified as Wnt co-receptors. It showed no effects on the canonical Wnt signalling by itself, but acted synergistically to Wnt. LRP5 binds axin and recruits axin to the membranes. Supposedly, the translocation of axin leads to stabilisation of β-catenin and target gene transcription . LRP5 binds 30 or more ligands with high affinity and delivers them to lysosomal compartments, where they are degraded. Known ligands include lipoproteins, proteinases, proteinase-inhibitor complexes, matrix-proteins, bacterial toxins and viruses, and various intracellular proteins . The presence of LRP5 correlates significantly with tumour metastasis in final stage osteosarcoma . Fujino et al. showed that LRP5 is also required in the normal cholesterol and glucose metabolism . Both LRP5 isoforms were found to be expressed in colon cancer cell lines and in patient tissues (Fig. 2F, Fig. 7). However, the expression of the alternative isoform in the tissue samples was low. Here, we describe a new, hitherto not published LRP5 isoform. Okubo et al.  has examined the human LRP5 gene in patients with Type III hyperlipoproteinemia by restriction enzyme fragment polymorphism (RFLP) and described seven other sequence polymorphisms, which means that a considerable heterogeneity exists in LRP5-splicing variants. On the protein level, the LRP5 exon skip would lead to the loss of part of the EGF-like -Domain (PF00008). The EGF precursor homology repeats in receptors are important for the dissociation of ligands from receptors in endocytic vesicles. The transmembrane domain is necessary for anchoring to membranes and the cytoplasmic domain is required for their targeting to coated pits and subsequent internalisation . The functional relevance of this loss has to be elucidated in future studies.
Axin's role as a scaffold protein in the Wnt signalling is very important for β-catenin degradation by the proteasome. So, axin acts as a negative regulator and as tumour suppressor . Axin mutations were found in human hepatocarcinomas with intact genes for β-catenin and APC . The expression of axin-1 was analysed by RT-PCR and quantitative real-time PCR. Exon 8 of the alternative isoform shows a 108 bp skip (Fig. 2E). The splice position was examined by RT-PCR with two different primer pairs (Fig. 6). The forward primer A (Table 1) was positioned within the skipped sequence. Therefore only the normal form is amplified by RT-PCR using this primer (Fig. 6, left part). Our results show that in RT-PCR with Axin-primer B (Fig. 6, right part) only the shorter alternative product was detected using different PCR conditions. We suppose that the PCR reaction in this case prefers the amplification of the short alternative product and used RT-PCR to examine this result. Using quantitative real-time PCR, we expected that primer T1 amplifies the normal form and primer T2 the alternative form (Fig. 8 line 1 and line 4). The results of the real-time PCR and the clinical data of patients are shown in table 4, 5. Both axin-1 splice isoforms were found in the tumour tissue as well as in healthy tissue. No differences in expression of either normal or alternative isoforms between tumour and mucosa were found (Fig. 9). However, in nearly all patient samples of normal mucosa and tumour, the alternative axin isoform was less expressed than the normal splice variant (Fig. 9). Axin expression was reduced in later stages (Fig 10). Similarly, axin expression is slightly reduced in samples from carcinomas with distant metastases (Fig. 11).
On the protein level, the alternative isoform is missing 14 aa within the DIX-domain. The DIX domain is a signalling module that is conserved in axin/conductin. In addition to this well-documented function in mediating protein-protein interactions, the DIX domain has been shown to target Dvl to actin stress fibres and vesicular membranes, and is indispensable for JNK activation by axin . Axin is incapable of dimerization if carrying mutations within the DIX-Domain. This leads to the loss of axin's function as a major scaffold protein for degradation of cytoplasmic β-catenin . The expression of axin was reduced in later stages (Fig. 10). Similarly, the expression of axin is slightly reduced in samples from carcinomas with distant metastases (Fig. 11). Expression of axin can be correlated inversely with depth of invasion, lymph node metastasis and lymphatic invasion in esophageal squamous cell carcinoma as shown by Nakajima .
In summary our results demonstrate the potential of bioinformatic prediction for the search of new splice sites (as shown for CTNNB1 and LRP5) and show the precision and reliability of EST based prediction of the existence of splice sites. Otherwise, the predicted cancer specificity suggested by the origin of the underlying ESTs was neither qualitatively nor significant quantitatively confirmed. That let us conclude that EST data can give adequate hints for the existence of alternative splicing in tumour tissues, which is significant information. That no difference in the expression of these splice forms between cancerous tissues and normal mucosa was found, may indicate that the existence of different splice forms is of less significance for cancer formation as previously suggested. One source of incorrect tumour specificity could be the lack of standard and quality control in annotation of the submitted ESTs in the past. Additionally, the differing experimental results from database analysis can be caused to still incomplete EST data and a bias of EST sequences towards some tissue types and pathological states. More EST data from colon cancer as well as from healthy colon is required to have a statistical power to make precise specificity predictions. Future work is necessary to reveal whether there is strict correlation between a loss of e.g. axin expression and the tumour progression or even the formation of distant metastasis.
Prediction of splice variants with bioinformatic methods
The Extended Alternatively Spliced EST Database (EASED)  project is establishing a comprehensive database of alternatively spliced human mRNAs. Moreover, EASED includes useful biological information and provides the possibility to search for biologically relevant data and for candidate genes for the origin of diseases.
All transcript sequences were extracted from the Ensembl  web site and the repetitive elements have been marked with RepeatMasker . The EST sequences were taken from dbEST  and were formatted for BLAST. All transcripts were aligned with the ESTs using WU-Blast  with an eValue of 0.001. The resulting high scoring pairs (HSPs) were analysed with an in-house written Perl-script. If one EST reveals more than one HSP with one transcript the following filtering criteria were used: (i) the gap between two HSPs has to be longer than 10 bp, (ii) both HSPs have to stretch over a sequence length of more than 30 bp, (iii) the overall identity for each HSP have to be at least 90%. The alternative splice event was named 'skip' if the gap between two HSPs was obtain within the transcript sequence. Otherwise it was termed as an 'insert'. All ESTs showing these HSPs with the mentioned filtering criteria were named as 'alternatively spliced ESTs'. The corresponding 'normally spliced ESTs' show one complete HSP without any gap over the entire EST sequence length. If possible, the HSPs from different ESTs that aligned to one transcript were clustered as follows: all HSP with start and end positions referring to the same transcript position (as an uncertainty a difference of not more than 5 bp is allowed) were merged to one single splice event. The resulting splice event positions arise from the mean of all HSP start positions and the mean of all HSP end positions.
The prediction of splice variants of investigated genes is based on Ensembl release 19 (February 2004) with 31609 human transcripts and dbEST release from December 2003 with 5,427,257 human ESTs. All processed and annotated alternative splice forms are stored in a MySQL database. It allows analysing each splice site separately instead of analysing complete transcripts. Each EST was categorised concerning its EST library annotation. At the moment of writing, EASED comprises four disease classes (cancer, normal, healthy, unknown), four development classes (adult, embryo, newborn, unknown) as well as 193 tissue classes. For the web version of EASED the 193 tissue classes were merged to 23 master classes.
In the following, we will indicate the transcript referring to the Ensembl database as the 'normal form'. In contrast, the EST based predicted isoform is named 'alternative form'.
Selection of genes to be analysed
In the first selection step we searched for genes playing a major role in the Wnt/β-catenin signalling pathway with a higher number of alternatively spliced ESTs in colon cancer compared to normally spliced ESTs. To increase the number of genes we further considered genes with a comparable high percentage of alternatively spliced ESTs from other gastrointestinal tissues (as e.g. stomach) or from cancer tissues apart from colon. To further increase the number of candidates we examined the ratio of cancer ESTs to the total number of annotated ESTs. Moreover, we were interested to find out unpublished splice sequences. The selection criteria for the several genes that were chosen to be analysed are given in Table 3.
Patients and tissue samples
Fourteen consecutive cases of colorectal cancers were examined. The patients had been treated in the Clinic for Surgery and Surgical Oncology, Robert-Roessle-Klinik in Berlin, Germany. Tissue samples were collected after obtaining informed consent. Carcinoma tissues were excised carefully from tumour specimens. Non-malignant mucosa was also scraped off in a distance of over 50 mm from tumour margin. Tissue was snap-frozen in liquid nitrogen by the Tumorbank service of the Robert-Roessle-Klinik. Carcinomas were classified and characterised according to the World Health Organization's UICC (Union Internationale Contre le Cancer) system: invasiveness (pT), lymph-node metastases (pN), primary distant metastases (M), grading (G), lymphatic vessel-invasion (L), venous vessel-invasion (V), type of resection (R0-, R1- or R2) (Table 4). Seven of these patients were women and seven were men. The age range was 51–78 years, with the median age of 64.5 years. The range of follow-up was 3–105 months, with a median follow-up of 46 months. The location of the carcinomas was as follows: rectum, four patients; colon ascendens, three patients; caecum, four patients; sigma, two patients, and one patient suffered from a carcinoma in colon transversum. Seven patients had metastases at the time of surgery. In the subsequent 2–4 years after surgery, six patients have died due to their disease. Information on gender, age, stage of disease, and histopathological factors was abstracted from the clinical and pathological records. Furthermore, material obtained during coloscopy from tumour free patients was examined for control purposes.
All cell lines were of human origin, colon adenocarcinoma cell lines HT-29 (ATCC: HTB-38), HCT-116 (ATCC:CCL-247), LS-174T (ATCC: CL-188), SW-480 (ATCC:CCL-228), and mammary carcinoma cells MDA-MB-435. The adenocarcinoma lines HT-29, LS-174T, MDA-MB-435, and SW-480 were grown in Dulbecco's modified Eagle's medium, HCT-116 was grown in Iscove's medium, always supplemented with 10 % FBS, 2 mM glutamine, 100 U penicillin and 0.1 mg streptomycin. Subconfluent adherent cells were harvested by a mixture of trypsin (0.05%) and EDTA (0.02%), rescued with their own medium and washed with phosphate buffered saline (PBS).
Tissue preparation and RNA extraction for RT-PCR
Cryo-sections from colorectal and normal tissue specimens were transferred into RLT-Buffer containing β-mercaptoethanol and frozen at -80°C. Only specimens containing more than 60% epithelial cells, no Peyer's patches and no necrotic areas were further processed. RNA was isolated using QIAshredder and RNeasy Kits. Integrity of isolated mRNA was checked by β-Actin One-Step RT PCR, and RNA concentration was determined photometrically. About 3 μg of RNA were used for cDNA synthesis. In the presence of 4 μM random hexamer primer RNA was incubated for 10 min at 70°C, then the mixture was held on ice for another 10 min. After addition of M-MLV-Buffer, 8 U/μL M-MLV Reverse Transcriptase, 0.8 U/μL RNasin, 0.1 μg/μL BSA and 1.25 mM dNTP and 1 h incubation at 37°C, twice the volume of absolute ethanol was added. Following 30 min at -40°C, samples were spun down in a centrifuge at 15,000 g at 4°C for 20 min and washed once with 70% ethanol. Finally, cDNA was reconstituted in 50 μL DEPC water and frozen at -20°C.
Primers (Table 1) used for the identification of splice variants were designed so that they bind to sequences before or after the predicted splice borders or only to the spliced sequence itself. As a result we would expect two reaction products which are different in length if a splice event takes place. All primers were produced by BioTez (Berlin). The sequences for the primer pairs were drafted using the Primer Express 2.0 software (Applied Biosystems).
RT-PCR was done mainly according to Petretti et al. . The optimal concentrations, temperatures and reaction times were determined experimentally for each individual PCR. The target sequences were amplified in 25 μl reaction volume, containing: cDNA, 2 μM dNTPs (Amersham-Pharmacia Bioteches), 50 μM of each primer (BioTez, Berlin), PCR-Buffer with 1.5 mM MgCl2 (Applied Biosystems), AmpliTaq DNA polymerase (Applied Biosystems) and DEPC-treated water (Ambion). The PCR conditions vary for each individual gene (Table 2). Reaction products obtained after 30 cycles were electrophoretically separated in 2 % agarose containing ethidium-bromide. After electrophoresis the gels were evaluated with the gel documentation system E.A.S.Y. Win32 (Herolab). All amplicons obtained were sequenced and compared with published data. For extraction of the interesting RT-PCR reaction products from the agarose gel, Qiaex II Kit (Qiagen) or Turbo Nucleic Acid Purification Kit (Geneclean) was used. Sequence determination was done by SEQLAB (Goettingen, Germany).
PCRs were carried out for each splice (normal/alternative) position with the cDNA of six cell lines (colonic: LS-174T, SW-480, HCT-116, HT-29; mammary: MDA-MB-435, MCF-7), and with RNA extracted from surgical specimens – normal mucosa and cancer tissue of 10 patients – and two samples of normal mucosa from healthy individuals which underwent a coloscopy.
Quantitative real-time PCR for Axin-1
For the TaqMan assays, 5 μL cDNA extracted from patient samples or LS174T cells, 100 nM T1 primers (normal axin splicing form), 200 nM T2 primers (alternative axin splicing isoform) and 150 nM β-Actin primers were used (primer sequences cf Table 1).
After mixing with the appropriate volume of TaqMan Universal SYBR®GREEN PCR Master Mix, quantitative real-time PCR was run in a MicroAmp Optical 96-Well Reaction Plate using the ABI 7000 Sequence Detection System. Thermal cycle conditions were as follows: 95°C for 10 min, initially, then 40 cycles of 95°C for 15 s and 62°C for 1 min. Dissociation curves were generated after amplification. Data was collected and analysed with ABI PRISM 7000 Sequence Detection System thermal cycler and software system (Applied Biosystems) according to comparative CT Method in the User Bulletin.
In order to analyse expression data according to the ΔCt method , Ct values were exported from the ABI Prism 7000 SDS Software into Microsoft Excel (Seattle, WA). Ct values of a cDNA stock from the cell line LS174T were used as calibrator. Thereby gene expression in a sample under investigation is reported as a multiple of cell culture expression. This is achieved by employing the equation (1) wherein ΔCt stands for the difference between Ct values of gene of interest and housekeeper β-Actin.
Relative Amount = Equation (1)
The SPSS software package version 10 was used for all statistical calculations. The Wilcoxon-signed-rank test was applied for calculating the difference between mRNA-expression in tumour and in mucosa of cancer patients. Associations between mRNA-expression and established histopathological parameters were calculated using Mann-Whitney- and Kruskal-Wallis-tests. A significance level of p < 0.05 was used throughout the analyses.
Dralyuk I, Brudno M, Gelfand MS, Zorn M, Dubchak I: ASDB: database of alternatively spliced genes. Nucleic Acids Res. 2000, 28: 296-297. 10.1093/nar/28.1.296.
Thanaraj TA, Stamm S, Clark F, Riethoven JJ, Le Texier V, Muilu J: ASD: the Alternative Splicing Database. Nucleic Acids Res. 2004, 32: D64-D69. 10.1093/nar/gkh030.
Kan Z, Rouchka EC, Gish WR, States DJ: Gene structure prediction and alternative splicing analysis using genomically aligned ESTs. Genome Res. 2001, 11: 889-900. 10.1101/gr.155001.
Xu Q, Modrek B, Lee C: Genome-wide detection of tissue-specific alternative splicing in the human transcriptome. Nucleic Acids Res. 2002, 30: 3754-3766. 10.1093/nar/gkf492.
Krause A, Haas SA, Coward E, Vingron M: SYSTERS, GeneNest, SpliceNest: exploring sequence space from genome to protein. Nucleic Acids Res. 2002, 30: 299-300. 10.1093/nar/30.1.299.
Lee C, Atanelov L, Modrek B, Xing Y: ASAP: the Alternative Splicing Annotation Project. Nucleic Acids Res. 2003, 31: 101-105. 10.1093/nar/gkg029.
Pospisil H, Herrmann A, Bortfeldt RH, Reich JG: EASED: Extended Alternatively Spliced EST Database. Nucleic Acids Res. 2004, 32: D70-D74. 10.1093/nar/gkh136.
Roberts GC, Smith CW: Alternative splicing: combinatorial output from the genome. Curr Opin Chem Biol. 2002, 6: 375-383. 10.1016/S1367-5931(02)00320-4.
Modrek B, Lee C: A genomic view of alternative splicing. Nat Genet. 2002, 30: 13-19. 10.1038/ng0102-13.
Venables JP: Aberrant and alternative splicing in cancer. Cancer Res. 2004, 64: 7647-7654. 10.1158/0008-5472.CAN-04-1910.
Lustig B, Behrens J: The Wnt signaling pathway and its role in tumor development. J Cancer Res Clin Oncol. 2003, 129: 199-221.
Thorstensen L, Lothe RA: The WNT Signaling Pathway and Its Role in Human Solid Tumors. Atlas Genet Cytogenet Oncol Haematol. 2003, [http://www.infobiogen.fr/services/chromcancer/Deep/WNTSignPathID20042.html]
Brinkman BM: Splice variants as cancer biomarkers. Clin Biochem. 2004, 37: 584-594. 10.1016/j.clinbiochem.2004.05.015.
Gupta S, Zink D, Korn B, Vingron M, Haas SA: Strengths and weaknesses of EST-based prediction of tissue-specific alternative splicing. BMC Genomics. 2004, 5: 72-10.1186/1471-2164-5-72.
Kikuchi A: Tumor formation by genetic mutations in the components of the Wnt signaling pathway. Cancer Sci. 2003, 94: 225-229. 10.1111/j.1349-7006.2003.tb01424.x.
Cong F, Zhang J, Pao W, Zhou P, Varmus H: A protein knowndown strategy to study the function of beta-catenin in tumorgenesis. BMC Mol Bio. 2003, 4: 10-10.1186/1471-2199-4-10.
Hlubek F, Spaderna S, Jung A, Kirchner T, Brabletz T: Beta-catenin activates a coordinated expression of the proinvasive factors laminin-5 gamma2 chain and MT1-MMP in colorectal carcinomas. Int J Cancer. 2004, 108: 321-326. 10.1002/ijc.11522.
Hulsken J, Birchmeier W, Behrens J: E-cadherin and APC compete for the interaction with beta-catenin and the cytoskeleton. J Cell Biol. 1994, 127: 2061-2069. 10.1083/jcb.127.6.2061.
Nollet F, Berx G, Molemans F, van Roy F: Genomic organization of the human beta-catenin gene (CTNNB1). Genomics. 1996, 32: 413-424. 10.1006/geno.1996.0136.
Valenta T, Lukas J, Korinek V: HMG box transcription factor TCF-4's interaction with CtBP1 controls the expression of the Wnt target Axin2/Conductin in human embryonic kidney cells. Nucleic Acids Res. 2003, 31: 2369-2380. 10.1093/nar/gkg346.
Grooteclaes M, Deveraux Q, Hildebrand J, Zhang Q, Goodman RH, Frisch SM: C-terminal-binding protein corepresses epithelial and proapoptotic gene expression programs. PNAS. 2003, 100: 4568-4573. 10.1073/pnas.0830998100.
Sparks AB, Morin PJ, Vogelstein B, Kinzler KW: Mutational analysis of the APC/beta-catenin/Tcf pathway in colorectal cancer. Cancer Res. 1998, 58: 1130-1134.
Mukai F, Ishiguro K, Sano Y, Fujita SC: Alternative splicing isoform of tau protein kinase I/glycogen synthase kinase 3beta. J Neurochem. 2002, 81: 1073-1083. 10.1046/j.1471-4159.2002.00918.x.
Mao J, Wang J, Liu B, Pan W, Farr GH, Flynn C, Yuan H, Takada S, Kimelman D, Li l, Wu D: Low-Density Lipoprotein Receptor-Related Protein-5 Binds to Axin and Regulates the Canonical Wnt Signaling Pathway. Molecular Cell. 2001, 7: 801-809. 10.1016/S1097-2765(01)00224-6.
Strickland DK, Ranganathan S: Diverse role of LDL receptor-related protein in the clearance of proteases and in signalling. J Thromb Haemost. 2003, 1: 1663-1670. 10.1046/j.1538-7836.2003.00330.x.
Hoang BH, Kubo T, Healey JH, Sowers R, Mazza B, Yang R, Huvos AG, Meyers PA, Gorlick R: Expression of LDL receptor-related protein 5 (LRP5) as a novel marker for disease progression in high-grade osteosarcoma. Int J Cancer. 2004, 109: 106-111. 10.1002/ijc.11677.
Fujino T, Asaba H, Kang MJ, Ikeda Y, Sone H, Takada S, Kim DH, Ioka RX, Ono M, Tomoyori H, Okubo M, Murase T, Kamataki A, Yamamoto J, Magoori K, Takahashi S, Miyamoto Y, Oishi H, Nose M, Okazaki M, Usui S, Imaizumi K, Yanagisawa M, Sakai J, Yamamoto TT: Low-density lipoprotein receptor-related protein 5 (LRP5) is essential for normal cholesterol metabolism and glucose-induced insulin secretion. PNAS. 2003, 100: 229-234. 10.1073/pnas.0133792100.
Okubo M, Horinishi A, Kim DH, Yamamoto TT, Murase T: Seven novel sequence variants in the human low density lipoprotein receptor related protein 5 (LRP5) gene. Hum Mut. 2002, 19: 186-10.1002/humu.9012.
Hussain MM: Structural, biochemical and signaling properties of the low-density lipoprotein receptor gene family. Front Biosci. 2001, 6: D417-D428.
Kishida S, Yamamoto H, Ikeda S, Kishida M, Sakamoto I, Koyama S, Kikuchi A: Axin, a negative regulator of the Wnt signaling pathway, directly interacts with adenomatous polyposis coli and regulates the stabilization of β-catenin. J Biol Chem. 1998, 273: 10823-10826. 10.1074/jbc.273.18.10823.
Satoh S, Daigo Y, Furukawa Y, Kato T, Miwa N, Nishiwaki T, Kawasoe T, Ishiguro H, Fujita M, Tokino T, Sasaki Y, Imaoka S, Murata M, Shimano T, Yamaoka Y, Nakamura Y: AXIN1 mutations in hepatocellular carcinomas, and growth suppression in cancer cells by virus-mediated transfer of AXIN1. Nat Genet. 2000, 24: 245-250. 10.1038/73448.
Wong CK, Luo W, Deng Y, Zou H, Ye Z, Lin SC: The DIX domain protein coiled-coil-DIX1 inhibits c-Jun N-terminal kinase activation by Axin and dishevelled through distinct mechanisms. J Biol Chem. 2004, 279: 39366-39373. 10.1074/jbc.M404598200.
Nakajima M, Fukuchi M, Miyazaki T, Masuda N, Kato H, Kuwano H: Reduced expression of Axin correlates with tumour progression of oesophageal squamous cell carcinoma. Br J Cancer. 2003, 88: 1734-1739. 10.1038/sj.bjc.6600941.
Luo W, Zou H, Jin L, Lin S, Li Q, Ye Z, Rui H, Lin SC: Axin contains three separable domains that confer intramolecular, homodimeric, and heterodimeric interactions involved in distinct functions. J Biol Chem. 2005, 280: 5054-5060. 10.1074/jbc.M412340200.
Birney E, Andrews D, Bevan P, Caccamo M, Cameron G, Chen Y, Clarke L, Coates G, Cox T, Cuff J, Curwen V, Cutts T, Down T, Durbin R, Eyras E, Fernandez-Suarez XM, Gane P, Gibbins B, Gilbert J, Hammond M, Hotz H, Iyer V, Kahari A, Jekosch K, Kasprzyk A, Keefe D, Keenan S, Lehvaslaiho H, McVicker G, Melsopp C, Meidl P, Mongin E, Pettett R, Potter S, Proctor G, Rae M, Searle S, Slater G, Smedley D, Smith J, Spooner W, Stabenau A, Stalker J, Storey R, Ureta-Vidal A, Woodwark C, Clamp M, Hubbard T: Ensembl 2004. Nucleic Acids Res. 2004, 32: D468-D470. 10.1093/nar/gkh038.
Benson DA, Karsch-Mizrachi I, Lipman DJ, Ostell J, Wheeler DL: GenBank: update. Nucl Acids Res. 2004, 32: D23-D26. 10.1093/nar/gkh045.
Gish W: 1996, [http://blast.wustl.edu]
Smit AFA, Hubley R, Green P: RepeatMasker Open-3.0. 2004, [http://www.repeatmasker.org]
Petretti T, Kemmner W, Schulze B, Schlag PM: Altered mRNA expression of glycosyltransferases in human colorectal carcinomas and liver metastases. GUT. 2000, 46: 359-366. 10.1136/gut.46.3.359.
Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262.
This work was partially supported from the Bundesministerium für Forschung und Bildung, Germany, through its contribution to the Helmholtz Network for Bioinformatics, and by the ATD project (funded by the European Commission within its FP6 Programme, LHSG-CT-2003-503329). We thank Ian F. Korf, Joseph A. Bedell, and Warren R. Gish for providing MaskerAid and WU Blast. The authors thank the reviewers for their useful comments.
AH wrote the prediction software and established the database EASED. SP participated in the selection of the Wnt components. HP and JR provided guidance for the computational work. KB performed the RT-PCR experiments under the guidance of WK. All Authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.