The wax gourd genomes offer insights into the genetic diversity and ancestral cucurbit karyotype

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.

Graph-Based Pangenome of Actinidia chinensis Reveals Structural Variations Mediating Fruit Degreening

Fruit ripening is associated with the degreening process (loss of chlorophyll) that occurs in most fruit species. Kiwifruit is one of the special species whose fruits may maintain green flesh by accumulating a large amount of chlorophyll even after ripening. However, little is known about the genetic variations related to the fruit degreening process. Here, a graph-based kiwifruit pangenome by analyzing 14 chromosome-scale haplotype-resolved genome assemblies from seven representative cultivars or lines in Actinidia chinensis is built. A total of 49,770 non-redundant gene families are identified, with core genes constituting 46.6%, and dispensable genes constituting 53.4%. A total of 84,591 non-redundant structural variations (SVs) are identified. The pangenome graph integrating both reference genome sequences and variant information facilitates the identification of SVs related to fruit color. The SV in the promoter of the AcBCM gene determines its high expression in the late developmental stage of fruits, which causes chlorophyll accumulation in the green-flesh fruits by post-translationally regulating AcSGR2, a key enzyme of chlorophyll catabolism. Taken together, a high-quality pangenome is constructed, unraveled numerous genetic variations, and identified a novel SV mediating fruit coloration and fruit quality, providing valuable information for further investigating genome evolution and domestication, QTL genes function, and genomics-assisted breeding.

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.

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.

Chemical compositions and cryo-adhesive probing of the epicuticular wax crystals on fruit surface of wax gourd (Benincasa hispida)

Surface wax crystals play important roles in protecting plants from pest and disease invasions, and UV irradiation. The wax crystals are less probed individually from the fruit surfaces. Herein the morphologies, chemicals and an efficient method to sample the wax blooms of white wax gourd were addressed. Various crystalloids such as rodlets, platelets, fragments, and granules were observed, which stacked as fine wax film covering on wax gourd fruit surface. The wax blooms were effectively removed by cryo-adhesive after consecutive manipulating set by a high-end device with cylinders. Wax crystals were dominated by triterpenols and triterpenol acetates over 61 % of total crystals, followed by vey-long-chain aliphatics. Accordingly, the high-end device with cryo-adhesive provides an efficient approach to selectively probe the wax crystals from those fruits covering wax blooms. The elucidation of morphologies and chemical compositions of wax crystals may help to better understand their regulations on fruit quality traits.

Identification and functional characterization of conserved cis-regulatory elements responsible for early fruit development in cucurbit crops

A number of cis-regulatory elements (CREs) conserved during evolution have been found to be responsible for phenotypic novelty and variation. Cucurbit crops such as cucumber (Cucumis sativus), watermelon (Citrullus lanatus), melon (Cucumis melo), and squash (Cucurbita maxima) develop fruits from an inferior ovary and share some similar biological processes during fruit development. Whether conserved regulatory sequences play critical roles in fruit development of cucurbit crops remains to be explored. In six well-studied cucurbit species, we identified 392,438 conserved noncoding sequences (CNSs), including 82,756 that are specific to cucurbits, by comparative genomics. Genome-wide profiling of accessible chromatin regions (ACRs) and gene expression patterns mapped 20,865 to 43,204 ACRs and their potential target genes for two fruit tissues at two key developmental stages in six cucurbits. Integrated analysis of CNSs and ACRs revealed 4,431 syntenic orthologous CNSs, including 1,687 cucurbit-specific CNSs that overlap with ACRs that are present in all six cucurbit crops and that may regulate the expression of 757 adjacent orthologous genes. CRISPR mutations targeting two CNSs present in the 1,687 cucurbit-specific sequences resulted in substantially altered fruit shape and gene expression patterns of adjacent NAC1 (NAM, ATAF1/2, and CUC2) and EXT-like (EXTENSIN-like) genes, validating the regulatory roles of these CNSs in fruit development. These results not only provide a number of target CREs for cucurbit crop improvement, but also provide insight into the roles of CREs in plant biology and during evolution.

A nonsynonymous mutation in BhLS, encoding an acyl-CoA N-acyltransferase leads to fruit and seed size variation in wax gourd (Benincasa hispida)

Wax gourd (Benincasa hispida (Thunb.) Cogn., 2n = 2x = 24) is an economically important vegetable crop cultivated widely in many tropical and subtropical regions, including China, India, and Japan. Both fruit and seeds are prized agronomic attributes in wax gourd breeding and production. However, the genetic mechanisms underlying these traits remain largely unexplored. In this study, we observed a strong correlation between fruit size and seed size variation in our mapping population, indicating genetic control by a single gene, BhLS, with large size being dominant over small. Through bulk segregant analysis sequencing and fine mapping with a large F2 population, we precisely located the BhLS gene within a 47.098-kb physical interval on Chromosome 10. Within this interval, only one gene, Bhi10M000649, was identified, showing homology to Arabidopsis HOOKLESS1. A nonsynonymous mutation (G to C) in the second exon of Bhi10M000649 was found to be significantly associated with both fruit and seed size variation in wax gourd. These findings collectively highlight the pleiotropic effect of the BhLS gene in regulating fruit and seed size in wax gourd. Our results offer molecular insights into the variation of fruit and seed size in wax gourd and establish a fundamental framework for breeding wax gourd cultivars with desired traits.

Karyotype and LTR-RTs analysis provide insights into oak genomic evolution

BACKGROUND

Whole-genome duplication and long terminal repeat retrotransposons (LTR-RTs) amplification in organisms are essential factors that affect speciation, local adaptation, and diversification of organisms. Understanding the karyotype projection and LTR-RTs amplification could contribute to untangling evolutionary history. This study compared the karyotype and LTR-RTs evolution in the genomes of eight oaks, a dominant lineage in Northern Hemisphere forests.

RESULTS

Karyotype projections showed that chromosomal evolution was relatively conservative in oaks, especially on chromosomes 1 and 7. Modern oak chromosomes formed through multiple fusions, fissions, and rearrangements after an ancestral triplication event. Species-specific chromosomal rearrangements revealed fragments preserved through natural selection and adaptive evolution. A total of 441,449 full-length LTR-RTs were identified from eight oak genomes, and the number of LTR-RTs for oaks from section Cyclobalanopsis was larger than in other sections. Recent amplification of the species-specific LTR-RTs lineages resulted in significant variation in the abundance and composition of LTR-RTs among oaks. The LTR-RTs insertion suppresses gene expression, and the suppressed intensity in gene regions was larger than in promoter regions. Some centromere and rearrangement regions indicated high-density peaks of LTR/Copia and LTR/Gypsy. Different centromeric regional repeat units (32, 78, 79 bp) were detected on different Q. glauca chromosomes.

CONCLUSION

Chromosome fusions and arm exchanges contribute to the formation of oak karyotypes. The composition and abundance of LTR-RTs are affected by its recent amplification. LTR-RTs random retrotransposition suppresses gene expression and is enriched in centromere and chromosomal rearrangement regions. This study provides novel insights into the evolutionary history of oak karyotypes and the organization, amplification, and function of LTR-RTs.

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.

Genome-wide identification of bZIP gene family and expression analysis of BhbZIP58 under heat stress in wax gourd

BACKGROUND

The basic leucine zipper (bZIP) transcription factor family is one of the most abundant and evolutionarily conserved gene families in plants. It assumes crucial functions in the life cycle of plants, including pathogen defense, secondary metabolism, stress response, seed maturation, and flower development. Although the genome of wax gourd has been published, little is known about the functions, evolutionary background, and gene expression patterns of the bZIP gene family, which limits its utilization.

RESULTS

A total of 61 bZIP genes (BhbZIPs) were identified from wax gourd (Benincasa hispida) genome and divided into 12 subgroups. Whole-genome duplication (WGD) and dispersed duplication (DSD) were the main driving forces of bZIP gene family expansion in wax gourd, and this family may have undergone intense purifying selection pressure during the evolutionary process. We selected BhbZIP58, only one in the member of subgroup B, to study its expression patterns under different stresses, including heat, salt, drought, cold stress, and ABA treatment. Surprisingly, BhbZIP58 had a dramatic response under heat stress. BhbZIP58 showed the highest expression level in the root compared with leaves, stem, stamen, pistil, and ovary. In addition, BhbZIP58 protein was located in the nucleus and had transcriptional activation activity. Overexpression of BhbZIP58 in Arabidopsis enhanced their heat tolerance.

CONCLUSIONS

In this study, bZIP gene family is systematically bioinformatically in wax gourd for the first time. Particularly, BhbZIP58 may have an important role in heat stress. It will facilitate further research on the bZIP gene family regarding their evolutionary history and biological functions.

Comparative transcriptome provides insights into gene regulation network associated with the resistance to Fusarium wilt in grafted wax gourd Benincasa hispida

Introduction

Wilt is a soil-borne disease caused by Fusarium that has become a serious threat to wax gourd production. Disease-resistant graft wax gourds are an effective treatment for Fusarium wilt. However, there are few reports on the defense mechanism of graft wax gourd against wilt diseases.

Methods

In the present study, disease and growth indices were compared between grafted and original wax gourds after infection with Fusarium. High level of disease resistance was observed in the grafted wax gourd, with a lower disease index and low impacts on growth after infection. RNA-seq was performed to identify the differentially expressed genes (DEGs) between the adjacent treatment time points in the grafted and original wax gourds, respectively. Then a comparative temporal analysis was performed and a total of 1,190 genes were identified to show different gene expression patterns between the two wax gourd groups during Fusarium infection.

Result and discussion

Here, high level of disease resistance was observed in the grafted wax gourd, with a lower disease index and low impacts on growth after infection. The DEG number was higher in grafted group than original group, and the enriched functional categories and pathways of DEGs were largely inconsistent between the two groups. These genes were enriched in multiple pathways, of which the mitogen-activated protein kinase (MAPK) signaling pathway enhanced the early defense response, and cutin, suberin, and wax biosynthesis signaling pathways enhanced surface resistance in grafted wax gourd in comparison to original group. Our study provides insights into the gene regulatory mechanisms underlying the resistance of grafted wax gourds to Fusarium wilt infection, which will facilitate the breeding and production of wilt-resistant rootstocks.

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

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

Haplotype-resolved genome assembly of Coriaria nepalensis a non-legume nitrogen-fixing shrub

Coriaria nepalensis Wall. (Coriariaceae) is a nitrogen-fixing shrub which forms root nodules with the actinomycete Frankia. Oils and extracts of C. nepalensis have been reported to be bacteriostatic and insecticidal, and C. nepalensis bark provides a valuable tannin resource. Here, by combining PacBio HiFi sequencing and Hi-C scaffolding techniques, we generated a haplotype-resolved chromosome-scale genome assembly for C. nepalensis. This genome assembly is approximately 620 Mb in size with a contig N50 of 11 Mb, with 99.9% of the total assembled sequences anchored to 40 pseudochromosomes. We predicted 60,862 protein-coding genes of which 99.5% were annotated from databases. We further identified 939 tRNAs, 7,297 rRNAs, and 982 ncRNAs. The chromosome-scale genome of C. nepalensis is expected to be a significant resource for understanding the genetic basis of root nodulation with Frankia, toxicity, and tannin biosynthesis.

Lateral Root Initiation in Cucumber (Cucumis sativus): What Does the Expression Pattern of Rapid Alkalinization Factor 34 (RALF34) Tell Us?

In Arabidopsis, the small signaling peptide (peptide hormone) RALF34 is involved in the gene regulatory network of lateral root initiation. In this study, we aimed to understand the nature of the signals induced by RALF34 in the non-model plant cucumber (Cucumis sativus), where lateral root primordia are induced in the apical meristem of the parental root. The RALF family members of cucumber were identified using phylogenetic analysis. The sequence of events involved in the initiation and development of lateral root primordia in cucumber was examined in detail. To elucidate the role of the small signaling peptide CsRALF34 and its receptor CsTHESEUS1 in the initial stages of lateral root formation in the parental root meristem in cucumber, we studied the expression patterns of both genes, as well as the localization and transport of the CsRALF34 peptide. CsRALF34 is expressed in all plant organs. CsRALF34 seems to differ from AtRALF34 in that its expression is not regulated by auxin. The expression of AtRALF34, as well as CsRALF34, is regulated in part by ethylene. CsTHESEUS1 is expressed constitutively in cucumber root tissues. Our data suggest that CsRALF34 acts in a non-cell-autonomous manner and is not involved in lateral root initiation in cucumber.

Identification of clade-wide putative cis-regulatory elements from conserved non-coding sequences in Cucurbitaceae genomes

Cis-regulatory elements regulate gene expression and play an essential role in the development and physiology of organisms. Many conserved non-coding sequences (CNSs) function as cis-regulatory elements. They control the development of various lineages. However, predicting clade-wide cis-regulatory elements across several closely related species remains challenging. Based on the relationship between CNSs and cis-regulatory elements, we present a computational approach that predicts the clade-wide putative cis-regulatory elements in 12 Cucurbitaceae genomes. Using 12-way whole-genome alignment, we first obtained 632 112 CNSs in Cucurbitaceae. Next, we identified 16 552 Cucurbitaceae-wide cis-regulatory elements based on collinearity among all 12 Cucurbitaceae plants. Furthermore, we predicted 3 271 potential regulatory pairs in the cucumber genome, of which 98 were verified using integrative RNA sequencing and ChIP sequencing datasets from samples collected during various fruit development stages. The CNSs, Cucurbitaceae-wide cis-regulatory elements, and their target genes are accessible at http://cmb.bnu.edu.cn/cisRCNEs_cucurbit/. These elements are valuable resources for functionally annotating CNSs and their regulatory roles in Cucurbitaceae genomes.

A comprehensive and conceptual overview of omics-based approaches for enhancing the resilience of vegetable crops against abiotic stresses

Abiotic stresses adversely affect the productivity and production of vegetable crops. The increasing number of crop genomes that have been sequenced or re-sequenced provides a set of computationally anticipated abiotic stress-related responsive genes on which further research may be focused. Knowledge of omics approaches and other advanced molecular tools have all been employed to understand the complex biology of these abiotic stresses. A vegetable can be defined as any component of a plant that is eaten for food. These plant parts may be celery stems, spinach leaves, radish roots, potato tubers, garlic bulbs, immature cauliflower flowers, cucumber fruits, and pea seeds. Abiotic stresses, such as deficient or excessive water, high temperature, cold, salinity, oxidative, heavy metals, and osmotic stress, are responsible for the adverse activity in plants and, ultimately major concern for decreasing yield in many vegetable crops. At the morphological level, altered leaf, shoot and root growth, altered life cycle duration and fewer or smaller organs can be observed. Likewise different physiological and biochemical/molecular processes are also affected in response to these abiotic stresses. In order to adapt and survive in a variety of stressful situations, plants have evolved physiological, biochemical, and molecular response mechanisms. A comprehensive understanding of the vegetable's response to different abiotic stresses and the identification of tolerant genotypes are essential to strengthening each vegetable's breeding program. The advances in genomics and next-generation sequencing have enabled the sequencing of many plant genomes over the last twenty years. A combination of modern genomics (MAS, GWAS, genomic selection, transgenic breeding, and gene editing), transcriptomics, and proteomics along with next-generation sequencing provides an array of new powerful approaches to the study of vegetable crops. This review examines the overall impact of major abiotic stresses on vegetables, adaptive mechanisms and functional genomic, transcriptomic, and proteomic processes used by researchers to minimize these challenges. The current status of genomics technologies for developing adaptable vegetable cultivars that will perform better in future climates is also examined.

A reference-grade genome assembly for Astragalus mongholicus and insights into the biosynthesis and high accumulation of triterpenoids and flavonoids in its roots

Astragalus membranaceus var. mongholicus (AMM), a member of the Leguminosae, is one of the most important medicinal plants worldwide. The dried roots of AMM have a wide range of pharmacological effects and are a traditional Chinese medicine. Here, we report the first chromosome-level reference genome of AMM, comprising nine pseudochromosomes with a total size of 1.47 Gb and 27 868 protein-encoding genes. Comparative genomic analysis reveals that AMM has not experienced an independent whole-genome duplication (WGD) event after the WGD event shared by the Papilionoideae species. Analysis of long terminal repeat retrotransposons suggests a recent burst of these elements at approximately 0.13 million years ago, which may explain the large size of the AMM genome. Multiple gene families involved in the biosynthesis of triterpenoids and flavonoids were expanded, and our data indicate that tandem duplication has been the main driver for expansion of these families. Among the expanded families, the phenylalanine ammonia-lyase gene family was primarily expressed in the roots of AMM, suggesting their roles in the biosynthesis of phenylpropanoid compounds. The functional versatility of 2,3-oxidosqualene cyclase genes in cluster III may play a critical role in the diversification of triterpenoids in AMM. Our findings provide novel insights into triterpenoid and flavonoid biosynthesis and can facilitate future research on the genetics and medical applications of AMM.

Comparative karyotype analysis of eight Cucurbitaceae crops using fluorochrome banding and 45S rDNA-FISH

To have an insight into the karyotype variation of eight Cucurbitaceae crops including Cucumissativus Linnaeus, 1753, Cucumismelo Linnaeus, 1753, Citrulluslanatus (Thunberg, 1794) Matsumura et Nakai, 1916, Benincasahispida (Thunberg, 1784) Cogniaux, 1881, Momordicacharantia Linnaeus, 1753, Luffacylindrica (Linnaeus, 1753) Roemer, 1846, Lagenariasicerariavar.hispida (Thunberg, 1783) Hara, 1948 and Cucurbitamoschata Duchesne ex Poiret, 1819, well morphologically differentiated mitotic metaphase chromosomes were prepared using the enzymatic maceration and flame-drying method, and the chromosomal distribution of heterochromatin and 18S-5.8S-26S rRNA genes (45S rDNA) was investigated using sequential combined PI and DAPI (CPD) staining and fluorescence in situ hybridization (FISH) with 45S rDNA probe. Detailed karyotypes were established using the dataset of chromosome measurements, fluorochrome bands and rDNA FISH signals. Four karyotype asymmetry indices, CVCI, CVCL, MCA and Stebbins' category, were measured to elucidate the karyological relationships among species. All the species studied had symmetrical karyotypes composed of metacentric and submetacentric or only metacentric chromosomes, but their karyotype structure can be discriminated by the scatter plot of MCA vs. CVCL. The karyological relationships among these species revealed by PCoA based on x, 2n, TCL, MCA, CVCL and CVCI was basically in agreement with the phylogenetic relationships revealed by DNA sequences. CPD staining revealed all 45S rDNA sites in all species, (peri)centromeric GC-rich heterochromatin in C.sativus, C.melo, C.lanatus, M.charantia and L.cylindrica, terminal GC-rich heterochromatin in C.sativus. DAPI counterstaining after FISH revealed pericentromeric DAPI+ heterochromatin in C.moschata. rDNA FISH detected two 45S loci in five species and five 45S loci in three species. Among these 45S loci, most were located at the terminals of chromosome arms, and a few in the proximal regions. In C.sativus, individual chromosomes can be precisely distinguished by the CPD band and 45S rDNA signal patterns, providing an easy method for chromosome identification of cucumber. The genome differentiation among these species was discussed in terms of genome size, heterochromatin, 45S rDNA site, and karyotype asymmetry based on the data of this study and previous reports.

A chromosome-level reference genome of the wax gourd (Benincasa hispida)

The wax gourd (Benincasa hispida), the only species in the genus Benincasa, is an important crop native to Asia that has been widely planted for multi-purpose uses. The first wax gourd draft genome was published three years ago, but it was incomplete and highly-fragmented due to data and technical limitations. Herein, we report a new chromosome-level genome assembly and annotation of B. hispida. We generated 974.87 Mb of unitigs with N50 size of 2.43 Mb via a hybrid assembly strategy by using PacBio long reads and Illumina short reads. We then joined them into scaffolds with Hi-C data, resulting 1862 scaffolds with a total length of 975.62 Mb, and 94.92% of the length (926.05 Mb) is contained in the 12 largest scaffolds corresponding to the 12 chromosomes of B. hispida. We predicted 37,092 protein-coding genes, and 85.05% of them were functionally annotated. This chromosome-level reference genome provides significant improvement to the earlier version of draft genome and would be valuable resource for research and molecular breeding of the wax gourd.

Vegetable biology and breeding in the genomics era

Vegetable crops provide a rich source of essential nutrients for humanity and represent critical economic values to global rural societies. However, genetic studies of vegetable crops have lagged behind major food crops, such as rice, wheat and maize, thereby limiting the application of molecular breeding. In the past decades, genome sequencing technologies have been increasingly applied in genetic studies and breeding of vegetables. In this review, we recapitulate recent progress on reference genome construction, population genomics and the exploitation of multi-omics datasets in vegetable crops. These advances have enabled an in-depth understanding of their domestication and evolution, and facilitated the genetic dissection of numerous agronomic traits, which jointly expedites the exploitation of state-of-the-art biotechnologies in vegetable breeding. We further provide perspectives of further directions for vegetable genomics and indicate how the ever-increasing omics data could accelerate genetic, biological studies and breeding in vegetable crops.

Chromosome-Level Genome Assembly of Herpetospermum pedunculosum (Cucurbitaceae)

This study presents a chromosome-level reference genome assembly of a traditional Tibetan medicinal plant, Herpetospermum pedunculosum belonging to the Cucurbitaceae family. Following a combined PacBio high-fidelity sequencing and Hi-C analysis, a final H. pedunculosum genome assembly, 804.11 Mb in length was obtained, 90.45% of which was anchored into ten pseudochromosomes with a contig N50 of 24.39 Mb. In addition, 579.55 Mb repetitive sequences and 23,924 high-confidence protein-coding genes were annotated. Phylogenetic analysis revealed that H. pedunculosum was sister to a clade formed by cucumber, zucchini, and wax gourd. Further whole-genome duplication analysis revealed no recent polyploidization event in the H. pedunculosum genome. The high-quality H. pedunculosum genome presented here will be highly useful in investigating the molecular mechanisms underlying the biosynthesis of its active compounds and adaptation strategies to the extreme environment. It will also provide great insights into comparative genomic studies of Cucurbitaceae and flowering plants.

Telomere-to-telomere genome assembly of bitter melon (Momordica charantia L. var. abbreviata Ser.) reveals fruit development, composition and ripening genetic characteristics

Momordica charantia L. var. abbreviata Ser. (Mca), known as bitter gourd or bitter melon, is a Momordica variety with medicinal value and belongs to the Cucurbitaceae family. In view of the lack of genomic information on bitter gourd and other Momordica species and to promote Mca genomic research, we assembled a 295.6-Mb telomere-to-telomere (T2T) high-quality Mca genome with six gap-free chromosomes after Hi-C correction. This genome is anchored to 11 chromosomes, which is consistent with the karyotype information, and comprises 98 contigs (N50 of 25.4 Mb) and 95 scaffolds (N50 of 25.4 Mb). The Mca genome harbors 19 895 protein-coding genes, of which 45.59% constitute predicted repeat sequences. Synteny analysis revealed variations involved in fruit quality during the divergence of bitter gourd. In addition, assay for transposase-accessible chromatin by high-throughput sequencing and metabolic analysis showed that momordicosides and other substances are characteristic of Mca fruit pulp. A combined transcriptomic and metabolomic analysis revealed the mechanisms of pigment accumulation and cucurbitacin biosynthesis in Mca fruit peels, providing fundamental molecular information for further research on Mca fruit ripening. This report provides a new genetic resource for Momordica genomic studies and contributes additional insights into Cucurbitaceae phylogeny.

Fine mapping and candidate gene analysis of gynoecy trait in chieh-qua (Benincasa hispida Cogn. var. chieh-qua How)

Gynoecy demonstrates an earlier production of hybrids and a higher yield and improves the efficiency of hybrid seed production. Therefore, the utilization of gynoecy is beneficial for the genetic breeding of chieh-qua. However, little knowledge of gynoecious-related genes in chieh-qua has been reported until now. Here, we used an F₂ population from the cross between the gynoecious line 'A36' and the monoecious line 'SX' for genetic mapping and revealed that chieh-qua gynoecy was regulated by a single recessive gene. We fine-mapped it into a 530-kb region flanked by the markers Indel-3 and KASP145 on Chr.8, which harbors eight candidate genes. One of the candidate genes, Bhi08G000345, encoding networked protein 4 (CqNET4), contained a non-synonymous SNP resulting in the amino acid substitution of isoleucine (ATA; I) to methionine (ATG; M). CqNET4 was prominently expressed in the female flower, and only three genes related to ethylene synthesis were significantly expressed between 'A36' and 'SX.' The results presented here provide support for the CqNET4 as the most likely candidate gene for chieh-qua gynoecy, which differed from the reported gynoecious genes.

The Current Developments in Medicinal Plant Genomics Enabled the Diversification of Secondary Metabolites' Biosynthesis

Medicinal plants produce important substrates for their adaptation and defenses against environmental factors and, at the same time, are used for traditional medicine and industrial additives. Plants have relatively little in the way of secondary metabolites via biosynthesis. Recently, the whole-genome sequencing of medicinal plants and the identification of secondary metabolite production were revolutionized by the rapid development and cheap cost of sequencing technology. Advances in functional genomics, such as transcriptomics, proteomics, and metabolomics, pave the way for discoveries in secondary metabolites and related key genes. The multi-omics approaches can offer tremendous insight into the variety, distribution, and development of biosynthetic gene clusters (BGCs). Although many reviews have reported on the plant and medicinal plant genome, chemistry, and pharmacology, there is no review giving a comprehensive report about the medicinal plant genome and multi-omics approaches to study the biosynthesis pathway of secondary metabolites. Here, we introduce the medicinal plant genome and the application of multi-omics tools for identifying genes related to the biosynthesis pathway of secondary metabolites. Moreover, we explore comparative genomics and polyploidy for gene family analysis in medicinal plants. This study promotes medicinal plant genomics, which contributes to the biosynthesis and screening of plant substrates and plant-based drugs and prompts the research efficiency of traditional medicine.

Comparative cytogenomics reveals genome reshuffling and centromere repositioning in the legume tribe Phaseoleae

The tribe Phaseoleae includes several legume crops with assembled genomes. Comparative genomic studies have evidenced the preservation of large genomic blocks among legumes, although chromosome dynamics during Phaseoleae evolution has not been investigated. We conducted a comparative genomic analysis to define an informative genomic block (GB) system and to reconstruct the ancestral Phaseoleae karyotype (APK). We identified GBs based on the orthologous genes between Phaseolus vulgaris and Vigna unguiculata and searched for GBs in different genomes of the Phaseolinae (P. lunatus) and Glycininae (Amphicarpaea edgeworthii) subtribes and Spatholobus suberectus (sister to Phaseolinae and Glycininae), using Medicago truncatula as the outgroup. We also used oligo-FISH probes of two P. vulgaris chromosomes to paint the orthologous chromosomes of two non-sequenced Phaseolinae species. We inferred the APK as having n = 11 and 19 GBs (A to S), hypothesizing five chromosome fusions that reduced the ancestral legume karyotype to n = 11. We identified the rearrangements among the APK and the subtribes and species, with extensive centromere repositioning in Phaseolus. We also reconstructed the chromosome number reduction in S. suberectus. The development of the GB system and the proposed APK provide useful approaches for future comparative genomic analyses of legume species.

A common whole-genome paleotetraploidization in Cucurbitales

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

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

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

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

CuGenDBv2: an updated database for cucurbit genomics

The Cucurbitaceae (cucurbit) family consists of about 1,000 species in 95 genera, including many economically important and popular fruit and vegetable crops. During the past several years, reference genomes have been generated for >20 cucurbit species, and variome and transcriptome profiling data have been rapidly accumulated for cucurbits. To efficiently mine, analyze and disseminate these large-scale datasets, we have developed an updated version of Cucurbit Genomics Database. The updated database, CuGenDBv2 (http://cucurbitgenomics.org/v2), currently hosts 34 reference genomes from 27 cucurbit species/subspecies belonging to 10 different genera. Protein-coding genes from these genomes have been comprehensively annotated by comparing their protein sequences to various public protein and domain databases. A novel 'Genotype' module has been implemented to facilitate mining and analysis of the functionally annotated variome data including SNPs and small indels from large-scale genome sequencing projects. An updated 'Expression' module has been developed to provide a comprehensive gene expression atlas for cucurbits. Furthermore, synteny blocks between any two and within each of the 34 genomes, representing a total of 595 pair-wise genome comparisons, have been identified and can be explored and visualized in the database.

Diploid chromosome-level reference genome and population genomic analyses provide insights into Gypenoside biosynthesis and demographic evolution of Gynostemma pentaphyllum (Cucurbitaceae)

Gynostemma pentaphyllum (Thunb.) Makino is a perennial creeping herbaceous plant in the family Cucurbitaceae, which has great medicinal value and commercial potential, but urgent conservation efforts are needed due to the gradual decreases and fragmented distribution of its wild populations. Here, we report the high-quality diploid chromosome-level genome of G. pentaphyllum obtained using a combination of next-generation sequencing short reads, Nanopore long reads, and Hi-C sequencing technologies. The genome is anchored to 11 pseudo-chromosomes with a total size of 608.95 Mb and 26 588 predicted genes. Comparative genomic analyses indicate that G. pentaphyllum is estimated to have diverged from Momordica charantia 60.7 million years ago, with no recent whole-genome duplication event. Genomic population analyses based on genotyping-by-sequencing and ecological niche analyses indicated low genetic diversity but a strong population structure within the species, which could classify 32 G. pentaphyllum populations into three geographical groups shaped jointly by geographic and climate factors. Furthermore, comparative transcriptome analyses showed that the genes encoding enzyme involved in gypenoside biosynthesis had higher expression levels in the leaves and tendrils. Overall, the findings obtained in this study provide an effective molecular basis for further studies of demographic genetics, ecological adaption, and systematic evolution in Cucurbitaceae species, as well as contributing to molecular breeding, and the biosynthesis and biotransformation of gypenoside.

Natural products of medicinal plants: biosynthesis and bioengineering in post-genomic era

Globally, medicinal plant natural products (PNPs) are a major source of substances used in traditional and modern medicine. As we human race face the tremendous public health challenge posed by emerging infectious diseases, antibiotic resistance and surging drug prices etc., harnessing the healing power of medicinal plants gifted from mother nature is more urgent than ever in helping us survive future challenge in a sustainable way. PNP research efforts in the pre-genomic era focus on discovering bioactive molecules with pharmaceutical activities, and identifying individual genes responsible for biosynthesis. Critically, systemic biological, multi- and inter-disciplinary approaches integrating and interrogating all accessible data from genomics, metabolomics, structural biology, and chemical informatics are necessary to accelerate the full characterization of biosynthetic and regulatory circuitry for producing PNPs in medicinal plants. In this review, we attempt to provide a brief update on the current research of PNPs in medicinal plants by focusing on how different state-of-the-art biotechnologies facilitate their discovery, the molecular basis of their biosynthesis, as well as synthetic biology. Finally, we humbly provide a foresight of the research trend for understanding the biology of medicinal plants in the coming decades.

Integrated analysis of multi-omics and fine-mapping reveals a candidate gene regulating pericarp color and flavonoids accumulation in wax gourd (Benincasa hispida)

Wax gourd (Benincasa hispida), a popular fruit of the Cucurbitaceae (cucurbits) family, contains many nutrients with health benefits and is widely grown in China and other tropical areas. In this study, a wax gourd mutant hfc12 with light-color pericarp was obtained through ethane methylsulfonate (EMS) mutagenesis. Integrative analysis of the metabolome and transcriptome identified 31 differentially accumulated flavonoids (DAFs; flavonoids or flavonoid glycosides) and 828 differentially expressed genes (DEGs) between the hfc12 mutant and wild-type 'BWT'. Furthermore, BSA-seq and kompetitive allele specific PCR (KASP) analysis suggested that the light-color pericarp and higher flavonoid content was controlled by a single gene BhiPRR6 (Bhi12M000742), a typical two-component system (TCS) pseudo-response regulator (PRR). Genetic analysis detected only one nonsynonymous mutation (C-T) in the second exon region of the BhiPRR6. Weighted correlation network analysis (WGCNA) identified the downstream target genes of BhiPRR6, probably regulated by light and were intermediated in the regulatory enzyme reaction of flavonoid biosynthetic pathway. Thus, these results speculated that the transcription factor BhiPRR6, interacting with multiple genes, regulates the absorption of light signals and thereby changes the pericarp color and synthesis of flavonoids in wax gourd.

A time-course transcriptome analysis of wax gourd fruit development reveals predominant genes regulating taste and nutrition

Wax gourd, which belongs to Cucurbitaceae, is an excellent plant resource with the concomitant function of both medicine and foodstuff. Its unique taste and rich nutrition are deeply accepted by consumers. However, the main flavor and nutrients are still unclear, which restricts the quality breeding process of wax gourd. Here, we discovered that monosaccharides, malic acid and citrulline affect the flavor and nutrition of wax gourd and clarified the dynamic accumulation process of these metabolites. To gain insights into the underlying predominant genes regulating accumulation of these metabolites, we performed a time-course transcriptome analysis using RNA-sequencing analysis and compared the expression of screened genes among twenty-four germplasms with different metabolites levels. In addition, the expression abundance among the homologous genes were also analyzed. Finally, a total of 8 genes related to sugar [AGA2 (Bhi03G001926), SUS (Bhi12G001032)], malic acid [MDH (Bhi12G001426, Bhi01G000427), PEPC (Bhi12G000721, Bhi09G002867), ME (Bhi01G002616)] and citrulline [ASS (Bhi02G000401)], respectively were determined. In summary, understanding the core genes influencing taste or nutrition will provide a theoretical basis for fruit quality improvement in wax gourd.

Genetic mapping and genome-wide association study identify BhYAB4 as the candidate gene regulating seed shape in wax gourd (Benincasa hispida)

Wax gourd is an important vegetable crop of the Cucurbitaceae family. According to the shape and structure of the seed coat, the seeds of the wax gourd can be divided into bilateral and unilateral. Bilateral seeds usually germinate quickly and have a high germination rate than unilateral seeds. Thereby, wax gourd varieties with bilateral seeds are more welcomed by seed companies and growers. However, the genetic basis and molecular mechanism regulating seed shape remain unclear in the wax gourd. In this study, the genetic analysis demonstrated that the seed shape of wax gourd was controlled by a single gene, with bilateral dominant to unilateral. Combined with genetic mapping and genome-wide association study, Bhi04G000544 (BhYAB4), encoding a YABBY transcription factor, was identified as the candidate gene for seed shape determination in the wax gourd. A G/A single nucleotide polymorphism variation of BhYAB4 was detected among different germplasm resources, with BhYAB4^G specifically enriched in bilateral seeds and BhYAB4^A in unilateral seeds. The G to A mutation caused intron retention and premature stop codon of BhYAB4. Expression analysis showed that both BhYAB4^G and BhYAB4^A were highly expressed in seeds, while the nuclear localization of BhYAB4^A protein was disturbed compared with that of BhYAB4^G protein. Finally, a derived cleaved amplified polymorphic sequence marker that could efficiently distinguish between bilateral and unilateral seeds was developed, thereby facilitating the molecular marker-assisted breeding of wax gourd cultivars.

Impact of LTR-Retrotransposons on Genome Structure, Evolution, and Function in Curcurbitaceae Species

Long terminal repeat (LTR)-retrotransposons (LTR-RTs) comprise a major portion of many plant genomes and may exert a profound impact on genome structure, function, and evolution. Although many studies have focused on these elements in an individual species, their dynamics on a family level remains elusive. Here, we investigated the abundance, evolutionary dynamics, and impact on associated genes of LTR-RTs in 16 species in an economically important plant family, Cucurbitaceae. Results showed that full-length LTR-RT numbers and LTR-RT content varied greatly among different species, and they were highly correlated with genome size. Most of the full-length LTR-RTs were amplified after the speciation event, reflecting the ongoing rapid evolution of these genomes. LTR-RTs highly contributed to genome size variation via species-specific distinct proliferations. The Angela and Tekay lineages with a greater evolutionary age were amplified in Trichosanthes anguina, whereas a recent activity burst of Reina and another ancient round of Tekay activity burst were examined in Sechium edule. In addition, Tekay and Retand lineages belonging to the Gypsy superfamily underwent a recent burst in Gynostemma pentaphyllum. Detailed investigation of genes with intronic and promoter LTR-RT insertion showed diverse functions, but the term of metabolism was enriched in most species. Further gene expression analysis in G.pentaphyllum revealed that the LTR-RTs within introns suppress the corresponding gene expression, whereas the LTR-RTs within promoters exert a complex influence on the downstream gene expression, with the main function of promoting gene expression. This study provides novel insights into the organization, evolution, and function of LTR-RTs in Cucurbitaceae genomes.

The slow-evolving Acorus tatarinowii genome sheds light on ancestral monocot evolution

Monocots are one of the most diverse groups of flowering plants, and tracing the evolution of their ancestral genome into modern species is essential for understanding their evolutionary success. Here, we report a high-quality assembly of the Acorus tatarinowii genome, a species that diverged early from all the other monocots. Genome-wide comparisons with a range of representative monocots characterized Acorus as a slowly evolved genome with one whole-genome duplication. Our inference of the ancestral monocot karyotypes provides new insights into the chromosomal evolutionary history assigned to modern species and reveals the probable molecular functions and processes related to the early adaptation of monocots to wetland or aquatic habitats (that is, low levels of inorganic phosphate, parallel leaf venation and ephemeral primary roots). The evolution of ancestral gene order in monocots is constrained by gene structural and functional features. The newly obtained Acorus genome offers crucial evidence for delineating the origin and diversification of monocots, including grasses.

Genome sequence and population genomics provide insights into chromosomal evolution and phytochemical innovation of Hippophae rhamnoides

Plants of the Elaeagnaceae family are widely used to treat various health disorders owing to their natural phytochemicals. Seabuckthorn (Hippophae rhamnoides L.) is an economically and ecologically important species within the family with richness of biologically and pharmacologically active substances. Here, we present a chromosome-level genome assembly of seabuckthorn (http://hipp.shengxin.ren/), the first genome sequence of Elaeagnaceae, which has a total length of 849.04 Mb with scaffold N50 of 69.52 Mb and 30 864 annotated genes. Two sequential tetraploidizations with one occurring ~36-41 million years ago (Mya) and the last ~24-27 Mya were inferred, resulting in expansion of genes related to ascorbate and aldarate metabolism, lipid biosynthesis, and fatty acid elongation. Comparative genomic analysis reconstructed the evolutionary trajectories of the seabuckthorn genome with the predicted ancestral genome of 14 proto-chromosomes. Comparative transcriptomic and metabonomic analyses identified some key genes contributing to high content of polyunsaturated fatty acids and ascorbic acid (AsA). Additionally, we generated and analysed 55 whole-genome sequences of diverse accessions, and identified 9.80 million genetic variants in the seabuckthorn germplasms. Intriguingly, genes in selective sweep regions identified through population genomic analysis appeared to contribute to the richness of AsA and fatty acid in seabuckthorn fruits, among which GalLDH, GMPase and ACC, TER were the potentially major-effect causative genes controlling AsA and fatty acid content of the fruit, respectively. Our research offers novel insights into the molecular basis underlying phytochemical innovation of seabuckthorn, and provides valuable resources for exploring the evolution of the Elaeagnaceae family and molecular breeding.

Genome evolution and diversity of wild and cultivated potatoes

Potato (Solanum tuberosum L.) is the world's most important non-cereal food crop, and the vast majority of commercially grown cultivars are highly heterozygous tetraploids. Advances in diploid hybrid breeding based on true seeds have the potential to revolutionize future potato breeding and production^1-4. So far, relatively few studies have examined the genome evolution and diversity of wild and cultivated landrace potatoes, which limits the application of their diversity in potato breeding. Here we assemble 44 high-quality diploid potato genomes from 24 wild and 20 cultivated accessions that are representative of Solanum section Petota, the tuber-bearing clade, as well as 2 genomes from the neighbouring section, Etuberosum. Extensive discordance of phylogenomic relationships suggests the complexity of potato evolution. We find that the potato genome substantially expanded its repertoire of disease-resistance genes when compared with closely related seed-propagated solanaceous crops, indicative of the effect of tuber-based propagation strategies on the evolution of the potato genome. We discover a transcription factor that determines tuber identity and interacts with the mobile tuberization inductive signal SP6A. We also identify 561,433 high-confidence structural variants and construct a map of large inversions, which provides insights for improving inbred lines and precluding potential linkage drag, as exemplified by a 5.8-Mb inversion that is associated with carotenoid content in tubers. This study will accelerate hybrid potato breeding and enrich our understanding of the evolution and biology of potato as a global staple food crop.

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.

Towards Balanced Fertilizer Management in South China: Enhancing Wax Gourd (Benincasa hispida) Yield and Produce Quality

Balanced fertilizer management promotes plant growth, enhances produce quality, minimizes inputs, and reduces negative environmental impacts. Wax gourd (Benincasa hispida) is an important vegetable crop species in China and in South Asia. Two crop nutrition options, NPK and the natural mineral polyhalite, were tested, separately and combined, with the aim of enhancing wax gourd yield and quality and simultaneously to increase nutrient use efficiency and reducing inputs. The experiments tested the optimization of NPK by reducing the proportion of phosphorus (P), and the effect of enriching the soil with essential macronutrients by the use of the supplementary mineral fertilizer polyhalite containing magnesium (Mg), calcium (Ca) and sulfur (S). Two experiments were carried out in Foshan County, Guangdong, China, in 2018 and 2019. Experiments included four treatments: (1) Conventional NPK (15:15:15); (2) Optimized NPK (16:8:18); (3) Conventional NPK + polyhalite; (4) Optimized NPK + polyhalite. Fertilizers were applied prior to planting. While optimized NPK alone had no effects on fruit yield and quality, supplementary polyhalite resulted in a 10–17% increase in yield and significantly improved produce quality due to increased nutrient uptake from polyhalite, resulting in better foliar biomass. We conclude that the combined crop nutrition options improved yield and quality, enhanced nutrient use efficiency, and reduced risks of nutrient pollution. Inclusion of polyhalite in balanced fertilization practices as a supplementary source of secondary macronutrients seems promising. Nevertheless, plenty of space remains open for further adjustments of NPK application management, focusing on reduced rates, optimized ratio, and accurate timing of application for each nutrient.

Identification, Characterization and Function of Orphan Genes Among the Current Cucurbitaceae Genomes

Orphan genes (OGs) that are missing identifiable homologs in other lineages may potentially make contributions to a variety of biological functions. The Cucurbitaceae family consists of a wide range of fruit crops of worldwide or local economic significance. To date, very few functional mechanisms of OGs in Cucurbitaceae are known. In this study, we systematically identified the OGs of eight Cucurbitaceae species using a comparative genomics approach. The content of OGs varied widely among the eight Cucurbitaceae species, ranging from 1.63% in chayote to 16.55% in wax gourd. Genetic structure analysis showed that OGs have significantly shorter protein lengths and fewer exons in Cucurbitaceae. The subcellular localizations of OGs were basically the same, with only subtle differences. Except for aggregation in some chromosomal regions, the distribution density of OGs was higher near the telomeres and relatively evenly distributed on the chromosomes. Gene expression analysis revealed that OGs had less abundantly and highly tissue-specific expression. Interestingly, the largest proportion of these OGs was significantly more tissue-specific expressed in the flower than in other tissues, and more detectable expression was found in the male flower. Functional prediction of OGs showed that (1) 18 OGs associated with male sterility in watermelon; (2) 182 OGs associated with flower development in cucumber; (3) 51 OGs associated with environmental adaptation in watermelon; (4) 520 OGs may help with the large fruit size in wax gourd. Our results provide the molecular basis and research direction for some important mechanisms in Cucurbitaceae species and domesticated crops.

Multiple Stressors in Vegetable Production: Insights for Trait-Based Crop Improvement in Cucurbits

Vegetable production is a key determinant of contribution from the agricultural sector toward national Gross Domestic Product in a country like India, the second largest producer of fresh vegetables in the world. This calls for a careful scrutiny of the threats to vegetable farming in the event of climate extremes, environmental degradation and incidence of plant pests/diseases. Cucurbits are a vast group of vegetables grown almost throughout the world, which contribute to the daily diet on a global scale. Increasing food supply to cater to the ever-increasing world population, calls for intensive, off-season and year-round cultivation of cucurbits. Current situation predisposes these crops to a multitude of stressors, often simultaneously, under field conditions. This scenario warrants a systematic understanding of the different stress specific traits/mechanisms/pathways and their crosstalk that have been examined in cucurbits and identification of gaps and formulation of perspectives on prospective research directions. The careful dissection of plant responses under specific production environments will help in trait identification for genotype selection, germplasm screens to identify superior donors or for direct genetic manipulation by modern tools for crop improvement. Cucurbits exhibit a wide range of acclimatory responses to both biotic and abiotic stresses, among which a few like morphological characters like waxiness of cuticle; primary and secondary metabolic adjustments; membrane thermostability, osmoregulation and, protein and reactive oxygen species homeostasis and turnover contributing to cellular tolerance, appear to be common and involved in cross talk under combinatorial stress exposures. This is assumed to have profound influence in triggering system level acclimation responses that safeguard growth and metabolism. The possible strategies attempted such as grafting initiatives, molecular breeding, novel genetic manipulation avenues like gene editing and ameliorative stress mitigation approaches, have paved way to unravel the prospects for combined stress tolerance. The advent of next generation sequencing technologies and big data management of the omics output generated have added to the mettle of such emanated concepts and ideas. In this review, we attempt to compile the progress made in deciphering the biotic and abiotic stress responses of cucurbits and their associated traits, both individually and in combination.

The vegetable SNP database: An integrated resource for plant breeders and scientists

Single nucleotide polymorphisms (SNPs) are widely used in genetic research and molecular breeding. To date, the genomes of many vegetable crops have been assembled, and hundreds of core germplasms for each vegetable have been sequenced. However, these data are not currently easily accessible because they are stored on different public databases. Therefore, a vegetable crop SNP database should be developed that hosts SNPs demonstrated to have a high success rate in genotyping for genetic research (herein, "alpha SNPs"). We constructed a database (VegSNPDB, http://www.vegsnpdb.cn/) containing the sequence data of 2032 germplasms from 16 vegetable crop species. VegSNPDB hosts 118,725,944 SNPs of which 4,877,305 were alpha SNPs. SNPs can be searched by chromosome number, position, SNP type, genetic population, or specific individuals, as well as the values of MAF, PIC, and heterozygosity. We hope that VegSNPDB will become an important SNP database for the vegetable research community.

A critical review on cytogenetics of Cucurbitaceae with updates on Indian taxa

The cytogenetic relationships in the species of Cucurbitaceae are becoming immensely important to answer questions pertaining to genome evolution. Here, a simplified and updated data resource on cytogenetics of Cucurbitaceae is presented on the basis of foundational parameters (basic, zygotic and gametic chromosome numbers, ploidy, genome size, karyotype) and molecular cytogenetics. We have revised and collated our own findings on seven agriculturally important Indian cucurbit species in a comparative account with the globally published reports. Chromosome count (of around 19% species) shows nearly three-fold differences while genome size (of nearly 5% species) shows 5.84-fold differences across the species. There is no significant correlation between chromosome numbers and nuclear genome sizes. The possible trend of evolution is discussed here based on molecular cytogenetics data, especially the types and distribution of nucleolus organizer regions (NORs). The review supersedes the scopes of general chromosome databases and invites scopes for continuous updates. The offline resource serves as an exclusive toolkit for research and breeding communities across the globe and also opens scope for future establishment of web-database on Cucurbitaceae cytogenetics.

Genome-Wide Identification of YABBY Gene Family in Cucurbitaceae and Expression Analysis in Cucumber (Cucumis sativus L.)

YABBY transcription factors play important roles in plant growth and development. However, little is known about YABBY genes in Cucurbitaceae. Here, we identified 59 YABBY genes from eight cucurbit species, including cucumber (C. sativus L.), melon (C. melon L.), watermelon (C. lanatus), wax gourd (B. hispida), pumpkin (C. maxima), zucchini (C. pepo L.), silver-seed gourd (C. argyrosperma), and bottle gourd (L. siceraria). The 59 YABBY genes were clustered into five subfamilies wherein the gene structures and motifs are conserved, suggesting similar functions within each subfamily. Different YABBY gene numbers in eight cucurbit species indicated that gene loss or duplication events exist in an evolutionary process across Cucurbitaceae. The cis-acting elements analysis implied that the YABBYs may be involved in plant development, and phytohormone, stress, and light responses. Importantly, YABBY genes exhibited organ-specific patterns in expression in cucumber. Furthermore, a gene CsaV3_6G038650 was constitutively expressed at higher levels at different fruit development stages and might play a crucial role in cucumber fruit development. Collectively, our work will provide a better understanding for further function identifications of YABBY genes in Cucurbitaceae.

Comparative Chloroplast Genome Analysis of Wax Gourd (Benincasa hispida) with Three Benincaseae Species, Revealing Evolutionary Dynamic Patterns and Phylogenetic Implications

Benincasa hispida (wax gourd) is an important Cucurbitaceae crop, with enormous economic and medicinal importance. Here, we report the de novo assembly and annotation of the complete chloroplast genome of wax gourd with 156,758 bp in total. The quadripartite structure of the chloroplast genome comprises a large single-copy (LSC) region with 86,538 bp and a small single-copy (SSC) region with 18,060 bp, separated by a pair of inverted repeats (IRa and IRb) with 26,080 bp each. Comparison analyses among B. hispida and three other species from Benincaseae presented a significant conversion regarding nucleotide content, genome structure, codon usage, synonymous and non-synonymous substitutions, putative RNA editing sites, microsatellites, and oligonucleotide repeats. The LSC and SSC regions were found to be much more varied than the IR regions through a divergent analysis of the species within Benincaseae. Notable IR contractions and expansions were observed, suggesting a difference in genome size, gene duplication and deletion, and the presence of pseudogenes. Intronic gene sequences, such as trnR-UCU–atpA and atpH–atpI, were observed as highly divergent regions. Two types of phylogenetic analysis based on the complete cp genome and 72 genes suggested sister relationships between B. hispida with the Citrullus, Lagenaria, and Cucumis. Variations and consistency with previous studies regarding phylogenetic relationships are discussed. The cp genome of B. hispida provides valuable genetic information for the detection of molecular markers, research on taxonomic discrepancies, and the inference of the phylogenetic relationships of Cucurbitaceae.

Graph-based pan-genome reveals structural and sequence variations related to agronomic traits and domestication in cucumber

Structural variants (SVs) represent a major source of genetic diversity and are related to numerous agronomic traits and evolutionary events; however, their comprehensive identification and characterization in cucumber (Cucumis sativus L.) have been hindered by the lack of a high-quality pan-genome. Here, we report a graph-based cucumber pan-genome by analyzing twelve chromosome-scale genome assemblies. Genotyping of seven large chromosomal rearrangements based on the pan-genome provides useful information for use of wild accessions in breeding and genetic studies. A total of ~4.3 million genetic variants including 56,214 SVs are identified leveraging the chromosome-level assemblies. The pan-genome graph integrating both variant information and reference genome sequences aids the identification of SVs associated with agronomic traits, including warty fruits, flowering times and root growth, and enhances the understanding of cucumber trait evolution. The graph-based cucumber pan-genome and the identified genetic variants provide rich resources for future biological research and genomics-assisted breeding.

Increasing studies have suggested that single reference genome is insufficient to capture all variations in the genome. Here, the authors report a graph-based cucumber pan-genome by analyzing 12 chromosome-scale assemblies and reveal variations associated with agronomic traits and domestication.