Application of Genome Wide Association and Genomic Prediction for Improvement of Cacao Productivity and Resistance to Black and Frosty Pod Diseases

Three de novo assembled wild cacao genomes from the Upper Amazon

Theobroma cacao, the chocolate tree, is indigenous to the Amazon basin, the greatest biodiversity hotspot on earth. Recent advancement in plant genomics highlights the importance of de novo sequencing of multiple reference genomes to capture the genome diversity present in different cacao populations. In this study, three high-quality chromosome-level genomes of wild cacao were constructed, de novo assembled with HiFi long reads sequencing, and scaffolded using a reference-free strategy. These genomes represent the three most important genetic clusters of cacao trees from the Upper Amazon region. The three wild cacao genomes were compared with two reference genomes of domesticated cacao. The five cacao genetic clusters were inferred to have diverged in the early and middle Pleistocene period, approximately 1.83-0.69 million years ago. The results shown here serve as an example of understanding how the Amazonian biodiversity was developed. The three wild cacao genomes provide valuable resources for studying genetic diversity and advancing genetic improvement of this species.

Genomic selection in tropical perennial crops and plantation trees: a review

To overcome the multiple challenges currently faced by agriculture, such as climate change and soil deterioration, more efficient plant breeding strategies are required. Genomic selection (GS) is crucial for the genetic improvement of quantitative traits, as it can increase selection intensity, shorten the generation interval, and improve selection accuracy for traits that are difficult to phenotype. Tropical perennial crops and plantation trees are of major economic importance and have consequently been the subject of many GS articles. In this review, we discuss the factors that affect GS accuracy (statistical models, linkage disequilibrium, information concerning markers, relatedness between training and target populations, the size of the training population, and trait heritability) and the genetic gain expected in these species. The impact of GS will be particularly strong in tropical perennial crops and plantation trees as they have long breeding cycles and constrained selection intensity. Future GS prospects are also discussed. High-throughput phenotyping will allow constructing of large training populations and implementing of phenomic selection. Optimized modeling is needed for longitudinal traits and multi-environment trials. The use of multi-omics, haploblocks, and structural variants will enable going beyond single-locus genotype data. Innovative statistical approaches, like artificial neural networks, are expected to efficiently handle the increasing amounts of heterogeneous multi-scale data. Targeted recombinations on sites identified from profiles of marker effects have the potential to further increase genetic gain. GS can also aid re-domestication and introgression breeding. Finally, GS consortia will play an important role in making the best of these opportunities.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01326-4.

Cacao (Theobroma cacao L.) Response to Water Stress: Physiological Characterization and Antioxidant Gene Expression Profiling in Commercial Clones

The increase in events associated with drought constraints plant growth and crop performance. Cacao (Theobroma cacao L.) is sensitive to water deficit stress (DS), which limits productivity. The aim of this research was to characterise the response of seven (CCN51, FEAR5, ICS1, ICS60, ICS95, EET8, and TSH565) commercially important cacao clones to severe and temporal water deficit stress. Ten-month-old cacao trees were submitted to two treatments: well-watered and water-stressed until the leaf water potential (Ψ leaf) reached values between -3.0 and -3.5 MPa. The effects of hydric stress on water relations, gas exchange, photochemical activity, membrane integrity and oxidative stress-related gene expression were evaluated. All clones showed decreases in Ψ leaf, but TSH565 had a higher capacity to maintain water homeostasis in leaves. An initial response phase consisted of stomatal closure, a general mechanism to limit water loss: as a consequence, the photosynthetic rate dropped by approximately 98% on average. In some clones, the photosynthetic rate reached negative values at the maximum stress level, evidencing photorespiration and was confirmed by increased intracellular CO2. A second and photosynthetically limited phase was characterized by a drop in PSII quantum efficiency, which affected all clones. On average, all clones were able to recover after 4 days of rewatering. Water deficit triggered oxidative stress at the early phase, as evidenced by the upregulation of oxidative stress markers and genes encoding ROS scavenging enzymes. The effects of water deficit stress on energy metabolism were deduced given the upregulation of fermentative enzyme-coding genes. Altogether, our results suggest that the EET8 clone was the highest performing under water deficit while the ICS-60 clone was more susceptible to water stress. Importantly, the activation of the antioxidant system and PSII repair mechanism seem to play key roles in the observed differences in tolerance to water deficit stress among clones.

Favorable Conditions for Genomic Evaluation to Outperform Classical Pedigree Evaluation Highlighted by a Proof-of-Concept Study in Poplar

Forest trees like poplar are particular in many ways compared to other domesticated species. They have long juvenile phases, ongoing crop-wild gene flow, extensive outcrossing, and slow growth. All these particularities tend to make the conduction of breeding programs and evaluation stages costly both in time and resources. Perennials like trees are therefore good candidates for the implementation of genomic selection (GS) which is a good way to accelerate the breeding process, by unchaining selection from phenotypic evaluation without affecting precision. In this study, we tried to compare GS to pedigree-based traditional evaluation, and evaluated under which conditions genomic evaluation outperforms classical pedigree evaluation. Several conditions were evaluated as the constitution of the training population by cross-validation, the implementation of multi-trait, single trait, additive and non-additive models with different estimation methods (G-BLUP or weighted G-BLUP). Finally, the impact of the marker densification was tested through four marker density sets. The population under study corresponds to a pedigree of 24 parents and 1,011 offspring, structured into 35 full-sib families. Four evaluation batches were planted in the same location and seven traits were evaluated on 1 and 2 years old trees. The quality of prediction was reported by the accuracy, the Spearman rank correlation and prediction bias and tested with a cross-validation and an independent individual test set. Our results show that genomic evaluation performance could be comparable to the already well-optimized pedigree-based evaluation under certain conditions. Genomic evaluation appeared to be advantageous when using an independent test set and a set of less precise phenotypes. Genome-based methods showed advantages over pedigree counterparts when ranking candidates at the within-family levels, for most of the families. Our study also showed that looking at ranking criteria as Spearman rank correlation can reveal benefits to genomic selection hidden by biased predictions.

Genome-Wide Association Study Reveals Novel Candidate Genes Associated with Productivity and Disease Resistance to Moniliophthora spp. in Cacao (Theobroma cacao L.)

Cacao (Theobroma cacao L.), the source of chocolate, is one of the most important commodity products worldwide that helps improve the economic livelihood of farmers. Diseases like frosty pod rot caused by Moniliophthora roreri and witches’ broom caused by Moniliophthora perniciosa limit the cacao productivity, this can be solved by using resistant varieties. In the current study, we sequenced 229 cacao accessions using genotyping-by-sequencing to examine the genetic diversity and population structure employing 9,003 and 8,131 single nucleotide polymorphisms recovered by mapping against two cacao genomes (Criollo B97-61/B2 v2 and Matina 1-6 v1.1). In the phenotypic evaluation, three promising accessions for productivity and 10 with good tolerance to the frosty pod rot and witches’ broom diseases were found. A genome-wide association study was performed on 102 accessions, discovering two genes associated with productivity and seven to disease resistance. The results enriched the knowledge of the genetic regions associated with important cacao traits that can have significant implications for conservation and breeding strategies like marker-assisted selection.

Identification of Climate and Genetic Factors That Control Fat Content and Fatty Acid Composition of Theobroma cacao L. Beans

The main ingredients of chocolate are usually cocoa powder, cocoa butter, and sugar. Both the powder and the butter are extracted from the beans of the cacao tree (Theobroma cacao L.). The cocoa butter represents the fat in the beans and possesses a unique fatty acid profile that results in chocolate's characteristic texture and mouthfeel. Here, we used a linkage mapping population and phenotypic data of 3,292 samples from 420 progeny which led to the identification of 27 quantitative trait loci (QTLs) for fatty acid composition and six QTLs for fat content. Progeny showed extensive variation in fat levels and composition, with the level of palmitic acid negatively correlated to the sum of stearic acid, oleic acid, and linoleic acid. A major QTL explaining 24% of the relative level of palmitic acid was mapped to the distal end of chromosome 4, and those higher levels of palmitic acid were associated with the presence of a haplotype from the "TSH 1188" parent in the progeny. Within this region of chromosome 4 is the Thecc1EG017405 gene, an orthologue and isoform of the stearoyl-acyl carrier protein (ACP) desaturase (SAD) gene in plants, which is involved in fatty acid biosynthesis. Besides allelic differences, we also show that climate factors can change the fatty acid composition in the beans, including a significant positive correlation between higher temperatures and the higher level of palmitic acid. Moreover, we found a significant pollen donor effect from the variety "SIAL 70" which was associated with decreased palmitic acid levels.

Low-cost assembly of a cacao crop genome is able to resolve complex heterozygous bubbles

Cacao (Theobroma cacao) is a tropical tree that produces the essential raw material for chocolate. Because yields have been stagnant, land use has expanded to provide for increasing chocolate demand. Assembled genomes of key parents could modernize breeding programs in the remote and under-resourced locations where cacao is grown. The MinION, a long read sequencer that runs off of a laptop computer, has the potential to facilitate the assembly of the complex genomes of high-yielding F₁ hybrids. Here, we validate the MinION's application to heterozygous crops by creating a de novo genome assembly of a key parent in breeding programs, the clone Pound 7. Our MinION-only assembly was 20% larger than the latest released cacao genome, with 10-fold greater contiguity, and the resolution of complex heterozygosity and repetitive elements. Polishing with Illumina short reads brought the predicted completeness of our assembly to similar levels to the previously released cacao genome assemblies. In contrast to previous cacao genome projects, our assembly required only a small scientific team and limited reagents. Our sequencing and assembly methods could easily be adopted by under-resourced breeding programs, speeding crop improvement in the developing world.

Association Between SSR Molecular Markers and Chemical and Sensory Traits of Cacao Samples Using Multiple Regression Analysis

The sensory properties of cacao beans are linked to the chemical composition of the seeds, and both characteristics are the partial results of its allelic composition. Therefore, it is useful to search for molecular markers associated with these traits. We perform multiple regression analysis to associate previously generated data of alleles generated with 12 SSR (of cultivated cacao trees) with data obtained from chemical and sensory characterization (of beans) of plants grown in the southern region from Mexico. When the association was significant, the mathematical models for predictive purposes were proposed. All phenotypic traits evaluated showed equations with setting values R² > 0.5. All chemical characters tested have a significant association with at least two alleles (P < 0.05). In addition, the fat content was associated with six molecular markers (mTcCIR03₂₀₉, mTcCIR12₁₈₈, mTcCIR19₂₈₆, mTcCIR07₁₅₀, mTcCIR19₃₁₀). The most common allele was mTcCIR12₁₈₈, which was associated with the contents of eicosanoic acid, moisture, fat and total polyphenols content. The mTcCIR28₃₆₂ allele is associated with sensory characters bitterness, musty odor, and roasted odor. These alleles could be useful as molecular markers of chemical and sensory characteristics of cacao samples.

Mapping of a Major QTL for Ceratocystis Wilt Disease in an F1 Population of Theobroma cacao

Cacao is an important crop, its beans are key raw materials for the chocolate and cosmetic industries. Ceratocystis wilt of cacao (CWC) caused by Ceratocystis cacaofunesta is a lethal disease for the crop. Therefore, the selection of resistant cacao varieties is one of the viable ways to minimize losses in cacao production. In this paper, we described the identification of a major QTL associated with CWC in an F1 mapping population from a cross between a resistant, "TSH 1188," and a susceptible genotype, "CCN 51." A set of 266 trees were genotyped using 3,526 single nucleotide polymorphic markers and then multiple QTL mapping analyses were performed. Two QTLs were identified on chromosomes IV and VI. The major QTL was located at 20 cM from the top position of chromosome VI, accounting for more than 60% of the phenotypic variation. The favorable allele T1, with haplotype GTT, came from the "TSH 1188" parent. It was evident that the haplotype combination T1C2 on chromosome VI was the most significant for resistance, since 93% of resistant trees had this haplotype. The major QTL converged to a genomic region of 739.4 kb that harbored nine candidate genes, including two major classes of resistance genes, which would make them the primary candidates involved in the resistance to CWC. The haplotypes detected are now used to improve the efficiency and precision of the selection of resistant trees in cacao breeding.

Prediction of Cacao (Theobroma cacao) Resistance to Moniliophthora spp. Diseases via Genome-Wide Association Analysis and Genomic Selection

Cacao (Theobroma cacao) is a globally important crop, and its yield is severely restricted by disease. Two of the most damaging diseases, witches' broom disease (WBD) and frosty pod rot disease (FPRD), are caused by a pair of related fungi: Moniliophthora perniciosa and Moniliophthora roreri, respectively. Resistant cultivars are the most effective long-term strategy to address Moniliophthora diseases, but efficiently generating resistant and productive new cultivars will require robust methods for screening germplasm before field testing. Marker-assisted selection (MAS) and genomic selection (GS) provide two potential avenues for predicting the performance of new genotypes, potentially increasing the selection gain per unit time. To test the effectiveness of these two approaches, we performed a genome-wide association study (GWAS) and GS on three related populations of cacao in Ecuador genotyped with a 15K single nucleotide polymorphism (SNP) microarray for three measures of WBD infection (vegetative broom, cushion broom, and chirimoya pod), one of FPRD (monilia pod) and two productivity traits (total fresh weight of pods and % healthy pods produced). GWAS yielded several SNPs associated with disease resistance in each population, but none were significantly correlated with the same trait in other populations. Genomic selection, using one population as a training set to estimate the phenotypes of the remaining two (composed of different families), varied among traits, from a mean prediction accuracy of 0.46 (vegetative broom) to 0.15 (monilia pod), and varied between training populations. Simulations demonstrated that selecting seedlings using GWAS markers alone generates no improvement over selecting at random, but that GS improves the selection process significantly. Our results suggest that the GWAS markers discovered here are not sufficiently predictive across diverse germplasm to be useful for MAS, but that using all markers in a GS framework holds substantial promise in accelerating disease-resistance in cacao.

Phenotypic Description of Theobroma cacao L. for Yield and Vigor Traits From 34 Hybrid Families in Costa Rica Based on the Genetic Basis of the Parental Population

A comprehensive understanding of the genetic basis of target traits in any crop is critical to design breeding strategies for the development and release of new improved varieties. In this study, 34 cacao families were evaluated for vigor and yield related traits over the course of 6 years in Costa Rica. Linear mixed models provided the variance components for the partitioning of additive and non-additive effects. Heritabilities of yield over time ranged from 0.085 to 0.576, from 0.127 to 0.399 for vigor, and 0.141 to 0.146 for disease resistance traits. Significant (p < 0.001) general combining abilities were observed for ICS-43 and LcTeen-37 with negative effect on average yield (-0.674, -0.690), respectively. Specific combining abilities for yield had significant (p < 0.001) positive effect from the cross GU-154-L x UF-273 Type 2 (0.703) and strong negative interaction between ICS-43 and LF-1 (-0.975). A weighted index was used to select the top performers while providing the corresponding genetic gains. At an 1% selection intensity, yield component gains ranged from 17.8 to 331.9%. Agronomic traits such as branch angle, trunk diameter and jorquette height had lower genetic gains and lower heritabilities. In addition, the parents in this study were genotyped with a 96-SNP marker off-typing set and a significant positive correlation of 0.39 (p = 0.019) was found between genetic distance and specific combining abilities for yield. Preliminary comparison of clonal parents vs. seedlings yield in the family with the highest SCA suggest for the first time presence of heterobeltiosis in cacao.

Genetic Parameters and the Impact of Off-Types for Theobroma cacao L. in a Breeding Program in Brazil

Breeding programs of cacao (Theobroma cacao L.) trees share the many challenges of breeding long-living perennial crops, and genetic progress is further constrained by both the limited understanding of the inheritance of complex traits and the prevalence of technical issues, such as mislabeled individuals (off-types). To better understand the genetic architecture of cacao, in this study, 13 years of phenotypic data collected from four progeny trials in Bahia, Brazil were analyzed jointly in a multisite analysis. Three separate analyses (multisite, single site with and without off-types) were performed to estimate genetic parameters from statistical models fitted on nine important agronomic traits (yield, seed index, pod index, % healthy pods, % pods infected with witches broom, % of pods other loss, vegetative brooms, diameter, and tree height). Genetic parameters were estimated along with variance components and heritabilities from the multisite analysis, and a trial was fingerprinted with low-density SNP markers to determine the impact of off-types on estimations. Heritabilities ranged from 0.37 to 0.64 for yield and its components and from 0.03 to 0.16 for disease resistance traits. A weighted index was used to make selections for clonal evaluation, and breeding values estimated for the parental selection and estimation of genetic gain. The impact of off-types to breeding progress in cacao was assessed for the first time. Even when present at <5% of the total population, off-types altered selections by 48%, and impacted heritability estimations for all nine of the traits analyzed, including a 41% difference in estimated heritability for yield. These results show that in a mixed model analysis, even a low level of pedigree error can significantly alter estimations of genetic parameters and selections in a breeding program.