浙江雁荡山山脉叶栖酵母菌资源与物种多样性

臧　威,刘春月,,谢广发,沈　赤,刘新展,邹彗君,汤访评,冯富娟,孙剑秋,*,白逢彦

1 绍兴文理学院生命科学学院,绍兴　312000 2 东北林业大学生命科学学院,哈尔滨　150040 3 浙江古越龙山绍兴酒股份有限公司,绍兴　312000 4 中国科学院微生物研究所真菌学国家重点实验室,北京　100101

酵母菌只是形态学术语而不是分类学名词,是一类主要以芽殖或裂殖进行无性繁殖、有性繁殖阶段不形成子实体的单细胞真菌[1]。在植物、土壤、淡水、海水等生境中广泛分布着各类酵母菌,与人类生活关系密切[2]。基于酵母菌在工业、农业、医学、食品等行业和基础研究、生态功能等领域的重要作用[1],考察和分离各种生境中存在的酵母菌资源是生物技术发展的重要基础[3],所以特定生境内酵母菌类群的分布情况一直受到广泛关注。对于酒类酿造[4- 7]、传统食品发酵[8- 10]、水环境[11]、土壤[12]、植物表面[3,13- 15]、人体[16]和患处[17- 18]、极端环境[19]及其他生境[20]中酵母菌资源的研究成果也充分说明,自然界中还有大量的酵母菌种类亟待发现。

雁荡山山脉位于浙江省东南部,绵延数百公里,按照地理位置可以分为北雁荡山、中雁荡山、南雁荡山、西雁荡山(泽雅)、东雁荡山(洞头半屏山),气候温和湿润,适合酵母菌栖息繁殖,但是迄今对该地区野生酵母菌资源研究未见报道。我们采用经典的酵母菌分离培养方法和现代酵母菌分子鉴定技术,对雁荡山山脉广阔区域内的叶栖酵母菌资源进行考察,旨在弄清该地区叶栖酵母菌群落的多样性、组成和结构,发现特殊酵母菌类群以进一步丰富我国酵母菌资源库。

1　材料与方法

1.1　采样地点

研究范围集中在浙江省东南部雁荡山地区,包括北雁荡山、中雁荡山、南雁荡山、西雁荡山(泽雅)、东雁荡山(洞头半屏山),该地区属亚热带海洋性气候,雨量充沛、气候温暖,年平均气温13.5℃,年平均降雨量1935.6 mm,年平均无霜期269 d。

1.2　采样方法

在2015—2016年间,自雁荡山山脉采集树木半枯叶片暂存于无菌自封袋内,尽快运回实验室进行叶栖酵母菌的分离和纯化。采样过程中,对宿主植物信息、周围环境、纬度和经度等进行详细记载。

1.3　酵母菌分离

参照Nakase和Takashima方法[21]进行。将采集的半枯叶片,用凡士林粘贴在PDA平板的皿盖内侧,在20℃下正置恒温培养7 d,期间每天更换新鲜培养基。更换下的PDA平板在20℃下倒置培养至渐见菌落形成,挑取形态不同的酵母菌菌落进行平板划线,即可获得叶栖酵母菌的纯培养。

1.4　酵母菌DNA提取

挑取少量在PDA斜面上生长旺盛的幼龄酵母菌菌体进行DNA微量提取,参照Makimura方法[22]完成。

1.5　酵母菌26S rDNA D1/D2区PCR扩增

参考Kurtzman和Robnett方法[23],采用引物NL4(5′-GGTCCGTGTTTCAAGACGG- 3′)和NL1(5′-GCATATCAATAAGCGGAGGAAAAG- 3′)对菌株26S rDNA D1/D2进行扩增,按照反应条件94℃ 1 min,52℃ 1 min,72℃ 1 min扩增30个循环。本研究中引物NL4和NL1的合成由生工生物工程(上海)股份有限公司完成,北京六合华大基因科技有限公司对扩增产物进行测序。

1.6　序列分析

将测得酵母菌26S rDNA D1/D2序列在GenBank核酸序列数据库中进行同源序列搜索(BLAST search),比较实验菌株与已知菌株序列的相似程度。采用Clustal X软件比对相应模式菌株序列与实验菌株26S rDNA D1/D2序列的碱基差异[24]。通过MEGA5软件和Neighbor-Joining方法,构建分子系统树以明确相关酵母菌的分类地位[25- 26]。根据系统树上标注的序列号(accession numbers),可以从GeneBank下载相关菌株26S rDNA D1/D2参考序列。

1.7　数据分析

相对频率(Relative frequency, RF)是指某雁荡山地区分离到的某种叶栖酵母菌占该地区叶栖酵母菌总菌株数的百分数,若≥1%则属于常见种类。

根据以下Shannon-Weiner公式计算多样性指数(Diversity index,H′),式中k是指雁荡山某地区分离到的常见叶栖酵母菌种类的总数,Pi是指某种叶栖酵母菌的菌株数量占该地区叶栖酵母菌全部菌株数的百分数。

丰度(Richness, R),是指某雁荡山地区分离到的常见叶栖酵母菌种类数。

相似性系数(Similarity coefficients,Cs),根据Sorenson系数公式进行计算,Cs=2j/(a+b),式中j是指雁荡山山脉两个地区之间都可以分离到的常见叶栖酵母菌种类数,a是指雁荡山山脉一个地区可以分离到的常见叶栖酵母菌种类数,b是指雁荡山山脉另一个地区可以分离到的常见叶栖酵母菌种类数。

采用bootstrap指数获得期望物种丰富度(S),在如下计算公式中So代表雁荡山某地区分离到的常见叶栖酵母菌种类数,Pk代表某种叶栖酵母菌的菌株数占该地区叶栖酵母菌全部菌株数的百分数,m是指全部样本数量。借助EstimateS软件完成分析和计算,通过Excel绘制雁荡山山脉叶栖酵母菌的物种累积曲线(Species accumulation curves)。

2　结果与分析

2.1　雁荡山山脉叶栖酵母菌的已知类群

从北雁荡山、中雁荡山、南雁荡山、西雁荡山和东雁荡山分别采集到树木半枯叶片标本30、32、30、31和30份,分离获得酵母菌144、209、257、337株和171株。参照1118株酵母菌26S rDNA D1/D2序列在GenBank核酸序列数据库中同源序列的搜索(BLAST search)结果,采用Clustal X软件将相应模式菌株与实验菌株26S rDNA D1/D2序列进行比较分析并选择典型菌株绘制分子系统树,其中1072个菌株与已知模式菌株26S rDNA D1/D2序列相似性在99%以上,被鉴定为37属56种(图1)。

图1　基于26S rDNA D1/D2序列构建雁荡山山脉典型酵母菌株与近缘种间的系统进化树Fig.1　The phylogenetic tree based on 26S rDNA D1/D2 sequences, showing the relationships among representative strains and their close relatives编号中E、S、W、N和M分别表示该菌株分离自东雁荡山、南雁荡山、西雁荡山、北雁荡山和中雁荡山

2.2　雁荡山山脉叶栖酵母菌的疑似新种

在雁荡山山脉分离得到1118个叶栖酵母菌株中,有46个菌株26S rDNA D1/D2序列较为特殊,与最近缘的已知种类26S rDNA D1/D2序列均存在1%以上的碱基差异,根据目前国际酵母菌分类标准[27- 29]疑似属于新的酵母菌类群。其中,分离自北雁荡山共有23株、中雁荡山7株、南雁荡山10株、西雁荡山2株和东雁荡山4株。通过Clustal X软件比较并与相应模式菌株比对,包括子囊菌酵母4属5种和担子菌酵母9属14种(表1)。

表1　疑似新种与最近缘已知酵母菌模式种26S rDNA D1/D2序列比较

续表疑似新种Proposednewspecies菌株Strains最近缘模式种(序列号)与相似性/%Closestrelatives(accessionnumber)andsequenceidentity宿主植物Hosts纬度和经度Latitudeandlongitude地点SitesColacogloeasp．2N0921C．foliorum(KY106943)96．15Cryptomeriafortunei28°21′47″121°3′59″北雁荡Colacogloeasp．3N2121C．eucalyptica(EU075184)97．11Schimasuperba28°22′2″121°5′26″北雁荡Derxomycessp．S0372D．simaoensis(KY107627)97．28Sloaneasinensis27°36′6″120°16′14″南雁荡S1412Elaeocarpusdecipiens27°35′57″120°16′22″Fellozymasp．N3051F．inositophila(KY107699)93．19Ligustrumsp．28°22′5″121°5′21″北雁荡S095293．01Phoebechekiangensis27°35′56″120°16′20″南雁荡Kockovaellasp．S0541K．barringtoniae(AB292854)98．04Camelliaoleifera27°36′4″120°16′11″南雁荡Kondoasp．1M151K．myxariophila(AF189904)96．10Rosalaevigata？中雁荡E1511Rhuschinensis27°49′52″121°7′11″东雁荡E2111Hoveniaacerba27°50′5″121°7′13″E2711Ficuserecta27°50′8″121°7′16″Kondoasp．2S1322K．myxariophila(AF189904)97．37Syzygiumjambos27°35′58″120°12′23″南雁荡S1512Euonymuscarnosus27°35′57″120°16′22″Saitozymasp．W0951S．flavus(AF075497)91．83Adinarubella28°2′11″120°22′20″西雁荡N0931Cryptomeriafortunei28°21′47″121°3′59″北雁荡N0941N0961N0991

2.3　雁荡山山脉叶栖酵母菌的分布特征

从雁荡山山脉分离到的1118株叶栖酵母菌被鉴定为37属56种已知种类和13属19种疑似新种(图1和表1)。其中,在东雁荡山有12种常见类群,占该地区叶栖酵母菌的95.4%；在南雁荡山有14种常见类群,占该地区叶栖酵母菌的86.8%；在西雁荡山有9种常见类群,占该地区叶栖酵母菌的94.4%；在北雁荡山有14种常见类群,占该地区叶栖酵母菌的93.3%；中雁荡山分离到13种常见种类,占该地区叶栖酵母菌的94%(表2)。雁荡山山脉常见的叶栖酵母菌种类中,只有Derxomycesmrakii和Derxomycespseudoschimicola两种酵母菌在东雁荡山、西雁荡山、南雁荡山、北雁荡山、中雁荡山都可以分离到；除北雁荡山的优势酵母为Derxomycespseudoschimicola之外,其他地区分离到的叶栖酵母菌中Derxomycesmrakii均表现出明显的分布优势；特别是在东雁荡山,Derxomycesmrakii在该地区常见叶栖酵母菌中的相对频率可以达到55% (表2)。

表2　雁荡山山脉常见叶栖酵母菌的相对频率/%

稀有种类; The rare taxa (RF﹤1%);

东雁荡山East Yandang mountain:Bulleribasidiumfoliicola,Bulleribasidiumsp.1,Ceraceosorusafricanus,Erythrobasidiumsp.,Hannaellasinensis,Microsporomycesmagnisporus,Naganishiaalbida,Tilletiopsissp.

南雁荡山South Yandang mountain:Bannoahahajimensis,Chrysozymasp.1,Chrysozymasp.2,Colacogloeafalcata,Colacogloeasp.1,Cystobasidiomycetessp.,Derxomycesboekhoutii,Derxomycessp.,Dioszegiazsoltii,Exobasidiumsymploci-japonicae,Fellozymasp.,Golubeviapallescens,Kockovaellasp.,Kockovaellavietnamensis,Kondoasp.2,Microsporomycesmagnisporus,Microstromaphylloplanum,Nielozymaformosana,Phyllozymacorallina,Rhodotorulasp.,Sporobolomyceskoalae,Sporobolomycespararoseus,Symmetrosporacoprosmae,Takashimellaformosensis.

西雁荡山West Yandang mountain:Bannoaogasawarensis,Bulleribasidiumsp.2,Cryptococcussp.,Derxomycesbambusicola,Dioszegiazsoltii,Hannaellasurugaensis,Myriangiumsp.,Oberwinklerozymayarrowii,Piskurozymacylindricus,Rhodotorulasp.,Saitozymasp.

北雁荡山North Yandang mountain:Carlosrosaeasp.,Colacogloeafalcate,Colacogloeasp.2,Colacogloeasp.3,Derxomycesboekhoutii,Fellozymasp.,Meiranashicola,Moesziomycesaphidis,Papiliotremaaurea,Sporobolomyceskoalae.

中雁荡山Middle Yandang mountain:Collophorasp.,Coniosporiumsp.,Cryptococcussp.,Derxomycesboekhoutii,Dioszegiazsoltii,Gjaerumiaminor,Jaminaeaangkorensis,Kondoasp.1,Oberwinklerozymayarrowii,Papiliotremaaurea,Rhodotorulasp.,Tilletiopsissp.

多样性指数(H′)和丰度(R)可以反映该地区叶栖酵母菌的物种多样性程度,在雁荡山山脉广阔地区分布的叶栖酵母菌具有丰富的多样性。从多样性指数(H′)来看,北雁荡山﹥中雁荡山﹥南雁荡山﹥西雁荡山﹥东雁荡山；根据酵母菌的丰度值,也可以看出分离自北雁荡山、中雁荡山和南雁荡山的叶栖酵母菌种类相对更多(表3)。

表3　雁荡山山脉常见叶栖酵母菌的多样性指数(H′)和丰度(R)

相似性系数可以比较雁荡山山脉两个地区之间叶栖酵母菌物种组成的相似性程度。根据表4发现,中雁荡山常见叶栖酵母菌组成与东、南、西、北雁荡山等周边地区的相似性系数较高,达到0.4—0.52；而南雁荡山与西雁荡山、中雁荡山常见叶栖酵母菌组成最为相近,相似性系数达到0.52的较高水平；东、西雁荡山和南、北雁荡山的叶栖酵母菌种类的相似程度则较低,相似性系数分别为0.19和0.21。

表4　雁荡山山脉常见叶栖酵母菌的相似性系数

2.4　雁荡山山脉叶栖酵母菌的序列差异

以东雁荡山、西雁荡山、南雁荡山、北雁荡山和中雁荡山都可以分离到的Derxomycesmrakii为代表,对雁荡山山脉属于同一种类的叶栖酵母菌物种资源D1/D2序列进行分析和比较。结果发现,分离自雁荡山山脉的广布种Derxomycesmrakii包括13种基因型,根据每种基因型与模式菌株CBS8288的D1/D2序列比较结果及其在雁荡山不同地区分布的相对频率(表5),其中与模式菌株D1/D2序列完全相同的基因型01最常见,在雁荡山各个地区可以占到Derxomycesmrakii菌株的87.4%—98.2%。

表5　Derxomyces mrakii各种基因型的相对频率及其与模式菌株的序列差异/%

2.5　雁荡山山脉叶栖酵母菌的物种累积曲线

图2　雁荡山山脉叶栖酵母菌的物种累积曲线Fig.2　The species accumulation curves for the yeasts from the Yandang Mountain area1：物种累积曲线；2、3：在95%置信区间内物种累积曲线的下限与上限95%

基于雁荡山山脉常见叶栖酵母菌在东雁荡山、西雁荡山、南雁荡山、北雁荡山和中雁荡山的分布数据(表2),采用EstimateS软件对该地区叶栖酵母菌多样性进行了分析和计算,并绘制出雁荡山山脉叶栖酵母菌的物种累积曲线(图2)。

根据物种累积曲线特征和EstimateS软件计算结果,其Bootstrap值为46.95,估计在本项研究中已经将雁荡山地区80.94%的叶栖酵母菌分离出来。

3　讨论

雁荡山是我国华东地区的重要山脉之一,幅员辽阔、资源丰富,在生物多样性研究方面受到广泛重视[30]。该地区温度、湿度等环境因子适合野生酵母菌栖息,系统、全面的研究雁荡山地区酵母菌资源情况具有一定代表性。在野生酵母菌研究方面,植物叶表是人们比较关注的酵母菌栖息环境[31- 32],而且每年秋季的成熟植物叶片上酵母菌的种类和数量最为丰富[33]。现在一般认为,同土壤等生境相比较,叶表可能更容易分离到特殊种类的酵母菌。

目前国际上通行的酵母菌鉴定依据主要是26S rRNA基因D1/D2区域的序列组成[23,27],种内不同菌株之间D1/D2序列差异应该在1%以内,如果株间该序列碱基差异在1%以上则基本可以判定为不同种类[28]。近年来酵母菌分类系统变化较大[34- 36],在进行雁荡山地区叶栖酵母菌的序列比对和系统发育分析时,参照最新的酵母菌分类系统对鉴定结果进行了必要的修正。

分离自雁荡山山脉的1072株叶栖酵母菌,经过26S rRNA基因D1/D2序列比对和系统发育分析,属于56个已知类群(图1)；另外还有46个菌株26S rDNA D1/D2序列与已知种类差异明显,根据26S rRNA基因D1/D2序列比对和系统发育分析结果,基本可以确定这些菌株代表19个新类群(表1)；尽管有些菌株,还需要详细研究其ITS序列特征和生理生化特征、细胞形态特征、培养特征以进一步证明其新种地位,但是疑似新种比较常见的事实也说明,研究雁荡山地区叶栖酵母菌资源具有重要的理论意义。

在雁荡山地区有38种常见叶栖酵母菌类群(表2),包括Derxomyces属5种,Sporobolomyces属3种,Bannoa、Bulleribasidium、Rhodosporidiobolus、Symmetrospora、Taphrina、Tilletiopsi和Udeniomyces属各2种,Bullera、Coniochaeta、Coniosporium、Cryptococcus、Elsinoe、Erythrobasidium、Fellozyma、Golubevia、Kockovaella、Kondoa、Leucosporidium、Moesziomyces、Oberwinklerozyma、Phyllozyma、Ruinenia和Saitozyma属各1种,其中Derxomycesmrakii在该地区分离到的叶栖酵母菌中具有明显分布优势。Hamamoto和Nakase[37]最早从新西兰的三种植物Leptospermumscoparium、Nothofagusfusca、Pseudowinteracolonata上分离到Derxomycesmrakii培养物,并将其鉴定为Bullera属新种,Wang和Bai[38]在建立Derxomyces属时将其重合组合在Derxomyces属中。在分析比较雁荡山地区Derxomycesmrakii的D1/D2序列时,发现大量培养物与模式酵母的序列并非完全一致,表现出一定程度的遗传多样性(表5)。

雁荡山山脉叶栖酵母菌资源丰富(表3),根据不同地区常见叶栖酵母菌组成(表2)和相似程度(表4),发现中雁荡山与东雁荡山、西雁荡山、南雁荡山和北雁荡山叶栖酵母菌相似性系数达到0.4—0.52,表明中雁荡山与其他4个地区的叶栖酵母菌种类组成更加相似；而东、西雁荡山之间和南、北雁荡山之间叶栖酵母菌组成相似性系数分别为0.19和0.21,酵母菌种类组成存在更大差异；可能地理上相近地区的酵母菌物种分布具有更大相似性,这个推测可以供今后野生酵母菌资源考察工作参考。

在生物多样性研究领域,物种累积曲线被广泛用于样本量充分性的判断以及物种丰富程度的估计[39]。根据图2,我们已经分离到雁荡山地区绝大多数的叶栖酵母菌种类,研究结果可以真实反映出该地区叶栖酵母菌资源和物种多样性的存在状况。