Genome-wide analysis reveals divergent patterns of gene expression during zygotic and somatic embryo maturation of Theobroma cacao L., the chocolate tree

Genome-wide identification, characterization, and expression analysis of the small auxin-up RNA gene family during zygotic and somatic embryo maturation of the cacao tree (Theobroma cacao)

Small auxin-up RNA (SAUR) proteins were known as a large family that supposedly participated in various biological processes in higher plant species. However, the SAUR family has been still not explored in cacao (Theobroma cacao L.), one of the most important industrial trees. The present work, as an in silico study, revealed comprehensive aspects of the structure, phylogeny, and expression of TcSAUR gene family in cacao. A total of 90 members of the TcSAUR gene family have been identified and annotated in the cacao genome. According to the physic-chemical features analysis, all TcSAUR proteins exhibited slightly similar characteristics. Phylogenetic analysis showed that these TcSAUR proteins could be categorized into seven distinct groups, with 10 sub-groups. Our results suggested that tandemly duplication events, segmental duplication events, and whole genome duplication events might be important in the growth of the TcSAUR gene family in cacao. By re-analyzing the available transcriptome databases, we found that a number of TcSAUR genes were exclusively expressed during the zygotic embryogenesis and somatic embryogenesis. Taken together, our study will be valuable to further functional characterizations of candidate TcSAUR genes for the genetic engineering of cacao.

Single-base resolution methylomes of somatic embryogenesis in Theobroma cacao L. reveal epigenome modifications associated with somatic embryo abnormalities

Propagation by somatic embryogenesis in Theobroma cacao has some issues to be solved, as many morphologically abnormal somatic embryos that do not germinate into plants are frequently observed, thus hampering plant production on a commercial scale. For the first time the methylome landscape of T. cacao somatic embryogenesis was examined, using whole-genome bisulfite sequencing technique, with the aim to understand the epigenetic basis of somatic embryo abnormalities. We identified 873 differentially methylated genes (DMGs) in the CpG context between zygotic embryos, normal and abnormal somatic embryos, with important roles in development, programmed cell death, oxidative stress, and hypoxia induction, which can help to explain the morphological abnormalities of somatic embryos. We also identified the role of ethylene and its precursor 1-aminocyclopropane-1-carboxylate in several biological processes, such as hypoxia induction, cell differentiation and cell polarity, that could be associated to the development of abnormal somatic embryos. The biological processes and the hypothesis of ethylene and its precursor involvement in the somatic embryo abnormalities in cacao are discussed.

The GASA Gene Family in Theobroma cacao: Genome wide Identification and Expression Analysis.. [DOI: 10.1101/2021.01.27.425041]

The gibberellic acid-stimulated Arabidopsis (GASA/GAST) gene family is widely distributed in plants. The role of the GASA gene family has been reported previously in various physiological and biological processes, such as cell division, root and seed development, stem growth, and fruit ripening. These genes also provide resistance to abiotic and biotic stresses including antimicrobial, antiviral, and antifungal. Here, we report 17 tcGASA genes in Theobroma cacao L. distributed on six chromosomes. The gene structure, promoter-region sequences, protein structure, and biochemical properties, expression, and phylogenetics of all tcGASAs were analyzed. Phylogenetic analyses divided tcGASA proteins into five groups. The nine segmentally duplicating genes form four pairs and cluster together in phylogenetic tree. Purifying selection pressure was recorded on tcGASA, including duplicated genes. Several stress/hormone-responsive cis-regulatory elements were also recognized in the promoter region of tcGASAs. Differential expression analyses revealed that most of the tcGASA genes showed elevated expression in the seeds (cacao food), implying their role in seed development. The black rod disease of genus Phytophthora caused up to 20–25% loss (700,000 metric tons) in world cacao production. The role of tcGASA genes in conferring fungal resistance was also explored based on RNAseq data against Phytophthora megakarya. The differential expression of tcGASA genes was recorded between the tolerant and susceptible cultivars of cacao plants, which were inoculated with the fungus for 24h and 72h. This differential expression indicating possible role of tcGASA genes to fungal resistant in cacao. Our findings provide new insight into the function, evolution, and regulatory system of the GASA family genes in T. cacao and provide new target genes for development of fungi-resistant cacao varieties in breeding programs.

Attractor Concepts to Evaluate the Transcriptome-wide Dynamics Guiding Anaerobic to Aerobic State Transition in Escherichia coli

For any dynamical system, like living organisms, an attractor state is a set of variables or mechanisms that converge towards a stable system behavior despite a wide variety of initial conditions. Here, using multi-dimensional statistics, we investigate the global gene expression attractor mechanisms shaping anaerobic to aerobic state transition (AAT) of Escherichia coli in a bioreactor at early times. Out of 3,389 RNA-Seq expression changes over time, we identified 100 sharply changing genes that are key for guiding 1700 genes into the AAT attractor basin. Collectively, these genes were named as attractor genes constituting of 6 dynamic clusters. Apart from the expected anaerobic (glycolysis), aerobic (TCA cycle) and fermentation (succinate pathways) processes, sulphur metabolism, ribosome assembly and amino acid transport mechanisms together with 332 uncharacterised genes are also key for AAT. Overall, our work highlights the importance of multi-dimensional statistical analyses for revealing novel processes shaping AAT.

Transcriptional profiling reveals differentially expressed genes involved in lipid biosynthesis during cacao seed development

Theobroma cacao is a plant of economic value due to the use of its seed lipid for chocolate, confectionery, and cosmetic industries. The seed lipid contains a stable ratio of saturated and unsaturated fatty acids, which determines its unique melting temperature. However, little is known about the molecular mechanism determining the fatty acid ratio and lipid content in cacao. To gain insight into the unique properties of lipid synthesis in cacao, biochemical and transcriptomic approaches were used to compare the lipid accumulation between high and low lipid content cacao accessions. Lipid accumulation rates and lipid content were different between the two accessions. Moreover, differentially expressed genes were detected between high and low lipid content cacao accessions. The data allowed the identification of distinct candidate genes and furthered our understanding of lipid accumulation, potentially explaining the differences in lipid content between various cacao accessions. The results might be used to develop molecular tools and engineer alternative pathways for cacao breeding with improved lipid production potentials.

De novo transcriptome analysis of Tibetan medicinal plant Dysphania schraderiana

Dysphania schraderiana is widely distributed in Lhasa (Tibet, China) and used as a traditional medicine. However, the lack of genetic information hinders the understanding of its physiological processes, such as the biosynthesis of secondary metabolites. Herein, we used Illumina Hiseq4000 platform to sequence the transcriptome of flower and leaf tissues from D. schraderiana for the first time. Totally, 40,142 unigenes were assembled from approximately 5.2 million clean reads. All unigenes underwent gene prediction and were subsequently annotated in a NR (NCBI non-redundant protein) database, COG (Clusters of Orthologous Groups of proteins) database, and KEGG (Kyoto Encyclopedia of Genes and Genomes) database. Among the 40,142 unigenes, 2,579 genes were identified as differentially expressed between flowers and leaves, and used in further enrichment analysis. Also, 2,156 unigenes were annotated as transcription factors. Furthermore, our transcriptome analysis resulted in the identification of candidate unigenes annotated to enzymes involved in terpenoid biosynthesis. Taken together, this work has laid the foundation for the investigation of secondary metabolite biosynthesis and other physiological processes of D. schraderiana.

Metabolome and Transcriptome Association Analysis Reveals Dynamic Regulation of Purine Metabolism and Flavonoid Synthesis in Transdifferentiation during Somatic Embryogenesis in Cotton

Plant regeneration via somatic embryogenesis (SE) is a key step during genetic engineering. In the current study, integrated widely targeted metabolomics and RNA sequencing were performed to investigate the dynamic metabolic and transcriptional profiling of cotton SE. Our data revealed that a total of 581 metabolites were present in nonembryogenic staged calli (NEC), primary embryogenic calli (PEC), and initiation staged globular embryos (GE). Of the differentially accumulated metabolites (DAMs), nucleotides, and lipids were specifically accumulated during embryogenic differentiation, whereas flavones and hydroxycinnamoyl derivatives were accumulated during somatic embryo development. Additionally, metabolites related to purine metabolism were significantly enriched in PEC vs. NEC, whereas in GE vs. PEC, DAMs were remarkably associated with flavonoid biosynthesis. An association analysis of the metabolome and transcriptome data indicated that purine metabolism and flavonoid biosynthesis were co-mapped based on the Kyoto encyclopedia of genes and genomes (KEGG) database. Moreover, purine metabolism-related genes associated with signal recognition, transcription, stress, and lipid binding were significantly upregulated. Moreover, several classic somatic embryogenesis (SE) genes were highly correlated with their corresponding metabolites that were involved in purine metabolism and flavonoid biosynthesis. The current study identified a series of potential metabolites and corresponding genes responsible for SE transdifferentiation, which provides a valuable foundation for a deeper understanding of the regulatory mechanisms underlying cell totipotency at the molecular and biochemical levels.

Cocoa Bean Proteins-Characterization, Changes and Modifications due to Ripening and Post-Harvest Processing

The protein fractions of cocoa have been implicated influencing both the bioactive potential and sensory properties of cocoa and cocoa products. The objective of the present review is to show the impact of different stages of cultivation and processing with regard to the changes induced in the protein fractions. Special focus has been laid on the major seed storage proteins throughout the different stages of processing. The study starts with classical introduction of the extraction and the characterization methods used, while addressing classification approaches of cocoa proteins evolved during the timeline. The changes in protein composition during ripening and maturation of cocoa seeds, together with the possible modifications during the post-harvest processing (fermentation, drying, and roasting), have been documented. Finally, the bioactive potential arising directly or indirectly from cocoa proteins has been elucidated. The “state of the art” suggests that exploration of other potentially bioactive components in cocoa needs to be undertaken, while considering the complexity of reaction products occurring during the roasting phase of the post-harvest processing. Finally, the utilization of partially processed cocoa beans (e.g., fermented, conciliatory thermal treatment) can be recommended, providing a large reservoir of bioactive potentials arising from the protein components that could be instrumented in functionalizing foods.

Novel Avenues for Plant Protection: Plant Propagation by Somatic Embryogenesis Enhances Resistance to Insect Feeding

Somatic embryogenesis (SE), a clonal propagation method utilizing somatic cells, occurs under conditions that activate plant stress adaptation mechanisms such as production of protective secondary metabolites. Surprisingly, possible differences in susceptibility to insect pests between SE-generated and conventionally cultivated plants have not been previously explored. Here, we recorded frequencies and levels of bark-feeding damage by pine weevils (Hylobius abietis) in two large field trials, consisting of emblings (SE-propagated plants) and seedlings from 50 half-sib Norway spruce (Picea abies) families. We found that emblings were less frequently attacked by pine weevils, and when attacked, they were damaged to a lesser extent than seedlings. Moreover, we detected significant additive genetic variation in damage levels received by plants, indicating a heritable component to differences in resistance to insect herbivory among half-sib families. We present first-time evidence that emblings can be more resistant than seedlings to herbivorous insect damage, thus, SE appears to confer a previously unknown plant protection advantage. This finding indicates novel avenues to explore mechanisms underlying plant resistance and new approaches to develop non-toxic measures against insect pests.

The GhmiR157a-GhSPL10 regulatory module controls initial cellular dedifferentiation and callus proliferation in cotton by modulating ethylene-mediated flavonoid biosynthesis

MicroRNAs (miRNAs) modulate many biological processes through inactivation of specific mRNA targets such as those encoding transcription factors. A delicate spatial/temporal balance between specific miRNAs and their targets is central to achieving the appropriate biological outcomes. Somatic embryogenesis in cotton (Gossypium hirsutum), which goes through initial cellular dedifferentiation, callus proliferation, and somatic embryo development, is of great importance for both fundamental research and biotechnological applications. In this study, we characterize the function of the GhmiR157a-GhSPL10 miRNA-transcription factor module during somatic embryogenesis in cotton. We show that overexpression of GhSPL10, a target of GhmiR157a, increases free auxin and ethylene content and expression of associated signaling pathways, activates the flavonoid biosynthesis pathway, and promotes initial cellular dedifferentiation and callus proliferation. Inhibition of expression of the flavonoid synthesis gene F3H in GhSPL10 overexpression lines (35S:rSPL10-7) blocked callus initiation, while exogenous application of several types of flavonol promoted callus proliferation, associated with cell cycle-related gene expression. Inhibition of ethylene synthesis by aminoethoxyvinylglycine treatment in the 35S:rSPL10-7 line severely inhibited callus initiation, while activation of ethylene signaling through 1-aminocyclopropane 1-carboxylic acid treatment, EIN2 overexpression, or inhibition of the ethylene negative regulator CTR1 by RNA interference promoted flavonoid-related gene expression and flavonol accumulation. These results show that an up-regulation of ethylene signaling and the activation of flavonoid biosynthesis in GhSPL10 overexpression lines were associated with initial cellular dedifferentiation and callus proliferation. Our results demonstrate the importance of a GhmiR157a-GhSPL10 gene module in regulating somatic embryogenesis via hormonal and flavonoid pathways.

Cacao biotechnology: current status and future prospects

Theobroma cacao-The Food of the Gods, provides the raw material for the multibillion dollar chocolate industry and is also the main source of income for about 6 million smallholders around the world. Additionally, cocoa beans have a number of other nonfood uses in the pharmaceutical and cosmetic industries. Specifically, the potential health benefits of cocoa have received increasing attention as it is rich in polyphenols, particularly flavonoids. At present, the demand for cocoa and cocoa-based products in Asia is growing particularly rapidly and chocolate manufacturers are increasing investment in this region. However, in many Asian countries, cocoa production is hampered due to many reasons including technological, political and socio-economic issues. This review provides an overview of the present status of global cocoa production and recent advances in biotechnological applications for cacao improvement, with special emphasis on genetics/genomics, in vitro embryogenesis and genetic transformation. In addition, in order to obtain an insight into the latest innovations in the commercial sector, a survey was conducted on granted patents relating to T. cacao biotechnology.

Indirect somatic embryogenesis of Theobroma cacao L. in liquid medium and improvement of embryo-to-plantlet conversion rate

The establishment of cocoa embryogenic cell lines in liquid medium starting from high frequency somatic embryogenesis (HFSE) callus is described. The growth kinetics of the cultures during the multiplication and the expression steps conducted in 250 mL Erlenmeyer flasks were described for three genotypes selected for their agronomical traits (EET95, EET96, and EET103). The glucose and dissolved oxygen concentrations and the absorption of Murashige and Skoog medium macronutrients (nitrate, ammonium, potassium, sulfate, calcium, phosphorus, and magnesium) were monitored. The multiplication of the embryogenic calluses in a medium containing 2,4,5-trichlorophenoxyacetic acid (2,4,5-T) at 1 mg L^-1, initiated with an inoculation density of 20 g L^-1 of callus, was achieved. The growth rate was characterized by two phases, with the second being concomitant with a depletion of phosphorus and magnesium, and a decrease in the embryogenic potential of the callus. The expression of the callus embryogenic capacity was conducted in an auxin-free medium. The embryo production starting from 1 and 5 g L^-1 inoculation densities was compared. When placed in the optimal expression conditions in flasks, 1 g of callus produced 1000 to 1500 embryos within 5 to 7 wk. Finally, two paths for improving the plantlet regenerative capacities of cocoa SE produced in liquid medium were identified. Supplementing the expression medium with myo-inositol used as an osmotic agent at a concentration of 50 g L^-1 increased the embryo-to-plantlet conversion rate from 13-16% to 40-48%. A 6-wk culture of the embryos on a maturation medium in Petri dishes optimized their subsequent development into plantlets.

Transcription Factors: Their Role in the Regulation of Somatic Embryogenesis in Theobroma cacao L. and Other Species

Transcription factors are proteins that help with the control and regulation in the transcription of the DNA to mRNA by binding to special DNA sequences. With the aim to understand more about gene transcription regulation in Theobroma cacao L., this review outlines the principal transcription factors that were reported in other plants especially Arabidopsis thaliana and attempts at looking for the homologies with transcription factors in T. cacao. The information cited in this work is about the initiation, development, and maturation of the cacao somatic embryos and other crops. It is important to underline that there are very few publications in T. cacao discussing transcription factors that control the somatic embryogenesis process, but there is some information about transcription factors in other crops that we have used as a guide to try to understand this process.

Transcriptome Analysis of mRNA and miRNA in Somatic Embryos of Larix leptolepis Subjected to Hydrogen Treatment

Hydrogen is a therapeutic antioxidant that has been used extensively in clinical trials. It also acts as a bioactive molecule that can alleviate abiotic stress in plants. However, the biological effects of hydrogen in somatic embryos and the underlying molecular basis remain largely unknown. In this study, the morphological and physiological influence of exogenous H₂ treatment during somatic embryogenesis was characterized in Larix leptolepis Gordon. The results showed that exposure to hydrogen increased the proportions of active pro-embryogenic cells and normal somatic embryos. We sequenced mRNA and microRNA (miRNA) libraries to identify global transcriptome changes at different time points during H₂ treatment of larch pro-embryogenic masses (PEMs). A total of 45,393 mRNAs and 315 miRNAs were obtained. Among them, 4253 genes and 96 miRNAs were differentially expressed in the hydrogen-treated libraries compared with the control. Further, a large number of the differentially expressed mRNAs and miRNAs were related to reactive oxygen species (ROS) homeostasis and cell cycle regulation. We also identified 4399 potential target genes for 285 of the miRNAs. The differential expression data and the mRNA-miRNA interaction network described here provide new insights into the molecular mechanisms that determine the performance of PEMs exposed to H₂ during somatic embryogenesis.

Protocorms and Protocorm-Like Bodies Are Molecularly Distinct from Zygotic Embryonic Tissues in Phalaenopsis aphrodite

The distinct reproductive program of orchids provides a unique evolutionary model with pollination-triggered ovule development and megasporogenesis, a modified embryogenesis program resulting in seeds with immature embryos, and mycorrhiza-induced seed germination. However, the molecular mechanisms that have evolved to establish these unparalleled developmental programs are largely unclear. Here, we conducted comparative studies of genome-wide gene expression of various reproductive tissues and captured the molecular events associated with distinct reproductive programs in Phalaenopsis aphrodite Importantly, our data provide evidence to demonstrate that protocorm-like body (PLB) regeneration (the clonal regeneration practice used in the orchid industry) does not follow the embryogenesis program. Instead, we propose that SHOOT MERISTEMLESS, a class I KNOTTED-LIKE HOMEOBOX gene, is likely to play a role in PLB regeneration. Our studies challenge the current understanding of the embryonic identity of PLBs. Taken together, the data obtained establish a fundamental framework for orchid reproductive development and provide a valuable new resource to enable the prediction of gene regulatory networks that is required for specialized developmental programs of orchid species.

Theobroma cacao L. pathogenesis-related gene tandem array members show diverse expression dynamics in response to pathogen colonization

Background

The pathogenesis-related (PR) group of proteins are operationally defined as polypeptides that increase in concentration in plant tissues upon contact with a pathogen. To date, 17 classes of highly divergent proteins have been described that act through multiple mechanisms of pathogen resistance. Characterizing these families in cacao, an economically important tree crop, and comparing the families to those in other species, is an important step in understanding cacao’s immune response.

Results

Using publically available resources, all members of the 17 recognized pathogenesis-related gene families in the genome of Theobroma cacao were identified and annotated resulting in a set of ~350 members in both published cacao genomes. Approximately 50 % of these genes are organized in tandem arrays scattered throughout the genome. This feature was observed in five additional plant taxa (three dicots and two monocots), suggesting that tandem duplication has played an important role in the evolution of the PR genes in higher plants. Expression profiling captured the dynamics and complexity of PR genes expression at basal levels and after induction by two cacao pathogens (the oomycete, Phytophthora palmivora, and the fungus, Colletotrichum theobromicola), identifying specific genes within families that are more responsive to pathogen challenge. Subsequent qRT-PCR validated the induction of several PR-1, PR-3, PR-4, and PR-10 family members, with greater than 1000 fold induction detected for specific genes.

Conclusions

We describe candidate genes that are likely to be involved in cacao’s defense against Phytophthora and Colletotrichum infection and could be potentially useful for marker-assisted selection for breeding of disease resistant cacao varieties. The data presented here, along with existing cacao—omics resources, will enable targeted functional genetic screening of defense genes likely to play critical functions in cacao’s defense against its pathogens.

Electronic supplementary material

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

Enhanced somatic embryogenesis in Theobroma cacao using the homologous BABY BOOM transcription factor

BACKGROUND

Theobroma cacao, the chocolate tree, is an important economic crop in East Africa, South East Asia, and South and Central America. Propagation of elite varieties has been achieved through somatic embryogenesis (SE) but low efficiencies and genotype dependence still presents a significant limitation for its propagation at commercial scales. Manipulation of transcription factors has been used to enhance the formation of SEs in several other plant species. This work describes the use of the transcription factor Baby Boom (BBM) to promote the transition of somatic cacao cells from the vegetative to embryonic state.

RESULTS

An ortholog of the Arabidopsis thaliana BBM gene (AtBBM) was characterized in T. cacao (TcBBM). TcBBM expression was observed throughout embryo development and was expressed at higher levels during SE as compared to zygotic embryogenesis (ZE). TcBBM overexpression in A. thaliana and T. cacao led to phenotypes associated with SE that did not require exogenous hormones. While transient ectopic expression of TcBBM provided only moderate enhancements in embryogenic potential, constitutive overexpression dramatically increased SE proliferation but also appeared to inhibit subsequent development.

CONCLUSION

Our work provides validation that TcBBM is an ortholog to AtBBM and has a specific role in both somatic and zygotic embryogenesis. Furthermore, our studies revealed that TcBBM transcript levels could serve as a biomarker for embryogenesis in cacao tissue. Results from transient expression of TcBBM provide confirmation that transcription factors can be used to enhance SE without compromising plant development and avoiding GMO plant production. This strategy could compliment a hormone-based method of reprogramming somatic cells and lead to more precise manipulation of SE at the regulatory level of transcription factors. The technology would benefit the propagation of elite varieties with low regeneration potential as well as the production of transgenic plants, which similarly requires somatic cell reprogramming.