Different horticultural types of lettuce exhibit tremendous morphological variation. However, the molecular basis for domestication and divergence among the different horticultural types of lettuce remains unknown. Here, we report the RNA sequencing of 240 lettuce accessions sampled from the major horticultural types and wild relatives, generating 1.1 million single-nucleotide polymorphisms (SNPs). Demographic modeling indicates that there was a single domestication event for lettuce. We identify a list of regions as putative selective sweeps that occurred during domestication and divergence, respectively. Genome-wide association studies (GWAS) identify 5311 expression quantitative trait loci (eQTL) regulating the expression of 4105 genes, including nine eQTLs regulating genes associated with flavonoid biosynthesis. GWAS for leaf color detects six candidate loci responsible for the variation of anthocyanins in lettuce leaves. Our study provides a comprehensive understanding of the domestication and the accumulation of anthocyanins in lettuce, and will facilitate the breeding of cultivars with improved nutritional value.

BackgroundMolecular breeding of pepper (Capsicum spp.) can be accelerated by developing DNA markers associated with transcriptomes in breeding germplasm. Before the advent of next generation sequencing (NGS) technologies, the majority of sequencing data were generated by the Sanger sequencing method. By leveraging Sanger EST data, we have generated a wealth of genetic information for pepper including thousands of SNPs and Single Position Polymorphic (SPP) markers. To complement and enhance these resources, we applied NGS to three pepper genotypes: Maor, Early Jalapeño and Criollo de Morelos-334 (CM334) to identify SNPs and SSRs in the assembly of these three genotypes.ResultsTwo pepper transcriptome assemblies were developed with different purposes. The first reference sequence, assembled by CAP3 software, comprises 31,196 contigs from >125,000 Sanger-EST sequences that were mainly derived from a Korean F1-hybrid line, Bukang. Overlapping probes were designed for 30,815 unigenes to construct a pepper Affymetrix GeneChip® microarray for whole genome analyses. In addition, custom Python scripts were used to identify 4,236 SNPs in contigs of the assembly. A total of 2,489 simple sequence repeats (SSRs) were identified from the assembly, and primers were designed for the SSRs. Annotation of contigs using Blast2GO software resulted in information for 60% of the unigenes in the assembly. The second transcriptome assembly was constructed from more than 200 million Illumina Genome Analyzer II reads (80–120 nt) using a combination of Velvet, CLC workbench and CAP3 software packages. BWA, SAMtools and in-house Perl scripts were used to identify SNPs among three pepper genotypes. The SNPs were filtered to be at least 50 bp from any intron-exon junctions as well as flanking SNPs. More than 22,000 high-quality putative SNPs were identified. Using the MISA software, 10,398 SSR markers were also identified within the Illumina transcriptome assembly and primers were designed for the identified markers. The assembly was annotated by Blast2GO and 14,740 (12%) of annotated contigs were associated with functional proteins.ConclusionsBefore availability of pepper genome sequence, assembling transcriptomes of this economically important crop was required to generate thousands of high-quality molecular markers that could be used in breeding programs. In order to have a better understanding of the assembled sequences and to identify candidate genes underlying QTLs, we annotated the contigs of Sanger-EST and Illumina transcriptome assemblies. These and other information have been curated in a database that we have dedicated for pepper project.

Our ability to assemble complex genomes and construct ultradense genetic maps now allows the determination of recombination rates, translocations, and the extent of genomic collinearity between populations, species, and genera. We developed two ultradense genetic linkage maps for pepper from single-position polymorphisms (SPPs) identified de novo with a 30,173 unigene pepper genotyping array. The Capsicum frutescens × C. annuum interspecific and the C. annuum intraspecific genetic maps were constructed comprising 16,167 and 3,878 unigene markers in 2108 and 783 genetic bins, respectively. Accuracies of marker groupings and orders are validated by the high degree of collinearity between the two maps. Marker density was sufficient to locate the chromosomal breakpoint resulting in the P1/P8 translocation between C. frutescens and C. annuum to a single bin. The two maps aligned to the pepper genome showed varying marker density along the chromosomes. There were extensive chromosomal regions with suppressed recombination and reduced intraspecific marker density. These regions corresponded to the pronounced nonrecombining pericentromeric regions in tomato, a related Solanaceous species. Similar to tomato, the extent of reduced recombination appears to be more pronounced in pepper than in other plant species. Alignment of maps with the tomato and potato genomes shows the presence of previously known translocations and a translocation event that was not observed in previous genetic maps of pepper.

Phytophthora capsici is an important pepper (Capsicum annuum) pathogen causing fruit and root rot, and foliar blight in field and greenhouse production. Previously, an F6 recombinant inbred line population was evaluated for fruit rot susceptibility. Continuous variation among lines and partial and isolate-specific resistance were found. In this study, Phytophthora fruit rot resistance was mapped in the same F6 population between Criollo del Morelos 334 (CM334), a landrace from Mexico, and 'Early Jalapeno' using a high-density genetic map. Isolate-specific resistance was mapped independently in 63 of the lines evaluated and the two parents. Heritability of the resistance for each isolate at 3 and 5 days postinoculation (dpi) was high (h(2) = 0.63 to 0.68 and 0.74 to 0.83, respectively). Significant additive and epistatic quantitative trait loci (QTL) were identified for resistance to isolates OP97 and 13709 (3 and 5 dpi) and 12889 (3 dpi only). Mapping of fruit traits showed potential linkage with few disease resistance QTL. The partial fruit rot resistance from CM334 suggests that this may not be an ideal source for fruit rot resistance in pepper.

The identification of genes involved in variation of peach fruit quality would assist breeders to create new cultivars with improved fruit quality. Peach is a genetic and genomic model within the Rosaceae. A large quantity of useful data suitable for fine mapping using Single Nucleotide Polymorphisms (SNPs) from the peach genome sequence was used in this study. A set of 94 individuals from a peach germplasm collection was phenotyped and genotyped, including local Spanish and modern cultivars maintained at the Experimental Station of Aula Dei, Spain. Phenotypic evaluation based on agronomical, pomological and fruit quality traits was performed at least 3 years. A set of 4,558 out of a total of 8,144 SNPs markers developed by the Illumina Infinium BeadArray (v1.0) technology platform, covering the peach genome, were analyzed for population structure analysis and genome-wide association studies (GWAS). Population structure analysis identified two subpopulations, with admixture within them. While one subpopulation contains only modern cultivars, the other one is formed by local Spanish and several modern cultivars from international breeding programs. To test the marker trait associations between markers and phenotypic traits, four models comprising both general linear model (GLM) and mixed linear model (MLM) were selected. The MLM approach using co-ancestry values from population structure and kinship estimates (K model) identified a maximum of 347 significant associations between markers and traits. The associations found appeared to map within the interval where many candidate genes involved in different pathways are predicted in the peach genome. These results represent a promising situation for GWAS in the identification of SNP variants associated to fruit quality traits, potentially applicable in peach breeding programs.