Analysis of differentially expressed genes during embryogenesis in ovary culture of cucumber (Cucumis sativus L.)

Analysis of differentially expressed genes during embryogenesis in ovary culture of cucumber (Cucumis sativus L.) | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research article Analysis of differentially expressed genes during embryogenesis in ovary culture of cucumber (Cucumis sativus L.) Ying Deng, Wenyuan Fu, Bing Tang, Lian Tao, Lu Zhang, Xia Zhou, and 2 more This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.2.10287/v1 This work is licensed under a CC BY 4.0 License Status: Posted Version 1 posted You are reading this latest preprint version Abstract Background: Ovary culture has been a useful way to generate double haploid (DH) plant in cucumber (Cucumis sativus L.). However, the rate of embryo induction is low, and the probability for the induced embryo to grow into normal embryo is low as well. This is largely due to unknown of the mechanism of embryogenesis in cucumber. In this study, the differentially expressed genes during embryogenesis, including the early stages of embryo formation, embryo maturation and shoot formation, was investigated with transcriptomic technique to set up basis for a more efficient ovary culture technology Results: Cytological observations led to suggestions that cell enlargement is the symbol that gametophytes had switched to the sporophyte development pathway during the early embryogenesis stage. In this stage, RNA-seq revealed 3468 up-regulated genes, 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. The total of 480 genes that function in protein complex binding, microtubule binding, tetrapyrrole binding, tubulin binding and other microtubule activities were continuously up-regulated during the embryo maturation stage, indicating that the cytoskeleton structure was continuously being built and maintained by the action of microtubule-binding proteins and enzyme modification during embryo development. In the shoot formation stage, 1383 genes were up-regulated, which were mainly enriched in phenylpropanoid biosynthesis, plant hormone signal transduction, phenylalanine metabolism, and starch and sucrose metabolism. The shoot formation stage might be regulated by 6 transcription factors that contained a B3 domain, 9 genes in the AP2/ERF family and 2 genes encoded WUS homologous domain proteins. Conclusions: Findings from this study offer a valuable framework for explaining the transcriptional regulatory mechanism underlying embryogenesis in cucumber ovary culture. Plant Molecular Biology and Genetics Cucumber Ovary culture Embryogenesis Transcriptome Differentially expressed genes Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Figure 6 Figure 7 Figure 8 Figure 9 Background Gametophyte culture involved a stress-induced reprogramming of male or female gametophytes to develop into embryo-like structures, which can be directly regenerated into completely homozygous, doubled haploid (DH). The study of male gametophyte (microspore) embryogenesis began in the 1960s [1]. Microspore embryogenesis in vitro can be easily monitored and provides a convenient experimental platform for large-scale physiological and biochemical analyses [2]. However, the study of gynogenesis is not simple; the embryo sac is small and embedded in surrounding tissues, making early embryos difficult to observe and isolate. Despite these difficulties, to inducing gynogenesis in vitro, especially when the haploid induced by the male gametophyte is not successful or the induction rate is too low and when self-incompatibility, male-sterile and dioecious plants occurs, is an effective way to obtain haploid or double haploid plants [3]. In vitro gynogenesis is mainly used in the Cucurbitaceae , and previous studies have revealed that this technique is still imperfect and that the embryogenesis rate is low [4-11]. The study of the mechanism of embryogenesis has therefore been delayed, and the understanding of the mechanisms underlying gynogenesis induction is extremely limited. In general, embryogenesis is a stress-induced phenomenon. Micropores are cultured and develop into embryos in vitro, and stress treatments involving cold or heat shock and hormones are used as trigger factors to induce gametocytes to follow the sporophyte development pathway [12]. The genes associated with the reprogramming phase and the early stage of embryogenesis were slowly characterized, such as AGL15 -related proteins and the AGL15 gene, the BBM gene and HSP gene [13-15]. More recently, the use of functional genomics tools has made it possible to identify more genes associated with different stages of microspore embryogenesis. In Brassica napus , the expression of embryogenesis-related genes such as BBM , LEC 1 and LEC2 was detected after culture under heat stress at 32 o C for three days; after incubation at 24 o C, embryo-specific expressed genes, such as LEC1 , LEC2 were detected [2]. The expression of genes associated with metabolism, chromosome remodeling, transcription, and translation signaling was up-regulated during the stress treatment stage in tobacco [16]. Similarly, genes related to metabolism, cell walls and the cell membrane, cell tissue control, cell communication and signal transduction were detected in the early stage of rape embryogenesis [17]. High levels of BBM and LEC gene expression have been confirmed in the early embryonic development of sweet pepper anther cultures [18]. The transition from microspores to developing embryos is mainly manifested in the induction of transcription factor genes that play an important role in early embryogenesis. Many genes are involved in hormone biosynthesis and plant hormone signal transduction and are involved in secondary metabolism [19]. Studies have shown that differentially expressed genes can be detected by gene chips in in vitro gynogenesis of cucumber in the early stage and have suggested that phenylalanine metabolism and phenylalanine synthesis may play an important role in the early development of cucumber in vitro [20]. Cucumber ( Cucumis sativus L.) is the fourth most important vegetable worldwide [21]. Cucumber production in China ranks first worldwide. To date, the cultivated area of cucumber in China has reached 1.258 Mha [22]. Cucumber is thought to be one of the oldest vegetable crops and has been domesticated in China for probably 2000 years [23]. To meet the needs of production, breeders are constantly looking for valuable germplasm resources of valuable characteristics, especially resistance to disease and environmental stresses such as cold, drought, or salt stress [24]. To facilitate these breeding efforts, we can collect extensive germplasm resources for variety improvements. However, cucumber is a cross-pollinated plant with obvious heterosis. Therefore, almost all the cucumber varieties in production are hybrids; however, the purification process of the hybrids is slow and long. In the breeding of parents, it takes 6-8 years of artificial self-breeding to develop a stable inbred line. To accelerate the purification of cucumber parents and improve breeding efficiency, haploid gametophyte cultures can be used to induce embryoids, and homozygous double haploids can be obtained in 1-2 years. DH technology is a powerful tool to speed up plant breeding. However, additional research is needed to improve our understanding of the genes or the roles of genes involved in embryogenesis and the mechanism of haploid induction in embryo sacs [25]. This study was based on the highly efficient in vitro ovary culture technology system of cucumber (unpublished). We divided the process from acquiring embryogenic potential to plant regeneration into three stages: early embryo development, embryo maturation (from pre-embryos to cotyledon embryos) and the shoot formation stage. We evaluated the embryo morphology and transcriptomes at different developmental stages during embryogenesis in cucumbers via ovary culture. The metabolism and biological process of embryogenesis in the gynogenesis of cucumber are discussed, and several key genes regulating embryogenesis were identified. Result Morphological and cytological characterization of embryogenesis in cucumber The embryogenesis early stage, embryo maturation stage and shoot formation stage were observed to discern the cytological changes during embryogenesis by ovary culture in cucumber. Ovaries at 6 h before anthesis were collected and pretreated at 4 o C. The ovaries were sterilized, sliced, inoculated into the medium, incubated for 2 d at 33 o C in the dark and then transferred to 25 o C until plants formed. Interestingly, based on the morphological and cytological observations, critical developmental transitions in embryos were discovered at 2 (T1), 10 (T2), 20 (T3), 30 (T4) and 60 (T5) d of culture. Fresh unpollinated ovaries before culture were named T0 (Fig. 1). At T0, fresh unpollinated ovaries were selected, in which the ovules were obvious (Fig. 1a) and the embryo sacs were mature (Fig. 1a1).At T1, heat shock stress induced the embryogenesis of cucumber (Fig. 1b). The most obvious sign that referred to the development of the embryo was the enlargement of the cells in the embryo sac (Fig. 1b1). The enlargement of microspores have been correlated with embryogenic potential acquisition during induction of androgenesis in many crop species [26~27]. After culturing for 2 d, 4 d and 6 d, the expanded cells formed 2 (Fig. 1b2), 4(Fig. 1b3) and multiple cell structures(Fig. 1b4), respectively. At T2, the explants were cultured for 10 d, the ovules that successfully protruded from the surrounding tissues, and a pro-embryo had formed (Fig. 1c). At T3, after 20 d of culture, the total number of embryos increased continuously in one ovary slice (Fig. 1d). From T3 to T4 (culture for 30 d), the whole embryo had formed, including globular embryos, heart-shaped embryos, torpedo-shaped embryos and cotyledon-shaped embryos (Fig. 1 g). At the same time, the color of the embryonic cells appeared green, indicating that the development of plants from the embryo would begin (Fig. 1e). At T5 (60 d after culture), shoots had formed (Fig. 1f). At T1, ovule enlargement and cell enlargement could be clearly observed, suggesting that the transformation of the development pathway from the gametophyte to the sporophyte was induced by stress and that the stage between T0 and T1 was the key to inducing embryogenesis. After the switching of cell development, embryogenic potential was acquired. Mitosis then occurred continuously to form pro-embryos, globular embryos, heart-shaped embryos, torpedo-shaped embryos and cotyledon-shaped embryos, similar to zygote development, and this process continued until T4. After 6 d of culture, the mature embryos developed into shoots. Therefore, we divided the six time points of ovary culture into three stages: T0 to T1 (embryogenesis early stage), T1 to T4 (embryo maturation stage), and T4 to T5 (shoot formation stage). Transcriptome data were generated via Illumina II HiSeq TM 2000 sequencing, and the statistics of RNA-seq alignment were shown in Table S1. Validation of differentially expressed genes Validation of the Illumina sequencing data and the expression patterns of the DEGs revealed by RNA-seq was performed to examine the expression patterns of 10 DEGs, including 7 genes involved in plant hormone signal transduction and 3 plant-pathogen interaction genes. The results showed that the relative expression levels revealed by RNA-seq and qRT-PCR were closely correlated (Pearson's r = 0.53, 0.77, 0.57, 0.74, 0.79, 0.69, 0.81, 0.81, 0.86, 0.62, 0.81). The fold changes in the qRT-PCR analysis were different from those in the RNA-seq analysis, which might be due to the difference in sensitivity between the qRT-PCR analysis and the RNA-seq technique. However, the qRT-PCR analysis showed that the up- and down-regulation trends of the differential gene expression were consistent with those of the RNA-seq analysis (Fig. 2). Expression of embryogenesis-related genes in cucumber ovary culture At the early stage of embryogenesis, after heat shock stress, ovule enlargement was obvious (Fig. 1b). Cytological observations also showed that one of the cells in the embryo sac expanded (Fig. 1b1), suggesting acquisition of embryogenic potential. The differentially expressed genes (DEGs) at this stage (T0 vs T1) were identified as comprising 3468 up-regulated genes and 3065 down-regulated genes (Fig. 3). We performed a GO enrichment analysis of these genes, and the results showed that the proteins encoded by these genes were assigned to 4 biological processes, 3 molecular functions and 9 cellular components (Table S2). Different genes were involved in different biological processes, and the same gene was involved in a variety of biological processes. The majority of DEGs were involved in the ‘single-organism process’ (GO:0044699) and the ‘oxidation-reduction process’ (GO:0055114). They were mainly distributed among the terms ‘membrane’ (GO:0016020), ‘membrane part’ (GO:0044425), ‘intrinsic component of membrane’ (GO:0031224) and ‘integral component of membrane’ (GO:0016021), having good molecular function in oxidoreductase activity (GO:0016491). Among the DEGs, the expression of related genes involved in the functional classification might play an important role in early embryogenesis. To compare and summarize the results of this stage, we performed a pathway analysis to identify potential target genes. Based on the KEGG database, the pathway enrichment analysis of these genes revealed significant involvement in 8 distinct pathways (Fig. 4): pathways involving plant hormone signal transduction, phenylpropanoid biosynthesis, plant-pathogen interaction, glutathione metabolism, cysteine and methionine metabolism, drug metabolism-cytochrome P450, phenylalanine metabolism and metabolism of xenobiotics by cytochrome P450. Plant hormone signal transduction genes The DEGs were significantly enriched in the plant hormone signal transduction pathway (Ko04075), and a large number of genes encoding cytokinin transfer proteins ( Csa2M373410.1 , Csa6M067360.1 , Csa7M452370.1 ), cytokinin response proteins ( Csa1M006300.1 , Csa3M822100.1 , Csa4M436980.1 , Csa5M223020.1 , Csa5M434550.1 , Csa5M603910.1 , Csa5M623800.1 , Csa6M383530.1 ), ethylene-binding proteins ( CsaUNG009930 ) and ethylene-responsive protein transcription factors ( Csa3M389850.1 ) were detected, which were significantly up-regulated after heat treatment. Moreover, the partially up-regulated DEGs encoding auxin influx carriers, auxin-responsive proteins, SAUR family proteins, abscisic acid receptor proteins, serine/threonine-protein kinases, ABA-responsive element binding factors, and brassinosteroid signal-responsive protein kinases were significantly enriched in 8 pathways, some of which involved tryptophan metabolism as well as zeatin, carotene, and brassinosteroid biosynthesis, suggesting that hormone responses and hormone signal transduction were closely related to stress-induced embryogenesis. The up-regulated genes encoding major enzymes and receptor proteins involved in the plant hormone signal transduction pathway were expressed during the whole stage of ovary culture (Fig. 5). Most of the genes were expressed dynamically after their up-regulated expression during embryogenesis. Only two genes ( Csa6M147590.1 , Csa3M389850.1 ) were up-regulated at the stage of embryo initiation; both almost no expression in the later stages. Csa6M147590.1 and Csa3M389850.1 encode auxin-induced proteins and ethylene-responsive transcription factors, respectively, indicating that they might be closely related to embryogenesis. Stress response proteins related to plant defens e m echanisms DEGs were significantly enriched in the plant-pathogen interaction pathway, which was related to plant defense mechanisms and complex physiological responses to heat shock-induced embryogenesis. In this regulatory pathway, the expression of Pti1 ( Csa7M420160.1 ), HSP90 ( Csa3M183950.1 ), EDS1 ( Csa1M006320.1 ) and other related protein genes were up-regulated. Heat shock protein 90 ( HSP90 ) was a protein that was up-regulated under heat stress. This gene had been highly conserved during evolution and had plant defense functions. Response of genes involved in microtubule organization in cucumber ovary culture In the stage of embryo maturation, the gametophyte development pathway switched to the sporophyte development pathway, in which the cells enlarged and initiating mitosis, forming 2-cell, 4-cell and multicellular pro-embryos after heat stress (Fig. 1b1~b4); afterward, globular embryos, heart-shaped embryos, torpedo-shaped embryos and cotyledon embryos were formed (Fig. 1g). The DEGs analyses of T1 vs T2, T2 vs T3 and T3 vs T4 revealed 480 continuously expressed genes. Additionally, the numbers of uniquely expressed genes between T1 and T2, T2 and T3 and T3 and T4 were 2236, 371 and 1558, respectively (Fig. 6). GO categories of the 480 continuously expressed genes were significantly assigned to 2 biological processes, including cell movement or subcellular component (GO:0006928) and microtubule-based movement (GO:0007018), which were mainly associated with microtubules (GO:0005874), supramolecular fibers (GO:0099512) and polymeric cytoskeletal fibers (GO:0099513). Subsequently, proteins encoded by these genes were classified into 11 functional categories, including protein complex binding (GO:0032403), microtubule binding (GO:0008017), tetrapyrrole binding (GO:0046906), tubulin binding (GO:0015631), microtubule motor activity (GO:0003777), motor activity (GO:0003774), macromolecular complex binding (GO:0044877), heme binding (GO:0020037), oxidoreductase activity (GO:0016491), iron ion binding (GO:0005506) and oxidoreductase activity(GO:0016491), acting on paired donors or the molecular incorporation of oxygen (GO:0016705) (Table S3). This findings indicated that, during certain stages, the cytoskeleton structure was continuously built and maintained by the action of microtubule binding proteins and enzyme modifications, providing the basis for positioning the various organelles and implementation functions, which ensured the orderly activities in various cells in time and space. To determine the involvement of these differentially expressed genes in embryo development, we performed a pathway analysis to identify the potential target genes, and 17 significant pathways were obtained by mapping to the KEGG database (Fig. 7): protein processing in the endoplasmic reticulum, phenylpropanoid biosynthesis, photosynthesis-antenna proteins, limonene and pinene degradation, meiosis-yeast, the estrogen signaling pathway, cell cycle, cell cycle-yeast, diterpenoid biosynthesis, chloroalkane and chloroalkene degradation, carotenoid biosynthesis, progesterone-mediated oocyte maturation, glycerolipid metabolism, histidine metabolism, fatty acid degradation, antigen processing and presentation, and ascorbate and aldarate metabolism. Among them, the pathways of protein processing in the endoplasmic reticulum and phenylpropanoid biosynthesis were the most significant. Expression of main oxidation-reduction and metabolic process-related genes in cucumber ovary culture In the stage of shoot formation, the embryos further differentiated into shoots (Fig. 1f). In total, 3320 genes were differentially expressed between T4 and T5, among which 1383 genes were up-regulated and 1837 were down-regulated in T5 (Fig. 3). Next, we performed an enrichment analysis of the genes using Gene Ontology (GO), which revealed 18 biological processes, 15 molecular functions and 5 cellular components (Table S4). The terms annotated under the biological process category mainly included oxidation-reduction processes (GO: 0055114), the regulation of primary metabolic processes (GO: 0080090), the regulation of cellular metabolic processes (GO: 0031323), and the regulation of macromolecule metabolic processes (GO: 0060255). The molecular function category mainly included the terms oxidoreductase activity (GO: 0016491) and DNA binding (GO: 0003677). The cellular component category mainly included extracellular region (GO: 0005576). The development of the ovary culture was mainly focused on the process of oxidation-reduction and metabolism. The DEGs were subsequently annotated using the KEGG database to identify pathway enrichments. A variety of pathways were found to be significantly enriched (Fig. 8), including those involved in phenylpropanoid biosynthesis; plant hormone signal transduction; stilbenoid, diarylheptanoid and gingerol biosynthesis; metabolism of xenobiotics by cytochrome P450; drug metabolism-cytochrome P450; flavonoid biosynthesis; phenylalanine metabolism; starch and sucrose metabolism; and zeatin biosynthesis. Most of the differentially expressed genes were mainly enriched in phenylpropanoid biosynthesis, plant hormone signal transduction, phenylalanine metabolism, and starch and sucrose metabolism. Pathways involving phenylpropanoid biosynthesis and phenylalanine were important for metabolizing secondary metabolites in plants and were closely related to cell differentiation and pigmentation in plant development. In the plant hormone signal transduction pathway, the expression of an auxin response factor, the gibberellin receptor GID1 , ethylene-responsive transcription factor 1 ( ERF1 ) and cyclin D3 were up-regulated, promoting plant development and maturation, cell enlargement and division as well as promoting the development of stems. In addition, the DEGs involved in starch and sugar metabolism were also significantly enriched and mainly encoding 3-βD glucosidase, T6Ps, GlgB, endoglucanase and alpha-trehalase, which satisfied the requirements of embryo growth and development. Discussion Many previous experiments have shown that early embryogenesis of microspore cultures can be divided into three main phases: acquisition of embryogenic potential, initiation of cell division, and pattern formation [27]. In the present study, the whole process was investigated from acquisition of embryogenic potential to plant regeneration. We divided the ovary culture into three stages. In the early stage of embryogenesis, the acquisition of embryogenic potential by stress (e.g., low or high temperature and hormone induction) was observed together with the repression of gametophytic development, leading to the dedifferentiation of cells. In the embryo maturation stage, cell divisions gave rise to the formation of pro-embryos (cell clumps), globular embryos, heart-shaped embryos, torpedo-shaped embryos and cotyledon embryos (mature embryos). In the stage of shoot formation, mature embryos developed into shoots. There are active molecular events that regulate embryo development at different stages of development. Embryogenesis-related genes are expressed in cucumber ovary culture The embryogenic ability and transformation of somatic cells were regulated by various plant hormones, such as auxin, abscisic acid, cytokinin and ethylene [28]. Moreover, the synthesis of auxin and ethylene increased significantly after ovary pollination and fertilization [29]. In normal growth, mitotic asymmetry of fertilized eggs was found in the gametophyte development pathway. Previous work using molecular and genetic approaches has identified auxin as a critical signal for the proper development of the asymmetric structure of the female gametophyte in Arabidopsis [30]. In in vitro culture, a certain stress condition is needed to block the development direction of the original gametophyte, followed by turning into the direction of the sporophyte and carrying on the symmetry splitting to eventually lead to embryogenesis [31]. In ovary culture of cucumber, heat shock pretreatment, silver nitrate, genotype and hormone combination factors could play key roles in embryo and callus production independently and simultaneously [32]. In the ovary culture of linseed, cultivar, combination of growth regulators, type of carbohydrates and their interaction significantly influenced callus induction and shoot formation frequency, and a relatively high shoot regeneration frequency was obtained when the medium was supplemented with TDZ and NAA [33]. These results implied that the expression of specific genes in embryogenesis was activated by exogenous hormone regulation, which laid the foundation for DNA replication in the cell division stage. In our study, a large number of DEGs were enriched in the plant hormone signal transduction pathway during the stress process of embryogenesis. The pathway annotation obtained by KEGG analysis showed that phytohormones such as cytokinin and ethylene were significantly up-regulated in multiple biosynthesis and metabolic pathways, while up-regulated and down-regulated genes related to auxin and abscisic acid both existed. Generally, high-level endogenous cytokinin and low-level endogenous ethylene are beneficial for the acquisition of embryogenic ability. Embryogenesis is the comprehensive performance of interactions among different hormones, and all kinds of endogenous hormones show dynamic changes. The up-regulated expression of the main enzymes and receptor proteins in the plant hormone signal transduction pathway could promote embryogenesis, indicating that the high expression of related genes may play important roles in the process of switching from the gametophyte to the sporophyte development pathway. In this pathway, we found two genes related to hormone regulation that were highly expressed only at the stage of early embryo development ( Csa6M147590.1 and Csa3M389850.1 );moreover, the functional categories of their specific functions and the location of the regulatory pathways were identified, thus laying an important theoretical basis for the elucidation of the mechanism. Heat shock protein 90 ( HSP90 ) is an evolutionarily conserved molecular chaperone induced by abiotic stress. HSP90 plays an important role in cell cycle control, genomic silencing, protein transduction and signal transduction[20]. The HSP90 chaperone is involved in maintaining phenotypic plasticity and developmental stability [34-35]. Previous studies have shown that HSP genes were specifically expressed in the spore stress-induction process of many crop species [36-38]. HSPs had been described as genes associated with the reprogramming phase and the early stages of embryogenesis [39]. After heat shock stress treatment, we found that 26 up-regulated genes encoded proteins and transcription factors related to heat shock, and HSP90 ( Csa3M183950.1 ) was involved in the regulation of the hypersensitive response in the plant and pathogen interaction pathway. The expression of these genes might play a regulatory role in the process of embryogenesis. The BBM gene was the first key gene to be isolated in the process of spore cell division; which was first expressed in zygotic embryogenesis and microspore embryogenesis [14, 40]. We found that two BBM genes ( Csa3M827310.1 , Csa3M827320.1 ) were up-regulated at the beginning of embryogenesis and were subsequently down-regulated. In addition, AGL15 a member of the MADS family, was also considered a regulatory protein that acts at the start of cell division [13]. We found that the AGL gene ( Csa3M258140.1 ) was up-regulated in our study, indicating that all three of these genes play important roles in embryogenesis. We considered that BBM, HSP90 and AGL might be the critical genes involved in the induction of embryogenesis by ovary culture in cucumber. Microtubule organization genes play an important role in the embryogenesis of cucumber ovary culture There is a similar process between gametophyte embryogenesis and zygotic embryogenesis, and the only difference is that the former occurs in sex cells, while the latter occurs in fertilized egg cells. The process of gametophyte embryogenesis development is similar to that of zygote differentiation into embryos. Microtubules play a transport role in cells, participating in the construction of cell walls and promoting cell differentiation and division. Electron microscope-based observations of carrot culture showed that the appearance of microtubules was accompanied with the formation of somatic embryos [41]. The different patterns of microtubule organization in the cells of the mature embryo sac reflect their structural adaptations for their future function [42]. In our study, after embryogenic potential was acquired, the different embryo shapes formed through active cell metabolism and rapid cell division. GO analysis showed that in different embryo periods continuously expressed genes participated in microtubule-based movement and cell or subcellular component processes, providing a good function of protein complex binding, microtubule binding, tetrapyrrole binding, tubulin binding and other microtubule activities; these genes mainly participated in protein processing of the endoplasmic reticulum and phenylpropanoid biosynthesis. We considered that the cellular reprogramming and morphological changes in embryos were controlled by microtubule organization genes. Expression of main oxidation-reduction and metabolic process-related genes in cucumber ovary culture Mature embryos are further cultured to form shoots, and at this stage, many biological processes are involved. GO analysis showed that differentially expressed genes participated in the oxidation-reduction process, a variety of metabolic processes, biosynthesis processes and the response to auxin, together with the molecular functions of oxidoreductase activity, binding capacity, catalytic activity, transporter activity, transcription factor activity and so on. Among them, the regulation of transcription, DNA template (GO: 0006355), the regulation of RNA biosynthetic processes (GO: 2001141), the regulation of RNA metabolic processes (GO: 0051252), DNA binding (GO: 0003677), transcription factor activity (GO: 0003700), and sequence-specific DNA binding (GO: 0043565) lay a foundation for cell division and differentiation. In addition, two processes, cellular response to auxin stimulus (GO: 0071365) and response to auxin (GO: 0009733), might promote the embryo to form shoots. At the stage of shoot formation, a large number of transcription factors were expressed. The first constitutes transcription factors that contain a B3 domain [43]. These transcription factors can encode regulatory proteins involved in the embryonic development process, maintaining embryo development during late embryonic development [43, 44]. In our study, total of 6 transcription factors( Csa2M359980.1 , Csa2M292240.1 , Csa5M608380.1 , Csa6M489980.1 , Csa6M489940.1 and CsaUNG031640 ) that contain a B3 domain were up-regulated . The second was the AIL gene from the AP2/ERF family, which is involved in key developmental processes throughout the whole plant life cycle. Some genes in the AP2/ERF gene family are expressed in many tissues and participate in many plant development processes, such as embryogenesis and shoot development [14, 45). In addition, we found that the expression levels of 9 genes ( Csa1M423190.1 , Csa3M114480.1 , Csa3M652380.1 , Csa3M114470.1 , Csa4M644740.1 , Csa5M175970.1 , Csa5M608380.1 , Csa6M496390.1 ,and CsaUNG031640 ) in the AP2/ERF family were up-regulated. Furthermore, WUS homologous domain proteins not only alter the cell fate of the shoot and flower meristem but also promote the development of somatic embryos into seedlings. The function of WUS proteins has no direct connection with the characteristics of the embryo but will alter the development state of the tissue by maintaining cells in an undifferentiated state in response to different stimuli [46]. The ability to transform the vegetative growth phase to the embryonic stage by WUS , as well as to eventually form somatic cells, indicates that this homologous domain protein also plays an important role in embryo maturation in addition to its role in the development of the meristem [47]. During this stage, two genes encoding WUS homologous domain proteins ( Csa6M301060.1 ) and ( Csa6M505860.1 ) were up-regulated. We considered that those homologous domain proteins might directly or indirectly regulate shoot formation. Conclusion Several studies have reported on the embryogenesis of ovary culture in cucumbers. However, the mechanism that drives the process from embryogenic acquisition to the formation of embryos and plant regeneration is not well understood. In this study, we explored embryogenesis mechanism of ovary culture in cucumber. Inducing cucumber embryogenesis could be divided into three stages: early embryogenesis, embryo maturation (from proembryos to cotyledon embryos) and shoot formation (Fig. 9a). The early stage of embryogenesis was the turning point for the formation of embryos, which had experienced dedifferentiation and loss of photosynthetic capacity, requiring the provision of exogenous nutrients and carbon sources such as plant hormones and sucrose in the media. Some physical factors, such as temperature, light quality, photoperiod and the presence of specific hormones, could affect the ability of the embryo sac to adapt to these conditions and survive the developmental transition. These dynamic changes were helpful for cell physiology reprogramming, metabolic alterations, revival of cell division and differentiation, morphogenesis and so on. Therefore, the cells of the embryo sac began to divide by dedifferentiation and started to form cell clumps (Fig. 1b1~b4). Although there was no obvious change in appearance, the metabolism of some macromolecules in the cells changed markedly. A large number of hormone-related genes, cell protection-related genes, and some unique protein kinases and transcription factors ( Csa6M147590.1 and Csa3M389850.1 ) were expressed only at the stage of early embryo development . The reported embryogenesis-related genes BBM , HSP90 and AGL were also actively expressed during this stage (Fig. 9b). In the stage of embryo maturation, GO analysis showed continuously expressed genes participated in microtubule-based movement, movement of the cell or subcellular component processes, giving a good function of protein complex binding, microtubule binding, tetrapyrrole binding, tubulin binding and other microtubule activities, which are involved in protein processing in the endoplasmic reticulum and phenylpropanoid biosynthesis(Fig. 9b). Shoot formation was regulated by 6 transcription factors that contain a B3 domain, total of 9 genes in the AP2/ERF family and 2 genes encoding WUS homologous domain proteins (Fig. 9b). Methods Plant materials ‘SG033’ (F 1 cultivar of Kunming Huaxing seed Industry Co., Ltd., Kunming, Yunnan) was used for ovary culture; this material had a high embryo rate after a large number of screenings. Cucumber plants from southern China were cultivated and stored in our lab. The plants were grown in a greenhouse with a 12 h photoperiod, a mean daily air temperature of 25/15 o C (day/night), a relative humidity of 85% and a photosynthetic photo flux density of 500 μmol·m -2 ·s -1 at the Horticulture Institute of Guizhou Academy of Agricultural Sciences, Guizhou, China. Ovaries were collected from ‘SG033’ for one month after the first anthesis flower appeared. Ovary culture Ovaries were harvested at 6 h before anthesis, after treated for 24 h in a refrigerator at 4 o C, which papilloma was removed. The ovaries were then sterilized in 75% ethanol for 30-60 s and rinsed in sterile distilled water 3 times, followed by soaking in a 0.5% sodium hypochlorite for 20 min and then rinsing in sterile distilled water 3 times. The ovaries were cut into small round slices of 2 mm under sterile conditions and then placed on 30 mL of solid media in cylindrical flasks. In all experiments, MS medium was used as a basal medium; it was supplemented with 0.06 mg·L -1 TDZ and 3% (w/v) sucrose and was solidified with 7% (w/v) agar. The pH of the medium was adjusted to 5.9 before autoclaving at 116 o C and 1.1 kg/cm 2 for 30 min. All materials were treated with high temperature of 33 o C for 2 d, after which recovery growth was allowed until shoot formation occurred at 25 o C under a 16/8 h (light/dark) photoperiod with a 4000 lx light intensity. RNA isolation The materials were harvested separately at the following time points: T0 (the ovules were cultured for 0 d, fresh unpollinated ovaries), 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 culture for 6 d). These materials were stored at -80 o C for RNA extraction, and the remaining materials were maintained in culture to observe plantlet regeneration. Furthermore, all experimental procedures, such as culture medium replacement and sample collection, were performed under similar conditions to minimize possible circadian effects. Ovules at T0 and T1 were extracted by hand under a stereoscopic microscope, and embryo-like structures with integuments were selected directly at the time points from T2 to T5. As described above, samples were collected from six time points of ovary culture for RNA-seq analysis. For each sample, the ovules or embryo-like structures were immersed in liquid nitrogen in a mortar and ground into powder (three biological replicates per sample). Total RNA was isolated using TRIzol (Invitrogen), and DNase I (Fermentas) digestion was performed for 30 min at 25 o C to remove DNA, according to the manufacturer’s instructions. The integrity and quality of the total RNA were checked using a NanoDrop 1000 spectrophotometer and formaldehyde-agarose gel electrophoresis. RNA was used only when the OD 260 : OD 280 ratio was above 1.8. Transcriptome analysis Transcriptomes were analyzed by an Illumina HiSeq TM 2000. The reads with more than 10% unknown bases and those of low quality were removed from the raw sequencing data using the software Illumina GA Pipeline (v1.6). The filtered reads were aligned to the cucumber Chinese long v2 genome using TopHat2 ( http://tophat.cbcb.umd.edu/ ) software. The reference genome and gene database information were obtained from the public website: http://cmb.bnu. edu.cn/Cucumis_sativus_v20/ [42]. Differential expression analysis was determined using DESeq (version 1.18.0). The resulting P values were adjusted by using Benjamini and Hochberg’s approach for controlling the false discovery rate. Genes with an adjusted P value of <0.05 and an absolute value of |log2fold change≥1| according to DESeq were considered differentially expressed. GO annotation was performed using Blast2GO software (GO association was performed by a BLASTX search against the NCBI NR database). GO enrichment analysis of differentially expressed genes (DEGs) was then performed via the BiNGO plugin for Cytoscape. Overrepresented GO terms were identified using a hypergeometric test with a significance threshold of 0.05 after the Benjamini-Hochberg FDR correction. KEGG enrichment analysis of the differentially expressed genes was performed using the KOBAS (2.0) [48] software. A ‘ P value ≤ 0.05’ was used as the threshold to judge the significantly enriched pathways of the differentially expressed genes. Real-time quantitative PCR Ten DEGs were selected from all the DEGs for analysis using real-time quantitative PCR (qRT-PCR). For this step, the total RNA was treated with the DNase I enzyme and subsequently converted to single-strand cDNA by applying the GoScript TM Reverse Transcription System (Promega) according to the manufacturer’s protocol. The specific gene primers were designed using Primer 5.0 according to the cDNA sequences in Table S5. Actin was used as the internal control to normalize small differences in template quantities. Quantitative real-time PCR was performed with GoTaq qPCR Master Mix (Promega) in a Bio-Rad iQ1 real-time PCR system (Bio-Rad), and the collected data were analyzed using the 2 -ΔΔCt method. Abbreviations DH, Double haploid; DEG, Differentially expressed genes; qRT-PCR, Quantitative real-time PCR; SSR, Simple sequence repeats; FDR, False discovery rate; AGL15 , Agamous-like15; AP2/ERF, Ethylene-responsive element-binding factor; BBM , Baby boom; HSP , Heat shock proteins; Pti 1 , Pto-interacting protein 1; AIL ,Aintegumenta-like; LEC, Leafy Cotyledon; GID1 , gibberellin receptor, ERF1 , ethylene-responsive transcription factor1; WUS , WUSCHEL; TDZ, thidiazuron; 3-βD glucosidase, 3-beta-D-glucosidase; T6Ps, Trehalose 6-phosphate synthase; GlgB, 1,4-alpha-glucan branching enzyme. Declarations Ethics approval and consent to participate Not applicable. Consent for publication Not applicable. Availability of data and materials The data sets supporting the results of this article are included within the article and its additional files. Competing Interests The authors declare that they have no conflict of interest. Funding The financial support from National Natural Science Foundation of China (No. 31460515), Jiangsu Agricultural Innovation of New Cultivars (No. PZCZ201719)，Guizhou Science and Technology Support Plan （NO.20192259）and Special Fund for Agro-Scientific Research of Guizhou Academy of Agricultural Sciences (No. 2015-19) are deeply appreciated. Guizhou Construction of Genetics and Breeding Laboratory of Vegetable Industry System(2019-0102). Author contributions JFC guide and designed the research. YD , BT, WYF, LT, LZ , XZ and JL performed the experiment, analyzed, interpreted date, and rote the manuscript. All authors reviewed and contributed to draft the manuscript. All authors read and approved the final manuscript. Acknowledgements The authors would like to acknowledge the Dr Tan (Institute of Horticulture, Guizhou Academy of Agricultural Sciences, Guiyang , China) for critical reading of the manuscript , Jiaji Zhang( Personal Biotechnology Corporation, Shanghai, China) for their efforts in bioinformatics analysis and Springer Nature Author Services for writing services. References Guha S, Maheshwari SC. In vitro production of embryos from anthers of datura. nature. 1964;204:497. Malik MR, Wang F, Dirpaul JM, Zhou N, Polowick PL, Ferrie, AMR, Krochko JE. Transcript profiling and identification of molecular markers for early microspore embryogenesis in Brassica napus . Plant Physiology. 2007;144:134-154. Bhat JG, Murthy HN. Factors affecting in-vitro gynogenic haploid production in niger ( Guizotia abyssinica (L. f.) Cass.). Plant Growth Regulation. 2007;52:241-248. Metwally EI , Moustafa SA , El-Sawy BI , Haroun SA, Shalaby TA. Production of haploid plants from in vitro culture of unpollinated ovules of Cucurbita pepo . Plant Cell Tissue and Organ Culture. 1998;52:117-121. Gémes-Juhász A, Balogh P, Ferenczy AZ, Kristóf. Effect of optimal stage of female gametophyte and heat treatment on in vitro gynogenesis induction in cucumber ( Cucumis sativus L.). Plant Cell Reports. 2002;21:105-111. Halit Y, Nebahat S. A new method for haploid muskmelon ( Cucumis melo L.) dihaploidization. Scientia Horticulturae. 2003;98:277–283. Shalaby TA. Factors affecting haploid induction through in vitro gynogenesis in summer squash ( Cucurbita pepo L.). Scientia Horticulturae. 2007;115:1-6. Diao WP, Jia YY, Song H, Zhang XQ, Lou QF, Chen JF. Efficient embryo induction in cucumber ovary culture and homozygous identification of the regenetants using SSR markers. Scientia Horticulturae. 2009;119:246-251. Li JW, Si SW, Cheng JY, Li JX, Liu JQ. Thidiazuron and silver nitrate enhanced gynogenesis of unfertilized ovule cultures of Cucumis sativus . Biologia Plantarum. 2013;57:164-168. Plapung P, Khumsukdee S, Smitamana P. Development of cucumber lines resistant to cucumber mosaic virus by ovule culture. International Journal of Agricultural Technology. 2014;733-741. Tantasawat PA, Sorntip A, Poolsawat O, Chaowiset W, Pornbungkerd P. Evaluation of factors affecting embryo-like structure and callus formation in unpollinated ovary culture of cucumber( Cucumis sativus ). International Journal of Agriculture & Biology. 2015;17:613-618. Touraev A, Vicente O, Heberle-Bors E. Initiation of microspore embryogenesis by stress. Trends in Plant Science. 1997;2:297–302. Perry SE, Fernandez DE. The MADS-domain protein AGAMOUS-like 15 accumulates in embryonic tissues with diverse origins. Plant Physiology. 1999;120:121-129. Boutilier K, Offringa R, Sharma VK, Kieft H, Ouellet T, Zhang L, Hattori J, Liu CM, van Lammeren AAM, Miki BLA, Custers JBM, van Looleren Campagne MM. ectopic expression of BABY BOOM triggers a conversion from vegetative to embryonic growth. Plant Cell. 2002;14:1737-1749. Maraschin SF, Priester WD, Spaink HP, Spaink, Wang M. Androgenic switch: an example of plant embryogenesis from the male gametophyte perspective. Journal of Experimental Botany. 2005;56:1711-1726. Hosp J, Tashpulatov A, Roessner U, Barsova E, Katholnigg H, Steinborn R, Melikant B, Lukyanov S, Heberle-Bors E, Touraev A. Transcriptional and metabolic profiles of stress-induced, embryogenic tobacco microspores. Plant Molecular Biology. 2007;63:137-149. Rosa Angélica SánchezDíaz, Ana María Castillo, María Pilar Vallés. Microspore embryogenesis in wheat: new marker genes for early, middle and late stages of embryo development. Plant Reproduction. 2013;26:287-296. Irikova T, Grozeva S, Denev I. Identification of BABY BOOM and LEAFY COTYLEDON genes in sweet pepper ( Capsicum annuum L.) genome by their partial gene sequences. Plant Growth Regulation. 2012; 67:191-198. Bélanger S, Marchand S, Jacques PÉ, Meyers B, Belzile F. Differential expression profiling of microspores during the early stages of isolated microspore culture using the responsive barley cultivar gobernadora. Genes Genomes Genetics.2018;8:1603-1614 Zhang J , Li J , Liu B , Zhang L, Chen J, Lu M. Genome-wide analysis of the populus HSP90 gene family reveals differential expression patterns localization and heat stress responses. BMC Genomics. 2013;14:532-546. Lv J , Qi J , Shi Q , Shen D, Zhang S, Shao G, Li H, Sun ZY, Weng Y, Gu XF, Li XX, Zhu XG, Zhang JZ, van Treuren R, van Dooijeweert W, Zhang ZH, Huang S. Genetic diversity and population structure of cucumber ( Cucumis sativus L.). PLoS ONE. 2012;7:e46919. Qing ZW, Xin M, Liu D, Zhang YJ, Wang XG, Zhou XY. A new cucumber F 1 hybrid 'dongnong 812'. Chinese Vegetables. 2018;5:71-73.(in Chinese) Golabadi M, Golkar P, Eghtedary AR. Assessment of genetic variation in cucumber ( Cucumis sativus L.) genotypes. European Journal of Experimental Biology. 2012;2 :1382-1388. Wang X, Bao K, Reddy UK, Bai Y, Hammar SA, Jiao C, Wehner TC, Ramírez-Madera AO, Weng Y, Grumet R, Fei Z. The USDA cucumber ( Cucumis sativus L.) collection: genetic diversity, population structure, genome-wide association studies, and core collection development. Horticulture Research. 2018;5:64-67 Chen JF, Cui L, Malik AA, Kere George Mbira. In vitro haploid and dihaploid production via unfertilized ovule culture. Plant Cell Tissue and Organ Culture. 2011; 104:311-319. Hoekstra S, Hoekstra IR, Hoekstra RA, Hoekstra E. The Interaction of 2, 4-D application and mannitol pretreatment in anther and microspore culture of Hordeum vulgare L. cv. Igri. Journal of Plant Physiology. 1996;148:696-700. Maraschin, SF. Androgenic switch: an example of plant embryogenesis from the male gametophyte perspective. Journal of Experimental Botany. 2005;56:1711-1726. Ikeda-Iwai M, Umehara M, Satoh S, Kamada H. Stress‐ induced somatic embryogenesis in vegetative tissues of Arabidopsis thaliana . Plant Journal. 2003;34:107-114. Mól R, Filek M, Machackova I, Matthys-Rochon E. Ethylene synthesis and auxin augmentation in pistil tissues are important for egg cell differentiation after pollination in maize. Plant & Cell Physiology. 2004;45:1396-1405. Pagnussat GC, Alandetesaez M, Bowman JL, Sundaresan v. Auxin-dependent patterning and gamete specification in the Arabidopsis female gametophyte. Science. 2009;324:1684-1689. Fan Z, Armstrong KC, Keller WA. Development of microspores in vivo , and in vitro , in Brassica napus L. Protoplasma. 1988;147:191-199. Golabadi M, Ghanbari Y, Ghanbari S, Keighobadi K. Embryo and callus induction by different factors in ovary culture of cucumber. Journal of Applied Botany & Food Quality. 2017;90:68-75. Aušra B, Natalija B, Ramunė M. Genotypic and exogenous factors affecting linseed ovary culture. Zemdirbyste. 2017;104:243-248. Hahn A, Bublak D, Schleiff E, Scharf KD. Crosstalk between Hsp90 and Hsp70 Chaperones and Heat Stress Transcription Factors in Tomato. Plant Cell. 2011; 23:741-755. Sangster TA, Queitsch C. The HSP90 chaperone complex, an emerging force in plant development and phenotypic plasticity. Current Opinion in Plant Biology. 2005;8:86-92. Zarsky V, Garrido D, Eller N, Eller N, Tupy J, Vicente O, Schoffl F. The expression of a small heat shock gene is activated during induction of tobacco pollen embryogenesis by starvation. Plant Cell & Environment. 2010;18:139-147. Binarova P, Hause G, Cenklová V, Cordewener JHG, Looderen Campagne MM. A short severe heat shock is required to induce embryogenesis in late bicellular pollen of Brassica napus L.. Sexual Plant Reproduction. 1997;10:200-208. Sabehat A, Weiss D. Expression of small heat-shock proteins at low temperatures: A possible role in protecting against chilling injuries. Plant Physiology. 1998;117:651-8. Maraschin SF. Androgenic switch: an example of plant embryogenesis from the male gametophyte perspective. Journal of Experimental Botany. 2005; 56:1711-1726. Pechan PM, Bartels D, Brown DC, Schell J. Messenger-RNA and protein changes associated with induction of Brassica , microspore embryogenesis. Planta. 1991;184:161-165. Halperin W, Jensen WA. Ultrastructural changes during growth and embryogenesis in carrot cell cultures. Journal of Ultrastruct Research. 1967;18:428-443. Huang BQ, Sheridan WF. Female Gametophyte Development in maize: microtubular organization and embryo sac polarity. Plant Cell. 1994;6:845-861. Stone SL, Kwong LW, Yee KM, Yee KM, Pelletier J, Lepiniec L, Fischer RL. LEAFY COTYLEDON2 encodes a B3 domain transcription factor that induces embryo development. Proceedings of the National Academy of Sciences. 2001;98:11806-11811. Lotan T, Ohto M, Yee K M, West MA, Lo R, Kwong RW, Yamagishi, Fischer RL, Goldberg RB, Harada JJ. Arabidopsis LEAFY COTYLEDON1 is sufficient to induce embryo development in vegetative cells. Cell. 1998;93:1195-1205. Riechmann JL, Meyerowitz EM. The AP2/EREBP family of plant transcription factors. Biological Chemistry. 1998;379:633-646. Mayer KF, Schoof H, Haecker A, Lenhard M, Laux T. Role of WUSCHEL in Regulating Stem Cell Fate in the Arabidopsis Shoot Meristem. Cell. 1998;95:805-815. Palovaara J, Hakman I. Conifer WOX-related homeodomain transcription factors, developmental consideration and expression dynamic of WOX2 during Picea abies somatic embryogenesis. Plant Molecular Biology. 2008;66:533-549. Mao X, Cai T, Olyarchuk JG, Wei L. Automated genome annotation and pathway identification using the KEGG Orthology (KO) as a controlled vocabulary. Bioinformatics. 2005;21:3787–3793. Additional file legend Table S1. Statistics of RNA-Seq alignment Table S2. GO classification of common expressed genes in gynogenesis stage. Table S3. GO classification of common expressed genes from embryogenesis to the mature stage. Table S4. GO classification of common expressed genes in organogenesis stage. Table S5. Primers for Real-time quantitative PCR Supplementary Files supplement1.docx Cite Share Download PDF Status: Posted Version 1 posted You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-1316","acceptedTermsAndConditions":true,"allowDirectSubmit":true,"archivedVersions":[],"articleType":"Research article","associatedPublications":[],"authors":[{"id":96200,"identity":"188b058d-ea55-47d0-9ec0-3e7e77ca2341","order_by":1,"name":"Ying Deng","email":"","orcid":"","institution":"Nanjing Agricultural University;Institute of Horticulture, Guizhou Academy of Agricultural Sciences","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Ying","middleName":"","lastName":"Deng","suffix":""},{"id":96201,"identity":"4e3a95ba-853a-4a88-b894-18c054f5b834","order_by":2,"name":"Wenyuan Fu","email":"","orcid":"","institution":"Institute of Horticulture,Guizhou Academy of Agricultural Sciences","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Wenyuan","middleName":"","lastName":"Fu","suffix":""},{"id":96202,"identity":"f8c28d8c-9c33-4ac5-aba2-d043d9517d45","order_by":3,"name":"Bing Tang","email":"","orcid":"","institution":"Institute of Horticulture,Guizhou Academy of Agricultural Sciences","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Bing","middleName":"","lastName":"Tang","suffix":""},{"id":96203,"identity":"49d538f5-24e8-4720-a603-fc5f3a2a5b7f","order_by":4,"name":"Lian Tao","email":"","orcid":"","institution":"Institute of Horticulture, Guizhou Academy of Agricultural Sciences","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Lian","middleName":"","lastName":"Tao","suffix":""},{"id":96204,"identity":"8f39b6c0-4646-46bd-b43e-3970fb3af0af","order_by":5,"name":"Lu Zhang","email":"","orcid":"","institution":"Nanjing Agricultural University","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Lu","middleName":"","lastName":"Zhang","suffix":""},{"id":96205,"identity":"3ff8f768-a94c-49d1-9b89-1f74a9631006","order_by":6,"name":"Xia Zhou","email":"","orcid":"","institution":"Nanjing Agricultural University","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Xia","middleName":"","lastName":"Zhou","suffix":""},{"id":96206,"identity":"889a606c-6b09-493e-aa60-f44b908f2172","order_by":7,"name":"Ji LI","email":"","orcid":"","institution":"Nanjing Agriculture University","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Ji","middleName":"","lastName":"LI","suffix":""},{"id":96207,"identity":"3bc6713f-eaf5-48d8-aa14-ea2be7a89cf3","order_by":8,"name":"Jinfeng Chen","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAA0UlEQVRIiWNgGAWjYBACxuYDbEDKhoGxAcRlI0ZLWwJIWRoJWhjYwFoOw3hE6GBuY3724OOO83nM7WcMGD6UHWbgn91AyGFs5oYzz9wuZuzJMWCcce4wg8SdAwS0zO9hk+Ztu53Y2JBjwMzbdpjBQCKBkC08bNJ/284lNva/MWD+S7QWxrYDiY0zgLYwEqeFzUyyty0ZqOVZwcGec+k8EjcIaDEEhpjEzza7xI39yRsf/CizluOfQUhLAxLjAJDmwa8eCOQxGKNgFIyCUTAK0AEAfP9CZBznivkAAAAASUVORK5CYII=","orcid":"https://orcid.org/0000-0003-0127-3070","institution":"Nanjing Agricultural University","correspondingAuthor":true,"submittingAuthor":false,"prefix":"","firstName":"Jinfeng","middleName":"","lastName":"Chen","suffix":""}],"badges":[],"createdAt":"2019-06-12 14:07:49","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.2.10287/v1","doiUrl":"https://doi.org/10.21203/rs.2.10287/v1","draftVersion":[],"editorialEvents":[],"editorialNote":"","failedWorkflow":false,"files":[{"id":2614188,"identity":"babb62d3-e431-4c23-be99-519691e7009f","added_by":"399a65fa-c8d4-440b-992c-2804352e6ac3","created_at":"2020-09-25 20:54:05","extension":"jpg","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":395068,"visible":true,"origin":"","legend":"Embryogenic process of ovary culture in cucumber \na-g observation under stereoscopic microscope. a selected fresh unpollinated ovaries(0 d, T0). b ovule were treated with the high temperature of 33 oC for 2 d; c the embryo cultured for10 d, the ovule enlargement and embryo initiation. d embryos of culturing 10 d, an average of 6 embryoids in each slice of ovary; e embryos of culturing 30 d, embryo maturation and cotyledon-embryo formation. f culturing for 60 d, the embryo differentiated into shoot. g globular-embryo switch into cotyledon shape-embryo. Histological observations: the histological sections were stained with Delafield’s hematoxylin. a1 multicellular embryo sac before culturing, b1 one of the cells expands in the embryo sac, 2 d of culture, b2 cell mitosis, 4 d of culture. b3 continue to be divided into four cells, 6 d of culture. b4 cell clumps structure, 8 d of culture. G-e globular embryo, H-se heart-shaped embryo, T-se torpedo-shaped embryo, C-se cotyledon-shaped embryo, E sembryo sac. The solid white lines represent 2mm, and the black lines represent 200μm. Ovule was marked with a white arrow.","description":"","filename":"figure1.jpg","url":"https://assets-eu.researchsquare.com/files/rs-1316/v1/figure_1.jpg"},{"id":2614182,"identity":"0d9c5d63-7559-4c68-b7f8-25e61f6ad1e2","added_by":"399a65fa-c8d4-440b-992c-2804352e6ac3","created_at":"2020-09-25 20:54:05","extension":"png","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":131083,"visible":true,"origin":"","legend":"Validation of differentially expressed genes and correlation between RNA-seq and qRT-PCR \nqRT-PCR of 7 up-regulated genes involved in plant hormone signal transduction and 3 plant-pathogen interaction genes in whole ovary culture were analyzed. Close correlations (Person’s r=0.53, 0.77, 0.57, 0.74, 0.79, 0.69, 0.81, 0.86, 0.62, 0.81) were observed between relative expression levels in RNA-seq and qRT-PCR, validating the RNA-Seq methodology described here for quantitative analysis of the cucumber transcriptome. Three independent experimental replicates were analyzed for each sample, and data were showed as mean ± SE (n = 3).","description":"","filename":"figure2.png","url":"https://assets-eu.researchsquare.com/files/rs-1316/v1/figure_2.png"},{"id":2614184,"identity":"0ce2da6f-34d8-425e-9035-fecbe454c9af","added_by":"399a65fa-c8d4-440b-992c-2804352e6ac3","created_at":"2020-09-25 20:54:05","extension":"jpg","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":32809,"visible":true,"origin":"","legend":"The total number of up-regulated and down-regulated genes\nRed color represents up-regulated, and green represents down-regulated","description":"","filename":"figure3.jpg","url":"https://assets-eu.researchsquare.com/files/rs-1316/v1/figure_3.jpg"},{"id":2614180,"identity":"fdc49307-7b8e-4635-a63a-4d10be0628e1","added_by":"399a65fa-c8d4-440b-992c-2804352e6ac3","created_at":"2020-09-25 20:54:05","extension":"jpg","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":166648,"visible":true,"origin":"","legend":"KEGG pathway enrichment analysis based on the differentially expressed genes in early embryo development \nThe longitudinal coordinates was KEGG Pathway, and the abscissa was richfactor; the size of the dots in the graph represented the number of differentially expressed genes annotated to the path, and the color represented the significant P value of the path. The 20 most significant paths were shown in the picture.","description":"","filename":"figure4.jpg","url":"https://assets-eu.researchsquare.com/files/rs-1316/v1/figure_4.jpg"},{"id":2614185,"identity":"eaa3d411-365c-4797-b1dd-8b42ad8123ba","added_by":"399a65fa-c8d4-440b-992c-2804352e6ac3","created_at":"2020-09-25 20:54:05","extension":"png","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":210424,"visible":true,"origin":"","legend":"The matrix graph of up-regulated expression genes in the plant hormone signal transduction pathway in six time points\nThe absolute values of the relative expression of up-regulated genes were indicated by the depth of color. The abscissa of the matrix graph showed the six time points. The vertical ordinate of the matrix graph showed the name of the up-regulated genes. The up-regulated genes were arrayed according to the hormonal signal transduction flows, which were classified into six groups: auxin, cytokinin, abscisic acid, ethylene, brassinosteroid, salicylic acid.","description":"","filename":"figure5.png","url":"https://assets-eu.researchsquare.com/files/rs-1316/v1/figure_5.png"},{"id":2614179,"identity":"67be1e1f-a9f4-4d6c-83e6-d7d3e67aa435","added_by":"399a65fa-c8d4-440b-992c-2804352e6ac3","created_at":"2020-09-25 20:54:05","extension":"jpg","order_by":6,"title":"Figure 6","display":"","copyAsset":false,"role":"figure","size":38575,"visible":true,"origin":"","legend":"Venn diagram showed differentially expressed genes in the stage from embryogenesis to the mature stage \nThe DEG sets (T1 vs T2, T2 vs T3 and T3 vs T4) described in Fig. 3 were analyzed by using the Venn method. The numbers marked in the diagram indicated the number of continue express genes significantly expressed among the three DEGs sets.","description":"","filename":"figure6.jpg","url":"https://assets-eu.researchsquare.com/files/rs-1316/v1/figure_6.jpg"},{"id":2614186,"identity":"a3d617fd-d6d9-44ac-a7db-ae63f4df129f","added_by":"399a65fa-c8d4-440b-992c-2804352e6ac3","created_at":"2020-09-25 20:54:05","extension":"jpg","order_by":7,"title":"Figure 7","display":"","copyAsset":false,"role":"figure","size":166404,"visible":true,"origin":"","legend":"KEGG pathway enrichment analysis based on the differentially expressed genes from embryogenesis to the mature stage\nThe longitudinal coordinates was KEGG pathway; the abscissa was rich factor; and the size of the dots in the graph represented the number of differentially expressed genes that were annotated to the path; the color represented the significant P value of the path. The 20 most significant paths are shown in the picture.","description":"","filename":"figure7.jpg","url":"https://assets-eu.researchsquare.com/files/rs-1316/v1/figure_7.jpg"},{"id":2614183,"identity":"db9c21ad-d6d8-4f3d-ba40-3c619721f0a2","added_by":"399a65fa-c8d4-440b-992c-2804352e6ac3","created_at":"2020-09-25 20:54:05","extension":"jpg","order_by":8,"title":"Figure 8","display":"","copyAsset":false,"role":"figure","size":161473,"visible":true,"origin":"","legend":"KEGG pathway enrichment analysis based on the differentially expressed genes at shoot formation stage\nThe longitudinal coordinates was KEGG pathway; the abscissa was richfactor; the size of the dots in the graph represented the number of differentially expressed genes that were annotated to the path, and the color represented the significant P value of the path. The 20 most significant paths were shown in the picture.","description":"","filename":"figure8.jpg","url":"https://assets-eu.researchsquare.com/files/rs-1316/v1/figure_8.jpg"},{"id":2614181,"identity":"6eb9def1-d68a-44e8-9810-ce8137f6787f","added_by":"399a65fa-c8d4-440b-992c-2804352e6ac3","created_at":"2020-09-25 20:54:05","extension":"jpg","order_by":9,"title":"Figure 9","display":"","copyAsset":false,"role":"figure","size":84653,"visible":true,"origin":"","legend":"Pattern diagram of embryogenesis and plant formation induced by cucumber ovary culture\na embryogenesis and plant formation process; b molecular events at various stages; ES embryo sac; Ge globular embryo; H-se heart-shaped embryo; T-se torpedo embryo; C-se cotyledon embryo; RP regenerated plant.","description":"","filename":"figure9.jpg","url":"https://assets-eu.researchsquare.com/files/rs-1316/v1/figure_9.jpg"},{"id":13467844,"identity":"d5cbbbe7-3ff6-4303-a1d2-7d94efa36d6e","added_by":"auto","created_at":"2021-09-16 20:56:34","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":1741095,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-1316/v1/502ca03b-e5b8-42d6-bc89-6ebf1520e1d2.pdf"},{"id":2614187,"identity":"7d67c40e-9c47-4213-af84-bb0d77ba3494","added_by":"399a65fa-c8d4-440b-992c-2804352e6ac3","created_at":"2020-09-25 20:54:05","extension":"docx","order_by":1,"title":"","display":"","copyAsset":false,"role":"supplement","size":36977,"visible":true,"origin":"","legend":"","description":"","filename":"supplement1.docx","url":"https://assets-eu.researchsquare.com/files/rs-1316/v1/supplement_1.docx"}],"financialInterests":"","formattedTitle":"Analysis of differentially expressed genes during embryogenesis in ovary culture of cucumber (Cucumis sativus L.)","fulltext":[{"header":"Background","content":" \n\u003cp\u003e\n Gametophyte culture involved a stress-induced reprogramming of male or female gametophytes\n to develop into embryo-like structures, which can be directly regenerated into completely\n homozygous, doubled haploid (DH). The study of male gametophyte (microspore) embryogenesis began in the 1960s [1]. Microspore embryogenesis in vitro can be easily monitored and provides a convenient\n experimental platform for large-scale physiological and biochemical analyses [2]. However, the study of gynogenesis is not simple; the embryo sac is small and embedded in surrounding\n tissues, making early embryos difficult to observe and isolate. Despite these difficulties, to inducing gynogenesis in vitro, especially when the\n haploid induced by the male gametophyte is not successful or the induction rate is\n too low and when self-incompatibility, male-sterile and dioecious plants occurs, is\n an effective way to obtain haploid or double haploid plants [3]. In vitro gynogenesis is mainly used in the \u003cem\u003eCucurbitaceae\u003c/em\u003e, and previous studies have revealed that this technique is still imperfect and that the embryogenesis rate is low [4-11]. The study of the mechanism of embryogenesis has therefore been delayed, and the understanding of the mechanisms underlying gynogenesis\n induction is extremely limited.\u003c/p\u003e\n \n\u003cp\u003e\n In general, embryogenesis is a stress-induced phenomenon. Micropores are cultured and develop into embryos in vitro, and stress\n treatments involving cold or heat shock and hormones are used as trigger factors to\n induce gametocytes to follow the sporophyte development pathway [12]. The genes associated with the reprogramming phase and the early stage of embryogenesis\n were slowly characterized, such as \u003cem\u003eAGL15\u003c/em\u003e-related proteins and the \u003cem\u003eAGL15\u003c/em\u003e gene, the \u003cem\u003eBBM\u003c/em\u003e gene and \u003cem\u003eHSP\u003c/em\u003egene [13-15]. More recently, the use of functional genomics tools has made it possible to identify\n more genes associated with different stages of microspore embryogenesis. In \u003cem\u003eBrassica napus\u003c/em\u003e, the expression of embryogenesis-related genes such as \u003cem\u003eBBM\u003c/em\u003e, \u003cem\u003eLEC\u003c/em\u003e\u003cem\u003e1\u003c/em\u003e and \u003cem\u003eLEC2\u003c/em\u003e was detected after culture under heat stress at 32 \u003csup\u003eo\u003c/sup\u003eC for three days; after incubation at 24 \u003csup\u003eo\u003c/sup\u003eC, embryo-specific expressed genes, such as \u003cem\u003eLEC1\u003c/em\u003e, \u003cem\u003eLEC2\u003c/em\u003ewere detected [2]. The expression of genes associated with metabolism, chromosome remodeling, transcription, and translation\n signaling was up-regulated during the stress treatment stage in tobacco [16]. Similarly, genes related to metabolism, cell walls and the cell membrane, cell tissue\n control, cell communication and signal transduction were detected in the early stage\n of rape embryogenesis [17]. High levels of \u003cem\u003eBBM\u003c/em\u003e and \u003cem\u003eLEC\u003c/em\u003e gene expression have been confirmed in the early embryonic development of sweet pepper\n anther cultures [18]. The transition from microspores to developing embryos is mainly manifested in the\n induction of transcription factor genes that play an important role in early embryogenesis.\n Many genes are involved in hormone biosynthesis and plant hormone signal transduction\n and are involved in secondary metabolism [19]. Studies have shown that differentially expressed genes can be detected by gene chips in in vitro gynogenesis of cucumber in the early\n stage and have suggested that phenylalanine metabolism and phenylalanine synthesis\n may play an important role in the early development of cucumber in vitro [20].\u003c/p\u003e\n \n\u003cp\u003e\n Cucumber (\u003cem\u003eCucumis sativus \u003c/em\u003eL.) is the fourth most important vegetable worldwide [21]. Cucumber production in China ranks first worldwide. To date, the cultivated area of\n cucumber in China has reached 1.258 Mha [22]. Cucumber is thought to be one of the oldest vegetable crops and has been domesticated\n in China for probably 2000 years [23]. To meet the needs of production, breeders are constantly looking for valuable germplasm\n resources of valuable characteristics, especially resistance to disease and environmental\n stresses such as cold, drought, or salt stress [24]. To facilitate these breeding efforts, we can collect extensive germplasm resources\n for variety improvements. However, cucumber is a cross-pollinated plant with obvious\n heterosis. Therefore, almost all the cucumber varieties in production are hybrids;\n however, the purification process of the hybrids is slow and long. In the breeding\n of parents, it takes 6-8 years of artificial self-breeding to develop a stable inbred\n line. To accelerate the purification of cucumber parents and improve breeding efficiency,\n haploid gametophyte cultures can be used to induce embryoids, and homozygous double\n haploids can be obtained in 1-2 years. DH technology is a powerful tool to speed up\n plant breeding. However, additional research is needed to improve our understanding\n of the genes or the roles of genes involved in embryogenesis and the mechanism of\n haploid induction in embryo sacs [25].\u003c/p\u003e\n \n\u003cp\u003e\n This study was based on the highly efficient in vitro ovary culture technology system\n of cucumber (unpublished). We divided the process from acquiring embryogenic potential\n to plant regeneration into three stages: early embryo development, embryo maturation\n (from pre-embryos to cotyledon embryos) and the shoot formation stage. We evaluated\n the embryo morphology and transcriptomes at different developmental stages during\n embryogenesis in cucumbers via ovary culture. The metabolism and biological process\n of embryogenesis in the gynogenesis of cucumber are discussed, and several key genes\n regulating embryogenesis were identified.\u003c/p\u003e\n "},{"header":"Result","content":"\n\u003ch2\u003eMorphological and cytological characterization of embryogenesis in cucumber\u003c/h2\u003e\n \n\u003cp\u003eThe embryogenesis early stage, embryo maturation stage and shoot formation stage were observed to discern the cytological changes during embryogenesis by ovary culture in cucumber. Ovaries at 6 h before anthesis were collected and pretreated at 4 \u003csup\u003eo\u003c/sup\u003eC. The ovaries were sterilized, sliced, inoculated into the medium, incubated for 2 d\n at 33 \u003csup\u003eo\u003c/sup\u003eC in the dark and then transferred to 25 \u003csup\u003eo\u003c/sup\u003eC until plants formed. Interestingly, based on the morphological and cytological observations, critical developmental transitions in embryos were discovered at 2 (T1), 10 (T2), 20 (T3), 30 (T4) and 60 (T5) d of culture. Fresh unpollinated ovaries before\n culture were named T0 (Fig. 1). At T0, fresh unpollinated ovaries were selected, in which the ovules were obvious (Fig. 1a) and the embryo sacs were mature (Fig. 1a1).At T1, heat shock stress induced the embryogenesis of cucumber (Fig. 1b). The most obvious sign that referred to the development of the embryo was the enlargement\n of the cells in the embryo sac (Fig. 1b1). The enlargement of microspores have been correlated with embryogenic potential acquisition during induction of androgenesis in many crop species [26~27]. After culturing for 2 d, 4 d and 6 d, the expanded cells formed 2 (Fig. 1b2), 4(Fig.\n 1b3) and multiple cell structures(Fig. 1b4), respectively. At T2, the explants were cultured for 10 d, the ovules that successfully protruded\n from the surrounding tissues, and a pro-embryo had formed (Fig. 1c). At T3, after\n 20 d of culture, the total number of embryos increased continuously in one ovary slice\n (Fig. 1d). From T3 to T4 (culture for 30 d), the whole embryo had formed, including\n globular embryos, heart-shaped embryos, torpedo-shaped embryos and cotyledon-shaped\n embryos (Fig. 1 g). At the same time, the color of the embryonic cells appeared green,\n indicating that the development of plants from the embryo would begin (Fig. 1e). At T5 (60 d after culture), shoots had formed (Fig. 1f).\u003c/p\u003e\n \n\u003cp\u003eAt T1, ovule enlargement and cell enlargement could be clearly observed, suggesting\n that the transformation of the development pathway from the gametophyte to the sporophyte\n was induced by stress and that the stage between T0 and T1 was the key to inducing\n embryogenesis. After the switching of cell development, embryogenic potential was\n acquired. Mitosis then occurred continuously to form pro-embryos, globular embryos,\n heart-shaped embryos, torpedo-shaped embryos and cotyledon-shaped embryos, similar\n to zygote development, and this process continued until T4. After 6 d of culture, the mature embryos developed into shoots. Therefore, we divided the six time points of ovary culture into three stages: T0 to T1 (embryogenesis early stage), T1 to T4 (embryo maturation stage), and T4 to T5 (shoot formation stage). Transcriptome data were generated via Illumina II HiSeq\u003csup\u003eTM\u003c/sup\u003e 2000 sequencing, and the statistics of RNA-seq alignment were shown in Table S1.\u003c/p\u003e\n \n\u003ch2\u003eValidation of differentially expressed genes\u003c/h2\u003e\n \n\u003cp\u003eValidation of the Illumina sequencing data and the expression patterns of the DEGs\n revealed by RNA-seq was performed to examine the expression patterns of 10 DEGs, including 7 genes involved\n in plant hormone signal transduction and 3 plant-pathogen interaction genes. The results\n showed that the relative expression levels revealed by RNA-seq and qRT-PCR were closely correlated (Pearson's r = 0.53, 0.77, 0.57, 0.74, 0.79, 0.69, 0.81,\n 0.81, 0.86, 0.62, 0.81). The fold changes in the qRT-PCR analysis were different from those in the RNA-seq analysis, which might be due to the difference in sensitivity between the qRT-PCR analysis and the RNA-seq technique. However, the qRT-PCR analysis showed that the\n up- and down-regulation trends of the differential gene expression were consistent\n with those of the RNA-seq analysis (Fig. 2).\u003c/p\u003e\n \n\u003ch2\u003eExpression of embryogenesis-related genes in cucumber ovary culture\u003c/h2\u003e\n \n\u003cp\u003eAt the early stage of embryogenesis, after heat shock stress, ovule enlargement was obvious (Fig. 1b). Cytological observations also showed\n that one of the cells in the embryo sac expanded (Fig. 1b1), suggesting acquisition of embryogenic potential.\u003c/p\u003e\n \n\u003cp\u003eThe differentially expressed genes (DEGs) at this stage (T0 vs T1) were identified as comprising 3468 up-regulated genes and 3065 down-regulated genes (Fig. 3). We performed a GO enrichment analysis of these genes, and the results showed that\n the proteins encoded by these genes were assigned to 4 biological processes, 3 molecular\n functions and 9 cellular components (Table S2). Different genes were involved in different biological processes, and the same gene was involved in a variety of biological processes. The majority of DEGs were involved in the ‘single-organism process’ (GO:0044699) and the ‘oxidation-reduction process’ (GO:0055114). They were mainly distributed among\n the terms ‘membrane’ (GO:0016020), ‘membrane part’ (GO:0044425), ‘intrinsic component\n of membrane’ (GO:0031224) and ‘integral component of membrane’ (GO:0016021), having good molecular function in oxidoreductase activity (GO:0016491).\n Among the DEGs, the expression of related genes involved in the functional classification\n might play an important role in early embryogenesis.\u003c/p\u003e\n \n\u003cp\u003eTo compare and summarize the results of this stage, we performed a pathway analysis\n to identify potential target genes. Based on the KEGG database, the pathway enrichment analysis of these genes revealed\n significant involvement in 8 distinct pathways (Fig. 4): pathways involving plant hormone signal transduction, phenylpropanoid biosynthesis, plant-pathogen interaction,\n glutathione metabolism, cysteine and methionine metabolism, drug metabolism-cytochrome P450, phenylalanine metabolism and metabolism of xenobiotics by cytochrome P450.\u003c/p\u003e\n \n\u003ch2\u003e\u003cem\u003ePlant hormone signal transduction genes\u003c/em\u003e\u003c/h2\u003e\n \n\u003cp\u003eThe DEGs were significantly enriched in the plant hormone signal transduction pathway\n (Ko04075), and a large number of genes encoding cytokinin transfer proteins (\u003cem\u003eCsa2M373410.1\u003c/em\u003e, \u003cem\u003eCsa6M067360.1\u003c/em\u003e, \u003cem\u003eCsa7M452370.1\u003c/em\u003e), cytokinin response proteins (\u003cem\u003eCsa1M006300.1\u003c/em\u003e, \u003cem\u003eCsa3M822100.1\u003c/em\u003e, \u003cem\u003eCsa4M436980.1\u003c/em\u003e, \u003cem\u003eCsa5M223020.1\u003c/em\u003e, \u003cem\u003eCsa5M434550.1\u003c/em\u003e, \u003cem\u003eCsa5M603910.1\u003c/em\u003e, \u003cem\u003eCsa5M623800.1\u003c/em\u003e, \u003cem\u003eCsa6M383530.1\u003c/em\u003e), ethylene-binding proteins (\u003cem\u003eCsaUNG009930\u003c/em\u003e) and ethylene-responsive protein transcription factors (\u003cem\u003eCsa3M389850.1\u003c/em\u003e) were detected, which were significantly up-regulated after heat treatment. Moreover, the partially\n up-regulated DEGs encoding auxin influx carriers, auxin-responsive proteins, SAUR\n family proteins, abscisic acid receptor proteins, serine/threonine-protein kinases,\n ABA-responsive element binding factors, and brassinosteroid signal-responsive protein\n kinases were significantly enriched in 8 pathways, some of which involved tryptophan\n metabolism as well as zeatin, carotene, and brassinosteroid biosynthesis, suggesting\n that hormone responses and hormone signal transduction were closely related to stress-induced embryogenesis. The up-regulated genes encoding major enzymes and receptor proteins involved in the\n plant hormone signal transduction pathway were expressed during the whole stage of\n ovary culture (Fig. 5). Most of the genes were expressed dynamically after their up-regulated expression\n during embryogenesis. Only two genes (\u003cem\u003eCsa6M147590.1\u003c/em\u003e,\u003cem\u003e Csa3M389850.1\u003c/em\u003e) were up-regulated at the stage of embryo initiation; both almost no expression in the later stages. \u003cem\u003eCsa6M147590.1 \u003c/em\u003eand \u003cem\u003eCsa3M389850.1 \u003c/em\u003eencode auxin-induced proteins and ethylene-responsive transcription factors, respectively,\n indicating that they might be closely related to embryogenesis.\u003c/p\u003e\n \n\u003ch2\u003e\u003cem\u003eStress response proteins related to plant defens\u003c/em\u003e\u003cem\u003ee m\u003c/em\u003e\u003cem\u003eechanisms\u003c/em\u003e\u003c/h2\u003e\n \n\u003cp\u003e\n DEGs were significantly enriched in the plant-pathogen interaction pathway, which\n was related to plant defense mechanisms and complex physiological responses to heat\n shock-induced embryogenesis. In this regulatory pathway, the expression of \u003cem\u003ePti1\u003c/em\u003e (\u003cem\u003eCsa7M420160.1\u003c/em\u003e), \u003cem\u003eHSP90 \u003c/em\u003e(\u003cem\u003eCsa3M183950.1\u003c/em\u003e), \u003cem\u003eEDS1\u003c/em\u003e (\u003cem\u003eCsa1M006320.1\u003c/em\u003e) and other related protein genes were up-regulated. \u003cem\u003eHeat shock protein 90\u003c/em\u003e (\u003cem\u003eHSP90\u003c/em\u003e) was a protein that was up-regulated under heat stress. This gene had been highly\n conserved during evolution and had plant defense functions.\u003c/p\u003e\n \n\u003ch2\u003eResponse of genes involved in microtubule organization in cucumber ovary culture\u003c/h2\u003e\n \n\u003cp\u003e\n In the stage of embryo maturation, the gametophyte development pathway switched to the sporophyte development\n pathway, in which the cells enlarged and initiating mitosis, forming 2-cell, 4-cell and multicellular pro-embryos after heat stress (Fig. 1b1~b4); afterward, globular embryos, heart-shaped embryos, torpedo-shaped embryos and cotyledon\n embryos were formed (Fig. 1g). The DEGs analyses of T1 vs T2, T2 vs T3 and T3 vs T4\n revealed 480 continuously expressed genes. Additionally, the numbers of uniquely expressed genes between T1 and T2, T2 and T3\n and T3 and T4 were 2236, 371 and 1558, respectively (Fig. 6).\u003c/p\u003e\n \n\u003cp\u003eGO categories of the 480 continuously expressed genes were significantly assigned to 2 biological processes,\n including cell movement or subcellular component (GO:0006928) and microtubule-based movement (GO:0007018), which were mainly associated\n with microtubules (GO:0005874), supramolecular fibers (GO:0099512) and polymeric cytoskeletal fibers (GO:0099513). Subsequently,\n proteins encoded by these genes were classified into 11 functional categories, including protein complex binding (GO:0032403), microtubule binding (GO:0008017),\n tetrapyrrole binding (GO:0046906), tubulin binding (GO:0015631), microtubule motor activity (GO:0003777), motor activity (GO:0003774), macromolecular complex binding (GO:0044877), heme binding (GO:0020037),\n oxidoreductase activity (GO:0016491), iron ion binding (GO:0005506) and oxidoreductase activity(GO:0016491),\n acting on paired donors or the molecular incorporation of oxygen (GO:0016705) (Table\n S3). This findings indicated that, during certain stages, the cytoskeleton structure was continuously built and maintained by the action of microtubule\n binding proteins and enzyme modifications, providing the basis for positioning the various organelles and implementation functions,\n which ensured the orderly activities in various cells in time and space.\u003c/p\u003e\n \n\u003cp\u003eTo determine the involvement of these differentially expressed genes in embryo development,\n we performed a pathway analysis to identify the potential target genes, and 17 significant\n pathways were obtained by mapping to the KEGG database (Fig. 7): protein processing in the endoplasmic reticulum, phenylpropanoid biosynthesis, photosynthesis-antenna proteins, limonene and pinene degradation, meiosis-yeast, the estrogen signaling pathway, cell cycle, cell cycle-yeast, diterpenoid biosynthesis, chloroalkane and chloroalkene degradation, carotenoid biosynthesis, progesterone-mediated oocyte maturation, glycerolipid metabolism, histidine metabolism, fatty acid degradation, antigen processing and presentation, and ascorbate and aldarate metabolism. Among them, the pathways of protein processing in the endoplasmic reticulum and\n phenylpropanoid biosynthesis were the most significant.\u003c/p\u003e\n \n\u003ch2\u003eExpression of main oxidation-reduction and metabolic process-related genes in cucumber\n ovary culture\u003c/h2\u003e\n \n\u003cp\u003eIn the stage of shoot formation, the embryos further differentiated into shoots (Fig. 1f). In total, 3320 genes were differentially expressed between T4 and T5, among which 1383 genes were up-regulated and 1837 were down-regulated in T5 (Fig. 3). Next, we performed an enrichment analysis\n of the genes using Gene Ontology (GO), which revealed 18 biological processes, 15\n molecular functions and 5 cellular components (Table S4). The terms annotated under the biological process category mainly included oxidation-reduction processes (GO: 0055114), the regulation of primary metabolic processes (GO: 0080090),\n the regulation of cellular metabolic processes (GO: 0031323), and the regulation of\n macromolecule metabolic processes (GO: 0060255). The molecular function category mainly\n included the terms oxidoreductase activity (GO: 0016491) and DNA binding (GO: 0003677).\n The cellular component category mainly included extracellular region (GO: 0005576).\n The development of the ovary culture was mainly focused on the process of oxidation-reduction\n and metabolism.\u003c/p\u003e\n \n\u003cp\u003eThe DEGs were subsequently annotated using the KEGG database to identify pathway enrichments. A variety of pathways were found to be\n significantly enriched (Fig. 8), including those involved in phenylpropanoid biosynthesis; plant hormone signal transduction; stilbenoid, diarylheptanoid and gingerol biosynthesis; metabolism of xenobiotics by cytochrome P450; drug metabolism-cytochrome P450; flavonoid biosynthesis; phenylalanine metabolism; starch and sucrose metabolism; and zeatin biosynthesis. Most of the differentially expressed genes were mainly enriched in phenylpropanoid\n biosynthesis, plant hormone signal transduction, phenylalanine metabolism, and starch\n and sucrose metabolism. Pathways involving phenylpropanoid biosynthesis and phenylalanine were important\n for metabolizing secondary metabolites in plants and were closely related to cell\n differentiation and pigmentation in plant development. In the plant hormone signal transduction pathway, the expression of an auxin response\n factor, the gibberellin receptor \u003cem\u003eGID1\u003c/em\u003e, ethylene-responsive transcription factor 1 (\u003cem\u003eERF1\u003c/em\u003e) and cyclin D3 were up-regulated, promoting plant development and maturation, cell\n enlargement and division as well as promoting the development of stems. In addition,\n the DEGs involved in starch and sugar metabolism were also significantly enriched\n and mainly encoding 3-βD glucosidase, T6Ps, GlgB, endoglucanase and alpha-trehalase, which satisfied the requirements of embryo growth and development.\u003c/p\u003e"},{"header":"Discussion","content":" \n\u003cp\u003eMany previous experiments have shown that early embryogenesis of microspore cultures\n can be divided into three main phases: acquisition of embryogenic potential, initiation of cell division, and pattern formation [27]. In the present study, the whole process was investigated from acquisition of embryogenic\n potential to plant regeneration. We divided the ovary culture into three stages. In the early stage of embryogenesis, the acquisition of embryogenic potential by stress (e.g., low or high temperature\n and hormone induction) was observed together with the repression of gametophytic development,\n leading to the dedifferentiation of cells. In the embryo maturation stage, cell divisions gave rise to the formation of pro-embryos\n (cell clumps), globular embryos, heart-shaped embryos, torpedo-shaped embryos and\n cotyledon embryos (mature embryos). In the stage of shoot formation, mature embryos developed into shoots. There are\n active molecular events that regulate embryo development at different stages of development.\u003c/p\u003e\n \n\u003ch2\u003eEmbryogenesis-related genes are expressed in cucumber ovary culture\u003c/h2\u003e\n \n\u003cp\u003e\n The embryogenic ability and transformation of somatic cells were regulated by various plant hormones, such as auxin, abscisic acid, cytokinin and ethylene [28]. Moreover, the synthesis of auxin and ethylene increased significantly after ovary\n pollination and fertilization [29]. In normal growth, mitotic asymmetry of fertilized eggs was found in the gametophyte\n development pathway. Previous work using molecular and genetic approaches has identified\n auxin as a critical signal for the proper development of the asymmetric structure\n of the female gametophyte in Arabidopsis [30]. In \u003cem\u003ein vitro\u003c/em\u003e culture, a certain stress condition is needed to block the development direction\n of the original gametophyte, followed by turning into the direction of the sporophyte\n and carrying on the symmetry splitting to eventually lead to embryogenesis [31]. In ovary culture of cucumber, heat shock pretreatment, silver nitrate, genotype\n and hormone combination factors could play key roles in embryo and callus production\n independently and simultaneously [32]. In the ovary culture of linseed, cultivar, combination of growth regulators, type\n of carbohydrates and their interaction significantly influenced callus induction and\n shoot formation frequency, and a relatively high shoot regeneration frequency was\n obtained when the medium was supplemented with TDZ and NAA [33]. These results implied that the expression of specific genes in embryogenesis was\n activated by exogenous hormone regulation, which laid the foundation for DNA replication\n in the cell division stage.\u003c/p\u003e\n \n\u003cp\u003eIn our study, a large number of DEGs were enriched in the plant hormone signal transduction\n pathway during the stress process of embryogenesis. The pathway annotation obtained by KEGG analysis showed that phytohormones such as cytokinin and ethylene were significantly up-regulated in multiple biosynthesis\n and metabolic pathways, while up-regulated and down-regulated genes related to auxin\n and abscisic acid both existed. Generally, high-level endogenous cytokinin and low-level\n endogenous ethylene are beneficial for the acquisition of embryogenic ability. Embryogenesis\n is the comprehensive performance of interactions among different hormones, and all\n kinds of endogenous hormones show dynamic changes. The up-regulated expression of\n the main enzymes and receptor proteins in the plant hormone signal transduction pathway\n could promote embryogenesis, indicating that the high expression of related genes\n may play important roles in the process of switching from the gametophyte to the sporophyte development pathway. In this pathway, we found two genes related to hormone regulation that were highly\n expressed only at the stage of early embryo development (\u003cem\u003eCsa6M147590.1\u003c/em\u003e and \u003cem\u003eCsa3M389850.1\u003c/em\u003e);moreover, the functional categories of their specific functions and the location of\n the regulatory pathways were identified, thus laying an important theoretical basis for the elucidation of the mechanism.\u003c/p\u003e\n \n\u003cp\u003eHeat shock protein 90 (\u003cem\u003eHSP90\u003c/em\u003e) is an evolutionarily conserved molecular chaperone induced by abiotic stress. \u003cem\u003eHSP90\u003c/em\u003e plays an important role in cell cycle control, genomic silencing, protein transduction\n and signal transduction[20]. The \u003cem\u003eHSP90\u003c/em\u003e chaperone is involved in maintaining phenotypic plasticity and developmental stability [34-35]. Previous studies have shown that \u003cem\u003eHSP\u003c/em\u003e genes were specifically expressed in the spore stress-induction process of many crop\n species [36-38]. \u003cem\u003eHSPs\u003c/em\u003e had been described as genes associated with the reprogramming phase and the early\n stages of embryogenesis [39]. After heat shock stress treatment, we found that 26 up-regulated genes encoded proteins and transcription\n factors related to heat shock, and \u003cem\u003eHSP90 \u003c/em\u003e(\u003cem\u003eCsa3M183950.1\u003c/em\u003e) was involved in the regulation of the hypersensitive response in the plant and pathogen\n interaction pathway. The expression of these genes might play a regulatory role in\n the process of embryogenesis.\u003c/p\u003e\n \n\u003cp\u003e\n The \u003cem\u003eBBM \u003c/em\u003egene was the first key gene to be isolated in the process of spore cell division;\n which was first expressed in zygotic embryogenesis and microspore embryogenesis [14, 40]. We found that two \u003cem\u003eBBM\u003c/em\u003e genes (\u003cem\u003eCsa3M827310.1\u003c/em\u003e, \u003cem\u003eCsa3M827320.1\u003c/em\u003e) were up-regulated at the beginning of embryogenesis and were subsequently down-regulated.\n In addition, \u003cem\u003eAGL15\u003c/em\u003e a member of the MADS family, was also considered a regulatory protein that acts at\n the start of cell division [13]. We found that the \u003cem\u003eAGL\u003c/em\u003e gene (\u003cem\u003eCsa3M258140.1\u003c/em\u003e) was up-regulated in our study, indicating that all three of these genes play important\n roles in embryogenesis. We considered that \u003cem\u003eBBM, HSP90 \u003c/em\u003eand \u003cem\u003eAGL \u003c/em\u003emight be the critical genes involved in the induction of embryogenesis by ovary culture\n in cucumber.\u003c/p\u003e\n \n\u003ch2\u003eMicrotubule organization genes play an important role in the embryogenesis of cucumber\n ovary culture\u003c/h2\u003e\n \n\u003cp\u003eThere is a similar process between gametophyte embryogenesis and zygotic embryogenesis,\n and the only difference is that the former occurs in sex cells, while the latter occurs in fertilized egg cells. The process of gametophyte embryogenesis development is similar to that of zygote differentiation\n into embryos. Microtubules play a transport role in cells, participating in the construction of cell walls and\n promoting cell differentiation and division. Electron microscope-based observations of carrot culture showed that the appearance of microtubules was accompanied with the formation of somatic embryos [41]. The different patterns of microtubule organization in the cells of the mature embryo\n sac reflect their structural adaptations for their future function [42].\u003c/p\u003e\n \n\u003cp\u003eIn our study, after embryogenic potential was acquired, the different embryo shapes\n formed through active cell metabolism and rapid cell division. GO analysis showed that in different embryo periods continuously expressed genes participated in microtubule-based movement and cell or subcellular component processes, providing a good function of protein complex binding, microtubule binding, tetrapyrrole binding, tubulin binding\n and other microtubule activities; these genes mainly participated in protein processing of the endoplasmic reticulum\n and phenylpropanoid biosynthesis. We considered that the cellular reprogramming and morphological changes in embryos\n were controlled by microtubule organization genes.\u003c/p\u003e\n \n\u003ch2\u003eExpression of main oxidation-reduction and metabolic process-related genes in cucumber\n ovary culture\u003c/h2\u003e\n \n\u003cp\u003eMature embryos are further cultured to form shoots, and at this stage, many biological\n processes are involved. GO analysis showed that differentially expressed genes participated\n in the oxidation-reduction process, a variety of metabolic processes, biosynthesis\n processes and the response to auxin, together with the molecular functions of oxidoreductase\n activity, binding capacity, catalytic activity, transporter activity, transcription\n factor activity and so on. Among them, the regulation of transcription, DNA template (GO: 0006355), the regulation of RNA biosynthetic processes (GO: 2001141), the regulation\n of RNA metabolic processes (GO: 0051252), DNA binding (GO: 0003677), transcription factor activity (GO: 0003700), and sequence-specific DNA binding (GO: 0043565) lay a foundation for\n cell division and differentiation. In addition, two processes, cellular response to\n auxin stimulus (GO: 0071365) and response to auxin (GO: 0009733), might promote the\n embryo to form shoots.\u003c/p\u003e\n \n\u003cp\u003eAt the stage of shoot formation, a large number of transcription factors were expressed.\n The first constitutes transcription factors that contain a B3 domain [43]. These transcription factors can encode regulatory proteins involved in the embryonic\n development process, maintaining embryo development during late embryonic development [43, 44]. In our study, total of 6 transcription factors(\u003cem\u003eCsa2M359980.1\u003c/em\u003e, \u003cem\u003eCsa2M292240.1\u003c/em\u003e, \u003cem\u003eCsa5M608380.1\u003c/em\u003e, \u003cem\u003eCsa6M489980.1\u003c/em\u003e, \u003cem\u003eCsa6M489940.1 \u003c/em\u003eand \u003cem\u003eCsaUNG031640\u003c/em\u003e) that contain a B3 domain were up-regulated\u003cem\u003e.\u003c/em\u003eThe second was the \u003cem\u003eAIL\u003c/em\u003e gene from the \u003cem\u003eAP2/ERF\u003c/em\u003e family, which is involved in key developmental processes throughout the whole plant life\n cycle. Some genes in the \u003cem\u003eAP2/ERF\u003c/em\u003e gene family are expressed in many tissues and participate in many plant development\n processes, such as embryogenesis and shoot development [14, 45). In addition, we found that the expression levels of 9 genes (\u003cem\u003eCsa1M423190.1\u003c/em\u003e,\u003cem\u003e Csa3M114480.1\u003c/em\u003e, \u003cem\u003eCsa3M652380.1\u003c/em\u003e,\u003cem\u003e Csa3M114470.1\u003c/em\u003e,\u003cem\u003e Csa4M644740.1\u003c/em\u003e, \u003cem\u003eCsa5M175970.1\u003c/em\u003e, \u003cem\u003eCsa5M608380.1\u003c/em\u003e,\u003cem\u003e Csa6M496390.1\u003c/em\u003e,and \u003cem\u003eCsaUNG031640\u003c/em\u003e) in the\u003cem\u003eAP2/ERF\u003c/em\u003e family were up-regulated. Furthermore, \u003cem\u003eWUS \u003c/em\u003ehomologous domain proteins not only alter the cell fate of the shoot and flower meristem but also promote the\n development of somatic embryos into seedlings. The function of \u003cem\u003eWUS\u003c/em\u003e proteins has no direct connection with the characteristics of the embryo but will\n alter the development state of the tissue by maintaining cells in an undifferentiated\n state in response to different stimuli [46]. The ability to transform the vegetative growth phase to the embryonic stage by \u003cem\u003eWUS\u003c/em\u003e, as well as to eventually form somatic cells, indicates that this homologous domain protein also plays an important role in embryo maturation in addition to its role in the development of the meristem [47]. During this stage, two genes encoding \u003cem\u003eWUS\u003c/em\u003e homologous domain proteins (\u003cem\u003eCsa6M301060.1\u003c/em\u003e) and (\u003cem\u003eCsa6M505860.1\u003c/em\u003e) were up-regulated. We considered that those homologous domain proteins might directly or indirectly regulate\n shoot formation.\u003c/p\u003e\n \n\u003ch1\u003eConclusion\u003c/h1\u003e\n \n\u003cp\u003eSeveral studies have reported on the embryogenesis of ovary culture in cucumbers. However, the mechanism that drives the process from embryogenic acquisition to the\n formation of embryos and plant regeneration is not well understood. In this study,\n we explored embryogenesis mechanism of ovary culture in cucumber. Inducing cucumber embryogenesis could be divided into three stages: early embryogenesis, embryo maturation (from proembryos to cotyledon embryos) and shoot formation (Fig. 9a).\u003c/p\u003e\n \n\u003cp\u003eThe early stage of embryogenesis was the turning point for the formation of embryos, which had experienced dedifferentiation and loss of photosynthetic capacity, requiring the provision of exogenous nutrients and carbon sources such as plant hormones and\n sucrose in the media. Some physical factors, such as temperature, light quality, photoperiod and the presence\n of specific hormones, could affect the ability of the embryo sac to adapt to these\n conditions and survive the developmental transition. These dynamic changes were helpful\n for cell physiology reprogramming, metabolic alterations, revival of cell division\n and differentiation, morphogenesis and so on. Therefore, the cells of the embryo sac\n began to divide by dedifferentiation and started to form cell clumps (Fig. 1b1~b4). Although there was no obvious change in appearance, the metabolism of some macromolecules\n in the cells changed markedly. A large number of hormone-related genes, cell protection-related\n genes, and some unique protein kinases and transcription factors (\u003cem\u003eCsa6M147590.1\u003c/em\u003e and\u003cem\u003e Csa3M389850.1\u003c/em\u003e) were expressed only at the stage of early embryo development\u003cem\u003e. \u003c/em\u003eThe reported embryogenesis-related genes\u003cem\u003e BBM\u003c/em\u003e,\u003cem\u003e HSP90\u003c/em\u003e and \u003cem\u003eAGL \u003c/em\u003ewere also actively expressed during this stage (Fig. 9b).\u003c/p\u003e\n \n\u003cp\u003eIn the stage of embryo maturation, GO analysis showed continuously expressed genes participated in microtubule-based movement, movement of the cell or subcellular component processes,\n giving a good function of protein complex binding, microtubule binding, tetrapyrrole\n binding, tubulin binding and other microtubule activities, which are involved in protein\n processing in the endoplasmic reticulum and phenylpropanoid biosynthesis(Fig. 9b).\u003c/p\u003e\n \n\u003cp\u003eShoot formation was regulated by 6 transcription factors that contain a B3 domain, total of 9 genes\n in the \u003cem\u003eAP2/ERF\u003c/em\u003e family and 2 genes encoding \u003cem\u003eWUS\u003c/em\u003e homologous domain proteins (Fig. 9b). \u003c/p\u003e"},{"header":"Methods","content":"\u003ch2\u003ePlant materials\u003c/h2\u003e\n \n\u003cp\u003e\n ‘SG033’ (F\u003csub\u003e1\u003c/sub\u003e cultivar of Kunming Huaxing seed Industry Co., Ltd., Kunming, Yunnan) was used for\n ovary culture; this material had a high embryo rate after a large number of screenings.\n Cucumber plants from southern China were cultivated and stored in our lab. The plants\n were grown in a greenhouse with a 12 h photoperiod, a mean daily air temperature of 25/15 \u003csup\u003eo\u003c/sup\u003eC (day/night), a relative humidity of 85% and a photosynthetic photo flux density of\n 500 μmol·m\u003csup\u003e-2\u003c/sup\u003e·s\u003csup\u003e -1\u003c/sup\u003e at the Horticulture Institute of Guizhou Academy of Agricultural Sciences, Guizhou,\n China. Ovaries were collected from ‘SG033’ for one month after the first anthesis\n flower appeared.\u003c/p\u003e\n \n\u003ch2\u003eOvary culture\u003c/h2\u003e\n \n\u003cp\u003e\n Ovaries were harvested at 6 h before anthesis, after treated for 24 h in a refrigerator\n at 4 \u003csup\u003eo\u003c/sup\u003eC, which papilloma was removed. The ovaries were then sterilized in 75% ethanol for 30-60 s and rinsed\n in sterile distilled water 3 times, followed by soaking in a 0.5% sodium hypochlorite\n for 20 min and then rinsing in sterile distilled water 3 times. The ovaries were cut into small round slices of 2 mm under sterile conditions and then placed on 30 mL of solid media in cylindrical flasks.\n In all experiments, MS medium was used as a basal medium; it was supplemented with\n 0.06 mg·L\u003csup\u003e-1\u003c/sup\u003e TDZ and 3% (w/v) sucrose and was solidified with 7% (w/v) agar. The pH of the medium\n was adjusted to 5.9 before autoclaving at 116 \u003csup\u003eo\u003c/sup\u003eC and 1.1 kg/cm\u003csup\u003e2\u003c/sup\u003e for 30 min. All materials were treated with high temperature of 33 \u003csup\u003eo\u003c/sup\u003eC for 2 d, after which recovery growth was allowed until shoot formation occurred at\n 25 \u003csup\u003eo\u003c/sup\u003eC under a 16/8 h (light/dark) photoperiod with a 4000 lx light intensity.\u003c/p\u003e\n \n\u003ch2\u003eRNA isolation\u003c/h2\u003e\n \n\u003cp\u003e\n The materials were harvested separately at the following time points: T0 (the ovules\n were cultured for 0 d, fresh unpollinated ovaries), T1 (the ovules were cultured for\n 2 d), T2 (the embryos were cultured for 10 d), T3 (the embryos were cultured for 20\n d), T4 (the embryos were cultured for 30 d), and T5 (the shoots after culture for\n 6 d). These materials were stored at -80 \u003csup\u003eo\u003c/sup\u003eC for RNA extraction, and the remaining materials were maintained in culture to observe\n plantlet regeneration. Furthermore, all experimental procedures, such as culture medium replacement and sample collection, were\n performed under similar conditions to minimize possible circadian effects. Ovules at T0 and T1 were extracted by hand under a stereoscopic microscope, and embryo-like\n structures with integuments were selected directly at the time points from T2 to T5.\u003c/p\u003e\n \n\u003cp\u003e\n As described above, samples were collected from six time points of ovary culture for\n RNA-seq analysis. For each sample, the ovules or embryo-like structures were immersed\n in liquid nitrogen in a mortar and ground into powder (three biological replicates per sample). Total RNA was isolated using TRIzol (Invitrogen),\n and DNase I (Fermentas) digestion was performed for 30 min at 25 \u003csup\u003eo\u003c/sup\u003eC to remove DNA, according to the manufacturer’s instructions. The integrity and quality\n of the total RNA were checked using a NanoDrop 1000 spectrophotometer and formaldehyde-agarose\n gel electrophoresis. RNA was used only when the OD\u003csub\u003e260\u003c/sub\u003e: OD\u003csub\u003e280\u003c/sub\u003e ratio was above 1.8.\u003c/p\u003e\n \n\u003ch2\u003eTranscriptome analysis\u003c/h2\u003e\n \n\u003cp\u003eTranscriptomes were analyzed by an Illumina HiSeq\u003csup\u003eTM\u003c/sup\u003e 2000. The reads with more than 10% unknown bases and those of low quality were removed\n from the raw sequencing data using the software Illumina GA Pipeline (v1.6). The filtered\n reads were aligned to the cucumber Chinese long v2 genome using TopHat2 (\u003ca href=\"http://tophat.cbcb.umd.edu/)\"\u003e\u003ca href=\"http://tophat.cbcb.umd.edu/\"\u003ehttp://tophat.cbcb.umd.edu/\u003c/a\u003e)\u003c/a\u003e software. The reference genome and gene database information were obtained from the\n public website: \u003ca href=\"http://cmb.bnu.\"\u003ehttp://cmb.bnu.\u003c/a\u003e\u003ca href=\"http://edu.cn/Cucumis_sativus_v20/\"\u003eedu.cn/Cucumis_sativus_v20/\u003c/a\u003e [42]. Differential expression analysis was determined using DESeq (version 1.18.0). The resulting \u003cem\u003eP\u003c/em\u003e values were adjusted by using Benjamini and Hochberg’s approach for controlling the\n false discovery rate. Genes with an adjusted \u003cem\u003eP \u003c/em\u003evalue of \u0026lt;0.05 and an absolute value of |log2fold change≥1| according to DESeq were considered differentially expressed. GO annotation was\n performed using Blast2GO software (GO association was performed by a BLASTX search\n against the NCBI NR database). GO enrichment analysis of differentially expressed\n genes (DEGs) was then performed via the BiNGO plugin for Cytoscape. Overrepresented\n GO terms were identified using a hypergeometric test with a significance threshold\n of 0.05 after the Benjamini-Hochberg FDR correction. KEGG enrichment analysis of the differentially expressed genes was performed\n using the KOBAS (2.0) [48] software. A ‘\u003cem\u003eP\u003c/em\u003e value ≤ 0.05’ was used as the threshold to judge the significantly enriched pathways of the\n differentially expressed genes.\u003c/p\u003e\n \n\u003ch2\u003eReal-time quantitative PCR\u003c/h2\u003e\n \n\u003cp\u003eTen DEGs were selected from all the DEGs for analysis using real-time quantitative\n PCR (qRT-PCR). For this step, the total RNA was treated with the DNase I enzyme and subsequently\n converted to single-strand cDNA by applying the GoScript\u003csup\u003eTM\u003c/sup\u003e Reverse Transcription System (Promega) according to the manufacturer’s protocol. The specific gene primers were designed\n using Primer 5.0 according to the cDNA sequences in Table S5. Actin was used as the\n internal control to normalize small differences in template quantities. Quantitative\n real-time PCR was performed with GoTaq qPCR Master Mix (Promega) in a Bio-Rad iQ1\n real-time PCR system (Bio-Rad), and the collected data were analyzed using the 2\u003csup\u003e-ΔΔCt\u003c/sup\u003e method. \u003c/p\u003e"},{"header":"Abbreviations","content":" \n\u003cp\u003eDH, Double haploid; DEG, Differentially expressed genes; qRT-PCR, Quantitative real-time\n PCR; SSR, Simple sequence repeats; FDR, False discovery rate; \u003cem\u003eAGL15\u003c/em\u003e, Agamous-like15; \u003cem\u003eAP2/ERF, \u003c/em\u003eEthylene-responsive element-binding factor; \u003cem\u003eBBM\u003c/em\u003e, Baby boom; \u003cem\u003eHSP\u003c/em\u003e, Heat shock proteins; \u003cem\u003ePti\u003c/em\u003e\u003cem\u003e1\u003c/em\u003e, Pto-interacting protein 1; \u003cem\u003eAIL\u003c/em\u003e,Aintegumenta-like; \u003cem\u003eLEC,\u003c/em\u003e Leafy Cotyledon; \u003cem\u003eGID1\u003c/em\u003e, gibberellin receptor, \u003cem\u003eERF1\u003c/em\u003e, ethylene-responsive transcription factor1;\u003cem\u003eWUS\u003c/em\u003e, WUSCHEL; TDZ, thidiazuron; \u003c/p\u003e\n \n\u003cp\u003e3-βD glucosidase, 3-beta-D-glucosidase; T6Ps, Trehalose 6-phosphate synthase; GlgB, 1,4-alpha-glucan branching enzyme.\u003c/p\u003e"},{"header":"Declarations","content":" \n\u003ch2\u003eEthics approval and consent to participate\u003c/h2\u003e\n \n\u003cp\u003eNot applicable.\u003c/p\u003e\n \n\u003ch2\u003eConsent for publication\u003c/h2\u003e\n \n\u003cp\u003eNot applicable.\u003c/p\u003e\n \n\u003ch2\u003eAvailability of data and materials\u003c/h2\u003e\n \n\u003cp\u003eThe data sets supporting the results of this article are included within the article\n and its additional files.\u003c/p\u003e\n \n\u003ch2\u003eCompeting Interests\u003c/h2\u003e\n \n\u003cp\u003eThe authors declare that they have no conflict of interest.\u003c/p\u003e\n \n\u003ch2\u003eFunding\u003c/h2\u003e\n \n\u003cp\u003eThe financial support from National Natural Science Foundation of China (No. 31460515),\n Jiangsu Agricultural Innovation of New Cultivars (No. PZCZ201719)，Guizhou Science\n and Technology Support Plan （NO.20192259）and Special Fund for Agro-Scientific Research\n of Guizhou Academy of Agricultural Sciences (No. 2015-19) are deeply appreciated.\n Guizhou Construction of Genetics and Breeding Laboratory of Vegetable Industry System(2019-0102).\u003c/p\u003e\n \n\u003ch2\u003eAuthor contributions \u003c/h2\u003e\n \n\u003cp\u003eJFC guide and designed the research. YD , BT, WYF, LT, LZ , XZ and JL performed the experiment, analyzed,\n interpreted date, and rote the manuscript. All authors reviewed and contributed to\n draft the manuscript. All authors read and approved the final manuscript.\u003c/p\u003e\n \n\u003ch2\u003eAcknowledgements \u003c/h2\u003e\n \n\u003cp\u003eThe authors would like to acknowledge the Dr Tan (Institute of Horticulture, Guizhou\n Academy of Agricultural Sciences, Guiyang , China) for critical reading of the manuscript\n , Jiaji Zhang( Personal Biotechnology Corporation, Shanghai, China) for their efforts in bioinformatics analysis and Springer Nature\n Author Services for writing services.\u003c/p\u003e\n "},{"header":"References","content":"\u003col\u003e \n\u003cli\u003eGuha S, Maheshwari SC. In vitro production of embryos from anthers of datura. nature. 1964;204:497.\u003c/li\u003e\n \n \n \n\u003cli\u003eMalik MR, Wang F, Dirpaul JM, Zhou N, Polowick PL, Ferrie, AMR, Krochko JE. Transcript\n profiling and identification of molecular markers for early microspore embryogenesis\n in \u003cem\u003eBrassica napus\u003c/em\u003e. Plant Physiology. 2007;144:134-154.\u003c/li\u003e\n \n \n \n\u003cli\u003eBhat JG, Murthy HN. Factors affecting \u003cem\u003ein-vitro\u003c/em\u003e gynogenic haploid production in niger (\u003cem\u003eGuizotia abyssinica\u003c/em\u003e(L. f.) Cass.). Plant Growth Regulation. 2007;52:241-248.\u003c/li\u003e\n \n \n \n\u003cli\u003eMetwally EI , Moustafa SA , El-Sawy BI , Haroun SA, Shalaby TA. Production of haploid\n plants from in vitro culture of unpollinated ovules of \u003cem\u003eCucurbita pepo\u003c/em\u003e. Plant Cell Tissue and Organ Culture. 1998;52:117-121.\u003c/li\u003e\n \n \n \n\u003cli\u003eGémes-Juhász A, Balogh P, Ferenczy AZ, Kristóf. Effect of optimal stage of female gametophyte and heat treatment on in vitro gynogenesis\n induction in cucumber (\u003cem\u003eCucumis sativus\u003c/em\u003e L.). Plant Cell Reports. 2002;21:105-111.\u003c/li\u003e\n \n \n \n\u003cli\u003eHalit Y, Nebahat S. A new method for haploid muskmelon (\u003cem\u003eCucumis melo \u003c/em\u003eL.) dihaploidization. Scientia Horticulturae. 2003;98:277–283.\u003c/li\u003e\n \n \n \n\u003cli\u003eShalaby TA. Factors affecting haploid induction through \u003cem\u003ein vitro\u003c/em\u003e gynogenesis in summer squash (\u003cem\u003eCucurbita pepo\u003c/em\u003e L.). Scientia Horticulturae. 2007;115:1-6.\u003c/li\u003e\n \n \n \n\u003cli\u003eDiao WP, Jia YY, Song H, Zhang XQ, Lou QF, Chen JF. Efficient embryo induction in\n cucumber ovary culture and homozygous identification of the regenetants using SSR\n markers. Scientia Horticulturae. 2009;119:246-251.\u003c/li\u003e\n \n \n \n\u003cli\u003eLi JW, Si SW, Cheng JY, Li JX, Liu JQ. Thidiazuron and silver nitrate enhanced gynogenesis\n of unfertilized ovule cultures of \u003cem\u003eCucumis sativus\u003c/em\u003e. Biologia Plantarum. 2013;57:164-168.\u003c/li\u003e\n \n \n \n\u003cli\u003ePlapung P, Khumsukdee S, Smitamana P. Development of cucumber lines resistant to cucumber mosaic\n virus by ovule culture. International Journal of Agricultural Technology. 2014;733-741.\u003c/li\u003e\n \n \n \n\u003cli\u003eTantasawat PA, Sorntip A, Poolsawat O, Chaowiset W, Pornbungkerd P. Evaluation of factors affecting\n embryo-like structure and callus formation in unpollinated ovary culture of cucumber(\u003cem\u003eCucumis sativus\u003c/em\u003e). International Journal of Agriculture \u0026amp; Biology. 2015;17:613-618.\u003c/li\u003e\n \n \n \n\u003cli\u003eTouraev A, Vicente O, Heberle-Bors E. Initiation of microspore embryogenesis by stress. Trends\n in Plant Science. 1997;2:297–302.\u003c/li\u003e\n \n \n \n\u003cli\u003ePerry SE, Fernandez DE. The MADS-domain protein \u003cem\u003eAGAMOUS-like 15\u003c/em\u003e accumulates in embryonic tissues with diverse origins. Plant Physiology. 1999;120:121-129.\u003c/li\u003e\n \n \n \n\u003cli\u003eBoutilier K, Offringa R, Sharma VK, Kieft H, \u0026nbsp;\u003ca href=\"http://med.wanfangdata.com.cn/Paper/Search?q=%e4%bd%9c%e8%80%85%3a(Th%c3%a9r%c3%a8se%2c+Ouellet)\"\u003eOuellet\u003c/a\u003e T, Zhang L, Hattori J, Liu CM, van Lammeren AAM, Miki BLA, Custers JBM, van Looleren\n Campagne MM. ectopic expression of \u003cem\u003eBABY BOOM\u003c/em\u003e triggers a conversion from vegetative to embryonic growth. Plant Cell. 2002;14:1737-1749.\u003c/li\u003e\n \n \n \n\u003cli\u003eMaraschin SF, Priester WD, Spaink HP, Spaink, Wang M. Androgenic switch: an example\n of plant embryogenesis from the male gametophyte perspective. Journal of Experimental\n Botany. 2005;56:1711-1726.\u003c/li\u003e\n \n \n \n\u003cli\u003eHosp J, Tashpulatov A, Roessner U, \u003ca href=\"http://med.wanfangdata.com.cn/Paper/Search?q=%e4%bd%9c%e8%80%85%3a(Ekaterina%2c+Barsova)\"\u003eBarsova\u003c/a\u003e E, \u003ca href=\"http://med.wanfangdata.com.cn/Paper/Search?q=%e4%bd%9c%e8%80%85%3a(Heidrun%2c+Katholnigg)\"\u003eKatholnigg\u003c/a\u003e H, \u003ca href=\"http://med.wanfangdata.com.cn/Paper/Search?q=%e4%bd%9c%e8%80%85%3a(Ralf%2c+Steinborn)\"\u003eSteinborn\u003c/a\u003e R, \u003ca href=\"http://med.wanfangdata.com.cn/Paper/Search?q=%e4%bd%9c%e8%80%85%3a(Bal%c3%a1zs%2c+Melikant)\"\u003eMelikant\u003c/a\u003e B, \u003ca href=\"http://med.wanfangdata.com.cn/Paper/Search?q=%e4%bd%9c%e8%80%85%3a(Sergey%2c+Lukyanov)\"\u003eLukyanov\u003c/a\u003e S, \u003ca href=\"http://med.wanfangdata.com.cn/Paper/Search?q=%e4%bd%9c%e8%80%85%3a(Erwin%2c+Heberle-Bors)\"\u003eHeberle-Bors\u003c/a\u003e E, \u003ca href=\"http://med.wanfangdata.com.cn/Paper/Search?q=%e4%bd%9c%e8%80%85%3a(Alisher%2c+Touraev)\"\u003eTouraev\u003c/a\u003e A. Transcriptional and metabolic profiles of stress-induced, embryogenic tobacco\n microspores. Plant Molecular Biology. 2007;63:137-149.\u003c/li\u003e\n \n \n \n\u003cli\u003eRosa Angélica SánchezDíaz, Ana María Castillo, María Pilar Vallés. Microspore embryogenesis\n in wheat: new marker genes for early, middle and late stages of embryo development.\n Plant Reproduction. 2013;26:287-296.\u003c/li\u003e\n \n \n \n\u003cli\u003eIrikova T, Grozeva S, Denev I. Identification of \u003cem\u003eBABY BOOM\u003c/em\u003e and\u003cem\u003e LEAFY COTYLEDON \u003c/em\u003egenes in sweet pepper (\u003cem\u003eCapsicum annuum\u003c/em\u003e L.) genome by their partial gene sequences. Plant Growth Regulation. 2012; 67:191-198.\u003c/li\u003e\n \n \n \n\u003cli\u003eBélanger S, Marchand S, Jacques PÉ, Meyers B, Belzile F. Differential expression profiling of microspores during the early stages of isolated\n microspore culture using the responsive barley cultivar gobernadora. Genes Genomes Genetics.2018;8:1603-1614\u003c/li\u003e\n \n \n \n\u003cli\u003eZhang J , Li J , Liu B , Zhang L, Chen J, Lu M. Genome-wide analysis of the populus \u003cem\u003eHSP90 \u003c/em\u003egene family reveals differential expression patterns localization and heat stress\n responses. BMC Genomics. 2013;14:532-546.\u003c/li\u003e\n \n \n \n\u003cli\u003eLv J , Qi J , Shi Q , Shen D, Zhang S, Shao G, Li H, Sun ZY, Weng Y, Gu XF, Li XX,\n Zhu XG, Zhang JZ, van Treuren R, van Dooijeweert W, Zhang ZH, Huang S. Genetic diversity\n and population structure of cucumber (\u003cem\u003eCucumis sativus\u003c/em\u003e L.). PLoS ONE. 2012;7:e46919.\u003c/li\u003e\n \n \n \n\u003cli\u003eQing ZW, Xin M, Liu D, Zhang YJ, Wang XG, Zhou XY. A new cucumber F\u003csub\u003e1\u003c/sub\u003e hybrid 'dongnong 812'. Chinese Vegetables. 2018;5:71-73.(in Chinese)\u003c/li\u003e\n \n \n \n\u003cli\u003eGolabadi M, Golkar P, Eghtedary AR. Assessment of genetic variation in cucumber (\u003cem\u003eCucumis sativus \u003c/em\u003eL.) genotypes. European Journal of Experimental Biology. 2012;2 :1382-1388.\u003c/li\u003e\n \n \n \n\u003cli\u003eWang X, Bao K, Reddy UK, Bai Y, Hammar SA, Jiao C, Wehner TC, Ramírez-Madera AO, Weng\n Y, Grumet R, Fei Z. The USDA cucumber (\u003cem\u003eCucumis sativus\u003c/em\u003e L.) collection: genetic diversity, population structure, genome-wide association\n studies, and core collection development. Horticulture Research. 2018;5:64-67\u003c/li\u003e\n \n \n \n\u003cli\u003eChen JF, Cui L, Malik AA, Kere George Mbira. \u003cem\u003eIn vitro\u003c/em\u003e haploid and dihaploid production \u003cem\u003evia\u003c/em\u003e unfertilized ovule culture. Plant Cell Tissue and Organ Culture. 2011; 104:311-319.\u003c/li\u003e\n \n \n \n\u003cli\u003eHoekstra S, Hoekstra IR, Hoekstra RA, Hoekstra E. The Interaction of 2, 4-D application\n and mannitol pretreatment in anther and microspore culture of \u003cem\u003eHordeum vulgare\u003c/em\u003e L. cv. Igri. Journal of Plant Physiology. 1996;148:696-700.\u003c/li\u003e\n \n \n \n\u003cli\u003eMaraschin, SF. Androgenic switch: an example of plant embryogenesis from the male gametophyte\n perspective. Journal of Experimental Botany. 2005;56:1711-1726.\u003c/li\u003e\n \n \n \n\u003cli\u003eIkeda-Iwai M, Umehara M, Satoh S, Kamada H. Stress‐ induced somatic embryogenesis in vegetative tissues of \u003cem\u003eArabidopsis thaliana\u003c/em\u003e. Plant Journal. 2003;34:107-114.\u003c/li\u003e\n \n \n \n\u003cli\u003eMól R, Filek M, Machackova I, Matthys-Rochon E. Ethylene synthesis and auxin augmentation\n in pistil tissues are important for egg cell differentiation after pollination in\n maize. Plant \u0026amp; Cell Physiology. 2004;45:1396-1405.\u003c/li\u003e\n \n \n \n\u003cli\u003ePagnussat GC, Alandetesaez M, Bowman JL, Sundaresan v. Auxin-dependent patterning\n and gamete specification in the \u003cem\u003eArabidopsis\u003c/em\u003e female gametophyte. Science. 2009;324:1684-1689.\u003c/li\u003e\n \n \n \n\u003cli\u003eFan Z, Armstrong KC, Keller WA. Development of microspores \u003cem\u003ein vivo\u003c/em\u003e, and \u003cem\u003ein vitro\u003c/em\u003e, in \u003cem\u003eBrassica napus\u003c/em\u003e L. Protoplasma. 1988;147:191-199.\u003c/li\u003e\n \n \n \n\u003cli\u003eGolabadi M, Ghanbari Y, Ghanbari S, Keighobadi K. Embryo and callus induction by different\n factors in ovary culture of cucumber. Journal of Applied Botany \u0026amp; Food Quality. 2017;90:68-75.\u003c/li\u003e\n \n \n \n\u003cli\u003eAušra B, Natalija B, Ramunė M. Genotypic and exogenous factors affecting linseed ovary culture. Zemdirbyste. 2017;104:243-248.\u003c/li\u003e\n \n \n \n\u003cli\u003eHahn A, Bublak D, Schleiff E, Scharf KD. Crosstalk between \u003cem\u003eHsp90\u003c/em\u003e and \u003cem\u003eHsp70\u003c/em\u003e Chaperones and Heat Stress Transcription Factors in Tomato. Plant Cell. 2011; 23:741-755.\u003c/li\u003e\n \n \n \n\u003cli\u003eSangster TA, Queitsch C. The \u003cem\u003eHSP90\u003c/em\u003e chaperone complex, an emerging force in plant development and phenotypic plasticity.\n Current Opinion in Plant Biology. 2005;8:86-92.\u003c/li\u003e\n \n \n \n\u003cli\u003eZarsky V, Garrido D, Eller N, Eller N, Tupy J, Vicente O, Schoffl F. The expression\n of a small heat shock gene is activated during induction of tobacco pollen embryogenesis\n by starvation. Plant Cell \u0026amp; Environment. 2010;18:139-147.\u003c/li\u003e\n \n \n \n\u003cli\u003eBinarova P, Hause G, Cenklová V, Cordewener JHG, Looderen Campagne MM. A short severe\n heat shock is required to induce embryogenesis in late bicellular pollen of \u003cem\u003eBrassica napus\u003c/em\u003e L.. Sexual Plant Reproduction. 1997;10:200-208.\u003c/li\u003e\n \n \n \n\u003cli\u003eSabehat A, Weiss D. Expression of small heat-shock proteins at low temperatures: A\n possible role in protecting against chilling injuries. Plant Physiology. 1998;117:651-8.\u003c/li\u003e\n \n \n \n\u003cli\u003eMaraschin SF. Androgenic switch: an example of plant embryogenesis from the male gametophyte\n perspective. Journal of Experimental Botany. 2005; 56:1711-1726.\u003c/li\u003e\n \n \n \n\u003cli\u003ePechan PM, Bartels D, Brown DC, Schell J. Messenger-RNA and protein changes associated\n with induction of \u003cem\u003eBrassica\u003c/em\u003e, microspore embryogenesis. Planta. 1991;184:161-165.\u003c/li\u003e\n \n \n \n\u003cli\u003eHalperin W, Jensen WA. Ultrastructural changes during growth and embryogenesis in\n carrot cell cultures. Journal of Ultrastruct Research. 1967;18:428-443.\u003c/li\u003e\n \n \n \n\u003cli\u003eHuang BQ, Sheridan WF. Female Gametophyte Development in maize: microtubular organization\n and embryo sac polarity. Plant Cell. 1994;6:845-861.\u003c/li\u003e\n \n \n \n\u003cli\u003eStone SL, Kwong LW, Yee KM, Yee KM, Pelletier J, Lepiniec L, Fischer RL. \u003cem\u003eLEAFY COTYLEDON2\u003c/em\u003e encodes a B3 domain transcription factor that induces embryo development. Proceedings of the National Academy of Sciences. 2001;98:11806-11811.\u003c/li\u003e\n \n \n \n\u003cli\u003eLotan T, Ohto M, Yee K M, West MA, Lo R, Kwong RW, Yamagishi, Fischer RL, Goldberg\n RB, Harada JJ. \u003cem\u003eArabidopsis\u003c/em\u003e \u003cem\u003eLEAFY COTYLEDON1\u003c/em\u003e is sufficient to induce embryo development in vegetative cells. Cell. 1998;93:1195-1205.\u003c/li\u003e\n \n \n \n\u003cli\u003eRiechmann JL, Meyerowitz EM. The \u003cem\u003eAP2/EREBP\u003c/em\u003e family of plant transcription factors. Biological Chemistry. 1998;379:633-646.\u003c/li\u003e\n \n \n \n\u003cli\u003eMayer KF, Schoof H, Haecker A, Lenhard M, Laux T. Role of \u003cem\u003eWUSCHEL\u003c/em\u003e in Regulating Stem Cell Fate in the \u003cem\u003eArabidopsis\u003c/em\u003e Shoot Meristem. Cell. 1998;95:805-815.\u003c/li\u003e\n \n \n \n\u003cli\u003ePalovaara J, Hakman I. Conifer WOX-related homeodomain transcription factors, developmental\n consideration and expression dynamic of \u003cem\u003eWOX2\u003c/em\u003e during Picea abies somatic embryogenesis. Plant Molecular Biology. 2008;66:533-549.\u003c/li\u003e\n \n \n \n\u003cli\u003eMao X, Cai T, Olyarchuk JG, Wei L. Automated genome annotation and pathway identification\n using the KEGG Orthology (KO) as a controlled vocabulary. Bioinformatics. 2005;21:3787–3793.\u003c/li\u003e\n\u003c/ol\u003e"},{"header":"Additional file legend","content":"\u003cp\u003eTable S1. Statistics of RNA-Seq alignment\u003c/p\u003e\n \n\u003cp\u003eTable S2. GO classification of common expressed genes in gynogenesis stage.\u003c/p\u003e\n \n\u003cp\u003eTable S3. GO classification of common expressed genes from embryogenesis to the mature stage.\u003c/p\u003e\n \n\u003cp\u003eTable S4. GO classification of common expressed genes in organogenesis stage.\u003c/p\u003e\n \n\u003cp\u003eTable S5. Primers for Real-time quantitative PCR\u003c/p\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":true,"highlight":"","institution":"","isAcceptedByJournal":false,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true},"keywords":"Cucumber, Ovary culture, Embryogenesis, Transcriptome, Differentially expressed genes","lastPublishedDoi":"10.21203/rs.2.10287/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.2.10287/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"Background: Ovary culture has been a useful way to generate double haploid (DH) plant in cucumber (Cucumis sativus L.). However, the rate of embryo induction is low, and the probability for the induced embryo to grow into normal embryo is low as well. This is largely due to unknown of the mechanism of embryogenesis in cucumber. In this study, the differentially expressed genes during embryogenesis, including the early stages of embryo formation, embryo maturation and shoot formation, was investigated with transcriptomic technique to set up basis for a more efficient ovary culture technology\nResults: Cytological observations led to suggestions that cell enlargement is the symbol that gametophytes had switched to the sporophyte development pathway during the early embryogenesis stage. In this stage, RNA-seq revealed 3468 up-regulated genes, 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.\nThe total of 480 genes that function in protein complex binding, microtubule binding, tetrapyrrole binding, tubulin binding and other microtubule activities were continuously up-regulated during the embryo maturation stage, indicating that the cytoskeleton structure was continuously being built and maintained by the action of microtubule-binding proteins and enzyme modification during embryo development.\nIn the shoot formation stage, 1383 genes were up-regulated, which were mainly enriched in phenylpropanoid biosynthesis, plant hormone signal transduction, phenylalanine metabolism, and starch and sucrose metabolism. The shoot formation stage might be regulated by 6 transcription factors that contained a B3 domain, 9 genes in the AP2/ERF family and 2 genes encoded WUS homologous domain proteins.\nConclusions: Findings from this study offer a valuable framework for explaining the transcriptional regulatory mechanism underlying embryogenesis in cucumber ovary culture.","manuscriptTitle":"Analysis of differentially expressed genes during embryogenesis in ovary culture of cucumber (Cucumis sativus L.)","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2019-06-13 18:12:48","doi":"10.21203/rs.2.10287/v1","editorialEvents":[{"type":"communityComments","content":1}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"582b17bf-30cd-46e9-9c30-5b9faf5c7455","owner":[],"postedDate":"June 13th, 2019","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"posted","subjectAreas":[{"id":13521,"name":"Plant Molecular Biology and Genetics"}],"tags":[],"updatedAt":"","versionOfRecord":[],"versionCreatedAt":"2019-06-13 18:12:48","video":"","vorDoi":"","vorDoiUrl":"","workflowStages":[]},"version":"v1","identity":"rs-1316","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"identity":"rs-1316","version":["v1"]},"buildId":"WrCJVZZCHTDjtuVLN7oU0","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}