【www.zhangdahai.com--其他范文】
王彤彤 张留所
(1. 中国科学院海洋研究所 中国科学院实验海洋生物学重点实验室 山东青岛 266071; 2. 青岛海洋科学与技术试点国家实验室海洋生物学与生物技术功能实验室 山东青岛 266237; 3. 中国科学院海洋大科学研究中心 山东青岛 266071; 4. 中国科学院大学 北京 100049)
全球气候变化影响海洋动物的生长发育、繁殖和代谢等生命过程, 变化的海洋环境如何影响海洋动物发育的分子机制还不清楚。本研究团队构建了潜在的海洋线虫Litoditismarina模式动物研究体系和平台(Xieet al, 2020; Zhaoet al, 2021), 动物早期发育阶段对环境变化非常敏感(丛岩懿等, 2020; Xieet al, 2021),但海洋线虫早期发育的分子调控网络还不清楚。
秀丽线虫Caenorhabditiselegans的发育起始于受精卵, 经过多轮的细胞分裂和分化, 孵化为幼虫, 依次经历L1、L2、L3 和L4 共四个幼虫阶段发育成为具有多种组织和细胞类型的成虫(Zhanget al, 2015)。在20 °C 的条件下, 秀丽线虫从L1 幼虫发育到成虫约需3 d (Zhanget al, 2015)。Sulston 等科学家首次绘制出了秀丽线虫从受精卵到成体的完整细胞谱系, 在20 °C下, 从受精卵孵化为L1 约需16 h, 刚孵化的L1 幼虫含558 个细胞, 幼虫通过细胞分裂和分化并最终发育成具有959 个体细胞的成体, 不同个体间细胞发育谱系高度一致(Sulstonet al, 1977, 1983)。
研究发现约3 千万年前分化的两种同属线虫C.elegans和C.briggsae在发育过程中mRNA 和蛋白质的变化高度保守, 但同一物种内mRNA 和蛋白质的变化关联较弱, 该报道也对秀丽线虫胚胎期和L1 幼虫间的转录组变化特征进行了详细描述(Grünet al,2014)。另一项研究报道胚胎时期几乎一半的基因表达参与DNA 复制和染色质形成, 以实现快速的细胞分裂和细胞类型增加(Boecket al, 2016)。通过单细胞转录组测序, 已有报道量化了16 细胞期前基因表达变化, 发现胚胎期转录本多样性随着时间推移而增加(Tintoriet al, 2016)。通过全长转录组测序分析, 研究报道秀丽线虫胚胎转录本比幼虫阶段的转录本长,进一步揭示了秀丽线虫发育过程中转录组的动态变化(Liet al, 2020)。
我们团队聚焦的海洋线虫L.marina与陆生模式动物秀丽线虫同属线虫动物门小杆线虫科, 属于海-陆线虫近缘物种(Xieet al, 2020; Zhaoet al, 2021)。秀丽线虫发育机制研究已经比较深入, 而海洋线虫早期发育阶段的研究基本处于未知阶段。海洋线虫L.marina为雌雄异体, 其发育起始于受精卵, 受精卵在线虫子宫内完成早期胚胎发育, 产出体外的胚胎(eggs)孵化后为L1 幼虫, 20 °C 下, 从L1 幼虫发育为成体约需4～5 d(Xieet al, 2020, 2021)。L.marina从L1 幼虫发育为成体依次经过L2、L3 和L4 幼虫阶段, 从L1 幼虫发育为L2 约需46 h, 从L2 幼虫发育为L3 约需19 h, 从L3 幼虫发育为L4 约需16 h, 而从L4 幼虫发育为成体约需15 h (Zhaoet al, 2021)。研究人员通过绘制L.marina的胚胎发育细胞谱系, 发现L.marina与秀丽线虫胚胎细胞谱系同源性高达95.5%, 但末端细胞分化命运的相似度只有76.4%, 表明了细胞分化和发育命运决定机制在海陆线虫近缘种间既有相似性又有物种特异性(Houthoofdet al, 2003)。为了初步探究海洋线虫与秀丽线虫早期发育的异同, 本研究对L.marina胚胎时期和孵化后L1 发育阶段2 h、4 h 和6 h 的样品进行了转录组测序和分析, 研究结果表明, 海洋线虫与秀丽线虫早期发育过程中的调控机制高度相似, 秀丽线虫作为一种优秀的模式生物, 具有丰富的研究资源,这些资源将极大地辅助海洋线虫的基础生物学研究。同时进一步通过基因编辑和转基因营救等技术和方法深入比较研究L.marina和C.elegans发育调控关键基因的功能将为海-陆近缘线虫间的发育进化机制、海洋线虫对潮间带环境适应以及全球气候变化应答的分子机制研究提供新认知。
1.1 品系的获得与培养
本研究所用的海洋线虫L.marina野生型品系HQ1 采集自青岛汇泉湾潮间带的表层泥沙, 已在实验室连续培养5 年, 在20 °C 的条件下, 将其培养于90 mm SW-NGM 平板(Xieet al, 2020)上, 以大肠杆菌E.coliOP50 为食物。
1.2 海洋线虫大规模同步化
将海洋线虫在20 °C 下扩大培养到30 个板以上,待其生长至成虫且交配产生了大量的卵时, 将平板上的全部线虫、卵以及菌用灭菌海水洗脱至15 mL离心管中, 静置5 min 使大虫自然沉降, 卵及细菌在上清中, 将上清转移至新的15 mL 离心管中, 3 000g离心1 min, 弃上清。加入灭菌水, 混匀, 1 300g离心1 min, 弃上清, 重复三次。加入3 mL 灭菌水 + 3 mL裂解液[4 mL 高乐式®漂白水(≥6.25%的次氯酸钠溶液) + 1 mL 10 mol/L NaOH + 5 mL 灭菌水], 充分震荡75～90 s, 再加入灭菌水到12 mL, 1 300g离心, 弃上清, 重复3 次。此时得到大量的卵, 将其转移至海水中孵化18～20 h 得到大量同步化的幼虫L1。
1.3 RNAseq 样本制备
将同步化后孵化前得到的卵通过400 目网筛过滤,转移至1.5 mL EP 管, 8 000g离心, 弃上清, 置于液氮中速冻30 min 后保存在–80 °C, 得到胚胎样品。将同步化的L1 依次通过450 目和500 目网筛过滤后, 置于OP50 的SW-NGM 板上分别培养2 h、4 h、6 h 后, 将平板上的L1 用灭菌海水洗脱至15 mL 离心管, 600g离心1 min, 弃上清。加入灭菌水, 600g离心2.5 min, 弃上清, 重复3 次。将L1 转移至1.5 mL 离心管中, 加入灭菌水, 3 000g离心1 min, 弃上清, 置于液氮中速冻30 min 后保存至–80 °C, 得到L1 期2 h、4 h 和6 h 的样品。
1.4 RNAseq 建库及分析
提取样本 RNA, 并使用 NanoDrop 2000 测量RNA 浓度和纯度, 使用Agilent 生物分析仪2100 系统的RNA Nano 6000 检测试剂盒对RNA 完整性进行评估。使用Illumina 的NEBNext UltraTMRNA Library Prep Kit (NEB, USA)试剂盒生成测序文库。使用TruSeq PE Cluster Kit v4-cBot-HS (Illumia)在cBot Cluster Generation System 上对指标编码的样本进行聚类, 然后在Illumina 平台上对文库进行测序生成双端reads, 最终获得文库。
原始fastq 数据过滤掉含接头和低质量的reads后, 得到高质量Clean Data。使用HISAT2 (Kimet al,2015)将其比对到参考基因组得到Mapped reads, 然后利用String Tie (Perteaet al, 2015)将比对上的reads进行组装和定量。对各样品中的Mapped reads 的数目和转录本长度进行归一化, 采用FPKM (Floreaet al,2013)作为衡量基因表达水平的指标, 对每个样本分别进行基因表达水平定量。使用DESeq2 (Loveet al,2014)对两组样品进行差异表达分析, 将Fold Change≥1.5 且FDR≤0.05 作为差异表达基因筛选标准[差异倍数(Fold Change)表示两样品(组)间表达量的比值;错误发现率(False Discovery Rate, FDR)是使用Benjamini 和Hochberg 的方法通过对差异显著性P值(P-value)进行校正得到的]。将差异表达基因注释到KEGG (Kanehisaet al, 2008)数据库(http://www.genome.jp/kegg/), 使用KOBAS (Maoet al, 2005)软件来检测KEGG 通路中差异表达基因的统计富集, 筛选出两个样本间FDR≤0.05 的显著性差异通路。
1.5 实时荧光定量PCR 分析
本研究以EVM0015995_esyt-2为内参基因, 随机选择11 个关键基因(表1)进行qPCR 验证。使用ReverTra AceTMqPCR RT Master Mix with gDNA Remover 试剂盒将RNA 样品反转录为cDNA。cDNA中加入SYBR Green Realtime PCR Master Mix 试剂和正反向引物, 在QuantStudioTM6 Flex 荧光实时定量PCR 系统中检测, 每个样品设三个生物学重复。算出每个比较组合的ΔΔCt 后取平均, 2–ΔΔCt即为该基因在两个比较组合之间的倍数变化。使用Graphpad 对qPCR 和RNAseq 数据进行相关性分析和T 检验, 得到两组数据的皮尔逊相关系数和P值。
表1 qPCR 验证基因的引物序列Tab.1 The primer sequences of genes verified by qPCR
2.1 海洋线虫L. marina 转录组分析
选取L.marina胚胎期(Embryo, EM)和L1 幼虫2 h、4 h、6 h 四个阶段, 每个阶段设三个生物学重复, 使用二代Illumina 技术对其进行转录组测序。通过PCA主成分分析, 发现相同发育阶段的生物学重复高度相关, 不同发育阶段能够清晰区分(图1)。
图1 海洋线虫L. marina 早期发育转录组主成分分析(PCA)Fig.1 Principal component analysis of marine nematode L.marina in early developmental stages
通过分析不同样品间的差异表达基因数目, 发现L1 期2 h、4 h 和6 h 三个发育时期相较于EM, 均有超过9 000 个基因的表达发生了显著变化, 而L1 2 h、4 h 和6 h 间的差异基因数目较少(表2)。通过KEGG通路富集分析, 发现L1 幼虫2 h、4 h 和6 h 相较于EM 有超过30 个通路发生显著上调(图2), 而超过20个通路发生了显著下调(图3)。
图2 海洋线虫L. marina 各比较组合中上调基因的KEGG 富集情况Fig.2 KEGG enrichment of up-regulated genes in the marine nematode L. marina in each combination
图3 海洋线虫L. marina 各比较组合中下调基因的KEGG 富集情况Fig.3 KEGG enrichment of down-regulated genes in the marine nematode L. marina in each combination
表2 不同组间差异基因数目Tab.2 The number of differentially expressed genes between groups
2.2 相较于EM, L1 幼虫核糖体和核糖体发生相关基因的转录水平显著上调
通过KEGG 富集分析发现, 与胚胎期相比, L1 发育时期2 h、4 h 和 6 h 的多个核糖体蛋白大亚基编码基因(如EVM0009601_rpl-1和EVM0014298_rpl-3)和多个核糖体蛋白小亚基编码基因(EVM0000892_rps-1和EVM0013684_rps-7等)的表达均发生了显著上调(图4a), 核糖体发生通路中的52 个基因(如核仁蛋白EVM0010210_nol-6, Rio 激酶 EVM0005680_riok-1,核糖体RNA 加工蛋白质EVM0015730_rrp-8)在L1幼虫期发生了显著上调(图4b), 氨酰基-rRNA 生物合成通路中多个氨酰基-tRNA 合成酶(如丙氨酰基tRNA合成酶EVM0006554_aars-2, 天冬氨酸氨酰基tRNA合成酶EVM0015786_dars-1, 精氨酸氨酰基tRNA 合成酶EVM0008023_rars-1)在L1 幼虫期发生了显著上调(图4c)。
图4 核糖体相关通路基因表达变化情况Fig.4 Changes in gene expression of ribosomal related pathway
2.3 糖酵解/糖原异生、TCA 循环和氧化磷酸化通路中相关基因转录水平在L1 期显著上调
通过KEGG 富集分析发现, 与胚胎期相比, 糖酵解和糖原异生通路中的多个基因(如醛脱氢酶EVM0001434_alh-4和二氢硫辛酰胺脱氢酶EVM0004862_dld-1)在L1 发生显著上调(图5a)。柠檬酸(TCA)循环通路中的多个基因(如顺乌头酸酶EVM0011486_aco-2、柠檬酸合成酶 EVM0011427_cts-1)在L1 期发生显著上调(图5b)。氧化磷酸化通路中的多个基因(如ATP 合成酶EVM0008198_atp-1、液泡型H+ATP 酶EVM0002567_vha-2)在L1 发生显著上调(图5c)。
2.4 多个神经递质受体基因在L1 期转录水平显著上调
通过KEGG 分析, 发现L1 期相较于EM, 超过100 个基因在神经活性配体-受体相互作用通路中上调(图6)。EVM0000190_dop-1和 EVM0009467_dop-2等多个多巴胺受体相关基因在L1 期表达量上调(图6a),GABA β 受体相关基因如 EVM0013292_gbb-2,EVM0015297_gbb-1和GABA 受体基因EVM0014833_gab-1( 图 6b), 多个神经肽受体家族基因如EVM0011355_npr-29和EVM0012055_npr-34等(图6c),FMRFamide 肽受体家族基因EVM0006577_frpr-3、EVM0008649_frpr-14等也在L1 期表达量上调(图6d)。
图6 不同神经递质的基因表达变化情况Fig.6 Changes in gene expression of different neurotransmitters
2.5 DNA 复制和修复相关通路基因的表达在L1 期显著下调
通过KEGG 分析, 发现DNA 复制、错配修复、碱基切除修复和核苷酸切除修复通路在L1 期下调。DNA 复制通路中多个相关基因(如DNA 解旋酶/核酸内切酶EVM0005762_dna-2、DNA 连接酶EVM0006442_pri-1)。错配修复通路中多个相关基因(DNA 复制因子EVM0005944_rfc-1、复制蛋白 A EVM0007970_rpa-1)。碱基切除修复通路中多个相关基因(如DNA聚合酶EVM0014593_pole-1、无嘌呤/无嘧啶核糖内切酶EVM0009978_apn-1)和核苷酸切除修复通路中多个相关基因(如复制蛋白A EVM0002773_rpa-2、DNA 复制因子EVM0013177_rfc-4)在L1 时期显著下调(图7)。
图7 DNA 复制(a)、错配修复(b)、碱基切除修复(c)和核苷酸切除(d)修复通路中相关基因表达变化情况Fig.7 Changes in expression of related genes in DNA replication (a), mismatch repair (b), base exicision repair (c) and nucleotide excision repair (d)
2.6 Notch、Hippo 和Hedgehog 等重要发育相关信号通路基因在L1 期显著下调
Notch 信号通路中30 个基因(如EVM0011417_lag-1, EVM0004706_bir-1, EVM0008121_hda-1)的转录水平在L1 期显著下调(图8a), Hippo 信号通路中10多个基因(如EVM0004098_yap-1, EVM0016569_zyx-1)的转录水平在L1 期显著下调(图8b), Hedgehog 信号通路中10 个基因(如EVM0004427_ptc-3, EVM0015926_wrt-8)的转录水平在L1 期显著下调(图8c)。
图8 重要发育通路-Notch 信号通路(a)、Hippo 信号通路(b)和Hedgehog 信号通路(c)的相关基因表达变化情况Fig.8 Changes in gene expression of important developmental pathways-Notch signaling pathway (a), Hippo signaling pathway (b),and Hedgehog signaling pathway (c)
2.7 剪切体通路中多个剪切因子基因在L1 期显著下调
剪切因子对于RNA 的加工至关重要(Mortonet al,2011)。本研究中, 多个剪切因子(如BCAS 剪切因子EVM0013266_bcas-2、DBI1 剪切因子EVM0005825_dib-1、RSP 蛋白EVM0003635_rsp-6、PRP 剪切因子EVM0016298_prp-6、剪切因子3B 亚基 EVM0010128_sftb-1和 EVM0009161_sftb-2、SLU7 剪 切 因 子EVM0000329_sluh-7在L1 期下调(图9)。
图9 剪切体相关基因表达变化情况Fig.9 Changes in expression of spliceosomal related genes
2.8 qPCR 验证组学分析结果
本研究随机选择11 个基因进行qPCR 验证, 得出倍数变化(图10a)。采用皮尔逊相关系数r作为衡量转录组数据和qPCR 数据相关性的指标,r越接近1,相关性越好。通过Graphpad 分析, 得出转录组测序结果和qPCR 的结果相关性r为0.961 6, 显著性P值小于0.000 1 (图10b), 结果显示两者相关性很高。
图10 qPCR 验证结果Fig.10 Verification results through qPCR analysis
3.1 海洋线虫L. marina L1 幼虫发育需要更多的蛋白合成、能量供应与代谢活力
核糖体是蛋白质合成的场所, 与细胞生长和分裂密不可分(Fromont-Racineet al, 2003; Thomsonet al, 2013), 对秀丽线虫全基因组RNAi 筛查发现核糖体编码基因如rpl-4、rpl-35、rps-7、rps-13等基因的敲降会导致发育速率减慢(Sönnichsenet al, 2005)。研究发现RNAi 敲降rps-14导致秀丽线虫生长缓慢和不育(Chanet al, 2009)。通过GFP 标记秀丽线虫nog-1(编码核仁GTP 酶)启动子, 发现nog-1在胚胎早期到成虫均有表达, 使用RNAi 抑制nog-1的功能导致秀丽线虫孵化数减少、生长减慢、寿命增加和脂肪储存增加(Kimet al, 2014)。秀丽线虫riok-1编码RIO 激酶,RNAi 抑制riok-1的功能会导致性腺形态异常, 从而使卵数量减少(Weinberget al, 2014)。通过EMS 化学诱变和核糖体图谱分析, 发现秀丽线虫rrp-8(编码核糖体RNA 加工蛋白)参与核糖体RNA 加工, 对于核糖体的生物发生具有重要作用(Wuet al, 2018)。氨酰基-tRNA 生物合成通路中的氨酰基-tRNA 合成酶催化氨基酸与其同源tRNA 的共价连接, 在翻译过程中具有至关重要的作用(Rajendranet al, 2018)。rrt-1(rars-1)编码精氨酰tRNA 合成酶, 通过RNAi 抑制其功能发现,rrt-1对于秀丽线虫的发育、运动速度、缺氧敏感性、寿命有至关重要的作用(Fröhlichet al,2017)。本研究发现多个核糖体、核糖体发生和氨酰基-tRNA 生物合成通路相关基因在L1 时期的表达相较于胚胎期发生了显著上调(图 4), 意味着在L.marina中, 相较于胚胎时期, L1 期的发育可能需要更多蛋白质合成。据报道, 在秀丽线虫中, 相较于胚胎时期, L1 期的多个核糖体编码基因rpl-和rps-发生了显著上调(Grünet al, 2014; Boecket al, 2016; Liet al,2020), 这种变化模式与本研究海洋线虫L.marina中发现的一致。然而在秀丽线虫中, 相比于胚胎期, 核糖体发生通路多个基因在L1 期下调, 氨酰基tRNA合成通路中的基因在秀丽线虫EM 和L1 期之间变化不明显(Grünet al, 2014), 总之, 我们发现秀丽线虫这两个通路相关基因的表达变化特征和L.marina不同, 推测可能和两种线虫细胞发育过程存在的差异相关。
核糖体的生物发生是所有细胞中对能量要求最高的生命过程之一(Kressleret al, 2010; Thomsonet al,2013)。本研究中, L1 期相较于胚胎期, 糖酵解/糖原异生通路、TCA 循环通路和氧化磷酸化通路等与能量相关的通路中多个基因在L1 期显著上调(图5)。在糖酵解/糖原异生通路中,alh-4编码醛脱氢酶, 通过正向遗传筛选发现其缺失会增加脂肪醛, 减少脂肪的储存(Zenget al, 2021)。RNAi 抑制秀丽线虫中dld-1(编码二氢硫辛酰胺脱氢酶)的功能, 可改变能量代谢,激活与蛋白酶体降解相关的通路, 改善细胞间信号传导, 延长寿命(Ahmadet al, 2021)。TCA 循环中,cts-1编码的柠檬酸合成酶是TCA 循环的起始酶和限速酶, RNAi 抑制秀丽线虫中cts-1的功能, 会导致早期发育停滞(Hadaet al, 2019)。秀丽线虫aco-2编码乌头酸酶,idha-1编码异柠檬酸脱氢酶, RNAi 抑制两者的功能会导致柠檬酸的积累, 并延迟秀丽线虫的生长和发育(Yanget al, 2022)。氧化磷酸化通路中:vha-8和vha-19编码液泡型H+ATP 酶,vha-8是PH 动态平衡和幼虫发育所必需的,vha-19是胚胎发生所必需的,RNAi 抑制秀丽线虫vha-8或vha-19的功能会致死(Jiet al, 2006; Knightet al, 2012)。因此, 我们推测在L1期, 这些与能量相关的基因转录水平的增加, 为核糖体的生物合成提供了必要的能量, 同时也为L1 幼虫的发育提供了充足的能量供应。研究发现秀丽线虫从胚胎期发育到L1 期, 上述三个能量通路中的多个基因的表达均发生显著上调(Grünet al, 2014; Boecket al, 2016; Liet al, 2020)。因此, 我们推测L.marina早期发育的能量供应机制与秀丽线虫高度相似。
本研究通过KEGG 富集分析, 发现在L1 期, 17个代谢通路在L1 期显著上调, 相关的通路是碳代谢,嘌呤代谢, 花生四烯酸代谢, 亚油酸代谢, 淀粉和蔗糖代谢, 甘氨酸、丝氨酸、苏氨酸代谢, D-谷氨酰胺和D-谷氨酸代谢, 谷胱甘肽代谢, 2-氧代羧酸代谢,甘油磷脂代谢代谢, 组氨酸代谢, 亚麻酸代谢, β-丙氨酸代谢, 丙酸酯代谢, 氮代谢(图2)。4 个生物合成通路在L1 期显著上调, 相关的通路是氨基酸的生物合成, 精氨酸的生物合成, 叶酸的生物合成, 类固醇的生物合成(图2)。与海洋线虫相似, 这些通路中的大部分基因在秀丽线虫L1 期表达量也显著高于胚胎期(Grünet al, 2014; Boecket al, 2016; Liet al, 2020),据此推测L.marina早期发育的代谢调控机制与秀丽线虫高度保守。
3.2 海洋线虫L. marina L1 幼虫神经受体基因表达显著上调
通过KEGG 富集分析, 发现L1 期相较于EM, 超过100 个基因在神经活性配体-受体相互作用通路中上调, 大多是神经递质或神经肽受体, 可能与L1 幼虫的觅食和运动等行为调控相关。多巴胺受体和GABA 受体在L1 表达量较胚胎期显著上调(图6a,6b)。多巴胺和谷氨酸-氨基丁酸(γ-aminobuyric acid,GABA)对于线虫的基本运动功能至关重要(Jorgensen,2005; Omuraet al, 2012), 研究发现编码多巴胺受体的dop-2在多巴胺能神经元中表达, 可以调节多巴胺的合成和释放(Suoet al, 2004), 多巴胺能神经元在感觉到细菌食物时释放多巴胺从而引发运动减慢反应(Suoet al, 2004; Nagashimaet al, 2016)。研究发现多巴胺受体编码基因dop-3缺失会导致秀丽线虫运动缺陷(Chaseet al, 2004)。GABA 是一种重要的抑制性神经递质, 在线虫运动和觅食时, 起到放松身体肌肉的作用(Jorgensen, 2005)。研究发现通过 RNAi 抑制gbb-1(编码代谢型GABAβ 受体)的功能增强抗氧化能力并延长秀丽线虫寿命(Yuanet al, 2019)。
多个神经肽受体家族基因在L1 期表达量显著上升(图6c, 6d)。在秀丽线虫中, 目前已鉴定出113 个神经肽基因, 对于感觉和运动至关重要(Liet al, 2008)。比如, 通过RNAi 敲降实验发现, 秀丽线虫运动和感觉唤醒需要frpr-3(编码FMRFamide 肽受体) (Chewet al, 2018)。在食物充足的情况下,npr-1(编码神经肽受体)可以导致秀丽线虫在进食过程中聚集, 并在细菌边缘聚集(Liet al, 2008)。因此, 随着线虫孵化后开始运动, 这些神经递质表达量显著升高。除此之外, 速激肽受体家族TKR, 胆囊收缩素受体CKR, 神经调节肽NMUR, NMDA 类谷氨酸受体NMR, 代谢型谷氨酸受体家族MGL 等神经递质在L1 期表达量较EM均发生显著上调。在秀丽线虫中, 相比于胚胎期, 多种神经递质和神经肽受体基因也均在L1 期转录水平上调(Grünet al, 2014; Boecket al, 2016; Liet al,2020)。据此推测, 海洋线虫和秀丽线虫在L1 幼虫发育阶段神经受体调控方面可能具有类似的机制。
3.3 海洋线虫L. marina 胚胎期复制与修复和细胞分化调控相关基因表达显著上调
DNA 复制是细胞分裂周期的一个重要方面, 发生在细胞周期的合成阶段(S 期), 准确无误的DNA 复制对于保持基因组的完整性必不可少(Bellet al, 2002;Ariaset al, 2007)。在海洋线虫中DNA 复制、错配修复、碱基切除修复和核苷酸切除修复通路中的多个基因, 与胚胎期相比, 在L1 幼虫期显著下调(图7)。研究发现秀丽线虫RPA 蛋白(编码复制蛋白A)参与有丝分裂复制、减数分裂、同源重组和修复等多个过程,对于DNA 复制过程至关重要, CRISPR/Cas9 敲除rpa-1会导致秀丽线虫早期幼虫致死(Hefelet al,2021)。在酵母和哺乳动物中, RFC 蛋白(编码DNA 复制因子)参与将某些DNA 聚合酶加工因子加载到单链DNA 上, 调控DNA 复制过程(Venclovaset al, 2000;Griffithet al, 2002; Bermudezet al, 2003)。秀丽线虫中,6 个mcm-基因(mcm-2、mcm-3、mcm-4、mcm-5、mcm-6和mcm-7)编码的6 个MCM 蛋白(微小染色体维持蛋白)形成六聚体复合体, 作为DNA 复制前起始复合体的关键组分, 在DNA 复制中起到DNA 解旋酶的作用(Korzeliuset al, 2011)。RNAi 抑制秀丽线虫dna-2的功能会降低胚胎发育速率、缩短寿命, 其编码的DNA解旋酶/核酸内切酶可以维持基因组的稳定性, 在细胞分裂和DNA 复制中都具有非常重要的作用(Leeet al, 2011)。秀丽线虫div-1、pola-1、pri-1和pri-2编码的DNA 聚合酶α-引发酶复合物启动增殖细胞的DNA 复制, RNAi 抑制这些基因的功能会使得胚胎致死率升高, 此外, 增殖细胞命运的维持需要pola-1、pri-1、pri-2(Yoonet al, 2018)。我们推测, 与之相关基因在胚胎期的高表达可能反映了胚胎期较高的DNA 复制率和细胞分裂, 因此, 相较于L1, 胚胎期细胞分裂速度可能更快, 从而需要复制与修复相关基因的高表达。这些基因在秀丽线虫中, 相较于胚胎期, 其转录水平在L1 期也均显著下降(Rünet al, 2014;Boecket al, 2016; Liet al, 2020), 因此L.marina和秀丽线虫类似, 胚胎时期复制和修复等多个通路基因表达显著上调可能用以支持快速、精准的细胞分裂和增殖。
KEGG 富集分析发现海洋线虫在L1 期, Notch、Hippo 和Hedgehog 信号通路中多个基因, 与胚胎期相比, 在L1 期转录水平显著下调。lag-1编码Notch信号通路中的重要转录因子LAG-1, 在发育过程中介导细胞间的相互作用, RNAi 敲降lag-1会导致孵化后的L1 秀丽线虫发育停滞(Maicaset al, 2021)。Hippo信号通路中,yap-1(编码Yes 相关蛋白)参与秀丽线虫Wnt 介导的神经元极化(Leeet al, 2018), 也与秀丽线虫的耐热性和衰老有关(Iwasaet al, 2013)。通过秀丽线虫EMS 诱变发现,zyx-1(编码含LIM 结构域的黏着斑蛋白)调节秀丽线虫突触维持(Luoet al, 2014)。Hedgehog 信号通路中,ptc-3(编码PTC 蛋白)在秀丽线虫中的表达是高度动态的, 通过RNAi 抑制秀丽线虫ptc-3基因的功能会导致致死、不完全蜕皮、体型缩小、vulva 发育缺陷等(Solovievet al, 2011)。秀丽线虫wrt-8(编码Hedgehog 相关蛋白)的过度表达会导致轴突缺陷(Riveiroet al, 2017)。三个通路中在秀丽线虫中的大多数同源基因, 相较于胚胎期, L1 期的转录表达也均发生了显著下调(Grünet al, 2014; Boecket al, 2016; Liet al, 2020)。因此, 我们推测海洋线虫和秀丽线虫在胚胎期均需要更高的发育和细胞分化调控因子表达水平。
KEGG 富集分析表明, 相比于胚胎期, 海洋线虫在L1 期多个剪切体通路相关基因的转录水平显著下降。剪切因子对于RNA 的加工过程至关重要(Mortonet al, 2011)。研究表明prp-40调控秀丽线虫可变剪切,其编码的PRP-40 mRNA 前体加工因子是剪接体的重要组成之一(Choudharyet al, 2021)。通过RNAi 抑制prp-8的功能导致秀丽线虫剪切异常, 其对于维持整个胚胎发生过程中的细胞分裂周期调控起重要作用(Hebeisenet al, 2008)。sftb-1编码剪切因子3B 亚基,mCherry 标记秀丽线虫sftb-1发现其是一种普遍表达的发育必需基因(Serratet al, 2019)。在秀丽线虫中,这些剪切因子的基因L1 期表达量也均低于胚胎期(Grünet al, 2014; Boecket al, 2016; Liet al, 2020)。最近报道发现秀丽线虫胚胎转录本的平均长度比胚胎后幼虫阶段的较长(Liet al, 2020), 因此推测海洋线虫与秀丽线虫在胚胎发育的过程中, 需要更高的剪切因子相关基因表达, 可能与胚胎期较长的转录本特征维持和细胞的快速分裂和分化等相关。
我们的研究结果表明, 与胚胎期相比, 孵化后早期L.marinaL1 幼虫的多个通路如核糖体、核糖体生物发生、糖酵解/糖原异生、TCA 循环和氧化磷酸化通路相关基因发生了显著上调; 多个神经递质和神经肽受体基因如dop-、npr-和frpr-等在L1 期转录水平也显著上调。与胚胎期相比, L1 幼虫DNA 复制和修复相关、Notch、Hippo 和Hedgehog 信号通路和剪切体调控相关基因发生了显著下调。本研究发现的海洋线虫早期发育转录组变化模式与陆生模式生物秀丽线虫高度相似, 但同一上调或下调通路中具体发生表达变化的基因有些不同。另外,L.marinaL1 幼虫显著上调的核糖体生物发生通路相关基因在秀丽线虫中发生了显著下调。因此,L.marina与秀丽线虫早期发育的表达调控机制高度保守, 秀丽线虫丰富的生物信息学和突变体库等资源和成熟的分子生物学和遗传学等研究手段将极大地辅助海洋线虫的基础生物学研究, 并推动模式海洋线虫体系和平台构建。结合秀丽线虫成熟的研究方法, 在L.marina中进一步优化基因编辑等遗传操作将会为海-陆近缘线虫间的发育进化机制, 海洋线虫对潮间带环境适应以及海洋动物如何应答和适应全球气候变化的分子机制研究提供新认知, 也可为蜕皮类经济物种如虾蟹的遗传育种与产业发展提供重要理论和技术支持。
猜你喜欢 核糖体线虫幼虫 核糖体成熟因子RimP、Era和RimJ的研究进展内蒙古民族大学学报(自然科学版)(2022年2期)2022-11-22夏季蔬菜换茬期线虫防治要注意今日农业(2022年14期)2022-09-15核糖体生物合成与肿瘤的研究进展肿瘤防治研究(2020年5期)2020-07-09例析翻译过程中核糖体移动方向的判断生物学教学(2019年9期)2019-09-23地黄花对秀丽线虫寿命的影响中成药(2018年2期)2018-05-09飞舞吧，凤蝶金色少年(奇趣科普)(2017年4期)2017-06-05朝鲜孢囊线虫——浙江省孢囊线虫新记录种浙江大学学报（农业与生命科学版）(2017年1期)2017-04-17柿树夏季咋防柿蒂虫农村农业农民·B版(2016年7期)2016-10-21线虫共生菌Xenorhabdus budapestensis SN19次生代谢产物的分离纯化与结构鉴定天然产物研究与开发(2016年6期)2016-06-05白带锯蛱蝶幼虫和蛹的营养成分比较分析天然产物研究与开发(2014年8期)2014-04-27本文来源：http://www.zhangdahai.com/shiyongfanwen/qitafanwen/2023/0809/637082.html
