• <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ù)挖掘方法分類對比研究
    精品久久久久久久人妻蜜臀av| 亚洲一级一片aⅴ在线观看| 国产精品电影一区二区三区| 麻豆国产av国片精品| 日韩欧美精品v在线| 边亲边吃奶的免费视频| 极品教师在线视频| 美女 人体艺术 gogo| 男女做爰动态图高潮gif福利片| 欧美一区二区精品小视频在线| 一级黄色大片毛片| 国产探花极品一区二区| 波多野结衣巨乳人妻| 高清毛片免费看| 女人十人毛片免费观看3o分钟| 26uuu在线亚洲综合色| 赤兔流量卡办理| 啦啦啦观看免费观看视频高清| 国产乱人偷精品视频| 日本三级黄在线观看| 天堂影院成人在线观看| 欧美成人精品欧美一级黄| 最近的中文字幕免费完整| 高清毛片免费观看视频网站| 69人妻影院| 两个人视频免费观看高清| 国产一区二区三区在线臀色熟女| 国内揄拍国产精品人妻在线| 久久99精品国语久久久| 男女啪啪激烈高潮av片| 亚洲电影在线观看av| 国产高清不卡午夜福利| 岛国毛片在线播放| 伊人久久精品亚洲午夜| 亚洲性久久影院| 国产午夜精品一二区理论片| 少妇人妻精品综合一区二区 | 搡老妇女老女人老熟妇| 日韩强制内射视频| 国产成年人精品一区二区| 亚洲欧美精品专区久久| 亚洲一区二区三区色噜噜| 日韩成人伦理影院| 国产色爽女视频免费观看| 日本一本二区三区精品| 青春草视频在线免费观看| 免费人成在线观看视频色| 亚洲丝袜综合中文字幕| 免费观看精品视频网站| 国产在视频线在精品| 亚洲欧美精品综合久久99| 亚洲一区高清亚洲精品| 麻豆av噜噜一区二区三区| 男人和女人高潮做爰伦理| 色播亚洲综合网| 国产成人a∨麻豆精品| 亚洲七黄色美女视频| 久久久精品欧美日韩精品| 亚洲av熟女| 一个人观看的视频www高清免费观看| 亚洲内射少妇av| 不卡一级毛片| 22中文网久久字幕| 少妇被粗大猛烈的视频| 亚洲国产欧美人成| 久久国内精品自在自线图片| 国产69精品久久久久777片| 校园春色视频在线观看| 免费观看精品视频网站| 国产精品野战在线观看| 欧美不卡视频在线免费观看| 亚洲欧美精品自产自拍| 搡女人真爽免费视频火全软件| 黑人高潮一二区| 亚洲欧美成人综合另类久久久 | 亚洲经典国产精华液单| 两性午夜刺激爽爽歪歪视频在线观看| 日日摸夜夜添夜夜添av毛片| 成人永久免费在线观看视频| 日韩,欧美,国产一区二区三区 | 国产精品99久久久久久久久| 亚洲av中文字字幕乱码综合| 岛国毛片在线播放| 亚洲中文字幕日韩| 此物有八面人人有两片| 国产精品无大码| 高清日韩中文字幕在线| 少妇的逼好多水| 成人漫画全彩无遮挡| 欧美zozozo另类| 亚洲欧美清纯卡通| 亚洲国产色片| 国产精品一二三区在线看| 国产91av在线免费观看| 乱码一卡2卡4卡精品| 国产日韩欧美在线精品| 18+在线观看网站| 久久欧美精品欧美久久欧美| 色综合色国产| 免费av不卡在线播放| av在线蜜桃| 国产极品天堂在线| 国产不卡一卡二| av又黄又爽大尺度在线免费看 | 熟妇人妻久久中文字幕3abv| 亚洲最大成人中文| 白带黄色成豆腐渣| 看十八女毛片水多多多| 变态另类成人亚洲欧美熟女| 欧美日韩一区二区视频在线观看视频在线 | 久久久精品大字幕| 国产精品永久免费网站| 国产探花在线观看一区二区| 国产精品野战在线观看| 成年女人看的毛片在线观看| 亚洲av男天堂| 一区二区三区四区激情视频 | 亚洲色图av天堂| 国产精品1区2区在线观看.| 99国产极品粉嫩在线观看| 91av网一区二区| 国产高清视频在线观看网站| 青青草视频在线视频观看| 午夜福利在线观看吧| 老司机影院成人| 成人亚洲精品av一区二区| 精品久久久久久久久久免费视频| 免费观看精品视频网站| 久久亚洲国产成人精品v| 国语自产精品视频在线第100页| 久久精品91蜜桃| 神马国产精品三级电影在线观看| 成人三级黄色视频| 99久久精品热视频| 亚洲国产欧洲综合997久久,| 中文字幕熟女人妻在线| 国产精品久久久久久av不卡| 国产精品久久视频播放| 在线观看午夜福利视频| 麻豆国产av国片精品| 亚洲欧美精品综合久久99| 国产在线精品亚洲第一网站| 成人午夜精彩视频在线观看| 国产伦理片在线播放av一区 | 中文亚洲av片在线观看爽| 成人毛片a级毛片在线播放| 波多野结衣高清无吗| 亚洲精品成人久久久久久| 高清午夜精品一区二区三区 | 五月伊人婷婷丁香| 久久久久久久久大av| 变态另类丝袜制服| 嫩草影院新地址| 卡戴珊不雅视频在线播放| 禁无遮挡网站| 国产成人aa在线观看| 99热这里只有是精品在线观看| 晚上一个人看的免费电影| 日本在线视频免费播放| 国产成人a区在线观看| 国产成人午夜福利电影在线观看| 中文资源天堂在线| 亚洲最大成人中文| 人妻夜夜爽99麻豆av| 久久99蜜桃精品久久| 黄片无遮挡物在线观看| 午夜a级毛片| 色吧在线观看| 91狼人影院| 午夜久久久久精精品| 国产激情偷乱视频一区二区| 欧美极品一区二区三区四区| a级一级毛片免费在线观看| 国产成人freesex在线| 国产毛片a区久久久久| 国产高清不卡午夜福利| 欧美+亚洲+日韩+国产| 日韩欧美一区二区三区在线观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 久久久色成人| 国产老妇伦熟女老妇高清| 尤物成人国产欧美一区二区三区| 亚洲三级黄色毛片| 国产亚洲精品av在线| 久久精品91蜜桃| 村上凉子中文字幕在线| 高清日韩中文字幕在线| 一个人观看的视频www高清免费观看| 午夜激情福利司机影院| 国产乱人偷精品视频| 免费一级毛片在线播放高清视频| 国产精品久久电影中文字幕| 在线国产一区二区在线| 色综合色国产| 久久久午夜欧美精品| 国产成人精品婷婷| 免费人成视频x8x8入口观看| 国产麻豆成人av免费视频| 国产成人精品一,二区 | 国产精品一区二区三区四区久久| 综合色丁香网| 美女cb高潮喷水在线观看| 少妇熟女欧美另类| 中文在线观看免费www的网站| 国产成人一区二区在线| 九九在线视频观看精品| 日产精品乱码卡一卡2卡三| 亚洲综合色惰| 少妇裸体淫交视频免费看高清| 六月丁香七月| 亚洲一区高清亚洲精品| 亚洲人成网站在线观看播放| 国产一级毛片七仙女欲春2| 免费观看a级毛片全部| 亚洲av成人av| 毛片一级片免费看久久久久| 麻豆av噜噜一区二区三区| 天堂av国产一区二区熟女人妻| av福利片在线观看| 深夜a级毛片| 少妇裸体淫交视频免费看高清| 国产美女午夜福利| 十八禁国产超污无遮挡网站| 人妻制服诱惑在线中文字幕| 91午夜精品亚洲一区二区三区| 久久精品国产清高在天天线| 男女视频在线观看网站免费| 婷婷色av中文字幕| 亚洲精华国产精华液的使用体验 | 成人特级av手机在线观看| 午夜老司机福利剧场| 免费看av在线观看网站| 可以在线观看毛片的网站| 女同久久另类99精品国产91| 国内揄拍国产精品人妻在线| 欧美性猛交╳xxx乱大交人| 午夜爱爱视频在线播放| 亚洲天堂国产精品一区在线| 国产午夜精品论理片| 日本黄大片高清| 婷婷六月久久综合丁香| 深爱激情五月婷婷| 成人性生交大片免费视频hd| 久久久a久久爽久久v久久| 伦理电影大哥的女人| av在线观看视频网站免费| 国产一区亚洲一区在线观看| 亚洲精品影视一区二区三区av| 秋霞在线观看毛片| 精品久久久久久久久久免费视频| 欧美最新免费一区二区三区| 国产亚洲5aaaaa淫片| 床上黄色一级片| 亚洲欧美清纯卡通| 乱码一卡2卡4卡精品| 日韩一区二区视频免费看| 亚洲欧美精品自产自拍| 久久99热这里只有精品18| 嫩草影院新地址| 夜夜夜夜夜久久久久| 亚洲av成人精品一区久久| 亚洲第一区二区三区不卡| 久久国产乱子免费精品| 国产黄色小视频在线观看| 91精品国产九色| 亚洲在久久综合| ponron亚洲| 免费不卡的大黄色大毛片视频在线观看 | 久久精品综合一区二区三区| 嫩草影院入口| 免费看av在线观看网站| 美女被艹到高潮喷水动态| 国内精品美女久久久久久| 搡女人真爽免费视频火全软件| 激情 狠狠 欧美| 午夜老司机福利剧场| 欧美极品一区二区三区四区| 欧美日韩国产亚洲二区| 国产高潮美女av| av国产免费在线观看| 日本av手机在线免费观看| 熟女人妻精品中文字幕| av在线老鸭窝| 亚洲国产精品合色在线| 亚洲乱码一区二区免费版| 欧美性猛交╳xxx乱大交人| 最近的中文字幕免费完整| 日韩一区二区视频免费看| 不卡一级毛片| 长腿黑丝高跟| 日韩亚洲欧美综合| 精品国内亚洲2022精品成人| 国产乱人视频| 成年女人永久免费观看视频| 欧美在线一区亚洲| 免费黄网站久久成人精品| 大香蕉久久网| 九九久久精品国产亚洲av麻豆| 免费人成视频x8x8入口观看| 亚洲精品乱码久久久v下载方式| 99视频精品全部免费 在线| 最好的美女福利视频网| 超碰av人人做人人爽久久| 夜夜夜夜夜久久久久| 中文资源天堂在线| 日韩中字成人| 99久久无色码亚洲精品果冻| 直男gayav资源| 免费无遮挡裸体视频| 精品久久久久久久人妻蜜臀av| 搞女人的毛片| 日韩欧美精品v在线| 免费看美女性在线毛片视频| 性欧美人与动物交配| av在线亚洲专区| 亚洲精品成人久久久久久| 观看美女的网站| 蜜桃亚洲精品一区二区三区| 波多野结衣巨乳人妻| 国产乱人视频| 国产精品日韩av在线免费观看| 亚洲精品国产av成人精品| 天堂影院成人在线观看| 小说图片视频综合网站| 成人国产麻豆网| videossex国产| 免费黄网站久久成人精品| 天天一区二区日本电影三级| 中文欧美无线码| a级毛片免费高清观看在线播放| 久久午夜亚洲精品久久| 欧美日本视频| 日本成人三级电影网站| 久久亚洲精品不卡| 岛国在线免费视频观看| 一卡2卡三卡四卡精品乱码亚洲| 欧美最新免费一区二区三区| 丰满的人妻完整版| 99久久无色码亚洲精品果冻| 国产精品永久免费网站| 在线观看66精品国产| 日韩精品有码人妻一区| 精品久久久久久久末码| 国产久久久一区二区三区| 亚洲国产高清在线一区二区三| 久久亚洲精品不卡| 18禁在线无遮挡免费观看视频| 亚洲久久久久久中文字幕| 欧美区成人在线视频| 国产精品免费一区二区三区在线| 好男人视频免费观看在线| 成人美女网站在线观看视频| 亚洲中文字幕一区二区三区有码在线看| 丰满的人妻完整版| 一级毛片我不卡| 国产色爽女视频免费观看| 12—13女人毛片做爰片一| 久久精品综合一区二区三区| 国产黄色视频一区二区在线观看 | 久久精品国产亚洲网站| 看十八女毛片水多多多| 国产精品久久久久久精品电影| 午夜福利视频1000在线观看| 免费人成在线观看视频色| 日韩人妻高清精品专区| 免费观看人在逋| 亚洲欧美精品综合久久99| 午夜福利在线观看免费完整高清在 | 国产av在哪里看| 校园人妻丝袜中文字幕| 一级av片app| 色哟哟哟哟哟哟| 一级毛片久久久久久久久女| 成人亚洲精品av一区二区| 色哟哟哟哟哟哟| 日韩强制内射视频| 国产精品久久久久久久久免| 婷婷精品国产亚洲av| av又黄又爽大尺度在线免费看 | 少妇丰满av| 国产成人影院久久av| 青春草视频在线免费观看| 久久精品人妻少妇| 成人漫画全彩无遮挡| 性欧美人与动物交配| 色综合色国产| 日日撸夜夜添| 日韩亚洲欧美综合| 国产成人午夜福利电影在线观看| 亚洲av中文av极速乱| 99九九线精品视频在线观看视频| 欧美丝袜亚洲另类| 国产成人一区二区在线| 久久精品国产亚洲av香蕉五月| 一边摸一边抽搐一进一小说| 老师上课跳d突然被开到最大视频| 国产午夜福利久久久久久| 亚洲无线在线观看| 精品久久久久久成人av| 国产午夜福利久久久久久| 国产精品日韩av在线免费观看| 麻豆乱淫一区二区| 天天一区二区日本电影三级| 99久久中文字幕三级久久日本| 一本一本综合久久| 免费一级毛片在线播放高清视频| 欧美日本亚洲视频在线播放| 毛片一级片免费看久久久久| 亚洲三级黄色毛片| 深夜a级毛片| 午夜免费激情av| 欧美日本视频| 国产高清激情床上av| 1024手机看黄色片| 尾随美女入室| www.色视频.com| 成人av在线播放网站| 色5月婷婷丁香| 国产私拍福利视频在线观看| 一进一出抽搐动态| 看十八女毛片水多多多| 国产一区二区激情短视频| 亚洲av电影不卡..在线观看| 欧美区成人在线视频| 九色成人免费人妻av| 欧美成人精品欧美一级黄| 亚洲中文字幕一区二区三区有码在线看| 精品久久国产蜜桃| 精品久久久久久成人av| 欧美精品一区二区大全| or卡值多少钱| 成人二区视频| 国产免费男女视频| 久久精品久久久久久噜噜老黄 | 国产一区二区在线观看日韩| 岛国在线免费视频观看| 99久久成人亚洲精品观看| 色哟哟哟哟哟哟| 综合色丁香网| 精品久久久久久久久久免费视频| 12—13女人毛片做爰片一| 午夜免费男女啪啪视频观看| 久久精品国产亚洲av涩爱 | 久久久a久久爽久久v久久| 国产精品av视频在线免费观看| 一级黄片播放器| 麻豆久久精品国产亚洲av| 亚洲精品影视一区二区三区av| 日日摸夜夜添夜夜添av毛片| 亚洲内射少妇av| 国产成人a区在线观看| 男女做爰动态图高潮gif福利片| 女人被狂操c到高潮| 国产亚洲5aaaaa淫片| 精品人妻熟女av久视频| 成年女人看的毛片在线观看| 丝袜美腿在线中文| 国产精品野战在线观看| 久久久成人免费电影| 麻豆av噜噜一区二区三区| 天天躁夜夜躁狠狠久久av| av卡一久久| 欧美激情久久久久久爽电影| 我要看日韩黄色一级片| 国产探花在线观看一区二区| 国产精品一及| 亚洲av电影不卡..在线观看| 亚洲av成人精品一区久久| 亚洲乱码一区二区免费版| 在线免费观看不下载黄p国产| 人体艺术视频欧美日本| 搡老妇女老女人老熟妇| 国产精品久久电影中文字幕| 寂寞人妻少妇视频99o| 高清毛片免费观看视频网站| 成人高潮视频无遮挡免费网站| 国产av不卡久久| 亚洲精品乱码久久久v下载方式| 亚洲乱码一区二区免费版| 人妻少妇偷人精品九色| 欧美日韩综合久久久久久| 欧美性感艳星| 国产色婷婷99| 菩萨蛮人人尽说江南好唐韦庄 | 久久人人爽人人爽人人片va| 尤物成人国产欧美一区二区三区| 免费观看a级毛片全部| 午夜激情欧美在线| 国内少妇人妻偷人精品xxx网站| 简卡轻食公司| 日韩国内少妇激情av| a级毛片免费高清观看在线播放| 国产成人午夜福利电影在线观看| 欧美一区二区国产精品久久精品| 麻豆一二三区av精品| 国产精品日韩av在线免费观看| 99精品在免费线老司机午夜| 哪里可以看免费的av片| 国产综合懂色| www.色视频.com| 久久草成人影院| 深夜a级毛片| 免费看日本二区| 一进一出抽搐gif免费好疼| 免费观看的影片在线观看| 国产午夜福利久久久久久| 成人无遮挡网站| 亚洲性久久影院| 国产成人精品婷婷| 亚洲美女搞黄在线观看| 日韩中字成人| 99久国产av精品| 日本色播在线视频| 精品日产1卡2卡| 亚洲精品久久久久久婷婷小说 | 丰满乱子伦码专区| 国产成人午夜福利电影在线观看| 久久久久久久午夜电影| 日韩一区二区视频免费看| 亚洲精品色激情综合| 国内久久婷婷六月综合欲色啪| 国产av在哪里看| 久99久视频精品免费| 亚洲电影在线观看av| 亚洲五月天丁香| 黄色一级大片看看| 日日摸夜夜添夜夜爱| 国产探花极品一区二区| 美女xxoo啪啪120秒动态图| 中文字幕av在线有码专区| 哪里可以看免费的av片| 精品久久久噜噜| 久久人人爽人人爽人人片va| 久久久久久久亚洲中文字幕| 在线免费观看的www视频| 亚洲天堂国产精品一区在线| 国产精品久久久久久精品电影| 国产日本99.免费观看| 在线观看av片永久免费下载| 麻豆国产av国片精品| 日本-黄色视频高清免费观看| 春色校园在线视频观看| 网址你懂的国产日韩在线| 99九九线精品视频在线观看视频| 天堂av国产一区二区熟女人妻| 高清日韩中文字幕在线| 久久久久久久久久久丰满| 搡女人真爽免费视频火全软件| 国产精品伦人一区二区| 少妇人妻一区二区三区视频| 免费在线观看成人毛片| 2022亚洲国产成人精品| 亚洲国产精品成人久久小说 | 黄色视频,在线免费观看| 久久久久久久久中文| 一级毛片电影观看 | 人妻少妇偷人精品九色| 人妻夜夜爽99麻豆av| 99在线视频只有这里精品首页| 成人特级黄色片久久久久久久| 高清毛片免费观看视频网站| 全区人妻精品视频| 亚洲va在线va天堂va国产| 人体艺术视频欧美日本| 一边摸一边抽搐一进一小说| 国产精品一区二区在线观看99 | 亚洲欧美精品综合久久99| 男女下面进入的视频免费午夜| 日本在线视频免费播放| 免费不卡的大黄色大毛片视频在线观看 | 国产av一区在线观看免费| 精品久久久久久久久亚洲| 免费看日本二区| 日本色播在线视频| 国产精品伦人一区二区| 国产午夜精品久久久久久一区二区三区| 欧美3d第一页| 亚洲成人中文字幕在线播放| 99热这里只有是精品在线观看| 亚洲精品自拍成人| 美女 人体艺术 gogo| 日韩精品青青久久久久久| 精品无人区乱码1区二区| 亚洲欧美精品专区久久| 极品教师在线视频| 草草在线视频免费看| 桃色一区二区三区在线观看| 国产日韩欧美在线精品| 亚洲精品粉嫩美女一区| 国产成人一区二区在线| 欧美在线一区亚洲| 在线观看免费视频日本深夜| 99国产精品一区二区蜜桃av| 亚洲色图av天堂| 天天躁日日操中文字幕| 国产精品国产高清国产av| 美女内射精品一级片tv| 欧美日韩一区二区视频在线观看视频在线 | 国产精品久久久久久久久免| www日本黄色视频网| 日本在线视频免费播放| 色吧在线观看| 看片在线看免费视频| 欧美一区二区亚洲| 搞女人的毛片| 少妇猛男粗大的猛烈进出视频 | 美女脱内裤让男人舔精品视频 | 国内久久婷婷六月综合欲色啪| 成人国产麻豆网| 日本在线视频免费播放| 久久人人爽人人爽人人片va| 可以在线观看毛片的网站| 我要看日韩黄色一级片| 不卡视频在线观看欧美| 欧美zozozo另类|