ABSTRACT
Bitter gourd is an important vegetable crop grown throughout the tropics mainly because of its high nutritional value. Sex expression and identification of gynoecious trait in cucurbitaceous vegetable crops has facilitated the hybrid breeding programme in a great way to improve productivity. In bitter gourd, gynoecious sex expression is poorly reported and detailed molecular pathways involve yet to be studied. The present experiment was conducted to study the inheritance, identify the genomic regions associated with gynoecious sex expression and to reveal possible candidate genes through QTL-seq. Segregation for the gynoecious and monoecious sex forms in the F2 progenies indicated single recessive gene controlling gynoecious sex expression in the genotype, PVGy-201. Gynoecious parent, PVGy-201, Monoecious parent, Pusa Do Mausami (PDM), and two contrasting bulks were constituted for deep-sequencing. A total of 10.56, 23.11, 15.07, and 19.38 Gb of clean reads from PVGy-201, PDM, gynoecious bulk and monoecious bulks were generated. Based on the ΔSNP index, 1.31 Mb regions on the chromosome 1 was identified to be associated with gynoecious sex expression in bitter gourd. In the QTL region 293,467 PVGy-201 unique variants, including SNPs and indels, were identified. In the identified QTL region, a total of 1019 homozygous variants were identified between PVGy1 and PDM genomes and 71 among them were non-synonymous variants (SNPS and INDELs), out of which 11 variants (7 INDELs, 4 SNPs) were classified as high impact variants with frame shift/stop gain effect. In total twelve genes associated with male and female gametophyte development were identified in the QTL-region. Ethylene-responsive transcription factor 12, Auxin response factor 6, Copper-transporting ATPase RAN1, CBL-interacting serine/threonine-protein kinase 23, ABC transporter C family member 2, DEAD-box ATP-dependent RNA helicase 1 isoform X2, Polygalacturonase QRT3-like isoform X2, Protein CHROMATIN REMODELING 4 were identified with possible role in gynoecious sex expression. Promoter region variation in 8 among the 12 genes indicated their role in determining gynoecious sex expression in bitter gourd genotype, DBGy-1. The findings in the study provides insight about sex expression in bitter gourd and will facilitate fine mapping and more precise identification of candidate genes through their functional validation.
ABSTRACT
Introduction: Momordica balsamina is the closest wild species that can be crossed with an important fruit vegetable crop, Momordica charantia, has immense medicinal value, and placed under II subclass of primary gene pool of bitter gourd. M. balsamina is tolerant to major biotic and abiotic stresses. Genome characterization of Momordica balsamina as a wild relative of bitter gourd will contribute to the knowledge of the gene pool available for improvement in bitter gourd. There is potential to transfer gene/s related to biotic resistance and medicinal importance from M. balsamina to M. charantia to produce high-quality, better yielding and stress tolerant bitter gourd genotypes. Methods: The present study provides the first and high-quality chromosome-level genome assembly of M. balsamina with size 384.90 Mb and N50 30.96 Mb using sequence data from 10x Genomics, Nanopore, and Hi-C platforms. Results: A total of 6,32,098 transposons elements; 2,15,379 simple sequence repeats; 5,67,483 transcription factor binding sites; 3,376 noncoding RNA genes; and 41,652 protein-coding genes were identified, and 4,347 disease resistance, 67 heat stress-related, 05 carotenoid-related, 15 salt stress-related, 229 cucurbitacin-related, 19 terpenes-related, 37 antioxidant activity, and 06 sex determination-related genes were characterized. Conclusion: Genome sequencing of M. balsamina will facilitate interspecific introgression of desirable traits. This information is cataloged in the form of webgenomic resource available at http://webtom.cabgrid.res.in/mbger/. Our finding of comparative genome analysis will be useful to get insights into the patterns and processes associated with genome evolution and to uncover functional regions of cucurbit genomes.
