Click on author name to view affiliation information
Head splitting resistance (HSR) in cabbage is an important trait closely related to appearance, yield, storability, and mechanical harvestability. In this study, a doubled haploid (DH) population derived from a cross between head splitting-susceptible inbred cabbage line 79-156 and resistant line 96-100 was used to analyze inheritance and detect quantitative trait loci (QTLs) for HSR during 2011–12 in Beijing, China. The analysis was performed using a mixed major gene/polygene inheritance method and QTL mapping. This approach, which uncovered no cytoplasmic effect, indicated that HSR can be attributed to additive-epistatic effects of three major gene pairs combined with those of polygenes. Major gene and polygene heritabilities were estimated to be 88.03% to 88.22% and 5.65% to 7.60%, respectively. Using the DH population, a genetic map was constructed with simple sequence repeat (SSR) markers anchored on nine linkage groups spanning 906.62 cM. Eight QTLs for HSR were located on chromosomes C4, C5, C7, and C9 based on 2 years of phenotypic data using both multiple-QTL mapping and inclusive composite interval mapping. The identified QTLs collectively explained 37.6% to 46.7% of phenotypic variation. Three or four major QTLs (Hsr 4.2, 7.2, 9.3, and/or 9.1) showing a relatively larger effect were robustly detected in different years or with different mapping methods. The HSR trait was shown to have a complex genetic basis. Results from QTL mapping and classical genetic analysis were consistent. Our results provide a foundation for further research on HSR genetic regulation and molecular marker-assisted selection (MAS) for HSR in cabbage.
Cabbage (Brassica oleracea L. var. capitata) is one of the world’s most widely cultivated vegetables. In addition to being a significant year-round component of the global vegetable supply, cabbage is of immense importance for human nutrition, providing dietary fiber, vitamins, and cancer-preventing substances (Fahey and Talalay, 1995; Masarirambi et al., 2011). During later stages of vegetative growth, cabbage heads are vulnerable to cracking. This phenomenon seriously affects appearance, yield, storability, and mechanical harvestability. In addition, susceptibility to head splitting hinders the prolongation of harvest time and thus the ability of growers to select harvest times for optimal selling price. Head splitting resistance is thus a very desirable property in cabbage (Holt and Schoorl, 1983). Previous studies (Liu et al., 2009; Zeng et al., 2009) about this trait mainly concentrated on cultivating conditions, physiological, cytological characters, and types of the head. It is proved that many factors, such as irrigation, fertilizer, planting density, endogenous hormones content, leaf surface microconfiguration, and tissue character could affect the cracking head. To alleviate late-stage head splitting, farmers generally reduce irrigation times during cultivation. Such practices affect normal cabbage head growth, however, and cannot completely prevent head splitting. Under water-sufficient conditions, head splitting in cabbage is mainly caused by genetic factors (Qin et al., 1994). Improvement of HSR in newly developed varieties has therefore become a priority in cabbage breeding programs. The complex genetic nature of the HSR trait has only been explored in a few studies. After performing a cross between early- and late-splitting cabbage lines and obtaining their F1, F2, and backcross progenies, Chiang (1972) concluded that HSR is controlled by at least three pairs of genes within an average of six generations. Gene action was found to be mostly additive. Partial dominance for early splitting was detected, and narrow sense heritability was estimated as 47%. Zhuang et al. (2009) further analyzed combining ability and heritability according to the Griffing IV method. Using 15 cross combinations involving six backbone parents, they revealed that HSR is controlled by both additive and nonadditive effects. Additive effects were found to be the most important, with the percentage contribution of additive effects increasing with prolongation of field duration. These classical genetic methods revealed that HSR is complex and controlled by many genes, and also estimated the total gene effect. Nevertheless, major gene effects and polygene effects were not clearly distinguished and the gene relationships were not identified.
Major gene/polygene genetic segregation analysis and QTL mapping are the main approaches used to clarify the genetic basis of quantitative traits (Gai et al., 2007; Zhang et al., 2003). These methods have been successfully applied to uncover inheritance patterns and QTLs of important traits in many crops, vegetables, and flowers (Anbessa et al., 2006; Cheng et al., 2011; Hao et al., 2008; Liang et al., 2014; Wang et al., 2012b; Zhang et al., 2010, 2013).
For QTL mapping, the use of DH populations is a powerful tool. Because each line is homozygous, they can be replicated between test sites, and trialed over multiple years. So the standard error of QTL genotype means is thus decreased, thereby allowing a better estimate of trait heritability and increasing QTL detection ability (Pink et al., 2008; Simon et al., 2008). At the same time, multiple mapping methods for detecting QTLs are needed to effectively identify and verify important QTLs (Su et al., 2010; Xing et al., 2012).
Recent progress in genetic mapping and molecular marker development in B. oleracea have laid an important foundation for genetic studies and QTL analysis of important B. oleracea traits (Brown et al., 2014; Gao et al., 2007; Walley et al., 2012; Wang et al., 2012a). QTL mapping studies in cabbage have mainly concentrated on plant size (Lan and Paterson, 2001), flowering time (Okazaki et al., 2007), and disease resistance (Kifuji et al., 2012; Nagaoka et al., 2010; Mei et al., 2013). No QTL mapping studies associated with cabbage HSR have been reported.
To thoroughly dissect the genetic architecture of HSR in cabbage, we first developed an intracrop ‘immortal’ DH mapping population derived from a cabbage x cabbage cross and used it in a segregation analysis over multiple generations (P1, P2, and DHs) to explain HSR inheritance. We then constructed a framework linkage map based on this DH population with SSR markers. Finally, we located the first-known cabbage QTLs for HSR on the generated linkage map and analyzed QTL stability across years with different mapping programs.
Two inbred cabbage lines (Fig. 1), 96-100 and 79-156, were used as parents to develop a DH population. The male parental line, 96-100 (P2, about 50 d form transplantation to maturity in autumn), is resistant to head splitting and was self-developed from germplasm of Indian origin introduced by the Bejo Sheetal Company; the female parental line, 79-156 (P1, about 50 d form transplantation to maturity in autumn), is susceptible to head splitting and self-developed from a germplasm introduced from Denmark. Both of the parents are early maturing inbred lines and were provided by the Cabbage and Broccoli Laboratory, Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences (CAAS). The DH population was derived from the F1 (79-156 × 96-100) by microspore culture (Yuan et al., 2011, 2012) during 2009–12. Seeds were obtained from 157 DH individuals in the summer of 2011 and from 181 individuals in 2012. All DH lines obtained in the 2 years were used for genetic segregation analysis. For linkage map construction and mapping, data from the first 157 lines were used.
Citation: HortScience horts 50, 7; 10.21273/HORTSCI.50.7.944
Parents, F1 (79-156 × 96-100), RF1 (96-100 × 79-156), and DH lines were planted in the field at the experimental farm of the Institute of Vegetables and Flowers, CAAS (Changping, Beijing, China) in 2011 and 2012. Seeds of each generation were sown in a greenhouse in mid-July, and the resulting seedlings were transplanted into 10-cm diameter plastic pots after 3 weeks. The plants were then transplanted into the field at ≈5 weeks of age (19 Aug. 2011 and 22 Aug. 2012). Parental, F1, and RF1 lines were planted in a randomized complete block design with three replicates. In each replicate, parents and F1 lines were planted in three rows; each row comprised 16 plants, spaced 40 cm apart within a row and 50 cm between rows. For the DH population, a block in replication design was adopted with three replicates (Gai, 2000; Supplemental Table 1). Each replicate consisted of single rows of 16 plants. The experimental plots were surrounded by two additional rows planted to serve as a protective buffer. Irrigation and insecticide application were consistent with commercial production conditions. In each replicate, 15 plants at the same growth level of each parental, F1, and RF1 line and 10 plants of each DH line were marked for phenotypic measurements.
Because of differences in maturation rates due to genotypic variation, maturation dates were recorded for each DH line. At maturity, the height and the circumference of each marked cabbage head at its widest point were measured with a flexible plastic ruler. The arc length and width of the largest split were also measured with a flexible plastic ruler 15 d after maturity. Splitting was assessed on a 6-point scale, assigned according to the number of split layers and the ratio of the split size to the entire surface area (calculated as S [%] = S1/S2 × 100%, where S1 = arc length × width of the largest split and S2 = height × half of the largest circumference of the head) (Fig. 2). Assigned splitting scores were as follows: 0 = no split; 1 = 1 split layer; 2 = 2 split layers and S < 50; 3 = 3–5 split layers and S < 50, or 2 split layers and S ≥ 50; 4 = 6–10 split layers and S < 50, or 3–5 split layers and S ≥ 50; 5 = more than 10 split layers, or 6–10 split layers and S ≥ 50. The head splitting index was calculated as: ∑ (splitting score × number of plants with that score)/(total number of plants × highest possible splitting score) × 100%. The data were analyzed with Excel 2003 and SAS 8.1 software.
Citation: HortScience horts 50, 7; 10.21273/HORTSCI.50.7.944
To determine an appropriate genetic model for the HSR trait, the cabbage DH population dataset was analyzed by major gene/polygene mixed inheritance analysis (Gai, 2006; Gai et al., 2003). Joint segregation analysis was performed following Hao et al. (2008) and Zhang et al. (2010) under the basic assumptions described by those authors. To select the genetic model best explaining the quantitative trait variation, 38 genetic models of 7 different types were considered (Supplemental Table 2). Maximum likelihood estimates of component parameters in each genetic model were generated using the iterated expectation and conditional maximization algorithm (Zhang et al., 2003). The proportion, mean, and variance of each component distribution in the likelihood function were included in the estimates. Best-fit models were chosen according to the Akaike information criterion (AIC) (Akaike, 1977) and a suite of goodness-of-fit tests. Finally, genetic parameters of major genes and polygenes were estimated based on the least squares principle.
Total DNA was isolated from expanding leaves of 3-week-old plants using the modified cetyltrimethylammonium bromide method (Saghai-Maroof et al., 1984). The genomic DNA samples were diluted to 100 ng/μl with Tris-EDTA (pH 8.0) and stored at −20 °C for use as polymerase chain reaction (PCR) templates. In addition, leaf tissue was lyophilized for use in future DNA extractions.
A set of 2,170 SSR markers developed from cabbage sequence scaffolds (Wang et al., 2012a) and 1,013 expressed sequence tag (EST)-SSR markers (Chen et al., 2010a) were used to scan for polymorphisms between the two parents. DNA amplification of SSR markers was carried out in 20-μl volumes containing 1 unit of Taq polymerase, 0.1 μM of each primer, 200 μM dNTPs, 2 μl of 10 × buffer (25 mm Mg2+), and 100 ng genomic DNA template. The PCR thermal profile was as follows: initial denaturation for 5 min at 94 °C, followed by 35 cycles of DNA denaturation for 30 s at 94 °C, annealing for 30 s at 55 °C, and extension for 45 s at 72 °C, and a final extension for 7 min at 72 °C. PCR amplifications were carried out in a GeneAmp PCR system 9700 (Applied Biosystems, Foster City, CA). The resulting products were subjected to 8% polyacrylamide gel electrophoresis at 160 V for 1.5 h followed by silver staining (Brant et al., 1991).
For map construction, the DH population was genotyped for all SSR markers that showed polymorphisms between the parental 79-156 and 96-100 lines. The genotyping data were coded as type “a” or “b”, corresponding to parental lines 79-156 and 96-100, respectively, with ambiguous and missing data indicated by “-”. A linkage map was constructed using JoinMap 4.0 software (Van Ooijen, 2006). For map distance calculations, recombination frequencies were converted to centiMorgans (cM) using Kosambi’s method, and linkage groups were assigned to chromosomes C1–C9 of B. oleracea based on markers in common with the reference (Lv et al., 2014; Wang et al., 2012a).
QTLs were estimated by multiple-QTL modeling (MQM) with MapQTL 4.0 (Van Ooijen et al., 2002) and by inclusive composite interval mapping (ICIM) using QTL IciMapping v3.0 software (Li et al., 2007). In MQM, a 1,000-permutation run was first performed to estimate the significance threshold of the test statistic for a QTL based upon a 5% experiment-wise error rate. Interval mapping at 1-cM intervals along the chromosomes was then used to scan for QTLs based on a logarithm of odds (LOD) threshold of 2.5. Markers closely linked to positions with the highest LOD score were taken as cofactors for MQM analysis. To select significant markers during the first step of ICIM stepwise regression, P values for entering and removing variables were set respectively at 0.001 and 0.002; in the second step, a minimum LOD threshold of 2.5 was used to declare a QTL significant.
Parents 79-156 and 96-100 differed significantly with respect to HSR (Fig. 1; Table 1). Line 96-100 was resistant to head splitting, with a head splitting index of 9.33 and 9.36 in 2011 and 2012, respectively; in contrast, 79-156 was highly susceptible to head splitting, with corresponding head splitting index values of 86.87 and 85.32. Head splitting indexes of F1 and RF1 were not significantly different from one another; the values were intermediate to those of the parents, indicating that there is no cytoplasmic effect on HSR inheritance.
Analysis of variance revealed significant differences (P < 0.01) in HSR among DHs in both 2011 and 2012 (Table 2), indicating the existence of heritable variation and is suitable for genetic analysis. Significant differences between years and among lines were observed at P < 0.05 and P < 0.01 levels, respectively; this result demonstrates that the HSR trait is mainly under genetic control, with climate also playing an important role.
In the DH population, the head splitting index showed continuous variation, suggesting that comprehensive HSR in cabbage is a typical quantitative genetic character (Fig. 3). Nevertheless, the frequency distribution of phenotypes deviated from a normal distribution, with skewness and kurtosis values of 0.76 and −0.73 in 2011 and 0.97 and −0.29 in 2012, respectively. These skewed and multipeak phenomena indicate the possible existence of major genes for HSR. The head splitting index of DHs ranged from 0.00 to 98.52 in 2011 and 0.00 to 100.00 in 2012, with maximum values of the index much greater than that in the parents. The observed transgressive segregation indicates that genes controlling HSR are scattered throughout the genome. Extreme phenotypes can produce in both positive and negative directions through gene recombination. In 2011 and 2012, the corresponding average head splitting index values were 35.16 and 30.46, with the distribution of the data in the DH population skewed toward the resistant parental type.
Citation: HortScience horts 50, 7; 10.21273/HORTSCI.50.7.944
Models F-1, F-2, G-0, and G-1 in 2011 and B-1-2, F-1, G-0, and G-1 in 2012 were chosen as candidates according to the smaller Akaike’s Information Criterion (AIC) (Table 3). Fitness evaluation including the uniform test (U12, U22, U32), Smirnov test (nW2) and Kolmogorov test (Dn) ruled out the applicability of F-2, F-1, and G-1 models because 5, 4, and 4 of the 15 parameters were observed to be significant in the 3 respective tests, while no statistic was significant under the G-0 model in 2011. The B-1-2 model could be eliminated in 2012 because five of the 15 parameters were found to be significant. Although only one parameter was significant under the other three models, the G-0 model was regarded as the best fit because it was associated with the lowest AIC score (Supplemental Table 3). We therefore deduced that HSR in cabbage can be described by a model corresponding to three pairs of additive-epistatic major genes plus additive-epistatic polygenes.
Maximum likelihood frequency values of eight component distributions were estimated under the G-0 model, with first- and second-order genetic parameters then calculated from the results (Table 4) using the least squares method. The calculated parameters showed similar tendencies in both years; additive effects (d) of the three major genes were estimated as 20.4, 8.59, and 15.89. All the additive effects were positive, indicating that the susceptible parental line 79-156 had important effects on total variability. The additive epistatic effect of the first two major genes, iab, was 3.64, while that between the first and third major genes, iac, was 10.94. The additive epistatic effect of the second and third major genes, ibc, was only −0.86, while that among the three major genes, iabc, was −5.81. Similar results were observed in 2012. Heritabilities of major genes and polygenes in 2011 were 88.03 and 7.60%, respectively, with corresponding values of 88.22 and 5.65% in 2012. The heritability values of the major genes were much larger than those of the polygenes. These observations demonstrate that HSR is mainly controlled by major genes, and that selection for this trait should be carried out in early generations.
A preliminary screening of 79-156 and 96-100 parental genotypes using SSR markers developed from cabbage sequence scaffolds and B. oleracea ESTs identified 252 polymorphic markers between the two parents. Of these polymorphic markers, 149 were used to genotype 157 members of the DH population. In addition to 12 EST-SSR markers, 123 scaffold-SSR markers were incorporated into the DH linkage map. The linkage analysis uncovered nine linkage groups. Based on the presence of reference SSR markers, the linkage groups were designated as C1–C9 in accordance with the nomenclature used by Wang et al. (2012a) and Lv et al. (2014).
The framework linkage map (Fig. 4) was 906.62 cM in length, with an average between-marker distance of 6.72 cM, a minimum between-marker distance of 0.46 cM, and a maximum distance of 34.90 cM. The largest linkage group, encompassing 189.46 cM, was C3; the smallest, C2, spanned 32.01 cM (Table 5). The maximum average distance between markers (10.67 cM) was that of C2, which featured the lowest number of markers (3); the minimum average distance (3.34 cM) was found on C7, which had the highest number of markers (24). The cabbage genome has been variously estimated to comprise 603 (Uzunova et al., 1995) or 630 Mbp (Liu et al., 2014). Using the latest estimate of Liu et al. (2014), the average physical distance between mapped markers was calculated to be 4.67 Mbp.
Citation: HortScience horts 50, 7; 10.21273/HORTSCI.50.7.944
Segregation distortion is a feature common to Brassica DH populations, indicating possible preferential selection of genotypes responsive to microspore culture and/or the ability to produce seed during regeneration and seed-bulking phases (Sebastian et al., 2000; Wang et al., 2012a). Based on a χ2 test for goodness of fit to the expected 1:1 Mendelian segregation ratio, 90 of 135 loci displayed varying degrees of distortion (P ≤ 0.01) (Table 5). Although slightly more, 79-156 alleles (63.3%), were present among the 90 distorted loci compared with 96-100 (36.7%) alleles (Table 5), there was not a significant departure from a 1:1 Mendelian ratio. Linkage groups C1, C4, and C9 had clusters of markers that were distorted toward 96-100 alleles, whereas linkage groups C3, C5, C6, and C7 contained clusters of markers distorted in favor of 79-156. Linkage group C8 included small clusters of both parental genotypes. No heterozygous loci were scored during genotyping of the molecular markers.
As shown in Table 6 and Fig. 4, eight QTLs for HSR during 2 years were detected using MQM and ICIM methods with Map QTL 4.0 and IciMapping 3.0 software. The QTLs were located on cabbage chromosomes C4, C6, C7, and C9 and individually explained 5.50% to 13.94% of observed phenotypic variation. Using the MQM program, six QTLs were detected on chromosomes C4, C7, and C9; these QTLs collectively accounted for 38.4% and 38.0% of the phenotypic variation in 2011 and 2012, respectively, with the effect of each QTL ranging from 7.1 to 12.1% over the 2 years. Hsr 4.2 and Hsr 7.2, which were robust QTLs showing a relatively large effect, could be detected in both years and were distributed between markers Scaffold12597a/Scaffold55516 and Scaffold195/Scaffold46873. In both years, the 96-100 allele at the Hsr 4.2 locus increased the tendency toward HSR, whereas this allele at the Hsr 7.2 locus decreased it. Using ICIM, five of the same QTLs were detected on chromosomes C4, C6, and C9 and together explained 37.6% and 46.7% of observed phenotypic variation in 2011 and 2012, respectively. Their individual effects ranged from 5.5% to 13.9%—the same trend observed using MQM. Three of the QTLs (Hsr 4.2, Hsr 9.2, and Hsr 9.3) either overlapped or were adjacent to the corresponding regions detected using MQM. Hsr 9.3 was located in the interval between BOE344–BOE975 on chromosome C9; it explained 10.1% to 13.9% of the phenotypic variation and showed the largest effect of the three QTLs. The locus Hsr 9.1, detected only by ICIM, was a major QTL explaining 9.2% to 10.5% of the phenotypic variation in both years. The 96-100 allele increased HSR at the Hsr 9.3 locus and decreased it at Hsr 9.1 in both years.
HSR is an important agronomic trait, being associated with appearance, yield, mechanical harvestability, marketability, and storability. Qin et al. (1994) have reported that head splitting in cabbage under water-sufficient conditions is mainly due to genetic factors. In the present study, head splitting is mainly controlled by genetic features, with environmental conditions having a lesser influence. Similar results have been recorded for chinese cabbage as well as watermelon, tomato, and various other fruit crops cracking (Budan, 1999; Cortés et al., 1983). The development of varieties having HSR is thus an essential strategy for breeders.
Few relevant genetic analyses have been conducted, however, and previous genetic studies (Chiang, 1972; Zhuang et al., 2009) of HSR inheritance have focused only on HSR as a polygenic system and ignored the effect of individual genes. The major gene/polygene genetic segregation analysis, used in the present study, has the ability of determining individual effects up to major genes and the collective epistatic effects of polygene (Wang et al., 2001). It is strongly recommended for plant breeders to apply it as a simple and useful technique to know the number of major genes, their kinds of genetic effects, heritability values and genetic information on all kinds of genetic effects, and heritability value of whole polygenes without any extra requirements on laboratory conditions except a precise experiment (Gai et al., 2007). It was used for head cracking in our previous study with six basic generations (Su et al., 2012). However, the present approach designed for the six basic populations was only capable to find the genetic mechanism up to two major genes plus polygenes (Gai et al., 2003) and the separation generations can’t be repeated, which may lead to large errors. In this study, we performed a segregation analysis of the HSR trait using a DH population comprising more than 150 lines with three replications over 2 years. HSR is primarily controlled by hereditable factors with the additive-epistatic effects of three major genes as well as polygenes, and no cytoplasmic effect was observed. Additive effects predominated overall other types of genetic effects, and, as in previous studies, higher heritabilities of the trait were recorded with respect to major genes than polygenes (Chiang, 1972; Su et al., 2012; Zhuang et al., 2009).
Compared with F2, DH populations are ideal for the genetic study of quantitative traits, because they are composed of genetically fixed DH lines, which can be replicated between test sites and trialed over years, decreasing the standard error of QTL genotype means and allowing a better estimate of trait heritability and increased power to detect QTLs (Soller and Beckmann, 1990). Although the number of available polymorphic loci may be reduced using intracrop crosses, such a strategy—as pointed out by Walley et al. (2012)—allows a direct relationship to be established between trait and crop type; the genetic variation captured in this fashion reduces the time required to incorporate important agronomic traits into elite breeding material. Here, we report a new cabbage linkage map using SSR markers in a DH population derived from an intrasubspecies cross of cabbage. Our newly constructed map and the generated DH lines are therefore not only important for research on related characteristics of cabbage, but will also contribute to the exchange of materials between laboratories and successive research in the future (Walley et al., 2012; Wang et al., 2012a).
With the advent of molecular markers, QTL mapping has become increasingly important in molecular breeding. Marker-assisted selection and gene discovery are now widely used for the breeding of field crops and vegetables (Cheng et al., 2011; Sabouri, 2009; Walley et al., 2012). Although some QTLs have been identified in cabbage, research is still at a preliminary stage: QTL cloning has not yet been reported and functional analysis studies are rare (Lv et al., 2014). For cabbage breeders, HSR is an important characteristic that affects both cabbage yield and quality, but no QTL studies have previously been reported for this trait.
The use of different genetic models, algorithms, and mapping procedures can produce different mapping results, even for the same set of data (Su et al., 2010). The results of different mapping methods can potentially complement one another. Whole-genome scanning with multiple mapping methods has thus been recommended for QTL mapping (Xing et al., 2012). Using data obtained over multiple years increases environmental heterogeneity and also allows improved estimates of QTLs that may not reach the genome-wide significance threshold in just one environment (Chen et al., 2010b; Piepho, 2005; Van Eeuwijk et al., 2010). We used two mapping methods to analyze agronomic trait data collected during 2 years in a DH population, and our study results, the first to identify QTLs associated with the HSR trait in cabbage, provide support for the above-mentioned assertions. For example, QTLs distributed on C7 and C6 were only identified by one mapping method; in addition, Hsr 9.2 and Hsr 9.3 could both be detected in 2011 and 2012 by ICIM, but only one was identified using MQM. On chromosome C9, QTL Hsr 9.1 with a larger effect is physically quite close to Hsr 9.2. Whether these QTLs are independent or instead both components of a larger QTL is an open question. In brief, three or four major QTLs (Hsr 4.2, Hsr 7.2, Hsr 9.3, and/or Hsr 9.1) and some minor QTLs were detected over 2 years using the two mapping methods. This result is consistent with the results of classical genetic analysis, both segregation analysis and QTL mapping with molecular markers can be used as a mutual control and supplement (Gai et al., 2007).
Because cabbage lines resistant to head splitting have traditionally been selected in the field during the mature stage, splitting severity depends on environmental factors and plant conditions. The use of DNA markers enables reliable selection of resistant plants, even at the seedling stage. The makers link to QTLs identified in our study will be helpful for the identification of genes related to HSR and for MAS in cabbage breeding programs.
Contributor Notes
This work was funded by the National High Technology Research and Development Program (863 Program, 2012AA100104, 2012AA100105), Chinese National Key Technology R&D Program (2013BAD01B04), the Earmarked Fund for Modern Agro-industry Technology Research System (CARS-25-A), and Key Laboratory of Biology and Genetic Improvement of Horticultural Crops, Ministry of Agriculture, People’s Republic of China.
To whom reprint requests should be addressed; e-mail liuyumei@caas.cn.
