Heterogeneity in polyadenylation cleavage sites in mammalian mRNA sequences: implications for SAGE analysis

Similarities and differences of polyadenylation signals in human and fly

BACKGROUND

Cleavage of messenger RNA (mRNA) precursors is an essential step in mRNA maturation. The signal recognized by the cleavage enzyme complex has been characterized as an A rich region upstream of the cleavage site containing a motif with consensus AAUAAA, followed by a U or UG rich region downstream of the cleavage site.

RESULTS

We studied these signals using exhaustive databases of cleavage sites obtained from aligning raw expressed sequence tags (EST) sequences to genomic sequences in Homo sapiens and Drosophila melanogaster. These data show that the polyadenylation signal is highly conserved in human and fly. In addition, de novo motif searches generated a refined description of the U-rich downstream sequence (DSE) element, which shows more divergence between the two species. These refined motifs are applied, within a Hidden Markov Model (HMM) framework, to predict mRNA cleavage sites.

CONCLUSION

We demonstrate that the DSE is a specific motif in both human and Drosophila. These findings shed light on the sequence correlates of a highly conserved biological process, and improve in silico prediction of 3' mRNA cleavage and polyadenylation sites.

Abundant transcripts of malting barley identified by serial analysis of gene expression (SAGE)

Serial analysis of gene expression (SAGE) was applied to the major cereal crop barley (Hordeum vulgare) to characterize the transcriptional profile of grain during the malting process. Seven SAGE libraries were generated from seed at different time points during malting, in addition to one library from dry mature seed. A total of 155,206 LongSAGE tags, representing 41,909 unique sequences, was generated. This study reports an in-depth analysis of the most abundant transcripts from each of eight specific time points in a malting barley time course. The 100 most abundant tags from each library were analysed to identify the putative functional role of highly abundant transcripts. The largest functional groups included transcripts coding for stress response and cell defence, ribosomal proteins and storage proteins. The most abundant tag represented B22EL8, a barley metallothionein, which showed significant up-regulation across the malting time course. Considerable changes in the abundance profiles of some of the highly abundant tags occurred at 24 h post-steeping, indicating that it may be an important time point for gene expression changes associated with barley seed germination.

A systematic analysis of disease-associated variants in the 3' regulatory regions of human protein-coding genes I: general principles and overview

The 3' regulatory regions (3' RRs) of human genes play an important role in regulating mRNA 3' end formation, stability/degradation, nuclear export, subcellular localization and translation and are consequently rich in regulatory elements. Although 3' RRs contain only approximately 0.2% of known disease-associated mutations, this is likely to represent a rather conservative estimate of their actual prevalence. In an attempt to catalogue 3' RR-mediated disease and also to gain a greater understanding of the functional role of regulatory elements within 3' RRs, we have performed a systematic analysis of disease-associated 3' RR variants; 121 3' RR variants in 94 human genes were collated. These included 17 mutations in the upstream core polyadenylation signal sequence (UCPAS), 81 in the upstream sequence (USS) between the translational termination codon and the UCPAS, 6 in the left arm of the 'spacer' sequence (LAS) between the UCPAS and the pre-mRNA cleavage site (CS), 3 in the right arm of the 'spacer' sequence (RAS) or downstream core polyadenylation signal sequence (DCPAS) and 7 in the downstream sequence (DSS) of the 3'-flanking region, with 7 further mutations being treated as isolated examples. All the UCPAS mutations and the rather unusual cases of the DMPK, SCA8, FCMD and GLA mutations exert a significant effect on the mRNA phenotype and are usually associated with monogenic disease. By contrast, most of the remaining variants are polymorphisms that exert a comparatively minor influence on mRNA expression, but which may nevertheless predispose to or otherwise modify complex clinical phenotypes. Considerable efforts have been made to validate/elucidate the mechanisms through which the 3' untranslated region (3' UTR) variants affect gene expression. It is hoped that the integrative approach employed here in the study of naturally occurring variants of actual or potential pathological significance will serve to complement ongoing efforts to identify all functional regulatory elements in the human genome.

AltTrans: transcript pattern variants annotated for both alternative splicing and alternative polyadenylation

BACKGROUND

The three major mechanisms that regulate transcript formation involve the selection of alternative sites for transcription start (TS), splicing, and polyadenylation. Currently there are efforts that collect data & annotation individually for each of these variants. It is important to take an integrated view of these data sets and to derive a data set of alternate transcripts along with consolidated annotation. We have been developing in the past computational pipelines that generate value-added data at genome-scale on individual variant types; these include AltSplice on splicing and AltPAS on polyadenylation. We now extend these pipelines and integrate the resultant data sets to facilitate an integrated view of the contributions from splicing and polyadenylation in the formation of transcript variants.

DESCRIPTION

The AltSplice pipeline examines gene-transcript alignments and delineates alternative splice events and splice patterns; this pipeline is extended as AltTrans to delineate isoform transcript patterns for each of which both introns/exons and 'terminating' polyA site are delineated; EST/mRNA sequences that qualify the transcript pattern confirm both the underlying splicing and polyadenylation. The AltPAS pipeline examines gene-transcript alignments and delineates all potential polyA sites irrespective of underlying splicing patterns. Resultant polyA sites from both AltTrans and AltPAS are merged. The generated database reports data on alternative splicing, alternative polyadenylation and the resultant alternate transcript patterns; the basal data is annotated for various biological features. The data (named as integrated AltTrans data) generated for both the organisms of human and mouse is made available through the Alternate Transcript Diversity web site at http://www.ebi.ac.uk/atd/.

CONCLUSION

The reported data set presents alternate transcript patterns that are annotated for both alternative splicing and alternative polyadenylation. Results based on current transcriptome data indicate that the contribution of alternative splicing is larger than that of alternative polyadenylation.

Characterization of the arom gene in Rhizoctonia solani, and transcription patterns under stable and induced hypovirulence conditions

The quinate pathway is induced by quinate in the wild-type virulent Rhizoctonia solani isolate Rhs 1AP but is constitutive in the hypovirulent, M2 dsRNA-containing isolate Rhs 1A1. Constitutive expression of the quinate pathway results in downregulation of the shikimate pathway, which includes the pentafunctional arom gene in Rhs 1A1. The arom gene has 5,323 bp including five introns as opposed to a single intron found in arom in ascomycetes. A 199-bp upstream sequence has a GC box, no TATAA box, but two GTATTAGA repeats. The largest arom transcript is 5,108 nucleotides long, excluding the poly(A) tail. It contains an open reading frame of 4,857 bases, coding for a putative 1,618-residue pentafunctional AROM protein. A Kozak sequence (GCGCCATGG) is present between +127 and +135. The 5'-end of the arom mRNA includes two nucleotides (UA) that are not found in the genomic sequence, and are probably added post-transcriptionally. Size and sequence heterogeneity were observed at both 5'- and 3'-end of the mRNA. Northern blot and suppression subtractive hybridization analyses showed that presence of a low amount of quinate, inducer of the quinate pathway, resulted in increased levels of arom mRNA, consistent with the compensation effect observed in ascomycetes.

SuperSAGE combined with PCR walking allows global gene expression profiling of banana (Musa acuminata), a non-model organism

Super-serial analysis of gene expression (SuperSAGE) was used to characterize, for the first time, the global gene expression pattern in banana (Musa acuminata). A total of 10,196 tags were generated from leaf tissue, representing 5,292 expressed genes. Forty-nine tags of the top 100 most abundantly expressed transcripts were annotated by homology to cDNA or EST sequences. Typically for leaf tissue, analysis of the transcript profiles showed that the majority of the abundant transcripts are involved in energy production, mainly photosynthesis. However, the most abundant tag was derived from a type 3 metallothionein transcript, which accounted for nearly 3% of total transcripts analysed. Furthermore, the 26-bp long SuperSAGE tags were applied in 3'-rapid amplification of cDNA ends (3'RACE) for the identification of unknown tags. In combination with thermal asymmetric interlaced PCR (TAIL-PCR), this allowed the recovery of a full gene sequence of a novel NADPH:protochlorophyllide oxidoreductase, the key enzyme in chlorophyll biosynthesis. SuperSAGE in conjunction with 3'RACE and TAIL-PCR will be a powerful tool for transcriptomics of non-model, but otherwise important organisms.

Evaluation of the similarity of gene expression data estimated with SAGE and Affymetrix GeneChips

Background

Serial Analysis of Gene Expression (SAGE) and microarrays have found awidespread application, but much ambiguity exists regarding the evaluation of these technologies. Cross-platform utilization of gene expression data from the SAGE and microarray technology could reduce the need for duplicate experiments and facilitate a more extensive exchange of data within the research community. This requires a measure for the correspondence of the different gene expression platforms. To date, a number of cross-platform evaluations (including a few studies using SAGE and Affymetrix GeneChips) have been conducted showing a variable, but overall low, concordance. This study evaluates these overall measures and introduces the between-ratio difference as a concordance measure pergene.

Results

In this study, gene expression measurements of Unigene clusters represented by both Affymetrix GeneChips HG-U133A and SAGE were compared using two independent RNA samples. After matching of the data sets the final comparison contains a small data set of 1094 unique Unigene clusters, which is unbiased with respect to expression level. Different overall correlation approaches, like Up/Down classification, contingency tables and correlation coefficients were used to compare both platforms. In addition, we introduce a novel approach to compare two platforms based on the calculation of differences between expression ratios observed in each platform for each individual transcript. This approach results in a concordance measure per gene (with statistical probability value), as opposed to the commonly used overall concordance measures between platforms.

Conclusion

We can conclude that intra-platform correlations are generally good, but that overall agreement between the two platforms is modest. This might be due to the binomially distributed sampling variation in SAGE tag counts, SAGE annotation errors and the intensity variation between probe sets of a single gene in Affymetrix GeneChips. We cannot identify or advice which platform performs better since both have their (dis)-advantages. Therefore it is strongly recommended to perform follow-up studies of interesting genes using additional techniques. The newly introduced between-ratio difference is a filtering-independent measure for between-platform concordance. Moreover, the between-ratio difference per gene can be used to detect transcripts with similar regulation on both platforms.

Pre-cholesterol precursors in gametogenesis

Meiosis activating sterols (MAS) are biologically active post-lanosterol intermediates of cholesterol biosynthesis that are synthetized primarily in the gonads, including the sperm. MAS reinitiate the meiosis of oocytes in vitro while in vivo they seem to contribute to the oocyte quality and the progression of meiosis. The mRNAs for the MAS-producing enzyme lanosterol 14alpha-demethylase (CYP51) arise by alternative poly (A) signal selection. Only signals with low cleavage activity are used in the testis. Translation of mammalian CYP51s starts at one of the tandem in-frame ATGs. CYP51 protein of the bull is shorter compared to the human due to the usage of a more downstream translation start site. CYP51 proteins are post-translationally modified by glycosylations in the Golgi and on acrosomal membranes of the sperm. Green fluorescence protein-based ex vivo system has been developed to aid studying the intracellular transport of the MAS-producing CYP51. The influence of the post-translational modifications on MAS-synthesizing capacity is under investigation.

PolyA_DB: a database for mammalian mRNA polyadenylation

Messenger RNA polyadenylation is one of the key post-transcriptional events in eukaryotic cells. A large number of genes in mammalian species can undergo alternative polyadenylation, which leads to mRNAs with variable 3' ends. As the 3' end of mRNAs often contains cis elements important for mRNA stability, mRNA localization and translation, the implications of the regulation of polyadenylation can be multifold. Alternative polyadenylation is controlled by cis elements and trans factors, and is believed to occur in a tissue- or disease-specific manner. Given the availability of many databases devoted to other aspects of mRNA metabolism, such as transcriptional initiation and splicing, systematic information on polyadenylation, including alternative polyadenylation and its regulation, is noticeably lacking. Here, we present a database named polyA_DB, through which we strive to provide several types of information regarding polyadenylation in mammalian species: (i) polyadenylation sites and their locations with respect to the genomic structure of genes; (ii) cis elements surrounding polyadenylation sites; (iii) comparison of polyadenylation configuration between orthologous genes; and (iv) tissue/organ information for alternative polyadenylation sites. Currently, polyA_DB contains 45,565 polyadenylation sites for 25,097 human and mouse genes, representing the most comprehensive polyadenylation database till date. The database is accessible via the website (http://polya.umdnj.edu/polyadb).

A comparative analysis of transcript abundance using SAGE and Affymetrix arrays

A number of methods are currently used for gene expression profiling. They differ in scale, economy and sensitivity. We present the results of a direct comparison between serial analysis of gene expression (SAGE) and the Barley1 Affymetrix GeneChip. Both technology platforms were used to obtain quantitative measurements of transcript abundance using identical RNA samples and assessed for their ability to quantify differential gene expression. For SAGE, a total of 82,122 tags were generated from two independent libraries representing whole developing barley caryopsis and dissected embryos. The Barley1 GeneChip contains 22,791 probe sets. Results obtained from both methods are generally comparable, indicating that both will lead to similar conclusions regarding transcript levels and differential gene expression. However, excluding singletons, 24.4% of the unique SAGE tags had no corresponding probe set on the Barley1 array indicating that a broader snapshot of gene expression was obtained by SAGE. Discrepancies were observed for a number of "genes" and these are discussed.

A large-scale analysis of mRNA polyadenylation of human and mouse genes

mRNA polyadenylation is a critical cellular process in eukaryotes. It involves 3′ end cleavage of nascent mRNAs and addition of the poly(A) tail, which plays important roles in many aspects of the cellular metabolism of mRNA. The process is controlled by various cis-acting elements surrounding the cleavage site, and their binding factors. In this study, we surveyed genome regions containing cleavage sites [herein called poly(A) sites], for 13 942 human and 11 155 mouse genes. We found that a great proportion of human and mouse genes have alternative polyadenylation (∼54 and 32%, respectively). The conservation of alternative polyadenylation type or polyadenylation configuration between human and mouse orthologs is statistically significant, indicating that alternative polyadenylation is widely employed by these two species to produce alternative gene transcripts. Genes belonging to several functional groups, indicated by their Gene Ontology annotations, are biased with respect to polyadenylation configuration. Many poly(A) sites harbor multiple cleavage sites (51.25% human and 46.97% mouse sites), leading to heterogeneous 3′ end formation for transcripts. This implies that the cleavage process of polyadenylation is largely imprecise. Different types of poly(A) sites, with regard to their relative locations in a gene, are found to have distinct nucleotide composition in surrounding genomic regions. This large-scale study provides important insights into the mechanism of polyadenylation in mammalian species and represents a genomic view of the regulation of gene expression by alternative polyadenylation.

Mining SAGE data allows large-scale, sensitive screening of antisense transcript expression

As a growing number of complementary transcripts, susceptible to exert various regulatory functions, are being found in eukaryotes, high throughput analytical methods are needed to investigate their expression in multiple biological samples. Serial Analysis of Gene Expression (SAGE), based on the enumeration of directionally reliable short cDNA sequences (tags), is capable of revealing antisense transcripts. We initially detected them by observing tags that mapped on to the reverse complement of known mRNAs. The presence of such tags in individual SAGE libraries suggested that SAGE datasets contain latent information on antisense transcripts. We raised a collection of virtual tags for mining these data. Tag pairs were assembled by searching for complementarities between 24-nt long sequences centered on the potential SAGE-anchoring sites of well-annotated human expressed sequences. An analysis of their presence in a large collection of published SAGE libraries revealed transcripts expressed at high levels from both strands of two adjacent, oppositely oriented, transcription units. In other cases, the respective transcripts of such cis-oriented genes displayed a mutually exclusive expression pattern or were co-expressed in a small number of libraries. Other tag pairs revealed overlapping transcripts of trans-encoded unique genes. Finally, we isolated a group of tags shared by multiple transcripts. Most of them mapped on to retroelements, essentially represented in humans by Alu sequences inserted in opposite orientations in the 3'UTR of otherwise different mRNAs. Registering these tags in separate files makes possible computational searches focused on unique sense-antisense pairs. The method developed in the present work shows that SAGE datasets constitute a major resource of rapidly investigating with high sensitivity the expression of antisense transcripts, so that a single tag may be detected in one library when screening a large number of biological samples.

HLA DQA1 genes generate multiple transcripts by alternative splicing and polyadenylation of the 3' untranslated region

Regulation of the human leucocyte antigen (HLA) class II genes expression is an important field in immunology, because these molecules play a crucial role in the function of the immune system. HLA DQ genes expression is a complex phenomenon regulated at both transcriptional and post-transcriptional levels. In this study, we have investigated the post-transcriptional mechanisms accounting for alleles-dependent length polymorphism of DQA1 mRNA. We have first sequenced the genomic DNA encoding the 3' untranslated region (UTR) of DQA1 *0101, *0102, *0103, *0201, *0301, *0401, and *0501 alleles. We have identified two competing splicing sites: a unique splicing donor site AG/GTA located 20 nucleotides downstream from the stop codon associated to two spliced acceptor sequences, approximately 165 and approximately 370 nucleotides downstream. In addition, three polyadenylation signals have been identified, respectively, at approximately 475, approximately 795, and approximately 855 nucleotides downstream from the stop codon. Subsequently, we have analyzed mRNAs derived from DQA1 alleles in homozygous B lymphoblastoid cell lines by reverse transcriptase-polymerase chain reaction. We show that allele-dependent length polymorphism of DQA1 mRNA-3' UTR results from a combination of differential splicing and alternative polyadenylations. Four mRNA isoforms (two spliced variant cleaved at two distinct polyadenylation sites) were detected in DQA1 *0101, *0102, and *0103 homozygous cell lines, and six mRNA species (three spliced variant cleaved at two polyadenylation-sequence signal) were generated by the other four alleles. Possible advantages for cells to generate multiple transcripts previously undetected are discussed.

Complex patterns of transcription at the insertion site of a retrotransposon in the mouse

Here we report that transcriptional effects of the insertion of a retrotransposon can occur simultaneously both upstream and downstream of the insertion site. We have identified an intra-cisternal A particle (IAP) retrotransposon in intron 6 of a gene that we have named Cabp (CDK5 activator binding protein). The presence of the IAP is associated with an aberrant transcript initiating from a cryptic promoter in the IAP, reading out into the adjacent Cabp gene sequence. The expression of this transcript is highly variable among isogenic mice within the C57BL/6J strain and so Cabp(IAP) can be classified as a metastable epiallele. As expected, the presence or absence of the transcript correlates with differential DNA methylation of the 5' LTR of the IAP. More surprisingly, in mice where the retrotransposon is unmethylated and presumably transcriptionally active, we find a number of short Cabp transcripts which initiate at the normal 5' end of the gene but terminate prematurely, just 5' of the retrotransposon. This is the first report of a retrotransposon having both upstream and downstream effects on transcription at the site of insertion and it suggests that alternative polyadenylation may sometimes be caused by a downstream convergent transcription unit.

Modified serial analysis of gene expression method for construction of gene expression profiles of microbial eukaryotic species

Serial analysis of gene expression (SAGE) is a powerful approach for the identification of differentially expressed genes, providing comprehensive and quantitative gene expression profiles in the form of short tag sequences. Each tag represents a unique transcript, and the relative frequencies of tags in the SAGE library are equal to the relative proportions of the transcripts they represent. One of the major obstacles in the preparation of SAGE libraries from microorganisms is the requirement for large amounts of starting material (i.e., mRNA). Here, we present a novel approach for the construction of SAGE libraries from small quantities of total RNA by using Y linkers to selectively amplify 3' cDNA fragments. To validate this method, we constructed comprehensive gene expression profiles of the toxic dinoflagellate Pfiesteria shumwayae. SAGE libraries were constructed from an actively toxic fish-fed culture of P. shumwayae and from a recently toxic alga-fed culture. P. shumwayae-specific gene transcripts were identified by comparison of tag sequences in the two libraries. Representative tags with frequencies ranging from 0.026 to 3.3% of the total number of tags in the libraries were chosen for further analysis. Expression of each transcript was confirmed in separate control cultures of toxic P. shumwayae. The modified SAGE method described here produces gene expression profiles that appear to be both comprehensive and quantitative, and it is directly applicable to the study of gene expression in other environmentally relevant microbial species.

5′-end SAGE for the analysis of transcriptional start sites

Identification of the mRNA start site is essential in establishing the full-length cDNA sequence of a gene and analyzing its promoter region, which regulates gene expression. Here we describe the development of a 5'-end serial analysis of gene expression (5' SAGE) that can be used to globally identify transcriptional start sites and the frequency of individual mRNAs. Of the 25,684 5' SAGE tags in the HEK293 human cell library, 19,893 matched to the human genome. Among 15,448 tags in one locus of the genome, 85.8%-96.1% of the 5' SAGE tags were assigned within -500 to +200 nt of mRNA start sites using the RefSeq, UniGene and DBTSS databases. This technique should facilitate 5'-end transcriptome analysis in a variety of cells and tissues.

Rna14-Rna15 assembly mediates the RNA-binding capability of Saccharomyces cerevisiae cleavage factor IA

The Rna14-Rna15 complex is a core component of the cleavage factor IA RNA-processing complex from Saccharomyces cerevisiae. To understand the assembly and RNA-binding properties, we have isolated and characterized the Rna14-Rna15 complex using a combination of biochemical and biophysical methods. Analysis of the purified complex, using transmission electron microscopy, reveals that the two proteins assemble into a kinked rod-shaped structure and that these rods are able to further self-associate. Analytical ultracentrifugation reveals that Rna14 mediates this association and facilitates assembly of an A2B2 tetramer (M(r) 230 000), where relatively compact Rna14-Rna15 heterodimers are in rapid equilibrium with tetramers that have a more extended shape. The Rna14-Rna15 complex, unlike the individual components, binds to an RNA oligonucleotide derived from the 3'-untranslated region of the S.cerevisiae GAL7 gene. Based on these structural and thermodynamic data, we propose that CFIA assembly regulates RNA-binding activity.

A preoperative diagnostic test that distinguishes benign from malignant thyroid carcinoma based on gene expression

Accurate diagnosis of thyroid tumors is challenging. A particular problem is distinguishing between follicular thyroid carcinoma (FTC) and benign follicular thyroid adenoma (FTA), where histology of fine-needle aspirates is not conclusive. It is often necessary to remove healthy thyroid to rule out carcinoma. In order to find markers to improve diagnosis, we quantified gene transcript expression from FTC, FTA, and normal thyroid, revealing 73 differentially expressed transcripts (P < or = 0.0001). Using an independent set of 23 FTCs, FTAs, and matched normal thyroids, 17 genes with large expression differences were tested by real-time RT-PCR. Four genes (DDIT3, ARG2, ITM1, and C1orf24) differed between the two classes FTC and FTA, and a linear combination of expression levels distinguished FTC from FTA with an estimated predictive accuracy of 0.83. Furthermore, immunohistochemistry for DDIT3 and ARG2 showed consistent staining for carcinoma in an independent set 59 follicular tumors (estimated concordance, 0.76; 95% confidence interval, [0.59, 0.93]). A simple test based on a combination of these markers might improve preoperative diagnosis of thyroid nodules, allowing better treatment decisions and reducing long-term health costs.

Genes differentially expressed in thyroid carcinoma identified by comparison of SAGE expression profiles

To identify transcripts that distinguish malignant from benign thyroid disease serial analysis of gene expression (SAGE) profiles of papillary thyroid carcinoma and of normal thyroid are compared. Of the 21,000 tags analyzed, 204 tags are differentially expressed with statistical significance in the tumor. Thyroid tumor specificity of these transcripts is determined in silico using the tissue preferential expression (TPE) algorithm. TPE values demonstrate that 42 tags of the 204 are thyroid tumor specific. BC013035, a cDNA encoding a novel protein, is up-regulated from 0 to 24 tags in the thyroid tumor SAGE library. In a tissue panel of 30 thyroid tumors and 12 controls, it has an expression pattern similar to thyroid peroxidase, indicating possible involvement of BC013035 in thyroid differentiation. A tag coding for extracellular matrix protein 1 (ECM1) is absent in the normal thyroid SAGE library and present 55 times in the tumor. ECM1, a protein recently associated with angiogenesis and expressed in metastatic breast carcinoma, is up-regulated in 50% of all thyroid carcinoma and absent in normal controls and follicular adenoma. In conclusion, SAGE analysis and subsequent determination of TPE values facilitates the rapid distinction of genes specifically expressed in cancer tissues.

Serial Analysis of Gene Expression: Applications in Human Studies

Serial analysis of gene expression (SAGE) is a powerful tool, which provides quantitative and comprehensive expression profile of genes in a given cell population. It works by isolating short fragments of genetic information from the expressed genes that are present in the cell being studied. These short sequences, called SAGE tags, are linked together for efficient sequencing. The frequency of each SAGE tag in the cloned multimers directly reflects the transcript abundance. Therefore, SAGE results in an accurate picture of gene expression at both the qualitative and the quantitative levels. It does not require a hybridization probe for each transcript and allows new genes to be discovered. This technique has been applied widely in human studies and various SAGE tags/SAGE libraries have been generated from different cells/tissues such as dendritic cells, lung fibroblast cells, oocytes, thyroid tissue, B-cell lymphoma, cultured keratinocytes, muscles, brain tissues, sciatic nerve, cultured Schwann cells, cord blood-derived mast cells, retina, macula, retinal pigment epithelial cells, skin cells, and so forth. In this review we present the updated information on the applications of SAGE technology mainly to human studies.

The intronic prothrombin 19911A>G polymorphism influences splicing efficiency and modulates effects of the 20210G>A polymorphism on mRNA amount and expression in a stable reporter gene assay system

The common prothrombin gene cleavage site mutation 20210G>A is associated with elevated prothrombin levels and thrombosis. The pathomechanism of the 20210G>A mutation was explained by increased mRNA formation and/or more efficient translation. Human studies also showed an influence of the intronic 19911A>G polymorphism on prothrombin activity. We established HepG2 cell lines stably transfected with prothrombin mini-genes containing the last 2 prothrombin exons, the last intron, 3' untranslated region (UTR), and flanking sequence. The highest mRNA expression and protein activity resulted from the mutant haplotype 19911A-20210A. Haplotypes with wild-type cleavage site (19911A-20210G, 19911G-20210G) also differed significantly as a consequence of the intronic 19911 mutation; the 19911G-20210G haplotype showed lower expression than the 19911A-20210G haplotype, whereas previous clinical studies have reported elevated prothrombin activity with the 19911G-20210G haplotype. The cleavage site pattern was homogeneous with 20210A, which may cause a favorable intracellular processing, and heterogeneous with 20210G. In an independent assay for splicing efficiency, 19911G showed about 30% higher efficiency than 19911A. We conclude that the intronic 19911A>G single nucleotide polymorphism is itself functional and changes splicing efficiency by altering a known functional pentamer motif. Further studies are needed to define the value of additional prothrombin 19911 genotyping for thrombophilia screening, especially in cases heterozygous for 20210G>A.

A sensitive transcriptome analysis method that can detect unknown transcripts

We have developed an AFLP-based gene expression profiling method called 'high coverage expression profiling' (HiCEP) analysis. By making improvements to the selective PCR technique we have reduced the rate of false positive peaks to approximately 4% and consequently the number of peaks, including overlapping peaks, has been markedly decreased. As a result we can determine the relationship between peaks and original transcripts unequivocally. This will make it practical to prepare a database of all peaks, allowing gene assignment without having to isolate individual peaks. This precise selection also enables us to easily clone peaks of interest and predict the corresponding gene for each peak in some species. The procedure is highly reproducible and sensitive enough to detect even a 1.2-fold difference in gene expression. Most importantly, the low false positive rate enables us to analyze gene expression with wide coverage by means of four instead of six nucleotide recognition site restriction enzymes for fingerprinting mRNAs. Therefore, the method detects 70-80% of all transcripts, including non-coding transcripts, unknown and known genes. Moreover, the method requires no sequence information and so is applicable even to eukaryotes for which there is no genome information available.

Assessment of SAGE in transcript identification

An essential step in Serial Analysis of Gene Expression (SAGE) is tag mapping, which refers to the unambiguous determination of the gene represented by a SAGE tag. Current resources for tag mapping are incomplete, and thus do not allow assessment of the efficacy of SAGE in transcript identification. A method of tag mapping is described here and applied to the Drosophila melanogaster and Caenorhabditis elegans genomes, which permits detailed SAGE assessment and provides tag-mapping resources that were unavailable previously for these organisms. In our method, a conceptual transcriptome is constructed using genomic sequence and annotation by extending predicted coding regions to include UTRs on the basis of EST and cDNA alignments, UTR length distributions, and polyadenylation signals. Analysis of extracted tags suggests that, using the standard SAGE procedure, expression of 8% of D. melanogaster and 15% of C. elegans genes cannot be detected unambiguously by SAGE due to shared sequence or lack of NlaIII-anchoring enzyme sites. Both increasing tag length by 2-3 bp and using Sau3A instead of NlaIII as the anchoring enzyme increases potential for transcript detection. This work identifies and quantifies genes not amenable to SAGE analysis, in addition to providing tag-to-gene mappings for two model organisms.

Transcript identification by analysis of short sequence tags--influence of tag length, restriction site and transcript database

There exist a number of gene expression profiling techniques that utilize restriction enzymes for generation of short expressed sequence tags. We have studied how the choice of restriction enzyme influences various characteristics of tags generated in an experiment. We have also investigated various aspects of in silico transcript identification that these profiling methods rely on. First, analysis of 14 248 mRNA sequences derived from the RefSeq transcript database showed that 1-30% of the sequences lack a given restriction enzyme recognition site. Moreover, 1-5% of the transcripts have recognition sites located less than 10 bases from the poly(A) tail. The uniqueness of 10 bp tags lies in the range 90-95%, which increases only slightly with longer tags, due to the existence of closely related transcripts. Furthermore, 3-30% of upstream 10 bp tags are identical to 3' tags, introducing a risk of misclassification if upstream tags are present in a sample. Second, we found that a sequence length of 16-17 bp, including the recognition site, is sufficient for unique transcript identification by BLAST based sequence alignment to the UniGene Human non-redundant database. Third, we constructed a tag-to-gene mapping for UniGene and compared it to an existing mapping database. The mappings agreed to 79-83%, where the selection of representative sequences in the UniGene clusters is the main cause of the disagreement. The results of this study may serve to improve the interpretation of sequence-based expression studies and the design of hybridization arrays, by identifying short tags that have a high reliability and separating them from tags that carry an inherent ambiguity in their capacity to discriminate between genes. To this end, supplementary information in the form of a web companion to this paper is located at http:// biobase.biotech.kth.se/tagseq.

Comprehensive sampling of gene expression in human cell lines with massively parallel signature sequencing

Whereas information is rapidly accumulating about the structure and position of genes encoded in the human genome, less is known about the complexity and relative abundance of their expression in individual human cells and tissues. Here, we describe the characteristics of the transcriptomes of two cultured cell lines, HB4a (normal breast epithelium) and HCT-116 (colon adenocarcinoma), using massively parallel signature sequencing (MPSS). We generated in excess of 10(7) short signature sequences per cell line, thus providing a comprehensive snapshot of gene expression, within the technical limitations of the method. The number of genes expressed at one copy per cell or more in either of the lines was estimated to be between 10,000 and 15,000. The vast majority of the transcripts found in these cells can be mapped to known genes and their polyadenylation variants. Among the genes that could be identified from their signature sequences, approximately 8,500 were expressed by both cell lines, whereas 6,000 showed cellular specificity. Taking into account sequence tags that map uniquely to the genome but not to known transcripts, overall the data are consistent with an upper limit of 17,000 for the total number of genes expressed at more than one copy per cell in one or both of the two cell lines examined.

Identifying novel transcripts and novel genes in the human genome by using novel SAGE tags

The number of genes in the human genome is still a controversial issue. Whereas most of the genes in the human genome are said to have been physically or computationally identified, many short cDNA sequences identified as tags by use of serial analysis of gene expression (SAGE) do not match these genes. By performing experimental verification of more than 1,000 SAGE tags and analyzing 4,285,923 SAGE tags of human origin in the current SAGE database, we examined the nature of the unmatched SAGE tags. Our study shows that most of the unmatched SAGE tags are truly novel SAGE tags that originated from novel transcripts not yet identified in the human genome, including alternatively spliced transcripts from known genes and potential novel genes. Our study indicates that by using novel SAGE tags as probes, we should be able to identify efficiently many novel transcripts/novel genes in the human genome that are difficult to identify by conventional methods.

Long-range heterogeneity at the 3' ends of human mRNAs

The publication of a draft of the human genome and of large collections of transcribed sequences has made it possible to study the complex relationship between the transcriptome and the genome. In the work presented here, we have focused on mapping mRNA 3' ends onto the genome by use of the raw data generated by the expressed sequence tag (EST) sequencing projects. We find that at least half of the human genes encode multiple transcripts whose polyadenylation is driven by multiple signals. The corresponding transcript 3' ends are spread over distances in the kilobase range. This finding has profound implications for our understanding of gene expression regulation and of the diversity of human transcripts, for the design of cDNA microarray probes, and for the interpretation of gene expression profiling experiments.

A novel root-specific gene, MIC-3, with increased expression in nematode-resistant cotton (Gossypium hirsutum L.) after root-knot nematode infection

A full-length cDNA, MIC-3, has been identified from a lambda ZAPII cDNA library constructed from the mRNA of nematode-resistant cotton (Gossypium hirsutum L.) roots after infection with root-knot nematode (Meloidogyne incognita). The putative open reading frame of MIC-3 encoded a protein of 141 amino acids with a calculated molecular mass of 15.3 kDa. Seven alternative polyadenylation sites have been identified for the MIC-3 transcripts, and the major transcripts are the longest ones. The MIC-3 gene contains a single intron within its coding region and belongs to a novel, multi-gene family containing up to six members. Expression of MIC-3 is root localized and specifically enhanced in the nematode induced, immature galls of resistant cotton line M-249, suggesting that MIC-3 may play a critical role in the resistance response to root-knot nematode.

A partial transcriptome of human epidermis

Serial analysis of gene expression (SAGE) is a powerful technique for global expression profiling without prior knowledge of the genes of interest. We carried out SAGE analysis of purified keratinocytes derived from human skin biopsy specimens, resulting in a partial transcriptome of human epidermis. We identified 7645 unique SAGE tags with quantitative information from 15,131 collected SAGE tags obtained from approximately 3 x 10(6) epidermal cells. This catalog contains a large number of genes that were not previously known to be expressed by human epidermis. Comparison with the databases of all known human SAGE tags allowed us to identify a number of keratinocyte-specific tags that putatively correspond to formerly unknown genes. Surprisingly, human epidermal keratinocytes in vivo show relatively low expression levels of genes typically associated with epidermal differentiation, whereas the expression levels of housekeeping genes are considerably higher than in cultured keratinocytes. This study provides a first step toward a transcriptome of human epidermis and, as such, harbors a wealth of information to identify genes involved in skin function, and candidate genes for genetic skin disorders.

Reexamining the polyadenylation signal: were we wrong about AAUAAA?

Polyadenylation is the process by which most eukaryotic mRNAs form their 3' ends. It was long held that polyadenylation required the sequence AAUAAA and that 90% of mRNAs had AAUAAA within 30 nucleotides of the site of poly(A) addition. More recent studies, aided by computer analysis of sequences made available in GenBank and expressed sequence tag (EST) databases, have suggested that the actual incidence of AAUAAA is much lower, perhaps as low as 50-60%. Reproductive biologists have long recognized that a large number of mRNAs in male germ cells of mammals lack AAUAAA but are otherwise normally polyadenylated. Recent research in our laboratory has uncovered a new form of an essential polyadenylation protein, tauCstF-64, that is most highly expressed in male germ cells, and to a smaller extent in the brain, and which we propose plays a significant role in AAUAAA-independent mRNA polyadenylation in germ cells.

Identification of alternate polyadenylation sites and analysis of their tissue distribution using EST data

Alternate polyadenylation affects a large fraction of higher eucaryote mRNAs, producing mature transcripts with 3' ends of variable length. This variation is poorly represented in the current transcript catalogs derived from whole genome sequences, mostly because such posttranscriptional events are not detectable directly at the DNA level. Alternate polyadenylation of an mRNA is better understood by comparison to EST databases. Comparing ESTs to mRNAs, however, is a difficult task subjected to the pitfalls of internal priming, presence of intron sequences, repeated elements, chimerical ESTs or matches with EST from paralogous genes. We present here a computer program that addresses these problems and displays ESTs matches to a query mRNA sequence to predict alternate polyadenylation and to suggest library-specific forms. The output highlights effective polyadenylation signals, possible sources of artifacts such as A-rich stretches in the mRNA sequences, and allows for a direct visualization of EST libraries using color codes. Statistical biases in the distribution of alternative mRNA forms among EST libraries were systematically sought. About 1450 human and 200 mouse mRNAs displayed such biases, suggesting in each case a tissue- or disease-specific regulation of polyadenylation.