A COSII genetic map of the pepper genome provides a detailed picture of synteny with tomato and new insights into recent chromosome evolution in the genus Capsicum

Unravelling inversions: Technological advances, challenges, and potential impact on crop breeding

Inversions, a type of chromosomal structural variation, significantly influence plant adaptation and gene functions by impacting gene expression and recombination rates. However, compared with other structural variations, their roles in functional biology and crop improvement remain largely unexplored. In this review, we highlight technological and methodological advancements that have allowed a comprehensive understanding of inversion variants through the pangenome framework and machine learning algorithms. Genome editing is an efficient method for inducing or reversing inversion mutations in plants, providing an effective mechanism to modify local recombination rates. Given the potential of inversions in crop breeding, we anticipate increasing attention on inversions from the scientific community in future research and breeding applications.

Genetic analysis and QTL mapping of domestication-related traits in chili pepper (Capsicum annuum L.)

Chili pepper (Capsicum annuum L.) is one of the oldest and most phenotypically diverse pre-Columbian crops of the Americas. Despite the abundance of genetic resources, the use of wild germplasm and landraces in chili pepper breeding is limited. A better understanding of the evolutionary history in chili peppers, particularly in the context of traits of agronomic interest, can contribute to future improvement and conservation of genetic resources. In this study, an F₂:₃ mapping population derived from a cross between a C. annuum wild accession (Chiltepin) and a cultivated variety (Puya) was used to identify genomic regions associated with 19 domestication and agronomic traits. A genetic map was constructed consisting of 1023 single nucleotide polymorphism (SNP) markers clustered into 12 linkage groups and spanning a total of 1,263.87 cM. A reciprocal translocation that differentiates the domesticated genome from its wild ancestor and other related species was identified between chromosomes 1 and 8. Quantitative trait locus (QTL) analysis detected 20 marker-trait associations for 13 phenotypes, from which 14 corresponded to previously identified loci, and six were novel genomic regions related to previously unexplored domestication-syndrome traits, including form of unripe fruit, seedlessness, deciduous fruit, and growth habit. Our results revealed that the genetic architecture of Capsicum domestication is similar to other domesticated species with few loci with large effects, the presence of QTLs clusters in different genomic regions, and the predominance of domesticated recessive alleles. Our analysis indicates the domestication process in chili pepper has also had an effect on traits not directly related to the domestication syndrome. The information obtained in this study provides a more complete understanding of the genetic basis of Capsicum domestication that can potentially guide strategies for the exploitation of wild alleles.

Induced mutation in ELONGATED HYPOCOTYL5 abolishes anthocyanin accumulation in the hypocotyl of pepper

The causal gene, CaHY5 of a chemical induced green-hypocotyl mutant was identified by molecular mapping. CaHY5 regulates anthocyanin accumulation by directly binding to the promoter of genes in anthocyanin pathway. Morphological markers at seedling stage are useful indicators for F₁ hybrid seeds screening. Pepper is a worldwide vegetable with diverse uses, and F₁ hybrids are popular in the pepper industry. Hypocotyl color is a useful marker to identify F₁ hybrid seeds. However, most pepper accessions have purple hypocotyl caused by anthocyanin accumulation, while green hypocotyl pepper accessions are rare. In this study, we identified a green hypocotyl mutant (e1898) from a pepper ethylmethanesulfonate (EMS) mutant library. By combining bulked segregant RNA-seq (BSR), genome resequencing and recombinant analysis, it was found that CaHY5 is the causal gene of this mutant. Virus-induced gene silencing (VIGS) of CaHY5 resulted in the decrease of anthocyanin accumulation in pepper hypocotyls. RNA-seq data showed that many genes related to anthocyanin biosynthesis and transport decreased significantly in the mutant. Yeast one-hybrid (Y1H) assays showed that CaHY5 can bind to the promoter of CaF3H, CaF3'5'H, CaDFR, CaANS and CaGST, which are important genes in anthocyanin biosynthesis or transport. Our results indicate that CaHY5 directly regulates anthocyanin biosynthesis and transport, thus governing anthocyanin accumulation in pepper hypocotyl. The mutant and gene identified in this work shall be valuable in the purity control of hybrid pepper seeds.

Identifying genomic regions targeted during eggplant domestication using transcriptome data

Identifying genes and traits that have diverged during domestication provides key information of importance for maintaining and even increasing yield and nutrients in existing crops. A “bottom-up” population genetics approach was used to identify signatures of selection across the eggplant genome, to better understand the process of domestication. RNA-seq data were obtained for 4 wild eggplants (Solanum insanum L.) and 16 domesticated eggplants (S. melongena L.) and mapped to the eggplant genome. Single-nucleotide polymorphism (SNPs) exhibiting signatures of selection in domesticates were identified as those exhibiting high F_ST between the 2 populations (evidence of significant divergence) and low π for the domesticated population (indicative of a selective sweep). Some of these regions appear to overlap with previously identified quantitative trait loci for domestication traits. Genes in regions of linkage disequilibrium surrounding these SNPs were searched against the Arabidopsis thaliana and tomato genomes to find orthologs. Subsequent gene ontology (GO) enrichment analysis identified over-representation of GO terms related to photosynthesis and response to the environment. This work reveals genomic changes involved in eggplant domestication and improvement, and how this compares to observed changes in the tomato genome, revealing shared chromosomal regions involved in the domestication of both species.

A reference chromosomal map of the hot chili pepper Capsicum pubescens cv. “locoto” (Solanaceae)

Abstract Capsicum pubescens is a cultivated hot chili pepper, consumed in Latin American cuisine as a distinctive ingredient, and popularly known as “locoto” or “rocoto”. This taxon is also an outstanding source of resistance to biotic and abiotic stresses as well as other valuable fruit traits for breeding of the worldwidely cultivated C. annuum and related species. In this study, the chromosome complement of C. pubescens cv. “locoto” (2n = 24) was deeply characterized through a sequential combination of conventional and molecular cytogenetics approaches comprising: Ag-NOR staining, heterochromatic fluorescent C-DAPI, DAPI/AMD-CMA/DA bandings, fluorescence in situ hybridization (FISH) of Capsicum-derived probes of the 5S and 18S-25S rRNA genes and different regions of spacers of the ribosomal unit, as well as telomeric probe. The markers identified were systematically combined with morphological karyotype parameters - number, size, centromeres, satellites - to produce a physical map which allowed the identification of several landmarks in each individual chromosome. The reference chromosomal map of C. pubescens here presented is the most comprehensively developed in Capsicum so far. It is envisioned that this chromosomal map will serve as a reference framework for the upcoming sequencing projects and as starting point to assist future genetic mapping of important agronomic traits.

Frequency, Origins, and Evolutionary Role of Chromosomal Inversions in Plants

Chromosomal inversions have the potential to play an important role in evolution by reducing recombination between favorable combinations of alleles. Until recently, however, most evidence for their likely importance derived from dipteran flies, whose giant larval salivary chromosomes aided early cytogenetic studies. The widespread application of new genomic technologies has revealed that inversions are ubiquitous across much of the plant and animal kingdoms. Here we review the rapidly accumulating literature on inversions in the plant kingdom and discuss what we have learned about their establishment and likely evolutionary role. We show that inversions are prevalent across a wide range of plant groups. We find that inversions are often associated with locally favored traits, as well as with traits that contribute to assortative mating, suggesting that they may be key to adaptation and speciation in the face of gene flow. We also discuss the role of inversions in sex chromosome formation, and explore possible parallels with inversion establishment on autosomes. The identification of inversion origins, as well as the causal variants within them, will advance our understanding of chromosomal evolution in plants.

Frequency, Origins, and Evolutionary Role of Chromosomal Inversions in Plants

Chromosomal inversions have the potential to play an important role in evolution by reducing recombination between favorable combinations of alleles. Until recently, however, most evidence for their likely importance derived from dipteran flies, whose giant larval salivary chromosomes aided early cytogenetic studies. The widespread application of new genomic technologies has revealed that inversions are ubiquitous across much of the plant and animal kingdoms. Here we review the rapidly accumulating literature on inversions in the plant kingdom and discuss what we have learned about their establishment and likely evolutionary role. We show that inversions are prevalent across a wide range of plant groups. We find that inversions are often associated with locally favored traits, as well as with traits that contribute to assortative mating, suggesting that they may be key to adaptation and speciation in the face of gene flow. We also discuss the role of inversions in sex chromosome formation, and explore possible parallels with inversion establishment on autosomes. The identification of inversion origins, as well as the causal variants within them, will advance our understanding of chromosomal evolution in plants.

Development and Application of InDel Markers for Capsicum spp. Based on Whole-Genome Re-Sequencing

Genome-wide identification of Insertion/Deletion polymorphisms (InDels) in Capsicum spp. was performed through comparing whole-genome re-sequencing data from two Capsicum accessions, C. annuum cv. G29 and C. frutescens cv. PBC688, with the reference genome sequence of C. annuum cv. CM334. In total, we identified 1,664,770 InDels between CM334 and PBC688, 533,523 between CM334 and G29, and 1,651,856 between PBC688 and G29. From these InDels, 1605 markers of 3-49 bp in length difference between PBC688 and G29 were selected for experimental validation: 1262 (78.6%) showed polymorphisms, 90 (5.6%) failed to amplify, and 298 (18.6%) were monomorphic. For further validation of these InDels, 288 markers were screened across five accessions representing five domesticated species. Of these assayed markers, 194 (67.4%) were polymorphic, 87 (30.2%) monomorphic and 7 (2.4%) failed. We developed three interspecific InDels, which associated with three genes and showed specific amplification in five domesticated species and clearly differentiated the interspecific hybrids. Thus, our novel PCR-based InDel markers provide high application value in germplasm classification, genetic research and marker-assisted breeding in Capsicum species.

Genome-wide analysis of genes encoding MBD domain-containing proteins from tomato suggest their role in fruit development and abiotic stress responses

In tomato, DNA methylation has an inhibitory effect on fruit ripening. The inhibition of DNA methyltransferase by 5-azacytidine results in premature fruit ripening. Methyl CpG binding domain (MBD) proteins are the readers of DNA methylation marks and help in the recruitment of chromatin-modifying enzymes which affect gene expression. Therefore, we investigate their contribution during fruit development. In this study, we identified and analyzed 18 putative genes of Solanum lycopersicum and Solanum pimpinellifolium encoding MBD proteins. We also identified tomato MBD syntelogs in Capsicum annum and Solanum tuberosum. Sixty-three MBD genes identified from four different species of solanaceae were classified into three groups. An analysis of the conserved domains in these proteins identified additional domains along with MBD motif. The transcript profiling of tomato MBDs in wild-type and two non-ripening mutants, rin and Nr, indicated constructive information regarding their involvement during fruit development. When we performed a stage-specific expression analysis during fruit ripening, a gradual decrease in transcript accumulation in the wild-type fruit was detected. However, a very low expression was observed in the ripening mutants. Furthermore, many ethylene-responsive cis-elements were found in SlMBD gene promoters, and some of them were induced in the presence of exogenous ethylene. Further, we detected the possible role of these MBDs in abiotic stresses. We found that few genes were differentially expressed under various abiotic stress conditions. Our results provide an evidence of the involvement of the tomato MBDs in fruit ripening and abiotic stress responses, which would be helpful in further studies on these genes in tomato fruit ripening.

Identification of Cucumber mosaic resistance 2 (cmr2) That Confers Resistance to a New Cucumber mosaic virus Isolate P1 (CMV-P1) in Pepper (Capsicum spp.)

Cucumber mosaic virus (CMV) is one of the most devastating phytopathogens of Capsicum. The single dominant resistance gene, Cucumber mosaic resistant 1 (Cmr1), that confers resistance to the CMV isolate P0 has been overcome by a new isolate (CMV-P1) after being deployed in pepper (Capsicum annuum) breeding for over 20 years. A recently identified Indian C. annuum cultivar, "Lam32," displays resistance to CMV-P1. In this study, we show that the resistance in "Lam32" is controlled by a single recessive gene, CMV resistance gene 2 (cmr2). We found that cmr2 conferred resistance to CMV strains including CMV-Korean, CMV-Fny, and CMV-P1, indicating that cmr2 provides a broad-spectrum type of resistance. We utilized two molecular mapping approaches to determine the chromosomal location of cmr2. Bulked segregant analysis (BSA) using amplified fragment-length polymorphism (AFLP) (BSA-AFLP) revealed one marker, cmvAFLP, located 16 cM from cmr2. BSA using the Affymetrix pepper array (BSA-Affy) identified a single-nucleotide polymorphism (SNP) marker (Affy4) located 2.3 cM from cmr2 on chromosome 8. We further screened a pepper germplasm collection of 4,197 accessions for additional CMV-P1 resistance sources and found that some accessions contained equivalent levels of resistance to that of "Lam32." Inheritance and allelism tests demonstrated that all the resistance sources examined contained cmr2. Our result thus provide genetic and molecular evidence that cmr2 is a single recessive gene that confers to pepper an unprecedented resistance to the dangerous new isolate CMV-P1 that had overcome Cmr1.

Comparative genome-wide characterization leading to simple sequence repeat marker development for Nicotiana

BACKGROUND

Simple sequence repeats (SSRs) are tandem repeats of DNA that have been used to develop robust genetic markers. These molecular markers are powerful tools for basic and applied studies such as molecular breeding. In the model plants in Nicotiana genus e.g. N. benthamiana, a comprehensive assessment of SSR content has become possible now because several Nicotiana genomes have been sequenced. We conducted a genome-wide SSR characterization and marker development across seven Nicotiana genomes.

RESULTS

Here, we initially characterized 2,483,032 SSRs (repeat units of 1-10 bp) from seven genomic sequences of Nicotiana and developed SSR markers using the GMATA® software package. Of investigated repeat units, mono-, di- and tri-nucleotide SSRs account for 98% of all SSRs in Nicotiana. More complex SSR motifs, although rare, are highly variable between Nicotiana genomes. A total of 1,224,048 non-redundant Nicotiana (NIX) markers were developed, of which 99.98% are novel. An efficient and uniform genotyping protocol for NIX markers was developed and validated. We created a web-based database of NIX marker information including amplicon sizes of alleles in each genome for downloading and online analysis.

CONCLUSIONS

The present work constitutes the first deep characterization of SSRs in seven genomes of Nicotiana, and the development of NIX markers for these SSRs. Our online marker database and an efficient genotyping protocol facilitate the application of these markers. The NIX markers greatly expand Nicotiana marker resources, thus providing a useful tool for future research and breeding. We demonstrate a novel protocol for SSR marker development and utilization at the whole genome scale that can be applied to any lineage of organisms. The Tobacco Markers & Primers Database (TMPD) is available at http://biodb.sdau.edu.cn/tmpd/index.html.

Mapping of fruit length related QTLs in interspecific cross (Capsicum annuum L. × Capsicum galapagoense Hunz.) of chilli

Fruit length in chilli is quantitatively inherited trait and selection based on phenotypic performance is tedious and time consuming. To detect QTLs determining fruit length in Capsicum spp., an interspecific F₂ mapping population was developed from the cross of C. annuum L. cv. 'FL 201' with C. galapagoense Hunz. accession 'TC 07245'. Fruit length in this cross showed a quantitative inheritance with the population depicting a symmetric distribution in histogram. To map quantitative trait loci (QTLs) for fruit length 400 SSR markers were surveyed on the parental genotypes but only 28 markers were observed to be polymorphic indicating less genetic diversity between the two Capsicum species. Polymorphic markers were then analyzed in F₂ population consisting of 210 plants and 24 of these markers were mapped on to three linkage groups (LGs): LG 1, LG 2 and LG 3. Two fruit length determining QTLs designated as paufl2.1 and paufl2.2 were identified and both the QTLs were mapped on to LG 2. The two QTLs together explained 21.78 per cent of the phenotypic variation. Apart from the two QTLs, positive alleles were detected in the small fruited parent 'TC 07245' which might be of potential use in chilli breeding programs.

Construction of a high-density genetic map and identification of QTLs for cucumber mosaic virus resistance in pepper (Capsicum annuum L.) using specific length amplified fragment sequencing (SLAF-seq)

Pepper (Capsicum) is one of the most important vegetable and spice crops. Aphid-transmitted cucumber mosaic virus (CMV) causes significant damage to pepper crops across the world. The genetic basis of CMV resistance in pepper is complex, and the mechanisms underlying resistance remain largely unknown. Here, we employed a SLAF-seq approach to generate a high-density genetic map of pepper. The map spanned 1,785.46 cM, containing 12,727 markers on 12 chromosomes, with a mean marker distance of 0.16 cM between adjacent markers. We used this map and the interval mapping (IM) and multiple QTL mapping (MQM) procedures to detect genetic regions associated with quantitative trait for CMV resistance. Three QTLs, qcmv11.1, qcmv11.2 and qcmv12.1, conferred resistance to CMV and showed trait variation of 10.2%, 19.2% and 7.3% respectively. Our results will help to develop markers linked to CMV-resistant QTLs to improve pepper resistance to CMV.

Reference quality assembly of the 3.5-Gb genome of Capsicum annuum from a single linked-read library

Linked-Read sequencing technology has recently been employed successfully for de novo assembly of human genomes, however, the utility of this technology for complex plant genomes is unproven. We evaluated the technology for this purpose by sequencing the 3.5-gigabase (Gb) diploid pepper (Capsicum annuum) genome with a single Linked-Read library. Plant genomes, including pepper, are characterized by long, highly similar repetitive sequences. Accordingly, significant effort is used to ensure that the sequenced plant is highly homozygous and the resulting assembly is a haploid consensus. With a phased assembly approach, we targeted a heterozygous F₁ derived from a wide cross to assess the ability to derive both haplotypes and characterize a pungency gene with a large insertion/deletion. The Supernova software generated a highly ordered, more contiguous sequence assembly than all currently available C. annuum reference genomes. Over 83% of the final assembly was anchored and oriented using four publicly available de novo linkage maps. A comparison of the annotation of conserved eukaryotic genes indicated the completeness of assembly. The validity of the phased assembly is further demonstrated with the complete recovery of both 2.5-Kb insertion/deletion haplotypes of the PUN1 locus in the F₁ sample that represents pungent and nonpungent peppers, as well as nearly full recovery of the BUSCO2 gene set within each of the two haplotypes. The most contiguous pepper genome assembly to date has been generated which demonstrates that Linked-Read library technology provides a tool to de novo assemble complex highly repetitive heterozygous plant genomes. This technology can provide an opportunity to cost-effectively develop high-quality genome assemblies for other complex plants and compare structural and gene differences through accurate haplotype reconstruction.

Conservation and Dispersion of Genes Conferring Resistance to Tomato Begomoviruses between Tomato and Pepper Genomes

In the present climate change scenario, controlling plant disease through exploitation of host plant resistance could contribute toward the sustainable crop production and global food security. In this respect, the identification of new sources of resistance and utilization of genetic diversity within the species may help in the generation of cultivars with improved disease resistance. Begomoviruses namely, Tomato yellow leaf curl virus (TYLCV) and Chilli leaf curl virus (ChLCV) are known to cause major yield losses in several economically important crop plants of the family Solanaceae. Though co-occurrence, association and synergistic interactions among these viruses in the host plants is reported, whether orthologous genetic loci in related host plants could be responsible for conferring resistance to these viruses has not been investigated yet. Several loci including Ty1, Ty2, Ty3, Ty4, and ty5 have been reported to confer resistance to leaf curl viruses in tomato. Here, we examined the pepper orthologous markers, corresponding to these QTL regions, for polymorphism between ChLCV susceptible and resistant genotypes of pepper. Further, to examine if the polymorphic markers are segregating with the disease resistance, Bulk Segregant Analysis (BSA) was performed on F₂ population derived from crosses between resistant and susceptible lines. However, none of the markers showed polymorphism in BSA suggesting that the tested markers are not linked to genes/QTLs responsible for conferring resistance to ChLCV in the selected genotypes. In silico analysis was performed to study the synteny and collinearity of genes located within these QTL regions in tomato and pepper genomes, which revealed that more than 60% genes located in Ty2 and Ty4, 13.71% genes in Ty1, 23.07% in Ty3, and 44.77% genes located within ty5 QTL region in tomato are conserved in pepper genome. However, despite such a high conservation in gene content, the linkage relationship in these regions seems to be greatly affected by gross rearrangements in both the species.

Expansion of sesquiterpene biosynthetic gene clusters in pepper confers nonhost resistance to the Irish potato famine pathogen

Chemical barriers contribute to nonhost resistance, which is defined as the resistance of an entire plant species to nonadapted pathogen species. However, the molecular basis of metabolic defense in nonhost resistance remains elusive. Here, we report genetic evidence for the essential role of phytoalexin capsidiol in nonhost resistance of pepper (Capsicum spp.) to potato late blight Phytophthora infestans using transcriptome and genome analyses. Two different genes for capsidiol biosynthesis, 5-epi-aristolochene synthase (EAS) and 5-epi-aristolochene-1,3-dihydroxylase (EAH), belong to multigene families. However, only a subset of EAS/EAH gene family members were highly induced upon P. infestans infection, which was associated with parallel accumulation of capsidiol in P. infestans-infected pepper. Silencing of EAS homologs in pepper resulted in a significant decrease in capsidiol accumulation and allowed the growth of nonadapted P. infestans that is highly sensitive to capsidiol. Phylogenetic and genomic analyses of EAS/EAH multigene families revealed that the emergence of pathogen-inducible EAS/EAH genes in Capsicum-specific genomic regions rendered pepper a nonhost of P. infestans. This study provides insights into evolutionary aspects of nonhost resistance based on the combination of a species-specific phytoalexin and sensitivity of nonadapted pathogens.

Genetic and biochemical analysis reveals linked QTLs determining natural variation for fruit post-harvest water loss in pepper (Capsicum)

Molecular markers linked to QTLs controlling post-harvest fruit water loss in pepper may be utilized to accelerate breeding for improved shelf life and inhibit over-ripening before harvest. Bell pepper (Capsicum annuum L.) is an important vegetable crop world-wide. However, marketing is limited by the relatively short shelf life of the fruit due to water loss and decay that occur during prolonged storage. Towards breeding pepper with reduced fruit post-harvest water loss (PWL), we studied the genetic, physiological and biochemical basis for natural variation of PWL. We performed quantitative trait locus (QTL) mapping of fruit PWL in multiple generations of an interspecific cross of pepper, which resulted in the identification of two linked QTLs on chromosome 10 that control the trait. We further developed near-isogenic lines (NILs) for characterization of the QTL effects. Transcriptome analysis of the NILs allowed the identification of candidate genes associated with fruit PWL-associated traits such as cuticle biosynthesis, cell wall metabolism and fruit ripening. Significant differences in PWL between the NILs in the immature fruit stage, differentially expressed cuticle-associated genes and differences in the content of specific chemical constituents of the fruit cuticle, indicated a likely influence of cuticle composition on the trait. Reduced PWL in the NILs was associated with delayed over-ripening before harvest, low total soluble solids before storage, and reduced fruit softening after storage. Our study enabled a better understanding of the genetic and biological processes controlling natural variation in fruit PWL in pepper. Furthermore, the genetic materials and molecular markers developed in this study may be utilized to breed peppers with improved shelf life and inhibited over-ripening before harvest.

Divergent evolution of multiple virus-resistance genes from a progenitor in Capsicum spp

Plants have evolved hundreds of nucleotide-binding and leucine-rich domain proteins (NLRs) as potential intracellular immune receptors, but the evolutionary mechanism leading to the ability to recognize specific pathogen effectors is elusive. Here, we cloned Pvr4 (a Potyvirus resistance gene in Capsicum annuum) and Tsw (a Tomato spotted wilt virus resistance gene in Capsicum chinense) via a genome-based approach using independent segregating populations. The genes both encode typical NLRs and are located at the same locus on pepper chromosome 10. Despite the fact that these two genes recognize completely different viral effectors, the genomic structures and coding sequences of the two genes are strikingly similar. Phylogenetic studies revealed that these two immune receptors diverged from a progenitor gene of a common ancestor. Our results suggest that sequence variations caused by gene duplication and neofunctionalization may underlie the evolution of the ability to specifically recognize different effectors. These findings thereby provide insight into the divergent evolution of plant immune receptors.

Genome-wide SNP discovery and population structure analysis in pepper (Capsicum annuum) using genotyping by sequencing

Background

Knowledge on population structure and genetic diversity in vegetable crops is essential for association mapping studies and genomic selection. Genotyping by sequencing (GBS) represents an innovative method for large scale SNP detection and genotyping of genetic resources. Herein we used the GBS approach for the genome-wide identification of SNPs in a collection of Capsicum spp. accessions and for the assessment of the level of genetic diversity in a subset of 222 cultivated pepper (Capsicum annum) genotypes.

Results

GBS analysis generated a total of 7,568,894 master tags, of which 43.4% uniquely aligned to the reference genome CM334. A total of 108,591 SNP markers were identified, of which 105,184 were in C. annuum accessions. In order to explore the genetic diversity of C. annuum and to select a minimal core set representing most of the total genetic variation with minimum redundancy, a subset of 222 C. annuum accessions were analysed using 32,950 high quality SNPs. Based on Bayesian and Hierarchical clustering it was possible to divide the collection into three clusters. Cluster I had the majority of varieties and landraces mainly from Southern and Northern Italy, and from Eastern Europe, whereas clusters II and III comprised accessions of different geographical origins. Considering the genome-wide genetic variation among the accessions included in cluster I, a second round of Bayesian (K = 3) and Hierarchical (K = 2) clustering was performed. These analysis showed that genotypes were grouped not only based on geographical origin, but also on fruit-related features.

Conclusions

GBS data has proven useful to assess the genetic diversity in a collection of C. annuum accessions. The high number of SNP markers, uniformly distributed on the 12 chromosomes, allowed the accessions to be distinguished according to geographical origin and fruit-related features. SNP markers and information on population structure developed in this study will undoubtedly support genome-wide association mapping studies and marker-assisted selection programs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3297-7) contains supplementary material, which is available to authorized users.

Development of a SNP array and its application to genetic mapping and diversity assessment in pepper (Capsicum spp.)

The development and application of single nucleotide polymorphisms (SNPs) is in its infancy for pepper. Here, a set of 15,000 SNPs were chosen from the resequencing data to develop an array for pepper with 12,720 loci being ultimately synthesized. Of these, 8,199 (~64.46%) SNPs were found to be scorable and covered ~81.18% of the whole genome. With this array, a high-density interspecific genetic map with 5,569 SNPs was constructed using 297 F2 individuals, and genetic diversity of a panel of 399 pepper elite/landrace lines was successfully characterized. Based on the genetic map, one major QTL, named Up12.1, was detected for the fruit orientation trait. A total of 65 protein-coding genes were predicted within this QTL region based on the current annotation of the Zunla-1 genome. In summary, the thousands of well-validated SNP markers, high-density genetic map and genetic diversity information will be useful for molecular genetics and innovative breeding in pepper. Furthermore, the mapping results lay foundation for isolating the genes underlying variation in fruit orientation of Capsicum.

Plant Genetic Background Increasing the Efficiency and Durability of Major Resistance Genes to Root-knot Nematodes Can Be Resolved into a Few Resistance QTLs

With the banning of most chemical nematicides, the control of root-knot nematodes (RKNs) in vegetable crops is now based essentially on the deployment of single, major resistance genes (R-genes). However, these genes are rare and their efficacy is threatened by the capacity of RKNs to adapt. In pepper, several dominant R-genes are effective against RKNs, and their efficacy and durability have been shown to be greater in a partially resistant genetic background. However, the genetic determinants of this partial resistance were unknown. Here, a quantitative trait loci (QTL) analysis was performed on the F2:3 population from the cross between Yolo Wonder, an accession considered partially resistant or resistant, depending on the RKN species, and Doux Long des Landes, a susceptible cultivar. A genetic linkage map was constructed from 130 F2 individuals, and the 130 F3 families were tested for resistance to the three main RKN species, Meloidogyne incognita, M. arenaria, and M. javanica. For the first time in the pepper-RKN pathosystem, four major QTLs were identified and mapped to two clusters. The cluster on chromosome P1 includes three tightly linked QTLs with specific effects against individual RKN species. The fourth QTL, providing specific resistance to M. javanica, mapped to pepper chromosome P9, which is known to carry multiple NBS-LRR repeats, together with major R-genes for resistance to nematodes and other pathogens. The newly discovered cluster on chromosome P1 has a broad spectrum of action with major additive effects on resistance. These data highlight the role of host QTLs involved in plant-RKN interactions and provide innovative potential for the breeding of new pepper cultivars or rootstocks combining quantitative resistance and major R-genes, to increase both the efficacy and durability of RKN control by resistance genes.

FISH-mapping of the 5S rDNA locus in chili peppers (Capsicum-Solanaceae)

We present here the physical mapping of the 5S rDNA locus in six wild and five cultivated taxa of Capsicum by means of a genus-specific FISH probe. In all taxa, a single 5S locus per haploid genome that persistently mapped onto the short arm of a unique metacentric chromosome pair at intercalar position, was found. 5S FISH signals of almost the same size and brightness intensity were observed in all the analyzed taxa. This is the first cytological characterization of the 5S in wild taxa of Capsicum by using a genus-derived probe, and the most exhaustive and comprehensive in the chili peppers up to now. The information provided here will aid the cytomolecular characterization of pepper germplasm to evaluate variability and can be instrumental to integrate physical, genetic and genomic maps already generated in the genus.

A HapMap leads to a Capsicum annuum SNP infinium array: a new tool for pepper breeding

The Capsicum genus (Pepper) is a part of the Solanacae family. It has been important in many cultures worldwide for its key nutritional components and uses as spices, medicines, ornamentals and vegetables. Worldwide population growth is associated with demand for more nutritionally valuable vegetables while contending with decreasing resources and available land. These conditions require increased efficiency in pepper breeding to deal with these imminent challenges. Through resequencing of inbred lines we have completed a valuable haplotype map (HapMap) for the pepper genome based on single-nucleotide polymorphisms (SNP). The identified SNPs were annotated and classified based on their gene annotation in the pepper draft genome sequence and phenotype of the sequenced inbred lines. A selection of one marker per gene model was utilized to create the PepperSNP16K array, which simultaneously genotyped 16 405 SNPs, of which 90.7% were found to be informative. A set of 84 inbred and hybrid lines and a mapping population of 90 interspecific F2 individuals were utilized to validate the array. Diversity analysis of the inbred lines shows a distinct separation of bell versus chile/hot pepper types and separates them into five distinct germplasm groups. The interspecific population created between Tabasco (C. frutescens chile type) and P4 (C. annuum blocky type) produced a linkage map with 5546 markers separated into 1361 bins on twelve 12 linkage groups representing 1392.3 cM. This publically available genotyping platform can be used to rapidly assess a large number of markers in a reproducible high-throughput manner for pepper. As a standardized tool for genetic analyses, the PepperSNP16K can be used worldwide to share findings and analyze QTLs for important traits leading to continued improvement of pepper for consumers. Data and information on the array are available through the Solanaceae Genomics Network.

New Insights on Eggplant/Tomato/Pepper Synteny and Identification of Eggplant and Pepper Orthologous QTL

Eggplant, pepper, and tomato are the most exploited berry-producing vegetables within the Solanaceae family. Their genomes differ in size, but each has 12 chromosomes which have undergone rearrangements causing a redistribution of loci. The genome sequences of all three species are available but differ in coverage, assembly quality and percentage of anchorage. Determining their syntenic relationship and QTL orthology will contribute to exploit genomic resources and genetic data for key agronomic traits. The syntenic analysis between tomato and pepper based on the alignment of 34,727 tomato CDS to the pepper genome sequence, identified 19,734 unique hits. The resulting synteny map confirmed the 14 inversions and 10 translocations previously documented, but also highlighted 3 new translocations and 4 major new inversions. Furthermore, each of the 12 chromosomes exhibited a number of rearrangements involving small regions of 0.5-0.7 Mbp. Due to high fragmentation of the publicly available eggplant genome sequence, physical localization of most eggplant QTL was not possible, thus, we compared the organization of the eggplant genetic map with the genome sequence of both tomato and pepper. The eggplant/tomato syntenic map confirmed all the 10 translocations but only 9 of the 14 known inversions; on the other hand, a newly detected inversion was recognized while another one was not confirmed. The eggplant/pepper syntenic map confirmed 10 translocations and 8 inversions already detected and suggested a putative new translocation. In order to perform the assessment of eggplant and pepper QTL orthology, the eggplant and pepper sequence-based markers located in their respective genetic map were aligned onto the pepper genome. GBrowse in pepper was used as reference platform for QTL positioning. A set of 151 pepper QTL were located as well as 212 eggplant QTL, including 76 major QTL (PVE ≥ 10%) affecting key agronomic traits. Most were confirmed to cluster in orthologous chromosomal regions. Our results highlight that the availability of genome sequences for an increasing number of crop species and the development of "ultra-dense" physical maps provide new and key tools for detailed syntenic and orthology studies between related plant species.

Combined Use of Molecular Markers and High-Resolution Melting (HRM) to Assess Chromosome Dosage in Potato Hybrids

In plants, the most widely used cytological techniques to assess parental genome contributions are based on in situ hybridization (FISH and GISH), but they are time-consuming and need specific expertise and equipment. Recent advances in genomics and molecular biology have made PCR-based markers a straightforward, affordable technique for chromosome typing. Here, we describe the development of a molecular assay that uses single-copy conserved ortholog set II (COSII)-based single nucleotide polymorphisms (SNPs) and the high-resolution melting (HRM) technique to assess the chromosome dosage of interspecific hybrids between a Solanum phureja-S. tuberosum diploid (2n = 2x = 24) hybrid and its wild relative S. commersonii. Screening and analysis of 45 COSII marker sequences allowed S. commersonii-specific SNPs to be identified for all 12 chromosomes. Combining the HRM technique with the establishment of synthetic DNA hybrids, SNP markers were successfully used to predict the expected parental chromosome ratio of 5 interspecific triploid hybrids. These results demonstrate the ability of this strategy to distinguish diverged genomes from each other, and to estimate chromosome dosage. The method could potentially be applied to any species as a tool to assess paternal to maternal ratios in the framework of a breeding program or following transformation techniques.

QTL mapping of thrips resistance in pepper

A QTL for thrips resistance on pepper chromosome 6 was identified and validated. This QTL affects thrips larval development and explains 50% of the variation. Thrips is one of the most damaging pests in pepper (Capsicum). Resistance to thrips was identified in Capsicum annuum. This study was aimed at the elucidation of the genetic background of thrips resistance in Capsicum through QTL mapping. The QTL analysis was carried out for Frankliniella occidentalis resistance in an F2 population consisting of 196 plants derived from an interspecific cross between the highly resistant C. annuum AC 1979 as female parent and the highly susceptible C. chinense 4661 as male parent. Fifty-seven SSR, 109 AFLP, and 5 SNP markers were used to construct a genetic map with a total length of 1636 cM. Damage caused by larvae and the survival of first and second instar larval stages observed in a no-choice test were used as parameters of resistance. Interval mapping detected one QTL for each of these parameters, all co-localizing near the same marker on chromosome 6. Use of this marker as co-factor in a multiple-QTL mapping analysis failed to uncover any additional QTLs. This QTL explained about 50% of the genetic variation, and the resistance allele of this QTL was inherited from the resistant parent. Thrips resistance was not linked to trichome density.

Ultra-High Density, Transcript-Based Genetic Maps of Pepper Define Recombination in the Genome and Synteny Among Related Species

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.

Molecular characterization of Pvr9 that confers a hypersensitive response to Pepper mottle virus (a potyvirus) in Nicotiana benthamiana

There are some R genes against potyviruses which were mapped in pepper. However, none of them has been characterized at the molecular level. In this study, we characterized Pvr9 which is an Rpi-blb2 ortholog from pepper and confers a hypersensitive response to Pepper mottle virus (PepMoV) in a transient expression system in Nicotiana benthamiana. This gene putatively encoded for 1298 amino acids and is located on pepper chromosome 6. PepMoV NIb was the elicitor of the Pvr9-mediated hypersensitive response. NIb from several other potyviruses also elicited the hypersensitive response. Inoculation of pepper with PepMoV resulted in a minor increase in Pvr9 transcription in the resistant cultivar CM334 and a slight down-regulation in the susceptible cultivar Floral Gem. The 5' upstream region of Pvr9 from cultivar CM334 had higher transcription activity than the region from cultivar Floral Gem. The cultivars CM334 and Floral Gem had non-functional Pvr9 homologs with loss-of-function mutations.

CaAP2 transcription factor is a candidate gene for a flowering repressor and a candidate for controlling natural variation of flowering time in Capsicum annuum

The APETALA2 transcription factor homolog CaAP2 is a candidate gene for a flowering repressor in pepper, as revealed by induced-mutation phenotype, and a candidate underlying a major QTL controlling natural variation in flowering time. To decipher the genetic control of transition to flowering in pepper (Capsicum spp.) and determine the extent of gene function conservation compared to model species, we isolated and characterized several ethyl methanesulfonate (EMS)-induced mutants that vary in their flowering time compared to the wild type. In the present study, we report on the isolation of an early-flowering mutant that flowers after four leaves on the primary stem compared to nine leaves in the wild-type 'Maor'. By genetic mapping and sequencing of putative candidate genes linked to the mutant phenotype, we identified a member of the APETALA2 (AP2) transcription factor family, CaAP2, which was disrupted in the early-flowering mutant. CaAP2 is a likely ortholog of AP2 that functions as a repressor of flowering in Arabidopsis. To test whether CaAP2 has an effect on controlling natural variation in the transition to flowering in pepper, we performed QTL mapping for flowering time in a cross between early and late-flowering C. annuum accessions. We identified a major QTL in a region of chromosome 2 in which CaAP2 was the most significant marker, explaining 52 % of the phenotypic variation of the trait. Sequence comparison of the CaAP2 open reading frames in the two parents used for QTL mapping did not reveal significant variation. In contrast, significant differences in expression level of CaAP2 were detected between near-isogenic lines that differ for the flowering time QTL, supporting the putative function of CaAP2 as a major repressor of flowering in pepper.

Construction of an interspecific genetic map based on InDel and SSR for mapping the QTLs affecting the initiation of flower primordia in pepper (Capsicum spp.)

Re-sequencing permits the mining of genome-wide variations on a large scale and provides excellent resources for the research community. To accelerate the development and application of molecular markers and identify the QTLs affecting the flowering time-related trait in pepper, a total of 1,038 pairs of InDel and 674 SSR primers from different sources were used for genetic mapping using the F2 population (n = 154) derived from a cross between BA3 (C. annuum) and YNXML (C. frutescens). Of these, a total of 224 simple PCR-based markers, including 129 InDels and 95 SSRs, were validated and integrated into a map, which was designated as the BY map. The BY map consisted of 13 linkage groups (LGs) and spanned a total genetic distance of 1,249.77 cM with an average marker distance of 5.60 cM. Comparative analysis of the genetic and physical map based on the anchored markers showed that the BY map covered nearly the whole pepper genome. Based on the BY map, one major and five minor QTLs affecting the number of leaves on the primary axis (Nle) were detected on chromosomes P2, P7, P10 and P11 in 2012. The major QTL on P2 was confirmed based on another subset of the same F2 population (n = 147) in 2014 with selective genotyping of markers from the BY map. With the accomplishment of pepper whole genome sequencing and annotations (release 2.0), 153 candidate genes were predicted to embed in the Nle2.2 region, of which 12 important flowering related genes were obtained. The InDel/SSR-based interspecific genetic map, QTLs and candidate genes obtained by the present study will be useful for the downstream isolation of flowering time-related gene and other genetic applications for pepper.

Development of molecular markers tightly linked to Pvr4 gene in pepper using next-generation sequencing

It is imperative to identify highly polymorphic and tightly linked markers of a known trait for molecular marker-assisted selection. Potyvirus resistance 4 (Pvr4) locus in pepper confers resistance to three pathotypes of potato virus Y and to pepper mottle virus. We describe the use of next-generation sequencing technology to generate molecular markers tightly linked to Pvr4. Initially, comparative genomics was carried out, and a syntenic region of tomato on chromosome ten was used to generate PCR-based markers and map Pvr4. Subsequently, the genomic sequence of pepper was used, and more than 5000 single-nucleotide variants (SNVs) were identified within the interval. In addition, we identified nucleotide binding site-leucine-rich repeat-type disease resistance genes within the interval. Several of these SNVs were converted to molecular markers desirable for large-scale molecular breeding programmes.

Genetic analysis to identify good combiners for ToLCV resistance and yield components in tomato using interspecific hybridization

The interspecific hybridization for tomato leaf curl virus (ToLCV) resistance was carried out among 10 genetically diverse tomato genotypes (diversified by 50 SSR markers). Among the 10 parents, four susceptible cultivars of Solanum lycopersicum were crossed with six resistant wilds, such as S. pimpinellifolium, S. habrochaites, S. chemielewskii, S. ceraseforme, S. peruvianum and S. chilense in a line x tester mating design. All the 24 hybrids and their parents were grown in the field and glasshouse conditions to determine the general-combining abilities (GCA) and specific-combining abilities (SCA). The variances due to SCA and GCA showed both additive and nonadditive gene effects. Based on GCA estimates, EC-520061 and WIR-5032 were good general combiners while based on SCA estimates, PBC x EC-520061 and PBC x EC-521080 were best specific combiners for coefficient of infection and fruit yield per plant in both the environments. These lines could be selected and utilized in ToLCV resistance and high yield breeding programme for improving the traits.

An InDel-based linkage map of hot pepper (Capsicum annuum)

Two independent pepper (Capsicum annuum) genomes were published recently, opening a new era of molecular genetics research on pepper. However, pepper molecular marker technologies are still mainly focusing on the simple sequence repeats derived from public database or genomic library. The development and application of the third generation marker system such as single nucleotide polymorphisms, structure variations as well as insertion/deletion polymorphisms (InDels) is still in its infancy. In the present study, we developed InDel markers for pepper genetic mapping with the convenience of two whole-genome re-sequenced inbred lines BA3 (C. annuum) and B702 (C. annuum). A total of 154,519 and 149,755 InDel (1-5 bp) sites were identified for BA3 and B702, respectively, by the alignment of re-sequencing reads to Zunla-1 reference genome. Then, 14,498 InDel sites (only 4 and 5 bp) that are different between BA3 and B702 were predicted. Finally, within a random set of 1,000 primer pairs, 251 InDel markers were validated and mapped onto a linkage map using F₂ population derived from the intraspecific cross BA3 × B702. The first InDel-based map, named as BB-InDel map, consisted of 12 linkage groups, covered a genetic distance of 1,178.01 cM and the average distance between bin markers was 5.01 cM. Compared to the Zunla-1 reference physical map, high consistency was observed on all 12 chromosomes, and the total length of scaffold anchored and physical distance covered by this map was 299.66 and 2,558.68 Mb, respectively, which accounted for 8.95 and 76.38 % of the Zunla-1 reference genome (3.35 Gb), respectively. Furthermore, 37 scaffolds (total length of 36.21 Mb) from the pseudo-chromosome (P0) of the current genome assembly were newly assigned to the corresponding chromosomes by 40 InDel markers. Thus, this map provided good genome coverage and would be useful for basic and applied research in pepper.

New insights into Capsicum spp relatedness and the diversification process of Capsicum annuum in Spain

The successful exploitation of germplasm banks, harbouring plant genetic resources indispensable for plant breeding, will depend on our ability to characterize their genetic diversity. The Vegetable Germplasm Bank of Zaragoza (BGHZ) (Spain) holds an important Capsicum annuum collection, where most of the Spanish pepper variability is represented, as well as several accessions of other domesticated and non-domesticated Capsicum spp from all over the five continents. In the present work, a total of 51 C. annuum landraces (mainly from Spain) and 51 accessions from nine Capsicum species maintained at the BGHZ were evaluated using 39 microsatellite (SSR) markers spanning the whole genome. The 39 polymorphic markers allowed the detection of 381 alleles, with an average of 9.8 alleles per locus. A sizeable proportion of alleles (41.2%) were recorded as specific alleles and the majority of these were present at very low frequencies (rare alleles). Multivariate and model-based analyses partitioned the collection in seven clusters comprising the ten different Capsicum spp analysed: C. annuum, C. chinense, C. frutescens, C. pubescens, C. bacatum, C. chacoense and C. eximium. The data clearly showed the close relationships between C. chinense and C. frutescens. C. cardenasii and C. eximium were indistinguishable as a single, morphologically variable species. Moreover, C. chacoense was placed between C. baccatum and C. pubescens complexes. The C. annuum group was structured into three main clusters, mostly according to the pepper fruit shape, size and potential pungency. Results suggest that the diversification of C. annuum in Spain may occur from a rather limited gene pool, still represented by few landraces with ancestral traits. This ancient population would suffer from local selection at the distinct geographical regions of Spain, giving way to pungent and elongated fruited peppers in the South and Center, while sweet blocky and triangular types in Northern Spain.

Comparative analysis of the radish genome based on a conserved ortholog set (COS) of Brassica

This manuscript provides a Brassica conserved ortholog set (COS) that can be used as diagnostic cross-species markers as well as tools for genetic mapping and genome comparison of the Brassicaceae. A conserved ortholog set (COS) is a collection of genes that are conserved in both sequence and copy number between closely related genomes. COS is a useful resource for developing gene-based markers and is suitable for comparative genome mapping. We developed a COS for Brassica based on proteome comparisons of Arabidopsis thaliana, B. rapa, and B. oleracea to establish a basis for comparative genome analysis of crop species in the Brassicaceae. A total of 1,194 conserved orthologous single-copy genes were identified from the genomes based on whole-genome BLASTP analysis. Gene ontology analysis showed that most of them encoded proteins with unknown function and chloroplast-related genes were enriched. In addition, 152 Brassica COS primer sets were applied to 16 crop and wild species of the Brassicaceae and 57.9-92.8 % of them were successfully amplified across the species representing that a Brassica COS can be used as diagnostic cross-species markers of diverse Brassica species. We constructed a genetic map of Raphanus sativus by analyzing the segregation of 322 COS genes in an F2 population (93 individuals) of Korean cultivars (WK10039 × WK10024). Comparative genome analysis based on the COS genes showed conserved genome structures between R. sativus and B. rapa with lineage-specific rearrangement and fractionation of triplicated subgenome blocks indicating close evolutionary relationship and differentiation of the genomes. The Brassica COS developed in this study will play an important role in genetic, genomic, and breeding studies of crop Brassicaceae species.

Chromosomal evolution and patterns of introgression in helianthus

Knowledge of the nature and extent of karyotypic differences between species provides insight into the evolutionary history of the genomes in question and, in the case of closely related species, the potential for genetic exchange between taxa. We constructed high-density genetic maps of the silverleaf sunflower (Helianthus argophyllus) and Algodones Dune sunflower (H. niveus ssp. tephrodes) genomes and compared them to a consensus map of cultivated sunflower (H. annuus) to identify chromosomal rearrangements between species. The genetic maps of H. argophyllus and H. niveus ssp. tephrodes included 17 linkage groups each and spanned 1337 and 1478 cM, respectively. Comparative analyses revealed greater divergence between H. annuus and H. niveus ssp. tephrodes (13 inverted segments, 18 translocated segments) than between H. annuus and H. argophyllus (10 inverted segments, 8 translocated segments), consistent with their known phylogenetic relationships. Marker order was conserved across much of the genome, with 83 and 64% of the H. argophyllus and H. niveus ssp. tephrodes genomes, respectively, being syntenic with H. annuus. Population genomic analyses between H. annuus and H. argophyllus, which are sympatric across a portion of the natural range of H. annuus, revealed significantly elevated genetic structure in rearranged portions of the genome, indicating that such rearrangements are associated with restricted gene flow between these two species.

Whole-genome sequencing of cultivated and wild peppers provides insights into Capsicum domestication and specialization

As an economic crop, pepper satisfies people's spicy taste and has medicinal uses worldwide. To gain a better understanding of Capsicum evolution, domestication, and specialization, we present here the genome sequence of the cultivated pepper Zunla-1 (C. annuum L.) and its wild progenitor Chiltepin (C. annuum var. glabriusculum). We estimate that the pepper genome expanded ∼0.3 Mya (with respect to the genome of other Solanaceae) by a rapid amplification of retrotransposons elements, resulting in a genome comprised of ∼81% repetitive sequences. Approximately 79% of 3.48-Gb scaffolds containing 34,476 protein-coding genes were anchored to chromosomes by a high-density genetic map. Comparison of cultivated and wild pepper genomes with 20 resequencing accessions revealed molecular footprints of artificial selection, providing us with a list of candidate domestication genes. We also found that dosage compensation effect of tandem duplication genes probably contributed to the pungent diversification in pepper. The Capsicum reference genome provides crucial information for the study of not only the evolution of the pepper genome but also, the Solanaceae family, and it will facilitate the establishment of more effective pepper breeding programs.

Quantitative trait loci from the host genetic background modulate the durability of a resistance gene: a rational basis for sustainable resistance breeding in plants

The combination of major resistance genes with quantitative resistance factors is hypothesized as a promising breeding strategy to preserve the durability of resistant cultivar, as recently observed in different pathosystems. Using the pepper (Capsicum annuum)/Potato virus Y (PVY, genus Potyvirus) pathosystem, we aimed at identifying plant genetic factors directly affecting the frequency of virus adaptation to the major resistance gene pvr2(3) and at comparing them with genetic factors affecting quantitative resistance. The resistance breakdown frequency was a highly heritable trait (h(2)=0.87). Four loci including additive quantitative trait loci (QTLs) and epistatic interactions explained together 70% of the variance of pvr2(3) breakdown frequency. Three of the four QTLs controlling pvr2(3) breakdown frequency were also involved in quantitative resistance, strongly suggesting that QTLs controlling quantitative resistance have a pleiotropic effect on the durability of the major resistance gene. With the first mapping of QTLs directly affecting resistance durability, this study provides a rationale for sustainable resistance breeding. Surprisingly, a genetic trade-off was observed between the durability of PVY resistance controlled by pvr2(3) and the spectrum of the resistance against different potyviruses. This trade-off seemed to have been resolved by the combination of minor-effect durability QTLs under long-term farmer selection.

Capturing flavors from Capsicum baccatum by introgression in sweet pepper

Biochemical characterization in combination with genetic analyses in BC 2 S 1 plants and near-isogenic lines led to the detection and validation of C. baccatum loci affecting flavor, terpenoid content and Brix level. The species Capsicum baccatum includes the most common hot peppers of the Andean cuisine, known for their rich variation in flavors and aromas. So far the C. baccatum genetic variation remained merely concealed for Capsicum annuum breeding, due to post-fertilization genetic barriers encountered in interspecific hybridization. However, to exploit the potential flavor wealth of C. baccatum we combined interspecific crossing with embryo rescue, resulting in a multi-parent BC2S1 population. Volatile and non-volatile compounds plus some physical characters were measured in mature fruits, in combination with taste evaluation by a sensory panel. An enormous variation in biochemical composition and sensory attributes was found, with almost all traits showing transgression. A population-specific genetic linkage map was developed for QTL mapping. BC2S1 QTLs were validated in an experiment with near-isogenic lines, resulting in confirmed genetic effects for physical, biochemical and sensory traits. Three findings are described in more detail: (1) A small C. baccatum LG3 introgression caused an extraordinary effect on flavor, resulting in significantly higher scores for the attributes aroma, flowers, spices, celery and chives. In an attempt to identify the responsible biochemical compounds few consistently up- and down-regulated metabolites were detected. (2) Two introgressions (LG10.1 and LG1) had major effects on terpenoid content of mature fruits, affecting at least 15 different monoterpenes. (3) A second LG3 fragment resulted in a strong increase in Brix without negative effects on fruit size. The mapping strategy, the potential application of studied traits and perspectives for breeding are discussed.

Conserved miRNAs and their response to salt stress in wild eggplant Solanum linnaeanum roots

The Solanaceae family includes some important vegetable crops, and they often suffer from salinity stress. Some miRNAs have been identified to regulate gene expression in plant response to salt stress; however, little is known about the involvement of miRNAs in Solanaceae species. To identify salt-responsive miRNAs, high-throughput sequencing was used to sequence libraries constructed from roots of the salt tolerant species, Solanum linnaeanum, treated with and without NaCl. The sequencing identified 98 conserved miRNAs corresponding to 37 families, and some of these miRNAs and their expression were verified by quantitative real-time PCR. Under the salt stress, 11 of the miRNAs were down-regulated, and 3 of the miRNAs were up-regulated. Potential targets of the salt-responsive miRNAs were predicted to be involved in diverse cellular processes in plants. This investigation provides valuable information for functional characterization of miRNAs in S. linnaeanum, and would be useful for developing strategies for the genetic improvement of the Solanaceae crops.

Genetic mapping of semi-polar metabolites in pepper fruits (Capsicum sp.): towards unravelling the molecular regulation of flavonoid quantitative trait loci

Untargeted LCMS profiling of semi-polar metabolites followed by metabolite quantitative trait locus (mQTL) analysis was performed in ripe pepper fruits of 113 F2 plants derived from a cross between Capsicum annuum AC1979 (no. 19) and Capsicum chinense No. 4661 Selection (no. 18). The parental accessions were selected based on their variation in fruit morphological characteristics and fruit content of some target phytonutrients. Clear segregation of fruit colour and fruit metabolite profiles was observed in the F2 population. The F2 plants formed three clusters based on their metabolite profiles. Of the total of 542 metabolites, 52 could be annotated, including a range of flavonoids, such as flavone C-glycosides, flavonol O-glycosides and naringenin chalcone, as well as several phenylpropanoids, a capsaicin analogue, fatty acid derivatives and amino acid derivatives. Interval mapping revealed 279 mQTLs in total. Two mQTL hotspots were found on chromosome 9. These two chromosomal regions regulated the relative levels of 35 and 103 metabolites, respectively. Analysis also revealed an mQTL for a capsaicin analogue, located on chromosome 7. Confirmation of flavonoid mQTLs using a set of six flavonoid candidate gene markers and their corresponding expression data (expression QTLs) indicated the Ca-MYB12 transcription factor gene on chromosome 1 and the gene encoding flavone synthase (FS-2) on chromosome 6 as likely causative genes determining the variation in naringenin chalcone and flavone C-glycosides, respectively, in this population. The combination of large-scale metabolite profiling and QTL analysis provided valuable insight into the genomic regions and genes important for the production of (secondary) metabolites in pepper fruit. This will impact breeding strategies aimed at optimising the content of specific metabolites in pepper fruit.

De novo transcriptome assembly and novel microsatellite marker information in Capsicum annuum varieties Saengryeg 211 and Saengryeg 213

BACKGROUND

Pepper, Capsicum annuum L., Solanaceae, is a major staple economically important vegetable crop worldwide. Limited functional genomics resources and whole genome association studies could be substantially improved through the application of molecular approach for the characterization of gene content and identification of molecular markers. The massive parallel pyrosequencing of two pepper varieties, the highly pungent, Saengryeg 211, and the non-pungent, Saengryeg 213, including de novo transcriptome assembly, functional annotation, and in silico discovery of potential molecular markers is described. We performed 454 GS-FLX Titanium sequencing of polyA-selected and normalized cDNA libraries generated from a single pool of transcripts obtained from mature fruits of two pepper varieties.

RESULTS

A single 454 pyrosequencing run generated 361,671 and 274,269 reads totaling 164.49 and 124.60 Mb of sequence data (average read length of 454 nucleotides), which assembled into 23,821 and 17,813 isotigs and 18,147 and 15,129 singletons for both varieties, respectively. These reads were organized into 20,352 and 15,781 'isogroups' for both varieties. Assembled sequences were functionally annotated based on homology to genes in multiple public databases and assigned with Gene Ontology (GO) terms. Sequence variants analyses identified a total of 3,766 and 2,431 potential (Simple Sequence Repeat) SSR motifs for microsatellite analysis for both varieties, where trinucleotide was the most common repeat unit (84%), followed by di (9.9%), hexa (4.1%) and pentanucleotide repeats (2.1%). GAA repeat (8.6%) was the most frequent repeat motif, followed by TGG (7.2%), TTC (6.5%), and CAG (6.2%).

CONCLUSIONS

High-throughput transcriptome assembly, annotation and large scale of SSR marker discovery has been achieved using next generation sequencing (NGS) of two pepper varieties. These valuable informations for functional genomics resource shall help to further improve the pepper breeding efforts with respect to genetic linkage maps, QTL mapping and marker-assisted trait selection.

CaHAM is autoregulated and regulates CaSTM expression and is required for shoot apical meristem organization in pepper

The angiosperm shoot apical meristem (SAM) is characterized by tightly organized cell layers and zones. The SAM's organization allows it to maintain its indeterminate nature while producing determinate lateral organs. Alterations in SAM gene expression partly account for the immense diversity in plant architecture. The GRAS protein family gene HAIRY MERISTEM (HAM) is an important regulator of SAM organization in Petunia and Arabidopsis. Here we describe CaHAM loss-of-function pepper mutants characterized by an arrested SAM following the formation of several leaves on the primary stem, complete inhibition of axillary meristem development, an expanded tunica domain and trichome formation on the SAM epidermis. CaHAM is expressed in the periphery of the SAM and in the vasculature of young leaves throughout plant development, reaching its highest level in the reproductive growth stage. Analysis of the effect of CaHAM loss-of-function on its own expression showed that CaHAM is negatively autoregulated. Furthermore, CaHAM negatively regulates the expression level and pattern of pepper SHOOT MERISTEMLESS (CaSTM), which is required to maintain the SAM in an undifferentiated state. We conclude that CaHAM is regulated to achieve adjusted functional levels and has a conserved role in controlling SAM maintenance, organization and axillary meristem formation.

Development of universal genetic markers based on single-copy orthologous (COSII) genes in Poaceae

KEY MESSAGE : We develop a set of universal genetic markers based on single-copy orthologous (COSII) genes in Poaceae. Being evolutionary conserved, single-copy orthologous (COSII) genes are particularly useful in comparative mapping and phylogenetic investigation among species. In this study, we identified 2,684 COSII genes based on five sequenced Poaceae genomes including rice, maize, sorghum, foxtail millet, and brachypodium, and then developed 1,072 COSII markers whose transferability and polymorphism among five bamboo species were further evaluated with 46 pairs of randomly selected primers. 91.3 % of the 46 primers obtained clear amplification in at least one bamboo species, and 65.2 % of them produced polymorphism in more than one species. We also used 42 of them to construct the phylogeny for the five bamboo species, and it might reflect more precise evolutionary relationship than the one based on the vegetative morphology. The results indicated a promising prospect of applying these markers to the investigation of genetic diversity and the classification of Poaceae. To ease and facilitate access of the information of common interest to readers, a web-based database of the COSII markers is provided ( http://www.sicau.edu.cn/web/yms/PCOSWeb/PCOS.html ).

DNA marker applications to molecular genetics and genomics in tomato

Tomato is an important crop and regarded as an experimental model of the Solanaceae family and of fruiting plants in general. To enhance breeding efficiency and advance the field of genetics, tomato has been subjected to DNA marker studies as one of the earliest targets in plants. The developed DNA markers have been applied to the construction of genetic linkage maps and the resultant maps have contributed to quantitative trait locus (QTL) and gene mappings for agronomically important traits, as well as to comparative genomics of Solanaceae. The recently released whole genome sequences of tomato enable us to develop large numbers of DNA markers comparatively easily, and even promote new genotyping methods without DNA markers. In addition, databases for genomes, DNA markers, genetic linkage maps and other omics data, e.g., transcriptome, proteome, metabolome and phenome information, will provide useful information for molecular breeding in tomatoes. The use of DNA marker technologies in conjunction with new breeding techniques will promise to advance tomato breeding.

Identification of QTLs for capsaicinoids, fruit quality, and plant architecture-related traits in an interspecific Capsicum RIL population

Quantitative trait loci (QTL) analyses in pepper are common for horticultural, disease resistance, and fruit quality traits; although none of the studies to date have used sequence-based markers associated with genes. In this study we measured plant architectural, phenological, and fruit quality traits in a pepper mapping population consisting of 92 recombinant inbred lines derived from a cross between Capsicum frutescens acc. 2814-6 and C. annuum var. NuMexRNAKY. Phenotypic measurements were correlated to loci in a high-density EST-based genetic map. In total, 96 QTL were identified for 38 traits, including 12 QTL associated with capsaicinoid levels. Twenty-one loci showed correlation among seemingly unrelated phenotypic categories, highlighting tight linkage or shared genetics between previously unassociated traits in pepper.

Development of Capsicum EST-SSR markers for species identification and in silico mapping onto the tomato genome sequence

Capsicum spp. are widely cultivated for use as vegetables and spices. The Kihara Institute for Biological Research, Yokohama City University, Japan, has stocks of approximately 800 lines of Capsicum spp. collected from various regions of Central and South America, the regions of origin for Capsicum spp. In this study, 5,751 primer pairs for simple sequence repeat markers, based on 118,060 publicly available sequences of expressed sequence tags of Capsicum annuum, were designed and subjected to a similarity search against the genomic sequence of tomato, a model Solanaceae species. Nucleotide sequences spanning 2,245 C. annuum markers were successfully mapped onto the tomato genome, and 96 of these, which spanned the entire tomato genome, were selected for further analysis. In genotyping analysis, 60 out of the 77 markers that produced specific DNA amplicons showed polymorphism among the Capsicum lines examined. On the basis of the resulting data, the 192 tested lines were grouped into five main clusters. The additional sequencing analysis of the plastid genes, matK and rbcL, divided the resources into three groups. As a result, 19 marker loci exhibited genotypes specific to species and cluster, suggesting that the DNA markers are useful for species identification. Information on the DNA markers will contribute to Capsicum genetics, genomics, and breeding. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1007/s11032-012-9774-z) contains supplementary material, which is available to authorized users.

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.