Mapping quantitative trait loci for bean traits and ovule number in Theobroma cacao L

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.

Molecular Network for Regulation of Ovule Number in Plants

In seed-bearing plants, the ovule ("small egg") is the organ within the gynoecium that develops into a seed after fertilization. The gynoecium located in the inner compartment of the flower turns into a fruit. The number of ovules in the ovary determines the upper limit or the potential of seed number per fruit in plants, greatly affecting the final seed yield. Ovule number is an important adaptive characteristic for plant evolution and an agronomic trait for crop improvement. Therefore, understanding the mechanism and pathways of ovule number regulation becomes a significant research aspect in plant science. This review summarizes the ovule number regulators and their regulatory mechanisms and pathways. Specially, an integrated molecular network for ovule number regulation is constructed, in which phytohormones played a central role, followed by transcription factors, enzymes, other protein and micro-RNA. Of them, AUX, BR and CK are positive regulator of ovule number, whereas GA acts negatively on it. Interestingly, many ovule number regulators have conserved functions across several plant taxa, which should be the targets of genetic improvement via breeding or gene editing. Many ovule number regulators identified to date are involved in the diverse biological process, such as ovule primordia formation, ovule initiation, patterning, and morphogenesis. The relations between ovule number and related characteristics/traits especially of gynoecium/fruit size, ovule fertility, and final seed number, as well as upcoming research questions, are also discussed. In summary, this review provides a general overview of the present finding in ovule number regulation, which represents a more comprehensive and in-depth cognition on it.

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.

Resistant and susceptible cacao genotypes exhibit defense gene polymorphism and unique early responses to Phytophthora megakarya inoculation

Key genes potentially involved in cacao disease resistance were identified by transcriptomic analysis of important cacao cultivars. Defense gene polymorphisms were identified which could contribute to pathogen recognition capacity. Cacao suffers significant annual losses to the water mold Phytophthora spp. (Oomycetes). In West Africa, P. megakarya poses a major threat to farmer livelihood and the stability of cocoa production. As part of a long-term goal to define key disease resistance genes in cacao, here we use a transcriptomic analysis of the disease-resistant cacao clone SCA6 and the susceptible clone NA32 to characterize basal differences in gene expression, early responses to infection, and polymorphisms in defense genes. Gene expression measurements by RNA-seq along a time course revealed the strongest transcriptomic response 24 h after inoculation in the resistant genotype. We observed strong regulation of several pathogenesis-related genes, pattern recognition receptors, and resistance genes, which could be critical for the ability of SCA6 to combat infection. These classes of genes also showed differences in basal expression between the two genotypes prior to infection, suggesting that prophylactic expression of defense-associated genes could contribute to SCA6's broad-spectrum disease resistance. Finally, we analyzed polymorphism in a set of defense-associated receptors, identifying coding variants between SCA6 and NA32 which could contribute to unique capacities for pathogen recognition. This work is an important step toward characterizing genetic differences underlying a successful defense response in cacao.

Cacao biotechnology: current status and future prospects

Theobroma cacao-The Food of the Gods, provides the raw material for the multibillion dollar chocolate industry and is also the main source of income for about 6 million smallholders around the world. Additionally, cocoa beans have a number of other nonfood uses in the pharmaceutical and cosmetic industries. Specifically, the potential health benefits of cocoa have received increasing attention as it is rich in polyphenols, particularly flavonoids. At present, the demand for cocoa and cocoa-based products in Asia is growing particularly rapidly and chocolate manufacturers are increasing investment in this region. However, in many Asian countries, cocoa production is hampered due to many reasons including technological, political and socio-economic issues. This review provides an overview of the present status of global cocoa production and recent advances in biotechnological applications for cacao improvement, with special emphasis on genetics/genomics, in vitro embryogenesis and genetic transformation. In addition, in order to obtain an insight into the latest innovations in the commercial sector, a survey was conducted on granted patents relating to T. cacao biotechnology.

Colombia a Source of Cacao Genetic Diversity As Revealed by the Population Structure Analysis of Germplasm Bank of Theobroma cacao L

Beans of the species Theobroma cacao L., also known as cacao, are the raw material to produce chocolate. Colombian cacao has been classified as a fine flavor cacao that represents the 5% of cacao world's production. Colombian genetic resources from this species are conserved in ex situ and in-field germplasm banks, since T. cacao has recalcitrant seeds to desication and long-term storage. Currently, the collection of T. cacao of the Colombian Corporation of Agricultural Research (CORPOICA) has approximately 700 germplasm accessions. We conducted a molecular analysis of Corpoica's cacao collection and a morphological characterization of some accessions with the goal to study its genetic diversity and population structure and, to select interesting accessions for the cacao's breeding program. Phenotypic evaluation was performed based on 18 morphological traits and 4 biochemical traits. PCA analysis of morphological traits explained 60.6% of the total variation in seven components and 100% of the total variation of biochemical traits in four components, grouping the collection in 4 clusters for both variables. We explored 565 accessions from Corpoica's germplasm and 252 accessions from reference populations using 96 single nucleotide polymorphism (SNP) molecular markers. Molecular patterns of cacao Corpoica's collection were obtained amplifying specific alleles in a Fluidigm platform that used integrated circuits of fluids. Corpoica's collection showed highest genetic diversity [Expected Heterozygosity (HE = 0.314), Observed Heterozygosity (HO = 0.353)] that is reduced when reference populations were included in the dataset (HE = 0.294, HO = 0.261). The collection was divided into four clusters based on population structure analysis. Cacao accessions from distinct groups showed some taxonomic concordance and reflected their geographic origins. For instance, accessions classified as Criollo were clearly differentiated in one group and we identified two new Colombian genetic groups. Using a number of allelic variations based on 87 SNP markers and 22 different morphological/biochemical traits, a core collection with a total of 232 accessions was selected as a primary genetic resource for cacao breeders.

Association mapping of seed and disease resistance traits in Theobroma cacao L

Microsatellite and single nucleotide polymorphism markers that could be used in marker assisted breeding of cacao were identified for number of filled seeds, black pod resistance and witches' broom disease resistance. An association mapping approach was employed to identify markers for seed number and resistance to black pod and witches' broom disease (WBD) in cacao (Theobroma cacao L.). Ninety-five microsatellites (SSRs) and 775 single nucleotide polymorphisms (SNPs) were assessed on 483 unique trees in the International Cocoa Genebank Trinidad (ICGT). Linkage disequilibrium (LD) and association mapping studies were conducted to identify markers to tag the phenotypic traits. Decay of LD occurred over an average 9.3 cM for chromosomes 1-9 and 2.5 cM for chromosome 10. Marker/trait associations were generally identified based on general linear models (GLMs) that incorporated principal components from molecular information on relatedness factor. Seven markers (mTcCIR 8, 66, 126, 212; TcSNP368, 697, 1370) on chromosomes 1 and 9 were identified for number of filled seeds (NSEED). A single marker was found for black pod resistance (mTcCIR280) on chromosome 3, whereas six markers on chromosomes 4, 5, 6, 8, and 10 were detected for WBD (mTcCIR91, 183; TcSNP375, 720, 1230 and 1374). It is expected that this association mapping study in cacao would contribute to the knowledge of the genetic determinism of cocoa traits and that the markers identified herein would prove useful in marker assisted breeding of cacao.

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.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2415-x) contains supplementary material, which is available to authorized users.

Sequencing of a QTL-rich region of the Theobroma cacao genome using pooled BACs and the identification of trait specific candidate genes

Background

BAC-based physical maps provide for sequencing across an entire genome or a selected sub-genomic region of biological interest. Such a region can be approached with next-generation whole-genome sequencing and assembly as if it were an independent small genome. Using the minimum tiling path as a guide, specific BAC clones representing the prioritized genomic interval are selected, pooled, and used to prepare a sequencing library.

Results

This pooled BAC approach was taken to sequence and assemble a QTL-rich region, of ~3 Mbp and represented by twenty-seven BACs, on linkage group 5 of the Theobroma cacao cv. Matina 1-6 genome. Using various mixtures of read coverages from paired-end and linear 454 libraries, multiple assemblies of varied quality were generated. Quality was assessed by comparing the assembly of 454 reads with a subset of ten BACs individually sequenced and assembled using Sanger reads. A mixture of reads optimal for assembly was identified. We found, furthermore, that a quality assembly suitable for serving as a reference genome template could be obtained even with a reduced depth of sequencing coverage. Annotation of the resulting assembly revealed several genes potentially responsible for three T. cacao traits: black pod disease resistance, bean shape index, and pod weight.

Conclusions

Our results, as with other pooled BAC sequencing reports, suggest that pooling portions of a minimum tiling path derived from a BAC-based physical map is an effective method to target sub-genomic regions for sequencing. While we focused on a single QTL region, other QTL regions of importance could be similarly sequenced allowing for biological discovery to take place before a high quality whole-genome assembly is completed.

Development of a marker assisted selection program for cacao

ABSTRACT Production of cacao in tropical America has been severely affected by fungal pathogens causing diseases known as witches' broom (WB, caused by Moniliophthora perniciosa), frosty pod (FP, caused by M. roreri) and black pod (BP, caused by Phytophthora spp.). BP is pan-tropical and causes losses in all producing areas. WB is found in South America and parts of the Caribbean, while FP is found in Central America and parts of South America. Together, these diseases were responsible for over 700 million US dollars in losses in 2001 (4). Commercial cacao production in West Africa and South Asia are not yet affected by WB and FP, but cacao grown in these regions is susceptible to both. With the goal of providing new disease resistant cultivars the USDA-ARS and Mars, Inc. have developed a marker assisted selection (MAS) program. Quantitative trait loci have been identified for resistance to WB, FP, and BP. The potential usefulness of these markers in identifying resistant individuals has been confirmed in an experimental F(1) family in Ecuador.

Adding value to cocoa (Theobroma cacao L.) germplasm information with domestication history and admixture mapping

A sound understanding of crop history can provide the basis for deriving novel genetic information through admixture mapping. We confirmed this, by using characterization data from an international collection of cocoa, collected 25 years ago, and from a contemporary plantation. We focus on the trees derived from three centuries of admixture between Meso-American Criollo and South American Forastero genomes. In both cacao sets of individuals, linkage disequilibrium extended over long genetic distances along chromosome regions, as expected in populations derived from recent admixture. Based on loose genome scans, genomic regions involved in useful traits were identified. Fifteen genomic regions involved in seed and fruit weight variation were highlighted. They correspond to ten previously identified QTLs and five novel ones. Admixture mapping can help to add value to genetic resources and thus, help to encourage investment in their conservation.

Creation of BAC genomic resources for cocoa ( Theobroma cacao L.) for physical mapping of RGA containing BAC clones

We have constructed and validated the first cocoa ( Theobroma cacao L.) BAC library, with the aim of developing molecular resources to study the structure and evolution of the genome of this perennial crop. This library contains 36,864 clones with an average insert size of 120 kb, representing approximately ten haploid genome equivalents. It was constructed from the genotype Scavina-6 (Sca-6), a Forastero clone highly resistant to cocoa pathogens and a parent of existing mapping populations. Validation of the BAC library was carried out with a set of 13 genetically-anchored single copy and one duplicated markers. An average of nine BAC clones per probe was identified, giving an initial experimental estimation of the genome coverage represented in the library. Screening of the library with a set of resistance gene analogues (RGAs), previously mapped in cocoa and co-localizing with QTL for resistance to Phytophthora traits, confirmed at the physical level the tight clustering of RGAs in the cocoa genome and provided the first insights into the relationships between genetic and physical distances in the cocoa genome. This library represents an available BAC resource for structural genomic studies or map-based cloning of genes corresponding to important QTLs for agronomic traits such as resistance genes to major cocoa pathogens like Phytophthora spp ( palmivora and megakarya), Crinipellis perniciosa and Moniliophthora roreri.

A new cacao linkage map based on codominant markers: development and integration of 201 new microsatellite markers

A linkage map of cacao based on codominant markers has been constructed by integrating 201 new simple sequence repeats (SSR) developed in this study with a number of isoenzymes, restriction fragment length polymorphisms (RFLP), microsatellite markers and resistance and defence gene analogs (Rgenes-RFLP) previously mapped in cacao. A genomic library enriched for (GA)(n) and (CA)(n) was constructed, and 201 new microsatellite loci were mapped on 135 individuals from the same mapping population used to establish the first reference maps. This progeny resulted from a cross between two heterozygous cacao clones: an Upper-Amazon Forastero (UPA 402) and a Trinitario (UF 676). The new map contains 465 markers (268 SSRs, 176 RFLPs, five isoenzymes and 16 Rgenes-RFLP) arranged in ten linkage groups corresponding to the haploid chromosome number of cacao. Its length is 782.8 cM, with an average interval distance between markers of 1.7 cM. The new microsatellite markers were distributed throughout all linkage groups of the map, but their distribution was not random. The length of the map established with only SSRs was 769.6 cM, representing 94.8% of the total map. The current level of genome coverage is approximately one microsatellite every 3 cM. This new reference map provides a set of useful markers that is transferable across different mapping populations and will allow the identification and comparison of the most important regions involved in the variation of the traits of interest and the development of marker-assisted selection strategies.