Takahiro Arakawa

Functional annotation of a full-length mouse cDNA collection

The RIKEN Mouse Gene Encyclopaedia Project, a systematic approach to determining the full coding potential of the mouse genome, involves collection and sequencing of full-length complementary DNAs and physical mapping of the corresponding genes to the mouse genome. We organized an international functional annotation meeting (FANTOM) to annotate the first 21,076 cDNAs to be analysed in this project. Here we describe the first RIKEN clone collection, which is one of the largest described for any organism. Analysis of these cDNAs extends known gene families and identifies new ones.The RIKEN Mouse Gene Encyclopaedia Project, a systematic approach to determining the full coding potential of the mouse genome, involves collection and sequencing of full-length complementary DNAs and physical mapping of the corresponding genes to the mouse genome. We organized an international functional annotation meeting (FANTOM) to annotate the first 21,076 cDNAs to be analysed in this project. Here we describe the first RIKEN clone collection, which is one of the largest described for any organism. Analysis of these cDNAs extends known gene families and identifies new ones.

Cap analysis gene expression for high-throughput analysis of transcriptional starting point and identification of promoter usage

We introduce cap analysis gene expression (CAGE), which is based on preparation and sequencing of concatamers of DNA tags deriving from the initial 20 nucleotides from 5′ end mRNAs. CAGE allows high-throughout gene expression analysis and the profiling of transcriptional start points (TSP), including promoter usage analysis. By analyzing four libraries (brain, cortex, hippocampus, and cerebellum), we redefined more accurately the TSPs of 11-27% of the analyzed transcriptional units that were hit. The frequency of CAGE tags correlates well with results from other analyses, such as serial analysis of gene expression, and furthermore maps the TSPs more accurately, including in tissue-specific cases. The high-throughput nature of this technology paves the way for understanding gene networks via correlation of promoter usage and gene transcriptional factor expression.

The regulated retrotransposon transcriptome of mammalian cells

Although repetitive elements pervade mammalian genomes, their overall contribution to transcriptional activity is poorly defined. Here, as part of the FANTOM4 project, we report that 6–30% of cap-selected mouse and human RNA transcripts initiate within repetitive elements. Analysis of approximately 250,000 retrotransposon-derived transcription start sites shows that the associated transcripts are generally tissue specific, coincide with gene-dense regions and form pronounced clusters when aligned to full-length retrotransposon sequences. Retrotransposons located immediately 5′ of protein-coding loci frequently function as alternative promoters and/or express noncoding RNAs. More than a quarter of RefSeqs possess a retrotransposon in their 3′ UTR, with strong evidence for the reduced expression of these transcripts relative to retrotransposon-free transcripts. Finally, a genome-wide screen identifies 23,000 candidate regulatory regions derived from retrotransposons, in addition to more than 2,000 examples of bidirectional transcription. We conclude that retrotransposon transcription has a key influence upon the transcriptional output of the mammalian genome.

Tiny RNAs associated with transcription start sites in animals

It has been reported that relatively short RNAs of heterogeneous sizes are derived from sequences near the promoters of eukaryotic genes. As part of the FANTOM4 project, we have identified tiny RNAs with a modal length of 18 nt that map within −60 to +120 nt of transcription start sites (TSSs) in human, chicken and Drosophila. These transcription initiation RNAs (tiRNAs) are derived from sequences on the same strand as the TSS and are preferentially associated with G+C-rich promoters. The 5′ ends of tiRNAs show peak density 10–30 nt downstream of TSSs, indicating that they are processed. tiRNAs are generally, although not exclusively, associated with highly expressed transcripts and sites of RNA polymerase II binding. We suggest that tiRNAs may be a general feature of transcription in metazoa and possibly all eukaryotes.

Rapid SNP diagnostics using asymmetric isothermal amplification and a new mismatch-suppression technology

We developed a rapid single nucleotide polymorphism (SNP) detection system named smart amplification process version 2 (SMAP 2). Because DNA amplification only occurred with a perfect primer match, amplification alone was sufficient to identify the target allele. To achieve the requisite fidelity to support this claim, we used two new and complementary approaches to suppress exponential background DNA amplification that resulted from mispriming events. SMAP 2 is isothermal and achieved SNP detection from whole human blood in 30 min when performed with a new DNA polymerase that was cloned and isolated from Alicyclobacillus acidocaldarius (Aac pol). Furthermore, to assist the scientific community in configuring SMAP 2 assays, we developed software specific for SMAP 2 primer design. With these new tools, a high-precision and rapid DNA amplification technology becomes available to aid in pharmacogenomic research and molecular-diagnostics applications.

Rapid Detection of Epidermal Growth Factor Receptor Mutations in Lung Cancer by the SMart-Amplification Process

Purpose: A positive response to gefitinib in non–small cell lung cancer (NSCLC) has been correlated to mutations in epidermal growth factor receptor (EGFR) gene. Previous reports have been based mainly on diagnostic screening by sequencing. However, sequencing is a time-consuming and complicated procedure, not suitable for routine clinical use. Experimental Design: We have developed rapid, simple, and sensitive mutation detection assays based on the SMart Amplification Process (SMAP) and applied it for analyzing EGFR gene mutations in clinical samples. By using SMAP, we can detect mutations within 30 min including sample preparation. To validate the assay system for potential use in clinical diagnostics, we examined 45 NSCLC patients for EGFR mutations using sequencing and SMAP. Results: The outcomes of the SMAP assay perfectly matched the sequencing results, except in one case where SMAP was able to identify a mutation that was not detected by sequencing. We also evaluated the sensitivity and specificity of SMAP in mutation detection for EGFR. In a serial dilution study, SMAP was able to find a mutation in a sample containing only 0.1% of the mutant allele in a mixture of wild-type genomic DNA. We also could show amplification of mutated DNA with only 30 copies per reaction. Conclusions: The SMAP method offers higher sensitivity and specificity than alternative technologies, while eliminating the need for sequencing to identify mutations in the EGFR gene of NSCLC. It provides a robust and point-of-care accessible approach for a rapid identification of most patients likely to respond to gefitinib.

Solution Structure of the SEA Domain from the Murine Homologue of Ovarian Cancer Antigen CA125 (MUC16)

Human CA125, encoded by the MUC16 gene, is an ovarian cancer antigen widely used for a serum assay. Its extracellular region consists of tandem repeats of SEA domains. In this study we determined the three-dimensional structure of the SEA domain from the murine MUC16 homologue using multidimensional NMR spectroscopy. The domain forms a unique α/β sandwich fold composed of two α helices and four antiparallel β strands and has a characteristic turn named the TY-turn between α1 and α2. The internal mobility of the main chain is low throughout the domain. The residues that form the hydrophobic core and the TY-turn are fully conserved in all SEA domain sequences, indicating that the fold is common in the family. Interestingly, no other residues are conserved throughout the family. Thus, the sequence alignment of the SEA domain family was refined on the basis of the three-dimensional structure, which allowed us to classify the SEA domains into several subfamilies. The residues on the surface differ between these subfamilies, suggesting that each subfamily has a different function. In the MUC16 SEA domains, the conserved surface residues, Asn-10, Thr-12, Arg-63, Asp-75, Asp-112, Ser-115, and Phe-117, are clustered on the β sheet surface, which may be functionally important. The putative epitope (residues 58-77) for anti-MUC16 antibodies is located around the β2 and β3 strands. On the other hand the tissue tumor marker MUC1 has a SEA domain belonging to another subfamily, and its GSVVV motif for proteolytic cleavage is located in the short loop connecting β2 and β3.

Solution structure of the RWD domain of the mouse GCN2 protein

GCN2 is the α‐subunit of the only translation initiation factor (eIF2α) kinase that appears in all eukaryotes. Its function requires an interaction with GCN1 via the domain at its N‐terminus, which is termed the RWD domain after three major RWD‐containing proteins: RING finger‐containing proteins, WD‐repeat‐containing proteins, and yeast DEAD (DEXD)‐like helicases. In this study, we determined the solution structure of the mouse GCN2 RWD domain using NMR spectroscopy. The structure forms an α + β sandwich fold consisting of two layers: a four‐stranded antiparallel β‐sheet, and three side‐by‐side α‐helices, with an αββββαα topology. A characteristic YPXXXP motif, which always occurs in RWD domains, forms a stable loop including three consecutive β‐turns that overlap with each other by two residues (triple β‐turn). As putative binding sites with GCN1, a structure‐based alignment allowed the identification of several surface residues in α‐helix 3 that are characteristic of the GCN2 RWD domains. Despite the apparent absence of sequence similarity, the RWD structure significantly resembles that of ubiquitin‐conjugating enzymes (E2s), with most of the structural differences in the region connecting β‐strand 4 and α‐helix 3. The structural architecture, including the triple β‐turn, is fundamentally common among various RWD domains and E2s, but most of the surface residues on the structure vary. Thus, it appears that the RWD domain is a novel structural domain for protein‐binding that plays specific roles in individual RWD‐containing proteins.

Gene organization in rice revealed by full-length cDNA mapping and gene expression analysis through microarray.

Rice (Oryza sativa L.) is a model organism for the functional genomics of monocotyledonous plants since the genome size is considerably smaller than those of other monocotyledonous plants. Although highly accurate genome sequences of indica and japonica rice are available, additional resources such as full-length complementary DNA (FL-cDNA) sequences are also indispensable for comprehensive analyses of gene structure and function. We cross-referenced 28.5K individual loci in the rice genome defined by mapping of 578K FL-cDNA clones with the 56K loci predicted in the TIGR genome assembly. Based on the annotation status and the presence of corresponding cDNA clones, genes were classified into 23K annotated expressed (AE) genes, 33K annotated non-expressed (ANE) genes, and 5.5K non-annotated expressed (NAE) genes. We developed a 60mer oligo-array for analysis of gene expression from each locus. Analysis of gene structures and expression levels revealed that the general features of gene structure and expression of NAE and ANE genes were considerably different from those of AE genes. The results also suggested that the cloning efficiency of rice FL-cDNA is associated with the transcription activity of the corresponding genetic locus, although other factors may also have an effect. Comparison of the coverage of FL-cDNA among gene families suggested that FL-cDNA from genes encoding rice- or eukaryote-specific domains, and those involved in regulatory functions were difficult to produce in bacterial cells. Collectively, these results indicate that rice genes can be divided into distinct groups based on transcription activity and gene structure, and that the coverage bias of FL-cDNA clones exists due to the incompatibility of certain eukaryotic genes in bacteria.

Libraries enriched for alternatively spliced exons reveal splicing patterns in melanocytes and melanomas

It is becoming increasingly clear that alternative splicing enables the complex development and homeostasis of higher organisms. To gain a better understanding of how splicing contributes to regulatory pathways, we have developed an alternative splicing library approach for the identification of alternatively spliced exons and their flanking regions by alternative splicing sequence enriched tags sequencing. Here, we have applied our approach to mouse melan-c melanocyte and B16-F10Y melanoma cell lines, in which 5,401 genes were found to be alternatively spliced. These genes include those encoding important regulatory factors such as cyclin D2, Ilk, MAPK12, MAPK14, RAB4, melastatin 1 and previously unidentified splicing events for 436 genes. Real-time PCR further identified cell line–specific exons for Tmc6, Abi1, Sorbs1, Ndel1 and Snx16. Thus, the ASL approach proved effective in identifying splicing events, which suggest that alternative splicing is important in melanoma development.