收稿日期：2021-04-30" " "第一作者簡介：鐘鎮(zhèn)濤（1995―），男，碩士研究生，21816016@zju.edu.cn。 * 通信作者：shjzhu@zju.edu.cn

基金項目：國家重點研發(fā)計劃（2016YFD0101400）；國家棉花產(chǎn)業(yè)技術(shù)體系（CARS-18-25）；江蘇省現(xiàn)代作物生產(chǎn)協(xié)同創(chuàng)新中心

摘要：【目的】棉花是1種較耐重金屬鎘的作物，培育與種植耐鎘脅迫棉花品種對于鎘污染土壤的修復(fù)具有重要意義。本文旨在研究陸地棉耐鎘性的遺傳并對相關(guān)數(shù)量性狀位點 （quantitative trait locus， QTL）進(jìn)行定位分析?！痉椒ā恳躁懙孛薹N質(zhì)HS46和MARCABUCAG8US-1-88為親本構(gòu)建的188個重組自交系（recombinant inbred line， RIL）為材料，調(diào)查RIL群體在鎘脅迫條件下的發(fā)芽率、發(fā)芽勢，以及苗期葉片的葉綠素含量、株高和干/鮮物質(zhì)質(zhì)量等表型數(shù)據(jù)，并計算各家系的鎘脅迫系數(shù)；采用完備區(qū)間作圖法對RIL群體的各性狀的鎘脅迫系數(shù)及鎘脅迫綜合系數(shù)進(jìn)行QTL定位和遺傳分析?！窘Y(jié)果】MARCABUCAG8US-1-88為耐鎘材料，HS46則為鎘敏感材料，188個RIL家系的耐鎘性相關(guān)性狀基本服從正態(tài)分布，包含強(qiáng)耐鎘型材料23份、耐鎘型材料49份、中等耐鎘型材料67份、鎘敏感型材料49份。連鎖分析共鑒定到28個與鎘脅迫系數(shù)相關(guān)的QTL，分布于14條染色體上，單個位點解釋1.58%～8.41%的表型變異，其中根鮮物質(zhì)質(zhì)量鎘脅迫系數(shù)位點qRFW-13-1在2個環(huán)境下均被檢測到，分別解釋7.66%和7.71%的表型變異?！窘Y(jié)論】棉花耐鎘性為多基因控制的數(shù)量遺傳，受環(huán)境影響較大，定位到的耐鎘QTL qRFW-13-1為穩(wěn)定QTL，該結(jié)果對于棉花耐鎘分子育種具有一定的應(yīng)用價值。

關(guān)鍵詞：棉花；鎘脅迫；QTL定位

QTL mapping and genetic analysis of cadmium stress tolerance in upland cotton

Zhong Zhentao， Zhao Tianlun， Chen Jinhong， Zhu Shuijin*

（College of Agriculture and Biotechnology/Zhejiang Key Laboratory of Crop Germplasm， Zhejiang University， Hangzhou 310058， China）

Abstract： [Objective] Cotton is a relatively tolerant crop to cadmium （Cd） stress. It is of great significance to breeding and extending the Cd tolerant cotton cultivars for remediation of Cd contaminated soil. This study aims to analysis the inheritance of cotton tolerance to Cd stress and to map the quantitative trait locus （QTL） related to Cd stress. [Methods] In this study， 188 recombinant inbred lines （RIL） called HM188 derived from an intraspecific hybrid between HS46 and MARCABUCAG8US-

1-88 were used to study the phenotypic data of germination rate， germination potential， chlorophyll content， plant height， dry and fresh weight at seedling stage under Cd stress， and the Cd stress coefficient of each material was calculated as well. QTL mapping was conducted for each trait and Cd stress coefficient of HM188 population by using complete interval mapping method. [Results] According to the performance and Cd stress coefficient of each trait under Cd stress， MARCABUCAG8US-

1-88 was the Cd tolerant parent， while HS46 was a sensitive one. The Cd tolerance of 188 RIL was basically normal distribution， including 23 strong Cd tolerant lines， 49 Cd tolerant line， 67 medium tolerant lines， and 49 sensitive lines. A total of 28 QTL for Cd stress coefficient was identified， which were distributed on 14 chromosomes with explaining of variation of 1.58%～8.41%. Among them， qRFW-13-1 related to root fresh weight was detected in both environments， which explained variation of 7.66% and 7.71%， respectively. [Conclusion] Cd tolerance in upland cotton is a quantitative trait controlled by multiple genes， which was greatly affected by the environment. However， qRFW-13-1 for Cd tolerance may be a stable QTL， which has some application values in molecular breeding for Cd tolerant cotton.

Keywords： cotton; cadmium stress; QTL mapping

隨著全球人口劇增，工農(nóng)業(yè)生產(chǎn)以及日常生活廢棄物對土壤產(chǎn)生了嚴(yán)重的污染。2014年我國原國土資源部和環(huán)境保護(hù)部聯(lián)合發(fā)布的《全國土壤污染狀況調(diào)查公報》顯示，我國土壤環(huán)境狀況總體不樂觀，有19.4%的受調(diào)查耕地被污染，影響耕地環(huán)境質(zhì)量的主要污染物是重金屬，其中鎘（cadmium，Cd）為首要污染物，在超過40%的受污染耕地中發(fā)現(xiàn)重金屬Cd超標(biāo)[1]。高濃度的Cd會對土壤生物產(chǎn)生毒害作用，并很容易轉(zhuǎn)移到植物中，最終進(jìn)入食物鏈[2]。攝入過多Cd會嚴(yán)重?fù)p害人體健康，導(dǎo)致癌癥、骨質(zhì)疏松和其他疾病[3]。因此，對Cd污染土壤的修復(fù)刻不容緩。目前針對重金屬污染土壤的修復(fù)技術(shù)包括物理修復(fù)、化學(xué)修復(fù)、生物修復(fù)和農(nóng)業(yè)生態(tài)修復(fù)技術(shù)等，生物修復(fù)法中的植物修復(fù)技術(shù)因成本低，可修復(fù)大面積污染土壤且具有一定綠化美化作用而備受關(guān)注。目前全世界已發(fā)現(xiàn)500多種富集重金屬的植物，其中部分植物表現(xiàn)出對Cd的選擇性吸收，如菠菜、煙草、向日葵等[4]。孫紅艷[5]在大麥單雙倍體（double haploid， DH）群體鑒定到11個與Cd含量相關(guān)的數(shù)量性狀位點（quantitative trait locus， QTL），并研究了ZIP3、ZIP5、ZIP7和ZIP8這4個基因?qū)Υ篼淐d轉(zhuǎn)運與積累的調(diào)控與影響。潘晨陽等[6]在水稻重組自交系（recombinant inbred line， RIL）群體中鑒定到了7個Cd積累相關(guān)QTL，并篩選出LOCOs02g50240、LOCOs02g52780、LOCOs09g31200、LOCOs09g35030和LOCOs09g-

37949為耐Cd基因。

棉花是我國重要的經(jīng)濟(jì)作物。研究表明中低濃度Cd處理對棉花生長發(fā)育影響較小，棉花對土壤Cd具有良好的富集能力，棉纖維為非食用或飼用產(chǎn)品且Cd吸附量較低，棉花是理想的Cd污染耕地的修復(fù)作物之一[7-11]。植棉修復(fù)重金屬污染耕地具有經(jīng)濟(jì)效益高、對土壤和環(huán)境友好、適用于重金屬復(fù)合污染等優(yōu)勢[10]。然而，土壤中Cd質(zhì)量濃度≥20 mg·kg－1時會嚴(yán)重抑制棉花種子萌發(fā)和幼苗生長發(fā)育，顯著降低單株鈴數(shù)、鈴重、單鈴種子數(shù)、籽棉產(chǎn)量和皮棉產(chǎn)量[12]。目前對于棉花耐Cd性和Cd積累的研究較少，尚未見與棉花耐Cd性和Cd積累相關(guān)的QTL報道，對于棉花耐Cd的遺傳表現(xiàn)了解甚少，強(qiáng)耐Cd、高Cd富集棉花品種不足以滿足Cd污染耕地修復(fù)的需要。本研究以陸地棉（Gossypium hirsutum）RIL群體為材料，研究Cd脅迫條件下棉花發(fā)芽和苗期生長情況，分析棉花苗期耐Cd性的遺傳特點，并對棉花耐Cd性進(jìn)行QTL定位，為耐Cd品種選育提供遺傳資源和材料資源。

1 材料與方法

1.1 試驗材料

以綜合性狀優(yōu)良的陸地棉品種HS46和對多種逆境都有抗性的陸地棉種質(zhì)系MARCABUCAG8US-1-88（簡稱MAR），以及兩者雜交獲得的188個RIL群體（簡稱HM188）為供試材料[13]。RIL通過單籽粒傳法自交8代構(gòu)建完成，自交繁殖保存。試驗用種為2018年冬種于海南三亞收獲的南繁自交種子，由浙江大學(xué)農(nóng)學(xué)院棉花基因工程實驗室提供。

1.2 試驗方法

發(fā)芽試驗：各材料脫絨種子曬干。發(fā)芽前用75%的酒精消毒，用無菌水沖洗干凈后浸種12 h。挑選50粒種子均勻平鋪到墊有雙層濾紙的培養(yǎng)皿中，每皿加入600 μmol·L－1 CdCl2溶液20 mL進(jìn)行Cd脅迫處理，以蒸餾水作對照，4次重復(fù)，放置于28 ℃、光暗周期為14 h/10 h的植物培養(yǎng)室中發(fā)芽，記錄每天發(fā)芽種子數(shù)。

苗期試驗：試驗設(shè)2個環(huán)境，分別于2019年6－7月在浙江省杭州市浙江大學(xué)華家池校區(qū)實驗站和2020年7－8月在浙江大學(xué)紫金港校區(qū)控溫光生長室進(jìn)行。采用裂區(qū)隨機(jī)試驗設(shè)計，主區(qū)為Cd脅迫處理（Cd），以非脅迫處理（CK）為對照；副區(qū)為RIL群體和雙親，對照為中棉所49。按試驗設(shè)計，硫酸脫絨的飽滿種子分別播于50（10×5）孔穴盤內(nèi)，每個株系種植1行（10孔），每孔播2粒飽滿的脫絨種子，加營養(yǎng)土至與孔口平。每穴盤中間1行播中棉所49，用以校正穴盤之間的試驗誤差。待幼苗1片真葉完全展開后，將處理組（Cd）用600 μmol·L－1 CdCl2溶液澆灌至飽和，保持10 d；對照組（CK）用水澆灌至飽和，保持10 d。處理10 d后，各處理每個株系隨機(jī)取10株幼苗，分別測定其葉片葉綠素含量（以soil and plant analyzer development即SPAD值反映葉片葉綠素相對含量）、株高（plant height， PH）、根鮮物質(zhì)質(zhì)量（root fresh weight， RFW）、根干物質(zhì)質(zhì)量（root dry weight， RDW）、地上部鮮物質(zhì)質(zhì)量（shoot fresh weight， SFW）和地上部干物質(zhì)質(zhì)量（shoot dry weight， SDW），取平均值，為1次重復(fù)。試驗設(shè)2次重復(fù)。2020年利用葉綠素分析儀（SPAD-502Plus，日本）測定幼苗1、2、3片真葉的SPAD值，取平均值為該株材料的SPAD值。

1.3 性狀統(tǒng)計

利用Excel 2014和SPSS 20.0軟件對親本及RIL群體的表型數(shù)據(jù)進(jìn)行統(tǒng)計分析和聚類分析，運用 Kolmogorov-Smirnov（D檢驗）和 Shapiro-

Wilk（W檢驗）進(jìn)行正態(tài)分布檢驗。參照朱協(xié)飛等[14]的方法計算Cd脅迫系數(shù)：

Cd脅迫系數(shù)＝×100%

參照岳松青等[15]的方法，計算各性狀Cd脅迫系數(shù)的隸屬函數(shù)值Xμ：

Xμ＝；

式中，X為供試材料某一性狀的Cd脅迫系數(shù)，Xmax為供試材料中某一性狀Cd脅迫系數(shù)的最大值，Xmin為供試材料中某一性狀Cd脅迫系數(shù)的最小值。將每個材料各性狀的隸屬函數(shù)值相加后取平均值，作為該品種的Cd脅迫綜合系數(shù)，綜合系數(shù)越小，該材料的耐Cd性越強(qiáng)。

1.4 連鎖圖譜與QTL定位

本研究所用RIL群體的連鎖圖譜由Li等[16]構(gòu)建。利用QTL IciMapping 4.2軟件中BIP模塊，采用完備區(qū)間作圖法（inclusive composite interval mapping，ICIM）分別對葉片SPAD值、株高、根干/鮮物質(zhì)質(zhì)量和地上部干/鮮物質(zhì)質(zhì)量進(jìn)行QTL定位，參數(shù)設(shè)置：步長為1 厘摩（centiMorgan， cM），PIN為0.000 1，以對數(shù)優(yōu)勢比（logarithm of odd score，LOD）≥2.5作為閾值判斷QTL存在。QTL命名規(guī)則[17]：q＋性狀英文名稱縮寫＋染色體序號＋QTL編號。

2 結(jié)果與分析

2.1 HM188群體及雙親的發(fā)芽情況及苗期耐Cd表現(xiàn)

2個環(huán)境下的HM188及其親本在Cd脅迫條件下的發(fā)芽情況、苗期生長情況及對應(yīng)的Cd脅迫系數(shù)統(tǒng)計結(jié)果見表1～3。

從表1來看，在對照和Cd脅迫條件下，親本MAR的發(fā)芽率、發(fā)芽勢均顯著高于HS46，且MAR發(fā)芽率和發(fā)芽勢的Cd脅迫系數(shù)顯著低于HS46（表3），說明2個親本的耐Cd性存在顯著差異，MAR較HS46耐Cd。在Cd脅迫條件下，HM188群體的平均發(fā)芽勢顯著低于對照，發(fā)芽率與對照無顯著差異。在對照條件下，HM188的發(fā)芽率為57.0%～95.0%，表現(xiàn)出超親分離現(xiàn)象。在Cd脅迫下，HM188的發(fā)芽率為37.6%～94.0%，變異幅度更大。在Cd脅迫條件下，HM188群體的發(fā)芽勢為1.2%～59.9%，表現(xiàn)出超親分離現(xiàn)象。這些結(jié)果表明，在Cd脅迫下棉花種子能發(fā)芽，但種子活力顯著下降。

脅迫系數(shù)結(jié)果（表3）表明，HM188群體的發(fā)芽勢在Cd脅迫條件下平均下降了41.18%，發(fā)芽勢的Cd脅迫系數(shù)（－18.80%～97.50%）呈連續(xù)分布，偏度與峰度絕對值均小于1，基本符合正態(tài)分布，且表現(xiàn)出超親分離。說明Cd脅迫下的發(fā)芽勢可作為棉花耐Cd性指標(biāo)之一。

在正常和Cd脅迫處理條件下，親本MAR和HS46間的SPAD值均無顯著差異，但在Cd脅迫條件下MAR和HS46的SPAD值較對照分別下降了17.5%和28.3%。HM188群體的SPAD平均值在Cd脅迫條件下顯著低于無脅迫的對照。脅迫系數(shù)結(jié)果表明，親本MAR葉片SPAD值的Cd脅迫系數(shù)顯著低于HS46，HM188群體的葉片SPAD值在Cd脅迫條件下平均下降了23.68%，SPAD值的Cd脅迫系數(shù)呈連續(xù)分布（18.52%～29.54%），且表現(xiàn)出超親分離，Cd脅迫系數(shù)分布的偏度與峰度絕對值均小于1，符合正態(tài)分布。盡管只有1年的測定數(shù)據(jù)，但也能夠說明Cd脅迫影響植株葉片的SPAD值并影響光合作用。

在正常條件下，親本MAR的株高顯著低于HS46，而Cd脅迫條件下MAR的株高顯著高于HS46（表2），盡管2個環(huán)境的結(jié)果不盡一致，但趨勢相同（結(jié)果未顯示）。親本MAR的株高Cd脅迫系數(shù)顯著低于HS46（表3）。Cd脅迫條件下HM188群體的株高為8.1～11.0 cm，平均為9.8 cm，顯著低于CK，平均下降了30.50%。親本MAR和HS46的株高Cd脅迫系數(shù)分別為22.97%和37.23%；HM188的株高Cd脅迫系數(shù)呈連續(xù)分布，變幅為21.47%～40.97%，具有超親表現(xiàn)，Cd脅迫系數(shù)分布的偏度與峰度絕對值均小于1，符合正態(tài)分布。Cd脅迫下的株高可作為棉花耐Cd性指標(biāo)之一。

在正常條件下，親本MAR的根鮮物質(zhì)質(zhì)量和根干物質(zhì)質(zhì)量均顯著低于HS46，而在Cd脅迫條件下2個親本的根鮮物質(zhì)質(zhì)量和根干物質(zhì)質(zhì)量均無顯著差異（表2）。在Cd脅迫條件下，HM188群體的根鮮物質(zhì)質(zhì)量和根干物質(zhì)質(zhì)量分別為0.462～0.823 g和0.032～0.052 g，平均分別為0.635 g和0.037 g，比CK分別顯著降低29.99%和35.09%（表2）。從表3來看，親本MAR根鮮物質(zhì)質(zhì)量、根干物質(zhì)質(zhì)量的Cd脅迫系數(shù)顯著低于HS46；在Cd脅迫條件下HM188的根鮮物質(zhì)質(zhì)量和干物質(zhì)質(zhì)量的Cd脅迫系數(shù)呈連續(xù)分布，分別為4.88%～44.44%和12.02%～49.15%，均表現(xiàn)出超親現(xiàn)象，偏度與峰度絕對值均小于1，符合正態(tài)分布。Cd脅迫下的根鮮物質(zhì)質(zhì)量和干物質(zhì)質(zhì)量可作為棉花耐Cd性指標(biāo)之一。

CK條件下，親本MAR的地上部鮮物質(zhì)質(zhì)量和干物質(zhì)質(zhì)量均顯著低于HS46，在Cd脅迫條件下，2個親本的地上部鮮物質(zhì)質(zhì)量和地上部干物質(zhì)質(zhì)量均無顯著差異（表2）。Cd脅迫條件下，HM188群體的地上部鮮物質(zhì)質(zhì)量和干物質(zhì)質(zhì)量變幅分別為1.421～2.097 g和0.161～0.245 g，平均分別為1.766 g和0.203 g，比CK顯著降低32.08%和36.36%（表2）。脅迫系數(shù)（表3）結(jié)果表明，親本MAR地上部鮮物質(zhì)質(zhì)量和干物質(zhì)質(zhì)量的Cd脅迫系數(shù)均顯著低于HS46，地上部鮮物質(zhì)質(zhì)量和干物質(zhì)質(zhì)量的Cd脅迫系數(shù)呈連續(xù)分布，并有超親表現(xiàn)，分布的偏度與峰度絕對值均小于1，符合正態(tài)分布。Cd脅迫下的地上部鮮物質(zhì)質(zhì)量和干物質(zhì)質(zhì)量可作為棉花耐Cd性指標(biāo)之一。

Cd脅迫條件下HM188群體的發(fā)芽勢、葉片SPAD值、株高、根干/鮮物質(zhì)質(zhì)量、地上部干/鮮物質(zhì)質(zhì)量比CK顯著降低。親本MAR的Cd脅迫綜合系數(shù)顯著低于HS46，表現(xiàn)出良好的耐Cd性。

2.2 HM188群體耐Cd性的聚類分析和遺傳力分析

采用最長距離法對各材料的Cd脅迫綜合系數(shù)進(jìn)行聚類，結(jié)果見圖1。聚類分析表明，HM188群體和雙親被劃分為4個類群：Ⅰ類為強(qiáng)耐Cd型材料，包含親本MAR、RIL-6、RIL-8和RIL-14等24份材料，占供試材料的12.6%；Ⅱ類為耐Cd型材料，包含RIL-5、RIL-7和RIL-13等49份材料，占供試材料的25.8%；III類為中等耐Cd型材料，包含RIL-1、RIL-2和RIL-4等67份材料，占供試材料的35.3%；IV類為Cd敏感型材料，包含親本HS46、RIL-9、RIL-10和RIL-32等50份材料，占供試材料的26.3%。

基于2019年和2020年的Cd脅迫系數(shù)，計算株高、根鮮物質(zhì)質(zhì)量、根干物質(zhì)質(zhì)量、地上部鮮物質(zhì)質(zhì)量和地上部干物質(zhì)質(zhì)量Cd脅迫系數(shù)的廣義遺傳率（H2），分別為19.709%、36.294%、21.909%、39.818%和37.048 %（表4），其中地上部鮮物質(zhì)質(zhì)量Cd脅迫系數(shù)的廣義遺傳率相對最高，地上部干物質(zhì)質(zhì)量的次之，株高的最低?？梢婈懙孛弈虲d性的遺傳率較低，受環(huán)境影響較大，遺傳改良較為困難。

2.3 Cd脅迫系數(shù)的QTL定位

HM188群體的發(fā)芽率、發(fā)芽勢、葉片SPAD值、株高、根干/鮮物質(zhì)質(zhì)量、地上部干/鮮物質(zhì)質(zhì)量這8個性狀的Cd脅迫系數(shù)及Cd脅迫綜合系數(shù)均呈連續(xù)性分布，偏度與峰度絕對值均小于1，基本服從正態(tài)分布，適合用于QTL定位（表3）。

通過QTL定位，在2個環(huán)境下共檢測到28個與這些性狀的Cd脅迫系數(shù)相關(guān)的QTL，5個與Cd脅迫綜合系數(shù)相關(guān)的QTL，分布在5、10、12、13、15、16、17、18、19、20、21、23、24和25號染色體上（表5、圖2）。

葉片SPAD值：共檢測到2個與葉片SPAD值 Cd脅迫系數(shù)相關(guān)的QTL，分布在19和21號染色體上，解釋2.64%～2.70%的表型變異，qSPAD-19-1的有利等位變異來自HS46，qSPAD-

21-1的有利等位變異來自MAR。

株高：共檢測到3個與株高Cd脅迫系數(shù)相關(guān)的QTL，分布在12、15和19號染色體上，解釋2.44%～7.26%的表型變異，qPH-15-1的有利等位變異來自HS46，qPH-12-1和qPH-19-1的有利等位變異來自MAR。

根鮮物質(zhì)質(zhì)量：共檢測到4個與根鮮物質(zhì)質(zhì)量Cd脅迫系數(shù)相關(guān)的QTL，5號染色體上2個、13和17號染色體上各1個，解釋6.62%～8.41%的表型變異。其中qRFW-13-1在2019和2020年2個環(huán)境下均被檢測到，分別解釋7.66%和7.71%的表型變異。

根干物質(zhì)質(zhì)量：共檢測到6個與根干物質(zhì)質(zhì)量Cd脅迫系數(shù)相關(guān)的QTL，5、13、17、19、23和25號染色體上各1個，解釋了1.58%～8.32%的表型變異。其中qRDW-5-1與根鮮物質(zhì)質(zhì)量Cd脅迫系數(shù)位點qRFW-5-1為同一位點，有利等位變異來自MAR，qRDW-13-1與根鮮物質(zhì)質(zhì)量Cd脅迫系數(shù)位點qRFW-13-1為同一位點，有利等位變異來自HS46，qRDW-17-1與根鮮物質(zhì)質(zhì)量Cd脅迫系數(shù)位點qRFW-17-1為同一位點，有利等位變異來自MAR。

地上部鮮物質(zhì)質(zhì)量：共檢測到3個與地上部鮮物質(zhì)質(zhì)量Cd脅迫系數(shù)相關(guān)的QTL，10、13和18號染色體上各1個，解釋1.92%～7.65%的表型變異，有利等位變異均來自MAR。

地上部干物質(zhì)質(zhì)量：共檢測到4個與地上部干物質(zhì)質(zhì)量Cd脅迫系數(shù)相關(guān)的QTL，16號染色體上2個，10和24號染色體上各1個，解釋4.93%～6.31%的表型變異，4個QTL的有利等位變異均來自MAR。其中qSDW-10-1和地上部鮮物質(zhì)質(zhì)量Cd脅迫系數(shù)位點qSFW-10-1為同一位點。

共檢測到1個與發(fā)芽率Cd脅迫系數(shù)相關(guān)的QTL，位于12號染色體，解釋2.47%的表型變異。共檢測到5個與發(fā)芽勢Cd脅迫系數(shù)相關(guān)的QTL，5、12、20、21和24號染色體上各1個，解釋2.07%～6.30%的表型變異。其中，qGP-12-1和發(fā)芽率Cd脅迫系數(shù)位點qGR-12-1、株高Cd脅迫系數(shù)位點qPH-12-1為同一位點，其有利等位變異來自HS46。

Cd脅迫綜合系數(shù)（comprehensive coefficient， CC）：共檢測到5個與Cd脅迫綜合系數(shù)相關(guān)的QTL，分別位于12、15、16、19和21號染色體上，解釋1.18%～9.79%的表型變異。其中，qCC-12-1和株高Cd脅迫系數(shù)位點qPH-12-1、發(fā)芽率Cd脅迫系數(shù)位點qGR-12-1、發(fā)芽勢Cd脅迫系數(shù)位點qGP-12-1為同一位點，qCC-15-1和株高Cd脅迫系數(shù)位點qPH-15-1為同一位點，qCC-16-1和地上部干物質(zhì)質(zhì)量Cd脅迫系數(shù)位點qSDW-16-1為同一位點，qCC-21-1和葉片SPAD值Cd脅迫系數(shù)位點qSPAD-21-1為同一位點。

3 討論

本研究構(gòu)建HM188群體的親本為陸地棉HS46和MARCABUCAGSUS-l-88，兩者親緣關(guān)系較遠(yuǎn)， HS46是高產(chǎn)優(yōu)質(zhì)陸地棉品種， MAR是多抗的陸地棉種質(zhì)資源。2個親本在農(nóng)藝性狀、纖維品質(zhì)和抗病性等方面都存在顯著的差異，但目前缺少關(guān)于該群體對低溫、高溫、干旱、澇漬和重金屬等非生物脅迫的耐性/抗性的研究。

前人已對HM188群體的抗病性、種子品質(zhì)、纖維品質(zhì)、產(chǎn)量等方面開展了遺傳效應(yīng)分析和QTL定位等工作。任愛霞等[17]調(diào)查了該群體枯萎病和黃萎病的田間發(fā)病情況，比較抗病、非抗病株系的比例，得出枯萎病抗性與黃萎病抗性分別為數(shù)量性狀和質(zhì)量性狀。欒啟福等[13]利用包括基因型與環(huán)境互作在內(nèi)的加性模型分析該群體，證實枯萎病抗病能力主要受加性效應(yīng)影響。孔廣超等[19]調(diào)查了該群體在4個環(huán)境下的產(chǎn)量及纖維品質(zhì)性狀的表型數(shù)據(jù)，通過連鎖分析檢測到109個與籽棉產(chǎn)量、皮棉產(chǎn)量、衣分、鈴重、籽指和衣指相關(guān)的QTL。秦利[20]將該群體株系間隨機(jī)組配，構(gòu)建了1套包括376個組合的永久F2群體，檢測到控制籽指、仁指、種仁率、蛋白質(zhì)含量、油分含量和棉酚含量相關(guān)的加性QTL 18個、上位性QTL 12個，解釋性狀2.42%～23.04%的表型變異。Li等[21]將2個親本HS46和MAR分別作為父本與HM188群體回交創(chuàng)建了2個陸地棉回交F1群體，定位到205個與產(chǎn)量性狀相關(guān)的QTL、195個與纖維品質(zhì)性狀相關(guān)的QTL、138個與產(chǎn)量性狀雜種優(yōu)勢相關(guān)的QTL和65個與纖維品質(zhì)性狀雜種優(yōu)勢相關(guān)的QTL。

耐Cd性苗期鑒定因環(huán)境問題而受到限制。本研究苗期試驗采用穴盤種植，通過在每個穴盤上種植1行中棉所49作為數(shù)據(jù)校準(zhǔn)依據(jù)，減少穴盤間差異導(dǎo)致的誤差。Cd脅迫綜合系數(shù)計算采用隸屬函數(shù)法，由于目前尚無文獻(xiàn)明確得出棉花各表型性狀的權(quán)重，參照王準(zhǔn)等[22]的方法，采用客觀賦權(quán)法計算各性狀的權(quán)重進(jìn)行綜合評價，為棉花耐重金屬脅迫鑒定提供借鑒。采用該方法對HM188群體的苗期耐Cd性進(jìn)行鑒定和QTL定位，其中qRFW-13-1位點在2019和2020年共2個環(huán)境下均被檢測到，qPH-12-1/qGR-12-1/qGP-

12-1、qPH-15-1、qSDW-16-1、qSPAD-21-1這4個位點在性狀Cd脅迫系數(shù)轉(zhuǎn)換為Cd脅迫綜合系數(shù)后同樣可以定位到，可以為陸地棉耐Cd相關(guān)候選基因的篩選和克隆提供參考。

4 結(jié)論

通過對HM188群體及雙親的發(fā)芽率、發(fā)芽勢以及5個苗期性狀的鑒定，發(fā)現(xiàn)HM188群體及雙親的耐Cd性變異豐富，其中強(qiáng)耐Cd型材料24份、耐Cd型材料49份、中等耐Cd型材料67份、Cd敏感型材料50份，為棉花耐Cd性的遺傳研究和品種選育提供了寶貴的種質(zhì)材料。通過連鎖分析檢測到28個與陸地棉苗期表型性狀Cd脅迫系數(shù)相關(guān)的QTL和5個與Cd脅迫綜合系數(shù)相關(guān)的QTL，填補(bǔ)了陸地棉耐Cd相關(guān)遺傳定位研究的空白，可為棉花耐Cd分子育種提供參考。但本研究僅針對陸地棉種子萌發(fā)期和苗期進(jìn)行脅迫研究，對于Cd脅迫對陸地棉產(chǎn)量、纖維品質(zhì)以及各組織Cd吸收量的研究需要在后續(xù)的試驗中進(jìn)行補(bǔ)充。

參考文獻(xiàn)：

[1]中華人民共和國環(huán)境保護(hù)部. 全國土壤污染狀況調(diào)查公報[EB/OL]. （2012-04-17） [2021-04-25]. http：//www.gov.cn/foot/

site1/20140417/782bcb88840814ba158d01.pdf.

Ministry of Environmental Protection of the People's Republic of China. Bulletin of national soil pollution survey[EB/OL]. （2012-

04-17） [2021-04-25]. http：//www.gov.cn/foot/site1/20140417/

782bcb88840814ba158d01.pdf.

[2] Li Y S， Tang H， Hu Y X， et al. Enrofloxacin at environmentally relevant concentrations enhances uptake and toxicity of cadmium in the earthworm Eisenia fetida in farm soils[J/OL]. Journal of Hazardous Materials， 2016， 308： 312-320[2021-04-25]. https：//

doi.org/10.1016/j.jhazmat.2016.01.057.

[3] Moynihan M， Peterson K E， Cantoral A， et al. Dietary predictors of urinary cadmium among pregnant women and children [J/OL]. Science of The Total Environment， 2017， 575： 1255-

1262[2021-04-25]. https：//doi.org/10.1016/j.scitotenv.2016.09.

204.

[4] 盧紅玲， 肖光輝， 劉青山， 等. 土壤Cd污染現(xiàn)狀及其治理措施研究進(jìn)展[J/OL]. 南方農(nóng)業(yè)學(xué)報， 2014， 45（11）： 1986-1993[2021-04-25]. https：//dio.org/10.3969/j.issn.2095-1191.2014.11.

1986.

Lu Hongling， Xiao Guanghui， Liu Qingshan， et al. Advances in soil Cd pollution and solution measures[J/OL]. Journal of Southern Agriculture， 2014， 45（11）： 1986-1993[2021-04-25]. https：//dio.org/10.3969/j.issn.2095-1191.2014.11.1986.

[5] 孫紅艷. 大麥耐Cd及Cd積累相關(guān)QTL定位與候選基因篩選及調(diào)控研究[D]. 杭州： 浙江大學(xué)， 2012.

Sun Hongyan. Identification of QTLs and candidate genes associated with cadmium tolerance and accumulation in barley and its regulation[D]. Hangzhou： Zhejiang University， 2012.

[6] 潘晨陽， 葉涵斐， 周維永， 等. 水稻籽粒Cd積累QTL定位及候選基因分析[J/OL]. 植物學(xué)報， 2021， 56（1）： 25-32[2021-04-25]. https：//doi.org/10.11983/CBB20148.

Pan Chenyang， Ye Hanfei， Zhou Weiyong， et al. QTL mapping of candidate genes involved in Cd accumulation in rice grain [J/OL]. Chinese Bulletin of Botany， 2021， 56（1）： 25-32[2021-

04-25]. https：//doi.org/10.11983/CBB20148.

[7] Kirungu J N， Magwanga R O， Shiraku M L， et al. Genome-wide identification and analysis of the CNGC gene family in upland cotton under multiple stress conditions[J/OL]. Journal of Cotton Research， 2023， 6： 17[2024-04-25]. https：//doi.org/10.1186/

s42397-023-00152-z.

[8] 郭利雙， 陳浩東， 賀云新， 等. Cd高積累棉花品種水培篩選結(jié)果初報[J/OL]. 中國棉花， 2015， 42（10）： 14-16[2021-04-25]. https：//

doi.org/10.11963/issn.1000-632X.201510005.

Guo Lishuang， Chen Haodong， He Yunxin， et al. Preliminary results on cotton cultivars screening for high accumulation of cadmium[J/OL]. China Cotton， 2015， 42（10）： 14-16[2021-04-

25]. https：//doi.org/10.11963/issn.1000-632X.201510005.

[9] 李玲， 陳進(jìn)紅， 何秋伶， 等. 3個陸地棉種質(zhì)系重金屬Cd的積累、轉(zhuǎn)運和富集特性分析[J/OL]. 棉花學(xué)報， 2012， 24（6）： 535-

540[2021-04-25]. https：//doi.org/10.3969/j.issn.1002-7807.2012.

06.010.

Li Ling， Chen Jinhong， He Qiuling， et al. Accumulation， transportation， and bioconcentration of cadmium in three upland cotton plants under cadmium stress[J/OL]. Cotton Science， 2012， 24（6）： 535-540[2021-04-25]. https：//doi.org/10.3969/j.issn.1002-

7807.2012.06.010.

[10] 樊霆， 葉文玲， 陳海燕， 等. 農(nóng)田土壤重金屬污染狀況及修復(fù)技術(shù)研究[J]. 生態(tài)環(huán)境學(xué)報， 2013， 22（10）： 1727-1736.

Fan Ting， Ye Wenling， Chen Haiyan， et al. Review on contamination and remediation technology of heavy metal in agricultural soil[J]. Ecology and Environmental Sciences， 2013， 22（10）： 1727-1736.

[11] 秦普豐， 鐵柏清， 周細(xì)紅， 等. 鉛與Cd對棉花和水稻萌發(fā)及生長的影響[J/OL]. 湖南農(nóng)業(yè)大學(xué)學(xué)報， 2000， 26（3）： 205-207[2021-04-25]. https：//doi.org/10.3321/j.issn：1007-1032.2000.

03.013.

Qin Pufeng， Tie Baiqing， Zhou Xihong， et al. Effects of cadmium and lead in soil on the germination and growth of rice and cotton[J/OL]. Journal of Hunan Agricultural University， 2000， 26（3）： 205-207[2021-04-25]. https：//doi.org/10.3321/j.issn：1007-

1032.2000.03.013.

[12] 劉連濤， 陳靜， 孫紅春， 等. Cd脅迫對棉花幼苗生長效應(yīng)及不同器官Cd積累的影響[J/OL]. 棉花學(xué)報， 2014， 26（5）： 466-

470[2021-04-25]. https：//doi.org/10.3969/j.issn.1002-7807.2014.

05.012.

Liu Liantao， Chen Jing， Sun Hongchun， et al. Effects of cadmium stress on growth and cadmium accumulation in cotton （Gossy-pium hirsutum L.） seedlings[J/OL]. Cotton Science， 2014， 26（5）： 466-470[2021-04-25]. https：//doi.org/10.3969/j.issn.1002-7807.

2014.05.012.

[13] 欒啟福， 祝水金. 陸地棉重組近交系HM188及其性狀表現(xiàn)[J]. 棉花學(xué)報， 2003， 15（4）： 252-253.

Luan Qifu， Zhu Shuijin. Introduction of a recombinant inbred line of upland cotton （G. hirsutum L.） HM188[J]. Cotton Science， 2003， 15（4）： 252-253.

[14] 朱協(xié)飛， 司占峰. 棉花導(dǎo)入系耐鹽性鑒定及耐鹽基因QTL定位[J/OL]. 棉花學(xué)報， 2019， 31（1）： 23-30[2021-04-25]. https：//

doi.org/10.11963/1002-7807.zxfzxf.20181226.

Zhu Xiefei， Si Zhanfeng. Evaluation and QTL mapping of tolerance to salinity using interspecific introgression lines from Gossypium barbadense in G. hirsutum[J/OL]. Cotton Science， 2019， 31（1）： 23-30[2021-04-25]. https：//doi.org/10.11963/1002-

7807.zxfzxf.20181226.

[15] 岳松青， 曹輝， 荀咪， 等. 基于指標(biāo)綜合分析的蘋果砧木耐Cd性評價[J]. 植物生理學(xué)報， 2019， 55（5）： 649-656.

Yue Songqing， Cao Hui， Xun Mi， et al. Evaluation of cadmium tolerance for apple rootstocks based on index comprehensive analysis[J]. Plant Physiology Journal， 2019， 55 （5）： 649-656.

[16] Li C， Dong Y T， Zhao T L， et al. Genome-wide SNP linkage mapping and QTL analysis for fiber quality and yield traits in the upland cotton recombinant inbred lines population[J/OL]. Frontiers in Plant Science， 2016， 7： 1356[2021-04-25]. https：//

doi.org/10.3389/fpls.2016.01356.

[17] McCouch S R， Chen X L， Panaud O， et al. Microsatellite marker development， mapping and applications in rice genetics and breeding[J]. Plant Molecular Biology， 1997， 35（1/2）： 89-99.

[18] 任愛霞， 胡家恕， 祝水金， 等. 棉花黃萎病抗性與過氧化物酶同工酶分析[J]. 棉花學(xué)報， 2002， 14（5）： 273-276.

Ren Aixia， Hu Jiasu， Zhu Shuijin， et al. Analysis of peroxidase isozymes and resistance to Verticillium dahliae in cotton[J]. Cotton Science， 2002， 14（5）： 273-276.

[19] 孔廣超， 秦利， 徐海明， 等. 棉花IF2群體構(gòu)建及其在纖維品質(zhì)遺傳和雜種優(yōu)勢研究中的應(yīng)用[J]. 作物學(xué)報， 2010， 36（6）： 940-944.

Kong Guangchao， Qin Li， Xu Haiming， et al. Construction of IF2 population and its application in studies on genetic effects and heterosis for fiber quality in upland cotton （G. hirsutum L.）[J]. Acta Agronomica Sinica， 2010， 36（6）： 940-944.

[20] 秦利. 陸地棉種子品質(zhì)性狀遺傳及其QTL定位研究[D]. 杭州： 浙江大學(xué)， 2009.

Qin Li. Inheritance of seed quality traits and their QTL mapping in upland cotton （G. hirsutum L.）[D]. Hangzhou： Zhejiang University， 2009.

[21] Li C， Yu H R， Li C， et al. QTL Mapping and heterosis analysis for fiber quality traits across multiple genetic populations and environments in upland cotton[J/OL]. Frontiers in Plant Science， 2018， 9： 1364[2021-04-25]. https：//doi.org/ 10.3389/fpls.2018.

01364.

[22] 王準(zhǔn)， 張恒恒， 董強(qiáng)， 等. 棉花耐低氮和氮敏感種質(zhì)篩選及驗證[J/OL]. 棉花學(xué)報， 2020， 32（6）： 538-551[2021-04-25]. https：//doi.org/10.11963/1002-7807.wzsmz.20201023.

Wang Zhun， Zhang Hengheng， Dong Qiang， et al. Screening and verification of low nitrogen tolerant and nitrogen sensitive cotton germplasm[J/OL]. Cotton Science， 2020， 32（6）： 538-551[2021-04-25]. https：//doi.org/10.11963/1002-7807.wzsmz.20201023.

（責(zé)任編輯：王國鑫" " 責(zé)任校對：王小璐）" ●