Genetic Diversity and Natural Population Structure of Cacao (Theobroma cacao L.) from the Brazilian Amazon Evaluated by Microsatellite Markers

Genetic diversity and verification of plant material compliance of Cocoa (Theobroma cacao L.) in the Barombi-Kang Regional variety trial

Cocoa (Theobroma cacao L.) is a pivotal agricultural commodity in Cameroon, which ranks as one of the top five global cocoa producers. This study focused on evaluating the genetic diversity and verifying the plant material compliance of cocoa genotypes in the Barombi-Kang Regional variety trial, employing 12 highly polymorphic SSR markers. A comprehensive analysis of 318 hybrid families and 15 parental genotypes was conducted, which revealed extensive genetic variability. The study found an average polymorphic information content (PIC) of 0.72 for hybrids and 0.68 for parents, alongside observed heterozygosity rates of 0.54 and 0.42, respectively, indicating a rich genetic reservoir. Significantly, an 18.55% labeling error rate was identified, underscoring prevalent issues in germplasm management that could impact the efficacy of breeding programs. These errors highlight the critical need for enhanced genetic verification protocols to ensure the accuracy and reliability of plant materials used in breeding. The genetic analysis also demonstrated substantial allelic richness with the hybrids showing an average of 72 alleles per locus, suggesting a high capacity for selection within the breeding pool. The data from this study not only reinforce the potential for genetic improvement of cocoa in Cameroon but also provide crucial insights into the genetic structure and population dynamics within the trial. Addressing the genetic and management challenges identified could lead to the development of superior cocoa varieties, enhancing yield, disease resistance, and environmental stress tolerance, thereby contributing to the sustainable advancement of the cocoa industry in Cameroon and beyond.

Local working collections as the foundation for an integrated conservation of Theobroma cacao L. in Latin America

The intraspecific diversity of cacao has traditionally been preserved in genebanks. However, these establishments face various challenges, notably insufficient funding, accession redundancy, misidentification and lack of wild cacao population samples. In natural environments, it is expected that unknown varieties of cacao may still be found, but wild populations of cacao are increasingly threatened by climate change, deforestation, habitat loss, land use changes and poor knowledge. Farmers also retain diversity, but on-farm conservation is affected by geopolitical, economic, management and cultural issues, that are influenced at multiple scales, from the household to the international market. Taking separately, ex situ, in situ and on-farm conservation have not achieved adequate conservation fostering the inclusion of all stakeholders and the broad use of cacao diversity. We analyze the use of the traditional conservation strategies (ex situ, in situ and on-farm) and propose an integrated approach based on local working collections to secure cacao diversity in the long term. We argue that national conservation networks should be implemented in countries of origin to simultaneously maximize alpha (diversity held in any given working collection), beta (the change in diversity between working collections in different regions) and gamma diversity (overall diversity in a country).

Genome-wide SNP genotyping as a simple and practical tool to accelerate the development of inbred lines in outbred tree species: An example in cacao (Theobroma cacao L.)

Cacao is a globally important crop with a long history of domestication and selective breeding. Despite the increased use of elite clones by cacao farmers, worldwide plantations are established mainly using hybrid progeny material derived from heterozygous parents, therefore displaying high tree-to-tree variability. The deliberate development of hybrids from advanced inbred lines produced by successive generations of self-pollination has not yet been fully considered in cacao breeding. This is largely due to the self-incompatibility of the species, the long generation cycles (3–5 years) and the extensive trial areas needed to accomplish the endeavor. We propose a simple and accessible approach to develop inbred lines based on accelerating the buildup of homozygosity based on regular selfing assisted by genome-wide SNP genotyping. In this study we genotyped 90 clones from the Brazilian CEPEC´s germplasm collection and 49 inbred offspring of six S₁ or S₂ cacao families derived from self-pollinating clones CCN-51, PS-13.19, TSH-1188 and SIAL-169. A set of 3,380 SNPs distributed across the cacao genome were interrogated on the EMBRAPA multi-species 65k Infinium chip. The 90 cacao clones showed considerable variation in genome-wide SNP homozygosity (mean 0.727± 0.182) and 19 of them with homozygosity ≥90%. By assessing the increase in homozygosity across two generations of self-pollinations, SNP data revealed the wide variability in homozygosity within and between S₁ and S₂ families. Even in small families (<10 sibs), individuals were identified with up to ~1.5 standard deviations above the family mean homozygosity. From baseline homozygosities of 0.476 and 0.454, offspring with homozygosities of 0.862 and 0.879 were recovered for clones TSH-1188 and CCN-51 respectively, in only two generations of selfing (81–93% increase). SNP marker assisted monitoring and selection of inbred individuals can be a practical tool to optimize and accelerate the development of inbred lines of outbred tree species. This approach will allow a faster and more accurate exploitation of hybrid breeding strategies in cacao improvement programs and potentially in other perennial fruit and forest trees.

Structural and Functional Genomics of the Resistance of Cacao to Phytophthora palmivora

Black pod disease, caused by Phytophthora spp., is one of the main diseases that attack cocoa plantations. This study validated, by association mapping, 29 SSR molecular markers flanking to QTL (Quantitative Trait Loci) associated with Phytophthora palmivora Butler (Butler) (PP) resistance, in three local ancient varieties of the Bahia (Comum, Pará, and Maranhão), varieties that have a high potential in the production of gourmet chocolate. Four SSR loci associated with resistance to PP were detected, two on chromosome 8, explaining 7.43% and 3.72% of the Phenotypic Variation (%PV), one on chromosome 2 explaining 2.71%PV and one on chromosome 3 explaining 1.93%PV. A functional domains-based annotation was carried out, in two Theobroma cacao (CRIOLLO and MATINA) reference genomes, of 20 QTL regions associated with cocoa resistance to the pathogen. It was identified 164 (genome CRIOLLO) and 160 (genome MATINA) candidate genes, hypothetically involved in the recognition and activation of responses in the interaction with the pathogen. Genomic regions rich in genes with Coiled-coils (CC), nucleotide binding sites (NBS) and Leucine-rich repeat (LRR) domains were identified on chromosomes 1, 3, 6, 8, and 10, likewise, regions rich in Receptor-like Kinase domain (RLK) and Ginkbilobin2 (GNK2) domains were identified in chromosomes 4 and 6.

Geographic Patterns of Genetic Variation among Cacao (Theobroma cacao L.) Populations Based on Chloroplast Markers

The cacao tree (Theobroma cacao L.) is native to the Amazon basin and widely cultivated in the tropics to produce seeds, the valuable raw material for the chocolate industry. Conservation of cacao genetic resources and their availability for breeding and production programs are vital for securing cacao supply. However, relatively little is still known about the phylogeographic structure of natural cacao populations. We studied the geographic distribution of cpDNA variation in different populations representing natural cacao stands, cacao farms in Ecuador, and breeding populations. We used six earlier published cacao chloroplast microsatellite markers to genotype 233 cacao samples. In total, 23 chloroplast haplotypes were identified. The highest variation of haplotypes was observed in western Amazonia including geographically restricted haplotypes. Two observed haplotypes were widespread across the Amazon basin suggesting long distance seed dispersal from west to east in Amazonia. Most cacao genetic groups identified earlier using nuclear SSRs are associated with specific chloroplast haplotypes. A single haplotype was common in selections representing cacao plantations in west Ecuador and reference Trinitario accessions. Our results can be used to determine the chloroplast diversity of accessions and in combination with phenotypic assessments can help to select geographically distinctive varieties for cacao breeding programs.

Genetic structuring in a Neotropical palm analyzed through an Andean orogenesis-scenario

Andean orogenesis has driven the development of very high plant diversity in the Neotropics through its impact on landscape evolution and climate. The analysis of the intraspecific patterns of genetic structure in plants would permit inferring the effects of Andean uplift on the evolution and diversification of Neotropical flora. In this study, using microsatellite markers and Bayesian clustering analyses, we report the presence of four genetic clusters for the palm Oenocarpus bataua var. bataua which are located within four biogeographic regions in northwestern South America: (a) Chocó rain forest, (b) Amotape-Huancabamba Zone, (c) northwestern Amazonian rain forest, and (d) southwestern Amazonian rain forest. We hypothesize that these clusters developed following three genetic diversification events mainly promoted by Andean orogenic events. Additionally, the distinct current climate dynamics among northwestern and southwestern Amazonia may maintain the genetic diversification detected in the western Amazon basin. Genetic exchange was identified between the clusters, including across the Andes region, discarding the possibility of any cluster to diversify as a distinct intraspecific variety. We identified a hot spot of genetic diversity in the northern Peruvian Amazon around the locality of Iquitos. We also detected a decrease in diversity with distance from this area in westward and southward direction within the Amazon basin and the eastern Andean foothills. Additionally, we confirmed the existence and divergence of O. bataua var. bataua from var. oligocarpus in northern South America, possibly expanding the distributional range of the latter variety beyond eastern Venezuela, to the central and eastern Andean cordilleras of Colombia. Based on our results, we suggest that Andean orogenesis is the main driver of genetic structuring and diversification in O. bataua within northwestern South America.

Cocoa agronomy, quality, nutritional, and health aspects

The history of cocoa and chocolate including the birth and the expansion of the chocolate industry was described. Recent developments in the industry and cocoa economy were briefly depicted. An overview of the classification of cacao as well as studies on phenotypic and genetic diversity was presented. Cocoa agronomic practices including traditional and modern propagation techniques were reviewed. Nutrition-related health benefits derived from cocoa consumption were listed and widely reviewed. The specific action of cocoa antioxidants was compared to those of teas and wines. Effects of adding milk to chocolate and chocolate drinks versus bioavailability of cocoa polyphenols were discussed. Finally, flavor, sensory, microbiological, and toxicological aspects of cocoa consumption were presented.

Genetic Structure and Molecular Diversity of Cacao Plants Established as Local Varieties for More than Two Centuries: The Genetic History of Cacao Plantations in Bahia, Brazil

Bahia is the most important cacao-producing state in Brazil, which is currently the sixth-largest country worldwide to produce cacao seeds. In the eighteenth century, the Comum, Pará and Maranhão varieties of cacao were introduced into southern Bahia, and their descendants, which are called ‘Bahian cacao’ or local Bahian varieties, have been cultivated for over 200 years. Comum plants have been used to start plantations in African countries and extended as far as countries in South Asia and Oceania. In Brazil, two sets of clones selected from Bahian varieties and their mutants, the Agronomic Institute of East (SIAL) and Bahian Cacao Institute (SIC) series, represent the diversity of Bahian cacao in germplasm banks. Because the genetic diversity of Bahian varieties, which is essential for breeding programs, remains unknown, the objective of this work was to assess the genetic structure and diversity of local Bahian varieties collected from farms and germplasm banks. To this end, 30 simple sequence repeat (SSR) markers were used to genotype 279 cacao plants from germplasm and local farms. The results facilitated the identification of 219 cacao plants of Bahian origin, and 51 of these were SIAL or SIC clones. Bahian cacao showed low genetic diversity. It could be verified that SIC and SIAL clones do not represent the true diversity of Bahian cacao, with the greatest amount of diversity found in cacao trees on the farms. Thus, a core collection to aid in prioritizing the plants to be sampled for Bahian cacao diversity is suggested. These results provide information that can be used to conserve Bahian cacao plants and applied in breeding programs to obtain more productive Bahian cacao with superior quality and tolerance to major diseases in tropical cacao plantations worldwide.

Extraction of high-quality DNA from ethanol-preserved tropical plant tissues

Background

Proper conservation of plant samples, especially during remote field collection, is essential to assure quality of extracted DNA. Tropical plant species contain considerable amounts of secondary compounds, such as polysaccharides, phenols, and latex, which affect DNA quality during extraction. The suitability of ethanol (96% v/v) as a preservative solution prior to DNA extraction was evaluated using leaves of Jatropha curcas and other tropical species.

Results

Total DNA extracted from leaf samples stored in liquid nitrogen or ethanol from J. curcas and other tropical species (Theobroma cacao, Coffea arabica, Ricinus communis, Saccharum spp., and Solanum lycopersicon) was similar in quality, with high-molecular-weight DNA visualized by gel electrophoresis. DNA quality was confirmed by digestion with EcoRI or HindIII and by amplification of the ribosomal gene internal transcribed spacer region. Leaf tissue of J. curcas was analyzed by light and transmission electron microscopy before and after exposure to ethanol. Our results indicate that leaf samples can be successfully preserved in ethanol for long periods (30 days) as a viable method for fixation and conservation of DNA from leaves. The success of this technique is likely due to reduction or inactivation of secondary metabolites that could contaminate or degrade genomic DNA.

Conclusions

Tissue conservation in 96% ethanol represents an attractive low-cost alternative to commonly used methods for preservation of samples for DNA extraction. This technique yields DNA of equivalent quality to that obtained from fresh or frozen tissue.

Gene flow and genetic diversity in cultivated and wild cacao (Theobroma cacao) in Bolivia

PREMISE OF THE STUDY

The role of pollen flow within and between cultivated and wild tropical crop species is little known. To study the pollen flow of cacao, we estimated the degree of self-pollination and pollen dispersal distances as well as gene flow between wild and cultivated cacao (Theobroma cacao L.).

METHODS

We studied pollen flow and genetic diversity of cultivated and wild cacao populations by genotyping 143 wild and 86 cultivated mature plants and 374 seedlings raised from 19 wild and 25 cultivated trees at nine microsatellite loci.

KEY RESULTS

A principal component analysis distinguished wild and cultivated cacao trees, supporting the notion that Bolivia harbors truly wild cacao populations. Cultivated cacao had a higher level of genetic diversity than wild cacao, presumably reflecting the varied origin of cultivated plants. Both cacao types had high outcrossing rates, but the paternity analysis revealed 7-14% self-pollination in wild and cultivated cacao. Despite the tiny size of the pollinators, pollen was transported distances up to 3 km; wild cacao showed longer distances (mean = 922 m) than cultivated cacao (826 m). Our data revealed that 16-20% of pollination events occurred between cultivated and wild populations.

CONCLUSIONS

We found evidence of self-pollination in both wild and cultivated cacao. Pollination distances are larger than those typically reported in tropical understory tree species. The relatively high pollen exchange from cultivated to wild cacao compromises genetic identity of wild populations, calling for the protection of extensive natural forest tracts to protect wild cacao in Bolivia.

Development of novel microsatellites from Moniliophthora perniciosa, causal agent of the witches' broom disease of Theobroma cacao

Moniliophthora perniciosa is the causal agent of the witches' broom disease of cacao. Based on available genomic sequences, we identified 30 new microsatellite loci, which were analysed using 50 isolates from four populations sampled over a wide geographical area in Brazil, including three populations from the Amazon, the fungal putative centre of diversity, plus one from Bahia. Nine loci were polymorphic, with an average of 2.9 alleles per locus. The level of polymorphism observed was low, but these markers may allow the evaluation of pathogen diversity and the establishment of molecular standards for isolate fingerprinting to support cacao breeding.

ARG-based genome-wide analysis of cacao cultivars

Background

Ancestral recombinations graph (ARG) is a topological structure that captures the relationship between the extant genomic sequences in terms of genetic events including recombinations. IRiS is a system that estimates the ARG on sequences of individuals, at genomic scales, capturing the relationship between these individuals of the species. Recently, this system was used to estimate the ARG of the recombining X Chromosome of a collection of human populations using relatively dense, bi-allelic SNP data.

Results

While the ARG is a natural model for capturing the inter-relationship between a single chromosome of the individuals of a species, it is not immediately apparent how the model can utilize whole-genome (across chromosomes) diploid data. Also, the sheer complexity of an ARG structure presents a challenge to graph visualization techniques. In this paper we examine the ARG reconstruction for (1) genome-wide or multiple chromosomes, (2) multi-allelic and (3) extremely sparse data. To aid in the visualization of the results of the reconstructed ARG, we additionally construct a much simplified topology, a classification tree, suggested by the ARG.

As the test case, we study the problem of extracting the relationship between populations of Theobroma cacao. The chocolate tree is an outcrossing species in the wild, due to self-incompatibility mechanisms at play. Thus a principled approach to understanding the inter-relationships between the different populations must take the shuffling of the genomic segments into account. The polymorphisms in the test data are short tandem repeats (STR) and are multi-allelic (sometimes as high as 30 distinct possible values at a locus). Each is at a genomic location that is bilaterally transmitted, hence the ARG is a natural model for this data. Another characteristic of this plant data set is that while it is genome-wide, across 10 linkage groups or chromosomes, it is very sparse, i.e., only 96 loci from a genome of approximately 400 megabases. The results are visualized both as MDS plots and as classification trees. To evaluate the accuracy of the ARG approach, we compare the results with those available in literature.

Conclusions

We have extended the ARG model to incorporate genome-wide (ensemble of multiple chromosomes) data in a natural way. We present a simple scheme to implement this in practice. Finally, this is the first time that a plant population data set is being studied by estimating its underlying ARG. We demonstrate an overall precision of 0.92 and an overall recall of 0.93 of the ARG-based classification, with respect to the gold standard. While we have corroborated the classification of the samples with that in literature, this opens the door to other potential studies that can be made on the ARG.

Insight into the wild origin, migration and domestication history of the fine flavour Nacional Theobroma cacao L. variety from Ecuador

Ecuador's economic history has been closely linked to Theobroma cacao L cultivation, and specifically to the native fine flavour Nacional cocoa variety. The original Nacional cocoa trees are presently in danger of extinction due to foreign germplasm introductions. In a previous work, a few non-introgressed Nacional types were identified as potential founders of the modern Ecuadorian cocoa population, but so far their origin could not be formally identified. In order to determine the putative centre of origin of Nacional and trace its domestication history, we used 80 simple sequence repeat (SSR) markers to analyse the relationships between these potential Nacional founders and 169 wild and cultivated cocoa accessions from South and Central America. The highest genetic similarity was observed between the Nacional pool and some wild genotypes from the southern Amazonian region of Ecuador, sampled along the Yacuambi, Nangaritza and Zamora rivers in Zamora Chinchipe province. This result was confirmed by a parentage analysis. Based on our results and on data about pre-Columbian civilization and Spanish colonization history of Ecuador, we determined, for the first time, the possible centre of origin and migration events of the Nacional variety from the Amazonian area until its arrival in the coastal provinces. As large unexplored forest areas still exist in the southern part of the Ecuadorian Amazonian region, our findings could provide clues as to where precious new genetic resources could be collected, and subsequently used to improve the flavour and disease resistance of modern Ecuadorian cocoa varieties.

Present spatial diversity patterns of Theobroma cacao L. in the neotropics reflect genetic differentiation in pleistocene refugia followed by human-influenced dispersal

Cacao (Theobroma cacao L.) is indigenous to the Amazon basin, but is generally believed to have been domesticated in Mesoamerica for the production of chocolate beverage. However, cacao’s distribution of genetic diversity in South America is also likely to reflect pre-Columbian human influences that were superimposed on natural processes of genetic differentiation. Here we present the results of a spatial analysis of the intra-specific diversity of cacao in Latin America, drawing on a dataset of 939 cacao trees genotypically characterized by means of 96 SSR markers. To assess continental diversity patterns we performed grid-based calculations of allelic richness, Shannon diversity and Nei gene diversity, and distinguished different spatially coherent genetic groups by means of cluster analysis. The highest levels of genetic diversity were observed in the Upper Amazon areas from southern Peru to the Ecuadorian Amazon and the border areas between Colombia, Peru and Brazil. On the assumption that the last glaciation (22,000–13,000 BP) had the greatest pre-human impact on the current distribution and diversity of cacao, we modeled the species’ Pleistocene niche suitability and overlaid this with present-day diversity maps. The results suggest that cacao was already widely distributed in the Western Amazon before the onset of glaciation. During glaciations, cacao populations were likely to have been restricted to several refugia where they probably underwent genetic differentiation, resulting in a number of genetic clusters which are representative for, or closest related to, the original wild cacao populations. The analyses also suggested that genetic differentiation and geographical distribution of a number of other clusters seem to have been significantly affected by processes of human management and accompanying genetic bottlenecks. We discuss the implications of these results for future germplasm collection and in situ, on farm and ex situ conservation of cacao.

Genetic diversity of Lagerstroemia (Lythraceae) species assessed by simple sequence repeat markers

Lagerstroemia (crape myrtle) are famous ornamental plants with large pyramidal racemes, long flower duration, and diverse colors. However, little is known about the genetic structure and diversity of germplasm in Lagerstroemia. We genotyped 81 L. indica cultivars, five other species of Lagerstroemia, and 10 interspecific hybrids using 30 simple sequence repeat markers; 275 alleles were generated with a mean of nine alleles per locus. The mean polymorphism information content value, a measure of gene diversity, was 0.63, with a range from 0.25 to 0.86. The mean observed heterozygosity (0.51) tended to be lower than the mean expected heterozygosity (0.67). The mean F-statistics (F(ST), F(IS), and F(IT)) were 0.05, 0.20, and 0.24, respectively, indicating a high level of genetic variation among cultivars. Clustering analysis based on genetic distance divided the 96 genotypes into three distinct groups, which corresponded with their genetic backgrounds and geographic regions. L. indica cultivars and the other five L. species were grouped into different sub-clusters. Chinese and North American cultivars were divided into different clusters. These data about the genetic relationship among cultivars demonstrated the potential value of L. indica cultivars and other Lagerstroemia species for widening the genetic basis of breeding programs for this ornamental flower.

Development of microsatellite primers for Jatropha curcas (Euphorbiaceae) and transferability to congeners

PREMISE OF THE STUDY

Microsatellite primers were developed for Jatropha curcas (Euphorbiaceae), a tree species with large potential for biofuel production, to investigate its natural genetic diversity and mating system to facilitate the establishment of tree improvement and conservation programs.

METHODS AND RESULTS

Using a protocol for genomic library enrichment, 104 clones containing 195 repeat motifs were identified. Primer pairs were developed for 40 microsatellite loci and validated in 41 accessions of J. curcas from six provenances. Nine loci were polymorphic revealing from two to eight alleles per locus, and six primers were able to amplify alleles in the congeners J. podagrica, J. pohliana, and J. gossypifolia, but not in other Euphorbiaceae species, such as Hevea brasiliensis, Manihot esculenta, or Ricinus communis.

CONCLUSIONS

The primers developed here revealed polymorphic loci that are suitable for genetic diversity and structure, mating system, and gene flow studies in J. curcas, and some congeners.

Understanding the genetic diversity, spatial genetic structure and mating system at the hierarchical levels of fruits and individuals of a continuous Theobroma cacao population from the Brazilian Amazon

Understanding the mating patterns of populations of tree species is a key component of ex situ genetic conservation. In this study, we analysed the genetic diversity, spatial genetic structure (SGS) and mating system at the hierarchical levels of fruits and individuals as well as pollen dispersal patterns in a continuous population of Theobroma cacao in Pará State, Brazil. A total of 156 individuals in a 0.56 ha plot were mapped and genotyped for nine microsatellite loci. For the mating system analyses, 50 seeds were collected from nine seed trees by sampling five fruits per tree (10 seeds per fruit). Among the 156 individuals, 127 had unique multilocus genotypes, and the remaining were clones. The population was spatially aggregated; it demonstrated a significant SGS up to 15 m that could be attributed primarily to the presence of clones. However, the short seed dispersal distance also contributed to this pattern. Population matings occurred mainly via outcrossing, but selfing was observed in some seed trees, which indicated the presence of individual variation for self-incompatibility. The matings were also correlated, especially within (Ρ(p(m))=0.607) rather than among the fruits (Ρ(p(m))=0.099), which suggested that a small number of pollen donors fertilised each fruit. The paternity analysis suggested a high proportion of pollen migration (61.3%), although within the plot, most of the pollen dispersal encompassed short distances (28 m). The determination of these novel parameters provides the fundamental information required to establish long-term ex situ conservation strategies for this important tropical species.

Geographic and genetic population differentiation of the Amazonian chocolate tree (Theobroma cacao L)

Numerous collecting expeditions of Theobroma cacao L. germplasm have been undertaken in Latin-America. However, most of this germplasm has not contributed to cacao improvement because its relationship to cultivated selections was poorly understood. Germplasm labeling errors have impeded breeding and confounded the interpretation of diversity analyses. To improve the understanding of the origin, classification, and population differentiation within the species, 1241 accessions covering a large geographic sampling were genotyped with 106 microsatellite markers. After discarding mislabeled samples, 10 genetic clusters, as opposed to the two genetic groups traditionally recognized within T. cacao, were found by applying Bayesian statistics. This leads us to propose a new classification of the cacao germplasm that will enhance its management. The results also provide new insights into the diversification of Amazon species in general, with the pattern of differentiation of the populations studied supporting the palaeoarches hypothesis of species diversification. The origin of the traditional cacao cultivars is also enlightened in this study.

Genes differentially expressed in Theobroma cacao associated with resistance to witches' broom disease caused by Crinipellis perniciosa

SUMMARY The basidiomycete Crinipellis perniciosa is the causal agent of witches' broom disease of Theobroma cacao (cocoa). Hypertrophic growth of infected buds ('brooms') is the most dramatic symptom, but the main economic losses derive from pod infection. To identify cocoa genes differentially expressed during the early stages of infection, two cDNA libraries were constructed using the suppression subtractive hybridization (SSH) approach. Subtraction hybridization was conducted between cDNAs from infected shoot-tips of the susceptible genotype 'ICS 39' and the resistant 'CAB 214', in both directions. A total of 187 unique sequences were obtained, with 83 from the library enriched for the susceptible 'ICS 39' sequences, and 104 for the resistant 'CAB 214'. By homology search and ontology analyses, the identified sequences were mainly putatively categorized as belonging to 'signal transduction', 'response to biotic and abiotic stress', 'metabolism', 'RNA and DNA metabolism', 'protein metabolism' and 'cellular maintenance' classes. Quantitative reverse transcription amplification (RT-qPCR) of 23 transcripts identified as differentially expressed between genotypes revealed distinct kinetics of gene up-regulation at the asymptomatic stage of the disease. Expression induction in the susceptible 'ICS 39' in response to C. perniciosa was delayed and limited, while in 'CAB 214' there was a quicker and more intense reaction, with two peaks of gene induction at 48 and 120 h after inoculation, corresponding to morphological and biochemical changes previously described during colonization. Similar differences in gene induction were validated for another resistant genotype ('CAB 208') in an independent experiment. Validation of these genes corroborated similar hypothetical mechanisms of resistance described in other pathosystems.