探秘银杏变异斑叶呈色：结构、生理与分子调控的多维解析

探秘银杏变异斑叶呈色：结构、生理与分子调控的多维解析一、引言1.1研究背景与意义银杏（Ginkgobiloba）作为世界著名的孑遗树种，享有“活化石”的美誉，是现存种子植物中最古老的孑遗植物之一，在植物演化历史中占据着独特而关键的地位。其不仅见证了漫长岁月中的地质变迁与气候变化，为植物进化研究提供了珍贵的“活标本”，而且具有多方面重要的经济价值。在药用领域，银杏的果实（白果）和叶片富含黄酮类、萜内酯等多种生物活性成分，具有敛肺定喘、祛痰止带、改善心脑血管循环、抗氧化、抗炎等药用功效，被广泛应用于心脑血管疾病的预防与治疗以及保健产品的开发；在食用方面，白果含有丰富的营养物质，如脂肪、胡萝卜素、钙、铁等，经过适当加工处理后可供食用，为人们提供独特的营养补充；在材用方面，银杏木材材质优良，纹理直，结构细，不翘不裂，可用于制作家具、雕刻工艺品等，具有较高的工艺价值。此外，银杏叶形独特，呈扇形且具有二叉状平行叶脉，在秋季叶色会由绿转金黄，满树金黄的景观极具视觉冲击力，为城市、公园、景区等增添了独特的景致，使其成为备受青睐的园林景观树种，具有极高的观赏价值，深受人们喜爱。植物叶色变异现象在自然界中较为普遍，叶色突变体是植物在自然生长过程中或受到外界环境因素影响，导致遗传物质发生改变而产生的特殊类型，其叶色表现出与正常植株明显不同的特征。这类突变体是研究植物叶绿体发育、光合作用以及光形态建成等生理过程的理想材料。叶绿体作为植物进行光合作用的重要细胞器，其正常发育对植物的生长和生存至关重要。叶色突变体往往伴随着叶绿体结构和功能的异常，通过对这些突变体的研究，可以深入了解叶绿体发育的分子机制和调控网络，揭示光合作用过程中光能吸收、传递和转化的奥秘，以及光信号如何调控植物的形态建成和生长发育。在过去，虽然对叶色突变体的研究取得了一定进展，但大多集中在突变体的发生频率、遗传规律以及突变基因的初步定位等方面，对于叶色变异呈色机理，尤其是在分子水平上的深入研究还相对较少，许多关键的调控机制和信号通路仍有待进一步探索和阐明。银杏在长期的自然生长和人工栽培过程中，由于立地条件的差异以及天然、人工选择的作用，出现了众多栽培变种，其中银杏变异斑叶便是一种独特的叶色突变类型。这些变异斑叶呈现出黄绿相间的条纹状或斑块状，其独特的叶色表型不仅丰富了银杏的观赏多样性，更为研究植物叶色变异的分子机制提供了难得的材料。研究银杏变异斑叶呈色机理在理论和应用方面均具有重要意义。在理论层面，通过深入探究银杏变异斑叶呈色的分子机制，可以进一步揭示植物叶色形成和调控的复杂过程，填补植物叶色变异研究领域在银杏这一古老树种上的空白，完善植物色素代谢、叶绿体发育以及基因表达调控等相关理论体系，为理解植物适应环境变化的分子机制提供新的视角和理论依据。在应用层面，明确银杏变异斑叶的呈色机理，有助于开发基于分子标记的银杏优良观赏品种选育技术，加快培育具有更加丰富叶色和更高观赏价值的银杏新品种，满足园林景观对多样化观赏植物的需求，推动银杏在园林景观中的广泛应用和推广；同时，对于利用基因工程技术改良其他植物的叶色，创造具有独特观赏性状的植物新材料也具有重要的借鉴意义，为园艺植物的遗传改良和品种创新提供新的思路和方法。1.2研究现状近年来，随着分子生物学技术的飞速发展，植物叶色变异的研究取得了显著进展，尤其是在模式植物拟南芥以及一些重要经济作物如水稻、玉米、小麦等方面。在拟南芥中，已鉴定出大量与叶色相关的突变基因，这些基因参与了叶绿素合成、叶绿体发育、光信号传导等多个关键过程。通过对这些突变基因的功能分析，揭示了许多叶色变异的分子调控机制，为植物叶色研究提供了重要的理论基础和研究范式。在水稻中，叶色突变体的研究不仅有助于解析水稻光合作用和叶绿体发育的分子机制，还为水稻遗传育种提供了重要的标记性状和种质资源。例如，一些叶色突变体由于光合效率的改变，表现出独特的生长特性和产量潜力，通过对其遗传机制的深入研究，可以为水稻高产优质品种的选育提供新的思路和方法。在银杏叶色变异研究方面，目前已有一些相关报道。扬州大学银杏研究团队对银杏金叶突变体进行了深入研究，利用细胞生物学、生理学、转录组学测序及分子生物学等多学科手段，系统地比较了突变体金叶与正常绿叶的差异。研究发现，突变体金叶中叶绿体数量少、体积较小，基粒片层结构松散、类囊体膜结构不清晰，质体小球数量明显增多。在色素含量方面，金叶中叶绿素合成代谢路径中关键产物发生变化，UrogenIII和CoprogenIII含量增加，而ProtoIX、Mg-Proto以及Pchlide含量下降，叶绿素b含量显著降低，类胡萝卜素中主要是叶黄素含量显著升高。通过高通量测序及RT-qPCR验证，进一步鉴定出叶绿素、类胡萝卜和类黄酮等色素合成代谢通路中的关键功能基因发生了改变，如PPO、CHLD、CHLH、CHLI的下调表达阻遏了镁原卟啉IX的合成，导致叶绿素生物合成效率降低；NYC/NOL上调表达致使叶绿素b加速降解为叶绿素a，是影响突变叶中叶绿素积累减少的主要原因；叶绿体发育分化过程中关键差异基因GbGLK、GbFtsz在突变叶中显著下调，解释了叶绿体发育不完整的原因，表明银杏突变体的呈色与叶绿体的发育密切相关；类胡萝卜素合成过程中的关键基因如ZDS、LCYE等均在不同程度上调表达，促进了金叶中类胡萝卜素的积累。研究最终鉴定出银杏金叶形成的关键功能基因，明确了叶绿素代谢通路中合成受阻和降解加速、以及类胡萝卜素中叶黄素的合成增加，导致类胡萝卜素与叶绿素比值的显著升高，是银杏叶色突变体金叶呈色的关键因素。另有研究以银杏斑叶变异植株的叶片和正常叶片为试验材料，通过叶色比较、色素含量测定、解剖结构观察以及高通量测序转录组分析，探究银杏斑叶叶色形成的变化规律及叶色变化的分子调控机制。结果表明，银杏斑叶中叶绿素色素含量相对较低，类胡萝卜素含量较高；超微结构显示斑叶的黄色组织中叶绿体发育不健全，基粒片层垛叠程度较低，质体小球积累较多。进一步的转录组测序分析发现，斑叶中编码叶绿素合成的大部分基因表达下调，而部分类胡萝卜素和花青素合成的基因表达上调。综合上述研究结果，认为斑叶银杏的形成主要是由于叶绿素合成障碍以及类胡萝卜和花青素的合成增加，使得叶片内色素含量及比例发生改变，从而导致银杏斑叶叶色发生变异。尽管目前在银杏叶色变异研究方面取得了一定成果，但对于银杏变异斑叶呈色机理的研究仍存在诸多不足。现有研究主要集中在金叶突变体和斑叶突变体的整体分析上，对于变异斑叶中不同颜色区域（如黄色斑块与绿色区域）在细胞结构、生理生化以及分子水平上的差异对比研究还不够深入和系统。在分子机制方面，虽然已经鉴定出一些与叶色变异相关的基因，但这些基因之间的相互作用关系以及它们如何协同调控银杏变异斑叶的呈色过程尚未完全明确。此外，环境因素对银杏变异斑叶呈色的影响机制也有待进一步探究，如光照强度、温度、土壤养分等环境因子如何通过影响基因表达和生理生化过程来调控叶色变化，目前还缺乏深入的研究。因此，深入开展银杏变异斑叶呈色机理的研究，对于全面揭示银杏叶色变异的分子机制，丰富植物叶色变异理论体系，具有重要的科学意义。二、材料与方法2.1实验材料本研究选取的银杏变异斑叶植株与正常绿叶植株，均来源于[具体种植地点]的银杏种质资源圃。该种质资源圃位于[地理位置描述]，地理位置优越，光照充足，年平均光照时长达到[X]小时，能够充分满足银杏生长对光照的需求；气候温暖湿润，年平均气温在[X]℃左右，年平均降水量约为[X]毫米，为银杏的生长提供了适宜的气候条件。土壤类型主要为[土壤类型名称]，土层深厚肥沃，土壤pH值维持在[X]之间，呈微酸性至中性，且富含多种矿物质元素和有机质，能够为银杏生长提供丰富的养分。圃内的银杏植株均由专业人员按照科学的栽培管理技术进行日常养护，定期进行浇水、施肥、修剪、病虫害防治等工作，确保植株生长健壮，无明显病虫害侵扰。在2023年5月中旬，于银杏种质资源圃中选取生长状况良好、树龄均为15年左右且生长环境基本一致的银杏变异斑叶植株和正常绿叶植株各5株作为实验材料。从每株实验植株的中部外围枝条上，选取生长发育正常、大小基本一致且无病虫害的叶片。对于银杏变异斑叶植株，分别采集其黄色斑块部分（变异黄叶）和绿色部分（正常绿叶）的叶片，每种叶片各采集30片，迅速用液氮速冻后，置于-80℃超低温冰箱中保存，用于后续的生理生化指标测定和分子生物学实验；同时，另取5片新鲜叶片用于叶绿体超微结构观察。对于正常绿叶植株，采集30片正常绿色叶片，同样经液氮速冻后保存于-80℃超低温冰箱，另取5片新鲜叶片用于叶绿体超微结构观察。在采集过程中，确保每片叶片都具有代表性，且采集时间尽量集中在上午9:00-11:00，以减少因时间差异导致的生理状态变化对实验结果的影响。2.2形态学观察从2023年3月银杏萌芽开始，至12月叶片自然落叶结束，对选取的银杏变异斑叶植株和正常绿叶植株进行定期观察记录。每隔10天，使用卷尺测量叶片的长度、宽度，并计算叶片的长宽比；用游标卡尺测量叶柄的长度和直径，详细记录数据。同时，仔细观察叶片的形状、颜色、斑纹分布等特征，如叶片是完整的扇形还是有缺刻、裂片，颜色是均匀的绿色还是黄绿相间的斑纹，斑纹是呈规则的条纹状还是无规则的斑块状，以及斑纹在叶片上的具体分布位置和比例等，并采用拍照的方式进行直观记录，每次拍照时均选择相同的拍摄角度和光线条件，以确保照片的可比性。在观察过程中，特别关注叶片从萌芽到完全展开以及叶色从最初的嫩绿逐渐变化的过程。对于银杏变异斑叶植株，重点记录黄色斑块出现的时间、部位以及随着时间推移其面积的变化情况；对于正常绿叶植株，记录其叶片颜色在整个生长周期内的变化，以及是否出现异常的颜色变化或形态特征。此外，还对叶片的生长速率进行分析，通过计算相邻两次测量时叶片长度、宽度等指标的差值，除以测量间隔天数，得到叶片在不同生长阶段的日平均生长速率，以此来比较银杏变异斑叶与正常绿叶在生长动态上的差异。2.3细胞学观察2.3.1叶片解剖结构观察在叶片解剖结构观察中，采用石蜡切片法。于2023年7月中旬，从银杏变异斑叶植株和正常绿叶植株上，分别选取生长状况一致的成熟叶片各5片。将选取的叶片沿中脉迅速切成大小约为8mm×6mm的小片，立即投入FAA固定液中，固定时间不少于3天，确保叶片组织充分固定。固定液配方为：50%或70%酒精90mL、冰醋酸5mL、福尔马林5mL，考虑到银杏叶片质地，选用70%酒精配制固定液。固定完成后，进行洗涤操作，由于采用FAA固定液，可直接进入脱水步骤。脱水过程使用不同浓度梯度的酒精溶液，依次为50%、70%、85%、95%、100%酒精，每级脱水时间为1小时，特别注意在95%和100%酒精中放置时间不宜过长，若需过夜则放置于85%以下低浓度酒精中。脱水后，使用二甲苯作为透明剂进行透明处理，步骤如下：100%酒精→2/3酒精+1/3二甲苯→1/2酒精+1/2二甲苯→1/3酒精+2/3二甲苯→二甲苯→二甲苯，在每个溶液中停留1小时，使叶片组织充分透明。完成透明后，进行透蜡操作。将熔化的石蜡倒入盛有透明后的样品和二甲苯的熔蜡杯中，熔蜡与二甲苯比例为1：1。将熔蜡杯置于熔蜡箱中，设定恒温使石蜡刚好熔化，保持3-4小时。同时在熔蜡箱中放置三套纯石蜡小杯，每隔3小时将样品轻移至纯石蜡小杯中，以彻底除去组织内的透明剂。透蜡后，进行包埋。将材料和石蜡一起倒入小纸盒中，用镊子在酒精灯上烧热后，调整材料位置，使材料之间保持一定间距，并注意材料方向。用烧热的拨针沿着纸盒边缘及材料周围轻轻移动，赶走蜡中的气泡，使蜡凝固均匀。待蜡块凝固后，贴上标签，注明样品信息。使用旋转切片机将包埋好的蜡块切成厚度为8μm的切片，将切片粘贴在预先涂有粘片剂的载玻片上。将载玻片置于电热温台上，40℃烘烤30分钟，使切片牢固粘贴。然后进行脱蜡，将载玻片放入二甲苯中浸泡30分钟，依次经过1/2酒精+1/2二甲苯、100%酒精、95%酒精、85%酒精、70%酒精、50%酒精，每个步骤停留3分钟。染色采用1%番红（50%酒精配制）染色12小时以上，再用50%酒精冲洗10-30秒，接着依次经过70%酒精（2分钟）、85%酒精（2分钟）、95%酒精（2分钟）、1%固绿（95%酒精配制，染色10-20秒）、95%酒精（2分钟）、100%酒精（2分钟）、100%酒精（2分钟）、1/2酒精+1/2二甲苯（2分钟）、二甲苯（2分钟）、二甲苯（2分钟），最后用中性树胶封片。在OlympusBX53光学显微镜下观察切片，使用显微镜自带的成像系统拍照记录，对比分析银杏变异斑叶与正常绿叶的叶片厚度、表皮细胞形态与厚度、叶肉组织厚度及细胞排列方式、叶脉结构等解剖结构差异。2.3.2叶绿体超微结构观察在叶绿体超微结构观察实验中，使用透射电子显微镜（TEM）。于2023年7月中旬，从银杏变异斑叶植株和正常绿叶植株上，分别选取新鲜叶片各5片，迅速将叶片切成1mm×1mm左右的小块。将切好的叶片小块立即投入2.5%戊二醛固定液（用0.1mol/L磷酸缓冲液配制，pH7.2-7.4）中，4℃固定4小时。固定完成后，用0.1mol/L磷酸缓冲液冲洗3次，每次15分钟，以去除多余的固定液。然后用1%锇酸固定液（同样用0.1mol/L磷酸缓冲液配制）4℃固定2小时，再次用0.1mol/L磷酸缓冲液冲洗3次，每次15分钟。接下来进行脱水操作，依次用30%、50%、70%、80%、95%、100%的乙醇溶液进行脱水，每个浓度停留15分钟，100%乙醇溶液需脱水2次。脱水完成后，用环氧丙烷置换乙醇2次，每次15分钟。将样品置于环氧丙烷与包埋剂（Epon812）按1：1混合液中渗透2小时，再放入纯包埋剂中渗透过夜。第二天，将样品置于模具中，加入新鲜包埋剂，70℃聚合8小时，制成包埋块。使用超薄切片机将包埋块切成厚度约70nm的超薄切片，将切片捞至铜网上。用2%醋酸铀和柠檬酸铅进行双重染色，在JEOLJEM-1400透射电子显微镜下观察叶绿体的超微结构，加速电压为80kV。仔细观察并拍照记录叶绿体的形状、大小、数量，叶绿体膜的完整性，基粒片层的垛叠程度、数量和排列方式，类囊体膜的结构，质体小球的数量、大小和分布情况，淀粉粒的有无及大小、数量等特征。通过ImageJ软件对拍摄的电镜照片进行分析，测量叶绿体的长径、短径、基粒片层厚度、质体小球直径等参数，统计单位面积内叶绿体、基粒片层、质体小球的数量，对比分析银杏变异斑叶与正常绿叶叶绿体超微结构的差异。2.4色素含量测定采用乙醇-丙酮混合提取法测定银杏叶片中叶绿素、类胡萝卜素的含量。精确称取剪碎并混匀的银杏叶片样品（包括银杏变异斑叶的黄色斑块部分和绿色部分、正常绿叶）各0.2g，放入研钵中，加入少量石英砂和碳酸钙粉，以保护色素不被破坏。再加入10mL体积比为1：1的95%乙醇和丙酮混合提取液，在冰浴条件下迅速研磨成匀浆，以减少研磨过程中色素的分解。将匀浆转移至25mL棕色容量瓶中，用少量混合提取液冲洗研钵及杵棒2-3次，将冲洗液一并倒入容量瓶中，然后用混合提取液定容至刻度线。充分摇匀后，将容量瓶置于黑暗环境中浸提24小时，使色素充分溶解于提取液中。24小时后，将浸提液用定性滤纸过滤至干净的棕色试剂瓶中，得到叶绿体色素提取液。使用UV-2600紫外可见分光光度计，以95%乙醇-丙酮混合提取液作为空白对照，分别在波长665nm、649nm和470nm处测定提取液的吸光度。根据Arnon公式计算叶绿素a（Chla）、叶绿素b（Chlb）以及类胡萝卜素（Car）的含量。计算公式如下：Chla(mg/L)=12.21\timesA_{665}-2.81\timesA_{649}Chlb(mg/L)=20.13\timesA_{649}-5.03\timesA_{665}Car(mg/L)=\frac{1000\timesA_{470}-3.27\timesChla-104\timesChlb}{229}式中，A_{665}、A_{649}、A_{470}分别为提取液在波长665nm、649nm和470nm处的吸光度。最后，根据公式计算单位质量叶片中各色素的含量，公式为：色ç´

含量(mg/g)=\frac{C\timesV}{W\times1000}式中，C为色素浓度（mg/L），V为提取液总体积（mL），W为叶片样品质量（g）。每个样品设置3次生物学重复，取平均值作为最终结果，并计算标准偏差，以评估数据的可靠性和重复性。2.5转录组测序分析取适量保存于-80℃超低温冰箱中的银杏变异斑叶的黄色斑块部分（变异黄叶）和正常绿叶植株的正常绿色叶片样品，使用TRIzol试剂（Invitrogen公司）按照其说明书提供的操作步骤进行总RNA的提取。提取过程中，严格控制操作环境，避免RNA酶的污染，确保提取的RNA质量和完整性。提取完成后，使用1%琼脂糖凝胶电泳检测RNA的完整性，观察28S和18SrRNA条带的清晰度和亮度，以判断RNA是否有降解；利用Nanodrop2000超微量分光光度计（ThermoScientific公司）测定RNA的浓度和纯度，确保OD260/OD280比值在1.8-2.0之间，OD260/OD230比值大于2.0，以保证RNA的纯度满足后续实验要求。将质量合格的RNA样品送至专业的生物测序公司（如华大基因）进行文库构建和高通量测序。文库构建采用IlluminaTruSeqStrandedmRNALTSamplePrepKit试剂盒，具体步骤如下：首先利用带有Oligo（dT）的磁珠富集真核生物mRNA，对于银杏这样的植物，由于其mRNA结构特点，该方法能有效富集目标mRNA；接着使用FragmentationBuffer将mRNA打断成短片段，在逆转录酶的作用下，以短片段mRNA为模板合成第一链cDNA，随后合成第二链cDNA，合成过程中使用dUTP代替dTTP，以便后续区分正反向链；对双链cDNA进行末端修复、加A尾，并连接测序接头，经过PCR扩增后得到最终的测序文库。文库构建完成后，使用Agilent2100生物分析仪（AgilentTechnologies公司）对文库的片段大小进行检测，确保文库片段分布符合预期；采用Qubit2.0荧光定量仪（Invitrogen公司）对文库浓度进行精确定量。将合格的文库在IlluminaHiSeqXTen测序平台上进行双端测序，测序读长为150bp。测序过程中，严格控制测序反应条件，确保测序数据的质量和准确性。测序完成后，对下机数据进行质量控制和预处理。利用FastQC软件对原始测序数据进行质量评估，查看数据的碱基质量分布、GC含量分布、测序接头污染等情况。使用Trimmomatic软件对原始数据进行过滤，去除低质量碱基（质量值低于20的碱基）、含测序接头的序列以及N含量超过10%的reads，得到高质量的cleanreads。将cleanreads与银杏参考基因组（可选用已公布的高质量银杏基因组序列）进行比对，使用Hisat2软件进行比对分析，统计比对到基因组上的reads数量和比例，评估测序数据的比对效率。通过StringTie软件对每个样品比对上的reads进行转录本组装，将组装得到的转录本与已知的基因注释文件进行合并，得到最终的基因表达定量结果。利用DESeq2软件对银杏变异黄叶和正常绿叶的基因表达数据进行差异表达分析，以|log2(FoldChange)|≥1且padj<0.05作为筛选差异表达基因的标准，筛选出在银杏变异斑叶的黄色斑块部分和正常绿叶中表达差异显著的基因。对筛选出的差异表达基因进行功能注释和富集分析，将差异表达基因映射到GO（GeneOntology）数据库和KEGG（KyotoEncyclopediaofGenesandGenomes）数据库，利用clusterProfiler软件进行GO功能富集分析和KEGG通路富集分析，找出差异表达基因显著富集的生物学过程、分子功能、细胞组成以及相关的代谢通路，为深入探究银杏变异斑叶呈色的分子机制提供线索。2.6实时荧光定量PCR验证为了验证转录组测序结果的可靠性，从转录组测序筛选出的差异表达基因中，选取10个与色素合成、叶绿体发育相关的关键基因进行实时荧光定量PCR（qRT-PCR）验证。这10个基因包括参与叶绿素合成的CHLH（镁离子螯合酶H亚基基因）、CHLI（镁离子螯合酶I亚基基因），参与类胡萝卜素合成的PSY（八氢番茄红素合成酶基因）、ZDS（ζ-胡萝卜素脱氢酶基因），以及与叶绿体发育密切相关的FtsZ（原核细胞分裂蛋白基因同源物）、GLK（黄金2样转录因子基因）等。使用PrimeScriptRTreagentKitwithgDNAEraser（TaKaRa公司）反转录试剂盒，将提取的总RNA反转录合成cDNA。具体操作步骤如下：在冰上配制反应体系，总RNA1μg，5×gDNAEraserBuffer2μL，gDNAEraser1μL，RNaseFreedH₂O补足至10μL。将上述反应体系轻轻混匀，42℃孵育2min，以去除基因组DNA污染。短暂离心后，继续在冰上加入5×PrimeScriptBuffer24μL，PrimeScriptRTEnzymeMixI1μL，RTPrimerMix1μL，RNaseFreedH₂O4μL，总体积为20μL。将反应管置于PCR仪中，37℃孵育15min，85℃孵育5s，得到的cDNA产物保存于-20℃备用。根据NCBI数据库中已公布的银杏基因序列，使用PrimerPremier5.0软件设计qRT-PCR引物。引物设计原则如下：引物长度为18-25bp，GC含量在40%-60%之间，Tm值在58-62℃之间，避免引物二聚体和发夹结构的形成。引物序列由生工生物工程（上海）股份有限公司合成。合成后的引物用ddH₂O稀释至10μmol/L的工作浓度，保存于-20℃。采用SYBRPremixExTaqII（TaKaRa公司）试剂盒进行qRT-PCR反应。在冰上配制20μL反应体系，包括2×SYBRPremixExTaqII10μL，上下游引物（10μmol/L）各0.8μL，cDNA模板1μL，ROXReferenceDyeII（50×）0.4μL，ddH₂O7μL。将反应体系充分混匀后，转移至96孔板中，使用CFX96TouchReal-TimePCRDetectionSystem（Bio-Rad公司）进行扩增。反应程序设置如下：95℃预变性30s；95℃变性5s，60℃退火30s，共40个循环；在每个循环的退火阶段收集荧光信号。扩增结束后，进行熔解曲线分析，以确定扩增产物的特异性，熔解曲线分析程序为95℃15s，60℃1min，95℃15s。以银杏的Actin基因作为内参基因，用于校正和标准化目的基因的表达水平。每个样品设置3次技术重复，采用2^(-ΔΔCt)法计算目的基因的相对表达量。其中，ΔCt=Ct目的基因-CtActin，ΔΔCt=ΔCt变异黄叶-ΔCt正常绿叶。将qRT-PCR得到的基因相对表达量与转录组测序得到的基因表达量进行相关性分析，使用SPSS软件计算Pearson相关系数，以评估转录组测序结果的可靠性。三、银杏变异斑叶与正常绿叶的形态发育特征3.1生长周期的形态变化在萌芽期，约3月中旬，银杏变异斑叶植株与正常绿叶植株的芽体均开始萌动。芽体呈圆锥状，外被褐色鳞片，鳞片表面有一层薄薄的绒毛，起到保护芽体的作用。随着气温逐渐升高，芽体开始膨大，鳞片逐渐张开，露出内部嫩绿的幼叶。此时，变异斑叶植株与正常绿叶植株的幼叶在形态上并无明显差异，均为淡黄绿色，叶片质地柔软，表面光滑，叶片边缘呈波浪状，主脉和侧脉尚未完全分化清晰。进入展叶期，大约在4月初，正常绿叶植株的叶片迅速展开，叶片形状逐渐呈现出典型的扇形，边缘具有浅裂或波状缺刻，长度一般在5-8厘米之间，宽度为4-7厘米。叶片颜色由淡黄绿色逐渐转变为深绿色，表面富有光泽，主脉明显，向叶片边缘延伸，侧脉以二叉状平行排列，清晰可见。而银杏变异斑叶植株的叶片在展开过程中，开始出现明显的性状分化。叶片上逐渐出现黄绿相间的条纹状或斑块状区域，黄色部分颜色较浅，呈浅黄色或淡金黄色，绿色部分则保持正常的绿色。随着叶片的进一步展开，斑纹的分布范围逐渐扩大，且斑纹的形状和大小各异，有的呈长条状贯穿叶片，有的则形成不规则的斑块分布在叶片上。5-8月为旺盛生长期，正常绿叶植株的叶片生长迅速，叶面积不断增大，长度可达8-12厘米，宽度为7-10厘米。叶片质地变厚，颜色深绿且富有光泽，叶片的光合作用能力较强，为植株的生长提供充足的能量和物质。在这一时期，正常绿叶的形态基本稳定，叶片的大小、形状和颜色变化较小。银杏变异斑叶植株在旺盛生长期，斑纹的特征更加明显。黄色区域的面积进一步扩大，部分叶片的黄色区域甚至超过绿色区域。黄色部分的颜色也逐渐加深，呈现出金黄色，与绿色部分形成鲜明的对比，极具观赏价值。同时，叶片的厚度和质地也有所增加，但相较于正常绿叶，变异斑叶的质地略显柔软，叶片的韧性稍差。9-12月是衰老期，随着气温逐渐降低，光照时间缩短，正常绿叶植株的叶片开始逐渐变黄。首先从叶片边缘开始，黄色逐渐向叶片中心蔓延，叶片的绿色逐渐褪去，直至整个叶片变为金黄色。在这一过程中，叶片的质地逐渐变脆，失去光泽，叶片边缘开始卷曲。到11月下旬至12月，叶片逐渐干枯，从树枝上脱落。银杏变异斑叶植株在衰老期，黄色斑纹部分的颜色变化不明显，但绿色部分也开始逐渐变黄，最终整个叶片变为金黄色。与正常绿叶相比，变异斑叶在衰老期的叶片脱落时间稍早，且脱落过程较为迅速。到12月上旬，大部分变异斑叶已经脱落，而正常绿叶则在12月中旬左右才基本全部脱落。3.2形态差异总结在整个生长周期中，银杏变异斑叶与正常绿叶呈现出明显且具有规律性的形态差异。在萌芽期，二者的芽体外观基本一致，均被褐色鳞片包裹，鳞片上的绒毛也无明显区别。随着生长进程推进，进入展叶期后，差异逐渐显现。正常绿叶迅速展开并呈现典型的扇形，边缘的浅裂或波状缺刻较为规整，叶片长宽比例相对稳定，在5-8厘米长、4-7厘米宽的范围内，长宽比约为1.2-1.4。而变异斑叶在展开时就出现黄绿相间的斑纹，斑纹的形状和分布随机性较大，导致叶片整体形状不如正常绿叶规则，长宽比也因斑纹分布的不均衡而有所波动。在旺盛生长期，正常绿叶生长稳定，叶面积增大，质地变厚，颜色深绿且富有光泽，此时其叶片的各项形态指标趋于稳定，长度在8-12厘米之间，宽度为7-10厘米。银杏变异斑叶的斑纹特征进一步强化，黄色区域面积扩大，部分叶片黄色区域甚至超过绿色区域。由于黄色区域细胞结构和生理功能与绿色区域存在差异，导致变异斑叶质地略显柔软，韧性较差，在叶片大小上，虽然整体也有所增大，但受斑纹影响，不同部位的生长速率不一致，使得叶片大小的变化范围相对较广。进入衰老期，正常绿叶从边缘开始逐渐变黄，直至整片叶子变为金黄色，质地变脆，边缘卷曲，最终脱落。银杏变异斑叶的绿色部分同样逐渐变黄，但由于其本身具有黄色斑纹，在衰老过程中颜色变化的视觉差异相对较小。并且，变异斑叶的脱落时间稍早于正常绿叶，脱落过程更为迅速。综上所述，银杏变异斑叶与正常绿叶在整个生长周期中的形态差异，从叶片形状、大小、质地到颜色变化和脱落时间等方面均有体现。这些差异不仅反映了二者在生长发育过程中的不同生理状态，也暗示了它们在细胞结构、色素含量及相关基因表达调控等方面存在显著不同，为深入研究银杏变异斑叶的呈色机理提供了重要的形态学依据。四、银杏变异斑叶呈色的细胞学基础4.1叶片解剖结构差异通过石蜡切片法对银杏变异斑叶与正常绿叶的叶片解剖结构进行观察分析，结果显示二者在多个结构层次上存在显著差异（图1）。在表皮结构方面，正常绿叶的上、下表皮细胞均呈扁平状，排列紧密且整齐，犹如紧密排列的砖块，形成了一道有效的保护屏障。表皮细胞外壁具有明显的角质层，厚度均匀，能够有效防止水分散失和外界有害物质的侵入。上表皮细胞厚度约为[X]μm，下表皮细胞厚度约为[X]μm。而银杏变异斑叶的上表皮细胞形态不规则，部分细胞出现膨大或凹陷的现象，排列较为疏松，细胞之间的间隙明显增大。下表皮细胞同样存在形态异常的情况，且角质层厚度不均匀，在某些区域明显变薄。变异斑叶上表皮细胞厚度在[X]-[X]μm之间波动，下表皮细胞厚度在[X]-[X]μm之间变化，与正常绿叶相比，变异斑叶表皮细胞厚度的波动范围更大。叶肉组织是叶片进行光合作用的主要部位，正常绿叶的叶肉组织分化明显，可清晰分为栅栏组织和海绵组织。栅栏组织位于上表皮下方，由1-2层紧密排列的柱状细胞组成，细胞长轴与叶片表面垂直，如同紧密排列的栅栏，细胞内叶绿体含量丰富，叶绿体呈长椭圆形，排列整齐，基粒片层结构完整，垛叠紧密，能够高效地捕获光能，为光合作用提供充足的能量。海绵组织位于栅栏组织下方，细胞形状不规则，排列疏松，细胞间隙较大，有利于气体交换。海绵组织细胞内叶绿体含量相对较少，但叶绿体同样具有完整的结构。正常绿叶栅栏组织厚度约为[X]μm，海绵组织厚度约为[X]μm。银杏变异斑叶的叶肉组织分化不明显，栅栏组织和海绵组织的界限模糊。栅栏组织细胞排列疏松，细胞形状不规则，部分细胞出现扭曲或变形的情况。细胞内叶绿体数量明显减少，且叶绿体形态异常，呈不规则形状，叶绿体膜模糊不清，基粒片层垛叠程度较低，片层结构疏松分散，分化程度低，这使得叶绿体的光合作用功能受到严重影响。海绵组织细胞排列更为松散，细胞间隙进一步增大，细胞内叶绿体含量极少，且叶绿体结构受损严重。变异斑叶栅栏组织厚度在[X]-[X]μm之间，海绵组织厚度在[X]-[X]μm之间，与正常绿叶相比，变异斑叶叶肉组织的厚度整体变薄，且各部分厚度的稳定性较差。叶脉是叶片的重要组成部分，主要负责水分和养分的运输以及对叶片的机械支持。正常绿叶的主脉和侧脉结构清晰，维管束发达。主脉维管束由木质部和韧皮部组成，木质部位于上方，由导管、管胞等组成，负责将根部吸收的水分和无机盐向上运输到叶片的各个部位；韧皮部位于下方，由筛管、伴胞等组成，主要负责将叶片光合作用产生的有机物质向下运输到植物的其他部位。维管束周围有一层厚壁细胞组成的维管束鞘，起到保护和支持维管束的作用。侧脉结构与主脉相似，但维管束的规模逐渐减小。正常绿叶主脉维管束直径约为[X]μm。银杏变异斑叶的叶脉结构相对不发达，主脉和侧脉的维管束较小。木质部和韧皮部的细胞数量减少，导管和筛管的口径变小，这可能会影响水分和养分的运输效率。维管束鞘细胞的厚度也有所变薄，对维管束的保护和支持作用减弱。变异斑叶主脉维管束直径约为[X]μm，明显小于正常绿叶。综上所述，银杏变异斑叶与正常绿叶在叶片解剖结构上存在显著差异，这些差异可能直接或间接地影响了叶片的生理功能，如光合作用、水分和养分的运输等，进而对叶色的形成和表现产生重要影响。表皮结构的差异可能影响叶片对环境的适应能力和水分散失速率；叶肉组织的异常导致光合作用效率降低，影响叶片中光合产物的积累；叶脉结构的不发达可能导致水分和养分供应不足，无法满足叶片正常生长和色素合成的需求。这些解剖结构上的差异为深入探究银杏变异斑叶呈色机理提供了重要的细胞学基础。结构正常绿叶变异斑叶表皮上、下表皮细胞扁平，排列紧密整齐，外壁角质层均匀，上表皮细胞厚约[X]μm，下表皮细胞厚约[X]μm上表皮细胞形态不规则、排列疏松、间隙大，下表皮细胞类似，角质层厚度不均，上表皮细胞厚[X]-[X]μm，下表皮细胞厚[X]-[X]μm叶肉叶肉组织分化明显，栅栏组织1-2层柱状细胞紧密排列，叶绿体丰富、结构完整，海绵组织细胞不规则、排列疏松、间隙大，叶绿体较少但结构完整，栅栏组织厚约[X]μm，海绵组织厚约[X]μm叶肉组织分化不明显，栅栏组织和海绵组织界限模糊，栅栏组织细胞排列疏松、形状不规则、叶绿体少且形态异常、片层结构差，海绵组织细胞排列更松散、间隙更大、叶绿体极少且结构受损严重，栅栏组织厚[X]-[X]μm，海绵组织厚[X]-[X]μm叶脉主脉和侧脉结构清晰，维管束发达，木质部和韧皮部组成合理，维管束鞘厚壁细胞起保护支持作用，主脉维管束直径约[X]μm叶脉结构不发达，主脉和侧脉维管束小，木质部和韧皮部细胞数量减少、导管筛管口径变小，维管束鞘细胞变薄，主脉维管束直径约[X]μm图1银杏变异斑叶与正常绿叶叶片解剖结构对比（标尺=50μm）：A为正常绿叶横切面；B为变异斑叶横切面；1为上表皮；2为栅栏组织；3为海绵组织；4为下表皮；5为主脉维管束。4.2叶绿体超微结构差异4.2.1叶绿体形态与结构通过透射电子显微镜对银杏变异斑叶与正常绿叶的叶绿体超微结构进行观察（图2），发现二者在叶绿体的形态与结构方面存在显著差异。正常绿叶中，叶绿体形状规则，呈典型的长椭圆形，长度约为[X]μm，宽度约为[X]μm。叶绿体被两层清晰完整的膜所包裹，内膜和外膜界限分明，膜结构连续且光滑，厚度均匀，有效维持了叶绿体内部环境的稳定性。基粒片层结构高度发达，由多个扁平的类囊体紧密垛叠而成，类囊体之间排列整齐、紧密有序，宛如一摞摞整齐摆放的硬币。每个基粒通常包含[X]-[X]个类囊体，基粒片层的厚度约为[X]nm。类囊体膜上分布着丰富的光合色素和光合蛋白复合体，如光系统Ⅰ（PSⅠ）、光系统Ⅱ（PSⅡ）、Cytb6f复合物以及CF0-CF1ATP合酶复合物等，这些物质在光合作用的光反应过程中发挥着关键作用，能够高效地捕获光能，并将光能转化为化学能。此外，叶绿体基质中还分布着少量的质体小球，质体小球呈球形，直径约为[X]nm，它们在叶绿体的脂质代谢和抗氧化防御等过程中具有重要作用。银杏变异斑叶的叶绿体则呈现出明显的异常形态与结构。叶绿体形状不规则，部分呈扭曲状、哑铃状或不规则多边形，与正常的长椭圆形差异显著。叶绿体的大小也存在较大差异，长度在[X]-[X]μm之间波动，宽度在[X]-[X]μm之间变化，整体大小不均匀。叶绿体膜结构模糊不清，内膜和外膜部分区域出现破损、断裂或融合现象，导致膜的完整性遭到破坏，这可能会影响叶绿体与外界环境之间的物质交换和信号传递。基粒片层垛叠程度相对较低，类囊体排列疏松，相互之间的连接不紧密，部分类囊体出现膨胀、变形或解体的情况。每个基粒所含的类囊体数量明显减少，仅为[X]-[X]个，基粒片层的厚度也变薄，约为[X]nm。由于类囊体结构的异常，导致光合色素和光合蛋白复合体的分布和功能受到影响，进而降低了光合作用的效率。在变异斑叶的叶绿体中，质体小球数量明显增多，且排列密集，质体小球直径在[X]-[X]nm之间，比正常绿叶中的质体小球更大。大量质体小球的积累可能与叶绿体的发育异常以及脂质代谢紊乱有关。结构正常绿叶变异斑叶叶绿体形状长椭圆形，长约[X]μm，宽约[X]μm不规则，扭曲状、哑铃状或多边形，长[X]-[X]μm，宽[X]-[X]μm叶绿体膜双层膜清晰完整，内膜和外膜界限分明、连续光滑、厚度均匀膜结构模糊，内膜和外膜部分破损、断裂或融合基粒片层高度发达，由多个类囊体紧密垛叠，每个基粒含[X]-[X]个类囊体，厚约[X]nm，类囊体排列整齐紧密，膜上光合色素和蛋白复合体丰富垛叠程度低，类囊体排列疏松、连接不紧密，部分膨胀、变形或解体，每个基粒含[X]-[X]个类囊体，厚约[X]nm，光合色素和蛋白复合体分布与功能受影响质体小球少量，球形，直径约[X]nm数量明显增多、排列密集，直径[X]-[X]nm图2银杏变异斑叶与正常绿叶叶绿体超微结构对比（标尺=1μm）：A为正常绿叶叶绿体；B为变异斑叶叶绿体；1为叶绿体膜；2为基粒片层；3为质体小球。4.2.2质体小球与其他结构质体小球作为叶绿体中的重要结构，在银杏变异斑叶与正常绿叶中表现出明显的差异。在正常绿叶中，质体小球数量较少，它们均匀地分散在叶绿体基质中，与基粒片层和其他细胞器保持着相对稳定的空间位置关系。质体小球主要由脂类物质组成，包含甘油三酯、磷脂、类胡萝卜素等成分，这些成分在叶绿体的正常生理功能中发挥着重要作用。甘油三酯是质体小球的主要储能物质，在植物生长发育的特定阶段，如种子萌发和幼苗生长初期，当植物需要额外的能量供应时，甘油三酯可以被分解代谢，为细胞提供能量；磷脂是生物膜的重要组成成分，质体小球中的磷脂可以参与叶绿体膜的修复和更新，维持叶绿体膜的完整性和稳定性；类胡萝卜素不仅是光合作用中的辅助色素，能够吸收和传递光能，拓宽植物对光的吸收范围，增强光合作用效率，还具有抗氧化功能，能够清除叶绿体在光合作用过程中产生的活性氧自由基，保护叶绿体免受氧化损伤。在银杏变异斑叶中，质体小球数量显著增加，它们在叶绿体基质中大量聚集，甚至出现相互融合的现象。质体小球的分布不再均匀，部分区域质体小球密集堆积，导致叶绿体内部结构紊乱。这种数量和分布的变化与叶绿体的发育异常密切相关。由于叶绿体膜结构受损和基粒片层发育不良，使得叶绿体的正常代谢过程受到干扰，从而影响了质体小球的合成、降解和运输平衡。一方面，可能是质体小球的合成过程增强，导致其数量增多；另一方面，质体小球的降解和运输途径受阻，使得它们在叶绿体中大量积累。质体小球的异常积累可能会进一步影响叶绿体的功能。过多的质体小球占据了叶绿体内部的空间，挤压了基粒片层和其他细胞器，影响了光合色素和光合蛋白复合体的正常分布和功能，进而降低了光合作用效率。此外，质体小球中富含的类胡萝卜素虽然具有一定的抗氧化作用，但过多的类胡萝卜素积累可能会改变叶绿体内部的氧化还原状态，影响叶绿体的正常代谢和信号传导。除了质体小球的变化，银杏变异斑叶的叶绿体中还出现了其他细胞结构的异常。部分叶绿体中淀粉粒的含量明显减少，甚至在一些叶绿体中几乎检测不到淀粉粒的存在。淀粉粒是植物光合作用的产物，它的积累和降解与光合作用的强度以及植物的生长发育阶段密切相关。在正常情况下，叶绿体通过光合作用产生的糖类物质会转化为淀粉粒进行储存，当植物需要能量时，淀粉粒再被分解为糖类供细胞利用。银杏变异斑叶中淀粉粒含量的减少，表明其光合作用产生的糖类物质可能无法正常积累，这可能是由于叶绿体结构和功能的异常，导致光合作用效率降低，糖类合成减少，或者是糖类的代谢和运输途径发生改变，使得合成的糖类不能及时转化为淀粉粒储存起来。此外，在变异斑叶的叶绿体周围，还观察到一些线粒体的形态和分布发生了变化。线粒体是细胞进行有氧呼吸的主要场所，为细胞提供能量。正常情况下，线粒体呈椭圆形或棒状，均匀地分布在细胞质中。而在银杏变异斑叶中，部分线粒体出现肿胀、变形的现象，线粒体嵴的数量减少且排列紊乱。线粒体的这种变化可能会影响其呼吸功能，导致细胞能量供应不足。叶绿体和线粒体在细胞的能量代谢中相互关联，叶绿体通过光合作用产生的ATP和糖类物质为线粒体的呼吸作用提供底物，而线粒体呼吸作用产生的ATP又为叶绿体的光合作用和其他生理过程提供能量。因此，叶绿体的发育异常可能会间接影响线粒体的功能，反之亦然。银杏变异斑叶中叶绿体和线粒体结构和功能的相互影响，进一步说明了叶色变异与细胞内能量代谢和物质代谢的密切关系。五、银杏变异斑叶呈色的生理机制5.1色素含量变化5.1.1叶绿素含量差异通过乙醇-丙酮混合提取法对银杏变异斑叶和正常绿叶中的叶绿素含量进行测定，结果显示二者存在显著差异（图3）。正常绿叶中叶绿素a含量较高，平均值为[X]mg/g，叶绿素b含量相对较低，平均值为[X]mg/g，总叶绿素含量为叶绿素a与叶绿素b含量之和，约为[X]mg/g。叶绿素a在光合作用中主要负责吸收和传递光能，将光能转化为化学能，其含量的高低直接影响光合作用的光反应效率。叶绿素b则主要起到辅助吸收光能的作用，拓宽植物对光的吸收范围，增强光合作用的稳定性。在正常绿叶中，叶绿素a与叶绿素b的比例相对稳定，约为[X]，这种比例关系有助于维持叶绿体中光合色素蛋白复合体的正常结构和功能，保证光合作用的高效进行。银杏变异斑叶的黄色斑块部分，叶绿素a含量显著降低，仅为[X]mg/g，约为正常绿叶中叶绿素a含量的[X]%；叶绿素b含量同样大幅下降，为[X]mg/g，约为正常绿叶中叶绿素b含量的[X]%；总叶绿素含量降至[X]mg/g，仅为正常绿叶总叶绿素含量的[X]%。叶绿素含量的显著降低使得变异斑叶的黄色部分对光能的吸收和利用能力大幅下降，光合作用受到严重抑制。这可能是由于变异斑叶中叶绿体发育异常，叶绿体数量减少、形态不规则，基粒片层垛叠程度低，类囊体膜结构受损，导致参与叶绿素合成的关键酶活性降低，叶绿素合成途径受阻。同时，叶绿素的降解过程可能也有所增强，进一步加剧了叶绿素含量的减少。例如，在叶绿素合成过程中，镁离子螯合酶是催化镁离子插入原卟啉IX形成镁原卟啉IX的关键酶，其活性的降低会直接影响叶绿素的合成效率。而在叶绿素降解过程中，叶绿素酶的活性升高可能会加速叶绿素的分解。此外，相关基因的表达变化也可能对叶绿素的合成和降解产生调控作用。叶片类型叶绿素a（mg/g）叶绿素b（mg/g）总叶绿素（mg/g）叶绿素a/b正常绿叶[X][X][X][X]变异斑叶（黄叶部分）[X][X][X][X]图3银杏变异斑叶与正常绿叶叶绿素含量对比：不同小写字母表示差异显著（P<0.05）。5.1.2类胡萝卜素含量差异类胡萝卜素作为植物体内重要的色素之一，在银杏变异斑叶和正常绿叶中的含量也存在明显差异（图4）。正常绿叶中类胡萝卜素含量相对较低，平均值为[X]mg/g。类胡萝卜素在植物光合作用中扮演着重要角色，它不仅能够吸收和传递光能，将吸收的光能传递给叶绿素，参与光合作用的光反应过程，还具有抗氧化功能，能够清除光合作用过程中产生的活性氧自由基，保护叶绿体免受氧化损伤，维持叶绿体的正常结构和功能。在正常生理状态下，类胡萝卜素与叶绿素的含量保持一定的比例关系，共同维持植物叶片的正常绿色和光合作用的稳定进行。银杏变异斑叶的黄色斑块部分，类胡萝卜素含量显著升高，平均值达到[X]mg/g，约为正常绿叶中类胡萝卜素含量的[X]倍。类胡萝卜素主要包括胡萝卜素和叶黄素两大类，其中胡萝卜素呈橙黄色，叶黄素呈黄色。在银杏变异斑叶的黄色区域，类胡萝卜素含量的增加，尤其是叶黄素含量的相对升高，使得叶片呈现出明显的黄色。这可能是由于在变异斑叶中，类胡萝卜素合成途径中的关键基因表达上调，促进了类胡萝卜素的合成。例如，八氢番茄红素合成酶（PSY）基因是类胡萝卜素合成途径的关键基因，其表达水平的升高会导致八氢番茄红素合成增加，进而促进类胡萝卜素的合成。同时，类胡萝卜素的降解过程可能受到抑制，使得类胡萝卜素在叶片中积累。此外，叶绿体发育异常导致的质体小球数量增多，也可能与类胡萝卜素含量的升高有关。质体小球是类胡萝卜素的储存场所，质体小球数量的增加为类胡萝卜素的积累提供了更多的空间。类胡萝卜素含量的升高以及其与叶绿素含量比例的改变，是导致银杏变异斑叶呈现黄色的重要生理原因之一。叶片类型类胡萝卜素（mg/g）正常绿叶[X]变异斑叶（黄叶部分）[X]图4银杏变异斑叶与正常绿叶类胡萝卜素含量对比：不同小写字母表示差异显著（P<0.05）。综上所述，银杏变异斑叶与正常绿叶在叶绿素和类胡萝卜素含量上存在显著差异。变异斑叶中叶绿素含量的显著降低和类胡萝卜素含量的显著升高，导致叶片中色素组成和比例发生改变，使得叶片呈现出黄绿相间的斑纹。这些色素含量的变化与叶绿体的发育异常密切相关，进一步深入研究色素合成和代谢相关基因的表达调控机制，对于揭示银杏变异斑叶呈色的分子机理具有重要意义。5.2色素含量与叶色的关系叶绿素作为植物叶片中最重要的光合色素，其含量的变化对叶色起着决定性作用。在正常绿叶中，叶绿素含量丰富，使得叶片呈现出浓郁的绿色。这是因为叶绿素能够吸收光能，参与光合作用的光反应过程，将光能转化为化学能，为植物的生长和发育提供能量。同时，叶绿素的绿色能够掩盖其他色素的颜色，使得叶片在外观上主要呈现出绿色。而在银杏变异斑叶的黄色斑块部分，叶绿素a和叶绿素b的含量均显著降低。叶绿素a是光合作用中光系统Ⅰ和光系统Ⅱ的核心组成部分，其含量的减少直接影响了光合作用中光能的捕获和转化效率。叶绿素b作为辅助色素，能够吸收不同波长的光能，并将能量传递给叶绿素a，其含量的降低也进一步削弱了叶片对光能的利用能力。叶绿素含量的大幅下降，使得叶片对绿光的吸收能力减弱，绿光被大量反射，从而导致叶片绿色变浅。当叶绿素含量降低到一定程度时，其他色素的颜色便得以显现，这是银杏变异斑叶黄色部分形成的重要原因之一。类胡萝卜素在银杏变异斑叶的呈色过程中也发挥着关键作用。正常情况下，类胡萝卜素在叶片中的含量相对较低，其颜色被叶绿素所掩盖。但在银杏变异斑叶的黄色区域，类胡萝卜素含量显著升高。类胡萝卜素主要包括胡萝卜素和叶黄素，其中叶黄素呈黄色。随着类胡萝卜素含量的增加，尤其是叶黄素含量的相对升高，叶片对蓝光和绿光的吸收能力增强，而对红光和橙光的吸收相对减少，使得叶片呈现出明显的黄色。类胡萝卜素不仅赋予了银杏变异斑叶独特的颜色，还在光合作用中起到保护叶绿素的作用。在强光条件下，类胡萝卜素能够吸收多余的光能，将其转化为热能散失掉，避免叶绿素因吸收过多光能而受到光氧化损伤。此外，类胡萝卜素还参与了植物的抗氧化防御系统，能够清除光合作用过程中产生的活性氧自由基，维持叶绿体的正常结构和功能。因此，在银杏变异斑叶中，类胡萝卜素含量的升高不仅改变了叶片的颜色，还可能对叶片的光合作用和抗氧化能力产生重要影响。综上所述，银杏变异斑叶独特的黄绿相间颜色是叶绿素和类胡萝卜素含量变化共同作用的结果。叶绿素含量的显著降低削弱了叶片对绿光的吸收能力，使得叶片绿色变浅；而类胡萝卜素含量的显著升高，尤其是叶黄素含量的相对增加，增强了叶片对蓝光和绿光的吸收，使得叶片呈现出黄色。这两种色素含量的动态变化，打破了正常叶片中色素的平衡，从而导致银杏变异斑叶呈现出独特的颜色。深入研究色素含量变化的调控机制，对于揭示银杏变异斑叶呈色的分子机理具有重要意义。六、银杏变异斑叶呈色的转录调控机制6.1转录组测序数据分析对银杏变异斑叶的黄色斑块部分（变异黄叶）和正常绿叶植株的正常绿色叶片进行转录组测序后，获得了大量的原始测序数据。经过严格的数据质量控制和预处理，利用FastQC软件对原始数据进行质量评估，结果显示碱基质量值分布良好，平均质量值均在30以上，表明测序数据的准确性较高；GC含量分布在45%-50%之间，处于正常范围，无明显的GC偏好性；测序接头污染率极低，低于0.1%，说明数据质量可靠，可用于后续分析。使用Trimmomatic软件过滤低质量碱基、含测序接头的序列以及N含量超过10%的reads后，得到高质量的cleanreads。变异黄叶样本获得的cleanreads数量为[X]条，正常绿叶样本获得的cleanreads数量为[X]条，两个样本的cleanreads占原始reads的比例均达到95%以上。将cleanreads与银杏参考基因组进行比对，使用Hisat2软件进行比对分析，结果显示变异黄叶样本的比对率为[X]%，正常绿叶样本的比对率为[X]%，表明大部分cleanreads能够成功比对到银杏参考基因组上，为后续的基因表达分析提供了可靠的数据基础。通过StringTie软件对每个样品比对上的reads进行转录本组装，将组装得到的转录本与已知的基因注释文件进行合并，得到最终的基因表达定量结果。利用DESeq2软件对银杏变异黄叶和正常绿叶的基因表达数据进行差异表达分析，以|log2(FoldChange)|≥1且padj<0.05作为筛选差异表达基因的标准，共筛选出[X]个差异表达基因，其中上调表达的基因有[X]个，下调表达的基因有[X]个。为了直观展示差异表达基因的分布情况，绘制了火山图（图5）。在火山图中，横坐标表示基因表达量的变化倍数（log2(FoldChange)），纵坐标表示差异表达的显著性水平（-log10(padj)）。图中红色点表示上调表达的差异基因，蓝色点表示下调表达的差异基因，黑色点表示表达无显著差异的基因。从火山图中可以清晰地看出，上调表达和下调表达的差异基因分布在两侧，且数量较多，表明银杏变异斑叶的黄色斑块部分与正常绿叶在基因表达水平上存在显著差异。图5银杏变异黄叶与正常绿叶差异表达基因火山图：红色点表示上调表达的差异基因，蓝色点表示下调表达的差异基因，黑色点表示表达无显著差异的基因。对差异表达基因进行层次聚类分析，绘制了聚类热图（图6）。聚类热图以颜色深浅来表示基因表达量的高低，红色表示高表达，蓝色表示低表达。通过聚类分析，可以将差异表达基因分为不同的簇，每个簇中的基因具有相似的表达模式。从聚类热图中可以看出，银杏变异黄叶和正常绿叶的基因表达模式明显不同，且差异表达基因在不同样本中呈现出较为明显的聚类特征，进一步验证了两者在基因表达水平上的显著差异。这也表明，这些差异表达基因可能在银杏变异斑叶的呈色过程中发挥着重要作用。图6银杏变异黄叶与正常绿叶差异表达基因聚类热图：每一行代表一个差异表达基因，每一列代表一个样本，红色表示高表达，蓝色表示低表达。6.2差异表达基因的功能注释与富集分析6.2.1功能注释将筛选出的[X]个差异表达基因映射到多个权威数据库中进行功能注释，以期全面、深入地了解这些基因的生物学功能。在GO（GeneOntology）数据库注释结果中，从生物过程（BiologicalProcess）角度分析，众多差异表达基因参与了光合作用相关过程。其中，有[X]个基因参与了光合生物中光能捕获过程，如编码捕光叶绿素a/b结合蛋白的基因，这些基因在正常绿叶中高表达，而在银杏变异斑叶的黄色斑块部分表达显著下调。捕光叶绿素a/b结合蛋白能够捕获光能，并将其传递给光系统中的反应中心，在光合作用的光反应过程中起着关键作用。此外，有[X]个基因参与了光合作用光反应中的电子传递链，如编码细胞色素b6/f复合体亚基的基因，其表达变化可能影响电子传递效率，进而影响光合作用的进行。在碳固定过程中，也发现了[X]个相关差异表达基因，如编码1,5-二磷酸核酮糖羧化酶/加氧酶（Rubisco）小亚基的基因，Rubisco是光合作用碳固定的关键酶，其亚基基因表达的改变可能对碳固定速率产生重要影响。从分子功能（MolecularFunction）层面来看，许多差异表达基因具有结合活性。其中，[X]个基因具有叶绿素结合活性，这些基因编码的蛋白能够与叶绿素特异性结合，维持叶绿素在叶绿体中的稳定存在，并参与叶绿素在光合作用中的光能捕获和传递过程。另外，[X]个基因具有辅酶结合活性，如编码参与光合作用辅酶合成或代谢相关酶的基因，辅酶在光合作用的酶促反应中起着不可或缺的作用，参与能量传递和物质转化。还有[X]个基因具有氧化还原酶活性，这些基因参与了光合作用过程中的氧化还原反应，如光系统Ⅰ和光系统Ⅱ中的电子传递和氧化还原过程，对维持光合作用中电子传递链的正常运行至关重要。在细胞组成（CellularComponent）方面，大量差异表达基因与叶绿体相关。[X]个基因参与了叶绿体类囊体膜的组成，类囊体膜是光合作用光反应的主要场所，这些基因编码的蛋白构成了类囊体膜的结构框架，并参与了光合色素和蛋白复合体的组装与定位。[X]个基因与叶绿体基质相关，叶绿体基质中进行着光合作用的暗反应，这些基因可能参与了暗反应中碳固定、糖类合成等过程的调控。此外，还有[X]个基因与叶绿体被膜有关，叶绿体被膜对维持叶绿体的结构完整性和物质交换起着重要作用，其相关基因表达的变化可能影响叶绿体与细胞质之间的物质和能量交流。在KEGG（KyotoEncyclopediaofGenesandGenomes）数据库注释中，发现差异表达基因显著富集在多个与光合作用和色素代谢相关的通路中。在光合作用通路中，涉及光系统Ⅰ、光系统Ⅱ、细胞色素b6/f复合体、ATP合酶等多个关键组件的基因表达发生显著变化。例如，编码光系统Ⅱ中D1蛋白的基因在银杏变异斑叶的黄色斑块部分表达下调，D1蛋白是光系统Ⅱ反应中心的核心蛋白，其表达降低可能导致光系统Ⅱ活性下降，影响光能的捕获和转化。在光合电子传递链中，编码细胞色素b6/f复合体亚基的基因表达异常，这可能干扰电子传递过程，进而影响ATP和NADPH的合成，而ATP和NADPH是光合作用暗反应所必需的能量和还原力。在色素代谢通路中，叶绿素代谢通路中的多个关键基因表达受到显著调控。如参与叶绿素合成的镁离子螯合酶（Mg-chelatase）基因家族成员CHLH、CHLI表达下调，镁离子螯合酶是催化镁离子插入原卟啉IX形成镁原卟啉IX的关键酶，其基因表达下调可能导致镁原卟啉IX合成受阻，进而影响叶绿素的合成。同时，叶绿素降解相关基因NYC1（Non-YellowColoring1）表达上调，NYC1基因编码的蛋白参与叶绿素b向叶绿素a的转化以及叶绿素的降解过程，其表达上调可能加速叶绿素的降解，导致银杏变异斑叶中叶绿素含量降低。在类胡萝卜素代谢通路中，八氢番茄红素合成酶（PSY）基因、ζ-胡萝卜素脱氢酶（ZDS）基因等表达上调，PSY是类胡萝卜素合成途径的第一个关键酶，催化两分子的牻牛儿基牻牛儿基焦磷酸（GGPP）缩合形成八氢番茄红素，ZDS则参与八氢番茄红素向番茄红素的转化过程，这些基因表达上调可能促进类胡萝卜素的合成，使得银杏变异斑叶中类胡萝卜素含量升高。6.2.2富集分析通过GO功能富集分析，进一步揭示了差异表达基因在生物学过程、分子功能和细胞组成等方面的显著富集情况（图7）。在生物学过程类别中，光合作用相关的生物过程显著富集。除了前面提到的光能捕获、光合电子传递链和碳固定过程外，还包括光合作用的调节过程。有多个差异表达基因参与了光合作用的调节，如编码光合调控蛋白的基因，这些基因可能通过调控光合色素蛋白复合体的组装、活性以及光合作用相关酶的表达和活性，来调节光合作用的速率和效率。此外，与氧化还原过程相关的生物过程也显著富集，这与光合作用过程中伴随着大量的氧化还原反应密切相关。在光合作用的光反应和暗反应中，都涉及到电子的传递和氧化还原反应，如光系统Ⅰ和光系统Ⅱ中的电子传递、NADP+的还原以及卡尔文循环中的氧化还原反应等，这些过程需要一系列氧化还原酶和电子载体的参与，而差异表达基因在这些相关过程中的富集，表明它们可能对光合作用中的氧化还原平衡和电子传递过程产生重要影响。在分子功能类别中，除了结合活性和氧化还原酶活性显著富集外，还发现与转运活性相关的功能显著富集。有多个差异表达基因编码的蛋白具有离子转运活性，如镁离子转运蛋白基因，镁离子是叶绿素分子的中心原子，其转运活性的改变可能影响镁离子在细胞内的分布和浓度，进而影响叶绿素的合成。此外，还有一些基因编码的蛋白具有小分子转运活性，可能参与了光合作用相关的辅酶、代谢产物等小分子物质的转运过程，对光合作用的正常进行起着重要的支持作用。在细胞组成类别中，除了叶绿体相关的细胞组成显著富集外，还发现与膜结构相关的细胞组成显著富集。除了叶绿体类囊体膜和被膜外，还包括其他细胞器膜以及细胞膜。这可能是因为叶色变异不仅影响了叶绿体的结构和功能，还可能对细胞内其他膜系统的结构和功能产生间接影响。例如，细胞膜的结构和功能变化可能影响细胞与外界环境之间的物质交换和信号传递，进而影响细胞的生理状态和代谢过程。细胞器膜的变化可能影响细胞器之间的物质运输和信息交流，对细胞内的能量代谢和物质合成等过程产生影响。在KEGG通路富集分析中，除了光合作用和色素代谢通路外，还发现差异表达基因在植物激素信号转导通路中显著富集（图8）。在植物激素信号转导通路中，涉及生长素、细胞分裂素、脱落酸、乙烯等多种植物激素的信号转导过程。例如，生长素信号转导通路中，编码生长素响应因子（ARF）的基因表达发生变化，ARF能够与生长素响应元件结合，调控下游基因的表达，参与植物的生长发育过程，如细胞伸长、分化和器官形成等。在细胞分裂素信号转导通路中，编码细胞分裂素响应调节因子（RR）的基因表达改变，细胞分裂素在植物的细胞分裂、分化、衰老等过程中发挥重要作用，RR基因表达的变化可能影响细胞分裂素信号的传递和响应，进而影响植物叶片的生长和发育。脱落酸信号转导通路中，编码脱落酸受体（PYR/PYL/RCAR）的基因表达异常，脱落酸在植物应对逆境胁迫以及种子休眠、萌发等过程中起重要作用，其受体基因表达的改变可能影响植物对逆境的响应以及叶片的衰老进程。乙烯信号转导通路中，编码乙烯响应因子（ERF）的基因表达变化，乙烯参与植物的生长发育、果实成熟、衰老等过程，ERF基因表达的改变可能对银杏叶片的生长和叶色变化产生影响。植物激素信号转导通路的异常可能通过调控相关基因的表达，影响叶绿体的发育、色素的合成与代谢以及光合作用等过程，进而参与银杏变异斑叶的呈色调控。综上所述，通过GO功能富集分析和KEGG通路富集分析，揭示了银杏变异斑叶与正常绿叶之间差异表达基因在生物学过程、分子功能、细胞组成以及代谢通路等方面的显著富集情况。这些富集结果为深入理解银杏变异斑叶呈色的分子机制提供了重要线索，表明银杏变异斑叶的呈色过程涉及光合作用、色素代谢、植物激素信号转导等多个生物学过程的协同调控，以及叶绿体结构和功能、膜系统等细胞组成的变化。图7银杏变异黄叶与正常绿叶差异表达基因GO富集分析柱状图：横坐标表示富集的GOterm，纵坐标表示富集到该GOterm的差异表达基因数量。不同颜色柱子分别表示生物过程、分子功能和细胞组成三个类别。图8银杏变异黄叶与正常绿叶差异表达基因KEGG通路富集分析气泡图：横坐标表示富集因子，即差异表达基因在某一通路中的富集程度，富集因子越大，说明该通路在差异表达基因中富集越显著；纵坐标表示KEGG通路名称；气泡大小表示富集到该通路的差异表达基因数量，气泡越大，基因数量越多；气泡颜色表示p-value值的大小，颜色越红，p-value值越小，富集越显著。6.3色素合成代谢相关基因的表达分析6.3.1叶绿素合成代谢基因叶绿素的合成是一个复杂且精细调控的过程，涉及多个关键基因的协同作用。在本研究中，通过转录组测序和实时荧光定量PCR验证，对银杏变异斑叶和正常绿叶中叶绿素合成代谢通路关键基因的表达进行了深入分析。在叶绿素合成的起始阶段，谷氨酸经过一系列反应转化为5-氨基乙酰丙酸（ALA），这一过程由谷氨酸-tRNA还原酶（GTR）和谷氨酸-1-半醛转氨酶（GSA）催化，编码这两个酶的基因在正常绿叶中表达水平相对稳定。然而，在银杏变异斑叶的黄色斑块部分，GTR基因的表达量显著下调，约为正常绿叶的[X]%，GSA基因的表达量也有所下降，降至正常绿叶的[X]%。这表明在变异斑叶中，ALA的合成受到抑制，可能是由于相关基因表达下调导致酶活性降低，进而影响了叶绿素合成的起始步骤。从ALA合成胆色素原（PBG）的过程由ALA脱水酶（ALAD）催化，在正常绿叶中，ALAD基因表达正常。但在银杏变异斑叶中，ALAD基因表达下调，其表达量仅为正常绿叶的[X]%。PBG进一步经过一系列反应生成原卟啉IX（ProtoIX），这一过程涉及多个酶的参与，如胆色素原脱氨酶（PBGD）、尿卟啉原III合酶（UROS）、尿卟啉原III脱羧酶（UROD）和粪卟啉原III氧化酶（CPOX）。在银杏变异斑叶中，PBGD基因表达下调，表达量为正常绿叶的[X]%；UROS基因表达量降至正常绿叶的[X]%；UROD基因表达下调更为明显，仅为正常绿叶的[X]%；CPOX基因表达量也显著降低，为正常绿叶的[X]%。这些基因表达的下调可能导致原卟啉IX合成受阻，影响叶绿素合成的后续步骤。原卟啉IX是叶绿素和血红素合成的共同前体，在叶绿素合成途径中，原卟啉IX在镁离子螯合酶（Mg-chelatase）的作用下，插入镁离子形成镁原卟啉IX（Mg-ProtoIX）。镁离子螯合酶由CHLH、CHLD和CHLI三个亚基组成，在正常绿叶中，这三个亚基基因均正常表达。但在银杏变异斑叶中，CHLH基因表达显著下调，表达量仅为正常绿叶的[X]%；CHLD基因表达量降至正常绿叶的[X]%；CHLI基因表达下调明显，为正常绿叶的[X]%。镁离子螯合酶亚基基因表达的下调，使得镁原卟啉IX合成减少，这是导致叶绿素合成受阻的关键环节之一。镁原卟啉IX经过一系列甲基化、环化等反应，最终形成叶绿素。在这一过程中，编码镁原卟啉IX甲基转移酶（CHLM）、镁原卟啉IX单甲酯环化酶（CRD1）等酶的基因在银杏变异斑叶中表达也出现不同程度的下调。CHLM基因表达量为正常绿叶的[X]%，CRD1基因表达量降至正常绿叶的[X]%。这些基因表达的改变进一步影响了叶绿素合成的后续反应，导致叶绿素合成量减少。综上所述，银杏变异斑叶中叶绿素合成代谢通路关键基因的表达显著下调，从ALA的合成到叶绿素合成的各个步骤均受到抑制，这是导致变异斑叶中叶绿素含量降低的重要分子机制。这些基因表达的变化可能与叶绿体发育异常、转录因子调控以及信号转导途径的改变等因素有关，进一步深入研究这些基因的调控机制，对于揭示银杏变异斑叶呈色的分子机理具有重要意义。6.3.2类胡萝卜素合成代谢基因类胡萝卜素的合成同样是一个复杂的生化过程，涉及多个关键基因的有序表达和协同调控。通过对银杏变异斑叶和正常绿叶中类胡萝卜素合成代谢相关基因的表达分析，发现这些基因的表达变化与变异斑叶中类胡萝卜素含量的升高密切相关。类胡萝卜素合成的起始阶段是由异戊烯基焦磷酸（IPP）和二甲基烯丙基焦磷酸（DMAPP）在牻牛儿基牻牛儿基焦磷酸合成酶（GGPS）的催化下，逐步缩合形成牻牛儿基牻牛儿基焦磷酸（GGPP）。在正常绿叶中，GGPS基因表达维持在一定水平。而在银杏变异斑叶的黄色斑块部分，GGPS基因表达显著上调，表达量约为正常绿叶的[X