Fine-mapping and candidate gene analysis of the Mcgy1 locus responsible for gynoecy in bitter gourd (Momordica spp.).

The Mcgy1 locus responsible for gynoecy was fine-mapped into a 296.94-kb region, in which four single-nucleotide variations and six genes adjacent to them might be associate with sex differentiation in bitter gourd. Gynoecy plays an important role in high-efficiency hybrid seed production, and gynoecious plants are excellent materials for dissecting sex differentiation in Cucurbitaceae crop species, including bitter gourd. However, the gene responsible for gynoecy in bitter gourd is unknown. Here, we first identified a gynoecy locus designated Mcgy1 using the F population (n = 291) crossed from the gynoecious line S156G and the monoecious line K8-201 via bulked segregant analysis with whole-genome resequencing (BSA-seq) and molecular marker linkage analysis. Then, a large S156G × K8-201 F population (n = 5,656) was used for fine-mapping to delimit the Mcgy1 locus into a 296.94-kb physical region on pseudochromosome MC01, where included 33 annotated genes different from any homologous gynoecy genes previously reported in Cucurbitaceae species. Within this region, four underlying single-nucleotide variations (SNVs) that might cause gynoecy were identified by multiple genomic sequence variation analysis, and their six neighbouring genes were considered as potential candidate genes for Mcgy1. Of these, only MC01g1681 showed a significant differential expression at two-leaf developmental stage between S156G and its monoecious near-isogenic line S156 based on RNA sequencing (RNA-seq) and qRT-PCR analyses. In addition, transcriptome analysis revealed 21 key differentially expressed genes (DEGs) and possible regulatory pathways of the formation of gynoecy in bitter gourd. Our findings provide a new clue for researching on gynoecious plants in Cucurbitaceae species and a theoretical basis for breeding gynoecious bitter gourd lines by the use of molecular markers-assisted selection.