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

Transcriptome and Gene Expression Analysis Revealed CeNA1: A Potential New Marker for Somatic Embryogenesis in Common Centaury (Centaurium erythraea Rafn.)

Centaurium erythraea Rafn. is a medicinal plant used as a model for studying plant developmental processes due to its developmental plasticity and ease of manipulation in vitro. Identifying the genes involved in its organogenesis and somatic embryogenesis (SE) is the first step toward unraveling the molecular mechanisms underlying its morphogenic plasticity. Although SE is the most common method of centaury regeneration, the genes involved in this have not yet been identified. The aim of this study was to identify the differentially expressed genes (DEGs) during key stages of SE and organogenesis using transcriptome data, with a focus on novel SE-related genes. The transcriptomic analysis revealed a total of 4040 DEGs during SE and 12,708 during organogenesis. Gene Ontology (GO) annotation showed that the highest number of SE-related genes was involved in defense responses. The expression of fifteen selected SE-related candidate genes was assessed by RT-qPCR across nine centaury developmental stages, including embryogenic tissues. Notably, a newly reported transcript, named CeNA1, was specifically activated during embryogenic callus (ec) induction, making it a potential novel marker for early SE. These findings provide, for the first time, insight into SE-related transcriptional patterns, representing a step closer to uncovering the molecular basis of centaury's developmental plasticity.

The WOX Genes from the Intermediate Clade: Influence on the Somatic Embryogenesis in Medicago truncatula

Transcription factors from the WOX family are well-known regulators of cell proliferation and differentiation in plants. Herein, we focused on several WOX genes from the intermediate clade and checked their impact on somatic embryogenesis using the model legume object Medicago truncatula. As a result, we show that MtWOX9-1 overexpression not only stimulates somatic embryogenesis in the embryogenic M. truncatula line, as it was shown previously, but can also induce somatic embryogenesis in the non-embryogenic line. Other intermediate clade WOX, including the close paralog of MtWOX9-1, as well as WOX11 homologs, did not have any significant impact on somatic embryogenesis in our in vitro cultivation system. Together, our results give new information about the diversity of the WOX family proteins and their specific functions. These data can be used for the search of new regeneration stimulators.

Genome-Wide Identification and Expression Analysis of GATA Family Genes in Dimocarpus longan Lour

GATA transcription factors, which are DNA-binding proteins with type IV zinc finger binding domains, have a role in transcriptional regulation in biological organisms. They have an indispensable role in the growth and development of plants, as well as in improvements in their ability to face various environmental stresses. To date, GATAs have been identified in many gene families, but the GATA gene in longan (Dimocarpus longan Lour) has not been studied in previous explorations. Various aspects of genes in the longan GATA family, including their identification and classification, the distribution of their positions on chromosomes, their exon/intron structures, a synteny analysis, their expression at different temperatures, concentration of PEG, early developmental stages of somatic embryos and their expression levels in different tissues, and concentrations of exogenous hormones, were investigated in this study. This study showed that the 22 DlGATAs could be divided into four subfamilies. There were 10 pairs of homologous GATA genes in the synteny analysis of DlGATA and AtGATA. Four segmental replication motifs and one pair of tandem duplication events were present among the DlGATA family members. The cis-acting elements located in promoter regions were also found to be enriched with light-responsive elements, which contained related hormone-responsive elements. In somatic embryos, DlGATA4 is upregulated for expression at the globular embryo (GE) stage. We also found that DlGATA expression was strongly up-regulated in roots and stems. The study demonstrated the expression of DlGATA under hormone (ABA and IAA) treatments in embryogenic callus of longan. Under ABA treatment, DlGATA4 was up-regulated and the other DlGATA genes did not respond significantly. Moreover, as demonstrated with qRT-PCR, the expression of DlGATA genes showed strong up-regulated expression levels under 100 μmol·L-1 concentration IAA treatment. This experiment further studied these and simulated their possible connections with a drought response mechanism, while correlating them with their expression under PEG treatment. Overall, this experiment explored the GATA genes and dug into their evolution, structure, function, and expression profile, thus providing more information for a more in-depth study of the characteristics of the GATA family of genes.

Transcriptomic Analysis During the Induction of Somatic Embryogenesis in Coffea canephora

Omic tools have changed the way of doing research in experimental biology. The somatic embryogenesis (SE) study has not been immune to this benefit. The transcriptomic tools have been used to compare the genes expressed during the induction of SE with the genes expressed in zygotic embryogenesis or to compare the development of the different stages embryos go through. It has also been used to compare the expression of genes during the development of calli from which SE is induced, as well as many other applications. The protocol described here is employed in our laboratory to extract RNA and generate several transcriptomes for the study of SE on Coffea canephora.

An Introduction to Plant Cell, Tissue, and Organ Culture: Current Status and Perspectives

Plant cell, tissue, and organ cultures (PCTOC) have been used as experimental systems in basic research, allowing gene function demonstration through gene overexpression or repression and investigating the processes involved in embryogenesis and organogenesis or those related to the potential production of secondary metabolites, among others. On the other hand, PCTOC has also been applied at the commercial level for the vegetative multiplication (micropropagation) of diverse plant species, mainly ornamentals but also horticultural crops such as potato or fruit and tree species, and to produce high-quality disease-free plants. Moreover, PCTOC protocols are important auxiliary systems in crop breeding crops to generate pure lines (homozygous) to produce hybrids for the obtention of polyploid plants with higher yields or better performance. PCTOC has been utilized to preserve and conserve the germplasm of different crops or threatened species. Plant genetic improvement through genetic engineering and genome editing has been only possible thanks to the establishment of efficient in vitro plant regeneration protocols. Different companies currently focus on commercializing plant secondary metabolites with interesting biological activities using in vitro PCTOC. The impact of omics on PCTOC is discussed.

MtWOX2 and MtWOX9-1 Effects on the Embryogenic Callus Transcriptome in Medicago truncatula

WOX family transcription factors are well-known regulators of plant development, controlling cell proliferation and differentiation in diverse organs and tissues. Several WOX genes have been shown to participate in regeneration processes which take place in plant cell cultures in vitro, but the effects of most of them on tissue culture development have not been discovered yet. In this study, we evaluated the effects of MtWOX2 gene overexpression on the embryogenic callus development and transcriptomic state in Medicago truncatula. According to our results, overexpression of MtWOX2 leads to an increase in callus weight. Furthermore, transcriptomic changes in MtWOX2 overexpressing calli are, to a large extent, opposite to the changes caused by overexpression of MtWOX9-1, a somatic embryogenesis stimulator. These results add new information about the mechanisms of interaction between different WOX genes and can be useful for the search of new regeneration regulators.

Proteomic Approach during the Induction of Somatic Embryogenesis in Coffea canephora

Plant growth regulators (PGR) are essential for somatic embryogenesis (SE) in different species, and Coffea canephora is no exception. In our study model, previously, we have been able to elucidate the participation of various genes involved in SE by using different strategies; however, until now, we have not used a proteomic approach. This research seeks to contribute to understanding the primary cellular pathways involved in developing SE in C. canephora. The process of our model consists of two stages: (1) preconditioning in MS medium with auxin (NAA) and cytokinin (KIN), and (2) induction in Yasuda liquid medium added with cytokinin (BA). Therefore, in this study, we analyzed different days of the SE induction process using shotgun label-free proteomics. An amount of 1630 proteins was found among different sampling days of the process, of which the majority were accumulated during the induction stage. We found that some of the most enriched pathways during this process were the biosynthesis of amino acids and secondary metabolites. Eighteen proteins were found related to auxin homeostasis and two to cytokinin metabolism, such as ABC, BIG, ILR, LOG, and ARR. Ten proteins and transcription factors related to SE were also identified, like SERK1, SKP1, nuclear transcription factor Y, MADS-box, and calreticulin, and 19 related to other processes of plant development, among which the 14-3-3 and PP2A proteins stand out. This is the first report on the proteomic approach to elucidate the mechanisms that operate during the induction of SE in C. canephora. So, our findings provide the groundwork for future, more in-depth research. Data are available via ProteomeXchange Consortium via the PRIDE partner repository with the dataset identifier PXD047172.

Integrated Metabolomic and Transcriptomic Analyses Reveal the Potential Regulation of Flavonoids in the Production of Embryogenic Cultures during Early Somatic Embryogenesis of Longan (Dimocarpus longan Lour.)

Embryogenic cultures of longan (Dimocarpus longan Lour.) contain various metabolites with pharmacological properties that may function in the regulation of somatic embryogenesis (SE). In this study, based on widely targeted metabolomics, 501 metabolites were obtained from the embryogenic calli, incomplete compact proembryogenic cultures, and globular embryos during early SE of longan, among which 41 flavonoids were differentially accumulated during the SE. Using RNA sequencing, 36 flavonoid-biosynthesis-related genes and 43 MYB and 52 bHLH transcription factors were identified as differentially expressed genes. Furthermore, Kyoto Encyclopedia of Genes and Genomes enrichment analysis revealed that the flavonoid metabolism-related pathways were significantly enriched during the early SE. These results suggested that the changes in flavonoid levels in the embryogenic cultures of longan were mediated by MYBs and bHLHs via regulating flavonoid-biosynthesis-related genes, thus potentially regulating early SE. The identified metabolites in the embryogenic cultures of longan can be used to develop pharmaceutical ingredients.

Single-cell RNA sequencing analysis of the embryogenic callus clarifies the spatiotemporal developmental trajectories of the early somatic embryo in Dimocarpus longan

Plant embryogenic calli (ECs) can undergo somatic embryogenesis to regenerate plants. This process is mediated by regulatory factors, such as transcription factors and specifically expressed genes, but the precise molecular mechanisms underlying somatic embryogenesis at the single-cell level remain unclear. In this study, we performed high-resolution single-cell RNA sequencing analysis to determine the cellular changes in the EC of the woody plant species Dimocarpus longan (longan) and clarify the continuous cell differentiation trajectories at the transcriptome level. The highly heterogeneous cells in the EC were divided into 12 putative clusters (e.g., proliferating, meristematic, vascular, and epidermal cell clusters). We determined cluster-enriched expression marker genes and found that overexpression of the epidermal cell marker gene GDSL ESTERASE/LIPASE-1 inhibited the hydrolysis of triacylglycerol. In addition, the stability of autophagy was critical for the somatic embryogenesis of longan. The pseudo-timeline analysis elucidated the continuous cell differentiation trajectories from early embryonic cell division to vascular and epidermal cell differentiation during the somatic embryogenesis of longan. Moreover, key transcriptional regulators associated with cell fates were revealed. We found that ETHYLENE RESPONSIVE FACTOR 6 was characterized as a heat-sensitive factor that negatively regulates longan somatic embryogenesis under high-temperature stress conditions. The results of this study provide new spatiotemporal insights into cell division and differentiation during longan somatic embryogenesis at single-cell resolution.

Genome-wide high-throughput chromosome conformation capture analysis reveals hierarchical chromatin interactions during early somatic embryogenesis

Somatic embryogenesis (SE), like zygotic embryo development, is a progressive process. Early SE is the beginning of a switch from a somatic to an embryogenic state and is an important stage for initiating chromatin reprogramming of SE. Previous studies suggest that changes in chromatin accessibility occur during early SE, although information on the 3D structure of chromatin is not yet available. Here, we present a chromosome-level genome assembly of longan (Dimocarpus longan) using PacBio combined with high-through chromosome conformation capture scaffolding, which resulted in a 446 Mb genome assembly anchored onto 15 scaffolds. During early SE, chromatin was concentrated and then decondensed, and a large number of long terminal repeat retrotransposons (LTR-RTs) were enriched in the local chromatin interaction region, suggesting LTR-RTs were involved in chromatin reorganization. Early SE was accompanied by the transformation from A to B compartments, and the interactions between B compartments were enhanced. Results from chromatin accessibility, monomethylation of histone H3 at lysine 4 (H3K4me1) modification, and transcription analyses further revealed a gene regulatory network for cell wall thickening during SE. Particularly, we found that the H3K4me1 differential peak binding motif showed abnormal activation of ethylene response factor transcription factors and participation in SE. The chromosome-level genomic and multiomics analyses revealed the 3D conformation of chromatin during early SE, providing insight into the molecular mechanisms underlying cell wall thickening and the potential regulatory networks of TFs during early SE in D. longan. These results provide additional clues for revealing the molecular mechanisms of plant SE.

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

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

Apoplastic and Symplasmic Markers of Somatic Embryogenesis

Somatic embryogenesis (SE) is a process that scientists have been trying to understand for many years because, on the one hand, it is a manifestation of the totipotency of plant cells, so it enables the study of the mechanisms regulating this process, and, on the other hand, it is an important method of plant propagation. Using SE in basic research and in practice is invaluable. This article describes the latest, but also historical, information on changes in the chemical composition of the cell wall during the transition of cells from the somatic to embryogenic state, and the importance of symplasmic communication during SE. Among wall chemical components, different pectic, AGP, extensin epitopes, and lipid transfer proteins have been discussed as potential apoplastic markers of explant cells during the acquisition of embryogenic competence. The role of symplasmic communication/isolation during SE has also been discussed, paying particular attention to the formation of symplasmic domains within and between cells that carry out different developmental processes. Information about the number and functionality of plasmodesmata (PD) and callose deposition as the main player in symplasmic isolation has also been presented.

Identification of WUSCHEL-related homeobox (WOX) gene family members and determination of their expression profiles during somatic embryogenesis in Phoebe bournei

WUSCHEL-related homeobox (WOX) transcription factor (TF)-encoding genes play crucial roles during embryo development. The function of WOX genes in embryonic development has been thoroughly studied in Arabidopsis thaliana, but little is known about their function in woody species, especially Phoebe bournei, an endemic and endangered species in China. In the present study, a total of 15 WOX genes were identified in P. bournei, and phylogenetic analysis resulted in their assignment to three typical clades: an ancient clade, an intermediate clade, and a modern/WUS clade. The gene structure and sequence characteristics and the physicochemical properties of WOX proteins were also analyzed. Promoter prediction indicated that WOX genes are likely involved in plant growth and development and hormone responses. Subsequently, we evaluated the expression patterns of WOX genes in response to auxin (IAA), abscisic acid (ABA), and methyl jasmonate (MeJA) treatments. According to tissue-specific expression patterns, we screened nine WOX genes that were present in embryonic calli and that might participate in the somatic embryogenesis (SE) of P. bournei. Furthermore, the expression profiles of these nine WOX genes during three phases of embryogenic calli development and three phases of somatic embryo development, namely, spheroid embryogenesis, immature cotyledon-producing embryogenesis and mature cotyledon-producing embryogenesis, were monitored. Overall, we systematically analyzed the expression patterns of WOX genes in P. bournei during SE, the information of which provides a basis for further elucidating the molecular mechanism through which WOX TFs function in P. bournei embryo development.

Global transcriptome profiling reveals differential regulatory, metabolic and hormonal networks during somatic embryogenesis in Coffea arabica

BACKGROUND

Somatic embryogenesis (SE) is one of the most promising processes for large-scale dissemination of elite varieties. However, for many plant species, optimizing SE protocols still relies on a trial and error approach. We report the first global scale transcriptome profiling performed at all developmental stages of SE in coffee to unravel the mechanisms that regulate cell fate and totipotency.

RESULTS

RNA-seq of 48 samples (12 developmental stages × 4 biological replicates) generated 90 million high quality reads per sample, approximately 74% of which were uniquely mapped to the Arabica genome. First, the statistical analysis of transcript data clearly grouped SE developmental stages into seven important phases (Leaf, Dedifferentiation, Primary callus, Embryogenic callus, Embryogenic cell clusters, Redifferentiation and Embryo) enabling the identification of six key developmental phase switches, which are strategic for the overall biological efficiency of embryo regeneration. Differential gene expression and functional analysis showed that genes encoding transcription factors, stress-related genes, metabolism-related genes and hormone signaling-related genes were significantly enriched. Second, the standard environmental drivers used to control SE, i.e. light, growth regulators and cell density, were clearly perceived at the molecular level at different developmental stages. Third, expression profiles of auxin-related genes, transcription factor-related genes and secondary metabolism-related genes were analyzed during SE. Gene co-expression networks were also inferred. Auxin-related genes were upregulated during dedifferentiation and redifferentiation while transcription factor-related genes were switched on from the embryogenic callus and onward. Secondary metabolism-related genes were switched off during dedifferentiation and switched back on at the onset of redifferentiation. Secondary metabolites and endogenous IAA content were tightly linked with their respective gene expression. Lastly, comparing Arabica embryogenic and non-embryogenic cell transcriptomes enabled the identification of biological processes involved in the acquisition of embryogenic capacity.

CONCLUSIONS

The present analysis showed that transcript fingerprints are discriminating signatures of cell fate and are under the direct influence of environmental drivers. A total of 23 molecular candidates were successfully identified overall the 12 developmental stages and can be tested in many plant species to optimize SE protocols in a rational way.

CRISPR in Modulating Antibiotic Resistance of ESKAPE Pathogens

The ESKAPE (Enterococcus spp., Staphylococcus aureus, Klebsiella pneumoniae, Acinetobacter baumannii, Pseudomonas aeruginosa, and Enterobacter spp.) isolates both from the clinical settings and food products are demonstrated to gain resistance to multiple antimicrobials. Therefore, the ESKAPE pathogens pose a serious threat to public health, which warrants specific attention to developing alternative novel therapeutics. The clustered regularly interspaced short palindromic repeats associated (CRISPR-Cas) system is one of the novel methods for managing antibiotic-resistant strains. Specific Cas nucleases can be programmed against bacterial genomic sequences to decrease bacterial resistance to antibiotics. Moreover, a few CRISPR-Cas nucleases have the ability to the sequence-specific killing of bacterial strains. However, some pathogens acquire antibiotic resistance due to the presence of the CRISPR-Cas system. In brief, there is a wide range of functional diversity of CRISPR-Cas systems in bacterial pathogens. Hence, to be an effective and safe infection treatment strategy, a comprehensive understanding of the role of CRISPR-Cas systems in modulating antibiotic resistance in ESKAPE pathogens is essential. The present review summarizes all the mechanisms by which CRISPR confers and prevents antibiotic resistance in ESKAPE. The review also emphasizes the relationship between CRISPR-Cas systems, biofilm formation, and antibiotic resistance in ESKAPE.

Changes and transcriptome regulation of endogenous hormones during somatic embryogenesis in Ormosia henryi Prain

Introduction

Ormosia henryi is a rare and endangered plant growing in southern China. Somatic embryo culture is an effective measure for the rapid propagation of O. henryi. It has not been reported how regulatory genes induce somatic embryogenesis by regulating endogenous hormone changes during the process of somatic embryogenesis in O. henryi.

Methods

In this study, we analysed the endogenous hormone levels and transcriptome data of nonembryogenic callus (NEC), embryogenic callus (EC), globular embryo (GE) and cotyledon embryo (CE) in O. henryi.

Results

The results showed that the indole-3-acetic acid (IAA) content was higher and the cytokinins (CKs) content was lower in EC than in NEC, and the gibberellins (GAs) and abscisic acid (ABA) contents were significantly higher in NEC than in EC. The contents of IAA, CKs, GAs and ABA increased significantly with EC development. The expression patterns of differentially expressed genes (DEGs) involved in the biosynthesis and signal transduction of auxin (AUX) (YUCCA and SAUR), CKs (B-ARR), GAs (GA3ox, GA20ox, GID1 and DELLA) and ABA (ZEP, ABA2, AAO3, CYP97A3, PYL and ABF) were consistent with the levels of endogenous hormones during somatic embryogenesis (SE). In this study, 316 different transcription factors (TFs) regulating phytohormones were detected during SE. AUX/IAA were downregulated in the process of EC formation and GE differentiation into CE, but other TFs were upregulated and downregulated.

Conclusion

Therefore, we believe that relatively high IAA content and low CKs, GAs and ABA contents contribute to EC formation. The differential expression of AUX, CKs, GAs and ABA biosynthesis and signal transduction genes affected the endogenous hormone levels at different stages of SE in O. henryi. The downregulated expression of AUX/IAA inhibited NEC induction, promoted EC formation and GE differentiation into CE.

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.

miR156 regulates somatic embryogenesis by modulating starch accumulation in citrus

Somatic embryogenesis (SE) is a major regeneration approach for in vitro cultured tissues of plants, including citrus. However, SE capability is difficult to maintain, and recalcitrance to SE has become a major obstacle to plant biotechnology. We previously reported that miR156-SPL modules regulate SE in citrus callus. However, the downstream regulatory pathway of the miR156-SPL module in SE remains unclear. In this study, we found that transcription factors CsAGL15 and CsFUS3 bind to the CsMIR156A promoter and activate its expression. Suppression of csi-miR156a function leads to up-regulation of four target genes, SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (CsSPL) genes, and reduction of SE efficiency. In the short tandem target mimic (STTM)-miR156a overexpression callus (MIM156), the number of amyloplasts and starch content were significantly reduced, and genes involved in starch synthesis and transport were down-regulated. csi-miR172d was down-regulated, whereas the target genes, CsTOE1.1 and CsTOE1.2, which inhibit the expression of starch biosynthesis genes, were up-regulated. In our working model, CsAGL15 and CsFUS3 activate csi-miR156a, which represses CsSPLs and further regulates csi-miR172d and CsTOEs, thus altering starch accumulation in callus cells and regulating SE in citrus. This study elucidates the pathway of miR156-SPLs and miR172-TOEs-mediated regulation of SE, and provides new insights into enhancing SE capability in citrus.

A Comparative Transcriptome Analysis Reveals the Molecular Mechanisms That Underlie Somatic Embryogenesis in Peaonia ostii ‘Fengdan’

Low propagation rate is the primary problem that limits industry development of tree peony. In this study, a highly efficient regeneration system for tree peony using somatic embryogenesis (SE) was established. The transcriptomes of zygotic embryo explants (S0), non-embryonic callus (S1), embryonic callus (S2), somatic embryos (S3), and regenerated shoots (S4) were analyzed to determine the regulatory mechanisms that underlie SE in tree peony. The differentially expressed genes (DEGs) were identified in the pairwise comparisons of S1-vs-S2 and S1-vs-S3, respectively. The enriched DEGs were primarily involved in hormone signal transduction, stress response and the nucleus (epigenetic modifications). The results indicated that cell division, particularly asymmetric cell division, was enhanced in S3. Moreover, the genes implicated in cell fate determination played central roles in S3. Hormone signal pathways work in concert with epigenetic modifications and stress responses to regulate SE. SERK, WOX9, BBM, FUS3, CUC, and WUS were characterized as the molecular markers for tree peony SE. To our knowledge, this is the first study of the SE of tree peony using transcriptome sequencing. These results will improve our understanding of the molecular mechanisms that underly SE in tree peony and will benefit the propagation and genetic engineering of this plant.

Genome-wide circular RNA profiling and competing endogenous RNA regulatory network analysis provide new insights into the molecular mechanisms underlying early somatic embryogenesis in Dimocarpus longan Lour

Circular RNAs (circRNAs) are widely involved in plant growth and development. However, the function of circRNAs in plant somatic embryogenesis (SE) remains elusive. Here, by using high-throughput sequencing, a total of 5029 circRNAs were identified in the three stages of longan (Dimocarpus longan Lour.) early SE. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed that differentially expressed (DE) circRNA host genes were enriched in the 'non-homologous end-joining' (NHEJ) and 'butanoate metabolism' pathways. In addition, the reactive oxygen species (ROS) content during longan early SE was determined. The results indicated that ROS-induced DNA double-strand breaks may not depend on the NHEJ repair pathway. Correlation analyses of the levels of related metabolites (glutamate, γ-aminobutyrate and pyruvate) and the expression levels of circRNAs and their host genes involved in butanoate metabolism were performed. The results suggested that circRNAs may act as regulators of the expression of cognate mRNAs, thereby affecting the accumulation of related compounds. A competing endogenous RNA (ceRNA) network of DE circRNAs, DE mRNAs, DE long noncoding RNAs (lncRNAs) and DE microRNAs (miRNAs) was constructed. The results showed that the putative targets of the noncoding RNA (ncRNAs) were significantly enriched in the KEGG pathways 'mitogen-activated protein kinase signaling' and 'nitrogen metabolism'. Furthermore, the expression patterns of the candidate circRNAs, lncRNAs, miRNAs and mRNAs confirmed the negative correlation between miRNAs and ceRNAs. In addition, two circRNA overexpression vectors were constructed to further verify the ceRNA network correlations in longan early SE. Our study revealed the potential role of circRNAs in longan early SE, providing new insights into the intricate regulatory mechanism underlying plant SE.

Genome-wide identification, evolution analysis of cytochrome P450 monooxygenase multigene family and their expression patterns during the early somatic embryogenesis in Dimocarpus longan Lour

Cytochrome P450 (CYP), a multi-gene superfamily, is involved in a broad range of physiological processes, including hormone responses and secondary metabolism throughout the plant life cycle. Longan (Dimocarpus longan), a subtropical and tropical evergreen fruit tree, its embryonic development is closely related to the yield and quality of fruits. And a large number of secondary metabolites, such as flavonoids and carotenoids, are also produced during the longan somatic embryogenesis (SE). It is important, therefore, to study potential functions of CYPs in longan. However, the knowledge of longan CYPs is still very limited. Here, a total of 327 DlCYPs were identified using the genome-search method, which could be classified into nine clans. The expansion of the DlCYP family was mainly caused by tandem duplication (TD) events. Promoter cis-acting elements analysis elucidated that DlCYPs played important roles in hormonal responses. A total of 246 DlCYPs exhibited six different expression patterns during the early SE based on longan transcriptomic data. Eight DlCYPs underwent alternative splicing (AS) events, and they might produce one to six isoforms. And the AS transcript of DlCYP97C1 might act as an alternative to the full-length transcript in ICpEC and GE stages. Finally, protein-protein interaction (PPI) networks and miRNA target prediction elucidated that DlCYPs might be involved in the phenylpropanoid metabolic pathway and primarily regulated and targeted by miR413. In summary, our results provided valuable inventory for understanding the classification and biological functions of DlCYPs and provided insight into further functional verification of DlCYPs during the longan early SE.

miR171 modulates induction of somatic embryogenesis in citrus callus

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

CRISPR/Cas9-Mediated Targeted Mutagenesis of GmUGT Enhanced Soybean Resistance Against Leaf-Chewing Insects Through Flavonoids Biosynthesis

Leaf-chewing insects are important pests that cause yield loss and reduce seed quality in soybeans (Glycine max). Breeding soybean varieties that are resistant to leaf-chewing insects can minimize the need for insecticide use and reduce yield loss. The marker gene for QTL-M, Glyma.07g110300 (LOC100775351) that encodes a UDP-glycosyltransferase (UGT) is the major determinant of resistance against leaf-chewing insects in soybean; it exhibits a loss of function in insect-resistant soybean germplasms. In this study, Agrobacterium-mediated transformation introduced the CRISPR/Cas9 expression vector into the soybean cultivar Tianlong No. 1 to generate Glyma.07g110300-gene mutants. We obtained two novel types of mutations, a 33-bp deletion and a single-bp insertion in the GmUGT coding region, which resulted in an enhanced resistance to Helicoverpa armigera and Spodoptera litura. Additionally, overexpressing GmUGT produced soybean varieties that were more sensitive to H. armigera and S. litura. Both mutant and overexpressing lines exhibited no obvious phenotypic changes. The difference in metabolites and gene expression suggested that GmUGT is involved in imparting resistance to leaf-chewing insects by altering the flavonoid content and expression patterns of genes related to flavonoid biosynthesis and defense. Furthermore, ectopic expression of the GmUGT gene in the ugt72b1 mutant of Arabidopsis substantially rescued the phenotype of H. armigera resistance in the atugt72b1 mutant. Our study presents a strategy for increasing resistance against leaf-chewing insects in soybean through CRISPR/Cas9-mediated targeted mutagenesis of the UGT genes.

Global Transcriptome and Coexpression Network Analyses Reveal New Insights Into Somatic Embryogenesis in Hybrid Sweetgum (Liquidambar styraciflua × Liquidambar formosana)

Somatic embryogenesis (SE) is a process of somatic cells that dedifferentiate to totipotent embryonic stem cells and generate embryos in vitro. Despite recent scientific headway in deciphering the difficulties of somatic embryogenesis, the overall picture of key genes, pathways, and co-expression networks regulating SE is still fragmented. Therefore, deciphering the molecular basis of somatic embryogenesis of hybrid sweetgum remains pertinent. In the present study, we analyzed the transcriptome profiles and gene expression regulation changes via RNA sequencing from three distinct developmental stages of hybrid sweetgum: non-embryogenic callus (NEC), embryogenic callus (EC), and redifferentiation. Comparative transcriptome analysis showed that 19,957 genes were differentially expressed in ten pairwise comparisons of SE. Among these, plant hormone signaling-related genes, especially the auxin and cytokinin signaling components, were significantly enriched in NEC and EC early. The K-means method was used to identify multiple transcription factors, including HB-WOX, B3-ARF, AP2/ERF, and GRFs (growth regulating factors). These transcription factors showed distinct stage- or tissue-specific expression patterns mirroring each of the 12 superclusters to which they belonged. For example, the WOX transcription factor family was expressed only at NEC and EC stages, ARF transcription factor was expressed in EC early, and GRFs was expressed in late SE. It was noteworthy that the AP2/ERF transcription factor family was expressed during the whole SE process, but almost not in roots, stems and leaves. A weighted gene co-expression network analysis (WGCNA) was used in conjunction with the gene expression profiles to recognize the genes and modules that may associate with specific tissues and stages. We constructed co-expression networks and revealed 22 gene modules. Four of these modules with properties relating to embryonic potential, early somatic embryogenesis, and somatic embryo development, as well as some hub genes, were identified for further functional studied. Through a combination analysis of WGCNA and K-means, SE-related genes including AUX22, ABI3, ARF3, ARF5, AIL1, AIL5, AGL15, WOX11, WOX9, IAA29, BBM1, MYB36, LEA6, SMR4 and others were obtained, indicating that these genes play an important role in the processes underlying the progression from EC to somatic embryos (SEs) morphogenesis. The transcriptome information provided here will form the foundation for future research on genetic transformation 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; Key genes obtained from coexpression network analysis at each critical stage of somatic embryo can be considered as potential candidate genes to verify these networks.

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.

Early nodulin 93 protein gene: essential for induction of somatic embryogenesis in oil palm

Transcript profiling during the early induction phase of oil palm tissue culture and RNAi studies in a model somatic embryogenesis system showed that EgENOD93 expression is essential for somatic embryogenesis. Micropropagation of oil palm through tissue culture is vital for the generation of superior and uniform elite planting materials. Studies were carried out to identify genes to distinguish between leaf explants with the potential to develop into embryogenic or non-embryogenic callus. Oil palm cDNA microarrays were co-hybridized with cDNA probes of reference tissue, separately with embryo forming (media T527) and non-embryo (media T694) forming leaf explants sampled at Day 7, Day 14 and Day 21. Analysis of the normalized datasets has identified 77, 115 and 127 significantly differentially expressed genes at Day 7, Day 14, and Day 21, respectively. An early nodulin 93 protein gene (ENOD93), was highly expressed at Day 7, Day 14, and Day 21 and in callus (media T527), as assessed by RT-qPCR. Validation of EgENOD93 across tissue culture lines of different genetic background and media composition showed the potential of this gene as an embryogenic marker. In situ RNA hybridization and functional characterization in Medicago truncatula provided additional evidence that ENOD93 is essential for somatic embryogenesis. This study supports the suitability of EgENOD93 as a marker to predict the potential of leaf explants to produce embryogenic callus. Crosstalk among stresses, auxin, and Nod-factor like signalling molecules likely induces the expression of EgENOD93 for embryogenic callus formation.

Genome-wide identification and expression analysis of the SR gene family in longan (Dimocarpus longan Lour.)

Longan (Dimocarpus longan Lour.) is an important commercial fruit tree in southern China. The embryogenesis of longan affects the quality and yield of fruit. A large number of alternative splicing events occurs during somatic embryogenesis (SE), which is regulated by serine/arginine-rich (SR) proteins. However, the functions of SR proteins in longan are poorly understood. In this study, 21 Dlo-SR gene family members belonging to six subfamilies were identified, among which Dlo-RSZ20a, Dlo-SR30, Dlo-SR17, Dlo-SR53 and Dlo-SR32 were localized in the nucleus, Dlo-RSZ20b, Dlo-RSZ20c, Dlo-RSZ20d, Dlo-SC18, Dlo-RS2Z29, Dlo-SCL41, and Dlo-SR33 were localized in chloroplasts, and Dlo-RS43, Dlo-SC33, Dlo-SC37, Dlo-RS2Z33, Dlo-RS2Z16, Dlo-RS2Z24, Dlo-SCL43, Dlo-SR112, and Dlo-SR59 were localized in the nucleus and chloroplasts. The Dlo-SR genes exhibited differential expression patterns in different tissues of longan. The transcript levels of Dlo-RSZ20a, Dlo-SC18, Dlo-RS2Z29, DLo-SR59, Dlo-SR53, and Dlo-SR17 were low in all analyzed tissues, whereas Dlo-RS43, Dlo-RS2Z16, Dlo-RS2Z24, and Dlo-SR30 were highly expressed in all tissues. To clarify their function during SE, the transcript levels of Dlo-SR genes were analyzed at different four stages of SE, comprising non-embryonic callus (NEC), friable-embryogenic callus (EC), incomplete compact pro-embryogenic culture (ICpEC) and globular embryo (GE). Interestingly, the transcript levels of Dlo-RS2Z29 and Dlo-SR112 were increased in embryogenic cells compared with the NEC stage, whereas transcript levels of Dlo-RSZ20a, Dlo-RS43, Dlo-SC37, and Dlo-RS2Z16 were especially increased at the GE stage compared with the other stages. Alternative splicing events of Dlo-SR mRNA precursors (pre-mRNAs) was detected during SE, with totals of 41, 29, 35, and 44 events detected during NEC, EC, ICpEC, and GE respectively. Protein–protein interaction analysis showed that SR proteins were capable of interaction with each other. The results indicate that the alternative splicing of Dlo-SR pre-mRNAs occurs during SE and that Dlo-SR proteins may interact to regulate embryogenesis of longan.

Genome-wide identification and expression analyses of Sm genes reveal their involvement in early somatic embryogenesis in Dimocarpus longan Lour

The Sm proteins are a conserved protein family with Sm motifs. The family includes Sm and Sm-like proteins, which play important roles in pre-mRNA splicing. Most research on the Sm proteins have been conducted in herbaceous plants, and less in woody plants such as Dimocarpus longan (longan). And the embryo development status significantly affects the quality and yield of longan. In this study, we conducted a genome-wide analysis of longan Sm genes (DlSm) to clarify their roles during somatic embryogenesis (SE) and identified 29 Sm genes. Phylogenetic analysis deduced longan Sm proteins clustered into 17 phylogenetic groups with the homologous Sm proteins of Arabidopsis thaliana. We also analyzed the gene structures, motif compositions, and conserved domains of the longan Sm proteins. The promoter sequences of the DlSm genes contained many light, endosperm development, hormone, and temperature response elements, which suggested their possible functions. In the non-embryogenic callus(NEC) and during early SE in longan, the alternative splicing(AS) events of DlSm genes indicated that these genes may influence SE development by changing gene structures and sequences. The kinetin(KT) hormone, and blue and white light treatments affected the differentiation and growth of longan embryonic callus(EC) probably by affecting the AS events of DlSm genes. Expression profiles showed the possible functional divergence among Sm genes, and different hormones and light qualities affected their expression levels. The expression trends of the DlSm genes determined by RNA sequencing as fragments per kilobase of exon model per million mapped reads (FPKM) and by real-time quantitative PCR(qRT-PCR) during early SE in longan showed that the expression of the DlSm genes was affected by the growth and differentiation of longan SE, and decreased as the somatic embryo differentiation progressed. The results will contributed to understanding the longan Sm gene family and provide a basis for future functional validation studies.

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

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