Japanese foxtail is a grass weed infesting in oilseed rape and wheat (Triticum aestivum L.) fields in eastern China (Bi et al. Reference Bi, Liu, Guo, Li, Yuan, Du and Wang2016; Cui et al. Reference Cui, Wang, Han, Chen and Li2015; Tang et al. Reference Tang, Li, Dong, Dong, Lü and Zhu2012; Yang et al. Reference Yang, Dong, Li and Moss2007). Fenoxaprop-P-ethyl, an acetyl-CoA carboxylase (ACCase)-inhibiting herbicide, was firstly registered to control grass weeds in wheat by Bayer AG in 1998 in China (www.chinapesticide.gov.cn). Since then, fenoxaprop-P-ethyl has been used continuously by farmers to control grass weeds by farmers continuously. Many Japanese foxtail populations have evolved fenoxaprop-P-ethyl resistance, as well as cross-resistance to other ACCase-inhibiting herbicides (Cui et al. Reference Cui, Wang, Han, Chen and Li2015; Pan et al. Reference Pan, Li, Xia, Zhang and Dong2015; Tang et al. Reference Tang, Li, Dong, Dong, Lü and Zhu2012; Xu et al. Reference Xu, Zhu, Wang, Li and Dong2013; Yang et al. Reference Yang, Dong, Li and Moss2007).
Target-site resistance (TSR) and non–target site resistance are the two categories of mechanisms for resistance to ACCase-inhibiting herbicides in grasses (Kaundun Reference Kaundun2014). ACCase-inhibiting herbicides bind to the site of action in the carboxyl transferase (CT) domain of ACCase and inhibit de novo fatty acid synthesis in sensitive grass weeds, leading to necrosis and plant death (Kaundun Reference Kaundun2014; Post-Beittenmiller et al. Reference Post-Beittenmiller, Roughan and Ohlrogge1992). There are three types of ACCase-inhibiting herbicides: aryloxyphenoxypropionate (FOP), cyclohexanedione (DIM), and phenylpyrazolin (DEN). TSR prevents an herbicide from binding to the target enzyme because of amino acid substitutions at the CT domain of ACCase (Devine Reference Devine1997). Fourteen target-site mutations (amino acid substitutions) referring to six positions that induce resistance to different ACCase-inhibiting herbicides have been identified within the CT domain of homomeric chloroplastic ACCase. These substitutions include Ile-1781 to Leu, Val, Thr, or Ala (Délye et al. Reference Délye, Pernin and Michel2011; Kaundun Reference Kaundun2014; Zagnitko et al. Reference Zagnitko, Jelenska, Tevzadze, Haselkorn and Gornicki2001; Zhang and Powles Reference Zhang and Powles2006); Trp-1999 to Cys, Leu, or Ser (Liu et al. Reference Liu, Harrison, Chalupska, Gornicki, O’Donnell, Adkins, Haselkorn and Williams2007); Trp-2027 to Cys (Gherekhloo et al. Reference Gherekhloo, Osuna and De Prado2012; Petit et al. Reference Petit, Bay, Pernin and Délye2010; Yu et al. Reference Yu, Collavo, Zheng, Owen, Sattin and Powles2007); Ile-2041 to Asn or Val (Délye et al. Reference Délye, Zhang, Chalopin, Michel and Powles2003; Scarabel et al. Reference Scarabel, Panozzo, Varotto and Sattin2011; Tang et al. Reference Tang, Zhou, Chen and Zhou2014); Asn-2078 to Gly (Collavo et al. Reference Collavo, Panozzo, Lucchesi, Scarabel and Sattin2011; Kaundun Reference Kaundun2010); Cys-2088 to Arg (Yu et al. Reference Yu, Collavo, Zheng, Owen, Sattin and Powles2007); and Gly-2096 to Ala or Ser (Cruz-Hipolito et al. Reference Cruz-Hipolito, Domínguez-Valenzuela, Osuna and De Prado2012; Li et al. Reference Li, Du, Liu, Yuan and Wang2014).
TSR often causes different cross-resistance to herbicides that possess the same mode of action (Beckie and Tardif Reference Beckie and Tardif2012). Since 2007, when Japanese foxtail resistance to ACCase-inhibiting herbicides was first reported, multiple amino acid mutations have been identified (Bi et al. Reference Bi, Liu, Guo, Li, Yuan, Du and Wang2016; Cui et al. Reference Cui, Wang, Han, Chen and Li2015; Xu et al. Reference Xu, Zhu, Wang, Li and Dong2013, Reference Xu, Li, Zhang, Cheng, Jiang and Dong2014a; Yang et al. Reference Yang, Dong, Li and Moss2007). Previous studies showed that the Ile-1781-Leu substitution in Japanese foxtail confers high resistance to fenoxaprop-P-ethyl and cross-resistance to clodinafop-propargyl, clethodim, and pinoxaden (Bi et al. Reference Bi, Liu, Guo, Li, Yuan, Du and Wang2016; Cui et al. Reference Cui, Wang, Han, Chen and Li2015; Mohamed et al. Reference Mohamed, Li, You and Li2012). The Ile-2041-Asn substitution confers resistance to haloxyfop-R-methyl and cross-resistance to other FOPs in Japanese foxtail (Tang et al. Reference Tang, Li, Dong, Dong, Lü and Zhu2012). Additionally, the mutation at position 1999 confers resistance to fenoxaprop-P-ethyl, but not to clodinafop-propargyl, haloxyfop-R-methyl, or clethodim (Xu et al. Reference Xu, Li, Zhang, Cheng, Jiang and Dong2014a). The Trp-2027-Cys substitution in Japanese foxtail confers resistance to fenoxaprop-P-ethyl and cross-resistance to other FOPs and pinoxaden but not to DIMs (Xu et al. Reference Xu, Zhu, Wang, Li and Dong2013).
Therefore, it is important to elucidate the mechanisms of resistance in order to design more effective weed management strategies. We previously studied the cross-resistance patterns of Japanese foxtail populations that harbored mutations at positions 1999, 2027, and 2041 for five ACCase-inhibiting herbicides (Tang et al. Reference Tang, Li, Dong, Dong, Lü and Zhu2012; Xu et al. Reference Xu, Zhu, Wang, Li and Dong2013, Reference Xu, Li, Zhang, Cheng, Jiang and Dong2014a). All the Japanese foxtail populations harboring those three mutations were collected from eastern China. In the present study we focused on three populations with high resistance to fenoxaprop-P-ethyl also collected from eastern China. Our objectives here were to (1) identify different mutations in fenoxaprop-P-ethyl–resistant Japanese foxtail, (2) develop an effective dCAPS protocol for detecting the mutation type and frequency in fenoxaprop-P-ethyl–resistant Japanese foxtail populations, and (3) determine cross-resistance patterns to ACCase inhibitors associated with different mutations in Japanese foxtail.
Materials and Methods
Plant Material and Herbicides
Three putative resistant Japanese foxtail populations were collected in 2013 from fields in Jiangsu province, China. A reference sensitive population was collected in 2011 from a field that was never treated with any herbicide in Jiangsu province, China (Figure 1). The collected seeds were kept under dry conditions at room temperature until use. The following herbicides were used for the dose–response tests: fenoxaprop-P-ethyl, clodinafop-propargyl, haloxyfop-R-methyl, clethodim, and pinoxaden. All the herbicides used in this study were commercial formulations and are listed in Table 1.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20170717141601-22851-mediumThumb-S0043174517000066_fig1g.jpg?pub-status=live)
Figure 1 Counties with confirmed resistance to ACCase-inhibitor herbicides in Japanese foxtail in eastern China (Bi et al. Reference Bi, Liu, Guo, Li, Yuan, Du and Wang2016; Cui et al. Reference Cui, Wang, Han, Chen and Li2015; Mohamed et al. Reference Mohamed, Li, You and Li2012; Xu et al. Reference Xu, Zhu, Wang, Li and Dong2013, 2014; Yang et al. Reference Yang, Dong, Li and Moss2007), and the locations of three resistant Japanese foxtail populations used in this study.
Table 1 Herbicides used in this study.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20170717140420347-0610:S0043174517000066:S0043174517000066_tab1.gif?pub-status=live)
a Abbreviations: EC, emulsifiable concentrate; EW, water emulsion; WP, wettable powder.
Sensitivity to Fenoxaprop- P -ethyl
Whole-plant dose–response experiments were conducted under greenhouse conditions at 20/15 C (day/night), as described in Xu et al. (Reference Xu, Zhu, Wang, Li and Dong2013). At the 3- to 4-leaf stage, the plants of the susceptible population (JNXW-2) were sprayed with fenoxaprop-P-ethyl at a dose of 0.00, 0.81, 3.23, 12.94, 51.75, 207.00, or 828.00 g ai ha−1, whereas those of the putative resistant populations (JCJT-2, JZJR-1, and JCWJ-3) were sprayed at a concentration of 0.00, 3.23, 12.94, 51.75, 207.00, 828.00, or 3312.00 g ha−1. The experiment was conducted as a completely randomized design with four replications (each pot=replication, 20 plants per pot), and the experiment was conducted twice. Three weeks after treatment, the aboveground fresh weight per pot and the effective dose of herbicide causing 50% reduction of fresh weight (GR50) were determined. The resistance index (RI) was calculated as the ratio of GR50 values for the resistant and sensitive control populations.
Total DNA Extraction and Plastidic ACCase CT Domain Cloning
Total DNA was extracted from 100 mg of young shoot tissue using the Plant Genomic DNA Kit (Tiangen Biotech, Beijing, China), according to the manufacturer’s instructions. Two primer pairs were designed to amplify a 1,230-bp DNA fragment containing the entire CT domain of Japanese foxtail ACCase (Xu et al. Reference Xu, Li, Zhang, Cheng, Jiang and Dong2014a): ACCp1F/ACCp1R to amplify a 553-bp region containing codon 1781, and ACCp2F/ACCp2R to amplify an 873-bp region containing codons 1999, 2027, 2041, 2078, 2088, and 2096 (Xu et al. Reference Xu, Li, Zhang, Cheng, Jiang and Dong2014a). PCR for generating and sequencing ACCase was carried out as described by Pan et al. (Reference Pan, Li, Xia, Zhang and Dong2015). DNA samples of eight plants from each population were sequenced, and at least five clones of the PCR product for each biological replicate were used to construct the ACCase consensus sequence. Fragments were sequenced by Invitrogen Biotechnology (Shanghai, China) in both forward and reverse directions to minimize sequencing errors. The sequence data of putative resistant populations (JCJT-2, JZJR-1, and JCWJ-3) and the sensitive population (JNXW-2) were thus compared.
dCAPS for Mutation Detection and Genotype Analysis
Three new dCAPS markers were developed for detecting three different mutations (Ile-1781-Leu, Ile-2041-Asn, and Asp-2078-Gly) according to methods reported by Délye et al. (Reference Délye, Pernin and Michel2011) and Xu et al. (Reference Xu, Li, Zhang, Cheng, Jiang and Dong2014a). The primers and corresponding restriction enzymes were designed using dCAPS Finder (http://helix.wustl.edu/dcaps/dcaps.html). They are shown in Tables 2 and 3, and restriction enzymes EcoO109I (Thermo Fisher Scientific, Waltham, MA), EcoRV (Thermo Fisher Scientific), and KpnI (Thermo Fisher Scientific) were used for Ile-1781-Leu, Ile-2041-Ala, and Asp-2078-Gly mutations, respectively. PCR was performed as described by Xu et al. (Reference Xu, Li, Zhang, Cheng, Jiang and Dong2014a), and PCR products were separated on 3% agarose gels, stained with ethidium bromide, and visualized using UV.
Table 2 dCAPS primers used in this study.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20170717141601-62990-mediumThumb-S0043174517000066_tab2.jpg?pub-status=live)
a The forced mismatches introduced to create a restriction site are underlined.
Table 3 dCAPS markers and restriction enzymes.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20170717140420347-0610:S0043174517000066:S0043174517000066_tab3.gif?pub-status=live)
Sensitivity to Other ACCase-inhibiting Herbicides
Four ACCase-inhibiting herbicides representing three different chemical classes were selected for whole-plant bioassays, including two FOPs, one DIM, and one DEN. The whole-plant pot bioassays were conducted to determine the sensitivity to different herbicides. Herbicide doses (Table 4) were determined based on preliminary experiments. The experiment was conducted as a completely randomized design with four replications, and the experiment was conducted twice.
Table 4 Herbicide treatments applied for the dose–response treatments to the resistant and susceptible Japanese foxtail populations.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20170717140420347-0610:S0043174517000066:S0043174517000066_tab4.gif?pub-status=live)
Statistical Analysis
The fresh weight data were expressed as percentages of the nontreated control. Regression analysis and the three-parameter logistic function (LL.3) were performed using the ‘drc’ add-on package (Ritz and Streibig Reference Ritz and Streibig2007) in R v. 3.1.3 (R Core Team 2015) to pool and fit the data to a nonlinear logistic model (Chen et al. Reference Chen, Wang, Yao, Zhu and Dong2016; Valverde et al. Reference Valverde, Boddy, Pedroso, Eckert and Fischer2014). Thus, the effective rate of herbicide causing 50% inhibition in fresh weight (GR50) was determined for each population referring to each herbicide treated. Resistance factors (RFs) were calculated as the ratio between the GR50 of the resistant population and the GR50 of the susceptible population: S, not resistant (<2); L, low resistance (2–5); M, moderate resistance (6–10); and H, high resistance (>10) (Beckie and Tardif Reference Beckie and Tardif2012).
Results and Discussion
Sensitivity to Fenoxaprop-P-ethyl
The GR50 values of JCJT-2, JZJR-1, and JCWJ-3 were higher than the recommended field dose for fenoxaprop-P-ethyl, whereas that of JNXW-2 was lower than the field dose (Table 5). The GR50 values of JCJT-2, JZJR-1, and JCWJ-3 were 55-fold, 52-fold, and 80-fold higher, respectively, than those of JNXW-2 (Table 5; Figure 2). In China, fenoxaprop-P-ethyl has been used to control grass weeds in wheat and rapeseed fields since 1998. Two counties in Anhui province and four counties in Jiangsu province have been confirmed to have fenoxaprop-P-ethyl–resistant Japanese foxtail (Figure 1). Meanwhile, a group of grass weeds in Chinese wheat lands have also evolved fenoxaprop-P-ethyl resistance, including American sloughgrass [Beckmannia syzigachne (Steud.) Fernald] (Li et al. Reference Li, Du, Liu, Yuan and Wang2014; Pan et al. Reference Pan, Li, Xia, Zhang and Dong2015), shortawn foxtail (Alopecurus aequalis Sobol.) (Guo et al. Reference Guo, Zhang, Wang, Li, Liu and Wang2017; Xia et al. Reference Xia, Pan, Li, Wang, Feng and Dong2015), Asia Minor bluegrass (Polypogon fugax Nees ex Steud.) (Tang et al. Reference Tang, Zhou, Chen and Zhou2014), and barnyardgrass [Echinochloa crus-galli (L.) Beauv.] (Chen et al. Reference Chen, Wang, Yao, Zhu and Dong2016). Therefore, some herbicides have been introduced as alternatives to fenoxaprop-P-ethyl, such as pyroxsulam, mesosulfuron-methyl, and tralkoxydim.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20170717141601-53360-mediumThumb-S0043174517000066_fig2g.jpg?pub-status=live)
Figure 2 Fenoxaprop-P-ethyl dose–response tests for the four Japanese foxtail populations.
Table 5 Summary of fenoxaprop-P-ethyl dose–response analyses and target-site resistance (TSR) mutations for resistant and susceptible Japanese foxtail populations.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20170717140420347-0610:S0043174517000066:S0043174517000066_tab5.gif?pub-status=live)
a RF, resistance factor.
Identification of ACCase Mutations
Sequences of the entire CT domain of Japanese foxtail ACCase were obtained from all four populations. Two copies of ACCase (Acc1-1 and Acc1-2) were identified, but mutations were detected only in Acc1-1. Sequence comparison revealed that an A (JNXW-2) to C (JCJT-2) substitution resulted in an Ile (ATA) to Leu (CTA) substitution at codon position 1781 (Ile-1781-Leu); a T (JNXW-2) to A (JZJR-1) substitution resulted in an Ile (ATT) to Asn (AAT) substitution at codon position 2041 (Ile-2041-Asn); and an A (JNXW-2) to G (JCWJ-3) substitution resulted in an Asp (GAT) to Gly (GGT) substitution at codon position 2078 (Asp-2078-Gly).
Multiple copies of genes encoding plastidic ACCase have been identified in rice barnyardgrass [Echinochloa phyllopogon (Stapf) Koso-Pol.], sterile oat (Avena sterilis L.), Japanese foxtail, and wild oat (Avena fatua L.) (Christoffers et al. Reference Christoffers, Berg and Messersmith2002; Iwakami et al. Reference Iwakami, Uchino, Watanabe, Yamasue and Inamura2012; Xu et al. Reference Xu, Zhang, Zhang, Li, Wu and Dong2014b; Yu et al. Reference Yu, Ahmad-Hamdani, Han, Christoffers and Powles2013), and it has been demonstrated that any of the three Acc1 homologues (Acc1-1, Acc1-2, and Acc1-3) in wild oat and two Acc1 homologues (Acc1-1 and Acc1-2) in Japanese foxtail, can also host ACCase resistance mutations (Christoffers et al. Reference Christoffers, Berg and Messersmith2002; Xu et al. Reference Xu, Li, Zhang, Cheng, Jiang and Dong2014a; Yu et al. Reference Yu, Ahmad-Hamdani, Han, Christoffers and Powles2013). Mutations Ile-1781-Leu and Ile-2041-Asn in Japanese foxtail have been reported (Bi et al. Reference Bi, Liu, Guo, Li, Yuan, Du and Wang2016; Cui et al. Reference Cui, Wang, Han, Chen and Li2015; Mohamed et al. Reference Mohamed, Li, You and Li2012; Tang et al. Reference Tang, Li, Dong, Dong, Lü and Zhu2012), but the Asp-2078-Gly mutation in Japanese foxtail has not been reported.
dCAPS for Mutation Detection and Genotype Analysis
The dCAPS method detected these mutations successfully in Japanese foxtail. Of the 96 plants tested in JCJT-2, two were wild type and 94 were mutants with the Ile-1781-Leu substitution (Table 6); of the 61 plants tested in JZJR-1, five were wild type and 56 were mutants with the Ile-2041-Asn substitution; and of the 23 plants tested in JCWJ-3, three were wild type and 20 were mutant with the Asp-2078-Gly substitution. The mutation frequencies in JCJT-2, JZJR-1, and JCWJ-3 were approximately 98%, 92%, and 87%, respectively. dCAPS is a high-throughput PCR-based method to detect known resistance-endowing mutations, and this method was developed to easily determine resistance caused by mutations of ACCase or ALS in different weed species (Délye et al. Reference Délye, Pernin and Michel2011; Pan et al. Reference Pan, Li, Xia, Zhang and Dong2015).
Table 6 The derived cleaved amplified polymorphic sequence (dCAPS) detection of mutations in the resistant populations of Japanese foxtail.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20170717140420347-0610:S0043174517000066:S0043174517000066_tab6.gif?pub-status=live)
Cross-resistance to Other ACCase-inhibiting Herbicides
JCJT-2 (Ile-1781-Leu) was highly resistant to clodinafop-propargyl, slightly resistant to haloxyfop-R-methyl and pinoxaden, but susceptible to clethodim (Table 7; Figure 3). JZJR-1 (Ile-2041-Asn) was highly resistant to clodinafop-propargyl and haloxyfop-R-methyl, slightly resistant to pinoxaden, but susceptible to clethodim (Table 7; Figure 3). Previous studies have shown that the Ile-1781-Leu substitution confers high resistance to FOPs, DIMs, and pinoxaden in blackgrass (Alopecurus myosuroides Huds.), wild oat, sterile oat, green foxtail [Setaria viridis (L.) Beauv.], and rigid ryegrass (Lolium rigidum Gaudin) (Beckie and Tardif Reference Beckie and Tardif2012; Brown et al. Reference Brown, Moss, Wilson and Field2002; Christoffers et al. Reference Christoffers, Berg and Messersmith2002; Délye et al. Reference Délye, Pernin and Michel2011; Tal and Rubin Reference Tal and Rubin2004). In the present study, JCJT-2, harboring the Ile-1781-Leu substitution, showed resistance to haloxyfop-R-methyl, clodinafop-propargyl, and pinoxaden, and was susceptible to clethodim. Our previous study found that American sloughgrass with the Ile-1781-Leu mutation showed resistance to all ACCase-inhibiting herbicides tested; while the cross-resistance varied for different DIMs: specifically, high resistance (RF=31) to sethoxydim and low resistance (RF=5) to clethodim (Pan et al. Reference Pan, Li, Xia, Zhang and Dong2015). Yu et al. (Reference Yu, Ahmad-Hamdani, Han, Christoffers and Powles2013) also found the Ile-1781-Leu mutation in wild oat ACCase confers high resistance to sethoxydim (RI=11) and low resistance to clethodim (RI=4). Moreover, two blackgrass populations harboring Ile-1781-Leu mutation collected in Denmark were found to be sensitive (RF=1.3) or to have low resistance (RF=2.4) to a DIM herbicide, cycloxydim (Keshtkar et al. Reference Keshtkar, Mathiassen, Moss and Kudsk2015). It is reported that blackgrass, hood canarygrass (Phalaris paradoxa L.), and rigid ryegrass with the Ile-2041-Asn substitution were highly resistant to FOPs, but not to DIMs, and showed either no resistance or moderate resistance to DENs (Beckie and Tardif Reference Beckie and Tardif2012; Hochberg et al. Reference Hochberg, Sibony and Rubin2009). These findings were consistent with those reported in the present study. JZJR-1, which harbors the Ile-2041-Asn substitution, was highly resistant to the two FOPs, slightly resistant to pinoxaden (DEN herbicide), but susceptible to clethodim (DIM herbicide), indicating that the mutation confers high resistance to FOPs, but not to DIMs or DENs.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20170717141601-75138-mediumThumb-S0043174517000066_fig3g.jpg?pub-status=live)
Figure 3 Dose-dependent response of the four Japanese foxtail populations to four ACCase-inhibiting herbicides.
Table 7 The sensitivities of the resistant and sensitive Japanese foxtail populations to ACCase-inhibiting herbicides.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20170717141601-34627-mediumThumb-S0043174517000066_tab7.jpg?pub-status=live)
a RF, resistance factor.
JCWJ-3 (Asp-2078-Gly) was highly resistant to all four ACCase-inhibiting herbicides, and the RFs for clodinafop-propargyl and haloxyfop-R-methyl, pinoxaden and clethodim were 46.85, 10.59, 46.50, and 12.37, respectively (Table 7; Figure 3). Previous studies have shown that blackgrass, rigid ryegrass, sterile oat, and hood canarygrass harboring the Asp-2078-Gly substitution are resistant to all three classes of ACCase-inhibiting herbicides (Beckie and Tardif Reference Beckie and Tardif2012; Liu et al. Reference Liu, Harrison, Chalupska, Gornicki, O’Donnell, Adkins, Haselkorn and Williams2007; Yu et al. Reference Yu, Collavo, Zheng, Owen, Sattin and Powles2007).
Conclusions
An effective dCAPS protocol was developed to detect the resistance-endowing mutations in the ACCase gene of Japanese foxtail. Ile-1781-Leu, Ile-2041-Asn, and Asp-2078-Gly were detected in the three Japanese foxtail populations studied. JCJT-2 (Ile-1781-Leu) and JZJR-1 (Ile-2041-Asn) showed cross-resistance to haloxyfop-R-methyl, clodinafop-propargyl, and pinoxaden, but not to clethodim. JCWJ-3 (Asp-2078-Gly) showed cross-resistance to all tested ACCase-inhibiting herbicides. Our results might allow farmers to select the most appropriate herbicides to manage fenoxaprop-P-ethyl–resistant Japanese foxtail populations.
Acknowledgments
This study was funded by the Natural Science Foundation for Young Scientists of Jiangsu province (BK20160724) and the Fundamental Research Funds for the Central Universities (KYZ201511). We thank Jun Li and Di Zhang from Nanjing Agricultural University for their help. Thanks are also due to the reviewers and editors for their helpful comments on earlier drafts of the article.