A revisited history of cacao domestication in pre-Columbian times revealed by archaeogenomic approaches.

Humans have a long history of transporting and trading plants, contributing to the evolution of domesticated plants. Theobroma cacao originated in the Neotropics from South America. However, little is known about its domestication and use in these regions. In this study, ceramic residues from a large sample of pre-Columbian cultures from South and Central America were analyzed using archaeogenomic and biochemical approaches. Here we show, for the first time, the widespread use of cacao in South America out of its native Amazonian area of origin, extending back 5000 years, likely supported by cultural interactions between the Amazon and the Pacific coast. We observed that strong genetic mixing between geographically distant cacao populations occurred as early as the middle Holocene, in South America, driven by humans, favoring the adaptation of T. cacao to new environments. This complex history of cacao domestication is the basis of today's cacao tree populations and its knowledge can help us better manage their genetic resources.

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

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

Diversity and determinants of bitterness, astringency, and fat content in cultivated Nacional and native Amazonian cocoa accessions from Ecuador.

Cocoa (Theobroma cacao L.) is the only tree that can produce cocoa. Cocoa beans are highly sought after by chocolate makers to produce chocolate. Cocoa can be fine aromatic, characterized by floral and fruity notes, or it can be described as standard cocoa with a more pronounced cocoa aroma and bitterness. In this study, the genetic and biochemical determinants of sensorial notes and nonvolatile compounds related to bitterness, astringency, fat content, and protein content will be investigated in two populations: a cultivated modern Nacional population and a population of cocoa accessions collected recently in the Ecuadorian South Amazonia area of origin of the Nacional ancestral variety. For this purpose, a genome-wide association study (GWAS) was carried out on both populations, with results of biochemical compounds evaluated by near-infrared spectroscopy (NIRS) assays and with sensory evaluations. Twenty areas of associations were detected for sensorial data especially bitterness and astringency. Fifty-three areas of associations were detected linked to nonvolatile compounds. A total of 81 candidate genes could be identified in the areas of the association.

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

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

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

MAIN CONCLUSION: 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.

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.