Tomato Genomic Resources Database: An Integrated Repository of Useful Tomato Genomic Information for Basic and Applied Research

PLOS ONE, Jan 2014

Tomato Genomic Resources Database (TGRD) allows interactive browsing of tomato genes, micro RNAs, simple sequence repeats (SSRs), important quantitative trait loci and Tomato-EXPEN 2000 genetic map altogether or separately along twelve chromosomes of tomato in a single window. The database is created using sequence of the cultivar Heinz 1706. High quality single nucleotide polymorphic (SNP) sites between the genes of Heinz 1706 and the wild tomato S. pimpinellifolium LA1589 are also included. Genes are classified into different families. 5′-upstream sequences (5′-US) of all the genes and their tissue-specific expression profiles are provided. Sequences of the microRNA loci and their putative target genes are catalogued. Genes and 5′-US show presence of SSRs and SNPs. SSRs located in the genomic, genic and 5′-US can be analysed separately for the presence of any particular motif. Primer sequences for all the SSRs and flanking sequences for all the genic SNPs have been provided. TGRD is a user-friendly web-accessible relational database and uses CMAP viewer for graphical scanning of all the features. Integration and graphical presentation of important genomic information will facilitate better and easier use of tomato genome. TGRD can be accessed as an open source repository at

A PDF file should load here. If you do not see its contents the file may be temporarily unavailable at the journal website or you do not have a PDF plug-in installed and enabled in your browser.

Alternatively, you can download the file locally and open with any standalone PDF reader:

Tomato Genomic Resources Database: An Integrated Repository of Useful Tomato Genomic Information for Basic and Applied Research

Chattopadhyay D (2014) Tomato Genomic Resources Database: An Integrated Repository of Useful Tomato Genomic Information for Basic and Applied Research. PLoS ONE 9(1): e86387. doi:10.1371/journal.pone.0086387 Tomato Genomic Resources Database: An Integrated Repository of Useful Tomato Genomic Information for Basic and Applied Research B. Venkata Suresh 0 Riti Roy 0 Kamlesh Sahu 0 Gopal Misra 0 Debasis Chattopadhyay 0 Tongming Yin, Nanjing Forestry University, China 0 National Institute of Plant Genome Research, Aruna Asaf Ali Marg , New Delhi , India Tomato Genomic Resources Database (TGRD) allows interactive browsing of tomato genes, micro RNAs, simple sequence repeats (SSRs), important quantitative trait loci and Tomato-EXPEN 2000 genetic map altogether or separately along twelve chromosomes of tomato in a single window. The database is created using sequence of the cultivar Heinz 1706. High quality single nucleotide polymorphic (SNP) sites between the genes of Heinz 1706 and the wild tomato S. pimpinellifolium LA1589 are also included. Genes are classified into different families. 59-upstream sequences (59-US) of all the genes and their tissuespecific expression profiles are provided. Sequences of the microRNA loci and their putative target genes are catalogued. Genes and 59-US show presence of SSRs and SNPs. SSRs located in the genomic, genic and 59-US can be analysed separately for the presence of any particular motif. Primer sequences for all the SSRs and flanking sequences for all the genic SNPs have been provided. TGRD is a user-friendly web-accessible relational database and uses CMAP viewer for graphical scanning of all the features. Integration and graphical presentation of important genomic information will facilitate better and easier use of tomato genome. TGRD can be accessed as an open source repository at - Funding: This work was financially supported by the core grant of National Institute of Plant Genome Research, New Delhi, India and the Department of Biotechnology, Government of India. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Competing Interests: The authors have declared that no competing interests exist. Tomato (Solanum lycopersicum) is model for a number of biological studies important to agriculture such as, fruit development and ripening, disease resistance and biochemical pathways of important nutrients [13]. Tomato, a Solanaceae member, has been used for pioneering research on cell wall and storage polysaccharide synthesis and degradation [46]. Tomato was the first plant from which a gene for gene class of R-gene for disease resistance was cloned [7]. Tomato fruits are a rich source of carotenoid pigment, which is the precursor of vitamin A. Tomato has been used extensively for genetic studies because of several reasons such as, its diploid genome, short generation time, availability of homozygous inbred lines, efficient transformation technology and its genes are largely sequestered in contiguous euchromatic regions [810]. Apart from basic research, tomato is the second mostconsumed vegetable in the world. Researchers from fourteen different countries together have published a reference genome sequence of an inbred tomato cultivar Heinz1706 and a draft sequence of wild tomato, Solanum pimpinellifolium [11]. These sequence information have provided the basic and applied researchers an opportunity to scout for gene function, genetic diversity and evolution not only in tomato, but also in other Solanaceae members, for studying basic biology and bio-diversity based breeding. The Sol genomics network (SGN) serves as the most referred host for storing and integrating most of the information on Solanaceae crops and is continuously developing [12]. There has been existence of a number of other tomato databases on specific aspects even before the release of tomato genome such as, SolEST database, MoTo DB, Tomato functional genomics database, miSolRNA, TOMATOMA, KaTomicsDB ( and several others [1318]. These databases serve extensively as resources for genomic and biochemical information and biological material for the Solanaceae research community. However, there is a need to amalgamate different genomic and trait information together graphically along all the chromosomes of tomato in a single window to enable the basic and applied researchers to integrate and utilize all the information in a better and easier way. Microsatellite markers (simple sequence repeats, SSRs) are important for various applications such as construction of highdensity linkage maps, comparative genome mapping, identification of variety, marker-assisted selection, studying genetic diversity and so on. They are reliable because of high experimental reproducibility, multiallelic nature and co-dominant inheritance [19]. Simple sequence repeats (SSR) are the most popular microsatellite markers. Genes are the most important functional part of a genome and define the trait of an organism. Integration of SSRs with genes in a single window will allow the researchers to look for presence of SSRs in the coding, non-coding and the upstream activating sequences of a gene and subsequently for intra- and inter-species polymorphism. Genetic variation of a trait within a species can be attributed to a single gene or joint action of many genes that can be mapped on the genome (QTL, Quantitative trait loci) by genetic markers [20]. Integration of wellknown QTLs with the genes and SSRs along chromosomal length will provide the users an easier access to the genes and SSRs colocalizing with a QTL. Changes in gene expression may lead to phenotypic differences between two individuals. Apart from transcriptional regulation, gene expression can be regulated at various stages including through microRNAs (miRNAs). Several studies demonstrated that development and metabolism of plants are regulated by miRNAs [2123]. miRNAs might also have potential role in transgressive phenotype [24]. Therefore, inclusion of tissue-specific expression of genes, miRNAs and their predicted target genes together with SSRs, genes and QTLs would be useful. Small and red-fruited Solanum pimpinellifolium is the closest wild relative of domesticated tomato. Introgression lines for S. pimpinellifolium in the background of cultivated tomato are used for mapping different traits. There is a 0.6% nucleotide divergence between the inbred cultivar Heinz1706 and the S. pimpinellifolium accession LA1589. This huge resource of SNPs, especially in the genic region, would immensely benefit the breeders to utilize the natural trait reservoir for crop improvement. In view of above, we have integrated SSRs, genes, miRNAs, known QTLs and SNPs between the sequenced cultivated and wild tomato accessions in a single window along the length of all the twelve chromosomes of tomato. Materials and Methods Sequence Retrieval S. lycopersicum Heinz 1706 reference genome sequence, raw sequence reads of S. pimpinellifolium LA1589, latest annotation (ITAG 2.3), and gff files were downloaded from Sol genomics network database (SGN) ( The RNAseq data was downloaded from short read archive of NCBI [accession number SRX118613 (leaf) SRX118614 (root) SRX118615 (flower) SRX118621 (mature green fruit)]. QTL and EXPEN2000 map coordinates were retrieved from SGN. Microsatellites/SSRs MIcroSAtellite identification tool, MISA ( was used to scan the tomato reference genome with default parameters to identify repeats. MISA program allows the user to specify the minimal length of the consecutive nucleotide repeat and reports the SSR type, SSR motif, motif repeat, length of repeat and coordinates of the SSRs in the genome. Mononucleotide repeats were not included. The minimum number of repeats allowed for the dimers was six and for tri-hexamers was five. Primer3 program was used to design the primers from the flanking regions within 100 bp of each of the identified SSR [25]. Coordinates for all the genes were retrieved from the gff file of ITAG2.3 annotation. Coordinates and sequence of 2 kilobase (kb)-long upstream sequences (59-US) of the genes were extracted with PERL script. Accordingly, SSRs were assigned to 59-US and different parts (exon, intron, and 39 untranslated region) of the genes in addition to other parts of the genome. All the SSRs, for which primers could be designed, were anchored on to all the twelve chromosomes and can be visualized interactively using CMAP viewer [26]. SNPs and miRNAs Raw Illumina reads of S pimpinellifolium were filtered for high quality through NGS QC Toolkit v2.3 with default parameter [27]. Paired end reads were mapped on tomato reference genome through Bowtie version 2.1.0 with default parameters [28]. SAMTools was used to convert SAM file to BAM file, remove duplicate reads and for SNP calling [29]. Stringent SNP filtering criteria of minimum read depth 5, minimum root-mean-square mapping quality 30 and all the mapped reads showing same nonallelic base were applied. SNPs in the protein-coding genes were extracted using gene coordinates and accordingly assigned to the different parts of the genes as mentioned above. MicroRNAs with their targets were retrieved based on the available literature and from miRBase database, and the co-ordinates were retrieved from.gff file of reference genome [30,31,11]. Genetic Markers and QTLs Genetic markers and QTLs were retrieved from EXPEN-2000 and QTL genetic maps, respectively, from SGN database and literature (Ty1/3, Pto, Bs4, Cf-9, Ve, Mi) [3237]. All the EXPEN2000 genetic markers were anchored to the chromosomes according to their physical locations. All the QTLs were assigned according to the physical locations of the nearest flanking sequence-characterized genetic markers mentioned in SGN. Where the sequences of the flanking markers were not available, physical coordinates of the corresponding EXPEN2000 genetic positions, if available, were retrieved from Kazusa Tomato Marker Database ( for assigning to chromosomal locations. Gene Expression RNA sequence data for four tissues, namely leaf (3 week-old), root (3 week-old), flower (unopened bud) and mature green fruit, of tomato cv. Heinz1706 was retrieved as mentioned above and used for the analysis. Illumina RNAseq reads were first aligned to ribosomal RNA sequences using Bowtie 2 to eliminate possible rRNA sequence contamination. To quantify the expressions of genes all the Illumina reads from 4 tissue samples were mapped onto the genes using RSEM [38]. The number of reads mapped was normalized by RPKM (reads per kilobase per million) method. The heat map showing gene specific expression was generated on the RPKM for each gene in all the tissue samples using TIGR MultiExperiment Viewer [39,40]. Differential gene expression analysis was performed using DESEQ [41]. The genes showing greater than two fold expression with p#0.05 were regarded as differentially expressed. Genes expressing with more than 5 RPKM in a tissue as opposed to less than 1 RPKM in other tissues were considered as preferentially expressed. Genes preferentially expressed in each tissue sample, as compared with others, in a tissue-by-tissue comparison is presented as heat map. To find out collinear gene blocks on tomato genome, all-versus-all BlastP match was performed using the tomato proteome (version: ITAG2.3) with cut off e-value of 1e25. The BLASTP result was fed into MCScanX toolkit ( mcscan2/) to identify collinear blocks with the following parameters; e-value of 1e25 and match size of 10. A total of 773 collinear blocks were detected with 20 genes as a minimum block size. Structure of TGRD The Tomato Genomic Resources Database (http://59.163.192. 91/tomato2/) is an online and interactive relational database developed using open source software, Apache (http://www., MySQL 5.0 (, JavaScript and PHP 5.4 ( and is hosted on a 64-bit Linux Server. The database was designed based on Three-Level Schema Architecture (figure S1). A flow chart explaining the database architecture, organization and workflow of TGRD has been presented in figure 1. The CMAP viewer and its schema were integrated with TGRD for graphical scanning and analysis of various genomic features. The user-friendly interface for TGRD has been developed using PHP 5.4, JavaScript and HTML to query and retrieve the data based on user needs. The interface is completely interactive and interlinked with each component of the genome. For example, the gene page contains information about the associated SSRs and SNPs; and the genic SSR page contains information about the gene, which contains the SSR and so on. The TGRD database contains a robust graphical tutorial to facilitate better use of the database. Results and Discussions Sequence Based Physical Mapping All the genomic features, eg. genes, SSRs, genic SNPs, QTL, miRNAs and genetic markers (EXPEN-2000) are mapped on twelve chromosomes and visualized by using CMAP viewer. The CMAP viewer can be accessed from MAPs tab. The MAP page is developed with rich interactivity using JavaScript and HTML. It shows all the twelve chromosomes in circular form with interactivity. The user can visualize all the twelve chromosomes together (All Chromosomes tab at the centre) or individually with all or individual genomic features to have better comprehension about tomato genome. On CMAP viewer, each genomic feature contains links to get further information. The CMAP viewer contains Feature, Display and Advanced options at the bottom to control and reconstruct the map to view the features on the physical map as per users option. For better utilization of space, the overlapping features were collapsed. Choosing No option for Collapse Overlapping Features in Feature option would separate the features. Different features are color-coded. The maps can be cropped and magnified for any particular region of the chromosome to see the detail features. EXPEN-2000 genetic markers are included because they are well known to the breeders and, therefore, the other feature can be viewed with respect to the positions of these markers. EXPEN-2000 markers are linked to SGN site for detail information. Examples of two CMAP viewer pages with all the genomic features are presented in figure 2. Simple Sequence Repeats We mined the reference genome of tomato for microsatellite repeats. Total 68,641 microsatellite repeat motifs spanning 781 Mb assembled genome were identified. Dinucleotide repeats (60.18%) are much more abundant than tri (19.56%) and other repeats. Scaffolds assigned to twelve chromosomes spanning 760 Mb contain 66,823 sequence repeats of which, 55,396 (82.90%) and 11,427 (17.10%) were simple and compound repeats, respectively. Of the 55,396 simple repeats, primer pairs for 38,659 SSRs (69.78%) could be successfully designed. Frequency and distribution of various SSRs were analyzed. Chromosome-wise repeat statistics is presented in figure 3 and table S1. All the SSRs, for which primer pairs were designed, were anchored to individual chromosomes according to their coordinates. Physical locations of the SSRs in the exon, intron or 2 kb 59upstream regions of the genes were assigned. The user can query the database from SSRs tab using motif sequence, chromosome number, two coordinates of a chromosome or combination of repeat length and minimum number of repeat. The SSRs present in the gene coding (Genic SSRs) or in the 59-upstream sequences of the genes (59-US SSRs) also can be searched. Total of 5841 and 4773 SSRs were present in 33,840 chromosome-assigned genes and their 59-upstream sequences, respectively, with average frequencies of 0.172 SSRs/gene and 0.14 SSRs/59US. The query results information on marker ID, motif sequence, motif length, motif repeat, repeat length, region, chromosomal location, start and end positions, functional information in case of genic SSRs, primer information and link to CMAP to view location on chromosomes for each entry present in the database. These SSRs, physically localized along the chromosomes, would serve as an immense resource for comparative genome mapping, genetic diversity study and identification of polymorphic marker for quantitative trait loci. Genes, Gene Families and microRNAs According to ITAG2.3 annotation, tomato reference genome contains 34,727 protein-coding genes, of which 33,840 are assigned to twelve chromosomes so far. Chromosome 1 has the highest estimated (108.0 Mb) and assembled (90.30 Mb) lengths and accordingly, it codes for the highest number (4293) of genes. Whereas, the assembled chromosome 11 possesses the least (2385) number of genes, although its estimated (64.7 Mb) and assembled (53.4 Mb) lengths are more than the estimated (53.8 Mb) and assembled (46.04 Mb) lengths of the smallest tomato chromosome 6, which codes for 2813 genes [42,11]. Each gene record provides information about chromosomal location, gene expression (RPKM values), sequence information (gene, CDS, protein, 59US) and associated SSRs and SNPs at different locations of a gene. Under the GENE SEARCH tab, the user can search for any gene by ID, chromosome number, gene function and two coordinates of any chromosome. This tab is facilitated with BLAST search to look for any homologous tomato gene. Genes coding for transcription factors, heat shock proteins, protein kinases and transporter proteins are generally highly studied. Tomato has been extensively used as a model for disease resistance and fruit development. Apart from various transcription factors, enzymes and receptors for ethylene synthesis and perception, and red-light photoreceptors influence fruit development, ripening and quality. In addition, several enzymes related to modification of cell wall architecture have agronomic importance for fruit quality. Cytochrome P-450 family of genes involved in toxic alkaloid production showed significant contraction in tomato [11]. Therefore, genes encoding these protein families are specifically mentioned under the tab GENE FAMILIEs. All these genes can be viewed on CMAP viewer as a whole family or individually on all the twelve chromosomes. Total 2505 transcription factors assigned to the chromosomes are further subdivided into 89 families according to the domains present in them. Chromosome 1 possesses the highest number (287) of genes encoding transcription factors, while chromosome 9 contains the least (150). Similarly, 510 R genes are subdivided into eight families including the receptor like protein kinases. Chromosome 4 contains the highest number (85) of R genes. Functional genomics involves tissue-specific and treatment-specific expression of genes at different developmental stages of plants to predict function of a gene. We analyzed RNAseq data of tomato to reveal tissue-specific expression profiles of tomato genes. Hierarchical expression of all the genes in a family in leaf, root, flower and matured green fruit is provided in the form of heat maps and in RPKM values for individual genes. Total 6404 genes did not show any expression in these four tissues at these developmental stages. 73 genes were found to express preferentially in leaf while 665, 442 and 43 genes were found to express preferentially in root, flower and mature green fruit, respectively, at the specific developmental stage as mentioned (Fig. 4 and table S2). Physical positions of all the genes in a family can be viewed under this tab and, therefore, paralogous genes can be identified. Further, functional divergence of the genes present in collinear blocks can be analyzed from expression data as shown for two collinear blocks on chromosome 1 and 2 (Fig. 5). Integrated with QTL map, physical map of the genes provide an opportunity to associate candidate genes with important agronomic traits. 887 genes yet to be assigned to the chromosomes were listed with their detail under chromosome zero separately. Out of 96 conserved miRNAs annotated in genome sequence, 93 have been assigned to twelve chromosomes and can be searched under microRNAs tab. They are categorized in 26 families. As mentioned in the materials and methods, target genes for these miRNAs have been predicted based on literature and miRBase database. Each miRNA family contains information about their target genes, chromosomal location, strand orientation, mature and pre-mature RNA and DNA sequences. Single Nucleotide Polymorphisms (SNPs) Cultivated tomato, which experienced a genetic bottleneck while carried from South America to Europe, has narrow genetic base with estimated ,5% genetic variation in their wild relatives [43]. The wide genetic variation present in the self-compatible and incompatible wild tomato species is being investigated and exploited in tomato improvement by introgression breeding. S. pimpinellifolium is used to study tomato traits such as fruit size and shape and to introduce bacterial spot resistance. High quality genomic sequence reads of S. pimpinellifolium were mapped to the tomato reference genome. Stringent filtering criteria were applied to extract highly probable SNPs between these two accessions. Total 2,971,025 high quality SNPs were obtained of which 290,445 are in protein-coding genes. Chromosome 7 and 12 possesses the highest (5361/Mb) and lowest (2551/Mb) frequencies of SNPs, respectively (table S3). Frequency distribution of genomic and genic SNPs and genes in each million base interval were plotted for each chromosome (figures S2 and S3) for better understanding of SNP distribution across chromosomes and especially in euchromatic regions. Except chromosome 12, heterochromatic regions of all the chromosomes showed higher sequence divergence than the euchromatic part (figure S2). Comparison of distributions of genes and genic SNPs (figure S3) identified ten regions of 1 Mb interval with twenty or more SNPs per gene and all these regions are in the low-gene heterochromatin regions of the respective chromosomes (table S4). A list of genes containing more than ten SNPs is listed in table S5. Most of them encode unknown proteins. Three genes encoding calciumtransporting ATPase, two encoding WD-repeat proteins and one encoding receptor-like protein kinase are among the annotated high SNP-containing protein-coding genes. All the genic SNPs are anchored to their respective locations in each chromosome and can be accessed directly from map or from the genes. Flanking sequences of 100-base length for each SNP are provided for designing primers to convert them as markers. The genic SNPs identified can be utilized in finding allelic variation of known genes of economic importance and may initiate new studies. As an example, fruit size locus fw2.2 is one of the QTLs selected during tomato domestication and accounts for as much as 30% of the difference between the fruit sizes of smallfruited wild S. pennellii and large-fruited domesticated S. lycopersicum [44]. fw2.2 corresponds to a gene (Solyc02g090730.2.1, Chr 02:46830407-46831197) that encodes a protein Cell number regulator 19 having structural similarity to human oncogene Ras. Of 42 SNPs between two alleles, only three SNPs change three amino acids within first nine residues of the protein. Therefore, the allelic variation in fw2.2 was thought to modulate fruit size by differentially regulating carpel cell number due to differential expression of two genes rather than structural differences in proteins [45,46]. Sequence comparison of Solyc02g090730.2.1 of domesticated Heinz1706 and another small-fruited wild S. pimpinellifolium LA1589 detected only two SNPs in the open reading frame. Of them, only one (SNP-Chr-02-23437), common for both S. pennellii and S. pimpinellifolium, changes the third amino acid. However, the 59-upstream sequence (2 kb from translation start site), which regulates transcriptional expression of a gene, showed twenty one structural variations between Heinz1706 and LA1589, of which 15 are SNPs and six are in-dels (table S6). Significance of these allelic variations in regulating fruit size awaits experimental validation. Quantitative Trait Loci Genetic mapping of a QTL allows researchers to look for candidate genes segregating around a trait locus, which ultimately results in QTL characterization. However, use of different mapping populations for QTL mapping hinders unified physical localization of multiple QTLs on a single physical map. We tried to anchor some QTLs to tomato chromosomes according to the physical locations of their flanking markers. We used QTLs described in SGN QTL map. This map is based on a mapping population derived from a cherry tomato line Cervil and a round large-fruited line Levovil, which was different from the mapping population (S. lycopersicum M82 X L. pennellii LA716) used for EXPEN2000 map [20]. Hence, the boundaries of QTLs could not be delineated accurately due to unavailability of sequencecharacterized flanking genetic markers. The nearest sequenceFigure 5. Differential expression profiles of genes present in two collinear blocks located on chromosome 1 and 2. The color scale (0 20) represents RPKM values. IDs of the genes in the collinear blocks in chromosome 1 and 2 are mentioned in the left and right sides, respectively. doi:10.1371/journal.pone.0086387.g005 characterized genetic markers or the corresponding EXPEN2000 genetic locations were, therefore, used to physically localize the QTLs, which resulted in assignment of a little wider physical length than the actual. Following this approach we could localize fifty-four QTLs to eight chromosomes. According to the recombination nodule frequency map the assigned QTLs are in high recombination frequency regions of the corresponding chromosomes [11]. Diseases caused by different bacteria, fungus and viruses tremendously attenuate yield potential of tomato. Resistance loci associated with diseases caused by Tomato yellow leaf curl virus (Ty1/3), Pseudomonas syringae (Pto), Xanthomonas campestris (Bs4), Cladosporium fulvum (Cf-9), Verticillium albo-atrum (Ve) and root knot nematode (Mi) were identified by map-based cloning (3237). Physical locations of these pathogen-resistance loci, as derived from their marker or gene sequences, were also shown on the map as references. Tomato genome sequence has created a huge impact on Solanaceae research. The reference genome is increasingly being used for sequence-based approaches to answer basic biological questions of plant development as well as for agricultural improvement of fruit quality and quantity. Although several databases for tomato in addition to SGN exist, we tried to provide the most sought-after basic information about a genome in TGRD in a user-friendly way and we feel that graphical presentation of this important genomic information will facilitate a better use of tomato genome. A tutorial with graphical and video description has been provided for better use of the site. Regulation of gene expression is equally important as the gene sequence in influencing traits of an organism. For this reason we provided 2 kb-long 59upstream sequences with all the genes. Incorporation of SSR- and SNP-information within this region will be useful in exploiting intra- and inter-specific polymorphic sequences in studying functional divergence of homologous genes. Focus of the breeding programs carried out in the last century has been on increase in yield while fruit quality traits did not receive proper importance [47]. Improvement in traits such as, flavor, taste and nutrient content is challenging because these are regulated by biochemically complex processes. Recent advances in mapping some of the QTLs for these traits will lead to our understanding of these complex traits. Availability of tomato genome sequence and effort to assign physical locations of some of the flavor-related QTLs on the reference genome integrated with SSRs and SNPs in TGRD may facilitate biochemical characterization of these traits. Transcriptome sequences of a few wild relatives of tomato are now available [48]. We will keep on updating our site by incorporating SNP data of those wild accessions once their genome sequence data are available to allow robust mining of SNPs with stringent filtering criteria. Supporting Information Figure S1 Description of schema used to construct tomato genome database. Table S2 Genes preferentially/specifically expressed in each tissue sample, as compared to others in a tissueby-tissue comparison. Genes showing at least two-fold change (upregulated above the blank cell and downregulated below the blank cell) as compared to the other tissue samples were given. (DOC) Table S4 Genomic regions showing high frequency of SNP per gene in 1 Mb interval between S. lycopersicum Heinz1706 and S. pimpinellifolium LA1589. (XLS) Acknowledgement is due to the Director, National Institute of Plant Genome Research, New Delhi, India (NIPGR) for providing facilities. RR and KS acknowledge NIPGR and University grants commission for research fellowships. Computational resources under BTISNET program at NIPGR are gratefully acknowledged. Conceived and designed the experiments: DC. Performed the experiments: BVS RR KS GM. Analyzed the data: DC BVS RR KS GM. Wrote the paper: DC. 1. Giovannoni J ( 2001 ) Molecular biology of fruit maturation and ripening . Ann Rev Plant Physiol Plant Mol Biol 52 : 725 - 749 . 2. Adams-Phillips L , Barry C , Giovannoni J ( 2004 ) Signal transduction systems regulating fruit ripening . Trends Plant Sci 9 : 331 - 8 . 3. Ercolano MR , Sanseverino W , Carli P , Ferriello F , Frusciante L ( 2012 ) Genetic and genomic approaches for R-gene mediated disease resistance in tomato: retrospects and prospects . Plant Cell Rep 31 : 973 - 85 . 4. Greve LC , Labavitch JM ( 1991 ) Cell Wall Metabolism in Ripening Fruit: V. Analysis of Cell Wall Synthesis in Ripening Tomato Pericarp Tissue Using a d [U-C] Glucose Tracer and Gas Chromatography-Mass Spectrometry . Plant Physiol 97 : 1456 - 1461 . 5. Brady C , MacAlpine G , McGlassen WB ( 1982 ) Polygalacturonase in tomato fruits and the induction of ripening . Aust J Plant Physiol 9 : 171 - 178 . 6. Buckeridge MS ( 2010 ) Seed cell wall storage polysaccharides: models to understand cell wall biosynthesis and degradation . Plant Physiol 154 : 1017 - 1023 . 7. Martin GB , Brommonschenkel SH , Chunwongse J , Frary A , Ganal MW , et al. ( 1993 ) Map-based cloning of a protein kinase gene conferring disease resistance in tomato . Science 262 : 1432 - 1436 . 8. Khush GS , Rick CM , Robinson RW ( 1964 ) Genetic Activity in a Heterochromatic Chromosome Segment of the Tomato . Science 145 : 1432 - 1434 . 9. Emmanuel E , Levy AA ( 2002 ) Tomato mutants as tools for functional genomics . Curr Opin Plant Biol 5 : 112 - 117 . 10. Zamir D ( 2001 ) Improving plant breeding with exotic genetic libraries . Nat Rev Genet 2 : 983 - 989 . 11. Tomato Genome Consortium ( 2012 ) The tomato genome sequence provides insights into fleshy fruit evolution . Nature 485 : 635 - 641 . 12. Bombarely A , Menda N , Tecle IY , Buels RM , Strickler S , et al. ( 2011 ) The Sol Genomics Network ( growing tomatoes using Perl . Nucleic Acids Res 39 : D1149 - 55 . 13. D'Agostino N , Traini A , Frusciante L , Chiusano ML ( 2009 ) SolEST database: a ''one-stop shop'' approach to the study of Solanaceae transcriptomes . BMC Plant Biol 9 : 142 . 14. Grennan AK ( 2009 ) MoTo DB: a metabolic database for tomato . Plant Physiol 151 : 1701 - 1702 . 15. Fei Z , Joung J , Tang X , Zheng Y , Huang M , et al. ( 2011 ) Tomato Functional Genomics Database: a comprehensive resource and analysis package for tomato functional genomics . Nucleic Acids Research 39 : D1156 - D1163 . 16. Bazzini AA , Ass R , Gonzalez V , Bassi S , Conte M , et al. ( 2010 ) miSolRNA: A tomato micro RNA relational database . BMC Plant Biology 10 : 240 . 17. Saito T , Ariizumi T , Okabe Y , Asamizu E , Hiwasa-Tanase K , et al. ( 2011 ) TOMATOMA: a novel tomato mutant database distributing Micro-Tom mutant collections . Plant Cell Physiol 52 : 283 - 296 . 18. Shirasawa K , Hirakawa H ( 2013 ) DNA marker applications to molecular genetics and genomics in tomato . Breed Sci 63 : 21 - 30 . 19. Kalia RK , Rai MK , Kalia S , Singh R , Dhawan AK ( 2011 ) Microsatellites markers: an overview of the recent progress in Plants . Euphytica 177 : 309 - 334 . 20. Causse M , Saliba-Colombani V , Lecomte L , Duffe P, Rousselle P , et al. ( 2002 ), QTL analysis of fruit quality in fresh market tomato: a few chromosome regions control the variation of sensory and instrumental traits . J Exp Bot 53 : 2089 - 2098 . 21. Chen K , Rajewsky N ( 2007 ) The evolution of gene regulation by transcription factors and microRNAs . Nat Rev Genet 8 : 93 - 103 . 22. Sunkar R , Li YF , Jagadeeswaran G ( 2012 ) Functions of microRNAs in plant stress responses . Trends Plant Sci 17 : 196 - 203 . 23. Rubio-Somoza I , Weigel D ( 2011 ) MicroRNA networks and developmental plasticity in plants . Trends Plant Sci 16 : 258 - 264 . 24. Shivaprasad PV , Chen HM , Patel K , Bond DM , Santos BA , et al. ( 2012 ) A microRNA superfamily regulates nucleotide binding site-leucine-rich repeats and other mRNAs . Plant Cell 24 : 859 - 874 . 25. Untergasser A , Cutcutache I , Koressaar T , Ye J , Faircloth BC , et al. ( 2012 ) Primer3-new capabilities and interfaces . Nucleic Acids Res 40 : e115 . 26. Youens-Clark K , Faga B , Yap IV , Stein L , Ware D ( 2009 ) CMap 1.01: a comparative mapping application for the Internet . Bioinformatics 25 : 3040 - 3042 . 27. Patel RK , Jain M ( 2012 ) NGS QC Toolkit: a toolkit for quality control of next generation sequencing data . PLoS One 7 : e30619 . 28. Langmead B , Salzberg SL ( 2012 ) Fast gapped-read alignment with Bowtie 2 . Nat Methods 9 : 357 - 3599 . 29. Li H , Handsaker B , Wysoker A , Fennell T , Ruan J , et al. ( 2009 ) The Sequence alignment/map (SAM) format and SAMtools . Bioinformatics 25 : 2078 - 2079 . 30. Karlova R , van Haarst JC , Maliepaard C , van de Geest H , Bovy AG , et al. ( 2013 ) Identification of microRNA targets in tomato fruit development using high-throughput sequencing and degradome analysis . J Exp Bot 64 : 1863 - 1878 . 31. Kozomara A , Griffiths-Jones S ( 2011 ) miRBase: integrating microRNA annotation and deep-sequencing data . Nucleic Acids Res 39 : D152 - 157 . 32. Verlaan MG , Hutton SF , Ibrahem RM , Kormelink R , Visser RG , et al. ( 2013 ) The tomato yellow leaf curl virus resistance genes Ty-1 and Ty-3 are allelic and code for DFDGD-class RNA-dependent RNA polymerases . PLoS Genet 9 : e1003399 . 33. Martin GB , Williams JGK , Tanksley SD ( 1991 ) Rapid identification of markers linked to a Pseudomonas resistance gene in tomato by using random primers and near-isogenic lines . Proc Natl Acad Sci 88 : 2336 - 2340 . 34. Schornack S , Ballvora A , Gurlebeck D , Gurlebeck D , Peart J , et al. ( 2004 ) The tomato resistance protein Bs4 is a predicted non-nuclear TIR-NB-LRR protein that mediates defense responses to severely truncated derivatives of AvrBs4 and overexpressed AvrBs3 . Plant Jr 37 : 46 - 60 . 35. Jones DA , Thomas CM , Hammondkosack KE , Balintkurti PJ , Jones JDG ( 1994 ) Isolation of the tomato Cf-9 gene for resistance to Cladosporium fulvum by transposon tagging . Science 266 : 789 - 793 . 36. Kawchuk LM , Hachey J , Lynch DR , Kulcsar F , van Rooijen G , et al. ( 2001 ) Tomato Ve disease resistance genes encode cell surface-like receptors . Proc Natl Acad Sci USA 98 : 6511 - 6515 . 37. Milligan SB , Bodeau J , Yaghoobi J , Kaloshian I , Zabel P , et al. ( 1998 ) The root knot nematode resistance gene Mi from tomato is a member of the leucine zipper, nucleotide binding, leucine-rich repeat family of plant genes . Plant Cell 10 : 1307 - 1319 . 38. Bo L , Dewey CN ( 2011 ) RSEM: accurate transcript quantification from RNASeq data with or without a reference genome . BMC Bioinformatics 12 : 323 . 39. Saeed AI , Bhagabati NK , Braisted JC , Liang W , Sharov V , et al ( 2006 ) TM4 microarray software suite . Methods in Enzymology 411 : 134 - 193 . 40. Saeed AI , Sharov V , White J , Li J , Liang W , et al. ( 2003 ) TM4: a free, opensource system for microarray data management and analysis . Biotechniques 34 ( 2 ): 374 - 8 . 41. Anders S , Huber W ( 2010 ) Differential expression analysis for sequence count data . Genome Biol 11 : R106 . 42. Mueller LA , Lankhorst RK , Tanksley SD , Giovannoni JJ , White R , et al. ( 2009 ) A snapshot of the emerging tomato genome sequence . Plant Genome 2 : 78 - 92 . 43. Miller JC , Tanksley SD ( 1990 ) RFLP analysis of phylogenetic relationships and genetic variation in the genus Lycopersicon . Theoretical and Applied Genetics 80 : 437 - 448 . 44. Tanksley SD ( 2004 ) The Genetic, Developmental, and Molecular Bases of Fruit Size and Shape Variation in Tomato, The Plant Cell 16 : S181 - S189 . 45. Frary A , Nesbitt TC , Grandillo S , Knaap E , Cong B et al. ( 2000 ) fw2.2: a quantitative trait locus key to the evolution of tomato fruit size . Science 289 : 85 - 88 . 46. Cong B , Liu J , Tanksley SD ( 2002 ) Natural alleles at a tomato fruit size quantitative trait locus differ by heterochronic regulatory mutations , Proc Natl Acad Sci USA 99 : 13606 - 13611 . 47. Klee HJ , Tieman DM ( 2013 ) Genetic challenges of flavor improvement in tomato . Trends in Genetics 29 : 257 - 262 . 48. Koenig D , Jimenez-Gomez JM , Kimura S , Fulop D , Chitwood DH et al. ( 2013 ) Comparative transcriptomics reveals patterns of selection in domesticated and wild tomato . Proc Natl Acad Sci USA 110 : E2655 - 62 .

This is a preview of a remote PDF:

B. Venkata Suresh, Riti Roy, Kamlesh Sahu, Gopal Misra, Debasis Chattopadhyay. Tomato Genomic Resources Database: An Integrated Repository of Useful Tomato Genomic Information for Basic and Applied Research, PLOS ONE, 2014, DOI: 10.1371/journal.pone.0086387