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

    玉米尾孢菌灰斑病抗性種質(zhì)鑒定及其抗性基因分析

    2020-02-22 02:54譚靜羅吉王文瑞王琨高佳琪
    關(guān)鍵詞:關(guān)聯(lián)分析玉米

    譚靜 羅吉 王文瑞 王琨 高佳琪

    摘要:灰斑病是全球玉米的重要病害,特別對中國西南地區(qū)玉米的生產(chǎn)造成了嚴(yán)重影響,玉米尾孢菌是中國西南地區(qū)玉米灰斑病的主要致病菌種。本研究采用來源廣泛的241份玉米自交系構(gòu)建關(guān)聯(lián)群體,在云南玉米灰斑病重病區(qū)進(jìn)行灰斑病的田間自然發(fā)病抗性鑒定;利用覆蓋玉米全基因組的單核苷酸多態(tài)性(SNP)標(biāo)記對由241份玉米自交系構(gòu)建的關(guān)聯(lián)群體進(jìn)行基因型分析,在此基礎(chǔ)上通過全基因組關(guān)聯(lián)分析進(jìn)行玉米灰斑病抗性基因的定位及其候選基因的預(yù)測,為玉米灰斑病抗性育種提供技術(shù)支撐。2017-2018年的玉米灰斑病抗性鑒定結(jié)果顯示,2年共有109份熱帶/亞熱帶玉米自交系之間和109份溫帶玉米自交系材料之間以及材料與年份互作間的差異均達(dá)極顯著水平;在2017-2018年的抗性鑒定結(jié)果中均表現(xiàn)為抗病的材料有44份,均表現(xiàn)為高抗的材料僅有5份。用20 586個(gè)高質(zhì)量SNP標(biāo)記、2年玉米灰斑病表型數(shù)據(jù)進(jìn)行全基因組關(guān)聯(lián)分析,結(jié)果顯示,在P≤0.001(-lgP≥3)條件下,共檢測到44個(gè)與玉米灰斑病抗性顯著關(guān)聯(lián)的SNP位點(diǎn),其中34個(gè)SNP位點(diǎn)對應(yīng)31個(gè)基因,分布在除10號染色體外的其他染色體上。結(jié)合已有的連鎖定位研究結(jié)果可知,在以上31個(gè)基因中,除了4個(gè)基因外,其余27個(gè)基因均分布在已知數(shù)量性狀座位(QTL)區(qū)間內(nèi),將這27個(gè)基因在生物信息數(shù)據(jù)庫中與已有信息進(jìn)行比對,共獲得8個(gè)與抗病相關(guān)的候選基因。研究結(jié)果可為由玉米尾孢菌引起的玉米灰斑病的抗性育種及其抗性基因的克隆、抗病機(jī)制的研究奠定基礎(chǔ)。

    關(guān)鍵詞:玉米;玉米尾孢菌;灰斑病;抗性種質(zhì);抗病基因;關(guān)聯(lián)分析

    中圖分類號:S435.131.4+9文獻(xiàn)標(biāo)識碼:A文章編號:1000-4440(2020)06-1373-09

    Abstract: Gray leaf spot (GLS) is an important disease of maize in the globe, which brings about severe impact on maize production especially in Southwestern China. Cercospora zeina is the main pathogen of GLS in Southwestern China. In this study, 241 maize inbred lines with wide sources were used to construct association population, and the resistance phenotype to GLS under natural condition in the severe disease area of Yunnan was evaluated. Genotypes of the association population constructed by 241 maize inbred lines were analyzed by single nucleotide polymorphism (SNP) markers covering the whole genome. Then the loci of GLS-resistant genes and the candidate genes were predicted to provide technical support for the GLS-resistant breeding of maize by genome-wide association analysis. The two-year (2017 to 2018) identification results of GLS-resistance indicated that there were significant differences among 109 tropical/subtropical maize inbred lines, 109 temperate maize inbred lines and the interaction between materials and years. There were 44 inbred lines with disease resistance and five inbred lines with high resistance in two years of resistant evaluation. The results of genome-wide association analysis by using 20 586 high-quality SNP markers and two-year phenotypic data of GLS indicated that, under the condition of P≤0.001 (-lgP≥3), a total of 44 SNP loci significantly associated with GLS-resistance were detected, among which 34 SNP loci corresponding to 31 genes were distributed on chromosomes except for chromosome 10. Compared with existing research results on linkage location, it could be seen that except for four genes, the remaining 27 genes of the 31 genes were located in the quantitative trait locus (QTL). Eight candidate genes related to diseaseresistance were obtained by comparing the 27 genes in bioinformatics database. The results will lay a foundation for the resistance breeding, the cloning of resistant genes and the resistance mechanism research on maize GLS (Cercospora zeina).

    Key words:maize;Cercospora zeina;gray leaf spot;resistant germplasm;resistant gene;correlation analysis

    玉米(Zea Mays L.)是全球第一大糧食作物,灰斑病是玉米葉部的主要病害,在世界各玉米產(chǎn)區(qū)均有發(fā)生,一般會造成玉米產(chǎn)量損失達(dá)5%~30%,嚴(yán)重時(shí)可達(dá)60%以上,甚至造成絕收[1-4]。在中國西南高海拔玉米種植區(qū),灰斑病已經(jīng)成為玉米生產(chǎn)中最重要的病害,在云南地區(qū)已經(jīng)超過玉米大斑病成為玉米生產(chǎn)中的第一大病害[2, 5-6]。長期以來,國內(nèi)研究者一致認(rèn)為,玉蜀黍尾孢菌是中國玉米灰斑病的致病菌種[7],直到2013年才有報(bào)道指出,玉米尾孢菌是中國西南地區(qū)玉米灰斑病的主要致病菌種[8-9]。

    培育和種植抗病品種是控制植物病害最經(jīng)濟(jì)有效的途徑,而抗病育種的成效取決于對獲得抗原、抗病遺傳機(jī)制的認(rèn)識。目前,國內(nèi)外研究者已經(jīng)開展了大量玉米灰斑病的抗源篩選工作,雖然有研究結(jié)果表明,玉米自交系、雜交種對灰斑病的抗性差異明顯,但均未發(fā)現(xiàn)對灰斑病免疫的玉米種質(zhì),大多數(shù)玉米種質(zhì)對灰斑病表現(xiàn)為感病[3,10-12]。多數(shù)研究發(fā)現(xiàn),玉米對灰斑病的抗性主要表現(xiàn)為加性遺傳效應(yīng),屬于數(shù)量性狀遺傳效應(yīng)[13-14]。研究者已經(jīng)用不同群體定位了多個(gè)抗玉米灰斑病QTL,這些QTL在玉米的10條染色體上均有分布[15-24]。但是這些研究大多是針對玉蜀黍尾孢菌玉米灰斑病的抗性,加上受連鎖分析自身的局限,不能全面挖掘玉米灰斑病抗性基因,且多數(shù)定位精度不高,難以在玉米灰斑病的抗性育種中得以應(yīng)用。

    關(guān)聯(lián)分析是一種以連鎖不平衡為基礎(chǔ),鑒定與表型變異緊密聯(lián)系的功能性等位基因/基因型的分析方法。具有研究周期短、分辨率高、可同時(shí)檢測數(shù)十個(gè)基因等優(yōu)點(diǎn),是復(fù)雜數(shù)量性狀基因挖掘和驗(yàn)證的首選方法[25-27]。Mammadov等[28]結(jié)合連鎖和關(guān)聯(lián)分析法,在玉米第1、6、7、8號染色體上獲得了抗玉米灰斑病QTL;Kuki等[29]通過全基因組關(guān)聯(lián)分析檢測到7個(gè)單核苷酸多態(tài)性(SNP)位點(diǎn)與玉米灰斑病抗性顯著相關(guān),并發(fā)現(xiàn)3個(gè)抗玉米灰斑病的候選基因。本研究利用自然發(fā)病原理對241份來源廣泛的玉米自交系進(jìn)行玉米尾孢菌灰斑病抗性鑒定,為西南地區(qū)玉米抗灰斑病育種提供材料支撐。同時(shí),本研究利用覆蓋全基因組的SNP標(biāo)記,通過全基因組關(guān)聯(lián)分析全面挖掘抗玉米尾孢菌灰斑病的基因及其位點(diǎn),并對檢測到的顯著性位點(diǎn)與前人研究得到的QTL區(qū)間進(jìn)行比較,結(jié)合生物信息學(xué)分析方法,從一致性位點(diǎn)中預(yù)測候選基因,為玉米尾孢菌灰斑病抗性基因的克隆和抗病機(jī)制的研究奠定基礎(chǔ)。

    1材料與方法

    1.1試驗(yàn)材料

    參考已有表型數(shù)據(jù)和系譜關(guān)系,從來源廣泛的630份玉米自交系中優(yōu)先選擇沒有直接親緣關(guān)系、抗灰斑病差異大的241份玉米自交系構(gòu)建關(guān)聯(lián)群體,其中111份來自CIMMYT(國際玉米小麥改良中心),30份來自美國,6份來自泰國,其余94份來自國內(nèi)不同玉米生態(tài)區(qū),其中包括在中國玉米產(chǎn)業(yè)中大面積推廣使用的一些重要雜交種的親本。

    1.2田間試驗(yàn)與抗性鑒定

    利用重病區(qū)自然發(fā)病原理進(jìn)行玉米灰斑病的抗性鑒定,于2017-2018年在云南保山市郊進(jìn)行。根據(jù)系譜將研究材料分為熱帶/亞熱帶、溫帶2組,組內(nèi)采用隨機(jī)區(qū)組設(shè)計(jì),設(shè)2次重復(fù),同時(shí)選取4個(gè)典型材料(抗病自交系Y32、感病自交系Ye478、抗病雜交種海禾2號、感病雜交種會單四號)作為對照,每40行設(shè)置1組對照(含2個(gè)抗病對照材料、2個(gè)感病對照材料),在試驗(yàn)地外圍種植2行感病材料作為保護(hù)行。種植規(guī)格參照當(dāng)?shù)亓?xí)慣,加強(qiáng)田間水肥管理,試驗(yàn)期間防蟲但不防病,在開花后第2周參照王曉鳴等[30]的玉米灰斑病分級標(biāo)準(zhǔn)對植株的抗病性進(jìn)行鑒定,按照抗病性由強(qiáng)到弱分為1級、3級、5級、7級、9級(其中1級為高抗,9級為高感),以后每隔1周調(diào)查1次,共調(diào)查3次,取病級最符合正態(tài)分布的1次調(diào)查結(jié)果。每次調(diào)查均由2人進(jìn)行獨(dú)立評價(jià),每個(gè)小區(qū)最終的病害級別為2人評價(jià)值的平均值。

    1.3基因型鑒定

    用Affymetrix SNP 56K芯片(包含56 110個(gè)SNP標(biāo)記)對上述自交系進(jìn)行基因型分析,對芯片檢測結(jié)果進(jìn)行嚴(yán)格的質(zhì)量控制,最終選取20 586個(gè)高質(zhì)量SNP標(biāo)記用于后續(xù)的數(shù)據(jù)分析。

    1.4數(shù)據(jù)分析

    用SPSS 25.0軟件對表型數(shù)據(jù)進(jìn)行描述性統(tǒng)計(jì)和方差分析,用Structure(V 2.3.4)軟件進(jìn)行群體結(jié)構(gòu)分析,再利用TASSEL 5.0軟件中的混合線性模型(Mixed-linear model,MLM),在考慮群體結(jié)構(gòu)和親緣關(guān)系的條件下進(jìn)行SNP標(biāo)記與玉米灰斑病表型之間的關(guān)聯(lián)分析,參照P≤0.001(-lgP≥3)的標(biāo)準(zhǔn)進(jìn)行顯著性位點(diǎn)檢測,并對檢測到的顯著性位點(diǎn)與前人研究得到的QTL區(qū)間進(jìn)行比較,結(jié)合生物信息學(xué)方法,預(yù)測與前人研究得到的QTL一致性位點(diǎn)中的候選基因。

    2結(jié)果與分析

    2.1抗性鑒定

    由于田間出苗、生長等因素而使部分材料缺失,2017年獲得了232份玉米自交系的灰斑病抗性鑒定結(jié)果,2018年獲得了225份玉米自交系的灰斑病抗性鑒定結(jié)果。2017年的玉米灰斑病抗性鑒定結(jié)果顯示,抗病自交系Y32的平均病級為1.50,抗病雜交種海禾2號的平均病級為2.19;感病自交系Ye478的平均病級為6.88,感病雜交種會單四號的平均病級為6.31。2018年的玉米灰斑病抗性鑒定結(jié)果顯示,抗病自交系Y32的平均病級為1.00,抗病雜交種海禾2號的平均病級為3.25;感病自交系Ye478的平均病級為7.63,感病雜交種會單四號的平均病級為7.94。由此可見,在鑒定試驗(yàn)中抗病對照、感病對照的發(fā)病病級適宜且二者間差異明顯,說明在本試驗(yàn)條件下得出的玉米灰斑病發(fā)病程度和鑒定效果可以用來區(qū)分不同玉米材料的抗性差異。

    對2017-2018年共有的218份玉米自交系的灰斑病抗性鑒定結(jié)果進(jìn)行方差分析,由表1可見,材料間以及年份與材料間互作的差異均達(dá)極顯著水平,說明玉米灰斑病表型與環(huán)境之間存在極顯著互作,109份熱帶/亞熱帶玉米自交系和109份溫帶玉米自交系對玉米灰斑病的抗性均存在極顯著差異。正態(tài)分布檢驗(yàn)結(jié)果表明,2017-2018年玉米灰斑病表型的偏度、峰度絕對值均小于1,符合正態(tài)分布,2017-2018年玉米灰斑病病級的頻率分布直方圖見圖1。

    在2017年,232份玉米自交系的灰斑病病級在1.00~9.00級均有分布,平均病級為4.79級,其中高感材料(病級>7.00)37份,感病材料(5.00<病級≤7.00)63份,中抗材料(3.00<病級≤5.00)59份,抗病材料(1.00<病級≤3.00)61份,高抗材料(病級≤1.00)12份。在2018年,225份玉米自交系的灰斑病病級在1.00~9.00級均有分布,平均病級為4.96級,其中高感材料33份,感病材料61份,中抗材料55份,抗病材料57份,高抗材料19份。綜合2017-2018年的玉米灰斑病抗性鑒定結(jié)果可知,2年都表現(xiàn)為抗病的材料有44份,2年都表現(xiàn)為高抗的材料有5份(分別為DTMA25、CIMBL32、CIMBL70、GEMS19、PS61-20),可以在玉米灰斑病抗病育種中應(yīng)用這5份高抗材料。

    2.2全基因組分析

    對241份玉米自交系的群體結(jié)構(gòu)分析可知,△K在K=2時(shí)出現(xiàn)第1個(gè)峰值(圖2a、圖2b),將241份玉米自交系劃分成2個(gè)亞群,與系譜信息進(jìn)行比對后發(fā)現(xiàn),第1個(gè)亞群對應(yīng)育種實(shí)踐中的熱帶/亞熱帶(Tropical and subtropical, TST)種質(zhì),第2個(gè)亞群對應(yīng)育種實(shí)踐中的溫帶(Temperate, TEM)種質(zhì);△K在K=4時(shí)出現(xiàn)第2個(gè)峰值(圖2a、圖2c),說明該群體能進(jìn)一步被劃分為4個(gè)亞群,對應(yīng)系譜中這4個(gè)亞群分別為TST群、溫帶堅(jiān)稈(Stiff stalk)群、溫帶非堅(jiān)稈(Non-stiff stalk)群和溫帶唐四平頭群。

    分析結(jié)果表明,在P≤0.001(-lgP≥3)條件下,2017年檢測到22個(gè)與抗性顯著關(guān)聯(lián)的SNP位點(diǎn),2018年檢測到23個(gè)與抗性顯著關(guān)聯(lián)的SNP位點(diǎn),2年均檢測到的SNP位點(diǎn)僅有1個(gè),對應(yīng)標(biāo)記為AX86259524(表2)。通過查閱注釋文件發(fā)現(xiàn),在2年內(nèi)檢測到的與玉米灰斑病抗性關(guān)聯(lián)的44個(gè)SNP位點(diǎn)中,有10個(gè)位點(diǎn)的注釋信息位于基因間區(qū)域,其余34個(gè)有注釋信息的SNP位點(diǎn)共對應(yīng)31個(gè)基因,這31個(gè)基因分布在除10號染色體外的其他染色體上,其中分布在6號染色體上的基因數(shù)最多,有6個(gè),分布在4號染色體上的基因數(shù)最少,僅有1個(gè),分布在其余染色體上的基因數(shù)為2~4個(gè)。

    2.3抗性基因與已知抗性QTL比較

    將本研究檢測到的31個(gè)基因與前人研究確定的玉米灰斑病抗性QTL區(qū)間進(jìn)行比較,發(fā)現(xiàn)除了位于bin 1.11的2個(gè)基因(登錄號:GRMZM2G163912、GRMZM2G156543)和位于bin 2.04的1個(gè)基因(登錄號:GRMZM2G161932)、位于bin 5.02的1個(gè)基因(登錄號:GRMZM2G057129)外,其余27個(gè)基因均分布在前人研究的QTL區(qū)間內(nèi)(表3)。雖然以上4個(gè)基因不位于前人研究的QTL區(qū)間內(nèi),但都與已定位的QTL相距不遠(yuǎn),它們是否屬于新的QTL還有待進(jìn)一步研究確定。其中1號染色體上登錄號為GRMZM2G163912的基因與Zhang等[20]定位的QTL區(qū)間的右標(biāo)記umc1500在遺傳圖譜上只相距54.20 cM,與閆偉等[31]通過QTL元分析定位的QTL區(qū)間右標(biāo)記ias13a的距離更近,為47.12 cM,登錄號為GRMZM2G156543的基因與閆偉等[31]通過QTL元分析定位的QTL區(qū)間右標(biāo)記ias13a相距99.52 cM;2號染色體上登錄號為GRMZM2G161932的基因與Balint-Kurti等[18]定位的mmp167~umc2249 QTL區(qū)間的物理距離最近,為432 kb;5號染色體上登錄號為GRMZM2G057129的基因與Benson等[32]定位的QTL區(qū)間右邊界的物理距離只有14 kb。

    3討論

    3.1玉米灰斑病抗性鑒定與抗性種質(zhì)培育

    準(zhǔn)確的抗性鑒定是抗病育種及抗病基因定位的根本前提。人工培養(yǎng)玉米灰斑病的病原菌較難,發(fā)病要求的環(huán)境較難控制,而利用重病區(qū)的自然發(fā)病原理是玉米灰斑病抗性鑒定的有效途徑[35]。本試驗(yàn)地點(diǎn)位于云南省保山市市郊(海拔1 600 m,地理坐標(biāo)為25°08′N,99°10′E),屬于玉米灰斑病的重災(zāi)區(qū),當(dāng)?shù)貧夂驕嘏睗瘢镩g相對濕度較高,并且土壤屬微酸性紅壤,有利于病原菌的存活、侵染和在田間發(fā)病流行。目前,許多有關(guān)玉米灰斑病的研究都將試驗(yàn)地點(diǎn)設(shè)在云南省保山市市郊[20-21,23],而且已有研究發(fā)現(xiàn),該地區(qū)的玉米灰斑病病原菌為玉米尾孢菌[2,8-9]。由此可見,本研究所選試驗(yàn)地點(diǎn)可以從環(huán)境、氣候、病原菌等方面為玉米尾孢菌灰斑病的發(fā)生與流行提供條件,從而為本研究得到真實(shí)、可靠的表型數(shù)據(jù)提供重要保障。

    為了準(zhǔn)確評估研究材料的抗病情況,本研究在田間試驗(yàn)設(shè)計(jì)的基礎(chǔ)上采用嚴(yán)格的措施控制試驗(yàn)誤差,一是將關(guān)聯(lián)群體按照熱帶/亞熱帶、溫帶種質(zhì)分為2組,組內(nèi)采用隨機(jī)區(qū)組設(shè)計(jì),2次重復(fù);二是設(shè)置大量對照,每40行放1組對照(2個(gè)抗病對照材料,2個(gè)感病對照材料),以提高試驗(yàn)數(shù)據(jù)的準(zhǔn)確性。從試驗(yàn)結(jié)果看,對照材料抗病自交系Y32在2017年、2018年這2年的抗性鑒定結(jié)果中表現(xiàn)為一致抗病,感病自交系Ye478表現(xiàn)為一致感病;2017年232份玉米自交系的平均病級為4.79級,2018年225份玉米自交系的平均病級為4.96級。方差分析結(jié)果顯示,玉米灰斑病表型不僅在玉米材料間表現(xiàn)出極顯著差異,在年份與材料互作間也表現(xiàn)出極顯著差異。分析其原因,一方面由于玉米灰斑病抗性是由許多微效基因控制的數(shù)量性狀,易受環(huán)境影響[13-14];另一方面,也與玉米灰斑病抗性評價(jià)標(biāo)準(zhǔn)有關(guān)。在本研究中,評價(jià)玉米灰斑病抗性采用的是中國病害鑒定通用的葉斑病病級劃分標(biāo)準(zhǔn),即根據(jù)發(fā)病葉片部位及發(fā)病面積進(jìn)行分級,目前大多數(shù)玉米灰斑病研究都采用這個(gè)方法進(jìn)行抗性評價(jià)[20-21,23],該方法易于掌握和在田間實(shí)施,但是通過目測葉片發(fā)病面積進(jìn)行分級可能會存在一定的主觀性,要求調(diào)查者能準(zhǔn)確識別玉米灰斑病并正確評價(jià)其病級。本研究在調(diào)查過程中采用2人獨(dú)立評價(jià)的方法,再取2人評價(jià)值的平均值代表該材料的病級,以盡量避免因個(gè)人主觀因素帶來的表型誤差。

    盡管目前國內(nèi)外在玉米灰斑病抗性種質(zhì)鑒定方面進(jìn)行了大量研究,但是大多數(shù)都是針對玉蜀黍尾孢菌灰斑病抗性種質(zhì)鑒定,而且目前在世界范圍內(nèi)缺少抗玉米尾孢菌灰斑病的種質(zhì)資源。本研究通過2年的玉米尾孢菌灰斑病抗性鑒定,篩選出44份抗病材料和5份高抗材料,可以在玉米尾孢菌灰斑病的抗性遺傳改良中加以利用。對這些材料的來源進(jìn)行分析發(fā)現(xiàn),這些抗病材料大都來源于CIMMYT(82%)。進(jìn)一步查看241份玉米自交系中所有來源于CIMMYT的材料發(fā)現(xiàn),中抗及以上的材料占比達(dá)到了66.7%以上,而CIMMYT種質(zhì)基本上屬于熱帶/亞熱帶種質(zhì),表明來源于熱帶/亞熱帶的玉米種質(zhì)對于玉米尾孢菌灰斑病普遍具有較好的抗性。吳紀(jì)昌等[10]對來源廣泛的915份玉米自交系進(jìn)行灰斑病抗性鑒定,發(fā)現(xiàn)其中的高抗材料多數(shù)是從國外引進(jìn)的具有熱帶/亞熱帶血緣的材料。Shi等[36-37]通過對玉米種質(zhì)進(jìn)行抗性鑒定發(fā)現(xiàn),在中國主要的玉米種質(zhì)中,具有熱帶/亞熱帶血緣的PB種質(zhì)對灰斑病的抗性表現(xiàn)較好。綜上,今后應(yīng)加強(qiáng)熱帶/亞熱帶玉米種質(zhì)的研究和應(yīng)用,培育優(yōu)良的抗灰斑病玉米品種。

    3.2玉米灰斑病抗性位點(diǎn)與候選基因預(yù)測

    雖然本研究檢測到的與玉米灰斑病抗性顯著關(guān)聯(lián)的44個(gè)SNP位點(diǎn)中只有1個(gè)位點(diǎn)在2年內(nèi)均被檢測到,但與他人的研究結(jié)果進(jìn)行比較可以發(fā)現(xiàn),本研究檢測到的44個(gè)SNP位點(diǎn)中有34個(gè)SNP位點(diǎn)對應(yīng)31個(gè)基因,分布在除10號染色體外的其余染色體上。由此可以推測,本研究通過關(guān)聯(lián)分析得到的玉米灰斑病抗性位點(diǎn)是真實(shí)有效的,同時(shí)也進(jìn)一步證明,玉米灰斑病抗性是由許多微效基因控制的數(shù)量性狀[13-14]。本研究檢測到的44個(gè)與玉米灰斑病抗性顯著關(guān)聯(lián)的SNP位點(diǎn)的表型貢獻(xiàn)率都較低(4.72%~8.46%),其中2017年、2018年重復(fù)檢測到的SNP標(biāo)記AX86259524的表型貢獻(xiàn)率分別僅為5.15%、5.40%。正是因?yàn)檫@些位點(diǎn)的表型貢獻(xiàn)率較小、極易受外界環(huán)境的影響,因此它們被重復(fù)檢驗(yàn)的難度較大,說明表型鑒定對于受微效基因控制的數(shù)量性狀研究極為重要。

    將本研究檢測到的位于已知玉米灰斑病QTL區(qū)間內(nèi)的27個(gè)基因在生物信息數(shù)據(jù)庫中進(jìn)行比對,共找到8個(gè)與植物抗病相關(guān)的基因。其中登錄號為GRMZM2G163426基因(TRAM/LAG/CRN8 5基因,簡稱TLC5)、登錄號為GRMZM2G130278基因(內(nèi)腔蛋白保留受體家族蛋白編碼基因)和登錄號為GRMZM2G069177基因(自噬相關(guān)蛋白18編碼基因)都可以通過調(diào)控細(xì)胞程序性死亡來參與植物免疫反應(yīng)[38-39]。Pei等[40]研究發(fā)現(xiàn),自噬相關(guān)蛋白編碼基因在應(yīng)對小麥白粉病的免疫反應(yīng)中起著積極作用。登錄號為GRMZM5G889999基因(激酶-相互作用激酶1編碼基因)和登錄號為GRMZM2G374309基因(富半胱氨酸受體樣蛋白激酶10編碼基因)都包含絲氨酸/蘇氨酸激酶結(jié)構(gòu)域(STK),在植物免疫系統(tǒng)中充當(dāng)細(xì)胞表面受體并轉(zhuǎn)導(dǎo)信號,以觸發(fā)植物的先天性免疫反應(yīng)[41-42]。此外,許多研究都發(fā)現(xiàn),絲氨酸/蘇氨酸蛋白激酶在植物抗病方面有著重要作用[43-44]。已有的研究結(jié)果也多次將絲氨酸/蘇氨酸激酶基因作為玉米灰斑病抗性的候選基因[36, 45-46]。登錄號為GRMZM2G068665基因(吡咯啉-5-羧酸還原酶基因)參與植物脯氨酸合成代謝途徑,而游離脯氨酸可能在脅迫耐受和滲透調(diào)節(jié)中起作用[47-48]。郭紅蓮等[49]研究發(fā)現(xiàn),游離脯氨酸含量與玉米品種對灰斑病的抗性間呈正相關(guān)。登錄號為GRMZM2G007372基因(五肽重復(fù)序列蛋白編碼基因)和登錄號為GRMZM2G100146基因(組蛋白脫乙酰酶HDT2基因)都通過調(diào)控基因表達(dá)進(jìn)而參與到植物免疫防御過程中[50-51],在擬南芥、水稻中都有這2個(gè)基因與抗病性相關(guān)的報(bào)道[52-53]。雖然有報(bào)道表明,以上8個(gè)基因在植物中都有抗病防御的功能,但是它們是否對玉米灰斑病具有抗性還需進(jìn)行相關(guān)基因的功能驗(yàn)證。這8個(gè)候選基因不僅分別位于不同染色體及其區(qū)段上,其抗病機(jī)制也各不相同,并且其基因功能涉及信號轉(zhuǎn)導(dǎo)、轉(zhuǎn)錄調(diào)控、物質(zhì)代謝和結(jié)構(gòu)蛋白功能,暗示玉米灰斑病的抗病反應(yīng)是一個(gè)非常復(fù)雜的過程,可能涉及多種調(diào)控途徑的相互作用,因此在研究抗玉米灰斑病的過程中,還需要挖掘更多抗性基因和抗性機(jī)制。

    參考文獻(xiàn):

    [1]盧宗志,李艷君,李海春,等. 玉米灰斑病對玉米產(chǎn)量及產(chǎn)量特性的影響研究[J]. 玉米科學(xué), 2008, 16(6): 126-129.

    [2]趙立萍,王曉鳴,段燦星,等. 中國玉米灰斑病發(fā)生現(xiàn)狀與未來擴(kuò)散趨勢分析[J]. 中國農(nóng)業(yè)科學(xué), 2015, 48(18): 3612-3626.

    [3]HILTY J W, HADDEN C H, GARDEN F T. Response of maize hybrids and inbred lines to gray leaf spot disease and the effects on yield in Tennessee[J]. Plant Disease Reporter, 1979, 63(6): 515-518.

    [4]KINYUA Z M, SMITH J J, KIBATA G N, et al. Status of grey leaf spot disease in Kenyan maize production ecosystems[J]. African Crop Science Journal, 2010, 18(4): 183-194.

    [5]李曉,張小飛,崔麗娜,等. 警惕四川玉米灰斑病的發(fā)生危害[J]. 四川農(nóng)業(yè)科技, 2011(12): 37.

    [6]周惠萍,吳景芝,李月秋,等. 云南省玉米灰斑病發(fā)生規(guī)律研究[J]. 西南農(nóng)業(yè)學(xué)報(bào), 2011, 24(6): 2207-2212.

    [7]張益先,呂國忠,梁景頤,等. 玉米灰斑病菌生物學(xué)特性研究[J]. 植物病理學(xué)報(bào), 2003, 33(4): 292-295.

    [8]劉慶奎,秦子惠,張小利,等. 中國玉米灰斑病病原菌的鑒定及其基本特征研究[J]. 中國農(nóng)業(yè)科學(xué), 2013, 46(19): 4044-4057.

    [9]LIU K J, XU X D. First report of gray leaf spot of maize caused by Cercospora zeina in China[J]. Plant Disease, 2013, 97(12): 1656.

    [10]吳紀(jì)昌,馬麗君,王作英. 玉米抗尾孢菌葉斑病鑒定與抗病材料利用[J]. 遼寧農(nóng)業(yè)科學(xué), 1997(5): 25-28.

    [11]董懷玉,姜鈺,王麗娟,等. 玉米種質(zhì)資源抗灰斑病鑒定與評價(jià)[J]. 植物遺傳資源學(xué)報(bào), 2005, 6(4): 441-443.

    [12]呂香玲,李新海,陳陽,等. 玉米種質(zhì)抗灰斑病鑒定與評價(jià)[J]. 玉米科學(xué), 2011, 19(6): 125-128.

    [13]李世強(qiáng),陳威,譚靜,等. 玉米灰斑病的抗性機(jī)理研究進(jìn)展[J]. 生物技術(shù)進(jìn)展, 2011, 1(2): 112-115.

    [14]ABEBE M, MARIA A. Genetic analysis of resistance to gray leaf spot of maize inbred lines[J]. Crop Science, 2005, 45(1): 163-170.

    [15]BUBECK D M, GOODMAN M M, BEAVIS W D, et al. Quantitative trait loci controlling resistance to gray leaf spot in maize[J]. Crop Science, 1993, 33: 838-847.

    [16]CLEMENTS M J, DUDLEY J W, WHITE D G. Quantitative trait loci associated with resistance to gray leaf spot of corn[J]. Phytopathology, 2000, 90(9):1018-1025.

    [17]LEHMENSIEK A, ESTERHUIZEN A M, VAN-STADEN D, et al. Genetic mapping of gray leaf spot (GLS) resistance genes in maize [J]. Theoretical and Applied Genetics, 2001, 103: 797-803.

    [18]BALINT-KURTI P J, WISSER R, ZWONITZER J C. Use of the IBM population for precise mapping of quantitative trait loci for gray leaf spot resistance in maize[J]. Crop Science, 2008, 48(5): 1696-1704.

    [19]POZAR G, BUTRUILLE D, SILVA H D, et al. Mapping and validation of quantitative trait loci for resistance to Cercospora zea-maydis infection in tropical maize (Zea mays L.)[J]. Theoretical and Applied Genetics, 2009, 118(3): 533-564.

    [20]ZHANG Y, XU L, FAN X M, et al. QTL mapping of resistance to gray leaf spot in maize[J]. Theoretical and Applied Genetics, 2012, 125(8): 1797-1808.

    [21]XU L, ZHANG Y, SHAO S Q, et al. High-resolution mapping and characterization of qRgls2, a major quantitative trait locus involved in maize resistance to gray leaf spot[J]. BMC Plant Biology, 2014, 14: 230-240.

    [22]BERGER D K, CARSTENS M, KORSMAN J N, et al. Mapping QTL conferring resistance in maize to gray leaf spot disease caused by Cercospora zeina[J]. BMC Genetics, 2014, 15: 60.

    [23]LIU L, TAN J, ZHANG Y D, et al. QTL mapping for gray leaf spot resistance in a tropical maize population[J]. Plant Disease, 2016, 100(2): 304-312.

    [24]HE W Z, YANG L, LENG Y F, et al. QTL mapping for resistance of maize to grey leaf spot[J]. Journal of Phytopathology, 2018, 166(3): 167-176.

    [25]WANG W Y, BARRATT B J, CLAYTON D G, et al. Genome-wide association studies: theoretical and practical concerns[J]. Nature Review Genetics, 2005, 6: 109-118.

    [26]WITTE J S. Genome-wide association studies and beyond[J]. Annual Review of Public Health, 2010, 31(1): 9-20.

    [27]ZHU C S, GORE M, BUCKLER E S, et al. Status and prospects of association mapping in plants[J]. The Plant Genome, 2008, 1(1): 5-20.

    [28]MAMMADOV J A, SUN X C, GAO Y X, et al. Combining powers of linkage and association mapping for precise dissection of QTL controlling resistance to gray leaf spot disease in maize (Zea mays L.)[J]. BMC Genomics, 2015, 16(1): 916-931.

    [29]KUKI M C, SCAPIM C A, ROSSI E S, et al. Genome wide association study for gray leaf spot resistance in tropical maize core[J]. PLoS One, 2018, 13(6): e0199539.

    [30]王曉鳴,石潔,晉齊鳴,等. 玉米病蟲害田間手冊——病蟲害鑒別與抗性鑒定[M]. 北京: 中國農(nóng)業(yè)科學(xué)技術(shù)出版社, 2010: 8.

    [31]閆偉,李元,宋茂興,等. 玉米抗灰斑病QTL元分析及其驗(yàn)證[J]. 作物學(xué)報(bào), 2016, 42(5): 758-767.

    [32]BENSON J M, POLAND J A, BENSON B M, et al. Resistance to gray leaf spot of maize: genetic architecture and mechanisms elucidated through nested association mapping and near-isogenic line analysis[J]. PLoS Genetics, 2015, 11(3): e1005045.

    [33]ZWONITZER J C, COLES N D, KRAKOWSKY M D, et al. Mapping resistance quantitative trait loci for three foliar diseases in a maize recombinant inbred line population-evidence for multiple disease resistance?[J]. Phytopathology, 2010, 100(1): 72-79.

    [34]JINES M P. Identification of quantitative trait loci (QTL) for gray leaf spot resistance, maturity, and grain yield in a semi-tropical recombinant inbred population of maize[D]. North Carolina, USA: North Carolina State University, 2004:55-69.

    [35]譚靜,羅吉,孫彩梅,等. 玉米灰斑病抗性種質(zhì)與基因位點(diǎn)的研究進(jìn)展[J]. 種子, 2019, 38(8): 57-60,64.

    [36]SHI L Y, LV X L, WENG J F, et al. Genetic characterization and linkage disequilibrium mapping of resistance to gray leaf spot in maize (Zea mays L.)[J].The Crop Journal, 2014(Z1): 132-143.

    [37]劉可杰,董懷玉,姜鈺,等. 300份玉米種質(zhì)對兩種灰斑病菌的抗性評價(jià)[J]. 玉米科學(xué), 2018, 26(4): 162-165.

    [38]TAKAHASHI Y, BERBERICH T, KANZAKI H, et al. Serine palmitoyltransferase, the first step enzyme in sphingolipid biosynthesis, is involved in non-host resistance[J]. Molecular Plant-Microbe Interactions, 2009, 22(1): 31-38.

    [39]XU G Y, LI S Z, XIE K, et al. Plant ERD2‐like proteins function as endoplasmic reticulum luminal protein receptors and participate in programmed cell death during innate immunity[J]. The Plant Journal, 2012, 72(1): 57-69.

    猜你喜歡
    關(guān)聯(lián)分析玉米
    收玉米啦!
    我的玉米送給你
    為什么玉米長“胡須”
    玉米骨干親本及其衍生系中基因的序列變異及與株高等性狀的關(guān)聯(lián)分析
    玉米骨干親本及其衍生系中基因的序列變異及與株高等性狀的關(guān)聯(lián)分析
    基于隨機(jī)函數(shù)Petri網(wǎng)的系統(tǒng)動(dòng)力學(xué)關(guān)聯(lián)分析模型
    關(guān)聯(lián)分析技術(shù)在學(xué)生成績分析中的應(yīng)用
    基于關(guān)聯(lián)分析的學(xué)生活動(dòng)參與度與高校社團(tuán)管理實(shí)證研究
    不同的數(shù)據(jù)挖掘方法分類對比研究
    老司机在亚洲福利影院| 母亲3免费完整高清在线观看| 亚洲18禁久久av| 熟女人妻精品中文字幕| 午夜a级毛片| 精品电影一区二区在线| 中文字幕人成人乱码亚洲影| 亚洲av免费在线观看| 日韩欧美 国产精品| 伊人久久大香线蕉亚洲五| 久久精品91蜜桃| 国产高清视频在线播放一区| 午夜影院日韩av| 黄色丝袜av网址大全| 天美传媒精品一区二区| 九色成人免费人妻av| 欧美三级亚洲精品| 亚洲美女黄片视频| 国产成+人综合+亚洲专区| h日本视频在线播放| 国产精品久久久久久人妻精品电影| 老司机福利观看| 精品一区二区三区视频在线 | 欧美日韩中文字幕国产精品一区二区三区| 午夜视频国产福利| 99久国产av精品| 国产一级毛片七仙女欲春2| 亚洲不卡免费看| 丰满人妻一区二区三区视频av | 国产淫片久久久久久久久 | 夜夜夜夜夜久久久久| 三级毛片av免费| a在线观看视频网站| 精品乱码久久久久久99久播| 露出奶头的视频| 一二三四社区在线视频社区8| 国产真人三级小视频在线观看| av国产免费在线观看| 国产高清激情床上av| 国产黄片美女视频| 最近最新中文字幕大全电影3| 很黄的视频免费| 美女cb高潮喷水在线观看| 欧美中文日本在线观看视频| 午夜视频国产福利| 怎么达到女性高潮| 免费搜索国产男女视频| 两个人视频免费观看高清| 免费大片18禁| 男女视频在线观看网站免费| 亚洲专区中文字幕在线| 啪啪无遮挡十八禁网站| 亚洲国产精品成人综合色| 欧美黑人欧美精品刺激| 特大巨黑吊av在线直播| 欧美黄色片欧美黄色片| 人妻夜夜爽99麻豆av| 中文字幕av在线有码专区| 又黄又粗又硬又大视频| 久久久国产成人免费| 国产日本99.免费观看| 久久久久精品国产欧美久久久| 亚洲成人久久爱视频| 男人和女人高潮做爰伦理| 国产黄色小视频在线观看| 欧美性猛交黑人性爽| 老鸭窝网址在线观看| av女优亚洲男人天堂| 欧美日韩一级在线毛片| 欧洲精品卡2卡3卡4卡5卡区| 国产色婷婷99| 欧美国产日韩亚洲一区| 国内精品一区二区在线观看| 99久久九九国产精品国产免费| 亚洲精品在线美女| 精品久久久久久,| 十八禁网站免费在线| 久久久国产精品麻豆| 最近视频中文字幕2019在线8| 美女高潮喷水抽搐中文字幕| 欧美最新免费一区二区三区 | 久久午夜亚洲精品久久| 国产v大片淫在线免费观看| 国产三级黄色录像| 中文字幕av成人在线电影| 老司机午夜十八禁免费视频| 国产 一区 欧美 日韩| 色老头精品视频在线观看| 中文字幕精品亚洲无线码一区| 免费在线观看影片大全网站| 午夜激情福利司机影院| 亚洲熟妇熟女久久| av片东京热男人的天堂| 久久精品国产自在天天线| 99热这里只有精品一区| 日韩有码中文字幕| 中文字幕熟女人妻在线| 国语自产精品视频在线第100页| 欧美精品啪啪一区二区三区| 一个人免费在线观看的高清视频| 香蕉av资源在线| 亚洲片人在线观看| 男人舔女人下体高潮全视频| 在线天堂最新版资源| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 一个人观看的视频www高清免费观看| av国产免费在线观看| netflix在线观看网站| 中文在线观看免费www的网站| 99久久99久久久精品蜜桃| 99国产精品一区二区三区| 日韩中文字幕欧美一区二区| 一区福利在线观看| 69人妻影院| 狂野欧美白嫩少妇大欣赏| 男人的好看免费观看在线视频| 亚洲成a人片在线一区二区| tocl精华| 国产成年人精品一区二区| 国产毛片a区久久久久| 国产精品98久久久久久宅男小说| 成年版毛片免费区| 男人舔女人下体高潮全视频| 我要搜黄色片| 欧美在线黄色| 中文资源天堂在线| 国产爱豆传媒在线观看| 老汉色av国产亚洲站长工具| 久久人人精品亚洲av| 男女午夜视频在线观看| 日韩人妻高清精品专区| 国产精品日韩av在线免费观看| 欧美成人一区二区免费高清观看| 十八禁网站免费在线| 久久久久久久午夜电影| 一级毛片女人18水好多| 美女被艹到高潮喷水动态| 久久天躁狠狠躁夜夜2o2o| 变态另类丝袜制服| 亚洲国产精品成人综合色| 色老头精品视频在线观看| 欧美成人a在线观看| 嫁个100分男人电影在线观看| 99久久综合精品五月天人人| 国产精品亚洲美女久久久| 欧美一区二区精品小视频在线| 精品人妻一区二区三区麻豆 | 日韩国内少妇激情av| 亚洲成人久久爱视频| 桃色一区二区三区在线观看| 免费一级毛片在线播放高清视频| 成人永久免费在线观看视频| 国产真实伦视频高清在线观看 | 日韩中文字幕欧美一区二区| 亚洲成人免费电影在线观看| 欧美性感艳星| 亚洲七黄色美女视频| 岛国在线观看网站| 精品欧美国产一区二区三| 搡女人真爽免费视频火全软件 | 伊人久久精品亚洲午夜| 噜噜噜噜噜久久久久久91| 99久久精品国产亚洲精品| 精品久久久久久成人av| 日本精品一区二区三区蜜桃| 国产私拍福利视频在线观看| 在线观看日韩欧美| 99视频精品全部免费 在线| 久久国产乱子伦精品免费另类| 成人午夜高清在线视频| 国产亚洲欧美在线一区二区| 国产熟女xx| 此物有八面人人有两片| 身体一侧抽搐| 久久香蕉国产精品| 亚洲国产高清在线一区二区三| 国产久久久一区二区三区| 麻豆成人av在线观看| 久久国产乱子伦精品免费另类| 在线观看66精品国产| 国产国拍精品亚洲av在线观看 | 香蕉丝袜av| 熟女电影av网| 国产真实乱freesex| 国模一区二区三区四区视频| 亚洲片人在线观看| 精品久久久久久久久久久久久| 免费大片18禁| 少妇人妻一区二区三区视频| 床上黄色一级片| 网址你懂的国产日韩在线| 精品久久久久久久久久免费视频| 欧美日本视频| 日本一本二区三区精品| 老司机午夜福利在线观看视频| 亚洲狠狠婷婷综合久久图片| 成年免费大片在线观看| 夜夜躁狠狠躁天天躁| svipshipincom国产片| 国模一区二区三区四区视频| 噜噜噜噜噜久久久久久91| 人人妻人人澡欧美一区二区| 黄色视频,在线免费观看| 给我免费播放毛片高清在线观看| 一个人观看的视频www高清免费观看| 18禁裸乳无遮挡免费网站照片| 欧美+日韩+精品| 国产伦人伦偷精品视频| 亚洲国产欧洲综合997久久,| 国产精品香港三级国产av潘金莲| 久久久久久久久中文| 欧美高清成人免费视频www| 日韩欧美在线乱码| 国产熟女xx| 午夜免费男女啪啪视频观看 | 少妇的逼好多水| 国产欧美日韩精品亚洲av| a级毛片a级免费在线| 欧美性猛交╳xxx乱大交人| 乱人视频在线观看| 最近最新中文字幕大全免费视频| 午夜免费激情av| 国产精华一区二区三区| 日韩亚洲欧美综合| 两人在一起打扑克的视频| 成年免费大片在线观看| 久久久久久久亚洲中文字幕 | 深爱激情五月婷婷| 精品人妻1区二区| 国产高清激情床上av| 久久精品国产亚洲av香蕉五月| 此物有八面人人有两片| 亚洲精品一卡2卡三卡4卡5卡| 又粗又爽又猛毛片免费看| 国产真实伦视频高清在线观看 | 国产一区二区亚洲精品在线观看| 少妇裸体淫交视频免费看高清| 舔av片在线| 男人的好看免费观看在线视频| 三级毛片av免费| 成熟少妇高潮喷水视频| 99久久精品一区二区三区| 国内少妇人妻偷人精品xxx网站| 88av欧美| 久9热在线精品视频| 免费在线观看成人毛片| 国产精品99久久久久久久久| 精品电影一区二区在线| 国产精品日韩av在线免费观看| 国内精品久久久久久久电影| 看片在线看免费视频| 好看av亚洲va欧美ⅴa在| 欧美+亚洲+日韩+国产| 此物有八面人人有两片| av国产免费在线观看| 欧美bdsm另类| 岛国在线观看网站| 欧美最黄视频在线播放免费| 露出奶头的视频| 一二三四社区在线视频社区8| 无人区码免费观看不卡| 女人高潮潮喷娇喘18禁视频| 亚洲成人久久性| 婷婷精品国产亚洲av在线| 色综合站精品国产| 性欧美人与动物交配| xxxwww97欧美| 国产中年淑女户外野战色| av黄色大香蕉| 亚洲精华国产精华精| 国产亚洲av嫩草精品影院| 欧美国产日韩亚洲一区| 亚洲avbb在线观看| 欧美午夜高清在线| 人人妻人人澡欧美一区二区| 婷婷精品国产亚洲av| 搞女人的毛片| 国产精华一区二区三区| 中文字幕精品亚洲无线码一区| 欧美3d第一页| 日韩大尺度精品在线看网址| 99热这里只有精品一区| 中文字幕人妻熟人妻熟丝袜美 | 又爽又黄无遮挡网站| 看免费av毛片| 国语自产精品视频在线第100页| 好男人电影高清在线观看| 国产一区二区在线观看日韩 | 听说在线观看完整版免费高清| 久久天躁狠狠躁夜夜2o2o| 男插女下体视频免费在线播放| 免费av不卡在线播放| 色哟哟哟哟哟哟| 欧美性猛交黑人性爽| 国产一区二区亚洲精品在线观看| 国产精品嫩草影院av在线观看 | 成人精品一区二区免费| 中文字幕av在线有码专区| 国产精品永久免费网站| 两个人看的免费小视频| 成人欧美大片| 在线十欧美十亚洲十日本专区| 久久精品亚洲精品国产色婷小说| 中文亚洲av片在线观看爽| 午夜福利在线在线| 深夜精品福利| 亚洲,欧美精品.| 国产高清视频在线观看网站| 亚洲av电影在线进入| 一夜夜www| 少妇的逼好多水| 国产激情偷乱视频一区二区| 亚洲成av人片免费观看| 欧美色欧美亚洲另类二区| 亚洲人成网站高清观看| 床上黄色一级片| 制服丝袜大香蕉在线| 成人性生交大片免费视频hd| 在线十欧美十亚洲十日本专区| 成人国产一区最新在线观看| 国产欧美日韩精品一区二区| 日本三级黄在线观看| 两个人看的免费小视频| 亚洲国产精品成人综合色| 性色avwww在线观看| 99国产精品一区二区蜜桃av| 久久草成人影院| 美女高潮的动态| 久久人妻av系列| 色综合婷婷激情| 成人av在线播放网站| 国产高清videossex| 国产日本99.免费观看| 亚洲精品乱码久久久v下载方式 | 波多野结衣高清作品| or卡值多少钱| 国产精品自产拍在线观看55亚洲| 亚洲乱码一区二区免费版| 久久精品国产亚洲av香蕉五月| 淫秽高清视频在线观看| 免费观看的影片在线观看| 欧美日韩福利视频一区二区| av黄色大香蕉| 国产主播在线观看一区二区| 一进一出抽搐动态| 亚洲专区国产一区二区| 国产一区二区亚洲精品在线观看| 欧美极品一区二区三区四区| 每晚都被弄得嗷嗷叫到高潮| av在线天堂中文字幕| 欧美大码av| 日韩精品青青久久久久久| 69av精品久久久久久| 啦啦啦韩国在线观看视频| 91av网一区二区| 精品欧美国产一区二区三| 国产 一区 欧美 日韩| 男女午夜视频在线观看| 黄片大片在线免费观看| 成人永久免费在线观看视频| 九九久久精品国产亚洲av麻豆| 欧美成狂野欧美在线观看| 国产乱人伦免费视频| 亚洲精品粉嫩美女一区| 免费看日本二区| 在线观看免费视频日本深夜| 久久亚洲精品不卡| 亚洲无线在线观看| 日本在线视频免费播放| 神马国产精品三级电影在线观看| 叶爱在线成人免费视频播放| 女警被强在线播放| 真人做人爱边吃奶动态| 日本 av在线| 麻豆成人午夜福利视频| 日韩国内少妇激情av| 老汉色av国产亚洲站长工具| 亚洲avbb在线观看| 午夜影院日韩av| 国产亚洲av嫩草精品影院| 18+在线观看网站| 午夜福利在线在线| 久9热在线精品视频| 老熟妇乱子伦视频在线观看| 嫩草影院精品99| 国产在线精品亚洲第一网站| 亚洲自拍偷在线| 日本黄色片子视频| 久久精品91蜜桃| av在线天堂中文字幕| 国产aⅴ精品一区二区三区波| 亚洲精品亚洲一区二区| 毛片女人毛片| 亚洲色图av天堂| 少妇的逼好多水| 亚洲第一电影网av| 国产野战对白在线观看| 麻豆国产av国片精品| 欧美日韩国产亚洲二区| 国产高清视频在线播放一区| 中国美女看黄片| 免费av观看视频| 亚洲18禁久久av| 欧美中文综合在线视频| av中文乱码字幕在线| www.www免费av| 亚洲中文日韩欧美视频| 乱人视频在线观看| 久久久久久久亚洲中文字幕 | 国产伦在线观看视频一区| 亚洲av免费在线观看| 美女高潮的动态| 午夜福利高清视频| 精品久久久久久久末码| 国产成+人综合+亚洲专区| 一进一出好大好爽视频| 黄色片一级片一级黄色片| 岛国在线观看网站| АⅤ资源中文在线天堂| 99热这里只有精品一区| 人人妻,人人澡人人爽秒播| a在线观看视频网站| 伊人久久精品亚洲午夜| 成人特级av手机在线观看| 最新美女视频免费是黄的| 成人特级黄色片久久久久久久| 51午夜福利影视在线观看| 高潮久久久久久久久久久不卡| 五月伊人婷婷丁香| 亚洲欧美激情综合另类| 国产免费一级a男人的天堂| 在线观看美女被高潮喷水网站 | 国产高清视频在线观看网站| 欧美中文综合在线视频| 午夜影院日韩av| 国产黄色小视频在线观看| 一级毛片女人18水好多| 亚洲片人在线观看| 99精品久久久久人妻精品| 女同久久另类99精品国产91| 一个人免费在线观看电影| 十八禁人妻一区二区| 欧美一级a爱片免费观看看| 伊人久久大香线蕉亚洲五| 久久精品91无色码中文字幕| 欧美激情久久久久久爽电影| 久久精品影院6| 国产日本99.免费观看| 久久久精品大字幕| 日韩中文字幕欧美一区二区| 中出人妻视频一区二区| 久久精品91蜜桃| 欧美一区二区亚洲| 国内精品久久久久精免费| 国产一区在线观看成人免费| 别揉我奶头~嗯~啊~动态视频| 婷婷精品国产亚洲av| 叶爱在线成人免费视频播放| 淫妇啪啪啪对白视频| 亚洲精品在线观看二区| 亚洲国产色片| 国产精品99久久久久久久久| 亚洲无线在线观看| 人妻丰满熟妇av一区二区三区| 国产蜜桃级精品一区二区三区| 欧美成人免费av一区二区三区| 日本与韩国留学比较| 别揉我奶头~嗯~啊~动态视频| 在线天堂最新版资源| 天天躁日日操中文字幕| 国产精品乱码一区二三区的特点| 乱人视频在线观看| 免费一级毛片在线播放高清视频| 亚洲国产欧洲综合997久久,| 色尼玛亚洲综合影院| 美女高潮喷水抽搐中文字幕| 在线国产一区二区在线| 99国产精品一区二区蜜桃av| 熟女电影av网| 中文字幕精品亚洲无线码一区| 97人妻精品一区二区三区麻豆| 亚洲国产中文字幕在线视频| 88av欧美| 欧美日本视频| 99国产精品一区二区蜜桃av| 亚洲精品一卡2卡三卡4卡5卡| 精品欧美国产一区二区三| 国产精品影院久久| 人妻久久中文字幕网| 丰满的人妻完整版| 国产激情偷乱视频一区二区| 欧美zozozo另类| 精品国产超薄肉色丝袜足j| 欧美日韩一级在线毛片| 美女免费视频网站| 黑人欧美特级aaaaaa片| 丰满乱子伦码专区| 757午夜福利合集在线观看| 嫁个100分男人电影在线观看| 亚洲 国产 在线| 精品午夜福利视频在线观看一区| 欧美黄色片欧美黄色片| 国产熟女xx| 久久天躁狠狠躁夜夜2o2o| 国产高潮美女av| 久久精品国产99精品国产亚洲性色| 国产免费av片在线观看野外av| 99精品在免费线老司机午夜| 黄色视频,在线免费观看| h日本视频在线播放| 久9热在线精品视频| 无遮挡黄片免费观看| 久久伊人香网站| 麻豆国产97在线/欧美| 国产成+人综合+亚洲专区| 亚洲精品久久国产高清桃花| 久久精品人妻少妇| 内地一区二区视频在线| 亚洲精品国产精品久久久不卡| 成人高潮视频无遮挡免费网站| 美女cb高潮喷水在线观看| 亚洲精品一区av在线观看| 精品人妻一区二区三区麻豆 | 精品欧美国产一区二区三| 日韩欧美精品v在线| www.熟女人妻精品国产| 一个人免费在线观看的高清视频| av福利片在线观看| 18禁裸乳无遮挡免费网站照片| 在线a可以看的网站| 亚洲国产精品成人综合色| 熟妇人妻久久中文字幕3abv| 90打野战视频偷拍视频| 看黄色毛片网站| 亚洲狠狠婷婷综合久久图片| 在线国产一区二区在线| 国产一区二区在线av高清观看| 精品久久久久久久久久免费视频| 欧美黑人欧美精品刺激| 夜夜躁狠狠躁天天躁| 久久精品影院6| 国产精品野战在线观看| 天堂√8在线中文| 看免费av毛片| av欧美777| 国产久久久一区二区三区| 久久久久久久久大av| 久久国产精品人妻蜜桃| 男女床上黄色一级片免费看| 免费在线观看日本一区| 国产精品日韩av在线免费观看| 18+在线观看网站| 日韩欧美免费精品| 国产黄片美女视频| 国产中年淑女户外野战色| 国产精品久久久久久精品电影| 窝窝影院91人妻| 老鸭窝网址在线观看| 日韩精品青青久久久久久| 99久久成人亚洲精品观看| 国产精品爽爽va在线观看网站| 在线观看免费视频日本深夜| 日本黄色片子视频| 中文字幕人妻丝袜一区二区| 国内精品一区二区在线观看| 好男人在线观看高清免费视频| 中文字幕人成人乱码亚洲影| 欧美最黄视频在线播放免费| 欧美黑人巨大hd| 亚洲精品在线观看二区| xxx96com| av视频在线观看入口| 亚洲专区中文字幕在线| 级片在线观看| 国内毛片毛片毛片毛片毛片| 色哟哟哟哟哟哟| 国产精品久久电影中文字幕| 淫秽高清视频在线观看| 全区人妻精品视频| 少妇的丰满在线观看| 国产精品一区二区三区四区久久| 国产精品乱码一区二三区的特点| 日韩有码中文字幕| 亚洲激情在线av| 午夜免费男女啪啪视频观看 | 九九在线视频观看精品| x7x7x7水蜜桃| 亚洲,欧美精品.| 欧美精品啪啪一区二区三区| 18禁黄网站禁片午夜丰满| 久久精品国产清高在天天线| 内射极品少妇av片p| 精品欧美国产一区二区三| 3wmmmm亚洲av在线观看| 在线看三级毛片| 村上凉子中文字幕在线| 国产97色在线日韩免费| 少妇人妻一区二区三区视频| 99热这里只有精品一区| 日韩欧美一区二区三区在线观看| 在线播放国产精品三级| 长腿黑丝高跟| 中文字幕av在线有码专区| 中文字幕熟女人妻在线| 国内精品一区二区在线观看| 亚洲精品久久国产高清桃花| 国产一区二区三区视频了| 99精品久久久久人妻精品| 亚洲精品在线观看二区| 蜜桃亚洲精品一区二区三区| 久久久成人免费电影| 窝窝影院91人妻| 欧美日韩乱码在线| 99视频精品全部免费 在线| 亚洲国产色片|