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

    玉米抗灰斑病QTL元分析及其驗證

    2016-07-15 05:59:37宋茂興張曠野孫銘澤李鳳海鐘雪梅杜萬里呂香玲
    作物學(xué)報 2016年5期
    關(guān)鍵詞:元分析

    閆 偉 李 元 宋茂興 張曠野 孫銘澤 瞿 會 李鳳海鐘雪梅 朱 敏 杜萬里 呂香玲

    沈陽農(nóng)業(yè)大學(xué)特種玉米研究所, 遼寧沈陽 110866

    ?

    玉米抗灰斑病QTL元分析及其驗證

    閆 偉 李 元 宋茂興 張曠野 孫銘澤 瞿 會 李鳳海鐘雪梅 朱 敏 杜萬里 呂香玲*

    沈陽農(nóng)業(yè)大學(xué)特種玉米研究所, 遼寧沈陽 110866

    摘 要:玉米灰斑病是危害玉米生產(chǎn)的主要病害之一, 目前對抗灰斑病基因數(shù)目、位置及作用方式仍然不清楚, 這嚴重制約著玉米抗灰斑病育種進展。本研究利用元分析方法分析并整理了14篇玉米抗灰斑病QTL文獻的信息, 共篩選確定了13個一致性QTL區(qū)間。利用以自交系81162為輪回親本、自交系CN165為非輪回親本構(gòu)建的回交導(dǎo)入群體根據(jù)連鎖不平衡原理對13個一致性QTL進行驗證, 在13個一致性QTL區(qū)間共獲得20多個偏分離位點。第1和第4染色體上偏分離最嚴重, 其他染色體上偏分離度較小。說明第1和第4染色體上存在著效應(yīng)較大的抗病QTL。第1染色體標(biāo)記umc2227、bnlg1832、umc1243、umc2025、umc1515、umc1297、umc1461處供體基因頻率均在50%以上, 可能存在幾個連鎖的抗病基因。第4染色體上基因位于標(biāo)記bnlg2291和umc1194之間。研究為精細定位供體CN165中第1和第4染色體上的抗灰斑病QTL奠定了基礎(chǔ)。

    關(guān)鍵詞:灰斑病; 元分析; 一致性QTL; 回交導(dǎo)入系; 主效QTL

    本研究由國家自然科學(xué)基金項目(31301322)和遼寧省科技特派項目(2014215031)資助。

    This study was supported by the National Natural Science Foundation of China (31301322) and the Special Talent Appointment in Technology of Liaoning Province (2014215031).

    第一作者聯(lián)系方式∶ E-mail∶ yanwei1452@hotmail.com, Tel∶ 18512468548

    URL∶ http∶//www.cnki.net/kcms/detail/11.1809.S.20160314.1444.006.html

    玉米灰斑病是由尾孢菌侵染玉米葉部引起的真菌病害, 又被稱為玉米尾孢菌葉斑?。?-2]。玉蜀黍尾孢菌(Cecrospora zeae-maydis Tehon & Daniels)為玉米灰斑病的真菌侵染病原, 主要殘存在土壤中的病植物殘體上越冬, 其分生孢子通過風(fēng)和雨水飛濺散布到玉米植株上。溫和潮濕的條件下灰斑病會進一步侵染, 導(dǎo)致嚴重的葉片衰老和莖稈腐爛。于20世紀(jì)20年代在美國首次發(fā)現(xiàn)灰斑病, 如今已成為危害玉米生產(chǎn)的主要病害[3]。長期的生產(chǎn)實踐證明, 選育和推廣抗病品種是防治玉米灰斑病的最有效的途徑之一[4-5]。

    研究表明, 玉米對灰斑病的抗性屬多基因控制的數(shù)量性狀, 其不同的基因作用方式和效應(yīng)不同, 抗性基因的累積可使品種產(chǎn)生較高的抗性。因此玉米抗灰斑病主效QTL的發(fā)掘?qū)τ衩卓共∮N具有重要意義[6]。自20世紀(jì)90年代以來, 國內(nèi)外學(xué)者開展了一些玉米抗灰斑病的基因定位研究, 發(fā)掘出數(shù)量較多的抗病QTL[7-20]。但由于不同研究間存在材料遺傳背景、發(fā)病環(huán)境條件、圖譜標(biāo)記密度等差異, 定位結(jié)果間一致性較差。目前對抗灰斑病基因數(shù)目、位置及作用方式仍然不清楚, 這嚴重制約著玉米抗灰斑病育種進展。

    元分析(Meta-analysis)統(tǒng)計方法是對眾多現(xiàn)有實證文獻的再次統(tǒng)計, 通過將不同文獻的QTL位點統(tǒng)一投射到參考圖譜上, 根據(jù)每個QTL的中心位置及其置信區(qū)間, 即可以根據(jù)最大概率的原理估計出真實QTL的位置、縮小QTL的置信區(qū)間。目前, 高密度的玉米IBM遺傳圖譜為不同試驗QTL整合和一致性圖譜構(gòu)建奠定了基礎(chǔ)[21]。Chardon等[22]提出利用overview和元分析的方法對玉米花期相關(guān)的313個QTL信息進行整合分析, 考慮QTL定位中作圖群體及其大小、QTL效應(yīng)等因素, 最終獲得62個一致性QTL。Truntzler等[23]利用元分析整合了11個不同材料的與青貯玉米品質(zhì)相關(guān)的QTL原始圖譜, 得到一致性QTL。這種基于性狀定位的比較研究不僅有助于發(fā)掘一致性主效基因位點及其連鎖標(biāo)記, 同時為基因/QTL的克隆和利用提供了新手段[24]。

    雜交、回交與自交方法相結(jié)合, 構(gòu)建高世代回交導(dǎo)入系(或近等基因系)基因定位的理論, 成為目前開展性狀理論研究和育種實踐相結(jié)合的一個有效途徑。因為構(gòu)建的回交導(dǎo)入系, 其每一導(dǎo)入系與導(dǎo)入系, 或與輪回親本僅在某幾個染色體片段上有差異, 這樣的植株適用于性狀分析,可作為目標(biāo)性狀遺傳與生理學(xué)基礎(chǔ)研究分析的材料; 同時導(dǎo)入系在遺傳背景上繼承了輪回親本的主要優(yōu)良基因,基本保持了原有輪回親本的一般配合力與特殊配合力,可以快速應(yīng)用于育種[25]。如徐建龍等[26]、Li等[27]、李芳等[28]、Xie等[29]分別利用回交導(dǎo)入系進行了水稻抗旱、抗病及產(chǎn)量品質(zhì)相關(guān)的QTL研究; Ho等[30]、呂香玲等[31]、Li等[32]、Salvi等[33]、Pea等[34]、Teng等[35]分別利用導(dǎo)入系進行了玉米抗旱、抗病和產(chǎn)量相關(guān)QTL的定位研究。

    本文利用元分析方法對不同玉米群體的抗灰斑病QTL進行整合分析, 并對分析結(jié)果進行有條件的篩選,確定玉米抗灰斑病一致性 QTL。然后利用回交導(dǎo)入系群體進行一致性QTL驗證, 發(fā)掘玉米抗灰斑病主效QTL。

    1 材料與方法

    1.1 玉米抗灰斑病QTL文獻下載及格式文件整理

    從NCBI、SpringerLink等網(wǎng)站下載已發(fā)表的玉米抗灰斑病QTL相關(guān)文獻, 根據(jù)元分析軟件格式要求按照QTL名稱、染色體位置、置信區(qū)間、臨近標(biāo)識、群體類型等信息整理數(shù)據(jù)。

    1.2 玉米抗灰斑病QTL信息整合及元分析

    通過閱讀文獻, 將整理好的每個QTL的相關(guān)信息按一定格式載入BioMercator V4.2.1軟件中, 以IBM2 2008 Neighbors作為參考圖譜對所有QTL整合并繪制成整合圖譜。根據(jù)整合結(jié)果, 分別對第1、第2、第3、……、第10染色體進行元分析, 選取模型AIC值最小的為最適合模型,以包括2個及以上群體和3個及以上QTL為篩選條件, 確定一致性QTL的最可能位置及置信區(qū)間。

    1.3 利用回交導(dǎo)入系進行一致性QTL驗證

    1.3.1 供試材料 以CN165為供體和81162為受體通過2代回交和2代自交構(gòu)建出BC2F3群體, 自交系CN165為熱帶血緣材料, 高抗多種玉米病害, 植株屬于平展型, 較高大; 81162是吉林農(nóng)業(yè)科學(xué)院選育的優(yōu)良自交系, 株型緊湊, 葉片夾角小, 但對玉米灰斑病敏感。

    1.3.2 灰斑病抗性鑒定 2012年在沈陽農(nóng)業(yè)大學(xué)試驗田種植親本自交系和BC2F3群體, 親本各種植2行, BC2F3群體種植60行, 行長4 m, 每行17株, 利用人工接種技術(shù)在該群體內(nèi)篩選抗玉米灰斑病單株并自交留種, 每一果穗構(gòu)建為一個回交導(dǎo)入系。2013年種植親本自交系和回交導(dǎo)入系并人工接種鑒定。采用完全隨機區(qū)組排列, 單行區(qū)種植每個回交導(dǎo)入系, 行長4 m, 每行17株, 2次重復(fù)。鑒定用的玉米灰斑病菌是從沈陽地區(qū)采集的玉米感病葉片上的病斑, 按常規(guī)分離法分離培養(yǎng)的純培養(yǎng)物。采用玉米葉粉碳酸鈣瓊脂培養(yǎng)基(MLPCA)進行分生孢子的培養(yǎng)。用無菌水配制成接種用分生孢子懸液, 將其濃度調(diào)至2.5×103個 mL-1。在玉米植株喇叭口期(第9~第11葉期),用噴嘴處裝有20 mL注射器針頭的手提式注射器, 從植株喇叭口處平行插入, 將每株10 mL的病菌孢子懸液注入植株心葉, 接種所有單株。

    按抗病性調(diào)查分級方法[36], 于玉米乳熟期田間目測調(diào)查接種部位上、下葉片的發(fā)病情況, 記載病情級別。當(dāng)鑒定圃中的感病對照丹340達到其相應(yīng)的感病程度(7級)以上時, 鑒定則視為有效。1級為葉片上無病斑或僅有零星病斑, 病斑占葉面積少于或等于5%, 屬高抗(HR); 3級為葉片上有少量病斑, 占葉面積6%~10%, 屬抗(R); 5級為葉片上病斑較多, 占葉面積11%~30%, 屬中抗(MR); 7級為葉片上有大量病斑, 病斑相連, 占葉面積31%~70%,屬感?。⊿); 9級為葉片基本被病斑覆蓋, 葉片枯死, 屬高感(HS)。

    1.3.3 回交導(dǎo)入系的基因型 采用CTAB法提取54個抗玉米灰斑病單株基因組DNA。在MaizeGDB數(shù)據(jù)庫中根據(jù)獲得的一致性QTL區(qū)間所在的染色體區(qū)域查找SSR標(biāo)記信息, 下載引物序列, 由生工生物工程(上海)股份有限公司合成引物。SSR標(biāo)記分析包括PCR擴增、聚丙烯酰胺凝膠電泳、銀染等程序, 參考LABORATARY PROTOCOLS分析方法[37], 并根據(jù)經(jīng)驗優(yōu)化檢測技術(shù)。首先在親本間進行多態(tài)性篩選, 獲得具有多態(tài)性的引物分析導(dǎo)入系基因型。在同等遷移速率下, 用“A”表示后代單株擴增條帶與親本CN165條帶相一致, “B”表示與親本81162條帶相一致,“H”表示雜合條帶, “U”表示缺失條帶。

    1.3.4 主效QTL分析方法 根據(jù)連鎖不平衡理論分析電泳檢測結(jié)果, 當(dāng)所選單株的某一標(biāo)記位點基因型頻率與理論期望值之間存在顯著差異時, 認為存在該選擇位點。通常選擇P≤0.0001的顯著水平來減小誤差。其中供體位點的基因頻率高, 表明存在對供體基因的有利選擇,反之則不利; 同時可以推斷在高頻率基因型的標(biāo)記處存在性狀相關(guān)的QTL。當(dāng)這種偏分離相關(guān)性是由供體純合子引起時, 表明基因間存在著加性效應(yīng); 如果偏分離是由供體雜合子引起, 表明基因間存在著超顯性效應(yīng); 若偏分離是由供體純合子和雜合子共同引起, 則表明存在著部分或者完全顯性。

    2 結(jié)果與分析

    表1 玉米抗灰斑病QTL的文獻信息Table 1 Reference informations of resistance to gray leaf spot of corn

    2.1 玉米抗灰斑病QTL文獻及整理

    通過瀏覽NCBI等網(wǎng)站, 本文整理了14篇玉米抗灰斑病QTL文獻的信息, 共得到98個QTL區(qū)間(表1)。

    2.2 玉米抗灰斑病QTL的整合與一致性QTL確定

    整合已經(jīng)定位的98個QTL。表明, 在玉米的10條染色體上均有玉米抗灰斑病QTL存在, 每條染色體上分布的數(shù)量不同, 主要集中在第1、第2、第4、第5、第7、第8、第9染色體上。第6染色體上整合的QTL數(shù)最少(3個QTL),第4條染色體上整合的QTL數(shù)最多(21個QTL)。各QTL整合到IBM2 2008 Neighbors參考圖譜后的區(qū)間大小差異較大,從最小的1.7 cM到最大的277.83 cM。

    對98個抗灰斑病QTL進行元分析。在軟件給出的5個分析模型中, 選擇AIC值最小的模型, 共得到34個QTL區(qū)間(表2)。本文以群體數(shù)≥2、整合QTL數(shù)≥3為篩選條件,最終確定13個一致性QTL。這些QTL分別位于第1染色體3個, 置信區(qū)間分別為405.29~409.74 cM、525.28~538.44 cM 和927.07~973.79 cM; 第2染色體4個, 置信區(qū)間分別為280.51~288.84 cM、372.29~383.96 cM、468.84~486.22 cM 和624.34~635.58 cM; 第4染色體4個, 置信區(qū)間分別為175.09~187.61 cM、288.88~292.71 cM、433.84~444.70 cM和524.51~525.48 cM; 第7染色體1個, 置信區(qū)間為186.17~259.98 cM;第8染色體1個, 置信區(qū)間為378.64~390.84 cM (圖1)。

    2.3 回交導(dǎo)入系抗病性表現(xiàn)

    2012年對BC2F3群體單株接種鑒定, 共獲得54個抗病表現(xiàn)3級以上的抗病單株并自交留種獲得其回交導(dǎo)入系。2013年對54個回交導(dǎo)入系重復(fù)鑒定, 從每一回交導(dǎo)入系每一重復(fù)選取連續(xù)10株調(diào)查病斑覆蓋度, 并由10個單株病斑覆蓋度平均水平確定該導(dǎo)入系的抗感級別。對一個導(dǎo)入系, 選取2次重復(fù)鑒定結(jié)果中發(fā)病級別較高的結(jié)果作為導(dǎo)入系的最終抗感級別。54個回交導(dǎo)入系的抗病表現(xiàn)鑒定結(jié)果見表3。其中抗病性表現(xiàn)為1級的家系有21個, 3級的家系有19個, 5級的家系有14個。

    2.4 標(biāo)記位點基因型分析

    根據(jù)獲得的13個一致性QTL結(jié)果及玉米MaizeGDB信息選取一致性 QTL區(qū)間的標(biāo)記進行親本多態(tài)性篩選,獲得74對在親本間表現(xiàn)多態(tài)性的引物。利用該74對親本多態(tài)性引物對 54個抗病單株進行標(biāo)記基因型分析, 并計算供體基因頻率。供體基因頻率最高達89%, 其中達50%以上的標(biāo)記有15個, 分別位于第1和第4染色體上(表4)。

    表2 玉米抗灰斑病QTL的元分析結(jié)果Table 2 QTL meta-analysis of resistance to gray leaf spot of corn

    2.5 主效QTL分析

    本試驗54個導(dǎo)入系DNA采自BC2F3世代, 根據(jù)計算可知, 在非選擇的情況下, 群體中某一單株在單一標(biāo)記位點上基因型的理論概率P(AA)∶P(Aa)∶P(aa)=3∶2∶27, 分別代表了來源于供體親本的基因純合、基因雜合和來源于輪回親本的基因純合的理論概率。根據(jù)此理論值,計算3種基因型的χ2值, 利用χ2值的顯著性檢驗判斷標(biāo)記位點附近是否存在抗病性QTL。

    表4表明, 74個標(biāo)記位點中20多個位點出現(xiàn)了嚴重偏分離現(xiàn)象。其中第 1染色體標(biāo)記位點 χ2值最大達到121.46, 第 2染色體達64.7, 第 4染色體達213.98, 第 7染色體達18.98, 第8染色體達24.76。由此說明第1和第4染色體上存在主效抗病 QTL, 第 2染色體上抗病 QTL效應(yīng)次之, 其他染色體上的抗病QTL效應(yīng)可能較小。

    表3 54個導(dǎo)入系抗病性表現(xiàn)Table 3 Resistance performance of 54 backcross introgression lines

    表4 標(biāo)記位點供體基因頻率及卡平方值Table 4 Frequency of donor parent's genotype and χ2

    (續(xù)表4)

    (續(xù)表4)

    圖2 卡平方顯著性檢驗Fig.2 Significance testing by χ2test

    根據(jù)IBM2 2008 Neighbors標(biāo)記順序, 分析染色體上χ2值變化。由圖2可見, 第1染色體上χ2值沒有顯著的一個高峰點, 而是在標(biāo)記 umc2227、bnlg1832、umc1243、umc2025、umc1515、umc1297、umc1461處出現(xiàn)一個平緩的峰線。第2染色體上χ2高峰點位于標(biāo)記phi435417處。第4染色體上χ2高峰點位于標(biāo)記bnlg2291和umc1194處。第7染色體上χ2高峰點位于標(biāo)記umc2333處。第8染色體上χ2高峰點位于標(biāo)記phi080處。這說明在這些位點附近可能存在抗病QTL。

    3 討論

    迄今, 已有較多玉米抗灰斑病 QTL定位的研究, 然而仍然未能找到通用性標(biāo)記被育種家應(yīng)用于輔助育種。近年來, 元分析和生物信息學(xué)方法在整合不同 QTL信息并進一步挖掘候選基因上得到了廣泛的應(yīng)用。這一基于性狀的比較定位途徑可以減少分子育種的基本投入, 為開發(fā)復(fù)雜數(shù)量性狀的主效基因及標(biāo)記, 進而為克隆基因/QTL提供了技術(shù)支持。本研究通過對前人研究的總結(jié)及元分析處理, 共得到13個一致性QTL。利用回交導(dǎo)入群體根據(jù)連鎖不平衡原理對13個一致性 QTL進行驗證, 在13個一致性QTL區(qū)間共獲得20多個偏分離位點。第1和第4染色體上偏分離最嚴重, 其他染色體上偏分離度較少。說明第1和第4染色體上存在著效應(yīng)較大的抗病QTL。從χ2高峰點來看, 第 1染色體比較復(fù)雜, 標(biāo)記 umc2227、bnlg1832、umc1243、umc2025、umc1515、umc1297、umc1461 處 χ2值基本接近, 均較高, 供體基因頻率也均在 50%以上,這說明第1染色體抗病基因表現(xiàn)十分復(fù)雜, 可能存在幾個緊密連鎖的抗病基因。這也解釋了為什么多數(shù)研究在第1染色體上均定位到QTL, 但位置差異較大。第4染色體上χ2值最高, 而且供體基因頻率達到最大值89%, 這說明在本研究供體親本中第 4染色體上存在一個效應(yīng)較大的抗病基因。在前人的研究中, 第 4條染色體上獲得的 QTL數(shù)最多(21個)。如Bubeck等[19]發(fā)現(xiàn)在10條染色體上均存在抗病QTL; Saghai Maroof等[12]發(fā)現(xiàn)在第1、第2、第4、第8、第10染色體上有QTL; Gordon等[8,38]發(fā)現(xiàn)抗病QTL主要分布于第2和第4染色體的長臂上; Zwonitzer等[16]定位的QTL分別位于染色體bins 1.05~1.06,1.08~1.09, 2.02~2.03, 3.04~3.05, 4.05, 7.02, 8.07~8.09,10.05; Benson等[39]發(fā)現(xiàn)在第1、第2、第4染色體上存在QTL。綜合玉米抗灰斑病基因(QTL)定位研究表明, 玉米10條染色體上均可能存在抗病 QTL, 而且不同材料控制灰斑病抗性的 QTL數(shù)目與位置不同, 多數(shù)研究都含有2~3個主效QTL。如Gordon等[8]分別在染色體bin2.09和bin 4.08定位到2個主效QTL; Lehmensiek等[9]研究發(fā)現(xiàn)染色體 bin1.04和 bin2.04各存在一個效應(yīng)值較高的主效QTL; Veiga等[13]得到位于bin4.03、bin4.04和bin4.05的3個主效QTL; Chung等[15]定位出2個主效QTL, 分別位于bin1.04和bin8.05; Zwonitzer等[16]檢測到3個主效QTL,分別位于bin1.05、bin1.08和bin2.02; Asea等[17]分別在染色體bin2.09和bin4.08檢測到主效QTL; Balint等[18]定位到3個與抗病有關(guān)的主效QTL, 分別位于bin1.05、bin2.04 和bin4.05。本研究中供體CN165中至少含有2個以上主效抗病QTL, 它們分別位于第1和第4染色體。

    References

    [1]王桂清, 陳捷.玉米灰斑病抗病性研究進展.沈陽農(nóng)業(yè)大學(xué)學(xué)報, 2000, 31∶ 418-422 Wang G Q, Chen J.The research progress of resistance to gray leaf spot of maize.J Shenyang Agric Univ, 2000, 31∶ 418-422 (in Chinese with English abstract)

    [2]Wang J, Levy M, Dunkle L D.Sibling species of Cercospora associated with gray leaf spot of maize.Phytopathology, 1998, 88∶1269-1275

    [3]Ward J M J, Stromberg E L, Nowell D C, Nutter F W J.Gray leaf spot∶a disease of global importance in maize production.Plant Dis, 1999, 83∶ 884-895

    [4]李世強, 陳威, 譚靜, 王坤, 番興明.玉米灰斑病的抗性機理研究進展.生物技術(shù)進展, 2011, (1)∶ 112-115 Li S Q, Chen W, Tan J, Wang K, Fan X M.The progress of resistance mechanism on maize gray leaf spot.Curr Biotechnol, 2011,(1)∶ 112-115 (in Chinese with English abstract)

    [5]曹國輝.玉米灰斑病及抗性研究.玉米科學(xué), 2009, 17∶152-155 Cao G H.The research advance on resistance to grey leaf spot in maize.J Maize Sci, 2009, 17∶ 152-155 (in Chinese with English abstract)

    [6]呂香玲, 李新海, 陳陽, 史振聲, 李鳳海, 傅俊范.玉米種質(zhì)抗灰斑病鑒定與評價.玉米科學(xué), 2011, 19(6)∶ 125-128 Lü X L, Li X H, Chen Y, Shi Z S, Li F H, Fu J F.Evaluation and identification of resistance to gray leaf spot (gls) in maize germplasm.J Maize Sci, 2011, 19(6)∶ 125-128 (in Chinese with English abstract)

    [7]Zhang Y, Xu L, Fan X M, Tan J, Chen W.QTL mapping of resistance to gray leaf spot in maize.Theor Appl Genet, 2012, 125∶1797-1808

    [8]Gordon S G, Bartsch M, Matthies I, Gevers H O, Lipps P E.Linkage of molecular markers to Cercosporazeae-maydis resistance in maize.Crop Sci, 2004, 44∶ 628-636

    [9]Lehmensiek A, Esterhuizen A M, Staden D V, Nelson S W, Retief A E.Genetic mapping of gray leaf spot(GLS) resistance genes in maize.Theor Appl Genet, 2001, 103∶ 797-803

    [10]Danson J, Lagat M, Kimani M, Kuria A.Quantitative trait loci (QTLs) for resistance to gray leaf spot and common rust diseases of maize.Afr J Biotechnol, 2008, 7∶ 3247-3254

    [11]Clements M J, Dudley J W, White D G.Quantitative trait loci associated with resistance to gray leaf spot of corn.Phytopathology, 2000, 90∶ 1018-1025

    [12]Saghai Maroof M A, Yue Y G, Xiang Z X, Stromberg E L,Rufener G K.Identification of quantitative trait loci controlling resistance to gray leaf spot disease in maize.Theor Appl Genet,1996, 93∶ 539-546

    [13]Veiga A D, Pinho R G V, Resende L V.Quantitative trait loci associated with resistance to gray leaf spot and grain yield in corn.Ciênc Agrotec, 2012, 36∶ 31-38

    [14]Pozar G, Butruille D, Silva H D, McCuddin Z P, Penna J C V.Mapping and validation of quantitative trait loci for resistance to Cercosporazeae-maydis infection in tropical maize.Theor Appl Genet, 2009, 118∶ 553-564

    [15]Chung C L, Poland J, Kump K, Benson J, Longfellow J, Walsh E,Balint K P, Nelson R.Targeted discovery of quantitative trait loci for resistance to northern leaf blight and other diseases of maize.Theor Appl Genet, 2011, 123∶ 307-326

    [16]Zwonitzer J C, Coles N D, Krakowsky M D, Arellano C, Holland J B, McMullen M D, Pratt R C, Balint-Kurti P J.Mapping resistance quantitative trait Loci for three foliar diseases in a maize recombinant inbred line population-evidence for multiple disease resistance? Phytopathology, 2010, 100∶ 72-79

    [17]Asea G, Vivek B S, Lipps P E, Pratt R C.Genetic gain and cost efficiency of marker-assisted selection of maize for improved resistance to multiple foliar pathogens.Mol Breed, 2012, 29∶515-527

    [18]Balint K P J, Wisser R, Zwonitzer J C.Use of an advanced intercross line population for precise mapping of quantitative trait loci for gray leaf spot resistance in maize.Crop Sci, 2008, 48∶1696-1704

    [19]Bubeck D M, Goodman M M, Beavis W D, Grant D.Quantitative trait loci controlling resistance to gray leaf spot in maize.Crop Sci, 1993, 33∶ 838-847

    [20]Berger D K, Carstens M, Korsman J N, Middleton F, Kloppers F J, Tongoona P, Myburg A A.Mapping QTL conferring resistance in maize to gray leaf spot disease caused by Cercospora zeina.BMC Genet, 2014, 15∶ 60, DOI∶ 10.1186/1471-2156-15-60

    [21]Lee M, Sharopova N, Beavis W D, Grant D, Katt M, Blair D,Hallauer A.Expanding the genetic map of maize with the intermated B73 × Mo17 (IBM) population.Plant Mol Biol, 2002, 48∶453-461

    [22]Chardon F, Virlon B, Moreau L, Falque M, Joets J, Decousset L,Murigneux A, Charcosset A.Genetic architecture of flowering time in maize as inferred from quantitative trait loci meta-analysis and synteny conservation with the rice genome.Genetics, 2004, 168∶ 2169-2185

    [23]Truntzler M, Barrière Y, Sawkins M C, Lespinasse D, Betran J,Charcosset A, Moreau L.Meta-analysis of QTL involved in silage quality of maize and comparison with the position of candidate genes.Theor Appl Genet, 2010, 121∶ 1465-1482

    [24]Lü X L, Li X H, Xie C X, Hao Z F, Ji H L, Shi L Y, Zhang S H.Comparative QTL mapping of resistance to sugarcane mosaic virus in maize based on bioinformatics.Front Agric China, 2008,24∶ 365-371

    [25]Bernacchi D, Beck-Bunn T, Eshed Y, Lopez J, Petiard V, Uhlig J,Zamir D, Tanksley S.Identification of QTLs for traits of agronomic importance from Lycopersicon hirsutum.Theor Appl Genet, 1998, 97∶ 381-397

    [26]徐建龍, 薛慶中, 羅利軍, 黎志康.近等基因?qū)胂刀ㄎ凰究沟厩?shù)量性狀位點的研究初報.浙江農(nóng)業(yè)學(xué)報, 2002,14(1)∶ 14-19 Xu J L, Xue Q Z, Luo L J, Li Z K.Preliminary report on quantitative trait loci mapping of false smut resistance using near-isogenic introgression lines in rice.Acta Agric Zhejiangensis,2002, 14(1)∶ 14-19 (in Chinese with English abstract)

    [27]Li Z K, Fu B Y, Gao Y M, Xu J L, Ali J, Lafitte H R, Jiang Y Z,Domingo R J, Vijayakumar C H M, Maghirang R, Zheng T Q,Zhu L H.Genome-wide introgression lines and their use in genetic and molecular dissection of complex phenotypes in rice (Oryza sativa L.).Plant Mol Biol, 2005, 59∶ 33-52

    [28]李芳, 程立銳, 許美容, 徐建龍, 黎志康.利用品質(zhì)性狀的回交選擇導(dǎo)入系挖掘水稻抗紋枯病 QTL.作物學(xué)報, 2009, 35∶1729-1737 Li F, Cheng L R, Xu M R, Xu J L, Li Z K.QTL mining for sheath blight resistance using the backcross selected introgression lines for grain quality in rice.Acta Agron Sin, 2009, 35∶ 1729-1737 (in Chinese with English abstract)

    [29]Xie X B, Jin F X, Song M H, Suh J P, Hwang H G, Kim Y G,McCouch S R, Ahn S N.Fine mapping of a yield-enhancing QTL cluster associated with transgressive variation in an Oryza sativa × O.rufipogon cross.Theor Appl Genet, 2008, 116∶ 613-622

    [30]Ho J C, McCouch S R, Smith M E.Improvement of hybrid yield by advanced backcross QTL analysis in elite maize.Theor Appl Genet, 2002, 105∶ 440-448

    [31]呂香玲, 李新海, 郝轉(zhuǎn)芳, 吉海蓮, 史利玉, 張世煌.基于近等基因?qū)胂蛋l(fā)掘玉米抗甘蔗花葉病毒主效基因.玉米科學(xué),2007, 15(3)∶ 9-14 Lü X L, Li X H, Hao Z F, Ji H L, Shi L Y, Zhang S H.Identification of major QTL for resistance to sugarcane mosaic virus in maize based on introgression lines analysis.J Maize Sci, 2007,15(3)∶ 9-14 (in Chinese with English abstract)

    [32]Li Y, Liu C, Shi Y, Song Y, Wang T, Li Y.Genetic location and evaluation of chromosomal segments for drought tolerance at flowering stage in maize using selected backcross populations.Plant Breed, 2009, 128∶ 342-347

    [33]Salvi S, Corneti S, Bellotti M, Carraro N, Sanguineti M C, Castelletti S, Tuberosa R.Genetic dissection of maize phenology using an intraspecific introgression library.BMC Plant Biol, 2011,11∶ 4

    [34]Pea G, Aung H H, Frascaroli E, Landi P, Pè M E.Extensive genomic characterization of a set of near-isogenic lines for heterotic QTL in maize (Zea mays L.).BMC Genomics, 2013, 14∶ 61

    [35]Teng F, Zhai L, Liu R, Bai W, Wang L, Huo D, Tao Y, Zheng Y,Zhang Z.ZmGA3ox2, a candidate gene for a major QTL, qPH3.1,for plant height in maize.Plant J, 2013, 73∶ 405-416

    [36]董懷玉, 姜鈺, 王麗娟, 徐秀德.玉米種質(zhì)資源抗灰斑病鑒定與評價.植物遺傳資源學(xué)報, 2005, 6∶ 441-443 Dong H Y, Jiang Y, Wang L J, Xu X D.Evaluation on maize germplasm resources for resistance to gray leaf spot.J Plant Genet Resour, 2005, 6∶ 441-443 (in Chinese with English abstract)

    [37]CIMMYT.Applied Molecular Genetics Laboratory.Laboratory Protocol.Mexico∶ CIMMYT, 1998

    [38]Gordon S G, Lipps P E, Pratt R C.Heritability and components of resistance to Cercospora zeae-maydis derived from maize inbred VO613Y.Phytopathology, 2006, 96∶ 593-598

    [39]Benson J M, Poland J A, Benson B M, Stromberg E L, Nelson R J.Resistance to gray leaf spot of maize∶ genetic architecture and mechanisms elucidated through nested association mapping and near-isogenic line analysis.PLoS Genet, 2015, 11∶e1005045

    DOI:10.3724/SP.J.1006.2016.00758

    *通訊作者(

    Corresponding author)∶ 呂香玲, E-mail∶ lvxiangling521@126.com, Tel∶ 13664105185

    收稿日期Received()∶ 2015-10-15; Accepted(接受日期)∶ 2016-03-02; Published online(網(wǎng)絡(luò)出版日期)∶ 2016-03-14.

    Meta-analysis and Validation of QTL for Resistance to Gray Leaf Spot in Maize

    YAN Wei, LI Yuan, SONG Mao-Xing, ZHANG Kuang-Ye, SUN Ming-Ze, QU Hui, LI Feng-Hai, ZHONG Xue-Mei, ZHU Min, DU Wan-Li, and Lü Xiang-Ling*
    Special Maize Institute, Shenyang Agricultural University, Shenyang 110866, China

    Abstract:Gray leaf spot (GLS) is one of the most severe leaf diseases of maize worldwide.The breeding progress for resistance to GLS has been seriously hindered by less knowledge about QTL number, QTL intervals and mechanism of GLS.Based on meta-analysis, we conformed 13 consensus QTL regions from 14 articles on resistance to gray leaf spot of maize.One backcross with inbred line 81162 as recurrent parent and inbred line CN165 as donor parent, was used to test the consensus of those 13 QTL regions on the basic of linkage disequilibrium getting more than 20 partial separated loci.High level of partial separation indicated that there were QTLs with high effects of resistance to GLS on chromosome 1 and chromosome 4.On chromosome 1, the donor genes closed to markers of umc2227, bnlg1832, umc1243, umc2025, umc1515, umc1297, and umc1461 showed the frequency over 50%.Therefore, we inferred that there were a few of highly linked QTLs on chromosome 1.The resistant QTL on chromosome 4 was located between markers bnlg2291 and umc1194.Consequently, this study could lay a foundation for the QTL fine mapping on chromosome 1 and chromosome 4 in the donor parent CN165.

    Keywords:Gray leaf spot; Meta-analysis; Consensus QTLs; Backcross introgression lines; Major QTL

    猜你喜歡
    元分析
    基于元分析的管理創(chuàng)新與企業(yè)績效的關(guān)系研究
    軟科學(xué)(2017年2期)2017-04-01 21:34:31
    高等職業(yè)教育領(lǐng)域中產(chǎn)教融合研究的元分析
    職教論壇(2017年3期)2017-03-07 22:45:59
    樂觀—悲觀傾向與主觀幸福感關(guān)系的元分析
    工作記憶容量與第二語言水平關(guān)系的元分析
    護理實踐教學(xué)中在線學(xué)習(xí)效果的元分析
    信任性別差異的元分析
    制度環(huán)境與跨區(qū)域并購效應(yīng)關(guān)系研究
    會計之友(2016年8期)2016-04-20 15:13:54
    大學(xué)生主觀幸福感變遷的元分析研究
    信任關(guān)鍵變化量元分析及其對中國信任文化釋義
    語音記憶和視覺記憶在兒童漢語識字能力中同樣重要嗎?一個元分析綜述
    国产精品 国内视频| 丝袜美腿诱惑在线| 91成人精品电影| 国产精品国产av在线观看| 大片电影免费在线观看免费| 亚洲男人天堂网一区| 三上悠亚av全集在线观看| 久久精品国产综合久久久| 精品一区二区三区av网在线观看 | 国语对白做爰xxxⅹ性视频网站| 国产精品亚洲av一区麻豆| 亚洲欧美一区二区三区黑人| 久久精品国产综合久久久| 狂野欧美激情性xxxx| 男的添女的下面高潮视频| 精品人妻在线不人妻| 中文字幕制服av| 国精品久久久久久国模美| 欧美日韩福利视频一区二区| 一级毛片女人18水好多 | 黑人猛操日本美女一级片| 久久久久久久大尺度免费视频| 母亲3免费完整高清在线观看| 18禁观看日本| 国产亚洲欧美精品永久| 一本综合久久免费| 久久久国产一区二区| 精品国产一区二区久久| 国产伦人伦偷精品视频| 中文字幕人妻丝袜制服| 亚洲熟女毛片儿| 国产在线免费精品| 国产精品久久久av美女十八| 一边摸一边抽搐一进一出视频| 精品熟女少妇八av免费久了| 欧美成狂野欧美在线观看| 电影成人av| 中国美女看黄片| 久久久精品免费免费高清| 午夜影院在线不卡| 一级毛片女人18水好多 | 少妇人妻久久综合中文| 女性生殖器流出的白浆| 日韩,欧美,国产一区二区三区| 妹子高潮喷水视频| 五月天丁香电影| 久久综合国产亚洲精品| 欧美精品一区二区免费开放| 国产成人精品久久二区二区91| 在线 av 中文字幕| 国产爽快片一区二区三区| cao死你这个sao货| 不卡av一区二区三区| 国产一区二区 视频在线| 天天躁夜夜躁狠狠久久av| 18禁观看日本| 国产欧美日韩一区二区三区在线| 美女福利国产在线| 欧美人与性动交α欧美精品济南到| 男女之事视频高清在线观看 | 亚洲午夜精品一区,二区,三区| 丰满饥渴人妻一区二区三| 手机成人av网站| 我要看黄色一级片免费的| 国产人伦9x9x在线观看| 午夜久久久在线观看| 国产片内射在线| 日本vs欧美在线观看视频| 麻豆av在线久日| 免费在线观看完整版高清| 精品一品国产午夜福利视频| 精品亚洲乱码少妇综合久久| 免费在线观看影片大全网站 | 免费av中文字幕在线| 国产黄色免费在线视频| 电影成人av| 视频区欧美日本亚洲| 蜜桃在线观看..| 国语对白做爰xxxⅹ性视频网站| 国产主播在线观看一区二区 | 免费久久久久久久精品成人欧美视频| 国产在视频线精品| 欧美日本中文国产一区发布| 在线观看人妻少妇| 永久免费av网站大全| 亚洲精品国产av成人精品| 日本一区二区免费在线视频| 亚洲成人国产一区在线观看 | 精品少妇内射三级| netflix在线观看网站| 国产高清videossex| 久久这里只有精品19| 久久久久久免费高清国产稀缺| 国产激情久久老熟女| 午夜福利在线免费观看网站| 日日爽夜夜爽网站| 国产成人系列免费观看| 久久99一区二区三区| 激情视频va一区二区三区| 激情五月婷婷亚洲| 国产成人系列免费观看| 亚洲欧美日韩高清在线视频 | 电影成人av| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲专区中文字幕在线| 亚洲,一卡二卡三卡| 丰满迷人的少妇在线观看| 国产熟女午夜一区二区三区| 国产一区二区三区综合在线观看| 日本欧美视频一区| 欧美精品亚洲一区二区| 久久久精品94久久精品| 亚洲免费av在线视频| 亚洲综合色网址| 中文字幕色久视频| 国产三级黄色录像| xxxhd国产人妻xxx| 十分钟在线观看高清视频www| 99国产精品免费福利视频| 成人亚洲欧美一区二区av| 男女免费视频国产| 叶爱在线成人免费视频播放| av线在线观看网站| 国产亚洲av高清不卡| 国产精品99久久99久久久不卡| 水蜜桃什么品种好| 999久久久国产精品视频| 一边摸一边抽搐一进一出视频| 国产免费一区二区三区四区乱码| 国产成人a∨麻豆精品| 青春草亚洲视频在线观看| 亚洲av电影在线观看一区二区三区| 欧美精品av麻豆av| 婷婷丁香在线五月| 一级毛片女人18水好多 | av国产精品久久久久影院| 国产日韩一区二区三区精品不卡| 亚洲 欧美一区二区三区| 老司机影院毛片| 欧美日本中文国产一区发布| 天天躁夜夜躁狠狠久久av| 久久久久久亚洲精品国产蜜桃av| 国产高清国产精品国产三级| 一区二区三区四区激情视频| 色综合欧美亚洲国产小说| 最黄视频免费看| 七月丁香在线播放| 亚洲色图综合在线观看| 9热在线视频观看99| 久久人人97超碰香蕉20202| 考比视频在线观看| 精品国产一区二区三区久久久樱花| 国产极品粉嫩免费观看在线| 欧美另类一区| 99国产精品一区二区三区| 人人妻,人人澡人人爽秒播 | 爱豆传媒免费全集在线观看| 久久毛片免费看一区二区三区| 国产高清视频在线播放一区 | 免费看十八禁软件| 大香蕉久久成人网| 十八禁人妻一区二区| 色网站视频免费| 亚洲av电影在线进入| 80岁老熟妇乱子伦牲交| 亚洲,欧美,日韩| 女人久久www免费人成看片| 大陆偷拍与自拍| 在线观看免费日韩欧美大片| 麻豆国产av国片精品| 国产极品粉嫩免费观看在线| 国产亚洲欧美精品永久| 各种免费的搞黄视频| 久久久精品国产亚洲av高清涩受| 成人国语在线视频| 久久久精品国产亚洲av高清涩受| 日韩大码丰满熟妇| 香蕉丝袜av| 亚洲国产日韩一区二区| 国产在线观看jvid| 亚洲国产精品国产精品| 久久精品国产a三级三级三级| 高清不卡的av网站| 免费一级毛片在线播放高清视频 | 咕卡用的链子| 咕卡用的链子| 夜夜骑夜夜射夜夜干| 19禁男女啪啪无遮挡网站| 97在线人人人人妻| a级片在线免费高清观看视频| 国产一卡二卡三卡精品| 亚洲国产看品久久| 久久国产精品大桥未久av| 亚洲人成电影观看| 777米奇影视久久| 看免费av毛片| 亚洲欧美一区二区三区国产| 国产高清国产精品国产三级| 老司机午夜十八禁免费视频| 超色免费av| 欧美日韩成人在线一区二区| 亚洲精品日韩在线中文字幕| a级毛片在线看网站| 男的添女的下面高潮视频| 色94色欧美一区二区| 亚洲成av片中文字幕在线观看| 美女午夜性视频免费| 久久久久国产精品人妻一区二区| 人人妻人人添人人爽欧美一区卜| 亚洲天堂av无毛| 看十八女毛片水多多多| 两个人免费观看高清视频| 午夜影院在线不卡| 我要看黄色一级片免费的| 国产成人免费无遮挡视频| 国产成人91sexporn| 精品人妻一区二区三区麻豆| 国产精品免费视频内射| 日日爽夜夜爽网站| 91精品伊人久久大香线蕉| 人人妻人人添人人爽欧美一区卜| 99久久99久久久精品蜜桃| 久久青草综合色| 国产亚洲午夜精品一区二区久久| 一本久久精品| 久久精品久久久久久久性| 亚洲 国产 在线| 在线观看一区二区三区激情| 三上悠亚av全集在线观看| 成在线人永久免费视频| 中文字幕人妻丝袜制服| 精品视频人人做人人爽| 80岁老熟妇乱子伦牲交| 国产又色又爽无遮挡免| 精品国产乱码久久久久久男人| av欧美777| 精品少妇黑人巨大在线播放| 人体艺术视频欧美日本| 精品一区二区三区四区五区乱码 | 高清欧美精品videossex| 狠狠婷婷综合久久久久久88av| av又黄又爽大尺度在线免费看| 久久久久久亚洲精品国产蜜桃av| 国产一区有黄有色的免费视频| 男女无遮挡免费网站观看| 中文字幕另类日韩欧美亚洲嫩草| 亚洲欧洲国产日韩| 日韩大码丰满熟妇| 国产一区亚洲一区在线观看| 大片免费播放器 马上看| 国产无遮挡羞羞视频在线观看| 欧美av亚洲av综合av国产av| 两性夫妻黄色片| 我要看黄色一级片免费的| 久久影院123| 丰满人妻熟妇乱又伦精品不卡| 精品国产超薄肉色丝袜足j| 亚洲av成人不卡在线观看播放网 | 老熟女久久久| 久久国产精品大桥未久av| 99久久精品国产亚洲精品| 这个男人来自地球电影免费观看| 日韩中文字幕视频在线看片| 中文字幕色久视频| 亚洲国产精品一区三区| 超碰97精品在线观看| 日韩,欧美,国产一区二区三区| 只有这里有精品99| 国产成人免费无遮挡视频| 国产亚洲精品第一综合不卡| 黄色片一级片一级黄色片| 亚洲专区国产一区二区| 精品免费久久久久久久清纯 | 成人亚洲精品一区在线观看| 欧美日韩成人在线一区二区| 成年人午夜在线观看视频| 国产在线免费精品| 久久99热这里只频精品6学生| 午夜免费观看性视频| 99热网站在线观看| av线在线观看网站| 99热国产这里只有精品6| 国产精品国产三级国产专区5o| 精品国产乱码久久久久久男人| av在线播放精品| 免费看不卡的av| 日韩大片免费观看网站| 黑人巨大精品欧美一区二区蜜桃| 国产成人影院久久av| 国产精品久久久久成人av| 国产一区二区三区综合在线观看| 国产精品一区二区免费欧美 | 色精品久久人妻99蜜桃| 五月天丁香电影| 午夜福利在线免费观看网站| 午夜福利视频精品| 欧美日本中文国产一区发布| 制服人妻中文乱码| 午夜福利一区二区在线看| 18禁黄网站禁片午夜丰满| 国产日韩欧美视频二区| 亚洲一区二区三区欧美精品| 欧美国产精品va在线观看不卡| 国产亚洲欧美在线一区二区| 中文字幕av电影在线播放| 久久av网站| 一边亲一边摸免费视频| 搡老岳熟女国产| 免费人妻精品一区二区三区视频| 亚洲欧美精品综合一区二区三区| 亚洲国产欧美一区二区综合| 手机成人av网站| 久久久精品94久久精品| h视频一区二区三区| 国产成人欧美在线观看 | 国产黄色免费在线视频| 97人妻天天添夜夜摸| 国产精品.久久久| 9色porny在线观看| 国产成人一区二区在线| 9色porny在线观看| 国产成人一区二区在线| 欧美精品av麻豆av| 国产91精品成人一区二区三区 | 免费不卡黄色视频| 极品少妇高潮喷水抽搐| 一级毛片 在线播放| 丝瓜视频免费看黄片| 国产色视频综合| 亚洲精品成人av观看孕妇| 大话2 男鬼变身卡| 亚洲精品日本国产第一区| 亚洲国产精品999| 精品国产乱码久久久久久小说| 亚洲国产欧美网| 精品亚洲乱码少妇综合久久| 欧美日本中文国产一区发布| 一级,二级,三级黄色视频| 亚洲精品久久久久久婷婷小说| 日本av免费视频播放| 免费少妇av软件| 制服人妻中文乱码| 好男人视频免费观看在线| 一级黄色大片毛片| 亚洲av在线观看美女高潮| 女人爽到高潮嗷嗷叫在线视频| 99国产精品免费福利视频| 黄色毛片三级朝国网站| 美女福利国产在线| 超碰97精品在线观看| 国产免费又黄又爽又色| 日韩大码丰满熟妇| 国产99久久九九免费精品| 国产日韩欧美在线精品| 99国产精品一区二区蜜桃av | 在线观看免费午夜福利视频| 午夜91福利影院| 成年人午夜在线观看视频| 在线观看一区二区三区激情| 亚洲激情五月婷婷啪啪| 精品少妇久久久久久888优播| 亚洲精品一区蜜桃| 欧美日韩综合久久久久久| bbb黄色大片| 欧美日韩国产mv在线观看视频| 飞空精品影院首页| 国产av精品麻豆| 别揉我奶头~嗯~啊~动态视频 | 丝瓜视频免费看黄片| 日韩av在线免费看完整版不卡| 国精品久久久久久国模美| 精品一品国产午夜福利视频| 精品久久久久久电影网| 天堂俺去俺来也www色官网| 亚洲一区二区三区欧美精品| 天堂俺去俺来也www色官网| 久久精品亚洲av国产电影网| 日韩一区二区三区影片| 国产午夜精品一二区理论片| 一本—道久久a久久精品蜜桃钙片| 欧美亚洲 丝袜 人妻 在线| 成人18禁高潮啪啪吃奶动态图| 亚洲av在线观看美女高潮| 视频在线观看一区二区三区| 精品久久久精品久久久| 国产在线观看jvid| 精品人妻熟女毛片av久久网站| 男女高潮啪啪啪动态图| 丁香六月天网| 一级毛片黄色毛片免费观看视频| 久久久久视频综合| 悠悠久久av| 国产色视频综合| 亚洲国产精品成人久久小说| 岛国毛片在线播放| 国产视频首页在线观看| 老司机影院成人| 精品久久久久久久毛片微露脸 | 男女之事视频高清在线观看 | 99久久人妻综合| 只有这里有精品99| 2021少妇久久久久久久久久久| 久久av网站| 国产精品亚洲av一区麻豆| 一区二区三区激情视频| 天天操日日干夜夜撸| 丰满迷人的少妇在线观看| 伊人久久大香线蕉亚洲五| 日韩av在线免费看完整版不卡| 亚洲精品日本国产第一区| 精品国产乱码久久久久久男人| 国产高清国产精品国产三级| 少妇粗大呻吟视频| 一级,二级,三级黄色视频| 真人做人爱边吃奶动态| 在线亚洲精品国产二区图片欧美| 人成视频在线观看免费观看| 又黄又粗又硬又大视频| 国产精品一区二区在线观看99| 国产日韩一区二区三区精品不卡| 欧美黑人精品巨大| 亚洲免费av在线视频| 欧美+亚洲+日韩+国产| 大香蕉久久成人网| 亚洲五月色婷婷综合| 欧美精品人与动牲交sv欧美| 少妇的丰满在线观看| 美女大奶头黄色视频| 欧美国产精品va在线观看不卡| 成在线人永久免费视频| 国精品久久久久久国模美| 久久av网站| 亚洲精品av麻豆狂野| 天天躁狠狠躁夜夜躁狠狠躁| 少妇的丰满在线观看| 久久精品熟女亚洲av麻豆精品| 伦理电影免费视频| 黄色视频不卡| 一级毛片电影观看| 国产精品国产av在线观看| 成人18禁高潮啪啪吃奶动态图| 国产欧美日韩精品亚洲av| 午夜久久久在线观看| 精品一品国产午夜福利视频| 青青草视频在线视频观看| 女性生殖器流出的白浆| 一本久久精品| 亚洲精品自拍成人| 欧美成人精品欧美一级黄| 亚洲色图 男人天堂 中文字幕| 亚洲精品在线美女| 老司机午夜十八禁免费视频| 天天躁狠狠躁夜夜躁狠狠躁| 久久久久精品国产欧美久久久 | 午夜福利,免费看| 亚洲精品久久午夜乱码| 久久中文字幕一级| 97精品久久久久久久久久精品| 一区二区三区精品91| 久久国产精品影院| 亚洲人成电影免费在线| 亚洲成人免费av在线播放| 久久久久视频综合| 18禁国产床啪视频网站| 久久99热这里只频精品6学生| av有码第一页| 色婷婷久久久亚洲欧美| 国产成人影院久久av| 在线亚洲精品国产二区图片欧美| 狂野欧美激情性xxxx| 免费看不卡的av| 亚洲精品久久午夜乱码| 肉色欧美久久久久久久蜜桃| 欧美成人午夜精品| 操出白浆在线播放| av在线老鸭窝| 国产激情久久老熟女| 欧美激情 高清一区二区三区| 国产成人av激情在线播放| av有码第一页| 亚洲一区中文字幕在线| 亚洲天堂av无毛| 国产亚洲av片在线观看秒播厂| 国产精品欧美亚洲77777| 久久精品国产亚洲av涩爱| 一级毛片我不卡| 色精品久久人妻99蜜桃| 国产爽快片一区二区三区| 亚洲国产av影院在线观看| 国产精品久久久av美女十八| 丰满迷人的少妇在线观看| videosex国产| 老司机在亚洲福利影院| 婷婷色综合大香蕉| 久久国产精品男人的天堂亚洲| 亚洲欧洲精品一区二区精品久久久| 啦啦啦 在线观看视频| 超碰成人久久| 亚洲人成电影观看| 在线亚洲精品国产二区图片欧美| 自拍欧美九色日韩亚洲蝌蚪91| 久久av网站| 天天躁日日躁夜夜躁夜夜| 人成视频在线观看免费观看| 免费看av在线观看网站| 在线av久久热| 赤兔流量卡办理| 久久久久精品国产欧美久久久 | 亚洲一码二码三码区别大吗| 婷婷成人精品国产| 久久久久国产一级毛片高清牌| 热99久久久久精品小说推荐| 欧美日韩亚洲高清精品| 欧美成人精品欧美一级黄| 午夜免费男女啪啪视频观看| 亚洲黑人精品在线| 午夜激情久久久久久久| 亚洲国产精品成人久久小说| 黄片小视频在线播放| 日韩av在线免费看完整版不卡| 精品国产一区二区久久| 精品久久久精品久久久| 国产野战对白在线观看| 女人久久www免费人成看片| 午夜福利在线免费观看网站| 久久中文字幕一级| 91国产中文字幕| 丰满人妻熟妇乱又伦精品不卡| 国精品久久久久久国模美| 国产欧美日韩一区二区三区在线| 亚洲人成网站在线观看播放| 青春草亚洲视频在线观看| 男女床上黄色一级片免费看| 国产成人精品久久二区二区免费| 美女大奶头黄色视频| 亚洲国产欧美日韩在线播放| 精品久久蜜臀av无| 中文字幕人妻丝袜一区二区| 无限看片的www在线观看| 色精品久久人妻99蜜桃| 久久人人97超碰香蕉20202| 2018国产大陆天天弄谢| 国产欧美亚洲国产| 男人爽女人下面视频在线观看| 免费人妻精品一区二区三区视频| 日本五十路高清| av网站在线播放免费| 人妻一区二区av| av电影中文网址| 欧美亚洲 丝袜 人妻 在线| 久久精品国产a三级三级三级| 亚洲综合色网址| 人人妻人人爽人人添夜夜欢视频| 国产免费一区二区三区四区乱码| 男的添女的下面高潮视频| 欧美亚洲日本最大视频资源| 日韩中文字幕视频在线看片| 亚洲第一青青草原| 美国免费a级毛片| 母亲3免费完整高清在线观看| 国产成人欧美| av网站在线播放免费| 在线天堂中文资源库| 中文字幕精品免费在线观看视频| 亚洲人成电影观看| 在线观看人妻少妇| 青青草视频在线视频观看| 狠狠婷婷综合久久久久久88av| 久久久久国产精品人妻一区二区| 免费在线观看黄色视频的| 日韩精品免费视频一区二区三区| 日日爽夜夜爽网站| 电影成人av| 纯流量卡能插随身wifi吗| 中文精品一卡2卡3卡4更新| 国产一区二区三区综合在线观看| 免费在线观看视频国产中文字幕亚洲 | 少妇的丰满在线观看| 另类亚洲欧美激情| 精品久久久精品久久久| 婷婷色麻豆天堂久久| 在线亚洲精品国产二区图片欧美| 亚洲国产中文字幕在线视频| 九草在线视频观看| 在线 av 中文字幕| 又大又黄又爽视频免费| 久久亚洲精品不卡| 99久久99久久久精品蜜桃| 视频在线观看一区二区三区| 日本91视频免费播放| 免费在线观看日本一区| 国产精品 国内视频| 99久久人妻综合| 久久热在线av| 男女之事视频高清在线观看 | 亚洲精品自拍成人| 丝袜喷水一区| 一级毛片女人18水好多 | 亚洲精品中文字幕在线视频| 久久久久国产一级毛片高清牌| 美女高潮到喷水免费观看| a级毛片黄视频| 国产国语露脸激情在线看| 熟女av电影| 国产老妇伦熟女老妇高清| 精品熟女少妇八av免费久了| 成人午夜精彩视频在线观看| 精品国产国语对白av| 欧美日韩视频精品一区| 国产淫语在线视频| 一区二区三区四区激情视频| 亚洲精品日韩在线中文字幕| 国产日韩欧美在线精品| 丝袜美足系列| 波多野结衣一区麻豆|