，，，，，，*

(1.石河子大學(xué),新疆 石河子832003；2.中國(guó)農(nóng)業(yè)科學(xué)院北京畜牧獸醫(yī)研究所,北京100193)

開花性狀是紫花苜蓿(Medicagosativa)重要的農(nóng)藝性狀，產(chǎn)量積累主要發(fā)生在開花之前，通過(guò)推遲始花期能夠提高產(chǎn)量[1]。開花是植物由營(yíng)養(yǎng)生長(zhǎng)向生殖生長(zhǎng)轉(zhuǎn)變的重要標(biāo)志之一，對(duì)植物的有性生殖起著決定性的作用[2]。始花期與苜蓿干草的產(chǎn)量和品質(zhì)密切相關(guān)，一直是苜蓿遺傳育種過(guò)程中的主要關(guān)注指標(biāo)之一。掌握始花期遺傳規(guī)律能夠幫助管理者選擇最適合的放牧?xí)r期和收獲期，從而獲得更高的經(jīng)濟(jì)效益[3]。與其他作物相比，苜蓿開花性狀的相關(guān)報(bào)道比較少，這主要與紫花苜蓿遺傳性狀的復(fù)雜性及相關(guān)分析技術(shù)的限制有關(guān)。苜蓿為無(wú)限花序，一般群體的開花時(shí)間可以延續(xù)40～60 d。將始花期作為一個(gè)苜蓿開花過(guò)程的基準(zhǔn)，有助于更加精細(xì)地描述、分析和預(yù)測(cè)植物的發(fā)育。始花期屬于多基因控制的數(shù)量性狀，由遺傳因素和環(huán)境因素共同決定[4]。研究始花期性狀的遺傳構(gòu)成對(duì)深入研究開花性狀的遺傳規(guī)律具有重要意義。前人利用紫花苜蓿雜交后代研究不同農(nóng)藝性狀之間的關(guān)系已有報(bào)道，通過(guò)調(diào)控始花期能夠獲得更高的產(chǎn)量[5]。將始花期具有差異的材料進(jìn)行雜交能夠獲得始花期性狀差異明顯的F1代材料[6]。但是從表型性狀進(jìn)行紫花苜蓿早熟遺傳的相關(guān)研究國(guó)內(nèi)還未見報(bào)道，而通過(guò)IECM算法計(jì)算主+多基因混合遺傳模型是分析數(shù)量性狀遺傳規(guī)律的有效方法[7]。在紫花苜蓿研究中，利用主多基因進(jìn)行表型性狀分析的研究還比較少，這主要是因?yàn)橹鞫嗷蛞话銘?yīng)用于二倍體聯(lián)合世代分析[8]。但是利用該方法分析F2代遺傳模型的研究也有報(bào)道[9]，這說(shuō)明利用主多基因模型進(jìn)行遺傳分析也可以應(yīng)用于單個(gè)世代。曹錫文[10]開發(fā)出主多基因SEA軟件，這也為主多基因分析提供了極大的便利條件。

紫花苜蓿為同源四倍體，其雜合水平高，導(dǎo)致多種構(gòu)圖群體存在偏差。由于多倍體復(fù)雜的遺傳特性和相關(guān)的輔助構(gòu)圖軟件的限制，使得遺傳圖譜的構(gòu)建相對(duì)困難。借助二倍體分析軟件對(duì)紫花苜蓿進(jìn)行研究是一種有效的方法[11]。分子遺傳連鎖圖譜的構(gòu)建是數(shù)量性狀基因定位、分子標(biāo)記輔助育種以及基因克隆的重要理論基礎(chǔ), 利用遺傳連鎖圖譜尋找分子標(biāo)記與控制優(yōu)良性狀基因之間的連鎖關(guān)系，對(duì)植物基礎(chǔ)遺傳育種研究具有重要的應(yīng)用價(jià)值[12]。王夢(mèng)穎等[13]對(duì)早熟材料產(chǎn)生的雜交F1群體進(jìn)行分析，初步構(gòu)建出紫花苜蓿遺傳連鎖圖譜。賈瑞等[14]也利用15個(gè)雜交組合產(chǎn)生的F1代單株進(jìn)行了生物學(xué)性狀測(cè)定。近年來(lái)植物遺傳圖譜的構(gòu)建研究相對(duì)來(lái)說(shuō)比較多，北柴胡(Bupleurumchinense)[15]、甘薯(Dioscoreaesculenta)[16]等植物的遺傳圖譜已經(jīng)被構(gòu)建。但是利用紫花苜蓿進(jìn)行遺傳連鎖圖譜研究的報(bào)道仍然較少，使得相關(guān)QTL定位分析相對(duì)困難，難以加快苜蓿的育種進(jìn)程。GBS技術(shù)的出現(xiàn)彌補(bǔ)了分子標(biāo)記數(shù)量的限制，較傳統(tǒng)的標(biāo)記檢測(cè)方法來(lái)說(shuō)能夠獲得更多的SNP標(biāo)記，提高遺傳連鎖圖譜的密度。本研究利用紫花苜?；ㄆ诓町惷黠@的F1代群體進(jìn)行相關(guān)研究，并結(jié)合GBS測(cè)序技術(shù)進(jìn)行遺傳分析，重點(diǎn)分析始花期性狀的遺傳特性，并進(jìn)行早熟性狀的QTL定位分析，以期利用植物遺傳連鎖圖以及QTL定位方法，為解析紫花苜蓿開花性狀的遺傳規(guī)律、培育高產(chǎn)苜蓿新品種提供理論依據(jù)。

1 材料與方法

1.1 供試材料

供試材料為中苜系列品種選育過(guò)程中的中間材料：低產(chǎn)早熟紫花苜蓿(父本，來(lái)源于滄州苜蓿)、高產(chǎn)晚熟紫花苜蓿(母本，來(lái)源于中苜1號(hào) )單株個(gè)體及二者雜交產(chǎn)生的152個(gè)F1代單株個(gè)體。

1.2 試驗(yàn)設(shè)計(jì)

雜交群體表型觀測(cè)試驗(yàn)于2014-2016年在中國(guó)農(nóng)業(yè)科學(xué)院國(guó)際農(nóng)業(yè)高新技術(shù)產(chǎn)業(yè)園基地(河北省廊坊市境內(nèi))進(jìn)行。該地位于河北省西北部，屬暖溫帶大陸性季風(fēng)氣候, 年平均氣溫11.9 ℃左右, 最冷月份(1月)平均氣溫為-4.7 ℃, 最熱月份(7月)平均氣溫為26.2 ℃。年降水量554.9 mm。土質(zhì)為中壤土，含有機(jī)質(zhì)1.69%，pH值7.37。

雜交種子在溫室培養(yǎng)獲得F1代幼苗，對(duì)幼苗進(jìn)行DNA提取并按照王夢(mèng)穎等[13]的方法利用SSR分子標(biāo)記鑒定雜交種，最終獲得152個(gè)F1代單株個(gè)體。通過(guò)扦插獲得親本和F1代單株個(gè)體的3個(gè)重復(fù)。采用隨機(jī)區(qū)組排列進(jìn)行試驗(yàn)，設(shè)3次重復(fù)，行距1 m，株距0.8 m。在2014年入冬前進(jìn)行一次刈割，從而保證不同單株間的一致性。于2015和2016年進(jìn)行始花期的調(diào)查。每年從5月1日起統(tǒng)計(jì)不同單株的開花情況，始花期計(jì)算從連續(xù)5 d日平均溫度大于10 ℃算起(2015年3月10日；2016年3月16日)，直至第一朵小花出現(xiàn)即表示開花所需天數(shù)。

1.3 連鎖圖譜的構(gòu)建

通過(guò)對(duì)測(cè)序數(shù)據(jù)進(jìn)行分析，獲得20334個(gè)親本間具有差異的SNP位點(diǎn)。按照Li等[11]對(duì)紫花苜蓿遺傳連鎖圖譜構(gòu)建過(guò)程中的分析，僅利用親本間的單雜合位點(diǎn)分析來(lái)解決同源四倍體分離復(fù)雜的簡(jiǎn)便算法，并且具有較高的準(zhǔn)確性[17]。本研究篩選出缺失值小于10%的單雜合SNP位點(diǎn)(AAAAXAAAB型)，得到7468個(gè)SNP位點(diǎn)(雜合父本5888個(gè)，雜合母本1580個(gè))。最后篩選出具有子代分離比為1∶1(P>0.001)的候選位點(diǎn)1320個(gè)，雜合父本位點(diǎn)1011個(gè)，雜合母本位點(diǎn)309個(gè)(表4，表5)。以這些SNP位點(diǎn)為遺傳標(biāo)記進(jìn)行連鎖圖譜構(gòu)建及QTL定位。

1.4 數(shù)據(jù)處理

將2015，2016年的154個(gè)單株的3個(gè)重復(fù)均值進(jìn)行基本統(tǒng)計(jì)分析，并利用3個(gè)重復(fù)的親本數(shù)據(jù)進(jìn)行親本差異性T檢驗(yàn)；運(yùn)用主多基因遺傳模型分析遺傳率。利用R繪制正態(tài)分布圖。根據(jù)植物數(shù)量性狀主多基因混合模型分析方法對(duì)不同表型性狀進(jìn)行分析。采用極大似然法和IECM算法估計(jì)各世代、各成分分布的參數(shù)，利用AIC值，U12、U22、U32檢驗(yàn)，Smimov檢驗(yàn)(nW2)，Kolmogorov檢驗(yàn)(Dn)和Heritability參數(shù),選擇最優(yōu)遺傳模型。使用SEA軟件包進(jìn)行植物數(shù)量性狀主基因+多基因遺傳體系分析。遺傳連鎖圖譜構(gòu)建主要利用GBS測(cè)序產(chǎn)生的SNP標(biāo)記，根據(jù)親本和F1雜交后代的分離類型進(jìn)行標(biāo)記分析。將篩選到的SNP標(biāo)記導(dǎo)入Joinmap軟件中，以LOD=3為標(biāo)準(zhǔn)構(gòu)建連鎖圖譜，最終構(gòu)建出紫花苜蓿的遺傳連鎖圖譜。將表型數(shù)據(jù)和基因型數(shù)據(jù)導(dǎo)入QTL IciMapping中進(jìn)行QTL定位分析，軟件參數(shù)為默認(rèn)參數(shù)。

2 結(jié)果與分析

2.1 親本及雜交F1代始花期數(shù)據(jù)的基本統(tǒng)計(jì)結(jié)果

由表1結(jié)果可知，兩年試驗(yàn)結(jié)果具有一致性，并且試驗(yàn)重現(xiàn)性較好。親本間始花期數(shù)據(jù)差異明顯，平均開花時(shí)間相差達(dá)13 d，T檢驗(yàn)表明親本間差異達(dá)到極顯著水平，同時(shí)Levene’s 檢驗(yàn)證明親本間方差具有同質(zhì)性。F1代均值介于雙親之間，但是變異系數(shù)不大。F1代變異范圍超過(guò)雙親差異范圍，這表現(xiàn)出超親分布的特點(diǎn)。

表1 親本和雜交F1代的始花期性狀基本統(tǒng)計(jì)分析Table 1 Summary statistics analysis of phenotypes for early flowering-time in the F1 progeny and parents

從圖1的概率密度分布圖中可知，大部分單株開花時(shí)間集中在一個(gè)區(qū)域，但是有少部分植株開花時(shí)間集中在另一個(gè)區(qū)域。父本和母本的始花期差異達(dá)到極顯著水平(P<0.01)。同時(shí)父本開花較早，母本開花較晚，親本差異明顯。另外152個(gè)F1代植株的開花期差別也很大，最早開花和最晚開花植株間差異能夠達(dá)到20 d，同時(shí)兩年始花期數(shù)據(jù)結(jié)果具有相似的分離特點(diǎn)。

正態(tài)性檢驗(yàn)(Kolmogorov-Smirnov)表明,F1代單株分布為非正態(tài)分布(P<0.05)，圖1表明F1代單株的分布呈現(xiàn)出雙峰分布的特點(diǎn)，黑色實(shí)線為正態(tài)分布曲線。

圖1 2015年和2016年始花期的概率密度分布圖Fig.1 Density estimate of initial time of flowering in 2015 and 2016 P1, P2分別代表父本，母本。P1, P2 represent paternal and maternal, respectively.

2.2 適宜遺傳模型選擇及遺傳參數(shù)估算

利用數(shù)量性狀主多基因分析法對(duì)F1代單株兩年始花期數(shù)據(jù)進(jìn)行基因聯(lián)合分析，通過(guò)極大似然法和IECM算法估算AIC值。根據(jù)最適遺傳模型篩選原則，篩選出AIC值最小和較小者，2015年篩選出2MG-A和2MG-ADI兩個(gè)最適遺傳模型(表2)。并對(duì)候選模型進(jìn)行適合性檢驗(yàn)，結(jié)果顯示2MG-A的nW2檢驗(yàn)達(dá)到顯著水平，其他檢驗(yàn)未達(dá)到顯著水平。而2MG-ADI的U32和nW2檢驗(yàn)達(dá)到顯著水平，這兩個(gè)模型都具有兩對(duì)主效基因，同時(shí)適合性檢驗(yàn)結(jié)果類似。因此2MG-A和2MG-ADI都可作為候選遺傳模型。2016年篩選出2MG-ADI和2MG-A兩個(gè)最適遺傳模型(表3)。并對(duì)候選模型進(jìn)行適合性檢驗(yàn)，結(jié)果顯示2MG-ADI的nW2檢驗(yàn)達(dá)到顯著水平，其他檢驗(yàn)未達(dá)到顯著水平。而2MG-A的U22和nW2檢驗(yàn)達(dá)到顯著水平，這兩個(gè)模型都具有兩對(duì)主效基因，同時(shí)適合性檢驗(yàn)結(jié)果類似。因此2MG-A和2MG-ADI都可作為候選遺傳模型，兩年試驗(yàn)結(jié)果具有一致性，始花期數(shù)據(jù)由兩對(duì)主效基因控制，并且兩年結(jié)果具有相同的遺傳模型。

表2 2015年始花期的遺傳參數(shù)估計(jì)Table 2 Estimates of genetic parameters for early flowering-time in 2015

MG:主基因;A:加性效應(yīng);D:顯性效應(yīng);I:上位性效應(yīng);E:相等；N：負(fù)向。U12、U22、U32、nW2、Dn 分別指均勻性U12、U22、U32檢驗(yàn), Smimov檢驗(yàn)和Kolmogorov檢驗(yàn)。Heritability：主基因遺傳率。AIC:Akaike信息準(zhǔn)則。下同。

MG: Major gene model; A: Additive effect; D: Dominance effect; I: Epistatic interaction; E: Equal; N：Negative. U12、U22、U32、nW2、Dn represent the uniform test U12, U22, U32, Smimov test and Kolmogorov test, respectively. Heritability: Heritability of major gene. AIC:Akaike’s information criterion. The same below.

表3 2016年始花期的遺傳參數(shù)估計(jì)Table 3 Estimates of genetic parameters for early flowering-time in 2016

根據(jù)不同遺傳模型的極大似然估計(jì)值，確定不同遺傳模型的遺傳參數(shù)(表2和表3)。同時(shí)根據(jù)AIC值和適合性檢驗(yàn)得到最適遺傳模型為2MG-A和2MG-ADI。2015年試驗(yàn)結(jié)果表明，二者主基因遺傳率都達(dá)到98%以上，但是2016年結(jié)果僅有2MG-A達(dá)到98%，因此本研究判斷候選模型為2MG-A，主基因的遺傳率為99%。始花期主要受兩對(duì)主效基因控制，同時(shí)具有加性作用。

2.3 遺傳圖譜的構(gòu)建

將上步篩選到的SNP標(biāo)記(20334個(gè))導(dǎo)入Joinmap軟件中，按照SNP標(biāo)記所在的染色體進(jìn)行連鎖群設(shè)定標(biāo)準(zhǔn)，以LOD=3為標(biāo)準(zhǔn)構(gòu)建連鎖圖譜，進(jìn)而以重組率小于0.4和LOD值大于1進(jìn)行SNP標(biāo)記的排序和篩選,最終得到遺傳圖譜。父本連鎖圖共有430個(gè)SNP標(biāo)記，覆蓋圖距1386 cM，平均標(biāo)記間圖距3.2 cM,連鎖群長(zhǎng)度在122.696～208.514 cM，其中g(shù)roup 8連鎖群最短(表4和圖2)。母本連鎖圖共有99個(gè)SNP標(biāo)記，覆蓋圖距798.73 cM, 平均圖距8.07 cM，連鎖群長(zhǎng)度在54.476～148.826 cM，其中g(shù)roup 6連鎖群最短(表5和圖3)。

表5 母本連鎖群基本信息Table 5 Maternal chain group basic information

圖2 父本遺傳連鎖圖Fig.2 Paternal genetic linkage map cM表示重組頻率的測(cè)量單位。左側(cè)表示SNP標(biāo)記在染色體上的相對(duì)位置。右側(cè)代表父本SNP編號(hào)。cM represents the unit of measurement of the recombination frequency. The left side represents the relative position of the SNP on the chromosome. The right side represents the paternal SNP number.

圖3 母本遺傳連鎖圖Fig.3 Maternal genetic linkage map cM表示重組頻率的測(cè)量單位。左側(cè)表示SNP標(biāo)記在染色體上的相對(duì)位置。右側(cè)代表母本SNP標(biāo)記的編號(hào)。cM represents the unit of measurement of the recombination frequency. The left side represents the relative position of the SNP on the chromosome. The right represents the maternal SNP number.

2.4 QTL定位

利用SNP標(biāo)記信息和2年表型信息進(jìn)行QTL定位，最終兩年數(shù)據(jù)分別定位到2個(gè)QTL位點(diǎn)。其中2015年結(jié)果定位到的QTL位點(diǎn)位于3號(hào)連鎖群上，2016年定位到的QTL位點(diǎn)位于2號(hào)連鎖群上，兩年結(jié)果并沒(méi)有相同的QTL位點(diǎn)(表6)。兩個(gè)QTL位點(diǎn)的LOD值分別為3.1299 和3.6756，貢獻(xiàn)率分別為12.1334%和11.0157%。

表6 父本QTL定位結(jié)果Table 6 Paternal QTL mapping results

PVE：Phenotypic variation explained by QTL at the current position.

3 討論

3.1 主多基因模型分析的優(yōu)勢(shì)

主多基因模型是用傳統(tǒng)方法預(yù)測(cè)表型變異信息的一種有效方法，它在許多植物中都有應(yīng)用。周清元等[18]采用主+多基因混合模型預(yù)測(cè)甘藍(lán)型油菜(Brassicanapus)6世代變異信息，得出果身長(zhǎng)、角果長(zhǎng)、果喙長(zhǎng)等性狀的遺傳模型，并根據(jù)不同遺傳模型的參數(shù)結(jié)果估計(jì)出加性效應(yīng)和顯性效應(yīng)值，進(jìn)而從不同模型間的相互作用中推斷出遺傳進(jìn)化過(guò)程。李忠南等[19]對(duì)玉米(Zeamays)6世代進(jìn)行聯(lián)合分析得出，葉綠素SPAD值主要受兩對(duì)主基因+多基因控制。同時(shí)該遺傳模型也可以應(yīng)用于動(dòng)物研究[20],這些模型適合性檢驗(yàn)結(jié)果良好，因此該方法非常適合應(yīng)用于遺傳分析研究。同時(shí)該方法已經(jīng)被沿用多年，在預(yù)測(cè)準(zhǔn)確性方面有一定的優(yōu)勢(shì)，且各種試驗(yàn)結(jié)果良好。因此該方法適合作為關(guān)聯(lián)分析的基礎(chǔ)步驟，通過(guò)表型預(yù)測(cè)初步估計(jì)遺傳變異信息，為進(jìn)一步的遺傳分析研究提供參考依據(jù)。

兩年試驗(yàn)結(jié)果，能夠通過(guò)不同年份的試驗(yàn)結(jié)果信息消除環(huán)境因素，進(jìn)一步通過(guò)不同環(huán)境中穩(wěn)定存在的主多基因模型預(yù)測(cè)始花期性狀的遺傳變異規(guī)律。兩年數(shù)據(jù)結(jié)果具有一致性，始花期性狀的最適遺傳模型為2MG-A和2MG-ADI，這兩個(gè)遺傳模型在不同年份分別具有最小的AIC值，同時(shí)主基因遺傳率能夠解釋98%以上的表型變異信息。此外兩個(gè)遺傳模型都是由2對(duì)主效基因控制，因此可推斷始花期性狀由兩對(duì)主效基因控制。但是2MG-ADI的主基因遺傳率在2016年為0.45%，因此并不是最優(yōu)遺傳模型。而2MG-A遺傳模型在兩年結(jié)果中都能解釋98%的表型變異信息，因此始花期性狀的遺傳信息可以被2MG-A遺傳模型較好解釋，同時(shí)始花期性狀受環(huán)境影響因素較小，該性狀主要被遺傳因素影響，相關(guān)遺傳信息能夠在不同環(huán)境中穩(wěn)定表達(dá)。

3.2 QTL定位

紫花苜蓿為同源四倍體，異花授粉植物，其雜合水平高，因此導(dǎo)致多種構(gòu)圖群體存在偏差。由于多倍體復(fù)雜的遺傳特性和相關(guān)的輔助構(gòu)圖軟件的限制，使得遺傳圖譜的構(gòu)建相對(duì)困難[21]。在四倍體苜蓿中已經(jīng)建立了一些基因連鎖圖譜。應(yīng)用QTL作圖，一些重要的農(nóng)藝性狀，如產(chǎn)量[22-23]、耐寒性等[24-25]QTL已鑒定出來(lái)。同時(shí)已有大量試驗(yàn)證實(shí)在蒺藜苜蓿(Medicagotruncatula)中鑒定出的QTL，對(duì)苜蓿的改良具有應(yīng)用價(jià)值。如蒺藜苜蓿的遺傳圖譜在抗春季黑莖病和葉斑病的QTL中被確定[26]，而且與氮素營(yíng)養(yǎng)[27-28]及形態(tài)特征[29-30]相關(guān)的QTL也被鑒定出。因此，蒺藜苜蓿QTL定位與遺傳圖譜構(gòu)建的研究為四倍體苜蓿的相關(guān)研究奠定了良好的基礎(chǔ)[26-31]。目前只有為數(shù)不多的幾個(gè)完整的紫花苜蓿的遺傳圖譜。最早的報(bào)道是用F1群體，分析32個(gè)RAPD標(biāo)記的分離結(jié)果，發(fā)現(xiàn)了9個(gè)連鎖群，并且構(gòu)建了四倍體苜蓿的遺傳圖譜[32]。Irwin等[33]利用RAPD和AFLP標(biāo)記對(duì)紫花苜蓿親本進(jìn)行抗炭疽病研究, 共得到了10個(gè)與抗炭疽顯著相關(guān)的標(biāo)記。Musial等[34]運(yùn)用RAPD、AFLP和SSR技術(shù)從紫花苜蓿的遺傳圖譜中得到與產(chǎn)量正相關(guān)的標(biāo)記及與產(chǎn)量負(fù)相關(guān)的基因。然而，大多數(shù)分子標(biāo)記，如SSR等不能夠滿足基因的精細(xì)定位和全基因組關(guān)聯(lián)研究[17]。SNP標(biāo)記構(gòu)建遺傳圖譜，具有共顯性，位點(diǎn)穩(wěn)定等優(yōu)點(diǎn)，是構(gòu)建遺傳圖譜較為理想的分子標(biāo)記之一。本試驗(yàn)通過(guò)對(duì)測(cè)序數(shù)據(jù)進(jìn)行分析，獲得20334個(gè)親本間具有差異的SNP位點(diǎn)。最后篩選出具有子代分離比為1∶1(P>0.001)的候選位點(diǎn)1320個(gè)，雜合父本位點(diǎn)1011個(gè)，雜合母本位點(diǎn)309個(gè)，最終構(gòu)建了包含8個(gè)連鎖群的遺傳圖譜。父本遺傳圖譜覆蓋圖距為1386 cM，標(biāo)記間平均圖距3.2 cM, 母本覆蓋圖距798.73 cM, 平均圖距8.07 cM。同時(shí)定位到早熟性狀相關(guān)的兩個(gè)QTL位點(diǎn)。本結(jié)論和前人研究結(jié)果存在一定差異[11]，這主要是分析方法的差別。本研究利用蒺藜苜蓿作為參考基因組進(jìn)行SNP分型，這樣進(jìn)行分析的優(yōu)點(diǎn)是SNP在染色體上的位置能夠確定，但是定位到染色體上的SNP標(biāo)記較少。結(jié)合性狀鑒定結(jié)果，本研究定位到早熟性狀QTL位點(diǎn)，并且表型解釋率較高，因此可以說(shuō)明本研究的分析方法具有現(xiàn)實(shí)可行性。同時(shí)也說(shuō)明這兩個(gè)標(biāo)記附近是控制開花性狀的重要區(qū)域，這些區(qū)域的發(fā)現(xiàn)為紫花苜蓿早熟育種的分子標(biāo)記輔助選擇提供了有利證據(jù)。

4 結(jié)論

1)利用數(shù)量性狀主多基因分析法對(duì)紫花苜蓿雜交F1代單株的兩年始花期數(shù)據(jù)進(jìn)行基因聯(lián)合分析，篩選出2MG-A為最適遺傳模型。2015年的主基因遺傳率為99%，2016年的主基因遺傳率為98.5%。始花期主要受兩對(duì)主效基因控制，同時(shí)具有加性作用。

2)構(gòu)建了紫花苜蓿遺傳連鎖圖譜，其中父本遺傳圖譜覆蓋圖距為1386 cM，標(biāo)記間平均圖距3.2 cM；母本覆蓋圖距798.73 cM, 平均圖距8.07 cM。QTL定位分析獲得2個(gè)早熟相關(guān)QTL位點(diǎn)，LOD值分別為3.1299 和3.6756，貢獻(xiàn)率分別為12.1334%和11.0157%。