• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于RIL和CSSL群體定位大豆脂肪酸組分QTL

    2021-08-20 01:36:10渠可心韓露謝建國潘文婧張澤鑫辛大偉劉春燕陳慶山齊照明
    中國農(nóng)業(yè)科學(xué) 2021年15期
    關(guān)鍵詞:油分連鎖組分

    渠可心,韓露,謝建國,潘文婧,張澤鑫,辛大偉,劉春燕,陳慶山,齊照明

    基于RIL和CSSL群體定位大豆脂肪酸組分QTL

    渠可心,韓露,謝建國,潘文婧,張澤鑫,辛大偉,劉春燕,陳慶山,齊照明

    東北農(nóng)業(yè)大學(xué)農(nóng)學(xué)院大豆遺傳改良實(shí)驗(yàn)室,哈爾濱 150030

    【】大豆()原產(chǎn)于中國,高品質(zhì)的大豆在食品、飼料、紡織品等多種加工業(yè)中廣泛應(yīng)用,因此,選育高品質(zhì)大豆已成為育種者和生產(chǎn)者的聚焦問題。通過對(duì)大豆脂肪酸各組分進(jìn)行QTL定位及候選基因的篩選,為大豆品質(zhì)改良奠定分子基礎(chǔ)。以美國大豆品種Charleston和東農(nóng)594為親本構(gòu)建重組自交系(RILs)、以栽培大豆綏農(nóng)14與野生大豆ZYD00006為親本構(gòu)建染色體片段代換系(CSSLs)為試驗(yàn)材料。利用氣相色譜法測(cè)定2個(gè)群體的脂肪酸含量,根據(jù)東北農(nóng)業(yè)大學(xué)農(nóng)學(xué)院大豆遺傳改良實(shí)驗(yàn)室已構(gòu)建的遺傳圖譜,通過Windows QTL Cartographer 2.5和ICIMapping軟件對(duì)2017—2018年RIL群體與CSSL群體中的大豆脂肪酸組分進(jìn)行QTL定位研究,并對(duì)所獲得的QTL置信區(qū)間進(jìn)行候選基因的挖掘。2017—2018年,RIL群體和CSSL群體分別定位到34和20個(gè)與脂肪酸組分相關(guān)的QTL,分布在除B2、C1、G、H、J、M和O以外的13個(gè)連鎖群上。比較2個(gè)群體的QTL定位結(jié)果,發(fā)現(xiàn)在2個(gè)群體中重復(fù)檢測(cè)到10對(duì)QTL,其中,分布在A1、C2、D1a、F、K和N連鎖群上的QTL與多種脂肪酸含量相關(guān),在A1連鎖群上檢測(cè)到亞油酸和油分含量重疊的QTL;在C2連鎖群上檢測(cè)到硬脂酸和油分含量重疊的QTL;在D1a連鎖群上檢測(cè)到硬脂酸和油分含量重疊的QTL;在F連鎖群上檢測(cè)到棕櫚酸、硬脂酸和油分含量重疊的QTL;在K連鎖群上檢測(cè)到亞油酸和亞麻酸含量重疊的QTL;在N連鎖群上檢測(cè)到棕櫚酸和油分含量重疊的QTL、油酸和亞油酸含量重疊的QTL。對(duì)QTL定位獲得的所有置信區(qū)間進(jìn)行候選基因的挖掘,從基因注釋數(shù)據(jù)集中共篩選出485個(gè)候選基因,其中,271個(gè)候選基因具有GO注釋,進(jìn)一步進(jìn)行GO富集數(shù)據(jù)分析,共有15個(gè)候選基因與脂肪酸相關(guān)。主要通過編碼植物?;??;d體蛋白(ACP)硫酯酶、脂肪酸去飽和酶、磷脂酶D1、脂肪酸-羥化酶、丙酮酸激酶和參與?;o酶A生物合成、調(diào)節(jié)脂肪酸鏈的延伸,從而影響脂肪酸的合成。共檢測(cè)到54個(gè)與大豆脂肪酸各組分相關(guān)的QTL,在2個(gè)群體中重復(fù)檢測(cè)到10對(duì)QTL,對(duì)QTL定位獲得的置信區(qū)間進(jìn)行候選基因的篩選,共有15個(gè)候選基因與脂肪酸相關(guān)。這些穩(wěn)定的脂肪酸相關(guān)的QTL和脂肪酸相關(guān)的候選基因可用于大豆脂肪酸改良的分子標(biāo)記輔助選擇。

    大豆;導(dǎo)入系;氣相色譜法;QTL定位;基因挖掘;富集分析

    0 引言

    【研究意義】大豆()原產(chǎn)于中國,是一種重要的經(jīng)濟(jì)作物,占世界油籽產(chǎn)量的61%,蛋白質(zhì)食物消費(fèi)量的68%[1]。大豆是植物油的主要供給原料,油分含量占籽粒干重的18%—23%,主要包含5種重要脂肪酸組分:棕櫚酸(palmitic acid)、硬脂酸(stearic acid)、油酸(oleic acid)、亞油酸(linoleic acid)和亞麻酸(linolenic acid)。大豆已經(jīng)成為人們?nèi)粘I钪猩攀诚牟豢扇鄙俚氖称?,與國民糧食安全的可持續(xù)發(fā)展緊密相連[2-4]。目前,中國處于大豆進(jìn)口最多的國家,而且高品質(zhì)的大豆在食品、飼料、紡織品等多種加工工業(yè)中的需求量日益劇增,出現(xiàn)了較為明顯的供需矛盾,鑒于這一問題,選育高品質(zhì)大豆新品種已成為育種者和生產(chǎn)者的聚焦問題[5-6]?!厩叭搜芯窟M(jìn)展】大豆脂肪酸中,亞油酸含量最高,油酸次之,亞麻酸、棕櫚酸、硬脂酸含量相對(duì)較低[7],其中,油酸已被證明可以降低血液中的膽固醇含量[4],亞油酸已被證明可以防止動(dòng)脈硬化的作用[8],亞麻酸具有降血壓、健腦抗衰老等作用[9]。因此,在衡量大豆油品質(zhì)的優(yōu)劣中,脂肪酸各組分的比例是一個(gè)非常重要的指標(biāo)。大豆脂肪酸組分含量的檢測(cè)方法主要包括氣相色譜法(gas chromatography,GC)[10-11]、近紅外反射光譜分析法(near infrared reflectance spectrometry,NIRS)[12-13]和高效液相色譜法(high performance liquid chromatography,HPLC)[14-15],而氣相色譜法是應(yīng)用最為廣泛的[16]。氣相色譜法是色譜法中的一種,主要是利用混合物各物質(zhì)之間的沸點(diǎn)、極性及吸附性的差異性來實(shí)現(xiàn)分離,目前,現(xiàn)代分離分析主要應(yīng)用該技術(shù)[17]。大豆脂肪酸的5種組分是由多基因控制的數(shù)量性狀,容易受不同環(huán)境條件的影響[18]。國內(nèi)外研究學(xué)者鑒定影響大豆籽粒貯藏物含量的QTL位點(diǎn),主要通過分子標(biāo)記技術(shù)結(jié)合區(qū)間作圖(interval mapping,IM)[19]、復(fù)合區(qū)間作圖(composite interval mapping,CIM)[20]、多區(qū)間作圖(multiple interval mapping,MIM)[21]和完備區(qū)間作圖(inclusive composite interval mapping,ICIM)[22]等常用的分析方法來進(jìn)行。QTL定位的實(shí)質(zhì)是分析分子標(biāo)記與目標(biāo)性狀之間的連鎖關(guān)系[23],隨著分子生物技術(shù)的飛速發(fā)展,國內(nèi)外研究學(xué)者已開發(fā)出適應(yīng)各種情況的QTL定位方法,大大地解決了分子輔助育種過程中所遇到的難題,為分子輔助育種的高效快速發(fā)展奠定了堅(jiān)實(shí)基礎(chǔ)[24]。目前,國內(nèi)外的研究學(xué)者對(duì)大豆脂肪酸含量進(jìn)行了大量的QTL定位研究[25-26]。Li等[27]共定位到26個(gè)與脂肪酸含量相關(guān)的穩(wěn)定QTL;XIA等[28]對(duì)2個(gè)群體在不同環(huán)境下進(jìn)行大豆脂肪酸各組分的QTL定位分析,共定位到16個(gè)相關(guān)的QTL;FAN等[29]共定位到52個(gè)相關(guān)的QTL;盛英華等[30]對(duì)2個(gè)群體的脂肪酸含量進(jìn)行QTL定位分析,共定位到57個(gè)相關(guān)的QTL?!颈狙芯壳腥朦c(diǎn)】目前,雖然與大豆脂肪酸相關(guān)的QTL已被大量報(bào)道,但是可以在不同環(huán)境或不同遺傳背景中被重復(fù)檢測(cè)到的QTL較少[31]?!緮M解決的關(guān)鍵問題】本研究于2017年和2018年測(cè)定190份染色體段代換系(chromosome segment substitution lines,CSSLs)群體及親本和147份重組自交系(recombinant inbred lines,RILs)群體及親本的脂肪酸各組分含量,并對(duì)大豆脂肪酸各組分進(jìn)行QTL定位,以期獲得能在不同遺傳背景下被重復(fù)檢測(cè)到的脂肪酸QTL。對(duì)檢測(cè)到QTL的定位區(qū)段進(jìn)行候選基因的挖掘,并選擇與脂肪酸相關(guān)的基因,為后續(xù)與脂肪酸各組分相關(guān)基因的精細(xì)定位及基因功能的研究奠定基礎(chǔ),并對(duì)大豆品質(zhì)改良方面具有一定的指導(dǎo)意義。

    1 材料與方法

    1.1 試驗(yàn)材料

    選用東北農(nóng)業(yè)大學(xué)農(nóng)學(xué)院大豆遺傳改良實(shí)驗(yàn)室歷經(jīng)12年構(gòu)建,以栽培大豆綏農(nóng)14與野生大豆ZYD00006為親本,并以綏農(nóng)14為輪回親本,采用分子輔助選擇段構(gòu)建的一套全基因組染色體片段代換系(CSSLs)的190份材料。以美國半矮稈大豆品種Charleston(♂,中國農(nóng)業(yè)科學(xué)院作物科學(xué)研究所提供)和高蛋白大豆品系東農(nóng)594(♀,東北農(nóng)業(yè)大學(xué)大豆研究所提供)親本構(gòu)建的重組自交系(RILs)群體的147份材料。2017—2018年播種在具有適當(dāng)土壤含水量(15%—20%)的中國哈爾濱(45.75°N,126.53°E)的向陽農(nóng)場(chǎng)試驗(yàn)田,采用完全隨機(jī)區(qū)組設(shè)計(jì),行長5 m,株距5 cm,壟距60 cm,3次重復(fù),田間管理遵循正常的生產(chǎn)實(shí)踐。

    1.2 脂肪酸的測(cè)定方法

    采用Agilent7890B氣相色譜儀測(cè)定大豆脂肪酸的5種組分,每份材料隨機(jī)選取3粒籽粒,于65℃烘箱中烘干(12—14 h),研磨成粉末,稱取5 mg左右樣品,置于儲(chǔ)存管中。在儲(chǔ)存管中分別加入1 mL 2.5%的濃硫酸甲醇溶液、5 μL BHT(2,6-二叔丁基-4-甲基苯酚50 mg·mL-1,甲醇溶解)還原劑和50 μL的內(nèi)標(biāo)十七烷酸(10 mg·mL-1,乙酸乙酯溶解),混勻后85℃水浴1.5 h,冷卻至室溫,加入160 μL 9%的NaCl溶液和700 μL正己烷,渦旋、震蕩,4 000 r/min常溫離心10 min,吸取上清液400—500 μL于新儲(chǔ)存管中,置于通風(fēng)櫥中吹干、過夜。最后加400 μL乙酸乙酯并迅速用氣相色譜儀進(jìn)行脂肪酸的測(cè)定,每個(gè)樣品測(cè)定3次,取平均值作為該脂肪酸組分的含量。

    Agilent7890B氣相色譜儀所用色譜柱型號(hào)為:30 m×320 μm×0.25 μm。載氣:氮?dú)?0 mL·min-1,氫氣60 mL·min-1,空氣450 mL·min-1。進(jìn)樣量:1 μL,分流進(jìn)樣方式,分流比為10﹕1;進(jìn)樣口溫度170℃。反應(yīng)程序?yàn)?80℃保持1 min,以25℃·min-1的速率上升至250℃保持7 min。

    氣相色譜分析過程中,首先試樣在汽化室內(nèi)被汽化,然后被惰性氣體(即流動(dòng)相)帶入色譜柱,柱內(nèi)含有固定相,由于試樣中各組分之間的氣相和固定相之間的分子作用力不同,使得在色譜柱中各個(gè)組分的運(yùn)行速度不同,從而產(chǎn)生了不同的峰值,最后分析色譜圖中各種組分的峰值,計(jì)算出不同組分的相對(duì)含量[21]。

    1.3 脂肪酸各組分含量計(jì)算

    氣相色譜儀的峰圖是以標(biāo)準(zhǔn)溶液濃度為橫坐標(biāo),色譜峰面積為縱坐標(biāo),從而繪制不同濃度梯度的標(biāo)準(zhǔn)曲線[16]。根據(jù)5種脂肪酸標(biāo)準(zhǔn)樣品的保留時(shí)間和峰面積進(jìn)行定性定量分析各個(gè)脂肪酸組分,以十七烷酸的峰值為內(nèi)標(biāo),計(jì)算5種脂肪酸的絕對(duì)含量(mg·g-1)和相對(duì)含量(%),油分含量為5種脂肪酸的總和。

    相對(duì)量計(jì)算公式:。

    1.4 脂肪酸含量表型及QTL分析

    采用Microsoft Excel 2010軟件對(duì)2個(gè)群體親本和后代單株脂肪酸各組分進(jìn)行基本統(tǒng)計(jì)量分析,包括最小值、最大值、平均值、變異系數(shù)、標(biāo)準(zhǔn)差、變異幅度、偏度和峰度。采用SPSS軟件對(duì)脂肪酸各組分進(jìn)行相關(guān)性分析。

    RIL群體的遺傳圖譜由QI等[32]構(gòu)建完成,圖譜總長度為2 655.68 cm,包括5 308個(gè)特異長度擴(kuò)增片段標(biāo)記,覆蓋20個(gè)連鎖群(linkage group,LGs),標(biāo)記間平均距離為0.5 cm。CSSL群體的遺傳圖譜由XIN等[33]構(gòu)建完成。使用Win-QTLCart2.5軟件通過CIM對(duì)RILs群體進(jìn)行QTL定位,使用1 000次排列檢驗(yàn)分析表型數(shù)據(jù),并設(shè)置顯著性=0.05,取LOD=2.5為閾值判斷QTL的存在。采用ICIM Mapping軟件中ICIM對(duì)CSSLs群體進(jìn)行與脂肪酸相關(guān)的QTL定位分析,使用1 000次排列檢驗(yàn)分析表型數(shù)據(jù),并設(shè)置顯著性=0.05。QTL命名基于前研究者的方法,簡(jiǎn)而言之,QTL名稱的構(gòu)造如下:+性狀名稱+連鎖群或連鎖群數(shù)+“-”+QTL數(shù)目[34]。

    2 結(jié)果

    2.1 2個(gè)群體的脂肪酸含量和油分總量分布情況

    運(yùn)用氣象色譜儀對(duì)2017—2018年的RILs和CSSLs群體籽粒中的5種脂肪酸含量和油分總量進(jìn)行測(cè)定,并計(jì)算5種脂肪酸和油分含量的頻率分布,以分析大豆籽粒在不同年份、不同群體脂肪酸積累的特點(diǎn)。結(jié)果表明,2017—2018年2個(gè)群體的脂肪酸各組分都呈現(xiàn)連續(xù)正態(tài)分布(圖1和圖2),表明脂肪酸是以數(shù)量方式遺傳的,并適用于脂肪酸各組分含量的QTL定位。油分含量變化范圍為19.66%—24.26%,符合正態(tài)分布(圖3),適用于油分含量的QTL定位,此外,在后代中也觀察到明顯的超親分離現(xiàn)象[27],表明父母雙方對(duì)脂肪酸組成都有貢獻(xiàn)。

    X軸代表脂肪酸各組分所占總油分的百分比,Y軸代表每個(gè)脂肪酸各組分所占總油分百分比的株系數(shù),黑色實(shí)線表示RILs群體的正常曲線。2018年RILs群體的5種脂肪酸含量頻率分布與2017年的結(jié)果相同,呈現(xiàn)連續(xù)正態(tài)分布,所以在這里僅展示了2017年的結(jié)果。下同The X-axis represents the percentage of each fatty acid component to the total oil content, the Y-axis represents the plant coefficient of the percentage of each fatty acid component in the total oil content, and the black solid line represents the normal curve of the RILs population. The frequency distribution of the five fatty acid content of the RILs population in 2018 is the same as the 2017 results, showing continuous normal distribution, so only the 2017 results are shown here. The same as below

    黑色實(shí)線表示CSSLs群體的正常曲線

    2.2 RILs和CSSLs群體脂肪酸含量表型分析

    對(duì)2017—2018年2個(gè)群體親本及后代單株的脂肪酸各組分含量和油分總量進(jìn)行基本統(tǒng)計(jì)分析(表1)。2年的RILs、CSSLs群體母本的棕櫚酸、硬脂酸和油酸含量均低于父本,而亞油酸和亞麻酸含量均高于父本,RILs群體2個(gè)親本的油酸含量差異較大,而CSSLs群體2個(gè)親本的油酸、亞油酸含量差異較大。后代群體的5種脂肪酸含量中除硬脂酸的平均值比親本高,其余性狀的平均值均位于兩親本之間,并且具有雙向超親的現(xiàn)象。這5種脂肪酸在RILs和CSSLs 2個(gè)群體中均存在較大變異,變異系數(shù)為1.13%—6.99%,其中,亞油酸的變異系數(shù)(coefficient of variation,CV)最小,為1.76%—2.31%;硬脂酸的變異系數(shù)(CV)最大,為3.91%—6.99%,2017年RILs群體的變異系數(shù)達(dá)到6.99%;油酸的變異系數(shù)僅次于硬脂酸,表明硬脂酸和油酸在后代群體具有較大差異,提高其在大豆5種脂肪酸組分中所占的比例具有一定意義。同一群體中5種脂肪酸各組分最大值和最小值差異較明顯,說明在后代分離群體中的表現(xiàn)分離程度較大,且變異連續(xù),是屬于多基因控制的數(shù)量性狀,受基因和環(huán)境的共同影響。

    黑色箭頭代表RILs群體和CSSLs群體的父本Charleston和野生豆(ZYD00006)油分總量所在位置,紅色箭頭代表RILs群體和CSSLs群體的母本東農(nóng)594和綏農(nóng)14油分總量所在位置

    對(duì)2017—2018年RILs和CSSLs群體的5種脂肪酸含量進(jìn)行相關(guān)性分析(表2和表3),表明2個(gè)群體中5種脂肪酸含量相關(guān)性的趨勢(shì)一致,其中,油酸與棕櫚酸、硬脂酸、亞油酸和亞麻酸呈顯著負(fù)相關(guān);棕櫚酸與硬脂酸、亞油酸和亞麻酸呈顯著正相關(guān),亞油酸與亞麻酸呈顯著正相關(guān)。其中,油酸和亞油酸的相關(guān)系數(shù)最大,2018年RILs群體中達(dá)到-0.939,2017年CSSLS群體中達(dá)到-0.991,說明可以通過對(duì)油酸或亞油酸的選擇,從而達(dá)到對(duì)另一種成分的選擇。

    2.3 大豆脂肪酸含量QTL定位

    基于所構(gòu)建的大豆遺傳連鎖圖譜,采用Win-QTLCart2.5和ICIM Mapping軟件分別對(duì)RILs和CSSLs群體進(jìn)行QTL定位(表4和表5),2017—2018年的2個(gè)群體共定位到54個(gè)與脂肪酸組分相關(guān)的QTL,RILs和CSSLs群體分別定位到34和20個(gè)QTL,分布在除B2、C1、G、H、J、M和O以外的13個(gè)連鎖群上,其中,以N連鎖群上最多,有10個(gè)QTL;其次是A1連鎖群上,有8個(gè)QTL;而A2和L連鎖群上均只有1個(gè)。

    檢測(cè)到與棕櫚酸相關(guān)的QTL共有12個(gè),RILs群體檢測(cè)到10個(gè)QTL分別位于A1、C2、D1b、E、F和N等6個(gè)連鎖群上,貢獻(xiàn)率為6.47%—13.50%,LOD值為2.64—5.60。其中,被檢測(cè)到位于E連鎖群上,QTL效應(yīng)值較大,貢獻(xiàn)率為13.50%。CSSLs群體檢測(cè)到2個(gè)QTL都位于N連鎖群上,表型貢獻(xiàn)率為4.62和7.75,LOD值為2.86和6.87。

    檢測(cè)到6個(gè)與硬脂酸相關(guān)的QTL,其中RILs群體有4個(gè)QTL,分布在D1b、D2、F和N等4個(gè)連鎖群上,貢獻(xiàn)率為6.90%—11.11%,LOD值為2.84—4.51。其中,被檢測(cè)到位于F連鎖群上,QTL效應(yīng)值較大,貢獻(xiàn)率為11.11%;CSSLs群體檢測(cè)到2個(gè)QTL分別位于C2和D1b連鎖群上,表型貢獻(xiàn)率為3.21和5.00,LOD值為1.78和2.95。

    表1 2個(gè)群體親本和后代單株脂肪酸含量的基本統(tǒng)計(jì)

    PA:棕櫚酸;SA:硬脂酸;OA:油酸;LA:亞油酸;LNA:亞麻酸。下同

    PA: Palmitic acid; SA: Stearic acid; OA: Oleic acid; LA: Linoleic acid; LNA: Linolenic acid. The same as below

    表2 RILs群體脂肪酸含量的相關(guān)分析

    對(duì)角線下方為2017年RILs和CSSLs群體各脂肪酸之間的相關(guān)系數(shù),上方為2018年RILs和CSSLs群體各脂肪酸之間的相關(guān)系數(shù);**:表示在<0.01水平上顯著相關(guān)。下同

    Below the diagonal line is the correlation coefficient between the fatty acids in the RILs and CSSLs populations in 2017, and above is the correlation coefficient between the fatty acids in the RILs and CSSLs populations in 2018; **: indicates a significant correlation at the<0.01 level. The same as below

    表3 CSSLs群體脂肪酸含量的相關(guān)分析

    表4 RIL群體脂肪酸含量QTL定位

    表5 CSSL群體脂肪酸含量QTL定位

    檢測(cè)到與油酸相關(guān)的QTL共有10個(gè),其中RILs群體檢測(cè)到7個(gè)QTL分別位于A1、D1a、D2、和N等4個(gè)連鎖群上,貢獻(xiàn)率為5.64%—9.43%,LOD值為2.65—3.94;CSSLs群體檢測(cè)到3個(gè)QTL分別位于D2、E和F連鎖群上,表型貢獻(xiàn)率為4.49%、5.22%和8.80%,LOD值為2.75、3.98和9.60。

    檢測(cè)到8個(gè)與亞油酸相關(guān)的QTL,其中RILs群體有6個(gè)QTL,分布在A1、I、K、L和N等5個(gè)連鎖群上,貢獻(xiàn)率為6.40%—8.34%,LOD值為2.54—3.26;CSSLs群體檢測(cè)到2個(gè)QTL都位于A1連鎖群上,表型貢獻(xiàn)率為22.09%和4.28%,LOD值為11.65和2.59。其中,被檢測(cè)到位于A1連鎖群上,QTL效應(yīng)值較大,貢獻(xiàn)率為22.09%。

    檢測(cè)到11個(gè)與亞麻酸相關(guān)的QTL,其中RILs群體有5個(gè)QTL,分布在B1、I、K和N等4個(gè)連鎖群上,貢獻(xiàn)率為7.74%—13.49%,LOD值為3.20—5.46,其中被檢測(cè)到位于K連鎖群上,QTL效應(yīng)值較大,貢獻(xiàn)率為13.49%;CSSLs群體檢測(cè)到6個(gè)QTL分別位于A1、A2、E、F和K等5個(gè)連鎖群上,表型貢獻(xiàn)率為2.92%—5.86%,LOD值為2.31—5.51。

    檢測(cè)到7個(gè)與油分相關(guān)的QTL,其中RILs群體有2個(gè)QTL,分布在D1b和F等2個(gè)連鎖群上,貢獻(xiàn)率為6.95%和10.04%,LOD值為2.67和3.79,其中被檢測(cè)到位于F連鎖群上,QTL效應(yīng)值較大,貢獻(xiàn)率為10.04%;CSSLs群體檢測(cè)到5個(gè)QTL分別位于A1、C2、D1a、F和N等5個(gè)連鎖群上,表型貢獻(xiàn)率為3.02%—8.14%,LOD值為1.84—3.67。

    2.4 RIL和CSSL群體QTL分析結(jié)果比較

    2017—2018年RIL和CSSL群體共定位到54個(gè)QTL,分布在13個(gè)連鎖群上,在同一群體中被重復(fù)檢測(cè)到的QTL有10對(duì),分別位于A1、B1、C2、D1a、D2、F、K和N連鎖群上(表4、表5和圖4)。在RILs群體中重復(fù)檢測(cè)到的QTL有6對(duì),其中,在A1 和D2連鎖群上有2對(duì)QTL,均連續(xù)2年被重復(fù)檢測(cè)到與油酸相關(guān),說明這是2對(duì)是油酸相關(guān)的穩(wěn)定QTL;在B1連鎖群上的1對(duì)QTL同樣連續(xù)2年被重復(fù)檢測(cè)到與亞麻酸相關(guān),分別解釋了7.74%和9.89%的貢獻(xiàn)率,說明該QTL是與亞麻酸相關(guān)的穩(wěn)定QTL;另外在F連鎖群上定位到的3個(gè)QTL,與棕櫚酸、硬脂酸和油分均相關(guān),定位在同一區(qū)間內(nèi),且位置相同,說明為棕櫚酸、硬脂酸和油分含量重疊的QTL,可能調(diào)控多種脂肪酸的合成;定位在連鎖群K上的1對(duì)QTL,檢測(cè)與亞油酸、亞麻酸相關(guān),說明該QTL可能會(huì)調(diào)控多種脂肪酸的合成;定位在連鎖群N上的1對(duì)QTL,檢測(cè)為油酸和亞油酸含量重疊的QTL,說明該QTL可能會(huì)調(diào)控多種脂肪酸的合成。

    :2017年RIL群體定位到的QTL;:2018年RIL群體定位到的QTL;:2017年CSSL群體定位到的QTL;:2018年CSSL群體定位到的QTL

    在CSSLs群體中被重復(fù)檢測(cè)的有4對(duì)。定位在N連鎖群上的3個(gè)QTL,其中和連續(xù)2年被重復(fù)檢測(cè)到并與棕櫚酸相關(guān),說明可能是棕櫚酸相關(guān)的穩(wěn)定QTL,同時(shí)又與油分相關(guān)的為同一個(gè)QTL,說明其為棕櫚酸和油分含量重疊的QTL,并且能夠調(diào)控多種脂肪酸的合成,值得關(guān)注;定位在C2連鎖群上的1對(duì)QTL和定位在D1a連鎖群上的1對(duì)QTL,均為硬脂酸和油分含量重疊的QTL,說明可能是調(diào)控多種脂肪酸合成的QTL;在A1連鎖群上有3個(gè)QTL,其中和連續(xù)2年被檢測(cè)到并與亞油酸相關(guān),并且在2017年表型貢獻(xiàn)率達(dá)到了22.09%,說明可能是亞油酸相關(guān)的穩(wěn)定QTL,同時(shí)又與油分相關(guān)的為同一個(gè)QTL,說明其為亞油酸和油分含量重疊的QTL,可能是控制大豆脂肪酸含量的重要位點(diǎn)。

    對(duì)不同群體定位到的脂肪酸組分的QTL進(jìn)行比較發(fā)現(xiàn),N連鎖群上檢測(cè)到的QTL最多,并且5種脂肪酸的QTL均被檢測(cè)到,說明可能較多的與脂肪酸相關(guān)的QTL定位在該連鎖群上。RILs群體的與CSSLs群體的被檢測(cè)定位在A1連鎖群上,且區(qū)間較近,可能為同一個(gè)QTL,可能是一個(gè)在不同群體中能夠穩(wěn)定遺傳的QTL,定位結(jié)果需進(jìn)一步驗(yàn)證。

    2.5 候選基因挖掘及富集分析

    通過對(duì)RILs和CSSLs群體的5種脂肪酸和油分總量進(jìn)行的QTL定位,定位到54個(gè)QTL,對(duì)QTL定位所獲得的區(qū)間進(jìn)行候選基因的挖掘及GO富集數(shù)據(jù)分析。從基因注釋數(shù)據(jù)集[35]中發(fā)現(xiàn)有9個(gè)區(qū)間不含基因,共篩選出485個(gè)候選基因,對(duì)具有GO注釋的候選基因進(jìn)行進(jìn)一步的GO富集數(shù)據(jù)分析,如表6所示,有15個(gè)基因與大豆脂肪酸的形成直接或間接相關(guān)。這15個(gè)候選基因都在脂肪酸的合成或分解過程中發(fā)揮重要作用。其中候選基因與酰基輔酶a/?;D(zhuǎn)移酶有關(guān),能夠調(diào)節(jié)脂肪酸鏈的延伸,影響脂肪酸的合成[36-37];編碼植物?;?酰基載體蛋白(acyl carrier protein,ACP)硫酯酶,ACP的特異性決定了大多數(shù)植物組織中飽和脂肪酸的水平[36,38],候選基因?qū)χ舅岬男纬煞浅V匾@兩個(gè)基因都與脂肪酸合成的調(diào)控有關(guān)。候選基因編碼磷脂酶D1,已被證實(shí)磷脂酶D在甘油磷脂的代謝中起重要作用[39];候選基因與二羥基丙酮激酶/FAD還原酶有關(guān),是甘油醛生產(chǎn)的重要介質(zhì),與脂類的甘油代謝過程有關(guān)[40],這兩個(gè)基因主要影響脂質(zhì)代謝,從而調(diào)節(jié)植物組織中的脂肪酸含量。候選基因、和都與乙醇乙?;D(zhuǎn)移酶有關(guān),是脂肪酸合成途徑中的重要催化劑[41]。候選基因編碼一種脂肪酸-羥化酶,通過介導(dǎo)長鏈脂肪酸轉(zhuǎn)化為ω-羥基脂肪酸,影響角質(zhì)、軟木脂和蠟的生物合成[42]。是丙酮酸激酶家族蛋白成員,丙酮酸激酶是糖酵解途徑中的一種限速酶,能夠催化磷酸烯醇式丙酮酸生成ATP和丙酮酸的相關(guān)反應(yīng),而丙酮酸又是三羧酸循環(huán)的前體物質(zhì),參與物質(zhì)代謝和氨基酸合成等多種生理生化反應(yīng)[43],有研究發(fā)現(xiàn),丙酮酸激酶有調(diào)節(jié)與控制光合作用的產(chǎn)物轉(zhuǎn)化為油脂的作用[44]、是催化光合產(chǎn)物生成油脂的關(guān)鍵酶[45]。另外編碼脂肪酸去飽和酶,而與酮酸脫羧酶有關(guān),、、和4個(gè)候選基因都與脂肪酶/?;饷赣嘘P(guān),這幾個(gè)基因編碼的蛋白質(zhì)在幾個(gè)方面影響脂肪酸含量。

    3 討論

    本研究進(jìn)行基因精細(xì)定位所利用的染色體片段代換系(CSSLs)群體,已成功應(yīng)用在包括水稻、玉米等作物中[46-48]。CIM和ICIM是QTL定位的重要方法,廣泛應(yīng)用于大豆性狀的QTL定位中,CIM能夠克服區(qū)間作圖法連鎖背景的影響;而ICIM能夠解決QTL算法上的缺陷,充分利用所有標(biāo)記,本研究同時(shí)運(yùn)用了這兩種QTL定位方法,以期能獲得較為準(zhǔn)確的QTL[19-21]。對(duì)2個(gè)群體的5脂肪酸的相關(guān)性分析結(jié)果:棕櫚酸與硬脂酸、亞油酸、亞麻酸呈顯著正相關(guān),而油酸與其余4種組分均呈顯著負(fù)相關(guān);與于福寬[49]研究結(jié)果相同,說明5種脂肪酸之間的相關(guān)性穩(wěn)定,可以在以后的大豆品質(zhì)改良研究中作為參考。

    表6 基因注釋候選基因

    本研究于2017—2018年在2個(gè)群體中共定位到54個(gè)QTL,與目前報(bào)道的脂肪酸相關(guān)的QTL結(jié)果比較,發(fā)現(xiàn)2個(gè)群體可能共有5個(gè)QTL區(qū)間與前人結(jié)果相近,其中,在D1a連鎖群上定位到的區(qū)間與盛英華等[30]相近,但表型解釋率較低,為6.95%;在D2連鎖群上定位到的區(qū)間與盛英華等[30]、朱明月[50]的定位結(jié)果相近,這個(gè)區(qū)間存在2個(gè)與油酸相關(guān)的QTL,其中的表型解釋率較高,為9.4%,可能是一個(gè)穩(wěn)定的QTL區(qū)間且與油酸相關(guān);在E連鎖群上定位到的區(qū)間與DIERS等[51]的相同,且位置相近可能是同一個(gè)QTL,其中和的表型解釋率較高,分別為13.5%和9.4%,可能是一個(gè)穩(wěn)定的QTL,并且該區(qū)間與多種脂肪酸相關(guān);在L連鎖群上定位的區(qū)間與HYTEN等[52]的相近,但表型解釋率較低,為7.26%;在I連鎖群上定位到的區(qū)間在FAN等[29]定位的區(qū)間內(nèi),推測(cè)有多種脂肪酸相關(guān)的QTL分布在這一區(qū)間內(nèi)。本研究定位到的QTL從貢獻(xiàn)率上來說,大部分是一些微效的QTL,可能是由于雙親與后代群體的脂肪酸各組分和油分的變異系數(shù)較小,最大變異系數(shù)僅為6.99%(表1),微效基因較多與雙親和后代群體的表型差異較小有關(guān)。除了部分QTL與已有的定位結(jié)果相近,同時(shí)發(fā)現(xiàn)了一些新的QTL區(qū)間,并且2個(gè)群體中共有10對(duì)QTL被多次重復(fù)檢測(cè)到,表現(xiàn)穩(wěn)定;同時(shí)還有部分的QTL區(qū)間存在重疊,并與多種脂肪酸組分相關(guān),可能是由于一個(gè)基因可以影響生物多種性狀的表型,也可能是由于脂肪酸各組分QTL之間的緊密連鎖[30];但能在2個(gè)群體重復(fù)檢測(cè)到的QTL較少,僅有1對(duì)??傊?,本研究檢測(cè)到的部分QTL與前人定位的QTL區(qū)間一致,進(jìn)一步驗(yàn)證了本研究的可靠性;另一方面能在不同年份、不同群體中被重復(fù)檢測(cè)到的QTL區(qū)間和新發(fā)現(xiàn)的區(qū)間均是與大豆脂肪酸含量相關(guān)的重要區(qū)段,具有更高的利用價(jià)值。

    進(jìn)一步對(duì)所定位到的這些區(qū)間進(jìn)行候選基因的挖掘及富集分析,一共有15個(gè)候選基因與脂肪酸相關(guān),這些候選基因主要分布在第2、3、5、6和17染色體上,其中,基因被許多學(xué)者研究報(bào)道[36],編碼植物?;??;d體蛋白(ACP)硫酯酶,說明該基因?qū)χ舅岬男纬煞浅V匾℅O:0006633)[36];是丙酮酸激酶家族蛋白成員,其家族成員已在甘藍(lán)和擬南芥中被報(bào)道與脂肪酸代謝相關(guān)[53];與酰基輔酶a/?;D(zhuǎn)移酶有關(guān),輔酶A已經(jīng)被多次報(bào)道參與飽和脂肪酸的生物合成[36-37,53]。此外,本研究在其他QTL的置信區(qū)間內(nèi)發(fā)現(xiàn)許多未知功能基因,可能存在與大豆脂肪酸相關(guān)的基因,有待進(jìn)一步探索。

    4 結(jié)論

    2017—2018年在RILs群體和CSSLs群體中分別定位到34和20個(gè)與脂肪酸含量相關(guān)的QTL。其中,有10對(duì)QTL在不同年份的2個(gè)群體中被重復(fù)檢測(cè)到,在RILs群體重復(fù)檢測(cè)到6對(duì),在CSSLs群體中重復(fù)檢測(cè)到4對(duì),并且與多種脂肪酸相關(guān)。檢測(cè)到15個(gè)候選基因與脂肪酸相關(guān)。

    [1] QI Z M, ZHANG Z G, WANG Z Y, YU J Y, Qin H T, MAO X R, JIANG H W, XIN D W, YIN Z G, ZHU R S, LIU C Y, YU W, HU Z B, WU X X, LIU J, CHEN Q S. Meta-analysis and transcriptome profiling reveal hub genes for soybean seed storage composition during seed development. Plant Cell & Environment, 2018, 41(9): 2109-2127.

    [2] Wilson R F. Soybean: market driven research needs//Genetics and Genomics of soybean. NewYork: Springer, 2008: 3-15.

    [3] Bellaloui N, Bruns H A, Abbas H K, Mengistu A, Fisher D K, Reddy K N. Agricultural practices altered soybean seed protein, oil, fatty acids, sugars and minerals in the Midsouth USA. Frontiers in Plant Science, 2015, 6(31): 31-44.

    [4] Spencer M, Pantalone V, Meyer E, Landau-Ellis D, Hyten D. Mapping the Fas locus controlling stearic acid content in soybean. Theoretical and Applied Genetics, 2003, 106(4): 615-619.

    [5] 左進(jìn)華, 董海洲, 侯漢學(xué). 大豆蛋白生產(chǎn)與應(yīng)用現(xiàn)狀. 糧食與油脂, 2007, 5(5): 12-15.

    ZUO J H, DONG H Z, HOU H X. Current status of soy protein production and application. Grains and fats, 2007, 5(5): 12-15. (in Chinese)

    [6] 王連錚. 國內(nèi)外大豆生產(chǎn)的現(xiàn)狀和大豆品種創(chuàng)新問題. 中國食物與營養(yǎng), 2006, 7(6): 6-9.

    WANG L Z. The status quo of soybean production at home and abroad and the innovation of soybean varieties. Chinese food and nutrition, 2006, 7(6): 6-9. (in Chinese)

    [7] 任波, 李毅. 大豆種子脂肪酸合成代謝的研究進(jìn)展. 分子植物育種, 2005(3): 301-306.

    REN B, LI Y. Research progress on fatty acid synthesis and metabolism of soybean seeds.Molecular Plant Breeding, 2005(3): 301-306. (in Chinese)

    [8] Whigham L D, Watras A C, Schoeller D A. Efficacy of conjugated linoleic acid for reducing fat mass: a meta-analysis in humans. American Journal of Clinical Nutrition, 2007(5): 1203-1211.

    [9] Tompkins C, Perkins E G. Frying performance of low-linolenic acid soybean oil. Journal of the American Oil Chemists’ Society, 2000, 77(3): 223-229.

    [10] Slover H T, Lanza E. Quantitative analysis of food fatty acids by capillary gas chromatography. Journal of the American Oil Chemists’ Society, 1979, 56(12): 933-943.

    [11] Stoffel W, Chu F, Ahrens Jr E H. Analysis of long-chain fatty acids by gas-liquid chromatography. Analytical Chemistry, 1959, 31(2): 307-308.

    [12] Pazdernik D L, Killam A S, Orf J H. Analysis of amino and fatty acid composition in soybean seed, using near infrared reflectance spectroscopy. Agronomy Journal, 1997, 89(4): 679-685.

    [13] Sato T, Kawano S, Iwamoto M. Near infrared spectral patterns of fatty acid analysis from fats and oils. Journal of the American Oil Chemists’ Society, 1991, 68(11): 827-833.

    [14] Wood R, Lee T. High-performance liquid chromatography of fatty acids: quantitative analysis of saturated, monoenoic, polyenoic and geometrical isomers. Journal of Chromatography A, 1983, 254(JAN): 237-246.

    [15] Aveldano M I, VanR ollins M, Horrocks L A. Separation and quantitation of free fatty acids and fatty acid methyl esters by reverse phase high pressure liquid chromatography. Journal of Lipid Research, 1983, 24(1): 83-93.

    [16] 范勝栩, 李斌, 孫君明, 韓粉霞, 閆淑榮, 王嵐. 氣相色譜方法定量檢測(cè)大豆5種脂肪酸. 中國油料作物學(xué)報(bào), 2015, 37(4): 548.

    FAN S X, LI B, SUN J M, HAN F X, YAN S R, WANG L. Gas chromatography method for quantitative detection of 5 fatty acids in soybean.Chinese Journal of Oil Crops, 2015, 37 (4): 548. (in Chinese)

    [17] 王芹, 馮景春, 馮開. 氣相色譜法及其應(yīng)用. 廣東化工, 2014, 41(12): 202-208.

    WANG Q, FENG J C, FENG K. Gas chromatography and its application. Guangdong Chemical Industry, 2014, 41(12): 202-208. (in Chinese)

    [18] LI H H, YE G Y, WANG J K. A modified algorithm for the improvement of composite interval mapping. Genetics, 2007, 175(1): 361-374.

    [19] Jansen R C. Interval mapping of multiple quantitative trait loci. Genetics, 1993, 135(1): 205-211.

    [20] KAO C H, ZENG Z B, Teasdale R D. Multiple interval mapping for quantitative trait loci. Genetics, 2004, 152(3): 1987-2002.

    [21] Rodolphe F, Lefort M. A multi-marker model for detecting chromosomal segments displaying QTL activity. Genetics, 1993, 134(4): 1277-1288.

    [22] LI H H, Ribaut J M, LI Z L, WANG J K. Inclusive composite interval mapping (ICIM) for digenic epistasis of quantitative traits in biparental populations. Theoretical & Applied Genetics, 2008, 116(2): 243-260.

    [23] 李慧慧. 數(shù)量性狀基因的完備區(qū)間作圖方法[D]. 北京: 北京師范大學(xué), 2009.

    LI H H. A complete interval mapping method for quantitative trait genes[D]. Beijing: Beijing Normal University, 2009. (in Chinese)

    [24] Weller J I. Maximum likelihood techniques for the mapping and analysis of quantitative trait loci with the aid of genetic markers. Biometrics, 1986, 42(3): 627-640.

    [25] Akond M, LIU S, Boney M, KANTARTZI S K, KASSEM M A. Identification of quantitative trait loci (QTL)underlying protein, oil, and five major fatty acids' contents in soybean. American Journal of Plant Sciences, 2014, 5(1): 158-167.

    [26] QIN H T, LIU Z X, WANG Y Y, XU M Y, QI Z M. Meta-analysis and overview analysis of quantitative trait locis associated with fatty acid content in soybean for candidate gene mining.Plant Breed, 2018, 137(2): 181-193.

    [27] LI B, FAN S X, YU F K, CHEN Y, ZHANG S R, HAN F X, YAN S R, WANG L Z, SUN J M. High-resolution mapping of QTL for fatty acid composition in soybean using specific locus amplified fragment sequencing. Theoretical & Applied Genetics, 2017, 130(7): 1467-1479.

    [28] XIA N, WU D P, ZHAN Y H, LIU Y, SUN M Y, ZHAO X, TENG W L, HAN Y P. Dissection of genetic architecture for oil content in soybean seed using two backcross populations. Plant Breed,2017, 136(7): 365-371.

    [29] FAN S X, LI B, YU F K, HAN F X, YAN S R, WANG L Z, SUN J M. Analysis of additive and epistatic quantitative trait loci underlying fatty acid concentrations in soybean seeds across multiple environments. Euphytica, 201, 206(3): 689-700.

    [30] 盛英華, 張延瑞, 戴亞楠, 昝光敏, 周凱, 王賢智. 不同群體中大豆脂肪酸組分QTL定位研究. 中國油料作物學(xué)報(bào), 2020(5): 796-806.

    SHENG Y H, ZHANG Y R, DAI Y N, ZAN G M, ZHOU K, WANG X Z. QTL mapping of soybean fatty acid components in different populations. Chinese Journal of Oil Crops, 2020(5): 796-806. (in Chinese)

    [31] CHEN Q S, ZHANG Z C, LIU C Y, XIN D W, QIU H M, SHAN D P, SHAN C Y, HU G H. QTL analysis of major agronomic traits in soybean. Scientia Agriculture Sinica, 2007, 6(4): 399-405.

    [32] QI Z M, HUANG L, ZHU R S, XIN D W, LIU C Y, HAN X, JIANG H W, HONG W G, HU G H, ZHENG H K, CHEN Q S. A high-density genetic map for soybean based on specific length amplified fragment sequencing.PLoS ONE, 2014, 9(8): e104871.

    [33] Xin D W, Qi Z M, Jiang H W, Zhang Z G, Zhu R S, Hu J H, HAN H Y, HU G H, LIU C Y, CHEN Q S. QTL locationand epistatic effect analysis of 100-seed weight using wild soybean (Sieb. & Zucc.) chromosome segment substitution lines. Plos ONE, 2016, 11(3): e0149380.

    [34] Mccouch S R, Cho Y G, Yano M, PAUL E, BLINSTRUB M, MORISHIMA H, KINOSITA T. Report on QTL nomenclature. Rice Genetics Newsletter, 1997, 14: 11-13.

    [35] JIANG H W, LI Y Y, QIN H T, LI Y L, QI H D, LI C D, WANG N N, LI R C, ZHAO Y Y, HAUNG S Y, YU J Y, WANG X Y, ZHU R S, LIU C Y, HU Z B, QI Z M, XIN D W, WU X X, CHEN Q S. Identification of major QTLs associated with first pod height and candidate gene mining in soybean. Frontiers in plant science, 2018, 9: 1280

    [36] WANG X Y, LI Q Y, ZHANG Q, YU J Y, QIN H T, QI H D, LI Y L, LI Y Y, YIN Z G, HAN X, WU X X, XIN D W, CHEN Q S, QI Z M. Identification of soybean genes related to fatty acid content based on a soybean genome collinearity analysis.Plant Breeding, 2019, 138(6): 696-707.

    [37] Baud S, Guyon V, Kronenberger J, WUILLEME S, MIQUEL M, CABOCHE M, LEPINIEC L, ROCHAT C. Multifunctional acetyl-CoA carboxylase 1 is essential for very long chain fatty acid elongation and embryo development in. The Plant Journal, 2010, 33(1): 75-86.

    [38] Goettel W, Ramirez M, Upchurch R G, CHARLES Y Q. Identification and characterization of large DNA deletions affecting oil quality traits in soybean seeds through transcriptome sequencing analysis. Theoretical and Applied Genetics, 2016, 129(8): 1577-1593.

    [39] Pettitt T R, Martin A, Horton T, LIOSSIS C, LORD J M, WAKELAM M. Diacylglycerol and phosphatidate generated by phospholipases C and D, respectively, have distinct fatty acid compositions and functions phospholipase d-derived diacylglycerol does not activate protein kinase c in porcine aortic endothelial cells. Journal of Biological Chemistry, 1997, 272(28): 17354-17359.

    [40] Wilson C H, Shalini S, Filipovska A, RICHMAN T R, KUMAR S. Age-related proteostasis and metabolic alterations in Caspase-2-deficient mice. Cell Death & Disease, 2015, 6(1): e1597.

    [41] Bartley I M, Stoker P G, Martin A D E, HATFIELD S G S, KNEE M. Synthesis of aroma compounds by apples supplied with alcohols and methyl esters of fatty acids. Journal of the Science of Food and Agriculture, 1985, 36: 567-574.

    [42] Koteles J. Fatty acid ω-hydroxylases in soybean[D].Canada Ontario: The University of Western Ontario, 2012.

    [43] Ambasht P K, Kayastha A M. Plant pyruvate kinase. Biologia Plantarum, 2002, 45(1): 1-10.

    [44] Andre C, Froehlich J E, Moll M R, BENNING C. A heteromeric plastidic pyruvate kinase complex involved in seed oil biosynthesis in. The Plant Cell, 2007, 19(6): 2006-2022.

    [45] HUANG P Y, LUO L J. Effect on pyruvate kinase in high plants.Journal of Anhui Agricultural Sciences, 2009, 37(20): 9352-9354.

    [46] 蔣洪蔚, 劉春燕, 高運(yùn)來, 李燦東, 張聞博, 胡國華, 陳慶山. 作物QTL定位常用作圖群體. 生物技術(shù)通報(bào), 2008, 1(20): 12-17.

    JIANG H W, LIU C Y, GAO Y L, LI C D, ZHANG W B, HU G H, CHEN Q S. Crop QTL mapping is often used as a map population. Biotechnology Bulletin, 2008, 1(20): 12-17. (in Chinese)

    [47] Ma X, Chen X P, Zhao J, WANG S S, TAN L B, SUN C Q, LIU F X. Identification of QTLs related to cadmium tolerance from wild rice () using a high-density genetic map for a set of introgression lines. Euphytica, 2019, 215(12): 1-12.

    [48] 李晶晶, 王利鋒, 馬娟, 曹言勇, 王浩, 王麗艷, 賈騰蛟, 董春林, 李會(huì)勇. 基于昌7-2導(dǎo)入系發(fā)掘干旱脅迫下玉米產(chǎn)量相關(guān)QTL位點(diǎn). 玉米科學(xué), 2019, 27(4): 64-70.

    LI J J, WANG L F, MA J, CAO Y Y, WANG H, WANG L Y, JIA T J, DONG C L, LI H Y. Discovery of QTLs related to maize yield under drought stress based on Chang 7-2 introduced line. Maize Science, 2019, 27(4): 64-70. (in Chinese)

    [49] 于福寬. 大豆種質(zhì)脂肪酸主要組分鑒定與QTL標(biāo)記定位[D]. 北京: 中國農(nóng)業(yè)科學(xué)院, 2011.

    YU F K. Identification of the main fatty acid components of soybean germplasm and QTL mapping[D].Beijing:Chinese Academy of Agricultural Sciences, 2011. (in Chinese)

    [50] 朱明月. 利用回交導(dǎo)入系群體定位大豆蛋白質(zhì)、脂肪含量及脂肪酸含量QTL[D]. 北京: 中國農(nóng)業(yè)科學(xué)院, 2017.

    ZHU M Y. Using backcross introduction line population to locate soybean protein, fat content and fatty acid content QTL[D]. Beijing:Chinese Academy of Agricultural Sciences, 2017. (in Chinese)

    [51] Diers B W, Keim P, Fehr W R, SHOEMAKER R C. RFLP analysis of soybean seed protein and oil content. Theoretical and Applied Genetic, 1992, 83(5): 608-612.

    [52] Hyten D L, Pantalone V R, Saxton A M, SCHMIDT M E, SAMS C E. Molecular mapping and identification of soybean fatty acid modifier quantitative trait loci. Journal of the American Oil Chemists Society, 2004, 81(12): 1115-1118.

    [53] 葉桑, 崔翠, 郜歡歡, 雷維, 王劉艷, 王瑞莉, 陳柳依, 曲存民, 唐章林, 李加納. 基于SNP遺傳圖譜對(duì)甘藍(lán)型油菜部分脂肪酸組成性狀的QTL定位. 中國農(nóng)業(yè)科學(xué), 2019, 52(21): 26-40.

    YE S, CUI C, GAO H H, LEI W, WANG L Y, WANG R L, CHEN L Y, QU C M, TANG Z L, LI J N. QTL mapping of some fatty acid composition traits inbased on SNP genetic map. China Agricultural Sciences, 2019, 52(21): 26-40. (in Chinese)

    Mapping QTL for Soybean Fatty Acid Composition Based on RIL and CSSL Population

    QU KeXin, HAN Lu, XIE JianGuo, PAN WenJing, ZHANG ZeXin, XIN DaWei, LIU ChunYan, CHEN QingShan, QI ZhaoMing

    Soybean Genetic Improvement Laboratory, College of Agriculture, Northeast Agricultural University, Harbin 150030

    【】Soybeans () originated from China. High-quality soybeans are widely used in various processing industries such as food, feeding, textiles, etc. Therefore, high-quality soybean breeding is a key point for soybean breeders and producers. This study conducted QTL mapping of each component of soybean fatty acid and screening of candidate genes, which would lay the molecular foundation for soybean quality improvement. 【】A recombinant inbred lines (RILs) population crossed by Charleston (American soybean varieties ) and Dongnong 594, and a chromosome segment substitution lines (CSSLs) population crossed by Suinong 14 (cultivated soybean) and ZYD00006 (wild soybean) were used for QTL mapping. We used gas chromatography to determine the fatty acid content of these two populations. As the genetic maps have been published by the soybean genetic improvement laboratory of the Agricultural College of Northeast Agricultural University before, QTL mapping of soybean fatty acid components in RIL and CSL populations were performed by the Windows QTL Cartographer 2.5 and ICIMapping software. And the candidate genes were screened from the QTL interval. 【】Based on 2017 to 2018 years data, 34 and 20 QTLs related to fatty acid components were mapped in the RIL population and the CSSL population, respectively. These QTLs distribute in 13 linkage groups except B2, C1, G, H, J, M, and O. QTL mapping of the two populations was compared that ten pairs of QTLs were detected in the two populations. We found that QTLs distributed in the A1, C2, D1a, F, K, and N linkage groups were related to the content of multiple fatty acids components. An overlapping QTL related to linoleic acid and oil content was detected on the A1 linkage group, QTL related to stearic acid and oil content on the C2, QTL related to stearic acid and oil content on the D1a, QTLs related to palmitic acid, stearic acid and oil content on the F, QTLs related to linoleic acid and linolenic acid content on the K, QTLs related to palmitic acid and oil content, and QTL related to oleic acid and linoleic acid content on the N. Candidate genes were screened out from QTL intervals. In total, 485 candidate genes were screened from the gene annotation data set and 271 of them annotated within GO annotations. GO enrichment analysis showed that 15 candidate genes involved in fatty acids pathway. These genes affect synthesis of fatty acids mainly through encoding plant acyl carrier protein (ACP) thioesterase, fatty acid desaturase, phospholipase D1, fatty acid-hydroxylase and pyruvate kinase, participating in the biosynthesis of acyl-CoA, and regulating the extension of fatty acid chain. 【】54 QTLs related to soybean fatty acid were detected, and 10 pairs QTLs were stable detected from the two mapping populations. We used the confidence intervals from QTL mapping to screen candidate genes, and 15 candidate genes related to fatty acids pathway were screened out. These stable QTLs and candidate genes can be used for molecular marker-assisted selection of soybean fatty acid improvement.

    soybean; introduction line; gas chromatography; QTL mapping; gene mining; enrichment analysis

    10.3864/j.issn.0578-1752.2021.15.003

    2021-02-03;

    2021-03-22

    國家自然科學(xué)基金(31701449)、東北農(nóng)業(yè)大學(xué)青年人才骨干項(xiàng)目(18XG01)

    渠可心,E-mail:18145648226@163.com。通信作者陳慶山,E-mail:qshchen@126.com。通信作者齊照明,E-mail:qizhaoming1860@126.com

    (責(zé)任編輯 李莉)

    猜你喜歡
    油分連鎖組分
    組分分發(fā)管理系統(tǒng)在天然氣計(jì)量的應(yīng)用
    油分芯導(dǎo)致TSA-2.4A型空氣壓縮機(jī)著火原因分析
    一種難溶難熔未知組分板材的定性分析
    大豆關(guān)聯(lián)重組自交系群體蛋白質(zhì)、油分含量的QTL分析
    專注零售連鎖空間打造
    烤煙油分判定標(biāo)準(zhǔn)的修訂與應(yīng)用探討
    黑順片不同組分對(duì)正常小鼠的急性毒性
    中成藥(2018年8期)2018-08-29 01:28:26
    庫里受傷的連鎖效應(yīng)
    NBA特刊(2018年7期)2018-06-08 05:48:32
    布拉格Burrito Loco連鎖快餐店
    金雀花中黃酮苷類組分鑒定及2種成分測(cè)定
    中成藥(2018年2期)2018-05-09 07:20:09
    最近最新中文字幕大全免费视频| svipshipincom国产片| 久久国产亚洲av麻豆专区| 激情在线观看视频在线高清| 欧美一区二区精品小视频在线| 久久婷婷人人爽人人干人人爱| 欧美激情高清一区二区三区| 国产亚洲av嫩草精品影院| 欧洲精品卡2卡3卡4卡5卡区| 亚洲专区字幕在线| 久久精品人妻少妇| av有码第一页| 一进一出抽搐动态| 高潮久久久久久久久久久不卡| 一级a爱视频在线免费观看| 国产野战对白在线观看| 制服丝袜大香蕉在线| 午夜老司机福利片| 99久久99久久久精品蜜桃| 老汉色∧v一级毛片| 亚洲激情在线av| 大型av网站在线播放| 桃色一区二区三区在线观看| 久久精品国产综合久久久| 12—13女人毛片做爰片一| 亚洲,欧美精品.| 亚洲七黄色美女视频| 久久九九热精品免费| 国产色视频综合| 久9热在线精品视频| 亚洲色图 男人天堂 中文字幕| 这个男人来自地球电影免费观看| 真人一进一出gif抽搐免费| 99久久国产精品久久久| www国产在线视频色| 久久伊人香网站| 男女那种视频在线观看| 真人做人爱边吃奶动态| 国产伦在线观看视频一区| 白带黄色成豆腐渣| 欧美三级亚洲精品| 日本免费一区二区三区高清不卡| 亚洲自拍偷在线| 久久天堂一区二区三区四区| 91国产中文字幕| 国产人伦9x9x在线观看| 成人特级黄色片久久久久久久| 亚洲精品av麻豆狂野| 色老头精品视频在线观看| 日韩 欧美 亚洲 中文字幕| 午夜亚洲福利在线播放| 美女扒开内裤让男人捅视频| 精品少妇一区二区三区视频日本电影| 亚洲,欧美精品.| 亚洲五月婷婷丁香| 国产av一区在线观看免费| 曰老女人黄片| 草草在线视频免费看| 亚洲aⅴ乱码一区二区在线播放 | 欧美 亚洲 国产 日韩一| 午夜免费成人在线视频| 午夜福利免费观看在线| 后天国语完整版免费观看| 亚洲人成网站高清观看| 欧美三级亚洲精品| 90打野战视频偷拍视频| 久久久国产成人精品二区| 国产伦在线观看视频一区| 久久国产乱子伦精品免费另类| 黄色视频,在线免费观看| 夜夜夜夜夜久久久久| 久久精品国产综合久久久| 午夜影院日韩av| 国产精品av久久久久免费| 97超级碰碰碰精品色视频在线观看| 亚洲av成人一区二区三| 亚洲欧美日韩无卡精品| 成人午夜高清在线视频 | 亚洲成人国产一区在线观看| 后天国语完整版免费观看| 亚洲 欧美一区二区三区| 国内少妇人妻偷人精品xxx网站 | 久久香蕉国产精品| aaaaa片日本免费| 欧美日本视频| 啦啦啦免费观看视频1| 免费人成视频x8x8入口观看| 一本精品99久久精品77| 国产午夜福利久久久久久| 欧美黄色片欧美黄色片| 国产国语露脸激情在线看| 真人做人爱边吃奶动态| 国产在线观看jvid| 色精品久久人妻99蜜桃| 人妻丰满熟妇av一区二区三区| 老汉色∧v一级毛片| 久久久久久国产a免费观看| bbb黄色大片| 一区二区三区国产精品乱码| 美女大奶头视频| 欧美成人性av电影在线观看| 欧美国产日韩亚洲一区| 精品一区二区三区视频在线观看免费| av有码第一页| 校园春色视频在线观看| 国产一级毛片七仙女欲春2 | 国产精品98久久久久久宅男小说| 美女高潮喷水抽搐中文字幕| 亚洲人成电影免费在线| 一本大道久久a久久精品| 国产在线观看jvid| 久久性视频一级片| 日日爽夜夜爽网站| 美女免费视频网站| 听说在线观看完整版免费高清| 午夜老司机福利片| 欧美成人一区二区免费高清观看 | 亚洲一区中文字幕在线| 欧美午夜高清在线| 麻豆国产av国片精品| 国产1区2区3区精品| www.自偷自拍.com| 69av精品久久久久久| 国产久久久一区二区三区| 亚洲 欧美 日韩 在线 免费| 国产一区二区三区在线臀色熟女| 欧美中文综合在线视频| 999精品在线视频| 美女高潮到喷水免费观看| x7x7x7水蜜桃| 成人国产综合亚洲| 国产精品免费一区二区三区在线| 18禁黄网站禁片午夜丰满| 91在线观看av| 丰满的人妻完整版| 好男人在线观看高清免费视频 | 久久久久久久精品吃奶| 国产真实乱freesex| 91老司机精品| 亚洲精品在线美女| 国产精品,欧美在线| 久久久国产成人免费| 日日爽夜夜爽网站| 欧美av亚洲av综合av国产av| 国产精品 国内视频| 制服人妻中文乱码| 欧美又色又爽又黄视频| 免费看十八禁软件| 中文字幕av电影在线播放| 国产男靠女视频免费网站| 日本a在线网址| 色尼玛亚洲综合影院| 少妇 在线观看| 波多野结衣高清作品| 日韩有码中文字幕| 黑人操中国人逼视频| 久久国产精品男人的天堂亚洲| 亚洲五月天丁香| 正在播放国产对白刺激| 少妇粗大呻吟视频| 成在线人永久免费视频| 级片在线观看| 亚洲精品在线观看二区| 欧美zozozo另类| 淫妇啪啪啪对白视频| 亚洲专区中文字幕在线| 免费人成视频x8x8入口观看| 久久热在线av| 国产精品久久久久久亚洲av鲁大| 国产伦人伦偷精品视频| 精品久久久久久久久久免费视频| 欧美av亚洲av综合av国产av| 成人手机av| 看黄色毛片网站| 国产精品香港三级国产av潘金莲| av福利片在线| 美女国产高潮福利片在线看| 国产av不卡久久| 国产97色在线日韩免费| 身体一侧抽搐| 男女视频在线观看网站免费 | 欧美成人午夜精品| 国产精品爽爽va在线观看网站 | 亚洲精品在线观看二区| 桃红色精品国产亚洲av| 日本熟妇午夜| 国产私拍福利视频在线观看| 女警被强在线播放| 国产亚洲av高清不卡| 女性生殖器流出的白浆| 国产精品永久免费网站| 99国产极品粉嫩在线观看| av视频在线观看入口| xxx96com| 成人三级黄色视频| 真人一进一出gif抽搐免费| 搞女人的毛片| 亚洲精品色激情综合| 老司机午夜十八禁免费视频| 在线观看66精品国产| 国产欧美日韩精品亚洲av| e午夜精品久久久久久久| 精品少妇一区二区三区视频日本电影| 亚洲性夜色夜夜综合| 欧洲精品卡2卡3卡4卡5卡区| 免费在线观看影片大全网站| 国产精品久久久久久人妻精品电影| 久久久久久久久免费视频了| 国产真人三级小视频在线观看| 久9热在线精品视频| 757午夜福利合集在线观看| 十八禁网站免费在线| 欧美黄色淫秽网站| 欧美一区二区精品小视频在线| 亚洲av熟女| 久久热在线av| 人人妻人人澡欧美一区二区| 国产亚洲av嫩草精品影院| 又黄又粗又硬又大视频| 久久精品影院6| 狂野欧美激情性xxxx| 色哟哟哟哟哟哟| 亚洲七黄色美女视频| 亚洲在线自拍视频| 国产主播在线观看一区二区| 午夜老司机福利片| 999久久久精品免费观看国产| 99国产综合亚洲精品| 亚洲天堂国产精品一区在线| 亚洲国产精品999在线| 成年人黄色毛片网站| 无限看片的www在线观看| 天堂动漫精品| 国产精品永久免费网站| 亚洲第一av免费看| 日本撒尿小便嘘嘘汇集6| 侵犯人妻中文字幕一二三四区| 欧美日韩中文字幕国产精品一区二区三区| 最近最新中文字幕大全电影3 | 黄色片一级片一级黄色片| 老熟妇乱子伦视频在线观看| 国产精品av久久久久免费| 久久久久久久久免费视频了| 久久精品国产亚洲av香蕉五月| 很黄的视频免费| 一区二区三区激情视频| 色精品久久人妻99蜜桃| 欧美日韩瑟瑟在线播放| 精品久久久久久久毛片微露脸| 午夜福利成人在线免费观看| 午夜免费成人在线视频| 亚洲成人久久爱视频| 国语自产精品视频在线第100页| av有码第一页| 人妻丰满熟妇av一区二区三区| 国产精品香港三级国产av潘金莲| 欧美黄色淫秽网站| 精品电影一区二区在线| 一区二区三区精品91| 日本熟妇午夜| 香蕉丝袜av| 啦啦啦观看免费观看视频高清| 中文字幕人成人乱码亚洲影| 欧美中文日本在线观看视频| 伊人久久大香线蕉亚洲五| 1024视频免费在线观看| 色老头精品视频在线观看| a级毛片a级免费在线| 一本精品99久久精品77| 久久精品国产亚洲av香蕉五月| x7x7x7水蜜桃| 不卡一级毛片| or卡值多少钱| 人妻久久中文字幕网| 久久亚洲真实| 97超级碰碰碰精品色视频在线观看| 又黄又爽又免费观看的视频| 黄色女人牲交| 在线观看午夜福利视频| 亚洲精品国产一区二区精华液| 99久久久亚洲精品蜜臀av| 琪琪午夜伦伦电影理论片6080| 老司机福利观看| 亚洲成人久久性| 一夜夜www| 最近最新中文字幕大全电影3 | 天堂√8在线中文| 狂野欧美激情性xxxx| 女同久久另类99精品国产91| 国产黄a三级三级三级人| 欧美成人免费av一区二区三区| 每晚都被弄得嗷嗷叫到高潮| 久久久久久久午夜电影| 欧美成人免费av一区二区三区| 少妇 在线观看| 久久香蕉国产精品| 欧美成人免费av一区二区三区| 麻豆成人av在线观看| 国产免费av片在线观看野外av| 午夜免费激情av| 日韩欧美在线二视频| 国产精品国产高清国产av| av欧美777| 99riav亚洲国产免费| 黄片播放在线免费| www.自偷自拍.com| 黑丝袜美女国产一区| 成年女人毛片免费观看观看9| 亚洲av熟女| 欧美成人免费av一区二区三区| 国产黄片美女视频| av欧美777| 视频在线观看一区二区三区| 高潮久久久久久久久久久不卡| 少妇被粗大的猛进出69影院| 成人18禁在线播放| 日韩欧美一区二区三区在线观看| 91九色精品人成在线观看| 91麻豆精品激情在线观看国产| 熟女少妇亚洲综合色aaa.| 欧美av亚洲av综合av国产av| 久久久国产成人精品二区| 婷婷精品国产亚洲av在线| 久久香蕉激情| 熟女电影av网| 麻豆av在线久日| 成人一区二区视频在线观看| 嫩草影视91久久| 天天添夜夜摸| 成人亚洲精品av一区二区| 日本成人三级电影网站| 好男人电影高清在线观看| 免费av毛片视频| 精品熟女少妇八av免费久了| 日本熟妇午夜| 欧美激情 高清一区二区三区| 亚洲成av人片免费观看| 国产成人精品久久二区二区免费| 久久99热这里只有精品18| 久久久久国产一级毛片高清牌| 国产精品影院久久| 夜夜夜夜夜久久久久| 日韩欧美国产一区二区入口| 97超级碰碰碰精品色视频在线观看| 韩国精品一区二区三区| 在线看三级毛片| 此物有八面人人有两片| xxxwww97欧美| 亚洲 国产 在线| 欧美人与性动交α欧美精品济南到| 亚洲精品中文字幕在线视频| 香蕉国产在线看| 亚洲电影在线观看av| 两个人免费观看高清视频| 俺也久久电影网| 曰老女人黄片| 1024视频免费在线观看| 成人特级黄色片久久久久久久| 亚洲精华国产精华精| 美女国产高潮福利片在线看| 色精品久久人妻99蜜桃| 夜夜看夜夜爽夜夜摸| 中国美女看黄片| 日本成人三级电影网站| 在线观看午夜福利视频| 久久久久久免费高清国产稀缺| 青草久久国产| 999精品在线视频| 后天国语完整版免费观看| 不卡一级毛片| 一进一出抽搐动态| www.自偷自拍.com| 久久久久久亚洲精品国产蜜桃av| 啦啦啦免费观看视频1| 老司机福利观看| 国产不卡一卡二| 色综合欧美亚洲国产小说| 亚洲国产看品久久| 午夜成年电影在线免费观看| 久久精品aⅴ一区二区三区四区| 成人亚洲精品av一区二区| 日韩 欧美 亚洲 中文字幕| 一卡2卡三卡四卡精品乱码亚洲| 波多野结衣av一区二区av| 国产亚洲欧美98| 男女午夜视频在线观看| 欧美成狂野欧美在线观看| 亚洲一区中文字幕在线| 黄色毛片三级朝国网站| 婷婷精品国产亚洲av在线| 男女视频在线观看网站免费 | 久久精品成人免费网站| 淫妇啪啪啪对白视频| 欧美在线黄色| 久久午夜亚洲精品久久| 老司机深夜福利视频在线观看| 天堂动漫精品| 色av中文字幕| 999久久久国产精品视频| 久久久久久久精品吃奶| 精品高清国产在线一区| 草草在线视频免费看| 老司机午夜十八禁免费视频| 亚洲电影在线观看av| 亚洲七黄色美女视频| 国产99久久九九免费精品| 啦啦啦免费观看视频1| 久99久视频精品免费| 日本成人三级电影网站| 91av网站免费观看| 曰老女人黄片| 国产精华一区二区三区| 18禁观看日本| av片东京热男人的天堂| 欧美在线黄色| 国产aⅴ精品一区二区三区波| 波多野结衣高清作品| 日日摸夜夜添夜夜添小说| av有码第一页| 超碰成人久久| 国产精品 国内视频| 亚洲中文字幕日韩| 午夜精品久久久久久毛片777| 韩国精品一区二区三区| 精品欧美一区二区三区在线| 一区二区日韩欧美中文字幕| 国产又黄又爽又无遮挡在线| 啦啦啦免费观看视频1| 热99re8久久精品国产| 老司机在亚洲福利影院| 变态另类丝袜制服| 亚洲avbb在线观看| ponron亚洲| 国产成人影院久久av| 国产黄片美女视频| 午夜精品在线福利| 伊人久久大香线蕉亚洲五| 成年免费大片在线观看| 久热这里只有精品99| 国产在线精品亚洲第一网站| 久久香蕉国产精品| 免费在线观看视频国产中文字幕亚洲| 久久精品91无色码中文字幕| 12—13女人毛片做爰片一| 日韩一卡2卡3卡4卡2021年| 久久人人精品亚洲av| 99热这里只有精品一区 | 欧美精品亚洲一区二区| 欧美色视频一区免费| 大型黄色视频在线免费观看| 精品欧美一区二区三区在线| 免费电影在线观看免费观看| 日日夜夜操网爽| 精品久久久久久久毛片微露脸| 不卡av一区二区三区| 啪啪无遮挡十八禁网站| 国产亚洲欧美精品永久| 国产精品久久久人人做人人爽| 亚洲在线自拍视频| 亚洲狠狠婷婷综合久久图片| 久久久久国产精品人妻aⅴ院| 久久久久国内视频| 国产高清视频在线播放一区| 久久精品国产99精品国产亚洲性色| 婷婷六月久久综合丁香| 99精品在免费线老司机午夜| 叶爱在线成人免费视频播放| 国产精品免费视频内射| 可以在线观看的亚洲视频| 最好的美女福利视频网| 高清在线国产一区| 久久久久久久午夜电影| 国产精品国产高清国产av| 国产精华一区二区三区| 午夜免费成人在线视频| 亚洲国产毛片av蜜桃av| 国产精品永久免费网站| 国产三级黄色录像| 免费在线观看视频国产中文字幕亚洲| 久久久久九九精品影院| 国产麻豆成人av免费视频| av天堂在线播放| 色综合欧美亚洲国产小说| 在线免费观看的www视频| 黄色毛片三级朝国网站| 国产精品精品国产色婷婷| 激情在线观看视频在线高清| 最近在线观看免费完整版| 午夜影院日韩av| 丰满人妻熟妇乱又伦精品不卡| 国产精品一区二区三区四区久久 | 久久精品aⅴ一区二区三区四区| 欧美一级a爱片免费观看看 | 婷婷丁香在线五月| 亚洲久久久国产精品| 精品国内亚洲2022精品成人| 搡老妇女老女人老熟妇| 午夜两性在线视频| 精品一区二区三区四区五区乱码| 日韩 欧美 亚洲 中文字幕| 精品久久久久久久久久免费视频| 久久人妻福利社区极品人妻图片| 在线观看66精品国产| bbb黄色大片| 日韩精品免费视频一区二区三区| 国产成人av激情在线播放| 国产精品av久久久久免费| 在线观看66精品国产| 女警被强在线播放| 国产高清激情床上av| 亚洲精品av麻豆狂野| 两个人视频免费观看高清| 不卡一级毛片| 精品欧美一区二区三区在线| 国产精品一区二区三区四区久久 | 亚洲无线在线观看| 国产精品98久久久久久宅男小说| 国产国语露脸激情在线看| 国内揄拍国产精品人妻在线 | 亚洲第一青青草原| netflix在线观看网站| 日韩av在线大香蕉| 999久久久国产精品视频| 亚洲精品国产一区二区精华液| 亚洲精品美女久久av网站| 欧美精品啪啪一区二区三区| 久9热在线精品视频| 欧美在线黄色| 成人欧美大片| 成人午夜高清在线视频 | 99久久综合精品五月天人人| 亚洲精品美女久久久久99蜜臀| 两性午夜刺激爽爽歪歪视频在线观看 | 午夜福利一区二区在线看| 久久国产亚洲av麻豆专区| 免费搜索国产男女视频| 午夜视频精品福利| 人妻久久中文字幕网| 黄色视频不卡| 激情在线观看视频在线高清| 一进一出抽搐动态| 少妇被粗大的猛进出69影院| 欧美日本视频| 国产亚洲av嫩草精品影院| 在线视频色国产色| 黄片播放在线免费| 最近最新中文字幕大全免费视频| 一进一出抽搐gif免费好疼| 亚洲男人的天堂狠狠| 久久久久久久午夜电影| 变态另类成人亚洲欧美熟女| 在线看三级毛片| 亚洲第一青青草原| 两性午夜刺激爽爽歪歪视频在线观看 | 黄色片一级片一级黄色片| 婷婷六月久久综合丁香| 9191精品国产免费久久| 久久精品成人免费网站| 亚洲美女黄片视频| 日韩大尺度精品在线看网址| 人成视频在线观看免费观看| 97碰自拍视频| 制服人妻中文乱码| 性色av乱码一区二区三区2| 亚洲国产欧美日韩在线播放| 特大巨黑吊av在线直播 | 欧美 亚洲 国产 日韩一| 1024手机看黄色片| 午夜a级毛片| 精品卡一卡二卡四卡免费| 俺也久久电影网| 亚洲国产毛片av蜜桃av| 日日夜夜操网爽| 中亚洲国语对白在线视频| 亚洲真实伦在线观看| 精品久久久久久久久久免费视频| 亚洲av熟女| 自线自在国产av| 国产片内射在线| 免费观看人在逋| 日韩欧美 国产精品| 欧美国产日韩亚洲一区| 丝袜美腿诱惑在线| 国产午夜福利久久久久久| 午夜福利成人在线免费观看| 亚洲无线在线观看| 亚洲电影在线观看av| 成熟少妇高潮喷水视频| 国产一区二区激情短视频| 精品不卡国产一区二区三区| 国产成人欧美在线观看| 女生性感内裤真人,穿戴方法视频| 在线观看66精品国产| 久久久国产成人精品二区| 黄色视频不卡| 国产视频一区二区在线看| 久久久久国产一级毛片高清牌| 女性生殖器流出的白浆| 黄网站色视频无遮挡免费观看| 色综合欧美亚洲国产小说| 90打野战视频偷拍视频| 久久久久久免费高清国产稀缺| 99热这里只有精品一区 | 亚洲欧美日韩无卡精品| 亚洲精品av麻豆狂野| 又大又爽又粗| 国产真人三级小视频在线观看| 看黄色毛片网站| 成年人黄色毛片网站| 一区二区日韩欧美中文字幕| 成人永久免费在线观看视频| 欧美色视频一区免费| 丰满的人妻完整版| 欧美性长视频在线观看| 午夜老司机福利片|