• <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
    精品卡一卡二卡四卡免费| 欧美日韩av久久| 日韩 亚洲 欧美在线| 亚洲国产精品国产精品| 免费黄频网站在线观看国产| 亚洲成av片中文字幕在线观看| 亚洲人成77777在线视频| 久久久亚洲精品成人影院| 99香蕉大伊视频| 亚洲av福利一区| 亚洲精品久久午夜乱码| 久久久久国产精品人妻一区二区| 在线观看一区二区三区激情| 亚洲成人av在线免费| 国产精品 国内视频| 日本色播在线视频| 亚洲av电影在线观看一区二区三区| 美国免费a级毛片| 啦啦啦啦在线视频资源| 考比视频在线观看| 永久免费av网站大全| 99热国产这里只有精品6| √禁漫天堂资源中文www| 99九九在线精品视频| 性色av一级| 男人添女人高潮全过程视频| 免费观看av网站的网址| 国产一区二区在线观看av| 9色porny在线观看| 精品少妇一区二区三区视频日本电影 | av一本久久久久| 久久精品国产综合久久久| 久久精品久久久久久久性| 精品一区二区三卡| 中文字幕另类日韩欧美亚洲嫩草| 色婷婷av一区二区三区视频| 男男h啪啪无遮挡| 亚洲欧美精品综合一区二区三区| 国产日韩欧美在线精品| 久久人妻熟女aⅴ| 免费在线观看黄色视频的| 欧美中文综合在线视频| av不卡在线播放| 黄色 视频免费看| 美女午夜性视频免费| 久久久久久久久久久免费av| 最黄视频免费看| 美女扒开内裤让男人捅视频| 2021少妇久久久久久久久久久| 悠悠久久av| 男女无遮挡免费网站观看| 亚洲精品美女久久av网站| a级毛片在线看网站| 日韩制服骚丝袜av| 美女福利国产在线| 九九爱精品视频在线观看| 人妻一区二区av| 美女大奶头黄色视频| 成人国语在线视频| 美国免费a级毛片| 久久久久久久国产电影| 十八禁网站网址无遮挡| 69精品国产乱码久久久| av在线播放精品| 麻豆乱淫一区二区| 黄片小视频在线播放| 亚洲精品aⅴ在线观看| 99热网站在线观看| 欧美乱码精品一区二区三区| 精品人妻一区二区三区麻豆| 色网站视频免费| 99香蕉大伊视频| 精品一区二区三区av网在线观看 | 亚洲国产日韩一区二区| 最黄视频免费看| 黄片播放在线免费| 欧美亚洲日本最大视频资源| 欧美日本中文国产一区发布| 中文字幕人妻熟女乱码| 一二三四中文在线观看免费高清| 999久久久国产精品视频| 一本一本久久a久久精品综合妖精| 久久久久久久久久久久大奶| av一本久久久久| 黄色视频在线播放观看不卡| 考比视频在线观看| 久久99精品国语久久久| 亚洲情色 制服丝袜| 国产精品免费视频内射| 免费看不卡的av| 中文天堂在线官网| 多毛熟女@视频| 成年人午夜在线观看视频| 少妇精品久久久久久久| 亚洲第一区二区三区不卡| 免费不卡黄色视频| 我要看黄色一级片免费的| 亚洲av成人不卡在线观看播放网 | 亚洲精品久久成人aⅴ小说| av.在线天堂| 免费不卡黄色视频| 女性被躁到高潮视频| 丁香六月天网| 肉色欧美久久久久久久蜜桃| 又粗又硬又长又爽又黄的视频| 最近手机中文字幕大全| 日本欧美视频一区| 美女大奶头黄色视频| 2018国产大陆天天弄谢| 久热爱精品视频在线9| 男女无遮挡免费网站观看| 亚洲av日韩在线播放| 一级毛片 在线播放| 午夜日韩欧美国产| 天天影视国产精品| 日韩av不卡免费在线播放| 夫妻午夜视频| 亚洲av综合色区一区| 熟女av电影| 汤姆久久久久久久影院中文字幕| 黄片播放在线免费| 自线自在国产av| 亚洲av日韩精品久久久久久密 | 久久久久精品久久久久真实原创| 丰满少妇做爰视频| 午夜影院在线不卡| 国产成人精品福利久久| 日本色播在线视频| 欧美日韩国产mv在线观看视频| 九九爱精品视频在线观看| 女性生殖器流出的白浆| 久久久久视频综合| 777久久人妻少妇嫩草av网站| 久久精品熟女亚洲av麻豆精品| 最近最新中文字幕大全免费视频 | 国产精品久久久久久精品电影小说| 美女午夜性视频免费| 亚洲精品aⅴ在线观看| 国产精品三级大全| 国产精品久久久久久精品古装| 一区二区日韩欧美中文字幕| 国产精品偷伦视频观看了| 青春草亚洲视频在线观看| 国产日韩欧美视频二区| 看免费成人av毛片| 国产成人91sexporn| 咕卡用的链子| 午夜福利一区二区在线看| 一级,二级,三级黄色视频| 久久影院123| 日韩一区二区三区影片| 看免费成人av毛片| 中文字幕色久视频| 青春草国产在线视频| 亚洲男人天堂网一区| 69精品国产乱码久久久| 一本久久精品| 午夜福利一区二区在线看| 日韩大码丰满熟妇| 一二三四中文在线观看免费高清| 国产成人系列免费观看| 亚洲av国产av综合av卡| 久久精品久久精品一区二区三区| av在线老鸭窝| 丝袜人妻中文字幕| 国产1区2区3区精品| 伊人久久大香线蕉亚洲五| 天天添夜夜摸| 熟女少妇亚洲综合色aaa.| 9热在线视频观看99| 国产一卡二卡三卡精品 | a 毛片基地| 欧美变态另类bdsm刘玥| 少妇被粗大猛烈的视频| 又黄又粗又硬又大视频| 久久天堂一区二区三区四区| 一区二区av电影网| 日韩 亚洲 欧美在线| 午夜福利,免费看| 一区在线观看完整版| 午夜日韩欧美国产| 欧美变态另类bdsm刘玥| 午夜老司机福利片| 国产精品嫩草影院av在线观看| 国产免费又黄又爽又色| 香蕉国产在线看| 丁香六月欧美| 亚洲第一av免费看| 黄片无遮挡物在线观看| 中文字幕精品免费在线观看视频| 热99久久久久精品小说推荐| 日韩制服骚丝袜av| 亚洲五月色婷婷综合| 免费在线观看完整版高清| 哪个播放器可以免费观看大片| 青春草亚洲视频在线观看| 亚洲四区av| 欧美在线黄色| 成人黄色视频免费在线看| 卡戴珊不雅视频在线播放| 最近的中文字幕免费完整| 日韩精品免费视频一区二区三区| 男女下面插进去视频免费观看| 欧美日韩福利视频一区二区| 国产精品久久久久久精品电影小说| 熟女av电影| 欧美精品一区二区免费开放| 中文乱码字字幕精品一区二区三区| 热99久久久久精品小说推荐| 性少妇av在线| 一区二区三区激情视频| 欧美日韩一区二区视频在线观看视频在线| 黄色怎么调成土黄色| 久久久国产欧美日韩av| 日本午夜av视频| 最近手机中文字幕大全| 精品少妇黑人巨大在线播放| 成人三级做爰电影| 纯流量卡能插随身wifi吗| 亚洲av日韩在线播放| 久久久久久久久免费视频了| 国产免费又黄又爽又色| 视频区图区小说| www日本在线高清视频| 尾随美女入室| 成人国产av品久久久| 国产av精品麻豆| 夫妻午夜视频| 波多野结衣av一区二区av| 视频在线观看一区二区三区| 男男h啪啪无遮挡| 成人午夜精彩视频在线观看| 亚洲成国产人片在线观看| 丰满饥渴人妻一区二区三| 夫妻性生交免费视频一级片| 国产精品三级大全| 无遮挡黄片免费观看| 男女边摸边吃奶| 男男h啪啪无遮挡| 精品人妻在线不人妻| 99香蕉大伊视频| 男人爽女人下面视频在线观看| 亚洲国产日韩一区二区| 免费久久久久久久精品成人欧美视频| 999久久久国产精品视频| 美女大奶头黄色视频| 美女福利国产在线| 国产精品二区激情视频| 日本猛色少妇xxxxx猛交久久| 纯流量卡能插随身wifi吗| 亚洲av欧美aⅴ国产| 亚洲精品av麻豆狂野| 日韩一区二区视频免费看| 久久这里只有精品19| 一级黄片播放器| 国产av码专区亚洲av| 久久综合国产亚洲精品| 国产av一区二区精品久久| 日韩伦理黄色片| 男人爽女人下面视频在线观看| 黄色一级大片看看| 黄片播放在线免费| 国产午夜精品一二区理论片| 成年女人毛片免费观看观看9 | 女人被躁到高潮嗷嗷叫费观| 波野结衣二区三区在线| 黄色怎么调成土黄色| 青春草视频在线免费观看| 亚洲熟女精品中文字幕| 97人妻天天添夜夜摸| 欧美激情极品国产一区二区三区| 国产精品欧美亚洲77777| 成年av动漫网址| 亚洲精品久久久久久婷婷小说| 丁香六月欧美| 亚洲在久久综合| xxx大片免费视频| 欧美成人精品欧美一级黄| 伦理电影大哥的女人| 午夜久久久在线观看| 色播在线永久视频| 精品人妻在线不人妻| 中文欧美无线码| 极品人妻少妇av视频| 亚洲精品日本国产第一区| 国产欧美亚洲国产| 国产视频首页在线观看| 999精品在线视频| 国产女主播在线喷水免费视频网站| 高清av免费在线| 九九爱精品视频在线观看| 一级毛片黄色毛片免费观看视频| www.精华液| 新久久久久国产一级毛片| 中文字幕人妻丝袜制服| 热re99久久国产66热| 日日啪夜夜爽| 国产精品av久久久久免费| 青青草视频在线视频观看| 老司机靠b影院| 一级片免费观看大全| av天堂久久9| 少妇精品久久久久久久| 视频在线观看一区二区三区| 综合色丁香网| 国产成人一区二区在线| bbb黄色大片| 下体分泌物呈黄色| 19禁男女啪啪无遮挡网站| 中文字幕最新亚洲高清| 熟妇人妻不卡中文字幕| 男人操女人黄网站| 老司机深夜福利视频在线观看 | 精品久久久精品久久久| 亚洲第一青青草原| 99热国产这里只有精品6| 亚洲精品第二区| 男女床上黄色一级片免费看| 久久精品久久久久久噜噜老黄| 激情五月婷婷亚洲| 99久国产av精品国产电影| 亚洲国产毛片av蜜桃av| 99精国产麻豆久久婷婷| 一区在线观看完整版| 精品卡一卡二卡四卡免费| 日韩一区二区三区影片| 久久狼人影院| 亚洲av日韩精品久久久久久密 | 亚洲人成网站在线观看播放| 99热国产这里只有精品6| 免费不卡黄色视频| 又黄又粗又硬又大视频| 免费少妇av软件| 精品少妇黑人巨大在线播放| 免费黄频网站在线观看国产| 搡老岳熟女国产| 在线亚洲精品国产二区图片欧美| √禁漫天堂资源中文www| 久久97久久精品| 韩国av在线不卡| 亚洲第一青青草原| 中文天堂在线官网| av又黄又爽大尺度在线免费看| 亚洲国产最新在线播放| 国产深夜福利视频在线观看| 亚洲,欧美,日韩| 精品国产露脸久久av麻豆| 精品一区二区免费观看| 国产一区二区 视频在线| a级毛片黄视频| 精品少妇一区二区三区视频日本电影 | 人妻人人澡人人爽人人| 巨乳人妻的诱惑在线观看| 国产精品秋霞免费鲁丝片| 亚洲四区av| 精品少妇黑人巨大在线播放| 国产一区亚洲一区在线观看| 日本91视频免费播放| 日韩av在线免费看完整版不卡| 黄色怎么调成土黄色| 涩涩av久久男人的天堂| 午夜老司机福利片| 国产精品免费视频内射| 美女高潮到喷水免费观看| 午夜免费鲁丝| 在线观看人妻少妇| 久久国产精品大桥未久av| 成人国产av品久久久| 国产精品嫩草影院av在线观看| 黄片小视频在线播放| 操出白浆在线播放| 女人被躁到高潮嗷嗷叫费观| 老鸭窝网址在线观看| 国产精品欧美亚洲77777| 男人添女人高潮全过程视频| 亚洲天堂av无毛| 国语对白做爰xxxⅹ性视频网站| 国产成人精品在线电影| 丰满乱子伦码专区| 老熟女久久久| 欧美精品一区二区大全| 亚洲欧美清纯卡通| 少妇精品久久久久久久| 波多野结衣一区麻豆| 国产在视频线精品| 国产高清不卡午夜福利| 久久久久久久久免费视频了| 巨乳人妻的诱惑在线观看| 啦啦啦在线观看免费高清www| 国产亚洲一区二区精品| 中文字幕另类日韩欧美亚洲嫩草| 99九九在线精品视频| 国产精品久久久人人做人人爽| 国产欧美日韩综合在线一区二区| 日韩免费高清中文字幕av| 国产精品av久久久久免费| 男女国产视频网站| 波多野结衣av一区二区av| av又黄又爽大尺度在线免费看| 天堂俺去俺来也www色官网| 美女脱内裤让男人舔精品视频| 天天躁日日躁夜夜躁夜夜| 亚洲第一青青草原| 99re6热这里在线精品视频| 高清黄色对白视频在线免费看| 成人国语在线视频| 国产深夜福利视频在线观看| av在线观看视频网站免费| 一区二区三区乱码不卡18| 国产一区二区激情短视频 | av在线app专区| 国产精品秋霞免费鲁丝片| 97人妻天天添夜夜摸| 久久国产精品大桥未久av| 一本久久精品| 精品少妇黑人巨大在线播放| 国产视频首页在线观看| 精品少妇黑人巨大在线播放| 在线观看三级黄色| 如何舔出高潮| 亚洲精品视频女| 欧美 亚洲 国产 日韩一| 大话2 男鬼变身卡| 国产免费视频播放在线视频| 久久久久久久久久久免费av| 色视频在线一区二区三区| 国产乱来视频区| 欧美日韩视频精品一区| 亚洲一区中文字幕在线| 大陆偷拍与自拍| 国产精品嫩草影院av在线观看| 岛国毛片在线播放| 亚洲国产精品一区三区| 19禁男女啪啪无遮挡网站| 美女扒开内裤让男人捅视频| 国产一区有黄有色的免费视频| 亚洲精品国产一区二区精华液| 国产欧美日韩一区二区三区在线| 这个男人来自地球电影免费观看 | 亚洲国产精品999| 又大又黄又爽视频免费| videos熟女内射| 哪个播放器可以免费观看大片| a 毛片基地| 国产在线免费精品| 精品一区二区三区av网在线观看 | 亚洲伊人久久精品综合| 秋霞在线观看毛片| 成人毛片60女人毛片免费| 国产精品免费视频内射| 国产成人精品无人区| 99久久综合免费| 免费观看人在逋| 美女扒开内裤让男人捅视频| 亚洲精品一区蜜桃| 视频区图区小说| 亚洲国产毛片av蜜桃av| 日本色播在线视频| 欧美黄色片欧美黄色片| 热99久久久久精品小说推荐| 捣出白浆h1v1| 日韩中文字幕视频在线看片| 熟妇人妻不卡中文字幕| 亚洲国产欧美一区二区综合| 女人爽到高潮嗷嗷叫在线视频| 国产在线一区二区三区精| 99香蕉大伊视频| 国产精品熟女久久久久浪| 国产xxxxx性猛交| 肉色欧美久久久久久久蜜桃| 18禁裸乳无遮挡动漫免费视频| 国产 精品1| 在线精品无人区一区二区三| 哪个播放器可以免费观看大片| 国产亚洲一区二区精品| 国产熟女欧美一区二区| 精品一区在线观看国产| a 毛片基地| 亚洲国产日韩一区二区| 国产一区亚洲一区在线观看| 国产精品嫩草影院av在线观看| 香蕉丝袜av| 青春草亚洲视频在线观看| 飞空精品影院首页| 欧美激情高清一区二区三区 | 大话2 男鬼变身卡| 久久97久久精品| av在线老鸭窝| 国产成人欧美在线观看 | 精品第一国产精品| 亚洲精品在线美女| 国产精品国产三级专区第一集| 国产免费视频播放在线视频| 最近中文字幕2019免费版| 人人妻人人添人人爽欧美一区卜| 亚洲国产精品999| 可以免费在线观看a视频的电影网站 | 色婷婷久久久亚洲欧美| 久久久精品免费免费高清| 熟女少妇亚洲综合色aaa.| 嫩草影视91久久| 哪个播放器可以免费观看大片| av国产久精品久网站免费入址| 免费看不卡的av| 一级毛片电影观看| 婷婷色av中文字幕| 狠狠婷婷综合久久久久久88av| 精品少妇内射三级| 乱人伦中国视频| 国产av码专区亚洲av| 国产精品人妻久久久影院| 精品亚洲成a人片在线观看| 亚洲av欧美aⅴ国产| 一级毛片电影观看| 精品国产超薄肉色丝袜足j| 亚洲国产欧美一区二区综合| 久久热在线av| 欧美精品亚洲一区二区| 熟女av电影| 国产精品欧美亚洲77777| 国产欧美日韩一区二区三区在线| 久久久久久久久久久久大奶| 美女国产高潮福利片在线看| 免费看不卡的av| bbb黄色大片| 成年av动漫网址| 日韩一本色道免费dvd| 黄色视频在线播放观看不卡| 国产xxxxx性猛交| 国产精品亚洲av一区麻豆 | 免费少妇av软件| 亚洲欧美成人综合另类久久久| 国产 精品1| 欧美黑人欧美精品刺激| 欧美 亚洲 国产 日韩一| 精品国产乱码久久久久久男人| 免费黄色在线免费观看| 18在线观看网站| av在线app专区| 女人精品久久久久毛片| 免费在线观看完整版高清| 欧美日韩亚洲国产一区二区在线观看 | 黄色视频不卡| 老司机影院成人| 人体艺术视频欧美日本| 久久久精品区二区三区| 九草在线视频观看| 一边摸一边抽搐一进一出视频| 天天添夜夜摸| 精品卡一卡二卡四卡免费| 丝袜人妻中文字幕| 亚洲图色成人| 国产极品粉嫩免费观看在线| 在线观看免费视频网站a站| 免费av中文字幕在线| 女人久久www免费人成看片| 久久99热这里只频精品6学生| 肉色欧美久久久久久久蜜桃| 免费日韩欧美在线观看| 免费观看a级毛片全部| 久久久久久久大尺度免费视频| 又粗又硬又长又爽又黄的视频| 丝袜在线中文字幕| 777米奇影视久久| 九九爱精品视频在线观看| 伦理电影免费视频| 国产毛片在线视频| 中文字幕精品免费在线观看视频| 亚洲,欧美,日韩| 国产老妇伦熟女老妇高清| 日本一区二区免费在线视频| 国产爽快片一区二区三区| 久久久精品94久久精品| av视频免费观看在线观看| 精品久久久精品久久久| 久久精品国产亚洲av高清一级| 国产成人欧美| 国产熟女午夜一区二区三区| 欧美黑人精品巨大| 午夜激情久久久久久久| 91成人精品电影| 看免费成人av毛片| 欧美日韩亚洲国产一区二区在线观看 | 免费观看性生交大片5| 看免费成人av毛片| 久热爱精品视频在线9| 免费女性裸体啪啪无遮挡网站| 欧美 日韩 精品 国产| 午夜日本视频在线| 国产乱人偷精品视频| 各种免费的搞黄视频| 国产亚洲av高清不卡| 狂野欧美激情性xxxx| 久久久国产欧美日韩av| av在线播放精品| 黄频高清免费视频| 国产精品熟女久久久久浪| 国产国语露脸激情在线看| 国产精品国产av在线观看| 国产人伦9x9x在线观看| 亚洲欧美精品自产自拍| 亚洲美女搞黄在线观看| 国产成人欧美在线观看 | av免费观看日本| 我要看黄色一级片免费的| 80岁老熟妇乱子伦牲交| 午夜福利影视在线免费观看| 国产 一区精品| 乱人伦中国视频| 午夜免费鲁丝| 99国产精品免费福利视频| 黑人欧美特级aaaaaa片| 国产xxxxx性猛交| 一区二区av电影网| 97人妻天天添夜夜摸|