Abstract
Powdery mildew (PM), caused by Podosphaera xanthii (Px), is one of the most devastating fungal diseases of melon worldwide. The use of resistant cultivars is considered to be the best and most effective approach to control this disease. In this study, an F2 segregating population derived from a cross between a resistant (wm-6) and a susceptible cultivar (12D-1) of melon was used to map major powdery mildew resistance genes using bulked segregant analysis (BSA), in combination with next-generation sequencing (NGS). A novel quantitative trait locus (QTL) named qCmPMR-12 for resistance to PM on chromosome 12 was identified, which ranged from 22.0 Mb to 22.9 Mb. RNA-Seq analysis indicated that the MELO3C002434 gene encoding an ankyrin repeat-containing protein was considered to be the most likely candidate gene that was associated with resistance to PM. Moreover, 15 polymorphic SNPs around the target area were successfully converted to Kompetitive Allele-Specific PCR (KASP) markers (P < 0.0001). The novel QTL and candidate gene identified from this study provide insights into the genetic mechanism of PM resistance in melon, and the tightly linked KASP markers developed in this research can be used for marker-assisted selection (MAS) to improve powdery mildew resistance in melon breeding programs.
Introduction
Melon (Cucumis melo L.), which belongs to the Cucurbitaceae family, is an important horticultural and economic crop worldwide (). In view of its unique biological characteristics, such as a sweet fruit, unique aroma and rich nutritional value, melon is highly favored by consumers. Powdery mildew (PM) is a fungal disease that occurs commonly on leaves, petioles, and stems of most cucurbit crops in both field and greenhouse conditions (). This disease can result in a decrease in weight-based productivity and a reduction in fruit quality (), thereby causing severe economic losses in many areas of the world ().
Melon PM is often caused by Podosphaera xanthii (Px) and Golovinomyces cichoracearum (Gc) (; ). These two fungi can be distinguished by observation of the pathogens’ cleistothecia, conidia germination, microscopic morphology, and host identification (; ). Both pathogens exist either as an asexual stage (conidia) or sexual stage (Cleistothecia); the morphological characteristics of the asexual stage are similar for the two pathogens, and the morphological distinction between them is based on whether fibrous bodies occur on the spore (Px spores possess fibrous bodies, whereas these are absent from those of Gc) (). In China, Px is generally considered to be the main causal agent of PM in melon (; Zhang et al., 2013). Recent reports have shown that melon PM is caused by Px in Shanghai (; ). The five phases of the Px life cycle are separately germinating conidia, the formation of a primary germ tube, hyphae, conidiophores, and colonization (). After a spore lands on the leaf cuticle, the interaction is initiated by the formation of a primary germ tube and is followed by the formation and maturation of an appressiorial germ tube; the next step involves haustorium formation within an attacked epidermal cell and fungal development proceeds via the formation of secondary hyphae and haustoria and terminates with sporulation (). Px can coexist with the diseased plants in the soil as hyphae, conidia and cleistothecia, or it can overwinter as hyphae on greenhouse plants, and when the temperature rises the subsequent year, the conidia are distributed by air or water to infect the plants (Zhang et al., 2007).
Currently, the main method applied to control powdery mildew is chemical prevention by the application of chemical fungicides, which is often time-consuming and labor-intensive (). Furthermore, chemical control not only leads to the appearance of resistance and to mutations in PM, which therefore means that this treatment becomes ineffective, but it also increases the input costs and has a negative impact on the environment (; ). Therefore, breeding for disease resistance is a safe alternative or a complement to the chemical control of this disease (). Marker-assisted selection (MAS) is a powerful genomic tool that assists phenotypic selection for the development of disease-resistant cultivars and can help breeders incorporate and pyramid resistance genes into breeding material, thereby reducing disease severity (; ; Zhu et al., 2019). At present, MAS has been extensively applied to search for the molecular markers that are linked to a specific trait during the development of disease-resistant cultivars (). To date, a variety of molecular markers have been developed, such as RFLPs, RAPDs, SSRs, InDels and SNPs, to detect allelic variation within different samples at the DNA level ().
Bulked segregant analysis (BSA) is an important technique used to map quantitative trait loci (QTLs) and identify DNA markers. Compared with the traditional QTL mapping method, which is time-consuming and involves screening polymorphic markers and genotyping, BSA provides a convenient and rapid method with which to identify resistance genes by generating two DNA bulks with a contrasting target trait (; ; ). Recently, due to the release of sequenced genomes and the significant reduction in the costs of next-generation sequencing (NGS), whole-genome resequencing has been coupled with BSA to map the genes of interest that are associated with a given phenotype. The combined application of BSA with NGS (BSA-Seq) has accelerated the identification of tightly linked markers for gene identification and QTL mapping (Zou et al., 2016). To date, BSA-Seq has been successfully used in mapping the traits of early flowering, flesh thickness and downy mildew resistance in cucumber (; ; ), cold tolerance and blast resistance in rice (Yang et al., 2013; Zheng et al., 2016), cotyledon color, and a high-sucrose and low-oil seed phenotype in soybean (; ).
At present, the availability of sequence information has facilitated the identification and development of single nucleotide polymorphism (SNP) markers, which have largely replaced simple sequence repeats (SSRs) as markers in crop species (). Because of the low assay cost, high genomic abundance, ease of documentation, locus specificity, co-dominant inheritance, the potential for high-throughput analysis, and relatively low genotyping error rates, the use of SNPs has emerged as a powerful approach for many genetic applications in areas such as germplasm characterization, quality control (QC) analysis, linkage mapping, linkage-based and linkage disequilibrium-based QTL mapping, allele mining, marker-assisted backcrossing (MABC), genomic selection (GS), and MAS (; ; ). Kompetitive Allele-Specific PCR (KASP) is a high-throughput SNP genotyping platform. Due to its low cost and genotyping error rates, and its high reliability and reproducibility, KASP has evolved to become a global benchmark technology and has been widely used for genetic mapping and trait-specific marker development (; ; ; ).
To date, several genes and QTLs that confer resistance to powdery mildew have been identified in melon, such as the genes of Pm-w from WMR 29 (), Pm-x from PI 414723 and Pm-y from VA 435 (), and Pm-1 from the AF125Pm–1 Cantalupensis Charentais-type breeding line (), and the QTLs of PmV.1 and PmXII.1 from PI 124112 (), Pm-R from TGR-1551 (Yustelisbona et al., 2011) and BPm12.1 from MR-1 (). In previous studies, many of these genes and QTLs have been found on chromosomes 2, 4, 5, and 12 (; ; ; Zhang et al., 2013; ). Differing views exist concerning the genetic basis of PM resistance in melon. Some studies have indicated that PM resistance in melon is controlled by a single dominant gene (; Zhang et al., 2008; ; ), whereas other research has reported that it is controlled by a recessive gene (), by two dominant genes (), or by one dominant and one recessive gene (; Yuste-Lisbona et al., 2010). Moreover, it is also reported that resistance to PM in melon is controlled by different sets of QTL ().
The melon cultivars wm-6 and 12D-1 are both high-generation inbred lines developed by our group and our previous study has shown that the melon cultivar wm-6 is highly resistant to PM, whereas 12D-1 is highly susceptible (Data not shown). In this study, we obtained an F2 population from a cross between wm-6 (female) and 12D-1 (male). A major QTL that confers PM resistance on chromosome 12 was identified by BSA-Seq analysis, and a most likely candidate gene was predicted from RNA-Seq data in wm-6 melon. In addition, 15 suitable KASP markers were developed by the KASP SNP genotyping method. This will facilitate the cloning and functional validation of the candidate resistance gene and the linked markers will further provide a useful tool for MAS in melon breeding programs.
Materials and Methods
Plant Materials, Growth Conditions and Inoculation With Powdery Mildew Fungus
Two inbred lines, wm-6 (P1, resistant to PM) and 12D-1 (P2, susceptible to PM), were used as parental lines to generate F1 and F2 populations for the QTL mapping of PM resistance in melon (Cucumis melo L.). The P1, P2, F1 and F2 individuals were all placed in a culture room at a temperature of 25/20°C (day/night) with a photoperiod of 14 h light and relative humidity of 50–75%. All the seeds used in this study were provided by the Shanghai Academy of Agricultural Sciences.
The PM fungus (P. xanthii) used in this study was isolated from leaves of diseased Cucurbitaceae plants according to the method of . The plants were grown on the experimental farm of the Shanghai Academy of Agricultural Sciences, and the PM fungus was maintained by infection of susceptible melon cultivar plants. When the three true seedling leaves of melon plants were fully expanded, the fungus was collected and suspended in sterile distilled water containing 0.01% Tween 20 and was then used to inoculate plants at a concentration of 1 × 106 as previously described (Zhang et al., 2011).
Disease Evaluation for Resistance to Powdery Mildew
Phenotyping for powdery mildew resistance of melon was performed according to Zhang et al. (2013) with some modifications at 12-d post-inoculation (dpi), and each infected leaf was analyzed individually. Briefly, the disease grade of powdery mildew was categorized on a scale of 0–5 as follows: Class 0, no infection; Class 1, infection of less than 30% of the leaf; Class 3, infection of less than 70% of the leaf; Class 5, infection of approximately the entire leaf and coverage with heavy sporulation. The disease severity index (DSI) was calculated from the disease-rating scale using the following formula: DSI = 100 × Σ[(disease grade × number of plants in that grade)/(Total number of plants × maximum disease grade)]. For the inheritance study, lines with a DSI ≤ 20 were considered resistant and lines with a DSI > 20 were considered susceptible. For the F2 population, the same protocol was followed to identify resistant and susceptible plants. The phenotype of the two parental lines was analyzed at least three independent experiments with >15 seedlings examined in each experiment. The DSI of the two parental lines and the F1 hybrid plants was separately measured three times for 20 seedlings in each measurement. The DSI of F2 generation was calculated based on the phenotype of 193 F2 plants. On the basis of the DSI scores, plants with a DSI of 0–1 were categorized as resistant, and those with a DSI of 3–5, as susceptible.
RNA Extraction, Library Construction and Sequencing
For RNA-Seq analysis, the two parental lines wm-6 and 12D-1 were separately treated with water (mock) or PM fungus. After 3 days, the leaves of wm-6 and 12D-1 plants (named wm-6K, wm-6P, 12D-1K, 12D-1P, respectively) were harvested, immediately frozen in liquid nitrogen and stored at −80°C before RNA extraction.
Total RNA was extracted using the Trizol reagent kit (Invitrogen, Carlsbad, CA, United States) according to the manufacturer’s protocol. RNA quality was assessed on an Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, United States) and analyzed using RNase-free agarose gel electrophoresis. Following total RNA extraction, eukaryotic mRNA was enriched by Oligo(dT) beads, whereas prokaryotic mRNA was enriched by removing rRNA with the Ribo-ZeroTM Magnetic Kit (Epicentre, Madison, WI, United States). The enriched mRNA was fragmented into short fragments using fragmentation buffer and reverse transcribed into cDNA with random primers. Second-strand cDNA was synthesized by DNA polymerase I, RNase H, dNTPs and buffer. The cDNA fragments were purified with a QiaQuick PCR extraction kit (Qiagen, Venlo, The Netherlands), and following end repair and the addition of poly(A), was ligated to Illumina sequencing adapters. The ligation products were selected according to size by agarose gel electrophoresis, amplified by PCR, and sequenced using Illumina HiSeq2500 by the Gene Denovo Biotechnology Company (Guangzhou, China).
Genomic DNA Extraction, Library Construction for Bulked Segregant Analysis and Whole-Genome Resequencing
Young leaves from the two parental lines, and from the F1 and the F2 populations were collected, and total genomic DNA was extracted using the CTAB method (). For bulked segregant analysis, four DNA pools were constructed, consisting of two parent bulks and two F2 segregating bulks. The parent bulks were separately constructed from the female parent (wm-6) and male parent (12D-1), and the two F2 segregating bulks were separately constructed by mixing an equal amount of DNA extracted from 25 extremely resistant (R-bulk) and 25 susceptible (S-bulk) F2 plants. After the four sequencing libraries were prepared according to the standard protocol of Illumina, they were sequenced on an Illumina HisSeq platform (Illumina, San Diego, CA, United States). Short reads obtained from the four bulks were aligned against the melon reference genome sequence to obtain the consensus sequence using BWA software. SNP calling was performed with GATK tools. The heterozygous alleles in both parents were filtered out during the process. The raw sequence reads are deposited in the NCBI Sequence Read Archive (SRA; Accession number: PRJNA655764).
Genetic Mapping
In this study, four methods (SNP-ratio, ED4, G value and LOD) were used to map QTLs that underlay resistance to PM. The SNP-ratio (resistant alleles/sensitive alleles) of the R-bulk and S-bulk were calculated as described by and the SNP-ratio of the R-bulk was then divided by the SNP-ratio of the S-bulk and plotted across the genomic regions that showed ratio peaks, which indicate the possible existence of the QTLs. The read depth for each allele at segregating allelic SNPs in 500-kb sliding windows was summed using a 100-kb step increment. The Euclidean distance of each SNP (ED-SNP) was calculated as described by and the ED was the sum of 100 ED-SNP values within a window of 100 consecutive SNPs. The ED4 was calculated by raising ED to the fourth power. The G value averaged across neighboring SNPs was calculated according to . The LOD (logarithm of the odds) score was calculated as described by Zhang et al. (2019).
Development and Analysis of KASP Markers
The polymorphic SNPs identified around the target regions that associated with powdery mildew resistance were converted into KASP markers using PolyMarker software1. For each SNP, two allele-specific forward primers and one common reverse primer based on the flanking sequences around the variant position (SNP) were designed using Primer 3 software. The polymorphic SNP primers were converted to KASP markers to test their ability to differentiate the polymorphism by genotyping the two parents, and the KASP markers were then verified with the entire F2 population. Each KASP reaction was carried out using a 3-μL reaction mixture consisting of 1.48 μL KASP 2 × reaction mix, 50 ng DNA template, 0.17 μM Hex forward primer, 0.17 μM FAM forward primer and 0.42 μM universal reverse primer. The cycling conditions were as follows: 94°C for 15 min followed by 10 touchdown cycles at 95°C for 20 s and 65°C for 60 s (dropping 0.8°C per cycle); after the final annealing temperature of 56°C was reached, 26 cycles were performed at 94°C for 20 s and at 57°C for 60 s. Thermocycling and fluorescence readings were performed on a Hydrocycler and PHERAstar of LGC SNPline platform. Genotyping data were viewed as a cluster plot by SNPviewer software supported from LGC Genomics2. The significance of the correlation coefficients between phenotype and genotype was determined with t-tests. Linkage groups were constructed using JoinMap 4.1.
Results
Evaluation of Resistance to Powdery Mildew in wm-6 × 12D-1
The artificial inoculation results showed that the parental line wm-6 was highly resistant to P. xanthii (DSI = 8.0), whereas the other parental line 12D-1 (DSI = 94.0) was susceptible to the fungus (Figure 1, Table 1 and Supplementary Table 1). The DSI of the F1 plants was 68.0, which meant that the F1 generation was susceptible to PM fungus (Table 1 and Supplementary Table 1). The evaluation of infection by PM indicated that 58 F2 plants showed PM resistance and 135 showed PM susceptibility, with a 1:3 segregation between resistant and susceptible individuals (χ2 = 2.63, P = 0.11) (Table 1). This indicated that resistance to P. xanthii in wm-6 was conferred by a single recessive gene.
FIGURE 1
TABLE 1
| Cross | Generation | Inoculation time | Disease rating1 | Total | DSI | |||
| 0 | 1 | 3 | 5 | |||||
| wm-6 | P1 | 2019 | 12 | 8 | 0 | 0 | 20 | 8 |
| 12D-1 | P2 | 2019 | 0 | 0 | 3 | 17 | 20 | 94 |
| wm-6 × 12D-1 | F1 | 2019 | 0 | 5 | 6 | 9 | 20 | 68 |
| wm-6 × 12D-1 | F2 | 2019 | 26 | 32 | 38 | 97 | 193 | 65.39 |
Disease evaluation of powdery mildew in wm-6, 12D-1, F1 and F2 plants at 12 days post inoculation with P. xanthii.
The disease severity index of the two parental lines and the F1 hybrid plants were separately measured with 20 seedlings. The measurements were repeated three times independently with similar results, and the data of one representative experiment are shown. 1The DSI of leaves was rated on a 0–5 scale to determine the response of melon genotypes to powdery mildew, where 0 = immune, no symptom; 1 = highly resistant, infection of less than 30% of the leaf with low sporulation; 3 = moderately susceptible, infection of less than 70% of the leaf with moderate to high sporulation, and 5 = highly susceptible, infection of approximately the entire leaf with heavy sporulation.
Sequencing Data Analysis of Four DNA Bulks
BSA-Seq analysis was performed with the DNAs of four libraries (wm-6, 12D-1, R-bulk and S-bulk) using the Illumina HiSeq 2500 platform. In total, 60,401,327 and 57,755,318 clean reads for wm-6 and 12D-1, respectively, and 92,327,580 and 101,809,441 short reads for the R-bulk and S-bulk libraries, respectively, were generated. The GC content ranged from 36.76% to 37.1% and the Q20 and Q30 of each pool were over 98% and 93%, respectively. Over 95% of the reads were mapped to the melon reference genome, and the coverage rates were 95.78% in wm-6, 95.85% in 12D-1, 97.02% in R-bulk, and 97.05% in S-bulk, approximately resulting in a 23 × coverage depth for parental bulks and at least 35 × coverage for two F2 progeny bulks (Table 2). These results indicated that the quantity and quality of the data were sufficient for further analysis.
TABLE 2
| Sample name1 | Clean reads2 | Clean reads (%) | GC content (%) | Q20 (%)3 | Q30 (%)4 | Mapped reads (%)5 | Coverage (%)6 | Sequencing depth (×) |
| wm-6 | 60401327 | 95.86 | 36.80 | 98.02 | 93.17 | 95.80 | 95.78 | 23.34 |
| 12D-1 | 57755318 | 95.86 | 36.76 | 98.06 | 93.29 | 97.86 | 95.85 | 22.81 |
| R-bulk | 92327580 | 96.3 | 37.02 | 98.21 | 93.71 | 96.55 | 97.02 | 35.54 |
| S-bulk | 101809441 | 95.98 | 37.1 | 98.10 | 93.37 | 95.95 | 97.05 | 38.31 |
Summary of sequencing data and the data aligned to the melon reference genome for the parental lines and the resistant and susceptible pools by BSA-Seq.
1wm-6, powdery mildew-resistant parent; 12D-1, powdery mildew-susceptible parent; R-bulk, powdery mildew resistance bulk; S-bulk, powdery mildew susceptible bulk. 2Number of reads after trimming and adapter removal. 3The percentage of base recognition accuracy above 99%. 4The percentage of base recognition accuracy above 99.9%. 5Alignment to the melon genome assembly V4.0 (www.melonomics.net/). 6Coverage (≥1 read).
QTL Mapping of PM Resistance
After alignment to the reference genome of melon, 2,624,079 SNPs were identified between the parental lines wm-6 and 12D-1, and the distribution of these SNPs on each chromosome is listed in Supplementary Table 2. The resistance-related candidate regions were identified by analyzing the resistant pool and the susceptible pool. Using four statistic methods, a 0.9-Mb region spanning 22.0–22.9 Mb on chromosome 12 was defined as the target region associated with PM resistance (Figures 2A–D), and the significant QTL was designated qCmPMR-12. Within the candidate region, a total of 4,033 SNPs showed polymorphisms, 115 of which were located within the 3′ UTR, 78 in the 5′ UTR, 547 in the downstream region, 397 in the exonic region, 2,243 in the intergenic region, 1,005 in the intronic region, and 735 in the upstream region (Table 3). In addition, the target region included 182 non-synonymous SNV, 3 stopgain, 1 stoploss and 211 synonymous SNV polymorphisms (Table 4). Moreover, the candidate region contained 476 small InDels, 29 of which were located within the 3′ UTR, 23 in the 5′ UTR, 199 in the downstream region, 28 in the exonic region, 892 in the intergenic region, 387 in the intronic region, and 298 in the upstream region (Table 3). The InDels included 9 frameshift deletions, 6 frameshift insertions, 8 non-frameshift deletions and 5 non-frameshift insertions (Table 4).
FIGURE 2
TABLE 3
| Chromosomal location | 3′ UTR | 5′ UTR | Downstream | Exonic region | Intergenic region | Intronic region | Upstream |
| SNP number | 115 | 78 | 547 | 397 | 2243 | 1005 | 735 |
| InDel number | 29 | 23 | 199 | 28 | 892 | 387 | 298 |
The chromosomal locations and numbers of SNPs and InDels in the target QTL region associated with powdery mildew resistance between the melon parental lines wm-6 and 12D-1.
TABLE 4
| SNP type | SNP number | InDel type | InDel number |
| non-synonymous SNV1 | 182 | frameshift deletion | 9 |
| stopgain | 3 | frameshift insertion | 6 |
| stoploss | 1 | non-frameshift deletion | 8 |
| synonymous SNV1 | 211 | non-frameshift insertion | 5 |
The types and numbers of SNPs and InDels in the target QTL region associated with powdery mildew resistance between the melon parental lines wm-6 and 12D-1.
1Single nucleotide variants.
According to the melon gene annotation database3, the target region contains approximately 126 annotated genes, and the specific information for these genes is listed in Supplementary Table 3.
RNA-Seq and Candidate Gene Analysis
To elucidate the changes in gene expression of resistant and susceptible melon lines after inoculation with P. xanthii, we performed an RNA-Seq analysis at 3 dpi. The dataset submitted to NCBI include the raw reads of the assembled transcriptome sequences from control and pathogen-treated melon plants. All transcriptome raw reads have been deposited in NCBI SRA4 under the accession number PRJNA670091.
Analysis of the RNA-Seq results for the 126 annotated genes obtained by BSA-Seq indicated that the transcript levels of the two genes MELO3C002434 and MELO3C002477 (| log2FC| > 1) were clearly up-regulated in the parental lines after inoculation with P. xanthii5. The two genes respectively encode an ankyrin repeat-containing (ANK) protein and a homeobox-leucine zipper protein (Supplementary Table 3). Further analysis showed that only the expression level of MELO3C002434 was significantly higher in the resistant line wm-6 than in the susceptible line 12D-1 following PM infection. Notably, BSA-Seq showed that this gene contained 10 SNPs and 1 InDel (Table 5). Previous studies have shown that ANK proteins play important roles in regulating immune responses against various pathogens (; ; ; ). Therefore, we speculated that MELO3C002434 may confer PM resistance, or at least be a PM resistance-related gene in melon.
TABLE 5
| Type | Position | Reference | Allele | wm-6 | 12D-1 |
| SNP | 22665882 | C | T | 1/1 | 0/0 |
| SNP | 22666490 | A | G | 1/1 | 0/0 |
| SNP | 22666491 | G | A | 1/1 | 0/0 |
| SNP | 22666786 | A | G | 0/0 | 1/1 |
| SNP | 22666869 | G | A | 0/0 | 1/1 |
| SNP | 22667081 | C | G | 0/0 | 1/1 |
| SNP | 22667159 | T | C | 0/0 | 1/1 |
| SNP | 22668006 | A | C | 0/0 | 1/1 |
| SNP | 22668100 | T | C | 0/0 | 1/1 |
| SNP | 22668300 | T | C | 0/0 | 1/1 |
| InDel | 22666219 | C | CATT | 1/1 | 0/0 |
Analysis of SNPs and InDels present within the MELO3C002434 gene in two parental melon lines.
0/0 means that the genotype of the inbred was the same as the reference, and 1/1 means that the genotype of the inbred was the same as the allele.
KASP Marker Development and Physical Map Construction
KASP assays were designed for SNPs across the 20.0–23.9 Mb region on chromosome 12 and were tested on the F2 population to determine which markers showed the highest association with PM resistance. In this region, thirty-seven chromosome-specific SNPs selected for conversion to KASP markers were used to screen the parents and bulks to confirm their polymorphisms, and 15 out of 37 markers successfully distinguished the parents and bulks (Figure 3), indicating that these KASP markers were suitable for use in MAS to improve the level of powdery mildew resistance in melon breeding. Specific information and the sequences of the polymorphic KASP markers are separately listed in Supplementary Tables 4, 5. The close-up view of QTL and KASP markers linked to the target regions is shown in Figure 4.
FIGURE 3
FIGURE 4
Using the 15 KASP markers, the genetic physical map and linkage map were constructed by JoinMap 4.1 (Figures 5A,B). As shown in the linkage map, the peak of the QTL was located in a 0.6-cM interval spanned by KASP markers KA002213 and KA002215 (Figure 5B).
FIGURE 5
Discussion
The melon cultivar wm-6 is a high-generation inbred line developed by our group, which is highly resistant to powdery mildew (Figure 1A, Table 1 and Supplementary Table 1). In this study, a genomic region ranging from 22.0 Mb to 22.9 Mb on chromosome 12 was identified using BSA-Seq technology. Similarly, detected a major QTL for PM resistance on chromosome 12 and suggested that resistance to P. xanthii in MR-1 was controlled by a single dominant gene. However, in this study, the segregation ratio of the F2 population indicated that the PM resistance in wm-6 was controlled by a single recessive gene, implying a different genetic basis of the PM resistance mechanisms between the melon cultivars wm-6 and MR-1. Other studies have also demonstrated that the genetic basis of resistance to PM differs depending on the tested melon material (; Zhang et al., 2008; ; Yuste-Lisbona et al., 2010; ; ).
To identify the PM resistance gene in wm-6, RNA-Seq analysis was performed with the melon lines wm-6 and 12D-1 in parallel. Combination of the results from BSA-Seq and RNA-Seq suggested that the At3g12360-like gene MELO3C002434, which encodes an ANK protein was the most likely candidate gene to confer PM resistance, because only this gene was significantly more highly expressed in the resistant line wm-6 than that in the susceptible line 12D-1 following infection with P. xanthii (see text footnote 5). Previous studies revealed that ANK proteins have critical functions in various biological processes of plant growth and development as well as in response to biotic and abiotic stresses (; ; ; ; ; ; ; ). It has been reported that the ANK protein NPR1 is important both in the SA-dependent immune response and in SA-independent resistance responses induced by the root-associated bacteria (; ); AKR2 functions in the oxidative metabolism of disease resistance and stress response in Arabidopsis (); a plasma membrane-localized ANK protein, ACD6, is involved in SA-dependent signaling in defense responses and programmed cell death (, ), and ectopic expression of the ANK protein OsBIANK1 of rice confers enhanced disease resistance to Botrytis cinerea and Pseudomonas syringae in Arabidopsis (). Interestingly, the gene MELO3C002434 contained 10 SNPs and 1 InDel (Table 5) and 5 out of 15 KASP markers that developed in this study were located within this gene (Supplementary Table 5). All of these results further implied that MELO3C002434 is the most likely candidate gene associated with PM resistance in melon wm-6. In addition, the report of suggested MELO3C002434 to be one of the genes related to PM resistance in spite of the different genetic basis of the PM resistance mechanism in melon wm-6 and MR-1. Therefore, further genetic studies and more detailed analyses are required to confirm the role and molecular mechanism of action of MELO3C002434 in the PM defense response.
In this study, 37 SNPs surrounding the candidate region were used to design KASP markers using the PolyMarker website and 15 of these were polymorphic among the bulk and parent populations. To the best of our knowledge, this is the first report that KASP markers have been developed and used in MAS to improve PM resistance in melon breeding, although several markers linked to resistance genes have been reported previously in melon (; ).
Conclusion
In this study, a major QTL that is associated with PM resistance was identified in a 0.9-Mb interval on chromosome 12 of melon using BSA-Seq technology. Additional RNA-Seq data suggested that an ankyrin repeat-containing gene within this region, MELO3C002434, was implicated to be the most important candidate gene. In addition, 15 suitable KASP markers that were tightly linked to the resistance phenotype were developed for the MAS of melon. These data can be used to improve PM resistance in breeding programs and to facilitate understanding of the molecular mechanisms that underlying PM resistance in melon.
Statements
Data availability statement
The raw sequence reads produced by BSA-Seq are deposited in the NCBI Sequence Read Archive (SRA; Accession number: PRJNA655764).
Author contributions
YZ, YC, and HZ designed the research. YC and QD prepared the plant materials. YC performed the experiments, analyzed the data, and wrote the manuscript. YZ, YC, and HJ revised the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This study was supported by the Shanghai Science and Technology Commission (20392000300), Shanghai Minhang Science and Technology Commission (2019MHC051), Excellent Team of Shanghai Academy of Agricultural Sciences (B-06), and Shanghai Melon and Fruit Industrial Technology System [(2020)1].
Conflict of interest
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2020.593207/full#supplementary-material
Supplementary Table 1Means, standard errors and ranges of DSI of the parents and F1 plants at 12 dpi with P. xanthii in the year of 2019.
Supplementary Table 2The number of SNPs distributed on each chromosome of melon detected by BSA-Seq.
Supplementary Table 3The positions and annotations of candidate genes in the target region conferring powdery mildew resistance from 22.0 Mb to 22.9 Mb based on the melon reference genome (CM3.6.1).
Supplementary Table 4KASP primer sequences used for MAS in melon breeding programs.
Supplementary Table 5Genotypes and positions of the 15 KASP markers linked to powdery mildew resistance gene of melon.
Footnotes
1.^http://polymarker.tgac.ac.uk
3.^https://www.melonomics.net/
References
1
AbeA.KosugiS.YoshidaK.NatsumeS.TakagiH.KanzakiH.et al (2012). Genome sequencing reveals agronomically important loci in rice using MutMap.Nat. Biotechnol.30174–178. 10.1038/nbt.2095
2
AlbertS.DespresB.GuilleminotJ.BechtoldN.PelletierG.DelsenyM.et al (1999). The EMB 506 gene encodes a novel ankyrin repeat containing protein that is essential for the normal development of Arabidopsis embryos.Plant J.17169–179. 10.1046/j.1365-313x.1999.00361.x
3
CandidoV.CampanelliG.ViggianiG.LazzeriL.CameleI. (2014). Melon yield response to the control of powdery mildew by environmentally friendly substances.Entia Horticult.16670–77. 10.1016/j.scienta.2013.12.008
4
CaoH.GlazebrookJ.ClarkeJ. D.VolkoS.DongX. (1997). The Arabidopsis NPR1 gene that controls systemic acquired resistance encodes a novel protein containing ankyrin repeats.Cell8857–63. 10.1016/s0092-8674(00)81858-9
5
ChenX. M. (2013). High-temperature adult-plant resistance, key for sustainable control of stripe rust.Am. J. Plant Sci.04608–627. 10.4236/ajps.2013.43080
6
ChenX. M. (2014). Integration of cultivar resistance and fungicide application for control of wheat stripe rust.Can. J. Plant Pathol.36311–326. 10.1080/07060661.2014.924560
7
ChengH.KongW.SuY.LiuD.HeQ. (2011). Identification on powdery mildew and physiological race of melon in china.J. Changjiang Vegetables181–5. 10.29405/j.bes/411-104725
8
ClementsJ. (2014). Powdery mildew resistance gene of ssr molecular markers linked to the ‘wmr-29’ strains of melon.Northern Horticult.61347–357.
9
DobbelsA. A.MichnoJ. M.CampbellB. W.VirdiK. S.StecA. O.MuehlbauerG. J.et al (2017). An Induced Chromosomal Translocation in Soybean Disrupts a KASI Ortholog and Is Associated with a High-Sucrose and Low-Oil Seed Phenotype.G371215–1223. 10.1534/g3.116.038596
10
DoyleJ. (1991). DNA protocols for plants-CTAB total DNA isolation.Mole. Techniq. Taxonomy43:25.
11
EpinatC. M.PitratM.BertrandF. (1992). Genetic analysis of resistance of 5 melon lines to powdery mildews.Acta Horticult.65135–144. 10.1007/bf00022575
12
ErtiroB. T.OgugoV.WorkuM. (2015). Comparison of Kompetitive Allele Specific PCR (KASP) and genotyping by sequencing (GBS) for quality control analysis in maize.BMC Genomics16:908. 10.1186/s12864-015-2180-2
13
FukinoN.OharaT.MonforteA. J.SugiyamaM.SakataY.KunihisaM.et al (2008). Identification of QTLs for resistance to powdery mildew and SSR markers diagnostic for powdery mildew resistance genes in melon (Cucumis melo L.).Theoret. Appl. Genet.118165–175. 10.1007/s00122-008-0885-1
14
Garcia-MasJ.BenjakA.SanseverinoW.BourgeoisM.BaulcombeD. C. (2012). The genome of melon (Cucumis melo L.).Proc. Natl. Acad. Sci.10911872–11877.
15
GarcionC.GuilleminotJ.KrojT.ParcyF.GiraudatJ.DevicM. (2006). AKRP and EMB506 are two ankyrin repeat proteins essential for plastid differentiation and plant development in Arabidopsis.Plant J.48895–906. 10.1111/j.1365-313x.2006.02922.x
16
GuH. F.ZhangX.ZhangW. F.JinC. Y.XuL.LuX. L.et al (2010). Identification of physiological races of watermelon and melon powdery mildew in shanghai.Acta Agricult. Shanghai26155–158.
17
HaC. M.JunJ. H.NamH. G.FletcherJ. C. (2004). BLADE-ON-PETIOLE1 encodes a BTB/POZ domain protein required for leaf morphogenesis in Arabidopsis thaliana.Plant Cell Physiol.451361–1370. 10.1093/pcp/pch201
18
HanY.LvP.HouS.LiS.JiG.MaX.et al (2015). Combining next generation sequencing with bulked segregant analysis to fine map a stem moisture locus in sorghum (Sorghum bicolor L. Moench).PLoS One10:e0127065. 10.1371/journal.pone.0127065
19
HeC.HolmeJ.AnthonyJ. (2014). SNP genotyping: the KASP assay.Methods Mol. Biol.114575–86. 10.1007/978-1-4939-0446-4_7
20
HemsleyP. A.KempA. C.GriersonC. S. (2005). The TIP GROWTH DEFECTIVE1 S-acyl transferase regulates plant cell growth in Arabidopsis.Plant Cell172554–2563. 10.1105/tpc.105.031237
21
HillJ. T.DemarestB. L.BisgroveB. W.GorsiB.SuY. C.YostH. J. (2013). Mmappr: Mutation Mapping Analysis Pipeline for Pooled Rna-Seq.Genome Res.23687–697. 10.1101/gr.146936.112
22
HollomonD. W.WheelerI. E.BélangerR. R.BushnellW. R.DikA. J.CarverT. L. W. (2002). “The Powdery Mildews,” in A Comprehensive Treatise, edsBélangerR. R.BushnellW. R.DikA. J.CarverT. L. W. (Saint Paul, MN: APS Press), 249–255.
23
KřístkováE.LebedaA.SedlákováB. (2009). Species spectra, distribution and host range of cucurbit powdery mildews in the Czech Republic, and in some other European and Middle Eastern countries.Phytoparasitica37337–350. 10.1007/s12600-009-0045-4
24
LiB.ZhaoY.ZhuQ.ZhangZ.FanC.AmanullahS.et al (2017). Mapping of powdery mildew resistance genes in melon (Cucumis melo L.) by bulked segregant analysis.Sci. Horticult.220160–167. 10.1016/j.scienta.2017.04.001
25
LiD.WangF.LiuB.ZhangY.HuangL.ZhangH.et al (2013). Ectopic Expression of Rice OsBIANK1, Encoding an Ankyrin Repeat-Containing Protein, in Arabidopsis Confers Enhanced Disease Resistance to Botrytis cinerea and Pseudomonas syringae.J. Phytopathol.16127–34. 10.1111/jph.12023
26
LiP. F.ZhuL. L.YangX. P.XuJ. H.LiuG.YaoX. F.et al (2015). Species and physiological race identification of melon powdery mildew fungus in Jiangsu, Zhejiang and Shanghai areas.China Cucurbits Vegetables2816–20.
27
LiangQ. L.XuB. L.YanH. X. (2010). Hull-less pumpkin powdery mildew and its host range.Mycosystema29636–643.
28
LiuL. Z.ChenY. Y.SuZ. H.ZhangH.ZhuW. M. (2010). A sequence-amplified characterized region marker for a single, dominant gene in melon PI 134198 that confers resistance to a unique race of podosphaera xanthii in china.Hortscience451407–1410. 10.21273/hortsci.45.9.1407
29
LuH.LinT.KleinJ.WangS.QiJ.ZhouQ.et al (2014). QTL-seq identifies an early flowering QTL located near Flowering Locus T in cucumber.Theor. Appl. Genet.1271491–1499. 10.1007/s00122-014-2313-z
30
LuH.LiuY.GreenbergJ. T. (2005). Structure-function analysis of the plasma membrane-localized Arabidopsis defense component ACD6.Plant J.44798–809. 10.1111/j.1365-313x.2005.02567.x
31
LuH.RateD. N.SongJ. T.GreenbergJ. T. (2003). ACD6, a novel ankyrin protein, is a regulator and an effector of salicylic acid signaling in the Arabidopsis defense response.Plant Cell152408–2420. 10.1105/tpc.015412
32
MagweneP. M.WillisJ. H.KellyJ. K. (2011). Te Statistics of Bulk Segregant Analysis Using Next Generation Sequencing.PLoS Comput. Biol.7:e1002255. 10.1371/journal.pcbi.1002255
33
McCreightJ. D.CoffeyM. D. (2011). Inheritance of resistance in melon PI 313970 to cucurbit powdery mildew incited by podosphaera xanthii race s.HortScience46838–840. 10.21273/hortsci.46.6.838
34
McGrathT. M. (2001). Fungicide resistance in cucurbit powdery mildew: experiences and challenges.Plant Dis.85236–245. 10.1094/pdis.2001.85.3.236
35
MichelmoreR. W.ParanI.KesseliR. V. (1991). Identification of markers linked to disease resistance genes by bulked segregant analysis: a rapid method to detect markers in specific genomic regions by using segregating populations Proc.Natl. Acad. Sci. U S A889828–9832. 10.1073/pnas.88.21.9828
36
NieJ.HeH.PengJ.YangX.BieB.ZhaoJ.et al (2015). Identification and fine mapping ofpm5.1: a recessive gene for powdery mildew resistance in cucumber (cucumis sativus L.).Mole. Breed.35:7.
37
NingX.WangX.GaoX.ZhangZ.ZhangL.YanW.et al (2014). Inheritances and location of powdery mildew resistance gene in melon Edisto47.Euphytica1951–9.
38
PerchepiedL.BardinM.DogimontC.PitratM. (2005). Relationship Between Loci Conferring Downy Mildew and Powdery Mildew Resistance in Melon Assessed by Quantitative Trait Loci Mapping.Phytopathology95556–565. 10.1094/phyto-95-0556
39
PerchepiedL.MarcB.CatherineD.PitratM. (2005). Relationship between loci conferring downy mildew and powdery mildew resistance in melon assessed by quantitative trait loci mapping.Phytopathology95556–565. 10.1094/phyto-95-0556
40
Perez-GarciaA.RomeroD.Fernandez-OrtunoD.Lopez-RuizF.De VicenteA.ToresJ. A. (2009). The powdery mildew fungus Podosphaera fusca (synonym Podosphaera xanthii), a constant threat to cucurbits.Mole. Plant Pathol.10153–160.
41
PérinC.HagenL.ContoV. D.KatzirN.Danin-PolegY.PortnoyV.et al (2002). A reference map of Cucumis melo based on two recombinant inbred line populations.Theor. Appl. Genet1041017–1034. 10.1007/s00122-002-0864-x
42
PitratM. (1991). Linkage groups in Cucumis melo L.Rep. Cucurbit Genetics Cooperative82406–411. 10.1093/oxfordjournals.jhered.a111112
43
RafalskiA. (2002). Applications of single nucleotide polymorphisms in crop genetics.Curr. Opin. Plant Biol.594–100. 10.1016/s1369-5266(02)00240-6
44
RasheedA.WenW.GaoF.ZhaiS.JinH.LiuJ.et al (2016). Development and validation of KASP assays for genes underpinning key economic traits in bread wheat.Theoret. Appl. Genet.1291843–1860. 10.1007/s00122-016-2743-x
45
RibautJ. M.BanzigerM.BetranJ.JiangC.EdmeadesG. O.DreherK.et al (2002). “Use of molecular markers in plant breeding: Drought tolerance improvement in tropical maize,” in itor, Quantitative genetics, genomics, and plant breeding, ed.KangM. S. (Wallingford, UK: CABI Publishing).
46
RomeroD.RiveraM. E.CazorlaF. M.CodinaJ. C.Fernandez-OrtunoD.ToresJ. A.et al (2008). Comparative histochemical analyses of oxidative burst and cell wall reinforcement in compatible and incompatible melon-powdery mildew (Podosphaera fusca) interactions.J. Plant Physiol.1651895–1905. 10.1016/j.jplph.2008.04.020
47
RyalsJ.WeymannK.LawtonK. (1997). The Arabidopsis NIM1 protein shows homology to the mammalian transcription factor inhibitor I kappa B.Plant Cell9425–439. 10.1105/tpc.9.3.425
48
SakamotoH.MatsudaO.IbaK. (2008). ITN1, a novel gene encoding an ankyrin-repeat protein that affects the ABA-mediated production of reactive oxygen species and is involved in salt-stress tolerance in Arabidopsis thaliana.Plant J.56411–422. 10.1111/j.1365-313X.2008.03614.x
49
SchlottererC. (2004). The evolution of molecular markers-just a matter of fashion?Nat. Rev. Genet563–69. 10.1038/nrg1249
50
SemagnK.BabuR.HearneS.OlsenM. (2014). Single nucleotide polymorphism genotyping using Kompetitive Allele Specific PCR (KASP): overview of the technology and its application in crop improvement.Mole. Breed.331–14. 10.1007/s11032-013-9917-x
51
ShenG.KuppuS.VenkataramaniS.WangJ.YanJ.QiuX.et al (2010). ANKYRIN REPEAT-CONTAINING PROTEIN 2A is an essential molecular chaperone for peroxisomal membranebound ASCORBATE PEROXIDASE3 in Arabidopsis.Plant Cell22811–831. 10.1105/tpc.109.065979
52
SongJ.LiZ.LiuZ.GuoY.QiuL. (2017). Next-Generation Sequencing from Bulked-Segregant Analysis Accelerates the Simultaneous Identification of Two Qualitative Genes in Soybean.Front. Plant Sci.8919–919. 10.3389/fpls.2017.00919
53
SoykS.LemmonZ. H.OvedM.FisherJ.LiberatoreK. L.ParkS. J.et al (2017). Bypassing Negative Epistasis on Yield in Tomato Imposed by a Domestication Gene.Cell1691142–1155 e1112.
54
SunX. L.LiuD.ZhangH. Q.HuoN. X.ZhouR. H.JiaJ. Z. (2010). Identification and mapping of two new genes conferring resistance to powdery mildew from aegilops tauschii (coss.) schmal.J. Integr. Plant Biol.481204–1209. 10.1111/j.1744-7909.2006.00328.x
55
TanC.AssangaS.ZhangG.RuddJ. C.HaleyS. D.XueQ.et al (2017). Development and Validation of KASP Markers for Wheat Streak Mosaic Virus Resistance Gene Wsm2.Crop Sci.57340–349. 10.2135/cropsci2016.04.0234
56
TeixeiraA.MatosobarretoA. P.AranhaL. E. (2008). An AFLP marker linked to the Pm-1 gene that confers resistance to Podosphaera xanthii race 1 in Cucumis melo.Genet. Mol. Biol.31547–550. 10.1590/s1415-47572008000300023
57
WangX. Z.XuB. Y.WangP.LuanF. S.MaH. Y.MaY. Y. (2013). Identification of Powdery Mildew Pathogen and Ribosomal DNA-ITS Sequence Analysis on Melon.J. Northeast Agricult. Univ.1018–97.
58
WangY. H.WuD. H.HuangJ. H.TsaoS. J.HwuK. K.LoH. F. (2016). Mapping quantitative trait loci for fruit traits and powdery mildew resistance in melon (cucumis melo).Botanical. Stud.57:19.
59
WangY.NingZ.HuY.ChenJ.ZhaoR.ChenH.et al (2015). Molecular Mapping of Restriction-Site Associated DNA Markers in Allotetraploid Upland Cotton.PLoS One10:e0124781. 10.1371/journal.pone.0124781
60
WinK. T.VegasJ.ZhangC.SongK.LeeS. (2017). QTL mapping for downy mildew resistance in cucumber via bulked segregant analysis using next-generation sequencing and conventional methods.Theoret. Appl. Genet.130199–211. 10.1007/s00122-016-2806-z
61
WolterM.HollricherK.SalaminiF.Schulze-LefertP. (1993). The mlo resistance alleles to powdery mildew infection in barley trigger a developmentally controlled defence mimic phenotype.MGG239122–128. 10.1007/bf00281610
62
XuX.LuL.ZhuB.XuQ.QiX.ChenX. (2015). QTL mapping of cucumber fruit flesh thickness by slaf-seq.Scient. Rep.5:15829.
63
YanJ.WangJ.ZhangH. (2002). An ankyrin repeat-containing protein plays a role in both disease resistance and antioxidation metabolism.Plant J.29193–202. 10.1046/j.0960-7412.2001.01205.x
64
YangZ.HuangD.TangW.ZhengY.LiangK.CutlerA. J.et al (2013). Mapping of Quantitative Trait Loci Underlying Cold Tolerance in Rice Seedlings via High-Throughput Sequencing of Pooled Extremes.PLoS One8:e68433. 10.1371/journal.pone.0068433
65
YustelisbonaF. J.CapelC.SarriaE.TorreblancaR.GomezguillamonM. L.CapelJ.et al (2011). Genetic linkage map of melon (Cucumis melo L.) and localization of a major QTL for powdery mildew resistance.Mole. Breed.27181–192. 10.1007/s11032-010-9421-5
66
Yuste-LisbonaF. J.López-SeséA. I.ómez-GuillamónM. L. G. (2010). Inheritance of resistance to races 1, 2 and 5 of powdery mildew in the melon tgr-1551.Plant Breed.12972–75. 10.1111/j.1439-0523.2009.01655.x
67
ZhangC.RenY.GuoS. (2013). Application of comparative genomics in developing markers tightly linked to the Pm-2f gene for powdery mildew resistance in melon (Cucumis Melo L.).Euphytica1901–12. 10.12816/0011695
68
ZhangH. Y.SuF.GuoS. G.GongG. Y.YiH. P.WuM. Z.et al (2008). Genetic analysis and specific fragments linked to powdery mildew resistant gene Pm-2F in melon.Acta Horticult. Sin.121773–1780. 10.16420/j.issn.0513-353x.2008.12.009
69
ZhangH.GuoS.GongG.RenY.DavisA. R.XuY. (2011). Sources of resistance to race 2WF powdery mildew in U.S. watermelon plant introductions.Hortscience461349–1352. 10.21273/hortsci.46.10.1349
70
ZhangH.WangX.PanQ.LiP.LiuY.LuX.et al (2019). QTG-Seq Accelerates QTL Fine Mapping through QTL Partitioning and Whole-Genome Sequencing of Bulked Segregant Samples.Mol. Plant12426–437. 10.1016/j.molp.2018.12.018
71
ZhangS. X.CuiJ.SongY. Y.LiuS. Q. (2007). Occurrence and control of the melon cultivar.Decid. Fruits3941–42.
72
ZhangX. J.JiJ.ZhaiW. Q.LiM. H.WangH. J.YiH. P. (2013). Identification of physiological races of powdery mildew on Cucumis melo ssp. melo in xinjiang.Xinjiang Agricult. Sci.501450–1455.
73
ZhengW.WangY.WangL.MaZ.ZhaoJ.WangP.et al (2016). Genetic mapping and molecular marker development for Pi65(t), a novel broad-spectrum resistance gene to rice blast using next-generation sequencing.Theoret. Appl. Genet.1291035–1044. 10.1007/s00122-016-2681-7
74
ZhuH.ZhaiW.LiX.ZhuY. (2019). Two QTLs controlling clubroot resistance identified from bulked segregant sequencing in pakchoi (brassica campestris ssp. chinensis makino).Scient. Rep.9:9228.
75
ZouC.WangP.XuY. (2016). Bulked sample analysis in genetics, genomics and crop improvement.Plant Biotechnol. J.141941–1955. 10.1111/pbi.12559
Summary
Keywords
bulked segregant analysis, KASP markers, melon, powdery mildew, resistance gene
Citation
Cao Y, Diao Q, Chen Y, Jin H, Zhang Y and Zhang H (2021) Development of KASP Markers and Identification of a QTL Underlying Powdery Mildew Resistance in Melon (Cucumis melo L.) by Bulked Segregant Analysis and RNA-Seq. Front. Plant Sci. 11:593207. doi: 10.3389/fpls.2020.593207
Received
10 August 2020
Accepted
28 October 2020
Published
05 February 2021
Volume
11 - 2020
Edited by
Magdalena Arasimowicz-Jelonek, Adam Mickiewicz University, Poland
Reviewed by
Christina Cowger, Plant Science Research Unit (USDA-ARS), United States; Zhiyong Liu, Institute of Genetics and Developmental Biology, Chinese Academy of Sciences, China
Updates
Copyright
© 2021 Cao, Diao, Chen, Jin, Zhang and Zhang.
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(s) 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: Yongping Zhang, zyp123944@126.comHongmei Zhang, zhanghongmei@saas.sh.cn
This article was submitted to Plant Pathogen Interactions, a section of the journal Frontiers in Plant Science
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.