Hea-Young Lee
Na-Young Ro2
Abhinandan Patil
Jin-Kyung Kwon
Byoung-Cheorl Kang- 1Department of Plant Science, Plant Genomics and Breeding Institute and Vegetable Breeding Research Center, College of Agriculture and Life Sciences, Seoul National University, Seoul, South Korea
- 2National Academy of Agricultural Science, National Agrobiodiversity Center, Rural Development Administration, Jeonju, South Korea
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.
Introduction
Pepper (Capsicum species), like other Solanaceae family members including tomato and potato, is a New World genus with a primary center of diversity in Bolivia and Peru (Nee et al., 2006). Capsicum comprises more than 30 species, and the domestication of five of these species in the Americas, including the economically important plants C. annuum, C. baccatum, C. chinense, C. frutescens, and C. pubescens, dates back to 6,000 BC (Moscone et al., 2007; Paran and Van Der Knaap, 2007; Cheng et al., 2016). Peppers are referred to as capsicum, pimento, sweet pepper, red pepper, cayenne pepper, bird’s eye pepper, jalapenos, or habaneros based on fruit shape and pungency (Moscone et al., 2007; Babu et al., 2011), and have various uses as vegetables, seasonings, ornamental plants, and medicinal crops. The easy cultivation of pepper has led to their widespread use worldwide, especially in tropical regions (Moscone et al., 2007; Babu et al., 2011). The majority of wild forms of Capsicum spp. display perennial herbaceous growth, with a small, erect, deciduous growth habit and red, pungent, soft-fleshed fruits (Paran and Van Der Knaap, 2007).
Among Solanaceae species, domestication-related traits have been described for tomato (Bai and Lindhout, 2007; Giovannoni 2018), potato (Li et al., 2018), and eggplant (Doganlar et al., 2002; Meyer et al., 2012). These traits are generally referred to as “domestication syndrome” because they can be used to distinguish cultivated crops from their progenitors (Doebley et al., 2006). The domestication syndrome traits are not fully elucidated for pepper. During domestication, Capsicum spp. might have been selected for fruit morphology and pungency (Babu et al., 2011; Che and Zhang, 2019). Other pepper domestication traits include a non-deciduous habit, fruit that remains on the plant until harvest, and pendent fruit orientation (Kaiser, 1935 Paran and Van Der Knaap, 2007). However, underlying genes are largely known.
Genetic and genomic analyses of cultivated crops and wild relatives have provided evidence for domestication by revealing selection footprints in the key genes controlling domestication traits (Zeder et al., 2006; Stitzer and Ross-Ibarra, 2018). Recent genetic and archaeological studies have revealed the spatiotemporal origins and processes underlying the domestication of these traits and have allowed domestication traits to be divided into two types based on the underlying genes. Some domestication traits are controlled by genes called ‘domestication genes’ that were subjected to early selection of major-effect QTLs, while other traits are controlled by genes that were selected later to produce diversified, improved crops; these genes are called ‘improvement genes’ (Pickersgill, 2007). Wang and Bosland (2006) published a comprehensive summary of genetic studies on Capsicum genes performed from 1912 to 2006 that lists 292 genes for various traits in pepper, including morphological and physiological traits, male sterility, and resistance to nematodes, diseases, and herbicides.
Most traditional QTL analyses in pepper have focused on fruit morphology-related traits. These studies have involved low-throughput genotyping or focused only on identifying the genes governing these traits. For instance, genetic mapping studies identified QTLs for fruit length, fruit width (FWd), fruit weight (FWg), pericarp thickness (PT), and fruit position (FP). Among these, fs2.1, FrSHP2.1, and fs3.1 are the major QTLs for fruit shape; these QTLs are located on chromosomes 2 and 3 (Chaim et al., 2001; Chaim et al., 2003; Rao et al., 2003; Zygier et al., 2005; Barchi et al., 2009; Borovsky and Paran, 2011; Mimura et al., 2012; Hill et al., 2017; Chunthawodtiporn et al., 2018).
By contrast, due to advancements in next-generation sequencing (NGS) techniques and the availability of newer populations of tomato, six representative gene families were identified to control fruit size in this crop, including the Cell Number Regulator (CNR), Cytochrome P450 A78 class (CYP78A), IQ domain, Ovate Family Protein (OFP), YABBY, and WOX gene families (Monforte et al., 2014; Lin et al., 2014; Sacco et al., 2015; Soyk et al., 2017). Candidate genes belonging to these families such as CNR, SlKLUH, SUN, OVATE, FAS, and LC have been cloned, and their roles in regulating fruit elongation, locule number, and fruit shape have been well characterized (Xiao et al., 2008; Guo and Simmons, 2011; Rodriguez et al., 2011; Chakrabarti et al., 2013). These findings from tomato were successfully utilized for QTL mapping and downstream gene analysis in pepper, shedding light on the complex genetic architecture and genomic regions that govern these quantitatively inherited traits (Ramchiary et al., 2014; Wang et al., 2015; Chunthawodtiporn et al., 2018; Colonna et al., 2019).
Since the release of the first reference genome of pepper (Kim et al., 2014), genome-wide association study (GWAS) has been used to analyze only a few traits in pepper such as fruit weight (Nimmakayala et al., 2016a), capsaicinoid contents (Nimmakayala et al., 2016a; Han et al., 2018), peduncle length (Nimmakayala et al., 2016b), and fruit size and shape (Colonna et al., 2019) using diverse pepper germplasms. Combined QTL mapping and GWAS has been utilized to avoid identifying false-positive QTLs or associations for major fruit-related traits in pepper.
The goals of the current study were to (1) explore the correlations among five important fruit-related traits in pepper and (2) determine the significant genetic regions or genes governing genetic variations in the major fruit-related traits with strong evidence for selection during domestication. We obtained high-quality SNPs via genotype-by-sequencing (GBS) and subjected them to GWAS. We obtained candidate genes underlying the QTLs detected by GWAS and characterized their functions, laying the foundation for further functional validation and cloning of candidate genes for major fruit-related traits in pepper.
Materials and Methods
Plant Materials
A collection of 351 Capsicum accessions known as the ‘pepper GWAS population’ was used for analysis, comprising four major domesticated species including C. annuum (230 accessions), C. baccatum (48 accessions), C. chinense (48 accessions), and C. frutescens (25 accessions) (Table S1). Among the accessions in the pepper GWAS population, 250 were previously selected as a core set representing the genetic diversity of more than 4,600 accessions from 97 countries (Lee et al., 2016).
Phenotypic Evaluation and Correlation Analysis
The pepper GWAS population was planted in a greenhouse at the Rural Development Administration (RDA)-Gene bank Jeonju, Republic of Korea (35°49′51.3” N, 127°03′47.1” E). Over a three-year period (2015–2017), six plants per accession were randomly planted, and the phenotypes of three plants per accession were evaluated. Five domestication traits were evaluated using a randomized block design, including four quantitative traits (fruit length [FL], width [FWd], weight [FWg], and pericarp thickness [PT]) and one qualitative trait (fruit position [FP]). All quantitative traits were measured using an electronic scale and a ruler. FP was scored as 1 to 3 (1 = erect, 2 = declining like a pendant, and 3 = intermediate). The correlation among the five traits was evaluated by Pearson correlation (r) analysis with SPSS software (IBM Corp. Released 2017. IBM SPSS Statistics for Windows, Version 25.0. Armonk, NY).
Genomic DNA Extraction and GBS
Genomic DNA was extracted from the samples using the CTAB method (Lee et al., 2017; Siddique et al., 2019; Solomon et al., 2019) and diluted to 80 ng/µl in distilled water. GBS libraries were constructed via double digestion with two sets of restriction enzymes (PstI/MseI and EcoRI/MseI) as previously described (Han et al., 2018; Siddique et al., 2019; Solomon et al., 2019). The digested DNA was ligated to adapters and amplified with ‘TA’ primers. The libraries were pooled in five tubes. The contents of the tubes were sequenced in separate lanes using the HiSeq 2000 platform (Illumina, San Diego, CA) at Macrogen (Seoul, Republic of Korea).
Reference-Based SNP Calling and Construction of the SNP Set
Raw 101-bp reads from the libraries were trimmed to a minimum length of 80 bp and filtered to a quality score >Q30. The filtered reads were aligned to the C. annuum ‘CM334’ reference genome v.1.6, http://peppergenome.snu.ac.kr (Kim et al., 2017) using the Burrows-Wheeler Aligner program v.0.7.12 (Li and Durbin, 2010). For SNP calling and filtering, the GATK Unified Genotyper v.3.3-0 was used (Depristo et al., 2011). The SNP set was constructed using three steps: pre-filtering, imputation, and major filtering. First, pre-imputed SNPs were filtered to removed mono and tri-allelic SNP types and SNPs with a call rate >0.1. After pre-imputation filtering, SNPs with missing data were imputed using the FILLIN method in TASSEL (Bradbury et al., 2007). To obtain SNPs of suitable quality, hapSize was applied to obtain sequences ranging from 100 to 8,000 with two minSites (25, 50) and two minPres (250, 500). The final selected imputation option was dependent on the best option of regression (R2) values and the imputed ratio of minor and major alleles. Finally, the major filtering step was performed under the following conditions: minor allele frequency >0.05, SNP coverage >0.6, and inbreeding coefficient (IF) >0.8.
Population Structure (Q) and Linkage Disequilibrium (LD) Estimations
To identify population stratification, principal component analysis (PCA) was performed using the ‘pcaMethods’ library in R software (Stacklies et al., 2007). The values of each PC were used as variables in the GWAS. The LD block of the GWAS population was estimated using PLINK v.1.9 (Chang et al., 2015) with the following settings: ‘–no-parents –no-sex –blocks no-pheno-req no-small-max-span –blocks-inform-frac 0.95 –blocks-max-kb 2000 –blocks-min-maf 0.05 –blocks-recomb-highci 0.9 –blocks-strong-highci 0.98 –blocks-strong-lowci 0.7’. The calculated LD block intervals were used to search for candidate genes for specific traits.
Genome-Wide Association Study (GWAS) and Candidate Gene Identification
GWAS based on the compressed mixed linear model (CMLM) was conducted using the R package of Genomic Association and Prediction Integrated Tool (GAPIT) (Lipka et al., 2012) with the forward model selection using the Bayesian information criterion (BIC). The significance threshold −log10 P-value of the GWAS was determined using Bonferroni (1936) correction (FDR P-value < 0.05) based on the number of independent SNPs in a population. The candidate genes inside LD regions with significant SNPs were investigated. Gene prediction was performed based on gff file v.2.0 of the CM334 v.1.6 reference genome (http://peppergenome.snu.ac.kr), and the function of each gene was predicted using Blast2GO (Götz et al., 2008) based on deduced protein sequences. To detect the physical positions of previous QTLs or orthologous genes from other species, BLAST searches of the pepper genome from the NCBI database (https://www.ncbi.nlm.nih.gov) were performed. To predict the molecular and biological functions of genes, the NCBI, Solanaceae (https://solgenomics.net), and Arabidopsis databases (https://www.arabidopsis.org) were used.
QTL Mapping Using Recombinant Inbred Lines
To validate the GWAS results, recombinant inbred lines (RILs) were used for QTL mapping as described by Han et al. (2016). All information about the plant materials and phenotypes in this study were described in Han et al. (2016), but the genotyping results were altered using a more recent version of the reference genome (C. annuum ‘CM334’ v.1.6, http://peppergenome.snu.ac.kr). In brief, 120 F7–F10 RILs derived from a cross between pungent C. annuum ‘Perennial’ and non-pungent C. annuum ‘Dempsey’ were grown for 3 years (2011, 2012, and 2014) in two locations: Anseong (2011 and 2012a) and Suwon, Korea (2012b and 2014). After SNP calling, sequencing reads were aligned to C. annuum ‘CM334’ v.1.6. A modified sliding window approach was used to investigate recombination breakpoints and to construct a bin map of the RILs. Bins were used as markers to construct a genetic map using the Carthagene program (De Givry et al., 2005) with default threshold values. All detailed options were adapted from Han et al. (2016) and Siddique et al. (2019). Of the 18 reported traits (Han et al., 2016), we utilized major four fruit domestication-related traits (FL, FWd, FWg, and FP) for QTL analysis. Composite interval mapping (CIM) was performed with Windows QTL Cartographer 2.5 (Wang et al., 2012). The phenotypic values of each trait in the respective years and locations were analyzed separately to detect QTLs. The log of odds (LOD) threshold was determined by performing 1,000 permutation tests with 5% probability (P) for each trait, and the proportion of phenotypic variation (R2) for each QTL was estimated. The 95% confidence interval was used to represent the location of each QTL.
Results
SNP Discovery and Population Structure (Q) of a Pepper GWAS Population
We aligned sequences derived from GBS to the C. annuum cv. CM334 reference genome v. 1.6, http://peppergenome.snu.ac.kr (Kim et al., 2017). GBS genotyping of 351 accessions (Table S1) with two sets of libraries constructed using double digestion with two sets of restriction enzymes generated 8,717,361 SNPs. The GBS generated data is available in National Agricultural Biotechnology Information Center (NABIC, https://nabic.rda.go.kr/, ID= NV-0630-000001) Trimming and filtering-out of SNPs with a quality score <30, call rate <10%, and mono or tri-allelic SNPs types resulted in 1,869,524 SNPs (Table S2). To avoid potential errors in the interpretation of the GWAS results, the missing genotypes were imputed using the FILLIN method in the TASSEL package. Accordingly, approximately 26% of genotypes were imputed to minor alleles, and 21 and 59% of genotypes were imputed to hetero and major alleles, respectively, with a regression (R2) value of 0.82. Using this imputed genotype, the final filtering step was performed under the following conditions: MAF >0.05, SNP coverage >0.6, and IF >0.8. This step resulted in a set of 507,713 high-quality SNPs, which were evenly distributed on 12 chromosomes (Figure 1A, Table S2). Each SNP marker generated from this SNP set was named according to its physical position in the pepper reference genome.
Figure 1 SNP distribution and population structure (Q) of the GWAS population. (A, B) SNP distribution across genotypes on 12 chromosomes within a 1 Mbp window size (A) SNP distribution across all accessions using 507,713 SNPs, (B) SNP distribution across all 230 C. annuum accessions using 187,966 SNPs). (C–E) Genetic distributions based on principal component analysis (C) Genetic distribution across all accessions using 507,713 SNPs, (D) Genetic distribution across 230 C. annuum accessions using 507,713 SNPs, (E) Genetic distribution across 230 C. annuum accessions using 187,966 SNPs).
As the genetic structure of a population can strongly affect the results of GWAS, we performed principal component analysis (PCA) to analyze population stratification. This analysis yielded four genetic clusters, with 22.26% (PC1) and 13.47% (PC2) of the genotypic variance in the first and second axes. Each cluster was well clustered by species, including C. annuum, C. baccatum, C. chinense, and C. frutescens. Although some accessions showed slight admixture, there were no conspicuous sub-clusters in the structure (Figure 1C).
Phenotypic Diversity of Major Fruit-Related Domestication Traits in the GWAS Population
We evaluated 351 Capsicum accessions from four species with maximum genetic diversity (Lee et al., 2016) for three years to assess the range of phenotype variation of the five fruit-related domestication traits (FL, FWd, FWg, PT, and FP). We detected a consecutive reduction in the mean values of the quantitative traits during the three years of the experiment, except for FWd, which had a higher value in 2016 (25.1 mm) compared to 2015 (24.3 mm) and 2017 (23.7 mm). While the maximum average FL value (72.4 mm) was obtained in 2015, the minimum value (66.9 mm) was obtained in 2017. Similarly, FWg and PT showed the highest average values in 2015 (22.9 g and 2.5 mm, respectively), followed by 2016 (19.9 g and 2.0 mm, respectively) (Table 1). Three species showed either erect (40%) or pendant (60%) FP, whereas C. frutescens showed all erect FP. Intermediate FP was also observed in all species except C. frutescens (Figure 2B). High broad-sense heritability (H2) values were recorded for FL (0.8), FWd (0.81), FWg (0.83), and PT (0.72) (Table 1).
Table 1 Phenotypic variation of four quantitative traits identified by GWAS population-based analysis over a three-year period.
Figure 2 Phenotype performance and correlations among five major fruit-related domestication traits. (A) Pearson correlation coefficients (r) among the five investigated traits. Numbers indicate the correlations between two traits. Red blocks indicate positive correlations. Asterisks (**) represent a significant difference at P-value 0.01. (B–F) morphological distributions of four domestication-related traits among species (B: FP, C: FL, D: FWd, E: FWg, F: PT). Except FP, each box depicts the upper and lower quantile, with the median represented by a horizontal solid line. Outliers are indicated by dots. Different letters indicate significant difference at P-value <0.05, as determined by one-way ANOVA with Scheffe multiple comparison post-hoc test.
All five fruit-related domestication traits showed significant positive correlations (P = 0.01). Specifically, highly strong positive correlations were detected between FWg and FWd (r = 0.91), followed by FWd and PT (r = 0.90) and FWg and PT (r = 0.88). Although FL had slightly lower positive correlations with these three traits (FWd, FWg, PT), it was the most highly correlated with FP (r = 0.46) compared to the three other traits (Figure 2A). As the GWAS population was clustered by species, we performed ANOVA, which validated the variability in the traits among species (Figures 2B–F, Tables S3 and S4). This analysis uncovered significant variation between species groups for the five traits, ranging from 9.63 to 22.68 (F, with p = 0.00), with a mean difference of 0.08 to 0.17 (η2), which also supported the differences among species (Table S4). Most of the traits showed the greatest mean values in the C. annuum accessions (FL: 78.8 mm, FWd: 27.6 mm, FWg: 26.8 g, PT: 2.4 mm), whereas the lowest mean values were detected in C. frutescens (Figures 2C–F).
Linkage Disequilibrium (LD) Pattern and GWAS of the C. annuum Cluster
The minimize the confounding effect of interspecies variation and the corresponding false-positive errors, among the entire population set used in the experiment, we selected the C. annuum cluster, as it contained a sufficient number of accessions with high levels of phenotypic and genotypic diversity without any interrelated population stratification (Figures 1D, E). Since the five fruit-related traits showed high broad-sense heritability (H2) values (Table 1), indicating that genetic factors were the major determinants of the observed phenotypic variability, we subjected all traits to GWAS using the C. annuum CM334 v.1.6, http://peppergenome.snu.ac.kr reference genome.
In the C. annuum cluster, 187,966 high-quality SNPs were filtered for use in GWAS following the criteria described above (Figure 1B, Table S2). Using this SNP set, we compared the common and unique patterns of genetic variation in adjacent marker pairs of each chromosome by performing LD analysis throughout the genomes of the C. annuum accessions. We identified 12,234 LD blocks, with an average of 1,020 per chromosome. The average block size was 149 kb, each containing an average of 14 SNPs (Table S5). The LD blocks were named based on their order on each chromosome.
We detected high variation for all fruit-related traits over the three years of analysis except for the qualitative trait (FP), which was observed for only one year. Common SNPs that were consistently correlated for at least two years of investigation exceeding the significance threshold (−log10 P >6.575) were used to describe our results (Figures S1 and S2, Table S6). Accordingly, a total of 178 common SNPs were identified, including 1 for FL, 148 for FWd, 28 for FWg, and 1 for PT. For FP, 52 significant SNPs from accessions with pendant, erect, and intermediate phenotypes were used for analysis (Table S6).
Of the 230 common SNPs, one SNP located on chromosome 4 (S04_227983120) was common both to FWg and PT. Furthermore, five SNPs on chromosome 9 (S09_100362495, S09_133144036, S09_136634514, S09_136634573, S09_143733895) and two SNPs on chromosome 12 (S12_14626660, S12_17471128) were detected for both FWd and FWg. Unlike these eight common SNPs, 222 SNPs were associated with 64 LD blocks distributed on chromosomes 2, 3, 4, 5, 6, 9, 10, and 12 (Figure S3, Table S7).
In detail, for marker–trait associations per year, eight significant SNPs were identified for FL; of these, three SNPs (S03_19218749, S03_19218759, and S03_19254384) on chromosome 3 were detected only in 2017. Two SNPs were identified on chromosome 4, including S04_211838587 only in 2016 and S04_211848210 in all three years of the study. Two and one additional SNP (S05_12080328, S07_214135504, and S11_72669050) each on chromosomes 5, 7, and 11 were identified only in 2015 and 2017, respectively (Figure S1). SNP S04_211848210, which is associated with LD block H04-0562 on chromosome 4 in the region between 211.8 Mb and 211.9 Mb (spanning an interval of 64,916 bp), was expected to be highly correlated with FL, as it was consistently detected throughout the experiment (Figures 3A, B). For FWd, we detected 281 significant SNPs throughout the experimental period, with an average –log10 P-value ranging from 6 to 10.57 (Figure S1). Most of the significant SNPs (98.9%) were located on chromosome 9 from 67.4 to 171 Mb and were linked to 26 LD blocks (Figure S3). In 2015, nine unique significant SNP positions (eight on chromosome 9 and one on chromosome 7) were detected. There were 123 unique SNPs in 2016, all on the middle and distal regions of chromosome 9. While 146 common SNPs were detected in 2015 and 2016, only two SNPs were identified on chromosome 12. Notably, two SNPs (S09_133144036 and S09_136634573) associated with H09-0745 and H09-0756 were consistently identified in all three years of the study (Figure 4, Table S6).
Figure 3 GWAS identifies significant SNPs and LD block regions containing genes controlling FL, PT, and FP in pepper. The left side of each Manhattan plot shows chromosome-wide associations. The most significant areas are indicated by arrows. The significant areas are shown in detail, with the three years of analysis represented by various shapes (circle: 2015, triangle: 2016, X: 2017). Under the Manhattan plot, LD blocks are represented by gray boxes. Candidate genes are indicated by red dotted lines, with the names and positions inside the LD block indicated by blue bars. (A, B): FL, (C–E): PT, (F–H): FP, (A) association with FL on chromosome 4. (B) Close-up view of the significant LD block regions (211.5–212.5 Mbp) on chromosome 4. (C) Association with PT on chromosome 4. (D) Close-up view of the significant LD block regions (227.5–228.5 Mbp) on chromosome 4. (E) Gene structure with DNA polymorphism. Below the gene structure, boxplots show PT based on allelic differences of significant SNP; the width of each box is proportional to the square root of the number of accessions. (F) Association with FP on chromosome 12. (G) Close-up view of the significant LD block regions (211.5–212.5 Mbp) on chromosome 12. (H) close-up view of the significant LD block regions (214–215 Mbp) on chromosome 12.
Figure 4 GWAS identifies significant LD block regions and correlated genes controlling FWd in pepper. The significant areas are shown in detail next to chromosome-wide Manhattan plots, in which all three years of analysis are represented by various shapes (circle: 2015, triangle: 2016, X: 2017). Under the Manhattan plot, LD blocks are represented by gray boxes. Candidate genes are represented by red dotted lines, with the names and positions indicated inside the LD blocks indicated by blue bars. (A, I): chromosome-wide Manhattan plots. The most significant areas are indicated by arrows. (B–D): close-up views of the significant LD block regions (100–139 Mbp) on chromosome 9. (E, H): gene structure with DNA polymorphism. Below the gene structure, boxplots of FWd based on allelic differences of significant SNPs are shown; the width of each box is proportional to the square root of the number of accessions. (F, G): close-up views of the significant LD block regions (142.5–170 Mbp) on chromosome 9. (J, K): close-up views of the significant LD block regions (14–18 Mbp) on chromosome 12.
Of the 101 significant SNPs detected for FWg, 58.4% were located on chromosome 9 (Figures 5F–J). Unlike the other traits examined in the study, significant SNPs for FWg were identified on all chromosomes except chromosomes 3 and 5. While chromosomes 1, 4, 7, and 8 contained one SNP each, chromosomes 2, 6, 10, 11, and 12 contained 3, 11, 9, 2, and 13 significant SNPs for this trait, respectively. Twenty-nine SNPs were commonly detected in at least two years on chromosomes 2, 4, 6, 9, 10, and 12. Seven unique SNPs on chromosomes 6, 7, 8, 9 10, and 12 were detected in 2015; 64 unique SNPs were detected on chromosomes 1, 6, 9, 10, 11, and 12 in 2016; and two unique SNPs were detected on chromosomes 10 and 11 in 2017 (Figure S1).
Figure 5 GWAS identifies significant LD block regions and correlated genes controlling FWg in pepper. (A, D, F, K): chromosome-wide Manhattan plots. The most significant areas are indicated by arrows. The significant areas are shown in detail next to the chromosome-wide Manhattan plot, in which all three years of analysis are represented by various shapes (circle: 2015, triangle: 2016, X: 2017). Under the Manhattan plot, LD blocks are represented by gray boxes. Candidate genes are represented by red dotted lines with the names and positions indicated inside the LD blocks indicated by blue bars. (B) Close-up view of the significant LD block regions (227.5–228.5 Mbp) on chromosome 4. (C) Gene structure with DNA polymorphism. Below the gene structure, boxplots of FWg based on allelic differences of significant SNPs are shown; the width of each box is proportional to square root of the number of accessions. (E) close-up view of the significant LD block regions (194.5–195.5 Mbp) on chromosome 6. (G–J): close-up views of the significant LD block regions (100–144.5 Mbp) on chromosome 9. (L, M): close-up views of the significant LD block regions (14–18 Mbp) on chromosome 12.
PT was associated with 9 SNPs, which were located in the 227, 197–199, 174, 12, and 62–243 Mbp regions of chromosomes 4, 6, 7, 11, and 12, respectively (Figure S1). Among these, seven SNPs (S06_197114855, S06_198980398, S06_199214893, S06_199214897, S07_174200667, S11_12782432, and S12_61862450) were detected only in 2015, S12_243181724 was detected only in 2017, and S04_227983120 was detected in both 2015 and 2016 (Figures 3C–E, Figure S1, Table S6).
A genome-wide association scan also revealed 52 significant SNPs associated with the variation in FP (Table S6). Most of the significant SNPs were located on chromosome 12, while three SNPs were detected in the 220 Mbp region on chromosome 3, and two SNPs were located at 161 and 198 Mbp on chromosome 5, respectively. Inside of chromosome 12, except for two SNPs detected in a 143.3 Mbp region, most significant SNPs were detected near the 211 to 219 Mbp region, with the highest association detected in the H12-0566 block area (Figures 3F–H).
QTLs of Major Fruit-Related Domestication Traits in the RILs
To confirm the GWAS results, we examined QTLs for four major fruit-related domestication traits (FL, FWd, FWg, FP) using 120 RILs derived from a cross between C. annuum ‘Perennial’ and C. annuum ‘Dempsey’ (PDRIL); in these lines, 86 QTLs for 17 horticultural traits were previously mapped (Han et al., 2016). The only difference in the technique used in the current compared to the previous study is that here, we used the genetic map developed from the more recent version of the reference genome (CM334 v.1.6, http://peppergenome.snu.ac.kr). Based on the reference genome, we used 444,405 SNPs from 120 RILs and both parental lines to construct a binmap. Using a sliding window approach, all SNPs were grouped into 2,050 bins (Figure S4, Table S8). The average length of the bins was 0.55 Mb, ranging from 100 kb to 83.5 Mbp. The total genetic distance of the bin map was 1,123.6 cM (Table S9).
Using the same phenotypic information, a total of 17 QTLs were identified (Table 2). Each QTL was named based on an abbreviation of the trait name and the chromosome number following ‘PD_’. For each trait, three to five QTLs were detected, which were distributed throughout chromosomes 2, 3, 4, 5, 7, 9, 10, and 12. The phenotypic variation (R2) explained by each QTL ranged from 8.3% (PD_FP9) to 38.5% (PD_FWd4).
Table 2 QTLs controlling FL, FWd, FWg, and FP detected in PDRIL based on the CM334 v.1.6 reference genome.
Based on the mapping results, five minor QTLs for FL were detected on chromosomes 2, 3, 5, and 9 in one environment. For FWd, four minor QTLs were identified on chromosomes 3, 4, and 7. All major and minor QTLs for FWg were detected on chromosome 7 at 29.6 to 30.7 cM, explaining more than 18.4% (LOD >7.5) of the phenotypic variation (R2) among four environments. For FP, three QTLs (PD_FP9, PD_FP10, PD_FP12.1) were commonly identified in two different environments but explained less than 10% of the phenotypic variation. However, the two remaining QTLs detected on chromosome 12 at 41.3 to 44.7 cM explained higher phenotype variation (>18%) in one environment. Except for PD_FWd3, all QTLs for FL, FWd, and FWg had positive additive effects, meaning that RILs with the maternal genotype had higher values, while all five QTLs for FP showed negative additive effects.
Using the same criteria as Han et al. (2016), QTLs detected in more than two environments with threshold R2 values of 10% were considered to be major QTLs. Only one QTL, PD_FWg7.1, was identified as a major QTL for FWg, with R2 (%) values ranging from 19.7 to 30.2. This large variation in R2 values indicates that FWg is highly affected by genotype × environment interactions in this population.
Candidate Genes Influencing Major Fruit-Related Domestication Traits Under Selection
Based on the GWAS results, we selected 64 significant LD blocks and 230 SNPs related to the five major fruit-related domestication traits to predict candidate genes using Blast2GO. Among the 111 genes identified in the significant LD blocks, 1, 70, 39, and 16 genes were correlated with FL, FWd, FWg, PT, and FP, respectively, with some duplication (Table S7). Based on their predicted functions and communality for two or more closely related traits, 16 genes appeared to have close correlations with major fruit-related domestication traits.
First, a gene (CA.PGA v.1.6.scaffold517.20) located in the 211 Mb region of chromosome 4 in H04-0562 was strongly associated with FL. This gene, which is annotated as low-affinity sulfate transporter 3-like, is located approximately 1.7 kb from S04_211848210. A gene in the same mapping region at a 47 kb distance from SNP_211848210, CA.PGA v.1.6.scaffold517.21, is annotated as Agamous-like MADS-box protein AGL104; this gene appears to be an important regulator of FL (Figures 3A, B).
In the 227 Mb region of chromosome 4, a single gene, CA.PGA v.1.6.scaffold1239.15, was detected for both PT and FWg. This gene, encoding growth-regulating factor 1-like, and is be closely linked with SNP S04_227983120, a significant SNP located inside the 2nd exon (Figures 3C–E and 5A–C). The varieties carrying the G allele had heavier fruits with thicker pericarps than varieties carrying the C allele (Figures 3E and 5C).
CA.PGAv.1.6.scaffold1368.1 (associated with LD block H12-0553) and CA.PGAv.1.6.scaffold1387.3 (associated with LD block H12-0570) were predicted to be very important for FP due to their known associations with this trait. These two genes, which are physically positioned between 211 and 215 Mbp on chromosome 12, encode auxin-binding protein ABP19a-like and the protein BIG GRAIN 1-like A, respectively (Figures 3F–H).
Ten genes were closely related to FWd, including eight genes encoding various transcription factors and hormone-regulated genes on chromosome 9 (CA.PGAv.1.6.scaffold3.11, CA.PGAv.1.6.scaffold3.10, CA.PGAv.1.6.scaffold5.32, CA.PGAv.1.6.scaffold5.22, CA.PGAv.1.6.scaffold5.16, CA.PGAv.1.6.scaffold5.14, CA.PGAv.1.6.scaffold283.11, and CA.PGAv.1.6.scaffold133.5) and two genes on chromosome 12 (CA.PGAv.1.6.scaffold730.39 and CA.PGAv.1.6.scaffold534.6) assembled in seven LD blocks (Figure 4). Among these, CA.PGAv.1.6.scaffold3.10 and CA.PGAv.1.6.scaffold5.16, annotated as transcription repressor OFP12-like and leucine-rich repeat and IQ domain-containing protein 1-like isoform X3, respectively, are homologous to gene family members involved in domestication in tomato (OFP, IQ domain family) (Figures 4B, D). Additionally, two genes (CA.PGAv.1.6.scaffold5.32 and CA.PGAv.1.6.scaffold5.22) associated with two stable SNPs in all three years of the experiment (S09_133144036, S09_136634573), which are annotated as elongation factor 1-beta-like and uncharacterized protein LOC107842678 isoform X1, respectively are also predicted to be important for FWd (Figures 4C, D). In addition, two genes (CA.PGAv.1.6.scaffold5.14, and CA.PGAv.1.6.scaffold133.5), which significant SNPs inside their coding regions, are annotated as mRNA cap guanine-N7 methyltransferase 1 and DNA-directed RNA polymerase II subunit 1, respectively (Figures 4D, E, G, H). Two SNPs (S09_138787607, S09_138787665) are located inside the 13th intron of CA.PGAv.1.6.scaffold5.14. Analysis of allelic frequency showed that plants with the T allele had wider fruits than plants with the C and G alleles (Figure 4E). Another SNP, S09_169434758, was located in the 4th exon of CA.PGAv.1.6.scaffold133.5. Among the 230 accessions, 186 accessions carrying the A allele had narrow fruits (average width of 17.95 mm), while 34 accessions carrying the G allele had relatively wide fruits (average width of 50.45 mm; Figure 4H).
Finally, four genes (CA.PGAv.1.6.scaffold3.11, CA.PGAv.1.6.scaffold283.11, CA.PGAv.1.6.scaffold730.39, and CA.PGAv.1.6.scaffold534.6) were closely related to significant SNPs commonly associated with FWg (Figures 4B, F, J, K).
Nine candidate genes are predicted to regulate FWg (Figure 5). Of these, CA.PGAv.1.6.scaffold1239.15, which regulates PT, as described above, is located on chromosome 4 and contains a significant SNP inside its coding region (Figures 5B, C). Moreover, CA.PGAv.1.6.scaffold422.15, which is annotated as peroxidase 41-like, is located 229.8 kb away from the significant SNP S06_194967541, which was consistently identified all three years of the experiment (Figure 5E). The seven remaining genes, which were commonly identified with FWd-associated genes, play roles in plant immunity and defense mechanisms (Figures 5F–M).
Discussion
GWAS is often used to explore the genetic basis of complex traits in field-grown and horticultural crops due to its efficient detection of many natural allelic variations underlying phenotypic diversity (Brachi et al., 2011). Despite its successful use, however, it is still difficult to link the trait of interest to causal genes due to the widespread existence of population structure inside the diversity panels (Pritchard et al., 2000; Zhang et al., 2009). Population stratification and cryptic relationships can generate spurious associations between phenotypes and unlinked SNPs, leading to false positives (Ioannidis, 2005; Moonesinghe et al., 2007). In the current study, we identified 64 significant LD blocks linked to fruit-related traits and uncovered 16 candidate genes as major genes related to pepper domestication.
Pepper germplasm accessions have been divided into sub-clusters based on species, geographical origin, fruit characteristics, or different routes of introduction (Nicolai et al., 2013; Lee et al., 2016). Similar to previous reports, we identified four distinct sub-populations of pepper based on Capsicum species classification. To improve the reliability and credibility of the association results, we performed GWAS using only the C. annuum sub-cluster, which contains a large number of accessions, with great phenotypic variability but without any strong population stratification. We generated 187,966 genome-wide high-quality SNP markers from the C. annuum sub-cluster of the GWAS population using the GBS method. The LD blocks had an average size of 149 kb, indicating that at least 23,490 genome-wide SNPs are required for GWAS in pepper. Based on the estimated LD block size, the number SNP markers generated in this study is sufficient for GWAS in pepper.
We analyzed marker–trait associations for five major fruit-related domestication traits (FL, FWd, FWg, PT, FP) by GWAS. As a result, we identified 111 candidate genes within the 65 LD blocks. Of these, we selected 16 genes as strong candidate causal genes regulating fruit morphology according to the following criteria: 1) developmental genes known to be related to domestication in other plants; 2) genes within LD blocks containing significant SNPs detected in all three years of the study; and 3) SNP-containing genes associated with more than two traits.
Three genes (CA.PGAv.1.6.scaffold517.21, CA.PGAv.1.6.scaffold3.10, CA.PGAv.1.6.scaffold5.16), which are annotated as Agamous-like MADS-box protein AGL104, transcription repressor OFP12-like, and leucine-rich repeat and IQ domain-containing protein 1-like isoform X3, respectively, satisfied the first criterion, as they belong to the MADS domain subfamily, Ovate Family Protein (OFP) family, and IQ domain family, respectively. The OFP and IQ domain gene families include the well-known ovate and sun genes in tomato (Rodriguez et al., 2011). A nonsense mutation in the ovate gene is responsible for the development of pear-shaped fruit instead of oval-shaped fruit in tomato (Wang et al., 2007; Schmitz et al., 2015). In Arabidopsis, this gene regulates the production of a gibberellic acid (GA) biosynthesis enzyme to control cell elongation (Wang et al., 2007). AGAMOUS-like (AGL) transcription factors, which belong to the plant type I MADS domain subfamily, regulate reproductive development. A number of AGL transcription factor genes are specifically expressed in the central cell of the female gametophyte and endosperm in Arabidopsis (Bemer et al., 2010). Two genes associated with FP (CA.PGAv.1.6.scaffold1368.1, CA.PGAv.1.6.scaffold1387.3) are thought to be important candidates due to their regulation by auxin. The gene CA.PGAv.1.6.scaffold1387.3, which is annotated as BIG GRAIN 1-like A, is homologous to an auxin transport protein gene in Arabidopsis. This gene controls the adaxial–abaxial polarity of the pedicel (Yamaguchi et al., 2007), making it a good candidate gene for FP.
Five genes (CA.PGAv.1.6.scaffold517.21, CA.PGAv.1.6.scaffold517.20, CA.PGAv.1.6.scaffold422.15, CA.PGAv.1.6.scaffold5.32, and CA.PGAv.1.6.scaffold5.22) were chosen as candidate causal genes based on the second criterion: these genes are annotated as Agamous-like MADS-box protein AGL104, low-affinity sulfate transporter 3-like, peroxidase 41-like, elongation factor 1-beta-like, and uncharacterized protein LOC107842678, respectively. The first two genes, which are closely related to FL, are located at 211 Mbp on chromosome 4. In detail, CA.PGAv.1.6.scaffold517.20, a member of sulfate transporter family group 2, might be involved in the internal transport of sulfate between cellular or subcellular compartments within the plant (Hawkesford, 2003). Although sulfate is essential nutrient required for the biosynthesis of a wide range of sulfur-containing compounds, the functions of these genes in plants are unclear (Saito, 2000). The homolog of CA.PGAv.1.6.scaffold422.15 (associated with FWg and located at 194 Mbp on chromosome 6) regulates pollen germination and pollen tube growth (Becker et al., 2003; Wang et al., 2008). The FWd-related genes include CA.PGAv.1.6.scaffold5.22 and CA.PGAv.1.6.scaffold5.32. CA.PGAv.1.6.scaffold5.32 is a homolog of Arabidopsis high amplitude circadian-regulating, which plays fundamental roles in nearly all aspects of plant growth and development (Covington et al., 2008). By contrast, the exact nature of the CA.PGAv.1.6.scaffold5.22 gene homolog has yet to be characterized.
Seven genes (CA.PGAv.1.6.scaffold3.11, CA.PGAv.1.6.scaffold3.10, CA.PGAv.1.6.scaffold5.32, CA.PGAv.1.6.scaffold5.22, CA.PGAv.1.6.scaffold283.11, CA.PGAv.1.6.scaffold730.39, and CA.PGAv.1.6.scaffold534.6) were commonly identified as candidates for both FWd and FWg. Three of these genes (CA.PGAv.1.6.scaffold3.11, CA.PGAv.1.6.scaffold730.39, CA.PGAv.1.6.scaffold534.6) are closely related to plant immune responses and are annotated probable serine/threonine-protein kinase Cx32, probable LRR receptor-like serine/threonine-protein kinase At3g47570, and flower-specific defensin-like, respectively. Besides their major roles, a few studies have focused on their roles in plant growth and development (Chevalier and Walker, 2005; Schulz et al., 2013). CA.PGAv.1.6.scaffold283.11, annotated as calcium-dependent protein kinase 13, is homologous to an Arabidopsis gene encoding a transcriptional regulator essential for Nod-factor-induced gene expression in response to elevated calcium levels, which regulate secondary growth and biomass accumulation (Sehr et al., 2010).
CA.PGAv.1.6.scaffold1239.15 encodes a Growth-regulating factor (GRF) that is correlated with both PT and FWg. GRFs are plant-specific transcription factors that were originally identified for their roles in stem and leaf development. Recent studies have highlighted their importance in other central developmental processes including flower and seed formation, root development, and the coordination of growth processes under adverse environmental conditions (Omidbakhshfard et al., 2015). We subjected the results of our phenotypic survey (conducted for three years to examine morphological traits in pepper) to Pearson correlation (r) analysis, which also supported the GWAS results.
A comparison of the QTLs mapped based on the PDRIL and GWAS results from the GWAS population revealed only one common genetic area associated with FP. Region 141.6 to 144.6 Mbp on chromosome 12 contains three QTLs (PD_FP12.1, PD_FP12.2, PD_FP12.3) and two significant SNPs (S12_143380249, S12_143380271). Inside this common area, two genes were identified (CA.PGAv.1.6.scaffold18.1, and CA.PGAv.1.6.scaffold172.10); these genes are annotated as L-ascorbate oxidase and ELKS/Rab6-interacting/CAST family member 1 isoform X1, respectively. Although the functional relevance of these candidate genes requires further validation, based on their putative functions, they represent strong candidate genes involved in pepper domestication. Among the 17 detected QTLs, one major QTL for FWg, PD_FWg7.1, spanning around 68.9 to 73.6 Mbp on chromosome 7 was identified. In this position, we were able to detect a relatively high peak than the surrounding area in GWAS. However, the P-values of those SNPs (–log10 P-value <2.9) did not pass a significant threshold. Some QTL positions for FWg and FWd were corresponding to QTLs reported by Wu et al. (2019). Unexpectedly, however, most QTLs or significant SNPs in QTL analysis and GWAS for fruit traits were not common. This may be due to several reasons including Beavis effect, differences in models or fundamental differences analysis as suggested Hansson et al. (2018).
In summary, we successfully used GWAS to identify genes responsible for major fruit-related traits in pepper. The significant haplotypes identified in this study provide unique molecular footprints for developing markers for pre-breeding or genomic selection. Future functional validation of the candidate genes identified in this study should provide additional targets for the improvement of major horticultural traits in pepper via breeding.
Data Availability Statement
The datasets generated for this study can be found in The National Agricultural Biotechnology Information Center http://nabic.rda.go.kr/, ID NV-0630-000001.
Author Contributions
Conceptualization: H-YL, B-CK. Data curation: H-YL, N-YR. Formal analysis: H-YL. Funding acquisition: B-CK, J-KK. Investigation: H-YL, N-YR. Methodology: H-YL. Project administration: H-YL, B-CK. Resources: B-CK, N-YR. Software: H-YL, J-HL. Validation: B-CK. Visualization: H-YL. Writing—original draft: H-YL. Writing—review and editing: B-CK, AP.
Funding
This work was carried out with the support of “Cooperative Research Program for Agriculture Science & Technology Development (Project No. PJ01322901)” Rural Development Administration, Republic of Korea.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2020.01100/full#supplementary-material
References
Babu, B. S., Pandravada, S. R., Prasada Rao, R. D. V. J., Anitha, K., Chakrabarty, S. K., Varaprasad, K. S. (2011). Global sources of pepper genetic resources against arthropods, nematodes and pathogens. Crop Prot. 30, 389–400. doi: 10.1016/j.cropro.2010.12.011
Bai, Y., Lindhout, P. (2007). Domestication and breeding of tomatoes: What have we gained and what can we gain in the future? Ann. Bot. 100, 1085–1094. doi: 10.1093/aob/mcm150
Barchi, L., Lefebvre, V., Sage-Palloix, A., Lanteri, S., Palloix, A. (2009). QTL analysis of plant development and fruit traits in pepper and performance of selective phenotyping. Theor. Appl. Genet. 118, 1157–1171. doi: 10.1007/s00122-009-0970-0
Becker, D., Boavida, L. C., Carneiro, J., Haury, M., Feijo, A. (2003). Transcriptional Profiling of Arabidopsis Tissues. Society 133, 713–725. doi: 10.1104/pp.103.028241
Bemer, M., Heijmans, K., Airoldi, C., Davies, B., Angenent, G. C. (2010). An atlas of type I MADS box gene expression during female gametophyte and seed development in Arabidopsis. Plant Physiol. 154, 287–300. doi: 10.1104/pp.110.160770
Bonferroni, C. E. (1936). Teoria statistica delle classi e calcolo delle probabilita. Pubblicazioni del R Istituto Superiore di Scienze Economiche e Commericiali di Firenze. 8, 3–62.
Borovsky, Y., Paran, I. (2011). Characterization of fs10.1, a major QTL controlling fruit elongation in Capsicum. Theor. Appl. Genet. 123 (4), 657–65. doi: 10.1007/s00122-011-1615-7
Brachi, B., Morris, G. P., Borevitz, J. O. (2011). Genome-wide association studies in plants: the missing heritability is in the field. Genom. Biol. 12 (10), 232. doi: 10.1186/gb-2011-12-10-232
Bradbury, P. J., Zhang, Z., Kroon, D. E., Casstevens, T. M., Ramdoss, Y., Buckler, E. S. (2007). TASSEL: software for association mapping of complex traits in diverse samples. Bioinformatics 23 (19), 2633–2635. doi: 10.1093/bioinformatics/btm308
Chaim, A., Ben, P., II, Grube, R. C., Jahn, M., Van Wijk, R., Peleman, J. (2001). QTL mapping of fruit-related traits in pepper (Capsicum annuum). Theor. Appl. Genet. 102 (6–7), 1016–1028. doi: 10.1007/s001220000461
Chaim, A. B., Borovsky, Y., De Jong, W., Paran, I. (2003). Linkage of the A locus for the presence of anthocyanin and fs10.1, a major fruit-shape QTL in pepper. Theor. Appl. Genet. 106 (5), 889–894. doi: 10.1007/s00122-002-1132-9
Chakrabarti, M., Zhang, N., Sauvage, C., Munos, S., Blanca, J., Canizares, J., et al. (2013). A cytochrome P450 regulates a domestication trait in cultivated tomato. Proc. Natl. Acad. Sci. 110 (42), 17125–17130. doi: 10.1073/pnas.1307313110
Chang, C. C., Chow, C. C., Tellier, L. C. A. M., Vattikuti, S., Purcell, S. M., Lee, J. J. (2015). Second-generation PLINK: Rising to the challenge of larger and richer datasets. Gigascience 4 (1), 1–16. doi: 10.1186/s13742-015-0047-8
Che, G., Zhang, X. (2019). Molecular basis of cucumber fruit domestication. Curr. Opin. Plant Biol. 47, 38–46. doi: 10.1016/j.pbi.2018.08.006
Cheng, J., Qin, C., Tang, X., Zhou, H., Hu, Y., Zhao, Z., et al. (2016). Development of a SNP array and its application to genetic mapping and diversity assessment in pepper (Capsicum spp.). Sci. Rep. 6, 1–11. doi: 10.1038/srep33293
Chevalier, D., Walker, J. C. (2005). Functional genomics of protein kinases in plants. Briefings Funct. Genomics Proteomics 3, 362–371. doi: 10.1093/bfgp/3.4.362
Chunthawodtiporn, J., Hill, T., Stoffel, K., Van Deynze, A. (2018). Quantitative Trait Loci Controlling Fruit Size and Other Horticultural Traits in Bell Pepper. Plant Genom. 11 (1), 0. doi: 10.3835/plantgenome2016.12.0125
Colonna, V., D’Agostino, N., Garrison, E., Albrechtsen, A., Meisner, J., Facchiano, A., et al. (2019). Genomic diversity and novel genome-wide association with fruit morphology in Capsicum, from 746k polymorphic sites. Sci. Rep. 9 (1), 10067. doi: 10.1038/s41598-019-46136-5
Covington, M. F., Maloof, J. N., Straume, M., Kay, S. A., Harmer, S. L. (2008). Global transcriptome analysis reveals circadian regulation of key pathways in plant growth and development. Genome Biol. 9. doi: 10.1186/gb-2008-9-8-r130
De Givry, S., Bouchez, M., Chabrier, P., Milan, D., Schiex, T. (2005). CARHTA GENE: multipopulation integrated genetic and radiation hybrid mapping. Bioinformatics 21, 1703–1704. doi: 10.1093/bioinformatics/bti222
Depristo, M. A., Banks, E., Poplin, R., Garimella, K. V., Maguire, J. R., Hartl, C., et al. (2011). A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat. Genet. 43 (5), 491–501. doi: 10.1038/ng.806
Doebley, J. F., Gaut, B. S., Smith, B. D. (2006). The molecular genetics of crop domestication. Cell 127, 1309–1321. doi: 10.1016/j.cell.2006.12.006
Doganlar, S., Frary, A., Daunay, M. C., Lester, R. N., Tanksley, S. D. (2002). Conservation of gene function in the Solanaceae as revealed by comparative mapping of domestication traits in eggplant. Genetics 161, 1713–1726.
Giovannoni, J. (2018). Tomato Multiomics Reveals Consequences of Crop Domestication and Improvement. Cell 172, 6–8. doi: 10.1016/j.cell.2017.12.036
Götz, S., García-Gómez, J. M., Terol, J., Williams, T. D., Nagaraj, S. H., Nueda, M. J., et al. (2008). High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res. 36 (10), 3420–3435. doi: 10.1093/nar/gkn176
Guo, M., Simmons, C. R. (2011). Cell number counts - The fw2.2 and CNR genes and implications for controlling plant fruit and organ size. Plant Sci. 181 (1), 1–7. doi: 10.1016/j.plantsci.2011.03.010
Han, K., Jeong, H. J., Yang, H. B., Kang, S. M., Kwon, J. K., Kim, S., et al. (2016). An ultra-high-density bin map facilitates high-throughput QTL mapping of horticultural traits in pepper (Capsicum annuum). DNA Res. 23 (2), 81–91. doi: 10.1093/dnares/dsv038
Han, K., Lee, H. Y., Ro, N. Y., Hur, O. S., Lee, J. H., Kwon, J. K., et al. (2018). QTL mapping and GWAS reveal candidate genes controlling capsaicinoid content in Capsicum. Plant Biotechnol. J. 16 (9), 1546–1558. doi: 10.1111/pbi.12894
Hansson, B., Sigeman, H., Stervander, M., Tarka, M., Ponnikas, S., Strandh, M., et al. (2018). Contrasting results from GWAS and QTL mapping on wing length in great reed warblers. Mol. Ecol. Resour. 18 (4), 867–876. doi: 10.1111/1755-0998.12785
Hawkesford, M. J. (2003). Transporter gene families in plants: the sulphate transporter gene family-redundancy or specialization? Physiol. Plant 117, 155–163. doi: 10.1034/j.1399-3054.2003.00034.x
Hill, T. A., Chunthawodtiporn, J., Ashrafi, H., Stoffel, K., Weir, A., Van Deynze, A. (2017). Regions underlying population structure and the genomics of organ size determination in Capsicum annuum. Plant Genome 10 (3), 1–14. doi: 10.3835/plantgenome2017.03.0026
Ioannidis, J. P. A. (2005). Why most published research findings are false. PloS Med. 2 (8), 696–701. doi: 10.1371/journal.pmed.0020124
Kaiser, S. (1935). The Inheritance of a Geotropic Response in Capsicum Fruits. Bulletin of the Torrey Botanical Club 62, 75–80. Published by: Torrey Botanical Society.
Kim, S., Park, M., Yeom, S. I., Kim, Y. M., Lee, J. M., Lee, H. A., et al. (2014). Genome sequence of the hot pepper provides insights into the evolution of pungency in Capsicum species. Nat. Genet. 46 (3), 270. doi: 10.1038/ng.2877
Kim, S., Park, J., Yeom, S. I., Kim, Y. M., Seo, E., Kim, K. T., et al. (2017). New reference genome sequences of hot pepper reveal the massive evolution of plant disease-resistance genes by retroduplication. Genome Biol. 18 (1), 1–11. doi: 10.1186/s13059-017-1341-9
Lee, H. Y., Ro, N. Y., Jeong, H. J., Kwon, J. K., Jo, J., Ha, Y., et al. (2016). Genetic diversity and population structure analysis to construct a core collection from a large Capsicum germplasm. BMC Genet. 17 (1), 1–13. doi: 10.1186/s12863-016-0452-8
Lee, J. H., An, J. T., Siddique, M. I., Han, K., Choi, S., Kwon, J. K., et al. (2017). Identification and molecular genetic mapping of Chili veinal mottle virus (ChiVMV) resistance genes in pepper (Capsicum annuum). Mol. Breed. 37 (10), 121. . doi: 10.1007/s11032-017-0717-6
Li, H., Durbin, R. (2010). Fast and accurate long-read alignment with Burrows-Wheeler transform. Bioinformatics 26 (5), 589–95. doi: 10.1093/bioinformatics/btp698
Li, Y., Colleoni, C., Zhang, J., Liang, Q., Hu, Y., Ruess, H., et al. (2018). Genomic Analyses Yield Markers for Identifying Agronomically Important Genes in Potato. Mol. Plant 11, 473–484. doi: 10.1016/j.molp.2018.01.009
Lin, T., Zhu, G., Zhang, J., Xu, X., Yu, Q., Zheng, Z., et al. (2014). Genomic analyses provide insights into the history of tomato breeding. Nat. Genet. 46 (11), 1220–1226. doi: 10.1038/ng.3117
Lipka, A. E., Tian, F., Wang, Q., Peiffer, J., Li, M., Bradbury, P. J., et al. (2012). GAPIT: Genome association and prediction integrated tool. Bioinformatics 28 (18), 2397–29. doi: 10.1093/bioinformatics/bts444
Meyer, R. S., Karol, K. G., Little, D. P., Nee, M. H., Litt, A. (2012). Phylogeographic relationships among Asian eggplants and new perspectives on eggplant domestication. Mol. Phylogenet. Evol. 63, 685–701. doi: 10.1016/j.ympev.2012.02.006
Mimura, Y., Inoue, T., Minamiyama, Y., Kubo, N. (2012). An SSR-based genetic map of pepper (Capsicum annuum L.) serves as an anchor for the alignment of major pepper maps. Breed. Sci. 62, 93–98. doi: 10.1270/jsbbs.62.93
Monforte, A. J., Diaz, A., Caño-Delgado, A., Van Der Knaap, E. (2014). The genetic basis of fruit morphology in horticultural crops: Lessons from tomato and melon. J. Exp. Bot. 65 (16), 4625–4637. doi: 10.1093/jxb/eru017
Moonesinghe, R., Khoury, M. J., Janssens, A. C. J. W. (2007). Most Published Research Findings Are False-But a Little Replication Goes a Long Way. PloS Med. 4 (2), e28. doi: 10.1371/journal.pmed.0040028
Moscone, E. A., Scaldaferro, M. A., Grabiele, M., Cecchini, N. M., Sánchez García, Y., Jarret, R., et al. (2007). The evolution of chili peppers (Capsicum – Solanaceae): A cytogenetic perspective. Acta Hortic. 745, 137–170. doi: 10.17660/ActaHortic.2007.745.5
Nee, M., Bohs, L., Knapp, S. (2006). New species of Solanum and Capsicum (Solanaceae) from Bolivia, with clarification of nomenclature in some Bolivian Solanum. Brittonia 58 (4), 322–356. doi: 10.1663/0007-196X(2006)58[322:NSOSAC]2.0.CO;2
Nicolaï, M., Cantet, M., Lefebvre, V., Sage-Palloix, A. M., Palloix, A. (2013). Genotyping a large collection of pepper (Capsicum spp.) with SSR loci brings new evidence for the wild origin of cultivated C. annuum and the structuring of genetic diversity by human selection of cultivar types. Genet. Resour. Crop Evol. 60, 2375–2390. doi: 10.1007/s10722-013-0006-0
Nimmakayala, P., Abburi, V. L., Saminathan, T., Alaparthi, S. B., Almeida, A., Davenport, B., et al. (2016a). Genome-wide Diversity and Association Mapping for Capsaicinoids and Fruit Weight in Capsicum annuum L. Sci. Rep. 6, 1–14. doi: 10.1038/srep38081
Nimmakayala, P., Abburi, V. L., Saminathan, T., Almeida, A., Davenport, B., Davidson, J., et al. (2016b). Genome-Wide Divergence and Linkage Disequilibrium Analyses for Capsicum baccatum Revealed by Genome-Anchored Single Nucleotide Polymorphisms. Front. Plant Sci. 7, 1–12. doi: 10.3389/fpls.2016.01646
Omidbakhshfard, M. A., Proost, S., Fujikura, U., Mueller-Roeber, B. (2015). Growth-Regulating Factors (GRFs): A Small Transcription Factor Family with Important Functions in Plant Biology. Mol. Plant 8, 998–1010. doi: 10.1016/j.molp.2015.01.013
Paran, I., Van Der Knaap, E. (2007). Genetic and molecular regulation of fruit and plant domestication traits in tomato and pepper. J. Exp. Bot. 58 (14), 3841–3852. doi: 10.1093/jxb/erm257
Pickersgill, B. (2007). Domestication of plants in the Americas: Insights from Mendelian and molecular genetics. Ann. Bot. 100, 925–940. doi: 10.1093/aob/mcm193
Pritchard, J. K., Stephens, M., Rosenberg, N. A., Donnelly, P. (2000). Association mapping in structured populations. Am. J. Hum. Genet. 67 (1), 170–181. doi: 10.1086/302959
Ramchiary, N., Kehie, M., Brahma, V., Kumaria, S., Tandon, P. (2014). Application of genetics and genomics towards Capsicum translational research. Plant Biotechnol. Rep. 8 (2), 101–123. doi: 10.1007/s11816-013-0306-z
Rao, G. U., Ben Chaim, A., Borovsky, Y., Paran, I. (2003). Mapping of yield-related QTLs in pepper in an interspecific cross of Capsicum annuum and C.frutescens. Theor. Appl. Genet. 106 (8), 1457–1466. doi: 10.1007/s00122-003-1204-5
Rodriguez, G. R., Munos, S., Anderson, C., Sim, S.-C., Michel, A., Causse, M., et al. (2011). Distribution of SUN, OVATE, LC, and FAS in the Tomato Germplasm and the Relationship to Fruit Shape Diversity. Plant Physiol. 156 (1), 275–285. doi: 10.1104/pp.110.167577
Sacco, A., Ruggieri, V., Parisi, M., Festa, G., Rigano, M. M., Picarella, M. E., et al. (2015). Exploring a tomato landraces collection for fruit-related traits by the aid of a high-throughput genomic platform. PloS One 10 (9), 1–20. doi: 10.1371/journal.pone.0137139
Saito, K. (2000). Regulation of sulfate transport and synthesis of sulfur-containing amino acids. Curr. Opin. Plant Biol. 3, 188–195. doi: 10.1016/S1369-5266(00)00063-7
Schmitz, A. J., Begcy, K., Sarath, G., Walia, H. (2015). Rice Ovate Family Protein 2 (OFP2) alters hormonal homeostasis and vasculature development. Plant Sci. 241, 177–188. doi: 10.1016/j.plantsci.2015.10.011
Schulz, P., Herde, M., Romeis, T. (2013). Calcium-dependent protein kinases: Hubs in plant stress signaling and development. Plant Physiol. 163, 523–530. doi: 10.1104/pp.113.222539
Sehr, E. M., Agusti, J., Lehner, R., Farmer, E. E., Schwarz, M., Greb, T. (2010). Analysis of secondary growth in the Arabidopsis shoot reveals a positive role of jasmonate signalling in cambium formation. Plant J. 63, 811–822. doi: 10.1111/j.1365-313X.2010.04283.x
Siddique, M. I., Lee, H. Y., Ro, N. Y., Han, K., Venkatesh, J., Solomon, A. M., et al. (2019). Identifying candidate genes for Phytophthora capsici resistance in pepper (Capsicum annuum) via genotyping-by-sequencing-based QTL mapping and genome-wide association study. Sci. Rep. 9 (1), 9962. doi: 10.1038/s41598-019-46342-1
Solomon, A. M., Han, K., Lee, J. H., Lee, H. Y., Jang, S., Kang, B. C. (2019). Genetic diversity and population structure of Ethiopian Capsicum germplasms. PloS One 14 (5), e0216886. doi: 10.1371/journal.pone.0216886
Soyk, S., Lemmon, Z. H., Oved, M., Fisher, J., Liberatore, K. L., Park, S. J., et al. (2017). Bypassing Negative Epistasis on Yield in Tomato Imposed by a Domestication Gene. Cell 169 (6), 1142–1155.e12. doi: 10.1016/j.cell.2017.04.032
Stacklies, W., Redestig, H., Scholz, M., Walther, D., Selbig, J. (2007). pcaMethods-a Bioconductor package providing PCA methods for incomplete Data. Bioinformatics 23 (9), 1164–1167. doi: 10.1093/bioinformatics/btm069
Stitzer, M. C., Ross-Ibarra, J. (2018). Maize domestication and gene interaction. New Phytol. 220, 395–408. doi: 10.1111/nph.15350
Wang, D., Bosland, P. W. (2006). The genes of Capsicum. HortScience 41, 1169–1187. doi: 10.21273/hortsci.41.5.1169
Wang, S., Basten, C. J., Zeng, Z. B. (2012). Windows QTL Cartographer 2.5. Department of Statistics (Raleigh, NC: North Carolina State University).
Wang, L., Li, J., Zhao, J., He, C. (2015). Evolutionary developmental genetics of fruit morphological variation within the Solanaceae. Front. Plant Sci. 6, 1–10. doi: 10.3389/fpls.2015.00248
Wang, S., Chang, Y., Guo, J., Chen, J. G. (2007). Arabidopsis Ovate Family Protein 1 is a transcriptional repressor that suppresses cell elongation. Plant J. 50, 858–872. doi: 10.1111/j.1365-313X.2007.03096.x
Wang, Y., Zhang, W. Z., Song, L. F., Zou, J. J., Su, Z., Wu, W. H. (2008). Transcriptome analyses show changes in gene expression to accompany pollen germination and tube growth in arabidopsis. Plant Physiol. 148, 1201–1211. doi: 10.1104/pp.108.126375
Wu, L., Wang, P., Wang, Y., Chen, Q., Lu, Q., Lou, J., et al. (2019). Genome-wide correlation of 36 agronomic traits in the 287 pepper (Capsicum) accessions obtained from the SLAF-se-based GWAS. Internal J. Mol. Sci. 20 (22), 5675. doi: 10.3390/ijms20225675
Xiao, H., Jiang, N., Schaffner, E., Stockinger, E. J., Van Der Knaap, E. (2008). Variation of Tomato Fruit. Science 319, 1527–30. doi: 10.1126/science.1153040
Yamaguchi, N., Suzuki, M., Fukaki, H., Morita-Terao, M., Tasaka, M., Komeda, Y. (2007). CRM1/BIG-Mediated Auxin Action Regulates Arabidopsis Inflorescence Development. Plant Cell Physiol. 48 (9), 1275–1290. doi: 10.1093/pcp/pcm094
Zeder, M. A., Emshwiller, E., Smith, B. D., Bradley, D. G. (2006). Documenting domestication: the intersection of genetics and archaeology. Trends Genet. 22, 139–155. doi: 10.1016/j.tig.2006.01.007
Zhang, Z., Buckler, E. S., Casstevens, T. M., Bradbury, P. J. (2009). Software engineering the mixed model for genome-wide association studies on large samples. Brief Bioinform. 10 (6), 664–675. doi: 10.1093/bib/bbp050
Keywords: Capsicum, domestication, fruit-related traits, genotype-by-sequencing, genome-wide association study, quantitative trait locus, linkage disequilibrium
Citation: Lee H-Y, Ro N-Y, Patil A, Lee J-H, Kwon J-K and Kang B-C (2020) Uncovering Candidate Genes Controlling Major Fruit-Related Traits in Pepper via Genotype-by-Sequencing Based QTL Mapping and Genome-Wide Association Study. Front. Plant Sci. 11:1100. doi: 10.3389/fpls.2020.01100
Received: 10 April 2020; Accepted: 03 July 2020;
Published: 23 July 2020.
Edited by:
José Antonio Fernández, University of Castilla La Mancha, SpainReviewed by:
Francesca Taranto, Italian National Research Council, ItalySanghyeob Lee, Sejong University, South Korea
Copyright © 2020 Lee, Ro, Patil, Lee, Kwon and Kang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Byoung-Cheorl Kang, bk54@snu.ac.kr