ORIGINAL RESEARCH article
Sec. Plant Breeding
Volume 9 - 2018 | https://doi.org/10.3389/fpls.2018.00350
Identifying Quantitative Trait Loci (QTLs) and Developing Diagnostic Markers Linked to Orange Rust Resistance in Sugarcane (Saccharum spp.)
- 1Department of Agronomy, Genetics Institute, University of Florida, Gainesville, FL, United States
- 2Sugarcane Field Station, United States Department of Agriculture-Agricultural Research Service, Canal Point, FL, United States
- 3Key Laboratory of Genetics, Breeding and Multiple Utilization of Corps, Center for Genomics and Biotechnology, Ministry of Education and Fujian Provincial Key Laboratory of Haixia Applied Plant Systems Biology, Fujian Agriculture and Forestry University, Fujian, China
Sugarcane (Saccharum spp.) is an important economic crop, contributing up to 80% of table sugar used in the world and has become a promising feedstock for biofuel production. Sugarcane production has been threatened by many diseases, and fungicide applications for disease control have been opted out for sustainable agriculture. Orange rust is one of the major diseases impacting sugarcane production worldwide. Identifying quantitative trait loci (QTLs) and developing diagnostic markers are valuable for breeding programs to expedite release of superior sugarcane cultivars for disease control. In this study, an F1 segregating population derived from a cross between two hybrid sugarcane clones, CP95-1039 and CP88-1762, was evaluated for orange rust resistance in replicated trails. Three QTLs controlling orange rust resistance in sugarcane (qORR109, qORR4 and qORR102) were identified for the first time ever, which can explain 58, 12 and 8% of the phenotypic variation, separately. We also characterized 1,574 sugarcane putative resistance (R) genes. These sugarcane putative R genes and simple sequence repeats in the QTL intervals were further used to develop diagnostic markers for marker-assisted selection of orange rust resistance. A PCR-based Resistance gene-derived maker, G1 was developed, which showed significant association with orange rust resistance. The putative QTLs and marker developed in this study can be effectively utilized in sugarcane breeding programs to facilitate the selection process, thus contributing to the sustainable agriculture for orange rust disease control.
Sugarcane (Saccharum spp.) is one of the most important economic crops, cultivated on ∼27.1 million hectares in over 100 countries with a worldwide harvest of 1.9 billion metric tons and a gross production value of $81.5 billion (FAO, 2013, 2014). As the most important sugar resource, sugarcane contributes up to 80% of sugar production in the world. Additionally, sugarcane accounts for ∼60% of global bio-ethanol production, with an energy output-to-input ratio five times higher than that of maize (Waclawovsky et al., 2010; Dahlquist, 2013). Nevertheless, like many other crops, sugarcane production has been vulnerable by many diseases.
Orange rust, caused by Puccinia kuehnii, was first reported in Australia more than 100 years ago; however, this disease only began to draw attention due to its devastating epidemic in 2000 on a major Australian sugarcane cultivar Q124 (Braithwaite et al., 2009). Orange rust was first identified in the Western Hemisphere, specifically in Florida, United States, in 2007 (Comstock et al., 2008) and shortly after in Guatemala (Ovalle et al., 2008), Mexico, El Salvador, Panama (Flores et al., 2009), Costa Rica, Nicaragua (Chavarria et al., 2009), Brazil (Barbasso et al., 2010), Colombia (Cadavid et al., 2012), and recently Ecuador (Garcés et al., 2014). The disease immediately became an emerging threat to sugarcane production and breeding programs since its discovery (Comstock et al., 2008, 2010), which can cause up to 50% cane yield loss in its susceptible hosts (Magarey et al., 2011).
In the United States, Florida is the top sugarcane producer for sugar production with ∼28.0 million tons of cane and a gross value of $909.7 million annually (FAO, 2014). The majority of commercial cultivars grown in Florida and the parental germplasm used in the sugarcane breeding program are susceptible to the orange rust (Comstock et al., 2010), marking this disease a major concern in sugarcane industry. Application of fungicides has proven to increase sugar production from 7.9 to 26% on CP80-1743 in commercial fields (Comstock et al., 2010). However, fungicide treatments can lead to input cost increasing as well as potential environmental problems, which is against the principles advocated in sustainable agriculture (Lichtfouse et al., 2009). Therefore, developing and utilizing resistant cultivars is undoubtedly the favorable approach for sugarcane breeders and farmers for the disease control.
Screening for disease resistant sugarcane breeding materials can be performed by evaluating plant reactions both in the field and/or in greenhouse after rust spore inoculation artificially. The phenotype-based selection is direct and efficient when plants are inoculated with viable inoculum at the right growth stages in an environment favoring disease development. Unfortunately, selection based on phenotyping is time-consuming, labor-intensive, and environment-dependent due to requirement of large space, inoculation process, and possible escapes if favorable conditions or plant developmental stages are not achieved during inoculation. Marker-assisted selection (MAS) of resistant materials is a desirable alternative method, in which resistant individuals can be selected in the laboratory in less than 24 h by testing the presence of molecular markers linked to disease resistance. The MAS depends on reliable markers that are tightly linked to genes or genomic regions controlling disease resistance through quantitative trait loci (QTL) mapping and association study. Though sugarcane is a complex auto-polyploid species with large genome size, several disease resistance loci have been identified, such as Bru1 (Daugrois et al., 1996) and Bru2 (Raboin et al., 2006) for brown rust resistance, four DNA markers for pachymetra root rot and leaf scald resistance, five for Fiji leaf gall resistance, and 11 for smut resistance (Wei et al., 2006), a major QTL for yellow spot resistance (Aljanabi et al., 2007), a major quantitative trait allele for the Sugarcane yellow leaf resistance (Costet et al., 2012b) were identified in sugarcane. Once identified and validated, these molecular markers could be used in MAS to quickly identify and select desirable resistant resources or progeny, though no disease resistance gene has been cloned in sugarcane yet. For example, since the discovery of Bru1 and the identification of its co-segregating markers (R12H16 and 9O20-F4) (Costet et al., 2012a), the diagnostic markers have widely been used in sugarcane breeding programs. In Florida sugarcane breeding program, the markers for Bru1 are being utilized to evaluate brown rust resistance in sugarcane germplasm and hybrid clones (Glynn et al., 2013). It was estimated that 27% of the clones used for crossing (a total of 1027) contained Bru1, and the frequency of Bru1 in sugarcane clones increased from 15% (1975–1985) to 47% (2002–2012) after brown rust was introduced to Florida.
Although reliable molecular markers are of great value for controlling disease effectively in a way fitting the goals of sustainable agriculture, so far none has been reported to be linked to orange rust resistance in sugarcane. The objectives of this study were (1) to phenotypically evaluate the orange rust resistance reaction in a mapping population; (2) to map QTLs controlling orange rust resistance in sugarcane; and (3) to develop diagnostic markers for MAS of orange rust resistance in sugarcane. The putative QTLs and markers developed in this study can be effectively utilized in sugarcane breeding programs to expedite release of resistant cultivars, thus contributing to sustainable agriculture for orange rust disease control.
Materials and Methods
The sugarcane mapping population comprised of 173 F1 progeny that were derived from a cross between sugarcane clones CP95-1039 and CP88-1762, which were developed by the United States Department of Agriculture (USDA), Agricultural Research Service (ARS) Sugarcane Field Station at Canal Point, FL, United States. Individuals of the progeny were clonally propagated for leaf sampling and phenotype evaluation. The whole population was genotyped using genotyping by sequencing (GBS), and sequence variations were called according to Yang et al. (2017a). The GBS libraries and sequencing were performed at the Institute of Genomic Diversity, Cornell University following the optimized protocol (Elshire et al., 2011) with a rare cutting restriction enzyme PstI (CTGCAG) for library construction and a 96-plex (95 DNA samples and one negative control) for sequencing on the Illumina HiSeq 2000 platform (Yang et al., 2017b). The cleaned reads were aligned to sorghum genome v3.0 (Paterson et al., 2009) for calling single nucleotide polymorphisms (SNPs) using seven different SNP callers described by Yang et al. (2017b). Only the single dose markers were used for genetic map construction using Joinmap 4.0 (Van Ooijen, 2006). Two high density genetic maps were constructed with one for CP95-1039 including 2,453 markers and a total length of 4224.4 cM and the other for CP88-1762 including 2,154 markers and a total length of 4373.2 cM (Yang et al., 2017b). For validating the diagnostic marker developed in this study, a diversity panel with 165 sugarcane clones derived from multiple crosses in multiple years (Supplementary Table S1) was used.
Evaluation of Orange Rust Resistance
The whole F1 mapping population along with the two parental clones, CP95-1039 and CP88-1762, were used in orange rust resistance evaluation. CP95-1039 is resistant to orange rust, whereas CP88-1762 is susceptible to the disease (Supplementary Figure S1). Rust resistance evaluation of this population was conducted using artificial inoculation methods under both field and greenhouse conditions.
The inoculum preparation and inoculation were performed following the protocol described by Sood et al. (2009). Briefly, sugarcane orange rust urediniospores were collected by vacuuming the abaxial side of symptomatic leaves since P. kuehnii is biotrophic and difficult to culture in medium. A large number of leaves diagnosed with orange rust disease symptoms from multiple highly susceptible cultivars were used to collect enough inoculum for this experiment. Freshly collected urediniospores were used as inoculum or stored at -70°C. Rust spores were suspended in sterile distilled water containing 0.1% (V/V) Tween-20 and 0.002% 1-nonanol with an adjusted concentration of 104 urediniospores/ml.
For the field experiment, seedlings from the mapping population grown in the greenhouse were transplanted in the field on August 2010 for population establishment at the USDA-ARS Sugarcane Field Station, Canal Point, FL, United States. The planting canes were trimmed above the soil in February 2011 for first ratoon regrowth. The inoculum was applied to three stalks of each genotype when the sugarcane seedlings were 2 months old after transplanting in October 2010, and when the regrown ratoon plants were 2 and 3 months old after cutting in April 2011 and May 2011, respectively. To infect plants, one-third of the tips of the uppermost leaves were cut off as a marker, and a 0.5-ml aliquot of the spore suspension was deposited inside the leaf whorl of the marked stalks. Three weeks post-inoculation, any rust signs or symptoms were recorded from the newly emerged leaf following a 0–4 scale (Sood et al., 2009). A score of 0 indicated that the plants were asymptomatic. A score of 1 corresponded to the presence of some chlorotic flecks without any pathogen spores on the leaves. A score of 2 indicated that leaves had some brown discoloration, but no pustules or sporulation. A score of 3 corresponded to the presence of ≤10 pustules and sporulation on the leaves. A score of 4 corresponded to >10 pustules with massive sporulation. The inoculation on plant cane was scored one time, while each inoculation on ratoons was scored by two researchers independently. Five observations were conducted after artificial inoculation in the field; one was on the plant cane and the other four were on ratoons.
For the greenhouse experiment, individuals of the mapping population were clonally propagated with four plants per sugarcane clone in 46 cm × 38 cm × 10 cm flats filled with Miracle-Gro® Potting Mix (Scotts Miracle-Gro® Company, Marysville, OH, United States). Scott’s Osmocote Classic 14-14-14 fertilizer (Scotts-Sierra, Marysville, OH, United States) and daily irrigation were applied to maintain healthy plants. Each flat contained four genotypes, with four stalks for each genotype. The inoculum (prepared in the same way as described above) was applied to three healthy stalks (4 months old) with a complete leaf whorl developed. Three weeks after inoculation, any rust signs or symptoms were scored from the newly emerged leaf following the 0 (no disease present) to 4 (highly susceptible) scale (Sood et al., 2009). Artificial inoculation experiments were conducted under greenhouse conditions in March 2014, November 2014 and May 2015, respectively. The first greenhouse inoculation experiment was performed on plant cane, whereas the other two inoculation experiments were on ratoons. The inoculation and scoring of rust symptoms in the greenhouse were the same as described above. Therefore, three observations were conducted after artificial inoculation in the greenhouse.
Data Analysis and QTL Mapping
A pairwise Pearson correlation among scores of eight observations of orange rust resistance was calculated using the ‘Hmisc’ package in R3.0.2 (R Development Core Team, 2013; Harrell, 2018). Analysis of variance was performed using the ‘nlme’ package in R (Pinheiro et al., 2009). Broad sense heritability (H2) was estimated using the formula below:
where and are the genetic variance of progeny lines, and the error variance, respectively.
The eight observations of orange rust resistance were manually checked for ‘escape’ inoculation (progeny with a score of two or more after inoculation were considered as ‘susceptible,’ and then only observations with scores higher than two were used to calculate average scores), and average scores were used for QTL analysis. A pseudo-testcross strategy was employed for QTL analysis using WinQTLCart 2.5 (Wang et al., 2012). Composite interval mapping (CIM) was performed using forward and backward stepwise regressions to select markers as cofactors with 10 cM window size and 1 cM walking speed in WinQTLCart 2.5 (Wang et al., 2012). The QTL nomenclature was used according to McCouch (2008).
Identification of Sugarcane Putative R Genes
Five sugarcane transcript sequence databases (Table 1), including a total of two million sugarcane transcripts, were used to identify sugarcane transcripts having homologous genes in the sorghum genome. Redundant transcripts were removed by using cdhit/4.6 (sequence identity threshold 0.95) to reduce subsequently computational and manual efforts (Li and Godzik, 2006; Fu et al., 2012). The non-redundant sugarcane transcripts and their corresponding sorghum orthologs (Paterson et al., 2009) were identified and verified by using reciprocal BLASTN. The sugarcane transcripts with corresponding sorghum orthologs were used for subsequent analysis.
TABLE 1. Sugarcane transcript database sources used for identification of sugarcane putative resistance genes.
Sugarcane putative R genes were identified by two complementary methods. First, BLASTN of sugarcane transcripts with corresponding sorghum orthologs was performed against known plant R genes, which included 112 reference resistance genes with experimental support from the Plant Resistance Genes database (PRGdb)1 (Sanseverino et al., 2013). Second, InterProScan 5.0 was used to characterize sugarcane nucleotide-binding site (NBS)–leucine-rich repeat (LRR) genes from the sugarcane transcripts with corresponding sorghum orthologs following the procedures described previously (Jones et al., 2014; Yang and Wang, 2016). All bioinformatics analyses were performed locally at the University of Florida High Performance Computing Center using command lines.
Development of PCR-Based Markers for Orange Rust Resistance Selection
Sugarcane simple sequence repeats (SSRs) sequences were collected from published literatures and aligned to sorghum genome 3.0 (Liu et al., unpublished data). Primers for putative R genes were designed using Primer3.2 Primers across introns based on the sorghum gene models were selected. Only primers uniquely aligned to the sorghum genome and aligned to sugarcane transcripts with no more than two mismatches were selected for synthesis (Invitrogen, Carlsbad, CA, United States). All the SSRs and gene markers were first screened using two parental clones, CP95-1039 (resistant to orange rust) and CP88-1762 (susceptible to orange rust). A touchdown PCR program was used as described by Li et al. (2011) for amplification. In brief, the reaction mixture was incubated at 95°C for 5 min, then five cycles of 60 s of denaturing at 96°C, 5 min of annealing at 68°C with a decrease of 2°C in each subsequent cycle, and 1 min of extension at 72°C; For another five cycles, the annealing temperature started at 58°C for 2 min with a decrease of 2°C for each subsequent cycle; PCR continued through an additional 25 cycles of 60 s at 96°C, 1 min at 50°C, and 1 min at 72°C with a final extension at 72°C for 5 min. PCR products were run on a 1% agarose gel in a horizontal electrophoresis apparatus (Bio-Rad Laboratories, Hercules, CA, United States). Polymorphic markers were used to genotype the mapping population of 173 F1 progeny derived from the cross between CP95-1039 and CP88-1762. Candidate diagnostic markers (PCR-based) were further genotyped in a diversity panel with 165 sugarcane clones derived from multiple crosses in multiple years (Supplementary Table S1). These 165 sugarcane clones were also screened phenotypically for orange rust resistance in the field following the protocol described above in the evaluation of orange rust resistance section. The Student t-tests of each marker were conducted for mean orange rust scores by using functions in Excel 2010 (Microsoft Corp., Redmond, WA, United States).
Based on artificial inoculation data, progeny with scores of 0, 1, and 2, that were not sporulating, were categorized as phenotypic ‘resistant,’ and progeny with scores of 3 and 4, that were sporulating, were categorized as phenotypic ‘susceptible.’ The 173 progeny were considered as ‘resistant’ or ‘susceptible’ based on sporulation for marker evaluation. Based on the closest flanking SNP markers to the significant QTLs, each individual in the whole population was characterized as genotypic ‘resistant’ or genotypic ‘susceptible.’ Selection accuracy was defined as the percentage of individuals which were both phenotypic and genotypic resistant among the total number of genotypic resistant individuals. Selection efficiency was defined as the percentage of individuals which were both phenotypic and genotypic resistant among the total number of phenotypic resistant individuals in the mapping population. Selection accuracy and selection efficiency were calculated following formula:
Selection accuracy = (the number of both phenotypic and genotypic resistant individuals/the number of genotypic resistant individuals) ∗ 100.
Selection efficiency = (the number of both phenotypic and genotypic resistant individuals/the number of phenotypic resistant individuals in the mapping population) ∗ 100.
Phenotypic Data Analysis
Orange rust resistance reactions of the mapping population showed a nearly normal distribution (Figure 1), supporting orange rust resistance segregating in this bi-parental mapping population is controlled by multiple genes/alleles. Orange rust resistance reactions of artificial inoculations were highly correlated (P < 0.0001) with an average coefficient of 0.64 (Table 2). Heritability of orange rust resistance was 0.67, indicating that the disease resistance in the mapping population was largely controlled by genetic factors. In this study, the mean scores of the eight rust resistance observations were used for subsequent QTL analysis, which can reduce discrepancy caused by experimental errors in single artificial inoculations.
FIGURE 1. Distribution of orange rust (OR) resistance reactions in the mapping population across multiple trials observations. F, Field; G, Greenhouse; AVR, average.
TABLE 2. Pearson’s correlation coefficients among rust resistance reactions on the F1 population derived from the cross CP95-1039 × CP88-1762.
High density genetic maps of both parental clones (Yang et al., 2017b) were used for detecting QTLs employing WinQTLCart 2.5 (Wang et al., 2012). Three QTLs controlling orange rust resistance in sugarcane, qORR109, qORR4 and qORR102, were detected in the mapping population (Figure 2 and Table 3). Of the three QTLs, one was located on linkage group (LG) 109 of the resistant parent CP95-1039 map, explaining 58% of the phenotypic variance with a logarithm of odds (LOD) of 19.3. The other two QTLs were located on LG 4 and LG 102 of the CP88-1762 map, explaining 12 and 8% of the phenotypic variance with a LOD of 5.4 and 3.4, respectively. Of the three QTLs, heterozygous genotypes of qORR109 and qORR102 positively contributed to orange rust resistance, whereas that of QTL qORR4 showed a negative contribution (Table 3). The three nearest markers, 2SNP350, 5SNPUN3354 and 3SNP3092, were 2 cM, 0.2 cM and 2.2 cM away from their corresponding QTL peak, qORR109, qORR4 and qORR102, respectively.
FIGURE 2. Logarithm of odds (LOD) score plots of the quantitative trait loci (QTLs) for sugarcane orange rust resistance in the regions harboring the three identified QTL clusters.
TABLE 3. Quantitative trait loci (QTL) associated with orange rust resistance of the F1 population derived from the cross CP95-1039 × CP88-1762.
To further study the genetic effect of the three QTLs, the phenotypic variation was dissected for each QTL and their combination (Table 4). Progeny containing the three susceptible QTL intervals based on the nearest marker alleles (Homozygote for QTL qORR109 and qORR102, and heterozygote for QTL qORR4) had a mean disease resistance score of 3.62, which served as the baseline for calculating orange rust disease reduction in the mapping population. With a single QTL, the disease severity can be reduced by 37.8 and 46.8% for QTLs qORR4 and qORR102, respectively. Since only one progeny existed in this category of single QTL qORR109 with heterozygote, the calculation of its contribution to rust resistance was ignored. With a combination of two QTLs, the disease severity could be reduced by 45.9, 34.3, and 19.6%, separately (Table 4). The combination of the three positive QTLs (Heterozygote for QTL qORR109 and qORR102, and homozygote for QTL qORR4) could reduce the disease severity by 56.4%.
TABLE 4. Phenotypic effects of quantitative trait loci (QTL) allele combinations on orange rust resistance in the F1 population from the cross CP95-1039 × CP88-1762.
Sugarcane SSRs and Putative R Genes
Out of 6,149 available SSR primers recruited from a literature search, a total of 1,095 sugarcane SSR primer pairs were uniquely aligned to the sorghum reference genome (Liu et al., unpublished data). Of all the mapped SSRs, 96 were found in the vicinity of three orange rust resistance QTL intervals and were close to the nearest markers of the QTLs according to the sorghum genome (Supplementary Table S2).
A custom-made sugarcane transcript database was formed including 700,194 sequences after removing redundant sequences. In total, 24,665 sugarcane transcripts were identified with orthologs in the sorghum genome. To identify sugarcane putative resistance genes, only the transcripts with orthologs in sorghum were utilized for subsequent analysis for a collinear comparison and application of the QTL intervals based on the sorghum genome. Out of the 24,665 sugarcane transcripts, 1,574 were characterized as putative R genes (Supplementary Table S3). Transcript sequences, corresponding sorghum orthologs, and reference R genes for sugarcane putative R genes were deposited in Supplementary Table S3.
Developing Diagnostic Markers for MAS of Orange Rust Resistance in Sugarcane
Sugarcane SSRs in vicinity of the three orange rust resistance QTL intervals were used for diagnostic marker development. Of the 96 SSRs, only two showed robust polymorphisms between the two parental clones, CP95-1039 and CP88-1762. In addition, 61 gene-derived primer pairs crossing introns of 18 sugarcane putative R gene sequences were designed (Supplementary Table S2). Of the 61 pairs of R gene-derived primers, three showed polymorphisms between the two parental clones. The five polymorphic markers (two SSRs and three gene-derived) were first genotyped in the F1 mapping population. Interestingly, two gene-derived markers, G1 (P < 0.001) and M16 (P < 0.05), were significantly associated with orange rust resistance based on single marker analysis (Table 5 and Supplementary Figure S2). The candidate diagnostic markers were further validated by genotyping 165 diverse sugarcane clones (Supplementary Table S1). The result further confirmed that G1 maker was associated with orange rust resistance (P < 0.05), while M16 was weakly associated. Based on the F1 mapping population, selection accuracy of G1 marker for orange rust resistance was 65.8% and the selection efficiency was 85.6%.
Previously restricted to Asia, and the Pacific regions, sugarcane orange rust has also been detected in multiple countries in the western hemisphere (Comstock et al., 2008; Ovalle et al., 2008; Chavarria et al., 2009; Flores et al., 2009; Barbasso et al., 2010; Cadavid et al., 2012; Garcés et al., 2014). The outbreak of this disease could result in devastating economic loss, with an estimation of Aus$150–210 million in 2000 in Australia (Braithwaite et al., 2009). In Florida, orange rust has gained full attention due to its sudden emergence and negative impact on sugarcane production. In one crop season, economic losses associated with orange rust were estimated at $40 million (Dixon et al., 2010). It continues to be a serious disease for growers for several reasons. From a chemical perspective, optimal fungicide selection, and the rate, timing, and frequency of applications to prevent or control this newly emerged disease are unknown. Though several fungicides (e.g., strobilurin class fungicides and triazole class fungicides) have recently been registered and used to manage orange rust, more research is needed to fine-tune control recommendations that are economically reasonable (Rott et al., 2014). Biologically, orange rust pathogen releases a high density of spores into the air, and tolerates high temperatures, which allows the disease to develop throughout most of the sugarcane growing season and production areas, assuming infection conditions are suitable (Rott et al., 2014). Finally, the use of resistant host plants is currently limited because 77.5% of commercial cultivars grown in Florida and 41% of parental germplasm are susceptible to orange rust (Comstock et al., 2010). Unfortunately, it takes 12–14 years to release a new sugarcane cultivar with disease resistance, and even longer to propagate the cultivar for commercial production.
Compared to fungicide treatment, developing and using resistant cultivars is undoubtedly the thumbs-up approach for disease control because of a low input cost and healthy environment for sustainable sugarcane production. With reliable diagnostic markers developed, MAS would make this approach achievable for sugarcane breeding programs by quickly identifying resistant germplasm for crossing and selecting resistance progeny. In this study, we evaluated an F1 segregating population following a pseudo-testcross strategy for orange rust resistance QTL identification. Since the sugarcane is an outcrossing species and the parental clones for the F1 population were heterozygous, the F1 population was segregating genetically in the way similar as the backcross population of self-pollinating species if only the loci with homozygous genotype in one parent and heterozygous genotype in the other parent were considered. With multiple evaluation of the phenotype and high stringent genotypes of this segregating F1 population, we identified three QTLs controlling orange rust resistance for the first time, which can explain 58, 12, and 8% of the phenotypic variation, separately. Moreover, a PCR-based R gene-derived maker, G1 was developed, which showed significant association with orange rust resistance. The results of this research would be valuable for expediting release of resistant cultivars, thus contributing to the sustainable agriculture for orange rust disease control.
QTL Analysis for Orange Rust Resistance
By performing QTL analysis using high density genetic maps and the mean scores of the eight rust resistance evaluations, for the first time we identified three QTLs controlling orange rust resistance. It’s noteworthy that GBS empowered QTL identification in sugarcane in this study. We used multiple SNP callers to maximize the detection of single dose markers, applied a robust filtering process on raw SNPs for each individual sample to improve SNP calling, and exploited a high LOD (≥10) in genetic map construction to avoid confusing linkage from different homoeologs (Yang et al., 2017b), which provide high density and reliable genetic maps for QTL analysis in this study. The three QTLs combined together could reduce 56.4% of disease severity compared with progeny containing the three susceptible QTL alleles (Table 3), indicating pyramiding these positive QTL alleles into sugarcane cultivars through MAS could effectively reduce the impact of the disease. So far, we have no evidence on whether these QTLs were linked to other disease resistance or agronomically important traits yet. The strategy for pyramiding multiple genes to improve disease resistance has been discussed in bean (Kelly et al., 1995), wheat (Liu et al., 2000), and rice (Singh et al., 2001; Jiang et al., 2004). Interestingly, two QTLs, qORR4 and qORR102, were identified from the same susceptible parental clone, CP88-1762, suggesting breeding materials that have robust agronomic traits but are susceptible to the disease should not directly be eliminated from sugarcane breeding programs as they may contain valuable minor resistance QTL alleles.
Although the molecular mechanisms of pathogen recognition and resistance to orange rust disease remain unclear, the resistance characterized in the current mapping population is less likely to be a single pathogen-race specific resistance since a mixture of orange rust spores collected from multiple highly susceptible cultivars from multiple fields were used to screen for resistance. Therefore, most likely a broad-based and durable resistance was identified within the segregating population. Although no genetic diversity has been reported in P. kuehnii so far, however, a change of orange rust pathogenicity has most likely been occurred in Florida. Earlier, major sugarcane cultivars CP88-1762 and CP89-2143 had been symptomless or resistant before 2010–2011, however, recently exhibited severe disease symptoms, suggesting a change in P. kuehnii pathogenicity (Comstock et al., 2010; Philippe Rott, Personal Communication). The resistance genes identified in the mapping population of this study may recognize multiple elicitors from orange rust pathogens, including pathogen-associated molecular patterns and/or weak effectors. Subsequent defense responses in sugarcane can be triggered to reduce the damage with contributing effects from each resistance gene. Since spores were artificially injected in the leaf, so lesions were also observed on leaves in the resistant type, but no sporulation was observed even after 4 weeks of infection (Supplementary Figure S1). The resistance in this mapping population may belong to a general non-specific resistance, in which sporulation of the fungus is inhibited to a low level to reduce the disease spread and severity (Kiraly et al., 2007). Therefore, this resistance should be fully exploited for breeding resistant sugarcane cultivars.
Development of Diagnostic Markers for MAS of Orange Rust Resistance in Sugarcane
Release of orange rust resistant cultivars could be expedited by quickly identifying parental resistant germplasm and selecting resistant progeny through MAS using reliable molecular markers, especially PCR-based markers, which are easy to run for breeders without relying on fancy equipment and complicated procedure. In this study, we constructed a sugarcane putative R gene database (Supplementary Table S3) by bioinformatics analysis of several transcript databases. A total of 1,574 sugarcane putative R genes were annotated as NBS-LRR genes. Information on these putative R genes, with corresponding sorghum gene models, corresponding reference R genes, annotation, sugarcane transcript sequences (Supplementary Table S3) will be an important source to sugarcane research community for sugarcane disease resistance genes mining, genetic map construction, and PCR-based marker development.
In this study, we designed 61 R gene-derived primers according to the orange rust QTL intervals and the sugarcane putative R database. Out of the 61 primers, three (4.9%) showed polymorphisms with different fragment sizes between the two parental clones, indicating this strategy of crossing introns could functionally generate markers for subsequent analysis. Even more interestingly, after single marker analysis in the bi-parental mapping population and the diversity sugarcane clones, G1 marker was significantly linked to orange rust resistance, while M16 showed weak association (Table 5). We also screened 96 SSR markers; however, none of them was linked to orange rust resistance. Lack of linkage between the two polymorphic SSRs to orange rust resistance could be due to that the marker were far away from the QTL; or most likely that the marker linked with QTLs in a repulsion phase. The sugarcane genome is complex, with up to 12 sets of homo(eo)logs. A single locus may have multiple alleles and/or alleles with different dosages. Theoretically, only markers in coupling linkage with target QTLs contributing alleles can be used for MAS. Therefore, a successful rate of identifying PCR-based markers in the expected linkage is relatively low using this strategy.
Furthermore, we found that the target gene of G1 marker, Sobic.002G166150, is homologous to RFO1, a dominant Arabidopsis disease resistance gene, which encoded a wall-associated receptor-like kinase and conferred resistance to a broad spectrum of Fusarium races (Diener and Ausubel, 2005), further indicating that the resistance QTLs identified in this study carried durable instead of race-specific resistant genes. In this study, we used a mixture of orange rust spores collected from infected leaf of multiple sugarcane cultivars as inoculum (hard to culture this pathogen in medium). Thus the inoculum was most likely a mixture of multiple races, though the Puccinia race variation was not reported in North America yet. If the inoculum contained multiple races (we speculated that multiple races of Puccinia existed due to observation of the changing of host resistance over years; Comstock et al., 2010; Philippe Rott, Personal Communication), then resistance to the most popular orange rust pathogen races would be detected in our experiments. In another word, the QTLs we identified should harbor durable orange rust resistance gene(s) instead of a single race-specific resistance gene. The resistance provided by RFO1 is quantitative and stronger resistance achieved if combined with other loci (additive effect). We also observed the same patterns in our results (Table 4). However, further experiments are needed to confirm whether sugarcane ortholog of Sobic.002G166150 is the gene contributing to the orange rust resistance in sugarcane, and to elucidate the interactions among the three QTLs in sugarcane.
QTL analysis revealed three QTLs controlling orange rust resistance, which can explain 58, 12, and 8% of the phenotypic variation, respectively. The three QTLs together can reduce 56.4% of disease severity compared with progeny containing the three susceptible QTL alleles. This is the first time that QTLs controlling orange rust resistance in sugarcane have been reported. To develop diagnostic markers for orange rust resistance, we identified 1,574 sugarcane putative R genes and aligned publicly available sugarcane SSRs to the sorghum genome. Out of 61 gene-derived markers and 96 SSRs screened, five were polymorphic with different fragment sizes between the two parental clones. The G1 marker was significantly linked to orange rust resistance based on single marker analysis. The putative QTLs and marker developed in this study can be utilized in sugarcane breeding programs to utilize MAS of orange rust resistance in sugarcane.
JW conceived and designed the experiments. XY, SS, MI, and SM performed the experiments. XY performed the data analysis. XY and JW drafted the manuscript. XY, JW, SS, SM, EH, MI, and JC revised the manuscript. All authors read and approved the final manuscript.
This research is financially supported by United States Department of Agriculture – Agricultural Research Service CRIS projects 6030-22000-010-00D, USDA National Institute of Food and Agriculture, Hatch Project 1011664, and the Florida Sugar Cane League. Publication of this article was funded in part by the University of Florida Open Access Publishing Fund.
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.
We gratefully thank Dr. Neil Glynn for developing and providing the population. We also express their gratitude to Kay McCorkle, Ken Peterkin and Brittany Read at Sugarcane Field Station, USDA, ARS for their help during disease screening and field trial.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2018.00350/full#supplementary-material
Aljanabi, S. M., Parmessur, Y., Kross, H., Dhayan, S., Saumtally, S., Ramdoyal, K., et al. (2007). Identification of a major quantitative trait locus (QTL) for yellow spot (Mycovellosiella koepkei) disease resistance in sugarcane. Mol. Breed. 19, 1–14. doi: 10.1007/s11032-006-9008-3
Barbasso, D., Jordo, H., Maccheroni, W., Boldini, J., Bressiani, J., and Sanguino, A. (2010). First report of Puccinia kuehnii, causal agent of orange rust of sugarcane, in Brazil. Plant Dis. 94:1170. doi: 10.1094/PDIS-94-9-1170C
Braithwaite, K. S., Croft, B. J., Magarey, R. C., and Scharaschkin, T. (2009). Phylogenetic placement of the sugarcane orange rust pathogen Puccinia kuehnii in a historical and regional context. Aust. Plant Pathol. 38, 380–388. doi: 10.1071/AP09012
Cadavid, M., Ángel, J. C., and Victoria, J. I. (2012). First report of orange rust of sugarcane caused by Puccinia kuehnii in Colombia. Plant Dis. 96, 143. doi: 10.1094/PDIS-05-11-0406
Chavarria, E., Subiros, F., Vega, J., Ralda, G., Glynn, N. C., Comstock, J. C., et al. (2009). First report of orange rust of sugarcane caused by Puccinia kuehnii in Costa Rica and Nicaragua. Plant Dis. 93, 425. doi: 10.1094/PDIS-93-4-0425C
Comstock, J., Glynn, N., and Davdison, R. (2010). Sugarcane rusts in Florida. Proc. Int. Soc. Sugar Cane Technol. 27, 1–9.
Comstock, J. C., Sood, S. G., Glynn, N. C., Shine, J. M. Jr., McKemy, J. M., and Castlebury, L. A. (2008). First report of Puccinia kuehnii, causal agent of orange rust of sugarcane, in the United States and Western Hemisphere. Plant Dis. 92, 175. doi: 10.1094/PDIS-92-1-0175A
Costet, L., Le Cunff, L., Royaert, S., Raboin, L., Hervouet, C., Toubi, L., et al. (2012a). Haplotype structure around Bru1 reveals a narrow genetic basis for brown rust resistance in modern sugarcane cultivars. Theor. Appl. Genet. 125, 825–836. doi: 10.1007/s00122-012-1875-x
Costet, L., Raboin, L., Payet, M., D’Hont, A., and Nibouche, S. (2012b). A major quantitative trait allele for resistance to the Sugarcane yellow leaf virus (Luteoviridae). Plant Breed. 131, 637–640. doi: 10.1111/j.1439-0523.2012.02003.x
Dahlquist, E. (2013). Biomass as Energy Source: Resources, Systems and Applications. London: CRC Press.
Daugrois, J. H., Grivet, L., Roques, D., Hoarau, J. Y., Lombard, H., Glaszmann, J. C., et al. (1996). A putative major gene for rust resistance linked with a RFLP marker in sugarcane cultivar ‘R570’. Theor. Appl. Genet. 92, 1059–1064. doi: 10.1007/BF00224049
Diener, A. C., and Ausubel, F. M. (2005). RESISTANCE TO FUSARIUM OXYSPORUM 1, a dominant Arabidopsis disease-resistance gene, is not race specific. Genetics 171, 305–321. doi: 10.1534/genetics.105.042218
Dixon, L. J., Castlebury, L. A., Aime, M. C., Glynn, N. C., and Comstock, J. C. (2010). Phylogenetic relationships of sugarcane rust fungi. Mycol. Prog. 9, 459–468. doi: 10.1007/s11557-009-0649-6
Elshire, R. J., Glaubitz, J. C., Sun, Q., Poland, J. A., Kawamoto, K., Buckler, E. S., et al. (2011). A robust, simple genotyping-by-sequencing (GBS) approach for high diversity species. PLoS One 6:e19379. doi: 10.1371/journal.pone.0019379
FAO (2013). FAOSTAT, FAO Statistical Databases. Available at: http://faostat3.fao.org/ [accessed December 18, 2015].
FAO (2014). FAOSTAT, FAO Statistical Databases. Available at: http://faostat3.fao.org/ [accessed June 07, 2016].
Flores, R. C., Loyo, J. R., Ojeda, R. A., Rangel, O., Cerón, F. A., Marquez, W., et al. (2009). First report of orange rust of sugarcane caused by Puccinia kuehnii in Mexico, El Salvador, and Panama. Plant Dis. 93:1347. doi: 10.1094/PDIS-93-12-1347B
Fu, L., Niu, B., Zhu, Z., Wu, S., and Li, W. (2012). CD-HIT: accelerated for clustering the next-generation sequencing data. Bioinformatics 28, 3150–3152. doi: 10.1093/bioinformatics/bts565
Garcés, F. F., Fiallos, F. F., Silva, E., Martinez, F., Aime, M. C., Comstock, J. C., et al. (2014). First report of orange rust of sugarcane caused by Puccinia kuehnii in Ecuador. Plant Dis. 98, 842. doi: 10.1094/PDIS-05-13-0574-PDN
Glynn, N. C., Laborde, C., Davidson, R. W., Irey, M. S., Glaz, B., D’Hont, A., et al. (2013). Utilization of a major brown rust resistance gene in sugarcane breeding. Mol. Breed. 31, 323–331. doi: 10.1007/s00122-017-2968-3
Harrell, F. E. Jr. (2018). Hmisc: Harrell Miscellaneous. Available at: http://cran.R-project.Org/web/packages/Hmisc
Jiang, G. H., Xu, C. G., Tu, J. M., Li, X. H., He, Y. Q., and Zhang, Q. F. (2004). Pyramiding of insect-and disease-resistance genes into an elite indica, cytoplasm male sterile restorer line of rice,‘Minghui 63’. Plant Breed. 123, 112–116. doi: 10.1046/j.1439-0523.2003.00917.x
Jones, P., Binns, D., Chang, H. Y., Fraser, M., Li, W., McAnulla, C., et al. (2014). InterProScan 5: genome-scale protein function classification. Bioinformatics 30, 1236–1240. doi: 10.1093/bioinformatics/btu031
Kelly, J. D., Afanador, L., and Haley, S. D. (1995). Pyramiding genes for resistance to bean common mosaic virus. Euphytica 82, 207–212. doi: 10.1007/s00122-011-1746-x
Kiraly, L., Barna, B., and Király, Z. (2007). Plant resistance to pathogen infection: forms and mechanisms of innate and acquired resistance. J. Phytopathol. 155, 385–396. doi: 10.1111/j.1439-0434.2007.01264.x
Li, T., Bai, G., Wu, S., and Gu, S. (2011). Quantitative trait loci for resistance to fusarium head blight in a Chinese wheat landrace Haiyanzhong. Theor. Appl. Genet. 122, 1497–1502. doi: 10.1007/s00122-011-1549-0
Li, W., and Godzik, A. (2006). Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics 22, 1658–1659. doi: 10.1093/bioinformatics/btl158
Lichtfouse, E., Navarrete, M., Debaeke, P., Souchère, V., Alberola, C., and Ménassieu, J. (2009). Agronomy for sustainable agriculture: a review. Agron. Sustain. Dev. 29, 1–6. doi: 10.1051/agro:2008054
Liu, J., Liu, D., Tao, W., Li, W., Wang, S., Chen, P., et al. (2000). Molecular marker-facilitated pyramiding of different genes for powdery mildew resistance in wheat. Plant Breed. 119, 21–24. doi: 10.1046/j.1439-0523.2000.00431.x
Magarey, R. C., Royal, A., Williams, D. J., and Bull, J. I. (2011). A brief history of disease epidemics in Queensland and of some economic outcomes. Proc. Aust. Soc. Sugarcane Technol. 33, 1–12.
McCouch, S. R. (2008). Gene nomenclature system for rice. Rice 1, 72–84. doi: 10.1007/s12284-008-9004-9
Ovalle, W., Comstock, J. C., Glynn, N. C., and Castlebury, L. A. (2008). First report of Puccinia kuehnii, causal agent of orange rust of sugarcane, in Guatemala. Plant Dis. 92:973. doi: 10.1094/PDIS-92-6-0973C
Paterson, A. H., Bowers, J. E., Bruggmann, R., Dubchak, I., Grimwood, J., Gundlach, H., et al. (2009). The Sorghum bicolor genome and the diversification of grasses. Nature 457, 551–556. doi: 10.1038/nature07723
Pinheiro, J., Bates, D., DebRoy, S., Sarkar, D., and R Core Team (2009). nlme: Linear and Nonlinear Mixed Effects Models. R Package Version. Available at: https://CRAN.R-project.org/package=nlme
R Development Core Team (2013). R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing.
Raboin, L., Oliveira, K. M., Lecunff, L., Telismart, H., Roques, D., Butterfield, M., et al. (2006). Genetic mapping in sugarcane, a high polyploid, using bi-parental progeny: identification of a gene controlling stalk colour and a new rust resistance gene. Theor. Appl. Genet. 112, 1382–1391. doi: 10.1007/s00122-006-0240-3
Rott, P. C., Sood, S., Comstock, J. C., Raid, R. N., Glynn, N. C., Gilbert, R. A., et al. (2014). Sugarcane Orange Rust. Available at: http://edis.ifas.ufl.edu/sc099
Sanseverino, W., Hermoso, A., D’Alessandro, R., Vlasova, A., Andolfo, G., Frusciante, L., et al. (2013). PRGdb 2.0: towards a community-based database model for the analysis of R-genes in plants. Nucleic Acids Res. 41, 1167. doi: 10.1093/nar/gks1183
Singh, S., Sidhu, J. S., Huang, N., Vikal, Y., Li, Z., Brar, D. S., et al. (2001). Pyramiding three bacterial blight resistance genes (xa5, xa13 and Xa21) using marker-assisted selection into indica rice cultivar PR106. Theor. Appl. Genet. 102, 1011–1015. doi: 10.1007/s001220000495
Sood, S. G., Comstock, J. C., and Glynn, N. C. (2009). Leaf whorl inoculation method for screening sugarcane rust resistance. Plant Dis. 93, 1335–1340. doi: 10.1094/PDIS-93-12-1335
Van Ooijen, J. W. (2006). JoinMap 4. Software for the Calculation of Genetic Linkage Maps in Experimental Populations. Wageningen: Kyazma BV.
Waclawovsky, A. J., Sato, P. M., Lembke, C. G., Moore, P. H., and Souza, G. M. (2010). Sugarcane for bioenergy production: an assessment of yield and regulation of sucrose content. Plant. Biotechnol. J. 8, 263–276. doi: 10.1111/j.1467-7652.2009.00491.x
Wang, S., Basten, C. J., and Zeng, Z. B. (2012). Windows QTL Cartographer 2.5. Raleigh, NC: North Carolina State University.
Wei, X., Jackson, P. A., McIntyre, C. L., Aitken, K. S., and Croft, B. (2006). Associations between DNA markers and resistance to diseases in sugarcane and effects of population substructure. Theor. Appl. Genet. 114, 155–164. doi: 10.1007/s00122-006-0418-8
Yang, X., Song, J., You, Q., Paudel, D. R., Zhang, J., and Wang, J. (2017a). Mining sequence variations in representative polyploid sugarcane germplasm accessions. BMC Genomics 18:594. doi: 10.1186/s12864-017-3980-3
Yang, X., Sood, S., Glynn, N., Islam, M. S., Comstock, J., and Wang, J. (2017b). Constructing high-density genetic maps for polyploid sugarcane (Saccharum spp.) and identifying quantitative trait loci controlling brown rust resistance. Mol Breed. 37, 116. doi: 10.1007/s11032-017-0716-7
Yang, X., and Wang, J. (2016). Genome-wide analysis of NBS-LRR genes in sorghum genome revealed several events contributing to NBS-LRR gene evolution in grass species. Evol. Bioinform. Online 12, 9–21. doi: 10.4137/EBO.S36433
Keywords: orange rust disease, sugarcane, Saccharum spp., quantitative trait locus (QTL), marker-assisted selection (MAS)
Citation: Yang X, Islam MS, Sood S, Maya S, Hanson EA, Comstock J and Wang J (2018) Identifying Quantitative Trait Loci (QTLs) and Developing Diagnostic Markers Linked to Orange Rust Resistance in Sugarcane (Saccharum spp.). Front. Plant Sci. 9:350. doi: 10.3389/fpls.2018.00350
Received: 25 September 2017; Accepted: 02 March 2018;
Published: 19 March 2018.
Edited by:Piergiorgio Stevanato, Università degli Studi di Padova, Italy
Reviewed by:Hongjun Liu, Shandong Agricultural University, China
Christopher John Lambrides, The University of Queensland, Australia
Copyright © 2018 Yang, Islam, Sood, Maya, Hanson, Comstock and Wang. 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) and the copyright owner 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: Jianping Wang, firstname.lastname@example.org