Stefan Royaert

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

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

Optimization of a SNP assay for genotyping Theobroma cacao under field conditions

The tropical tree crop Theobroma cacao L. is grown commercially for its beans, which are used in the production of cocoa butter and chocolate. Although the upper Amazon region is the center of origin for cacao, 70% of the world’s supply of cacao beans currently comes from small farms in West Africa. While cacao breeding programs in producer nations are the source of improved planting material, modern marker-based breeding is difficult to perform due to the lack of genotyping facilities in these countries. While DNA extraction can be routinely performed, the equipment needed to analyze simple sequence repeats (SSRs) is seldom available, forcing the outsourcing of genotyping to foreign laboratories and delaying the breeding process. We describe a 5′ nuclease (TaqMan)-based single nucleotide polymorphism (SNP) assay for genotyping cacao plants under conditions similar to those found in most cacao-producing areas. The assay was tested under field conditions by planting open pollinated seeds of seven pods from four different maternal plants. The resulting 171 seedlings were successfully genotyped with 18 SNP markers representing 12 loci. The ability to use temperature-stable reagents and rapid DNA extraction methods is also explored. Additionally, by examining the seedling genotypes for the SNP markers and 14 additional SSR markers, we investigated whether seeds in a pod are the result of single or multiple pollination events. This simple, effective method of genotyping cacao seedlings in the field should allow for more efficient resource management of seed gardens and is currently being implemented in Ghana.

Assessing microsatellite linkage disequilibrium in wild, cultivated, and mapping populations of Theobroma cacao L. and its impact on association mapping

Linkage disequilibrium (LD) measured over the genomes of a species can provide important indications for how future association analyses should proceed. This information can be advantageous especially for slow-growing, perennial crops such as Theobroma cacao, where experimental crosses are inherently time-consuming and logistically expensive. While LD has been evaluated in cacao, previous work has been focused on relatively narrow genetic bases. We use microsatellite marker data collected from a uniquely diverse sample of individuals broadly covering both wild and cultivated varieties to gauge the LD present in the different cacao diversity groups and populations. We find that genome-wide LD decays far more rapidly in the wild and primitive diversity groups of cacao as compared to those representing cultivated varieties. The impact that such differences can have on association analyses is demonstrated using phenotypic data on pod color and genotypic data from two cacao populations with contrasting patterns of LD decay. Our results indicate that the more rapid LD decay in wild and primitive germplasm can lead to higher-resolution mapping intervals when compared to results from cultivated germplasm. Through simulations, we demonstrate how future association mapping analyses, comprising of cacao samples with a wild or primitive background, will likely exhibit lower LD and would be more suitable for fine-scale association mapping analyses. As many traits targeted by cacao breeders are found exclusively in wild and primitive germplasm, association mapping in wild cacao populations holds significant promise for cacao improvement through marker-assisted breeding and emphasize the need to further explore the natural diversity of Amazonian cacao.

iXora: exact haplotype inferencing and trait association.

BackgroundWe address the task of extracting accurate haplotypes from genotype data of individuals of large F1 populations for mapping studies. While methods for inferring parental haplotype assignments on large F1 populations exist in theory, these approaches do not work in practice at high levels of accuracy.ResultsWe have designed iXora (Identifying crossovers and recombining alleles), a robust method for extracting reliable haplotypes of a mapping population, as well as parental haplotypes, that runs in linear time. Each allele in the progeny is assigned not just to a parent, but more precisely to a haplotype inherited from the parent. iXora shows an improvement of at least 15% in accuracy over similar systems in literature. Furthermore, iXora provides an easy-to-use, comprehensive environment for association studies and hypothesis checking in populations of related individuals.ConclusionsiXora provides detailed resolution in parental inheritance, along with the capability of handling very large populations, which allows for accurate haplotype extraction and trait association. iXora is available for non-commercial use from http://researcher.ibm.com/project/3430.

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.

Population genomic analyses of the chocolate tree, Theobroma cacao L., provide insights into its domestication process

Domestication has had a strong impact on the development of modern societies. We sequenced 200 genomes of the chocolate plant Theobroma cacao L. to show for the first time to our knowledge that a single population, the Criollo population, underwent strong domestication ~3600 years ago (95% CI: 2481–13,806 years ago). We also show that during the process of domestication, there was strong selection for genes involved in the metabolism of the colored protectants anthocyanins and the stimulant theobromine, as well as disease resistance genes. Our analyses show that domesticated populations of T. cacao (Criollo) maintain a higher proportion of high-frequency deleterious mutations. We also show for the first time the negative consequences of the increased accumulation of deleterious mutations during domestication on the fitness of individuals (significant reduction in kilograms of beans per hectare per year as Criollo ancestry increases, as estimated from a GLM, P = 0.000425).Omar Cornejo et al. report a genomic analysis of 200 cacao plants (Theobroma cacao L.) representing more than 10 genetically distinct populations. They identify metabolic and disease resistance genes as contributing to the domestication of cacao and show that domesticated populations maintain a high proportion of deleterious mutations.

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.

Genomic insights into the domestication of the chocolate tree, Theobroma cacao L.

Domestication has had a strong impact on the development of modern societies. We sequenced 200 genomes of the chocolate plant Theobroma cacao L. to show for the first time that a single population underwent strong domestication approximately 3,600 years (95% CI: 2481 – 10,903 years ago) ago, the Criollo population. We also show that during the process of domestication, there was strong selection for genes involved in the metabolism of the colored protectants anthocyanins and the stimulant theobromine, as well as disease resistance genes. Our analyses show that domesticated populations of T. cacao (Criollo) maintain a higher proportion of high frequency deleterious mutations. We also show for the first time the negative consequences the increase accumulation of deleterious mutations during domestication on the fitness of individuals (significant negative correlation between Criollo ancestry and Kg of beans per hectare per year, P = 0.000425).