The genome of melon (Cucumis melo L.)

Identification of powdery mildew resistance quantitative trait loci in melon and development of resistant near-isogenic lines through marker-assisted backcrossing

BACKGROUND

Melon (Cucumis melo L.), an important cucurbit crop, faces production limitations due to powdery mildew (PM). Developing resistant varieties offers a sustainable, genetics-based alternative to chemical treatments. Therefore, identifying PM resistance quantitative trait loci (QTL) and creating trait-associated markers are essential for efficient melon PM resistance improvement through marker-assisted backcrossing (MABC).

RESULTS

Three F2 populations, A6, B2, and C4, were generated for QTL mapping of PM resistance. Major QTL were identified on chromosome 2 in A6, chromosome 5 in B2, and chromosomes 5 and 12 in C4. A series of TaqMan® assays targeting regions on chromosomes 2, 5, and 12 were developed and validated for foreground and recombinant selection, complemented by the double digest restriction-site associated DNA genotyping system to evaluate the recurrent parent genome recovery. Three MABC programs using resistant donor parents from A6 and C4 crossed with elite susceptible recurrent parents with green and orange fruit flesh were implemented. After two to three cycles of MABC, individual QTL was successfully introgressed into elite genetic backgrounds, giving six PM resistance lines in each green- and orange-fleshed background. PM inoculation on the twelve near-isogenic lines confirmed their resistance to PM.

CONCLUSIONS

We have identified major PM resistance QTL for melon on chromosomes 2, 5, and 12 and have introgressed individual QTL to elite genetic backgrounds using MABC in three and a half years. This study demonstrates the power of combining high-throughput genotyping with breeding efforts and showcases the efficiency of molecular breeding.

Genome-Wide Identification and Expression Analysis of the Melon Aldehyde Dehydrogenase (ALDH) Gene Family in Response to Abiotic and Biotic Stresses

Through the integration of genomic information, transcriptome sequencing data, and bioinformatics methods, we conducted a comprehensive identification of the ALDH gene family in melon. We explored the impact of this gene family on melon growth, development, and their expression patterns in various tissues and under different stress conditions. Our study discovered a total of 17 ALDH genes spread across chromosomes 1, 2, 3, 4, 5, 7, 8, 11, and 12 in the melon genome. Through a phylogenetic analysis, these genes were classified into 10 distinct subfamilies. Notably, genes within the same subfamily exhibited consistent gene structures and conserved motifs. Our study discovered a pair of fragmental duplications within the melon ALDH gene. Furthermore, there was a noticeable collinearity relationship between the melon's ALDH gene and that of Arabidopsis (12 times), and rice (3 times). Transcriptome data reanalysis revealed that some ALDH genes consistently expressed highly across all tissues and developmental stages, while others were tissue- or stage-specific. We analyzed the ALDH gene's expression patterns under six stress types, namely salt, cold, waterlogged, powdery mildew, Fusarium wilt, and gummy stem blight. The results showed differential expression of CmALDH2C4 and CmALDH11A3 under all stress conditions, signifying their crucial roles in melon growth and stress response. RT-qPCR (quantitative reverse transcription PCR) analysis further corroborated these findings. This study paves the way for future genetic improvements in melon molecular breeding.

Genome-Wide Identification of the Cyclic Nucleotide-Gated Ion Channel Gene Family and Expression Profiles Under Low-Temperature Stress in Luffa cylindrica L

Cyclic nucleotide-gated ion channels (CNGCs) are cell membrane channel proteins for calcium ions. They have been reported to play important roles in survival and in the responses to environmental factors in various plants. However, little is known about the CNGC family and its functions in luffa (Luffa cylindrica L.). In this study, a bioinformatics-based method was used to identify members of the CNGC gene family in L. cylindrica. In total, 20 LcCNGCs were detected, and they were grouped into five subfamilies (I, II, Ⅲ, IV-a, and IV-b) in a phylogenetic analysis with CNGCs from Arabidopsis thaliana (20 AtCNGCs) and Momordica charantia (17 McCNGCs). The 20 LcCNGC genes were unevenly distributed on 11 of the 13 chromosomes in luffa, with none on Chromosomes 1 and 5. The members of each subfamily encoded proteins with highly conserved functional domains. An evolutionary analysis of CNGCs in luffa revealed three gene losses and a motif deletion. An examination of gene replication events during evolution indicated that two tandemly duplicated gene pairs were the primary driving force behind the evolution of the LcCNGC gene family. PlantCARE analyses of the LcCNGC promoter regions revealed various cis-regulatory elements, including those responsive to plant hormones (abscisic acid, methyl jasmonate, and salicylic acid) and abiotic stresses (light, drought, and low temperature). The presence of these cis-acting elements suggested that the encoded CNGC proteins may be involved in stress responses, as well as growth and development. Transcriptome sequencing (RNA-seq) analyses revealed tissue-specific expression patterns of LcCNGCs in various plant parts (roots, stems, leaves, flowers, and fruit) and the upregulation of some LcCNGCs under low-temperature stress. To confirm the accuracy of the RNA-seq data, 10 cold-responsive LcCNGC genes were selected for verification by quantitative real-time polymerase chain reaction (RT-qPCR) analysis. Under cold conditions, LcCNGC4 was highly upregulated (>50-fold increase in its transcript levels), and LcCNGC3, LcCNGC6, and LcCNGC13 were upregulated approximately 10-fold. Our findings provide new information about the evolution of the CNGC family in L. cylindrica and provide insights into the functions of the encoded CNGC proteins.

FoodAtlas: Automated knowledge extraction of food and chemicals from literature

Automated generation of knowledge graphs that accurately capture published information can help with knowledge organization and access, which have the potential to accelerate discovery and innovation. Here, we present an integrated pipeline to construct a large-scale knowledge graph using large language models in an active learning setting. We apply our pipeline to the association of raw food, ingredients, and chemicals, a domain that lacks such knowledge resources. By using an iterative active learning approach of 4120 manually curated premise-hypothesis pairs as training data for ten consecutive cycles, the entailment model extracted 230,848 food-chemical composition relationships from 155,260 scientific papers, with 106,082 (46.0 %) of them never been reported in any published database. To augment the knowledge incorporated in the knowledge graph, we further incorporated information from 5 external databases and ontology sources. We then applied a link prediction model to identify putative food-chemical relationships that were not part of the constructed knowledge graph. Validation of the 443 hypotheses generated by the link prediction model resulted in 355 new food-chemical relationships, while results show that the model score correlates well (R2 = 0.70) with the probability of a novel finding. This work demonstrates how automated learning from literature at scale can accelerate discovery and support practical applications through reproducible, evidence-based capture of latent interactions of diverse entities, such as food and chemicals.

Loss of Lateral suppressor gene is associated with evolution of root nodule symbiosis in Leguminosae

BACKGROUND

Root nodule symbiosis (RNS) is a fascinating evolutionary event. Given that limited genes conferring the evolution of RNS in Leguminosae have been functionally validated, the genetic basis of the evolution of RNS remains largely unknown. Identifying the genes involved in the evolution of RNS will help to reveal the mystery.

RESULTS

Here, we investigate the gene loss event during the evolution of RNS in Leguminosae through phylogenomic and synteny analyses in 48 species including 16 Leguminosae species. We reveal that loss of the Lateral suppressor gene, a member of the GRAS-domain protein family, is associated with the evolution of RNS in Leguminosae. Ectopic expression of the Lateral suppressor (Ls) gene from tomato and its homolog MONOCULM 1 (MOC1) and Os7 from rice in soybean and Medicago truncatula result in almost completely lost nodulation capability. Further investigation shows that Lateral suppressor protein, Ls, MOC1, and Os7 might function through an interaction with NODULATION SIGNALING PATHWAY 2 (NSP2) and CYCLOPS to repress the transcription of NODULE INCEPTION (NIN) to inhibit the nodulation in Leguminosae. Additionally, we find that the cathepsin H (CTSH), a conserved protein, could interact with Lateral suppressor protein, Ls, MOC1, and Os7 and affect the nodulation.

CONCLUSIONS

This study sheds light on uncovering the genetic basis of the evolution of RNS in Leguminosae and suggests that gene loss plays an essential role.

Superstatistics Applied to Cucurbitaceae DNA Sequences

The short and long statistical correlations are essential in the genomic sequence. Such correlations are long-range for introns, whereas, for exons, these are short. In this study, we employed superstatistics to investigate correlations and fluctuations in the distribution of nucleotide sequence lengths of the Cucurbitaceae family. We established a time series for exon sizes to probe these correlations and fluctuations. We used data from the National Center for Biotechnology Information (NCBI) gene database to extract the temporal evolution of exon sizes, measured in terms of the number of base pairs (bp). To assess the model's viability, we utilized a timescale extraction method to determine the statistical properties of our time series, including the local distribution and fluctuations, which provide the exon size distributions based on the q-Gamma and inverse q-Gamma distributions. From the Bayesian statistics standpoint, both distributions are excellent for capturing the correlations and fluctuations from the data.

Identification of MYC genes in four Cucurbitaceae species and their roles in the response to temperature stress

BACKGROUND

Myelocytomatosis (MYC) transcription factors are crucial mediators of the response of plants to environmental stresses through via binding to DNA regulatory regions. However, few systematic characterizations of MYC genes are available in Cucurbitaceae species.

RESULTS

In this study, we identified 10, 8, 12, and 10 MYC genes in Cucumis sativus, Cucumis melo, Citrullus lanatus, and Benincasa hispida, respectively. Characterization revealed that all of the MYC proteins contain a highly conserved H4-V5-E6-E8-R9-R11-R12 sequence, which is essential for the binding of DNA regulatory regions. Evolutionary analysis enabled us to categorize 40 predicted MYC proteins from seven species into five distinct groups and revealed that the expansion of the MYC genes occurred before the divergence of monocots and dicots. The upstream promoter regions of the MYC genes contain a variety of developmental, stress, and hormone-responsive regulatory elements. The expression of cucumber MYC genes varies significantly across organs, with particularly high expression of CsaV3_3G001710 observed across all organs. Transcriptomic analysis revealed that certain cucumber MYC genes undergo specific upregulation or downregulation in response to both biotic and abiotic stressors. In particular, under temperature stress, the cucumber genes CsaV3_3G007980 and CsaV3_3G001710 were significantly upregulated. Interestingly, the homologs of these two genes in C. lanatus presented a similar expression pattern to that in C. sativus, whereas in B. hispida, they presented the opposite pattern, i.e., significant downregulation. These findings indicated that these two genes indeed respond to temperature stress but with different expression patterns, highlighting the divergent functions of homologous genes across different species.

CONCLUSIONS

This study analyzed the size and composition of the MYC gene family in four Cucurbitaceae species and investigated stress-responsive expression profiles, especially under temperature stress. All the results showed that MYC genes play important roles in development and stress responses, laying a theoretical foundation for further investigations of these response mechanisms.

Genome-wide identification of the GATA gene family in melon (Cucumis melo) and analysis of their expression characteristics under biotic and abiotic stresses

GATA transcription factors are an important class of transcription factors in plants, known for their roles in tissue development, signal transduction, and responses to biotic and abiotic stresses. To date, there have been no reports on the GATA gene family in melon (Cucumis melo). In this study, 24 CmGATA genes were identified from the melon genome. These family members exhibit significant differences in protein length, molecular weight, and theoretical isoelectric point and are primarily located in the nucleus. Based on the classification of Arabidopsis thaliana GATA members, the phylogenetic tree divided them into four groups: group I, group II, group III, and group IV, containing 10, 8, 4, and 2 genes, respectively. Notably, CmGATA genes within the same group have highly conserved protein motifs and similar exon-intron structures. The CmGATA family members are unevenly distributed across 10 chromosomes, with six pairs of segmentally duplicated genes and one pair of tandemly duplicated genes, suggesting that gene duplication may be the primary factor in the expansion of the CmGATA family. Melon shares 21, 4, 38, and 34 pairs of homologous genes with A. thaliana, Oryza sativa, Cucumis sativus, and Citrullus lanatus, respectively. The promoter regions are enriched with various cis-acting elements related to growth and development (eight types), hormone regulation (nine types), and stress responses (six types). Expression patterns indicate that different CmGATA family members are significantly expressed in seeds, roots, stems, leaves, tendrils, mesocarp, and epicarp, exhibiting distinct tissue-specific expression characteristics. Quantitative fluorescence analysis revealed that five genes, CmGATA3, CmGATA7, CmGATA16, CmGATA22, and CmGATA24, may be highly active under 48-h drought stress, while CmGATA1 and CmGATA22 may enhance melon resistance to heavy metal lead stress. Additionally, CmGATA22 and CmGATA24 are suggested to regulate melon resistance to Fusarium wilt infection. CmGATA22 appears to comprehensively regulate melon responses to both biotic and abiotic stresses. Lastly, potential protein interaction networks were predicted for the CmGATA family members, identifying CmGATA8 as a potential hub gene and predicting 2,230 target genes with enriched GO functions. This study preliminarily explores the expression characteristics of CmGATA genes under drought stress, heavy metal lead stress, and Fusarium wilt infection, providing a theoretical foundation for molecular mechanisms in melon improvement and stress resistance.

Genetic analysis and preliminary mapping by BSA-seq of the CmSR gene regulating the spotted rind trait in melon (Cucumis melo L.)

Melon (Cucumis melo L.) is an economically important horticultural crop. Spotted rind at maturity is an important appearance quality trait in melons. However, the gene controlling this trait remains unknown. In this study, the inheritance pattern of this trait was explored, and the candidate gene underlying this trait was also successfully identified. Genetic analysis showed that a single dominant gene, Cucumis melo Spotted Rind (CmSR), regulates the spotted rind trait. A preliminary genetic mapping analysis was conducted based on a BSA-seq approach. The CmAPRR2 gene was identified to be linked with the spotted rind trait and was located on the short arm of chromosome 4. It harbored two single-nucleotide mutations (chr4: 687014 G/A and chr4: 687244 C/A) in the non-spotted line 'Yellow 2', which may result in the alternative splicing of the transcript and an amino acid change in the respective protein, from proline to glutamine, respectively. Moreover, marker SNP687014-G/A was developed and co-segregated with the spotted rind trait. Therefore, it is speculated that the CmAPRR2 gene may be involved in the regulation of the spotted rind trait in melon. This study provides a theoretical foundation for further research on the gene regulatory mechanism of the rind color in melon.

Chromosome-level genome assembly and annotation of Flueggea virosa (Phyllanthaceae)

Flueggea virosa (Roxb. ex Willd.) Royle, an evergreen shrub and small tree in the Phyllanthaceae family, holds significant potential in garden landscaping and pharmacological applications. However, the lack of genomic data has hindered further scientific understanding of its horticultural and medicinal values. In this study, we have assembled a haplotype-resolved genome of F. virosa for the first time. The two haploid genomes, named haplotype A genome and haplotype B genome, are 487.33 Mb and 477.53 Mb in size, respectively, with contig N50 lengths of 31.45 Mb and 32.81 Mb. More than 99% of the assembled sequences were anchored to 13 pairs of pseudo-chromosomes. Furthermore, 21,587 and 21,533 protein-coding genes were predicted in haplotype A and haplotype B genomes, respectively. The availability of this chromosome-level genome fills the gap in genomic data for F. virosa and provides valuable resources for molecular studies of this species, supporting future research on speciation, functional genomics, and comparative genomics within the Phyllanthaceae family.

Chromosome-level and haplotype-resolved genome assembly of Dracaena cambodiana (Asparagaceae)

Dracaena cambodiana Pierre ex Gagn. (Asparagaceae) is the source plant of Dragon's blood and has high ornamental values in gardening. Currently, this species is classified as the second-class state-protected species in the National Key Protected Wild Plants (NKPWP) of China. However, limited genomic data has hindered a more comprehensive scientific understanding of the processes involved in the production of Dragon's blood and the related conservation genomics research. In this study, we assembled a haplotype-resolved genome of D. cambodiana. The haploid genomes, haplotype A and haplotype B, are 1,015.22 Mb and 1,003.13 Mb in size, respectively. The completeness of haplotype A and haplotype B genomes was 98.60% and 98.20%, respectively, using the "embryophyta_10" dataset. Haplotype A and haplotype B genomes contained 27,361 and 27,066 protein-coding genes, respectively, with nearly all being functionally annotated. These findings provide new insights into the genomic characteristics of D. cambodiana and will offer additional genomic resources for studying the biosynthesis mechanism of Dragon's blood and the horticultural application of Dragon trees.

A deletion in FLS2 and its expansion after domestication caused global dissemination of melon cultivars defective in flagellin recognition

FLAGELLIN SENSING 2 (FLS2) encodes a pattern recognition receptor that perceives bacterial flagellin. While putative FLS2 orthologs are broadly conserved in plants, their functional characterization remains limited. Here, we report the identification of orthologs in cucumber (Cucumis sativus) and melon (C. melo), named CsFLS2 and CmFLS2, respectively. Homology searching identified CsFLS2, and virus-induced gene silencing (VIGS) demonstrated that CsFLS2 is required for flg22-triggered ROS generation. Interestingly, genome re-sequencing of melon cv. Lennon and subsequent genomic PCR revealed that Lennon has two CmFLS2 haplotypes, haplotype I encoding full-length CmFLS2 and haplotype II encoding a truncated form. We show that VIGS-mediated knockdown of CmFLS2 haplotype I resulted in a significant reduction in both flg22-triggered ROS generation and immunity to a bacterial pathogen in melon cv. Lennon. Remarkably, genomic PCR of CmFLS2 revealed that 68% of tested commercial melon cultivars possess only CmFLS2 haplotype II: these cultivars thus lack functional CmFLS2. To explore evolutionary aspects of CmFLS2 haplotype II occurrence, we genotyped the CmFLS2 locus in 142 melon accessions by genomic PCR and analyzed 437 released sequences. The results suggest that CmFLS2 haplotype II is derived from C. melo subsp. melo. Furthermore, we suggest that the proportion of CmFLS2 haplotype II increased among the improved melo group compared with the primitive melo group. Collectively, these findings suggest that the deleted FLS2 locus generated in the primitive melo subspecies expanded after domestication, resulting in the spread of commercial melon cultivars defective in flagellin recognition, which is critical for bacterial immunity.

Characterization and stress-responsive regulation of CmPHT1 genes involved in phosphate uptake and transport in Melon (Cucumis melo L.)

BACKGROUND

Phosphorus (P) deficiency, a major nutrient stress, greatly hinders plant growth. Phosphate (Pi) uptake in plant roots relies on PHT1 family transporters. However, melon (Cucumis melo L.) lacks comprehensive identification and characterization of PHT1 genes, particularly their response patterns under diverse stresses.

RESULTS

This study identified and analyzed seven putative CmPHT1 genes on chromosomes 3, 4, 5, 6, and 7 using the melon genome. Phylogenetic analysis revealed shared motifs, domain compositions, and evolutionary relationships among genes with close histories. Exon number varied from 1 to 3. Collinearity analysis suggested segmental and tandem duplications as the primary mechanisms for CmPHT1 gene family expansion. CmPHT1;4 and CmPHT1;5 emerged as a tandemly duplicated pair. Analysis of cis-elements in CmPHT1 promoters identified 14 functional categories, including putative PHR1-binding sites (P1BS) in CmPHT1;4, CmPHT1;6, and CmPHT1;7. We identified that three WRKY transcription factors regulated CmPHT1;5 expression by binding to its W-box element. Notably, CmPHT1 promoters harbored cis-elements responsive to hormones and abiotic factors. Different stresses regulated CmPHT1 expression differently, suggesting that the adjusted expression patterns might contribute to plant adaptation.

CONCLUSIONS

This study unveils the characteristics, evolutionary diversity, and stress responsiveness of CmPHT1 genes in melon. These findings lay the foundation for in-depth investigations into their functional mechanisms in Cucurbitaceae crops.

Genome Assembly and Structural Variation Analysis of Luffa acutangula Provide Insights on Flowering Time and Ridge Development

Luffa spp. is an important worldwide cultivated vegetable and medicinal plant from the Cucurbitaceae family. In this study, we report a high-quality chromosome-level genome of the high-generation inbred line SG261 of Luffa acutangula. The genomic sequence was determined by PacBio long reads, Hi-C sequencing reads, and 10× Genomics sequencing, with an assembly size of 739.82 Mb, contig N50 of 18.38 Mb, and scaffold N50 of 56.08 Mb. The genome of L. acutangula SG261 was predicted to contain 27,312 protein-coding genes and 72.56% repetitive sequences, of which long terminal repeats (LTRs) were an important form of repetitive sequences, accounting for 67.84% of the genome. Phylogenetic analysis reveals that L. acutangula evolved later than Luffa cylindrica, and Luffa is closely related to Momodica charantia. Comparing the genome of L. acutangula SG261 and L. cylindrica with PacBio data, 67,128 high-quality structural variations (SVs) and 55,978 presence-absence variations (PAVs) were identified in SG261, resulting in 2424 and 1094 genes with variation in the CDS region, respectively, and there are 287 identical genes affected by two different structural variation analyses. In addition, we found that the transcription factor FY (FLOWERING LOCUS Y) families had a large expansion in L. acutangula SG261 (flowering in the morning) compared to L. cylindrica (flowering in the afternoon), which may result in the early flowering time in L. acutangula SG261. This study provides valuable reference for the breeding of and pan-genome research into Luffa species.

Melon CmRLCK VII-8 kinase genes CmRLCK27, CmRLCK30 and CmRLCK34 modulate resistance against bacterial and fungal diseases in Arabidopsis

Receptor-like cytoplasmic kinases (RLCKs) represent a distinct class of receptor-like kinases crucial for various aspects of plant biology, including growth, development, and stress responses. This study delves into the characterization of RLCK VII-8 members within cucurbits, particularly in melon, examining both structural features and the phylogenetic relationships of these genes/proteins. The investigation extends to their potential involvement in disease resistance by employing ectopic overexpression in Arabidopsis. The promoters of CmRLCK VII-8 genes harbor multiple phytohormone- and stress-responsive cis-acting elements, with the majority (excluding CmRLCK39) displaying upregulated expression in response to defense hormones and fungal infection. Subcellular localization studies reveal that CmRLCK VII-8 proteins predominantly reside on the plasma membrane, with CmRLCK29 and CmRLCK30 exhibiting additional nuclear distribution. Notably, Arabidopsis plants overexpressing CmRLCK30 manifest dwarfing and delayed flowering phenotypes. Overexpression of CmRLCK27, CmRLCK30, and CmRLCK34 in Arabidopsis imparts enhanced resistance against Botrytis cinerea and Pseudomonas syringae pv. tomato DC3000, concomitant with the strengthened expression of defense genes and reactive oxygen species accumulation. The CmRLCK VII-8 members actively participate in chitin- and flg22-triggered immune responses. Furthermore, CmRLCK30 interacts with CmMAPKKK1 and CmARFGAP, adding a layer of complexity to the regulatory network. In summary, this functional characterization underscores the regulatory roles of CmRLCK27, CmRLCK30, and CmRLCK34 in immune responses by influencing pathogen-induced defense gene expression and ROS accumulation.

Deciphering recursive polyploidization in Lamiales and reconstructing their chromosome evolutionary trajectories

Lamiales is an order of core eudicots with abundant diversity, and many Lamiales plants have important medicinal and ornamental values. Here, we comparatively reanalyzed 11 Lamiales species with well-assembled genome sequences and found evidence that Lamiales plants, in addition to a hexaploidization or whole-genome triplication (WGT) shared by core eudicots, experienced further polyploidization events, establishing new groups in the order. Notably, we identified a whole-genome duplication (WGD) occurred just before the split of Scrophulariaceae from the other Lamiales families, such as Acanthaceae, Bignoniaceae, and Lamiaceae, suggesting its likely being the causal reason for the establishment and fast divergence of these families. We also found that a WGT occurred ∼68 to 78 million years ago (Mya), near the split of Oleaceae from the other Lamiales families, implying that it may have caused their fast divergence and the establishment of the Oleaceae family. Then, by exploring and distinguishing intra- and intergenomic chromosomal homology due to recursive polyploidization and speciation, respectively, we inferred that the Lamiales ancestral cell karyotype had 11 proto-chromosomes. We reconstructed the evolutionary trajectories from these proto-chromosomes to form the extant chromosomes in each Lamiales plant under study. We must note that most of the inferred 11 proto-chromosomes, duplicated during a WGD thereafter, have been well preserved in jacaranda (Jacaranda mimosifolia) genome, showing the credibility of the present inference implementing a telomere-centric chromosome repatterning model. These efforts are important to understand genome repatterning after recursive polyploidization, especially shedding light on the origin of new plant groups and angiosperm cell karyotype evolution.

CmERF1 acts as a positive regulator of fruits and leaves growth in melon (Cucumis melo L.)

Melon (Cucumis melo L.) is an important horticultural and economic crop. ETHYLENE RESPONSE FACTOR1 (ERF1) plays an important role in regulating plant development, and the resistance to multiple biotic and abiotic stresses. In this study, developmental biology, molecular biology and biochemical assays were performed to explore the biological function of CmERF1 in melon. Abundant transcripts of CmERF1 were found in ovary at green-yellow bud (GYB) and rapid enlargement (ORE) stages. In CmERF1 promoter, the cis-regulatory elements for indoleacetic acid (IAA), methyl jasmonate (MeJA), salicylic acid (SA), abscisic acid (ABA), gibberellic acid (GA), light and low temperature responses were found. CmERF1 could be significantly induced by ethylene, IAA, MeJA, SA, ABA, and respond to continuous light and low temperature stresses in melon. Ectopic expression of CmERF1 increased the length of siliqua and carpopodium, and expanded the size of leaves in Arabidopsis. Knockdown of CmERF1 led to smaller ovary at anthesis, mature fruit and leaves in melon. In CmERF1-RNAi #2 plants, 75 genes were differently expressed compared with control, and the promoter regions of 28 differential expression genes (DEGs) contained the GCC-box (AGCCGCC) or DRE (A/GCCGAC) cis-acting elements of CmERF1. A homolog of cell division cycle protein 48 (CmCDC48) was proved to be the direct target of CmERF1 by the yeast one-hybrid assay and dual-luciferase (LUC) reporter (DLR) system. These results indicated that CmERF1 was able to promote the growth of fruits and leaves, and involved in multiple hormones and environmental signaling pathways in melon.

Chromosome-level genome assembly and population genomics reveals crucial selection for subgynoecy development in chieh-qua

Chieh-qua is an important cucurbit crop and very popular in South China and Southeast Asia. Despite its significance, its genetic basis and domestication history are unclear. In this study, we have successfully generated a chromosome-level reference genome assembly for the chieh-qua 'A36' using a hybrid assembly strategy that combines PacBio long reads and Illumina short reads. The assembled genome of chieh-qua is approximately 953.3 Mb in size and is organized into 12 chromosomes, with contig N50 of 6.9 Mb and scaffold N50 of 68.2 Mb. Notably, the chieh-qua genome is comparable in size to the wax gourd genome. Through gene prediction analysis, we have identified a total of 24 593 protein-coding genes in the A36 genome. Additionally, approximately 56.6% (539.3 Mb) of the chieh-qua genome consists of repetitive sequences. Comparative genome analysis revealed that chieh-qua and wax gourd are closely related, indicating a close evolutionary relationship between the two species. Population genomic analysis, employing 129 chieh-qua accessions and 146 wax gourd accessions, demonstrated that chieh-qua exhibits greater genetic diversity compared to wax gourd. We also employed the GWAS method to identify related QTLs associated with subgynoecy, an interested and important trait in chieh-qua. The MYB59 (BhiCQ0880026447) exhibited relatively high expression levels in the shoot apex of four subgynoecious varieties compared with monoecious varieties. Overall, this research provides insights into the domestication history of chieh-qua and offers valuable genomic resources for further molecular research.

Degenerate oligonucleotide primer MIG-seq: an effective PCR-based method for high-throughput genotyping

Next-generation sequencing (NGS) library construction often involves using restriction enzymes to decrease genome complexity, enabling versatile polymorphism detection in plants. However, plant leaves frequently contain impurities, such as polyphenols, necessitating DNA purification before enzymatic reactions. To overcome this problem, we developed a PCR-based method for expeditious NGS library preparation, offering flexibility in number of detected polymorphisms. By substituting a segment of the simple sequence repeat sequence in the MIG-seq primer set (MIG-seq being a PCR method enabling library construction with low-quality DNA) with degenerate oligonucleotides, we introduced variability in detectable polymorphisms across various crops. This innovation, named degenerate oligonucleotide primer MIG-seq (dpMIG-seq), enabled a streamlined protocol for constructing dpMIG-seq libraries from unpurified DNA, which was implemented stably in several crop species, including fruit trees. Furthermore, dpMIG-seq facilitated efficient lineage selection in wheat and enabled linkage map construction and quantitative trait loci analysis in tomato, rice, and soybean without necessitating DNA concentration adjustments. These findings underscore the potential of the dpMIG-seq protocol for advancing genetic analyses across diverse plant species.

Complete genome assembly provides a high-quality skeleton for pan-NLRome construction in melon

Melon (Cucumis melo L.), being under intensive domestication and selective breeding, displays an abundant phenotypic diversity. Wild germplasm with tolerance to stress represents an untapped genetic resource for discovery of disease-resistance genes. To comprehensively characterize resistance genes in melon, we generate a telomere-to-telomere (T2T) and gap-free genome of wild melon accession PI511890 (C. melo var. chito) with a total length of 375.0 Mb and a contig N50 of 31.24 Mb. The complete genome allows us to dissect genome architecture and identify resistance gene analogs. We construct a pan-NLRome using seven melon genomes, which include 208 variable and 18 core nucleotide-binding leucine-rich repeat receptors (NLRs). Multiple disease-related transcriptome analyses indicate that most up-regulated NLRs induced by pathogens are shell or cloud NLRs. The T2T gap-free assembly and the pan-NLRome not only serve as essential resources for genomic studies and molecular breeding of melon but also provide insights into the genome architecture and NLR diversity.

A nanopore-based cucumber genome assembly reveals structural variations at two QTLs controlling hypocotyl elongation

The Xishuangbanna (XIS) cucumber (Cucumis sativus var. xishuangbannanesis) is a semiwild variety that has many distinct agronomic traits. Here, long reads generated by Nanopore sequencing technology helped assembling a high-quality genome (contig N50 = 8.7 Mb) of landrace XIS49. A total of 10,036 structural/sequence variations (SVs) were identified when comparing with Chinese Long (CL), and known SVs controlling spines, tubercles, and carpel number were confirmed in XIS49 genome. Two QTLs of hypocotyl elongation under low light, SH3.1 and SH6.1, were fine-mapped using introgression lines (donor parent, XIS49; recurrent parent, CL). SH3.1 encodes a red-light receptor Phytochrome B (PhyB, CsaV3_3G015190). A ∼4 kb region with large deletion and highly divergent regions (HDRs) were identified in the promoter of the PhyB gene in XIS49. Loss of function of this PhyB caused a super-long hypocotyl phenotype. SH6.1 encodes a CCCH-type zinc finger protein FRIGIDA-ESSENTIAL LIKE (FEL, CsaV3_6G050300). FEL negatively regulated hypocotyl elongation but it was transcriptionally suppressed by long terminal repeats retrotransposon insertion in CL cucumber. Mechanistically, FEL physically binds to the promoter of CONSTITUTIVE PHOTOMORPHOGENIC 1a (COP1a), regulating the expression of COP1a and the downstream hypocotyl elongation. These above results demonstrate the genetic mechanism of cucumber hypocotyl elongation under low light.

Identification and Characterization of Evolutionary Conserved Muskmelon Non-coding miRNAs and Their Target Proteins

MicroRNAs (miRNAs) are short, endogenously encoded small RNAs, 18-26 nucleotides in length, which can posttranscriptionally regulate gene expression through translation inhibition or endonucleolytic cleavage. The muskmelon is one of the most widely cultivated fruits in the Cucurbitaceae family. Despite its significance, only 120 miRNAs from different families have been reported in muskmelon. In this study, we aimed to expand this knowledge base by predicting 40 new miRNAs in muskmelon using a spectrum of genomic-based tools. Precursor and mature sequences were obtained from microRNA registry database as reference and analyzed via the basic local alignment search tool (Blastn) for ESTs identification. After removing the non-coding sequences, the remaining candidate sequences were analyzed using MFOLD to generate secondary structures for the newly predicted miRNAs. Additionally, the predicted muskmelon miRNAs were validated using a set of five randomly chosen primers and RT-PCR. Through gene ontology (GO) analysis, we identified 89 targets associated with newly predicted muskmelon miRNAs. Transcription factor-coding genes play a crucial role in plant growth and development. Additionally, the miR4249 has been found to have the same targets in muskmelon that have been linked to cell signaling and transcription factors. The identified targets are integral for diverse biological processes including plant growth, development, metabolism, aging, disease resistance, and resistance to environmental stresses, such as salt, cold, and oxidative stress. As a result, the outcomes of this study demonstrate that this mechanism not only contributes to the production of a higher quality crop but also enhances overall production.

Novel Allelic Gene Variations in CmCLAVATA3 (CmCLV3) Were Identified in a Genetic Population of Melon (Cucumis melo L.)

Carpel number (CN) is an important trait affecting the fruit size and shape of melon, which plays a crucial role in determining the overall appearance and market value. A unique non-synonymous single nucleotide polymorphism (SNP) in CmCLAVATA3 (CmCLV3) is responsible for the variation of CN in C. melo ssp. agrestis (hereafter agrestis), but it has been unclear in C. melo ssp. melo (hereafter melo). In this study, one major locus controlling the polymorphism of 5-CN (multi-CN) and 3-CN (normal-CN) in melo was identified using bulked segregant analysis (BSA-seq). This locus was then fine-mapped to an interval of 1.8 Mb on chromosome 12 using a segregating population containing 1451 progeny. CmCLV3 is still present in the candidate region. A new allele of CmCLV3, which contains five other nucleotide polymorphisms, including a non-synonymous SNP in coding sequence (CDS), except the SNP reported in agrestis, was identified in melo. A cis-trans test confirmed that the candidate gene, CmCLV3, contributes to the variation of CNs in melo. The qRT-PCR results indicate that there is no significant difference in the expression level of CmCLV3 in the apical stem between the multi-CN plants and the normal-CN plants. Overall, this study provides a genetic resource for melon fruit development research and molecular breeding. Additionally, it suggests that melo has undergone similar genetic selection but evolved into an independent allele.

Deciphering the Genomic Characterization of the GGP Gene Family and Expression Verification of CmGGP1 Modulating Ascorbic Acid Biosynthesis in Melon Plants

Ascorbic acid (AsA), also known as vitamin C, is a well-known antioxidant found in living entities that plays an essential role in growth and development, as well as in defensive mechanisms. GDP-L-galactose phosphorylase (GGP) is a candidate gene regulating AsA biosynthesis at the translational and transcriptional levels in plants. In the current study, we conducted genome-wide bioinformatic analysis and pinpointed a single AsA synthesis rate-limiting enzyme gene in melon (CmGGP1). The protein prediction analysis depicted that the CmGGP1 protein does not have a signaling peptide or transmembrane structure and mainly functions in the chloroplast or nucleus. The constructed phylogenetic tree analysis in multispecies showed that the CmGGP1 protein has a highly conserved motif in cucurbit crops. The structural variation analysis of the CmGGP1 gene in different domesticated melon germplasms showed a single non-synonymous type-base mutation and indicated that this gene was selected by domestication during evolution. Wild-type (WT) and landrace (LDR) germplasms of melon depicted close relationships to each other, and improved-type (IMP) varieties showed modern domestication selection. The endogenous quantification of AsA content in both the young and old leaves of nine melon varieties exhibited the major differentiations for AsA synthesis and metabolism. The real-time quantitative polymerase chain reaction (qRT-PCR) analysis of gene co-expression showed that AsA biosynthesis in leaves was greater than AsA metabolic consumption, and four putative interactive genes (MELO3C025552.2, MELO3C007440.2, MELO3C023324.2, and MELO3C018576.2) associated with the CmGGP1 gene were revealed. Meanwhile, the CmGGP1 gene expression pattern was noticed to be up-regulated to varying degrees in different acclimated melons. We believe that the obtained results would provide useful insights for an in-depth genetic understanding of the AsA biosynthesis mechanism, aimed at the development of improving crop plants for melon.

Crop Landraces and Indigenous Varieties: A Valuable Source of Genes for Plant Breeding

Landraces and indigenous varieties comprise valuable sources of crop species diversity. Their utilization in plant breeding may lead to increased yield and enhanced quality traits, as well as resilience to various abiotic and biotic stresses. Recently, new approaches based on the rapid advancement of genomic technologies such as deciphering of pangenomes, multi-omics tools, marker-assisted selection (MAS), genome-wide association studies (GWAS), and CRISPR/Cas9 gene editing greatly facilitated the exploitation of landraces in modern plant breeding. In this paper, we present a comprehensive overview of the implementation of new genomic technologies and highlight their importance in pinpointing the genetic basis of desirable traits in landraces and indigenous varieties of annual, perennial herbaceous, and woody crop species cultivated in the Mediterranean region. The need for further employment of advanced -omic technologies to unravel the full potential of landraces and indigenous varieties underutilized genetic diversity is also indicated. Ultimately, the large amount of genomic data emerging from the investigation of landraces and indigenous varieties reveals their potential as a source of valuable genes and traits for breeding. The role of landraces and indigenous varieties in mitigating the ongoing risks posed by climate change in agriculture and food security is also highlighted.

Possibility of genome editing for melon breeding

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

Identification of WRKY Family Members and Characterization of the Low-Temperature-Stress-Responsive WRKY Genes in Luffa (Luffa cylindrica L.)

The plant-specific WRKY transcription factor family members have diverse regulatory effects on the genes associated with many plant processes. Although the WRKY proteins in Arabidopsis thaliana and other species have been thoroughly investigated, there has been relatively little research on the WRKY family in Luffa cylindrica, which is one of the most widely grown vegetables in China. In this study, we performed a genome-wide analysis to identify L. cylindrica WRKY genes, which were subsequently classified and examined in terms of their gene structures, chromosomal locations, promoter cis-acting elements, and responses to abiotic stress. A total of 62 LcWRKY genes (471-2238 bp) were identified and divided into three phylogenetic groups (I, II, and III), with group II further divided into five subgroups (IIa, IIb, IIc, IId, and IIe) in accordance with the classification in other plants. The LcWRKY genes were unevenly distributed across 13 chromosomes. The gene structure analysis indicated that the LcWRKY genes contained 0-11 introns (average of 4.4). Moreover, 20 motifs were detected in the LcWRKY proteins with conserved motifs among the different phylogenetic groups. Two subgroup IIc members (LcWRKY16 and LcWRKY31) contained the WRKY sequence variant WRKYGKK. Additionally, nine cis-acting elements related to diverse responses to environmental stimuli were identified in the LcWRKY promoters. The subcellular localization analysis indicated that three LcWRKY proteins (LcWRKY43, LcWRKY7, and LcWRKY23) are localized in the nucleus. The tissue-specific LcWRKY expression profiles reflected the diversity in LcWRKY expression. The RNA-seq data revealed the effects of low-temperature stress on LcWRKY expression. The cold-induced changes in expression were verified via a qRT-PCR analysis of 24 differentially expressed WRKY genes. Both LcWRKY7 and LcWRKY12 were highly responsive to the low-temperature treatment (approximately 110-fold increase in expression). Furthermore, the LcWRKY8, LcWRKY12, and LcWRKY59 expression levels increased by more than 25-fold under cold conditions. Our findings will help clarify the evolution of the luffa WRKY family while also providing valuable insights for future studies on WRKY functions.

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.

Whole-genome sequencing in medicinal plants: current progress and prospect

Advancements in genomics have dramatically accelerated the research on medicinal plants, and the development of herbgenomics has promoted the "Project of 1K Medicinal Plant Genome" to decipher their genetic code. However, it is difficult to obtain their high-quality whole genomes because of the prevalence of polyploidy and/or high genomic heterozygosity. Whole genomes of 123 medicinal plants were published until September 2022. These published genome sequences were investigated in this review, covering their classification, research teams, ploidy, medicinal functions, and sequencing strategies. More than 1,000 institutes or universities around the world and 50 countries are conducting research on medicinal plant genomes. Diploid species account for a majority of sequenced medicinal plants. The whole genomes of plants in the Poaceae family are the most studied. Almost 40% of the published papers studied species with tonifying, replenishing, and heat-cleaning medicinal effects. Medicinal plants are still in the process of domestication as compared with crops, thereby resulting in unclear genetic backgrounds and the lack of pure lines, thus making their genomes more difficult to complete. In addition, there is still no clear routine framework for a medicinal plant to obtain a high-quality whole genome. Herein, a clear and complete strategy has been originally proposed for creating a high-quality whole genome of medicinal plants. Moreover, whole genome-based biological studies of medicinal plants, including breeding and biosynthesis, were reviewed. We also advocate that a research platform of model medicinal plants should be established to promote the genomics research of medicinal plants.

High-quality Momordica balsamina genome elucidates its potential use in improving stress resilience and therapeutic properties of bitter gourd

Introduction

Momordica balsamina is the closest wild species that can be crossed with an important fruit vegetable crop, Momordica charantia, has immense medicinal value, and placed under II subclass of primary gene pool of bitter gourd. M. balsamina is tolerant to major biotic and abiotic stresses. Genome characterization of Momordica balsamina as a wild relative of bitter gourd will contribute to the knowledge of the gene pool available for improvement in bitter gourd. There is potential to transfer gene/s related to biotic resistance and medicinal importance from M. balsamina to M. charantia to produce high-quality, better yielding and stress tolerant bitter gourd genotypes.

Methods

The present study provides the first and high-quality chromosome-level genome assembly of M. balsamina with size 384.90 Mb and N50 30.96 Mb using sequence data from 10x Genomics, Nanopore, and Hi-C platforms.

Results

A total of 6,32,098 transposons elements; 2,15,379 simple sequence repeats; 5,67,483 transcription factor binding sites; 3,376 noncoding RNA genes; and 41,652 protein-coding genes were identified, and 4,347 disease resistance, 67 heat stress-related, 05 carotenoid-related, 15 salt stress-related, 229 cucurbitacin-related, 19 terpenes-related, 37 antioxidant activity, and 06 sex determination-related genes were characterized.

Conclusion

Genome sequencing of M. balsamina will facilitate interspecific introgression of desirable traits. This information is cataloged in the form of webgenomic resource available at http://webtom.cabgrid.res.in/mbger/. Our finding of comparative genome analysis will be useful to get insights into the patterns and processes associated with genome evolution and to uncover functional regions of cucurbit genomes.

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.

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.

Identification of Gene Responsible for Conferring Resistance against Race KN2 of Podosphaera xanthii in Melon

Powdery mildew caused by Podosphaera xanthii is a serious fungal disease which causes severe damage to melon production. Unlike with chemical fungicides, managing this disease with resistance varieties is cost effective and ecofriendly. But, the occurrence of new races and a breakdown of the existing resistance genes poses a great threat. Therefore, this study aimed to identify the resistance locus responsible for conferring resistance against P. xanthii race KN2 in melon line IML107. A bi-parental F2 population was used in this study to uncover the resistance against race KN2. Genetic analysis revealed the resistance to be monogenic and controlled by a single dominant gene in IML107. Initial marker analysis revealed the position of the gene to be located on chromosome 2 where many of the resistance gene against P. xanthii have been previously reported. Availability of the whole genome of melon and its R gene analysis facilitated the identification of a F-box type Leucine Rich Repeats (LRR) to be accountable for the resistance against race KN2 in IML107. The molecular marker developed in this study can be used for marker assisted breeding programs.

A highly diversified NLR cluster in melon contains homologs that confer powdery mildew and aphid resistance

Podosphaera xanthii is the main causal agent of powdery mildew (PM) on Cucurbitaceae. In Cucumis melo, the Pm-w resistance gene, which confers resistance to P. xanthii, is located on chromosome 5 in a cluster of nucleotide-binding leucine-rich repeat receptors (NLRs). We used positional cloning and transgenesis, to isolate the Pm-wWMR 29 gene encoding a coiled-coil NLR (CC-NLR). Pm-wWMR 29 conferred high level of resistance to race 1 of PM and intermediate level of resistance to race 3 of PM. Pm-wWMR 29 turned out to be a homolog of the Aphis gossypii resistance gene Vat-1PI 161375. We confirmed that Pm-wWMR 29 did not confer resistance to aphids, while Vat-1PI 161375 did not confer resistance to PM. We showed that both homologs were included in a highly diversified cluster of NLRs, the Vat cluster. Specific Vat-1PI 161375 and Pm-wWMR 29 markers were present in 10% to 13% of 678 accessions representative of wild and cultivated melon types worldwide. Phylogenic reconstruction of 34 protein homologs of Vat-1PI 161375 and Pm-wWMR 29 identified in 24 melon accessions revealed an ancestor with four R65aa-a specific motif in the LRR domain, evolved towards aphid and virus resistance, while an ancestor with five R65aa evolved towards PM resistance. The complexity of the cluster comprising the Vat/Pm-w genes and its diversity in melon suggest that Vat homologs may contribute to the recognition of a broad range of yet to be identified pests and pathogens.

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.

Automatic annotation of the bHLH gene family in plants

BACKGROUND

The bHLH transcription factor family is named after the basic helix-loop-helix (bHLH) domain that is a characteristic element of their members. Understanding the function and characteristics of this family is important for the examination of a wide range of functions. As the availability of genome sequences and transcriptome assemblies has increased significantly, the need for automated solutions that provide reliable functional annotations is emphasised.

RESULTS

A phylogenetic approach was adapted for the automatic identification and functional annotation of the bHLH transcription factor family. The bHLH_annotator, designed for the automated functional annotation of bHLHs, was implemented in Python3. Sequences of bHLHs described in literature were collected to represent the full diversity of bHLH sequences. Previously described orthologs form the basis for the functional annotation assignment to candidates which are also screened for bHLH-specific motifs. The pipeline was successfully deployed on the two Arabidopsis thaliana accessions Col-0 and Nd-1, the monocot species Dioscorea dumetorum, and a transcriptome assembly of Croton tiglium. Depending on the applied search parameters for the initial candidates in the pipeline, species-specific candidates or members of the bHLH family which experienced domain loss can be identified.

CONCLUSIONS

The bHLH_annotator allows a detailed and systematic investigation of the bHLH family in land plant species and classifies candidates based on bHLH-specific characteristics, which distinguishes the pipeline from other established functional annotation tools. This provides the basis for the functional annotation of the bHLH family in land plants and the systematic examination of a wide range of functions regulated by this transcription factor family.

Insights into Sucrose Metabolism and Its Ethylene-Dependent Regulation in Cucumis melo L

The melon (Cucumis melo L.), a fruit crop of significant economic importance, is prized for its sweet and succulent fruits. Among variations of soluble sugars, sucrose, a disaccharide composed of glucose and fructose, is a key carbohydrate present in melon fruits. The sucrose content also determines the quality and value of melon fruits. However, the accumulation of sucrose is a complex process involving the coordinated actions of multiple enzymes and pathways. In melon species, there are two types of fruit ripening modes including climacteric and non-climacteric. Due to this biological characteristic, melon is emerging as a good model for studying the ripening process. Ethylene is a well-known phytohormone regulating the ripening of climacteric fruits. Recently, a few studies have elucidated a primary ethylene-dependent signaling pathway of sucrose accumulation in melon fruits. This review aims to provide a careful overview of the sucrose biosynthesis pathways in melon. It is essential to understand the molecular mechanisms of sucrose metabolism as well as its regulation mode. The information will be useful for developing molecular marker-assisted breeding as well as genetic engineering strategies aiming to improve the sucrose content and quality of melon fruits. In addition, even though limited, the impacts of genetic background and environmental factors on sucrose accumulation in melon fruits are also discussed. These are useful for practical applications in melon cultivation and quality management.

Analyzing genomic variation in cultivated pumpkins and identification of candidate genes controlling seed traits

Pumpkins are important vegetable crops widely grown worldwide, and seeds are considered a popular nutraceutical food and an excellent source of protein, oil, and vitamins. Seed size is one of the most important targets for commercial breeding in Cucurbita species; studies have shown that pumpkin seed size variation has a similar trend with fruit size, shape, and seed yield. However, few studies have been conducted to identify genetic loci controlling seed-related traits in cultivated pumpkins. This study analyzed the genomic characteristics of pumpkin breeding materials of 321 Cucurbita accessions collected worldwide, including Cucurbita moschata, Cucurbita maxima, and Cucurbita pepo, using extensive single nucleotide polymorphisms obtained from the genotyping-by-sequencing method, significant genetic variations were identified within and between Cucurbita species. Four major cultivar fruit types were further revealed in C. moschata species, and significant differentiation patterns were detected in several chromosomal regions. A total of 15 significant loci associated with pumpkin seed traits were mapped through a genome-wide association approach; 32 genes previously reported to be associated with seed size regulation in Arabidopsis and Oryza sativa were located in the intervals defined by linkage disequilibrium. Through this study, we gained a deep understanding of the genomic variation distribution across Cucurbita species. The available genetic resources and the associated genetic contents could be used in commercial pumpkin breeding and will facilitate molecular marker-assisted selection in pumpkin seed trait improvement.

Genomic analysis and characterization of new loci associated with seed protein and oil content in soybeans

Breeding for increased protein without a reduction in oil content in soybeans [Glycine max (L.) Merr.] is a challenge for soybean breeders but an expected goal. Many efforts have been made to develop new soybean varieties with high yield in combination with desirable protein and/or oil traits. An elite line, R05-1415, was reported to be high yielding, high protein, and low oil. Several significant quantitative trait loci (QTL) for protein and oil were reported in this line, but many of them were unstable across environments or genetic backgrounds. Thus, a new study under multiple field environments using the Infinium BARCSoySNP6K BeadChips was conducted to detect and confirm stable genomic loci for these traits. Genetic analyses consistently detected a single major genomic locus conveying these two traits with remarkably high phenotypic variation explained (R2 ), varying between 24.2% and 43.5%. This new genomic locus is located between 25.0 and 26.7 Mb, distant from the previously reported QTL and did not overlap with other commonly reported QTL and the recently cloned gene Glyma.20G085100. Homolog analysis indicated that this QTL did not result from the paracentric chromosome inversion with an adjacent genomic fragment that harbors the reported QTL. The pleiotropic effect of this QTL could be a challenge for improving protein and oil simultaneously; however, a further study of four candidate genes with significant expressions in the seed developmental stages coupled with haplotype analysis may be able to pinpoint causative genes. The functionality and roles of these genes can be determined and characterized, which lay a solid foundation for the improvement of protein and oil content in soybeans.

Multi-Omics Analysis Reveals Intricate Gene Networks Involved in Female Development in Melon

Sexual differentiation is an important developmental phenomenon in cucurbits that directly affects fruit yield. The natural existence of multiple flower types in melon offers an inclusive structure for studying the molecular basis of sexual differentiation. The current study aimed to identify and characterize the molecular network involved in sex determination and female development in melon. Male and female pools separated by the F2 segregated generation were used for sequencing. The comparative multi-omics data revealed 551 DAPs and 594 DEGs involved in multiple pathways of melon growth and development, and based on functional annotation and enrichment analysis, we summarized four biological process modules, including ethylene biosynthesis, flower organ development, plant hormone signaling, and ubiquitinated protein metabolism, that are related to female development. Furthermore, the detailed analysis of the female developmental regulatory pathway model of ethylene biosynthesis, signal transduction, and target gene regulation identified some important candidates that might have a crucial role in female development. Two CMTs ((cytosine-5)-methyltransferase), one AdoHS (adenosylhomocysteinase), four ACSs (1-aminocyclopropane-1-carboxylic acid synthase), three ACOs (ACC oxidase), two ARFs (auxin response factor), four ARPs (auxin-responsive protein), and six ERFs (Ethylene responsive factor) were identified based on various female developmental regulatory models. Our data offer new and valuable insights into female development and hold the potential to offer a deeper comprehension of sex differentiation mechanisms in melon.

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.

Telomere-to-telomere genome assembly of melon (Cucumis melo L. var. inodorus) provides a high-quality reference for meta-QTL analysis of important traits

Melon is an important horticultural crop with extensive diversity in many horticultural groups. To explore its genomic diversity, it is necessary to assemble more high-quality complete genomes from different melon accessions. Meanwhile, a large number of QTLs have been mapped in several studies. Integration of the published QTLs onto a complete genome can provide more accurate information for candidate gene cloning. To address these problems, a telomere-to-telomere (T2T) genome of the elite melon landrace Kuizilikjiz (Cucumis melo L. var. inodorus) was de novo assembled and all the published QTLs were projected onto it in this study. The results showed that a high-quality Kuizilikjiz genome with the size of 379.2 Mb and N50 of 31.7 Mb was de novo assembled using the combination of short reads, PacBio high-fidelity long reads, Hi-C data, and a high-density genetic map. Each chromosome contained the centromere and telomeres at both ends. A large number of structural variations were observed between Kuizilikjiz and the other published genomes. A total of 1294 QTLs published in 67 studies were collected and projected onto the T2T genome. Several clustered, co-localized, and overlapped QTLs were determined. Furthermore, 20 stable meta-QTLs were identified, which significantly reduced the mapping intervals of the initial QTLs and greatly facilitated identification of the candidate genes. Collectively, the T2T genome assembly together with the numerous projected QTLs will not only broaden the high-quality genome resources but also provide valuable and abundant QTL information for cloning the genes controlling important traits in melon.

The haplotype-resolved T2T reference genome highlights structural variation underlying agronomic traits of melon

Melon (Cucumis melo L.) is an important vegetable crop that has an extensive history of cultivation. However, the genome of wild and semi-wild melon types that can be used for the analysis of agronomic traits is not yet available. Here we report a chromosome-level T2T genome assembly for 821 (C. melo ssp. agrestis var. acidulus), a semi-wild melon with two haplotypes of ~373 Mb and ~364 Mb, respectively. Comparative genome analysis discovered a significant number of structural variants (SVs) between melo (C. melo ssp. melo) and agrestis (C. melo ssp. agrestis) genomes, including a copy number variation located in the ToLCNDV resistance locus on chromosome 11. Genome-wide association studies detected a significant signal associated with climacteric ripening and identified one candidate gene CM_ac12g14720.1 (CmABA2), encoding a cytoplasmic short chain dehydrogenase/reductase, which controls the biosynthesis of abscisic acid. This study provides valuable genetic resources for future research on melon breeding.

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.

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.

Genomics of plant speciation

Studies of plants have been instrumental for revealing how new species originate. For several decades, botanical research has complemented and, in some cases, challenged concepts on speciation developed via the study of other organisms while also revealing additional ways in which species can form. Now, the ability to sequence genomes at an unprecedented pace and scale has allowed biologists to settle decades-long debates and tackle other emerging challenges in speciation research. Here, we review these recent genome-enabled developments in plant speciation. We discuss complications related to identification of reproductive isolation (RI) loci using analyses of the landscape of genomic divergence and highlight the important role that structural variants have in speciation, as increasingly revealed by new sequencing technologies. Further, we review how genomics has advanced what we know of some routes to new species formation, like hybridization or whole-genome duplication, while casting doubt on others, like population bottlenecks and genetic drift. While genomics can fast-track identification of genes and mutations that confer RI, we emphasize that follow-up molecular and field experiments remain critical. Nonetheless, genomics has clarified the outsized role of ancient variants rather than new mutations, particularly early during speciation. We conclude by highlighting promising avenues of future study. These include expanding what we know so far about the role of epigenetic and structural changes during speciation, broadening the scope and taxonomic breadth of plant speciation genomics studies, and synthesizing information from extensive genomic data that have already been generated by the plant speciation community.

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 characterization of the PLATZ gene family in watermelon (Citrullus lanatus L.) with putative functions in biotic and abiotic stress response

Plant AT-rich sequence and zinc-binding (PLATZ) proteins are plant-specific transcription factors involved in growth, development, and stress responses. Here, we conducted a genome-wide characterization of the watermelon ClPLATZ family and examined its expression responsiveness to defense hormones and pathogen infection along with putative functions in biotic and abiotic stress responses. The watermelon genome contains 12 putative ClPLATZ genes, encoding proteins with a characteristic PLATZ domain, and their promoters contain various cis-elements related to plant growth, development, phytohormones and stress response. The ClPLATZ genes, except ClPLATZ6, are differentially expressed in response to defense hormones (e.g., salicylic acid and methyl jasmonate) and fungal infections caused by Fusarium oxysporum f. sp. niveum and Stagonosporopsis cucurbitacearum. Most ClPLATZ proteins interact with other proteins (viz., ClDP, ClRPT2a, and ClRPC53). Among ClPLATZ proteins, ClPLATZ8, 9, 10, and 11 are predominately localized in the nucleus. ClPLATZ3 and 8 positively, but ClPLATZ11 negatively regulate resistance against Pseudomonas syringe pv. tomato DC3000 in transgenic Arabidopsis lines. ClPLATZ8 and 11 positively regulate stress tolerance to NaCl and mannitol during seed germination in transgenic Arabidopsis. In conclusion, the characterization of the ClPLATZ family provides insights into the biological functions of ClPLATZ genes in growth, development, and stress response in watermelon. Further, the involvement of certain ClPLATZ genes in biotic and abiotic stress response in transgenic Arabidopsis suggests their potential application in engineering stress-tolerant crops.

The Progression in Developing Genomic Resources for Crop Improvement

Sequencing technologies have rapidly evolved over the past two decades, and new technologies are being continually developed and commercialized. The emerging sequencing technologies target generating more data with fewer inputs and at lower costs. This has also translated to an increase in the number and type of corresponding applications in genomics besides enhanced computational capacities (both hardware and software). Alongside the evolving DNA sequencing landscape, bioinformatics research teams have also evolved to accommodate the increasingly demanding techniques used to combine and interpret data, leading to many researchers moving from the lab to the computer. The rich history of DNA sequencing has paved the way for new insights and the development of new analysis methods. Understanding and learning from past technologies can help with the progress of future applications. This review focuses on the evolution of sequencing technologies, their significant enabling role in generating plant genome assemblies and downstream applications, and the parallel development of bioinformatics tools and skills, filling the gap in data analysis techniques.