Sequencing and Analysis of Approximately 40 000 Soybean cDNA Clones from a Full-Length-Enriched cDNA Library

A large collection of full-length cDNAs is essential for the correct annotation of genomic sequences and for the functional analysis of genes and their products. We obtained a total of 39 936 soybean cDNA clones (GMFL01 and GMFL02 clone sets) in a full-length-enriched cDNA library which was constructed from soybean plants that were grown under various developmental and environmental conditions. Sequencing from 5′ and 3′ ends of the clones generated 68 661 expressed sequence tags (ESTs). The EST sequences were clustered into 22 674 scaffolds involving 2580 full-length sequences. In addition, we sequenced 4712 full-length cDNAs. After removing overlaps, we obtained 6570 new full-length sequences of soybean cDNAs so far. Our data indicated that 87.7% of the soybean cDNA clones contain complete coding sequences in addition to 5′- and 3′-untranslated regions. All of the obtained data confirmed that our collection of soybean full-length cDNAs covers a wide variety of genes. Comparative analysis between the derived sequences from soybean and Arabidopsis, rice or other legumes data revealed that some specific genes were involved in our collection and a large part of them could be annotated to unknown functions. A large set of soybean full-length cDNA clones reported in this study will serve as a useful resource for gene discovery from soybean and will also aid a precise annotation of the soybean genome.Key words: EST, full-length cDNA, functional annotation, legume, soybean     

Construction of expression-ready cDNA clones for KIAA genes: manual curation of 330 KIAA cDNA clones.

We have accumulated information on protein-coding sequences of uncharacterized human genes, which are known as KIAA genes, through cDNA sequencing. For comprehensive functional analysis of the KIAA genes, it is necessary to prepare a set of cDNA clones which direct the synthesis of functional KIAA gene products. However, since the KIAA cDNAs were derived from long mRNAs (> 4 kb), it was not expected that all of them were full-length. Thus, as the first step toward preparing these clones, we evaluated the integrity of protein-coding sequences of KIAA cDNA clones through comparison with homologous protein entries in the public database. As a result, 1141 KIAA cDNAs had at least one homologous entry in the database, and 619 of them (54%) were found to be truncated at the 5' and/or 3' ends. In this study, 290 KIAA cDNA clones were tailored to be full-length or have considerably longer sequences than the original clones by isolating additional cDNA clones and/or connected parts of additional cDNAs or PCR products of the missing portion to the original cDNA clone. Consequently, 265, 8, and 17 predicted CDSs of KIAA cDNA clones were increased in the amino-, carboxy-, and both terminal sequences, respectively. In addition, 40 cDNA clones were modified to remove spurious interruption of protein-coding sequences. The total length of the resultant extensions at amino- and carboxy-terminals of KIAA gene products reached 97,000 and 7,216 amino acid residues, respectively, and various protein domains were found in these extended portions.     

Cloning and functional annotation of rare mRNA species from drought-stressed hot pepper (<Emphasis Type="Italic">Capsicum annuum</Emphasis>)

To better understand gene expression at very low levels, we have designed a method to eliminate cDNA clones representing abundant mRNAs. A cDNA library for drought-stressed hot pepper (Capsicum annuum) (Choi et al., 2002) underwent double-negative screening, once with probes made from a drought-stressed plant, the second time, with probes from a non-stressed plant. The cDNA clones that showed very weak or negative signals were isolated for further analysis, which resulted in 1399 cDNA clones from about 20,000 screened clones. When nucleotide sequences were determined, we obtained 1142 tentative unique genes, with a redundancy rate of 20.41%. An homology database search for the deduced amino acid sequences revealed that about 79% of the cDNA clones could not be matched for functioning with previously characterized sequences. However, when these uncategorized clones were subjected to classification based on functional domains, most could be cited. Notably, clones with possible functions in RNA transport, protein synthesis, and regulation of protein activity showed a dramatic increase in appearance while those coding for transposable elements, viral proteins, and plasmid proteins occupied a much smaller portion compared with those in theArabidopsis thaliana genome. In addition, those coding for proteins targeted to the endoplasmic reticulum were dramatically more abundant in our clones compared with theArabidopsis database.     

Development of Full-Length cDNAs from Chinese Cabbage (Brassica rapa Subsp. pekinensis) and Identification of Marker Genes for Defence Response

Arabidopsis belongs to the Brassicaceae family and plays an important role as a model plant for which researchers have developed fine-tuned genome resources. Genome sequencing projects have been initiated for other members of the Brassicaceae family. Among these projects, research on Chinese cabbage (Brassica rapa subsp. pekinensis) started early because of strong interest in this species. Here, we report the development of a library of Chinese cabbage full-length cDNA clones, the RIKEN BRC B. rapa full-length cDNA (BBRAF) resource, to accelerate research on Brassica species. We sequenced 10 000 BBRAF clones and confirmed 5476 independent clones. Most of these cDNAs showed high homology to Arabidopsis genes, but we also obtained more than 200 cDNA clones that lacked any sequence homology to Arabidopsis genes. We also successfully identified several possible candidate marker genes for plant defence responses from our analysis of the expression of the Brassica counterparts of Arabidopsis marker genes in response to salicylic acid and jasmonic acid. We compared gene expression of these markers in several Chinese cabbage cultivars. Our BBRAF cDNA resource will be publicly available from the RIKEN Bioresource Center and will help researchers to transfer Arabidopsis-related knowledge to Brassica crops.     

High-efficiency cloning of Arabidopsis full-length cDNA by biotinylated CAP trapper

Full-length cDNAs are essential for functional analysis of plant genes. We constructed high-content, full-length cDNA libraries from Arabidopsis thaliana plants based on chemical introduction of a biotin group into the diol residue of the CAP structure of eukaryotic mRNA, followed by RNase I treatment, to select full-length cDNA. More than 90% of the total clones obtained were of full length; recombinant clones were obtained with high efficiency (2.2 × 10⁶/9 μg starting mRNA). Sequence analysis of 111 randomly picked clones indicated that 32 isolated cDNA groups were derived from novel genes in the A. thaliana genome.     

Molecular characterization of multiple cDNA clones for ADP-glucose pyrophosphorylase from Arabidopsis thaliana

PCR amplification of cDNA prepared from poly(A)⁺ RNA from aerial parts of Arabidopsis thaliana, using degenerate nucleotide primers based on conserved regions between the large and small subunits of ADP-glucose pyrophosphorylase (AGP), yielded four different cDNAs of ca. 550 nucleotides each. Based on derived amino acid sequences, the identities between the clones varied from 49 to 69%. Sequence comparison to previously published cDNAs for AGP from various species and tissues has revealed that three of the amplified cDNAs (ApL1, ApL2 and ApL3) correspond to the large subunit of AGP, and one cDNA (ApS) encodes the small subunit of AGP. Both ApL1 and ApS were subsequently found to be present in a cDNA library made from Arabidopsis leaves. All four PCR products are encoded by single genes, as found by genomic Southern analysis.     

Efficient and high-throughput vector construction and Agrobacterium-mediated transformation of Arabidopsis thaliana suspension-cultured cells for functional genomics

We established a large-scale, high-throughput protocol to construct Arabidopsis thaliana suspension-cultured cell lines, each of which carries a single transgene, using Agrobacterium-mediated transformation. We took advantage of RIKEN Arabidopsis full-length (RAFL) cDNA clones and the Gateway cloning system for high-throughput preparation of binary vectors carrying individual full-length cDNA sequences. Throughout all cloning steps, multiple-well plates were used to treat 96 samples simultaneously in a high-throughput manner. The optimal conditions for Agrobacterium-mediated transformation of 96 independent binary vector constructs were established to obtain transgenic cell lines efficiently. We evaluated the protocol by generating transgenic Arabidopsis T87 cell lines carrying individual 96 metabolism-related RAFL cDNA fragments, and showed that the protocol was useful for high-throughput and large-scale production of gain-of-function lines for functional genomics.     

A genome-wide gain-of-function analysis of rice genes using the FOX-hunting system

The latest report has estimated the number of rice genes to be approximately 32,000. To elucidate the functions of a large population of rice genes and to search efficiently for agriculturally useful genes, we have been taking advantage of the Full-length cDNA Over-eXpresser (FOX) gene-hunting system. This system is very useful for analyzing various gain-of-function phenotypes from large populations of transgenic plants overexpressing cDNAs of interest and others with unknown or important functions. We collected the plasmid DNAs of 13,980 independent full-length cDNA (FL-cDNA) clones to produce a FOX library by placing individual cDNAs under the control of the maize Ubiquitin-1 promoter. The FOX library was transformed into rice by Agrobacterium-mediated high-speed transformation. So far, we have generated approximately 12,000 FOX-rice lines. Genomic PCR analysis indicated that the average number of FL-cDNAs introduced into individual lines was 1.04. Sequencing analysis of the PCR fragments carrying FL-cDNAs from 8615 FOX-rice lines identified FL-cDNAs in 8225 lines, and a database search classified the cDNAs into 5462 independent ones. Approximately 16.6% of FOX-rice lines examined showed altered growth or morphological characteristics. Three super-dwarf mutants overexpressed a novel gibberellin 2-oxidase gene,confirming the importance of this system. We also show here the other morphological alterations caused by individual FL-cDNA expression. These dominant phenotypes should be valuable indicators for gene discovery and functional analysis.     

Development of 5006 Full-Length CDNAs in Barley: A Tool for Accessing Cereal Genomics Resources

A collection of 5006 full-length (FL) cDNA sequences was developed in barley. Fifteen mRNA samples from various organs and treatments were pooled to develop a cDNA library using the CAP trapper method. More than 60% of the clones were confirmed to have complete coding sequences, based on comparison with rice amino acid and UniProt sequences. Blastn homologies (E<1E-5) to rice genes and Arabidopsis genes were 89 and 47%, respectively. Of the 5028 possible amino acid sequences derived from the 5006 FLcDNAs, 4032 (80.2%) were classified into 1678 GreenPhyl multigenic families. There were 555 cDNAs showing low homology to both rice and Arabidopsis. Gene ontology annotation by InterProScan indicated that many of these cDNAs (71%) have no known molecular functions and may be unique to barley. The cDNAs showed high homology to Barley 1 GeneChip oligo probes (81%) and the wheat gene index (84%). The high homology between FLcDNAs (27%) and mapped barley expressed sequence tag enabled assigning linkage map positions to 151–233 FLcDNAs on each of the seven barley chromosomes. These comprehensive barley FLcDNAs provide strong platform to connect pre-existing genomic and genetic resources and accelerate gene identification and genome analysis in barley and related species.Key words: full-length cDNA, Hordeum vulgare, mRNA, gene ontology     

A General Bacterial Expression System for Functional Analysis of cDNA-Encoded Proteins

A general system for functional analysis of cDNA-encoded proteins is described. The basic concept involves the expression inEscherichia coliof selected portions of cDNAs in an approach toward the understanding of the function of the corresponding proteins. A selected cDNA is expressed as part of a fusion protein used for immunization to elicit antibodies, and a corresponding fusion protein, having the cDNA-encoded portion in common, for purification of target protein-specific antibodies. This antiserum could be used for functional analysis of the cDNA-encoded protein, e.g., by immunohistology. Two general expression vector systems forE. colihave been constructed, both (i) designed with multiple cloning sites in three different reading frames, (ii) having their protein production controlled by the tightly regulated T7 promoter, and (iii) enabling affinity purification of the expressed target proteins by fusions to IgG-binding domains derived from staphylococcal protein A or a serum albumin-binding protein derived from streptococcal protein G, respectively. This novel system has been evaluated by expressing five cDNAs, isolated from pre- pubertal mouse testis by a differential cDNA library screening strategy. All five clones could be expressed intracellularly inE. colias fusion proteins with high production levels, ranging from 4 to 500 mg/liter, and affinity purification yielded essentially full-length products. Characterization of affinity-purified antibodies revealed that there exists no cross-reactivity between the two fusion systems and that such antibodies indeed could be used for immunohistology. The implications for the described system for large-scale functional analysis of cDNA libraries are discussed.     

Cloning of cDNAs for genes that are early-responsive to dehydration stress (ERDs) inArabidopsis thaliana L.: identification of three ERDs as HSP cognate genes

InArabidopsis thaliana L., accumulation of abscisic acid (ABA) began to increase 2 h after plants had been subjected to dehydration stress and reached maximum levels after 10h. Differential hybridization was used to isolate 26Arabidopsis cDNAs with gene expression induced by a 1 h dehydration treatment. The cDNA clones were classified into 16 groups based on Southern blot hybridization, and named ERD (early-responsive todehydration) clones. Partial sequencing of the cDNA clones revealed that three ERDs were identical to those of HSP cognates (Athsp70-1, Athsp81-2, and ubiquitin extension protein). Dehydration stress strongly induced the expression of genes for the three ERDs, while application of ABA, which is known to act as a signal transmitter in dehydration-stressed plants, did not significantly affect the ERD gene expression. This result suggests that these HSP cognates are preferentially responsive to dehydration stress inA. thaliana, and that signaling pathways for the expression of these genes under conditions of dehydration stress are not mainly mediated by ABA. We also discuss the possible functions of these three ERD gene products against dehydration stress.