Original Research ARTICLE
Uncovering the genetic architecture of Colletotrichum lindemuthianum resistance through QTL mapping and epistatic interaction analysis in common bean
- 1Grupo de Biología de Agrosistemas, Misión Biológica de Galicia-CSIC, Pontevedra, Spain
- 2Departamento de Biología y Geología (Genética), Centro de Investigación en Biotecnología Agroalimentaria (BITAL), Universidad de Almería, Almería, Spain
Colletotrichum lindemuthianum is a hemibiotrophic fungal pathogen that causes anthracnose disease in common bean. Despite the genetics of anthracnose resistance has been studied for a long time, few quantitative trait loci (QTLs) studies have been conducted on this species. The present work examines the genetic basis of quantitative resistance to races 23 and 1545 of C. lindemuthianum in different organs (stem, leaf and petiole). A population of 185 recombinant inbred lines (RIL) derived from the cross PMB0225 × PHA1037 was evaluated for anthracnose resistance under natural and artificial photoperiod growth conditions. Using multi-environment QTL mapping approach, 10 and 16 main effect QTLs were identified for resistance to anthracnose races 23 and 1545, respectively. The homologous genomic regions corresponding to 17 of the 26 main effect QTLs detected were positive for the presence of resistance-associated gene cluster encoding nucleotide-binding and leucine-rich repeat (NL) proteins. Among them, it is worth noting that the main effect QTLs detected on linkage group 05 for resistance to race 1545 in stem, petiole and leaf were located within a 1.2 Mb region. The NL gene Phvul.005G117900 is located in this region, which can be considered an important candidate gene for the non-organ-specific QTL identified here. Furthermore, a total of 39 epistatic QTL (E-QTLs) (21 for resistance to race 23 and 18 for resistance to race 1545) involved in 20 epistatic interactions (eleven and nine interactions for resistance to races 23 and 1545, respectively) were identified. None of the main and epistatic QTLs detected displayed significant environment interaction effects. The present research provides essential information not only for the better understanding of the plant-pathogen interaction but also for the application of genomic assisted breeding for anthracnose resistance improvement in common bean through application of marker-assisted selection (MAS).
Members of the ascomycete genus Colletotrichum cause devastating anthracnose diseases in many agronomically important crops in temperate, tropical and subtropical regions (Bailey and Jeger, 1992). The specialized hemibiotrophic fungus C. lindemuthianum [(Sacc. & Magnus) Lams. - Scrib]. has a reduced number of plant hosts, mainly common bean (Phaseolus vulgaris L.), although in less extent and severity, it can also colonize P. acutifolius var. lactifolius, P. coccineus, P. aureus, P. lunatus, P. limensis, Medicago sativa, and Vicia faba (Sicard et al., 1997; Mahuku et al., 2002). The pathogen has a sequential biotrophic- and necrotrophic-infection process to invade and colonize the plant hosts, that involves the transition from an asymptomatic biotrophic phase (characterized by intracellular thick primary hyphae) to a destructive necrotrophic phase (characterized by thin filamentous secondary hyphae) referred to as the biotrophynecrotrophy switch, which is essential for anthracnose disease development (Bhadauria et al., 2011). The remarkable resistance of C. lindemuthianum and its capacity for survivability in any environmental condition renders its presence responsible for losses in crops. In fact, the damage caused by this fungus in bean crops is so great that it has produced an economical loss in productive countries (Vigidal-Filho et al., 2007). Besides, C. lindemuthianum causes a hypersensitive response in bean resistant plants - groups of red-brownish wounds of different sizes that are produced by the plant to delimit the spread of the pathogenic fungus (Martínez-Pacheco et al., 2009). The process of co-evolution between the fungus and bean resistant plants has led this fungal species to produce new pathogenic variants, which can be detected on the basis of the phenotypic response to anthracnose infection shown by different varieties of common bean (Melotto et al., 2000; Rodríguez-Guerra et al., 2003). Thus, more than 100 races have been described for C. lindemuthianum (Rodríguez-Guerra et al., 2003) and new pathotypes are reported every day, indicating a large pathogenic variability of this fungus.
In common bean, up to 40 genes conferring resistance to specific races (designated as Co-) have been described, mainly due to C. lindemuthianum pathogenic variability. Anthracnose resistance is related to the presence of closely linked race-specific loci, which comprise different single, duplicate or complementary dominant genes, except for the recessive co-8 gene (Kelly and Vallejo, 2004; Ferreira et al., 2013; Campa et al., 2014). Based on the hypothesis that the same gene confers the resistance to different races in a bean genotype, most classical studies considered that different resistance spectra in genotypes were due to different alleles of the same gene. As a result, different alleles were described for genes Co-1, Co-3, and Co-4 (Kelly and Vallejo, 2004; Ferreira et al., 2013). Most identified resistance genes have been mapped on the different linkage groups (LG) of the common bean genetic map: genes Co-1, Co-12, Co-13, Co-14, Co-15, Co-165-X, Co-173-X, Co-x, and Co-w were mapped on LG01 (Barrus, 1915; Melotto and Kelly, 2000; Méndez-Vigo, 2001; Gonçalves-Vidigal and Kelly, 2006; Rodríguez-Suárez et al., 2007; Geffroy et al., 2008; Vallejo and Kelly, 2008; Gonçalves-Vidigal et al., 2011; Campa et al., 2014); CoPv02c3-X, CoPv02c7-X, CoPv02c19-X, CoPv0c2449-X, and Co-u on LG02 (Kelly et al., 2003; Campa et al., 2014); Co-13 on LG03 (Lacanallo et al., 2010); Co-3, Co-3c3-X, Co-3c7-X, Co-3c19-X, Co-3c449-X, Co3c453-X, Co-9, Co-y, Co-z, Co-10, and Co-15 on LG04 (Geffroy et al., 1999; Alzate-Marín et al., 2003; Méndez-Vigo et al., 2005; Rodríguez-Suárez et al., 2007, 2008; Gonçalves-Vidigal et al., 2013; Sousa et al., 2013; Campa et al., 2014); Co-5, Co-6, and Co-v on LG07 (Fouilloux, 1976; Campa et al., 2009); Co-4 on LG08 (Melotto et al., 2004; Rodríguez-Suárez et al., 2007; Campa et al., 2014); CoPv09c453-C on LG09 (Campa et al., 2014); and Co-2, Co-26-C, Co-239-C, Co-238-C, Co-2357-C, Co-265-C, Co-27-C, Co-23-C, Co-219-C, Co-2449-C, and Co-2453-C on LG11 (Adam-Blondon et al., 1994; Campa et al., 2014). Although genetic analyses support that the Co-genes behave as major Mendelian factors, they most likely exist as resistance gene clusters in which individual gene(s) confers resistance to one specific race. Most of the resistance-associated genes encode nucleotide-binding and leucine-rich repeat proteins, which are known as NB-LRR (NL) genes (Meyers et al., 2005). The presence of these clusters is widespread among higher plants, and clusters of NL genes have also been described in the common bean genome (Schmutz et al., 2014). In particular, two large clusters identified on chromosomes 4 and 11 could co-localize with previously mapped Co-3 and Co-2 anthracnose resistance genes, respectively.
Resistance to anthracnose in common bean generally follows a qualitative mode of inheritance where resistant and susceptible reactions are clearly differentiated. The specific resistance genes follow the classic gene-for-gene model (Flor, 1955), and the qualitative resistance provided by them is often less durable than quantitative resistance, since pathogens can more easily adapt to single gene-mediated resistance (St. Clair, 2010). In contrast, quantitative resistance usually confers broad-spectrum protection toward different races of biotrophic or necrotrophic pathogens (Oliver and Ipcho, 2004). The genetic regulation of quantitative traits is often complex due to their polygenic nature. However, trait dissection through Quantitative Trait Loci (QTL) analysis is a useful approach to identify chromosomal regions harboring genes that control these quantitative traits. Yet, in addition to mapping main effect QTLs, epistatic interactions between QTLs are important. Identification of quantitative disease resistance main and epistatic effects from multiple environments does not only help to extend the applicability of results, but is also essential for the development of an efficient marker-assisted selection (MAS) program aimed at improving breeding efficiency.
Despite the fact that genetics of anthracnose resistance in common bean has been studied for a long time, few QTL studies have been conducted on this species. The present work studies the genetic basis of quantitative resistance to two races of C. lindemuthianum in different organs of a segregating common bean recombinant inbred line population (RIL) from the cross PMB0225 × PHA1037. Using multi-environment QTL mapping approach, race specific anthracnose resistance QTLs were identified showing significant main additive effects in stem, petiole and leaf organs, which were co-localized with NL genes. In addition to identifying main effect QTLs, this analysis revealed epistatic interactions that explained phenotypic variation beyond those controlled by main effects of individual loci. Thus, markers associated with QTLs reported here constitute useful tools for MAS breeding programs directed toward improved anthracnose resistance.
Materials and Methods
A RIL population consisting of 185 F7 lines was developed by single-seed descent from an F2 population from the cross between PMB0225 (a common bean line as P1) and PHA1037 (nuña bean line abbreviated as P2) accessions belonging to the Andean gene pool. Mesoamerican (17, 73, 448, and 1545) and Andean (7, 23, 39, 55, and 102) anthracnose races were inoculated to the two parents. The twelve differential cultivars (Michelite, MDRK, Perry Marrow, Cornell 49242, Widusa, Kaboon, Mexico 222, PI207262, TO, TU, AB136, and G2333) were used to confirm the identity of the C. lindemuthianum races. Only races 23 and 1545 were pathogenic on PMB0225 parent and chosen for the present study.
Plant Growth Conditions, Inoculation, and Disease Evaluation
Plants were grown in plastic pots containing a mixture of clay soil and organic compound (1:1; v/v), under natural and artificial (12-h photoperiod, 166 μE s−1 m−2) photoperiod growth conditions with average day and night temperatures of 25 and 20°C, respectively. Plants were irrigated according to water needs. The anthracnose races were kept on potato-dextrose agar (PDA) at 19–21°C in darkness. To obtain conidia, fungus was grown for sporulation for about 15 days, and medium plates were flooded with 10 mL of 0.01% Tween 80 in distilled water. The conidial suspension was collected and filtered twice to remove mycelial fragments. The number of conidia was estimated using a haemocytometer and inoculum concentration was adjusted to 3 × 106 conidia mL−1 with distilled water. Spore suspension was sprayed-inoculated onto 2-week-old bean plants showing the fully expanded primary leaves using an atomizer. Inoculated plants were sealed in order to increase 95–100% humidity for 48 h.
The infected phenotypes were assessed on the basis of symptom severity on the primary leaves (L), stems (S), and petioles (P) at intervals of 7, 14, and 21 days post-inoculation (dpi). Numerical disease scores (DC) were assigned based on visual appreciation of the percentage of the organ presenting symptoms. A score of 1 represented no observed symptoms, while 9 corresponded to 100% of the organ covered by brown typical lesions of anthracnose (Figure 1). The Area Under the Disease Progress Curve (AUDPC) was calculated according to Shaner and Finney (1977) as: AUDPC = , where xi is the disease score on date i, n the number of evaluations made and the time in days between evaluations xiand xi + 1. The use of AUDPC is an effective method to take both duration and severity of disease into account.
Experimental Design and Statistical Data Analysis
The experiment was set up as a randomized complete block design with four replicates in artificial (named A-Ppd) and natural (named N-Ppd) photoperiod conditions, respectively. Each RIL genotype was represented by one plant in each block. Independent four-block experiments were carried out for each race, and the parental lines PMB0225 and PHA1037 were included.
Descriptive statistical parameters (mean value, standard deviation and range of variation) and normality (Kolmogorov–Smirnov test) were obtained for each quantitative trait and environment. Variation in the expression of traits through the environments was analyzed using PROC MIXED (SAS Institute Inc V. 9.04, Cary, NC, USA). Variance components and broad-sense heritabilities with their standard errors were estimated by restricted maximum likelihood (REML) option of the PROC MIXED and IML (SAS Institute Inc. v. 9.04, Cary, NC, USA) for the phenotypic traits (Holland et al., 2003; Holland, 2006). Phenotypic Pearson correlation coefficients among traits were implemented using PROC CORR across the environments (SAS Institute Inc. v. 9.04, Cary, NC, USA).
The genetic linkage map described by Yuste-Lisbona et al. (2012) was used for QTL analysis. The SCAR SW13 and SW12 (Fourie et al., 2004; Rodríguez-Suárez et al., 2008), 1 AFLP, 3 SSR, 29 SNP, and the seed coat color gene (P) were added to this map, which finally consisted of 229 loci (86 AFLP, 98 SSR, 42 SNP, 2 SCAR, and P locus) distributed on 11 LGs. The map spanned 858.4 cM, with an average distance of 3.7 cM between adjacent markers. Marker data were analyzed by JoinMap® 4.0 software (Van Ooijen, 2006). A minimum logarithm of odds ratio (LOD) score of 6.0 and a recombination frequency value of 0.3 were set as the linkage threshold for grouping markers. The Kosambi map function (Kosambi, 1944) was used to calculate the genetic distance between markers. The LGs were designated according to Pedrosa-Harand et al. (2008). QTLNetwork 2.0 software (Yang et al., 2008) was used to identify single-locus QTL, epistatic QTL (E-QTL) and their environment interaction effects (QTL × Environment, QE; and E-QTL × Environment, E-QE). The mixed-model based on composite interval mapping method (MCIM) was carried out for one-dimensional genome scan to detect putative single-locus QTL (defined as those showing significant main additive effects) and their environment interactions. In addition, a two-dimensional genome scan was carried out to identify epistatic interaction effects. An experimental-wise significance level of 0.05 was designated for candidate interval selection, putative QTL detection, and QTL effect. Both testing and filtration window size were set at 10 cM, with a walk speed of 1 cM. The critical F-value to declare putative QTLs was determined by a 1000 permutation test at 95% confidence level. The effects of QTL and environment interactions were estimated by the Markov Chain Monte Carlo method (Wang et al., 1994). QTL with only genetic effects indicated that these were expressed in the same way across environments. In addition, QTL with environment interaction effects suggested that their expressions were environmentally dependent. The detected QTLs were designated as recommended by Miklas and Porch (2010). The genetic map and the QTL detected were drawn using the MapChart 2.2 software (Voorrips, 2002).
Identifying Location of QTL in Common Bean Genome
Nucleotide sequences of the markers flanking the main effect QTLs were used as queries for BLASTN search (Altschul et al., 1997) against the first chromosome scale version of P. vulgaris genome (Schmutz et al., 2014) available in the Phytozome database (http://www.phytozome.net/).
Resistance Variation in the RIL Population
The bean accession PMB0225 was fully susceptible to anthracnose infection to race 1545 in all tested organs, and displayed susceptibility in leaf, intermediate resistance in stem, and full resistance in petiole to race 23. The PHA-1037 accession was fully resistant to both races in all organs. Table 1 shows the mean values and standard errors of the parental genotypes and the RIL population, as well as the ranges of variation of the RIL population for the resistance traits for each environment. In the RIL population, a continuous but bimodal distribution skewed toward the resistant parent PHA1037 was found regardless of the organ and race tested (Figure S1). The relative skewedness toward the resistant PHA1037 parent would imply that multiple genes with complementary additive effects are conferring resistance to anthracnose. Variance analysis was conducted for each environment and difference between blocks was not significant for most of the environments and resistant traits (Table 1). PMB0225 and PHA1037 parents and RIL progeny were significantly different for resistance traits in each environment (P = 0.001), demonstrating a genetic origin for the different levels of resistance in the RIL population.
Table 1. Estimates of means, standard errors, range of variation, and variance analysis results for anthracnose resistance to races 23 and 1545 of the two common bean parents, PMB0225 and PHA1037, and the RIL population, grown in two environments (Env).
The estimated broad-sense heritability estimates for resistance traits between organs for a given race were high, with values ≥0.70 (Table 2). These values are in accordance with those reported by Geffroy et al. (2000) for anthracnose resistance between aerial organs of the plant to isolates 45 and A7, with heritability values ranging from 0.90 to 0.98. There was a strong correlation (P ≤ 0.001) for disease resistance scores between the different organs for a given race (Table 2). This was in sharp contrast to the absence of correlation when the data regarding both races were compared. These findings suggest that different genes could be involved in resistance for a given race, while identical genes condition resistance against the same race in different organs. In addition and in order to determine if the same set of genes encode resistance against both races, co-segregations for resistance to races 23 and 1545 were considered. Twenty-nine and thirty-four RILs were resistant to each race 23 and 1545, respectively, while evidence of co-segregation of disease response to both races was observed in 114 RILs (14 and 74 RILs were susceptible and resistant to both races in all organs, and 26 RILs presented resistance to both races but not in all organs). Accordingly, these results suggest that different genes determined specific resistance to races 23 and 1545.
Table 2. Phenotypic correlation coefficients and heritability estimates with their standard errors among resistance traits for anthracnose races 23 and 1545.
Mapping of Main Effect QTL
The evaluation of the RIL population developed from the cross PMB0225 x PHA1037 under different environments has led to the identification of 10 and 16 main effect QTLs involved in resistance against anthracnose races 23 and 1545, respectively. These QTLs were mapped on eight LGs, with the exception of LGs 02, 10, and 11 (Figure 2). All QTLs detected had significant additive effects and did not display significant additive-by-environment interaction effects. A complete report of the single-locus QTL detected for anthracnose resistance traits is given in Table 3.
Figure 2. Location of main effect QTLs for resistance to anthracnose races 23 and 1545 on a genetic linkage map of common bean based on the RIL population developed from the cross PMB0225 × PHA1037. Distances among markers are indicated in cM to the left of the linkage groups (LG); names of markers are shown on the right. QTLs are depicted as vertical bars to the right of the LG. QTL detected for resistance to race 23 are indicated in black. QTLs identified for resistance to race 1545 are shown in gray. Co-location of resistance-associated genes encoding NB-LRR (NL) proteins with anthracnose QTLs are represented to the right of the QTLs.
Ten main effect QTLs were identified for resistance to race 23: one on each LGs 01 and 03, four on LG04, and two on each LGs 05 and 09 (Figure 2). Three of them had significant effects on stem resistance (SDC and SAUDPC traits), positioned on LGs 01 and 04, and explaining a phenotypic variance from 2.64 to 7.90%. The remaining QTLs were involved in leaf resistance (LDC and LAUDPC traits) on LGs 03, 04, 05, and 09, explaining a phenotypic variance from 3.06 to 11.43%. Two QTLs for resistance in leaf, LDC23-9 and LAUDPC23-9, were co-localized on LG09, explaining 11.43 and 10.19% of the phenotypic variation, respectively. Likewise, the QTLs LDC23-5.1 and LAUDPC23-5.1 for resistance in leaf were co-localized on LG05, explaining 5.98 and 4.36% of the phenotypic variation, respectively. On LG04, four QTLs (SDC23-4, SAUDPC23-4, LDC23-4.1, and LAUDPC23-4) for resistance in stem and leaf were co-localized or nearly co-localized, indicating that this genomic region could confer a non-organ-specific resistance to anthracnose race 23. The total phenotypic variation explained by the main effect QTLs detected varied from 2.64 to 10.93% for stem traits (SAUDPC and SDC, respectively) and from 17.61 to 25.16% for leaf traits (LAUDPC and LDC, respectively). All these QTLs had negative additive values, which indicate that the increase in resistance is due to the presence of the alleles from PHA1037.
For resistance to race 1545, 16 main effect QTLs were found: two on each LG 01, 03, and 07, six on LG05, and four on LG08 (Figure 2). Five of them had significant effects on stem resistance (SDC and SAUDPC traits), positioned on LGs 03, 05, and 08, and explaining a phenotypic variance from 4.33 to 14.24%. For petiole resistance (PDC and PAUDPC traits), five QTLs were detected on LGs 01, 03, and 05, which explain a phenotypic variance from 4.08 to 10.22%. The remaining six QTLs were involved in leaf resistance (LDC and LAUDPC traits) on LGs 05, 07, and 08, with phenotypic variance explained ranging from 1.49 to 16.75%. Organ-specific QTLs were identified for petiole resistance on LG01 (PDC1545-1 and PAUDPC1545-1) and for leaf resistance on LG07 (LDC1545-7 and LAUDPC1545-7). Six main effect QTLs (SDC1545-5, SAUDPC1545-5, PDC1545-5, PAUDPC1545-5, LDC1545-5, and LAUDPC1545-5) were co-localized or nearly co-localized on LG05 for resistance to the three evaluated organs, showing the existence of non-organ-specific resistance to race 1545 in this genomic region. In addition, non-organ-specific QTLs were also found for stem and leaf on LG08 (SDC1545-8, SAUDPC1545-8, LDC1545-8, and LAUDPC1545-8), and for stem and petiole on LG03 (SDC1545-3.1 and PAUDPC1545-3). The total phenotypic variation explained by the additive effects of all putative QTLs identified varied from 17.65 to 23.89% for stem traits (SAUDPC and SDC, respectively), from 15.81 to 17.39% for petiole traits (PDC and PAUDPC, respectively), and from 22.20 to 22.22% for leaf traits (LAUDPC and LDC, respectively). PHA1037 alleles were associated with resistance for most of the QTLs detected, with the exception of the QTLs SDC1545-3.1, PDC1545-1, PAUDPC1545-1, and PAUDPC1545-3, which had positive additive values, indicating that resistance alleles came from PMB0225.
The location of the anthracnose resistance main effect QTLs to races 23 and 1545 was different in most cases, which is in agreement with the absence of correlation between the resistance reactions against both races. Only the QTLs LDC23-3, SDC1545-3.1, and PAUDPC1545-3 were co-localized or nearly co-localized on LG03 for resistance to race 23 in leaf, and to race 1545 in stem and petiole, although with opposite additive values (Table 3).
Detection of Epistatic QTL
A total of 39 E-QTLs (twenty-one for resistance to race 23 and 18 for resistance to race 1545) involved in 20 epistatic interactions (eleven and nine interactions for resistance to races 23 and 1545, respectively) were detected by the combined analysis of the multi-environment phenotypic values. None of the epistatic interactions detected displayed significant environment interaction effects. These E-QTLs were mapped on seven LGs, with the exception of LGs 02, 07, 10, and 11. The positive and negative additive-by-additive epistatic effect values obtained for these epistatic interactions indicate that both parent alleles could contribute to increasing the resistance to anthracnose races 23 and 1545. A complete description of digenic epistatic interaction analysis for anthracnose resistance traits toward both races is shown in Table 4.
Table 4. Epistatic QTLs (E-QTLs) and E-QTL × Environment (E-QE) interaction effects for resistance to anthracnose races 23 and 1545.
For resistance to race 23, 6 of the 21 E-QTLs identified were previously detected as main effect QTLs. Thus, not only did the E-QTLs E-SDC23-1, E-SDC23-4, E-LDC23-3, E-LAUDPC23-4, E-LAUDPC23-5.1, and E-LAUDPC23-9 participate in epistatic interactions, but they also had an individual genetic effect. The analysis revealed novel loci on LGs 01, 03, 05, and 08 interacting so as to influence resistance to race 23. The percentage of phenotypic variance explained by the interaction of the E-QTLs ranged from 1.13 to 7.92%. Among the E-QTLs detected, it is noteworthy that the genomic region located between markers BMc316 and E45M50-69 on LG08 bears E-QTLs (E-SAUDPC23-8, E-LDC23-8, and E-LAUDPC23-8) involved in epistatic interactions for resistance in stem and leaf, which indicates that this region could participate in non-organ-specific resistance to anthracnose race 23. Collectively, the percentage of phenotypic variance explained by the interaction of the E-QTLs varied from 6.85 to 16.06% for stem traits (SDC and SAUDPC, respectively), and from 7.92 to 15.29% for leaf traits (LAUDPC and LDC, respectively).
Novel loci on LGs 01, 03, 04, 06, 08, and 09 were detected as involved in resistance to race 1545. Among the 18 E-QTLs identified, two E-QTLs were previously identified as main effect QTLs (PAUDPC1545-1 and PAUDPC1545-5). The percentage of phenotypic variance explained by the interaction of the E-QTLs varied from 2.25 to 7.05%. All the E-QTLs detected for resistance to race 1545 were organ-specific. The total phenotypic variation explained by the additive-by-additive epistatic effects of all E-QTLs detected ranged from 6.38 to 7.05% for stem traits (SAUDPC and SDC, respectively), from 3.34 to 7.46% for petiole traits (PDC and PAUDPC, respectively), and from 9.45 to 9.89% for leaf traits (LAUDPC and LDC, respectively).
Identifiying Location of QTL in Common Bean Genome
BLASTN analysis of the nucleotide sequences of the markers flanking the main effect QTLs with common bean genome revealed that the homologous regions spanning 17 of the 26 QTLs identified were positive for the presence of NL genes (Figure 2). Thus, the main effect QTLs SDC23-1, PDC1545-1, and PAUDPC1545-1 covered 19.32 cM (30.57–49.89 cM) on LG01, while the corresponding genomic region covered 5.2 Mb on Chromosome (Chr) 1. Within this region, there is a cluster consisting of 17 NL genes. Likewise, the QTLs SDC23-4, SAUDPC23-4, LDC23-4, and LAUDPC23-4 detected on LG04 (13.54–30.95 cM) were located within an important cluster of 41 NL genes on Chr4 (0.4–12.6 Mb). Meanwhile, the QTLs SDC1545-5, SAUDPC1545-5, PDC1545-5, PAUDPC1545-5, LDC1545-5, and LAUDPC1545-5 covered 4.22 cM (45.77–49.99 cM) on LG05, whereas the homologous genomic regions spanned 1.2 Mb on Chr5 (32.7–33.9 Mb). The NL gene Phvul.005G117900 is located in this region, which encodes for a leucine-rich repeat (LRR) protein. In addition, the QTLs SDC1545-8, SAUDPC1545-8, LDC1545-8, and LAUDPC1545-8 identified on LG08 (13.63–30.65 cM) were placed in a cluster of 17 NL genes on Chr8 (1.3–6.6 Mb).
To gain insight in basic knowledge dealing with resistance, a pathosystem that involves C. Lindemuthianum and P. vulgaris model legume has been characterized. The gene action governing anthracnose resistance was studied in a broad set of RILs generated from a cross between susceptible and resistant Andean accessions. Thus, insights into the number of quantitative resistance loci involved in anthracnose resistance to races 23 and 1545 were provided, as well as their epistatic interactions. The pathogen infected and colonized PMB0225 line, although symptoms and pathogen development were significantly reduced in race 23 as compared to race 1545. Genetic and molecular analysis revealed different features associated with the resistance of PHA1037.
Genetic Architecture of Resistance
The phenotypic dissection of anthracnose resistance carried out in the RIL population has led to the detection of different kinds of resistance components. The pathogen resistance response was consistent between the testing environments, which evidenced that anthracnose resistance is mostly influenced by genes rather than environmental conditions. The occurrence of major resistance factors was found across both races and in the different organs tested. Moreover, genomic regions controlling anthracnose resistance displayed additive main effects, epistatic effects or both. This architecture has been frequently reported for other quantitative resistances (Young, 1996). Hence, in addition to main effect QTLs, significant epistatic interactions between QTLs have previously been reported in quantitative resistance against other fungus, such as Phytophthora capsici in pepper (Lefebvre and Palloix, 1996; Thabuis et al., 2003), Rhizoctonia solani in rice (Liu et al., 2014) or Puccinia triticina in wheat (Singh et al., 2014). However, genetic mapping studies considering epistatic interaction effects have not been performed so far in common bean. In this work, depending on the race and organ tested, the total phenotypic variation explained by main effect QTLs ranged from 2.64 to 23.89% (SAUDPC and SDC traits, respectively), whereas epistatic interactions explained a total phenotypic variation from 3.34 to 15.29% (PDC and LDC traits, respectively). Most of the epistatic interactions detected were due to loci without detectable QTL additive main effects, which show the importance of the epistatic effects in genetic resistance to anthracnose. Furthermore, resistant alleles came from the resistant parent PHA1037 more frequently, but they occasionally originated from the susceptible parent PMB0225, as observed in petiole resistance to race 1545. This result suggests that the susceptible parent also develops defense mechanisms, even though their activity could be insufficient to stop fungal progression.
The dissection of resistance into distinct phenotypic resistance components allows for a more precise QTL detection and facilitates the exhaustive selection of resistance factors in breeding programs. This type of approach has led to the identification of organ-specific defense mechanisms for resistance to Phytophthora infestans in potato (Gao et al., 2013) and Colletotrichum graminicola in maize (Balmer et al., 2013). This study has identified QTLs located in the same genomic region for resistance to different organs (SDC23-4, SAUDPC23-4, LDC23-4.1, LAUDPC23-4, SDC1545-5, SAUDPC1545-5, PDC1545-5, PAUDPC1545-5, LDC1545-5, LAUDPC1545-5, SDC1545-8, SAUDPC1545-8, LDC1545-8, and LAUDPC1545-8) or for different races of infection (LDC23-3, SDC1545-3.1, and PAUDPC1545-3), which are usually described as “generalist QTLs” (Lefebvre and Palloix, 1996; Thabuis et al., 2003). In addition, “specialist QTLs” have also been identified, which were involved in organ- or race-specific resistance (SDC23-1, PDC1545-1, PAUDPC1545-1, LDC23-5.1, LAUDPC23-5.1, LDC1545-7, LAUDPC1545-7, LDC23-9 and LAUDPC23-9). However, it is not possible to conclude whether those genomic regions containing “generalist QTLs” resulted from the clustering of “specialist QTLs” or from the pleiotropic effect of a single gene. The presence gene clusters acting on the same trait is widespread among higher plants, and it has also been described in the common bean genome (Schmutz et al., 2014). Geffroy et al. (1999) pointed out that the origin of these clusters preceded the geographic separation of the wild common bean gene pools, as well as the role of selection in the emergence of such clusters, which might confer a selective advantage to the genotype that possesses them.
Co-Localization of QTL with Resistance Genes
The association between NL genes and QTLs conferring resistance to Colletotrichum species has been reported in several plant species. In maize, a gene conferring resistance to anthracnose stalk rot, caused by C. graminicola, encoded a CC-NB-LRR protein (Abad et al., 2006). The physical mapping of RCT1, a host resistance gene against C. trifolii in M. truncatula showed that RCT1 was part of a complex locus containing numerous genes homologous to previously characterized TIR-NB-LRR resistance genes (Yang et al., 2007). In common bean, QTLs associated with anthracnose resistance were mapped in a cluster on the LG04, which was composed by CC-NB-LRR genes (Ferrier-Cana et al., 2003; Geffroy et al., 2009). Therefore, there is strong evidence that NB-LRR genes confer gene-for-gene resistance to Colletotrichum species in diverse plant hosts. Schmutz et al. (2014) stated that the majority of NL genes were physically organized in complex clusters in the common bean genome. In the present study, the homologous regions spanning 17 of the 26 main effect QTLs detected were positive for the presence of NL genes. The main effect QTLs detected on LG01 were co-localized with a cluster of 17 NL genes at the bottom of Chr1. The Co-1 anthracnose resistance cluster is also located in this position of the Chr1, which includes Co-1, Co-12, Co-13 (Melotto and Kelly, 2000), Co-14 (Gonçalves-Vidigal et al., 2011), Co-15 (Gonçalves-Vidigal and Kelly, 2006), Co-165-X, and Co-173-X (Campa et al., 2014). Regarding the linked genes Co-w and Co-x (Geffroy et al., 2008), Richard et al. (2014) have recently positioned Co-x at the end of Chr1 to a 58 kb region that comprises eight genes: three phosphoinositide-specific phospholipases C (PI-PlC), one zinc finger protein, and four kinases, which suggests that Co-x is not a classical NL gene. Moreover, genes for resistance to angular leaf spot, common bacterial blight, Fusarium root rot, and white mold have been located at the bottom of Chr1 (Miklas and Singh, 2007), as well as Phg-1 and Ur-9 genes, which confer resistance against Pseudocercospora griseola and Uromyces appendiculatus, respectively (Kelly and Vallejo, 2004; Gonçalves-Vidigal et al., 2011).
Resistance main effect QTLs to race 23 acting in stem and leaf organs were positioned in an important cluster of 41 NL genes on Chr4. Geffroy et al. (1999) stated that three specific resistant genes for anthracnose were clustered in this region, which originated either from the Mesoamerican BAT93 parent (Co-9) or the Andean Jalo-EEP558 parent (Co-x, Co-y). Furthermore, a major-effect QTL for resistance to isolate 45 (for leaf, stem, and petiole resistance) and a reverse-effect QTL (for leaf resistance) for resistance to isolate A7 (Geffroy et al., 2000), as well as the Co-3 anthracnose resistance cluster (Ferreira et al., 2013; Campa et al., 2014) were also located in this genomic region of Chr4. On the other hand, the main effect QTLs detected at the top of LG08 were placed within a cluster of 17 NL, where the Co-4 anthracnose resistance cluster is located (Melotto et al., 2004; Rodríguez-Suárez et al., 2007; Campa et al., 2014).
In addition, it is worth noting that the main effect QTLs detected on LG05 for resistance to race 1545 in stem, petiole and leaf were positioned within a 1.2 Mb region where the LRR gene Phvul.005G117900 is located. Based on genome sequence analysis, it can be considered an important candidate gene for the non-organ-specific QTL identified here. However, given that regions containing fast evolving genes, such NL genes, that are susceptible to chromosomal rearrangement and transposition or genomic duplication (Meyers et al., 2005), it is not possible to determine if the non-organ-specific resistance resulted from the pleiotropic effect of the Phvul.005G117900 gene or from the clustering of different genes, which could not be present in the reference genome sequence. Thereby, further studies on fine mapping of the target genomic regions would be necessary to draw definitive conclusions.
The results stated herein provide essential information not only for a better understanding of the plant-pathogen interaction but also for the application of genomic assisted breeding for anthracnose resistance improvement in common bean. This research has also shown the importance of the epistatic effects in genetic resistance to anthracnose, which has never been studied so far. Thereby, both main and epistatic interaction effects of genes or QTLs should be considered for a successful application of MAS, which provides an opportunity to use a pyramiding strategy for durable resistance. As well as providing useful tools for MAS of anthracnose resistance in common bean, this work also offers valuable clues for further study on cloning the candidate gene corresponding to the non-organ-specific QTL for resistance to race 1545 located on Chr5.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
This work was financially supported by the Ministerio de Economía y Competitividad (AGL2011-25562) and UE-FEDER Program. The authors would like to thank Junta de Andalucía (grant number P10-AGR-06931) and Campus de Excelencia Internacional Agroalimentario-CeiA3 for partially supporting this work financially. The authors gratefully acknowledge the technical support for this research from María Lores, Luis Godoy and José Cruz. The authors would also like to thank Dr James D. Kelly from Department of Plant, Soil and Microbial Sciences of the Michigan State University (USA) for providing strains of the fungus.
The Supplementary Material for this article can be found online at: http://www.frontiersin.org/journal/10.3389/fpls.2015.00141/abstract
Figure S1. Disease score distribution of the RIL population: (A) stem for race 23, (B) stem for race 1545, (C) leaf for race 23, (D) leaf for race 1545, and (E) petiole for race 1545. Red line curve: normal probability distribution; blue arrows indicate the mean value of both parental lines in each environment.
Abad, L., Wolters, P., Stucker, D., and Davis, P. (2006). “Advances in anthracnose stalk rot resistance,” in The Fifth National IPM Symposium “Delivering on a Promise.” Available online at: http://www.ipmcenters.org/ipmsymposiumv/posters/037.pdf
Altschul, S. F., Madden, T. L., Schaffer, A. A., Zhang, J., Zhang, Z., Miller, W., et al. (1997). Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 25, 3389–3402. doi: 10.1093/nar/25.17.3389
Alzate-Marín, A. L., Costa, M. R., Arruda, K. M., Barros, E. G., and Moreira, M. A. (2003). Characterization of the anthracnose resistance gene present in Ouro Negro (Honduras 35) common bean cultivar. Euphytica 133, 165–169. doi: 10.1023/A:1025551006444
Balmer, D., de Papajewski, D. V., Planchamp, C., Glauser, G., and Mauch-Mani, B. (2013). Induced resistance in maize is based on organ-specific defence responses. Plant J. 74, 213–225. doi: 10.1111/tpj.12114
Bhadauria, V., Banniza, S., Vandenberg, A., Selvaraj, G., and Wei, Y. (2011). EST analysis identifies proteins putatively secreted by the anthracnose pathogen Colletotrichum truncatum. BMC Genomics 12:e327. doi: 10.1186/1471-2164-12-327
Campa, A., Giraldez, R., and Ferreira, J. J. (2009). Genetic dissection of the resistance to nine anthracnose races in the common bean differential cultivars MDRK and TU. Theor. Appl. Genet. 119, 1–11. doi: 10.1007/s00122-009-1011-8
Campa, A., Rodríguez-Suárez, C., Giraldez, R., and Ferreira, J. J. (2014). Genetic analysis of the response to eleven Colletotrichum lindemuthianum races in a RIL population of common bean (Phaseolus vulgaris L.). BMC Plant Biol. 14:115. doi: 10.1186/1471-2229-14-115
Ferreira, J. J., Campa, A., and Kelly, J. D. (2013). “Organization of genes conferring resistance to anthracnose in common bean,” in Translational Genomic for Crop Breeding: Biotic Stress, eds R. K. Varshney and R. Tuberosa (Chichester: John Wiley & Sons), 151–182. doi: 10.1002/9781118728475.ch9
Ferrier-Cana, E., Geffroy, V., Macadre, C., Creusot, F., Imbert-Bollore, P., Sevicnac, M., et al. (2003). Characterization of expressed NBS-LRR resistance gene candidates from common bean. Theor. Appl. Genet. 106, 251–261. doi: 10.1007/s00122-002-1032-z
Gao, L., Tu, Z. J., Millet, B. P., and Bradeen, J. M. (2013). Insights into organ-specific pathogen defence responses in plants: RNA-seq analysis pf potato tuber-Phytophthora infestans interactions. BMC Genomics 14:340. doi: 10.1186/1471-2164-14-340
Geffroy, V., Macadre, C., David, P., Pedrosa-Harand, A., Sevignac, M., Dauga, C., et al. (2009). Molecular analysis of a large subtelomeric nucleotide-binding-site-leucine-rich-repeat family in two representative genotypes of the major gene pools of Phaseolus vulgaris. Genetics 181, 405–419. doi: 10.1534/genetics.108.093583
Geffroy, V., Sévignac, M., Billant, P., Dron, M., and Langin, T. (2008). Resistance to Colletotrichum lindemuthianum in Phaseolus vulgaris: a case study for mapping two independent genes. Theor. Appl. Genet. 116, 407–415. doi: 10.1007/s00122-007-0678-y
Geffroy, V., Sévignac, M., De Oliveira, J. C., Fouilloux, G., Skroch, P., Thoquet, P., et al. (2000). Inheritance of partial resistance against Colletotrichum lindemuthianum in Phaseolus vulgaris and co-localization of quantitative trait loci with genes involved in specific resistance. Mol. Plant Microbe. Interact. 13, 287–296. doi: 10.1094/MPMI.2000.13.3.287
Geffroy, V., Sicard, D., de Oliveira, J. C. F., Sévignac, M., Cohen, S., Gepts, P., et al. (1999). Identification of an ancestral resistance gene cluster involved in the coevolution process between Phaseolus vulgaris and is fungal pathogen Colletotrichum lindemuthianum. Mol. Plant Microbe. Interact. 12, 774–784. doi: 10.1094/MPMI.19188.8.131.524
Gonçalves-Vidigal, M. C., Cruz, A. S., Garcia, A., Kami, J., Vigidal Filho, P. S., and Sousa, L. L. (2011). Linkage mapping of the Phg-1 and Co-14 genes for resistance to angular leaf spot and anthracnose in the common bean cultivar AND 277. Theor. Appl. Genet. 122, 893–903. doi: 10.1007/s00122-010-1496-1
Gonçalves-Vidigal, M. C., Cruz, A. S., Lacanallo, G. F., Vidigal Filho, P. S., Sousa, L. L., Pacheco, C. M. N. A., et al. (2013). Co-segregation analysis and mapping of the anthracnose Co-10 and angular leaf spot Phg-ON disease-resistance genes in the common bean cultivar Ouro Negro. Theor. Appl. Genet. 126, 2245–2255. doi: 10.1007/s00122-013-2131-8
Holland, J. B. (2006). Estimating genotypic correlations and their standard errors using multivariate restricted maximum likelihood estimation with SAS Proc MIXED. Crop Sci. 46, 642–654. doi: 10.2135/cropsci2005.0191
Holland, J. B., Nyquist, W. E., and Cervantes-Martínez, C. T. (2003). “Estimating and interpreting heritability for plant breeding: an update,” in Plant Breeding Review, ed J. Janick (New Jersey: John Wiley & Sons), 9–112.
Kelly, J. D., Gepts, P., Miklas, P. N., and Coyne, D. P. (2003). Tagging and mapping of genes and QTL and molecular marker-assisted selection for traits of economic importance in bean and cowpea. Field Crops Res. 82, 135–154. doi: 10.1016/S0378-4290(03)00034-0
Lacanallo, G. F., Gonçalves-Vidigal, M. C., Vidigal Filho, P. S., Kami, J., and Gonela, A. (2010). Mapping of an Andean gene for resistance to anthracnose in the landrace Jalo Listras Pretas. Annu. Rep. Bean Improv. Coop. 53, 96–97.
Lefebvre, V., and Palloix, A. (1996). Both epistatic and additive effects of QTLs are involved in polygenic induced resistance to disease: a case study, the interaction pepper–Phytophthora capsici Leonian. Theor. Appl. Genet. 93, 503–511. doi: 10.1007/BF00417941
Liu, Y., Chen, L., Fu, D., Lou, Q., Mei, H., Xiong, L., et al. (2014). Dissection of additive, epistatic effect and QTL environment interaction of quantitative trait loci for sheath blight resistance in rice. Hereditas 151, 28–37. doi: 10.1111/hrd2.00026
Mahuku, G. S., Jara, C. E., Cajiao, C., and Beebe, S. (2002). Sources of resistance to Colletotrichum lindemuthianum in the secondary gene pool of Phaseolus vulgaris and in crosses of primary and secondary gene pools. Plant Dis. 86, 1383–1387. doi: 10.1094/PDIS.2002.86.12.1383
Martínez-Pacheco, M. M., Saucedo-Luna, J., Flores-García, A., Martínez-Muñoz, R. E., and Campos-García, J. (2009). Colletotrichum lindemuthianum (Sacc. & Magn.) Scrib. is a potential cellulases producer microorganism. Rev. Lat. Am. Microbiol. 51, 23–31.
Melotto, M., Balardin, R. S., and Kelly, J. D. (2000). “Host-pathogen interaction and variability of Colletotrichum lindemuthianum,” in Colletotrichum: Host Specific City, Pathology, and Host-Pathogen Interaction, eds D. Prusky, S. Freeman, and M. B. Dickman (St. Paul, MN: APS Press), 346–361.
Melotto, M., Coelho, M. F., Pedrosa-Harand, A., Kelly, J. D., and Camargo, L. E. A. (2004). The anthracnose resistance locus Co-4 of common bean is located on chromosome 3 and contains putative disease resistance-related genes. Theor. Appl. Genet. 109, 690–699. doi: 10.1007/s00122-004-1697-6
Méndez-Vigo, B., Rodríguez-Suárez, C., Pañeda, A., Ferreira, J. J., and Giraldez, R. (2005). Molecular markers and allelic relationships of anthracnose resistance gene cluster B4 in common bean. Euphytica 141, 237–245. doi: 10.1007/s10681-005-7075-x
Miklas, P. N., and Singh, S. P. (2007). “Common bean,” in Genome Mapping and Molecular Breeding in Plants, Sugar and Tuber Crops Pulses, Vol. 3, ed C. Kole (Berlin: Springer), 1–31. doi: 10.1007/978-3-540-34516-9_1
Oliver, R. P., and Ipcho, S. V. S. (2004). Arabidopsis pathology breathes new life into the necrotrophs-vs.-biotrophs classification of fungal pathogens. Mol. Plant Pathol. 5, 347–352. doi: 10.1111/j.1364-3703.2004.00228.x
Richard, M. M. S., Pflieger, S., Sévignac, M., Thareau, V., Blanchet, S., Li, Y., et al. (2014). Fine mapping of Co-x an athracnose resistance gene to a highly virulent strain of Colletotrichum lindemuthianum in common bean. Theor. Appl. Genet. 127, 1653–1166. doi: 10.1007/s00122-014-2328-5
Rodríguez-Guerra, R., Ramírez, M. T., Martínez de la Vega, O., and Simpson, J. (2003). Variation in genotype, pathotype and anastomosis groups of Colletotrichum lindemuthianum isolates from Mexico. Plant Pathol. 52, 228–235. doi: 10.1046/j.1365-3059.2003.00808.x
Rodríguez-Suárez, C., Ferreira, J. J., Campa, A., Pañeda, A., and Giraldez, R. (2008). Molecular mapping and intra-cluster recombination between anthracnose race-specific resistance genes in the common bean differential cultivars Mexico 222 and Widusa. Theor. Appl. Genet. 116, 807–814. doi: 10.1007/s00122-008-0714-6
Rodríguez-Suárez, C., Méndez-Vigo, B., Pañeda, A., Ferreira, J. J., and Giraldez, R. (2007). A genetic linkage map of Phaseolus vulgaris L. and localization of genes for specific resistance to six races of anthracnose (Colletotrichum lindemuthianum). Theor. Appl. Genet. 114, 713–722. doi: 10.1007/s00122-006-0471-3
Schmutz, J., McClean, P. E., Mamidi, S., Wu, G. A., Cannon, S. B., Grimwood, J., et al. (2014). A reference genome for common bean and genome-wide analysis of dual domestications. Nat. Genet. 46, 707–713. doi: 10.1038/ng.3008
Sicard, D., Michalakis, Y., Dron, M., and Neema, C. (1997). Genetic diversity and pathogenicity variation of Colletotrichum lindemuthianum in the three centers of diversity of its host, Phaseolus vulgaris. Phytopathology 87, 807–813. doi: 10.1094/PHYTO.19184.108.40.2067
Singh, A., Knox, R. E., DePauw, R. M., Singh, A. K., Cuthbert, R. D., Campbell, H. L., et al. (2014). Stripe rust and leaf rust resistance QTL mapping, epistatic interactions, and co-localization with stem rust resistance loci in spring wheat evaluated over three continents. Theor. Appl. Genet. 127, 2465–2477. doi: 10.1007/s00122-014-2390-z
Sousa, L. L., Gonçalves-Vidigal, M. C., Gonçalves, A. O., Vidigal Filho, P. S., Awale, H., and Kelly, J. D. (2013). Molecular mapping of the anthracnose resistance gene Co-15 in the common bean cultivars Corinthiano. Ann. Rep. Bean Improv. Coop. 56, 45–46.
Thabuis, A., Palloix, A., Pflieger, S., Daubeze, A. M., Caranta, C., and Lefebvre, V. (2003). Comparative mapping of Phytophthora resistance loci in pepper germplasm: evidence for conserved resistance loci across Solanaceae and for a large genetic diversity. Theor. Appl. Genet. 106, 1473–1485. doi: 10.1007/s00122-003-1206-3
Vigidal-Filho, P. S., Gonçalves-Vigida, M. C., Kelly, J. D., and Kirk, W. W. (2007). Sources of resistance to anthracnose in traditional common bean cultivars from Parana – Brazil. J. Phytopatholy. 155, 108–113. doi: 10.1111/j.1439-0434.2007.01203.x
Wang, C. S., Rutledge, J. J., and Gianola, D. (1994). Bayesian analysis of mixed linear models via Gibbs sampling with an application to litter size in Iberian pigs. Genet. Sel. Evol. 26, 91–115. doi: 10.1186/1297-9686-26-2-91
Yang, J., Hu, C. C., Hu, H., Yu, R. D., Xia, Z., Ye, X., et al. (2008). QTLNetwork: mapping and visualizing genetic architecture of complex traits in experimental populations. Bioinformatics 24, 721–723. doi: 10.1093/bioinformatics/btm494
Yang, S., Gao, M., Deshpande, S., Lin, S., Roe, B. A., and Zhu, H. (2007). Genetic and physical localization of an anthracnose resistance gene in Medicago truncatula Theor. Appl. Genet. 116, 45–52. doi: 10.1007/s00122-007-0645-7
Yuste-Lisbona, F. J., Santalla, M., Capel, C., García-Alcázar, M., De La Fuente, M., Capel, J., et al. (2012). Marker-based linkage map of Andean common bean (Phaseolus vulgaris L.) and mapping of QTLs underlying popping ability traits. BMC Plant Biol. 12:136. doi: 10.1186/1471-2229-12-136
Keywords: anthracnose, epistasis, Phaseolus vulgaris L., QTL, resistance gene clusters, NL genes
Citation: González AM, Yuste-Lisbona FJ, Rodiño AP, De Ron AM, Capel C, García-Alcázar M, Lozano R and Santalla M (2015) Uncovering the genetic architecture of Colletotrichum lindemuthianum resistance through QTL mapping and epistatic interaction analysis in common bean. Front. Plant Sci. 6:141. doi: 10.3389/fpls.2015.00141
Received: 22 December 2014; Paper pending published: 17 February 2015;
Accepted: 22 February 2015; Published: 17 March 2015.
Edited by:Paul Gepts, University of California, Davis, USA
Reviewed by:Eduard Akhunov, Kansas State University, USA
Swarup Kumar Parida, National Institute of Plant Genome Research, India
Copyright © 2015 González, Yuste-Lisbona, Rodiño, De Ron, Capel, García-Alcázar, Lozano and Santalla. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Marta Santalla, Grupo de Biología de Agrosistemas, Misión Biológica de Galicia-CSIC, PO Box 28, 36080 Pontevedra, Spain firstname.lastname@example.org
†These authors have contributed equally to this work.