Towards the understanding of the cocoa transcriptome: Production and analysis of an exhaustive dataset of ESTs of Theobroma cacao L. generated from various tissues and under various conditions

Using ddRADseq to assess the genetic diversity of in-farm and gene bank cacao resources in the Baracoa region, eastern Cuba, for use and conservation purposes

The Baracoa region, eastern Cuba, hosts around 80 % of the country cacao (Theobroma cacao L.) plantations. Cacao plants in farms are diverse in origin and propagation, with grafted and hybrid plants being the more common ones. Less frequent are plants from cuttings, TSH progeny, and traditional Cuban cacao. A national cacao gene bank is also present in Baracoa, with 282 accessions either prospected in Cuba or introduced from other countries. A breeding program associated with the gene bank started in the 1990s based on agro-morphological descriptors. The genetic diversity of cacao resources in Baracoa has been poorly described, except for traditional Cuban cacao, affecting the proper development of the breeding program and the cacao planting policies in the region. To assess the population structure and genetic diversity of cacao resources in Baracoa region, we genotyped plants from both cacao gene bank (CG) and cacao farms (CF) applying a new ddRADseq protocol for cacao. After data processing, two SNPs datasets containing 11,425 and 6,481 high-quality SNPs were generated with 238 CG and 135 CF plants, respectively. SNPs were unevenly distributed along the 10 cacao chromosomes and laid mainly in noncoding regions of the genome. Population structure analysis with these SNP datasets identified seven and four genetic groups in CG and CF samples, respectively. Clustering using UPGMA and principal component analysis mostly agree with population structure results. Amelonado was the predominant cacao ancestry, accounting for 49.22 % (CG) and 57.73 % (CF) of the total. Criollo, Contamana, Iquitos, and Nanay ancestries were detected in both CG and CF samples, while Nacional and Marañon backgrounds were only identified in CG. Genetic differentiation among CG (FST ranging from 0.071 to 0.407) was higher than among CF genetic groups (FST: 0.093-0.282). Genetic diversity parameters showed similar values for CG and CF samples. The CG and CF genetic groups with the lowest genetic diversity parameters had the highest proportion of Amelonado ancestry. These results should contribute to reinforcing the ongoing breeding program and updating the planting policies on cacao farms, with an impact on the social and economic life of the region.

Cacao pod transcriptome profiling of seven genotypes identifies features associated with post-penetration resistance to Phytophthora palmivora

The oomycete Phytophthora palmivora infects the fruit of cacao trees (Theobroma cacao) causing black pod rot and reducing yields. Cacao genotypes vary in their resistance levels to P. palmivora, yet our understanding of how cacao fruit respond to the pathogen at the molecular level during disease establishment is limited. To address this issue, disease development and RNA-Seq studies were conducted on pods of seven cacao genotypes (ICS1, WFT, Gu133, Spa9, CCN51, Sca6 and Pound7) to better understand their reactions to the post-penetration stage of P. palmivora infection. The pod tissue-P. palmivora pathogen assay resulted in the genotypes being classified as susceptible (ICS1, WFT, Gu133 and Spa9) or resistant (CCN51, Sca6 and Pound7). The number of differentially expressed genes (DEGs) ranged from 1625 to 6957 depending on genotype. A custom gene correlation approach identified 34 correlation groups. De novo motif analysis was conducted on upstream promoter sequences of differentially expressed genes, identifying 76 novel motifs, 31 of which were over-represented in the upstream sequences of correlation groups and associated with gene ontology terms related to oxidative stress response, defense against fungal pathogens, general metabolism and cell function. Genes in one correlation group (Group 6) were strongly induced in all genotypes and enriched in genes annotated with defense-responsive terms. Expression pattern profiling revealed that genes in Group 6 were induced to higher levels in the resistant genotypes. An additional analysis allowed the identification of 17 candidate cis-regulatory modules likely to be involved in cacao defense against P. palmivora. This study is a comprehensive exploration of the cacao pod transcriptional response to P. palmivora spread after infection. We identified cacao genes, promoter motifs, and promoter motif combinations associated with post-penetration resistance to P. palmivora in cacao pods and provide this information as a resource to support future and ongoing efforts to breed P. palmivora-resistant cacao.

Genome-wide association studies and genomic selection assays made in a large sample of cacao (Theobroma cacao L.) germplasm reveal significant marker-trait associations and good predictive value for improving yield potential

A genome-wide association study (GWAS) was undertaken to unravel marker-trait associations (MTAs) between SNP markers and phenotypic traits. It involved a subset of 421 cacao accessions from the large and diverse collection conserved ex situ at the International Cocoa Genebank Trinidad. A Mixed Linear Model (MLM) in TASSEL was used for the GWAS and followed by confirmatory analyses using GAPIT FarmCPU. An average linkage disequilibrium (r2) of 0.10 at 5.2 Mb was found across several chromosomes. Seventeen significant (P ≤ 8.17 × 10-5 (-log10 (p) = 4.088)) MTAs of interest, including six that pertained to yield-related traits, were identified using TASSEL MLM. The latter accounted for 5 to 17% of the phenotypic variation expressed. The highly significant association (P ≤ 8.17 × 10-5) between seed length to width ratio and TcSNP 733 on chromosome 5 was verified with FarmCPU (P ≤ 1.12 × 10-8). Fourteen MTAs were common to both the TASSEL and FarmCPU models at P ≤ 0.003. The most significant yield-related MTAs involved seed number and seed length on chromosome 7 (P ≤ 1.15 × 10-14 and P ≤ 6.75 × 10-05, respectively) and seed number on chromosome 1 (P ≤ 2.38 × 10-05), based on the TASSEL MLM. It was noteworthy that seed length, seed length to width ratio and seed number were associated with markers at different loci, indicating their polygenic nature. Approximately 40 candidate genes that encode embryo and seed development, protein synthesis, carbohydrate transport and lipid biosynthesis and transport were identified in the flanking regions of the significantly associated SNPs and in linkage disequilibrium with them. A significant association of fruit surface anthocyanin intensity co-localised with MYB-related protein 308 on chromosome 4. Testing of a genomic selection approach revealed good predictive value (genomic estimated breeding values (GEBV)) for economic traits such as seed number (GEBV = 0.611), seed length (0.6199), seed width (0.5435), seed length to width ratio (0.5503), seed/cotyledon mass (0.6014) and ovule number (0.6325). The findings of this study could facilitate genomic selection and marker-assisted breeding of cacao thereby expediting improvement in the yield potential of cacao planting material.

TcTI, a Kunitz-type trypsin inhibitor from cocoa associated with defense against pathogens

Protease inhibitors (PIs) are important biotechnological tools of interest in agriculture. Usually they are the first proteins to be activated in plant-induced resistance against pathogens. Therefore, the aim of this study was to characterize a Theobroma cacao trypsin inhibitor called TcTI. The ORF has 740 bp encoding a protein with 219 amino acids, molecular weight of approximately 23 kDa. rTcTI was expressed in the soluble fraction of Escherichia coli strain Rosetta [DE3]. The purified His-Tag rTcTI showed inhibitory activity against commercial porcine trypsin. The kinetic model demonstrated that rTcTI is a competitive inhibitor, with a Ki value of 4.08 × 10-7 mol L-1. The thermostability analysis of rTcTI showed that 100% inhibitory activity was retained up to 60 °C and that at 70-80 °C, inhibitory activity remained above 50%. Circular dichroism analysis indicated that the protein is rich in loop structures and β-conformations. Furthermore, in vivo assays against Helicoverpa armigera larvae were also performed with rTcTI in 0.1 mg mL-1 spray solutions on leaf surfaces, which reduced larval growth by 70% compared to the control treatment. Trials with cocoa plants infected with Mp showed a greater accumulation of TcTI in resistant varieties of T. cacao, so this regulation may be associated with different isoforms of TcTI. This inhibitor has biochemical characteristics suitable for biotechnological applications as well as in resistance studies of T. cacao and other crops.

A Non-Shedding Fruit Elaeis oleifera Palm Reveals Perturbations to Hormone Signaling, ROS Homeostasis, and Hemicellulose Metabolism

The developmentally programmed loss of a plant organ is called abscission. This process is characterized by the ultimate separation of adjacent cells in the abscission zone (AZ). The discovery of an American oil palm (Elaeis oleifera) variant that does not shed its has allowed for the study of the mechanisms of ripe fruit abscission in this species. A comparative transcriptome analysis was performed to compare the fruit AZs of the non-shedding E. oleifera variant to an individual of the same progeny that sheds its ripe fruit normally. The study provides evidence for widespread perturbation to gene expression in the AZ of the non-shedding variant, compared to the normal fruit-shedding control, and offers insight into abscission-related functions. Beyond the genes with known or suspected roles during organ abscission or indehiscence that were identified, a list of genes with hormone-related functions, including ethylene, jasmonic acid, abscisic acid, cytokinin and salicylic acid, in addition to reactive oxygen species (ROS) metabolism, transcriptional responses and signaling pathways, was compiled. The results also allowed a comparison between the ripe fruit abscission processes of the African and American oil palm species at the molecular level and revealed commonalities with environmental stress pathways.

RNA-Seq analysis reveals potential regulators of programmed cell death and leaf remodelling in lace plant (Aponogeton madagascariensis)

BACKGROUND

The lace plant (Aponogeton madagascariensis) is an aquatic monocot that develops leaves with uniquely formed perforations through the use of a developmentally regulated process called programmed cell death (PCD). The process of perforation formation in lace plant leaves is subdivided into several developmental stages: pre-perforation, window, perforation formation, perforation expansion and mature. The first three emerging "imperforate leaves" do not form perforations, while all subsequent leaves form perforations via developmentally regulated PCD. PCD is active in cells called "PCD cells" that do not retain the antioxidant anthocyanin in spaces called areoles framed by the leaf veins of window stage leaves. Cells near the veins called "NPCD cells" retain a red pigmentation from anthocyanin and do not undergo PCD. While the cellular changes that occur during PCD are well studied, the gene expression patterns underlying these changes and driving PCD during leaf morphogenesis are mostly unknown. We sought to characterize differentially expressed genes (DEGs) that mediate lace plant leaf remodelling and PCD. This was achieved performing gene expression analysis using transcriptomics and comparing DEGs among different stages of leaf development, and between NPCD and PCD cells isolated by laser capture microdissection.

RESULTS

Transcriptomes were sequenced from imperforate, pre-perforation, window, and mature leaf stages, as well as PCD and NPCD cells isolated from window stage leaves. Differential expression analysis of the data revealed distinct gene expression profiles: pre-perforation and window stage leaves were characterized by higher expression of genes involved in anthocyanin biosynthesis, plant proteases, expansins, and autophagy-related genes. Mature and imperforate leaves upregulated genes associated with chlorophyll development, photosynthesis, and negative regulators of PCD. PCD cells were found to have a higher expression of genes involved with ethylene biosynthesis, brassinosteroid biosynthesis, and hydrolase activity whereas NPCD cells possessed higher expression of auxin transport, auxin signalling, aspartyl proteases, cysteine protease, Bag5, and anthocyanin biosynthesis enzymes.

CONCLUSIONS

RNA sequencing was used to generate a de novo transcriptome for A. madagascariensis leaves and revealed numerous DEGs potentially involved in PCD and leaf remodelling. The data generated from this investigation will be useful for future experiments on lace plant leaf development and PCD in planta.

Convergent Biochemical Pathways for Xanthine Alkaloid Production in Plants Evolved from Ancestral Enzymes with Different Catalytic Properties

Convergent evolution is widespread but the extent to which common ancestral conditions are necessary to facilitate the independent acquisition of similar traits remains unclear. In order to better understand how ancestral biosynthetic catalytic capabilities might lead to convergent evolution of similar modern-day biochemical pathways, we resurrected ancient enzymes of the caffeine synthase (CS) methyltransferases that are responsible for theobromine and caffeine production in flowering plants. Ancestral CS enzymes of Theobroma, Paullinia, and Camellia exhibited similar substrate preferences but these resulted in the formation of different sets of products. From these ancestral enzymes, descendants with similar substrate preference and product formation independently evolved after gene duplication events in Theobroma and Paullinia. Thus, it appears that the convergent modern-day pathways likely originated from ancestral pathways with different inferred flux. Subsequently, the modern-day enzymes originated independently via gene duplication and their convergent catalytic characteristics evolved to partition the multiple ancestral activities by different mutations that occurred in homologous regions of the ancestral proteins. These results show that even when modern-day pathways and recruited genes are similar, the antecedent conditions may be distinctive such that different evolutionary steps are required to generate convergence.

Pathogenicity of Moniliophthora roreri isolates from selected morphology groups in harvested cacao pods and in vitro sensitivity to compost tea

Moniliopthora roreri is the frosty pod rot disease (FPD) and one of the most devastating cacao pathogens worldwide. However, M. roreri pathogenicity on harvested cacao pods and sensitivity to compost tea have not been fully described. Monosporic cultures of M. roreri from different morphology groups were obtained. The isolates’ pathogenicity was tested by inoculation onto harvested cacao pods, and symptoms were evaluated at 3-day intervals during 16 days before estimating the area under the disease progress curve (AUDPC). The sensitivity of M. roreri to compost tea was evaluated on potato dextrose agar (PDA) amended with 1 to 5 % compost tea. All morphology groups could infect harvested cacao pods during the 16 days with a disease severity index abode 75 %. Compost tea completely inhibited the growth of M. roreri when used at 4.5 % or higher. Results suggest a shortened biotrophic phase during the infection in harvested pods and a medium to high sensitivity of M. roreri to compost tea.

Multi-scale comparative transcriptome analysis reveals key genes and metabolic reprogramming processes associated with oil palm fruit abscission

BACKGROUND

Fruit abscission depends on cell separation that occurs within specialized cell layers that constitute an abscission zone (AZ). To determine the mechanisms of fleshy fruit abscission of the monocot oil palm (Elaeis guineensis Jacq.) compared with other abscission systems, we performed multi-scale comparative transcriptome analyses on fruit targeting the developing primary AZ and adjacent tissues.

RESULTS

Combining between-tissue developmental comparisons with exogenous ethylene treatments, and naturally occurring abscission in the field, RNAseq analysis revealed a robust core set of 168 genes with differentially regulated expression, spatially associated with the ripe fruit AZ, and temporally restricted to the abscission timing. The expression of a set of candidate genes was validated by qRT-PCR in the fruit AZ of a natural oil palm variant with blocked fruit abscission, which provides evidence for their functions during abscission. Our results substantiate the conservation of gene function between dicot dry fruit dehiscence and monocot fleshy fruit abscission. The study also revealed major metabolic transitions occur in the AZ during abscission, including key senescence marker genes and transcriptional regulators, in addition to genes involved in nutrient recycling and reallocation, alternative routes for energy supply and adaptation to oxidative stress.

CONCLUSIONS

The study provides the first reference transcriptome of a monocot fleshy fruit abscission zone and provides insight into the mechanisms underlying abscission by identifying key genes with functional roles and processes, including metabolic transitions, cell wall modifications, signalling, stress adaptations and transcriptional regulation, that occur during ripe fruit abscission of the monocot oil palm. The transcriptome data comprises an original reference and resource useful towards understanding the evolutionary basis of this fundamental plant process.

A Combined RNA Preservation and Extraction Protocol for Gene Expression Studies in Cacao Beans

Despite the high economic importance of cacao beans, few RNA-based studies have been conducted on this plant material and hence no optimal RNA-extraction has been reported. Moreover, extraction of high-quality RNA from recalcitrant cacao bean tissue has shown many difficulties and requires optimization. Furthermore, cacao beans are mostly found at remote and under-resourced locations, which pressures the outsourcing of such analysis and thereby demands RNA-stable preservation and transportation of cacao beans. This study aims to select an appropriate RNA extraction and preservation/transportation method for cacao beans. For this purpose, three sample homogenization and five extraction protocols on cacao beans were compared. In addition, 13 preservation conditions-differing in tissue crushing degree, preservation method, duration, and temperature-were compared and evaluated. A comparative analysis revealed that CTAB-based homogenization and extraction outcompeted all tested commercial protocols in RNA yield and integrity, respectively. Preservation at -80°C affected RNA quality the least, whereas freeze-drying was most suitable for transportation at room temperature for maximum 1 week. The cacao bean RNA obtained from the selected methods were compatible for downstream applications. The results of this study will facilitate on-field sampling and transportation of genetically sensitive cacao material prior to cacao bean transcriptomic studies. In addition, valuable insights on sample homogenization, extraction, preservation, and transportation have been provided, which is of interest to every plant geneticist.

Adaptation of hydroxymethylbutenyl diphosphate reductase enables volatile isoprenoid production

Volatile isoprenoids produced by plants are emitted in vast quantities into the atmosphere, with substantial effects on global carbon cycling. Yet, the molecular mechanisms regulating the balance between volatile and non-volatile isoprenoid production remain unknown. Isoprenoids are synthesised via sequential condensation of isopentenyl pyrophosphate (IPP) to dimethylallyl pyrophosphate (DMAPP), with volatile isoprenoids containing fewer isopentenyl subunits. The DMAPP:IPP ratio could affect the balance between volatile and non-volatile isoprenoids, but the plastidic DMAPP:IPP ratio is generally believed to be similar across different species. Here we demonstrate that the ratio of DMAPP:IPP produced by hydroxymethylbutenyl diphosphate reductase (HDR/IspH), the final step of the plastidic isoprenoid production pathway, is not fixed. Instead, this ratio varies greatly across HDRs from phylogenetically distinct plants, correlating with isoprenoid production patterns. Our findings suggest that adaptation of HDR plays a previously unrecognised role in determining in vivo carbon availability for isoprenoid emissions, directly shaping global biosphere-atmosphere interactions.

The pathogen Moniliophthora perniciosa promotes differential proteomic modulation of cacao genotypes with contrasting resistance to witches´ broom disease

BACKGROUND

Witches' broom disease (WBD) of cacao (Theobroma cacao L.), caused by Moniliophthora perniciosa, is the most important limiting factor for the cacao production in Brazil. Hence, the development of cacao genotypes with durable resistance is the key challenge for control the disease. Proteomic methods are often used to study the interactions between hosts and pathogens, therefore helping classical plant breeding projects on the development of resistant genotypes. The present study compared the proteomic alterations between two cacao genotypes standard for WBD resistance and susceptibility, in response to M. perniciosa infection at 72 h and 45 days post-inoculation; respectively the very early stages of the biotrophic and necrotrophic stages of the cacao x M. perniciosa interaction.

RESULTS

A total of 554 proteins were identified, being 246 in the susceptible Catongo and 308 in the resistant TSH1188 genotypes. The identified proteins were involved mainly in metabolism, energy, defense and oxidative stress. The resistant genotype showed more expressed proteins with more variability associated with stress and defense, while the susceptible genotype exhibited more repressed proteins. Among these proteins, stand out pathogenesis related proteins (PRs), oxidative stress regulation related proteins, and trypsin inhibitors. Interaction networks were predicted, and a complex protein-protein interaction was observed. Some proteins showed a high number of interactions, suggesting that those proteins may function as cross-talkers between these biological functions.

CONCLUSIONS

We present the first study reporting the proteomic alterations of resistant and susceptible genotypes in the T. cacao x M. perniciosa pathosystem. The important altered proteins identified in the present study are related to key biologic functions in resistance, such as oxidative stress, especially in the resistant genotype TSH1188, that showed a strong mechanism of detoxification. Also, the positive regulation of defense and stress proteins were more evident in this genotype. Proteins with significant roles against fungal plant pathogens, such as chitinases, trypsin inhibitors and PR 5 were also identified, and they may be good resistance markers. Finally, important biological functions, such as stress and defense, photosynthesis, oxidative stress and carbohydrate metabolism were differentially impacted with M. perniciosa infection in each genotype.

Development of a core SNP panel for cacao (Theobroma cacao L.) identity analysis

Single nucleotide polymorphisms (SNPs) are preferred markers for DNA fingerprinting and diversity studies in cacao (Theobroma cacao L.). Yet, a consensus SNP panel with a minimum number of SNPs for optimal identity analysis is unavailable for cacao. An initial set of 146 SNP panels of varying sizes were assembled based on heterozygosity, linkage disequilibrium (LD), linkage group (LG) distribution, major allele frequency, minor allele frequency (M_iAF), polymorphism information content (PIC), and random distribution. These panels were assessed to determine their ability to distinguish among a training set of 155 accessions. The panels with the best separation ability were supplemented with additional SNPs to create 16 designer panels, which separated all 155 accessions. The 16 designer SNP panels were then assessed on a dataset of 1220 accessions coming from 10 ancestral groups. Increasing the number of SNPs generally yielded improved resolution of genetic identities with concomitant reduction of synonymous groups. The number and choice of SNPs were critical factors with LD, M_iAF, and PIC being important selection attributes but an even LG distribution was unnecessary. A robust set of 96 SNPs is recommended as a minimal core SNP panel for cacao DNA fingerprinting to the international cacao community.

Draft genome sequence of fastidious pathogen Ceratobasidium theobromae, which causes vascular-streak dieback in Theobroma cacao

Background

Ceratobasidium theobromae, a member of the Ceratobasidiaceae family, is the causal agent of vascular-streak dieback (VSD) of cacao, a major threat to the chocolate industry in the South-East Asia. The fastidious pathogen is very hard to isolate and maintain in pure culture, which is a major bottleneck in the study of its genetic diversity and genome.

Result

This study describes for the first time, a 33.90 Mbp de novo assembled genome of a putative C. theobromae isolate from cacao. Ab initio gene prediction identified 9264 protein-coding genes, of which 800 are unique to C. theobromae when compared to Rhizoctonia spp., a closely related group. Transcriptome analysis using RNA isolated from 4 independent VSD symptomatic cacao stems identified 3550 transcriptionally active genes when compared to the assembled C. theobromae genome while transcripts for only 4 C. theobromae genes were detected in 2 asymptomatic stems. De novo assembly of the non-cacao associated reads from the VSD symptomatic stems uniformly produced genes with high identity to predicted genes in the C. theobromae genome as compared to Rhizoctonia spp. or genes found in Genbank. Further analysis of the predicted C. theobromae transcriptome was carried out identifying CAZy gene classes, KEGG-pathway associated genes, and 138 putative effector proteins.

Conclusion

These findings put forth, for the first time, a predicted genome for the fastidious basidiomycete C. theobromae causing VSD on cacao providing a model for testing and comparison in the future. The C. theobromae genome predicts a pathogenesis model involving secreted effector proteins to suppress plant defense mechanisms and plant cell wall degrading enzymes.

Assessing hidden parentage and genetic integrity of the "United Fruit Clones" of cacao (Theobroma cacao) from Costa Rica using SNP markers

The international cacao collection in CATIE, Costa Rica contains nearly 1200 accessions of cacao, mainly from the center of genetic diversity of this species. Among these accessions, the United Fruit clones (UF clones) were developed by the United Fruit Company in Costa Rica, and they represent one of the earliest groups of improved cacao germplasm in the world. Some of these UF clones have been used as key progenitors for breeding resistance/tolerance to Frosty Pod and Black Pod diseases in the Americas. Accurate information on the identity and background of these clones is important for their effective use in breeding. Using Single Nucleotide Polymorphism (SNP) markers, we genotyped 273 cacao germplasm accessions including 44 UF clones and 229 reference accessions. We verified the true-to-type identity of UF clones in the CATIE cacao collection and analyzed their population memberships using maximum-likelihood-based approaches. Three duplicate groups, representing approximately 30% of the UF clones, were identified. Both distance- and model-based clustering methods showed that the UF clones were mainly composed of Trinitario, ancient Nacional and hybrids between ancient Nacional and Amelonado. This result filled the information gap about the UF clones thus will improve their utilization for cacao breeding.

The use and domestication of Theobroma cacao during the mid-Holocene in the upper Amazon

Cacao (Theobroma cacao L.) is an important economic crop, yet studies of its domestication history and early uses are limited. Traditionally, cacao is thought to have been first domesticated in Mesoamerica. However, genomic research shows that T. cacao's greatest diversity is in the upper Amazon region of northwest South America, pointing to this region as its centre of origin. Here, we report cacao use identified by three independent lines of archaeological evidence-cacao starch grains, absorbed theobromine residues and ancient DNA-dating from approximately 5,300 years ago recovered from the Santa Ana-La Florida (SALF) site in southeast Ecuador. To our knowledge, these findings constitute the earliest evidence of T. cacao use in the Americas and the first unequivocal archaeological example of its pre-Columbian use in South America. They also reveal the upper Amazon region as the oldest centre of cacao domestication yet identified.

Transcriptome Profiles of Nod Factor-independent Symbiosis in the Tropical Legume Aeschynomene evenia

Nod factors (NF) were assumed to be indispensable for the establishment of a rhizobium-legume symbiosis until the discovery that certain Bradyrhizobium strains interacting with certain Aeschynomene species lack the canonical nodABC genes required for their synthesis. So far, the molecular dialogue between Aeschynomene and its symbionts remains an open question. Here we report a time course transcriptional analysis of Aeschynomene evenia in response to inoculation with Bradyrhizobium ORS278. The NF-independent symbiotic process was monitored at five time points between bacterial infection and nodule maturity. The five time points correspond to three specific events, root infection by crack entry, nodule organogenesis, and the establishment of the nitrogen fixing process. During the third stage, about 80 NCR-like genes and eight symbiotic genes known to be involved in signaling, bacterial infection or nodulation regulation were highly expressed. Comparative gene expression analyses at the five time points also enabled the selection of genes with an expression profile that makes them promising markers to monitor early plant responses to bacteria. Such markers could be used in bioassays to identify the nature of the bacterial signal(s). Our data represent valuable resources for investigation of this Nod factor-independent symbiosis.

Transcriptome analysis of embryonic domains in Norway spruce reveals potential regulators of suspensor cell death

The terminal differentiation and elimination of the embryo-suspensor is the earliest manifestation of programmed cell death (PCD) during plant ontogenesis. Molecular regulation of suspensor PCD remains poorly understood. Norway spruce (Picea abies) embryos provide a powerful model for studying embryo development because of their large size, sequenced genome, and the possibility to obtain a large number of embryos at a specific developmental stage through somatic embryogenesis. Here, we have carried out global gene expression analysis of the Norway spruce embryo-suspensor versus embryonal mass (a gymnosperm analogue of embryo proper) using RNA sequencing. We have identified that suspensors have enhanced expression of the NAC domain-containing transcription factors, XND1 and ANAC075, previously shown to be involved in the initiation of developmental PCD in Arabidiopsis. The analysis has also revealed enhanced expression of Norway spruce homologues of the known executioners of both developmental and stress-induced cell deaths, such as metacaspase 9 (MC9), cysteine endopeptidase-1 (CEP1) and ribonuclease 3 (RNS3). Interestingly, a spruce homologue of bax inhibitor-1 (PaBI-1, for Picea abies BI-1), an evolutionarily conserved cell death suppressor, was likewise up-regulated in the embryo-suspensor. Since Arabidopsis BI-1 so far has been implicated only in the endoplasmic reticulum (ER)-stress induced cell death, we investigated its role in embryogenesis and suspensor PCD using RNA interference (RNAi). We have found that PaBI-1-deficient lines formed a large number of abnormal embryos with suppressed suspensor elongation and disturbed polarity. Cytochemical staining of suspensor cells has revealed that PaBI-1 deficiency suppresses vacuolar cell death and induces necrotic type of cell death previously shown to compromise embryo development. This study demonstrates that a large number of cell-death components are conserved between angiosperms and gymnosperms and establishes a new role for BI-1 in the progression of vacuolar cell death.

Convergent evolution of caffeine in plants by co-option of exapted ancestral enzymes

Convergent evolution is a process that has occurred throughout the tree of life, but the historical genetic and biochemical context promoting the repeated independent origins of a trait is rarely understood. The well-known stimulant caffeine, and its xanthine alkaloid precursors, has evolved multiple times in flowering plant history for various roles in plant defense and pollination. We have shown that convergent caffeine production, surprisingly, has evolved by two previously unknown biochemical pathways in chocolate, citrus, and guaraná plants using either caffeine synthase- or xanthine methyltransferase-like enzymes. However, the pathway and enzyme lineage used by any given plant species is not predictable from phylogenetic relatedness alone. Ancestral sequence resurrection reveals that this convergence was facilitated by co-option of genes maintained over 100 million y for alternative biochemical roles. The ancient enzymes of the Citrus lineage were exapted for reactions currently used for various steps of caffeine biosynthesis and required very few mutations to acquire modern-day enzymatic characteristics, allowing for the evolution of a complete pathway. Future studies aimed at manipulating caffeine content of plants will require the use of different approaches given the metabolic and genetic diversity revealed by this study.

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.

Colombia a Source of Cacao Genetic Diversity As Revealed by the Population Structure Analysis of Germplasm Bank of Theobroma cacao L

Beans of the species Theobroma cacao L., also known as cacao, are the raw material to produce chocolate. Colombian cacao has been classified as a fine flavor cacao that represents the 5% of cacao world's production. Colombian genetic resources from this species are conserved in ex situ and in-field germplasm banks, since T. cacao has recalcitrant seeds to desication and long-term storage. Currently, the collection of T. cacao of the Colombian Corporation of Agricultural Research (CORPOICA) has approximately 700 germplasm accessions. We conducted a molecular analysis of Corpoica's cacao collection and a morphological characterization of some accessions with the goal to study its genetic diversity and population structure and, to select interesting accessions for the cacao's breeding program. Phenotypic evaluation was performed based on 18 morphological traits and 4 biochemical traits. PCA analysis of morphological traits explained 60.6% of the total variation in seven components and 100% of the total variation of biochemical traits in four components, grouping the collection in 4 clusters for both variables. We explored 565 accessions from Corpoica's germplasm and 252 accessions from reference populations using 96 single nucleotide polymorphism (SNP) molecular markers. Molecular patterns of cacao Corpoica's collection were obtained amplifying specific alleles in a Fluidigm platform that used integrated circuits of fluids. Corpoica's collection showed highest genetic diversity [Expected Heterozygosity (HE = 0.314), Observed Heterozygosity (HO = 0.353)] that is reduced when reference populations were included in the dataset (HE = 0.294, HO = 0.261). The collection was divided into four clusters based on population structure analysis. Cacao accessions from distinct groups showed some taxonomic concordance and reflected their geographic origins. For instance, accessions classified as Criollo were clearly differentiated in one group and we identified two new Colombian genetic groups. Using a number of allelic variations based on 87 SNP markers and 22 different morphological/biochemical traits, a core collection with a total of 232 accessions was selected as a primary genetic resource for cacao breeders.

Expression of the Theobroma cacao Bax-inhibitor-1 gene in tomato reduces infection by the hemibiotrophic pathogen Moniliophthora perniciosa

Programmed cell death (PCD) plays a key role in plant responses to pathogens, determining the success of infection depending on the pathogen lifestyle and on which participant of the interaction triggers cell death. The hemibiotrophic basidiomycete Moniliophthora perniciosa is the causal agent of witches' broom disease of Theobroma cacao L. (cacao), a serious constraint for production in South America and the Caribbean. It has been hypothesized that M. perniciosa pathogenesis involves PCD, initially as a plant defence mechanism, which is diverted by the fungus to induce necrosis during the dikaryotic phase of the mycelia. Here, we evaluated whether the expression of a cacao anti-apoptotic gene would affect the incidence and severity of M. perniciosa infection using the 'Micro-Tom' (MT) tomato as a model. The cacao Bax-inhibitor-1 (TcBI-1) gene, encoding a putative basal attenuator of PCD, was constitutively expressed in MT to evaluate function. Transformants expressing TcBI-1, when treated with tunicamycin, an inducer of endoplasmic reticulum stress, showed a decrease in cell peroxidation. When the same transformants were inoculated with the necrotrophic fungal pathogens Sclerotinia sclerotiorum, Sclerotium rolfsii and Botrytis cinerea, a significant reduction in infection severity was observed, confirming TcBI-1 function. After inoculation with M. perniciosa, TcBI-1 transformant lines showed a significant reduction in disease incidence compared with MT. The overexpression of TcBI-1 appears to affect the ability of germinating spores to penetrate susceptible tissues, restoring part of the non-host resistance in MT against the S-biotype of M. perniciosa.

The cacao Criollo genome v2.0: an improved version of the genome for genetic and functional genomic studies

Background

Theobroma cacao L., native to the Amazonian basin of South America, is an economically important fruit tree crop for tropical countries as a source of chocolate. The first draft genome of the species, from a Criollo cultivar, was published in 2011. Although a useful resource, some improvements are possible, including identifying misassemblies, reducing the number of scaffolds and gaps, and anchoring un-anchored sequences to the 10 chromosomes.

Methods

We used a NGS-based approach to significantly improve the assembly of the Belizian Criollo B97-61/B2 genome. We combined four Illumina large insert size mate paired libraries with 52x of Pacific Biosciences long reads to correct misassembled regions and reduced the number of scaffolds. We then used genotyping by sequencing (GBS) methods to increase the proportion of the assembly anchored to chromosomes.

Results

The scaffold number decreased from 4,792 in assembly V1 to 554 in V2 while the scaffold N50 size has increased from 0.47 Mb in V1 to 6.5 Mb in V2. A total of 96.7% of the assembly was anchored to the 10 chromosomes compared to 66.8% in the previous version. Unknown sites (Ns) were reduced from 10.8% to 5.7%. In addition, we updated the functional annotations and performed a new RefSeq structural annotation based on RNAseq evidence.

Conclusion

Theobroma cacao Criollo genome version 2 will be a valuable resource for the investigation of complex traits at the genomic level and for future comparative genomics and genetics studies in cacao tree. New functional tools and annotations are available on the Cocoa Genome Hub (http://cocoa-genome-hub.southgreen.fr).

Electronic supplementary material

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

Phytophthora megakarya and P. palmivora, Causal Agents of Black Pod Rot, Induce Similar Plant Defense Responses Late during Infection of Susceptible Cacao Pods

Phytophthora megakarya (Pmeg) and Phytophthora palmivora (Ppal) cause black pod rot of Theobroma cacao L. (cacao). Of these two clade 4 species, Pmeg is more virulent and is displacing Ppal in many cacao production areas in Africa. Symptoms and species specific sporangia production were compared when the two species were co-inoculated onto pod pieces in staggered 24 h time intervals. Pmeg sporangia were predominantly recovered from pod pieces with unwounded surfaces even when inoculated 24 h after Ppal. On wounded surfaces, sporangia of Ppal were predominantly recovered if the two species were simultaneously applied or Ppal was applied first but not if Pmeg was applied first. Pmeg demonstrated an advantage over Ppal when infecting un-wounded surfaces while Ppal had the advantage when infecting wounded surfaces. RNA-Seq was carried out on RNA isolated from control and Pmeg and Ppal infected pod pieces 3 days post inoculation to assess their abilities to alter/suppress cacao defense. Expression of 4,482 and 5,264 cacao genes was altered after Pmeg and Ppal infection, respectively, with most genes responding to both species. Neural network self-organizing map analyses separated the cacao RNA-Seq gene expression profiles into 24 classes, 6 of which were largely induced in response to infection. Using KEGG analysis, subsets of genes composing interrelated pathways leading to phenylpropanoid biosynthesis, ethylene and jasmonic acid biosynthesis and action, plant defense signal transduction, and endocytosis showed induction in response to infection. A large subset of genes encoding putative Pr-proteins also showed differential expression in response to infection. A subset of 36 cacao genes was used to validate the RNA-Seq expression data and compare infection induced gene expression patterns in leaves and wounded and unwounded pod husks. Expression patterns between RNA-Seq and RT-qPCR were generally reproducible. The level and timing of altered gene expression was influenced by the tissues studied and by wounding. Although, in these susceptible interactions gene expression patterns were similar, some genes did show differential expression in a Phytophthora species dependent manner. The biggest difference was the more intense changes in expression in Ppal inoculated wounded pod pieces further demonstrating its rapid progression when penetrating through wounds.

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.

Antimicrobial properties of two novel peptides derived from Theobroma cacao osmotin

The osmotin proteins of several plants display antifungal activity, which can play an important role in plant defense against diseases. Thus, this protein can be useful as a source for biotechnological strategies aiming to combat fungal diseases. In this work, we analyzed the antifungal activity of a cacao osmotin-like protein (TcOsm1) and of two osmotin-derived synthetic peptides with antimicrobial features, differing by five amino acids residues at the N-terminus. Antimicrobial tests showed that TcOsm1 expressed in Escherichia coli inhibits the growth of Moniliophthora perniciosa mycelium and Pichia pastoris X-33 in vitro. The TcOsm1-derived peptides, named Osm-pepA (H-RRLDRGGVWNLNVNPGTTGARVWARTK-NH2), located at R23-K49, and Osm-pepB (H-GGVWNLNVNPGTTGARVWARTK-NH2), located at G28-K49, inhibited growth of yeasts (Saccharomyces cerevisiae S288C and Pichia pastoris X-33) and spore germination of the phytopathogenic fungi Fusarium f. sp. glycines and Colletotrichum gossypi. Osm-pepA was more efficient than Osm-pepB for S. cerevisiae (MIC=40μM and MIC=127μM, respectively), as well as for P. pastoris (MIC=20μM and MIC=127μM, respectively). Furthermore, the peptides presented a biphasic performance, promoting S. cerevisiae growth in doses around 5μM and inhibiting it at higher doses. The structural model for these peptides showed that the five amino acids residues, RRLDR at Osm-pepA N-terminus, significantly affect the tertiary structure, indicating that this structure is important for the peptide antimicrobial potency. This is the first report of development of antimicrobial peptides from T. cacao. Taken together, the results indicate that the cacao osmotin and its derived peptides, herein studied, are good candidates for developing biotechnological tools aiming to control phytopathogenic fungi.

Cloning, Expression Profiling and Functional Analysis of CnHMGS, a Gene Encoding 3-hydroxy-3-Methylglutaryl Coenzyme A Synthase from Chamaemelum nobile

Roman chamomile (Chamaemelum nobile L.) is renowned for its production of essential oils, which major components are sesquiterpenoids. As the important enzyme in the sesquiterpenoid biosynthesis pathway, 3-hydroxy-3-methylglutaryl coenzyme A synthase (HMGS) catalyze the crucial step in the mevalonate pathway in plants. To isolate and identify the functional genes involved in the sesquiterpene biosynthesis of C. nobile L., a HMGS gene designated as CnHMGS (GenBank Accession No. KU529969) was cloned from C. nobile. The cDNA sequence of CnHMGS contained a 1377 bp open reading frame encoding a 458-amino-acid protein. The sequence of the CnHMGS protein was highly homologous to those of HMGS proteins from other plant species. Phylogenetic tree analysis revealed that CnHMGS clustered with the HMGS of Asteraceae in the dicotyledon clade. Further functional complementation of CnHMGS in the mutant yeast strain YSC6274 lacking HMGS activity demonstrated that the cloned CnHMGS cDNA encodes a functional HMGS. Transcript profile analysis indicated that CnHMGS was preferentially expressed in flowers and roots of C. nobile. The expression of CnHMGS could be upregulated by exogenous elicitors, including methyl jasmonate and salicylic acid, suggesting that CnHMGS was elicitor-responsive. The characterization and expression analysis of CnHMGS is helpful to understand the biosynthesis of sesquiterpenoid in C. nobile at the molecular level and also provides molecular wealth for the biotechnological improvement of this important medicinal plant.

First Microsatellite Markers Developed from Cupuassu ESTs: Application in Diversity Analysis and Cross-Species Transferability to Cacao

The cupuassu tree (Theobroma grandiflorum) (Willd. ex Spreng.) Schum. is a fruitful species from the Amazon with great economical potential, due to the multiple uses of its fruit´s pulp and seeds in the food and cosmetic industries, including the production of cupulate, an alternative to chocolate. In order to support the cupuassu breeding program and to select plants presenting both pulp/seed quality and fungal disease resistance, SSRs from Next Generation Sequencing ESTs were obtained and used in diversity analysis. From 8,330 ESTs, 1,517 contained one or more SSRs (1,899 SSRs identified). The most abundant motifs identified in the EST-SSRs were hepta- and trinucleotides, and they were found with a minimum and maximum of 2 and 19 repeats, respectively. From the 1,517 ESTs containing SSRs, 70 ESTs were selected based on their functional annotation, focusing on pulp and seed quality, as well as resistance to pathogens. The 70 ESTs selected contained 77 SSRs, and among which, 11 were polymorphic in cupuassu genotypes. These EST-SSRs were able to discriminate the cupuassu genotype in relation to resistance/susceptibility to witches’ broom disease, as well as to pulp quality (SST/ATT values). Finally, we showed that these markers were transferable to cacao genotypes, and that genome availability might be used as a predictive tool for polymorphism detection and primer design useful for both Theobroma species. To our knowledge, this is the first report involving EST-SSRs from cupuassu and is also a pioneer in the analysis of marker transferability from cupuassu to cacao. Moreover, these markers might contribute to develop or saturate the cupuassu and cacao genetic maps, respectively.

TcCYPR04, a Cacao Papain-Like Cysteine-Protease Detected in Senescent and Necrotic Tissues Interacts with a Cystatin TcCYS4

The interaction amongst papain-like cysteine-proteases (PLCP) and their substrates and inhibitors, such as cystatins, can be perceived as part of the molecular battlefield in plant-pathogen interaction. In cacao, four cystatins were identified and characterized by our group. We identified 448 proteases in cacao genome, whereof 134 were cysteine-proteases. We expressed in Escherichia coli a PLCP from cacao, named TcCYSPR04. Immunoblottings with anti-TcCYSPR04 exhibited protein increases during leaf development. Additional isoforms of TcCYSPR04 appeared in senescent leaves and cacao tissues infected by Moniliophthora perniciosa during the transition from the biotrophic to the saprophytic phase. TcCYSPR04 was induced in the apoplastic fluid of Catongo and TSH1188 cacao genotypes, susceptible and resistant to M. perniciosa, respectively, but greater intensity and additional isoforms were observed in TSH1188. The fungal protein MpNEP induced PLCP isoform expression in tobacco leaves, according to the cross reaction with anti-TcCYSPR04. Several protein isoforms were detected at 72 hours after treatment with MpNEP. We captured an active PLCP from cacao tissues, using a recombinant cacao cystatin immobilized in CNBr-Sepharose. Mass spectrometry showed that this protein corresponds to TcCYSPR04. A homology modeling was obtained for both proteins. In order to become active, TcCYSPR04 needs to lose its inhibitory domain. Molecular docking showed the physical-chemical complementarities of the interaction between the cacao enzyme and its inhibitor. We propose that TcCYSPR04 and its interactions with cacao cystatins are involved in the senescence and necrosis events related to witches' broom symptoms. This molecular interaction may be the target for future interventions to control witches' broom disease.

Two Theobroma cacao genotypes with contrasting pathogen tolerance show aberrant transcriptional and ROS responses after salicylic acid treatment

Understanding the genetic basis of pathogen susceptibility in various crop plants is crucial to increasing the stability of food, feed, and fuel production. Varietal differences in defence responses provide insights into the mechanisms of resistance and are a key resource for plant breeders. To explore the role of salicylic acid in the regulation of defence in cacao, we demonstrated that SA treatment decreased susceptibility to a pod rot pathogen, Phytophthora tropicalis in two genotypes, Scavina 6 and Imperial College Selection 1, which differ in their resistance to several agriculturally important pathogens. Transient overexpression of TcNPR1, a major transcriptional regulator of the SA-dependent plant immune system, also increased pathogen tolerance in cacao leaves. To explore further the genetic basis of resistance in cacao, we used microarrays to measure gene expression profiles after salicylic acid (SA) treatment in these two cacao genotypes. The two genotypes displayed distinct transcriptional responses to SA. Unexpectedly, the expression profile of the susceptible genotype ICS1 included a larger number of pathogenesis-related genes that were induced by SA at 24h after treatment, whereas genes encoding many chloroplast and mitochondrial proteins implicated in reactive oxygen species production were up-regulated in the resistant genotype, Sca6. Sca6 accumulated significantly more superoxide at 24h after treatment of leaves with SA. These experiments revealed critical insights regarding the molecular differences between cacao varieties, which will allow a better understanding of defence mechanisms to help guide breeding programmes.

Tc-MYBPA an Arabidopsis TT2-like transcription factor and functions in the regulation of proanthocyanidin synthesis in Theobroma cacao

BACKGROUND

The flavan-3-ols catechin and epicatechin, and their polymerized oligomers, the proanthocyanidins (PAs, also called condensed tannins), accumulate to levels of up to 15 % of the total weight of dry seeds of Theobroma cacao L. These compounds have been associated with several health benefits in humans. They also play important roles in pest and disease defense throughout the plant. In Arabidopsis, the R2R3 type MYB transcription factor TT2 regulates the major genes leading to the synthesis of PA.

RESULTS

To explore the transcriptional regulation of the PA synthesis pathway in cacao, we isolated and characterized an R2R3 type MYB transcription factor MYBPA from cacao. We examined the spatial and temporal gene expression patterns of the Tc-MYBPA gene and found it to be developmentally expressed in a manner consistent with its involvement in PAs and anthocyanin synthesis. Functional complementation of an Arabidopsis tt2 mutant with Tc-MYBPA suggested that it can functionally substitute the Arabidopsis TT2 gene. Interestingly, in addition to PA accumulation in seeds of the Tc-MYBPA expressing plants, we also observed an obvious increase of anthocyanidin accumulation in hypocotyls. We observed that overexpression of the Tc-MYBPA gene resulted in increased expression of several key genes encoding the major structural enzymes of the PA and anthocyanidin pathway, including DFR (dihydroflavanol reductase), LDOX (leucoanthocyanidin dioxygenase) and BAN (ANR, anthocyanidin reductase).

CONCLUSION

We conclude that the Tc-MYBPA gene that encodes an R2R3 type MYB transcription factor is an Arabidopsis TT2 like transcription factor, and may be involved in the regulation of both anthocyanin and PA synthesis in cacao. This research may provide molecular tools for breeding of cacao varieties with improved disease resistance and enhanced flavonoid profiles for nutritional and pharmaceutical applications.

Making a chocolate chip: development and evaluation of a 6K SNP array for Theobroma cacao

Theobroma cacao, the key ingredient in chocolate production, is one of the world's most important tree fruit crops, with ∼4,000,000 metric tons produced across 50 countries. To move towards gene discovery and marker-assisted breeding in cacao, a single-nucleotide polymorphism (SNP) identification project was undertaken using RNAseq data from 16 diverse cacao cultivars. RNA sequences were aligned to the assembled transcriptome of the cultivar Matina 1-6, and 330,000 SNPs within coding regions were identified. From these SNPs, a subset of 6,000 high-quality SNPs were selected for inclusion on an Illumina Infinium SNP array: the Cacao6kSNP array. Using Cacao6KSNP array data from over 1,000 cacao samples, we demonstrate that our custom array produces a saturated genetic map and can be used to distinguish among even closely related genotypes. Our study enhances and expands the genetic resources available to the cacao research community, and provides the genome-scale set of tools that are critical for advancing breeding with molecular markers in an agricultural species with high genetic diversity.

Analysis of gene expression and proteomic profiles of clonal genotypes from Theobroma cacao subjected to soil flooding

Soil flooding causes changes in gene transcription, synthesis and degradation of proteins and cell metabolism. The main objective of this study was to understand the biological events of Theobroma cacao during soil flooding-induced stress, using the analyses of gene expression and activity of key enzymes involved in fermentation, as well as the identification of differentially expressed proteins by mass spectrometry in two contrasting genotypes for flooding tolerance (tolerant - TSA-792 and susceptible - TSH-774). Soil anoxia caused by flooding has led to changes in the expression pattern of genes associated with the biosynthesis of alcohol dehydrogenase (ADH), pyruvate decarboxylase (PDC) and lactate dehydrogenase (LDH) in leaves and roots of the two evaluated genotypes. Significant differences were observed between the enzyme activities of the two genotypes. Leaves and roots of the TSA-792 genotype showed higher ADH activity as compared to the TSH-774 genotype, whereas the activities of PDC and LDH have varied over the 96 h of soil flooding, being higher for TSA-792 genotype, at the initial stage, and TSH-774 genotype, at the final stage. Some of the identified proteins are those typical of the anaerobic metabolism-involved in glycolysis and alcoholic fermentation-and different proteins associated with photosynthesis, protein metabolism and oxidative stress. The ability to maintain glycolysis and induce fermentation was observed to play an important role in anoxia tolerance in cacao and may also serve to distinguish tolerant and susceptible genotypes in relation to this stressor.

Pervasive effects of a dominant foliar endophytic fungus on host genetic and phenotypic expression in a tropical tree

It is increasingly recognized that macro-organisms (corals, insects, plants, vertebrates) consist of both host tissues and multiple microbial symbionts that play essential roles in their host's ecological and evolutionary success. Consequently, identifying benefits and costs of symbioses, as well as mechanisms underlying them are research priorities. All plants surveyed under natural conditions harbor foliar endophytic fungi (FEF) in their leaf tissues, often at high densities. Despite producing no visible effects on their hosts, experiments have nonetheless shown that FEF reduce pathogen and herbivore damage. Here, combining results from three genomic, and two physiological experiments, we demonstrate pervasive genetic and phenotypic effects of the apparently asymptomatic endophytes on their hosts. Specifically, inoculation of endophyte-free (E−) Theobroma cacao leaves with Colletotrichum tropicale (E+), the dominant FEF species in healthy T. cacao, induces consistent changes in the expression of hundreds of host genes, including many with known defensive functions. Further, E+ plants exhibited increased lignin and cellulose content, reduced maximum rates of photosynthesis (A_max), and enrichment of nitrogen-15 and carbon-13 isotopes. These phenotypic changes observed in E+ plants correspond to changes in expression of specific functional genes in related pathways. Moreover, a cacao gene (Tc00g04254) highly up-regulated by C. tropicale also confers resistance to pathogen damage in the absence of endophytes or their products in host tissues. Thus, the benefits of increased pathogen resistance in E+ plants are derived in part from up-regulation of intrinsic host defense responses, and appear to be offset by potential costs including reduced photosynthesis, altered host nitrogen metabolism, and endophyte heterotrophy of host tissues. Similar effects are likely in most plant-endophyte interactions, and should be recognized in the design and interpretation of genetic and phenotypic studies of plants.

Sequence and expression analyses of ethylene response factors highly expressed in latex cells from Hevea brasiliensis

The AP2/ERF superfamily encodes transcription factors that play a key role in plant development and responses to abiotic and biotic stress. In Hevea brasiliensis, ERF genes have been identified by RNA sequencing. This study set out to validate the number of HbERF genes, and identify ERF genes involved in the regulation of latex cell metabolism. A comprehensive Hevea transcriptome was improved using additional RNA reads from reproductive tissues. Newly assembled contigs were annotated in the Gene Ontology database and were assigned to 3 main categories. The AP2/ERF superfamily is the third most represented compared with other transcription factor families. A comparison with genomic scaffolds led to an estimation of 114 AP2/ERF genes and 1 soloist in Hevea brasiliensis. Based on a phylogenetic analysis, functions were predicted for 26 HbERF genes. A relative transcript abundance analysis was performed by real-time RT-PCR in various tissues. Transcripts of ERFs from group I and VIII were very abundant in all tissues while those of group VII were highly accumulated in latex cells. Seven of the thirty-five ERF expression marker genes were highly expressed in latex. Subcellular localization and transactivation analyses suggested that HbERF-VII candidate genes encoded functional transcription factors.

Expansion of banana (Musa acuminata) gene families involved in ethylene biosynthesis and signalling after lineage-specific whole-genome duplications

Whole-genome duplications (WGDs) are widespread in plants, and three lineage-specific WGDs occurred in the banana (Musa acuminata) genome. Here, we analysed the impact of WGDs on the evolution of banana gene families involved in ethylene biosynthesis and signalling, a key pathway for banana fruit ripening. Banana ethylene pathway genes were identified using comparative genomics approaches and their duplication modes and expression profiles were analysed. Seven out of 10 banana ethylene gene families evolved through WGD and four of them (1-aminocyclopropane-1-carboxylate synthase (ACS), ethylene-insensitive 3-like (EIL), ethylene-insensitive 3-binding F-box (EBF) and ethylene response factor (ERF)) were preferentially retained. Banana orthologues of AtEIN3 and AtEIL1, two major genes for ethylene signalling in Arabidopsis, were particularly expanded. This expansion was paralleled by that of EBF genes which are responsible for control of EIL protein levels. Gene expression profiles in banana fruits suggested functional redundancy for several MaEBF and MaEIL genes derived from WGD and subfunctionalization for some of them. We propose that EIL and EBF genes were co-retained after WGD in banana to maintain balanced control of EIL protein levels and thus avoid detrimental effects of constitutive ethylene signalling. In the course of evolution, subfunctionalization was favoured to promote finer control of ethylene signalling.

Improving palm oil quality through identification and mapping of the lipase gene causing oil deterioration

The oil palm fruit mesocarp contains high lipase activity that increases free fatty acids and necessitates post-harvest inactivation by heat treatment of fruit bunches. Even before heat treatment the mesocarp lipase activity causes consequential oil losses and requires costly measures to limit free fatty acids quantities. Here we demonstrate that elite low-lipase lines yield oil with substantially less free fatty acids than standard genotypes, allowing more flexibility for post-harvest fruit processing and extended ripening for increased yields. We identify the lipase and its gene cosegregates with the low-/high-lipase trait, providing breeders a marker to rapidly identify potent elite genitors and introgress the trait into major cultivars. Overall, economic gains brought by wide adoption of this material could represent up to one billion dollars per year. Expected benefits concern all planters but are likely to be highest for African smallholders who would be more able to produce oil that meets international quality standards.

Lipase in the mesocarp of oil palm fruits reduces the quality of oil and leads to decreased yield. In this study, the authors identify the gene encoding the mesocarp lipase, raising the possibility of breeding palm genotypes that produce better quality oil.

Accurate determination of genetic identity for a single cacao bean, using molecular markers with a nanofluidic system, ensures cocoa authentication

Cacao (Theobroma cacao L.), the source of cocoa, is an economically important tropical crop. One problem with the premium cacao market is contamination with off-types adulterating raw premium material. Accurate determination of the genetic identity of single cacao beans is essential for ensuring cocoa authentication. Using nanofluidic single nucleotide polymorphism (SNP) genotyping with 48 SNP markers, we generated SNP fingerprints for small quantities of DNA extracted from the seed coat of single cacao beans. On the basis of the SNP profiles, we identified an assumed adulterant variety, which was unambiguously distinguished from the authentic beans by multilocus matching. Assignment tests based on both Bayesian clustering analysis and allele frequency clearly separated all 30 authentic samples from the non-authentic samples. Distance-based principle coordinate analysis further supported these results. The nanofluidic SNP protocol, together with forensic statistical tools, is sufficiently robust to establish authentication and to verify gourmet cacao varieties. This method shows significant potential for practical application.

TcNPR3 from Theobroma cacao functions as a repressor of the pathogen defense response

BACKGROUND

Arabidopsis thaliana (Arabidopsis) NON-EXPRESSOR OF PR1 (NPR1) is a transcription coactivator that plays a central role in regulating the transcriptional response to plant pathogens. Developing flowers of homozygous npr3 mutants are dramatically more resistant to infection by the pathogenic bacterium Pseudomonas syringae, suggesting a role of NPR3 as a repressor of NPR1-mediated defense response with a novel role in flower development.

RESULTS

We report here the characterization of a putative NPR3 gene from the tropical tree species Theobroma cacao (TcNPR3). Like in Arabidopsis, TcNPR3 was constitutively expressed across a wide range of tissue types and developmental stages but with some differences in relative levels compared to Arabidopsis. To test the function of TcNPR3, we performed transgenic complementation analysis by introducing a constitutively expressing putative TcNPR3 transgene into an Arabidopsis npr3 mutant. TcNPR3 expressing Arabidopsis plants were partially restored to the WT pathogen phenotype (immature flowers susceptible to bacterial infection). To test TcNPR3 function directly in cacao tissues, a synthetic microRNA targeting TcNPR3 mRNA was transiently expressed in cacao leaves using an Agrobacterium-infiltration method. TcNPR3 knock down leaf tissues were dramatically more resistance to infection with Phytophthora capsici in a leaf bioassay, showing smaller lesion sizes and reduced pathogen replication.

CONCLUSIONS

We conclude that TcNPR3 functions similar to the Arabidopsis NPR3 gene in the regulation of the cacao defense response. Since TcNPR3 did not show a perfect complementation of the Arabidopsis NPR3 mutation, the possibility remains that other functions of TcNPR3 remain to be found. This novel knowledge can contribute to the breeding of resistant cacao varieties against pathogens through molecular markers based approaches or biotechnological strategies.

Proanthocyanidin synthesis in Theobroma cacao: genes encoding anthocyanidin synthase, anthocyanidin reductase, and leucoanthocyanidin reductase

BACKGROUND

The proanthocyanidins (PAs), a subgroup of flavonoids, accumulate to levels of approximately 10% total dry weight of cacao seeds. PAs have been associated with human health benefits and also play important roles in pest and disease defense throughout the plant.

RESULTS

To dissect the genetic basis of PA biosynthetic pathway in cacao (Theobroma cacao), we have isolated three genes encoding key PA synthesis enzymes, anthocyanidin synthase (ANS), anthocyanidin reductase (ANR) and leucoanthocyanidin reductase (LAR). We measured the expression levels of TcANR, TcANS and TcLAR and PA content in cacao leaves, flowers, pod exocarp and seeds. In all tissues examined, all three genes were abundantly expressed and well correlated with PA accumulation levels, suggesting their active roles in PA synthesis. Overexpression of TcANR in an Arabidopsis ban mutant complemented the PA deficient phenotype in seeds and resulted in reduced anthocyanidin levels in hypocotyls. Overexpression of TcANS in tobacco resulted in increased content of both anthocyanidins and PAs in flower petals. Overexpression of TcANS in an Arabidopsis ldox mutant complemented its PA deficient phenotype in seeds. Recombinant TcLAR protein converted leucoanthocyanidin to catechin in vitro. Transgenic tobacco overexpressing TcLAR had decreased amounts of anthocyanidins and increased PAs. Overexpressing TcLAR in Arabidopsis ldox mutant also resulted in elevated synthesis of not only catechin but also epicatechin.

CONCLUSION

Our results confirm the in vivo function of cacao ANS and ANR predicted based on sequence homology to previously characterized enzymes from other species. In addition, our results provide a clear functional analysis of a LAR gene in vivo.

Tc-cAPX, a cytosolic ascorbate peroxidase of Theobroma cacao L. engaged in the interaction with Moniliophthora perniciosa, the causing agent of witches' broom disease

The level of hydrogen peroxide (H2O2) in plants signalizes the induction of several genes, including that of ascorbate peroxidase (APX-EC 1.11.1.11). APX isoenzymes play a central role in the elimination of intracellular H2O2 and contribute to plant responses to diverse stresses. During the infection process in Theobroma cacao by Moniliophthora perniciosa oxidative stress is generated and the APX action recruited from the plant. The present work aimed to characterize the T. cacao APX involved in the molecular interaction of T. cacao-M. perniciosa. The peroxidase activity was analyzed in protein extracts from cocoa plants infected by M. perniciosa and showed the induction of peroxidases like APX in resistant cocoa plants. The cytosolic protein of T. cacao (GenBank: ABR68691.2) was phylogenetically analyzed in relation to other peroxidases from the cocoa genome and eight genes encoding APX proteins with conserved domains were also analyzed. The cDNA from cytosolic APX was cloned in pET28a and the recombinant protein expressed and purified (rTc-cAPX). The secondary structure of the protein was analyzed by Circular Dichroism (CD) displaying high proportion of α-helices when folded. The enzymatic assay shows stable activity using ascorbate and guaiacol as an electron donor for H2O2 reduction. The pH 7.5 is the optimum for enzyme activity. Chromatographic analysis suggests that rTc-cAPX is a homodimer in solution. Results indicate that the rTc-cAPX is correctly folded, stable and biochemically active. The purified rTc-cAPX presented biotechnological potential and is adequate for future structural and functional studies.

The small RNA profile in latex from Hevea brasiliensis trees is affected by tapping panel dryness

Natural rubber is harvested by tapping Hevea brasiliensis (Willd. ex A. Juss.) Müll. Arg. Harvesting stress can lead to tapping panel dryness (TPD). MicroRNAs (miRNAs) are induced by abiotic stress and regulate gene expression by targeting the cleavage or translational inhibition of target messenger RNAs. This study set out to sequence miRNAs expressed in latex cells and to identify TPD-related putative targets. Deep sequencing of small RNAs was carried out on latex from trees affected by TPD using Solexa technology. The most abundant small RNA class size was 21 nucleotides for TPD trees compared with 24 nucleotides in healthy trees. By combining the LeARN pipeline, data from the Plant MicroRNA database and Hevea EST sequences, we identified 19 additional conserved and four putative species-specific miRNA families not found in previous studies on rubber. The relative transcript abundance of the Hbpre-MIR159b gene increased with TPD. This study revealed a small RNA-specific signature of TPD-affected trees. Both RNA degradation and a shift in miRNA biogenesis are suggested to explain the general decline in small RNAs and, particularly, in miRNAs.

Comparative transcriptome analysis of three oil palm fruit and seed tissues that differ in oil content and fatty acid composition

Oil palm (Elaeis guineensis) produces two oils of major economic importance, commonly referred to as palm oil and palm kernel oil, extracted from the mesocarp and the endosperm, respectively. While lauric acid predominates in endosperm oil, the major fatty acids (FAs) of mesocarp oil are palmitic and oleic acids. The oil palm embryo also stores oil, which contains a significant proportion of linoleic acid. In addition, the three tissues display high variation for oil content at maturity. To gain insight into the mechanisms that govern such differences in oil content and FA composition, tissue transcriptome and lipid composition were compared during development. The contribution of the cytosolic and plastidial glycolytic routes differed markedly between the mesocarp and seed tissues, but transcriptional patterns of genes involved in the conversion of sucrose to pyruvate were not related to variations for oil content. Accumulation of lauric acid relied on the dramatic up-regulation of a specialized acyl-acyl carrier protein thioesterase paralog and the concerted recruitment of specific isoforms of triacylglycerol assembly enzymes. Three paralogs of the WRINKLED1 (WRI1) transcription factor were identified, of which EgWRI1-1 and EgWRI1-2 were massively transcribed during oil deposition in the mesocarp and the endosperm, respectively. None of the three WRI1 paralogs were detected in the embryo. The transcription level of FA synthesis genes correlated with the amount of WRI1 transcripts and oil content. Changes in triacylglycerol content and FA composition of Nicotiana benthamiana leaves infiltrated with various combinations of WRI1 and FatB paralogs from oil palm validated functions inferred from transcriptome analysis.

Comparative transcriptome analysis of three oil palm fruit and seed tissues that differ in oil content and fatty acid composition

Oil palm (Elaeis guineensis) produces two oils of major economic importance, commonly referred to as palm oil and palm kernel oil, extracted from the mesocarp and the endosperm, respectively. While lauric acid predominates in endosperm oil, the major fatty acids (FAs) of mesocarp oil are palmitic and oleic acids. The oil palm embryo also stores oil, which contains a significant proportion of linoleic acid. In addition, the three tissues display high variation for oil content at maturity. To gain insight into the mechanisms that govern such differences in oil content and FA composition, tissue transcriptome and lipid composition were compared during development. The contribution of the cytosolic and plastidial glycolytic routes differed markedly between the mesocarp and seed tissues, but transcriptional patterns of genes involved in the conversion of sucrose to pyruvate were not related to variations for oil content. Accumulation of lauric acid relied on the dramatic up-regulation of a specialized acyl-acyl carrier protein thioesterase paralog and the concerted recruitment of specific isoforms of triacylglycerol assembly enzymes. Three paralogs of the WRINKLED1 (WRI1) transcription factor were identified, of which EgWRI1-1 and EgWRI1-2 were massively transcribed during oil deposition in the mesocarp and the endosperm, respectively. None of the three WRI1 paralogs were detected in the embryo. The transcription level of FA synthesis genes correlated with the amount of WRI1 transcripts and oil content. Changes in triacylglycerol content and FA composition of Nicotiana benthamiana leaves infiltrated with various combinations of WRI1 and FatB paralogs from oil palm validated functions inferred from transcriptome analysis.

The banana genome hub

Banana is one of the world’s favorite fruits and one of the most important crops for developing countries. The banana reference genome sequence (Musa acuminata) was recently released. Given the taxonomic position of Musa, the completed genomic sequence has particular comparative value to provide fresh insights about the evolution of the monocotyledons. The study of the banana genome has been enhanced by a number of tools and resources that allows harnessing its sequence. First, we set up essential tools such as a Community Annotation System, phylogenomics resources and metabolic pathways. Then, to support post-genomic efforts, we improved banana existing systems (e.g. web front end, query builder), we integrated available Musa data into generic systems (e.g. markers and genetic maps, synteny blocks), we have made interoperable with the banana hub, other existing systems containing Musa data (e.g. transcriptomics, rice reference genome, workflow manager) and finally, we generated new results from sequence analyses (e.g. SNP and polymorphism analysis). Several uses cases illustrate how the Banana Genome Hub can be used to study gene families. Overall, with this collaborative effort, we discuss the importance of the interoperability toward data integration between existing information systems.

Database URL: http://banana-genome.cirad.fr/

Identification of the Hevea brasiliensis AP2/ERF superfamily by RNA sequencing

Background

Rubber tree (Hevea brasiliensis) laticifers are the source of natural rubber. Rubber production depends on endogenous and exogenous ethylene (ethephon). AP2/ERF transcription factors, and especially Ethylene-Response Factors, play a crucial role in plant development and response to biotic and abiotic stresses. This study set out to sequence transcript expressed in various tissues using next-generation sequencing and to identify AP2/ERF superfamily in the rubber tree.

Results

The 454 sequencing technique was used to produce five tissue-type transcript libraries (leaf, bark, latex, embryogenic tissues and root). Reads from all libraries were pooled and reassembled to improve mRNA lengths and produce a global library. One hundred and seventy-three AP2/ERF contigs were identified by in silico analysis based on the amino acid sequence of the conserved AP2 domain from the global library. The 142 contigs with the full AP2 domain were classified into three main families (20 AP2 members, 115 ERF members divided into 11 groups, and 4 RAV members) and 3 soloist members. Fifty-nine AP2/ERF transcripts were found in latex. Alongside the microRNA172 already described in plants, eleven additional microRNAs were predicted to inhibit Hevea AP2/ERF transcripts.

Conclusions

Hevea has a similar number of AP2/ERF genes to that of other dicot species. We adapted the alignment and classification methods to data from next-generation sequencing techniques to provide reliable information. We observed several specific features for the ERF family. Three HbSoloist members form a group in Hevea. Several AP2/ERF genes highly expressed in latex suggest they have a specific function in Hevea. The analysis of AP2/ERF transcripts in Hevea presented here provides the basis for studying the molecular regulation of latex production in response to abiotic stresses and latex cell differentiation.

The Physalis peruviana leaf transcriptome: assembly, annotation and gene model prediction

Background

Physalis peruviana commonly known as Cape gooseberry is a member of the Solanaceae family that has an increasing popularity due to its nutritional and medicinal values. A broad range of genomic tools is available for other Solanaceae, including tomato and potato. However, limited genomic resources are currently available for Cape gooseberry.

Results

We report the generation of a total of 652,614 P. peruviana Expressed Sequence Tags (ESTs), using 454 GS FLX Titanium technology. ESTs, with an average length of 371 bp, were obtained from a normalized leaf cDNA library prepared using a Colombian commercial variety. De novo assembling was performed to generate a collection of 24,014 isotigs and 110,921 singletons, with an average length of 1,638 bp and 354 bp, respectively. Functional annotation was performed using NCBI’s BLAST tools and Blast2GO, which identified putative functions for 21,191 assembled sequences, including gene families involved in all the major biological processes and molecular functions as well as defense response and amino acid metabolism pathways. Gene model predictions in P. peruviana were obtained by using the genomes of Solanum lycopersicum (tomato) and Solanum tuberosum (potato). We predict 9,436 P. peruviana sequences with multiple-exon models and conserved intron positions with respect to the potato and tomato genomes. Additionally, to study species diversity we developed 5,971 SSR markers from assembled ESTs.

Conclusions

We present the first comprehensive analysis of the Physalis peruviana leaf transcriptome, which will provide valuable resources for development of genetic tools in the species. Assembled transcripts with gene models could serve as potential candidates for marker discovery with a variety of applications including: functional diversity, conservation and improvement to increase productivity and fruit quality. P. peruviana was estimated to be phylogenetically branched out before the divergence of five other Solanaceae family members, S. lycopersicum, S. tuberosum, Capsicum spp, S. melongena and Petunia spp.

cDNA-AFLP-based genetical genomics in cotton fibers

Genetical genomics, or genetic analysis applied to gene expression data, has not been widely used in plants. We used quantitative cDNA-AFLP to monitor the variation in the expression level of cotton fiber transcripts among a population of inter-specific Gossypium hirsutum × G. barbadense recombinant inbred lines (RILs). Two key fiber developmental stages, elongation (10 days post anthesis, dpa), and secondary cell wall thickening (22 dpa), were studied. Normalized intensity ratios of 3,263 and 1,201 transcript-derived fragments (TDFs) segregating over 88 RILs were analyzed for quantitative trait loci (QTL) mapping for the 10 and 22 dpa fibers, respectively. Two-thirds of all TDFs mapped between 1 and 6 eQTLs (LOD > 3.5). Chromosome 21 had a higher density of eQTLs than other chromosomes in both data sets and, within chromosomes, hotspots of presumably trans-acting eQTLs were identified. The eQTL hotspots were compared to the location of phenotypic QTLs for fiber characteristics among the RILs, and several cases of co-localization were detected. Quantitative RT-PCR for 15 sequenced TDFs showed that 3 TDFs had at least one eQTL at a similar location to those identified by cDNA-AFLP, while 3 other TDFs mapped an eQTL at a similar location but with opposite additive effect. In conclusion, cDNA-AFLP proved to be a cost-effective and highly transferable platform for genome-wide and population-wide gene expression profiling. Because TDFs are anonymous, further validation and interpretation (in silico analysis, qPCR gene profiling) of the eQTL and eQTL hotspots will be facilitated by the increasing availability of cDNA and genomic sequence resources in cotton.

Analysis of expressed sequence tags (ESTs) from cocoa (Theobroma cacao L) upon infection with Phytophthora megakarya

Phytophthora megakarya, the causative agent of cacao black pod disease in West African countries causes an extensive loss of yield. In this study we have analyzed 4 libraries of ESTs derived from Phytophthora megakarya infected cocoa leaf and pod tissues. Totally 6379 redundant sequences were retrieved from ESTtik database and EST processing was performed using seqclean tool. Clustering and assembling using CAP3 generated 3333 non-redundant (907 contigs and 2426 singletons) sequences. The primary sequence analysis of 3333 non-redundant sequences showed that the GC percentage was 42.7 and the sequence length ranged from 101 - 2576 nucleotides. Further, functional analysis (Blast, Interproscan, Gene ontology and KEGG search) were executed and 1230 orthologous genes were annotated. Totally 272 enzymes corresponding to 114 metabolic pathways were identified. Functional annotation revealed that most of the sequences are related to molecular function, stress response and biological processes. The annotated enzymes are aldehyde dehydrogenase (E.C: 1.2.1.3), catalase (E.C: 1.11.1.6), acetyl-CoA C-acetyltransferase (E.C: 2.3.1.9), threonine ammonia-lyase (E.C: 4.3.1.19), acetolactate synthase (E.C: 2.2.1.6), O-methyltransferase (E.C: 2.1.1.68) which play an important role in amino acid biosynthesis and phenyl propanoid biosynthesis. All this information was stored in MySQL database management system to be used in future for reconstruction of biotic stress response pathway in cocoa.