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

Development of a 45K pepper GBTS liquid-phase gene chip and its application in genome-wide association studies

Introduction

Pepper (Capsicum spp.) is a vegetable that is cultivated globally and has undergone extensive domestication, leading to a significant diversification in its agronomic traits. With the advancement of genomics in pepper and the reduction in sequencing costs, the high-throughput detection of single nucleotide polymorphisms (SNPs) and small insertions-deletions (indels) has become increasingly critical for analyzing pepper germplasms and improving breeding programs. As a result, there is a pressing need for a cost-effective, high-throughput, and versatile technique suitable for both foreground and background selection in pepper breeding.

Methods

In the present study, Python-based web scraping scripts were utilized to systematically extract data from published literatures and relevant sequence databases focusing on pepper genomes. Subsequent to data extraction, SNPs and indels were meticulously identified and filtered. This process culminated in the delineation of core polymorphic sites, which were instrumental in the development of specific probes. Following this, comprehensive phenotypic and genotypic analyses were conducted on a diverse collection of 420 pepper germplasms. Concurrently, a genome-wide association study (GWAS) was conducted to elucidate the genetic determinants of helical fruit shape in peppers.

Results

In this study, a 45K pepper Genotyping-By-Target-Sequencing (GBTS) liquid-phase gene chip was developed on the GenoBaits platform. This chip is composed of 45,389 probes, of which 42,535 are derived from core polymorphic sites (CPS) in the background genetic landscape, while 2,854 are associated with foreground agronomic traits, spanning across 43 traits. The CPS probes are spaced at an average interval of 68 Kb. We have assessed the performance of this chip on 420 pepper germplasms, with successful capture of target DNA fragments by 45,387 probes. Furthermore, the probe capture ratio surpassed 70% in 410 of the 420 germplasms tested. Using this chip, we have efficiently genotyped 273 germplasms for spiciness levels and elucidated the genetic relationships among 410 pepper germplasms. Our results allowed for precise clustering of sister lines and C. chinense germplasms. In addition, through a GWAS for helical fruit shape, we identified three quantitative trait loci (QTLs): heli2.1, heli11.1, and heli11.2. Within the heli11.1 QTL, a gene encoding the tubulin alpha chain was identified, suggesting its potential role in the helical growth pattern of pepper fruits.

Discussion

In summary, the 45K pepper GBTS liquid-phase gene chip offers robust detection of polymorphic sites and is a promising tool for advancing research into pepper germplasm and the breeding of new pepper varieties.

Revitalizing agriculture: next-generation genotyping and -omics technologies enabling molecular prediction of resilient traits in the Solanaceae family

This review highlights -omics research in Solanaceae family, with a particular focus on resilient traits. Extensive research has enriched our understanding of Solanaceae genomics and genetics, with historical varietal development mainly focusing on disease resistance and cultivar improvement but shifting the emphasis towards unveiling resilience mechanisms in genebank-preserved germplasm is nowadays crucial. Collecting such information, might help researchers and breeders developing new experimental design, providing an overview of the state of the art of the most advanced approaches for the identification of the genetic elements laying behind resilience. Building this starting point, we aim at providing a useful tool for tackling the global agricultural resilience goals in these crops.

Analysis of codon usage patterns in 48 Aconitum species

BACKGROUND

The Aconitum genus is a crucial member of the Ranunculaceae family. There are 350 Aconitum species worldwide, with about 170 species found in China. These species are known for their various pharmacological effects and are commonly used to treat joint pain, cold abdominal pain, and other ailments. Codon usage bias (CUB) analysis contributes to evolutionary relationships and phylogeny. Based on protein-coding sequences (PCGs), we selected 48 species of Aconitum for CUB analysis.

RESULTS

The results revealed that Aconitum species had less than 50% GC content. Furthermore, the distribution of GC content was irregular and followed a trend of GC1 > GC2 > GC3, indicating a bias towards A/T bases. The relative synonymous codon usage (RSCU) heat map revealed the presence of conservative codons with slight variations within the genus. The effective number of codons (ENC)-Plot and the parity rule 2 (PR2)-bias plot analysis indicate that natural selection is the primary factor influencing the variation in codon usage. As a result, we screened various optimal codons and found that A/T bases were preferred as the last codon. Furthermore, our Maximum Likelihood (ML) analysis based on PCGs among 48 Aconitum species yielded results consistent with those obtained from complete chloroplast (cp.) genome data. This suggests that analyzing mutation in PCGs is an efficient method for demonstrating the phylogeny of species at the genus level.

CONCLUSIONS

The CUB analysis of 48 species of Aconitum was mainly influenced by natural selection. This study reveals the CUB pattern of Aconitum and lays the foundation for future genetic modification and phylogenetic analyses.

Construction of a high-density genetic map for faba bean (Vicia faba L.) and quantitative trait loci mapping of seed-related traits

Faba bean (Vicia faba L.) is a valuable legume crop and data on its seed-related traits is required for yield and quality improvements. However, basic research on faba bean is lagging compared to that of other major crops. In this study, an F₂ faba bean population, including 121 plants derived from the cross WY7×TCX7, was genotyped using the Faba_bean_130 K targeted next-generation sequencing genotyping platform. The data were used to construct the first ultra-dense faba bean genetic map consisting of 12,023 single nucleotide polymorphisms markers covering 1,182.65 cM with an average distance of 0.098 cM. The map consisted of 6 linkage groups, which is consistent with the 6 faba bean chromosome pairs. A total of 65 quantitative trait loci (QTL) for seed-related traits were identified (3 for 100-seed weight, 28 for seed shape, 12 for seed coat color, and 22 for nutritional quality). Furthermore, 333 candidate genes that are likely to participate in the regulation of seed-related traits were also identified. Our research findings can provide a basis for future faba bean marker-assisted breeding and be helpful to further modify and improve the reference genome.

Fine mapping and gene silencing pinpoint Capana10g002229 as a strong candidate gene regulating the deciduous character of ripe pepper fruit (Capsicum spp.)

The pepper S locus, which controls the deciduous character of ripe fruit, was first fine mapped into an interval with a physical length of ~ 38.03 kb on chromosome P10. Capana10g002229, encoding a polygalacturonase, was proposed as a strong candidate gene based on sequence comparison, expression pattern analysis and virus-induced gene silencing (VIGS). The deciduous character of ripe fruit, which is controlled by the dominant S locus, is a domesticated trait with potential value in the pepper processing industry (Capsicum spp.). However, the gene associated with the S locus has not been identified. Here, one major QTL designated S10.1 was detected by using the F₂ population (n = 155) derived from BA3 (Capsicum annuum) × YNXML (Capsicum frutescens) and was further verified in an intraspecific backcross population (n = 254) derived from the cross between BB3 (C. annuum) and its wild relative Chiltepin (C. annuum var. glabriusculum) with BB3 as the recurrent parent. Then, a large BC₁F₂ population derived from the self-pollination of BB3 × (BB3 × Chiltepin) individuals and comprising 4217 individuals was used to screen the recombinants, and the S locus was ultimately delimited into a 38.03-kb region on chromosome P10 harbouring four annotated genes. Capana10g002229, encoding a polygalacturonase (PG), was proposed as the best candidate gene for S based on sequence comparison and expression pattern analyses. Downregulation of Capana10g002229 in fruits through VIGS significantly delayed fruit softening and abscission from the fruit-receptacle junction. Taken together, the results show that Capana10g002229 could be regarded as a strong candidate gene associated with the S locus in pepper. These findings not only lay a foundation for deciphering the molecular mechanisms underlying pepper domestication but also provide a strategy for genetic improvement of the deciduous character of ripe fruit using a marker-assisted selection approach.

Genetic diversity between local landraces and current breeding lines of pepper in China

Based on 22 qualitative traits, 13 quantitative traits, and 27 molecular markers (26 SSR and 1 InDel), in the current study we compared the diversity and population structure of 94 local landraces and 85 current breeding lines of pepper in China. The results showed that the Shannon Diversity indices of 9 qualitative traits and 8 quantitative traits in current breeding lines were greater than those of landraces, of which 11 were fruit organ-related traits. Compared with current breeding lines, the mean values of Gene Diversity index and Polymorphism Information content of local landraces were higher by 0.08 and 0.09, respectively. Population structure and phylogenetic tree analysis showed that the 179 germplasm resources could be divided into two taxa, dominated by local landraces and current breeding lines, respectively. The above results indicated that the diversity of quantitative traits of current breeding lines were higher than that of local landraces, especially traits related to fruit organs, but the genetic diversity based on molecular markers was lower than that of local landraces. Therefore, in the future breeding process, we should not only focus on the selection of target traits, but also strengthen the background selection based on molecular markers. Moreover, the genetic information of other domesticated species and wild species will be transferred to the breeding lines through interspecific crosses to expand the genetic background of the breeding material.

A 163-bp insertion in the Capana10g000198 encoding a MYB transcription factor causes male sterility in pepper (Capsicum annuum L.)

Male sterility provides an efficient approach for commercial exploitation of heterosis. Despite more than 20 genic male sterile (GMS) mutants documented in pepper (Capsicum annuum L.), only two causal genes have been successfully identified. Here, a novel spontaneous recessive GMS mutant, designated msc-3, is identified and characterized at both phenotypic and histological levels. Pollen abortion of msc-3 mutant may be due to the delayed tapetum degradation, leading to the non-degeneration of tetrads callosic wall. Then, a modified MutMap method and molecular marker linkage analysis were employed to fine mapping the msc-3 locus, which was delimited to the ~139.91-kb region harboring 10 annotated genes. Gene expression and structure variation analyses indicate the Capana10g000198, encoding a R2R3-MYB transcription factor, is the best candidate gene for the msc-3 locus. Expression profiling analysis shows the Capana10g000198 is an anther-specific gene, and a 163-bp insertion in the Capana10g000198 is highly correlated with the male sterile (MS) phenotype. Additionally, downregulation of Capana10g000198 in male fertile plants through virus-induced gene silencing resulted in male sterility. Finally, possible regulatory relationships of the msc-3 gene with the other two reported pepper GMS genes, msc-1 and msc-2, have been studied, and comparative transcriptome analysis reveals the expression of 16 GMS homologs are significantly downregulated in the MS anthers. Overall, our results reveal that Capana10g000198 is the causal gene underlying the msc-3 locus, providing important theoretical clues and basis for further in-depth study on the regulatory mechanisms of pollen development in pepper.

Challenges in quantifying genome erosion for conservation

Massive defaunation and high extinction rates have become characteristic of the Anthropocene. Genetic effects of population decline can lead populations into an extinction vortex, where declining populations show lower genetic fitness, in turn leading to lower populations still. The lower genetic fitness in a declining population due to a shrinking gene pool is known as genetic erosion. Three different types of genetic erosion are highlighted in this review: overall homozygosity, genetic load and runs of homozygosity (ROH), which are indicative of inbreeding. The ability to quantify genetic erosion could be a very helpful tool for conservationists, as it can provide them with an objective, quantifiable measure to use in the assessment of species at risk of extinction. The link between conservation status and genetic erosion should become more apparent. Currently, no clear correlation can be observed between the current conservation status and genetic erosion. However, the high quantities of genetic erosion in wild populations, especially in those species dealing with habitat fragmentation and habitat decline, may be early signs of deteriorating populations. Whole genome sequencing data is the way forward to quantify genetic erosion. Extra screening steps for genetic load and hybridization can be included, since they could potentially have great impact on population fitness. This way, the information yielded from genetic sequence data can provide conservationists with an objective genetic method in the assessment of species at risk of extinction. However, the great complexity of genome erosion quantification asks for consensus and bridging science and its applications, which remains challenging.

Determining Genetic Diversity and Population Structure of Common Bean (Phaseolus vulgaris L.) Landraces from Türkiye Using SSR Markers

Assessment of genetic diversity among different varieties helps to improve desired characteristics of crops, including disease resistance, early maturity, high yield, and resistance to drought. Molecular markers are one of the most effective tools for discovering genetic diversity that can increase reproductive efficiency. Simple sequence repeats (SSRs), which are codominant markers, are preferred for the determination of genetic diversity because they are highly polymorphic, multi-allelic, highly reproducible, and have good genome coverage. This study aimed to determine the genetic diversity of 40 common bean (Phaseolus vulgaris L.) landraces collected from the Ispir district located in the Northeast Anatolia region of Türkiye and five commercial varieties using SSR markers. The Twenty-seven SSR markers produced a total of 142 polymorphic bands, ranging from 2 (GATS91 and PVTT001) to 12 (BM153) alleles per marker, with an average number of 5.26 alleles. The gene diversity per marker varied between 0.37 and 0.87 for BM053 and BM153 markers, respectively. When heterozygous individuals are calculated proportional to the population, the heterozygosity ranged from 0.00 to 1.00, with an average of 0.30. The expected heterozygosity of the SSR locus ranged from 0.37 (BM053) to 0.88 (BM153), with an average of 0.69. Nei's gene diversity scored an average of 0.69. The polymorphic information content (PIC) values of SSR markers varied from 0.33 (BM053) to 0.86 (BM153), with an average of 0.63 per locus. The greatest genetic distance (0.83) was between lines 49, 50, 53, and cultivar Karacaşehir-90, while the shortest (0.08) was between lines 6 and 26. In cluster analysis using Nei's genetic distance, 45 common bean genotypes were divided into three groups and very little relationship was found between the genotypes and the geographical distances. In genetic structure analysis, three subgroups were formed, including local landraces and commercial varieties. The result confirmed that the rich diversity existing in Ispir bean landraces could be used as a genetic resource in designing breeding programs and may also contribute to Türkiye bean breeding programs.

Fine mapping and candidate gene analysis of the up locus determining fruit orientation in pepper (Capsicum spp.)

The up locus determining fruit orientation was fine-mapped into a region with a physical length of ~169.51 kb on chromosome P12 in pepper. Capana12g000958, encoding a developmentally regulated G protein 2, was proposed as the strongest candidate via sequence comparison and expression analysis. Fruit orientation is an important horticultural and domesticated trait, which is controlled by a single semi-dominant gene (up) in pepper. However, the gene underlying up locus has not yet been identified. In this study, the previously detected major QTL UP12.1 was firstly verified using a backcross population (n = 225) stem from the cross of BB3 (C. annuum) and its wild relative Chiltepin (C. annuum var. glabriusculum) using BB3 as the recurrent parent. Then, a large BC₁F₂ population (n = 1827) was used for recombinant screening to delimit the up locus into an interval with ~ 169.51 kb in length. Sequence comparison and expression analysis suggested that Capana12g000958, encoding a developmentally regulated G protein 2, was the most likely candidate gene for the up locus. There is no difference within the coding sequences of Capana12g000958 between BB3 and Chiltepin, while a SNP in the upstream of Capana12g000958 showed a complete correlation with the fruit orientation among a panel of 40 diverse pepper inbred lines. These findings will form a basis for gene isolation and reveal of genetic mechanism underlying the fruit orientation domestication in pepper.

Fine Mapping and Candidate Gene Identification for the CapUp Locus Controlling Fruit Orientation in Pepper (Capsicum spp.)

The orientation of fruits is a distinguishing morphological feature of pepper (Capsicum spp.) varieties. The pendent (downward curved) growth of the fruit stalks, known as pedicels, is highly correlated with fruit weight and pedicel length. A previous genetic analysis revealed that the pendent fruit orientation is governed by a dominant gene, and incomplete inheritance is also observed in some Capsicum accessions. To identify and localize this gene, a single quantitative trait locus (QTL) analysis was performed on one F₂ and two recombinant inbred line (RIL) populations, and a genome-wide association study (GWAS) was performed using a core collection. Common QTL regions associated with fruit orientation were detected on chromosome 12. A total of 187,966 SNPs were identified in a genotyping-by-sequencing (GBS) for GWAS analysis of 196 Capsicum annuum, 25 Capsicum baccatum, 21 Capsicum chinense, and 14 Capsicum frutescens accessions, representing the germplasm collection of South Korea. The results of these analyses enabled us to narrow down the CapUp region of interest to 200-250 Mbp on chromosome 12. Seven candidate genes were found to be located between two markers that were completely cosegregated with the fruit orientation phenotype. The findings and markers developed in this study will be helpful for additional understanding of pepper fruit development and breeding for fruit orientation.

Single nucleotide polymorphisms reveal genetic diversity in New Mexican chile peppers (Capsicum spp.)

BACKGROUND

Chile peppers (Capsicum spp.) are among the most important horticultural crops in the world due to their number of uses. They are considered a major cultural and economic crop in the state of New Mexico in the United States. Evaluating genetic diversity in current New Mexican germplasm would facilitate genetic improvement for different traits. This study assessed genetic diversity, population structure, and linkage disequilibrium (LD) among 165 chile pepper genotypes using single nucleotide polymorphism (SNP) markers derived from genotyping-by-sequencing (GBS).

RESULTS

A GBS approach identified 66,750 high-quality SNP markers with known map positions distributed across the 12 chromosomes of Capsicum. Principal components analysis revealed four distinct clusters based on species. Neighbor-joining phylogenetic analysis among New Mexico State University (NMSU) chile pepper cultivars showed two main clusters, where the C. annuum genotypes grouped together based on fruit or pod type. A Bayesian clustering approach for the Capsicum population inferred K = 2 as the optimal number of clusters, where the C. chinense and C. frutescens grouped in a single cluster. Analysis of molecular variance revealed majority of variation to be between the Capsicum species (76.08 %). Extensive LD decay (~ 5.59 Mb) across the whole Capsicum population was observed, demonstrating that a lower number of markers would be required for implementing genome wide association studies for different traits in New Mexican type chile peppers. Tajima's D values demonstrated positive selection, population bottleneck, and balancing selection for the New Mexico Capsicum population. Genetic diversity for the New Mexican chile peppers was relatively low, indicating the need to introduce new alleles in the breeding program to broaden the genetic base of current germplasm.

CONCLUSIONS

Genetic diversity among New Mexican chile peppers was evaluated using GBS-derived SNP markers and genetic relatedness on the species level was observed. Introducing novel alleles from other breeding programs or from wild species could help increase diversity in current germplasm. We present valuable information for future association mapping and genomic selection for different traits for New Mexican chile peppers for genetic improvement through marker-assisted breeding.

Population structure, genetic diversity and genomic selection signatures among a Brazilian common bean germplasm

Brazil is the world's largest producer of common bean. Knowledge of the genetic diversity and relatedness of accessions adapted to Brazilian conditions is of great importance for the conservation of germplasm and for directing breeding programs aimed at the development of new cultivars. In this context, the objective of this study was to analyze the genetic diversity, population structure, and linkage disequilibrium (LD) of a diversity panel consisting of 219 common bean accessions, most of which belonging to the Mesoamerican gene pool. Genotyping by sequencing (GBS) of these accessions allowed the identification of 49,817 SNPs with minor allele frequency > 0.05. Of these, 17,149 and 12,876 were exclusive to the Mesoamerican and Andean pools, respectively, and 11,805 SNPs could differentiate the two gene pools. Further the separation according to the gene pool, bayesian analysis of the population structure showed a subdivision of the Mesoamerican accessions based on the origin and color of the seed tegument. LD analysis revealed the occurrence of long linkage blocks and low LD decay with physical distance between SNPs (LD half decay in 249 kb, corrected for population structure and relatedness). The GBS technique could effectively characterize the Brazilian common bean germplasms, and the diversity panel used in this study may be of great use in future genome-wide association studies.

Fine mapping of restorer-of-fertility gene based on high-density genetic mapping and collinearity analysis in pepper (Capsicum annuum L.)

The pepper restorer-of-fertility (CaRf) gene was fine mapped based on conjoint analysis of recombinants and collinearity between the two pepper reference genomes. Capana06g003028, encoding an Rf-like PPR protein, was proposed as the strongest candidate for pepper CaRf based on sequence comparison and expression analysis. The cytoplasmic male sterility (CMS)/restorer-of-fertility (Rf) system not only provides an excellent model to dissect genetic interactions between the mitochondria and nucleus but also plays a vital role in high-efficiency hybrid seed production in crops including pepper (Capsicum spp.). Although two important CMS candidate genes, orf507 and Ψatp6-2, have previously been suggested, the pepper Rf gene (CaRf) has not yet been isolated. In this study, a high-density genetic map comprising 7566 SNP markers in 1944 bins was first constructed with the array genotyping results from 317 F₂ individuals. Based on this map, the CaRf gene was preliminarily mapped to a region of 1.15 Mb in length at the end of chromosome P6. Then, by means of a conjoint analysis of recombinants and collinearity between the two pepper reference genomes, an important candidate interval with 270.10 kb in length was delimited for CaRf. Finally, Capana06g003028, which encodes an Rf-like PPR protein, was proposed as the strongest candidate for CaRf based on sequence analysis and expression characteristics in sterile and fertile plants. The high-density genetic map and fine mapping results provided here lay a foundation for the application of molecular breeding, as well as cloning and functional analysis of CaRf, in pepper.

Uncovering Candidate Genes Controlling Major Fruit-Related Traits in Pepper via Genotype-by-Sequencing Based QTL Mapping and Genome-Wide Association Study

All modern pepper accessions are products of the domestication of wild Capsicum species. However, due to the limited availability of genome-wide association study (GWAS) data and selection signatures for various traits, domestication-related genes have not been identified in pepper. Here, to address this problem, we obtained data for major fruit-related domestication traits (fruit length, width, weight, pericarp thickness, and fruit position) using a highly diverse panel of 351 pepper accessions representing the worldwide Capsicum germplasm. Using a genotype-by-sequencing (GBS) method, we developed 187,966 genome-wide high-quality SNP markers across 230 C. annuum accessions. Linkage disequilibrium (LD) analysis revealed that the average length of the LD blocks was 149 kb. Using GWAS, we identified 111 genes that were linked to 64 significant LD blocks. We cross-validated the GWAS results using 17 fruit-related QTLs and identified 16 causal genes thought to be associated with fruit morphology-related domestication traits, with molecular functions such as cell division and expansion. The significant LD blocks and candidate genes identified in this study provide unique molecular footprints for deciphering the domestication history of Capsicum. Further functional validation of these candidate genes should accelerate the cloning of genes for major fruit-related traits in pepper.

High-density genetic map construction and QTL mapping of first flower node in pepper (Capsicum annuum L.)

BACKGROUND

First flower node (FFN) is an important trait for evaluating fruit earliness in pepper (Capsicum annuum L.). The trait is controlled by quantitative trait loci (QTL); however, studies have been limited on QTL mapping and genes contributing to the trait.

RESULTS

In this study, we developed a high density genetic map using specific-locus amplified fragment sequencing (SLAF-seq), a high-throughput strategy for de novo single nucleotide polymorphism discovery, based on 146 recombinant inbred lines (RILs) derived from an intraspecific cross between PM702 and FS871. The map contained 9328 SLAF markers on 12 linkage groups (LGs), and spanned a total genetic distance of 2009.69 centimorgan (cM) with an average distance of 0.22 cM. The sequencing depth for the map was 72.39-fold in the male parent, 57.04-fold in the female parent, and 15.65-fold in offspring. Using the genetic map, two major QTLs, named Ffn2.1 and Ffn2.2, identified on LG02 were strongly associated with FFN, with a phenotypic variance explanation of 28.62 and 19.56%, respectively. On the basis of the current annotation of C. annuum cv. Criollo de Morelos (CM334), 59 candidate genes were found within the Ffn2.1 and Ffn2.2 region, but only 3 of 59 genes were differentially expressed according to the RNA-seq results. Eventually we identified one gene associated with the FFN based on the function through GO, KEGG, and Swiss-prot analysis.

CONCLUSIONS

Our research showed that the construction of high-density genetic map using SLAF-seq is a valuable tool for fine QTL mapping. The map we constructed is by far the most saturated complete genetic map of pepper, and using it we conducted fine QTL mapping for the important trait, FFN. QTLs and candidate genes obtained in this study lay a good foundation for the further research on FFN-related genes and other genetic applications in pepper.

Mapping and identification of CsUp, a gene encoding an Auxilin-like protein, as a putative candidate gene for the upward-pedicel mutation (up) in cucumber

BACKGROUND

Pedicel orientation can affect the female flower orientation and seed yield in cucumber. A spontaneous mutant possessing upward growth of pedicels was identified in the wild type inbred strain 9930 and named upward-pedicel (up). The morphological and genetic analyses of up were performed in this study. In order to clone the up gene, 933 F₂ individuals and 524 BC₁ individuals derived from C-8-6 (WT) and up were used for map-based cloning.

RESULTS

up was mapped to a 35.2 kb physical interval on chromosome 1, which contains three predicted genes. Sequencing analysis revealed that a 5-bp deletion was found in the second exon of Csa1G535800, and it led to a frameshift mutation resulting in a premature stop codon. The candidate gene of CsUp (Csa1G535800) was further confirmed via genomic and cDNA sequencing in biparental and natural cucumber populations. Sequencing data showed that a 4-bp deletion was found in the sixth exon of Csa1G535800 in CGN19839, another inbred line, and there was also a mutation of an amino acid in Csa1G535800 that could contribute to the upward growth of pedicels in CGN19839. Moreover, it was found that Csa1G535800 exhibited strong expression in the pedicel of WT, suggesting its important role in development of pedicel orientation. Thus, Csa1G535800 was considered to be the candidate gene of CsUp.

CONCLUSIONS

CsUp encodes an Auxilin-like protein and controls pedicel orientation in cucumber. The identification of CsUp may help us to understand the mechanism of pedicel orientation development and allow for investigation of novel functions of Auxilin-like proteins in cucumber.

A MYB transcription factor is a candidate to control pungency in Capsicum annuum

Identification of a novel pungency-controlling gene Pun3, which acts as a master regulator of capsaicinoid biosynthetic genes in Capsicum annuum. Capsaicinoid is a unique compound that gives hot peppers (Capsicum spp.) their spicy taste. The Pun1 and Pun2 loci are known to control pungency in Capsicum species. Whereas Pun1 encodes an acyltransferase, the identity of Pun2 is currently unknown. Here, we used recombinant inbred lines and F₂ plants derived from a cross between the non-pungent C. annuum accession 'YCM334' and the pungent C. annuum cultivar 'Tean' to identify a novel non-pungency locus. Inheritance studies showed that non-pungency in C. annuum 'YCM334' is controlled by a single recessive gene, which we named Pun3. Using a high-density SNP map derived from genotyping-by-sequencing, Pun3 was mapped to chromosome 7. By comparing physical information about the Pun3 region in the C. annuum 'Zunla-1' and C. chinense 'PI159236' reference genomes, we identified candidate genes in this target region. One cDNA sequence from 'PI159236' was homologous to an unannotated gene in 'Zunla-1.' This sequence was also homologous to CaMYB31, which is expressed only in 'Tean' and harbors one stop codon in the non-pungent accession 'YCM334.' RNA-Seq analysis showed that major structural genes in the capsaicinoid biosynthetic pathway were significantly downregulated in 'YCM334' compared to pungent pepper. Therefore, CaMYB31 is a candidate gene for Pun3, which may act as a master regulator of capsaicinoid biosynthetic genes in pepper.

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.

Genetic diversity, population structure, and relationships in a collection of pepper (Capsicum spp.) landraces from the Spanish centre of diversity revealed by genotyping-by-sequencing (GBS)

Pepper (Capsicum spp.) is one of the most important vegetable crops; however, pepper genomic studies lag behind those of other important Solanaceae. Here we present the results of a high-throughput genotyping-by-sequencing (GBS) study of a collection of 190 Capsicum spp. accessions, including 183 of five cultivated species (C. annuum, C. chinense, C. frutescens, C. baccatum, and C. pubescens) and seven of the wild form C. annuum var. glabriusculum. Sequencing generated 6,766,231 high-quality read tags, of which 40.7% were successfully aligned to the reference genome. SNP calling yielded 4083 highly informative segregating SNPs. Genetic diversity and relationships of a subset of 148 accessions, of which a complete passport information was available, was studied using principal components analysis (PCA), discriminant analysis of principal components (DAPC), and phylogeny approaches. C. annuum, C. baccatum, and C. chinense were successfully separated by all methods. Our population was divided into seven clusters by DAPC, where C. frutescens accessions were clustered together with C. chinense. C. annuum var. glabriusculum accessions were spread into two distinct genetic pools, while European accessions were admixed and closely related. Separation of accessions was mainly associated to differences in fruit characteristics and origin. Phylogeny studies showed a close relation between Spanish and Mexican accessions, supporting the hypothesis that the first arose from a main genetic flow from the latter. Tajima's D statistic values were consistent with positive selection in the C. annuum clusters, possibly related to domestication or selection towards traits of interest. This work provides comprehensive and relevant information on the origin and relationships of Spanish landraces and for future association mapping studies in pepper.

Genetic mapping of anthocyanin accumulation-related genes in pepper fruits using a combination of SLAF-seq and BSA

Anthocyanins have significant functions in stress tolerance in pepper (Capsicum annuum L.) and also benefit human health. Nevertheless, the key structural genes and regulatory genes involved in anthocyanin accumulation in pepper fruits are still not well understood and fine mapped. For the present study, 383 F₂ plants from a cross between the green-fruited C. annuum line Z5 and the purple-fruited line Z6 was developed. Two separate bulked DNA pools were constructed with DNAs extracted from either 37 plants with high anthocyanin content or from 18 plants with no anthocyanin. A combination of specific-locus amplified fragment sequencing (SLAF-seq) and bulked segregant analysis (BSA) was used to identify candidates for regions associated with anthocyanin accumulation. We identified a total of 127,004 high-quality single nucleotide polymorphism (SNP) markers, and detected 1674 high-quality SNP markers associated with anthocyanin accumulation. Three candidate anthocyanin-associated regions including the intervals from 12.48 to 20.00 Mb, from 54.67 to 56.59 Mb, and from 192.17 to 196.82 Mb were identified within a 14.10-Mb interval on chromosome 10 containing 126 candidate genes. Based on their annotations, we identified 12 candidate genes that are predicted to be associated with anthocyanin expression. The present results provide an efficient strategy for genetic mapping of and valuable insights into the genetic mechanisms of anthocyanin accumulation in pepper fruit, and allow us to clone and functionally analyze the genes that influence anthocyanin accumulation in this species.

Candidate genes for first flower node identified in pepper using combined SLAF-seq and BSA

First flower node (FFN) is an important trait for evaluating fruit earliness in pepper (Capsicum annuum L.), but the genetic mechanisms that control FFN are still poorly understood. In the present study, we developed 249 F₂ plants derived from an intraspecific cross between the inbred pepper lines Z4 and Z5. Thirty plants with the highest FFN and 30 plants with the lowest FFN were chosen and their DNAs were pooled according to phenotype to construct two bulked DNA pools. Specific-locus amplified fragment sequencing (SLAF-seq) was combined with bulked segregant analysis (BSA) to identify candidate regions related to FFN. According to our genetic analysis, the FFN trait is quantitatively inherited. A total of 106,848 high-quality single nucleotide polymorphism (SNP) markers were obtained, and 393 high-quality SNP markers associated with FFN were detected. Ten candidate regions within an interval of 3.98 Mb on chromosome 12 harboring 23 candidate genes were identified as closely correlated with FFN. Five genes (CA12g15130, CA12g15160, CA12g15370, CA12g15360, and CA12g15390) are predicted based on their annotations to be associated with expression of the FFN trait. The present study demonstrates an efficient genetic mapping strategy and lays a good foundation for molecular marker-assisted breeding using SNP markers linked to FFN and for cloning and functional analysis of the key genes controlling FFN.