Comparative analysis of expressed genes from cacao meristems infected by Moniliophthora perniciosa

State of the Art of the Molecular Biology of the Interaction between Cocoa and Witches’ Broom Disease: A Systematic Review

Significant scientific advances to elucidate the Moniliophthora perniciosa pathosystem have been achieved in recent years, but the molecular biology of this pathogen-host interaction is still a field with many unanswered questions. In order to present insights at the molecular level, we present the first systematic review on the theme. All told, 1118 studies were extracted from public databases. Of these, 109 were eligible for the review, based on the inclusion and exclusion criteria. The results indicated that understanding the transition from the biotrophic-necrotrophic phase of the fungus is crucial for control of the disease. Proteins with strong biotechnological potential or that can be targets for pathosystem intervention were identified, but studies regarding possible applications are still limited. The studies identified revealed important genes in the M. perniciosa-host interaction and efficient molecular markers in the search for genetic variability and sources of resistance, with Theobroma cacao being the most common host. An arsenal of effectors already identified and not explored in the pathosystem were highlighted. This systematic review contributes to the understanding of the pathosystem at the molecular level, offering new insights and proposing different paths for the development of new strategies to control witches’ broom disease.

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.

The glutathione peroxidase family of Theobroma cacao: Involvement in the oxidative stress during witches' broom disease

The glutathione peroxidases (GPXs) are enzymes which are part of the cell antioxidant system inhibiting the ROS-induced damages of membranes and proteins. In cacao (Theobroma cacao L.) genome, five GPX genes were identified. Cysteine insertion codons (UGU) were found in TcPHGPX, TcGPX2, TcGPX4, TcGPX6 and tryptophan insertion codon (UGG) in TcGPX8. Multiple alignments revealed conserved domains between TcGPXs and other plants and human GPXs. Homology modeling was performed using the Populus trichocarpa GPX5 structure as template, and the molecular modeling showed that TcGPXs have affinity with selenometionine in their active site. In silico analysis of the TcGPXs promoter region revealed the presence of conserved cis-elements related to biotic stresses and hormone responsiveness. The expression analysis of TcGPXs in cacao plantlet meristems infected by M. perniciosa showed that TcGPXs are most expressed in susceptible variety than in resistant one, mainly in disease stages in which oxidative stress and programmed cell death occurred. This data, associated with phylogenetic and location analysis suggested that TcGPXs may play a role in protecting cells from oxidative stress as a try of disease progression reduction. To our knowledge, this is the first study of the overall GPX family from T. cacao.

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.

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.

The selenium-binding protein of Theobroma cacao: A thermostable protein involved in the witches' broom disease resistance

The selenium-binding proteins are known to be inducers of apoptosis in human and animals, and have been studied as target for the treatment of various types of cancer. In plants, SBP expression has been related to abiotic and biotic stress resistance. The SBP from Theobroma cacao (TcSBP) was first identified from a cocoa-Moniliophthora perniciosa cDNA library. The present study provides details on the TcSBP gene and protein structure. Multiple alignments revealed conserved domains between SBP from plants, human and archea. Homology modeling and molecular docking were performed and showed that the TcSBP has affinity to selenite in the active CSSC site. This result was confirmed by circular dichroism of the recombinant TcSBP, which also presented thermostable behavior. RT-qPCR analysis showed that TcSBP was differentially expressed in resistant vs susceptible cacao varieties inoculated by M. perniciosa and its expression was probably due to hormone induction via cis-regulating elements present in its promotor. The presence of the CSSC domain suggested that TcSBP acted by altering oxidation/reduction of proteins during H₂O₂ production and programmed cell death in the final stages of the witches' broom disease. To our knowledge, this is the first in silico and in vitro analysis of the SBP from cacao.

Developing Blight-Tolerant American Chestnut Trees

An invasive fungal pathogen has reduced the American chestnut (Castanea dentata), once a keystone tree species within its natural range in the eastern United States and Canada, to functional extinction. To help restore this important canopy tree, blight-tolerant American chestnut trees have been developed using an oxalate oxidase-encoding gene from wheat. This enzyme breaks down oxalate, which is produced by the pathogen and forms killing cankers. Expressing oxalate oxidase results in blight tolerance, where the tree and the fungus can coexist, which is a more evolutionarily stable relationship than direct pathogen resistance. Genetic engineering (GE) typically makes a very small change in the tree's genome, potentially avoiding incompatible gene interactions that have been detected in some chestnut hybrids. The GE American chestnut also retains all the wild American chestnut's alleles for habitat adaptation, which are important for a forest ecosystem restoration program.

Sensitivity of Yeast Mutants Deficient in Mitochondrial or Vacuolar ABC Transporters to Pathogenesis-Related Protein TcPR-10 of Theobroma cacao

Pathogenesis-related proteins (PRs) are induced in plants after infection by pathogens and/or abiotic stress. Among these proteins, the family 10 (PR-10) influences the biosynthesis of secondary metabolites and shows antimicrobial ribonuclease activity. TcPR-10p (Pathogenesis-related Protein 10 of Theobroma cacao) was isolated from resistant and susceptible Moniliophthora perniciosa cacao cultivars. Cell survival with Saccharomyces cerevisiae mutant lines deficient in ATP-binding cassette (ABC) transporter proteins indicated the influence on resistance to TcPR-10p. Proteins of the ABC transport type are considered important in the process of resistance to antimicrobials and toxins. Thus, the objective of this work was to observe the sensitivity of ABC transporter yeast mutants in the presence of the TcPR-10p. Chronic exposure of S. cerevisiae mitochondrial (BYatm1Δ and BYmdl1Δ) and vacuole (BYnft1Δ, BYvmr1Δ, BYybt1Δ, BYycf1Δ and BYbpt1Δ) ABC transporter mutants to TcPR-10p (3 μg/mL, 0, 6, 12 and 24 h) was performed. Two TcPR-10p sensitive strains (BYmdl1Δ and BYnft1Δ) were submitted to a fluorescence test with the fluorogenic dihydroethidium (DHE), to visualize the presence of oxidative stress in the cells. Oxidative stress-increased sensitivity was confirmed by flow cytometry indicating induced cell death either via apoptosis or necrosis. This yeast data combined with previous data of literature (of M. perniciosa sensitivity to TcPR-10p) show that increased sensitivity to TcPR-10p in these mutants could be due to the TcPR10p-generated higher levels of intracellular reactive oxygen species (ROS), leading to increased cell death either via necrosis or apoptosis.

Genome-wide identification and characterization of cacao WRKY transcription factors and analysis of their expression in response to witches' broom disease

Transcriptional regulation, led by transcription factors (TFs) such as those of the WRKY family, is a mechanism used by the organism to enhance or repress gene expression in response to stimuli. Here, we report on the genome-wide analysis of the Theobroma cacao WRKY TF family and also investigate the expression of WRKY genes in cacao infected by the fungus Moniliophthora perniciosa. In the cacao genome, 61 non-redundant WRKY sequences were found and classified in three groups (I to III) according to the WRKY and zinc-finger motif types. The 61 putative WRKY sequences were distributed on the 10 cacao chromosomes and 24 of them came from duplication events. The sequences were phylogenetically organized according to the general WRKY groups. The phylogenetic analysis revealed that subgroups IIa and IIb are sister groups and share a common ancestor, as well as subgroups IId and IIe. The most divergent groups according to the plant origin were IIc and III. According to the phylogenetic analysis, 7 TcWRKY genes were selected and analyzed by RT-qPCR in susceptible and resistant cacao plants infected (or not) with M. perniciosa. Some TcWRKY genes presented interesting responses to M. perniciosa such as Tc01_p014750/Tc06_p013130/AtWRKY28, Tc09_p001530/Tc06_p004420/AtWRKY40, Tc04_p016130/AtWRKY54 and Tc10_p016570/ AtWRKY70. Our results can help to select appropriate candidate genes for further characterization in cacao or in other Theobroma species.

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.

Genome size, cytogenetic data and transferability of EST-SSRs markers in wild and cultivated species of the genus Theobroma L. (Byttnerioideae, Malvaceae)

The genus Theobroma comprises several trees species native to the Amazon. Theobroma cacao L. plays a key economic role mainly in the chocolate industry. Both cultivated and wild forms are described within the genus. Variations in genome size and chromosome number have been used for prediction purposes including the frequency of interspecific hybridization or inference about evolutionary relationships. In this study, the nuclear DNA content, karyotype and genetic diversity using functional microsatellites (EST-SSR) of seven Theobroma species were characterized. The nuclear content of DNA for all analyzed Theobroma species was 1C = ~ 0.46 pg. These species presented 2n = 20 with small chromosomes and only one pair of terminal heterochromatic bands positively stained (CMA⁺/DAPI⁻ bands). The small size of Theobroma ssp. genomes was equivalent to other Byttnerioideae species, suggesting that the basal lineage of Malvaceae have smaller genomes and that there was an expansion of 2C values in the more specialized family clades. A set of 20 EST-SSR primers were characterized for related species of Theobroma, in which 12 loci were polymorphic. The polymorphism information content (PIC) ranged from 0.23 to 0.65, indicating a high level of information per locus. Combined results of flow cytometry, cytogenetic data and EST-SSRs markers will contribute to better describe the species and infer about the evolutionary relationships among Theobroma species. In addition, the importance of a core collection for conservation purposes is highlighted.

Association mapping of seed and disease resistance traits in Theobroma cacao L

Microsatellite and single nucleotide polymorphism markers that could be used in marker assisted breeding of cacao were identified for number of filled seeds, black pod resistance and witches' broom disease resistance. An association mapping approach was employed to identify markers for seed number and resistance to black pod and witches' broom disease (WBD) in cacao (Theobroma cacao L.). Ninety-five microsatellites (SSRs) and 775 single nucleotide polymorphisms (SNPs) were assessed on 483 unique trees in the International Cocoa Genebank Trinidad (ICGT). Linkage disequilibrium (LD) and association mapping studies were conducted to identify markers to tag the phenotypic traits. Decay of LD occurred over an average 9.3 cM for chromosomes 1-9 and 2.5 cM for chromosome 10. Marker/trait associations were generally identified based on general linear models (GLMs) that incorporated principal components from molecular information on relatedness factor. Seven markers (mTcCIR 8, 66, 126, 212; TcSNP368, 697, 1370) on chromosomes 1 and 9 were identified for number of filled seeds (NSEED). A single marker was found for black pod resistance (mTcCIR280) on chromosome 3, whereas six markers on chromosomes 4, 5, 6, 8, and 10 were detected for WBD (mTcCIR91, 183; TcSNP375, 720, 1230 and 1374). It is expected that this association mapping study in cacao would contribute to the knowledge of the genetic determinism of cocoa traits and that the markers identified herein would prove useful in marker assisted breeding of cacao.

Protein profile and protein interaction network of Moniliophthora perniciosa basidiospores

BACKGROUND

Witches' broom, a disease caused by the basidiomycete Moniliophthora perniciosa, is considered to be the most important disease of the cocoa crop in Bahia, an area in the Brazilian Amazon, and also in the other countries where it is found. M. perniciosa germ tubes may penetrate into the host through intact or natural openings in the cuticle surface, in epidermis cell junctions, at the base of trichomes, or through the stomata. Despite its relevance to the fungal life cycle, basidiospore biology has not been extensively investigated. In this study, our goal was to optimize techniques for producing basidiospores for protein extraction, and to produce the first proteomics analysis map of ungerminated basidiospores. We then presented a protein interaction network by using Ustilago maydis as a model.

RESULTS

The average pileus area ranged from 17.35 to 211.24 mm(2). The minimum and maximum productivity were 23,200 and 6,666,667 basidiospores per basidiome, respectively. The protein yield in micrograms per million basidiospores were approximately 0.161; 2.307, and 3.582 for germination times of 0, 2, and 4 h after germination, respectively. A total of 178 proteins were identified through mass spectrometry. These proteins were classified according to their molecular function and their involvement in biological processes such as cellular energy production, oxidative metabolism, stress, protein synthesis, and protein folding. Furthermore, to better understand the expression pattern, signaling, and interaction events of spore proteins, we presented an interaction network using orthologous proteins from Ustilago maydis as a model. Most of the orthologous proteins that were identified in this study were not clustered in the network, but several of them play a very important role in hypha development and branching.

CONCLUSIONS

The quantities of basidiospores 7 × 10(9); 5.2 × 10(8), and 6.7 × 10(8) were sufficient to obtain enough protein mass for the three 2D-PAGE replicates, for the 0, 2, and 4 h-treatments, respectively. The protein extraction method that is based on sedimentation, followed by sonication with SDS-dense buffer, and phenolic extraction, which was utilized in this study, was effective, presenting a satisfactory resolution and reproducibility for M. perniciosa basidiospores. This report constitutes the first comprehensive study of protein expression during the ungerminated stage of the M. perniciosa basidiospore. Identification of the spots observed in the reference gel enabled us to know the main molecular interactions involved in the initial metabolic processes of fungal development.

Changes in cellular infrastructure after induced endoplasmic reticulum stress in Moniliophthora perniciosa

Moniliophthora perniciosa is a basidiomycete fungus that causes witches' broom disease in Theobroma cacao We analyzed the morphology and survival of fungal hyphae and endoplasmic reticulum (ER) remodeling in either glucose- or glycerol-grown M. perniciosa after treatment with ER stress-inducing chemicals dithiothreitol (DTT) or tunicamycin (TM). Changes in intracellular redox potential can cause endoplasmic reticulum (ER) stress due to diminished efficiency in protein folding that could in turn reduce cell survival. Such stress diminishes protein-folding efficiency that could in turn reduce cell survival. Light microscopy revealed morphological changes in hyphae after TM but not after DTT treatment, regardless of the media carbon source. Decrease in fungal survival, after both TM and DTT treatments, was dose-dependent and glycerol-grown cells showed a higher resistance to both chemicals compared to glucose-grown cells. Electron microscopy showed TM and DDT-induced ER stress in M. perniciosa as evidenced by structural alterations of the organelle. The volume of ER structures increased as a typical consequence of unfolded protein stress, and the number of autophagosomes was higher. In glycerol-grown fungus DTT treatment slightly induced expression of molecular chaperone BiP. The TM exposure-induced expression of gene MpIRE1, involved in signaling of the unfolded protein response, was higher in glycerol than glucose-grown cells. Such difference was not observable with expression of gene MpATG8, encoding a key protein in autosome formation, that was induced 1.4-fold and 1.2-fold in glucose or glycerol-grown cells, respectively. DHE-based fluorescence assay showed M. perniciosa oxidative stress induced by H₂O₂, and treated cells had a higher level of oxidative stress compared to control. A comprehensive study of remodeling of ER is important in understanding M. perniciosa fungus resistance to oxidative stress and its ability to implement a successful infection in T. cacao.

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.

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.

Differential expression of jasmonate biosynthesis genes in cacao genotypes contrasting for resistance against Moniliophthora perniciosa

The resistance mechanism of cacao against M. perniciosa is likely to be mediated by JA/ET-signaling pathways due to the preferential TcAOS and TcSAM induction in a resistant genotype. The basidiomycete Moniliophthora perniciosa causes a serious disease in cacao (Theobroma cacao L.), and the use of resistant varieties is the only sustainable long-term solution. Cacao resistance against M. perniciosa is characterized by pathogen growth inhibition with reduced colonization and an attenuation of disease symptoms, suggesting a regulation by jasmonate (JA)/ethylene (ET) signaling pathways. The hypothesis that genes involved in JA biosynthesis would be active in the interaction of T. cacao and M. perniciosa was tested here. The cacao JA-related genes were evaluated for their relative quantitative expression in susceptible and resistant genotypes upon the exogenous application of ET, methyl-jasmonate (MJ), and salicylic acid (SA), or after M. perniciosa inoculation. MJ treatment triggered changes in the expression of genes involved in JA biosynthesis, indicating that the mechanism of positive regulation by exogenous MJ application occurs in cacao. However, a higher induction of these genes was observed in the susceptible genotype. Further, a contrast in JA-related transcriptional expression was detected between susceptible and resistant plants under M. perniciosa infection, with the induction of the allene oxide synthase gene (TcAOS), which encodes a key enzyme in the JA biosynthesis pathway in the resistant genotype. Altogether, this work provides additional evidences that the JA-dependent signaling pathway is modulating the defense response against M. perniciosa in a cacao-resistant genotype.

Altered Phenotypes in Saccharomyces cerevisiae by Heterologous Expression of Basidiomycete Moniliophthora perniciosa SOD2 Gene

Heterologous expression of a putative manganese superoxide dismutase gene (SOD2) of the basidiomycete Moniliophthora perniciosa complemented the phenotypes of a Saccharomyces cerevisiae sod2Δ mutant. Sequence analysis of the cloned M. perniciosa cDNA revealed an open reading frame (ORF) coding for a 176 amino acid polypeptide with the typical metal-binding motifs of a SOD2 gene, named MpSOD2. Phylogenetic comparison with known manganese superoxide dismutases (MnSODs) located the protein of M. perniciosa (MpSod2p) in a clade with the basidiomycete fungi Coprinopsis cinerea and Laccaria bicolor. Haploid wild-type yeast transformants containing a single copy of MpSOD2 showed increased resistance phenotypes against oxidative stress-inducing hydrogen peroxide and paraquat, but had unaltered phenotype against ultraviolet-C (UVC) radiation. The same transformants exhibited high sensitivity against treatment with the pro-mutagen diethylnitrosamine (DEN) that requires oxidation to become an active mutagen/carcinogen. Absence of MpSOD2 in the yeast sod2Δ mutant led to DEN hyper-resistance while introduction of a single copy of this gene restored the yeast wild-type phenotype. The haploid yeast wild-type transformant containing two SOD2 gene copies, one from M. perniciosa and one from its own, exhibited DEN super-sensitivity. This transformant also showed enhanced growth at 37 °C on the non-fermentable carbon source lactate, indicating functional expression of MpSod2p. The pro-mutagen dihydroethidium (DHE)-based fluorescence assay monitored basal level of yeast cell oxidative stress. Compared to the wild type, the yeast sod2Δ mutant had a much higher level of intrinsic oxidative stress, which was reduced to wild type (WT) level by introduction of one copy of the MpSOD2 gene. Taken together our data indicates functional expression of MpSod2 protein in the yeast S. cerevisiae.

The activity of TcCYS4 modified by variations in pH and temperature can affect symptoms of witches' broom disease of cocoa, caused by the fungus Moniliophthora perniciosa

The phytocystatins regulate various physiological processes in plants, including responses to biotic and abiotic stresses, mainly because they act as inhibitors of cysteine proteases. In this study, we have analyzed four cystatins from Theobroma cacao L. previously identified in ESTs libraries of the interaction with the fungus Moniliophthora perniciosa and named TcCYS1, TcCYS2, TcCYS3 and TcCYS4. The recombinant cystatins were purified and subjected to the heat treatment, at different temperatures, and their thermostabilities were monitored using their ability to inhibit papain protease. TcCYS1 was sensitive to temperatures above 50°C, while TcCYS2, TcCYS3, and TcCYS4 were thermostable. TcCYS4 presented a decrease of inhibitory activity when it was treated at temperatures between 60 and 70°C, with the greater decrease occurring at 65°C. Analyses by native gel electrophoresis and size-exclusion chromatography showed that TcCYS4 forms oligomers at temperatures between 60 and 70°C, condition where reduction of inhibitory activity was observed. TcCYS4 oligomers remain stable for up to 20 days after heat treatment and are undone after treatment at 80°C. TcCYS4 presented approximately 90% of inhibitory activity at pH values between 5 and 9. This protein treated at temperatures above 45°C and pH 5 presented reduced inhibitory activity against papain, suggesting that the pH 5 enhances the formation of TcCYS4 oligomers. A variation in the titratable acidity was observed in tissues of T. cacao during the symptoms of witches’ broom disease. Our findings suggest that the oligomerization of TcCYS4, favored by variations in pH, is an endergonic process. We speculate that this process can be involved in the development of the symptoms of witches’ broom disease in cocoa.

High-resolution transcript profiling of the atypical biotrophic interaction between Theobroma cacao and the fungal pathogen Moniliophthora perniciosa

Witches' broom disease (WBD), caused by the hemibiotrophic fungus Moniliophthora perniciosa, is one of the most devastating diseases of Theobroma cacao, the chocolate tree. In contrast to other hemibiotrophic interactions, the WBD biotrophic stage lasts for months and is responsible for the most distinctive symptoms of the disease, which comprise drastic morphological changes in the infected shoots. Here, we used the dual RNA-seq approach to simultaneously assess the transcriptomes of cacao and M. perniciosa during their peculiar biotrophic interaction. Infection with M. perniciosa triggers massive metabolic reprogramming in the diseased tissues. Although apparently vigorous, the infected shoots are energetically expensive structures characterized by the induction of ineffective defense responses and by a clear carbon deprivation signature. Remarkably, the infection culminates in the establishment of a senescence process in the host, which signals the end of the WBD biotrophic stage. We analyzed the pathogen's transcriptome in unprecedented detail and thereby characterized the fungal nutritional and infection strategies during WBD and identified putative virulence effectors. Interestingly, M. perniciosa biotrophic mycelia develop as long-term parasites that orchestrate changes in plant metabolism to increase the availability of soluble nutrients before plant death. Collectively, our results provide unique insight into an intriguing tropical disease and advance our understanding of the development of (hemi)biotrophic plant-pathogen interactions.

The pathogenesis-related protein PR-4b from Theobroma cacao presents RNase activity, Ca(2+) and Mg(2+) dependent-DNase activity and antifungal action on Moniliophthora perniciosa

BACKGROUND

The production and accumulation of pathogenesis-related proteins (PR proteins) in plants in response to biotic or abiotic stresses is well known and is considered as a crucial mechanism for plant defense. A pathogenesis-related protein 4 cDNA was identified from a cacao-Moniliophthora perniciosa interaction cDNA library and named TcPR-4b.

RESULTS

TcPR-4b presents a Barwin domain with six conserved cysteine residues, but lacks the chitin-binding site. Molecular modeling of TcPR-4b confirmed the importance of the cysteine residues to maintain the protein structure, and of several conserved amino acids for the catalytic activity. In the cacao genome, TcPR-4b belonged to a small multigene family organized mainly on chromosome 5. TcPR-4b RT-qPCR analysis in resistant and susceptible cacao plants infected by M. perniciosa showed an increase of expression at 48 hours after infection (hai) in both cacao genotypes. After the initial stage (24-72 hai), the TcPR-4b expression was observed at all times in the resistant genotypes, while in the susceptible one the expression was concentrated at the final stages of infection (45-90 days after infection). The recombinant TcPR-4b protein showed RNase, and bivalent ions dependent-DNase activity, but no chitinase activity. Moreover, TcPR-4b presented antifungal action against M. perniciosa, and the reduction of M. perniciosa survival was related to ROS production in fungal hyphae.

CONCLUSION

To our knowledge, this is the first report of a PR-4 showing simultaneously RNase, DNase and antifungal properties, but no chitinase activity. Moreover, we showed that the antifungal activity of TcPR-4b is directly related to RNase function. In cacao, TcPR-4b nuclease activities may be related to the establishment and maintenance of resistance, and to the PCD mechanism, in resistant and susceptible cacao genotypes, respectively.

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.

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.

Evaluation of the allergenicity potential of TcPR-10 protein from Theobroma cacao

BACKGROUND

The pathogenesis related protein PR10 (TcPR-10), obtained from the Theobroma cacao-Moniliophthora perniciosa interaction library, presents antifungal activity against M. perniciosa and acts in vitro as a ribonuclease. However, despite its biotechnological potential, the TcPR-10 has the P-loop motif similar to those of some allergenic proteins such as Bet v 1 (Betula verrucosa) and Pru av 1 (Prunus avium). The insertion of mutations in this motif can produce proteins with reduced allergenic power. The objective of the present work was to evaluate the allergenic potential of the wild type and mutant recombinant TcPR-10 using bioinformatics tools and immunological assays.

METHODOLOGY/PRINCIPAL FINDINGS

Mutant substitutions (T10P, I30V, H45S) were inserted in the TcPR-10 gene by site-directed mutagenesis, cloned into pET28a and expressed in Escherichia coli BL21(DE3) cells. Changes in molecular surface caused by the mutant substitutions was evaluated by comparative protein modeling using the three-dimensional structure of the major cherry allergen, Pru av 1 as a template. The immunological assays were carried out in 8-12 week old female BALB/c mice. The mice were sensitized with the proteins (wild type and mutants) via subcutaneous and challenged intranasal for induction of allergic airway inflammation.

CONCLUSIONS/SIGNIFICANCE

We showed that the wild TcPR-10 protein has allergenic potential, whereas the insertion of mutations produced proteins with reduced capacity of IgE production and cellular infiltration in the lungs. On the other hand, in vitro assays show that the TcPR-10 mutants still present antifungal and ribonuclease activity against M. perniciosa RNA. In conclusion, the mutant proteins present less allergenic potential than the wild TcPR-10, without the loss of interesting biotechnological properties.

Transcriptomics and systems biology analysis in identification of specific pathways involved in cacao resistance and susceptibility to witches' broom disease

This study reports on expression analysis associated with molecular systems biology of cacao-Moniliophthora perniciosa interaction. Gene expression data were obtained for two cacao genotypes (TSH1188, resistant; Catongo, susceptible) challenged or not with the fungus M. perniciosa and collected at three time points through disease. Using expression analysis, we identified 154 and 227 genes that are differentially expressed in TSH1188 and Catongo, respectively. The expression of some of these genes was confirmed by RT-qPCR. Physical protein-protein interaction (PPPI) networks of Arabidopsis thaliana orthologous proteins corresponding to resistant and susceptible interactions were obtained followed by cluster and gene ontology analyses. The integrated analysis of gene expression and systems biology allowed designing a general scheme of major mechanisms associated with witches' broom disease resistance/susceptibility. In this sense, the TSH1188 cultivar shows strong production of ROS and elicitors at the beginning of the interaction with M. perniciosa followed by resistance signal propagation and ROS detoxification. On the other hand, the Catongo genotype displays defense mechanisms that include the synthesis of some defense molecules but without success in regards to elimination of the fungus. This phase is followed by the activation of protein metabolism which is achieved with the production of proteasome associated with autophagy as a precursor mechanism of PCD. This work also identifies candidate genes for further functional studies and for genetic mapping and marker assisted selection.

The mitochondrial genome of Moniliophthora roreri, the frosty pod rot pathogen of cacao

In this study, we report the sequence of the mitochondrial (mt) genome of the Basidiomycete fungus Moniliophthora roreri, which is the etiologic agent of frosty pod rot of cacao (Theobroma cacao L.). We also compare it to the mtDNA from the closely-related species Moniliophthora perniciosa, which causes witches' broom disease of cacao. The 94 Kb mtDNA genome of M. roreri has a circular topology and codes for the typical 14 mt genes involved in oxidative phosphorylation. It also codes for both rRNA genes, a ribosomal protein subunit, 13 intronic open reading frames (ORFs), and a full complement of 27 tRNA genes. The conserved genes of M. roreri mtDNA are completely syntenic with homologous genes of the 109 Kb mtDNA of M. perniciosa. As in M. perniciosa, M. roreri mtDNA contains a high number of hypothetical ORFs (28), a remarkable feature that make Moniliophthoras the largest reservoir of hypothetical ORFs among sequenced fungal mtDNA. Additionally, the mt genome of M. roreri has three free invertron-like linear mt plasmids, one of which is very similar to that previously described as integrated into the main M. perniciosa mtDNA molecule. Moniliophthora roreri mtDNA also has a region of suspected plasmid origin containing 15 hypothetical ORFs distributed in both strands. One of these ORFs is similar to an ORF in the mtDNA gene encoding DNA polymerase in Pleurotus ostreatus. The comparison to M. perniciosa showed that the 15 Kb difference in mtDNA sizes is mainly attributed to a lower abundance of repetitive regions in M. roreri (5.8 Kb vs 20.7 Kb). The most notable differences between M. roreri and M. perniciosa mtDNA are attributed to repeats and regions of plasmid origin. These elements might have contributed to the rapid evolution of mtDNA. Since M. roreri is the second species of the genus Moniliophthora whose mtDNA genome has been sequenced, the data presented here contribute valuable information for understanding the evolution of fungal mt genomes among closely-related species.

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.

Hydrogen peroxide formation in cacao tissues infected by the hemibiotrophic fungus Moniliophthora perniciosa

In plant-pathogen interaction, the hydrogen peroxide (H₂O₂) may play a dual role: its accumulation inhibits the growth of biotrophic pathogens, while it could help the infection/colonization process of plant by necrotrophic pathogens. One of the possible pathways of H₂O production involves oxalic acid (Oxa) degradation by apoplastic oxalate oxidase. Here, we analyzed the production of H₂O₂, the presence of calcium oxalate (CaOx) crystals and the content of Oxa and ascorbic acid (Asa)--the main precursor of Oxa in plants--in susceptible and resistant cacao (Theobroma cacao L.) infected by the hemibiotrophic fungus Moniliophthora perniciosa. We also quantified the transcript level of ascorbate peroxidase (Apx), germin-like oxalate oxidase (Glp) and dehydroascorbate reductase (Dhar) by RT-qPCR. We report that the CaOx crystal amount and the H₂O₂ levels in the two varieties present distinct temporal and genotype-dependent patterns. Susceptible variety accumulated more CaOx crystals than the resistant one, and the dissolution of these crystals occurred in the early infection steps and in the final stage of the disease in the resistant and the susceptible variety, respectively. High expression of the Glp and accumulation of Oxa were observed in the resistant variety. The content of Asa increased in the inoculated susceptible variety, but remained constant in the resistant one. The susceptible variety presented reduced Dhar expression. The role of H₂O₂ and its formation from Oxa via Apx and Glp in resistant and susceptible variety infected by M. perniciosa were discussed.

Theobroma cacao cystatins impair Moniliophthora perniciosa mycelial growth and are involved in postponing cell death symptoms

Three cystatin open reading frames named TcCys1, TcCys2 and TcCys3 were identified in cDNA libraries from compatible interactions between Theobroma cacao (cacao) and Moniliophthora perniciosa. In addition, an ORF named TcCys4 was identified in the cDNA library of the incompatible interaction. The cDNAs encoded conceptual proteins with 209, 127, 124, and 205 amino acid residues, with a deduced molecular weight of 24.3, 14.1, 14.3 and 22.8 kDa, respectively. His-tagged recombinant proteins were purified from Escherichia coli expression, and showed inhibitory activities against M. perniciosa. The four recombinant cystatins exhibited K(i) values against papain in the range of 152-221 nM. Recombinant TcCYS3 and TcCYS4 immobilized in CNBr-Sepharose were efficient to capture M. perniciosa proteases from culture media. Polyclonal antibodies raised against the recombinant TcCYS4 detected that the endogenous protein was more abundant in young cacao tissues, when compared with mature tissues. A ~85 kDa cacao multicystatin induced by M. perniciosa inoculation, MpNEP (necrosis and ethylene-inducing protein) and M. perniciosa culture supernatant infiltration were detected by anti-TcCYS4 antibodies in cacao young tissues. A direct role of the cacao cystatins in the defense against this phytopathogen was proposed, as well as its involvement in the development of symptoms of programmed cell death.

Expression analysis of transcription factors from the interaction between cacao and Moniliophthora perniciosa (Tricholomataceae)

Cacao (Theobroma cacao) is one of the most important tropical crops; however, production is threatened by numerous pathogens, including the hemibiotrophic fungus Moniliophthora perniciosa, which causes witches' broom disease. To understand the mechanisms that lead to the development of this disease in cacao, we focused our attention on cacao transcription factors (TFs), which act as master regulators of cellular processes and are important for the fine-tuning of plant defense responses. We developed a macroarray with 88 TF cDNA from previously obtained cacao-M. perniciosa interaction libraries. Seventy-two TFs were found differentially expressed between the susceptible (Catongo) and resistant (TSH1188) genotypes and/or during the disease time course--from 24 h to 30 days after infection. Most of the differentially expressed TFs belonged to the bZIP, MYB and WRKY families and presented opposite expression patterns in susceptible and resistant cacao-M. perniciosa interactions (i.e., up-regulated in Catongo and down-regulated in TSH1188). The results of the macroarray were confirmed for bZIP and WRKY TFs by real-time PCR. These differentially expressed TFs are good candidates for subsequent functional analysis as well as for plant engineering. Some of these TFs could also be localized on the cacao reference map related to witches' broom resistance, facilitating the breeding and selection of resistant cacao trees.

Comparative gene expression analysis of susceptible and resistant near-isogenic lines in common wheat infected by Puccinia triticina

Gene expression after leaf rust infection was compared in near-isogenic wheat lines differing in the Lr10 leaf rust resistance gene. RNA from susceptible and resistant plants was used for cDNA library construction. In total, 55 008 ESTs were sequenced from the two libraries, then combined and assembled into 14 268 unigenes for further analysis. Of these ESTs, 89% encoded proteins similar to (E value of < or =10(-5)) characterized or annotated proteins from the NCBI non-redundant database representing diverse molecular functions, cellular localization and biological processes based on gene ontology classification. Further, the unigenes were classified into susceptible and resistant classes based on the EST members assembled from the respective libraries. Several genes from the resistant sample (14-3-3 protein, wali5 protein, actin-depolymerization factor and ADP-ribosylation factor) and the susceptible sample (brown plant hopper resistance protein, caffeic acid O-methyltransferase, pathogenesis-related protein and senescence-associated protein) were selected and their differential expression in the resistant and susceptible samples collected at different time points after leaf rust infection was confirmed by RT-PCR analysis. The molecular pathogenicity of leaf rust in wheat was studied and the EST data generated made a foundation for future studies.

Early development of Moniliophthora perniciosa basidiomata and developmentally regulated genes

BACKGROUND

The hemibiotrophic fungus Moniliophthora perniciosa is the causal agent of Witches' broom, a disease of Theobroma cacao. The pathogen life cycle ends with the production of basidiocarps in dead tissues of the infected host. This structure generates millions of basidiospores that reinfect young tissues of the same or other plants. A deeper understanding of the mechanisms underlying the sexual phase of this fungus may help develop chemical, biological or genetic strategies to control the disease.

RESULTS

Mycelium was morphologically analyzed prior to emergence of basidiomata by stereomicroscopy, light microscopy and scanning electron microscopy. The morphological changes in the mycelium before fructification show a pattern similar to other members of the order Agaricales. Changes and appearance of hyphae forming a surface layer by fusion were correlated with primordia emergence. The stages of hyphal nodules, aggregation, initial primordium and differentiated primordium were detected. The morphological analysis also allowed conclusions on morphogenetic aspects. To analyze the genes involved in basidiomata development, the expression of some selected EST genes from a non-normalized cDNA library, representative of the fruiting stage of M. perniciosa, was evaluated. A macroarray analysis was performed with 192 selected clones and hybridized with two distinct RNA pools extracted from mycelium in different phases of basidiomata formation. This analysis showed two groups of up and down-regulated genes in primordial phases of mycelia. Hydrophobin coding, glucose transporter, Rho-GEF, Rheb, extensin precursor and cytochrome p450 monooxygenase genes were grouped among the up-regulated. In the down-regulated group relevant genes clustered coding calmodulin, lanosterol 14 alpha demethylase and PIM1. In addition, 12 genes with more detailed expression profiles were analyzed by RT-qPCR. One aegerolysin gene had a peak of expression in mycelium with primordia and a second in basidiomata, confirming their distinctiveness. The number of transcripts of the gene for plerototolysin B increased in reddish-pink mycelium and indicated an activation of the initial basidiomata production even at this culturing stage. Expression of the glucose transporter gene increased in mycelium after the stress, coinciding with a decrease of adenylate cyclase gene transcription. This indicated that nutrient uptake can be an important signal to trigger fruiting in this fungus.

CONCLUSION

The identification of genes with increased expression in this phase of the life cycle of M. perniciosa opens up new possibilities of controlling fungus spread as well as of genetic studies of biological processes that lead to basidiomycete fruiting. This is the first comparative morphologic study of the early development both in vivo and in vitro of M. perniciosa basidiomata and the first description of genes expressed at this stage of the fungal life cycle.

dsRNA-induced gene silencing in Moniliophthora perniciosa, the causal agent of witches' broom disease of cacao

The genome sequence of the hemibiotrophic fungus Moniliophthora perniciosa revealed genes possibly participating in the RNAi machinery. Therefore, studies were performed in order to investigate the efficiency of gene silencing by dsRNA. We showed that the reporter gfp gene stably introduced into the fungus genome can be silenced by transfection of in vitro synthesized gfpdsRNA. In addition, successful dsRNA-induced silencing of endogenous genes coding for hydrophobins and a peroxiredoxin were also achieved. All genes showed a silencing efficiency ranging from 18% to 98% when compared to controls even 28d after dsRNA treatment, suggesting systemic silencing. Reduction of GFP fluorescence, peroxidase activity levels and survival responses to H(2)O(2) were consistent with the reduction of GFP and peroxidase mRNA levels, respectively. dsRNA transformation of M. perniciosa is shown here to efficiently promote genetic knockdown and can thus be used to assess gene function in this pathogen.

High-affinity copper transport and Snq2 export permease of saccharomyces cerevisiae modulate cytotoxicity of PR-10 from Theobroma cacao

A pathogenesis-related (PR) protein from Theobroma cacao (TcPR-10) was identified from a cacao-Moniliophthora perniciosa interaction cDNA library. Nucleotide and amino acid sequences showed homology with other PR-10 proteins having P loop motif and Betv1 domain. Recombinant TcPR-10 showed in vitro and in vivo ribonuclease activity, and antifungal activity against the basidiomycete cacao pathogen M. perniciosa and the yeast Saccharomyces cerevisiae. Fluorescein isothiocyanate-labeled TcPR-10 was internalized by M. perniciosa hyphae and S. cerevisiae cells and inhibited growth of both fungi. Energy and temperature-dependent internalization of the TcPR-10 suggested an active importation into the fungal cells. Chronical exposure to TcPR-10 of 29 yeast mutants with single gene defects in DNA repair, general membrane transport, metal transport, and antioxidant defenses was tested. Two yeast mutants were hyperresistant compared with their respective isogenic wild type: ctr3Delta mutant, lacking the high-affinity plasma membrane copper transporter and mac1Delta, the copper-sensing transcription factor involved in regulation of high-affinity copper transport. Acute exposure of exponentially growing yeast cells revealed that TcPR-10 resistance is also enhanced in the Snq2 export permease-lacking mutant which has reduced intracellular presence of TcPR-10.

A genome survey of Moniliophthora perniciosa gives new insights into Witches' Broom Disease of cacao

Background

The basidiomycete fungus Moniliophthora perniciosa is the causal agent of Witches' Broom Disease (WBD) in cacao (Theobroma cacao). It is a hemibiotrophic pathogen that colonizes the apoplast of cacao's meristematic tissues as a biotrophic pathogen, switching to a saprotrophic lifestyle during later stages of infection. M. perniciosa, together with the related species M. roreri, are pathogens of aerial parts of the plant, an uncommon characteristic in the order Agaricales. A genome survey (1.9× coverage) of M. perniciosa was analyzed to evaluate the overall gene content of this phytopathogen.

Results

Genes encoding proteins involved in retrotransposition, reactive oxygen species (ROS) resistance, drug efflux transport and cell wall degradation were identified. The great number of genes encoding cytochrome P450 monooxygenases (1.15% of gene models) indicates that M. perniciosa has a great potential for detoxification, production of toxins and hormones; which may confer a high adaptive ability to the fungus. We have also discovered new genes encoding putative secreted polypeptides rich in cysteine, as well as genes related to methylotrophy and plant hormone biosynthesis (gibberellin and auxin). Analysis of gene families indicated that M. perniciosa have similar amounts of carboxylesterases and repertoires of plant cell wall degrading enzymes as other hemibiotrophic fungi. In addition, an approach for normalization of gene family data using incomplete genome data was developed and applied in M. perniciosa genome survey.

Conclusion

This genome survey gives an overview of the M. perniciosa genome, and reveals that a significant portion is involved in stress adaptation and plant necrosis, two necessary characteristics for a hemibiotrophic fungus to fulfill its infection cycle. Our analysis provides new evidence revealing potential adaptive traits that may play major roles in the mechanisms of pathogenicity in the M. perniciosa/cacao pathosystem.

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

BACKGROUND

Theobroma cacao L., is a tree originated from the tropical rainforest of South America. It is one of the major cash crops for many tropical countries. T. cacao is mainly produced on smallholdings, providing resources for 14 million farmers. Disease resistance and T. cacao quality improvement are two important challenges for all actors of cocoa and chocolate production. T. cacao is seriously affected by pests and fungal diseases, responsible for more than 40% yield losses and quality improvement, nutritional and organoleptic, is also important for consumers. An international collaboration was formed to develop an EST genomic resource database for cacao.

RESULTS

Fifty-six cDNA libraries were constructed from different organs, different genotypes and different environmental conditions. A total of 149,650 valid EST sequences were generated corresponding to 48,594 unigenes, 12,692 contigs and 35,902 singletons. A total of 29,849 unigenes shared significant homology with public sequences from other species.Gene Ontology (GO) annotation was applied to distribute the ESTs among the main GO categories.A specific information system (ESTtik) was constructed to process, store and manage this EST collection allowing the user to query a database.To check the representativeness of our EST collection, we looked for the genes known to be involved in two different metabolic pathways extensively studied in other plant species and important for T. cacao qualities: the flavonoid and the terpene pathways. Most of the enzymes described in other crops for these two metabolic pathways were found in our EST collection.A large collection of new genetic markers was provided by this ESTs collection.

CONCLUSION

This EST collection displays a good representation of the T. cacao transcriptome, suitable for analysis of biochemical pathways based on oligonucleotide microarrays derived from these ESTs. It will provide numerous genetic markers that will allow the construction of a high density gene map of T. cacao. This EST collection represents a unique and important molecular resource for T. cacao study and improvement, facilitating the discovery of candidate genes for important T. cacao trait variation.

Protein extraction for proteome analysis from cacao leaves and meristems, organs infected by Moniliophthora perniciosa, the causal agent of the witches' broom disease

Preparation of high-quality proteins from cacao vegetative organs is difficult due to very high endogenous levels of polysaccharides and polyphenols. In order to establish a routine procedure for the application of proteomic and biochemical analysis to cacao tissues, three new protocols were developed; one for apoplastic washing fluid (AWF) extraction, and two for protein extraction--under denaturing and nondenaturing conditions. The first described method allows a quick and easy collection of AWF--using infiltration-centrifugation procedure--that is representative of its composition in intact leaves according to the smaller symplastic contamination detected by the use of the hexose phosphate isomerase marker. Protein extraction under denaturing conditions for 2-DE was remarkably improved by the combination of chemically and physically modified processes including phenol, SDS dense buffer and sonication steps. With this protocol, high-quality proteins from cacao leaves and meristems were isolated, and for the first time well-resolved 1-DE and 2-DE protein patterns of cacao vegetative organs are shown. It also appears that sonication associated with polysaccharide precipitation using tert-butanol was a crucial step for the nondenaturing protein extraction and subsequent enzymatic activity detection. It is expected that the protocols described here could help to develop high-level proteomic and biochemical studies in cacao also being applicable to other recalcitrant plant tissues.

Identification of a second family of genes in Moniliophthora perniciosa, the causal agent of witches' broom disease in cacao, encoding necrosis-inducing proteins similar to cerato-platanins

The hemibiotrophic basidiomycete Moniliophthora perniciosa is the causal agent of witches' broom disease in cacao. This is a dimorphic species, with monokaryotic hyphae during the biotrophic phase, which is converted to dikaryotic mycelia during the saprophytic phase. The infection in pod is characterized by the formation of hypertrophic and hyperplasic tissues in the biotrophic phase, which is followed by necrosis and complete degradation of the organ. We found at least five sequences in the fungal genome encoding putative proteins similar to cerato-platanin (CP)-like proteins, a novel class of proteins initially found in the phytopathogen Ceratocystis fimbriata. One M. perniciosa CP gene (MpCP1) was expressed in vitro and proved to have necrosis-inducing ability in tobacco and cacao leaves. The protein is present in solution as dimers and is able to recover necrosis activity after heat treatment. Transcription analysis ex planta showed that MpCP1 is more expressed in biotrophic-like mycelia than saprotrophic mycelia. The necrosis profile presented is different from that caused by M. perniciosa necrosis and ethylene-inducing proteins (MpNEPs), another family of elicitors expressed by M. perniciosa. Remarkably, a mixture of MpCP1 with MpNEP2 led to a synergistic necrosis effect very similar to that found in naturally infected plants. This is the first report of a basidiomycete presenting both NEP1-like proteins (NLPs) and CPs in its genome.

Carbon source-dependent variation of acquired mutagen resistance of Moniliophthora perniciosa: similarities in natural and artificial systems

The basidiomycete Moniliophthora perniciosa causes Witches' Broom disease in Theobroma cacao. We studied the influence of carbon source on conditioning hyphae to oxidative stress agents (H(2)O(2), paraquat, 4NQO) and to UVC, toward the goal of assessing the ability of this pathogen to avoid plant defenses involving ROS. Cells exhibited increased resistance to H(2)O(2) when shifted from glucose to glycerol and from glycerol to glycerol. When exposed to paraquat, cells grown in fresh medium were always more resistant. Apparently glycerol and/or fresh media, but not old glucose media, up-regulate oxidative stress defenses in this fungus. For the mutagens UVC and 4NQO, whose prime action on DNA is not via ROS, change of carbon source did not elicit a clear change in sensitivity/resistance. These results correlate with expression of fungal genes that protect against ROS and with biochemical changes observed in infected cacao tissues, where glycerol and high amounts of ROS have been detected in green brooms.

Production of calcium oxalate crystals by the basidiomycete Moniliophthora perniciosa, the causal agent of witches' broom disease of Cacao

Oxalic acid has been shown as a virulence factor for some phytopathogenic fungi, removing calcium from pectin and favoring plant cell wall degradation. Recently, it was published that calcium oxalate accumulates in infected cacao tissues during the progression of Witches' Broom disease (WBD). In the present work we report that the hemibiotrophic basidiomycete Moniliophthora perniciosa, the causal agent of WBD, produces calcium oxalate crystals. These crystals were initially observed by polarized light microscopy of hyphae growing on a glass slide, apparently being secreted from the cells. The analysis was refined by Scanning electron microscopy and the compositon of the crystals was confirmed by energy-dispersive x-ray spectrometry. The production of oxalate by M. perniciosa was reinforced by the identification of a putative gene coding for oxaloacetate acetylhydrolase, which catalyzes the hydrolysis of oxaloacetate to oxalate and acetate. This gene was shown to be expressed in the biotrophic-like mycelia, which in planta occupy the intercellular middle-lamella space, a region filled with pectin. Taken together, our results suggest that oxalate production by M. perniciosa may play a role in the WBD pathogenesis mechanism.