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.

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.

Time-Series Transcriptome of Cucumis melo Reveals Extensive Transcriptomic Differences with Different Maturity

As the most important melon cultivar grown in the north-western provinces of China, Hami melon (Cucumis melo) produces large edible fruits that serve as an important dietary component in the world. In general, as a climacteric plant, melon harvested at 60% maturity results in a product with bad quality, while the highest-quality product can be guaranteed when harvesting at 90% maturity. In order to clarify the genetic basis of their distinct profiles of metabolite accumulation, we performed systematic transcriptome analyses between 60% and 90% maturity melons. A total of 36 samples were sequenced and over 1.7 billion reads were generated. Differentially expressed genes in 60% and 90% maturity melons were detected. Hundreds of these genes were functionally enriched in the sucrose and citric acid accumulation process of C. melo. We also detected a number of distinct splicing events between 60% and 90% maturity melons. Many genes associated with sucrose and citric acid accumulation displayed as differentially expressed or differentially spliced between different degrees of maturity of Hami melons, including CmCIN2, CmSPS2, CmBGAL3, and CmSPS2. These results demonstrate that the phenotype pattern differences between 60% and 90% maturity melons may be largely resulted from the significant transcriptome regulation.

CmCML11 interacts with CmCAMTA5 to enhance γ-aminobutyric acid (GABA) accumulation by regulating GABA shunt in fresh-cut cantaloupe

The effect of calcium chloride (CaCl2) treatment on γ-aminobutyric acid (GABA) accumulation in fresh-cut cantaloupe and the involved mechanisms were investigated. The result showed that 1% (w/v) CaCl2 treatment increased GABA content and activities of glutamate decarboxylase (GAD) and succinate semialdehyde dehydrogenase (SSADH), while decreased glutamate (Glu) content and GABA transaminase (GABA-T) activities in fresh-cut cantaloupe. CmCML11 and CmCAMTA5 expressions of CaCl2-treated fruit increased by 187.4% and 165.6% than control fruit in the initial 6 h. Besides, expressions of GABA shunt genes, including CmGAD1, CmGAD2, CmGABA-T and CmSSADH were also up-regulated by CaCl2 treatment during early storage. Moreover, acting as a transcriptional activator, CmCAMTA5 could bind to the CG-box in promoters of CmGAD1, CmGABA-T and CmSSADH and activate their transcription. Furthermore, the interaction between CmCML11 and CmCAMTA5 could enhance the transcriptional activation on GABA shunt genes which were regulated by CmCAMTA5. Collectively, our findings revealed that CaCl2 treatment promoted GABA accumulation in fresh-cut cantaloupe via the combined effect of CmCML11 and CmCAMTA5 in the regulation of expressions of CmGAD1, CmGABA-T, and CmSSADH in GABA shunt.

Genome-wide identification of the melon (Cucumis melo L.) response regulator gene family and functional analysis of CmRR6 and CmPRR3 in response to cold stress

The response regulator (RR) gene family play crucial roles in cytokinin signal transduction, plant development, and resistance to abiotic stress. However, there are no reports on the identification and functional characterization of RR genes in melon. In this study, a total of 18 CmRRs were identified and classified into type A, type B, and clock PRRs, based on phylogenetic analysis. Most of the CmRRs displayed tissue-specific expression patterns, and some were induced by cold stress according to two RNA-seq datasets. The expression patterns of CmRR2/6/11/15 and CmPRR2/3 under cold treatment were confirmed by qRT-PCR. Subcellular localization assays indicated that CmRR6 and CmPRR3 were primarily localized in the nucleus and chloroplast. Furthermore, when either CmRR6 or CmPRR3 were silenced using tobacco ringspot virus (TRSV), the cold tolerance of the virus-induced gene silencing (VIGS) melon plants were significantly enhanced, as evidenced by measurements of chlorophyll fluorescence, ion leakage, reactive oxygen, proline, and malondialdehyde levels. Additionally, the expression levels of CmCBF1, CmCBF2, and CmCBF3 were significantly increased in CmRR6-silenced and CmPRR3-silenced plants under cold treatment. Our findings suggest that CmRRs contribute to cold stress responses and provide new insights for further pursuing the molecular mechanisms underlying CmRRs-mediated cold tolerance in melon.

Genome-wide identification of the OMT gene family in Cucumis melo L. and expression analysis under abiotic and biotic stress

Background

O-methyltransferase (OMT)-mediated O-methylation is a frequent modification that occurs during natural product biosynthesis, and it increases the diversity and stability of secondary metabolites. However, detailed genome-wide identification and expression analyses of OMT gene family members have not been performed in melons. In this study, we aimed to perform the genome-wide identification of OMT gene family members in melon to identify and clarify their actions during stress.

Methods

Genome-wide identification of OMT gene family members was performed using data from the melon genome database. The Cucumis melo OMT genes (CmOMTs) were then compared with the genes from two representative monocotyledons and three representative dicotyledons. The basic information, cis-regulatory elements in the promoter, predicted 3-D-structures, and GO enrichment results of the 21 CmOMTs were analyzed.

Results

In our study, 21 CmOMTs (named CmOMT1-21) were obtained by analyzing the melon genome. These genes were located on six chromosomes and divided into three groups composed of nine, six, and six CmOMTs based on phylogenetic analysis. Gene structure and motif descriptions were similar within the same classes. Each CmOMT gene contains at least one cis-acting element associated with hormone transport regulation. Analysis of cis-acting elements illustrated the potential role of CmOMTs in developmental regulation and adaptations to various abiotic and biotic stresses. The RNA-seq and quantitative real-time PCR (qRT-PCR) results indicated that NaCl stress significantly induced CmOMT6/9/14/18 and chilling and high temperature and humidity (HTH) stresses significantly upregulated CmOMT14/18. Furthermore, the expression pattern of CmOMT18 may be associated with Fusarium oxysporum f. sp. melonis race 1.2 (FOM1.2) and powdery mildew resistance. Our study tentatively explored the biological functions of CmOMT genes in various stress regulation pathways and provided a conceptual basis for further detailed studies of the molecular mechanisms.

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.

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.

Genome-wide identification and analysis of Lateral Organ Boundaries Domain (LBD) transcription factor gene family in melon (Cucumis melo L.)

Background

Lateral Organ Boundaries Domain (LBD) transcription factor (TF) gene family members play very critical roles in several biological processes like plant-spesific development and growth process, tissue regeneration, different biotic and abiotic stress responses in plant tissues and organs. The LBD genes have been analyzed in various species. Melon (Cucumis melo L.), a member of the Cucurbitaceae family, is economically important and contains important molecules for nutrition and human health such as vitamins A and C, β-carotenes, phenolic acids, phenolic acids, minerals and folic acid. However, no studies have been reported so far about LBD genes in melon hence this is the first study for LBD genes in this plant.

Results

In this study, 40 melon CmLBD TF genes were identified, which were separated into seven groups through phylogenetic analysis. Cis-acting elements showed that these genes were associated with plant growth and development, phytohormone and abiotic stress responses. Gene Ontology (GO) analysis revealed that of CmLBD genes especially function in regulation and developmental processes. The in silico and qRT-PCR expression patterns demonstrated that CmLBD01 and CmLBD18 are highly expressed in root and leaf tissues, CmLBD03 and CmLBD14 displayed a high expression in male-female flower and ovary tissues.

Conclusions

These results may provide important contributions for future research on the functional characterization of the melon LBD gene family and the outputs of this study can provide information about the evolution and characteristics of melon LBD gene family for next studies.

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.

Natural allelic variation in the EamA-like transporter, CmSN, is associated with fruit skin netting in melon

KEY MESSAGE

A SNP mutation in CmSN, encoding an EamA-like transporter, is responsible for fruit skin netting in melon. In maturing melon (Cucumis melo L.), the rind becomes reticulated or netted, a unique characteristic that dramatically changes the appearance of the fruit. However, little is known about the molecular basis of fruit skin netting formation in this important cucurbit crop. Here, we conducted map-based cloning of a skin netting (CmSN) locus using segregating populations derived from the cross between the smooth-fruit line H906 and the netted-fruit line H581. The results showed that CmSN was controlled by a single dominant gene and was primarily positioned on melon chromosome 2, within a physical interval of ~ 351 kb. Further fine mapping in a large F2 population narrowed this region to a 71-kb region harboring 5 genes. MELO3C010288, which encodes a protein in the EamA-like transporter family, is the best possible candidate gene for the netted phenotype. Two nonsynonymous single nucleotide polymorphisms (SNPs) were identified in the third and sixth exons of the CmSN gene and co-segregated with the skin netting (SN) phenotype among the genetic population. A genome-wide association study (GWAS) determined that CmSN is probably a domestication gene under selective pressure during the subspecies C. melo subsp. melo differentiation. The SNP in the third exon of CmSN (the leading SNP in GWAS) revealed a bi-allelic diversity in natural accessions with SN traits. Our results lay a foundation for deciphering the molecular mechanism underlying the formation of fruit skin netting in melon, as well as provide a strategy for genetic improvement of netted fruit using a marker-assisted selection approach.

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.

Physiological analysis and transcriptome profiling reveals the impact of microplastic on melon (Cucumis melo L.) seed germination and seedling growth

The wide application of agricultural plastics leads to microplastic (MP) accumulation in the soil and inevitably result in MP pollution. Melon is an economically important horticultural crop that is widely cultivated with plastic film mulching. However, the impact of MP pollution on plant growth remains largely unclear. Here we reported the morphological, physiological, biochemical responses and transcriptome re-programing of melon responses to MP on seed germination and seedling growth. Polyvinyl chloride particles were added to potting mix to simulate MP exposure environment (MEE). The results showed that low and medium concentrations (1-4 g kg-1) of MEE had a significant adverse effect on seed germination and seedling growth. In both cases, the germination potential was decreased, young root forks increased, and tips decreased; and the dry weight of seedlings, the total length, surface area, forks and tips of root were also decreased. However, the root activity was increased. The concentration of MEE to give the best parameters was at 2 g kg-1. Catalase enzymatic activity and reactive oxygen species (ROS) in roots were decreased continuously with increased MEE concentrations. The peak values of peroxidase activity, O2.- content and generation rate, ROS enrichment and malondialdehyde content all reached the highest at 2 g kg-1. MEE also increased the proline content and decreased the contents of ascorbic acid, soluble sugar and soluble protein in these seedlings. Medium and high concentrations of MEE (4-8 g kg-1) also increased the chlorophyll b content. Low concentrations MEE (1-2 g kg-1) inhibited actual photochemical efficiency of photosystem II and photochemical quenching, two key chlorophyll fluorescence parameters. Transcriptome analysis showed that the differentially expressed genes caused by the MEE were mainly belonged to defense response, signal transduction, hormone metabolism, plant-pathogen interaction, and phenylpropanoid biosynthesis. The results of this study will help to understand the ecotoxicological effects of MEE on melons and provide data for ecological risk assessment of Cucurbitaceae vegetable cultivation.

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.

CmMYB44 might interact with CmAPS2-2 to regulate starch metabolism in oriental melon fruit

Sugar content is one of the determining factors for melon fruit maturity. Studies have shown that starch gradually degrades during fruit ripening, resulting in sugar accumulation. But the specific relationship between starch metabolism and sucrose accumulation was still unknown. Here, the starch and sugar contents, the activities of key enzymes and the expression patterns of genes related to starch-sucrose metabolism were determined in the fruit of high sugar and starch variety 'HS' and low sugar and starch variety 'LW'. It was found that starch accumulated during fruit development process, and then degraded at 30 days after anthesis (DAA), which was synchronized with sucrose accumulation in 'HS' fruit, while starch and sucrose contents were always at a lower level during 'LW' fruit maturation. Furthermore, starch metabolism-related enzymes (Adenine dinucleotide phosphate -glucose pyrophosphorylase (AGPase), α-amylase (AMY), β-amylase (BMY)) and the key enzymes for sucrose accumulation (sucrose phosphate synthase (SPS) and sucrose synthase (SS)) were significantly increased at ripening stage of 'HS' fruit, and their activities were consistent with the expressions of CmAPS2-2, CmAMY2, CmBAM1, CmBAM9 and CmSPS1. However, the contents of starch and sucrose and the activities of AGPase and SPS in 'LW' fruit didn't change significantly. We discovered an R2R3-type MYB transcription factor, CmMYB44, screened from yeast one hybrid library, could directly bind to the promoter of CmAPS2-2 to inhibit its transcription. These results revealed that the targeted down-regulation of CmAPS2-2 by CmMYB44 might be involved in the starch accumulation process, which affect the flavor quality of oriental melon fruit.

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.

Novel Bisexual Flower Control Gene Regulates Sex Differentiation in Melon (Cucumis melo L.)

The sex-control system involves several mechanisms in melon. The present study identified a novel bisexual flower control gene from the hermaphroditic melon germplasm, different from the previously recognized one. Genetic analysis showed that a single recessive gene in the newly identified locus b controlled the bisexual flower phenotype in melons. We generated 1431 F₂ segregating individuals for genetic mapping of locus b, which was delimited to a 47.94 kb region. Six candidate genes were identified in the delimited interval, and candidate No. 4 encoding melon CPR5 protein was selected as the suitable one for locus b and was denoted CmCPR5. CPR5 reportedly interacted with ethylene receptor ETR1 to regulate ethylene signal transduction. Moreover, the ethephon assays showed that the parental lines (unisexual line and bisexual line) had contrasting expression patterns of CmCPR5. The BiFC and LCI assays also confirmed that CmCPR5 interacted with CmETR1 in 0426 but not in Y101. However, crossover tests showed that CmETR1 functioned normally in both parental lines, suggesting CPR5 malfunction in Y101. This study proposed a corollary mechanism of bisexual flower regulation during stamen primordium development in which the inhibition of stamen primordia development was prevented by the malfunctioning CmCPR5, resulting in bisexual flowers.

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.

Physiological and Transcriptomic Analyses Uncover the Reason for the Inhibition of Photosynthesis by Phosphate Deficiency in Cucumis melo L

Phosphate (Pi) deficiency is a common phenomenon in agricultural production and limits plant growth. Recent work showed that long-term Pi deficiency caused the inhibition of photosynthesis and inefficient electron transport. However, the underlying mechanisms are still unknown. In this study, we used the physiological, histochemical, and transcriptomic methods to investigate the effect of low-Pi stress on photosynthetic gas exchange parameters, cell membrane lipid, chloroplast ultrastructure, and transcriptional regulation of key genes in melon seedlings. The results showed that Pi deficiency significantly downregulated the expression of aquaporin genes, induced an increase in ABA levels, and reduced the water content and free water content of melon leaves, which caused physiological drought in melon leaves. Therefore, gas exchange was disturbed. Pi deficiency also reduced the phospholipid contents in leaf cell membranes, caused the peroxidation of membrane lipids, and destroyed the ultrastructure of chloroplasts. The transcriptomic analysis showed that 822 differentially expressed genes (DEGs) were upregulated and 1254 downregulated by Pi deficiency in leaves. GO and KEGG enrichment analysis showed that DEGs significantly enriched in chloroplast thylakoid membrane composition (GO:0009535), photosynthesis-antenna proteins (map00196), and photosynthesis pathways (map00195) were downregulated by Pi deficiency. It indicated that Pi deficiency regulated photosynthesis-related genes at the transcriptional level, thereby affecting the histochemical properties and physiological functions, and consequently causing the reduced light assimilation ability and photosynthesis efficiency. It enriches the mechanism of photosynthesis inhibition by Pi deficiency.

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.

Nano-selenium enhances the antioxidant capacity, organic acids and cucurbitacin B in melon (Cucumis melo L.) plants

Pesticides are widely used in melon production causing safety issues around the consumption of melon and increasing pathogen and insect tolerance to pesticides. This study investigated whether a nano-selenium (Nano-Se) spray treatment can improve resistance to biological stress in melon plants, reducing the need for pesticides, and how this mechanism is activated. To achieve this, we examine the ultrastructure and physio-biochemical responses of two melon cultivars after foliar spraying with Nano-Se. Nano-Se treatment reduced plastoglobulins in leaf mesophyll cells, thylakoid films were left intact, and compound starch granules increased. Nano-Se treatment also increased root mitochondria and left nucleoli intact. Nano-Se treatment enhanced ascorbate peroxidase, peroxidase, phenylalanine ammonia lyase, β-1,3-glucanase, chitinase activities and their mRNA levels in treated melon plants compared to control plants (without Nano-Se treatments). Exogenous application of Nano-Se improved fructose, glucose, galactitol, stachyose, lactic acid, tartaric acid, fumaric acid, malic acid and succinic acid in treated plants compared to control plants. In addition, Nano-Se treatment enhanced cucurbitacin B and up-regulated eight cucurbitacin B synthesis-related genes. We conclude that Nano-Se treatment of melon plants triggered antioxidant capacity, photosynthesis, organic acids, and up-regulated cucurbitacin B synthesis-related genes, which plays a comprehensive role in stress resistance in melon plants.

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.

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.

Biparental genetic mapping reveals that CmCLAVATA3 (CmCLV3) is responsible for the variation in carpel number in melon (Cucumis melo L.)

Genetic analysis revealed that CmCLV3 is a candidate gene for the variation in melon carpel number. Carpel number (CN) is an important trait in melon. Three-CN melon fruit is oval, while 5-CN melon fruit has a round or flat shape. Herein, a genetic analysis of a population in which the CN locus was segregated indicated that 3-CN is controlled by a major dominant effective gene. Bulked segregant analysis and initial linkage mapping placed the CN locus in a 6.67 Mb region on chromosome 12, and it was narrowed to 882.19 kb with molecular markers and recombinant plants. Fine mapping with a large F₂ population containing 1026 individuals further narrowed the locus to an 83.98 kb region harboring five annotated genes. Gene structure alignment between the parental lines revealed MELO3C035640.2 (annotated as CLAVATA3, CmCLV3) as the best candidate gene for the CN trait. CmCLV3 was more highly expressed in 3- than 5-CN lines and specifically expressed in terminal buds rather than in young leaves, hypocotyls, and roots. The CmCLV3 coding region was cloned from eight 3- or 5-CN melon accessions, and a nonsynonymous SNP site was highly correlated with CN variation. This SNP site was also related to CN variations among 40 melon lines according to their resequencing data, causing a helix alteration in the CmCLV3 protein. Promoter region sequence alignment and activity analysis showed that, unlike in cucumber and tomato, CmCLV3 promoter variation and activity were not the main reasons for CN alteration. Overall, this study provides a genetic resource for melon fruit development research and molecular breeding tools for melon CN improvement.

Genome-wide characterization and expression analysis of SAUR gene family in Melon (Cucumis melo L.)

We identified 66 melon SAUR genes by bioinformatic analyses. CmSAUR19, 38, 58, 62 genes are specifically expressed in different stages of fruit growth, suggesting their participation in regulating fruit development. Auxin plays a crucial role in plant growth by regulating the multiple auxin response genes. However, in melon (Cucumis melo L.), the functions of the auxin early response gene family SAUR (Small auxin up RNA) genes in fruit development are still poorly understood. Through genome-wide characterization of CmSAUR family in melon, we identified a total of 66 CmSAUR genes. The open reading frames of the CmSAUR genes ranged from 234 to 525 bp, containing only one exon and lacking introns. Chromosomal position and phylogenetic tree analyses found that the two gene clusters in the melon chromosome are highly homologous in the Cucurbitaceae plants. Among the four conserved motifs in CmSAUR proteins, motif 1, motif 2, and motif 3 located in most of the family protein sequences, and motif 4 showed a close correlation with the two gene clusters. The CmSAUR28 and CmSAUR58 genes have auxin response elements located in the promoters, suggesting they may be involved in the auxin signaling pathway to regulate fruit development. Through transcriptomic profiling in the four developmental stages of fruit and different lateral organs, we selected 16 differentially-expressed SAUR genes for performing further expression analyses. qRT-PCR results showed that five SAUR genes are specifically expressed in flower organs and ovaries. CmSAUR19 and CmSAUR58 were significantly accumulated in the early developmental stage of the fruit. CmSAUR38 and CmAUR62 showed high expression in the climacteric and post-climacteric stages, suggesting their specific role in controlling fruit ripening. This work provides a foundation for further exploring the function of the SAUR gene in fruit development.

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.

CmPMRl and CmPMrs are responsible for resistance to powdery mildew caused by Podosphaera xanthii race 1 in Melon

Two genes for resistance to Podosphaera xanthii race 1 in melon were identified on chromosomes 10 and 12 of the Cucumis melo cultivar MR-1. Cucumis melo L. is an economically important crop, the production of which is threatened by the prevalence of melon powdery mildew (PM) infections. We herein utilized the MR-1 (P₁; resistant to PM) and M4-7 (P₂; susceptible to PM) accessions to assess the heritability of PM (race 1) resistance in these melon plants. PM resistance in MR-1 leaves was linked to a dominant gene (CmPMRl), whereas stem resistance was under the control of a recessive gene (CmPMrs), with the dominant gene having an epistatic effect on the recessive gene. The CmPMRl gene was mapped to a 50 Kb interval on chromosome 12, while CmPMrs was mapped to an 89 Kb interval on chromosome 10. The CmPMRl candidate gene MELO3C002441 and the CmPMrs candidate gene MELO3C012438 were identified through sequence alignment, functional annotation, and expression pattern analyzes of all genes within these respective intervals. MELO3C002441 and MELO3C012438 were both localized to the cellular membrane and were contained conserved NPR gene-like and MLO domains, respectively, which were linked to PM resistance. In summary, we identified patterns of PM resistance in the disease-resistant MR-1 melon cultivar and identified two putative genes linked to resistance. Our results offer new genetic resources and markers to guide future marker-assisted breeding for PM resistance in melon.

CmFSI8/CmOFP13 encoding an OVATE family protein controls fruit shape in melon

Fruit shape is an important quality and yield trait in melon (Cucumis melo). Although some quantitative trait loci for fruit shape have been reported in in this species, the genes responsible and the underlying mechanisms remain poorly understood. Here, we identified and characterized a gene controlling fruit shape from two melon inbred lines, B8 with long-horn fruit and HP22 with flat-round fruit. Genetic analysis suggested that the shape was controlled by a single and incompletely dominant locus, which we designate as CmFSI8/CmOFP13. This gene was finely mapped to a 53.7-kb interval on chromosome 8 based on bulked-segregant analysis sequencing and map-based cloning strategies. CmFSI8/CmOFP13 encodes an OVATE family protein (OFP) and is orthologous to AtOFP1 and SlOFP20. The transcription level of CmFSI8/CmOFP13 in the ovary of HP22 was significantly higher than that in B8, and sequence analysis showed that a 12.5-kb genomic variation with a retrotransposon insertion identified in the promoter was responsible for elevating the expression, and this ultimately caused the differences in fruit shape. Ectopic overexpression of CmFSI8/CmOFP13 in Arabidopsis led to multiple phenotypic changes, including kidney-shaped leaves and shortened siliques. Taken together, our results demonstrate the involvement of an OFP in regulating fruit shape in melon, and our improved understanding of the molecular mechanisms will enable us to better manipulate fruit shape in breeding.

Fine-mapping and identification of a candidate gene controlling seed coat color in melon (Cucumis melo L. var. chinensis Pangalo)

MELO3C019554 encoding a homeobox protein (PHD transcription factor) is a candidate gene that involved in the formation of seed coat color in melon. Seed coat color is related to flavonoid content which is closely related to seed dormancy. According to the genetic analysis of a six-generation population derived from two parents (IC2508 with a yellow seed coat and IC2518 with a brown seed coat), we discovered that the yellow seed coat trait in melon is controlled by a single dominant gene, named CmBS-1. Bulked segregant analysis sequencing (BSA-Seq) revealed that the gene is located at 11,860,000-15,890,000 bp (4.03 Mb) on Chr 6. The F₂ population was genotyped using insertion-deletions (InDels), from which cleaved amplified polymorphic sequence (dCAPS) markers were derived to construct a genetic map. The gene was then fine-mapped to a 233.98 kb region containing 12 genes. Based on gene sequence analysis with two parents, we found that the MELO3C019554 gene encoding a homeobox protein (PHD transcription factor) had a nonsynonymous single nucleotide polymorphism (SNP) mutation in the coding sequence (CDS), and the SNP mutation resulted in the conversion of an amino acid (A → T) at residue 534. In addition, MELO3C019554 exhibited lower relative expression levels in the yellow seed coat than in the brown seed coat. Furthermore, we found that MELO3C019554 is related to 12 flavonoid metabolites. Thus, we predicted that MELO3C019554 is a candidate gene controlling seed coat color in melon. The study lays a foundation for further cloning projects and functional analysis of this gene, as well as marker-assisted selection breeding.

CmLHP1 proteins play a key role in plant development and sex determination in melon (Cucumis melo)

In monoecious melon (Cucumis melo), sex is determined by the differential expression of sex determination genes (SDGs) and adoption of sex-specific transcriptional programs. Histone modifications such as H3K27me3 have been previously shown to be a hallmark associated to unisexual flower development in melon; yet, no genetic approaches have been conducted for elucidating the roles of H3K27me3 writers, readers, and erasers in this process. Here we show that melon homologs to Arabidopsis LHP1, CmLHP1A and B, redundantly control several aspects of plant development, including sex expression. Cmlhp1ab double mutants displayed an overall loss and redistribution of H3K27me3, leading to a deregulation of genes involved in hormone responses, plant architecture, and flower development. Consequently, double mutants display pleiotropic phenotypes and, interestingly, a general increase of the male:female ratio. We associated this phenomenon with a general deregulation of some hormonal response genes and a local activation of male-promoting SDGs and MADS-box transcription factors. Altogether, these results reveal a novel function for CmLHP1 proteins in maintenance of monoecy and provide novel insights into the polycomb-mediated epigenomic regulation of sex lability in plants.

A novel insight into transcriptional and epigenetic regulation underlying sex expression and flower development in melon (Cucumis melo L.)

Melon (Cucumis melo L.) is an important cucurbit and has been considered as a model plant for studying sex determination. The four most common sexual morphotypes in melon are monoecious (A-G-M), gynoecious (--ggM-), andromonoecious (A-G-mm), and hermaphrodite (--ggmm). Sex expression in melons is complex, as the genes and associated networks that govern the sex expression are not fully explored. Recently, RNA-seq transcriptomic profiling, ChIP-qPCR analysis integrated with gene ontology annotation and Kyoto Encyclopedia of Genes and Genomes pathways predicted the differentially expressed genes including sex-specific ACS and ACO genes, in regulating the sex-expression, phytohormonal cross-talk, signal transduction, and secondary metabolism in melons. Integration of transcriptional control through genetic interaction in between the ACS7, ACS11, and WIP1 in epistatic or hypostatic manner, along with the recruitment of H3K9ac and H3K27me3, epigenetically, overall determine sex expression. Alignment of protein sequences for establishing phylogenetic evolution, motif comparison, and protein-protein interaction supported the structural conservation while presence of the conserved hydrophilic and charged residues across the diverged evolutionary group predicted the functional conservation of the ACS protein. Presence of the putative cis-binding elements or DNA motifs, and its further comparison with DAP-seq-based cistrome and epicistrome of Arabidopsis, unraveled strong ancestry of melons with Arabidopsis. Motif comparison analysis also characterized putative genes and transcription factors involved in ethylene biosynthesis, signal transduction, and hormonal cross-talk related to sex expression. Overall, we have comprehensively reviewed research findings for a deeper insight into transcriptional and epigenetic regulation of sex expression and flower development in melons.

Reconstruction of ancestral karyotype illuminates chromosome evolution in the genus Cucumis

Karyotype dynamics driven by complex chromosome rearrangements constitute a fundamental issue in evolutionary genetics. The evolutionary events underlying karyotype diversity within plant genera, however, have rarely been reconstructed from a computed ancestral progenitor. Here, we developed a method to rapidly and accurately represent extant karyotypes with the genus, Cucumis, using highly customizable comparative oligo-painting (COP) allowing visualization of fine-scale genome structures of eight Cucumis species from both African-origin and Asian-origin clades. Based on COP data, an evolutionary framework containing a genus-level ancestral karyotype was reconstructed, allowing elucidation of the evolutionary events that account for the origin of these diverse genomes within Cucumis. Our results characterize the cryptic rearrangement hotspots on ancestral chromosomes, and demonstrate that the ancestral Cucumis karyotype (n = 12) evolved to extant Cucumis genomes by hybridizations and frequent lineage- and species-specific genome reshuffling. Relative to the African species, the Asian species, including melon (Cucumis melo, n = 12), Cucumis hystrix (n = 12) and cucumber (Cucumis sativus, n = 7), had highly shuffled genomes caused by large-scale inversions, centromere repositioning and chromothripsis-like rearrangement. The deduced reconstructed ancestral karyotype for the genus allowed us to propose evolutionary trajectories and specific events underlying the origin of these Cucumis species. Our findings highlight that the partitioned evolutionary plasticity of Cucumis karyotype is primarily located in the centromere-proximal regions marked by rearrangement hotspots, which can potentially serve as a reservoir for chromosome evolution due to their fragility.

QTL mapping of resistance to Pseudoperonospora cubensis clade 1, mating type A2, in Cucumis melo

This is the first identification of QTLs underlying resistance to Pseudoperonospora cubensis in Cucumis melo using a genetically characterized isolate. Pseudoperonospora cubensis, causal organism of cucurbit downy mildew (CDM), is one of the largest threats to cucurbit production in the eastern USA. Currently, no Cucumis melo (melon) cultivars have significant levels of resistance. Additionally, little is understood about the genetic basis of resistance in C. melo. Recombinant inbred lines (RILs; N = 169) generated from a cross between the resistant melon breeding line MR-1 and susceptible cultivar Ananas Yok'neam were phenotyped for CDM resistance in both greenhouse and growth chamber studies. A high-density genetic linkage map with 5,663 binned SNPs created from the RIL population was utilized for QTL mapping. Nine QTLs, including two major QTLs, were associated with CDM resistance. Of the major QTLs, qPcub-10.1 was stable across growth chamber and greenhouse tests, whereas qPcub-8.2 was detected only in growth chamber tests. qPcub-10.1 co-located with an MLO-like protein coding gene, which has been shown to confer resistance to powdery mildew and Phytophthora in other plants. This is the first screening of C. melo germplasm with a genetically characterized P. cubensis isolate.

High-quality chromosome-level genomes of Cucumis metuliferus and Cucumis melo provide insight into Cucumis genome evolution

Cucumis metuliferus (African horned cucumber), a wild relative of Cucumis sativus (cucumber) and Cucumis melo (melon), displays high-level resistance to several important plant pathogens (e.g., root-knot nematodes and several viruses). Here, we report a chromosome-level genome assembly for C. metuliferus, with a 316 Mb genome sequence comprising 29 039 genes. Phylogenetic analysis of related species in family Cucurbitaceae indicated that the divergence time between C. metuliferus and melon was 17.8 million years ago. Comparisons between the C. metuliferus and melon genomes revealed large structural variations (inversions and translocations >1 Mb) in eight chromosomes of these two species. Gene family comparison showed that C. metuliferus has the largest number of resistance-related nucleotide-binding site leucine-rich repeat (NBS-LRR) genes in Cucurbitaceae. The loss of NBS-LRR loci caused by large insertions or deletions (indels) and pseudogenization caused by small indels explained the loss of NBS-LRR genes in Cucurbitaceae. Population structure analysis suggested that C. metuliferus originated in Zimbabwe, then spread to other southern African regions where it likely underwent similar domestic selection as melon. This C. metuliferus reference sequence will accelerate the understanding of the molecular evolution of resistance-related genes and enhance cucurbit crop improvement efforts.

Genome-Wide SNP Markers for Genotypic and Phenotypic Differentiation of Melon (Cucumis melo L.) Varieties Using Genotyping-by-Sequencing

Melon (Cucumis melo L.) is an economically important horticultural crop with abundant morphological and genetic variability. Complex genetic variations exist even among melon varieties and remain unclear to date. Therefore, unraveling the genetic variability among the three different melon varieties, muskmelon (C. melo subsp. melo), makuwa (C. melo L. var. makuwa), and cantaloupes (C. melo subsp. melo var. cantalupensis), could provide a basis for evolutionary research. In this study, we attempted a systematic approach with genotyping-by-sequencing (GBS)-derived single nucleotide polymorphisms (SNPs) to reveal the genetic structure and diversity, haplotype differences, and marker-based varieties differentiation. A total of 6406 GBS-derived SNPs were selected for the diversity analysis, in which the muskmelon varieties showed higher heterozygote SNPs. Linkage disequilibrium (LD) decay varied significantly among the three melon varieties, in which more rapid LD decay was observed in muskmelon (r² = 0.25) varieties. The Bayesian phylogenetic tree provided the intraspecific relationships among the three melon varieties that formed, as expected, individual clusters exhibiting the greatest genetic distance based on the posterior probability. The haplotype analysis also supported the phylogeny result by generating three major networks for 48 haplotypes. Further investigation for varieties discrimination allowed us to detect a total of 52 SNP markers that discriminated muskmelon from makuwa varieties, of which two SNPs were converted into cleaved amplified polymorphic sequence markers for practical use. In addition to these markers, the genome-wide association study identified two SNPs located in the genes on chromosome 6, which were significantly associated with the phenotypic traits of melon seed. This study demonstrated that a systematic approach using GBS-derived SNPs could serve to efficiently classify and manage the melon varieties in the genebank.

QTLs and candidate genes analyses for fruit size under domestication and differentiation in melon (Cucumis melo L.) based on high resolution maps

BACKGROUND

Melon is a very important horticultural crop produced worldwide with high phenotypic diversity. Fruit size is among the most important domestication and differentiation traits in melon. The molecular mechanisms of fruit size in melon are largely unknown.

RESULTS

Two high-density genetic maps were constructed by whole-genome resequencing with two F₂ segregating populations (WAP and MAP) derived from two crosses (cultivated agrestis × wild agrestis and cultivated melo × cultivated agrestis). We obtained 1,871,671 and 1,976,589 high quality SNPs that show differences between parents in WAP and MAP. A total of 5138 and 5839 recombination events generated 954 bins in WAP and 1027 bins in MAP with the average size of 321.3 Kb and 301.4 Kb respectively. All bins were mapped onto 12 linkage groups in WAP and MAP. The total lengths of two linkage maps were 904.4 cM (WAP) and 874.5 cM (MAP), covering 86.6% and 87.4% of the melon genome. Two loci for fruit size were identified on chromosome 11 in WAP and chromosome 5 in MAP, respectively. An auxin response factor and a YABBY transcription factor were inferred to be the candidate genes for both loci.

CONCLUSION

The high-resolution genetic maps and QTLs analyses for fruit size described here will provide a better understanding the genetic basis of domestication and differentiation, and provide a valuable tool for map-based cloning and molecular marker assisted breeding.

Drought-induced ABA, H2O2 and JA positively regulate CmCAD genes and lignin synthesis in melon stems

BACKGROUND

Cinnamyl alcohol dehydrogenase (CAD) is an important enzyme functions at the last step in lignin monomer synthesis pathway. Our previous work found that drought induced the expressions of CmCAD genes and promoted lignin biosynthesis in melon stems.

RESULTS

Here we studied the effects of abscisic acid (ABA), hydrogen peroxide (H2O2) and jasmonic acid (JA) to CmCADs under drought stress. Results discovered that drought-induced ABA, H2O2 and MeJA were prevented efficiently from increasing in melon stems pretreated with fluridone (Flu, ABA inhibitor), imidazole (Imi, H2O2 scavenger) and ibuprofen (Ibu, JA inhibitor). ABA and H2O2 are involved in the positive regulations to CmCAD1, 2, 3, and 5, and JA is involved in the positive regulations to CmCAD2, 3, and 5. According to the expression profiles of lignin biosynthesis genes, ABA, H2O2 and MeJA all showed positive regulations to CmPAL2-like, CmPOD1-like, CmPOD2-like and CmLAC4-like. In addition, positive regulations were also observed with ABA to CmPAL1-like, CmC4H and CmCOMT, with H2O2 to CmPAL1-like, CmC4H, CmCCR and CmLAC17-like, and with JA to CmCCR, CmCOMT, CmLAC11-like and CmLAC17-like. As expected, the signal molecules positively regulated CAD activity and lignin biosynthesis under drought stress. Promoter::GUS assays not only further confirmed the regulations of the signal molecules to CmCAD1~3, but also revealed the important role of CmCAD3 in lignin synthesis due to the strongest staining of CmCAD3 promoter::GUS.

CONCLUSIONS

CmCADs but CmCAD4 are positively regulated by ABA, H2O2 and JA under drought stress and participate in lignin synthesis.

Genetic Mapping and Nucleotide Diversity of Two Powdery Mildew Resistance Loci in Melon (Cucumis melo)

Powdery mildew (PM) significantly and negatively affects the yield and quality of melon (Cucumis melo) worldwide. Race 2F is the predominant physiological race of the pathogen Podosphaera xanthii in many regions. We used accessions PMR 6 (P₁; resistant to PM) and M1-7 (P₂; susceptible to PM) to analyze the inheritance of resistance to PM (race 2F). The ratio between resistant and susceptible individuals fits a Mendelian segregation ratio of 13:3 in a total of 256 F₂ individuals and 1:1 in BC₁P₂. The resistance to PM in PMR 6 was governed by two genes: a dominant (AA) gene with an epistatic effect and a recessive gene (bb). Only individuals with aaBB or aaBb genotypes were susceptible to PM. Two PM resistance loci, Pm2.1 and pm12.1, were mapped on chromosomes 2 and 12 by bulked segregant analysis and secondary mapping by quantitative trait loci analysis with 18 markers. A new marker-assisted selection system to identify melon genotypes resistant or susceptible to PM was developed and tested in 93 melon accessions. Nucleotide diversity (π) and fixation index (Fst) for the two PM resistance loci were estimated using resequencing data of 336 melons from three groups: C. melo subsp. agrestis, Cucumis melo subsp. melo, and the intermediate type. The lowest π was observed in C. melo ssp. agrestis, and the highest Fst value was between C. melo ssp. agrestis and C. melo ssp. melo. The findings provide a promising tool that can be used to accelerate breeding for durable resistance to PM.

Transcriptome profiling reveals the occurrence mechanism of bisexual flowers in melon (Cucumis melo L.)

Most cultivated melons are andromonoecies in which male flowers arose both in main stem and lateral branches but bisexual flowers only emerged from the leaf axils of lateral branches. However, bisexual flowers emerged in leaf axils of main stem after ethephon treatment. Therefore, the mechanism regulating the occurrence of bisexual flowers were investigated by performing transcriptome analysis in two comparison sets: shoot apex of main stem (MA) versus that of lateral branches (LA), and shoot apex of main stem after ethephon treatment (Eth) versus control (Cont). KEGG results showed that genes involved in "plant hormone signal transduction", "MAPK signaling pathway" and "carbon metabolism" were significantly upregulated both in LA and Eth. Further, details of DEGs involved in ethylene signaling pathway were surveyed and six genes were co-upregulated in two comparison sets. Among these, CmERF1, downstream in ethylene signaling pathway, showed the most significantly difference and expressed higher in bisexual buds than that in male buds. Furthermore, fifteen DEGs were found to contain GCC box or CRT/DRE cis-element for CmERF1 in their putative promoter region, and these DEGs involved in several plant hormones signaling pathway, camalexin synthesis, carbon metabolism and plant pathogen interaction.

Genome-wide identification, characterization, and expression analysis related to autotoxicity of the GST gene family in Cucumis melo L

Glutathione S-transferase (GST) plays an important role in plant resistance to biotic and abiotic stresses. In this paper, the characteristics of melon GST gene family members were analyzed from a genome-wide perspective. Forty-nine GSTs were identified in melon genome, belonging to eight classes. Through the phylogenetic analysis of GST proteins in melon and other plants, it was found that members from the same subfamily in different species clustered together, indicating that the subfamilies of GST have diversified before the divergence within these species. The results of chromosome mapping showed that GSTs were present in all chromosomes except for chromosome 5. Gene replication events played an important role in the expansion and evolution of melon GST gene family. Ten GSTs with significant differential expression were screened in the transcriptome database related to melon autotoxicity stress. The differential expression of these 10 GSTs was detected in roots and leaves of melon seedlings treated with cinnamic acid. The relative expression level of CmGSTU7, CmGSTU10, CmGSTU18, CmGSTF2 and CmGSTL1 in roots of melon seedlings was significantly higher than that in control group. It suggested that the five GSTs might play an important role in cinnamic acid mediated autotoxicity stress in melon. The results of this paper were helpful to reveal the evolution and functional succession of GST family and further understand the response of GST to autotoxicity stress in melon.

A calmodulin-like CmCML13 from Cucumis melo improved transgenic Arabidopsis salt tolerance through reduced shoot's Na+, and also improved drought resistance

The calmodulin-like proteins (CMLs) are a large family involved in plant biological processes. A calmodulin-like gene CmCML13 (GenBank accession number: MT340534) from melon (Cucumis melo L.) was isolated and functionally analyzed. CmCML13 was predicted to possess 3 EF-hands in which only the first EF-hand could bind with Ca²⁺. Subcellular localization assay revealed that CmCML13 was localized in nucleus, cell membrane, vacuolar membrane and cytoplasmic strand. The transcript level of CmCML13 was temporally and spatially regulated under salt stress. Constitutive expression of CmCML13 in the Arabidopsis thaliana enhanced salt tolerance at seeds germination. CmCML13 improved the transgenic Arabidopsis plants salt tolerance by significantly reducing Na⁺ content of shoots, which was unrelated to HKT1-involving pathway. Moreover, overexpressing of CmCML13 in Arabidopsis showed stronger drought tolerance. This study demonstrates that the CmCML13 is an important multifunctional protein associated with salt and drought stress, which may play a key role in stress signaling pathway.

Specific response mechanism to autotoxicity in melon (Cucumis melo L.) root revealed by physiological analyses combined with transcriptome profiling

Melon is of great value in food, medicine and industry. In recent years, the continuous cropping obstacles of melon is increasingly prominent, which seriously affects the cultivation. Autotoxicity is the key factor for the obstacles. Root is the first line against autotoxicity and main organs for autotoxins secretion. Some physiological responses and differentially expressed genes (DEGs) related to autotoxicity are only limited to root system. Considering the lack of relevant research, physiological researches combined with transcriptome sequencing of melon seedling after autotoxicity stress mediated by root exudates (RE) was performed to help characterize the response mechanism to autotoxicity in melon roots. The results showed that autotoxicity inhibited root morphogenesis of melon seedlings, induced the excessive accumulation of reactive oxygen species (ROS) and lipid peroxidation in roots, and activated most antioxidant enzymes. Compared with the control group, the osmoregulation substance content was always at a high level. DEGs response to autotoxicity in roots were distinguished from that in leaves. Functional annotation of these DEGs suggested that autotoxicity affected biological regulation in a negative manner. DEGs were mainly involved in the synthesis of antioxidants, DNA damage and metabolism, and stress response. These setbacks were associated with the deterioration of root morphogenesis, generation of dwarf and slender roots, and ultimately leading to plant death. The results may provide important information for revealing the response mechanism of root to autotoxicity, and provide theoretical basis for solving the continuous cropping obstacles in melon.

Supportive role of the Na+ transporter CmHKT1;1 from Cucumis melo in transgenic Arabidopsis salt tolerance through improved K+/Na+ balance

KEY MESSAGE

CmHKT1;1 selectively exports Na+ from plant cells. Upon NaCl stress, its expression increased in a salt-tolerant melon cultivar. Overexpression of CmHKT1;1 increased transgenic Arabidopsis salt tolerance through improved K+/Na+ balance. High-affinity K+ transporters (HKTs) are thought to be involved in reducing Na+ in plant shoots under salt stress and modulating salt tolerance, but their function in a moderately salt-tolerant species of melon (Cucumis melo L.) remains unclear. In this study, a Na+ transporter gene, CmHKT1;1 (GenBank accession number: MK986658), was isolated from melons based on genome data. The transcript of CmHKT1;1 was relatively more abundant in roots than in stems or leaves from melon seedlings. The tobacco transient expression system showed that CmHKT1;1 was plasma-membrane localized. Upon salt stress, CmHKT1;1 expression was more strongly upregulated in a salt-tolerant melon cultivar, 'Bingxuecui' (BXC) compared with a salt-sensitive cultivar, 'Yulu' (YL). Electrophysiological evidence demonstrated that CmHKT1;1 only transported Na+, rather than K+, when expressed in Xenopus laevis oocytes. Overexpression of CmHKT1;1 increased salt sensitivity in Saccharomyces cerevisiae and salt tolerance in Arabidopsis thaliana. Under NaCl treatments, transgenic Arabidopsis plants accumulated significantly lower concentrations of Na+ in shoots than wild type plants and showed a better K+/Na+ balance, leading to better Fv/Fm, root length, biomass, and enhanced plant growth. The CmHKT1;1 gene may serve as a useful candidate for improving crop salt tolerance.

Genome-wide characterization of two-component system (TCS) genes in melon (Cucumis melo L.)

To better understand cytokinin signaling in melon (Cucumis melo L.), one of the most important fruit crops in the Cucurbitaceae family, we identified and characterized melon two-component system (TCS) genes in this study. The results showed that there were 51 genes encoding putative TCS proteins in melon, and these TCS genes were classified into 3 subgroups, with 17 HK(L)s (histidine kinase/histidine-kinase like; 9 HKs and 8 HKLs), 9 HPs (histidine phosphotransfer proteins; 6 authentic and 3 pseudo), and 25 RRs (response regulators; 8 Type-A, 11 Type-B and 6 pseudo). The identity values of these cytokinin signaling proteins were revealed by analyzing their conserved motifs, domains and amino acid sequences. By analyzing TCS genes in different plant species, we found that melon HK(L)s, HPs and RRs had closer phylogenetic relationships with cucumber genes than with the genes of other plants, and the expansion of melon cytokinin signaling genes might be attributed to segmental duplication events. Analysis of the putative promoter regions (2-kb upstream regions of the start codon) revealed the enrichment of stress- and hormone-response cis-elements. The involvement of these putative TCS genes in melon cytokinin signaling was further supported by qRT-PCR data.

High-density NGS-based map construction and genetic dissection of fruit shape and rind netting in Cucumis melo

Melon is an important crop that exhibits broad variation for fruit morphology traits that are the substrate for genetic mapping efforts. In the post-genomic era, the link between genetic maps and physical genome assemblies is key for leveraging QTL mapping results for gene cloning and breeding purposes. Here, using a population of 164 melon recombinant inbred lines (RILs) that were subjected to genotyping-by-sequencing, we constructed and compared high-density sequence- and linkage-based recombination maps that were aligned to the reference melon genome. These analyses reveal the genome-wide variation in recombination frequency and highlight regions of disrupted collinearity between our population and the reference genome. The population was phenotyped over 3 years for fruit size and shape as well as rind netting. Four QTLs were detected for fruit size, and they act in an additive manner, while significant epistatic interaction was found between two neutral loci for this trait. Fruit shape displayed transgressive segregation that was explained by the action of four QTLs, contributed by alleles from both parents. The complexity of rind netting was demonstrated on a collection of 177 diverse accessions. Further dissection of netting in our RILs population, which is derived from a cross of smooth and densely netted parents, confirmed the intricacy of this trait and the involvement of major locus and several other interacting QTLs. A major netting QTL on chromosome 2 co-localized with results from two additional populations, paving the way for future study toward identification of a causative gene for this trait.

Transcriptome profiling of non-climacteric 'yellow' melon during ripening: insights on sugar metabolism

BACKGROUND

The non-climacteric 'Yellow' melon (Cucumis melo, inodorus group) is an economically important crop and its quality is mainly determined by the sugar content. Thus, knowledge of sugar metabolism and its related pathways can contribute to the development of new field management and post-harvest practices, making it possible to deliver better quality fruits to consumers.

RESULTS

The RNA-seq associated with RT-qPCR analyses of four maturation stages were performed to identify important enzymes and pathways that are involved in the ripening profile of non-climacteric 'Yellow' melon fruit focusing on sugar metabolism. We identified 895 genes 10 days after pollination (DAP)-biased and 909 genes 40 DAP-biased. The KEGG pathway enrichment analysis of these differentially expressed (DE) genes revealed that 'hormone signal transduction', 'carbon metabolism', 'sucrose metabolism', 'protein processing in endoplasmic reticulum' and 'spliceosome' were the most differentially regulated processes occurring during melon development. In the sucrose metabolism, five DE genes are up-regulated and 12 are down-regulated during fruit ripening.

CONCLUSIONS

The results demonstrated important enzymes in the sugar pathway that are responsible for the sucrose content and maturation profile in non-climacteric 'Yellow' melon. New DE genes were first detected for melon in this study such as invertase inhibitor LIKE 3 (CmINH3), trehalose phosphate phosphatase (CmTPP1) and trehalose phosphate synthases (CmTPS5, CmTPS7, CmTPS9). Furthermore, the results of the protein-protein network interaction demonstrated general characteristics of the transcriptome of young and full-ripe melon and provide new perspectives for the understanding of ripening.

Centromeres of Cucumis melo L. comprise Cmcent and two novel repeats, CmSat162 and CmSat189

Centromeres are prerequisite for accurate segregation and are landmarks of primary constrictions of metaphase chromosomes in eukaryotes. In melon, high-copy-number satellite DNAs (SatDNAs) were found at various chromosomal locations such as centromeric, pericentromeric, and subtelomeric regions. In the present study, utilizing the published draft genome sequence of melon, two new SatDNAs (CmSat162 and CmSat189) of melon were identified and their chromosomal distributions were confirmed using fluorescence in situ hybridization. DNA probes prepared from these SatDNAs were successfully hybridized to melon somatic and meiotic chromosomes. CmSat162 was located on 12 pairs of melon chromosomes and co-localized with the centromeric repeat, Cmcent, at the centromeric regions. In contrast, CmSat189 was found to be located not only on centromeric regions but also on specific regions of the chromosomes, allowing the characterization of individual chromosomes of melon. It was also shown that these SatDNAs were transcribed in melon. These results suggest that CmSat162 and CmSat189 might have some functions at the centromeric regions.

A tunable filter for high molecular weight DNA selection and linked-read sequencing

In third generation sequencing, the production of quality data requires the selection of molecules longer than ∼20 kbp, but the size selection threshold of most purification technologies is smaller than this target. Here, we describe a technology operated in a capillary with a tunable selection threshold in the range of 3 to 40 kbp controlled by an electric field. We demonstrate that the selection cut-off is sharp, the purification yield is high, and the purification throughput is scalable. We also provide an analytical model that the actuation settings of the filter. The selection of high molecular weight genomic DNA from the melon Cucumis melo L., a diploid organism of ∼0.45 Gbp, is then reported. Linked-read sequencing data show that the N50 phase block size, which scores the correct representation of two chromosomes, is enhanced by a factor of 2 after size selection, establishing the relevance and versatility of our technology.

Molecular cloning and functional characterization of CmFT (FLOWERING LOCUS T) from Cucumis melo L

CmFT homologous gene in muskmelon was obtained by homologous cloning, introducing CmFT gene by Agrobacterium mediated transformation. The results of subcellular localization showed that CmFT protein was expressed in cytoplasm and nucleus. qRTPCR results showed that the expression levels of AtLFY, AtFT, AtCO, AtFLC, AtSOC1 and AtAP1 were upregulated in the 35S::MeFT Arabidopsis line. The CmFT gene was introduced into wild-type Arabidopsis by Agrobacterium-mediated transformation, and the growth status of T2 transgenic Arabidopsis thaliana and wild-type A. thaliana was observed. The results showed that wild-type Arabidopsis began to bolt on the 25th day after sowing, we can initially confirm that the FT gene of melon can promote the early flowering of melon in the growth and development of melon.

Genome-wide identification and characterization of long non-coding RNAs involved in fruit ripening and the climacteric in Cucumis melo

BACKGROUND

Cucumis melo is a suitable study material for investigation of fruit ripening owing to its climacteric nature. Long non-coding RNAs have been linked to many important biological processes, such as fruit ripening, flowering time regulation, and abiotic stress responses in plants. However, knowledge of the regulatory roles of lncRNAs underlying the ripening process in C. melo are largely unknown. In this study the complete transcriptome of Cucumis melo L. cv. Hetao fruit at four developmental stages was sequenced and analyzed. The potential role of lncRNAs was predicted based on the function of differentially expressed target genes and correlated genes.

RESULTS

In total, 3857 lncRNAs were assembled and annotated, of which 1601 were differentially expressed between developmental stages. The target genes of these lncRNAs and the regulatory relationship (cis- or trans-acting) were predicted. The target genes were enriched with GO terms for biological process, such as response to auxin stimulus and hormone biosynthetic process. Enriched KEGG pathways included plant hormone signal transduction and carotenoid biosynthesis. Co-expression network construction showed that LNC_002345 and LNC_000154, which were highly expressed, might co-regulate with mutiple genes associated with auxin signal transduction and acted in the same pathways. We identified lncRNAs (LNC_000987, LNC_000693, LNC_001323, LNC_003610, LNC_001263 and LNC_003380) that were correlated with fruit ripening and the climacteric, and may participate in the regulation of ethylene biosynthesis and metabolism and the ABA signaling pathway. A number of crucial transcription factors, such as ERFs, WRKY70, NAC56, and NAC72, may also play important roles in the regulation of fruit ripening in C. melo.

CONCLUSIONS

Our results predict the regulatory functions of the lncRNAs during melon fruit development and ripening, and 142 highly expressed lncRNAs (average FPKM > 100) were identified. These lncRNAs participate in the regulation of auxin signal transduction, ethylene, sucrose biosynthesis and metabolism, the ABA signaling pathway, and transcription factors, thus regulating fruit development and ripening.