Genome-Wide Correlation of 36 Agronomic Traits in the 287 Pepper (Capsicum) Accessions Obtained from the SLAF-seq-Based GWAS

Genome-wide association study and candidate gene identification for agronomic traits in 182 upward-growing fruits of C. frutescens and C. annuum

Pepper agronomic traits serve as pivotal indicators for characterizing germplasm attributes and correlations. It is important to study differential genotypic variation through phenotypic differences of target traits. Whole genome resequencing was used to sequence the whole genome among different individuals of species with known reference genomes and annotations, and based on this, differential analyses of individuals or populations were carried out to identify SNPs for agronomic traits related to pepper. This study conducted a genome-wide association study encompassing 26 key agronomic traits in 182 upward-growing fruits of C. frutescens and C. annuum. The population structure (phylogenetics, population structure, population principal component analysis, genetic relationship) and linkage disequilibrium analysis were realized to ensure the accuracy and reliability of GWAS results, and the optimal statistical model was determined. A total of 929 SNPs significantly associated with 26 agronomic traits, were identified, alongside the detection of 519 candidate genes within 100 kb region adjacent to these SNPs. Additionally, through gene annotation and expression pattern scrutiny, genes such as GAUT1, COP10, and DDB1 correlated with fruit traits in Capsicum frutescens and Capsicum annuum were validated via qRT-PCR. In the CH20 (Capsicum annuum) and YB-4 (Capsicum frutescens) cultivars, GAUT1 and COP10 were cloned with cDNA lengths of 1065 bp and 561 bp, respectively, exhibiting only a small number of single nucleotide variations and nucleotide deletions. This validation provides a robust reference for molecular marker-assisted breeding of pepper agronomic traits, offering both genetic resources and theoretical foundations for future endeavors in molecular marker-assisted breeding for pepper.

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.

A Genome-Wide Association Study Reveals the Genetic Mechanisms of Nutrient Accumulation in Spinach

Spinach is a significant source of vitamins, minerals, and antioxidants. These nutrients make it delicious and beneficial for human health. However, the genetic mechanism underlying the accumulation of nutrients in spinach remains unclear. In this study, we analyzed the content of chlorophyll a, chlorophyll b, oxalate, nitrate, crude fiber, soluble sugars, manganese, copper, and iron in 62 different spinach accessions. Additionally, 3,356,182 high-quality, single-nucleotide polymorphisms were found using resequencing and used in a genome-wide association study. A total of 2077 loci were discovered that significantly correlated with the concentrations of the nutritional elements. Data mining identified key genes in these intervals for four traits: chlorophyll, oxalate, soluble sugar, and Fe. Our study provides insights into the genetic architecture of nutrient variation and facilitates spinach breeding for good nutrition.

Multi-environment association study highlights candidate genes for robust agronomic quantitative trait loci in a novel worldwide Capsicum core collection

Investigating crop diversity through genome-wide association studies (GWAS) on core collections helps in deciphering the genetic determinants of complex quantitative traits. Using the G2P-SOL project world collection of 10 038 wild and cultivated Capsicum accessions from 10 major genebanks, we assembled a core collection of 423 accessions representing the known genetic diversity. Since complex traits are often highly dependent upon environmental variables and genotype-by-environment (G × E) interactions, multi-environment GWAS with a 10 195-marker genotypic matrix were conducted on a highly diverse subset of 350 Capsicum annuum accessions, extensively phenotyped in up to six independent trials from five climatically differing countries. Environment-specific and multi-environment quantitative trait loci (QTLs) were detected for 23 diverse agronomic traits. We identified 97 candidate genes potentially implicated in 53 of the most robust and high-confidence QTLs for fruit flavor, color, size, and shape traits, and for plant productivity, vigor, and earliness traits. Investigating the genetic architecture of agronomic traits in this way will assist the development of genetic markers and pave the way for marker-assisted selection. The G2P-SOL pepper core collection will be available upon request as a unique and universal resource for further exploitation in future gene discovery and marker-assisted breeding efforts by the pepper community.

The genetic architecture of prolificacy in maize revealed by association mapping and bulk segregant analysis

KEY MESSAGE

Here, we revealed maize prolificacy highly correlated with domestication and identified a causal gene ZmEN1 located in one novel QTL qGEN261 that regulating maize prolificacy by using multiple-mapping methods. The development of maize prolificacy (EN) is crucial for enhancing yield and breeding specialty varieties. To achieve this goal, we employed a genome-wide association study (GWAS) to analyze the genetic architecture of EN in maize. Using 492 inbred lines with a wide range of EN variability, our results demonstrated significant differences in genetic, environmental, and interaction effects. The broad-sense heritability (H2) of EN was 0.60. Through GWAS, we identified 527 significant single nucleotide polymorphisms (SNPs), involved 290 quantitative trait loci (QTL) and 806 genes. Of these SNPs, 18 and 509 were classified as major effect loci and minor loci, respectively. In addition, we performed a bulk segregant analysis (BSA) in an F2 population constructed by a few-ears line Zheng58 and a multi-ears line 647. Our BSA results identified one significant QTL, qBEN1. Importantly, combining the GWAS and BSA, four co-located QTL, involving six genes, were identified. Three of them were expressed in vegetative meristem, shoot tip, internode and tip of ear primordium, with ZmEN1, encodes an unknown auxin-like protein, having the highest expression level in these tissues. It suggested that ZmEN1 plays a crucial role in promoting axillary bud and tillering to encourage the formation of prolificacy. Haplotype analysis of ZmEN1 revealed significant differences between different haplotypes, with inbred lines carrying hap6 having more EN. Overall, this is the first report about using GWAS and BSA to dissect the genetic architecture of EN in maize, which can be valuable for breeding specialty maize varieties and improving maize yield.

Fine mapping of Rf2, a minor Restorer-of-fertility (Rf) gene for cytoplasmic male sterility in chili pepper G164 (Capsicum annuum L.)

Genome re-sequencing and recombination analyses identified Capana06g000193 as a strong candidate for the minor male fertility restoration locus Rf2 in chili pepper G164 harboring two dominant male fertility restoration genes. Male fertility restoration genes of chili pepper restorer line G164 (Capsicum annuum L.) were studied using molecular marker genotypes of an F₂ population (7G) of G164 crossed with the cytoplasmic male sterility line 77013A. The ratio of sterile to fertile single plants in the F₂ population was 1:15. This result indicates that chili pepper G164 has two dominant restoration genes, which we designated as Rf1 and Rf2. An individual plant recessive for Rf1 and heterozygous for Rf2, 7G-112 (rf1rf1Rf2rf2), was identified by molecular marker selection and genetic analysis, and a single Rf2 gene-segregating population with a 3:1 ratio of fertile to sterile plants was developed from the self-pollination of male fertile individuals of 77013A and 7G-112 hybrid progeny. Bulk segregant analysis of fertile and sterile pools from the segregating populations was used to genetically map Rf2 to a 3.1-Mb region on chromosome 6. Rf2 was further narrowed to a 179.3-kb interval through recombination analysis of molecular markers and obtained the most likely candidate gene, Capana06g000193.

Fine Mapping and Gene Analysis of restorer-of-fertility Gene CaRfHZ in Pepper (Capsicum annuum L.)

Cytoplasmic male sterility (CMS) is a common biological phenomenon used in hybrid production of peppers (Capsicum annuum L.). Although several restorer-of-fertility (Rf) genes of pepper CMS lines have been mapped, there is no report that the Rf gene with clear gene function has been isolated. Here, pepper CMS line HZ1A and its restorer line HZ1C were used to construct (HZ1A × HZ1C) F₂ populations and map the Rf gene. A single dominant gene CaRfHZ conferred male fertility according to inheritance analysis. Using sterile plants from (HZ1A × HZ1C) F₂ populations and bulked segregant analysis (BSA), the CaRf_HZ gene was mapped between P06gInDel-66 and P06gInDel-89 on chromosome 6. This region spans 533.81 kb, where four genes are annotated according to Zunla-1 V2.0 gene models. Based on the analysis of genomic DNA sequences, gene expressions, and protein structures, Capana06g002968 was proposed as the strongest candidate for the CaRf_HZ gene. Our results may help with hybrid pepper breeding and to elucidate the mechanism of male fertility restoration in peppers.

Chile Pepper (Capsicum) Breeding and Improvement in the "Multi-Omics" Era

Chile pepper (Capsicum spp.) is a major culinary, medicinal, and economic crop in most areas of the world. For more than hundreds of years, chile peppers have "defined" the state of New Mexico, USA. The official state question, "Red or Green?" refers to the preference for either red or the green stage of chile pepper, respectively, reflects the value of these important commodities. The presence of major diseases, low yields, decreased acreages, and costs associated with manual labor limit production in all growing regions of the world. The New Mexico State University (NMSU) Chile Pepper Breeding Program continues to serve as a key player in the development of improved chile pepper varieties for growers and in discoveries that assist plant breeders worldwide. Among the traits of interest for genetic improvement include yield, disease resistance, flavor, and mechanical harvestability. While progress has been made, the use of conventional breeding approaches has yet to fully address producer and consumer demand for these traits in available cultivars. Recent developments in "multi-omics," that is, the simultaneous application of multiple omics approaches to study biological systems, have allowed the genetic dissection of important phenotypes. Given the current needs and production constraints, and the availability of multi-omics tools, it would be relevant to examine the application of these approaches in chile pepper breeding and improvement. In this review, we summarize the major developments in chile pepper breeding and present novel tools that can be implemented to facilitate genetic improvement. In the future, chile pepper improvement is anticipated to be more data and multi-omics driven as more advanced genetics, breeding, and phenotyping tools are developed.

Mapping of CaPP2C35 involved in the formation of light-green immature pepper (Capsicum annuum L.) fruits via GWAS and BSA

Genome-wide association study, bulked segregant analysis, and genetic analysis delimited the LG locus controlling light-green immature pepper fruits into a 35.07 kbp region on chromosome 10. A strong candidate gene, CaPP2C35, was identified in this region. In pepper (Capsicum annuum L.), the common colors of immature fruits are yellowish white, milky yellow, green, purple, and purplish black. Genes related to dark green, white, and purple immature fruits have been cloned; however, only a few studies have investigated light-green immature fruits. Here, we performed a genetic study using light-green (17C827) and green (17C658) immature fruits. The light-green color of immature fruits was controlled by a single locus-dominant genetic trait compared with the green color of immature fruits. We also performed a genome-wide association study and bulked segregant analysis of immature-fruit color and mapped the LG locus to a 35.07 kbp region on chromosome 10. Only one gene, Capana10g001710, was found in this region. A G-A substitution occurred at the 313th base of the Capana10g001710 coding sequence in 17C827, resulting in the conversion of the α-helix of its encoded PP2C35 protein into a β-fold. The expression of Capana10g001710 (termed CaPP2C35) in 17C827 was significantly higher than in 17C658. Silencing CaPP2C35 in 17C827 resulted in an increase in chlorophyll content in the exocarp and the appearance of green stripes on the surface of the fruit. These results indicate that CaPP2C35 may be involved in the formation of light-green immature fruits by regulating the accumulation of chlorophyll content in the exocarp. Thus, these findings lay the foundation for further studies and genetic improvement of immature-fruit color in pepper.

Hairiness Gene Regulated Multicellular, Non-Glandular Trichome Formation in Pepper Species

Trichomes are unicellular or multicellular epidermal structures that play a defensive role against environmental stresses. Although unicellular trichomes have been extensively studied as a mechanistic model, the genes involved in multicellular trichome formation are not well understood. In this study, we first classified the trichome morphology structures in Capsicum species using 280 diverse peppers. We cloned a key gene (Hairiness) on chromosome 10, which mainly controlled the formation of multicellular non-glandular trichomes (types II, III, and V). Hairiness encodes a Cys2-His2 zinc-finger protein, and virus-induced gene silencing of the gene resulted in a hairless phenotype. Differential expression of Hairiness between the hairiness and hairless lines was due to variations in promoter sequences. Transgenic experiments verified the hypothesis that the promoter of Hairiness in the hairless line had extremely low activity causing a hairless phenotype. Hair controlled the formation of type I glandular trichomes in tomatoes, which was due to nucleotide differences. Taken together, our findings suggest that the regulation of multicellular trichome formation might have similar pathways, but the gene could perform slightly different functions in crops.

Comparative Transcriptome Analysis of the Anthers from the Cytoplasmic Male-Sterile Pepper Line HZ1A and Its Maintainer Line HZ1B

Cytoplasmic male-sterility (CMS) is important for the utilization of crop heterosis and study of the molecular mechanisms involved in CMS could improve breeding programs. In the present study, anthers of the pepper CMS line HZ1A and its maintainer line HZ1B were collected from stages S1, S2, and S3 for transcriptome sequencing. A total of 47.95 million clean reads were obtained, and the reads were assembled into 31,603 unigenes. We obtained 42 (27 up-regulated and 15 down-regulated), 691 (346 up-regulated and 345 down-regulated), and 709 (281 up-regulated and 428 down-regulated) differentially expressed genes (DEGs) in stages S1, S2, and S3, respectively. Through Gene Ontology (GO) analysis, the DEGs were found to be composed of 46 functional groups. Two GO terms involved in photosynthesis, photosynthesis (GO:0015986) and photosystem I (GO:0009522), may be related to CMS. Through Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, oxidative phosphorylation (ko00190) and phenylpropanoid biosynthesis (ko00940) were significantly enriched in the S1 and S2 stages, respectively. Through the analysis of 104 lipid metabolism-related DEGs, four significantly enriched KEGG pathways may help to regulate male sterility during anther development. The mitochondrial genes orf470 and atp6 were identified as candidate genes of male sterility for the CMS line HZ1A. Overall, the results will provide insights into the molecular mechanisms of pepper CMS.

Development and application of KASP markers associated with Restorer-of-fertility gene in Capsicum annuum L

Fertility restoration of cytoplasmic male sterility (CMS) in Capsicum annuum is controlled by multiple alleles of Restorer-of-fertility (Rf) genes. The isolation of additional Rf genes should therefore enrich the knowledge of CMS/Rf systems and accelerate their exploitation in hybrid seed production. In this study, the fertility restorer gene CaRfm of '0601 M', a non-pungent bell pepper, was genetically mapped to a 1.2-cM region flanked by KASP markers S761 and S183. CaRfm was then physically mapped to a 128.96-Kb interval predicted from 24 recombinants with two co-segregated markers, S423 and S424. CaPPR6 encoding a pentatricopeptide repeat (PPR) protein was suggested as the most likely candidate gene for the CaRfm locus on the basis of sequence alignment as well as genotyping of tightly linked markers. In addition, molecular markers S1597 and S1609, which are immediately adjacent to CaRfm at 15.7 and 57.8-Kb respectively, were developed and applied to marker-assisted selection. The results provided friendly markers for breeding pepper restorer lines and laid the foundation for elucidating the male fertility restoration mechanism.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01109-9.

The APETALA2 homolog CaFFN regulates flowering time in pepper

Flowering time is an important agronomic trait that contributes to fitness in plants. However, the genetic basis of flowering time has not been extensively studied in pepper. To understand the genetics underlying flowering time, we constructed an F₂ population by crossing a spontaneous early flowering mutant and a late-flowering pepper line. Using bulked segregant RNA-seq, a major locus controlling flowering time in this population was mapped to the end of chromosome 2. An APETALA2 (AP2) homolog (CaFFN) cosegregated with flowering time in 297 individuals of the F₂ population. A comparison between the parents revealed a naturally occurring rare SNP (SNP2T > C) that resulted in the loss of a start codon in CaFFN in the early flowering mutant. Transgenic Nicotiana benthamiana plants with high CaFFN expression exhibited a delay in flowering time and floral patterning defects. On the other hand, pepper plants with CaFFN silencing flowered early. Therefore, the CaFFN gene acts as a flowering repressor in pepper. CaFFN may function as a transcriptional activator to activate the expression of CaAGL15 and miR156e and as a transcriptional repressor to repress the expression of CaAG, CaAP1, CaSEP3, CaSOC1, and miR172b based on a qRT-PCR assay. Direct activation of CaAGL15 by CaFFN was detected using yeast one-hybrid and dual-luciferase reporter assays, consistent with the hypothesis that CaFFN regulates flowering time. Moreover, the CaFFN gene association analysis revealed a significant association with flowering time in a natural pepper population, indicating that the CaFFN gene has a broad effect on flowering time in pepper. Finally, the phylogeny, evolutionary expansion and expression patterns of CaFFN/AP2 homologs were analyzed to provide valuable insight into CaFFN. This study increases our understanding of the involvement of CaFFN in controlling flowering time in pepper, thus making CaFFN a target gene for breeding early maturing pepper.

Molecular Breeding Strategy and Challenges Towards Improvement of Downy Mildew Resistance in Cauliflower (Brassica oleracea var. botrytis L.)

Cauliflower (Brassica oleracea var. botrytis L.) is one of the important, nutritious and healthy vegetable crops grown and consumed worldwide. But its production is constrained by several destructive fungal diseases and most importantly, downy mildew leading to severe yield and quality losses. For sustainable cauliflower production, developing resistant varieties/hybrids with durable resistance against broad-spectrum of pathogens is the best strategy for a long term and reliable solution. Identification of novel resistant resources, knowledge of the genetics of resistance, mapping and cloning of resistance QTLs and identification of candidate genes would facilitate molecular breeding for disease resistance in cauliflower. Advent of next-generation sequencing technologies (NGS) and publishing of draft genome sequence of cauliflower has opened the flood gate for new possibilities to develop enormous amount of genomic resources leading to mapping and cloning of resistance QTLs. In cauliflower, several molecular breeding approaches such as QTL mapping, marker-assisted backcrossing, gene pyramiding have been carried out to develop new resistant cultivars. Marker-assisted selection (MAS) would be beneficial in improving the precision in the selection of improved cultivars against multiple pathogens. This comprehensive review emphasizes the fascinating recent advances made in the application of molecular breeding approach for resistance against an important pathogen; Downy Mildew (Hyaloperonospora parasitica) affecting cauliflower and Brassica oleracea crops and highlights the QTLs identified imparting resistance against this pathogen. We have also emphasized the critical research areas as future perspectives to bridge the gap between availability of genomic resources and its utility in identifying resistance genes/QTLs to breed downy mildew resistant cultivars. Additionally, we have also discussed the challenges and the way forward to realize the full potential of molecular breeding for downy mildew resistance by integrating marker technology with conventional breeding in the post-genomics era. All this information will undoubtedly provide new insights to the researchers in formulating future breeding strategies in cauliflower to develop durable resistant cultivars against the major pathogens in general and downy mildew in particular.

A Comprehensive Transcriptional Profiling of Pepper Responses to Root-Knot Nematode

Genetic resistance remains a key component in integrated pest management systems. The cosmopolitan root-knot nematode (RKN; Meloidogyne spp.) proves a significant management challenge as virulence and pathogenicity vary among and within species. RKN greatly reduces commercial bell pepper yield, and breeding programs continuously develop cultivars to emerging nematode threats. However, there is a lack of knowledge concerning the nature and forms of nematode resistance. Defining how resistant and susceptible pepper cultivars mount defenses against RKN attacks can help inform breeding programs. Here, we characterized the transcriptional responses of the highly related resistant (Charleston Belle) and susceptible (Keystone Resistance Giant) pepper cultivars throughout early nematode infection stages. Comprehensive transcriptomic sequencing of resistant and susceptible cultivar roots with or without Meloidogyneincognita infection over three-time points; covering early penetration (1-day), through feeding site maintenance (7-days post-inoculation), produced > 300 million high quality reads. Close examination of chromosome P9, on which nematode resistance hotspots are located, showed more differentially expressed genes were upregulated in resistant cultivar at day 1 when compared to the susceptible cultivar. Our comprehensive approach to transcriptomic profiling of pepper resistance revealed novel insights into how RKN causes disease and the plant responses mounted to counter nematode attack. This work broadens the definition of resistance from a single loci concept to a more complex array of interrelated pathways. Focus on these pathways in breeding programs may provide more sustainable and enduring forms of resistance.

Genetic diversity and population structure of eddoe taro in China using genome-wide SNP markers

Taro (Colocasia esculenta) is an important root and tuber crop cultivated worldwide. There are two main types of taro that vary in morphology of corm and cormel, ‘dasheen’ and ‘eddoe’. The eddoe type (Colocasia esculenta var. antiquorium) is predominantly distributed throughout China. Characterizing the genetic diversity present in the germplasm bank of taro is fundamental to better manage, conserve and utilize the genetic resources of this species. In this study, the genetic diversity of 234 taro accessions from 16 provinces of China was assessed using 132,869 single nucleotide polymorphism (SNP) markers identified by specific length amplified fragment-sequencing (SLAF-seq). Population structure and principal component analysis permitted the accessions to be categorized into eight groups. The genetic diversity and population differentiation of the eight groups were evaluated using the characterized SNPs. Analysis of molecular variance showed that the variation among eight inferred groups was higher than that within groups, while a relatively small variance was found among the two morphological types and 16 collection regions. Further, a core germplasm set comprising 41 taro accessions that maintained the genetic diversity of the entire collection was developed based on the genotype. This research is expected to be valuable for genetic characterization, germplasm conservation, and breeding of taro.

Construction of a High-Density Genetic Map Based on SLAF Markers and QTL Analysis of Leaf Size in Rice

Leaf shape is an important agronomic trait for constructing an ideal plant type in rice, and high-density genetic map is facilitative in improving accuracy and efficiency for quantitative trait loci (QTL) analysis of leaf trait. In this study, a high-density genetic map contained 10,760 specific length amplified fragment sequencing (SLAF) markers was established based on 149 recombinant inbred lines (RILs) derived from the cross between Rekuangeng (RKG) and Taizhong1 (TN1), which exhibited 1,613.59 cM map distance with an average interval of 0.17 cM. A total of 24 QTLs were detected and explained the phenotypic variance ranged from 9% to 33.8% related to the leaf morphology across two areas. Among them, one uncloned major QTL qTLLW1 (qTLL1 and qTLLW1) involved in regulating leaf length and leaf width with max 33.8% and 22.5% phenotypic variance respectively was located on chromosome 1, and another major locus qTLW4 affecting leaf width accounted for max 25.3% phenotypic variance was mapped on chromosome 4. Fine mapping and qRT-PCR expression analysis indicated that qTLW4 may be allelic to NAL1 (Narrow leaf 1) gene.

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.