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

反相離子對色譜法對寡核苷酸藥物分析的進展

2015-07-07 15:15:58楊洪淼范慧紅
中國生化藥物雜志 2015年7期
關鍵詞:分析

楊洪淼,范慧紅

(中國食品藥品檢定研究院,北京 100050)

?

反相離子對色譜法對寡核苷酸藥物分析的進展

楊洪淼,范慧紅Δ

(中國食品藥品檢定研究院,北京 100050)

寡核苷酸藥物是近年來藥物研發領域中的熱點之一,相應分析方法的研究也被廣泛開展,其中反相離子對色譜法以分離度高、易與MS對接等優勢成為寡核苷酸分析中最常用的方法。本文綜述了反相離子對色譜法在寡核苷酸分析方面的新進展,以下將從色譜柱、流動相,分離規律的探索,MS檢測器方面進行系統介紹。

離子對試劑;高效液相色譜;寡核苷酸藥物;綜述

寡核苷酸藥物是指含有10~30個堿基的RNA或DNA分子,及其類似物,可以在基因水平上干擾致病蛋白質的產生,從而達到治療疾病的目的。目前藥物研發的方向根據作用機理主要分為:反義寡核苷酸(antisense oligonucleotide ,asON),是指含有15~30個堿基的核酸分子或其類似物,能與特定的DNA、RNA以堿基互補配對的方式結合并阻止其轉錄與翻譯,這是一類經典的寡核苷酸藥物,通常會進行結構修飾以增加對體內核酸酶的穩定性,最常見的修飾為將磷酸二酯鍵中的一個非橋氧原子用硫原子取代,被稱為硫代寡核苷酸[1];另一類是小分子干擾RNA(small-interfering RNAs,siRNAs),通常為含有21~23 個堿基的單鏈或雙鏈RNA,可以引發RNA 干擾現象,即指由于目標mRNA降解或翻譯抑制所引起的轉錄后基因沉默現象[2];此外也有核酸適配體[3]等研究方向。至今FDA已批準了數個寡核苷酸藥物的上市,如2013年獲批上市的米泊美生(mipomersen),而處于臨床研發階段的藥物數量更多。由此可見,開展對寡核苷酸分析方法的研究工作非常必要。

寡核苷酸藥物的分子量較大,通常為6000~8000,同時在水中幾乎以陰離子形式存在,會形成一個負電荷聚集的分子體[4]。根據這樣的結構特征,常用的分析方法有反相離子對色譜、親水色譜、離子交換色譜、毛細管電泳和PCR等方法[5-8]。其中反相離子對色譜以分離度高、結果穩定、操作簡便和易與MS檢測器對接等優勢,在近年的文獻中越來越多被使用到。本文通過列舉一些寡核苷酸藥物分析方法開發的實例,和對寡核苷酸類物質分析方法研究的新進展等內容的介紹,希望對今后這一類新型藥物分析方法的開發提供參考。

1 色譜柱

使用離子對色譜方法對寡核苷酸進行分析時,最常用的色譜柱即C18柱。 Fountain等[9]使用小粒徑的C18柱(4.6 mm×75 mm,2.5 μm)可同時將5~30個堿基的寡核苷酸實現N與N-1序列的分離。Studzińska等[10]考察了6根UPLC色譜柱對含有細微差別的多個寡核苷酸的分離效果,以C18柱和苯基柱為最優,并同時考察了離子對試劑的種類和濃度對分離的影響。小粒徑色譜柱可以縮短寡核苷酸分子在固定相內擴散的路徑,從而抑制峰形展寬,更有利于分離。此外采用高分子骨架的色譜柱可以對偏堿性的流動相及較高的柱溫條件更加耐受。

近年來一些新型的色譜柱也被用于寡核苷酸的分析中。Sharma等[11],Xiong等[12]嘗試采用聚苯乙烯-二乙烯基苯(PS-DVB)共聚物為固定相的整體柱(monolithic column)可以將修飾基團連接到不同堿基位點的寡核苷酸同分異構體進行成功分離。Biba等[13],Zimmermann等[14]使用一種結合離子交換與反相2種分離模式的色譜柱對20個堿基左右的寡核苷酸進行分析,分離過程同時利用了寡核苷酸的負電性和疏水性。文獻中發現若使用離子對試劑作為流動相,分離仍以反相模式為主導;若使用離子交換色譜的條件,如NaCl溶液梯度洗脫,所得色譜圖與只使用離子交換色譜柱的結果差異非常大,說明此時分離過程為混合模式,并且使用NaCl或NaBr作為流動相可以顯著改善寡核苷酸的峰形。Studzińska等[15]在實驗室制備中作為新型固定相將硅膠基質上同時鍵合1個疏水基團和1個親水基團,根據疏水基團的不同分為酰胺型和膽固醇型,此種類型固定相首次應用于寡核苷酸的色譜分析中,結果表明它們的分離效果與C18柱相當,但分離時間會更短。

2 流動相

由于寡核苷酸具有強烈的負電性,在常規的反相色譜條件下保留較弱,所以需要在流動相中添加離子對試劑。強離子對試劑如季銨鹽會使寡核苷酸產生強保留,但同時會削弱不同堿基疏水性差異對保留的影響,反而不利于寡核苷酸藥物與雜質的分離[16]。而三乙胺(triethylamine,TEA)體系為較弱的陽離子對體系,在寡核苷酸分析中最為常用。三乙胺/醋酸(triethylamine/acetic acid,TEAA)系統是最先被應用的,通常30個堿基以內未修飾的寡核苷酸均可以達到N與N-1序列的分離[17]。但硫代寡核苷酸由于硫原子取代了氧原子,使得每一個磷原子形成為一個手性中心,理論上N個堿基長度的硫代寡核苷酸可以形成2n-1個非對映異構體,這些異構體在色譜保留行為上的細小差異最終造成嚴重的峰形展寬[4]。Apffel等[18]開創性地引入六氟異丙醇(HFIP)至離子對系統中,該試劑可以有效地抑制硫代寡核苷酸的峰形展寬,從而提高分離度。因為流動相中HFIP會顯著降低TEA在水中的溶解度,所以推測HFIP會迫使TEA更多地貼合在硅膠柱的表面,使硅膠表面形成一層穩定的TEA+膜,此時的分離模式更傾向于離子交換,從而抑制對映異構體由于疏水性差異造成的峰形展寬[4]。

Cantide、CT102和流感泰得[19-21]均為研發中的硫代反義寡核苷酸藥物,分別為20、20和13個堿基長度。它們均可以應用C18柱、離子對試劑TEA/HFIP的方法使主成分與雜質達到良好地分離。硫代寡核苷酸的合成雜質主要由短序列和硫代失敗序列構成,其中以5’N-1序列和一個氧未被硫代(P=S/O)的N序列最多,由于雜質與主成分結構極為相似,分離難度很大。之前通常采用陰離子交換色譜法分析P=S/O雜質及其他雜質,但離子交換色譜對N-1雜質無能為力,需再采用毛細管電泳方法分析N-1雜質。而反相離子對色譜法的優勢是可以同時分離N-1雜質和P=S/O雜質,簡化了分析步驟。根據文獻中優化的結果,通常TEA與HFIP的濃度范圍分別為10~40 mmol/L和100~400 mmol/L的時候,更利于20個堿基長度左右的硫代寡核苷酸與雜質的分離。但不同序列的硫代寡核苷酸適用的TEA/HFIP的比例各不相同,需要通過實驗優化。此外由于HFIP在乙腈中的溶解度很小,所以流動相中一般使用甲醇作為有機相。

除此之外,N,N-二甲基丁胺/醋酸(DMBAA)、四乙基溴化銨/醋酸(TEAB)、醋酸正己基銨(HAA)等離子對系統也被應用于寡核苷酸的反相色譜分析中[10,22-23]。Studzińska等[10]比較了DMBAA與TEAA對未修飾寡核苷酸分析的差別,結果顯示DMBAA雖增加保留時間但會降低分離度。Levin等[24]嘗試開發多種離子對試劑聯合應用于siRNA的分析,結果顯示TEA/丙胺(PA)的組合可使整體分離度提高,并可明顯改善峰形。

3 分離規律的探索

3.1 保留時間的預測 由于組成寡核苷酸的核苷的疏水性各不相同,A和T的疏水性要大于C和G,所以長度或序列不相同的寡核苷酸在反相色譜中的保留時間也不同。Gilar等[25]建立的數學模型可以預測10~60個堿基內不同長度與序列的寡核苷酸的保留行為,采用的色譜條件為100 mmol/L TEAA離子對試劑,pH為7。并進一步設計39個寡核苷酸對此模型進行驗證,結果顯示預測值與實測值的相關性很好,r2為0.946。此研究可用于合成寡核苷酸的純化中,通過該數學模型即可得到最佳的梯度洗脫的初始有機相比例,節省優化色譜條件的時間。在此基礎上,Sturm等[26]也建立了對寡核苷酸保留時間的預測模型,并且有2個進步,一是使用支持向量回歸(SVR)的方法代替單一的線性回歸或對數回歸,SVR可以模擬非線性關系;二是將寡核苷酸二級結構的信息引入模型中,從而提高預測的準確性。Lei等[27]建立的定量結構-保留時間關系(QSRR)模型可以將寡核苷酸的結構信息進行量化的輸入,并且可以很好地預測多個柱溫下的保留時間。而Studzińska等[28]進一步將QSRR模型成功應用于寡核苷酸在酰胺柱和膽固醇柱上的保留行為的預測。

3.2 流速與柱溫對分離的影響 由于寡核苷酸的傳質速率較低,根據Van Deemter方程,影響塔板高度最主要的參數項為傳質阻抗項,所以影響傳質速度的因素都將顯著影響寡核苷酸的分離。基于此理論,降低流速和提高柱溫被認為是提高寡核苷酸分離度的有效手段,如流速0.5 mL/min和柱溫50 ℃為常見色譜條件。駱雪芳等[29-30]對此進行了深入的研究,設計19、20、21個堿基的不同長度的3條寡核苷酸鏈,和相同長度(20個堿基)但序列上有細微差異的3條寡核苷酸鏈。通過實驗考察流速和柱溫的變化對分離度的影響。結果顯示,對于堿基數目較大、以長度差異為主要特征的寡核苷酸序列,增加溫度是改善分離的有效途徑。而對于同長序列,提高溫度的影響各不相同,由于疏水性強的堿基處于末端對保留貢獻大,若處于中間位置則取決于其在空間結構中的位置,所以疏水堿基處于末端等優勢位置時,隨溫度的改變保留影響不大,處于中間等劣勢位置容易受到溫度的影響,產生較大的保留改變。對于流速的影響,所有寡核苷酸鏈隨流速增加,保留均迅速下降,分離度也呈下降趨勢,但基于長度差異的樣品的分離度下降得更明顯,此類樣品應注重流速的優化,但過低的流速會使樣品縱向擴散增加,同樣發生分離度下降情況。

4 MS檢測器

反相離子對色譜(IP-HPLC)通過添加揮發性的離子對試劑可以與質譜更好地對接,而質譜由于其高靈敏度和快速定性等特性成為寡核苷酸藥物分析中的重要分析手段。

反相離子對色譜與ESI-MS聯用在藥物代謝研究中發揮著重要作用[31-33]。Wei等[31]采用IP-HPLC-MS/MS的方法對硫代寡核苷酸藥物GTI-2040的藥代產物進行鑒定。IP-HPLC的條件可以得到N與N-1序列的分離,高分離度是藥代產物鑒定的前提。文中使用SOS寡核苷酸序列分析軟件(the simple oligonucleotide sequencer software),共鑒別了5種代謝產物。Dai等[32]也采用上述IP-HPLC-MS/MS的方法對硫代寡核苷酸藥物G3139的藥代產物進行鑒定,并建立了G3139和3個代謝產物的定量分析方法,最小定量限(LOQ)為17.6 nmol/mL。

使用MS/MS對短鏈的寡核苷酸進行測序已是較為成熟的方法[34],但經過化學修飾的寡核苷酸,裂解行為發生改變,測序工作則更具挑戰。Kretschmer等[35]采用IP-HPLC-MS/MS的方法對化學修飾的siRNA進行測序,待測物包含2’-O-甲基化和硫代2種修飾方式,該方法的建立可以很好地對寡核苷酸合成的過程進行實時監測。

HFIP可以顯著提高寡核苷酸的離子化效率,提高質譜檢測的靈敏度。近年來對適用于質譜的離子對試劑的研究也很多[36-39]。McGinnis等[36]發現小于100 mmol/L的HFIP對質譜的信號增強作用更明顯,比較了13種烷基胺分別與HFIP組成離子對試劑對質譜信號的增強作用,得出待測物寡核苷酸的疏水性與離子對試劑的選擇關系很大,10 mmol/L二異丙胺(DIPA)和25 mmol/L HFIP適用于大部分親水性的siRNA,15 mmol/L三丙胺(TPA)和50 mmol/L HFIP適用于中等疏水性的未修飾DNA片段,10 mmol/L二異丙基乙胺(DIEA)和50 mmol/L HFIP適用于疏水性更強的硫代寡核苷酸。

此外也有多篇文獻是與ICP-MS聯用或與MALDI-TOF-MS聯用的應用實例[40-42]。Vallone等[40]采用IP-HPLC的方法將短鏈DNA進行分離和收集,后采用MALDI-TOF-MS進行鑒定,液相步驟中所采用的離子對試劑為TEAA。Studzińska等[42]通過IP-HPLC與ICP-MS聯用的方式建立了血漿中硫代寡核酸的含量測定方法,并進行了方法學驗證。離子對試劑采用TEA/HFIP,但發現高濃度的HFIP會產生較高的背景干擾。

5 展望

綜上所述,反相離子對色譜法是寡核苷酸分析的一種非常有效的手段,今后的研究方向一方面為更多新型的固定相被應用,如兼具反相與離子交換模式的固定相,以及更多的離子對試劑被開發,尤其是可以提高MS檢測靈敏度的離子對試劑;另一方面反相離子對色譜法分析寡核苷酸的分離機理及分析規律也將更多地被揭示。

[1] 張嵬.硫代反義寡核苷酸質量控制方法的研究進展[J].國外醫學藥學分冊,2005,32(2):132-135.

[2] 于鐵成,王巖,劉建國,等.小分子干擾RNA和反義寡核苷酸技術的比較及其應用意義[J].中國組織工程研究與臨床康復,2007,2(36):7262-7266.

[3] Lin CY,Huang Z,Jaremko W,et al.High-performance liquid chromatography purification of chemically modified RNA aptamers[J].Anal Biochem,2014(449):106-108.

[4] McGinnis AC,Chen B,Bartlett MG.Chromatographic methods for the determination of therapeutic oligonucleotides[J].J Chromatogr B,2012(883-884):76-94.

[6] Cook K,Thayer J.Advantages of ion-exchange chromatography for oligonucleotide analysis[J].Bioanalysis,2011,3(10):1109-1120.

[7] Chen B,Bartlett M.G.Determination of therapeutic oligonucleotides using capillary gel electrophoresis[J].Biomed Chromatogr,2012,26 (4):409-418.

[8] Yu RZ,Geary RS,Levin AA.Application of novel quantitative bioanalytical methods for pharmacokinetic and pharmacokinetic / pharmacodynamic assessments of antisense oligonucleotides[J].Curr Opin Drug Discov,2004,7(2):195-203.

[9] Fountain KJ,Gilar M,Gebler JC.Analysis of native and chemically modified oligonucleotides by tandem ion-pair reversed-phase high-performance liquid chromatography/electrospray ionization mass spectrometry[J].Rapid Commun Mass Spectrom,2003,17(7):646-653.

[10] Studzińska S,Buszewski B.Evaluation of ultrahigh-performance liquid chromatography columns for the analysis of unmodified and antisense oligonucleotides[J].Anal Bioanal Chem,2014,406(28):7127-7136.

[11] Sharma VK,Glick J,Vouros P.Reversed-phase ion-pair liquid chromatography electrospray ionization tandem mass spectrometry for separation,sequencing and mapping of sites of base modification of isomeric oligonucleotide adducts using monolithic column[J].J Chromatogr A,2012,1245:65-74.

[12] Xiong W,Glick J,Vouros P,et al.Separation and sequencing of isomeric oligonucleotide adducts using monolithic columns by ion-pair reversed-phase nano-HPLC coupled to ion trap mass spectrometry[J].Anal Chem,2007,79(14):5312-5321.

[13] Biba M,Jiang E,Welch CJ,et al.Factors influencing the separation of oligonucleotides using reversed-phase /ion-exchange mixed-mode high performance liquid chromatography columns[J].J Chromatogr A,2013,1304:69-77.

[14] Zimmermann A,Greco R,L?mmerhofer M,et al.Synthetic oligonucleotide separations by mixed-mode reversed phase /weak anion-exchange liquid chromatography[J].J Chromatogr A,2014(1354):43-55.

[15] Studzińska S,Pietrzak L,Buszewski B.The Effects of Stationary Phases on Retention and Selectivity of Oligonucleotides in IP-RP-HPLC[J].Chromatographia,2014,77(23-24):1589-1596.

[16] Lloyd LL,Millichip MI,Mapp KJ.Rigid polymerics:the future of oligonucleotide analysis and purification[J].J Chromatogr A,2003,1009(1-2):223-230.

[17] Gilar M.Analysis and purification of synthetic oligonucleotides by reversed-phase high-performance liquid chromatography with photodiode array and mass spectrometry detection[J].Anal Biochem,2001,298(2):196-206.

[18] Apffel A,Chakel JA,Fisher S,et al.Analysis of Oligonucleotides by HPLC-Electrospray Ionization Mass Spectrometry[J].Anal Chem,1997,69(7):1320-1325.

[19] Shu L,Lu DD,Wang SQ,et.al.An Improved Ion-Pair Reversed Phase LC Method for Analysis of Major Impurities of Phosphorothioate Oligonucleotide Cantide[J].Chromatographia,2010,72(3-4):215-223.

[20] 賈巧云,鄧新秀,王升啟,等.反相離子對高效液相色譜法分析硫代反義寡核苷酸CT102相關雜質[J].藥物分析雜質,2013,33(7):1157-1162.

[21] 王穎,魯丹丹,王升啟,等.離子對反相液相色譜法分析硫代寡核苷酸藥物流感泰得的有關物質[J].分析化學,2011,39(12):1811-1816.

[22] Huber C,Oberacher H.Analysis of nucleic acids by on-line liquid chromatography mass spectrometry[J].Mass Spectrom Rev,2001,20 (5):310-343.

[23] McCarthy SM,Gilar M,Gebler J.Reversed-phase ion-pair liquid chromatography analysis and purification of small interfering RNA[J].Anal Biochem,2009,390 (2):181-188.

[24] Levin DS,Shepperd BT,Gruenloh CJ.Combining ion pairing agents for enhanced analysis of oligonucleotide therapeutics by reversed phase ion pairing ultra performance liquid chromatography (UPLC)[J].J Chromatogr B Analyt Technol Biomed Life Sci,2011,879(19):1587-1595.

[25] Gilar M,Fountain KJ,Gebler JC.Ion-pair reversed-phase high-performance liquid chromatography analysis of oligonucleotides:Retention prediction[J].J Chromatogr A,2002,958(1-2):167-182.

[26] Sturm M,Quinten S,Kohlbacher O,et al.A statistical learning approach to the modeling of chromatographic retention of oligonucleotides incorporating sequence and secondary structure data[J].Nucleic Acids Research,2007,35(12):4195-4202.

[27] Lei B,Li S,Yao X,et al.Novel approaches for retention time prediction of oligonucleotides in ion-pair reversed-phase high-performance liquid chromatography[J].J Chromatogr A,2009,1216(20):4434-4439.

[28] Studzińska S,Buszewski B.Different approaches to quantitative structure-retention relationships in the prediction of oligonucleotide retention[J].J Sep Sci,2015,38(12):2076-2084.

[29] 駱雪芳,尹小英,胡育筑,等.柱溫和流速對反相離子對色譜法分析寡核苷酸的影響[J].分析化學,2008,36(5):627-631.

[30] 駱雪芳,陳蓉,胡育筑,等.反相離子對色譜法分析3’,5’反向寡核苷酸[J].色譜,2007,25(6):814-819.

[31] Wei XH,Dai GW,Chan KK,et al.Metabolism of GTI-2040,a Phosphorothioate Oligonucleotide Antisense,Using Ion-pair Reversed Phase High Performance Liquid Chromatography (HPLC) Coupled With Electrospray Ion-trap Mass Spectrometry[J].AAPS J,2006,8(4):743-755.

[32] Dai GW,Wei XH,Chan KK,et al.Characterization and quantification of Bcl-2 antisense G3139 and metabolites in plasma and urine by ion-pair reversed phase HPLC coupled with electrospray ion-trap mass spectrometry[J].J Chromatogr B Analyt Technol Biomed Life Sci,2005,825(2):201-213.

[33] Erb R,Leithner K,Oberacher H.Phosphorothioate oligonucleotide quantification byμ-liquid chromatography-mass spectrometry[J].AAPS J,2012,14(4):728-737.

[34] Oberacher H,Pitterl F.On the use of ESI-Qq-TOF-MS/MS for the comparative sequencing of nucleic acids[J].Biopolymers,2009,91(6):401-409.

[35] Kretschmer M,Lavine G,Manoharan M,et al.An automated algorithm for sequence confirmation of chemically modified oligonucleotides by tandem mass spectrometry[J].Anal Biochem,2010,405(2):213-223.

[36] McGinnis AC,Grubb EC,Bartlett MG.Systematic optimization of ion-pairing agents and hexafluoroisopropanol for enhanced electrospray ionization mass spectrometry of oligonucleotides[J].Rapid Commun Mass Spectrom,2013,27(33):2655-2664.

[37] Gong L,McCullagh JS.Comparing ion-pairing reagents and sample dissolution solvents for ion-pairing reversed-phase liquid chromatography/electrospray ionization mass spectrometry analysis of oligonucleotides[J].Rapid Commun Mass Spectrom,2014,28(4):339-350.

[38] Erb R,Oberacher H.Comparison of mobile-phase systems commonly applied in liquid chromatography-mass spectrometry of nucleic acids[J].Electrophoresis,2014,35(9):1226-1235.

[39] Chen B,Bartlett MG.Evaluation of mobile phase composition for enhancing sensitivity of targeted quantification of oligonucleotides using ultra-high performance liquid chromatography and mass spectrometry:application to phosphorothioate deoxyribonucleic acid[J].J Chromatogr A,2013,1288:73-81.

[40] Vallone PM,Devaney JM,Butler JM.A strategy for examining complex mixtures of deoxyoligonucleotides using ion-pair-reverse-phase high-performance liquid chromatography,matrix-assisted laser desorption ionization time-of-flight mass spectrometry,and informatics[J].Anal Biochem,2002,304(2):257-265.

[41] Cai HL,Lu DD,Wang SQ,et al.The analysis of major impurities of lipophilic-conjugated phosphorothioate oligonucleotides by ion-pair reversed-phase HPLC combined with MALDI-TOF-MS[J].Anal Bioanal Chem,2012,403(5):1333-1342.

[42] Studzińska S,Mounicou S,Buszewski B,et al.New approach to the determination phosphorothioate oligonucleotides by ultra high performance liquid chromatography coupled with inductively coupled plasma mass spectrometry[J].Anal Chim Acta,2015(855):13-20.

(編校:王儼儼)

Progress of ion-pairing reverse liquid chromatography for analysis of oligonucleotides

YANG Hong-miao, FAN Hui-hongΔ

(National Institutes for Food and Drug Control, Beijing 100050)

Therapeutic oligonucleotides have been one of the hottest areas in R&D of drugs in recent years.The research of analysis method has been widely carried out.Ion-pairing reverse liquid chromatography is the most commonly used method for good resolution and easy docking with MS.A review is presented here on the progress of oligonucleotides analysis with ion-pairing reverse liquid chromatography.The column, mobile phase, MS detection, mechanism of analysis will be discussed.

ion-pairing agents; HPLC; therapeutic oligonucleotides; review

楊洪淼,女,碩士,助理研究員,研究方向:生化藥品的質量分析和質量標準的研究,E-mail:yanghongmiao@163.com;范慧紅,通訊作者,女,博士,研究員,研究方向:生化藥品的質量控制和質量標準的研究,E-mail:shenghuayaoshi@126.com。

R917

A

1005-1678(2015)07-0161-04

猜你喜歡
分析
禽大腸桿菌病的分析、診斷和防治
隱蔽失效適航要求符合性驗證分析
電力系統不平衡分析
電子制作(2018年18期)2018-11-14 01:48:24
電力系統及其自動化發展趨勢分析
經濟危機下的均衡與非均衡分析
對計劃生育必要性以及其貫徹實施的分析
現代農業(2016年5期)2016-02-28 18:42:46
GB/T 7714-2015 與GB/T 7714-2005對比分析
出版與印刷(2016年3期)2016-02-02 01:20:11
中西醫結合治療抑郁癥100例分析
偽造有價證券罪立法比較分析
在線教育與MOOC的比較分析
主站蜘蛛池模板: 国产精品私拍在线爆乳| 香蕉99国内自产自拍视频| 国产精品久久自在自线观看| 精品国产aⅴ一区二区三区| 狂欢视频在线观看不卡| 一级全免费视频播放| 中文字幕伦视频| 日韩一级毛一欧美一国产| 久久无码av一区二区三区| 欧美日韩亚洲综合在线观看 | 国产视频入口| 毛片三级在线观看| 老司机aⅴ在线精品导航| 女人av社区男人的天堂| 久久性视频| 五月激激激综合网色播免费| 久热re国产手机在线观看| 91小视频在线观看| aⅴ免费在线观看| 中文无码伦av中文字幕| 免费国产好深啊好涨好硬视频| 久久久国产精品无码专区| 伊人久久影视| 毛片久久久| 亚洲天堂福利视频| 在线国产91| 五月婷婷欧美| 亚洲精品在线影院| 国产视频你懂得| 精品国产成人国产在线| 成人免费网站在线观看| 国产极品嫩模在线观看91| 国产成人福利在线| 沈阳少妇高潮在线| 亚洲天堂久久新| 97免费在线观看视频| 日韩国产一区二区三区无码| 亚洲an第二区国产精品| 久久网综合| 无码人妻热线精品视频| 天天综合色网| 亚洲无码37.| av在线手机播放| 99热6这里只有精品| 成人欧美在线观看| 国产成人调教在线视频| 91精品国产自产在线观看| 天天综合天天综合| 国产高清不卡视频| 尤物精品视频一区二区三区| 亚洲综合一区国产精品| 97综合久久| 少妇精品久久久一区二区三区| 日韩不卡高清视频| 青青草国产免费国产| 亚洲男人天堂网址| 国产91透明丝袜美腿在线| 97久久免费视频| 国产交换配偶在线视频| 免费国产福利| 影音先锋亚洲无码| 精品国产毛片| 国产成人亚洲精品蜜芽影院| 婷婷激情五月网| 91精品aⅴ无码中文字字幕蜜桃| 欧美日韩免费在线视频| 激情在线网| 国产视频入口| 国产精品自在在线午夜| 男女男精品视频| 中文字幕 欧美日韩| 国产精品污污在线观看网站| 一本久道热中字伊人| 国产精品女主播| 天天色综网| 欧美自拍另类欧美综合图区| 国产精品视频第一专区| 国产又爽又黄无遮挡免费观看| 欧美自慰一级看片免费| 麻豆a级片| 性视频久久| 国产色婷婷|