Identification of candidate genes involved in Witches’ broom disease resistance in a segregating mapping population of Theobroma cacao L. in Brazil
- Stefan Royaert†1,
- Johannes Jansen†2,
- Daniela Viana da Silva3,
- Samuel Martins de Jesus Branco3,
- Donald S. LivingstoneIII4,
- Guiliana Mustiga4,
- Jean-Philippe Marelli1,
- Ioná Santos Araújo5,
- Ronan Xavier Corrêa3 and
- Juan Carlos Motamayor4Email authorView ORCID ID profile
© Royaert et al. 2016
Received: 14 May 2015
Accepted: 26 January 2016
Published: 11 February 2016
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.
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’.
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.
Theobroma cacao L. (cacao), which is native to the Amazonian basin of South America [1, 2], is a perennial tree that belongs to the Malvaceae sensu lato  and has a diploid chromosome number of 2n = 2x = 20 . Recently, two cacao genomes have been published. ‘B97-61/B2’, a member of the Criollo genetic group, has an estimated genome size of 409 Mbp  whereas the ‘Matina 1–6’ genome, a member of the Amelonado genetic group and more representative of the cacao that is cultivated worldwide, has an estimated genome size of 445 Mbp .
Witches’ broom disease (WBD), which causes significant losses, is only present in South America and in the Caribbean  and was first reported in 1785 by the explorer Alexandre Rodrigues in the Brazilian Amazon basin . The fungus Moniliophthora perniciosa (Stahel) Aime & Phillips-Mora is a basidiomycete and produces basidiocarps, which emerge from dry brooms in the form of mushroom-like structures. Direct production losses are due to infected mature pods and flower cushions causing parthenocarpic fruits, thereby reducing significantly the potential production .
A major source of resistance was found in the clones ‘SCAVINA 6’ (‘SCA 6’) and ‘SCAVINA 12’ (‘SCA 12’)  and in two other clones, ‘CAB 208’ and ‘CAB 214’, which are genetically different from ‘SCA 6’ . QTL associated with resistance to vegetative brooms (VB) have been identified in three mapping populations. The first population was developed by selfing one tree, ‘TSH 516’, which was obtained from a cross between the resistant ‘SCA 6’ and the susceptible ‘ICS 1’. A major QTL was identified on chromosome IX [12–14], whereas a minor QTL was identified on chromosome I . The two other populations consist of the progenies of a cross between the susceptible ‘ICS 39’ and the resistant ‘CAB 208’ and of a cross between ‘ICS 39’ and the resistant ‘CAB 214’ . One QTL was identified on chromosome IX and is located at a slightly different position than the QTL that was identified in the F2 mapping population of ‘SCA 6’ × ‘ICS 1’. Two other QTL were identified on chromosomes IV and VIII .
One of the main requirements of an efficient marker-assisted selection (MAS) program is to have reliable markers associated with traits of interest. The availability of dense genetic linkage maps is one of the prerequisites. The first linkage map in cacao, which contained 202 markers, was published in 1995 . This reference map was improved over time through the addition of more trees and markers [16–19]. The most recent consensus map, which also has the highest marker density, combined the data of the cross ‘UPA402’ × ‘UF676’ and of the F2 mapping population resulting from selfing a single tree of the cross ‘SCA 6’ × ‘ICS 1’. This map contains 1,262 co-segregating markers and has a length of 733.6 cM .
Here, we describe the development of a dense genetic linkage map using single nucleotide polymorphisms (SNP) based on 459 trees of a cross between the WBD-resistant ‘TSH 1188’ and the tolerant ‘CCN 51’ at the Mars Center for Cocoa Science (MCCS) in Barro Preto, Bahia, Brazil. Furthermore, we identified seven QTL on five different chromosomes that are associated with WBD using a multi-trait QTL analysis for outbreeders. Various potential candidate genes that were identified within the QTL may be involved in different steps in disease resistance. Potential SNP markers that can be applied in the marker-assisted selection of WBD-resistant material are proposed.
Plant material and design of the field experiment
For this study, a segregating population ‘MP01’ consisting of 459 offspring from a cross between ‘TSH 1188’ and ‘CCN 51’ was used. Both parents are very different for various traits of interest and the offspring segregates for pod color , self-compatibility, resistance to WBD, black pod and Ceratocystis wilt, fat content and fat composition, methylxanthines (theobromine and caffeine), flavanol content and some other agronomic traits. ‘TSH 1188’ is considered resistant to WBD, whereas ‘CCN 51’ is considered moderately resistant. A recent study focused on the number of cushion brooms showed that ‘TSH 1188’ had about 3 % infected flower cushions, whereas ‘CCN 51’ has close to 30 % infected flower cushions . Individual trees were planted in 2000 in a 3 × 3 m grid in a hydromorphic and typical tropudalf (Itabuna modal) mixed soil type. Every fifth row in the field contained 5–19 trees (depending on the length of the rows) of the variety ‘Comum’, which is susceptible to WBD and serves as a natural and permanent inoculum source. The layout of the field experiment is shown in Additional file 1. Shade was provided using the traditional ‘cabruca’ system in which the trees are grown amongst the Atlantic Forest's native canopies.
SNP identification, SNP array evaluation and genotyping of trees using the Illumina Infinium SNP chip
DNA for genotyping was isolated from leaf samples according to the protocol described by Motamayor et al. . The SNPs that were used on the cacao Illumina Infinium SNP chip were identified as part of the work to identify conserved ortholog set II (COSII) SNP markers for MAS in cacao as described by Kuhn et al.  and Livingstone et al. . The strategy was to sequence leaf RNA samples of 15 members of a panel of diverse cacao genotypes and to align the sequences against the Theobroma cacao Matina 1–6 leaf transcriptome. After alignment, a variant report was generated for all of the identified SNPs . Further filtering of the SNPs was performed to select the final 6,000 SNPs to produce the 6 K SNP chip . Of the 6,000 SNPs that were submitted for inclusion on the bead array, 5,388 were successfully synthesized. Failing to produce genotypic data for any sample, 174 SNPs were removed before analysis. Additionally, to serve as a quality control for genotype calling, only the markers with GenTrain and 10 % GenCall scores greater than 0.4 and 0.2, respectively, were included for data analysis. This process resulted in a final set of 5,149 SNPs. The Cacao6kSNP array was originally run with 1,152 DNA samples (including duplicated controls), but seven DNA samples failed completely. In total, 5,149 SNPs were genotyped in 1,145 DNA samples for a total of 5,895,605 data points or SNP genotypes .
SNP genotype data analysis and construction of a genetic map
JoinMap®4.1  was used to create the genetic map. SNPs with more than 5 % missing data were discarded. SNPs were also checked for strongly deviating segregation ratios. For each linkage group, the marker data were analyzed per segregation type (markers segregating only in the mother ‘TSH 1188’, markers segregating only in the father ‘CCN 51’, and cosegregating markers). Finally, integrated genetic maps were obtained using the Maximum Likelihood (ML) mapping algorithm . Markers with a nearest-neighbor stress greater than 2 were discarded, followed by recalculation of the marker order. Map distances were calculated using the Haldane mapping function and the resulting maps were drawn using MapChart 2.1 .
Phenotypic evaluation of WBD resistance
VB and cushion brooms (CB) were counted once or twice a year over a period of 4 years (April 2008, March 2009, March 2010, December 2010 and December 2011) and removed after they had been counted. The data from the 459 genotyped trees were recorded as the number of VB per year and the number of CB per year.
The distributions of the numbers of VB and CB per year were extremely skewed with many zeros. Therefore, the presence/absence of brooms, being the main feature of the phenotypic variation, was used for further analysis. To get an indication of the genetic component of variation for the number of VB and CB over time, determination coefficients were calculated from the variance components for genotypes and the genotype-by-year interaction using the REML facilities of GenStat 15 .
The number of years in which each tree carried VB and the number of years in which it carried CB were used for further analysis; the values for each of these traits varied between zero and four. To remove the effect of possible spatial patterns of infestation by WBD, the data was adjusted for row and column effects using the REML facilities of GenStat 15 . The adjusted data of the two traits was used for a multi-trait QTL analysis , which was carried out using the facilities for QTL analysis in GenStat 15. The QTL profiles were obtained using all of the available markers using the integrated linkage map. The values of the test statistic (expressed as -log10(P)) are based on the additive effect of ‘TSH 1188’, the additive effect of ‘CCN 51’, and the interaction effect. The initial selection of QTL was performed by simple interval mapping using a threshold of three for the -log10(P), and the final selection of QTL was performed by backward elimination using the default parameters provided by GenStat.
To obtain a better sense of the precision of the QTL positions and their usefulness for cacao breeding, a simulation study was carried out in which random samples of 200 individuals were drawn from the original 459 individuals. For each sample, a QTL analysis was carried out as described above. The sample size of 200 was chosen as a reasonable size compared to the population size of many QTL experiments. For each sample, the marker positions that were selected as QTL were recorded as ‘hit’. The results of 500 random samples were stored, and a graphical display of the results was obtained using the procedure DQSQTLSCAN of GenStat 15.
Identification of associations between haplotypes and WBD resistance
Haplotypes of the offspring and the parental haplotypes surrounding the SNP markers on each chromosome with a significant association with WBD resistance in the QTL analysis were identified using version 1.1 of the Matina 1–6 genome and iXora [6, 28]. To test the significance of the haplotype-phenotype associations, a chi-squared test was applied independently for each identified marker, based on the number of total VB as a trait, where the trees with more than 10 brooms were considered susceptible. ‘TSH 1188’ has two parental haplotypes that are designated as T1 and T2, whereas ‘CCN 51’ has the two parental haplotypes, C1 and C2. Chi-squared tests were also used to test whether the parental haplotype combinations that were present in the offspring (T1-C1, T1-C2, T2-C1 and T2-C2) were non-randomly associated with WBD resistance.
Origin of disease resistance alleles
The phasing of 1162 individuals including three mapping populations and a set of diversity panel members  from distinct T. cacao structural groups was run with fastPHASE . Sixty-two markers (from Tcm009s07051488 to Tcm009s09263083) in the WBD QTL9.1 region were used for the phasing. To improve the accuracy of phasing and to account for haplotype structure, a subpopulation label was used in the form of an integer. Known members that were sampled from the same population (i.e., mapping populations) were assigned the same population, while unrelated individuals were assigned a distinct integer. The expectation-maximization (EM) algorithm for computing the maximum likelihoods was controlled by the following options: 20 random starts, 25 iterations, and 200 haplotypes that were sampled from the posterior distribution from a particular random start. The default allelic two-parameter error model for inferring true genotypes was also used to scan for genotype errors. To create the phylogenetic neighbor joining (NJ) tree, thirty diversity panel members were chosen, including the parents from MP01, CATIE type 1, and PNG mapping populations . The distance matrix for phylogeny estimation was created with the MEGA 6 phylogeny reconstruction setting, using the method of Maximum Composite Likelihood using 1000 bootstrap samples .
Candidate gene identification
Using the location of the identified SNP markers with the highest association in the QTL analysis and the gene annotations in the version 1.1 Matina 1–6 genome , a search was performed manually for possible candidate genes involved in disease resistance. The specific keywords that were searched for were resistance, senescence, jasmonic acid, salicylic acid (SA), ethylene (ET), gibberellic acid, auxin, reactive oxygen species (ROS), mitogen-activated protein kinases (MAPKs), F-box, LRR-receptor, cysteine, serine/threonine protein kinase, phosphatase, ubiquitin, NPR1, and WRKY. All of the genes within an interval of approximately 2 cM surrounding the most significant marker for each QTL were investigated in more detail using BLASTX and additional literature study. All details about the experimental set-up and data analysis for the differential gene expression of a subset of candidate genes are given in the Additional file 2 and in Additional file 3.
Results and discussion
Construction of the integrated genetic map and relationship with the physical map
An overview of the mapping process
Number of SNPs
Map length (cM)
ab × aa
aa × ab
ab × ab
ab × aa
aa × ab
ab × ab
Compared to the most recent map , the MP01 map contains three times more markers (3,526 vs. 1,043), has a map length that is 12 % longer (852.8 cM vs. 751.7 cM) and has a marker density that is three times denser (0.24 cM vs. 0.72 cM). The shortest chromosome was chromosome VII (56.5 cM and 156 SNPs), and the longest chromosome was chromosome IX (114.0 cM and 522 SNPs). The average gap without SNP markers was 2.3 cM, with the smallest gap of 1.8 cM on chromosomes III and VIII and the largest gap of 3.0 cM on chromosome X. The average distance between two markers varied from 0.19 cM on chromosome IV to 0.37 cM on chromosome X (Table 1 , Additional file 5). In contrast to the previously published studies [15–19, 31], no distorted segregation ratios were observed; all SNPs of the type ‘ab × aa’ and ‘aa × ab’ showed a segregation ratio between 0.4 and 0.6 (Additional file 6). The relationship between the positions of the SNPs on the integrated linkage map and the corresponding position on the physical map is shown in Additional file 7. Apart from a few deviating points, a nearly perfect correspondence exists between the linkage map and the physical map. All ten chromosomes show physical regions with reduced recombination rates, presumably the centromeric regions, which agrees with the karyotyping results and demonstrates the lack of a long arm for some of the chromosomes and that chromosomes VI, VII, VIII and X are shorter than the other chromosomes .
Phenotypic evaluation of WBD resistance
‘TSH 1188’ is considered resistant to WBD, whereas ‘CCN 51’ is considered moderately resistant. A recent study in MP01 focused on the number of cushion brooms showed that ‘TSH 1188’ had about 3 % infected flower cushions, whereas ‘CCN 51’ had close to 30 % infected flower cushions . Nowadays we notice that ‘CCN 51’ has many more infected flower cushions located on the stem and the lower part of the main branches, compared to ‘TSH 1188’, but the latter one seems to have more infected flower cushions higher up in the canopy. Nevertheless, the number of vegetative brooms is much lower in ‘TSH 1188 than in ‘CCN 51’. It is also noteworthy that both clones are as good as resistant for WB infection of the pods. There also exist region-specific differences, which are most likely associated with differences in existing pathogen strains. It is therefore important to mention that the outcome of a plant-pathogen interaction does not depend solely on the plant genotype and it may change completely depending on the strain of the pathogen.
Proportions of trees in MP01 with vegetative brooms and cushion brooms in the years 2008–2011
Identification of QTL for WBD resistance
Effects of selected QTL, values of the significance level -log10(P) and % variance accounted for by including markers in the final QTL model
TSH 1188 × CCN 51
The major QTL (for both VB and CB) on chromosome IX most likely corresponds to the QTL identified in the F2 population obtained by selfing an individual of a cross between the WBD-resistant ‘SCA 6’ and the WBD-susceptible ‘ICS 1’ [12–14]. This result might not be surprising because ‘SCA 6’ is present in the parental lineage of ‘TSH 1188’. ‘TSH 1188’ results from a cross between ‘POUND 18’ and ‘TSH 753’ [32, 33]; ‘TSH 753’ is the result of an open-pollination of ‘TSA 641’ [33, 34], and ‘TSA 641’ is the result of a cross between ‘SCA 6’ and ‘IMC 67’ [33, 35]. The second QTL on chromosome IX had an additive effect in both ‘TSH 1188’ and ‘CCN 51’ for VB, but had an effect only in ‘CCN 51’ for CB. Other ‘CCN 51’ specific effects for CB were identified for the QTL on chromosome III and on chromosome IV. Similar effects in ‘CCN 51’ for various other QTL were noticed, indicating that these newly discovered QTL might contribute to the known moderate resistance/tolerance of ‘CCN 51’ to WBD and can be used as alternative sources of resistance. ‘CCN 51’ is derived from crosses involving ‘IMC 67’ and ‘ICS 95’ . ‘IMC 67’ shows resistance or susceptibility to WBD according to different studies [33, 37, 38], whereas ‘ICS 95’ has intermediate resistance [33, 39, 40]. At MCCS, the two clones show an intermediate level of resistance to WBD. The fact that about 10–17 % of the trees in the offspring doesn’t have any CB or VB (Additional file 8), suggests that they may have a combination of favorable alleles originating from ‘SCA 6’, ‘IMC 67’ and ‘ICS 95’. The higher number of resistant trees (about 65 %) in MP01 might indicate that they have a combination of favorable alleles from any of these three clones. A preliminary study performed by Santos et al.  used 18 SSRs (some linked to resistance in the F2 mapping population [12–14]) on 16 trees (eight resistant and eight susceptible) of a subset of 50 trees of MP01 under natural and artificial infection pressure with M. perniciosa. A simple regression analysis indicated that in addition to the known markers on chromosome IX, there were also markers that showed an association with resistance on chromosomes I, III and IV. The markers on chromosomes III and IX were from the ‘TSH 1188’ parent, whereas the markers on chromosomes I and IV were from the ‘CCN 51’ parent, most likely indicating an additional source of resistance . In the present study, however, no QTL were located at chromosome I. The identified QTL on chromosomes III and IV were located on the opposite side of the chromosomes compared to our results. Alternative QTL were identified in ‘CAB 208’ and ‘CAB 214’, which were originally collected from the Brazilian Amazon and did not correspond to the QTL that were identified in this study . These results suggest that there are multiple sources of resistance in the different mapping populations for both VB and CB. To better understand the precision of the QTL positions and their usefulness for cacao breeding, a resampling study was carried out in which random samples of 200 individuals were drawn 500 times from the original 459 individuals, and the results are shown in Additional file 10. The sample size of 200 individuals was chosen because most of the published QTL results in cacao are based on a population size ranging from 95 to 264 individuals [12, 15–19, 31, 42–44]. For the region between 40.8 cM and 45.1 cM of chromosome IX, the average percentage of samples in which a marker position is selected as a QTL (= %hit) was 84.0 %; nearly half of these samples (41.2 %) were at 44.4 cM. For a slightly larger region between 37.9 cM and 47.5 cM, the %hit was 94.0 %. For the region between 5.2 and 18.3 cM of chromosome IX, the %hit was 15.8 %. The hits on the other chromosomes were more dispersed (Additional file 10). From the simulation study, it became clear, however, that the more trees are used in the study, the more precise the QTL analysis becomes, smaller QTL effects can detected, but also the positions of the QTL with small effects are not very accurate.
Identification of associations between the parental haplotypes and WBD resistance
Alleles for each of the parental haplotypes (T1, T2, C1 and C2) of ‘TSH 1188’ and ‘CCN 51’ for each QTL and representative SNP marker associated with WBD
Chi-squared test to calculate the associations between the parental haplotypes (T1, T2, C1 and C2) and the observed and expected number of resistant/susceptible trees
QTL and haplotypes
Chi-squared test to calculate the associations between the parental haplotype combinations (T1-C1, T1-C2, T2-C1 and T2-C2) and the number of resistant/susceptible trees
To identify the best marker combinations of two QTL for the selection of resistant trees, each haplotype combination of QTL9.1 was combined with each haplotype combination for the other QTL. All the associated P-values are presented in Additional file 11, but the five highest associations for a combination of two QTL are presented in Table 6. The highest association was found for QTL9.1_T2-C1 with QTL9.2_T2-C1 (P = 1.33 × 10−9), followed by QTL9.1_T2-C2 with QTL6.1_T2-C2 (P = 5.79 × 10−8).
To identify the best marker combinations of three QTL for the selection of resistant trees, two haplotype combinations of QTL9.1 (T2-C1 and T2-C2) were combined with the more-significant haplotype combinations for two other QTL. All of the associated P-values are presented in Additional file 11, but the five highest associations for a combination of three QTL are presented in Table 6. The highest association was found for the combination of QTL9.1_T2-C2 with QTL6.1_T2-C2 and QTL6.2_T2-C2 (P = 0.000057), followed by the combination of QTL9.1_T2-C2 with QTL9.2_T2-C2 and QTL6.1_T2-C2 (P = 0.000074). The percentage of resistant trees for the top five of each QTL combination varied between 76.9 and 95.0 % depending on the number of trees that were available with that specific combination (Table 6).
Furthermore, we used the trees that showed a recombination event between the different haplotypes within a small region around the most significant marker and extended this region by adding flanking SNP markers. This extension was designed to help us to localize more precisely the most likely candidate genes, but we only managed to do this successfully for the most significant QTL9.1 on chromosome IX (Additional file 12). Between the markers Tcm009s08010009 and Tcm009s08051015, we identified three trees (MP01-212, 555 and 795) with a recombination event between the maternal haplotypes. Between Tcm009s08051015 and the most significant marker Tcm009s08066239, we identified a single tree (MP01-218) with a recombination event between the maternal haplotypes. One recombination event (MP01-462) was observed between Tcm009s08066239 and Tcm009s08105371. Furthermore, we identified two recombination events (MP01-243 and 593) between the markers Tcm009s08105371 and Tcm009s08172795 (Additional file 12). Based on the recombinants, the size of QTL9.1 between the first and the last SNP marker is 162,786 bp and contains five candidate genes possibly involved in disease resistance (see below).
For the most significant marker (Tcm009s08066239) and its neighbor (Tcm009s08051015), all of the trees that have the T2 haplotype, i.e., the G allele or the T allele, respectively, had ten or fewer total VB over the 4 years and can be considered resistant (Additional file 12). The resistant trees that carry the A allele (in orange, Additional file 12) possess one or more associated alleles of the other identified QTL (data not shown). The susceptible trees (more than ten total VB, in red, Additional file 12), with the exception of MP01-114, do not possess any of the associated alleles of the other identified QTL. MP01-114 contains the favorable alleles of QTL6.1_T2-C2, QTL6.2_T2-C2 and QTL4.1_T2-C1. Sixty-four trees with a single recombination event were identified in a larger interval from 7,293,462 to 9,668,548 bp. Of the 18 susceptible trees within this interval, only MP01-281 has the G allele, i.e., only 6 % of the trees carrying the G allele are at risk of being misidentified as being resistant based on the allele information (data not shown).
This result also indicates that these allele-specific markers, either by themselves or in combination with some of the other allele-specific markers of the other QTL, can be used for the diagnostic disease resistance screening of other populations and germplasm collections. The associated SNP markers can also be relatively easily converted into PCR-based markers using the dCAPS method .
Origin of disease resistance alleles
Potential candidate genes possibly involved in Witches’ broom resistance through programmed cell death
All of the genes within an interval of 2–3 cM (approximately 300,000 to 700,000 bp) surrounding the most significant marker for each QTL were investigated. The number of genes within each QTL ranged from 23 in QTL9.1 to 90 in QTL6.1, and the number of candidate genes that were potentially involved in disease resistance ranged from 6 in QTL9.1 to 29 in QTL6.1 (Additional file 13). The potential biological processes for each of the genes were determined by GO Slim annotation (Additional file 14).
One of the first steps in the recognition of pathogens is the activation of plant nucleotide-binding (NB)-leucine-rich repeat (LRR) receptors by direct or indirect mechanisms . On QTL6.1, three putative disease resistance-related proteins were identified that contain the NB and LRR domains. Thecc1EG028968 and Thecc1EG028973 encode a CC-NBS-LRR resistance protein, whereas Thecc1EG028972 encodes an NB-ARC domain-containing disease resistance protein. A fourth gene, Thecc1EG028979 encodes a disease resistance-responsive family protein. The SNP Tcm006s19715703, which is associated with QTL6.1, is located within Thecc1EG028962, a gene encoding a hydroxyproline-rich glycoprotein family protein (HRGP). Cross-linking of HRGPs is an important process by which to strengthen the cell walls, contributing to plant defense reactions . Another gene, Thecc1EG028959, in QTL6.1 is the leaf senescence-associated receptor-like protein kinase, which has homology to the senescence-induced receptor-like serine/threonine protein kinase (SIRK). This gene is specifically expressed in the leaves at senescence and after pathogen infection and is a putative target of WRKY6 , which is a transcription factor that is activated through the MAPK pathway . On QTL6.2, the most associated SNP, Tcm006s25375496, is within Thecc1EG030009, a MYB-like 102 gene that does not seem to function in disease resistance.
The next stages of the plant immune response include a rapid influx of calcium ions, an oxidative burst of ROS and the activation of MAPKs [50–52]. Four genes that were identified in QTL7.1 have been annotated as 2-oxoglutarate and Fe(II)-dependent oxygenase superfamily proteins (Thecc1EG032367, Thecc1EG032371, Thecc1EG032376 and Thecc1EG032382), which are involved in the formation of hydrogen peroxide at different stages in cacao depending on the disease status of the trees (resistant vs. susceptible) [53, 54]. On QTL7.1, the most associated SNP, Tcm007s10302466, is within Thecc1EG032345, a nucleotide/sugar transporter protein, which does not seem to function in disease resistance. Another neighboring gene, Thecc1EG032344, is annotated as an aldolase-type TIM barrel family protein, which in Arabidopsis thaliana seems to be involved in defense response to bacteria and in the hydrogen peroxide biosynthetic process .
The MAPK signaling cascade is involved in the biosynthesis and signaling of plant stress hormones, ROS generation, stomatal closure, defense gene activation, phytoalexin biosynthesis, cell wall strengthening and programmed cell death (PCD) . In A. thaliana, one of the genes involved in this cascade is MAPK kinase 5, which is also one of the genes (Thecc1EG016564) that we identified in QTL3.1 [46, 48, 51]. Another disease resistance-related gene, Thecc1EG016546, in QTL3.1 encodes a glycosyl hydrolase superfamily protein. In crucifers, these enzymes break down glucosinolates and convert them into biologically active compounds. These compounds restrain insect feeding activity on plants and have potent antimicrobial activity, mainly against necrotrophic fungi [56, 57]. NDR1 (Thecc1EG016558), which is in the same QTL, was shown to be required for disease resistance in A. thaliana to both bacterial and fungal pathogens . The SNP marker Tcm003s33466269, which is associated with QTL3.1, is located within Thecc1EG016568, a gene encoding an RNA-binding (RRM/RBD/RNP motifs) family protein.
Another important pathway in disease response and PCD is the ubiquitin/26S proteasome-mediated pathway, which is involved in the selective degradation of proteins in cells of eukaryotic organisms [59–61]. At least five U-box E3 ligases have been identified that are involved in disease response to various pathogens in different plants species [62–67]. For example, disruption of the A. thaliana Plant U-box protein 13 (PUB13) by T-DNA insertion results in spontaneous cell death, the accumulation of hydrogen peroxide and SA, and elevated resistance to biotrophic pathogens, as well as increased susceptibility to necrotrophic pathogens . Tcm009s08066239, the SNP in QTL9.1 with the best association with WBD, is located into Thecc1EG038262 that has been annotated as a RING/U-box superfamily protein. Another RING/U-box superfamily protein (Thecc1EG038259) has been identified within the same QTL. Three other potential candidate genes in this QTL have been annotated as follows: DNA-binding protein phosphatase 1 (Thecc1EG038249, dephosphorylation and inactivation of MAPK ), uveal autoantigen with coiled-coil domains and ankyrin repeats (UACA, Thecc1EG038261), and serine/threonine protein kinase CTR1/EDR1 with an octicosapeptide/Phox/Bem1p domain (Thecc1EG038267). In humans, UACA is involved in a series of molecular signals in which an intracellular signal is conveyed to trigger apoptotic cell death . This pathway is induced by the detection of DNA damage, which can be caused by ROS, thereby linking to one of the primary symptoms of WBD. CTR1 is a Raf-like MAPKKK that is involved in ET signal transduction [69, 70], and EDR1 was found to encode a MAPKKK that is similar to CTR1. This MAPKKK is involved in the induction of several defense responses, including host cell death, and it negatively regulates SA-dependent defense responses, abscisic acid (ABA) signaling, and ET-induced senescence . Tcm009s02031341, the SNP in QTL9.2 with the best association with WBD, is located within Thecc1EG047089 that has been annotated as encoding an RNA recognition motif (RRM)-containing protein. It is interesting that, according to STRING v9.1 , this protein in A. thaliana may interact with a TIR-NBS-LRR class disease resistance protein. Among the other candidate genes, two were annotated as BURP domain-containing proteins. These proteins are involved in abiotic and biotic stresses in Brassica .
Ca2+ and MAPKs also control the synthesis and/or signaling of the hormones that are involved in the later steps of defense gene activation and PCD. The SNP marker Tcm004s00110232 that is associated with QTL4.1 is located within Thecc1EG016777 encoding a UDP-glucose pyrophosphorylase 2. Overexpression of this enzyme in poplar showed extremely high quantities of the glycoside of salicylic acid (SAG). SA accumulation, and therefore SAG, could catalyze the increase in other plant defense metabolites . Also in QTL4.1 Thecc1EG016815, the senescence-associated gene 101, is part of the Enhanced Disease Susceptibility 1 (EDS1) pathway and is responsible for SA accumulation after infection with biotrophic pathogens of A. thaliana .
Of the 11 candidate genes selected for differential gene expression (Additional file 3), only six showed considerable expression differences between infected and mock-infected plants for both ‘CCN 51’ and ‘TSH 1188’. The most remarkable expression differences were seen for the UACA gene. When comparing the infected ‘CCN 51’ plants, it was clear that this gene was approximately 4-fold upregulated at 3 and 12 h after infection (HAI), whereas it had a basal level of expression at the other time points. In the resistant ‘TSH 1188’, this gene was approximately 8-fold downregulated at 3 HAI, whereas it had a similar expression profile for the other time points (see Additional files 15 and 16). According to the study of Sena et al.  this early response may correspond with the early infection stage of the fungus, where the basidiospores start to germinate around 2 HAI in the susceptible ‘CATONGO’, and around 4 HAI in the resistant ‘TSH 1188’, but both reach their maximum at 6 HAI. Penetration into the plant tissues was observed at 6 HAI for both susceptible and resistant genotypes, and in susceptible genotypes, primary hyphae were observed in the cortex under the epidermis at 48 HAI .
A study performed by Teixeira et al.  showed a detailed transcriptome analysis of the interaction between M. perniciosa and cacao at 30 days after artificial infection (which corresponds with the green broom stage). They identified 1,269 upregulated and 698 downregulated differentially expressed genes in cacao. Of those differentially expressed genes 26 were located within our identified QTL. Three of those genes were pathogen-related (CDG0008741 in QTL6.2) or receptor-like kinases involved in signal transduction (CDG00002586 and CDG00002587 in QTL6.1). The later one showed a positive fold change of 25.37 compared to the untreated control. This gene is also located between our best SNP marker on QTL6.1 (Tcm006s19715703) and the closest neighboring SNP (Tcm006s19387431). Three other genes, CGD0026033 in QTL6.2, CDG0031650 and CDG0031682 in QTL9.2, are involved in photosynthetic, carbon and nitrogen metabolism, more specifically in sugar transport, beta-oxidation of fatty acids and amino acid metabolism, respectively . It is interesting that many of these identified genes are actually very close to the most significant SNPs in each QTL identified in our study (Additional file 13). It would be of interest as well to investigate all these genes identified by Teixeira et al.  during the very early stages of infection.
One of the main requirements of an efficient MAS program in cacao is to have reliable markers associated with the trait(s) of interest. These markers and ultimate knowledge of the underlying genes must be identified first. This knowledge demands a combination of large mapping populations with very distinct parents that segregate for the trait(s) of interest, a dense set of mapped markers, and rigorous phenotyping capabilities. The availability of a sequenced genome helps in the identification and annotation of potential candidate genes, although it is important to mention here that annotations are based on the Matina 1–6 genome, which is susceptible to WBD. It is therefore not unlikely that one or more genes responsible for the resistance phenotype occur specifically in the resistant cultivars, which might render the use of the Matina 1–6 genome inappropriate. We combined these requirements to identify several genomic regions that are associated with WBD resistance, and 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.
We would like to thank Dr. Karina Peres Gramacho (CEPLAC/CEPEC) for preparing the spore inoculation and for her help with the inoculation of the plants for the gene expression experiment. Thanks are due to Dayana Rodezno (Mars, Inc.) for her help with the gene expression analysis and to Dr. J. Conrad Stack (Mars, Inc.) for helpful comments on the manuscript. Thanks are also due to Valdevino Santana do Carmo (MCCS) for his help with the field data collection. We would also like to thank Mars, Inc. and their funding of Trust Agreement # 58-6631-6-123: Genetic Improvement of Theobroma cacao, to the United States Department of Agriculture - Agricultural Research Service, Subtropical Horticulture Research Station (USDA-ARS SHRS) cacao program, which supported, in part, this research. This research was conducted in Brazil under the CGEN permit for access to genetic resources number 02000.000070/2015-05.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Motamayor JC, Lachenaud P, da Silva E Mota JW, Loor R, Kuhn DN, Brown JS, et al. Geographic and genetic population differentiation of the amazonian chocolate tree (Theobroma cacao L.). PLoS One. 2008;3(10):e3311.PubMed CentralView ArticlePubMedGoogle Scholar
- Motamayor JC, Risterucci AM, Lopez PA, Ortiz CF, Moreno A, Lanaud C. Cacao domestication I: the origin of the cacao cultivated by the Mayas. Heredity. 2002;89(5):380–6.View ArticlePubMedGoogle Scholar
- Alverson WS, Whitlock BA, Nyffeler R, Bayer C, Baum DA. Phylogeny of the core Malvales: evidence from ndhF sequence data. Am J Bot. 1999;86(10):1474–86.View ArticlePubMedGoogle Scholar
- Dantas LG, Guerra M. Chromatin differentiation between Theobroma cacao L. and T. grandiflorum Schum. Genet Mol Biol. 2010;33:94–8.PubMed CentralView ArticlePubMedGoogle Scholar
- Argout X, Salse J, Aury J-M, Guiltinan MJ, Droc G, Gouzy J, et al. The genome of Theobroma cacao. Nat Genet. 2011;43(2):101–8.View ArticlePubMedGoogle Scholar
- Motamayor JC, Mockaitis K, Schmutz J, Haiminen N, Livingstone D, Cornejo O, et al. The genome sequence of the most widely cultivated cacao type and its use to identify candidate genes regulating pod color. Genome Biol. 2013;14(6):R53.PubMed CentralView ArticlePubMedGoogle Scholar
- Bowers JH, Bailey BA, Hebbar PK, Sanogo S, Lumsden RD. The impact of plant diseases on world chocolate production. Online. Plant Health Progress. doi:101094/PHP-2001-0709-01-RV 2001.Google Scholar
- Teixeira PJPL, Thomazella DPT, Pereira GAG. Time for chocolate: current understanding and New perspectives on cacao Witches’ broom disease research. PLoS Pathog. 2015;11(10):e1005130.PubMed CentralView ArticlePubMedGoogle Scholar
- Aime MC, Phillips-Mora W. The causal agents of witches’ broom and frosty pod rot of cacao (chocolate, Theobroma cacao) form a new lineage of Marasmiaceae. Mycologia. 2005;97(5):1012–22.View ArticlePubMedGoogle Scholar
- Pound FJ. Cacao and witches’ broom disease (Marasmius perniciosa). Report on a recent visit to the Amazon territory of Peru, September, 1942 - February, 1943. Port of Spain: Government Print; 1943.Google Scholar
- Figueira A, Albuquerque P, Leal-Jr G. Genetic mapping and differential gene expression of Brazilian alternative resistance sources to Witches’ Broom (causal agent Crinipellis perniciosa), 15th Proc Intl Cocoa Res Conf IS. 2006. p. 241–55.Google Scholar
- Brown JS, Schnell II RJ, Motamayor JC, Lopes U, Kuhn DN. Resistance gene mapping for witches’ broom disease in Theobroma cacao L. in an F2 population using SSR markers and candidate genes. J Am Soc Hortic Sci. 2005;130(3):366–73.Google Scholar
- Faleiro FG, Queiroz VT, Lopes UV, Guimarães CT, Pires JL, Yamada MM, et al. Mapping QTLs for Witches’ Broom (Crinipellis Perniciosa) Resistance in Cacao (Theobroma Cacao L.). Euphytica. 2006;149(1–2):227–35.View ArticleGoogle Scholar
- Queiroz VT, Guimarães CT, Anhert D, Schuster I, Daher RT, Pereira MG, et al. Identification of a major QTL in cocoa (Theobroma cacao L.) associated with resistance to witches’ broom disease. Plant Breed. 2003;122(3):268–72.View ArticleGoogle Scholar
- Lanaud C, Risterucci AM, N’Goran AKJ, Clement D, Flament MH, Laurent V, et al. A genetic linkage map of Theobroma cacao L. Theor Appl Genet. 1995;91(6–7):987–93.PubMedGoogle Scholar
- Risterucci AM, Grivet L, N’Goran JAK, Pieretti I, Flament MH, Lanaud C. A high-density linkage map of Theobroma cacao L. Theor Appl Genet. 2000;101(5–6):948–55.View ArticleGoogle Scholar
- Pugh T, Fouet O, Risterucci AM, Brottier P, Abouladze M, Deletrez C, et al. A new cacao linkage map based on codominant markers: development and integration of 201 new microsatellite markers. Theor Appl Genet. 2004;108(6):1151–61.View ArticlePubMedGoogle Scholar
- Fouet O, Allegre M, Argout X, Jeanneau M, Lemainque A, Pavek S, et al. Structural characterization and mapping of functional EST-SSR markers in Theobroma cacao. Tree Genet Genomes. 2011;7(4):799–817.View ArticleGoogle Scholar
- Allegre M, Argout X, Boccara M, Fouet O, Roguet Y, Berard A, et al. Discovery and mapping of a new expressed sequence tag-single nucleotide polymorphism and simple sequence repeat panel for large-scale genetic studies and breeding of Theobroma cacao L. DNA Res. 2012;19(1):23–35.PubMed CentralView ArticlePubMedGoogle Scholar
- Silva DV, Araújo IS, Branco SMJ, Aguilar-Vildoso CI, Lopes UV, Marelli JP, et al. Analysis of resistance to witches’ broom disease (Moniliophthora perniciosa) in flower cushions of Theobroma cacao in a segregating population. Plant Pathol. 2014;63(6):1264–71.View ArticleGoogle Scholar
- Kuhn D, Livingstone III D, Main D, Zheng P, Saski C, Feltus FA, et al. Identification and mapping of conserved ortholog set (COS) II sequences of cacao and their conversion to SNP markers for marker-assisted selection in Theobroma cacao and comparative genomics studies. Tree Genet Genomes. 2012;8(1):97–111.View ArticleGoogle Scholar
- Livingstone D, Royaert S, Stack C, Mockaitis K, May G, Farmer A, Saski C, Schnell R, Kuhn D, Motamayor JC. Making a Chocolate Chip: Development and Evaluation of a 6K SNP Array for Theobroma cacao. DNA research. 2015. in press.Google Scholar
- Van Ooijen JW. JoinMap®4, software for the calculation of genetic linkage maps in experimental populations. Kyazma B.V.: Wageningen; 2006.Google Scholar
- Van Ooijen JW. Multipoint maximum likelihood mapping in a full-sib family of an outbreeding species. Genet Res (Camb). 2011;93(5):343–9.View ArticleGoogle Scholar
- Voorrips RE. MapChart: software for the graphical presentation of linkage maps and QTLs. J Hered. 2002;93(1):77–8.View ArticlePubMedGoogle Scholar
- Payne R. A Guide to ANOVA and Design in GenStat® (15th Edition). 5 The Waterhouse, Waterhouse Street, Hemel Hempstead, Hertfordshire HP1 1ES, UK: VSN International; 2012.Google Scholar
- Malosetti M, Ribaut J, Vargas M, Crossa J, Eeuwijk F. A multi-trait multi-environment QTL mixed model with an application to drought and nitrogen stress trials in maize (Zea mays L.). Euphytica. 2008;161(1–2):241–57.View ArticleGoogle Scholar
- Utro F, Haiminen N, Livingstone D, Cornejo O, Royaert S, Schnell R, et al. iXora: exact haplotype inferencing and trait association. BMC Genet. 2013;14(1):48.PubMed CentralView ArticlePubMedGoogle Scholar
- Scheet P, Stephens M. A fast and flexible statistical model for large-scale population genotype data: applications to inferring missing genotypes and haplotypic phase. Am J Hum Genet. 2006;78(4):629–44.PubMed CentralView ArticlePubMedGoogle Scholar
- Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: Molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30(12):2725–9.PubMed CentralView ArticlePubMedGoogle Scholar
- Clement D, Risterucci AM, Motamayor JC, N’Goran J, Lanaud C. Mapping quantitative trait loci for bean traits and ovule number in Theobroma cacao L. Genome. 2003;46(1):103–11.View ArticlePubMedGoogle Scholar
- Maharaj K, Maharaj P, Bekele FL, Ramnath D, Bidaisee GG, Bekele I, et al. Trinidad selected hybrids: an investigation of the phenotypic and agro-economic traits of 20 selected cacao cultivars. Trop Agric (Trinidad). 2011;88(4):175–85.Google Scholar
- Turnbull CJ, Hadley P. International Cocoa Germplasm Database (ICGD). [Online Database]. CRA Ltd/NYSE Liffe/University of Reading, UK Available: http://wwwicgdreadingacuk (30 September 2013) 2013.Google Scholar
- Bartley BGD. Additional Notes with the translation of Coral’s field notes from R. Ucayali, Peru (1987–88) expedition. B.G.D Bartley. Personal Communications. 1996.Google Scholar
- Lockwood G, Gyamfi MMO. The CRIG cocoa germplasm collection with notes on codes used in the breeding programme at Tafo and elsewhere. Technical bulletin, no. 10. Tafo: Cocoa Research Institute; 1979.Google Scholar
- Boza EJ, Motamayor JC, Amores FM, Cedeño-Amador S, Tondo CL, Livingstone DS, et al. Genetic characterization of the cacao cultivar CCN 51: its impact and significance on global cacao improvement and production. J Am Soc Hortic Sci. 2014;139(2):219–29.Google Scholar
- Enriquez GA, Soria VJ. Cacao cultivars register. Costa Rica: IICA Teaching and Research Centre; 1967.Google Scholar
- Morera J, Mora A. Banco de Germoplasma de Cacao del CATIE. Turrialba: CATIR|E; 1990.Google Scholar
- Mejía LA, Argüello CO. Tecnología para el Mejoramiento del Sistema de Producción de Cacao. Bucaramanga: CORPOICA; 2000.Google Scholar
- Trinidad-CRU. Cocoa research unit report for 1994. Trinidad: The University of the West Indies; 1994.Google Scholar
- Santos RMF, Lopes UV, Bahia RC, Machado RCR, Ahnert D, Corrêa RX. Marcadores microssatélites relacionados com a resistência à vassoura-de-bruxa do cacaueiro. Pesq Agropec Bras. 2007;42:1137–42.View ArticleGoogle Scholar
- Flament MH, Kebe I, Clément D, Pieretti I, Risterucci AM, N’Goran JAK, et al. Genetic mapping of resistance factors to Phytophthora palmivora in cocoa. Genome. 2001;44(1):79–85.View ArticlePubMedGoogle Scholar
- Brown JS, Phillips-Mora W, Power EJ, Krol C, Cervantes-Martinez C, Motamayor JC, et al. Mapping QTLs for resistance to frosty pod and black pod diseases and horticultural traits in Theobroma cacao L. Crop Sci. 2007;47(5):1851.View ArticleGoogle Scholar
- Crouzillat D, Lerceteau E, Petiard V, Morera J, Rodriguez H, Walker D, et al. Theobroma cacao L.: a genetic linkage map and quantitative trait loci analysis. Theor Appl Genet. 1996;93(1):205–14.View ArticlePubMedGoogle Scholar
- Neff MM, Neff JD, Chory J, Pepper AE. dCAPS, a simple technique for the genetic analysis of single nucleotide polymorphisms: experimental applications in Arabidopsis thaliana genetics. Plant J. 1998;14(3):387–92.View ArticlePubMedGoogle Scholar
- Dodds PN, Rathjen JP. Plant immunity: towards an integrated view of plant–pathogen interactions. Nat Rev Genet. 2010;11(8):539–48.View ArticlePubMedGoogle Scholar
- Deepak S, Shailasree S, Kini RK, Muck A, Mithöfer A, Shetty SH. Hydroxyproline-rich glycoproteins and plant defence. J Phytopathol. 2010;158(9):585–93.Google Scholar
- Asai T, Tena G, Plotnikova J, Willmann MR, Chiu W-L, Gomez-Gomez L, et al. MAP kinase signalling cascade in Arabidopsis innate immunity. Nature. 2002;415(6875):977–83.View ArticlePubMedGoogle Scholar
- Robatzek S, Somssich IE. Targets of AtWRKY6 regulation during plant senescence and pathogen defense. Genes Dev. 2002;16(9):1139–49.PubMed CentralView ArticlePubMedGoogle Scholar
- Tena G, Boudsocq M, Sheen J. Protein kinase signaling networks in plant innate immunity. Curr Opin Plant Biol. 2011;14(5):519–29.PubMed CentralView ArticlePubMedGoogle Scholar
- Meng X, Zhang S. MAPK cascades in plant disease resistance signaling. Annu Rev Phytopathol. 2013;51(1):245–66.View ArticlePubMedGoogle Scholar
- Nanda AK, Andrio E, Marino D, Pauly N, Dunand C. Reactive oxygen species during plant-microorganism early interactions. J Integr Plant Biol. 2010;52(2):195–204.View ArticlePubMedGoogle Scholar
- Ceita GO, Macêdo JNA, Santos TB, Alemanno L, da Silva Gesteira A, Micheli F, et al. Involvement of calcium oxalate degradation during programmed cell death in Theobroma cacao tissues triggered by the hemibiotrophic fungus Moniliophthora perniciosa. Plant Sci. 2007;173(2):106–17.View ArticleGoogle Scholar
- Dias CV, Mendes JS, dos Santos AC, Pirovani CP, da Silva GA, Micheli F, et al. Hydrogen peroxide formation in cacao tissues infected by the hemibiotrophic fungus Moniliophthora perniciosa. Plant Physiol Biochem. 2011;49(8):917–22.View ArticlePubMedGoogle Scholar
- Rojas CM, Senthil-Kumar M, Wang K, Ryu C-M, Kaundal A, Mysore KS. Glycolate oxidase modulates reactive oxygen species–mediated signal transduction during nonhost resistance in nicotiana benthamiana and arabidopsis. The Plant Cell Online. 2012;24(1):336–52.View ArticleGoogle Scholar
- Bednarek P, Osbourn A. Plant-microbe interactions: chemical diversity in plant defense. Science. 2009;324(5928):746–8.View ArticlePubMedGoogle Scholar
- Wang H, Wu J, Sun S, Liu B, Cheng F, Sun R, et al. Glucosinolate biosynthetic genes in Brassica rapa. Gene. 2011;487(2):135–42.View ArticlePubMedGoogle Scholar
- Century KS, Holub EB, Staskawicz BJ. NDR1, a locus of Arabidopsis thaliana that is required for disease resistance to both a bacterial and a fungal pathogen. Proc Natl Acad Sci U S A. 1995;92(14):6597–601.PubMed CentralView ArticlePubMedGoogle Scholar
- Smalle J, Vierstra RD. The ubiquitin 26S proteasome proteolytic pathway. Annu Rev Plant Biol. 2004;55:555–90.View ArticlePubMedGoogle Scholar
- Vierstra RD. The ubiquitin/26S proteasome pathway, the complex last chapter in the life of many plant proteins. Trends Plant Sci. 2003;8(3):135–42.View ArticlePubMedGoogle Scholar
- Zeng L-R, Vega-Sanchez ME, Zhu T, Wang G-L. Ubiquitination-mediated protein degradation and modification: an emerging theme in plant-microbe interactions. Cell Res. 2006;16(5):413–26.View ArticlePubMedGoogle Scholar
- Kirsch C, Logemann E, Lippok B, Schmelzer E, Hahlbrock K. A highly specific pathogen-responsive promoter element from the immediate-early activated CMPG1 gene in Petroselinum crispum. Plant J. 2001;26(2):217–27.View ArticlePubMedGoogle Scholar
- Rowland O, Ludwig AA, Merrick CJ, Baillieul F, Tracy FE, Durrant WE, et al. Functional analysis of Avr9/Cf-9 rapidly elicited genes identifies a protein kinase, ACIK1, that is essential for full Cf-9-dependent disease resistance in tomato. Plant Cell. 2005;17(1):295–310.PubMed CentralView ArticlePubMedGoogle Scholar
- Heise A, Lippok B, Kirsch C, Hahlbrock K. Two immediate-early pathogen-responsive members of the AtCMPG gene family in Arabidopsis thaliana and the W-box-containing elicitor-response element of AtCMPG1. Proc Natl Acad Sci U S A. 2002;99(13):9049–54.PubMed CentralView ArticlePubMedGoogle Scholar
- Li W, Ahn IP, Ning Y, Park CH, Zeng L, Whitehill JG, et al. The U-Box/ARM E3 ligase PUB13 regulates cell death, defense, and flowering time in Arabidopsis. Plant Physiol. 2012;159(1):239–50.PubMed CentralView ArticlePubMedGoogle Scholar
- Zeng LR, Park CH, Venu RC, Gough J, Wang GL. Classification, expression pattern, and E3 ligase activity assay of rice U-box-containing proteins. Mol Plant. 2008;1(5):800–15.View ArticlePubMedGoogle Scholar
- Zeng L-R, Qu S, Bordeos A, Yang C, Baraoidan M, Yan H, et al. Spotted leaf11, a negative regulator of plant cell death and defense, encodes a U-Box/Armadillo repeat protein endowed with E3 ubiquitin ligase activity. The Plant Cell Online. 2004;16(10):2795–808.View ArticleGoogle Scholar
- Nagase T, Kikuno R, Nakayama M, Hirosawa M, Ohara O. Prediction of the coding sequences of unidentified human genes. XVIII. The complete sequences of 100 new cDNA clones from brain which code for large proteins in vitro. DNA Res. 2000;7(4):273–81.View ArticlePubMedGoogle Scholar
- Hahn A, Harter K. Mitogen-activated protein kinase cascades and ethylene: signaling, biosynthesis, or both? Plant Physiol. 2009;149(3):1207–10.PubMed CentralView ArticlePubMedGoogle Scholar
- Kieber JJ, Rothenberg M, Roman G, Feldmann KA, Ecker JR. CTR1, a negative regulator of the ethylene response pathway in Arabidopsis, encodes a member of the raf family of protein kinases. Cell. 1993;72(3):427–41.View ArticlePubMedGoogle Scholar
- Frye CA, Tang D, Innes RW. Negative regulation of defense responses in plants by a conserved MAPKK kinase. Proc Natl Acad Sci U S A. 2001;98(1):373–8.PubMed CentralView ArticlePubMedGoogle Scholar
- Franceschini A, Szklarczyk D, Frankild S, Kuhn M, Simonovic M, Roth A, et al. STRING v9.1: protein-protein interaction networks, with increased coverage and integration. Nucleic Acids Res. 2013;41(D1):D808–15.PubMed CentralView ArticlePubMedGoogle Scholar
- Shunwu Y, Zhang L, Zuo K, Li Z, Tang K. Isolation and characterization of a BURP domain-containing gene BnBDC1 from Brassica napus involved in abiotic and biotic stress. Physiol Plant. 2004;122(2):210–8.View ArticleGoogle Scholar
- Coleman HD, Canam T, Kang K-Y, Ellis DD, Mansfield SD. Over-expression of UDP-glucose pyrophosphorylase in hybrid poplar affects carbon allocation. J Exp Bot. 2007;58(15–16):4257–68.View ArticlePubMedGoogle Scholar
- Feys BJ, Wiermer M, Bhat RA, Moisan LJ, Medina-Escobar N, Neu C, et al. Arabidopsis SENESCENCE-ASSOCIATED GENE101 stabilizes and signals within an ENHANCED DISEASE SUSCEPTIBILITY1 complex in plant innate immunity. Plant Cell. 2005;17(9):2601–13.PubMed CentralView ArticlePubMedGoogle Scholar
- Sena K, Alemanno L, Gramacho KP. The infection process of Moniliophthora perniciosa in cacao. Plant Pathol. 2014;63(6):1272–81.View ArticleGoogle Scholar
- Teixeira PJPL, Thomazella DPdT, Reis O, do Prado PFV, do Rio MCS, Fiorin GL, José J, Costa GGL, Negri VA, Mondego JMC et al. High-Resolution Transcript Profiling of the Atypical Biotrophic Interaction between Theobroma cacao and the Fungal Pathogen Moniliophthora perniciosa. The Plant Cell Online. 2014;26(11):4245-4269.Google Scholar