Skip to main content

Exogenous spraying of IAA improved the efficiency of microspore embryogenesis in Wucai (Brassica campestris L.) by affecting the balance of endogenous hormones, energy metabolism, and cell wall degradation



Microspore embryogenesis is an extraordinarily complicated process, comprehensively regulated by a composite network of physiological and molecular factors, among which hormone is one of the most crucial factors. Auxin is required for stress-induced microspore reprogramming, however, the mechanism of its regulation of microspore embryogenesis is still unclear.


In this study, we found exogenously spraying 100 mg·L− 1 IAA on the buds of Wucai significantly increased the rate of microspore embryogenesis, and moreover accelerated the process of embryogenesis. Physiological and biochemical tests showed that the contents of amino acids, soluble total sugar, soluble protein, and starch were significantly increased after IAA treatment. Furthermore, exogenously spraying 100 mg·L− 1 IAA significantly enhanced IAA, GA4, and GA9 content, increased catalase (CAT) and malondialdehyde (MDA) activity, and reduced abscisic acid (ABA), MDA and soluble protopectin content, H2O2 and O2· production rate in the bud with the largest population of late-uninucleate-stage microspores. Transcriptome sequencing was performed on buds respectively treated with 100 mg·L− 1 IAA and fresh water. A total of 2004 DEGs were identified, of which 79 were involved in micropores development, embryonic development and cell wall formation and modification, most of which were upregulated. KEGG and GO analysis revealed that 9.52% of DEGs were enriched in plant hormone synthesis and signal transduction pathways, pentose and glucuronic acid exchange pathways, and oxidative phosphorylation pathways.


These findings indicated that exogenous IAA altered the contents of endogenous hormone content, total soluble sugar, amino acid, starch, soluble protein, MDA and protopectin, the activities of CAT and peroxidase (POD), and the production rate of H2O2 and O2·. Combined with transcriptome analysis, it was found that most genes related to gibberellin (GA) and Auxin (IAA) synthesis and signal transduction, pectin methylase (PME) and polygalacturonase (PGs) genes and genes related to ATP synthesis and electron transport chain were upregulated, and genes related to ABA synthesis and signal transduction were downregulated. These results indicated that exogenous IAA treatment could change the balance of endogenous hormones, accelerate cell wall degradation, promote ATP synthesis and nutrient accumulation, inhibit ROS accumulation, which ultimately promote microspore embryogenesis.

Peer Review reports


Wucai (Brassica campestris L. ssp. chinensis var. rosularis Tsen) is a variant of non-heading Chinese cabbage (Brassica campestris L.), which is an important species in the Brassicaceae family [1]. As an important autumn and winter vegetable crop, this crop is cultured widely in the Yangtze-Huaihe River Basin. With a high nutritional value, Wucai is becoming an increasingly popular crop in other countries [2]. Breeding homozygous parental lines using traditional plant-breeding methods is difficult because Brassicaceae family crops are cross-pollinated. Microspore embryogenesis is the most effective tool of obtaining doubled haploids (DH) plants that are utilized in plant breeding for accelerating production of new cultivars [3].

Microspore embryogenesis involve a complex network of factors [4]. There are a number of intrinsic and extrinsic factors that affect the embryogenesis frequency, including the genotype and physiological status of donor plants, the developmental stage of the microspore, cold pretreatment and heat shock, culture medium (pH, carbohydrates, nitrogen source and other additives) and plant hormones [5]. Although the majority of microspore embryogenesis do not require exogenous hormones as inducers, a transient stress is essential to reprogram of microspores, probably as a response mediated by endogenous hormones.

Among plant growth regulators, auxin is the most significant (p < 0.05) hormone in plant development and it is the key regulator of cell division and differentiation [6]. During embryogenesis in planta, the major form of auxin, indoleacetic acid (IAA), has demonstrated key functions, particularly in embryo patterning [7]. The synthesis, activity, and polar transport of indoleacetic acid (IAA) were found to be necessary for the initiation and continuation of microspore embryogenesis in monocotyledonous barley and dicotyledonous Brassica napus [8]. Furthermore, gibberellin is also a key regulator of plant growth and development. The two hormones usually work together to control cell division and differentiation. In microspore cultures of B. napus, GA3 improved plantlet regeneration, mainly via elongation of the embryo axis and acceleration of its maturation [9]. Inhibitor of GA-biosynthesis, uniconazole, applied to B. napus embryo at the globular stage of development significantly reduced axis elongation [10]. Besides auxins and gibberellins, abscisic acid (ABA), known as a ubiquitous plant stress hormone, has been claimed to play a role in microspore embryogenesis (ME)-inducing signal transduction system [11]. In recent studies, the function of ABA synthesis and signaling during the in vitro embryogenic process has been analyzed, with Abscisic Acid Insensitive 3 (ABI3) and ABI4 transcription factors shown to play a relevant role in embryo formation [12]. High frequency of plant development could be obtained by abscisic acid (ABA) treatment of embryos [13]. On the contrary, higher level of ABA significantly inhibited embryogenesis.

Plant cell wall is complex and dynamic structures that involved in a series of functions such as interaction with the environment, intercellular communication, regulation of plant growth and development, and the determination of cell size and shape [14]. Previous studies indicated that cell wall remodeling plays an indispensable role in microspore embryogenesis [15]. Pectin are the family of oligosaccharides and polysaccharides with common characteristics and the major component of the primary wall [16], which play a vital role in microspore embrogenesis [15]. Pectin exist in the cell wall with highly methylated form, and its methylation level is mainly regulated by pectin methyl esterases (PMEs) [17]. PMEs catalyze the removal of methyl esters from homogalacturonan backbone of pectins in plant cell wall. The de-methylesterified residues can either form Ca2+ bonds or turn into a target for polygalacturonases, one of pectin-degrading enzymes, which affect the texture and rigidity of the cell wall [18]. In Brassica napus, Changes in pectin esterification levels are relevant to temporal expression patterns of PME genes in microspore embryogenesis [14]. De-methylesterified pectins are enriched in cell and PME genes are highly expressed in torpedo and cotyledon embryos of Brassica napus. Furthermore, inhibition of PME activity led to impair somatic embryogenesis in Quercus suber[19], indicating that pectins esterification play a role in inducing embryogenesis. Testillano et al. considered that PME-mediated configuration variation of pectins could be a crucial factor for microspore mebryogenesis, acting by promoting the cell wall remodeling in the process [15].

Intracellular reactive oxygen species (ROS) levels have different degrees of influence on microspore development. In previous studies, ROS were involved in the programmed cell death (PCD) of tapetal cells [20, 21]. The abnormal metabolism of ROS during plant PCD leads to delayed tapetal degradation, insufficient nutrition supply for microspore development, and microspore abortion [22,23,24]. To maintain the balance of intracellular ROS, plants have evolved a series of clear mechanisms for enzymatic and non-enzymatic ROS [25]. The activity of antioxidant enzymes plays an important role in clarifying free radicals and regulating cell tissue activity during embryogenesis [26].

In this study, the optimal concentration of exogenous IAA was identified in Wucai, which significantly (p < 0.05) increased the rate of microspore embryogenesis and moreover accelerated the process of embryogenesis. RNA-seq analysis suggested that plant hormone synthesis and signal transduction pathways, pentose and glucuronic acid exchange pathways, and oxidative phosphorylation pathways were involved in the regulation of microspores embryogenesis. Combined with analysis of changes in endogenous hormones and physiological indicators and transcriptome analysis, it was found that exogenous IAA treatment could change the balance of endogenous hormones, accelerate cell wall degradation, promote ATP synthesis and nutrient accumulation, inhibit ROS accumulation, which ultimately promote microspore embryogenesis.


Effects of IAA on the rate of microspore embryogenesis

The rate of isolated microspore-derive embryogenesis in Wucai was significantly (p < 0.05) increased treated with IAA (Fig. 1). Statistical analysis showed that IAA treatment enhanced the occurrence of embryoid bodies in a dose dependent manner (Table 1). Low concentration of exogenous IAA had little or no effect on the incidence of microspore embryogenesis, while high concentration even inhibited microspore embryogenesis. This experiment found that 100 mg. L − 1 IAA treatment had the best effect on induction of embryogenesis, the embryo yield reached 16.4 embryos. bud − 1, 3.15 times higher than that in the control group. At the same time, 150 mg. L − 1 IAA inhibited the development of isolated microspore embryos. These results indicate that treatment with an appropriate IAA concentration is conducive to the embryogenesis of isolated microspores in Wucai.

Fig. 1
figure 1

Images showing the embryogenic potential of isolated microspores of Wucai treated with different concentrations of IAA: (a) 0 mg·L− 1 IAA; (b) 50 mg·L− 1 IAA; (c) 75 mg·L− 1 IAA; (d) 100 mg·L− 1 IAA; (e) 125 mg·L− 1 IAA; and (f) 150 mg·L− 1 IAA. Bar = 1 cm

Cytological observation of microspore embryonic development

According to the first division of microspore, the embryonic development pathways were divided into two types: Route A (unequal division) and Route B (equal division). Route A can be further subdivided according to the second division and subsequent situations. Studies had shown that there were two pathways, A and B, for embryonic development in isolated microspore cultures of non-heading Chinese cabbage, and pathway A was more common [27].

In this study, the microspore mother cells underwent meiosis to form a tetrad (Fig. 2), and four parenchyma cells were surrounded by a thick layer of callose (Fig. 2A, a). The tetrad was separated to form isolated microspores. The microspore just separated from the tetrad was in the mononuclear metaphase, when the nuclear to cytoplasm ratio was large (Fig. 2B, b). Subsequently, the nucleus was pushed aside by vacuoles, i.e., the late uni-nucleated stage (Fig. 2C, c), the most sensitive stage for the induction of microspore embryogenesis. The microspore undergoes mitosis and forms two cell nucleus, one large and one small (Fig. 2D, d), i.e. unequal division. Therefore, we speculated that the developmental pathway of isolated microspores in Wucai was pathway A.

Fig. 2
figure 2

Cytological observation of isolated microspore development between the control (A–F) and IAA-treated (a–f) groups. A(a) tetrad stage (0 d). B(b) early uni-nucleated microspore stage (0 d). C(c) late uni-nucleated microspore stage (0 d). D(d) bi-nucleated microspore stage (3 d). E(e) Multicellular proembryo (5 d). F(f) Embryoid structure (10 d). G(g) Embryoid (15 d). msp, microspore, cy, cytoplasm, n, nucleus

Compared to the control, the IAA treated microspores had a dispersed arrangement of four microspores during the tetrad stage (Fig. 2b). At the mid-mononucleate stage, the nucleus was significantly larger, and at the mononucleate leaning stage, it was triclinic with a distinct budding furrow, in preparation for subsequent divisions (Fig. 2c). The late uni-nucleate stage (Fig. 2d) was the spmpling stage of the buds in our experiment, as it was the best stage for inducing microspore embryogenesis [28, 29]. In the 3rd day of isolated microspore culture, the microspores were observed to undergo unequal division, forming a large vegetative nucleus and a small reproductive nucleus (Fig. 2d). On the 5th day of culture, the multicellular clusters formed were regular and round. (Fig. 2e). On the 10th day, in the control group, the early embryos were mostly torpedo shaped (Fig. 2F, G), while in the IAA treated group, the early embryos first formed a long stalk and eventually developed into cotyledons (Fig. 2f, g). The multicellular proto embryos divided further in culture and eventually formed microspore embryos. In addition, the related research indicated that the microspore embryo would pass through the globular embryo, the heart shaped embryo, the torpedo shaped embryo and finally develop into the cotyledonary embryo [28]. In this study, microscopic observation of microspores after 5 days of culture and statistical analysis of 10 randomly selected fields of view showed that the number of multicellular proembryos treated with IAA was significantly higher than that of the control group (fresh water treatment) (Additional file 1: Fig. s1) The initial time of microspore embryogenesis was counted after spraying with different concentrations of IAA (Additional file 2: Table s1). Compared with the control group, the time of microspore embryogenesis was 4 d earlier after 100 mg. L− 1 IAA treatment.

Comparison of changes in the content of endogenous plant hormones

After exogenous IAA treatment, 0.3 g of fresh flower bud samples were used to determine the endogenous levels of plant hormones such as auxin, CKs, ABA, GAs, JA and SLs. Compared to the control, the content of GAs in monocotyledonous leaning stage flower buds was dramatically increased, and the ABA content was significantly (p < 0.05) decreased after exogenous IAA treatment, while the content of other phytohormones did not show obvious (p < 0.05) changes (Fig. 3).

Fig. 3
figure 3

Determination of endogenous plant hormone content between buds of the control a nd IAA treatment groups. (A) Auxin content, (B) cytokinin content, (C) gibberellin content, (D) abscisic acid content, (E) jasmonic acid content, and (F) monolayer content. Data were derived from three independent experiments. Error bars represent ± SD. Different letters indicate significant differences (p < 0.05)

In the plant system, the physiological effects of different categories of plant hormones have a mutually reinforcing or antagonistic effect, and the regulation of plant growth and development is usually the result of the combined action of multiple hormones. Therefore, the ratio of several hormones was analyzed in this experiment. The ratios of (auxin + GAs)/ABA, (auxin + GAs + CKs)/ABA, GAs/ABA, and auxin/ABA were all higher in the treatment group than in the control group (Addition File 3: Table s2), while the ratios of (auxin + GAs)/CKs and auxin/CKs were not obviously different from those in the control group. These results indicate that microspore development may also be related to the relative content and the dynamic balance between hormones.

Transcriptome sequencing data and quality control

To explore the molecular mechanism of exogenous IAA in enhancing the embryogenesis of isolated microspores in Wucai, we performed transcriptomic analysis using the buds of the largest population of late-uninucleate-stage were analyzed after treatment with 100 mg. L− 1 IAA and the same dose of water for 24 h. After sequencing quality control and data analysis, 256,144,226 clean reads were obtained, and the clean data of each sample reached more than 5 GB. The percentage of Q30 was more than 92% (Table 2). In a direct comparison of the density and discrete distribution of expression levels for different samples, the sequencing quality and gene expression levels were the same (Additional file 4: Fig s2). These indicated that the sequencing quality was good, and could be analyzed in the next step.

Table 1 Effect of exogenous IAA on embryo rate of isolated microspores
Table 2 Illumina sequencing data and results of de novo assembly

Identification of differentially expressed genes

To better understand the expression patterns of differentially expressed genes (DEGs), DESeq was used for differential expression analysis [30, 31]. DEGs between the buds treated with 100 mg. L− 1 IAA and control group were identified. The screening conditions of differential genes were as follows: absolute value of log2(FoldChang) greater than or equal to 1, and the p-value less than 0.05. A total of 2004 DEGs were identified, comprising 1002 upregulated genes and 1002 downregulated genes. In addition, genes that were not compared to the reference genome were considered new, and we identified 345 new genes, 52 of which could be not annotated to any database.

GO and KEGG analysis of differentially expressed genes

Detailed GO term annotations of the DEGs were categorized according to biological processes, molecular functions, and cellular components. The GO categories for the set of DEGs (Fig. 4) revealed that most of the encoded products were associated with “metabolic process”, “cellular process”, “multicellular BPs”, “developmental process”, “cellular anatomical entity”, “catalytic activity” and “binding”.

Fig. 4
figure 4

GO enrichment of flower bud DEGs between the control and IAA-treated groups

KEGG [32] ( ) enrichment analysis showed that 1334 differential genes were enriched in 121 metabolic pathways (Additional file5, Fig s3). The metabolic pathways were annotated with “metabolic pathway”, “pentose and glucuronate interconversions”, “plant-pathogen interaction”, “biosynthesis of secondary metabolism”, “plant hormone signal transduction” and “oxidative phosphorylation” (Fig. 5). These results indicated that a variety of complex metabolic pathways were involved in microspore development.

Fig. 5
figure 5

First 50 pathway entries for KEGG enrichment of bud differentially expressed genes (DEGs) between the control and IAA-treated groups

Pollen and embryo development related genes

Pollen development was a complex process involved many events, including microspore embryonic development. Many genes related to pollen and embryonic development have been identified in the model plant Arabidopsis thaliana [33, 34]. A total of 79 DEGs were identified as being associated with pollen and embryonic development, which 46 were upregulated and 23 were downregulated (Additional file 6; Table s3), including GDSL lipase genes, BHLH transcription factor family genes, AP2/ERF transcription factor family genes, pectinesterase genes, WUS-related genes, GST (glutathione S-transferase) related genes, and RALF (protein RALF-like) genes. These results indicate that these DEGs participate in pollen and embryonic development, as well as reproductive development, and play a certain role in promoting the subsequent embryogenesis of isolated microspores.

Based on GO and KEGG enrichment analysis and identification of development-related genes, we screened out 139 genes under the condition of |log2FC|>=2, P < 0.05, and presented them in the heatmap (Additional file 7: Fig s4).

Transcriptional and metabolic regulatory networks for plant hormone synthesis and signal transduction

To explore the role of the regulatory network of genes and metabolites in plant hormone synthesis and signal transduction, all detected DEGs and metabolites were mapped to the related pathways (Fig. 6A, B).

Fig. 6
figure 6

Metabolic pathway of plant hormone synthesis and signal transduction between the control and IAA-treated groups. (A) Plant hormone synthesis and signal transduction metabolic pathway (ko04075, ko00904, ko00308, ko00906). (B) Heat map of differentially expressed genes (DEGs) related to plant hormone synthesis and signal transduction. (C) IAA content. (D) ABA content. (E) GA4 content. (F) GA9 content. Data were derived from three independent experiments. Error bars represent ± SD. Different letters indicate significant differences (p < 0.05)

In the process of gibberellin synthesis (ko00904), non-bioactive GA24 was converted into active GA4 and GA9 under the catalysis of the GA20OX enzyme. GA2OX, as a key enzyme in the degradation process of GAs in plants, oxidizes bioactive GA4 to inactive GA34, and oxidizes GA9, the precursor of GA4, to inactive GA51 and GA29, thereby maintaining the dynamic balance between bioactive GAs and their intermediates in plants. During phytohormone signaling (ko04075), the transcription factors PIF3 and PIF4 were activated to mediate the gibberellin signal response. GA4 (Fig. 6E) and GA9 (Fig. 6F) contents were remarkably higher than that of the control group.

During auxin biosynthesis (ko00380), tryptophan is catalyzed by enzymes to synthesize indole 3 acetamide (IAM) and indole 3 acetaldehyde (IAAld), prerequisite substances of IAA. IAA1d is synthesized into indoleacetic acid through the action of aldehyde dehydrogenase. In our study, CYP71AB was downregulated, and DAO was upregulated. There were four genes encoding aldehyde dehydrogenase (ALDH), of which one was upregulated and three were downregulated. During auxin signal transduction (ko04075), more than 60% of DEGs were upregulated. At the same time, we found that the content of IAA was significantly (p < 0.05) higher than that in the control group (Fig. 6C). Changes in these DEGs may have contributed to the accumulation of IAA.

During ABA synthesis (ko00906), β-cyclohydroxylase (LUT5), β-carotene isomerase (DWARF27), and abscisic β-glucosyltransferase (AOG) were all downregulated, while abscisic acid 8’-hydroxylase (CYP707A) was upregulated. Leading to a decrease in ABA accumulation in vivo. Both the abscisic acid receptor (PYL) and protein phosphatase (PP2C) were downregulated during signal transduction (ko04075). This was consistent with our detection of a significant reduction (p < 0.05) in ABA content (Fig. 6D).

Oxidative phosphorylation pathway

Normal microspore development requires an adequate energy supply, and the energy produced by mitochondrial oxidative phosphorylation (ko00190) in plants supplies the needs for plant growth and breeding. We plotted a model of oxidative phosphorylatipn metabolism (Fig. 7A). A total of 17 DEGs were detected in the oxidative phosphorylation pathway, of which 14 were upregulated and 3 were downregulated (Fig. 7B). DEGs were enriched in NADH dehydrogenase, succinate dehydrogenase and ATP synthase (complex V), which were involved in electron transfer and ATP synthesis. Exogenous IAA treatment significantly (p < 0.05) decreased H2O2 (Fig. 7C), MDA (Fig. 7D), and O2· (Fig. 7E) contents, and enhanced POD (Fig. 7F), CAT (Fig. 7G) activity, which maintained the balance of ROS levels during microspore development and prevented ROS accumulation, reducing membrane damage.

Fig. 7
figure 7

Metabolic pathways of oxidative phosphorylation between the control and IAA-treated groups. (A) Oxidative phosphorylation metabolic pathway (ko00190). (B) Heat map of oxidative phosphorylation metabolism-related differentially expressed genes (DEGs). (C) H2O2 content. (D) MDA content. (E) O2. content. (F) POD content. (G) CAT content. Data were derived from three independent experiments. Error bars represent ± SD. Different letters indicate significant differences (p < 0.05)

To further prove that the development of microspores was related to the energy supply, the energy substances and nutrient substances were measured. Compared with the control group, the total soluble sugar (Fig. 8A), amino acid (Fig. 8B), starch (Fig. 8C) and soluble protein content (Fig. 8D) in buds dramatically increased (p < 0.05) after treated with IAA.

Fig. 8
figure 8

Determination of physiological indices. (A) Total soluble sugar content, (B) glucose content, (C) starch content, and (D) soluble protein content. Data were derived from three independent experiments. Error bars represent ± SD. Different letters indicate significant differences (p < 0.05)

Pentose and glucuronate interconversions pathway

According to the KEGG analyses, the pentose and glucuronic acid exchange pathway (ko00040) was significantly enriched. The regulatory pattern of related differential genes in the metabolic pathway was shown in Fig. 9A. In the pentose and glucuronate interconversion pathway, more than 80% upregulated expression of DEGs were detected, including the pectin methylase (PME) and polygalacturonase (PGs), key enzymes in pectin metabolism, encoding polygalacturonases (GSVIVT) and UDP-sugar pyrophosphorylase (USP) (Fig. 9C). The pectins content were significantly (p < 0.05) lower than that of the control (Fig. 9B).

Fig. 9
figure 9

Pentose and glucuronide exchange pathways between the control and IAA-treated groups. (A) Pentose and glucose exchange pathway (ko00040). (B) Heat map of the differentially expressed genes (DEGs) associated with the pentose and glucuronic acid interchange pathway. (C) Pectin content. Data were derived from three independent experiments. Error bars represent ± SD. Different letters indicate significant differences (p < 0.05)

The expression patterns of DEGs are verified by RT-qPCR

16 DEGs were selected from the important metabolic pathways for RT-qPCR to verify the reliability of RNA-Seq (Fig. 10). These DEGs were involved in plant hormone synthesis and signal transduction (IAA26, D27, DAO, GA2OX2, CYP71A13, ABAH3, C7352, GH33, CED3), pentose and glucuronic acid exchange pathways (PME5, PME1, PGLR1, PLY, BP19), oxidative phosphorylation metabolic pathways (SDHB, SDB19). The gene expression patterns obtained from RT-qPCR and RNA-Seq data showed similar trends, which confirmed the accuracy of the RNA-Seq results obtained in this study.

Fig. 10
figure 10

Expression levels of DEGs according to RT-qPCR (histogram) and RNA-Seq (black line chart). The FPKM values are based on the RNA-Seq data. The data obtained by RT-qPCR represent the means based on three replicates


Treatment with exogenous IAA to initiate microspore embryogenesis by affecting the balance of endogenous hormones

Hormones are an important factor affecting the development of microspore embryos. Auxin biosynthesis, transport, and response severely affect embryo development [35] IAA can be synthesized via the IPA pathway, the IAOx pathway, and other pathways [36]. The IAOx pathway may be a species-specific pathway for Brassicaceae species [37,38,39]. Premature degeneration of the tapetum during microspore development leads to cytoplasmic male sterility, and the ALDH gene has the function of overcoming this phenomenon [40]. The primary function of CYP71A13 is to catalyze the conversion of IAOx to IAN [41]. The regulation of active plant hormone levels is possible at β-glucosidases (BGLU) [42]. In this study, CYP71A13 was downregulated, and BGLU was upregulated; one gene encoding ALDH was upregulated. During plant hormone signal transduction, SAUR, and AUX/IAA, which are annotated as auxin-responsive proteins, are upregulated with high signal transduction levels. Related research shows that AUX prevents auxin leakage by transporting extracellular auxin back to the cytoplasm [43]. This may be the reason why the endogenous IAA content is significantly higher than that of the control group.

Studies have shown that exogenous IAA applications can improve endogenous GA biosynthesis [44]. The production of embryonic cells increased as GA content increased, and the amount of bioactive GA positively correlated with the formation of somatic embryos [45]. After treatment with exogenous 100 mg. L− 1 IAA, the contents of endogenous GAs, GA4and GA9 were significantly increased. In the GA biosynthesis process, the expression of GA20OX and GA2OX regulates the level of GA with biological activity in vivo. The molecular mechanism of GA on microspore embryogenesis needs to be further studied.

ABA accumulation inhibits the occurrence of microspore embryos, and low concentration of ABA promotes the occurrence of microspore embryos [46]. The absence of LUT5 results in a decrease in beta-carotene-derived lutein [47]. DWARF27 is a protein that converts trans-carotene to 9-cis-carotene [48]. 9-Cis-epoxycarotenoid dioxygenase (NCED) plays a key role in the development of ABA synthesis in plants [49]. Abscisic acid 8’-hydroxyase (CYP707A) is a key enzyme for ABA catabolism [50]. ABA seems to play a predominant role in the conversion of environmental signals into changes in the gene expression of plants. In this study, DWARF27, LUT5, NCED, and abscisic beta-glucosyltransferase (AOG) were downregulated, while CYP707A was upregulated. Both abscisic acid receptors (PYR/PYL) and protein phosphatase 2 C (PP2C) were downregulated, i.e., so ABA synthesis and transcription levels were low. This was consistent with significant decrease in ABA content.

The balanced relationship between plant hormones affects microspore development and embryogenesis [51]. In this study, the ratio of endogenous hormone contents in flower buds at the uni-nuclear stage after treatment with 100 mg. L− 1 IAA was statistically analyzed (Additional file 3: Table s2), and changes in hormone balance were found, which were consistent with the results of previous studies [51, 52]. Based on the analysis of plant hormone synthesis and signal transduction pathways, as well as the results of endogenous plant hormone content determination, we speculated that exogenous IAA affected the occurrence of isolated microspore embryos by affecting the level of endogenous IAA, the content of bioactive GAs, and the level of ABA.

Sufficient energy supply is key to microspore development

The energy supply is crucial for microspore development. Mitochondria are important sites for energy production [53]. The electron transport chain in the mitochondrial membrane is the main site of ROS production [54]. ROS are involved in the regulation of reproductive growth in plants. During microspore development, abnormal ROS accumulation leads to programmed cell death. H2O2, MDA, O2·, CAT, and POD are important indicators of the response to ROS levels. H2O2 stimulates the formation of embryonic cells [55]. Several studies have shown that H2O2 accumulation can produce a stress response, which has a promotional effect on embryonic development. The level of antioxidant enzyme activity plays an important role in scavenging free radicals and regulating cellular tissue activity during embryogenesis[26]. In this study, the results showed that the contents of MDA, O2·, and H2O2 were lower than those of the control group, and the activities of CAT and POD were higher than those of the control group, which were consistent with the results of previous studies on crops such as platycodon grandiflorum [56, 57]. This indicated that production and removal of reactive oxygen species in the bud were balanced and the damage caused by membrane lipid peroxidation was prevented.

NADH dehydrogenase and succinate dehydrogenase (SDH) play important roles in mitochondrial electron transport. NADH dehydrogenase complex 1 was an enzyme that catalyzes the transfer of electrons from NADH to CoQ at the inner mitochondrial membrane, and was the “entry enzyme” for mitochondrial oxidative phosphorylation [58]. Succinate dehydrogenase (SDH) was one of the hubs that connects electron transport and oxidative phosphorylation, providing electrons for the respiratory chain [59]. Higher plant cells had two F1-F0 ATPase complexes that were involved in mitochondrial oxidative phosphorylation. The F1F0-ATP complex enzyme can catalyze the hydrolysis or synthesis of ATP by adjusting the electrochemical potential gradient through proton penetration [60]. The accumulation of cytoplasmic PPI was toxic and can lead to serious growth defects and even cell death [61]. PPAs were enzymes that mainly participate in the hydrolysis of PPI into two inorganic phosphates (PI), which keep plant cells at a low level of cytoplasmic PPI [62]. In this study, most of the DEGs related to electron transport in the oxidative phosphorylation pathway were upregulated genes, including SDH, SAH, PPA, AHA, and VHA. NADH and succinate are dehydrogenated by NADH dehydrogenase and succinate dehydrogenase, respectively, which release protons and accelerate the efficiency of the electron transport chain, respectively, Catalyzed by PPA enzymes, PPi was hydrolyzed to Pi and ADP, while the upregulated expression of VHA and AHA accelerated the rate of electron transfer inside and outside the mitochondria and promoted the ATP synthesis [63]. We speculated that the upregulation of SDH, SAH, PPA, AHA, and VHA accelerated the electron transport rate inside and outside of the mitochondria, inhibited the generation of ROS, improved the efficiency of mitochondrial oxidative phosphorylation and the production efficiency of ATP, thus providing sufficient energy for the development of microspores, which might be the important reason for the accelerated rate of isolated microspore embryogenesis. It should be noted that carbohydrates such as sugars, as energy donors, played an important role in the growth and development of plants. In this study, after treatment with 100 mg. L− 1 IAA, the contents of soluble total sugar, amino acids, starch and soluble protein in flower buds at the late uni-nucleus stage were markedly higher than those in the control group. This might further explain that the adequate supply of energy was the important reason for the increase of the incidence of isolated microspore embryogenesis.

Degradation of the cell wall may have accelerated microspore development

Cell proliferation, expansion, and differentiation are the essential cytological basis for plant growth and development. There is a close link between this series of cytological events and cell wall metabolism. As a significant component of the primary cell wall of higher plants, the metabolism of pectin modification, synthesis, and degradation was involved in numerous aspects of plant growth and development by affecting the mobility and malleability of the cell wall [63]. Pectin methyl esterase (PME) demethylates pectin to form pectic acid and methanol, which alters the pH value of the cell wall and relax it [64]. Several PME genes are involved in pollen development in Arabidopsis and cabbage[65]. Polygalacturonase (PGs) uses polygalacturonic acid as the substrate. It acts on alpha-1, 4-D- galactosides bond to hydrolyze the pectic acid on the polygalacturonic acid main chain into galacturonic acid, galacturonic acid, the main component of pectin in the plant cell wall, is degraded into a single polygalactose residue which exists in tapetum, sporophyte, and pollen [66, 67]. In this study, almost all PME and PGs were upregulated, and in addition, the pectin content was significantly lower than that in the control group. Therefore, we speculated that upregulation of key genes involved in pectinase accelerated pectin degradation, resulting in a decrease in pectin content, degradation of the cell wall of pollen mother cells, and accelerated release of tetrad microspores. Pectate lyase (PLY) is present in anthers and pollen and can loosen cell walls [68]. USP catalyzes the reversible conversion of D-glucuronate-1P to UDP-D-glucuronate and ultimately UDP-glucose synthesis. Studies have shown that a deficiency in USP activity leads to male sterility [69,70,71]. USP is involved in the rapid growth of plant cells [72].

In this study, 75% of PLY and USP were upregulated genes, and we speculate that upregulation of PLY may lead to partial autolysis of cell wall and loosening of cell wall. The upregulation of USP may have accelerated cell growth, which may further lead to an increased incidence of isolated microspore embryogenesis.

The significant expression of DEGs may be related to microspore embryogenesis

Exploring the genes involved in the initiation of microspore embryos will help to better understand the molecular mechanism of the molecular mechanism of microspore embryogenesis. In this study, we further screened 139 DEGs from genes related to metabolic pathways and development, 97 of which were upregulated. WUS, GSTU12, PEX3, PRK3, BGL13, BFRUCT3, PME5, CYP707A3, AHA8, AHA9, PPA3, and PPA6 were highly expressed in the IAA treated group. WUS plays an important role in the embryonic development of Arabidopsis thaliana [73], and a study has shown that GST is related to the embryonic development of sorghum [74]. We speculated that these genes were related to the occurrence and development of microspore embryos, and further research was needed.


In this study, we first determined that 100 mg·L− 1 IAA was the optimal hormone concentration for inducing isolated microspore embryogenesis by counting the embryo yield. Cytological observation and the statistics of initial time of embryogenesis indicated that IAA treatment accelerated the development of microspores. Treatment with exogenous IAA caused changes in the content and dynamic balance of endogenous plant hormones. We have focused on the transcript and pathway, particularly affecting embryonic development. Transcriptome analysis, it was found that most genes related to GA and Auxin (IAA) synthesis and signal transduction, PMEs and PGs genes, and genes related to ATP synthesis and electron transport chain were upregulated, and genes related to ABA synthesis and signal transduction were downregulated. Physiological and biochemical indexes showed that exogenous IAA altered the contents of soluble total sugars, amino acids, starch, soluble protein, malondialdehyde and pectin, CAT and POD activities, and the production rates of H2O2 and O2·. These results indicated that exogenous IAA treatment could change the balance of endogenous hormones, accelerate cell wall degradation, promote ATP synthesis and nutrient accumulation, and inhibit ROS accumulation, which ultimately promotes microspore embryogenesis.

Materials and methods

Plant material and microspore culture

Line Huiwu (HW) was screened by Chenggang Wang, and all seed materials were kept in the Vegetable Breeding Laboratory of the College of Horticulture, Anhui Agricultural University (Hefei, Anhui Province, China). No specific permits were required for the described field studies. The location is not privately-owned or protected in any way, and the field studies did not involve endangered or protected species.

The Wucai varieties HW grew under 25℃ and light intensity was 20,000 lx conditions in the plastic greenhouse of the Wanjiang Vegetable Industry and Technology Research Institute (Ma’anshan, Anhui Province, China) in the spring.

At the full bloom stage (25-75% of the whole plant flowers open period), the buds at the mononuclear stage were treated with 0 mg·L− 1, 50 mg·L− 1, 75 mg·L− 1, 100 mg·L− 1, 125 mg·L− 1, 150 mg·L− 1 IAA and collected after 24 h. Bud pre-treatment was performed according to method of Zhang et al. [75]. In short, the collected buds were first sterilized with 75% alcohol for 30 s, then with 0.1% mercuric chloride for 6 min, and finally rinsed with sterile fresh water three times for 5 min each time. Microspores were isolated from the anthers by grinding the buds in B5 liquid medium using sterile glass rods and filtering the suspension through a 30 μm nylon mesh. Centrifugation was performed three times (800×g, 3 min each), and the microspores were resuspended in isolation medium. The centrifuged microspores were collected and resuspended in NLN-13 medium, with 13% sucrose and buffered at pH 5.8.The microspore density was adjusted to 1 × 105 microspores per milliliter using a hemocytometer (AP-0650010,MARIENFELD, Germany) [76]. The microspore suspension was dispensed into 60 mm × 15 mm Petri dishes with 5 mL per plate, incubated in the dark at 33℃ for 24 h and then transferred to the dark at 25 ℃ for embryonic development.

Cytological observation of microspore development

Cytological examination of the microspores for the first two weeks of culture after treatment with 100 mg·L− 1 IAA and the control group. Representative microspores cultured for 0,3,5,10, and 15 days were selected from the IAA treatment and the control group respectively, with three replicates. The microspores in 3 petri dishes were randomly selected together with the medium and inhaled into a 20 mL centrifuge tube, centrifuged at 800 rpm for 4 min, and the supernatant was discarded. After collection, the microspores were sucked onto the slide, covered with the cover slide, cut into pieces, and observed under an inverted microscope.

Determination of endogenous plant hormone content

Six hormones were selected and determined in flower buds treated with 100 mg. L− 1 IAA or the fresh water respectively by using LC-MS/MS method, including Auxin, CKs, GA, JA and ABA. 0.3 g fresh buds were collected, immediately frozen in liquid nitrogen, ground into powder (30 Hz, 1 min), and stored at -80 °C. 50 mg of plant sample was weighed into a 2 mL plastic microtube and frozen in liquid nitrogen, dissolved in 1 mL methanol/water/formic acid (15:4:1, V/V/V). 10 µL internal standard mixed solution (100 ng/mL) was added into the extract as internal standards (IS) for the quantication. The mixture was vortexed for 10 min, then centrifugation for 5 min (12,000 r/min, and 4 °C), the supernatant was transferred to clean plastic microtubes, followed by evaporation to dryness and dissolved in 100 µL 80% methanol (V/V), and filtered through a 0.22 μm membrane filter for further LC-MS/MS analysis [77,78,79].

The chromatographic column was a Waters ACQUITY UPLC HSS T3 C18, 1.8 μm, 2.1 mm × 100 mm. The mobile phase (phase A) was ultrapure water with 0.04% acetic acid, and the organic phase (phase B) was acetonitrile with 0.04% acetic acid. The elution gradient was set as follows: 0 min, A: B = 95:5 (V: V); 1 min, A: B = 95:5 (V: V); 8 min, 5:95 (V: V); 9 min, 5:95 (V: V); 9.1 min, 95:5 (V: V); and 12 min, 95:5 (V: V). The flow rate was 0.35 mL/min, the column temperature was 40 ◦C, and the injection volume was 2 µL.

Mass detection using Q-Trap 6500+ (SCIEX, USA). The mass spectrometry conditions were as follows: electrospray ionization (ESI) temperature, 550℃; mass spectrum voltage, 5500 V; curtain gas (CUR), 35 psi. Linear ion trap (LIT) and triple quadrupole (QQQ) scans were acquired on a triple quadrupole-linear ion trap mass spectrometer (QTRAP), QTRAP® 6500 + LC-MS/MS System, equipped with an ESI Turbo Ion-Spray interface, operating in both positive and negative ion mode and controlled by Analyst 1.6.3 software (Sciex). Multiquant 3.0.3 software (Sciex) was used to quantify all metabolites [80,81,82].

Transcriptome analysis and gene annotation

Total RNA from buds of treatment with 100 mg·L− 1 IAA and the control group were extracted using RNAprep Pure Plant Plus kit (Tiangen Biotech, Beijing, China) according to the manufacturer’s instructions. RNA degradation and contamination was monitored on 1% agarose, RNA purity was checked using the NanoPhotometer® spectrophotometer (IMPLEN, CA, USA). RNA concentration was measured using Qubit® RNA Assay Kit in Qubit®2.0 Flurometer (Life Technologies, CA, USA). RNA integrity was assessed using the RNA Nano 6000 Assay Kit of the Bioanalyzer 2100 system (Agilent Technologies, CA. USA). Sequencing libraries were generated using NEBNext® UltraTM RNA Library Prep Kit for Illumina® (NEB, USA) following manufacturer’s recommendations and index codes were added to attribute sequences to each sample. After the completion of library construction, the Qubit2.0 Fluorometer was used for initial quantification, followed by Agilent 2100 bioanalyzer for insert size of the library. After the insert size met the expectation, qRT-PCR accurately quantified the effective concentration of the library (the effective concentration of the library was higher than 2nM) to ensure the quality of the library. The clustering of the index-coded samples was performed on a Bot Cluster Generation System using TruSeq PE Cluster Kit v3-cBot-HS (Illumia) according to the manufacturer’s instructions. After cluster generation, the library preparations were sequenced on an Illumina platform and 150 bp paired-end. Use fastp v 0.19.3 to filter the original data, mainly to remove reads with adapters, all subsequent analyses are based on clean reads. Use HISAT v2.1.0 to compare clean reads with Brassica rapa v 3.0 reference genome. Use String Tie v1.3.4d for new gene (genes not compared to reference genome) prediction. Use featureCounts v1.6.2 /StringTie v1.3.4d to calculate the gene alignment and FPKM. DESeq2 v1.22.1 /edgeR v3.24.3 was used to analyze the differential expression between the two groups, and the P-value was corrected using the Benjamini & Hochberg method. The corrected P-value and |log2foldchange| are used as the threshold for significant difference expression. The enrichment analysis was performed based on the hypergeometric test. For KEGG, the hypergeometric distribution test is performed with the unit of pathway; for GO, it was performed based on the GO term.

Bud physiological indicators and enzyme activity determination

The buds were treated with 100 mg. L− 1 of IAA between 8: 00 a.m. and 9: 00 a.m. at the full bloom stage (25-75% flowering period of the whole plant). The control treatment was sprayed with fresh water under the same conditions. After 24 h, buds with the largest population of late-uninucleate-stage microspores stage were selected for determination of soluble sugar, soluble protein, starch, ROS related enzymes and Oxidative damage index. The experiment was repeated three times. The soluble sugar content was determined by an anthrone colorimetric method, soluble protein content was determined by Coomassie brilliant blue G-250 method [83], and starch content was determined by perchloric acid extraction-continuous flow method [84].

A centrifuge tube was used to collect 50 mg dry flower bud sample, which was mixed with 4 mL of 80% ethanol. The sample was treated with a fresh water bath at 80 ℃ for 30 min, cooled and centrifuged at 5000×g for 3 min. The supernatant was collected, and the residue was washed with 80% ethanol three times. After each washing, the supernatants were extracted and combined. Then 10 mg of activated carbon was added, decolorized at 80 ℃ for 30 min, diluted to 10 mL, and then filtered to obtain the extract.

Determination of soluble total sugar

20 µL of the above extractive solution was combined with 480 µL of fresh water and 2.5 mL of anthrone, mixed, reacted in a 90℃ fresh water bath for 15 min, taken out, cooled rapidly, and develop for 10 min before being measured at the wavelength of 620 nm.

Determination of starch content

10 mL of 52% perchloric acid was added into the above residue, the residue was centrifuged at 4000×g for 3 min, and the supernatant was collected. Then 10 mg of activated carbon was added, and the residue was treated with a fresh water bath at 80℃ for 30 min. After decolorization, filtration was performed, 5 mL of anthrone was added into 2 mL of filtrate, and the filtrate was shaken up immediately. Then the filtrate was treated with a fresh water bath at 80℃ for 10 min, cooled, developed in the dark for 10 min, and then the absorbance value was measured at a wavelength of 620 nm.

Determination of soluble protein

After washing, 0.2 g of fresh buds were placed in a pre-cooled mortar, 1.6 mL of 5 mmol. L− 1 pre-cooled phosphate buffer solution (pH 7.8) was added and ground on an ice bath to homogenate, then transferred to a centrifuge tube for centrifugation at 4℃ and 12,000×g for 20 min. To 20 uL of the extract, 80 uL of fresh water was added, followed by 2.9 mL of Coomassie brilliant blue solution. After reacting for 2 min, the absorbance was measured at 595 nm.

Determination of amino acid, CAT, POD, MDA, O2·, H2O2, and Protopectin content

Amino acid (AA), MDA content, superoxide anion (O2·), hydrogen peroxide (H2O2) content and catalase (CAT), peroxidase (POD) activities, Protopectin content were determined using Solarbio reagent kits (BC1570, BC0020, BC0290, BC3595, BC0205, BC0090, BC3685; Solarbio, Beijing, China). Each sample (fresh buds) about 50–200 mg (fresh weight) thoroughly ground for the determination of the above indicators.

RT-qPCR identification

RNA was isolated from the flower buds with single nucleus sidelined treated with 100 mg. L− 1 IAA and treated with the same dose of fresh water, and DNase-treated RNA (1 mg) was reverse transcribed into cDNA using a PrimeScript™ RT kit (TaKaRa, Japan). RT-qPCR was performed using SYBR® Premix Ex Taq™ II kit (TaKaRa, Japan). PCR amplification was performed in a Bio-Rad CFX96 instrument according to the manufacturer’s instructions. Three biological replicates were performed for each sample, and three technical replicates were performed for each gene. The data were standardized with the expression level of BnaActin as the control. The relative expression levels were calculated using the 2 −ΔΔCt method. Primer Software Version v6.0 (Premier Biosoft International, Palo Alto, CA, USA) was used to design primers, and primer sequences were listed in additional file 8: Table s4.

Data analysis

All data are expressed as the mean ± standard deviation, using at least three biological replicates. Data were statistically analyzed using spass.26 (SPSS Institute, Chicago, IL, USA) software, and means were compared using the t tests with a significance level of 0.05. Graphs were drawn using GraphPad Prism software (GraphPad Software Inc., USA), Origin 2018 64 Bit, and Adobe Illustrator CC 2019.

Data Availability

The raw RNA-Seq data used in this study have been deposited in the Nation Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA) database under the accession number PRJNA933255 ( PRJNA933255).



Indoleacetic acid


Differentially Expressed Genes


Kyoto Encyclopedia of Genes


Gene Ontology


Quantitative real-time polymerase chain reaction


RNA sequencing technology


The expected number of fragments per kilobase of transcripts per million mapped fragments


Doubled haploids




non-specific lipid transferase protein


glutathione S-transferase






Abscisic acid


Programmed cell death


Reactive Oxygen Species


Jasmonic acid


Salicylic acid




protein RALF-like


indole 3 acetamide


indole 3 acetaldehyde


aldehyde dehydrogenase




β-carotene isomerase


abscisic β-glucosyltransferase


abscisic acid 8’-hydroxylase


abscisic acid receptor


protein phosphatase


succinic dehydrogenase


V-type proton ATPase




UDP- sugar pyrophosphorylase


9-Cis-epoxycarotenoid dioxygenase


inorganic phosphates


Pectin methyl esterase




Pectate lyase


Amino acid





H2O2 :

hydrogen peroxide


  1. Zou M, Yuan L, Zhu S, et al. Effects of heat stress on photosynthetic characteristics and chloroplast ultrastructure of a heat-sensitive and heat-tolerant cultivar of Wucai (Brassica campestris L). Acta Physiol Plant. 2016;39(1):1–10.

    Google Scholar 

  2. Zhao M, Yuan L, Wang J, et al. Transcriptome analysis reveals a positive effect of brassinosteroids on the photosynthetic capacity of Wucai under low temperature[J]. BMC Genomics. 2019;20(1):1–19.

    Article  Google Scholar 

  3. M DJ. Haploids in flowering plants: origins and exploitation. Plant Biotechnol J. 2010;8(4):377–424.

    Article  Google Scholar 

  4. Testillano PS. Microspore embryogenesis: targeting the determinant factors of stress-induced cell reprogramming for crop improvement[J]. J Exp Bot. 2019;70(11):2965–78.

    Article  CAS  PubMed  Google Scholar 

  5. Chen Y, Wang Y, Xu L, et al. Effects of genotype and culture conditions on microspore embryogenesis in radish (Raphanus sativus L.) [J]. Mol Breeding. 2022;42(8):43.

    Article  Google Scholar 

  6. Weijers D, Nemhauser J, Yang Z. Auxin: small molecule, big impact[J]. J Exp Bot. 2018;69(2):133–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Rademacher EH, Lokerse AS, Schlereth A, et al. Different auxin response machineries control distinct cell fates in the early plant embryo[J]. Dev Cell. 2012;22(1):211–22.

    Article  CAS  PubMed  Google Scholar 

  8. Pérez-Pérez Y, El-Tantawy AA, Solís MT, et al. Stress-induced microspore embryogenesis requires endogenous auxin synthesis and polar transport in barley[J]. Front Plant Sci. 2019;10:1200.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Haddadi P, Moieni A, Karimzadeh G, et al. Effects of gibberellin, abscisic acid and embryo desiccation on normal plantlet regeneration, secondary embryogenesis and callogenesis in microspore culture of Brassica napus L. cv. PF704[J]. Int J Plant Prod. 2008;2(2):153–62.

    CAS  Google Scholar 

  10. Hays DB, Yeung EC, Pharis RP. The role of gibberellins in embryo axis development[J]. J Exp Bot. 2002;53(375):1747–51.

    Article  CAS  PubMed  Google Scholar 

  11. Żur I, Dubas E, Krzewska M, et al. Hormonal requirements for effective induction of microspore embryogenesis in triticale (× Triticosecale Wittm.) Anther cultures[J]. Plant Cell Rep. 2015;34:47–62.

    Article  PubMed  Google Scholar 

  12. Chen B, Fiers M, Dekkers BJW, et al. ABA signalling promotes cell totipotency in the shoot apex of germinating embryos[J]. J Exp Bot. 2021;72(18):6418–36.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Hansen M. ABA treatment and desiccation of microspore-derived embryos of cabbage (Brassica oleracea ssp. capitata L.) improves plant development[J]. J Plant Physiol. 2000;156(2):164–7.

    Article  CAS  Google Scholar 

  14. Solís M, Berenguer E, Risueño M, et al. BnPME is progressively induced after microspore reprogramming to embryogenesis, correlating with pectin de-esterification and cell differentiation in Brassica napus. BMC Plant Biol. 2016;16:1–15.

    Article  Google Scholar 

  15. Testillano P. Microspore embryogenesis: targeting the determinant factors of stress-induced cell reprogramming for crop improvement. J Exp Bot. 2019;70(11):2965–78.

    Article  CAS  PubMed  Google Scholar 

  16. Ridley BL, O’Neill MA, Mohnen D. Pectins: structure, biosynthesis, and oligogalacturonide-related signaling[J]. Phytochemistry. 2001;57(6):929–67.

    Article  CAS  PubMed  Google Scholar 

  17. Pelloux J, Rustérucci C, Mellerowicz E. New insights into pectin methylesterase structure and function. Trends Plant Sci. 2007;12:267–77.

    Article  CAS  PubMed  Google Scholar 

  18. Gómez M, Renau-Morata B, Roque E, Polaina J, Beltran J, et al. PsPMEP, a pollen-specific pectin methylesterase of pea (Pisum sativum L). Plant Reprod. 2013;26:245–54.

    Article  PubMed  Google Scholar 

  19. Pérez-Pérez Y, Bárány I, Berenguer E, et al. Modulation of autophagy and protease activities by small bioactive compounds to reduce cell death and improve stress-induced microspore embryogenesis initiation in rapeseed and barley[J]. Volume 14. Plant signaling & behavior; 2019. p. 1559577. 2.

  20. De Pinto MC, Locato V, De Gara L. Redox regulation in plant programmed cell death[J]. Plant Cell Environ. 2012;35(2):234–44.

    Article  PubMed  Google Scholar 

  21. Gechev TS, Van Breusegem F, Stone JM, et al. Reactive oxygen species as signals that modulate plant stress responses and programmed cell death[J]. BioEssays. 2006;28(11):1091–101.

    Article  CAS  PubMed  Google Scholar 

  22. Li Y, Qin T, Wei C, et al. Using transcriptome analysis to screen for key genes and pathways related to cytoplasmic male sterility in cotton (Gossypium hirsutum L.)[J]. Int J Mol Sci. 2019;20(20):5120.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Yan M, Xie D, Cao J, et al. Brassinosteroid-mediated reactive oxygen species are essential for tapetum degradation and pollen fertility in tomato[J]. Plant J. 2020;102(5):931–47.

    Article  CAS  PubMed  Google Scholar 

  24. Zhou B, Liu Y, Chen Z, et al. Comparative transcriptome analysis reveals the cause for Accumulation of reactive oxygen species during Pollen abortion in cytoplasmic male-sterile Kenaf Line 722HA[J]. Int J Mol Sci. 2019;20(21):5515.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Mittler R. ROS are good. Trends Plant Sci. 2017;22(1):11–9.

    Article  CAS  PubMed  Google Scholar 

  26. Wu G, Wei X, Wang X, et al. Changes in biochemistry and histochemical characteristics during somatic embryogenesis in Ormosia henryi Prain[J]. Volume 144. Plant Cell, Tissue and Organ Culture (PCTOC); 2021. pp. 505–17.

  27. Wang Q. Studies on microspore-derived embryogenesis and developmental pathway in non-heading chinese cabbage[J]. Journal of Nanjing Agricultural University; 2012.

  28. Dong Y, Gao Y, Zhao T, et al. Influencing factors and physiochemical changes of embryogenesis through in vitro isolated microspore culture in Brassica species[J]. Biologia. 2021;76(9):2629–54.

    Article  CAS  Google Scholar 

  29. Gu H, Zhou W, Hagberg P. High frequency spontaneous production of doubled haploid plants in microspore cultures of Brassica rapa ssp. chinensis[J] Euphytica. 2003;134:239–45.

    Article  Google Scholar 

  30. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2[J]. Genome Biol. 2014;15(12):1–21.

    Article  Google Scholar 

  31. Varet H, Brillet-Guéguen L, Coppée JY, et al. SARTools: a DESeq2-and EdgeR-based R pipeline for comprehensive differential analysis of RNA-Seq data[J]. PLoS ONE. 2016;11(6):e0157022.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Kanehisa Minoru,Furumichi Miho,Sato Yoko,Kawashima Masayuki,IshiguroWatanabe Mari. KEGG for taxonomy-based analysis of pathways and genomes [J]. Nucleic acids research,2022.

  33. Bayer M, Nawy T, Giglione C, et al. Paternal control of embryonic patterning in Arabidopsis thaliana[J]. Science. 2009;323(5920):1485–8.

    Article  CAS  PubMed  Google Scholar 

  34. Yan L, Tan Z, Zhu L, et al. MYB97, MYB101 and MYB120 function as male factors that control pollen tube-synergid interaction in Arabidopsis thaliana fertilization[J]. PLoS Genet. 2013;9(11):e1003933.

    Article  Google Scholar 

  35. Yao X, Tian L, Yang J, et al. Auxin production in diploid microsporocytes is necessary and sufficient for early stages of pollen development[J]. PLoS Genet. 2018;14(5):e1007397.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Song S, Liu J, Tang C, et al. Metabolism and signaling of auxins and their roles in regulating seed dormancy and germination[J]. Chin Sci Bull. 2020;65(34):3924–43.

    Article  Google Scholar 

  37. Bender J, Celenza JL. Indolic glucosinolates at the crossroads of tryptophan metabolism[J]. Phytochem Rev. 2009;8:25–37.

    Article  CAS  Google Scholar 

  38. Kasahara H. Current aspects of auxin biosynthesis in plants[J]. Bioscience, Biotechnology, and Biochemistry, 2016, 80(1): 34–42.

  39. Sugawara S, Hishiyama S, Jikumaru Y, et al. Biochemical analyses of indole-3-acetaldoxime-dependent auxin biosynthesis in Arabidopsis. Proc Natl Acad Sci. 2009;106(13):5430–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Skibbe D, Liu F, Wen T, et al. Characterization of the aldehyde dehydrogenase gene families of Zea mays and Arabidopsis[J]. Plant Mol Biol. 2002;48:751–64.

    Article  CAS  PubMed  Google Scholar 

  41. Nafisi M, Goregaoker S, Botanga C, et al. Arabidopsis cytochrome P450 monooxygenase 71A13 catalyzes the conversion of indole-3-acetaldoxime in camalexin synthesis[J]. Plant Cell. 2007;19(6):2039–52.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Xu Z, Escamilla-Treviño L, Zeng L, et al. Functional genomic analysis of Arabidopsis thaliana glycoside hydrolase family 1[J]. Plant Mol Biol. 2004;55:343–67.

    Article  CAS  PubMed  Google Scholar 

  43. Shi G, Huang F, Gong Y, et al. RNA-Seq analysis reveals that multiple phytohormone biosynthesis and signal transduction pathways are reprogrammed in curled-cotyledons mutant of soybean [Glycine max (L.) Merr.][J]. BMC Genomics. 2014;15(1):1–13.

    Article  Google Scholar 

  44. Zhu L, Jiang B, Zhu J, et al. Auxin promotes fiber elongation by enhancing gibberellic acid biosynthesis in cotton[J]. Plant Biotechnol J. 2022;20(3):423.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Zhao J, Yao X. Mechanism of Microspore Embryogenesis in Plant[J]. Biotechnol Bull, 2009(06): 7–11.

  46. Gao Y, Qu G, Huang S, et al. Comparison between germinated seed and isolated microspore EMS mutagenesis in chinese cabbage (Brassica rapa L. ssp. pekinensis)[J]. Horticulturae. 2022;8(3):232.

    Article  Google Scholar 

  47. Kim J, DellaPenna D. Defining the primary route for lutein synthesis in plants: the role of Arabidopsis carotenoid β-ring hydroxylase CYP97A3. Proc Natl Acad Sci. 2006;103(9):3474–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Alder A, Jamil M, Marzorati M, et al. The path from β-carotene to carlactone, a strigolactone-like plant hormone[J]. Science. 2012;335(6074):1348–51.

    Article  CAS  PubMed  Google Scholar 

  49. Tan B, Schwartz S, Zeevaart J et al. Genetic control of abscisic acid biosynthesis in maize[J]. Proceedings of the National Academy of Sciences, 1997, 94(22): 12235–12240.

  50. Okamoto M, Kuwahara A, Seo M, et al. CYP707A1 and CYP707A2, which encode abscisic acid 8′-hydroxylases, are indispensable for proper control of seed dormancy and germination in Arabidopsis[J]. Plant Physiol. 2006;141(1):97–107.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Zhang Z, Zhang H. Relationship between fertility transition of thermo-photo-sensitive wheat male sterile line BNS and endogenous hormone contents in its developing ear[J]. J Plant Genetic Resour. 2016;17:913–9.

    CAS  Google Scholar 

  52. Jiang H, Sun Y, Xu Y, et al. Physiological and biochemical changes during bud development in male sterile plant of Camellia crassicolumna[J]. Plant Physiol J. 2020;56(9):1807–17.

    Google Scholar 

  53. Hu L, Liang W, Yin C, et al. Rice MADS3 regulates ROS homeostasis during late anther development[J]. Plant Cell. 2011;23(2):515–33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Møller IM. Plant mitochondria and oxidative stress: electron transport, NADPH turnover, and metabolism of reactive oxygen species[J]. Annu Rev Plant Biol. 2001;52(1):561–91.

    Article  Google Scholar 

  55. Cui KR, Xing GS, Liu XM, et al. Effect of hydrogen peroxide on somatic embryogenesis of Lycium barbarum L[J]. Plant Sci. 1999;146(1):9–16.

    Article  CAS  Google Scholar 

  56. Wang Y. Studies on the physiological, biochemical and molecular mechanisms of the sterility characteristics of fresh watermelon genic male sterile line Se18. D. Northwest A & F University; 2020.

  57. Liu Z, Yang Y, Hu T, et al. Antioxidant and osmoregulation substance properties of flower buds in Plant grandiflorus male sterile lines[J]. Acta Prataculturae Sinica. 2012;21(02):77–82.

    Google Scholar 

  58. Serrano-Bueno G, Hernández A, López-Lluch G, et al. Inorganic pyrophosphatase defects lead to cell cycle arrest and autophagic cell death through NAD + depletion in fermenting yeast[J]. J Biol Chem. 2013;288(18):13082–9267.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Aggarwal RK, Luchtel RA, Machha V, et al. Functional succinate dehydrogenase deficiency is a common adverse feature of clear cell renal cancer. Proc Natl Acad Sci. 2021;118(39):e2106947118.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Stitt M. Pyrophosphate as an Energy Donor in the cytosol of Plant cells: an enigmatic alternative to ATP. Bot Acta. 1998;111(3):167–75.

    Article  CAS  Google Scholar 

  61. Serrano-Bueno G, Hernández A, López-Lluch G, et al. Inorganic pyrophosphatase defects lead to cell cycle arrest and autophagic cell death through NAD + depletion in fermenting yeast[J]. J Biol Chem. 2013;288(18):13082–92.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Hoffmann RD, Portes MT, Olsen LI, et al. Plasma membrane H+-ATPases sustain pollen tube growth and fertilization[J]. Nat Commun. 2020;11(1):2395.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Ren C, Kermode AR. An increase in pectin methyl esterase activity accompanies dormancy breakage and germination of yellow cedar seeds[J]. Plant Physiol. 2000;124(1):231–42.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  64. Rhee SY, Somerville CR. Tetrad pollen formation in quartet mutants of Arabidopsis thaliana is associated with persistence of pectic polysaccharides of the pollen mother cell wall[J]. Plant J. 1998;15(1):79–88.

    Article  CAS  PubMed  Google Scholar 

  65. Xiong X, Zhou D, Xu L, et al. BcPME37c is involved in pollen intine formation in Brassica campestris[J]. Biochem Biophys Res Commun. 2019;517(1):63–8.

    Article  CAS  PubMed  Google Scholar 

  66. Dubald M, Barakate A, Mandaron P, et al. The ubiquitous presence of exopolygalacturonase in maize suggests a fundamental cellular function for this enzyme[J]. Plant J. 1993;4(5):781–91.

    Article  CAS  PubMed  Google Scholar 

  67. Barakate A, Martin W, Quigley F, et al. Characterization of a multigene family encoding an exopolygalacturonase in maize[J]. J Mol Biol. 1993;229(3):797–801.

    Article  CAS  PubMed  Google Scholar 

  68. Marín-Rodríguez MC, Orchard J, Seymour GB. Pectate lyases, cell wall degradation and fruit softening[J]. J Exp Bot. 2002;53(377):2115–9.

    Article  PubMed  Google Scholar 

  69. Chen R, Zhao X, Shao Z, et al. Rice UDP-glucose pyrophosphorylase1 is essential for pollen callose deposition and its cosuppression results in a new type of thermosensitive genic male sterility[J]. Plant Cell. 2007;19(3):847–61.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  70. Mu H, Ke J, Liu W, et al. UDP-glucose pyrophosphorylase2 (OsUgp2), a pollen-preferential gene in rice, plays a critical role in starch accumulation during pollen maturation[J]. Chin Sci Bull. 2009;54(2):234–43.

    CAS  Google Scholar 

  71. Park J, Ishimizu T, Suwabe K, et al. UDP-glucose pyrophosphorylase is rate limiting in vegetative and reproductive phases in Arabidopsis thaliana[J]. Plant Cell Physiol. 2010;51(6):981–96.

    Article  CAS  PubMed  Google Scholar 

  72. Kotake T, Hojo S, Yamaguchi D, et al. Properties and physiological functions of UDP-sugar pyrophosphorylase in Arabidopsis[J]. Biosci Biotechnol Biochem. 2007;71(3):761–71.

    Article  CAS  PubMed  Google Scholar 

  73. Wójcikowska B, Gaj MD. Somatic embryogenesis in Arabidopsis[J]. Somatic Embryogenesis: Fundamental Aspects and Applications, 2016: 185–99.

  74. Liu H, Micic N, Miller S, et al. Species-specific dynamics of specialized metabolism in germinating sorghum grain revealed by temporal and tissue-resolved transcriptomics and metabolomics[J]. Plant Physiol Biochem. 2023;196:807–20.

    Article  CAS  PubMed  Google Scholar 

  75. Zhang K, Wang Y, Wang Q. Study on the key factors influencing isolated Microspore Culture in Pakchoi[J]. Seed. 2022;41(05):144–8.

    Google Scholar 

  76. Camacho-Fernández C, Hervás D, Rivas-Sendra A, et al. Comparison of six different methods to calculate cell densities[J]. Plant Methods. 2018;14:1–15.

    Article  Google Scholar 

  77. Sheng Y, Chang W, Jiao S, et al. Stamen structure development and physiological and biochemical characteristics in male sterile melon[J]. Plant Physiol J. 2016;52(7):1028–34.

    Google Scholar 

  78. Li Y, Zhou C, Yan X, et al. Simultaneous analysis of ten phytohormones in Sargassum horneri by high-performance liquid chromatography with electrospray ionization tandem mass spectrometry[J]. J Sep Sci. 2016;39(10):1804–13.

    Article  PubMed  Google Scholar 

  79. Floková K, Tarkowská D, Miersch O, et al. UHPLC–MS/MS based target profiling of stress-induced phytohormones[J]. Phytochemistry. 2014;105:147–57.

    Article  PubMed  Google Scholar 

  80. Pan X, Welti R, Wang X. Quantitative analysis of major plant hormones in crude plant extracts by high-performance liquid chromatography–mass spectrometry[J]. Nat Protoc. 2010;5(6):986–92.

    Article  CAS  PubMed  Google Scholar 

  81. Cui K, Lin Y, Zhou X, et al. Comparison of sample pretreatment methods for the determination of multiple phytohormones in plant samples by liquid chromatography–electrospray ionization-tandem mass spectrometry[J]. Microchem J. 2015;121:25–31.

    Article  CAS  Google Scholar 

  82. Šimura J, Antoniadi I, Široká J, et al. Plant hormonomics: multiple phytohormone profiling by targeted metabolomics[J]. Plant Physiol. 2018;177(2):476–89.

    Article  PubMed  PubMed Central  Google Scholar 

  83. Zhu LC, Tang SJ, Zhou L. Teaching practice and methodological investigation of protein content determination using coomassie brilliant blue G250[C]//Educ. Teach Forum. 2020;23:266–9.

    Google Scholar 

  84. Li R, Guo P, Ning Z, et al. Comparison of detection methods for starch in tobacco[J]. Mod Food Sci Technol. 2012;28(02):229–32.

    Google Scholar 

Download references


We thank SZ for his work and would like to thank the reviewers for their valuable comments and suggestions. We thank LetPub (www.letpub. com) for its linguistic assistance during the preparation of this manuscript.


This work was supported by the Key Project Foundation of Natural Science Research in Universities of Anhui Province (KJ2019A0189), Scientific Research Project of Anhui Education Department (GXXT-2022-072).

Author information

Authors and Affiliations



CW and PZ conceived and designed the experiments. XT assisted in transcriptome data analysis and interpretation. PZ and YH performed molecular, physiological and biochemical analysis. FH, XW and HL critically analyzed the data and contributed to the writing of the manuscript. YH made a feasible plan for the revision of the paper. LY, HJ and GC directed the study, including experimental design and manuscript revision. WW and JW supervised the study. All authors have read and approved the final manuscript.

Corresponding author

Correspondence to Xiaoyan Tang.

Ethics declarations

Ethics approval and consent to participate

Line HW was identified by the corresponding author, and all seed materials were kept in the Vegetable Breeding Laboratory of the College of Horticulture, Anhui Agricultural University. No specific permits were required for the described field studies. The location is not privately-owned or protected in any way, and the field studies did not involve endangered or protected species. We complied with the IUCN Policy Statement on Research Involving Species at Risk of Extinction and the Convention on the Trade in Endangered Species of Wild Fauna and Flora.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wang, C., Zhang, P., He, Y. et al. Exogenous spraying of IAA improved the efficiency of microspore embryogenesis in Wucai (Brassica campestris L.) by affecting the balance of endogenous hormones, energy metabolism, and cell wall degradation. BMC Genomics 24, 380 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: