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

    AA染色體組野生種花生Ty3-gypsy類逆轉(zhuǎn)座子逆轉(zhuǎn)錄酶序列的多樣性

    2022-06-09 22:36:23劉俊仙劉菁唐榮華陽太億韓柱強(qiáng)唐秀梅賀梁瓊鐘瑞春蔣菁黃志鵬吳海寧韋榮昌熊發(fā)前
    江蘇農(nóng)業(yè)科學(xué) 2022年9期
    關(guān)鍵詞:逆轉(zhuǎn)錄酶多樣性

    劉俊仙 劉菁 唐榮華 陽太億 韓柱強(qiáng) 唐秀梅 賀梁瓊 鐘瑞春 蔣菁 黃志鵬 吳海寧 韋榮昌 熊發(fā)前

    摘要:以2份AA染色體組野生種花生為試驗(yàn)材料,擴(kuò)增分離Ty3-gypsy類逆轉(zhuǎn)座子逆轉(zhuǎn)錄酶序列,分析其序列特征、多樣性及進(jìn)化關(guān)系。利用根據(jù)Ty3-gypsy類逆轉(zhuǎn)座子逆轉(zhuǎn)錄酶的保守區(qū)設(shè)計(jì)的簡并引物進(jìn)行PCR擴(kuò)增,對目的條帶進(jìn)行回收、克隆和測序后,對逆轉(zhuǎn)錄酶序列進(jìn)行生物信息學(xué)分析。目的條帶大小均約為430 bp,分別從2份AA染色體組野生種花生材料中分離到32、33條逆轉(zhuǎn)錄酶序列,對于65條逆轉(zhuǎn)錄酶序列,其長度變化范圍為 397~440 bp,A+T所占比例范圍為56.48%~68.14%,(A+T)/(G+C)為1.3~2.13,核苷酸序列間相似性范圍為 62.6%~97.9%;65條逆轉(zhuǎn)錄酶序列被劃分為9個(gè)家族,其中家族Ⅰ為主要成分;翻譯成氨基酸后,序列間相似性范圍為12.4%~98.6%,相比核苷酸序列,氨基酸序列表現(xiàn)出更高的異質(zhì)性;65條逆轉(zhuǎn)錄酶序列中有26條發(fā)生了無義突變,2份野生種花生材料的無義突變發(fā)生率相當(dāng);序列間保守基序大體一致,但也發(fā)生了一定的變異,呈現(xiàn)出一定的異質(zhì)性;9個(gè)家族所選代表序列的蛋白質(zhì)結(jié)構(gòu)總體類似,但也在螺旋結(jié)構(gòu)數(shù)、折疊結(jié)構(gòu)數(shù)、轉(zhuǎn)角數(shù)、氫鍵數(shù)、α-螺旋數(shù)、β-折疊數(shù)上存在差別;系統(tǒng)進(jìn)化樹顯示,所有逆轉(zhuǎn)錄酶序列被分為13類,A類和B類中包含大部分逆轉(zhuǎn)錄酶序列,E類中的5條AA染色體組野生種花生與其他物種植物的逆轉(zhuǎn)錄酶序列具有較高相似性,表明這些序列在進(jìn)化過程中有可能發(fā)生了橫向傳遞;通過比對花生EST數(shù)據(jù)庫,發(fā)現(xiàn)3條來自Archis duranensis(PI262133)和2條來自A.duranensis(PI219823)的Ty3-gypsy類逆轉(zhuǎn)座子具有轉(zhuǎn)錄活性,本研究為下一步分離Ty3-gypsy類逆轉(zhuǎn)座子全長序列、研究其轉(zhuǎn)錄轉(zhuǎn)座活性和功能提供序列基礎(chǔ),也為基于Ty3-gypsy類逆轉(zhuǎn)座子的花生屬分子標(biāo)記開發(fā)奠定基礎(chǔ)。

    關(guān)鍵詞:野生種花生;Ty3-gypsy類逆轉(zhuǎn)座子;逆轉(zhuǎn)錄酶;多樣性

    中圖分類號:S565.201 ??文獻(xiàn)標(biāo)志碼: A

    文章編號:1002-1302(2022)09-0043-12

    花生被譽(yù)為“長生果”,是世界上最重要的油料作物和經(jīng)濟(jì)作物之一,也是最重要的植物食用油來源之一[1]。我國是世界花生生產(chǎn)第一大國,種植面積達(dá)500萬hm2/年以上,約占全球的20%;我國花生產(chǎn)量達(dá)1 700萬t/年,約占全球的40%,居全球首位,花生占我國油料作物總產(chǎn)(不包括大豆)的46%;花生油約占國產(chǎn)植物油的25%,僅次于菜籽油,是國產(chǎn)植物油的第二大來源[2]。

    轉(zhuǎn)座子是指能從同一條染色體或不同染色體上的一個(gè)位點(diǎn)轉(zhuǎn)移到另外一個(gè)位點(diǎn)的可移動DNA序列[3]。轉(zhuǎn)座子分為逆轉(zhuǎn)座子和DNA轉(zhuǎn)座子,反轉(zhuǎn)錄轉(zhuǎn)座子分為長末端重復(fù)序列(LTR)逆轉(zhuǎn)座子和非LTR逆轉(zhuǎn)座子,Ty1-copia和Ty3-gypsy是LTR逆轉(zhuǎn)座子的主要兩大類型[4-6],這兩大類LTR逆轉(zhuǎn)座子的逆轉(zhuǎn)錄酶序列可以利用簡并PCR技術(shù)擴(kuò)增出來。LTR逆轉(zhuǎn)座子非常適合用來開發(fā)成分子標(biāo)記,序列特異擴(kuò)增多態(tài)性(S-SAP)、逆轉(zhuǎn)座子位點(diǎn)間擴(kuò)增多態(tài)性(IRAP)和逆轉(zhuǎn)座子-微衛(wèi)星擴(kuò)增多態(tài)性(REMAP)是目前基于LTR逆轉(zhuǎn)座子的主要分子標(biāo)記技術(shù)[7-8]。

    MITEs轉(zhuǎn)座子已在花生上被較廣泛地研究和應(yīng)用[9-23],而花生LTR逆轉(zhuǎn)座子的研究報(bào)道較少[24-25]。Nielen等先后分離出花生Ty3-gypsy類逆轉(zhuǎn)座子的FIDEL和Ty1-copia類逆轉(zhuǎn)座子的Matita,對其特性和作用進(jìn)行分析[24-25]。

    筆者曾系統(tǒng)對花生LTR逆轉(zhuǎn)座子和MITE轉(zhuǎn)座子的分離及其應(yīng)用的研究進(jìn)展進(jìn)行了歸納[1],也對四倍體野生種花生(Arachis monticola)的Ty1-copia類逆轉(zhuǎn)座子逆轉(zhuǎn)錄酶基因進(jìn)行克隆與分析[26],但尚未見對AA染色體組野生種花生Ty3-gypsy類逆轉(zhuǎn)座子逆轉(zhuǎn)錄酶序列進(jìn)行分離和多樣性分析的報(bào)道。

    本研究擬分離AA染色體組野生種花生A.duranensis的Ty3-gypsy類逆轉(zhuǎn)座子逆轉(zhuǎn)錄酶序列,分析其序列特征和多樣性,弄清其序列組成和變異模式及其與其他物種植物之間的系統(tǒng)發(fā)育進(jìn)化關(guān)系,為進(jìn)一步分離其全長序列、研究其轉(zhuǎn)錄轉(zhuǎn)座活性和功能提供序列基礎(chǔ),也為基于Ty3-gypsy類逆轉(zhuǎn)座子的花生屬分子標(biāo)記開發(fā)和應(yīng)用奠定基礎(chǔ)。

    1 材料與方法

    1.1 試驗(yàn)材料

    2份AA染色體組野生種花生材料分別為A.duranensis(PI262133)和A.duranensis(PI219823)。

    1.2 方法

    1.2.1 DNA提取 采用改良過的CTAB法[27]提取高質(zhì)量花生基因組DNA。

    1.2.2 逆轉(zhuǎn)錄酶序列的PCR擴(kuò)增 Ty3-gypsy類逆轉(zhuǎn)座子逆轉(zhuǎn)錄酶序列的PCR擴(kuò)增采用前人設(shè)計(jì)的簡并引物,上下游引物序列分別為:Gyrt1,5′-AGMGRTATGTGYGTSGAYTAT-3′;Gyrt2,5′-CAMCCMRAAMWCACAMTT-3′。其中,R=A/G,Y=C/T,M=A/C,S=C/G,W=A/T,N=A/T/C/G[27]。PCR擴(kuò)增反應(yīng)體系、擴(kuò)增程序以及PCR產(chǎn)物的分離檢測參考文獻(xiàn)[26]。

    1.2.3 目的條帶的克隆和測序 具體操作參考文獻(xiàn)[26]。

    1.2.4 逆轉(zhuǎn)錄酶序列分析 序列相似性檢索、序列統(tǒng)計(jì)分析、序列圖及Logo圖生成、蛋白質(zhì)的二級結(jié)構(gòu)和三級結(jié)構(gòu)預(yù)測、蛋白質(zhì)三級結(jié)構(gòu)的轉(zhuǎn)角數(shù)和氫鍵數(shù)統(tǒng)計(jì)、保守基序預(yù)測等參考文獻(xiàn)[26]。運(yùn)用MEGA 6.0軟件的鄰接法(No. of differences模型)構(gòu)建系統(tǒng)進(jìn)化樹,自展值設(shè)置為1 000。與花生EST數(shù)據(jù)庫比對,鑒定具有轉(zhuǎn)錄活性的Ty3-gypsy類反轉(zhuǎn)錄轉(zhuǎn)座子。B5F2FA2E-8CD4-41A0-93B1-4A7CB397A425

    2 結(jié)果與分析

    2.1 逆轉(zhuǎn)錄酶序列的PCR擴(kuò)增及測序

    對2份AA染色體組野生種花生材料的基因組DNA進(jìn)行PCR擴(kuò)增,擴(kuò)增結(jié)果(圖1)顯示,2份材料均擴(kuò)增出了1條約430 bp大小的目的條帶。對目的條帶進(jìn)行回收、克隆和測序,從A.duranensis(PI262133)和A.duranensis(PI219823)中均獲得了34條序列。去除重復(fù)序列和非目標(biāo)序列后,分別從A.duranensis(PI262133)和A.duranensis(PI219823)中獲得了32、33條逆轉(zhuǎn)錄酶序列(分別命名為AdRT3-X和AdRT4-X),對65條序列進(jìn)行多重比對分析并生成序列l(wèi)ogo(圖2、圖3)。

    2.2 逆轉(zhuǎn)錄酶序列分析

    所有序列長度都在397~440 bp之間。在A.duranensis(PI262133)的32條序列中,AdRT3-23的序列長度最短,為397 bp,AdRT3-27的序列長度最長,為440 bp,有26條序列的長度均為432 bp,占所克隆序列的81.25%(表1);A、T、C、G數(shù)量變化范圍分別為109~140、107~153、56~86、71~114個(gè),A+T所占比例范圍為56.48%~65.97%,A+T與G+C比例為1.3~1.94(表1);核苷酸序列間相似性范圍為63.7%~99.3%,其中AdRT3-15與AdRT3-24以及AdRT3-20與AdRT3-30之間的相似性最高,達(dá)99.3%,AdRT3-12 與AdRT3-14之間的相似性最低,為63.7%;氨基酸序列間相似性范圍為10.9%~99.3%(表2)。在A.duranensis(PI219823)的33條序列中,AdRT4-27的序列長度最短,為428 bp,AdRT4-28的序列長度最長,為440 bp,有27條序列的長度均為432 bp,占所克隆序列的81.82%(1);A、T、C、G數(shù)量變化范圍分別為117~139、108~157、51~85、77~114個(gè),A+T所占比例范圍為56.48%~68.14%,(A+T)/(G+C)為1.3~2.13(表1);核苷酸序列間相似性范圍為62.8%~99.3%,其中,AdRT4-5與AdRT4-28之間的相似性最高,達(dá)99.3%,AdRT4-1 與AdRT4-23之間的相似性最低,為62.8%,氨基酸序列間相似性范圍為31.5%~99.3%(表2);將2份花生材料的65條序列合并進(jìn)行分析,核苷酸序列間相似性范圍為62.6%~97.9%,氨基酸序列間相似性范圍為12.4%~98.6%。

    2.3 逆轉(zhuǎn)錄酶序列聚類分析

    遺傳進(jìn)化樹(圖4)顯示,65條序列被劃分為9個(gè)家族。家族Ⅰ包含20條序列,占總序列數(shù)的30.80%,該家族只有2條序列來自A.duranensis(PI262133);家族Ⅱ包含7條序列,只有1條來自A.duranensis(PI219823);家族Ⅲ和家族Ⅳ均只包含1條來自A.duranensis(PI262133)的序列,這2條序列與家族Ⅰ親緣關(guān)系較遠(yuǎn)而單獨(dú)聚為一類;家族Ⅴ包含7條序列,有1條來自A.duranensis(PI219823);家族Ⅵ包含9條序列,有2條來自A.duranensis(PI262133);家族Ⅶ包含7條來自A.duranensis(PI262133)和1條來自A.duranensis(PI219823)的序列;家族Ⅷ包含5條序列;家族Ⅸ包含7條序列;除家族Ⅲ和家族Ⅳ外,其他每個(gè)家族中都存在來自這2份野生種花生材料的序列。

    2.4 逆轉(zhuǎn)錄酶氨基酸序列分析

    65條序列中有26條存在無義突變,其中有12條序列來自A.duranensis(PI262133),占該材料序列總數(shù)的37.5%,有14條序列來自A.duranensis(PI219823),占該材料序列總數(shù)的42.42%(圖5)。無義突變在2份野生種花生材料序列中的具體表現(xiàn)為AdRT3-7突變數(shù)最多,存在14個(gè),分別在第36、第38、第49、第51、第63、第76、第87、第95、第103、第109、第113、第128、第136、第138個(gè)氨基酸處;AdRT4-23存在12個(gè)無義突變,分別在第32、第59、第72、第77、第87、第88、第96、第100、第108、第121、第128、第138個(gè)氨基酸處;AdRT3-3存在9個(gè)無義突變,分別在第30、第77、第78、第88、第100、第104、第108、第121、第128個(gè)氨基酸處;AdRT4-27存在9個(gè)無義突變,分別在第32、第38、第42、第99、第103、第113、第120、第122、第127個(gè)氨基酸處;AdRT4-11存在8個(gè)無義突變,分別在第8、第77、第88、第101、第104、第108、第121、第128個(gè)氨基酸處;AdRT4-30存在8個(gè)無義突變,分別在第8、第77、第88、第101、第104、第108、第121、第128個(gè)氨基酸處;AdRT3-26存在7個(gè)無義突變,分別在第第52、第54、第80、第91、第111、第124、第131個(gè)氨基酸處;AdRT3-33存在5個(gè)無義突變,分別在第96、第101、第104、第121、第128個(gè)氨基酸處;AdRT3-23存在4個(gè)無義突變,分別在第86、第93、第102、第127個(gè)氨基酸處;AdRT3-6存在4個(gè)無義突變,分別在第98、第105、第114、第139個(gè)氨基酸處;AdRT4-5(第106、第112、第136個(gè)氨基酸處)、AdRT4-17(第62、第105、第106個(gè)氨基酸處)和AdRT4-28(第106、第112、第136個(gè)氨基酸處)各存在3個(gè)無義突變;AdRT3-14(第9、第115個(gè)氨基酸處)和AdRT4-14(第32、第129個(gè)氨基酸處)各存在2個(gè)無義突變;AdRT3-11(第48個(gè)氨基酸處)、AdRT3-13(第43個(gè)氨基酸處)、AdRT3-17(第130個(gè)氨基酸處)、AdRT3-21(第32個(gè)氨基酸處)、AdRT3-32(第30個(gè)氨基酸處)、AdRT4-12(第131個(gè)氨基酸處)、AdRT4-15(第30個(gè)氨基酸處)、AdRT4-16(第48個(gè)氨基酸處)、AdRT4-18(第9個(gè)氨基酸處)、AdRT4-29(第43個(gè)氨基酸處)、AdRT4-36(第45個(gè)氨基酸處)各存在1個(gè)無義突變;部分序列存在連續(xù)無義突變的現(xiàn)象。無義突變導(dǎo)致逆轉(zhuǎn)座子失去轉(zhuǎn)錄活性。B5F2FA2E-8CD4-41A0-93B1-4A7CB397A425

    2.5 逆轉(zhuǎn)錄酶的蛋白結(jié)構(gòu)預(yù)測

    翻譯成氨基酸后,根據(jù)核苷酸聚類結(jié)果,選擇2份AA染色體組野生種花生材料中每個(gè)家族中的代表序列各1條,家族Ⅲ和家族Ⅳ均只包含1條序列,其選擇均來自A.duranensis(PI262133)。利用在線軟件Phyre2預(yù)測AA染色體組野生種花生Ty3-gypsy類逆轉(zhuǎn)座子逆轉(zhuǎn)錄酶蛋白質(zhì)的二級結(jié)構(gòu)和三級結(jié)構(gòu),代表序列蛋白三級結(jié)構(gòu)匹配覆蓋度最高的模板為c2opqA、c5dmqA、d2zd1b1和c3kk1B,置信度均為100,都屬于逆轉(zhuǎn)錄酶家族。二級結(jié)構(gòu)包含4~6個(gè)α-螺旋和6~9個(gè)β-折疊;三級結(jié)構(gòu)包含 14~19個(gè)轉(zhuǎn)角、70~96個(gè)氫鍵,還存在6個(gè)明顯的螺旋結(jié)構(gòu)和7個(gè)明顯的折疊結(jié)構(gòu)(表3)。其中,家族Ⅰ代表序列AdRT3-10和AdRT4-18,其蛋白質(zhì)二級結(jié)構(gòu)和三級結(jié)構(gòu)見圖6。

    2.6 逆轉(zhuǎn)錄酶保守基序預(yù)測

    65條序列共存在10種保守基序(圖7),其中有57條序列包含motif 1,占序列總數(shù)的87.69%;除AdRT3-26外,剩余的64條序列均包含motif 2,占序列總數(shù)的98.46%;52條序列包含motif 3,占序列總數(shù)的80%;62條序列包含motif 4,占序列總數(shù)的95.38%;4條序列包含motif 5,5條序列包含motif 6,4條序列包含motif 7,3條序列包含motif 8,7條序列包含motif 9,3條序列包含motif 10。AdRT3-26序列存在7個(gè)終止密碼子突變,只包含motif 5和motif 9,該序列上游部分無保守基序且與其他序列差異較大,在系統(tǒng)進(jìn)化樹中與其他序列遺傳距離最遠(yuǎn),單獨(dú)歸為一類;AdRT3-3、AdRT4-11和AdRT4-30同時(shí)包含motif 5和motif 7,在系統(tǒng)進(jìn)化樹中為同一類;AdRT3-6、AdRT3-7、AdRT3-23同時(shí)包含motif 6和motif 8,在系統(tǒng)進(jìn)化樹中為同一類。

    2.7 逆轉(zhuǎn)錄酶序列系統(tǒng)進(jìn)化樹構(gòu)建

    從NCBI數(shù)據(jù)庫中下載已登錄的來源于其他物種植物的Ty3-gypsy類逆轉(zhuǎn)座子逆轉(zhuǎn)錄酶氨基酸序列(表4),與本研究中克隆所獲得的65條逆轉(zhuǎn)錄酶序列一起構(gòu)建系統(tǒng)進(jìn)化樹,分析物種以及序列之間的系統(tǒng)進(jìn)化關(guān)系(圖8)。進(jìn)化樹顯示,所有逆轉(zhuǎn)錄酶序列可分為13類。A類和B類中包含大部分AA染色體組野生種花生Ty3-gypsy類逆轉(zhuǎn)座子逆轉(zhuǎn)錄酶序列,不包含來自其他物種植物的逆轉(zhuǎn)錄酶序列;C類只有AdRT3-17,D類也只有AdRT4-7。E類包含5條來自A.duranensis(PI262133)和1條來自A.duranensis(PI219823)的逆轉(zhuǎn)錄酶序列,其與來自綠豆(Vigna radiata,AAT85841.1)、黃瓜(Cucumis sativus,ADD83121.1)、山桑(Morus bombycis,BAB40830.1)、落葉松(Larix gmelinii,BAQ22332.1)、牡丹(Paeonia suffruticosa,AFQ94056.1)、火龍果(Hylocereus undatus,AOS58468.1)、李子(Prunus salicina,AGX45501.1)、油棕(Elaeis guineensis,CAD45567.1)、擬南芥(Arabidopsis thaliana,BAB40828.1)、棗(Ziziphus jujuba,AFR43612.1)、白皮松(Phelipanche bungeana,ABD43118.1)、多花黑麥草(Lolium multiflorum,BAB40827.1)、荸薺(Eleocharis uniglumis,ADF46121.1)、蘋果(Malus domestica,ABS11067.1)和銀杏(Ginkgo biloba,CAA12930.1)等15個(gè)物種植物的逆轉(zhuǎn)錄酶序列之間具有較高的相似性,親緣關(guān)系較近。Ⅰ類包含來自大豆和菠菜的2條逆轉(zhuǎn)錄酶序列,除Ⅰ類外,F(xiàn)類~M類中均是來自2份AA染色體組野生種花生的逆轉(zhuǎn)錄酶序列。

    2.8 逆轉(zhuǎn)座子轉(zhuǎn)錄活性分析

    將65條AA染色體組野生種花生Ty3-gypsy類逆轉(zhuǎn)座子逆轉(zhuǎn)錄酶序列提交到NCBI與花生EST數(shù)據(jù)庫進(jìn)行比對, 以檢測AA染色體組野生種花生Ty3-gypsy類逆轉(zhuǎn)座子的轉(zhuǎn)錄活性。結(jié)果顯示,當(dāng)查詢覆蓋度都為97%時(shí),AdRT3-7、AdRT3-13、AdRT3-36與花生EST數(shù)據(jù)庫中GO261148.1之間的一致性分別為85.48%、86.26%、85.99%;當(dāng)查詢覆蓋度分別為96%和93%時(shí),AdRT4-7、AdRT4-28與花生EST數(shù)據(jù)庫中GO266033.1之間的一致性分別為88.22%和85.68%(表5)。說明GO261148.1和GO266033.1這2條序列為AA染色體組野生種花生Ty3-gypsy類逆轉(zhuǎn)座子的部分轉(zhuǎn)錄序列。

    3 討論與結(jié)論

    本研究首次PCR擴(kuò)增AA染色體組野生種花生中Ty3-gypsy類逆轉(zhuǎn)座子逆轉(zhuǎn)錄酶序列,結(jié)果從2份AA染色體組野生種花生材料中均擴(kuò)增出大小約430 bp的目的條帶,這與前人的研究結(jié)果[28-33]一致,最終也從2份AA染色體組野生種花生材料中分別獲得了32、33條逆轉(zhuǎn)錄酶序列。以上說明 Ty3-gypsy 類逆轉(zhuǎn)座子廣泛存在于本研究AA染色體組野生種花生材料基因組中,當(dāng)然也證實(shí)了采用簡并引物從本研究AA染色體組野生種花生材料中擴(kuò)增和克隆Ty3-gypsy類逆轉(zhuǎn)座子逆轉(zhuǎn)錄酶序列的策略是行之有效的[28]。

    所有逆轉(zhuǎn)錄酶序列長度在397~440 bp之間,存在缺失或插入突變;2份AA染色體組野生種花生逆轉(zhuǎn)錄酶序列的A+T所占比例范圍分別為56.48%~65.97%和56.48%~68.14%,(A+T)/(G+C)分別為1.3~1.94和1.3~2.13,均富含 A+T 堿基,A/T堿基含量的增加使序列呈現(xiàn)較高異質(zhì)性;2份AA染色體組野生種花生逆轉(zhuǎn)錄酶核苷酸序列間相似性范圍分別為63.7%~99.3%和62.8%~99.3%,呈現(xiàn)較低異質(zhì)性;翻譯后65條逆轉(zhuǎn)錄酶序列中有26條發(fā)生了無義突變,其中有12條來自A.duranensis(PI262133),有14條來自A.duranensis(PI219823),2份材料的無義突變發(fā)生率相當(dāng),無義突變導(dǎo)致產(chǎn)生較高的異質(zhì)性,無義突變也會使基因功能發(fā)生改變或喪失;2份AA染色體組野生種花生所有逆轉(zhuǎn)錄酶氨基酸序列間相似性范圍分別為10.9%~99.3%和31.5%~99.3%,呈現(xiàn)出高度的異質(zhì)性,相比核苷酸序列,氨基酸序列表現(xiàn)出更高的異質(zhì)性。B5F2FA2E-8CD4-41A0-93B1-4A7CB397A425

    9個(gè)家族所選代表序列的蛋白質(zhì)結(jié)構(gòu)總體類似,但也在螺旋結(jié)構(gòu)數(shù)、折疊結(jié)構(gòu)數(shù)、轉(zhuǎn)角數(shù)、氫鍵數(shù)、α-螺旋數(shù)、β-折疊數(shù)上存在著差別,如家族Ⅵ中AdRT4-5和Ⅸ中AdRT4-23的α-螺旋數(shù)、β- 折疊數(shù)和氫鍵數(shù)少于其他代表序列,這些差異可能會影響AA染色體組野生種花生Ty3-gypsy類逆轉(zhuǎn)座子逆轉(zhuǎn)錄酶的轉(zhuǎn)錄活性、轉(zhuǎn)座效率及拷貝數(shù),也是AA染色體組野生種花生Ty3-gypsy類逆轉(zhuǎn)座子存在較高異質(zhì)性和多態(tài)性的體現(xiàn)。AA染色體組野生種花生Ty3-gypsy類逆轉(zhuǎn)座子逆轉(zhuǎn)錄酶序列的主要保守基序是4種保守基序(motif 1~motif 4),motif 5~motif 10這6種保守基序在所克隆序列中出現(xiàn)的頻率較低,可能是這些序列在進(jìn)化過程中發(fā)生了突變,各序列間保守基序的差異反映了AA染色體組野生種花生Ty3-gypsy類逆轉(zhuǎn)座子逆轉(zhuǎn)錄酶序列具有較高異質(zhì)性和多態(tài)性。

    根據(jù)聚類分析,所有逆轉(zhuǎn)錄酶核苷酸序列被劃分為9個(gè)家族,其中構(gòu)成AA染色體組野生種花生Ty3-gypsy類逆轉(zhuǎn)座子的主要成分是家族Ⅰ,表明AA染色體組野生種花生Ty3-gypsy類逆轉(zhuǎn)座子逆轉(zhuǎn)錄酶序列具有較高的保守性與相似性。系統(tǒng)進(jìn)化樹顯示,E類中的5條AA染色體組野生種花生逆轉(zhuǎn)錄酶序列與綠豆、黃瓜、山桑、落葉松、牡丹、火龍果、李子、油棕、擬南芥、棗、白皮松、多花黑麥草、荸薺、蘋果和銀杏的逆轉(zhuǎn)錄酶序列之間具有較高的相似性,親緣關(guān)系較近,表明這些逆轉(zhuǎn)錄酶序列在進(jìn)化過程中有可能發(fā)生了橫向傳遞。A類和B類中包含大部分AA染色體組野生種花生Ty3-gypsy類逆轉(zhuǎn)座子逆轉(zhuǎn)錄酶序列,不包含來自其他物種植物的逆轉(zhuǎn)錄酶序列,C類只有AdRT3-17,D類也只有AdRT4-7,這4類可能較為特異,有可能只存在于AA染色體組野生種花生中。F類~M類中的逆轉(zhuǎn)錄酶序列均來自AA染色體組野生種花生,這些逆轉(zhuǎn)錄酶序列不僅與其他物種植物的同源序列遺傳距離較遠(yuǎn),也與本研究中AA染色體組野生種花生材料的逆轉(zhuǎn)錄酶序列遺傳距離較遠(yuǎn),表明這些逆轉(zhuǎn)錄酶序列在起源和進(jìn)化上可能較為古老,特異性比較強(qiáng),有可能是本研究AA染色體組野生種花生材料所特有。

    通過比對花生EST數(shù)據(jù)庫,發(fā)現(xiàn)3條來自A.duranensis(PI262133)的逆轉(zhuǎn)錄酶序列與GO261148.1一致性高,發(fā)現(xiàn)2條來自A.duranensis(PI219823)的逆轉(zhuǎn)錄酶序列與GO266033.1一致性高,表明AA染色體組野生種花生材料中存在可轉(zhuǎn)錄的Ty3-gypsy類逆轉(zhuǎn)座子。結(jié)合核苷酸序列的聚類分析,發(fā)現(xiàn)這2種來源的序列分別屬于家族Ⅴ和家族Ⅵ。

    參考文獻(xiàn):

    [1]熊發(fā)前,劉俊仙,賀梁瓊,等. 花生LTR和MITE轉(zhuǎn)座子及其分子標(biāo)記開發(fā)利用研究進(jìn)展[J]. 分子植物育種,2017,15(2):640-647.

    [2]Zhuang W,Chen H,Yang M,et al. The genome of cultivated peanut provides insight into legume karyotypes,polyploid evolution and crop domestication[J]. Nature Genetics,2019,51(5):865-876.

    [3]Kumar A,Pearce S R,McLean K,et al. The Ty1-copia group of retrotransposons in plants:genomic organisation,evolution,and use as molecular markers[J]. Genetica,1997,100(1/2/3):205-217.

    [4]Kumar A,Bennetzen J L. Plant retrotransposons[J]. Annual Review of Genetics,1999,33:479-532.

    [5]Feschotte C,Jiang N,Wessler S R. Plant transposable elements:where genetics meets genomics[J]. Nature Reviews Genetics,2002,3(5):329-341.

    [6]Bonchev G,Parisod C. Transposable elements and microevolutionary changes in natural populations[J]. Molecular Ecology Resources,2013,13(5):765-775.

    [7]Waugh R,McLean K,F(xiàn)lavell A J,et al. Genetic distribution of Bare-1-like retrotransposable elements in the barley genome revealed by sequence-specific amplification polymorphisms (S-SAP)[J]. Molecular & General Genetics,1997,253(6):687-694.

    [8]Kalendar R,Grob T,Regina M,et al. IRAP and REMAP:two new retrotransposon-based DNA fingerprinting techniques[J]. Theoretical and Applied Genetics,1999,98(5):704-711.

    [9]Patel M,Jung S,Moore K,et al. High-oleate peanut mutants result from a MITE insertion into the FAD2 gene[J]. Theoretical and Applied Genetics,2004,108(8):1492-1502.B5F2FA2E-8CD4-41A0-93B1-4A7CB397A425

    [10]Gowda M V C,Bhat R S,Motagi B N,et al. Association of high-frequency origin of late leaf spot resistant mutants with AhMITE1 transposition in peanut[J]. Plant Breeding,2010,129(5):567-569.

    [11]Gowda M V C,Bhat R S,Sujay V,et al. Characterization of AhMITE1 transposition and its association with the mutational and evolutionary origin of botanical types in peanut (Arachis spp.)[J]. Plant Systematics and Evolution,2011,291(3/4):153-158.

    [12]Shirasawa K,Hirakawa H,Tabata S,et al. Characterization of active miniature inverted-repeat transposable elements in the peanut genome[J]. Theoretical and Applied Genetics,2012,124(8):1429-1438.

    [13]Shirasawa K,Koilkonda P,Aoki K,et al. In silico polymorphism analysis for the development of simple sequence repeat and transposon markers and construction of linkage map in cultivated peanut[J]. BMC Plant Biology,2012,12:80.

    [14]Shirasawa K,Bertioli D J,Varshney R K,et al. Integrated consensus map of cultivated peanut and wild relatives reveals structures of the A and B genomes of Arachis and divergence of the legume genomes[J]. DNA Research,2013,20(2):173-184.

    [15]Mondal S,Hande P,Badigannavar A M.Identification of transposable element markers for a rust (Puccinia arachidis Speg.) resistance gene in cultivated peanut[J]. Journal of Phytopathology,2014,162(7/8):548-552.

    [16]Hake A A,Shirasawa K,Yadawad A,et al. Mapping of important taxonomic and productivity traits using genic and non-genic transposable element markers in peanut (Arachis hypogaea L.)[J]. PLoS One,2017,12(10):e0186113.

    [17]Gayathri M,Shirasawa K,Varshney R K,et al. Development of AhMITE1 markers through genome-wide analysis in peanut (Arachis hypogaea L.)[J]. BMC Research Notes,2018,11(1):10.

    [18]王 潔,李雙鈴,王 輝,等. 利用AhMITE1轉(zhuǎn)座子分子標(biāo)記鑒定花生F 1代雜種[J]. 花生學(xué)報(bào),2012,41(2):8-12.

    [19]王 輝,李雙鈴,任 艷,等. 利用AhMITE轉(zhuǎn)座子分子標(biāo)記研究花生栽培種及高世代材料的親緣關(guān)系[J]. 農(nóng)業(yè)生物技術(shù)學(xué)報(bào),2013,21(10):1176-1184.

    [20]尹 亮,任 艷,石延茂,等. 利用AhMITE1轉(zhuǎn)座子分子標(biāo)記鑒定栽培花生雜交F 1代種子真?zhèn)蝃J]. 山東農(nóng)業(yè)科學(xué),2015,47(12):1-5.

    [21]吳 琪,曹廣英,尹 亮,等. 利用AhMITE轉(zhuǎn)座子分子標(biāo)記鑒定花生雜交F 1代真假雜種[J]. 花生學(xué)報(bào),2017,46(3):48-53,60.

    [22]劉 婷,王傳堂,唐月異,等. 利用近紅外技術(shù)和轉(zhuǎn)座子標(biāo)記鑒定花生雜交F 1真雜種[J]. 分子植物育種,2017,15(9):3592-3598.

    [23]許夢琦,李雙鈴,任 艷,等. 花生作圖親本間分子標(biāo)記的多態(tài)性分析[J]. 湖北農(nóng)業(yè)科學(xué),2015,54(11):2763-2766.

    [24]Nielen S,Campos-Fonseca F,Leal-Bertioli S,et al. FIDEL-a retrovirus-like retrotransposon and its distinct evolutionary histories in the A-and B-genome components of cultivated peanut[J]. Chromosome Research,2010,18(2):227-246.B5F2FA2E-8CD4-41A0-93B1-4A7CB397A425

    [25]Nielen S,Vidigal B S,Leal-Bertioli S C M,et al. Matita,a new retroelement from peanut:characterization and evolutionary context in the light of the Arachis A-B genome divergence[J]. Molecular Genetics and Genomics,2012,287(1):21-38.

    [26]陽太億,劉俊仙,劉 菁,等. 四倍體野生種花生Ty1-copia類逆轉(zhuǎn)座子逆轉(zhuǎn)錄酶基因的克隆與分析[J]. 山東農(nóng)業(yè)科學(xué),2019,51(9):9-20.

    [27]熊發(fā)前,劉俊仙,劉 菁,等. 花生DNA的五種改良CTAB提取方法的比較分析及其應(yīng)用[J]. 分子植物育種,2019,17(7):2207-2216.

    [28]Kumekawa N,Ohtsubo E,Ohtsubo H.Identification and phylogenetic analysis of Gypsy-type retrotransposons in the plant kingdom[J]. Genes & Genetic Systems,1999,74(6):299-307.

    [29]王慶竹,李慧平,文曉鵬,等. 桂花LTR類反轉(zhuǎn)錄轉(zhuǎn)座子RT序列的克隆及分析[J]. 園藝學(xué)報(bào),2018,45(2):309-320.

    [30]侯小改,郭大龍,黃燕梅,等. 牡丹Ty3-gypsy類反轉(zhuǎn)錄轉(zhuǎn)座子反轉(zhuǎn)錄酶序列的克隆及分析[J]. 園藝學(xué)報(bào),2013,40(1):98-106.

    [31]彭 磊,吳 艷,劉小翠,等. 火龍果Ty3-gypsy類反轉(zhuǎn)錄轉(zhuǎn)座子反轉(zhuǎn)錄酶序列的克隆及分析[J]. 果樹學(xué)報(bào),2017,34(2):186-195.

    [32]張文波,陳 凌,李雪輝,等. 興安落葉松Ty3-gypsy類逆轉(zhuǎn)座子逆轉(zhuǎn)錄酶的多樣性分析[J]. 分子植物育種,2016,14(5):1098-1106.

    [33]白 楊,林曉飛,張文波.雜交構(gòu)樹Ty3-gypsy類逆轉(zhuǎn)座子逆轉(zhuǎn)錄酶序列的多樣性分析[J]. 分子植物育種,2018,16(22):7429-7437.B5F2FA2E-8CD4-41A0-93B1-4A7CB397A425

    猜你喜歡
    逆轉(zhuǎn)錄酶多樣性
    四倍體野生種花生Ty1-copia類逆轉(zhuǎn)座子逆轉(zhuǎn)錄酶基因的克隆與分析
    核苷類逆轉(zhuǎn)錄酶抑制劑(NRTIs)對大鼠卵母細(xì)胞功能的影響
    由古典戲曲看“代言體”在中國的前世今生
    戲劇之家(2016年22期)2016-11-30 15:13:39
    淺談新時(shí)期群文輔導(dǎo)工作的特征
    新時(shí)期群文輔導(dǎo)工作的特征
    海洋微生物次生代謝的生物合成機(jī)制
    科技資訊(2016年19期)2016-11-15 10:39:12
    乙型肝炎病毒逆轉(zhuǎn)錄酶基因突變的臨床意義
    舞蹈表演的表現(xiàn)形式多樣性研究
    人間(2016年27期)2016-11-11 16:27:23
    水磨地區(qū)蕨類植物多樣性調(diào)查分析
    端粒酶逆轉(zhuǎn)錄酶與轉(zhuǎn)錄激活蛋白-1在喉癌組織中的表達(dá)及相關(guān)性研究
    国产精品国产高清国产av| aaaaa片日本免费| 伊人久久大香线蕉亚洲五| 丝袜人妻中文字幕| 国产精品综合久久久久久久免费| 很黄的视频免费| 国产又色又爽无遮挡免费看| 久久精品国产99精品国产亚洲性色| 日本免费一区二区三区高清不卡| 欧美黑人巨大hd| 久久久久久九九精品二区国产 | 嫁个100分男人电影在线观看| 中文在线观看免费www的网站 | a级毛片a级免费在线| 国产精品免费视频内射| 一卡2卡三卡四卡精品乱码亚洲| 最近最新中文字幕大全免费视频| 免费在线观看视频国产中文字幕亚洲| 亚洲欧美激情综合另类| 国产精品免费一区二区三区在线| 欧美久久黑人一区二区| 十八禁人妻一区二区| 丁香欧美五月| 日韩有码中文字幕| 亚洲九九香蕉| 宅男免费午夜| 一边摸一边做爽爽视频免费| 精品午夜福利视频在线观看一区| 成人午夜高清在线视频| 久久精品夜夜夜夜夜久久蜜豆 | 黄色 视频免费看| 色综合欧美亚洲国产小说| 中亚洲国语对白在线视频| 久久精品国产亚洲av香蕉五月| 在线a可以看的网站| 人妻丰满熟妇av一区二区三区| 久久性视频一级片| 黑人巨大精品欧美一区二区mp4| 亚洲av熟女| 五月伊人婷婷丁香| 国内精品一区二区在线观看| 最新在线观看一区二区三区| 精品国产乱码久久久久久男人| 丰满的人妻完整版| 又粗又爽又猛毛片免费看| 国语自产精品视频在线第100页| 婷婷丁香在线五月| 99国产综合亚洲精品| 美女午夜性视频免费| 成人特级黄色片久久久久久久| 91老司机精品| 九色成人免费人妻av| 91字幕亚洲| 久久国产精品影院| 日日摸夜夜添夜夜添小说| 国产成人一区二区三区免费视频网站| 1024视频免费在线观看| 黄色女人牲交| 久久国产乱子伦精品免费另类| 国产精品亚洲一级av第二区| 精品久久久久久久人妻蜜臀av| 很黄的视频免费| 一边摸一边抽搐一进一小说| 五月玫瑰六月丁香| 亚洲 国产 在线| 国产av一区在线观看免费| 亚洲国产欧美一区二区综合| 中文字幕人妻丝袜一区二区| 少妇粗大呻吟视频| 亚洲自偷自拍图片 自拍| 国产成人一区二区三区免费视频网站| 50天的宝宝边吃奶边哭怎么回事| 禁无遮挡网站| 首页视频小说图片口味搜索| 1024香蕉在线观看| 999精品在线视频| 搡老熟女国产l中国老女人| 亚洲第一电影网av| 亚洲成人久久爱视频| 免费搜索国产男女视频| 日韩成人在线观看一区二区三区| 在线观看舔阴道视频| 丰满人妻熟妇乱又伦精品不卡| 久久久精品欧美日韩精品| 一进一出抽搐gif免费好疼| 亚洲国产欧美网| 狂野欧美激情性xxxx| 午夜福利在线观看吧| 久久精品国产清高在天天线| 人人妻,人人澡人人爽秒播| 黄色成人免费大全| 99久久99久久久精品蜜桃| 三级毛片av免费| 午夜免费观看网址| 国产精品亚洲一级av第二区| 在线看三级毛片| 国产成人精品久久二区二区91| 99在线人妻在线中文字幕| avwww免费| 精品国产亚洲在线| 两性午夜刺激爽爽歪歪视频在线观看 | 午夜精品在线福利| 国产成人精品无人区| 免费av毛片视频| 婷婷亚洲欧美| 一区二区三区国产精品乱码| 国产久久久一区二区三区| 国产69精品久久久久777片 | 国产91精品成人一区二区三区| 中文字幕最新亚洲高清| 99国产精品一区二区蜜桃av| 久久精品aⅴ一区二区三区四区| 91在线观看av| 欧美大码av| 精品日产1卡2卡| 亚洲人成伊人成综合网2020| 成年人黄色毛片网站| 国产精品久久久久久人妻精品电影| 一二三四社区在线视频社区8| 夜夜看夜夜爽夜夜摸| 国产成人精品久久二区二区免费| 黄片大片在线免费观看| 丰满的人妻完整版| 久久久国产欧美日韩av| 午夜两性在线视频| 国产亚洲精品一区二区www| 脱女人内裤的视频| 日本在线视频免费播放| 欧美中文日本在线观看视频| 国产熟女午夜一区二区三区| 91成年电影在线观看| 丰满的人妻完整版| 一二三四在线观看免费中文在| 久久热在线av| 欧美中文日本在线观看视频| 午夜精品久久久久久毛片777| 亚洲18禁久久av| 中文字幕人成人乱码亚洲影| 我的老师免费观看完整版| 在线a可以看的网站| 久久久久免费精品人妻一区二区| 国产久久久一区二区三区| 国产亚洲欧美在线一区二区| 久久精品夜夜夜夜夜久久蜜豆 | 又黄又粗又硬又大视频| 亚洲av成人精品一区久久| 亚洲男人的天堂狠狠| 久久久久久亚洲精品国产蜜桃av| 亚洲精品国产精品久久久不卡| 女同久久另类99精品国产91| 欧美色视频一区免费| 亚洲五月天丁香| 我要搜黄色片| 黄色视频不卡| 国产爱豆传媒在线观看 | 成人特级黄色片久久久久久久| 又黄又粗又硬又大视频| 99久久99久久久精品蜜桃| 亚洲天堂国产精品一区在线| 欧美日韩亚洲国产一区二区在线观看| 别揉我奶头~嗯~啊~动态视频| 91麻豆精品激情在线观看国产| 精品人妻1区二区| 色综合婷婷激情| 国产不卡一卡二| av片东京热男人的天堂| 日本a在线网址| 天天添夜夜摸| 天堂动漫精品| 18禁国产床啪视频网站| 国产精品综合久久久久久久免费| 999久久久精品免费观看国产| 国产黄片美女视频| 亚洲色图av天堂| 欧美激情久久久久久爽电影| 国产亚洲精品久久久久久毛片| 天天一区二区日本电影三级| 18禁裸乳无遮挡免费网站照片| 国产亚洲欧美在线一区二区| 久久久久久久久久黄片| 男女做爰动态图高潮gif福利片| 在线永久观看黄色视频| 丁香欧美五月| 神马国产精品三级电影在线观看 | 免费在线观看视频国产中文字幕亚洲| 淫妇啪啪啪对白视频| 午夜精品一区二区三区免费看| 两个人看的免费小视频| 色av中文字幕| 亚洲色图 男人天堂 中文字幕| 免费在线观看黄色视频的| 久久亚洲真实| 日本免费a在线| 亚洲国产精品999在线| 亚洲乱码一区二区免费版| 黑人巨大精品欧美一区二区mp4| 无限看片的www在线观看| 久久精品国产综合久久久| 99久久99久久久精品蜜桃| www国产在线视频色| tocl精华| 好男人在线观看高清免费视频| 国内毛片毛片毛片毛片毛片| 亚洲中文字幕一区二区三区有码在线看 | 欧美丝袜亚洲另类 | 中文字幕熟女人妻在线| 久久国产乱子伦精品免费另类| 久久精品人妻少妇| 国产激情欧美一区二区| 成年人黄色毛片网站| 日韩欧美国产一区二区入口| 久99久视频精品免费| 久久人妻福利社区极品人妻图片| 国产三级黄色录像| 亚洲一区二区三区色噜噜| www日本黄色视频网| 欧美 亚洲 国产 日韩一| 免费观看精品视频网站| 国产野战对白在线观看| 亚洲精品一卡2卡三卡4卡5卡| 久久国产精品影院| 亚洲成人免费电影在线观看| 国产一区二区三区视频了| 三级国产精品欧美在线观看 | 禁无遮挡网站| 黄色视频,在线免费观看| 99国产精品一区二区蜜桃av| 亚洲一码二码三码区别大吗| 日韩精品免费视频一区二区三区| 欧美av亚洲av综合av国产av| 18禁黄网站禁片免费观看直播| 每晚都被弄得嗷嗷叫到高潮| 午夜久久久久精精品| 亚洲精品国产一区二区精华液| 男插女下体视频免费在线播放| 亚洲精品久久国产高清桃花| 啪啪无遮挡十八禁网站| 丰满的人妻完整版| 久久这里只有精品中国| 亚洲成人国产一区在线观看| 无遮挡黄片免费观看| 国产精品综合久久久久久久免费| 亚洲中文字幕一区二区三区有码在线看 | 99国产极品粉嫩在线观看| 亚洲国产高清在线一区二区三| 国产探花在线观看一区二区| 伦理电影免费视频| 日本 欧美在线| 久久久久九九精品影院| 亚洲国产精品成人综合色| 亚洲,欧美精品.| 一级黄色大片毛片| 国产精品爽爽va在线观看网站| 国产不卡一卡二| 99国产精品99久久久久| 91国产中文字幕| 美女午夜性视频免费| 在线观看舔阴道视频| 国产午夜精品论理片| 一本大道久久a久久精品| 欧美日韩国产亚洲二区| 婷婷亚洲欧美| 老司机靠b影院| 欧美又色又爽又黄视频| 国产爱豆传媒在线观看 | 欧美大码av| xxxwww97欧美| 精品日产1卡2卡| 欧美大码av| 99国产精品一区二区蜜桃av| 久久99热这里只有精品18| 中文亚洲av片在线观看爽| 给我免费播放毛片高清在线观看| 在线观看日韩欧美| 99国产极品粉嫩在线观看| 国产爱豆传媒在线观看 | 久久久久久人人人人人| 成人永久免费在线观看视频| 老司机深夜福利视频在线观看| 欧美精品啪啪一区二区三区| 激情在线观看视频在线高清| √禁漫天堂资源中文www| 又爽又黄无遮挡网站| 日本一区二区免费在线视频| 热99re8久久精品国产| 99久久国产精品久久久| 亚洲人成伊人成综合网2020| 精品久久蜜臀av无| 亚洲在线自拍视频| netflix在线观看网站| 色综合站精品国产| 久久人妻福利社区极品人妻图片| 51午夜福利影视在线观看| 日本精品一区二区三区蜜桃| 国产精品日韩av在线免费观看| 中文字幕熟女人妻在线| 一夜夜www| 男插女下体视频免费在线播放| 黄色片一级片一级黄色片| 亚洲色图 男人天堂 中文字幕| 国产av又大| 国产区一区二久久| 日韩欧美在线二视频| 欧美在线一区亚洲| 久久中文字幕人妻熟女| 美女高潮喷水抽搐中文字幕| 宅男免费午夜| 一个人免费在线观看电影 | 老司机靠b影院| 国产伦在线观看视频一区| 啦啦啦观看免费观看视频高清| 国产99久久九九免费精品| 国产激情欧美一区二区| 99久久无色码亚洲精品果冻| 午夜福利在线在线| 99久久久亚洲精品蜜臀av| 国产午夜精品论理片| 久久国产精品影院| 免费看a级黄色片| 欧美不卡视频在线免费观看 | 在线看三级毛片| 久久久久九九精品影院| 久久久久国内视频| 丝袜人妻中文字幕| 久久中文字幕一级| 久久久精品欧美日韩精品| 黄色视频,在线免费观看| 精品国产美女av久久久久小说| 五月伊人婷婷丁香| 日本五十路高清| 最新美女视频免费是黄的| 亚洲人与动物交配视频| 禁无遮挡网站| 欧美一级毛片孕妇| cao死你这个sao货| 欧美日韩亚洲综合一区二区三区_| 日韩欧美精品v在线| 久久午夜综合久久蜜桃| 日日爽夜夜爽网站| 欧美激情久久久久久爽电影| 欧美+亚洲+日韩+国产| 午夜免费成人在线视频| av免费在线观看网站| 久久精品国产亚洲av高清一级| 午夜免费激情av| 嫁个100分男人电影在线观看| av片东京热男人的天堂| 制服诱惑二区| √禁漫天堂资源中文www| 欧美绝顶高潮抽搐喷水| 精品欧美国产一区二区三| 精品电影一区二区在线| 正在播放国产对白刺激| 亚洲中文av在线| 久久久久久久精品吃奶| 2021天堂中文幕一二区在线观| 免费在线观看完整版高清| 99国产综合亚洲精品| 青草久久国产| 国产精华一区二区三区| 成人国语在线视频| 成人18禁高潮啪啪吃奶动态图| 一级毛片女人18水好多| 国产成人影院久久av| 色在线成人网| 国产高清视频在线观看网站| 99国产综合亚洲精品| 久久久水蜜桃国产精品网| 日本精品一区二区三区蜜桃| 日本五十路高清| 日韩大尺度精品在线看网址| 国产aⅴ精品一区二区三区波| 成人三级做爰电影| 亚洲aⅴ乱码一区二区在线播放 | 成人高潮视频无遮挡免费网站| 日本五十路高清| 特级一级黄色大片| 好男人电影高清在线观看| 国产又色又爽无遮挡免费看| 国产91精品成人一区二区三区| АⅤ资源中文在线天堂| 国产精品久久久久久亚洲av鲁大| 精品高清国产在线一区| 麻豆久久精品国产亚洲av| 日韩精品免费视频一区二区三区| 大型av网站在线播放| 99久久久亚洲精品蜜臀av| 久久婷婷成人综合色麻豆| 亚洲精品久久成人aⅴ小说| videosex国产| 久久久久久久精品吃奶| 法律面前人人平等表现在哪些方面| 日本熟妇午夜| 黑人巨大精品欧美一区二区mp4| 午夜a级毛片| 中国美女看黄片| 国产人伦9x9x在线观看| 啦啦啦韩国在线观看视频| 麻豆久久精品国产亚洲av| 久久久水蜜桃国产精品网| 久9热在线精品视频| 最新在线观看一区二区三区| 桃红色精品国产亚洲av| 老司机福利观看| 日本一区二区免费在线视频| 日日摸夜夜添夜夜添小说| 国产精品综合久久久久久久免费| 波多野结衣高清作品| 一级作爱视频免费观看| 午夜精品一区二区三区免费看| 国产av麻豆久久久久久久| 中文亚洲av片在线观看爽| 999精品在线视频| 久久香蕉精品热| 欧美av亚洲av综合av国产av| 丝袜美腿诱惑在线| 叶爱在线成人免费视频播放| 18禁裸乳无遮挡免费网站照片| 特大巨黑吊av在线直播| 亚洲午夜理论影院| 久久香蕉精品热| 大型黄色视频在线免费观看| 长腿黑丝高跟| 最近在线观看免费完整版| 亚洲 欧美一区二区三区| 久久天躁狠狠躁夜夜2o2o| 免费看美女性在线毛片视频| 国产成人欧美在线观看| 日本熟妇午夜| 黄色a级毛片大全视频| 日本黄色视频三级网站网址| 丰满人妻一区二区三区视频av | 精品国产亚洲在线| 久久精品国产清高在天天线| 制服丝袜大香蕉在线| 免费在线观看亚洲国产| 日本一本二区三区精品| 国产午夜福利久久久久久| 亚洲av美国av| 欧美乱妇无乱码| 欧美一级毛片孕妇| 色综合婷婷激情| 波多野结衣高清无吗| 每晚都被弄得嗷嗷叫到高潮| 亚洲欧美精品综合一区二区三区| 亚洲美女视频黄频| av国产免费在线观看| 亚洲人成77777在线视频| 亚洲人与动物交配视频| 黑人操中国人逼视频| 好男人电影高清在线观看| 久久久久久大精品| 免费看十八禁软件| 国产69精品久久久久777片 | av片东京热男人的天堂| www国产在线视频色| 99国产精品99久久久久| 国产av在哪里看| 色播亚洲综合网| 黄色丝袜av网址大全| 国产精品国产高清国产av| 欧美在线黄色| 男女之事视频高清在线观看| 91麻豆av在线| 国产精品98久久久久久宅男小说| 18禁黄网站禁片免费观看直播| 香蕉丝袜av| 中文字幕熟女人妻在线| 高清毛片免费观看视频网站| 国产精品亚洲av一区麻豆| 在线国产一区二区在线| 亚洲成人精品中文字幕电影| 日韩av在线大香蕉| 热99re8久久精品国产| 精品国产乱子伦一区二区三区| 欧美激情久久久久久爽电影| 国产真人三级小视频在线观看| 亚洲精品久久国产高清桃花| 亚洲av电影不卡..在线观看| 日本黄色视频三级网站网址| 国产精品日韩av在线免费观看| 麻豆国产av国片精品| 黑人欧美特级aaaaaa片| 999久久久精品免费观看国产| 久久久久久国产a免费观看| 久久久久性生活片| 亚洲最大成人中文| 免费在线观看完整版高清| 色综合亚洲欧美另类图片| 午夜福利在线在线| 欧美日韩瑟瑟在线播放| 黄色丝袜av网址大全| 老司机福利观看| 少妇的丰满在线观看| 别揉我奶头~嗯~啊~动态视频| 三级毛片av免费| 国产v大片淫在线免费观看| 脱女人内裤的视频| 国产片内射在线| 国产片内射在线| 麻豆久久精品国产亚洲av| a级毛片a级免费在线| 女人被狂操c到高潮| 一本大道久久a久久精品| 一a级毛片在线观看| 人成视频在线观看免费观看| www.熟女人妻精品国产| 亚洲一区二区三区色噜噜| 亚洲精品美女久久久久99蜜臀| 狠狠狠狠99中文字幕| 国内久久婷婷六月综合欲色啪| 国产精品一区二区三区四区免费观看 | 夜夜躁狠狠躁天天躁| 免费在线观看影片大全网站| 俄罗斯特黄特色一大片| 免费在线观看成人毛片| 免费一级毛片在线播放高清视频| 精品电影一区二区在线| av视频在线观看入口| 91老司机精品| 成人三级黄色视频| 免费观看人在逋| 男女之事视频高清在线观看| 亚洲男人的天堂狠狠| 老汉色∧v一级毛片| 搡老妇女老女人老熟妇| 国产成人av激情在线播放| 又粗又爽又猛毛片免费看| 国产精品自产拍在线观看55亚洲| 久久久久性生活片| 一级黄色大片毛片| 91成年电影在线观看| 成人三级做爰电影| 老司机午夜福利在线观看视频| 国产精品亚洲一级av第二区| 国产一区二区在线av高清观看| 欧美zozozo另类| 国产成人av激情在线播放| 午夜免费成人在线视频| a级毛片在线看网站| xxx96com| 国产成人系列免费观看| 国产成人影院久久av| 日日爽夜夜爽网站| 女人高潮潮喷娇喘18禁视频| 亚洲av电影在线进入| 欧美日韩亚洲国产一区二区在线观看| 一区二区三区激情视频| 久久99热这里只有精品18| 久久精品国产亚洲av高清一级| 成人一区二区视频在线观看| 日韩中文字幕欧美一区二区| 久久久国产欧美日韩av| 亚洲熟妇中文字幕五十中出| 免费观看人在逋| 久久午夜综合久久蜜桃| 亚洲成av人片在线播放无| 两个人免费观看高清视频| 久9热在线精品视频| 欧美激情久久久久久爽电影| 麻豆国产97在线/欧美 | 午夜视频精品福利| 成在线人永久免费视频| 欧美日韩一级在线毛片| 国产亚洲精品一区二区www| 男人舔奶头视频| 俺也久久电影网| 日韩大尺度精品在线看网址| 亚洲欧美日韩无卡精品| 99精品久久久久人妻精品| 91大片在线观看| 91字幕亚洲| 久久欧美精品欧美久久欧美| 亚洲精品美女久久av网站| 欧美人与性动交α欧美精品济南到| 一卡2卡三卡四卡精品乱码亚洲| 亚洲精品在线美女| 精华霜和精华液先用哪个| 嫁个100分男人电影在线观看| 成人一区二区视频在线观看| 999久久久国产精品视频| 草草在线视频免费看| 久久这里只有精品中国| 色精品久久人妻99蜜桃| 热99re8久久精品国产| 夜夜看夜夜爽夜夜摸| 欧美黄色片欧美黄色片| 亚洲欧美精品综合久久99| 特级一级黄色大片| 免费一级毛片在线播放高清视频| 黄色女人牲交| 欧美色视频一区免费| 国产亚洲欧美98| 丁香欧美五月| 首页视频小说图片口味搜索| 亚洲欧美激情综合另类| 久久国产乱子伦精品免费另类| 又爽又黄无遮挡网站| 成人三级黄色视频| 国产亚洲欧美98| 国产野战对白在线观看| 色av中文字幕| 他把我摸到了高潮在线观看| 50天的宝宝边吃奶边哭怎么回事| 亚洲熟妇熟女久久| 99久久99久久久精品蜜桃| 亚洲av第一区精品v没综合| 老熟妇乱子伦视频在线观看| 99riav亚洲国产免费| 亚洲中文字幕一区二区三区有码在线看 | 亚洲精品中文字幕一二三四区| 熟女少妇亚洲综合色aaa.| 一进一出抽搐gif免费好疼| 青草久久国产| 婷婷亚洲欧美|