999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

產氣莢膜梭菌α-毒素基因數字PCR 方法的建立及其在標準物質研制方面的應用

2024-04-01 11:05:19張銘洋張喜悅趙格曲志娜高宏偉孫敏程慧敏徐佳微王君瑋鄒明
山東農業科學 2024年1期

張銘洋 張喜悅 趙格 曲志娜 高宏偉 孫敏 程慧敏 徐佳微 王君瑋 鄒明

摘要:產氣莢膜梭菌(Clostridium perfringens)是一種重要的人獸共患病原菌,α-毒素是其分子生物學檢測的主要靶標。 試驗選取Clper 引物及探針建立微滴式數字PCR(droplet digtial PCR,ddPCR)方法,最佳引物濃度0.85 μmol/ L,探針濃度為0.3 μmol/ L,退火溫度為60 ℃。 利用建立的ddPCR 方法檢測單核細胞增生李斯特菌、沙門氏菌、大腸桿菌、腸球菌、金黃色葡萄球菌、彎曲桿菌的DNA,結果均為陰性,表明該方法具有較好的特異性。 合成產氣莢膜梭菌α-毒素基因,制成DNA 標準品,選取3 個稀釋梯度的標準品進行重復試驗,組內和組間變異系數均小于5%,證明該方法具有較好的重復性。 采用7 個稀釋梯度的標準品進行重復檢測,計算出該方法的檢測限為12.39 copies/ μL,定量限為33.97 copies/ μL。 通過對68 份臨床樣品進行檢測,證實該方法可用于臨床,且用于檢測質粒DNA 時無明顯的基質效應。 應用ddPCR 方法對標準品進行均勻性和穩定性檢驗,并聯合9 家實驗室進行定值,最終獲得了國家市場監督管理總局頒發的標準物質證書(GBW(E)091237)。 α-毒素基因ddPCR 方法的建立和DNA 標準物的研制,為產氣莢膜梭菌分子生物學檢測試劑的標準化奠定了基礎。

關鍵詞:產氣莢膜梭菌;ddPCR;α-毒素基因;標準物質;特異性;靈敏度;可重復性;定值

中圖分類號:S852.61文獻標識號:A文章編號:1001-4942(2024)01-0156-08

微滴式數字PCR ( droplet digital PCR,ddPCR)是近年來迅速發展起來的一項定量分析技術,能夠實現對靶基因的精準定量檢測。 數字PCR 的概念由Vogelstein 等[1] 首次提出,他們在檢測大腸癌患者的K-RAS 基因突變時,進行了微升級的PCR 擴增,分別用不同熒光檢測突變基因和正常基因,通過計算突變基因和正常基因的比例得出突變率。 隨著生物技術的不斷發展與優化,PCR 技術從普通PCR 逐漸發展為熒光定量PCR,目前的ddPCR 屬于第三代PCR。 ddPCR 技術通過對樣品進行微滴化處理,經PCR 擴增后對所有微滴進行逐個檢測,有熒光信號的判讀為1,無熒光信號的判讀為0,再根據泊松分布原理和陽性微滴的數量與比例可得出目的分子的起始拷貝數。 與傳統PCR 及熒光定量PCR 相比,ddPCR具有樣本需求量低、絕對定量以及高靈敏度、高重復性、高耐受性等優勢[2-4] 。 目前,該技術已廣泛應用于轉基因生物、致病菌、病毒、食品安全等領域。

產氣莢膜梭菌(Clostridium perfringens) 廣泛分布于土壤、污水、飼料和動物腸道中,可導致動物的壞死性腸炎、腸毒血癥,以及人類的腹瀉和食物中毒等疾病[1-2] 。 近年來,關于產氣莢膜梭菌分子生物學診斷方法的報道較多,但多為傳統的PCR 和熒光定量PCR 方法,使用的靶基因均為α-毒素基因。 隨著產氣莢膜梭菌分子生物學檢測試劑的廣泛應用,對標準物質研發的需求也越來越強烈。 標準物質可以滿足檢測結果標準化的需求,因此,本研究旨在建立產氣莢膜梭菌α-毒素基因ddPCR 檢測方法,并應用于質粒標準物質的研發。

1 材料與方法

1.1 供試菌株及樣品

產氣莢膜梭菌(ATCC13124)、單核細胞增生李斯特菌(ATCC19111)、沙門氏菌(ATCC14028)、大腸桿菌(ATCC25922)、腸球菌(ATCC29212)、金黃色葡萄球菌(ATCC29213)、彎曲桿菌(ATCC33559)均由國家獸醫微生物菌(毒)種保藏中心動物衛生與流行病學分中心保存。 產氣莢膜梭菌陽性病料由中國動物衛生與流行病學中心致病微生物監測室保存。

1.2 試劑和儀器

ddPCR 探針法超預混液(不含dUTP)、微滴發生油購自北京匯力宏業生物科技有限公司。α-毒素基因、引物和探針序列由生工生物工程(上海)股份有限公司合成。 微滴式數字PCR 檢測系統(伯樂QX200,美國)。

1.3 試驗設計與方法

1.3.1 ddPCR 方法的建立 參照Grau-Roma[5] 、Xu[6] 等的方法分別合成2 對引物和探針(表1)。

分別使用上述2 對引物建立ddPCR 方法并進行反應條件優化。 引物和探針分別設置4 個濃度梯度,其中10 μmol/ L 引物加樣量設置為:1.3、1.5、1.7、1.9 μL;10 μmol/ L 探針加樣量設置為:0.3、0.6、0.9、1.2 μL。 退火溫度設置為56、58、60、62 ℃。通過比較ddPCR 的檢測拷貝數、微滴生成數、微滴分布狀態、微滴熒光信號的強度確定最佳反應條件。

1.3.2 特異性、重復性、檢出限和定量限試驗 分別以產氣莢膜梭菌、單核細胞增生李斯特菌、沙門氏菌、大腸桿菌、腸球菌和金黃色葡萄球菌的DNA 為模板, 用產氣莢膜梭菌α - 毒素基因ddPCR 方法檢測以評價其特異性。 選取3 個濃度梯度(101、102、103 copies/ μL) 的標準品DNA為模板,在間隔第7 天和第14 天分別進行重復試驗,每個濃度設3 個重復,記錄各次試驗的拷貝數,并分析組內和組間的變異系數,以檢測該方法的重復性。 參考?分析方法檢出限和定量限的評估?(GB/ T 27415—2013)[7] 中的方法進行檢出限(IDE)和定量限(IQE)的測定。 以7 個濃度梯度(10-1、100、101、102、103、104、105 copies/ μL)的產氣莢膜梭菌α-毒素基因質粒標準品DNA 為模板,每個梯度重復檢測9 次,用產氣莢膜梭菌α-毒素基因的ddPCR 方法對其擴增,評估方法的檢出限、定量限和動態范圍。

1.3.3 ddPCR 方法的臨床應用及基質效應評估

取實驗室保存的68 份臨床陽性病料,其中產氣莢膜梭菌陽性病料38 份,沙門氏菌、大腸桿菌、腸球菌、單核細胞增生李斯特菌、金黃色葡萄球菌、彎曲桿菌陽性病料各5 份,提取DNA 作為臨床樣本進行產氣莢膜梭菌α-毒素基因ddPCR 方法鑒定。 依據?基質效應與互通性評估指南?(WS/ T356—2011)[8] 中的方法對制備的標準物質進行基質效應與互通性評估驗證。 選取產氣莢膜梭菌臨床陽性樣本20 份,濃度在102 ~ 104 copies/ μL間隨機分布。 準備5 份制備樣本,將制備樣本與臨床樣本隨機穿插排列,使用ddPCR 方法測定所有樣本。

1.3.4 候選標準物質的制備 通過NCBI 網站Blast 比對,合成具有代表性的產氣莢膜梭菌α-toxin 基因(GenBank:AY277724.1) 片段,克隆至pUC57-α-toxin 質粒并轉化E. coli Top 10 感受態細胞。 經擴大培養后提取質粒備用。 應用KpnⅠ(NEB,美國)對質粒進行線性化處理,反應體系為2 μL cutsmart、3 μL KpnⅠ、5 μL ddH2O、10 μL質粒(53.65 ng/ μL),37 ℃水浴3 h。 回收后的線性化質粒,稀釋后分裝作為候選標準物質,-80 ℃保存備用。

1.3.5 候選標準物質的鑒定及檢驗 對產氣莢膜梭菌α-毒素基因質粒進行常規PCR 鑒定、雙酶切和測序鑒定,并對其進行均勻性和穩定性檢驗,鑒定方法參照文獻[9]進行。

1.3.6 定值 選擇9 家具備認可資質的實驗室,對候選標準物質進行聯合定值。 對測得的數據進行正態分布檢驗、組內可疑值檢驗、組間等精度檢驗以及各組平均值之間的t 檢驗,并計算樣品的不確定度,對候選標準物質定值,并通過有證標準物質評估實驗室定值的穩定性[7,9] 。

2 結果與分析

2.1 ddPCR 方法的建立

分別使用2 對引物建立ddPCR 方法,結果Clper 引物探針具有更好的特異性,故在本研究中,使用該引物探針進行目標基因的擴增。 對引物和探針的濃度進行優化,結果10 μmol/ L 引物的用量為1.7 μL、10 μmol/ L 探針的用量為0.6μL,即反應濃度分別為0.85 μmol/ L 和0.3 μmol/ L時,ddPCR 檢測拷貝數最多,微滴分布狀態、微滴熒光信號的強度也均較好(圖1A)。 退火溫度優化結果顯示,60 ℃ 時特異性最好,且信號強度較強(圖1B)。最終建立的ddPCR 方法,主要包括微滴生成、PCR 擴增和微滴讀取3 個步驟。 ①微滴生成。 反應體系為supermix 10 μL、上、下游引物(10 μmol/ L) 各1.7 μL、探針(10 μmol/ L) 0.6μL、模板1 μL、ddH2O 5.0 μL,加入至微滴生成卡的sample 孔;在微滴生成卡的oil 孔中加入70 μL微滴生成油;將生成卡放入微滴生成儀,生成微滴。 ②PCR 擴增。 將生成的微滴轉移至96 孔板,蓋膜后放入熱封儀封膜;將封好膜的96 孔板放入PCR 儀進行擴增;PCR 擴增程序為:50 ℃ 2min;95 ℃ 10 min;95 ℃ 30 s,60 ℃ 1 min,40 個循環;98 ℃ 10 min,4 ℃ 60 min。 ③微滴讀取。 擴增完畢后將96 板轉移至QX200 Droplet Reader 微滴讀取儀上進行微滴讀取。 在以上確定的反應條件下,陰、陽性微滴劃分閾值劃定在熒光強度為2 000的位置。

加樣說明參見表2。

2.2 特異性、重復性、檢出限和定量限結果

分別以單核細胞增生李斯特菌、沙門氏菌、大腸桿菌、腸球菌、彎曲桿菌、金黃色葡萄球菌的DNA 為模板,用建立的產氣莢膜梭菌α-毒素基因ddPCR 方法檢測,結果均為陰性(圖2),表明特異性較好。

選取3 個濃度梯度(101、102、103 copies/ μL)的標準品DNA 為模板,分別在第1、7、14 天進行重復試驗,每個梯度設3 個平行,記錄各次試驗的拷貝數,并統計分析組內和組間的變異系數。 結果3個梯度標準品的組內變異系數為0.73% ~3.17%,組間變異系數為3.36% ~4.27%(表3),重復性和重現性較好。

采用7 個濃度梯度(10-1、100、101、102、103、104copies/ μL 和105 copies/ μL)的產氣莢膜梭菌α-毒素基因質粒DNA 進行檢測。 每個稀釋度3個重復。 結果(圖3)表明,ddPCR 反應產生的微滴數量超過10 000 個,當稀釋度為101 ~104時,產氣莢膜梭菌α-毒素基因拷貝數與濃度之間呈良好的線性關系,可繪制標準曲線(圖4)。 標準曲線方程為y= 4.319x+14.53,R2 = 0.999,計算得出檢出限為12. 39 copies/ μL, 定量限為33. 97copies/ μL。

A08、B08、C08 的DNA 濃度為105 copies/ μL,D 08、E08、F08的DNA濃度為104 copies / μL,G08、H08 的DNA 濃度為103 copies/ μL, B09、C09、D09 的DNA 濃度為102 copies/ μL, E09、F09、G09 的DNA 濃度為101 copies/ μL, H09、D10、E10 的DNA 濃度為100 copies/ μL, F10、G10、H10 的DNA 濃度為10-1 copies/ μL。

2.3 候選標準物質的制備與鑒定

將提取的pUC57-α-toxin 質粒進行雙酶切和PCR 鑒定,結果均可見大小為1 146 bp 的片段,符合預期(圖5A、C)。 電泳鑒定線性化質粒,僅見1 條條帶且大小正確,證實酶切完全(圖5B)。將質粒交由公司進行測序,結果其序列與大腸產氣莢膜梭菌α-toxin 完全一致。 證明產氣莢膜梭菌α-毒素基因質粒DNA 候選標準物質基因序列正確。

2.4 臨床樣品的驗證及基質效應評價

對68 份臨床陽性病料進行檢測,結果38 份產氣莢膜梭菌陽性病料的拷貝數均高于1 000(熒光值均高于2 000),其他非特異性菌株陽性病料的拷貝數均為零(熒光值均低于2 000)。 證明建立的ddPCR 方法可用于臨床樣品的檢測。測定的結果采用線性回歸分析方法進行分析,參考?基質效應與互通性評估指南?[8] 中的方法,計算樣本測定值雙側95%置信區間。 結果表明,標準物質的測定均值在臨床樣本95%置信區間內(圖6),表明標準物質無明顯的基質效應。

2.5 候選標準物質的檢驗

對20 管候選標準物質進行均勻性檢驗,采用單因素方差分析法(F 檢驗法)對結果進行分析,結果F<Fα,表明樣品是均勻的。 經穩定性檢驗,樣品可在-20 ℃條件下穩定儲存18 個月,基因含量未發生顯著變化;將樣品置于4 ℃保存8 d,結果樣品的檢測值無顯著變化,表明標準物質在4 ℃運輸條件下可穩定保持。 模擬反復凍融,結果樣品凍融30 次,基因含量未發生顯著變化。

2.6 定值

聯合9 家實驗室應用ddPCR 方法進行定值。應用SPSS 軟件,對定值數據進行正態分布檢驗,結果顯示數據符合正態分布;應用Microsoft Excel做出正態分布的直方圖,擬合正態曲線(圖7)。對組內數據可疑值檢驗,結果各組所有最大值和最小值之差的絕對值|vi |均小于l(α,n)s,證實組內沒有可疑值。 對各組數據進行等精度檢驗,計算出C 值為0.18,實驗室取α 為0.05,數據組數9,重復測定次數9,查科克倫檢驗臨界值C(α,m,n) 為0.27;結果C≤C(α,m,n) ,表明各組數據為等精度。 對平均值最大組及平均值最小組之間的數據進行t 檢驗,t 值為2.74,查t 值表tα為2.78,t<tα,證明各組平均值之間無顯著差異[7] 。

2.7 不確定度評定

用uA,rel表示相對A 類不確定度,經計算uA,rel為2.1%。 用uB,rel 表示B 類相對不確定度,經計算uB,rel為3.2%。 用ubb,rel表示均勻性差異帶來的不確定度,經計算為2.3%。 用Ults和Usts分別表示長期和短期不確定度,經計算分別為6.60%和3.42%。將各分量合成得到標準物質的合成相對不確定度uCRM,rel 為8.67%,擴展相對不確定度uCRM,rel為17.34%,擴展不確定度uCRM為0.59×103copies/ μL。 綜合確定制備的標準物質的定值為(3.43±0.43)×103 copies/ μL。

2.8 比對與驗證

購置中國計量科學研究院制備的豬繁殖與呼吸障礙綜合征標準物質,用ddPCR 方法對其進行測定。 3 次的結果分別為9.41×108、9.38×108、9.42×108 copies/ μL,經分析,該值符合標準物質說明書上的定值(標準值9.4×108 copies/ μL,擴展不確定度0.9×108 copies/ μL)。 證明本研究建立的ddPCR 方法對標準物質的定值正確。

3 討論與結論

產氣莢膜梭菌分布廣泛,對人類和動物的健康造成嚴重危害。 該菌可產生多種毒素,其中A型產氣莢膜梭菌分泌的α-毒素是動物和人食源性感染的主要毒素型。 研究表明,α-毒素能夠突破血腦屏障、侵害延髓、破壞呼吸系統及心臟的傳導系統,最終引發動物“猝死癥”[10] 。 產氣莢膜梭菌可跟據其產生的4 種主要毒素(α、β、ε、ι)分為5 種毒素型(A、B、C、D、E)。 其中α-毒素是產氣莢膜梭菌分泌的主要致死毒素,是一種多功能磷脂酶,幾乎所有產氣莢膜梭菌都能產生α-毒素。Yonogi 等[11] 建立了用于檢測產氣莢膜梭菌各種毒素的多重PCR 檢測方法;van Asten 等[12] 針對α、β、ε 等毒素建立了多重PCR 檢測方法;Shan ̄non 等[13] 建立了用于檢測城市污水中特定病原菌的熒光定量PCR 檢測方法;Hanabara 等[14] 建立了用于檢測人類糞便中食源性致病菌的熒光定量PCR 檢測方法;另有其他眾多關于檢測產氣莢膜梭菌的報道[15-22] ,他們使用的靶基因均為α-毒素基因。

ddPCR 方法是一種對核酸分子進行絕對定量的檢測技術,已廣泛應用于病原檢測領域,尤其適用于對低拷貝、低含量病原早期的精準檢測,可有效彌補普通熒光定量PCR 方法的不足。ddPCR 將反應體系分散為上萬個微滴,微滴經PCR 擴增后,利用微滴檢測儀對每個微滴進行檢測,最后根據泊松分布原理,直接計算陽性微滴的個數,從而實現直接定量核酸濃度。 本研究建立的產氣莢膜梭菌α-毒素基因ddPCR 方法,具有較好的特異性、靈敏性和重復性,與Li[23-24] 、Cor ̄bisier[25] 、Dupas[26] 等的報道一致。ddPCR 檢測中,DNA 濃度過高會超出儀器的檢測范圍,無法直接用于檢測;而DNA 濃度過低則會對標準物質的穩定性產生不利影響。 為保證標準物質的穩定性,以更好地應用于臨床,本研究制備了濃度相對較低(103 copies/ μL) 的標準物質。 ddPCR 為直接定量方法,其對標準物質的定值結果是全國標準物質管理委員會評價其可信度的主要指標。 本研究利用建立的ddPCR 方法對標準物質進行定值,并獲得了國家市場監督管理總局頒發的標準物質證書(GBW(E) 091237)。該標準物質的研發成功,可促進產氣莢膜梭菌診斷試劑的標準化,減少檢測誤差,確保檢驗人員的檢測能力,規范試劑盒的管理。

相較于其他定值方法,ddPCR 法對技術人員和試劑等要求更高,但隨著技術進步和試劑成本的逐步降低,本研究建立的產氣莢膜梭菌α-毒素基因ddPCR 檢測方法,必將具有廣闊的應用前景。

參 考 文 獻:

[1] Vogelstein B, Kinzler K W. Digital PCR[J]. Proceedings ofthe National Academy of Sciences,1999,96: 9236-9241.

[2] 黃瑾,梁濤波,許恒毅. 數字PCR 在生物學檢測中應用的研究進展[J]. 生命科學,2021,33(2):255-264.

[3] Sacher A G, Paweletz C, Dahlberg S E, et al. Prospective val ̄idation of rapid plasma genotyping for the detection of EGFRand KRAS mutations in advanced lung cancer[J]. JAMA On ̄cology, 2016, 2(8): 1014-1022.

[4] van Ginkel J H, Huibers M M H, van Es R J J, et al. Dropletdigital PCR for detection and quantification of circulating tumorDNA in plasma of head and neck cancer patients[J]. BMCCancer, 2017, 17(1): 1-8.

[5] Grau ̄Roma L, Navarro M, Blatter S, et al. Clostridium per ̄fringens ̄associated necrotic enteritis ̄like disease in coconut lor ̄ikeets (Trichoglossus haematodus)[J]. Veterinary Pathology,2021, 58(2): 423-427.

[6] Xu W P, Wang H R, Liu L X, et al. Prevalence and charac ̄terization of Clostridium perfringens isolated from differentchicken farms in China[J]. Anaerobe, 2021, 72: 102467.

[7] 中華人民共和國國家質量監督檢驗檢疫總局. 分析方法檢出限和定量限的評估:GB/ T 27415—2013[S]. 北京:中國標準出版社,2013.

[8] 中華人民共和國衛生部. 基質效應與互通性評估指南:WS/T 356—2011[S]. 北京:中國標準出版社,2012.

[9] 張喜悅,張銘洋,趙格,等. 彎曲桿菌16S rRNA 基因質粒DNA 標準物質的研制[J]. 中國動物檢疫,2023,40(7):103-109.

[10] 王苗利,蔡玉梅,柴同杰,等. 感染A 型產氣莢膜梭菌小鼠體內α 毒素的分布及定位[J]. 動物醫學進展,2011,32(10):19-23.

[11] Yonogi S, Kanki M, Ohnishi T, et al. Development and appli ̄cation of a multiplex PCR assay for detection of the Clostridiumperfringens enterotoxin ̄encoding genes cpe and becAB [ J].Journal of Microbiological Methods, 2016, 127: 172-175.

[12] van Asten A J A M, van der Wiel C W, Nikolaou G, et al. Amultiplex PCR for toxin typing of Clostridium perfringens iso ̄lates[J]. Veterinary Microbiology (Amsterdam), 2009, 136(3/4): 411-412.v[13] Shannon K E, Lee D Y, Trevors J T, et al. Application of real ̄time quantitative PCR for the detection of selected bacterialpathogens during municipal wastewater treatment[J]. Scienceof The Total Environment,2007, 382(1): 121-129.

[14] Hanabara Y, Ueda Y. A Rapid and simple real ̄time PCR assayfor detecting foodborne pathogenic bacteria in human feces[J].Japanese Journal of Infectious Diseases, 2016, 69(6): 471-476.

[15] 任宏榮,李苗云,朱瑤迪,等. 產氣莢膜梭菌在食品中的危害及其控制研究進展[J]. 食品科學,2021,42(7):352-359.

[16] 孫雨,王曉英,董浩,等. 產氣莢膜梭菌外毒素基因與相關疫苗的研究進展[J]. 中國草食動物科學,2016,36(2):58-62.

[17] 劉立兵,李睿文,陳志敏,等. 產氣莢膜梭菌實時熒光PCR和實時熒光RPA 檢測方法的建立和比較[J]. 食品科學,2020,41(4):268-272.

[18] 石玉玲,曾蘭蘭,陳麗丹,等. 產氣莢膜梭菌實時熒光PCR方法的建立[J]. 生物技術通訊,2010,21(3):389-392.

[19] Brynestad S, Granum P E. Clostridium perfringens and food ̄borne infections[J]. International Journal of Food Microbiolo ̄gy, 2002, 74(3): 195-202.

[20] Maikanov B, Mustafina R, Auteleyeva L, et al. Clostridiumbotulinum and Clostridium perfringens occurrence in Kazakhhoney samples[J]. Toxins, 2019, 11(8): 472.

[21] Brady J, Hernandez ̄Doria J D, Bennett C, et al. Toxinotypingof necrotic enteritis ̄producing and commensal isolates of Clos ̄tridium perfringens from chickens fed organic diets[J]. AvianPathology, 2010, 39(6): 475-481.

[22] 門佩璇,肖宇鋒,張玢. 基于計量學方法分析數字PCR 技術的臨床應用現狀與技術熱點[J]. 生物技術進展,2022,12(4):606-613.

[23] Li J, Li L, Zhang L, et al. Development of a certified genomicDNA reference material for detection and quantification of ge ̄netically modified rice KMD[J]. Analytical and BioanalyticalChemistry, 2020, 412(25): 7007-7016.

[24] Li J, Zhai S S, Gao H F, et al. Development and assessment ofa duplex droplet digital PCR method for quantification of GMrice Kemingdao[J]. Analytical and Bioanalytical Chemistry, 2021,413(16): 4341-4351.

[25] Corbisier P,Pinheiro L,Mazoua S,et al. DNA copy numberconcentration measured by digital and droplet digital quantita ̄tive PCR using certified reference materials[J]. Analytical andBioanalytical Chemistry,2015,407(7): 1831-1840.

[26] Dupas E,Legendre B,Olivier V,et al. Comparison of real ̄time PCR and droplet digital PCR for the detection of Xylellafastidiosa in plants[ J]. Journal of Microbiological Methods,2019, 162: 86-95.

主站蜘蛛池模板: 一级毛片在线播放免费| 无码福利日韩神码福利片| 日日噜噜夜夜狠狠视频| 91精品久久久无码中文字幕vr| 国产精品国产三级国产专业不 | 欧美日韩精品一区二区在线线| 久久精品电影| 成人一级免费视频| a欧美在线| 98超碰在线观看| 国产在线一区视频| 亚卅精品无码久久毛片乌克兰| 欧美一级大片在线观看| 日韩激情成人| 高潮爽到爆的喷水女主播视频| 国产成人高清精品免费| 呦女精品网站| 无码中文AⅤ在线观看| 国产极品嫩模在线观看91| www.91中文字幕| av免费在线观看美女叉开腿| 亚洲娇小与黑人巨大交| 99re免费视频| 美女被操91视频| 成人午夜视频免费看欧美| www.av男人.com| 美女被操黄色视频网站| 中文字幕 欧美日韩| 国产99久久亚洲综合精品西瓜tv| 久久国语对白| 国产毛片久久国产| 欧美中文字幕第一页线路一| 九九热精品在线视频| 中文字幕啪啪| 91精品情国产情侣高潮对白蜜| 少妇精品在线| 欧美午夜性视频| 毛片一区二区在线看| 免费黄色国产视频| 欧美一级爱操视频| 国产91小视频| 国产99精品久久| 99尹人香蕉国产免费天天拍| 精品成人一区二区三区电影| 动漫精品中文字幕无码| 特级aaaaaaaaa毛片免费视频 | 婷婷色婷婷| 欧美日韩高清| 欧美精品伊人久久| 日韩毛片基地| 国产特一级毛片| 青草免费在线观看| 色综合狠狠操| 成人午夜精品一级毛片| 久久久久久久97| 看你懂的巨臀中文字幕一区二区| 精品伊人久久久久7777人| 国产成人亚洲精品色欲AV| 亚洲国语自产一区第二页| 欧美在线网| 日本尹人综合香蕉在线观看| 中文字幕亚洲精品2页| 99久久精品国产综合婷婷| 在线播放国产一区| 日韩视频免费| 午夜精品国产自在| h网站在线播放| 国产女同自拍视频| 国产在线啪| 亚洲无限乱码一二三四区| 国产第二十一页| 亚洲第一福利视频导航| 青青草a国产免费观看| 综合五月天网| 欧美成人午夜视频免看| 国产一区二区精品福利| 无码免费试看| 99热亚洲精品6码| 蜜桃视频一区二区| 亚洲大尺码专区影院| 美女免费精品高清毛片在线视| 国产成人综合在线视频|