Single-base resolution methylomes of somatic embryogenesis in Theobroma cacao L. reveal epigenome modifications associated with somatic embryo abnormalities.

Propagation by somatic embryogenesis in Theobroma cacao has some issues to be solved, as many morphologically abnormal somatic embryos that do not germinate into plants are frequently observed, thus hampering plant production on a commercial scale. For the first time the methylome landscape of T. cacao somatic embryogenesis was examined, using whole-genome bisulfite sequencing technique, with the aim to understand the epigenetic basis of somatic embryo abnormalities. We identified 873 differentially methylated genes (DMGs) in the CpG context between zygotic embryos, normal and abnormal somatic embryos, with important roles in development, programmed cell death, oxidative stress, and hypoxia induction, which can help to explain the morphological abnormalities of somatic embryos. We also identified the role of ethylene and its precursor 1-aminocyclopropane-1-carboxylate in several biological processes, such as hypoxia induction, cell differentiation and cell polarity, that could be associated to the development of abnormal somatic embryos. The biological processes and the hypothesis of ethylene and its precursor involvement in the somatic embryo abnormalities in cacao are discussed.

QTL mapping and identification of SNP-haplotypes affecting yield components of Theobroma cacao L.

Cacao is a crop of global relevance that faces constant demands for improved bean yield. However, little is known about the genomic regions controlling the crop yield and genes involved in cacao bean filling. Hence, to identify the quantitative trait loci (QTL) associated with cacao yield and bean filling, we performed a QTL mapping in a segregating mapping population comprising 459 trees of a cross between 'TSH 1188' and 'CCN 51'. All variables showed considerable phenotypic variation and had moderate to high heritability values. We identified 24 QTLs using a genetic linkage map that contains 3526 single nucleotide polymorphism (SNP) markers. Haplotype analysis at the significant QTL region on chromosome IV pointed to the alleles from the maternal parent, 'TSH 1188', as the ones that affect the cacao yield components the most. The recombination events identified within these QTL regions allowed us to identify candidate genes that may take part in the different steps of pod growth and bean filling. Such candidate genes seem to play a significant role in the source-to-sink transport of sugars and amino acids, and lipid metabolism, such as fatty acid production. The SNP markers mapped in our study are now being used to select potential high-yielding cacao varieties through marker-assisted selection in our existing cacao-breeding experiments.

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.

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

Climate change could threaten cocoa production: Effects of 2015-16 El Niño-related drought on cocoa agroforests in Bahia, Brazil.

Climate models predict a possible increase in the frequency of strong climate events such as El Niño-Southern Oscillation (ENSO), which in parts of the tropics are the cause of exceptional droughts, these threaten global food production. Agroforestry systems are often suggested as promising diversification options to increase farmers' resilience to extreme climatic events. In the Northeastern state of Bahia, where most Brazilian cocoa is grown in wildlife-friendly agroforests, ENSOs cause severe droughts which negatively affect forest and agriculture. Cocoa (Theobroma cacao) is described as being sensitive to drought but there are no field-studies of the effect of ENSO-related drought on adult cocoa trees in the America's; there is one study of an experimentally-imposed drought in Indonesia which resulted in 10 to 46% yield loss. In our study, in randomly chosen farms in Bahia, Brazil, we measured the effect of the 2015-16 severe ENSO, which caused an unprecedented drought in cocoa agroforests. We show that drought caused high cocoa tree mortality (15%) and severely decreased cocoa yield (89%); the drought also increased infection rate of the chronic fungal disease witches' broom (Moniliophthora perniciosa). Ours findings showed that Brazilian cocoa agroforests are at risk and that increasing frequency of strong droughts are likely to cause decreased cocoa yields in the coming decades. Furthermore, because cocoa, like many crops, is grown somewhat beyond its climatic limits, it and other crops could be the 'canaries in the coalmine' warning of forthcoming major drought effects on semi-natural and natural vegetation.

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.

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.

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.

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.

Analysis of the Petunia TM6 MADS box gene reveals functional divergence within the DEF/AP3 lineage.

Antirrhinum majus DEFICIENS (DEF) and Arabidopsis thaliana APETALA3 (AP3) MADS box proteins are required to specify petal and stamen identity. Sampling of DEF/AP3 homologs revealed two types of DEF/AP3 proteins, euAP3 and TOMATO MADS BOX GENE6 (TM6), within core eudicots, and we show functional divergence in Petunia hybrida euAP3 and TM6 proteins. Petunia DEF (also known as GREEN PETALS [GP]) is expressed mainly in whorls 2 and 3, and its expression pattern remains unchanged in a blind (bl) mutant background, in which the cadastral C-repression function in the perianth is impaired. Petunia TM6 functions as a B-class organ identity protein only in the determination of stamen identity. Atypically, Petunia TM6 is regulated like a C-class rather than a B-class gene, is expressed mainly in whorls 3 and 4, and is repressed by BL in the perianth, thereby preventing involvement in petal development. A promoter comparison between DEF and TM6 indicates an important change in regulatory elements during or after the duplication that resulted in euAP3- and TM6-type genes. Surprisingly, although TM6 normally is not involved in petal development, 35S-driven TM6 expression can restore petal development in a def (gp) mutant background. Finally, we isolated both euAP3 and TM6 genes from seven solanaceous species, suggesting that a dual euAP3/TM6 B-function system might be the rule in the Solanaceae.

The duplicated B-class heterodimer model: whorl-specific effects and complex genetic interactions in Petunia hybrida flower development.

In both Antirrhinum (Antirrhinum majus) and Arabidopsis (Arabidopsis thaliana), the floral B-function, which specifies petal and stamen development, is embedded in a heterodimer consisting of one DEFICIENS (DEF)/APETALA3 (AP3)-like and one GLOBOSA (GLO)/PISTILLATA (PI)-like MADS box protein. Here, we demonstrate that gene duplications in both the DEF/AP3 and GLO/PI lineages in Petunia hybrida (petunia) have led to a functional diversification of their respective members, which is reflected by partner specificity and whorl-specific functions among these proteins. Previously, it has been shown that mutations in PhDEF (formerly known as GREEN PETALS) only affect petal development. We have isolated insertion alleles for PhGLO1 (FLORAL BINDING PROTEIN1) and PhGLO2 (PETUNIA MADS BOX GENE2) and demonstrate unique and redundant properties of PhDEF, PhGLO1, and PhGLO2. Besides a full homeotic conversion of petals to sepals and of stamens to carpels as observed in phglo1 phglo2 and phdef phglo2 flowers, we found that gene dosage effects for several mutant combinations cause qualitative and quantitative changes in whorl 2 and 3 meristem fate, and we show that the PHDEF/PHGLO1 heterodimer controls the fusion of the stamen filaments with the petal tube. Nevertheless, when the activity of PhDEF, PhGLO1, and PhGLO2 are considered jointly, they basically appear to function as DEF/GLO does in Antirrhinum and to a lesser extent as AP3/PI in Arabidopsis. By contrast, our data suggest that the function of the fourth B-class MADS box member, the paleoAP3-type PETUNIA HYBRIDA TM6 (PhTM6) gene, differs significantly from the known euAP3-type DEF/AP3-like proteins; PhTM6 is mainly expressed in the developing stamens and ovary of wild-type flowers, whereas its expression level is upregulated in whorls 1 and 2 of an A-function floral mutant; PhTM6 is most likely not involved in petal development. The latter is consistent with the hypothesis that the evolutionary origin of the higher eudicot petal structure coincided with the appearance of the euAP3-type MADS box genes.