Comparative transcriptome analysis between somatic embryos (SEs) and zygotic embryos in cotton: evidence for stress response functions in SE development

Proteomic and Phosphoproteomic Analyses during Plant Regeneration Initiation in Cotton (Gossypium hirsutum L.)

Somatic embryogenesis (SE) is a biotechnological tool used to generate new individuals and is the preferred method for rapid plant regeneration. However, the molecular basis underlying somatic cell regeneration through SE is not yet fully understood, particularly regarding interactions between the proteome and post-translational modifications. Here, we performed association analysis of high-throughput proteomics and phosphoproteomics in three representative samples (non-embryogenic calli, NEC; primary embryogenic calli, PEC; globular embryos, GE) during the initiation of plant regeneration in cotton, a pioneer crop for genetic biotechnology applications. Our results showed that protein accumulation is positively regulated by phosphorylation during SE, as revealed by correlation analyses. Of the 1418 proteins that were differentially accumulated in the proteome and the 1106 phosphoproteins that were differentially regulated in the phosphoproteome, 115 proteins with 229 phosphorylation sites overlapped (co-differential). Furthermore, seven dynamic trajectory patterns of differentially accumulated proteins (DAPs) and the correlated differentially regulated phosphoproteins (DRPPs) pairs with enrichment features were observed. During the initiation of plant regeneration, functional enrichment analysis revealed that the overlapping proteins (DAPs-DRPPs) were considerably enriched in cellular nitrogen metabolism, spliceosome formation, and reproductive structure development. Moreover, 198 DRPPs (387 phosphorylation sites) were specifically regulated at the phosphorylation level and showed four patterns of stage-enriched phosphorylation susceptibility. Furthermore, enrichment annotation analysis revealed that these phosphoproteins were significantly enriched in endosomal transport and nucleus organization processes. During embryogenic differentiation, we identified five DAPs-DRPPs with significantly enriched characteristic patterns. These proteins may play essential roles in transcriptional regulation and signaling events that initiate plant regeneration through protein accumulation and/or phosphorylation modification. This study enriched the understanding of key proteins and their correlated phosphorylation patterns during plant regeneration, and also provided a reference for improving plant regeneration efficiency.

Single-cell transcriptome atlas reveals somatic cell embryogenic differentiation features during regeneration

Understanding somatic cell totipotency remains a challenge facing scientific inquiry today. Plants display remarkable cell totipotency expression, illustrated by single-cell differentiation during somatic embryogenesis (SE) for plant regeneration. Determining cell identity and exploring gene regulation in such complex heterogeneous somatic cell differentiation have been major challenges. Here, we performed high-throughput single-cell sequencing assays to define the precise cellular landscape and revealed the modulation mode of marker genes during embryogenic differentiation in cotton (Gossypium hirsutum L.) as the crop for biotechnology application. We demonstrated that nonembryogenic calli (NEC) and primary embryogenic calli (PEC) tissues were composed of heterogeneous cells that could be partitioned into four broad populations with six distinct cell clusters. Enriched cell clusters and cell states were identified in NEC and PEC samples, respectively. Moreover, a broad repertoire of new cluster-specific genes and associated expression modules were identified. The energy metabolism, signal transduction, environmental adaptation, membrane transport pathways, and a series of transcription factors were preferentially enriched in cell embryogenic totipotency expression. Notably, the SE-ASSOCIATED LIPID TRANSFER PROTEIN (SELTP) gene dose-dependently marked cell types with distinct embryogenic states and exhibited a parabolic curve pattern along the somatic cell embryogenic differentiation trajectory, suggesting that SELTP could serve as a favorable quantitative cellular marker for detecting embryogenic expression at the single-cell level. In addition, RNA velocity and Scissor analysis confirmed the pseudo-temporal model and validated the accuracy of the scRNA-seq data, respectively. This work provides valuable marker-genes resources and defines precise cellular taxonomy and trajectory atlases for somatic cell embryogenic differentiation in plant regeneration.

Transcriptome analysis reveals the effect of propyl gallate on kiwifruit callus formation

KEY MESSAGE

Exploring the potential action mechanisms of reactive oxygen species during the callus inducing, they can activate specific metabolic pathways in explants to regulate callus development. Reactive oxygen species (ROS) play an important role in the regulation of plant growth and development, but the mechanism of their action on plant callus formation remains to be elucidated. To address this question, kiwifruit was selected as the explant for callus induction, and the influence of ROS on callus formation was investigated by introducing propyl gallate (PG) as an antioxidant into the medium used for inducing callus. The results have unveiled that the inclusion of PG in the medium has disturbed the equilibrium of ROS during the formation of the kiwifruit callus. We selected the callus that was induced by the addition of 0.05 mmol/L PG to the MS medium. The callus exhibited a significant difference in the amount compared to the control medium without PG. The callus induced by the MS medium without PG was used as the control for comparison. KEGG enrichment indicated that PG exposure resulted in significant differences in gene expression in related pathways, such as phytohormone signaling and glutathione in kiwifruit callus. Weighted gene co-expression analysis indicated that the pertinent regulatory networks of both ROS and phytohormone signaling were critical for the establishment of callus in kiwifruit leaves. In addition, during the process of callus establishment, the ROS level of the explants was also closely related to the genes for transmembrane transport of substances, cell wall formation, and plant organ establishment. This investigation expands the theory of ROS-regulated callus formation and presents a new concept for the expeditious propagation of callus in kiwifruit.

Differential expression of genes potentially related to the callogenesis and in situ hybridization of SERK gene in macaw palm (Acrocomia aculeata Jacq.) Lodd. ex Mart

For the purpose of understanding the molecular processes triggered during callus formation in macaw palm, the expression of seven genes potentially involved in this process, identified in previous studies and from the literature, was investigated by RT-qPCR. In addition, in situ hybridization of the SERK gene was performed. Leaf tissues from adult plants from two macaw palm accession were inoculated in a medium combined with Picloram at a concentration of 450 μM to induce callus. The expression analysis was performed from leaf samples from two accessions of different origins (Municipalities of Tiros, MG, and Buriti Vermelho, DF, Brazil), which are characterized as non-responsive (NR) and responsive (R), respectively. The material was collected before callus induction (0 DAI, initial day) and 120 days after callus induction (120 DAI). Genes related to development (SERK, OASA, EF1, ANN1) and stress (LEA, CAT2, and MDAR5) were evaluated. The results obtained showed that all the genes involved with the development had their expressions downregulated at 0 DAI when the accession R was compared with the accession NR. On the other hand, it was possible to observe that these genes were upregulated at 120 DAI. The LEA stress gene showed a tendency to increase expression in the NR accession, while the R accession showed decreased expression and the CAT2 and MDAR5 genes showed upregulation in both accessions. In situ hybridization showed SERK transcripts in the vascular bundles, indicating the expression of SERK in this region, in addition to its expression in calluses. The results obtained in this study support our hypothesis that the regulation of genes involved in the control of oxidative stress and development is crucial for the formation of calluses in macaw palm.

Structure-activity relationship of 2,4-D correlates auxinic activity with the induction of somatic embryogenesis in Arabidopsis thaliana

2,4-dichlorophenoxyacetic acid (2,4-D) is a synthetic analogue of the plant hormone auxin that is commonly used in many in vitro plant regeneration systems, such as somatic embryogenesis (SE). Its effectiveness in inducing SE, compared to the natural auxin indole-3-acetic acid (IAA), has been attributed to the stress triggered by this compound rather than its auxinic activity. However, this hypothesis has never been thoroughly tested. Here we used a library of forty 2,4-D analogues to test the structure-activity relationship with respect to the capacity to induce SE and auxinic activity in Arabidopsis thaliana. Four analogues induced SE as effectively as 2,4-D and 13 analogues induced SE but were less effective. Based on root growth inhibition and auxin response reporter expression, the 2,4-D analogues were classified into different groups, ranging from very active to not active auxin analogues. A halogen at the 4-position of the aromatic ring was important for auxinic activity, whereas a halogen at the 3-position resulted in reduced activity. Moreover, a small substitution at the carboxylate chain was tolerated, as was extending the carboxylate chain with an even number of carbons. The auxinic activity of most 2,4-D analogues was consistent with their simulated TIR1-Aux/IAA coreceptor binding characteristics. A strong correlation was observed between SE induction efficiency and auxinic activity, which is in line with our observation that 2,4-D-induced SE and stress both require TIR1/AFB auxin co-receptor function. Our data indicate that the stress-related effects triggered by 2,4-D and considered important for SE induction are downstream of auxin signalling.

Molecular mechanism of somatic embryogenesis in paeonia ostii 'Fengdan' based on transcriptome analysis combined histomorphological observation and metabolite determination

BACKGROUND

Tree peony (Paeonia sect. Moutan DC.) is a famous flower native to China with high ornamental, medicinal, and oil value. However, the low regeneration rate of callus is one of the main constraints for the establishment of a genetic transformation system in tree peony. By histomorphological observation, transcriptomic analysis and metabolite determination, we investigated the molecular mechanism of somatic embryogenesis after the establishment of a culture system and the induction of somatic embryo(SE) formation.

RESULTS

We found that SE formation was successfully induced when cotyledons were used as explants. A total of 3185 differentially expressed genes were screened by comparative transcriptomic analysis of embryogenic callus (EC), SE, and non-embryogenic callus (NEC). Compared to NEC, the auxin synthesis-related genes GH3.6 and PCO2 were up-regulated, whereas cytokinin dehydrogenase (CKX6) and CYP450 family genes were down-regulated in somatic embryogenesis. In SE, the auxin content was significantly higher than the cytokinin content. The methyltransferase-related gene S-adenosylmethionine synthase (SAMS) and the flavonoid biosynthesis-related gene (ANS and F3'5'H) were down-regulated in somatic embryogenesis. The determination of flavonoids showed that rhoifolin and hyperoside had the highest content in SE. The results of transcriptome analysis were consistent with the relative expression of 8 candidate genes by quantitative polymerase chain reaction analysis.

CONCLUSION

The results revealed that auxin and cytokinin may play a key role in 'Fengdan' somatic embryogenesis. The genes related to somatic embryogenesis were revealed, which has partly elucidated the molecular mechanism of somatic embryogenesis in 'Fengdan'.

Comparative transcriptome analysis of high- and low-embryogenic Hevea brasiliensis genotypes reveals involvement of phytohormones in somatic embryogenesis

BACKGROUND

Rubber plant (Hevea brasiliensis) is one of the major sources of latex. Somatic embryogenesis (SE) is a promising alterative to its propagation by grafting and seed. Phytohormones have been shown to influence SE in different plant species. However, limited knowledge is available on the role of phytohormones in SE in Hevea. The anther cultures of two Hevea genotypes (Yunyan 73477-YT and Reken 628-RT) with contrasting SE rate were established and four stages i.e., anthers (h), anther induced callus (y), callus differentiation state (f), and somatic embryos (p) were studied. UPLC-ESI-MS/MS and transcriptome analyses were used to study phytohormone accumulation and related expression changes in biosynthesis and signaling genes.

RESULTS

YT showed higher callus induction rate than RT. Of the two genotypes, only YT exhibited successful SE. Auxins, cytokinins (CKs), abscisic acid (ABA), jasmonic acid (JA), salicylic acid (SA), gibberellins (GAs), and ethylene (ETH) were detected in the two genotypes. Indole-3-acetic acid (IAA), CKs, ABA, and ETH had notable differences in the studied stages of the two genotypes. The differentially expressed genes identified in treatment comparisons were majorly enriched in MAPK and phytohormone signaling, biosynthesis of secondary metabolites, and metabolic pathways. The expression changes in IAA, CK, ABA, and ETH biosynthesis and signaling genes confirmed the differential accumulation of respective phytohormones in the two genotypes.

CONCLUSION

These results suggest potential roles of phytohormones in SE in Hevea.

Shoot differentiation from Dendrocalamus brandisii callus and the related physiological roles of sugar and hormones during shoot differentiation

Only a few calli regeneration systems of bamboos were successfully established, which limited the research on physiological mechanism of callus differentiation. In this study, we successfully established the callus differentiation systems of Dendrocalamus brandisii via seeds. The results showed that the best medium for callus induction of D. brandisii seeds was basal MS media amended with 5.0 mg L-1 2,4-D and 0.5 mg L-1 KT, and the optimal medium for shoot differentiation was the basal MS media supplemented with 4.0 mg L-1 BA and 0.5 mg L-1 NAA. Callus tissues had apparent polarity in cell arrangement, and developed their own meristematic cell layers. α-amylase, STP and SUSY played a dominant role in carbohydrates degradation in callus during shoot differentiation. PPP and TCA pathways up-regulated in the shoot-differentiated calli. The dynamics of BA and KT contents in calli was consistent with their concentrations applied in medium. IAA synthesis and the related signal transduction were down-regulated, while the endogenous CTKs contents were up-regulated by the exogenous CTKs application in shoot-differentiated calli, and their related synthesis, transport and signal transduction pathways were also up-regulated. The downregulated signal transduction pathways of IAA and ABA revealed that they did not play the key role in shoot differentiation of bamboos. GAs also played a role in shoot differentiation based on the down-regulation of DELLA and the up-regulation of PIF4 genes. The overexpression of DbSNRK2 and DbFIF4 genes further confirmed the negative role of ABA and the positive role of GAs in shoot differentiation.

Proteome profiles during early stage of somatic embryogenesis of two Eucalyptus species

BACKGROUND

Somatic embryogenesis (SE) was recognized as an important tool for plants to propagate. However, our knowledge about the proteins involved in early SE including the callus dedifferentiation is still limited, especially in the economic woody tree - Eucalyptus.

RESULTS

We used the data-independent acquisition mass-spectrometry to study the different proteome profiles of early SE of two Eucalyptus species-E. camaldulensis (high regeneratively potential) and E. grandis x urophylla (low regenerative potential). Initially, 35,207 peptides and 7,077 proteins were identified in the stem and tissue-culture induced callus of the two Eucalyptus species. MSstat identified 2,078 and 2,807 differentially expressed proteins (DEPs) in early SE of E. camaldulensis and E. grandis x urophylla, respectively. They shared 760 upregulated and 420 downregulated proteins, including 4 transcription factors, 31 ribosomal proteins, 1 histone, 3 zinc finger proteins (ZFPs), 16 glutathione transferases, 10 glucosyltransferases, ARF19, WOX8 and PIN1. These proteins might be involved in the early SE of Eucalyptus. By combining the miRNA and RNA-Seq results, some miRNA ~ gene/protein regulatory networks were identified in early SE of Eucalyptus, such as miR160 ~ TPP2, miR164 ~ UXS2, miR169 ~ COX11 and miR535 ~ Eucgr.E01067. Further, we found SERK, WRKY, ZFP and ABC transporter might be related with high SE potential.

CONCLUSIONS

Overall, our study identified proteins involved in the early SE and related to the high regeneration potential of Eucalyptus. It greatly enhanced our understanding of the early SE and the SE capacity of Eucalyptus.

Integrated transcriptomic and metabolic analyses provide insights into the maintenance of embryogenic potential and the biosynthesis of phenolic acids and flavonoids involving transcription factors in Larix kaempferi (Lamb.) Carr

Somatic embryogenesis (SE) techniques have been established for micropropagation or basic research related to plant development in many conifer species. The frequent occurrence of non-embryogenic callus (NEC) during SE has impose constraints on the application of somatic embryogenesis SE in Larix kaempferi (Lamb.) Carr, but the potential regulatory mechanisms are poorly understood. In this study, integrated transcriptomic and metabolomic analyses were performed in embryogenic callus (EC) and NEC originating from a single immature zygotic embryo to better decipher the key molecular and metabolic mechanisms required for embryogenic potential maintenance. The results showed that a total of 13,842 differentially expressed genes (DEGs) were found in EC and NEC, among which many were enriched in plant hormone signal transduction, starch and sucrose metabolism, phenylpropanoid biosynthesis, flavonoid biosynthesis, and the biosynthesis of amino acids pathways. Metabolite profiling showed that 441 differentially accumulated metabolites (DAMs) were identified in EC and NEC. Both EC and NEC had vigorous primary metabolic activities, while most secondary metabolites were upregulated in NEC. Many totipotency-related transcription factor (TF) genes such as BBMs, WUSs, and LEC1 showed higher expression levels in EC compared with NEC, which may result in the higher accumulation of indole 3-acetic acid (IAA) in EC. NEC was characterized by upregulation of genes and metabolites associated with stress responses, such as DEGs involved in jasmonic acid (JA) and ethylene (ETH) biosynthesis and signal transduction pathways, and DEGs and DAMs related to phenylpropanoid and flavonoid biosynthesis. We predicted and analyzed TFs that could target several key co-expressed structural DEGs including two C4H genes, two CcoAOMT genes and three HCT genes involved in phenylpropanoid and flavonoid biosynthesis. Based on the targeted relationship and the co-expression network, two ERFs (Lk23436 and Lk458687), one MYB (Lk34626) and one C2C2-dof (Lk37167) may play an important role in regulating phenolic acid and flavonoid biosynthesis by transcriptionally regulating the expression of these structural genes. This study shows an approach involving integrated transcriptomic and metabolic analyses to obtain insights into molecular events underlying embryogenic potential maintenance and the biosynthesis mechanisms of key metabolites involving TF regulation, which provides valuable information for the improvement of SE efficiency in L. kaempferi.

Transcriptomic and Metabolomic Analysis of Korean Pine Cell Lines with Different Somatic Embryogenic Potential

The embryogenesis capacity of conifer callus is not only highly genotype-dependent, but also gradually lost after long-term proliferation. These problems have seriously limited the commercialization of conifer somatic embryogenesis (SE) technology. In this study, the responsive SE cell line (R-EC), the blocked SE cell line (B-EC), and the loss of SE cell line (L-EC) were studied. The morphological, physiological, transcriptomic, and metabolomic profiles of these three types of cells were analyzed. We found that R-EC had higher water content, total sugar content, and putrescine (Put) content, as well as lower superoxide dismutase (SOD) activity and H2O2 content compared to B-EC and L-EC. A total of 2566, 13,768, and 13,900 differentially expressed genes (DEGs) and 219, 253, and 341 differentially expressed metabolites (DEMs) were found in the comparisons of R-EC versus B-EC, R-EC versus B-EC, and B-EC versus L-EC, respectively. These DEGs and DEMs were mainly found to be involved in plant signal transduction, starch and sugar metabolism, phenylpropane metabolism, and flavonoid metabolism. We found that the AUX1 and AUX/IAA families of genes were significantly up-regulated after the long-term proliferation of callus, resulting in higher auxin content. Most phenylpropane and flavonoid metabolites, which act as antioxidants to protect cells from damage, were found to be significantly up-regulated in R-EC.

GhTCE1-GhTCEE1 dimers regulate transcriptional reprogramming during wound-induced callus formation in cotton

Wounded plant cells can form callus to seal the wound site. Alternatively, wounding can cause adventitious organogenesis or somatic embryogenesis. These distinct developmental pathways require specific cell fate decisions. Here, we identify GhTCE1, a basic helix-loop-helix family transcription factor, and its interacting partners as a central regulatory module of early cell fate transition during in vitro dedifferentiation of cotton (Gossypium hirsutum). RNAi- or CRISPR/Cas9-mediated loss of GhTCE1 function resulted in excessive accumulation of reactive oxygen species (ROS), arrested callus cell elongation, and increased adventitious organogenesis. In contrast, GhTCE1-overexpressing tissues underwent callus cell growth, but organogenesis was repressed. Transcriptome analysis revealed that several pathways depend on proper regulation of GhTCE1 expression, including lipid transfer pathway components, ROS homeostasis, and cell expansion. GhTCE1 bound to the promoters of the target genes GhLTP2 and GhLTP3, activating their expression synergistically, and the heterodimer TCE1-TCEE1 enhances this activity. GhLTP2- and GhLTP3-deficient tissues accumulated ROS and had arrested callus cell elongation, which was restored by ROS scavengers. These results reveal a unique regulatory network involving ROS and lipid transfer proteins, which act as potential ROS scavengers. This network acts as a switch between unorganized callus growth and organized development during in vitro dedifferentiation of cotton cells.

An overview of salinity stress, mechanism of salinity tolerance and strategies for its management in cotton

Salinity stress is one of the primary threats to agricultural crops resulting in impaired crop growth and development. Although cotton is considered as reasonably salt tolerant, it is sensitive to salt stress at some critical stages like germination, flowering, boll formation, resulting in reduced biomass and fiber production. The mechanism of partial ion exclusion (exclusion of Na+ and/or Cl-) in cotton appears to be responsible for the pattern of uptake and accumulation of harmful ions (Na+ and Cl) in tissues of plants exposed to saline conditions. Maintaining high tissue K+/Na+ and Ca2+/Na+ ratios has been proposed as a key selection factor for salt tolerance in cotton. The key adaptation mechanism in cotton under salt stress is excessive sodium exclusion or compartmentation. Among the cultivated species of cotton, Egyptian cotton (Gossypium barbadense L.) exhibit better salt tolerance with good fiber quality traits as compared to most cultivated cotton and it can be used to improve five quality traits and transfer salt tolerance into Upland or American cotton (Gossypium hirsutum L.) by interspecific introgression. Cotton genetic studies on salt tolerance revealed that the majority of growth, yield, and fiber traits are genetically determined, and controlled by quantitative trait loci (QTLs). Molecular markers linked to genes or QTLs affecting key traits have been identified, and they could be utilized as an indirect selection criterion to enhance breeding efficiency through marker-assisted selection (MAS). Transfer of genes for compatible solute, which are an important aspect of ion compartmentation, into salt-sensitive species is, theoretically, a simple strategy to improve tolerance. The expression of particular stress-related genes is involved in plant adaptation to environmental stressors. As a result, enhancing tolerance to salt stress can be achieved by marker assisted selection added with modern gene editing tools can boost the breeding strategies that defend and uphold the structure and function of cellular components. The intent of this review was to recapitulate the advancements in salt screening methods, tolerant germplasm sources and their inheritance, biochemical, morpho-physiological, and molecular characteristics, transgenic approaches, and QTLs for salt tolerance in cotton.

Weighted Gene Correlation Network Analysis (WGCNA) of Arabidopsis Somatic Embryogenesis (SE) and Identification of Key Gene Modules to Uncover SE-Associated Hub Genes

Somatic embryogenesis (SE), which occurs naturally in many plant species, serves as a model to elucidate cellular and molecular mechanisms of embryo patterning in plants. Decoding the regulatory landscape of SE is essential for its further application. Hence, the present study was aimed at employing Weighted Gene Correlation Network Analysis (WGCNA) to construct a gene coexpression network (GCN) for Arabidopsis SE and then identifying highly correlated gene modules to uncover the hub genes associated with SE that may serve as potential molecular targets. A total of 17,059 genes were filtered from a microarray dataset comprising four stages of SE, i.e., stage I (zygotic embryos), stage II (proliferating tissues at 7 days of induction), stage III (proliferating tissues at 14 days of induction), and stage IV (mature somatic embryos). This included 1,711 transcription factors and 445 EMBRYO DEFECTIVE genes. GCN analysis identified a total of 26 gene modules with the module size ranging from 35 to 3,418 genes using a dynamic cut tree algorithm. The module-trait analysis revealed that four, four, seven, and four modules were associated with stages I, II, III, and IV, respectively. Further, we identified a total of 260 hub genes based on the degree of intramodular connectivity. Validation of the hub genes using publicly available expression datasets demonstrated that at least 78 hub genes are potentially associated with embryogenesis; of these, many genes remain functionally uncharacterized thus far. In silico promoter analysis of these genes revealed the presence of cis-acting regulatory elements, "soybean embryo factor 4 (SEF4) binding site," and "E-box" of the napA storage-protein gene of Brassica napus; this suggests that these genes may play important roles in plant embryo development. The present study successfully applied WGCNA to construct a GCN for SE in Arabidopsis and identified hub genes involved in the development of somatic embryos. These hub genes could be used as molecular targets to further elucidate the molecular mechanisms underlying SE in plants.

miR171 modulates induction of somatic embryogenesis in citrus callus

Overexpression of miR171 restored SE competence in the recalcitrant citrus callus, and inhibition of miR171 function weakened SE competence in the strongly embryogenic citrus callus. Somatic embryogenesis (SE) is an important way of in vitro regeneration for plants. For perennial woody crops such as citrus, embryogenic callus is usually induced from unfertilized aborted ovules and widely used in biotechnology aided breeding. However, SE capacity always declines in callus during subculture, which makes regeneration difficult and hinders the application of biotechnology. We previously found that miR171 may be a regulator of SE in citrus, based on the abundant expression of csi-miR171c in the embryogenic callus and during SE of citrus. Here, we report that miR171 promotes SE and is required for SE in citrus. Overexpression of miR171 restored SE competence in the recalcitrant callus of 'Guoqing No.1' Satsuma mandarin (G1), whereas inhibition of miR171 function by Short Tandem Target Mimic (STTM) weakened SE competence in the strongly embryogenic callus of 'Valencia' sweet orange (V). The comparative transcriptomic analysis in miR171 overexpressed callus line (OE) and the wild type callus (WT) indicated that overexpression of miR171 decreased the expression level of its SCARECROW-LIKE (CsSCL) targets, and activated stress response related biological processes and metabolic processes that are required for cell differentiation. However, CsSCLs were up-regulated in the OE callus during SE induction process, which activated the cell division and developmental processes that are required for embryogenesis progress. Our results validate the function of miR171 in regulation of SE and reveal the biological responses provoked by miR171 in citrus that may promote SE.

Uniconazole Augments Abscisic Acid in Promoting Somatic Embryogenesis in Cotton (Gossypium hirsutum L.)

During somatic embryogenesis (SE), somatic cells initiate embryogenic development under appropriate conditions. Uniconazole, a plant growth regulator, was found to inhibit the proliferation of callus but promoted the conversion of callus into an embryogenic callus (EC) in cotton. The supplementation of uniconazole in the culture medium significantly suppressed the endogenous auxin [indole acetic acid (IAA)] level in callus tissues in both the callus initiation and proliferation stage but enhanced the abscisic acid (ABA) level only in the callus proliferation stage. Exogenous ABA and uniconazole showed cooperative effects on promoting the differentiation rate of callus into EC. These findings were verified by RNA-seq analysis, which elucidated that the genes involved in the IAA biosynthesis, metabolism, and signaling, and ABA metabolism pathways were regulated by uniconazole during the callus development and SE. Overall, the results suggest that uniconazole could modulate callus proliferation and callus differentiation rate by regulating the endogenous levels of IAA and ABA.

Transcriptome analysis reveals key developmental and metabolic regulatory aspects of oil palm (Elaeis guineensis Jacq.) during zygotic embryo development

BACKGROUND

Oil palm is the most efficient oil-producing crop in the world, and the yield of palm oil is associated with embryonic development. However, a comprehensive understanding of zygotic embryo development at the molecular level remains elusive. In order to address this issue, we report the transcriptomic analysis of zygotic embryo development in oil palm, specifically focusing on regulatory genes involved in important biological pathways.

RESULTS

In this study, three cDNA libraries were prepared from embryos at S1 (early-stage), S2 (middle-stage), and S3 (late-stage). There were 16,367, 16,500, and 18,012 genes characterized at the S1, S2, and S3 stages of embryonic development, respectively. A total of 1522, 2698, and 142 genes were differentially expressed in S1 vs S2, S1 vs S3, and S2 vs S3, respectively. Using Gene Ontology (GO) term enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis to identify key genes and pathways. In the hormone signaling pathway, genes related to auxin antagonize the output of cytokinin which regulates the development of embryo meristem. The genes related to abscisic acid negatively regulating the synthesis of gibberellin were strongly up-regulated in the mid-late stage of embryonic development. The results were reported the early synthesis and mid-late degradation of sucrose, as well as the activation of the continuous degradation pathway of temporary starch, providing the nutrients needed for differentiation of the embryonic cell. Moreover, the transcripts of genes involved in fatty acid synthesis were also abundantly accumulated in the zygotic embryos.

CONCLUSION

Taken together, our research provides a new perspective on the developmental and metabolic regulation of zygotic embryo development at the transcriptional level in oil palm.

PAs Regulate Early Somatic Embryo Development by Changing the Gene Expression Level and the Hormonal Balance in Dimocarpus longan Lour

Polyamines (PAs) play an important regulatory role in many basic cellular processes and physiological and biochemical processes. However, there are few studies on the identification of PA biosynthesis and metabolism family members and the role of PAs in the transition of plant embryogenic calli (EC) into globular embryos (GE), especially in perennial woody plants. We identified 20 genes involved in PA biosynthesis and metabolism from the third-generation genome of longan (Dimocarpus longan Lour.). There were no significant differences between longan and other species regarding the number of members, and they had high similarity with Citrus sinensis. Light, plant hormones and a variety of stress cis-acting elements were found in these family members. The biosynthesis and metabolism of PAs in longan were mainly completed by DlADC2, DlSAMDC2, DlSAMDC3, DlSPDS1A, DlSPMS, DlCuAOB, DlCuAO3A, DlPAO2 and DlPAO4B. In addition, 0.01 mmol∙L⁻¹ 1-aminocyclopropane-1-carboxylic acid (ACC), putrescine (Put) and spermine (Spm), could promote the transformation of EC into GE, and Spm treatment had the best effect, while 0.01 mmol∙L⁻¹ D-arginine (D-arg) treatment inhibited the process. The period between the 9th and 11th days was key for the transformation of EC into GE in longan. There were higher levels of gibberellin (GA), salicylic acid (SA) and abscisic acid (ABA) and lower levels of indole-3-acetic acid (IAA), ethylene and hydrogen peroxide (H₂O₂) in this key period. The expression levels in this period of DlADC2, DlODC, DlSPDS1A, DlCuAOB and DlPAO4B were upregulated, while those of DlSAMDC2 and DlSPMS were downregulated. These results showed that the exogenous ACC, D-arg and PAs could regulate the transformation of EC into GE in longan by changing the content of endogenous hormones and the expression levels of PA biosynthesis and metabolism genes. This study provided a foundation for further determining the physicochemical properties and molecular evolution characteristics of the PA biosynthesis and metabolism gene families, and explored the mechanism of PAs and ethylene for regulating the transformation of plant EC into GE.

Stress responses and epigenomic instability mark the loss of somatic embryogenesis competence in grapevine

Somatic embryogenesis (SE) represents the most appropriate tool for next-generation breeding methods in woody plants such as grapevine (Vitis vinifera L.). However, in this species, the SE competence is strongly genotype-dependent and the molecular basis of this phenomenon is poorly understood. We explored the genetic and epigenetic basis of SE in grapevine by profiling the transcriptome, epigenome, and small RNAome of undifferentiated, embryogenic, and non-embryogenic callus tissues derived from two genotypes differing in competence for SE, Sangiovese and Cabernet Sauvignon. During the successful formation of embryonic callus, we observed the upregulation of epigenetic-related transcripts and short interfering RNAs in association with DNA hypermethylation at transposable elements in both varieties. Nevertheless, the switch to nonembryonic development matched the incomplete reinforcement of transposon silencing, and the evidence of such effect was more apparent in the recalcitrant Cabernet Sauvignon. Transcriptomic differences between the two genotypes were maximized already at early stage of culture where the recalcitrant variety expressed a broad panel of genes related to stress responses and secondary metabolism. Our data provide a different angle on the SE molecular dynamics that can be exploited to leverage SE as a biotechnological tool for fruit crop breeding.

Somatic Embryogenesis of Norway Spruce and Scots Pine: Possibility of Application in Modern Forestry

Somatic embryogenesis (SE) is an important method for the vegetative propagation of trees. SE is the developmental in vitro process in which embryos are produced from somatic cells. This method can be integrated with other biotechnological techniques, genomic breeding and cryopreservation, which enables commercial-scale sapling production of selected high-yielding genotypes in wood production combined with fast breeding cycles. The SE is potential tool to improve plant stock in comparison with seed orchards. It can be useful for ecologically and economically important species, such as Norway spruce (Picea abies L. Karst.) and Scots pine (Pinus sylvestris L.), ensuring stable production in the era of climate change and biodiversity crisis. In this review, we summarize the current state of research on problems associated with somatic embryogenesis in P. abies and P. sylvestris.

Transduction of Signals during Somatic Embryogenesis

Somatic embryogenesis (SE) is an in vitro biological process in which bipolar structures (somatic embryos) can be induced to form from somatic cells and regenerate into whole plants. Acquisition of the embryogenic potential in culture is initiated when some competent cells within the explants respond to inductive signals (mostly plant growth regulators, PRGs), and de-differentiate into embryogenic cells. Such cells, "canalized" into the embryogenic developmental pathway, are able to generate embryos comparable in structure and physiology to their in vivo counterparts. Genomic and transcriptomic studies have identified several pathways governing the initial stages of the embryogenic process. In this review, the authors emphasize the importance of the developmental signals required for the progression of embryo development, starting with the de-differentiation of somatic cells and culminating with tissue patterning during the formation of the embryo body. The action and interaction of PGRs are highlighted, along with the participation of master regulators, mostly transcription factors (TFs), and proteins involved in stress responses and the signal transduction required for the initiation of the embryogenic process.

Synthetic Strigolactone GR24 Improves Arabidopsis Somatic Embryogenesis through Changes in Auxin Responses

Somatic embryogenesis in Arabidopsis encompasses an induction phase requiring auxin as the inductive signal to promote cellular dedifferentiation and formation of the embryogenic tissue, and a developmental phase favoring the maturation of the embryos. Strigolactones (SLs) have been categorized as a novel group of plant hormones based on their ability to affect physiological phenomena in plants. The study analyzed the effects of synthetic strigolactone GR24, applied during the induction phase, on auxin response and formation of somatic embryos. The expression level of two SL biosynthetic genes, MOREAXILLARY GROWTH 3 and 4 (MAX3 and MAX4), which are responsible for the conversion of carotene to carotenal, increased during the induction phase of embryogenesis. Arabidopsis mutant studies indicated that the somatic embryo number was inhibited in max3 and max4 mutants, and this effect was reversed by applications of GR24, a synthetic strigolactone, and exacerbated by TIS108, a SL biosynthetic inhibitor. The transcriptional studies revealed that the regulation of GR24 and TIS108 on somatic embryogenesis correlated with changes in expression of AUXIN RESPONSIVE FACTORs 5, 8, 10, and 16, known to be required for the production of the embryogenic tissue, as well as the expression of WUSCHEL (WUS) and Somatic Embryogenesis Receptor-like Kinase 1 (SERK1), which are markers of cell dedifferentiation and embryogenic tissue formation. Collectively, this work demonstrated the novel role of SL in enhancing the embryogenic process in Arabidopsis and its requirement for inducing the expression of genes related to auxin signaling and production of embryogenic tissue.

The transcriptional dynamics during de novo shoot organogenesis of Ma bamboo (Dendrocalamus latiflorus Munro): implication of the contributions of the abiotic stress response in this process

De novo shoot organogenesis is an important biotechnological tool for fundamental studies in plant. However, it is difficult in most bamboo species, and the genetic control of this highly dynamic and complicated regeneration process remains unclear. In this study, based on an in-depth analysis at the cellular level, the shoot organogenesis from calli of Ma bamboo (Dendrocalamus latiflorus Munro) was divided into five stages. Subsequently, single-molecule long-read isoform sequencing of tissue samples pooled from all five stages was performed to generate a full-length transcript landscape. A total of 83 971 transcripts, including 73 209 high-quality full-length transcripts, were captured, which served as an annotation reference for the subsequent RNA sequencing analysis. Time-course transcriptome analysis of samples at the abovementioned five stages was conducted to investigate the global gene expression atlas showing genome-wide expression of transcripts during the course of bamboo shoot organogenesis. K-means clustering analysis and stage-specific transcript identification revealed important dynamically expressed transcription regulators that function in bamboo shoot organogenesis. The majority of abiotic stress-responsive genes altered their expression levels during this process, and further experiments demonstrated that exogenous application of moderate but not severe abiotic stress increased the shoot regeneration efficiency. In summary, our study provides an overview of the genetic flow dynamics during bamboo shoot organogenesis. Full-length cDNA sequences generated in this study can serve as a valuable resource for fundamental and applied research in bamboo in the future.

Genetic and Molecular Control of Somatic Embryogenesis

Somatic embryogenesis is a method of asexual reproduction that can occur naturally in various plant species and is widely used for clonal propagation, transformation and regeneration of different crops. Somatic embryogenesis shares some developmental and physiological similarities with zygotic embryogenesis as it involves common actors of hormonal, transcriptional, developmental and epigenetic controls. Here, we provide an overview of the main signaling pathways involved in the induction and regulation of somatic embryogenesis with a focus on the master regulators of seed development, LEAFY COTYLEDON 1 and 2, ABSCISIC ACID INSENSITIVE 3 and FUSCA 3 transcription factors whose precise role during both zygotic and somatic embryogenesis remains to be fully elucidated.

Leaf transcriptomic signatures for somatic embryogenesis potential of Elaeis guineensis

Potentially embryogenic oil palms can be identified through leaf transcriptomic signatures. Differential expression of genes involved in flowering time, and stress and light responses may associate with somatic embryogenesis potential. Clonal propagation is an attractive approach for the mass propagation of high yielding oil palms. A major issue hampering the effectiveness of oil palm tissue culture is the low somatic embryogenesis rate. Previous studies have identified numerous genes involved in oil palm somatic embryogenesis, but their association with embryogenic potential has not been determined. In this study, differential expression analysis of leaf transcriptomes from embryogenic and non-embryogenic mother palms revealed that transcriptome profiles from non- and poor embryogenic mother palms were more similar than highly embryogenic palms. A total of 171 genes exhibiting differential expression in non- and low embryogenesis groups could also discriminate high from poor embryogenesis groups of another tissue culture agency. Genes related to flowering time or transition such as FTIP, FRIGIDA-LIKE, and NF-YA were up-regulated in embryogenic ortets, suggesting that reproduction timing of the plant may associate with somatic embryogenesis potential. Several light response or photosynthesis-related genes were down-regulated in embryogenic ortets, suggesting a link between photosynthesis activity and embryogenic potential. As expression profiles of the differentially expressed genes are very similar between non- and low embryogenic groups, machine learning approaches with several candidate genes may generate a more sensitive model to better discriminate non-embryogenic from embryogenic ortets.

Similarities and Differences in the GFP Movement in the Zygotic and Somatic Embryos of Arabidopsis

Intercellular signaling during embryo patterning is not well understood and the role of symplasmic communication has been poorly considered. The correlation between the symplasmic domains and the development of the embryo organs/tissues during zygotic embryogenesis has only been described for a few examples, including Arabidopsis. How this process occurs during the development of somatic embryos (SEs) is still unknown. The aim of these studies was to answer the question: do SEs have a restriction in symplasmic transport depending on the developmental stage that is similar to their zygotic counterparts? The studies included an analysis of the GFP distribution pattern as expressed under diverse promoters in zygotic embryos (ZEs) and SEs. The results of the GFP distribution in the ZEs and SEs showed that 1/the symplasmic domains between the embryo organs and tissues in the SEs was similar to those in the ZEs and 2/the restriction in symplasmic transport in the SEs was correlated with the developmental stage and was similar to the one in their zygotic counterparts, however, with the spatio-temporal differences and different PDs SEL value between these two types of embryos.

Transcriptomes analysis reveals novel insight into the molecular mechanisms of somatic embryogenesis in Hevea brasiliensis

BACKGROUND

Somatic embryogenesis (SE) is a promising technology for plant vegetative propagation, which has an important role in tree breeding. Though rubber tree (Hevea brasiliensis Muell. Arg.) SE has been founded, few late SE-related genes have been identified and the molecular regulation mechanisms of late SE are still not well understood.

RESULTS

In this study, the transcriptomes of embryogenic callus (EC), primary embryo (PE), cotyledonary embryo (CE), abnormal embryo (AE), mature cotyledonary embryo (MCE) and withered abnormal embryo (WAE) were analyzed. A total of 887,852,416 clean reads were generated, 85.92% of them were mapped to the rubber tree genome. The de novo assembly generated 36,937 unigenes. The differentially expressed genes (DEGs) were identified in the pairwise comparisons of CE vs. AE and MCE vs. WAE, respectively. The specific common DEGs were mainly involved in the phytohormones signaling pathway, biosynthesis of phenylpropanoid and starch and sucrose metabolism. Among them, hormone signal transduction related genes were significantly enriched, especially the auxin signaling factors (AUX-like1, GH3.1, SAUR32-like, IAA9-like, IAA14-like, IAA27-like, IAA28-like and ARF5-like). The transcription factors including WRKY40, WRKY70, MYBS3-like, MYB1R1-like, AIL6 and bHLH93-like were characterized as molecular markers for rubber tree late SE. CML13, CML36, CAM-7, SERK1 and LEAD-29-like were also related to rubber tree late SE. In addition, histone modification had crucial roles during rubber tree late SE.

CONCLUSIONS

This study provides important information to elucidate the molecular regulation during rubber tree late SE.

Shotgun label-free proteomic and biochemical study of somatic embryos (cotyledonary and maturation stage) in Catharanthus roseus (L.) G. Don

Somatic embryogenesis is an important and wonderful biotechnological tool used to develop whole plant from a single or a group of somatic cells. The differentiated somatic cells become totipotent stem cells by drastic reprogramming of a wide range of cellular activities, leading to the acquisition of embryogenic competence. After acquiring competence, the cells pass through globular, heart, torpedo and cotyledonary stages of embryo; however, all advanced embryos do not convert into full plant, produce adventive embryos or callus instead, thus reverses the programming. This is a big limitation in propagation of many plants. Understanding and unraveling the proteins at this 'embryo to plantlet' transition stage will help to get more numbers of plants. Thus, our study was aimed at an identification of differentially abundant proteins between two important advanced stages, i.e. cotyledonary-(T1) and maturation stage (T2) of somatic embryos in Catharanthus roseus. A total of 2949 and 3030 proteins were identified in cotyledonary and maturation stage, respectively. Of these, 1129 proteins were common to both. Several proteins were found to be differentially accumulated in two different embryo stages in which over 60 proteins were most accumulated during somatic embryo maturation time. More chlorophyll accumulation was noted at this time under the influence of gibberellic acid (GA3). Proteins like Mg-protoporphyrin IX chelatase, chlorophyll a-b-binding protein, photosystem I iron-sulfur center, photosystem II Psb, photosystem II subunit P-1, P-II domain-containing protein, RuBisCO large chain, RuBisCO small chain, RuBisCO activase, RuBisCO large subunit-binding proteins were synthesized. Some of the identified proteins are linked to chlorophyll synthesis, carbohydrate metabolism and stress. The identified proteins are categorized into different groups on the basis of their cellular location, role and other metabolic processes. Biochemical attributes like protein, sugar, proline, antioxidant enzyme (APX, SOD and CAT) activities were high in T2 as compared to T1. The proteins like peroxidases, pathogenesis-related proteins, the late-embryogenesis abundant proteins, argonaute, germin and others have been discussed in C. roseus somatic embryo maturation process.

Taking the Wheel - de novo DNA Methylation as a Driving Force of Plant Embryonic Development

During plant embryogenesis, regardless of whether it begins with a fertilized egg cell (zygotic embryogenesis) or an induced somatic cell (somatic embryogenesis), significant epigenetic reprogramming occurs with the purpose of parental or vegetative transcript silencing and establishment of a next-generation epigenetic patterning. To ensure genome stability of a developing embryo, large-scale transposon silencing occurs by an RNA-directed DNA methylation (RdDM) pathway, which introduces methylation patterns de novo and as such potentially serves as a global mechanism of transcription control during developmental transitions. RdDM is controlled by a two-armed mechanism based around the activity of two RNA polymerases. While PolIV produces siRNAs accompanied by protein complexes comprising the methylation machinery, PolV produces lncRNA which guides the methylation machinery toward specific genomic locations. Recently, RdDM has been proposed as a dominant methylation mechanism during gamete formation and early embryo development in Arabidopsis thaliana, overshadowing all other methylation mechanisms. Here, we bring an overview of current knowledge about different roles of DNA methylation with emphasis on RdDM during plant zygotic and somatic embryogenesis. Based on published chromatin immunoprecipitation data on PolV binding sites within the A. thaliana genome, we uncover groups of auxin metabolism, reproductive development and embryogenesis-related genes, and discuss possible roles of RdDM at the onset of early embryonic development via targeted methylation at sites involved in different embryogenesis-related developmental mechanisms.

Changes in abscisic acid metabolism in relation to the maturation of grapevine (Vitis vinifera L., cv. Mencía) somatic embryos

BACKGROUND

Somatic embryogenesis in grapevines is a complex process that depends on many physiological and genetic factors. One of its main limitations is the process of precocious germination of the somatic embryos in differentiation medium. This process lowers plant conversion rates from the somatic embryos, and it is probably caused by a low endogenous abscisic acid (ABA) content.

RESULTS

Precocious germination of the somatic embryos was successfully avoided by culturing grapevine cv. Mencía embryogenic aggregates over a semipermeable membrane extended on top of the differentiation medium. The weekly analysis of the endogenous ABA and ABA-glucosyl ester (ABA-GE) contents in the aggregates showed their rapid accumulation. The expression profiles of 9-cis-epoxycarotenoid dioxygenase (VvNCED1), 8'-hydroxylase (VvHyd2), UDP-glucosyltransferase (VvUGT) and β-glucosidase (VvBG2) genes in grapevine revealed that the occurrence of a first accumulation peak of endogenous ABA in the second week of culture over the semipermeable membrane was mainly dependent on the expression of the VvNCED1 gene. A second increase in the endogenous ABA content was observed in the fourth week of culture. At this point in the culture, our results suggest that of those genes involved in ABA accumulation, one (VvNCED1) was repressed, while another (VvBG2) was activated. Similarly, of those genes related to a reduction in ABA levels, one (VvUGT) was repressed while another (VvHyd2) was activated. The relative expression level of the VvNCED1 gene in embryogenic aggregates cultured under the same conditions and treated with exogenous ABA revealed the significant downregulation of this gene.

CONCLUSIONS

Our results demonstrated the involvement of ABA metabolism in the control of the maturation of grapevine somatic embryos cultured over a semipermeable membrane and two important control points for their endogenous ABA levels. Thus, subtle differences in the expression of the antagonistic genes that control ABA synthesis and degradation could be responsible for the final level of ABA during the maturation of grapevine somatic embryos in vitro. In addition, the treatment of somatic embryos with exogenous ABA suggested the feedback-based control of the expression of the VvNCED1 gene by ABA during the maturation of grapevine somatic embryos.

Genetic activity during early plant embryogenesis

Seeds are essential for human civilization, so understanding the molecular events underpinning seed development and the zygotic embryo it contains is important. In addition, the approach of somatic embryogenesis is a critical propagation and regeneration strategy to increase desirable genotypes, to develop new genetically modified plants to meet agricultural challenges, and at a basic science level, to test gene function. We briefly review some of the transcription factors (TFs) involved in establishing primary and apical meristems during zygotic embryogenesis, as well as TFs necessary and/or sufficient to drive somatic embryo programs. We focus on the model plant Arabidopsis for which many tools are available, and review as well as speculate about comparisons and contrasts between zygotic and somatic embryo processes.

Heterologous Expression of Poplar WRKY18/35 Paralogs in Arabidopsis Reveals Their Antagonistic Regulation on Pathogen Resistance and Abiotic Stress Tolerance via Variable Hormonal Pathways

WRKY transcription factors (WRKY TFs) are one of the largest protein families in plants, and most of them play vital roles in response to biotic and abiotic stresses by regulating related signaling pathways. In this study, we isolated two WRKY TF genes PtrWRKY18 and PtrWRKY35 from Populustrichocarpa and overexpressed them in Arabidopsis. Expression pattern analyses showed that PtrWRKY18 and PtrWRKY35 respond to salicylic acid (SA), methyl JA (MeJA), abscisic acid (ABA), B. cinereal, and P. syringae treatment. The transgenic plants conferred higher B. cinerea tolerance than wild-type (WT) plants, and real-time quantitative (qRT)-PCR assays showed that PR3 and PDF1.2 had higher expression levels in transgenic plants, which was consistent with their tolerance to B. cinereal. The transgenic plants showed lower P. syringae tolerance than WT plants, and qRT-PCR analysis (PR1, PR2, and NPR1) also corresponded to this phenotype. Germination rate and root analysis showed that the transgenic plants are less sensitive to ABA, which leads to the reduced tolerance to osmotic stress and the increase of the death ratio and stomatal aperture. Compared with WT plants, a series of ABA-related genes (RD29A, ABO3, ABI4, ABI5, and DREB1A) were significantly down-regulated in PtrWRKY18 and PtrWRKY35 overexpression plants. All of these results demonstrated that the two WRKY TFs are multifunctional transcription factors in plant resistance.

Research Tools for the Functional Genomics of Plant miRNAs During Zygotic and Somatic Embryogenesis

During early plant embryogenesis, some of the most fundamental decisions on fate and identity are taken making it a fascinating process to study. It is no surprise that higher plant embryogenesis was intensively analysed during the last century, while somatic embryogenesis is probably the most studied regeneration model. Encoded by the MIRNA, short, single-stranded, non-coding miRNAs, are commonly present in all Eukaryotic genomes and are involved in the regulation of the gene expression during the essential developmental processes such as plant morphogenesis, hormone signaling, and developmental phase transition. During the last few years dedicated to miRNAs, analytical methods and tools have been developed, which have afforded new opportunities in functional analyses of plant miRNAs, including (i) databases for in silico analysis; (ii) miRNAs detection and expression approaches; (iii) reporter and sensor lines for a spatio-temporal analysis of the miRNA-target interactions; (iv) in situ hybridisation protocols; (v) artificial miRNAs; (vi) MIM and STTM lines to inhibit miRNA activity, and (vii) the target genes resistant to miRNA. Here, we attempted to summarise the toolbox for functional analysis of miRNAs during plant embryogenesis. In addition to characterising the described tools/methods, examples of the applications have been presented.

Dynamic Transcriptome Analysis Reveals Uncharacterized Complex Regulatory Pathway Underlying Genotype-Recalcitrant Somatic Embryogenesis Transdifferentiation in Cotton

As a notable illustration of totipotency and plant regeneration, somatic embryogenesis (SE) is the developmental reprogramming of somatic cells toward the embryogenesis pathway, the key step for genetic engineering. Investigations examining the totipotency process are of great fundamental and practical importance in crop biotechnology. However, high-frequency regeneration of cotton via SE has been limited due to genotype-dependent response. The molecular basis deciphering SE genotype recalcitrance remains largely unexplored in cotton. In the current study, to comprehensively investigate the dynamic transcriptional profiling and gene regulatory patterns involved in SE process, a genome-wide RNA sequencing analysis was performed in two cotton genotypes with distinct embryogenic abilities, the highly embryogenic genotype Yuzao 1 (YZ) and the recalcitrant genotype Lumian 1 (LM). Three typical developmental staged cultures of early SE-hypocotyls (HY), nonembryogenic calli (NEC) and primary embryogenic calli (PEC)-were selected to establish the transcriptional profiles. Our data revealed that a total of 62,562 transcripts were present amongst different developmental stages in the two genotypes. Of these, 18,394 and 26,514 differentially expressed genes (DEGs) were identified during callus dedifferentiation (NEC-VS-HY) and embryogenic transdifferentiation (PEC-VS-NEC), respectively in the recalcitrant genotype, 21,842 and 22,343 DEGs in the highly embryogenic genotype. Furthermore, DEGs were clustered into six expression patterns during cotton SE process in the two genotypes. Moreover, functional enrichment analysis revealed that DEGs were significantly enriched in fatty acid, tryptophan and pyruvate metabolism in the highly embryogenic genotype and in DNA conformation change otherwise in the recalcitrant genotype. In addition, critical SE-associated expressed transcription factors, as well as alternative splicing events, were notably and preferentially activated during embryogenic transdifferentiation in the highly embryogenic genotype compared with the recalcitrant genotype. Taken together, by systematically comparing two genotypes with distinct embryogenic abilities, the findings in our study revealed a comprehensive overview of the dynamic gene regulatory patterns and uncharacterized complex regulatory pathways during cotton SE genotype-dependent response. Our work provides insights into the molecular basis and important gene resources for understanding the underlying genotype recalcitrance during SE process and plant regeneration, thereby holding great promise for accelerating the application of biotechnology to cotton for improving its breeding efficiency.

Genome-wide identification of miRNAs and their targets during early somatic embryogenesis in Dimocarpus longan Lour

miRNAs are endogenous regulatory factors that play pivotal roles in post-transcriptional regulation. However, their specific roles in early somatic embryogenesis (SE) remain unclear. Study of the SE system is fundamental for clarifying the molecular mechanisms in Dimocarpus longan. We identified 289 known miRNAs from 106 different miRNA families and 1087 novel miRNAs during early longan SE, including embryogenic callus (EC), incomplete pro-embryogenic culture (ICpEC), globular embryo (GE), and non-embryogenic callus (NEC). The abundances of known miRNAs were concentrated in GE. The differentially expression (DE) miRNAs showed five expression patterns during early SE. Largely miRNAs were expressed highly and specially in EC, ICpEC, and GE, respectively. Some miRNAs and putative target genes were enriched in lignin metabolism. Most potential targets were related to the pathways of plant hormone signal transduction, alternative splicing, tyrosine metabolism and sulfur metabolism in early longan SE. The regulatory relationships between dlo-miR166a-3p and DlHD-zip8, dlo-miR397a and DlLAC7, dlo-miR408-3p and DlLAC12 were confirmed by RNA ligase-mediated rapid amplification of cDNA ends. The expression patterns of eight DE miRNAs detected by qRT-PCR were consistent with RNA-seq. Finally, the miRNA regulatory network in early SE was constructed, which provided new insight into molecular mechanism of early SE in longan.

Characterization of a novel cotton MYB gene, GhMYB108-like responsive to abiotic stresses

Transcriptional factors are the major regulators of plant signaling pathways in response to environmental stresses i.e., drought, salinity and cold. Hereby, the GhMYB108-like was characterized to determine whether it regulate these stresses. The GhMYB108-like cDNA consisted of 1107 base pairs (bp) with 807 open reading frame encoded a protein of 268 amino acids. Its isoelectric point and molecular weight are 5.51 and 30.3 kDa respectively. Phylogenetic analysis and online databases revealed that GhMYB108-like proteins are closely related with the Arabidopsis thaliana MYB2. Important cis-elements were detected in the promotor region of GhMYB108-like responding to stresses and phytohormones. The 3D structure of GhMYB108-like protein has been predicted. In addition, various physico-chemical properties of GhMYB108-like have been determined. Subcellular localization confirmed that GhMYB108-like are nuclear localized protein. Quantitative expression analysis showed that polyethylene glycol and salt treatments significantly induced the expression of GhMYB108-like. Overall, our findings suggest that GhMYB108-like is an important gene that would plays important regulatory role in response to drought and salt stresses.

Full-Length Transcriptome Analysis of the ABCB, PIN/PIN-LIKES, and AUX/LAX Families Involved in Somatic Embryogenesis of Lilium pumilum DC. Fisch

Plant cell totipotency is one of the 25 major topics in current scientific research, and somatic embryos are good experimental material for studying cell totipotency. Polar auxin transport plays an important regulatory role in somatic embryogenesis (SE). However, little is known about the auxin transport genes and their regulatory mechanisms in Lilium SE. In this study, we applied single-molecule real-time (SMRT) sequencing to Lilium pumilum DC. Fisch. for the first time and obtained a total of 119,649 transcripts, of which 14 encoded auxin transport genes. Correlation analyses between somatic embryo induction and gene expression under different treatments revealed that auxin transport genes, especially ATP-binding cassette (ABC) transporter B family member 21 (ABCB21) and PIN-FORMED (PIN) LIKES 7 (PILS7), may be key players in SE, and the necessary duration of picloram (PIC) treatment to induce SE is as short as 3 days. Our research provides valuable genetic information on Lilium pumilum, elucidating the candidate auxin transport genes involved in SE and their influencing factors. This study lays a foundation for elucidating the regulatory mechanism of auxin transport in SE.

Global scale transcriptome analysis reveals differentially expressed genes involve in early somatic embryogenesis in Dimocarpus longan Lour

BACKGROUND

Somatic embryogenesis (SE) is a process of somatic cells that dedifferentiate to totipotent embryonic stem cells and generate embryos in vitro. Longan SE has been established and wildly used as model system for studying embryogenesis in woody plants, SE-related genes had been characterized. In spite of that, a comprehensive overview of SE at a molecular level is still absent. To understand the molecular mechanisms during longan SE, we examined the transcriptome changes by using Illumina HiSeq from the four distinct developmental stages, including non-embryogenic callus (NEC), embryogenic callus (EC), incomplete compact pro-embryogenic cultures (ICpEC), globular embryos (GE).

RESULTS

RNA-seq of the four samples generated a total of 243.78 million high quality reads, approximately 81.5% of the data were mapped to longan genome. The cDNA libraries of NEC, EC, ICpEC and GE, generated 22,743, 19,745, 21,144, 21,102 expressed transcripts, 1935, 1710, 1816, 1732 novel transcripts, 2645, 366, 505, 588 unique genes, respectively. Comparative transcriptome analysis showed that a total of 10,642, 4180, 5846 and 1785 genes were differentially expressed in the pairwise comparisons of NEC_vs_EC, EC_vs_ICpEC, EC_vs_GE, ICpEC_vs_GE, respectively. Among them, plant hormones signalling related genes were significantly enriched, especially the auxin and cytokinin signalling components. The transcripts of flavonoid biosynthesis related genes were mainly expressed in NEC, while fatty acid biosynthesis related genes mainly accumulated in early SE. In addition, the extracelluar protein encoding genes LTP, CHI, GLP, AGP, EP1 were related to longan SE. Combined with the FPKM value of longan nine tissues transcription, 27 SE specific or preferential genes (LEC1, LEC1-like, PDF1.3, GH3.6, AGL80, PIN1, BBM, WOX9, WOX2, ABI3, et al.) and 28 NEC preferential genes (LEA5, CNOT3, DC2.15, PR1-1, NsLTP2, DIR1, PIP1, PIP2.1, TIP2-1, POD-P7 and POD5 et al.) were characterized as molecular markers for longan early SE. qRT-PCR validation of SE-related genes showed a high correlation between RNA-seq and qRT-PCR data.

CONCLUSION

This study provides new insights into the role of the transcriptome during early SE in longan. Differentially expressed genes reveal that plant hormones signalling, flavonoid and fatty acid biosynthesis, and extracelluar protein related genes were involved in longan early SE. It could serve as a valuable platform resource for further functional studies addressing embryogenesis in woody plants.

Insights into Drought Stress Signaling in Plants and the Molecular Genetic Basis of Cotton Drought Tolerance

Drought stress restricts plant growth and development by altering metabolic activity and biological functions. However, plants have evolved several cellular and molecular mechanisms to overcome drought stress. Drought tolerance is a multiplex trait involving the activation of signaling mechanisms and differentially expressed molecular responses. Broadly, drought tolerance comprises two steps: stress sensing/signaling and activation of various parallel stress responses (including physiological, molecular, and biochemical mechanisms) in plants. At the cellular level, drought induces oxidative stress by overproduction of reactive oxygen species (ROS), ultimately causing the cell membrane to rupture and stimulating various stress signaling pathways (ROS, mitogen-activated-protein-kinase, Ca2+, and hormone-mediated signaling). Drought-induced transcription factors activation and abscisic acid concentration co-ordinate the stress signaling and responses in cotton. The key responses against drought stress, are root development, stomatal closure, photosynthesis, hormone production, and ROS scavenging. The genetic basis, quantitative trait loci and genes of cotton drought tolerance are presented as examples of genetic resources in plants. Sustainable genetic improvements could be achieved through functional genomic approaches and genome modification techniques such as the CRISPR/Cas9 system aid the characterization of genes, sorted out from stress-related candidate single nucleotide polymorphisms, quantitative trait loci, and genes. Exploration of the genetic basis for superior candidate genes linked to stress physiology can be facilitated by integrated functional genomic approaches. We propose a third-generation sequencing approach coupled with genome-wide studies and functional genomic tools, including a comparative sequenced data (transcriptomics, proteomics, and epigenomic) analysis, which offer a platform to identify and characterize novel genes. This will provide information for better understanding the complex stress cellular biology of plants.

Genome-Wide Analysis of Cotton Auxin Early Response Gene Families and Their Roles in Somatic Embryogenesis

Auxin is well known to regulate growth and development processes. Auxin early response genes serve as a critical component of auxin signaling and mediate auxin regulation of diverse physiological processes. In the present study, a genome-wide identification and comprehensive analysis of auxin early response genes were conducted in upland cotton. A total of 71 auxin response factor (ARF), 86 Auxin/Indole-3-Acetic Acid (Aux/IAA), 63 Gretchen Hagen3 (GH3), and 194 small auxin upregulated RNA (SAUR) genes were identified in upland cotton, respectively. Phylogenetic analysis revealed that the ARF, GH3, and SAUR families were likely subject to extensive evolutionary divergence between Arabidopsis and upland cotton, while the Aux/IAA family was evolutionary conserved. Expression profiles showed that the ARF, Aux/IAA, GH3, and SAUR family genes were extensively involved in embryogenic competence acquisition of upland cotton callus. The Aux/IAA family genes generally showed a higher expression level in the non-embryogenic callus (NEC) of highly embryogenic cultivar CCRI24 than that of recalcitrant cultivar CCRI12, which may be conducive to initializing the embryogenic transformation. Auxin early response genes were tightly co-expressed with most of the known somatic embryogenesis (SE) related genes, indicating that these genes may regulate upland cotton SE by interacting with auxin early response genes.

Metabolome and Transcriptome Association Analysis Reveals Dynamic Regulation of Purine Metabolism and Flavonoid Synthesis in Transdifferentiation during Somatic Embryogenesis in Cotton

Plant regeneration via somatic embryogenesis (SE) is a key step during genetic engineering. In the current study, integrated widely targeted metabolomics and RNA sequencing were performed to investigate the dynamic metabolic and transcriptional profiling of cotton SE. Our data revealed that a total of 581 metabolites were present in nonembryogenic staged calli (NEC), primary embryogenic calli (PEC), and initiation staged globular embryos (GE). Of the differentially accumulated metabolites (DAMs), nucleotides, and lipids were specifically accumulated during embryogenic differentiation, whereas flavones and hydroxycinnamoyl derivatives were accumulated during somatic embryo development. Additionally, metabolites related to purine metabolism were significantly enriched in PEC vs. NEC, whereas in GE vs. PEC, DAMs were remarkably associated with flavonoid biosynthesis. An association analysis of the metabolome and transcriptome data indicated that purine metabolism and flavonoid biosynthesis were co-mapped based on the Kyoto encyclopedia of genes and genomes (KEGG) database. Moreover, purine metabolism-related genes associated with signal recognition, transcription, stress, and lipid binding were significantly upregulated. Moreover, several classic somatic embryogenesis (SE) genes were highly correlated with their corresponding metabolites that were involved in purine metabolism and flavonoid biosynthesis. The current study identified a series of potential metabolites and corresponding genes responsible for SE transdifferentiation, which provides a valuable foundation for a deeper understanding of the regulatory mechanisms underlying cell totipotency at the molecular and biochemical levels.

Signaling Overview of Plant Somatic Embryogenesis

Somatic embryogenesis (SE) is a means by which plants can regenerate bipolar structures from a somatic cell. During the process of cell differentiation, the explant responds to endogenous stimuli, which trigger the induction of a signaling response and, consequently, modify the gene program of the cell. SE is probably the most studied plant regeneration model, but to date it is the least understood due to the unclear mechanisms that occur at a cellular level. In this review, the authors seek to emphasize the importance of signaling on plant SE, highlighting the interactions between the different plant growth regulators (PGR), mainly auxins, cytokinins (CKs), ethylene and abscisic acid (ABA), during the induction of SE. The role of signaling is examined from the start of cell differentiation through the early steps on the embryogenic pathway, as well as its relation to a plant's tolerance of different types of stress. Furthermore, the role of genes encoded to transcription factors (TFs) during the embryogenic process such as the LEAFY COTYLEDON (LEC), WUSCHEL (WUS), BABY BOOM (BBM) and CLAVATA (CLV) genes, Arabinogalactan-proteins (AGPs), APETALA 2 (AP2) and epigenetic factors is discussed.

Somatic Embryogenesis Induction in Woody Species: The Future After OMICs Data Assessment

Very early somatic embryogenesis has been recognized as a powerful method to propagate plants in vitro. For some woody species and in particular for some coniferous trees, somatic embryogenesis induction has become a routine procedure. For the majority, the application of this technology presents yet many limitations especially due to the genotype, the induction conditions, the number of embryos produced, maturation, and conversion, among other factors that compromise the systematic use of somatic embryogenesis for commercial purposes especially of woody species and forest trees in particular. The advancements obtained on somatic embryogenesis in Arabidopsis and the development of OMIC technologies allowed the characterization of genes and the corresponding proteins that are conserved in woody species. This knowledge will help in understanding the molecular mechanisms underlying the complex regulatory networks that control somatic embryogenesis in woody plants. In this revision, we report on developments of OMICs (genomics, transcriptomics, metabolomics, and proteomics) applied to somatic embryogenesis induction and its contribution for understanding the change of fate giving rise to the expression of somatic embryogenesis competence.

GhL1L1 affects cell fate specification by regulating GhPIN1-mediated auxin distribution

Auxin is as an efficient initiator and regulator of cell fate during somatic embryogenesis (SE), but the molecular mechanisms and regulating networks of this process are not well understood. In this report, we analysed SE process induced by Leafy cotyledon1-like 1 (GhL1L1), a NF-YB subfamily gene specifically expressed in embryonic tissues in cotton. We also identified the target gene of GhL1L1, and its role in auxin distribution and cell fate specification during embryonic development was analysed. Overexpression of GhL1L1 accelerated embryonic cell formation, associated with an increased concentration of IAA in embryogenic calluses (ECs) and in the shoot apical meristem, corresponding to altered expression of the auxin transport gene GhPIN1. By contrast, GhL1L1-deficient explants showed retarded embryonic cell formation, and the concentration of IAA was decreased in GhL1L1-deficient ECs. Disruption of auxin distribution accelerated the specification of embryonic cell fate together with regulation of GhPIN1. Furthermore, we showed that PHOSPHATASE 2AA2 (GhPP2AA2) was activated by GhL1L1 through targeting the G-box of its promoter, hence regulating the activity of GhPIN1 protein. Our results indicate that GhL1L1 functions as a key regulator in auxin distribution to regulate cell fate specification in cotton and contribute to the understanding of the complex process of SE in plant species.

Temperature and Water Availability During Maturation Affect the Cytokinins and Auxins Profile of Radiata Pine Somatic Embryos

Somatic embryogenesis (SE) provides us a potent biotechnological tool to manipulate the physical and chemical conditions (water availability) along the process and to study their effect in the final success in terms of quantity of somatic embryos produced. In the last years, our research team has been focused on the study of different aspects of the SE in Pinus spp. One of the main aspects affecting SE is the composition of culture media; in this sense, phytohormones play one of the most crucial roles in this propagation system. Many studies in conifers have shown that different stages of SE and somatic embryo development are correlated with distinct endogenous phytohormone profiles under the stress conditions needed for the process (i.e., cytokinins play a regulatory role in stress signaling, which it is essential for radiata pine SE). Based on this knowledge, the aim of this study was to test the effect of different temperatures (18, 23, and 28°C) and gelling agent concentrations (8, 9, and 10 gL-1) during the maturation stage of Pinus radiata SE in maturation and germination rates. Parallel, phytohormone profile of somatic embryos developed was evaluated. In this sense, the highest gellan gum concentration led to significantly lower water availability. At this gellan gum concentration and 23°C a significantly higher number of somatic embryos was obtained and the overall success of the process increased with respect to other treatments assayed. The somatic embryos produced in these conditions showed the highest concentration of iP-type cytokinins and total ribosides. Although, the different conditions applied during maturation of somatic embryos led to different hormonal profiles, they did not affect the ex vitro survival of the resulting somatic plants, where no significant differences were observed.

Physical state of the culture medium triggers shift in morphogenetic pattern from shoot bud formation to somatic embryo in Solanum khasianum

Solanum khasianum is a rich source of steroidal alkaloids that are important secondary metabolites with enormous pharmaceutical uses. Development of plantlets from somatic tissues, under in vitro conditions, takes place both through adventitious shoot bud differentiation or somatic embryogenesis (SE) pathway. We observed that the physical state of medium, solid or liquid, determined the regenerant differentiation patterns from root segment explants in S. khasianum. In the solidified medium, the root segments developed adventitious shoot buds whereas somatic embryos were regenerated in the liquid medium. Varying gradients from liquid to solid medium were further used to confirm the effect of solidified condition on regeneration pathway. Histological analysis of developing shoot buds and somatic embryos was also performed to confirm their development and differentiation patterns. In order to further confirm the developmental pathways, qRT-PCR analysis of the marker genes of SE and shoot regeneration was also performed. While SOMATIC EMBRYOGENESIS RECEPTOR KINASE1 (SkSERK1) expression was significantly up-regulated during the early embryogenic stage, the LATE EMBRYOGENESIS ABUNDANT (SkLEA) protein was found to be highly expressed in the mature embryos. Expression of the HISTONE DEACETYLASE (HDA6), a repressor of SE related genes, was highly decreased during embryogenesis in the liquid culture. Furthermore, expression of the ENHANCER OF SHOOT REGENERATION  (ESR) gene was comparatively increased during shoot regeneration in the culture using solid medium. Our results point out that the physical state of the medium in S. khasianum plays a decisive role in differentiation pattern which was independent of hormonal supplements.

Comparative transcriptome analysis of the wild-type model apomict Hieracium praealtum and its loss of parthenogenesis (lop) mutant

BACKGROUND

Asexual seed formation (apomixis) has been observed in diverse plant families but is rare in crop plants. The generation of apomictic crops would revolutionize agriculture, as clonal seed production provides a low cost and efficient way to produce hybrid seed. Hieracium (Asteraceae) is a model system for studying the molecular components of gametophytic apomixis (asexual seed reproduction).

RESULTS

In this study, a reference transcriptome was produced from apomictic Hieracium undergoing the key apomictic events of apomeiosis, parthenogenesis and autonomous endosperm development. In addition, transcriptome sequences from pre-pollination and post-pollination stages were generated from a loss of parthenogenesis (lop) mutant accession that exhibits loss of parthenogenesis and autonomous endosperm development. The transcriptome is composed of 147,632 contigs, 50% of which were annotated with orthologous genes and their probable function. The transcriptome was used to identify transcripts differentially expressed during apomictic and pollination dependent (lop) seed development. Gene Ontology enrichment analysis of differentially expressed transcripts showed that an important difference between apomictic and pollination dependent seed development was the expression of genes relating to epigenetic gene regulation. Genes that mark key developmental stages, i.e. aposporous embryo sac development and seed development, were also identified through their enhanced expression at those stages.

CONCLUSION

The production of a comprehensive floral reference transcriptome for Hieracium provides a valuable resource for research into the molecular basis of apomixis and the identification of the genes underlying the LOP locus.

Jasmonic Acid-Mediated Aliphatic Glucosinolate Metabolism Is Involved in Clubroot Disease Development in Brassica napus L

Glucosinolate (GSL) is associated with clubroot disease, which is caused by the obligate biotrophic protist Plasmodiophora brassicae. Due to the complicated composition of GSLs, their exact role in clubroot disease development remains unclear. By investigating clubroot disease resistance in cruciferous plants and characterizing the GSL content in seeds, we can determine if clubroot disease development is related to the components of GSLs. The difference in the infection process between Matthiola incana L. (resistant) and Brassica napus L. (susceptible) was determined. Root hair infection was definitely observed in both resistant and susceptible hosts, but no infection was observed during the cortical infection stage in resistant roots; this finding was verified by molecular detection of P. brassicae via PCR amplification at various times after inoculation. Based on the time course detection of the contents and compositions of GSLs after P. brassicae inoculation, susceptible roots exhibited increased accumulation of aliphatic, indolic, and aromatic GSLs in B. napus, but only aromatic GSLs were significantly increased in M. incana. Gluconapin, which was the main aliphatic GSL in B. napus and present only in B. napus, was significantly increased during the secondary infection stage. Quantification of the internal jasmonic acid (JA) concentration showed that both resistant and susceptible plants exhibited an enhanced level of JA, particularly in susceptible roots. The exogenous JA treatment induced aliphatic GSLs in B. napus and aromatic GSLs in M. incana. JA-induced aromatic GSLs may be involved in the defense against P. brassicae, whereas aliphatic GSLs induced by JA in B. napus likely play a role during the secondary infection stage. Three candidate MYB28 genes regulate the content of aliphatic GSLs identified in B. napus; one such gene was BnMYB28.1, which was significantly increased following both the treatment with exogenous JA and P. brassicae inoculation. In summary, the increased content of JA during the secondary infection stage may induce the expression of BnMYB28.1, which caused the accumulation of aliphatic GSLs in clubroot disease development.

miR156-SPL modules regulate induction of somatic embryogenesis in citrus callus

miR156 is a highly conserved plant miRNA and has been extensively studied because of its versatile roles in plant development. Here, we report a novel role of miR156 in regulating somatic embryogenesis (SE) in citrus, one of the most widely cultivated fruit crops in the world. SE is an important means of in vitro regeneration, but over the course of long-term sub-culturing there is always a decline in the SE potential of the preserved citrus embryogenic callus, and this represents a key obstacle for citrus biotechnology. In this study, the SE competence of citrus callus of wild kumquat (Fortunella hindsii) was significantly enhanced by either overexpression of csi-miR156a or by individual knock-down of the two target genes, CsSPL3 and CsSPL14, indicating that the effect of miR156-SPL modules was established during the initial phases of SE induction. Biological processes that might promote SE in response to miR156 overexpression were explored using RNA-seq, and mainly included hormone signaling pathways, stress responses, DNA methylation, and the cell cycle. CsAKIN10 was identified as interacting protein of CsSPL14. Our results provide insights into the regulatory pathway through which miR156-SPL modules enhance the SE potential of citrus callus, and provide a theoretical basis for improvement of plant SE competence.