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     

Calmodulin isoforms in Arabidopsis encoded by multiple divergent mRNAs

Three new, unique cDNA sequences encoding isoforms of calmodulin (CaM) were isolated from an Arabidopsis cDNA library cloned in gt10. These sequences (ACaM-4, -5, and -6) represent members of the Arabidopsis CaM gene family distinct from the three DNA sequences previously reported. ACaM-4 and -6 encode full-length copies of CaM mRNAs of ca. 0.75 kb. The ACaM-5 sequence encodes a partial length copy of CaM mRNA that is lacking sequences encoding the amino-terminal 10 amino acids of mature CaM and the initiator methionine. The derived amino acid sequence of ACaM-5 is identical to the sequences encoded by two of the previously characterized ACaM cDNAs, and is identical to TCH-1 mRNA, whose accumulation was increased by touch stimulation. The polypeptides encoded by ACaM-4 and -6 differ from that encoded by ACaM-5 by six and two amino acid substititions, respectively. Most of the deduced amino acid sequence substitutions in the Arabidopsis CaM isoforms occurred in the fourth Ca²⁺-binding domain. Polymerase chain reaction amplification assays of ACaM-4, -5 and -6 mRNA sequences indicated that each accumulated in Arabidopsis leaf RNA fractions, but only ACaM-4 and -5 mRNAs were detected in silique total RNA. The six different CaM cDNA sequences each hybridize with unique Eco RI restriction fragments in genomic Southern blots of Arabidopsis DNA, indicating that these sequences were derived from distinct structural genes. Our results suggest that CaM isoforms in Arabidopsis may have evolved to optimize the interaction of this Ca²⁺-receptor protein with specific subsets of response elements.     

新基因水稻 OsLSD1 的克隆及拟南芥和水稻类 LSD1 基因家族的生物信息学分析

类 LSD1 (LSD1-like) 基因家族是一类特殊的 C2C2 型锌指蛋白基因,编码植物特有的转录因子 . 目前已经研究的 2 个成员拟南芥 LSD1 (lesions stimulating disease resistance 1) 和 LOL1 (LSD-One-Like 1) 基因均参与植物细胞程序化死亡 (programmed cell death, PCD) 的调控 . 从水稻 cDNA 文库中克隆到 1 个类 LSD1 基因,命名为 OsLSD1. 该基因长 988 bp ,包含一个 432 bp 的开放阅读框,推导的氨基酸序列 (143 个氨基酸 ) 含有 3 个内部保守的锌指结构域 . DNA 印迹结果表明 OsLSD1 基因在水稻基因组中为单拷贝,且在根、茎和叶中表达 . 借助于生物信息学分析技术,从拟南芥和水稻数据库中各识别出 5 个和 7 个 ( 包括 OsLSD1) 类 LSD1 基因 . 分析了这些类 LSD1 基因的结构,蛋白质结构域组成 . 系统进化分析表明,无论基于编码区的核苷酸或氨基酸序列都可以将这些类 LSD1 基因分为 2 类 . 虽然不存在拟南芥或水稻特有的类 LSD1 蛋白,但有些结构域是水稻所特有的,也有些基因是来源于复制事件 .     

Molecular cloning,characterization and expression analysis of two catalase isozyme genes in barley

Clones representing two distinct barley catalase genes, Cat1 and Cat2, were found in a cDNA library prepared from seedling polysomal mRNA. Both clones were sequenced, and their deduced amino acid sequences were found to have high homology with maize and rice catalase genes. Cat1 had a 91% deduced amino acid sequence identity to CAT-1 of maize and 92% to CAT B of rice. Cat2 had 72 and 79% amino acid sequence identities to maize CAT-2 and-3 and 89% to CAT A of rice. Barley, maize or rice isozymes could be divided into two distinct groups by amino acid homologies, with one group homologous to the mitochondria-associated CAT-3 of maize and the other homologous to the maize peroxisomal/glyoxysomal CAT-1. Both barley CATs contained possible peroxisomal targeting signals, but neither had favorable mitochondrial targeting sequences. Cat1 mRNA occurred in whole endosperms (aleurones plus starchy endosperm), in isolated aleurones and in developing seeds, but Cat2 mRNA was virtually absent. Both mRNAs displayed different developmental expression patterns in scutella of germinating seeds. Cat2 mRNA predominated in etiolated seedling shoots and leaf blades. Barley genomic DNA contained two genes for Cat1 and one gene for Cat2. The Cat2 gene was mapped to the long arm of chromosome 4, 2.9 cM in telomeric orientation from the mlo locus conferring resistance to the powdery mildew fungus (Erysiphe graminis f.sp. hordei).     

Plant aldolase: cDNA and deduced amino-acid sequences of the chloroplast and cytosol enzyme from spinach

We report the sequences of full-length cDNAs for the nuclear genes encoding the chloroplastic and cytosolic fructose-1,6-bisphosphate aldolase (EC 4.1.2.13) from spinach. A comparison of the deduced amino-acid sequences with one another and with published cytosolic aldolase sequences of other plants revealed that the two enzymes from spinach share only 54% homology on their amino acid level whereas the homology of the cytosolic enzyme of spinach with the known sequences of cytosolic aldolases of maize, rice and Arabidopsis range from 67 to 92%. The sequence of the chloroplastic enzyme includes a stroma-targeting N-terminal transit peptide of 46 amino acid residues for import into the chloroplast. The transit peptide exhibits essential features similar to other chloroplast transit peptides. Southern blot analysis implies that both spinach enzymes are encoded by single genes.     

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.     

Comprehensive sequence analysis of 24,783 barley full-length cDNAs derived from 12 clone libraries

Full-length cDNA (FLcDNA) libraries consisting of 172,000 clones were constructed from a two-row malting barley cultivar (Hordeum vulgare 'Haruna Nijo') under normal and stressed conditions. After sequencing the clones from both ends and clustering the sequences, a total of 24,783 complete sequences were produced. By removing duplicates between these and publicly available sequences, 22,651 representative sequences were obtained: 17,773 were novel barley FLcDNAs, and 1,699 were barley specific. Highly conserved genes were found in the barley FLcDNA sequences for 721 of 881 rice (Oryza sativa) trait genes with 50% or greater identity. These FLcDNA resources from our Haruna Nijo cDNA libraries and the full-length sequences of representative clones will improve our understanding of the biological functions of genes in barley, which is the cereal crop with the fourth highest production in the world, and will provide a powerful tool for annotating the barley genome sequences that will become available in the near future.     

Cloning of cDNA Sequences Encoding the Calcium-Binding Protein, Calmodulin, from Barley (Hordeum vulgare L.)

Full- and partial-length cDNAs encoding calmodulin mRNA have been cloned and sequenced from barley (Hordeum vulgare L.). Barley leaf mRNA, size-fractionated in sucrose density gradients, was used to synthesize double-stranded cDNA. The cDNA was cloned in λgt10 and screened with a synthetic, 14-nucleotide oligonucleotide probe, which was designed using the predicted coding sequences of the carboxy termini of spinach and wheat calmodulin proteins. The primary structure of barley calmodulin, predicted from DNA sequencing experiments, consists of 148 amino acids and differs from that of wheat calmodulin in only three positions. In two of the three positions, the amino acid changes are conservative, while the third change consists of an apparent deletion/insertion. The overall nucleotide sequence similarity between the amino acid coding regions of barley and vertebrate calmodulin mRNAs is approximately 77%. However, a region encoding 11 amino acids of the second Ca²⁺-binding domain is very highly conserved at the nucleotide level compared with the rest of the coding sequences (94% sequence identity between barley and chicken calmodulin mRNAs). Genomic Southern blots reveal that barley calmodulin is encoded by a single copy gene. This gene is expressed as a single size class of mRNA in all tissues of 7-day-old barley seedlings. In addition, these analyses indicate that a barley calmodulin cDNA coding region subclone is suitable as a probe for isolating calmodulin genes from other plants.     

Isolation and annotation of 10828 putative full length cDNAs from <Emphasis Type="Italic">indica</Emphasis> rice

We reported the isolation and identification of 10828 putative full-length cDNAs (FL-cDNA) from an indica rice cultivar, Minghui 63, with the long-term goal to isolate all full-length cDNAs from indica genome. Comparison with the databases showed that 780 of them are new rice cDNAs with no match in japonica cDNA database. Totally, 9078 of the FL-cDNAs contained predicted ORFs matching with japonica FL-cDNAs and 6543 could find homologous proteins with complete ORFs. 53% of the matched FL-cDNAs isolated in this study had longer 5′UTR than japonica FL-cDNAs. In silico mapping showed that 9776 (90.28%) of the FL-cDNAs had matched genomic sequences in the japonica genome and 10046 (92.78%) had matched genomic sequences in the indica genome. The average nucleotide sequence identity between the two subspecies is 99.2%. A majority of FL-cDNAs (90%) could be classified with GO (gene ontology) terms based on homology proteins. More than 60% of the new cDNAs isolated in this study had no homology to the known proteins. This set of FL-cDNAs should be useful for functional genomics and proteomics studies.     

A root-specific wall-associated kinase gene, HvWAK1, regulates root growth and is highly divergent in barley and other cereals

Wall-associated receptor-like kinases (WAKs) are important candidates for directly linking the extracellular matrix with intracellular compartments and are involved in developmental processes and stress response. WAK gene family has been identified in plants such as Arabidopsis and rice. Here, we present a detailed analysis of the WAK1 gene from barley cv. Golden Promise, mapped to chromosome 5H. Three BAC clones corresponding to the WAK fragment were sequenced and the full-length WAK1 gene was characterized. The gene has three exons and two short introns with a coding region of 2,178 bp encoding a protein of 725 amino acids. A regulatory region was analyzed in ?1,000 bp sequence upstream to start codon. Using conserved domains database and SMART, various conserved domains such as GUB WAK Bind, epidermal growth factor CA, and protein kinase C as well as other regions like signal peptides, active sites, and transmembrane domains were identified. The gene organization of HvWAK1 was compared with wheat (TaWAK1) and Arabidopsis (AtWAK1), suggesting that the WAK1 gene organization has remained highly conserved. Nonetheless, WAK1 was found to be highly divergent when compared with sequences available from barley cv. Haruna Nijo (50 %), rice (46 %), wheat (21 %), Arabidopsis (25 %), and maize (19 %). This divergence may have facilitated a better adaptation to surrounding environments due to its role in communication between the extracellular matrix, cell, and outer environment. Semiquantitative RT-PCR-based expression analysis indicates HvWAK1 expression is specific to roots. Significant differences in root growth between GP wild type and GP-Ds mutant seedlings were observed under control and salt stress conditions.     

Codon usage in Xanthophyllomyces dendrorhous (formerly Phaffia rhodozyma)

By sequence analysis of 96 randomly selected clones in a cDNA library of Xanthophyllomyces dendrorhous, ten novel, full-length clones encoding cytoplasmic ribosomal proteins (rp) were found. The deduced amino acid sequences showed significant homology to their counterparts from eukaryotic origin including mammals, fungi and plants. Some ribosmal protein encoding cDNAs appeared several times, but by Southern blot analysis it was shown they are encoded by a single copy gene. The nucleotide sequences of ten full length cDNAs were used to investigate the codon usage in X. dendrorhous.     

Characterization of an HSP70 Cognate Gene Family in Arabidopsis

Analysis of the polypeptide composition of extracts from heat-shocked leaves of Arabidopsis indicated the presence of at least 12 HSP70-related polypeptides, most of which were constitutively expressed. In vitro translation of mRNA from heat-shocked and control leaves indicated that the amount of mRNA encoding four HSP70 polypeptides was increased strongly by heat-shock. Three Arabidopsis genes which exhibit homology to a Drosophila HSP70 gene were cloned. Two of the three genes are arranged in direct orientation approximately 1.5 kilobases apart. The third gene is not closely linked to the other two. Nucleotide sequence analysis of the 5′ regions of the two linked genes revealed that both contain a TATA box, the CAAT motif, and several short sequences which are homologous to the Drosophila heat-shock consensus sequence. The deduced partial amino acid sequence of the open reading frames were 79 and 72% homologous to the corresponding regions of the Drosophila HSP70-cognate and HSP70 sequences, respectively. As with the two maize HSP70 genes which have been characterized, and the Drosophila HSP70-cognate genes, the Arabidopsis genes contained a putative intron in the codon specifying amino acid 72. Analysis of mRNA levels with gene-specific oligonucleotide probes indicated that two of the genes were not expressed or were expressed at very low levels in leaves during normal growth or after heat-shock, whereas the other gene was constitutively expressed. By analogy with the results of similar studies of other organisms, it appears that the three cloned genes are members of a small family which are most closely related to the HSP70-cognate genes found in other species.     

Cloning and partial characteristics of the barley D-hordein sequence]

A number of clones containing major endosperm-specifically transcribed gene copies were selected from a cDNA library developed on the basis of barley endosperm mRNA. Approx. 30% of the recombinant clones carried sequences homologous to mRNA of various cereal storage proteins. Some of them appeared to be related to cDNA clones of wheat and barley storage proteins. The typical insert length ranged from 0.3 to 1.7 kB. A couple of clones among them were selected which revealed positive hybridization with all probes used. The positive signals disappeared after stringent washing of the filters. The nucleotide sequences of two representatives of the group were determined and corresponding amino acid sequence deduced after subsequent computer analysis. The comparison with known cereal storage protein genes revealed relatively high homology level with the central part of wheat high molecular weight (HMW) glutenine subunit genes. The fact suggests the cloned gene to belong to barley D-hordein family.