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

    鰻弧菌(Vibrio anguillarum)核酸適配體篩選中高頻序列的進(jìn)化多樣性研究*

    2021-10-11 09:53:18劉慧敏林筱鈞江興龍鄢慶枇范云庭湯學(xué)敏
    海洋與湖沼 2021年5期
    關(guān)鍵詞:優(yōu)勢(shì)

    許 凈 劉慧敏 林筱鈞 鄭 江,,3① 江興龍,鄢慶枇 范云庭 湯學(xué)敏

    (1. 集美大學(xué)水產(chǎn)學(xué)院 廈門 361021; 2. 鰻鱺現(xiàn)代產(chǎn)業(yè)技術(shù)教育部工程研究中心 廈門 361021; 3. 福建省水產(chǎn)生物育種與健康養(yǎng)殖工程研究中心 廈門 361021; 4. 福建省特種水產(chǎn)配合飼料重點(diǎn)實(shí)驗(yàn)室 福清 350308)

    核酸適配體是通過(guò)指數(shù)富集配體的系統(tǒng)進(jìn)化技術(shù)(Systematic Evolution of Ligands by Exponential Enrichment), 即SELEX 篩選技術(shù), 獲得的對(duì)靶目標(biāo)有較好親和特異性的寡核苷酸分子(Ellingtonet al,1990)。核酸適配體具有靶目標(biāo)范圍廣、體積小、易合成、親和力高、特異性強(qiáng)等優(yōu)點(diǎn), 在癌癥治療、微生物檢測(cè)、食品安全等領(lǐng)域中呈現(xiàn)出廣闊的應(yīng)用前景(Hamulaet al, 2011; 于寒松等, 2015; Yazdian-Robatiet al, 2017; Liuet al, 2020)。SELEX 篩選中, 需要對(duì)獲得的篩選產(chǎn)物進(jìn)行測(cè)序, 然后從測(cè)序結(jié)果中選擇部分序列進(jìn)行親和特異性的驗(yàn)證, 親和特異性較好的序列則被確認(rèn)為核酸適配體序列。然而, 通過(guò)傳統(tǒng)克隆測(cè)序獲得的序列有成百上千條, 通過(guò)高通量測(cè)序獲得的序列更是高達(dá)數(shù)萬(wàn)條, 因此, 如何從如此眾多的序列中精準(zhǔn)、高效的選擇出相應(yīng)的候選序列進(jìn)行核酸適配體的驗(yàn)證, 是一個(gè)亟待解決的問(wèn)題。

    高頻序列是測(cè)序結(jié)果中出現(xiàn)的頻率大于等于2次的序列, 研究表明, 相當(dāng)多的高頻序列都是在篩選進(jìn)化中具有競(jìng)爭(zhēng)優(yōu)勢(shì)的序列, 且大概率是對(duì)靶目標(biāo)有較好親和力的核酸適配體序列(鄭江等, 2014), 而高頻序列的數(shù)量在總序列中的占比常常很低, 因此,對(duì)高頻序列進(jìn)行分析研究, 研究其在篩選進(jìn)化過(guò)程中的優(yōu)勢(shì)序列及其多樣性變化, 將能有針對(duì)性地尋找到相應(yīng)的核酸適配體序列, 從而能夠降低篩選的盲目性, 提高篩選的效率。

    鰻弧菌(Vibrio anguillarum)是多種淡海水養(yǎng)殖品種的致病菌(饒穎竹等, 2016; 王鳳青等, 2018;Mohamadet al, 2019), 由其導(dǎo)致的疾病給水產(chǎn)養(yǎng)殖業(yè)造成巨大的經(jīng)濟(jì)損失, 而對(duì)鰻弧菌進(jìn)行快速的分析檢測(cè)則是其病害防治的前提和基礎(chǔ)。核酸適配體因其具有親和力高、特異性強(qiáng)等優(yōu)點(diǎn), 已廣泛應(yīng)用于微生物檢測(cè)(Tanget al, 2013; Zhenget al, 2015; Yuet al,2019; Sadsriet al, 2020), 因此篩選鰻弧菌的核酸適配體, 對(duì)于鰻弧菌的快速檢測(cè)及其病害防治都具有積極意義。而研究其篩選過(guò)程中高頻序列的進(jìn)化特點(diǎn)及其多樣性變化, 不僅能快速準(zhǔn)確地獲得鰻弧菌核酸適配體的候選序列, 對(duì)于理解掌握篩選文庫(kù)的進(jìn)化特點(diǎn), 提高相關(guān)水產(chǎn)病原菌核酸適配體的篩選效率都具有積極意義。

    本文選擇了13 輪鰻弧菌核酸適配體的篩選產(chǎn)物進(jìn)行高通量測(cè)序, 采用生態(tài)學(xué)中常用的豐富度指數(shù)、均勻度指數(shù)、多樣性指數(shù)和優(yōu)勢(shì)度指數(shù)對(duì)其中的高頻序列進(jìn)行了多樣性分析, 并采用相對(duì)重要性指數(shù)(IRI)對(duì)高頻序列進(jìn)行了歸納分類, 最后利用在線網(wǎng)站對(duì)IRI 較高的7 種核酸適配體的二級(jí)結(jié)構(gòu)進(jìn)行了預(yù)測(cè)模擬, 相關(guān)研究對(duì)于核酸適配體的篩選以及鰻弧菌的檢測(cè)防治都具有重要意義。

    1 材料與方法

    1.1 序列的來(lái)源及其種類、個(gè)體數(shù)的定義

    本文中的高頻序列來(lái)自于以鰻弧菌為靶目標(biāo)的SELEX 篩選, 具體流程如下: (1) 構(gòu)建并合成初始隨機(jī)寡核苷酸文庫(kù): 5′-TCAGTCGCTTCGCCGTCTCC TTC(N35)GCACAAGAG GGAGACCCCAGAGGG-3′,N35 為含35 個(gè)隨機(jī)堿基的寡核苷酸序列。(2) 孵育:將上述寡核苷酸文庫(kù)與鰻弧菌共同孵育30 min—2 h,使寡核苷酸分子與鰻弧菌結(jié)合; (3) 分離: 利用離心、洗滌、熱變性等分離技術(shù), 將未結(jié)合的或結(jié)合較弱的寡核苷酸分子與鰻弧菌分離, 獲得能與鰻弧菌較好結(jié)合的寡核苷酸分子; (4) 擴(kuò)增: 將能與鰻弧菌結(jié)合的寡核苷酸分子進(jìn)行不對(duì)稱PCR 擴(kuò)增, 此為完成一輪篩選, 該P(yáng)CR 產(chǎn)物即為該輪的篩選產(chǎn)物, 也是下一輪篩選的次級(jí)文庫(kù); (5) 重復(fù): 重復(fù)步驟(1)到(4), 即次級(jí)文庫(kù)再次與鰻弧菌進(jìn)行孵育、分離、擴(kuò)增, 獲得相應(yīng)的篩選產(chǎn)物; (6) 測(cè)序: 對(duì)篩選產(chǎn)物進(jìn)行高通量測(cè)序, 從中選擇核酸適配體的候選序列進(jìn)行后續(xù)的驗(yàn)證和鑒定。我們共進(jìn)行了15 輪的篩選, 選擇了其中的13 輪篩選產(chǎn)物, 即第1、2、3、4、5、6、7、9、11、12、13、14、15 輪的PCR 產(chǎn)物, 進(jìn)行高通量測(cè)序, 選擇測(cè)序結(jié)果中的高頻序列進(jìn)行分析研究。高通量測(cè)序由生工生物工程有限公司(上海)完成。

    為了方便理解, 我們對(duì)本文中核酸序列的種類和個(gè)體數(shù)作如下定義或表述: 不同種類的序列是指堿基的組成、數(shù)量或排列順序不同的序列, 通常被簡(jiǎn)稱為不同序列; 同一種序列在測(cè)序結(jié)果中如果出現(xiàn)多次, 它的出現(xiàn)頻率可以看作是這種序列的個(gè)體數(shù),而每輪測(cè)序結(jié)果中的總序列數(shù)則可看作該輪測(cè)序獲得的總個(gè)體數(shù)。

    1.2 高頻序列在總序列中的占比

    統(tǒng)計(jì)每輪測(cè)序結(jié)果中高頻序列的個(gè)體數(shù)和每種高頻序列出現(xiàn)的頻率(即該高頻序列的個(gè)體數(shù)), 則該輪高頻序列個(gè)體數(shù)的占比(%)=(該輪高頻序列個(gè)體數(shù)/該輪測(cè)序序列總個(gè)體數(shù))×100=(該輪所有高頻序列的頻率之和/該輪測(cè)序序列總個(gè)體數(shù))×100, 分別計(jì)算出13 個(gè)測(cè)序輪中高頻序列的占比, 比較分析其變化規(guī)律。

    1.3 高頻序列的多樣性分析

    參考王旭娜等(2018)和黃雅琴等(2020)的方法,分別選用豐富度指數(shù)、均勻度指數(shù)、多樣性指數(shù)和優(yōu)勢(shì)度指數(shù), 來(lái)研究高頻序列的多樣性, 這些指數(shù)的原計(jì)算公式如下:

    原公式中,S為所有物種的種類數(shù),N為每輪物種的個(gè)體數(shù)量的總和,Ni為第i個(gè)物種的個(gè)體數(shù),Pi=Ni/N,Nmax為優(yōu)勢(shì)種的總個(gè)體數(shù),Nt為各物種的總個(gè)體數(shù)。將上述公式應(yīng)用于本文的高頻序列, 則公式中的S就是某測(cè)序輪中高頻序列的種類數(shù),N為該輪中高頻序列的總個(gè)體數(shù),Ni是該輪第i個(gè)高頻序列的個(gè)體數(shù), 即第i個(gè)高頻序列的出現(xiàn)次數(shù)或頻率,Pi=Ni/N,Nmax為該輪中優(yōu)勢(shì)高頻序列的個(gè)體數(shù),Nt是該輪中高頻序列的總出現(xiàn)次數(shù)。因此, 根據(jù)上述公式, 采用Primer 6.0 軟件可計(jì)算出每一個(gè)測(cè)序輪中的上述四種指數(shù)。

    1.4 高頻序列的相對(duì)重要性分析

    優(yōu)勢(shì)種的研究通常是先選擇一系列的站點(diǎn), 然后計(jì)算每個(gè)站點(diǎn)中的各個(gè)物種的相對(duì)重要性指數(shù)(Index of relative importance, IRI), 再根據(jù)該指數(shù)對(duì)相應(yīng)的物種進(jìn)行判定分類(王雪輝等, 2010; 劉曉霞等,2016; 張鵬弛等, 2017; 粟麗等, 2018)。IRI 的原始計(jì)算公式為: IRI = (N+W)×F×10 000 (沃佳等, 2017; 丁朋朋等, 2019), 式中,N為某一種類的個(gè)體數(shù)占總個(gè)體數(shù)的百分率;W為某一種類的重量占總重量的百分率;F為某一種類出現(xiàn)的站數(shù)占調(diào)查總站數(shù)的百分率。在本文高頻序列的研究中, 我們將各個(gè)篩選輪或測(cè)序輪看作是原公式中的各個(gè)站點(diǎn), 則總共13 個(gè)測(cè)序輪就對(duì)應(yīng)著13 個(gè)站點(diǎn), 某種高頻序列出現(xiàn)的次數(shù)或頻率就相當(dāng)于該高頻序列的個(gè)體數(shù), 則公式中N就是某種高頻序列在13 個(gè)測(cè)序輪中出現(xiàn)的總次數(shù)占13輪中全部高頻序列的出現(xiàn)總次數(shù)的百分比,F則為該高頻序列出現(xiàn)的輪數(shù)占全部輪數(shù)(即13 輪)的百分比。另外, 由于核酸序列重量與其序列長(zhǎng)度成正比, 本文中的高頻序列的長(zhǎng)度大多在82 nt 左右, 差別不大,因此公式中的重量百分比W就近似等于其數(shù)量百分比N, 因此IRI 公式就無(wú)需加和, 可只取數(shù)量百分比N, 由此相對(duì)重要性指數(shù)IRI 的公式可修正為IRI =N×F×10 000。然后參考相關(guān)文獻(xiàn)的判定方法(陳國(guó)寶等, 2007; 李顯森等, 2013; 粟麗等, 2018), 依據(jù)IRI的值對(duì)不同序列進(jìn)行判斷分類, 具體如下: 當(dāng)該序列的IRI≥500時(shí), 該序列為優(yōu)勢(shì)序列; 100≤IRI<500時(shí),該序列為重要序列; 10≤IRI<100 時(shí), 該序列為常見(jiàn)序列, 1≤IRI<10 時(shí), 該序列為一般序列; IRI<1 時(shí), 該序列為少見(jiàn)序列。

    1.5 序列的二級(jí)結(jié)構(gòu)模擬

    利用RNA structure 網(wǎng)站(http://rna.urmc.rochester.edu/RNAstructureWeb/)進(jìn)行在線預(yù)測(cè), 核酸類型選擇為DNA, 其他參數(shù)為默認(rèn)值, 最后選擇最大可能性的二級(jí)結(jié)構(gòu), 即MaxExpect 二級(jí)結(jié)構(gòu)作為各序列的二級(jí)結(jié)構(gòu)。

    2 結(jié)果

    2.1 篩選過(guò)程中高頻序列的種類數(shù)及其在總序列中的占比變化

    高頻序列種類數(shù)的變化(圖1)與其占比的變化(圖2)相類似, 都是在第4、5 輪出現(xiàn)突然性的大幅上升,在第4 輪達(dá)到最大值, 之后快速回落后出現(xiàn)穩(wěn)步增加,在最后一輪達(dá)到次高點(diǎn)。從總趨勢(shì)上看, 高頻序列的種類數(shù)在前6 輪出現(xiàn)較劇烈的波動(dòng), 之后隨著篩選的進(jìn)行呈穩(wěn)步增加的趨勢(shì), 高頻序列的占比也呈現(xiàn)類似變化, 不過(guò)高頻序列占比的波動(dòng)要更劇烈些, 后續(xù)的增加趨勢(shì)也要更明顯。

    圖1 篩選過(guò)程中高頻序列的種類數(shù)變化Fig.1 Changes of high frequency sequences in the species during the selection

    圖2 篩選過(guò)程中高頻序列在總序列中的占比變化Fig.2 The proportion changes of high frequency sequence in the total sequence during the selection

    2.2 高頻序列的多樣性分析

    反映多樣性的Shannon-Wiener 指數(shù)和豐富度指數(shù)與高頻序列的種類數(shù)和占比變化(如圖1 和圖2)也類似, 也是在篩選的第4 輪達(dá)到最高, 之后快速回落,然后再逐步上升, 說(shuō)明隨著篩選進(jìn)行, 高頻序列也在進(jìn)化, 其豐富度和多樣性在經(jīng)歷了前6 輪的劇烈波動(dòng)后, 都隨著篩選進(jìn)化過(guò)程而逐步提高。

    Pielou 均勻度指數(shù)反映的是各物種個(gè)體數(shù)目分配的均勻程度, 指數(shù)越高, 說(shuō)明個(gè)體間數(shù)目越均勻。圖3 顯示, 該指數(shù)在前6 輪在一定區(qū)間內(nèi)呈現(xiàn)較大波動(dòng), 之后則趨于穩(wěn)定, 并呈現(xiàn)逐漸降低趨勢(shì), 說(shuō)明隨著篩選進(jìn)行, 均勻度在經(jīng)歷了前6 輪的劇烈波動(dòng)后,呈現(xiàn)越來(lái)越低的趨勢(shì), 這也進(jìn)一步說(shuō)明篩選進(jìn)化是向著某一方向進(jìn)行的, 適應(yīng)這一進(jìn)化方向的高頻序列, 其個(gè)體數(shù)就會(huì)越來(lái)越多, 不適應(yīng)這一方向的高頻序列, 其個(gè)體數(shù)就會(huì)越來(lái)越少, 從而呈現(xiàn)出越來(lái)越大的差異和不均勻性。

    優(yōu)勢(shì)度指數(shù)反映的是優(yōu)勢(shì)種的優(yōu)勢(shì)程度, 優(yōu)勢(shì)度指數(shù)越高, 說(shuō)明該輪優(yōu)勢(shì)種的優(yōu)勢(shì)越突出。圖3 顯示, 該指數(shù)在前6 輪波動(dòng)較為劇烈, 說(shuō)明優(yōu)勢(shì)種的優(yōu)勢(shì)度并不穩(wěn)定, 其他種的高頻序列可能對(duì)優(yōu)勢(shì)種呈現(xiàn)較大的競(jìng)爭(zhēng)壓力, 之后隨著篩選進(jìn)行, 優(yōu)勢(shì)度指數(shù)也逐漸走高, 說(shuō)明優(yōu)勢(shì)種在篩選進(jìn)化中的競(jìng)爭(zhēng)力越來(lái)越強(qiáng), 其在進(jìn)化競(jìng)爭(zhēng)中的優(yōu)勢(shì)越來(lái)越明顯, 這也反映了篩選效果越來(lái)越好。

    圖3 篩選過(guò)程中Margalef 豐富度指數(shù)、Shannon-Wiener 多樣性指數(shù)、Pielou 均勻度指數(shù)、Berger-Parker 優(yōu)勢(shì)度指數(shù)的變化Fig.3 Changes of the Margalef richness index, Shannon-Wiener diversity index, Pielou evenness index, and Berger-Parker dominance index during the selection

    2.3 高頻序列的相對(duì)重要性分析及其分類

    對(duì)13 輪中的79 個(gè)高頻序列進(jìn)行了相對(duì)重要性指數(shù)(IRI)的計(jì)算, 結(jié)果如表1。依據(jù)IRI 的大小可將這些高頻序列分為以下五種類型, 即優(yōu)勢(shì)序列(IRI≥500)、重要序列(100 ≤IRI<500)、常見(jiàn)序列(10 ≤IRI<100)、一般序列(1≤IRI<10)和少見(jiàn)序列(IRI<1)。從表中可以看出, 優(yōu)勢(shì)序列和重要序列分別只有1 種,約占全部高頻序列的1.27%, 常見(jiàn)序列、一般序列和少見(jiàn)序列, 則分別約占全部高頻序列的 8.86%、18.99%和69.62%。另外, 優(yōu)勢(shì)序列的IRI 值高達(dá)8 707.25, 是重要序列的19 倍多, 這表明優(yōu)勢(shì)序列在篩選進(jìn)化中的競(jìng)爭(zhēng)優(yōu)勢(shì)是較為突出的。

    通過(guò)表1 的分類可看出, 優(yōu)勢(shì)序列、重要序列和常見(jiàn)序列應(yīng)該是篩選進(jìn)化中的相對(duì)優(yōu)勢(shì)種群, 它們對(duì)靶目標(biāo)鰻弧菌的競(jìng)爭(zhēng)力大概率要高于一般序列和少見(jiàn)序列, 因此, 按照IRI 指數(shù)從高到低, 優(yōu)先從優(yōu)勢(shì)序列、重要序列中挑選相應(yīng)的序列進(jìn)行核酸適配體的驗(yàn)證, 將能大大降低驗(yàn)證的盲目性和隨機(jī)性, 提高篩選效率。后續(xù)選擇了IRI 最高的7 個(gè)高頻序列(H5、H1、H6、H12、H17、H28、H21)進(jìn)行了親和特異性驗(yàn)證, 發(fā)現(xiàn)這7 個(gè)高頻序列均是親和特異性較好的核酸適配體(相關(guān)數(shù)據(jù)另文發(fā)表)。

    表1 高頻序列的相對(duì)重要性指數(shù)及其分類Tab.1 Index of relative importance of high frequency sequences and their classification

    2.4 核酸適配體序列的二級(jí)結(jié)構(gòu)

    選擇已驗(yàn)證確認(rèn)的7 個(gè)核酸適配體(H5、H1、H6、H12、H17、H28、H21)進(jìn)行了二級(jí)結(jié)構(gòu)的模擬(圖4), 從圖中可看出, 這些序列均形成了多個(gè)大小不一的閉合環(huán)狀結(jié)構(gòu), 表2 對(duì)這些核酸適配體序列中的閉合環(huán)進(jìn)行總結(jié), 可以發(fā)現(xiàn), 每個(gè)序列都含有2—6 個(gè)閉合環(huán), 每個(gè)環(huán)由5—25 核苷酸組成, 根據(jù)每個(gè)環(huán)中的核苷酸數(shù)量可以把這些環(huán)大致分為三種:小型環(huán)(由5—8 個(gè)核苷酸構(gòu)成)、中型環(huán)(由9—13 個(gè)核苷酸構(gòu)成)、大型環(huán)(由17 個(gè)以上的核苷酸構(gòu)成)。這些環(huán)的大小和數(shù)量很可能與核酸適配體的結(jié)合位點(diǎn)相匹配。

    表2 核酸適配體序列的IRI 值及其二級(jí)結(jié)構(gòu)中環(huán)的數(shù)量和大小Tab.2 Values of index of relative importance of aptamers and the numbers and sizes of their rings in their secondary structures

    圖4 核酸適配體二級(jí)結(jié)構(gòu)Fig.4 Secondary structures of the aptamers

    3 討論

    研究鰻弧菌核酸適配體篩選中高頻序列的進(jìn)化多樣性對(duì)于揭示核酸適配體的篩選進(jìn)化規(guī)律、提高核酸適配體的篩選效率都具有重要意義。本文分別從高頻序列的種類數(shù)、多樣性、相對(duì)重要性指數(shù)(IRI)和二級(jí)結(jié)構(gòu)等四個(gè)方面對(duì)高頻序列的篩選進(jìn)化特點(diǎn)進(jìn)行了研究, 下面也從這四個(gè)方面對(duì)其進(jìn)行分析討論。

    3.1 篩選過(guò)程中高頻序列的種類數(shù)和占比變化

    理論上, 高頻序列是在篩選競(jìng)爭(zhēng)中占優(yōu)勢(shì)的序列, 其種類數(shù)量和在總序列中的占比應(yīng)該隨著篩選輪數(shù)的增加而逐步增加。但我們的研究卻發(fā)現(xiàn), 在篩選的中間階段出現(xiàn)了高頻序列大量爆增、劇烈波動(dòng)的現(xiàn)象, 其種類和占比都超過(guò)了最后一輪。出現(xiàn)這種劇烈波動(dòng)現(xiàn)象的原因很可能和篩選條件有關(guān)。前6 輪的篩選是在較為寬松的條件下進(jìn)行的(結(jié)合1—2 h, 洗滌1—2 次), 所以, 雖然第4、5 輪出現(xiàn)了高頻序列的大爆發(fā), 但大部分高頻序列的親和力可能并不高, 第7 輪后為了去除親和力較弱的序列, 篩選條件逐漸嚴(yán)苛, 結(jié)合時(shí)間縮短為0.5 h, 洗滌次數(shù)增加為3 次, 從而導(dǎo)致相當(dāng)多親和力不高的高頻序列被淘汰, 最后雖然高頻序列的數(shù)量和占比沒(méi)有達(dá)到最高, 但其中高親和力序列的數(shù)量大大增加, 文庫(kù)中高頻序列的質(zhì)量得到了較大的改善。

    3.2 篩選過(guò)程中高頻序列的進(jìn)化多樣性

    Margalef 豐富度指數(shù)、Shannon-Wiener 多樣性指數(shù)、Pielou 均勻度指數(shù)和Berger-Parker 優(yōu)勢(shì)度指數(shù),經(jīng)常被用于動(dòng)植物的多樣性研究。作為一種有效的分析和描述方法, 這些指數(shù)還被用于手足口病的流行性特征分析(Yanget al, 2020)、口腔微生物菌群的組成分析(Wolffet al, 2019)、以及土壤微生物多樣性的研究(夏圍圍等, 2014)。本文則采用這四個(gè)指數(shù)研究了核酸適配體篩選過(guò)程中高頻寡核苷酸序列的進(jìn)化多樣性, 結(jié)果表明, 隨著篩選進(jìn)行, 高頻序列的豐富度、多樣性和優(yōu)勢(shì)度都呈現(xiàn)逐步增加趨勢(shì), 而均勻度出現(xiàn)降低。這一結(jié)果很可能是因?yàn)樵诤Y選競(jìng)爭(zhēng)中, 在總序列中占比90%以上的、大量親和力較弱的非高頻序列逐漸被淘汰, 而親和力較強(qiáng)的非高頻序列則逐漸進(jìn)化為高頻序列, 從而導(dǎo)致含有大量核酸適配體的、有較高親和力的高頻序列在篩選競(jìng)爭(zhēng)中逐漸占據(jù)優(yōu)勢(shì), 其種類數(shù)量不斷增加, 競(jìng)爭(zhēng)優(yōu)勢(shì)也越來(lái)越強(qiáng),并呈現(xiàn)出強(qiáng)者越強(qiáng)的定向進(jìn)化趨勢(shì)。這些變化趨勢(shì)與核酸適配體的篩選進(jìn)化特點(diǎn)也是一致的。

    3.3 相對(duì)重要性指數(shù)的意義

    核酸適配體的篩選中, 需要對(duì)篩選產(chǎn)物進(jìn)行測(cè)序, 然后再?gòu)谋姸嗟臏y(cè)序序列中挑選出可能的候選序列, 進(jìn)行親和特異性的驗(yàn)證, 只有親和特異性較好的序列才被認(rèn)定為核酸適配體。如何挑選候選序列,目前尚沒(méi)有統(tǒng)一的方法。早期是采用隨機(jī)的方法從測(cè)序結(jié)果中隨機(jī)挑選數(shù)個(gè)序列去進(jìn)行親和特異性驗(yàn)證(李衛(wèi)濱等, 2007), 或者利用同源性將測(cè)序的序列分為數(shù)個(gè)家族, 然后從數(shù)個(gè)家族中隨機(jī)挑選序列進(jìn)行驗(yàn)證(Liuet al, 2014), 這種隨機(jī)方法要進(jìn)行大量的驗(yàn)證實(shí)驗(yàn), 工作量較大。另外, 也有學(xué)者是按照序列的頻率和占比來(lái)挑選候選序列的(Zhuet al, 2021)。而我們則綜合考慮了頻率和出現(xiàn)輪數(shù), 采用相對(duì)重要性指數(shù)(IRI)這個(gè)指標(biāo), 按照IRI 從大到小來(lái)挑選相應(yīng)的候選序列, 挑選出的7 個(gè)高頻序列后續(xù)均被確認(rèn)是親和特異性較好的核酸適配體, 說(shuō)明采用IRI 指數(shù)來(lái)挑選候選核酸適配體序列也是可行的。

    3.4 高頻序列的二級(jí)結(jié)構(gòu)

    核酸適配體通常借助其結(jié)構(gòu)中的堆積作用、形狀互補(bǔ)、靜電作用和氫鍵等相互作用, 通過(guò)頸環(huán)等結(jié)構(gòu),與靶目標(biāo)或嵌合或包被, 從而形成穩(wěn)定的復(fù)合物(Ellingtonet al, 1992; Hermannet al, 2000)。有研究發(fā)現(xiàn), 適配體二級(jí)結(jié)構(gòu)中的不同環(huán)區(qū)可以聯(lián)合參與靶標(biāo)的識(shí)別與結(jié)合(孫羽菡, 2019), 還有研究表明, 兩個(gè)頸環(huán)的適配體要比一個(gè)頸環(huán)的適配體有更好的穩(wěn)定性和親和特異性(Bhardwajet al, 2019), 我們的研究也表明二級(jí)結(jié)構(gòu)中環(huán)的數(shù)量和大小很可能是與其結(jié)合位點(diǎn)相匹配的, 進(jìn)而影響其與靶目標(biāo)的結(jié)合和穩(wěn)定性。

    4 結(jié)論

    利用Margalef 豐富度指數(shù)、Shannon-Wiener 多樣性指數(shù)、Pielou 均勻度指數(shù)和Berger-Parker 優(yōu)勢(shì)度指數(shù), 研究了鰻弧菌核酸適配體的篩選中高頻序列的進(jìn)化多樣性, 結(jié)果發(fā)現(xiàn), 在篩選進(jìn)化過(guò)程中, 高頻序列的豐富度、多樣性、優(yōu)勢(shì)度和均勻度都在前6 輪出現(xiàn)劇烈的波動(dòng), 之后隨著篩選的進(jìn)行, 豐富度、多樣性和優(yōu)勢(shì)度都呈逐步增加趨勢(shì), 而均勻度則呈現(xiàn)逐步降低的趨勢(shì), 說(shuō)明含有大量核酸適配體的高頻序列在篩選競(jìng)爭(zhēng)中逐漸占據(jù)優(yōu)勢(shì), 其種類數(shù)量不斷增加, 競(jìng)爭(zhēng)優(yōu)勢(shì)越來(lái)越強(qiáng); 利用相對(duì)重要性指數(shù)(IRI)可將高頻序列分成優(yōu)勢(shì)序列、重要序列、常見(jiàn)序列、一般序列和少見(jiàn)序列, IRI 越高的序列在篩選進(jìn)化中的競(jìng)爭(zhēng)優(yōu)勢(shì)也越大, 也越有可能成為高親和力適配體;最后, 對(duì)IRI 較高的7 個(gè)核酸適配體的二級(jí)結(jié)構(gòu)進(jìn)行了模擬分析, 發(fā)現(xiàn)它們的二級(jí)結(jié)構(gòu)中都存在有大、中、小三種不同尺寸的環(huán)狀結(jié)構(gòu), 推測(cè)這些環(huán)的大小和數(shù)量是與靶目標(biāo)上的結(jié)合位點(diǎn)相對(duì)應(yīng)的。

    猜你喜歡
    優(yōu)勢(shì)
    優(yōu)勢(shì) 等
    創(chuàng)新發(fā)揮僑務(wù)優(yōu)勢(shì) 拓展海外統(tǒng)戰(zhàn)工作
    矮的優(yōu)勢(shì)
    老父親的優(yōu)勢(shì)
    畫與話
    發(fā)揚(yáng)優(yōu)勢(shì) 有所作為
    談“五老”的五大特殊優(yōu)勢(shì)
    第二優(yōu)勢(shì)
    從優(yōu)勢(shì)到勝勢(shì)
    棋藝(2001年19期)2001-11-25 19:55:34
    從優(yōu)勢(shì)到勝勢(shì)
    棋藝(2001年23期)2001-01-06 19:08:36
    久久精品人妻少妇| 成人国产一区最新在线观看| 久久中文看片网| 狠狠狠狠99中文字幕| 国产中年淑女户外野战色| 国产又黄又爽又无遮挡在线| 亚洲国产高清在线一区二区三| 久久精品夜夜夜夜夜久久蜜豆| 我要看日韩黄色一级片| 日本撒尿小便嘘嘘汇集6| 免费观看精品视频网站| 国产精品一区二区性色av| 我的女老师完整版在线观看| 99riav亚洲国产免费| 欧美高清性xxxxhd video| 尾随美女入室| 精华霜和精华液先用哪个| 18禁在线播放成人免费| 国产亚洲欧美98| 国产亚洲精品综合一区在线观看| 亚洲精品亚洲一区二区| 男女边吃奶边做爰视频| 国产麻豆成人av免费视频| 啦啦啦观看免费观看视频高清| 熟女电影av网| 亚洲国产欧洲综合997久久,| 亚洲国产色片| 看片在线看免费视频| 麻豆久久精品国产亚洲av| 色吧在线观看| 在现免费观看毛片| 精品久久久久久久久久免费视频| 国产免费一级a男人的天堂| 国产精品免费一区二区三区在线| 国产精品久久视频播放| 国产色婷婷99| 91午夜精品亚洲一区二区三区 | 亚洲人成网站在线播| 欧美绝顶高潮抽搐喷水| 亚洲乱码一区二区免费版| 露出奶头的视频| 国产成人aa在线观看| 国产成人一区二区在线| 久久欧美精品欧美久久欧美| 乱人视频在线观看| 国产v大片淫在线免费观看| 成人永久免费在线观看视频| 欧美zozozo另类| 国产成人av教育| 欧美最黄视频在线播放免费| 99国产极品粉嫩在线观看| 综合色av麻豆| 91麻豆av在线| 99久久九九国产精品国产免费| 精华霜和精华液先用哪个| а√天堂www在线а√下载| 免费看日本二区| 老司机深夜福利视频在线观看| 午夜影院日韩av| 成人特级av手机在线观看| 欧美成人性av电影在线观看| 男人的好看免费观看在线视频| 成年女人毛片免费观看观看9| 国产伦一二天堂av在线观看| 最新中文字幕久久久久| 成人无遮挡网站| 高清在线国产一区| 麻豆成人av在线观看| 久久久久久久久大av| 一边摸一边抽搐一进一小说| 最近最新中文字幕大全电影3| 日本成人三级电影网站| 欧美黑人欧美精品刺激| 亚洲四区av| 亚洲avbb在线观看| 变态另类成人亚洲欧美熟女| 国产av不卡久久| 最近在线观看免费完整版| 日韩大尺度精品在线看网址| 97碰自拍视频| 69人妻影院| 色av中文字幕| 成人综合一区亚洲| 亚洲第一区二区三区不卡| 99热只有精品国产| 可以在线观看的亚洲视频| 小蜜桃在线观看免费完整版高清| 久久国产乱子免费精品| 色精品久久人妻99蜜桃| 免费在线观看日本一区| 亚洲精品一区av在线观看| 国产精品女同一区二区软件 | 九九久久精品国产亚洲av麻豆| 欧美成人a在线观看| 久久精品国产鲁丝片午夜精品 | 99热6这里只有精品| 91在线观看av| 久久亚洲精品不卡| 在现免费观看毛片| 小蜜桃在线观看免费完整版高清| 午夜福利欧美成人| 99久久无色码亚洲精品果冻| 日本一本二区三区精品| 九九在线视频观看精品| 深夜精品福利| 久久久久性生活片| 婷婷亚洲欧美| 色吧在线观看| 久久久国产成人免费| 成人精品一区二区免费| 国产真实伦视频高清在线观看 | 在线观看66精品国产| 丰满人妻一区二区三区视频av| 午夜福利视频1000在线观看| 亚洲精华国产精华液的使用体验 | 国产极品精品免费视频能看的| 精品午夜福利在线看| 久99久视频精品免费| 中文资源天堂在线| 草草在线视频免费看| 搡老妇女老女人老熟妇| 午夜精品在线福利| x7x7x7水蜜桃| 又黄又爽又刺激的免费视频.| 亚洲欧美日韩高清在线视频| 91在线观看av| 国产毛片a区久久久久| 欧美日韩黄片免| 日韩欧美精品v在线| 久久精品影院6| 嫩草影院入口| 韩国av在线不卡| av在线老鸭窝| 成人国产综合亚洲| 国产极品精品免费视频能看的| 国产精品福利在线免费观看| 亚洲内射少妇av| 日本欧美国产在线视频| 亚洲无线观看免费| 我的女老师完整版在线观看| 久久99热6这里只有精品| 黄色配什么色好看| 国产av不卡久久| xxxwww97欧美| 亚洲成人免费电影在线观看| 国产精品爽爽va在线观看网站| 99热6这里只有精品| 欧美日韩瑟瑟在线播放| 有码 亚洲区| 最新在线观看一区二区三区| 变态另类成人亚洲欧美熟女| 乱人视频在线观看| 乱码一卡2卡4卡精品| 精品人妻偷拍中文字幕| 一区二区三区四区激情视频 | 最新在线观看一区二区三区| 亚洲中文字幕一区二区三区有码在线看| 久久精品综合一区二区三区| 国产精品伦人一区二区| 永久网站在线| 嫩草影院新地址| 精品人妻一区二区三区麻豆 | 美女免费视频网站| 精品欧美国产一区二区三| 日韩欧美在线乱码| 婷婷六月久久综合丁香| 别揉我奶头~嗯~啊~动态视频| 在线观看舔阴道视频| 婷婷六月久久综合丁香| 人妻少妇偷人精品九色| 亚洲va日本ⅴa欧美va伊人久久| 色吧在线观看| 国产v大片淫在线免费观看| 国产精品久久久久久久电影| 少妇裸体淫交视频免费看高清| 久久精品国产鲁丝片午夜精品 | 成人亚洲精品av一区二区| 亚洲va日本ⅴa欧美va伊人久久| 美女大奶头视频| 中国美女看黄片| 国产精品野战在线观看| 黄片wwwwww| 国产亚洲av嫩草精品影院| 日日摸夜夜添夜夜添av毛片 | 窝窝影院91人妻| 国产精品无大码| 一本久久中文字幕| 午夜免费男女啪啪视频观看 | 一边摸一边抽搐一进一小说| 熟妇人妻久久中文字幕3abv| 丰满人妻一区二区三区视频av| 久久精品国产亚洲网站| 永久网站在线| 亚洲精品影视一区二区三区av| 国产伦精品一区二区三区四那| 97人妻精品一区二区三区麻豆| 嫁个100分男人电影在线观看| 久久久久久久久中文| 在线观看一区二区三区| 在线看三级毛片| 小说图片视频综合网站| 一夜夜www| 国产高清三级在线| 久久精品久久久久久噜噜老黄 | 五月伊人婷婷丁香| 可以在线观看的亚洲视频| 婷婷亚洲欧美| 偷拍熟女少妇极品色| 老女人水多毛片| 人人妻,人人澡人人爽秒播| 中国美女看黄片| 国产亚洲精品av在线| 草草在线视频免费看| 日韩 亚洲 欧美在线| 成年版毛片免费区| 极品教师在线免费播放| 一本精品99久久精品77| 波多野结衣高清无吗| 女人十人毛片免费观看3o分钟| 国产成人a区在线观看| 国产一区二区三区av在线 | 久久午夜福利片| 国产精品人妻久久久久久| 日韩欧美一区二区三区在线观看| 干丝袜人妻中文字幕| 动漫黄色视频在线观看| 69av精品久久久久久| 国产精品野战在线观看| 深夜精品福利| 久久久久精品国产欧美久久久| 人人妻,人人澡人人爽秒播| 最后的刺客免费高清国语| 尾随美女入室| 欧美日韩中文字幕国产精品一区二区三区| 久久人妻av系列| 午夜福利18| 99在线视频只有这里精品首页| 亚洲国产精品合色在线| 亚洲内射少妇av| 九九爱精品视频在线观看| 嫩草影视91久久| 一边摸一边抽搐一进一小说| 国产精品国产高清国产av| 久久久久久久午夜电影| 精品久久久久久久久久久久久| 人人妻人人澡欧美一区二区| 99精品久久久久人妻精品| 色综合色国产| 免费观看人在逋| 联通29元200g的流量卡| 99久久中文字幕三级久久日本| 成人鲁丝片一二三区免费| 悠悠久久av| 性插视频无遮挡在线免费观看| 国产亚洲精品av在线| 深夜精品福利| 欧洲精品卡2卡3卡4卡5卡区| 亚洲自拍偷在线| 别揉我奶头~嗯~啊~动态视频| 黄色日韩在线| 久久6这里有精品| 国产欧美日韩精品一区二区| 热99在线观看视频| 99久国产av精品| 欧美+亚洲+日韩+国产| 嫁个100分男人电影在线观看| 成人高潮视频无遮挡免费网站| 有码 亚洲区| 亚洲天堂国产精品一区在线| 热99在线观看视频| 精品午夜福利视频在线观看一区| av在线观看视频网站免费| 国产欧美日韩一区二区精品| 女同久久另类99精品国产91| 久久天躁狠狠躁夜夜2o2o| 伦精品一区二区三区| 久久精品影院6| 在线播放国产精品三级| 国产乱人视频| 日日干狠狠操夜夜爽| 哪里可以看免费的av片| 69av精品久久久久久| 啦啦啦啦在线视频资源| 中国美白少妇内射xxxbb| 日韩精品青青久久久久久| a级毛片a级免费在线| 久久99热6这里只有精品| 欧美不卡视频在线免费观看| 99久久中文字幕三级久久日本| 国产精品一区www在线观看 | 色综合亚洲欧美另类图片| 国产亚洲精品久久久com| a级一级毛片免费在线观看| 精品久久久久久,| 观看美女的网站| 99久久成人亚洲精品观看| 国产精品一区二区性色av| 久久久精品欧美日韩精品| 精品日产1卡2卡| 国产视频一区二区在线看| 国产亚洲91精品色在线| 两人在一起打扑克的视频| 亚洲色图av天堂| 我的女老师完整版在线观看| 午夜爱爱视频在线播放| 欧美激情国产日韩精品一区| 国产伦精品一区二区三区视频9| 美女cb高潮喷水在线观看| 亚洲午夜理论影院| 精品国产三级普通话版| 性色avwww在线观看| 久久精品国产鲁丝片午夜精品 | a在线观看视频网站| 亚洲天堂国产精品一区在线| 久久热精品热| 男女之事视频高清在线观看| 欧美国产日韩亚洲一区| 一进一出抽搐gif免费好疼| 超碰av人人做人人爽久久| 久久久久久久久久久丰满 | 亚洲狠狠婷婷综合久久图片| 亚洲av不卡在线观看| 日本黄色视频三级网站网址| 亚洲av日韩精品久久久久久密| 如何舔出高潮| 身体一侧抽搐| 18禁在线播放成人免费| 久久精品国产鲁丝片午夜精品 | 哪里可以看免费的av片| 干丝袜人妻中文字幕| 永久网站在线| or卡值多少钱| 国产精品久久电影中文字幕| 亚洲精品一卡2卡三卡4卡5卡| 亚洲av成人av| 欧美3d第一页| 中文亚洲av片在线观看爽| 听说在线观看完整版免费高清| 色av中文字幕| 麻豆久久精品国产亚洲av| 日本免费a在线| 国产大屁股一区二区在线视频| 欧美又色又爽又黄视频| 色精品久久人妻99蜜桃| 午夜免费男女啪啪视频观看 | 国产一区二区激情短视频| 麻豆av噜噜一区二区三区| 老熟妇乱子伦视频在线观看| 变态另类丝袜制服| 精品午夜福利视频在线观看一区| 91久久精品电影网| 午夜亚洲福利在线播放| 久久久久久久亚洲中文字幕| 亚洲乱码一区二区免费版| 禁无遮挡网站| 日韩欧美在线乱码| 国产精品久久久久久精品电影| 精品人妻熟女av久视频| 国产毛片a区久久久久| АⅤ资源中文在线天堂| 蜜桃亚洲精品一区二区三区| 动漫黄色视频在线观看| 校园人妻丝袜中文字幕| 成人美女网站在线观看视频| 亚洲成人中文字幕在线播放| 久久久国产成人免费| 欧美不卡视频在线免费观看| 国产国拍精品亚洲av在线观看| 精品午夜福利视频在线观看一区| 悠悠久久av| 99久久无色码亚洲精品果冻| 赤兔流量卡办理| 九九久久精品国产亚洲av麻豆| 在线免费观看不下载黄p国产 | www.www免费av| 色尼玛亚洲综合影院| 日本黄大片高清| av在线天堂中文字幕| 最近最新中文字幕大全电影3| 无遮挡黄片免费观看| 国产淫片久久久久久久久| 免费无遮挡裸体视频| 国国产精品蜜臀av免费| 九九爱精品视频在线观看| 国产三级中文精品| a在线观看视频网站| 伦精品一区二区三区| 十八禁网站免费在线| 精品日产1卡2卡| 久久精品人妻少妇| 一本久久中文字幕| 伊人久久精品亚洲午夜| 国产av不卡久久| 亚洲图色成人| 国产男靠女视频免费网站| 99久久精品国产国产毛片| 又黄又爽又免费观看的视频| 成人特级av手机在线观看| 中文字幕人妻熟人妻熟丝袜美| 美女免费视频网站| 一级黄片播放器| 国产老妇女一区| 免费观看人在逋| 国产欧美日韩精品一区二区| 精品一区二区免费观看| 欧美3d第一页| 一进一出抽搐gif免费好疼| 嫩草影院入口| 国产伦在线观看视频一区| 久久人人爽人人爽人人片va| 日本三级黄在线观看| 免费看a级黄色片| 成人性生交大片免费视频hd| 欧美区成人在线视频| 天美传媒精品一区二区| 真人做人爱边吃奶动态| 久久精品国产清高在天天线| 国产v大片淫在线免费观看| 国产乱人视频| 日本一二三区视频观看| 国内精品一区二区在线观看| 亚洲最大成人av| 精品久久久噜噜| 草草在线视频免费看| 欧美中文日本在线观看视频| 91久久精品电影网| 美女cb高潮喷水在线观看| 欧美不卡视频在线免费观看| 黄色视频,在线免费观看| 精品国内亚洲2022精品成人| 丰满的人妻完整版| 欧美成人免费av一区二区三区| 此物有八面人人有两片| 久久午夜亚洲精品久久| 亚洲av二区三区四区| 男女之事视频高清在线观看| 中文字幕av在线有码专区| 美女高潮喷水抽搐中文字幕| 日韩欧美在线二视频| 久久国产乱子免费精品| 国产麻豆成人av免费视频| 精品一区二区三区视频在线观看免费| 丝袜美腿在线中文| 伦精品一区二区三区| 91久久精品国产一区二区三区| 黄色欧美视频在线观看| 国产蜜桃级精品一区二区三区| 中文字幕免费在线视频6| 久久久午夜欧美精品| 亚洲国产色片| 亚洲第一区二区三区不卡| 午夜福利在线观看吧| 黄色一级大片看看| a级毛片免费高清观看在线播放| 国内精品久久久久精免费| 国产91精品成人一区二区三区| 色av中文字幕| 美女大奶头视频| 日本黄色片子视频| av天堂在线播放| 欧美xxxx性猛交bbbb| 日本黄大片高清| 给我免费播放毛片高清在线观看| av福利片在线观看| 欧美在线一区亚洲| 亚洲最大成人中文| 在线观看美女被高潮喷水网站| 小蜜桃在线观看免费完整版高清| 国产亚洲91精品色在线| 午夜老司机福利剧场| 伊人久久精品亚洲午夜| 欧美性感艳星| 在线播放无遮挡| 露出奶头的视频| 看片在线看免费视频| 国产国拍精品亚洲av在线观看| 男女那种视频在线观看| 国产私拍福利视频在线观看| 黄色丝袜av网址大全| 草草在线视频免费看| 亚洲男人的天堂狠狠| 精品久久久久久久久亚洲 | 最好的美女福利视频网| 18禁黄网站禁片午夜丰满| 身体一侧抽搐| 成年人黄色毛片网站| 夜夜看夜夜爽夜夜摸| 日本一二三区视频观看| 夜夜看夜夜爽夜夜摸| 日韩精品有码人妻一区| 国产精品久久电影中文字幕| 欧美bdsm另类| 亚洲av不卡在线观看| 国产一区二区在线观看日韩| 99热这里只有是精品在线观看| 麻豆久久精品国产亚洲av| 久久九九热精品免费| 精品不卡国产一区二区三区| 亚洲黑人精品在线| 国产精品一区二区三区四区免费观看 | 亚洲欧美日韩无卡精品| 国产成人av教育| 91午夜精品亚洲一区二区三区 | 此物有八面人人有两片| 午夜激情福利司机影院| 22中文网久久字幕| 简卡轻食公司| 国产高清视频在线播放一区| 国产欧美日韩一区二区精品| 999久久久精品免费观看国产| 男人和女人高潮做爰伦理| 不卡一级毛片| 1024手机看黄色片| 12—13女人毛片做爰片一| 好男人在线观看高清免费视频| 两个人的视频大全免费| 国产淫片久久久久久久久| 一区二区三区高清视频在线| 精品欧美国产一区二区三| 1000部很黄的大片| 午夜爱爱视频在线播放| 又粗又爽又猛毛片免费看| 中文字幕精品亚洲无线码一区| 亚洲av中文av极速乱 | 国产精品无大码| 久久精品国产99精品国产亚洲性色| 国产 一区精品| 婷婷丁香在线五月| 午夜a级毛片| 尤物成人国产欧美一区二区三区| 国产私拍福利视频在线观看| 国产又黄又爽又无遮挡在线| 国产高清视频在线观看网站| 成人一区二区视频在线观看| 亚洲国产欧美人成| a级一级毛片免费在线观看| a级毛片a级免费在线| 欧美日韩国产亚洲二区| 搡老妇女老女人老熟妇| 日韩精品中文字幕看吧| 久久天躁狠狠躁夜夜2o2o| 中文亚洲av片在线观看爽| 国内精品宾馆在线| 91久久精品电影网| 非洲黑人性xxxx精品又粗又长| 亚洲精品一区av在线观看| 久久精品国产99精品国产亚洲性色| 搡老岳熟女国产| 九九爱精品视频在线观看| 色噜噜av男人的天堂激情| 亚洲人成网站在线播| 久久精品国产清高在天天线| 国产成人a区在线观看| 国产一区二区在线av高清观看| 免费av观看视频| 亚洲av日韩精品久久久久久密| 在线播放无遮挡| 亚洲av第一区精品v没综合| 亚洲av五月六月丁香网| 亚洲五月天丁香| 成年女人看的毛片在线观看| 精品一区二区三区视频在线| 亚洲精品在线观看二区| 国产亚洲欧美98| 久久人人精品亚洲av| 色视频www国产| 3wmmmm亚洲av在线观看| 尤物成人国产欧美一区二区三区| 精品午夜福利在线看| a级一级毛片免费在线观看| 极品教师在线免费播放| 成人高潮视频无遮挡免费网站| 最近最新免费中文字幕在线| 免费观看在线日韩| 亚洲人与动物交配视频| 成人三级黄色视频| 黄色丝袜av网址大全| 免费高清视频大片| 欧美日本视频| 成人午夜高清在线视频| 国产亚洲精品久久久久久毛片| 亚洲av免费高清在线观看| 此物有八面人人有两片| 久久精品国产99精品国产亚洲性色| 成人国产麻豆网| 99久久成人亚洲精品观看| АⅤ资源中文在线天堂| 色5月婷婷丁香| 亚洲欧美日韩卡通动漫| 在线天堂最新版资源| 99精品久久久久人妻精品| 亚洲乱码一区二区免费版| 亚洲久久久久久中文字幕| 97超视频在线观看视频| 精品久久久久久久人妻蜜臀av| 一个人免费在线观看电影| 国产白丝娇喘喷水9色精品| 给我免费播放毛片高清在线观看| 亚洲 国产 在线| 婷婷丁香在线五月| 色吧在线观看| 亚洲av中文字字幕乱码综合| 校园人妻丝袜中文字幕| 精品不卡国产一区二区三区| 可以在线观看的亚洲视频| 成人三级黄色视频| 天堂影院成人在线观看| 欧美性猛交黑人性爽| 少妇熟女aⅴ在线视频| 亚洲五月天丁香| 亚洲成人免费电影在线观看| 51国产日韩欧美| 婷婷色综合大香蕉| 国产高清三级在线| 狂野欧美激情性xxxx在线观看|