• <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ù)挖掘方法分類對比研究
    国产精品98久久久久久宅男小说| 亚洲国产欧洲综合997久久,| 波多野结衣高清无吗| 国产精品国产高清国产av| 丁香六月欧美| 久久久精品大字幕| 日韩欧美在线二视频| 两人在一起打扑克的视频| 成人性生交大片免费视频hd| 男插女下体视频免费在线播放| 观看美女的网站| 欧美日韩中文字幕国产精品一区二区三区| 最新中文字幕久久久久 | 亚洲,欧美精品.| 久久久成人免费电影| 丁香欧美五月| 此物有八面人人有两片| 夜夜看夜夜爽夜夜摸| 国产亚洲av高清不卡| 757午夜福利合集在线观看| 香蕉久久夜色| 波多野结衣巨乳人妻| 国产精品电影一区二区三区| 日韩欧美国产在线观看| 黄色丝袜av网址大全| www.精华液| 国产单亲对白刺激| 三级毛片av免费| 精品一区二区三区视频在线 | 18禁美女被吸乳视频| 精品久久久久久久毛片微露脸| 色吧在线观看| 一进一出好大好爽视频| 最新中文字幕久久久久 | 偷拍熟女少妇极品色| 久久久久久久午夜电影| 亚洲熟妇熟女久久| 日韩欧美三级三区| 国产私拍福利视频在线观看| 精品国产三级普通话版| 中文字幕人成人乱码亚洲影| avwww免费| 国产久久久一区二区三区| 成人鲁丝片一二三区免费| 中文字幕精品亚洲无线码一区| 色精品久久人妻99蜜桃| 久久久久国内视频| 国产精品久久视频播放| 亚洲黑人精品在线| 一个人观看的视频www高清免费观看 | 国产成+人综合+亚洲专区| 麻豆一二三区av精品| 免费一级毛片在线播放高清视频| 脱女人内裤的视频| 午夜精品一区二区三区免费看| 亚洲欧洲精品一区二区精品久久久| 亚洲精品乱码久久久v下载方式 | 桃色一区二区三区在线观看| 国产日本99.免费观看| 美女 人体艺术 gogo| av片东京热男人的天堂| 欧美黄色淫秽网站| 成年免费大片在线观看| 国产视频内射| 色av中文字幕| svipshipincom国产片| 欧洲精品卡2卡3卡4卡5卡区| 久久久久精品国产欧美久久久| 中文在线观看免费www的网站| 国产高清videossex| 九九久久精品国产亚洲av麻豆 | 国产三级黄色录像| 欧美性猛交黑人性爽| 久久久久久大精品| 亚洲av电影在线进入| 久久精品综合一区二区三区| 亚洲狠狠婷婷综合久久图片| 久久中文字幕一级| 国产伦人伦偷精品视频| а√天堂www在线а√下载| 偷拍熟女少妇极品色| 男女之事视频高清在线观看| 搡老妇女老女人老熟妇| 亚洲国产精品成人综合色| 亚洲五月天丁香| 丝袜人妻中文字幕| 禁无遮挡网站| 1024手机看黄色片| 国产美女午夜福利| e午夜精品久久久久久久| 日本撒尿小便嘘嘘汇集6| 欧美3d第一页| 精品99又大又爽又粗少妇毛片 | 国产久久久一区二区三区| 国产av一区在线观看免费| 中文字幕最新亚洲高清| 91老司机精品| 亚洲中文日韩欧美视频| 亚洲欧美日韩东京热| 手机成人av网站| 亚洲专区国产一区二区| 999久久久精品免费观看国产| 午夜福利高清视频| 嫩草影视91久久| 国产真人三级小视频在线观看| 日本与韩国留学比较| 亚洲乱码一区二区免费版| 97人妻精品一区二区三区麻豆| 少妇丰满av| 黄色女人牲交| 麻豆国产97在线/欧美| 欧美性猛交╳xxx乱大交人| 亚洲av第一区精品v没综合| 亚洲精华国产精华精| 亚洲成av人片免费观看| 国内精品久久久久久久电影| 校园春色视频在线观看| 亚洲第一电影网av| 99热精品在线国产| 亚洲午夜理论影院| 亚洲欧美日韩东京热| 又紧又爽又黄一区二区| 身体一侧抽搐| 亚洲av成人不卡在线观看播放网| 色在线成人网| 国产精品av视频在线免费观看| 国产精品免费一区二区三区在线| 日韩三级视频一区二区三区| 久久久精品欧美日韩精品| 日韩欧美免费精品| 淫秽高清视频在线观看| 国产伦人伦偷精品视频| 狠狠狠狠99中文字幕| 麻豆av在线久日| 黄片大片在线免费观看| 青草久久国产| 麻豆国产av国片精品| 18禁黄网站禁片午夜丰满| 国产精品美女特级片免费视频播放器 | 亚洲欧美一区二区三区黑人| 狠狠狠狠99中文字幕| 亚洲五月天丁香| 18美女黄网站色大片免费观看| 欧美黑人巨大hd| 母亲3免费完整高清在线观看| a级毛片在线看网站| 国产成人影院久久av| 成人av在线播放网站| 三级毛片av免费| 在线免费观看不下载黄p国产 | 亚洲成人中文字幕在线播放| 黄色女人牲交| 免费高清视频大片| 久久久久久久午夜电影| 色综合站精品国产| 色综合站精品国产| 黄色视频,在线免费观看| 黄色视频,在线免费观看| 亚洲精品中文字幕一二三四区| cao死你这个sao货| 亚洲五月天丁香| 操出白浆在线播放| 两人在一起打扑克的视频| 窝窝影院91人妻| 久久精品亚洲精品国产色婷小说| 免费av不卡在线播放| 久久国产精品影院| 嫁个100分男人电影在线观看| 欧美日韩一级在线毛片| av片东京热男人的天堂| 91麻豆av在线| 亚洲第一电影网av| 精品国产乱子伦一区二区三区| 精品国内亚洲2022精品成人| 国产激情偷乱视频一区二区| www日本在线高清视频| 黑人巨大精品欧美一区二区mp4| 午夜福利高清视频| av欧美777| 噜噜噜噜噜久久久久久91| 麻豆av在线久日| 99热6这里只有精品| 国产久久久一区二区三区| 99热精品在线国产| h日本视频在线播放| 此物有八面人人有两片| 久久久成人免费电影| 国产精品女同一区二区软件 | 久久亚洲真实| 可以在线观看的亚洲视频| 校园春色视频在线观看| 窝窝影院91人妻| 国产乱人伦免费视频| 最好的美女福利视频网| 国产精品,欧美在线| 嫩草影院精品99| 一本久久中文字幕| 最新在线观看一区二区三区| 国产av一区在线观看免费| 好看av亚洲va欧美ⅴa在| 久久香蕉精品热| 亚洲成人精品中文字幕电影| 99热只有精品国产| av在线天堂中文字幕| tocl精华| 一卡2卡三卡四卡精品乱码亚洲| 精品国内亚洲2022精品成人| 国产精品一区二区精品视频观看| 国产高清视频在线播放一区| 一a级毛片在线观看| or卡值多少钱| 午夜福利视频1000在线观看| 首页视频小说图片口味搜索| 亚洲 欧美一区二区三区| 成人国产综合亚洲| 国产av在哪里看| 精品乱码久久久久久99久播| 亚洲av第一区精品v没综合| 国产麻豆成人av免费视频| 最近最新中文字幕大全电影3| 久久精品国产综合久久久| 999久久久精品免费观看国产| 在线观看一区二区三区| 国产亚洲欧美在线一区二区| 一本久久中文字幕| 午夜福利成人在线免费观看| 99久久99久久久精品蜜桃| 夜夜看夜夜爽夜夜摸| 一a级毛片在线观看| 最近视频中文字幕2019在线8| 丝袜人妻中文字幕| 成人亚洲精品av一区二区| 国产伦精品一区二区三区四那| 久久香蕉精品热| 亚洲成人中文字幕在线播放| or卡值多少钱| 精品欧美国产一区二区三| 欧美中文综合在线视频| 欧美日韩国产亚洲二区| 成人av一区二区三区在线看| 婷婷六月久久综合丁香| 亚洲av熟女| 国产亚洲精品久久久久久毛片| 免费在线观看日本一区| 桃色一区二区三区在线观看| 欧美日韩福利视频一区二区| 真实男女啪啪啪动态图| 欧美不卡视频在线免费观看| 丁香欧美五月| 亚洲av电影在线进入| 国产精品美女特级片免费视频播放器 | 亚洲av成人精品一区久久| 激情在线观看视频在线高清| av中文乱码字幕在线| 亚洲精品粉嫩美女一区| 日本一二三区视频观看| 特大巨黑吊av在线直播| 亚洲中文日韩欧美视频| 色哟哟哟哟哟哟| 亚洲激情在线av| 麻豆av在线久日| 丰满人妻熟妇乱又伦精品不卡| 国产免费男女视频| 欧美+亚洲+日韩+国产| 十八禁网站免费在线| 69av精品久久久久久| 99久久综合精品五月天人人| 99久久99久久久精品蜜桃| 国产午夜福利久久久久久| 久久精品夜夜夜夜夜久久蜜豆| 搡老妇女老女人老熟妇| 中出人妻视频一区二区| 色综合婷婷激情| www.自偷自拍.com| 亚洲精品色激情综合| 久久久久国内视频| 亚洲欧洲精品一区二区精品久久久| 亚洲成av人片免费观看| av黄色大香蕉| 免费av不卡在线播放| 久久中文看片网| 高清毛片免费观看视频网站| 真实男女啪啪啪动态图| 亚洲中文字幕一区二区三区有码在线看 | 国产一区二区在线av高清观看| 欧美乱妇无乱码| 高潮久久久久久久久久久不卡| 亚洲中文字幕日韩| 18禁国产床啪视频网站| a级毛片在线看网站| svipshipincom国产片| 老熟妇乱子伦视频在线观看| 亚洲国产精品sss在线观看| 嫩草影院精品99| 亚洲成a人片在线一区二区| 亚洲国产欧美网| 国产三级在线视频| 久久中文字幕一级| 国产单亲对白刺激| 国产伦在线观看视频一区| 国产三级中文精品| av在线天堂中文字幕| xxxwww97欧美| 国产熟女xx| 在线视频色国产色| 亚洲精品中文字幕一二三四区| 欧美中文日本在线观看视频| 韩国av一区二区三区四区| 久久精品人妻少妇| 国产精品爽爽va在线观看网站| 黄色片一级片一级黄色片| 小蜜桃在线观看免费完整版高清| 日韩欧美 国产精品| 久久精品综合一区二区三区| 美女 人体艺术 gogo| 亚洲国产欧美网| 最好的美女福利视频网| 国内毛片毛片毛片毛片毛片| 日本撒尿小便嘘嘘汇集6| 国产午夜精品论理片| 精品久久久久久,| 性色av乱码一区二区三区2| 精品午夜福利视频在线观看一区| 国产伦在线观看视频一区| 人人妻人人澡欧美一区二区| 中文字幕熟女人妻在线| 国产亚洲欧美98| 一个人看视频在线观看www免费 | 国产精品综合久久久久久久免费| 成年女人永久免费观看视频| 久99久视频精品免费| 少妇的逼水好多| 男人舔女人的私密视频| 九九热线精品视视频播放| e午夜精品久久久久久久| 日韩国内少妇激情av| av欧美777| 18美女黄网站色大片免费观看| 久久久久久久久免费视频了| 国产一级毛片七仙女欲春2| 久久精品影院6| 可以在线观看的亚洲视频| 在线国产一区二区在线| 婷婷精品国产亚洲av在线| 国内毛片毛片毛片毛片毛片| 欧美一级毛片孕妇| 精品福利观看| 宅男免费午夜| 99在线视频只有这里精品首页| 日日夜夜操网爽| 午夜福利在线在线| 听说在线观看完整版免费高清| 老司机午夜十八禁免费视频| 久久精品国产综合久久久| 亚洲乱码一区二区免费版| 亚洲av第一区精品v没综合| 精品99又大又爽又粗少妇毛片 | 国产视频一区二区在线看| 丰满人妻一区二区三区视频av | 久久久久久久久久黄片| 别揉我奶头~嗯~啊~动态视频| 中文字幕高清在线视频| 人人妻,人人澡人人爽秒播| 日本a在线网址| 一本精品99久久精品77| 99久久久亚洲精品蜜臀av| 最近视频中文字幕2019在线8| 精品久久蜜臀av无| 亚洲专区国产一区二区| 中文字幕人妻丝袜一区二区| 熟女人妻精品中文字幕| 亚洲国产日韩欧美精品在线观看 | 国产av在哪里看| cao死你这个sao货| 真实男女啪啪啪动态图| 亚洲av成人不卡在线观看播放网| 在线观看美女被高潮喷水网站 | 国产精品亚洲一级av第二区| 日本精品一区二区三区蜜桃| 精品无人区乱码1区二区| 成人高潮视频无遮挡免费网站| 精品国产三级普通话版| 欧美成狂野欧美在线观看| 99久久99久久久精品蜜桃| av国产免费在线观看| 久久久久国内视频| 一级毛片女人18水好多| 人人妻人人澡欧美一区二区| 69av精品久久久久久| 九色国产91popny在线| 男人舔女人的私密视频| 日本免费a在线| 天堂影院成人在线观看| 亚洲成人久久性| 黄色视频,在线免费观看| 成熟少妇高潮喷水视频| 久久香蕉精品热| 亚洲av日韩精品久久久久久密| 免费在线观看亚洲国产| 老熟妇乱子伦视频在线观看| 欧美黑人欧美精品刺激| aaaaa片日本免费| 亚洲专区国产一区二区| 琪琪午夜伦伦电影理论片6080| 精品久久蜜臀av无| 亚洲中文字幕日韩| 在线观看美女被高潮喷水网站 | 男女之事视频高清在线观看| av黄色大香蕉| 一进一出抽搐gif免费好疼| 免费观看人在逋| 国产综合懂色| 人人妻人人澡欧美一区二区| 老鸭窝网址在线观看| 一级毛片精品| av国产免费在线观看| 嫁个100分男人电影在线观看| 午夜日韩欧美国产| 日本 av在线| 日韩有码中文字幕| 免费无遮挡裸体视频| 老汉色av国产亚洲站长工具| 亚洲成a人片在线一区二区| 在线观看66精品国产| 久久久久久久午夜电影| 美女扒开内裤让男人捅视频| 欧美三级亚洲精品| 极品教师在线免费播放| 成人18禁在线播放| 国产精品久久久久久亚洲av鲁大| 99国产精品一区二区蜜桃av| 黑人欧美特级aaaaaa片| 午夜免费成人在线视频| 一边摸一边抽搐一进一小说| 亚洲av成人一区二区三| 亚洲美女视频黄频| 三级毛片av免费| 在线观看66精品国产| 亚洲国产欧美人成| 欧美国产日韩亚洲一区| 久久久久亚洲av毛片大全| 免费看十八禁软件| 国产熟女xx| 午夜免费观看网址| 亚洲一区高清亚洲精品| 欧美日韩黄片免| 狠狠狠狠99中文字幕| 日日夜夜操网爽| 热99在线观看视频| 亚洲人成电影免费在线| h日本视频在线播放| 国产精品日韩av在线免费观看| 欧美国产日韩亚洲一区| 午夜福利在线在线| 99riav亚洲国产免费| 熟妇人妻久久中文字幕3abv| 黄片小视频在线播放| 性欧美人与动物交配| 1000部很黄的大片| а√天堂www在线а√下载| 国产激情欧美一区二区| 国产精品99久久久久久久久| 99re在线观看精品视频| 免费观看的影片在线观看| 欧美日韩国产亚洲二区| 免费高清视频大片| 法律面前人人平等表现在哪些方面| 午夜福利在线观看免费完整高清在 | 午夜福利在线在线| 嫩草影院入口| 好看av亚洲va欧美ⅴa在| 美女高潮喷水抽搐中文字幕| 亚洲欧美精品综合久久99| 国产精品综合久久久久久久免费| 岛国在线免费视频观看| 亚洲午夜理论影院| 国产精品免费一区二区三区在线| 亚洲第一欧美日韩一区二区三区| 两个人的视频大全免费| 欧美xxxx黑人xx丫x性爽| 韩国av一区二区三区四区| 欧美黑人欧美精品刺激| 精品日产1卡2卡| 国内少妇人妻偷人精品xxx网站 | 亚洲最大成人中文| a级毛片a级免费在线| 夜夜躁狠狠躁天天躁| 一区福利在线观看| 最新在线观看一区二区三区| 国产一级毛片七仙女欲春2| 国产精品久久久久久人妻精品电影| 99国产精品99久久久久| 久久精品国产清高在天天线| 久久精品国产99精品国产亚洲性色| 国产高潮美女av| a在线观看视频网站| 亚洲最大成人中文| 久久久久久大精品| 天堂网av新在线| 热99在线观看视频| 精品久久蜜臀av无| 伊人久久大香线蕉亚洲五| 两性夫妻黄色片| 999久久久精品免费观看国产| 91久久精品国产一区二区成人 | 热99re8久久精品国产| 欧美三级亚洲精品| 精品久久久久久久末码| 国产伦在线观看视频一区| 久久精品综合一区二区三区| xxx96com| 一区福利在线观看| 99久国产av精品| 波多野结衣高清无吗| 国产亚洲精品久久久com| 黑人巨大精品欧美一区二区mp4| 观看美女的网站| 99久久99久久久精品蜜桃| 国产欧美日韩精品亚洲av| 伦理电影免费视频| 国产高清激情床上av| 欧美一级毛片孕妇| 婷婷精品国产亚洲av在线| 在线观看日韩欧美| 欧美三级亚洲精品| 麻豆av在线久日| 欧美日韩瑟瑟在线播放| 精品国产乱码久久久久久男人| 免费一级毛片在线播放高清视频| 在线免费观看的www视频| 国产高清videossex| 久久精品夜夜夜夜夜久久蜜豆| 亚洲av电影不卡..在线观看| 美女免费视频网站| 美女高潮的动态| 亚洲真实伦在线观看| 国产成人av教育| 午夜视频精品福利| 在线观看午夜福利视频| 黄色日韩在线| 男女视频在线观看网站免费| 国产精品国产高清国产av| 91在线精品国自产拍蜜月 | 天堂av国产一区二区熟女人妻| 亚洲精品美女久久av网站| 夜夜躁狠狠躁天天躁| 99热只有精品国产| 久久久久精品国产欧美久久久| 最新中文字幕久久久久 | 久久人妻av系列| 色综合亚洲欧美另类图片| 高清在线国产一区| 久久这里只有精品中国| 在线播放国产精品三级| 久久精品国产亚洲av香蕉五月| 欧洲精品卡2卡3卡4卡5卡区| 久久精品夜夜夜夜夜久久蜜豆| 久久久久国产一级毛片高清牌| 精品电影一区二区在线| 99精品久久久久人妻精品| 欧美日韩综合久久久久久 | 一区二区三区国产精品乱码| 国产一区二区在线观看日韩 | 国产三级黄色录像| 天堂av国产一区二区熟女人妻| 19禁男女啪啪无遮挡网站| 精品一区二区三区av网在线观看| 91麻豆精品激情在线观看国产| av黄色大香蕉| 国产精品,欧美在线| 国产 一区 欧美 日韩| 国产私拍福利视频在线观看| 国产又黄又爽又无遮挡在线| 久久婷婷人人爽人人干人人爱| 校园春色视频在线观看| 五月伊人婷婷丁香| 色av中文字幕| av欧美777| 色尼玛亚洲综合影院| 国产精品久久久人人做人人爽| 欧美国产日韩亚洲一区| 国产蜜桃级精品一区二区三区| 成人18禁在线播放| 热99在线观看视频| 色尼玛亚洲综合影院| 久久天堂一区二区三区四区| 国产 一区 欧美 日韩| 欧美黑人欧美精品刺激| 亚洲在线自拍视频| 国产欧美日韩精品亚洲av| 脱女人内裤的视频| 日韩三级视频一区二区三区| 啦啦啦观看免费观看视频高清| 变态另类丝袜制服| 小蜜桃在线观看免费完整版高清| 中文字幕熟女人妻在线| 国产精品亚洲美女久久久| 亚洲七黄色美女视频| 白带黄色成豆腐渣| 欧美xxxx黑人xx丫x性爽| 久久国产精品影院| 国产精华一区二区三区| 三级国产精品欧美在线观看 | 亚洲熟女毛片儿| 国产亚洲av嫩草精品影院| 我的老师免费观看完整版| 日韩有码中文字幕| av天堂在线播放| 婷婷丁香在线五月| 热99在线观看视频| 制服人妻中文乱码| 国产精品一区二区三区四区久久| 亚洲 欧美一区二区三区| 九九热线精品视视频播放| 欧美日韩国产亚洲二区|