Transcriptome analysis of callus from Picea balfouriana

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

BACKGROUND

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

RESULTS

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

CONCLUSION

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

Evaluating Tartary Buckwheat Genotypes with High Callus Induction Rates and the Transcriptomic Profiling during Callus Formation

Due to their complex genotypes, low in vitro regeneration rates, and difficulty in obtaining transgenic plants, studies concerning basic biological research and molecular breeding in Tartary buckwheat (TB) are greatly limited. In this study, the hypocotyls of 60 genotypes of TB (TBC1~60) were used as explants. Of these, TBC14 was selected due to a high callus induction rate of 97.78% under dark and a proliferation coefficient (PC) of 28.2 when cultured on MS medium supplemented with 2.0 mg/L of 2,4-D and 1.5 mg/L of 6-BA. Subsequently, the samples of the calli obtained from TBC14 were collected at 0, 10, 20, and 30 d, and their transcriptomes were sequenced where identified. GO enrichment led to the detection of the most significant active gene set, which was the DNA binding transcription factor activity. The DEGs related to the pathways concerning metabolism, the biosynthesis of secondary metabolites, and hormone signal transduction were the most enriched in the KEGG database. The sets of MYB, AP2/ERF, and bHLH TFs exhibited the highest number of DEGs. Using this enrichment analysis, 421 genes encoding TFs, 47 auxin- and cytokinin-related genes, and 6 signal transduction-associated genes were screened that may play significant roles in callus formation (CF) in TB. Furthermore, FtPinG0008123200.01 (bZIP), a key gene promoting CF, was screened in terms of the weighted gene co-expression network associated with the various stages of CF. Our study not only provides valuable information about the molecular mechanism of CF but also reveals new genes involved in this process.

Transcriptomic Profiling of Embryogenic and Non-Embryogenic Callus Provides New Insight into the Nature of Recalcitrance in Cannabis

Differential gene expression profiles of various cannabis calli including non-embryogenic and embryogenic (i.e., rooty and embryonic callus) were examined in this study to enhance our understanding of callus development in cannabis and facilitate the development of improved strategies for plant regeneration and biotechnological applications in this economically valuable crop. A total of 6118 genes displayed significant differential expression, with 1850 genes downregulated and 1873 genes upregulated in embryogenic callus compared to non-embryogenic callus. Notably, 196 phytohormone-related genes exhibited distinctly different expression patterns in the calli types, highlighting the crucial role of plant growth regulator (PGRs) signaling in callus development. Furthermore, 42 classes of transcription factors demonstrated differential expressions among the callus types, suggesting their involvement in the regulation of callus development. The evaluation of epigenetic-related genes revealed the differential expression of 247 genes in all callus types. Notably, histone deacetylases, chromatin remodeling factors, and EMBRYONIC FLOWER 2 emerged as key epigenetic-related genes, displaying upregulation in embryogenic calli compared to non-embryogenic calli. Their upregulation correlated with the repression of embryogenesis-related genes, including LEC2, AGL15, and BBM, presumably inhibiting the transition from embryogenic callus to somatic embryogenesis. These findings underscore the significance of epigenetic regulation in determining the developmental fate of cannabis callus. Generally, our results provide comprehensive insights into gene expression dynamics and molecular mechanisms underlying the development of diverse cannabis calli. The observed repression of auxin-dependent pathway-related genes may contribute to the recalcitrant nature of cannabis, shedding light on the challenges associated with efficient cannabis tissue culture and regeneration protocols.

Integrated transcriptomic and metabolomic analyses revealed the molecular mechanism of terpenoid formation for salicylic acid resistance in Pulsatilla chinensis callus

As a kind of traditional Chinese medicine, Pulsatilla chinensis (Bunge) Regel is well known for its anti-inflammation and anti-cancer activities, which are attributed to its active components including total saponins and monomers. To clarify the synthesis and metabolism mechanisms of class components in callus terpenes of P. chinensis, a certain concentration of salicylic acid (SA) hormone elicitor was added to the callus before being analysed by transcriptomic and metabolomic techniques. Results showed that the content of Pulsatilla saponin B4 in the callus suspension culture was significantly increased up to 1.99% with the addition of SA. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis showed that the differentially expressed genes were mainly enriched in 122 metabolic pathways, such as terpenoid metabolism-related pathways: terpenoid skeleton synthesis pathway, monoterpenoid biosynthesis pathways, diterpenoid biosynthesis pathways, and ubiquinone and other terpenoid-quinone biosynthesis pathways. A total of 31 differentially accumulated metabolites were obtained from four differential groups. Amongst 21 kinds of known chemical components in P. chinensis, deoxyloganic acid was the only monoterpenoid; the others are triterpenoids. In summary, this study found that SA elicitors can affect the metabolic changes of terpenoids in P. chinensis callus, which provided a basis for analysing the genetic regulation of terpenoid components of leucons.

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

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

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

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

Integrated transcriptome and endogenous hormone analysis provides new insights into callus proliferation in Osmanthus fragrans

Osmanthus fragrans is an important evergreen species with both medicinal and ornamental value in China. Given the low efficiency of callus proliferation and the difficulty of adventitious bud differentiation, tissue culture and regeneration systems have not been successfully established for this species. To understand the mechanism of callus proliferation, transcriptome sequencing and endogenous hormone content determination were performed from the initial growth stages to the early stages of senescence on O. fragrans calli. In total, 47,340 genes were identified by transcriptome sequencing, including 1798 previously unidentified genes specifically involved in callus development. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of differentially expressed genes (DEGs) was significantly enriched in plant hormone signal transduction pathways. Furthermore, our results from the orthogonal projections to latent structures discrimination analysis (OPLS-DA) of six typical hormones in five development stages of O. fragrans calli showed jasmonic acid (JA) could play important role in the initial stages of calli growth, whereas JA and auxin (IAA) were dominant in the early stages of calli senescence. Based on the weighted gene co-expression network analysis, OfSRC2, OfPP2CD5 and OfARR1, OfPYL3, OfEIL3b were selected as hub genes from the modules with the significant relevance to JA and IAA respectively. The gene regulation network and quantitative real-time PCR implied that during the initial stages of callus growth, the transcription factors (TFs) OfERF4 and OfMYC2a could down-regulate the expression of hub genes OfSRC2 and OfPP2CD5, resulting in decreased JA content and rapid callus growth; during the late stage of callus growth, the TFs OfERF4, OfMYC2a and OfTGA21c, OfHSFA1 could positively regulate the expression of hub genes OfSRC2, OfPP2CD5 and OfARR1, OfPYL3, OfEIL3b, respectively, leading to increased JA and IAA contents and inducing the senescence of O. fragrans calli. Hopefully, our results could provide new insights into the molecular mechanism of the proliferation of O. fragrans calli.

Protuberances are organized distinct regions of long-term callus: histological and transcriptomic analyses in kiwifruit

KEY MESSAGE

Macroscopic, ultrastructural, and molecular features-like a ball shape, the presence of starch granules, and the up-regulation of genes involved in carbohydrate metabolism and secondary metabolite biosynthesis-distinguish PT regions within a callus. The modification of the mass of pluripotent cells into de novo shoot bud regeneration is highly relevant to developmental biology and for agriculture and biotechnology. This study deals with protuberances (PT), structures that appear during the organogenic long-term culturing of callus (OC) in kiwifruit. These ball-shaped regions of callus might be considered the first morphological sign of the subsequent shoot bud development. Sections of PT show the regular arrangement of some cells, especially on the surface, in contrast to the regions of OC beyond the PT. The cells of OC possess chloroplasts; however, starch granules were observed only in PTs' plastids. Transcriptomic data revealed unique gene expression for each kind of sample: OC, PT, and PT with visible shoot buds (PT-SH). Higher expression of the gene involved in lipid (glycerol-3-phosphate acyltransferase 5 [GPAT5]), carbohydrate (granule-bound starch synthase 1 [GBSS1]), and secondary metabolite (beta-glucosidase 45 [BGL45]) pathways were detected in PT and could be proposed as the markers of these structures. The up-regulation of the regulatory associated protein of TOR (RAPTOR1) was found in PT-SH. The highest expression of the actinidain gene in leaves from two-year-old regenerated plants suggests that the synthesis of this protein takes place in fully developed organs. The findings indicate that PT and PT-SH are specific structures within OC but have more features in common with callus tissue than with organs.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Genotype and transcriptome effects on somatic embryogenesis in Cryptomeria japonica

Somatic embryogenesis (SE), which is in vitro regeneration of plant bodies from somatic cells, represents a useful means of clonal propagation and genetic engineering of forest trees. While protocols to obtain calluses and induce regeneration in somatic embryos have been reported for many tree species, the knowledge of molecular mechanisms of SE development is still insufficient to achieve an efficient supply of somatic embryos required for the industrial application. Cryptomeria japonica, a conifer species widely used for plantation forestry in Japan, is one of the tree species waiting for a secure SE protocol; the probability of normal embryo development appears to depend on genotype. To discriminate the embryogenic potential of embryonal masses (EMs) and efficiently obtain normal somatic embryos of C. japonica, we investigated the effects of genotype and transcriptome on the variation in embryogenic potential. Using an induction experiment with 12 EMs each from six genotypes, we showed that embryogenic potential differs between/within genotypes. Comparisons of gene expression profiles among EMs with different embryogenic potentials revealed that 742 differently expressed genes were mainly associated with pattern forming and metabolism. Thus, we suggest that not only genotype but also gene expression profiles can determine success in SE development. Consistent with previous findings for other conifer species, genes encoding leafy cotyledon, wuschel, germin-like proteins, and glutathione-S-transferases are likely to be involved in SE development in C. japonica and indeed highly expressed in EMs with high-embryogenic potential; therefore, these proteins represent candidate markers for distinguishing embryogenic potential.

Microspore embryogenesis: targeting the determinant factors of stress-induced cell reprogramming for crop improvement

Under stress, isolated microspores are reprogrammed in vitro towards embryogenesis, producing doubled haploid plants that are useful biotechnological tools in plant breeding as a source of new genetic variability, fixed in homozygous plants in only one generation. Stress-induced cell death and low rates of cell reprogramming are major factors that reduce yield. Knowledge gained in recent years has revealed that initiation and progression of microspore embryogenesis involve a complex network of factors, whose roles are not yet well understood. Here, I review recent findings on the determinant factors underlying stress-induced microspore embryogenesis, focusing on the role of autophagy, cell death, auxin, chromatin modifications, and the cell wall. Autophagy and cell death proteases are crucial players in the response to stress, while cell reprogramming and acquisition of totipotency are regulated by hormonal and epigenetic mechanisms. Auxin biosynthesis, transport, and action are required for microspore embryogenesis. Initial stages involve DNA hypomethylation, H3K9 demethylation, and H3/H4 acetylation. Cell wall remodelling, with pectin de-methylesterification and arabinogalactan protein expression, is necessary for embryo development. Recent reports show that treatments with small modulators of autophagy, proteases, and epigenetic marks reduce cell death and enhance embryogenesis initiation in several crops, opening up new possibilities for improving in vitro embryo production in breeding programmes.

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.

PICEAdatabase: a web database for Picea omics and phenotypic information

Picea belongs to the Pinaceae family and is a famous commercial tree species because of its straight trunk and excellent timber traits. Recently, omics have been widely used for fundamental and mechanism studies on Picea plants. To improve the accessibility to omics and phenotypic data and facilitate further studies, we compiled the sequences of 2 chloroplast genomes (Picea crassifolia and Picea asperata) and 32 complete omics data sets, including 20 transcriptomes, 4 proteomes, 2 degradomes and 6 microRNAs from P. crassifolia, P. asperata, Picea balfouriana and Picea abies tissues under different treatments, in PICEAdatabase. In addition, phenotypic data on plant growth and wood property traits were collected from two field trials of P. crassifolia. PICEAdatabase also includes useful analysis tools, such as BLAST, DESeq2 and JBrowse, to assist with analyses.

Comprehensive assembly and analysis of the transcriptome of maritime pine developing embryos

BACKGROUND

There are clear differences in embryo development between angiosperm and gymnosperm species. Most of the current knowledge on gene expression and regulation during plant embryo development has derived from studies on angiosperms species, in particular from the model plant Arabidopsis thaliana. The few published studies on transcript profiling of conifer embryogenesis show the existence of many putative embryo-specific transcripts without an assigned function. In order to extend the knowledge on the transcriptomic expression during conifer embryogenesis, we sequenced the transcriptome of zygotic embryos for several developmental stages that cover most of Pinus pinaster (maritime pine) embryogenesis.

RESULTS

Total RNA samples collected from five zygotic embryo developmental stages were sequenced with Illumina technology. A de novo transcriptome was assembled as no genome sequence is yet published for Pinus pinaster. The transcriptome of reference for the period of zygotic embryogenesis in maritime pine contains 67,429 transcripts, which likely encode 58,527 proteins. The annotation shows a significant percentage, 31%, of predicted proteins exclusively present in pine embryogenesis. Functional categories and enrichment analysis of the differentially expressed transcripts evidenced carbohydrate transport and metabolism over-representation in early embryo stages, as highlighted by the identification of many putative glycoside hydrolases, possibly associated with cell wall modification, and carbohydrate transport transcripts. Moreover, the predominance of chromatin remodelling events was detected in early to middle embryogenesis, associated with an active synthesis of histones and their post-translational modifiers related to increased transcription, as well as silencing of transposons.

CONCLUSIONS

Our results extend the understanding of gene expression and regulation during zygotic embryogenesis in conifers and are a valuable resource to support further improvements in somatic embryogenesis for vegetative propagation of conifer species. Specific transcripts associated with carbohydrate metabolism, monosaccharide transport and epigenetic regulation seem to play an important role in pine early embryogenesis and may be a source of reliable molecular markers for early embryogenesis.

Cell Wall Epitopes and Endoploidy as Reporters of Embryogenic Potential in Brachypodium Distachyon Callus Culture

Effective regeneration of callus tissue into embryos and then into whole plants is essential for plant biotechnology. The embryonic potential is often low and can further decrease with time in culture, which limits the utilisation of calli for transformation procedures and in vitro propagation. In this study, we show that the loss of embryogenic potential in callus cultures of Brachypodium distachyon is progressive over time. Flow cytometry analyses indicated endoploidy levels increased in 60- and 90-day-old calli with effective loss of the 2C DNA content peak in the latter. Analysis of indolic compounds content revealed a decrease in 60- and 90-day-old calli compared to either freshly isolated explants or 30-day-old calli. Immunohistochemical analysis revealed a decrease in arabinogalactan proteins (AGP) signal with the time of culture, but extensin (EXT) epitopes either increased (JIM12 epitopes) or decreased (JIM11 epitopes). The transcript accumulation levels of AGPs and EXTs confirmed these results, with most of AGP and EXT transcripts gradually decreasing. Some chimeric EXT transcripts significantly increased on the 30th day of culture, perhaps because of an increased embryogenic potential. Selected somatic embryogenesis-related genes and cyclins demonstrated a gradual decrease of transcript accumulation for YUCCA (YUC), AINTEGUMENTA-LIKE (AIL), BABY BOOM (BBM), and CLAVATA (CLV3) genes, as well as for most of the cyclins, starting from the 30th day of culture. Notably, WUSCHEL (WUS) transcript was detectable only on the 30th and 60th day and was not detectable in the zygotic embryos and in 90-day-old calli.

The role of miRNA in somatic embryogenesis

Somatic embryogenesis (SEG) is one of the best techniques for mass production of economically important plants. It is also used for the study of morphology, anatomy, physiology, genetics and molecular mechanism of embryo development. Somatic Embryos (SE) are bipolar structures that develop from a cell other than a gamete or zygote. SEG reflects the unique developmental potential of plant somatic cells, resulting in the transition of the differentiated somatic cells to embryogenic cells to follow the zygotic embryo stages. There are several biochemical and physiological processes that transformed a single somatic cell to a whole plant. SE studies provide insight into cell mechanisms governing the totipotency process in plants. Previously, in vitro studies have suggested the role of various regulatory genes in embryogenic transition that are triggered by plant hormones in response to stress. The omic studies identify the specific genes, transcripts, and proteins required for somatic embryogenesis development. MicroRNAs (miRNAs) are small, 19-24 nucleotides (nt), non-coding small RNA regulatory molecules controlling a large number of biological processes. In addition to their role in SEG, miRNAs play vital role in plant development, secondary metabolite synthesis and metabolism of macromolecules, hormone signal transduction, and tolerance of plants to biotic and abiotic stresses. During last decade several types of miRNAs involved in SEG have been reported. Among these miRNAs, miR156, miR162, miR166a, miR167, miR168, miR171a/b, miR171c, miR393, miR397 and miR398 played very active role during various stages of SEG. In this review, we highlighted the role of these as well as other miRNAs in some economically important plants.

Comprehensive transcriptome profiling of soybean leaves in response to simulated acid rain

As a source of edible oil and protein, soybean is a major globally important economic crop; Improving its production has been an important objective of soybean breeding. Acid rain has been shown to influence soybean growth and productivity, with consequent adverse impacts on its production for use by human populations. In this study, RNA sequencing technology was utilized to examine changes in gene expression when soybean was exposed to simulated acid rain (SAR). We sampled soybean leaves at five time intervals (0, 6, 30, 54, 78, and 102 h), and built the cDNA library. In total, 54,175 expression genes were found, including 2016 genes with differential expression. A total of 416 genes were considered, as they were closely related to the response to SAR. Genes related to the regulation of sulfur and nitrogen metabolism, carbohydrate metabolism, photosynthesis, and reactive oxygen species were among those differentially expressed in response to SAR. In this study, we examined the response mechanisms of soybean under SAR exposure. Our findings will improve our understanding of the molecular mechanisms employed by soybean in responding to abiotic stress, and therefore provides important information in developing soybean breeding to improve tolerance to these stresses.

Molecular Dissection of the Regenerative Capacity of Forest Tree Species: Special Focus on Conifers

Somatic embryogenesis (SE) and organogenesis have become leading biotechnologies for forest tree improvement and the implementation of multi-varietal forestry. Despite major advances in clonal propagation using these technologies, many forest tree species, such as conifers, show a low regeneration capacity. Developmental factors such as genotype, the type and age of the explant or tissue, and the age and maturity of the mother tree are limiting factors for the success of propagation programs. This review summarizes recent research on the molecular pathways involved in the regulation of key steps in SE and organogenesis of forest tree species, mainly conifers. The interaction between auxin and stress conditions, the induction of cell identity regulators and the role of cell wall remodeling are reviewed. This information is essential to develop tools and strategies to improve clonal propagation programs for forest tree species.

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.

Identification of novel miRNAs and miRNA expression profiling in embryogenic tissues of Picea balfouriana treated by 6-benzylaminopurine

Here, we compared miRNA expression profiles in embryonic cell cultures of the conifer Picea balfouriana following application of the synthetic cytokinin 6-benzylaminopurine (6-BAP). We used next-generation sequencing to analyze three libraries of small RNAs from the treated embryogenic cell cultures and generated 24,000,000 raw reads from each of the libraries. Over 70 differentially regulated micro RNA (miRNA) families (≥2 fold change in expression) were identified between pairs of treatments. A quantitative analysis showed that miR3633 and miR1026 were upregulated in tissues with the highest embryogenic ability. These two miRNAs were predicted to target genes encoding receptor-like protein kinase and GAMYB transcription factors, respectively. In one library, miR1160, miR5638, miR1315, and miR5225 were downregulated. These four miRNAs were predicted to target genes encoding APETALA2, calmodulin-binding protein, and calcium-dependent protein kinase transcription factors. The expression patterns of the miRNAs and their targets were negatively correlated. Approximately 181 potentially novel P. balfouriana miRNAs were predicted from the three libraries, and seven were validated during the quantitative analysis. This study is the first report of differential miRNA regulation in tissues treated with 6-BAP during somatic embryogenesis. The differentially expressed miRNAs will be of value for investigating the mechanisms of embryogenic processes that are responsive to 6-BAP in P. balfouriana.

Transcriptome Characterization of Gnetum parvifolium Reveals Candidate Genes Involved in Important Secondary Metabolic Pathways of Flavonoids and Stilbenoids

Gnetum is a small, unique group of Gnetophyta with a controversial phylogenetic position. Gnetum parvifolium is an important Chinese traditional medicinal plant, which is rich in bioactive compounds such as flavonoids and stilbenoids. These compounds provide significant medicinal effects, mostly as antioxidant, anticancer, and antibacterial agents. However, the mechanisms involved in the biosynthesis and regulation of these compounds in G. parvifolium are still unknown. In this study, we found that flavonoids and stilbene compounds accumulated at different levels in various tissues of G. parvifolium. We further obtained and analyzed massive sequence information from pooled samples of G. parvifolium by transcriptome sequencing, which generated 94,816 unigenes with an average length of 724 bp. Functional annotation of all these unigenes revealed that many of them were associated with several important secondary metabolism pathways including flavonoids and stilbenoids. In particular, several candidate unigenes (PAL-, C4H-, 4CL-, and STS-like genes) involved in stilbenoids biosynthesis were highly expressed in leaves and mature fruits. Furthermore, high temperature and UV-C strongly induced the expression of these genes and enhanced stilbene production (i.e., resveratrol and piceatannol) in leaves of young seedlings. Our present transcriptomic and biochemical data on secondary metabolites in G. parvifolium should encourage further investigation on evolution, ecology, functional genomics, and breeding of this plant with strong pharmaceutical potential.

Changes in the Metabolome of Picea balfouriana Embryogenic Tissues That Were Linked to Different Levels of 6-BAP by Gas Chromatography-Mass Spectrometry Approach

Embryogenic cultures of Picea balfouriana, which is an important commercial species for reforestation in Southern China, easily lose their embryogenic ability during long-term culture. Embryogenic tissue that proliferated at lower concentrations (3.6 μM and 2.5 μM) of 6-benzylaminopurine (6-BAP) were more productive, and generated 113 ± 6 and 89 ± 3 mature embryos per 100 mg embryogenic tissue, respectively. A metabolomic approach was used to study the changes in metabolites linked to embryogenic competence related to three different 6-BAP concentrations (2.5 μM, 3.6 μM, and 5 μM). A total of 309 compounds were obtained, among which 123 metabolites mapped to Kyoto Encyclopedia of Genes and genomes (KEGG) pathways. The levels of 35 metabolites were significantly differentially regulated among the three 6-BAP treatments, and 32 metabolites differed between the 2.5 μM and 5 μM treatments. A total of 17 metabolites appeared only once among the three comparisons. The combination of a score plot and a loading plot showed that in the samples with higher embryogenic ability (3.6 μM and 2.5 μM), up-regulated metabolites were mostly amino acids and down-regulated metabolites were mostly primary carbohydrates (especially sugars). These results suggested that 6-BAP may influence embryogenic competence by nitrogen metabolism, which could cause an increase in amino acid levels and higher amounts of aspartate, isoleucine, and leucine in tissues with higher embryogenic ability. Furthermore, we speculated that 6-BAP may affect the amount of tryptophan in tissues, which would change the indole-3-acetic acid levels and influence the embryogenic ability.

Global scale transcriptome analysis of Arabidopsis embryogenesis in vitro

Background

Somatic embryogenesis (SE) in plants is a process by which embryos are generated directly from somatic cells, rather than from the fused products of male and female gametes. Despite the detailed expression analysis of several somatic-to-embryonic marker genes, a comprehensive understanding of SE at a molecular level is still lacking. The present study was designed to generate high resolution transcriptome datasets for early SE providing the way for future research to understand the underlying molecular mechanisms that regulate this process. We sequenced Arabidopsis thaliana somatic embryos collected from three distinct developmental time-points (5, 10 and 15 d after in vitro culture) using the Illumina HiSeq 2000 platform.

Results

This study yielded a total of 426,001,826 sequence reads mapped to 26,520 genes in the A. thaliana reference genome. Analysis of embryonic cultures after 5 and 10 d showed differential expression of 1,195 genes; these included 778 genes that were more highly expressed after 5 d as compared to 10 d. Moreover, 1,718 genes were differentially expressed in embryonic cultures between 10 and 15 d. Our data also showed at least eight different expression patterns during early SE; the majority of genes are transcriptionally more active in embryos after 5 d. Comparison of transcriptomes derived from somatic embryos and leaf tissues revealed that at least 4,951 genes are transcriptionally more active in embryos than in the leaf; increased expression of genes involved in DNA cytosine methylation and histone deacetylation were noted in embryogenic tissues. In silico expression analysis based on microarray data found that approximately 5% of these genes are transcriptionally more active in somatic embryos than in actively dividing callus and non-dividing leaf tissues. Moreover, this identified 49 genes expressed at a higher level in somatic embryos than in other tissues. This included several genes with unknown function, as well as others related to oxidative and osmotic stress, and auxin signalling.

Conclusions

The transcriptome information provided here will form the foundation for future research on genetic and epigenetic control of plant embryogenesis at a molecular level. In follow-up studies, these data could be used to construct a regulatory network for SE; the genes more highly expressed in somatic embryos than in vegetative tissues can be considered as potential candidates to validate these networks.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1504-6) contains supplementary material, which is available to authorized users.

Transition from somatic embryo to friable embryogenic callus in cassava: dynamic changes in cellular structure, physiological status, and gene expression profiles

Friable embryogenic callus (FEC) is considered as the most suitable material for efficient genetic transformation of cassava. Heavy genotype dependence of FEC induction and amenability to somaclonal variation limits the production and maintenance of reliable FEC. Identifying key elements involved in biological processes from somatic embryos (SEs) to FEC at different stages provides critical insights for FEC improvement. Cytological observation showed a dramatic change of subcellular structures among SEs, fresh FEC (FFEC), and old FEC (OFEC). Decrease of sucrose and increase of fructose and glucose were detected in OFEC. A total of 6871 differentially expressed genes (DEGs) were identified from SEs, FFEC, and OFEC by RNA-seq. Analysis of the DEGs showed that FEC induction was accompanied by the process of dedifferentiation, whereas the epigenetics modification occurred during the continuous subculturing process. The cell structure was reconstructed, mainly including the GO terms of "cell periphery" and "external encapsulating structure"; in parallel, the internal mechanisms changed correspondingly, including the biological process of glycolysis and metabolisms of alanine, aspartate, and glutamate. The significant reduction of genomic DNA methylation in OFEC indicated altered gene expression via chromatin modification. These results indicate that the induction and long-term subculture of FEC is a complicated biological process involving changes of genome modification, gene expression, and subcellular reconstruction. The findings will be useful for improving FEC induction and maintenance from farmer-preferred cassava cultivars recalcitrant to genetic transformation, hence improving cassava through genetic engineering.