RNA-Seq improves annotation of protein-coding genes in the cucumber genome

A nuclease-dead Cas9-derived tool represses target gene expression

Manipulation of gene expression is central to understanding gene function, engineering cell behavior, and altering biological traits according to production demands. Nuclease-dead Cas9 (dCas9), a variant of active Cas9, offers a versatile platform for the precise control of genome function without DNA cleavage. Notably, however, an effective and universal dCas9-based transcriptional repression system remains unavailable in plants. The noncanonical histone acetyltransferase TENDRIL-LESS (CsTEN) is responsible for chromatin loosening and histone modification in cucumber (Cucumis sativus). In this study, we engineered a gene regulation tool by fusing TEN and its truncated proteins with dCas9. The full-length dCas9-TEN protein substantially repressed gene expression, with the N-terminal domain identified as the core repression domain. We subsequently validated the specificity and efficacy of this system through both transient infection and genetic transformation in cucumber and Arabidopsis (Arabidopsis thaliana). The electrophoretic mobility shift assay (EMSA) revealed the ability of the N-terminal domain of TEN to bind to chromatin, which may promote target binding of the dCas9 complex and enhance the transcriptional repression effect. Our tool enriches the arsenal of genetic regulation tools available for precision breeding in crops.

Updated Gene Prediction of the Cucumber (9930) Genome through Manual Annotation

Thorough and precise gene structure annotations are essential for maximizing the benefits of genomic data and unveiling valuable genetic insights. The cucumber genome was first released in 2009 and updated in 2019. To increase the accuracy of the predicted gene models, 64 published RNA-seq data and 9 new strand-specific RNA-seq data from multiple tissues were used for manual comparison with the gene models. The updated annotation file (V3.1) contains an increased number (24,145) of predicted genes compared to the previous version (24,317 genes), with a higher BUSCO value of 96.9%. A total of 6231 and 1490 transcripts were adjusted and newly added, respectively, accounting for 31.99% of the overall gene tally. These newly added and adjusted genes were renamed (CsaV3.1_XGXXXXX), while genes remaining unaltered preserved their original designations. A random selection of 21 modified/added genes were validated using RT-PCR analyses. Additionally, tissue-specific patterns of gene expression were examined using the newly obtained transcriptome data with the revised gene prediction model. This improved annotation of the cucumber genome will provide essential and accurate resources for studies in cucumber.

Comparative Genomic Analysis of PEBP Genes in Cucurbits Explores the Interactors of Cucumber CsPEBPs Related to Flowering Time

The family of phosphatidylethanolamine-binding proteins (PEBPs) participates in various plant biological processes, mainly flowering regulation and seed germination. In cucurbit crops, several PEBP genes have been recognized to be responsible for flowering time. However, the investigation of PEBP family members across the genomes of cucurbit species has not been reported, and their conservation and divergence in structure and function remain largely unclear. Herein, PEBP genes were identified from seven cucurbit crops and were used to perform a comparative genomics analysis. The cucurbit PEBP proteins could be classified into MFT, FT, TFL, and PEBP clades, and further, the TFL clade was divided into BFT-like, CEN-like, and TFL1-like subclades. The MFT-like, FT-like, and TFL-like proteins were clearly distinguished by a critical amino acid residue at the 85th position of the Arabidopsis FT protein. In gene expression analysis, CsaPEBP1 was highly expressed in flowers, and its expression levels in females and males were 70.5 and 89.2 times higher, respectively, than those in leaves. CsaPEBP5, CsaPEBP6, and CsaPEBP7 were specifically expressed in male flowers, with expression levels 58.1, 17.3, and 15.7 times higher, respectively, than those of leaves. At least five CsaPEBP genes exhibited the highest expression during the later stages of corolla opening. Through clustering of time-series-based RNA-seq data, several potential transcription factors (TFs) interacting with four CsaPEBPs were identified during cucumber corolla opening. Because of the tandem repeats of binding sites in promoters, NF-YB (Csa4G037610) and GATA (Csa7G64580) TFs appeared to be better able to regulate the CsaPEBP2 and CsaPEBP5 genes, respectively. This study would provide helpful information for further investigating the roles of PEBP genes and their interacting TFs in growth and development processes, such as flowering time regulation in cucurbit crops.

Do DEEPER ROOTING 1 Homologs Regulate the Lateral Root Slope Angle in Cucumber (Cucumis sativus)?

The architecture of the root system is fundamental to plant productivity. The rate of root growth, the density of lateral roots, and the spatial structure of lateral and adventitious roots determine the developmental plasticity of the root system in response to changes in environmental conditions. One of the genes involved in the regulation of the slope angle of lateral roots is DEEPER ROOTING 1 (DRO1). Its orthologs and paralogs have been identified in rice, Arabidopsis, and several other species. However, nothing is known about the formation of the slope angle of lateral roots in species with the initiation of lateral root primordia within the parental root meristem. To address this knowledge gap, we identified orthologs and paralogs of the DRO1 gene in cucumber (Cucumis sativus) using a phylogenetic analysis of IGT protein family members. Differences in the transcriptional response of CsDRO1, CsDRO1-LIKE1 (CsDRO1L1), and CsDRO1-LIKE2 (CsDRO1L2) to exogenous auxin were analyzed. The results showed that only CsDRO1L1 is auxin-responsive. An analysis of promoter-reporter fusions demonstrated that the CsDRO1, CsDRO1L1, and CsDRO1L2 genes were expressed in the meristem in cell files of the central cylinder, endodermis, and cortex; the three genes displayed different expression patterns in cucumber roots with only partial overlap. A knockout of individual CsDRO1, CsDRO1L1, and CsDRO1L2 genes was performed via CRISPR/Cas9 gene editing. Our study suggests that the knockout of individual genes does not affect the slope angle formation during lateral root primordia development in the cucumber parental root.

Pan-genome Analysis of WOX Gene Family and Function Exploration of CsWOX9 in Cucumber

Cucumber is an economically important vegetable crop, and the warts (composed of spines and Tubercules) of cucumber fruit are an important quality trait that influences its commercial value. WOX transcription factors are known to have pivotal roles in regulating various aspects of plant growth and development, but their studies in cucumber are limited. Here, genome-wide identification of cucumber WOX genes was performed using the pan-genome analysis of 12 cucumber varieties. Our findings revealed diverse CsWOX genes in different cucumber varieties, with variations observed in protein sequences and lengths, gene structure, and conserved protein domains, possibly resulting from the divergent evolution of CsWOX genes as they adapt to diverse cultivation and environmental conditions. Expression profiles of the CsWOX genes demonstrated that CsWOX9 was significantly expressed in unexpanded ovaries, especially in the epidermis. Additionally, analysis of the CsWOX9 promoter revealed two binding sites for the C2H2 zinc finger protein. We successfully executed a yeast one-hybrid assay (Y1H) and a dual-luciferase (LUC) transaction assay to demonstrate that CsWOX9 can be transcriptionally activated by the C2H2 zinc finger protein Tu, which is crucial for fruit Tubercule formation in cucumber. Overall, our results indicated that CsWOX9 is a key component of the molecular network that regulates wart formation in cucumber fruits, and provide further insight into the function of CsWOX genes in cucumber.

Identification of QTL associated with resistance to Phytophthora fruit rot in cucumber (Cucumis sativus L.)

Phytophthora fruit rot (PFR) caused by the soilborne oomycete pathogen, Phytophthora capsici, can cause severe yield loss in cucumber. With no resistant variety available, genetic resources are needed to develop resistant varieties. The goal of this work was to identify quantitative trait loci (QTL) associated with resistance to PFR using multiple genomic approaches and populations. Two types of resistances have been identified: age-related resistance (ARR) and young fruit resistance. ARR occurs at 12-16 days post pollination (dpp), coinciding with the end of exponential fruit growth. A major QTL for ARR was discovered on chromosome 3 and a candidate gene identified based on comparative transcriptomic analysis. Young fruit resistance, which is observed during the state of rapid fruit growth prior to commercial harvest, is a quantitative trait for which multiple QTL were identified. The largest effect QTL, qPFR5.1, located on chromosome 5 was fine mapped to a 1-Mb region. Genome-wide association studies (GWAS) and extreme-phenotype genome-wide association study (XP-GWAS) for young fruit resistance were also performed on a cucumber core collection representing > 96% of the genetic diversity of the USDA cucumber germplasm. Several SNPs overlapped with the QTL identified from QTL-seq analysis on biparental populations. In addition, novel SNPs associated with the resistance were identified from the germplasm. The resistant alleles were found mostly in accessions from India and South Asia, the center of diversity for cucumber. The results from this work can be applied to future disease resistance studies and marker-assisted selection in breeding programs.

Genome-Wide Identification and Expression Analysis of the Cucumber FKBP Gene Family in Response to Abiotic and Biotic Stresses

The FKBP (FK506-binding protein) gene family is an important member of the PPlase protease family and plays a vital role during the processes of plant growth and development. However, no studies of the FKBP gene family have been reported in cucumber. In this study, 19 FKBP genes were identified in cucumber, which were located on chromosomes 1, 3, 4, 6, and 7. Phylogenetic analysis divided the cucumber FKBP genes into three subgroups. The FKBP genes in the same subgroup exhibited similar structures and conserved motifs. The cis-acting elements analysis revealed that the promoters of cucumber FKBP genes contained hormone-, stress-, and development-related cis-acting elements. Synteny analysis of the FKBP genes among cucumber, Arabidopsis, and rice showed that 12 kinds of syntenic relationships were detected between cucumber and Arabidopsis FKBP genes, and 3 kinds of syntenic relationships were observed between cucumber and rice FKBP genes. The tissue-specific expression analysis showed that some FKBP genes were expressed in all tissues, while others were only highly expressed in part of the 10 types of tissues. The expression profile analysis of cucumber FKBP genes under 13 types of stresses showed that the CsaV3_1G007080 gene was differentially expressed under abiotic stresses (high temperature, NaCl, silicon, and photoperiod) and biotic stresses (downy mildew, green mottle mosaic virus, Fusarium wilt, phytophthora capsica, angular leaf spot, and root-knot nematode), which indicated that the CsaV3_1G007080 gene plays an important role in the growth and development of cucumber. The interaction protein analysis showed that most of the proteins in the FKBP gene family interacted with each other. The results of this study will lay the foundation for further research on the molecular biological functions of the cucumber FKBP gene family.

Genome-wide characterization, phylogenetic and expression analysis of ABCG gene subfamily in cucumber (Cucumis sativus L.)

The ABCG is the largest subfamily of the ABC family with extensive functions, and only a few members have been identified in detail. However, more and more studies have shown that the members of this family are very important and are involved in many life processes such as plant development and response to various stresses. Cucumber is an important vegetable crops around the world. The cucumber development is essential for its production and quality. Meanwhile, various stresses have caused serious losses of cucumber. However, the ABCG genes were not well characterized and functioned in cucumber. In this study, the cucumber CsABCG gene family were identified and characterized, and their evolutionary relationship and functions were analyzed. The cis-acting elements and expression analysis showed that they played important role in development and responding to various biotic and abiotic stresses in cucumber. Phylogenetic analysis, sequence alignment and MEME (Multiple Em for Motif Elicitation) analysis indicated that the functions of ABCG proteins in different plants are evolutionarily conserved. Collinear analysis revealed that the ABCG gene family was highly conserved during the evolution. In addition, the potential binding sites of the CsABCG genes targeted by miRNA were predicted. These results will lay a foundation for further research on the function of the CsABCG genes in cucumber.

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.

Genome-wide identification and expression analysis of BZR gene family and associated responses to abiotic stresses in cucumber (Cucumis sativus L.)

BACKGROUND

BRASSINAZOLE-RESISTANT (BZR) is a class of specific transcription factor (TFs) involved in brassinosteroid (BR) signal transduction. The regulatory mechanism of target genes mediated by BZR has become one of the key research areas in plant BR signaling networks. However, the functions of the BZR gene family in cucumber have not been well characterized.

RESULTS

In this study, six CsBZR gene family members were identified by analyzing the conserved domain of BES1 N in the cucumber genome. The size of CsBZR proteins ranges from 311 to 698 amino acids and are mostly located in the nucleus. Phylogenetic analysis divided CsBZR genes into three subgroups. The gene structure and conserved domain showed that the BZR genes domain in the same group was conserved. Cis-acting element analysis showed that cucumber BZR genes were mainly involved in hormone response, stress response and growth regulation. The qRT-PCR results also confirmed CsBZR response to hormones and abiotic stress.

CONCLUSION

Collectively, the CsBZR gene is involved in regulating cucumber growth and development, particularly in hormone response and response to abiotic stress. These findings provide valuable information for understanding the structure and expression patterns of BZR genes.

Integrative Proteomics and Metabolomics Analysis Reveals the Role of Small Signaling Peptide Rapid Alkalinization Factor 34 (RALF34) in Cucumber Roots

The main role of RALF small signaling peptides was reported to be the alkalization control of the apoplast for improvement of nutrient absorption; however, the exact function of individual RALF peptides such as RALF34 remains unknown. The Arabidopsis RALF34 (AtRALF34) peptide was proposed to be part of the gene regulatory network of lateral root initiation. Cucumber is an excellent model for studying a special form of lateral root initiation taking place in the meristem of the parental root. We attempted to elucidate the role of the regulatory pathway in which RALF34 is a participant using cucumber transgenic hairy roots overexpressing CsRALF34 for comprehensive, integrated metabolomics and proteomics studies, focusing on the analysis of stress response markers. CsRALF34 overexpression resulted in the inhibition of root growth and regulation of cell proliferation, specifically in blocking the G2/M transition in cucumber roots. Based on these results, we propose that CsRALF34 is not part of the gene regulatory networks involved in the early steps of lateral root initiation. Instead, we suggest that CsRALF34 modulates ROS homeostasis and triggers the controlled production of hydroxyl radicals in root cells, possibly associated with intracellular signal transduction. Altogether, our results support the role of RALF peptides as ROS regulators.

Nanomaterial Size and Surface Modification Mediate Disease Resistance Activation in Cucumber (Cucumis sativus)

Crop disease represents a serious and increasing threat to global food security. Lanthanum oxide nanomaterials (La₂O₃ NMs) with different sizes (10 and 20 nm) and surface modifications (citrate, polyvinylpyrrolidone [PVP], and poly(ethylene glycol)) were investigated for their control of the fungal pathogen Fusarium oxysporum (Schl.) f. sp cucumerinum Owen on six-week-old cucumber (Cucumis sativus) in soil. Seed treatment and foliar application of the La₂O₃ NMs at 20-200 mg/kg (mg/L) significantly suppressed cucumber wilt (decreased by 12.50-52.11%), although the disease control efficacy was concentration-, size-, and surface modification-dependent. The best pathogen control was achieved by foliar application of 200 mg/L PVP-coated La₂O₃ NMs (10 nm); disease severity was decreased by 67.6%, and fresh shoot biomass was increased by 49.9% as compared with pathogen-infected control. Importantly, disease control efficacy was 1.97- and 3.61-fold greater than that of La₂O₃ bulk particles and a commercial fungicide (Hymexazol), respectively. Additionally, La₂O₃ NMs application enhanced cucumber yield by 350-461%, increased fruit total amino acids by 295-344%, and improved fruit vitamin content by 65-169% as compared with infected controls. Transcriptomic and metabolomic analyses revealed that La₂O₃ NMs: (1) interacted with calmodulin, subsequently activating salicylic acid-dependent systemic acquired resistance; (2) increased the activity and expression of antioxidant and related genes, thereby alleviating pathogen-induced oxidative stress; and (3) directly inhibited in vivo pathogen growth. The findings highlight the significant potential of La₂O₃ NMs for suppressing plant disease in sustainable agriculture.

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.

Cytological, genetic and transcriptomic characterization of a cucumber albino mutant

Photosynthesis, a fundamental process for plant growth and development, is dependent on chloroplast formation and chlorophyll synthesis. Severe disruption of chloroplast structure results in albinism of higher plants. In the present study, we report a cucumber albino alc mutant that presented white cotyledons under normal light conditions and was unable to produce first true leaf. Meanwhile, alc mutant could grow creamy green cotyledons under dim light conditions but died after exposure to normal light irradiation. No chlorophyll and carotenoid were detected in the alc mutant grown under normal light conditions. Using transmission electron microscopy, impaired chloroplasts were observed in this mutant. The genetic analysis indicated that the albino phenotype was recessively controlled by a single locus. Comparative transcriptomic analysis between the alc mutant and wild type revealed that genes involved in chlorophyll metabolism and the methylerythritol 4-phosphate pathway were affected in the alc mutant. In addition, three genes involved in chloroplast development, including two FtsH genes and one PPR gene, were found to have negligible expression in this mutant. The quality of RNA sequencing results was further confirmed by real-time quantitative PCR analysis. We also examined 12 homologous genes from alc mutant in other plant species, but no genetic variation in the coding sequences of these genes was found between alc mutant and wild type. Taken together, we characterized a cucumber albino mutant with albinism phenotype caused by chloroplast development deficiency and this mutant can pave way for future studies on plastid development.

Identification and Expression Profile Analysis of the OSCA Gene Family Related to Abiotic and Biotic Stress Response in Cucumber

Simple Summary

Hyperosmolality-gated calcium-permeable channels (OSCAs) are calcium nonselective cation channel proteins involved in multiple biological processes. In this work, the members of the OSCA family in cucumber were systematically analyzed, including their sequence characteristics, phylogenetic relationships, conserved motifs, gene structures, promoter regions, and tissue expression patterns. In addition, the effects of different osmotic-related abiotic stresses [salt (NaCl), drought (PEG), and abscisic acid (ABA)] and three biotic stresses [powdery mildew (PM), downy mildew (DM), and root-knot nematode (RKN)] on OSCA family genes were also determined. The results indicated that cucumber OSCA genes play important roles in response to osmotic-related abiotic stresses and pathogen invasion. Overall, this study lays a foundation for research on the biological function and evolutionary process of OSCA family genes in cucumber.

Abstract

Calcium ions are important second messengers, playing an important role in the signal transduction pathways. Hyperosmolality gated calcium-permeable channels (OSCA) gene family members play critical modulating roles in response to osmotic-related abiotic stress as well as other abiotic and biotic stresses, which has been reported in many plant species such as Arabidopsis, rice, maize, and wheat. However, there has been no report about the identification and expression profile of the OSCA genes in cucumber. In this study, a total of nine OSCA genes were identified, which are unevenly distributed on the six chromosomes of cucumber. Phylogenetic analysis revealed that the OSCAs of cucumber, Arabidopsis, rice and maize were clustered into four clades. The motif arrangement of CsOSCAs was strongly conserved, and the CsOSCA genes in each group had similar genetic structure. A total of 11 and 10 types of cis-elements related to hormone and stress, respectively, were identified in the promoter regions of CsOSCA genes. Gene expression analysis indicated that the CsOSCA genes have different expression patterns in various tissues, and some of them were regulated by three osmotic-related abiotic stresses (salt, drought and ABA) and three biotic stresses (powdery mildew, downy mildew, and root-knot nematode). As the first genome-wide identification and characterization of the OSCA gene family in cucumber, this study lays a foundation for research on the biological function and evolutionary process of this gene family, which is of great significance for exploiting stress resistant cucumber varieties.

Genome-Wide Identification and Characterization of the CC-NBS-LRR Gene Family in Cucumber (Cucumis sativus L.)

The NBS-LRR (NLR) gene family plays a pivotal role in regulating disease defense response in plants. Cucumber is one of the most important vegetable crops in the world, and various plant diseases, including powdery mildew (PM), cause severe losses in both cucumber productivity and quality annually. To characterize and understand the role of the CC-NBS-LRR(CNL) family of genes in disease defense response in cucumber plants, we performed bioinformatical analysis to characterize these genes systematically. We identified 33 members of the CNL gene family in cucumber plants, and they are distributed on each chromosome with chromosome 4 harboring the largest cluster of five different genes. The corresponding CNL family member varies in the number of amino acids and exons, molecular weight, theoretical isoelectric point (pI) and subcellular localization. Cis-acting element analysis of the CNL genes reveals the presence of multiple phytohormone, abiotic and biotic responsive elements in their promoters, suggesting that these genes might be responsive to plant hormones and stress. Phylogenetic and synteny analysis indicated that the CNL proteins are conserved evolutionarily in different plant species, and they can be divided into four subfamilies based on their conserved domains. MEME analysis and multiple sequence alignment showed that conserved motifs exist in the sequence of CNLs. Further DNA sequence analysis suggests that CsCNL genes might be subject to the regulation of different miRNAs upon PM infection. By mining available RNA-seq data followed by real-time quantitative PCR (qRT-PCR) analysis, we characterized expression patterns of the CNL genes, and found that those genes exhibit a temporospatial expression pattern, and their expression is also responsive to PM infection, ethylene, salicylic acid, and methyl jasmonate treatment in cucumber plants. Finally, the CNL genes targeted by miRNAs were predicted in cucumber plants. Our results in this study provided some basic information for further study of the functions of the CNL gene family in cucumber plants.

Better Agreement of Human Transcriptomic and Proteomic Cancer Expression Data at the Molecular Pathway Activation Level

Previously, we have shown that the aggregation of RNA-level gene expression profiles into quantitative molecular pathway activation metrics results in lesser batch effects and better agreement between different experimental platforms. Here, we investigate whether pathway level of data analysis provides any advantage when comparing transcriptomic and proteomic data. We compare the paired proteomic and transcriptomic gene expression and pathway activation profiles obtained for the same human cancer biosamples in The Cancer Genome Atlas (TCGA) and the NCI Clinical Proteomic Tumor Analysis Consortium (CPTAC) projects, for a total of 755 samples of glioblastoma, breast, liver, lung, ovarian, pancreatic, and uterine cancers. In a CPTAC assay, expression levels of 15,112 protein-coding genes were profiled using the Thermo QE series of mass spectrometers. In TCGA, RNA expression levels of the same genes were obtained using the Illumina HiSeq 4000 engine for the same biosamples. At the gene level, absolute gene expression values are compared, whereas pathway-grade comparisons are made between the pathway activation levels (PALs) calculated using average sample-normalized transcriptomic and proteomic profiles. We observed remarkably different average correlations between the primary RNA- and protein expression data for different cancer types: Spearman Rho between 0.017 (p = 1.7 × 10−13) and 0.27 (p < 2.2 × 10−16). However, at the pathway level we detected overall statistically significantly higher correlations: averaged Rho between 0.022 (p < 2.2 × 10−16) and 0.56 (p < 2.2 × 10−16). Thus, we conclude that data analysis at the PAL-level yields results of a greater similarity when comparing high-throughput RNA and protein expression profiles.

Comparative genomics analysis of bHLH genes in cucurbits identifies a novel gene regulating cucurbitacin biosynthesis

The basic helix-loop-helix (bHLH) family of transcription factors (TFs) participate in a variety of biological regulatory processes in plants, and have undergone significant expansion during land plant evolution by gene duplications. In cucurbit crops, several bHLH genes have been found to be responsible for the agronomic traits such as bitterness. However, the characterization of bHLH genes across the genomes of cucurbit species has not been reported, and how they have evolved and diverged remains largely unanswered. Here we identified 1160 bHLH genes in seven cucurbit crops and performed a comprehensive comparative genomics analysis. We determined orthologous and paralogous bHLH genes across cucurbit crops by syntenic analysis between or within species. Orthology and phylogenetic analysis of the tandem-duplicated bHLH genes in the Bt cluster which regulate the biosynthesis of cucurbitacins suggest that this cluster is derived from three ancestral genes after the cucurbit-common tetraploidization event. Interestingly, we identified a new conserved cluster paralogous to the Bt cluster that includes two tandem bHLH genes, and the evolutionary history and expression profiles of these two genes in the new cluster suggest the involvement of one gene (Brp) in the regulation of cucurbitacin biosynthesis in roots. Further biochemical and transgenic assays in melon hairy roots support the function of Brp. This study provides useful information for further investigating the functions of bHLH TFs and novel insights into the regulation of cucurbitacin biosynthesis in cucurbit crops and other plants.

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.

Genome-wide identification and expression profile of YABBY genes in Averrhoa carambola

Background

Members of the plant-specific YABBY gene family are thought to play an important role in the development of leaf, flower, and fruit. The YABBY genes have been characterized and regarded as vital contributors to fruit development in Arabidopsis thaliana and tomato, in contrast to that in the important tropical economic fruit star fruit (Averrhoa carambola), even though its genome is available.

Methods

In the present study, a total of eight YABBY family genes (named from AcYABBY1 to AcYABBY8) were identified from the genome of star fruit, and their phylogenetic relationships, functional domains and motif compositions, physicochemical properties, chromosome locations, gene structures, protomer elements, collinear analysis, selective pressure, and expression profiles were further analyzed.

Results

Eight AcYABBY genes (AcYABBYs) were clustered into five clades and were distributed on five chromosomes, and all of them had undergone negative selection. Tandem and fragment duplications rather than WGD contributed to YABBY gene number in the star fruit. Expression profiles of AcYABBYs from different organs and developmental stages of fleshy fruit indicated that AcYABBY4 may play a specific role in regulating fruit size. These results emphasize the need for further studies on the functions of AcYABBYs in fruit development.

Genome-wide identification and expression analysis of the cucumber PYL gene family

Abscisic acid (ABA) is a very important hormone in plants. It regulates growth and development of plants and plays an important role in biotic and abiotic stresses. The Pyrabactin resistance 1-like (PYR/PYL) proteins play a central role in ABA signal transduction pathways. The working system of PYL genes in cucumber, an important economical vegetable (Cucumis sativus L.), has not been fully studied yet. Through bioinformatics, a total of 14 individual PYL genes were identified in Chinese long '9930' cucumber. Fourteen PYL genes were distributed on six chromosomes of cucumber, and their encoded proteins predicted to be distributed in cytoplasm and nucleus. Based on the phylogenetic analysis, the PYL genes of cucumber, Arabidopsis, rice, apple, Brachypodium distachyon and soybeancould be classified into three groups. Genetic structures and conserved domains analysis revealed that CsPYL genes in the same group have similar exons and conserved domains. By predicting cis-elements in the promoters, we found that all CsPYL members contained hormone and stress-related elements. Additionally, the expression patterns of CsPYL genes were specific in tissues. Finally, we further examined the expression of 14 CsPYL genes under ABA, PEG, salt stress. The qRT-PCR results showed that most PYL gene expression levels were up-regulated. Furthermore, with different treatments about 3h, the relative expression of PYL8 was up-regulated and more than 20 times higher than 0h. It indicated that this gene may play an important role in abiotic stress.

Terpene synthases in cucumber (Cucumis sativus) and their contribution to herbivore-induced volatile terpenoid emission

Terpenoids play important roles in flavour, pollinator attraction and defence of plants. In cucumber (Cucumis sativus) they are important components of the herbivore-induced plant volatile blend that attracts natural enemies of herbivores. We annotated the cucumber TERPENE SYNTHASE gene (CsTPS) family and characterized their involvement in the response towards herbivores with different feeding guilds using a combined molecular and biochemical approach. Transcripts of multiple CsTPS genes were upregulated in leaves upon herbivory and the products generated by the expressed proteins match the terpenoids recorded in the volatile blend released by herbivore-damaged leaves. Spatial and temporal analysis of the promoter activity of CsTPS genes showed that cell content-feeding spider mites (Tetranychus urticae) and thrips (Frankliniella occidentalis) induced promoter activity of CsTPS9 and CsTPS19 within hours after initiation of infestation, while phloem-feeding aphids (Myzus persicae) induced CsTPS2 promoter activity. Our findings offer detailed insights into the involvement of the TPS gene family in the dynamics and fine-tuning of the emission of herbivore-induced plant volatiles in cucumber, and open a new avenue to understand molecular mechanisms that affect plant-herbivore interactions.

UNAGI: Yeast Transcriptome Reconstruction and Gene Discovery Using Nanopore Sequencing

Computational approaches are the main approaches used in genome annotation. However, accuracy is low. Untranslated regions are not identified, complex isoforms are not predicted correctly and discovery rate of noncoding RNA is low. RNA-seq has revolutionized transcriptome reconstruction over the last decade. However, fragmentation included in cDNA sequencing leads to information loss, requiring transcripts to be assembled and reconstructed, thus affecting the accuracy of reconstructed transcriptome. Recently, long-read sequencing has been introduced with technologies such as Oxford Nanopore sequencing. cDNA is sequenced directly without fragmentation producing long reads that don't need to be assembled keeping the transcript structure intact and increasing the accuracy of transcriptome reconstruction.Here we present a protocol and a pipeline to reconstruct the transcriptome of compact genomes including yeasts. It involves generating full-length cDNA and using Oxford Nanopore ligation-based sequencing kit to sequence multiple samples in the same run. The pipeline (1) strands the generated long reads, (2) corrects the reads by mapping them to the reference genome, (3) identifies transcripts including 5'UTR and 3'UTR, (4) profiles the isoforms, filtering out artifacts resulting from low accuracy in sequencing, and (5) improves accuracy of provided annotations. Using long reads improves the accuracy of transcriptome reconstruction and helps in discovering a significant number of novel RNAs.

Transcriptome analysis of ovary culture-induced embryogenesis in cucumber (Cucumis sativus L.)

Background.

Ovary culture is a useful technique used to generate double haploid (DH) cucumber (Cucumis sativus L.) plants. However, cucumber ovary culture have a low rate of embryo induction and plant regeneration. Moreover, the cucumber embryogenesis mechanism remains unclear. In this study, we explored the molecular basis of cucumber embryogenesis in order to establish a foundation for a more efficient ovary culture method. Using transcriptome sequencing, we also investigated the differential expression of genes during the embryogenesis process.

Methods.

Cytological and morphological observations have divided cucumber ovary culture into three stages: early embryo development (T0), embryo morphogenesis (T1, T2, T3 and T4), and shoot formation (T5). We selected six key time points for transcriptome sequencing and analysis: T0 (the ovules were cultured for 0 d), T1 (the ovules were cultured for 2 d), T2 (the embryos were cultured for 10 d), T3 (the embryos were cultured for 20 d), T4 (the embryos were cultured for 30 d), and T5 (the shoots after 60 d culture).

Results.

We used cytology and morphology to observe the characteristics of the cucumber’s developmental transformation during embryogenesis and plant regeneration. The differentially expressed genes(DEGs) at developmental transition points were analyzed using transcriptome sequencing. In the early embryogenesis stage, the cells expanded, which was the signal for gametophytes to switch to the sporophyte development pathway. RNA-seq revealed that when compared to the fresh unpollinated ovaries, there were 3,468 up-regulated genes in the embryos, including hormone signal transduction genes, hormone response genes, and stress-induced genes. The reported embryogenesis-related genes BBM, HSP90 and AGL were also actively expressed during this stage. In the embryo morphogenesis stage (from cell division to cotyledon-embryo formation), 480 genes that functioned in protein complex binding, microtubule binding, tetrapyrrole binding, tubulin binding and other microtubule activities were continuously up-regulated during the T1, T2, T3 and T4 time points. This indicated that the cytoskeleton structure was continuously being built and maintained by the action of microtubule-binding proteins and enzyme modification. In the shoot formation stage, 1,383 genes were up-regulated that were mainly enriched in phenylpropanoid biosynthesis, plant hormone signal transduction, phenylalanine metabolism, and starch and sucrose metabolism. These up-regualted genes included six transcription factors that contained a B3 domain, nine genes in the AP2/ERF family, and two genes encoding WUS homologous domain proteins.

Conclusions.

Evaluation of molecular gynogenesis events may contribute to a better understanding of the molecular mechanism of cucumber ovarian culture.

HECATE2 acts with GLABROUS3 and Tu to boost cytokinin biosynthesis and regulate cucumber fruit wart formation

Warty fruit in cucumber (Cucumis sativus L.) is an important quality trait that greatly affects fruit appearance and market value. The cucumber wart consists of fruit trichomes (spines) and underlying tubercules, in which the existence of spines is prerequisite for tubercule formation. Although several regulators have been reported to mediate spine or tubercule formation, the direct link between spine and tubercule development remains unknown. Here, we found that the basic Helix-Loop-Helix (bHLH) gene HECATE2 (CsHEC2) was highly expressed in cucumber fruit peels including spines and tubercules. Knockout of CsHEC2 by the CRISPR/Cas9 system resulted in reduced wart density and decreased cytokinin (CTK) accumulation in the fruit peel, whereas overexpression of CsHEC2 led to elevated wart density and CTK level. CsHEC2 is directly bound to the promoter of the CTK hydroxylase-like1 gene (CsCHL1) that catalyzes CTK biosynthesis, and activated CsCHL1 expression. Moreover, CsHEC2 physically interacted with GLABROUS3 (CsGL3, a key spine regulator) and Tuberculate fruit (CsTu, a core tubercule formation factor), and such interactions further enhanced CsHEC2-mediated CsCHL1 expression. These data suggested that CsHEC2 promotes wart formation by acting as an important cofactor for CsGL3 and CsTu to directly stimulate CTK biosynthesis in cucumber. Thus, CsHEC2 can serve as a valuable target for molecular breeding of cucumber varieties with different wart density requirements.

Comprehensive analysis of 14-3-3 family genes and their responses to cold and drought stress in cucumber

The 14-3-3 proteins play essential roles in regulating various biological processes and abiotic stress responses in plants. However, there have been few studies of 14-3-3 family members in cucumber. Here, we identified a total of ten 14-3-3 genes (named as CsGF14a-j) in the cucumber genome. These genes are unevenly distributed across six cucumber chromosomes, and six of them were found to be segmentally duplicated. A phylogenetic analysis of 14-3-3 proteins in cucumber and other plant species showed that they could be divided into two distinct groups (ε and non-ε). Members in the same group tend to have similar exon-intron structure and conserved motif patterns. Several hormone-, stress- and development-related cis-elements associated with transcriptional regulation were found in the promoters of CsGF14 genes. RNA-seq data showed that most CsGF14 genes have broad expression in different tissues, and some had preferential expression in specific tissues and variable expression at certain developmental stages during fruit development. Quantitative real-time PCR (qRT-PCR) results revealed that nearly all tested CsGF14 genes were significantly up-regulated under cold and drought stress at certain time points. These results provide important information about the functions of CsGF14 genes in cucumber.

Genome-wide analysis of valine-glutamine motif-containing proteins related to abiotic stress response in cucumber (Cucumis sativus L.)

BACKGROUND

Cucumber (Cucumis sativus L.) is one of the most important economic crops and is susceptible to various abiotic stresses. The valine-glutamine (VQ) motif-containing proteins are plant-specific proteins with a conserved "FxxhVQxhTG" amino acid sequence that regulates plant growth and development. However, little is known about the function of VQ proteins in cucumber.

RESULTS

In this study, a total of 32 CsVQ proteins from cucumber were confirmed and characterized using comprehensive genome-wide analysis, and they all contain a conserved motif with 10 variations. Phylogenetic tree analysis revealed that these CsVQ proteins were classified into nine groups by comparing the CsVQ proteins with those of Arabidopsis thaliana, melon and rice. CsVQ genes were distributed on seven chromosomes. Most of these genes were predicted to be localized in the nucleus. In addition, cis-elements in response to different stresses and hormones were observed in the promoters of the CsVQ genes. A network of CsVQ proteins interacting with WRKY transcription factors (CsWRKYs) was proposed. Moreover, the transcripts of CsVQ gene were spatio-temporal specific and were induced by abiotic adversities. CsVQ4, CsVQ6, CsVQ16-2, CsVQ19, CsVQ24, CsVQ30, CsVQ32, CsVQ33, and CsVQ34 were expressed in the range of organs and tissues at higher levels and could respond to multiple hormones and different stresses, indicating that these genes were involved in the response to stimuli.

CONCLUSIONS

Together, our results reveal novel VQ resistance gene resources, and provide critical information on CsVQ genes and their encoded proteins, which supplies important genetic basis for VQ resistance breeding of cucumber plants.

Characterization of Germin-like Proteins (GLPs) and Their Expression in Response to Abiotic and Biotic Stresses in Cucumber

Germins and germin-like proteins (GLPs) are glycoproteins closely associated with plant development and stress response in the plant kingdom. Here, we carried out genome-wide identification and expression analysis of the GLP gene family in cucumber to study their possible functions. A total of 38 GLP genes were identified in cucumber, which could be mapped to six out of the seven cucumber chromosomes. A phylogenetic analysis of the GLP members from cucumber, Arabidopsis and rice showed that these GLPs could be divided into six groups, and cucumber GLPs in the same group had highly similar conserved motif distribution and gene structure. Gene duplication analysis revealed that six cucumber GLP genes were located in the segmental duplication regions of cucumber chromosomes, while 14 genes were associated with tandem duplications. Tissue expression profiles of cucumber GLP genes showed that many genes were preferentially expressed in specific tissues. In addition, some cucumber GLP genes were differentially expressed under salt, drought and ABA treatments, as well as under DM inoculation. Our results provide important information for the functional identification of GLP genes in the growth, development and stress response of cucumber.

Photosynthetic contribution and characteristics of cucumber stems and petioles

BACKGROUND

Photosynthesis in the green leafless blade tissues or organs of plants has been studied in some plants, but the photosynthetic characteristics of stems and petioles are poorly understood. Cucurbitaceous plants are climbing plants that have substantial stem and petiole biomass. Understanding the photosynthetic contribution of cucumber stems and petioles to their growth and the underlying molecular mechanisms are important for the regulating of growth in cucumber production.

RESULTS

In this study, the photosynthetic capacity of cucumber stems and petioles were determined by ¹⁴CO₂ uptake. The total carbon fixed by the stems and petioles was approximately 4% of that fixed by one leaf blade in the cucumber seedling stage, while the proportion of the carbon accumulated in the stems and petioles that redistributed to sink organs (roots and shoot apexes) obviously increased under leafless conditions. The photosynthetic properties of cucumber stems and petioles were studied using a combination of electron microscopy and isotope tracers to compare these properties of stems and petioles with those of leaf blade using two genotypes of cucumber (dark green and light green). Compared with those of the leaf blades, the chlorophyll contents of the cucumber stems and petioles were lower, and the stems and petioles had lower chloroplast numbers and lower stoma numbers but higher thylakoid grana lamella numbers and larger stoma sizes. The Chl a/b ratios were also decreased in the petioles and stems compared with those in the leaf blades. The total photosynthetic rates of the stems and petioles were equivalent to 6 ~ 8% of that of one leaf blade, but the respiration rates were similar in all the three organs, with an almost net 0 photosynthetic rate in the stems and petioles. Transcriptome analysis showed that compared with the leaf blades, the stems and petioles has significantly different gene expression levels in photosynthesis, porphyrin and chlorophyll metabolism; photosynthetic antenna proteins; and carbon fixation. PEPC enzyme activities were higher in the stems and petioles than in the leaf blades, suggesting that the photosynthetic and respiratory mechanisms in stems and petioles are different from those in leaf blade, and these results are consistent with the gene expression data.

CONCLUSIONS

In this study, we confirmed the photosynthetic contribution to the growth of cucumber stems and petioles, and showed their similar photosynthetic patterns in the terms of anatomy, molecular biology and physiology, which were different from those of cucumber leaf blades.

Genome-Wide Identification, Phylogenetic and Expression Pattern Analysis of GATA Family Genes in Cucumber (Cucumis sativus L.)

GATA transcription factors are a class of transcriptional regulatory proteins that contain a characteristic type-IV zinc finger DNA-binding domain, which play important roles in plant growth and development. The GATA gene family has been characterized in various plant species. However, GATA family genes have not been identified in cucumber. In this study, 26 GATA family genes were identified in cucumber genome, whose physicochemical characteristics, chromosomal distributions, phylogenetic tree, gene structures conserved motifs, cis-regulatory elements in promoters, homologous gene pairs, downstream target genes were analyzed. Tissue expression profiles of cucumber GATA family genes exhibited that 17 GATA genes showed constitutive expression, and some GATA genes showed tissue-specific expression patterns. RNA-seq analysis of green and virescent leaves revealed that seven GATA genes might be involved in the chloroplast development and chlorophyll biosynthesis. Importantly, expression patterns analysis of GATA genes in response to abiotic and biotic stresses indicated that some GATA genes respond to either abiotic stress or biotic stress, some GATA genes such as Csa2G162660, Csa3G017200, Csa3G165640, Csa4G646060, Csa5G622830 and Csa6G312540 were simultaneously functional in resistance to abiotic and biotic stresses. Overall, this study will provide useful information for further analysis of the biological functions of GATA factors in cucumber.

The amino acid permease (AAP) genes CsAAP2A and SlAAP5A/B are required for oomycete susceptibility in cucumber and tomato

Cucurbit downy mildew (DM), caused by the obligate biotroph Pseudoperonospora cubensis, is a destructive disease in cucumber. A valuable source of DM resistance is the Indian cucumber accession PI 197088, which harbours several quantitative trait loci (QTLs) contributing to quantitatively inherited DM resistance. With a combination of fine-mapping and transcriptomics, we identified Amino Acid Permease 2A (CsAAP2A) as a candidate gene for QTL DM4.1.3. Whole-genome and Sanger sequencing revealed the insertion of a Cucumis Mu-like element (CUMULE) transposon in the allele of the resistant near-isogenic line DM4.1.3. To confirm whether loss of CsAAP2A contributes to partial DM resistance, we performed targeting induced local lesions in genomes on a DM-susceptible cucumber genotype to identify an additional csaap2a mutant, which indeed was partially DM resistant. In view of the loss of the putative function as amino acid transporter, we measured amino acids in leaves. We found that DM-inoculated leaves of line DM4.1.3 (with the csaap2a mutation) contained significantly fewer amino acids than wild-type cucumber. The decreased flow of amino acids towards infected leaves in csaap2a plants compared to the wild type might explain the resistant phenotype of the mutant, as this would limit the available nutrients for the pathogen and thereby its fitness. To examine whether AAP genes play a conserved role as susceptibility factors in plant-oomycete interactions, we made targeted mutations in two AAP genes from tomato and studied the effect on susceptibility to Phytophthora infestans. We conclude that not only CsAAP2A but also SlAAP5A/SlAAP5B are susceptibility genes for oomycete pathogens.

Complexity untwined: The structure and function of cucumber (Cucumis sativus L.) shoot phloem

Cucurbit phloem is complex, with large sieve tubes on both sides of the xylem (bicollateral phloem), and extrafascicular elements that form an intricate web linking the rest of the vasculature. Little is known of the physical interconnections between these networks or their functional specialization, largely because the extrafascicular phloem strands branch and turn at irregular angles. Here, export in the phloem from specific regions of the lamina of cucumber (Cucumis sativus L.) was mapped using carboxyfluorescein and ¹⁴ C as mobile tracers. We also mapped vascular architecture by conventional microscopy and X-ray computed tomography using optimized whole-tissue staining procedures. Differential gene expression in the internal (IP) and external phloem (EP) was analyzed by laser-capture microdissection followed by RNA-sequencing. The vascular bundles of the lamina form a nexus at the petiole junction, emerging in a predictable pattern, each bundle conducting photoassimilate from a specific region of the blade. The vascular bundles of the stem interconnect at the node, facilitating lateral transport around the stem. Elements of the extrafascicular phloem traverse the stem and petiole obliquely, joining the IP and EP of adjacent bundles. Using pairwise comparisons and weighted gene coexpression network analysis, we found differences in gene expression patterns between the petiole and stem and between IP and EP, and we identified hub genes of tissue-specific modules. Genes related to transport were expressed primarily in the EP while those involved in cell differentiation and development as well as amino acid transport and metabolism were expressed mainly in the IP.

Identification and Transcriptional Analysis of Zinc Finger-Homeodomain (ZF-HD) Family Genes in Cucumber

Zinc finger-homeodomain (ZF-HD) proteins encode a family of plant-specific transcription factors that play essential roles in regulating plant growth and development as well as responses to abiotic/biotic stresses by activating or repressing the target genes. In this study, genome-wide characterization and expression profiling of the ZF-HD gene family in cucumber (Cucumis sativus) were performed for the first time. By using bioinformatics approaches, a total of 13 ZF-HD genes (designated as CsMIF1-CsMIF3 and CsZHD1-CsZHD10) were identified in the cucumber genome, which were unevenly distributed on six chromosomes. According to the phylogenetic analysis of cucumber and other species, they were divided into two distinct families, MINI ZINC FINGER (MIF) and zinc finger-homeodomain (ZHD), and the ZHD family was further divided into six subfamilies (ZHDI-ZHDVI). CsZF-HD members were mostly conserved in each subfamily with minor variations in motif distribution, and gene structure analysis showed that the CsZF-HD genes had only one intron or no intron at all. Expression analysis showed that most CsZF-HD genes had tissue-specific expression patterns, and some of them exhibited highly variable expression during fruit development. qRT-PCR results indicated that the selected CsZF-HD genes were responsive to drought stress, and some of them were differentially expressed in response to the inoculation of powdery mildew (PM) and downy mildew (DM) based on publicly available RNA-seq data. The results lay the foundation for further functional analysis of the ZF-HD genes and explore their potential application to the improvement of stress tolerance in cucumber.

Genome-Wide Characterization of HSP90 Gene Family in Cucumber and Their Potential Roles in Response to Abiotic and Biotic Stresses

Heat shock protein 90 (HSP90) possesses critical functions in plant developmental control and defense reactions. The HSP90 gene family has been studied in various plant species. However, the HSP90 gene family in cucumber has not been characterized in detail. In this study, a total of six HSP90 genes were identified from the cucumber genome, which were distributed to five chromosomes. Phylogenetic analysis divided the cucumber HSP90 genes into two groups. The structural characteristics of cucumber HSP90 members in the same group were similar but varied among different groups. Synteny analysis showed that only one cucumber HSP90 gene, Csa1G569290, was conservative, which was not collinear with any HSP90 gene in Arabidopsis and rice. The other five cucumber HSP90 genes were collinear with five Arabidopsis HSP90 genes and six rice HSP90 genes. Only one pair of paralogous genes in the cucumber HSP90 gene family, namely one pair of tandem duplication genes (Csa1G569270/Csa1G569290), was detected. The promoter analysis showed that the promoters of cucumber HSP90 genes contained hormone, stress, and development-related cis-elements. Tissue-specific expression analysis revealed that only one cucumber HSP90 gene Csa3G183950 was highly expressed in tendril but low or not expressed in other tissues, while the other five HSP90 genes were expressed in all tissues. Furthermore, the expression levels of cucumber HSP90 genes were differentially induced by temperature and photoperiod, gibberellin (GA), downy mildew, and powdery mildew stimuli. Two cucumber HSP90 genes, Csa1G569270 and Csa1G569290, were both differentially expressed in response to abiotic and biotic stresses, which means that these two HSP90 genes play important roles in the process of cucumber growth and development. These findings improve our understanding of cucumber HSP90 family genes and provide preliminary information for further studies of cucumber HSP90 gene functions in plant growth and development.

Identification and Expression Analysis of the CsMYB Gene Family in Root Knot Nematode-Resistant and Susceptible Cucumbers

MYB (myeloblastosis) transcription factors (TFs) play important roles in controlling various physiological processes in plants, such as responses to biotic and abiotic stress, metabolism, and defense. A previous study identified a gene, Csa6G410090, encoding a plant lipid transfer protein (LTP), as a possible regulator in cucumber (Cucumis sativus L.) of the resistance response to root-knot nematode (RKN) [Meloidogyne incognita Kofoid and White (Chitwood)]. Myb-type DNA-binding TFs were presumed to regulate downstream genes expression, including LTPs, however, the regulation mechanism remained unclear. To elucidate whether and which MYB TFs may be involved in regulation of the resistance response, this study identified 112 genes as candidate members of the CsMYB gene family by combining CDD and SMART databases, using the Hidden Markov Model (HMM) and manual calibration. Within this group, ten phylogenetic subgroups were resolved according to sequence-based classification, consistent with results from comprehensive investigation of gene structure, conserved motifs, chromosome locations, and cis-element analysis. Distribution and collinearity analysis indicated that amplification of the CsMYB gene family in cucumber has occurred mainly through tandem repeat events. Spatial gene expression analysis showed that 8 CsMYB genes were highly expressed at differing levels in ten different tissues or organs. The roots of RKN-resistant and susceptible cucumbers were inoculated with M. incognita, finding that CsMYB (Csa6G538700, Csa1G021940, and Csa5G641610) genes showed up-regulation coincident with upregulation of the "hub" gene LTP (Csa6G410090) previously implicated as a major gene in the resistance response to RKN in cucumber. Results of this study suggest hypotheses regarding the elements and regulation of the resistant response as well as possible RKN resistance-enhancing strategies in cucumber and perhaps more broadly in plants.

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

Snake gourd (Trichosanthes anguina L.), which belongs to the Cucurbitaceae family, is a popular ornamental and food crop species with medicinal value and is grown in many parts of the world. Although progress has been made in its genetic improvement, the organization, composition, and evolution of the snake gourd genome remain largely unknown. Here, we report a high-quality genome assembly for snake gourd, comprising 202 contigs, with a total size of 919.8 Mb and an N50 size of 20.1 Mb. These findings indicate that snake gourd has one of the largest genomes of Cucurbitaceae species sequenced to date. The snake gourd genome assembly harbors 22,874 protein-coding genes and 80.0% of the genome consists of repetitive sequences. Phylogenetic analysis reveals that snake gourd is closely related to sponge gourd but diverged from their common ancestor ~33-47 million years ago. The genome sequence reported here serves as a valuable resource for snake gourd genetic research and comparative genomic studies in Cucurbitaceae and other plant species. In addition, fruit transcriptome analysis reveals the candidate genes related to quality traits during snake gourd fruit development and provides a basis for future research on snake gourd fruit development and ripening at the transcript level.

Analysis of QTL DM4.1 for Downy Mildew Resistance in Cucumber Reveals Multiple subQTL: A Novel RLK as Candidate Gene for the Most Important subQTL

One of the biggest problems in cucumber cultivation is cucurbit downy mildew (DM), caused by the obligate biotroph Pseudoperonospora cubensis. Whereas DM in cucumber was previously efficiently controlled by the dm-1 gene from Indian cucumber accession PI 197087, this resistance was broken by new DM strains, prompting the search for novel sources of resistance. A promising source of resistance is the wild cucumber accession PI 197088. It was previously shown that DM resistance in this genotype inherits polygenically. In this paper, we put the focus on one of the QTL, DM4.1 that is located on chromosome 4. QTL DM4.1 was shown to consist of three subQTL: DM4.1.1 affected pathogen-induced necrosis, DM4.1.2 was shown to have an additive effect on sporulation, and DM4.1.3 had a recessive effect on chlorosis as well as an effect on sporulation. Near-isogenic lines (NILs) were produced by introgressing the subQTLs into a susceptible cucumber line (HS279) with good horticultural traits. Transcriptomic analysis revealed that many genes in general, and defense pathway genes in particular, were differentially expressed in NIL DM4.1.1/.2 compared to NIL DM4.1.3 and the susceptible parent HS279. This indicates that the resistance from subQTL DM4.1.1 and/or subQTL DM4.1.2 likely involves defense signaling pathways, whereas resistance due to subQTL DM4.1.3 is more likely to be independent of known defense pathways. Based on fine-mapping data, we identified the RLK gene CsLRK10L2 as a likely candidate for subQTL DM4.1.2, as this gene was found to have a loss-of-function mutation in the susceptible parent HS279, and was strongly upregulated by P. cubensis inoculation in NIL DM4.1.1/.2. Heterologous expression of this gene triggered necrosis, providing further evidence that this gene is indeed causal for subQTL DM4.1.2.

Genome-wide analysis of the WRKY gene family in the cucumber genome and transcriptome-wide identification of WRKY transcription factors that respond to biotic and abiotic stresses

BACKGROUND

Cucumber (Cucumis sativus L.) is an economically important vegetable crop species. However, it is susceptible to various abiotic and biotic stresses. WRKY transcription factors play important roles in plant growth and development, particularly in the plant response to biotic and abiotic stresses. However, little is known about the expression pattern of WRKY genes under different stresses in cucumber.

RESULTS

In the present study, an analysis of the new assembly of the cucumber genome (v3.0) allowed the identification of 61 cucumber WRKY genes. Phylogenetic and synteny analyses were performed using related species to investigate the evolution of the cucumber WRKY genes. The 61 CsWRKYs were classified into three main groups, within which the gene structure and motif compositions were conserved. Tissue expression profiles of the WRKY genes demonstrated that 24 CsWRKY genes showed constitutive expression (FPKM > 1 in all samples), and some WRKY genes showed organ-specific expression, suggesting that these WRKYs might be important for plant growth and organ development in cucumber. Importantly, analysis of the CsWRKY gene expression patterns revealed that five CsWRKY genes strongly responded to both salt and heat stresses, 12 genes were observed to be expressed in response to infection from downy mildew and powdery mildew, and three CsWRKY genes simultaneously responded to all treatments analysed. Some CsWRKY genes were observed to be induced/repressed at different times after abiotic or biotic stress treatment, demonstrating that cucumber WRKY genes might play different roles during different stress responses and that their expression patterns vary in response to stresses.

CONCLUSIONS

Sixty-one WRKY genes were identified in cucumber, and insight into their classification, evolution, and expression patterns was gained in this study. Responses to different abiotic and biotic stresses in cucumber were also investigated. Our results provide a better understanding of the function of CsWRKY genes in improving abiotic and biotic stress resistance in cucumber.

Developmentally regulated activation of defense allows for rapid inhibition of infection in age-related resistance to Phytophthora capsici in cucumber fruit

Background

Age-related resistance (ARR) is a developmentally regulated phenomenon conferring resistance to pathogens or pests. Although ARR has been observed in several host-pathogen systems, the underlying mechanisms are largely uncharacterized. In cucumber, rapidly growing fruit are highly susceptible to Phytophthora capsici but become resistant as they complete exponential growth. We previously demonstrated that ARR is associated with the fruit peel and identified gene expression and metabolomic changes potentially functioning as preformed defenses.

Results

Here, we compare the response to infection in fruit at resistant and susceptible ages using microscopy, quantitative bioassays, and weighted gene co-expression analyses. We observed strong transcriptional changes unique to resistant aged fruit 2–4 h post inoculation (hpi). Microscopy and bioassays confirmed this early response, with evidence of pathogen death and infection failure as early as 4 hpi and cessation of pathogen growth by 8–10 hpi. Expression analyses identified candidate genes involved in conferring the rapid response including those encoding transcription factors, hormone signaling pathways, and defenses such as reactive oxygen species metabolism and phenylpropanoid biosynthesis.

Conclusion

The early pathogen death and rapid defense response in resistant-aged fruit provide insight into potential mechanisms for ARR, implicating both pre-formed biochemical defenses and developmentally regulated capacity for pathogen recognition as key factors shaping age-related resistance.

De novo assemblies of Luffa acutangula and Luffa cylindrica genomes reveal an expansion associated with substantial accumulation of transposable elements

Luffa spp. (sponge gourd or ridge gourd) is an economically important vegetable crop widely cultivated in China, India and Southeast Asia. Here, we employed PacBio long-read single-molecule real-time (SMRT) sequencing to perform de novo genome assemblies of two commonly cultivated Luffa species, L. acutangula and L. cylindrica. We obtained preliminary draft genomes of 734.6 Mb and 689.8 Mb with scaffold N50 of 786,130 and 578,616 bases for L. acutangula and L. cylindrica, respectively. We also applied long-range Chicago and HiC techniques to obtain the first chromosome-scale whole-genome assembly of L. acutangula. The final assembly contained 13 pseudomolecules, corresponding to the haploid chromosome number in Luffa spp. (1n = 13, 2n = 26). The sizes of the assembled Luffa genomes are approximately twice as large as the genome assemblies of related Cucurbitaceae. A large proportion of L. acutangula (62.17%; 456.69 Mb) and L. cylindrica (56.78%; 391.65 Mb) genome assemblies contained repetitive elements. Phylogenetic analyses revealed that the substantial accumulation of transposable elements likely contributed to the expansion of the Luffa genomes. We also investigated alternative splicing events in Luffa using full-length transcript sequences obtained from PacBio Isoform Sequencing (Iso-seq). While the predominant form of alternative splicing in most plant species examined was intron retention, alternative 3' acceptor site selection appeared to be a major event observed in Luffa. High-quality genome assemblies for L. acutangula and L. cylindrica reported here provide valuable resources for Luffa breeding and future genetics and comparative genomics studies in Cucurbitaceae.

Genome-Wide Analysis Reveals Transcription Factors Regulated by Spider-Mite Feeding in Cucumber (Cucumis sativus)

To gain insight into the regulatory networks that underlie the induced defense in cucumber against spider mites, genes encoding transcription factors (TFs) were identified in the cucumber (Cucumissativus) genome and their regulation by two-spotted spider mite (Tetranychusurticae) herbivory was analyzed using RNA-seq. Of the total 1212 annotated TF genes in the cucumber genome, 119 were differentially regulated upon spider-mite herbivory during a period of 3 days. These TF genes belong to different categories but the MYB, bHLH, AP2/ERF and WRKY families had the highest relative numbers of differentially expressed genes. Correlation analysis of the expression of TF genes with defense-associated genes during herbivory and pathogen infestation, and in different organs resulted in the putative identification of regulators of herbivore-induced terpenoid and green-leaf-volatile biosynthesis. Analysis of the cis-acting regulatory elements (CAREs) present in the promoter regions of the genes responsive to spider-mite feeding revealed potential TF regulators. This study describes the TF genes in cucumber that are potentially involved in the regulation of induced defense against herbivory by spider mites.

IC4R-2.0: Rice Genome Reannotation Using Massive RNA-seq Data

Genome reannotation aims for complete and accurate characterization of gene models and thus is of critical significance for in-depth exploration of gene function. Although the availability of massive RNA-seq data provides great opportunities for gene model refinement, few efforts have been made to adopt these precious data in rice genome reannotation. Here we reannotate the rice (Oryza sativa L. ssp. japonica) genome based on integration of large-scale RNA-seq data and release a new annotation system IC4R-2.0. In general, IC4R-2.0 significantly improves the completeness of gene structure, identifies a number of novel genes, and integrates a variety of functional annotations. Furthermore, long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs) are systematically characterized in the rice genome. Performance evaluation shows that compared to previous annotation systems, IC4R-2.0 achieves higher integrity and quality, primarily attributable to massive RNA-seq data applied in genome annotation. Consequently, we incorporate the improved annotations into the Information Commons for Rice (IC4R), a database integrating multiple omics data of rice, and accordingly update IC4R by providing more user-friendly web interfaces and implementing a series of practical online tools. Together, the updated IC4R, which is equipped with the improved annotations, bears great promise for comparative and functional genomic studies in rice and other monocotyledonous species. The IC4R-2.0 annotation system and related resources are freely accessible at http://ic4r.org/.

Transcriptional and metabolite analysis reveal a shift in direct and indirect defences in response to spider-mite infestation in cucumber (Cucumis sativus)

KEY MESSAGE

Cucumber plants adapt their transcriptome and metabolome as result of spider mite infestation with opposite consequences for direct and indirect defences in two genotypes. Plants respond to arthropod attack with the rearrangement of their transcriptome which lead to subsequent phenotypic changes in the plants' metabolome. Here, we analysed transcriptomic and metabolite responses of two cucumber (Cucumis sativus) genotypes to chelicerate spider mites (Tetranychus urticae) during the first 3 days of infestation. Genes associated with the metabolism of jasmonates, phenylpropanoids, terpenoids and L-phenylalanine were most strongly upregulated. Also, genes involved in the biosynthesis of precursors for indirect defence-related terpenoids were upregulated while those involved in the biosynthesis of direct defence-related cucurbitacin C were downregulated. Consistent with the observed transcriptional changes, terpenoid emission increased and cucurbitacin C content decreased during early spider-mite herbivory. To further study the regulatory network that underlies induced defence to spider mites, differentially expressed genes that encode transcription factors (TFs) were analysed. Correlation analysis of the expression of TF genes with metabolism-associated genes resulted in putative identification of regulators of herbivore-induced terpenoid, green-leaf volatiles and cucurbitacin biosynthesis. Our data provide a global image of the transcriptional changes in cucumber leaves in response to spider-mite herbivory and that of metabolites that are potentially involved in the regulation of induced direct and indirect defences against spider-mite herbivory.

Genotyping by RAD Sequencing Analysis Assessed the Genetic Distinctiveness of Experimental Lines and Narrowed Down the Genomic Region Responsible for Leaf Shape in Endive (Cichorium endivia L.)

The characterization of genetic diversity in elite breeding stocks is crucial for the registration and protection of new varieties. Moreover, experimental population structure analysis and information about the genetic distinctiveness of commercial materials are essential for crop breeding programs. The purpose of our research was to assess the genetic relationships of 32 endive (Cichorium endivia L.) breeding lines, 18 from var. latifolium (escarole) and 14 from var. crispum (curly), using heterologous Cichorium intybus-derived simple sequence repeats (SSR) markers and single-nucleotide polymorphisms (SNP) markers. We found that 14 out of 29 SSR markers were successfully amplified, but only 8 of them were related to polymorphic loci. To overcome the limitation of the low number of informative SSR marker loci, an alternative SNP-based approach was employed. The 4621 SNPs produced by a restriction site-associated DNA marker sequencing approach were able to fully discriminate the 32 endive accessions; most importantly, as many as 50 marker loci were found to distinguish the curly group from the escarole group. Interestingly, 24 of the marker loci mapped within a peripheral segment of chromosome 8 of lettuce (Lactuca sativa L.), spanning a chromosomal region of 49.6 Mb. Following Sanger sequencing-based validation, three genes were determined to carry nonsynonymous SNPs, and one of them matched a putative ortholog of AtELP1, subunit 1 of the Elongator complex. Considering that several previously characterized Elongator complex subunit mutants exhibited elongated and/or curly leaf phenotypes, this gene should be taken into consideration for a better understanding of the underlying mechanism controlling leaf shape in endive.

Identification and Expression Analysis of Stress-Associated Proteins (SAPs) Containing A20/AN1 Zinc Finger in Cucumber

Stress-associated proteins (SAPs) are a class of zinc finger proteins that confer tolerance to a variety of abiotic and biotic stresses in diverse plant species. However, in cucumber (Cucumis sativus L.), very little is known about the roles of SAP gene family members in regulating plant growth, development, and stress responses. In this study, a total of 12 SAP genes (named as CsSAP1-CsSAP12) were identified in the cucumber genome, which were unevenly distributed on six chromosomes. Gene duplication analysis detected one tandem duplication and two segmental duplication events. Phylogenetic analysis of SAP proteins from cucumber and other plants suggested that they could be divided into seven groups (sub-families), and proteins in the same group generally had the same arrangement of AN1 (ZnF-AN1) and A20 (ZnF-A20) domains. Most of the CsSAP genes were intronless and harbored a number of stress- and hormone-responsive cis-elements in their promoter regions. Tissue expression analysis showed that the CsSAP genes had a broad spectrum of expression in different tissues, and some of them displayed remarkable alteration in expression during fruit development. RT-qPCR results indicated that all the selected CsSAP genes displayed transcriptional responses to cold, drought, and salt stresses. These results enable the first comprehensive description of the SAP gene family in cucumber and lay a solid foundation for future research on the biological functions of CsSAP genes.

Genome-wide analysis of a putative lipid transfer protein LTP_2 gene family reveals CsLTP_2 genes involved in response of cucumber against root-knot nematode (Meloidogyne incognita)

Plant lipid transfer proteins (LTPs) are small basic proteins that play important roles in the regulation of various plant biological processes as well as the response to biotic and abiotic stresses. However, knowledge is limited on how this family of proteins is regulated in response to nematode infection in cucumber. In the present study, a total of 39 CsLTP_2 genes were identified by querying databases for cucumber-specific LTP_2 using a Hidden Markov Model approach and manual curation. The family has a five-cysteine motif (5CM) with the basic form CC-Xn-CXC-Xn-C, which differentiates it from typical nsLTPs. The members of CsLTP_2 were grouped into six families according to their structure and their phylogenetic relationships. Expression data of CsLTP_2 genes in 10 cucumber tissues indicated that they were tissue-specific genes. Two genes showed significant expression change in roots of resistant and susceptible lines during nematode infection, indicating their involvement in response to Meloidogyne incognita. This systematic analysis provides a foundation of knowledge for future studies of the biological roles of CsLTP_2 genes in cucumber in response to nematode infection and may help in the efforts to improve M. incognita-resistance breeding in cucumber.

Transcriptome Analysis Reveals the Different Response to Toxic Stress in Rootstock Grafted and Non-Grafted Cucumber Seedlings

Autotoxicity of root exudates is one of the main reasons for consecutive monoculture problem (CMP) in cucumber under greenhouse cultivation. Rootstock grafting may improve the tolerance of cucumber plants to autotoxic stress. To verify the enhanced tolerance to autotoxic stress and illuminate relevant molecular mechanism, a transcriptomic comparative analysis was performed between rootstock grafted (RG) and non-grafted (NG) cucumber plants by a simulation of exogenous cinnamic acid (CA). The present study confirmed that relatively stable plant growth, biomass accumulation, chlorophyll content, and photosynthesis was observed in RG than NG under CA stress. We identified 3647 and 2691 differentially expressed genes (DEGs) in NG and RG cucumber plants when compared to respective control, and gene expression patterns of RNA-seq was confirmed by qRT-PCR. Functional annotations revealed that DEGs response to CA stress were enriched in pathways of plant hormone signal transduction, MAPK signaling pathway, phenylalanine metabolism, and plant-pathogen interaction. Interestingly, the significantly enriched pathway of photosynthesis-related, carbon and nitrogen metabolism only identified in NG, and most of DEGs were down-regulated. However, most of photosynthesis, Calvin cycle, glycolysis, TCA cycle, and nitrogen metabolism-related DEGs exhibited not or slightly down-regulated in RG. In addition, several stress-related transcription factor families of AP2/ERF, bHLH, bZIP, MYB. and NAC were uniquely triggered in the grafted cucumbers. Overall, the results of this study suggest that rootstock grafting improve the tolerance of cucumber plants to autotoxic stress by mediating down-regulation of photosynthesis, carbon, and nitrogen metabolism-related DEGs and activating the function of stress-related transcription factor. The transcriptome dataset provides an extensive sequence resource for further studies of autotoxic mechanism at molecular level.

Phylogenetic informativeness analyses to clarify past diversification processes in Cucurbitaceae

Phylogenomic studies have so far mostly relied on genome skimming or target sequence capture, which suffer from representation bias and can fail to resolve relationships even with hundreds of loci. Here, we explored the potential of phylogenetic informativeness and tree confidence analyses to interpret phylogenomic datasets. We studied Cucurbitaceae because their small genome size allows cost-efficient genome skimming, and many relationships in the family remain controversial, preventing inferences on the evolution of characters such as sexual system or floral morphology. Genome skimming and PCR allowed us to retrieve the plastome, 57 single copy nuclear genes, and the nuclear ribosomal ITS from 29 species representing all but one tribe of Cucurbitaceae. Node support analyses revealed few inter-locus conflicts but a pervasive lack of phylogenetic signal among plastid loci, suggesting a fast divergence of Cucurbitaceae tribes. Data filtering based on phylogenetic informativeness and risk of homoplasy clarified tribe-level relationships, which support two independent evolutions of fringed petals in the family. Our study illustrates how formal analysis of phylogenomic data can increase our understanding of past diversification processes. Our data and results will facilitate the design of well-sampled phylogenomic studies in Cucurbitaceae and related families.