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

    采用粒子群優(yōu)化的基因轉(zhuǎn)錄差異分析模型

    2010-09-11 01:46:10馮偉興王科俊
    關(guān)鍵詞:高通個數(shù)測序

    馮偉興王科俊

    (哈爾濱工程大學(xué)模式識別與智能系統(tǒng)研究所,哈爾濱 150001)

    采用粒子群優(yōu)化的基因轉(zhuǎn)錄差異分析模型

    馮偉興*王科俊

    (哈爾濱工程大學(xué)模式識別與智能系統(tǒng)研究所,哈爾濱 150001)

    高通測序技術(shù)可以更加全面地檢測實(shí)驗樣本中基因轉(zhuǎn)錄水平,這使得對樣本間基因轉(zhuǎn)錄差異的分析精度越來越準(zhǔn)確。根據(jù)通過DNA高通測序技術(shù)獲得的基因轉(zhuǎn)錄區(qū)內(nèi)PolⅡ蛋白個數(shù),提出了兩個樣本間基因轉(zhuǎn)錄差異分析模型。該模型在考慮同一基因間轉(zhuǎn)錄差異的同時,引入了全基因的轉(zhuǎn)錄分布特性以提高模型分析精度。模型采用預(yù)測概率最大化統(tǒng)計算法進(jìn)行參數(shù)求取。在概率最大化步驟中,由于難以得到模型參數(shù)解析解和提高算法效率,采用粒子群優(yōu)化算法直接進(jìn)行求取模型參數(shù)數(shù)值解。乳腺癌實(shí)驗數(shù)據(jù)測試表明,該模型可有效分析不同樣本間基因轉(zhuǎn)錄差異。

    基因;轉(zhuǎn)錄差異;粒子群;分析模型

    Abstract:Gene transcriptional level in an experiment can be more widely measured with high-throughput sequencing technology,which can obviously improve the analysis accuracy of gene transcriptional change between two experimental samples.Using numbers of PolⅡproteins contained inside each gene transcription region checked with DNA sequencing technology,a novel model was proposed to analyze gene transcription change considering gene transcription levels and the whole genome transcriptional distribution factor as well.Expectation-maximization(EM)statistics algorithm was utilized to resolve the model parameters.For obtaining analytical solution of model parameters and improving the efficiency of algorithm,particle swarm optimization(PSO)was used in maximization step to directly search parameters optimization numerical values.The test to analyze gene transcription change between normal and breast cancer samples verified the effectiveness of the proposed model.

    Key words:gene;transcription difference;particle swarm;analysis model

    引言

    國際基因組測序計劃的順利完成,形成了由幾億個堿基代碼組成的DNA序列集,并在其上發(fā)現(xiàn)了25 000個以上的人體基因。在此基礎(chǔ)上,隨著研究的逐步深入,發(fā)現(xiàn)在眾多調(diào)控因子高度非線性的綜合調(diào)控下,以各種功能蛋白為橋梁,基因間形成著錯綜復(fù)雜的調(diào)控關(guān)系。一個基因的表達(dá)發(fā)生變異,就有可能導(dǎo)致整個細(xì)胞功能異常,甚至癌病變。因此,在正常和異常兩個樣本間高精度地分析出表達(dá)發(fā)生變異的基因,可以為相關(guān)疾病診治提供有意義的參考依據(jù)。

    隨著計算機(jī)、控制和信息技術(shù)的飛速發(fā)展,人們獲取細(xì)胞樣本內(nèi)基因表達(dá)數(shù)據(jù)的方法越來越成熟[1]。20世紀(jì)末出現(xiàn)并逐步發(fā)展起來的基因芯片技術(shù),提供了測量樣本中基因表達(dá)程度的有效方法,基于通過基因芯片技術(shù)獲得的基因表達(dá)數(shù)據(jù),Kerr和Wolfinger設(shè)計了方差分析模型進(jìn)行不同樣本間基因表達(dá)差異測試[2-3]。Dudoit則采用統(tǒng)計學(xué)中的T測試法組合比較多樣本間的基因表達(dá)差異[4]。由于基因芯片僅能做到大體上的全基因表達(dá)測量,因此采用簡單的統(tǒng)計學(xué)方法進(jìn)行樣本間基因表達(dá)差異分析是足夠的。但隨著近幾年高通測序技術(shù)的出現(xiàn)和成熟,真正意義的全基因表達(dá)測量成為可能。與已有的不同樣本間基因表達(dá)差異分析方法相比,基于高通測序技術(shù),不僅可以參考樣本間基因表達(dá)測得數(shù)據(jù),還可借鑒樣本中基因表達(dá)的分布特性,從而可以更準(zhǔn)確地對基因表達(dá)差異進(jìn)行分析。

    相比較而言,DNA高通測序技術(shù)已經(jīng)比較成熟,可用來獲取DNA上所結(jié)合轉(zhuǎn)錄因子、PolⅡ等蛋白信息,該信息因可被用于分析基因前轉(zhuǎn)錄調(diào)控機(jī)制,并間接反映基因表達(dá)水平及其基因間調(diào)控關(guān)系而得到廣泛應(yīng)用[4-7]。在此基礎(chǔ)上,RNA高通測序技術(shù)也已逐漸成熟。該技術(shù)可直接反映基因表達(dá)水平并可用于基因變異分析。為了更全面地獲取有用信息,相比于一個基因芯片僅能獲得幾十萬條有用信息,目前一次高通測序周期可同時處理30億個堿基并生成接近1億條讀數(shù)。據(jù)估計,2009年年底該數(shù)字將達(dá)到1 000億個堿基和30億條讀數(shù)。面對更加全面的數(shù)據(jù),通過設(shè)計新方法對其進(jìn)行有效的生物信息處理,無疑將得到與原有技術(shù)相比精度更高的分析成果[8-11]。由于DNA高通測序技術(shù)發(fā)展得比較成熟,測得數(shù)據(jù)精度也較高,本文利用該技術(shù)測得的數(shù)據(jù)進(jìn)行不同樣本間基因表達(dá)差異分析。具體就是利用DNA高通測序技術(shù)測量基因轉(zhuǎn)錄區(qū)內(nèi)PolⅡ蛋白數(shù)目,進(jìn)而分析不同樣本間基因轉(zhuǎn)錄差異。盡管基因轉(zhuǎn)錄差異是基因表達(dá)差異的有效間接反映,但對于分析不同樣本間致病基因及其致病機(jī)理仍然具有很大的參考價值?;诖?,本文利用DNA高通測序技術(shù)測得的數(shù)據(jù)提出一種分析不同樣本間基因轉(zhuǎn)錄差異的分析模型,該模型不僅參考樣本間基因轉(zhuǎn)錄測得數(shù)據(jù),還借鑒了樣本中基因轉(zhuǎn)錄的分布特性。鑒于模型的復(fù)雜性,在模型參數(shù)的求取過程中采用智能尋優(yōu)算法。

    1 基因轉(zhuǎn)錄差異分析模型

    在基因轉(zhuǎn)錄過程中,PolⅡ蛋白結(jié)合在基因轉(zhuǎn)錄區(qū)內(nèi)并對基因進(jìn)行轉(zhuǎn)錄,因此,基于基因轉(zhuǎn)錄區(qū)內(nèi)PolⅡ蛋白數(shù)量可以進(jìn)行不同樣本間基因轉(zhuǎn)錄差異分析[12]。

    1.1 基因轉(zhuǎn)錄數(shù)據(jù)獲取

    利用DNA高通測序中的 ChIP-seq技術(shù)可以對基因轉(zhuǎn)錄區(qū)內(nèi)PolⅡ蛋白數(shù)量進(jìn)行測量和統(tǒng)計,以直接反映基因轉(zhuǎn)錄水平。該技術(shù)首先利用超聲波將DNA鏈降解為DNA片段,然后利用特制的抗體俘獲結(jié)合在DNA片段上的PolⅡ蛋白,再利用沉淀技術(shù)(IP)將含有抗體的DNA片段濾出,隨后通過測序技術(shù)(seq)對所有濾出的DNA片段測序并通過序列比對映射回DNA上,最后根據(jù)基因轉(zhuǎn)錄區(qū)在DNA上的位置定義即可實(shí)現(xiàn)對基因轉(zhuǎn)錄區(qū)內(nèi)PolⅡ蛋白數(shù)量的測量和統(tǒng)計[13]。

    1.2 基因轉(zhuǎn)錄差異分析

    在兩個樣本之間,例如正常樣本和癌癥樣本,為了分析各個基因的轉(zhuǎn)錄是否發(fā)生變化以及變化程度,這里分為基因轉(zhuǎn)錄變化和基因轉(zhuǎn)錄不變化兩類。然后,依據(jù)不同樣本中分布在各基因轉(zhuǎn)錄區(qū)內(nèi)PolⅡ蛋白的個數(shù)及其個數(shù)變化程度來判斷各個基因?qū)儆谀囊活愐约半`屬程度。為此,提出基因轉(zhuǎn)錄模型及其轉(zhuǎn)錄差異模型并對此進(jìn)行分析。

    1.2.1 基因轉(zhuǎn)錄模型

    由于細(xì)胞中對全部n個基因的調(diào)控強(qiáng)度不同,因此一個PolⅡ蛋白結(jié)合在不同基因轉(zhuǎn)錄區(qū)內(nèi)以對其進(jìn)行轉(zhuǎn)錄的概率是不同的,設(shè)該概率為Pi(i=1,…,n)。假設(shè)細(xì)胞中PolⅡ蛋白的個數(shù)處于較低的水平,設(shè)其個數(shù)為m,則在第i個基因的轉(zhuǎn)錄區(qū)內(nèi)結(jié)合k個PolⅡ蛋白的概率可用二項分布表示:

    按照這一概率分布,隨著細(xì)胞中PolⅡ蛋白的個數(shù)增多,結(jié)合在各基因轉(zhuǎn)錄區(qū)內(nèi)的PolⅡ蛋白個數(shù)將隨之增加。但隨著細(xì)胞中PolⅡ蛋白的個數(shù)繼續(xù)增多,由于基因轉(zhuǎn)錄區(qū)上可容納的PolⅡ蛋白的個數(shù)有限,新的PolⅡ蛋白結(jié)合在基因轉(zhuǎn)錄區(qū)內(nèi)的概率將隨之下降,最終,基因轉(zhuǎn)錄區(qū)內(nèi)結(jié)合的PolⅡ蛋白個數(shù)將出現(xiàn)飽和,即mPi(i=1,…,n)將保持不變。事實(shí)上,為了保證正常生理功能,細(xì)胞中PolⅡ蛋白的個數(shù)多處于飽和狀態(tài)。此時,各基因轉(zhuǎn)錄區(qū)內(nèi)結(jié)合PolⅡ蛋白的個數(shù)應(yīng)遵循泊松分布:

    式中,yi,j是第 j(j=1,2)個樣本中第 i個基因轉(zhuǎn)錄區(qū)內(nèi)結(jié)合 Pol Ⅱ蛋白的個數(shù),而λi,j=mjPi,j為分布常數(shù)。按照泊松分布特性,實(shí)際測量的PolⅡ蛋白個數(shù) yi,j是基因 i在樣本 j中實(shí)際轉(zhuǎn)錄程度 λi,j的表象,但如何通過 yi,j準(zhǔn)確 計算出 λi,j,僅依靠yi,j是無法做到的。此外,即使計算出 λi,j,如何評估其樣本間變化程度也是很困難的。一個基因轉(zhuǎn)錄從100變?yōu)?00和另一個基因轉(zhuǎn)錄從1 000變?yōu)? 000,其表示的轉(zhuǎn)錄差異程度是不同的,前者更劇烈一些。為此,結(jié)合樣本中所有基因的表達(dá)分布特性,設(shè)計了基因轉(zhuǎn)錄差異模型分析樣本間基因的轉(zhuǎn)錄差異及其差異程度。

    1.2.2 基因轉(zhuǎn)錄差異模型

    據(jù)Newton實(shí)驗分析,在同一樣本中,所有基因的表達(dá)水平將遵循伽瑪分布[14~15]。基于此,可認(rèn)為:同一樣本內(nèi),基因轉(zhuǎn)錄實(shí)際值 λi(i=1,…,n)應(yīng)遵循伽瑪分布:

    其中,α,β為模型常數(shù)。

    兩個樣本之間,如果第i個基因的轉(zhuǎn)錄未發(fā)生變化,則在兩個樣本PolⅡ總數(shù)歸一化前提下,其在兩個樣本內(nèi)均遵循同一伽瑪分布,即:

    否則:

    為了計算基因轉(zhuǎn)錄差異程度,本研究定義變量Zi來描述兩個樣本之間第i個基因的轉(zhuǎn)錄是否發(fā)生變化,如果完全變化則Zi=1;完全不變化則Zi=0。

    基于以上分析,本研究最終構(gòu)建了一個概率模型PA來綜合描述細(xì)胞中所有基因出現(xiàn)當(dāng)前轉(zhuǎn)錄變化情形的概率:

    其中,yi1,yi2為基因 i在兩個樣本中測得的轉(zhuǎn)錄值。λi1,λi2為基因 i在兩個樣本中的實(shí)際轉(zhuǎn)錄值。α1,β1,α2,β2,Zi為模型參數(shù)。P 為任一個基因在兩個樣本之間發(fā)生轉(zhuǎn)錄變化概率。

    如果概率模型PA的模型參數(shù)得到準(zhǔn)確求取,則通過Zi值就可以分析出第i個基因在不同樣本間,其轉(zhuǎn)錄是否發(fā)生變化以及變化程度。Zi越大則變化程度越大,反之亦然。

    在上述概率模型中,如果Zi已知,即第i個基因在兩個樣本間轉(zhuǎn)錄變化程度是已知的,則可以通過調(diào)整 α1,β1,α2,β2的取值,直接計算出 λi,λi1,λi2值并使得PA最大化,從而實(shí)現(xiàn)對兩個樣本中基因轉(zhuǎn)錄情形的透徹分析。但由于Zi的取值是未知的,無法一次對所有模型參數(shù)在PA最大化前提下進(jìn)行求取,因此,采用統(tǒng)計學(xué)中預(yù)測概率最大化迭代算法來對 α1,β1,α2,β2進(jìn)行 PA概率最大化條件下的求取,并最終確定 λi,λi1,λi2及其 Zi值,以確定每個基因是否在不同樣本間發(fā)生轉(zhuǎn)錄變化及其變化程度。

    這里采用的預(yù)測概率最大化算法分為反復(fù)迭代的兩步。

    第一步是概率預(yù)測,即假設(shè) α1,β1,α2,β2已知的情況下求取Zi值。

    首先對 λi,λi1,λi2進(jìn)行估算,即在 yi1,yi2,α1,β1,α2,β2已知的情況下尋找 λi,λi1,λi2值使 PA最大。

    若Zi=0:

    若Zi=1:

    隨后,Zi可估算為:第二步是概率最大化,即在計算出的情況下尋找最優(yōu)的 α1,β1,α2,β2使 PA最大。通過對 PA表達(dá)式(6)的分析發(fā)現(xiàn),這相當(dāng)于通過尋找 α1,β1,α2,β2值使下兩式最大

    最后,對基因在兩個樣本之間發(fā)生轉(zhuǎn)錄變化的概率P值進(jìn)行更新:

    以上過程迭代進(jìn)行,直至收斂。

    1.3 粒子群優(yōu)化

    在應(yīng)用預(yù)測概率最大化算法進(jìn)行概率模型PA的模型參數(shù)求取過程中發(fā)現(xiàn),針對第一步,即概率預(yù)測,其計算比較容易實(shí)現(xiàn)。但對于第二步,即概率最大化,如果進(jìn)行精確的解析求解將很困難且計算量相當(dāng)大。與之相對應(yīng),在復(fù)雜解空間內(nèi)搜尋最優(yōu)解恰好是智能尋優(yōu)算法的優(yōu)勢[16],因此,這里采用智能尋優(yōu)算法中的粒子群算法,在概率最大化步驟中對模型參數(shù)直接進(jìn)行數(shù)值解尋優(yōu)求取。通過初步計算,發(fā)現(xiàn)模型參數(shù) α1,β1和 α2,β2的解空間比較簡單,因此采用基本粒子群算法進(jìn)行 α1,β1和α2,β2的數(shù)值解求取。另外,在式(10)中,由于 α1,β1和 α2,β2不具有相關(guān)性,因此這里可對 α1,β1和α2,β2分別進(jìn)行尋優(yōu)。

    具體粒子群算法設(shè)計如下:

    算法初始化時,取h個粒子,其初始值均勻分布在二維解空間內(nèi)。對于每個粒子,其解空間內(nèi)當(dāng)前位置矢量的評估函數(shù)采用式(10)。對于第k(k=1,…,h)個粒子,為了在解空間內(nèi)搜索更優(yōu)的位置矢量,其移動速度V′k采用下式進(jìn)行更新

    式中,Pk為第k個粒子在解空間內(nèi)當(dāng)前位置矢量。Pglobal為所有粒子中解空間內(nèi)當(dāng)前最優(yōu)的位置矢量。Pk-local為第 k個粒子解空間內(nèi)當(dāng)前最優(yōu)的位置矢量。Vk為第 k個粒子當(dāng)前移動速度。C0,C1,C2為各部分權(quán)值以控制第k個粒子在解空間內(nèi)的移動方向和速度。

    最后,第k個粒子在解空間內(nèi)的新位置矢量按下式計算:

    以上過程重復(fù)進(jìn)行,直至無法在解空間內(nèi)找到更優(yōu)的全局位置矢量。

    2 實(shí)驗測試結(jié)果和討論

    2.1 實(shí)驗設(shè)計

    為了對本研究所提出基因轉(zhuǎn)錄差異分析模型的有效性進(jìn)行測試,這里將選取美國俄亥俄州立大學(xué)公開提供的兩種MCF-7乳腺癌樣本PolⅡ測量數(shù)據(jù)。一種樣本是普通的MCF-7乳腺癌樣本,另一種是具有抗藥性的 MCF-7乳腺癌樣本。該數(shù)據(jù)共分為四組,分別是加藥前普通乳腺癌、加藥后普通乳腺癌、加藥前抗藥乳腺癌和加藥后抗藥乳腺癌樣本中利用DNA高通測序技術(shù)測得的25470個基因轉(zhuǎn)錄區(qū)內(nèi)結(jié)合的PolⅡ蛋白數(shù)量。

    按預(yù)計,在加藥的情況下,普通乳腺癌樣本內(nèi)基因的轉(zhuǎn)錄變化將比抗藥乳腺癌樣本更劇烈,這構(gòu)造出測試樣本分析模型有效性的良好環(huán)境。

    對這四組實(shí)驗數(shù)據(jù)在進(jìn)行PolⅡ總數(shù)歸一化后的統(tǒng)計分析結(jié)果如表1所示。

    表1 測試數(shù)據(jù)統(tǒng)計分析(單位:個)Tab.1 Statistical analysis of experiment data(Unit:Number )

    由表1可以分析出,在4組測試數(shù)據(jù)歸一化后,不同樣本間的測試數(shù)據(jù)在均值一致的前提下方差變化有著明顯的差異。其中,普通乳腺癌樣本在加藥后基因整體的轉(zhuǎn)錄差異明顯降低,這可能是由于部分轉(zhuǎn)錄水平高的基因加藥后轉(zhuǎn)錄水平得到抑制導(dǎo)致的;而抗藥乳腺癌樣本加藥后基因的整體轉(zhuǎn)錄水平變化則不明顯,這說明其對藥物不敏感。這一現(xiàn)象是符合兩種樣本的自然特性的,因此,可證明實(shí)驗測得數(shù)據(jù)是有意義的。

    在獲取所有基因整體的轉(zhuǎn)錄差異變化后,具體到每一個基因,如何評估其轉(zhuǎn)錄變化則將用本研究所提出分析模型進(jìn)行分析。

    2.2 基因轉(zhuǎn)錄差異分析

    首先,如果某個基因所在DNA區(qū)域內(nèi)組蛋白修飾或DNA甲基化等原因?qū)е翽olⅡ蛋白無法在其轉(zhuǎn)錄區(qū)內(nèi)結(jié)合,則該基因轉(zhuǎn)錄將不遵循泊松分布。因此,為保證計算精度,在測試數(shù)據(jù)中將去除加藥前后轉(zhuǎn)錄區(qū)內(nèi)均沒有結(jié)合PolⅡ蛋白的基因。這樣處理后,在全部25 470個基因中,普通乳腺癌樣本測試數(shù)據(jù)中將有24 020個基因參與后續(xù)分析,而抗藥乳腺癌樣本測試數(shù)據(jù)中有24 230個基因參與后續(xù)分析。

    隨后,將4組轉(zhuǎn)錄數(shù)據(jù)分成兩個大組。每個大組含有一個乳腺癌樣本加藥前后的兩組基因轉(zhuǎn)錄數(shù)據(jù)。然后分別采用本文所提出的基因轉(zhuǎn)錄差異分析模型對其進(jìn)行個體基因轉(zhuǎn)錄差異分析。

    在預(yù)測概率最大化算法的迭代參數(shù)優(yōu)化過程中,在概率預(yù)測步驟中,由于針對n個個體基因構(gòu)造的泊松模型參數(shù) λi,λi1,λi2(i=1,…,n)的最優(yōu)參數(shù)計算式(7)和式(8)也沒有解析解,依據(jù)其所遵循的伽瑪函數(shù)凸特性設(shè)計參數(shù)尋優(yōu)算法如下:λi,λi1,λi2分別由1開始以1為步長進(jìn)行增加并計算其對應(yīng)的式(7)和式(8)中伽瑪函數(shù)值,當(dāng)函數(shù)值由小變大過程中首次出現(xiàn)下降時,在其對應(yīng)的 λi,λi1,λi2值附近進(jìn)行尋優(yōu),并找到 λi,λi1,λi2的最優(yōu)值。在概率最大化步驟中,將采用粒子群優(yōu)化算法。粒子群優(yōu)化算法的參數(shù)設(shè)置為:粒子個數(shù)為100,C0=1,C1=2,C2=2。

    對于普通乳腺癌和抗藥乳腺癌兩大組基因轉(zhuǎn)錄數(shù)據(jù),預(yù)測概率最大化算法的參數(shù)尋優(yōu)均在迭代4~5次后收斂,收斂后的 Zi(i=1,2,…,n)值將反映第i個基因轉(zhuǎn)錄變化的程度。其最終分析結(jié)果見表2和表3。

    表2 普通乳腺癌分析結(jié)果Tab.2 Analytical results of normal breast cancer

    表3 抗藥乳腺癌分析結(jié)果Tab.3 Analytical results of drug resist brea st cancer

    結(jié)合表1和表2,普通乳腺癌樣本加藥后,盡管基因間整體轉(zhuǎn)錄差異有所降低,但就個體基因而言,在所考慮的24 020個基因中,仍有6505個基因的轉(zhuǎn)錄發(fā)生了明顯變化,占27.1%。其中,轉(zhuǎn)錄提高的基因共有4 346(18.1%)個,提高幅度大的(Z≥0.9)有2 474(10.3%)個,提高幅度小的(0.1≤Z<0.9)有1 872(7.8%)個。而轉(zhuǎn)錄降低的基因僅為轉(zhuǎn)錄提高基因數(shù)目的一半,有2 159(9.0%)個,其中降低大的為1 373(5.7%),降低小的為786(3.3%)??梢娂铀幒?,普通乳腺癌樣本的反應(yīng)還是較劇烈的,更多的反應(yīng)是基因的轉(zhuǎn)錄水平發(fā)生了提高。

    結(jié)合表1和表3,抗藥乳腺癌樣本加藥后,盡管基因間整體轉(zhuǎn)錄差異變化不大,但就個體基因而言,仍有6 505(8.4%)個基因的轉(zhuǎn)錄發(fā)生了明顯變化。其中,轉(zhuǎn)錄提高的基因較少,僅有489(2.0%)個,提高幅度大的有123(0.5%)個,提高幅度小的有366(1.5%)個。而轉(zhuǎn)錄降低的基因有1 543(6.4%)個,其中降低大的為1 273(5.3%),降低小的為270(1.1%)??梢娂铀幒?,抗藥乳腺癌樣本的反應(yīng)相比普通乳腺癌樣本明顯減弱。而這種減弱主要體現(xiàn)在加藥后基因的轉(zhuǎn)錄反而受到了抑制,這從基因轉(zhuǎn)錄降低的數(shù)目明顯多于基因轉(zhuǎn)錄提高數(shù)目中可以看出。

    從以上分析可以看出,加藥對普通乳腺癌樣本和抗藥乳腺癌樣本的影響有著明顯的不同,相比較而言,普通乳腺癌樣本比抗藥乳腺癌樣本對加藥的反應(yīng)要強(qiáng)烈得多。而這種差異與實(shí)驗設(shè)計的預(yù)期結(jié)果是相吻合的。這充分反映出本研究所提出基因轉(zhuǎn)錄差異分析模型的有效性。

    3 結(jié)論

    利用可直接反映基因轉(zhuǎn)錄水平的基因轉(zhuǎn)錄區(qū)內(nèi)所結(jié)合PolⅡ 蛋白的個數(shù),提出一個分析模型,分析每一個基因在兩個相關(guān)樣本間轉(zhuǎn)錄水平的變化程度,通過實(shí)際數(shù)據(jù)測試證明了該模型的有效性。

    由于不同樣本不同實(shí)驗下所測得PolⅡ 蛋白整體個數(shù)差異較大,因此,即使通過歸一化進(jìn)行了校正,仍然會對本研究所提出模型的分析精度造成影響。目前,隨著DNA測序技術(shù)的發(fā)展,一次實(shí)驗所測得的數(shù)據(jù)越來越大。據(jù)估計,2009年年底該數(shù)字將達(dá)到1000億個堿基和30億條讀數(shù)。隨著測試數(shù)據(jù)量的提高,不同實(shí)驗下所測得PolⅡ 蛋白整體個數(shù)差異對本模型分析精度的影響將會越來越小。

    另外,為了使分析結(jié)果有明顯的生物含義,所提出的模型是基于這樣一個前提設(shè)計的:即如果兩個樣本間基因轉(zhuǎn)錄水平未發(fā)生明顯變化,則其轉(zhuǎn)錄分布特性是基本不變的。這樣處理的好處是如果找到兩個樣本間轉(zhuǎn)錄分布特性發(fā)生明顯變化的基因,該基因就很可能與兩個樣本的自然差異有關(guān)(如加藥反應(yīng),癌病變等)。因此,只要遵從這一前提,本研究所提出模型即可用于對兩個樣本間基因轉(zhuǎn)錄水平的差異進(jìn)行分析,也可以用于多樣本分析。具體方法可以以一個樣本為基準(zhǔn),其它樣本對其分別利用本研究所提出模型進(jìn)行計算,并將所有計算結(jié)果集中在一起進(jìn)行綜合分析即可。

    此外,由于PolⅡ 蛋白的個數(shù)只能反映基因轉(zhuǎn)錄水平,不能直接反映更有生物意義的基因表達(dá)水平,因此本模型還有待進(jìn)一步發(fā)展,以最終用于基因表達(dá)差異分析。目前,最新發(fā)展的生物技術(shù)已使RNA高通測序成為可能。相比于DNA高通測序僅在DNA層次反映樣本組織的基因轉(zhuǎn)錄模式,RNA高通測序則可直接測量樣本組織中所有基因是否表達(dá)以及表達(dá)的程度。這就為直接對基因表達(dá)差異分析提供了可能。針對這一新技術(shù),通過將分析對象由基因轉(zhuǎn)錄區(qū)內(nèi)所結(jié)合 PolⅡ 蛋白個數(shù)替換為反映該基因表達(dá)水平的RNA高通測序序列讀數(shù),本模型即可用于分析兩個相關(guān)樣本間基因表達(dá)差異分析,但具體實(shí)現(xiàn)方式還有待針對RNA高通測序技術(shù)進(jìn)行進(jìn)一步研究。

    [1]張曉龍,楊艷霞.機(jī)器學(xué)習(xí)在生物信息學(xué)中的應(yīng)用[J].武漢科技大學(xué)學(xué)報,2005,28(2):201-204

    [2]Kerr MK,Martin M,Churchill GA.Analysis of variance for gene expression microarray data[J].Journal of Computational Biology,2000,7:819-837.

    [3]Wolfinger RD,Gibson G,Wolfinger ED.et al.Assessing Gene Significance from cDNA Microarray Expression Data via Mixed Models[J].Journal of Computational Biology,2001,8:625-637.

    [4]Dudoit S,Yang YH,Callow MJ,et al.Statistical methods for identifying differentially expressed genes in replicated cDNA microarray experiments[J].Statistica Sinica,2002,12:111-139.

    [5]Schmidt D,Wilson MD,Spyrou C,et al.ChIP-seq:Using high-throughput sequencing to discover protein-DNA interactions[J].Methods,2009,48(3):240-247.

    [6]Johnson DS,Mortazavi A,Myers RM,et al.Genome-wide mapping of in vivo protein-DNA interactions[J].Science,2007,316:1497-1499.

    [7]Ansorge WJ.Next-generation DNA sequencing techniques[J].N Biotechnol,2009,25(4):195-203.

    [8]Fullwood MJ,Wei CL,Liu ETB,et al.Next-generation DNA sequencing of paired-end tags(PET)for transcriptome and genome analyses[J].Genome Res,2009,19(4):521-32.

    [9]Cole T,Lior P,Steven LS.TopHat:discovering splice junctions with RNA-Seq[J].Bioinformatics,2009,25(9):1105-1111.

    [10]Daniel RZ,Ewan B.Velvet:Algorithms for de novo short read assembly using de Bruijn graphs[J].Genome Res,2008,18(5):821-829.

    [11]David JS,William GR.Transcriptome sequencing of malignant pleural mesothelioma tumors[J].PNAS,2008,105(9):3521-3526.

    [12]特納PC著.分子生物學(xué)[M].(第2版).北京:科學(xué)出版社,2001.

    [13]Raja J,Suresh C.Genome-wide identification of in vivo protein-DNA binding sites from ChIP-Seq data[J].Nucleic Acids Research,2008,36(16):5221-5231.

    [14]Newton MA,KendziorskiCM,RichmondCS,etal.On differential variability of expression ratios:Improving statistical inference about gene expression changes from microarray data[J].Journal of Computational Biology,2001,8:37-52.

    [15]Kendziorski C,Newton M,Lan H,et al.On parametric empirical bayes methods for comparing multiple groups using replicated gene expression profiles[J].Stat Med,2003,22(24):3899-3914.

    [16]Parsopoulos KE,Vrahatis MN.Recent approaches to global optimization problems through Particle Swarm Optimization[J].Natural Computing,2002,1:235-306.

    A Model for Analyzing Gene Transcription Difference Using Particle Swarm Optimization

    FENG Wei-Xing*WANG Ke-Jun
    (Pattern Recognition and Intelligent System Institute of Harbin Engineering University,Harbin,Heilongjiang 150001,China)

    Q332

    A

    0258-8021(2010)02-0229-06

    10.3969/j.issn.0258-8021.2010.02.012

    2009-07-10,

    2009-11-04

    國家高技術(shù)研究發(fā)展(863)計劃(2008AA01Z148)

    *通訊作者。 E-mail:fengweixing@hrbeu.edu.cn

    猜你喜歡
    高通個數(shù)測序
    杰 Sir 帶你認(rèn)識宏基因二代測序(mNGS)
    新民周刊(2022年27期)2022-08-01 07:04:49
    怎樣數(shù)出小正方體的個數(shù)
    二代測序協(xié)助診斷AIDS合并馬爾尼菲籃狀菌腦膜炎1例
    傳染病信息(2021年6期)2021-02-12 01:52:58
    等腰三角形個數(shù)探索
    怎樣數(shù)出小木塊的個數(shù)
    怎樣數(shù)出小正方體的個數(shù)
    歷史轉(zhuǎn)折后的高通前執(zhí)行董事長
    外資力薦高通并購賽靈思對抗英特爾
    電子世界(2015年22期)2015-12-29 02:49:41
    高通24億美元收購芯片制造商CSR
    IT時代周刊(2015年9期)2015-11-11 05:51:53
    高通市場份額跌至47%
    IT時代周刊(2015年8期)2015-11-11 05:50:48
    亚洲欧美日韩卡通动漫| 日本猛色少妇xxxxx猛交久久| 九草在线视频观看| 色5月婷婷丁香| 欧美xxⅹ黑人| 亚洲一级一片aⅴ在线观看| 少妇猛男粗大的猛烈进出视频| 少妇人妻精品综合一区二区| 美女中出高潮动态图| 99国产综合亚洲精品| 99热国产这里只有精品6| 国产精品免费大片| av女优亚洲男人天堂| 十八禁高潮呻吟视频| 免费大片黄手机在线观看| 91aial.com中文字幕在线观看| 欧美国产精品一级二级三级| 最近2019中文字幕mv第一页| 欧美变态另类bdsm刘玥| 国产精品国产三级专区第一集| 亚洲精品色激情综合| 国模一区二区三区四区视频| 国产亚洲精品第一综合不卡 | 亚洲一区二区三区欧美精品| 亚洲欧洲精品一区二区精品久久久 | 三上悠亚av全集在线观看| 午夜免费观看性视频| 制服丝袜香蕉在线| 高清黄色对白视频在线免费看| 亚洲国产精品专区欧美| av卡一久久| 18禁动态无遮挡网站| 日韩一区二区视频免费看| 欧美xxⅹ黑人| 免费av中文字幕在线| 成年人午夜在线观看视频| 高清毛片免费看| 免费看光身美女| 亚洲av国产av综合av卡| 青青草视频在线视频观看| 久久精品国产亚洲网站| 伦精品一区二区三区| 蜜臀久久99精品久久宅男| 欧美另类一区| 欧美3d第一页| 91精品三级在线观看| 2022亚洲国产成人精品| 久久青草综合色| 国产一级毛片在线| videosex国产| 免费播放大片免费观看视频在线观看| 亚洲国产最新在线播放| 国产成人精品一,二区| 高清午夜精品一区二区三区| 成人国产av品久久久| 亚洲欧美色中文字幕在线| 亚洲第一av免费看| 寂寞人妻少妇视频99o| 爱豆传媒免费全集在线观看| 欧美 亚洲 国产 日韩一| 91成人精品电影| 99精国产麻豆久久婷婷| 一级毛片 在线播放| 成人无遮挡网站| 伊人久久国产一区二区| 交换朋友夫妻互换小说| 狂野欧美激情性xxxx在线观看| 人妻少妇偷人精品九色| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 成年女人在线观看亚洲视频| 汤姆久久久久久久影院中文字幕| 日本爱情动作片www.在线观看| 亚洲av综合色区一区| 亚洲久久久国产精品| 亚洲美女黄色视频免费看| 国产综合精华液| 欧美亚洲日本最大视频资源| 久久女婷五月综合色啪小说| 日韩中文字幕视频在线看片| 国产黄色视频一区二区在线观看| 欧美激情国产日韩精品一区| 一区二区av电影网| 国产精品久久久久成人av| 涩涩av久久男人的天堂| 一级毛片aaaaaa免费看小| 日日爽夜夜爽网站| 老司机亚洲免费影院| a级毛片黄视频| 王馨瑶露胸无遮挡在线观看| 中文乱码字字幕精品一区二区三区| a 毛片基地| 国产无遮挡羞羞视频在线观看| 欧美激情国产日韩精品一区| 视频区图区小说| 久久午夜福利片| av黄色大香蕉| 国产永久视频网站| av国产精品久久久久影院| 特大巨黑吊av在线直播| 你懂的网址亚洲精品在线观看| 国产欧美另类精品又又久久亚洲欧美| 97精品久久久久久久久久精品| 日日摸夜夜添夜夜添av毛片| 日本av免费视频播放| 一区二区三区四区激情视频| 午夜福利视频在线观看免费| 亚洲精品亚洲一区二区| 免费观看性生交大片5| 国产爽快片一区二区三区| 免费av不卡在线播放| 亚洲av男天堂| 熟女人妻精品中文字幕| 黄色配什么色好看| 亚洲av福利一区| 黄色视频在线播放观看不卡| 国产精品不卡视频一区二区| 国产精品一国产av| 性色avwww在线观看| 日韩视频在线欧美| 男人操女人黄网站| 99久久精品国产国产毛片| 久热这里只有精品99| 亚洲av国产av综合av卡| av国产久精品久网站免费入址| 午夜福利视频在线观看免费| 国产永久视频网站| 国产黄片视频在线免费观看| 国产精品麻豆人妻色哟哟久久| 免费观看在线日韩| 成人18禁高潮啪啪吃奶动态图 | 看十八女毛片水多多多| 欧美日韩视频精品一区| 在线观看免费日韩欧美大片 | 久久99精品国语久久久| 成人毛片60女人毛片免费| 男人操女人黄网站| 18禁在线播放成人免费| 久久国产亚洲av麻豆专区| 成年人午夜在线观看视频| 99久久人妻综合| a级片在线免费高清观看视频| 久久精品人人爽人人爽视色| 草草在线视频免费看| 色吧在线观看| xxxhd国产人妻xxx| 老司机亚洲免费影院| 亚洲少妇的诱惑av| 少妇的逼水好多| 日日摸夜夜添夜夜添av毛片| 黄片无遮挡物在线观看| 色5月婷婷丁香| 不卡视频在线观看欧美| 天天操日日干夜夜撸| 精品国产一区二区久久| a级毛片免费高清观看在线播放| 中文字幕精品免费在线观看视频 | 天堂俺去俺来也www色官网| 国产av一区二区精品久久| 国产精品久久久久久久电影| 亚洲图色成人| 尾随美女入室| 国产一区二区三区av在线| 成人影院久久| 精品少妇内射三级| 午夜激情福利司机影院| 丁香六月天网| 亚洲中文av在线| 成年美女黄网站色视频大全免费 | 欧美日韩综合久久久久久| 亚洲av欧美aⅴ国产| 国产 一区精品| 精品午夜福利在线看| 丰满少妇做爰视频| tube8黄色片| 亚洲精华国产精华液的使用体验| 日本av手机在线免费观看| 亚洲国产最新在线播放| 亚洲欧美一区二区三区国产| 欧美一级a爱片免费观看看| 99国产综合亚洲精品| 欧美国产精品一级二级三级| 新久久久久国产一级毛片| 久久久久久久久久久久大奶| 亚洲高清免费不卡视频| 另类亚洲欧美激情| 我的女老师完整版在线观看| 热re99久久精品国产66热6| 99国产综合亚洲精品| 久久亚洲国产成人精品v| av福利片在线| 一边摸一边做爽爽视频免费| 国产高清不卡午夜福利| 精品少妇久久久久久888优播| 丝袜喷水一区| 欧美日韩精品成人综合77777| 夜夜骑夜夜射夜夜干| 蜜桃在线观看..| 搡女人真爽免费视频火全软件| 热99国产精品久久久久久7| 国产免费又黄又爽又色| 一个人看视频在线观看www免费| 久久99热6这里只有精品| av在线app专区| 色网站视频免费| 亚洲精品久久久久久婷婷小说| 国产亚洲欧美精品永久| 久久久久国产精品人妻一区二区| 精品酒店卫生间| 又粗又硬又长又爽又黄的视频| 久热久热在线精品观看| 久久久亚洲精品成人影院| 久久韩国三级中文字幕| 两个人免费观看高清视频| 日韩伦理黄色片| 亚洲国产精品成人久久小说| 国产高清不卡午夜福利| 国产高清有码在线观看视频| 少妇被粗大的猛进出69影院 | 日产精品乱码卡一卡2卡三| 精品少妇内射三级| 久久精品国产亚洲av涩爱| 97超视频在线观看视频| 高清毛片免费看| 午夜日本视频在线| 国产精品熟女久久久久浪| 伊人久久精品亚洲午夜| 亚洲欧洲精品一区二区精品久久久 | www.av在线官网国产| 国产男人的电影天堂91| 久久精品国产a三级三级三级| 亚洲av国产av综合av卡| 亚洲国产色片| 免费观看av网站的网址| 午夜免费鲁丝| 久久毛片免费看一区二区三区| 亚洲国产精品专区欧美| 中文字幕制服av| 久久热精品热| 女性被躁到高潮视频| 黄色一级大片看看| 9色porny在线观看| 午夜免费男女啪啪视频观看| 久久人人爽人人爽人人片va| 免费黄频网站在线观看国产| 亚洲av成人精品一二三区| 一边摸一边做爽爽视频免费| 人成视频在线观看免费观看| 性色avwww在线观看| 99国产精品免费福利视频| 国产国拍精品亚洲av在线观看| 在线亚洲精品国产二区图片欧美 | 中文字幕制服av| 亚洲精品自拍成人| 日韩亚洲欧美综合| 国产熟女欧美一区二区| 99久久综合免费| 久久亚洲国产成人精品v| 亚洲精品乱久久久久久| 亚洲精品国产av蜜桃| 久久久午夜欧美精品| 一本大道久久a久久精品| 久久精品国产亚洲av涩爱| 精品久久国产蜜桃| 一级爰片在线观看| 极品人妻少妇av视频| 另类精品久久| 国产精品国产三级专区第一集| 久久国内精品自在自线图片| 免费av不卡在线播放| 激情五月婷婷亚洲| 亚洲av综合色区一区| 亚洲经典国产精华液单| 国产男人的电影天堂91| 国产亚洲精品第一综合不卡 | 卡戴珊不雅视频在线播放| 亚洲精品一二三| 日韩中文字幕视频在线看片| 久久精品国产a三级三级三级| 亚洲美女黄色视频免费看| 熟女人妻精品中文字幕| 欧美+日韩+精品| 久久久久精品久久久久真实原创| 久久精品久久久久久久性| 三级国产精品欧美在线观看| av国产久精品久网站免费入址| 一个人看视频在线观看www免费| 欧美国产精品一级二级三级| 秋霞伦理黄片| 大片电影免费在线观看免费| 色婷婷av一区二区三区视频| 亚洲不卡免费看| 国产精品一国产av| 国产无遮挡羞羞视频在线观看| 久久久亚洲精品成人影院| 国产精品久久久久成人av| 精品酒店卫生间| 亚洲国产精品一区三区| 99久久精品一区二区三区| 七月丁香在线播放| 久久久久久久久久久丰满| 免费观看a级毛片全部| 少妇 在线观看| 成年美女黄网站色视频大全免费 | 国产女主播在线喷水免费视频网站| 王馨瑶露胸无遮挡在线观看| 久久女婷五月综合色啪小说| 热re99久久国产66热| 伊人久久精品亚洲午夜| 免费观看性生交大片5| 熟女av电影| 国产亚洲午夜精品一区二区久久| 亚洲丝袜综合中文字幕| 一本—道久久a久久精品蜜桃钙片| 免费黄网站久久成人精品| 男人操女人黄网站| 亚洲国产欧美日韩在线播放| 国产 一区精品| 三级国产精品片| 天天躁夜夜躁狠狠久久av| av不卡在线播放| 久久国产精品大桥未久av| 黄色毛片三级朝国网站| 亚洲中文av在线| 亚洲精品久久午夜乱码| 看免费成人av毛片| 日日爽夜夜爽网站| 久久久久久久久久久免费av| 国产成人a∨麻豆精品| 一个人看视频在线观看www免费| 中国国产av一级| av专区在线播放| 爱豆传媒免费全集在线观看| 精品久久久噜噜| 十分钟在线观看高清视频www| 黑人巨大精品欧美一区二区蜜桃 | 日日摸夜夜添夜夜爱| 人妻一区二区av| 中国国产av一级| 另类亚洲欧美激情| 制服人妻中文乱码| 国产不卡av网站在线观看| 日本猛色少妇xxxxx猛交久久| 精品久久久噜噜| 久久久久视频综合| 在线观看一区二区三区激情| 777米奇影视久久| 一级爰片在线观看| 午夜激情久久久久久久| 在线观看免费视频网站a站| 国产成人精品婷婷| 欧美97在线视频| 啦啦啦在线观看免费高清www| 亚洲欧美色中文字幕在线| 一区二区av电影网| 丰满饥渴人妻一区二区三| 亚洲欧美成人精品一区二区| 精品人妻一区二区三区麻豆| 大话2 男鬼变身卡| 99九九线精品视频在线观看视频| 国产精品麻豆人妻色哟哟久久| 亚洲精品aⅴ在线观看| 丰满迷人的少妇在线观看| 亚洲国产欧美日韩在线播放| 少妇被粗大猛烈的视频| 老熟女久久久| 中国三级夫妇交换| 男女无遮挡免费网站观看| 午夜激情久久久久久久| 80岁老熟妇乱子伦牲交| 国产高清国产精品国产三级| 国产伦理片在线播放av一区| 在线观看美女被高潮喷水网站| 欧美人与善性xxx| 黄色怎么调成土黄色| 亚洲性久久影院| 日韩成人av中文字幕在线观看| 一个人看视频在线观看www免费| 综合色丁香网| 女性生殖器流出的白浆| 一区二区三区乱码不卡18| 一级毛片黄色毛片免费观看视频| 在现免费观看毛片| 久久久亚洲精品成人影院| 少妇人妻精品综合一区二区| 欧美精品一区二区免费开放| 亚洲人成网站在线播| 国产国语露脸激情在线看| 亚洲不卡免费看| 亚洲伊人久久精品综合| 国产成人精品一,二区| 欧美精品国产亚洲| 热99久久久久精品小说推荐| 午夜免费鲁丝| 人人妻人人澡人人爽人人夜夜| videossex国产| 国产av码专区亚洲av| 插阴视频在线观看视频| 男女免费视频国产| 国产精品99久久99久久久不卡 | 亚洲欧美中文字幕日韩二区| 午夜日本视频在线| 91精品国产九色| 欧美日韩精品成人综合77777| 成年av动漫网址| 色婷婷av一区二区三区视频| 精品一区二区三卡| 中文字幕最新亚洲高清| 久久久久久久久久成人| 一二三四中文在线观看免费高清| .国产精品久久| 春色校园在线视频观看| 国产精品 国内视频| 日韩一区二区三区影片| av福利片在线| 国产亚洲av片在线观看秒播厂| 国产男人的电影天堂91| 日韩在线高清观看一区二区三区| 日韩成人av中文字幕在线观看| 精品酒店卫生间| videossex国产| 天堂中文最新版在线下载| a 毛片基地| 日产精品乱码卡一卡2卡三| 蜜桃国产av成人99| 亚洲美女视频黄频| 日韩视频在线欧美| 成人18禁高潮啪啪吃奶动态图 | 欧美xxⅹ黑人| 亚洲av免费高清在线观看| 亚洲av欧美aⅴ国产| 一边摸一边做爽爽视频免费| 看非洲黑人一级黄片| 亚洲av在线观看美女高潮| 黑人高潮一二区| 成人午夜精彩视频在线观看| 日本av手机在线免费观看| 久久精品人人爽人人爽视色| 老司机亚洲免费影院| 国产精品 国内视频| 热re99久久精品国产66热6| 国产极品粉嫩免费观看在线 | 日韩中字成人| 极品少妇高潮喷水抽搐| 国产亚洲午夜精品一区二区久久| 99国产精品免费福利视频| 丝袜美足系列| 边亲边吃奶的免费视频| 免费看av在线观看网站| 九九久久精品国产亚洲av麻豆| 亚洲av成人精品一区久久| 欧美少妇被猛烈插入视频| 亚洲久久久国产精品| 热re99久久国产66热| 制服丝袜香蕉在线| 天堂俺去俺来也www色官网| 最新中文字幕久久久久| 日本91视频免费播放| 亚洲欧洲日产国产| 精品少妇黑人巨大在线播放| av国产久精品久网站免费入址| 久久99一区二区三区| 国产乱来视频区| 少妇丰满av| 国产精品麻豆人妻色哟哟久久| 黑人猛操日本美女一级片| 免费看光身美女| 性色avwww在线观看| 22中文网久久字幕| 中文字幕亚洲精品专区| 精品午夜福利在线看| 亚洲欧美成人精品一区二区| 精品少妇内射三级| 纵有疾风起免费观看全集完整版| 久久99热6这里只有精品| 看非洲黑人一级黄片| 亚洲,欧美,日韩| 亚洲色图综合在线观看| 全区人妻精品视频| 国产乱来视频区| 亚洲欧美成人综合另类久久久| 国产色爽女视频免费观看| av一本久久久久| 久久鲁丝午夜福利片| 我要看黄色一级片免费的| 日韩av不卡免费在线播放| 欧美 日韩 精品 国产| 国产精品国产三级国产av玫瑰| 国产午夜精品一二区理论片| 日本爱情动作片www.在线观看| 在线天堂最新版资源| 水蜜桃什么品种好| 亚洲精品日韩av片在线观看| tube8黄色片| 中国国产av一级| 岛国毛片在线播放| 一本一本综合久久| av播播在线观看一区| 高清在线视频一区二区三区| 国产精品欧美亚洲77777| 精品人妻熟女毛片av久久网站| 日韩免费高清中文字幕av| 亚洲激情五月婷婷啪啪| 亚洲精品乱码久久久v下载方式| 中文精品一卡2卡3卡4更新| 熟女av电影| 欧美日韩国产mv在线观看视频| 国产永久视频网站| 久久久久久久大尺度免费视频| av又黄又爽大尺度在线免费看| 色5月婷婷丁香| av电影中文网址| 日产精品乱码卡一卡2卡三| 老熟女久久久| 国产亚洲最大av| 人人妻人人澡人人爽人人夜夜| 午夜激情av网站| 久久久久精品久久久久真实原创| h视频一区二区三区| 久久久久久久精品精品| 男女免费视频国产| 飞空精品影院首页| 在线观看三级黄色| 日日摸夜夜添夜夜添av毛片| 午夜福利影视在线免费观看| 成人手机av| 国产精品麻豆人妻色哟哟久久| 久久人人爽av亚洲精品天堂| 亚洲第一区二区三区不卡| 亚洲成人一二三区av| 97超视频在线观看视频| 亚洲在久久综合| 老司机影院成人| 婷婷色综合大香蕉| 亚洲精品久久成人aⅴ小说 | 亚洲国产精品999| 日韩欧美精品免费久久| 日产精品乱码卡一卡2卡三| 极品少妇高潮喷水抽搐| 丰满迷人的少妇在线观看| 高清午夜精品一区二区三区| 亚洲成人一二三区av| 精品久久蜜臀av无| 人妻制服诱惑在线中文字幕| 国产精品一区www在线观看| 免费观看a级毛片全部| 青春草视频在线免费观看| 啦啦啦中文免费视频观看日本| 制服丝袜香蕉在线| 女的被弄到高潮叫床怎么办| 大片电影免费在线观看免费| 国产精品.久久久| 久久精品夜色国产| 日韩制服骚丝袜av| a级片在线免费高清观看视频| 国产在视频线精品| 亚洲色图 男人天堂 中文字幕 | 麻豆成人av视频| 亚洲av免费高清在线观看| 欧美人与性动交α欧美精品济南到 | 亚洲欧洲精品一区二区精品久久久 | 我的女老师完整版在线观看| 欧美人与善性xxx| 色婷婷久久久亚洲欧美| 国产欧美日韩综合在线一区二区| 日韩一区二区视频免费看| 99久久综合免费| av在线观看视频网站免费| 一级爰片在线观看| 中文字幕av电影在线播放| 男女边摸边吃奶| 丁香六月天网| 日本爱情动作片www.在线观看| 91久久精品电影网| 亚洲人成网站在线播| 久久ye,这里只有精品| 中文字幕人妻丝袜制服| 久久久久久久大尺度免费视频| 国国产精品蜜臀av免费| 欧美最新免费一区二区三区| 人妻系列 视频| av又黄又爽大尺度在线免费看| 九草在线视频观看| 99热6这里只有精品| 久久久久视频综合| 精品酒店卫生间| 久久精品国产亚洲网站| 精品久久蜜臀av无| 久久久国产精品麻豆| 国产精品一二三区在线看| 亚洲精品乱久久久久久| 大香蕉久久网| 三级国产精品片| 中文精品一卡2卡3卡4更新| 久久久国产精品麻豆| 丰满饥渴人妻一区二区三| 一级片'在线观看视频| 亚洲经典国产精华液单| 免费看av在线观看网站| 国产精品三级大全| 久久国产精品大桥未久av| 亚洲av综合色区一区| √禁漫天堂资源中文www| 国产高清不卡午夜福利| 嘟嘟电影网在线观看| 男女边吃奶边做爰视频| 国产精品久久久久久久电影| 国产欧美日韩一区二区三区在线 | 欧美最新免费一区二区三区| 国产无遮挡羞羞视频在线观看| 一级a做视频免费观看| 在线天堂最新版资源| 亚洲精品国产av蜜桃| 丰满迷人的少妇在线观看| 人成视频在线观看免费观看| 精品亚洲成a人片在线观看|