An ultra-high-density bin map facilitates high-throughput QTL mapping of horticultural traits in pepper (Capsicum annuum)

Identification of QTLs involved in destemming and fruit quality for mechanical harvesting of New Mexico pod-type green chile

Introduction

Domestic production of pepper (Capsicum spp.) is shrinking while demand within the US is growing. Lack of availability and cost of labor often present an obstacle for domestic producers both practically and economically. As a result, switching to harvesting peppers mechanically is anticipated as a key strategy to help domestic producers compete in the international market. Mechanical harvest efficiency can be improved through breeding. One important trait that mechanical harvest compatible material should have is an easy destemming trait: low force separation of the pedicel and calyx from the fruit.

Methods

To detect the genetic sources underlying a novel easy destemming trait for the purpose of future breeding efforts in New Mexico pod-type green chile, we performed QTL analysis on three F2:F3 populations, coming from three New Mexico pod-type varieties: 'NuMex Odyssey,' 'NuMex Iliad,' and 'NuMex Joe E. Parker,' each crossed with a parent with an easy destemming trait: MUC14. Genotyping was done through genotyping by sequencing (GBS) and phenotyping was done for destemming and fruit trait measurements. Correlations between measurements were found through the R package hmisc and QTL analysis was done through R/qtl.

Results

A strong relationship was seen between destemming and aspects of fruit morphology, particularly, destemming force and fruit width (Pearson's correlation coefficient r=0.75). Major QTLs for destemming and fruit size were discovered. Of these, the largest destemming force QTLs for all populations (PVE=34.5-69.9%) were on chromosome 10, and in two populations QTLs for destemming force were found on chromosome 3 (Percent Variance Explained (PVE)=10.7-18.8%). Fruit size-related QTLs in all populations colocalized in these same areas on chromosomes 3 and 10.

Discussion

This suggests that fruit shape may be genetically linked to destemming, and breeders interested in selecting for easy destemming pepper will also have to pay attention to fruit size and shape.

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.

Insights into the genetic architecture of Phytophthora capsici root rot resistance in chile pepper (Capsicum spp.) from multi-locus genome-wide association study

BACKGROUND

Phytophthora root rot, a major constraint in chile pepper production worldwide, is caused by the soil-borne oomycete, Phytophthora capsici. This study aimed to detect significant regions in the Capsicum genome linked to Phytophthora root rot resistance using a panel consisting of 157 Capsicum spp. genotypes. Multi-locus genome wide association study (GWAS) was conducted using single nucleotide polymorphism (SNP) markers derived from genotyping-by-sequencing (GBS). Individual plants were separately inoculated with P. capsici isolates, 'PWB-185', 'PWB-186', and '6347', at the 4-8 leaf stage and were scored for disease symptoms up to 14-days post-inoculation. Disease scores were used to calculate disease parameters including disease severity index percentage, percent of resistant plants, area under disease progress curve, and estimated marginal means for each genotype.

RESULTS

Most of the genotypes displayed root rot symptoms, whereas five accessions were completely resistant to all the isolates and displayed no symptoms of infection. A total of 55,117 SNP markers derived from GBS were used to perform multi-locus GWAS which identified 330 significant SNP markers associated with disease resistance. Of these, 56 SNP markers distributed across all the 12 chromosomes were common across the isolates, indicating association with more durable resistance. Candidate genes including nucleotide-binding site leucine-rich repeat (NBS-LRR), systemic acquired resistance (SAR8.2), and receptor-like kinase (RLKs), were identified within 0.5 Mb of the associated markers.

CONCLUSIONS

Results will be used to improve resistance to Phytophthora root rot in chile pepper by the development of Kompetitive allele-specific markers (KASP®) for marker validation, genomewide selection, and marker-assisted breeding.

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.

Identification of quantitative trait loci associated with bacterial spot race T4 resistance in intra-specific populations of tomato (Solanum lycopersicum L.)

Bacterial spot of tomato is a serious disease caused by at least four species and four races of Xanthomonas- X. euvesicatoria (race T1), X. vesicatoria (race T2), X. perforans (race T3 and T4), and X. gardneri, with X. perforans race T4 being predominant in the southeast USA. Practical management of this disease is challenging because of the need for more effective chemicals and commercially resistant cultivars. Identification of genetic resistance is the first step to developing a disease-resistant variety. The objective of this study was to identify quantitative trait loci (QTL) conferring resistance to race T4 in two independent recombinant inbred lines (RILs) populations NC 10204 (intra-specific) and NC 13666 (interspecific) developed by crossing NC 30P x NC22L-1(2008) and NC 1CELBR x PI 270443, respectively. Seven QTLs on chromosomes 2, 6, 7, 11, and 12 were identified in NC 10204. The QTL on chromosome 6 explained the highest percentage of phenotypic variance (up to 21.3%), followed by the QTL on chromosome 12 (up to 8.2%). On the other hand, the QTLs on chromosomes 1, 3, 4, 6, 7, 8, 9, and 11 were detected in NC 13666. The QTLs on chromosomes 6, 7, and 11 were co-located in NC 10204 and NC 13666 populations. The donor of the resistance associated with these QTL in NC 10204 is a released breeding line with superior horticultural traits. Therefore, both the donor parent and the QTL information will be useful in tomato breeding programs as there will be minimal linkage drag associated with the bacterial spot resistance.

Mining the candidate genes of rice panicle traits via a genome-wide association study

Panicle traits are important for improving the panicle architecture and grain yield of rice. Therefore, we performed a genome-wide association study (GWAS) to analyze and determine the genetic determinants of five panicle traits. A total of 1.29 million single nucleotide polymorphism (SNP) loci were detected in 162 rice materials. We carried out a GWAS of panicle length (PL), total grain number per panicle (TGP), filled grain number per panicle (FGP), seed setting rate (SSR) and grain weight per panicle (GWP) in 2019, 2020 and 2021. Four quantitative trait loci (QTLs) for PL were detected on chromosomes 1, 6, and 9; one QTL for TGP, FGP, and GWP was detected on chromosome 4; two QTLs for FGP were detected on chromosomes 4 and 7; and one QTL for SSR was detected on chromosome 1. These QTLs were detected via a general linear model (GLM) and mixed linear model (MLM) in both years of the study period. In this study, the genomic best linear unbiased prediction (BLUP) method was used to verify the accuracy of the GWAS results. There are nine QTLs were both detected by the multi-environment GWAS method and the BLUP method. Moreover, further analysis revealed that three candidate genes, LOC_Os01g43700, LOC_Os09g25784, and LOC_Os04g47890, may be significantly related to panicle traits of rice. Haplotype analysis indicated that LOC_Os01g43700 and LOC_Os09g25784 are highly associated with PL and that LOC_Os04g47890 is highly associated with TGP, FGP, and GWP. Our results offer essential genetic information for the molecular improvement of panicle traits. The identified candidate genes and elite haplotypes could be used in marker-assisted selection to improve rice yield through pyramid breeding.

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

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

Pepper variome reveals the history and key loci associated with fruit domestication and diversification

Pepper (Capsicum spp.) is an important vegetable crop that provides a unique pungent sensation when eaten. Through construction of a pepper variome map, we examined the main groups that emerged during domestication and breeding of C. annuum, their relationships and temporal succession, and the molecular events underlying the main transitions. The results showed that the initial differentiation in fruit shape and pungency, increase in fruit weight, and transition from erect to pendent fruits, as well as the recent appearance of large, blocky, sweet fruits (bell peppers), were accompanied by strong selection/fixation of key alleles and introgressions in two large genomic regions. Furthermore, we identified Up, which encodes a BIG GRAIN protein involved in auxin transport, as a key domestication gene that controls erect vs pendent fruit orientation. The up mutation gained increased expression especially in the fruit pedicel through a 579-bp sequence deletion in its 5' upstream region, resulting in the phenotype of pendent fruit. The function of Up was confirmed by virus-induced gene silencing. Taken together, these findings constitute a cornerstone for understanding the domestication and differentiation of a key horticultural crop.

Genetic mapping revealed that the Pun2 gene in Capsicum chacoense encodes a putative aminotransferase

Several genes regulating capsaicinoid biosynthesis including Pun1 (also known as CS), Pun3, pAMT, and CaKR1 have been studied. However, the gene encoded by Pun2 in the non-pungent Capsicum chacoense is unknown. This study aimed to identify the Pun2 gene by genetic mapping using interspecific (C. chacoense × Capsicum annuum) and intraspecific (C. chacoense × C. chacoense) populations. QTL mapping using the interspecific F₂ population revealed two major QTLs on chromosomes 3 and 9. Two bin markers within the QTL regions on two chromosomes were highly correlated with the capsaicinoid content in the interspecific population. The major QTL, Pun2_PJ_Gibbs_3.11 on chromosome 3, contained the pAMT gene, indicating that the non-pungency of C. chacoense may be attributed to a mutation in the pAMT gene. Sequence analysis revealed a 7 bp nucleotide insertion in the 8^th exon of pAMT of the non-pungent C. chacoense. This mutation resulted in the generation of an early stop codon, resulting in a truncated mutant lacking the PLP binding site, which is critical for pAMT enzymatic activity. This insertion co-segregated with the pungency phenotype in the intraspecific F₂ population. We named this novel pAMT allele pamt¹¹ . Taken together, these data indicate that the non-pungency of C. chacoense is due to the non-functional pAMT allele, and Pun2 encodes the pAMT gene.

Teosinte confers specific alleles and yield potential to maize improvement

Teosinte improves maize grain yield and broadens the maize germplasm. Seventy-one quantitative trait loci associated with 24 differential traits between maize and teosinte were identified. Maize is a major cereal crop with a narrow germplasm that has limited its production and breeding progress. Teosinte, an ancestor of maize, provides valuable genetic resources for maize breeding. To identify the favorable alien alleles in teosinte and its yield potential for maize breeding, 4 backcrossed maize-teosinte recombinant inbred line (RIL) populations were cultivated under five conditions. A North Carolina mating design II experiment was conducted on inbred lines with B73 and Mo17 pedigree backgrounds to analyze their combining ability. Abundant phenotypic variation on 26 traits of four RIL populations were found, of which barren tip length, kernel height, and test weight showed positive genetic improvement potential. The hybrid FM132 (BD138/MP116) showed a superior grain yield to that of the check, with an average yield gain of 4.86%. Moreover, inbred lines BD138 and MP048 showed a higher general grain yield combining ability than those of their corresponding checks. We screened 4,964,439 high-quality single-nucleotide polymorphisms in the BD (B73/Zea diploperennis) RIL population for bin construction and used 2322 bin markers for genetic map construction and quantitative trait loci (QTL) mapping. Via inclusive composite interval mapping, 71 QTL associated with 24 differential traits were identified. Gene annotation and transcriptional expression suggested that Zm00001eb352570 and Zm00001eb352580, both annotated as ethylene-responsive transcription factors, were key candidate genes that regulate ear height and the ratio of ear to plant height. Our results indicate that teosinte could broaden the narrow maize germplasm, improve yield potential, and provide desirable alleles for maize breeding.

High-quality chromosome-scale genomes facilitate effective identification of large structural variations in hot and sweet peppers

Pepper (Capsicum annuum) is an important vegetable crop that has been subjected to intensive breeding, resulting in limited genetic diversity, especially for sweet peppers. Previous studies have reported pepper draft genome assemblies using short read sequencing, but their capture of the extent of large structural variants (SVs), such as presence-absence variants (PAVs), inversions, and copy-number variants (CNVs) in the complex pepper genome falls short. In this study, we sequenced the genomes of representative sweet and hot pepper accessions by long-read and/or linked-read methods and advanced scaffolding technologies. First, we developed a high-quality reference genome for the sweet pepper cultivar 'Dempsey' and then used the reference genome to identify SVs in 11 other pepper accessions and constructed a graph-based pan-genome for pepper. We annotated an average of 42 972 gene families in each pepper accession, defining a set of 19 662 core and 23 115 non-core gene families. The new pepper pan-genome includes informative variants, 222 159 PAVs, 12 322 CNVs, and 16 032 inversions. Pan-genome analysis revealed PAVs associated with important agricultural traits, including potyvirus resistance, fruit color, pungency, and pepper fruit orientation. Comparatively, a large number of genes are affected by PAVs, which is positively correlated with the high frequency of transposable elements (TEs), indicating TEs play a key role in shaping the genomic landscape of peppers. The datasets presented herein provide a powerful new genomic resource for genetic analysis and genome-assisted breeding for pepper improvement.

Construction of a High-Density Recombination Bin-Based Genetic Map Facilitates High-Resolution Mapping of a Major QTL Underlying Anthocyanin Pigmentation in Eggplant

High-density genetic maps can significantly improve the resolution of QTL mapping. We constructed a high-density recombination bin-based genetic map of eggplant based on 200 F₂ plants from an interspecific cross (Solanum melongena × S. incanum) using the whole genome resequencing strategy. The map was 2022.8 cM long, covering near 99% of the eggplant genome. The map contained 3776 bins, with 3644 (96.5%) being effective (position non-redundant) ones, giving a nominal average distance of 0.54 cM and an effective average distance of 0.56 cM between adjacent bins, respectively. Using this map and 172 F_2:3 lines, a major QTL with pleiotropic effects on two anthocyanin pigmentation-related traits, leaf vein color (LVC) and fruit pericarp color (FPC), was steadily detected in a bin interval of 2.28 cM (or 1.68 Mb) on chromosome E10 in two cropping seasons, explaining ~65% and 55% of the phenotypic variation in LVC and FPC, respectively. Genome-wide association analysis in this population validated the QTL and demonstrated the correctness of mapping two bins of chromosome E02 onto E10. Bioinformatics analysis suggested that a WDR protein gene inside the bin interval with reliable effective variation between the two parents could be a possible candidate gene of the QTL.

Identification of Genetic Factors Controlling the Formation of Multiple Flowers Per Node in Pepper (Capsicum spp.)

Flower production provides the foundation for crop yield and increased profits. Capsicum annuum is a pepper species with a sympodial shoot structure with solitary flowers. By contrast, C. chinense produces multiple flowers per node. C. annuum accounts for 80% of pepper production worldwide. The identification of C. chinense genes that control multiple flowers and their transfer into C. annuum may open the way to increasing fruit yield. In this study, we dissected the genetic factors were dissected controlling the multiple-flower-per-node trait in Capsicum. 85 recombinant inbred lines (RILs) between the contrasting C. annuum 'TF68' and C. chinense 'Habanero' accessions were phenotyped and genotyped. Quantitative Trait Loci (QTL) analysis identified four novel QTLs on chromosomes 1, 2, 7, and 11 that accounted for 65% of the total phenotypic variation. Genome-wide association study was also performed on a panel of 276 genotyped and phenotyped C. annuum accessions, which revealed 28 regions significantly associated with the multiple-flower trait, of which three overlapped the identified QTLs. Five candidate genes involved in the development of the shoot and flower meristems were identified and these genes could cause multiple flowers per node in pepper. These results contribute to our understanding of multiple flower formation in Capsicum and will be useful to develop high-yielding cultivars.

Molecular Bases of Heat Stress Responses in Vegetable Crops With Focusing on Heat Shock Factors and Heat Shock Proteins

The effects of the climate change including an increase in the average global temperatures, and abnormal weather events such as frequent and severe heatwaves are emerging as a worldwide ecological concern due to their impacts on plant vegetation and crop productivity. In this review, the molecular processes of plants in response to heat stress-from the sensing of heat stress, the subsequent molecular cascades associated with the activation of heat shock factors and their primary targets (heat shock proteins), to the cellular responses-have been summarized with an emphasis on the classification and functions of heat shock proteins. Vegetables contain many essential vitamins, minerals, antioxidants, and fibers that provide many critical health benefits to humans. The adverse effects of heat stress on vegetable growth can be alleviated by developing vegetable crops with enhanced thermotolerance with the aid of various genetic tools. To achieve this goal, a solid understanding of the molecular and/or cellular mechanisms underlying various responses of vegetables to high temperature is imperative. Therefore, efforts to identify heat stress-responsive genes including those that code for heat shock factors and heat shock proteins, their functional roles in vegetable crops, and also their application to developing vegetables tolerant to heat stress are discussed.

Genotyping-by-sequencing-based QTL mapping reveals novel loci for Pepper yellow leaf curl virus (PepYLCV) resistance in Capsicum annuum

Disease caused by Pepper yellow leaf curl virus (PepYLCV) is one of the greatest threats to pepper (Capsicum spp.) cultivation in the tropics and subtropics. Resistance to PepYLCV was previously identified in a few Capsicum accessions, but no resistance QTLs have been mapped. This study aimed to elucidate the genetics of PepYLCV resistance in C. annuum L. Augmented inoculation by the viruliferous whitefly Bemisia tabaci was used to evaluate parental lines and an F₂ segregating population derived from a cross between resistant C. annuum line LP97 and susceptible C. annuum line ECW30R. Final evaluation was performed six weeks after inoculation using a standardized 5-point scale (0 = no symptoms to 4 = very severe symptoms). A high-density linkage map was constructed using genotyping-by-sequencing (GBS) to identify single-nucleotide polymorphism (SNP) markers associated with PepYLCV resistance in the F₂ population. QTL analysis revealed three QTLs, peplcv-1, peplcv-7, and peplcv-12, on chromosomes P1, P7, and P12, respectively. Candidate genes associated with PepYLCV resistance in the QTL regions were inferred. In addition, single markers Chr7-LCV-7 and Chr12-LCV-12 derived from the QTLs were developed and validated in another F₂ population and in commercial varieties. This work thus provides not only information for mapping PepYLCV resistance loci in pepper but also forms the basis for future molecular analysis of genes involved in PepYLCV resistance.

Identification of Fruit Traits Related QTLs and a Candidate Gene, CaBRX, Controlling Locule Number in Pepper (Capsicum annuum L.)

Fruit traits are important in pepper (Capsicum annuum L.) and affect its quality and yield. These traits are controlled by quantitative trait loci (QTLs). In this study, we identified many major QTLs that control fruit length (Ftl), fruit diameter (Ftd), fruit shape (Fts), fruit weight (Ftw) and locule number (Lcn) in the F2 and F2:3 populations developed from the QTL mapping of GS6 (P1) and Qiemen (P2). A total of 111 simple sequence repeats and insertion/deletion markers were utilized to construct a linkage map with 12 linkage groups over a length of 1320.72 cM. An inclusive composite interval mapping analysis indicated that many QTLs were detected and included ftl2.1, ftd2.1, fts1.1, ftw2.1 and lcn1.1. As a novel QTL, lcn1.1 was located between HM1112 and EPMS709, and the genetic distance was 3.18 cM covering 60 predicted genes. Within the region, we identified Capana01g004285 as a candidate gene by functional annotation and expression analysis and found that it encodes the BREVIS RADIX (BRX) protein. Knockdown of CaBRX through the virus-induced gene silencing approach in GS6 reduced the number of locules and influenced the expressions of genes related to flower and locule development, suggesting that CaBRX plays an important function in the development of locules.

Detection of Chromosomal Segments Introgressed from Wild Species of Carrot into Cultivars: Quantitative Trait Loci Mapping for Morphological Features in Backcross Inbred Lines

Cultivated carrot is thought to have been domesticated from a wild species, and various phenotypes developed through human domestication and selection over the past several centuries. Little is known about the genomic contribution of wild species to the phenotypes of present-day cultivars, although several studies have focused on identifying genetic loci that contribute to the morphology of storage roots. A backcross inbred line (BIL) population derived from a cross between the wild species Daucus carota ssp. carota "Songzi" and the orange cultivar "Amsterdam forcing" was developed. The morphological features in the BIL population became more diverse after several generations of selfing BC₂F₁ plants. Only few lines retained features of wild parent. Genomic resequencing of the two parental lines and the BILs resulted in 3,223,651 single nucleotide polymorphisms (SNPs), and 13,445 bin markers were generated using a sliding window approach. We constructed a genetic map with 2027 bins containing 154,776 SNPs; the total genetic distance was 1436.43 cM and the average interval between the bins was 0.71 cm. Five stable QTLs related to root length, root shoulder width, dry material content of root, and ratio of root shoulder width to root middle width were consistently detected on chromosome 2 in both years and explained 23.4-66.9% of the phenotypic variance. The effects of introgressed genomic segments from the wild species on the storage root are reported and will enable the identification of functional genes that control root morphological traits in carrot.

Pepper Fruit Elongation Is Controlled by Capsicum annuum Ovate Family Protein 20

Fruit shape is one of the most important quality traits of pepper (Capsicum spp.) and is used as a major attribute for the classification of fruit types. Wide natural variation in fruit shape exists among the major cultivated species Capsicum annuum, allowing the identification of several QTLs controlling the trait. However, to date, no genes underlying fruit shape QTLs have been conclusively identified, nor has their function been verified in pepper. We constructed a mapping population from a cross of round- and elongated-fruited C. annuum parents and identified a single major QTL on chromosome 10, termed fs10, explaining 68 and 70% of the phenotypic variation for fruit shape index and for distal fruit end angle, respectively. The QTL was mapped in several generations and was localized to a 5 Mbp region containing the ortholog of SlOFP20 that suppresses fruit elongation in tomato. Virus-induced gene silencing of the pepper ortholog CaOFP20 resulted in increased fruit elongation on two independent backgrounds. Furthermore, CaOFP20 exhibited differential expression in fs10 near-isogenic lines, as well as in an association panel of elongated- and round-fruited accessions. A 42-bp deletion in the upstream region of CaOFP20 was most strongly associated with fruit shape variation within the locus. Histological observations in ovaries and fruit pericarps indicated that fs10 exerts its effect on fruit elongation by controlling cell expansion and replication. Our results indicate that CaOFP20 functions as a suppressor of fruit elongation in C. annuum and is the most likely candidate gene underlying fs10.

SNP-Based Genetic Linkage Map and Quantitative Trait Locus Mapping Associated with the Agronomically Important Traits of Hypsizygus marmoreus

White strains of Hypsizygus marmoreus are more difficult to cultivate than are brown strains; therefore, new white strain breeding strategies are required. Accordingly, we constructed the genetic map of H. marmoreus with 1996 SNP markers on 11 linkage groups (LGs) spanning 1380.49 cM. Prior to analysis, 82 backcrossed strains (HM8 lines) were generated by mating between KMCC03106-31 and the progenies of the F1 hybrid (Hami-18 × KMCC03106-93). Using HM8, the first 23 quantitative trait loci (QTLs) of yield-related traits were detected with high limit of detection (LOD) scores (1.98-9.86). The length, thickness, and hardness of the stipe were colocated on LG 1. Especially, length of stipe and thickness of stipe were highly correlated given that the correlation coefficients were negative (-0.39, p value ≤ .01). And a typical biomodal distribution was observed for lightness of the pileus and the lightness of the pileus trait belonged to the LG 8, as did traits of earliness and mycelial growth in potato dextrose agar (PDA) medium. Therefore, results for color traits can be suggested that color is controlled by a multi-gene of one locus. The yield trait was highly negatively correlated with the traits for thickness of the stipe (-0.45, p value ≤ .01). Based on additive effects, the white strain was confirmed as recessive; however, traits of mycelial growth, lightness, and quality were inherited by backcrossed HM8 lines. This new genetic map, finely mapped QTLs, and the strong selection markers could be used in molecular breeding of H. marmoreus.

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

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

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

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

A mutation in Zeaxanthin epoxidase contributes to orange coloration and alters carotenoid contents in pepper fruit (Capsicum annuum)

Phytoene synthase (PSY1), capsanthin-capsorubin synthase (CCS), and pseudo-response regulator 2 (PRR2) are three major genes controlling fruit color in pepper (Capsicum spp.). However, the diversity of fruit color in pepper cannot be completely explained by these three genes. Here, we used an F₂ population derived from Capsicum annuum 'SNU-mini Orange' (SO) and C. annuum 'SNU-mini Yellow' (SY), both harboring functional PSY1 and mutated CCS, and observed that yellow color was dominant over orange color. We performed genotyping-by-sequencing and mapped the genetic locus to a 6.8-Mb region on chromosome 2, which we named CaOr. We discovered a splicing mutation in the zeaxanthin epoxidase (ZEP) gene within this region leading to a premature stop codon. HPLC analysis showed that SO contained higher amounts of zeaxanthin and total carotenoids in mature fruits than SY. A color complementation assay using Escherichia coli harboring carotenoid biosynthetic genes showed that the mutant ZEP protein had reduced enzymatic activity. Transmission electron microscopy of plastids revealed that the ZEP mutation affected plastid development with more rod-shaped inner membrane structures in chromoplasts of mature SO fruits. To validate the role of ZEP in fruit color formation, we performed virus-induced gene silencing of ZEP in the yellow-fruit cultivar C. annuum 'Micropep Yellow' (MY). The silencing of ZEP caused significant changes in the ratios of zeaxanthin to its downstream products and increased total carotenoid contents. Thus, we conclude that the ZEP genotype can determine orange or yellow mature fruit color in pepper.

RNA-Seq-based high-resolution linkage map reveals the genetic architecture of fruiting body development in shiitake mushroom, Lentinula edodes

•

We constructed a reference genetic map of Lentinula edodes.

•

We re-assembled a chromosome-level genome of L. edodes.

•

We disclosed three hotspots regions for fruiting body-related traits in shiitake.

•

We scanned candidate genes for fruiting body-related traits.

Fruiting body development (FBD) of mushroom-forming fungi has attracted tremendous interest. However, the genetic and molecular basis of FBD is poorly known. Here, using Lentinula edodes (shiitake) as a model, we deciphered the genetic architecture underlying fruiting body-related traits (FBRTs) by combined genomic, genetic and phenotypic data. Using RNA-Seq of fruiting bodies from 110 dikaryons in a bi-parental mapping population, we constructed an ultra-high-density genetic map of L. edodes (Lemap2.0) with a total length of 810.14 cM, which covered 81.7% of the shiitake genome. A total of 94 scaffolds of the shiitake genome were aligned to Lemap2.0 and re-anchored into nine pseudo-chromosomes. Then via quantitative trait locus (QTL) analysis, we disclosed an outline of the genetic architecture of FBD in shiitake. Twenty-nine QTLs and three main genomic regions associated with FBD of shiitake were identified. Using meta-QTL analysis, seven pleiotropic QTLs for multiple traits were detected, which contributed to the correlations of FBRTs. In the mapped QTLs, the expression of 246 genes were found to significantly correlate with the phenotypic traits. Thirty-three of them were involved in FBD and could represent candidate genes controlling the shape and size of fruiting bodies. Collectively, our findings have advanced our understanding of the genetic regulation of FBD in shiitake and mushroom-forming fungi at large.

Mapping QTLs for 1000-grain weight and genes controlling hull type using SNP marker in Tartary buckwheat (Fagopyrum tataricum)

BACKGROUND

Tartary buckwheat (Fagopyrum tataricum), an important pseudocereal crop, has high economic value due to its nutritional and medicinal properties. However, dehulling of Tartary buckwheat is difficult owing to its thick and tough hull, which has greatly limited the development of the Tartary buckwheat processing industry. The construction of high-resolution genetic maps serves as a basis for identifying quantitative trait loci (QTLs) and qualitative trait genes for agronomic traits. In this study, a recombinant inbred lines (XJ-RILs) population derived from a cross between the easily dehulled Rice-Tartary type and Tartary buckwheat type was genotyped using restriction site-associated DNA (RAD) sequencing to construct a high-density SNP genetic map. Furthermore, QTLs for 1000-grain weight (TGW) and genes controlling hull type were mapped in multiple environments.

RESULTS

In total, 4151 bin markers comprising 122,185 SNPs were used to construct the genetic linkage map. The map consisted of 8 linkage groups and covered 1444.15 cM, with an average distance of 0.35 cM between adjacent bin markers. Nine QTLs for TGW were detected and distributed on four loci on chromosome 1 and 4. A major locus detected in all three trials was mapped in 38.2-39.8 cM region on chromosome 1, with an LOD score of 18.1-37.0, and explained for 23.6-47.5% of the phenotypic variation. The genes controlling hull type were mapped to chromosome 1 between marker Block330 and Block331, which was closely followed by the major locus for TGW. The expression levels of the seven candidate genes controlling hull type present in the region between Block330 and Block336 was low during grain development, and no significant difference was observed between the parental lines. Six non-synonymous coding SNPs were found between the two parents in the region.

CONCLUSIONS

We constructed a high-density SNP genetic map for the first time in Tartary buckwheat. The mapped major loci controlling TGW and hull type will be valuable for gene cloning and revealing the mechanism underlying grain development and easy dehulling, and marker-assisted selection in Tartary buckwheat.

Chromosome Level Assembly of Homozygous Inbred Line 'Wongyo 3115' Facilitates the Construction of a High-Density Linkage Map and Identification of QTLs Associated With Fruit Firmness in Octoploid Strawberry (Fragaria × ananassa)

Strawberry is an allo-octoploid crop with high genome heterozygosity and complexity, which hinders the sequencing and the assembly of the genome. However, in the present study, we have generated a chromosome level assembly of octoploid strawberry sourced from a highly homozygous inbred line 'Wongyo 3115', using long- and short-read sequencing technologies. The assembly of 'Wongyo 3115' produced 805.6 Mb of the genome with 323 contigs scaffolded into 208 scaffolds with an N50 of 27.3 Mb after further gap filling. The whole genome annotation resulted in 151,892 genes with a gene density of 188.52 (genes/Mb) and validation of a genome, using BUSCO analysis resulted in 94.10% complete BUSCOs. Firmness is one of the vital traits in strawberry, which facilitate the postharvest shelf-life qualities. The molecular and genetic mechanisms that contribute the firmness in strawberry remain unclear. We have constructed a high-density genetic map based on the 'Wongyo 3115' reference genome to identify loci associated with firmness in the present study. For the quantitative trait locus (QTL) identification, the 'BS F₂' populations developed from two inbred lines were genotyped, using an Axiom 35K strawberry chip, and marker positions were analyzed based on the 'Wongyo 3115' genome. Genetic maps were constructed with 1,049 bin markers, spanning the 3,861 cM. Using firmness data of 'BS F₂' obtained from 2 consecutive years, five QTLs were identified on chromosomes 3-3, 5-1, 6-1, and 6-4. Furthermore, we predicted the candidate genes associated with firmness in strawberries by utilizing transcriptome data and QTL information. Overall, we present the chromosome-level assembly and annotation of a homozygous octoploid strawberry inbred line and a linkage map constructed to identify QTLs associated with fruit firmness.

Bin-based genome-wide association studies reveal superior alleles for improvement of appearance quality using a 4-way MAGIC population in rice

•

4-way Multiparental population covered the limitations of the biparental structure.

•

The combination of SNP and bin-GWAS showed a powerful tool for QTL mapping.

•

qPGWC8.2 harbored a novel predicted gene for rice chalkiness quality.

Introduction

The multiparental population provides us the chance to identify superior alleles controlling a trait for genetic improvement. Genome wide association studies at bin level (bin-GWAS) are expected to be more power in QTL mapping than GWAS at SNP level (SNP-GWAS).

Objectives

This study is to estimate genetic effects of QTL conferring grain appearance quality in rice by SNP-GWAS and bin-GWAS, compare their power in QTL mapping and identify the superior alleles of all detected QTL from 4 parents for genetic improvement.

Methods

A 4-way MAGIC population and its four founders were cultivated in two environments to dissect the genetic basis of rice grain appearance quality. Both SNP-GWAS and bin-GWAS were conducted for QTL mapping. Multiple comparison among 4 parental bin/alleles was used to identify the superior alleles.

Results

A total of 16 and 20 QTL associated with grain appearance quality were identified by SNP- and bin-GWAS, respectively. A minor chalkiness QTL qPGWC8.2/qDEC8 was assigned to a 30-kb genomic region, in which OsMH_08T0121900 is the potential candidate gene because its encoded protein, glucan endo-1,3-beta-glucosidase precursor is involved in the starch and sucrose metabolism pathway. The superior parental alleles for GS3, GL3.1, GW5, GW7, and Chalk5 and two QTLs were almost carried by the high-quality parents Cypress and Yuejingsimiao (YJSM), while the poor-quality parent Guichao-2 (GC2) always carried the inferior alleles. The top five recombinant inbred lines with the highest quality of grain shape and chalkiness traits all carried gene combinations of superior alleles.

Conclusions

Both SNP- and bin-GWAS methods are encouraged for joint QTL mapping with MAGIC population. qPGWC8.2/qDEC8 is a novel candidate gene strongly associated with chalkiness. The superior alleles of GS3, GW5, GL3.1, GW7, Chalk5 and qPGWC8.2 were identified, and the pyramiding of these superior alleles is helpful to improve rice appearance quality.

Genetic dissection of an allotetraploid interspecific CSSLs guides interspecific genetics and breeding in cotton

Background

The low genetic diversity of Upland cotton limits the potential for genetic improvement. Making full use of the genetic resources of Sea-island cotton will facilitate genetic improvement of widely cultivated Upland cotton varieties. The chromosome segments substitution lines (CSSLs) provide an ideal strategy for mapping quantitative trait loci (QTL) in interspecific hybridization.

Results

In this study, a CSSL population was developed by PCR-based markers assisted selection (MAS), derived from the crossing and backcrossing of Gossypium hirsutum (Gh) and G. barbadense (Gb), firstly. Then, by whole genome re-sequencing, 11,653,661 high-quality single nucleotide polymorphisms (SNPs) were identified which ultimately constructed 1211 recombination chromosome introgression segments from Gb. The sequencing-based physical map provided more accurate introgressions than the PCR-based markers. By exploiting CSSLs with mutant morphological traits, the genes responding for leaf shape and fuzz-less mutation in the Gb were identified. Based on a high-resolution recombination bin map to uncover genetic loci determining the phenotypic variance between Gh and Gb, 64 QTLs were identified for 14 agronomic traits with an interval length of 158 kb to 27 Mb. Surprisingly, multiple alleles of Gb showed extremely high value in enhancing cottonseed oil content (SOC).

Conclusions

This study provides guidance for studying interspecific inheritance, especially breeding researchers, for future studies using the traditional PCR-based molecular markers and high-throughput re-sequencing technology in the study of CSSLs. Available resources include candidate position for controlling cotton quality and quantitative traits, and excellent breeding materials. Collectively, our results provide insights into the genetic effects of Gb alleles on the Gh, and provide guidance for the utilization of Gb alleles in interspecific breeding.

Whole genome resequencing of four Italian sweet pepper landraces provides insights on sequence variation in genes of agronomic value

Sweet pepper (Capsicum annuum L.) is a high value crop and one of the most widely grown vegetables belonging to the Solanaceae family. In addition to commercial varieties and F₁ hybrids, a multitude of landraces are grown, whose genetic combination is the result of hundreds of years of random, environmental, and farmer selection. High genetic diversity exists in the landrace gene pool which however has scarcely been studied, thus bounding their cultivation. We re-sequenced four pepper inbred lines, within as many Italian landraces, which representative of as many fruit types: big sized blocky with sunken apex ('Quadrato') and protruding apex or heart shaped ('Cuneo'), elongated ('Corno') and smaller sized sub-spherical ('Tumaticot'). Each genomic sequence was obtained through Illumina platform at coverage ranging from 39 to 44×, and reconstructed at a chromosome scale. About 35.5k genes were predicted in each inbred line, of which 22,017 were shared among them and the reference genome (accession 'CM334'). Distinctive variations in miRNAs, resistance gene analogues (RGAs) and susceptibility genes (S-genes) were detected. A detailed survey of the SNP/Indels occurring in genes affecting fruit size, shape and quality identified the highest frequencies of variation in regulatory regions. Many structural variations were identified as presence/absence variations (PAVs), notably in resistance gene analogues (RGAs) and in the capsanthin/capsorubin synthase (CCS) gene. The large allelic diversity observed in the four inbred lines suggests their potential use as a pre-breeding resource and represents a one-stop resource for C. annuum genomics and a key tool for dissecting the path from sequence variation to phenotype.

Multi-dimensional machine learning approaches for fruit shape phenotyping in strawberry

BACKGROUND

Shape is a critical element of the visual appeal of strawberry fruit and is influenced by both genetic and non-genetic determinants. Current fruit phenotyping approaches for external characteristics in strawberry often rely on the human eye to make categorical assessments. However, fruit shape is an inherently multi-dimensional, continuously variable trait and not adequately described by a single categorical or quantitative feature. Morphometric approaches enable the study of complex, multi-dimensional forms but are often abstract and difficult to interpret. In this study, we developed a mathematical approach for transforming fruit shape classifications from digital images onto an ordinal scale called the Principal Progression of k Clusters (PPKC). We use these human-recognizable shape categories to select quantitative features extracted from multiple morphometric analyses that are best fit for genetic dissection and analysis.

RESULTS

We transformed images of strawberry fruit into human-recognizable categories using unsupervised machine learning, discovered 4 principal shape categories, and inferred progression using PPKC. We extracted 68 quantitative features from digital images of strawberries using a suite of morphometric analyses and multivariate statistical approaches. These analyses defined informative feature sets that effectively captured quantitative differences between shape classes. Classification accuracy ranged from 68% to 99% for the newly created phenotypic variables for describing a shape.

CONCLUSIONS

Our results demonstrated that strawberry fruit shapes could be robustly quantified, accurately classified, and empirically ordered using image analyses, machine learning, and PPKC. We generated a dictionary of quantitative traits for studying and predicting shape classes and identifying genetic factors underlying phenotypic variability for fruit shape in strawberry. The methods and approaches that we applied in strawberry should apply to other fruits, vegetables, and specialty crops.

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

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

Genetic Dissection of Germinability under Low Temperature by Building a Resequencing Linkage Map in japonica Rice

Among all cereals, rice is highly sensitive to cold stress, especially at the germination stage, which adversely impacts its germination ability, seed vigor, crop stand establishment, and, ultimately, grain yield. The dissection of novel quantitative trait loci (QTLs) or genes conferring a low-temperature germination (LTG) ability can significantly accelerate cold-tolerant rice breeding to ensure the wide application of rice cultivation through the direct seeding method. In this study, we identified 11 QTLs for LTG using 144 recombinant inbred lines (RILs) derived from a cross between a cold-tolerant variety, Lijiangxintuanheigu (LTH), and a cold-sensitive variety, Shennong265 (SN265). By resequencing two parents and RIL lines, a high-density bin map, including 2,828 bin markers, was constructed using 123,859 single-nucleotide polymorphisms (SNPs) between two parents. The total genetic distance corresponding to all 12 chromosome linkage maps was 2,840.12 cm. Adjacent markers were marked by an average genetic distance of 1.01 cm, corresponding to a 128.80 kb physical distance. Eight and three QTL alleles had positive effects inherited from LTH and SN265, respectively. Moreover, a pleiotropic QTL was identified for a higher number of erected panicles and a higher grain number on Chr-9 near the previously cloned DEP1 gene. Among the LTG QTLs, qLTG3 and qLTG7b were also located at relatively small genetic intervals that define two known LTG genes, qLTG3-1 and OsSAP16. Sequencing comparisons between the two parents demonstrated that LTH possesses qLTG3-1 and OsSAP16 genes, and SN-265 owns the DEP1 gene. These comparison results strengthen the accuracy and mapping resolution power of the bin map and population. Later, fine mapping was done for qLTG6 at 45.80 kb through four key homozygous recombinant lines derived from a population with 1569 segregating plants. Finally, LOC_Os06g01320 was identified as the most possible candidate gene for qLTG6, which contains a missense mutation and a 32-bp deletion/insertion at the promoter between the two parents. LTH was observed to have lower expression levels in comparison with SN265 and was commonly detected at low temperatures. In conclusion, these results strengthen our understanding of the impacts of cold temperature stress on seed vigor and germination abilities and help improve the mechanisms of rice breeding programs to breed cold-tolerant varieties.

A Resequencing-Based Ultradense Genetic Map of Hericium erinaceus for Anchoring Genome Sequences and Identifying Genetic Loci Associated With Monokaryon Growth

Hericium erinaceus has attracted tremendous interest owing to its compelling health-promoting properties. However, breeding of elite cultivars of H. erinaceus is hindered by the lack of a genetic and molecular toolbox. Here, we performed resequencing analysis of 127 F1 single-spore isolates and constructed the first high-resolution genetic map of H. erinaceus. With the use of recombination bins as markers, an ultradense genetic map consisting of 1,174 bins (including 37,082 single-nucleotide polymorphisms) was generated. This newly developed genetic map covered 1,096.5 cM, with an average bin spacing of 0.95 cM. High collinearity between genetic map and H. erinaceus genome assembly was revealed by aligning scaffolds to this genetic map using bin markers as anchors. The application of this newly developed genetic map in quantitative trait locus (QTL) mapping was also elucidated, and four QTLs for monokaryon growth were recovered. One QTL, mgr1, which contributes 12.1% of growth variations, was located near the mating type A (MAT-A) loci. Overall, this newly constructed high-resolution genetic map (or bin map) could be used as reference in future genetic, genomic, and breeding studies on H. erinaceus.

Genomic Selection for Prediction of Fruit-Related Traits in Pepper (Capsicum spp.)

Pepper (Capsicum spp.) fruit-related traits are critical determinants of quality. These traits are controlled by quantitatively inherited genes for which marker-assisted selection (MAS) has proven insufficiently effective. Here, we evaluated the potential of genomic selection, in which genotype and phenotype data for a training population are used to predict phenotypes of a test population with only genotype data, for predicting fruit-related traits in pepper. We measured five fruit traits (fruit length, fruit shape, fruit width, fruit weight, and pericarp thickness) in 351 accessions from the pepper core collection, including 229 Capsicum annuum, 48 Capsicum baccatum, 48 Capsicum chinense, 25 Capsicum frutescens, and 1 Capsicum chacoense in 4 years at two different locations and genotyped these accessions using genotyping-by-sequencing. Among the whole core collection, considering its genetic distance and sexual incompatibility, we only included 302 C. annum complex (229 C. annuum, 48 C. chinense, and 25 C. frutescens) into further analysis. We used phenotypic and genotypic data to investigate genomic prediction models, marker density, and effects of population structure. Among 10 genomic prediction methods tested, Reproducing Kernel Hilbert Space (RKHS) produced the highest prediction accuracies (measured as correlation between predicted values and observed values) across the traits, with accuracies of 0.75, 0.73, 0.84, 0.83, and 0.82 for fruit length, fruit shape, fruit width, fruit weight, and pericarp thickness, respectively. Overall, prediction accuracies were positively correlated with the number of markers for fruit traits. We tested our genomic selection models in a separate population of recombinant inbred lines derived from two parental lines from the core collection. Despite the large difference in genetic diversity between the training population and the test population, we obtained moderate prediction accuracies of 0.32, 0.34, 0.50, and 0.48 for fruit length, fruit shape, fruit width, and fruit weight, respectively. This use of genomic selection for fruit-related traits demonstrates the potential use of core collections and genomic selection as tools for crop improvement.

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.

All-in-one sequencing: an improved library preparation method for cost-effective and high-throughput next-generation sequencing

BACKGROUND

Next generation sequencing (NGS) has been widely used in biological research, due to its rapid decrease in cost and increasing ability to generate data. However, while the sequence generation step has seen many improvements over time, the library preparation step has not, resulting in low-efficiency library preparation methods, especially for the most time-consuming and labor-intensive steps: size-selection and quantification. Consequently, there can be bottlenecks in projects with large sample cohorts.

RESULTS

We have described the all-in-one sequencing (AIO-seq) method, where instead of performing size-selection and quantification for samples individually, one sample one tube, up to 116 samples are pooled and analyzed in a single tube, 'All-In-One'. The AIO-seq method pools libraries based on the samples' expected data yields and the calculated concentrations of the size selected regions (target region), which can easily be obtained with the Agilent 2100 Bioanalyzer and Qubit Fluorometer. AIO-seq was applied to whole genome sequencing and RNA-seq libraries successfully, and it is envisaged that it could be applied to any type of NGS library, such as chromatin immunoprecipitation coupled with massively parallel sequencing, assays for transposase-accessible chromatin with high-throughput sequencing, and high-throughput chromosome conformation capture. We also demonstrated that for genetic population samples with low coverage sequences, like recombinant inbred lines (RIL), AIO-seq could be further simplified, by mixing the libraries immediately after PCR, without calculating the target region concentrations.

CONCLUSIONS

The AIO-seq method is thus labor saving and cost effective, and suitable for projects with large sample cohorts, like those used in plant breeding or population genetics research.

Target sequencing reveals genetic diversity, population structure, core-SNP markers, and fruit shape-associated loci in pepper varieties

BACKGROUND

The widely cultivated pepper (Capsicum spp.) is one of the most diverse vegetables; however, little research has focused on characterizing the genetic diversity and relatedness of commercial varieties grown in China. In this study, a panel of 92 perfect single-nucleotide polymorphisms (SNPs) was identified using re-sequencing data from 35 different C. annuum lines. Based on this panel, a Target SNP-seq genotyping method was designed, which combined multiplex amplification of perfect SNPs with Illumina sequencing, to detect polymorphisms across 271 commercial pepper varieties.

RESULTS

The perfect SNPs panel had a high discriminating capacity due to the average value of polymorphism information content, observed heterozygosity, expected heterozygosity, and minor allele frequency, which were 0.31, 0.28, 0.4, and 0.31, respectively. Notably, the studied pepper varieties were morphologically categorized based on fruit shape as blocky-, long horn-, short horn-, and linear-fruited. The long horn-fruited population exhibited the most genetic diversity followed by the short horn-, linear-, and blocky-fruited populations. A set of 35 core SNPs were then used as kompetitive allele-specific PCR (KASPar) markers, another robust genotyping technique for variety identification. Analysis of genetic relatedness using principal component analysis and phylogenetic tree construction indicated that the four fruit shape populations clustered separately with limited overlaps. Based on STRUCTURE clustering, it was possible to divide the varieties into five subpopulations, which correlated with fruit shape. Further, the subpopulations were statistically different according to a randomization test and F_st statistics. Nine loci, located on chromosomes 1, 2, 3, 4, 6, and 12, were identified to be significantly associated with the fruit shape index (p < 0.0001).

CONCLUSIONS

Target SNP-seq developed in this study appears as an efficient power tool to detect the genetic diversity, population relatedness and molecular breeding in pepper. Moreover, this study demonstrates that the genetic structure of Chinese pepper varieties is significantly influenced by breeding programs focused on fruit shape.

A recombination bin-map identified a major QTL for resistance to Tomato Spotted Wilt Virus in peanut (Arachis hypogaea)

Tomato spotted wilt virus (TSWV) is a devastating disease to peanut growers in the South-eastern region of the United States. Newly released peanut cultivars in recent years are crucial as they have some levels of resistance to TSWV. One mapping population of recombinant inbred line (RIL) used in this study was derived from peanut lines of SunOleic 97R and NC94022. A whole genome re-sequencing approach was used to sequence these two parents and 140 RILs. A recombination bin-based genetic map was constructed, with 5,816 bins and 20 linkage groups covering a total length of 2004 cM. Using this map, we identified three QTLs which were colocalized on chromosome A01. One QTL had the largest effect of 36.51% to the phenotypic variation and encompassed 89.5 Kb genomic region. This genome region had a cluster of genes, which code for chitinases, strictosidine synthase-like, and NBS-LRR proteins. SNPs linked to this QTL were used to develop Kompetitive allele specific PCR (KASP) markers, and the validated KASP markers showed expected segregation of alleles coming from resistant and susceptible parents within the population. Therefore, this bin-map and QTL associated with TSWV resistance made it possible for functional gene mapping, map-based cloning, and marker-assisted breeding. This study identified the highest number of SNP makers and demonstrated recombination bin-based map for QTL identification in peanut. The chitinase gene clusters and NBS-LRR disease resistance genes in this region suggest the possible involvement in peanut resistance to TSWV.

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

There are many agronomic traits of pepper (Capsicum L.) with abundant phenotypes that can benefit pepper growth. Using specific-locus amplified fragment sequencing (SLAF-seq), a genome-wide association study (GWAS) of 36 agronomic traits was carried out for 287 representative pepper accessions. To ensure the accuracy and reliability of the GWAS results, we analyzed the genetic diversity, distribution of labels (SLAF tags and single nucleotide polymorphisms (SNPs)) and population differentiation and determined the optimal statistical model. In our study, 1487 SNPs were highly significantly associated with 26 agronomic traits, and 2126 candidate genes were detected in the 100-kb region up- and down-stream near these SNPs. Furthermore, 13 major association peaks were identified for 11 key agronomic traits. Then we examined the correlations among the 36 agronomic traits and analyzed SNP distribution and found 37 SNP polymerization regions (total size: 264.69 Mbp) that could be selected areas in pepper breeding. We found that the stronger the correlation between the two traits, the greater the possibility of them being in more than one polymerization region, suggesting that they may be linked or that one pleiotropic gene controls them. These results provide a theoretical foundation for future multi-trait pyramid breeding of pepper. Finally, we found that the GWAS signals were highly consistent with those from the nuclear restorer-of-fertility (Rf) gene for cytoplasmic male sterility (CMS), verifying their reliability. We further identified Capana06g002967 and Capana06g002969 as Rf candidate genes by functional annotation and expression analysis, which provided a reference for the study of cytoplasmic male sterility in Capsicum.

Genomic diversity and novel genome-wide association with fruit morphology in Capsicum, from 746k polymorphic sites

Capsicum is one of the major vegetable crops grown worldwide. Current subdivision in clades and species is based on morphological traits and coarse sets of genetic markers. Broad variability of fruits has been driven by breeding programs and has been mainly studied by linkage analysis. We discovered 746k variable sites by sequencing 1.8% of the genome in a collection of 373 accessions belonging to 11 Capsicum species from 51 countries. We describe genomic variation at population-level, confirm major subdivision in clades and species, and show that the known major subdivision of C. annuum separates large and bulky fruits from small ones. In C. annuum, we identify four novel loci associated with phenotypes determining the fruit shape, including a non-synonymous mutation in the gene Longifolia 1-like (CA03g16080). Our collection covers all the economically important species of Capsicum widely used in breeding programs and represent the widest and largest study so far in terms of the number of species and number of genetic variants analyzed. We identified a large set of markers that can be used for population genetic studies and genetic association analyses. Our results provide a comprehensive and precise perspective on genomic variability in Capsicum at population-level and suggest that future fine genetic association studies will yield useful results for breeding.

Identifying candidate genes for Phytophthora capsici resistance in pepper (Capsicum annuum) via genotyping-by-sequencing-based QTL mapping and genome-wide association study

Phytophthora capsici (Leon.) is a globally prevalent, devastating oomycete pathogen that causes root rot in pepper (Capsicum annuum). Several studies have identified quantitative trait loci (QTL) underlying resistance to P. capsici root rot (PcRR). However, breeding for pepper cultivars resistant to PcRR remains challenging due to the complexity of PcRR resistance. Here, we combined traditional QTL mapping with GWAS to broaden our understanding of PcRR resistance in pepper. Three major-effect loci (5.1, 5.2, and 5.3) conferring broad-spectrum resistance to three isolates of P. capsici were mapped to pepper chromosome P5. In addition, QTLs with epistatic interactions and minor effects specific to isolate and environment were detected on other chromosomes. GWAS detected 117 significant SNPs across the genome associated with PcRR resistance, including SNPs on chromosomes P5, P7, and P11 that colocalized with the QTLs identified here and in previous studies. Clusters of candidate nucleotide-binding site-leucine-rich repeat (NBS-LRR) and receptor-like kinase (RLK) genes were predicted within the QTL and GWAS regions; such genes often function in disease resistance. These candidate genes lay the foundation for the molecular dissection of PcRR resistance. SNP markers associated with QTLs for PcRR resistance will be useful for marker-assisted breeding and genomic selection in pepper breeding.

Genetic diversity and population structure of Ethiopian Capsicum germplasms

We established a collection of 142 Capsicum genotypes from different geographical areas of Ethiopia with the aim of capturing genetic diversity. Morphological traits and high-resolution melting analysis distinguished one Capsicum baccatum, nine Capsicum frutescens and 132 Capsicum annuum accessions in the collection. Measurement of plant growth parameters revealed variation between germplasms in parameters including plant height, stem thickness, internode length, number of side branches, fruit width, and fruit length. Broad-sense heritability was maximum for fruit weight, followed by length and width of leaves. We used genotyping by sequencing (GBS) to identify single-nucleotide polymorphisms (SNPs) in the panel of 142 Capsicum germplasms and found 2,831,791 genome-wide SNP markers. Among these, we selected 53,284 high-quality SNPs and used them to estimate the level of genetic diversity, population structure, and phylogenetic relationships. From model-based ancestry analysis, the phylogenetic tree and principal-coordinate analysis (PCoA), we identified two distinct genetic populations: one comprising 132 C. annuum accessions and the other comprising the nine C. frutescens accessions. GWAS analysis detected 509 SNP markers that were significantly associated with fruit-, stem- and leaf-related traits. This is the first comprehensive report of the analysis of genetic variation in Ethiopian Capsicum species involving a large number of accessions. The results will help breeders utilize the germplasm collection to improve existing commercial cultivars.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

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

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

Sequencing-Based Bin Map Construction of a Tomato Mapping Population, Facilitating High-Resolution Quantitative Trait Loci Detection

Genotyping-by-sequencing (GBS) was employed to construct a highly saturated genetic linkage map of a tomato ( L.) recombinant inbred line (RIL) population, derived from a cross between cultivar NC EBR-1 and the wild tomato L. accession LA2093. A pipeline was developed to convert single nucleotide polymorphism (SNP) data into genomic bins, which could be used for fine mapping of quantitative trait loci (QTL) and identification of candidate genes. The pipeline, implemented in a python script named SNPbinner, adopts a hidden Markov model approach for calculation of recombination breakpoints followed by genomic bins construction. The total length of the newly developed high-resolution genetic map was 1.2-fold larger than previously estimated based on restriction fragment length polymorphism (RFLP) and polymerase chain reaction (PCR)-based markers. The map was used to verify and refine QTL previously identified for two fruit quality traits in the RIL population, fruit weight (FW) and fruit lycopene content (LYC). Two well-described FW QTL ( and ) were localized precisely at their known underlying causative genes, and the QTL intervals were decreased by two- to tenfold. A major QTL for LYC content () was verified at high resolution and its underlying causative gene was determined to be ζ (). The RIL population, the high resolution genetic map, and the easy-to-use genotyping pipeline, SNPbinner, are made publicly available.

A major QTL and candidate genes for capsaicinoid biosynthesis in the pericarp of Capsicum chinense revealed using QTL-seq and RNA-seq

A major QTL and candidate genes controlling capsaicinoid content in the pericarp were identified by QTL-seq and RNA-seq in Capsicum chinense. Capsaicinoid biosynthesis was previously thought to be restricted to the placental tissue; however, the recent discovery of their biosynthesis in the pericarp provides new opportunities to increase the capsaicinoid content in pepper fruits. Currently, the genetic mechanisms regulating capsaicinoid biosynthesis in the pericarp remain unknown. Here, we performed quantitative trait loci (QTL) mapping and RNA sequencing (RNA-seq) to reveal the genes controlling capsaicinoid biosynthesis in the pericarp. A whole-genome sequencing-based QTL-seq strategy was employed, identifying a major QTL on chromosome 6. To validate the QTL on chromosome 6, we performed traditional QTL mapping using the same population in QTL-seq with an additional biparental population. A total of 15 QTLs for capsaicinoid content distributed on chromosomes 3, 6, and 11 were newly identified. Among these QTLs, the genetic loci on the lower arm of chromosome 6 were commonly detected in the two mapping populations, corresponding to the location of the major QTL detected using whole-genome sequencing-based QTL-seq. Our RNA-seq analysis identified candidate genes within the common QTL that were differentially expressed in the pungent and non-pungent pericarp tissues. Our results are expected to contribute to the elucidation of the regulation of capsaicinoid biosynthesis. We also demonstrated that a combination of QTL mapping and RNA-seq is helpful for refining the candidate genes of a complicated trait of interest.

Construction of a high density genetic map of an interspecific cross of Capsicum chinense and Capsicum annuum and QTL analysis of floral traits

The yield of pepper plants (Capsicum spp.) is their most important trait and is affected by the flower number and flowering time. Capsicum annuum produces a single flower per node and has an early flowering habit. By contrast, Capsicum chinense yields multiple flowers per node and has a late flowering character. However, the genetic mechanism underlying the control of these floral traits remains largely unknown. In this study, 150 F₂ populations from an interspecific cross between the inbred lines 740 (C. chinense) and CA1 (C. annuum) and their parents were used to construct a molecular genetic linkage map using the specific length amplified fragment sequencing (SLAF-seq) technique. This linkage map, spanning 1,586.78 cM in length, contained 9,038 markers on 12 chromosomes, with a mean marker distance of 0.18 cM. Phenotypic data on the flowering time and flower number per node were collected over multiple years, and QTL analysis identified 6 QTLs for the flowering time and flower number per node by composite interval mapping (CIM) and genome-wide composite interval mapping (GCIM) methods at least in two environments. The candidate genes within the major QTL were predicted. In the major flowering time QTL, the candidate gene Capana02g000700, which encodes the homeotic protein APETALA2, was identified. Quantitative reverse-transcription PCR (qRT-PCR) analysis indicated that its expression level in 740 was higher than that in CA1. Gene expression analysis indicated that the expression of Capana02g000700 was significantly upregulated in flowers, and many floral development-related genes were found to be coexpressed with Capana02g000700, supporting the function of this gene in association with flowering time in C. chinense and C. annuum species.

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

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