• <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ì)的灰度圖彩色化算法研究
    97超级碰碰碰精品色视频在线观看| 香蕉丝袜av| 无限看片的www在线观看| 听说在线观看完整版免费高清| bbb黄色大片| 露出奶头的视频| 麻豆一二三区av精品| 亚洲精品中文字幕一二三四区| 精品电影一区二区在线| 日韩欧美免费精品| 国产三级在线视频| 中文字幕最新亚洲高清| 1024香蕉在线观看| 淫妇啪啪啪对白视频| 嫁个100分男人电影在线观看| 18禁黄网站禁片免费观看直播| 麻豆国产av国片精品| 欧美+亚洲+日韩+国产| 丁香六月欧美| 亚洲国产精品成人综合色| 最新美女视频免费是黄的| 岛国在线免费视频观看| 中文字幕人成人乱码亚洲影| 国产免费男女视频| 亚洲成a人片在线一区二区| 一二三四社区在线视频社区8| 国产熟女xx| 亚洲国产精品成人综合色| 国产三级黄色录像| 亚洲 欧美 日韩 在线 免费| 色综合欧美亚洲国产小说| 午夜亚洲福利在线播放| 国产亚洲精品久久久久久毛片| 狠狠狠狠99中文字幕| 欧美在线黄色| 久久精品国产99精品国产亚洲性色| av中文乱码字幕在线| 小说图片视频综合网站| 高清在线国产一区| av女优亚洲男人天堂 | 一个人免费在线观看电影 | 久久久国产精品麻豆| netflix在线观看网站| 午夜亚洲福利在线播放| 18禁黄网站禁片免费观看直播| 人妻丰满熟妇av一区二区三区| 很黄的视频免费| 欧美在线一区亚洲| 成人av一区二区三区在线看| 国产1区2区3区精品| 好男人在线观看高清免费视频| 国产精品免费一区二区三区在线| 三级国产精品欧美在线观看 | 夜夜爽天天搞| 亚洲专区中文字幕在线| 一本久久中文字幕| 村上凉子中文字幕在线| 免费看光身美女| 亚洲国产中文字幕在线视频| 麻豆av在线久日| 色综合亚洲欧美另类图片| 黑人操中国人逼视频| 亚洲av熟女| 国产视频一区二区在线看| 丰满人妻熟妇乱又伦精品不卡| 成人午夜高清在线视频| 国产精品99久久99久久久不卡| 午夜成年电影在线免费观看| 琪琪午夜伦伦电影理论片6080| 美女午夜性视频免费| 1024手机看黄色片| 成人av一区二区三区在线看| 91老司机精品| 特大巨黑吊av在线直播| 国产精品精品国产色婷婷| 午夜福利18| 日本黄大片高清| 久99久视频精品免费| 三级男女做爰猛烈吃奶摸视频| 欧美成人一区二区免费高清观看 | 国内精品美女久久久久久| 久久午夜综合久久蜜桃| 嫩草影院入口| 日韩欧美三级三区| 久久欧美精品欧美久久欧美| 看免费av毛片| 好男人在线观看高清免费视频| 97碰自拍视频| 国产人伦9x9x在线观看| 美女高潮喷水抽搐中文字幕| 熟女电影av网| 精品久久久久久成人av| 国产成人系列免费观看| 亚洲国产中文字幕在线视频| 男人的好看免费观看在线视频| 国产av在哪里看| bbb黄色大片| 日韩av在线大香蕉| 99精品欧美一区二区三区四区| 男人舔女人的私密视频| 欧美日韩瑟瑟在线播放| 老司机午夜福利在线观看视频| 中文字幕人成人乱码亚洲影| 国产99白浆流出| 国产高清三级在线| 亚洲欧美激情综合另类| 亚洲aⅴ乱码一区二区在线播放| 丰满人妻一区二区三区视频av | 久久久久精品国产欧美久久久| 天堂网av新在线| 非洲黑人性xxxx精品又粗又长| 国产精品精品国产色婷婷| 男女做爰动态图高潮gif福利片| 欧美日韩黄片免| 日本熟妇午夜| 成年免费大片在线观看| 欧美黑人巨大hd| 午夜免费成人在线视频| www日本在线高清视频| 久99久视频精品免费| 男女做爰动态图高潮gif福利片| 国产伦精品一区二区三区四那| 日本三级黄在线观看| 国产精品,欧美在线| 少妇熟女aⅴ在线视频| 国内精品久久久久精免费| 黄色片一级片一级黄色片| 亚洲中文字幕一区二区三区有码在线看 | av在线蜜桃| 国产乱人伦免费视频| 欧美中文日本在线观看视频| 19禁男女啪啪无遮挡网站| 麻豆成人av在线观看| 免费在线观看日本一区| 国产人伦9x9x在线观看| 久久性视频一级片| 两个人视频免费观看高清| 一区二区三区高清视频在线| 国产亚洲精品综合一区在线观看| 人妻夜夜爽99麻豆av| 亚洲片人在线观看| 黄色 视频免费看| 欧美乱码精品一区二区三区| 久久亚洲真实| 久久精品综合一区二区三区| 欧美黄色片欧美黄色片| 一级a爱片免费观看的视频| 色老头精品视频在线观看| 国产探花在线观看一区二区| 老司机福利观看| 男插女下体视频免费在线播放| 国产伦精品一区二区三区视频9 | 欧美一级a爱片免费观看看| 成人鲁丝片一二三区免费| 成人特级av手机在线观看| 88av欧美| 久久久成人免费电影| 中文字幕av在线有码专区| 成人性生交大片免费视频hd| 国产亚洲精品久久久com| 又大又爽又粗| 国产淫片久久久久久久久 | 国产精品98久久久久久宅男小说| 日韩欧美在线二视频| 欧美成狂野欧美在线观看| 无限看片的www在线观看| 成人亚洲精品av一区二区| 亚洲熟妇中文字幕五十中出| 亚洲国产精品sss在线观看| 国产高清videossex| 中文字幕熟女人妻在线| 国产精品免费一区二区三区在线| 波多野结衣高清无吗| 在线观看免费视频日本深夜| 嫩草影院入口| 亚洲熟妇熟女久久| 国产亚洲av高清不卡| 国产亚洲欧美在线一区二区| 午夜福利高清视频| av天堂中文字幕网| 99国产精品99久久久久| 国产高清三级在线| 亚洲人成伊人成综合网2020| 欧美色视频一区免费| 欧美性猛交╳xxx乱大交人| 97超视频在线观看视频| 波多野结衣高清无吗| 少妇熟女aⅴ在线视频| 热99re8久久精品国产| 日日干狠狠操夜夜爽| 成人特级av手机在线观看| 波多野结衣高清作品| 国产欧美日韩一区二区三| 国内少妇人妻偷人精品xxx网站 | 亚洲中文字幕日韩| 在线十欧美十亚洲十日本专区| 看黄色毛片网站| 99精品在免费线老司机午夜| 人妻久久中文字幕网| 亚洲欧美激情综合另类| 别揉我奶头~嗯~啊~动态视频| 成人特级av手机在线观看| 亚洲美女视频黄频| 久久久久久久久久黄片| 国产三级中文精品| 一进一出抽搐gif免费好疼| 国产视频一区二区在线看| 身体一侧抽搐| 亚洲一区二区三区色噜噜| 又紧又爽又黄一区二区| 老熟妇乱子伦视频在线观看| 在线a可以看的网站| 欧美不卡视频在线免费观看| 久久午夜综合久久蜜桃| 这个男人来自地球电影免费观看| 久久久久久大精品| 国产亚洲av嫩草精品影院| 18禁黄网站禁片免费观看直播| 夜夜看夜夜爽夜夜摸| 日本黄色视频三级网站网址| 麻豆久久精品国产亚洲av| 每晚都被弄得嗷嗷叫到高潮| 亚洲午夜精品一区,二区,三区| 白带黄色成豆腐渣| 一个人看视频在线观看www免费 | 久久久久久久午夜电影| 国产毛片a区久久久久| 亚洲一区高清亚洲精品| 日本黄色片子视频| 亚洲成人久久爱视频| 18禁裸乳无遮挡免费网站照片| 欧美在线一区亚洲| 性色avwww在线观看| 成年人黄色毛片网站| 嫩草影院入口| 午夜福利视频1000在线观看| 午夜亚洲福利在线播放| 精品国产乱码久久久久久男人| 国产精品爽爽va在线观看网站| 久久精品影院6| 首页视频小说图片口味搜索| 国产久久久一区二区三区| 精品久久久久久成人av| 亚洲av电影不卡..在线观看| 一卡2卡三卡四卡精品乱码亚洲| 国产成人欧美在线观看| 久久午夜亚洲精品久久| 成人无遮挡网站| 亚洲欧美日韩卡通动漫| 69av精品久久久久久| 欧美在线一区亚洲| 99riav亚洲国产免费| 免费看日本二区| 精品久久久久久成人av| 欧洲精品卡2卡3卡4卡5卡区| svipshipincom国产片| 午夜久久久久精精品| 亚洲天堂国产精品一区在线| 中文字幕最新亚洲高清| 亚洲欧洲精品一区二区精品久久久| 精品一区二区三区av网在线观看| 18禁黄网站禁片午夜丰满| 国产成人一区二区三区免费视频网站| 免费在线观看亚洲国产| 日本三级黄在线观看| 亚洲欧美精品综合久久99| 欧美黄色淫秽网站| 国产一区二区在线av高清观看| 久久精品综合一区二区三区| 中亚洲国语对白在线视频| 丝袜人妻中文字幕| tocl精华| aaaaa片日本免费| 欧美3d第一页| 欧美国产日韩亚洲一区| 亚洲精品456在线播放app | 日日摸夜夜添夜夜添小说| 国产精品亚洲一级av第二区| 国产精品综合久久久久久久免费| 国产主播在线观看一区二区| 亚洲中文字幕日韩| 久久精品人妻少妇| 丰满人妻一区二区三区视频av | 国产亚洲精品久久久久久毛片| av在线蜜桃| 国产又色又爽无遮挡免费看| 香蕉丝袜av| 黄色片一级片一级黄色片| 久久久久久久久免费视频了| 国产av一区在线观看免费| 欧美中文综合在线视频| 日本一本二区三区精品| 亚洲国产日韩欧美精品在线观看 | 亚洲国产日韩欧美精品在线观看 | 中文字幕人妻丝袜一区二区| 亚洲精品美女久久久久99蜜臀| 变态另类丝袜制服| 搡老妇女老女人老熟妇| 国产成人精品久久二区二区91| 最近最新免费中文字幕在线| 午夜精品在线福利| 很黄的视频免费| 精品久久蜜臀av无| 成人三级做爰电影| 两个人看的免费小视频| 亚洲中文日韩欧美视频| 欧美精品啪啪一区二区三区| 小说图片视频综合网站| 久久久久久九九精品二区国产| 精品久久久久久,| 一a级毛片在线观看| 深夜精品福利| a级毛片a级免费在线| 观看美女的网站| 日本一二三区视频观看| 亚洲人成伊人成综合网2020| 激情在线观看视频在线高清| 制服丝袜大香蕉在线| 国产精品爽爽va在线观看网站| 国产野战对白在线观看| 婷婷精品国产亚洲av在线| 少妇丰满av| 天堂√8在线中文| 18禁美女被吸乳视频| 亚洲欧美一区二区三区黑人| 男人和女人高潮做爰伦理| 色av中文字幕| 久久久久久久久免费视频了| 91在线精品国自产拍蜜月 | 欧美极品一区二区三区四区| 色老头精品视频在线观看| 人人妻人人看人人澡| 免费看光身美女| 免费av不卡在线播放| 神马国产精品三级电影在线观看| 热99在线观看视频| 国产蜜桃级精品一区二区三区| 舔av片在线| 午夜精品久久久久久毛片777| 欧美一级毛片孕妇| 小蜜桃在线观看免费完整版高清| 亚洲国产精品合色在线| 一进一出抽搐gif免费好疼| 亚洲av成人av| 免费搜索国产男女视频| 一本久久中文字幕| 夜夜夜夜夜久久久久| 免费观看精品视频网站| 国产精品98久久久久久宅男小说| 成年女人毛片免费观看观看9| 亚洲欧美激情综合另类| 男女那种视频在线观看| 一级毛片精品| 亚洲av五月六月丁香网| 神马国产精品三级电影在线观看| 香蕉av资源在线| 51午夜福利影视在线观看| 一级作爱视频免费观看| 国产人伦9x9x在线观看| 免费在线观看视频国产中文字幕亚洲| 国产av麻豆久久久久久久| 亚洲最大成人中文| 国产高清视频在线播放一区| 久久婷婷人人爽人人干人人爱| 久久精品综合一区二区三区| 国产精品久久视频播放| 国产视频内射| 亚洲在线自拍视频| 婷婷亚洲欧美| 国产精品国产高清国产av| 每晚都被弄得嗷嗷叫到高潮| 久久婷婷人人爽人人干人人爱| 亚洲av日韩精品久久久久久密| 欧美xxxx黑人xx丫x性爽| 无人区码免费观看不卡| 国产综合懂色| 亚洲国产精品合色在线| 日韩精品青青久久久久久| 欧美黄色片欧美黄色片| 久久亚洲真实| av片东京热男人的天堂| 亚洲欧美日韩高清在线视频| 国语自产精品视频在线第100页| 日韩欧美 国产精品| 黄色女人牲交| 黄频高清免费视频| 亚洲人成电影免费在线| 日本 av在线| 精品久久久久久久久久免费视频| 亚洲欧美精品综合久久99| 久久久久久久久免费视频了| 在线观看舔阴道视频| 国产99白浆流出| 18禁国产床啪视频网站| 制服丝袜大香蕉在线| 男人和女人高潮做爰伦理| 亚洲第一电影网av| 国产人伦9x9x在线观看| 特级一级黄色大片| av福利片在线观看| 国内精品美女久久久久久| 嫩草影视91久久| 久久精品国产亚洲av香蕉五月| а√天堂www在线а√下载| 真人一进一出gif抽搐免费| 婷婷丁香在线五月| 成人高潮视频无遮挡免费网站| 亚洲中文字幕一区二区三区有码在线看 | 真人做人爱边吃奶动态| 人妻夜夜爽99麻豆av| 日韩免费av在线播放| 亚洲熟妇熟女久久| 亚洲av第一区精品v没综合| 日本撒尿小便嘘嘘汇集6| 国产精品野战在线观看| 一级黄色大片毛片| 热99在线观看视频| 欧美xxxx黑人xx丫x性爽| x7x7x7水蜜桃| 人人妻人人看人人澡| 久久欧美精品欧美久久欧美| 又粗又爽又猛毛片免费看| 男女做爰动态图高潮gif福利片| 久久香蕉国产精品| 欧美成狂野欧美在线观看| 欧美另类亚洲清纯唯美| 白带黄色成豆腐渣| 久久香蕉国产精品| 国产精品久久久久久人妻精品电影| 美女扒开内裤让男人捅视频| 精品不卡国产一区二区三区| 五月玫瑰六月丁香| 亚洲成a人片在线一区二区| 国产精品一区二区精品视频观看| 97超视频在线观看视频| 18禁黄网站禁片免费观看直播| 叶爱在线成人免费视频播放| 午夜激情欧美在线| 久久精品亚洲精品国产色婷小说| 午夜福利欧美成人| 美女cb高潮喷水在线观看 | 亚洲熟妇中文字幕五十中出| 久久午夜亚洲精品久久| 国产伦一二天堂av在线观看| 亚洲精品美女久久av网站| 精品无人区乱码1区二区| 网址你懂的国产日韩在线| 九色成人免费人妻av| 色尼玛亚洲综合影院| 老熟妇仑乱视频hdxx| 偷拍熟女少妇极品色| 国产 一区 欧美 日韩| 国产男靠女视频免费网站| 久久亚洲真实| 窝窝影院91人妻| 色在线成人网| 精品一区二区三区视频在线 | 最近最新中文字幕大全电影3| 国内少妇人妻偷人精品xxx网站 | 少妇的逼水好多| 一进一出好大好爽视频| 在线观看午夜福利视频| 免费观看精品视频网站| 日韩成人在线观看一区二区三区| 亚洲欧美日韩卡通动漫| 国产高清三级在线| 99久久久亚洲精品蜜臀av| 国产精品久久久久久人妻精品电影| 日日夜夜操网爽| 国产不卡一卡二| 热99re8久久精品国产| 精品福利观看| 亚洲熟妇中文字幕五十中出| 精品日产1卡2卡| 亚洲 国产 在线| www.www免费av| 国产精品99久久久久久久久| 99热只有精品国产| 非洲黑人性xxxx精品又粗又长| 亚洲精品一卡2卡三卡4卡5卡| 黄色视频,在线免费观看| 日本一本二区三区精品| 午夜福利免费观看在线| 99国产极品粉嫩在线观看| cao死你这个sao货| 少妇的逼水好多| 色吧在线观看| 人妻久久中文字幕网| 神马国产精品三级电影在线观看| 精品久久久久久久末码| 国产精品电影一区二区三区| 日韩成人在线观看一区二区三区| 中国美女看黄片| 国产伦在线观看视频一区| 天堂√8在线中文| 国产精品久久电影中文字幕| 国产野战对白在线观看| 99在线视频只有这里精品首页| 国产精品乱码一区二三区的特点| 国产av不卡久久| 欧美日韩亚洲国产一区二区在线观看| 97碰自拍视频| 成年人黄色毛片网站| 精品福利观看| 国产精品女同一区二区软件 | 国产精品亚洲av一区麻豆| 999久久久精品免费观看国产| 国产伦人伦偷精品视频| 久久人妻av系列| 一区福利在线观看| 国产97色在线日韩免费| 成人特级黄色片久久久久久久| 香蕉丝袜av| 国产欧美日韩精品一区二区| 男人舔女人下体高潮全视频| 国产精品美女特级片免费视频播放器 | 亚洲中文字幕日韩| 国产高清有码在线观看视频| 亚洲av中文字字幕乱码综合| 18禁国产床啪视频网站| 中出人妻视频一区二区| 亚洲 国产 在线| 午夜福利视频1000在线观看| 欧美激情久久久久久爽电影| 九九热线精品视视频播放| 九色国产91popny在线| 在线观看一区二区三区| aaaaa片日本免费| 人妻夜夜爽99麻豆av| 麻豆av在线久日| 不卡av一区二区三区| 90打野战视频偷拍视频| 欧美一级a爱片免费观看看| 最新中文字幕久久久久 | 国产探花在线观看一区二区| 丰满的人妻完整版| 精品午夜福利视频在线观看一区| 日韩欧美国产一区二区入口| 1024手机看黄色片| 国产激情偷乱视频一区二区| 男女做爰动态图高潮gif福利片| 最新在线观看一区二区三区| av在线蜜桃| 一本一本综合久久| 亚洲熟女毛片儿| 欧美zozozo另类| 国产伦一二天堂av在线观看| 国产精品免费一区二区三区在线| 国产精品亚洲av一区麻豆| 亚洲熟女毛片儿| 成人午夜高清在线视频| 搡老熟女国产l中国老女人| 九九热线精品视视频播放| 日日夜夜操网爽| АⅤ资源中文在线天堂| 久久久久久久久免费视频了| 亚洲美女黄片视频| 99热6这里只有精品| 三级男女做爰猛烈吃奶摸视频| 性色av乱码一区二区三区2| 搡老熟女国产l中国老女人| 亚洲av电影在线进入| 真人做人爱边吃奶动态| 女生性感内裤真人,穿戴方法视频| 亚洲人成网站高清观看| 久久久久久久久久黄片| 国产成人一区二区三区免费视频网站| 在线观看免费视频日本深夜| 色播亚洲综合网| 757午夜福利合集在线观看| 国产精品九九99| 国产激情欧美一区二区| 19禁男女啪啪无遮挡网站| 国产成人av激情在线播放| 天堂动漫精品| 久久久精品大字幕| 亚洲精品美女久久久久99蜜臀| 蜜桃久久精品国产亚洲av| 在线观看一区二区三区| 日韩欧美在线二视频| 美女高潮喷水抽搐中文字幕| 一级毛片精品| 久久香蕉国产精品| 欧美乱色亚洲激情| 最近最新中文字幕大全免费视频| 亚洲色图av天堂| 最近最新免费中文字幕在线| 又紧又爽又黄一区二区| 欧美成人一区二区免费高清观看 | 亚洲熟妇中文字幕五十中出| 国内精品久久久久精免费| 天堂√8在线中文| 丁香六月欧美| 久久久久久久午夜电影| 久久久久久久久中文| 他把我摸到了高潮在线观看| 国产精品永久免费网站| 最好的美女福利视频网| 99热6这里只有精品| 欧美+亚洲+日韩+国产| 一级作爱视频免费观看| 成人av一区二区三区在线看| 在线观看免费午夜福利视频| 成人国产一区最新在线观看| 很黄的视频免费| 久久国产乱子伦精品免费另类| 成人鲁丝片一二三区免费| 美女免费视频网站| 亚洲国产精品久久男人天堂| 国产蜜桃级精品一区二区三区| 国产av一区在线观看免费| 精品乱码久久久久久99久播| 国产精品一区二区免费欧美| 99在线人妻在线中文字幕| 中文字幕高清在线视频|