Transcript profiling reveals complex auxin signalling pathway and transcription regulation involved in dedifferentiation and redifferentiation during somatic embryogenesis in cotton

Genetic Dissection of Maize Embryonic Callus Regenerative Capacity Using Multi-Locus Genome-Wide Association Studies

The regenerative capacity of the embryonic callus, a complex quantitative trait, is one of the main limiting factors for maize transformation. This trait was decomposed into five traits, namely, green callus rate (GCR), callus differentiating rate (CDR), callus plantlet number (CPN), callus rooting rate (CRR), and callus browning rate (CBR). To dissect the genetic foundation of maize transformation, in this study multi-locus genome-wide association studies (GWAS) for the five traits were performed in a population of 144 inbred lines genotyped with 43,427 SNPs. Using the phenotypic values in three environments and best linear unbiased prediction (BLUP) values, as a result, a total of 127, 56, 160, and 130 significant quantitative trait nucleotides (QTNs) were identified by mrMLM, FASTmrEMMA, ISIS EM-BLASSO, and pLARmEB, respectively. Of these QTNs, 63 QTNs were commonly detected, including 15 across multiple environments and 58 across multiple methods. Allele distribution analysis showed that the proportion of superior alleles for 36 QTNs was <50% in 31 elite inbred lines. Meanwhile, these superior alleles had obviously additive effect on the regenerative capacity. This indicates that the regenerative capacity-related traits can be improved by proper integration of the superior alleles using marker-assisted selection. Moreover, a total of 40 candidate genes were found based on these common QTNs. Some annotated genes were previously reported to relate with auxin transport, cell fate, seed germination, or embryo development, especially, GRMZM2G108933 (WOX2) was found to promote maize transgenic embryonic callus regeneration. These identified candidate genes will contribute to a further understanding of the genetic foundation of maize embryonic callus regeneration.

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.

Advanced Proteomic Approaches to Elucidate Somatic Embryogenesis

Somatic embryogenesis (SE) is a cell differentiation process by which a somatic cell changes its genetic program and develops into an embryonic cell. Investigating this process with various explant sources in vitro has allowed us to trace somatic embryo development from germination to plantlets and has led to the generation of new technologies, including genetic transformation, endangered species conservation, and synthetic seed production. A transcriptome data comparison from different stages of the developing somatic embryo has revealed a complex network controlling the somatic cell's fate, suggesting that an interconnected network acts at the protein level. Here, we discuss the current progress on SE using proteomic-based data, focusing on changing patterns of proteins during the establishment of the somatic embryo. Despite the advanced proteomic approaches available so far, deciphering how the somatic embryo is induced is still in its infancy. The new proteomics techniques that lead to the quantification of proteins with different abundances during the induction of SE are opening this area of study for the first time. These quantitative differences can elucidate the different pathways involved in SE induction. We envisage that the application of these proteomic technologies can be pivotal to identifying proteins critical to the process of SE, demonstrating the cellular localization, posttranslational modifications, and turnover protein events required to switch from a somatic cell to a somatic embryo cell and providing new insights into the molecular mechanisms underlying SE. This work will help to develop biotechnological strategies for mass production of quality crop material.

Elaboration of Transcriptome During the Induction of Somatic Embryogenesis

Somatic embryogenesis (SE) is one of the most studied developmental processes due to its applications, such as plant micropropagation, transformation, and germplasm conservation. The use of massive techniques of sequencing, as well as the use of subtractive hybridization and macroarrays, has led to the identification of hundreds of genes involved in the SE process. These have been important developments to study the molecular aspects of the progress of SE. With the advent of the new massive techniques for sequencing RNA, it has been possible to see a more complete picture of whole processes. In this chapter we present a technique to handle the elaboration of the transcriptome from the extraction of RNA until the assembly of the complete transcriptome.

De novo assembly and comparative analysis of the transcriptome of embryogenic callus formation in bread wheat (Triticum aestivum L.)

BACKGROUND

During asexual reproduction the embryogenic callus can differentiate into a new plantlet, offering great potential for fostering in vitro culture efficiency in plants. The immature embryos (IMEs) of wheat (Triticum aestivum L.) are more easily able to generate embryogenic callus than mature embryos (MEs). To understand the molecular process of embryogenic callus formation in wheat, de novo transcriptome sequencing was used to generate transcriptome sequences from calli derived from IMEs and MEs after 3d, 6d, or 15d of culture (DC).

RESULTS

In total, 155 million high quality paired-end reads were obtained from the 6 cDNA libraries. Our de novo assembly generated 142,221 unigenes, of which 59,976 (42.17%) were annotated with a significant Blastx against nr, Pfam, Swissprot, KOG, KEGG, GO and COG/KOG databases. Comparative transcriptome analysis indicated that a total of 5194 differentially expressed genes (DEGs) were identified in the comparisons of IME vs. ME at the three stages, including 3181, 2085 and 1468 DEGs at 3, 6 and 15 DC, respectively. Of them, 283 overlapped in all the three comparisons. Furthermore, 4731 DEGs were identified in the comparisons between stages in IMEs and MEs. Functional analysis revealed that 271transcription factor (TF) genes (10 overlapped in all 3 comparisons of IME vs. ME) and 346 somatic embryogenesis related genes (SSEGs; 35 overlapped in all 3 comparisons of IME vs. ME) were differentially expressed in at least one comparison of IME vs. ME. In addition, of the 283 overlapped DEGs in the 3 comparisons of IME vs. ME, excluding the SSEGs and TFs, 39 possessed a higher rate of involvement in biological processes relating to response to stimuli, in multi-organism processes, reproductive processes and reproduction. Furthermore, 7 were simultaneously differentially expressed in the 2 comparisons between the stages in IMEs, but not MEs, suggesting that they may be related to embryogenic callus formation. The expression levels of genes, which were validated by qRT-PCR, showed a high correlation with the RNA-seq value.

CONCLUSIONS

This study provides new insights into the role of the transcriptome in embryogenic callus formation in wheat, and will serve as a valuable resource for further studies addressing embryogenic callus formation in plants.

A transcriptional view on somatic embryogenesis

Somatic embryogenesis is a form of induced plant cell totipotency where embryos develop from somatic or vegetative cells in the absence of fertilization. Somatic embryogenesis can be induced in vitro by exposing explants to stress or growth regulator treatments. Molecular genetics studies have also shown that ectopic expression of specific embryo‐ and meristem‐expressed transcription factors or loss of certain chromatin‐modifying proteins induces spontaneous somatic embryogenesis. We begin this review with a general description of the major developmental events that define plant somatic embryogenesis and then focus on the transcriptional regulation of this process in the model plant Arabidopsis thaliana (arabidopsis). We describe the different somatic embryogenesis systems developed for arabidopsis and discuss the roles of transcription factors and chromatin modifications in this process. We describe how these somatic embryogenesis factors are interconnected and how their pathways converge at the level of hormones. Furthermore, the similarities between the developmental pathways in hormone‐ and transcription‐factor‐induced tissue culture systems are reviewed in the light of our recent findings on the somatic embryo‐inducing transcription factor BABY BOOM.

A genome-wide analysis of the small auxin-up RNA (SAUR) gene family in cotton

Background

Small auxin-up RNA (SAUR) gene family is the largest family of early auxin response genes in higher plants, which have been implicated in the regulation of multiple biological processes. However, no comprehensive analysis of SAUR genes has been reported in cotton (Gossypium spp.).

Results

In the study, we identified 145, 97, 214, and 176 SAUR homologous genes in the sequenced genomes of G. raimondii, G. arboreum, G. hirsutum, and G. barbadense, respectively. A phylogenetic analysis revealed that the SAUR genes can be classified into 10 groups. A further analysis of chromosomal locations and gene duplications showed that tandem duplication and segmental duplication events contributed to the expansion of the SAUR gene family in cotton. An exon-intron organization and motif analysis revealed the conservation of SAUR-specific domains, and the auxin responsive elements existed in most of the upstream sequences. The expression levels of 16 GhSAUR genes in response to an exogenous application of IAA were determined by a quantitative RT-PCR analysis. The genome-wide RNA-seq data and qRT-PCR analysis of selected SAUR genes in developing fibers revealed their differential expressions. The physical mapping showed that 20 SAUR genes were co-localized with fiber length quantitative trait locus (QTL) hotspots. Single nucleotide polymorphisms (SNPs) were detected for 12 of these 20 genes between G. hirsutum and G. barbadense, but no SNPs were identified between two backcross inbred lines with differing fiber lengths derived from a cross between the two cultivated tetraploids.

Conclusions

This study provides an important piece of genomic information for the SAUR genes in cotton and lays a solid foundation for elucidating the functions of SAUR genes in auxin signaling pathways to regulate cotton growth.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4224-2) contains supplementary material, which is available to authorized users.

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.

Endogenous small interfering RNAs associated with maize embryonic callus formation

The induction efficiency of maize embryonic callus is highly dependent on the genotype, and only a few lines possess a high capacity for callus formation. Although certain genes and pathways have been reported to contribute to the regulation of callus induction, to the best of our knowledge, the functions of the small interfering RNAs (siRNAs) involved in this process remain unknown. In this study, we identified 861 differentially expressed siRNAs and 576 target genes in the callus induction process. These target genes were classified into 3 clusters, and their functions involve controlling metalloexopeptidase activity, catalase activity, transcription regulation, and O-methyltransferase activity. In addition, certain genes related to auxin transport and stem cell or meristem development (e.g., PLT5-like, ARF15, SAUR-like, FAS1-like, Fea3, SCL5, and Zmwox2A) were regulated by the differentially expressed siRNAs. Moreover, zma-siR004119-2 directly cleaves the 5' UTR of Homeobox-transcription factor 25, which further leads to the down-regulation of its expression. Twelve 24-nt-siRNAs led to the hyper-methylation of GRMZM2G013465, which further decreases its expression. These results suggest that differentially expressed siRNAs regulate callus formation by controlling the expression of their target genes.

Auxin and nutritional stress coupled somatic embryogenesis in Oldenlandia umbellata L

Somatic embryos were induced from internodal segment derived callus of Oldenlandia umbellata L., in MS medium supplemented with different concentrations of 2,4-Dichlorophenoxy acetic acid (2,4-D). Initially calli were developed from internodes of microshoots inoculated in 2.5 µM NAA supplemented medium. Then calli were transferred to 2,4-D added medium for somatic embryogenesis. Nutritional stress coupled with higher concentration of 2,4-D triggered somatic embryogenesis. Nutritional stress was induced by culturing callus in a fixed amount of medium for a period up to 20 weeks without any external supply of nutrients. Addition of 2.5 µM 2,4-D gave 100% embryogenesis within 16 weeks of incubation. Callus mass bearing somatic embryos were transferred to germination medium facilitated production of in vitro plantlets. MS medium supplemented with 2.5 µM benzyl adenine and 0.5 µM α-naphthalene acetic acid produced 15.33 plants per culture within 4 weeks of culture. Somatic embryo germinated plants were then hardened and transferred to green house.

Differential proteome analysis during early somatic embryogenesis in Musa spp. AAA cv. Grand Naine

Endogenous hormone secretion proteins along with stress and defense proteins play predominant role in banana embryogenesis. This study reveals the underlying molecular mechanism during transition from vegetative to embryogenic state. Banana (Musa spp.) is well known globally as a food fruit crop for millions. The requirement of quality planting material of banana is enormous. Although mass multiplication through tissue culture is in vogue, high-throughput techniques like somatic embryogenesis (SE) as a mass multiplication tool needs to be improved. Apart from clonal propagation, SE has extensive applications in genetic improvement and mutation. SE in banana is completely genome-dependent and most of the commercial cultivars exhibit recalcitrance. Thus, understanding the molecular basis of embryogenesis in Musa will help to develop strategies for mass production of quality planting material. In this study, differentially expressed proteins between embryogenic calli (EC) and non-embryogenic calli (NEC) with respect to the explant, immature male flower buds (IMFB), of cv. Grand Naine (AAA) were determined using two-dimensional gel electrophoresis (2DE). The 2DE results were validated through qRT-PCR. In total, 65 proteins were identified: 42 were highly expressed and 23 were less expressed in EC compared to NEC and IMFB. qRT-PCR analysis of five candidate proteins, upregulated in EC, were well correlated with expression at transcript level. Further analysis of proteins showed that embryogenesis in banana is associated with the control of oxidative stress. The regulation of ROS scavenging system and protection of protein structure occurred in the presence of heat shock proteins. Alongside, high accumulation of stress-related cationic peroxidase and plant growth hormone-related proteins like indole-3-pyruvate monooxygenase and adenylate isopentenyltransferase in EC revealed the association with the induction of SE.

Transcript profiling for early stages during embryo development in Scots pine

BACKGROUND

Characterization of the expression and function of genes regulating embryo development in conifers is interesting from an evolutionary point of view. However, our knowledge about the regulation of embryo development in conifers is limited. During early embryo development in Pinus species the proembyo goes through a cleavage process, named cleavage polyembryony, giving rise to four embryos. One of these embryos develops to a dominant embryo, which will develop further into a mature, cotyledonary embryo, while the other embryos, the subordinate embryos, are degraded. The main goal of this study has been to identify processes that might be important for regulating the cleavage process and for the development of a dominant embryo.

RESULTS

RNA samples from embryos and megagametophytes at four early developmental stages during seed development in Pinus sylvestris were subjected to high-throughput sequencing. A total of 6.6 million raw reads was generated, resulting in 121,938 transcripts, out of which 36.106 contained ORFs. 18,638 transcripts were differentially expressed (DETs) in embryos and megagametophytes. GO enrichment analysis of transcripts up-regulated in embryos showed enrichment for different cellular processes, while those up-regulated in megagametophytes were enriched for accumulation of storage material and responses to stress. The highest number of DETs was detected during the initiation of the cleavage process. Transcripts related to embryogenic competence, cell wall modifications, cell division pattern, axis specification and response to hormones and stress were highly abundant and differentially expressed during early embryo development. The abundance of representative DETs was confirmed by qRT-PCR analyses.

CONCLUSION

Based on the processes identified in the GO enrichment analyses and the expression of the selected transcripts we suggest that (i) processes related to embryogenic competence and cell wall loosening are involved in activating the cleavage process; (ii) apical-basal polarization is strictly regulated in dominant embryos but not in the subordinate embryos; (iii) the transition from the morphogenic phase to the maturation phase is not completed in subordinate embryos. This is the first genome-wide transcript expression profiling of the earliest stages during embryo development in a Pinus species. Our results can serve as a framework for future studies to reveal the functions of identified genes.

Genome-wide analysis and expression characteristics of small auxin-up RNA (SAUR) genes in moso bamboo (Phyllostachys edulis)

Moso bamboo (Phyllostachys edulis) is well known for its rapid shoot growth. Auxin exerts pleiotropic effects on plant growth. The small auxin-up RNA (SAUR) genes are early auxin-responsive genes involved in plant growth. In total, 38 SAUR genes were identified in P. edulis (PheSAUR). A comprehensive overview of the PheSAUR gene family is presented, including the gene structures, phylogeny, and subcellular location predictions. A transcriptome analysis indicated that 37 (except PheSAUR18) of the PheSAUR genes were expressed during shoot growth process and that the PheSAUR genes were differentially expressed. Furthermore, quantitative real-time PCR analysis indicated that all of the PheSAUR genes could be induced in different tissues of seedlings and that 37 (except PheSAUR41) of the PheSAUR genes were up-regulated after indole-3-acetic acid (IAA) treatment. These results reveal a comprehensive overview of the PheSAUR gene family and may pave the way for deciphering their functions during bamboo development.

The effect of polar auxin transport on adventitious branches formation in Gracilaria lichenoides in vitro

Seaweed tissue culture (STC) is an important micropropagation tool that has been applied for strain improvement, micropropagation and genetic engineering. Because the mechanisms associated with STC are poorly understood, its application to these organisms lags far behind that of tissue culture propagation of higher plants. Auxin, calcium (Ca²⁺ ) and hydrogen peroxide (H₂ O₂ ) fluxes all play key roles during plant growth and development. In this study, we therefore measured indole-3-acetic acid, Ca²⁺ and H₂ O₂ fluxes of Gracilaria lichenoides explants during adventitious branches (ABs) formation for the first time using noninvasive micro-test technology. We confirmed that polar auxin transport (PAT) also occurs in the marine red alga G. lichenoides. We additionally found that N-1-naphthylphthalamic acid may suppress auxin efflux via ABCB1 transporters and then inhibit ABs formation from the apical region of G. lichenoides segments. The involvement of Ca²⁺ and H₂ O₂ fluxes in PAT-mediated AB formation in G. lichenoides was also investigated. We propose that complex feedback among Ca²⁺ , H₂ O₂ and auxin signaling and response systems may occur during ABs polar formation in G. lichenoides explants, similar to that in higher plants. Our results provide innovative insights that should aid future elucidation of mechanisms operative during STC.

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

Plant regeneration via somatic embryogenesis (SE) is the key step for genetic improvement of cotton (Gossypium hirsutum L.) through genetic engineering mediated by Agrobacteria, but the molecular mechanisms underlying SE in cotton is still unclear. Here, RNA-Sequencing was used to analyze the genes expressed during SE and their expression dynamics using RNAs isolated from non-embryogenic callus (NEC), embryogenic callus (EC) and somatic embryos (SEs). A total of 101, 670 unigenes were de novo assembled. The genes differentially expressed (DEGs) amongst NEC, EC and SEs were identified, annotated and classified. More DEGs were found between SEs and EC than between EC and NEC. A significant number of DEGs were related to hormone homeostasis, stress and ROS responses, and metabolism of polyamines. To confirm the expression dynamics of selected DEGs involved in various pathways, experiments were set up to investigate the effects of hormones (Indole-3-butytric acid, IBA; Kinetin, KT), polyamines, H2O2 and stresses on SE. Our results showed that exogenous application of IBA and KT positively regulated the development of EC and SEs, and that polyamines and H2O2 promoted the conversion of EC into SEs. Furthermore, we found that low and moderate stress is beneficial for proliferation of EC and SEs formation. Together, our global analysis of transcriptomic dynamics reveals that hormone homeostasis, polyamines, and stress response synergistically regulating SE in cotton.

miR393 contributes to the embryogenic transition induced in vitro in Arabidopsis via the modification of the tissue sensitivity to auxin treatment

miR393 was found to control embryogenic transition in somatic cells in Arabidopsis via control of the TIR1 and AFB2 auxin receptors genes of the F-box family. miR393 molecules are believed to regulate the expression of the auxin receptors of the TAAR clade. Considering the central role of auxin in the induction of somatic embryogenesis (SE) in plant explants cultured in vitro, the involvement of miR393 in the embryogenic transition of somatic cells has been hypothesised. To verify this assumption, the reporter, overexpressor and mutant lines in genes encoded MIR393 and TIR1/AFB proteins of the F-box family were analysed during SE in Arabidopsis. Expression profiling of MIR393a and MIR393b, mature miR393 and the target genes (TIR1, AFB1, AFB2, AFB3) were investigated in explants undergoing SE. In addition, the embryogenic potential of various genotypes with a modified activity of the MIR393 and TIR1/AFB targets was evaluated. The distinct increase in the accumulation of miR393 that was coupled with a notable down-regulation of TIR1 and AFB2 targets was observed at the early phase of SE induction. Relevant to this observation, the GUS/GFP monitored expression of MIR393, TIR1 and AFB2 transcripts was localised in explant tissue undergoing SE induction. The results suggest the miR393-mediated regulation of TIR1 and AFB2 during embryogenic transition induced in Arabidopsis and a modification of the explant sensitivity to auxin treatment is proposed as underlying this regulatory pathway.

ROS Homeostasis Regulates Somatic Embryogenesis via the Regulation of Auxin Signaling in Cotton

Somatic embryogenesis (S.E.) is a versatile model for understanding the mechanisms of plant embryogenesis and a useful tool for plant propagation. To decipher the intricate molecular program and potentially to control the parameters affecting the frequency of S.E., a proteomics approach based on two-dimensional gel electrophoresis (2-DE) combined with MALDI-TOF/TOF was used. A total of 149 unique differentially expressed proteins (DEPs) were identified at different stages of cotton S.E. compared with the initial control (0 h explants). The expression profile and functional annotation of these DEPs revealed that S.E. activated stress-related proteins, including several reactive oxygen species (ROS)-scavenging enzymes. Proteins implicated in metabolic, developmental, and reproductive processes were also identified. Further experiments were performed to confirm the role of ROS-scavenging enzymes, suggesting the involvement of ROS homeostasis during S.E. in cotton. Suppressing the expression of specifically identified GhAPX proteins resulted in the inhibition of dedifferentiation. Accelerated redifferentiation was observed in the suppression lines of GhAPXs or GhGSTL3 in parallel with the alteration of endogenous ascorbate metabolism and accumulation of endogenous H₂O₂ content. Moreover, disrupting endogenous redox homeostasis through the application of high concentrations of DPI, H₂O₂, BSO, or GSH inhibited the dedifferentiation of cotton explants. Mild oxidation induced through BSO treatment facilitated the transition from embryogenic calluses (ECs) to somatic embryos. Meanwhile, auxin homeostasis was altered through the perturbation of ROS homeostasis by chemical treatments or suppression of ROS-scavenging proteins, along with the activating/suppressing the transcription of genes related to auxin transportation and signaling. These results show that stress responses are activated during S.E. and may regulate the ROS homeostasis by interacting with auxin signaling.

Transcriptional identification and characterization of differentially expressed genes associated with embryogenesis in radish (Raphanus sativus L.)

Embryogenesis is an important component in the life cycle of most plant species. Due to the difficulty in embryo isolation, the global gene expression involved in plant embryogenesis, especially the early events following fertilization are largely unknown in radish. In this study, three cDNA libraries from ovules of radish before and after fertilization were sequenced using the Digital Gene Expression (DGE) tag profiling strategy. A total of 5,777 differentially expressed transcripts were detected based on pairwise comparison in the three libraries (0_DAP, 7_DAP and 15_DAP). Results from Gene Ontology (GO) and pathway enrichment analysis revealed that these differentially expressed genes (DEGs) were implicated in numerous life processes including embryo development and phytohormones biosynthesis. Notably, some genes encoding auxin response factor (ARF ), Leafy cotyledon1 (LEC1) and somatic embryogenesis receptor-like kinase (SERK ) known to be involved in radish embryogenesis were differentially expressed. The expression patterns of 30 genes including LEC1-2, AGL9, LRR, PKL and ARF8-1 were validated by qRT-PCR. Furthermore, the cooperation between miRNA and mRNA may play a pivotal role in the radish embryogenesis process. This is the first report on identification of DEGs profiles related to radish embryogenesis and seed development. These results could facilitate further dissection of the molecular mechanisms underlying embryogenesis and seed development in radish.

De novo comparative transcriptome analysis provides new insights into sucrose induced somatic embryogenesis in camphor tree (Cinnamomum camphora L.)

BACKGROUND

Somatic embryogenesis is a notable illustration of cell totipotency, by which somatic cells undergo dedifferentiation and then differentiate into somatic embryos. Our previous work demonstrated that pretreatment of immature zygotic embryos with 0.5 M sucrose solution for 72 h efficiently induced somatic embryo initiation in camphor tree. To better understand the molecular basis of somatic embryogenesis induced by osmotic stress, de novo transcriptome sequencing of three tissues of camphor tree (immature zygotic embryos, sucrose-pretreated immature zygotic embryos, and somatic embryos induced from sucrose-pretreated zygotic embryos) were conducted using Illumina Hiseq 2000 platform.

RESULTS

A total of 30.70 G high quality clean reads were obtained from cDNA libraries of the three samples. The overall de novo assembly of cDNA sequence data generated 205592 transcripts, with an average length of 998 bp. 114229 unigenes (55.56 % of all transcripts) with an average length of 680 bp were annotated with gene descriptions, gene ontology terms or metabolic pathways based on Blastx search against Nr, Nt, Swissprot, GO, COG/KOG, and KEGG databases. CEGMA software identified 237 out of 248 ultra-conserved core proteins as 'complete' in the transcriptome assembly, showing a completeness of 95.6 %. A total of 897 genes previously annotated to be potentially involved in somatic embryogenesis were identified. Comparative transcriptome analysis showed that a total of 3335 genes were differentially expressed in the three samples. The differentially expressed genes were divided into six groups based on K-means clustering. Expression level analysis of 52 somatic embryogenesis-related genes indicated a high correlation between RNA-seq and qRT-PCR data. Gene enrichment analysis showed significantly differential expression of genes responding to stress and stimulus.

CONCLUSIONS

The present work reported a de novo transcriptome assembly and global analysis focused on gene expression changes during initiation and formation of somatic embryos in camphor tree. Differential expression of somatic embryogenesis-related genes indicates that sucrose induced somatic embryogenesis may share or partly share the mechanisms of somatic embryogenesis induced by plant hormones. This study provides comprehensive transcript information and gene expression data for camphor tree. It could also serve as an important platform resource for further functional studies in plant embryogenesis.

Genome-Wide Analysis and Expression Profiling of the SUC and SWEET Gene Families of Sucrose Transporters in Oilseed Rape (Brassica napus L.)

Sucrose is the principal transported product of photosynthesis from source leaves to sink organs. SUTs/SUCs (sucrose transporters or sucrose carriers) and SWEETs (Sugars Will Eventually be Exported Transporters) play significant central roles in phloem loading and unloading. SUTs/SUCs and SWEETs are key players in sucrose translocation and are associated with crop yields. The SUT/SUC and SWEET genes have been characterized in several plant species, but a comprehensive analysis of these two gene families in oilseed rape has not yet been reported. In our study, 22 and 68 members of the SUT/SUCs and SWEET gene families, respectively, were identified in the oilseed rape (Brassica napus) genome through homology searches. An analysis of the chromosomal distribution, phylogenetic relationships, gene structures, motifs and the cis-acting regulatory elements in the promoters of BnSUC and BnSWEET genes were analyzed. Furthermore, we examined the expression of the 18 BnSUC and 16 BnSWEET genes in different tissues of "ZS11" and the expression of 9 BnSUC and 7 BnSWEET genes in "ZS11" under various conditions, including biotic stress (Sclerotinia sclerotiorum), abiotic stresses (drought, salt and heat), and hormone treatments (abscisic acid, auxin, cytokinin, brassinolide, gibberellin, and salicylic acid). In conclusion, our study provides the first comprehensive analysis of the oilseed rape SUC and SWEET gene families. Information regarding the phylogenetic relationships, gene structure and expression profiles of the SUC and SWEET genes in the different tissues of oilseed rape helps to identify candidates with potential roles in specific developmental processes. Our study advances our understanding of the important roles of sucrose transport in oilseed rape.

LEAFY COTYLEDON1-CASEIN KINASE I-TCP15-PHYTOCHROME INTERACTING FACTOR4 Network Regulates Somatic Embryogenesis by Regulating Auxin Homeostasis

Somatic embryogenesis (SE) is an efficient tool for the propagation of plant species and also, a useful model for studying the regulatory networks in embryo development. However, the regulatory networks underlying the transition from nonembryogenic callus to somatic embryos during SE remain poorly understood. Here, we describe an upland cotton (Gossypium hirsutum) CASEIN KINASE I gene, GhCKI, which is a unique key regulatory factor that strongly affects SE. Overexpressing GhCKI halted the formation of embryoids and plant regeneration because of a block in the transition from nonembryogenic callus to somatic embryos. In contrast, defective GhCKI in plants facilitated SE. To better understand the mechanism by which GhCKI regulates SE, the regulatory network was analyzed. A direct upstream negative regulator protein, cotton LEAFY COTYLEDON1, was identified to be targeted to a cis-element, CTTTTC, in the promoter of GhCKI. Moreover, GhCKI interacted with and phosphorylated cotton CINCINNATA-like TEOSINTE BRANCHED1-CYCLOIDEA-PCF transcription factor15 by coordinately regulating the expression of cotton PHYTOCHROME INTERACTING FACTOR4, finally disrupting auxin homeostasis, which led to increased cell proliferation and aborted somatic embryo formation in GhCKI-overexpressing somatic cells. Our results show a complex process of SE that is negatively regulated by GhCKI through a complex regulatory network.

SAUR Proteins as Effectors of Hormonal and Environmental Signals in Plant Growth

The plant hormone auxin regulates numerous aspects of plant growth and development. Early auxin response genes mediate its genomic effects on plant growth and development. Discovered in 1987, small auxin up RNAs (SAURs) are the largest family of early auxin response genes. SAUR functions have remained elusive, however, presumably due to extensive genetic redundancy. However, recent molecular, genetic, biochemical, and genomic studies have implicated SAURs in the regulation of a wide range of cellular, physiological, and developmental processes. Recently, crucial mechanistic insight into SAUR function was provided by the demonstration that SAURs inhibit PP2C.D phosphatases to activate plasma membrane (PM) H(+)-ATPases and promote cell expansion. In addition to auxin, several other hormones and environmental factors also regulate SAUR gene expression. We propose that SAURs are key effector outputs of hormonal and environmental signals that regulate plant growth and development.

The GRAS gene family in pine: transcript expression patterns associated with the maturation-related decline of competence to form adventitious roots

BACKGROUND

Adventitious rooting is an organogenic process by which roots are induced from differentiated cells other than those specified to develop roots. In forest tree species, age and maturation are barriers to adventitious root formation by stem cuttings. The mechanisms behind the respecification of fully differentiated progenitor cells, which underlies adventitious root formation, are unknown.

RESULTS

Here, the GRAS gene family in pine is characterized and the expression of a subset of these genes during adventitious rooting is reported. Comparative analyses of protein structures showed that pine GRAS members are conserved compared with their relatives in angiosperms. Relatively high GRAS mRNA levels were measured in non-differentiated proliferating embryogenic cultures and during embryo development. The mRNA levels of putative GRAS family transcription factors, including Pinus radiata's SCARECROW (SCR), PrSCR, and SCARECROW-LIKE (SCL) 6, PrSCL6, were significantly reduced or non-existent in adult tissues that no longer had the capacity to form adventitious roots, but were maintained or induced after the reprogramming of adult cells in rooting-competent tissues. A subset of genes, SHORT-ROOT (PrSHR), PrSCL1, PrSCL2, PrSCL10 and PrSCL12, was also expressed in an auxin-, age- or developmental-dependent manner during adventitious root formation.

CONCLUSIONS

The GRAS family of pine has been characterized by analyzing protein structures, phylogenetic relationships, conserved motifs and gene expression patterns. Individual genes within each group have acquired different and specialized functions, some of which could be related to the competence and reprogramming of adult cells to form adventitious roots.

Analysis of global gene expression profiles to identify differentially expressed genes critical for embryo development in Brassica rapa

Embryo development represents a crucial developmental period in the life cycle of flowering plants. To gain insights into the genetic programs that control embryo development in Brassica rapa L., RNA sequencing technology was used to perform transcriptome profiling analysis of B. rapa developing embryos. The results generated 42,906,229 sequence reads aligned with 32,941 genes. In total, 27,760, 28,871, 28,384, and 25,653 genes were identified from embryos at globular, heart, early cotyledon, and mature developmental stages, respectively, and analysis between stages revealed a subset of stage-specific genes. We next investigated 9,884 differentially expressed genes with more than fivefold changes in expression and false discovery rate ≤ 0.001 from three adjacent-stage comparisons; 1,514, 3,831, and 6,633 genes were detected between globular and heart stage embryo libraries, heart stage and early cotyledon stage, and early cotyledon and mature stage, respectively. Large numbers of genes related to cellular process, metabolism process, response to stimulus, and biological process were expressed during the early and middle stages of embryo development. Fatty acid biosynthesis, biosynthesis of secondary metabolites, and photosynthesis-related genes were expressed predominantly in embryos at the middle stage. Genes for lipid metabolism and storage proteins were highly expressed in the middle and late stages of embryo development. We also identified 911 transcription factor genes that show differential expression across embryo developmental stages. These results increase our understanding of the complex molecular and cellular events during embryo development in B. rapa and provide a foundation for future studies on other oilseed crops.

Microarray analysis of siberian ginseng cyclic somatic embryogenesis culture systems provides insight into molecular mechanisms of embryogenic cell cluster generation

Four systems of cyclic somatic embryogenesis of Siberian ginseng (Eleutherococcus senticosus Maxim) were used to study the mechanism of embryonic cell cluster generation. The first, direct somatic embryo induction (DSEI), generates secondary embryos directly from the primary somatic embryos; the second, direct embryogenic cell cluster induction (DEC)), induces embryogenic cell clusters directly from somatic embryos in agar medium. Subsequently, we found that when DEC-derived somatic embryos are transferred to suspension culture or a bioreactor culture, only somatic embryos are induced, and embryogenic cell clusters cannot form. Therefore, these new lines were named DEC cultured by liquid medium (ECS) and DEC cultured by bioreactor (ECB), respectively. Transmission electron microscopy showed that DEC epidermal cells contained a variety of inclusions, distinct from other lines. A cDNA library of DEC was constructed, and 1,948 gene clusters were obtained and used as probes. RNA was prepared from somatic embryos from each of the four lines and hybridized to a microarray. In DEC, 7 genes were specifically upregulated compared with the other three lines, and 4 genes were downregulated. EsXTH1 and EsPLT1, which were among the genes upregulated in DEC, were cloned using the rapid amplification of cDNA ends (RACE). Real-time quantitative PCR showed EsXTH1 was more highly expressed in DEC than in other lines throughout the culture cycle, and EsPLT1 expression in DEC increased as culture duration increased, but remained at a low expression level in other lines. These results suggest that EsXTH1 and EsPLT1 may be the essential genes that play important roles during the induction of embryogenic cell clusters.

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

As a product of asexual reproduction in plants, the somatic embryo (SE) differentiates into a new plantlet via a zygotic embryogenesis-like process. Here, we present the phenotypic and cellular differences between SEs and zygotic embryos (ZEs) revealed by histological section scanning using three parallel development stages of the two types of embryos of cotton (Gossypium hirsutum cv. YZ1), including globular, torpedo and cotyledonary-stages. To identify the molecular characteristics of SE development in cotton, the digital gene expression system was used to profile the genes active during SE and ZE development. A total of 4242 differentially expressed genes (DEGs) were identified in at least one developmental stage. Expression pattern and functional classification analysis based on these DEGs reveals that SE development exhibits a transcriptional activation of stress responses. RT-PCR analysis further confirmed enhanced expression levels of stress-related genes in SEs than in ZEs. Experimental stress treatment, induced by NaCl and ABA, accelerated SE development and increased the transcription of genes related to stress response, in parallel with decelerated proliferation of embryogenic calluses under stress treatment. Our data reveal that SE development involves the activation of stress responses, which we suggest may regulate the balance between cell proliferation and differentiation. These results provide new insight into the molecular mechanisms of SE development and suggest strategies that can be used for regulating the developmental processes of somatic embryogenesis.

AtWuschel promotes formation of the embryogenic callus in Gossypium hirsutum

Upland cotton (Gossypium hirsutum) is one of the most recalcitrant species for in vitro plant regeneration through somatic embryogenesis. Callus from only a few cultivars can produce embryogenic callus (EC), but the mechanism is not well elucidated. Here we screened a cultivar, CRI24, with high efficiency of EC produce. The expression of genes relevant to EC production was analyzed between the materials easy to or difficult to produce EC. Quantitative PCR showed that CRI24, which had a 100% EC differentiation rate, had the highest expression of the genes GhLEC1, GhLEC2, and GhFUS3. Three other cultivars, CRI12, CRI41, and Lu28 that formed few ECs expressed these genes only at low levels. Each of the genes involved in auxin transport (GhPIN7) and signaling (GhSHY2) was most highly expressed in CRI24, with low levels in the other three cultivars. WUSCHEL (WUS) is a homeodomain transcription factor that promotes the vegetative-to-embryogenic transition. We thus obtained the calli that ectopically expressed Arabidopsis thaliana Wus (AtWus) in G. hirsutum cultivar CRI12, with a consequent increase of 47.75% in EC differentiation rate compared with 0.61% for the control. Ectopic expression of AtWus in CRI12 resulted in upregulation of GhPIN7, GhSHY2, GhLEC1, GhLEC2, and GhFUS3. AtWus may therefore increase the differentiation potential of cotton callus by triggering the auxin transport and signaling pathways.

Direct reprogramming of adult somatic cells toward adventitious root formation in forest tree species: the effect of the juvenile-adult transition

Cellular plasticity refers, among others, to the capability of differentiated cells to switch the differentiation process and acquire new fates. One way by which plant cell plasticity is manifested is through de novo regeneration of organs from somatic differentiated cells in an ectopic location. However, switching the developmental program of adult cells prior to organ regeneration is difficult in many plant species, especially in forest tree species. In these species, a decline in the capacity to regenerate shoots, roots, or embryos from somatic differentiated cells is associated with tree age and maturation. The decline in the ability to form adventitious roots from stem cuttings is one of the most dramatic effects of maturation, and has been the subject of investigations on the basic nature of the process. Cell fate switches, both in plants and animals, are characterized by remarkable changes in the pattern of gene expression, as cells switch from the characteristic expression pattern of a somatic cell to a new one directing a new developmental pathway. Therefore, determining the way by which cells reset their gene expression pattern is crucial to understand cellular plasticity. The presence of specific cellular signaling pathways or tissue-specific factors underlying the establishment, maintenance, and redirection of gene expression patterns in the tissues involved in adventitious root formation could be crucial for cell fate switch and for the control of age-dependent cellular plasticity.

Combined small RNA and degradome sequencing reveals microRNA regulation during immature maize embryo dedifferentiation

Genetic transformation of maize is highly dependent on the development of embryonic calli from the dedifferentiated immature embryo. To better understand the regulatory mechanism of immature embryo dedifferentiation, we generated four small RNA and degradome libraries from samples representing the major stages of dedifferentiation. More than 186 million raw reads of small RNA and degradome sequence data were generated. We detected 102 known miRNAs belonging to 23 miRNA families. In total, we identified 51, 70 and 63 differentially expressed miRNAs (DEMs) in the stage I, II, III samples, respectively, compared to the control. However, only 6 miRNAs were continually up-regulated by more than fivefold throughout the process of dedifferentiation. A total of 87 genes were identified as the targets of 21 DEM families. This group of targets was enriched in members of four significant pathways including plant hormone signal transduction, antigen processing and presentation, ECM-receptor interaction, and alpha-linolenic acid metabolism. The hormone signal transduction pathway appeared to be particularly significant, involving 21 of the targets. While the targets of the most significant DEMs have been proved to play essential roles in cell dedifferentiation. Our results provide important information regarding the regulatory networks that control immature embryo dedifferentiation in maize.

Transcriptome profiling reveals auxin and cytokinin regulating somatic embryogenesis in different sister lines of cotton cultivar CCRI24

To get a broader view on the molecular mechanisms underlying somatic embryogenesis (SE) in cotton (Gossypium hirsutum L.), global analysis of cotton transcriptome dynamics during SE in different sister lines was performed using RNA-Seq. A total of 204 349 unigenes were detected by de novo assembly of the 214 977 462 Illumina reads. The quantitative reverse transcription-polymerase chain reaction (qRT-PCR) measurements were positively correlated with the RNA-Seq results for almost all the tested genes (R(2)  = 0.841, correlation was significant at the 0.01 level). Different phytohormone (auxin and cytokinin) concentration ratios in medium and the endogenous content changes of these two phytohormones at two stages in different sister lines suggested the roles of auxin and cytokinin during cotton SE. On the basis of global gene regulation of phytohormone-related genes, numerous genes from all the differentially expressed transcripts were involved in auxin and cytokinin biosynthesis and signal transduction pathways. Analyses of differentially expressed genes that were involved in these pathways revealed the substantial changes in gene type and abundance between two sister lines. Isolation, cloning and silencing/overexpressing the genes that revealed remarkable up- or down-expression during cotton SE were important. Furthermore, auxin and cytokinin play a primary role in SE, but potential cross-talk with each other or other factors remains unclear.

Molecular cloning and functional characterization of a novel cotton CBL-interacting protein kinase gene (GhCIPK6) reveals its involvement in multiple abiotic stress tolerance in transgenic plants

Plant CIPKs were specific Ser/Thr protein kinases, which were activated through interaction with calcineurin B-like protein (CBL) containing four EF hands for Ca(2+) binding. The CBL/CIPK complexes play an important role in signal transduction in biotic and abiotic stresses, as well as developmental processes. Here a Ser/Thr protein kinase gene (defined as GhCIPK6), which was isolated from RNA-Seq profile during cotton somatic embryogenesis in our previous research was characterized. The GhCIPK6 gene contains an ORF of 1296 bp that putatively encodes a polypeptide of 431 amino acids with a predicted molecular mass of 48.46 kDa and isoelectric point of 9.12. Sequence alignment analysis confirmed that GhCIPK6 has no intron, and it was homologous to AtCIPK6. Expression analysis of the GhCIPK6 suggested that they might function in diverse tissues, including styles and anthers but not fibers. In addition, expression of the GhCIPK6 gene was induced by salt, drought and ABA treatments. Overexpression of GhCIPK6 significantly enhances the tolerance to salt, drought and ABA stresses in transgenic Arabidopsis, indicating that GhCIPK6 acts as a positive regulator in response to salt and drought stress, and is supposed to be a potential candidate gene to improve stress tolerance by genetic manipulation in cotton and other crops.

Wuschel overexpression promotes somatic embryogenesis and induces organogenesis in cotton (Gossypium hirsutum L.) tissues cultured in vitro

This work shows that overexpression of the WUS gene from Arabidopsis enhanced the expression of embryogenic competence and triggered organogenesis from some cells of the regenerated embryo-like structures. Agrobacterium-mediated genetic transformation of cotton was described in the late 1980s, but is still time consuming and largely genotype dependant due to poor regeneration. To help solve this bottleneck, we over-expressed the WUSCHEL (WUS) gene, a homeobox transcription factor cloned in Arabidopsis thaliana, known to stimulate organogenesis and/or somatic embryogenesis in Arabidopsis tissues cultured in vitro. The AtWUS gene alone, and AtWUS gene fused to the GFP marker were compared to the GFP gene alone and to an empty construct used as a control. Somatic embryogenesis was improved in WUS expressed calli, as the percentage of explants giving rise to embryogenic tissues was significantly higher (×3) when WUS gene was over-expressed than in the control. An interesting result was that WUS embryogenic lines evolved in green embryo-like structures giving rise to ectopic organogenesis never observed in any of our previous transformation experiments. Using our standard in vitro culture protocol, the overexpression of AtWUS in tissues of a recalcitrant variety did not result in the production of regenerated plants. This achievement will still require the optimization of other non-genetic factors, such as the balance of exogenous phytohormones. However, our results suggest that targeted expression of the WUS gene is a promising strategy to improve gene transfer in recalcitrant cotton cultivars.

Small RNA and degradome sequencing reveal complex miRNA regulation during cotton somatic embryogenesis

MicroRNAs (miRNAs) are endogenous non-coding ~21 nucleotide RNAs that regulate gene expression at the transcriptional and post-transcriptional levels in plants and animals. They play an important role in development, abiotic stress, and pathogen responses. miRNAs with their targets have been widely studied in model plants, but limited knowledge is available on the small RNA population of cotton (Gossypium hirsutum)-an important economic crop, and global identification of related targets through degradome sequencing has not been developed previously. In this study, small RNAs and their targets were identified during cotton somatic embryogenesis (SE) through high-throughput small RNA and degradome sequencing, comparing seedling hypocotyl and embryogenic callus (EC) of G. hirsutum YZ1. A total of 36 known miRNA families were found to be differentially expressed, of which 19 miRNA families were represented by 29 precursors. Twenty-five novel miRNAs were identified. A total of 234 transcripts in EC and 322 transcripts in control (CK) were found to be the targets of 23 and 30 known miRNA families, respectively, and 16 transcripts were targeted by eight novel miRNAs. Interestingly, four trans-acting small interfering RNAs (tas3-siRNAs) were also found in degradome libraries, three of which perfectly matched their precursors. Several targets were further validated via RNA ligase-mediated rapid amplification of 5' cDNA ends (RLM 5'-RACE). The profiling of the miRNAs and their target genes provides new information on the miRNAs network during cotton SE.