K-seq, an affordable, reliable, and open Klenow NGS-based genotyping technology

[EN] Background: K-seq, a new genotyping methodology based on the amplification of genomic regions using two steps of Klenow amplification with short oligonucleotides, followed by standard PCR and Illumina sequencing, is presented. The protocol was accompanied by software developed to aid with primer set design. Results: As the first examples, K-seq in species as diverse as tomato, dog and wheat was developed. K-seq provided genetic distances similar to those based on WGS in dogs. Experiments comparing K-seq and GBS in tomato showed similar genetic results, although K-seq had the advantage of finding more SNPs for the same number of Illumina reads. The technology reproducibility was tested with two independent runs of the tomato samples, and the correlation coefficient of the SNP coverages between samples was 0.8 and the genotype match was above 94%. K-seq also proved to be useful in polyploid species. The wheat samples generated specific markers for all subgenomes, and the SNPs generated from the diploid ancestors were located in the expected subgenome with accuracies greater than 80%. Conclusion: K-seq is an open, patent-unencumbered, easy-to-set-up, cost-effective and reliable technology ready to be used by any molecular biology laboratory without special equipment in many genetic studies.This work was supported by the University Polytechnic of Valencia, Grant Number 20180051 "Desarrollo de herramientas para la identificacion de genes y loci de interes en la mejora genetica del tomate y otras horticolas".Ziarsolo, P.; Hasing, T.; Hilario, R.; García-Carpintero, V.; Blanca Postigo, JM.; Bombarely, A.; Cañizares Sales, J. (2021). K-seq, an affordable, reliable, and open Klenow NGS-based genotyping technology. Plant Methods. 17(1):1-11. https://doi.org/10.1186/s13007-021-00733-6S11117

Haplotype analyses reveal novel insights into tomato history and domestication driven by long-distance migrations and latitudinal adaptations

A novel haplotype-based approach that uses Procrustes analysis and automatic classification was used to provide further insights into tomato history and domestication. Agrarian societies domesticated species of interest by introducing complex genetic modifications. For tomatoes, two species, one of which had two botanical varieties, are thought to be involved in its domestication: the fully wild Solanum pimpinellifolium (SP), the wild and semi-domesticated Solanum lycopersicum var. cerasiforme (SLC) and the cultivated S. l. var. lycopersicum (SLL). The Procrustes approach showed that SP evolved into SLC during a gradual migration from the Peruvian deserts to the Mexican rainforests and that Peruvian and Ecuadorian SLC populations were the result of more recent hybridizations. Our model was supported by independent evidence, including ecological data from the accession collection site and morphological data. Furthermore, we showed that photosynthesis-, and f lowering time-related genes were selected during the latitudinal migrations

European traditional tomatoes galore: a result of farmers' selection of a few diversity-rich loci

A comprehensive collection of 1254 tomato accessions, corresponding to European traditional and modern varieties, early domesticated varieties, and wild relatives, was analyzed by genotyping by sequencing. A continuous genetic gradient between the traditional and modern varieties was observed. European traditional tomatoes displayed very low genetic diversity, with only 298 polymorphic loci (95% threshold) out of 64 943 total variants. European traditional tomatoes could be classified into several genetic groups. Two main clusters consisting of Spanish and Italian accessions showed higher genetic diversity than the remaining varieties, suggesting that these regions might be independent secondary centers of diversity with a different history. Other varieties seem to be the result of a more recent complex pattern of migrations and hybridizations among the European regions. Several polymorphic loci were associated in a genome-wide association study with fruit morphological traits in the European traditional collection. The corresponding alleles were found to contribute to the distinctive phenotypic characteristic of the genetic varietal groups. The few highly polymorphic loci associated with morphological traits in an otherwise a low-diversity population suggests a history of balancing selection, in which tomato farmers likely maintained the morphological variation by inadvertently applying a high selective pressure within different varietal types

Quantitative genetic analysis of floral traits shows current limits but potential evolution in the wild

The vast variation in floral traits across angiosperms is often interpreted as the result of adaptation to pollinators. However, studies in wild populations often find no evidence of pollinator-mediated selection on flowers. Evolutionary theory predicts this could be the outcome of periods of stasis under stable conditions, followed by shorter periods of pollinator change that provide selection for innovative phenotypes. We asked if periods of stasis are caused by stabilizing selection, absence of other forms of selection or by low trait ability to respond even if selection is present. We studied a plant predominantly pollinated by one bee species across its range. We measured heritability and evolvability of traits, using genome-wide relatedness in a large wild population, and combined this with estimates of selection on the same individuals. We found evidence for both stabilizing selection and low trait heritability as potential explanations for stasis in flowers. The area of the standard petal is under stabilizing selection, but the variability is not heritable. A separate trait, floral weight, presents high heritability, but is not currently under selection. We show how a simple pollination environment coincides with the absence of current prerequisites for adaptive evolutionary change, while heritable variation remains to respond to future selection pressures

Discovery of a major QTL controlling trichome IV density in tomato using K-seq genotyping

Trichomes are a common morphological defense against pests, in particular, type IV glandular trichomes have been associated with resistance against different invertebrates. Cultivated tomatoes usually lack or have a very low density of type IV trichomes. Therefore, for sustainable management of this crop, breeding programs could incorporate some natural defense mechanisms, such as those afforded by trichomes, present in certain Solanum species. We have identified a S. pimpinellifolium accession with very high density of this type of trichomes. This accession was crossed with a S. lycopersicum var. cerasiforme and a S. lycopersicum var. lycopersicum accessions, and the two resulting F2 populations have been characterized and genotyped using a new genotyping methodology, K-seq. We have been able to build an ultra-dense genetic map with 147,326 SNP markers with an average distance between markers of 0.2 cm that has allowed us to perform a detailed mapping. We have used two different families and two different approaches, QTL mapping and QTL-seq, to identify several QTLs implicated in the control of trichome type IV developed in this accession on the chromosomes 5, 6, 9 and 11. The QTL located on chromosome 9 is a major QTL that has not been previously reported in S. pimpinellifolium. This QTL could be easily introgressed in cultivated tomato due to the close genetic relationship between both species

Exploiting the diversity of tomato: the development of a phenotypically and genetically detailed germplasm collection

A collection of 163 accessions, including Solanum pimpinellifolium, Solanum lycopersicum var. cerasiforme and Solanum lycopersicum var. lycopersicum, was selected to represent the genetic and morphological variability of tomato at its centers of origin and domestication: Andean regions of Peru and Ecuador and Mesoamerica. The collection is enriched with S. lycopersicum var. cerasiforme from the Amazonian region that has not been analyzed previously nor used extensively. The collection has been morphologically characterized showing diversity for fruit, flower and vegetative traits. Their genomes were sequenced in the Varitome project and are publicly available (solgenomics.net/projects/varitome). The identified SNPs have been annotated with respect to their impact and a total number of 37,974 out of 19,364,146 SNPs have been described as high impact by the SnpEeff analysis. GWAS has shown associations for different traits, demonstrating the potential of this collection for this kind of analysis. We have not only identified known QTLs and genes, but also new regions associated with traits such as fruit color, number of flowers per inflorescence or inflorescence architecture. To speed up and facilitate the use of this information, F2 populations were constructed by crossing the whole collection with three different parents. This F2 collection is useful for testing SNPs identified by GWAs, selection sweeps or any other candidate gene. All data is available on Solanaceae Genomics Network and the accession and F2 seeds are freely available at COMAV and at TGRC genebanks. All these resources together make this collection a good candidate for genetic studies

GB3.0: a platform for plant bio-design that connects functional DNA elements with associated biological data

This is a pre-copyedited, author-produced version of an article accepted for publication in Nucleic Acids Research following peer review. The version of record Vázquez-Vilar, M.; Quijano-Rubio, A.; Fernandez Del Carmen, MA.; Sarrion-Perdigones, A.; Ochoa-Fernández, R.; Ziarsolo Areitioaurtena, P.; Blanca Postigo, JM.... (2017). GB3.0: a platform for plant bio-design that connects functional DNA elements with associated biological data. Nucleic Acids Research. 45(4):2196-2209. doi:10.1093/nar/gkw1326 is available online at: http://doi.org/10.1093/nar/gkw1326.[EN] Modular DNA assembly simplifies multigene engineering in Plant Synthetic Biology. Furthermore, the recent adoption of a common syntax to facilitate the exchange of plant DNA parts (phytobricks) is a promising strategy to speed up genetic engineering. Following this lead, here, we present a platform for plant biodesign that incorporates functional descriptions of phytobricks obtained under pre-defined experimental conditions, and systematically registers the resulting information as metadata for documentation. To facilitate the handling of functional descriptions, we developed a new version (v3.0) of the GoldenBraid (GB) webtool that integrates the experimental data and displays it in the form of datasheets. We report the use of the Luciferase/Renilla (Luc/Ren) transient agroinfiltration assay in Nicotiana benthamiana as a standard to estimate relative transcriptional activities conferred by regulatory phytobricks, and show the consistency and reproducibility of this method in the characterization of a synthetic phytobrick based on the CaMV35S promoter. Furthermore, we illustrate the potential for combinatorial optimization and incremental innovation of the GB3.0 platform in two separate examples, (i) the development of a collection of orthogonal transcriptional regulators based on phiC31 integrase and (ii) the design of a small genetic circuit that connects a glucocorticoid switch to a MYB/bHLH transcriptional activation module.Spanish Ministry of Economy and Competitiveness [BIO2013-42193-R and BIO2016-78601-R projects to A.G. and D.O.]. Funding for open access charge: Spanish Ministry of Economy and Competitiveness [BIO2013-42193-R and BIO2016-78601-R projects to A.G. and D.O.].Vázquez-Vilar, M.; Quijano-Rubio, A.; Fernández Del Carmen, MA.; Sarrion-Perdigones, A.; Ochoa-Fernández, R.; Ziarsolo Areitioaurtena, P.; Blanca Postigo, JM.... (2017). GB3.0: a platform for plant bio-design that connects functional DNA elements with associated biological data. Nucleic Acids Research. 45(4):2196-2209. https://doi.org/10.1093/nar/gkw1326S21962209454Dalal, J., Yalamanchili, R., La Hovary, C., Ji, M., Rodriguez-Welsh, M., Aslett, D., … Qu, R. (2015). A novel gateway-compatible binary vector series (PC-GW) for flexible cloning of multiple genes for genetic transformation of plants. Plasmid, 81, 55-62. doi:10.1016/j.plasmid.2015.06.003Cha-aim, K., Fukunaga, T., Hoshida, H., & Akada, R. (2009). Reliable fusion PCR mediated by GC-rich overlap sequences. Gene, 434(1-2), 43-49. doi:10.1016/j.gene.2008.12.014Nour-Eldin, H. H., Geu-Flores, F., & Halkier, B. A. (2010). USER Cloning and USER Fusion: The Ideal Cloning Techniques for Small and Big Laboratories. Methods in Molecular Biology, 185-200. doi:10.1007/978-1-60761-723-5_13Engler, C., Kandzia, R., & Marillonnet, S. (2008). A One Pot, One Step, Precision Cloning Method with High Throughput Capability. PLoS ONE, 3(11), e3647. doi:10.1371/journal.pone.0003647Blake, W. J., Chapman, B. A., Zindal, A., Lee, M. E., Lippow, S. M., & Baynes, B. M. (2010). Pairwise selection assembly for sequence-independent construction of long-length DNA. Nucleic Acids Research, 38(8), 2594-2602. doi:10.1093/nar/gkq123De Paoli, H. C., Tuskan, G. A., & Yang, X. (2016). An innovative platform for quick and flexible joining of assorted DNA fragments. Scientific Reports, 6(1). doi:10.1038/srep19278Engler, C., Gruetzner, R., Kandzia, R., & Marillonnet, S. (2009). Golden Gate Shuffling: A One-Pot DNA Shuffling Method Based on Type IIs Restriction Enzymes. PLoS ONE, 4(5), e5553. doi:10.1371/journal.pone.0005553Weber, E., Engler, C., Gruetzner, R., Werner, S., & Marillonnet, S. (2011). A Modular Cloning System for Standardized Assembly of Multigene Constructs. PLoS ONE, 6(2), e16765. doi:10.1371/journal.pone.0016765Sarrion-Perdigones, A., Falconi, E. E., Zandalinas, S. I., Juárez, P., Fernández-del-Carmen, A., Granell, A., & Orzaez, D. (2011). GoldenBraid: An Iterative Cloning System for Standardized Assembly of Reusable Genetic Modules. PLoS ONE, 6(7), e21622. doi:10.1371/journal.pone.0021622Patron, N. J., Orzaez, D., Marillonnet, S., Warzecha, H., Matthewman, C., Youles, M., … Rogers, C. (2015). Standards for plant synthetic biology: a common syntax for exchange ofDNAparts. New Phytologist, 208(1), 13-19. doi:10.1111/nph.13532Liu, W., & Stewart, C. N. (2015). Plant synthetic biology. Trends in Plant Science, 20(5), 309-317. doi:10.1016/j.tplants.2015.02.004Wang, Y.-H., Wei, K. Y., & Smolke, C. D. (2013). Synthetic Biology: Advancing the Design of Diverse Genetic Systems. Annual Review of Chemical and Biomolecular Engineering, 4(1), 69-102. doi:10.1146/annurev-chembioeng-061312-103351Engler, C., Youles, M., Gruetzner, R., Ehnert, T.-M., Werner, S., Jones, J. D. G., … Marillonnet, S. (2014). A Golden Gate Modular Cloning Toolbox for Plants. ACS Synthetic Biology, 3(11), 839-843. doi:10.1021/sb4001504Canton, B., Labno, A., & Endy, D. (2008). Refinement and standardization of synthetic biological parts and devices. Nature Biotechnology, 26(7), 787-793. doi:10.1038/nbt1413Sarrion-Perdigones, A., Vazquez-Vilar, M., Palaci, J., Castelijns, B., Forment, J., Ziarsolo, P., … Orzaez, D. (2013). GoldenBraid 2.0: A Comprehensive DNA Assembly Framework for Plant Synthetic Biology. PLANT PHYSIOLOGY, 162(3), 1618-1631. doi:10.1104/pp.113.217661Orzaez, D., Medina, A., Torre, S., Fernández-Moreno, J. P., Rambla, J. L., Fernández-del-Carmen, A., … Granell, A. (2009). A Visual Reporter System for Virus-Induced Gene Silencing in Tomato Fruit Based on Anthocyanin Accumulation. Plant Physiology, 150(3), 1122-1134. doi:10.1104/pp.109.139006Yoo, S.-D., Cho, Y.-H., & Sheen, J. (2007). Arabidopsis mesophyll protoplasts: a versatile cell system for transient gene expression analysis. Nature Protocols, 2(7), 1565-1572. doi:10.1038/nprot.2007.199Schneider, C. A., Rasband, W. S., & Eliceiri, K. W. (2012). NIH Image to ImageJ: 25 years of image analysis. Nature Methods, 9(7), 671-675. doi:10.1038/nmeth.2089Gutiérrez, S., Yvon, M., Thébaud, G., Monsion, B., Michalakis, Y., & Blanc, S. (2010). Dynamics of the Multiplicity of Cellular Infection in a Plant Virus. PLoS Pathogens, 6(9), e1001113. doi:10.1371/journal.ppat.1001113Vazquez-Vilar, M., Sarrion-Perdigones, A., Ziarsolo, P., Blanca, J., Granell, A., & Orzaez, D. (2015). Software-Assisted Stacking of Gene Modules Using GoldenBraid 2.0 DNA-Assembly Framework. Plant Functional Genomics, 399-420. doi:10.1007/978-1-4939-2444-8_20Vazquez-Vilar, M., Bernabé-Orts, J. M., Fernandez-del-Carmen, A., Ziarsolo, P., Blanca, J., Granell, A., & Orzaez, D. (2016). A modular toolbox for gRNA–Cas9 genome engineering in plants based on the GoldenBraid standard. Plant Methods, 12(1). doi:10.1186/s13007-016-0101-2Quinn, J. Y., Cox, R. S., Adler, A., Beal, J., Bhatia, S., Cai, Y., … Sauro, H. M. (2015). SBOL Visual: A Graphical Language for Genetic Designs. PLOS Biology, 13(12), e1002310. doi:10.1371/journal.pbio.1002310Thorpe, H. M., Wilson, S. E., & Smith, M. C. M. (2000). Control of directionality in the site-specific recombination system of the Streptomyces phage phiC31. Molecular Microbiology, 38(2), 232-241. doi:10.1046/j.1365-2958.2000.02142.xKhaleel, T., Younger, E., McEwan, A. R., Varghese, A. S., & Smith, M. C. M. (2011). A phage protein that binds φC31 integrase to switch its directionality. Molecular Microbiology, 80(6), 1450-1463. doi:10.1111/j.1365-2958.2011.07696.xOrt, D. R., Merchant, S. S., Alric, J., Barkan, A., Blankenship, R. E., Bock, R., … Zhu, X. G. (2015). Redesigning photosynthesis to sustainably meet global food and bioenergy demand. Proceedings of the National Academy of Sciences, 112(28), 8529-8536. doi:10.1073/pnas.1424031112Smolke, C. D. (2009). Building outside of the box: iGEM and the BioBricks Foundation. Nature Biotechnology, 27(12), 1099-1102. doi:10.1038/nbt1209-1099Cambray, G., Mutalik, V. K., & Arkin, A. P. (2011). Toward rational design of bacterial genomes. Current Opinion in Microbiology, 14(5), 624-630. doi:10.1016/j.mib.2011.08.001Mutalik, V. K., Guimaraes, J. C., Cambray, G., Mai, Q.-A., Christoffersen, M. J., Martin, L., … Arkin, A. P. (2013). Quantitative estimation of activity and quality for collections of functional genetic elements. Nature Methods, 10(4), 347-353. doi:10.1038/nmeth.2403Cambray, G., Guimaraes, J. C., Mutalik, V. K., Lam, C., Mai, Q.-A., Thimmaiah, T., … Endy, D. (2013). Measurement and modeling of intrinsic transcription terminators. Nucleic Acids Research, 41(9), 5139-5148. doi:10.1093/nar/gkt163Mutalik, V. K., Guimaraes, J. C., Cambray, G., Lam, C., Christoffersen, M. J., Mai, Q.-A., … Endy, D. (2013). Precise and reliable gene expression via standard transcription and translation initiation elements. Nature Methods, 10(4), 354-360. doi:10.1038/nmeth.2404Bartley, B., Beal, J., Clancy, K., Misirli, G., Roehner, N., Oberortner, E., … Sauro, H. (2015). Synthetic Biology Open Language (SBOL) Version 2.0.0. Journal of Integrative Bioinformatics, 12(2). doi:10.1515/jib-2015-272Angov, E. (2011). Codon usage: Nature’s roadmap to expression and folding of proteins. Biotechnology Journal, 6(6), 650-659. doi:10.1002/biot.201000332Hillson, N. J., Rosengarten, R. D., & Keasling, J. D. (2011). j5 DNA Assembly Design Automation Software. ACS Synthetic Biology, 1(1), 14-21. doi:10.1021/sb2000116Brazma, A., Hingamp, P., Quackenbush, J., Sherlock, G., Spellman, P., Stoeckert, C., … Vingron, M. (2001). Minimum information about a microarray experiment (MIAME)—toward standards for microarray data. Nature Genetics, 29(4), 365-371. doi:10.1038/ng1201-365Taylor, C. F., Paton, N. W., Lilley, K. S., Binz, P.-A., Julian, R. K., Jones, A. R., … Hermjakob, H. (2007). The minimum information about a proteomics experiment (MIAPE). Nature Biotechnology, 25(8), 887-893. doi:10.1038/nbt1329Yang, Y., Li, R., & Qi, M. (2000). In vivo analysis of plant promoters and transcription factors by agroinfiltration of tobacco leaves. The Plant Journal, 22(6), 543-551. doi:10.1046/j.1365-313x.2000.00760.xKapila, J., De Rycke, R., Van Montagu, M., & Angenon, G. (1997). An Agrobacterium-mediated transient gene expression system for intact leaves. Plant Science, 122(1), 101-108. doi:10.1016/s0168-9452(96)04541-4Schaumberg, K. A., Antunes, M. S., Kassaw, T. K., Xu, W., Zalewski, C. S., Medford, J. I., & Prasad, A. (2015). Quantitative characterization of genetic parts and circuits for plant synthetic biology. Nature Methods, 13(1), 94-100. doi:10.1038/nmeth.3659Dafny-Yelin, M., & Tzfira, T. (2007). Delivery of Multiple Transgenes to Plant Cells. Plant Physiology, 145(4), 1118-1128. doi:10.1104/pp.107.106104TZFIRA, T. (2004). Agrobacterium T-DNA integration: molecules and models. Trends in Genetics, 20(8), 375-383. doi:10.1016/j.tig.2004.06.004De Buck, S., De Wilde, C., Van Montagu, M., & Depicker, A. (2000). Determination of the T-DNA Transfer and the T-DNA Integration Frequencies upon Cocultivation ofArabidopsis thalianaRoot Explants. Molecular Plant-Microbe Interactions, 13(6), 658-665. doi:10.1094/mpmi.2000.13.6.658Konermann, S., Brigham, M. D., Trevino, A. E., Joung, J., Abudayyeh, O. O., Barcena, C., … Zhang, F. (2014). Genome-scale transcriptional activation by an engineered CRISPR-Cas9 complex. Nature, 517(7536), 583-588. doi:10.1038/nature14136Gilbert, L. A., Horlbeck, M. A., Adamson, B., Villalta, J. E., Chen, Y., Whitehead, E. H., … Weissman, J. S. (2014). Genome-Scale CRISPR-Mediated Control of Gene Repression and Activation. Cell, 159(3), 647-661. doi:10.1016/j.cell.2014.09.029Moore, I., Samalova, M., & Kurup, S. (2006). Transactivated and chemically inducible gene expression in plants. The Plant Journal, 45(4), 651-683. doi:10.1111/j.1365-313x.2006.02660.xButelli, E., Titta, L., Giorgio, M., Mock, H.-P., Matros, A., Peterek, S., … Martin, C. (2008). Enrichment of tomato fruit with health-promoting anthocyanins by expression of select transcription factors. Nature Biotechnology, 26(11), 1301-1308. doi:10.1038/nbt.1506Müller, K., Siegel, D., Rodriguez Jahnke, F., Gerrer, K., Wend, S., Decker, E. L., … Zurbriggen, M. D. (2014). A red light-controlled synthetic gene expression switch for plant systems. Mol. BioSyst., 10(7), 1679-1688. doi:10.1039/c3mb70579jDavuluri, G. R., van Tuinen, A., Fraser, P. D., Manfredonia, A., Newman, R., Burgess, D., … Bowler, C. (2005). Fruit-specific RNAi-mediated suppression of DET1 enhances carotenoid and flavonoid content in tomatoes. Nature Biotechnology, 23(7), 890-895. doi:10.1038/nbt1108Zhang, Y., Butelli, E., Alseekh, S., Tohge, T., Rallapalli, G., Luo, J., … Martin, C. (2015). Multi-level engineering facilitates the production of phenylpropanoid compounds in tomato. Nature Communications, 6(1). doi:10.1038/ncomms9635Ori, N., Juarez, M. T., Jackson, D., Yamaguchi, J., Banowetz, G. M., & Hake, S. (1999). Leaf Senescence Is Delayed in Tobacco Plants Expressing the Maize Homeobox Gene knotted1 under the Control of a Senescence-Activated Promoter. The Plant Cell, 11(6), 1073-1080. doi:10.1105/tpc.11.6.1073Azuma, M., Morimoto, R., Hirose, M., Morita, Y., Hoshino, A., Iida, S., … Shiratake, K. (2015). A petal‐specific In MYB 1 promoter from Japanese morning glory: a useful tool for molecular breeding of floricultural crops. Plant Biotechnology Journal, 14(1), 354-363. doi:10.1111/pbi.12389Paine, J. A., Shipton, C. A., Chaggar, S., Howells, R. M., Kennedy, M. J., Vernon, G., … Drake, R. (2005). Improving the nutritional value of Golden Rice through increased pro-vitamin A content. Nature Biotechnology, 23(4), 482-487. doi:10.1038/nbt1082Houmard, N. M., Mainville, J. L., Bonin, C. P., Huang, S., Luethy, M. H., & Malvar, T. M. (2007). High-lysine corn generated by endosperm-specific suppression of lysine catabolism using RNAi. Plant Biotechnology Journal, 5(5), 605-614. doi:10.1111/j.1467-7652.2007.00265.xPadidam, M. (2003). Chemically regulated gene expression in plants. Current Opinion in Plant Biology, 6(2), 169-177. doi:10.1016/s1369-5266(03)00005-0Kinkema, M., Geijskes, R. J., Shand, K., Coleman, H. D., De Lucca, P. C., Palupe, A., … Sainz, M. B. (2013). An improved chemically inducible gene switch that functions in the monocotyledonous plant sugar cane. Plant Molecular Biology, 84(4-5), 443-454. doi:10.1007/s11103-013-0140-

ngs_backbone: a pipeline for read cleaning, mapping and SNP calling using Next Generation Sequence

Background: The possibilities offered by next generation sequencing (NGS) platforms are revolutionizing biotechnological laboratories. Moreover, the combination of NGS sequencing and affordable high-throughput genotyping technologies is facilitating the rapid discovery and use of SNPs in non-model species. However, this abundance of sequences and polymorphisms creates new software needs. To fulfill these needs, we have developed a powerful, yet easy-to-use application. Results: The ngs_backbone software is a parallel pipeline capable of analyzing Sanger, 454, Illumina and SOLiD (Sequencing by Oligonucleotide Ligation and Detection) sequence reads. Its main supported analyses are: read cleaning, transcriptome assembly and annotation, read mapping and single nucleotide polymorphism (SNP) calling and selection. In order to build a truly useful tool, the software development was paired with a laboratory experiment. All public tomato Sanger EST reads plus 14.2 million Illumina reads were employed to test the tool and predict polymorphism in tomato. The cleaned reads were mapped to the SGN tomato transcriptome obtaining a coverage of 4.2 for Sanger and 8.5 for Illumina. 23,360 single nucleotide variations (SNVs) were predicted. A total of 76 SNVs were experimentally validated, and 85% were found to be real. Conclusions: ngs_backbone is a new software package capable of analyzing sequences produced by NGS technologies and predicting SNVs with great accuracy. In our tomato example, we created a highly polymorphic collection of SNVs that will be a useful resource for tomato researchers and breeders. The software developed along with its documentation is freely available under the AGPL license and can be downloaded from http://bioinf. comav.upv.es/ngs_backbone/ or http://github.com/JoseBlanca/franklin.Blanca Postigo, JM.; Pascual Bañuls, L.; Ziarsolo Areitioaurtena, P.; Nuez Viñals, F.; Cañizares Sales, J. (2011). Ngs_backbone: a pipeline for read cleaning, mapping and SNP calling using Next Generation Sequence. BMC Genomics. 12:1-8. doi:10.1186/1471-2164-12-285S1812Metzker ML: Sequencing technologies - the next generation. Nature Reviews Genetics. 2010, 11 (1): 31-46. 10.1038/nrg2626.454 sequencing. [ http://www.454.com/ ]Illumina Inc. [ http://www.illumina.com/ ]Flicek P, Birney E: Sense from sequence reads: methods for alignment and assembly (vol 6, pg S6, 2009). Nature Methods. 2010, 7 (6): 479-479.Chevreux B, Pfisterer T, Drescher B, Driesel AJ, Muller WEG, Wetter T, Suhai S: Using the miraEST assembler for reliable and automated mRNA transcript assembly and SNP detection in sequenced ESTs. Genome Research. 2004, 14 (6): 1147-1159. 10.1101/gr.1917404.Li H, Durbin R: Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009, 25 (14): 1754-1760. 10.1093/bioinformatics/btp324.Langmead B, Trapnell C, Pop M, Salzberg SL: Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biology. 2009, 10 (3):Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R, Genome Project Data P: The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009, 25 (16): 2078-2079. 10.1093/bioinformatics/btp352.1000 Genomes. A deep Catalog of Human Genetic Variation. [ http://1000genomes.org/wiki/doku.php?id=1000_genomes:analysis:vcf4.0 ]The seqanswers internet forum. [ http://seqanswers.com/ ]Blankenberg D, Taylor J, Schenck I, He JB, Zhang Y, Ghent M, Veeraraghavan N, Albert I, Miller W, Makova KD, Ross CH, Nekrutenko A: A framework for collaborative analysis of ENCODE data: Making large-scale analyses biologist-friendly. Genome Research. 2007, 17 (6): 960-964. 10.1101/gr.5578007.CloVR Automated Sequence Analysis from Your Desktop. [ http://clovr.org/ ]Papanicolaou A, Stierli R, Ffrench-Constant RH, Heckel DG: Next generation transcriptomes for next generation genomes using est2assembly. Bmc Bioinformatics. 2009, 10:Applied Biosystems by life technologies. [ http://www.appliedbiosystems.com/absite/us/en/home/applications-technologies/solid-next-generation-sequencing.html ]Wall PK, Leebens-Mack J, Chanderbali AS, Barakat A, Wolcott E, Liang HY, Landherr L, Tomsho LP, Hu Y, Carlson JE, Ma H, Schuster SC, Soltis DE, Soltis PS, Altman N, dePamphilis CW: Comparison of next generation sequencing technologies for transcriptome characterization. Bmc Genomics. 2009, 10:Murchison EP, Tovar C, Hsu A, Bender HS, Kheradpour P, Rebbeck CA, Obendorf D, Conlan C, Bahlo M, Blizzard CA, Pyecroft S, Kreiss A, Kellis M, Stark A, Harkins TT, Marshall Graves JA, Woods GM, Hanon GJ, Papenfuss AT: The Tasmanian Devil Transcriptome Reveals Schwann Cell Origins of a Clonally Transmissible Cancer. Science. 2010, 327 (5961): 84-87. 10.1126/science.1180616.Parchman TL, Geist KS, Grahnen JA, Benkman CW, Buerkle CA: Transcriptome sequencing in an ecologically important tree species: assembly, annotation, and marker discovery. Bmc Genomics. 2010, 11:Babik W, Stuglik M, Qi W, Kuenzli M, Kuduk K, Koteja P, Radwan J: Heart transcriptome of the bank vole (Myodes glareolus): towards understanding the evolutionary variation in metabolic rate. BMC Genomics. 2010, 11: 390-10.1186/1471-2164-11-390.Miller JC, Tanksley SD: RFLP analysis of phylogenetic-relationships and genetic-variation in the genus Lycopersicon. Theoretical and Applied Genetics. 1990, 80 (4): 437-448.Williams CE, Stclair DA: Phenetic relationships and levels of variability detected by restriction-fragment-length-polymorphism and random amplified polymorphic DNA analysis of cultivated and wild accessions of Lycopersicon-esculentum. Genome. 1993, 36 (3): 619-630. 10.1139/g93-083.Rick CM: Tomato, Lycopersicon esculentum (Solanaceae). Evolution of crop plants. Edited by: Simmonds NW. 1976, London: Longman Group, 268-273.Labate JA, Baldo AM: Tomato SNP discovery by EST mining and resequencing. Molecular Breeding. 2005, 16 (4): 343-349. 10.1007/s11032-005-1911-5.Yano K, Watanabe M, Yamamoto N, Maeda F, Tsugane T, Shibata D: Expressed sequence tags (EST) database of a miniature tomato cultivar, Micro-Tom. Plant and Cell Physiology. 2005, 46: S139-S139.Jimenez-Gomez JM, Maloof JN: Sequence diversity in three tomato species: SNPs, markers, and molecular evolution. Bmc Plant Biology. 2009, 9:Yang WC, Bai XD, Kabelka E, Eaton C, Kamoun S, van der Knaap E, Francis D: Discovery of single nucleotide polymorphisms in Lycopersicon esculentum by computer aided analysis of expressed sequence tags. Molecular Breeding. 2004, 14 (1): 21-34.Van Deynze A, Stoffel K, Buell CR, Kozik A, Liu J, van der Knaap E, Francis D: Diversity in conserved genes in tomato. Bmc Genomics. 2007, 8:Sim SC, Robbins MD, Chilcott C, Zhu T, Francis DM: Oligonucleotide array discovery of polymorphisms in cultivated tomato (Solanum lycopersicum L.) reveals patterns of SNP variation associated with breeding. Bmc Genomics. 2009, 10:Bioinformatics at COMAV. [ http://bioinf.comav.upv.es/ngs_backbone/index.html ]Broad institute. [ http://www.broadinstitute.org/igv ]Bioinformatics at COMAV. [ http://bioinf.comav.upv.es/ngs_backbone/install.html ]Github social coding. [ http://github.com/JoseBlanca/franklin ]Chou HH, Holmes MH: DNA sequence quality trimming and vector removal. Bioinformatics. 2001, 17 (12): 1093-1104. 10.1093/bioinformatics/17.12.1093.Picard. [ http://picard.sourceforge.net/index.shtml ]McKenna A, Hanna M, Banks E, Sivachenko A, Citulskis K, Kernytsky A, Garimella K, Altshuler D, Gabriel S, Daly M, DePristo MA: The Genome Analysis Toolkit: A MapReduce framework for analyzing next-generation DNA sequencing data. Genome Research. 2010, 20: 1297-1303. 10.1101/gr.107524.110.Sol Genomics Network. [ ftp://ftp.solgenomics.net/ ]NCBI Genbank. [ http://www.ncbi.nlm.nih.gov/genbank/ ]Gundry CN, Vandersteen JG, Reed GH, Pryor RJ, Chen J, Wittwer CT: Amplicon melting analysis with labeled primers: A closed-tube method for differentiating homozygotes and heterozygotes. Clinical Chemistry. 2003, 49 (3): 396-406. 10.1373/49.3.396

A consensus linkage map for molecular markers and Quantitative Trait Loci associated with economically important traits in melon (Cucumis melo L.)

Background A number of molecular marker linkage maps have been developed for melon (Cucumis melo L.) over the last two decades. However, these maps were constructed using different marker sets, thus, making comparative analysis among maps difficult. In order to solve this problem, a consensus genetic map in melon was constructed using primarily highly transferable anchor markers that have broad potential use for mapping, synteny, and comparative quantitative trait loci (QTL) analysis, increasing breeding effectiveness and efficiency via marker-assisted selection (MAS). Results Under the framework of the International Cucurbit Genomics Initiative (ICuGI, http://www.icugi.org webcite), an integrated genetic map has been constructed by merging data from eight independent mapping experiments using a genetically diverse array of parental lines. The consensus map spans 1150 cM across the 12 melon linkage groups and is composed of 1592 markers (640 SSRs, 330 SNPs, 252 AFLPs, 239 RFLPs, 89 RAPDs, 15 IMAs, 16 indels and 11 morphological traits) with a mean marker density of 0.72 cM/marker. One hundred and ninety-six of these markers (157 SSRs, 32 SNPs, 6 indels and 1 RAPD) were newly developed, mapped or provided by industry representatives as released markers, including 27 SNPs and 5 indels from genes involved in the organic acid metabolism and transport, and 58 EST-SSRs. Additionally, 85 of 822 SSR markers contributed by Syngenta Seeds were included in the integrated map. In addition, 370 QTL controlling 62 traits from 18 previously reported mapping experiments using genetically diverse parental genotypes were also integrated into the consensus map. Some QTL associated with economically important traits detected in separate studies mapped to similar genomic positions. For example, independently identified QTL controlling fruit shape were mapped on similar genomic positions, suggesting that such QTL are possibly responsible for the phenotypic variability observed for this trait in a broad array of melon germplasm. Conclusions Even though relatively unsaturated genetic maps in a diverse set of melon market types have been published, the integrated saturated map presented herein should be considered the initial reference map for melon. Most of the mapped markers contained in the reference map are polymorphic in diverse collection of germplasm, and thus are potentially transferrable to a broad array of genetic experimentation (e.g., integration of physical and genetic maps, colinearity analysis, map-based gene cloning, epistasis dissection, and marker-assisted selection).This work was supported in part by SNC Laboratoire ASL, Ruiter Seeds B.V., Enza Zaden B.V., Gautier Semences S.A., Nunhems B.V., Rijk Zwaan B.V., Sakata Seed Inc, Semillas Fito S. A., Seminis Vegetable Seeds Inc, Syngenta Seeds B. V., Takii and Company Ltd, Vilmorin & Cie S. A., and Zeraim Gedera Ltd (all of them as part of the support to the ICuGI); the grants AGL2009-12698-C02-02 from the Spanish "Ministerio de Ciencia e Innovacion" to AJM. NK lab was supported in part by Research Grant Award No. IS-4223-09C from BARD, the United States - Israel Binational Agricultural Research and Development Fund, and in part by Israel Science Foundation Grant No. 38606, De Ruiter Seeds, Enza Zaden, Keygene, Rijk Zwaan, Sakata Seed Corporation, Semillas Fito, Syngenta Seeds and Vilmorin Clause & Cie. AD was supported by a JAE-Doc contract from "Consejo Superior de Investigaciones Cientificas" (CSIC-Spain). MF was supported by a postdoctoral contract from CRAG. The research carried out at YX's laboratory was supported by Chinese funds (Grant No. 2008-Z42(3), 5100001, 2010AA101907).Díaz Bermúdez, A.; Fergany, M.; Formisano, G.; Ziarsolo, P.; Blanca Postigo, JM.; Fei, Z.; Staub, JE.... (2011). A consensus linkage map for molecular markers and Quantitative Trait Loci associated with economically important traits in melon. BMC Plant Biology. 11. https://doi.org/10.1186/1471-2229-11-111S1