De novo transcriptome analysis reveals insights into dynamic homeostasis regulation of somatic embryogenesis in upland cotton (G. hirsutum L.)

Physiological, Metabolic, and Transcriptomic Analyses Reveal Mechanisms of Proliferation and Somatic Embryogenesis of Litchi (Litchi chinensis Sonn.) Embryogenic Callus Promoted by D-Arginine Treatment

D-arginine (D-Arg) can promote embryogenic callus (EC) proliferation and increase the rate of somatic embryo induction of litchi (Litchi chinensis Sonn.), yet the mechanism underlying the processes is incompletely understood. To investigate the mechanism, physiological responses of polyamines (PAs) [putrescine (Put), spermidine (Spd), and spermine (Spm)] were investigated for D-Arg-treated litchi EC and enzyme activity related to polyamine metabolism, plant endogenous hormones, and polyamine- and embryogenic-related genes were explored. Results showed that the exogenous addition of D-Arg reduces the activity of diamine oxidase (DAO) and polyamine oxidase (PAO) in EC, reduces the production of H2O2, promotes EC proliferation, and increases the (Spd + Spm)/Put ratio to promote somatic embryo induction. Exogenous D-Arg application promoted somatic embryogenesis (SE) by increasing indole-3-acetyl glycine (IAA-Gly), kinetin-9-glucoside (K9G), and dihydrozeatin-7-glucoside (DHZ7G) levels and decreasing trans-zeatin riboside (tZR), N-[(-)-jasmonoyl]-(L)-valine (JA-Val), jasmonic acid (JA), and jasmonoyl-L-isoleucine (Ja-ILE) levels on 18 d, as well as promoting cell division and differentiation. The application of exogenous D-Arg regulated EC proliferation and somatic embryo induction by altering gene expression levels of the WRKY family, AP2/ERF family, C3H family, and C2H2 family. These results indicate that exogenous D-Arg could regulate the proliferation of EC and the SE induction of litchi by changing the biosynthesis of PAs through the alteration of gene expression pattern and endogenous hormone metabolism.

Transcriptome Analysis of Embryogenic and Non-Embryogenic Callus of Picea Mongolica

Picea mongolica is a rare tree species in China, which is of great significance in combating desertification and improving the harsh ecological environment. Due to the low rate of natural regeneration, high mortality, and susceptibility to pests and cold springs, Picea mongolica has gradually become extinct. At present, somatic embryogenesis (SE) is the most effective method of micro-proliferation in conifers, but the induction rate of embryogenic callus (EC) is low, and EC is difficult to differentiate from non-embryonic callus (NEC). Therefore, the EC and NEC of Picea mongolica were compared from the morphology, histological, physiological, and transcriptional levels, respectively. Morphological observation showed that the EC was white and transparent filamentous, while the NEC was compact and brownish-brown lumpy. Histological analyses showed that the NEC cells were large and loosely arranged; the nuclei attached to the edge of the cells were small; the cytoplasm was low; and the cell gap was large and irregular. In the EC, small cells, closely arranged cells, and a large nucleus and nucleolus were observed. Physiological studies showed significant differences in ROS-scavenging enzymes between the EC and NEC. Transcriptome profiling revealed that 13,267 differentially expressed genes (DEGs) were identified, 3682 were up-regulated, and 9585 were down-regulated. In total, 63 GO terms had significant enrichment, 32 DEGs in plant hormone signal transduction pathway were identified, and 502 different transcription factors (TFs) were characterized into 38 TF families. Meanwhile, we identified significant gene expression trends associated with somatic embryo development in plant hormones (AUX/IAA, YUCCA, LEA, etc.), stress (GST, HSP, GLP, etc.), phenylpropanoid metabolism (4CL, HCT, PAL, etc.), and transcription factors (AP2/ERF, MYB, WOX, etc.). In addition, nine genes were chosen for RT-qPCR, and the results were consistent with RNA-Seq data. This study revealed the changes in morphology, histology, physiology, and gene expression in the differentiation of NEC into EC and laid the foundation for finding the key genes to promote EC formation.

Endogenous auxin maintains embryonic cell identity and promotes somatic embryo development in Arabidopsis

Somatic embryogenesis (SE), or embryo development from in vitro cultured vegetative explants, can be induced in Arabidopsis by the synthetic auxin 2,4-dichlorophenoxyacetic acid (2,4-D) or by overexpression of specific transcription factors, such as AT-HOOK MOTIF NUCLEAR LOCALIZED 15 (AHL15). Here, we explored the role of endogenous auxin [indole-3-acetic acid (IAA)] during 2,4-D and AHL15-induced SE. Using the pWOX2:NLS-YFP reporter, we identified three distinct developmental stages for 2,4-D and AHL15-induced SE in Arabidopsis, with these being (i) acquisition of embryo identity; (ii) formation of pro-embryos; and (iii) somatic embryo patterning and development. The acquisition of embryo identity coincided with enhanced expression of the indole-3-pyruvic acid auxin biosynthesis YUCCA genes, resulting in an enhanced pDR5:GFP-reported auxin response in the embryo-forming tissues. Chemical inhibition of the indole-3-pyruvic acid pathway did not affect the acquisition of embryo identity, but significantly reduced or completely inhibited the formation of pro-embryos. Co-application of IAA with auxin biosynthesis inhibitors in the AHL15-induced SE system rescued differentiated somatic embryo formation, confirming that increased IAA levels are important during the last two stages of SE. Our analyses also showed that polar auxin transport, with AUXIN/LIKE-AUX influx and PIN-FORMED1 efflux carriers as important drivers, is required for the transition of embryonic cells to proembryos and, later, for correct cell fate specification and differentiation. Taken together, our results indicate that endogenous IAA biosynthesis and its polar transport are not required for the acquisition of embryo identity, but rather to maintain embryonic cell identity and for the formation of multicellular proembryos and their development into histodifferentiated embryos.

Dynamic Transcriptome Analysis Reveals Complex Regulatory Pathway Underlying Induction and Dose Effect by Different Exogenous Auxin IAA and 2,4-D During in vitro Embryogenic Redifferentiation in Cotton

Plant somatic cells can reprogram into differentiated embryos through somatic embryogenesis (SE) on the condition of plant growth regulators (PGRs). RNA sequencing analysis was performed to investigate transcriptional profiling on cotton redifferentiated callus that was induced by different auxin types (IAA and 2,4-D), different concentrations (0, 0.025, and 0.05 mg L^-1), and different incubation times (0, 5, and 20 days). Under the 2,4-D induction effect, signal transduction pathways of plant hormones were significantly enriched in the embryogenic response stage (5 days). These results indicated that auxin signal transduction genes were necessary for the initial response of embryogenic differentiation. In the pre-embryonic initial period (20 days), the photosynthetic pathway was significantly enriched. Most differentially expressed genes (DEGs) were downregulated under the induction of 2,4-D. Upon the dose effect of IAA and 2,4-D, respectively, pathways were significantly enriched in phenylpropanoid biosynthesis, fatty acid metabolism, and carbon metabolic pathways. Therefore, primary and secondary metabolism pathways were critical in cotton SE. These results showed that complex synergistic mechanisms involving multiple cellular pathways were the causes of the induction and dose effect of auxin-induced SE. This study reveals a systematic molecular response to auxin signals and reveals the way that regulates embryogenic redifferentiation during cotton SE.

One-Week Scutellar Somatic Embryogenesis in the Monocot Brachypodium distachyon

Plant somatic embryogenesis (SE) is a natural process of vegetative propagation. It can be induced in tissue cultures to investigate developmental transitions, to create transgenic or edited lines, or to multiply valuable crops. We studied the induction of SE in the scutellum of monocots with Brachypodium distachyon as a model system. Towards the in-depth analysis of SE initiation, we determined the earliest stages at which somatic scutellar cells acquired an embryogenic fate, then switched to a morphogenetic mode in a regeneration sequence involving treatments with exogenous hormones: first an auxin (2,4-D) then a cytokinin (kinetin). Our observations indicated that secondary somatic embryos could already develop in the proliferative calli derived from immature zygotic embryo tissues within one week from the start of in vitro culture. Cell states and tissue identity were deduced from detailed histological examination, and in situ hybridization was performed to map the expression of key developmental genes. The fast SE induction method we describe here facilitates the mechanistic study of the processes involved and may significantly shorten the production of transgenic or gene-edited plants.

Transcriptional Landscape of Cotton Fiber Development and Its Alliance With Fiber-Associated Traits

Cotton fiber development is still an intriguing question to understand fiber commitment and development. At different fiber developmental stages, many genes change their expression pattern and have a pivotal role in fiber quality and yield. Recently, numerous studies have been conducted for transcriptional regulation of fiber, and raw data were deposited to the public repository for comprehensive integrative analysis. Here, we remapped > 380 cotton RNAseq data with uniform mapping strategies that span ∼400 fold coverage to the genome. We identified stage-specific features related to fiber cell commitment, initiation, elongation, and Secondary Cell Wall (SCW) synthesis and their putative cis-regulatory elements for the specific regulation in fiber development. We also mined Exclusively Expressed Transcripts (EETs) that were positively selected during cotton fiber evolution and domestication. Furthermore, the expression of EETs was validated in 100 cotton genotypes through the nCounter assay and correlated with different fiber-related traits. Thus, our data mining study reveals several important features related to cotton fiber development and improvement, which were consolidated in the "CottonExpress-omics" database.

Embryogenic Calli Induction and Salt Stress Response Revealed by RNA-Seq in Diploid Wild Species Gossypium sturtianum and Gossypium raimondii

Wild cotton species can contribute to a valuable gene pool for genetic improvement, such as genes related to salt tolerance. However, reproductive isolation of different species poses an obstacle to produce hybrids through conventional breeding. Protoplast fusion technology for somatic cell hybridization provides an opportunity for genetic manipulation and targeting of agronomic traits. Transcriptome sequencing analysis of callus under salt stress is conducive to study salt tolerance genes. In this study, calli were induced to provide materials for extracting protoplasts and also for screening salt tolerance genes. Calli were successfully induced from leaves of Gossypium sturtianum (C₁ genome) and hypocotyls of G. raimondii (D₅ genome), and embryogenic calli of G. sturtianum and G. raimondii were induced on a differentiation medium with different concentrations of 2, 4-D, KT, and IBA, respectively. In addition, embryogenic calli were also induced successfully from G. raimondii through suspension cultivation. Transcriptome sequencing analysis was performed on the calli of G. raimondii and G. sturtianum, which were treated with 200 mM NaCl at 0, 6, 12, 24, and 48 h, and a total of 12,524 genes were detected with different expression patterns under salt stress. Functional analysis showed that 3,482 genes, which were differentially expressed in calli of G. raimondii and G. sturtianum, were associated with biological processes of nucleic acid binding, plant hormone (such as ABA) biosynthesis, and signal transduction. We demonstrated that DEGs or TFs which related to ABA metabolism were involved in the response to salt stress, including xanthoxin dehydrogenase genes (ABA2), sucrose non-fermenting 1-related protein kinases (SnRK2), NAM, ATAT1/2, and CUC2 transcription factors (NAC), and WRKY class of zinc-finger proteins (WRKY). This research has successfully induced calli from two diploid cotton species and revealed new genes responding to salt stress in callus tissue, which will lay the foundation for protoplast fusion for further understanding of salt stress responses in cotton.

miRNAs as key regulators via targeting the phytohormone signaling pathways during somatic embryogenesis of plants

Somatic embryogenesis is the regeneration of embryos from the somatic cell via dedifferentiation and redifferentiation without the occurrence of fertilization. A complex network of genes regulates the somatic embryogenesis process. Especially, microRNAs (miRNAs) have emerged as key regulators by affecting phytohormone biosynthesis, transport and signal transduction pathways. miRNAs are small, non-coding small RNA regulatory molecules involved in various developmental processes including somatic embryogenesis. Several types of miRNAs such as miR156, miR157, miR 159, miR 160, miR165, miR166, miR167, miR390, miR393 and miR396 have been reported to intricate in regulating somatic embryogenesis via targeting the phytohormone signaling pathways. Here we review current research progress on the miRNA-mediated regulation involved in somatic embryogenesis via regulating auxin, ethylene, abscisic acid and cytokinin signaling pathways. Further, we also discussed the possible role of other phytohormone signaling pathways such as gibberellins, jasmonates, nitric oxide, polyamines and brassinosteroids. Finally, we conclude by discussing the expression of miRNAs and their targets involved in somatic embryogenesis and possible regulatory mechanisms cross talk with phytohormones during somatic embryogenesis.

Transcriptome analyses provide insights into the homeostatic regulation of axillary buds in upland cotton (G. hirsutum L.)

BACKGROUND

The axillary bud is an important index of cotton plant-type traits, and the molecular mechanism of axillary bud development in upland cotton has not yet been reported. We obtained a mutant (designated mZ571) with a high-budding phenotype in axillary bud development from the low-budding phenotype variety G. hirsutum Z571 (CCRI 9A02), which provided ideal materials for the study of complex regulatory networks of axillary bud development. In this study, RNA sequencing was carried out to detect gene expression levels during three stages of axillary buds in Z571 (LB, low budding) and mZ571 mutant (HB, high budding).

RESULTS

A total of 7162 DEGs were identified in the three groups (HB-E vs. LB-E, HB-G1 vs. LB-G1, HB-G2 vs. LB-G2), including 4014 downregulated and 3184 upregulated DEGs. Additionally, 221 DEGs were commonly identified in all three groups, accounting for approximately 3.09% of the total DEGs. These DEGs were identified, annotated and classified. A significant number of DEGs were related to hormone metabolism, hormone signal transduction, and starch and sucrose metabolism. In addition, 45, 22 and 9 DEGs involved in hormone metabolic pathways and 67, 22 and 19 DEGs involved in hormone signal transduction pathwayspathway were identified in HB-E vs. LB-E, HB-G1 vs. LB-G1, and HB-G2 vs. LB-G2, respectively, suggesting that endogenous hormones are the primary factors influencing cotton axillary bud growth. Hormone and soluble sugar content measurements revealed that mZ571 exhibited higher concentrations of zeatin, gibberellins and soluble sugar in all three stages, which confirmed that these hormone metabolism-, hormone signal transduction- and starch metabolism-related genes showed interaction effects contributing to the divergence of axillary bud growth between mZ571 and Z571.

CONCLUSIONS

Our results confirmed the importance of endogenous hormones and sugars in the development of axillary buds, and we found that mZ571 plants, with a high-budding phenotype of axillary buds, exhibited higher endogenous hormone and sugar concentrations. Overall, we present a model for the emergence and development of cotton axillary buds that provides insights into the complexity and dynamic nature of the regulatory network during axillary bud emergence and development.

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.

Current Perspectives on the Auxin-Mediated Genetic Network that Controls the Induction of Somatic Embryogenesis in Plants

Auxin contributes to almost every aspect of plant development and metabolism as well as the transport and signalling of auxin-shaped plant growth and morphogenesis in response to endo- and exogenous signals including stress conditions. Consistently with the common belief that auxin is a central trigger of developmental changes in plants, the auxin treatment of explants was reported to be an indispensable inducer of somatic embryogenesis (SE) in a large number of plant species. Treating in vitro-cultured tissue with auxins (primarily 2,4-dichlorophenoxyacetic acid, which is a synthetic auxin-like plant growth regulator) results in the extensive reprogramming of the somatic cell transcriptome, which involves the modulation of numerous SE-associated transcription factor genes (TFs). A number of SE-modulated TFs that control auxin metabolism and signalling have been identified, and conversely, the regulators of the auxin-signalling pathway seem to control the SE-involved TFs. In turn, the different expression of the genes encoding the core components of the auxin-signalling pathway, the AUXIN/INDOLE-3-ACETIC ACIDs (Aux/IAAs) and AUXIN RESPONSE FACTORs (ARFs), was demonstrated to accompany SE induction. Thus, the extensive crosstalk between the hormones, in particular, auxin and the TFs, was revealed to play a central role in the SE-regulatory network. Accordingly, LEAFY COTYLEDON (LEC1 and LEC2), BABY BOOM (BBM), AGAMOUS-LIKE15 (AGL15) and WUSCHEL (WUS) were found to constitute the central part of the complex regulatory network that directs the somatic plant cell towards embryogenic development in response to auxin. The revealing picture shows a high degree of complexity of the regulatory relationships between the TFs of the SE-regulatory network, which involve direct and indirect interactions and regulatory feedback loops. This review examines the recent advances in studies on the auxin-controlled genetic network, which is involved in the mechanism of SE induction and focuses on the complex regulatory relationships between the down- and up-stream targets of the SE-regulatory TFs. In particular, the outcomes from investigations on Arabidopsis, which became a model plant in research on genetic control of SE, are presented.

Dynamic Transcriptome Analysis Reveals Uncharacterized Complex Regulatory Pathway Underlying Dose IBA-Induced Embryogenic Redifferentiation in Cotton

The somatic embryogenesis (SE) process of plants is regulated by exogenous hormones. During the SE, different genes sensitively respond to hormone signals through complex regulatory networks to exhibit plant totipotency. When cultured in indole-3-butyric acid (IBA) concentration gradient medium supplemented with 0 mg dm^-3, 0.025 mg dm^-3, and 0.05 mg dm^-3 IBA, the callus differentiation rate first increased then decreased in cotton. To characterize the molecular basis of IBA-induced regulating SE, transcriptome analysis was conducted on embryogenic redifferentiation. Upon the examination of the IBA's embryogenic inductive effect, it was revealed that pathways related to plant hormone signal transduction and alcohol degradation were significantly enriched in the embryogenic responsive stage (5 days). The photosynthesis, alcohol metabolism and cell cycle pathways were specifically regulated in the pre-embryonic initial period (20 days). Upon the effect of the IBA dose, in the embryogenic responsive stage (5 days), the metabolism of xenobiotics by the cytochrome P450 pathway and secondary metabolism pathways of steroid, flavonoid, and anthocyanin biosynthesis were significantly enriched. The phenylpropanoid, brassinosteroid, and anthocyanin biosynthesis pathways were specifically associated in the pre-embryonic initial period (20 days). At different developmental stages of embryogenic induction, photosynthesis, flavonoid biosynthesis, phenylpropanoid biosynthesis, mitogen-activated protein kinase (MAPK) signaling, xenobiotics metabolism by cytochrome P450, and brassinosteroid biosynthesis pathways were enriched at low a IBA concentration. Meanwhile, at high IBA concentration, the carbon metabolism, alcohol degradation, circadian rhythm and biosynthesis of amino acids pathways were significantly enriched. The results reveal that complex regulating pathways participate in the process of IBA-induced redifferentiation in cotton somatic embryogenesis. In addition, collections of potential essential signaling and regulatory genes responsible for dose IBA-induced efficient embryogenic redifferentiation were identified. Quantitative real-time PCR (qRT-PCR) was performed on the candidate genes with different expression patterns, and the results are basically consistent with the RNA-seq data. The results suggest that the complicated and concerted IBA-induced mechanisms involving multiple cellular pathways are responsible for dose-dependent plant growth regulator-induced SE. This report represents a systematic study and provides new insight into molecular signaling and regulatory basis underlying the process of dose IBA-induced embryogenic redifferentiation during SE.

Transcriptome analysis of the induction of somatic embryogenesis in Coffea canephora and the participation of ARF and Aux/IAA genes

Background

Somatic embryogenesis (SE) is a useful biotechnological tool to study the morpho-physiological, biochemical and molecular processes during the development of Coffea canephora. Plant growth regulators (PGR) play a key role during cell differentiation in SE. The Auxin-response-factor (ARF) and Auxin/Indole-3-acetic acid (Aux/IAA) are fundamental components involved in the signaling of the IAA. The IAA signaling pathway activates or represses the expression of genes responsive to auxins during the embryogenic transition of the somatic cells. The growing development of new generation sequencing technologies (NGS), as well as bioinformatics tools, has allowed us to broaden the landscape of SE study of various plant species and identify the genes directly involved.

Methods

Analysis of transcriptome expression profiles of the C. canephora genome and the identification of a particular set of differentially expressed genes (DEG) during SE are described in this study.

Results

A total of eight ARF and seven Aux/IAA differentially expressed genes were identified during the different stages of the SE induction process. The quantitative expression analysis showed that ARF18 and ARF5 genes are highly expressed after 21 days of the SE induction, while Aux/IAA7 and Aux/IAA12 genes are repressed.

Discussion

The results of this study allow a better understanding of the genes involved in the auxin signaling pathway as well as their expression profiles during the SE process.

Screening the Reference Genes for Quantitative Gene Expression by RT-qPCR During SE Initial Dedifferentiation in Four Gossypium hirsutum Cultivars that Have Different SE Capability

RNA sequencing (RNA-Seq)-based gene expression analysis is applicable to a wide range of biological purposes in various species. Reverse transcription quantitative PCR (RT-qPCR) is also used to assess target gene expression utilizing stably expressed reference genes as internal control under a given set of conditions. However, investigations of the reference genes for RT-qPCR normalization in the process of somatic embryogenesis (SE) initial dedifferentiation in Gossypium hirsutum are rarely reported. In this study, on the basis of our previous transcriptome data of three different induction stages during SE initial dedifferentiation process in four G. hirsutum cultivars that have different SE capability, 15 candidate genes were selected during SE initial dedifferentiation process, and their expression stability was evaluated by geNorm, NormFinder, and BestKeeper. The results indicated that the two genes of endonuclease 4 (ENDO4) and 18S ribosomal RNA (18S rRNA) showed stable expression in the four different G. hirsutum cultivars, endowing them to be appropriate reference genes during three induction stages in the four cotton cultivars. In addition, the stability and reliability of the two reference genes of ENDO4 and 18S rRNA were further verified by comparing the expressions of auxin-responsive protein 22 (AUX22) and ethylene-responsive transcription factor 17 (ERF17) between RT-qPCR results and the RNA-seq data, which showed strong positive correlation coefficient (R² = 0.8396–0.9984), validating again the steady expression of ENDO4 and 18S rRNA as the reliable reference genes. Our results provide effective reference genes for RT-qPCR normalization during SE process in different G. hirsutum cultivars.

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.

Dynamic TMT-Based Quantitative Proteomics Analysis of Critical Initiation Process of Totipotency during Cotton Somatic Embryogenesis Transdifferentiation

The somatic embryogenesis (SE) process of plants, as one of the typical responses to abiotic stresses with hormone, occurs through the dynamic expression of different proteins that constitute a complex regulatory network in biological activities and promotes plant totipotency. Plant SE includes two critical stages: primary embryogenic calli redifferentiation and somatic embryos development initiation, which leads to totipotency. The isobaric labels tandem mass tags (TMT) large-scale and quantitative proteomics technique was used to identify the dynamic protein expression changes in nonembryogenic calli (NEC), primary embryogenic calli (PEC) and globular embryos (GEs) of cotton. A total of 9369 proteins (6730 quantified) were identified; 805, 295 and 1242 differentially accumulated proteins (DAPs) were identified in PEC versus NEC, GEs versus PEC and GEs versus NEC, respectively. Eight hundred and five differentially abundant proteins were identified, 309 of which were upregulated and 496 down regulated in PEC compared with NEC. Of the 295 DAPs identified between GEs and PEC, 174 and 121 proteins were up- and down regulated, respectively. Of 1242 differentially abundant proteins, 584 and 658 proteins were up- and down regulated, respectively, in GEs versus NEC. We have also complemented the authenticity and accuracy of the proteomic analysis. Systematic analysis indicated that peroxidase, photosynthesis, environment stresses response processes, nitrogen metabolism, phytohormone response/signal transduction, transcription/posttranscription and modification were involved in somatic embryogenesis. The results generated in this study demonstrate a proteomic molecular basis and provide a valuable foundation for further investigation of the roles of DAPs in the process of SE transdifferentiation during cotton totipotency.

Transcriptome analysis of callus from melon

OBJECTIVE

To identify the key genes promoting the differentiation of melon non-embryogenic callus into embryogenic callus.

METHODS

The transcriptome sequencing analysis was used to analyze the mRNA sequence in embryogenic callus (Z) and non-embryogenic callus (F); transcript mapping, gene expression analysis, cluster analysis, classification analysis and enrichment analysis were then used to detect the differentially expressed genes and enriched pathways.

RESULTS

The correlation coefficient between sample Z and sample F was 0.929 after transcript mapping. The overall gene expression levels in sample Z were higher as compared with sample F. Furthermore, cluster analysis showed that the expression of genes involved in photosynthesis was increased in sample Z when comparing to F. Besides, the classification of differential Gene Ontology (GO) showed that many metabolic processes were affected with the metabolism enhanced in embryogenic callus. Interestingly, Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis further demonstrated the high metabolic activity and active secondary metabolite formation in the embryogenic callus.

CONCLUSION

The genes associated with photosynthesis, metabolic pathways and biosynthesis of secondary metabolites may promote the differentiation of callus into embryogenic callus.

Overexpression of the CsFUS3 gene encoding a B3 transcription factor promotes somatic embryogenesis in Citrus

In citrus, genetic improvement via biotechnology is challenging due to insufficient understanding of molecular barriers that prevent regeneration by somatic embryogenesis (SE). Our previous study indicated that LEC genes were involved in SE in citrus, but their regulatory roles remain to be elucidated. Here, we cloned one of the LEC genes, CsFUS3, and show that it is preferentially expressed during SE and in the embryogenic callus (EC) derived from citrus varieties with strong embryogenic competence. The overexpression of CsFUS3 in recalcitrant citrus callus restored embryogenic competence. Complementation of the loss-of-function Arabidopsis fus3 mutant with the CsFUS3 gene restored normal late embryogenesis, which is consistent with the CsFUS3 and AtFUS3 proteins contributing to the same regulatory network in Arabidopsis. Transcriptome profiling revealed that the expression of particular TFs that promote SE was up-regulated in the citrus overexpression (OE) line. The 104 differentially expressed genes associated with hormone biosynthesis, catabolism, and signaling are particularly noteworthy. The dynamic change in the ratio of ABA to GA during SE in wild-type callus mirrored the expression pattern of CsFUS3. In contrast, in the OE line, the ratio of ABA to GA was higher and the capacity for SE was greater when the OE line was separately treated with ABA and GA biosynthesis inhibitors. Taken together, our results demonstrate that the overexpression of CsFUS3 appears to establish a cellular environment favorable to SE, at least in part by promoting a high ABA to GA ratio and by regulating the expression of TFs that promote SE.

Trichostatin A Triggers an Embryogenic Transition in Arabidopsis Explants via an Auxin-Related Pathway

Auxin is an important regulator of plant ontogenies including embryo development and the exogenous application of this phytohormone has been found to be necessary for the induction of the embryogenic response in plant explants that have been cultured in vitro. However, in the present study, we show that treatment of Arabidopsis explants with trichostatin A (TSA), which is a chemical inhibitor of histone deacetylases, induces somatic embryogenesis (SE) without the exogenous application of auxin. We found that the TSA-treated explants generated somatic embryos that developed efficiently on the adaxial side of the cotyledons, which are the parts of an explant that are involved in auxin-induced SE. A substantial reduction in the activity of histone deacetylase (HDAC) was observed in the TSA-treated explants, thus confirming a histone acetylation-related mechanism of the TSA-promoted embryogenic response. Unexpectedly, the embryogenic effect of TSA was lower on the auxin-supplemented media and this finding further suggests an auxin-related mechanism of TSA-induced SE. Congruently, we found a significantly increased content of indolic compounds, which is indicative of IAA and an enhanced DR5::GUS signal in the TSA-treated explants. In line with these results, two of the YUCCA genes (YUC1 and YUC10), which are involved in auxin biosynthesis, were found to be distinctly up-regulated during TSA-induced SE and their expression was colocalised with the explant sites that are involved in SE. Beside auxin, ROS were extensively accumulated in response to TSA, thereby indicating that a stress-response is involved in TSA-triggered SE. Relevantly, we showed that the genes encoding the transcription factors (TFs) that have a regulatory function in auxin biosynthesis including LEC1, LEC2, BBM, and stress responses (MYB118) were highly up-regulated in the TSA-treated explants. Collectively, the results provide several pieces of evidence about the similarities between the molecular pathways of SE induction that are triggered by TSA and 2,4-D that involve the activation of the auxin-responsive TF genes that have a regulatory function in auxin biosynthesis and stress responses. The study suggests the involvement of histone acetylation in the auxin-mediated release of the embryogenic program of development in the somatic cells of Arabidopsis.

iTRAQ-based comparative proteomic analysis provides insights into somatic embryogenesis in Gossypium hirsutum L

KEY MESSAGE: iTRAQ based proteomic identified key proteins and provided new insights into the molecular mechanisms underlying somatic embryogenesis in cotton. Somatic embryogenesis, which involves cell dedifferentiation and redifferentiation, has been used as a model system for understanding molecular events of plant embryo development in vitro. In this study, we performed comparative proteomics analysis using samples of non-embryogenic callus (NEC), embryogenic callus (EC) and somatic embryo (SE) using the isobaric tags for relative and absolute quantitation (iTRAQ) technology. In total, 5892 proteins were identified amongst the three samples. The majority of these proteins (93.4%) were found to have catalytic activity, binding activity, transporter activity or structural molecular activity. Of these proteins, 1024 and 858 were differentially expressed in NEC versus EC and EC versus SE, respectively. Compared to NEC, EC had 452 and 572 down- and up-regulated proteins, respectively, and compared to EC, SE had 647 and 221 down- and up-regulated proteins, respectively. KEGG (Kyoto Encyclopedia of Genes and Genomes) analysis indicated that genetic information transmission, plant hormone transduction, glycolysis, fatty acid biosynthesis and metabolism, galactose metabolism were the top pathways involved in somatic embryogenesis. Our proteomics results not only confirmed our previous transcriptomic results on the role of the polyamine metabolic pathways and stress responses in cotton somatic embryogenesis, but identified key proteins important for cotton somatic embryogenesis and provided new insights into the molecular mechanisms underlying somatic embryogenesis in cotton.

Comparative Transcriptome Analysis of SE initial dedifferentiation in cotton of different SE capability

Somatic embryogenesis (SE) is a critical transition from vegetative to embryogenic growth in higher plants; however, few studies have investigated the mechanism that regulates SE initial differentiation. Most cotton varieties have not undergone regeneration by SE, so only a few varieties can be used in genetic engineering. Here, two varieties of cotton with different SE capabilities (HD, higher differentiation and LD, lower differentiation) were analyzed by high throughout RNA-Seq at the pre-induction stage (0h) and two induction stages (3h and 3d) under callus-induction medium (CIM). About 1150 million clean reads were obtained from 98.21% raw data. Transcriptomic analysis revealed that "protein kinase activity" and "oxidoreductase activity" were highly represented GO terms during the same and different treatment stages among HD and LD. Moreover, several stress-related transcription factors might play important roles in SE initiation. The SE-related regulation genes (SERKs) showed different expression patterns between HD and LD. Furthermore, the complex auxin and ethylene signaling pathway contributes to initiation of differentiation in SE. Thus, our RNA-sequencing of comparative transcriptome analysis will lay a foundation for future studies to better define early somatic formation in cotton with different SE capabilities.

Small RNA and Transcriptome Sequencing Reveal a Potential miRNA-Mediated Interaction Network That Functions during Somatic Embryogenesis in Lilium pumilum DC. Fisch

Plant somatic embryos are widely used in the fields of germplasm conservation, breeding for genetic engineering and artificial seed production. MicroRNAs (miRNAs) play pivotal roles in somatic embryogenesis (SE) regulation. However, their regulatory roles during various stages of SE remain unclear. In this study, six types of embryogenic samples of Lilium pumilum DC. Fisch., including organogenic callus, embryogenic callus induced for 4 weeks, embryogenic callus induced for 6 weeks, globular embryos, torpedo embryos and cotyledon embryos, were prepared for small RNA sequencing. The results revealed a total of 2,378,760 small RNA reads, among which the most common size was 24 nt. Four hundred and fifty-two known miRNAs, belonging to more than 86 families, 57 novel miRNAs and 40 miRNA*s were identified. The 86 known miRNA families were sorted according to an alignment with their homologs across 24 land plants into the following four categories: 23 highly conserved, 4 moderately conserved, 15 less conserved and 44 species-specific miRNAs. Differentially expressed known miRNAs were identified during various stages of SE. Subsequently, the expression levels of 12 differentially expressed miRNAs and 4 targets were validated using qRT-PCR. In addition, six samples were mixed in equal amounts for transcript sequencing, and the sequencing data were used as transcripts for miRNA target prediction. A total of 66,422 unigenes with an average length of 800 bp were assembled from 56,258,974 raw reads. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment indicated that 38,004 and 15,497 unigenes were successfully assigned to GO terms and KEGG pathways, respectively. Among the unigenes, 2,182 transcripts were predicted to be targets for 396 known miRNAs. The potential targets of the identified miRNAs were mostly classified into the following GO terms: cell, binding and metabolic process. Enriched KEGG analysis demonstrated that carbohydrate metabolism was the predominant pathway in Lilium SE. Thus, we performed systemic characterization, homology comparisons and profiling of miRNA expression, and we constructed an miRNA-target network during Lilium SE for the first time. Our findings establish a foundation for the further exploration of critical genes and elucidation of SE in Lilium.