Whole-genome sequencing of cultivated and wild peppers provides insights into Capsicum domestication and specialization

As an economic crop, pepper satisfies people's spicy taste and has medicinal uses worldwide. To gain a better understanding of Capsicum evolution, domestication, and specialization, we present here the genome sequence of the cultivated pepper Zunla-1 (C. annuum L.) and its wild progenitor Chiltepin (C. annuum var. glabriusculum). We estimate that the pepper genome expanded similar to 0.3 Mya (with respect to the genome of other Solanaceae) by a rapid amplification of retrotransposons elements, resulting in a genome comprised of similar to 81% repetitive sequences. Approximately 79% of 3.48-Gb scaffolds containing 34,476 protein-coding genes were anchored to chromosomes by a high-density genetic map. Comparison of cultivated and wild pepper genomes with 20 resequencing accessions revealed molecular footprints of artificial selection, providing us with a list of candidate domestication genes. We also found that dosage compensation effect of tandem duplication genes probably contributed to the pungent diversification in pepper. The Capsicum reference genome provides crucial information for the study of not only the evolution of the pepper genome but also, the Solanaceae family, and it will facilitate the establishment of more effective pepper breeding programs

Differential Proteomic Analysis of Anthers between Cytoplasmic Male Sterile and Maintainer Lines in Capsicum annuum L.

Cytoplasmic male sterility (CMS), widely used in the production of hybrid seeds, is a maternally inherited trait resulting in a failure to produce functional pollen. In order to identify some specific proteins associated with CMS in pepper, two-dimensional gel electrophoresis (2-DE) was applied to proteomic analysis of anthers/buds between a CMS line (designated NA3) and its maintainer (designated NB3) in Capsicum annuum L. Thirty-three spots showed more than 1.5-fold in either CMS or its maintainer. Based on mass spectrometry, 27 spots representing 23 distinct proteins in these 33 spots were identified. Proteins down-regulated in CMS anthers/buds includes ATP synthase D chain, formate dehydrogenase, alpha-mannosidas, RuBisCO large subunit-binding protein subunit beta, chloroplast manganese stabilizing protein-II, glutathione S-transferase, adenosine kinase isoform 1T-like protein, putative DNA repair protein RAD23-4, putative caffeoyl-CoA 3-O-methyltransferase, glutamine synthetase (GS), annexin Cap32, glutelin, allene oxide cyclase, etc. In CMS anthers/buds, polyphenol oxidase, ATP synthase subunit beta, and actin are up-regulated. It was predicted that male sterility in NA3 might be related to energy metabolism turbulence, excessive ethylene synthesis, and suffocation of starch synthesis. The present study lays a foundation for future investigations of gene functions associated with pollen development and cytoplasmic male sterility, and explores the molecular mechanism of CMS in pepper

A RAD-Based Genetic Map for Anchoring Scaffold Sequences and Identifying QTLs in Bitter Gourd (Momordica charantia)

Genetic mapping is a basic tool necessary for anchoring assembled scaffold sequences and for identifying QTLs controlling important traits. Though bitter gourd (Momordica charantia) is both consumed and used as a medicinal, research on its genomics and genetic mapping is severely limited. Here, we report the construction of a restriction site associated DNA (RAD)-based genetic map for bitter gourd using an F2 mapping population comprising 423 individuals derived from two cultivated inbred lines, the gynoecious line ‘K44’ and the monoecious line ‘Dali-11.’ This map comprised 1,009 SNP markers and spanned a total genetic distance of 2,203.95 cM across the 11 linkage groups. It anchored a total of 113 assembled scaffolds that covered about 251.32 Mb (85.48%) of the 294.01 Mb assembled genome. In addition, three horticulturally important traits including sex expression, fruit epidermal structure, and immature fruit color were evaluated using a combination of qualitative and quantitative data. As a result, we identified three QTL/gene loci responsible for these traits in three environments. The QTL/gene gy/fffn/ffn, controlling sex expression involved in gynoecy, first female flower node, and female flower number was detected in the reported region. Particularly, two QTLs/genes, Fwa/Wr and w, were found to be responsible for fruit epidermal structure and white immature fruit color, respectively. This RAD-based genetic map promotes the assembly of the bitter gourd genome and the identified genetic loci will accelerate the cloning of relevant genes in the future

Genome-wide identification and analysis of highly specific CRISPR/Cas9 editing sites in pepper (Capsicum annuum L.).

The CRISPR/Cas9 system is an efficient genome editing tool that possesses the outstanding advantages of simplicity and high efficiency. Genome-wide identification and specificity analysis of editing sites is an effective approach for mitigating the risk of off-target effects of CRISPR/Cas9 and has been applied in several plant species but has not yet been reported in pepper. In present study, we first identified genome-wide CRISPR/Cas9 editing sites based on the 'Zunla-1' reference genome and then evaluated the specificity of CRISPR/Cas9 editing sites through whole-genome alignment. Results showed that a total of 603,202,314 CRISPR/Cas9 editing sites, including 229,909,837 (~38.11%) NGG-PAM sites and 373,292,477 (~61.89%) NAG-PAM sites, were detectable in the pepper genome, and the systematic characterization of their composition and distribution was performed. Furthermore, 29,623,855 highly specific NGG-PAM sites were identified through whole-genome alignment analysis. There were 26,699,38 (~90.13%) highly specific NGG-PAM sites located in intergenic regions, which was 9.13 times of the number in genic regions, but the average density in genic regions was higher than that in intergenic regions. More importantly, 34,251 (~96.93%) out of 35,336 annotated genes exhibited at least one highly specific NGG-PAM site in their exons, and 90.50% of the annotated genes exhibited at least 4 highly specific NGG- PAM sites, indicating that the set of highly specific CRISPR/Cas9 editing sites identified in this study was widely applicable and conducive to the minimization of the off-target effects of CRISPR/Cas9 in pepper

Genome-Wide Analysis of Simple Sequence Repeats in Bitter Gourd (Momordica charantia)

Bitter gourd (Momordica charantia) is widely cultivated as a vegetable and medicinal herb in many Asian and African countries. After the sequencing of the cucumber (Cucumis sativus), watermelon (Citrullus lanatus), and melon (Cucumis melo) genomes, bitter gourd became the fourth cucurbit species whose whole genome was sequenced. However, a comprehensive analysis of simple sequence repeats (SSRs) in bitter gourd, including a comparison with the three aforementioned cucurbit species has not yet been published. Here, we identified a total of 188,091 and 167,160 SSR motifs in the genomes of the bitter gourd lines ‘Dali-11’ and ‘OHB3-1,’ respectively. Subsequently, the SSR content, motif lengths, and classified motif types were characterized for the bitter gourd genomes and compared among all the cucurbit genomes. Lastly, a large set of 138,727 unique in silico SSR primer pairs were designed for bitter gourd. Among these, 71 primers were selected, all of which successfully amplified SSRs from the two bitter gourd lines ‘Dali-11’ and ‘K44’. To further examine the utilization of unique SSR primers, 21 SSR markers were used to genotype a collection of 211 bitter gourd lines from all over the world. A model-based clustering method and phylogenetic analysis indicated a clear separation among the geographic groups. The genomic SSR markers developed in this study have considerable potential value in advancing bitter gourd research

Table_1_A point mutation in MC06g1112 encoding FLOWERING LOCUS T decreases the first flower node in bitter gourd (Momordica charantia L.).xlsx

In Cucurbitaceae crops, the first flower node (FFN) is an important agronomic trait which can impact the onset of maturity, the production of female flowers, and yield. However, the gene responsible for regulating FFN in bitter gourd is unknown. Here, we used a gynoecious line (S156G) with low FFN as the female parent and a monoecious line (K8-201) with high FFN as the male parent to obtain F1 and F2 generations. Genetic analysis indicated that the low FFN trait was incompletely dominant over the high FFN trait. A major quantitative trait locus (QTL)-Mcffn and four minor effect QTLs-Mcffn1.1, Mcffn1.2, Mcffn1.3, and Mcffn1.4 were detected by whole-genome re-sequencing-based QTL mapping in the S156G×K8-201 F2 population (n=234) cultivated in autumn 2019. The Mcffn locus was further supported by molecular marker-based QTL mapping in three S156G×K8-201 F2 populations planted in autumn 2019 (n=234), autumn 2020 (n=192), and spring 2022 (n=205). Then, the Mcffn locus was fine-mapped into a 77.98-kb physical region on pseudochromosome MC06 using a large S156G×K8-201 F2 population (n=2,402). MC06g1112, which is a homolog of FLOWERING LOCUS T (FT), was considered as the most likely Mcffn candidate gene according to both expression and sequence variation analyses between parental lines. A point mutation (C277T) in MC06g1112, which results in a P93S amino acid mutation between parental lines, may be responsible for decreasing FFN in bitter gourd. Our findings provide a helpful resource for the molecular marker-assisted selective breeding of bitter gourd.</p

DataSheet_1_A point mutation in MC06g1112 encoding FLOWERING LOCUS T decreases the first flower node in bitter gourd (Momordica charantia L.).docx

In Cucurbitaceae crops, the first flower node (FFN) is an important agronomic trait which can impact the onset of maturity, the production of female flowers, and yield. However, the gene responsible for regulating FFN in bitter gourd is unknown. Here, we used a gynoecious line (S156G) with low FFN as the female parent and a monoecious line (K8-201) with high FFN as the male parent to obtain F1 and F2 generations. Genetic analysis indicated that the low FFN trait was incompletely dominant over the high FFN trait. A major quantitative trait locus (QTL)-Mcffn and four minor effect QTLs-Mcffn1.1, Mcffn1.2, Mcffn1.3, and Mcffn1.4 were detected by whole-genome re-sequencing-based QTL mapping in the S156G×K8-201 F2 population (n=234) cultivated in autumn 2019. The Mcffn locus was further supported by molecular marker-based QTL mapping in three S156G×K8-201 F2 populations planted in autumn 2019 (n=234), autumn 2020 (n=192), and spring 2022 (n=205). Then, the Mcffn locus was fine-mapped into a 77.98-kb physical region on pseudochromosome MC06 using a large S156G×K8-201 F2 population (n=2,402). MC06g1112, which is a homolog of FLOWERING LOCUS T (FT), was considered as the most likely Mcffn candidate gene according to both expression and sequence variation analyses between parental lines. A point mutation (C277T) in MC06g1112, which results in a P93S amino acid mutation between parental lines, may be responsible for decreasing FFN in bitter gourd. Our findings provide a helpful resource for the molecular marker-assisted selective breeding of bitter gourd.</p

Image_2.tif

<p>Genetic mapping is a basic tool necessary for anchoring assembled scaffold sequences and for identifying QTLs controlling important traits. Though bitter gourd (Momordica charantia) is both consumed and used as a medicinal, research on its genomics and genetic mapping is severely limited. Here, we report the construction of a restriction site associated DNA (RAD)-based genetic map for bitter gourd using an F₂ mapping population comprising 423 individuals derived from two cultivated inbred lines, the gynoecious line ‘K44’ and the monoecious line ‘Dali-11.’ This map comprised 1,009 SNP markers and spanned a total genetic distance of 2,203.95 cM across the 11 linkage groups. It anchored a total of 113 assembled scaffolds that covered about 251.32 Mb (85.48%) of the 294.01 Mb assembled genome. In addition, three horticulturally important traits including sex expression, fruit epidermal structure, and immature fruit color were evaluated using a combination of qualitative and quantitative data. As a result, we identified three QTL/gene loci responsible for these traits in three environments. The QTL/gene gy/fffn/ffn, controlling sex expression involved in gynoecy, first female flower node, and female flower number was detected in the reported region. Particularly, two QTLs/genes, Fwa/Wr and w, were found to be responsible for fruit epidermal structure and white immature fruit color, respectively. This RAD-based genetic map promotes the assembly of the bitter gourd genome and the identified genetic loci will accelerate the cloning of relevant genes in the future.</p