Association Mapping Analysis for Fruit Quality Traits in Prunus persica Using SNP Markers

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.

RNA sequencing provides insights into the evolution of lettuce and the regulation of flavonoid biosynthesis

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.

De novo assembly of the pepper transcriptome (Capsicum annuum): a benchmark for in silico discovery of SNPs, SSRs and candidate genes

Background

Molecular 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.

Results

Two 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 F₁-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.

Conclusions

Before 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.