Mapping and QTL Analysis of Gynoecy and Earliness in Bitter Gourd (Momordica charantia L.) Using Genotyping-by-Sequencing (GBS) Technology

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

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.

Genetic Diversity and Population Structure Analyses in Bitter Gourd (Momordica charantia L.) Based on Agro-Morphological and Microsatellite Markers

Bitter gourd (Momordica charantia L.) is an important vine crop of the Cucurbitaceae family and is well known for its high nutritional and medicinal values. However, the genetic variation remains largely unknown. Herein, 96 diverse bitter gourd genotypes were undertaken for diversity analysis using 10 quantitative traits, and 82 simple sequence repeat (SSR) markers. Out of 82 SSRs, 33 were polymorphic and the mean polymorphism information content (PIC) value was 0.38. Marker, JY-003 revealed a maximum (0.81) PIC value and, the number of alleles per locus ranged from 2 to 7 (average 3.46). The value of gene diversity showed the presence of a significant level of polymorphism among these genotypes. The unweighted pair group method (UPGMA) cluster analysis grouped the genotypes into two major clusters of which Cluster I comprised mostly small and medium-fruited genotypes of both M. charantia var. charantia and M. charantia var. muricata, whereas Cluster II included mostly long and extra-long fruited genotypes. Furthermore, these genotypes were divided into six distinct groups based on population structure analysis. The diversity analysis based on 10 quantitative traits revealed that earliness and high-yielding ability were exhibited by the predominantly gynoecious line DBGS-21-06 followed by DBGS-48-00. The principal component analysis (PCA) revealed that the first two components exhibited more than 50% of the total genetic variation. The present study deciphered a higher magnitude of agro-morphological and genetic diversity in 96 bitter gourd genotypes. Therefore, trait-specific genotypes identified in this study could be utilized in breeding programmes directed towards the development of improved cultivars and hybrids of bitter gourd.

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.

Molecular mapping of genomic regions and identification of possible candidate genes associated with gynoecious sex expression in bitter gourd

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 F₂ 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.

Candidate gene based SSR and SNP markers for gynoecy in bitter gourd (Momordica charantia L.)

BACKGROUND

Even though the bitter gourd hybrids are shown to have significant heterosis for many of the economic traits, processes such as manual bagging and hand pollination make the hybrid seed production labour-intensive. Use of gynoecious line as female parent makes hybrid seed production more economical. This work was performed with the objective to identify the candidate gene based molecular markers for gynoecy in bitter gourd.

METHODS AND RESULTS

Seven putative genes for flowering and sex expression, isolated from the monoecious (MC-136) and gynoecious (KAU-MCGy-101) bitter gourd accessions, were sequence characterized. MADS-box transcription factor genes AG6 and McAG2 had nucleotide polymorphisms at five sites each and were potential candidates for marker development. An In/Del polymorphism of 48 bp ([TC]₂₄) in AG6 gene was used to develop an SSR marker and a transition mutation of [A/G] in this gene was used to develop a set of SNP markers. These markers have developed distinct polymorphism between the monoecious and gynoecious genotypes and were found suited for the marker assisted selection.

CONCLUSIONS

MADS box transcription factor genes AG6 and McAG2 are identified as candidates for sex expression in bitter gourd. Based on the InDels and transition in the intronic region of AG6, SSR marker BGAG6 and an SNP marker set segregating with the sex forms were developed. The markers have been validated using four other monoecious lines and are routinely used in our bitter gourd hybrid seed production programmes.

Construction of a density mutant collection in bitter gourd via new germplasms innovation and gene functional study

Although a few studies have elucidated the creation of bitter gourd mutants, the suitable concentration and duration of ethyl methanesulfonate (EMS) mutagenesis have not been determined. In this study, mutant collection was conducted to create new germplasms and widen genetic diversity. By employing the seeds of the inbred line Y52 as the mutagenic material, EMS as the mutagen, and the suitable mutagenic conditions for bitter gourd seeds (EMS concentration 0.2%, mutagenic time 10 h), we mutated 10,000 seeds and acquired 3223 independent M₁ lines. For the randomly selected 1000 M₂ lines, 199 M₂ lines with visible phenotypes were found, and 167 M₂ lines were mutants of fruit shape, size, and tubercles. Furthermore, fourteen dwarf, eleven leaf color, five leaf shape, and eight meristem defect mutants were discovered in this mutant collection. In addition, three lines of 1253, 2284, and 3269 represented recessive mutants crossed with Y52. Furthermore, the yellow leaf lines of 2284 and 3269 were not mutated at the same gene locus. This study constructed a mutant collection through innovative new germplasms and provided valuable resources for bitter gourd breeding and functional gene research.

Epidermal Patterning Factor 2-like (McEPFL2): A Putative Candidate for the Continuous Ridge (cr) Fruit Skin Locus in Bitter Gourd (Momordica charantia L.)

Bitter gourd (Momordica charantia L.) is an economically important vegetable and medicinal crop in many Asian countries. Limited work has been conducted in understanding the genetic basis of horticulturally important traits in bitter gourd. Bitter gourd is consumed primarily for its young, immature fruit, and fruit appearance plays an important role in market acceptability. One such trait is the ridges on the fruit skin. In the present study, molecular mapping of a locus underlying fruit ridge continuity was conducted. Genetic analysis in segregating populations, derived from the crosses between two inbred lines Y1 with continuous ridges (CR) and Z-1-4 with discontinuous ridges (DCR), suggested that CR was controlled by a single recessive gene (cr). High-throughput genome sequencing of CR and DCR bulks combined with high-resolution genetic mapping in an F₂ population delimited cr into a 108 kb region with 16 predicted genes. Sequence variation analysis and expression profiling supported the epidermal patterning factor 2-like (McEPFL2) gene as the best candidate of the cr locus. A 1 bp deletion in the first exon of McEPFL2 in Y1 which would result in a truncated McEPFL2 protein may be the causal polymorphism for the phenotypic difference between Y1 and Z-1-4. The association of this 1 bp deletion with CR was further supported by gDNA sequencing of McEPFL2 among 31 bitter gourd accessions. This work provides a foundation for understanding the genetic and molecular control of fruit epidermal pattering and development, which also facilitates marker-assisted selection in bitter melon breeding.

Analysis of Genetic Diversity and Population Structure in Bitter Gourd (Momordica charantia L.) Using Morphological and SSR Markers

The present investigation was carried out using 51 diverse bitter gourd accessions as material for studying genetic diversity and relatedness using morphological and SSR markers. A wide variation was observed for morphological traits like the number of days to the first female flower anthesis (37.33–60.67), the number of days to the first fruit harvest (47.67–72.00), the number of fruits/plant (12.00–46.67), fruit length (5.00–22.23 cm), fruit diameter (1.05–6.38 cm), average fruit weight (20.71–77.67 g) and yield per plant (513.3–1976 g). Cluster analysis for 10 quantitative traits grouped the 51 accessions into 6 clusters. Out of 61 SSR primers screened, 30 were polymorphic and highly informative as a means to differentiate these accessions. Based on genotyping, a high level of genetic diversity was observed, with a total of 99 alleles. The polymorphic information content (PIC) values ranged from 0.038 for marker BG_SSR-8 to 0.721 for S-24, with an average of 0.429. The numbers of alleles ranged from 2 to 5, with an average of 3.3 alleles per locus. Gene diversity ranged from 0.04 for BG_SSR-8 to 0.76 for S-24, showing a wide variation among 51 accessions. The UPGMA cluster analysis grouped these accessions into 3 major clusters. Cluster I comprised 4 small, fruited accessions that are commercially cultivated in central and eastern India. Cluster II comprised 35 medium- to long-sized fruited accessions, which made up an abundant and diverse group. Cluster III comprised 11 long and extra-long fruited accessions. The polymorphic SSR markers of the study will be highly useful in genetic fingerprinting and mapping, and for association analysis in Momordica regarding several economic traits.

Progress on Genetic Improvement and Analysis of Bottle Gourd [Lagenaria siceraria (Molina) Standl.] for Agronomic Traits, Nutrient Compositions, and Stress Tolerance: A Review

Bottle gourd [Lagenaria siceraria (Molina) Standl.] is an important multi-purpose cucurbit crop grown for its leaf, fruit, and seed. It is widely cultivated and used for human consumption in sub-Saharan Africa (SSA) providing vital human nutrition and serving as food security crop. There is wide genetic variation among bottle gourd genetic resources in Africa for diverse qualitative and quantitative attributes for effective variety design, product development, and marketing. However, the crop is under- researched and -utilized, and improved varieties are yet to be developed and commercialized in the region. Therefore, the objective of this review is to provide the progress on bottle gourd genetic improvement and genetic analysis targeting agronomic and horticultural attributes, nutritional composition, biotic, and abiotic stress tolerance to guide current and future cultivar development, germplasm access, and conservation in SSA. The first section of the paper presents progress on breeding of bottle gourd for horticultural traits, agronomic performance, nutritional and anti-nutritional composition, and biotic and abiotic stress tolerance. This is followed by important highlights on key genetic resources of cultivated and wild bottle gourd for demand driven breeding. Lastly, the review summaries advances in bottle gourd genomics, genetic engineering and genome editing. Information presented in this paper should aid bottle gourd breeders and agronomists to develop and deploy new generation and promising varieties with farmer- and market -preferred attributes.

Genetic analysis and QTL mapping of yield and fruit traits in bitter gourd (Momordica charantia L.)

Bitter gourd (Momordica charantia L.) is an economically important vegetable crop grown in tropical parts of the world. In this study, a high-density linkage map of M. charantia was constructed through genotyping-by-sequencing (GBS) technology using F2:3 mapping population generated from the cross DBGy-201 × Pusa Do Mausami. About 2013 high-quality SNPs were assigned on a total of 20 linkage groups (LGs) spanning over 2329.2 CM with an average genetic distance of 1.16 CM. QTL analysis was performed for six major yield-contributing traits such as fruit length, fruit diameter, fruit weight, fruit flesh thickness, number of fruits per plant and yield per plant. These six quantitative traits were mapped with 19 QTLs (9 QTLs with LOD > 3) using composite interval mapping (CIM). Among 19 QTLs, 12 QTLs derived from 'Pusa Do Mausami' revealed a negative additive effect when its allele increased trait score whereas 7 QTLs derived from 'DBGy-201' revealed a positive additive effect when its allele trait score increased. The phenotypic variation (R2%) elucidated by these QTLs ranged from 0.09% (fruit flesh thickness) on LG 14 to 32.65% (fruit diameter) on LG 16 and a total of six major QTLs detected. Most QTLs detected in the present study were located relatively very close, maybe due to the high correlation among the traits. This information will serve as a significant basis for marker-assisted selection and molecular breeding in bitter gourd crop improvement.

To be a male or a female flower, a question of ethylene in cucurbits

Within the Cucurbitaceae family, most of its species develop unisexual female and male flowers, either on the same plant (monoecy) or on different plants (dioecy). As in other plant families, these two sex morphotypes have evolved from hermaphrodite species; however, many evolutionary events have occurred in cucurbits allowing easy conversion from dioecy to monoecy and vice versa. The variability in sex morphotypes is higher in the domesticated species of the family, which together with recent advances in genomics, make cucurbits an ideal model to study the genetic and molecular mechanisms that control sex determination in plants. Conventional studies demonstrated that ethylene was the master regulator of sex determination in cucurbits, although some cultivated species may respond differently to ethylene action. In this article, we survey the new advances in hormonal and genetic control of sex determination in cucurbit species, control which establishes the ethylene biosynthesis and signaling genes as being those that determine the floral meristem towards a male, female or hermaphrodite flower. The interactions between these genes are integrated into a model that explains the occurrence and distribution of unisexal and hermaphrodite flowers within the different sex morphotypes. We underline the significance of this scientific progress with regard to breeding programs for agronomically-important sex-associated traits.

Quantitative Trait Loci Mapping for Earliness, Fruit, and Seed Related Traits Using High Density Genotyping-by-Sequencing-Based Genetic Map in Bitter Gourd (Momordica charantia L.)

Bitter gourd (Momordica charantia L.) is an important vegetable crop having numerous medicinal properties. Earliness and yield related traits are main aims of bitter gourd breeding program. High resolution quantitative trait loci (QTLs) mapping can help in understanding the molecular basis of phenotypic variation of these traits and thus facilitate marker-assisted breeding. The aim of present study was to identify genetic loci controlling earliness, fruit, and seed related traits. To achieve this, genotyping-by-sequencing (GBS) approach was used to genotype 101 individuals of F₄ population derived from a cross between an elite cultivar Punjab-14 and PAUBG-6. This population was phenotyped under net-house conditions for three years 2018, 2019, and 2021. The linkage map consisting of 15 linkage groups comprising 3,144 single nucleotide polymorphism (SNP) markers was used to detect the QTLs for nine traits. A total of 50 QTLs for these traits were detected which were distributed on 11 chromosomes. The QTLs explained 5.09-29.82% of the phenotypic variance. The highest logarithm of the odds (LOD) score for a single QTL was 8.68 and the lowest was 2.50. For the earliness related traits, a total of 22 QTLs were detected. For the fruit related traits, a total of 16 QTLs and for seed related traits, a total of 12 QTLs were detected. Out of 50 QTLs, 20 QTLs were considered as frequent QTLs (FQ-QTLs). The information generated in this study is very useful in the future for fine-mapping and marker-assisted selection for these traits in bitter gourd improvement program.

High-Density GBS-Based Genetic Linkage Map Construction and QTL Identification Associated With Yellow Mosaic Disease Resistance in Bitter Gourd (Momordica charantia L.)

Yellow mosaic disease (YMD) in bitter gourd (Momordica charantia) is a devastating disease that seriously affects its yield. Although there is currently no effective method to control the disease, breeding of resistant varieties is the most effective and economic option. Moreover, quantitative trait locus (QTL) associated with resistance to YMD has not yet been reported. With the objective of mapping YMD resistance in bitter gourd, the susceptible parent "Punjab-14" and the resistant parent "PAUBG-6" were crossed to obtain F₄ mapping population comprising 101 individuals. In the present study, the genotyping by sequencing (GBS) approach was used to develop the genetic linkage map. The map contained 3,144 single nucleotide polymorphism (SNP) markers, consisted of 15 linkage groups, and it spanned 2415.2 cM with an average marker distance of 0.7 cM. By adopting the artificial and field inoculation techniques, F_4:5 individuals were phenotyped for disease resistance in Nethouse (2019), Rainy (2019), and Spring season (2020). The QTL analysis using the genetic map and phenotyping data identified three QTLs qYMD.pau_3.1, qYMD.pau_4.1, and qYMD.pau_5.1 on chromosome 3, 4, and 5 respectively. Among these, qYMD.pau_3.1, qYMD.pau_4.1 QTLs were identified during the rainy season, explaining the 13.5 and 21.6% phenotypic variance respectively, whereas, during the spring season, qYMD.pau_4.1 and qYMD.pau_5.1 QTLs were observed with 17.5 and 22.1% phenotypic variance respectively. Only one QTL qYMD.pau_5.1 was identified for disease resistance under nethouse conditions with 15.6% phenotypic variance. To our knowledge, this is the first report on the identification of QTLs associated with YMD resistance in bitter gourd using SNP markers. The information generated in this study is very useful in the future for fine-mapping and marker-assisted selection for disease resistance.