- Open Access
Gene expression programs during callus development in tissue culture of two Eucalyptus species
BMC Plant Biology volume 22, Article number: 1 (2022)
Eucalyptus is a highly diverse genus of the Myrtaceae family and widely planted in the world for timber and pulp production. Tissue culture induced callus has become a common tool for Eucalyptus breeding, however, our knowledge about the genes related to the callus maturation and shoot regeneration is still poor.
We set up an experiment to monitor the callus induction and callus development of two Eucalyptus species - E. camaldulensis (high embryogenic potential) and E. grandis x urophylla (low embryogenic potential). Then, we performed transcriptome sequencing for primary callus, mature callus, shoot regeneration stage callus and senescence callus. We identified 707 upregulated and 694 downregulated genes during the maturation process of the two Eucalyptus species and most of them were involved in the signaling pathways like plant hormone and MAPK. Next, we identified 135 and 142 genes that might play important roles during the callus development of E. camaldulensis and E. grandis x urophylla, respectively. Further, we found 15 DEGs shared by these two Eucalyptus species during the callus development, including Eucgr.D00640 (stem-specific protein TSJT1), Eucgr.B00171 (BTB/POZ and TAZ domain-containing protein 1), Eucgr.C00948 (zinc finger CCCH domain-containing protein 20), Eucgr.K01667 (stomatal closure-related actinbinding protein 3), Eucgr.C00663 (glutaredoxin-C10) and Eucgr.C00419 (UPF0481 protein At3g47200). Interestingly, the expression patterns of these genes displayed “N” shape in the samples. Further, we found 51 genes that were dysregulated during the callus development of E. camaldulensis but without changes in E. grandis x urophylla, such as Eucgr.B02127 (GRF1-interacting factor 1), Eucgr.C00947 (transcription factor MYB36), Eucgr.B02752 (laccase-7), Eucgr.B03985 (transcription factor MYB108), Eucgr.D00536 (GDSL esterase/lipase At5g45920) and Eucgr.B02347 (scarecrow-like protein 34). These 51 genes might be associated with the high propagation ability of Eucalyptus and 22 might be induced after the dedifferentiation. Last, we performed WGCNA to identify the co-expressed genes during the callus development of Eucalyptus and qRT-PCR experiment to validate the gene expression patterns.
This is the first time to globally study the gene profiles during the callus development of Eucalyptus. The results will improve our understanding of gene regulation and molecular mechanisms in the callus maturation and shoot regeneration.
The regeneration of plant tissue or organs under culture conditions has been extensively used for decades in plants. Plants generate callus in response to stresses like wounding and pathogen infection, and callus cells are totipotent and are able to regenerate the whole plant body . Depending on the organs they generate, calli with some degrees of organ regeneration are called with different names, such as rooty, shooty, embryogenic callus, and compact callus [2, 3]. In vitro, exogenous application of auxin and cytokine has been proved to induce callus in plants. For example, Skoog and Miller showed that a high ratio of auxin-to-cytokinin can induce root regeneration, while a high ratio of cytokinin-to-auxin can induce shoot regeneration . In nature, wounding, pathogens and interspecific hybrids are common ways to induce callus and tumors . Some molecules have been identified to play key roles during the callus induction and development. For example, lateral organ boundaries domain (LBD) family of transcription factors (e.g., LBD16, LBD17, LBD18, LBD29) can mediate the expression of auxin response factors ARF7 and ARF19 [5, 6]. ARR1 and ARR21 have been identified to induce callus formation in Arabidopsis [7, 8]. RWP-RK domain transcription factors like RKD1, RKD2 and RKD4 have been found to mediate the gametogenesis and embryogenesis [9, 10]. In addition, some studies have been demonstrated to uncover the genes and proteins involved in the callus development. Tan et al. identified 73 proteins significantly differentially expressed during the callus development in Vanilla planifolia Andrews . Che et al. identified RAP2.6L as a key factor for shoot regeneration in Arabidopsis because the T-DNA knockdown mutations in RAP2.6L reduced the expression of many genes that are normally up-regulated during shoot development . However, our knowledge about the genes involved in the callus development and tissue regeneration process in plants is still poor.
Eucalyptus, a highly diverse genus of the Myrtaceae family, is widely planted in the world due to its significant economic values for timber and pulp . As we know, all commonly recognized methods of vegetative propagation have been applied with Eucalyptus, but most have resulted in failure especially when applied to adult tissues . Successful regeneration of plants or organs from selected Eucalyptus has never been reported until 1981 when callus was induced on embryos and sterile seedlings of selected trees of E. leichow . Then, this technique was applied with many Eucalyptus species like E. polybractea, E. cama1dulensis, E. gomphocephala and E. viminalis . Another technique that has been successfully applied with Eucalyptus is organ culture, in which differentiated tissues such as leaves, stems and roots are placed in a controlled system of nutrients and environment. Roots and/or buds can be induced on the explant either directly or after the formation of a callus and many tree species, including Eucalyptus, which can be propagated by organ culture techniques have applied with this method. Although organ culture techniques are often used in preference to the traditional methods of vegetative propagation due to its high multiplication rates, there are some problems usually happened in developing the organ culture techniques, such as obtaining aseptic tissue from field-grown plants, brown exudate, rooting and bud inhibitors . Also, some factors have been reported to affect the root initiation in nodes of Eucalyptus trees, such as the culture medium, the incubation conditions, the physiological state of the parent plant, and the position on the parent plan t[1, 14].
Some vegetative propagation associated studies have been demonstrated in Eucalyptus. Grattapaglia et al. identified some QTLs controlling the ability to form shoots of E. grandi s. Marques et al. identified QTLs related to adventitious rooting, sprouting ability and the stability of adventitious rooting . In plants, some genes have been reported to play key roles during the vegetative propagation, such as ARF19, SERK, LEC and WUS [17, 18]. Previously, our lab reported the transcriptome profiles of two Eucalyptus species during somatic embryogenesis and dedifferentiation . We identified genes encoding somatic embryogenesis receptor kinase, ethylene, auxin, ribosomal protein, zinc finger protein, heat shock protein, histone, cell wall related protein and multiple transcription factors that might control the ability of somatic embryogenesis and dedifferentiation. However, large is unknown about the gene regulations during the developmental process after the callus is induced in Eucalyptus.
In the present study, we aimed to investigate the transcriptome profiles of dedifferentiated callus tissues incubated on the culture medium. We also aimed to identify genes involved in the callus development process and controlling the ability of vegetative propagation. This is the first time to study the gene profiles of dedifferentiated callus tissues of Eucalyptus and our results will provide new insights of understanding the molecular mechanisms in the callus development and differentiation processes. More importantly, our results will improve our knowledge about the genes associated with the vegetative propagation ability of Eucalyptus.
Callus induction and incubation
To understand the gene expression profiles during the callus development of Eucalyptus, we obtained the stem tissues of two Eucalyptus species – E. camaldulensis (high embryogenic potential, A1) and E. grandis x urophylla (low embryogenic potential, B1). We performed the in vitro tissue-culture experiments on these stem samples and obtained callus tissues from different developmental stages (Figure 1A). Initially, we observed that the incubation time on CIM (callus-inducing medium) of the stems has a great impact on their regeneration ability. As shown in Figure 1B, the regeneration rates of the tissue culture induced callus by stem peaked at 21 days of incubation on CIM. It is notable that the regeneration rate of E. camaldulensis callus was much higher than E. grandis x urophylla. We obtained the callus (also called primary callus, pri-callus) at 10 days (A2 and B2) and mature callus at 21 days (A3 and B3). Then, the mature callus tissues were transferred to SIM (shooting-inducing medium) for further incubation. They were shown to start generating buds after 7 days incubation and 80% of the callus generated buds after 10 days incubation. We obtained the tissues of callus tissues incubated on SIM for 10 days for the two Eucalyptus species (A4 and B4). It is interesting that the callus tissues turned brown intensively, the bead-like protrusions also turned brown to black, and the callus lost the regeneration ability after incubation on CIM for 30 days. The callus tissues incubated on CIM for 35 days were obtained for the two species of Eucalyptus (A5 and B5).
Transcriptome sequencing and gene expression profiles
We employed the transcriptome sequencing for the callus tissues of E. camaldulensis (A2~A5) and E. grandis x urophylla (B2~B5). After data cleaning, we obtained 20.99 to 23.58 million reads for these samples and found 69.27% to 84.58% of the clean reads mapped to the reference Eucalyptus genome. Next, we used StringTie to identify genes expressed in the callus tissues of the two Eucalyptus species. After the average TPM (transcripts per million reads) values of all genes were calculated and lowly expressed genes (TPM < 5) were filtered, we identified 12,229 to 14,075 genes for all the samples. It showed in Figure 1C that 14,579 genes were identified in the callus tissues of both E. camaldulensis and E. grandis x urophylla and that 1,706 and 1,516 genes were expressed specifically in E. camaldulensis and E. grandis x urophylla, respectively. Then, we compared the genes identified in different stages of callus tissues. Figure 1D showed that 10,254 and 10,316 genes were commonly identified in all the callus tissues of E. camaldulensis and E. grandis x urophylla, respectively. Next, we analyzed the sample correlation during the callus development using the gene expression profiles. As expected, the replicates were performed well, and the samples can be distinguished from each other based on the gene expression profiles (Figure 1E). Further, we analyzed the gene expression profiles across species. Notably, the callus tissues showed similarities between E. camaldulensis and E. grandis x urophylla before mature callus developmental stage (Figure 1F), which indicates that the callus differentiation process varies in these two Eucalyptus species on molecular level. Based on the developmental stages, we divided the whole process into three parts to investigate the gene changes during the differentiation process, including pri-callus to mature callus (mat-callus), mature callus to shoot regeneration stage callus (SRS-callus), and mature callus to senescence callus (sen-callus).
DEGs in callus maturation
We compared the gene expression profiles of primary and mature callus tissues in the two Eucalyptus species. Initially, we identified 3,790 (1,834 upregulated and 1,956 downregulated) and 3,740 (1,834 upregulated and 1,956 downregulated) DEGs in E. camaldulensis (A3 compared to A2) and E. grandis x urophylla (B3 compared to B2), respectively (Figure 2A, additional file 1). In this process, the two Eucalyptus species shared 707 upregulated and 694 downregulated genes (Figure 2A). Notably, 106 genes were found with adverse regulation during the callus maturation process in the two Eucalyptus species (Figure 2A), including Eucgr.I01667 (kelch repeat-containing protein At3g27220), Eucgr.B00093 (HVA22-like protein e), Eucgr.G02764 (glutathione S-transferase DHAR2), Eucgr.F00590 (snakin-2), Eucgr.D00272 (CBL-interacting serine/threonine-protein kinase 23) and Eucgr.F00184 (low-temperature-induced cysteine proteinase). We next analyzed the signal transduction pathways involved by the DEGs. Notably, 127 and 125 DEGs were enriched in the plant hormone signaling transduction pathway of E. camaldulensis and E. grandis x urophylla, respectively (Figure 2B). In addition, we found 77 and 75 DEGs enriched in the plant MAPK signaling pathway of E. camaldulensis and E. grandis x urophylla, respectively (Figure 2B). Gene Ontology enrichment analysis identified that 165 and 136 DEGs were involved in the protein phosphorylation (GO:0006468) during the callus maturation process of E. camaldulensis and E. grandis x urophylla, respectively. Further, we compared the dysregulated genes during the process of stem to mature callus. It showed in Table 1 that 40 upregulated and 34 downregulated genes were identified in this process, and they might be involved in the dedifferentiation and callus development.
DEGs in callus development
We next identified DEGs in the differentiation process of callus in the two Eucalyptus species. Compared to A3, we identified 3,111 (1,793 upregulated and 1,318 downregulated) and 4,377 (2,300 upregulated and 2,077 downregulated) genes differentially expressed in A4 and A5, respectively (Figure 2C, additional file 1). As A4 has the propagation ability while A5 lost it, we next investigated the genes with diverse regulations in A4 and A5. It showed that a total of 135 genes with such characteristics (additional file 2). These genes were predicted to be enriched in the biological processes like “GO:0048827~phyllome development” (1 gene, Eucgr.B02127), “GO:0055114~oxidation-reduction process” (14 genes) and “GO:0009791~post-embryonic development” (1 gene, Eucgr.I02367). Likewise, compared to B3 we identified 4,295 (2,390 upregulated and 1,905 downregulated) and 2,358 (1,072 upregulated and 1,286 downregulated) genes differentially expressed in B4 and B5, respectively (Figure 2D). There were 142 genes with diverse regulations in B4 and B5 (additional file 2). These genes were enriched in the biological processes including “GO:0006542~glutamine biosynthetic process” (1 gene, Eucgr.G02570), “GO:0031145~anaphase-promoting complex-dependent catabolic process” (1 gene, Eucgr.J00733), “GO:0044765~single-organism transport” (1 gene, Eucgr.A00992), “GO:0006952~defense response” (1 gene, Eucgr.F03332), “GO:0016569~covalent chromatin modification” (1 gene, Eucgr.J03029), “GO:0008643~carbohydrate transport” (1 gene, Eucgr.C02790) and “GO:0006012~galactose metabolic process” (1 gene, Eucgr.C02197).
Further, we compared the DEGs with diverse regulations in SRS-callus and sen-callus, relative to mat-callus, in the two Eucalyptus species. We found 15 genes with such regulations in these samples (Table 2), including Eucgr.D00640 (stem-specific protein TSJT1), Eucgr.B00171 (BTB/POZ and TAZ domain-containing protein 1), Eucgr.C00948 (zinc finger CCCH domain-containing protein 20), Eucgr.K01667 (stomatal closure-related actin-binding protein 3), Eucgr.C00663 (glutaredoxin-C10) and Eucgr.C00419 (UPF0481 protein At3g47200). We showed the expression levels of six of these genes in Figure 2E. It is interesting that the expression patterns of all these genes showed an “N” shape during the callus development of both Eucalyptus species. For example, the expression level of Eucgr.D00640 decreased from pri-callus to mat-callus, went up from mat-callus to SRS-callus, and then went down in sen-callus.
Another group of genes involved in the callus development of Eucalyptus include the DEGs during the process from pri-callus to mat-callus to SRS-callus. In E. camaldulensis, a Venn diagram of DEGs identified in this process revealed that A4 vs A3 shared 14 upregulated and 146 downregulated genes with A3 vs A2 (Figure 2F, additional file 1). In E. grandis x urophylla, we identified 50 upregulated and 145 downregulated genes shared by B3 vs B2 and B4 vs B3 (Figure 2F, additional file 1). Next, we compared the shared DEGs in the two Eucalyptus species during the callus development process. It showed 4 upregulated and 41 downregulated genes shared by the two Eucalyptus species in this process (additional file 1), including Eucgr.B00168 (probable pectate lyase 18), Eucgr.E01615 (putative expansin-B2), Eucgr.K03562 (transcription factor MYB108), Eucgr.H03379 (GEM-like protein 5) and Eucgr.C03297 (ethylene-responsive transcription factor ERF017). Notably, Eucgr.B00168 and Eucgr.E01615 were downregulated in A5 vs A3 and B5 vs B3, respectively (additional file 1); no downregulated genes in the process of pri-callus to SRS-callus were found to be upregulated in sen-callus in both Eucalyptus species; and 17 downregulated genes were also downregulated in the sen-callus in both Eucalyptus species (additional file 1), including Eucgr.K03562, Eucgr.H03379 and Eucgr.C03297.
Key genes associated with high regeneration ability
We next investigated the genes related to high propagation ability of Eucalyptus. Due to the missing annotation of some genes in Eucalyptus genome, we first examined the expression profiles of previously reported vegetative propagation ability related genes, including ARF19, SERK, LEC and WUS. As shown in additional file 3, two genes (Eucgr.C02178, Eucgr.C03293) were found to encode ARF19 and only Eucgr.C02178 was downregulated between SRS-callus and mat-callus in both Eucalyptus species. We also found two genes (Eucgr.B01219, Eucgr.K03421) encoding LEC14B and none of them were changed in the callus tissues of E. grandis x urophylla(additional file 3). Among the three genes (Eucgr.H03383, Eucgr.I01078, Eucgr.F04151) encoding SERK, we found that Eucgr.I01078 was the dominant one. Eucgr.I01078 was not changed in the callus tissues of E. camaldulensis and peaked in B3 (additional file 3). It is interesting that Eucgr.F02320 encoding WUSCHEL-related homeobox 4 was upregulated in A3 vs A2, downregulated in B3 vs B2, and upregulated in B4/B5 compared to B3 (additional file 3).
Then, we analysed the 116 genes upregulated in A4 (compared to A3) and downregulated in A5 (compared to A3), and found 43 genes (Table 3) downregulated or with no changed in B4 (compared to B3), including Eucgr.F00590 (snakin-2), Eucgr.F02674 (putative laccase-9), Eucgr.A02259 (two-component response regulator ARR9), Eucgr.G01769(auxin transporter-like protein 5), Eucgr.B02127 (GRF1-interacting factor 1) and Eucgr.C00947 (transcription factor MYB36). The second group may contain genes downregulated in A4 (compared to A3), upregulated in A5 (compared to A3), but upregulated or not changed in B4 (compared to B3). Using these filters, we identified 8 genes (Table 3), such as Eucgr.B02752 (laccase-7), Eucgr.B03985 (transcription factor MYB108), Eucgr.D00536 (GDSL esterase/lipase At5g45920) and Eucgr.B02347 (scarecrow-like protein 34). A heat map (Figure 3A) showed the expression patterns of these 51 genes in the callus tissues of the two Eucalyptus species. It is interesting that there were 12 genes highly expressed in B3 and B4 but not changed between them, and they were upregulated in A4 compared to A3, such as Eucgr.C02990 (zinc finger CCCH domain-containing protein 2), Eucgr.A01269, Eucgr.B03374 (two-component response regulator ARR6) and Eucgr.B02127 (GRF1-interacting factor 1).
We next performed the weighted genes co-expression network analysis to identify co-expressed genes during the callus development process of the two Eucalyptus species. As shown in Figure 3B (left panel), the grey (8 genes) and blue (1,539 genes), yellow (820 genes), brown (961 genes), and turquoise (2,630 genes) modules of genes were identified to be correlated with A2, A3, A4 and A5, respectively. Notably, the yellow module of genes (co-expressed in A3) contained some transcription factor genes, such as Eucgr.C02208 (transcription factor bHLH35), Eucgr.I00291 (ethylene-responsive transcription factor ABR1), Eucgr.K01542 (transcription factor MYB44) and Eucgr.C01943 (probable WRKY transcription factor 40). Next, we showed the co-expressed genes in the callus development process of E. grandis x urophylla (Figure 3B, right panel). It is notable that there were 891 genes from the yellow module co-expressed in B3. Then, we compared the co-expressed genes at the same developmental stages of the two Eucalyptus species. It showed that 421, 166, 516 and 482 genes were co-expressed in pri-callus, mat-callus, SRS-callus and sen-callus of the two Eucalyptus species, respectively (Fig. 3C).
We next selected 12 genes and performed qRT-PCR to validate their expression patterns in the callus development of the two Eucalyptus species. The H2B gene was used as the internal control gene. The primer sequences of them can be accessed in additional file 4. We performed three reactions for each gene in one biological sample and a total of 9 reactions were used for one gene (n=9). For the qRT-PCR experiment we used log2RNE values to present the gene changes in the comparisons (mat-callus used as the control), and for transcriptome we used log2FC to show the gene changes. For the comparison of mat-callus and pri-callus (A3 vs A2 and B3 vs B2), we used -log2RNE to show the expression changes. Thus, in total there were 72 events (12 × 3 × 2) that to be validated and we have 61 (84.72%) events were agreed by both qRT-PCR and deep sequencing (Figure 4). It is notable that the dysregulation of some genes was confirmed by both experiments in E. camaldulensis and E. grandis x urophylla, such as Eucgr.B03816 (transcription factor LHW), Eucgr.C00948 (zinc finger CCCH domain-containing protein 20), Eucgr.C03301 (protein TIFY 10a), Eucgr.D00640 (stem-specific protein TSJT1) and Eucgr.J00388 (E3 ubiquitin-protein ligase MIEL1). High agreement of gene expression patterns in transcriptome sequencing and qRT-PCR indicate that the genes we found in this study might be associated with the callus development and differentiation of Eucalyptus. Their functions require more experiments to be explored.
In this study we analysed the transcriptome profiles during the callus maturation and development processes in two Eucalyptus species with different vegetative propagation capacity. The utilization of successful (SRS-callus) and failed (sen-callus) vegetative propagation callus tissues as contrast strongly support that the genes identified in this study might play important roles in the callus development and might be associated with the vegetative propagation ability in Eucalyptus.
In Eucalyptus, three major stages of regeneration from tissue culture have been reported – co-cultivation, callus induction and shoot regeneration . With our observation, after the dedifferentiation of stem tissue, pri-callus tissues still have low regeneration ability and require further incubation on CIM (Figure 1A and B), which agrees with the callus culture in maize . Not many studies have focused on the callus maturation process in plant and very little is known about the molecular changes in this process. Prior to transfer to the MS medium, fine chopping and partial desiccation of embryogenic calli can simulate the rapid maturation of somatic embryos in date palm . In Zea mays, Maturation of somatic embryos was enhanced by transferring the embryogenic callus after 3 weeks to medium containing 6% sucrose and lacking2,4-D . Many genes have been reported to be dysregulated during the embryo maturation in maize, such as histone and ribosomal protein genes, and genes encoding hydrolytic enzymes (nucleases, glucosidases and proteases) and a few storage genes (an α-zein and caleosin) . In this study, we also identified many dysregulated genes during the callus maturation (Figure 2A, additional file 1), including genes encoding histone, ribosomal proteins, nuclease, glucosidases and proteases. Notably, some genes were found to be continually up or down regulated from stem tissue to mature callus (Table 1), such as Eucgr.A02688 (ABC transporter G family member 4), Eucgr.B02604 (G-type lectin S-receptor-like serine/threonine-protein kinase LECRK2), Eucgr.E01615 (putative expansin-B2), Eucgr.A00514 (auxin transporter-like protein 4), Eucgr.C03807 (transcription factor bHLH68) and Eucgr.F01151 (early nodulin-like protein 2). ABCG14, a homologue of ABCG4, has been proved to be a positive regulator of plant growth and play an important role in the major root-to-shoot (acropetal) long-distance cytokinin (CK) transport via the xylem sap [24,25,26,27]. Another ABCG4 homologue ABCB4 is an auxin influx transporter which mediates the transport of auxin in roots and contributes to the basipetal transport in hypocotyls and root tips by establishing an auxin uptake sink in the root cap [28, 29]. Eucgr.E01615, which encodes the EXPB2 protein, was upregulated during the callus maturation process (additional file 1). Interestingly, four beta-expansin genes were found to be induced by treatment with gibberellin and by wounding in rice and correlated with rapid elongation of deep-water rice internodes . Further, EXPB2 was found to be a root-predominant gene and play a key role in the root-hair formation in rice . The expression pattern of Eucgr.A00514 (auxin transporter-like protein 4, LAX4) was downregulated in mat-callus compared to pri-callus in the two Eucalyptus species (additional file 1). It is consistent with the discovery of LAX3 gene in rice after grafting . The dysregulation of these genes may suggest that they might play an essential role from the beginning of wounding to encourage healing and preparation for downstream rapid development.
After maturation, shoot regeneration is another important stage indicating the successful propagation. In maize, genes encoding photosynthetic and other chloroplast components (e.g., chlorophyll a/b binding protein) were upregulated as shoots began to green . In the present study, we found that various genes encoding chloroplast components were dysregulated in SRS-callus and mat-callus in the two Eucalyptus species, including 9 genes encoding chlorophyll a-b binding proteins (additional file 1). Interestingly, these 9 genes were downregulated during the callus maturation but upregulated in the shoot regeneration process in Eucalyptus. We also found two genes (Eucgr.H00220 and Eucgr.F03055) encoding photosynthetic NDH subunit of subcomplex B 5, chloroplastic upregulated in B4 vs B3 only (additional file 1). Considering our study used the sen-callus as a contrast of SRS-callus, which may provide a deep insight into the genes associated with the shoot regeneration process of callus, there were 15 genes with diverse regulations in SRS-callus and sen-callus compared to mat-callus (Table 2), such as Eucgr.D00640, Eucgr.B00171, Eucgr.C00948, Eucgr.K01667, Eucgr.C00663 and Eucgr.C00419. Most of these 15 genes have not been reported to be associated with callus development or shoot regeneration. However, some of them have been proved to function in plant development. For example, the protein product of Eucgr.B00171 is BTB/POZ and TAZ domain-containing protein 1, which is a substrate-specific adapter of an E3 ubiquitin-protein ligase complex and involved in gametophyte development . Eucgr.C00419 encodes the UPF0481 protein At3g47200, which has been reported to be upregulated during the early flower development in Prunus mume . Eucgr.D00640 encodes the stem-specific protein TSJT1, which has been found in other plants like tobacco, grape, and soybean. The functions of these newly callus development associated genes require more experiments to investigate their roles in somatic embryogenesis, dedifferentiation, differentiation, and development.
Known vegetative propagation ability associated genes (e.g., ARF19, SERK, LEC and WUS) have been reported to play key roles during the dedifferentiation process in Eucalyptus , however, we did not find strong association between these genes and the callus development as we cannot determine their expression patterns in this process (additional file 3). However, because the two Eucalyptus species used in this study have distinct ability of vegetative propagation, we identified 51 genes that might be related to their somatic embryogenesis potential (Table 3) and 29 of them have been reported in the dedifferentiation process of Eucalyptus , including Eucgr.F02674 (putative laccase-9), Eucgr.B03374 (two-component response regulator ARR6), Eucgr.C00947 (transcription factor MYB36). These 29 genes might be triggered during the early dedifferentiation process and the other 22 genes might be induced at the callus maturation and shoot regeneration processes, such as Eucgr.A01788 (SRF3, protein STRUBBELIG-RECEPTOR FAMILY 3) and Eucgr.B02127 (GIF1, GRF1-interacting factor 1). SRF3 has been proved to be involved in the plant immunity . Recently, it was showed to be associated with the cell proliferation during the switch development from the apical buds to leaf marginal tissues . Interestingly, its homologue SRF4 was shown to play an important role in making plants display enlarged leaves through affecting cell wall formation . As a transcription coactivator, GIF1 has also been shown to control cell proliferation. Being a target of miR396, GIF1 and other GIFs act in the regulation of meristem function, at least partially through the control of cell proliferation . In addition, together with GRF5 GIF1 controls the development of appropriate leaf size and shape through the promotion and/or maintenance of cell proliferation activity in leaf primordia, GIF1 plays a role in adaxial/abaxial patterning and growth in leaf morphogenesis, and together with GATA18/HAN, GIF1 mediates the cotyledon identity by preventing ectopic root formation through the repression of PLT1 expression . We assume that the dysregulation of the 22 genes might be regulated by some mechanisms and the interaction network requires more experiments to be explored. When and how the vegetative propagation ability associated genes are expressed during the callus development are also valued research areas and will be focused in the future.
In conclusion, we analyzed the transcriptome profiles of callus tissues during the maturation and shoot regeneration processes of two Eucalyptus species which have distinct vegetative propagation ability. We observed that the regeneration rates of the tissue culture induced callus by stem peaked at 21 days of incubation on CIM. In the callus maturation process we identified 3,790 and 3,740 DEGs in E. camaldulensis and E. grandis x urophylla, respectively, including genes encoding histone/ribosomal proteins and genes involved in the plant hormone signalling transduction pathway. Then, using SRS-callus and sen-callus as contrast we identified 15 genes (e.g., Eucgr.D00640, Eucgr.B00171, Eucgr.C00948, Eucgr.K01667, Eucgr.C00663) which might play important roles during the development of mat-callus. They were annotated to encode the stem-specific protein TSJT1, zinc finger proteins, stomatal closure-related actin-binding proteins and glutaredoxin-C10 proteins. Further, 51 genes were identified to be associated with the ability of somatic embryogenesis of Eucalyptus, of which 22 genes (e.g., SRF3 and GIF1) might be induced after the dedifferentiation. This is the first time to study the transcriptome profiles of callus development in Eucalyptus. The results will improve our understanding of gene regulations and molecular mechanisms in the callus development and vegetative propagation of Eucalyptus. More importantly, the output of this study may benefit the Eucalyptus breeding program.
Plant material and culture conditions
The original seeds of E. camaldulensis (voucher ID: c0009) and E. grandis x urophylla (voucher ID: j0017) were obtained from the wild in 1984 without any restrictions. Then, the seeds and plants were confirmed by a senior botanist Prof. Dongyun Xiang and they were maintained in the experimental fields of Guangxi Forestry Research Institute. The stem tissues were obtained from the in vitro tissue-culture induced seedlings of E. camaldulensis (voucher ID: c0009, A1) and E. grandis x urophylla (B1) trees, and maintained on the callus induction MS medium (CIM, supplemented with 20mg/L Ca(NO3)2, 1 mg/L KT and 0.5 mg/L 2,4-D) for 10 days to get the pri-callus (A2, B2). Then, the pri-callus was continually incubated on the CIM for another 11 days to get the mature callus (A3, B3), which was transferred onto the shooting-inducing medium (SIM, MS medium supplemented with 20 mg/L Ca (NO3)2 +2.0 mg/L 6-BA + 0.2 mg/L NAA) for incubation. The callus was incubated on the SIM for 10 days to get the shoot regeneration stage callus (A4, B4), which were developed with some buds as successful propagation. While the mature callus was incubated on the CIM for another 14 days, we get the senescence callus (A5, B5) which totally lose the embryogenic capacity.
Total RNA extraction, cDNA library preparation and transcriptome sequencing
Total RNA was extracted from the plant tissues (A2~A5, B2~B5) using the TRIzol reagent, as previously described [19, 40]. Then, Agilent 2100 Bioanalyzer was used to evaluate the quantity and quality of the total RNA samples. Equal amount of total RNA (1 μg) was used for the cDNA library construction, as described . In brief, the poly-A mRNAs were enriched using the magnetic oligo (dT) beads and then were fragmented into 200 bp pieces. Next, random hexamer (N6) primers were used to build double strand cDNA libraries for all the samples. After the libraries were end-repaired by using phosphate at the 5’ end and sticky ‘A’ at the 3’ end, they were ligated with sequencing primers for BGISEQ-500. The libraries were sequenced on the BGISEQ-500 RS platform with paired-end 150 strategy.
Genome alignment, gene expression profiles and differential expression analysis
Raw reads were processed by SOAPnuke to remove sequencing adaptors, low quality reads and contamination reads . Then, the clean reads were aligned to the Eucalyptus genome (v2.0, https://plantgenie.org) using hisat2 . Stringtie  and Subread  were used to profile gene expression for each sample and count the read counts aligned to each gene, respectively. We next used the TPM (transcripts per million reads mapped) method to normalize gene expression in each sample and filtered lowly expressed genes (average TPM < 5). To perform differential expression analysis, we first calculated the coefficient of variation (CV) of each gene and genes with CV > 0.5 were filtered. Then, we employed edgeR with some cut-offs, including log2 fold change (log2FC) > 1 or < -1, p-value < 0.05, false discovery rate (FDR) < 0.1, to identify differentially expressed genes in two samples.
We next annotated the Eucalyptus genes by mapping them to the Gene Ontology (GO) and KEGG pathway databases, as previously described . Then, the enriched GO terms and KEGG pathways by differentially expressed genes were identified by p-value (< 0.05), calculated by Fisher’s exact test, and q-value (< 0.05), calculated by the R package ‘qvalue’.
We used the R package “WGCNA” to identify co-expressed genes during the callus development of the two Eucalyptus species , according to the manufacturer’s protocol.
We selected 12 genes for qPCR-PCR validation and used the H2B gene as the internal control. The primers of these 13 genes were predicted by Primer3 and synthesized at BGI-Shenzhen. The procedure of qRT-PCR was same as a previous study [19, 46]. For each gene in a biological replicate, we performed three qRT-PCR reactions and in total we have 9 replicates for one gene (n=9). After the Ct values were calculated and averaged, we used the ΔCt value to present the gene expression in each sample. Then, using the mature callus as the control sample we calculated the ΔΔCt value to show the difference of a gene in the callus development. Last, relative normalized expression (RNE) was used to show the gene expression change: RNE = 2−ΔΔCt and log2RNE was used to match the transcriptome sequencing method. For pri-callus and mature callus comparison (A3 vs A2 and B3 vs B2), we used -log2RNE to present the gene changes.
Availability of data and materials
The raw sequencing data can be accessed from the NCBI Sequence Read Archive (SRA) platform (http://trace.ncbi.nlm.nih.gov/Traces/sra/) under the accession number PRJNA761197.
Ikeuchi M, Sugimoto K, Iwase A. Plant callus: mechanisms of induction and repression. Plant Cell. 2013;25(9):3159–73.
Zimmerman JL. Somatic Embryogenesis: A Model for Early Development in Higher Plants. Plant Cell. 1993;5(10):1411–23.
Frank M, Rupp HM, Prinsen E, Motyka V, Van Onckelen H, Schmulling T. Hormone autotrophic growth and differentiation identifies mutant lines of Arabidopsis with altered cytokinin and auxin content or signaling. Plant Physiol. 2000;122(3):721–9.
Skoog F, Miller CO. Chemical regulation of growth and organ formation in plant tissues cultured in vitro. Symp Soc Exp Biol. 1957;11:118–30.
Lee DK, Geisler M, Springer PS. LATERAL ORGAN FUSION1 and LATERAL ORGAN FUSION2 function in lateral organ separation and axillary meristem formation in Arabidopsis. Development. 2009;136(14):2423–32.
Okushima Y, Fukaki H, Onoda M, Theologis A, Tasaka M. ARF7 and ARF19 regulate lateral root formation via direct activation of LBD/ASL genes in Arabidopsis. Plant Cell. 2007;19(1):118–30.
Sakai H, Honma T, Aoyama T, Sato S, Kato T, Tabata S, et al. ARR1, a transcription factor for genes immediately responsive to cytokinins. Science. 2001;294(5546):1519–21.
Tajima Y, Imamura A, Kiba T, Amano Y, Yamashino T, Mizuno T. Comparative studies on the type-B response regulators revealing their distinctive properties in the His-to-Asp phosphorelay signal transduction of Arabidopsis thaliana. Plant Cell Physiol. 2004;45(1):28–39.
Koszegi D, Johnston AJ, Rutten T, Czihal A, Altschmied L, Kumlehn J, et al. Members of the RKD transcription factor family induce an egg cell-like gene expression program. Plant J. 2011;67(2):280–91.
Waki T, Hiki T, Watanabe R, Hashimoto T, Nakajima K. The Arabidopsis RWP-RK protein RKD4 triggers gene expression and pattern formation in early embryogenesis. Curr Biol. 2011;21(15):1277–81.
Tan BC, Chin CF, Liddell S, Alderson P. Proteomic analysis of callus development in Vanilla planifolia Andrews. Plant molecular biology reporter. 2013;31(6):1220–9.
Che P, Lall S, Nettleton D, Howell SH. Gene expression programs during shoot, root, and callus development in Arabidopsis tissue culture. Plant Physiol. 2006;141(2):620–37.
Durand-Cresswell R, Boulay M, Franclet A: Vegetative Propagation of Eucalyptus. In: Tissue Culture in Forestry. Edited by Bonga JM, Durzan DJ. Dordrecht: Springer Netherlands; 1982: 150-181.
De Fossard R, Bennett M, Gorst J, Bourne R. Tissue culture propagation of Eucalyptus ficifolia F. In: Muell. International Plant Propagators' Society: Combined proceedings, International Plant Propagators', vol. 28: Society; 1978. p. 427–35.
Grattapaglia D, Bertolucci FL, Sederoff RR. Genetic mapping of QTLs controlling vegetative propagation in Eucalyptus grandis and E. urophylla using a pseudo-testcross strategy and RAPD markers. Theor Appl Genet. 1995;90(7-8):933–47.
Marques C, Vasquez-Kool J, Carocha V, Ferreira J, O’Malley D, Liu B-H, et al. Genetic dissection of vegetative propagation traits in Eucalyptus tereticornis and E. globulus. Theoretical and applied genetics. 1999;99(6):936–46.
Mohamed E, Claudio S, Aiming W. Molecular regulation of plant somatic embryogenesis. In Vitro Cellular & Developmental Biology-Plant. 2013.
Yang X, Zhang X. Regulation of Somatic Embryogenesis in Higher Plants. Critical Reviews in Plant Sciences. 2010;29(1):36–57.
Xiao Y, Li J, Zhang Y, Zhang X, Liu H, Qin Z, et al. Transcriptome analysis identifies genes involved in the somatic embryogenesis of Eucalyptus. BMC Genomics. 2020;21(1):803.
Chen ZZ, Ho CK, Ahn IS, Chiang VL. Eucalyptus. Methods Mol Biol. 2006;344:125–34.
Che P, Love TM, Frame BR, Wang K, Carriquiry AL, Howell SH. Gene expression patterns during somatic embryo development and germination in maize Hi II callus cultures. Plant Mol Biol. 2006;62(1-2):1–14.
Othmani A, Bayoudh C, Drira N, Marrakchi M, Trifi M. Somatic embryogenesis and plant regeneration in date palm Phœnix dactylifera L., cv. Boufeggous is significantly improved by fine chopping and partial desiccation of embryogenic callus. Plant Cell, Tissue and Organ Culture (PCTOC). 2009;97(1):71–9.
Kamo KK, Becwar MR, Hodges TK. Regeneration of Zea mays L. from embryogenic callus. Botanical gazette. 1985;146(3):327–34.
Zhang K, Novak O, Wei Z, Gou M, Zhang X, Yu Y, et al. Arabidopsis ABCG14 protein controls the acropetal translocation of root-synthesized cytokinins. Nat Commun. 2014;5:3274.
Ko D, Kang J, Kiba T, Park J, Kojima M, Do J, et al. Arabidopsis ABCG14 is essential for the root-to-shoot translocation of cytokinin. Proc Natl Acad Sci U S A. 2014;111(19):7150–5.
Borghi L, Kang J, Ko D, Lee Y, Martinoia E. The role of ABCG-type ABC transporters in phytohormone transport. Biochem Soc Trans. 2015;43(5):924–30.
Wang S, Wang S, Sun Q, Yang L, Zhu Y, Yuan Y, et al. A Role of Cytokinin Transporter in Arabidopsis Immunity. Mol Plant Microbe Interact. 2017;30(4):325–33.
Terasaka K, Blakeslee JJ, Titapiwatanakun B, Peer WA, Bandyopadhyay A, Makam SN, et al. PGP4, an ATP binding cassette P-glycoprotein, catalyzes auxin transport in Arabidopsis thaliana roots. Plant Cell. 2005;17(11):2922–39.
Santelia D, Vincenzetti V, Azzarello E, Bovet L, Fukao Y, Duchtig P, et al. MDR-like ABC transporter AtPGP4 is involved in auxin-mediated lateral root and root hair development. FEBS Lett. 2005;579(24):5399–406.
Lee Y, Kende H. Expression of beta-expansins is correlated with internodal elongation in deepwater rice. Plant Physiol. 2001;127(2):645–54.
Zou H, Wenwen Y, Zang G, Kang Z, Zhang Z, Huang J, et al. OsEXPB2, a β-expansin gene, is involved in rice root system architecture. Molecular Breeding. 2015;35(1):41.
Yousef AF, Ali MM, Rizwan HM, Gad AG, Liang D, Binqi L, et al. Light quality and quantity affect graft union formation of tomato plants. Sci Rep. 2021;11(1):9870.
Robert HS, Quint A, Brand D, Vivian-Smith A, Offringa R. BTB and TAZ domain scaffold proteins perform a crucial function in Arabidopsis development. Plant J. 2009;58(1):109–21.
Wu X, Shi T, Iqbal S, Zhang Y, Liu L, Gao Z. Genome-wide discovery and characterization of flower development related long non-coding RNAs in Prunus mume. BMC Plant Biol. 2019;19(1):64.
Alcazar R, Garcia AV, Kronholm I, de Meaux J, Koornneef M, Parker JE, et al. Natural variation at Strubbelig Receptor Kinase 3 drives immune-triggered incompatibilities between Arabidopsis thaliana accessions. Nat Genet. 2010;42(12):1135–9.
Bian X, Qu C, Zhang M, Li D, Jiang J, Liu G. Transcriptome analysis provides new insights into leaf shape variation in birch. Trees. 2019;33(5):1265–81.
Eyuboglu B, Pfister K, Haberer G, Chevalier D, Fuchs A, Mayer KF, et al. Molecular characterisation of the STRUBBELIG-RECEPTOR FAMILY of genes encoding putative leucine-rich repeat receptor-like kinases in Arabidopsis thaliana. BMC Plant Biol. 2007;7:16.
Fujikura U, Horiguchi G, Ponce MR, Micol JL, Tsukaya H. Coordination of cell proliferation and cell expansion mediated by ribosome-related processes in the leaves of Arabidopsis thaliana. Plant J. 2009;59(3):499–508.
Kanei M, Horiguchi G, Tsukaya H. Stable establishment of cotyledon identity during embryogenesis in Arabidopsis by ANGUSTIFOLIA3 and HANABA TARANU. Development. 2012;139(13):2436–46.
Chen M, Xu R, Rai A, Suwakulsiri W, Izumikawa K, Ishikawa H, et al. Distinct shed microvesicle and exosome microRNA signatures reveal diagnostic markers for colorectal cancer. PLOS ONE. 2019;14(1):e0210003.
Chen Y, Chen Y, Shi C, Huang Z, Zhang Y, Li S, et al. SOAPnuke: a MapReduce acceleration-supported software for integrated quality control and preprocessing of high-throughput sequencing data. Gigascience. 2018;7(1):1–6.
Pertea M, Kim D, Pertea GM, Leek JT, Salzberg SL: Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and Ballgown. Nat Protoc 2016, 11(9):1650-1667.
Liao Y, Smyth GK, Shi W. The Subread aligner: fast, accurate and scalable read mapping by seed-and-vote. Nucleic Acids Res. 2013;41(10):e108.
Wei S, Ma X, Pan L, Miao J, Fu J, Bai L, et al. Transcriptome Analysis of Taxillusi chinensis (DC.) Danser Seeds in Response to Water Loss. PLoS One. 2017;12(1):e0169177.
Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics. 2008;9:559.
Wei S, Wan L, He L, Wei Y, Long H, Ji X, et al. De novo Transcriptome Reveals Gene Changes in the Development of the Endosperm Chalazal Haustorium in Taxillus chinensis (DC.) Danser. Biomed Res Int. 2020;2020(7871918).
We acknowledge Mr. Dajie Zhou from Beijing Genomics Institute at Shenzhen for his help on the deep sequencing.
This work was financially supported by the National Natural Science Foundation of China (32160391) and the Guangxi Science and Technology Program (GuiKeAD18281086). The authors declare that the funding bodies have no role in the research design, the data collection and analysis, and the manuscript preparation.
Ethics approval and consent to participate
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.
Consent for publication
The authors declare no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Zhang, Y., Li, J., Li, C. et al. Gene expression programs during callus development in tissue culture of two Eucalyptus species. BMC Plant Biol 22, 1 (2022). https://doi.org/10.1186/s12870-021-03391-x
- vegetative propagation
- somatic embryogenesis
- transcription factor