Targeting a complex transcriptome: The construction of the mouse full-length cDNA encyclopedia

We report the construction of the mouse full-length cDNA encyclopedia, the most extensive view of a complex transcriptome, on the basis of preparing and sequencing 246 libraries. Before cloning, cDNAs were enriched in full-length by Cap-Trapper, and in most cases, aggressively subtracted/normalized. We have produced 1,442,236 successful 3'-end sequences clustered into 171,144 groups, from which 60,770 clones were fully sequenced cDNAs annotated in the FANTOM-2 annotation. We have also produced 547,149 5' end reads, which clustered into 124,258 groups. Altogether, these cDNAs were further grouped in 70,000 transcriptional units (TU), which represent the best coverage of a transcriptome so far. By monitoring the extent of normalization/subtraction, we define the tentative equivalent coverage (TEC), which was estimated to be equivalent to >12,000,000 ESTs derived from standard libraries. High coverage explains discrepancies between the very large. numbers of clusters (and TUs) of this project, which also include non-protein-coding RNAs, and the lower gene number estimation of genome annotations. Altogether, S'-end clusters identify regions that are potential promoters for 8637 known genes and S'-end clusters suggest the presence of almost 63,000 transcriptional starting points. An estimate of the frequency of polyadenylation signals suggests that at least half of the singletons in the EST set represent real mRNAs. Clones accounting for about half of the predicted TUs await further sequencing. The continued high-discovery rate suggests that the task of transcriptome discovery is not yet complete.

Selective loss of NMDA receptor NR1 subunit isoforms in Alzheimer's disease

Previous work had shown that the ratio of NMDA receptor NR1 subunit mRNA transcripts containing an N-terminal splice cassette to those that do not is markedly lower in regions of the Alzheimer's disease (AD) brain that are susceptible to pathological damage, compared with spared regions in the same cases or homotropic regions in controls. To elucidate the origins of this difference in proportionate expression, we measured the absolute levels of each of the eight NR1 transcripts by quantitative internally standardized RT-PCR assay. Expression of transcripts with the cassette was strongly attenuated in susceptible regions of Alzheimer's brain, whereas expression of non-cassette transcripts differed little from that in controls. The expression of other NR1 splice variants was not associated with pathology relevant to disease status, although some combinations of splice cassettes were well maintained in AD cases. The population profile of NR1 transcripts in occipital cortex differed from the profiles in other brain regions studied. Western analysis confirmed that the expression of protein isoforms containing the N-terminal peptide was very low in susceptible areas of the Alzheimer's brain. Cells that express NR1 subunits with the N-terminal cassette may be selectively vulnerable to toxicity in AD.

Differential expression of N-methyl-D-aspartate receptor NR2 isoforms in Alzheimer's disease

We have previously shown that the expression of NMDA receptor NR1 subunit mRNA splice variants in Alzheimer's disease (AD) brain varies according to regional susceptibility to pathological damage. Here we investigated the expression of the modulatory NR2 subunits of the NMDA receptor using quantitative RT-PCR to assay all NR2 isoforms. Significantly lower expression of NR2A and NR2B transcripts was found in susceptible regions of AD brain, whereas expression of NR2C and NR2D transcripts did not differ from that in controls. Western blot analysis confirmed a lower expression of the NR2A and NR2B isoforms at the protein level. The results suggest that NR2 subunit composition may modulate NMDA receptor-mediated excitotoxicity. NMDA receptor dysfunction might give rise to the regionally selective pattern of neuronal loss that is characteristic of AD.

Whole genome analysis reveals a high incidence of non-optimal codons in secretory signal sequences of Escherichia coli

Translational pausing may occur due to a number of mechanisms, including the presence of non-optimal codons, and it is thought to play a role in the folding of specific polypeptide domains during translation and in the facilitation of signal peptide recognition during see-dependent protein targeting. In this whole genome analysis of Escherichia coli we have found that non-optimal codons in the signal peptide-encoding sequences of secretory genes are overrepresented relative to the mature portions of these genes; this is in addition to their overrepresentation in the 5'-regions of genes encoding non-secretory proteins. We also find increased non-optimal codon usage at the 3' ends of most E. coli genes, in both non-secretory and secretory sequences. Whereas presumptive translational pausing at the 5' and 3' ends of E. coli messenger RNAs may clearly have a general role in translation, we suggest that it also has a specific role in sec-dependent protein export, possibly in facilitating signal peptide recognition. This finding may have important implications for our understanding of how the majority of non-cytoplasmic proteins are targeted, a process that is essential to all biological cells. (C) 2004 Elsevier Inc. All rights reserved.

Expression of the cell surface mucin gene family in adenocarcinomas

Cell surface mucins are complex glycoproteins expressed on the apical membrane surface of mucosal epithelial cells. In malignant epithelial cells they are thought to influence cell adhesion, and are clinical targets for tumor immunotherapy and serum tumor marker assays. We have compared expression of MUC1, MUC3, MUC4, MUC11, MUC12 and MUC13 mRNA in epithelial cancers and/or cell lines with non-malignant tissues. In non-malignant tissues, MUC3, 4, 11, 12 and 13 were expressed at highest levels in gastrointestinal tissues, whereas MUC1 was more widely distributed. Significant down-regulation of the MUC4, MUC12 and MUC13 genes was observed in colonic cancers compared with normal tissue, whereas MUC1 was upregulated. In rectal cancers, levels of all six mucin genes were not significantly different to those in normal rectal tissues. Both MUC1 and MUC4 were down-regulated in gastric cancers, whereas cancer and normal tissue levels were similar for MUC3, 11, 12 and 13. In esophageal cancers there was a general trend toward higher levels than in normal tissue for MUC1, 3, 12 and 13. In ovarian cancers MUC1 levels were very high, whereas only low levels of all other mucins were observed. We also report expression in renal cell carcinomas, bladder carcinomas and breast cancer cell lines. The reported expression profiles of the cell surface mucin gene family will help direct biological and clinical studies of these molecules in mucosal biology, and in malignant and inflammatory diseases of epithelial tissues.

5 '-Untranslated regions with multiple upstream AUG codons can support low-level translation via leaky scanning and reinitiation

Upstream AUGs (uAUGs) and upstream open reading frames (uORFs) are common features of mRNAs that encode regulatory proteins and have been shown to profoundly influence translation of the main ORF. In this study, we employed a series of artificial 5'-untranslated regions (5'-UTRs) containing one or more uAUGs/uORFs to systematically assess translation initiation at the main AUG by leaky scanning and reinitiation mechanisms. Constructs containing either one or two uAUGs in varying contexts but without an in-frame stop codon upstream of the main AUG were used to analyse the leaky scanning mechanism. This analysis largely confirmed the ranking of different AUG contextual sequences that was determined previously by Kozak. In addition, this ranking was the same for both the first and second uAUGs, although the magnitude of initiation efficiency differed. Moreover, similar to10% of ribosomes exhibited leaky scanning at uAUGs in the most favourable context and initiated at a downstream AUG. A second group of constructs containing different numbers of uORFs, each with optimal uAUGs, were used to measure the capacity for reinitiation. We found significant levels of initiation at the main ORF even in constructs containing four uORFs, with nearly 10% of ribosomes capable of reinitiating five times. This study shows that for mRNAs containing multiple uORFs/uAUGs, ribosome reinitiation and leaky scanning are efficient mechanisms for initiation at their main AUGs.

Gene codon composition determines differentiation-dependent expression of a viral capsid gene in keratinocytes in vitro and in vivo

By establishing mouse primary keratinocytes (KCs) in culture, we were able, for the first time, to express papillomavirus major capsid (L1) proteins by transient transfection of authentic or codon-modified L1 gene expression plasmids. We demonstrate in vitro and in vivo that gene codon composition is in part responsible for differentiation-dependent expression of L1 protein in KCs. L1 mRNA was present in similar amounts in differentiated and undifferentiated KCs transfected with authentic or codon-modified L1 genes and had a similar half-life, demonstrating that L1 protein production is posttranscriptionally regulated. We demonstrate further that KCs substantially change their tRNA profiles upon differentiation. Aminoacyl-tRNAs from differentiated KCs but not undifferentiated KCs enhanced the translation of authentic L1 mRNA, suggesting that differentiation-associated change to tRNA profiles enhances L1 expression in differentiated KCs. Thus, our data reveal a novel mechanism for regulation of gene expression utilized by a virus to direct viral capsid protein expression to the site of virion assembly in mature KCs. Analysis of two structural proteins of KCs, involucrin and keratin 14, suggests that translation of their mRNAs is also regulated, in association with KC differentiation in vitro, by a similar mechanism

Antisense transcription in the mammalian transcriptome

Antisense transcription (transcription from the opposite strand to a protein-coding or sense strand) has been ascribed roles in gene regulation involving degradation of the corresponding sense transcripts (RNA interference), as well as gene silencing at the chromatin level. Global transcriptome analysis provides evidence that a large proportion of the genome can produce transcripts from both strands, and that antisense transcripts commonly link neighboring genes in complex loci into chains of linked transcriptional units. Expression profiling reveals frequent concordant regulation of sense/antisense pairs. We present experimental evidence that perturbation of an antisense RNA can alter the expression of sense messenger RNAs, suggesting that antisense transcription contributes to control of transcriptional outputs in mammals.

Differential expression of the GABA transporters GAT-1 and GAT-3 in brains of rats, cats, monkeys and humans

The homeostasis of GABA is critical to normal brain function. Extracellular levels of GABA are regulated mainly by plasmalemmal gamma-aminobutyric acid (GABA) transporters. Whereas the expression of GABA transporters has been extensively studied in rodents, validation of this data in other species, including humans, has been limited. As this information is crucial for our understanding of therapeutic options in human diseases such as epilepsy, we have compared, by immunocytochemistry, the distributions of the GABA transporters GAT-1 and GAT-3 in rats, cats, monkeys and humans. We demonstrate subtle differences between the results reported in the literature and our results, such as the predominance of GAT-1 labelling in neurons rather than astrocytes in the rat cortex. We note that the optimal localisation of GAT-1 in cats, monkeys and humans requires the use of an antibody against the human sequence carboxyl terminal region of GAT-1 rather than against the slightly different rat sequence. We demonstrate that GAT-3 is localised mainly to astrocytes in hindbrain and midbrain regions of rat brains. However, in species such as cats, monkeys and humans, additional strong immunolabelling of oligodendrocytes has also been observed. We suggest that differences in GAT distribution, especially the expression of GAT-3 by oligodendrocytes in humans, must be accommodated in extrapolating rodent models of GABA homeostasis to humans.

Block of voltage-gated potassium channels by Pacific ciguatoxin-1 contributes to increased neuronal excitability in rat sensory neurons

The present study investigated the actions of the polyether marine toxin Pacific ciguatoxin-1 (P-CTX-1) on neuronal excitability in rat dorsal root ganglion (DRG) neurons using patch-clamp recording techniques. Under current-clamp conditions, bath application of 2-20 nM P-CTX-1 caused a rapid, concentration-dependent depolarization of the resting membrane potential in neurons expressing tetrodotoxin (TTX)-sensitive voltage-gated sodium (Na-v,.) channels. This action was completely suppressed by the addition of 200 nM TTX to the external solution, indicating that this effect was mediated through TTX-sensitive Na-v channels. In addition, P-CTX-1 also prolonged action potential and afterhyperpolarization (AHP) duration. In a subpopulation of neurons, P-CTX-1 also produced tonic action potential firing, an effect that was not accompanied by significant oscillation of the resting membrane potential. Conversely, in neurons expressing TTX-resistant Na-v currents, P-CTX-1 failed to alter any parameter of neuronal excitability examined in this study. Under voltage-clamp conditions in rat DRG neurons, P-CTX-1 inhibited both delayed-rectifier and 'A-type' potassium currents in a dose-dependent manner, actions that Occurred in the absence of alterations to the voltage dependence of activation. These actions appear to underlie the prolongation of the action potential and AHP. and contribute to repetitive firing. These data indicate that a block of potassium channels contributes to the increase in neuronal excitability, associated with a modulation of Na-v. channel gating, observed clinically in response to ciguatera poisoning. (c) 2004 Elsevier Inc. All rights reserved.

Developing pineapple fruit has a small transcriptome dominated by metallothionein

In a first step toward understanding the molecular basis of pineapple fruit development, a sequencing project was initiated to survey a range of expressed sequences from green unripe and yellow ripe fruit tissue. A highly abundant metallothionein transcript was identified during library construction, and was estimated to account for up to 50% of all EST library clones. Library clones with metallothionein subtracted were sequenced, and 408 unripe green and 1140 ripe yellow edited EST clone sequences were retrieved. Clone redundancy was high, with the combined 1548 clone sequences clustering into just 634 contigs comprising 191 consensus sequences and 443 singletons. Half of the EST clone sequences clustered within 13.5% and 9.3% of contigs from green unripe and yellow ripe libraries, respectively, indicating that a small subset of genes dominate the majority of the transcriptome. Furthermore, sequence cluster analysis, northern analysis, and functional classification revealed major differences between genes expressed in the unripe green and ripe yellow fruit tissues. Abundant genes identified from the green fruit include a fruit bromelain and a bromelain inhibitor. Abundant genes identified in the yellow fruit library include a MADS box gene, and several genes normally associated with protein synthesis, including homologues of ribosomal L10 and the translation factors SUI1 and eIF5A. Both the green unripe and yellow ripe libraries contained high proportions of clones associated with oxidative stress responses and the detoxification of free radicals.

Length suppression in histone messenger RNA 3′-end maturation: Processing defects of insertion mutant premessenger RNAs can be compensated by insertions into the U7 small nuclear RNA

Efficient 3′-end processing of cell cycle-regulated mammalian histone premessenger RNAs (pre-mRNAs) requires an upstream stem–loop and a histone downstream element (HDE) that base pairs with the U7 small ribonuclearprotein. Insertions between these elements have two effects: the site of cleavage moves in concert with the HDE and processing efficiency declines. We used Xenopus oocytes to ask whether compensatory length insertions in the human U7 RNA could restore the fidelity and efficiency of processing of mouse histone insertion pre-mRNAs. An insertion of 5 nt into U7 RNA that extends its complementary to the HDE compensated for both defects in processing of a 5-nt insertion substrate; a noncomplementary insertion into U7 did not. Yet, the noncomplementary insertion mutant U7 was shown to be active on insertion substrates further mutated to allow base pairing. Our results suggest that the histone pre-mRNA becomes rigidified upstream of its HDE, allowing the bound U7 small ribonucleoprotein to measure from the HDE to the cleavage site. Such a mechanism may be common to other RNA measuring systems. To our knowledge, this is the first demonstration of length suppression in an RNA processing system.

The evolution of controlled multitasked gene networks: The role of introns and other noncoding RNAs in the development of complex organisms

Eukaryotic phenotypic diversity arises from multitasking of a core proteome of limited size. Multitasking is routine in computers, as well as in other sophisticated information systems, and requires multiple inputs and outputs to control and integrate network activity. Higher eukaryotes have a mosaic gene structure with a dual output, mRNA (protein-coding) sequences and introns, which are released from the pre-mRNA by posttranscriptional processing. Introns have been enormously successful as a class of sequences and comprise up to 95% of the primary transcripts of protein-coding genes in mammals. In addition, many other transcripts (perhaps more than half) do not encode proteins at all, but appear both to be developmentally regulated and to have genetic function. We suggest that these RNAs (eRNAs) have evolved to function as endogenous network control molecules which enable direct gene-gene communication and multitasking of eukaryotic genomes. Analysis of a range of complex genetic phenomena in which RNA is involved or implicated, including co-suppression, transgene silencing, RNA interference, imprinting, methylation, and transvection, suggests that a higher-order regulatory system based on RNA signals operates in the higher eukaryotes and involves chromatin remodeling as well as other RNA-DNA, RNA-RNA, and RNA-protein interactions. The evolution of densely connected gene networks would be expected to result in a relatively stable core proteome due to the multiple reuse of components, implying,that cellular differentiation and phenotypic variation in the higher eukaryotes results primarily from variation in the control architecture. Thus, network integration and multitasking using trans-acting RNA molecules produced in parallel with protein-coding sequences may underpin both the evolution of developmentally sophisticated multicellular organisms and the rapid expansion of phenotypic complexity into uncontested environments such as those initiated in the Cambrian radiation and those seen after major extinction events.

Structured RNAs in the ENCODE selected regions of the human genome.

Functional RNA structures play an important role both in the context of noncoding RNA transcripts as well as regulatory elements in mRNAs. Here we present a computational study to detect functional RNA structures within the ENCODE regions of the human genome. Since structural RNAs in general lack characteristic signals in primary sequence, comparative approaches evaluating evolutionary conservation of structures are most promising. We have used three recently introduced programs based on either phylogenetic-stochastic context-free grammar (EvoFold) or energy directed folding (RNAz and AlifoldZ), yielding several thousand candidate structures (corresponding to approximately 2.7% of the ENCODE regions). EvoFold has its highest sensitivity in highly conserved and relatively AU-rich regions, while RNAz favors slightly GC-rich regions, resulting in a relatively small overlap between methods. Comparison with the GENCODE annotation points to functional RNAs in all genomic contexts, with a slightly increased density in 3'-UTRs. While we estimate a significant false discovery rate of approximately 50%-70% many of the predictions can be further substantiated by additional criteria: 248 loci are predicted by both RNAz and EvoFold, and an additional 239 RNAz or EvoFold predictions are supported by the (more stringent) AlifoldZ algorithm. Five hundred seventy RNAz structure predictions fall into regions that show signs of selection pressure also on the sequence level (i.e., conserved elements). More than 700 predictions overlap with noncoding transcripts detected by oligonucleotide tiling arrays. One hundred seventy-five selected candidates were tested by RT-PCR in six tissues, and expression could be verified in 43 cases (24.6%).

Small RNAs as regulators of primary and secondary metabolism in Pseudomonas species.

Small RNAs (sRNAs) exert important functions in pseudomonads. Classical sRNAs comprise the 4.5S, 6S, 10Sa and 10Sb RNAs, which are known in enteric bacteria as part of the signal recognition particle, a regulatory component of RNA polymerase, transfer-messenger RNA (tmRNA) and the RNA component of RNase P, respectively. Their homologues in pseudomonads are presumed to have analogous functions. Other sRNAs of pseudomonads generally have little or no sequence similarity with sRNAs of enteric bacteria. Numerous sRNAs repress or activate the translation of target mRNAs by a base-pairing mechanism. Examples of this group in Pseudomonas aeruginosa are the iron-repressible PrrF1 and PrrF2 sRNAs, which repress the translation of genes encoding iron-containing proteins, and PhrS, an anaerobically inducible sRNA, which activates the expression of PqsR, a regulator of the Pseudomonas quinolone signal. Other sRNAs sequester RNA-binding proteins that act as translational repressors. Examples of this group in P. aeruginosa include RsmY and RsmZ, which are central regulatory elements in the GacS/GacA signal transduction pathway, and CrcZ, which is a key regulator in the CbrA/CbrB signal transduction pathway. These pathways largely control the extracellular activities (including virulence traits) and the selection of the energetically most favourable carbon sources, respectively, in pseudomonads.