In silico screening for tumour-specific expressed sequences in human genome

Carcino-Evo-Devo, A Theory of the Evolutionary Role of Hereditary Tumors

A theory of the evolutionary role of hereditary tumors, or the carcino-evo-devo theory, is being developed. The main hypothesis of the theory, the hypothesis of evolution by tumor neofunctionalization, posits that hereditary tumors provided additional cell masses during the evolution of multicellular organisms for the expression of evolutionarily novel genes. The carcino-evo-devo theory has formulated several nontrivial predictions that have been confirmed in the laboratory of the author. It also suggests several nontrivial explanations of biological phenomena previously unexplained by the existing theories or incompletely understood. By considering three major types of biological development-individual, evolutionary, and neoplastic development-within one theoretical framework, the carcino-evo-devo theory has the potential to become a unifying biological theory.

lncRNA ELFN1-AS1 promotes proliferation, migration and invasion and suppresses apoptosis in colorectal cancer cells by enhancing G6PD activity

Tumour cells change their metabolic patterns to support high proliferation rates and cope with oxidative stress. The lncRNA ELFN1-AS1 is highly expressed in a wide range of cancers and is essential to the proliferation and apoptosis of tumour cells. Nevertheless, its function in the metabolic reprogramming of tumour cells is unclear. Here we show that ELFN1-AS1 promotes glucose consumption as well as lactate and NADPH production. Database searching, bioinformatics analysis, RNA immunoprecipitation (RIP) and RNA pull-down assays show that ELFN1-AS1 enhances glucose-6-phosphate dehydrogenase ( G6PD) expression and activates the pentose phosphate pathway (PPP) by promoting TP53 degradation. In addition, luciferase reporter assay and chromatin immunoprecipitation (ChIP) show that YY1 binds to the ELFN1-AS1 promoter to promote transcriptional activation of ELFN1-AS1. Consistent with the in vitro experiments, knockdown of ELFN1-AS1 impedes the growth of tumours transplanted into mice by inhibiting the expression of G6PD. In conclusion, this study reveals that ELFN1-AS1 activates the PPP, and validates the regulatory role of the YY1/ ELFN1-AS1/ TP53/ G6PD axis in colorectal cancer.

The Theory of Carcino-Evo-Devo and Its Non-Trivial Predictions

To explain the sources of additional cell masses in the evolution of multicellular organisms, the theory of carcino-evo-devo, or evolution by tumor neofunctionalization, has been developed. The important demand for a new theory in experimental science is the capability to formulate non-trivial predictions which can be experimentally confirmed. Several non-trivial predictions were formulated using carcino-evo-devo theory, four of which are discussed in the present paper: (1) The number of cellular oncogenes should correspond to the number of cell types in the organism. The evolution of oncogenes, tumor suppressor and differentiation gene classes should proceed concurrently. (2) Evolutionarily new and evolving genes should be specifically expressed in tumors (TSEEN genes). (3) Human orthologs of fish TSEEN genes should acquire progressive functions connected with new cell types, tissues and organs. (4) Selection of tumors for new functions in the organism is possible. Evolutionarily novel organs should recapitulate tumor features in their development. As shown in this paper, these predictions have been confirmed by the laboratory of the author. Thus, we have shown that carcino-evo-devo theory has predictive power, fulfilling a fundamental requirement for a new theory.

ELFN1-AS1 accelerates the proliferation and migration of colorectal cancer via regulation of miR-4644/TRIM44 axis

Faced with the increasing colorectal cancer (CRC) cases, the interrogation of pivotal molecules in CRC appears to be vitally important. Long non-coding RNAs (lncRNAs) are well-known regulators of gene expression at transcriptional, post-transcriptional or epigenetic level, among which the competing endogenous RNA (ceRNA) network is a common way that lncRNAs exert their properties. The current study aimed to provide a new insight into improving the outcomes of CRC patients. Our study detected that ELFN1-AS1 expression was elevated in CRC tissues and cells, and ELFN1-AS1 upregulation was correlated with poor prognosis of CRC sufferers. Besides, it was viewed that ELFN1-AS1 knockdown impeded the proliferation and migration abilities as well as activated the apoptosis ability of CRC cells. In subsequence, mechanism assays also displayed that ELFN1-AS1 targeted miR-4644 to augment TRIM44 level. Finally, rescue experiments confirmed that TRIM44 took part in the ELFN1-AS1-medatied promotional influences on CRC cells proliferation and migration. In conclusion, ELFN1-AS1 exerted pro-proliferation, anti-apoptosis and pro-migration functions on CRC cells by acting as a sponge of miR-4644 to increase TRIM44 expression at mRNA and protein level, providing an additional molecule responsible for the carcinogenesis and progression for CRC.

Evolutionarily novel genes are expressed in transgenic fish tumors and their orthologs are involved in development of progressive traits in humans

Abstract

Earlier we suggested a new hypothesis of the possible evolutionary role of hereditary tumors (Kozlov, Evolution by tumor Neofunctionalization, 2014), and described a new class of genes – tumor specifically expressed, evolutionarily novel (TSEEN) genes - that are predicted by this hypothesis (Kozlov, Infect Agents Cancer 11:34, 2016). In this paper we studied evolutionarily novel genes expressed in fish tumors after regression, as a model of evolving organs. As evolutionarily novel genes may not yet have organismal functions, we studied the acquisition of new gene functions by comparing fish evolutionarily novel genes with their human orthologs. We found that many genes involved in development of progressive traits in humans (lung, mammary gland, placenta, ventricular septum, etc.) originated in fish and are expressed in fish tumors and tumors after regression. These findings support a possible evolutionary role of hereditary tumors, and in particular the hypothesis of evolution by tumor neofunctionalization.

Research highlights

Earlier we described a new class of genes that are tumor-specifically expressed and evolutionarily novel (TSEEN). As the functions of TSEEN genes are often uncertain, we decided to study TSEEN genes of fishes so that we could trace the appearance of their new functions in higher vertebrates. We found that many human genes which are involved in development of progressive traits (placenta development, mammary gland and lung development etc.,) originated in fishes and are expressed in fish tumors.

Oncogenes, tumor suppressor and differentiation genes represent the oldest human gene classes and evolve concurrently

Earlier we showed that human genome contains many evolutionarily young or novel genes with tumor-specific or tumor-predominant expression. We suggest calling such genes Tumor Specifically Expressed, Evolutionarily New (TSEEN) genes. In this paper we performed a study of the evolutionary ages of different classes of human genes, using homology searches in genomes of different taxa in human lineage. We discovered that different classes of human genes have different evolutionary ages and confirmed the existence of TSEEN gene classes. On the other hand, we found that oncogenes, tumor-suppressor genes and differentiation genes are among the oldest gene classes in humans and their evolution occurs concurrently. These findings confirm non-trivial predictions made by our hypothesis of the possible evolutionary role of hereditary tumors. The results may be important for better understanding of tumor biology. TSEEN genes may become the best tumor markers.

Expression of evolutionarily novel genes in tumors

The evolutionarily novel genes originated through different molecular mechanisms are expressed in tumors. Sometimes the expression of evolutionarily novel genes in tumors is highly specific. Moreover positive selection of many human tumor-related genes in primate lineage suggests their involvement in the origin of new functions beneficial to organisms. It is suggested to consider the expression of evolutionarily young or novel genes in tumors as a new biological phenomenon, a phenomenon of TSEEN (tumor specifically expressed, evolutionarily novel) genes.

ELFN1-AS1: a novel primate gene with possible microRNA function expressed predominantly in human tumors

Human gene LOC100505644 uncharacterized LOC100505644 [Homo sapiens] (Entrez Gene ID 100505644) is abundantly expressed in tumors but weakly expressed in few normal tissues. Till now the function of this gene remains unknown. Here we identified the chromosomal borders of the transcribed region and the major splice form of the LOC100505644-specific transcript. We characterised the major regulatory motifs of the gene and its splice sites. Analysis of the secondary structure of the major transcript variant revealed a hairpin-like structure characteristic for precursor microRNAs. Comparative genomic analysis of the locus showed that it originated in primates de novo. Taken together, our data indicate that human gene LOC100505644 encodes some non-protein coding RNA, likely a microRNA. It was assigned a gene symbol ELFN1-AS1 (ELFN1 antisense RNA 1 (non-protein coding)). This gene combines features of evolutionary novelty and predominant expression in tumors.

Binding properties and dynamic localization of an alternative isoform of the cap-binding complex subunit CBP20

The nuclear cap-binding complex (CBC) is a heterodimer composed of CBP20 and CBP80 subunits and has roles in the biogenesis of messenger RNAs (mRNAs), small nuclear RNAs (snRNAs) and microRNAs. CBP20 is a phylogenetically conserved protein that interacts with the 7-methyl guanosine (m7G) cap added to the 5' end of all RNA polymerase II transcripts. CBP80 ensures high affinity binding of the cap by CBP20 and provides a platform for interactions with other factors. Here we characterize an alternative splice variant of CBP20, termed CBP20S. The CBP20S transcript has an in-frame deletion, leading to the translation of a protein lacking most of the RNA recognition motif (RRM). We show that CBP20S is conserved among mammalian species and is expressed in human cell lines and bone marrow cells. Unlike the full-length CBP20, CBP20S does not bind CBP80 or the m7G cap. Nevertheless, CBP20S does bind mRNA, is localized to an active transcription site and redistributed to nucleolar caps upon transcription inhibition. Our results suggest that this novel form CBP20S plays a role in transcription and/or RNA processing independent of CBP80 or the cap.

IN SILICO SEARCH FOR NATURAL ANTISENSE TRANSCRIPTS REVEALS THEIR DIFFERENTIAL EXPRESSION IN HUMAN TUMORS

We created an algorithm that allows high-throughput mapping of sense-antisense (SA) pairs of transcripts. By this method we mapped approximately 32 000 SA pairs of human mRNAs. Collected SA pairs were divided into three groups: SA pairs based on two or more UniGene clusters (17% of all sense-antisense pairs), SA pairs based on ESTs that belong to the same UniGene cluster (42%), and SA pairs formed by UniGene cluster and non-unique unclustered transcripts (41%). To study expression patterns of natural SA pairs we created a software application "Antisense Cluster Filter". We retrieved tissue expression data for all the transcripts forming identified SA pairs, including clustered and unclustered ones. After that, we selected 108 SA pairs represented by transcripts differentially regulated in human tumors. For each of these SA pairs one of the transcripts was expressed only in tumors, another one was expressed both in non-malignant and malignant tissues. Indicated SA pairs may represent a new class of tumor markers. An example of the tumor-specific natural antisense to C3orf4 mRNA is detailed.

Quantitative gene expression profiles in real time from expressed sequence tag databases

An accumulation of expressed sequence tag (EST) data in the public domain and the availability of bioinformatic programs have made EST gene expression profiling a common practice. However, the utility and validity of using EST databases (e.g., dbEST) has been criticized, particularly for quantitative assessment of gene expression. Problems with EST sequencing errors, library construction, EST annotation, and multiple paralogs make generation of specific and sensitive qualitative arid quantitative expression profiles a concern. In addition, most EST-derived expression data exists in previously assembled databases. The Virtual Northern Blot (VNB) (http: //tlab.bu.edu/vnb.html) allows generation, evaluation, and optimization of expression profiles in real time, which is especially important for alternatively spliced, novel, or poorly characterized genes. Representative gene families with variable nucleotide sequence identity, tissue specificity, and levels of expression (bcl-xl, aldoA, and cyp2d9) are used to assess the quality of VNB's output. The profiles generated by VNB are more sensitive and specific than those constructed with ESTs listed in preindexed databases at UCSC and NCBI. Moreover, quantitative expression profiles produced by VNB are comparable to quantization obtained from Northern blots and qPCR. The VNB pipeline generates real-time gene expression profiles for single-gene queries that are both qualitatively and quantitatively reliable.

The possible evolutionary role of tumors in the origin of new cell types

The ability of tumor cells to differentiate in combination with their ability to express genes that are not expressed in normal tissues, may result in the emergence of new cell types in evolution. Tumors may play an evolutionary role by providing conditions (space and resources) for the expression of newly evolving genes. Genetically or epigenetically predetermined tumors at the early stages of progression, benign tumors, and some tumor-like processes in invertebrates and plants, all of which are modes of excess cell growth which provide evolving multicellular organisms with extra cell masses, are considered as potentially evolutionarily meaningful. Malignant tumors at the late stages of progression, however, are not. The preexisting cell types of multicellular organisms had restricted potential for the expression of newly evolving genes. Because of regulation and gene competition, some of the newly evolving genes may stay silent. Multicellular organisms would need excess cell masses for the expression of newly evolving genes. The preexisting cell types cannot provide such excess cell masses because of limitations imposed on the number of possible cell divisions. Tumors could provide the evolving multicellular organisms with the excess cell masses for the expression of newly evolving genes. We suggest that tumors could be a sort of proving ground (or reservoir) for the expression of newly evolving genes that originate in the course of genome evolution in the DNA of germ cells (i.e., not in tumor cells themselves). The case in which the expression of a newly evolving gene in tumors results in the origin of a new function would be associated with the origin of new feedback and regulatory circuits, as in root nodules in legumes and macromelanophores in Xiphophorus fishes. Tumor cells would differentiate, resulting in a new cell type for the given multicellular species. This cell type would be inherited because of epigenomic mechanisms similar to those in preexisting cell types. Populations of tumor-bearing organisms with genetically or epigenetically programmed tumors could represent the transition between established species of organisms at different stages of progressive evolution. Experimental confirmation of the prediction of the hypothesis of evolution by tumor cells differentiation concerning the expression of evolutionarily new genes and/or silent (neutrally evolving) sequences in tumor cells is presented.

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.

The human T-box mesodermal transcription factor Brachyury is a candidate target for T-cell-mediated cancer immunotherapy

PURPOSE

Identification of tumor antigens is essential in advancing immune-based therapeutic interventions in cancer. Particularly attractive targets are those molecules that are selectively expressed by malignant cells and that are also essential for tumor progression.

EXPERIMENTAL DESIGN AND RESULTS

We have used a computer-based differential display analysis tool for mining of expressed sequence tag clusters in the human Unigene database and identified Brachyury as a novel tumor antigen. Brachyury, a member of the T-box transcription factor family, is a key player in mesoderm specification during embryonic development. Moreover, transcription factors that control mesoderm have been implicated in the epithelial-mesenchymal transition (EMT), which has been postulated to be a key step during tumor progression to metastasis. Reverse transcription-PCR analysis validated the in silico predictions and showed Brachyury expression in tumors of the small intestine, stomach, kidney, bladder, uterus, ovary, and testis, as well as in cell lines derived from lung, colon, and prostate carcinomas, but not in the vast majority of the normal tissues tested. An HLA-A0201 epitope of human Brachyury was identified that was able to expand T lymphocytes from blood of cancer patients and normal donors with the ability to lyse Brachyury-expressing tumor cells.

CONCLUSIONS

To our knowledge, this is the first demonstration that (a) a T-box transcription factor and (b) a molecule implicated in mesodermal development, i.e., EMT, can be a potential target for human T-cell-mediated cancer immunotherapy.

Functional and evolutionary analysis of alternatively spliced genes is consistent with an early eukaryotic origin of alternative splicing

Background

Alternative splicing has been reported in various eukaryotic groups including plants, apicomplexans, diatoms, amoebae, animals and fungi. However, whether widespread alternative splicing has evolved independently in the different eukaryotic groups or was inherited from their last common ancestor, and may therefore predate multicellularity, is still unknown. To better understand the origin and evolution of alternative splicing and its usage in diverse organisms, we studied alternative splicing in 12 eukaryotic species, comparing rates of alternative splicing across genes of different functional classes, cellular locations, intron/exon structures and evolutionary origins.

Results

For each species, we find that genes from most functional categories are alternatively spliced. Ancient genes (shared between animals, fungi and plants) show high levels of alternative splicing. Genes with products expressed in the nucleus or plasma membrane are generally more alternatively spliced while those expressed in extracellular location show less alternative splicing. We find a clear correspondence between incidence of alternative splicing and intron number per gene both within and between genomes. In general, we find several similarities in patterns of alternative splicing across these diverse eukaryotes.

Conclusion

Along with previous studies indicating intron-rich genes with weak intron boundary consensus and complex spliceosomes in ancestral organisms, our results suggest that at least a simple form of alternative splicing may already have been present in the unicellular ancestor of plants, fungi and animals. A role for alternative splicing in the evolution of multicellularity then would largely have arisen by co-opting the preexisting process.

High qualitative and quantitative conservation of alternative splicing in Caenorhabditis elegans and Caenorhabditis briggsae

Alternative splicing (AS) is an important contributor to proteome diversity and is regarded as an explanatory factor for the relatively low number of human genes compared with less complex animals. To assess the evolutionary conservation of AS and its developmental regulation, we have investigated the qualitative and quantitative expression of 21 orthologous alternative splice events through the development of 2 nematode species separated by 85-110 Myr of evolutionary time. We demonstrate that most of these alternative splice events present in Caenorhabditis elegans are conserved in Caenorhabditis briggsae. Moreover, we find that relative isoform expression levels vary significantly during development for 78% of the AS events and that this quantitative variation is highly conserved between the 2 species. Our results suggest that AS is generally tightly regulated through development and that the regulatory mechanisms controlling AS are to a large extent conserved during the evolution of Caenorhabditis. This strong conservation indicates that both major and minor splice forms have important functional roles and that the relative quantities in which they are expressed are crucial. Our results therefore suggest that the quantitative regulation of isoform expression levels is an intrinsic part of most AS events. Moreover, our results indicate that AS contributes little to transcript variation in Caenorhabditis genes and that gene duplication may be the major evolutionary mechanism for the origin of novel transcripts in these 2 species.

In silico whole-genome screening for cancer-related single-nucleotide polymorphisms located in human mRNA untranslated regions

BACKGROUND

A promising application of the huge amounts of genetic data currently available lies in developing a better understanding of complex diseases, such as cancer. Analysis of publicly available databases can help identify potential candidates for genes or mutations specifically related to the cancer phenotype. In spite of their huge potential to affect gene function, no systematic attention has been paid so far to the changes that occur in untranslated regions of mRNA.

RESULTS

In this study, we used Expressed Sequence Tag (EST) databases as a source for cancer-related sequence polymorphism discovery at the whole-genome level. Using a novel computational procedure, we focused on the identification of untranslated region (UTR)-localized non-coding Single Nucleotide Polymorphisms (UTR-SNPs) significantly associated with the tumoral state. To explore possible relationships between genetic mutation and phenotypic variation, bioinformatic tools were used to predict the potential impact of cancer-associated UTR-SNPs on mRNA secondary structure and UTR regulatory elements. We provide a comprehensive and unbiased description of cancer-associated UTR-SNPs that may be useful to define genotypic markers or to propose polymorphisms that can act to alter gene expression levels. Our results suggest that a fraction of cancer-associated UTR-SNPs may have functional consequences on mRNA stability and/or expression.

CONCLUSION

We have undertaken a comprehensive effort to identify cancer-associated polymorphisms in untranslated regions of mRNA and to characterize putative functional UTR-SNPs. Alteration of translational control can change the expression of genes in tumor cells, causing an increase or decrease in the concentration of specific proteins. Through the description of testable candidates and the experimental validation of a number of UTR-SNPs discovered on the secreted protein acidic and rich in cysteine (SPARC) gene, this report illustrates the utility of a cross-talk between in silico transcriptomics and cancer genetics.

Evolutionarily new sequences expressed in tumors

Background

Earlier we suggested the concept of the positive evolutionary role of tumors. According to this concept, tumors provide conditions for the expression of evolutionarily new and/or sleeping genes in their cells. Thus, tumors are considered as evolutionary proving ground or reservoir of expression. To support this concept we have previously characterized in silico and experimentally a new class of human tumor-related transcribed sequences.

Results

In this article we describe results of further studies of previously described tumor-related sequences. The results of molecular phylogeny studies, Southern hybridization experiments and computational comparison with genomes of other species are presented.

Conclusion

These results suggest that these previously described tumor-related human transcripts are also relatively evolutionarily new.

Maximum likelihood inference of imprinting and allele-specific expression from EST data

MOTIVATION

In a diploid organism the proportion of transcripts that are produced from the two parental alleles can differ substantially due, for example to epigenetic modification that causes complete or partial silencing of one parental allele or to cis acting polymorphisms that affect transcriptional regulation. Counts of SNP alleles derived from EST sequences have been used to identify both novel candidates for genomic imprinting as well as examples of genes with allelic differences in expression.

RESULTS

We have developed a set of statistical models in a maximum likelihood framework that can make highly efficient use of public transcript data to identify genes with unequal representation of alternative alleles in cDNA libraries. We modelled both imprinting and allele-specific expression and applied the models to a large dataset of SNPs mapped to EST sequences. Using simulations, matched closely to real data, we demonstrate significantly improved performance over existing methods that have been applied to the same data. We further validated the power of this approach to detect imprinting using a set of known imprinted genes and inferred a set of candidate imprinted genes, several of which are in close proximity to known imprinted genes. We report evidence that there are undiscovered imprinted genes in known imprinted regions. Overall, more than half of the genes for which the most data are available show some evidence of allele-specific expression.

AVAILABILITY

Software is available from the authors on request.

Template switching within exons 3 and 4 of KV11.1 (HERG) gives rise to a 5' truncated cDNA

K(V)11.1 (HERG) channels contribute to membrane potential in a number of excitable cell types. We cloned a variant of K(V)11.1 from human jejunum containing a 171 bp deletion spanning exons 3 and 4. Expression of a full-length cDNA clone containing this deletion gave rise to protein that trafficked to the cell membrane and generated robust currents. The deletion occurred in a G/C-rich region and identical sequence elements of UGGUGG were located at the deletion boundaries. In recent studies these features have been implicated to cause deletions via template switching during cDNA synthesis. To examine this possibility we compared cDNAs from human brain, heart, and jejunum synthesized at lower (42 degrees C) and higher temperatures (70 degrees C). The 171 bp deletion was absent at the higher temperature. Our results suggest that the sequence and secondary structure of mRNA in the G/C rich region leads to template switching producing a cDNA product with a 171 bp deletion.

Bioinformatic screening of human ESTs for differentially expressed genes in normal and tumor tissues

Background

Owing to the explosion of information generated by human genomics, analysis of publicly available databases can help identify potential candidate genes relevant to the cancerous phenotype. The aim of this study was to scan for such genes by whole-genome in silico subtraction using Expressed Sequence Tag (EST) data.

Methods

Genes differentially expressed in normal versus tumor tissues were identified using a computer-based differential display strategy. Bcl-xL, an anti-apoptotic member of the Bcl-2 family, was selected for confirmation by western blot analysis.

Results

Our genome-wide expression analysis identified a set of genes whose differential expression may be attributed to the genetic alterations associated with tumor formation and malignant growth. We propose complete lists of genes that may serve as targets for projects seeking novel candidates for cancer diagnosis and therapy. Our validation result showed increased protein levels of Bcl-xL in two different liver cancer specimens compared to normal liver. Notably, our EST-based data mining procedure indicated that most of the changes in gene expression observed in cancer cells corresponded to gene inactivation patterns. Chromosomes and chromosomal regions most frequently associated with aberrant expression changes in cancer libraries were also determined.

Conclusion

Through the description of several candidates (including genes encoding extracellular matrix and ribosomal components, cytoskeletal proteins, apoptotic regulators, and novel tissue-specific biomarkers), our study illustrates the utility of in silico transcriptomics to identify tumor cell signatures, tumor-related genes and chromosomal regions frequently associated with aberrant expression in cancer.

Computational and experimental analysis identifies Arabidopsis genes specifically expressed during early seed development

Background

Plant seeds are complex organs in which maternal tissues, embryo and endosperm, follow distinct but coordinated developmental programs. Some morphogenetic and metabolic processes are exclusively associated with seed development. The goal of this study was to explore the feasibility of incorporating the available online bioinformatics databases to discover Arabidopsis genes specifically expressed in certain organs, in our case immature seeds.

Results

A total of 11,032 EST sequences obtained from isolated immature seeds were used as the initial dataset (178 of them newly described here). A pilot study was performed using EST virtual subtraction followed by microarray data analysis, using the Genevestigator tool. These techniques led to the identification of 49 immature seed-specific genes. The findings were validated by RT-PCR analysis and in situ hybridization.

Conclusion

We conclude that the combined in silico data analysis is an effective data mining strategy for the identification of tissue-specific gene expression.

A genome-wide survey demonstrates widespread non-linear mRNA in expressed sequences from multiple species

We describe here the results of the first genome-wide survey of candidate exon repetition events in expressed sequences from human, mouse, rat, chicken, zebrafish and fly. Exon repetition is a rare event, reported in <10 genes, in which one or more exons is tandemly duplicated in mRNA but not in the gene. To identify candidates, we analysed database sequences for mRNA transcripts in which the order of the spliced exons does not follow the linear genomic order of the individual gene [events we term rearrangements or repetition in exon order (RREO)]. Using a computational approach, we have identified 245 genes in mammals that produce RREO events. RREO in mRNA occurs predominantly in the coding regions of genes. However, exon 1 is never involved. Analysis of the open reading frames suggests that this process may increase protein diversity and regulate protein expression via nonsense-mediated RNA decay. The sizes of the exons and introns involved around these events suggest a gene model structure that may facilitate non-linear splicing. These findings imply that RREO affects a significant subset of genes within a genome and suggests that non-linear information encoded within the genomes of complex organisms could contribute to phenotypic variation.

The importance of being divisible by three in alternative splicing

Alternative splicing events that are conserved in orthologous genes in different species are commonly viewed as reliable evidence of authentic, functionally significant alternative splicing events. Several recent bioinformatic analyses have shown that conserved alternative exons possess several features that distinguish them from alternative exons that are species-specific. One of the most striking differences between conserved and species-specific alternative exons is the high percentage of exons that preserve the reading frame (exons whose length is an exact multiple of 3, termed symmetrical exons) among the conserved alternative exons. Here, we examined conserved alternative exons and found several features that differentiate between symmetrical and non-symmetrical alternative exons. We show that symmetrical alternative exons have a strong tendency not to disrupt protein domain structures, whereas the tendency of non-symmetrical alternative exons to overlap with different fractions of protein domains is similar to that of constitutive exons. Additionally, skipping isoforms of non-symmetrical alternative exons are strongly underrepresented, compared with their including isoforms, suggesting that skipping of a large fraction of non-symmetrical alternative exons produces transcripts that are degraded by the nonsense-mediated mRNA decay mechanism. Non-symmetrical alternative exons also show a tendency to reside in the 5' half of the CDS. These findings suggest that alternative splicing of symmetrical and non-symmetrical exons is governed by different selective pressures and serves different purposes.

Sequence conservation, relative isoform frequencies, and nonsense-mediated decay in evolutionarily conserved alternative splicing

Studies of expressed sequence tag data sets have revealed large numbers of splicing variants for human genes, but it remains challenging to distinguish functionally important variants from aberrant splicing, clarify the nature of the alternative functions, and understand the signals that regulate splicing choices. To help address these issues, we have constructed and analyzed a large data set of 1,478 exon-skipping alternative splicing (AS) variants evolutionarily conserved in human and mouse. In about one-fifth of cases, one isoform appears subject to nonsense-mediated mRNA decay (NMD), supporting the idea that a major role of AS is to regulate gene expression; one-quarter of these NMD-inducing cases involve a conserved exon whose apparent sole purpose is to mediate destruction of the message when included. We explore sequence conservation likely related to splicing regulation, using in part a measure of the overall amount of conserved information in a sequence, and find that the increased conservation that has been observed within AS exons primarily affects synonymous sites, suggesting that regulatory signals significantly constrain synonymous substitution rates. We show that a lower frequency of the inclusion isoform relative to the exclusion isoform tends to be associated with weaker splice site signals, smaller exon size, and higher intronic sequence conservation, and provide evidence that all of these factors are under selection to control relative isoform frequencies. Some conserved instances of AS appear to represent aberrant splicing events that by chance have occurred in both species, and we develop a nonparametric likelihood approach to identify these.

GBA server: EST-based digital gene expression profiling

Expressed Sequence Tag-based gene expression profiling can be used to discover functionally associated genes on a large scale. Currently available web servers and tools focus on finding differentially expressed genes in different samples or tissues rather than finding co-expressed genes. To fill this gap, we have developed a web server that implements the GBA (Guilt-by-Association) co-expression algorithm, which has been successfully used in finding disease-related genes. We have also annotated UniGene clusters with links to several important databases such as GO, KEGG, OMIM, Gene, IPI and HomoloGene. The GBA server can be accessed and downloaded at .

Experimental Study of Human Expressed Sequences Newly Identified in Silico as Tumor Specific

In our previous studies, we used global computational differential display of ESTs that belong to UNIGENE clusters and identified human sequences differentially expressed in human tumors, as well as a considerable amount of transcripts represented only in tumor-derived cDNA libraries. Most of the tumor-specific EST clusters are derived from the plurality of the tumor types originated in tissues of both ectodermal and mesodermal origin. We found that many of such tumor-specific ESTs do not contain long open reading frames and cannot be classified as protein-encoding genes. To experimentally assess patterns of expression of these EST clusters, we studied four of them in PCR experiments on Clontech MTC panels. The experimental data confirm the results obtained by in silico screening, i.e. tumor specificity of their expression. We suggest that a significant increase in the expression of non-coding RNA is a fundamental feature of cancer cells, and that such transcripts could serve as markers for the diagnosis or monitoring of human malignancies.

Identification of alternatively spliced mRNA variants related to cancers by genome-wide ESTs alignment

Several databases have been published to predict alternative splicing of mRNAs by analysing the exon linkage relationship by alignment of expressed sequence tags (ESTs) to the genome sequence; however, little effort has been made to investigate the relationship between cancers and alternative splicing. We developed a program, Alternative Splicing Assembler (ASA), to look for splicing variants of human gene transcripts by genome-wide ESTs alignment. Using ASA, we constructed the biosino alternative splicing database (BASD), which predicted splicing variants for reference sequences from the reference sequence database (RefSeq) and presented them in both graph and text formats. EST clusters that differ from the reference sequences in at least one splicing site were counted as splicing variants. Of 4322 genes screened, 3490 (81%) were observed with at least one alternative splicing variants. To discover the variants associated with cancers, tissue sources of EST sequences were extracted from the UniLib database and ESTs from the same tissue type were counted. These were regarded as the indicators for gene expression level. Using Fisher's exact test, alternative splicing variants, of which EST counts were significantly different between cancer tissues and their counterpart normal tissues, were identified. It was predicted that 2149 variants, or 383 variants after Bonferroni correction, of 26 812 variants were likely tumor-associated. By reverse transcription-PCR, 11 of 13 novel alternative splicing variants and eight of nine variants' tissue specificity were confirmed in hepatocellular carcinoma and in lung cancer. The possible involvement of alternative splicing in cancer is discussed.

How prevalent is functional alternative splicing in the human genome?

Comparative analyses of ESTs and cDNAs with genomic DNA predict a high frequency of alternative splicing in human genes. However, there is an ongoing debate as to how many of these predicted splice variants are functional and how many are the result of aberrant splicing (or 'noise'). To address this question, we compared alternatively spliced cassette exons that are conserved between human and mouse with EST-predicted cassette exons that are not conserved in the mouse genome. Presumably, conserved exon-skipping events represent functional alternative splicing. We show that conserved (functional) cassette exons possess unique characteristics in size, repeat content and in their influence on the protein. By contrast, most non-conserved cassette exons do not share these characteristics. We conclude that a significant portion of cassette exons evident in EST databases is not functional, and might result from aberrant rather than regulated splicing.