Application of genomic tools in plant breeding

Trends in plant research using molecular markers

A deep bibliometric analysis has been carried out, obtaining valuable parameters that facilitate the understanding around the research in plant using molecular markers. The evolution of the improvement in the field of agronomy is fundamental for its adaptation to the new exigencies that the current world context raises. In addition, within these improvements, this article focuses on those related to the biotechnology sector. More specifically, the use of DNA markers that allow the researcher to know the set of genes associated with a particular quantitative trait or QTL. The use of molecular markers is widely extended, including: restriction fragment length polymorphism, random-amplified polymorphic DNA, amplified fragment length polymorphism, microsatellites, and single-nucleotide polymorphisms. In addition to classical methodology, new approaches based on the next generation sequencing are proving to be fundamental. In this article, a historical review of the molecular markers traditionally used in plants, since its birth and how the new molecular tools facilitate the work of plant breeders is carried out. The evolution of the most studied cultures from the point of view of molecular markers is also reviewed and other parameters whose prior knowledge can facilitate the approach of researchers to this field of research are analyzed. The bibliometric analysis of molecular markers in plants shows that top five countries in this research are: US, China, India, France, and Germany, and from 2013, this research is led by China. On the other hand, the basic research using Arabidopsis is deeper in France and Germany, while other countries focused its efforts in their main crops as the US for wheat or maize, while China and India for wheat and rice.

Molecular mechanisms underlying the diverse array of petal colors in chrysanthemum flowers

Chrysanthemum (Chrysanthemum morifolium Ramat.) is one of the most important floricultural crops in the world. Although the origin of modern chrysanthemum cultivars is uncertain, several species belonging to the family Asteraceae are considered to have been integrated during the long history of breeding. The flower color of ancestral species is limited to yellow, pink, and white, and is derived from carotenoids, anthocyanins, and the absence of both pigments, respectively. A wide range of flower colors, including purplish-red, orange, red, and dark red, has been developed by increasing the range of pigment content or the combination of both pigments. Recently, green-flowered cultivars containing chlorophylls in their ray petals have been produced, and have gained popularity. In addition, blue/violet flowers have been developed using a transgenic approach. Flower color is an important trait that influences the commercial value of chrysanthemum cultivars. Understanding the molecular mechanisms that regulate flower pigmentation may provide important implications for the rationale manipulation of flower color. This review describes the pigment composition, genetics, and molecular basis of ray petal color formation in chrysanthemum cultivars.

Parent-progeny imputation from pooled samples for cost-efficient genotyping in plant breeding

The increased usage of whole-genome selection (WGS) and other molecular evaluation methods in plant breeding relies on the ability to genotype a very large number of untested individuals in each breeding cycle. Many plant breeding programs evaluate large biparental populations of homozygous individuals derived from homozygous parent inbred lines. This structure lends itself to parent-progeny imputation, which transfers the genotype scores of the parents to progeny individuals that are genotyped for a much smaller number of loci. Here we introduce a parent-progeny imputation method that infers individual genotypes from non-barcoded pooled samples of DNA of multiple individuals using a Hidden Markov Model (HMM). We demonstrate the method for pools of simulated maize double haploids (DH) from biparental populations, genotyped using a genotyping by sequencing (GBS) approach for 3,000 loci at 0.125x to 4x coverage. We observed high concordance between true and imputed marker scores and the HMM produced well-calibrated genotype probabilities that correctly reflected the uncertainty of the imputed scores. Genomic estimated breeding values (GEBV) calculated from the imputed scores closely matched GEBV calculated from the true marker scores. The within-population correlation between these sets of GEBV approached 0.95 at 1x and 4x coverage when pooling two or four individuals, respectively. Our approach can reduce the genotyping cost per individual by a factor up to the number of pooled individuals in GBS applications without the need for extra sequencing coverage, thereby enabling cost-effective large scale genotyping for applications such as WGS in plant breeding.

GWAS for plant growth stages and yield components in spring wheat (Triticum aestivum L.) harvested in three regions of Kazakhstan

BACKGROUND

Spring wheat is the largest agricultural crop grown in Kazakhstan with an annual sowing area of 12 million hectares in 2016. Annually, the country harvests around 15 million tons of high quality grain. Despite environmental stress factors it is predicted that the use of new technologies may lead to increases in productivity from current levels of 1.5 to up to 3 tons per hectare. One way of improving wheat productivity is by the application of new genomic oriented approaches in plant breeding projects. Genome wide association studies (GWAS) are emerging as powerful tools for the understanding of the inheritance of complex traits via utilization of high throughput genotyping technologies and phenotypic assessments of plant collections. In this study, phenotyping and genotyping data on 194 spring wheat accessions from Kazakhstan, Russia, Europe, and CIMMYT were assessed for the identification of marker-trait associations (MTA) of agronomic traits by using GWAS.

RESULTS

Field trials in Northern, Central and Southern regions of Kazakhstan using 194 spring wheat accessions revealed strong correlations of yield with booting date, plant height, biomass, number of spikes per plant, and number of kernels per spike. The accessions from Europe and CIMMYT showed high breeding potential for Southern and Central regions of the country in comparison with the performance of the local varieties. The GGE biplot method, using average yield per plant, suggested a clear separation of accessions into their three breeding origins in relationship to the three environments in which they were evaluated. The genetic variation in the three groups of accessions was further studied using 3245 polymorphic SNP (single nucleotide polymorphism) markers. The application of Principal Coordinate analysis clearly grouped the 194 accessions into three clades according to their breeding origins. GWAS on data from nine field trials allowed the identification of 114 MTAs for 12 different agronomic traits.

CONCLUSIONS

Field evaluation of foreign germplasm revealed its poor yield performance in Northern Kazakhstan, which is the main wheat growing region in the country. However, it was found that EU and CIMMYT germplasm has high breeding potential to improve yield performance in Central and Southern regions. The use of Principal Coordinate analysis clearly separated the panel into three distinct groups according to their breeding origin. GWAS based on use of the TASSEL 5.0 package allowed the identification of 114 MTAs for twelve agronomic traits. The study identifies a network of key genes for improvement of yield productivity in wheat growing regions of Kazakhstan.

Data on identification of conserved and novel miRNAs in Elettaria cardamomum

Elettaria cardamomum (L.) Maton, or small cardamom referred as 'queen of spices', is a perennial herbaceous rhizomatous monocot of the family Zingiberaceae. Cardamom seeds and fruits are the economically significant parts and effectively used as a traditional medicine, food additive and flavoring agent. In the present study, using Ion Proton next generation sequencing technology we performed the small RNA sequencing, conserved and novel miRNA predictions of a wild and five cultivar genotypes of cardamom. Small RNA sequencing generated a total of 5,451,328 and 2,756,250 raw reads for wild and cultivar cardamom respectively. The raw data was submitted to SRA database of NCBI under the accession numbers and SRX2273863 (wild) and SRX2273862 (cultivars). The raw reads were quality filtered and predicted conserved and novel miRNAs for wild and cultivar cardamom. The predicted miRNAs, miRNA-targets and functional annotations might provide valuable insights into differences between wild progenitor and cultivated cardamom.

Genomic-based-breeding tools for tropical maize improvement

Maize has traditionally been the main staple diet in the Southern Asia and Sub-Saharan Africa and widely grown by millions of resource poor small scale farmers. Approximately, 35.4 million hectares are sown to tropical maize, constituting around 59% of the developing worlds. Tropical maize encounters tremendous challenges besides poor agro-climatic situations with average yields recorded <3 tones/hectare that is far less than the average of developed countries. On the contrary to poor yields, the demand for maize as food, feed, and fuel is continuously increasing in these regions. Heterosis breeding introduced in early 90 s improved maize yields significantly, but genetic gains is still a mirage, particularly for crop growing under marginal environments. Application of molecular markers has accelerated the pace of maize breeding to some extent. The availability of array of sequencing and genotyping technologies offers unrivalled service to improve precision in maize-breeding programs through modern approaches such as genomic selection, genome-wide association studies, bulk segregant analysis-based sequencing approaches, etc. Superior alleles underlying complex traits can easily be identified and introgressed efficiently using these sequence-based approaches. Integration of genomic tools and techniques with advanced genetic resources such as nested association mapping and backcross nested association mapping could certainly address the genetic issues in maize improvement programs in developing countries. Huge diversity in tropical maize and its inherent capacity for doubled haploid technology offers advantage to apply the next generation genomic tools for accelerating production in marginal environments of tropical and subtropical world. Precision in phenotyping is the key for success of any molecular-breeding approach. This article reviews genomic technologies and their application to improve agronomic traits in tropical maize breeding has been reviewed in detail.

Toward Genomics-Based Breeding in C3 Cool-Season Perennial Grasses

Most important food and feed crops in the world belong to the C3 grass family. The future of food security is highly reliant on achieving genetic gains of those grasses. Conventional breeding methods have already reached a plateau for improving major crops. Genomics tools and resources have opened an avenue to explore genome-wide variability and make use of the variation for enhancing genetic gains in breeding programs. Major C3 annual cereal breeding programs are well equipped with genomic tools; however, genomic research of C3 cool-season perennial grasses is lagging behind. In this review, we discuss the currently available genomics tools and approaches useful for C3 cool-season perennial grass breeding. Along with a general review, we emphasize the discussion focusing on forage grasses that were considered orphan and have little or no genetic information available. Transcriptome sequencing and genotype-by-sequencing technology for genome-wide marker detection using next-generation sequencing (NGS) are very promising as genomics tools. Most C3 cool-season perennial grass members have no prior genetic information; thus NGS technology will enhance collinear study with other C3 model grasses like Brachypodium and rice. Transcriptomics data can be used for identification of functional genes and molecular markers, i.e., polymorphism markers and simple sequence repeats (SSRs). Genome-wide association study with NGS-based markers will facilitate marker identification for marker-assisted selection. With limited genetic information, genomic selection holds great promise to breeders for attaining maximum genetic gain of the cool-season C3 perennial grasses. Application of all these tools can ensure better genetic gains, reduce length of selection cycles, and facilitate cultivar development to meet the future demand for food and fodder.

Emerging Avenues for Utilization of Exotic Germplasm

Breeders have been successful in increasing crop performance by exploiting genetic diversity over time. However, the reported annual yield increases are not sufficient in view of rapid human population growth and global environmental changes. Exotic germplasm possesses high levels of genetic diversity for valuable traits. However, only a small fraction of naturally occurring genetic diversity is utilized. Moreover, the yield gap between elite and exotic germplasm widens, which increases the effort needed to use exotic germplasm and to identify beneficial alleles and for their introgression. The advent of high-throughput genotyping and phenotyping technologies together with emerging biotechnologies provide new opportunities to explore exotic genetic variation. This review will summarize potential challenges for utilization of exotic germplasm and provide solutions.

Transcriptome analysis of lentil (Lens culinaris Medikus) in response to seedling drought stress

BACKGROUND

Drought stress is one of the most harmful abiotic stresses in crop plants. As a moderately drought tolerant crop, lentil is a major crop in rainfed areas and a suitable candidate for drought stress tolerance research work. Screening for drought tolerance stress under hydroponic conditions at seedling stage with air exposure is an efficient technique to select genotypes with contrasting traits. Transcriptome analysis provides valuable resources, especially for lentil, as here the information on complete genome sequence is not available. Hence, the present studies were carried out.

RESULTS

This study was undertaken to understand the biochemical mechanisms and transcriptome changes involved in imparting adaptation to drought stress at seedling stage in drought-tolerant (PDL-2) and drought-sensitive (JL-3) cultivars. Among different physiological and biochemical parameters, a significant increase was recorded in proline, glycine betaine contents and activities of SOD, APX and GPX in PDL-2 compared to JL-3while chlorophyll, RWC and catalase activity decreased significantly in JL-3. Transcriptome changes between the PDL-2 and JL-3 under drought stress were evaluated using Illumina HiSeq 2500 platform. Total number of bases ranged from 5.1 to 6.7 Gb. Sequence analysis of control and drought treated cDNA libraries of PDL-2 and JL-3 produced 74032, 75500, 78328 and 81523 contigs, respectively with respective N50 value of 2011, 2008, 2000 and 1991. Differential gene expression of drought treated genotypes along with their controls revealed a total of 11,435 upregulated and 6,934 downregulated transcripts. For functional classification of DEGs, KEGG pathway annotation analysis extracted a total of 413 GO annotation terms where 176 were within molecular process, 128 in cellular and 109 in biological process groups.

CONCLUSION

The transcriptional profiles provide a foundation for deciphering the underlying mechanism for drought tolerance in lentil. Transcriptional regulation, signal transduction and secondary metabolism in two genotypes revealed significant differences at seedling stage under severe drought. Our finding suggests role of candidate genes for improving drought tolerance in lentil.

Rice grain nutritional traits and their enhancement using relevant genes and QTLs through advanced approaches

BACKGROUND

Rice breeding program needs to focus on development of nutrient dense rice for value addition and helping in reducing malnutrition. Mineral and vitamin deficiency related problems are common in the majority of the population and more specific to developing countries as their staple food is rice.

RESULTS

Genes and QTLs are recently known for the nutritional quality of rice. By comprehensive literature survey and public domain database, we provided a critical review on nutritional aspects like grain protein and amino acid content, vitamins and minerals, glycemic index value, phenolic and flavonoid compounds, phytic acid, zinc and iron content along with QTLs linked to these traits. In addition, achievements through transgenic and advanced genomic approaches have been discussed. The information available on genes and/or QTLs involved in enhancement of micronutrient element and amino acids are summarized with graphical representation.

CONCLUSION

Compatible QTLs/genes may be combined together to design a desirable genotype with superior in multiple grain quality traits. The comprehensive review will be helpful to develop nutrient dense rice cultivars by integrating molecular markers and transgenic assisted breeding approaches with classical breeding.

Reverse Genetics and High Throughput Sequencing Methodologies for Plant Functional Genomics

In the post-genomic era, increasingly sophisticated genetic tools are being developed with the long-term goal of understanding how the coordinated activity of genes gives rise to a complex organism. With the advent of the next generation sequencing associated with effective computational approaches, wide variety of plant species have been fully sequenced giving a wealth of data sequence information on structure and organization of plant genomes. Since thousands of gene sequences are already known, recently developed functional genomics approaches provide powerful tools to analyze plant gene functions through various gene manipulation technologies. Integration of different omics platforms along with gene annotation and computational analysis may elucidate a complete view in a system biology level. Extensive investigations on reverse genetics methodologies were deployed for assigning biological function to a specific gene or gene product. We provide here an updated overview of these high throughout strategies highlighting recent advances in the knowledge of functional genomics in plants.

Target Amplicon Sequencing for Genotyping Genome-Wide Single Nucleotide Polymorphisms Identified by Whole-Genome Resequencing in Peanut

Genome-wide genotyping data regarding breeding materials are essential resources for improving breeding efficiency, especially in plants with complex genomes with a high degree of polyploidy. Several current breeding efforts in cultivated peanut ( L.), which has a tetraploid genome, are devoted to developing high oleic acid cultivars. Genetic maps for such breeding programs have been developed using simple-sequence repeat (SSR) markers, the use of which requires time-consuming electrophoretic analyses. Next-generation sequencing (NGS) technology can overcome this technical hurdle. Initially, we attempted double-digest restriction site-associated DNA sequencing on peanut breeding materials used to develop high oleic acid cultivars. However, this method was not effective because few single nucleotide polymorphism (SNPs) were available because of low genetic diversity of the lines. The genome sequences of the probable diploid ancestors of cultivated peanut, Krapov. & W. C. Greg. and Krapov. & W. C. Greg., are available. Therefore, we next employed whole-genome resequencing analysis to obtain genome-wide SNP data. In this analysis, we observed large biases in the numbers and genomic positions of interspecific and intraspecific SNPs. For genome-wide genotyping, we selected a subset of SNPs covering the peanut genome as the targets of amplicon sequencing analysis. Using this technique, genome-wide genotypes of the breeding materials were easily and rapidly determined. The SNP information and analytic methods developed in this study should accelerate genetics, genomics, and breeding in peanut.

Genomics Approaches For Improving Salinity Stress Tolerance in Crop Plants

Salinity is one of the major factors which reduces crop production worldwide. Plant responses to salinity are highly complex and involve a plethora of genes. Due to its multigenicity, it has been difficult to attain a complete understanding of how plants respond to salinity. Genomics has progressed tremendously over the past decade and has played a crucial role towards providing necessary knowledge for crop improvement. Through genomics, we have been able to identify and characterize the genes involved in salinity stress response, map out signaling pathways and ultimately utilize this information for improving the salinity tolerance of existing crops. The use of new tools, such as gene pyramiding, in genetic engineering and marker assisted breeding has tremendously enhanced our ability to generate stress tolerant crops. Genome editing technologies such as Zinc finger nucleases, TALENs and CRISPR/Cas9 also provide newer and faster avenues for plant biologists to generate precisely engineered crops.

Identification of SNP and SSR Markers in Finger Millet Using Next Generation Sequencing Technologies

Finger millet is an important cereal crop in eastern Africa and southern India with excellent grain storage quality and unique ability to thrive in extreme environmental conditions. Since negligible attention has been paid to improving this crop to date, the current study used Next Generation Sequencing (NGS) technologies to develop both Simple Sequence Repeat (SSR) and Single Nucleotide Polymorphism (SNP) markers. Genomic DNA from cultivated finger millet genotypes KNE755 and KNE796 was sequenced using both Roche 454 and Illumina technologies. Non-organelle sequencing reads were assembled into 207 Mbp representing approximately 13% of the finger millet genome. We identified 10,327 SSRs and 23,285 non-homeologous SNPs and tested 101 of each for polymorphism across a diverse set of wild and cultivated finger millet germplasm. For the 49 polymorphic SSRs, the mean polymorphism information content (PIC) was 0.42, ranging from 0.16 to 0.77. We also validated 92 SNP markers, 80 of which were polymorphic with a mean PIC of 0.29 across 30 wild and 59 cultivated accessions. Seventy-six of the 80 SNPs were polymorphic across 30 wild germplasm with a mean PIC of 0.30 while only 22 of the SNP markers showed polymorphism among the 59 cultivated accessions with an average PIC value of 0.15. Genetic diversity analysis using the polymorphic SNP markers revealed two major clusters; one of wild and another of cultivated accessions. Detailed STRUCTURE analysis confirmed this grouping pattern and further revealed 2 sub-populations within wild E. coracana subsp. africana. Both STRUCTURE and genetic diversity analysis assisted with the correct identification of the new germplasm collections. These polymorphic SSR and SNP markers are a significant addition to the existing 82 published SSRs, especially with regard to the previously reported low polymorphism levels in finger millet. Our results also reveal an unexploited finger millet genetic resource that can be included in the regional breeding programs in order to efficiently optimize productivity.

Molecular Approaches to Genetically Improve the Accumulation of Health-Promoting Secondary Metabolites in Staple Crops-A Case Study: The Lipoxygenase-B1 Genes and Regulation of the Carotenoid Content in Pasta Products

Secondary metabolites, also known as phytochemicals, represent a large subset of plant molecules that include compounds with health-promoting effects. Indeed, a number of epidemiological studies have shown that, when taken regularly and in adequate amounts, these molecules can have long-term beneficial effects on human health, through reduction of the incidence of degenerative diseases, such as cardiovascular diseases, obesity, diabetes, and cancer. As the dietary intake of these phytochemicals is often inadequate, various strategies are in use to improve their content in staple crops, and the end-products thereof. One of the most effective strategies is crop improvement through genetic approaches, as this is the only way to generate new cultivars in which the high accumulation of a given phytochemical is stably fixed. Efforts to genetically improve quality traits are rapidly evolving, from classical breeding to molecular-assisted approaches; these require sound understanding of the molecular bases underlying the traits, to identify the genes/alleles that control them. This can be achieved through global analysis of the metabolic pathway responsible for phytochemical accumulation, to identify the link between phytochemical content and the activities of key enzymes that regulate the metabolic pathway, and between the key enzymes and their encoding genes/alleles. Once these have been identified, they can be used as markers for selection of new improved genotypes through biotechnological approaches. This review provides an overview of the major health-promoting properties shown to be associated with the dietary intake of phytochemicals, and describes how molecular approaches provide means for improving the health quality of edible crops. Finally, a case study is illustrated, of the identification in durum wheat of the Lipoxygenase-B1 genes that control the final carotenoid content in semolina-based foods, such as pasta products.

Genome engineering and plant breeding: impact on trait discovery and development

KEY MESSAGE

New tools for the precise modification of crops genes are now available for the engineering of new ideotypes. A future challenge in this emerging field of genome engineering is to develop efficient methods for allele mining. Genome engineering tools are now available in plants, including major crops, to modify in a predictable manner a given gene. These new techniques have a tremendous potential for a spectacular acceleration of the plant breeding process. Here, we discuss how genetic diversity has always been the raw material for breeders and how they have always taken advantage of the best available science to use, and when possible, increase, this genetic diversity. We will present why the advent of these new techniques gives to the breeders extremely powerful tools for crop breeding, but also why this will require the breeders and researchers to characterize the genes underlying this genetic diversity more precisely. Tackling these challenges should permit the engineering of optimized alleles assortments in an unprecedented and controlled way.

The first de novo transcriptome of pepino (Solanum muricatum): assembly, comprehensive analysis and comparison with the closely related species S. caripense, potato and tomato

BACKGROUND

Solanum sect. Basarthrum is phylogenetically very close to potatoes (Solanum sect. Petota) and tomatoes (Solanum sect. Lycopersicon), two groups with great economic importance, and for which Solanum sect. Basarthrum represents a tertiary gene pool for breeding. This section includes the important regional cultigen, the pepino (Solanum muricatum), and several wild species. Among the wild species, S. caripense is prominent due to its major involvement in the origin of pepino and its wide geographical distribution. Despite the value of the pepino as an emerging crop, and the potential for gene transfer from both the pepino and S. caripense to potatoes and tomatoes, there has been virtually no genomic study of these species.

RESULTS

Using Illumina HiSeq 2000, RNA-Seq was performed with a pool of three tissues (young leaf, flowers in pre-anthesis and mature fruits) from S. muricatum and S. caripense, generating almost 111,000,000 reads among the two species. A high quality de novo transcriptome was assembled from S. muricatum clean reads resulting in 75,832 unigenes with an average length of 704 bp. These unigenes were functionally annotated based on similarity of public databases. We used Blast2GO, to conduct an exhaustive study of the gene ontology, including GO terms, EC numbers and KEGG pathways. Pepino unigenes were compared to both potato and tomato genomes in order to determine their estimated relative position, and to infer gene prediction models. Candidate genes related to traits of interest in other Solanaceae were evaluated by presence or absence and compared with S. caripense transcripts. In addition, by studying five genes, the phylogeny of pepino and five other members of the family, Solanaceae, were studied. The comparison of S. caripense reads against S. muricatum assembled transcripts resulted in thousands of intra- and interspecific nucleotide-level variants. In addition, more than 1000 SSRs were identified in the pepino transcriptome.

CONCLUSIONS

This study represents the first genomic resource for the pepino. We suggest that the data will be useful not only for improvement of the pepino, but also for potato and tomato breeding and gene transfer. The high quality of the transcriptome presented here also facilitates comparative studies in the genus Solanum. The accurate transcript annotation will enable us to figure out the gene function of particular traits of interest. The high number of markers (SSR and nucleotide-level variants) obtained will be useful for breeding programs, as well as studies of synteny, diversity evolution, and phylogeny.

Genomic Bayesian Prediction Model for Count Data with Genotype × Environment Interaction

Genomic tools allow the study of the whole genome, and facilitate the study of genotype-environment combinations and their relationship with phenotype. However, most genomic prediction models developed so far are appropriate for Gaussian phenotypes. For this reason, appropriate genomic prediction models are needed for count data, since the conventional regression models used on count data with a large sample size (nT) and a small number of parameters (p) cannot be used for genomic-enabled prediction where the number of parameters (p) is larger than the sample size (nT). Here, we propose a Bayesian mixed-negative binomial (BMNB) genomic regression model for counts that takes into account genotype by environment (G×E) interaction. We also provide all the full conditional distributions to implement a Gibbs sampler. We evaluated the proposed model using a simulated data set, and a real wheat data set from the International Maize and Wheat Improvement Center (CIMMYT) and collaborators. Results indicate that our BMNB model provides a viable option for analyzing count data.

Sequencing Crop Genomes: A Gateway to Improve Tropical Agriculture

Agricultural development in the tropics lags behind development in the temperate latitudes due to the lack of advanced technology, and various biotic and abiotic factors. To cope with the increasing demand for food and other plant-based products, improved crop varieties have to be developed. To breed improved varieties, a better understanding of crop genetics is necessary. With the advent of next-generation DNA sequencing technologies, many important crop genomes have been sequenced. Primary importance has been given to food crops, including cereals, tuber crops, vegetables, and fruits. The DNA sequence information is extremely valuable for identifying key genes controlling important agronomic traits and for identifying genetic variability among the cultivars. However, massive DNA re-sequencing and gene expression studies have to be performed to substantially improve our understanding of crop genetics. Application of the knowledge obtained from the genomes, transcriptomes, expression studies, and epigenetic studies would enable the development of improved varieties and may lead to a second green revolution. The applications of next generation DNA sequencing technologies in crop improvement, its limitations, future prospects, and the features of important crop genome projects are reviewed herein.

The genetic aspects of berries: from field to health

Berries are a relevant source of micronutrients and nonessential phytochemicals, such as polyphenol compounds, that play a synergistic and cumulative role in human health promotion. Several systematic analyses showed that berry phenolics are able to detoxify reactive oxygen and nitrogen species, blocking their production, to intervene in the cell cycle, participating in the transduction and expression of genes involved in apoptosis, and to repair oxidative DNA damage. As a consequence, the improvement of the nutritional quality of berries has become a new quality target of breeding and biotechnological strategies, to control or to increase the content of specific health-related compounds in fruits. This work reviews, on the basis of the in vitro and in vivo evidence, the main berries' phytochemical compounds and their possible mechanisms of action on pathways involved in several type of diseases, with particular attention to cancer, inflammation, neurodegeneration, diabetes and cardiovascular diseases.

GBS-SNP-CROP: a reference-optional pipeline for SNP discovery and plant germplasm characterization using variable length, paired-end genotyping-by-sequencing data

BACKGROUND

With its simple library preparation and robust approach to genome reduction, genotyping-by-sequencing (GBS) is a flexible and cost-effective strategy for SNP discovery and genotyping, provided an appropriate reference genome is available. For resource-limited curation, research, and breeding programs of underutilized plant genetic resources, however, even low-depth references may not be within reach, despite declining sequencing costs. Such programs would find value in an open-source bioinformatics pipeline that can maximize GBS data usage and perform high-density SNP genotyping in the absence of a reference.

RESULTS

The GBS SNP-Calling Reference Optional Pipeline (GBS-SNP-CROP) developed and presented here adopts a clustering strategy to build a population-tailored "Mock Reference" from the same GBS data used for downstream SNP calling and genotyping. Designed for libraries of paired-end (PE) reads, GBS-SNP-CROP maximizes data usage by eliminating unnecessary data culling due to imposed read-length uniformity requirements. Using 150 bp PE reads from a GBS library of 48 accessions of tetraploid kiwiberry (Actinidia arguta), GBS-SNP-CROP yielded on average three times as many SNPs as TASSEL-GBS analyses (32 and 64 bp tag lengths) and over 18 times as many as TASSEL-UNEAK, with fewer genotyping errors in all cases, as evidenced by comparing the genotypic characterizations of biological replicates. Using the published reference genome of a related diploid species (A. chinensis), the reference-based version of GBS-SNP-CROP behaved similarly to TASSEL-GBS in terms of the number of SNPs called but had an improved read depth distribution and fewer genotyping errors. Our results also indicate that the sets of SNPs detected by the different pipelines above are largely orthogonal to one another; thus GBS-SNP-CROP may be used to augment the results of alternative analyses, whether or not a reference is available.

CONCLUSIONS

By achieving high-density SNP genotyping in populations for which no reference genome is available, GBS-SNP-CROP is worth consideration by curators, researchers, and breeders of under-researched plant genetic resources. In cases where a reference is available, especially if from a related species or when the target population is particularly diverse, GBS-SNP-CROP may complement other reference-based pipelines by extracting more information per sequencing dollar spent. The current version of GBS-SNP-CROP is available at https://github.com/halelab/GBS-SNP-CROP.git.

Next generation sequencing and omics in cucumber (Cucumis sativus L.) breeding directed research

In the post-genomic era the availability of genomic tools and resources is leading us to novel generation methods in plant breeding, as they facilitate the study of the genotype and its relationship with the phenotype, in particular for complex traits. In this study we have mainly concentrated on the Cucumis sativus and (but much less) Cucurbitaceae family several important vegetable crops. There are many reports on research conducted in Cucurbitaceae plant breeding programs on the ripening process, phloem transport, disease resistance, cold tolerance and fruit quality traits. This paper presents the role played by new omic technologies in the creation of knowledge on the mechanisms of the formation of the breeding features. The analysis of NGS (NGS-next generation sequencing) data allows the discovery of new genes and regulatory sequences, their positions, and makes available large collections of molecular markers. Genome-wide expression studies provide breeders with an understanding of the molecular basis of complex traits. Firstly a high density map should be created for the reference genome, then each re-sequencing data could be mapped and new markers brought out into breeding populations. The paper also presents methods that could be used in the future for the creation of variability and genomic modification of the species in question. It has been shown also the state and usefulness in breeding the chloroplastomic and mitochondriomic study.

Molecular Breeding Strategy and Challenges Towards Improvement of Blast Disease Resistance in Rice Crop

Rice is a staple and most important security food crop consumed by almost half of the world's population. More rice production is needed due to the rapid population growth in the world. Rice blast caused by the fungus, Magnaporthe oryzae is one of the most destructive diseases of this crop in different part of the world. Breakdown of blast resistance is the major cause of yield instability in several rice growing areas. There is a need to develop strategies providing long-lasting disease resistance against a broad spectrum of pathogens, giving protection for a long time over a broad geographic area, promising for sustainable rice production in the future. So far, molecular breeding approaches involving DNA markers, such as QTL mapping, marker-aided selection, gene pyramiding, allele mining and genetic transformation have been used to develop new resistant rice cultivars. Such techniques now are used as a low-cost, high-throughput alternative to conventional methods allowing rapid introgression of disease resistance genes into susceptible varieties as well as the incorporation of multiple genes into individual lines for more durable blast resistance. The paper briefly reviewed the progress of studies on this aspect to provide the interest information for rice disease resistance breeding. This review includes examples of how advanced molecular method have been used in breeding programs for improving blast resistance. New information and knowledge gained from previous research on the recent strategy and challenges towards improvement of blast disease such as pyramiding disease resistance gene for creating new rice varieties with high resistance against multiple diseases will undoubtedly provide new insights into the rice disease control.

Can functional hologenomics aid tackling current challenges in plant breeding?

Molecular plant breeding usually overlooks the genetic variability that arises from the association of plants with endophytic microorganisms, when looking at agronomic interesting target traits. This source of variability can have crucial effects on the functionality of the organism considered as a whole (the holobiont), and therefore can be selectable in breeding programs. However, seeing the holobiont as a unit for selection and improvement in breeding programs requires novel approaches for genotyping and phenotyping. These should not focus just at the plant level, but also include the associated endophytes and their functional effects on the plant, to make effective desirable trait screenings. The present review intends to draw attention to a new research field on functional hologenomics that if associated with adequate phenotyping tools could greatly increase the efficiency of breeding programs.

Utilizing next-generation sequencing to study homeologous polymorphisms and herbicide-resistance-endowing mutations in Poa annua acetolactate synthase genes

BACKGROUND

Detection of single nucleotide polymorphisms (SNPs) related to herbicide resistance in non-model polyploid weed species is fraught with difficulty owing to the gene duplication and lack of reference sequences. Our research seeks to overcome these obstacles by Illumina HiSeq read mapping, SNP calling and allele frequency determinations. Our focus is on the acetolactate synthase (ALS) gene, the target site of ALS-inhibiting herbicides, in Poa annua, an allotetraploid weed species originating from two diploid parents, P. supina and P. infirma.

RESULTS

ALS contigs with complete coding regions of P. supina, P. infirma and P. annua were assembled and compared with ALS genes from other plant species. The ALS infirma-homeolog of P. annua showed higher levels of nucleotide sequence variability than the supina-homeolog. Comparisons of read mappings of P. annua and a simulated P. supina × P. infirma hybrid showed high resemblance. Two homeolog-specific primer pairs were designed and used to amplify a 1860 bp region covering all resistance-conferring codons in the ALS gene. Four P. annua populations, GN, RB, GW and LG, showed high resistance to two ALS inhibitors, bispyribac-sodium and foramsulfuron, and two populations, HD and RS, showed lower resistance in the rate-response trial. Mutations conferring Trp-574-Leu substitution were observed in the infirma-homeolog of GN and RB and in the supina-homeolog of GW and LG, but no resistance-conferring mutation was observed in the two populations of lower resistance, HD and RS.

CONCLUSION

In this study we have demonstrated the use of NGS data to study homeologous polymorphisms, parentage and herbicide resistance in an allotetraploid weed species, P. annua. Complete coding sequences of the ALS gene were assembled for P. infirma, P. supina, infirma-homeolog and supina-homeolog in P. annua. A pipeline consisting of read mapping, SNP calling and allele frequency calculation was developed to study the parentage of P. annua, which provided a new perspective on this topic besides the views of morphology, karyotype and phylogeny. Our two homeolog-specific primer pairs can be utilized in future research to separate the homeologs of the ALS gene in P. annua and cover all the codons that have been reported to confer herbicide resistance.

Transcriptional Profiling and Molecular Characterization of Astragalosides, Calycosin, and Calycosin-7-O-β-D-glucoside Biosynthesis in the Hairy Roots of Astragalus membranaceus in Response to Methyl Jasmonate

We used the next-generation Illumina/Solexa HiSeq2000 platform on RNA analysis to investigate the transcriptome of Astragalus membranaceus hairy roots in response to 100 μM methyl jasmonate (MeJA). In total, 77,758,230 clean reads were assembled into 48,636 transcripts (average length of 1398 bp), which were clustered into 23,658 loci (genes). Of these, 19,940 genes were annotated by BLASTx searches. In addition, DESeq analysis showed that 2127 genes were up-regulated, while 1247 genes were down-regulated by MeJA. Seventeen novel astragaloside (AST) biosynthetic genes and seven novel calycosin and calycosin-7-O-β-D-glucoside (CG) biosynthetic genes were isolated. The accumulation of ASTs, calycosin, and CG increased significantly in MeJA-treated hairy roots compared with control hairy roots. Our findings will provide a valuable resource for molecular characterization of AST, calycosin, and CG biosynthetic pathways and may lead to new approaches to maximize their production and biomass productivity in the hairy roots of A. membranaceus.

RNA-Seq analysis identifies genes associated with differential reproductive success under drought-stress in accessions of wild barley Hordeum spontaneum

BACKGROUND

The evolutionary basis of reproductive success in different environments is of major interest in the study of plant adaptation. Since the reproductive stage is particularly sensitive to drought, genes affecting reproductive success during this stage are key players in the evolution of adaptive mechanisms. We used an ecological genomics approach to investigate the reproductive response of drought-tolerant and sensitive wild barley accessions originating from different habitats in the Levant.

RESULTS

We sequenced mRNA extracted from spikelets at the flowering stage in drought-treated and control plants. The barley genome was used for a reference-guided assembly and differential expression analysis. Our approach enabled to detect biological processes affecting grain production under drought stress. We detected novel candidate genes and differentially expressed alleles associated with drought tolerance. Drought associated genes were shown to be more conserved than non-associated genes, and drought-tolerance genes were found to evolve more rapidly than other drought associated genes.

CONCLUSIONS

We show that reproductive success under drought stress is not a habitat-specific trait but a shared physiological adaptation that appeared to evolve recently in the evolutionary history of wild barley. Exploring the genomic basis of reproductive success under stress in crop wild progenitors is expected to have considerable ecological and economical applications.

Sequencing consolidates molecular markers with plant breeding practice

Plenty of molecular markers have been developed by contemporary sequencing technologies, whereas few of them are successfully applied in breeding, thus we present a review on how sequencing can facilitate marker-assisted selection in plant breeding. The growing global population and shrinking arable land area require efficient plant breeding. Novel strategies assisted by certain markers have proven effective for genetic gains. Fortunately, cutting-edge sequencing technologies bring us a deluge of genomes and genetic variations, enlightening the potential of marker development. However, a large gap still exists between the potential of molecular markers and actual plant breeding practices. In this review, we discuss marker-assisted breeding from a historical perspective, describe the road from crop sequencing to breeding, and highlight how sequencing facilitates the application of markers in breeding practice.

Application of in silico bulked segregant analysis for rapid development of markers linked to Bean common mosaic virus resistance in common bean

BACKGROUND

Common bean was one of the first crops that benefited from the development and utilization of molecular marker-assisted selection (MAS) for major disease resistance genes. Efficiency of MAS for breeding common bean is still hampered, however, due to the dominance, linkage phase, and loose linkage of previously developed markers. Here we applied in silico bulked segregant analysis (BSA) to the BeanCAP diversity panel, composed of over 500 lines and genotyped with the BARCBEAN_3 6K SNP BeadChip, to develop codominant and tightly linked markers to the I gene controlling resistance to Bean common mosaic virus (BCMV).

RESULTS

We physically mapped the genomic region underlying the I gene. This locus, in the distal arm of chromosome Pv02, contains seven putative NBS-LRR-type disease resistance genes. Two contrasting bulks, containing BCMV host differentials and ten BeanCAP lines with known disease reaction to BCMV, were subjected to in silico BSA for targeting the I gene and flanking sequences. Two distinct haplotypes, containing a cluster of six single nucleotide polymorphisms (SNP), were associated with resistance or susceptibility to BCMV. One-hundred and twenty-two lines, including 115 of the BeanCAP panel, were screened for BCMV resistance in the greenhouse, and all of the resistant or susceptible plants displayed distinct SNP haplotypes as those found in the two bulks. The resistant/susceptible haplotypes were validated in 98 recombinant inbred lines segregating for BCMV resistance. The closest SNP (~25-32 kb) to the distal NBS-LRR gene model for the I gene locus was targeted for conversion to codominant KASP (Kompetitive Allele Specific PCR) and CAPS (Cleaved Amplified Polymorphic Sequence) markers. Both marker systems accurately predicted the disease reaction to BCMV conferred by the I gene in all screened lines of this study.

CONCLUSIONS

We demonstrated the utility of the in silico BSA approach using genetically diverse germplasm, genotyped with a high-density SNP chip array, to discover SNP variation at a specific targeted genomic region. In common bean, many disease resistance genes are mapped and their physical genomic position can now be determined, thus the application of this approach will facilitate further development of codominant and tightly linked markers for use in MAS.

Application of in silico bulked segregant analysis for rapid development of markers linked to Bean common mosaic virus resistance in common bean

Background

Common bean was one of the first crops that benefited from the development and utilization of molecular marker-assisted selection (MAS) for major disease resistance genes. Efficiency of MAS for breeding common bean is still hampered, however, due to the dominance, linkage phase, and loose linkage of previously developed markers. Here we applied in silico bulked segregant analysis (BSA) to the BeanCAP diversity panel, composed of over 500 lines and genotyped with the BARCBEAN_3 6K SNP BeadChip, to develop codominant and tightly linked markers to the I gene controlling resistance to Bean common mosaic virus (BCMV).

Results

We physically mapped the genomic region underlying the I gene. This locus, in the distal arm of chromosome Pv02, contains seven putative NBS-LRR-type disease resistance genes. Two contrasting bulks, containing BCMV host differentials and ten BeanCAP lines with known disease reaction to BCMV, were subjected to in silico BSA for targeting the I gene and flanking sequences. Two distinct haplotypes, containing a cluster of six single nucleotide polymorphisms (SNP), were associated with resistance or susceptibility to BCMV. One-hundred and twenty-two lines, including 115 of the BeanCAP panel, were screened for BCMV resistance in the greenhouse, and all of the resistant or susceptible plants displayed distinct SNP haplotypes as those found in the two bulks. The resistant/susceptible haplotypes were validated in 98 recombinant inbred lines segregating for BCMV resistance. The closest SNP (~25-32 kb) to the distal NBS-LRR gene model for the I gene locus was targeted for conversion to codominant KASP (Kompetitive Allele Specific PCR) and CAPS (Cleaved Amplified Polymorphic Sequence) markers. Both marker systems accurately predicted the disease reaction to BCMV conferred by the I gene in all screened lines of this study.

Conclusions

We demonstrated the utility of the in silico BSA approach using genetically diverse germplasm, genotyped with a high-density SNP chip array, to discover SNP variation at a specific targeted genomic region. In common bean, many disease resistance genes are mapped and their physical genomic position can now be determined, thus the application of this approach will facilitate further development of codominant and tightly linked markers for use in MAS.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-903) contains supplementary material, which is available to authorized users.

Reducing capacity, chlorogenic acid content and biological activity in a collection of scarlet (Solanum aethiopicum) and Gboma (S. macrocarpon) eggplants

Scarlet (Solanum aethiopicum) and gboma (S. macrocarpon) eggplants are important vegetables in Sub-Saharan Africa. Few studies have been made on these crops regarding the diversity of phenolic content and their biological activity. We have studied the reducing activity, the chlorogenic acid and other phenolic acid contents in a collection of 56 accessions of scarlet eggplant, including the four cultivated groups (Aculeatum, Gilo, Kumba, Shum) and the weedy intermediate S. aethiopicum-S. anguivi types, as well as in eight accessions of gboma eggplant, including the cultivated S. macrocarpon and its wild ancestor, S. dasyphyllum. A sample of the accessions evaluated in this collection has been tested for inhibition of nitric oxide (NO) using macrophage cell cultures. The results show that there is a great diversity in both crops for reducing activity, chlorogenic acid content and chlorogenic acid peak area (% of total phenolic acids). Heritability (H²) for these traits was intermediate to high in both crops. In all samples, chlorogenic acid was the major phenolic acid and accounted for more than 50% of the chromatogram peak area. Considerable differences were found among and within groups for these traits, but the greatest values for total phenolics and chlorogenic acid content were found in S. dasyphyllum. In most groups, reducing activity was positively correlated (with values of up to 0.904 in the Aculeatum group) with chlorogenic acid content. Inhibition of NO was greatest in samples having a high chlorogenic acid content. The results show that both crops are a relevant source of chlorogenic acid and other phenolic acids. The high diversity found also indicates that there are good prospects for breeding new scarlet and gboma eggplant cultivars with improved content in phenolics and bioactive properties.

Natural allelic diversity, genetic structure and linkage disequilibrium pattern in wild chickpea

Characterization of natural allelic diversity and understanding the genetic structure and linkage disequilibrium (LD) pattern in wild germplasm accessions by large-scale genotyping of informative microsatellite and single nucleotide polymorphism (SNP) markers is requisite to facilitate chickpea genetic improvement. Large-scale validation and high-throughput genotyping of genome-wide physically mapped 478 genic and genomic microsatellite markers and 380 transcription factor gene-derived SNP markers using gel-based assay, fluorescent dye-labelled automated fragment analyser and matrix-assisted laser desorption ionization-time of flight (MALDI-TOF) mass array have been performed. Outcome revealed their high genotyping success rate (97.5%) and existence of a high level of natural allelic diversity among 94 wild and cultivated Cicer accessions. High intra- and inter-specific polymorphic potential and wider molecular diversity (11-94%) along with a broader genetic base (13-78%) specifically in the functional genic regions of wild accessions was assayed by mapped markers. It suggested their utility in monitoring introgression and transferring target trait-specific genomic (gene) regions from wild to cultivated gene pool for the genetic enhancement. Distinct species/gene pool-wise differentiation, admixed domestication pattern, and differential genome-wide recombination and LD estimates/decay observed in a six structured population of wild and cultivated accessions using mapped markers further signifies their usefulness in chickpea genetics, genomics and breeding.

Specific requirement for translation initiation factor 4E or its isoform drives plant host susceptibility to Tobacco etch virus

BACKGROUND

In plants, eIF4E translation initiation factors and their eIFiso4E isoforms are essential susceptibility factors for many RNA viruses, including potyviruses. Mutations altering these factors are a major source of resistance to the viruses. The eIF4E allelic series is associated with specific resistance spectra in crops such as Capsicum annum. Genetic evidence shows that potyviruses have a specific requirement for a given 4E isoform that depends on the host plant. For example, Tobacco etch virus (TEV) uses eIF4E1 to infect Capsicum annuum but uses eIFiso4E to infect Arabidopsis thaliana. Here, we investigated how TEV exploits different translation initiation factor isoforms to infect these two plant species.

RESULTS

A complementation system was set up in Arabidopsis to test the restoration of systemic infection by TEV. Using this system, Arabidopsis susceptibility to TEV was complemented with a susceptible pepper eIF4E1 allele but not with a resistant allele. Therefore, in Arabidopsis, TEV can use the pepper eIF4E1 instead of the endogenous eIFiso4E isoform so is able to switch between translation initiation factor 4E isoform to infect the same host. Moreover, we show that overexpressing the pepper eIF4E1 alleles is sufficient to make Arabidopsis susceptible to an otherwise incompatible TEV strain. Lastly, we show that the resistant eIF4E1 allele is similarly overcome by a resistance-breaking TEV strain as in pepper, confirming that this Arabidopsis TEV-susceptibility complementation system is allele-specific.

CONCLUSION

We report here a complementation system in Arabidopsis that makes it possible to assess the role of pepper pvr2-eIF4E alleles in susceptibility to TEV. Heterologous complementation experiments showed that the idiosyncratic properties of the 4E and iso4E proteins create a major checkpoint for viral infection of different hosts. This system could be used to screen natural or induced eIF4E alleles to find and study alleles of interest for plant breeding.

Genome-wide distribution of genetic diversity and linkage disequilibrium in a mass-selected population of maritime pine

Background

The accessibility of high-throughput genotyping technologies has contributed greatly to the development of genomic resources in non-model organisms. High-density genotyping arrays have only recently been developed for some economically important species such as conifers. The potential for using genomic technologies in association mapping and breeding depends largely on the genome wide patterns of diversity and linkage disequilibrium in current breeding populations. This study aims to deepen our knowledge regarding these issues in maritime pine, the first species used for reforestation in south western Europe.

Results

Using a new map merging algorithm, we first established a 1,712 cM composite linkage map (comprising 1,838 SNP markers in 12 linkage groups) by bringing together three already available genetic maps. Using rigorous statistical testing based on kernel density estimation and resampling we identified cold and hot spots of recombination. In parallel, 186 unrelated trees of a mass-selected population were genotyped using a 12k-SNP array. A total of 2,600 informative SNPs allowed to describe historical recombination, genetic diversity and genetic structure of this recently domesticated breeding pool that forms the basis of much of the current and future breeding of this species. We observe very low levels of population genetic structure and find no evidence that artificial selection has caused a reduction in genetic diversity. By combining these two pieces of information, we provided the map position of 1,671 SNPs corresponding to 1,192 different loci. This made it possible to analyze the spatial pattern of genetic diversity (H_e) and long distance linkage disequilibrium (LD) along the chromosomes. We found no particular pattern in the empirical variogram of H_e across the 12 linkage groups and, as expected for an outcrossing species with large effective population size, we observed an almost complete lack of long distance LD.

Conclusions

These results are a stepping stone for the development of strategies for studies in population genomics, association mapping and genomic prediction in this economical and ecologically important forest tree species.

A reverse genetics approach identifies novel mutants in light responses and anthocyanin metabolism in petunia

Flower color and plant architecture are important commercially valuable features for ornamental petunias (Petunia x hybrida Vilm.). Photoperception and light signaling are the major environmental factors controlling anthocyanin and chlorophyll biosynthesis and shade-avoidance responses in higher plants. The genetic regulators of these processes were investigated in petunia by in silico analyses and the sequence information was used to devise a reverse genetics approach to probe mutant populations. Petunia orthologs of photoreceptor, light-signaling components and anthocyanin metabolism genes were identified and investigated for functional conservation by phylogenetic and protein motif analyses. The expression profiles of photoreceptor gene families and of transcription factors regulating anthocyanin biosynthesis were obtained by bioinformatic tools. Two mutant populations, generated by an alkalyting agent and by gamma irradiation, were screened using a phenotype-independent, sequence-based method by high-throughput PCR-based assay. The strategy allowed the identification of novel mutant alleles for anthocyanin biosynthesis (CHALCONE SYNTHASE) and regulation (PH4), and for light signaling (CONSTANS) genes.

SNP genotyping reveals genetic diversity between cultivated landraces and contemporary varieties of tomato

Background

The tomato (Solanum lycopersium L.) is the most widely grown vegetable in the world. It was domesticated in Latin America and Italy and Spain are considered secondary centers of diversification. This food crop has experienced severe genetic bottlenecks and modern breeding activities have been characterized by trait introgression from wild species and divergence in different market classes.

Results

With the aim to examine patterns of polymorphism, characterize population structure and identify putative loci under positive selection, we genotyped 214 tomato accessions (which include cultivated landraces, commercial varieties and wild relatives) using a custom-made Illumina SNP-panel. Most of the 175 successfully scored SNP loci were found to be polymorphic. Population structure analysis and estimates of genetic differentiation indicated that landraces constitute distinct sub-populations. Furthermore, contemporary varieties could be separated in groups (processing, fresh and cherry) that are consistent with the recent breeding aimed at market-class specialization. In addition, at the 95% confidence level, we identified 30, 34 and 37 loci under positive selection between landraces and each of the groups of commercial variety (cherry, processing and fresh market, respectively). Their number and genomic locations imply the presence of some extended regions with high genetic variation between landraces and contemporary varieties.

Conclusions

Our work provides knowledge concerning the level and distribution of genetic variation within cultivated tomato landraces and increases our understanding of the genetic subdivision of contemporary varieties. The data indicate that adaptation and selection have led to a genomic signature in cultivated landraces and that the subpopulation structure of contemporary varieties is shaped by directed breeding and largely of recent origin. The genomic characterization presented here is an essential step towards a future exploitation of the available tomato genetic resources in research and breeding programs.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-14-835) contains supplementary material, which is available to authorized users.