Comparative transcriptome analysis reveals new molecular pathways for cucumber genes related to sex determination

Revisit and explore the ethylene-independent mechanism of sex expression in cucumber (Cucumis sativus)

KEY MESSAGE

This review provides a thorough and comprehensive perspective on the topic of cucumber sexual expression. Specifically, insights into sex expression mediated by pathways other than ethylene are highlighted. Cucumber (Cucumis sativus L.) is a common and important commercial crop that is cultivated and consumed worldwide. Additionally, this species is commonly used as a model for investigating plant sex expression. Cucumbers exhibit a variety of floral arrangements, comprising male, female, and hermaphroditic (bisexual) flowers. Generally, cucumber plants that produce female flowers are typically preferred due to their significant impact on the overall output. Various environmental conditions, such as temperature, light quality, and photoperiod, have been also shown to influence the sex expression in this species. Multiple lines of evidence indicate that ethylene and its biosynthesis genes are crucial in regulating cucumber sex expression. Gibberellins, another well-known phytohormone, can similarly influence cucumber sex expression via an ethylene-independent route. Further studies employing the next-generation sequencing technology also visualized a deeper slice of the molecular mechanism such as the role of the cell cycle program in the cucumber sex expression. This review aims to provide an overview of the sex expression of cucumber including its underlying molecular mechanism and regulatory aspects based on recent investigations.

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

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

Comparative transcriptome analysis reveals candidate genes related to the sex differentiation of Schisandra chinensis

Schisandra chinensis is a monoecious plant with unisex flowers. The fruit of S. chinensis is of high medical with economic value. The yield of S. chinensis fruit is related to the ratio of its female and male flowers. However, there is little research on its floral development and sex differentiation. To elucidate the possible mechanism for the sex differentiation of S. chinensis, we collected 18 samples of female and male flowers from three developmental stages and performed a comparative RNA-seq analysis aimed at identifying differentially expressed genes (DEGs) that may be related to sex differentiation. The results showed 936, 7179, and 6890 differentially expressed genes between female and male flowers at three developmental stages, respectively, and 466 candidate genes may play roles in sex differentiation. KEGG analysis showed genes involved in the flavonoid biosynthesis pathway and DNA replication pathway were essential for the development of female flowers. 51 MADS-box genes and 10 YABBY genes were identified in S. chinensis. The DEGs analysis indicated that MADS-box and YABBY genes were strongly related to the sex determination of S. chinensis. RT-qPCR confirmed the RNA-seq results of 20 differentially expressed genes, including three male-biased genes and 17 female-biased genes. A possible regulatory model of sex differentiation in S. chinensis was proposed according to our results. This study helps reveal the sex-differentiation mechanism of S. chinensis and lays the foundation for regulating the male-female ratio of S. chinensis in the future.

Understanding Transcription Factors and How They Affect Processes in Cucumber Sex Determination

Plant reproduction is a fundamental process on Earth from the perspective of biodiversity, biomass gain, and crop productivity. It is therefore important to understand the sex determination process, and many researchers are investigating the molecular basis of this phenomenon. However, information on the influence of transcription factors (TFs), genes that encode DNA-binding proteins, on this process is limited, although cucumber is a model plant in this regard. In the present study, based on RNA-seq data for differentially expressed genes (DEGs), we aimed to investigate the regulatory TFs that may influence the metabolic processes in the shoot apex containing the forming flower buds. Therefore, the annotation of the genome of the B10 cucumber line was supplemented with the assigned families of transcription factors. By performing ontology analyses of the DEGs, the processes they participate in were identified, and TFs were located among the results. In addition, TFs that have significantly overrepresented targets among DEGs were detected, and sex-specific interactome network maps were generated, indicating the regulatory TFs based on their effects on DEGs and furthermore, on the processes leading to the formation of different-sex flowers. Among the most overrepresented TF families in the sex comparisons were the NAC, bHLH, MYB, and bZIP families. An interaction network analysis indicated the most abundant families among DEGs' regulatory TFs were MYB, AP2/ERF, NAC, and bZIP, and those with the most significant impact on developmental processes were identified, namely the AP/ERF family, followed by DOF, MYB, MADS, and others. Thus, the networks' central nodes and key regulators were identified with respect to male, female, and hermaphrodite forms. Here, we proposed the first model of the regulatory network of TFs that influences the metabolism of sex development in cucumber. These findings may help us to understand the molecular genetics and functional mechanisms underlying sex determination processes.

Advance in sex differentiation in cucumber

Cucumber belongs to the family Cucurbitaceae (melon genus) and is an annual herbaceous vegetable crop. Cucumber is an important cash crop that is grown all over the world. From morphology to cytology, from canonical genetics to molecular biology, researchers have performed much research on sex differentiation and its regulatory mechanism in cucumber, mainly in terms of cucumber sex determination genes, environmental conditions, and the effects of plant hormones, revealing its genetic basis to improve the number of female flowers in cucumber, thus greatly improving the yield of cucumber. This paper reviews the research progress of sex differentiation in cucumber in recent years, mainly focusing on sex-determining genes, environmental conditions, and the influence of phytohormones in cucumber, and provides a theoretical basis and technical support for the realization of high and stable yield cultivation and molecular breeding of cucumber crop traits.

Identifying Genes Associated with Female Flower Development of Phellodendron amurense Rupr. Using a Transcriptomics Approach

Phellodendron amurense Rupr., a species of Rutaceae, is a nationally protected and valuable medicinal plant. It is generally considered to be dioecious. With the discovery of monoecious P. amurense, the phenomenon that its sex development is regulated by epigenetics has been revealed, but the way epigenetics affects the sex differentiation of P. amurense is still unclear. In this study, we investigated the effect of DNA methylation on the sexual development of P. amurense. The young inflorescences of male plants were treated with the demethylation agent 5-azaC, and the induced female flowers were obtained. The induced female flowers’ morphological functions and transcriptome levels were close to those of normally developed plants. Genes associated with the development of female flowers were studied by comparing the differences in transcriptome levels between the male and female flowers. Referring to sex-related genes reported in other plants, 188 candidate genes related to the development of female flowers were obtained, including sex-regulating genes, genes related to the formation and development of sexual organs, genes related to biochemical pathways, and hormone-related genes. RPP0W, PAL3, MCM2, MCM6, SUP, PIN1, AINTEGUMENTA, AINTEGUMENTA-LIKE6, AGL11, SEUSS, SHI-RELATED SEQUENCE 5, and ESR2 were preliminarily considered the key genes for female flower development. This study has demonstrated that epigenetics was involved in the sex regulation of P. amurense, with DNA methylation as one of its regulatory modes. Moreover, some candidate genes related to the sexual differentiation of P. amurense were obtained with analysis. These results are of great significance for further exploring the mechanism of sex differentiation of P. amurense and studying of sex differentiation of plants.

Molecular mapping of genomic regions and identification of possible candidate genes associated with gynoecious sex expression in bitter gourd

Bitter gourd is an important vegetable crop grown throughout the tropics mainly because of its high nutritional value. Sex expression and identification of gynoecious trait in cucurbitaceous vegetable crops has facilitated the hybrid breeding programme in a great way to improve productivity. In bitter gourd, gynoecious sex expression is poorly reported and detailed molecular pathways involve yet to be studied. The present experiment was conducted to study the inheritance, identify the genomic regions associated with gynoecious sex expression and to reveal possible candidate genes through QTL-seq. Segregation for the gynoecious and monoecious sex forms in the F₂ progenies indicated single recessive gene controlling gynoecious sex expression in the genotype, PVGy-201. Gynoecious parent, PVGy-201, Monoecious parent, Pusa Do Mausami (PDM), and two contrasting bulks were constituted for deep-sequencing. A total of 10.56, 23.11, 15.07, and 19.38 Gb of clean reads from PVGy-201, PDM, gynoecious bulk and monoecious bulks were generated. Based on the ΔSNP index, 1.31 Mb regions on the chromosome 1 was identified to be associated with gynoecious sex expression in bitter gourd. In the QTL region 293,467 PVGy-201 unique variants, including SNPs and indels, were identified. In the identified QTL region, a total of 1019 homozygous variants were identified between PVGy1 and PDM genomes and 71 among them were non-synonymous variants (SNPS and INDELs), out of which 11 variants (7 INDELs, 4 SNPs) were classified as high impact variants with frame shift/stop gain effect. In total twelve genes associated with male and female gametophyte development were identified in the QTL-region. Ethylene-responsive transcription factor 12, Auxin response factor 6, Copper-transporting ATPase RAN1, CBL-interacting serine/threonine-protein kinase 23, ABC transporter C family member 2, DEAD-box ATP-dependent RNA helicase 1 isoform X2, Polygalacturonase QRT3-like isoform X2, Protein CHROMATIN REMODELING 4 were identified with possible role in gynoecious sex expression. Promoter region variation in 8 among the 12 genes indicated their role in determining gynoecious sex expression in bitter gourd genotype, DBGy-1. The findings in the study provides insight about sex expression in bitter gourd and will facilitate fine mapping and more precise identification of candidate genes through their functional validation.

Integrative Analyses of Transcriptomics and Metabolomics in Sex Differentiation of Mulberry Flowers

Sex determination and sex differentiation of plants are important physiological processes of plant development. Mulberry (Morus indica L.) is an important economic tree being cultivated in sericulture countries, and mulberry leaf is commonly used for sericulture. The transcriptomic and metabolomic differences between the staminate flowers (SFs) and pistillate flowers (PFs) of mulberry were investigated by RNA sequencing and ultra-performance liquid chromatography-tandem mass spectrometry (UPLC-MS/MS). Overall, we uncovered 4,230 genes and 209 metabolites are significantly differentially expressed between the SFs and PFs of mulberry. The combined transcriptomic and metabolomic analysis revealed these differentially expressed genes (DEGs) and differentially expressed metabolites (DEMs) are involved in flavonoid biosynthesis, galactose metabolism, plant–pathogen interaction, and starch and sucrose metabolism, and these detected DEGs and DEMs may be associated with sex differentiation of mulberry through the regulation of the enrichment pathways, such as the MAPK pathway, flavonoid biosynthesis, galactose metabolism, plant–pathogen interaction, and starch and sucrose metabolism. This study will provide a rich source for the analysis of the molecular mechanism of mulberry sex differentiation processes.

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

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

Fruit transcriptional profiling of the contrasting genotypes for shelf life reveals the key candidate genes and molecular pathways regulating post-harvest biology in cucumber

Cucumber fruits are perishable in nature and become unfit for market within 2-3 days of harvesting. A natural variant, DC-48 with exceptionally high shelf life was developed and used to dissect the genetic architecture and molecular mechanism for extended shelf life through RNA-seq for first time. A total of 1364 DEGs were identified and cell wall degradation, chlorophyll and ethylene metabolism related genes played key role. Polygalacturunase (PG), Expansin (EXP) and xyloglucan were down regulated determining fruit firmness and retention of fresh green colour was mainly attributed to the low expression level of the chlorophyll catalytic enzymes (CCEs). Gene regulatory networks revealed the hub genes and cross-talk associated with wide variety of the biological processes. Large number of SSRs (21524), SNPs (545173) and InDels (126252) identified will be instrumental in cucumber improvement. A web genomic resource, CsExSLDb developed will provide a platform for future investigation on cucumber post-harvest biology.

Genome-Wide Association Studies for Sex Determination and Cross-Compatibility in Water Yam (Dioscorea alata L.)

Yam (Dioscorea spp.) species are predominantly dioecious, with male and female flowers borne on separate individuals. Cross-pollination is, therefore, essential for gene flow among and within yam species to achieve breeding objectives. Understanding genetic mechanisms underlying sex determination and cross-compatibility is crucial for planning a successful hybridization program. This study used the genome-wide association study (GWAS) approach for identifying genomic regions linked to sex and cross-compatibility in water yam (Dioscorea alata L.). We identified 54 markers linked to flower sex determination, among which 53 markers were on chromosome 6 and one on chromosome 11. Our result ascertained that D. alata is characterized by the male heterogametic sex determination system (XX/XY). The cross-compatibility indices, average crossability rate (ACR) and percentage high crossability (PHC), were controlled by loci on chromosomes 1, 6 and 17. Of the significant loci, SNPs located on chromosomes 1 and 17 were the most promising for ACR and PHC, respectively, and should be validated for use in D. alata hybridization activities to predict cross-compatibility success. A total of 61 putative gene/protein families with direct or indirect influence on plant reproduction were annotated in chromosomic regions controlling the target traits. This study provides valuable insights into the genetic control of D. alata sexual reproduction. It opens an avenue for developing genomic tools for predicting hybridization success in water yam breeding programs.

Using RNA-seq to characterize pollen-stigma interactions for pollination studies

Insects are essential for the reproduction of pollinator-dependent crops and contribute to the pollination of 87% of wild plants and 75% of the world’s food crops. Understanding pollen flow dynamics between plants and pollinators is thus essential to manage and conserve wild plants and ensure yields are maximized in food crops. However, the determination of pollen transfer in the field is complex and laborious. We developed a field experiment in a pollinator-dependent crop and used high throughput RNA sequencing (RNA-seq) to quantify pollen flow by measuring changes in gene expression between pollination treatments across different apple (Malus domestica Borkh.) cultivars. We tested three potential molecular indicators of successful pollination and validated these results with field data by observing single and multiple visits by honey bees (Apis mellifera) to apple flowers and measured fruit set in a commercial apple orchard. The first indicator of successful outcrossing was revealed via differential gene expression in the cross-pollination treatments after 6 h. The second indicator of successful outcrossing was revealed by the expression of specific genes related to pollen tube formation and defense response at three different time intervals in the stigma and the style following cross-pollination (i.e. after 6, 24, and 48 h). Finally, genotyping variants specific to donor pollen could be detected in cross-pollination treatments, providing a third indicator of successful outcrossing. Field data indicated that one or five flower visits by honey bees were insufficient and at least 10 honey bee flower visits were required to achieve a 25% probability of fruit set under orchard conditions. By combining the genotyping data, the differential expression analysis, and the traditional fruit set field experiments, it was possible to evaluate the pollination effectiveness of honey bee visits under orchards conditions. This is the first time that pollen-stigma-style mRNA expression analysis has been conducted after a pollinator visit (honey bee) to a plant (in vivo apple flowers). This study provides evidence that mRNA sequencing can be used to address complex questions related to stigma–pollen interactions over time in pollination ecology.

Comparative RNA-Seq analysis reveals genes associated with masculinization in female Cannabis sativa

Using RNA profiling, we identified several silver thiosulfate-induced genes that potentially control the masculinization of female Cannabis sativa plants. Genetically female Cannabis sativa plants normally bear female flowers, but can develop male flowers in response to environmental and developmental cues. In an attempt to elucidate the molecular elements responsible for sex expression in C. sativa plants, we developed genetically female lines producing both female and chemically-induced male flowers. Furthermore, we carried out RNA-Seq assays aimed at identifying differentially expressed genes responsible for male flower development in female plants. The results revealed over 10,500 differentially expressed genes, of which around 200 potentially control masculinization of female cannabis plants. These genes include transcription factors and other genes involved in male organ (i.e., anther and pollen) development, as well as genes involved in phytohormone signalling and male-biased phenotypes. The expressions of 15 of these genes were further validated by qPCR assay confirming similar expression patterns to that of RNA-Seq data. These genes would be useful for understanding predisposed plants producing flowers of both sex types in the same plant, and help breeders to regulate the masculinization of female plants through targeted breeding and plant biotechnology.

Transcriptional and Hormonal Responses in Ethephon-Induced Promotion of Femaleness in Pumpkin

The number and proportion of female flowers per plant can directly influence the yield and economic benefits of cucurbit crops. Ethephon is often used to induce female flowers in cucurbits. However, the mechanism through which it affects floral sex differentiation in pumpkin is unknown. We found that the application of ethephon on shoot apical meristem of pumpkin at seedling stage significantly increased the number of female flowers and expedited the appearance of the first female flower. These effects were further investigated by transcriptome and hormone analyses of plants sprayed with ethephon. A total of 647 differentially expressed genes (DEGs) were identified, among which 522 were upregulated and 125 were downregulated. Gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) analysis indicated that these genes were mainly enriched in plant hormone signal transduction and 1-aminocyclopropane-1-carboxylate oxidase (ACO). The results suggests that ethylene is a trigger for multiple hormone signaling, with approximately 4.2% of the identified DEGs involved in ethylene synthesis and multiple hormone signaling. Moreover, ethephon significantly reduced the levels of jasmonic acid (JA), jasmonoyl-L-isoleucine (JA-ILE), and para-topolin riboside (pTR) but increased the levels of 3-indoleacetamide (IAM). Although the level of 1-aminocyclopropanecarboxylic acid was not changed, the expression of ACO genes, which code for the enzyme catalyzing the key rate-limiting step in ethylene production, was significantly upregulated after ethephon treatment. The results indicate that the ethephon affects the transcription of ethylene synthesis and signaling genes, and other hormone signaling genes, especially auxin responsive genes, and modulates the levels of auxin, jasmonic acid, and cytokinin (CK), which may together contribute to femaleness.

Gibberellic Acid Induced Changes on Growth, Yield, Superoxide Dismutase, Catalase and Peroxidase in Fruits of Bitter Gourd (Momordica charantia L.)

Bitter gourd is one of the important cucurbits and highly liked among both farmers and consumers due to its high net return and nutritional value. However, being monoecious, it exhibits substantial variation in flower bearing pattern. Plant growth regulators (PGRs) are known to influence crop phenology while gibberellic acid (GA3) is one of the most prominent PGRs that influence cucurbits phenology. Therefore, a field trial was conducted at University of Agriculture Faisalabad to evaluate the impact of a commercial product of gibberellic acid (GA3) on growth, yield and quality attributes of two bitter gourd (Momordica charantiaL.) cultivars. We used five different concentrations (0.4 g, 0.6 g, 0.8 g, 1.0 g, and 1.2 g per litre) of commercial GA3 product (Gibberex, 10% Gibberellic acid). Results showed that a higher concentration of gibberex (1.0 and 1.20 g L−1 water) enhanced the petiole length, intermodal length, and yield of bitter gourd cultivars over control in Golu hybrid and Faisalabad Long. A significant decrease in the enzyme superoxidase dismutase, peroxidase, and catalase activities were observed with an increasing concentration of gibberex (1.0 and 1.20 gL−1 water) as compared to control. These results indicate that the exogenous application of gibberex at a higher concentration (1.2 g L−1) has a dual action in bitter gourd plant: i) it enhances the plant growth and yield, and ii) it also influenced the antioxidant enzyme activities in fruits. These findings may have a meaningful, practical use for farmers involved in agriculture and horticulture.

Transcriptome sequencing and screening of genes related to sex determination of Trichosanthes kirilowii Maxim

Trichosanthes kirilowii Maxim. (TK) is a dioecious plant in the Cucurbitaceae for which different sexes have separate medicinal uses. In order to study the genes related to sex determination, transcriptome sequencing was performed on flower buds of male and female plants using the high-throughput sequencing technology. A total of 145,975 unigenes and 7110 DEGs were obtained. There were 6776 DEGs annotated to 1234 GO terms and enriched to 18 functional groups, including five biological processes related to sugar metabolism. KEGG pathway analysis indicated genes involved in hormone transduction, hormone synthesis and carbohydrate metabolism. Many DEGs of TK are involved in reproductive organ formation, hormone signal transduction and regulatory networks. Combining the results of GO, KEGG and qRT-PCR, 11 sex determining candidate genes of TK were selected, including MYB80, MYB108, CER1, CBL9, ABCB19, SERK1, HSP81-3, ACS9, SEP3, AUX1 and YUC6. The results provide a foundation for the study of sex differentiation in TK.

Gene Interactions Regulating Sex Determination in Cucurbits

The family Cucurbitaceae includes many economically important crops, such as cucumber (Cucumis sativus), melon (Cucumis melo), watermelon (Citrullus lanatus), and zucchini (Cucurbita pepo), which share homologous gene pathways that control similar phenotypes. Sex determination is a research hotspot associated with yield and quality, and the genes involved are highly orthologous and conserved in cucurbits. In the field, six normal sex types have been categorized according to the distribution of female, male, or bisexual flowers in a given plant. To date, five orthologous genes involved in sex determination have been cloned, and their various combinations and expression patterns can explain all the identified sex types. In addition to genetic mechanisms, ethylene controls sex expression in this family. Two ethylene signaling components have been identified recently, which will help us to explore the ethylene signaling-mediated interactions among sex-related genes. This review discusses recent advances relating to the mechanism of sex determination in cucurbits and the prospects for research in this area.