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

    異尾次目寄居蟹總科線粒體基因組比較分析及系統(tǒng)發(fā)育研究

    2024-02-07 07:29:56江婷琪胡璧麟張楠楠汪凱欣呂振明
    水生生物學(xué)報(bào) 2024年2期
    關(guān)鍵詞:寄居蟹重排線粒體

    江婷琪 胡璧麟 張楠楠 汪凱欣 呂振明 龔 理

    (浙江海洋大學(xué)海洋科學(xué)與技術(shù)學(xué)院, 海洋生物種質(zhì)發(fā)掘與利用國家地方聯(lián)合工程研究中心, 舟山 316022)

    異尾次目(Anomura)隸屬于節(jié)肢動(dòng)物門(Arthropoda)軟甲綱(Malacostraca)十足目(Decapoda), 是一類體型介于蝦類和蟹類之間高度特化的甲殼動(dòng)物。根據(jù)形態(tài)特征, 異尾次目可分為鎧甲蝦型(Squat lobster form)、蟹型(Crab-like form)和寄居蟹型(Hermit form)三個(gè)主要類群; 其中, 寄居蟹型又可以分為對(duì)稱寄居蟹型(Symmetrical hermit form)和不對(duì)稱寄居蟹型(Asymmetrical hermit form)。對(duì)稱寄居蟹型常見于門螯寄居蟹科(Pylochelidae)種類, 不對(duì)稱寄居蟹型多見于陸寄居蟹科(Coenobitidae)、擬寄居蟹科(Parapaguridae)和寄居蟹科(Paguridae)絕大多數(shù)種類。寄居蟹獨(dú)特的形態(tài)特征代表了甲殼動(dòng)物從長尾類(Macrura)到短尾類(Brachyura)的過渡狀態(tài), 因而在甲殼動(dòng)物進(jìn)化過程中占據(jù)十分重要的地位[1]。

    雖然1802年Latreille就建立了寄居蟹總科(Paguridea), 但是該總科分類系統(tǒng)長期以來未有定論[2,3]。最初的研究把所有寄居蟹都?xì)w于寄居蟹總科, 直至1957年, MacDonald等[4]提出將寄居蟹分為寄居蟹總科和陸寄居蟹總科, 其中寄居蟹總科包括石蟹科(Lithodidae)、擬寄居蟹科和寄居蟹科, 陸寄居蟹總科包括陸寄居蟹科和活額寄居蟹科(Diogenidae),該分類體系并未考慮澳洲寄居蟹科(Lomisidae) 的歸屬。McLaughlin[5]則認(rèn)為澳洲寄居蟹科應(yīng)該歸屬獨(dú)立總科, 同時(shí)提議廢除陸寄居蟹總科, 僅保留寄居蟹總科。Martin和Davis[6]支持僅設(shè)寄居蟹總科的觀點(diǎn), 并新增門螯寄居蟹科。McLaughlin等[7]基于形態(tài)特征重建了寄居蟹總科系統(tǒng)發(fā)育樹, 提議將石蟹科獨(dú)立出來, 并將其歸于石蟹總科(Lithodoidea)。當(dāng)前, 石蟹科到底歸于寄居蟹總科還是石蟹總科仍有爭議。由于螯蓋寄居蟹科(Pylojacquesidae)物種稀少(目前僅發(fā)現(xiàn)2屬2種), 直至2001年該科才被建立。至此, 寄居蟹總科包括擬寄居蟹科、寄居蟹科、陸寄居蟹科、活額寄居蟹科、門螯寄居蟹科和螯蓋寄居蟹科6科的分類系統(tǒng)正式確立。

    雖然, 目前寄居蟹總科的分類系統(tǒng)被絕大多數(shù)學(xué)者所認(rèn)可[8—10], 但是該類群內(nèi)部分類及親緣關(guān)系仍具爭議。Richter和Scholz[3]首次基于形態(tài)特征分析了寄居蟹類的進(jìn)化關(guān)系, 認(rèn)為寄居蟹總科中除了最原始的門螯寄居蟹是并系群外, 其余種類均為單系群; 而當(dāng)前絕大多數(shù)分子數(shù)據(jù)都支持活額寄居蟹科為并系群[11,12]。隨著寄居蟹科分子數(shù)據(jù)不斷增加, 越來越多的研究發(fā)現(xiàn)該科也為非單系群, 如寄居蟹科中的Pagurus longicarpus不與同屬物種聚為一支, 而是與寄居蟹科和石蟹科組成的姐妹支聚支[10,13],支持寄居蟹科的非單系性。此外, 寄居蟹和石蟹之間的進(jìn)化關(guān)系也是甲殼動(dòng)物系統(tǒng)進(jìn)化學(xué)家關(guān)注的焦點(diǎn)[8,14,15]。因此, 開展寄居蟹總科系統(tǒng)分類學(xué)研究對(duì)異尾次目乃至整個(gè)甲殼類分類及系統(tǒng)演化具有重要意義。

    線粒體基因組被廣泛應(yīng)用于分子系統(tǒng)學(xué)研究,解決了很多長期以來備受爭議的分類鑒定及系統(tǒng)進(jìn)化等問題。如傳統(tǒng)形態(tài)學(xué)認(rèn)為沙蟹總科和方蟹總科是單系群, 但是越來越多的分子結(jié)果(包括線粒體基因組數(shù)據(jù))顯示這兩個(gè)總科為并系群[16—19]。除了線粒體基因組序列本身外, 也有研究發(fā)現(xiàn)線粒體基因排序也可以為分子系統(tǒng)學(xué)提供有用信息。如Tan等[20]比較了已有異尾次目和短尾次目線粒體全序列后提出基因重排可以用于異尾次目系統(tǒng)發(fā)育研究, 如前期基于線粒體基因組序列構(gòu)建的系統(tǒng)發(fā)育樹無法確定輝蝦屬(Aegla)(輝蝦總科)和澳洲寄居蟹屬(Lomis)(澳洲寄居蟹總科)之間的親緣關(guān)系,但是基因排序信息則顯示澳洲寄居蟹總科和柱螯蝦總科有著最近的親緣關(guān)系, 兩者再與輝蝦總科形成姐妹群關(guān)系, 隨后不斷有研究證實(shí)了這一觀點(diǎn)[12,13]。關(guān)于線粒體基因重排現(xiàn)象, 目前主要用串聯(lián)復(fù)制-隨機(jī)丟失(Tandem duplication and random loss)[21]、串聯(lián)復(fù)制-非隨機(jī)丟失(Tandem duplication and nonrandom loss)[22]和線粒體內(nèi)重組(Intramitochondrial recombination)[23]等模型假說來解釋。

    然而, 相比甲殼類其他類群, 研究人員對(duì)異尾次目線粒體基因組關(guān)注度明顯不足; 尤其是寄居蟹總科, 截2023年4月, GenBank數(shù)據(jù)庫中僅有16個(gè)物種的線粒體基因組全序列, 這極大阻礙了寄居蟹分子系統(tǒng)學(xué)研究發(fā)展。因此, 本研究以我國海域分布最廣泛的活額寄居蟹科物種作為研究對(duì)象, 以刺足真寄居蟹(Dardanus hessii)為代表種, 首次測(cè)定其線粒體基因組全序列, 并結(jié)合已公布的16個(gè)寄居蟹物種的基因組全序列, 對(duì)寄居蟹總科線粒體基因組序列特征及基因排序進(jìn)行了比較分析, 同時(shí)構(gòu)建了寄居蟹總科系統(tǒng)發(fā)育樹, 探究線粒體基因重排在系統(tǒng)進(jìn)化研究中的適用性。本研究結(jié)果不僅豐富了寄居蟹總科的線粒體基因組數(shù)據(jù), 同時(shí)為異尾次目分類及系統(tǒng)進(jìn)化提供新的觀點(diǎn)和思路。

    1 材料與方法

    1.1 樣品采集及測(cè)序

    刺足真寄居蟹樣品采自浙江省舟山市桃花島,參考《中國海活額寄居蟹科分類學(xué)研究》[24]進(jìn)行形態(tài)學(xué)鑒定。取約100 g肌肉組織利用二代測(cè)序技術(shù)進(jìn)行線粒體基因組全序列測(cè)定(上海元莘生物醫(yī)藥科技有限公司)。測(cè)序前先使用插入片段大小為400 bp的VAHTS Universal Plus DNA文庫制備試劑盒構(gòu)建基因文庫, 在Illumina NovaSeq 6000高通量測(cè)序平臺(tái)進(jìn)行PE150測(cè)序; 使用NOVO Plasty[25]軟件對(duì)去除接頭的有效數(shù)據(jù)(Clean data)進(jìn)行從頭組裝; 將鱗紋真寄居蟹線粒體基因組(GenBank登錄號(hào): MW147148)作為種子序列進(jìn)行序列延伸。為了評(píng)估二代測(cè)序的準(zhǔn)確性, 研究還利用一代測(cè)序技術(shù)(Sanger測(cè)序)測(cè)定了線粒體COⅠ基因序列, 比較兩種測(cè)序方法得到的序列相似性。

    1.2 線粒體基因組注釋和序列分析

    基于MITOS Web Server[26]及tRNAscan-SE 1.21[27]預(yù)測(cè)結(jié)果, 利用Sequin軟件(version 15.10, http://www.ncbi.nlm.nih.gov/Sequin/)對(duì)新組裝的刺足真寄居蟹線粒體基因組進(jìn)行注釋; 使用NCBI-BLAST對(duì)蛋白質(zhì)編碼基因和核糖體RNA基因進(jìn)行邊界分析; 控制區(qū)由于序列變異較大, 根據(jù)近緣種序列及其相鄰基因位置判斷。利用以下公式計(jì)算堿基的鏈不對(duì)稱性: GC-skew=(G-C)/(G+C); AT-skew=(A-T)/(A+T)[28]。使用 MEGA X軟件[29]計(jì)算核苷酸含量。從GenBank 數(shù)據(jù)庫中下載另外16種寄居蟹總科物種的線粒體COⅠ基因序列, 結(jié)合本研究測(cè)定的刺足真寄居蟹物種COⅠ序列, 利用軟件Clustal X 2.0[30]對(duì)這些序列進(jìn)行比對(duì), 并基于Kimura 2-parameter (K2P)參數(shù)替代模型計(jì)算17個(gè)物種COⅠ的遺傳距離, 利用軟件MatGAT 2.02[31]進(jìn)行多序列相似性比較分析。

    1.3 線粒體基因組比較分析

    以刺足真寄居蟹線粒體基因組作為參考序列,用CGView Comparison Tool (CCT)[32]軟件對(duì)所有17種寄居蟹總科物種線粒體基因組序列進(jìn)行比較,同時(shí)用Mauve v2.4.0[33]軟件對(duì)其進(jìn)行共線性分析。采用CREx[34]軟件預(yù)測(cè)基因重排事件, 基因重排包括倒置(R)、移位(T)和倒置移位(RT)三種類型。

    1.4 系統(tǒng)發(fā)育分析

    從GenBank數(shù)據(jù)庫下載目前已有16種寄居蟹總科物種線粒體基因組全序列, 加上本研究新測(cè)定的刺足真寄居蟹基因組, 以方蟹總科的側(cè)足厚蟹(Helice latimera)和伍氏擬厚蟹(Helicana wuana)為外類群, 基于13個(gè)蛋白質(zhì)編碼基因, 使用PhyloSuite[35]同時(shí)構(gòu)建寄居蟹總科最大似然樹(ML)和貝葉斯樹(BI)。最大似然樹使用IQ-TREE構(gòu)建, 其中bootstrap設(shè)置為Ultrafast, Num of bootstrap設(shè)置為100000,抽樣次數(shù)為1000。貝葉斯樹使用MrBayes構(gòu)建, 使用4條蒙特卡羅馬爾可夫鏈(Markov Chain Monte Carlo, MCMC)同時(shí)運(yùn)行2000000代, 抽樣頻率為1000, 摒棄25%的老化樣本。

    2 結(jié)果

    2.1 刺足真寄居蟹線粒體基因組序列特征

    刺足真寄居蟹線粒體基因組全序列長度為17821 bp(GenBank登錄號(hào)OP930882), 包含37個(gè)基因(13個(gè)蛋白質(zhì)編碼基因, 22個(gè)tRNA和2個(gè)rRNA基因)和1個(gè)控制區(qū)。13個(gè)蛋白質(zhì)編碼基因總長為11183 bp, 共編碼3716個(gè)氨基酸。所有蛋白質(zhì)編碼基因均以典型的ATN作為起始密碼子; 除ND5和ND6基因以T作為終止密碼子外, 其余蛋白質(zhì)編碼基因以TAA或TAG作為終止密碼子(表1)。tRNA基因分散在整個(gè)線粒體基因組中, 22個(gè)tRNA的總長度為1471 bp。16S rRNA位于ND1和tRNA-Val之間, 長度為1393 bp;12S rRNA位于tRNA-Val和CR之間, 長度為798 bp。雖然線粒體基因組各基因排列非常緊湊, 但是在刺足真寄居蟹基因組中23處共發(fā)現(xiàn)270 bp的基因間隔區(qū), 其中最長的為75 bp, 介于tRNA-His和ND4之間。當(dāng)然, 同時(shí)也在5處發(fā)現(xiàn)了累計(jì)15 bp的重疊區(qū)域, 最長的重疊區(qū)為7 bp, 位于ND4和ND4L之間(表1)。

    表1 刺足真寄居蟹線粒體基因組特征Tab.1 Features of the mitochondrial genome of D.hessii

    2.2 寄居蟹總科物種線粒體基因組比較分析

    用CCT軟件對(duì)17種寄居蟹總科的線粒體基因組進(jìn)行了比較分析, 結(jié)果顯示, 17個(gè)物種的線粒體基因組都發(fā)生了基因重排, 且基因排序和序列相似性都相對(duì)保守。除控制區(qū)和16S rRNA部分序列有較大變異外, 其余基因相似性普遍高于82%, 尤其是COⅠ和COⅢ基因在所有寄居蟹物種線粒體基因組中十分保守(圖1)。

    COⅠ基因作為最常用的DNA條形碼, 在物種分類鑒定有著廣泛的應(yīng)用。因此本研究計(jì)算了17種寄居蟹總科物種的COⅠ基因遺傳距離及其序列相似性。結(jié)果顯示, 日本寄居蟹(Pagurus japonicus)和長腕寄居蟹(Pagurus filholi)間的遺傳距離最小(0.005), 二者對(duì)應(yīng)的序列相似性也最高(99.5%), 表明兩者很有可能為同一物種; 最大的遺傳距離出現(xiàn)在Pagurus longicarpus和灰白陸寄居蟹(Coenobita rugosus)之間(0.292), 對(duì)應(yīng)的序列相似性為75.7%。本研究新測(cè)定的刺足真寄居蟹和同屬的紅星真寄居蟹(Dardanus asperses)及鱗紋真寄居蟹(Dardanus arrosor)的遺傳距離最小, 分別為0.203和0.204, 序列相似性分別為82.1%和82.0%。

    研究還對(duì)17種寄居蟹總科物種的線粒體基因組共線性進(jìn)行了比較分析, 軟件分析顯示, 所有寄居蟹物種線粒體基因組均可分為5個(gè)大的保守區(qū)域(圖2, A—E), 其中區(qū)域A(P-ND1-16S rRNA) 的位置和長度在所有寄居蟹物種線粒體基因組中都是最保守的, 區(qū)域B的位置和長度變異最大, 其余3個(gè)區(qū)域則相對(duì)保守。值得注意的是毛氏門螯寄居蟹(Pylocheles mortensenii) 基因組中的區(qū)域B(12S rRNA)發(fā)生了倒置, 由輕鏈編碼變成了重鏈編碼, 而在其余16種寄居蟹總科物種線粒體基因組均仍由輕鏈編碼; 此外, 區(qū)域E (T-ND6-Cytb-S2) 在紅星真寄居蟹基因組中發(fā)生了較明顯的移位, 而其在其余16種寄居蟹總科物種線粒體基因組中位置十分保守。

    圖2 17種寄居蟹總科物種線粒體基因組共線性分析Fig.2 Collinearity analysis of 17 Paguroidea mitogenomes

    進(jìn)一步分析發(fā)現(xiàn)17種寄居蟹總科物種的線粒體基因組共分為4種共線性類型, 其中7種寄居蟹科物種共享一種共線類型(1—7); 5種陸寄居蟹科和3種活額寄居蟹科物種共享一種共線類型(8—15);活額寄居蟹科的紅星真寄居蟹獨(dú)享一種共線類型(16), 與該科其他3個(gè)物種線粒體基因排序明顯不同, 表現(xiàn)在紅星真寄居蟹基因組中T-ND6-Cytb-S2移位至R-N-E基因簇下游, 而在5種陸寄居蟹和3種活額寄居蟹基因組中, 該區(qū)域移位至ND4L下游;門螯寄居蟹科的毛氏門螯寄居蟹獨(dú)享一種共線類型(17), 該類型與其他3種類型最大的區(qū)別是12S rRNA發(fā)生了倒置。

    2.3 線粒體基因重排分析

    與泛甲殼動(dòng)物(Pancrustacea) 線粒體基因組原始基因排序(Type 0)相比, 刺足真寄居蟹基因順序經(jīng)歷了大規(guī)模的重排, 共有6處基因或基因簇位置發(fā)生了顯著變化, 涉及到12個(gè)tRNA基因(L2、G、A、S1、P、L1、I、Q、M、W、C和Y)和2個(gè)蛋白質(zhì)編碼基因(ND3和ND2)。在這6處重排中, 包含了倒置移位(RT)和移位(T)兩種重排類型。單個(gè)的L2和L1基因、G-ND3-A-S1和W-C-Y基因簇發(fā)生了倒置移位, 單個(gè)的P基因和I-Q-M-ND2基因簇僅發(fā)生移位,未發(fā)生倒置; 其中L2基因由重鏈編碼倒轉(zhuǎn)至輕鏈,同時(shí)從COⅠ基因下游移至G基因的下游(圖3 ①);G-ND3-A-S1基因簇也是由重鏈編碼倒轉(zhuǎn)至輕鏈, 變成了S1-A-G-ND3排序, 位置從COⅢ基因的下游移至了控制區(qū)(CR)的下游(圖3 ②);P基因未發(fā)生倒置, 僅發(fā)生位置上的變化, 從T基因的下游移至S2基因的下游(圖3 ③);L1基因則由輕鏈編碼倒轉(zhuǎn)至重鏈, 位置從ND1基因的下游移位至COⅠ基因的下游(圖3 ④);I-Q-M-ND2基因簇僅位置發(fā)生變化, 其被分為兩部分, 其中I-M-ND2整體移至K基因的下游, 且I和M的前后位置發(fā)生了互換, 形成了M-IND2排序,Q基因則移至線粒體基因組的線性末端位置(圖3 ⑤); 最后一處重排發(fā)生在W-C-Y基因簇,其中W和Y基因發(fā)生了倒置移位, 形成了新的Y-WC排序(圖3 ⑥)。

    圖3 寄居蟹總科物種線粒體重排類型(A) 及寄居蟹總科系統(tǒng)發(fā)育關(guān)系(B)Fig.3 Mitochondrial rearrangement types of species of Paguridea (A) and phylogenetic relationship of Paguridea (B)

    研究進(jìn)一步比較了所有寄居蟹總科物種線粒體基因排序, 結(jié)果顯示17種寄居蟹總科物種線粒體基因組表現(xiàn)出7種不同類型的基因重排, 其中寄居蟹科物種線粒體基因組貢獻(xiàn)了4種基因重排類型(圖3Ⅰ—Ⅳ); 所有陸寄居蟹科物種和3種活額寄居蟹科共享一種重排類型(圖3Ⅴ); 活額寄居蟹科的紅星真寄居蟹和門螯寄居蟹科的毛氏門螯寄居蟹各自獨(dú)享一種重排類型(圖3Ⅵ和Ⅶ)。相比其余3種基因排序(圖3Ⅴ—Ⅶ) 之間的差異, 寄居蟹科中4種基因重排類型(圖3Ⅰ—Ⅳ) 內(nèi)部差異顯得很小,僅在AIMD基因簇排序、H及V的位置上有差異?;铑~寄居蟹科具有兩種不同的基因重排類型(圖3Ⅴ和Ⅵ), 兩者之間差異僅體現(xiàn)在T-ND6-Cytb-S2基因簇的位置, 包括大部分活額寄居蟹科在內(nèi)的寄居蟹總科物種T-ND6-Cytb-S2基因簇位于ND4L基因下游, 僅在活額寄居蟹科的紅星真寄居蟹基因組中位于R-N-E基因簇下游。門螯寄居蟹科基因排序(圖3Ⅶ)最大的不同在于12S rRNA的編碼鏈, 在其余6種基因重排類型中, 12S rRNA均由輕鏈編碼, 而該基因在門螯寄居蟹科線粒體基因組中發(fā)生了倒置, 由輕鏈編碼轉(zhuǎn)為重鏈編碼, 與線粒體基因組共線性分析結(jié)果一致。

    2.4 系統(tǒng)發(fā)育分析

    本研究基于13個(gè)蛋白質(zhì)編碼基因的核苷酸序列同時(shí)構(gòu)建了寄居蟹總科最大似然樹(ML tree)和貝葉斯樹(BI tree), 結(jié)果顯示兩種方法構(gòu)建的系統(tǒng)發(fā)育樹具有相同的拓?fù)浣Y(jié)構(gòu); 因此, 本研究只顯示其中一種樹, 節(jié)點(diǎn)處同時(shí)標(biāo)記出貝葉斯后驗(yàn)概率和最大似然樹自展值。系統(tǒng)樹顯示本研究新測(cè)定的刺足真寄居蟹跟鱗紋真寄居蟹及紅星真寄居蟹聚為一支, 表明這三者親緣關(guān)系很近, 組成了活額寄居蟹科的一部分; 然而, 活額寄居蟹科的另一個(gè)物種, 下齒細(xì)螯寄居蟹并沒有與同科的3個(gè)物種聚在一起, 而是先與陸寄居蟹科聚為一支, 再與3種活額寄居蟹科物種形成姐妹群, 表明活額寄居蟹科為非單系群。寄居蟹總科內(nèi)部親緣關(guān)系為活額寄居蟹科與陸寄居蟹科形成姐妹群后再與寄居蟹科聚在一起, 最后與門螯寄居蟹科聚支。需要提醒的是,NCBI數(shù)據(jù)庫中長腕寄居蟹(LC222528)與日本寄居蟹(LC222532)的線粒體基因組序列可能來自同一物種, 與COⅠ遺傳距離和序列相似性結(jié)果一致。兩者形態(tài)上有較大差異(如楯部比例、螯上小刺或短剛毛等), 暗示其中至少有一個(gè)物種鑒定有誤。

    3 討論

    異尾次目分類及演化是進(jìn)化系統(tǒng)學(xué)備受關(guān)注的科學(xué)問題。然而, 該類群分類及內(nèi)部親緣關(guān)系長期以來一直缺乏定論[8,15,36,37]。隨著測(cè)序技術(shù)的飛速發(fā)展及測(cè)序成本的不斷降低, 線粒體基因組全序列在分類鑒定、系統(tǒng)進(jìn)化和適應(yīng)性進(jìn)化等領(lǐng)域解決了很多備受爭議的難題[16,20,38]。本研究基于17種寄居蟹總科物種13個(gè)線粒體蛋白質(zhì)編碼基因序列重建了寄居蟹總科系統(tǒng)發(fā)育關(guān)系, 結(jié)果顯示活額寄居蟹科與陸寄居蟹科親緣關(guān)系最近, 兩者形成姐妹群后再與寄居蟹科聚為一支, 最后與門螯寄居蟹科聚在一起; 該聚類關(guān)系與絕大多數(shù)分子系統(tǒng)樹結(jié)果一致[12,20,39,40]。雖然, 本研究也支持了活額寄居蟹科的非單系性, 但是值得注意的是造成該科非單系的關(guān)鍵物種, 下齒細(xì)螯寄居蟹與陸寄居蟹科的關(guān)系與前人研究結(jié)果略有區(qū)別。在近期的異尾次目系統(tǒng)發(fā)育研究中, 陸寄居蟹科先與大多數(shù)活額寄居蟹科聚支后再與下齒細(xì)螯寄居蟹聚在一起[11], 而本研究并入刺足真寄居蟹后則顯示陸寄居蟹科先與下齒細(xì)螯寄居蟹聚為一支, 再與大多數(shù)活額寄居蟹科聚在一起, 該結(jié)果與Tsang等[41]基于5個(gè)核蛋白編碼基因構(gòu)建的寄居蟹總科系統(tǒng)發(fā)育關(guān)系一致。不難發(fā)現(xiàn)下齒細(xì)螯寄居蟹在系統(tǒng)樹中的位置并不完全穩(wěn)定, 今后可增加活額寄居蟹科和寄居蟹科代表種數(shù)量進(jìn)一步厘清這兩個(gè)科的親緣關(guān)系。

    脊椎動(dòng)物線粒體基因組排序相對(duì)穩(wěn)定, 發(fā)生重排的概率較低。而在無脊椎動(dòng)物線粒體基因組中,線粒體基因重排現(xiàn)象相對(duì)較多。以甲殼動(dòng)物十足目為例, 目前已發(fā)現(xiàn)相手蟹科(Sesarmidae)、溪蟹科(Potamonidae)、寶石蟹科(Mithracidae)、寄居蟹科、陸寄居蟹科和活額寄居蟹科等20余科物種線粒體基因組發(fā)生顯著重排[16,42,43]。用于解釋線粒體基因重排的機(jī)制主要有串聯(lián)復(fù)制隨機(jī)丟失模型、串聯(lián)復(fù)制非隨機(jī)丟失模型和線粒體重組等[23,44]。本研究中寄居蟹總科線粒體基因組呈現(xiàn)7種不同的基因排序, 涉及倒置移位和移位兩種重排類型。線粒體重組是最常用于解釋倒置現(xiàn)象的模型, 已廣泛用于解釋包括魚類[45]、爬行類[46]、昆蟲[47]和甲殼類[10]等在內(nèi)多種生物類群基因倒置現(xiàn)象。如半滑舌鰨(Cynoglossus semilaevis)線粒體基因組中Q基因由輕鏈編碼轉(zhuǎn)為重鏈編碼, 這是魚類線粒體基因組中首例基因倒置現(xiàn)象, 研究采用線粒體重組模型對(duì)Q基因倒置進(jìn)行了重排推演[40,42]。本研究中的倒置現(xiàn)象同樣也可以用重組模型來解釋。如在泛甲殼動(dòng)物線粒體基因組(Type 0)中為V-12S排序, 兩者均由輕鏈編碼; 而在毛氏門螯寄居蟹基因組(TypeⅦ)中變?yōu)?2S-V排序, 且由輕鏈編碼轉(zhuǎn)為重鏈編碼??赡艿闹嘏胚^程為:V的5′端和12S的3′端同時(shí)斷裂, 形成一個(gè)自由的V-12S片段, 該片段在重新連接時(shí)發(fā)生鏈倒置, 變成12S-V排序, 且由原來的輕鏈編碼變?yōu)橹劓?。線粒體重組模型在其他物種類似的倒置現(xiàn)象中也得到廣泛應(yīng)用[38,39]。線粒體基因移位是十分常見的基因重排現(xiàn)象, 常用串聯(lián)復(fù)制隨機(jī)丟失模型和串聯(lián)復(fù)制非隨機(jī)丟失模型來解釋。兩種模型最大的區(qū)別在于復(fù)制后的基因被刪除的方式, 前者為隨機(jī)刪除, 后者則取決于基因的轉(zhuǎn)錄極性和位置, 其刪除后的結(jié)果是具有相同極性的基因(由同一鏈編碼的基因)聚在一起。本研究中寄居蟹總科重排基因是散狀分布在基因組中, 而不是按相同極性聚在一起, 因而綜合重排基因序列特征及CREx軟件預(yù)測(cè)結(jié)果, 我們認(rèn)為串聯(lián)復(fù)制隨機(jī)丟失模型是用來解釋寄居蟹總科線粒體基因組產(chǎn)生大規(guī)?;蛑嘏努F(xiàn)象最合理的假說。

    研究還分析了寄居蟹總科線粒體基因排序跟系統(tǒng)進(jìn)化之間的關(guān)系。寄居蟹科雖然擁有4種不同的基因重排類型(Type Ⅰ—Ⅳ), 但是整科物種還是聚為一支, 形成一個(gè)單系群(Clade Ⅰ)。陸寄居蟹科和活額寄居蟹科中3個(gè)物種共享同一種重排類型(Type Ⅴ), 它們最先聚為一支, 暗示基因排序一定程度上能夠反映系統(tǒng)進(jìn)化信息[40,43]。但是值得注意的是活額寄居蟹科中鱗紋真寄居蟹先與基因排序具有較大區(qū)別的紅星真寄居蟹聚類, 兩者形成姐妹支后再與具有相同基因排序的刺足真寄居蟹聚在一起, 這又與相同基因排序優(yōu)先聚類的現(xiàn)象相矛盾。造成該“矛盾”現(xiàn)象的原因是紅星真寄居蟹基因重排并未對(duì)線粒體基因序列產(chǎn)生大的影響, 因而在基于線粒體序列構(gòu)建系統(tǒng)發(fā)育關(guān)系時(shí), 物種親緣關(guān)系并未按基因重排類型進(jìn)行聚類。因此, 基因排序是否適用于系統(tǒng)進(jìn)化分析應(yīng)當(dāng)視具體情況而定;即當(dāng)基因重排會(huì)對(duì)序列本身產(chǎn)生較大影響時(shí), 重排信息可能會(huì)成為一種理想的分子標(biāo)記; 反之基因重排可能不適用于系統(tǒng)進(jìn)化分析。

    猜你喜歡
    寄居蟹重排線粒體
    愛“搬家”的寄居蟹
    寄居蟹
    大學(xué)有機(jī)化學(xué)中的重排反應(yīng)及其歸納教學(xué)實(shí)踐
    棘皮動(dòng)物線粒體基因組研究進(jìn)展
    線粒體自噬與帕金森病的研究進(jìn)展
    重排濾波器的實(shí)現(xiàn)結(jié)構(gòu)*
    EGFR突變和EML4-ALK重排雙陽性非小細(xì)胞肺癌研究進(jìn)展
    寄居蟹
    寄居蟹
    基于像素重排比對(duì)的灰度圖彩色化算法研究
    最近最新免费中文字幕在线| 久久精品国产99精品国产亚洲性色 | 国产xxxxx性猛交| 国产精品香港三级国产av潘金莲| 99re6热这里在线精品视频| 一夜夜www| 波多野结衣一区麻豆| 欧美大码av| 成年版毛片免费区| 人妻 亚洲 视频| 亚洲精品国产精品久久久不卡| 国内毛片毛片毛片毛片毛片| 一区在线观看完整版| av一本久久久久| 美女国产高潮福利片在线看| 日韩精品免费视频一区二区三区| 欧美日韩中文字幕国产精品一区二区三区 | 日本撒尿小便嘘嘘汇集6| 国产av国产精品国产| 久久久久网色| 80岁老熟妇乱子伦牲交| 国产精品 国内视频| 高清av免费在线| 欧美 亚洲 国产 日韩一| 自线自在国产av| 国产成人系列免费观看| 亚洲熟妇熟女久久| 三级毛片av免费| 99精品在免费线老司机午夜| 大型黄色视频在线免费观看| 国产精品自产拍在线观看55亚洲 | 国产日韩一区二区三区精品不卡| 久久天躁狠狠躁夜夜2o2o| 国产精品亚洲av一区麻豆| netflix在线观看网站| 男人舔女人的私密视频| 国产男女超爽视频在线观看| 夜夜夜夜夜久久久久| 中文字幕av电影在线播放| 另类精品久久| 国产精品.久久久| 99香蕉大伊视频| 亚洲欧美日韩高清在线视频 | 欧美人与性动交α欧美精品济南到| 一级毛片电影观看| 日本黄色视频三级网站网址 | 又大又爽又粗| 中文字幕最新亚洲高清| 王馨瑶露胸无遮挡在线观看| 亚洲午夜理论影院| 国产精品国产高清国产av | 日日摸夜夜添夜夜添小说| 国产亚洲欧美精品永久| 国产熟女午夜一区二区三区| 免费观看a级毛片全部| 18禁裸乳无遮挡动漫免费视频| 国产黄色免费在线视频| 欧美日韩视频精品一区| 老鸭窝网址在线观看| 久久精品亚洲av国产电影网| tocl精华| 国产麻豆69| 欧美一级毛片孕妇| 可以免费在线观看a视频的电影网站| 狂野欧美激情性xxxx| 国产精品久久电影中文字幕 | 一本—道久久a久久精品蜜桃钙片| 国产一区二区 视频在线| 国产老妇伦熟女老妇高清| 国产精品久久久av美女十八| 最近最新免费中文字幕在线| 精品第一国产精品| 在线观看www视频免费| 久久av网站| 两性夫妻黄色片| 成人18禁高潮啪啪吃奶动态图| 亚洲av片天天在线观看| 丝袜喷水一区| 亚洲精品自拍成人| 两人在一起打扑克的视频| 精品少妇久久久久久888优播| 咕卡用的链子| 久久久久久人人人人人| 国内毛片毛片毛片毛片毛片| 大码成人一级视频| 中文亚洲av片在线观看爽 | 久久人妻福利社区极品人妻图片| 夜夜夜夜夜久久久久| 大片免费播放器 马上看| 女人被躁到高潮嗷嗷叫费观| 久久久水蜜桃国产精品网| 操出白浆在线播放| 在线 av 中文字幕| 高清视频免费观看一区二区| 国产成人一区二区三区免费视频网站| netflix在线观看网站| 色在线成人网| 亚洲人成电影观看| 午夜福利视频精品| 97在线人人人人妻| 国产精品久久久久久人妻精品电影 | 老司机靠b影院| 岛国毛片在线播放| 久久久久视频综合| 一级a爱视频在线免费观看| 国产不卡av网站在线观看| 一区二区三区精品91| 后天国语完整版免费观看| 极品教师在线免费播放| 欧美国产精品一级二级三级| 精品国产超薄肉色丝袜足j| 99国产精品免费福利视频| 久久久国产成人免费| 色播在线永久视频| 午夜福利欧美成人| 国产av又大| 男男h啪啪无遮挡| 日本五十路高清| 叶爱在线成人免费视频播放| 国产高清国产精品国产三级| 欧美 日韩 精品 国产| 国产精品九九99| 性少妇av在线| av不卡在线播放| 日韩欧美一区二区三区在线观看 | 美女主播在线视频| 国产精品熟女久久久久浪| 考比视频在线观看| 无人区码免费观看不卡 | 黄频高清免费视频| 成人18禁高潮啪啪吃奶动态图| 久久久久网色| 欧美激情极品国产一区二区三区| 女人精品久久久久毛片| 国产精品久久久久久精品电影小说| 精品国产一区二区三区四区第35| 亚洲国产欧美一区二区综合| 日日摸夜夜添夜夜添小说| 精品亚洲成a人片在线观看| 人成视频在线观看免费观看| 天堂8中文在线网| 91字幕亚洲| 免费观看人在逋| 亚洲国产中文字幕在线视频| 国产高清视频在线播放一区| 中文字幕人妻丝袜制服| 久久热在线av| 99国产极品粉嫩在线观看| 在线看a的网站| 国产高清videossex| 建设人人有责人人尽责人人享有的| 亚洲第一青青草原| 国产av又大| 咕卡用的链子| 国产精品国产av在线观看| 香蕉国产在线看| 9191精品国产免费久久| 久久国产精品男人的天堂亚洲| 久久99一区二区三区| 成人影院久久| 老熟女久久久| 亚洲精品自拍成人| 国产在视频线精品| 欧美午夜高清在线| 在线av久久热| 欧美亚洲日本最大视频资源| 亚洲九九香蕉| 午夜免费鲁丝| 国产xxxxx性猛交| 国产精品成人在线| netflix在线观看网站| 久久亚洲精品不卡| 中文欧美无线码| 99精品在免费线老司机午夜| 天天躁狠狠躁夜夜躁狠狠躁| 黑丝袜美女国产一区| 美女主播在线视频| 大码成人一级视频| 9191精品国产免费久久| 不卡一级毛片| 一个人免费看片子| 十八禁网站免费在线| 亚洲人成77777在线视频| 一区二区av电影网| 亚洲欧洲日产国产| 夜夜骑夜夜射夜夜干| 久热爱精品视频在线9| 国产一区有黄有色的免费视频| 巨乳人妻的诱惑在线观看| 国产精品.久久久| 99精品久久久久人妻精品| www.熟女人妻精品国产| 水蜜桃什么品种好| 男女午夜视频在线观看| 成人亚洲精品一区在线观看| 99香蕉大伊视频| 亚洲国产欧美网| 亚洲美女黄片视频| 亚洲精品自拍成人| 免费在线观看视频国产中文字幕亚洲| 一进一出好大好爽视频| 亚洲第一欧美日韩一区二区三区 | 亚洲人成77777在线视频| 欧美国产精品一级二级三级| 久久人人97超碰香蕉20202| 高清黄色对白视频在线免费看| 黑人猛操日本美女一级片| 女人被躁到高潮嗷嗷叫费观| 又大又爽又粗| a级毛片黄视频| 亚洲精品国产精品久久久不卡| 国产单亲对白刺激| 12—13女人毛片做爰片一| 一夜夜www| 国产99久久九九免费精品| 久久国产精品男人的天堂亚洲| 高清毛片免费观看视频网站 | 热99re8久久精品国产| 一边摸一边抽搐一进一出视频| 国产伦理片在线播放av一区| 久久久久久亚洲精品国产蜜桃av| 精品国产亚洲在线| 久久精品人人爽人人爽视色| 国产成人av激情在线播放| 国产又爽黄色视频| 亚洲国产av新网站| 麻豆成人av在线观看| 久久久国产精品麻豆| 男女高潮啪啪啪动态图| 乱人伦中国视频| 汤姆久久久久久久影院中文字幕| 亚洲专区字幕在线| 成人免费观看视频高清| 欧美国产精品一级二级三级| 天天影视国产精品| 日本wwww免费看| 欧美精品啪啪一区二区三区| 日韩熟女老妇一区二区性免费视频| 777米奇影视久久| 成人av一区二区三区在线看| 又大又爽又粗| 777久久人妻少妇嫩草av网站| 午夜福利乱码中文字幕| 国产精品亚洲av一区麻豆| 老司机靠b影院| 中亚洲国语对白在线视频| 9191精品国产免费久久| 大香蕉久久网| 日本av手机在线免费观看| 啦啦啦免费观看视频1| 99热国产这里只有精品6| 岛国在线观看网站| 久久中文字幕一级| 真人做人爱边吃奶动态| 人妻 亚洲 视频| 两人在一起打扑克的视频| 亚洲精品国产区一区二| 久久久久久久久久久久大奶| 午夜福利影视在线免费观看| 国产一区二区三区综合在线观看| 日韩中文字幕视频在线看片| 999精品在线视频| 一级片免费观看大全| 中文字幕av电影在线播放| 精品一品国产午夜福利视频| 欧美黄色淫秽网站| 国产精品影院久久| 麻豆国产av国片精品| 曰老女人黄片| 国产精品一区二区在线观看99| 国产精品 欧美亚洲| 丝瓜视频免费看黄片| 丝袜人妻中文字幕| 9热在线视频观看99| 90打野战视频偷拍视频| 丰满迷人的少妇在线观看| 午夜福利欧美成人| 亚洲国产欧美在线一区| 伊人久久大香线蕉亚洲五| 18禁美女被吸乳视频| 国产99久久九九免费精品| 另类亚洲欧美激情| 最近最新中文字幕大全电影3 | 国产精品久久久av美女十八| 操美女的视频在线观看| 黑人巨大精品欧美一区二区蜜桃| 亚洲性夜色夜夜综合| 黄频高清免费视频| 国产成人av激情在线播放| 丁香欧美五月| 十分钟在线观看高清视频www| 国产精品国产高清国产av | 极品人妻少妇av视频| 久久久精品区二区三区| 大型av网站在线播放| 日韩大码丰满熟妇| 欧美日韩精品网址| 人成视频在线观看免费观看| 欧美激情久久久久久爽电影 | 757午夜福利合集在线观看| 性高湖久久久久久久久免费观看| 淫妇啪啪啪对白视频| av欧美777| 桃花免费在线播放| 黑人巨大精品欧美一区二区蜜桃| 男人操女人黄网站| 久久久久久免费高清国产稀缺| 男女免费视频国产| 亚洲专区字幕在线| 精品久久久久久久毛片微露脸| 国产精品免费大片| 亚洲第一av免费看| 欧美一级毛片孕妇| 国产真人三级小视频在线观看| 国产极品粉嫩免费观看在线| 亚洲国产欧美日韩在线播放| 免费观看a级毛片全部| 亚洲三区欧美一区| 男女之事视频高清在线观看| 国产aⅴ精品一区二区三区波| 热re99久久精品国产66热6| 国产一区二区三区综合在线观看| 黄片大片在线免费观看| 国产成人av激情在线播放| 精品国产亚洲在线| 黄色视频不卡| 国产区一区二久久| 欧美人与性动交α欧美软件| 欧美老熟妇乱子伦牲交| 精品国产乱码久久久久久小说| 啦啦啦视频在线资源免费观看| 日韩免费av在线播放| 18禁美女被吸乳视频| 人人澡人人妻人| 精品亚洲成a人片在线观看| 老鸭窝网址在线观看| 黄色毛片三级朝国网站| 国产真人三级小视频在线观看| 这个男人来自地球电影免费观看| 精品国产乱子伦一区二区三区| 中文字幕制服av| 人成视频在线观看免费观看| 国产在线精品亚洲第一网站| 国产在线观看jvid| 91国产中文字幕| 热re99久久国产66热| 亚洲视频免费观看视频| 中文字幕色久视频| 中文字幕制服av| 亚洲精华国产精华精| 亚洲一区中文字幕在线| 日本a在线网址| 亚洲自偷自拍图片 自拍| 精品久久久久久久毛片微露脸| 高潮久久久久久久久久久不卡| 欧美变态另类bdsm刘玥| 搡老乐熟女国产| 他把我摸到了高潮在线观看 | 国产精品二区激情视频| 人人妻人人添人人爽欧美一区卜| a级片在线免费高清观看视频| 777米奇影视久久| 国产在线观看jvid| 国产精品免费大片| 国产免费av片在线观看野外av| 亚洲综合色网址| 18禁观看日本| 日本撒尿小便嘘嘘汇集6| 女人高潮潮喷娇喘18禁视频| e午夜精品久久久久久久| 国产日韩欧美亚洲二区| 最新在线观看一区二区三区| 欧美在线黄色| 国产成人啪精品午夜网站| 午夜福利欧美成人| av一本久久久久| 亚洲精品一二三| 久久国产精品人妻蜜桃| 国产欧美日韩精品亚洲av| 精品国产乱子伦一区二区三区| 黄片小视频在线播放| 最近最新中文字幕大全电影3 | 国产成人精品在线电影| 成人永久免费在线观看视频 | 777久久人妻少妇嫩草av网站| 狠狠狠狠99中文字幕| 日韩免费高清中文字幕av| 极品教师在线免费播放| 五月天丁香电影| 99国产精品一区二区蜜桃av | 一区福利在线观看| www.熟女人妻精品国产| 99国产精品免费福利视频| 日本av手机在线免费观看| 国产欧美亚洲国产| 王馨瑶露胸无遮挡在线观看| 亚洲欧美一区二区三区黑人| 国产精品99久久99久久久不卡| 亚洲国产精品一区二区三区在线| 国产精品久久久久成人av| 亚洲av电影在线进入| 国产成人精品久久二区二区免费| 国产男女超爽视频在线观看| 不卡av一区二区三区| 高清黄色对白视频在线免费看| 热99国产精品久久久久久7| 亚洲av日韩在线播放| 国产老妇伦熟女老妇高清| 一级片免费观看大全| 国产视频一区二区在线看| 中文欧美无线码| 老熟妇仑乱视频hdxx| 午夜老司机福利片| 精品亚洲成国产av| √禁漫天堂资源中文www| 国产精品久久久人人做人人爽| 久久久精品94久久精品| 99re6热这里在线精品视频| av国产精品久久久久影院| 国产在线免费精品| 又黄又粗又硬又大视频| 亚洲黑人精品在线| 一二三四在线观看免费中文在| 欧美精品人与动牲交sv欧美| 亚洲熟女毛片儿| 亚洲熟女精品中文字幕| 国产在线一区二区三区精| 欧美黄色淫秽网站| 久久精品国产亚洲av香蕉五月 | 午夜两性在线视频| 99热国产这里只有精品6| 午夜免费成人在线视频| 777米奇影视久久| 欧美变态另类bdsm刘玥| 脱女人内裤的视频| 99国产综合亚洲精品| 国产老妇伦熟女老妇高清| 视频在线观看一区二区三区| 国产无遮挡羞羞视频在线观看| 男人操女人黄网站| 免费看a级黄色片| 久久久久国产一级毛片高清牌| 国产xxxxx性猛交| 伦理电影免费视频| 亚洲久久久国产精品| 男男h啪啪无遮挡| 黄频高清免费视频| 中文亚洲av片在线观看爽 | 日日夜夜操网爽| 成在线人永久免费视频| 老司机亚洲免费影院| 亚洲一区二区三区欧美精品| 成人av一区二区三区在线看| 18禁裸乳无遮挡动漫免费视频| 欧美成狂野欧美在线观看| 天堂俺去俺来也www色官网| 欧美日韩一级在线毛片| 欧美av亚洲av综合av国产av| 成人影院久久| 久久久国产精品麻豆| 久久精品亚洲熟妇少妇任你| 天天操日日干夜夜撸| 九色亚洲精品在线播放| 男人操女人黄网站| 免费黄频网站在线观看国产| 国产高清videossex| 啦啦啦 在线观看视频| 国产在视频线精品| 又大又爽又粗| 日韩制服丝袜自拍偷拍| 美女扒开内裤让男人捅视频| 日本精品一区二区三区蜜桃| 一区二区三区乱码不卡18| av超薄肉色丝袜交足视频| 亚洲精品中文字幕在线视频| 欧美黄色淫秽网站| 欧美老熟妇乱子伦牲交| 午夜日韩欧美国产| 夜夜爽天天搞| 每晚都被弄得嗷嗷叫到高潮| 久久久久视频综合| 亚洲av美国av| 99九九在线精品视频| 欧美精品人与动牲交sv欧美| 欧美 日韩 精品 国产| 在线av久久热| 久久中文看片网| 国产成人精品久久二区二区免费| 亚洲七黄色美女视频| 黄网站色视频无遮挡免费观看| 免费观看a级毛片全部| 午夜福利视频在线观看免费| 久久ye,这里只有精品| 五月开心婷婷网| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲精品成人av观看孕妇| 美女主播在线视频| 日韩视频在线欧美| 久久99热这里只频精品6学生| 90打野战视频偷拍视频| 777久久人妻少妇嫩草av网站| 一区二区三区精品91| 久久精品国产a三级三级三级| 成年女人毛片免费观看观看9 | h视频一区二区三区| 国产精品 欧美亚洲| 黄片小视频在线播放| 9色porny在线观看| 精品国产国语对白av| 国产亚洲欧美在线一区二区| 国产免费视频播放在线视频| 国产伦理片在线播放av一区| 熟女少妇亚洲综合色aaa.| 午夜福利乱码中文字幕| 亚洲欧洲日产国产| 精品人妻在线不人妻| 久久久久视频综合| 亚洲欧美色中文字幕在线| 一本色道久久久久久精品综合| 亚洲 国产 在线| 国产国语露脸激情在线看| 天天躁夜夜躁狠狠躁躁| 一进一出抽搐动态| 亚洲av成人不卡在线观看播放网| 老鸭窝网址在线观看| 久久国产精品影院| 久久免费观看电影| 国产视频一区二区在线看| 久久性视频一级片| 久久久精品区二区三区| 香蕉国产在线看| 97人妻天天添夜夜摸| 免费看十八禁软件| 中文亚洲av片在线观看爽 | 国产精品 国内视频| 一边摸一边抽搐一进一出视频| 性色av乱码一区二区三区2| 99国产精品一区二区三区| 在线看a的网站| 老熟妇乱子伦视频在线观看| 亚洲国产欧美一区二区综合| 亚洲五月色婷婷综合| 这个男人来自地球电影免费观看| 久久国产精品影院| 青草久久国产| 亚洲少妇的诱惑av| 国产精品一区二区免费欧美| 香蕉丝袜av| 精品国产超薄肉色丝袜足j| 亚洲成国产人片在线观看| 午夜福利,免费看| 欧美人与性动交α欧美软件| 国产欧美日韩一区二区精品| 国产淫语在线视频| 99国产精品99久久久久| 男女边摸边吃奶| 老熟女久久久| 久久国产精品人妻蜜桃| 国产欧美日韩精品亚洲av| 丝袜美足系列| 日韩大码丰满熟妇| 少妇精品久久久久久久| 国产精品电影一区二区三区 | 久久久国产一区二区| 国产午夜精品久久久久久| 97人妻天天添夜夜摸| 精品久久久久久电影网| 亚洲第一青青草原| 国产精品熟女久久久久浪| 欧美日韩福利视频一区二区| 久久精品国产亚洲av香蕉五月 | 国产在线一区二区三区精| 国产高清国产精品国产三级| 免费观看av网站的网址| 欧美日韩亚洲高清精品| 一级毛片电影观看| 91av网站免费观看| 欧美 日韩 精品 国产| 国产精品欧美亚洲77777| 亚洲性夜色夜夜综合| 另类精品久久| 国产精品免费大片| 国产男女超爽视频在线观看| 热re99久久精品国产66热6| 久久久国产一区二区| 90打野战视频偷拍视频| 丰满饥渴人妻一区二区三| 老司机午夜福利在线观看视频 | 人妻久久中文字幕网| 女人久久www免费人成看片| 国产精品免费视频内射| 成人国产av品久久久| 91老司机精品| 精品福利永久在线观看| 十八禁高潮呻吟视频| 黑人欧美特级aaaaaa片| 国产成人精品无人区| 男人舔女人的私密视频| 在线播放国产精品三级| 中文字幕另类日韩欧美亚洲嫩草| a级毛片黄视频| 性少妇av在线| 国产欧美日韩精品亚洲av| 一二三四在线观看免费中文在| 精品国产一区二区久久| av视频免费观看在线观看| 日韩免费高清中文字幕av| 欧美日韩福利视频一区二区| 1024香蕉在线观看| 岛国在线观看网站| 日本vs欧美在线观看视频| 捣出白浆h1v1| 少妇 在线观看| 亚洲国产欧美在线一区| 午夜福利在线免费观看网站|