Genome-wide analysis, transcription factor network approach and gene expression profile of GH3 genes over early somatic embryogenesis in Coffea spp

Integrated proteome and acetylome analyses provide novel insights into early somatic embryogenesis of Dimocarpus longan

Longan (Dimocarpus longan) is a precious subtropical fruit with high nutritional value. The somatic embryogenesis (SE) affects the quality and yield of fruit. Apart from clonal propagation, SE has extensive applications in genetic improvement and mutation. Thus, understanding the molecular basis of embryogenesis in longan will help to develop strategies for mass production of quality planting material. Lysine acetylation (Kac) plays an important role in diverse cellular processes, but limited knowledge is available regarding acetylation modifications in plant early SE. In this study, the proteome and acetylome of longan embryogenic callus (ECs) and globular embryos (GEs) were investigated. In total, 7232 proteins and 14,597 Kac sites were identified, and this resulted in the discovery of 1178 differentially expressed proteins and 669 differentially expressed acetylated proteins. KEGG and GO analysis showed that glucose metabolism, carbon metabolism, fatty acid degradation, and oxidative phosphorylation pathways were influenced by Kac modification. Furthermore, sodium butyrate (Sb, a deacetylase inhibitor) led to reduced the proliferation and delayed the differentiation of ECs by regulating the homeostasis of reactive oxygen species (ROS) andindole-3-acetic acid (IAA). Our study provides a comprehensive proteomic and acetylomic analysis to aid in understanding the molecular mechanisms involved in early SE, representing a potential tool for genetic improvement of longan.

Whole Genome Analysis of SLs Pathway Genes and Functional Characterization of DlSMXL6 in Longan Early Somatic Embryo Development

Strigolactones (SLs), a new class of plant hormones, are implicated in the regulation of various biological processes. However, the related family members and functions are not identified in longan (Dimocarpus longan Lour.). In this study, 23 genes in the CCD, D27, and SMXL family were identified in the longan genome. The phylogenetic relationships, gene structure, conserved motifs, promoter elements, and transcription factor-binding site predictions were comprehensively analysed. The expression profiles indicated that these genes may play important roles in longan organ development and abiotic stress responses, especially during early somatic embryogenesis (SE). Furthermore, GR24 (synthetic SL analogue) and Tis108 (SL biosynthesis inhibitor) could affect longan early SE by regulating the levels of endogenous IAA (indole-3-acetic acid), JA (jasmonic acid), GA (gibberellin), and ABA (abscisic acid). Overexpression of SMXL6 resulted in inhibition of longan SE by regulating the synthesis of SLs, carotenoids, and IAA levels. This study establishes a foundation for further investigation of SL genes and provides novel insights into their biological functions.

Embryogenic Stem Cell Identity after Protoplast Isolation from Daucus carota and Recovery of Regeneration Ability through Protoplast Culture

Protoplasts are single cells isolated from tissues or organs and are considered a suitable system for cell studies in plants. Embryogenic cells are totipotent stem cells, but their regeneration ability decreases or becomes lost altogether with extension of the culture period. In this study, we isolated and cultured EC-derived protoplasts (EC-pts) from carrots and compared them with non-EC-derived protoplasts (NEC-pts) with respect to their totipotency. The protoplast isolation conditions were optimized, and the EC-pts and NEC-pts were characterized by their cell size and types. Both types of protoplasts were then embedded using the alginate layer (TAL) method, and the resulting EC-pt-TALs and NEC-pt-TALs were cultured for further regeneration. The expression of the EC-specific genes SERK1, WUS, BBM, LEC1, and DRN was analyzed to confirm whether EC identity was maintained after protoplast isolation. The protoplast isolation efficiency for EC-pts was 2.4-fold higher than for NEC-pts (3.5 × 106 protoplasts·g−1 FW). In the EC-pt group, protoplasts < 20 µm accounted for 58% of the total protoplasts, whereas in the NEC-pt group, small protoplasts accounted for only 26%. In protoplast culture, the number of protoplasts that divided was 2.6-fold higher for EC-pts than for NEC-pts (7.7 × 104 protoplasts·g−1 FW), with a high number of plants regenerated for EC-pt-TALs, whereas no plants were induced by NEC-pt-TAL. Five times more plants were regenerated from EC-pts than from ECs. Regarding the expression of EC-specific genes, WUS and SERK1 expression increased 12-fold, and LEC1 and BBM expression increased 3.6−6.4-fold in isolated protoplasts compared with ECs prior to protoplast isolation (control). These results reveal that the protoplast isolation process did not affect the embryogenic cell identity; rather, it increased the plant regeneration rate, confirming that EC-derived protoplast culture may be an efficient system for increasing the regeneration ability of old EC cultures through the elimination of old and inactivate cells. EC-derived protoplasts may also represent an efficient single-cell system for application in new breeding technologies such as genome editing.

Insight into VvGH3 genes evolutional relationship from monocotyledons and dicotyledons reveals that VvGH3-9 negatively regulates the drought tolerance in transgenic Arabidopsis

The Gretchen Hagen3 (GH3) gene family is necessary for growth and development in plants and is regulated by osmotic stress and various hormones. Although it has been reported in many plants, the evolutionary relationship of GH3 in grape has not been systematically analyzed from the perspective of monocotyledonous and dicotyledonous. This study identified and analyzed 188 GH3 genes, which were distinctly divided into 9 subgroups, and found these subgroups have obviously been clustered between monocotyledonous and dicotyledonous. VvGH3-x genes had higher synteny with apple and Arabidopsis than that of rice, and the average Ka/Ks value in monocotyledons was higher than that of dicotyledons. The codon usage index showed that monocotyledons preferred to use G3s, C3s, and GC3s, while dicotyledons preferred to use A3s and T3s. The GH3 genes of grape exhibited different expression patterns in various tissues, different abiotic stresses, and hormonal treatments. The subcellular localization showed that VvGH3-9 was expressed in the nucleus and cytoplasm. Additionally, under 20% PEG treatment, the IAA and ABA contents, relative expression levels of VvGH3-9, relative electrical conductivity (REC), as well as MDA were obviously increased in VvGH3-9 overexpression lines at 72 h. In contrast, compared to WT, the contents of proline and H₂O₂, the activities of POD, SOD, and CAT, and the relative expression levels of drought responsive genes were significantly decreased in overexpressing lines. Collectively, this study provided helpful insight for the evolution of GH3 genes and presented some possibilities to study the functions of GH3 genes in monocotyledons and dicotyledons.

The SAUR gene family in coffee: genome-wide identification and gene expression analysis during somatic embryogenesis

BACKGROUND

Small auxin-up RNA (SAUR) genes form a wide family supposedly involved in different physiological and developmental processes in plants such as leaf senescence, auxin signaling and transport, hypocotyl development and tolerance to abiotic stresses. The transcription of SAUR genes is quickly induced by auxins, a group of phytohormones of major importance on embryo development. To better understand the distribution and expression profile of such still not explored family in Coffea sp., especially during the development of somatic embryogenesis (SE), SAUR members were characterized in silico using the available Coffea canephora genome data and analyzed for gene expression by RT-qPCR in C. arabica embryogenic samples.

METHODS AND RESULTS

Over C. canephora genome 31 CcSAURs were distributed by 11 chromosomes. Out of these 31 gene members, 5 SAURs were selected for gene expression analysis in C. arabica embryogenic materials. CaSAUR12 and CaSAUR18 were the members highly expressed through almost all plant materials. The other genes had more expression in at least one of the developing embryo stages or plantlets. The CaSAUR12 was the only member to exhibit an increased expression in both non-embryogenic calli and the developing embryo stages.

CONCLUSION

The identification of SAUR family on C. canephora genome followed by the analysis of gene expression profile across coffee somatic embryogenesis process on C. arabica represents a further additional step towards a better comprehension of molecular components acting on SE. Along with new research about this gene family such knowledge may support studies about clonal propagation methods via somatic embryogenesis to help the scientific community towards improvements into coffee crop.

Genome-Wide Analysis, Modeling, and Identification of Amino Acid Binding Motifs Suggest the Involvement of GH3 Genes during Somatic Embryogenesis of Coffea canephora

Auxin plays a central role in growth and plant development. To maintain auxin homeostasis, biological processes such as biosynthesis, transport, degradation, and reversible conjugation are essential. The Gretchen Hagen 3 (GH3) family genes codify for the enzymes that esterify indole-3-acetic acid (IAA) to various amino acids, which is a key process in the induction of somatic embryogenesis (SE). The GH3 family is one of the principal families of early response to auxin genes, exhibiting IAA-amido synthetase activity to maintain optimal levels of free auxin in the cell. In this study, we carried out a systematic identification of the GH3 gene family in the genome of Coffea canephora, determining a total of 18 CcGH3 genes. Analysis of the genetic structures and phylogenetic relationships of CcGH3 genes with GH3 genes from other plant species revealed that they could be clustered in two major categories with groups 1 and 2 of the GH3 family of Arabidopsis. We analyzed the transcriptome expression profiles of the 18 CcGH3 genes using RNA-Seq analysis-based data and qRT-PCR during the different points of somatic embryogenesis induction. Furthermore, the endogenous quantification of free and conjugated indole-3-acetic acid (IAA) suggests that the various members of the CcGH3 genes play a crucial role during the embryogenic process of C. canephora. Three-dimensional modeling of the selected CcGH3 proteins showed that they consist of two domains: an extensive N-terminal domain and a smaller C-terminal domain. All proteins analyzed in the present study shared a unique conserved structural topology. Additionally, we identified conserved regions that could function to bind nucleotides and specific amino acids for the conjugation of IAA during SE in C. canephora. These results provide a better understanding of the C. canephora GH3 gene family for further exploration and possible genetic manipulation.

Osmotic stress-induced somatic embryo maturation of coffee Coffea arabica L., shoot and root apical meristems development and robustness

Somatic embryogenesis (SE) is the most important plant biotechnology process for plant regeneration, propagation, genetic transformation and genome editing of coffee, Coffea arabica L. Somatic embryo (SEs) conversion to plantlets is the principal bottleneck for basic and applied use of this process. In this study we focus on the maturation of SEs of C. arabica var. Typica. SEs conversion to plantlet up to 95.9% was achieved under osmotic stress, using 9 g/L gelrite, as compared with only 39.34% in non-osmotic stress. Mature SEs induced in osmotic stress developed shoot and root apical meristems, while untreated SEs were unable to do it. C. arabica regenerated plants from osmotic stress were robust, with higher leaf and root area and internode length. To understand a possible regulatory mechanism, gene expression of key genes of C. arabica, homologous to sequences in the Arabidopsis thaliana genome, were analyzed. A set of two component system and cytokinin signaling-related coding genes (AHK1, AHK3, AHP4 and ARR1) which interact with WUSCHEL and WOX5 homedomains and morphogenic genes, BABY-BOOM, LEC1, FUS3 and AGL15, underwent significant changes during maturation of SEs of C. arabica var. Typica. This protocol is currently being applied in genetic transformation with high rate of success.

Transcriptome analysis identifies genes involved in the somatic embryogenesis of Eucalyptus

Background

Eucalyptus, a highly diverse genus of the Myrtaceae family, is the most widely planted hardwood in the world due to its increasing importance for fiber and energy. Somatic embryogenesis (SE) is one large-scale method to provide commercial use of the vegetative propagation of Eucalyptus and dedifferentiation is a key step for plant cells to become meristematic. However, little is known about the molecular changes during the Eucalyptus SE.

Results

We compared the transcriptome profiles of the differentiated and dedifferentiated tissues of two Eucalyptus species – E. camaldulensis (high embryogenetic potential) and E. grandis x urophylla (low embryogenetic potential). Initially, we identified 18,777 to 20,240 genes in all samples. Compared to the differentiated tissues, we identified 9229 and 8989 differentially expressed genes (DEGs) in the dedifferentiated tissues of E. camaldulensis and E. grandis x urophylla, respectively, and 2687 up-regulated and 2581 down-regulated genes shared. Next, we identified 2003 up-regulated and 1958 down-regulated genes only in E. camaldulensis, including 6 somatic embryogenesis receptor kinase, 17 ethylene, 12 auxin, 83 ribosomal protein, 28 zinc finger protein, 10 heat shock protein, 9 histone, 122 cell wall related and 98 transcription factor genes. Genes from other families like ABA, arabinogalactan protein and late embryogenesis abundant protein were also found to be specifically dysregulated in the dedifferentiation process of E. camaldulensis. Further, we identified 48,447 variants (SNPs and small indels) specific to E. camaldulensis, including 13,434 exonic variants from 4723 genes (e.g., annexin, GN, ARF and AP2-like ethylene-responsive transcription factor). qRT-PCR was used to confirm the gene expression patterns in both E. camaldulensis and E. grandis x urophylla.

Conclusions

This is the first time to study the somatic embryogenesis of Eucalyptus using transcriptome sequencing. It will improve our understanding of the molecular mechanisms of somatic embryogenesis and dedifferentiation in Eucalyptus. Our results provide a valuable resource for future studies in the field of Eucalyptus and will benefit the Eucalyptus breeding program.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07214-5.