Genetic characterization of a locus responsible for low pungency using EMS-induced mutants in Capsicum annuum L

KEY MESSAGE

The pepper mutants ('221-2-1a' and '1559-1-2h') with very low pungency were genetically characterized. The Pun4 locus, responsible for the reduced pungency of the mutant fruits, was localized to a 208 Mb region on chromosome 6. DEMF06G16460, encoding 3-ketoacyl-CoA synthase, was proposed as a strong candidate gene based on the genetic analyses of bulked segregants, DEG, and expression analyses. Capsaicinoids are unique alkaloids present in pepper (Capsicum spp.), synthesized through the condensation of by-products from the phenylpropanoid and branched-chain fatty acid pathways, and accumulating in the placenta. In this study, we characterized two allelic ethyl methanesulfonate-induced mutant lines with extremely low pungency ('221-2-1a' and '1559-1-2h'). These mutants, derived from the pungent Korean landrace 'Yuwolcho,' exhibited lower capsaicinoid content than Yuwolcho but still contained a small amount of capsaicinoid with functional capsaicinoid biosynthetic genes. Genetic crosses between the mutants and Yuwolcho or pungent lines indicated that a single recessive mutation was responsible for the low-pungency phenotype of mutant 221-2-1a; we named the causal locus Pungency 4 (Pun4). To identify Pun4, we combined genome-wide polymorphism analysis and transcriptome analysis with bulked-segregant analysis. We narrowed down the location of Pun4 to a 208-Mb region on chromosome 6 containing five candidate genes, of which DEMF06G16460, encoding a 3-ketoacyl-CoA synthase associated with branched-chain fatty acid biosynthesis, is the most likely candidate for Pun4. The expression of capsaicinoid biosynthetic genes in placental tissues in Yuwolcho and the mutant was consistent with the branched-chain fatty acid pathway playing a pivotal role in the lower pungency observed in the mutant. We also obtained a list of differentially expressed genes in placental tissues between the mutant and Yuwolcho, from which we selected candidate genes using gene co-expression analysis. In summary, we characterized the capsaicinoid biosynthesis-related locus Pun4 through integrated of genetic, genomic, and transcriptome analyses. These findings will contribute to our understanding of capsaicinoid biosynthesis in pepper.

Molecular mapping of the broad bean wilt virus 2 resistance locus bwvr in Capsicum annuum using BSR-seq

KEY MESSAGE

Bulked segregant RNA seq of pools of pepper accessions that are susceptible or resistant to Broad bean wilt virus 2 identifies a gene that might confer resistance to this devastating pathogen. The single-stranded positive-sense RNA virus Broad bean wilt virus 2 (BBWV2) causes substantial damage to pepper (Capsicum annuum) cultivation. Here, we describe mapping the BBWV2 resistance locus bwvr using a F7:8 recombinant inbred line (RIL) population constructed by crossing the BBWV2-resistant pepper accession 'SNU-C' with the susceptible pepper accession 'ECW30R.' All F1 plants infected with the BBWV2 strain PAP1 were susceptible to the virus, and the RIL population showed a 1:1 ratio of resistance to susceptibility, indicating that this trait is controlled by a single recessive gene. To map bwvr, we performed bulked segregant RNA-seq (BSR-seq). We sequenced pools of resistant and susceptible lines from the RILs and aligned the reads to the high-quality 'Dempsey' reference genome to identify variants between the pools. This analysis identified 519,887 variants and selected the region from 245.9-250.8 Mb of the Dempsey reference genome as the quantitative trait locus region for bwvr. To finely map bwvr, we used newly designed high-resolution melting (HRM) and Kompetitive allele specific PCR (KASP) markers based on variants obtained from the BSR-seq reads and the PepperSNP16K array. Comparative analysis identified 11 SNU-C-specific SNPs within the bwvr locus. Using markers derived from these variants, we mapped the candidate bwvr locus to the region from 246.833-246.949 kb. SNU-C-specific variants clustered near DEM.v1.00035533 within the bwvr locus. DEM.v1.00035533 encodes the nitrate transporter NPF1.2 and contains a SNP within its 5' untranslated region. The bwvr locus, which contains four genes including DEM.v1.00035533, could represent a valuable resource for global pepper breeding programs.

Leaky mutations in the zeaxanthin epoxidase in Capsicum annuum result in bright-red fruit containing a high amount of zeaxanthin

Fruit color is one of the most important traits in peppers due to its esthetic value and nutritional benefits and is determined by carotenoid composition, resulting from diverse mutations of carotenoid biosynthetic genes. The EMS204 line, derived from an EMS mutant population, presents bright-red color, compared with the wild type Yuwolcho cultivar. HPLC analysis indicates that EMS204 fruit contains more zeaxanthin and less capsanthin and capsorubin than Yuwolcho. MutMap was used to reveal the color variation of EMS204 using an F3 population derived from a cross of EMS204 and Yuwolcho, and the locus was mapped to a 2.5-Mbp region on chromosome 2. Among the genes in the region, a missense mutation was found in ZEP (zeaxanthin epoxidase) that results in an amino acid sequence alteration (V291 → I). A color complementation experiment with Escherichia coli and ZEP in vitro assay using thylakoid membranes revealed decreased enzymatic activity of EMS204 ZEP. Analysis of endogenous plant hormones revealed a significant reduction in abscisic acid content in EMS204. Germination assays and salinity stress experiments corroborated the lower ABA levels in the seeds. Virus-induced gene silencing showed that ZEP silencing also results in bright-red fruit containing less capsanthin but more zeaxanthin than control. A germplasm survey of red color accessions revealed no similar carotenoid profiles to EMS204. However, a breeding line containing a ZEP mutation showed a very similar carotenoid profile to EMS204. Our results provide a novel breeding strategy to develop red pepper cultivars containing high zeaxanthin contents using ZEP mutations.

Development of virus-induced genome editing methods in Solanaceous crops

Genome editing (GE) using CRISPR/Cas systems has revolutionized plant mutagenesis. However, conventional transgene-mediated GE methods have limitations due to the time-consuming generation of stable transgenic lines expressing the Cas9/single guide RNA (sgRNA) module through tissue cultures. Virus-induced genome editing (VIGE) systems have been successfully employed in model plants, such as Arabidopsis thaliana and Nicotiana spp. In this study, we developed two VIGE methods for Solanaceous plants. First, we used the tobacco rattle virus (TRV) vector to deliver sgRNAs into a transgenic tomato (Solanum lycopersicum) line of cultivar Micro-Tom expressing Cas9. Second, we devised a transgene-free GE method based on a potato virus X (PVX) vector to deliver Cas9 and sgRNAs. We designed and cloned sgRNAs targeting Phytoene desaturase in the VIGE vectors and determined optimal conditions for VIGE. We evaluated VIGE efficiency through deep sequencing of the target gene after viral vector inoculation, detecting 40.3% and 36.5% mutation rates for TRV- and PVX-mediated GE, respectively. To improve editing efficiency, we applied a 37°C heat treatment, which increased the editing efficiency by 33% to 46% and 56% to 76% for TRV- and PVX-mediated VIGE, respectively. To obtain edited plants, we subjected inoculated cotyledons to tissue culture, yielding successful editing events. We also demonstrated that PVX-mediated GE can be applied to other Solanaceous crops, such as potato (Solanum tuberosum) and eggplant (Solanum melongena). These simple and highly efficient VIGE methods have great potential for generating genome-edited plants in Solanaceous crops.

Exploring horticultural traits and disease resistance in Capsicum baccatum through segmental introgression lines

KEY MESSAGE

Segmental introgression and advanced backcross lines were developed and validated as important tools for improving agronomically important traits in pepper, offering improved sensitivity in detecting quantitative trait loci for breeding. Segmental introgression lines (SILs) and advanced backcross lines (ABs) can accelerate genetics and genomics research and breeding in crop plants. This study presents the development of a complete collection of SILs and ABs in pepper using Capsicum annuum cv. 'CM334' as the recipient parent and Capsicum baccatum 'PBC81', which displays various agronomically important traits including powdery mildew and anthracnose resistance, as donor parent. Using embryo rescue to overcome abortion in interspecific crosses, and marker-assisted selection with genotyping-in-thousands by sequencing (GT-seq) to develop SILs and ABs containing different segments of the C. baccatum genome, we obtained 63 SILs and 44 ABs, covering 94.8% of the C. baccatum genome. We characterized them for traits including powdery mildew resistance, anthracnose resistance, anthocyanin accumulation, trichome density, plant architecture, and fruit morphology. We validated previously known loci for these traits and discovered new sources of variation and quantitative trait loci (QTLs). A total of 15 QTLs were identified, including four for anthracnose resistance with three novel loci, seven for plant architecture, and four for fruit morphology. This is the first complete collection of pepper SILs and ABs validated for agronomic traits and will enhance QTL detection and serve as valuable breeding resources. Further, these SILs and ABs will be useful for comparative genomics and to better understand the genetic mechanisms underlying important agronomic traits in pepper, ultimately leading to improved crop productivity and sustainability.

Development of a speed breeding protocol with flowering gene investigation in pepper (Capsicum annuum)

Pepper (Capsicum spp.) is a vegetable and spice crop in the Solanaceae family with many nutritional benefits for human health. During several decades, horticultural traits, including disease resistance, yield, and fruit quality, have been improved through conventional breeding methods. Nevertheless, cultivar development is a time-consuming process because of the long generation time of pepper. Recently, speed breeding has been introduced as a solution for shorting the breeding cycle in long-day or day-neutral field crops, but there have been only a few studies on speed breeding in vegetable crops. In this study, a speed breeding protocol for pepper was developed by controlling the photoperiod and light quality. Under the condition of a low red (R) to far-red (FR) ratio of 0.3 with an extended photoperiod (Epp) of 20 h (95 ± 0 DAT), the time to first harvest was shortened by 75 days after transplant (DAT) compared to that of the control treatment (170 ± 2 DAT), suggesting that Epp with FR light is an essential factor for flowering in pepper. In addition, we established the speed breeding system in a greenhouse with a 20 h photoperiod and a 3.8 R:FR ratio and promoted the breeding cycle of C. annuum for 110 days from seed to seed. To explain the accelerated flowering response to the Epp and supplemented FR light, genome-wide association study (GWAS) and gene expression analysis were performed. As a result of the GWAS, we identified a new flowering gene locus for pepper and suggested four candidate genes for flowering (APETALA2 (AP2), WUSCHEL-RELATED HOMEOBOX4 (WOX4), FLOWERING LOCUS T (FT), and GIGANTEA (GI)). Through expression analysis with the candidate genes, it appeared that Epp and FR induced flowering by up-regulating the flowering-promoting gene GI and down-regulating FT. The results demonstrate the effect of a combination of Epp and FR light by genetic analysis of flowering gene expression. This is the first study that verifies gene expression patterns associated with the flowering responses of pepper in a speed breeding system. Overall, this study demonstrates that speed breeding can shorten the breeding cycle and accelerate genetic research in pepper through reduced generation time.

Genome-wide association study of resistance to anthracnose in pepper (Capsicum chinense) germplasm

BACKGROUND

Anthracnose is a fungal disease caused by Colletotrichum spp. that has a significant impact on worldwide pepper production. Colletotrichum scovillei is the most common pathogenic anthracnose-causing species in the Republic of Korea.

RESULTS

The resistances of 197 pepper (Capsicum chinense) accessions deposited in Korea's National Agrobiodiversity Center were evaluated for their response against the virulent pathogens Colletotrichum acutatum isolate 'KSCa-1' and C. scovillei isolate 'Hana') in the field and in vitro methods for three consecutive years (2018 to 2020). The severity of the disease was recorded and compared between inoculation methods. Six phenotypically resistant pepper accessions were selected based on three years of disease data. All of the selected resistant pepper accessions outperformed the control resistant pepper in terms of resistance (PI 594,137). A genome-wide association study (GWAS) was carried out to identify single nucleotide polymorphisms (SNPs) associated with anthracnose resistance. An association analysis was performed using 53,518 SNPs and the disease score of the 2020 field and in vitro experiment results. Both field and in vitro experiments revealed 25 and 32 significantly associated SNPs, respectively. These SNPs were found on all chromosomes except Ch06 and Ch07 in the field experiment, whereas in the in vitro experiment they were found on all chromosomes except Ch04 and Ch11.

CONCLUSION

In this study, six resistant C. chinense accessions were selected. Additionally, in this study, significantly associated SNPs were found in a gene that codes for a protein kinase receptor, such as serine/threonine-protein kinase, and other genes that are known to be involved in disease resistance. This may strengthen the role of these genes in the development of anthracnose resistance in Capsicum spp. As a result, the SNPs discovered to be strongly linked in this study can be used to identify a potential marker for selecting pepper material resistant to anthracnose, which will assist in the development of resistant varieties.

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.

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.

Investigation of genetic factors regulating chlorophyll and carotenoid biosynthesis in red pepper fruit

Chlorophylls and carotenoids are synthesized in the chloroplast and chromoplast, respectively. Even though the two pigments are generated from the same precursor, the genetic correlation between chlorophyll and carotenoid biosynthesis has not yet been fully understood. We investigated the genetic correlation of chlorophyll and carotenoid biosynthesis during fruit ripening. Two recombinant inbred lines populations, "Long Sweet" × "AC2212" ("LA") RILs derived from a cross between Capsicum annuum "Long Sweet" with light-green and light-red fruit and C. annuum "AC2212" with dark-green and brown-fruit and "3501 (F)" × "3509 (C)" ("FC") RILs from C. annuum "3501" with dark-green and dark-red fruit and C. annuum "3509" with intermediate green and light-red fruit, were used. As the fruit ripened, three accessions produced high levels of xanthophyll. The dark-green immature fruit accumulated more total carotenoids than the light-green fruit. This trend corresponded to the expression pattern of 1-deoxy-d-xylulose 5-phosphate synthase (DXS) and CaGLK2 genes during fruit development. The expression levels of DXS and CaGLK2 in the dark-green accession "3501" were significantly higher than those of "3509" and "Long Sweet" during the early stages of fruit development. Furthermore, the genotype analysis of the transcription factor controlling chloroplast development (CaGLK2) in LA RILs revealed that CaGLK2 expression affected both carotenoid and chlorophyll contents. The single nucleotide polymorphism (SNP) linkage maps were constructed using genotyping-by-sequencing (GBS) for the two populations, and QTL analysis was performed for green fruit color intensity and carotenoid content. The QTL (LA_BG-CST10) for capsanthin content in LA RILs located at 24.4 to 100.4 Mbp on chromosome 10 was overlapped with the QTL (FC15-Cap10) for capsanthin content in FC RILs. Three QTLs for capsanthin content, American spice trade association (ASTA) value, and immature green fruit color intensity were also overlapped from 178.2 to 204 Mbp on chromosome 10. At the location, 151.6 to 165 Mbp on chromosome 8, QTLs (FC15-tcar8, FC17-ASTA8.1, and FC17-ASTA8.2) for total carotenoid content and ASTA value were discovered, and this region contained 2-C-methyl-d-erythritol 4-phosphate cytidylyltransferase (MCT), which is involved in the MEP pathway. This result is the first report to show the correlation between carotenoid and chlorophyll biosynthesis in pepper. This research will expand our understanding of the mechanism of the chloroplast-to-chromoplast transition and the development of high pigment pepper varieties.

Genomic selection with fixed-effect markers improves the prediction accuracy for Capsaicinoid contents in Capsicum annuum

Capsaicinoids provide chili peppers (Capsicum spp.) with their characteristic pungency. Several structural and transcription factor genes are known to control capsaicinoid contents in pepper. However, many other genes also regulating capsaicinoid contents remain unknown, making it difficult to develop pepper cultivars with different levels of capsaicinoids. Genomic selection (GS) uses genome-wide random markers (including many in undiscovered genes) for a trait to improve selection efficiency. In this study, we predicted the capsaicinoid contents of pepper breeding lines using several GS models trained with genotypic and phenotypic data from a training population. We used a core collection of 351 Capsicum accessions and 96 breeding lines as training and testing populations, respectively. To obtain the optimal number of single nucleotide polymorphism (SNP) markers for GS, we tested various numbers of genome-wide SNP markers based on linkage disequilibrium. We obtained the highest mean prediction accuracy (0.550) for different models using 3294 SNP markers. Using this marker set, we conducted GWAS and selected 25 markers that were associated with capsaicinoid biosynthesis genes and quantitative trait loci for capsaicinoid contents. Finally, to develop more accurate prediction models, we obtained SNP markers from GWAS as fixed-effect markers for GS, where 3294 genome-wide SNPs were employed. When four to five fixed-effect markers from GWAS were used as fixed effects, the RKHS and RR-BLUP models showed accuracies of 0.696 and 0.689, respectively. Our results lay the foundation for developing pepper cultivars with various capsaicinoid levels using GS for capsaicinoid contents.

Identification of CaAN3 as a fruit-specific regulator of anthocyanin biosynthesis in pepper (Capsicum annuum)

The novel gene CaAN3 encodes an R2R3 MYB transcription factor that regulates fruit-specific anthocyanin accumulation. The key regulatory gene CaAN2 encodes an R2R3 MYB transcription factor that regulates anthocyanin biosynthesis in various tissues in pepper (Capsicum annuum). However, CaAN2 is not expressed in certain pepper accessions showing fruit-specific anthocyanin accumulation. In this study, we identified the novel locus CaAN3 as a regulator of fruit-specific anthocyanin biosynthesis, using an F₂ population derived from a hybrid cultivar with purple immature fruits and segregating for CaAN3. We extracted total RNA, assembled two RNA pools according to fruit color, and carried out bulked segregant RNA sequencing. We aligned the raw reads to the pepper reference genome Dempsey and identified 6,672 significant single nucleotide polymorphisms (SNPs) by calculating the Δ(SNP-index) between the two pools. We then conducted molecular mapping to delimit the target region of CaAN3 to the interval 184.6-186.4 Mbp on chromosome 10. We focused on Dem.v1.00043895, encoding an R2R3 MYB transcription factor, as the strongest candidate gene. Sequence analysis revealed four insertion/deletion polymorphisms in the promoter region of the green CaAN3 allele. We employed virus-induced gene silencing and transient overexpression assays to characterize the function of the candidate gene. When Dem.v1.00043895 was silenced in pepper, anthocyanin accumulation decreased in the pericarp, while the transient overexpression of Dem.v1.00043895 in Nicotiana benthamiana leaves resulted in the accumulation of anthocyanins around the infiltration sites. These results showed that Dem.v1.00043895 is CaAN3, an activator of anthocyanin biosynthesis in pepper fruits.

The Unstable Restorer-of-fertility locus in pepper (Capsicum annuum. L) is delimited to a genomic region containing PPR genes

Unstable Restorer-of-fertility (Rf_u), conferring unstable fertility restoration in the pepper CGMS system, was delimited to a genomic region near Rf and is syntenic to the PPR-like gene-rich region in tomato. The use of cytoplasmic-genic male sterility (CGMS) systems greatly increases the efficiency of hybrid seed production. Although marker development and candidate gene isolation have been performed for the Restorer-of-fertility (Rf) gene in pepper (Capsicum annuum L.), the broad use of CGMS systems has been hampered by the instability of fertility restoration among pepper accessions, especially sweet peppers, due to the widespread presence of the Unstable Restorer-of-fertility (Rf_u) locus. Therefore, to investigate the genetic factors controlling unstable fertility restoration in sweet peppers, we developed a segregation population (BC₄F₅) from crosses using a male-sterile line and an Rf_u-containing line. Segregation did not significantly deviate from a 3:1 ratio for unstable fertility restoration to sterility, indicating single dominant locus control for unstable fertility restoration in this population. Genetic mapping delimited the Rf_u locus to a 398 kb genomic region on chromosome 6, which is close to but different from the previously identified Rf-containing region. The Rf_u-containing region harbors a pentatricopeptide repeat (PPR) gene, along with 10 other candidate genes. In addition, this region is syntenic to the genomic region containing the largest number of Rf-like PPR genes in tomato. Therefore, the dynamic evolution of PPR genes might be responsible for both the restoration and instability of fertility in pepper. During genetic mapping, we developed various molecular markers, including one that co-segregated with Rf_u. These markers showed higher accuracy for genotyping than previously developed markers, pointing to their possible use in marker-assisted breeding of sweet peppers.

Genome-Wide Association Study of Resistance to Phytophthora capsici in the Pepper (Capsicum spp.) Collection

One of the most serious pepper diseases is Phytophthora blight, which is caused by Phytophthora capsici. It is crucial to assess the resistance of pepper genetic resources to Phytophthora blight, understand the genetic resistances, and develop markers for selecting resistant pepper materials in breeding programs. In this study, the resistance of 342 pepper accessions to P. capsici was evaluated. The disease severity score method was used to evaluate the phenotypic responses of pepper accessions inoculated with the KCP7 isolate. A genome-wide association study (GWAS) was performed to identify single nucleotide polymorphisms (SNPs) linked to P. capsici (isolate KCP7) resistance. The pepper population was genotyped using the genotype-by-sequencing (GBS) method, and 45,481 SNPs were obtained. A GWAS analysis was performed using resistance evaluation data and SNP markers. Significantly associated SNPs for P. capsici resistance at 4 weeks after inoculation of the GWAS pepper population were selected. These SNPs for Phytophthora blight resistance were found on all chromosomes except Chr.05, Chr.09, and Chr.11. One of the SNPs found on Chr.02 was converted into a high-resolution melting (HRM) marker, and another marker (QTL5-1) from the previous study was applied to pepper accessions and breeding lines for validation and comparison. This SNP marker was selected because the resistance phenotype and the HRM marker genotype matched well. The selected SNP was named Chr02-1126 and was located at 112 Mb on Chr.02. The Chr02-1126 marker predicted P. capsici resistance with 78.5% accuracy, while the QTL5-1 marker predicted resistance with 80.2% accuracy. Along with the marker for major quantitative traits loci (QTLs) on Chr.05, this Chr02-1126 marker could be used to accurately predict Phytophthora blight resistance in pepper genetic resources. Therefore, this study will assist in the selection of resistant pepper plants in order to breed new phytophthora blight-resistant varieties.

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.

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 an Allelic Variant of the CsOr Gene Controlling Fruit Endocarp Color in Cucumber (Cucumis sativus L.) Using Genotyping-By-Sequencing (GBS) and Whole-Genome Sequencing

The cucumber is a major vegetable crop around the world. Fruit flesh color is an important quality trait in cucumber and flesh color mainly depends on the relative content of β-carotene in the fruits. The β-carotene serves as a precursor of vitamin A, which has dietary benefits for human health. Cucumbers with orange flesh contain a higher amount of β-carotene than white fruit flesh. Therefore, development of orange-fleshed cucumber varieties is gaining attention for improved nutritional benefits. In this study, we performed genotyping-by-sequencing (GBS) based on genetic mapping and whole-genome sequencing to identify the orange endocarp color gene in the cucumber breeding line, CS-B. Genetic mapping, genetic sequencing, and genetic segregation analyses showed that a single recessive gene (CsaV3_6G040750) encodes a chaperone DnaJ protein (DnaJ) protein at the Cucumis sativus(CsOr) locus was responsible for the orange endocarp phenotype in the CS-B line. The Or gene harbored point mutations T13G and T17C in the first exon of the coding region, resulting in serine to alanine at position 13 and isoleucine to threonine at position 17, respectively. CS-B line displayed increased β-carotene content in the endocarp tissue, corresponding to elevated expression of CsOr gene at fruit developmental stages. Identifying novel missense mutations in the CsOr gene could provide new insights into the role of Or mechanism of action for orange fruit flesh in cucumber and serve as a valuable resource for developing β-carotene-rich cucumbers varieties with increased nutritional benefits.

Comparative Genomic Analysis Reveals Genetic Variation and Adaptive Evolution in the Pathogenicity-Related Genes of Phytophthora capsici

Phytophthora capsici is an oomycete pathogen responsible for damping off, root rot, fruit rot, and foliar blight in popular vegetable and legume crops. The existence of distinct aggressiveness levels and physiological races among the P. capsici population is a major constraint to developing resistant varieties of host crops. In the present study, we compared the genomes of three P. capsici isolates with different aggressiveness levels to reveal their genomic differences. We obtained genome sequences using short-read and long-read technologies, which yielded an average genome size of 76 Mbp comprising 514 contigs and 15,076 predicted genes. A comparative genomic analysis uncovered the signatures of accelerated evolution, gene family expansions in the pathogenicity-related genes among the three isolates. Resequencing two additional P. capsici isolates enabled the identification of average 1,023,437 SNPs, revealing the frequent accumulation of non-synonymous substitutions in pathogenicity-related gene families. Furthermore, pathogenicity-related gene families, cytoplasmic effectors and ATP binding cassette (ABC) transporters, showed expansion signals in the more aggressive isolates, with a greater number of non-synonymous SNPs. This genomic information explains the plasticity, difference in aggressiveness levels, and genome structural variation among the P. capsici isolates, providing insight into the genomic features related to the evolution and pathogenicity of this oomycete pathogen.

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.

Identification of the determinant of tomato yellow leaf curl Kanchanaburi virus infectivity in tomato

Geminiviruses cause devastating diseases in solanaceous crops, with the bipartite begomoviruses tomato yellow leaf curl Kanchanaburi virus (TYLCKaV) and pepper yellow leaf curl Thailand virus (PYLCThV) major threats in Southeast Asia. To determine the molecular mechanism of geminivirus infection, we constructed infectious clones of TYLCKaV and PYLCThV. Both constructs infected Nicotiana benthamiana, but only TYLCKaV could infect Solanum lycopersicum 'A39'. A genome-swapping of TYLCKaV with PYLCThV revealed the TYLCKaV-B genome segment as the determinant of TYLCKaV infectivity in tomato. We constructed five geminivirus clones with chimeric TYLCKaV-B and PYLCThV-B genome segments to narrow down the region determining TYLCKaV infectivity in tomato. Only chimeric clones carrying the TYLCKaV intergenic region (IR) showed infectivity in S. lycopersicum 'A39', indicating that the IR of TYLCKaV-B is essential for TYLCKaV infectivity in tomato. Our results provide a foundation for elucidating the molecular mechanism of geminivirus infection in plants.

Development of a Panel of Genotyping-in-Thousands by Sequencing in Capsicum

Genotyping by sequencing (GBS) enables genotyping of multiple loci at low cost. However, the single nucleotide polymorphisms (SNPs) revealed by GBS tend to be randomly distributed between individuals, limiting their direct comparisons without applying the various filter options to obtain a comparable dataset of SNPs. Here, we developed a panel of a multiplex targeted sequencing method, genotyping-in-thousands by sequencing (GT-seq), to genotype SNPs in Capsicum spp. Previously developed Fluidigm^® SNP markers were converted to GT-seq markers and combined with new GT-seq markers developed using SNP information obtained through GBS. We then optimized multiplex PCR conditions: we obtained the highest genotyping rate when the first PCR consisted of 25 cycles. In addition, we determined that 101 primer pairs performed best when amplifying target sequences of 79 bp. We minimized interference of multiplex PCR by primer dimer formation using the PrimerPooler program. Using our GT-seq pipeline on Illumina Miseq and Nextseq platforms, we genotyped up to 1,500 (Miseq) and 1,300 (Nextseq) samples for the optimum panel size of 100 loci. To allow the genotyping of Capsicum species, we designed 332 informative GT-seq markers from Fluidigm SNP markers and GBS-derived SNPs. This study illustrates the first application of GT-seq in crop plants. The GT-seq marker set developed here will be a useful tool for molecular breeding of peppers in the future.

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.

Population Structure and Genetic Diversity of Cucurbita moschata Based on Genome-Wide High-Quality SNPs

Pumpkins (Cucurbita moschata) are one of the most important economic crops in genus Cucurbita worldwide. They are a popular food resource and an important rootstock resource for various Cucurbitaceae. Especially, C. moschata is widely used as a rootstock for the commercial production of bloomless cucumbers in East Asia. Since the genetic diversity of the commercial rootstock varieties is narrow, there has been an increasing demand for the trait development of abiotic and biotic stress tolerance breeding. In this study, 2071 high-quality SNPs that were distributed evenly across 20 chromosomes of pumpkins were discovered through the genotyping-by-sequencing (GBS) analysis of 610 accessions of C. moschata germplasm with a global origin. Using these SNPs, various analyses of the genetic diversity and the population structure were performed. Three subgroups were clustered from the germplasm collection, which included East Asia, Africa, and America, and these areas were included the most in each subgroup. Among those groups, accessions from Africa and South Asia showed the highest genetic diversity, which was followed by the Mexico accessions. This result reflected that large gene pools that consist of various native landraces have been conserved in those of countries. Based on the genetic diversity, we finally constructed the C. moschata core collection, which included 67 representative accessions from the 610 germplasms. Five morphological traits that are important in commercial grafting and rootstock seed production, which include the cotyledon length, the cotyledon width, the hypocotyl length, the internode length, and the number of female flowers, were investigated for three years and used to confirm the validity of the core collection selection. The results are expected to provide valuable information about the genetic structure of the worldwide C. moschata germplasm and help to create new gene pools to develop genetically diverse rootstock breeding materials.

Jasmonic acid and ERF family genes are involved in chilling sensitivity and seed browning of pepper fruit after harvest

Pepper (Capsicum annuum L.) fruit is sensitive to temperatures below 10 °C, which severely diminish fruit quality during cold chain distribution. Seed browning was a major chilling symptom in 36 genotypes of C. annuum fruit screened after storage at 2 °C for 3 weeks. Among them, pepper fruits of chilling-insensitive 'UZB-GJG-1999-51' and -sensitive 'C00562' were treated at 2 °C for 0 or 24 h, respectively. Analyses of integrated transcriptome-metabolome and relative gene expression in seeds obtained from the two genotypes were conducted to identify key factors involved in the seed browning induced by chilling. The relative contents of branched-chain amino acids such as leucine, isoleucine, and valine were significantly increased after chilling. Transcriptome identification showed 3,140 differentially expressed genes (log twofold change > 1.0 and FDR-corrected p value < 0.05) affected by chilling between the two genotypes. Particularly, genes related to jasmonic acid synthesis and signaling were differentially expressed. A regulatory network of jasmonic acid synthesis and signaling, and regulation of ERF family genes might contribute to chilling response in pepper fruit. The results of this study may help facilitate further studies to develop chilling-insensitive peppers and could be a basis for improving postharvest fruit quality.

F-Box Family Genes, LTSF1 and LTSF2, Regulate Low-Temperature Stress Tolerance in Pepper (Capsicum chinense)

The F-box proteins belong to a family of regulatory proteins that play key roles in the proteasomal degradation of other proteins. Plant F-box proteins are functionally diverse, and the precise roles of many such proteins in growth and development are not known. Previously, two low-temperature-sensitive F-box protein family genes (LTSF1 and LTSF2) were identified as candidates responsible for the sensitivity to low temperatures in the pepper (Capsicum chinense) cultivar ‘sy-2’. In the present study, we showed that the virus-induced gene silencing of these genes stunted plant growth and caused abnormal leaf development under low-temperature conditions, similar to what was observed in the low-temperature-sensitive ‘sy-2’ line. Protein–protein interaction analyses revealed that the LTSF1 and LTSF2 proteins interacted with S-phase kinase-associated protein 1 (SKP1), part of the Skp, Cullin, F-box-containing (SCF) complex that catalyzes the ubiquitination of proteins for degradation, suggesting a role for LTSF1 and LTSF2 in protein degradation. Furthermore, transgenic Nicotiana benthamiana plants overexpressing the pepper LTSF1 gene showed an increased tolerance to low-temperature stress and a higher expression of the genes encoding antioxidant enzymes. Taken together, these results suggest that the LTSF1 and LTSF2 F-box proteins are a functional component of the SCF complex and may positively regulate low-temperature stress tolerance by activating antioxidant-enzyme activities.

Phytoene synthase 2 can compensate for the absence of PSY1 in the control of color in Capsicum fruit

Phytoene synthase 1 (PSY1) and capsanthin-capsorubin synthase (CCS) are two major genes responsible for fruit color variation in pepper (Capsicum spp.). However, the role of PSY2 remains unknown. We used a systemic approach to examine the genetic factors responsible for the yellow fruit color of C. annuum 'MicroPep Yellow' (MY) and to determine the role of PSY2 in fruit color. We detected complete deletion of PSY1 and a retrotransposon insertion in CCS. Despite the loss of PSY1 and CCS function, both MY and mutant F2 plants from a cross between MY and the 'MicroPep Red' (MR) accumulated basal levels of carotenoids, indicating that other PSY genes may complement the loss of PSY1. qRT-PCR analysis indicated that PSY2 was constitutively expressed in both MR and MY fruits, and a color complementation assay using Escherichia coli revealed that PSY2 was capable of biosynthesizing a carotenoid. Virus-induced gene silencing of PSY2 in MY resulted in white fruits. These findings indicate that PSY2 can compensate for the absence of PSY1 in pepper fruit, resulting in the yellow color of MY fruits.

Transcriptome analysis to identify candidate genes associated with the yellow-leaf phenotype of a Cymbidium mutant generated by γ-irradiation

Leaf color is an important agronomic trait in flowering plants, including orchids. However, factors underlying leaf phenotypes in plants remain largely unclear. A mutant displaying yellow leaves was obtained by the γ-ray-based mutagenesis of a Cymbidium orchid and characterized using RNA sequencing. A total of 144,918 unigenes obtained from over 25 million reads were assigned to 22 metabolic pathways in the Kyoto Encyclopedia of Genes and Genomes database. In addition, gene ontology was used to classify the predicted functions of transcripts into 73 functional groups. The RNA sequencing analysis identified 2,267 differentially expressed genes between wild-type and mutant Cymbidium sp. Genes involved in the chlorophyll biosynthesis and degradation, as well as ion transport, were identified and assayed for their expression levels in wild-type and mutant plants using quantitative real-time profiling. No critical expression changes were detected in genes involved in chlorophyll biosynthesis. In contrast, seven genes involved in ion transport, including two metal ion transporters, were down-regulated, and chlorophyllase 2, associated with chlorophyll degradation, was up-regulated. Together, these results suggest that alterations in chlorophyll metabolism and/or ion transport might contribute to leaf color in Cymbidium orchids.

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.

Genome Editing of eIF4E1 in Tomato Confers Resistance to Pepper Mottle Virus

Many of the recessive virus-resistance genes in plants encode eukaryotic translation initiation factors (eIFs), including eIF4E, eIF4G, and related proteins. Notably, eIF4E and its isoform eIF(iso)4E are pivotal for viral infection and act as recessive resistance genes against various potyviruses in a wide range of plants. In this study, we used Clustered Regularly Interspaced Palindromic Repeats/CRISPR-associated protein 9 (CRISPR/Cas9)-mediated targeted mutagenesis to test whether novel sequence-specific mutations at eIF4E1 in Solanum lycopersicum (tomato) cv. Micro-Tom could confer enhanced resistance to potyviruses. This approach produced heritable homozygous mutations in the transgene-free E₁ generation. Sequence analysis of eIF4E1 from E₀ transgenic plants expressing Cas9 and eIF4E-sgRNA transcripts identified chimeric deletions ranging from 11 to 43 bp. Genotype analysis of the eIF4E1-edited lines in E₀, E₁, and E₂ transgenic tomato plants showed that the mutations were transmitted to subsequent generations. When homozygous mutant lines were tested for resistance to potyviruses, they exhibited no resistance to tobacco etch virus (TEV). Notably, however, several mutant lines showed no accumulation of viral particles upon infection with pepper mottle virus (PepMoV). These results indicate that site-specific mutation of tomato eIF4E1 successfully conferred enhanced resistance to PepMoV. Thus, this study demonstrates the feasibility of the use of CRISPR/Cas9 approach to accelerate breeding for trait improvement in tomato plants.

Candidate Gene Analysis Reveals That the Fruit Color Locus C1 Corresponds to PRR2 in Pepper (Capsicum frutescens)

The diverse fruit colors of peppers (Capsicum spp.) are due to variations in carotenoid composition and content. Mature fruit color in peppers is regulated by three independent loci, C1, C2, and Y. C2 and Y encode phytoene synthase (PSY1) and capsanthin-capsorubin synthase (CCS), respectively; however, the identity of the C1 gene has been unknown. With the aim of identifying C1, we analyzed two pepper accessions with different fruit colors: Capsicum frutescens AC08-045 and AC08-201, whose fruits are light yellow and white, respectively. Ultra-performance liquid chromatography showed that the total carotenoid content was six times higher in AC08-045 than in AC08-201 fruits, with similar composition of main carotenoids and slight difference in minor components. These results suggest that a genetic factor in AC08-201 may down-regulate overall carotenoid biosynthesis. Analyses of candidate genes related to carotenoid biosynthesis and plastid abundance revealed that both accessions carry non-functional alleles of CCS, golden2-like transcription factor (GLK2), and PSY1. However, a nonsense mutation (C2571T) in PRR2, a homolog of Arabidopsis pseudo response regulator2-like (APRR2), was present in only AC08-201. In a population derived from a cross between AC08-045 and AC08-201, a SNP marker based on the nonsense mutation co-segregated fully with fruit color, implying that the mutation in PRR2 may cause the white color of AC08-201 fruits. Transmission electron microscopy (TEM) of AC08-201 fruit pericarp also showed less developed granum structure in chloroplast and smaller plastoglobule in chromoplast compared to those of AC08-045. Virus-induced gene silencing (VIGS) of PRR2 significantly reduced carotenoid accumulation in Capsicum annuum 'Micropep Yellow', which carries non-functional mutations in both PSY1 and CCS. Furthermore, sequence analysis of PSY1, CCS, and PRR2 in other white pepper accessions of C. annuum and Capsicum chinense showed that they commonly have non-functional alleles in PSY1, CCS, and PRR2. Thus, our data demonstrate that the fruit color locus C1 in Capsicum spp. corresponds to the gene PRR2.

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.

Genome Editing of eIF4E1 in Tomato Confers Resistance to Pepper Mottle Virus

Many of the recessive virus-resistance genes in plants encode eukaryotic translation initiation factors (eIFs), including eIF4E, eIF4G, and related proteins. Notably, eIF4E and its isoform eIF(iso)4E are pivotal for viral infection and act as recessive resistance genes against various potyviruses in a wide range of plants. In this study, we used Clustered Regularly Interspaced Palindromic Repeats/CRISPR-associated protein 9 (CRISPR/Cas9)-mediated targeted mutagenesis to test whether novel sequence-specific mutations at eIF4E1 in Solanum lycopersicum (tomato) cv. Micro-Tom could confer enhanced resistance to potyviruses. This approach produced heritable homozygous mutations in the transgene-free E₁ generation. Sequence analysis of eIF4E1 from E₀ transgenic plants expressing Cas9 and eIF4E-sgRNA transcripts identified chimeric deletions ranging from 11 to 43 bp. Genotype analysis of the eIF4E1-edited lines in E₀, E₁, and E₂ transgenic tomato plants showed that the mutations were transmitted to subsequent generations. When homozygous mutant lines were tested for resistance to potyviruses, they exhibited no resistance to tobacco etch virus (TEV). Notably, however, several mutant lines showed no accumulation of viral particles upon infection with pepper mottle virus (PepMoV). These results indicate that site-specific mutation of tomato eIF4E1 successfully conferred enhanced resistance to PepMoV. Thus, this study demonstrates the feasibility of the use of CRISPR/Cas9 approach to accelerate breeding for trait improvement in tomato plants.

Effects of the total dose and duration of γ-irradiation on the growth responses and induced SNPs of a Cymbidium hybrid

Purpose: Ionizing radiation has been used for developing new cultivars of diverse plant species, including Cymbidium orchid species. The effects of the total dose on mutation induction have been investigated; however, there is relatively little research on the influence of the dose rate or irradiation duration.Materials and methods: Thus, we analyzed the effects of the total dose and irradiation duration on the growth of Cymbidium hybrid RB001 protocorm-like bodies (PLBs). We completed a genotyping-by-sequencing analysis to compare the induced SNPs among five γ-irradiated populations with similar growth responses (LD₅₀) to γ-rays.Results: The optimal time to assess the effects of the γ-irradiation was at 6 months after the treatment. On the basis of the survival rate of γ-irradiated PLBs, the optimal doses (LD₅₀) for each irradiation duration were estimated: 1 h, 16.1 Gy; 4 h, 23.6 Gy; 8 h, 37.9 Gy; 16 h, 37.9 Gy; and 24 h, 40.0 Gy. The estimated optimal doses were duration-dependent at irradiation durations shorter than 8 h, but not at irradiation durations exceeding 8 h. A SNP comparison revealed a lack of significant differences among the mutations induced by γ-irradiations.Conclusions: These results indicate the irradiation duration affects PLB growth in response to γ-rays. Moreover, the mutations induced by a short-term treatment may be similar to those induced by a treatment over a longer period.

Mitotypes Based on Structural Variation of Mitochondrial Genomes Imply Relationships With Morphological Phenotypes and Cytoplasmic Male Sterility in Peppers

Plant mitochondrial genomes characteristically contain extensive structural variation that can be used to define and classify cytoplasm types. We developed markers based on structural variation in the mitochondrial genomes of fertile and cytoplasmic male sterility (CMS) pepper lines and applied them to a panel of Capsicum accessions. We designed a total of 20 sequence characterized amplified region (SCAR) markers based on DNA rearrangement junctions or cytoplasm-specific segments that did not show high similarity to any nuclear mitochondrial DNA segments. We used those markers to classify the mitotypes of 96 C. annuum accessions into 15 groups. Precise genotyping of other Capsicum species (C. frutescens, C. chinense, and C. baccatum) was hampered because of various stoichiometric levels of marker amplicons. We developed a multiplex PCR system based on four of the markers that efficiently classified the C. annuum accessions into five mitotype groups. Close relationships between specific mitotypes and morphological phenotypes implied that diversification or domestication of C. annuum germplasm might have been accompanied by structural rearrangements of mitochondrial DNA or the selection of germplasms with specific mitotypes. Meanwhile, CMS lines shared the same amplification profile of markers with another mitotype. Further analysis using mitochondrial DNA (mtDNA) markers based on single-nucleotide polymorphisms (SNPs) or insertions and deletions (InDels) and CMS-specific open reading frames (orfs) provided new information about the origin of the CMS-specific mitotype and evaluation of candidates for CMS-associated genes, respectively.

A non-LTR retrotransposon activates anthocyanin biosynthesis by regulating a MYB transcription factor in Capsicum annuum

The flavonoid compound anthocyanin is an important plant metabolite with nutritional and aesthetic value as well as anti-oxidative capacity. MYB transcription factors are key regulators of anthocyanin biosynthesis in plants. In pepper (Capsicum annuum), the CaAn2 gene, encoding an R2R3 MYB transcription factor, regulates anthocyanin biosynthesis. However, no functional study or structural analysis of functional and dysfunctional CaAn2 alleles has been performed. Here, to elucidate the function of CaAn2, we generated transgenic Nicotiana benthamiana and Arabidopsis thaliana plants expressing CaAn2. All of the tissues in these plants were purple. Promoter analysis of CaAn2 in purple C. annuum 'KC00134' plants revealed the insertion of a non-long terminal repeat (LTR) retrotransposon designated Ca-nLTR-A. To determine the promoter activity and functional domain of Ca-nLTR-A, various constructs carrying different domains of Ca-nLTR-A fused with GUS were transformed into N. benthamiana. Promoter analysis showed that the 3' untranslated region (UTR) of the second open reading frame of Ca-nLTR-A is responsible for CaAn2 expression in 'KC00134'. Sequence analysis of Ca-nLTR-A identified transcription factor binding sites known to regulate anthocyanin biosynthesis. This study indicates that insertion of a non-LTR retrotransposon in the promoter may activate expression of CaAn2 by recruiting transcription factors at the 3' UTR and thus provides the first example of exaptation of a non-LTR retrotransposon into a new promoter in plants.

Current views on temperature-modulated R gene-mediated plant defense responses and tradeoffs between plant growth and immunity

Elevated ambient temperatures will likely be a key consequence of climate change over the next few decades. Adverse climatic changes could make crop plants more vulnerable to a number of biotic and abiotic stresses, which would have a major impact on worldwide food production in the future. Recent studies have indicated that elevated temperatures directly and/or indirectly affect plant-pathogen interactions. Elevated temperatures alter multiple signal transduction pathways related to stress responses in the host plant. High temperatures can also influence plant pathogenesis, but little is known about the molecular mechanisms associated with such effects. An improved understanding of the molecular genetic mechanisms involved in plant immune responses under elevated temperatures will be essential to mitigate the adverse effects of climate change to ensure future food security. In this review, we discuss recent advances in our understanding of the effects of temperature on resistance (R) gene and/or regulators of R genes in plant defense responses and summarize current evidence for tradeoffs between plant growth and immunity.

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.

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.

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.

Tandem gene duplication and recombination at the AT3 locus in the Solanaceae, a gene essential for capsaicinoid biosynthesis in Capsicum

Capsaicinoids are compounds synthesized exclusively in the genus Capsicum and are responsible for the burning sensation experienced when consuming hot pepper fruits. To date, only one gene, AT3, a member of the BAHD family of acyltransferases, is currently known to have a measurable quantitative effect on capsaicinoid biosynthesis. Multiple AT3 paralogs exist in the Capsicum genome, but their evolutionary relationships have not been characterized well. Recessive alleles at this locus result in absence of capsaicinoids in pepper fruit. To explore the evolution of AT3 in Capsicum and the Solanaceae, we sequenced this gene from diverse Capsicum genotypes and species, along with a number of representative solanaceous taxa. Our results revealed that the coding region of AT3 is highly conserved throughout the family. Further, we uncovered a tandem duplication that predates the diversification of the Solanaceae taxa sampled in this study. This pair of tandem duplications were designated AT3-1 and AT3-2. Sequence alignments showed that the AT3-2 locus, a pseudogene, retains regions of amino acid conservation relative to AT3-1. Gene tree estimation demonstrated that AT3-1 and AT3-2 form well supported, distinct clades. In C. rhomboideum, a non-pungent basal Capsicum species, we describe a recombination event between AT3-1 and AT3-2 that modified the putative active site of AT3-1, also resulting in a frame-shift mutation in the second exon. Our data suggest that duplication of the original AT3 representative, in combination with divergence and pseudogene degeneration, may account for the patterns of sequence divergence and punctuated amino acid conservation observed in this study. Further, an early rearrangement in C. rhomboidium could account for the absence of pungency in this Capsicum species.

Physical Localization of the Root-Knot Nematode (Meloidogyne incognita) Resistance Locus Me7 in Pepper (Capsicum annuum)

The root-knot nematode (RKN) Meloidogyne incognita severely reduces yields of pepper (Capsicum annuum) worldwide. A single dominant locus, Me7, conferring RKN resistance was previously mapped on the long arm of pepper chromosome P9. In the present study, the Me7 locus was fine mapped using an F₂ population of 714 plants derived from a cross between the RKN-susceptible parent C. annuum ECW30R and the RKN-resistant parent C. annuum CM334. CM334 exhibits suppressed RKN juvenile movement, suppressed feeding site enlargement and significant reduction in gall formation compared with ECW30R. RKN resistance screening in the F₂ population identified 558 resistant and 156 susceptible plants, which fit a 3:1 ratio confirming that this RKN resistance was controlled by a single dominant gene. Using the C. annuum CM334 reference genome and BAC library sequencing, fine mapping of Me7 markers was performed. The Me7 locus was delimited between two markers G21U3 and G43U3 covering a physical interval of approximately 394.7 kb on the CM334 chromosome P9. Nine markers co-segregated with the Me7 gene. A cluster of 25 putative nucleotide-binding site and leucine-rich repeat (NBS-LRR)-type disease resistance genes were predicted in the delimited Me7 region. We propose that RKN resistance in CM334 is mediated by one or more of these NBS-LRR class R genes. The Me7-linked markers identified here will facilitate marker-assisted selection (MAS) for RKN resistance in pepper breeding programs, as well as functional analysis of Me7 candidate genes in C. annuum.

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

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

HR-Mediated Defense Response is Overcome at High Temperatures in Capsicum Species

Resistance to Tomato spotted wilt virus isolated from paprika (TSWV-Pap) was overcome at high temperatures (30 ± 2°C) in both accessions of Capsicum annuum S3669 (Hana Seed Company) and C. chinense PI15225 (AVRDC Vegetable Genetic Resources). S3669 and PI15225, which carrying the Tsw gene, were mechanically inoculated with TSWV-Pap, and then maintained in growth chambers at temperatures ranging from 15 ± 2°C to 30 ± 2°C (in 5°C increments). Seven days post inoculation (dpi), a hypersensitivity reaction (HR) was induced in inoculated leaves of PI152225 and S3669 plants maintained at 25°C ± 2°C. Meanwhile, necrotic spots were formed in upper leaves of 33% of PI15225 plants maintained at 30 ± 2°C, while systemic mottle symptoms developed in 50% of S3669 plants inoculated. By 15 dpi, 25% of S3669 plants had recovered from systemic mottling induced at 30 ± 2°C. These results demonstrated that resistance to TSWV-Pap can be overcome at higher temperatures in both C. chinense and C. annuum. This is the first study reporting the determination of temperatures at which TSWV resistance is overcome in a C. annuum genetic resource expressing the Tsw gene. Our results indicated that TSWV resistance shown from pepper plants possess the Tsw gene could be overcome at high temperature. Thus, breeders should conduct evaluation of TSWV resistance in pepper cultivars at higher temperature than 30°C (constant temperature).

Fine Mapping of the Dominant Potyvirus Resistance Gene Pvr7 Reveals a Relationship with Pvr4 in Capsicum annuum

Pepper mottle virus (PepMoV) is the most common potyvirus infection of pepper plants and causes significant yield losses. The Pvr7 gene from Capsicum chinense PI159236 and the Pvr4 gene from C. annuum CM334 both have been reported to confer dominant resistance to PepMoV. The Pvr7 locus conferring resistance to PepMoV in C. annuum '9093' was previously mapped to chromosome 10. To develop a high-resolution map of the Pvr7 locus in 9093, we constructed an intraspecific F₂ mapping population consisting of 916 individuals by crossing PepMoV-resistant C. annuum '9093' and the PepMoV-susceptible C. annuum 'Jeju'. To delimit the Pvr7 target region, single-nucleotide polymorphism (SNP) markers derived from the Pvr4 region were used for genotyping the F₂ population. Molecular mapping delimited the Pvr7 locus to a physical interval of 258 kb, which was the same region as Pvr4 on chromosome 10. Three SNP markers derived from Pvr4 mapping perfectly cosegregated with PepMoV resistance. Sequencing analyses of the Pvr7 flanking markers and the Pvr4-specific gene indicated that Pvr7 and Pvr4 are the same gene. Resistance spectrum analysis of 9093 against pepper potyviruses showed that 9093 has a resistance spectrum similar to that of cultivar CM334. These combined results demonstrate that, unlike previously thought, the dominant PepMoV resistance in 9093 could be derived from C. annuum 'CM334', and that Pvr4 and Pvr7 should be considered as the same locus.

Large Sternocleidomastoid Haematoma after Minor Trauma in a Postradiation Neck

Sternocleidomastoid (SCM) muscle haematoma after minor blunt injury is rare. Radiation therapy as a treatment for head and neck tumours e.g. nasopharyngeal cancer can result in complications, including fibrosis, scarring, and atrophy. We report a 51-year-old man suffering from massive haematoma of SCM muscle after minor blunt injury to a postradiation neck. Signs of haematoma were absent during his first visit to the emergency department. Emergency physicians need to be aware that bleeding from a postradiation neck, even after minor trauma, may be life threatening.

New reference genome sequences of hot pepper reveal the massive evolution of plant disease-resistance genes by retroduplication

Background

Transposable elements are major evolutionary forces which can cause new genome structure and species diversification. The role of transposable elements in the expansion of nucleotide-binding and leucine-rich-repeat proteins (NLRs), the major disease-resistance gene families, has been unexplored in plants.

Results

We report two high-quality de novo genomes (Capsicum baccatum and C. chinense) and an improved reference genome (C. annuum) for peppers. Dynamic genome rearrangements involving translocations among chromosomes 3, 5, and 9 were detected in comparison between C. baccatum and the two other peppers. The amplification of athila LTR-retrotransposons, members of the gypsy superfamily, led to genome expansion in C. baccatum. In-depth genome-wide comparison of genes and repeats unveiled that the copy numbers of NLRs were greatly increased by LTR-retrotransposon-mediated retroduplication. Moreover, retroduplicated NLRs are abundant across the angiosperms and, in most cases, are lineage-specific.

Conclusions

Our study reveals that retroduplication has played key roles for the massive emergence of NLR genes including functional disease-resistance genes in pepper plants.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-017-1341-9) contains supplementary material, which is available to authorized users.

New reference genome sequences of hot pepper reveal the massive evolution of plant disease-resistance genes by retroduplication

BACKGROUND

Transposable elements are major evolutionary forces which can cause new genome structure and species diversification. The role of transposable elements in the expansion of nucleotide-binding and leucine-rich-repeat proteins (NLRs), the major disease-resistance gene families, has been unexplored in plants.

RESULTS

We report two high-quality de novo genomes (Capsicum baccatum and C. chinense) and an improved reference genome (C. annuum) for peppers. Dynamic genome rearrangements involving translocations among chromosomes 3, 5, and 9 were detected in comparison between C. baccatum and the two other peppers. The amplification of athila LTR-retrotransposons, members of the gypsy superfamily, led to genome expansion in C. baccatum. In-depth genome-wide comparison of genes and repeats unveiled that the copy numbers of NLRs were greatly increased by LTR-retrotransposon-mediated retroduplication. Moreover, retroduplicated NLRs are abundant across the angiosperms and, in most cases, are lineage-specific.

CONCLUSIONS

Our study reveals that retroduplication has played key roles for the massive emergence of NLR genes including functional disease-resistance genes in pepper plants.

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

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

Molecular Mapping of PMR1, a Novel Locus Conferring Resistance to Powdery Mildew in Pepper (Capsicum annuum)

Powdery mildew, caused by Leveillula taurica, is a major fungal disease affecting greenhouse-grown pepper (Capsicum annuum). Powdery mildew resistance has a complex mode of inheritance. In the present study, we investigated a novel powdery mildew resistance locus, PMR1, using two mapping populations: 102 'VK515' F_2:3 families (derived from a cross between resistant parental line 'VK515R' and susceptible parental line 'VK515S') and 80 'PM Singang' F₂ plants (derived from the F₁ 'PM Singang' commercial hybrid). Genetic analysis of the F_2:3 'VK515' and F₂ 'PM Singang' populations revealed a single dominant locus for inheritance of the powdery mildew resistance trait. Genetic mapping showed that the PMR1 locus is located on syntenic regions of pepper chromosome 4 in a 4-Mb region between markers CZ2_11628 and HRM4.1.6 in 'VK515R'. Six molecular markers including one SCAR marker and five SNP markers were localized to a region 0 cM from the PMR1 locus. Two putative nucleotide-binding site leucine-rich repeat (NBS-LRR)-type disease resistance genes were identified in this PMR1 region. Genotyping-by-sequencing (GBS) and genetic mapping analysis revealed suppressed recombination in the PMR1 region, perhaps due to alien introgression. In addition, a comparison of species-specific InDel markers as well as GBS-derived SNP markers indicated that C. baccatum represents a possible source of such alien introgression of powdery mildew resistance into 'VK515R'. The molecular markers developed in this study will be especially helpful for marker-assisted selection in pepper breeding programs for powdery mildew resistance.