Integrated full-length transcriptome and metabolome analysis reveals the defence response of melon to gummy stem blight

Gummy stem blight (GSB), a widespread disease causing great loss to cucurbit production, has become a major threat to melon cultivation. However, the melon-GSB interaction remains largely unknown. Here, full-length transcriptome and widely targeted metabolome were used to investigate the defence responses of resistant (PI511089) and susceptible (Payzawat) melon accessions to GSB pathogen infection at 24 h. The biosynthesis of secondary metabolites and MAPK signalling pathway were specifically enriched for differentially expressed genes in PI511890, while carbohydrate metabolism and amino acid metabolism were specifically enriched in Payzawat. More than 1000 novel genes were identified and MAPK signalling pathway was specifically enriched for them in PI511890. There were 11 793 alternative splicing events involving in the defence response to GSB. Totally, 910 metabolites were identified in Payzawat and PI511890, and flavonoids were the dominant metabolites. Integrated full-length transcriptome and metabolome analysis showed eriodictyol and oxalic acid were the potential marker metabolites for GSB resistance in melon. Moreover, posttranscription regulation was widely involved in the defence response of melon to GSB pathogen infection. These results not only improve our understanding on the interaction between melon and GSB, but also facilitate the genetic improvement of melon with GSB resistance.

Sequence characteristics, genetic diversity and phylogenetic analysis of the Cucurbita ficifolia (Cucurbitaceae) chloroplasts genome

BACKGROUND

Curcubita ficifolia Bouché (Cucurbitaceae) has high value as a food crop and medicinal plant, and also has horticultural value as rootstock for other melon species. China is home to many different cultivars, but the genetic diversity of these resources and the evolutionary relationships among them, as well as the differences between C. ficifolia and other Cucurbita species, remain unclear.

RESULTS

We investigated the chloroplast (cp) genomes of 160 C. ficifolia individuals from 31 populations in Yunnan, a major C. ficifolia production area in China. We found that the cp genome of C. ficifolia is ~151 kb and contains 128 genes, of which 86 are protein coding genes, 34 encode tRNA, and eight encode rRNAs. We also identified 64 SSRs, mainly AT repeats. The cp genome was found to contain a total of 204 SNP and 57 indels, and a total of 21 haplotypes were found in the 160 study individuals. The reverse repeat (IR) region of C. ficifolia contained a few differences compared with this region in the six other Cucurbita species. Sequence difference analysis demonstrated that most of the variable regions were concentrated in the single copy (SC) region. Moreover, the sequences of the coding regions were found to be more similar among species than those of the non-coding regions. The phylogenies reconstructed from the cp genomes of 61 representative species of Cucurbitaceae reflected the currently accepted classification, in which C. ficifolia is sister to the other Cucurbita species, however, different interspecific relationships were found between Cucurbita species.

CONCLUSIONS

These results will be valuable in the classification of C. ficifolia genetic resources and will contribute to our understanding of evolutionary relationships within the genus Cucurbita.

Developmental stages and episode-specific regulatory genes in andromonoecious melon flower development

BACKGROUND AND AIMS

Given the lack of specific studies on floral development in melon (Cucumis melo L.), we carried out an extensive study involving morphological and transcriptomic analyses to characterize floral development in this species.

METHODS

Using an andromonoecious line, we analysed the development of floral buds in male and hermaphrodite flowers with both light microscopy and scanning electron microscopy. Based on flower lengths, we established a correlation between the developmental stages and four main episodes of floral development and conducted an extensive RNA sequencing analysis of these episodes.

KEY RESULTS

We identified 12 stages of floral development, from the appearance of the floral meristems to anthesis. The main structural differences between male and hermaphrodite flowers appeared between stages 6 and 7; later stages of development leading to the formation of organs and structures in both types of flowers were also described. We analysed the gene expression patterns of the four episodes in flower development to find the genes that were specific to each given episode. Among others, we identified genes that defined the passage from one episode to the next according to the ABCDE model of floral development.

CONCLUSIONS

This work combines a detailed morphological analysis and a comprehensive transcriptomic study to enable characterization of the structural and molecular mechanisms that determine the floral development of an andromonoecious genotype in melon. Taken together, our results provide a first insight into gene regulation networks in melon floral development that are crucial for flowering and pollen formation, highlighting potential targets for genetic manipulation to improve crop yield of melon in the future.

Transcriptomic analysis provides an insight into the function of CmGH9B3, a key gene of β-1, 4-glucanase, during the graft union healing of oriental melon scion grafted onto squash rootstock

The melon (Cucumis melo L.) is a globally cherished and economically significant crop. The grafting technique has been widely used in the vegetative propagation of melon to promote environmental tolerance and disease resistance. However, mechanisms governing graft healing and potential incompatibilities in melons following the grafting process remain unknown. To uncover the molecular mechanism of healing of grafted melon seedlings, melon wild type (Control) and TRV-CmGH9B3 lines were obtained and grafted onto the squash rootstocks (C. moschata). Anatomical differences indicated that the healing process of the TRV-CmGH9B3 plants was slower than that of the control. A total of 335 significantly differentially expressed genes (DEGs) were detected between two transcriptomes. Most of these DEGs were down-regulated in TRV-CmGH9B3 grafted seedlings. GO and KEGG analysis showed that many metabolic, physiological, and hormonal responses were involved in graft healing, including metabolic processes, plant hormone signaling, plant MAPK pathway, and sucrose starch pathway. During the healing process of TRV-CmGH9B3 grafted seedlings, gene synthesis related to hormone signal transduction (auxin, cytokinin, gibberellin, brassinolide) was delayed. At the same time, it was found that most of the DEGs related to the sucrose pathway were down-regulated in TRV-CmGH9B3 grafted seedlings. The results showed that sugar was also involved in the healing process of melon grafted onto squash. These results deepened our understanding of the molecular mechanism of GH9B3, a key gene of β-1, 4-glucanase. It also provided a reference for elucidating the gene mechanism and function analysis of CmGH9B3 in the process of graft union healing.

Genome-Wide Analysis of the LBD Gene Family in Melon and Expression Analysis in Response to Wilt Disease Infection

LBD transcription factors are a class of transcription factors that regulate the formation of lateral organs, establish boundaries, and control secondary metabolism in plants. In this study, we identified 37 melon LBD transcription factors using bioinformatics methods and analyzed their basic information, chromosomal location, collinearity, evolutionary tree, gene structure, and expression patterns. The results showed that the genes were unevenly distributed across the 13 chromosomes of melon plants, with tandem repeats appearing on chromosomes 11 and 12. These 37 transcription factors can be divided into two major categories, Class I and Class II, and seven subfamilies: Ia, Ib, Ic, Id, Ie, IIa, and IIb. Of the 37 included transcription factors, 25 genes each contained between one to three introns, while the other 12 genes did not contain introns. Through cis-acting element analysis, we identified response elements such as salicylic acid, MeJA, abscisic acid, and auxin, gibberellic acid, as well as light response, stress response, and MYB-specific binding sites. Expression pattern analysis showed that genes in the IIb subfamilies play important roles in the growth and development of various organs in melon plants. Expression analysis found that the majority of melon LBD genes were significantly upregulated after infection with wilt disease, with the strongest response observed in the stem.

Reliable reference gene selection for quantitative real-time PCR (qRT-PCR) in floral developmental phases of dioecious species Coccinia grandis

The flowering process is intricate and regulated by a combination of external and internal factors. Delving into gene expression research has the potential to enhance our comprehension of the molecular foundations underlying floral development. Because of its accuracy, specificity, reproducibility, and efficiency, qRT-PCR is now a biological research tool for studying expression pattern of desired genes. The gene expression investigations using qRT-PCR required a reference gene with relatively uniform expression levels in multiple biological samples, including different developmental stages, tissues, and experimental conditions. In this study, experimental sets offloral and floral organ development in the male and female plants of C. grandis, a dioecious Cucurbitaceae species, qRT-PCR profiling was performed using six reference genes as internal control with B-class floral identity gene, PISTILLATA (PI). To analyse the data, algorithms such as geNorm, NormFinder, RefFinder, and BestKeeper were used to pick out the best internal controls from a group of candidates. The optimal reference gene for qRT-PCR studies with floral samples has been recommended as β-actin combined with β-tubulin. This is the first report on the validation of candidate reference genes across flower developmental stages in the dioecious species C. grandis, which will provide basic data for research on the molecular mechanism underlying flower development in this species and lay the groundwork for similar studies in other related species.

Identification and characterization analysis of the HDM gene family in melon (Cucumis melo L.)

Histone demethylase (HDM) play crucial roles in regulating plant growth and environmental adaptation. In this study, the HDM gene family in melon was identified by bioinformatics methods and the expression patterns of the CmHDM family members in different melon tissues were analyzed using transcriptome data. The results showed that 20 CmHDM genes were identified in the melon genome, which were unevenly distributed across each chromosome. These members fall into two major categories: LSD1 and JmjC. The JmjC group could be further divided into five subgroups with different numbers. The results of collinearity analysis of intraspecific and interspecific relationships showed that there were only one pair of segmental duplication in melon HDM genes, and more collinearity in genetic relationship of HDM genes between melon and tomato. The numbers of conserved domains, exons and introns in each member vary and various cis-acting elements responding to hormones and environmental signals existed in the respective promoter regions. Expression analysis showed that the respective gene members were expressed at different levels in male flowers, female flowers, roots, stems, leaves, ovary, and mature fruits of melon. These results will contribute to the understanding on the potential functions of the HDM genes and their potential functions in regulating melon growth and environmental adaptation.

The maturation profile triggers differential expression of sugar metabolism genes in melon fruits

Melons are commercially important crops that requires specific quality attributes for successful commercialization, including accumulation of sugars, particularly sucrose. This trait can be influenced by various factors, such as the type of ripening. Cucumis melo L. is an ideal species for studying sugar metabolism because it has both climacteric and non-climacteric cultivars. Thus, this study aimed to examine the gene expression of sucrose metabolism candidates using RT-qPCR, in conjunction with postharvest physiological analyzes and high-performance liquid chromatography-based sugar quantification, in the melon cultivars 'Gaúcho' (climacteric) and 'Eldorado' (non-climacteric). The results showed that sucrose synthase 1 played a role in the synthesis and accumulation of sucrose in both cultivars, whereas sucrose synthase 2 was more highly expressed in 'Gaúcho', contributing to lower hexose content. Invertase inhibitor 1 was more highly expressed in 'Eldorado' and may be involved in sugar-induced maturation. Neutral α-galactosidase had distinct functions, playing a role in substrate synthesis for the growth of young 'Eldorado' fruits, whereas in mature 'Gaúcho' fruits it participated in the metabolism of raffinose family oligosaccharides for sucrose accumulation. The expression of trehalose-6-phosphate synthase genes indicated a greater involvement of these enzymes in the sugar regulation in 'Gaúcho' melons. These findings shed light on the intraspecific differences related to fruit quality attributes in different types of maturation and contribute to a deeper understanding of the underlying molecular mechanisms involved in the metabolism of sugars in melons, which can inform breeding programs aimed at improving fruit quality attributes in this crop.

Insights into the early transcriptomic response against watermelon mosaic virus in melon

BACKGROUND

Watermelon mosaic virus (WMV) is one of the most prevalent viruses affecting melon worldwide. Recessive resistance to WMV in melon has previously been reported in the African accession TGR-1551. Moreover, the genomic regions associated to the resistance have also been described. Nevertheless, the transcriptomic response that might infer the resistance to this potyvirus has not been explored.

RESULTS

We have performed a comparative transcriptomic analysis using mock and WMV-inoculated plants of the susceptible cultivar "Bola de oro" (BO) and a resistant RIL (Recombinant inbred line) derived from the initial cross between "TGR-1551" and BO. In total, 616 genes were identified as differentially expressed and the weighted gene co-expression network analysis (WGCNA) detected 19 gene clusters (GCs), of which 7 were differentially expressed for the genotype x treatment interaction term. SNPs with a predicted high impact on the protein function were detected within the coding regions of most of the detected DEGs. Moreover, 3 and 16 DEGs were detected within the QTL regions previously described in chromosomes 11 and 5, respectively. In addition to these two specific genomic regions, we also observde large transcriptomic changes from genes spread across the genome in the resistant plants in response to the virus infection. This early response against WMV implied genes involved in plant-pathogen interaction, plant hormone signal transduction, the MAPK signaling pathway or ubiquitin mediated proteolysis, in detriment to the photosynthetic and basal metabolites pathways. Moreover, the gene MELO3C021395, which coded a mediator of RNA polymerase II transcription subunit 33A (MED33A), has been proposed as the candidate gene located on chromosome 11 conferring resistance to WMV.

CONCLUSIONS

The comparative transcriptomic analysis presented here showed that, even though the resistance to WMV in TGR-1551 has a recessive nature, it triggers an active defense response at a transcriptomic level, which involves broad-spectrum resistance mechanisms. Thus, this study represents a step forward on our understanding of the mechanisms underlaying WMV resistance in melon. In addition, it sheds light into a broader topic on the mechanisms of recessive resistances.

Comparative Analysis of Transcriptomics and Metabolomics Reveals Defense Mechanisms in Melon Cultivars against Pseudoperonospora cubensis Infection

Melon (Cucumis melo L.) represents an agriculturally significant horticultural crop that is widely grown for its flavorful fruits. Downy mildew (DM), a pervasive foliar disease, poses a significant threat to global melon production. Although several quantitative trait loci related to DM resistance have been identified, the comprehensive genetic underpinnings of this resistance remain largely uncharted. In this study, we utilized integrative transcriptomics and metabolomics approaches to identify potential resistance-associated genes and delineate the strategies involved in the defense against DM in two melon cultivars: the resistant 'PI442177' ('K10-1') and the susceptible 'Huangdanzi' ('K10-9'), post-P. cubensis infection. Even in the absence of the pathogen, there were distinctive differentially expressed genes (DEGs) between 'K10-1' and 'K10-9'. When P. cubensis was infected, certain genes, including flavin-containing monooxygenase (FMO), receptor-like protein kinase FERONIA (FER), and the HD-ZIP transcription factor member, AtHB7, displayed pronounced expression differences between the cultivars. Notably, our data suggest that following P. cubensis infection, both cultivars suppressed flavonoid biosynthesis via the down-regulation of associated genes whilst concurrently promoting lignin production. The complex interplay of transcriptomic and metabolic responses elucidated by this study provides foundational insights into melon's defense mechanisms against DM. The robust resilience of 'K10-1' to DM is attributed to the synergistic interaction of its inherent transcriptomic and metabolic reactions.

The effect of environmental factors on the genetic differentiation of Cucurbita ficifolia populations based on whole-genome resequencing

BACKGROUND

Cucurbita ficifolia is one of the squash species most resistant to fungal pathogens, and has especially high resistance to melon Fusarium wilt. This species is therefore an important germplasm resource for the breeding of squash and melon cultivars.

RESULTS

Whole-genome resequencing of 223 individuals from 32 populations in Yunnan Province, the main cucurbit production area in China, was performed and 3,855,120 single-nucleotide polymorphisms (SNPs) and 1,361,000 InDels were obtained. SNP analysis suggested that levels of genetic diversity in C. ficifolia were high, but that different populations showed no significant genetic differentiation or geographical structure, and that individual C. ficifolia plants with fruit rinds of a similar color did not form independent clusters. A Mantel test conducted in combination with geographical distance and environmental factors suggested that genetic distance was not correlated with geographical distance, but had a significant correlation with environmental distance. Further associations between the genetic data and five environmental factors were analyzed using whole-genome association analysis. SNPs associated with each environmental factor were investigated and genes 250 kb upstream and downstream from associated SNPs were annotated. Overall, 15 marker-trait-associated SNPs (MTAs) and 293 genes under environmental selection were identified. The identified genes were involved in cell membrane lipid metabolism, macromolecular complexes, catalytic activity and other related aspects. Ecological niche modeling was used to simulate the distribution of C. ficifolia across time, from the present and into the future. We found that the area suitable for C. ficifolia changed with the changing climate in different periods.

CONCLUSIONS

Resequencing of the C. ficifolia accessions has allowed identification of genetic markers, such as SNPs and InDels. The SNPs identified in this study suggest that environmental factors mediated the formation of the population structure of C. ficifolia in China. These SNPs and Indels might also contribute to the variation in important pathways of genes for important agronomic traits such as yield, disease resistance and stress tolerance. Moreover, the genome resequencing data and the genetic markers identified from 223 accessions provide insight into the genetic variation of the C. ficifolia germplasm and will facilitate a broad range of genetic studies.

Genome-wide identification and expression analysis of the JMJ-C gene family in melon (Cucumis melo L.) reveals their potential role in fruit development

BACKGROUND

Proteins with the jumonji (JMJ)-C domain belong to the histone demethylase family and contribute to reverse histone methylation. Although JMJ-C family genes have an essential role in regulating plant growth and development, the characterization of the JMJ-C family genes in melon has not been uncovered.

RESULTS

In this study, a total of 17 JMJ-C proteins were identified in melon (Cucumis melo L.). CmJMJs were categorized into five subfamilies based on the specific conserved domain: KDM4/JHDM3, KDM5/JARID1, JMJD6, KDM3/JHDM2, and JMJ-C domain-only. The chromosome localization analyses showed that 17 CmJMJs were distributed on nine chromosomes. Cis-acting element analyses of the 17 CmJMJ genes showed numerous hormone, light, and stress response elements distributed in the promoter region. Covariance analysis revealed one pair of replicated fragments (CmJMJ3a and CmJMJ3b) in 17 CmJMJ genes. We investigated the expression profile of 17 CmJMJ genes in different lateral organs and four developmental stages of fruit by RNA-seq transcriptome analysis and RT-qPCR. The results revealed that most CmJMJ genes were prominently expressed in female flowers, ovaries, and developing fruits, suggesting their active role in melon fruit development. Subcellular localization showed that the fruit-related CmJMJ5a protein is specifically localized in the cell nucleus.

CONCLUSIONS

This study provides a comprehensive understanding of the gene structure, classification, and evolution of JMJ-C in melon and supports the clarification of the JMJ-C functions in further research.

Resolving the phylogeny of Thladiantha (Cucurbitaceae) with three different target capture pipelines

BACKGROUND

Despite recent advances, reliable tools to simultaneously handle different types of sequencing data (e.g., target capture, genome skimming) for phylogenomics are still scarce. Here, we evaluate the performance of the recently developed pipeline Captus in comparison with the well-known target capture pipelines HybPiper and SECAPR. As test data, we analyzed newly generated sequences for the genus Thladiantha (Cucurbitaceae) for which no well-resolved phylogeny estimate has been available so far, as well as simulated reads derived from the genome of Arabidopsis thaliana.

RESULTS

Our pipeline comparisons are based on (1) the time needed for data assembly and locus extraction, (2) locus recovery per sample, (3) the number of informative sites in nucleotide alignments, and (4) the topology of the nuclear and plastid phylogenies. Additionally, the simulated reads derived from the genome of Arabidopsis thaliana were used to evaluate the accuracy and completeness of the recovered loci. In terms of computation time, locus recovery per sample, and informative sites, Captus outperforms HybPiper and SECAPR. The resulting topologies of Captus and SECAPR are identical for coalescent trees but differ when trees are inferred from concatenated alignments. The HybPiper phylogeny is similar to Captus in both methods. The nuclear genes recover a deep split of Thladiantha in two clades, but this is not supported by the plastid data.

CONCLUSIONS

Captus is the best choice among the three pipelines in terms of computation time and locus recovery. Even though there is no significant topological difference between the Thladiantha species trees produced by the three pipelines, Captus yields a higher number of gene trees in agreement with the topology of the species tree (i.e., fewer genes in conflict with the species tree topology).

Simple, fast and inexpensive hot sodium hydroxide and tris DNA extraction method for genotyping tomato and melon seeds

Seed commerce is a highly profitable global market. Most commercialized seeds are hybrid seeds originating from a controlled cross between two selected parental lines. The market value of hybrid seeds depends on their hybrid genetic purity. DNA molecular markers are a reliable and widespread tool to genotype plant materials; however, DNA extraction from seeds is challenging, often laborious and expensive. With the ultimate goal of creating a tomato and melon hybrid seeds purity test, various challenges arise. To overcome these problems and with the purpose of crude DNA extraction, a simple, fast, inexpensive and easily scalable adaptation of the hot sodium hydroxide and tris method coupled to a competitive allele-specific PCR genotyping method is proposed.

Embryonic microinjection of ribonucleoprotein complex (Cas9+sgRNA) of white gene in melon fly, Zeugodacus cucurbitae (Coquillett) (Diptera: Tephritidae) produced white eye phenotype

Melon fly, Zeugodacus cucurbitae (Coquillett) is a major pest of cucurbitaceous crops, and causes substantial yield losses and economic costs. CRISPR/Cas9 is a rapid and effective site-specific genome editing tool for the generation of genetic changes that are stable and heritable. The CRISPR/Cas9 tool uses synthetically designed single guide RNA (sgRNA) that is complementary to the target gene and guides the Cas9 enzyme to perform nuclease activity by making double-strand breaks in the target DNA sequences. This tool can be effectively exploited to improve traits critical for the management of insect pests by targeting specific genes encoding these traits without the need of extensive genetic information. The white gene is an important gene responsible for the transport of body pigment precursor molecules. In this study, we produced effective mutagenesis of the white gene of Z. cucurbitae using the CRISPR/Cas9 tool with double sgRNA to target multiple sites of white to increase the efficiency in the generation of frame-shift mutations resulting in the white eye phenotype in adults. This was achieved through embryonic microinjection of the ribonucleoprotein (RNP) complex in the pre-blastoderm embryo stage 1 h after embryo laying. Our success with the production of a white eye mutant fly by CRISPR/Cas9 mutagenesis is important for the research on gene function and protein-level modifications in melon fly and forms the basis for the development of new genetic control strategies such as precision guided sterile insect technique (pgSIT) for this pest of economic significance.

Molecular Marker-Assisted Mapping, Candidate Gene Identification, and Breeding in Melon (Cucumis melo L.): A Review

Melon (Cucumis melo L.) is an important crop that is cultivated worldwide for its fleshy fruit. Understanding the genetic basis of a plant's qualitative and quantitative traits is essential for developing consumer-favored varieties. This review presents genetic and molecular advances related to qualitative and quantitative phenotypic traits and biochemical compounds in melons. This information guides trait incorporation and the production of novel varieties with desirable horticultural and economic characteristics and yield performance. This review summarizes the quantitative trait loci, candidate genes, and development of molecular markers related to plant architecture, branching patterns, floral attributes (sex expression and male sterility), fruit attributes (shape, rind and flesh color, yield, biochemical compounds, sugar content, and netting), and seed attributes (seed coat color and size). The findings discussed in this review will enhance demand-driven breeding to produce cultivars that benefit consumers and melon breeders.

Ethylene produced in carpel primordia controls CmHB40 expression to inhibit stamen development

Sex determination evolved to control the development of unisexual flowers. In agriculture, it conditions how plants are cultivated and bred. We investigated how female flowers develop in monoecious cucurbits. We discovered in melon, Cucumis melo, a mechanism in which ethylene produced in the carpel is perceived in the stamen primordia through spatially differentially expressed ethylene receptors. Subsequently, the CmEIN3/CmEIL1 ethylene signalling module, in stamen primordia, activates the expression of CmHB40, a transcription factor that downregulates genes required for stamen development and upregulates genes associated with organ senescence. Investigation of melon genetic biodiversity revealed a haplotype, originating in Africa, altered in EIN3/EIL1 binding to CmHB40 promoter and associated with bisexual flower development. In contrast to other bisexual mutants in cucurbits, CmHB40 mutations do not alter fruit shape. By disentangling fruit shape and sex-determination pathways, our work opens up new avenues in plant breeding.

A Genome-Wide Identification and Comparative Analysis of the Heavy-Metal-Associated Gene Family in Cucurbitaceae Species and Their Role in Cucurbita pepo under Arsenic Stress

The heavy-metal-associated (HMA) proteins are a class of PB1-type ATPases related to the intracellular transport and detoxification of metals. However, due to a lack of information regarding the HMA gene family in the Cucurbitaceae family, a comprehensive genome-wide analysis of the HMA family was performed in ten Cucurbitaceae species: Citrullus amarus, Citrullus colocynthis, Citrullus lanatus, Citrullus mucosospermus, Cucumis melo, Cucumis sativus, Cucurbita maxima, Cucurbita moschata, Cucurbita pepo, and Legenaria siceraria. We identified 103 Cucurbit HMA proteins with various members, ranging from 8 (Legenaria siceraria) to 14 (Cucurbita pepo) across species. The phylogenetic and structural analysis confirmed that the Cucurbitaceae HMA protein family could be further classified into two major clades: Zn/Co/Cd/Pb and Cu/Ag. The GO-annotation-based subcellular localization analysis predicted that all HMA gene family members were localized on membranes. Moreover, the analysis of conserved motifs and gene structure (intron/exon) revealed the functional divergence between clades. The interspecies microsynteny analysis demonstrated that maximum orthologous genes were found between species of the Citrullus genera. Finally, nine candidate HMA genes were selected, and their expression analysis was carried out via qRT-PCR in root, leaf, flower, and fruit tissues of C. pepo under arsenic stress. The expression pattern of the CpeHMA genes showed a distinct pattern of expression in root and shoot tissues, with a remarkable expression of CpeHMA6 and CpeHMA3 genes from the Cu/Ag clade. Overall, this study provides insights into the functional analysis of the HMA gene family in Cucurbitaceae species and lays down the basic knowledge to explore the role and mechanism of the HMA gene family to cope with arsenic stress conditions.

Pan-genome analysis sheds light on structural variation-based dissection of agronomic traits in melon crops

Sweetness and appearance of fresh fruits are key palatable and preference attributes for consumers and are often controlled by multiple genes. However, fine-mapping the key loci or genes of interest by single genome-based genetic analysis is challenging. Herein, we present the chromosome-level genome assembly of 1 landrace melon accession (Cucumis melo ssp. agrestis) with wild morphologic features and thus construct a melon pan-genome atlas via integrating sequenced melon genome datasets. Our comparative genomic analysis reveals a total of 3.4 million genetic variations, of which the presence/absence variations (PAVs) are mainly involved in regulating the function of genes for sucrose metabolism during melon domestication and improvement. We further resolved several loci that are accountable for sucrose contents, flesh color, rind stripe, and suture using a structural variation (SV)-based genome-wide association study. Furthermore, via bulked segregation analysis (BSA)-seq and map-based cloning, we uncovered that a single gene, (CmPIRL6), determines the edible or inedible characteristics of melon fruit exocarp. These findings provide important melon pan-genome information and provide a powerful toolkit for future pan-genome-informed cultivar breeding of melon.

Comprehensive analyses of microtubule-associated protein MAP65 family genes in Cucurbitaceae and CsaMAP65s expression profiles in cucumber

MAP65 is a microtubule-binding protein family in plants and plays crucial roles in regulating cell growth and development, intercellular communication, and plant responses to various environmental stresses. However, MAP65s in Cucurbitaceae are still less understood. In this study, a total of 40 MAP65s were identified from six Cucurbitaceae species (Cucumis sativus L., Citrullus lanatus, Cucumis melo L., Cucurbita moschata, Lagenaria siceraria, and Benincasa hispida) and classified into five groups by phylogenetic analysis according to gene structures and conserved domains. A conserved domain (MAP65_ASE1) was found in all MAP65 proteins. In cucumber, we isolated six CsaMAP65s with different expression patterns in tissues including root, stem, leaf, female flower, male flower, and fruit. Subcellular localizations of CsaMAP65s verified that all CsaMAP65s were localized in microtubule and microfilament. Analyses of the promoter regions of CsaMAP65s have screened different cis-acting regulatory elements involved in growth and development and responses to hormone and stresses. In addition, CsaMAP65-5 in leaves was significantly upregulated by salt stress, and this promotion effect was higher in cucumber cultivars with salt tolerant than that without salt tolerant. CsaMAP65-1 in leaves was significantly upregulated by cold stress, and this promotion was higher in cold-tolerant cultivar than intolerant cultivar. With the genome-wide characterization and phylogenetic analysis of Cucurbitaceae MAP65s, and the expression profile of CsaMAP65s in cucumber, this study laid a foundation for further study on MAP65 functions in developmental processes and responses to abiotic stress in Cucurbitaceae species.

Mutagenesis of the melon Prv gene by CRISPR/Cas9 breaks papaya ringspot virus resistance and generates an autoimmune allele with constitutive defense responses

The majority of plant disease resistance (R) genes encode nucleotide binding-leucine-rich repeat (NLR) proteins. In melon, two closely linked NLR genes, Fom-1 and Prv, were mapped and identified as candidate genes that control resistance to Fusarium oxysporum f.sp. melonis races 0 and 2, and to papaya ringspot virus (PRSV), respectively. In this study, we validated the function of Prv and showed that it is essential for providing resistance against PRSV infection. We generated CRISPR/Cas9 [clustered regularly interspaced palindromic repeats (CRISPR)/CRISPR-associated protein 9] mutants using Agrobacterium-mediated transformation of a PRSV-resistant melon genotype, and the T1 progeny proved susceptible to PRSV, showing strong disease symptoms and viral spread upon infection. Three alleles having 144, 154, and ~3 kb deletions, respectively, were obtained, all of which caused loss of resistance. Interestingly, one of the Prv mutant alleles, prvΔ154, encoding a truncated product, caused an extreme dwarf phenotype, accompanied by leaf lesions, high salicylic acid levels, and defense gene expression. The autoimmune phenotype observed at 25 °C proved to be temperature dependent, being suppressed at 32 °C. This is a first report on the successful application of CRISPR/Cas9 to confirm R gene function in melon. Such validation opens up new opportunities for molecular breeding of disease resistance in this important vegetable crop.

Genome-Wide Identification and Characterization Analysis of WUSCHEL-Related Homeobox Family in Melon (Cucumis melo L.)

WUSCHEL-related homeobox (WOX) proteins are very important in controlling plant development and stress responses. However, the WOX family members and their role in response to abiotic stresses are largely unknown in melon (Cucumis melo L.). In this study, 11 WOX (CmWOX) transcript factors with conserved WUS and homeobox motif were identified and characterized, and subdivided into modern clade, ancient clade and intermediate clade based on bioinformatic and phylogenetic analysis. Evolutionary analysis revealed that the CmWOX family showed protein variations in Arabidopsis, tomato, cucumber, melon and rice. Alignment of protein sequences uncovered that all CmWOXs had the typical homeodomain, which consisted of conserved amino acids. Cis-element analysis showed that CmWOX genes may response to abiotic stress. RNA-seq and qRT-PCR results further revealed that the expression of partially CmWOX genes are associated with cold and drought. CmWOX13a and CmWOX13b were constitutively expressed under abiotic stresses, CmWOX4 may play a role in abiotic processes during plant development. Taken together, this study offers new perspectives on the CmWOX family's interaction and provides the framework for research on the molecular functions of CmWOX genes.

Lessons learned from metabolic engineering in hairy roots: Transcriptome and metabolic profile changes caused by Rhizobium-mediated plant transformation in Cucurbitaceae species

Cucurbitaceae species are used in traditional medicine around the world. Cucurbitacins are highly oxygenated triterpenoids found in Cucurbitaceae species and exhibit potent anticancer activity alone and in combination with other existing chemotherapeutic drugs. Therefore, increasing production of these specialized metabolites is of great relevance. We recently showed that hairy roots of Cucurbita pepo can be used as a platform for metabolic engineering of cucurbitacins to modify their structure and increase their production. To study the changes in cucurbitacin accumulation upon formation of hairy roots, an empty vector (EV) control and Cucurbitacin inducing bHLH transcription factor 1 (CpCUCbH1)-overexpressing hairy roots of C. pepo were compared to untransformed (WT) roots. Whilst CpCUCbH1-overexpression increased production of cucurbitacins I and B by 5-fold, and cucurbitacin E by 3-fold when compared to EV lines, this increase was not significantly different when compared to WT roots. This indicated that Rhizobium rhizogenes transformation lowered the cucurbitacins levels in hairy roots, but that increasing expression of cucurbitacin biosynthetic genes by CpCUCbH1-overexpression restored cucurbitacin production to WT levels. Subsequent metabolomic and RNA-seq analysis indicated that the metabolic profile and transcriptome of hairy roots was significantly changed when compared to WT roots. Interestingly, it was observed that 11% of the differentially expressed genes were transcription factors. It was noteworthy that the majority of transcripts showing highest Pearson correlation coefficients to the Rhizobium rhizogenes genes rolB, rolC and ORF13a, were predicted to be transcription factors. In summary, hairy roots are an excellent platform for metabolic engineering of plant specialized metabolites, but these extensive transcriptome and metabolic profile changes should be considered in subsequent studies.

The transcription factor CmERFI-2 represses CmMYB44 expression to increase sucrose levels in oriental melon fruit

Soluble sugar accumulation in fruit ripening determines fleshy fruit quality. However, the molecular mechanism for this process is not yet understood. Here, we showed a transcriptional repressor, CmMYB44 regulates sucrose accumulation and ethylene synthesis in oriental melon (Cucumis. melo var. makuwa Makino) fruit. Overexpressing CmMYB44 suppressed sucrose accumulation and ethylene production. Furthermore, CmMYB44 repressed the transcriptional activation of CmSPS1 (sucrose phosphate synthase 1) and CmACO1 (ACC oxidase 1), two key genes in sucrose and ethylene accumulation, respectively. During the later stages of fruit ripening, the repressive effect of CmMYB44 on CmSPS1 and CmACO1 could be released by overexpressing CmERFI-2 (ethylene response factor I-2) and exogenous ethylene in "HS" fruit (high sucrose accumulation fruit). CmERFI-2 acted upstream of CmMYB44 as a repressor by directly binding the CmMYB44 promoter region, indirectly stimulating the expression level of CmSPS1 and CmACO1. Taken together, we provided a molecular regulatory pathway mediated by CmMYB44, which determines the degree of sucrose and ethylene accumulation in oriental melon fruit and sheds light on transcriptional responses triggered by ethylene sensing that enable the process of fruit ripening.

Selective deployment of virulence effectors correlates with host specificity in a fungal plant pathogen

The hemibiotrophic fungal plant pathogen Colletotrichum orbiculare is predicted to secrete hundreds of effector proteins when the pathogen infects cucurbit crops, such as cucumber and melon, and tobacco (Nicotiana benthamiana), a distantly related Solanaceae species. Here, we report the identification of sets of C. orbiculare effector genes that are differentially required for fungal virulence to two phylogenetically distant host species. Through targeted gene knockout screening of C. orbiculare 'core' effector candidates defined based on in planta gene expression, we identified: four host-specific virulence effectors (named effector proteins for cucurbit infection, or EPCs) that are required for full virulence of C. orbiculare to cucurbit hosts, but not to the Solanaceae host N. benthamiana; and five host-nonspecific virulence effectors, which collectively contribute to fungal virulence to both hosts. During host infection, only a small subset of genes, including the host-specific EPC effector genes, showed preferential expression on one of the hosts, while gene expression profiles of the majority of other genes, including the five host-nonspecific effector genes, were common to both hosts. This work suggests that C. orbiculare adopts a host-specific effector deployment strategy, in addition to general host-blind virulence mechanisms, for adaptation to cucurbit hosts.

PIF8-WRKY42-mediated salicylic acid synthesis modulates red light induced powdery mildew resistance in oriental melon

Light signals and plant hormones are involved in regulating the growth, development and stress resistance of plants; however, it remains unclear whether light affects hormones and thus pathogen resistance in oriental melon. Here, we found that red light promoted salicylic acid (SA) accumulation and powdery mildew resistance by activating the transcription of CmICS, the key gene for SA biosynthesis, and silencing CmICS seriously weakened the induction effect of red light on powdery mildew resistance in oriental melon leaves. Further studies showed that red light induced the expression of CmWRKY42 under powdery mildew stress, and CmWRKY42 directly bound to the CmICS promoter to activate its expression and promote the accumulation of SA under red light. Furthermore, we found that PHYTOCHROME INTERACTING FACTOR 8 (PIF8), as a negative regulator of SA biosynthesis, inhibits CmWRKY42 transcriptional activation by binding to the CmWRKY42 promoter, and thus inhibits transcriptional activation of CmICS by CmWRKY42. Also, CmPIF8 binds to the CmICS promoter and directly inhibits its transcription. In conclusion, our study revealed a new molecular mechanism of the relationship between red light-SA-powdery mildew resistance and provided a theoretical basis for resistance breeding of oriental melon.

QTL mapping of resistance to Pseudoperonospora cubensis clade 2, mating type A1, in Cucumis melo and dual-clade marker development

This is the first identification of QTLs underlying resistance in Cucumis melo to an isolate of Pseudoperonospora cubensis identified as Clade 2/mating type A1. Pseudoperonospora cubensis, causal organism of cucurbit downy mildew (CDM), causes severe necrosis and defoliation on Cucumis melo (melon). A recombinant inbred line population (N = 169) was screened against an isolate of P. cubensis (Clade 2/mating type A1) in replicated greenhouse and growth chamber experiments. SNPs (n = 5633 bins) identified in the RIL population were used for quantitative trait loci (QTL) mapping. A single major QTL on chromosome 10 (qPcub-10.3-10.4) was consistently associated with resistance across all experiments, while a second major QTL on chromosome 8 (qPcub-8.3) was identified only in greenhouse experiments. These two major QTLs were identified on the same chromosomes (8 and 10) but in different locations as two major QTLs (qPcub-8.2 and qPcub-10.1) previously identified for resistance to P. cubensis Clade 1/mating type A2. Kompetitive allele-specific PCR (KASP) markers were developed for these four major QTLs and validated in the RIL population through QTL mapping. These markers will provide melon breeders a high-throughput genotyping toolkit for development of melon cultivars with broad tolerance to CDM.

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.

Comparative bioinformatics analysis and abiotic stress responses of expansin proteins in Cucurbitaceae members: watermelon and melon

Watermelon and melon are members of the Cucurbitaceae family including economically significant crops in the world. The expansin protein family, which is one of the members of the cell wall, breaks down the non-covalent bonds between cell wall polysaccharides, causing pressure-dependent cell expansion. Comparative bioinformatics and molecular characterization analysis of the expansin protein family were carried out in the watermelon (Citrullus lanatus) and melon (Cucumis melo) plants in the study. Gene expression levels of expansin family members were analyzed in leaf and root tissues of watermelon and melon under ABA, drought, heat, cold, and salt stress conditions by quantitative real-time PCR analysis. After comprehensive searches, 40 expansin proteins (22 ClaEXPA, 14 ClaEXPLA, and 4 ClaEXPB) in watermelon and 43 expansin proteins (19 CmEXPA, 15 CmEXPLA, 3 CmEXPB, and 6 CmEXPLB) in melon were identified. The greatest orthologous genes were identified with soybean expansin genes for watermelon and melon. However, the latest divergence time between orthologous genes was determined with poplar expansin genes for watermelon and melon expansin genes. ClaEXPA-04, ClaEXPA-09, ClaEXPB-01, ClaEXPB-03, and ClaEXPLA-13 genes in watermelon and CmEXPA-12, CmEXPA-10, and CmEXPLA-01 genes in melon can be involved in tissue development and abiotic stress response of the plant. The current study combining bioinformatics and experimental analysis can provide a detailed characterization of the expansin superfamily which has roles in growth and reaction to the stress of the plant. The study ensures detailed data for future studies examining gene functions including the roles in plant growth and stress conditions.

Genome-Wide Identification of BES1 Gene Family in Six Cucurbitaceae Species and Its Expression Analysis in Cucurbita moschata

The BES1 (BRI1-EMSSUPPRESSOR1) gene family play a vital role in the BR (brassinosteroid) signaling pathway, which is involved in the growth and development, biotic, abiotic, and hormone stress response in many plants. However, there are few reports of BES1 in Cucurbita moschata. In this study, 50 BES1 genes were identified in six Cucurbitaceae species by genome-wide analysis, which could be classified into 3 groups according to their gene structural features and motif compositions, and 13 CmoBES1 genes in Cucurbita moschata were mapped on 10 chromosomes. Quantitative real-time PCR analysis showed that the CmoBES1 genes displayed differential expression under different abiotic stress and hormone treatments. Subcellular localization showed that the most of CmoBES1 proteins localized in nucleus and cytoplasm, and transactivation assay indicated 9 CmoBES1 proteins played roles as transcription factors. Our analysis of BES1s diversity, localization, and expression in Curcubitaceae contributes to the better understanding of the essential roles of these transcription factors in plants.

Genetic diversity analysis and variety identification using SSR and SNP markers in melon

Melon is an important horticultural crop with a pleasant aromatic flavor and abundance of health-promoting substances. Numerous melon varieties have been cultivated worldwide in recent years, but the high number of varieties and the high similarity between them poses a major challenge for variety evaluation, discrimination, as well as innovation in breeding. Recently, simple sequence repeats (SSRs) and single nucleotide polymorphisms (SNPs), two robust molecular markers, have been utilized as a rapid and reliable method for variety identification. To elucidate the genetic structure and diversity of melon varieties, we screened out 136 perfect SSRs and 164 perfect SNPs from the resequencing data of 149 accessions, including the most representative lines worldwide. This study established the DNA fingerprint of 259 widely-cultivated melon varieties in China using Target-seq technology. All melon varieties were classified into five subgruops, including ssp. agrestis, ssp. melo, muskmelon and two subgroups of foreign individuals. Compared with ssp. melo, the ssp. agrestis varieties might be exposed to a high risk of genetic erosion due to their extremely narrow genetic background. Increasing the gene exchange between ssp. melo and ssp. agrestis is therefore necessary in the breeding procedure. In addition, analysis of the DNA fingerprints of the 259 melon varieties showed a good linear correlation (R² = 0.9722) between the SSR genotyping and SNP genotyping methods in variety identification. The pedigree analysis based on the DNA fingerprint of 'Jingyu' and 'Jingmi' series melon varieties was consistent with their breeding history. Based on the SNP index analysis, ssp. agrestis had low gene exchange with ssp. melo in chromosome 4, 7, 10, 11and 12, two specific SNP loci were verified to distinguish ssp. agrestis and ssp. melon varieties. Finally, 23 SSRs and 40 SNPs were selected as the core sets of markers for application in variety identification, which could be efficiently applied to variety authentication, variety monitoring, as well as the protection of intellectual property rights in melon.

A recessive gene Cmpmr2F confers powdery mildew resistance in melon (Cucumis melo L.)

Identified a recessive gene (Cmpmr2F) associated with resistance to infection by the powdery mildew causing agent Podosphaera xanthii race 2F. Powdery mildew (PM) is one of the most destructive fungal diseases of melon, which significantly reduces the crop yield and quality. Multiple studies are being performed for in-depth genetic understandings of PM-susceptibility or -resistance mechanisms in melon plants, but the holistic knowledge of the precise genetic basis of PM-resistance is unexplored. In this study, we characterized the recessive gene "Cmpmr2F" and found its association with resistance against the PM causative agent "Podosphaera xanthii race 2F." Fine genetic mapping revealed the major-effect region of a 26.25-kb interval on chromosome 12, which harbored the Cmpmr2F gene corresponding to the MELO3C002403, encoding allantoate amidohydrolase. The functional gene annotation, expression pattern, and sequence alignment analyses were carried out using two contrast parent lines of melon "X055" PM-susceptible and "PI 124112" PM-resistant. Further, gene silencing of Cmpmr2F using virus-induced gene silencing (VIGS) significantly increased PM-resistance in the susceptible plant. In contrast to the previously reported studies, we identified that Cmpmr2F-silenced plants showed no impairment in growth due to less apparent negative effects in silenced melon plants. So, it is believed that the Cmpmr2F gene has great potential for further breeding studies to increase the P. xanthii race 2F resistance in melon. In short, our study provides new genetic resources and a solid foundation for further functional analysis of PM-resistance genes in melon, as well as powerful molecular markers for marker-assisted breeding aimed at developing new melon varieties resistant to PM infection.

Identification of Bottle Gourd (Lagenaria siceraria) OVATE Family Genes and Functional Characterization of LsOVATE1

The OVATE gene family is a class of conserved transcription factors that play significant roles in plant growth, development, and abiotic stress, and also affect fruit shape in vegetable crops. Bottle gourd (Lagenaria siceraria), commonly known as calabash or gourd, is an annual climber belonging to the Cucurbitaceae family. Studies on bottle gourd OVATE genes are limited. In this study, we performed genome-wide identification of the OVATE gene family in bottle gourd, and identified a total of 20 OVATE family genes. The identified genes were unevenly distributed across 11 bottle gourd chromosomes. We also analyzed the gene homology, amino acid sequence conservation, and three-dimensional protein structure (via prediction) of the 20 OVATE family genes. We used RNA-seq data to perform expression analysis, which found 20 OVATE family genes to be differentially expressed based on spatial and temporal characteristics, suggesting that they have varying functions in the growth and development of bottle gourd. In situ hybridization and subcellular localization analysis showed that the expression characteristics of the LsOVATE1 gene, located on chromosome 7 homologous to OVATE, is a candidate gene for affecting the fruit shape of bottle gourd. In addition, RT-qPCR data from bottle gourd roots, stems, leaves, and flowers showed different spatial expression of the LsOVATE1 gene. The ectopic expression of LsOVATE1 in tomato generated a phenotype with a distinct fruit shape and development. Transgenic-positive plants that overexpressed LsOVATE1 had cone-shaped fruit, calyx hypertrophy, petal degeneration, and petal retention after flowering. Our results indicate that LsOVATE1 could serve important roles in bottle gourd development and fruit shape determination, and provide a basis for future research into the function of LsOVATE1.

Pan-genome and multi-parental framework for high-resolution trait dissection in melon (Cucumis melo)

Linking genotype with phenotype is a fundamental goal in biology and requires robust data for both. Recent advances in plant-genome sequencing have expedited comparisons among multiple-related individuals. The abundance of structural genomic within-species variation that has been discovered indicates that a single reference genome cannot represent the complete sequence diversity of a species, leading to the expansion of the pan-genome concept. For high-resolution forward genetics, this unprecedented access to genomic variation should be paralleled and integrated with phenotypic characterization of genetic diversity. We developed a multi-parental framework for trait dissection in melon (Cucumis melo), leveraging a novel pan-genome constructed for this highly variable cucurbit crop. A core subset of 25 diverse founders (MelonCore25), consisting of 24 accessions from the two widely cultivated subspecies of C. melo, encompassing 12 horticultural groups, and 1 feral accession was sequenced using a combination of short- and long-read technologies, and their genomes were assembled de novo. The construction of this melon pan-genome exposed substantial variation in genome size and structure, including detection of ~300 000 structural variants and ~9 million SNPs. A half-diallel derived set of 300 F₂ populations, representing all possible MelonCore25 parental combinations, was constructed as a framework for trait dissection through integration with the pan-genome. We demonstrate the potential of this unified framework for genetic analysis of various melon traits, including rind color intensity and pattern, fruit sugar content, and resistance to fungal diseases. We anticipate that utilization of this integrated resource will enhance genetic dissection of important traits and accelerate melon breeding.

A common whole-genome paleotetraploidization in Cucurbitales

Cucurbitales are an important order of flowering plants known for encompassing edible plants of economic and medicinal value and numerous ornamental plants of horticultural value. By reanalyzing the genomes of two representative families (Cucurbitaceae and Begoniaceae) in Cucurbitales, we found that the previously identified Cucurbitaceae common paleotetraploidization that occurred shortly after the core-eudicot-common hexaploidization event is shared by Cucurbitales, including Begoniaceae. We built a multigenome alignment framework for Cucurbitales by identifying orthologs and paralogs and systematically redating key evolutionary events in Cucurbitales. Notably, characterizing the gene retention levels and genomic fractionation patterns between subgenomes generated from different polyploidizations in Cucurbitales suggested the autopolyploid nature of the Begoniaceae common tetraploidization and the allopolyploid nature of the Cucurbitales common tetraploidization and the Cucurbita-specific tetraploidization. Moreover, we constructed the ancestral Cucurbitales karyotype comprising 17 proto-chromosomes, confirming that the most recent common ancestor of Cucurbitaceae contained 15 proto-chromosomes and rejecting the previous hypothesis for an ancestral Cucurbitaceae karyotype with 12 proto-chromosomes. In addition, we found that the polyploidization and tandem duplication events promoted the expansion of gene families involved in the cucurbitacin biosynthesis pathway; however, gene loss and chromosomal rearrangements likely limited the expansion of these gene families.

Genome-wide association study and genetic mapping of BhWAX conferring mature fruit cuticular wax in wax gourd

BACKGROUND

Wax gourd [Benincasa hispida (Thunb) Cogn. (2n = 2x = 24)] is an economically important vegetable crop of genus Benincasa in the Cucurbitaceae family. Fruit is the main consumption organ of wax gourd. The mature fruit cuticular wax (MFCW) is an important trait in breeding programs, which is also of evolutionary significance in wax gourd. However, the genetic architecture of this valuable trait remains unrevealed.

RESULTS

In this study, genetic analysis revealed that the inheritance of MFCW was controlled by a single gene, with MFCW dominant over non-MFCW, and the gene was primarily named as BhWAX. Genome-wide association study (GWAS) highlighted a 1.1 Mb interval on chromosome 9 associated with MFCW in wax gourd germplasm resources. Traditional fine genetic mapping delimited BhWAX to a 0.5 Mb region containing 12 genes. Based on the gene annotation, expression analysis and co-segregation analysis, Bhi09G001428 that encodes a membrane bound O-acyltransferase (MBOAT) was proposed as the candidate gene for BhWAX. Moreover, it was demonstrated that the efficiency of a cleaved amplified polymorphic sequences (CAPS) marker in the determination of MFCW in wax gourd reached 80%.

CONCLUSIONS

In closing, the study identified the candidate gene controlling MFCW and provided an efficient molecular marker for the trait in wax gourd for the first time, which will not only be beneficial for functional validation of the gene and marker-assisted breeding of wax gourd, but also lay a foundation for analysis of its evolutionary meaning among cucurbits.

Genome-Wide Identification of the SAUR Gene Family in Wax Gourd (Benincasa hispida) and Functional Characterization of BhSAUR60 during Fruit Development

The wax gourd (Benincasa hispida) is an important vegetable crop whose fruits contain nutrients and metabolites. Small auxin upregulated RNA (SAUR) genes constitute the largest early auxin-responsive gene family and regulate various biological processes in plants, although this gene family has not been studied in the wax gourd. Here, we performed genome-wide identification of the SAUR gene family in wax gourds and analyzed their syntenic and phylogenetic relationships, gene structures, conserved motifs, cis-acting elements, and expression patterns. A total of 68 SAUR (BhSAUR) genes were identified, which were distributed on nine chromosomes with 41 genes in two clusters. More than half of the BhSAUR genes were derived from tandem duplication events. The BhSAUR proteins were classified into seven subfamilies. BhSAUR gene promoters contained cis-acting elements involved in plant hormone and environmental signal responses. Further expression profiles showed that BhSAUR genes displayed different expression patterns. BhSAUR60 was highly expressed in fruits, and overexpression led to longer fruits in Arabidopsis. In addition, the plants with overexpression displayed longer floral organs and wavy stems. In conclusion, our results provide a systematic analysis of the wax gourd SAUR gene family and facilitate the functional study of BhSAUR60 during wax gourd fruit development.

Editing melon eIF4E associates with virus resistance and male sterility

The cap-binding protein eIF4E, through its interaction with eIF4G, constitutes the core of the eIF4F complex, which plays a key role in the circularization of mRNAs and their subsequent cap-dependent translation. In addition to its fundamental role in mRNA translation initiation, other functions have been described or suggested for eIF4E, including acting as a proviral factor and participating in sexual development. We used CRISPR/Cas9 genome editing to generate melon eif4e knockout mutant lines. Editing worked efficiently in melon, as we obtained transformed plants with a single-nucleotide deletion in homozygosis in the first eIF4E exon already in a T0 generation. Edited and non-transgenic plants of a segregating F2 generation were inoculated with Moroccan watermelon mosaic virus (MWMV); homozygous mutant plants showed virus resistance, while heterozygous and non-mutant plants were infected, in agreement with our previous results with plants silenced in eIF4E. Interestingly, all homozygous edited plants of the T0 and F2 generations showed a male sterility phenotype, while crossing with wild-type plants restored fertility, displaying a perfect correlation between the segregation of the male sterility phenotype and the segregation of the eif4e mutation. Morphological comparative analysis of melon male flowers along consecutive developmental stages showed postmeiotic abnormal development for both microsporocytes and tapetum, with clear differences in the timing of tapetum degradation in the mutant versus wild-type. An RNA-Seq analysis identified critical genes in pollen development that were down-regulated in flowers of eif4e/eif4e plants, and suggested that eIF4E-specific mRNA translation initiation is a limiting factor for male gametes formation in melon.

Karyotyping of commercial cultivars of melon (Cucumis melo L.)

BACKGROUND

This study on cultivars of melon (Cucumis melo L.) marketed in Brazil was conducted to obtain information to be used in breeding programs of this species. Little is known about the karyotype variability among C. melo L. cultivars targeted at the consumer market. The objective of the present study was to verify the karyotype variability in eight commercial melon cultivars used in the Brazilian market.

METHODS AND RESULTS

Slides were stained with 2% Giemsa and assembled with Neomount to perform chromosomal morphometry. GC-rich heterochromatin was observed by CMA₃/DAPI staining. 5 S rDNA, centromeric satellite DNA (SatDNA), and telomeric sites were visualized using fluorescence in situ hybridization. All images were captured on an Olympus BX41 microscope equipped with a 5 M Olympus DP25 digital camera and DP2-BSW software. The cultivars showed symmetrical karyotypes with significant differences in total chromosome length and average chromosome size. Heterochromatic CMA₃⁺ blocks were observed in terminal regions related to satellites (secondary constrictions), as well as in centromeric and pericentromeric regions. A single chromosomal pair of 5 S rDNA sites was observed in all cultivars, but at distinct locations. Centromeric satellite sequences, tested for the first time in melon, revealed only centromeric sites. Telomeric sites were observed in all the chromosomes of the cultivars.

CONCLUSIONS

Karyotype variation was observed in cultivars of melon, which were analyzed for chromosomal morphology and localization of GC-rich heterochromatin, as well centromeric SatDNA, rDNA, and telomeric chromosomal markers. Hence, these cultivars can be used in future breeding programs.

Plant Metabolic Engineering by Multigene Stacking: Synthesis of Diverse Mogrosides

Mogrosides are a group of health-promoting natural products that extracted from Siraitia grosvenorii fruit (Luo-han-guo or monk fruit), which exhibited a promising practical application in natural sweeteners and pharmaceutical development. However, the production of mogrosides is inadequate to meet the need worldwide, and uneconomical synthetic chemistry methods are not generally recommended for structural complexity. To address this issue, an in-fusion based gene stacking strategy (IGS) for multigene stacking has been developed to assemble 6 mogrosides synthase genes in pCAMBIA1300. Metabolic engineering of Nicotiana benthamiana and Arabidopsis thaliana to produce mogrosides from 2,3-oxidosqualene was carried out. Moreover, a validated HPLC-MS/MS method was used for the quantitative analysis of mogrosides in transgenic plants. Herein, engineered Arabidopsis thaliana produced siamenoside I ranging from 29.65 to 1036.96 ng/g FW, and the content of mogroside III at 202.75 ng/g FW, respectively. The production of mogroside III was from 148.30 to 252.73 ng/g FW, and mogroside II-E with concentration between 339.27 and 5663.55 ng/g FW in the engineered tobacco, respectively. This study provides information potentially applicable to develop a powerful and green toolkit for the production of mogrosides.

Melon shoot organization 1, encoding an AGRONAUTE7 protein, plays a crucial role in plant development

A melon gene MSO1 located on chromosome 10 by map-based cloning strategy, which encodes an ARGONAUTE 7 protein, is responsible for the development of shoot organization. Plant endogenous small RNAs (sRNAs) are involved in various plant developmental processes. In Arabidopsis, sRNAs combined with ARGONAUTE (AGO) proteins are incorporated into the RNA-induced silencing complex (RISC), which functions in RNA silencing or biogenesis of trans-acting siRNAs (ta-siRNAs). However, their roles in melon (Cucumis melo L.) are still unclear. Here, the melon shoot organization 1 (mso1) mutant was identified and shown to exhibit pleiotropic phenotypes in leaf morphology and plant architecture. Positional cloning of MSO1 revealed that it encodes a homologue of Arabidopsis AGO7/ZIPPY, which is required for the production of ta-siRNAs. The AG-to-C mutation in the second exon of MSO1 caused a frameshift mutation and significantly reduced its expression. Ectopic expression of MSO1 rescued the Arabidopsis ago7 phenotype. RNA-seq analysis showed that several genes involved in transcriptional regulation and plant hormones were significantly altered in mso1 compared to WT. A total of 304 and 231 miRNAs were identified in mso1 and WT by sRNA sequencing, respectively, and among them, 42 known and ten novel miRNAs were differentially expressed. cme-miR390a significantly accumulated, and the expression levels of the two ta-siRNAs were almost completely abolished in mso1. Correspondingly, their targets, the ARF3 and ARF4 genes, showed dramatically upregulated expression, indicating that the miR390-TAS3-ARF pathway has conserved roles in melon. These findings will help us better understand the molecular mechanisms of MSO1 in plant development in melon.

Natural variation in the NAC transcription factor NONRIPENING contributes to melon fruit ripening

The NAC transcription factor NONRIPENING (NOR) is a master regulator of climacteric fruit ripening. Melon (Cucumis melo L.) has climacteric and non-climacteric fruit ripening varieties and is an ideal model to study fruit ripening. Two natural CmNAC-NOR variants, the climacteric haplotype CmNAC-NOR^S,N and the non-climacteric haplotype CmNAC-NOR^A,S , have effects on fruit ripening; however, their regulatory mechanisms have not been elucidated. Here, we report that a natural mutation in the transcriptional activation domain of CmNAC-NOR^S,N contributes to climacteric melon fruit ripening. CmNAC-NOR knockout in the climacteric-type melon cultivar "BYJH" completely inhibited fruit ripening, while ripening was delayed by 5-8 d in heterozygous cmnac-nor mutant fruits. CmNAC-NOR directly activated carotenoid, ethylene, and abscisic acid biosynthetic genes to promote fruit coloration and ripening. Furthermore, CmNAC-NOR mediated the transcription of the "CmNAC-NOR-CmNAC73-CmCWINV2" module to enhance flesh sweetness. The transcriptional activation activity of the climacteric haplotype CmNAC-NOR^S,N on these target genes was significantly higher than that of the non-climacteric haplotype CmNAC-NOR^A,S . Moreover, CmNAC-NOR^S,N complementation fully rescued the non-ripening phenotype of the tomato (Solanum lycopersicum) cr-nor mutant, while CmNAC-NOR^A,S did not. Our results provide insight into the molecular mechanism of climacteric and non-climacteric fruit ripening in melon.

Regulation of ethylene biosynthesis and signal transduction by nitric oxide leading to resistance against Alternaria alternata in Hami melon

BACKGROUND

Hami melons are tasty and nutritive, but susceptibility to the fungus Alternaria alternata is one of the main problems leading to the postharvest loss of this fruit. The purpose of this research was to evaluate the effectiveness of nitric oxide (NO) on regulation of ethylene biosynthesis as well as signal transduction against black spot disease caused by A. alternata in the Hami melon.

RESULTS

Nitric oxide reduced the growth of lesion diameter and lesion depth in melons inoculated with A. alternata. Ethylene production was significantly inhibited by NO, which was supported by the reduction of 1-aminocyclopropene-1-carboxylate (ACC) synthase (ACS) activity and the deferment of ACC content and ACC oxidase (ACO) activity. Nitric oxide treatment also significantly regulated the expression of four ethylene biosynthesis genes CmACS1, CmACS2, CmACO1, and CmACO2, and eight signal ethylene transduction genes CmETR1, CmETR2, CmCTR1, CmEIN2, CmEIL1, CmEBF1, CmERF1B and CmERF2. The modes of NO regulating these genes can be divided into five categories: promotion (CmEIN2, and CmEIL1), delay (CmACS1, CmETR2, CmCTR1 and CmERF2), up-regulation (CmETR1, CmEBF1 and CmERF1B), down-regulation (CmACS2), and first inhibition and then induction (CmACO1 and CmACO2).

CONCLUSION

The NO treatment enhanced the postharvest disease resistance of Hami melon attacked by A. alternata, possibly by postponing ethylene biosynthesis and signal transduction. © 2021 Society of Chemical Industry.

Spanish Melon Landraces: Revealing Useful Diversity by Genomic, Morphological, and Metabolomic Analysis

Spain is a secondary centre of the diversification of the melon (Cucumis melo L.), with high diversity represented in highly appreciated landraces belonging to the Flexuosus and Ibericus groups. A collection of 47 accessions of Flexuosus, Chate, Piel de Sapo, Tendral, Amarillo, Blanco, and Rochet was analysed using a genotyping-by-sequencing (GBS) approach. A total of 66,971 quality SNPs were identified. Genetic analysis differentiated Ibericus accessions and exotic materials (Ameri, Momordica, Kachri, and Agrestis), while Flexuous accessions shared ancestry between them. Within the Ibericus group, no clear genomic distinction could be identified for the different landraces evaluated, with accessions of different landraces showing high genetic similarity. The morphological characterization confirmed that the external colour and fruit shape had been used as recognition patterns for Spanish melon landraces, but variability within a landrace exists. Differences were found in the sugars and acid and volatile profiles of the materials. Flexuosus and Chate melons at the immature commercial stage accumulated malic acid and low levels of hexoses, while Ibericus melons accumulated high contents of sucrose and citric acid. Specific trends could be identified in the Ibericus landraces. Tendral accumulated low levels of sugars and citric acid and high of malic acid, maintaining higher firmness, Rochet reached higher levels of sugars, and Amarillo tended to lower malic acid contents. Interestingly, high variability was found within landraces for the acidic profile, offering possibilities to alter taste tinges. The main volatile organic compounds (VOCs) in Flexuosus and Chate were aldehydes and alcohols, with clear differences between both groups. In the Ibericus landraces, general trends for VOC accumulation could be identified, but, again, a high level of variation exists. This situation highlights the necessity to develop depuration programs to promote on-farm in situ conservation and, at the same time, offers opportunities to establish new breeding program targets and to take advantage of these sources of variation.

A single substitution in Vacuolar protein sorting 4 is responsible for resistance to Watermelon mosaic virus in melon

In plants, introgression of genetic resistance is a proven strategy for developing new resistant lines. While host proteins involved in genome replication and cell to cell movement are widely studied, other cell mechanisms responsible for virus infection remain under investigated. Endosomal sorting complexes required for transport (ESCRT) play a key role in membrane trafficking in plants and are involved in the replication of several plant RNA viruses. In this work, we describe the role of the ESCRT protein CmVPS4 as a new susceptibility factor to the Potyvirus Watermelon mosaic virus (WMV) in melon. Using a worldwide collection of melons, we identified three different alleles carrying non-synonymous substitutions in CmVps4. Two of these alleles were shown to be associated with WMV resistance. Using a complementation approach, we demonstrated that resistance is due to a single non-synonymous substitution in the allele CmVps4P30R. This work opens up new avenues of research on a new family of host factors required for virus infection and new targets for resistance.

CRISPR/Cas9 gene editing uncovers the roles of CONSTITUTIVE TRIPLE RESPONSE 1 and REPRESSOR OF SILENCING 1 in melon fruit ripening and epigenetic regulation

Melon (Cucumis melo) has emerged as an alternative model to tomato for studying fruit ripening due to the coexistence of climacteric and non-climacteric varieties. Previous characterization of a major quantitative trait locus (QTL), ETHQV8.1, that is able to trigger climacteric ripening in a non-climacteric background resulted in the identification of a negative regulator of ripening CTR1-like (MELO3C024518) and a putative DNA demethylase ROS1 (MELO3C024516) that is the orthologue of DML2, a DNA demethylase that regulates fruit ripening in tomato. To understand the role of these genes in climacteric ripening, in this study we generated homozygous CRISPR knockout mutants of CTR1-like and ROS1 in a climacteric genetic background. The climacteric behavior was altered in both loss-of-function mutants in two growing seasons with an earlier ethylene production profile being observed compared to the climacteric wild type, suggesting a role of both genes in climacteric ripening in melon. Single-cytosine methylome analyses of the ROS1-knockout mutant revealed changes in DNA methylation in the promoter regions of the key ripening genes such as ACS1, ETR1, and ACO1, and in transcription factors associated with ripening including NAC-NOR, RIN, and CNR, suggesting the importance of ROS1-mediated DNA demethylation for triggering fruit ripening in melon.

Comparative Transcriptome Analysis Identified Key Pathways and Genes Regulating Differentiated Stigma Color in Melon (Cucumis melo L.)

Stigma color is an important morphological trait in many flowering plants. Visual observations in different field experiments have shown that a green stigma in melons is more attractive to natural pollinators than a yellow one. In the current study, we evaluated the characterization of two contrasted melon lines (MR-1 with a green stigma and M4-7 with a yellow stigma). Endogenous quantification showed that the chlorophyll and carotenoid content in the MR-1 stigmas was higher compared to the M4-7 stigmas. The primary differences in the chloroplast ultrastructure at different developmental stages depicted that the stigmas of both melon lines were mainly enriched with granum, plastoglobulus, and starch grains. Further, comparative transcriptomic analysis was performed to identify the candidate pathways and genes regulating melon stigma color during key developmental stages (S1–S3). The obtained results indicated similar biological processes involved in the three stages, but major differences were observed in light reactions and chloroplast pathways. The weighted gene co-expression network analysis (WGCNA) of differentially expressed genes (DEGs) uncovered a “black” network module (655 out of 5302 genes), mainly corresponding to light reactions, light harvesting, the chlorophyll metabolic process, and the chlorophyll biosynthetic process, and exhibited a significant contribution to stigma color. Overall, the expression of five key genes of the chlorophyll synthesis pathway—CAO (MELO03C010624), CHLH (MELO03C007233), CRD (MELO03C026802), HEMA (MELO03C011113), POR (MELO03C016714)—were checked at different stages of stigma development in both melon lines using quantitative real time polymerase chain reaction (qRT-PCR). The results exhibited that the expression of these genes gradually increased during the stigma development of the MR-1 line but decreased in the M4-7 line at S2. In addition, the expression trends in different stages were the same as RNA-seq, indicating data accuracy. To sum up, our research reveals an in-depth molecular mechanism of stigma coloration and suggests that chlorophyll and related biological activity play an important role in differentiating melon stigma color.

Transcriptome and Metabolome Analysis Provide New Insights into the Process of Tuberization of Sechium edule Roots

Chayote (Sechium edule) produces edible tubers with high starch content after 1 year of growth but the mechanism of chayote tuberization remains unknown. ‘Tuershao’, a chayote cultivar lacking edible fruits but showing higher tuber yield than traditional chayote cultivars, was used to study tuber formation through integrative analysis of the metabolome and transcriptome profiles at three tuber-growth stages. Starch biosynthesis- and galactose metabolism-related genes and metabolites were significantly upregulated during tuber bulking, whereas genes encoding sugars will eventually be exported transporter (SWEET) and sugar transporter (SUT) were highly expressed during tuber formation. Auxin precursor (indole-3-acetamide) and ethylene precursor, 1-aminocyclopropane-1-carboxylic acid, were upregulated, suggesting that both hormones play pivotal roles in tuber development and maturation. Our data revealed a similar tuber-formation signaling pathway in chayote as in potatoes, including complexes BEL1/KNOX and SP6A/14-3-3/FDL. Down-regulation of the BEL1/KNOX complex and upregulation of 14-3-3 protein implied that these two complexes might have distinct functions in tuber formation. Finally, gene expression and microscopic analysis indicated active cell division during the initial stages of tuber formation. Altogether, the integration of transcriptome and metabolome analyses unraveled an overall molecular network of chayote tuberization that might facilitate its utilization.

Specificity of Resistance and Tolerance to Cucumber Vein Yellowing Virus in Melon Accessions and Resistance Breaking with a Single Mutation in VPg

Cucumber vein yellowing virus (CVYV) is an emerging virus on cucurbits in the Mediterranean Basin, against which few resistance sources are available, particularly in melon. The melon accession PI 164323 displays complete resistance to isolate CVYV-Esp, and accession HSD 2458 presents a tolerance, i.e., very mild symptoms despite virus accumulation in inoculated plants. The resistance is controlled by a dominant allele Cvy-1¹, while the tolerance is controlled by a recessive allele cvy-2, independent from Cvy-1¹. Before introducing the resistance or tolerance in commercial cultivars through a long breeding process, it is important to estimate their specificity and durability. Upon inoculation with eight molecularly diverse CVYV isolates, the resistance was found to be isolate-specific because many CVYV isolates induced necrosis on PI 164323, whereas the tolerance presented a broader range. A resistance-breaking isolate inducing severe mosaic on PI 164323 was obtained. This isolate differed from the parental strain by a single amino acid change in the VPg coding region. An infectious CVYV cDNA clone was obtained, and the effect of the mutation in the VPg cistron on resistance to PI 164323 was confirmed by reverse genetics. This represents the first determinant for resistance-breaking in an ipomovirus. Our results indicate that the use of the Cvy-1¹ allele alone will not provide durable resistance to CVYV and that, if used in the field, it should be combined with other control methods such as cultural practices and pyramiding of resistance genes to achieve long-lasting resistance against CVYV.

Comparative transcriptomic analysis of two Cucumis melo var. saccharinus germplasms differing in fruit physical and chemical characteristics

BACKGROUND

Hami melon (Cucumis melo var. saccharinus) is a popular fruit in China because of its excellent taste, which is largely determined by its physicochemical characteristics, including flesh texture, sugar content, aroma, and nutrient composition. However, the mechanisms by which these characteristics are regulated have not yet been determined. In this study, we monitored changes in the fruits of two germplasms that differed in physicochemical characteristics throughout the fruit development period.

RESULTS

Ripe fruit of the bred variety 'Guimi' had significantly higher soluble sugar contents than the fruit of the common variety 'Yaolong.' Additionally, differences in fruit shape and color between these two germplasms were observed during development. Comparative transcriptome analysis, conducted to identify regulators and pathways underlying the observed differences at corresponding stages of development, revealed a higher number of differentially expressed genes (DEGs) in Guimi than in Yaolong. Moreover, most DEGs detected during early fruit development in Guimi were associated with cell wall biogenesis. Temporal analysis of the identified DEGs revealed similar trends in the enrichment of downregulated genes in both germplasms, although there were differences in the enrichment trends of upregulated genes. Further analyses revealed trends in differential changes in multiple genes involved in cell wall biogenesis and sugar metabolism during fruit ripening.

CONCLUSIONS

We identified several genes associated with the ripening of Hami melons, which will provide novel insights into the molecular mechanisms underlying the development of fruit characteristics in these melons.