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.

Identification of putative noncoding RNAs among the RIKEN mouse full-length cDNA collection

With the sequencing and annotation of genomes and transcriptomes of several eukaryotes, the importance of noncoding RNA (ncRNA)-RNA molecules that are not translated to protein products-has become more evident. A subclass of ncRNA transcripts are encoded by highly regulated, multi-exon, transcriptional units, are processed like typical protein-coding mRNAs and are increasingly implicated in regulation of many cellular functions in eukaryotes. This study describes the identification of candidate functional ncRNAs from among the RIKEN mouse full-length cDNA collection, which contains 60,770 sequences, by using a systematic computational filtering approach. We initially searched for previously reported ncRNAs and found nine murine ncRNAs and homologs of several previously described nonmouse ncRNAs. Through our computational approach to filter artifact-free clones that lack protein coding potential, we extracted 4280 transcripts as the largest-candidate set. Many clones in the set had EST hits, potential CpG islands surrounding the transcription start sites, and homologies with the human genome. This implies that many candidates are indeed transcribed in a regulated manner. Our results demonstrate that ncRNAs are a major functional subclass of processed transcripts in mammals.

Growth hormone as a cytokine

1. The growth hormone (GH) receptor was the first of the class 1 cytokine receptors to be cloned. It shares a number of structural characteristics with other family members and common signalling mechanisms based on common usage of the Janus kinase 2 (JAK2). 2. Growth hormone receptor activation is initiated by GH-induced homodimerization of receptor molecules. This has enabled the creation of specific hormone antagonists that block receptor dimerization. 3. The details of the transcription factors used by the activated receptor are being revealed as a result of promoter analyses and electrophoretic mobility gelshift analysis. 4. Growth hormone receptors are widespread and their discovery in certain tissues has led to the assignment of new physiological roles for GH, Some of these involve local or paracrine roles for GH, as befits its cytokine status. 5. Four examples of such novel roles are discussed, These are: (i) the brain GH axis; (ii) GH and the vitamin B-12 axis; (iii) GH in early pre-implantation development; and (iv) GH in development of the tooth. 6. We propose that the view that GH acts through the intermediacy of insulin-like growth factor-1 is simplistic; rather, GH acts to induce an array of growth factors and their receptors and the composition of this array varies with tissue type and, probably, stage of development.

Model for growth hormone receptor activation based on subunit rotation within a receptor dimer

Growth hormone is believed to activate the growth hormone receptor (GHR) by dimerizing two identical receptor subunits, leading to activation of JAK2 kinase associated with the cytoplasmic domain. However, we have reported previously that dimerization alone is insufficient to activate full-length GHR. By comparing the crystal structure of the liganded and unliganded human GHR extracellular domain, we show here that there is no substantial change in its conformation on ligand binding. However, the receptor can be activated by rotation without ligand by inserting a defined number of alanine residues within the transmembrane domain. Fluorescence resonance energy transfer ( FRET), bioluminescence resonance energy transfer (BRET) and coimmunoprecipitation studies suggest that receptor subunits undergo specific transmembrane interactions independent of hormone binding. We propose an activation mechanism involving a relative rotation of subunits within a dimeric receptor as a result of asymmetric placement of the receptor-binding sites on the ligand.

Expression of the CD44v2-10 isoform confers a metastatic phenotype: Importance of the heparan sulfate attachment site CD44v3

We expressed the full-length CD44v2-10 isoform in SKHep1 cells, a nonmetastatic human hepatocellular carcinoma cell line that does not express any endogenous CD44v isoforms. In SCID mice, expression of CD44v2-10 by SKHep1 cells had no effect on s.c. primary tumor development but caused pulmonary metastases in 41% (7 of 17) of animals compared with control SKHep1 cells (0 of 16; P < 0.01). CD44v2-10 expression by SKHep1 cells resulted in enhanced heparan sulfate (HS) attachment and an enhanced capacity to bind heparin-binding growth factors. Mutation of the v3 domain to prevent HS attachment and growth factor binding abolished the metastatic phenotype, demonstrating that HS modification of CD44v2-10 plays a critical role in the development of metastases in this model. However, in vitro proliferation, motility, and invasion were not altered by CD44v2-10 expression.

Sorting nexin 5 is localized to a subdomain of the early endosomes and is recruited to the plasma imembrane following EGF stimulation

Sorting nexins are a large family of proteins that contain the phosphoinositide-binding Phox homology (PX) domain. A number of sorting nexins are known to bind to PtdIns(3)P, which mediates their localization to membranes of the endocytic pathway. We show here that sorting nexin 5 (SNX5) can be recruited to two distinct membrane compartments. In non-stimulated cells, the PX domain was independently targeted to endosomal structures and colocalized with full-length SNX5. The membrane binding of the PX domain was inhibited by the PI 3-kinase inhibitor, wortmannin. Although SNX5 colocalized with a fluid-phase marker and was found predominantly within a PtdIns(3)P-rich endosomal domain, very little colocalization was observed between SNX5 and the PtdIns(3)P-binding protein, EEA1. Using liposome-based binding assays, we have shown that the PX domain of SNX5 interacts not only with PtdIns(3)P but also with PtdIns(3,4)P-2. In response to EGF stimulation, either the SNX5-PX domain or full-length SNX5 was rapidly recruited to the plasma membrane. The localization of SNX1, which does not bind PtdIns(3,4)P-2, was unaffected by EGF signalling. Therefore, SNX5 is localized to a subdomain of the early endosome distinct from EEA1 and, following EGF stimulation and elevation of PtdIns(3,4)P-2, is also transiently recruited to the plasma membrane. These results indicate that SNX5 may have functions not only associated with endosomal sorting but also with the phosphoinositide-signalling pathway.

Cloning and functional expression of venom prothrombin activator protease from Pseudonaja textilis with whole blood procoagulant activity

The snake venom group C prothrombin activators contain a number of components that enhance the rate of prothrombin activation. The cloning and expression of full-length cDNA for one of these components, an activated factor X (factor Xa)-like protease from Pseudonaja textilis as well as the generation of functional chimeric constructs with procoagulant activity were described. The complete cDNA codes for a propeptide, light chain, activation peptide (AP) and heavy chain related in sequence to mammalian factor X. Efficient expression of the protease was achieved with constructs where the AP was deleted and the cleavage sites between the heavy and light chains modified, or where the AP was replaced with a peptide involved in insulin receptor processing. In human kidney cells (H293F) transfected with these constructs, up to 80% of the pro-form was processed to heavy and light chains. Binding of the protease to barium citrate and use of specific antibodies demonstrated that gamma-carboxylation of glutamic acid residues had occurred on the light chain in both cases, as observed in human factor Xa and the native P. textilis protease. The recombinant protease caused efficient coagulation of whole citrated blood and citrated plasma that was enhanced by the presence of Ca2+. This study identified the complete cDNA sequence of a factor Xa-like protease from P. textilis and demonstrated for the first time the expression of a recombinant form of P. textilis protease capable of blood coagulation.

Identification and analysis of chromodomain-containing proteins encoded in the mouse transcriptome

The chromodomain is 40-50 amino acids in length and is conserved in a wide range of chromatic and regulatory proteins involved in chromatin remodeling. Chromodomain-containing proteins can be classified into families based on their broader characteristics, in particular the presence of other types of domains, and which correlate with different subclasses of the chromodomains themselves. Hidden Markov model (HMM)-generated profiles of different subclasses of chromodomains were used here to identify sequences encoding chromodomain-containing proteins in the mouse transcriptome and genome. A total of 36 different loci encoding proteins containing chromodomains, including 17 novel loci, were identified. Six of these loci (including three apparent pseudogenes, a novel HP1 ortholog, and two novel Msl-3 transcription factor-like proteins) are not present in the human genome, whereas the human genome contains four loci (two CDY orthologs and two apparent CDY pseuclogenes) that are not present in mouse. A number of these loci exhibit alternative splicing to produce different isoforms, including 43 novel variants, some of which lack the chromodomain. The likely functions of these proteins are discussed in relation to the known functions of other chromodomain-containing proteins within the same family.

Dynamics of motor nerve terminal remodeling unveiled using SNARE-cleaving botulinum toxins: the extent and duration are dictated by the sites of SNAP-25 truncation

Nerve sprouts emerge from motor nerve terminals following blockade of exo-endocytosis for more than 3 days by botulinum neurotoxin (BoNT), and form functional synapses, albeit temporary. Upon restoration of synaptic activity to the parent terminal 7 and 90 days after exposure to BoNT/F or A respectively, a concomitant retraction of the outgrowths was observed. BoNT/E caused short-term neuroparalysis, and dramatically accelerated the recovery of BoNT/A-paralyzed muscle by further truncation of SNAP-25 and its replenishment with functional full-length SNARE. The removal of 9 C-terminal residues from SNAP-25 by BoNT/A leads to persistence of the inhibitory product due to the formation of a nonproductive SNARE complex(es) at release sites, whereas deletion of a further 17 amino acids permits replenishment and a speedy recovery. (C) 2003 Elsevier Science (USA). All rights reserved.

Development and evaluation of an automated annotation pipeline and cDNA annotation system

Manual curation has long been held to be the gold standard for functional annotation of DNA sequence. Our experience with the annotation of more than 20,000 full-length cDNA sequences revealed problems with this approach, including inaccurate and inconsistent assignment of gene names, as well as many good assignments that were difficult to reproduce using only computational methods. For the FANTOM2 annotation of more than 60,000 cDNA clones, we developed a number of methods and tools to circumvent some of these problems, including an automated annotation pipeline that provides high-quality preliminary annotation for each sequence by introducing an uninformative filter that eliminates uninformative annotations, controlled vocabularies to accurately reflect both the functional assignments and the evidence supporting them, and a highly refined, Web-based manual annotation tool that allows users to view a wide array of sequence analyses and to assign gene names and putative functions using a consistent nomenclature. The ultimate utility of our approach is reflected in the low rate of reassignment of automated assignments by manual curation. Based on these results, we propose a new standard for large-scale annotation, in which the initial automated annotations are manually investigated and then computational methods are iteratively modified and improved based on the results of manual curation.

Systematic characterization of the zinc-finger-containing proteins in the mouse transcriptome

Zinc-finger-containing proteins can be classified into evolutionary and functionally divergent protein families that share one or more domains in which a zinc ion is tetrahedrally coordinated by cysteines and histidines. The zinc finger domain defines one of the largest protein superfamilies in mammalian genomes; 46 different conserved zinc finger domains are listed in InterPro (http://www.ebi.ac.uk/InterPro). Zinc finger proteins can bind to DNA, RNA, other proteins, or lipids as a modular domain in combination with other conserved structures. Owing to this combinatorial diversity, different members of zinc finger superfamilies contribute to many distinct cellular processes, including transcriptional regulation, mRNA stability and processing, and protein turnover. Accordingly, mutations of zinc finger genes lead to aberrations in a broad spectrum of biological processes such as development, differentiation, apoptosis, and immunological responses. This study provides the first comprehensive classification of zinc finger proteins in a mammalian transcriptome. Specific detailed analysis of the SP/Kruppel-like factors and the E3 ubiquitin-ligase RING-H2 families illustrates the importance of such an analysis for a more comprehensive functional classification of large protein families. We describe the characterization of a new family of C2H2 zinc-finger-containing proteins and a new conserved domain characteristic of this family, the identification and characterization of Sp8, a new member of the Sp family of transcriptional regulators, and the identification of five new RING-H2 proteins.

Impact of alternative initiation, splicing, and termination on the diversity of the mRNA transcripts encoded by the mouse transcriptome

We analyzed the FANTOM2 clone set of 60,770 RIKEN full-length mouse cDNA sequences and 44,122 public mRNA sequences. We developed a new computational procedure to identify and classify the forms of splice variation evident in this data set and organized the results into a publicly accessible database that can be used for future expression array construction, structural genomics, and analyses of the mechanism and regulation of alternative splicing. Statistical analysis shows that at least 41% and possibly as much as 60% of multiexon genes in mouse have multiple splice forms. Of the transcription units with multiple splice forms, 49% contain transcripts in which the apparent use of an alternative transcription start (stop) is accompanied by alternative splicing of the initial (terminal) exon. This implies that alternative transcription may frequently induce alternative splicing. The fact that 73% of all exons with splice variation fall within the annotated coding region indicates that most splice variation is likely to affect the protein form. Finally, we compared the set of constitutive (present in all transcripts) exons with the set of cryptic (present only in some transcripts) exons and found statistically significant differences in their length distributions, the nucleoticle distributions around their splice junctions, and the frequencies of occurrence of several short sequence motifs.

Alternative polyadenylation and splicing of mRNAs transcribed from the human Sin1 gene

Limited but significant sequence similarity has been observed between an uncharacterized human protein, SIN1, and the S. pombe SIN1, Dictyostelium RIP3 and S. cerevisiae AVO1 proteins. The human Sin1 gene has been automatically predicted (MAPKAP1; GenBank accession number NM_024117); however, this sequence appears to be incomplete. In this study, we have cloned and characterized the full-length human Sin1 mRNA and identified a highly conserved domain that defines the family of SIN1 orthologues, members of which are widely distributed in the fungal and metazoan kingdoms. We demonstrate that Sin1 transcripts can use alternative polyadenylation signals and describe a number of Sin1 splice variants that potentially encode functionally different isoforms. (C) 2004 Elsevier B.V. All rights reserved.

Endogenous Presentation of CD8+ T Cell Epitopes from Epstein-Barr Virus–encoded Nuclear Antigen 1

Epstein-Barr virus (EBV)-encoded nuclear antigen (EBNA)1 is thought to escape cytotoxic T lymphocyte (CTL) recognition through either self-inhibition of synthesis or by blockade of proteasomal degradation by the glycine-alanine repeat (GAr) domain. Here we show that EBNA1 has a remarkably varied cell type-dependent stability. However, these different degradation rates do not correspond to the level of major histocompatibility complex class I-restricted presentation of EBNA1 epitopes. In spite of the highly stable expression of EBNA1 in B cells, CTL epitopes derived from this protein are efficiently processed and presented to CD8(+) T cells. Furthermore, we show that EBV-infected B cells can readily activate EBNA1-specific memory T cell responses from healthy virus carriers. Functional assays revealed that processing of these EBNA1 epitopes is proteasome and transporter associated with antigen processing dependent. We also show that the endogenous presentation of these epitopes is dependent on the newly synthesized protein rather than the long-lived stable EBNA1. Based on these observations, we propose that defective ribosomal products, not the full-length antigen, are the primary source of endogenously processed CD8(+) T cell epitopes front EBNA1.

Aprataxin, a novel protein that protects against genotoxic stress

Ataxia-oculomotor apraxia (AOA1) is a neurological disorder with symptoms that overlap those of ataxia-telangiectasia, a syndrome characterized by abnormal responses to double-strand DNA breaks and genome instability. The gene mutated in AOA1, APTX, is predicted to code for a protein called aprataxin that contains domains of homology with proteins involved in DNA damage signalling and repair. We demonstrate that aprataxin is a nuclear protein, present in both the nucleoplasm and the nucleolus. Mutations in the APTX gene destabilize the aprataxin protein, and fusion constructs of enhanced green fluorescent protein and aprataxin, representing deletions of putative functional domains, generate highly unstable products. Cells from AOA1 patients are characterized by enhanced sensitivity to agents that cause single-strand breaks in DNA but there is no evidence for a gross defect in single-strand break repair. Sensitivity to hydrogen peroxide and the resulting genome instability are corrected by transfection with full-length aprataxin cDNA. We also demonstrate that aprataxin interacts with the repair proteins XRCC1, PARP-1 and p53 and that it co-localizes with XRCC1 along charged particle tracks on chromatin. These results demonstrate that aprataxin influences the cellular response to genotoxic stress very likely by its capacity to interact with a number of proteins involved in DNA repair.