亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        麝鼠香腺Elovl1基因生物信息學分析及表達研究

        2022-06-01 03:59:46李天峰李嘉昊白素英
        中國畜牧獸醫(yī) 2022年4期
        關鍵詞:脂肪酸測序物種

        雷 行,李天峰,李嘉昊,滿 欣,白素英,2

        (1.東北林業(yè)大學野生動物與自然保護地學院,哈爾濱 150040;2.國家林業(yè)和草原局野生動植物檢測中心,哈爾濱 150040)

        麝鼠(Ondatrazibethicus)是一種資源分布廣、適應性強、繁殖快的倉鼠科鼠類,其成年雄性麝鼠香腺在繁殖期會分泌麝鼠香,與麝香腺非分泌期也保有泌香狀態(tài)不同[1]。麝鼠香腺發(fā)育及泌香呈現明顯的周期性變化,分為香腺發(fā)育期、泌香盛期、泌香持續(xù)期和萎縮期。2~3月為麝鼠香腺發(fā)育期,4月為發(fā)育盛期進入泌香期,5~6月為泌香盛期,7~8月為泌香維持期,腺體在9月開始萎縮,并在10月萎縮到線狀或豆狀大小,到11月結締組織逐漸取代了腺體[2]。對麝鼠香化學成分進行氣相色譜-質譜分析顯示,其含有烷類、醇類、大環(huán)酮類、甾體類化合物及脂肪酸類化合物,主要成分為大環(huán)酮類和脂肪酸類化合物[3-4],且脂肪酸類化合物中長鏈脂肪酸(long chain fatty acids,LCFA)及其化合物占比偏高,主要集中在16C~24C之間。

        生物體內LCFA的合成是一個復雜且受多方面調控的過程,主要經過脂肪酸從頭合成形成不超過16C的飽和脂肪酸,再經過脂肪酸碳鏈的延長和不飽和鍵的形成,其中脂肪酸碳鏈的延伸主要發(fā)生于內質網,以脂酰CoA(fatty acyl-CoA)作為起點合成大于16C的LCFA[5-6]。目前已在哺乳動物中發(fā)現7個超長鏈脂肪酸延伸酶(elongase of verylong chain fatty acids,Elovl)家族成員,分別為Elovl1~Elovl7,其催化構成長鏈脂肪酸延伸循環(huán)的4個反應中的第1個反應,即限速反應。Elovl1延長飽和C18∶0~C26∶0及單不飽和C20∶1、C22∶1、C24∶1的脂酰CoA,是主要的脂肪酸延伸酶[7],廣泛存在于哺乳動物各個組織器官中。Elovl1為提供細胞自身LCFA需求起重要作用,分別對細胞中Elovl1進行敲除或抑制處理,檢測到LCFA含量明顯下調;反之對Elovl1進行過表達處理,檢測到LCFA含量明顯上調[8-9]。因此,Elovl家族成員可能與麝鼠香腺中的LCFA密不可分。鑒于此,本試驗對Elovl1基因生物信息學及轉錄水平展開分析,以期為深入了解麝鼠泌香的規(guī)律及LCFA相關合成通路在麝鼠香腺泌香周期中起到的調控作用提供參考。

        1 材料與方法

        1.1 試驗動物

        2019年2月末,在黑龍江省哈爾濱市五常麝鼠養(yǎng)殖場購入體況相近的1歲麝鼠40只,分別在2、3、4、5、6、7、8、9月末隨機挑選3只麝鼠,適應性飼養(yǎng)2 d后取其位于肌膜上腹部兩側的香腺,于液氮中速凍后放入凍存管中,-80 ℃保存?zhèn)溆谩?/p>

        1.2 主要試劑及儀器

        TRIzol試劑、氯仿、無水乙醇均購自哈爾濱欣科瑞經貿有限公司;總RNA提取試劑盒、RNA酶抑制劑均購自Thermo Fisher Scientific公司;反轉錄試劑盒購自哈爾濱擎科生物有限公司;熒光染料SYBR Green購自湖南艾科瑞生物工程有限公司;高保真Mix酶購自南京諾唯贊生物科技股份有限公司;瓊脂糖凝膠回收試劑盒(D2500 Gel Extraction Kit)購自Omega公司。

        PCR擴增儀購自ABI公司;NanoPhotometer-N50紫外分光光度計購自Implen公司;JY300C通用型電泳儀購自北京君意東方電泳設備有限公司;熒光定量PCR儀購自Bio-Rad公司。

        1.3 RNA提取與反轉錄

        將麝鼠香腺組織放入研缽中加入液氮快速研磨,研磨充分后使用TRIzol法提取其總RNA,使用超微量分光光度計測定RNA的濃度與純度,加入RNA酶抑制劑保存RNA。根據反轉錄試劑盒說明書對麝鼠香腺的總RNA進行反轉錄,-80 ℃保存?zhèn)溆谩?/p>

        1.4 PCR擴增及測序

        麝鼠Elovl1基因序列根據轉錄組測序獲得,由深圳華大基因公司完成,并進行PCR擴增測序檢驗。使用Primer Premier 6.0軟件設計引物,引物序列為:片段Ⅰ,F:5′-GCCCAGCAGATGAGG-AAGT-3′,R:5′-GCCTTGACCTTGGTGATAGC-3′;片段Ⅱ,F:5′-CCCAGCAGATGAGGAAGTG-3′,R:5′-TTTATCATGGCATGGAAAGAGC-3′。

        PCR反應體系50 μL:2×Phanta Max Buffer 25 μL,dNTP Mix 1 μL,10 μmol/L上、下游引物各2 μL,Phanta Max Super-Fidelity DNA Polymerase 1 μL,模板cDNA 5 μL,ddH2O補足體系。PCR反應程序:95 ℃預變性30 s;95 ℃變性15 s,56 ℃退火15 s,72 ℃延伸30 s,共30個循環(huán);72 ℃延伸5 min。將PCR產物進行1.0%瓊脂糖凝膠電泳,85 V恒壓電泳50~60 min。鑒定正確的PCR產物使用瓊脂糖凝膠回收試劑盒(D2500 Gel Extraction Kit)回收純化,純化后產物送哈爾濱擎科生物科技有限公司和吉林省庫美生物科技有限公司測序。

        1.5 系統(tǒng)進化樹構建

        將測序結果比對麝鼠全基因組庫獲得Elovl1基因CDS區(qū)序列,代碼設為Ondz X1和Ondz X2。從GenBank數據庫中下載嚙齒目不同物種的Elovl1基因及其多轉錄本CDS序列(表1),以家犬Elovl1基因作為外群,使用MegaX軟件進行多序列比對及系統(tǒng)進化樹構建。

        表1 各物種Elovl1基因CDS序列信息Table 1 CDS sequence information of Elovl1 gene in each species

        續(xù)表

        1.6 生物信息學分析

        使用ProtParam(https:∥web.expasy.org/protparam/)和ProtScale(https:∥web.expasy.org/protscale/)在線軟件分析蛋白理化性質和疏水性;使用TMHMM Server(https:∥services.healthtech.dtu.dk/service.php?TMHMM-2.0)預測蛋白跨膜結構;使用SignalP-5.0(http:∥www.cbs.dtu.dk/services/SignalP-5.0/)預測蛋白信號肽;使用SOPMA(https:∥npsa-prabi.ibcp.fr/cgi-bin/npsa_automat.pl?page=npsa_sopma.html)和Phyre2(http:∥www.sbg.bio.ic.ac.uk/phyre2/html/page.cgi?id=index)在線軟件預測蛋白二級結構和三級結構;使用Hum-mPLoc 3.0(http:∥www.csbio.sjtu.edu.cn/bioinf/Hum-mPLoc3/)進行蛋白亞細胞定位。

        1.7 Elovl1基因轉錄水平檢測

        使用實時熒光定量PCR方法檢測麝鼠香腺Elovl1基因的mRNA表達水平,利用Primer Premier 6.0軟件分別設計麝鼠不同CDS區(qū)的Elovl1基因通用引物和X2特異引物,引物序列為:通用引物,F:5′-CTCACTGCACATCAGCCAA-TACTACT-3′,R:5′-TGGAGAAGAGCACGAAG-AAGAAGGTA-3′;X2特異引物,F:5′-GACTGA-AGGTGAGCACTGGGAAG-3′,R:5′-TGATAC-AAGTTCACAACAGCCTCCAT-3′;內參基因β-actin引物[10],F:5′-TTGCTGATCCACATCTG-CT-3′,R:5′-GACAGGATGCAGAAGGAGAT-3′。 以Elovl1基因通用引物同時擴增X1+X2相同序列,以Elovl1基因X2引物特異性擴增X2序列。

        PCR反應體系20 μL:SYBR Green 10 μL,模板cDNA 2 μL,10 μmol/L上、下游引物各1 μL,ddH2O補足體系。PCR反應程序:95 ℃ 10 min;95 ℃ 10 s,60 ℃ 30 s,共40個循環(huán),退火時收集熒光信號。采用2-ΔΔCt法分析Elovl1基因在不同月份麝鼠香腺中的相對表達量。

        2 結 果

        2.1 麝鼠Elovl1基因序列分析

        經轉錄組測序,并經PCR擴增及測序驗證,麝鼠Elovl1基因2個轉錄本X1和X2的CDS區(qū)分別為942和840 bp,X2的5′-端缺失102 bp(圖1),其余序列完全相同。

        圖1 麝鼠Elovl1基因轉錄本X1和X2 CDS序列差異部分Fig.1 Differences in the CDS sequence of transcripts X1 and X2 of Elovl1 gene in muskrat

        2.2 Elovl1基因系統(tǒng)進化樹構建

        系統(tǒng)進化樹結果顯示,倉鼠科與鼠科聚在一起,二者與鼴形鼠科關系較近;八齒鼠科、毛絲鼠科、豚鼠科、濱鼠科親緣關系較近聚為一支,且這些科間分支的置信值普遍較高;松鼠科、河貍科則與其他科間關系較遠(圖2)。同屬物種間,親緣關系較近各自聚為一支,倉鼠科中,麝鼠與水鼠平、緩行田鼠、草原田鼠親緣關系最近;同一物種Elovl1基因的不同轉錄本間相對保守,均按物種聚在一起且置信值普遍較高,但旱獺例外,轉錄本間已出現進化分歧,其X3與黃腹旱獺聚在一起。

        麝鼠用粗線記號標出Muskrat is marked with a bold line圖2 Elovl1基因系統(tǒng)進化樹Fig.2 Phylogenetic tree of Elovl1 gene

        2.3 麝鼠Elovl1蛋白理化性質及親/疏水性預測

        Elovl1基因2個轉錄本X1和X2編碼蛋白分別由313和279個氨基酸組成,X2比X1缺失了N-端34個氨基酸,其他序列相同。利用在線軟件對麝鼠Elovl1基因編碼蛋白X1和X2的理化性質進行分析,結果顯示,蛋白分子質量分別為36.31和32.69 ku,分子式分別為C1732H2587N415O413S16和 C1558H2327N371O372S15,等電點分別為9.79和9.72,半衰期均為30 h,脂肪族氨基酸指數分別為102.43和105.09,不穩(wěn)定指數分別為34.50和31.12,推測該蛋白可能為穩(wěn)定堿性蛋白。蛋白包含20種氨基酸,帶負電荷氨基酸殘基(Asp+Glu)分別有10和9個,帶正電荷氨基酸殘基(Arg+Lys)分別有25和22個(表2),推測Elovl1蛋白整體帶正電。

        表2 麝鼠Elovl1蛋白X1和X2的氨基酸組成Table 2 Amino acid composition of X1 and X2 of Elovl1 protein in muskrat

        親/疏水性預測發(fā)現,麝鼠Elovl1蛋白X1、X2親水性最強的分別是第290和256位氨基酸處的蘇氨酸(-2.256),疏水性最強的分別是第107和73位氨基酸處的絲氨酸(2.856),肽鏈上的疏水性氨基酸多于親水性氨基酸,提示整條肽鏈表現為疏水性(圖3)。

        A,X1蛋白;B,X2蛋白。下同A,X1 protein;B,X2 protein.The same as below圖3 麝鼠Elovl1蛋白親/疏水性預測Fig.3 Hydrophilicity and hydrophobicity prediction of Elovl1 protein in muskrat

        2.4 麝鼠Elovl1蛋白信號肽分析

        利用SignalP-5.0軟件對麝鼠Elovl1蛋白X1和X2中信號肽進行分析,結果顯示,C、S、Y值均未達到信號肽的要求,信號肽存在可能性分別為0.0114和0.0005,推斷該蛋白不存在信號肽,為非分泌蛋白。

        2.5 麝鼠Elovl1蛋白二級結構和三級結構

        SOMPA在線軟件預測結果顯示,麝鼠Elovl1蛋白X1和X2的二級結構均由α-螺旋(46.45%和50.54%)、β-轉角(3.83%和2.51%)、延伸鏈(15.97%和16.85%)及無規(guī)則卷曲(36.74%和30.11%)構成(圖4),且Elovl1蛋白主要由α-螺旋和無規(guī)則卷曲構成。Phyre2在線軟件預測結果顯示,在Normol模式(同源建模法)下因缺少相關物種的該蛋白模型,均以人Elovl7為模板,預測的三級結構完全相同,由7個螺旋結構組成。遂采用Intensive模式(綜合同源建模法和從頭預測法)構建,結果顯示,Elovl1蛋白X1和X2的三級結構與二級結構分析結果基本一致,主要由7個螺旋結構組成(圖5),二者主要結構無明顯區(qū)別。

        e,延伸鏈;h,α-螺旋;c,無規(guī)則卷曲;t,β-轉角e,Extended chain;h,Alpha helix;c,Random coil;t,Beta turn圖4 麝鼠Elovl1蛋白二級結構預測Fig.4 Secondary structure prediction of Elovl1 protein in muskrat

        圖5 麝鼠Elovl1蛋白三級結構預測Fig.5 Tertiary structure prediction of Elovl1 protein in muskrat

        2.6 麝鼠Elovl1蛋白的跨膜結構與亞細胞定位

        利用TMHMM Server v 2.0分析發(fā)現,麝鼠Elovl1蛋白X1和X2跨膜螺旋數量均為7個(圖6),與三級結構預測一致。運用Hum-mPLoc 3.0預測麝鼠Elovl1蛋白X1和X2的亞細胞定位發(fā)現,該蛋白分布于內質網。

        圖6 Elovl1蛋白跨膜結構預測Fig.6 Transmembrane structure prediction of Elovl1 protein in muskrat

        2.7 麝鼠Elovl1基因轉錄水平檢測

        實時熒光定量PCR結果表明,Elovl1基因2個轉錄本在2~9月末的麝鼠香腺中均有表達,分別以X1+X2和X2的9月末表達量作為對照,結果顯示,X1+X2表達變化規(guī)律與X2基本一致,即2~9月末Elovl1基因在麝鼠香腺中的表達量整體均呈先上調后下調的趨勢(圖7)。由圖7可知,X1+X2在2、3、4、5、6、7、8月末的麝鼠香腺中的相對表達量分別為9月末的51.77、101.25、98.42、24.82、24.11、14.64和1.67倍,其中2~7月末的相對表達量均極顯著高于9月末(P<0.01);X2在2、3、4、5、6、7、8月末的麝鼠香腺中的相對表達量分別為9月末的13.66、265.39、157.61、106.91、106.02、83.34和3.83倍,其中3~7月末的相對表達量均極顯著高于9月末(P<0.01),2、8月末的相對表達量均顯著高于9月末(P<0.05)。

        ①A,X1+X2總擴增;B,X2擴增。②與9月末相比,*,差異顯著(P<0.05);**,差異極顯著(P<0.01)①A,X1+X2 total amplification;B,X2 amplification.②Compared with the end of September,*,Significant difference (P<0.05);**,Extremely significant difference (P<0.01)圖7 不同月份麝鼠香腺Elovl1基因相對表達量Fig.7 The relative expression of Elovl1 gene in muskrat scent gland in different months

        以X2各月表達量為對照,X1+X2與X2各月表達量倍數關系見圖8。由圖8可知,2和9月末的表達量最高,達131.33和134.01,各月平均倍數為74.6。說明X1+X2的表達量遠遠高于X2。

        圖8 X1+X2與X2不同月份相對表達量倍數關系Fig.8 The relationship of relative expression between X1+X2 and X2 in different months

        3 討 論

        3.1 麝鼠Elovl1基因的可變剪接

        Elovl1基因在麝鼠香腺中存在2個轉錄本,2個不同長度的CDS序列使用不同的5′-UTR,且轉錄本X2的CDS序列相比X1其5′-端缺失102 bp。在動物體中絕大多數基因存在可變剪切,通過產生不同的轉錄本極大程度地豐富了體內蛋白質的數量,且讓同一蛋白各亞型表現出不同的生物學特性[11-12]。如Bcl-x基因,具有Bcl-xl和Bcl-xs 2種轉錄本,較長的Bcl-xl發(fā)揮抗凋亡的作用,而較短的Bcl-xs則能夠促進凋亡[13]。不同轉錄本存在于不同組織器官中,能構成蛋白在不同組織器官中的功能差異[14];不同轉錄本也能表達在同一組織中,并因生長環(huán)境、發(fā)育階段、病理狀況等不同而發(fā)生改變[15-17]。通過對麝鼠Elovl1基因X1和X2的蛋白理化性質與結構分析發(fā)現,兩者差異并不顯著,均編碼穩(wěn)定的親脂性蛋白,這是由于Elovl1蛋白有著相對含量較高的脂肪族氨基酸側鏈(脂肪族氨基酸指數分別為102.43和105.09),增強了蛋白的疏水性和熱穩(wěn)定性。亞細胞定位均主要位于內質網,這與超長鏈脂肪酸及麝鼠香合成場所在內質網中相吻合[18]。X1和X2均在香腺組織中同時表達,暗示2個蛋白可能具有相同或相似功能,起到提高表達量的作用,但二者的具體功能還有待進一步研究驗證。

        可變剪接在物種的進化過程中也扮演著極為重要的角色,對從不同動物的多種器官中提取的RNA進行RNA-Seq測序分析后發(fā)現,即使在同種器官中,物種間大部分的可變剪接也并不保守,存在各自獨特的剪接模式,甚至在極為相近的物種之間也存在明顯的差別,如人和黑猩猩[19-20]。從進化的角度來看,可變剪接的進化速度遠在基因本身的改變之上[21]。Elovl1基因多轉錄本編碼不同CDS的現象在嚙齒目其他物種中普遍存在,同一物種的Elovl1基因不同轉錄本間相對保守并聚在一起,但也存在不同轉錄本間已出現進化分歧的物種,如旱獺。麝鼠Elovl1基因2個轉錄本X1和X2具有物種特異性,但其是否像其他嚙齒目物種(旱獺、非洲跳鼠和豚鼠)具有3個轉錄本仍有待后續(xù)試驗檢驗。

        3.2 Elovl1基因在麝鼠香腺中的作用

        Elovl1是脂肪酸延伸過程中的關鍵酶,對多種飽和及單不飽和?;鵆oA底物表現出高活性,參與不同LCFA的生產[22-23]。細胞中LCFA的含量隨Elovl1表達變化呈正相關變化,隨著妊娠后期向泌乳期轉變,與香腺同為分泌腺的乳腺細胞中Elovl1基因表達量顯著增高,且乳汁中LCFA含量顯著增高[24]。因此,Elovl1基因在麝鼠香腺發(fā)育期和泌香盛期顯著高表達,這與乳腺泌乳期的高表達相似,提示Elovl1基因參與了麝鼠香的分泌,這與LCFA及其化合物作為麝鼠香重要組成成分相符。

        LCFA是各種脂質分子中的結構成分,如甘油磷脂、鞘脂及蠟和甾醇酯,根據它們的鏈長和不飽和度差異,其對膜的流動性和其他物理、化學性質有重要貢獻,這些LCFA作為膜脂和脂質介質的前體也參與多種生物進程。研究發(fā)現,Elovl1基因以尾對尾雙向方式與編碼細胞周期蛋白CDC20基因相連[25]。這2個基因在不同類型細胞中都會非常緊密的連接,在特定的發(fā)育階段,某些細胞可能同步表達Elovl1和CDC20,其在細胞周期中存在復雜的共同調節(jié)。Elovl1可將脂肪酸延長至26個碳原子,在正常情況下,大量LCFA與長鏈鞘氨醇相連形成神經酰胺,構成了鞘磷脂和其他鞘脂。鞘脂是真核細胞質膜的主要脂質成分,鞘脂對細胞增殖至關重要,鞘脂合成受損導致酵母和哺乳動物細胞中的細胞生長停止[26-27]。CDC20作為有絲分裂后期促進復合體的正調控因子,在細胞分裂周期中是必需的。同時,CDC20作為原癌因子,對細胞生長及腫瘤發(fā)育有促進作用[28]。因此,可能存在Elovl1和CDC20共同表達調控細胞生長發(fā)育。本研究中轉錄表達分析顯示,香腺發(fā)育期(2、3、4月末)Elovl1基因的相對表達量明顯高于其他時期,說明Elovl1參與了香腺細胞的生長發(fā)育。

        研究表明,敲除Elovl1基因后,因HeLa細胞抗刺激性下降而增加細胞凋亡易感性[29],其原因是敲除Elovl1基因后C16以上LCFA生成途徑會減少,C16的鞘脂補償性增加。而鞘脂酰基鏈長度組成的變化可能影響細胞膜環(huán)境,導致細胞凋亡(Caspase-3/7)信號通路的活化,因此,鞘脂鏈的長度改變可能會影響質膜上調節(jié)的抗凋亡信號[30];線粒體中的鞘脂可能會影響線粒體外膜通透性,從而調節(jié)細胞色素釋放以響應細胞凋亡信號[31]。麝鼠作為季節(jié)性多發(fā)情動物,香腺細胞在發(fā)育期大量增殖生長。研究發(fā)現,在香腺發(fā)育期至泌香持續(xù)初期,香腺中多種介導細胞生存或增殖的基因表達量上調;而泌香持續(xù)末期至香腺萎縮期,則多種介導細胞凋亡或自噬的基因表達量上調[32-33]。因此,Elovl1或許通過調控膜脂成分從而間接參與香腺細胞凋亡通路。

        4 結 論

        麝鼠Elovl1基因CDS編碼的蛋白是穩(wěn)定的親脂性蛋白,主要定位于內質網,且具有7個跨膜螺旋區(qū),無信號肽。Elovl1基因在麝鼠香腺中表達2個轉錄本,其表達量隨月份增加整體呈先上調再下調的趨勢。

        猜你喜歡
        脂肪酸測序物種
        吃光入侵物種真的是解決之道嗎?
        英語世界(2023年10期)2023-11-17 09:18:18
        杰 Sir 帶你認識宏基因二代測序(mNGS)
        新民周刊(2022年27期)2022-08-01 07:04:49
        二代測序協(xié)助診斷AIDS合并馬爾尼菲籃狀菌腦膜炎1例
        傳染病信息(2021年6期)2021-02-12 01:52:58
        揭開反式脂肪酸的真面目
        回首2018,這些新物種值得關注
        揭開反式脂肪酸的真面目
        電咖再造新物種
        汽車觀察(2018年10期)2018-11-06 07:05:26
        瘋狂的外來入侵物種
        鱷梨油脂肪酸組成分析
        基因捕獲測序診斷血癌
        51国产偷自视频区视频| 久久精品视频按摩| 国产性感主播一区二区| 人妻av有码中文字幕| 国产真实乱对白精彩| 男人天堂网在线视频| 中文字幕色视频在线播放| 日本一区二区三区光视频| 人妻少妇精品视频专区| 亚洲中文字幕第一页在线| 国产高清女人对白av在在线| 蜜桃高清视频在线看免费1| 久久精品国产亚洲av麻豆图片| 亚洲欧美日韩另类精品一区| 日日碰狠狠躁久久躁96avv| 国产成社区在线视频观看| 一区二区激情偷拍老牛视频av| 青青河边草免费在线看的视频| 内射人妻无套中出无码| 一本大道久久东京热无码av| 日韩av一区二区毛片| 国产自拍一区二区三区| 胸大美女又黄的网站| 97久久久久人妻精品专区| 亚洲女同一区二区久久| 国产不卡在线视频观看| 亚洲午夜无码av毛片久久| 久久亚洲国产成人亚| 亚洲精品国产成人久久av盗摄| 亚洲国产中文字幕视频| 全部孕妇毛片丰满孕妇孕交| 香蕉视频免费在线| 蜜桃在线视频一区二区| 国产av夜夜欢一区二区三区| 少妇高潮喷水正在播放 | 亚洲中文字幕精品久久久久久直播 | 少妇高潮久久蜜柚av| 日本最新免费二区| av中文字幕综合在线| 国产99久久久国产精品免费| 欧美丰满少妇xxxx性|