- School of Horticulture and Plant Protection, Yangzhou University, Yangzhou, China
Fusarium wilt (FW) is a very serious soil-borne disease worldwide, which usually results in huge yield losses in cucumber production. However, the inheritance and molecular mechanism of the response to FW are still unknown in cucumber (Cucumis sativus L.). In this study, two inbred cucumber lines Superina (P1) and Rijiecheng (P2) were used as the sensitive and resistant lines, respectively. A mixed major gene plus polygene inheritance model was used to analyze the resistance to FW in different generations of cucumber, namely, P1, P2, F1 (P1×P2), B1, and B2, obtained by backcrossing F1 plants with Superina (B1) or Rijiecheng (B2), and F2, obtained by self-crossing the F1 plants. After screening 18 genetic models, we chose the E-1 model, which included two pairs of additive-dominance-epistatic major genes and additive-dominance polygenes, as the optimal model for resistance to FW on the basis of fitness tests. The major effect quantitative trait locus (QTL) fw2.1 was detected in a 1.91-Mb-long region of chromosome 2 by bulked-segregant analysis. We used five insertion/deletion markers to fine-map the fw2.1 to a 0.60 Mb interval from 1,248,093 to 1,817,308 bp on chromosome 2 that contained 80 candidate genes. We also used the transcriptome data of Rijiecheng inoculated with Fusarium oxysporum f. sp. cucumerinum (Foc) to screen the candidate genes. Twelve differentially expressed genes were detected in fw2.1, and five of them were significantly induced by FW. The expression levels of the five genes were higher in FW-resistant Rijiecheng inoculated with Foc than in the control inoculated with water. Our results will contribute to a better understanding of the genetic basis of FW resistance in cucumber, which may help in breeding FW-resistant cucumber lines in the future.
Introduction
Cucumber (Cucumis sativus L.) is an economically important vegetable crop that ranks fourth in vegetable production worldwide. In 2017, the cultivated area in China was 1.24 million hectares with an annual output of 64.88 million tons, respectively accounting for 54.4% and 77.5% of the world total 2017 (http://faostat3.fao.org/). In the cucumber production season, Fusarium wilt (FW), downy mildew, and powdery mildew are the three main diseases found in China (Cao et al., 2007). FW of cucumber occurs more readily and seriously under a continuous cropping system, with incidences ranging from 30% to 90%, leading to huge loss of cucumber output (Pu et al., 2011; Zhou and Wu, 2012).
FW caused by Fusarium oxysporum f. sp. cucumerinum Owen (Owen, 1955), which is a forma specialis that infects the vascular bundles of cucumber, leads to necrotic lesions on the stem base, foliar wilting and eventually whole plant wilt, and even death (Zhou et al., 2008; Ye et al., 2004). Until now, there are no methods to effectively control the occurrence and harm of cucumber FW (Ye et al., 2004; Yuan et al., 2003).
Understanding the inheritance of FW resistance is an important step in developing resistant breeding resources and in breeding resistant varieties. The inheritance of cucumber resistance to FW has been studied for a long time, but the results are inconsistent (Pierce and Wehner, 1990; Vakalounakis, 1993). Vakalounakis, (2015) found that cucumber resistance to FW was controlled by a single gene, whereas others have suggested that cucumber FW resistance was regulated by multiple genes (Liu et al., 2003; Wang, 2005). To date, the inheritance of cucumber FW resistance remains poorly understood. Disease-resistant traits are quantitative traits that are controlled by multiple genes, which are generally located in multiple effect quantitative trait loci (QTLs; Hiroki et al., 2013; Falconer and Mackay, 1996). Bulked-segregant analysis (BSA) is an effective method for identifying DNA markers tightly linked to FW resistance (Giovannoni et al., 1991; Michelmore et al., 1991). For example, Li et al. (2018) mapped powdery mildew resistance genes to a region of chromosome 1A of wheat, and Zhang et al. (2018) confirmed that resistance to downy mildew and powdery mildew shared common candidate intervals on chromosome 5 of cucumber (Zhang et al., 2018) by BSA. DNA samples from soybean with opposite phenotypes were subjected to BSA to detect DNA markers that exhibited differences between the two different samples to identify the QTL (Mansur et al., 1993). Genome-wide insertion/deletion (InDel) markers have been used for fine mapping of important economical traits in rice (Feng et al., 2005; Li et al., 2017), wheat (Chen et al., 2005; Shang, 2009), and tomato (Su et al., 2019). In this study, a mixed major gene plus polygene inheritance model was used to analyze FW resistance in cucumber, and the major effect QTL of FW was investigated by BSA. We combined transcriptome data and mapping data to detect the key candidate resistance genes in cucumber after inoculation with Foc.
Materials and Methods
Cucumber Materials and Treatments
The two cucumber inbred lines, Superina and Rijiecheng, that showed susceptibility and resistance to FW, respectively, were used in this study. Superina (female parent, P1) and Rijiecheng (male parent, P2) were used to construct different generations, namely, F1, B1 (F1×Superina), B2 (F1×Rijiecheng), and F2. The F1 generation (P1×P2) was planted in Autumn 2016 in a greenhouse at the experimental farm of the Department of Horticulture in Yangzhou University, and the B1, B2, and F2 generations were obtained by backcrossing the F1 plants with Superina (B1) or Rijiecheng (B2) and self-crossing the F1 plants (F2) in Spring 2017. The 200 P1, 200 P2, 200 F1, 200 B1, 200 B2, and 1500 F2 plants were planted in 36-cavity plates filled with aseptic organic substrates (N, P, and K = 40–60 g/kg, humus ≥350 g/kg, pH = 6.5–7.5) in the greenhouse in Autumn 2017. Another 500 F2 and 200 F2 were planted in Autumn 2018 and Spring 2019, respectively. Seedlings at the second true leaf stage were inoculated with the Foc suspension (concentration 106 conidia/ml) using the root irrigate method. The disease classification scale was as follows: grade 0, asymptomatic; grade 1, slightly discolored stem base and cotyledon; grade 2, necrotic patches at the stem base and slight wilting of the seedling; grade 3, color of necrotic patches at the stem was deep with longitudinal fissure and visible wilting of the seedling; and grade 4, the seedling was dead (modified from Zhou et al., 2010). The disease index was calculated as follows:
Phenotypic Statistics and Experimental Design
Ten days after Foc inoculation, the disease grade of the seedlings was recorded and analyzed. The extremely resistant and extremely sensitive seedlings were separately pooled into an R-pool and S-pool, respectively. Genomic regions that showed signatures of resistance to FW were detected by whole genome resequencing of the DNA from the two parents and the two pools. Regions that underwent specific selection in the opposite directions were selected (Hiroki et al., 2013). The reads in the R- and S-pools were mapped to the Cucumber (Chinese Long) v2 Genome (http://cucurbitgenomics.org/organism/2). The single nucleotide polymorphisms (SNP)-index and Euclidean distance of the two pools were calculated and compared. The Euclidean distance was estimated as described by Hill et al. (2013). The SNP-indexes of the two pools were different because of genotype selection and knock-on effects. QTLs related to the FW resistance trait were roughly located by taking the intersection of the above two results (SNP-index and Euclidean distance).
Genetic Analysis
The FW resistance in the six generations of cucumber (P1, P2, F1, B1, B2, and F2) was analyzed using the mixed major gene plus polygene inheritance model (Gai and Wang, 1998). The maximum log-likelihood value and Akaike information criterion (AIC) were obtained by estimating the parameters of each generation and component distributions using the iterated EMC (IEMC) algorithm (Gai and Wang, 1998). After repeated iteration, the algorithm converged to a relatively stable and consistent result. The optimal model was selected according to the AIC, and the corresponding component distribution parameters were obtained. The first- and second-order parameters of the optimal model were estimated using a least squares method.
DNA Extraction, Identification of InDel Markers, and Genotyping
Total genomic DNA was extracted from cotyledons of the parent seedlings using the Cetyltrimethyl Ammonium Bromide (CTAB) -acidic phenol extraction method (Xu et al., 2017). The concentration and quality of the extracted DNA were determined by ultraviolet spectrophotometer (Thermo Fisher, USA) and 1% agarose gel electrophoresis. The primer pairs used for screening for InDel markers were designed using Primer Premier 5.0 with the Cucumber v2 Genome sequence as the reference. For each sample, the PCR mixture (20 µl total volume) contained 2 µl 10× buffer, 1 µl dNTPs (10 mM), 1 µl primer F (50 ng/µl), 1 µl primer R (50 ng/µl), 1 µl DNA, 0.4 µl Taq DNA polymerase (10 U/µl), and 13.6 µl diethyl pyrocarbonate water (DEPC water). The touchdown PCR amplifications were performed using an Eppendorf Mastercycler Pro (Eppendorf). Subsequently, 1 µl PCR product was detected and analyzed by polyacrylamide gel electrophoresis.
PCR and polyacrylamide gel electrophoresis were used to identify polymorphisms in the F2 generation using the parent Superina and Rijiecheng plants, and the F1 plants are the reference. By combining the genotype and phenotype of the F2 seedlings, we determined the location of the interval associated with resistance on the cucumber genome sequence. We designed new InDel marker primers in the resistance interval to fine-map the interval associated with resistance and repeated this until the interval was small or until it contained no new markers.
Validation of Gene Expression by Real-Time PCR (qPCR)
We have transcriptome data for Rijiecheng, a relatively resistant line, at 0, 24, 48, 96, and 192 h after inoculation with Foc (unpublished data). We used the Foc-inoculated and water-inoculated Rijiecheng seedlings to verify the expression levels of candidate genes by qPCR. Total RNA of each sample was isolated using a Mini BEST Plant RNA Extraction Kit (TaKaRa, China) and then dissolved in Ultra Pure™ DNase/RNase-free distilled water (Invitrogen, USA). The total RNA was reverse-transcribed using a PrimeScript™ RT reagent kit with genomic DNA (gDNA) eraser (TaKaRa, China). Primer sequences were designed using Beacon Designer 7.0 and screened using SeqHunter 1.0 The qPCRs were performed using SYBR® Premix Ex Taq™ (TaKaRa, China), according to the manufacturer’s instructions. SYBR Green PCR cycling was performed on an IQTM5 Multicolor Real-Time PCR Detection System (Bio-Rad, USA) using 20 µl samples as follows: 95 for 3 min, followed by 39 cycles of 95℃ for 10 s, 60 for 20 s, and 72 for 20 s. The relative quantization of gene expression was calculated and normalized to tubulin alpha chain (Csa4G000580). Three biological replicates from each condition were used for qPCRs.
Results
Variations in FW Phenotypes Among the Cucumber Parents and Segregated Populations
Ten days after Foc inoculation, the disease symptom grades and disease indexes of the seedlings of different generation were recorded and calculated. The disease indexes of P1, P2, F1, B1, B2, and F2 were 73.21, 18.13, 51.43, 48.24, 37.68, and 51.71, respectively (Figure 1). The average disease index of F1 plants was higher than the mid parent value, indicating that the resistance traits of F1 plants tended to be from the male parent Rijiecheng. The disease indexes of the two backcross generations (B1 and B2) shifted towards the backcross parent. The F2 population showed positive correlation between the distribution of disease grades and obvious quantitative genetic characteristics using three replicates (Figure 2), which corresponds to the genetic characteristics of a mixed major gene plus polygene inheritance model.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/478559/fpls-10-01425-HTML-r1/image_m/fpls-10-01425-g001.jpg)
Figure 1 Disease index of each generation. F2 plants were obtained from the cross between Superina (female parent, P1) and Rijiecheng (male parent, P2). Superina (female parent, P1) and Rijiecheng (male parent, P2) were used to construct different generations, namely F1, B1 (F1× Superina), B2 (F1× Rijiecheng), and F2.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/478559/fpls-10-01425-HTML-r1/image_m/fpls-10-01425-g002.jpg)
Figure 2 Frequency distribution of the disease grades among the F2 population. The disease grade was recorded 10 days after inoculation in Autumn 2017 and 2018, and Spring 2019.
Best-Fitting Genetic Model for FW Resistance and the Affecting Factors
To obtain a genetic model of cucumber FW resistance, we analyzed the segregation of the resistant phenotype among the six generations using the mixed major gene plus polygene inheritance model. We analyzed and screened a total of 18 genetic models including five categories by combining the maximum log-likelihood value and AIC. The five categories were as follows: one major gene model, two major genes model, polygene model, one major gene plus polygene model, and two major genes plus polygene model (Table 1). We selected the two major genes plus polygene models (E, E-1, E-2, and E-3) that had the smallest AIC values as the candidate models. Fitness tests of these four models, including equal distribution (U12, U22, and U32), Smirnov (nW2), and Kolmogorov tests (Dn) (Xu et al., 2017), indicated the significance levels for the E, E-1, E-2, and E-3 models were 7, 1, 11, and 11, respectively (Table 2). Finally, we selected the E-1 model (two pairs of additive-dominance-epistatic major genes and additive-dominance polygenes) as the optimal model for resistance to FW on the basis of the combined AIC values and goodness-of-fit test.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/478559/fpls-10-01425-HTML-r1/image_m/fpls-10-01425-t001.jpg)
Table 1 Maximum log-likelihood values (MLVs) and Akaike information criterion (AIC) values under various genetic models estimated using the IECM algorithm.
To determine the genetic effects of major genes, we estimated first- and second-order distribution parameters for resistance to FW in the F1 generation using the optimal E-1 model (Table 3). The additive effect values of the major gene and polygene were both 0.05, and the dominant effect values of the major gene and polygene were 0.99 and 0.57, respectively, which indicated that the additive effect of the major gene was consistent with that of the polygene, whereas the dominant effect of the major gene was greater than that of the polygene. This suggested that the contribution rates of major gene and polygene were consistent with inheritance of resistance to FW. The heritability of major gene from B1, B2, and F2 generations was 22.91%, 29.70%, and 59.73%; the heritability of polygene from B1, B2, and F2 was 29.92%, 21.78%, and 2.43%; and the environmental variance accounted for 47.62%, 48.52%5, and 37.84% of the phenotypic variance for B1, B2 and F2, respectively. These results indicated that environmental factors had a large effect on cucumber resistance to FW.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/478559/fpls-10-01425-HTML-r1/image_m/fpls-10-01425-t003.jpg)
Table 3 Estimates of genetic parameters for resistance to FW in the F1 plants (Superina×Rijecheng) cross under the E-1 model
Identification of the Major Effect FW QTL by BSA
A total of 214,270,191 clean reads from the transcriptomes of Superina, Rijiecheng, S-pool, and R-pool were aligned to the Cucumber v2 Genome sequence; the Q30-level was >95%, and the GC-contents was >35%. The detailed statistics for each transcriptome including numbers of clean reads, numbers of bases, and average depth are given in Table 4. We identified a 3.78 Mb region on chromosome 2 (Chr2) that contained 625 genes that had an association threshold (Euclidean distance) of 0.12. One region on Chr2 with a total length of 1.91 Mb and containing 319 genes had an association threshold of 0.25 after fitting the △SNP-index. We used the intersection of these two results to identify the region of the genome associated FW resistance. The identified region was on Chr2, was 1.91 Mb long, and contained 319 genes (Figure 3), and was designated as the major effect QTL (fw2.1) related to the FW resistance trait.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/478559/fpls-10-01425-HTML-r1/image_m/fpls-10-01425-g003.jpg)
Figure 3 Identification of the major effect Fusarium wilt (FW) quantitative trait locus (QTL). Colored points indicate the Euclidean distance values of every SNP on each of the chromosomes. The black lines indicate the Euclidean distance values after matching. The red dotted line is the threshold line. Values above the threshold line were used to select the candidate FW trait-related interval as shown in the red box.
Fine Mapping QTL fw2.1 to 0.6 Mb on Chr2
To narrow down the major effect QTL identified by BSA, whole genome resequencing of the two parents was performed to confirm the high-quality InDels. We genotyped 500 F2 plants using two polymorphic InDel markers (InDel89749 and InDel1817308). Three other InDel markers that were evenly distributed in the major effect QTL interval were designed. Six different genotypes were detected among the F2 plants; however, two of them had disease grades 0 or 2 (Figure 4). By combining the genotypes and phenotypes of the F2 seedlings, the candidate interval was located between InDel1248093 and InDel1817308, which are at a physical distance of 569,215 bp in a region that contains 80 genes. Thus, the major effect QTL was fine-mapped to an approximately 0.60 Mb interval (from 1,248,093 to 1,817,308 bp) on Chr2 of the Cucumber v2 Genome assembly.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/478559/fpls-10-01425-HTML-r1/image_m/fpls-10-01425-g004.jpg)
Figure 4 Fine mapping of the major effect QTL fw2.1 on chromosome 2. The fw2.1 QTL was located between markers InDel1248093 and InDel1817308, shown in the black box. White bars indicate Fusarium oxysporum f. sp. cucumerinum (Foc)–sensitive Superina genotype; black bars indicate Foc-resistant Rijiecheng genotype; grey bars indicate F1 genotypes. DG is the disease grade. (A–F) indicate the six different genotypes that were detected among the F2 plants.
Analysis of Candidate Genes by RNA-Seq
To confirm the FW resistance genes, we analyzed the 80 candidate genes in the major effect QTL using RNA-sequencing (RNA-seq). We combined the RNA-seq data with the previously obtained transcriptome data [national center for biotechnology information’s (NCBI’s) sequence readarchive (SRA) database: PRJNA472169] to detect differentially expressed genes (DEGs) in the 0.60 Mb interval on Chr2. All the DEGs were positively regulated, except for Csa2G008760. DEGs were defined as having a false discovery rate ≤0.01 and fold change ≥2 or ≤−2. Twelve genes were found to be differentially induced by FW (Figure 5), namely, Csa2G007990 (calmodulin), Csa2G008030 [probable Guanosine diphosphate (GDP)-mannose transporter 2], Csa2G008110 (monogalactosyldiacylglycerol synthase), Csa2G008760 (chitinase), Csa2G008770 (adenylate kinase), Csa2G008780 and Csa2G009330 (unknown proteins), Csa2G009300 (DNA replication licensing factor MCM5), Csa2G009360 (RING-type finger protein 126), Csa2G009430 (transmembrane protein), Csa2G009440 (serine-rich protein), and Csa2G009470 (betaine aldehyde dehydrogenase). We also analyzed the expression patterns of the 12 candidate genes in Rijiecheng by qPCR (Figure 6; Supplementary Figure 1). The primers of these 12 genes are shown in Supplementary Table 1. We found that 5 of the 12 genes were significantly induced by FW, and their relative expression levels were higher in Foc-inoculated Rijiecheng than in the water-inoculated plants.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/478559/fpls-10-01425-HTML-r1/image_m/fpls-10-01425-g005.jpg)
Figure 5 Transcriptional profiles of 12 candidate genes in the major effect QTL fw2.1 on chromosome 2 at different times after inoculation. The expression data were converted to Log2 fragments per kilobase million (FPKM) to measure the expression levels of the FW resistance genes. The color scale indicates fold changes in the expression levels of the genes induced by FW.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/478559/fpls-10-01425-HTML-r1/image_m/fpls-10-01425-g006.jpg)
Figure 6 Relative expression patterns of five candidate genes in Superina and Rijiecheng plants inoculated with Foc and water control check (CK) . Each bar represents the average expression level of three independent biological replicates. Error bars show standard errors of the average values. * P ≤ 0.01–0.05 and ** P < 0.01 relative to the expression prior to inoculation with water.
Discussion
Several studies on the inheritance of resistance in cucumber investigated quantitative resistance traits controlled by multiple genes (Pierce and Wehner, 1990; Liu et al., 2003; Wang et al., 2005) and qualitative resistance traits controlled by a single gene (Netzer et al., 1977; Vakalounakis 2015). Contradictory results were obtained, possibly because there are different physiological races of Foc with different pathogenicities in different countries (Weng et al., 1989; Armsrong et al., 1987). We investigated the inheritance of cucumber FW resistance using Rijiecheng and Superina as the parents and concluded that the resistance in Rijiecheng was quantitative and the inheritance of FW resistance was controlled by multiple genes. This result was confirmed for F2 plants grown in different years and seasons, indicating that resistance to FW was a quantitative trait in cucumber. In dry pea, an F2 population segregating for high levels of resistance to Fusarium solani f. sp. Pisi also had a disease reaction phenotype in repeated greenhouse trials (Coyne et al., 2019).
Many studies on molecular linkage markers and genetic mapping of cucumber powdery mildew and downy mildew have been reported (Zhang et al., 2011; Zhang et al., 2013a; Zhang et al., 2013b; Xu et al., 2016). However, the genetic mapping of cucumber FW has rarely been reported (Zhang et al., 2014). One simple sequence repeat (SSR) marker linked to cucumber Foc2.1 was identified in a genetic interval of 5.98 cM (Wang, 2005) and validated among 46 germplasms. The accuracy rate of this SSR marker for selecting resistant germplasm was 87.88%. Zhang et al. (2014) found one major QTL on cucumber Chr2 (SSR03084–SSR17631) for FW resistance in a genetic distance of 2.4 cM. Zhou et al. (2015) concluded that cucumber Foc4 resistance to FW was located between SSR17631 and SSR00684 on Chr2. In this study, we found five InDel markers with polymorphisms in P1, P2, and F1 at different physical positions on Chr2 by BSA and successfully identified the major effect QTL on Chr2 with a physical distance of 0.60 Mb (InDel1248093–InDel1817308). Notably, the physical positions of all these QTLs are different: 2,526,888–3,262,528 (Zhang et al., 2014), 3,255,681–3,262,528 (Zhou et al., 2015), and 1,248,093–1,817,308 in the present study.
We found 12 DEGs in fw2.1 by combining BSA with RNA-seq after inoculating the plants with Foc infection and taking samples at different times. In previous studies, RNA-seq data were used to explore key susceptible genes in cucumber that responded to foliage diseases (Zheng et al., 2019). Liu et al. (2019) suggested that the differences in gene expression following cucumber mosaic virus (CMW) infection might explain the different resistance levels of two lines on the basis of RNA-seq data. Shi et al. (2018) analyzed the RNA-seq data of a recombinant inbred line and located a QTL on Chr5 that contained nine genes in cucumber infected with CMW. We found five candidate genes by combining QTL mapping and RNA-seq for resistance to FW in cucumber. Csa2G007990 encodes calmodulin and is highly expressed in cucumber. Naveed et al. (2019) found that when calmodulin bound to the effector peptide protein Avrblb2, it incited a hypersensitive response. Csa2G009430 encodes a transmembrane protein that has been widely studied in animal disease resistance, but not in plants. Csa2G009440 encodes a serine-rich protein, which has been linked to the ability of bacteria to attach to the hosts. For example, Santiago et al. (2007) suggested that serine-rich region may be involved in the negative regulation of phytochrome signal transduction, which is known to yield a hyperactive photoreceptor. Csa2G008780 and Csa2G009330 are highly expressed in cucumber, and their expression levels were significantly higher in Foc-inoculated plants than that in water-inoculated plants. The functions of these two genes are still unknown and require further study to determine their characteristics and functions. We plan to verify the functions of the five candidate genes using other methods such as virus-induced gene silencing and overexpression. Our results will help to better understand the genetic mechanisms and provide a strong basis for fine mapping of the major effect QTL and for cloning the candidate genes for resistance to FW in cucumber.
Data Availability Statement
The transcriptome data associated with this study can be found in NCBI using accession number PRJNA472169 (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA472169).
Author Contributions
XC and XX conceived the experiment. JD and JX performed the research and wrote the manuscript. JD and QX analyzed the data. All authors reviewed and approved this submission.
Funding
This research was supported financially by the Jiangsu Agriculture Science and Technology Innovation Fund (CX (17) 2004) and SCX (19) 3029, Special Funds for Three New Agricultural Projects in Jiangsu Province (SXGC [2017] 303), the National Natural Science Foundation of China (31902015), and Project of breeding for the major new agricultural variety of Jiangsu Province (PZCZ201720).
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.
Acknowledgments
We also thank Margaret Biswas, PhD, from Liwen Bianji, Edanz Group China (http://www.liwenbianji.cn/ac), for editing the English text of a draft of this manuscript.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2019.01425/full#supplementary-material
SUPPLEMENTARY FIGURE 1 | Relative expression patterns of the seven rejected candidate genes in the Superina and Rijiecheng plants inoculated with Foc and water (CK). Each bar represents the average expression level of three independent biological replicates. Error bars show standard errors of the average values. *P ≤0.01–0.05 and **P <0.01 relative to the expression prior to inoculation by water.
References
Armsrong, G. M., Armstrong, J. K., Neter, U. (1987). Pathogenic races of the cucumber-wilt Fusarium. Plant Dis. Rep. 62, 824–828.
Cao, Q. H., Wan, J., Chen, J. F. (2007). Research advances on downy mildew resistance in cucumber. China Cucurbits Veg 1, 27–30. doi: 10.16861/j.cnki.zggc.2007.01.013
Chen, J. B., Jing, R. L., Yuan, H. Y., Wei, B., Chang, X. P. (2005). Single nucleotide polymorphism of TaDREB1 gene in Wheat germplasm. Sci Agric Sin. 38, 2387–2394. doi: 10.1109/CIMCA.2005.1631314
Coyne, C. J., Porter, L. D., Boutet, G., Ma, Y., Mcgee, R. J., Lesne, A., et al. (2019). Confirmation of Fusarium root rot resistance QTL Fsp-Ps 2.1 of pea under controlled conditions. BMC Plant Biol. 98, 1–8. doi: 10.1186/s12870-019-1699-9
Falconer, D. S., Mackay, T. F. C. (1996). Introduction to Quantitative Genetics, 4th ed. Am. J. Hum. Genet. 46, 1231. doi: 10.1016/0168-9525(96)81458-2
Feng, F. J., Luo, L. J., Li, Y., Zhou, L. G., Xu, X. Y., Wu, X. Y., et al. (2005). Comparative analysis of polymorphism of InDel and SSR markers in Rice. Mol. Plant Breed. 3, 725–730. doi: 10.1080/00207177508921984
Gai, J. Y., Wang, J. K. (1998). Identification and estimation of a QTL model and its effects. Theor. Appl. Genet. 97, 1162–1168. doi: 10.1007/s001220051005
Giovannoni, J. J., Wing, R. A., Ganal, M. W., Tanksley, S. D. (1991). Isolation of molecular markers from specific chromosome intervals using DNA pools from existing mapping populations. Nucleic Acids Res. 19, 6553–6568. doi: 10.1093/nar/19.23.6553
Hill, J. T., Demarest, B. L., Bisgrove, B. W., Gorsi, B., Su, Y. C., Yost, H. J. (2013). MMAPPR: mutation mapping analysis pipeline for pooled RNA-seq. Genome Res. 23, 687–697. doi: 10.1101/gr.146936.112
Hiroki, T., Akira, A., Kentaro, Y., Shunichi, K., Satoshi, N., Chikako, M., et al. (2013). QTL-seq: rapid mapping of quantitative trait loci in rice by whole genome resequencing of DNA from two bulked populations. Plant J. 74, 174–183. doi: 10.1111/tpj.12105
Li, G., Carver, B. F., Cowqer, C., Bai, G., Xu, X. (2018). Pm223899, a new recessive powdery mildew resistance gene identified in Afghanistan landrace PI 223899. Theor. Appl. Genet. 131 (12), 2775–2783. doi: 10.1007/s00122-018-3199-y
Li, Z. Z., Xu, Q. K., Yu, H. P., Zhou, T. T., Xue, D. W., Zeng, D. L., et al (2017). Genetic analysis and gene mapping of yellow leaf and dwarf (yld) mutant in rice. Acta Agronomica Sinica. 4, 522–529. doi: 10.3724/SP.J.1006.2017.00522
Liu, D. L., Yang, R. H., Ha, Y. H. (2003). Research on genetic characteristics of cucumbers resisting to Fusarium wilt. Tianjin Agric. Sci. 9, 33–35. doi: 10.3969/j.issn.1006-6500.2003.02.011.
Liu, D., Cheng, Y. Z., Gong, M., Zhao, Q., Jiang, C. H., Cheng, L. R., et al. (2019). Comparative transcriptome analysis reveals differential gene expression in resistant and susceptible tobacco cultivars in response to infection by cucumber mosaic virus. Crop J. 7, 307–321. doi: 10.1016/j.cj.2018.11.008
Mansur, L. M., Orf, J., Lark, K. G. (1993). Determining the linkage of quantitative trait loci to RFLP markers using extreme phenotypes of recombinant inbreds of soybean (Glycine max L. Merr.). Theor. Appl. Genet. 86, 914–918. doi: 10.1007/BF00211041
Martínez, R., Aguilar, M. I., Guirado, M. L., Álvarez, A., Gómez, J. (2003). First report of fusarium wilt of cucumber caused by Fusarium oxysporum in Spain. Plant Pathol. 52, 410. doi: 10.1046/j.1365-3059.2003.00832.x
Michelmore, R. W., Paran, I., Kesseli, R. 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.A. 88, 9828–9832. doi: 10.1073/pnas.88.21.9828
Naveed, Z. A., Bibi, S., Ali, G. S. (2019). The Phytophthora RXLR effector Avrblb2 modulates plant immunity by interfering with Ca2+ signaling pathway. Front Plant Sci. 10. doi: 10.3389/fpls.2019.00374
Netzer, D., Niegro, S., Galun, E. (1977). A dominant gene conferring resistance to Fusarium wilt in cucumber. Phytopathology 77, 525–527. doi: 10.1094/Phyto-67-525
Owen, J.,.H. (1955). Fusarium wilt of cucumber. Phytopathology 45, 435–443. doi: 10.1111/j.1399-3054.1955.tb07793.x
Pierce, L. K., Wehner, T. C. (1990). Review of genes and linkage groups in cucumber. Hort Sci. 25, 605–615. doi: 10.21273/HORTSCI.25.6.605
Pu, Z. J., Zhang, Y. J., Liu, D., Dai, L. T., Wang, W. B. (2011). Research progress in biological control strategies for Fusarium wilt of cucumber. China Vegetables 6, 9–14.
Santiago, A. T., Dimitry, D., Andreas, H., Christian, F., Jorge, J. C. (2007). The serine-rich N-terminal region of Arabidopsis phytochrome A is required for protein stability. Plant Mol. Biol. 63, 669–678. doi: 10.1007/s11103-006-9115-x
Shang, S. J. D. (2009). Development and application of RGA markers in wheat. Chin. Acad. Agric. Sci.
Shi, L. X., Yang, Y. H., Xie, Q., Miao, H., Bo, K. L., Song, Z. C., et al. (2018). Inheritance and QTL mapping of cucumber mosaic virus resistance in cucumber (Cucumis Sativus L.). PloS One 13, 1–12. doi: 10.1371/journal.pone.0200571
Su, X., Zhu, G., Huang, Z., Wang, X., Guo, Y., Li, B., et al. (2019). Fine mapping and molecular marker development of the Sm gene conferring resistance to gray leaf spot (Stemphylium spp.) in tomato. Theor. Appl. Genet. 132 (4), 871–882. doi: 10.1007/s00122-018-3242-z
Takagi, H., Abe, A., Yoshida, K., Kosugi, S., Natsume, S., Mitsuoka, C., et al. (2013). QTL-seq: rapid mapping of quantitative trait loci in rice by whole genome resequencing of DNA from two bulked populations. Plant J. 74, 174–183. doi: 10.1111/tpj.12105
Vakalounakis, D. J. (2015). Inheritance and linkage of resistance in cucumber line SMR-18 to races 1 and 2 of Fusarium oxysporum f. sp. cucumerium. Plant Pathol. 44, 169–172.
Vakalounakis, D. J. (1993). Inheritance and genetic linkage of fusarium wilt (Fusarium oxysporum f. sp. cucumerinum race 1) and scab (Cladosporium cucumerinum) resistance genes in cucumber (Cucumis sativus L.). Ann. Appl. Biol. 123, 359–365. doi: 10.1111/j.1744-7348.1993.tb04098.x
Wang, Y. J. (2005). Studies on molecular marker of fusarium wilt resistance related gene in cucumber (Cucumis sativus L.). Northwest A&F University. 1–55.
Weng, Z. X., Xu, X. B., Feng, D. X. (1989). Study on the physiological race of Fusarium oxysporum f. sp. cucumerinum Owen. Chin. Vegetables 1, 19–21.
Xu, X. W., Ji, J., Xu, Q., Qi, X., Chen, X. (2017). Inheritance and quantitative trail loci mapping of adventitious root numbers in cucumber seedlings under waterlogging conditions. Mol. Genet. Genomics 292, 353–364. doi: 10.1007/s00438-016-1280-2
Xu, X., Yu, T., Xu, R., Shi, Y., Lin, X. Y., Xu, Q., et al. (2016). Fine mapping of a dominantly inherited powdery mildew resistance major-effect QTL, Pm1.1, in cucumber identifies a 41.1kb region containing two tandemly arrayed cysteine-rich receptor-like protein kinase genes. Theor. Appl. Genet. 129 (3), 507–516. doi: 10.1007/s00122-015-2644-4
Ye, S. F., Yu, J. Q., Peng, Y. H., Zheng, J. H., Zou, L. Y. (2004). Incidence of Fusarium wilt in Cucumis sativus L. is promoted by cinnamic acid, an autotoxin in root exudates. Plant Soil 263, 143–150. doi: 10.1023/B:PLSO.0000047721.78555.dc
Yuan, F., Zhang, C. L., Shen, Q. R. (2003). Alleviating effect of phenol compounds on cucumber fusarium wilt and mechanism. Agric. Sci. China 2, 647–652.
Zhang, K., Wang, X., Zhu, W., Qin, X., Xu, J., Cheng, C., et al. (2018). Complete resistance to powdery mildew and partial resistance to downy mildew in a Cucumis hystrix introgression line of cucumber were controlled by a co-localized locus. Theor. Appl. Genet. 131 (10), 2229–2243. doi: 10.1007/s00122-018-3150-2
Zhang, S. P., Liu, M. M., Miao, H., Zhang, S. Q., Yang, Y. H., Xie, B. Y., et al. (2011). QTL mapping of resistance genes to powdery mildew in cucumber (Cucumis sativus L.). Sci Agric Sin. 129 (3), 507–516. doi: 10.3864/j.issn.0578-1752.2011.17.011
Zhang, S. P., Liu, M. M., Miao, H., Zhang, S. Q., Yang, Y. H., Xie, B. Y., et al. (2013a). Chromosomal mapping and QTL analysis of resistance to downy mildew in Cucumis sativus. Plant Dis. 97, 245–251. doi: 10.1094/PDIS-11-11-0941-RE
Zhang, S. P., Miao, H., Yang, Y. H., Xie, B. Y., Wang, Y., Gu, X. F. (2014). A major quantitative trait locus conferring resistance to fusarium wilt was detected in cucumber by using recombinant inbred lines. Mol. Breed. 34, 1805–1815. doi: 10.1007/s11032-014-0140-1
Zhang, S., Miao, H., Sun, R., Wang, X., Huang, S. W., Wehner, T. C., et al. (2013b). Localization of a new gene for bitterness in cucumber. J. heredity 104, 134–139. doi: 10.1093/jhered/ess075
Zheng, X. Y., Yang, J. J., Lou, T. X., Zhang, J., Yu, W. J., Wen, C. L. (2019). Transcriptome profile analysis reveals that CsTCP14 induces susceptibility to Foliage Disease in cucumber. Int. J. Mol. Sci. 20, 2582. doi: 10.3390/ijms20102582
Zhou, H. M., Mao, A. J., Zhang, L. R., Zhang, F., Wang, Y. J., Yang, W. C. (2010). Research on inoculation method and the inheritance of resistance to fusarium oxysporum f. sp. cucumerinum on cucumber. Acta Agric. Boreali-Sinica 25, 186–190. doi: 10.3724/SP.J.1238.2010.00569
Zhou, X. G., Wu, F. Z., Wang, X. Z., Yuan, Y. (2008). Progresses in the mechanism of resistance to fusarium wilt in cucumber. J. Northeast Agric. Univ. 15, 1–6.
Zhou, H. M., Dong, C. J., Zhang, H. Y., Ren, Y., Guo, S. G., Yang, W. C., et al. (2015). SSR molecular markers and location of the gene Foc-4 linked to the resistance to Fusarium Wilt of cucumber. Mol. Plant Breed. 13, 1980–1986. doi: 10.13417/j.gab.034.001980
Keywords: cucumber, Fusarium wilt, inheritance, quantitative trait locus fine mapping, resistance genes
Citation: Dong J, Xu J, Xu X, Xu Q and Chen X (2019) Inheritance and Quantitative Trait Locus Mapping of Fusarium Wilt Resistance in Cucumber. Front. Plant Sci. 10:1425. doi: 10.3389/fpls.2019.01425
Received: 14 June 2019; Accepted: 14 October 2019;
Published: 02 December 2019.
Edited by:
Feishi Luan, Northeast Agricultural University, ChinaReviewed by:
Luming Yang, Henan Agricultural University, ChinaQunfeng Lou, Nanjing Agricultural University, China
Meiling Gao, Qiqihar University, China
Copyright © 2019 Dong, Xu, Xu, Xu and Chen. 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: Xuehao Chen, xhchen@yzu.edu.cn