Genetic Resources and Vulnerabilities of Major Cucurbit Crops

Quantitative Trait Locus Mapping Combined with RNA Sequencing Identified Candidate Genes for Resistance to Powdery Mildew in Bitter Gourd (Momordica charantia L.)

Improving the powdery mildew resistance of bitter gourd is highly important for achieving high yield and high quality. To better understand the genetic basis of powdery mildew resistance in bitter gourd, this study analyzed 300 lines of recombinant inbred lines (RILs) formed by hybridizing the powdery mildew-resistant material MC18 and the powdery mildew-susceptible material MC402. A high-density genetic map of 1222.04 cM was constructed via incorporating 1,996,505 SNPs generated by resequencing data from 180 lines, and quantitative trait locus (QTL) positioning was performed using phenotypic data at different inoculation stages. A total of seven QTLs related to powdery mildew resistance were identified on four chromosomes, among which qPm-3-1 was detected multiple times and at multiple stages after inoculation. By selecting 18 KASP markers that were evenly distributed throughout the region, 250 lines and parents were genotyped, and the interval was narrowed to 207.22 kb, which explained 13.91% of the phenotypic variation. Through RNA-seq analysis of the parents, 11,868 differentially expressed genes (DEGs) were screened. By combining genetic analysis, gene coexpression, and sequence comparison analysis of extreme materials, two candidate genes controlling powdery mildew resistance in bitter gourd were identified (evm.TU.chr3.2934 (C3H) and evm.TU.chr3.2946 (F-box-LRR)). These results represent a step forward in understanding the genetic regulatory network of powdery mildew resistance in bitter gourd and lay a molecular foundation for the genetic improvement in powdery mildew resistance.

The mutant STAY-GREEN (Cssgr) in cucumber interacts with the CSEP30 protein to elicit a defense response against Podosphaera xanthii

Disease-resistant plants activate immune responses by specifically recognition Candidate Secreted Effector Proteins (CSEPs) through resistance (R) proteins. In research on cucumber powdery mildew resistance breeding, several R genes and CSEPs have been identified; however, the specific interactions between R proteins and CSEPs are still largely unexplored. In this study, we used a luciferase reporter assay to identify six CSEPs from Podosphaera xanthii that potentially induce cell death in cucumber. Subsequent yeast two-hybrid analysis revealed that only the mature form of CSEP30 (CSEP30∆SP) interacted with the cucumber mutant STAY-GREEN (Cssgr), a gene previously recognized for its broad-spectrum resistance in genetic studies. This interaction was confirmed using pull-down and co-immunoprecipitation assays. Additionally, to determine if the interaction leads to phenotypic changes, Cssgr and CSEP30∆SP were transiently expressed in tobacco leaves. The infiltration of Cssgr in tobacco resulted in reduced chlorosis compared to the wild-type CsSGR. Co-infiltration of Cssgr with CSEP30∆SP induced distinct dry necrotic lesions, contrasting the effects observed when Cssgr and CSEP30∆SP were infiltrated separately. Additionally, after P. xanthii infection in moderately powdery mildew-resistant Gy14 cucumber, similar necrotic lesions and specific expression of Cssgr, as along with defense response-related genes (CsPR1 and CsLecRK6.1), were observed. This study suggests that the interaction between Cssgr and CSEP30∆SP could trigger cell death and defense response, offering new insights into the molecular function of Cssgr in disease resistance in Gy14 cucumber.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-024-01504-6.

Fine genetic mapping and transcriptomic analysis revealed major gene modulating the clear stripe margin pattern of watermelon peel

The peel stripe margin pattern is one of the most important quality traits of watermelon. In this study, two contrasted watermelon lines [slb line (P1) with a clear peel stripe margin pattern and GWAS-38 line (P2) with a blurred peel stripe margin pattern] were crossed, and biparental F2 mapping populations were developed. Genetic segregation analysis revealed that a single recessive gene is modulating the main-effect genetic locus (Clcsm) of the clear stripe margin pattern of peel. Bulked segregant analysis-based sequencing (BSA-Seq) and fine genetic mapping exposed the delimited Clcsm locus to a 19.686-kb interval on chromosome 6, and the Cla97C06G126680 gene encoding the MYB transcription factor family was identified. The gene mutation analysis showed that two non-synonymous single-nucleotide polymorphism (nsSNP) sites [Chr6:28438793 (A-T) and Chr6:28438845 (A-C)] contribute to the clear peel stripe margin pattern, and quantitative real-time polymerase chain reaction (qRT-PCR) also showed a higher expression trend in the slb line than in the GWAS-38 line. Further, comparative transcriptomic analysis identified major differentially expressed genes (DEGs) in three developmental periods [4, 12, and 20 days after pollination (DAP)] of both parental lines. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment analyses indicated highly enriched DEGs involved in metabolic processes and catalytic activity. A total of 44 transcription factor families and candidate genes belonging to the ARR-B transcription factor family are believed to regulate the clear stripe margin trait of watermelon peel. The gene structure, sequence polymorphism, and expression trends depicted significant differences in the peel stripe margin pattern of both parental lines. The ClMYB36 gene showed a higher expression trend for regulating the clear peel stripe margin of the slb line, and the ClAPRR5 gene depicted a higher expression for modulating the blurred peel stripe margin in the GWAS-38 line. Overall, our fine genetic mapping and transcriptomic analysis revealed candidate genes differentiating the clear and blurred peel stripe patterns of watermelon fruit.

Widening Genetic Diversity Using Embryo Rescue in Cucurbit Crops: A Review

Embryo rescue is a vital technique in cucurbit breeding and propagation, addressing challenges such as embryo abortion, poor seed viability, and incompatibility barriers. This method involves the excision of immature embryos from seeds followed by their in vitro culture on a nutrient medium, providing an environment conducive to their growth and development. In cucurbits, embryo rescue has been extensively utilized to overcome barriers to hybridization, enabling the production of interspecific and intergeneric hybrids with desired traits. Various factors, including genotype, developmental stage of embryos, and culture conditions, influence the success of embryo rescue in cucurbits. Optimal nutrient formulations, growth regulators, and culture techniques are critical for promoting embryo germination, shoot elongation, and subsequent plantlet establishment. Additionally, embryo rescue facilitates the recovery of valuable genetic material from wild and exotic cucurbit species, expanding genetic diversity and developing novel cultivars with improved traits such as disease resistance, yield, and quality. This review highlights the principles, applications, and advancements in embryo rescue technology in cucurbits, emphasizing its significance in cucurbit breeding programs and crop improvement efforts.

Possibility of genome editing for melon breeding

Genome editing technologies are promising for conventional mutagenesis breeding, which takes a long time to remove unnecessary mutations through backcrossing and create new lines because they directly modify the target genes of elite strains. In particular, this technology has advantages for traits caused by the loss of function. Many efforts have been made to utilize this technique to introduce valuable features into crops, including maize, soybeans, and tomatoes. Several genome-edited crops have already been commercialized in the US and Japan. Melons are an important vegetable crop worldwide, produced and used in various areas. Therefore, many breeding efforts have been made to improve its fruit quality, resistance to plant diseases, and stress tolerance. Quantitative trait loci (QTL) analysis was performed, and various genes related to important traits were identified. Recently, several studies have shown that the CRISPR/Cas9 system can be applied to melons, resulting in its possible utilization as a breeding technique. Focusing on two productivity-related traits, disease resistance, and fruit quality, this review introduces the progress in genetics, examples of melon breeding through genome editing, improvements required for breeding applications, and the possibilities of genome editing in melon breeding.

Whitefly-Transmitted Viruses of Cucurbits in the Southern United States

Cucurbits are economically important crops that are widely cultivated in many parts of the world, including the southern US. In recent years, higher temperatures have favored the rapid build-up of whiteflies in the fall-grown cucurbits in this region. As a result, whitefly-transmitted viruses (WTVs) have severely impacted the marketable yield of cucurbits. In this review, we discuss three major groups of WTVs negatively impacting cucurbit cultivation in the southern US, including begomoviruses, criniviruses, and ipomoviruses. Here, we discuss the available information on the biology, epidemiology and advances made toward detecting and managing these viruses, including sources of resistance and cultural practices.

Identification of QTL associated with resistance to Phytophthora fruit rot in cucumber (Cucumis sativus L.)

Phytophthora fruit rot (PFR) caused by the soilborne oomycete pathogen, Phytophthora capsici, can cause severe yield loss in cucumber. With no resistant variety available, genetic resources are needed to develop resistant varieties. The goal of this work was to identify quantitative trait loci (QTL) associated with resistance to PFR using multiple genomic approaches and populations. Two types of resistances have been identified: age-related resistance (ARR) and young fruit resistance. ARR occurs at 12-16 days post pollination (dpp), coinciding with the end of exponential fruit growth. A major QTL for ARR was discovered on chromosome 3 and a candidate gene identified based on comparative transcriptomic analysis. Young fruit resistance, which is observed during the state of rapid fruit growth prior to commercial harvest, is a quantitative trait for which multiple QTL were identified. The largest effect QTL, qPFR5.1, located on chromosome 5 was fine mapped to a 1-Mb region. Genome-wide association studies (GWAS) and extreme-phenotype genome-wide association study (XP-GWAS) for young fruit resistance were also performed on a cucumber core collection representing > 96% of the genetic diversity of the USDA cucumber germplasm. Several SNPs overlapped with the QTL identified from QTL-seq analysis on biparental populations. In addition, novel SNPs associated with the resistance were identified from the germplasm. The resistant alleles were found mostly in accessions from India and South Asia, the center of diversity for cucumber. The results from this work can be applied to future disease resistance studies and marker-assisted selection in breeding programs.

Molecular Mapping of Putative Genomic Regions Controlling Fruit and Seed Morphology of Watermelon

The genetic regulatory basis of qualitative and quantitative phenotypes of watermelon is being investigated in different types of molecular and genetic breeding studies around the world. In this study, biparental F2 mapping populations were developed over two experimental years, and the collected datasets of fruit and seed traits exhibited highly significant correlations. Whole-genome resequencing of comparative parental lines was performed and detected single nucleotide polymorphism (SNP) loci were converted into cleaved amplified polymorphic sequence (CAPS) markers. The screened polymorphic markers were genotyped in segregating populations and two genetic linkage maps were constructed, which covered a total of 2834.28 and 2721.45 centimorgan (cM) genetic lengths, respectively. A total of 22 quantitative trait loci (QTLs) for seven phenotypic traits were mapped; among them, five stable and major-effect QTLs (PC-8-1, SL-9-1, SWi-9-1, SSi-9-1, and SW-6-1) and four minor-effect QTLs (PC-2-1 and PC-2-2; PT-2-1 and PT-2-2; SL-6-1 and SSi-6-2; and SWi-6-1 and SWi-6-2) were observed with 3.77-38.98% PVE. The adjacent QTL markers showed a good fit marker-trait association, and a significant allele-specific contribution was also noticed for genetic inheritance of traits. Further, a total of four candidate genes (Cla97C09G179150, Cla97C09G179350, Cla97C09G180040, and Cla97C09G180100) were spotted in the stable colocalized QTLs of seed size linked traits (SL-9-1 and SWi-9-1) that showed non-synonymous type mutations. The gene expression trends indicated that the seed morphology had been formed in the early developmental stage and showed the genetic regulation of seed shape formation. Hence, we think that our identified QTLs and genes would provide powerful genetic insights for marker-assisted breeding aimed at improving the quality traits of watermelon.

Classification of multi-year and multi-variety pumpkin seeds using hyperspectral imaging technology and three-dimensional convolutional neural network

BACKGROUND

Pumpkin seeds are major oil crops with high nutritional value and high oil content. The collection and identification of different pumpkin germplasm resources play a significant role in the realization of precision breeding and variety improvement. In this research, we collected 75 species of pumpkin from the Zhejiang Province of China. 35,927 near-infrared hyperspectral images of 75 types of pumpkin seeds were used as the research object.

RESULTS

To realize the rapid classification of pumpkin seed varieties, position attention embedded three-dimensional convolutional neural network (PA-3DCNN) was designed based on hyperspectral image technology. The experimental results showed that PA-3DCNN had the best classification effect than other classical machine learning technology. The classification accuracy of 99.14% and 95.20% were severally reached on the training and test sets. We also demonstrated that the PA-3DCNN model performed well in next year's classification with fine-tuning and met with 94.8% accuracy.

CONCLUSIONS

The model performance improved by introducing double convolution and pooling structure and position attention module. Meanwhile, the generalization performance of the model was verified, which can be adopted for the classification of pumpkin seeds in multiple years. This study provided a new strategy and a feasible technical approach for identifying germplasm resources of pumpkin seeds.

Development, Validation, and Application of Reverse Transcription Real-Time and Droplet Digital PCR Assays for the Detection of the Potyviruses Watermelon Mosaic Virus and Zucchini Yellow Mosaic Virus in Cucurbits

Among the cucurbit-infecting viruses, watermelon mosaic virus (WMV) and zucchini yellow mosaic virus (ZYMV) (Potyvirus: Potyviridae) are responsible for severe symptoms on cucumber, melon, watermelon, and zucchini cultivations worldwide. In this study, reverse transcription real-time PCR (real-time RT-PCR) and droplet-digital PCR (RT-ddPCR) assays targeting the coat protein (CP) genes of WMV and ZYMV were developed and validated according to the international standards of plant pest diagnosis (EPPO PM 7/98 (5)). First, the diagnostic performance of WMV-CP and ZYMV-CP real-time RT-PCRs was evaluated, and the assays displayed an analytical sensitivity of 10^-5 and 10^-3, respectively. The tests also showed an optimal repeatability, reproducibility and analytical specificity, and were reliable for the virus detection in naturally infected samples and across a wide range of cucurbit hosts. Based on these results, the real-time RT-PCR reactions were adapted to set up RT-ddPCR assays. These were the first RT-ddPCR assays aiming at the detection and quantification of WMV and ZYMV and showed a high sensitivity, being able to detect until 9 and 8 copies/µL of WMV or ZYMV, respectively. The RT-ddPCRs allowed the direct estimation of the virus concentrations and opened to a broad range of applications in disease management, such as the evaluation of partial resistance in breeding processes, identification of antagonistic/synergistic events, and studies on the implementation of natural compounds in the integrated management strategies.

Genome-Wide Characteristics of GH9B Family Members in Melon and Their Expression Profiles under Exogenous Hormone and Far-Red Light Treatment during the Grafting Healing Process

β-1,4-glucanase can not only promote the wound healing of grafted seedlings but can also have a positive effect on a plant's cell wall construction. As a critical gene of β-1,4-glucanase, GH9B is involved in cell wall remodeling and intercellular adhesion and plays a vital role in grafting healing. However, the GH9B family members have not yet been characterized for melons. In this study, 18 CmGH9Bs were identified from the melon genome, and these CmGH9Bs were located on 15 chromosomes. Our phylogenetic analysis of these CmGH9B genes and GH9B genes from other species divided them into three clusters. The gene structure and conserved functional domains of CmGH9Bs in different populations differed significantly. However, CmGH9Bs responded to cis elements such as low temperature, exogenous hormones, drought, and injury induction. The expression profiles of CmGH9Bs were different. During the graft healing process of the melon scion grafted onto the squash rootstock, both exogenous naphthyl acetic acid (NAA) and far-red light treatment significantly induced the upregulated expression of CmGH9B14 related to the graft healing process. The results provided a technical possibility for managing the graft healing of melon grafted onto squash by regulating CmGH9B14 expression.

Characterization of the USDA Cucurbita pepo, C. moschata, and C. maxima germplasm collections

The Cucurbita genus is home to a number of economically and culturally important species. We present the analysis of genotype data generated through genotyping-by-sequencing of the USDA germplasm collections of Cucurbita pepo, C. moschata, and C. maxima. These collections include a mixture of wild, landrace, and cultivated specimens from all over the world. Roughly 1,500 - 32,000 high-quality single nucleotide polymorphisms (SNPs) were called in each of the collections, which ranged in size from 314 to 829 accessions. Genomic analyses were conducted to characterize the diversity in each of the species. Analysis revealed extensive structure corresponding to a combination of geographical origin and morphotype/market class. Genome-wide associate studies (GWAS) were conducted using both historical and contemporary data. Signals were observed for several traits, but the strongest was for the bush (Bu) gene in C. pepo. Analysis of genomic heritability, together with population structure and GWAS results, was used to demonstrate a close alignment of seed size in C. pepo, maturity in C. moschata, and plant habit in C. maxima with genetic subgroups. These data represent a large, valuable collection of sequenced Cucurbita that can be used to direct the maintenance of genetic diversity, for developing breeding resources, and to help prioritize whole-genome re-sequencing.

Primary mapping of quantitative trait loci regulating multivariate horticultural phenotypes of watermelon (Citrullus lanatus L.)

Watermelon fruits exhibit a remarkable diversity of important horticultural phenotypes. In this study, we initiated a primary quantitative trait loci (QTL) mapping to identify the candidate regions controlling the ovary, fruit, and seed phenotypes. Whole genome sequencing (WGS) was carried out for two differentiated watermelon lines, and 350 Mb (96%) and 354 Mb (97%) of re-sequenced reads covered the reference de novo genome assembly, individually. A total of 45.53% non-synonymous single nucleotide polymorphism (nsSNPs) and 54.47% synonymous SNPs (sSNPs) were spotted, which produced 210 sets of novel SNP-based cleaved amplified polymorphism sequence (CAPS) markers by depicting 46.25% co-dominant polymorphism among parent lines and offspring. A biparental F_2:3 mapping population comprised of 100 families was used for trait phenotyping and CAPS genotyping, respectively. The constructed genetic map spanned a total of 2,398.40 centimorgans (cM) in length and averaged 11.42 cM, with 95.99% genome collinearity. A total of 33 QTLs were identified at different genetic positions across the eight chromosomes of watermelon (Chr-01, Chr-02, Chr-04, Chr-05, Chr-06, Chr-07, Chr-10, and Chr-11); among them, eight QTLs of the ovary, sixteen QTLs of the fruit, and nine QTLs of the seed related phenotypes were classified with 5.32-25.99% phenotypic variance explained (PVE). However, twenty-four QTLs were identified as major-effect and nine QTLs were mapped as minor-effect QTLs across the flanking regions of CAPS markers. Some QTLs were exhibited as tightly localized across the nearby genetic regions and explained the pleiotropic effects of multigenic nature. The flanking QTL markers also depicted significant allele specific contributions and accountable genes were predicted for respective traits. Gene Ontology (GO) functional enrichment was categorized in molecular function (MF), cellular components (CC), and biological process (BP); however, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were classified into three main classes of metabolism, genetic information processing, and brite hierarchies. The principal component analysis (PCA) of multivariate phenotypes widely demonstrated the major variability, consistent with the identified QTL regions. In short, we assumed that our identified QTL regions provide valuable genetic insights regarding the watermelon phenotypes and fine genetic mapping could be used to confirm them.

Identification of an emerging cucumber virus in Taiwan using Oxford nanopore sequencing technology

BACKGROUND

In June 2020, severe symptoms of leaf mosaic and fruit malformation were observed on greenhouse-grown cucumber plants in Xizhou Township of Changhua County, Taiwan. An unknown virus, designated CX-2, was isolated from a diseased cucumber sample by single lesion isolation on Chenopodium quinoa leaves. Identification of CX-2 was performed. Moreover, the incidence of cucumber viruses in Taiwan was also investigated.

METHODS

Transmission electron microscopy was performed to examine virion morphology. The portable MinION sequencer released by Oxford Nanopore Technologies was used to detect viral sequences in dsRNA of CX-2-infected leaf tissue. The whole genome sequence of CX-2 was completed by Sanger sequencing and analyzed. Reverse transcription-polymerase chain reaction (RT-PCR) with species-specific primers and indirect enzyme-linked immunosorbent assay (ELISA) with anti-coat protein antisera were developed for virus detection in the field [see Additional file 1].

RESULTS

Icosahedral particles about 30 nm in diameter were observed in the crud leaf sap of CX-2-infected C. quinoa plant. The complete genome sequence of CX-2 was determined as 4577 nt long and shared 97.0-97.2% of nucleotide identity with that of two cucumber Bulgarian latent virus (CBLV) isolates in Iran and Bulgaria. Therefore, CX-2 was renamed CBLV-TW. In 2020-2022 field surveys, melon yellow spot virus (MYSV) had the highest detection rate of 74.7%, followed by cucurbit chlorotic yellows virus (CCYV) (32.0%), papaya ringspot virus virus watermelon type (PRSV-W) (10.7%), squash leaf curl Philippines virus (SLCuPV) (9.3%), CBLV (8.0%) and watermelon silver mottle virus (WSMoV) (4.0%). Co-infection of CBLV and MYSV could be detected in field cucumbers.

CONCLUSION

The emerging CBLV-TW was identified by nanopore sequencing. Whole genome sequence analysis revealed that CBLV-TW is closely related, but phylogenetically distinct, to two known CBLV isolates in Bulgaria and Iran. Detection methods including RT-PCR and indirect ELISA have been developed to detect CBLV and to investigate cucumber viruses in central Taiwan. The 2020-2022 field survey results showed that MYSV and CCYV were the main threats to cucumbers, with CBLV, SLCuPV and WSMoV were occasionally occurring.

Edible Vitalmelon Fruit Extract Inhibits Adipogenesis and Ameliorates High-Fat Diet-Induced Obesity

Conventional breeding of wild (Cucumis melo var. makuwa Makino (CM)) and cultivated (Cucumis melo var. reticulatus (CR)) melons is aimed at improving their biological traits. Here, we prepared a nontoxic, bioactive extract of vitalmelon (F1 hybrid) and evaluated its antiadipogenic and antiobesity effects in fully differentiated 3T3-L1 adipocytes and high-fat diet- (HFD-) induced obese C57BL/6 mice. In fully differentiated 3T3-L1 adipocytes, the vitalmelon extract reduced the DMI- (dexamethasone, 3-isobutyl-1-methylxanthine, and insulin-) induced increases in lipid droplet number and intracellular glucose and triglyceride levels. In addition, the extract inhibited 3T3-L1 preadipocyte differentiation by downregulating PPAR-γ and target genes LPL, CD36, HMGCR, and L-FABP. To investigate the inhibitory effects of the vitalmelon extract on lipid metabolism, we measured serum lipid, hormone, and cytokine concentrations; lipolytic activity; lipid accumulation; and adipogenesis in HFD-fed mice treated with the extract. The HFD+vitalmelon-fed mice showed lower blood cholesterol, free fatty acid, sugar, leptin, and insulin concentrations but higher blood adiponectin concentrations than the HFD-fed mice. Moreover, the HFD+vitalmelon-fed mice showed lower abdominal fat levels, smaller fat cells, lower weight, and fewer lipid droplets in the liver tissue than the HFD-fed mice. Therefore, in HFD-fed mice, vitalmelon regulated lipid metabolism through PPAR-γ, highlighting its potential as a promising antiobesity functional food.

Development of Whole Genome SNP-CAPS Markers and Preliminary QTL Mapping of Fruit Pedicel Traits in Watermelon

Fruit pedicel (FP) is an important determinant of premium fruit quality that directly affects commercial market value. However, in-depth molecular and genetic basis of pedicel-related traits has not been identified in watermelon. Herein, a quantitative trait locus (QTL) mapping strategy was used to identify the potential genetic regions controlling FP traits based on newly derived whole-genome single nucleotide polymorphism based cleaved amplified polymorphism sequence (SNP-CAPS) markers. Next-generation sequencing based whole-genome re-sequencing of two watermelon parent lines revealed 98.30 and 98.40% of average coverage, 4,989,869 SNP variants, and 182,949 CAPS loci pairs across the reference genome, respectively. A total of 221 sets of codominant markers exhibited 46.42% polymorphism rate and were effectively genotyped within 100-F_2:3 derived mapping population. The developed linkage map covered a total of 2,630.49 cM genetic length with averaged 11.90 cM, and depicted a valid marker-trait association. In total, 6 QTLs (qFPL4.1, qFPW4.1, qFPD2.1, qFPD2.2, qFPD8.1, qFPD10.1) were mapped with five major effects and one minor effect between the whole genome adjacent markers positioned over distinct chromosomes (02, 04, 08, 10), based on the ICIM-ADD mapping approach. These significant QTLs were similarly mapped in delimited flanking regions of 675.10, 751.38, 859.24, 948.39, and 947.51 kb, which collectively explained 8.64-13.60% PVE, respectively. A highly significant and positive correlation was found among the observed variables. To our knowledge, we first time reported the mapped QTLs/genes affecting FP traits of watermelon, and our illustrated outcomes will deliver the potential insights for fine genetic mapping as well as functional gene analysis through MAS-based breeding approaches.

Transcriptomic Analysis of Fusarium oxysporum Stress-Induced Pathosystem and Screening of Fom-2 Interaction Factors in Contrasted Melon Plants

Fusarium wilt is one of the most destructive and less controllable diseases in melon, which is usually caused by fusarium oxysporum. In this study, transcriptome sequencing and Yeast Two-Hybrid (Y2H) methods were used for quantification of differentially expressed genes (DEGs) involved in fusarium oxysporum (f. sp. melonis race 1) stress-induced mechanisms in contrasted melon varieties (M4-45 "susceptible" and MR-1 "resistant"). The interaction factors of Fom-2 resistance genes were also explored in response to the plant-pathogen infection mechanism. Transcriptomic analysis exhibited total 1,904 new genes; however, candidate DEGs analysis revealed a total of 144 specific genes (50 upregulated and 94 downregulated) for M4-45 variety and 104 specific genes (71 upregulated and 33 downregulated) for MR-1 variety, respectively. The analysis of Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway depicted some candidate DEGs, including Phenylalanine metabolism, phenylpropane biosynthesis, plants-pathogen interaction, and signal transduction of plant hormones, which were mainly involved in disease resistance metabolic pathways. The weighted gene co-expression network analysis (WGCNA) analysis revealed a strong correlation module and exhibited the disease resistance-related genes encoding course proteins, transcription factors, protein kinase, benzene propane biosynthesis path, plants-pathogen interaction pathway, and glutathione S-transferase. Meanwhile, the resistance-related specific genes expression was relatively abundant in MR-1 compared to the M4-45, and cell wall-associated receptor kinases (MELO3C008452 and MELO3C008453), heat shock protein (Cucumis_melo_newGene_172), defensin-like protein (Cucumis_melo_newGene_5490), and disease resistance response protein (MELO3C016325), activator response protein (MELO3C021623), leucine-rich repeat receptor protein kinase (MELO3C024412), lactyl glutathione ligase (Cucumis_melo_newGene_36), and unknown protein (MELO3C007588) were persisted by exhibiting the upregulated expressions. At the transcription level, the interaction factors between the candidate genes in response to the fusarium oxysporum induced stress, and Y2H screening signified the main contribution of MYB transcription factors (MELO3C009678 and MELO3C014597), BZIP (MELO3C011839 and MELO3C019349), unknown proteins, and key enzymes in the ubiquitination process (4XM334FK014). The candidate genes were further verified in exogenously treated melon plants with f. oxysporum (Fom-2, Race 1), Abscisic acid (ABA), Methyl Jasmonite (MeJA), and Salicylic acid (SA), using the fluorescence quantitative polymerase chain reaction (qRT-PCR) analysis. The overall expression results indicated that the SA signal pathway is involved in effective regulation of the Fom-2 gene activity.