The genome and transcriptome analysis of snake gourd provide insights into its evolution and fruit development and ripening

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.

Transcriptomics reveal useful resources for examining fruit development and variation in fruit size in Coccinia grandis

Introduction

The Cucurbitaceae family comprises many agronomically important members, that bear nutritious fruits and vegetables of great economic importance. Coccinia grandis, commonly known as Ivy gourd, belongs to this family and is widely consumed as a vegetable. Members of this family are known to display an impressive range of variation in fruit morphology. Although there have been studies on flower development in Ivy gourd, fruit development remains unexplored in this crop.

Methods

In this study, comparative transcriptomics of two Ivy gourd cultivars namely "Arka Neelachal Kunkhi" (larger fruit size) and "Arka Neelachal Sabuja" (smaller fruit size) differing in their average fruit size was performed. A de novo transcriptome assembly for Ivy gourd was developed by collecting fruits at different stages of development (5, 10, 15, and 20 days after anthesis i.e. DAA) from these two varieties. The transcriptome was analyzed to identify differentially expressed genes, transcription factors, and molecular markers.

Results

The transcriptome of Ivy gourd consisted of 155205 unigenes having an average contig size of 1472bp. Unigenes were annotated on publicly available databases to categorize them into different biological functions. Out of these, 7635 unigenes were classified into 38 transcription factor (TF) families, of which Trihelix TFs were most abundant. A total of 11,165 unigenes were found to be differentially expressed in both the varieties and the in silico expression results were validated through real-time PCR. Also, 98768 simple sequence repeats (SSRs) were identified in the transcriptome of Ivy gourd.

Discussion

This study has identified a number of genes, including transcription factors, that could play a crucial role in the determination of fruit shape and size in Ivy gourd. The presence of polymorphic SSRs indicated a possibility for marker-assisted selection for crop breeding in Ivy gourd. The information obtained can help select candidate genes that may be implicated in regulating fruit development and size in other fruit crops.

Fine mapping of TFL, a major gene regulating fruit length in snake gourd (Trichosanthes anguina L)

Fruit length is a crucial agronomic trait of snake gourd (Trichosanthes anguina L); however, genes associated with fruit length have not been characterised. In this study, F2 snake gourd populations were generated by crossing the inbred lines, S1 and S2 (fruit lengths: 110 and 20 cm, respectively). Subsequently, bulk segregant analysis, sequencing, and fine-mapping were performed on the F2 population to identify target genes. Our findings suggest that the fruit length of snake gourd is regulated by a major-effect regulatory gene. Mining of genes regulating fruit length in snake gourd to provide a basis for subsequent selection and breeding of new varieties. Genotype-phenotype association analysis was performed on the segregating F2 population comprising 6,000 plants; the results indicate that the target gene is located on Chr4 (61,846,126-61,865,087 bp, 18.9-kb interval), which only carries the annotated candidate gene, Tan0010544 (designated TFL). TFL belongs to the MADS-box family, one of the largest transcription factor families. Sequence analysis revealed a non-synonymous mutation of base C to G at position 202 in the coding sequence of TFL, resulting in the substitution of amino acid Gln to Glu at position 68 in the protein sequence. Subsequently, an InDel marker was developed to aid the marker-assisted selection of TFL. The TFL in the expression parents within the same period was analysed using quantitative real-time PCR; the TFL expression was significantly higher in short fruits than long fruits. Therefore, TFL can be a candidate gene for determining the fruit length in snake gourd. Collectively, these findings improve our understanding of the genetic components associated with fruit length in snake gourds, which could aid the development of enhanced breeding strategies for plant species.

Positive selection and relaxed purifying selection contribute to rapid evolution of male-biased genes in a dioecious flowering plant

Sex-biased genes offer insights into the evolution of sexual dimorphism. Sex-biased genes, especially those with male bias, show elevated evolutionary rates of protein sequences driven by positive selection and relaxed purifying selection in animals. Although rapid sequence evolution of sex-biased genes and evolutionary forces have been investigated in animals and brown algae, less is known about evolutionary forces in dioecious angiosperms. In this study, we separately compared the expression of sex-biased genes between female and male floral buds and between female and male flowers at anthesis in dioecious Trichosanthes pilosa (Cucurbitaceae). In floral buds, sex-biased gene expression was pervasive, and had significantly different roles in sexual dimorphism such as physiology. We observed higher rates of sequence evolution for male-biased genes in floral buds compared to female-biased and unbiased genes. Male-biased genes under positive selection were mainly associated with functions to abiotic stress and immune responses, suggesting that high evolutionary rates are driven by adaptive evolution. Additionally, relaxed purifying selection may contribute to accelerated evolution in male-biased genes generated by gene duplication. Our findings, for the first time in angiosperms, suggest evident rapid evolution of male-biased genes, advance our understanding of the patterns and forces driving the evolution of sexual dimorphism in dioecious plants.

Semiochemicals from Trichosanthes anguina (Cucurbitaceae) plants influence behavior in Diaphania indica

BACKGROUND

First to third instars of Diaphania indica (Saunders) (Lepidoptera: Crambidae) feed on the lower surface of leaves, while fourth and fifth instars gregariously consume leaves of Trichosanthes anguina L. After defoliating, the caterpillar also attacks flowers and fruits of the plant and finally, results in loss of crop yield. Therefore, behavioral responses of D. indica adults were investigated to volatiles from undamaged (UD), insect-damaged (ID, plants after feeding by D. indica larvae) and jasmonic acid (JA) treated T. anguina plants.

RESULTS

Females showed attraction to volatiles of UD and ID plants of three T. anguina cultivars [MNSR-1 (MNS), Baruipur Long (BAR) and Polo No. 1 (POLO)] in Y-tube olfactometer bioassays. Females did not show significant negative responses from volatiles of JA treated plants. Females were more attracted to volatiles of ID plants than UD plants. Females showed attraction to volatiles of UD or ID plants compared to JA treated plants. Females were attracted to certain synthetic blends resembling volatiles of insect-damaged MNS, BAR and POLO plants in olfactometer bioassays. Females could not distinguish among these three certain synthetic blends in olfactometer bioassays. A synthetic blend of 3Z-hexen-1-ol, α-pinene, hexyl acetate, benzyl alcohol and 6Z-nonenal at mole ratios of 1.47:1.20:1:1.82:1.21 was prepared at 20 mg/mL dichloromethane and 100 μL when used as lure in funnel traps resulted in the capture of the highest number of D. indica adults in field trails.

CONCLUSION

The earlier five-component chemical lure could be used in traps in an integrated pest management program of the insect pest, D. indica. © 2023 Society of Chemical Industry.

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.

Transcriptomic, metabolomic, and ATAC-seq analysis reveal the regulatory mechanism of senescence of post-harvest tomato fruit

Several physiological changes occur during fruit storage, which include the regulation of genes, metabolisms and transcription factors. In this study, we compared 'JF308' (a normal tomato cultivar) and 'YS006' (a storable tomato cultivar) to determine the difference in accumulated metabolites, gene expression, and accessible chromatin regions through metabolome, transcriptome, and ATAC-seq analysis. A total of 1006 metabolites were identified in two cultivars. During storage time, sugars, alcohols and flavonoids were found to be more abundant in 'YS006' compared to 'JF308' on day 7, 14, and 21, respectively. Differentially expressed genes, which involved in starch and sucrose biosynthesis were observed higher in 'YS006'. 'YS006' had lower expression levels of CesA (cellulose synthase), PL (pectate lyase), EXPA (expansin) and XTH (xyglucan endoglutransglucosylase/hydrolase) than 'JF308'. The results showed that phenylpropanoid pathway, carbohydrate metabolism and cell wall metabolism play important roles in prolonging the shelf life of tomato (Solanum lycopersicum) fruit. The ATAC-seq analysis revealed that the most significantly up-regulated transcription factors during storage were TCP 2,3,4,5, and 24 in 'YS006' compared to 'JF308' on day 21. This information on the molecular regulatory mechanisms and metabolic pathways of post-harvest quality changes in tomato fruit provides a theoretical foundation for slowing post-harvest decay and loss, and has theoretical importance and application value in breeding for longer shelf life cultivars.

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.

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.

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.

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.

Cucurbitaceae genome evolution, gene function and molecular breeding

The Cucurbitaceae is one of the most genetically diverse plant families in the world. Many of them are important vegetables or medicinal plants and are widely distributed worldwide. The rapid development of sequencing technologies and bioinformatic algorithms has enabled the generation of genome sequences of numerous important Cucurbitaceae species. This has greatly facilitated research on gene identification, genome evolution, genetic variation and molecular breeding of cucurbit crops. So far, genome sequences of 18 different cucurbit species belonging to tribes Benincaseae, Cucurbiteae, Sicyoeae, Momordiceae and Siraitieae have been deciphered. This review summarizes the genome sequence information, evolutionary relationship, and functional genes associated with important agronomic traits (e.g., fruit quality). The progress of molecular breeding in cucurbit crops and prospects for future applications of Cucurbitaceae genome information are also discussed.

Different scales of gene duplications occurring at different times have jointly shaped the NBS-LRR genes in Prunus species

In this study, genome-wide identification, phylogenetic relationships, duplication time and selective pressure of the NBS-LRR genes, an important group of plant disease-resistance genes (R genes), were performed to uncover their genetic evolutionary patterns in the six Prunus species. A total of 1946 NBS-LRR genes were identified; specifically, 589, 361, 284, 281, 318, and 113 were identified in Prunus yedoensis, P. domestica, P. avium, P. dulcis, P. persica and P. yedoensis var. nudiflora, respectively. Two NBS-LRR gene subclasses, TIR-NBS-LRR (TNL) and non-TIR-NBS-LRR (non-TNL), were also discovered. In total, 435 TNL and 1511 non-TNL genes were identified and could be classified into 30/55/75 and 103/158/191 multi-gene families, respectively, according to three different criteria. Higher Ks and Ka/Ks values were detected in TNL gene families than in non-TNL gene families. These results indicated that the TNL genes had more members involved in relatively ancient duplications and were affected by stronger selection pressure than the non-TNL genes. In general, the NBS-LRR genes were shaped by species-specific duplications, and lineage-specific duplications occurred at recent and relatively ancient periods among the six Prunus species. Therefore, different duplicated copies of NBS-LRRs can resist specific pathogens and will provide an R-gene library for resistance breeding in Prunus species.

De novo transcriptome analysis and identification of candidate genes associated with triterpenoid biosynthesis in Trichosanthes cucumerina L

De novo transcriptome analysis from callus, leaf, and fruit of Trichosanthes cucumerina L. for the identification of genes associated with triterpenoid biosynthesis, especially bryonolic acid and cucurbitacin B. Trichosanthes cucumerina L. (TC) has been used as a medicinal plant in Thailand with various potential functions. Two major triterpenoids found in this plant, bryonolic acid and cucurbitacin B, are receiving increased attention for their activities. Here, we provide TC transcriptome data to identify genes involved in the triterpenoid biosynthetic pathway through callus, where was previously suggested as a novel source for bryonolic acid production as opposed to leaf and fruit. A de novo assembly of approximately 290-thousand transcripts generated from these tissues led to two putative oxidosqualene cyclases: isomultiflorenol synthase (IMS) and cucurbitadienol synthase (CBS). TcIMS and TcCBS, genes that encode substrates for two characteristic triterpenoids in cucurbitaceous plants, were identified as isomultiflorenol synthase and cucurbitadienol synthase, respectively. These two genes were functionally characterised in mutant yeast Gil77 systems, which led to the productions of isomultiflorenol and cucurbitadienol. Moreover, the callus-specific gene expression profiles were also presented. These obtained information showed candidate cytochrome P450s with predicted full-length sequences, which were most likely associated with triterpenoid biosynthesis, especially bryonolic acid. Our study provides useful information and a valuable reference for the further studies on cucurbitaceous triterpenoids.

Combined genomic, transcriptomic, and metabolomic analyses provide insights into chayote (Sechium edule) evolution and fruit development

Chayote (Sechium edule) is an agricultural crop in the Cucurbitaceae family that is rich in bioactive components. To enhance genetic research on chayote, we used Nanopore third-generation sequencing combined with Hi-C data to assemble a draft chayote genome. A chromosome-level assembly anchored on 14 chromosomes (N50 contig and scaffold sizes of 8.40 and 46.56 Mb, respectively) estimated the genome size as 606.42 Mb, which is large for the Cucurbitaceae, with 65.94% (401.08 Mb) of the genome comprising repetitive sequences; 28,237 protein-coding genes were predicted. Comparative genome analysis indicated that chayote and snake gourd diverged from sponge gourd and that a whole-genome duplication (WGD) event occurred in chayote at 25 ± 4 Mya. Transcriptional and metabolic analysis revealed genes involved in fruit texture, pigment, flavor, flavonoids, antioxidants, and plant hormones during chayote fruit development. The analysis of the genome, transcriptome, and metabolome provides insights into chayote evolution and lays the groundwork for future research on fruit and tuber development and genetic improvements in chayote.

Review on the Development and Applications of Medicinal Plant Genomes

With the development of sequencing technology, the research on medicinal plants is no longer limited to the aspects of chemistry, pharmacology, and pharmacodynamics, but reveals them from the genetic level. As the price of next-generation sequencing technology becomes affordable, and the long-read sequencing technology is established, the medicinal plant genomes with large sizes have been sequenced and assembled more easily. Although the review of plant genomes has been reported several times, there is no review giving a systematic and comprehensive introduction about the development and application of medicinal plant genomes that have been reported until now. Here, we provide a historical perspective on the current situation of genomes in medicinal plant biology, highlight the use of the rapidly developing sequencing technologies, and conduct a comprehensive summary on how the genomes apply to solve the practical problems in medicinal plants, like genomics-assisted herb breeding, evolution history revelation, herbal synthetic biology study, and geoherbal research, which are important for effective utilization, rational use and sustainable protection of medicinal plants.