蛇龍珠等釀酒葡萄品種羧酸酯酶基因片段的SNP分析

陳小宇， 劉林德*， 葛宜和， 張 莉，李記明， 姜文廣， 裴廣仁

(1.魯東大學(xué)生命科學(xué)學(xué)院，山東煙臺(tái) 264025;2.煙臺(tái)張?jiān)＜瘓F(tuán)有限公司，山東省葡萄酒微生物發(fā)酵技術(shù)重點(diǎn)實(shí)驗(yàn)室，山東煙臺(tái) 264001)

?

蛇龍珠等釀酒葡萄品種羧酸酯酶基因片段的SNP分析

陳小宇1， 劉林德1*， 葛宜和1， 張 莉1，李記明2*， 姜文廣2， 裴廣仁2

(1.魯東大學(xué)生命科學(xué)學(xué)院，山東煙臺(tái) 264025;2.煙臺(tái)張?jiān)＜瘓F(tuán)有限公司，山東省葡萄酒微生物發(fā)酵技術(shù)重點(diǎn)實(shí)驗(yàn)室，山東煙臺(tái) 264001)

[目的]利用單核苷酸多態(tài)笥(SNP)分析蛇龍珠等釀酒葡萄的單核苷酸多態(tài)性及遺傳親緣關(guān)系，為今后研究釀酒葡萄的抗性機(jī)制以及培育抗病性葡萄提供參考。[方法]以釀酒葡蛇龍珠8個(gè)新株系(E01～E08)、品麗珠、 赤霞珠為材料，以歐亞種黑比諾為對(duì)照，對(duì)羧酸酯酶基因(CXEs)片段進(jìn)行克隆和序列分析，研究蛇龍珠等釀酒葡萄之間的遺傳差異。[結(jié)果]序列對(duì)比顯示，該基因片段在編碼區(qū)和非編碼區(qū)均存在不同程度的堿基替換。11份葡萄材料的基因片段共發(fā)現(xiàn)87個(gè)多態(tài)性位點(diǎn)，平均74 bp左右的序列長(zhǎng)度上能夠檢測(cè)到一個(gè)SNP位點(diǎn)，核苷酸多樣度(Pi)0.050 76±0.011 49，單倍型多樣度(Hd)1.000±0.039，表現(xiàn)出豐富的遺傳多樣性。3種中性檢驗(yàn)方法比較序列變異模式，結(jié)果符合中性生化模型。通過聚類分析可將蛇龍珠E02、E03聚為一類，E07、E08聚為一類，其遺傳距離和地理距離一致。該序列片段存在一些突變，這些突變能將蛇龍珠8個(gè)新株系與其他釀酒葡萄品種明顯區(qū)分。[結(jié)論]CXEs基因片段的突變位點(diǎn)可作為遺傳標(biāo)記，利用SNP方法可分析葡萄的遺傳多樣性和親緣關(guān)系。

蛇龍珠；羧酸酯酶；基因克??；序列分析

羧酸酯酶(Carboxylesterases，CXEs)是一種多功能酶，廣泛存在于生物界中，能夠催化各種水解反應(yīng)[1]。羧酸酯酶的功能較多，在哺乳動(dòng)物和昆蟲體內(nèi)分別參與藥物的代謝和殺蟲劑的代謝，在微生物中參與降解作用[2-3]。目前關(guān)于植物羧酸酯酶的研究較少。植物中的羧酸酯酶在代謝殺蟲劑[4]和病原菌的致病過程中[5]起重要作用。單核苷酸多態(tài)性(Single Nucleotide Polymorphism，SNP)作為第三代分子標(biāo)記[6]，在遺傳多樣性分析[7]、遺傳連鎖圖譜構(gòu)建[8]、品種鑒定[9]以及性狀的基因定位[10]等研究中應(yīng)用廣泛。近年來，歐洲釀酒葡萄的基因組框架圖和全基因組測(cè)序的完成[11-12]，以及大量EST序列和cDNA序列的公布[13]，為從分子水平研究葡萄的遺傳多樣性提供了大量信息，同時(shí)也為葡萄CXEs基因的深入研究提供了條件。

葡萄是我國重要的栽培產(chǎn)物之一，葡萄病害嚴(yán)重影響葡萄的產(chǎn)量和品質(zhì)[14]，培育抗病葡萄品種是葡萄育種的重要工作，羧酸酯酶在病原菌的致病過程中起重要作用。筆者通過對(duì)釀酒葡萄羧酸酯酶基因片段進(jìn)行克隆和測(cè)序，利用SNP分析蛇龍珠等釀酒葡萄的單核苷酸多態(tài)性及遺傳親緣關(guān)系，為今后研究釀酒葡萄的抗性機(jī)制以及培育抗病性葡萄提供參考。

1 材料與方法

1.1 試驗(yàn)材料 釀酒葡萄8個(gè)蛇龍珠新株系(E01：蓬萊，E02：棲霞北，E03：棲霞南，E04：煙臺(tái)市幸福街道；E05：萊山；E06：龍口；E07、E08：煙臺(tái)開發(fā)區(qū))，品麗珠、赤霞珠、黑比諾材料，由煙臺(tái)張?jiān)Ｆ咸逊N質(zhì)園提供。

1.2 葡萄總DNA的提取 采用Ezup Column Plant Genomic DNA Purification Kit提取基因組總DNA(Sangon,Biotech.Shanghai.Cat.:SK8262)，然后用1%瓊脂糖凝膠電泳檢測(cè)，超微量分光光度計(jì)(P 300)檢測(cè)DNA純度和濃度。

1.3CXEs基因片段的克隆 擴(kuò)增引物C4來自文獻(xiàn)[15]，PCR 反應(yīng)體系50 μL：10×Buffer 5 μL，Mg2+4 μL(25 mmol/L)，dNTPs 4 μL (2.5 mmol/L)，引物1.25 μL(10 μmol/L)，1 UTaq聚合酶，50 ng模板DNA。PCR反應(yīng)程序：94 ℃預(yù)變性2 min；94 ℃變性30 s，55 ℃退火30 s，72 ℃延伸45 s，35個(gè)循環(huán)；72 ℃延伸10 min，4 ℃保存。PCR產(chǎn)物在1.0%瓊脂糖凝膠上電泳分離，凝膠成像系統(tǒng)下觀察。PCR產(chǎn)物純化回收后送至上海生工生物工程有限公司進(jìn)行單向測(cè)序。

1.4CXEs基因片段的序列分析 將測(cè)序結(jié)果在NCBI上的BLAST程序與己登錄的葡萄屬植物進(jìn)行同源性比對(duì)。利用BioEdit version 7.0.5軟件進(jìn)行人工調(diào)整，利用軟件ClustalX 2.0將測(cè)序的11個(gè)樣品進(jìn)行序列比對(duì)，MEGA 5.0進(jìn)行系統(tǒng)發(fā)育分析，采用鄰接法(Neighbour Joining，NJ)構(gòu)建系統(tǒng)發(fā)育樹，以Bootstrap進(jìn)行1 000次可信度分析。利用SNiPlay在線工具(http://sniplay.cirad.fr)對(duì)序列SNP出現(xiàn)的頻率、突變位點(diǎn)以及多態(tài)性[16]進(jìn)行分析，核苷酸多樣度(Pi)、單倍型多樣度(Hd)、中性檢測(cè)運(yùn)用DnaSP 5.0軟件進(jìn)行分析。

2 結(jié)果與分析

2.1 PCR擴(kuò)增結(jié)果 利用引物對(duì)11份葡萄材料的基因組總DNA進(jìn)行PCR擴(kuò)增，在600 bp附近得到1條明亮清晰的條帶(圖1)，測(cè)序后得到葡萄的CXEs基因序列。

2.2CXEs基因片段的序列比對(duì)及氨基酸序列比對(duì) 陽性克隆測(cè)序成功后，利用NCBI上的BLAST對(duì)目的片段進(jìn)行同源性比對(duì)，結(jié)果顯示其與已登錄的葡萄屬CXEs基因序列(基因登錄號(hào)：XM_0106626666.1)的同源性高達(dá)99%，將目的片段產(chǎn)物的氨基酸序列進(jìn)行同源性比對(duì)，結(jié)果顯示其與已登錄的葡萄屬CXEs氨基酸序列(基因登錄號(hào)：XP_010660968.1)的同源率高達(dá)94%，因此確定該段序列為CXEs基因中一段序列。

注：M.DL 10000 DNA marker；1～8.蛇龍珠E01～E08；9.品麗珠；10.赤霞珠；11.黑比諾。Note：M.DL 10000 DNA marker; 1-8.Cabernet Gernischt E01-E08;9.Cabernet Franc;10.Cabernet Sauvignon;11.Pinot Noir.圖1 PCR純化產(chǎn)物凝膠檢測(cè)圖譜 Fig 1 Gel electrophoresis of PCR purification products

注： *表示相似性100%。Note:*：100% identity on sequence.圖2 CXEs基因片段序列比對(duì)Fig.2 The sequences alignment of CXEs gene segments

注：*表示該列殘基完全一致，：表示該列僅包含非常保守的殘基替換，.表示該列含有比較保守的殘基替換。Note：* indicates exactly the column residues ,: indicates that the column contains very conservative residues，.indicates that the column contains relatively conservative residues.圖3 CXEs基因片段氨基酸序列比對(duì)Fig.3 The amino acid sequence alignment of CXEs gene segments

CXEs基因片段的序列比對(duì)結(jié)果見圖2，核苷酸序列有差異，將供試材料與NCBI上黑比諾的CXEs基因進(jìn)行序列比對(duì)，在編碼區(qū)和非編碼區(qū)均存在不同程度的堿基突變。編碼區(qū)在51～319 bp位置，非編碼區(qū)的單核苷酸多態(tài)性位點(diǎn)明顯多于編碼區(qū)，包括堿基的顛換和轉(zhuǎn)換，說明該區(qū)域的遺傳多樣性較豐富。

表1 CXEs基因序列的SNP特征

利用MEGA 5.2軟件將測(cè)序所得的核苷酸序列轉(zhuǎn)換成相應(yīng)的氨基酸序列(圖3)，經(jīng)BioEdit version 7.0.5軟件分析發(fā)現(xiàn)該基因片段起始密碼子是ATG，位于51～53位核苷酸，終止密碼子是TGA，位于317～319位核苷酸，開放讀碼框267 bp，編碼89個(gè)氨基酸。蛇龍珠CXEs基因片段的非編碼區(qū)氨基酸序列表現(xiàn)出較高的多態(tài)性。

2.3CXEs基因片段的SNP特征分析 11個(gè)樣品的克隆片段大小在608～615 bp，其中，A含量在25.0%～26.6%，C含量在22.1%～23.3%，G含量在18.2%～19.2%，T含量在31.9%～33.9%。該基因片段共檢測(cè)出11個(gè)二等位SNPs(表1)?；蛲蛔冾愋桶ɑ蝾崜Q和基因轉(zhuǎn)換，突變位點(diǎn)的等位基因大部分為純合子。在這些二等位SNPs中，內(nèi)含子區(qū)域7個(gè)，外顯子區(qū)域4個(gè)，在外顯子中有1處堿基突變使編碼的氨基酸單核苷酸發(fā)生改變(Q/K)。

2.4CXEs基因片段的遺傳多樣性分析和中性檢驗(yàn) 利用DnaSP 5.0軟件將測(cè)序得到的核苷酸序列進(jìn)行遺傳多樣性分析，結(jié)果表明，其多態(tài)性位點(diǎn)數(shù)87，突變位點(diǎn)數(shù)87，單倍型數(shù)量11，單倍型多樣度(Hd)為1.000±0.039，核苷酸多樣度(Pi)為0.050 76±0.011 49， Tajima′s D、Fu and Li′s D和Fu and Li′s F 3種中性檢驗(yàn)均為正值(0.006 79、0.388 96、0.330 08)，且不顯著，符合中性進(jìn)化模型，說明這些葡萄樣品未發(fā)生種群擴(kuò)張。

2.5 進(jìn)化樹分析 采用MEGA 5.0軟件中的NJ方法構(gòu)建系統(tǒng)進(jìn)化樹(圖4)，11個(gè)樣品被分成2個(gè)分支，4個(gè)亞組。第一個(gè)分支包括黑比諾和赤霞珠。第二個(gè)分支包括2個(gè)亞組，品麗珠為第一亞組，來自不同地區(qū)的8個(gè)蛇龍珠聚為第二亞組。蛇龍珠與品麗珠關(guān)系較近，來自棲霞北的E02和來自棲霞南的E03聚為第一小分支，來自煙臺(tái)開發(fā)區(qū)的E07和E08聚為第二小分支，來自龍口的E06、來自蓬萊的E01、來自煙臺(tái)市幸福街道的E04、來自萊山的E05各為一小分支。來自相同地區(qū)的蛇龍珠E03和E02，蛇龍珠E07和E08各聚為一類，說明供試材料蛇龍珠8個(gè)新株系存在遺傳變異和地理分化，遺傳距離與地理距離一致。

圖4 不同材料CXEs基因片段的系統(tǒng)進(jìn)化樹Fig.4 Phylogenetic tree of CXEs gene from different samples

3 結(jié)論與討論

基因序列數(shù)據(jù)庫的獲得促進(jìn)除農(nóng)作物SNP標(biāo)記模型外其他植物的發(fā)展[17]。大量植物抗病基因已經(jīng)被確定，其中大部分已經(jīng)標(biāo)記了SNP位點(diǎn)[18]。該研究中，該基因片段共發(fā)現(xiàn)87個(gè)突變位點(diǎn)，平均74 bp左右有1個(gè)SNP位點(diǎn)，SNP頻率1.70%～29.00%，Velasco等[12]研究葡萄品種雜合子基因組的高質(zhì)量共識(shí)序列，結(jié)果表明，在葡萄遺傳連鎖圖譜中SNP頻率為4.0%，黑比諾每100 bp有1個(gè)SNP位點(diǎn)。這可能是由于蛇龍珠經(jīng)過長(zhǎng)期的人工培育和自然選擇，其基因組產(chǎn)生更多的SNP位點(diǎn)[19]。這些SNP位點(diǎn)對(duì)于遺傳鑒定、遺傳多樣性分析和親本鑒定具有重要作用。

該試驗(yàn)測(cè)序得到的CXEs基因片段的核苷酸序列有較大差異，共有11個(gè)二等位SNP位點(diǎn)，其中，非編碼區(qū)有7個(gè)SNP位點(diǎn)，編碼區(qū)有4個(gè)SNP位點(diǎn)，非編碼區(qū)的單核苷酸多態(tài)性位點(diǎn)多于編碼區(qū)，而在自花授粉的玉米和棉花的單核苷酸多態(tài)性研究中也得出了“非編碼區(qū)的單核苷酸多態(tài)性位點(diǎn)多于編碼區(qū)”的結(jié)論[20-21]。Dong等[17]研究11個(gè)歐亞種和5個(gè)歐美種的SNP特征,結(jié)果表明，基因轉(zhuǎn)換70%，基因顛換20%。Yang等[22]通過計(jì)算機(jī)輔助分析表達(dá)序列標(biāo)簽(EST)研究番茄的SNP，結(jié)果表明，基因轉(zhuǎn)換58.1%，基因顛換39.5%。而該試驗(yàn)得出基因轉(zhuǎn)換的頻率(71.8%)也高于基因顛換(20.7%)，與前人的研究結(jié)果一致。

通過CXEs基因片段序列差異性構(gòu)建的系統(tǒng)進(jìn)化樹可以有效地將蛇龍珠8個(gè)新株系及品麗珠、赤霞珠和黑比諾進(jìn)行聚類分析，且聚類分析和地理距離具有相關(guān)性；這些葡萄樣品的CXEs基因片段Pi(0.050 76±0.011 49)高于王佳媛等[23]得出的凹葉木蘭Pi(0.017 8)以及Zhu等[24]得出的大豆Pi(0.001 2)，Pi和Hd的數(shù)值越大，其多樣性程度越高，遺傳多樣性越豐富，說明該基因片段具有較高的遺傳多樣性。綜上所述，該基因片段可作為一種遺傳標(biāo)記用于鑒別葡萄的親緣關(guān)系。植物在進(jìn)化過程和多年栽培過程中發(fā)生變異引起個(gè)體之間的差異，這些差異決定植物的品種、特性等[25]。酯酶基因擴(kuò)增及突變是酯酶參與抗性的主要方式[26]，今后可對(duì)抗性基因的單核苷酸進(jìn)行關(guān)聯(lián)研究，為分子標(biāo)記輔助育種開發(fā)可利用的單核苷酸標(biāo)記提供理論依據(jù)。

[1] KO M，CHO J H，SEO H H，et al.Constitutive expression of a fungus-inducible carboxylesterase improves disease resistance in transgenic pepper plants[J].Planta，2016，244(2)：379-392.

[2] 周帥.羧酸酯酶基因的克隆表達(dá)及其應(yīng)用[D].長(zhǎng)春：長(zhǎng)春理工大學(xué)，2013.

[3] 吳麗麗.雌激素對(duì)小鼠肝臟羧酸酯酶的作用及其機(jī)制[D].南京：南京醫(yī)科大學(xué)，2015.

[4] GERSHATER M C，EDWARDS R.Regulating biological activity in plants with carboxylesterases[J].Plant Sci，2008，173(6)：579-588.

[5] SALERNO M I，GIANINAZZI S，ARNOULD C，et al.Ultrastructural and cell wall modifications during infection ofEucalyptusviminalisroots by a pathogenicFusariumoxysporumstrain[J].Journal of general plant pathology，2004，70(3)：145-152.

[6] 劉穎，朱方何，洪彥彬，等.單核苷酸多態(tài)性的研究進(jìn)展與應(yīng)用[J].廣東農(nóng)業(yè)科學(xué)，2011(S1)：50-53.

[7] INGHELANDT D V，MELCHINGER A E，LEBRETON C，et al.Population structure and genetic diversity in a commercial maize breeding program assessed with SSR and SNP markers[J].Theoretical & applied genetics，2010，120(7)：1289-1299.

[8] TREBBI D,MACCAFERRI M,HEER P D，et al.High-throughput SNP discovery and genotyping in durum wheat (TriticumdurumDesf.)[J].Theoretical & applied genetics,2011，123(4)：555-569.

[9] CIARMIELLO L F,PICCIRILLO P,PONTECORVO G，et al.A PCR based SNPs marker for specific characterization of English walnut (JuglansregiaL.) cultivars[J].Mol Biol Rep，2011，38(2)：1237-1249.

[10] SINGH A,SINGH P K,SINGH R，et al.SNP haplotypes of theBADH1 gene and their association with aroma in rice (OryzasativaL.)[J].Molecular breeding，2010，26(2)：325-338.

[11] JAILLON O,AURY J M,NOEL B，et al.The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla[J].Nature，2007，449(7161)：463-467.

[12] VELASCO R,ZHARKIKH A,TROGGIO M，et al.A high quality draft consensus sequence of the genome of a heterozygous grapevine variety[J].PLoS One，2007，2(12)：1326.

[13] BERGER R,HOFFMANN M，KELLER U.Molecular analysis of a gene encoding a cell-bound esterase fromStreptomyceschrysomallus[J].J Bacteriol，1998，180(23)：6396-6399.

[14] 齊慧霞，張京政，齊永順,等.釀酒葡萄園病蟲害調(diào)查及綜合防治研究[J].北方園藝，2007(8)：202-204.

[15] 郭大龍，李猛，張國海,等.葡萄SNP標(biāo)記的CAPS和TSP分析[J].園藝學(xué)報(bào)，2013，40(11)：2307-2315.

[16] DEREEPER A，NICOLAS S，CUNFF L L，et al.SNiPlay：A web-based tool for detection，management and analysis of SNPs.Application to grapevine diversity projects[J].BMC Bioinformatics，2011，12(1)：406-411.

[17] DONG Q H,CAO X,YANG G，et al.Discovery and characterization of SNPs inVitisviniferaand genetic assessment of some grapevine cultivars[J].Scientia horticulturae，2010，125(3)：233-238.

[18] GASPERO G D,CIPRIANI G,ADAM-BLONDON A F，et al.Linkage maps of grapevine displaying the chromosomal locations of 420 microsatellite markers and 82 markers for R-gene candidates[J].Theoretical & applied genetics，2007，114(7)：1249-1263.

[19] 李紅娟，周新明，劉勇強(qiáng)，等.‘蛇龍珠’親緣關(guān)系鑒定在營養(yǎng)系選種中的應(yīng)用[J].中國農(nóng)學(xué)通報(bào)，2012，28(28)：153-157.

[20] CHING A，CALDWELL K S，JUNG M，et al.SNP frequency，haplotype structure and linkage disequilibrium in elite maize inbred lines[J].BMC Genet，2002，3(7):19.

[21] AN C，SAHA S，JENKINS J N，et al.Cotton (Gossypiumspp.) R2R3-MYB transcription factors SNP identification，phylogenomic characterization，chromosome localization，and linkage mapping[J].Theoretical & applied genetics，2008，116(7)：1015-1026.

[22] YANG W，BAI X，KABELKA E，et al.Discovery of single nucleotide polymorphisms inLycopersiconesculentumby computer aided analysis of expressed sequence tags[J].Molecular breeding，2004，14(1)：21-34.

[23] 王佳媛，吳傳芳，唐亞.基于SNP分子標(biāo)記的凹葉木蘭遺傳多樣性初步研究[J].廣西植物，2012，32(4)：542-547.

[24] ZHU Y L，SONG Q J，HYTEN D L，et al.Single-nucleotide polymorphisms in soybean[J].Genetics，2003，163(3)：1123-1134.

[25] 陳吉寶，趙麗英，褚學(xué)英，等.植物中單核苷酸多態(tài)性的分布及其應(yīng)用[J].安徽農(nóng)業(yè)科學(xué)，2009，37(31)：15153-15154.

[26] 曲明靜，許新軍，韓召軍,等.二化螟酯酶基因克隆及序列分析[J].江西農(nóng)業(yè)學(xué)報(bào)，2006，18(3)：52-54.

The SNP Analysis of Carboxylesterases in Cabernet Gernischt (Vitis vinifera L.)and Other Wine Grape Cultivars

CHEN Xiao-yu1，LIU Lin-de1*，GE Yi-he,LI Ji-ming2*et al

(1.School of Life Science，Ludong University，Yantai，Shandong 264025;2.Shandong key Laboratory of Microbial Fermentation Technology of wine,Yantai Zhangyu Group Co. Ltd.,Yantai Shandong 264001)

[Objective] The aim was to analyze mononucleotide polymorphism and genetic relationships of Cabernet Gernischt (VitisviniferaL.)and other wine grape cultivars，provide reference for study on resistance mechanism of wine grape and cultivation of disease resistant grape.[Method] The difference among Carboxylesterase gene fragments was analyzed to investigate the genetic diversity of Cabernet Gernischt and other wine grape cultivars.Eight new strains of Cabernet Gernischt (E01-E08)，Cabernet Sauvignon，Cabernet Franc were used as experimental materials and Pinot Noir (V.vinifera) were taken as control materials，Carboxylesterase gene fragments were cloned and analyzed.[Result] The analysis of Carboxylesterase gene fragments sequence alignment showed that there were different degrees of base substitutions between the samples in the coding and non-coding regions.A total of 87 SNPs were detected in the sequenced fragments，a SNP locus was found on every 4 bp sequence on average.Nucleotide diversity (Pi) =0.050 76 (±0.011 49)，Haplotype diversity (Hd)=1.000 (±0.039) which give a high leve1 of polymorphism.Three neutrality tests were used to compare the patterns of sequence variations，the result indicates conform to neutral mutation.The clustering dendrogram was constructed by Neighbour Joining (NJ)，Cabernet Gernischt E02 (north of qixia) and E03 (south of qixia) classified into one group，E07(development zone) and E08 (development zone) classified into one group，genetic distance in accord with geographic distance.The results of the experiment showed that there were several mutations in the Carboxylesterase gene fragments might be selected to be the molecular marker to distinguish Cabernet Gernischt eight new strains and other wine grape varieties.[Conclusion]The mutation site of CXEs gene fragment can be used as a kind of genetic marker,SNP method can identify genetic relationships of grape.

Cabernet Gernischt; Carboxylesterase; Gene cloning; Sequence analysis

煙臺(tái)市科技發(fā)展計(jì)劃項(xiàng)目(2012JH204)；泰山產(chǎn)業(yè)領(lǐng)軍人才工程專項(xiàng)；國家高技術(shù)研究發(fā)展計(jì)劃項(xiàng)目(2013AA102108)。

陳小宇(1991- )，女，山東德州人，碩士研究生，研究方向:種群群落生態(tài)學(xué)。*通訊作者，劉林德，教授，博士，碩士生導(dǎo)師，從事種群生態(tài)學(xué)和遺傳學(xué)研究；*通訊作者，李記明，研究員，博士，博士生導(dǎo)師，從事葡萄和葡萄酒方面的研究。

2016-08-11

S 188

A

0517-6611(2016)28-0109-05