胰腺癌血清學標記物在胰腺癌初期診斷中的研究進展

曾 強，鞏 燕

解放軍總醫(yī)院 健康管理研究院，北京 100853

胰腺癌血清學標記物在胰腺癌初期診斷中的研究進展

曾 強，鞏 燕

解放軍總醫(yī)院 健康管理研究院，北京 100853

胰腺癌是目前最難診斷和治療的惡性腫瘤之一，研究開發(fā)方便實用的血清標記物篩查早期胰腺癌有重要的臨床意義。目前，胰腺癌血清學標記物的研究主要包括胰腺癌特異表達的分子、基因突變、異常甲基化DNA、microRNA和循環(huán)腫瘤細胞。聯(lián)合多種血清學標記物有助于其早期診斷，未來的研究應側(cè)重于篩選胰腺癌高危人群，形成高度敏感性和特異性的標記物，提高胰腺癌早期診斷水平。

胰腺癌；血清標記物；早期診斷

胰腺癌是目前最難診斷和治療的惡性腫瘤之一，80%以上的胰腺癌確診時已處于晚期，其中位生存期在8個月左右[1-4]。開發(fā)方便實用的血清標記物篩查早期胰腺癌有重要的臨床意義[5]。目前，胰腺癌血清學標記物的研究包括胰腺癌特異表達的蛋白質(zhì)分子、基因突變、異常甲基化DNA、microRNA和循環(huán)腫瘤細胞等，本文就其研究進展做一綜述[6-9]。

1 胰腺癌特異的蛋白質(zhì)分子或蛋白譜

1.1 抗黏蛋白1抗體(PAM4-based MUC1) PAM4是一種由RIP1小鼠胰腺癌移植瘤中獲得的抗黏蛋白1(mucin-1，MUC-1)的單克隆抗體，MUC-1是一種高度糖基化的黏蛋白，可以促進腫瘤細胞的侵襲和轉(zhuǎn)移，進而影響胰腺癌患者的生存時間[10]。Gold等[11-13]采用EILSA法測定胰腺癌、慢性胰腺炎及健康人群血清MUC-1表達水平，發(fā)現(xiàn)PAM4診斷胰腺癌的敏感度為81%，特異度為95%，對1期和2期的敏感度分別為62%和86%，并且診斷特異性優(yōu)于糖類抗原19-9(carbohydrate antigen19-9，CA19-9)(85% vs 68%)。通過免疫組化方法檢測非侵襲性胰腺腫瘤，如胰管上皮內(nèi)瘤變及導管內(nèi)乳頭狀黏液瘤的敏感性為87.3%，特異性為92.1%，提示該抗體所針對的抗原通常在胰腺癌發(fā)生的早期表達，對其進行檢測有助于胰腺癌的早期發(fā)現(xiàn)和診斷。

1.2 高流動蛋白1(extracellular high mobility group box-1，HMGB1) HMGB1被認為是一種核蛋白，其細胞外形式類似細胞因子發(fā)揮作用，參與如炎癥、細胞遷移、組織再生等許多重要的病理過程，有助于腫瘤的生長和侵襲[14-15]。比較正常、慢性胰腺炎和胰腺癌患者血清HMGB1水平，發(fā)現(xiàn)它與胰腺癌及分期有關(guān)，ROC曲線和Logistic回歸分析證實HMGB1作為預測胰腺癌單一或多重生物標記物之一，其靈敏度/特異性優(yōu)于CA 19-9或癌胚抗原(carcinoembryonic antigen，CEA)，Kaplan-Meier生存分析顯示高水平HMGB1的胰腺癌患者(＞30 ng/ml，中位生存期192 d)較低水平HMGB1患者(≤30 ng/ml，中位生存期514 d)預后差，Cox比例風險模型顯示，高血清HMGB1組較低血清HMGB1組患胰腺癌相對危險比為3.077。與CA19-9和CEA進行比較，HMGB1是胰腺癌診斷和預后評價更為理想的標記物[16]。

1.3 胰島素樣生長因子結(jié)合蛋白2(insulin-like growth factor binding protein 2，IGFBP2)和間皮素(mesotelin，MSLN) 采用EILSA法在84例胰腺癌患者、40例慢性胰腺炎和84名健康對照中，測定血清IGFBP2和MSLN水平，發(fā)現(xiàn)血清IGFPB2和MSLN診斷胰腺癌的靈敏度分別為22%(P=0.032)和17%(P=0.002)，較CA19-9低。然而，IGFBP2和(或) MSLN兩者聯(lián)合確診了CA19-9漏診的28例胰腺癌中的18例，這提示，血清IGFBP2和MSLN單獨診斷胰腺癌敏感性低，但在聯(lián)合診斷時具參考價值[17]。

1.4 蛋白譜 基于質(zhì)譜的血清肽譜是一種有前途的發(fā)現(xiàn)新的疾病相關(guān)蛋白質(zhì)的工具，Zapico-Muniz等[18]通過基質(zhì)輔助激光解吸/電離-飛行時間質(zhì)譜比較健康對照、慢性胰腺炎和胰腺癌患者的血清肽譜，忽略性別和年齡的影響，聯(lián)合2種血清肽譜和CA19-9診斷胰腺癌靈敏度為89.9%，特異度為92.7%，聯(lián)合3種血清肽譜區(qū)分胰腺癌與健康對照及慢性胰腺炎，其敏感度98.2%，特異度97.1%。Ehmann等[19]通過陰離子交換色譜法 確定載脂蛋白A-Ⅱ，轉(zhuǎn)甲狀腺素蛋白和載脂蛋白AI與胰腺癌相關(guān)，與健康對照組相比，它們在胰腺癌血清中水平均下降至少2倍。同樣經(jīng)此種方法確定血小板第4因子與胰腺癌有關(guān)[20]。采用2-維-凝膠電泳比較胰腺癌及對照血清標本，進行差異蛋白質(zhì)分析顯示，甘露糖結(jié)合凝集素2和肌球蛋白輕鏈激酶2蛋白表達上調(diào)，免疫印跡證實它們在胰腺癌組織表達升高[21]。聯(lián)合蛋白質(zhì)芯片技術(shù)及2維凝膠電泳發(fā)現(xiàn)熱休克蛋白27(heat shock protein 27，HSP27)在胰腺癌與健康對照組表達差異，通過ELISA法檢測35例胰腺癌患者和37名健康個體血清樣本，HSP27診斷胰腺癌靈敏度為100%。特異度為84%[22]。以上研究提示，血清差異蛋白有望作為胰腺癌早期診斷的潛在標記物。

2 K-ras基因突變

已知胰腺癌K-ras基因突變的發(fā)生率很高。Chen等[23]研究發(fā)現(xiàn)在30例(共91例)患者血漿中K-ras12號密碼子突變率為33%(30/91)，其中17個突變?yōu)閏.35G＞A(p.G12D)，11個突變?yōu)閏.35G＞T(p.G12V)，只有2個為c.34G＞C(p. G12R)。K-ras12號密碼子突變與腫瘤TNM分期(P=0.033)和肝轉(zhuǎn)移(P=0.014)有關(guān)。K-ras突變患者的中位生存時間少于野生型K-ras基因患者(3.9個月vs 10.2個月，P＜0.001)，另一項研究在30例胰腺癌患者的血漿和腫瘤組織，40例慢性胰腺炎患者血漿中，采用突變等位基因特異性擴增法和直接測序法檢測12號密碼子突變，同時檢測腫瘤標記物如CA19-9、CEA、CA50、CA242。在70%腫瘤組織樣品中檢測到K-ras突變，但在血漿DNA樣品中未檢測到K-ras突變，60%具有K-ras突變的患者CA199、CEA、CA50和CA242升高，但其靈敏度(20.00% ～ 56.67%)和特異度均偏低(56.67% ～ 77.5%)。這些數(shù)據(jù)表明，在血清中的K-ras基因突變用于篩選胰腺癌患者以及胰腺腫瘤高危進展人群的方法尚需進一步驗證[24]。

3 抑癌基因的甲基化

除基因突變外，抑癌基因的甲基化在腫瘤包括胰腺癌在內(nèi)其發(fā)生發(fā)展中發(fā)揮重要作用，研究它有望找到胰腺癌早期的標記物。Yi等[25]經(jīng)全基因轉(zhuǎn)錄組在胰腺癌細胞系鑒定出新的發(fā)生DNA甲基化的腫瘤特異性基因后，使用納米粒子處理的甲基化磁珠技術(shù)分析患者血清中BNC1和ADAMTS1基因啟動子區(qū)域DNA甲基化，結(jié)果提示143例胰腺癌中BNC1和ADAMTS1基因甲基化率分別為92%和68%，在Ⅲ期胰腺上皮內(nèi)瘤變兩者甲基化率接近100%，在Ⅰ期浸潤性癌其甲基化率為97%，42例胰腺癌患者血清中，BNC1敏感度為79%，ADAMTS1敏感度為48%，BNC1特異度89%，ADAMTS1特異度92%。聯(lián)合兩種標記物的靈敏度提高為81%，特異度僅為85%。Park等[26]研究發(fā)現(xiàn)，胰腺癌患者NPTX2甲基化率與慢性胰腺炎相比顯著升高(P=0.016)，其敏感度和特異度分別為80%和76%，NPTX2基因甲基化水平隨腫瘤進展顯著增高。

除單個基因甲基化的研究外，有報道采用甲基化譜進行胰腺癌的早期診斷。Melnikov等[27]納入30例胰腺癌患者，根據(jù)性別和年齡匹配原則入選30例健康對照組，比較兩者血漿甲基化譜，確定了包括CCND2、SOCS1、THBS1、PLAU、VHL在內(nèi)的基因組，聯(lián)合診斷胰腺癌敏感度為76%，特異度為59%，可應用于胰腺癌早期診斷。

4 MicroRNA

MicroRNA是小的非編碼轉(zhuǎn)錄體，參與包括腫瘤發(fā)生在內(nèi)的多種細胞生理病理過程。最近的研究表明，血漿/血清中穩(wěn)定地檢測到微RNA(miRNA)。與健康對照組相比，miR-200家族的兩個成員中miR-200a和miR-200b在胰腺癌和慢性胰腺炎患者的血清中均顯著升高(P＜0.000 1)，受試者工作曲線下面積分別為0.861和0.85[8]。由于胰腺癌顯示極端低氧特征，Ho等[28]推測miR-210(在低氧下誘導，與腫瘤不良反應有關(guān))可作為胰腺癌篩選或監(jiān)控的診斷性標記物，從血漿直接提取miRNA并逆轉(zhuǎn)錄為cDNA，celmiR-54作為正常對照，采用qRT-PCR檢測miR-210和cel-miR-54，與正常對照相比，胰腺癌患者miR-210表達超過4倍以上(P＜0.001)。有研究報道，胰腺癌患者血漿miR-221濃度顯著高于胰腺良性腫瘤(P=0.016)和健康對照組(P＜0.001)；在術(shù)后樣品中顯著減少(P=0.018)，與遠處轉(zhuǎn)移(P=0.041)和非切除狀態(tài)有關(guān)(P=0.021)[29]。

以上研究提示，胰腺癌相關(guān)血漿miR的測定可作為腫瘤檢測、監(jiān)測腫瘤動力學和評價預后的生物標記物，并可能有助于胰腺癌的臨床治療決策。

5 循環(huán)腫瘤細胞

胰腺癌常常經(jīng)血液轉(zhuǎn)移至遠處器官，如肝、肺和骨骼系統(tǒng)。循環(huán)腫瘤細胞(circulating tumor cells，CTCs)是具有進入循環(huán)系統(tǒng)能力的腫瘤細胞，這些細胞群在腫瘤遠處轉(zhuǎn)移中發(fā)揮舉足輕重的作用。臨床研究顯示，血液中CTCs的存在與乳腺癌、結(jié)直腸癌、前列腺癌的進展有關(guān)，并用來預測轉(zhuǎn)移性腫瘤的生存期[30]。Kurihara等[31]檢測胰腺癌患者的CTCs水平，預測胰腺癌的生存期，26例胰腺癌患者中11例CTC陽性，15例CTC陰性(Ⅱ期1例，Ⅲ期1例，Ⅳa期7例，Ⅳb期6例)，CTC陽性和陰性患者的中位生存時間(median survival time，MSTS)分別為110.5 d和375.8 d (P＜0.001)。對于Ⅳb患者，CTC陽性和陰性患者的MSTS分別為52.5 d和308.3 d (P＜0.01)。這項研究表明檢測外周血的CTCs可能有助于胰腺癌患者的預后判斷。

Han等[32]通過BM7(靶向黏蛋白1，mucin1)和VU1D9抗體(靶向上皮細胞黏附分子，EpCAM)的免疫磁性富集，評價34例全身治療前胰腺癌患者和40例健康對照者CTCs的水平，通過實時定量RT-PCR分析KRT19、MUC1、EPCAM、CEACAM5和BIRC5基因的表達，結(jié)果顯示47.1%胰腺癌CTCs陽性患者，與CTC陰性患者比較，具有較短的中位無進展生存期(progression free survival，PFS)，CTCs陽性患者的中位PFS時間66.0 d，CTC陰性患者的中位PFS 138.0 d (P=0.01)。

一項國際多中心隨機研究納入79例患者，化療開始前及治療2個月后進行CTC篩選，5%患者檢測到一個或多個CTCs，治療2個月后，9%的患者檢測到CTCs(總檢出率為11%)。CTC陽性與較差的腫瘤分化程度有關(guān)(P=0.04)，CTC陽性總生存期更短(RR=2.5，P=0.01)。CTC檢測可作為一種有前途的判斷胰腺癌患者預后的工具[33]。一篇薈萃分析結(jié)果也顯示，CTC陽性與較差無進展生存期顯著相關(guān)(HR=1.89，95% CI=1.25 ～ 4.00，P＜0.001)，與CTC陰性患者相比，CTC陽性胰腺癌患者表現(xiàn)出較差的整體存活率，(HR=1.23，95% CI=0.88 ～ 2.08，P＜0.001)。按種族亞組分析表明，亞洲人和白種人群中CTC陽性患者具有更差的生存期(均P＜0.05)。以上研究提示，在外周血中檢測CTC可能是一種很有前途的用于胰腺癌檢測和預后的生物標記物[34]。

6 結(jié)語

胰腺癌的發(fā)生涉及多個基因、信號途徑、轉(zhuǎn)錄因子、miRNA等相互作用，聯(lián)合多種血清學標志物有助于其早期診斷；Ginesta等[35]在61例細針抽吸胰腺腫塊(43胰腺癌和18例慢性胰腺炎)中，進行HRH2、EN1、SPARC、CDH13和APC基因的甲基化狀態(tài)和K-ras基因突變分析，結(jié)果提示基因甲基化對胰腺癌敏感度73%，特異度100%；K-ras突變敏感度77%，特異度100%；聯(lián)合分析其敏感度升高到84%，特異度100%。已有研究報道采用HMGB1與CA19-9、PAM4和CA 19-9、miR-16和CA19-9聯(lián)合后提高了胰腺癌檢測的靈敏性和特異性[36-38]。新型血清標記物的研究及標記物的聯(lián)合使用提高了胰腺癌診斷診斷和預后評估的敏感性和特異性，已在多個臨床試驗中證實，由于為較小的樣本量，研究結(jié)果可能仍處初步階段，對于胰腺癌早期診斷的臨床試驗仍需進一步完善。

1 Warshaw AL， Fernández-del Castillo C. Pancreatic carcinoma［J］. N Engl J Med， 1992， 326（7）：455-465.

2 Hu J， Zhao G， Wang HX， et al. A meta-analysis of gemcitabine containing chemotherapy for locally advanced and metastatic pancreatic adenocarcinoma［J］. J Hematol Oncol， 2011， 4： 11.

3 Carpelan-Holmstr?m M， Nordling S， Pukkala E， et al. Does anyone survive pancreatic ductal adenocarcinoma? A nationwide study reevaluating the data of the Finnish Cancer Registry［J］. Gut， 2005，54（3）：385-387.

4 Benson AB 3rd. Adjuvant therapy for pancreatic cancer： one small step forward［J］. JAMA， 2007， 297（3）：311-313.

5 Herreros-Villanueva M， Gironella M， Castells A， et al. Molecular markers in pancreatic cancer diagnosis［J］. Clinica Chimica Acta，2013， 418： 22-29.

6 Harsha HC， Kandasamy K， Ranganathan PA， et al. A compendium of potential biomarkers of pancreatic cancer［J］. PLoS Med， 2009，6（4）： e1000046.

7 Van Kampen JG， Marijnissen-van Zanten MA2， Simmer F3， et al. Epigenetic targeting in pancreatic cancer［J］. Cancer Treat Rev，2014， 40（5）：656-664.

8 Li A， Omura N， Hong SM， et al. Pancreatic cancers epigenetically silence SIP1 and hypomethylate and overexpress miR-200a/200b in association with elevated circulating miR-200a and miR-200b levels［J］. Cancer Res， 2010， 70（13）： 5226-5237.

9 Cen PT， Ni XL， Yang JX， et al. Circulating tumor cells in the diagnosis and management of pancreatic cancer［J］. Biochim Biophys Acta， 2012， 1826（2）： 350-356.

10 Gold DV， Lew K， Maliniak R， Hernandez M， et al. Characterization of monoclonal antibody PAM4 reactive with a pancreatic cancer mucin［J］. Int J Cancer， 1994， 57（2）：204-210.

11 Gold DV， Goggins M， Modrak DE， et al. Detection of Early-Stage pancreatic adenocarcinoma［J］. Cancer Epidemiol Biomarkers Prev， 2010， 19（11）： 2786-2794.

12 Gold DV， Karanjawala Z， Modrak DE， et al. PAM4-reactive MUC1 is a biomarker for early pancreatic adenocarcinoma［J］. Clin Cancer Res， 2007， 13（24）： 7380-7387.

13 Gold DV， Gaedcke J， Ghadimi BM， et al. PAM4 enzyme immunoassay alone and in combination with CA 19-9 for the detection of pancreatic adenocarcinoma［J］. Cancer， 2013， 119（3）： 522-528.

14 Palumbo R， Sampaolesi M， De Marchis F， et al. Extracellular HMGB1， a signal of tissue damage， induces mesoangioblast migration and proliferation［J］. J Cell Biol， 2004， 164（3）： 441-449.

15 Nehil M， Paquette J， Tokuyasu T， et al. High mobility group box 1 promotes tumor cell migration through epigenetic silencing of semaphorin 3A［J/OL］. http：//www.nature.com/onc/journal/vaop/ ncurrent/full/onc2013459a.html

16 Chung HW， Lim JB， Jang S， et al. Serum high mobility group box-1 is a powerful diagnostic and prognostic biomarker for pancreatic ductal adenocarcinoma［J］. Cancer Sci， 2012， 103（9）： 1714-1721.

17 Kendrick ZW， Firpo MA， Repko RC， et al. Serum IGFBP2 and MSLN as diagnostic and prognostic biomarkers for pancreatic cancer［J］. HPB （Oxford）， 2014， 16（7）：670-676.

18 Zapico-Muniz E， Farre-Viladrich A， Rico-Santana N， et al. Standardized peptidome profiling of human serum for the detection of pancreatic cancer［J］. Pancreas， 2010， 39（8）： 1293-1298.

19 Ehmann M， Felix K， Hartmann D， et al. Identification of potential markers for the detection of pancreatic cancer through comparative serum protein expression profiling［J］. Pancreas， 2007， 34（2）：205-214.

20 Fiedler GM， Leichtle AB， Kase J， et al. Serum peptidome profiling revealed platelet factor 4 as a potential discriminating peptide associated with pancreatic cancer［J］. Clin Cancer Res， 2009， 15（11）： 3812-3819.

21 Rong YF， Jin DY， Hou CR， et al. Proteomics analysis of serum protein profiling in pancreatic cancer patients by DIGE： upregulation of mannose-binding lectin 2 and myosin light chain kinase 2［J］. BMC Gastroenterol， 2010， 10： 68.

22 Melle C， Ernst G， Escher N， et al. Protein profiling of microdissected pancreas carcinoma and identification of HSP27 as a potential serum marker［J］. Clin Chem， 2007， 53（4）： 629-635.

23 Chen H， Tu H， Meng ZQ， et al. K-ras mutational status predicts poor prognosis in unresectable pancreatic cancer［J］. Eur J Surg Oncol，2010， 36（7）： 657-662.

24 Marchese R， Muleti A， Pasqualetti PA， et al. Low correspondence between K-ras mutations in pancreatic cancer tissue and detection of K-ras mutations in circulating DNA［J］. Pancreas， 2006， 32（2）：171-177.

25 Yi JM， Guzzetta AA， Bailey VJ， et al. Novel methylation biomarker panel for the early detection of pancreatic cancer［J］. Clin Cancer Res， 2013， 19（23）： 6544-6555.

26 Park JK， Ryu JK， Yoon WJ， et al. The role of quantitative NPTX2 hypermethylation as a novel serum diagnostic marker in pancreatic cancer［J］. Pancreas， 2012， 41（1）： 95-101.

27 Melnikov AA， Scholtens D， Talamonti MS， et al. Methylation profile of circulating plasma DNA in patients with pancreatic cancer［J］. J Surg Oncol， 2009， 99（2）： 119-122.

28 Ho AS， Huang X， Cao HB， et al. Circulating miR-210 as a novel hypoxia marker in pancreatic cancer［J］. Transl Oncol， 2010， 3（2）：109-113.

29 Kawaguchi T， Komatsu S， Ichikawa D， et al. Clinical impact of circulating miR-221 in plasma of patients with pancreatic cancer［J］. Br J Cancer， 2013， 108（2）： 361-369.

30 Tjensvoll K， Nordgard O， Smaaland R. Circulating tumor cells in pancreatic cancer patients： Methods of detection and clinical implications［J］. Int J Cancer， 2014， 134（1）： 1-8.

31 Kurihara T， Itoi T， Sofuni A， et al. Detection of circulating tumor cells in patients with pancreatic cancer： a preliminary result［J］. J Hepatobiliary Pancreat Surg， 2008， 15（2）： 189-195.

32 Han L， Chen W， Zhao QC. Prognostic value of circulating tumor cells in patients with pancreatic cancer： a meta-analysis［J］. Tumour Biol， 2014， 35（3）： 2473-2480.

33 Bidard FC， Huguet F， Louvet C， et al. Circulating tumor cells in locally advanced pancreatic adenocarcinoma： the ancillary Circe 07 study to the LAP 07 trial［J］. Ann Oncol， 2013， 24（8）： 2057-2061.

34 De Albuquerque A， Kubisch I， Breier G， et al. Multimarker gene analysis of circulating tumor cells in pancreatic cancer patients： a feasibility study［J］. Oncology， 2012， 82（1）： 3-10.

35 Ginesta MM， Mora J， Mayor R， et al. Genetic and epigenetic markers in the evaluation of pancreatic masses［J］. J Clin Pathol， 2013， 66（3）： 192-197.

36 王玉彬，孫麗麗，王曉雪．miRNAs聯(lián)合CA19-9對胰腺癌的診斷效能［J］.軍醫(yī)進修學院學報，2012（7）：772-774．

37 Chung HW， Lee SG， Kim H， et al. Serum high mobility group box-1（HMGB1） is closely associated with the clinical and pathologic features of gastric cancer［J］. J Transl Med， 2009， 7： 38-48.

38 Cheng BQ， Jia CQ， Liu CT， et al. Serum high mobility group box chromosomal protein 1 is associated with clinicopathologic features in patients with hepatocellular carcinoma［J］. Dig Liver Dis， 2008，40（6）： 446-452.

Serological markers of pancreatic cancer in early diagnosis of pancreatic cancer

ZENG Qiang, GONG Yan
Health Management and Research Institute, Chinese PLA General Hospital, Beijing 100853, China

Pancreatic cancer is a kind of malignant tumor which is very diff i cult to diagnose and cure. Researches in developing convenient serological markers in screening the early curable stage of pancreatic cancer have important clinical signif i cance. At present, the study of serological marker of pancreatic cancer includes the molecular of pancreatic cancer-specif i c expression, gene mutations, abnormal methylation of gene DNA, microRNA and circulating tumor cells. Combinations of different serological markers can contribute to the early diagnosis of pancreatic cancer, and further researches should focus on the screening of high-risk groups with pancreatic cancer and the development of highly sensitive and specif i c markers, in order to improve the diagnostic level of pancreatic cancer.

pancreatic cancer; serum markers; early diagnosis

R 730.3

A

2095-5227(2015)02-0184-04

10.3969/j.issn.2095-5227.2015.02.025

時間：2014-10-11 10:23

http://www.cnki.net/kcms/detail/11.3275.R.20141011.1023.002.html

2014-08-07

全軍十一五計劃保健專項(10BJZ18)

Supported by the 11th Five Years Programs for Health Care Plan of Chinese PLA (10BJZ18)

曾強，男，博士，研究員，主任，教授。研究方向：健康管理。Email: zq301@126.com

The fi rst author: ZENG Qiang. Email: zq301@126.com