Molecular Surveillance, Prevalence, and Distribution of Cacao Infecting Badnavirus Species in Côte d'Ivoire and Ghana.

The cacao swollen shoot disease (CSSD) caused by a complex of badnavirus species presents a major challenge for cacao production in West Africa, especially Ghana and Côte d'Ivoire. In this study, CSSD species detection efficiency, diversity, and geographic distribution patterns in cacao plantations in Ghana and Côte d'Ivoire were investigated through field surveillance, PCR detection assays, sequencing of positive amplicons, and phylogeographic clustering. Cumulatively, the detection efficiency of the tested CSSD primer sets that were targeting the movement protein domain of the virus ranged from 0.15% (CSSD-3 primer) to 66.91% (CSSD-1 primer) on all the symptomatic cacao leaf samples assessed. The identified CSSD species differed phylogenetically and overlapped in distribution, with the cacao swollen shoot Togo B virus (CSSTBV) (n = 588 sequences) being the most prevalent and widely distributed compared to the other CSSD species that were encountered in both countries. Geographically, the cacao swollen shoot CE virus (CSSCEV) species (n = 124 sequences) that was identified was largely restricted to the bordering regions of Ghana and Côte d'Ivoire. These results provide updated knowledge of the geographic distribution of the key CSSD species and their diagnostic efficiency and, thus, provide guidance in identifying locations for structured testing of cacao germplasm and optimal diagnostics for the predominant CSSD species in Ghana and Côte d'Ivoire.

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).

A Larger Chocolate Chip-Development of a 15K Theobroma cacao L. SNP Array to Create High-Density Linkage Maps.

Cacao (Theobroma cacao L.) is an important cash crop in tropical regions around the world and has a rich agronomic history in South America. As a key component in the cosmetic and confectionary industries, millions of people worldwide use products made from cacao, ranging from shampoo to chocolate. An Illumina Infinity II array was created using 13,530 SNPs identified within a small diversity panel of cacao. Of these SNPs, 12,643 derive from variation within annotated cacao genes. The genotypes of 3,072 trees were obtained, including two mapping populations from Ecuador. High-density linkage maps for these two populations were generated and compared to the cacao genome assembly. Phenotypic data from these populations were combined with the linkage maps to identify the QTLs for yield and disease resistance.

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.

Identification of candidate genes involved in Witches' broom disease resistance in a segregating mapping population of Theobroma cacao L. in Brazil.

BACKGROUND: Witches' broom disease (WBD) caused by the fungus Moniliophthora perniciosa is responsible for considerable economic losses for cacao producers. One of the ways to combat WBD is to plant resistant cultivars. Resistance may be governed by a few genetic factors, mainly found in wild germplasm. RESULTS: We developed a dense genetic linkage map with a length of 852.8 cM that contains 3,526 SNPs and is based on the MP01 mapping population, which counts 459 trees from a cross between the resistant 'TSH 1188' and the tolerant 'CCN 51' at the Mars Center for Cocoa Science in Barro Preto, Bahia, Brazil. Seven quantitative trait loci (QTL) that are associated with WBD were identified on five different chromosomes using a multi-trait QTL analysis for outbreeders. Phasing of the haplotypes at the major QTL region on chromosome IX on a diversity panel of genotypes clearly indicates that the major resistance locus comes from a well-known source of WBD resistance, the clone 'SCAVINA 6'. Various potential candidate genes identified within all QTL may be involved in different steps leading to disease resistance. Preliminary expression data indicate that at least three of these candidate genes may play a role during the first 12 h after infection, with clear differences between 'CCN 51' and 'TSH 1188'. CONCLUSIONS: We combined the information from a large mapping population with very distinct parents that segregate for WBD, a dense set of mapped markers, rigorous phenotyping capabilities and the availability of a sequenced genome to identify several genomic regions that are involved in WBD resistance. We also identified a novel source of resistance that most likely comes from the 'CCN 51' parent. Thanks to the large population size of the MP01 population, we were able to pick up QTL and markers with relatively small effects that can contribute to the creation and selection of more tolerant/resistant plant material.

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 world's 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.

The genome sequence of the most widely cultivated cacao type and its use to identify candidate genes regulating pod color.

BACKGROUND: Theobroma cacao L. cultivar Matina 1-6 belongs to the most cultivated cacao type. The availability of its genome sequence and methods for identifying genes responsible for important cacao traits will aid cacao researchers and breeders. RESULTS: We describe the sequencing and assembly of the genome of Theobroma cacao L. cultivar Matina 1-6. The genome of the Matina 1-6 cultivar is 445 Mbp, which is significantly larger than a sequenced Criollo cultivar, and more typical of other cultivars. The chromosome-scale assembly, version 1.1, contains 711 scaffolds covering 346.0 Mbp, with a contig N50 of 84.4 kbp, a scaffold N50 of 34.4 Mbp, and an evidence-based gene set of 29,408 loci. Version 1.1 has 10x the scaffold N50 and 4x the contig N50 as Criollo, and includes 111 Mb more anchored sequence. The version 1.1 assembly has 4.4% gap sequence, while Criollo has 10.9%. Through a combination of haplotype, association mapping and gene expression analyses, we leverage this robust reference genome to identify a promising candidate gene responsible for pod color variation. We demonstrate that green/red pod color in cacao is likely regulated by the R2R3 MYB transcription factor TcMYB113, homologs of which determine pigmentation in Rosaceae, Solanaceae, and Brassicaceae. One SNP within the target site for a highly conserved trans-acting siRNA in dicots, found within TcMYB113, seems to affect transcript levels of this gene and therefore pod color variation. CONCLUSIONS: We report a high-quality sequence and annotation of Theobroma cacao L. and demonstrate its utility in identifying candidate genes regulating traits.

iXora: exact haplotype inferencing and trait association.

BACKGROUND: We 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. RESULTS: We 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. CONCLUSIONS: iXora 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.

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.

Reduction of protease inhibitor activity by expression of a mutant Bowman-Birk gene in soybean seed.

A mutant Bowman-Birk gene was created that encoded an inactive high-sulfur product. It was used to transform soybean line Asgrow 3237. Transformants bearing the mutant gene were identified by GUS expression, PCR analysis, and Southern analysis. The amount of steady state mRNA from the mutant gene in the transformed plants showed that the gene was highly expressed, but the amount of message from the unmodified Bowman-Birk gene did not change detectably. Proteins synthesized at the direction of the mutant Bowman-Birk gene accumulated in seeds of the transformed plants, and there was a marked decrease in the ability of extracts prepared from these seeds to inhibit trypsin and chymotrypsin despite the presence of Kunitz trypsin inhibitor. The more prevalent mRNA from the mutant gene was considered to out-compete message from the native genes to decrease the amount of active Bowman-Birk inhibitor.

Two subtilisin-like proteases from soybean.

Two subtilisin-like proteases (SLP) were identified in soybean (Glycine max [L.] Merr.). The first, SLP-1, was localized in seed coats early in seed development, but became undetectable with anti-SLP-1 antibodies as seed fill progressed. A partial purification of SLP-1 was achieved using a two step chromatographic procedure. NH2-terminal sequence analysis of the partially purified enzyme permitted primers to be designed that were used to amplify cDNA encoding SLP-1. A genomic clone encoding SLP-1 was also obtained. Characterization of the cDNA and partially purified SLP-1 revealed the initial translation product was an 82 694 MW precursor. After removal of a signal peptide, the mature protein was formed by removal of an NH2-terminal propeptide. A COOH-terminal peptide also appeared to be removed from some of the protease molecules. DNA blot analysis suggested that at least one additional SLP gene was present in soybean. The second gene, SLP-2, was subsequently cloned and characterized. Although the coding regions for SLP-1 and SLP-2 were homologous, their promoters were quite divergent. RT-PCR revealed that SLP-2 message was found in the mature plant and in cotyledons of germinating seeds. Although SLP-2 mRNA could be identified in developing seeds, the message was at least an order of magnitude less abundant than that for SLP-1, and it was mis-spliced such that a chain termination event would preclude obtaining a product. As with SLPs from other organisms, the functions of the soybean proteases are unknown. However, SLP-1 is one of only a few proteins from soybean seed coats that have been described.