基于分子标记技术解析竹黄及其分离真菌的遗传多样性

基于分子标记技术解析竹黄及其分离真菌的遗传多样性一、引言1.1竹黄研究背景竹黄（ShiraiabambusicolaEHenn.），又名竹花、竹赤团子、竹三七、竹茧等，是一种专性寄生于某些竹子嫩枝上的珍稀寄生性药用真菌。其主要分布于中国和日本，在我国，主要集中在长江以南的四川、江西、江苏、浙江、湖北、湖南、安徽、贵州、福建和云南等省份。竹黄多生长于衰败或即将衰败的竹林中，常见寄主为箣竹属（BambusaSieb.etFranch）等竹子。竹黄在传统医学中应用历史悠久，性温味淡，具备止咳祛痛、舒筋活络、祛风利湿、补中益气、活血补血、散瘀通经等功效。民间常用竹黄浸酒，用以治疗风湿性关节炎、坐骨神经痛、腰肌劳损、跌打损伤、虚寒胃痛、小儿惊风及急性肝炎等病症。现代研究进一步揭示了竹黄丰富的药用价值，其有效成分包括竹红菌素、竹黄多糖、11,11'-二去氧沃替西林和蒽酮类化合物等。其中，竹红菌素作为目前已知的在可见光区内的优良光敏剂，具有显著的抗肿瘤、抗菌、抗炎等作用；竹黄多糖则表现出护肝等功效。此外，竹黄还具备抗病毒、镇痛、局部麻醉、护心血管、利尿等作用。近年来，随着对竹黄药用价值研究的不断深入，其市场需求日益增长。然而，由于大规模人工采集以及原生地生境遭到破坏，竹黄的生物多样性资源迅速减少，保护竹黄资源已成为当务之急。深入研究竹黄及其分离真菌的遗传多样性，不仅有助于我们更好地了解竹黄的遗传背景、进化历程和生态适应性，还能为竹黄的资源保护、品种改良以及可持续开发利用提供坚实的理论基础。同时，对于揭示竹黄与寄主竹子之间的相互作用关系，以及探索新的药用活性成分和生物防治策略等方面，也具有重要的科学意义。1.2研究目的与意义本研究旨在运用先进的分子生物学技术，全面、系统地分析竹黄及其分离真菌的遗传多样性，深入探究其群体结构和遗传变异规律。具体目标如下：其一，通过对不同地理居群竹黄样本的遗传多样性分析，明确竹黄的遗传分化情况，揭示地理因素对其遗传结构的影响，为竹黄种质资源的保护和管理提供科学依据；其二，从竹黄子实体中分离纯化真菌，并对其进行遗传多样性分析和系统学鉴定，明确与竹黄相关的真菌种类及其亲缘关系，为深入研究竹黄的生物学特性和生态功能提供参考；其三，筛选出具有优良性状的竹黄分离真菌菌株，为竹黄的人工栽培和产业化开发提供优质菌种资源。本研究具有重要的理论和实践意义。在理论方面，丰富了竹黄及其相关真菌的遗传多样性研究内容，有助于完善对竹黄这一珍稀药用真菌的生物学认知，推动真菌遗传学和生态学的发展。在实践方面，为竹黄资源的合理保护和可持续利用提供科学指导，通过筛选优质菌株，有望提高竹黄的人工栽培效率和质量，满足市场对竹黄日益增长的需求，促进竹黄产业的健康发展。同时，对于开发新型药物、生物防治剂等也具有潜在的应用价值，为人类健康和生态环境保护做出贡献。二、竹黄及其分离真菌概述2.1竹黄的生物学特性2.1.1形态特征竹黄的子座呈现出不规则的瘤状，在生长初期，其颜色为白色，质地较为柔软，表面平滑且富有光泽，宛如新生的嫩芽，充满生机。随着生长的推进，子座逐渐变为粉红色，此时质地也从柔软的肉质转变为稍硬的木栓质，表面开始出现龟裂，如同岁月在其身上留下的痕迹，这一变化通常发生在子座生长的中后期。其大小一般长1.5-4cm，宽1-2.5cm，形态虽不规则，但在自然环境中却有着独特的美感。子囊壳近球形，如同微小的球体，紧密地埋生于子座内部，直径在480-580μm之间。子囊则呈长圆柱状，长度为（280-340）μm，宽度为(22-35)μm，宛如细长的容器，整齐地排列在子座之中。子囊孢子单行排列，形状为长方形至梭形，两端大多尖锐，仿佛细长的针，且具有纵横隔膜，大小为（42-92）μm×(13-35)μm，在未成熟时无色或近无色，当它们成堆聚集时，则呈现出柿黄色，在显微镜下观察，这些孢子犹如排列有序的微小宝石，散发着独特的魅力。在不同的生长阶段，竹黄的形态变化显著。生长初期的白色肉质子座，是竹黄生长的起点，此时它对环境的要求较为苛刻，需要适宜的温度、湿度和光照条件，以保证其正常的生长和发育。随着时间的推移，子座颜色的转变和质地的变化，是其生长成熟的标志，这一过程不仅受到自身生长规律的调控，还与外界环境因素密切相关。例如，在温度适宜、湿度较高的环境中，子座的生长速度会加快，颜色转变也会更加迅速；而在光照过强或过弱的情况下，子座的形态和颜色可能会受到影响，出现生长不良或颜色异常的情况。竹黄的形态与环境之间存在着紧密的关联。竹黄多生长于衰败或即将衰败的竹林中，这种特殊的生境为竹黄的生长提供了独特的条件。衰败的竹林中，竹子的生理状态发生改变，其营养成分和代谢产物的分布也会相应变化，这些变化为竹黄的寄生和生长创造了适宜的环境。竹林的湿度、温度和光照等微环境因素，对竹黄的形态和生长发育也有着重要的影响。在阴暗潮湿、水沟旁或者低洼地段的竹林，由于湿度较高，竹黄的子座生长更为茂盛，颜色也更加鲜艳；而在阳光直射较多的地方，竹黄的子实体发生较少，这是因为过强的光照可能会影响竹黄的生理代谢，导致其生长受到抑制。2.1.2分布与寄主竹黄在地理分布上具有一定的局限性，主要集中在亚洲地区，其中中国和日本是其主要的分布区域。在中国，竹黄广泛分布于长江以南的诸多省份，如四川、江西、江苏、浙江、湖北、湖南、安徽、贵州、福建和云南等。这些地区的气候条件较为温暖湿润，年平均气温在15℃-25℃之间，年降水量在1000mm-2000mm之间，为竹黄的生长提供了适宜的环境。长江以南地区丰富的竹子资源，也为竹黄的寄生提供了充足的寄主。竹黄的常见寄主主要为箣竹属（BambusaSieb.etFranch）等竹子。其中，短穗竹、箭竹、毛环短穗竹等竹子品种是竹黄较为常见的寄主。不同寄主对竹黄的生长和遗传特性可能会产生一定的影响。寄主竹子的生理特性，如竹子的营养成分含量、代谢产物种类和含量等，会影响竹黄的生长速度和发育状况。研究表明，在营养成分丰富、代谢活性较高的寄主竹子上，竹黄的生长速度更快，子座更为饱满；而在营养成分匮乏、生长状态不佳的寄主竹子上，竹黄的生长可能会受到抑制，子座较小且发育不完全。寄主竹子的遗传背景也可能与竹黄的遗传特性相互作用，影响竹黄的遗传多样性和适应性。不同寄主竹子的基因表达差异，可能会导致竹黄在寄生过程中发生基因水平的变化，从而影响竹黄的遗传特征和生态适应性。2.1.3药用价值竹黄富含多种活性成分，这些成分赋予了竹黄丰富的药用功效。其主要活性成分包括竹红菌素、竹黄多糖、11,11'-二去氧沃替西林和蒽酮类化合物等。竹红菌素是竹黄中最为重要的活性成分之一，作为目前已知的在可见光区内的优良光敏剂，具有显著的抗肿瘤活性。其作用机制主要是通过光动力作用，在光照条件下，竹红菌素能够吸收光能，产生单线态氧等活性氧物质，这些活性氧物质可以破坏肿瘤细胞的细胞膜、线粒体等细胞器，诱导肿瘤细胞凋亡，从而达到抑制肿瘤生长的目的。竹红菌素还具有抗菌作用，能够抑制多种细菌和真菌的生长繁殖，其抗菌机制可能与破坏微生物的细胞膜结构和干扰其代谢过程有关；同时，它还具有抗炎作用，能够抑制炎症介质的释放，减轻炎症反应，对多种炎症相关疾病具有潜在的治疗作用。竹黄多糖也是竹黄的重要活性成分之一，具有护肝等功效。研究表明，竹黄多糖可以通过调节肝脏细胞的代谢功能，增强肝脏的抗氧化能力，减轻氧化应激对肝脏的损伤，从而起到保护肝脏的作用。竹黄多糖还可能参与调节免疫系统，增强机体的免疫力，对维持机体的健康状态具有重要意义。在传统医学中，竹黄性温味淡，被广泛应用于多种病症的治疗。它具有止咳祛痛、舒筋活络、祛风利湿、补中益气、活血补血、散瘀通经等功效。民间常用竹黄浸酒，用以治疗风湿性关节炎、坐骨神经痛、腰肌劳损、跌打损伤、虚寒胃痛、小儿惊风及急性肝炎等病症。现代医学研究进一步证实了竹黄在这些方面的药用价值，为其临床应用提供了科学依据。2.2竹黄分离真菌的研究现状从竹黄中分离真菌的常用方法主要包括组织分离法、稀释涂布平板法和富集培养法等。组织分离法是将竹黄子实体的组织块直接接种到培养基上，通过培养使组织块上的真菌生长繁殖，从而分离得到纯菌株。这种方法操作相对简单，能够直接从竹黄组织中获取真菌，但容易受到杂菌污染，需要在无菌条件下严格操作。稀释涂布平板法是将竹黄样品进行梯度稀释，然后将稀释液涂布在培养基平板上，使单个真菌细胞在平板上生长形成菌落，再通过挑取单菌落进行纯化培养，获得纯菌株。该方法能够有效地分离出不同种类的真菌，并且可以通过菌落形态和特征初步判断真菌的种类，但操作过程较为繁琐，需要一定的实验技巧。富集培养法是利用特定的培养基和培养条件，对竹黄样品中的目标真菌进行富集培养，增加目标真菌的数量，然后再通过其他分离方法进行分离。这种方法适用于分离一些生长缓慢或数量较少的真菌，但需要对目标真菌的生长特性有一定的了解，以便设计合适的富集培养条件。目前，从竹黄中已分离出多种真菌，包括细菌、酵母菌和霉菌等。其中，一些真菌与竹黄的生长发育密切相关。有些真菌可能与竹黄形成共生关系，它们能够帮助竹黄吸收营养物质，增强竹黄的抗逆性，促进竹黄的生长和发育；而有些真菌则可能是竹黄的病原菌，会导致竹黄发生病害，影响竹黄的品质和产量。对这些已分离真菌的特性研究，包括它们的形态特征、生理生化特性、遗传特性等，有助于深入了解竹黄与真菌之间的相互作用关系。通过对真菌的形态观察，可以初步判断其分类地位；通过生理生化实验，可以了解真菌的代谢类型和对营养物质的需求；通过遗传分析，可以揭示真菌的遗传多样性和进化关系，为进一步研究竹黄的生物学特性和生态功能提供重要参考。三、遗传多样性分析方法3.1分子标记技术原理与应用3.1.1ISSR标记ISSR（Inter-SimpleSequenceRepeat）标记，即简单序列重复区间标记，是1994年由Zietkeiwitcz等发展起来的一种基于微卫星的分子标记技术。其原理是利用在真核生物基因组中广泛存在的SSR（简单序列重复），在SSR的5’或3’末端锚定2-4个随机核苷酸作为引物。在PCR反应中，锚定引物能与互补的SSR末端退火，从而引发位于反向排列、间隔不太大的重复序列间的DNA片段进行扩增。所扩增的interSSR区域的多个条带可通过聚丙烯酰胺凝胶电泳或琼脂糖凝胶电泳进行分辨，扩增谱带多表现为显性。例如，以（AG）₈T为引物，其中（AG）₈是SSR重复序列，T为锚定碱基，它可以与基因组中具有（AG）n反向排列且间隔合适的DNA区域结合并扩增。在竹黄遗传多样性分析中，ISSR标记具有诸多优势。它比RAPD（随机扩增多态性DNA）、SSR多态性更为丰富，能检测到更多的遗传变异位点。ISSR引物开发相对简单，无需像SSR引物那样进行繁琐的测序以获得SSR两侧的单拷贝序列，开发成本较低，且引物可以在不同物种间通用，不像SSR标记具有较强的种特异性。ISSR标记的重复性好、稳定性高，实验操作相对简便，成本也较为适中。在对江苏、安徽和浙江3省的8个居群共107个竹黄样本进行遗传多样性分析时，通过ISSR分子标记技术，筛选得到的11条引物共检测出241个位点，其中多态性位点240个，占99.59%，居群水平的多态位点百分率平均为64.21%，充分展示了其在揭示竹黄遗传多样性方面的强大能力。ISSR标记的实验操作流程如下：首先，提取竹黄样本的基因组DNA，确保DNA的质量和浓度符合要求，可采用CTAB法等经典方法进行提取。然后，根据加拿大BritishColumbia大学设计的ISSR引物序列购买成品引物或自行合成。在25μl反应体系中，依次加入模板DNA1μl（30-50ng）、ISSR引物1μl（约5pmol）、10xPCRBuffer2.5μl、MgCl₂2μl、dNTP2μl、Taq酶1单位（U），最后加ddH₂O至25μl，充分混匀后稍离心，若PCR仪无热盖，需加一滴（约20μl）矿物油。接着，在PCR仪中进行扩增反应，一般预变性94℃2分钟，然后进入循环：94℃1分钟，45℃-68℃40秒（具体退火温度需根据引物的Tm值优化确定），72℃1-2分钟，共40轮循环，循环结束后，72℃延伸10分钟，4℃保存。PCR产物经1.6%-1.8%琼脂糖凝胶电泳分离，稳压50-100V，电泳结束后用溴化乙锭染色20分钟，最后用凝胶成像仪观察、拍照记录条带。在实验操作过程中，有以下注意事项：要保证DNA模板的质量，其质量和数量会直接影响PCR效果，提取的DNA应避免降解和杂质污染；不同引物的退火温度需根据引物的Tm值进行适当调整，ISSR引物长度一般在15-24bp，反应退火温度通常在52-55℃，比RAPD高；当ISSR-PCR扩增的带型背景较强时，可在PCR反应体系中加入2%甲酰胺降低由于引物滑动引起的背景模糊，或加入2%-4%DMSO提高反应的特异性。3.1.2ITS序列分析ITS（InternalTranscribedSpacer）序列分析是基于真菌核糖体RNA（rRNA）基因内转录间隔区的DNA序列进行分析的技术。真菌核糖体基因由小的亚单元（18S）、ITS1区、5.8S区、ITS2区和大的亚单元（28S）构成，头尾串联形成重复序列，一个基因组内有60-200个拷贝。其中，ITS1和ITS2是中度保守区域，位于18S和5.8SrRNA基因之间的是ITS1，位于5.8S和28SrRNA基因之间的是ITS2，它们常被合称为ITS，并且5.8SRNA基因也被包括在ITS之内。rDNA上的5.8、18和28SrRNA基因具有极大的保守性，存在广泛的异种同源性，而ITS区不加入成熟核糖体，在进化过程中承受的自然选择压力非常小，能容忍更多的变异，在绝大多数真核生物中表现出极为广泛的序列多态性，即使是亲缘关系非常接近的两个种也能在ITS序列上表现出差异，其进化速率是18SrDNA的10倍，这使得ITS序列适合用于真菌物种的分子鉴定以及属内物种间或种内差异较明显的菌群间的系统发育关系分析。在竹黄分离真菌的研究中，ITS序列分析可用于准确鉴定分离真菌的种类。通过对从江苏宜兴、安徽广德以及浙江的安吉和丽水采集的竹黄子实体上分离纯化得到的36株真菌进行ITS-rDNA序列的系统学分析，利用邻接法和最大简约法分别对这36株菌的ITS序列和其他15条来自GenBank的相关序列进行分析，能够明确这些真菌的分类地位。研究结果表明，这些真菌中29株属于子囊菌，3株属于担子菌，4株属于半知菌，其中子囊菌中又以格孢腔菌目和梨孢假壳科成员为优势菌群，其次为肉座菌目和粪壳菌目，为深入了解竹黄与这些真菌的关系提供了重要依据。ITS序列分析的序列测定与分析方法如下：首先提取竹黄分离真菌的基因组DNA，可采用适合真菌的DNA提取试剂盒进行操作。然后，利用真菌ITS区域的通用引物对ITS区域进行PCR扩增，常用的引物如ITS1（5’-TCCGTAGGTGAACCTGCGG-3’）和ITS4（5’-TCCTCCGCTTATTGATATGC-3’）。PCR反应体系和条件可根据实际情况进行优化，一般反应体系包含模板DNA、引物、dNTP、PCR缓冲液、Taq酶等，反应条件包括预变性、变性、退火、延伸等步骤，退火温度通常在50℃-60℃之间。PCR扩增产物经琼脂糖凝胶电泳检测后，将目的条带切胶回收，进行测序，可选择专业的测序公司进行。测序得到的ITS序列利用生物信息学软件进行分析，首先将序列在NCBI（美国国立生物技术信息中心）的GenBank数据库中进行BLAST比对，查找与之相似的已知序列，初步确定真菌的分类地位。再使用MEGA（MolecularEvolutionaryGeneticsAnalysis）等软件构建系统发育树，常用的建树方法有邻接法、最大简约法等，通过系统发育树可以直观地展示竹黄分离真菌与其他相关真菌之间的亲缘关系。3.1.3其他分子标记技术AFLP（AmplifiedFragmentLengthPolymorphism）即扩增片段长度多态性，是RFLP（限制性片段长度多态性）与PCR两种技术的优点相结合的产物。其原理是先利用限制性内切酶水解基因组DNA产生不同大小的DNA片段，再使双链人工接头与酶切片段相连接，作为扩增反应的模板DNA，然后以人工接头的互补链为引物进行预扩增，最后在接头互补链的基础上添加1-3个选择性核苷酸作引物对模板DNA基因再进行选择性扩增，扩增的酶切片段在高分辨率的顺序分析胶上进行电泳，多态性即以扩增片段的长度不同被检测出来。AFLP具有分辨率高、稳定性好、效率高的优点，分析所需DNA量少，仅需0.5μg，可重复性好，多态性强，不需要Southern杂交，样品适用性广，但技术费用昂贵，对DNA的纯度和内切酶的质量要求很高。SRAP（Sequence-relatedAmplifiedPolymorphism）即相关序列扩增多态性，其原理是基于基因外显子中GC含量丰富，而启动子和内含子中AT含量丰富的特点设计引物。正向引物针对外显子区域设计，反向引物针对内含子区域和启动子区域设计，通过对开放阅读框（ORF）进行扩增来检测多态性。SRAP具有简便、稳定、产率高、成本低等优点，引物具有通用性，可在不同物种间使用，扩增的片段与基因相关，有利于基因的克隆和功能分析，但该技术对引物设计要求较高，需要根据不同物种进行优化。虽然AFLP和SRAP等分子标记技术在竹黄及其分离真菌遗传多样性分析中应用相对较少，但它们各自独特的优势使其在某些研究中具有潜在的应用价值。在需要高分辨率和稳定性的研究中，AFLP可作为一种有效的分析手段；而SRAP由于其与基因的相关性，对于研究竹黄及其分离真菌的基因功能和进化关系可能具有重要意义，为竹黄及其分离真菌遗传多样性研究提供了更多的技术选择。3.2数据分析方法在竹黄及其分离真菌的遗传多样性分析中，常用的数据统计指标包括多态位点百分率、基因多样性指数等。多态位点百分率（PPB）是指多态性位点的数量占总检测位点数量的百分比，它直观地反映了群体中遗传变异位点的丰富程度。公式为：PPB=（多态性位点数量/总检测位点数量）×100%。在对竹黄的研究中，若筛选得到的11条引物共检测出241个位点，其中多态性位点240个，则多态位点百分率为（240/241）×100%≈99.59%，表明竹黄群体具有丰富的遗传变异位点。基因多样性指数（Nei’sgenediversityindex），又称Nei’s指数，用于衡量群体内基因的多样性程度。其计算公式为：H=1-\sum_{i=1}^{n}p_{i}^{2}，其中p_{i}是第i个等位基因的频率，n是等位基因的数量。Nei’s指数的值介于0-1之间，值越大表示基因多样性越高。例如，当一个群体中所有个体的基因完全相同，即只有一种等位基因时，p_{i}=1，H=0，表示基因多样性最低；而当群体中存在多种等位基因且频率相等时，H的值接近1，基因多样性较高。Shannon多态性信息指数（Shannon’sinformationindex）也是常用的遗传多样性指标，它综合考虑了群体中基因的丰富度和均匀度。公式为：I=-\sum_{i=1}^{n}p_{i}ln(p_{i})，其中p_{i}和n的含义与Nei’s指数公式中相同。I值越大，说明群体的遗传多样性越高，该指数不仅反映了基因的多态性，还考虑了不同等位基因在群体中的分布情况。聚类分析是揭示竹黄及其分离真菌遗传关系的重要数据处理方法之一，常用的聚类方法有非加权组平均法（UPGMA，UnweightedPair-GroupMethodwithArithmeticMean）。UPGMA根据遗传距离矩阵，将遗传距离最近的两个类群合并为一个新类群，然后重新计算新类群与其他类群之间的遗传距离，再将距离最近的类群合并，如此反复，直至所有类群合并为一个大类群，最终形成一个聚类树状图。在竹黄地理居群的研究中，通过UPGMA聚类分析发现，地理上距离较近的居群多聚在一起，表现出一定的地域性，说明各居群的亲缘关系与地理分布有着密切关系。主坐标分析（PCA，PrincipalCoordinatesAnalysis）是一种基于遗传距离矩阵的多元统计分析方法，它将高维的数据通过线性变换投影到低维空间中，使得数据在低维空间中的分布能够最大程度地反映原始数据的变异信息。PCA可以将竹黄及其分离真菌的遗传数据转化为几个主坐标，通过绘制主坐标图，直观地展示不同样本或居群之间的遗传关系。在主坐标图中，距离较近的点表示遗传关系较近，距离较远的点表示遗传关系较远，从而可以清晰地看出竹黄及其分离真菌的遗传结构和分化情况，有助于深入理解它们的遗传多样性和进化历程。四、竹黄遗传多样性分析实例4.1实验材料与方法本研究的竹黄样本分别采自江苏、安徽和浙江3省的8个自然居群，共收集了107个样本。在江苏宜兴采集了20个样本，该地位于太湖西岸，属北亚热带季风气候，竹林资源丰富，竹黄多寄生于短穗竹上；安徽广德采集18个样本，广德地处皖东南，气候温和湿润，其境内的竹林为竹黄提供了适宜的生长环境，主要寄主为箭竹；浙江安吉采集20个样本，安吉素有“中国竹乡”之称，生态环境优良，竹黄常见于毛环短穗竹上；浙江天目山采集15个样本，天目山森林覆盖率高，生态系统完整，为竹黄的生长创造了独特的条件，寄主竹子种类多样。此外，还在浙江丽水、临安等地采集了剩余样本，每个居群的样本采集点尽量覆盖该地区有竹黄分布的不同区域，以保证样本的代表性。采集时，选择生长良好、无明显病虫害的竹黄子实体，用无菌剪刀从竹子上剪下，放入无菌塑料袋中，做好标记，记录采集地点、时间、寄主等信息，迅速带回实验室，置于4℃冰箱保存备用。DNA提取采用CTAB法。具体步骤如下：取约0.1g竹黄组织，置于液氮中充分研磨成粉末状，迅速转移至含有700μl预热至65℃的2×CTAB提取缓冲液（含2%CTAB、100mmol/LTris-HCl，pH8.0、20mmol/LEDTA，pH8.0、1.4mol/LNaCl、0.2%β-巯基乙醇）的1.5ml离心管中，轻轻颠倒混匀，65℃水浴保温30-60分钟，期间每隔10分钟轻轻颠倒混匀一次。水浴结束后，冷却至室温，加入等体积的氯仿：异戊醇（24:1），轻轻颠倒混匀10-15分钟，使溶液充分乳化，12000r/min离心10分钟。将上清液转移至新的1.5ml离心管中，加入2/3体积预冷的异丙醇，轻轻颠倒混匀，-20℃放置30分钟，使DNA沉淀。12000r/min离心10分钟，弃上清液，用70%乙醇洗涤沉淀2-3次，每次洗涤后12000r/min离心5分钟，弃上清液，室温晾干或真空干燥DNA沉淀。加入50-100μlTE缓冲液（含10mmol/LTris-HCl，pH8.0、1mmol/LEDTA，pH8.0）溶解DNA，置于4℃冰箱保存备用，使用前通过紫外分光光度计检测DNA的浓度和纯度，要求OD₂₆₀/OD₂₈₀比值在1.8-2.0之间，OD₂₆₀/OD₂₃₀比值大于2.0。ISSR-PCR扩增在25μl反应体系中进行，各成分的优化过程如下：通过单因子试验，分别设置不同的Mg²⁺浓度（1.5mmol/L、2.0mmol/L、2.5mmol/L、3.0mmol/L、3.5mmol/L）、dNTPs浓度（0.1mmol/L、0.15mmol/L、0.2mmol/L、0.25mmol/L、0.3mmol/L）、TaqDNA聚合酶用量（0.5U、1.0U、1.5U、2.0U、2.5U）、引物浓度（0.2μmol/L、0.3μmol/L、0.4μmol/L、0.5μmol/L、0.6μmol/L）和模板DNA用量（20ng、30ng、40ng、50ng、60ng），进行PCR扩增，观察扩增条带的清晰度、亮度和特异性，初步确定各成分的适宜浓度范围。在此基础上，采用L₁₆（4⁵）正交设计方法，对这5个因素进行优化，每个因素设置4个水平，共进行16组实验，通过分析正交实验结果，确定最佳反应体系为：模板DNA30ng、ISSR引物0.4μmol/L、10×PCRBuffer2.5μl、MgCl₂2.5mmol/L、dNTP0.2mmol/L、Taq酶1.0U，加ddH₂O至25μl。扩增程序为：预变性94℃5分钟；然后进行35个循环，每个循环包括94℃变性45秒，52℃-55℃退火45秒（根据不同引物的Tm值进行调整），72℃延伸1.5分钟；循环结束后，72℃延伸10分钟，4℃保存。电泳检测采用1.6%-1.8%的琼脂糖凝胶，电泳缓冲液为1×TAE。将PCR扩增产物与6×LoadingBuffer按5:1的比例混合，取5-10μl混合液加入凝胶的加样孔中，同时加入DNAMarker作为分子量标准。稳压50-100V进行电泳，电泳时间根据凝胶长度和电压适当调整，一般为1-2小时，使DNA条带充分分离。电泳结束后，将凝胶置于含有0.5μg/ml溴化乙锭的溶液中染色20分钟，然后用凝胶成像仪观察、拍照记录条带，条带清晰、重复性好的样本用于后续分析。4.2实验结果与分析经过ISSR扩增，筛选得到的11条引物共检测出241个位点，其中多态性位点240个，多态性位点百分率（PPB）高达99.59%，这表明竹黄在物种水平上具有极为丰富的遗传变异。在居群水平上，多态位点百分率平均为64.21%，其中浙江天目山居群的多态位点百分率最高，达到74.27%，这可能与天目山独特的生态环境和丰富的竹子寄主资源有关，多样化的生态条件和寄主种类为竹黄的遗传变异提供了更多的机会；江苏宜兴居群的多态位点百分率最低，为45.23%，可能是由于该地区的竹黄分布相对集中，生态环境较为单一，遗传交流相对较少，导致遗传变异相对匮乏。通过数据分析，Nei’s基因多样性指数为0.3314，Shannon多态性信息指数为0.4996，进一步表明竹黄具有较高的遗传多样性水平。较高的基因多样性指数意味着竹黄群体中存在多种不同的等位基因，基因的丰富度较高；而Shannon多态性信息指数不仅考虑了基因的丰富度，还考虑了基因在群体中的分布均匀度，该指数较高说明竹黄群体中基因的分布相对较为均匀，没有出现某一种或几种基因占据绝对优势的情况。利用Arlequin软件进行AMOVA分析，结果显示居群间变异占39.96%，居群内变异占60.04%，这表明竹黄的遗传变异主要存在于居群内。这可能是因为竹黄的繁殖方式以无性繁殖为主，无性繁殖使得居群内个体之间的遗传差异相对较小，但由于长期的进化和适应不同的微环境，居群内仍然积累了一定的遗传变异；而居群间由于地理隔离、寄主差异等因素，虽然存在一定的遗传分化，但相对居群内变异而言，居群间变异所占比例较小。8个居群间的遗传距离变化范围在0.1115-0.2453之间，平均为0.1766，这表明不同居群的竹黄之间存在一定的遗传差异。遗传距离反映了不同居群间遗传物质的差异程度，遗传距离越大，说明两个居群之间的遗传差异越明显。居群间的基因流Nm=1.0086，说明居群间存在一定的基因交流，但基因交流相对较少。基因流是指由于个体的迁移或花粉、种子的传播等导致的基因在不同居群间的流动，基因流可以促进居群间的遗传交流，减少遗传分化。采用UPGMA聚类分析和主坐标分析（PCA）对竹黄居群的遗传关系进行研究。UPGMA聚类分析结果显示，地理上距离较近的居群多聚在一起，表现出一定的地域性，说明各居群的亲缘关系与地理分布有着密切关系。例如，江苏宜兴和浙江安吉的居群在聚类图中距离较近，它们在地理位置上相邻，生态环境也较为相似，可能存在一定程度的基因交流，从而导致亲缘关系较近；而贵州地区来源的菌株标本与其他地区来源的菌株标本亲缘关系较远，可能是由于贵州与其他地区地理距离较远，地理隔离限制了基因交流，使得贵州地区的竹黄在遗传上逐渐分化，形成了独特的遗传特征。PCA分析结果也进一步验证了这一结论，在主坐标图中，不同居群的竹黄样本按照地理分布呈现出一定的聚集趋势，进一步表明地理因素对竹黄的遗传结构具有重要影响。4.3影响竹黄遗传多样性的因素地理分布是影响竹黄遗传多样性的重要因素之一。不同地理区域的竹黄居群由于受到地理隔离的影响，基因交流受到限制，逐渐积累了不同的遗传变异，从而导致遗传分化。如前文所述，贵州地区的竹黄菌株标本与其他地区的亲缘关系较远，这是因为贵州与江苏、浙江等地地理距离较远，竹黄的传播主要依赖于风力、动物等自然因素，长距离的地理隔离使得贵州地区的竹黄难以与其他地区的竹黄进行有效的基因交流，在长期的进化过程中，形成了独特的遗传特征。不同的地理环境条件，包括气候、土壤、海拔等，也会对竹黄的遗传多样性产生影响。在气候温暖湿润、土壤肥沃、海拔适中的地区，竹黄的生长环境较为适宜，有利于其繁殖和生存，可能会出现更多的遗传变异；而在环境条件较为恶劣的地区，竹黄的生长受到限制，遗传多样性可能相对较低。海拔高度的变化会导致温度、光照、湿度等环境因素的改变，这些因素可能会影响竹黄的生理代谢和基因表达，从而影响其遗传多样性。在高海拔地区，温度较低，光照较强，竹黄可能会通过调整自身的基因表达来适应这种环境，导致遗传变异的积累。寄主种类对竹黄的遗传多样性也有显著影响。不同寄主竹子的生理特性、营养成分和代谢产物存在差异，这些差异会影响竹黄在寄主体内的生长、繁殖和遗传变异。在营养成分丰富、代谢活性较高的寄主竹子上，竹黄可能会获得更多的营养物质，生长更加旺盛，繁殖能力增强，从而增加遗传变异的机会；而在营养成分匮乏、生长状态不佳的寄主竹子上，竹黄的生长可能会受到抑制，遗传多样性也可能受到影响。寄主竹子的遗传背景也可能与竹黄的遗传特性相互作用，影响竹黄的遗传多样性。不同寄主竹子的基因表达差异，可能会导致竹黄在寄生过程中发生基因水平的变化，从而影响竹黄的遗传特征和生态适应性。环境因素，如温度、湿度、光照等，对竹黄的遗传变异有着重要的作用。温度是影响竹黄生长和代谢的关键环境因素之一，适宜的温度有利于竹黄的生长和繁殖，而极端温度可能会导致竹黄的生理功能紊乱，影响其遗传稳定性。在高温环境下，竹黄的酶活性可能会受到影响，导致代谢异常，进而影响基因的表达和遗传变异；在低温环境下，竹黄的生长速度可能会减缓，繁殖能力下降，也可能会影响遗传多样性。湿度对竹黄的影响也不容忽视，过高或过低的湿度都可能对竹黄的生长和发育产生不利影响。在高湿度环境下，竹黄容易受到病虫害的侵袭，可能会导致遗传变异的发生；而在低湿度环境下，竹黄可能会面临水分胁迫，影响其生理代谢和遗传稳定性。光照条件同样会影响竹黄的遗传变异，不同的光照强度和光周期可能会影响竹黄的光合作用和激素水平，进而影响其基因表达和遗传多样性。人为因素对竹黄遗传多样性的破坏作用不容忽视。过度采集是导致竹黄遗传多样性下降的主要人为因素之一。由于竹黄具有较高的药用价值，市场需求不断增加，人们对竹黄的采集强度也日益加大。过度采集使得竹黄的种群数量急剧减少，许多居群面临灭绝的危险，这不仅直接导致竹黄遗传资源的丧失，还会破坏其生态环境，影响其繁殖和生存，进一步降低遗传多样性。不合理的采集方式，如连根拔起、大面积砍伐竹林等，也会对竹黄的遗传多样性造成严重破坏。这些采集方式不仅会直接损害竹黄的子实体，还会破坏其寄主竹子的生长环境，影响竹黄与寄主之间的共生关系，导致竹黄的遗传多样性降低。此外，人类活动导致的生态环境破坏，如森林砍伐、环境污染等，也会对竹黄的遗传多样性产生负面影响。森林砍伐使得竹黄的栖息地减少，寄主竹子的数量下降，竹黄的生存空间受到挤压；环境污染，如土壤污染、水污染和空气污染等，可能会影响竹黄的生长和发育，导致其遗传变异异常，甚至引发基因突变，从而破坏竹黄的遗传多样性。五、竹黄分离真菌遗传多样性分析实例5.1分离真菌的获取与鉴定从竹黄子实体分离真菌时，选用马铃薯葡萄糖琼脂（PDA）培养基，该培养基富含多种营养成分，能够满足大多数真菌的生长需求。其配方为：马铃薯200g，葡萄糖20g，琼脂15-20g，水1000ml。先将马铃薯去皮，切成小块，煮沸30分钟后过滤，取滤液加入葡萄糖和琼脂，加热搅拌使其完全溶解，调节pH值至自然状态，分装到三角瓶中，用棉塞塞紧瓶口，包扎后于121℃高压蒸汽灭菌20-30分钟，待冷却至50℃-60℃时，在无菌条件下倒入无菌培养皿中，制成平板备用。采集江苏宜兴、安徽广德以及浙江的安吉和丽水等地的竹黄子实体，将采集到的子实体带回实验室后，先用无菌水冲洗表面，去除杂质，然后用75%乙醇消毒30-60秒，再用无菌水冲洗3-5次，以确保表面消毒彻底。在无菌操作台上，用无菌镊子和剪刀将竹黄子实体切成约0.5cm×0.5cm的小块，将小块接种到PDA平板上，每个平板接种3-5块，置于25℃恒温培养箱中培养，定期观察真菌的生长情况。待平板上长出真菌菌落时，进行形态学观察。不同真菌的菌落形态具有一定的特征，如颜色、形状、质地、边缘等。有些真菌的菌落呈白色，绒毛状，边缘整齐；有些则呈黑色，颗粒状，边缘不规则。通过仔细观察这些形态特征，可以初步判断真菌的类别。对菌丝和孢子的形态进行显微镜观察，进一步辅助分类。使用无菌牙签挑取少量菌丝或孢子，置于载玻片上，滴加一滴蒸馏水，盖上盖玻片，在显微镜下观察其形态、大小、结构等特征，如菌丝是否有分隔、孢子的形状和颜色等。为了更准确地鉴定分离真菌的种类，进行ITS序列分析。采用真菌DNA提取试剂盒提取分离真菌的基因组DNA，严格按照试剂盒说明书的步骤进行操作，确保提取的DNA纯度和浓度符合要求。利用通用引物ITS1（5’-TCCGTAGGTGAACCTGCGG-3’）和ITS4（5’-TCCTCCGCTTATTGATATGC-3’）对ITS区域进行PCR扩增。在25μl反应体系中，加入模板DNA1μl（约50ng）、引物ITS1和ITS4各1μl（10μmol/L）、10×PCRBuffer2.5μl、dNTP2μl（2.5mmol/L）、Taq酶1U，加ddH₂O至25μl。PCR扩增程序为：94℃预变性5分钟；然后94℃变性30秒，55℃退火30秒，72℃延伸1分钟，共35个循环；最后72℃延伸10分钟。PCR扩增产物经1.5%琼脂糖凝胶电泳检测，在紫外凝胶成像系统下观察并拍照，确认扩增条带的大小和特异性。将特异性扩增条带切胶回收，送测序公司进行测序。测序得到的ITS序列在NCBI的GenBank数据库中进行BLAST比对，查找与之相似的已知序列，根据比对结果初步确定分离真菌的分类地位。利用MEGA软件，采用邻接法（NJ）构建系统发育树，分析分离真菌与其他相关真菌之间的亲缘关系，进一步明确其分类地位。5.2遗传多样性分析结果基于ITS序列分析，对从竹黄子实体分离得到的36株真菌进行遗传多样性分析。通过序列比对和数据分析，发现这些真菌具有丰富的遗传多样性，在系统发育树上呈现出多个不同的分支，表明它们属于不同的真菌类群。系统发育树分析结果显示，这些真菌中29株属于子囊菌，3株属于担子菌，4株属于半知菌。在子囊菌中，格孢腔菌目（Pleosporales）和梨孢假壳科（Apiosporaceae）的成员构成了优势菌群。这可能是因为竹黄的生长环境为这些真菌提供了适宜的生存条件，它们与竹黄之间可能存在着某种共生或协同进化关系。这些优势菌群在竹黄的生态系统中可能扮演着重要角色，如参与竹黄的物质代谢、影响竹黄的生长发育等。半知菌和担子菌虽然数量相对较少，但它们在竹黄的生态系统中也具有不可忽视的作用。半知菌在自然界中广泛存在，它们可能参与竹黄的分解和转化过程，对竹黄的物质循环起到一定的促进作用；担子菌则可能与竹黄在营养获取、生态位竞争等方面存在相互作用，它们的存在丰富了竹黄生态系统的生物多样性。不同类群的真菌在竹黄子实体中的分布存在一定差异。在江苏宜兴采集的竹黄子实体中，分离得到的真菌以子囊菌中的格孢腔菌目成员为主；而在浙江安吉的样本中，梨孢假壳科的真菌相对较多。这种分布差异可能与不同地区的地理环境、气候条件以及竹子寄主的差异有关。地理环境的差异会影响土壤的酸碱度、养分含量等因素，从而影响真菌的生长和分布；气候条件的不同，如温度、湿度等，也会对真菌的生存和繁殖产生影响；竹子寄主的差异，包括竹子的品种、生理状态等，可能会导致竹黄子实体的化学成分和微生态环境发生变化，进而影响真菌的种类和分布。5.3分离真菌与竹黄的关系探讨从生态角度来看，分离真菌与竹黄之间可能存在共生或寄生关系。共生关系下，真菌与竹黄相互依存、互利共赢。某些真菌可能能够帮助竹黄吸收土壤中的养分，如磷、钾等，提高竹黄对营养物质的利用效率；真菌还可能分泌一些生长调节物质，促进竹黄的生长和发育，增强竹黄的抗逆性，使其能够更好地适应环境变化。寄生关系中，真菌则可能从竹黄中获取营养，对竹黄的生长产生一定的负面影响。一些病原菌可能会侵染竹黄，导致竹黄出现病害症状，如子座变色、腐烂等，影响竹黄的品质和产量。在遗传角度，分离真菌可能对竹黄的遗传多样性产生影响。真菌与竹黄在长期的相互作用过程中，可能会发生基因交流。真菌的某些基因可能会转移到竹黄基因组中，或者竹黄的基因可能会影响真菌的基因表达，从而导致两者的遗传物质发生改变，增加遗传多样性。这种基因交流可能会使竹黄获得一些新的性状，如对某些病原菌的抗性增强、对环境的适应性提高等，有利于竹黄在不同的生态环境中生存和繁衍。在进化过程中，竹黄与分离真菌之间存在着相互作用。它们可能在长期的进化过程中共同适应环境变化，逐渐形成了稳定的生态关系。随着环境的改变，竹黄和真菌可能会发生适应性进化，通过基因突变、基因重组等方式，调整自身的遗传特性，以更好地适应新的环境条件。这种相互作用推动了竹黄和分离真菌的共同进化，使得它们在生态系统中保持着相对稳定的状态，共同维持着生态平衡。六、结论与展望6.1研究主要结论本研究通过对竹黄及其分离真菌的遗传多样性分析，取得了一系列重要成果。在竹黄遗传多样性方面，利用ISSR分子标记技术对江苏、安徽和浙江3省8个居群的107个竹黄样本进行分析，结果显示竹黄在物种水平上具有极高的遗传多样性，多态位点百分率达99.59%。居群水平上，多态位点百分率平均为64.21%，其中浙江天目山居群最高，江苏宜兴居群最低。Nei’s基因多样性指数为0.3314，Shannon多态性信息指数为0.4996，进一步表明竹黄遗传多样性丰富。遗传结构分析表明，竹黄居群间变异占39.96%，居群内变异占60.04%，遗传变异主要存在于居群内。8个居群间遗传距离变化范围在0.1115-0.2453之间，平均为0.1766，基因流Nm=1.0086，说明居群间存在一定基因交流。UPGMA聚类分析和主坐标分析显示，地理距离相近的居群多聚在一起，表明各居群亲缘关系与地理分布密切相关。在竹黄分离真菌遗传多样性分析中，从江苏宜兴、安徽广德以及浙江的安吉和丽水等地采集的竹黄子实体上分离纯化得到36株真菌，经ITS-rDNA序列系统学分析，这些真菌中29株属于子囊菌，3株属于担子菌，4株属于半知菌。子囊菌中格孢腔菌目和梨孢假壳科成员为优势菌群，其次为肉座菌目和粪壳菌目，不同类群真菌在竹黄子实体中的分布存在差异。6.2研究的创新点与不足本研究的创新点主要体现在研究方法和研究角度上。在方法应用方面，运用ISSR分子标记技术对竹黄地理居群进行遗传多样性分析，该技术多态性丰富、稳定性高，能有效揭示竹黄的遗传变异情况；同时，采用ITS序列分析对竹黄分离真菌进行系统学鉴定，为准确识别与竹黄相关的真菌种类提供了可靠手段。在研究角度上，不仅关注竹黄自身的遗传多样性，还深入探究了竹黄分离真菌的遗传多样性及其与竹黄的关系，为全面了解竹黄生态系统提供了新的视角。然而，本研究也存在一些不足之处。样本数量有限，仅采集了江苏、安徽和浙江3省的竹黄样本，未能涵盖竹黄在我国的所有分布区域，可能无法全面反映竹黄的遗传多样性；对于竹黄分离真菌，虽然进行了ITS序列分析，但未能进一步深入研究其生理生化特性和功能，对它们与竹黄相互作用的分子机制了解有限；在遗传多样性分析中，仅采用了ISSR和ITS两种分子标记技术，分析方法相对单一，可能会遗漏一些遗传信息。6.3未来研究方向未来在竹黄及其分离真菌遗传多样性研究领域，可从以下几个方面展开工作。扩大样本范围，广泛采集我国其他地区以及国外竹黄分布区域的样本，增加样本数量，全面揭示竹黄的遗传多样性及其地理分布格局，深入研究不同地理区域竹黄的遗传差异和进化关系。采用多种分子标记技术联合分析，如将AFLP、SRAP等技术与ISSR、ITS相结合，从多个角度分析竹黄及其分离真菌的遗传多样性，提高分析结果的准确性和全面性。深入研究遗传多样性与环境适应性的关系，分析环境因素对竹黄及其分离真菌遗传变异的影响，揭示它们在不同环境条件下的适应性进化机制，为竹黄的资源保护和人工栽培提供科学依据。加强对竹黄分离真菌生理生化特性和功能的研究，探索它们与竹黄之间的共生、寄生等关系，明确它们在竹黄生长发育和生态系统中的作用，为竹黄的可持续开发利用提供理论支持。七、参考文献[1]赵宁，陈双林。竹黄遗传多样性与竹红菌素合成和抗癌作用的研究进展[J].生物技术通报，2018,34(04):1-6.[2]吴疆，侯文彬，高欣。药用真菌竹黄的研究进展[J].中草药，2009,40(S1):91-92.[3]胡晓，沈联德。竹黄化学成分的分离和结构鉴定[J].华西药学杂志，1992,(01):1-4.[4]范英。竹黄及其分离真菌的遗传多样性分析[D].南京师范大学，2011.[5]李向敏，陈双林，陈辉，李泰辉。广西产竹黄菌的分子鉴定[J].广西农业科学，2009,40(09):1113-1116.[6]张红，王芳，丁莉，赵天增。竹黄研究进展[J].河南科学，2006,(03):361-364.[7]李加友，李兆兰，焦庆才，张立新。竹黄菌发酵液萃取物的抗菌活性研究[J].南京中医药大学学报，2003,(03):159-160.[8]沈云修，荣先国，高宗华。竹黄的化学成分研究[J].中国中药杂志，2002,(09):674-676.[9]林海萍，陈虹，叶勇，吴林森。竹黄竹红菌甲素含量测定方法[J].浙江林学院学报，2002,(02):157-160.[10]陈佳佳，扶教龙，胡翠英，李兆兰，焦庆才。竹黄多糖的提取工艺及清除自由基活性研究[J].江苏农业科学，2008,(02):196-198.[11]刘雅峰，潘勤。真菌竹黄中的过氧麦角甾醇的分离[J].天津中医学院学报，2004,(01):15-16.[12]蒋丽金。竹红菌素的结构、性质、光化学反应及反应机制(Ⅰ)——竹红菌素的结构和性质[J].科学通报，1990,(21):1608-1616.[13]陈远腾，万象义。竹黄(ShiraiabambusicolaP.Henn.)主要光敏有效成分的初步探讨[J].云南大学学报(自然科学版),1981,(02):104-107.[14]钟琪，周雷罡，袁铭铭，廖圆月，周国平.HPLC法同时测定真菌竹黄胶囊中3个醌类成分的含量[J].药物分析杂志，2020,40(01):139-144.[15]钟琪，周雷罡，袁铭铭，廖圆月，周国平。真菌竹黄化学成分及药理活性研究进展[J].世界最新医学信息文摘，2020,20(44):42-43.[16]李加友。竹黄菌的液体发酵研究[J].食药用菌，2008,16(05):24-26.[17]蔡宇杰，丁彦蕊，张大兵，楼志华，石贵阳。竹黄菌固态发酵竹红菌素条件的研究[J].生物技术，2004,(04):46-47.[18]夏盛钦，陈景荣，王雪松，张宝文。一种新型葡萄糖基化竹红菌素的合成和ESR研究[J].感光科学与光化学，2004,(05):363-369.[19]石贵阳，张大兵，楼志华，蔡宇杰。竹黄菌液体培养条件下生成竹红菌素的研究[J].药物生物技术，2004,(05):299-301.[20]乐加昌，李美芬，庞素珍。激发态色氨酸与竹红菌乙素之间光诱导电子转移的ESR研究[J].化学通报，1993,(02):41-42.[21]王大辉，张志义，张过江，赫荣乔。竹红菌甲素激发态结构——激基缔合物和激基复合物[J].科学通报，1993,(11):1032-1036.[