Acceleration of Soybean Genomics Using Large Collections of DNA Markers for Gene Discovery

Legumes are known to provide nutritious proteins and vegetableoils while at the same time providing industrial products suchas biodiesel. It is estimated that approximately 25% of worldcrop production is derived from legumes. Recently, knowledgeof the molecular biology and genomics of legumes have been extendedsignificantly using two model species, Lotus japonicus (http://www.kazusa.or.jp/lotus/)and Medicago     

Sequence Analysis of the Genome of an Oil-Bearing Tree, Jatropha curcas L.

The whole genome of Jatropha curcas was sequenced, using a combination of the conventional Sanger method and new-generation multiplex sequencing methods. Total length of the non-redundant sequences thus obtained was 285 858 490 bp consisting of 120 586 contigs and 29 831 singlets. They accounted for ∼95% of the gene-containing regions with the average G + C content was 34.3%. A total of 40 929 complete and partial structures of protein encoding genes have been deduced. Comparison with genes of other plant species indicated that 1529 (4%) of the putative protein-encoding genes are specific to the Euphorbiaceae family. A high degree of microsynteny was observed with the genome of castor bean and, to a lesser extent, with those of soybean and Arabidopsis thaliana. In parallel with genome sequencing, cDNAs derived from leaf and callus tissues were subjected to pyrosequencing, and a total of 21 225 unigene data have been generated. Polymorphism analysis using microsatellite markers developed from the genomic sequence data obtained was performed with 12 J. curcas lines collected from various parts of the world to estimate their genetic diversity. The genomic sequence and accompanying information presented here are expected to serve as valuable resources for the acceleration of fundamental and applied research with J. curcas, especially in the fields of environment-related research such as biofuel production. Further information on the genomic sequences and DNA markers is available at http://www.kazusa.or.jp/jatropha/.     

Comparative Genomic Analysis of Two-Component Signal Transduction Systems in Probiotic Lactobacillus casei

Lactobacillus casei has traditionally been recognized as a probiotic, thus needing to survive the industrial production processes and transit through the gastrointestinal tract before providing benefit to human health. The two-component signal transduction system (TCS) plays important roles in sensing and reacting to environmental changes, which consists of a histidine kinase (HK) and a response regulator (RR). In this study we identified HKs and RRs of six sequenced L. casei strains. Ortholog analysis revealed 15 TCS clusters (HK–RR pairs), one orphan HKs and three orphan RRs, of which 12 TCS clusters were common to all six strains, three were absent in one strain. Further classification of the predicted HKs and RRs revealed interesting aspects of their putative functions. Some TCS clusters are involved with the response under the stress of the bile salts, acid, or oxidative, which contribute to survive the difficult journey through the human gastrointestinal tract. Computational predictions of 15 TCSs were verified by PCR experiments. This genomic level study of TCSs should provide valuable insights into the conservation and divergence of TCS proteins in the L. casei strains.     

SNP Discovery and Linkage Map Construction in Cultivated Tomato

Few intraspecific genetic linkage maps have been reported for cultivated tomato, mainly because genetic diversity within Solanum lycopersicum is much less than that between tomato species. Single nucleotide polymorphisms (SNPs), the most abundant source of genomic variation, are the most promising source of polymorphisms for the construction of linkage maps for closely related intraspecific lines. In this study, we developed SNP markers based on expressed sequence tags for the construction of intraspecific linkage maps in tomato. Out of the 5607 SNP positions detected through in silico analysis, 1536 were selected for high-throughput genotyping of two mapping populations derived from crosses between ‘Micro-Tom’ and either ‘Ailsa Craig’ or ‘M82’. A total of 1137 markers, including 793 out of the 1338 successfully genotyped SNPs, along with 344 simple sequence repeat and intronic polymorphism markers, were mapped onto two linkage maps, which covered 1467.8 and 1422.7 cM, respectively. The SNP markers developed were then screened against cultivated tomato lines in order to estimate the transferability of these SNPs to other breeding materials. The molecular markers and linkage maps represent a milestone in the genomics and genetics, and are the first step toward molecular breeding of cultivated tomato. Information on the DNA markers, linkage maps, and SNP genotypes for these tomato lines is available at http://www.kazusa.or.jp/tomato/.     

Sequence Analysis of the Genome of Carnation (Dianthus caryophyllus L.)

The whole-genome sequence of carnation (Dianthus caryophyllus L.) cv. ‘Francesco’ was determined using a combination of different new-generation multiplex sequencing platforms. The total length of the non-redundant sequences was 568 887 315 bp, consisting of 45 088 scaffolds, which covered 91% of the 622 Mb carnation genome estimated by k-mer analysis. The N50 values of contigs and scaffolds were 16 644 bp and 60 737 bp, respectively, and the longest scaffold was 1 287 144 bp. The average GC content of the contig sequences was 36%. A total of 1050, 13, 92 and 143 genes for tRNAs, rRNAs, snoRNA and miRNA, respectively, were identified in the assembled genomic sequences. For protein-encoding genes, 43 266 complete and partial gene structures excluding those in transposable elements were deduced. Gene coverage was ∼98%, as deduced from the coverage of the core eukaryotic genes. Intensive characterization of the assigned carnation genes and comparison with those of other plant species revealed characteristic features of the carnation genome. The results of this study will serve as a valuable resource for fundamental and applied research of carnation, especially for breeding new carnation varieties. Further information on the genomic sequences is available at http://carnation.kazusa.or.jp.     

Editor's choice: Genome-Wide SNP Genotyping to Infer the Effects on Gene Functions in Tomato

The genotype data of 7054 single nucleotide polymorphism (SNP) loci in 40 tomato lines, including inbred lines, F₁ hybrids, and wild relatives, were collected using Illumina''s Infinium and GoldenGate assay platforms, the latter of which was utilized in our previous study. The dendrogram based on the genotype data corresponded well to the breeding types of tomato and wild relatives. The SNPs were classified into six categories according to their positions in the genes predicted on the tomato genome sequence. The genes with SNPs were annotated by homology searches against the nucleotide and protein databases, as well as by domain searches, and they were classified into the functional categories defined by the NCBI''s eukaryotic orthologous groups (KOG). To infer the SNPs'' effects on the gene functions, the three-dimensional structures of the 843 proteins that were encoded by the genes with SNPs causing missense mutations were constructed by homology modelling, and 200 of these proteins were considered to carry non-synonymous amino acid substitutions in the predicted functional sites. The SNP information obtained in this study is available at the Kazusa Tomato Genomics Database (http://plant1.kazusa.or.jp/tomato/).     

Generation of 7137 non-redundant expressed sequence tags from a legume, Lotus japonicus.

For comprehensive analysis of genes expressed in a model legume, Lotus japonicus, a total of 22,983 5' end expressed sequence tags (ESTs) were accumulated from normalized and size-selected cDNA libraries constructed from young (2 weeks old) plants. The EST sequences were clustered into 7137 non-redundant groups. Similarity search against public non-redundant protein database indicated that 3302 groups showed similarity to genes of known function, 1143 groups to hypothetical genes, and 2692 were novel sequences. Homologues of 5 nodule-specific genes which have been reported in other legume species were contained in the collected ESTs, suggesting that the EST source generated in this study will become a useful tool for identification of genes related to legume-specific biological processes. The sequence data of individual ESTs are available at the web site: http://www.kazusa.or.jp/en/plant/lotus/EST/.     

Jatropha curcas hemagglutinin is similar to a 2S albumin allergen from the same source and has unique sugar affinities

We have previously reported the purification and preliminary X-ray characterization of a hemagglutinin from the seeds of Jatropha curcas and, with the detailed sequencing information available now, we find that it is similar to a 2S albumin allergen isolated from the same source. Through a search of Jatropha genome database (http://www.kazusa.or.jp/jatropha/), we map it to the sequence id JcCA0234191 (now referred to as Jcr4S00619.70 in the new version, release 4.5) which has a conserved alpha amylase inhibitor/seed storage protein domain found in the 2S albumin allergens. The putative sequence of the small and large chains of the protein is assigned and the total mass of the two subunits matches with the intact mass 10?kDa determined through MALDI. The protein retains hemagglutination activity between pH 6–9 and up to 60?°C on heat treatment and its hemagglutination activity is inhibited by sialic acid and fetuin. Bioinformatics studies show that the isolated protein sequence clusters in close association with a 2S albumin from Ricinus communis in phylogeny analysis and has a conservation of the characteristic four disulfide linkage pattern. Hemagglutinins and lectins are known to have allergenic effects through their interaction with immunoglobulin E and histamine release and earlier studies have shown that this interaction can be inhibited by lectin-specific sugars. We hope this report bridges the plant allergens and hemagglutinins further for exploring possible mediation of allergenic activity through sialic acid and complex sugar interactions and generates further interest in the area.     

An interspecific linkage map of SSR and intronic polymorphism markers in tomato

Despite the collection and availability of abundant tomato genome sequences, PCR-based markers adapted to large scale analysis have not been developed in tomato species. Therefore, using public genome sequence data in tomato, we developed three types of DNA markers: expressed sequence tag (EST)-derived simple sequence repeat (SSR) markers (TES markers), genome-derived SSR markers (TGS markers) and EST-derived intronic polymorphism markers (TEI markers). A total of 2,047 TES, 3,510 TGS and 674 TEI markers were established and used in the polymorphic analysis of a cultivated tomato (Solanum lycopersicum) ‘LA925’ and its wild relative Solanum pennellii ‘LA716’, parents of the Tomato-EXPEN 2000 mapping population. The polymorphic ratios between parents revealed by the TES, TGS and TEI markers were 37.3, 22.6 and 80.0%, respectively. Those showing polymorphisms were used to genotype the Tomato-EXPEN 2000 mapping population, and a high-density genetic linkage map composed of 1,433 new and 683 existing marker loci was constructed on 12 chromosomes, covering 1,503.1 cM. In the present map, 48% of the mapped TGS loci were located within heterochromatic regions, while 18 and 21% of TES and TEI loci, respectively, were located in heterochromatin. The large number of SSR and SNP markers developed in this study provide easily handling genomic tools for molecular breeding in tomato. Information on the DNA markers developed in this study is available at http://www.kazusa.or.jp/tomato/.     

Evolutionary history of mitogen-activated protein kinase (MAPK) genes in Lotus,Medicago, and Phaseolus

Mitogen-Activated Protein Kinase (MAPK) genes encode proteins that mediate various signaling pathways associated with biotic and abiotic stress responses in eukaryotes. The MAPK genes form a 3-tier signal transduction cascade between cellular stimuli and physiological responses. Recent identification of soybean MAPKs and availability of genome sequences from other legume species allowed us to identify their MAPK genes. The main objectives of this study were to identify MAPKs in 3 legume species, Lotus japonicus, Medicago truncatula, and Phaseolus vulgaris, and to assess their phylogenetic relationships. We used approaches in comparative genomics for MAPK gene identification and named the newly identified genes following Arabidopsis MAPK nomenclature model. We identified 19, 18, and 15 MAPKs and 7, 4, and 9 MAPKKs in the genome of Lotus japonicus, Medicago truncatula, and Phaseolus vulgaris, respectively. Within clade placement of MAPKs and MAPKKs in the 3 legume species were consistent with those in soybean and Arabidopsis. Among 5 clades of MAPKs, 4 founder clades were consistent to MAPKs of other plant species and orthologs of MAPK genes in the fifth clade-"Clade E" were consistent with those in soybean. Our results also indicated that some gene duplication events might have occurred prior to eudicot-monocot divergence. Highly diversified MAPKs in soybean relative to those in 3 other legume species are attributable to the polyploidization events in soybean. The identification of the MAPK genes in the legume species is important for the legume crop improvement; and evolutionary relationships and functional divergence of these gene members provide insights into plant genome evolution.     

Structural analysis of a Lotus japonicus genome. V. Sequence features and mapping of sixty-four TAC clones which cover the 6.4 mb regions of the genome.

We determined the nucleotide sequences of 64 TAC (transformation-competent artificial chromosome) clones selected from genomic libraries of Lotus japonicus accession Miyakojima MG-20 based on the sequence information of expressed sequence tags (ESTs), cDNAs, genes and DNA markers from L. japonicus and other legumes. The length of the DNA regions sequenced in this study was 6,370,255 bp, and the total length of the L. japonicus genome sequenced so far is 32,537,698 bp together with the nucleotide sequences of 256 TAC clones previously reported. Five hundred forty-eight potential protein-encoding genes with known or predicted functions, 127 gene segments and 224 pseudogenes were assigned to the newly sequenced regions by computer prediction and similarity searches against the sequences in protein and EST databases. Based on the nucleotide sequences of the clones, simple sequence repeat length polymorphism (SSLP) or derived cleaved amplified polymorphic sequence (dCAPS) markers were generated, and each clone was genetically localized onto the linkage map of two accessions of L. japonicus, MG-20 and Gifu B-129. The sequence data, gene information and mapping information are available through the World Wide Web at http://www.kazusa.or.jp/lotus/.     

Structural analysis of a Lotus japonicus genome. I. Sequence features and mapping of fifty-six TAC clones which cover the 5.4 mb regions of the genome.

A total of 56 TAC clones with an average insert size of 100 kb were isolated from a TAC library of the Lotus japonicus genome based on the expressed sequences tags (ESTs), cDNA and gene information, and their nucleotide sequences were determined according to the shot-gun based strategy. The total length of the sequenced regions is 5,473,195 bp. By comparison with the sequences in protein and EST databases and analysis with computer programs for gene modeling, a total of 605 potential protein-encoding genes with known or predicted functions, 69 gene segments, and 172 pseudogenes were identified. The average density of the genes assigned so far is 1 gene/8120 bp. Introns were identified in approximately 78% of the potential genes. There was an average of 3.8 introns per gene and the average length of the introns was 375 bp. DNA markers were generated based on the nucleotide sequences obtained, and each clone was mapped onto the linkage map using the F2 mapping population derived from a cross of L. japonicus Gifu B-129 and Miyakojima MG-20. The sequence data, gene information and mapping information are available through the World Wide Web at http://www.kazusa.or.jp/lotus/.     

Structural analysis of a Lotus japonicus genome. IV. Sequence features and mapping of seventy-three TAC clones which cover the 7.5 mb regions of the genome.

Using the sequence information of expressed sequences tags (ESTs), cDNAs and genes from Lotus japonicus and other legumes, 73 TAC (transformation-competent artificial chromosomes) clones were selected from a genomic library of L. japonicus accession MG-20, and their nucleotide sequences were determined. The length of the DNA sequenced in this study was 7,455,959 bp, and the total length of the DNA regions sequenced so far is 26,167,443 bp together with the nucleotide sequences of 183 TAC clones previously reported. By similarity searches against the sequences in protein and EST databases and prediction by computer programs, a total of 699 potential protein-encoding genes with known or predicted functions, 163 gene segments and 267 pseudogenes were assigned to the newly sequenced regions. Based oil the nucleotide sequences of the clones, simple sequence repeat length polymorphism (SSLP) or derived cleaved amplified polymorphic sequence (dCAPS) markers were generated, and each clone was located onto the linkage map of two accessions of L. japonicus, Gifu B-129 and Miyakojima MG-20. The sequence data, gene information and mapping information are available through the World Wide Web at http://www.kazusa.or.jp/lotus/.     

Genome-wide identification and expression analysis of two component system genes in Cicer arietinum

The two-component system (TCS) plays an important role in signal transduction pathways, cytokinin signaling and stress resistance of prokaryotes and eukaryotes. It is comprised of three types of proteins in plants; histidine kinases (HKs), histidine phosphotransfer proteins (HPs) and response regulators (RRs). Chickpea (Cicer arietinum L.) is one of the most important legume crops worldwide with special economic value in semi-arid tropics. Availability of complete genome sequence of chickpea presents a valuable resource for comparative analysis among angiosperms. In current study, Arabidopsis thaliana and Oryza sativa were used as reference plant species for comparative genomics analysis with C. arietinum. A genome-wide computational survey enabled us to identify putative members of TCS protein family including 18HKs, 26 RRs (7 type-A, 7 type-B, 2 type C and 10 pseudo) and 7 HPs (5 true and 2pseudo) genes in chickpea. The predicted TCS genes displayed family specific intron/exon organization and were randomly distributed across all the eight chromosomes. Comparative phylogenetic and evolutionary analysis suggested a variable conservation of TCS genes in relation to mono/dicot model plants and segmental duplication was the principal route of expansion for this family in chickpea. The promoter regions of TCS genes exhibited several abiotic stress-related cis-elements indicating their involvement in abiotic stress response. The expression analysis of TCS genes demonstrated stress (drought, heat, osmotic and salt) specific differential expression. Current study provides insight into TCS genes in C. arietinum, which will be helpful for further functional analysis of these genes in response to different abiotic stresses.     

Codon usage tabulated from international DNA sequence databases: status for the year 2000

The frequencies of each of the 257 468 complete protein coding sequences (CDSs) have been compiled from the taxonomical divisions of the GenBank DNA sequence database. The sum of the codons used by 8792 organisms has also been calculated. The data files can be obtained from the anonymous ftp sites of DDBJ, Kazusa and EBI. A list of the codon usage of genes and the sum of the codons used by each organism can be obtained through the web site http://www.kazusa.or.jp/codon/ . The present study also reports recent developments on the WWW site. The new web interface provides data in the CodonFrequency-compatible format as well as in the traditional table format. The use of the database is facilitated by keyword based search analysis and the availability of codon usage tables for selected genes from each species. These new tools will provide users with the ability to further analyze for variations in codon usage among different genomes.     

A Genome-Wide RNAi Screen for Enhancers of par Mutants Reveals New Contributors to Early Embryonic Polarity in Caenorhabditis elegans

The par genes of Caenorhabditis elegans are essential for establishment and maintenance of early embryo polarity and their homologs in other organisms are crucial polarity regulators in diverse cell types. Forward genetic screens and simple RNAi depletion screens have identified additional conserved regulators of polarity in C. elegans; genes with redundant functions, however, will be missed by these approaches. To identify such genes, we have performed a genome-wide RNAi screen for enhancers of lethality in conditional par-1 and par-4 mutants. We have identified 18 genes for which depletion is synthetically lethal with par-1 or par-4, or both, but produces little embryo lethality in wild type. Fifteen of the 18 genes identified in our screen are not previously known to function in C. elegans embryo polarity and 11 of them also increase lethality in a par-2 mutant. Among the strongest synthetic lethal genes, polarity defects are more apparent in par-2 early embryos than in par-1 or par-4, except for strd-1(RNAi), which enhances early polarity phenotypes in all three mutants. One strong enhancer of par-1 and par-2 lethality, F25B5.2, corresponds to nop-1, a regulator of actomyosin contractility for which the molecular identity was previously unknown. Other putative polarity enhancers identified in our screen encode cytoskeletal and membrane proteins, kinases, chaperones, and sumoylation and deubiquitylation proteins. Further studies of these genes should give mechanistic insight into pathways regulating establishment and maintenance of cell polarity.     

Structural analysis of a Lotus japonicus genome. III. Sequence features and mapping of sixty-two TAC clones which cover the 6.7 Mb regions of the genome.

A total of sixty-two clones were selected from a TAC (transformation-competent artificial chromosome) genomic library of the Lotus japonicus accession MG-20 based on the sequence information of expressed sequence tags (ESTs), cDNA and gene information, and their nucleotide sequences were determined. The length of the sequenced regions in this study is 6,682,189 bp, and the total length of the regions sequenced so far is 18,711,484 bp together with the nucleotide sequences of 121 TAC clones previously reported. By comparison with the sequences in protein and EST databases and analysis with computer programs for gene modeling, a total of 573 potential protein-coding genes with known or predicted functions, 91 gene segments and 272 pseudogenes were identified in the newly sequenced regions. Each of the sequenced clones was localized onto the linkage map of two accessions of L. japonicus, Gifu B-129 and Miyakojima MG-20, using simple sequence repeat length polymorphism (SSLP) or derived cleaved amplified polymorphic sequence (dCAPS) markers generated based on the nucleotide sequences of the clones. The sequence data, gene information and mapping information are available through the World Wide Web at http://www.kazusa.or.jp/lotus/.     

Genomic database resources for Dictyostelium discoideum

Dictyostelium is an attractive model system for the study of mechanisms basic to cellular function or complex multicellular developmental processes. Recent advances in Dictyostelium genomics have generated a wide spectrum of resources. However, much of the current genomic sequence information is still not currently available through GenBank or related databases. Thus, many investigators are unaware that extensive sequence data from Dictyostelium has been compiled, or of its availability and access. Here, we discuss progress in Dictyostelium genomics and gene annotation, and highlight the primary portals for sequence access, manipulation and analysis (http://genome.imb-jena.de/dictyostelium/; http://dictygenome.bcm.tmc.edu/; http://www.sanger. ac.uk/Projects/D_discoideum/; http://www.csm.biol. tsukuba.ac.jp/cDNAproject.html).