0 引言
【研究意义】棉花黄萎病是危害棉花最大的世界性病害之一,被称为棉花的“癌症”。实践证明,防治棉花黄病害最经济、安全、有效的方法就是培育和推广抗病品种[1]。抗病品种的选育已成为目前棉花育种的主要方向。然而,近年来,中国棉花抗黄萎病育种进展缓慢,主要原因是中国现有棉花品种资源的遗传基因狭窄,缺乏高抗黄萎病的棉花资源,与病原菌小种的快速变异不相适应[2-3]。在棉花抗黄萎病性状的筛选上,目前,有多种方法对黄萎病抗性进行鉴定[4],其中,最常用的鉴定方法是在大田环境条件下利用黄萎病病圃对育种材料进行鉴定。然而,由于大田鉴定条件下存在植株成熟早晚和接种量的差异,以及试验小区的干扰,导致目前棉花育种中对抗黄萎病性状的选择效率十分有限,表现为不同年份、不同棉区筛选的抗黄萎病材料的抗病性表现存在较大差异,难以筛选到真正的抗黄萎病材料。因此,如何克服环境条件对抗病性鉴定的影响,实现对抗病基因型的直接鉴定与选择,对棉花抗黄萎病品种的培育具有重要意义。【前人研究进展】近年来,随着分子标记技术的发展,利用分子标记筛选与目标性状基因紧密连锁的分子标记,可以大大提高性状选择的准确性,并且不受环境因子影响。目前,国内外已经定位的与棉花不同生长阶段的抗黄萎病性状相关的分子标记或QTL已超过了100个[5-17],利用这些分子标记进行辅助选择,可以有效地提高性状改良的选育效率,加快育种进程。在棉花抗病基因分子标记辅助选择育种上,王芙蓉等[18]利用与检测到的3个抗黄萎病QTL较紧密连锁的SSR标记对杂交后代辅助鉴定,标记BNL3590在纯合的父、母本基因型后代中,平均病指无显著差异;而标记NAU751和BNL1395的抗性基因型均能显著增加后代的黄萎病抗性,2个标记的抗性基因型聚合后,后代抗性水平极显著提高。孔祥瑞等[19]将与陆地棉黄萎病抗性相关的18个SSR标记用于大田辅助选择育种,结果表明,18个标记中的BNL1721、BNL2733和BNL3452标记在有/无标记基因型个体间病指差异显著,且选择效应能够在不同世代间稳定遗传。同时,随着遗传学研究的发展,基于连锁不平衡（LD）的关联分析在植物基因组学研究领域得到了成功的应用[20],为植物数量性状的研究开辟了新的途径。与QTL作图相比,关联分析有以下优点：以自然群体为材料,无需构建作图群体,因此,可以检测到更多的等位基因变异,从而挖掘到更多的与性状关联的分子标记,为基于分子标记的性状筛选奠定基础。【本研究切入点】由于目前在棉花中,利用分子标记对棉花抗黄萎病性状进行选择大多是对来自抗、感亲本杂交的后代分离群体进行选择,这些分子标记仅仅是作为一种辅助手段,仅能用于区分抗病基因型和感病基因型,无法用于多个材料间抗病性大小的比较。【拟解决的关键问题】本研究在前期对棉花抗黄萎病性状进行关联分析的基础上,进一步发掘与棉花抗黄萎病性状显著关联的优异等位基因位点,研究了利用优异等位基因位点进行棉花抗黄萎病性状预测的可行性,并分析了基于优异等位基因位点的抗病性分子鉴定与传统抗病性表型鉴定的相关性。该研究将为棉花抗黄萎病性状的分子检测及标记辅助选择育种提供重要的信息,有助于实现对抗病性的快速鉴定和对抗病基因型的直接选择,有效解决当前抗病性鉴定周期长、抗病性选择效率低的技术难题。1 材料与方法
1.1 棉花种质材料
棉花材料为125份陆地棉优异种质,这些材料均来自于笔者早期构建的由158份陆地棉优异种质组成的关联群体[21],在前期研究基础上,通过排除掉基因型数据或表型鉴定数据存在缺失的部分材料,最终确定了125材料用于本研究,这些材料的名称、编号及系谱来源见电子附表1。其中,材料的编号与前期研究所用的编号一致。1.2 抗病性鉴定
对125份材料的抗病性鉴定分别采用人工黄萎病病圃鉴定和温室接种鉴定2种方法,其中,人工黄萎病病圃鉴定采用国家棉花品种区试中抗黄萎病的鉴定方法[22],温室鉴定采用营养钵底部菌液蘸根法接种鉴定[23],对鉴定方法的详细描述请参考笔者以前的研究[21]。抗病性鉴定所用的黄萎病病菌均为中等致病力的安阳菌系—VD080,所用的感病对照为冀棉11,采用相对病情指数（relative disease index,RDI）表示各个材料的抗病性强弱,相对病情指数的计算方法参考文献[22]。1.3 优异等位基因的筛选及其效应值计算
在前期研究中,利用均匀分布于棉花全基因组的212对分子标记对158份陆地棉优异种质进行基因型分析,同时利用人工黄萎病病圃鉴定和温室鉴定2种方法对群体进行抗病表型鉴定,在此基础上对棉花抗黄萎病性状进行了全基因组关联分析研究,获得了40个与抗黄萎病表型显著关联的SSR标记位点[21]。将这些SSR标记序列与陆地棉基因组序列（测序版本号：BGI_Gossypium_hirsutum_v1.0,https://www. cottongen.org/species/Gossypium_hirsutum/bgi-AD1_genome_v1.0）进行BLAST比对,根据最好的比对结果（E值最低）确定SSR标记的物理位置,从而将这些SSR标记位点锚定在陆地棉染色体上。进一步分析了每个标记位点的各个等位基因在关联群体中的分布情况,根据ZHANG等[24]的描述,利用公式ai=∑xij/ni–∑Nk/nk计算等位基因的效应值,其中,ɑi是第i个等位基因的效应值,xij是含有第i个等位基因的第j个材料的表型值,ni是含有第i个等位基因的材料数目,Nk是所有材料的表型值,nk是所有材料的数目。将效应值最低的等位基因（黄萎病病指最低）确定为该标记位点的优异等位基因。1.4 基因扩增与检测
利用优异等位基因所对应的引物序列对待测的棉花材料进行PCR扩增,PCR反应体系10 µL,其中,含模板DNA 0.1 µg、10?PCR buffer（含Mg2+）1 µL、dNTP（10 nmol·µL-1）0.2 µL、Tag酶（5 U·µL-1）0.1 µL、上下游引物各1 µM。PCR反应程序为95℃ 3 min;94℃ 45 s,56℃ 45 s,72℃ 1 min,30个循环;72℃ 5 min。扩增产物用浓度10%的垂直板非变性聚丙烯酰胺凝胶电泳分离,电泳结束参考张军等[25]的方法银染检测。1.5 数据分析
根据基因扩增结果,对待测的棉花材料进行分子标记基因型分析,根据基因型分析结果统计每个材料含有的优异等位基因及其数目,根据每个优异等位基因的效应值,计算每份材料含有的所有优异等位基因的效应值之和;利用每份材料含有的优异等位基因的数目及其效应值的和来评价材料的抗病性强弱,实现不同棉花材料间抗黄萎病性状强弱的判断。2 结果
2.1 棉花优异种质的黄萎病抗性
对125份陆地棉材料进行抗病性鉴定表明,陆地棉在黄萎病抗性上表现出广泛的变异。其中,在田间病圃鉴定条件下,125份种质材料的抗黄萎病相对病指变化范围为10.10—76.6,平均为37.69;而在温室鉴定条件下,125份种质材料的抗黄萎病相对病指变化范围为17.01—72.63,平均为38.51（图1）。方差分析表明,不同材料的相对病情指数间存在极显著差异（P<0.001）,说明了这125份材料的抗病性变异极为丰富;而基于温室鉴定条件和基于田间病圃鉴定条件的相对病情指数（RDI）间不存在显著差异（P=0.5054）,即环境效应不显著,说明了温室鉴定和田间病圃鉴定的一致性。进一步比较发现,在温室鉴定条件下表现为耐病的材料（相对病指在20—35）有56份,在这56份材料中,在田间病圃条件下仍然表现为耐病的有36份,表现为感病的材料（相对病指在大于35）的有19份,表现为抗病的材料（相对病指在小于20）仅有1份;另外,在温室鉴定条件下相对病情指数最高的7份材料（相对病情指数在57.23—72.63）,在田间病圃鉴定条件下仍然表现为高度感病（相对病指在44.8—66）;然而,温室鉴定为抗病的材料仅有1份（相对病指17.1）,该材料在病圃条件下表现为耐病（相对病指29.2）。2.2 抗黄萎病优异等位基因的效应值
前期研究中,在田间病圃鉴定条件下获得30个与抗黄萎病性状显著关联的SSR标记位点,在温室鉴定条件下获得14个与抗黄萎病性状显著关联的SSR标记位点,同时获得在2个鉴定环境中均存在的抗黄萎病性位点4个。通过统计每个标记位点的各个等位基因在关联群体中的分布情况,利用公式ai=∑xij/ni–∑Nk/nk计算各个等位基因的效应值,将效应值最低的等位基因（黄萎病病指最低）确定为该标记位点的优异等位基因。40个标记位点共获得了40个优异等位基因,效应值的变化范围为-8.20—-0.39。其中,田间病圃鉴定条件下获得30个优异等位基因,效应值的变化范围为-8.2—-0.86;温室鉴定条件下获得14个优异等位基因,效应值的变化范围为-6.0—-0.39;同时获得在2个鉴定环境中均存在的优异等位基因4个。这些优异等位基因的基因组位置信息及其效应值如表1所示。
显示原图|下载原图ZIP|生成PPT图1田间病圃鉴定和温室鉴定条件下供试材料病指变化频率分布图
-->Fig. 1Histogram of relative disease index of tested lines identified in field disease nursery and in greenhouse
-->
Table 1
表1
表1优异等位基因的基因组位置信息及其效应值
Table 1Physical position and effective values of elite alleles
| 优势等位基因名称 Name of elite allele | 扩增产物大小 Size in PCR (bp) | 所在染色体Chromosome | 物理位置c Physical position | 基因效应值 Gene effective value | 代表性材料d Representatives | ||
|---|---|---|---|---|---|---|---|
| 序列起点 Start (bp) | 序列终点 End (bp) | 病圃a Disease nursery | 温室b Greenhouse | ||||
| NAU5233-2 | 590 | At_chr10 | 68085564 | 68086389 | -0.86 | 169, 317 | |
| NAU3592-2 | 185 | N | N | N | -3.68 | 317, 243 | |
| NAU3828-2e | 310 | Dt_chr9 | 66326715 | 66327125 | -2.73 | -1.56 | 169, 243 |
| NAU3212-1 | 360 | Dt_chr9 | 57900120 | 57900337 | -2.11 | 317, 243 | |
| DPL0222-1 | 580 | At_chr6 | 5277130 | 5276426 | -7.16 | -4.60 | 317, 243 |
| NAU3074-2 | 190 | Dt_chr7 | 3881156 | 3880922 | -2.35 | 169, 243 | |
| CIR196-1 | 300 | N | N | N | -2.59 | 169, 264 | |
| NAU980-2 | 180 | Dt_chr9 | 42307868 | 42307475 | -4.68 | 169, 317 | |
| NAU5428-3 | 490 | At_chr1 | 3327910 | 3328399 | -8.20 | 44, 209 | |
| BNL1034-2 | 400 | At_chr7 | 16189113 | 16188837 | -2.75 | 169, 317 | |
| NAU5064-2 | 310 | Dt_chr9 | 42242900 | 42242244 | -4.95 | 264, 185 | |
| BNL2646-1 | 130 | Dt_chr2 | 4194117 | 4194526 | -6.77 | 317, 85 | |
| BNL2441-1 | 250 | Dt_chr1 | 61687093 | 61687615 | -1.56 | 169, 243 | |
| BNL3319-2 | 210 | Dt_chr1 | 28243618 | 28243007 | -5.51 | 169, 185 | |
| TMB1114-2 | 220 | Dt_chr1 | 83322762 | 83322223 | -6.17 | 264, 44 | |
| NAU2887-2 | 460 | Dt_chr1 | 45872109 | 45871866 | -7.42 | 169, 317 | |
| NAU5120-2 | 170 | Dt_chr1 | 40145540 | 40145267 | -2.53 | 317, 264 | |
| BNL1606-1 | 220 | At_chr5 | 4066248 | 4065910 | -6.14 | -5.28 | 85, 167 |
| NAU2859-1 | 280 | Dt_chr5 | 17031760 | 17031155 | -3.89 | 169, 317 | |
| JESPR101-2 | 140 | Dt_chr3 | 38672549 | 38672780 | -5.06 | 167, 165 | |
| CIR364-2 | 250 | Dt_chr9 | 55697292 | 55697037 | -1.24 | 85, 137 | |
| NAU2894-1 | 190 | At_chr10 | 25818746 | 25818852 | -2.41 | 169, 317 | |
| BNL3649-2 | 200 | At_chr1 | 11226368 | 11225811 | -1.82 | -2.13 | 169, 317 |
| JESPR274-2 | 190 | Dt_chr6 | 21342913 | 21343016 | -0.90 | 72, 168 | |
| NAU1047-1 | 175 | Dt_chr6 | 6035428 | 6035210 | -1.96 | 169, 317 | |
| NAU4912-1 | 360 | Dt_chr8 | 41697924 | 41697169 | -4.27 | 169, 185 | |
| NAU5463-2 | 190 | At_chr10 | 70854329 | 70854641 | -2.32 | 169, 317 | |
| Gh268-1 | 190 | Dt_chr11 | 41841150 | 41840459 | -3.83 | 169, 185 | |
| Gh454-1 | 180 | Dt_chr11 | 41840459 | 41841121 | -3.99 | 169, 264 | |
| NAU3563-2 | 220 | Dt_chr1 | 105137054 | 105137323 | -1.68 | 317, 243 | |
| BNL2599-1 | 120 | N | N | N | -1.20 | 138, 132 | |
| BNL3255-2 | 240 | At_chr10 | 31849135 | 31849432 | -1.29 | 138, 132 | |
| NAU3201-2 | 240 | N | N | N | -0.64 | 138, 132 | |
| NAU3499-3 | 550 | Dt_chr10 | 6778862 | 6779369 | -2.28 | 138, 132 | |
| NAU2627-1 | 460 | Dt_chr1 | 46597406 | 46597039 | -3.31 | 124, 91 | |
| BNL4069-1 | 100 | Dt_chr9 | 26272316 | 26271919 | -1.42 | 138, 132 | |
| JESPR0001-1 | 500 | Dt_chr9 | 32274761 | 32274282 | -0.39 | 138, 287 | |
| BNL3646-1 | 230 | N | N | N | -6.00 | 287, 73 | |
| NAU3574-2 | 210 | Dt_chr11 | 31218783 | 31219254 | -1.34 | 138, 132 | |
| NAU2954-1 | 460 | Dt_chr6 | 7764989 | 7765806 | -1.55 | 138, 132 | |
新窗口打开
2.3 基于优异等位基因的抗病性强弱判断
统计每份材料含有的优异等位基因及其数目, 根据每个优异等位基因的效应值,计算每份材料含有的所有优异等位基因的效应值之和,该效应值之和代表了所有优异等位基因聚合情况下的基因效应。利用40个优异等位基因的分子标记对125份材料进行分析表明,每份材料含有的优异等位基因数目为1—24个,125份材料优异等位基因效应值之和的变化范围为-92.37—-0.86（表2）。其中,编号为191的材料（澳152）含有24个优异等位基因,其优异等位基因效应值之和为-92.37,在病圃条件下鉴定出的黄萎病相对病指为26.9,为耐病材料;编号为287的材料（新研96-48）仅含有1个优异等位基因,其优异等位基因效应值之和为-0.86,在病圃条件下鉴定出的黄萎病相对病指为46.5,为感病材料。在病圃鉴定条件下,优异等位基因效应值之和小于-32.6（平均值）的材料共56个,其中有47份材料（83.9%）的抗黄萎病病指在35以下,其中包括5份病指20以下的抗病材料(表3),这5份材料的名称分别为M-8124-1159（编号169）、中Arc-185（编号185）、GK99-1（编号243）、鲁458（编号264）和99633（编号317）,所含有的优异等位基因数目分别为16、14、12、11和15个,优异等位基因效应值之和分别为-50.27、-42.06、-34.10、-35.63和-52.43,在病圃条件下鉴定出的黄萎病相对病指分别为10.1、19.8、17.3、17.6和14.7 （表2）;优异等位基因效应值之和大于-32.6（平均值）的材料共69份,其中有49份材料（71%）的抗黄萎病病指在35以上（表2）,例如,高度感病的5份材料大铃棉（编号67）、常德184（编号5）、邯8944（编号123）、苏联21（编号46）和淮棉4（编号9）中,优异等位基因数目分别为6、4、7、8和6个,优异等位基因效应值之和分别为-12.69、-7.56、-13.97、-20.54和-12.21,在病圃条件下鉴定出的黄萎病相对病指分别为76.6、73.3、66、60.4和59.6。这说明利用优异等位基因的效应值之和可以有效的将抗（耐）病材料（黄萎病指35以下）与感病材料（黄萎病指35以上）区分开,从而实现比较不同材料抗病性强弱的目的。Table 2
表2
表2125材料的优异等位基因数目、效应值之和及抗黄萎病相对病指
Table 2The number and sum of effective values of elite alleles, and the relative disease index of Verticillium wilt resistance of 125 lines
| 材料编号 Line ID | 病圃鉴定 Tested in disease nursery | 温室鉴定 Tested in greenhouse | |||||
|---|---|---|---|---|---|---|---|
| 优异等位基因数目 Number of elite alleles | 优异等位基因效应值之和 Sum of effective value of elite alleles | 抗黄萎病相对病指 Relative disease index | 优异等位基因数目 Number of elite alleles | 优异等位基因效应值之和 Sum of effective value of elite alleles | 抗黄萎病相对病指 Relative disease index | ||
| 5 | 4 | -7.56 | 73.30 | 8 | -10.34 | 32.09 | |
| 8 | 9 | -25.11 | 53 | 10 | -17.52 | 28.85 | |
| 9 | 6 | -12.21 | 59.60 | 7 | -7.83 | 34.20 | |
| 12 | 13 | -41.65 | 35.2 | 7 | -8.93 | 47.00 | |
| 13 | 10 | -19.77 | 51.5 | 8 | -10.95 | 46.17 | |
| 15 | 11 | -30.80 | 49.2 | 9 | -12.24 | 38.53 | |
| 16 | 7 | -15.52 | 54.9 | 8 | -10.95 | 48.83 | |
| 22 | 12 | -28.38 | 53.7 | 2 | -1.96 | 71.97 | |
| 24 | 9 | -27.04 | 59.0 | 5 | -5.87 | 51.03 | |
| 26 | 11 | -31.11 | 26.2 | 7 | -8.97 | 47.15 | |
| 27 | 8 | -21.84 | 50.5 | 5 | -5.31 | 43.01 | |
| 32 | 16 | -43.75 | 33.1 | 9 | -12.24 | 32.66 | |
| 44 | 13 | -44.27 | 21.3 | 8 | -9.97 | 32.83 | |
| 45 | 14 | -45.24 | 25.7 | 9 | -16.39 | 55.77 | |
| 46 | 8 | -20.54 | 60.4 | 9 | -15.42 | 24.64 | |
| 48 | 12 | -32.04 | 27.5 | 5 | -5.87 | 52.51 | |
| 51 | 16 | -57.44 | 26.4 | 10 | -18.51 | 26.83 | |
| 54 | 8 | -22.26 | 44.8 | 5 | -6.15 | 67.36 | |
| 55 | 8 | -16.44 | 46.5 | 7 | -13.63 | 32.32 | |
| 59 | 7 | -13.94 | 52.0 | 9 | -13.44 | 47.64 | |
| 62 | 9 | -25.52 | 51.0 | 8 | -15.86 | 48.97 | |
| 63 | 9 | -23.80 | 46.7 | 7 | -12.50 | 53.95 | |
| 66 | 9 | -26.20 | 51.4 | 6 | -9.20 | 48.39 | |
| 67 | 6 | -12.69 | 76.6 | 6 | -7.45 | 41.73 | |
| 69 | 3 | -5.34 | 50.9 | 9 | -17.85 | 31.61 | |
| 70 | 10 | -23.38 | 52.8 | 7 | -9.01 | 56.39 | |
| 73 | 15 | -45.82 | 23.1 | 10 | -18.51 | 21.59 | |
| 76 | 17 | -52.47 | 26.5 | 10 | -18.25 | 35.97 | |
| 77 | 8 | -22.08 | 26.6 | 7 | -9.32 | 56.75 | |
| 84 | 10 | -32.34 | 26.9 | 7 | -9.80 | 42.45 | |
| 89 | 10 | -29.65 | 51.3 | 8 | -10.12 | 31.18 | |
| 90 | 8 | -18.58 | 50.7 | 6 | -8.43 | 35.49 | |
| 91 | 6 | -13.34 | 46.4 | 7 | -11.70 | 23.09 | |
| 94 | 11 | -31.56 | 51.8 | 8 | -11.86 | 36.96 | |
| 96 | 15 | -40.41 | 46.7 | 9 | -12.25 | 34.29 | |
| 98 | 4 | -12.54 | 58.5 | 6 | -6.50 | 39.89 | |
| 99 | 11 | -33.66 | 46.4 | 8 | -10.82 | 46.33 | |
| 103 | 6 | -11.11 | 48.0 | 9 | -15.96 | 30.80 | |
| 106 | 12 | -40.87 | 37.1 | 10 | -19.88 | 23.25 | |
| 107 | 14 | -46.35 | 32.2 | 10 | -20.58 | 41.07 | |
| 108 | 13 | -43.68 | 27.6 | 7 | -13.94 | 41.41 | |
| 续表2 Continued table 2 | |||||||
| 材料编号 Line ID | 病圃鉴定 Tested in disease nursery | 温室鉴定 Tested in greenhouse | |||||
| 优异等位基因数目 Number of elite alleles | 优异等位基因效应值之和 Sum of effective value of elite alleles | 抗黄萎病相对病指 Relative disease index | 优异等位基因数目 Number of elite alleles | 优异等位基因效应值之和 Sum of effective value of elite alleles | 抗黄萎病相对病指 Relative disease index | ||
| 111 | 15 | -47.94 | 32.1 | 7 | -10.43 | 53.86 | |
| 113 | 12 | -29.98 | 27.4 | 8 | -11.99 | 32.66 | |
| 116 | 12 | -35.87 | — | 7 | -9.27 | 30.04 | |
| 117 | 14 | -54.73 | 25.8 | 9 | -15.65 | 30.92 | |
| 121 | 7 | -13.24 | 57.7 | 9 | -12.24 | 50.44 | |
| 123 | 7 | -13.97 | 66.0 | 6 | -8.56 | 57.23 | |
| 124 | 13 | -37.50 | 26.9 | 11 | -21.11 | 21.69 | |
| 127 | 7 | -17.34 | 48.2 | 6 | -10.50 | 41.46 | |
| 128 | 10 | -29.95 | 49.1 | 7 | -10.75 | 42.75 | |
| 130 | 15 | -47.96 | 45.2 | 11 | -20.12 | 27.58 | |
| 132 | 6 | -12.44 | 29.9 | 9 | -13.42 | 20.42 | |
| 135 | 9 | -21.81 | 27.8 | 7 | -9.15 | 26.38 | |
| 136 | 15 | -50.69 | 28.7 | 10 | -13.81 | 29.87 | |
| 137 | 17 | -62.83 | 23.8 | 11 | -20.27 | 39.83 | |
| 138 | 17 | -53.45 | 29.2 | 11 | -19.08 | 17.01 | |
| 139 | 15 | -48.46 | 27.4 | 12 | -22.40 | 40.66 | |
| 140 | 10 | -28.74 | 27.8 | 10 | -16.48 | 37.73 | |
| 141 | 11 | -30.41 | 26.7 | 10 | -17.52 | 34.02 | |
| 143 | 10 | -32.40 | 29.3 | 8 | -10.91 | 39.01 | |
| 144 | 12 | -40.50 | 32.7 | 9 | -13.43 | 34.84 | |
| 150 | 12 | -38.84 | 23.6 | 8 | -10.91 | 35.95 | |
| 151 | 13 | -36.29 | 22.7 | 9 | -16.32 | 25.61 | |
| 158 | 15 | -50.69 | 30.2 | 9 | -12.47 | 34.89 | |
| 161 | 11 | -24.18 | 27.5 | 6 | -7.47 | 54.25 | |
| 163 | 3 | -8.11 | 27.6 | 3 | -3.82 | 31.83 | |
| 164 | 14 | -46.98 | 26.4 | 7 | -8.63 | 30.06 | |
| 166 | 10 | -29.44 | 31.9 | 9 | -12.25 | 28.66 | |
| 168 | 13 | -41.11 | 23.8 | 4 | -3.65 | 46.08 | |
| 169 | 16 | -50.27 | 10.1 | 8 | -11.32 | 42.91 | |
| 171 | 16 | -52.72 | 26.2 | 8 | -14.45 | 38.65 | |
| 172 | 11 | -32.97 | 27.2 | 9 | -15.57 | 36.29 | |
| 173 | 8 | -21.44 | 30.1 | 10 | -18.24 | 29.43 | |
| 174 | 13 | -39.59 | 22.7 | 10 | -13.81 | 27.42 | |
| 179 | 13 | -38.66 | 30.6 | 9 | -15.53 | 38.45 | |
| 180 | 8 | -26.70 | 51.3 | 8 | -12.44 | 46.68 | |
| 182 | 11 | -25.60 | 32.5 | 9 | -13.55 | 24.48 | |
| 185 | 14 | -42.06 | 19.8 | 10 | -18.51 | 43.54 | |
| 186 | 8 | -20.79 | 47.5 | 9 | -13.17 | 32.72 | |
| 189 | 10 | -23.41 | 45.5 | 9 | -14.36 | 30.24 | |
| 191 | 24 | -92.37 | 26.9 | 10 | -16.85 | 24.63 | |
| 199 | 5 | -10.10 | 29.1 | 10 | -18.24 | 34.06 | |
| 201 | 5 | -16.13 | 52.4 | 10 | -17.52 | 30.93 | |
| 207 | 7 | -20.63 | 29.7 | 6 | -7.14 | 56.39 | |
| 续表2 Continued table 2 | |||||||
| 材料编号 Line ID | 病圃鉴定 Tested in disease nursery | 温室鉴定 Tested in greenhouse | |||||
| 优异等位基因数目 Number of elite alleles | 优异等位基因效应值之和 Sum of effective value of elite alleles | 抗黄萎病相对病指 Relative disease index | 优异等位基因数目 Number of elite alleles | 优异等位基因效应值之和 Sum of effective value of elite alleles | 抗黄萎病相对病指 Relative disease index | ||
| 209 | 12 | -39.66 | 22.3 | 10 | -16.85 | 40.14 | |
| 211 | 15 | -50.60 | 28.2 | 9 | -15.65 | 29.64 | |
| 213 | 13 | -30.99 | 30.4 | 8 | -11.75 | 31.06 | |
| 216 | 10 | -27.87 | 48.2 | 8 | -10.92 | 72.63 | |
| 218 | 5 | -13.11 | 51.3 | 9 | -12.25 | 34.86 | |
| 222 | 14 | -48.64 | 29.7 | 9 | -15.65 | 28.40 | |
| 227 | 6 | -18.10 | 51.3 | 7 | -8.41 | 36.89 | |
| 229 | 7 | -20.15 | 55.9 | 5 | -5.21 | 64.63 | |
| 235 | 16 | -50.42 | 24.2 | 12 | -22.40 | 30.79 | |
| 237 | 17 | -48.30 | 25.1 | 10 | -13.81 | 36.49 | |
| 238 | 9 | -29.40 | 49.2 | 9 | -12.47 | 50.00 | |
| 240 | 16 | -50.31 | 36.1 | 9 | -12.38 | 40.17 | |
| 241 | 10 | -24.84 | 27.4 | 8 | -11.99 | 23.62 | |
| 243 | 12 | -34.10 | 17.3 | 10 | -16.99 | 37.56 | |
| 245 | 18 | -62.29 | 25.4 | 11 | -18.41 | 29.31 | |
| 253 | 13 | -37.08 | 27.6 | 10 | -17.53 | 27.78 | |
| 255 | 14 | -48.34 | 27.5 | 9 | -18.08 | 41.03 | |
| 261 | 6 | -14.24 | 53.2 | 8 | -10.19 | 36.80 | |
| 262 | 9 | -25.16 | 43.4 | 10 | -15.83 | 34.87 | |
| 263 | 13 | -43.20 | 31.7 | 8 | -14.36 | 30.23 | |
| 264 | 11 | -35.63 | 17.6 | 6 | -7.49 | 50.89 | |
| 265 | 13 | -37.84 | 50.4 | 8 | -9.96 | 45.26 | |
| 266 | 12 | -39.09 | 48.5 | 1 | -1.20 | 67.89 | |
| 269 | 8 | -24.52 | 52.3 | 6 | -9.09 | 41.79 | |
| 270 | 10 | -24.58 | 29.7 | 6 | -9.36 | 37.81 | |
| 276 | 5 | -10.19 | 50.5 | 7 | -9.71 | 53.79 | |
| 277 | 8 | -18.80 | 47.4 | 8 | -10.24 | 58.88 | |
| 280 | 12 | -28.41 | 45.1 | 0 | 0.00 | 43.59 | |
| 281 | 4 | -7.76 | 45.4 | 6 | -7.40 | 42.15 | |
| 282 | 12 | -44.00 | 46.2 | 6 | -8.20 | 38.62 | |
| 287 | 1 | -0.86 | 46.5 | 8 | -14.78 | 21.01 | |
| 288 | 8 | -25.90 | 46.6 | 8 | -14.82 | 26.50 | |
| 291 | 17 | -61.79 | 36.9 | 2 | -5.67 | 48.16 | |
| 295 | 6 | -18.34 | 45.3 | 5 | -7.78 | 43.14 | |
| 299 | 16 | -60.93 | 25.1 | 9 | -15.74 | 26.85 | |
| 307 | 16 | -63.68 | 25.3 | 8 | -11.99 | 50.20 | |
| 308 | 15 | -45.82 | 25.1 | 8 | -16.70 | 25.01 | |
| 311 | 15 | -50.60 | 27.9 | 7 | -13.70 | 30.47 | |
| 312 | 8 | -17.42 | 28.4 | 10 | -17.52 | 31.71 | |
| 317 | 15 | -52.43 | 14.7 | 9 | -16.21 | 24.23 | |
| 320 | 12 | -42.48 | 23.1 | 9 | -16.21 | 33.23 | |
| 平均 | 10.89 | -32.6 | 37.69 | 7.97 | -12.56 | 38.51 | |
新窗口打开
为了进一步判断利用优异等位基因进行抗病性强弱分析的可行性,将125份材料的优异等位基因数目、优异等位基因效应值之和与材料的黄萎病相对病指进行相关分析,结果表明,优异等位基因数目、效应值之和与材料的黄萎病相对病指间存在显著的相关性,其中,材料的优异等位基因效应值之和与材料的黄萎病病指间存在极显著的正相关,病圃鉴定与温室鉴定条件下的相关系数分别为0.616和0.566;同时,材料的优异等位基因数目与材料的黄萎病相对病指间存在极显著的负相关,相关系数分别为-0.618和-0.535（表3）,说明材料所含有的优异等位基因数目和优异等位基因效应值之和均可以用来判断材料的抗病性。
Table 3
表3
表3优异等位基因数目、效应值之和与黄萎病病指的相关分析
Table 3The correlation coefficients of number of elite alleles and sum of effective values of elite alleles with relative disease index of Verticillium wilt resistance
| 类别 Category | 优异等位基因数目 Number of elite alleles | 效应值之和 Sum of effective value | 相对病情指数 Relative disease index | |||||
|---|---|---|---|---|---|---|---|---|
| 病圃 Disease nursery | 温室 Greenhouse | 病圃 Disease nursery | 温室 Greenhouse | 病圃 Disease nursery | 温室 Greenhouse | |||
| 优势等位基因数目 Number of elite alleles | 1 | 1 | ||||||
| 效应值之和 Sum of effective value | -0.963** | -0.91** | 1 | 1 | ||||
| 相对病情指数 Relative disease index | -0.618** | -0.535** | 0.616** | 0.566** | 1 | 1 | ||
新窗口打开
3 讨论
目前,对于棉花黄萎病抗性的鉴定,主要有田间病圃鉴定、营养钵苗期鉴定和黄萎病菌毒素法鉴定[4]。其中,田间病圃鉴定能真实反映鉴定材料抗病性的强弱,是目前最常用的方法,然而该方法存在大田环境下的地力条件难以控制、需要每年补接菌量及鉴定周期长等限制因素;营养钵苗期鉴定通常在温室可控条件下进行,具有快速、准确、易操作等优点,目前被广泛使用。然而,由于黄萎病菌致病力的分化,常常导致不同批次、不同年度间同一品种的病情指数变异较大。近年来,通过设置合适的感病对照品种,利用相对病情指数来划分棉花品种的抗病性证明是一种减少环境变异的有效措施[23,26]。本研究中利用田间病圃鉴定和温室营养钵苗期鉴定研究了125份棉花优异种质的抗病性,田间病圃鉴定反映的是整个生育期的发病情况,而温室营养钵苗期鉴定主要反映苗期的发病情况,我们的结果表明,在对于耐病材料和感病的材料的鉴定方面,温室鉴定和田间病圃鉴定具有一定程度的一致性;尽管对于抗病材料的鉴定,由于主效抗病基因的表达具有时空性,导致少部分在温室条件下鉴定为抗病的材料在病圃鉴定条件下表现为耐病,或者反之,然而通过方差分析表明不同鉴定条件下获得的相对病情指数间不存在显著差异,即环境效应不显著,进一步说明了使用相对抗性指数进行抗病性鉴定可消除许多不确定因素对鉴定结果的影响。目前,国内外****通过构建双亲间抗病性存在显著差异的遗传分离群体,所定位的与棉花不同生长阶段的抗黄萎病性状相关的分子标记或QTL已超过了100个[5-17],这些分子标记或QTL为开展标记辅助选择抗病育种奠定了基础。本研究的分子标记是通过对大量抗病性存在显著差异的陆地棉优异种质进行关联分析获得的,通过与上述QTL结果进行比较,发现在本研究所获得的与抗黄萎病性状显著关联的SSR标记位点中,有8个SSR标记位点位于上述抗黄萎病QTL区间内,其中,NAU980-2、NAU5064-2、JESPR0001-1 3个优势等位基因位点（表1）均位于蒋锋等[15]定位的棉花抗黄萎病QTL qV-VD8M-D9-1基因组区域内（Dt_chr9染色体的29 224 283—55 190 112 bp）。这些结果证明了本研究与前人研究的一致性。祁伟彦等[3]利用和黄萎病抗性紧密连锁的SSR 标记对棉花黄萎病抗性进行了辅助筛选,证明了基于人工病圃筛选和分子标记辅助育种相结合是选育棉花抗黄萎病材料可行、高效的育种方法。同上述基于遗传分离群体的分子标记定位相比,本研究检测到了更多的优异种质材料中存在的抗黄萎病优异等位基因变异,有助于实现基于不同优势等位基因的分子标记辅助选择育种。
棉花抗黄萎病鉴定的一个主要限制因素是病原菌小种的易变性导致品种抗病性的变异。培育具有广谱抗病性的品种是克服病原菌小种变异的一种有效手段。而广谱抗病性产生的分子基础是多个抗黄萎病相关基因的累加。本研究发现种质材料中含有的优异等位基因数目与材料的黄萎病相对病指间存在极显著的负相关,说明了抗黄萎病优异等位基因的累加具有明显提高抗病性的作用。
近来,王铭等[27]通过比较多种鉴定棉花抗黄萎病的方法,提出一种使用黄萎病菌毒素进行联合鉴定的方法,以替代传统病圃鉴定法。尽管基于该方法的鉴定结果与田间病圃法鉴定结果的相关系数为0.94（P<0.01）,但是该方法仍然没有有效克服病原菌小种的易变性对鉴定结果的影响,也没有实现对抗病基因型的直接鉴定。本研究通过将棉花抗黄萎病分子标记进行等位基因分解,将抗病性分子标记转化为优异等位基因,并通过计算优异等位基因的效应值实现了分子标记的数量化,从而将抗病性表型与基因型联系起来。从理论上讲,优异等位基因的效应值在一定程度上反映了基因对表型的贡献,多个优异等位基因的聚合必将对表型产生积极的影响。我们的研究表明,不同材料中含有的优异等位基因数目、优异等位基因效应值之和与材料的黄萎病相对病指间存在显著的相关关系,其中,材料的优异等位基因效应值之和与材料的黄萎病相对病指间存在极显著的正相关,材料的优异等位基因数目与材料的黄萎病相对病指间存在极显著的负相关,在一定程度上证明了多个优异等位基因的聚合有提高材料抗病性的作用,说明材料所含有的优异等位基因数目和优异等位基因效应值之和在一定程度上可以反映材料抗病性的强弱,从而实现对抗病性的分子鉴定。
研究表明,陆地棉黄萎病抗性属数量性状遗传[28],因此容易受到环境条件变化的影响。由于田间病圃鉴定反映的是棉花整个生育期的发病情况,而温室营养钵苗期鉴定主要反映苗期的发病情况,再加上抗病相关基因的表达具有时空性,从而导致田间病圃鉴定条件下获得的优异等位基因位点与温室营养钵苗期接种鉴定条件下获得的优异等位基因位点存在差异。尽管如此,本研究中我们同时获得了在田间黄萎病病圃鉴定和温室接种鉴定2个环境中均显著关联的抗黄萎病优异等位基因及其效应值,可以推测,通过采用新的研究方法,发掘大量的多环境稳定的抗黄萎病优势等位基因,建立基于抗病优异等位基因效应值之和为标准的抗黄萎病实时检测技术,可以实现对抗病性的快速鉴定和对抗病基因型的直接选择,有效解决当前抗病性鉴定周期长、抗病性选择效率低的技术难题。
4 结论
棉花种质材料中抗黄萎病优异等位基因的累加具有明显提高抗病性的作用,材料所含有的优异等位基因数目和优异等位基因效应值之和在一定程度上可以反映材料抗病性的强弱,从而实现对抗病性的分子鉴定。Table 1
附表1
附表1本研究所选用的棉花材料
Table 1 Cotton lines used in this study
| 材料编号 Line ID | 材料名称 Line name | 系谱/来源 Pedigree | 材料编号 Line ID | 材料名称 Line name | 系谱/来源 Pedigree | |
|---|---|---|---|---|---|---|
| 5 | 常德184 Changde184 | N | 158 | Coker139 | 美国 USA | |
| 8 | 沪749513 Hu749513 | 巴州农业科学研究所 Bazhou Agricultural Research Institute | 161 | GP83 | 美国 USA | |
| 9 | 淮棉4 Huaimian4 | 安徽萧县 Xiaoxian, Anhui | 163 | 中89-1 Zhong89-1 | N | |
| 12 | 库车T94-6 KucheT94-6 | 新疆Xinjiang | 164 | Miscot78-27 | 美国 USA | |
| 13 | 沙帽长桃 Shamaochangtao | 湖北武汉 Wuhan, Hubei | 166 | 美D#1 Mei D#1 | 美国 USA | |
| 15 | 苏7036远缘 Su7036 | 86-1/异常棉//陆地棉 86-1/G.anomalum/G.hirsutum L. | 168 | AC239 | Line A/C6-5 | |
| 16 | 苏农6号 Sunong6 | N | 169 | M-8124-1159 | 美国 USA | |
| 22 | 邯8959 Han8959 | PD6186选 PD6186 | 171 | 美8123 Mei8123 | 美国 USA | |
| 24 | 冀91-31 | 邯88-S92 | 172 | 中资04184 | Arcot-1/SGK9708 | |
| Ji91-31 | Han88-S92 | Zhongzi04184 | ||||
| 26 | 荆55263 | 3247/临汾4032+中9 | 173 | Acala SJ-4 | C6TE/NMB3080 | |
| Jing55263 | 3247/Linfen4032+Zhong9 | |||||
| 27 | 奎屯80-2056W(新陆早3号) Kuitun80-2056W | 66-241/美国爱字无毒棉//荆州4588 66-241/Acala non-toxic cotton// Jingzhou4588 | 174 | Arcot-1 | 美国 USA | |
| 179 | LAPAR45 | 美国 USA | ||||
| 32 | 苏联8908 Soviet Union8908 | 前苏联 Former Soviet Union | 180 | PAR-51 | 美国 USA | |
| 44 | 山农6号 Shannong6 | 山东泰安 Taian, Shandong | 182 | 非洲E40 AfricaE40 | 非洲 Africa | |
| 45 | 石远638 Shiyuan638 | 陆地棉/比克氏棉 G.hirsutum L./G.bickii Prokh. | 185 | 中Arc-185 ZhongArc-185 | Arcot-1原子能诱变 Induced mutagenesis of Arcot-1 | |
| 46 | 苏联21 Soviet Union 21 | 前苏联 Former Soviet Union | 186 | L142-9 | L142-9 | |
| 48 | 苏徐138(早熟) Suxu138 | N | 189 | RNX189 | 光子Nn/Xu142的014HF GuangziNn/014HF of Xu142 | |
| 51 | 中ARR40682 ZhongARR40682 | Arcot-1原子能诱变 Induced mutagenesis of Arcot-1 | 191 | 澳152 Ao152 | 澳大利亚 Australia | |
| 54 | S-050031 | 美棉723446（岱46）/早熟品系202 Meimian723446（Deltapine46）/Zaoshu202 | 199 | 苏TKH-1 SuTKH-1 | 江苏南京 Nangjing, Jiangsu | |
| 55 | 冀邯2号选系 Jihan2 | 河北邯郸 Handan, Hebei | 201 | 莘棉5号 Shenmian5 | 斯517/8763依 Si517/8763Yi | |
| 59 | 农救 Nongjiu | 美国 USA | 207 | UA887 | 美国USA | |
| 62 | 肖县大铃 Xiaoxian big ball | N | 209 | DP33B | 岱字棉90\岱字棉50 转Bt基因 Deltapine90\Deltapine50.Transgene(Bt) | |
| 63 | Acala(大铃)B Acala(big boll)B | 美国 USA | 211 | 豫棉9号(豫早1109) Yumian9(Yuzao1109) | 中抗5/中棉所10号//中棉所14 Zhongkang5/Zhongmian10//Zhongmian14 | |
| 续表1 Continued table 1 | ||||||
| 材料编号 Line ID | 材料名称 Line name | 系谱/来源 Pedigree | 材料编号 Line ID | 材料名称 Line name | 系谱/来源 Pedigree | |
| 66 | 大铃棉 Big boll cotton | 新疆 Xinjiang | 213 | ST474 | 美国 USA | |
| 67 | 大铃棉69 Big boll cotton 69 | 新疆Xinjiang | 216 | 贝尔斯诺 Beiersinuo | 美国 USA | |
| 69 | 启丰大铃 Qifeng big boll | N | 218 | 泗棉2号 Simian2 | 泗阳437/墨西哥910 Siyang437/Mexico910 | |
| 70 | 新大铃 New big boll | N | 222 | 中9708-222 Zhong9708-222 | SGK9708原子能诱变 Induced mutagenesis of SGK9708 | |
| 73 | 中ARR40681 ZhongARR40681 | Arcot-1原子能诱变 Induced mutagenesis of Arcot-1 | 227 | 邯郸568 Handan568 | 河北邯郸 Handan, Hebei | |
| 76 | 中Arc-76 ZhongArc-76 | Arcot-1原子能诱变 Induced mutagenesis of Arcot-1 | 229 | 中资9103 Zhongzi9103 | Ct5R2Lc选 Ct5R2Lc | |
| 77 | 中AR683-77 ZhongAR683-77 | Arcot-1辐射诱变 Induced mutagenesis of Arcot-1 | 235 | 冀91-12 Ji91-12 | 辽1038选 Liao1038 | |
| 84 | GP70 | 美国 USA | 237 | 冀85-3 Ji85-3 | 海陆杂交85-5673 85-5673(G.hirsutum L.×G. barbadense L.) | |
| 89 | 鲁棉11号(3389) Lumian11(3389) | 鲁2/鲁1抗 Lu2/Lu1Kang | 238 | 冀91-32 Ji91-32 | 邯88-S92 Han88-S92 | |
| 90 | 豫284 Yu284 | N | 240 | FB20 | 巴基斯坦 Pakistan | |
| 91 | 中5913-2 Zhong5913-2 | N | 241 | 辽96-23-30 Liao96-23-30 | 辽宁 Liaoning | |
| 94 | 山农3号 Shannong3 | 108夫/5476依 108Fu/5476Yi | 243 | GK99-1 | GK12选系\82系 GK12 \82 | |
| 96 | 邢台79-11 Xingtai79-11 | 河北邢台 Xingtai, Hebei | 245 | DP410B | DP5409\531系 DP5409\531 | |
| 98 | 徐州半半棉 Xuzhoubanbanmian | 江苏经作所 Jiangsu Agricultural Research Institute | 253 | DP2156 | 美国 USA | |
| 99 | 盐城1115 Yancheng1115 | 江苏盐城 Yancheng, Jiangsu | 255 | r4136 | N | |
| 103 | 中2108 Zhong2108 | 乌干达4号/冀棉1号-6429//86-1 Uganda4/Jimian1-6429//86-1 | 261 | RT白絮 RT white fiber | TM-1/G.richmondi白絮 TM-1/G.richmondi | |
| 106 | CZA(70)33 | 美国 USA | 262 | MSCO-11 | Miscot 的选系 Selection from Miscot | |
| 107 | DES926 | 美国 USA | 263 | upland | 美国 USA | |
| 108 | GP93 | 美国 USA | 264 | 鲁458 Lu458 | N | |
| 111 | N74-250 | 美国 USA | 265 | 豫棉19 Yumian19 | 辐射矮源多亲复交 Composite cross | |
| 113 | Qik | 美国 USA | 266 | 望江长绒棉 Wangjiang long lint cotton | 安徽望江 Wangjiang, Anhui | |
| 116 | GK22 | N | 269 | 徐州261 Xuzhou261 | 江苏徐州 Xuzhou, Jiangsu | |
| 117 | 邯109 Han109 | N | 270 | 邯郸333 Handan333 | 陆海杂交 G.hirsutum L.×G. barbadense L. | |
| 121 | JCG59 | N | 276 | 郑州长绒棉 Zhengzhou long lint cotton | N | |
| 续表1 Continued table 1 | ||||||
| 材料编号 Line ID | 材料名称 Line name | 系谱/来源 Pedigree | 材料编号 Line ID | 材料名称 Line name | 系谱/来源 Pedigree | |
| 123 | 邯8944 Han8944 | N | 277 | 苏棉2号 Sumian2 | N | |
| 124 | 抗黄萎164 Verticillium wilt-resistance164 | {[Acala5658×陕5254]×4067}×中棉所10号}×冀8 {[Acala5658×Shan5254]×4067}×Zhongmian10}×Ji8 | 280 | PD6186 | SC-1/PD8619(SC-1//CKR421/PD2164) | |
| 281 | 高衣棉 Gaoyimian | N | ||||
| 127 | 陕3184 Shan3184 | N | 282 | MSCO-12 | Miscot 的选系 Selection from Miscot | |
| 128 | 陕960329-2远3 Shan960329-2yuan3 | N | 287 | 新研96-48 Xinyan96-48 | 新七8704-7/新七90-1,钴60辐射处理 Xinqi8704-7/Xinqi90-1,Induced mutagenesis and trangene（Bt+CPTI） | |
| 130 | 中870203 Zhong870203 | N | 288 | 豫棉2号 Yumian2 | 河南郑州 Zhengzhou, Henan | |
| 132 | 中植BD13 ZhongzhiBD13 | 75-23/辽7-10 75-23/Liao7-10 | 291 | LineF | N | |
| 135 | 中资640 Zhongzi640 | 鲁棉1号选系 Selection from Lumian1 | 295 | 中棉所12 Zhongmian12 | 乌干达4号/邢台6871 Uganda4/Xingtai6871 | |
| 136 | 辽4835 Liao4835 | N | 299 | 中棉所50 Zhongmian50 | H109/中662 H109/Zhong662 | |
| 137 | 辽4853 Liao4853 | N | 307 | 中507145 Zhong507145 | 中221 Zhong221 | |
| 138 | 辽棉16 Liaomian16 | 辽4086\辽棉6号+6496+4229）混粉 Liao4086\Liaomian6+6496+4229.Mixed pollination | 308 | 中Arc-308 ZhongArc-308 | Arcot-1原子能诱变 Induced mutagenesis of Arcot-1 | |
| 139 | 辽棉17 Liaomian17 | 辽4084\辽4258 Liao4084\Liao4258 | 311 | 中AR-RNX311 ZhongAR-RNX311 | 中ARR40683-4/RILNnXu0082 F6 ZhongARR40683-4/RILNnXu0082.F6 | |
| 140 | 辽棉18 Liaomian18 | 辽87-7109\辽1038-6 Liao87-7109\Liao1038-6 | 312 | 中G5 ZhongG5 | {[（黑山棉1号×简中1号）F1×斯特提棉]×中8}BC×中7 {[（Heishanmian1×Jianzhong1）F1×G．sturtianum]×Zhong8}×Zhong7 | |
| 141 | 永济1号 Yongji1 | N | ||||
| 143 | 中远9111 Zhongyuan9111 | 381//[亚洲棉/斯特提棉//陆地棉]F8 381//[G.arboreum L./G．sturtianum//G.hirsutum L]F8 | 317 | 99633 | Sicala996/DPL33B | |
| 144 | 中远9114 Zhongyuan9114 | 中10/墨西哥棉 Zhong10/Mexico cotton | 320 | Ari-320 | Ari971激光诱变 Induced mutagenesis of Ari971 | |
| 150 | 中远9112 Zhongyuan9112 | 亚洲棉/斯特提棉//中10/中10 G.arboreum.L./G．sturtianum//Zhong10/Zhong10 | ||||
| 151 | 豫棉2067 Yumian2067 | 中棉所12\豫植177 Zhongmian12\Yuzhi177 | ||||
新窗口打开
The authors have declared that no competing interests exist.
