First whole genome based microsatellite DNA marker database of tomato for mapping and variety identification
© Iquebal et al.; licensee BioMed Central Ltd. 2013
Received: 4 March 2013
Accepted: 26 November 2013
Published: 4 December 2013
The cultivated tomato is second most consumed vegetable of the world and is an important part of a diverse and balanced diet as a rich source of vitamins, minerals, phenolic antioxidants and antioxidant lycopene having anti-cancer properties. To reap benefit of genomics of the domestic tomato (Solanum lycopersicum L.) unravelled by Tomato Genome Consortium (The Tomato Genome Consortium, 2012), the bulk mining of its markers in totality is imperative and critically required. The solgenomics has limited number of microsatellite DNA markers (2867) pertaining to solanaceae family. As these markers are of linkage map having relative distance, the choice of selected markers based on absolute distance as of physical map is missing. Only limited microsatellite markers with limitations are reported for variety identification thus there is a need for more markers supplementing DUS test and also for traceability of product in global market.
We present here the first whole genome based microsatellite DNA marker database of tomato, TomSatDB (Tomato MicroSatellite Database) with more than 1.4 million markers mined in-silico, using MIcroSAtellite (MISA) tool. To cater the customized needs of wet lab, features with a novelty of an automated primer designing tool is added. TomSatDB (http://cabindb.iasri.res.in/tomsatdb), a user-friendly and freely accessible tool offers chromosome wise as well as location wise search of primers. It is an online relational database based on “three-tier architecture” that catalogues information of microsatellites in MySQL and user-friendly interface developed using PHP (Hypertext Pre Processor).
Besides abiotic stress, tomato is known to have biotic stress due to its susceptibility over 200 diseases caused by pathogenic fungi, bacteria, viruses and nematodes. These markers are expected to pave the way of germplasm management over abiotic and biotic stress as well as improvement through molecular breeding, leading to increased tomato productivity in India as well as other parts of the world. In era of IPR the new variety can be identified based on allelic variation among varieties supplementing DUS test and product traceability.
Tomato (Solanum lycopersicon L.), a new world solanaceous plant is an excellent model for plant genomic research. The genus Solanum is one of the largest angiosperm genera and the genome has 35,000 genes spread over 12 chromosomes and has few high copy number long terminal repeat (LTR) retrotransposons and largely comprised of low-copy DNA . The genome of tomato has been sequenced by The Tomato Genome Consortium in 2012 . A draft sequence of its closest wild relative, i.e. Solanum pimpinellifolium and the potato genome (Solanum tuberosum L.) has also been reported depicting the extent and pattern of similarities and dissimilarities among the three genomes. The tomato genome sequence will have implications on other plant species viz., strawberries, melons, apple etc., which share some characteristics with tomato. Especially common information related to gene and pathway involved in fruit ripening can be potentially applied to other crops also leading to improved fruit quality and reduced cold chain management costs .
The leading tomato producing countries of the world are China, USA, India, Turkey, Egypt, and Italy . Globally, USA (23.2%), Germany (16.5%), Russian Federation (9.6%) and United Kingdom (8.3%) are top countries importing tomatoes while Netherlands is the biggest exporter of tomatoes, exporting over 910 346 tons a year and accounting for 20.6% of world export market in tomatoes .
Worldwide, tomatoes are an important part of a diverse and balanced diet as a rich source of vitamins, minerals, phenolic antioxidants and anti-oxidant lycopene having anti-cancer properties [5, 6]. A major constraint in tomato production is the loss incurred due to several diseases. Tomato is known to be susceptible to over 200 diseases caused by pathogenic fungi, bacteria, viruses and nematodes . To accelerate conventional plant breeding, Marker Assisted Selection (MAS) with more than 80 resistance genes to major classes of pathogens (fungal, bacterial, virus and nematode) have been used extensively for pyramiding resistance genes [8, 9]. The other essential characteristics for tomato improvement are development of cultivars with broad adaptability, earliness and fruit quality. Over 75,000 accessions of the cultivated and wild species of tomato are maintained in Genebank around the world  but relative and absolute differentiation of these accessions need microsatellite (STR) markers. Thus, these bulk whole genome based STR markers are needed for mapping, variety identification and product traceability.
Earlier in silico works for STR mining were not based on whole genome and thus invariably yielded low/very less number of markers, for example just 80 STR were found in entire gene bank search over 2000 sequences. Though the markers were minimum but are highly potential in distinguishing closely related cultivars of tomato . Reported marker density based on in vivo method is relatively less on every chromosome, for example 12th chromosome has just 37 . For plant variety identification along with degree of admixture, STR is always preferred if they are in multiplex mode (for example, Basmati and non-Basmati rice can be differentiated by 8 plex/ single cocktail based PCR) . Such multiplex designing needs much more number of markers to design multiplex with thermodynamic compatibility, which can be accomplished from our large marker dataset. SGN database is having various classes of markers including STR which are chromosome wise and distance wise based on LOD (Logarithm of the Odds) score . Uniformly distributed markers over genome with an average spacing of 10.0 cM  and 1.2 cM (ca. 900 KB)  are reported but further higher marker density with average spacing of less than 10 KB has not been reported so far.
Our present work aims at development of such first microsatellite marker database based on whole genome based STR mining which is very user-friendly and freely accessible. Also, the feature of user defined primer designing has great advantage in terms of precise selection from each chromosome, from defined location, size of amplicons for ease of rapid genotyping in simple and low cost agarose gel.
Construction and content
Database processing pipeline
The chromosome wise tomato whole genome data available in public domain  was downloaded in FASTA format. All the 12 available chromosomes of the genome were chopped into manageable range using PERL script to be put into MIcroSAtellite identification (MISA) tool . The information on STR numbers, motifs, repeat number, length and size of the repeat, repeat type, GC content, start and end position of the repeat and STR sequence were compiled. Around 1.4 million STRs were generated from tomato genome. Scripts in PERL were written to arrange the output of MISA in proper format in order to create the data file to be further imported to MYSQL database.
Option to find STRs from tomato genome specifically based on chromosome location, type of motif, repeat motif and repeat kind are available. The advance option for search is also available for STRs in desired range of GC content, number of base pairs and copy number. Further, selected STRs can be used in wet lab by generating primers with the integrated Primer3 standalone tool . This obviates the need of manual primer designing using tool/ server.
TomSatDB has seven tabs (Home, About, Database, Analysis, Tutorial, Links and Team) where general information of the developed tomato microsatellite database, tomato, microsatellite markers, analysis of the tomato genome have been described. Tutorial contains the guidelines for users and terminologies used in the database contents. TomSatDB houses other useful links related to solanaceae family.
The TomSatDB is very flexible and easy to handle where the user may query for microsatellites over 12 chromosomes, either one or more chromosomes being selected at a time from tomato genome. The searches may further be customized based on various microsatellite characteristics like motif type (mono, di, tri, tetra, penta, hexa), repeat motif and repeat kind (simple and composite).
Primer3 tool  has been integrated in TomSatDB to generate primers of selected STRs. The STRs may be selected with the help of radio button for generation of primers. Provision to design primer for selected STR locus is provided with a template of approximately 1000 base pairs by selecting upto 500 base pairs of both flanking regions. The provided flexibilities would enable researchers to select markers at known location over the desired chromosomes.
Further, each individual STR of a targeted region over chromosome may be used to narrow down location of gene of interest or linked QTL. The users are given flexibility to replace degenerate bases with any of the alternative bases (A,T,G,C) in TomSatDB.
Analysis of tomato genome and relative abundance
Chromosome wise distribution of STRs
STR validation result of primers from Solgenomics
STR Markers fromhttp://solgenomics.net/
Total no. of primers reported
No. of positive primers
Utility and discussion
Limited attempt of STR development was reported for example, 12th chromosome has just 37 markers by earlier in vivo method, but in silico method has extracted 11245 from very same chromosome giving much more marker density, which is highly desirable . Though the genome size range of solanaceae is varying from 0.950 GB (tomato) to 2.70 GB (capsicum) but this family has fixed 12 chromosomes. Although the gene repertoire and gene order of solanaceae species are well conserved, the cause of the genome-size difference is not known .
STR plays important role in mapping, trait improvement, variety development, variety identification and product traceability. Traditionally, characterization of varieties is based on phenotypic observation but it is very difficult to distinguish varieties with very similar morphological characteristics and identification of the cultivars accurately is essential for maintaining cultivar integrity and Plant Breeders’ Rights.
Limited studies have been reported in variety identification of tomato using STR DNA markers. In one study, out of 20 STR markers, only 11 were able to discriminate 47 varieties  and in another study, 12 markers could differentiate 34 varieties . Studies based on 6000 SNP markers over 93 varieties have demonstrated that SNP based variety differentiation is also possible . However, in such SNP based studies, the genotyping data of "Moneymaker" and "Moneyberg" varieties were completely identical leading to no differentiation at all. So, more STR markers from tomato genome are warranted to address varietal differentiation and product tractability in the food chain. Also, DNA fingerprinting is an appropriate tool to track and trace the tomato supply chain, ensuring not only authenticity and integrity of the products but also the absence of any possible genetic contamination by other species or unwanted components [30–32].
Such use of STR in plant variety identification is well reported in other crops like barley varieties , S. tuberosum ssp. tuberosum, sugarcane , capsicum , eggplant  and identification of Basmati rice from that of non-Basmati rice . Also, the microsatellite STR markers are the method of first choice to complement the DUS (Distinctness, Uniformity and Stability) testing procedure [38, 39].
Our STR database can be a useful tool in MAS programme of tomato improvement. Such use of STR in crop improvement is already reported in sorghum , tagging stem rust resistance gene Sr35 in wheat , Fusarium head blight resistance in wheat , leaf rust resistance gene Lr35 in wheat  and mapping of resistance gene effective against Karnal bunt pathogen of wheat . Wheat improvement programs to enhance leaf rust resistance using STR markers has been attempted . STR markers are also used for introgression programs for trait improvement, for example Soltol QTLs in rice. The location of the Saltol QTL on chromosome 1 and identification of additional QTLs associated with salt tolerance is well identified .
The relative density of the tomato whole genome reported in the study is 154 markers per MBp. This is almost in range with Arabidopsis (157 MBp), the other crops reported are with higher number of markers like cucumber (367 MBp), rice (370–490 MBp), popular (485 MBp), grape (487 MBp), sorghum (818 MBp), soybean (1115 MBp), maize (2365 MBp), wheat (1000 MBp) and pigeon pea (833 MBp). Though a general negative correlation between genome size and STR density in plants has been reported  but we found distance between markers are not proportion to size of genome thus small genome size has enough marker density for mapping purpose.
TomSatDB is of great use to tomato breeders in molecular breeding. The customization of this tool for search based on chromosome may be used by breeders for mapping of gene by markers. It is likely to be accessed by biologists engaged in research with diverse objectives in the crop primarily to develop molecular markers and also to understand the functional significance of microsatellites in regulating gene expression and genome evolution. The comprehensive options to search for simple and compound microsatellites repeats in the genic regions allow users to explore new avenues of investigations on these repeats. The primer designing for PCR amplification of desired motifs will facilitate studies on mutability, microsatellite abundance etc. Association of microsatellites with a particular disease or phenotype may also be explored. Microsatellite data can also be used to investigate various anomalies using candidate gene approach. This microsatellite database will serve as an important application for extracting information in order to design experiments in new directions elucidating novel roles and functions of microsatellites. The STR markers (>1.4 m) reported here is not only relevant for tomato germplasm management using MAS against 200 biotic and abiotic stress but also to other crops. This database is expected to be of immense use across globe by respective statutory authorities for variety identification and varietal dispute resolution supplementing DUS test and product traceability.
A total of 146602 STR markers are reported for the first time using whole genome in the database. Though we have reported small attempt of in silico validation in our studies but extensive wet lab validation of these markers is warranted. These markers are expected to pave the way of germplasm management over abiotic and biotic stress as well as improvement through molecular breeding, leading to increased tomato productivity in various parts of the world. The marker reported in our database are ready to use for mapping as well as also for variety identification and product traceability, paving the pathway of best use of genomics and computational tool in endeavor of tomato improvement and variety management at global level.
Availability and requirement
TomSatDB, the tomato microsatellite marker database is freely accessible for research purposes for non-profit and academic organizations at http://cabindb.iasri.res.in/tomsatdb.
Financial assistance received by National Agricultural Innovation Project, Indian Council of Agricultural Research, New Delhi entitled “Establishment of National Agriculture Bioinformatics Grid in ICAR” is gratefully acknowledged. The technical assistance of Jai Bhagwan in maintaining the web server and AR Paul in designing the logo of TomSatDB are thankfully acknowledged. Authors acknowledge the critical input of all the anonymous reviewers and editor in improvement of the manuscript.
- The Tomato Genome Consortium: The tomato genome sequence provides insights into fleshy fruit evolution. Nature Lett. 2012, 485: 635-641. 10.1038/nature11119.View ArticleGoogle Scholar
- Ranjan A, Ichihashi Y, Sinha NR: The tomato genome: implications for plant breeding, genomics and evolution. Genome Biol. 2012, 13: 16. 10.1186/gb-2012-13-3-r16.View ArticleGoogle Scholar
- Food and agriculture organization. http://faostat.fao.org.
- Directorate Marketing: A profile of the South African tomato market value chain. Department - Agriculture, forestry and fisheries: Republic of South Africa; 2012.Google Scholar
- Rick CM: Tomato. Hybridization of crop plants. Edited by: Fehr WR, Hadley HH. Madison, USA: American Society of Agronomy-Crop Science Society of America;1980, 669-680.Google Scholar
- Willcox JK, Catignani GL, Lazarus S: Tomatoes and cardiovascular health. Crit Rev Food Sci Nutri. 2003, 43 (1): 1-18. 10.1080/10408690390826437.View ArticleGoogle Scholar
- Lukyanenko AN: Disease resistance in tomato. Genetic Improvement of Tomato. Edited by: Berlin KG. Germany: Springer; 1991, 99-119. Monographs on Theor Appl Genet, vol 14View ArticleGoogle Scholar
- Grube RC, Radwanski E, Jahn M: Comparative genetics of disease resistance within the Solanaceae. Genetics. 2000, 155: 873-88.PubMedPubMed CentralGoogle Scholar
- Foolad MR, Panthee DR: Marker-assisted selection in tomato breeding. Crit Rev Plant Sci. 2012, 31 (2): 93-123. 10.1080/07352689.2011.616057.View ArticleGoogle Scholar
- Larry R, Joanne L: Genetic resources of tomato (Lycopersicon esculentum Mill.) and Wild Relatives. Genetic improvement of solanaceous crops, Volume 2-Tomato. Edited by: Razdan MK, Mattoo AK. New Hampshire, USA: Science Publishers; 2007, 25-75.Google Scholar
- Smulders MJM, Bredemeijer G, Rus-Kortekaas W, Arens P, Vosman B: Use of short microsatellites from database sequences to generate polymorphisms among Lycopersicon esculentum cultivars and accessions of other Lycopersicon species. Theor Appl Genet. 1997, 94 (2): 264-272. 10.1007/s001220050409.View ArticleGoogle Scholar
- Subramaniam G, Palchamy P, Robert P, Eguru R, Jaw-Fen W: Development of tomato SSR markers from anchored BAC clones of chromosome 12 and their application for genetic diversity analysis and linkage mapping. Euphytica. 2011, 178 (2): 283-295. 10.1007/s10681-010-0331-8.View ArticleGoogle Scholar
- Archak S, Lakshminarayanareddy V, Nagaraju J: High-throughput multiplex microsatellite marker assay for detection and quantification of adulteration in Basmati rice (Oryza sativa). Electrophoresis. 2007, 28: 2396-2405. 10.1002/elps.200600646.PubMedView ArticleGoogle Scholar
- Sol genomics network. http://solgenomics.net/.
- Frary A, Xu Y, Liu J, Mitchell S, Tedeschi E, Tanksley S: Development of a set of PCR-based anchor markers encompassing the tomato genome and evaluation of their usefulness for genetics and breeding experiments. Theor Appl Genet. 2005, 111: 291-312. 10.1007/s00122-005-2023-7.PubMedView ArticleGoogle Scholar
- Tanksley SD, Ganal MW, Prince JP, de Vicente MC, Bonierbale MW, Broun P, Fulton TM, Giovannoni JJ, Grandillo S, Martin GB, Messeguer R, Miller J, Miller L, Paterson AH, Pineda O, Riider MS, Wing RA, Wu W, Young ND: High density molecular linkage maps of the tomato and potato genomes. Genetics. 1992, 132: 1141-1160.PubMedPubMed CentralGoogle Scholar
- Solanum lycopersicum project. http://mips.helmholtz-muenchen.de/plant/tomato/index.jsp.
- MicroSatellite analysis tool. http://pgrc.ipk-gatersleben.de/misa/.
- Rozen S, Skaletsky HJ: Primer3 on the WWW for general users and for biologist programmers. Bioinformatics Methods and Protocols: Methods in Molecular Biology. Edited by: Krawetz S, Misener S. 2000, Totowa, NJ: Humana Press, 365-386.Google Scholar
- Sharma PC, Grover A, Kahl G: Mining microsatellites in eukaryotic genomes. Trends Biotechnol. 2007, 25 (11): 490-498. 10.1016/j.tibtech.2007.07.013.PubMedView ArticleGoogle Scholar
- Kariin S, Burge C: Dinucleotide relative abundance extremes: a genomic signature. Trends Genet. 1995, 11 (7): 283-290. 10.1016/S0168-9525(00)89076-9.View ArticleGoogle Scholar
- Shioiri C, Takahata N: Skew of mononucleotide frequencies, relative abundance of dinucleotides and DNA strand asymmetry. J Mol Evol. 2001, 53: 364-376. 10.1007/s002390010226.PubMedView ArticleGoogle Scholar
- Haseneyer G, Schmutzer T, Seidel M: From RNA-seq to large-scale genotyping-genomics resources for rye (Secale cereale L.). BMC Plant Biol. 2011, 11: 131-143. 10.1186/1471-2229-11-131.PubMedPubMed CentralView ArticleGoogle Scholar
- Zietkiewicz E, Rafalski A, Labuda D: Genome finger printing by simple sequence repeat (SSR) - anchored polymerase chain reaction amplification. Genomics. 1994, 20: 176-183. 10.1006/geno.1994.1151.PubMedView ArticleGoogle Scholar
- Kim KS, Min MS, An JH, Lee H: Cross-species amplification of bovidae microsatellites and Low diversity of the endangered Korean goral. J Hered. 2004, 95 (6): 521-525. 10.1093/jhered/esh082.PubMedView ArticleGoogle Scholar
- Park M, Jo SH, Kwon J, Park J, Ahn JH, Kim S, Lee YH, Yang TJ, Hur CG, Kang BC, Kim BD, Choi D: Comparative analysis of pepper and tomato reveals euchromatin expansion of pepper genome caused by differential accumulation of Ty3/Gypsy-like elements. BMC Genomics. 2011, 12: 85-97. 10.1186/1471-2164-12-85.PubMedPubMed CentralView ArticleGoogle Scholar
- Sardaro MLS, Marmiroli M, Maestri E, Marmiroli N: Genetic characterization of Italian tomato varieties and their traceability in tomato food products. Food Sci Nutri. 2013, 1 (1): 54-62. 10.1002/fsn3.8.View ArticleGoogle Scholar
- Srivastava DK, Shoemaker LS, Franks CE, Sussman MD: Single laboratory validation of a microsatellite marker-based method in tomato variety identification. J AOAC Int. 2011, 94 (1): 251-258.PubMedGoogle Scholar
- Viquez-Zamora M, Vosman B, van de Geest H, Bovy H, Visser RGF, Finkers R, van Heusden AW: Tomato breeding in the genomics era: insights from a SNP array. BMC Genomics. 2013, 14: 354. 10.1186/1471-2164-14-354.PubMedPubMed CentralView ArticleGoogle Scholar
- Marmiroli N, Peano C, Maestri E: Advanced PCR techniques in identifying food components. Food authenticity and traceability. Edited by: Lees M. Cambridge UK: Woodhead Publishing; 2003, 3-33.View ArticleGoogle Scholar
- Marmiroli N, Maestri E, Pafundo S, Vietina M: Molecular traceability of olive oil: From plant genomics to food genomics. Advances in olive resources. Edited by: Berti L, Maury J. Trivandrum, India: Transworld Research Network; 2009, 157-172.Google Scholar
- Agrimonti C, Vietina M, Pafundo S, Marmiroli N: The use of food genomics to ensure the traceability of olive oil. Trends Food Sci Tech. 2011, 22: 237-244. 10.1016/j.tifs.2011.02.002.View ArticleGoogle Scholar
- Karakousis A, Barr AR, Chalmers KJ, Ablett GA, Holton TA, Henry RJ, Lim P, Langridge P: Potential of SSR markers for plant breeding and variety identification in Australian barley germplasm. Aust J of Agr Res. 2003, 54: 1197-1210. 10.1071/AR02178.View ArticleGoogle Scholar
- Kawchuk LM, Martin RF, Mcpherson J: Resistance in transgenic potato expressing the potato leafroll virus coat protein gene. Mol Plant Microbe In. 1990, 3: 301-307. 10.1094/MPMI-3-301.View ArticleGoogle Scholar
- Manigbas NL, Villegas LC: Microsatellite Markers in Hybridity tests to identify true hybrids of sugarcane. Philipp J Crop Sci. 2004, 29 (2): 23-32.Google Scholar
- Shirasawa K, Ishii K, Kim C, Ban T, Suzuki M, Ito T, Muranaka T, Kobayashi M, Nagata N, Isobe S, Tabata S: Development of Capsicum EST-SSR markers for species identification and in silico mapping onto the tomato genome sequence. Mol Breeding. 2013, 31 (1): 101-110. 10.1007/s11032-012-9774-z.View ArticleGoogle Scholar
- Stagel A, Portis E, Toppino L, Rotino GL, Lanteri S: Gene-based microsatellite development for mapping and phylogeny studies in eggplant. BMC Genomics. 2008, 9: 357-370. 10.1186/1471-2164-9-357.PubMedPubMed CentralView ArticleGoogle Scholar
- McCouch SR, Chen X, Panaud O, Temnykh S, Xu Y, Cho YG, Huang N, Ishii T, Blair M: Microsatellite marker development, mapping and applications in rice genetics and breeding. Plant Mol Biol. 1997, 35: 89-99. 10.1023/A:1005711431474.PubMedView ArticleGoogle Scholar
- Becher SA, Steinmetz K, Weising K, Boury S, Peltier D, Renou JP, Kahl G, Wolff K: Microsatellites for variety identification in Pelargonium. Theor Appl Genet. 2000, 101: 643-651. 10.1007/s001220051526.View ArticleGoogle Scholar
- Wang Y, Bible P, Loganantharaj R, Upadhyaya HD: Identification of SSR markers associated with height using pool-based genome-wide association mapping in sorghum. Mol Breeding. 2012, 30 (1): 281-292. 10.1007/s11032-011-9617-3.View ArticleGoogle Scholar
- Babiker E, Ibrahim AMH, Yen Y, Stein J: Identification of a microsatellite marker associated with stem rust resistance gene Sr35 in wheat. Aust J Crop Sci. 2009, 3: 195-200.Google Scholar
- Liu S, Anderson JA: Marker assisted evaluation of fusarium head blight resistant wheat germplasm. Crop Sci. 2003, 43: 760-766. 10.2135/cropsci2003.0760.View ArticleGoogle Scholar
- Seyfarth R, Feuillet C, Schachermayr G, Winzeler M, Keller B: Development of a molecular marker for the adult plant leaf rust resistance gene Lr35 in wheat. Theor Appl Genet. 1999, 99 (3–4): 554-560.PubMedView ArticleGoogle Scholar
- Singh R, Datta D, Priyamvada , Singh S, Tiwari R: Marker-assisted selection for leaf rust resistance genes Lr19 and Lr24 in wheat (Triticum aestivum L.). Appl Genet. 2004, 45 (4): 399-403.Google Scholar
- Kolmer JA, Anderson JA, Flor JM: Chromosome location, linkage with simple sequence repeat markers, and leaf rust resistance conditioned by gene in wheat. Crop Sci. 2010, 50: 2392-2395. 10.2135/cropsci2010.01.0005.View ArticleGoogle Scholar
- Thomson MJ, de Ocampo M, Egdane J, Rahman MA, Sajise AG, Adorada DL, Tumimbang-Raiz E, Blumwald E, Seraj ZI, Singh RK, Gregorio GB, Ismail AM: Characterizing the Saltol quantitative trait locus for salinity tolerance in rice. Rice. 2010, 3: 148-160. 10.1007/s12284-010-9053-8.View ArticleGoogle Scholar
- Morgante M, Olivieri AM: PCR-amplified microsatellites as markers in plant genetics. Plant J. 1993, 3: 175-182. 10.1111/j.1365-313X.1993.tb00020.x.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.