- Open Access
solQTL: a tool for QTL analysis, visualization and linking to genomes at SGN database
BMC Bioinformaticsvolume 11, Article number: 525 (2010)
A common approach to understanding the genetic basis of complex traits is through identification of associated quantitative trait loci (QTL). Fine mapping QTLs requires several generations of backcrosses and analysis of large populations, which is time-consuming and costly effort. Furthermore, as entire genomes are being sequenced and an increasing amount of genetic and expression data are being generated, a challenge remains: linking phenotypic variation to the underlying genomic variation. To identify candidate genes and understand the molecular basis underlying the phenotypic variation of traits, bioinformatic approaches are needed to exploit information such as genetic map, expression and whole genome sequence data of organisms in biological databases.
The Sol Genomics Network (SGN, http://solgenomics.net) is a primary repository for phenotypic, genetic, genomic, expression and metabolic data for the Solanaceae family and other related Asterids species and houses a variety of bioinformatics tools. SGN has implemented a new approach to QTL data organization, storage, analysis, and cross-links with other relevant data in internal and external databases. The new QTL module, solQTL, http://solgenomics.net/qtl/, employs a user-friendly web interface for uploading raw phenotype and genotype data to the database, R/QTL mapping software for on-the-fly QTL analysis and algorithms for online visualization and cross-referencing of QTLs to relevant datasets and tools such as the SGN Comparative Map Viewer and Genome Browser. Here, we describe the development of the solQTL module and demonstrate its application.
solQTL allows Solanaceae researchers to upload raw genotype and phenotype data to SGN, perform QTL analysis and dynamically cross-link to relevant genetic, expression and genome annotations. Exploration and synthesis of the relevant data is expected to help facilitate identification of candidate genes underlying phenotypic variation and markers more closely linked to QTLs. solQTL is freely available on SGN and can be used in private or public mode.
Many economically important traits are quantitatively inherited, influenced by the environment and controlled by many genes of small and large effect [1–4]. A powerful approach to understanding the genetic basis of such complex traits is through identification of quantitative trait loci (QTL) associated with the trait [1, 5, 6]. Once QTLs are identified, generally, confirmation and fine mapping of the QTLs is required, followed by map-based cloning or a candidate gene approach that, in turn, necessitates further functional analysis of the candidate genes through expression profiling, genetic complementation and/or insertional mutagenesis [2, 3, 7, 8]. Fine mapping of QTL requires mendelization of the locus  and large populations. This effectively means several generations of backcrosses to create nearly isogenic lines (NILs), which is a costly and time consuming effort [10, 11]. Furthermore, dissecting the quantitative trait nucleotide (QTN) variation and understanding the molecular and physiological basis underlying the phenotypic variation remains a great challenge [3, 12]. A step forward in identifying candidate genes and dissecting the QTN(s) would be to develop tools to harness the genetic resources as well as expression and whole genome sequence data of organisms in biological databases through bioinformatics approaches [3, 11–13].
The Sol Genomics Network (SGN, http://solgenomics.net; ) is a repository for phenotypic, genetic, genomic, expression and metabolic data for the Solanaceae and other related Asterids families and houses a variety of bioinformatic tools. SGN has added a new QTL module, solQTL http://solgenomics.net/qtl/, to its repertoire of bioinformatics resources. The solQTL module presents a new approach to QTL data organization, storage, and analysis. Additionally, the module allows dynamic cross-linking with other relevant data including to the annotated genome sequence, expression data and genetic maps found in internal and external databases. The solQTL module employs a user-friendly web interface for uploading raw phenotype and genotype data, R/QTL mapping software (http://www.rqtl.org/), an add-on package to the R statistical software, for on-the-fly QTL analysis and algorithms to cross-reference QTLs to relevant datasets.
In primary mapping populations, QTL regions may span up to 30 cM  containing thousands of genes [3, 16]. However, some studies have revealed that QTLs and genes of quantitative traits can be identified from within 0 to 1.9 cM for major QTLs, which explain >25% of the phenotypic variation, (e.g. fw2.2, ; ovate,  in tomato) and 0 to 3 cM for small QTLs (<25%) (Hd1, Hd2, and Hd3 in rice ) from the peak marker . Price (2006), in his review article, implied that searches for candidate genes of small QTLs can be done within 1 - 2 cM on either side of the mean QTL position in a primary QTL population, especially when fine mapping is problematic but datasets from multiple studies can be used. However, it is known that genes or QTNs that affect the QTL may be tens of kilobases (kbs) away from the mean QTL position or in genes with unrelated or regulatory functions [1, 3, 12, 20]. Nonetheless, exploration and synthesis of genetic, genomic and expression data on corresponding regions around the peak QTL and within region delimited by flanking markers is useful, and can yield data leading to identification of candidate genes for the QTL of interest.
SGN, as a primary data and bioinformatics center for the Solanaceae, has the database structure, information and comparative analysis tools at hand to facilitate the cross-referencing of QTLs with other internal and external datasets. The SGN database houses the tomato genome sequence, Solanum lycopersicum cv Heinz 1706, including annotations provided by the International Tomato Genome Annotation Group (ITAG) . SGN also hosts whole genome sequences of S. tuberosum cv DM1-3 516R44 (currently BAC-based sequencing only) and S. pimpinellifolium cv LA1589. With the on-going efforts to sequence 100 more Solanaceae species (SOL-100 project, ) and the SNP database under development at SGN (unpublished), we envision that anchoring QTLs to the genomes will open additional windows in facilitating the elucidation of the molecular basis of complex traits and may lead in identification of QTL markers for efficient marker-assisted breeding.
Currently, major plant model organism databases such as Gramene (http://www.gramene.org[23, 24], MaizeGDB (http://www.maizegdb.org, and animal databases such as Animal QTLdb (http://www.animalgenome.org, [26, 27]) rely on in-house curators or the research community to curate QTL mapping outputs from published literature. A problem with QTL data is that historically, a large part of the original data has been lost, as the raw data (phenotypic and genotypic scores of the populations) were not recorded in the publications. It is therefore impossible to evaluate many of these described QTLs independently and retroactively. The solQTL module is unique in that (1) it uses a user-friendly, guided, step-wise, web interface to allow researchers to upload their raw phenotype and genotype data; (2) researchers can choose the statistical parameters for the on-the-fly, real-time, QTL mapping; (3) users can opt to make their raw and analyzed data either private or public; and (4) as growing number of journals are requiring authors to deposit their data in public databases prior to publication of their work, solQTL could facilitate this for the Solanaceae research community, if it was made a requirement for QTL data. Importantly, this system stores all the primary data for the identification of a QTL, such that the data can be re-evaluated in the future.
Here, we describe the development of the solQTL tool and demonstrate its application.
Construction and content
The researchers upload their raw phenotype and genotype data into the solQTL database and select the statistical parameters for the QTL mapping and visualization. A step-by-step, user-friendly, web interface (Additional file 1) allows users to upload their QTL population details such as name, description, and parents, trait data (description, definition and units of measurement), phenotypic and genotypic data separately in a tab-delimited files. To allow cross-referencing QTL markers to relevant datasets in SGN, markers used in the mapping population must exist in the SGN marker database. Users are prompted to set the statistical parameters they would like to apply for the QTL analysis. The statistical options available at this stage are described in the 'QTL Statistical Analysis' section below. Users have the option to save their data in the database, for either public or private use. The public option allows all users to see the dataset, whereas the private option only lets the data owner and SGN curators access the data. For private data, the population name, description, and contact information of the data owner are publicly accessible.
Built-in validation algorithms check for mandatory fields that need to be filled in and verify data formatting; for the statistical parameters, the consistency of the selection of the parameters is also validated. For example, the combination of QTL mapping method and QTL genotype probability estimation method selected is checked. Appropriate informative messages are displayed to the user when a required data field is empty or filled with incompatible information. To assist users in the use of the tool, a link to the guideline document is available at each step in the uploading and analysis process http://solgenomics.net/qtl/guide.pl.
Currently, the solQTL database holds raw QTL data for three F2 intercross tomato QTL populations [28–30] submitted by Esther van der Knaap of Ohio State University. For each population, 18 - 46 fruit morphology related traits are evaluated for QTLs.
QTL statistical analysis
The solQTL module uses R/QTL for the QTL mapping analysis. Among the QTL mapping functions presently implemented using the R/QTL in the solQTL tool include a one-dimension single-QTL scan model along with the single marker analysis and interval mapping methods. With the single-QTL model, the genome is scanned at a self-determined genome size, one at a time, for the presence of a QTL . Genome scan sizes of 0 (marker regression), 1, 2.5, 5 or 10 cM can be selected. Users can choose between QTL mapping methods such as standard interval mapping (Maximum Likelihood with EM algorithm [32, 33]), Haley-Knott regression  and Multiple Imputation . For the Multiple Imputation QTL mapping method and the QTL genotype probability estimation using Simulation, one can select either 5, 10, 15, or 20 imputations to run. The underlying algorithm for calculating the QTL genotype probabilities and for dealing with missing or partially missing genotype data is hidden Markov model (HMM) technology . A permutation test  for determining the LOD threshold significance level is also implemented. Users have the option to run 100 or 1000 permutation tests or none at all. QTL confidence interval (flanking markers) is determined based on a 95% Bayesian Credible Interval method with the interval expanded to the nearest markers [35, 38]. Based on the highest LOD score in the linkage group, a peak marker for a QTL is identified using R/QTL's 'find.marker' method. When a marker is not located at the peak QTL, the nearest possible marker is identified .
Currently, solQTL uses genetic map data from a priori computation and thus users need to provide linkage map data with the genotype data file, as indicated in documentation provided with the tool http://solgenomics.net/qtl/guide.pl. The linkage map plotting and visualization, however, is done using Perl modules developed both in-house and from the open sources community via the Comprehensive Perl Archive Network (CPAN, http://cpan.org).
At this stage, the solQTL is available only for backcross and F2 intercross diploid populations for traits with continuous phenotypic variation. As the primary focus of solQTL is to link the predicted QTLs to annotated whole genome sequences, dense reference genetic maps and other relevant expression data at SGN using tools such as the SGN Comparative Map Viewer  and genome browser (GBrowse), the tool is designed for use at SGN and hence, applicable and most useful for organisms that are supported at SGN.
SGN currently uses PostgreSQL http://www.postgresql.org version 8.3 for the back-end database. The website is implemented in Perl http://www.perl.org version 5.10.0. In-house and externally developed Perl modules are used for organizing, formatting and communicating data to R for the statistical analysis, the visualization of the QTL mapping output on the browser, and the cross-linking of QTLs to relevant datasets on SGN. A 1.08-56 version of the R/QTL (http://www.rqtl.org, ) is implemented on a 2.7.1 version of the R statistical software (http://www.r-project.org, [40, 41]).
QTL mapping and visualization
Once a user successfully uploads the prerequisite data in the supported format and sets the statistical parameters, the respective QTL population detail page is displayed (Additional file 2). For the purpose of demonstrating the QTL analysis, output visualization and cross-referencing to relevant genetic and genomic datasets with solQTL, we use raw data from the F2 population of Solanum lycopersicum cv. Howard German x S. pimpinellifolium accession LA1589 http://solgenomics.net/phenome/population.pl?population_id=12, and specifically the trait "Fruit Area" . The statistical parameters used for the particular QTL mapping output are in the corresponding legend. However, it should be noted that validation of the predicted QTL for the trait is beyond the scope of this work. On the population webpage, the name and description of the population is displayed along with a list of the traits phenotyped with the corresponding descriptive statistics. Also on the same page are any publications that document the QTL data, if available, and links to download the public raw phenotype and genotypic data of the QTL population.
The QTL analysis is performed dynamically on a trait-by-trait basis when the user visits the trait of interest on the QTL population's webpage. On the QTL population's page, under the QTL(s) heading, clicking the graph icon corresponding to the trait of interest initiates the QTL mapping analysis. The analysis output is displayed on a new trait page where one can view its QTL analysis output (Figure 1) and frequency distribution of the phenotype data for the trait. QTL mapping locations from the output of the QTL analysis are viewable for the whole genome on a per-linkage group basis. The details of the analysis are displayed in a legend including genome-wide LOD threshold significance level, if the user opted for permutation test and probability level, size of genome scan and method used for the calculation of the QTL genotype probabilities, QTL mapping method and QTL model used. Clicking a chromosome with a QTL, after visually determining the chromosome with a significant QTL relative to the LOD threshold value in the legend, takes the user to a QTL detail page (Figure 2).
Linking QTLs to the genome and genetic maps
A crucial feature of the solQTL tool is that it links the putative QTL regions to relevant genetic and genomic data in the SGN database. In this subsection, we describe the options for linking the predicted QTLs with other datasets both at SGN and at external databases using bioinformatic tools such as the SGN Comparative Map  Viewer and GBrowse  tools.
On the QTL detail page (Figure 2) for each QTL, the peak and flanking markers with their corresponding positions on the respective genetic map for the population are shown. The 95% QTL confidence interval is linked to the corresponding linkage group page where one can view all the markers in the linkage group with the QTL segment zoomed out. On the linkage group webpage, users can perform comparative map analysis between the QTL region of interest and physical and genetic maps from several Solanaceae species in the database using the SGN Comparative Map Viewer, as long as the QTL flanking markers are mapped to the other maps. Marker-dense genetic maps available in SGN for comparative analysis include the S. lycopersicum LA925 x S. pennellii LA716 F2.2000 map  which contains more than 2500 Restriction Fragment Length Polymorphism (RFLP), Cleaved Amplified Polymorphic Sequences (CAPS) and Conserved Ortholog Set (COS) markers (Figure 3). The SGN Comparative Map Viewer is described in  and additional help and documentation for this tool can be found at http://solgenomics.net/help/cview.pl.
From the QTL detail page, the markers are linked to genomic regions visible in SGN's installation of GBrowse, provided that there is a match, where one can view the corresponding genomic sequences and annotations such as gene, mRNA, protein, and CDS models, ESTs, and unigene alignments and experimentally characterized genes from tomato or other organisms (Figure 4). Documentation using GBowse can be found in http://solgenomics.net/gbrowse/static/general_help.html and [42, 44].
The peak and flanking markers are also cross-referenced to their respective SGN detail pages where one can find, among other data, the primer sequences, map positions on any other genetic maps and matching BAC clones. For COS markers, one can also find matching Arabidopsis BAC clones and protein hits from NCBI http://www.ncbi.nlm.nih.gov/.
solQTL tool and data accessibility
solQTL is accessible at http://solgenomics.net/qtl/ or from the 'search' and 'tools' pull down menus on the SGN toolbar, located at the top of every SGN page. From the solQTL database frontpage, users interested in submitting new QTL data and in performing QTL mapping and linking to relevant datasets in SGN can follow the self-guiding 'Submit new QTL data' link. From the same solQTL front page, users intending to query the database for QTLs for traits of interest can also: (1) browse the alphabetically indexed list of traits; (2) use the search box and query the database using relevant keywords; or (3) browse through the list of QTL populations for traits of interest. In the cases of option 1 and 2, the user is directed to the trait's detail page where, among other data, QTL populations assayed for the trait are listed. Following the link to the population of interest initiates solQTL's on-the-fly QTL mapping and generates the QTL detail page for the trait of interest in that QTL population. In the case of option 3, users can select a trait of their interest evaluated in any of the populations and run the QTL mapping analysis.
The population genetic map can be accessed either by following a genetic map link at the population or QTL detail pages or by browsing for it from the 'maps' pull down menu on any SGN page.
Beyond QTL: comparative analysis and identification of candidate genes
Linking QTLs to corresponding regions in genomes and their annotations is pursued by several biological databases [24, 27]. Most databases store pre-computed QTL coordinates, extracted from published literature, and are mostly populated manually by in-house curators. In contrast, SGN has implemented a new approach that dynamically links phenotypes to genomes by implementing the R/QTL mapping software  and thus, allowing researchers to upload their raw phenotype and genotype QTL data to the SGN database and set the statistical parameters for the analysis using a user-friendly web interface. To our knowledge, none of these features are available at other databases.
R/QTL is a robust and widely used (with about 450 citations as of this publication, GoogleScholar) open source QTL mapping software. We have implemented its analytical functions such as single marker analysis and interval mapping methods (maximum likelihood [32, 33], Haley-Knott regression , and Multiple Imputation ) for diploid F2 intercross and backcross populations. solQTL is publicly and freely available, and SGN maintains and upgrades the software so that the latest versions are available to run on SGN web servers. In this sense, from the users perspective, the data are stored and analyzed in the 'cloud', as such internet-based services are sometimes referred to, although in this case the 'cloud' is currently limited to the SGN cluster infrastructure. The user-friendly web interface for uploading raw data and setting the statistical parameters combined with the on the-fly QTL analysis by users, in private or public mode, further reinforces the "community annotation" philosophy SGN has pioneered for its locus and phenotype databases .
SGN is a clade-oriented database with a computational infrastructure and capacity to store and make accessible data, in a comparative manner, for the Solanaceae family and related Asterids. It is a primary repository for genetic, phenotype, genome and expression data for the Solanaceae family . These datasets are deeply and extensively cross-referenced. As a result, provided that the QTL flanking and/or peak markers exist with their sequences in the SGN marker database, the solQTL tool links dynamically predicted QTLs to corresponding relevant datasets and builds upon the constantly growing networking of emerging datasets in SGN and external databases.
Predicted QTLs are seamlessly cross-linked to the SGN Comparative Map Viewer  and SGN's instance of GBrowse . The comparative map viewer enables comparative genetic analysis between QTLs of interest and corresponding regions in genetic and physical maps of the same or different Solanaceae species and Arabidopsis. Corresponding regions in other maps, for example the tomato F2.2000 map , with more dense and informative markers such as COS, help identify more tightly linked genetic markers and possibly matching genomic clones or regions. QTLs are currently linked to the GBrowse on a marker-by-marker basis. The GBrowse is also searchable using other genetic features such as unigenes which might enable investigators to infer some missing links from genetic markers to the genome. We anticipate, in the near future, further progress in the assembly of the genomes SGN hosts and search features under development will allow searching for genomic regions matching an entire QTL region. At this stage once a genome region matches a flanking or peak QTL marker, one can zoom in/out (bi-directionally) or scroll the genomic segment one-directionally to explore genomic annotations such as gene, mRNA, CDS and protein models, unigenes, ESTs and experimentally characterized loci. Currently, SGN's GBrowse database houses the Solanum lycopersicum cv Heinz 1706 genome sequence and annotations by ITAG . More genomes are in the pipeline to be sequenced and annotated and/or added to SGN's GBrowse and cross-referenced to solQTL.
The solQTL tool allows users to use non-controlled vocabularies trait names. However, users are highly encouraged to use controlled vocabularies for their trait names to facilitate cross-study comparisons for QTLs of the same traits and to fully exploit expression and phenotype data annotated with controlled vocabulary such as the Gene Ontology (GO)  and Plant Ontology (PO)  terms. SGN in collaboration with the Solanaceae breeders continues to develop and maintain a controlled vocabulary for Solanaceae traits, called Solanaceae Phenotype Ontology, accessible at http://solgenomics.net/tools/onto/.
Looking toward a longer time horizon, the SOL-100 project has an ambitious plan to sequence and make publicly available genomes of 100 or more Solanaceae species . A single nucleotide polymorphism (SNP) database and viewer is also under development at SGN (unpublished). As these projects bear fruit, solQTL will be a window linking phenotypes to this rich genomic data and analytical tools. We envision this integration and analysis of a wide variety of data in SGN and cross-referenced external databases will have an important role to play in identification of candidate genes for further functional analysis and gene cloning, identification of more closely linked QTL markers for marker-assisted breeding and understanding the molecular and physiological underlying phenotypic variation for traits of interest.
Conclusion and future plans
As genome sequencing becomes less costly, the number of organisms with sequenced and annotated genomes will increase. Linking the emerging wealth of genome data to phenotypes in order to help elucidate the molecular basis for phenotypic variation of quantitative traits will be a major challenge. We believe tools such as solQTL will help facilitate the dynamic linking of phenotypes to genomes for the Solanceae species.
At this stage, solQTL implements a limited number of QTL mapping capabilities of the R/QTL software. In the near future, we plan to implement the full capacity of the R/QTL tool such as genetic map estimation, two-dimensional two-QTL scan, QTL mapping for out-crossing and recombinant inbred lines, and account for covariates such as environment. We plan to extend its analysis capabilities to support multiple QTL mapping, QTL effects calculation, and interactive statistical analysis where users can predict QTLs using several combinations of statistical parameters.
We also plan to implement tools to allow meta-QTL analysis across multiple QTL populations. To enable communication between QTL databases and access to the raw and analyzed QTL data at SGN programmatically, we will generate QTL output in XML format as specified by the emerging standards of Minimum Information for QTLs and Association Studies (MIQAS).
Availability and requirements
All SGN source code and database schemas are open source and available for download at http://github.com/solgenomics. Raw phenotype and genotype data can also be freely downloaded from the SGN web page of the QTL population of interest, provided that the owner has released the data for the public use. solQTL works with all web browsers and good internet connection.
Fridman E, Pleban T, Zamir D: A recombination hotspot delimits a wild-species quantitative trait locus for tomato sugar content to 484 bp within an invertase gene. Proc Natl Acad Sci USA 2000, 97(9):4718–4723. 10.1073/pnas.97.9.4718
Causse M, Duffe P, Gomez MC, Buret M, Damidaux R, Zamir D, Gur A, Chevalier C, Lemaire-Chamley M, Rothan C: A genetic map of candidate genes and QTLs involved in tomato fruit size and composition. J Exp Bot 2004, 55(403):1671–1685. 10.1093/jxb/erh207
Pflieger S, Lefebvre V, Causse M: The candidate gene approach in plant genetics: a review. Molecular Breeding 2001, 7: 275–291. 10.1023/A:1011605013259
Zou F: QTL mapping in intercross and backcross populations. Methods Mol Biol 2009, 573: 157–173. full_text
Schauer N, Semel Y, Roessner U, Gur A, Balbo I, Carrari F, Pleban T, Perez-Melis A, Bruedigam C, Kopka J, Willmitzer L, Zamir D, Fernie AR: Comprehensive metabolic profiling and phenotyping of interspecific introgression lines for tomato improvement. Nat Biotechnol 2006, 24(4):447–454. 10.1038/nbt1192
Alpert KB, Tanksley SD: High-resolution mapping and isolation of a yeast artificial chromosome contig containing fw2.2: a major fruit weight quantitative trait locus in tomato. Proc Natl Acad Sci USA 1996, 93(26):15503–15507. 10.1073/pnas.93.26.15503
Xu X, Martin B, Comstock JP, Vision TJ, Tauer CG, Zhao B, Pausch RC, Knapp S: Fine mapping a QTL for carbon isotope composition in tomato. Theor Appl Genet 2008, 117(2):221–233. 10.1007/s00122-008-0767-6
Eshed Y, Zamir D: An introgression line population of Lycopersicon pennellii in the cultivated tomato enables the identification and fine mapping of yield associated QTL. Genetics 1995, 141: 1147–1162.
Alonso-Blanco C, Koornneef M: Naturally occurring variation in Arabidopsis: an underexploited resource for plant genetics. Trends Plant Sci 2000, 5(1):22–29. 10.1016/S1360-1385(99)01510-1
Price AH: Believe it or not, QTLs are accurate! Trends Plant Sci 2006, 11(5):213–216. 10.1016/j.tplants.2006.03.006
Bermudez L, Urias U, Milstein D, Kamenetzky L, Asis R, Fernie AR, Van Sluys MA, Carrari F, Rossi M: A candidate gene survey of quantitative trait loci affecting chemical composition in tomato fruit. J Exp Bot 2008, 59(10):2875–2890. 10.1093/jxb/ern146
Flint J, Mott R: Finding the molecular basis of quantitative traits: successes and pitfalls. Nat Rev Genet 2001, 2(6):437–445. 10.1038/35076585
Sharma A, Zhang L, Nino-Liu D, Ashrafi H, Foolad MR: A Solanum lycopersicum x Solanum pimpinellifolium Linkage Map of Tomato Displaying Genomic Locations of R-Genes, RGAs, and Candidate Resistance/Defense-Response ESTs. Int J Plant Genomics 2008, 2008: 926090.
Mueller LA, Solow TH, Taylor N, Skwarecki B, Buels R, Binns J, Lin C, Wright MH, Ahrens R, Wang Y, Herbst EV, Keyder ER, Menda N, Zamir D, Tanksley SD: The SOL Genomics Network: a comparative resource for Solanaceae biology and beyond. Plant Physiol 2005, 138(3):1310–1317. 10.1104/pp.105.060707
Broman KW, Wu H, Sen S, Churchill GA: R/qtl: QTL mapping in experimental crosses. Bioinformatics 2003, 19(7):889–890. 10.1093/bioinformatics/btg112
Arcade A, Labourdette A, Falque M, Mangin B, Chardon F, Charcosset A, Joets J: BioMercator: integrating genetic maps and QTL towards discovery of candidate genes. Bioinformatics 2004, 20(14):2324–2326. 10.1093/bioinformatics/bth230
Frary A, Nesbitt TC, Grandillo S, Knaap E, Cong B, Liu J, Meller J, Elber R, Alpert KB, Tanksley SD: Fw2.2: a Quantitative Trait Locus Key to the Evolution of Tomato Fruit Size. Science 2000, 289(5476):85–88. 10.1126/science.289.5476.85
Liu J, Van Eck J, Cong B, Tanksley SD: A new class of regulatory genes underlying the cause of pear-shaped tomato fruit. Proc Natl Acad Sci USA 2002, 99(20):13302–13306. 10.1073/pnas.162485999
Yamamoto T, Kuboki Y, Lin SY, Sasaki T, Yano M: Fine mapping of quantitative trait loci Hd-1, Hd-2 and Hd-3, controlling heading date of rice, as single Mendelian factors. TAG 1998, 97(1–2):37–44. 10.1007/s001220050864
Tanksley SD, Young ND, Paterson AH, Bonierbale MW: RFLP Mapping in Plant Breeding: New Tools for an Old Science. Nat Biotechnol 1989, 7: 257–264. 10.1038/nbt0389-257
Mueller LA, Lankhorst RK, Tanksley SD, Giovannoni JJ: A Snapshot of the Emerging Tomato Genome Sequence. The Plant Genome 2009, 2: 78–92. 10.3835/plantgenome2008.08.0005
Klein Lankhorst R: Sol 100. Sol Newsletter March 2010, 26: 1–3.
Liang C, Jaiswal P, Hebbard C, Avraham S, Buckler ES, Casstevens T, Hurwitz B, McCouch S, Ni J, Pujar A, Ravenscroft D, Ren L, Spooner W, Tecle I, Thomason J, Tung CW, Wei X, Yap I, Youens-Clark K, Ware D, Stein L: Gramene: a growing plant comparative genomics resource. Nucleic Acids Res 2008, 36(Database issue):D947–53.
Ni J, Pujar A, Youens-Clark K, Yap I, Jaiswal P, Tecle I, Tung CW, Ren L, Spooner W, Wei X, Avraham S, Ware D, Stein L, McCouch S: Gramene QTL database: development, content and applications. Database (Oxford) 2009, 2009: bap005.
Lawrence CJ, Harper LC, Schaeffer ML, Sen TZ, Seigfried TE, Campbell DA: MaizeGDB: The Maize Model Organism Database for Basic, Translational, and Applied Research. Int J Plant Genomics 2008, 2008: 496957.
Hu ZL, Fritz ER, Reecy JM: AnimalQTLdb: a livestock QTL database tool set for positional QTL information mining and beyond. Nucleic Acids Res 2007, 35(Database issue):D604–9. 10.1093/nar/gkl946
Hu ZL, Reecy JM: Animal QTLdb: beyond a repository. A public platform for QTL comparisons and integration with diverse types of structural genomic information. Mamm Genome 2007, 18(1):1–4. 10.1007/s00335-006-0105-8
Gonzalo MJ, van der Knaap E: A comparative analysis into the genetic bases of morphology in tomato varieties exhibiting elongated fruit shape. Theor Appl Genet 2008, 116(5):647–656. 10.1007/s00122-007-0698-7
Brewer MT, Moyseenko JB, Monforte AJ, van der Knaap E: Morphological variation in tomato: a comprehensive study of quantitative trait loci controlling fruit shape and development. J Exp Bot 2007, 58(6):1339–1349. 10.1093/jxb/erl301
van der Knaap E, Tanksley SD: The making of a bell pepper-shaped tomato fruit: identification of loci controlling fruit morphology in Yellow Stuffer tomato. Theor Appl Genet 2003, 107(1):139–147.
Broman KW, Sen S: Single-QTL Analysis. In A Guide to QTL mapping with R/Qtl, Statistics for Biology and Health. Edited by: Broman KW, Sen S. Springer; 2009:75–133. full_text
Dempster AP, Laird NM, Rubin DB: Maximum likelihood from incomplete data via the EM algorithm. J Roy Statist Soc 1977, 39(B):1–38.
Lander ES, Botstein D: Mapping mendelian factors underlying quantitative traits using RFLP linkage maps. Genetics 1989, 121(1):185–199.
Haley CS, Knott SA: A simple regression method for mapping quantitative trait loci in line crosses using flanking markers. Heredity 1992, 69(4):315–324.
Sen S, Churchill GA: A statistical framework for quantitative trait mapping. Genetics 2001, 159(1):371–87.
Baum LE, Petrie G, Weiss N: A maximization technique occurring in the statistical analysis of probabilistic functions of Markov chains. Ann Math Stat 1970, 41: 164–171. 10.1214/aoms/1177697196
Churchill GA, Doerge RW: Empirical threshold values for quantitative trait mapping. Genetics 1994, 138(3):963–971.
Manichaikul A, Dupuis J, Sen S, Broman KW: Poor performance of bootstrap confidence intervals for the location of a quantitative trait locus. Genetics 2006, 174(1):481–489. 10.1534/genetics.106.061549
Mueller LA, Mills AA, Skwarecki B, Buels RM, Menda N, Tanksley SD: The SGN comparative map viewer. Bioinformatics 2008, 24(3):422–423. 10.1093/bioinformatics/btm597
R Development Core Team: R: A language and environment for statistical computing. 2008., 2.7.1:
Zhao JH, Tan Q: Integrated analysis of genetic data with R. Hum Genomics 2006, 2(4):258–265.
Stein LD, Mungall C, Shu S, Caudy M, Mangone M, Day A, Nickerson E, Stajich JE, Harris TW, Arva A, Lewis S: The generic genome browser: a building block for a model organism system database. Genome Res 2002, 12(10):1599–1610. 10.1101/gr.403602
Fulton TM, Van der Hoeven R, Eannetta NT, Tanksley SD: Identification, analysis, and utilization of conserved ortholog set markers for comparative genomics in higher plants. Plant Cell 2002, 14(7):1457–1467. 10.1105/tpc.010479
Donlin MJ: Using the Generic Genome Browser (GBrowse). Curr Protoc Bioinformatics 2009., Chapter 9: Unit 9.9 Unit 9.9
Menda N, Buels RM, Tecle I, Mueller LA: A community-based annotation framework for linking solanaceae genomes with phenomes. Plant Physiol 2008, 147(4):1788–1799. 10.1104/pp.108.119560
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet 2000, 25(1):25–29. 10.1038/75556
Avraham S, Tung CW, Ilic K, Jaiswal P, Kellogg EA, McCouch S, Pujar A, Reiser L, Rhee SY, Sachs MM, Schaeffer M, Stein L, Stevens P, Vincent L, Zapata F, Ware D: The Plant Ontology Database: a community resource for plant structure and developmental stages controlled vocabulary and annotations. Nucleic Acids Res 2008, 36(Database issue):D449–54.
The Minimum Information for QTLs and Association Studies (MIQAS)[http://miqas.sourceforge.net/]
The project was funded by the USDA CSREES. For part of the duration of this project the first author of this manuscript has been a visiting scientist at Wageningen University, Plant Breeding and Genetics Department, The Netherlands. Many thanks to Dr. Richard Finkers and his lab for hosting the first author.
IYT designed the workflow, implemented the R/QTL code and wrote most of the code and part of the solQTL database schema or implemented modules contributed by NM, RMB, LAM and CPAN authors. NM, RMB, and LAM contributed in code writing, database schema development and discussion. LAM oversaw the development of the project. EVDK contributed the QTL data. The manuscript was drafted by IYT and all authors reviewed and contributed up to the final version.
Electronic supplementary material
About this article
- Quantitative Trait Locus
- Quantitative Trait Locus Mapping
- Quantitative Trait Locus Region
- Peak Quantitative Trait Locus
- Quantitative Trait Nucleotide