A major quantitative trait locus conferring subgynoecy in cucumber

Fine mapping and identification of CqMYB62 as the subgynoecy gene in chieh-qua (Benincasa hispida Cogn. var. Chieh-qua How)

KEY MESSAGE

It was hypothesized that Bch08G003160 (CqMYB62), located in the 51.08 Kb region on chromosome 08, might be an important candidate gene for the subgynoecy trait in chieh-qua, based on BSA-seq and linkage mapping approaches. In cucurbit crops, the use of female lines can greatly increase the yield of a single plant and is especially important for the production of hybrid generation seeds, thus being of great interest to breeders. To identify genes regulating sex differentiation in chieh-qua, genetic analysis of the subgynoecy trait was conducted using a chieh-qua F2 population. Initial localization of the locus was done using BSA-seq, followed by fine mapping with a large F2 population (n = 2,741). The locus was ultimately narrowed down to a 51.08 Kb region on Chr08, revealing a single gene Bch08G003160 (CqMYB62) in this region. Further analysis revealed that the presence of two variant loci (SNP_416 and SNP_317) in the coding region resulted in premature termination of the codon and loss of function of the structural domain of the PLN03212 superfamily. Moreover, our research indicated that the subgynoecy trait mediated by CqMYB62 in chieh-qua is potentially regulated by gibberellic acid (GA). Two efficient dCAPS markers were developed to distinguish subgynoecy. In summary, these findings highlight the critical role of CqMYB62 in subgynoecy trait regulation, offering potential implications for chieh-qua breeding programs.

The calcium-dependent protein kinase CmaCPK4 regulates sex determination in pumpkin (Cucurbita maxima D.)

Pumpkin (Cucurbita maxima D.) is typically monoecious with individual male and female flowers, and its yield is associated with the degree of femaleness, i.e. the ratio of female to male flowers produced by the plant. Subgynoecy represents a sex form with a high degree of femaleness, but the regulatory mechanisms in pumpkin remain poorly understood. In this study, using the F2 population crossed from the subgynoecious line 2013-12 and the monoecious line 9-6, we initially identified a recessive locus to control the subgynoecious trait and named it sg1. After bulked segregant analysis with whole-genome resequencing and molecular marker linkage analysis, the sg1 locus was mapped to pumpkin Chromosome 2. Genetic sequence analysis found a pumpkin calcium-dependent protein kinase (CPK) gene, CmaCPK4, in the mapping interval as the candidate gene. A retrotransposon insertion identified within the promoter elevated CmaCPK4 expression in 2013-12. Morphological characterization of near-isogenic lines containing the sg1 allele showed increases in the ratio of female flowers and high ethylene contents in terminal buds compared with the receptor parent. Heterologous overexpression of CmaCPK4 significantly increased the ratio of female flowers in cucumber (Cucumis sativus). Furthermore, CmaCPK4 directly interacts with and phosphorylates 1-aminocyclopropane-1-carboxylate synthase 5 (CmaACS5) and 1-aminocyclopropane-1-carboxylate synthase 7 (CmaACS7), resulting in increased ethylene content in 2013-12, which affected pumpkin sex determination. These findings provide insights into the role of the CmaCPK4-CmaACS5/CmaACS7 module in ethylene-induced sex determination in pumpkin.

Chromosome-level genome assembly and population genomics reveals crucial selection for subgynoecy development in chieh-qua

Chieh-qua is an important cucurbit crop and very popular in South China and Southeast Asia. Despite its significance, its genetic basis and domestication history are unclear. In this study, we have successfully generated a chromosome-level reference genome assembly for the chieh-qua 'A36' using a hybrid assembly strategy that combines PacBio long reads and Illumina short reads. The assembled genome of chieh-qua is approximately 953.3 Mb in size and is organized into 12 chromosomes, with contig N50 of 6.9 Mb and scaffold N50 of 68.2 Mb. Notably, the chieh-qua genome is comparable in size to the wax gourd genome. Through gene prediction analysis, we have identified a total of 24 593 protein-coding genes in the A36 genome. Additionally, approximately 56.6% (539.3 Mb) of the chieh-qua genome consists of repetitive sequences. Comparative genome analysis revealed that chieh-qua and wax gourd are closely related, indicating a close evolutionary relationship between the two species. Population genomic analysis, employing 129 chieh-qua accessions and 146 wax gourd accessions, demonstrated that chieh-qua exhibits greater genetic diversity compared to wax gourd. We also employed the GWAS method to identify related QTLs associated with subgynoecy, an interested and important trait in chieh-qua. The MYB59 (BhiCQ0880026447) exhibited relatively high expression levels in the shoot apex of four subgynoecious varieties compared with monoecious varieties. Overall, this research provides insights into the domestication history of chieh-qua and offers valuable genomic resources for further molecular research.

Genetic mapping and quantitative trait loci analysis for pistillate flowers per node and multi-pistillate flower traits in the F2 cucumber population

This study focused on cucumbers' multi-pistillate flower (MPF) trait, which is essential for high yields. A genetic linkage map was constructed using a population of 219 F2 plants to analyze quantitative trait loci (QTL) associated with MPF traits. Crossbreeding of EWSCU-809 (MPF) with EWSCU-989 (single pistillate flower: SPF) generated an F1 hybrid that self-pollinated to form an F2 population. Based on 244 single nucleotide polymorphic markers across seven cucumber chromosomes, a linkage map facilitated QTL analysis considering average pistillate flowers (PFs) per node and nodes with MPF traits. The results indicated a 9:6:1 epistatic ratio in the F2 populations, revealing recessive allele control of the MPF trait in gynoecious plants. Three QTLs (qMP2.1, qMP3, qMP7) on chromosomes 2, 3, and 7 were associated with average PFs per node, explaining 5.6 to 10.3% of phenotypic variation. Four QTLs (qMP2.2, qMP3, qMP4, qMP7) on chromosomes 2, 3, 4, and 7 were linked to the presence of nodes with MPF traits, explaining 5.8 to 10.6% of phenotypic variance. Notably, QTL regions overlapped between the two datasets, suggesting pleiotropic effects, particularly on chromosomes 3 and 7. These reliable QTLs have the potential to improve breeding programs, enhance PF development, and increase cucumber yields.

Identification of novel sex determination loci in Japanese weedy melon

KEY MESSAGE

Japanese weedy melon exhibits unique sex expression with interactions between previously reported sex determination genes and two novel loci. Sex expression contributes to fruit quality and yield in the Cucurbitaceae. In melon, orchestrated regulation by sex determination genes explains the mechanism of sex expression, resulting in a great variety of sexual morphologies. In this study, we examined the Japanese weedy melon UT1, which does not follow the reported model of sex expression. We conducted QTL analysis using F₂ plants for flower sex on the main stem and the lateral branch and mapped "occurrence of pistil-bearing flower on the main stem" locus on Chr. 3 (Opbf3.1) and "type of pistil-bearing flower" (female or bisexual) loci on Chr. 2 (tpbf2.1) and Chr. 8 (tpbf8.1). The Opbf3.1 included the known sex determination gene CmACS11. Sequence comparison of CmACS11 between parental lines revealed three nonsynonymous SNPs. A CAPS marker developed from one of the SNPs was closely linked to the occurrence of pistil-bearing flowers on the main stem in two F₂ populations with different genetic backgrounds. The UT1 allele on Opbf3.1 was dominant in F₁ lines from crosses between UT1 and diverse cultivars and breeding lines. This study suggests that Opbf3.1 and tpbf8.1 may promote the development of pistil and stamen primordia by inhibiting CmWIP1 and CmACS-7 functions, respectively, making the UT1 plants hermaphrodite. The results of this study provide new insights into the molecular mechanisms of sex determination in melons and considerations for the application of femaleness in melon breeding.

Recent progress on the molecular breeding of Cucumis sativus L. in China

Molecular breeding of Cucumis sativus L. is based on traditional breeding techniques and modern biological breeding in China. There are opportunities for further breeding improvement by molecular design breeding and the automation of phenotyping technology using untapped sources of genetic diversity. Cucumber (Cucumis sativus L.) is an important vegetable cultivated worldwide. It bears fruits of light fragrance, and crisp texture with high nutrition. China is the largest producer and consumer of cucumber, accounting for 70% of the world's total production. With increasing consumption demand, the production of Cucurbitaceae crops has been increasing yearly. Thus, new cultivars that can produce high-quality cucumber with high yield and easy cultivation are in need. Conventional genetic breeding has played an essential role in cucumber cultivar innovation over the past decades. However, its progress is slow due to the long breeding period, and difficulty in selecting stable genetic characters or genotypes, prompting researchers to apply molecular biotechnologies in cucumber breeding. Here, we first summarize the achievements of conventional cucumber breeding such as crossing and mutagenesis, and then focus on the current status of molecular breeding of cucumber in China, including the progress and achievements on cucumber genomics, molecular mechanism underlying important agronomic traits, and also on the creation of high-quality multi-resistant germplasm resources, new variety breeding and ecological breeding. Future development trends and prospects of cucumber molecular breeding in China are also discussed.

Gynoecy instability in cucumber (Cucumis sativus L.) is due to unequal crossover at the copy number variation-dependent Femaleness (F) locus

Cucumber, Cucumis sativus is an important vegetable crop, and gynoecy has played a critical role in yield increase of hybrid cucumber production. Cucumber has a unique genetic system for gynoecious sex expression, which is determined by the copy number variation (CNV)-based, dominant, and dosage-dependent femaleness (F) locus. However, this gynoecy expression system seems unstable since monecious plants could often be found in F-dependent gynoecious cucumber inbreds. We hypothesized that gynoecy instability (gynoecy loss) may be due to unequal crossing over (UCO) during meiosis among repeat units of the CNV. In this study, using high throughput genome resequencing, fiber-FISH and genomic qPCR analyses, we first confirmed and refined the structure of the F locus, which was a CNV of a 30.2-kb tandem repeat. Gynoecious plants contained three genes: CsACS1, CsACS1G, and CsMYB, of which CsACS1G is a duplication of CsACS1 but with a recombinant distal promoter that may contribute to gynoecy sex expression. In two large populations from self-pollinated gynoecious inbred lines, 'gynoecy loss' mutants were identified with similar mutation rates (~0.12%). We show that these monecious mutants have lost CsACS1G. In addition, we identified gynoecious lines in natural populations that carry two copies of CSACS1G. We proposed a model to explain gynoecy instability in F-dependent cucumbers, which is caused by UCO among CSACS1/G units during meiosis. The findings present a convincing case that the phenotypic variation of an economically important trait is associated with the dynamic changes of copy numbers at the F locus. This work also has important implications in cucumber breeding.

Molecularly tagged genes and quantitative trait loci in cucumber with recommendations for QTL nomenclature

Cucumber, Cucumis sativus L. (2n = 2x = 14), is an important vegetable crop worldwide. It was the first specialty crop with a publicly available draft genome. Its relatively small, diploid genome, short life cycle, and self-compatible mating system offers advantages for genetic studies. In recent years, significant progress has been made in molecular mapping, and identification of genes and QTL responsible for key phenotypic traits, but a systematic review of the work is lacking. Here, we conducted an extensive literature review on mutants, genes and QTL that have been molecularly mapped or characterized in cucumber. We documented 81 simply inherited trait genes or major-effect QTL that have been cloned or fine mapped. For each gene, detailed information was compiled including chromosome locations, allelic variants and associated polymorphisms, predicted functions, and diagnostic markers that could be used for marker-assisted selection in cucumber breeding. We also documented 322 QTL for 42 quantitative traits, including 109 for disease resistances against seven pathogens. By alignment of these QTL on the latest version of cucumber draft genomes, consensus QTL across multiple studies were inferred, which provided insights into heritable correlations among different traits. Through collaborative efforts among public and private cucumber researchers, we identified 130 quantitative traits and developed a set of recommendations for QTL nomenclature in cucumber. This is the first attempt to systematically summarize, analyze and inventory cucumber mutants, cloned or mapped genes and QTL, which should be a useful resource for the cucurbit research community.

Identification of quantitative trait loci governing subgynoecy in cucumber

QTL-seq analysis identified three major QTLs conferring subgynoecy in cucumbers. Furthermore, sequence and expression analyses predicted candidate genes controlling subgynoecy. The cucumber (Cucumis sativus L.) is a typical monoecious having individual male and female flowers, and sex differentiation is an important developmental process that directly affects its fruit yield. Subgynoecy represents a sex form with a high degree of femaleness and would have alternative use as gynoecy under limited resource conditions. Recently, many studies have been reported that QTL-seq, which integrates the advantages of bulked segregant analysis and high-throughput whole-genome resequencing, can be a rapid and cost-effective way of mapping QTLs. Segregation analysis in the F₂ and BC₁ populations derived from a cross between subgynoecious LOSUAS and monoecious BMB suggested the quantitative nature of subgynoecy in cucumbers. Both genome-wide SNP profiling of subgynoecious and monoecious bulks constructed from F₂ and BC₁ plants consistently identified three significant genomic regions, one on chromosome 3 (sg3.1) and another two on short and long arms of chromosome 1 (sg1.1 and sg1.2). Classical QTL analysis using the F₂ confirmed sg3.1 (R² = 42%), sg1.1 (R² = 29%) and sg1.2 (R² = 18%) as major QTLs. These results revealed the unique genetic inheritance of subgynoecious line LOSUAS through two distinct major QTLs, sg3.1 and sg1.1, which mainly increase degree of femaleness, while another QTL, sg1.2, contributes to decrease it. This study demonstrated that QTL-seq allows rapid and powerful detection of QTLs using preliminary generation mapping populations such as F₂ or BC₁ population and further that the identified QTLs could be useful for molecular breeding of cucumber lines with high yield potential.

Gene Interactions Regulating Sex Determination in Cucurbits

The family Cucurbitaceae includes many economically important crops, such as cucumber (Cucumis sativus), melon (Cucumis melo), watermelon (Citrullus lanatus), and zucchini (Cucurbita pepo), which share homologous gene pathways that control similar phenotypes. Sex determination is a research hotspot associated with yield and quality, and the genes involved are highly orthologous and conserved in cucurbits. In the field, six normal sex types have been categorized according to the distribution of female, male, or bisexual flowers in a given plant. To date, five orthologous genes involved in sex determination have been cloned, and their various combinations and expression patterns can explain all the identified sex types. In addition to genetic mechanisms, ethylene controls sex expression in this family. Two ethylene signaling components have been identified recently, which will help us to explore the ethylene signaling-mediated interactions among sex-related genes. This review discusses recent advances relating to the mechanism of sex determination in cucurbits and the prospects for research in this area.

Physiological, genomic and transcriptional diversity in responses to boron deficiency in rapeseed genotypes

Allotetraploid rapeseed (Brassica napus L. A_nA_nC_nC_n, 2n=4x=38) is highly susceptible to boron (B) deficiency, a widespread limiting factor that causes severe losses in seed yield. The genetic variation in the sensitivity to B deficiency found in rapeseed genotypes emphasizes the complex response architecture. In this research, a B-inefficient genotype, 'Westar 10' ('W10'), responded to B deficiencies during vegetative and reproductive development with an over-accumulation of reactive oxygen species, severe lipid peroxidation, evident plasmolysis, abnormal floral organogenesis, and widespread sterility compared to a B-efficient genotype, 'Qingyou 10' ('QY10'). Whole-genome re-sequencing (WGS) of 'QY10' and 'W10' revealed a total of 1 605 747 single nucleotide polymorphisms and 218 755 insertions/deletions unevenly distributed across the allotetraploid rapeseed genome (~1130Mb). Digital gene expression (DGE) profiling identified more genes related to B transporters, antioxidant enzymes, and the maintenance of cell walls and membranes with higher transcript levels in the roots of 'QY10' than in 'W10' under B deficiency. Furthermore, based on WGS and bulked segregant analysis of the doubled haploid (DH) line pools derived from 'QY10' and 'W10', two significant quantitative trait loci (QTLs) for B efficiency were characterized on chromosome C2, and DGE-assisted QTL-seq analyses then identified a nodulin 26-like intrinsic protein gene and an ATP-binding cassette (ABC) transporter gene as the corresponding candidates regulating B efficiency. This research facilitates a more comprehensive understanding of the differential physiological and transcriptional responses to B deficiency and abundant genetic diversity in rapeseed genotypes, and the DGE-assisted QTL-seq analyses provide novel insights regarding the rapid dissection of quantitative trait genes in plant species with complex genomes.