Deciphering the Theobroma cacao self-incompatibility system: from genomics to diagnostic markers for self-compatibility

Genomic decoding of Theobroma grandiflorum (cupuassu) at chromosomal scale: evolutionary insights for horticultural innovation

BACKGROUND

Theobroma grandiflorum (Malvaceae), known as cupuassu, is a tree indigenous to the Amazon basin, valued for its large fruits and seed pulp, contributing notably to the Amazonian bioeconomy. The seed pulp is utilized in desserts and beverages, and its seed butter is used in cosmetics. Here, we present the sequenced telomere-to-telomere genome of cupuassu, disclosing its genomic structure, evolutionary features, and phylogenetic relationships within the Malvaceae family.

FINDINGS

The cupuassu genome spans 423 Mb, encodes 31,381 genes distributed in 10 chromosomes, and exhibits approximately 65% gene synteny with the Theobroma cacao genome, reflecting a conserved evolutionary history, albeit punctuated with unique genomic variations. The main changes are pronounced by bursts of long-terminal repeat retrotransposons at postspecies divergence, retrocopied and singleton genes, and gene families displaying distinctive patterns of expansion and contraction. Furthermore, positively selected genes are evident, particularly among retained and dispersed tandem and proximal duplicated genes associated with general fruit and seed traits and defense mechanisms, supporting the hypothesis of potential episodes of subfunctionalization and neofunctionalization following duplication, as well as impact from distinct domestication process. These genomic variations may underpin the differences observed in fruit and seed morphology, ripening, and disease resistance between cupuassu and the other Malvaceae species.

CONCLUSIONS

The cupuassu genome offers a foundational resource for both breeding improvement and conservation biology, yielding insights into the evolution and diversity within the genus Theobroma.

3D pollination biology using micro-computed tomography and geometric morphometrics in Theobroma cacao

Premise

Imaging technologies that capture three-dimensional (3D) variation in floral morphology at micro- and nano-resolutions are increasingly accessible. In herkogamous flowers, such as those of Theobroma cacao, structural barriers between anthers and stigmas represent bottlenecks that restrict pollinator size and access to reproductive organs. To study the unresolved pollination biology of cacao, we present a novel application of micro-computed tomography (micro-CT) using floral dimensions to quantify pollinator functional size limits.

Methods

We generated micro-CT data sets from field-collected flowers and museum specimens of potential pollinators. To compare floral variation, we used 3D Slicer to place landmarks on the surface models and performed a geometric morphometric (GMM) analysis using geomorph R. We identified the petal side door (an opening between the petal hoods and filament) as the main bottleneck for pollinator access. We compared its mean dimensions with proposed pollinators to identify viable candidates.

Results

We identified three levels of likelihood for putative pollinators based on the number of morphological (body) dimensions that fit through the petal side door. We also found floral reward microstructures whose presence and location were previously unclear.

Discussion

Using micro-CT and GMM to study the 3D pollination biology of cacao provides new evidence for predicting unknown pollinators. Incorporating geometry and floral rewards will strengthen plant-pollinator trait matching models for cacao and other species.

Heterosis patterns and sources of self-compatibility, cross-compatibility and key nut traits within single and double hybrid crosses of kola [Cola nitida (Vent) Schott and Endl.]

Sexual incompatibility among kola genotypes accounted for over 50% yield loss. Compatible and high yielding varieties are in demand to develop commercial orchards. The objective of this study was to assess self-compatibility and cross-compatibility of kola (C. nitida) genotypes within self, single and double hybrid crosses and to determine heterosis pattern in the resulting hybrids for sexual compatibility and key nut yield and quality traits. Crosses among kola genotypes from three field gene banks (JX1, GX1, MX2) and one advanced germplasm (Bunso progeny) in Ghana were evaluated along their parents for sexual compatibility, nut yield and nut quality. Data were collected on pod set, pseudo-pod set, pod weight, number of nuts per pod, nut weight, brix, potential alcohol and nut firmness. Significant (P < 0.001) differential pod set was observed within Bunso progeny, JX1, GX1 and MX2 crosses; while pseudo-pod set differed only within JX1 and MX2 crosses (P < 0.001). Very large prevalence of mid-parent, heterobeltiosis, and economic heterosis was observed for sexual compatibility, outturn and brix for the single and double hybrid crosses. Heterosis was prominent among the double hybrid crosses as compared to the single hybrid crosses suggesting that recurrent selection of compatible varieties from advanced generations could result in genetic gain in kola improvement. The top five crosses with best heterosis for sexual compatibility and an appreciable positive heterosis for outturn and brix were B1/11 × B1/71 × B1/157 × B1/149, B1/11 × B1/71 × B1/296 × B1/177, GX1/46 × GX1/33 × B1/212 × B1/236, JX1/90 × JX1/51 and JX1/51 × JX1/36. These materials could serve as sources of beneficial alleles for improving Ghanaian kola hybrids and populations for yield and sexual compatibility.

Genome-wide SNP genotyping as a simple and practical tool to accelerate the development of inbred lines in outbred tree species: An example in cacao (Theobroma cacao L.)

Cacao is a globally important crop with a long history of domestication and selective breeding. Despite the increased use of elite clones by cacao farmers, worldwide plantations are established mainly using hybrid progeny material derived from heterozygous parents, therefore displaying high tree-to-tree variability. The deliberate development of hybrids from advanced inbred lines produced by successive generations of self-pollination has not yet been fully considered in cacao breeding. This is largely due to the self-incompatibility of the species, the long generation cycles (3–5 years) and the extensive trial areas needed to accomplish the endeavor. We propose a simple and accessible approach to develop inbred lines based on accelerating the buildup of homozygosity based on regular selfing assisted by genome-wide SNP genotyping. In this study we genotyped 90 clones from the Brazilian CEPEC´s germplasm collection and 49 inbred offspring of six S₁ or S₂ cacao families derived from self-pollinating clones CCN-51, PS-13.19, TSH-1188 and SIAL-169. A set of 3,380 SNPs distributed across the cacao genome were interrogated on the EMBRAPA multi-species 65k Infinium chip. The 90 cacao clones showed considerable variation in genome-wide SNP homozygosity (mean 0.727± 0.182) and 19 of them with homozygosity ≥90%. By assessing the increase in homozygosity across two generations of self-pollinations, SNP data revealed the wide variability in homozygosity within and between S₁ and S₂ families. Even in small families (<10 sibs), individuals were identified with up to ~1.5 standard deviations above the family mean homozygosity. From baseline homozygosities of 0.476 and 0.454, offspring with homozygosities of 0.862 and 0.879 were recovered for clones TSH-1188 and CCN-51 respectively, in only two generations of selfing (81–93% increase). SNP marker assisted monitoring and selection of inbred individuals can be a practical tool to optimize and accelerate the development of inbred lines of outbred tree species. This approach will allow a faster and more accurate exploitation of hybrid breeding strategies in cacao improvement programs and potentially in other perennial fruit and forest trees.

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.

Inflorescence Transcriptome Sequencing and Development of New EST-SSR Markers in Common Buckwheat (Fagopyrum esculentum)

Common buckwheat (Fagopyrum esculentum M.) is known for its adaptability, good nutrition, and medicinal and health care value. However, genetic studies of buckwheat have been hindered by limited genomic resources and genetic markers. In this study, Illumina HiSeq 4000 high-throughput sequencing technology was used to sequence the transcriptome of green-flower common buckwheat (Gr) with coarse pedicels and white-flower Ukrainian daliqiao (UD) with fine pedicels. A total of 118,448 unigenes were obtained, with an average length of 1248 bp and an N50 of 1850 bp. A total of 39,432 differentially expressed genes (DEGs) were identified, and the DEGs of the porphyrins and chlorophyll metabolic pathway had significantly upregulated expression in Gr. Then, a total of 17,579 sequences containing SSR loci were detected, and 20,756 EST-SSR loci were found. The distribution frequency of EST-SSR in the transcriptome was 17.52%, and the average distribution density was 8.21 kb. A total of 224 pairs of primers were randomly selected for synthesis; 35 varieties of common buckwheat and 13 varieties of Tartary buckwheat were verified through these primers. The clustering results well verified the previous conclusion that common buckwheat and Tartary buckwheat had a distant genetic relationship. The EST-SSR markers identified and developed in this study will be helpful to enrich the transcriptome information and marker-assisted selection breeding of buckwheat.

Pollen-Pistil Interactions as Reproductive Barriers

Pollen-pistil interactions serve as important prezygotic reproductive barriers that play a critical role in mate selection in plants. Here, we highlight recent progress toward understanding the molecular basis of pollen-pistil interactions as reproductive isolating barriers. These barriers can be active systems of pollen rejection, or they can result from a mismatch of required male and female factors. In some cases, the barriers are mechanistically linked to self-incompatibility systems, while others represent completely independent processes. Pollen-pistil reproductive barriers can act as soon as pollen is deposited on a stigma, where penetration of heterospecific pollen tubes is blocked by the stigma papillae. As pollen tubes extend, the female transmitting tissue can selectively limit growth by producing cell wall-modifying enzymes and cytotoxins that interact with the growing pollen tube. At ovules, differential pollen tube attraction and inhibition of sperm cell release can act as barriers to heterospecific pollen tubes.

Do s genes or deleterious recessives control late-acting self-incompatibility in Handroanthus heptaphyllus (Bignoniaceae)? A diallel study with four full-sib progeny arrays

BACKGROUND AND AIMS

Genetically controlled self-incompatibility (SI) mechanisms constrain selfing and thus have contributed to the evolutionary diversity of flowering plants. In homomorphic gametophytic SI (GSI) and homomorphic sporophytic SI (SSI), genetic control is usually by the single multi-allelic locus S. Both GSI and SSI prevent self pollen tubes reaching the ovary and so are pre-zygotic in action. In contrast, in taxa with late-acting self-incompatibility (LSI), rejection is often post-zygotic, since self pollen tubes grow to the ovary, where fertilization may occur prior to floral abscission. Alternatively, lack of self fruit set could be due to early-acting inbreeding depression (EID). The aim of our study was to investigate mechanisms underlying the lack of selfed fruit set in Handroanthus heptaphyllus in order to assess the likelihood of LSI versus EID.

METHODS

We employed four full-sib diallels to study the genetic control of LSI in H. heptaphyllus using a precociously flowering variant. We also used fluorescence microscopy to study the incidence of ovule penetration by pollen tubes in pistils that abscised following pollination or initiated fruits.

KEY RESULTS

All diallels showed reciprocally cross-incompatible full sibs (RCIs), reciprocally cross-compatible full sibs (RCCs) and non-reciprocally compatible full sibs (NRCs) in almost equal proportions. There was no significant difference between the incidences of ovule penetrations in abscised pistils following self- and cross-incompatible pollinations, but those in successful cross-pollinations were around 2-fold greater.

CONCLUSIONS

A genetic model postulating a single S locus with four S alleles, one of which, in the maternal parent, is dominant to the other three, will produce RCI, RCC and NRC full sib situations each at 33 %, consistent with our diallel results. We favour this simple genetic control over an EID explanation since none of our pollinations, successful or unsuccessful, resulted in partial embryo development, as would be expected under a whole-genome EID effect.

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.

In-depth Understanding of Camellia oleifera Self-incompatibility by Comparative Transcriptome, Proteome and Metabolome

Oil-tea tree (Camellia oleifera) is the most important edible oil tree species in China with late-acting self-incompatibility (LSI) properties. The mechanism of LSI is uncertain, which seriously hinders the research on its genetic characteristics, construction of genetic map, selection of cross breeding parents and cultivar arrangement. To gain insights into the LSI mechanism, we performed cytological, transcriptomic, proteomic and metabolomic studies on self- and cross-pollinated pistils. The studies identified 166,591 transcripts, 6851 proteins and 6455 metabolites. Transcriptomic analysis revealed 1197 differentially expressed transcripts between self- and cross-pollinated pistils and 47 programmed cell death (PCD)-control transcripts. Trend analysis by Pearson correlation categorized nine trend graphs linked to 226 differentially expressed proteins and 38 differentially expressed metabolites. Functional enrichment analysis revealed that the LSI was closely associated with PCD-related genes, mitogen-activated protein kinase (MAPK) signaling pathway, plant hormone signal transduction, ATP-binding cassette (ABC) transporters and ubiquitin-mediated proteolysis. These particular trends in transcripts, proteins and metabolites suggested the involvement of PCD in LSI. The results provide a solid genetic foundation for elucidating the regulatory network of PCD-mediated self-incompatibility in C. oleifera.

Clovamide, a Hydroxycinnamic Acid Amide, Is a Resistance Factor Against Phytophthora spp. in Theobroma cacao

The hydroxycinnamic acid amides (HCAAs) are a diverse group of plant-specialized phenylpropanoid metabolites distributed widely in the plant kingdom and are known to be involved in tolerance to abiotic and biotic stress. The HCAA clovamide is reported in a small number of distantly related species. To explore the contribution of specialized metabolites to disease resistance in cacao (Theobroma cacao L., chocolate tree), we performed untargeted metabolomics using liquid chromatography - tandem mass spectrometry (LC-MS/MS) and compared the basal metabolite profiles in leaves of two cacao genotypes with contrasting levels of susceptibility to Phytophthora spp. Leaves of the tolerant genotype 'Scavina 6' ('Sca6') were found to accumulate dramatically higher levels of clovamide and several other HCAAs compared to the susceptible 'Imperial College Selection 1' ('ICS1'). Clovamide was the most abundant metabolite in 'Sca6' leaf extracts based on MS signal, and was up to 58-fold higher in 'Sca6' than in 'ICS1'. In vitro assays demonstrated that clovamide inhibits growth of three pathogens of cacao in the genus Phytophthora, is a substrate for cacao polyphenol oxidase, and is a contributor to enzymatic browning. Furthermore, clovamide inhibited proteinase and pectinase in vitro, activities associated with defense in plant-pathogen interactions. Fruit epidermal peels from both genotypes contained substantial amounts of clovamide, but two sulfated HCAAs were present at high abundance exclusively in 'Sca6' suggesting a potential functional role of these compounds. The potential to breed cacao with increased HCAAs for improved agricultural performance is discussed.

Self-(In)compatibility Systems: Target Traits for Crop-Production, Plant Breeding, and Biotechnology

Self-incompatibility (SI) mechanisms prevent self-fertilization in flowering plants based on specific discrimination between self- and non-self pollen. Since this trait promotes outcrossing and avoids inbreeding it is a widespread mechanism of controlling sexual plant reproduction. Growers and breeders have effectively exploited SI as a tool for manipulating domesticated crops for thousands of years. However, only within the past thirty years have studies begun to elucidate the underlying molecular features of SI. The specific S-determinants and some modifier factors controlling SI have been identified in the sporophytic system exhibited by Brassica species and in the two very distinct gametophytic systems present in Papaveraceae on one side and in Solanaceae, Rosaceae, and Plantaginaceae on the other. Molecular level studies have enabled SI to SC transitions (and vice versa) to be intentionally manipulated using marker assisted breeding and targeted approaches based on transgene integration, silencing, and more recently CRISPR knock-out of SI-related factors. These scientific advances have, in turn, provided a solid basis to implement new crop production and plant breeding practices. Applications of self-(in)compatibility include widely differing objectives such as crop yield and quality improvement, marker-assisted breeding through SI genotyping, and development of hybrids for overcoming intra- and interspecific reproductive barriers. Here, we review scientific progress as well as patented applications of SI, and also highlight future prospects including further elucidation of SI systems, deepening our understanding of SI-environment relationships, and new perspectives on plant self/non-self recognition.

Genetic insights into the modification of the pre-fertilization mechanisms during plant domestication

Plant domestication is the process of adapting plants to human use by selecting specific traits. The selection process often involves the modification of some components of the plant reproductive mechanisms. Allelic variants of genes associated with flowering time, vernalization, and the circadian clock are responsible for the adaptation of crops, such as rice, maize, barley, wheat, and tomato, to non-native latitudes. Modifications in the plant architecture and branching have been selected for higher yields and easier harvests. These phenotypes are often produced by alterations in the regulation of the transition of shoot apical meristems to inflorescences, and then to floral meristems. Floral homeotic mutants are responsible for popular double-flower phenotypes in Japanese cherries, roses, camellias, and lilies. The rise of peloric flowers in ornamentals such as snapdragon and florists' gloxinia is associated with non-functional alleles that control the relative expansion of lateral and ventral petals. Mechanisms to force outcrossing such as self-incompatibility have been removed in some tree crops cultivars such as almonds and peaches. In this review, we revisit some of these important concepts from the plant domestication perspective, focusing on four topics related to the pre-fertilization mechanisms: flowering time, inflorescence architecture, flower development, and pre-fertilization self-incompatibility mechanisms.

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.