許燎原，趙麗穩(wěn)，胡宇峰2，邱海萍3，柴榮耀3，張 震（.寧波市種植業(yè)管理總站，浙江寧波3502；2.浙江省寧?？h農(nóng)業(yè)技術(shù)推廣總站，浙江寧海35600；3.浙江省農(nóng)業(yè)科學(xué)院植物保護(hù)與微生物研究所，浙江杭州3002）

稻瘟病菌孢子qPCR方法的建立及用于監(jiān)測氣傳菌源的研究

許燎原1，趙麗穩(wěn)1，胡宇峰2，邱海萍3，柴榮耀3，張 震3，*
（1.寧波市種植業(yè)管理總站，浙江寧波315012；2.浙江省寧?？h農(nóng)業(yè)技術(shù)推廣總站，浙江寧海315600；3.浙江省農(nóng)業(yè)科學(xué)院植物保護(hù)與微生物研究所，浙江杭州310021）

稻瘟病是水稻生產(chǎn)中最具毀滅性的真菌病害之一，監(jiān)測田間空氣中稻瘟病菌分生孢子數(shù)量動態(tài)變化，對于預(yù)測病情發(fā)生和指導(dǎo)及時防治具有重要意義。該研究以稻瘟病菌MHP1基因為靶序列，建立了1種稻瘟病菌分生孢子qPCR定量方法，并采用此方法對田間氣傳稻瘟病菌源進(jìn)行了監(jiān)測。監(jiān)測數(shù)據(jù)顯示，監(jiān)測期間空氣中稻瘟病菌分生孢子數(shù)量波動幅度較大。結(jié)合氣象條件進(jìn)行分析，菌源量的波動與氣象條件相關(guān)。該研究建立的孢子qPCR方法可應(yīng)用于監(jiān)測氣傳稻瘟病菌源。

稻瘟?。环稚咦恿?；孢子捕捉；qPCR；氣象因子

由子囊菌Magnaporthe oryzae侵染引起的稻瘟病是水稻生產(chǎn)中最具毀滅性的真菌病害之一［1］，該病屬于單年多循環(huán)病害，病稻草上的越冬菌源和當(dāng)年病斑上產(chǎn)生的分生孢子經(jīng)由氣流傳播是稻瘟病發(fā)生的初侵染源和再侵染源。因此，及時監(jiān)控稻田空氣中的稻瘟病菌源量，可以為有效進(jìn)行病害的控制和管理提供依據(jù)。

長期以來，空氣中稻瘟病菌源監(jiān)測以凡士林玻片或旋轉(zhuǎn)式孢子捕捉器等捕捉樣品的顯微計數(shù)為主［2-3］。顯微計數(shù)不僅工作量大，而且鑒于環(huán)境樣品的復(fù)雜性，包含固體顆粒物和形態(tài)近似孢子等，顯微計數(shù)的結(jié)果可能與實際情況存在較大偏差?；赥aqman的qPCR技術(shù)，以其特異性強、靈敏度高、重復(fù)性好、速度快等優(yōu)點已廣泛應(yīng)用于植物病原菌的定量監(jiān)測［4］。這些菌源定量監(jiān)測方法通常以病原真菌的rDNA區(qū)域保守序列為靶標(biāo)設(shè)計特異性的Taqman探針［5-6］。但由于這些序列與病原菌的致病性無關(guān)，其用于監(jiān)測環(huán)境樣本中致病菌源量，所得數(shù)據(jù)中包含大量的非致病菌，不利于病害的準(zhǔn)確預(yù)測。如將與病菌致病性相關(guān)的特異性標(biāo)記結(jié)合應(yīng)用于病害的流行預(yù)警，所得數(shù)據(jù)更具可靠性。近來，Su'udi等［7］根據(jù)稻瘟病菌致病相關(guān)基因MHP1設(shè)計了用于稻瘟病菌菌量檢測的Taqman探針，并成功實現(xiàn)了對病害嚴(yán)重度的比較?；谒麄兊难芯拷Y(jié)果，本研究研發(fā)空氣傳播的稻瘟病分生孢子定量檢測方法，并開展水稻生育期內(nèi)稻瘟病菌空氣傳播分生孢子數(shù)量的動態(tài)監(jiān)測，對及時預(yù)測田間病害發(fā)生情況具有重要的現(xiàn)實意義。

1 材料與方法

1.1 菌株

稻瘟病菌菌株12-10-1和小麥赤霉病菌12-02均由浙江省農(nóng)業(yè)科學(xué)院植物保護(hù)與微生物研究所糧食作物病害研究室分離保存。

稻瘟病菌分生孢子獲得：濾紙片上保存的菌株于CM平板28℃活化后，切取適量菌塊接種到CM平板，28℃、12 h（L）/12 h（D）光暗交替培養(yǎng)8 d。平皿內(nèi)倒入適量無菌水，用高壓滅菌的毛筆輕刷菌落，吸取菌懸液經(jīng)4層滅菌擦鏡紙過濾，濾液經(jīng)2 500 g離心10 min收集分生孢子，加入適量TE溶液重懸分生孢子，用血球計數(shù)板計數(shù)后調(diào)整孢子濃度至1.5×107個·mL-1，-20℃保存待用。

小麥赤霉病菌菌株于PDA平板28℃活化后，切取適量菌塊于100 mL 4%綠豆湯培養(yǎng)液中，25℃150 r·min-1振蕩培養(yǎng)5 d。吸取菌液經(jīng)4層滅菌擦鏡紙過濾，收集濾液經(jīng)2 500 g離心收集分生孢子，加入適量TE溶液懸浮，用血球計數(shù)板計數(shù)后調(diào)整孢子濃度至1.5×107個·mL-1，-20℃保存待用。

1.2 主要儀器與試劑

7500型實時熒光定量PCR儀和Premix Ex TaqTM熒光定量試劑盒，用于核酸的定量檢測，分別為美國Applied Biosystems公司和TaKaRa寶生物工程（大連）有限公司產(chǎn)品。Fastprep-24和酸洗玻璃珠（425—600 μm）用于分生孢子的機(jī)械破壁，分別購自美國MP Biomedicals公司和Sigma-Aldrich公司。7 d孢子捕捉儀，用于田間空氣傳播的稻瘟病菌分生孢子捕捉，是英國Burkard Scientific公司產(chǎn)品，屬吸入式7 d連續(xù)孢子捕捉儀，進(jìn)氣量10 L·min-1。HOBO u30小型自動氣象站，用于田間小氣候自動監(jiān)測與記錄，為美國Onset公司的產(chǎn)品。

1.3 稻瘟病菌分生孢子基因組DNA的提取

用TE溶液對上述稻瘟病菌孢子液進(jìn)行10倍梯度稀釋，配制1.5×107、1.5×106、1.5×105、1.5×104、1.5×103和1.5×102個·mL-1稻瘟病菌孢子懸浮液。對4個懸浮液中稻瘟病菌孢子濃度不足1.5×107個·mL-1，用赤霉病菌孢子母液進(jìn)行補充，使得各懸浮液中真菌孢子量保持一致。取上述不同濃度的稻瘟病菌孢子懸浮液400 μL用于提取基因組DNA。DNA提取相關(guān)試劑及過程參見Ward［8］的方法。

1.4 引物及熒光探針設(shè)計

引物和 Taqman探針參照 Su'udi等［7］的序列。引物序列為（MHP1F1）5′-TCGATGCCGACAACTTCTCCG-3′和（MHP1R1）5′-ACCCTGGTCAAGCTGTTCGATTGT-3′；探針序列：5′-TGCCCATCGTAAGTTCCTTCTTCGCA-3′，其中5′和3′分別采用FAM和TAMRA進(jìn)行修飾。上述引物探針均由美國Invitrogen公司合成。

1.5 qPCR擴(kuò)增體系

qPCR擴(kuò)增體系和反應(yīng)程序參照TaKaRa寶生物工程（大連）有限公司Premix Ex TaqTM熒光定量試劑盒操作說明，略有改動。反應(yīng)體系為2× Premix Ex TaqTM10 μL、10 μmol·L-1上下游引物及探針溶液各0.4 μL、50×ROX Reference DyeⅡ0.4 μL、DNA模板2 μL、ddH2O 6.4 μL，總體積為20 μL。反應(yīng)程序采用2步法。qPCR結(jié)果利用7500熒光定量PCR儀自帶軟件SDS（版本v2.0.4）進(jìn)行分析。

1.6 田間空氣中稻瘟病菌孢子的捕捉

孢子捕捉點位于浙江省寧?？h茶院鄉(xiāng)稻作區(qū)（121.42°E，29.3°N）。周邊種植的水稻品種以甬優(yōu)系列秈粳雜交水稻為主。7 d孢子捕捉儀進(jìn)氣口距地面1.2 m。捕捉儀整年開啟，每隔7 d更換用于收集樣品的聚酯片。取下的聚酯片進(jìn)行7等分，每份即為1 d的孢子捕捉樣品。將每份聚酯片剪成2 mm寬的窄條裝于2 mL離心管內(nèi)，并用Ward［8］的方法提取基因組DNA，最終獲得50 μL的DNA溶液。

1.7 田間小氣候條件的監(jiān)測

HOBO u30小型氣象站安裝于7 d孢子捕捉儀旁，整年開啟，實時監(jiān)測并記錄取樣點溫度、降雨量、相對濕度和風(fēng)速等氣象數(shù)據(jù)。

圖1 稻瘟病菌分生孢子定量標(biāo)準(zhǔn)曲線Fig.1 Standard curve used for quantifying rice blast spore

2 結(jié)果與分析

2.1 稻瘟病菌分生孢子定量標(biāo)準(zhǔn)曲線的建立

取上述1.5×107、1.5×106、1.5×105、1.5× 104、1.5×103和1.5×102個·mL-1的稻瘟病菌分生孢子懸浮液400 μL，分別相當(dāng)于取6×106、6× 105、6×104、6×103、6×102和60個稻瘟病菌分生孢子，用于提取基因組DNA，最終各樣品均獲得50 μL的DNA溶液。分別取2 μL作為模板進(jìn)行qPCR分析，結(jié)果顯示，除孢子濃度為1.5×102個·mL-1的樣品外，其他5個樣品擴(kuò)增曲線明顯，平均Ct值分別為21.68、25.82、30.27、33.09和35.14。據(jù)Su'udi等［7］報道，檢測下限為1 pg稻瘟病菌基因組DNA，按1 pg基因組DNA=0.978× 109bp計算［9］，稻瘟病菌基因組大小約40 Mbp，分生孢子具有3個細(xì)胞，其理論最低可檢出8個稻瘟病菌分生孢子。本研究中，孢子濃度為1.5× 102個·mL-1的樣品沒有擴(kuò)增曲線，這與其孢子數(shù)量低于Su'udi等［7］報道的最低檢測下限相符合。對5個有擴(kuò)增結(jié)果的樣品，按孢子數(shù)量的lg值與平均Ct值進(jìn)行線性回歸，結(jié)果顯示具有良好的相關(guān)性，R2=0.977 9（圖1）。

2.2 空氣中孢子捕捉樣品的檢測

對2012年6月7日—10月28日采集的21批次空氣中稻瘟病菌分生孢子捕捉樣品進(jìn)行了檢測。收集的聚酯片按天進(jìn)行等距分割，并提取基因組DNA，每天樣品獲得50 μL的DNA溶液。每個樣品取2 μL，進(jìn)行qPCR檢測，獲得各樣品的Ct值。根據(jù)上述稻瘟病菌分生孢子定量的標(biāo)準(zhǔn)曲線，計算各PCR樣品所對應(yīng)的稻瘟病菌分生孢子數(shù)量，然后根據(jù)7 d捕捉儀每天吸氣14 400 L，計算每立方空氣中所含的稻瘟病菌分生孢子數(shù)量。結(jié)果顯示，2012年6月7日—10月28日，稻田空氣中稻瘟病菌分生孢子數(shù)量呈現(xiàn)明顯的動態(tài)變化（圖2），在孢子監(jiān)測的144 d中有64 d未檢測到稻瘟病菌分生孢子，其中有53 d集中在8月25日至10月28日；而稻瘟病菌分生孢子釋放主要集中在6月7日至8月24日，檢出的分生孢子數(shù)量為16.9～825.3 cfu·m-3空氣，其中，8 月9日檢出孢子數(shù)量最多。

2.3 空氣中稻瘟病菌孢子數(shù)量與氣象因子的關(guān)系

生物學(xué)研究表明，稻瘟病菌菌絲可以在8～37℃生長，適宜溫度為26～28℃；分生孢子在10～35℃均可形成，最適溫度為25～28℃。根據(jù)HOBO u30氣象站數(shù)據(jù)，2012年6月7日至10 月28日，日平均氣溫為14.4～30.9℃，日最高氣溫為37.4℃，日最低溫度除了10月18日、19日、23日和24日外均在10℃以上（圖2），此溫度條件下稻瘟病菌均能夠生長和產(chǎn)孢。自然條件下，稻瘟病菌分生孢子大多數(shù)在夜間產(chǎn)生和釋放，特別是在2：00—6：00［8］。定量監(jiān)測結(jié)果顯示，稻瘟病菌分生孢子釋放主要集中在6月7日至8月24日，8月25日至10月28日監(jiān)測到的孢子數(shù)量較少，僅有12 d可監(jiān)測到稻瘟病菌孢子的存在。對每天夜間2：00—6：00的平均氣溫進(jìn)行分析發(fā)現(xiàn)，6月7日至8月24日的氣溫為24～28℃，是稻瘟病菌的最適產(chǎn)孢溫度；而8月25日至10月28日的絕大部分時間氣溫都低于20℃，不利于稻瘟病菌分生孢子的產(chǎn)生和釋放。

圖2 稻瘟病菌分生孢子數(shù)量與氣溫的動態(tài)變化Fig.2 Dynamic changes of temperature and number of rice blast spore

圖3 稻瘟病菌分生孢子數(shù)量與空氣相對濕度的動態(tài)變化Fig.3 Dynamic changes of relative air humidity and number of rice blast spore

稻瘟病菌分生孢子產(chǎn)生需要較高的相對濕度。本研究中，6月7日至10月28日中有79 d日平均相對濕度在60.5%以上，且絕大部分時間在85%以上（圖3）；2：00—6：00的空氣相對濕度絕大部分在90%以上，葉面濕度在100%?？傮w上，監(jiān)測點空氣相對濕度有利于稻瘟病病菌分生孢子的產(chǎn)生。而實際監(jiān)測數(shù)據(jù)顯示，9月和10月孢子數(shù)量較少，說明空氣相對濕度符合的條件下，當(dāng)年監(jiān)測點稻瘟病菌孢子產(chǎn)生主要受氣溫的影響，這與稻瘟病菌生物學(xué)特性的研究結(jié)果相符合［10-11］。

降雨可以增加空氣相對濕度，有利于滿足稻瘟病菌孢子形成所需的高濕環(huán)境條件，但雨水的沖刷不利于孢子向空中傳播。本研究的數(shù)據(jù)顯示，降雨對稻瘟病菌孢子經(jīng)由空氣傳播的影響與此相符（圖4）。6月7日至8月24日孢子集中釋放期間，降雨后往往伴隨著孢子的大量釋放，典型的如7月16日（孢子數(shù)量達(dá)到507.9個·m-3空氣）、8月5日（孢子數(shù)量達(dá)到523.8個·m-3空氣）和8月9日（孢子數(shù)量達(dá)到852.3個·m-3空氣），這3個日期之前均有連續(xù)2～3 d的降雨，而降雨當(dāng)天孢子檢出量相對較少，甚至未有孢子檢出。

風(fēng)力不僅影響稻瘟病菌分生孢子的機(jī)械脫落，也影響孢子在空氣中的傳播。本研究中監(jiān)測點鄰近浙江省三門灣，受海洋氣候的影響，取樣時間內(nèi)絕大部分日期最大瞬時風(fēng)速在5 m·s-1以上（圖5），其有利于引起分生孢子的機(jī)械脫落和傳播［12］。8月5日和8月9日孢子大量釋放的時間，都伴隨有7 m·s-1以上的最大瞬時風(fēng)速。

圖4 稻瘟病菌分生孢子數(shù)量與降雨量的動態(tài)變化Fig.4 Dynamic changes of rainfall and number of rice blast spore

圖5 稻瘟病菌分生孢子數(shù)量與最大風(fēng)速的動態(tài)變化Fig.5 Dynamic changes of wind velocity and number of rice blast spore

3 結(jié)論與討論

監(jiān)測空氣中稻瘟病菌分生孢子數(shù)量的動態(tài)變化，對于稻瘟病發(fā)生為害的預(yù)測預(yù)報及指導(dǎo)防治具有重要的現(xiàn)實意義。本研究利用 Su'udi等［7］設(shè)計的稻瘟病菌特異性Taqman探針，成功建立了1種稻瘟病菌分生孢子的定量方法，并應(yīng)用其對稻田空氣中稻瘟病菌分生孢子數(shù)量進(jìn)行了監(jiān)測。Su'udi等［7］的方法可檢測出至少8個分生孢子的基因組DNA，本研究建立稻瘟病菌分生孢子定量標(biāo)準(zhǔn)曲線時，發(fā)現(xiàn)1.5×102個·mL-1的樣品未有擴(kuò)增曲線，最低檢測限為2.4～24個分生孢子基因組DNA?？傮w上，本研究建立的稻瘟病菌分生孢子定量方法的靈敏度已足夠用于田間菌源量的動態(tài)監(jiān)測。

對孢子田間捕捉樣品的定量分析發(fā)現(xiàn)，水稻生育期空氣中稻瘟病菌分生孢子數(shù)量波動幅度較大，稻瘟病菌分生孢子釋放集中在6月7日至8月24日，基本與這段時間溫、濕度條件比較適宜稻瘟病菌分生孢子產(chǎn)生和釋放有關(guān)［10］。寧海絕大部分時間風(fēng)速較大，有利于孢子分散。據(jù)當(dāng)?shù)剞r(nóng)業(yè)技術(shù)部門調(diào)查，附近鄉(xiāng)鎮(zhèn)當(dāng)年種植的部分品種有葉瘟中度發(fā)生（個別品種重度發(fā)生）。早期的研究證明，稻瘟病菌株能夠侵染牛筋草、狗尾草、馬唐、稗、李氏禾、羅氏草和鋪地黍等禾本科雜草［13-14］，這些雜草在該稻區(qū)分布廣泛。采樣點鄰近水庫和溪流，多山地，6月7日至8月24日有利的溫、濕度條件使得稻瘟病菌在這些雜草上迅速繁殖。

田間稻瘟病的發(fā)生程度受菌源量、寄主抗性水平和氣象條件的共同決定。在明確品種抗性水平的前提下，利用本研究方法開展空氣中稻瘟病菌孢子數(shù)量的監(jiān)測，結(jié)合氣象條件，及時預(yù)測預(yù)報是否有利于病害發(fā)生、發(fā)展和流行，能夠?qū)Σ『皶r進(jìn)行防控。

（References）：

［1］ SKAMNIOTI P，GURR S J.Against the grain：safeguarding rice from rice blast disease［J］.Trends in Biotechnology，2009，27（3）：141-150.

［2］ 陳浩，康曉慧，張梅，等.稻瘟病菌空中孢子量與氣象因子的關(guān)系［J］.湖北農(nóng)業(yè)科學(xué)，2009，48（8）：1869-1871. CHEN H，KANG X H，ZHANG M，et al.Relationship between spatial spore amount of Pyricularia oryzae Cavara and meteorological factors［J］.Hubei Agricultural Sciences，2009，48（8）：1869-1871.（in Chinese with English abstract）

［3］ 宋成艷，王桂玲，劉乃生，等.寒地稻瘟病菌空中捕捉量與田間稻瘟病發(fā)病和氣象因子的關(guān)系［J］.黑龍江農(nóng)業(yè)科學(xué)，2011（9）：39-41. SONG C Y，WANG G L，LIU N S，et al.Relationship of the amounts of germs of rice blast caught in the air with the incidence of rice blast in the field and the meteorological factors in cold region［J］.Heilongjiang Agricultural Sciences，2011 （9）：39-41.（in Chinese with English abstract）

［4］ WEST J S，ATKINS S D，EMBERLIN J，et al.PCR to predict risk of airborne disease［J］.Cell，2008，16（8）：380 -387.

［5］ BORNEMAN J，HARTIN R J.PCR primers that amplify fungal rRNA genes from environmental samples［J］.Applied and Environmental Microbiology，2000，66（10）：4356-4360.

［6］ MARTIN K J，RYGIEWICZ P T.Fungal-specific PCR primers developed for analysis of the ITS region of environmental DNA extracts［J］.BMC Microbiology，2005，5（1）：1-11.

［7］ SU'UDI M，KIM J，PARK J M，et al.Quantification of rice blast disease progressions through Taqman real-time PCR［J］. Molecular Biotechnology，2013，55（1）：43-48.

［8］ WARD E.Use of molecular methods for the detection of fungal spores［M］.BURNS R.Plant pathology.New York：Humana Press，2009：147-159.

［9］ DOLEZEL J，BARTOS J，VOGLMAYR H，et al.Nuclear DNA content and genome size of trout and human［J］.Cytometry.Part A：the Journal of the International Society for Analytical Cytology，2003，51（2）：127-128.

［10］ OU S H.Rice diseases［M］.2nd ed，Surrey：Commonwealth Mycological Institute，1989：109-200.

［11］ 唐力瓊，張玲，胡運高，等.稻瘟病菌生物學(xué)特性及產(chǎn)孢條件研究［J］.河南農(nóng)業(yè)科學(xué)，2012，41（6）：97-100. TANG L Q，ZHANG L，HU Y G，et al.Biological characteristics and sporulation condition of Magnaporthe grisea［J］. Journal of Henan Agricultural Sciences，2012，41（6）：97-100.（in Chinese with English abstract）

［12］ 曹青，房輝，何有才，等.稻瘟病菌孢子傳播的影響因素與田間捕捉方法［J］.信陽農(nóng)業(yè)高等?？茖W(xué)校學(xué)報，2004，14（3）：7-9. CAO Q，F(xiàn)ANG H，HE Y C，et al.The influencing factors of Pyricularia oryzae Cavara spore distribution and spore trapping methods in fields［J］.Journal of Xingyang Agricultural College，2004，14（3）：7-9.（in Chinese with English abstract）

［13］ 杜新法，孫漱沅，陶榮祥，等.稻區(qū)雜草上的梨孢菌與稻瘟病發(fā)生的關(guān)系［J］.植物病理學(xué)報，1997，27（4）：33-37. DU X F，SUN S Y，TAO R X，et al.Effects of weed-hosts of Pyricularia on incidence of rice blast disease［J］.Acta Phytopathologica Sinica，1997，27（4）：33-37.（in Chinese with English abstract）

［14］ 沈瑛，袁筱萍，王艷麗.分子探針在稻瘟病流行病學(xué)中的應(yīng)用研究［J］.西南農(nóng)業(yè)大學(xué)學(xué)報，1998，20（5）：401 -408. SHEN Y，YUAN X P，WANG Y L.Molecular probes for epidemiological studies of the rice blast fungus［J］.Journal of Southwest Agricultural University，1998，20（5）：401-408. （in Chinese with English abstract）

（責(zé)任編輯 侯春曉）

Development of a qPCR detection method for monitoring conidial density of rice blast fungus in the air

XU Liao-yuan1，ZHAO Li-wen1，HU Yu-feng2，QIU Hai-ping3，CHAI Rong-yao3，ZHANG Zhen3，*
（1.Ningbo Plant Cultivation Management Station，Ningbo 315012，China；2.Ninghai Agricultural Technology Extension and Service Center，Ninghai 315600，China；3.Institute of Plant Protection and Microbiology，Zhejiang A-cademy of Agricultural Sciences，Hangzhou 310021，China）

Rice blast disease，caused by Magnaporthe oryzae，was one of the most devastating fungal diseases of rice.It was very important to monitor airborne conidial density of M.oryzae in fields for forecasting and control of blast disease.In present study，a qPCR detection method for conidia of M.oryzae was developed，based on M. oryzae MHP1 gene as the target sequence.Airborne conidia of M.oryzae in fields were trapped and quantified using this method.The results showed that airborne conidia population of M.oryzae was largely fluctuant，and indicated that the conidia release of M.oryzae was obviously affected by meteorological conditions during the monitoring periods.Generally，this qPCR method could be used to detect and quantified airborne conidia of rice blast fungus in fields.

rice blast；amount of conidia；spore trap；qPCR；meteorological factor

S435.121.4+5

A

1004-1524（2016）08-1368-06

10.3969/j.issn.1004-1524.2016.08.14

2016-04-11

國家公益性行業(yè)（農(nóng)業(yè)）科研專項（201203014）；浙江省公益技術(shù)研究農(nóng)業(yè)項目（2014C32018）

許燎原（1972—），男，浙江寧波人，高級農(nóng)藝師，主要從事植物保護(hù)等研究。E-mail：ddd0574@126.com
*

，張震，E-mail：zzhangcn@126.com