趙小珍，趙衛(wèi)國，張春，余坤江，彭門路，陳鋒，張維，孫程明，李保軍，王灝，王曉東*，張潔夫*

（1.南京農(nóng)業(yè)大學(xué)作物遺傳與種質(zhì)創(chuàng)新國家重點實驗室，江蘇 省現(xiàn)代作物生產(chǎn)協(xié)同創(chuàng)新中心，江蘇 南京，210095；2.江蘇省農(nóng)業(yè)科學(xué)院經(jīng)濟作物研究所，農(nóng)業(yè) 農(nóng)村部長江下游棉花與油菜重點實驗室，江蘇 南京，210014；3.陜西省雜交油菜研究中心，陜西 楊凌，712100）

油菜是我國第一大油料作物，菜籽油占國產(chǎn)食用植物油的47%以上[1]。隨著生活水平的提高，人們對油脂的需求不斷增強，進一步提高油菜產(chǎn)量變得尤為重要。分枝角度指的是有效分枝與主莖形成的夾角，是影響油菜產(chǎn)量的重要株型性狀[2]。研究表明，適度緊湊的株型有利于油菜中下部的通風(fēng)透光、降低病蟲害、提高群體光能利用率、提高葉面積系數(shù)和抗倒伏性，從而實現(xiàn)高產(chǎn)[3,4]。

植物分枝角度的形成是其感受器官受到重力刺激后，迫使各器官與重力形成最適角度的過程[5]，重力反應(yīng)的異常會導(dǎo)致分枝角度的異常[6]。水稻（Oryza sativaL.）的分蘗角度[7]及玉米（Zea maysL.）的莖葉夾角[8]等均受到重力反應(yīng)影響。當(dāng)前，人們普遍認(rèn)為重力反應(yīng)分為4 個步驟，即重力信號的感受、轉(zhuǎn)導(dǎo)、重力信號造成的生長素差異性分布及生長角度的形成[9]。目前就植物重力感受機制得到普遍認(rèn)可的是淀粉-平衡石假說[10]，該假說認(rèn)為平衡石細(xì)胞的淀粉體能夠感受重力刺激的改變并作出應(yīng)激反應(yīng)[11]，這里的平衡石細(xì)胞是一種富含淀粉體的高度極化細(xì)胞。植物在感受到重力刺激后隨之引起體內(nèi)信號轉(zhuǎn)導(dǎo)，研究表明，Ca2+、InsP3（1,4,5-三磷酸肌醇）、pH 等信號分子獨自或協(xié)作參與調(diào)控了該過程[12]。最后，通過將不同信號分子傳遞到反應(yīng)器官來激活下游底物發(fā)揮功能，進而完成重力反應(yīng)。Cholodny-Went 假說[13]認(rèn)為重力刺激引起生長素不對稱性分布，從而引起重力反應(yīng)部位差異性生長，其在擬南芥重力實驗中得到了證實[14]。除了生長素，乙烯[15]、油菜素內(nèi)酯[16]等激素也參與了重力反應(yīng)調(diào)控。

1938 年水稻lazy1（la1）突變體的研究首次分析了其分枝（分蘗）角度調(diào)控機理，該突變體重力反應(yīng)缺陷導(dǎo)致后期分蘗平臥生長[17]；2007年la1突變體的基因LAZY1被成功克隆[18]。除了LAZY1基因，控制水稻分枝（分蘗）角度的基因TAC1（TILLER ANGLE CONTROL 1）[19]、PROG1（PROSTRATE GROWTH 1）[20]、LPA1（LOOSE PLANT ARCHITECTURE 1）[21]及PAY1（PLANT ARCHITECTURE AND YIELD 1）[22]等相繼被克隆出來?；騆AZY1在玉米和擬南芥中的同源基因ZmLA1[8]和AtLAZY1[23]，基因TAC1在玉米和桃樹同源拷貝ZmTAC1[24]和PpeTAC1[25]等，也相繼被克隆。截止目前，部分分枝角度、分蘗角度、莖葉夾角候選基因已在水稻、擬南芥和玉米中被克隆和鑒定，但在油菜上研究較少，尚未有基因克隆的報道。

近年來，QTL 定位和GWAS 分析被廣泛應(yīng)用于識別復(fù)雜的數(shù)量性狀位點，并成功應(yīng)用于甘藍型油菜分枝角度性狀的研究中。張倩[26]利用F2群體定位到1 個分枝角度的QTL，位于LG1 連鎖群，解釋14.16%的表型變異。Liu 等[27]利用包含143 份株系的自然群體進行分枝角度的全基因組關(guān)聯(lián)分析，在A2、A3、A7、C3、C5 和C7 染色體上定位到25 個顯著QTL，并在A3 和C3 關(guān)聯(lián)位點附近找到了分枝角度基因LAZY1的同源拷貝。Sun 等[28]利用從國內(nèi)油菜主產(chǎn)區(qū)及歐洲等其它國家收集來的520份自然群體進行全基因組關(guān)聯(lián)分析，定位到56 個QTL，共解釋表型變異的51.1%，并在這些區(qū)域預(yù)測了多個候選基因，包括LAZY1、SGR2（SHOOT GRAVITROPISM 2）、PIN3（PIN-FORMED 3）、TIR1（TRANSPORT INHIBITOR RESPONSE 1）等。Wang 等[29]利用BSA 的方法結(jié)合二代測序技術(shù)對雙親分枝角度差異較大的F2群體進行分枝角度精細(xì)定位，最終在A6染色體上檢測到一個分枝角度主效QTL，解釋17.17%的表型變異，并在該QTL 的置信區(qū)間內(nèi)篩選到生長素合成相關(guān)基因BnaYUCCA6為候選基因。Li 等[30]利用60K SNP 芯片對472 份油菜核心種質(zhì)進行全基因組關(guān)聯(lián)分析，在6 個不同環(huán)境下進行分枝角度表型鑒定，MLM 和MRMLM 模型各檢測出46 個和38 個顯著相關(guān)的位點，分別解釋表型變異的62.2%和66.2%，并根據(jù)同源擬南芥的注釋，首次發(fā)現(xiàn)了定位在油菜中的分枝角度基因TAC1和SGR1等。Shen等[31]利用60K SNP 芯片對含有208 份株系的DH 群體進行基因分型，構(gòu)建遺傳圖譜并進行QTL 定位，共定位到17 個QTL。汪文祥等[32]利用由分枝角度差異顯著親本構(gòu)建的包含163個株系的DH群體，在2 個環(huán)境中共定位到17 個分枝角度QTL，并推測基因VAMP714，TAC1，IAR和ARF16等為油菜分枝角度的候選基因。以上研究初步獲得了一些油菜分枝角度性狀位點信息，但還遠遠不能揭示該性狀的遺傳機制。

前期利用油菜60K SNP 芯片，對包含189 個株系的甘藍型油菜重組自交系（RIL）群體進行基因分型，構(gòu)建了高密度SNP 遺傳圖譜[33]。本研究在此基礎(chǔ)上，對RIL群體的分枝角度性狀進行考察，結(jié)合高密度遺傳圖譜進行QTL 定位與分析，并在穩(wěn)定表達的QTL 區(qū)間內(nèi)篩選候選基因，為深入解析甘藍型油菜分枝角度的遺傳機理，精細(xì)定位和克隆分枝角度相關(guān)基因提供參考。

1 材料與方法

1.1 材料

前期利用父本Holly和母本APL01，構(gòu)建了一個包含189 個株系的RIL 群體并命名為AH。利用油菜60K SNP芯片，構(gòu)建了高密度的遺傳連鎖圖譜，其與甘藍型油菜參考基因組Darmor-bzh 具有良好的共線性[33]。該圖譜包含2755 bin（含11 458 個SNP）和57個SSR，覆蓋甘藍型油菜全基因組19個連鎖群的2027.53 cM，標(biāo)記間的平均距離為0.72 cM[33]。利用該群體，前期對無花瓣性狀[33]、脂肪酸組分[34]、初花期[35]及耐漬性狀[36]進行了QTL 定位及分析。本研究中，AH 群體用于分枝角度的QTL 定位及候選基因分析。

1.2 田間試驗

將RIL 群體及其親本于2015 年9 月至2016 年5月種植于江蘇南京和陜西楊凌兩個試驗點，兩個環(huán)境的實驗分別記錄為15NJ 和15YL。田間播種均采用隨機區(qū)組設(shè)計，每個小區(qū)每個株系種植2行，設(shè)置兩次重復(fù)，行間距40 cm，單株間平均間距20 cm。試驗材料的種植、管理和收獲等均按照當(dāng)?shù)卮筇锷a(chǎn)管理方式進行。

1.3 表型考察

參考孫程明[37]對分枝角度測量的方法，本實驗在油菜成熟后剪取帶有油菜頂部5個分枝的莖段作為單株的分枝角度，并將其裝袋編號。在15NJ 和15YL 兩試驗點，每個株系選取5 株進行考察，每個環(huán)境包括兩次重復(fù)。取樣回來的樣品將其放置在樣品置放板上，利用數(shù)碼相機拍下能清晰呈現(xiàn)主莖與每一個分枝節(jié)點夾角的部分，每株5個分枝共得5張照片，利用計算機Photoshop 軟件中的直線工具和角度工具測量照片中主莖和分枝之間的夾角，并將其記錄到Microsoft Excel文檔中，兩個重復(fù)的平均值作為株系在一個環(huán)境的表型值。利用R 軟件[38]計算其廣義遺傳力。

1.4 QTL定位分析

使用軟件Windows QTL Cartographer 2.5[39]的復(fù)合區(qū)間作圖（composite interval mapping,CIM）法，對2 個環(huán)境下分枝角度進行QTL 的檢測。窗口大小設(shè)為10 cM，選用2 cM 的步長（walking speed），同時設(shè)置5個標(biāo)記作為輔助因子。每個性狀的顯著性水平設(shè)為0.01，進行1000 次排列檢驗（permutation test），以此來判斷是否存在QTL，并將其稱為identified QTL。對于在2 個環(huán)境下都有表達且具有重疊置信區(qū)間（confidence intervals,CIs）的QTL 進行整合并稱為一致性QTL（consensus QTL）；對于在2 個環(huán)境中只檢測到一次的QTL同樣稱為一致性QTL。QTL的整合應(yīng)用軟件BioMercator V4.2[40]元分析（metaanalysis）法實現(xiàn)。參考McCouch 等[41]和Wang 等[42]對QTL的命名的介紹，本研究對QTL命名，以檢測到的QTLidentified QTL 英文縮寫iq 作為開頭，依次加上環(huán)境縮寫（15NJ、15YL）和連鎖群編號（A1-A10、C1-C9）。若在同一個連鎖群上出現(xiàn)了不止一個QTL，就依據(jù)順序?qū)ζ溥M行編號。如QTLiq15NJ. C4-1，表示2015 年南京環(huán)境中檢測到的C4 染色體上的第一個QTL。對于經(jīng)元分析整合后的一致性QTL，其命名法則類似上述identified QTL 命名原則，以一致性QTL 英文縮寫cq 開頭，再依次加上分枝角度性狀縮寫（branch angle,BA）和連鎖群編號。若在同一個連鎖群上出現(xiàn)了不止一個一致性QTL，就依據(jù)順序?qū)ζ溥M行編號。例如QTLcqBA. C4-2表示C4 染色體上分枝角度的第二個一致性QTL。

1.5 候選基因預(yù)測

候選基因的定位已經(jīng)實現(xiàn)了從區(qū)間定位到通過分析基因組數(shù)據(jù)庫信息篩選候選基因的階段[43]。本研究利用Illumina HiSeq 2500 平臺（Illumina，Inc；San Diego，CA，USA），對雙親APL01 和Holly 基因組進行了重測序，測序深度為30×，共獲得約249.95 Mb Clean Reads，74.89 Gb Clean Data，Q30平均達到92.00%。

候選基因的預(yù)測分為3步：首先，利用重測序的結(jié)果，在兩親本的全基因組內(nèi)鑒定變異位點；其次，利用遺傳圖譜上SNP 的探針序列，通過Blast（E 值≤1e-10）將QTL 區(qū)間兩側(cè)的標(biāo)記映射到參考基因組Darmor-bzh 的物理位置[44]，并根據(jù)第一步的結(jié)果，查找QTL 區(qū)間內(nèi)的SNP 位點和InDel 位點及造成非同義突變的SNP 和發(fā)生移碼突變的InDel；第三，以非同義突變的SNP和移碼突變的InDel為目標(biāo)，根據(jù)擬南芥基因的功能注釋篩選分枝角度可能的候選基因。依據(jù)對現(xiàn)有分枝角度相關(guān)文獻的參考，分枝角度候選基因的選擇偏向于功能注釋中關(guān)于向重力性，生長素的生物合成、運輸及代謝相關(guān)基因。測序工作由諾禾致源生物信息科技公司（北京）完成，數(shù)據(jù)分析工作由中國農(nóng)業(yè)科學(xué)院油料作物研究所油料作物分子改良理論與技術(shù)團隊完成。

2 結(jié)果與分析

2.1 RIL群體及其雙親分枝角度表型分析

對RIL群體及親本的分枝角度進行了兩個環(huán)境的考察（圖1、表1）。在15NJ 和15YL 2 個試驗點，APL01 的分枝角度分別為25.2°和32.8°，Holly 的分枝角度分別為33.5°和41.9°，兩親本的分枝角度在兩個環(huán)境中均存在極顯著差異（P＜0.01）（表1）。RIL 群體的分枝角度在兩個環(huán)境中具有廣泛的變異，平均變異幅度各為24.74°～37.76°和31.05°～46.67°。單個環(huán)境的最大值大于高值親本，最小值小于低值親本，表現(xiàn)出超親分離現(xiàn)象（圖1、表1），暗示增效等位基因在雙親中都有分布。兩環(huán)境下的分枝角度表現(xiàn)出連續(xù)性分布，表現(xiàn)出明顯的數(shù)量性狀特征，適合用于QTL檢測（圖1）。R軟件計算出分枝角度在2 個環(huán)境下的廣義遺傳力為86%，說明油菜的分枝角度主要受遺傳因素控制。

表1 兩環(huán)境下親本及RIL群體分枝角度性狀的統(tǒng)計分析Table 1 Statistical analysis of branch angle of parents and RIL populations under 2 environments

圖1 AH RIL群體在兩環(huán)境的分枝角度表型頻率分布Fig.1 Phenotypic frequency distributions of branch angle in AH RIL population under 2 environments

2.2 分枝角度QTL檢測

利用軟件Windows QTL Cartographer 2.5 中的復(fù)合區(qū)間作圖法進行QTL 定位，在15NJ 和15YL 兩環(huán)境中共檢測到11 個identified QTL，分布在A9、C3、C4 和C7 染色體上（表2），單個QTL 解釋表型變異的4.05%～9.60%，加性效應(yīng)范圍為-1.52～1.96。其中有7 個identified QTL 分布在C4 染色體，每個QTL 對表型變異的貢獻率在4.60%～9.60%。利用軟件BioMercator V4.2 將其中5 個具有重疊置信區(qū)間QTL 整合成2 個一致性QTLcqBA. C4-2和cqBA.C4-3（表2），這兩個QTL 均在15NJ 和15YL 被反復(fù)檢測到。其它6 個QTL 只在1 個環(huán)境中檢測到（表2），表現(xiàn)出環(huán)境特異性，說明環(huán)境對于這些QTL 位點影響的程度較大，利用它們對分枝角度進行遺傳改良時，環(huán)境因素不可忽略。經(jīng)過元分析后，本研究共獲得8個分枝角度的一致性QTL，其中7個一致性QTL 的加性效應(yīng)為正值，說明其增效基因來自于母本APL01；另一個一致性QTLcqBA. A9，具有-1.52 的負(fù)加性效應(yīng)，表明其增效基因來自于父本Holly（表2）。

表2 甘藍型油菜分枝角度QTL結(jié)果Table 2 QTL results for branch angle of Brassica napus

2.3 與已報道的油菜分枝角度QTL的比較

油菜中分枝角度的文獻報道較少，對已經(jīng)報道的分枝角度性狀定位結(jié)果進行整理[26～32]，發(fā)現(xiàn)除了C1染色體外，分枝角度位點在其它染色體上都有分布，其中A3，A7 和C3 染色體分布數(shù)目最多，各有10個以上位點，其余染色體各有2～7 個位點。為比較本研究定位的8 個QTL 和前人報道的QTL 是一致，將其連鎖/關(guān)聯(lián)的標(biāo)記/探針序列比對到參考基因組Darmor-bzh 以確定QTL 的物理位置。本研究在A9染色體上定位到一個QTLcqBA.A9，其位于22.15～23.05 Mb 區(qū)間（表2），通過與已研究報道的QTL[28,30]比較，置信區(qū)間沒有重疊（表3）；在C3 染色體上也定位到一個QTLcqBA.C3，其位于4.21～4.87 Mb 區(qū)間（表2），通過與已研究報道的QTL[27,28,30～32]比較，置信區(qū)間沒有重疊（表3）；在C4 染色體上定位到四個QTLcqBA. C4-1（0.082～1.64 Mb），cqBA. C4-2（10.32～14.3 Mb），cqBA.C4-3（44.39～44.46 Mb）和cqBA. C4-4（44.74～46.33 Mb）（表2）。經(jīng)比較，除cqBA. C4-1與汪文祥等[32]的區(qū)間有重疊外，其它3個QTL 與前人的定位結(jié)果不同[28,30,32]（表3）；在C7 染色體上定位到兩個QTLcqBA. C7-1（5.73～5.86 Mb）和cqBA.C7-2（30.72～32.18 Mb）與前人的定位結(jié)果均不同[27,28]（表3）。因此，本研究中定位的8 個QTL，除cqBA.C4-1外，其它7 個QTL 與前人的定位結(jié)果均不同，推測為新的控制油菜分枝夾角的QTL位點。

表3 本研究與前人報道的分枝角度QTL比較Table 3 Comparison of the QTL for branch angle from present and previous reports

2.4 分枝角度候選基因的預(yù)測

QTLcqBA. C4-2和cqBA. C4-3能在2 個環(huán)境中穩(wěn)定表達（表2），表明這2個QTL區(qū)間內(nèi)可能存在控制甘藍型油菜分枝角度的基因。這2個QTL均位于C4 連鎖群上，對應(yīng)物理圖譜上的位置為10.3～14.3 Mb 和44.4～44.5 Mb。根據(jù)Holly 和APL01 的重測序結(jié)果，在QTLcqBA.C4-2區(qū)間內(nèi)共發(fā)現(xiàn)有7226 個SNP 和829 個InDel，其中7226 個SNP 中包含220 個非同義突變和5 個stopgain 突變，829 個InDel 中有15個發(fā)生移碼突變，這些突變位點共對應(yīng)55個基因（附表1，見首頁OSID 碼）；在QTLcqBA. C4-3區(qū)間內(nèi)共發(fā)現(xiàn)有416 個SNP 和65 個InDel，其中416 個SNP中包含50個非同義突變，65個InDel中有1個發(fā)生移碼突變，這些突變位點共對應(yīng)15 個基因（附表2，見首頁OSID 碼）。在此基礎(chǔ)上，以這70 個基因為目標(biāo)，借助擬南芥同源基因的功能注釋，共篩選到4個分枝角度相關(guān)的候選基因，BnaC04g13100D，BnaC04g15900D，BnaC04g16280D和BnaC04g44330D（表4）。

表4 QTL cqBA.C4-2和cqBA.C4-3置信區(qū)間獲得的候選基因Table 4 Candidate genes in QTL confidence interval of cqBA.C4-2 and cqBA.C4-3

在QTLcqBA.C4-2的置信區(qū)間內(nèi)（10 323 625～14 297 777 bp），共篩選到3 個油菜分枝角度相關(guān)基因 ，BnaC04g13100D（AT2G31500,CPK24） 、BnaC04g15900D（AT2G28350,ARF10）和BnaC04g16 280D（AT2G28070,ABCG3）。根據(jù)擬南芥基因的功能注釋，CPK24（BnaC04g13100D）為鈣依賴性蛋白激酶編碼基因，CPK基因廣泛分布在植物體的不同組織，其產(chǎn)物CPK 蛋白通過解讀由各種發(fā)育和環(huán)境刺激觸發(fā)的Ca2+信號來參與調(diào)節(jié)植物的生長[45]。Ca2+作為重要的信號物質(zhì)，在植物的生長發(fā)育中起著關(guān)鍵作用，已有研究表明，Ca2+參與調(diào)控植物的重力反應(yīng)[12]，當(dāng)重力刺激改變時，會導(dǎo)致細(xì)胞內(nèi)Ca2+外流，這些外流Ca2+被CPK 和CaM（鈣調(diào)素）等感知，并將信號傳遞到下游效應(yīng)器，最終導(dǎo)致植物重力響應(yīng)[46]，其中CPK 參與了重力信號的傳遞[47]。推測鈣依賴性蛋白激酶成員基因CPK24（BnaC04g13100D）可能通過參與重力信號的傳遞來調(diào)控分枝角度。Roychoudhry 等[48]研究發(fā)現(xiàn)植物非垂直分枝的生長不僅僅由重力性因子決定，提出了反重力性抵消因子（antigravitropic offset, AGO），其與重力性因子共同調(diào)控側(cè)枝的生長角度。這種反重力抵消機制發(fā)揮作用依賴于生長素轉(zhuǎn)運，擬南芥生長素缺陷突變體，在生長素信號恢復(fù)之后，分枝角度也隨之減小，生長素通過一個復(fù)雜的調(diào)控網(wǎng)絡(luò)參與調(diào)控植物分枝角度的形成[49]。本研究在QTLcqBA.C4-2置信區(qū)間內(nèi)篩選到生長素響應(yīng)因子（auxin response factor,ARF）中的ARF10（BnaC04g15900D），該基因也存在于前人定位的分枝角度區(qū)間內(nèi)[28]?；駻BCG3（BnaC04g16280D）為ABC 轉(zhuǎn)運蛋白（ATP-binding cassette transporter）家族一員，ABC 轉(zhuǎn)運蛋白是植物界普遍存在的一類跨膜運輸?shù)鞍?，其通過消耗ATP產(chǎn)生的能量參與生物活動。Sun 等[28]和汪文祥等[32]在其置信區(qū)間相繼篩選出多個ABC 轉(zhuǎn)運家族為分枝角度基因，推測本研究ABC 轉(zhuǎn)運家族基因ABCG3可能為分枝角度候選基因。在QTLcqBA.C4-3的置信區(qū)間內(nèi)（44 386 986～44 459 092 bp）篩選到1 個油菜分枝角度相關(guān)基因BnaC04g44330D（AT2G3399 0,IQD9），是IQ67結(jié)構(gòu)域家族的成員，該家族成員只與CaM 特異性結(jié)合。CaM 作為Ca2+另一個重要的傳感蛋白，同樣通過傳遞鈣信號，繼而調(diào)控植物的生長發(fā)育。Bürstenbinder 等[50]認(rèn)為多個IQD基因為ARF 的潛在靶標(biāo)，其IQD 蛋白作為鈣信號和生長素信號的樞紐來調(diào)控植物的形態(tài)發(fā)生。因此推測IQD9可能參與調(diào)控分枝角度。

3 討論

分枝角度是決定作物株型的重要因素，對于作物合理密植及機械化收獲有重要意義。迄今為止，在水稻（分蘗角度）、玉米（莖葉夾角）和擬南芥（分枝角度）中已有較多研究，并對相關(guān)基因進行了克隆和鑒定，但在油菜上的研究較少且尚處于定位階段。甘藍型油菜中，除C1 染色體外，其余染色體均有分枝角度QTL 位點的分布[26～32]，其中A3、A7 和C3上的分布數(shù)目最多，各有10個以上位點。本研究在兩個環(huán)境中，共定位到8 個QTL，分布在A9、C3、C4和C7 染色體上。通過和前人研究中報道的分枝角度QTL 的物理位置比較發(fā)現(xiàn)，定位的8 個QTL 除cq-BA. C4-1外，其它7 個QTL 與前人報道的油菜分枝角度QTL 均沒有重疊，推測是新的油菜分枝角度QTL位點。

本研究對在2 個環(huán)境中均穩(wěn)定表達的QTLcq-BA.C4-2和cqBA.C4-3進行了候選基因的預(yù)測，共篩選到4 個油菜分枝角度相關(guān)基因，BnaC04g13100 D（CPK24）、BnaC04g15900D（ARF10）、BnaC04g1628 0D（ABCG3）和BnaC04g44330D（IQD9）。CPK24僅依賴于Ca2+信號來調(diào)節(jié)植物生長發(fā)育[51]，而重力刺激會引起Ca2+的變化[52]，同時有研究表明CPK 參與了重力信號的傳遞[47]。此外，有報道[53]稱CPK 可以通過磷酸化其載體進而影響植物體內(nèi)生長素水平，如定位于質(zhì)膜上的AtCRK5（CPK-related kinase 5）通過磷酸化生長素輸出載體AtPIN2，導(dǎo)致其生長素水平降低，重力反應(yīng)延遲。因此，推測重力因素的改變引起Ca2+信號發(fā)生變化繼而被CPK24識別，通過影響生長素運輸進而控制分枝角度的形成。ARF10為生長素響應(yīng)因子之一，能夠與生長素特異性地結(jié)合，以此影響植物的分枝。Sun 等[28]、Shen等[31]和汪文祥等[32]也在其置信區(qū)間相繼篩選出ARF10和ARF16為分枝角度基因。ABCG3為ABC轉(zhuǎn)運家族蛋白一員，通過消耗ATP 產(chǎn)生的能量來進行跨膜運輸，并主要運輸植物激素。有研究者表明[54]，ABCG36 為一種定位于質(zhì)膜上的ABC 轉(zhuǎn)運蛋白，其突變體pdr8（ABCG36 缺陷型）在生長素前體IBA（吲哚-3-丁酸）轉(zhuǎn)運上表現(xiàn)出缺陷。此外，ABC超家族PDR 轉(zhuǎn)運蛋白直接參與分枝角度的形成[55]。因此，推測ABCG3可能通過影響生長素運輸直接或者間接影響油菜分枝角度。IQD9 屬于IQD（IQ67-domain containing protein）基因家族蛋白，只與CaM特異性結(jié)合。研究表明，IQD 蛋白作為鈣信號和生長素信號的中介來調(diào)控植物的形態(tài)發(fā)生[50]，擬南芥和番茄中已證實了IQD家族基因的過表達導(dǎo)致莖葉等形態(tài)的改變[56,57]。因此，推測IQD9可能通過影響生長素信號和鈣信號來影響油菜分枝角度。

本研究利用油菜基因組序列信息及親本的重測序結(jié)果，從穩(wěn)定表達的QTLcqBA. C4-2和cqBA.C4-3區(qū)段內(nèi)，篩選出4 個油菜分枝角度的候選基因，為進一步理解分枝角度形成的分子機理提供了有價值的信息。