银杏小孢子母细胞响应秋水仙碱低频多倍化的细胞学和转录组分析.pdf

Journal of NuclearAgricultural Sciences1497核农学报2 0 2 3,37(8)：1497 150 6文章编号：10 0 0-8 551(2 0 2 3)0 8-1497-10银杏小孢子母细胞响应秋水仙碱低频多倍化的细胞学和转录组分析柯梦1邓厚银1郭文斌2李艳星？米跃骐1李云1孙宇涵1，（北京林业大学生物科学与技术学院/国家林业和草原局刺槐工程技术研究中心/林木育种与生态修复国家工程研究中心/林木、花卉遗传育种教育部重点实验室，北京100083；襄垣县森防检疫站，山西长治0 46 2 0 0）摘要：银杏（Ginkgo bilobaL.）三倍体具有很高的综合价值，但秋水仙碱诱导银杏有性多倍化仍较为困难，2 n雄配子最高得率仅有7%。为初步探究秋水仙碱诱导银杏2 n雄配子的形成机制，本研究使用0.6%的秋水仙碱对进入减数分裂粗线期的银杏小孢子母细胞进行避光瓶浸没处理，对不同处理时间下小孢子囊细胞中的微管骨架、生理生化变化以及转录本进行对比分析。结果显示，秋水仙碱处理的银杏小孢子母细胞游离脯氨酸、丙二醛含量及过氧化物酶活性在处理48 h时较高，可溶性蛋白含量和超氧化物酶活性低于对照组。与对照组相比，处理至2 4h时，细胞微管蛋白开始减少，48 h呈弥散状，厚度变薄，荧光显微下几乎无微管蛋白显现。转录组分析共筛选到46 90 个差异表达基因，其中热激基因HSP11在48 h高差异表达。另外，在48 h鉴定到的41个差异转录因子中，有2 3个转录因子显著上调表达，其中最大的三个家族是AP2/ERF、M YB-r e l a t e d 和NAC蛋白家族。上述结果表明，48 h是秋水仙碱诱导银杏小孢子母细胞加倍的关键时间点，HSP11基因、AP2/ERF、M YB-r e l a t e d 和NAC转录因子家族成员在此过程中均起着关键的调控作用。本研究结果为探索与调控秋水仙碱诱导银杏2 n雄配子形成有关的通路提供了参考依据，并为后续银杏有性多倍化分子机理研究奠定了理论基础。关键词：银杏；倍性育种；转录组；秋水仙碱；花粉母细胞D0I:10.11869/j.issn.1000-8551.2023.08.1497多倍化是指生物体中多套染色体在同一细胞核内稳定遗传，产生拥有3套或以上完整染色体组个体的过程,是植物进化的重要因素1。自然界中至少一半的被子植物都经历过多倍化事件2 ,而在裸子植物中，银杏的多倍化研究也一直是重点关注内容3。通过形成2 n配子参与受精的有性多倍化是自然界多倍体产生的主要途径4，但自然发生多倍体的频率并不高。人工诱导能极大提高多倍体得率，目前已有银灰杨(Populus canescens)5 杜仲(Eucommia ulmoides)6等多个树种通过诱导2 n配子成功获得多倍体植株。人工诱导多倍体的方法包括生物学诱导、物理诱导和化学诱导3类，其中化学诱导类的秋水仙碱处理一直被认为是多倍体育种中最常用、最经济的手段。银杏（Ginkgo biloba L.）为银杏科（Ginkgoaceae)落叶乔木，是现存裸子植物门中最古老的子遗植物，我国重要的经济林树种，具有广泛的栽培价值。然而由于市场需求旺盛，目前已获得的银杏新品种资源仍然供不应求7 。研究表明，三倍体在生长速度、果实产量、抗逆性等方面发挥巨大作用8 因此培育银杏三倍体对于提高其综合利用价值、满足市场需求至关重要。秋水仙碱处理能够诱导银杏小孢子母细胞产生2 n雄配子，但2 n雄配子得率低，最高只能达到7%9，不足以用于银杏多倍体培育。然而，即使已有研究初步推测此现象可能与银杏小孢子母细胞中的蛋白质有关10 ,但秋水仙碱诱导银杏2 n雄配子低得率的具体原因和分子调控模式仍有待明确。鉴于此，本研究利用免疫荧光染色技术,结合考马斯亮蓝G-250法、三酮显色法、四氮唑蓝法、愈创木酚法和硫代巴比妥酸法等细胞生理生化检测技术，从细胞骨架和生理生化角度出发，连续跟踪观察和分析收稿日期：2 0 2 2-10-2 8 接受日期：2 0 2 3-0 2-0 7基金项目：国家自然科学基金项目（3197 16 7 5），国家重点研发课题（2 0 2 1YFD2200302），中央高校基本科研业务费专项资金资助（PTYX202236）作者简介：柯梦，主要从事经济林木良种繁育研究。E-mail：k e m e n g 12 16 16 3.c o m*通讯作者：孙宇涵，副教授，主要从事林木遗传育种研究。E-mail：s y h 8 310 0 8 16 3.c o m149837卷报核农学秋水仙碱诱导的银杏小孢子母细胞在各时间点的不同变化，最终经过筛选并确定可能影响银杏小孢子母细胞受秋水仙碱加倍诱导过程中的关键时期。然后通过转录组测序技术挖掘响应秋水仙碱诱导的相关基因，以期为解决秋水仙碱诱导银杏2 n雄配子低得率这一问题提供数据参考，同时为从分子水平探索秋水仙碱诱导生殖细胞加倍过程中的调控机制，从而进一步提高林木有性多倍化个体得率奠定理论基础1材料与方法1.1试验材料以北京林业大学鹫峰林场（40.0 6 7 8 7 N，116.0 8 134E）的一株30 年生银杏为试验材料，初春时取带雄球花的枝条于温室中水培（2 0 30），观察其生殖细胞发育状况。当其小孢子母细胞即将进人减数分裂的粗线期时10)，对雄球花进行0.6%秋水仙碱（colchicine）避光瓶浸法处理2 d(处理组，简称为CC)和蒸馏水处理2 d(对照组,简称为CK)1.2生理生化变化检测收集秋水仙碱和蒸馏水处理0、12、2 4、36、48、6 0、72、8 4、96 h 后的小孢子囊，-8 0 冷冻保存后用于生化指标测定，每项指标重复3次。其中，采用考马斯亮蓝G-250法测定可溶性蛋白含量12 ,节三酮显色法测定游离脯氨酸(proline，PRO)含量13),四氮唑蓝法测定超氧化物歧化酶(superoxidedismutase，SO D)活性13,硫代巴比妥酸法测定丙二醛（malondialdehyde，M D A）含量13，愈创木酚法测定过氧化物酶（peroxidase，POD)活性13 1.3碳蜡切片定期随机摘取处理组和对照组发育良好的3 5个银杏花芽，用福尔马林-乙酸-乙醇（formalin-aceto-alcohol,FAA)固定液V（乙醇）:V（甲醛）：V（冰醋酸)=18:1:1固定6 0 min后，使用聚乙烯二醇（polyethleneglycol,PEG)包埋银杏样本,切片观察141.4免疫荧光观察选取秋水仙碱和蒸馏水处理0、12、2 4、48 和96 h后的银杏小孢子囊，固定剂处理6 0 min，用双抗免疫荧光一抗：anti-tubulin（T-90 2 6,Si g ma 公司，美国），二抗:anti-mouse IgG FITC conjugated(F-0257,Sigma公司，美国）染色，在激光共聚焦扫描显微镜下观察小孢子微管蛋白的变化情况141.5RNA-seq测序及差异表达基因分析根据免疫荧光观察和生化检测结果，选择不同时间点（0、2 4、48、96 h)对照组和处理组的材料作为测序样本，3次生物学重复。首先利用Sequel测序仪第三代（PacBio公司，美国对样本全长cDNA进行测序，将获得的全长转录本序列与已发表的银杏全基因组测序结果15 比对,并将三代转录组注释结果和原基因组结果进行合并，作为进行二代转录组检测差异表达基因分析的参考序列比对注释，然后利用二代高通量测序技术Illumina HiSeqT(Illumina公司，美国）对随机抽取的处理组和对照组各3份重复样本材料进行文库构建和测序分析，最后经过滤得到cleanreads比对到前期合并后所得参考序列，经过差异表达分析筛选出样品间差异表达基因，并利用基因本体富集分析（GeneOntology，G O)和京都基因和基因组百科全书富集分析（KyotoEncyclopedia ofGenes and Genomes,KEGG)对筛选的差异基因进行进一步的功能富集，从而确定其在秋水仙碱加倍诱导下所行使的主要生物学功能。1.6数据处理采用GraphPadPrism9.0进行数据统计和生理生化指标变化分析，并应用SPSS25.0软件的T检验进行统计学检验，统计的显著性水平为P0.05。2结果与分析2.1银杏小孢子母细胞生理生化变化差异以蒸馏水处理为对照，选择小孢子母细胞发育的0、0.5、1、3、6、12、2 4、36、48、6 0、7 2、8 4和96 h进行生理生化检测。结果发现，对照组和处理组的生化指标随时间增加均有不同程度的变化。与对照组相比,除12h以外，48 h前处理组可溶性蛋白含量始终高于对照组，而48 h后，处理组可溶性蛋白含量呈下降趋势且低于对照组（图1-A）。相较于可溶性蛋白，对照组和处理组中游离脯氨酸含量变化趋势均保持一致，且处理组始终高于对照组，但在36 48 h对照组含量持续下降时，处理组游离脯氨酸含量开始出现缓慢上升，并在48 h下降（图1-B）。处理组的超氧化物歧化酶活性在0 12 h呈下降趋势，对照组则缓慢上升，延至12 h后，两种处理方式的超氧化物歧化酶活性均在极速上升(12 2 4h)后趋于平缓，且二者间差异较小(图1-C)。另外,对照组和处理组的过氧化物酶活性在0 48 h均表现为下降趋势，但处理间差异较小，48 h后，处理组的过氧化物酶活性开始缓慢上升且始终高于变化不明显的对照组（图1-D）。秋水仙碱诱导的丙二醛含量在096h均比对照组高，且在6 12 h间，处理组的丙二醛1499银杏小孢子母细胞响应秋水仙碱低频多倍化的细胞学和转录组分析8期CKCCAB600300*550602250500200(-1.8)/450150400*水350100300502500012243648 6072849601224364860728496CD1400400*350200300000250*800200600150400*100eaod2005000122436486072849601224364860728496E处理时间Time/h231898916151401224364860728496处理时间Time/h注：*、*和*分别代表在相同时间下，CK和CC之间在P0.05、P 0.0 1和P0.001水平差异显著。Note:*,*and*represent the significant difference between CK and CC at the same time at 0.05.0.01 and 0.001 level,respectively.图1银杏小孢子母细胞内不同生理生化指标随秋水仙碱处理时间的变化情况Fig.1Changes of different physiological and biochemical indexes in G.biloba microspore mother cellstreated with colchicine at different time含量整体呈增加趋势，并在36 48 h后呈稳定的上升趋势(图1-E)。上述结果表明，随时间增加,经秋水仙碱诱导的银杏小孢子母细胞生理生化指标表现出不同程度的差异，尤其在处理后12、2 4、48 和96 h差异较明显。2.2秋水仙碱对银杏小孢子母细胞微管骨架的影响2.2.1银杏小孢子母细胞微管骨架结构对秋水仙碱的响应用荧光显微镜分别观察处理组和对照组0、12、2 4、48 和96 h的小孢子母细胞。结果发现，与对照组相比，经秋水仙碱处理的小孢子母细胞微管结构在012h变化不明显（图2-A、B、F、G）但在2 4h,处理组微管结构开始出现减弱现象（图2-H），48 h 时，指示微管骨架的荧光强度减到最弱，甚至消失（图2-1)，而对照组的微管结构在2 4 48 h无明显变化（图2-C、D）。值得注意的是，至96 h时，对照组中小孢子母细胞正常进行减数分裂，处理组的绿色荧光也在逐渐显像，但依旧呈弥散状分布在细胞核内（图2-E、G)。上述结果表明，48 h可能是秋水仙碱诱导银杏小孢子母细胞多倍化形成的关键时间点。2.2.2银杏小孢子母细胞微管蛋白分布和微管厚度差异通过对秋水仙碱处理48 h的银杏小孢子母细胞微管进行高倍显微三维扫描发现，对照组银杏小孢子母细胞内的微管结构呈聚合状（图3-A）,而处理组小孢子母细胞的微管结构则呈弥散状（图3-B）。分析其微管骨架厚度发现，48 h对照组微管结构完整，呈聚合状，由蓝到红，微管蛋白层次由少到多（图3-C），处理组则呈弥散状，微管骨架厚度较薄（图3-D）。上述结果表明，秋水仙碱处理银杏小孢子母细胞48 h时，其微管蛋白的分布和厚度均出现了明显变化,因此可以进一步推测48 h可能是秋水仙碱诱导银杏小孢子母细胞多倍化的关键节点。2.3转录组测序挖掘与银杏小孢子母细胞多倍化相关的差异基因选取秋水仙碱处理组和对照组0、2 4、48 和96 h的小孢子母细胞进行转录组分析，相关性分析表明样品间相关性较好（图4）。转录组测序的Q30（质量值30)碱基百分比至少为8 8%，表明测序结果可靠，可用于后150037卷报核农学注：AE分别表示空白对照0、12、2 4、48、96 h。F J分别表示秋水仙碱诱导处理0、12、2 4、48、96 h。红色荧光表示细胞核内的染色质或染色体，绿色荧光表示微管蛋白。其中A、F标尺为10 0 m，B、G 标尺为2 5m，C、D、H 标尺为5m，E、I标尺为7.5m，J标尺为10 m。Note:A-E indicate the negative control at 0,12,24,48 and 96 h,respectively.F-J indicate the treated group at 0,12,24,48 and 96 h,respectively.Red fluorescence means chromatin or chromosomes in the nucleus and green fluorescence represents microtubule.The scare bar of A,F was 100 m,B,Gwas25m,C,D,Hwas5m,E,Iwas7.5m,Jwas10m.图2 秋水仙碱处理银杏小孢子母细胞细胞骨架显微结构Fig.2Cytoskeletal microstructure of G.biloba microspore mother cells treated with colchicine after different time3注：A、B分别代表对照组和处理组的微管结构;C、D 分别代表秋水仙碱处理组微管骨架厚度。Note:A,B represent the microtubule structure of the control and thetreatment group,respectively.C,D indicates thethickness of the microtubule skeleton,respectively.图348 h银杏小孢子母细胞微管高倍显微扫描结果Fig.3 High magnification microscopy of microtubules aftercolchicine treatment of G.biloba microspore mother cells for 48 h续分析。另外，在GO、K EG G 等数据库中分别注释到的基因数目为16 8 0 7、110 6 3，注释率达56.36%。以矫正后的P值错误发现率（falsediscoveryrate，FDR)0.05,差异倍数(fold-change,FC)2为筛选条件，鉴定到46 90 个差异表达基因（图5）。经统计分析发现：与0 h相比，蒸馏水处理银杏小孢子母细胞2 4、48和96 h后，差异表达基因的数量变化不明显；但相比之下，秋水仙碱诱导的差异表达基因个数在2 4、48和96 h分别增加了2 2.1%、8 3.9%、6 0.6%，且随着秋水仙碱处理时间延长，差异表达的基因数目也随之增多,96 h筛选到的差异表达基因数目（92 0 1）比2 4h1.00CK0-1CK0-20.80CK0-30.60CK24-1CK24-20.40CK24-3CC24-10.20CC24-20.00CC24-3CK48-1CK48-2CK48-3CC48-1CC48-2CC48-3CK96-1CK96-2CK96-3CC96-1CC96-2CC96-34484888图4样品间皮尔逊相关性检验Fig.4Pearson correlation among different examples(5144)多了7 8.9%（图5-A）。分析同一时间发现，处理24h时,共筛选到了193个差异表达基因，其中上调基因135个，下调基因58 个；至48 h时，鉴定到上调基因684个，下调基因7 41个，共142 5个差异表达基因；到96h时，挖掘到的差异基因明显增多,共340 0 个,其中2165个上调基因，12 35个下调基因（图5-B）。通过Venn图分析发现，在2 4、48 和96 h存在2 3个重叠差异表达基因（图5-C)。2.4差异表达基因功能注释将所筛选的差异表达基因注释到GO数据库中，并对不同时间的前2 0 个GO富集结果进行统计分析。结果显示，随着秋水仙碱处理时间延长，差异表达基因富集的数量明显增多（图6）。具体来看，秋水仙碱处1501银杏小孢子母细胞响应秋水仙碱低频多倍化的细胞学和转录组分析8期A10000上调下调总计9201800075756000573056365144421441194095413240003.48035653121266126852.023200015535981.4340CK24vsCK0CC24vs CK0CK48vsCK0CC48vsCKOCK96vsCK0CC96vsCKOB35003400C口上调下调总计CC96 vs.CK9630002500314521652.0001500142523209123523100068474197501143500135193580CC24vsCK24CC48vsCK48BCC96vsCK96CC24vsCK24CC48vsCK48图5差差异表达基因统计图Fig.5Statistics of number of differentially expressed geneACC24vsCK24BCC48vsCK48CCC96vsCK96转录，DNA模版小分子结合翻译转录调控，DNA模版应激反应核糖体初级代谢过程的调控嘌岭核糖核苷酸结合有机物生物合成过程含核碱机基化合物代谢过程的调控嘌岭核糖核苷三磷酸结合细胞器核酸模板转录的调控嘌岭核糖核苷结合非膜结合细胞器Q值氮化合物代谢过程的调控1.00嘌岭核苷结合Q值大分子生物合成过程Q值氮化合物代谢过程的调控0.75蛋白磷酸化1.00细胞内成分1.000.75大分子代谢过程的调控蛋白质修饰过程细胞内细胞器0.750.50磷酸化0.500.50大分子生物合成过程的调控0.25细胞内非-膜结合细胞器0.250.25基因表达的调控0.00细胞或亚细胞组分的运动细胞内0.000.00细胞代谢过程的调控基因数量基于微管的运动基因数量基因表达基因数量细胞大分子生物合成过程的调控25微管结合50细胞质部分40.50细胞生物合成过程的调控细胞的应激反应100细胞质8075生物合成过程的调控150120Q100细胞蛋白质修饰过程细胞组分160RNA代谢过程的调控125细胞交流细胞蛋白质代谢过程RNA生物合成过程的调控催化活性细胞大分子生物合成过程含核碱基化合物的生物合成过程阴离子结合细胞生物合成过程核酸模板化转录腺苷酸核糖核苷酸结合细胞部分芳香族化合物生物合成过程腺苷酸核苷酸结合细胞RNA生物合成过程ATP结合生物合成过程0.0350.0400.0450.0500.10.20.30.40.20.3富集因子富集因子富集因子EnrichmentfactorsEnrichmentfactorsEnrichmentfactors图6 差异表达基因GO富集前2 0 条目Fig.6Top 20 terms for GO enrichment analysis of differentially expressed genes理2 4h时，富集到的差异基因主要参与DNA模板转录、细胞内大分子的生物合成、RNA代谢、RNA生物合成、含碱基化合物代谢、含碱基化合物生物合成等生物过程（图6-A）。48 h 时，虽然有少数基因富集到“催化活性”这一分子功能条目上，但是大多还是参与核苷酸结合相关过程、蛋白修饰、ATP结合等生物过程。值得注意的是，48 h富集到与细胞或亚细胞成分运动、微管运动、微管结合等生物过程相关基因，且大分子生物合成过程在前2 0 个条目中没有富集（图6-B）。而9 6 h时，大部分差异基因富集到细胞内成分、细胞器、细胞质等细胞成分一类上,却没有与微管相关的条目,并且细胞大分子生物合成、细胞生物合成、大分子生物合成等生物过程又重新在96 h富集（图6-C)。对筛选的差异基因进行KEGG富集分析,并统计前2 0 个代谢途径。结果表明，仅“植物-病原菌相互作用”在3个时间都有富集；植物昼夜节律、酮化物生物150237卷报核农学合成、类黄酮生物合成、戊糖和葡萄糖醛酸相互转化在24、48 h 的前2 0 条通路中存在；另外,2 4和96 h中都注释到蛋白酶体、GTP结合蛋白质通路和植物激素信号转导途径相关基因；类化合物生物合成、谷胱甘肽代谢、转运蛋白通路相关基因在48 和96 h前2 0 条通路中均有富集。而内质网蛋白质加工相关的基因只在48 h时富集（图7）。2.5银杏小孢子母细胞中响应秋水仙碱诱导的转录因子对各个时间段的差异表达基因进行统计分析，发现在2 4、48 和96 h分别鉴定到了12、41、54个，共2 1类转录因子，其中ERF、D O F、b H LH 和WRKY转录因子家族广泛存在于各个时间点（图8）。具体来看，在秋水仙碱处理2 4h后，发现有4类转录因子家族共12 个基因，包括ERF（9个）、DOF（1个）、bHLH（1个）以及WRKY(1个）参与表达调控；而在48 h，除ERF5个）、bHLH(9个）、DOF(1个）、WRKY(3个）外,MYB、A P2、HD-ZIP、NA C、B3、H SF、NF-YC、RA V、T r i h e l i x、WO X、YABBY等13类转录因子也高差异表达；在96 h,ERF、bHLH转录因子数目明显增加，分别为2 0 和11个,而MYB、H SF、A P2、H D-ZIP、H SF、T r i h e l i x 等转录因子依然参与调控,NF-YB、b ZIP、M IK C_M A D S、VO Z等3类转录因子仅高差异表达（图8-A）。为了准确锚定到小孢子母细胞中响应秋水仙碱诱导的关键转录因子，用log（FC)值对48 h鉴定到的41个转录因子进行了热图分析。结果如图8-B所示,处理组中共有2 3个转录因A04712植物昼夜节律01008聚酮类生物合成蛋白04626植物病原菌互作04075植物激素信号转导00941类黄酮生物合成04016MAPK信号通路-植物00908玉米素生物合成01001蛋白激酶00520氨基糖和核苷酸糖代谢00360苯丙酸代谢00970氨酰_IRNA生物合成B.09131膜转运02010ABC转运蛋白03009核糖体生物合成00600鞘脂代谢00350酪氨酸代谢03050蛋白酶侠04031GTP奇蛋白04090个CD分子00040戊糖和葡萄糖醛酸的相互转化0.000.501.001.502.002.503.003.504004.50B00941类黄酮生物合成00940苯丙素生物合成01008聚酮类生物合成蛋白00130醒等酷类醒类尘物各成04712套夜节律-植物00480谷胱甘肽代谢B.09132信转号03110分子伴信和折叠催化剂00040戊糖和葡萄糖醛酸的相互转化00250丙氨酸、天冬氢酸和容氨酸代谢04626植物病原菌互作00052半乳糖代谢03032DNA复制蛋白04141内质网蛋白质加工02000载体B09106基他氢基酸的代谢03000转录因子B09101碳水化合物代谢B09183蛋百家族：信号和细胞过程01003糖基转移酶0.001.002.003.004.005.006.00C03010核糖体00900粘类骨架生物合成00480谷胱食肽代谢04016MAPK信号通路一植物03050蛋自酶体00270半胱氨酸和甲硫氨酸代谢B09106其他氨基酸的代谢04031个GTP结合蛋白04626植物-病原菌互作02044分泌系统B09123折叠，分类和降解00537糖基磷脂酰肌醇（GPI)锚定蛋白00650宁凌代谢04070磷脂酰肌醇信号系统02000载体04120泛素介导的善自水解04075植物邀素贺转笑00254黄苗雀莓素生物合成00190氧化磷酸花03020RNA聚合酶0.002.004.006.008.0010.0012.0014.0016.00-loglo(Pvalue)图7差异表达基因KEGG富集前2 0 通路Fig.7Top 20 pathways of KEGG enrichment analysis ofdifferentially expressed genesbHLH308.00ABMYBP-26.00DOF2.44.00MYBP-3BHLH352.00ATHB-10.00LHY-2.00MYB3-4.00ERF(9)bHLH(9)TGT2.-6.00GAMIDOF(I)MYB(6)WOX8bHLH(I)ERF(5)HOX32WRKY(1)6HLH8!CC24vSWRKY(3)6HLH49ERF(20)AP2(2)MYBP-1CK24bHLH(I1)HD-ZIP(2)DHLH113MYB(3)其他(2)YABS6HLH169NF-YB(3)MYB-related(2)LHY其他（3)INAC(2)ERF79bZIP(2)RAVICC96VSCK96B3(1)WRKY15CC48VsCK48HD-ZIP(2)HSF(O)AILHSF(2)NF-YC(1)ODO1PAR2Trihelix(2)RAV(I)NAMB2AP2(1)Trihelix(1)NAC25DOF()WOX(1)ERF18NFYCIMIKC_MADS(I)YABBY(I)EAMAMYB-related(1)DOF(1)ERFT6VOZ(I)HSFA2eMYBP-4WRKY(I)RA212WRKY49bHLH1i2AB13SWRKY35TINYBBM2CC48vsCK48注：括号中的数字代表差异转录因子个数。Note:The numbers in brackets indicate the number of transcription factors.图8各时间点差异表达转录因子类别（A）和48 h差异表达转录因子热图（B）分析Fig.8Classification of transcription factor families in different treated time(A)and the heatmap ofdifferentially expressed transcription factors(B)1503银杏小孢子母细胞响应秋水仙碱低频多倍化的细胞学和转录组分析8期子上调表达,综合FPKM(fragments per kilobase of exonmodel per million mapped fragments）值（5），其中log（FC)值大于2 的转录因子有3个（BBM2、ERF18 和RA212），大于1的有12 个，分别是AP2/ERF转录因子家族的AIL1、ERF16、ERF7 9,NA C 转录因子家族的NAC25和NAMB2,MYB-related中的LHY-2和ODO1,NF-YC中的NFYC1,bHLH中的PAR2,WRKY家族的WRKY15,Trihelix家族中的RAV1以及热胁迫响应蛋白家族HSF中的HSFA2e。这些转录因子可能为与秋水仙碱诱导银杏小孢子母细胞相关的重要候选基因。2.6银杏小孢子母细胞中持续响应秋水仙碱的差异表达基因在筛选到的2 3个差异表达基因中，仅4个基因得到功能注释，分别为热应激蛋白HSP11、转录因子DOF2.4、跨膜运动相关蛋白STC和磷饥饿响应相关蛋白PHL1。其中,HSP11在秋水仙碱处理48 和96 h表达量较高，FPKM值分别为110.32 45和98.90 56;DOF2.4在对照组中的表达量始终高于处理组，2 4h两者相差10倍，最高96 h表达量相差59倍；而STC只在96 h高表达，且CK96的表达量（6 4.147 0 FPKM）高于CC96(22.3137FPKM）;PH L1在0 h的CK中表达量最高(0.2872FPKM），然后随着秋水仙碱处理时间的延长而逐渐降低（图9)。另外，两个未知基因（unknown11和unknown19)在48 h时高差异表达，可能是该过程中关unknownll2.50unknown62.00unknown71.50unknown51.00unknown150.50unknown120.00unknownlo-0.50unknown17-1.00unknown14-1.50unknown8-2.00unknowngunknownlunknown2HSP11unknown16unknown13DOF2.4unknown19STCunknown4unknown18unknown3PHLICK24图92 3个重叠差异表达基因热图分析Fig.9Heatmap analysis of 23 overlapped differentiallyexpressed gene键的候选基因。总体来看，筛选到的2 3个差异表达基因并未全部注释到Swissport数据库中，可能由于所参考的基因组信息量较少。3讨论秋水仙碱作为抗有丝分裂药物，其浓度、处理开始时间和持续时间均会影响染色体的加倍16 。微管是几乎所有细胞中都存在的细胞器，在形态发生、细胞运动和信号传导等方面发挥重要作用，微管结构的变化还会影响细胞周期的正常进行。微管的结构亚基是100kDa的微管蛋白17 ,微管蛋白通过参与纺锤体、成膜体的形成影响染色体分离和胞质分裂，从而参与真核细胞的细胞周期18 。微管蛋白的聚合-解聚受温度、化学药物等影响19-2 1。用秋水仙碱诱导2 n花粉时，花粉母细胞减数分裂和胞质分离异常与2 n花粉的产生密切相关，微管蛋白则调控减数分裂和胞质分离的正常进行2 2 。本研究发现，在处理48 h时，秋水仙碱显著抑制了微管蛋白合成，严重影响了微管蛋白的聚合-解聚，迟滞减数分裂过程中染色体分离而进一步发育产生2 n花粉，这与李赞等2 3 对银白杨(P.alba）花粉诱导加倍的结果一致。表明秋水仙碱诱导银杏2n花粉的原理与其他树种类似,均是通过影响细胞内的微管蛋白合成，进而干扰染色体的正常分离活性氧（reactive oxygene species,ROS）和丙二醛（M D A)是影响生活细胞的重要因子2 4-2 5,当植物受到胁迫时，细胞内ROS和MDA积累，植物细胞会诱导抗氧化酶(POD、SO D 等)活性2 6 ,清除体内自由基。本研究生化指标结果显示，48 h时，处理组的总蛋白含量、SOD、PO D 和MDA含量都出现了明显变化，这与郭鹏2 7 的研究结果基本相似，说明48 h是银杏小孢子母细胞响应秋水仙碱的一个重要时间节点。秋水仙碱诱导生殖细胞多倍化困难可能是细胞受胁迫损伤后相关基因表达调控所致2 8 。本研究在48 h特异性地富集到了与微管、细胞或亚细胞运动相关下调表达基因的同时，也富集到了上调表达的内质网蛋白质加工相关基因，这表明在48 h时，秋水仙碱与微管蛋白二聚体结合可能抑制了微管的组装，但未减少微管的解聚，导致微管减少2 9,进而减少了花粉细胞内染色体的聚集，降低2 n花粉得率。进一步寻找影响上述过程的关键因子至关重要。本研究在秋水仙碱诱导48 h的花粉母细胞中鉴定到了41个差异表达的转录因子，其中最大的三个家族是150437卷报农学核AP2/ERF(7)、M YB-r e l a t e d(2)和 NAC(2)。A y a 等30 研究发现，编码AP2类型的转录因子SMOS1可引起细胞内微管取向异常。AP2还能参与微管乙酰化,影响细胞定向运动和趋化性31。此外，细胞骨架解体是花粉绒毡层凋亡的重要表现32 。研究发现，MYB80转录因子可以通过调节细胞程序性死亡时间调控绒毡层发育33。而细胞周期的研究结果表明，MYB-related转录因子家族成员Ceflp蛋白功能的缺失会导致-微管蛋白合成基因TUB1的mRNA低效剪接，导致细胞周期停滞在G2/M期34。这说明本研究筛选到的AP2/ERF、M YB-r e l a t e d 家族的转录因子可能在秋水仙碱诱导银杏小孢子母细胞的过程中发挥着关键作用。本研究中NAC转录因子和热应激蛋白HSPs在秋水仙碱处理过程中均出现特异性表达。NAC转录因子作为植物中最大的转录因子家族之一35,已被证明可能是花粉减数分裂和四分体形成的基因调控模块中的关键基因36 。研究发现，在水稻花粉中,Osj10gBTF3（一种NAC蛋白）可以通过与Hsp90家族的分子伴侣OsHSP82协调作用，影响下游蛋白质的折叠和组装，从而影响花粉发育37 。而以往研究也发现,在性器官中高表达的Hsf11,在胁迫条件下会负调控应激转录因子HsfA2e,从而调节植物对胁迫的响应38 。另外,HsfA2和Hsp17-CI在番茄花药发育过程中受到精细调控，并在热应激条件下得到进一步诱导39,同时，HsfA2的抑制也会降低花粉在减数分裂阶段受到胁迫时的生活力40 。本研究通过对比分析，在48 h的处理组中筛选到了高差异表达的热应激响应基因HSP11、热应激转录因子HSFA2e以及NAC转录因子家族中的NAC25和NAMB2,综合已有研究，推测这些基因可能在银杏小孢子母细胞响应秋水仙碱的基因调控模块中发挥着重要作用，但这一假设仍有待进一步研究证实。4结论本研究通过避光瓶浸法，用0.6%的秋水仙碱处理进人减数分裂粗线期的银杏小孢子母细胞，在处理的096h内对其微管骨架进行显微观察，并分析小孢子囊内可溶性蛋白、游离脯氨酸、丙二醛含量以及超氧化物酶和过氧化物酶活性等生理生化指标变化趋势，发现48h是秋水仙碱诱导银杏小孢子母细胞2 n雄配子的关键时间节点，并利用转录组测序技术筛选到了46 90 个差异表达基因,其中热应激响应基因HSP11在48 h处理组中表达量增加，另外通过对比分析，在48 h筛选到了2 3个上调表达的转录因子，其中最大的三个家族是AP2/ERF、M YB-r e la t e d 和NAC蛋白家族，这些均可作为后续进一步研究中的关键候选调控因子。参考文献：【1李玉岭，闫少波，毛秀红，王翠艳，王金娜，刘翠兰，张谦，乔艳辉.秋水仙素诱导林木多倍体研究进展J.农学学报，2 0 2 2，12(8):55-61 2 Wu S D,Han B C,Jiao Y N.Genetic contribution of paleopolyploidyto adaptive evolution in angiosperms J.Molecular Plant,2020,13(1):59-713Li u H L,Wa n g X B,Wa n g G B,Cu i P,Wu S G,A i C,H u N,LiA,HeB,Shao X J,WuZ C,Feng H,Chang YX,MuDS,HouJ,Dai X G,Yin T M