Detection of Blackleg Resistance Gene Rlm1 in Double-Low Rapeseed Accessions from Sichuan Province, by Kompetitive Allele-Specific PCR

Article information

Plant Pathol J. 2021;37(2):194-199
Publication date (electronic) : 2021 April 1
doi : https://doi.org/10.5423/PPJ.OA.10.2020.0204
1Crop Research Institute, Sichuan Academy of Agricultural Sciences, Chengdu, Sichuan 610066, China
2Department of Plant Science, University of Manitoba, Winnipeg, Manitoba R3T 2N2, Canada
3Institute of Plant Protection, Sichuan Academy of Agricultural Sciences, Chengdu, Sichuan 610066, China
*Co-corresponding authors: L. Jiang, Tel) +86-28-84504235, E-mail) jlcrape@163.com. W. G. Dilantha Fernando, Tel) 204-474-8577, E-mail) dilantha.fernando@umanitoba.ca
Received 2020 November 1; Revised 2021 January 26; Accepted 2021 February 1.

Abstract

Blackleg is a serious disease in Brassica plants, causing moderate to severe yield losses in rapeseed worldwide. Although China has not suffered from this disease yet (more aggressive Leptosphaeria maculans is not present yet), it is crucial to take provisions in breeding for disease resistance to have excellent blackleg-resistant cultivars already in the fields or in the breeding pipeline. The most efficient strategy for controlling this disease is breeding plants with identified resistance genes. We selected 135 rapeseed accessions in Sichuan, including 30 parental materials and 105 hybrids, and we determined their glucosinolate and erucic acid content and confirmed 17 double-low materials. A recently developed single-nucleotide polymorphism (SNP) marker, SNP_208, was used to genotype allelic Rlm1/rlm1 on chromosome A07, and 87 AvrLm1-resistant materials. Combined with the above-mentioned seed quality data, we identified 11 AvrLm1-resistant double-low rapeseed accessions, including nine parental materials and two hybrids. This study lays the foundation of specific R gene-oriented breeding, in the case that the aggressive Leptosphaeria maculans invades and establishes in China in the future and a robust and less labor consuming method to identify resistance in canola germplasm.

Rapeseed (Brassica napus) is the second most important oilseed crop globally, after soybean (Glycine max) (Basunanda et al., 2010; Sun et al., 2018). Apart from edible oil and animal feed, rapeseed also provides raw material for biodiesel or industrial erucic acid production, as well as the landscape for tourism.

Blackleg, caused by two fungal pathogens, aggressive Leptosphaeria maculans or less (or non-) aggressive L. biglobosa (Fernando et al., 2016), is a severe disease in rapeseed and a major constraint to its production for many producing areas in the world, including Canada, Australia, and Europe (Fitt et al., 2006; Fu et al., 2019; Howlett, 2004). Blackleg causes a global economic loss of over $900 million per growing season (Fitt et al., 2006, 2008). China, as another high rapeseed producing country, of which oilseed production accounts for approximately 30% of the total global yield, has not been plagued by blackleg to date (Liu et al., 2014), although there are few reports of non-aggressive L. biglobosa (Zhang et al., 2016b). However, as the international seed trades are becoming increasingly frequent, it is necessary for Chinese breeders to take provisions to safeguard the crop. To the best of our knowledge, blackleg is difficult to control with fungicides (Kutcher et al., 2010); thus, the most efficient and environmentally sound approach for controlling this disease is deployment of plant resistance genes prior to L. maculans incursions into China (Delourme et al., 2006; Zhang et al., 2014).

There are two types of resistance against blackleg in rapeseed: seedling resistance, which is a qualitative trait conferred by a major gene (R gene), and adult plant resistance (APR), which is a quantitative trait conferred by multiple minor genes (Zhang et al., 2016a). Researchers have suggested that the rotation of R genes could be an efficient strategy to control blackleg (Marcroft et al., 2012), and the foundation of this is to validate the presence of R genes in plants and to specify them. So far, at least 18 R genes in Brassica species have been identified (Zhang et al., 2016a); however, only three of them, i.e., Rlm2 (Larkan et al., 2015), LepR3 (Larkan et al., 2013), and Rlm9 (Larkan et al., 2020), were cloned. Rlm1 is a single dominant R gene (Fu et al., 2019), and its interaction with corresponding avirulence gene AvrLm1 in L. maculans provided the first genetic evidence for the “gene-to-gene” interactions (Ansan-Melayah et al., 1998). Fu et al. (2019) used BC1 and F2 populations and fine-mapped Rlm1 gene in rapeseed by using single-nucleotide polymorphism (SNP) markers through kompetitive allele-specific PCR (KASP), a high throughput genotyping technology based on fluorescence (Semagn et al., 2014); then, an SNP marker (SNP_208) was developed that could distinguish Rlm1 from allele rlm1 by KASP. This allele SNP site (G for resistance while T for susceptibility) was located at the position of 19,989,284 bp on chromosome A07 and resulted in a missense variant in candidate gene BnaA07g27460D. That research provides a foundation for genotyping the allele Rlm1/rml1 in canola and related germplasm.

Double-low (low erucic acid and glucosinolates content) rapeseed oil is considered healthier for humans. In the last few decades, great gains in double-low rapeseed breeding have been achieved in China. In this study, we selected 135 rapeseed accessions, including 30 parental materials and 105 hybrids, to determine their glucosinolate and erucic acid contents; then, we performed KASP with the marker SNP_208 mentioned above to genotype these alleles from these rapeseed materials, screening out materials with a resistant Rlm1 locus. Thus, we identified many double-low and AvrLm1-resistant germplasm that would be potential sources of germplasm to move forward in our breeding program.

Plant materials

The 135 rapeseed germplasms, consisting of 30 parental materials and 105 hybrids, were all bred by and kept in the Crop Research Institute, Sichuan Academy of Agricultural Sciences. None of the hybrids were derived from these parental plants. The seeds were sowed in September 2018, in experimental fields in the Xindu District (30°47′10″N, 104°12′12″E) in the Sichuan Basin, with an annual average temperature of 16.2°C in humid subtropical monsoon climate. Each accession was planted in row of 8 plants, with 30 cm between plants within each row and between rows. Three biological replications were applied.

Seed quality determination

The seed quality characters, oil content, erucic acid content, and glucosinolate content were determined by near infrared spectroscopy (FOSS model 5000, NIR System Inc., Silver Spring, MD, USA), with the room temperature at 20–25°C and air humidity under 60%. The parameters were described by Gan et al. (2003). Three biological repetitions were applied and standard deviation was calculated by MS Excel. “Double-low” was defined as erucic acid content <2% of oil content and glucosinolate content <30 μmol/g for parental plants or <40 μmol/g for hybrids. For each accession, a specimen-cup of seeds (usually 6–7 g) were used for each time.

Genotyping of rapeseed germplasm

Leaves from each accession were sampled at the seedling stage. The total DNA was extracted by the cetyltrimethyl ammonium bromide method (Doyle and Doyle, 1990). The DNA quality was validated using a NanoDrop 2000 spectrophotometer (Thermo Scientific, Wilmington, DE, USA). Each DNA sample was then diluted to 20 ng/μl.

The 4 μl KASP reaction volume contained 2 μl of 2× Taq DNA polymerase mix, 1 μl of SNP primer mix (Table 1), and 4 ng DNA. KASP comprised a touchdown profile consisting of an initial denaturation at 94°C for 10 min; followed by 94°C for 20 s, 61°C (with a 0.6°C drop per cycle) for 45 s, and 72°C for 30 s for 10 cycles; then, 37 cycles of 94°C for 20 s, 55°C for 45 s, and 72°C for 30 s; with a final extension of 72°C for 10 min. The PCR was performed on a Bio-Rad MyCycler (Hercules, CA, USA), and then the fluorescence information was collected by KASP Mix (LGC, type 1536) with Omega F plate reader (LGC, Teddington, UK). The KlusterCaller software (LGC) was used to read the fluorescence information and analyze it.

Primers of marker SNP_208 for KASP

The analysis of the physical and chemical properties for the protein of candidate gene was performed using DNA-MAN 6.0 (Lynnon Biosoft, Quebec, Canada).

Oil content, erucic acid, and glucosinolate content in the seed

We measured three important characteristics of the seeds from the 135 accessions: the oil content, erucic acid, and glucosinolate content. The oil content of the seeds from 135 accessions varied from 38.29% to 55.61% (Fig. 1A), and pin-130 represented the highest value (Supplementary Table 1). The 42–45% group showed the highest number of accessions, at 60; 12 accessions had oil content over 50%. The erucic acid content varied significantly among accessions. Accession pin-112 showed the lowest erucic acid content at only 0.39%, while the highest value was 45.62% in accession pin-121 (Supplementary Table 1). In total, 42 accessions represented erucic acid contents lower than 2% (Fig. 1B), which was the defined as “low erucic acid content”. We found that 22 accessions had glucosinolate contents lower than 30 μmol/g, while 10 accessions represented contents ranging 30 to 40 μmol/g (Fig. 1C). The “low glucosinolate content” for rapeseed is usually defined as lower than 30 μmol/g in the seeds of parental plants and lower than 40 μmol/g for hybrids; the definition of “low erucic acid content” adopts the same standard for both parental plants and hybrids, at 2% of the oil content. According to this, fourteen parental plants (pin-101, pin-102, pin-103, pin-104, pin-105, pin-109, pin-111, pin-113, pin-115, pin-116, pin-128, pin-129, pin-130, and pin-131) and three hybrids (pin-10, pin-80, and pin-119) were identified as “double-low” (Table 2). In China, for edible oil, only the registered double-low rapeseed are currently allowed to be commercially released; while the hybrids, usually F1 hybrids derived from specific parental materials, are candidates for commercial use, which requires quality testing before registration.

Fig. 1

Distribution of seed oil quality in the 150 selected accessions. (A) The oil content in seed of 150 accessions. (B) Erucic acid content in seed of 150 accessions. (C) Glucosinolate content in seed of 150 accessions.

The characters of the seeds from the 17 double-low accessions

Genotyping

The allelic SNP was genotyped using KASP, and the positive samples were categorized mainly into three clusters (Supplementary Table 2): homozygous resistance (G/G), heterozygous (G/T), and homozygous susceptible (T/T). Among the 30 parental plants, 18 (60%) were identified as resistant, 8 were identified as susceptible, and the other 4 were negative, which meant it was neither G nor T on that particular locus; additionally, all 18 resistant samples were homozygous. In the 105 hybrids, 69 materials (65.7%) were clustered into a resistance group, including 52 homozygous samples and 17 heterozygous ones. We found 29 samples identified as susceptible, and the other 7 samples were negative (Table 3, Supplementary Table 2). This Rlm1-containing frequency seems higher than earlier reported by Zhang et al. (2016b). This is to be expected, as they used germplasm from mainly Anhui Province; herein, the materials tested were to breed double-low rapeseed suitable for Sichuan climate. It is common to have different parental germplasm in each breeding program, and China’s breeding material is no different from other countries in that. Some of them are progenies of common main parents (also known as backbone parental plants). This is a long-term selection process, so it is unavoidable to result in similar genetic background to some extent; thus, within such a range, it can result in higher Rlm1-containing frequency, which is good news for Sichuan’s rapeseed/canola breeders.

Number of lines for each genotype at the Rlm1 locus

Combined with the above-mentioned seed quality data, nine double-low parental materials (pin-101, pin-105, pin-109, pin-111, pin-115, pin-116, pin-128, pin-129, and pin-130) as well as two double-low hybrids (pin-80 and pin-119) were identified harboring allelic Rlm1. Rlm1 conferred resistance to AvrLm1 of L. maculans. As it is a dominant allele, it is beneficial and practical for hybrid utilization, because as long as one of the parental plants is homozygous in the Rlm1/Rlm1 site, the hybrid seeds are consequently resistant. Therefore, validating the presence of Rlm1 in parental material facilitates the resistant hybrid breeding.

The overuse or consecutive use of one single R gene tends to cause resistance breakdown, such as the ‘sylvestris’ resistance breakdown in Australia in 2003 (Van de Wouw et al., 2010). One solution to this is R gene rotation, which should be built on identification of these R genes. This study provides a foundation for that.

Notes

Handling Editor: Cecile Segonzac

Conflicts of Interest

No potential conflict of interest relevant to this article was reported.

Acknowledgements

This work was supported by the National Key Research and Development Plan- the International Cooperation Project for China-Canada (2018YFE0108200); Modern Agroindustry Technology Research System of China (CARS-12); Research and Public Service on Germplasm Resources of Main Crops.

Electronic Supplementary Material

Supplementary materials are available at The Plant Pathology Journal website (http://www.ppjonline.org/).

Supplementary Table 1.

Oil content, erucic acid and glucosinolate contents in the seeds of all the 135 accessions

Supplementary Table 2.

The allelic SNP information from the 135 rapeseed accessions

References

Ansan-Melayah D., Balesdent M. H., Delourme R., Pilet M. L., Tanguy X., Renard M., Rouxel T.. 1998;Genes for race-specific resistance against blackleg disease in Brassica napus L. Plant Breed 117:373–378.
Basunanda P., Radoev M., Ecke W., Friedt W., Becker H. C., Snowdon R. J.. 2010;Comparative mapping of quantitative trait loci involved in heterosis for seedling and yield traits in oilseed rape (Brassica napus L.). Theor. Appl. Genet 120:271–281.
Delourme R., Chèvre A. M., Brun H., Rouxel T., Balesdent M. H., Dias J. S., Salisbury P., Renard M., Rimmer S. R.. 2006;Major gene and polygenic resistance to Leptosphaeria maculans in oilseed rape (Brassica napus). Eur. J. Plant Pathol 114:41–52.
Doyle J. J., Doyle J. L.. 1990;Isolation of plant DNA from fresh tissues. Focus 12:13–15.
Fernando W. G. D., Zhang X., Amarasinghe C. C.. 2016;Detection of Leptosphaeria maculans and Leptosphaeria biglobosa causing blackleg disease in canola from Canadian canola seed lots and dockage. Plants 5:12.
Fitt B. D. L., Brun H., Barbetti M. J., Rimmer S. R.. 2006;World-wide importance of phoma stem canker (Leptosphaeria maculans and L. biglobosa) on oilseed rape (Brassica napus). Eur. J. Plant Pathol 114:3–15.
Fitt B. D. L., Hu B. C., Li Z. Q., Liu S. Y., Lange R. M., Kharbanda P. D., Butterworth M. H., White R. P.. 2008;Strategies to prevent spread of Leptosphaeria maculans (phoma stem canker) onto oilseed rape crops in China; costs and benefits. Plant Pathol 57:652–664.
Fu F., Liu X., Wang R., Zhai C., Peng G., Yu F., Fernando W. G. D.. 2019;Fine mapping of Brassica napus blackleg resistance gene Rlm1 through bulked segregant RNA sequencing. Sci. Rep 9:14600.
Gan L., Sun X., Jin L., Wang G., Xu J., Wei Z., Fu T.. 2003;Establishment of math models of NIRS analysis for oil and protein contents in seed of Brassica napus. Sci Agric Sin 36:1609–1613. (in Chinese).
Howlett B. J.. 2004;Current knowledge of the interaction between Brassica napus and Leptosphaeria maculans. Can. J. Plant Pathol 26:245–252.
Kutcher H. R., Yu F., Brun H.. 2010;Improving blackleg disease management of Brassica napus from knowledge of genetic interactions with Leptosphaeria maculans. Can. J. Plant Pathol 32:29–34.
Larkan N. J., Lydiate D. J., Parkin I. A. P., Nelson M. N., Epp D. J., Cowling W. A., Rimmer S. R., Borhan M. H.. 2013;The Brassica napus blackleg resistance gene LepR3 encodes a receptor-like protein triggered by the Leptosphaeria maculans effector AVRLM1. New Phytol 197:595–605.
Larkan N. J., Ma L., Borhan M. H.. 2015;The Brassica napus receptor-like protein RLM2 is encoded by a second allele of the LepR3/Rlm2 blackleg resistance locus. Plant Biotechnol J 13:983–992.
Larkan N. J., Ma L., Haddadi P., Buchwaldt M., Parkin I. A. P., Djavaheri M., Borhan M. H.. 2020;The Brassica napus wall-associated kinase-like (WAKL) gene Rlm9 provides race-specific blackleg resistance. Plant J 104:892–900.
Liu Z., Latunde-Dada A. O., Hall A. M., Fitt B. D. L.. 2014;Phoma stem canker disease on oilseed rape (Brassica napus) in China is caused by Leptosphaeria biglobosa ‘brassicae’. Eur. J. Plant Pathol 140:841–857.
Marcroft S. J., Van de Wouw A. P., Salisbury P. A., Potter T. D., Howlett B. J.. 2012;Effect of rotation of canola (Brassica napus) cultivars with different complements of blackleg resistance genes on disease severity. Plant Pathol 61:934–944.
Semagn K., Babu R., Hearne S., Olsen M.. 2014;Single nucleotide polymorphism genotyping using kompetitive allele specific PCR (KASP): overview of the technology and its application in crop improvement. Mol. Breed 33:1–14.
Sun L., Wang X., Yu K., Li W., Peng Q., Chen F., Zhang W., Fu S., Xiong D., Chu P., Guan R., Zhang J.. 2018;Mapping of QTLs controlling seed weight and seed-shape traits in Brassica napus L. using a high-density SNP map. Euphytica 214:228.
Van de Wouw A. P., Cozijnsen A. J., Hane J. K., Brunner P. C., McDonald B. A., Oliver R. P., Howlett B. J.. 2010;Evolution of linked avirulence effectors in Leptosphaeria maculans is affected by genomic environment and exposure to resistance genes in host plants. PLoS Pathog 6:e1001180.
Zhang X., Peng G., Kutcher H. R., Balesdent M.-H., Delourme R., Fernando W. G. D.. 2016a;Breakdown of Rlm3 resistance in the Brassica napus–Leptosphaeria maculans pathosystem in western Canada. Eur. J. Plant Pathol 145:659–674.
Zhang X., Peng G., Parks P., Hu B., Li Q., Jiang L., Niu Y., Huang S., Fernando W. G. D.. 2016b;Identifying seedling and adult plant resistance of Chinese Brassica napus germplasm to Leptosphaeria maculans. Plant Pathol 66:752–762.
Zhang X., White R. P., Demir E., Jedryczka M., Lange R. M., Islam M., Li Z. Q., Huang Y. J., Hall A. M., Zhou G., Wang Z., Cai X., Skelsey P., Fitt B. D. L.. 2014;Leptosphaeria spp., phoma stem canker and potential spread of L. maculans on oilseed rape crops in China. Plant Pathol 63:598–612.

Article information Continued

Fig. 1

Distribution of seed oil quality in the 150 selected accessions. (A) The oil content in seed of 150 accessions. (B) Erucic acid content in seed of 150 accessions. (C) Glucosinolate content in seed of 150 accessions.

Table 1

Primers of marker SNP_208 for KASP

Primer Sequence
VIC 5′-GAAGGTCGGAGTCAACGGATTCGTTAACAAACTCTGTCGCG-3′
FAM 5′-GAAGGTGACCAAGTTCATGCTCCGTTAACAAACTCTGTCGCT-3′
COM 5′-GGCGTCGATCAGGTGAGCC-3′

This marker is slight modified according to Fu et al. (2019).

SNP, single-nucleotide polymorphism; KASP, kompetitive allele specific PCR; FAM and VIC, two SNP alleles labelled with fluorescent tags Fam and Vic, respectively; COM, common reverse primer.

Table 2

The characters of the seeds from the 17 double-low accessions

Accession type Sample ID Oil content in seed (% of seed weight) Erucic acid content (% of oil content) Glucosinolate content (μmol/g in seed)
Parental plant pin-101 48.24 ± 1.28 1.39 ± 0.31 26.13 ± 1.17
pin-102 46.64 ± 1.24 1.18 ± 0.41 24.96 ± 1.05
pin-103 45.8 ± 2.14 1.65 ± 0.78 26.92 ± 2.01
pin-104 49.1 ± 1.06 0.55 ± 0.26 28.84 ± 1.49
pin-105 47.8 ± 0.93 1.35 ± 0.53 27.11 ± 1.06
pin-109 44.95 ± 2.17 1.39 ± 0.49 23.49 ± 3.49
pin-111 44.61 ± 0.59 0.88 ± 0.66 22.75 ± 0.53
pin-113 48.64 ± 1.63 1.48 ± 0.76 29.93 ± 1.66
pin-115 49.15 ± 0.82 0.65 ± 0.49 27.85 ± 1.8
pin-116 47.13 ± 2.58 1.68 ± 1.61 23.26 ± 0.91
pin-128 49.89 ± 0.29 0.68 ± 0.43 27.85 ± 2.29
pin-129 52.9 ± 1.26 0.42 ± 0.42 28.39 ± 0.97
pin-130 55.61 ± 1.15 0.65 ± 0.65 27.21 ± 3.03
pin-131 52.82 ± 1.5 1.36 ± 1.36 25.21 ± 1.17
Hybrid pin-10 41.87 ± 1.57 0.79 ± 0.12 39.48 ± 1.9
pin-80 42.03 ± 1.3 1.5 ± 0.73 39.97 ± 0.73
pin-119 49.04 ± 0.47 1.37 ± 0.7 28.13 ± 2.22

Values are presented as mean ± SE.

Table 3

Number of lines for each genotype at the Rlm1 locus

G/G G/T T/T Negative
Parental plant 18 0 8 4
Hybrid 52 17 29 7

Negative means it was neither G nor T on this locus.