Effect of high temperature stress on intestinal tissues morphology and transcriptome of Procambarus clarkii
-
摘要:
为探究克氏原螯虾 (Procambarus clarkii) 肠道在接近其生理极限的高温胁迫下的作用机制,对其进行32 和37 ℃的高温胁迫实验,分别于胁迫24、72 h后收集肠道组织,进行组织切片观察和转录组学分析。结果表明,高温对克氏原螯虾的肠道组织结构造成了明显损伤,且随着胁迫温度和胁迫时间的增加,损伤加剧。与对照组 (26 ℃) 相比,在32 ℃胁迫24和72 h条件下,分别鉴定到2 462和4 619个差异基因;在37 ℃胁迫24和72 h条件下,分别鉴定到1 825和7 298个差异基因。KEGG富集分析结果显示,差异基因在内质网蛋白加工、代谢途径等通路中显著富集,其中在同一胁迫时间内胁迫温度从32 ℃升至37 ℃时,内质网蛋白加工通路富集频率升高,表明其在应对高温胁迫中发挥重要作用;GO功能注释显示,差异基因主要在碳水化合物代谢、蛋白质折叠等相关过程中显著富集。
Abstract:To explore the mechanism of the intestinal tract of red swamp crayfish (Procambarus clarkii) under high temperature stress close to its physiological limit, we conducted the high-temperature stress experiments at 32 ℃ and 37 ℃, and collected the intestinal tissues after 24 h and 72 h of stress, for tissue sectioning observation and transcriptomic analysis. The results show that high temperature caused significant damages to the intestinal tissue structure of crayfish, and the damage increased with the increase of stress temperature and stress time. Compared with the control group (26 ℃), 2 462 and 4 619 differentially expressed genes (DEGs) were identified during 32 ℃ stress for 24 h and 72 h, respectively. During 37 ℃ stress for 24 h and 72 h, 1 825 and 7 298 DEGs were identified, respectively. KEGG enrichment analysis shows that the DEGs were significantly enriched in the endoplasmic reticulum protein processing and metabolic pathways, and the enrichment frequency of endoplasmic reticulum protein processing pathways increased when the stress temperature increased from 32 ℃ to 37 ℃ during the same stress time, which indicates that they played an important role in coping with high temperature stress. According to the GO functional annotation, DEGs were mainly enriched in the energy metabolism-related processes such as carbohydrate metabolism and protein folding.
-
Keywords:
- Procambarus clarkii /
- High-temperature stress /
- Transcriptome /
- Intestine /
- Paraffin section
-
角鲨烯 (Squalene, C30H50) 是一种含有6个双键的开链线性三萜,主要存在于深海鲨鱼肝油中,少部分存在于细菌、真菌、微藻和植物中[1-2]。角鲨烯具有抗肿瘤、耐缺氧、抗氧化、抗辐射和细胞保护等生物活性[3],已在水产、医药和化妆品领域获得广泛应用[4-5]。在水产领域,提升水产动物体内的角鲨烯含量可明显增强其非特异性免疫细胞水平和抗病力。给罗非鱼 (Oreochromis spp.) 注射角鲨烯可提高其巨噬细胞增殖量,获得了更多的非特异性免疫细胞[6]。此外,角鲨烯被认为是优良的黏膜免疫佐剂而应用于鱼类口服疫苗开发中[7]。
目前,鲨鱼肝油仍然是角鲨烯的主要来源,但不同种类鲨鱼肝油中的角鲨烯含量差异显著[8],随着海洋污染的加剧和鲨鱼资源的过度捕捞,亟需寻找其他可替代的天然资源。微生物发酵生产角鲨烯是一条绿色可循环、周期短、可规模化生产的途径,是工业化生产角鲨烯的优良选择[9]。裂殖壶菌 (Schizochitrium sp.)、大肠杆菌 (Escherichia coli) 和酿酒酵母 (Saccharomyces cerevisiae) 作为经典的工业发酵微生物,通过基因工程手段提高角鲨烯产量已有报道[9-11],然而这些微生物自身的角鲨烯产量并不高,选育更具发酵潜力的微生物菌株[12-13],利用基因工程手段进一步提高角鲨烯产量将是一种更有潜力的手段。Pseudozyma sp. SD301是本项目组自主筛选到的单细胞海洋类酵母,摇瓶培养时角鲨烯产量达0.79 g·L−1,具有较高的角鲨烯产量[12]。碳离子 (12C6+) 束辐照诱变育种具有诱变率高、诱变谱宽、突变体稳定等优势,已广泛应用于微生物育种领域[14-16]。本研究以Pseudozyma sp. SD301为出发菌株,采用12C6+束辐照诱变技术,选育更加高产角鲨烯的突变株。在此基础上,以高产突变株为例,开发高效的电转化方法,用于实现外源基因的表达,为进一步遗传改良提供了可能,为角鲨烯的微生物发酵生产奠定基础,并为深海来源的角鲨烯资源提供可能的替代选择。
1. 材料与方法
1.1 菌株
以实验室自主选育的野生型类酵母Pseudozyma sp. SD301 (CGMCC no. 9687) 作为初始菌株,该野生菌株具有较高的产角鲨烯性能[12]。
1.2 主要试剂
酵母提取物购自英国Oxoid;过氧化氢 (H2O2) 和葡萄糖购自上海国药集团化学试剂有限公司;琼脂粉和潮霉素B购自北京索莱宝科技有限公司;角鲨烯购自北京百灵威科技有限公司;质粒提取试剂盒、DNA片段纯化试剂盒和胶回收试剂盒均购自OMEGA生物公司;限制性内切酶购自Thermo Scientific公司;2×ES Taq MasterMix购自康为世纪生物科技有限公司;2×PCR Besttaq MasterMix购自中国南京爱必梦生物材料有限公司;其余常用生化试剂均为国产分析纯。
1.3 培养基与培养条件
GY液体培养基 (1 L):葡萄糖60 g,酵母提取物20 g,海水晶15 g。
GY固体培养基 (1 L):葡萄糖60 g,酵母提取物20 g,海水晶15 g,琼脂粉20 g。
培养条件:菌株在25 ℃、200 r·min−1条件下,使用GY液体培养基活化2次,每次活化2 d (批次培养时,2 d达指数生长期,3 d后达平台期,此时胞内角鲨烯质量最高)[12],然后按照10%接种量转接到摇瓶中培养,在相同条件下继续发酵3 d,之后收集菌体用于菌体干质量及角鲨烯质量分析。
1.4 12C6+束辐照、突变株的筛选及稳定性分析
12C6+束辐照:以SD301为出发菌株,培养2 d后,将菌液分装成每管2 mL,用于辐照实验。12C6+束辐照实验在中国科学院近代物理研究所重离子加速器国家实验室 (Heavy Ion Research Facility in Lanzhou, HIRFL) 进行,辐照条件为:采用12C6+束,能量80 MeV·µ−1,LET 31 keV·µm−1,辐照剂量为0、20、60、90、120、150和180 Gy,每个辐照剂量设置3个平行样品。经过辐照后的菌液稀释到1×10−4和1×10−5,取100 μL菌液稀释液涂布在GY固体平板上,之后25 ℃倒置培养,用于计数不同辐照剂量对应的存活菌落数。
突变株筛选条件:为了确定产角鲨烯类酵母SD301对H2O2的敏感性,将类酵母SD301菌液稀释到1×10−5后,涂布在H2O2浓度分别为10、100和1 000 μmol·L−1的平板上培养并记录菌落数。
突变株筛选:选择120、150和180 Gy的12C6+束辐照剂量处理后的菌液进行适当稀释,进行H2O2压力筛选 (稀释后的菌液涂布在含有100 μmol·L−1 H2O2的GY平板上,培养2 d),挑取单菌落在25 ℃、200 r·min−1培养3 d,收集菌体提取角鲨烯,应用气相色谱分析角鲨烯质量,每个样品平行3次测定。
突变株稳定性:将获得的突变株在25 ℃、200 r·min−1条件下,使用50 mL GY液体培养基传代培养,培养时间为2 d,连续传代15次后,测定突变株的角鲨烯质量,每个样品平行3次测定。
1.5 突变株的角鲨烯产量和质量分析
将对照菌株SD301和突变株PS120、PS150和PS180分别接种到GY培养液中25 ℃、200 r·min−1活化2次,之后分别在20 mL培养液和100 mL培养液中批次培养3 d,11 000×g离心2 min收集一定量菌体进行冷冻干燥,干燥后称质量计算生物量。称取25 mg干菌粉,加入300 μL 6 mol·L−1盐酸 (HCl),在−80 ℃和100 ℃下反复冻融3次,冷却后加入450 μL 氯仿-甲醇溶液[V(氯仿)∶V(甲醇)=2∶1]、500 μL饱和氯化钠 (NaCl) 和500 μL正己烷,混匀,10 000×g离心5 min,取上清过有机滤膜后用于气相色谱分析角鲨烯质量。角鲨烯质量根据角鲨烯标准曲线计算得到。每个样品平行3次测定。角鲨烯产量根据角鲨烯质量及生物量计算得到。
1.6 质粒pYES2-hph和pYES2-hph-egfp的构建
以潮霉素B作为遗传筛选标记,潮霉素B的质量浓度为300 mg·L−1。以pYES2 (Invitrogen)为初始质粒,以HPH (潮霉素B抗性基因) 的DNA为模板,用引物Hph-F/Hph-R扩增HPH基因,用BamHI和XhoI双酶切的方式连接到载体pYES2,构建质粒pYES2-hph,用引物P-C-F/P-C-R进行菌落PCR验证。以EGFP (增强型GFP基因) 的DNA为模板,用引物P-egfp-F/P-egfp-R扩增EGFP基因。用XhoI单酶切质粒进行线性化,将EGFP基因的PCR产物和线性化载体pYES2-hph (XhoI) 在Exnase的催化下进行重组反应,之后将反应产物加入大肠杆菌Top10的感受态细胞中进行转化,转化完成后,将细胞液涂布到含有100 mg·L−1氨苄青霉素的LB平板上,37 ℃倒置过夜培养。用引物P-egfp-C-F/P-C-R对平板上长出的克隆进行菌落PCR验证,在PCR验证目的片段大小正确后,使用OMEGA公司的胶回收试剂盒进行PCR产物回收并送到北京擎科生物科技有限公司进行测序。本研究中质粒构建所用引物见表1。
表 1 本研究中质粒构建所用引物Table 1. Oligonucleotides used in vector construction引物
Primer序列 (5'—3')
Sequence (5'—3')Hph-F CGGGATCCATGAAAAAGCCTGAACTCACCGC Hph-R CCGCTCGAGTTCCTTTGCCCTCGGACGAG P-C-F GGTCGCGTTCCTGAAACGCAG P-C-R CAAGCAAGGTTTTCAGTATAAT P-hph-C-F ATATGTCCTGCGGGTAAATAGC P-hph-C-R ATGCCTCCGCTCGAAGTAGCGCGTC P-egfp-F ACTCGTCCGAGGGCAAAGGAAATGGTGAGCAAGGGCGAGGAG P-egfp-R CCCTCTAGATGCATGCTCGAGCTTGTACAGCTCGTCCATGCC P-egfp-C-F ATCGTCCGATCCGGAGCCGGGACTGT 1.7 电转化方法
已有研究证实使用醋酸锂和山梨醇处理Pseudozyma属细胞制备的感受态细胞有利于提高其电转化效率[17-18],本研究中使用的电转化方法在Marchand等[18]、Masaaki等[19]和Morita等[20]报道的方法基础上进行了改进和电压优化。
YC预处理液 (1 L):山梨醇1 mol·L−1,Tris-HCl 10 mmol·L−1,LiAc 100 mmol·L−1,DTT 10 mmol·L−1,pH 6.4。
CX细胞重悬液 (1 L):1 mol·L−1山梨醇-GY液体培养基 (V∶V=1∶1)。
YP液体培养基 (1 L):蔗糖60 g,酵母提取物20 g,海水晶15 g,甘油100 g。
GP固体培养基 (1 L):半乳糖10 g,蔗糖60 g,酵母提取物20 g,海水晶15 g,琼脂粉20 g,甘油100 g,潮霉素B 300 mg。
感受态细胞制备:菌株在25 ℃、200 r·min−1条件下,使用YP液体培养基活化2次,然后取出1 mL菌液稀释10倍后测定其在600 nm的光密度 (OD),OD600在1左右即可用于感受态细胞制备。取1 mL菌液,在4 ℃、4 000×g离心2 min收集菌体,菌体用1 mL预冷的无菌蒸馏水洗涤2次,之后用1 mL预冷的YC预处理液漂洗2次,再用1 mL预冷的CX细胞重悬液洗涤2次,最后用预冷的CX细胞重悬液将菌液稀释250倍,使得细胞个数约为1×108个·mL−1。
电击条件:将50 μL感受态细胞液和5 μg质粒混匀后加入预冷的0.2 cm电转杯 (Bio-Rad),冰上静置15 min,之后置于自主研发的细胞电穿孔高压脉冲发生器 (SW11-3K) 中电击。参数设置:电压650、900、1 150和1 400 V,频率2 000 Hz,占空比10%和脉冲数40个,电击后立即加入1 mL的CX细胞重悬液到电转杯中,并将其转至1.5 mL试管,冰上放置15 min,之后于25 ℃、200 r·min−1复苏培养1 h,之后取出100 μL涂布GP固体培养板,之后25 ℃倒置培养。
转化子验证:挑取GP固体培养板上长出的单菌落,使用引物P-C-F和P-C-R对pYES2-hph-egfp的转化子进行PCR验证,将扩增出的大小正确的目的条带用OMEGA公司的胶回收试剂盒进行回收,并将回收到的PCR产物送到北京擎科生物科技有限公司进行测序,将测序结果与目的序列进行比对以确证克隆是否正确。选择验证正确的含有质粒pYES2-hph-egfp的阳性转化子进行激光共聚焦显微镜观察。
1.8 激光共聚焦显微镜分析
在载玻片上滴加40 μL无菌水,挑取验证正确的阳性转化子均匀分散到无菌水中,加上盖玻片后静置片刻,之后置于激光共聚焦显微镜下观察,选择480 nm激发波长用于增强型荧光蛋白EGFP的观察[21]。
2. 结果与分析
2.1 12C6+束辐照后的细胞存活曲线出现马鞍形波动
应用12C6+束对类酵母SD301菌株进行辐照,辐照剂量为20、60、90、120、150和180 Gy,细胞存活数分别为74、74、65、69、77和33。在辐照剂量为120、150和180 Gy时,细胞存活数出现了类似马鞍形的波动。12C6+束诱变的马鞍形效应在30多个属的微生物中被观测到,虽然具体的形成原因还不清楚,但马鞍形波动可能对应着较好的突变效果[15],因此选择马鞍形波动对应的经过120、150和180 Gy辐照的菌株进行突变效应筛选 (图1)。
2.2 H2O2可作为筛选压力筛选高产角鲨烯SD301突变株
研究发现,在本实验条件下所有H2O2处理组平板上均有明显的菌落生成,H2O2浓度达到100 μmol·L−1时,突变株菌落成活率保持在20%左右,所以选择100 μmol·L−1 H2O2作为筛选压力 (图2)。后续培养实验表明,该筛选压力下实验组表现出更高的角鲨烯质量 (22.53 mg·g−1 vs. 12.52 mg·g−1) 和产量 (267.31 mg·L−1 vs. 100.96 mg·L−1)。
2.3 高产角鲨烯SD301的突变株具有较好的遗传稳定性
经12C6+束对类酵母SD301菌株进行辐照,最终筛选到3株比出发菌株具有更高角鲨烯质量的突变株PS120、PS150和PS180。这3株菌株连续传代15次后,再次测定了角鲨烯质量,数据表明传代前后细胞内角鲨烯的质量无显著变化,获得的3株突变株经过传代后角鲨烯合成性能较为稳定 (图3)。
2.4 突变株比出发菌株具有更高的角鲨烯生产性能
摇瓶实验证实,突变株PS120的角鲨烯产量和质量均高于出发菌株SD301 (1.33 g·L−1 vs. 0.82 g·L−1,41.31 mg·g−1 vs. 26.82 mg·g−1,图4)。突变株PS120的角鲨烯产量明显高于同属菌株Pseudozyma sp. JCC 207 (0.34 g·L−1)[22],甚至高于经过代谢工程改造后的某些解脂耶氏酵母 (Yarrowia lipolytica ) (0.30 g·L−1)[23]、酿酒酵母 (0.21 g·L−1)[24]和大肠杆菌 (0.17 g·L−1)[25]。有研究发现,通过细胞质和过氧化物酶体双重调节,可进一步提高角鲨烯产量,最终通过优化发酵方法,可将角鲨烯产量提至11.00 g·L−1[17]。本文筛选到的突变株PS120作为遗传改造的出发菌株具有角鲨烯产量较高的起始优势,遗传改造手段的应用有望进一步大幅提高其角鲨烯产量,因此,本文将优先探索并建立稳定高效的电转化方法用于菌株的后续遗传改造。
2.5 电转化方法可以成功将质粒转入突变株PS120
潮霉素B通常对Pseudozyma属细胞具有抑制作用,因此本研究引入潮霉素B抗性基因作为构建重组质粒的抗性筛选标记[18-19]。本研究以半乳糖诱导型质粒pYES2为模板,将潮霉素B抗性基因HPH和编码增强型绿色荧光蛋白的EGFP基因引入载体,构建重组质粒pYES2-hph-egfp。通过菌落PCR扩增重组质粒中的EGFP基因 (EGFP基因片段的大小为1 014 bp),并对PCR产物进行琼脂糖凝胶验证、回收和测序研究,测序结果表明EGFP基因正确插入载体,重组质粒pYES2-hph-egfp构建成功 (图5-a)。
以突变株为出发菌株进行电转化,测试电压分别为650、900、1 150和1 400 V。通过菌落PCR扩增重组质粒中的HPH和EGFP基因表达盒 (该基因表达盒的大小为2 430 bp),并对PCR产物进行琼脂糖凝胶验证、回收和测序,测序结果表明重组质粒pYES2-hph-egfp成功转入到细胞中,成功获得含有质粒pYES2-hph-egfp的阳性转化子 (图5-b)。900 V电压条件下获得的阳性转化率 (30%,6个阳性克隆/20个克隆) 较1 400 V电压更高 (15%,3个阳性克隆/20个克隆)。此外,激光共聚焦显微镜观察也显示在900 V电压下获得的阳性转化子比1 400 V电压下有更多可以正常表达绿色荧光蛋白的细胞 (图6)。本研究使用的电压仅为900 V,远低于已报道的Pseudozyma属的电转化电压 (1 850 V[19]和5 kV[20]),操作上更具安全性。
3. 结论
本研究使用12C6+束辐照诱变产角鲨烯类酵母Pseudozyma sp. SD301,通过筛选初步获得了高产角鲨烯突变株PS120,摇瓶培养条件下其角鲨烯产量达到1.33 g·L−1,高于出发菌株SD301,是目前已报道过的使用诱变选育获得的最佳突变株之一[12,26]。PS120以葡萄糖为碳源进行1 L摇瓶培养3 d后菌体干质量 (干燥后细胞质量) 达22.10 g·L−1,高于产角鲨烯的南极酵母、酿酒酵母和裂殖壶菌[17,22,26-29],具备通过高密度发酵生产角鲨烯的潜力。除了通过诱变选育获得高产角鲨烯的突变株以外,通过代谢工程改造也可以有效提高角鲨烯的产量[17,30-34]。为了进一步提高突变株的角鲨烯产量,本研究改进了Pseudozyma属的电转化条件并优化了电压条件,结果显示在900 V电压下进行电转化可获得较高的阳性转化率 (30%),下一步工作重点是提高电转化效率和完善基因编辑系统。本研究结果可为产角鲨烯类酵母的遗传改造提供重要的技术支撑。
-
图 1 对照组及高温处理组中克氏原螯虾的肠道组织(200×)
注:a. C (26 ℃);b. M1 (32 ℃, 24 h);c. H1 (37 ℃, 24 h);d. M2 (32 ℃, 72 h);e. H2 (37 ℃, 72 h)。
Figure 1. Intestinal tissues of P. clarkii in control group and high-temperature treatment group (200×)
Note: a. C (26 ℃); b. M1 (32 ℃, 24 h); c. H1 (37 ℃, 24 h); d. M2 (32 ℃, 72 h); e. H2 (37 ℃, 72 h).
图 2 差异基因火山图
注:对照组 (26 ℃) 以C 表示,高温处理组 (32 ℃ 和37 ℃ 组) 分别以M 和H 表示;胁迫后的第24 小时以1 表示,第72 小时以2 表示。
Figure 2. Volcano map of differentially expressed genes
Note: C represents the control group (26 ℃); M and H represent the high-temperature treatment groups (32 ℃ and 37 ℃ groups), respectively; "1" represents the 24th hour after the stress, and "2" represents the 72nd hour after the stress.
图 3 差异基因聚类图
注:对照组 (26 ℃) 以C表示,高温处理组 (32 ℃ 和37 ℃ 组) 分别以M 和H 表示;胁迫后的第24 小时以1 表示,第72 小时以2 表示。横坐标表示样品名称及层次聚类结果,纵坐标表示差异基因及层次聚类结果,红色表示高表达,蓝色表示低表达。
Figure 3. Cluster map of differentially expressed genes
Note: In the samples, C represents the control group (26 ℃); M and H represent the high-temperature treatment groups (32 ℃ and 37 ℃ groups), respectively; "1" represents the 24th hour after the stress, and "2" represents the 72nd hour after the stress. The horizontal axis represents the sample name and hierarchical clustering results; the vertical axis represents the differentially expressed genes and hierarchical clustering results; red represents high expression, and blue represents low expression.
图 4 差异基因 GO 富集分析 (前15)
注:对照组 (26 ℃) 以C 表示,高温处理组 (32 ℃ 和37 ℃ 组) 分别以M 和H 表示;胁迫后的第24 小时以1 表示,第72 小时以2 表示。横坐标表示注释到该条目的差异基因数与差异基因总数的比值,纵坐标表示GO条目的名称。
Figure 4. Go enrichment analysis (Top 15)
Note: C represents the control group (26 ℃); M and H represent the high-temperature treatment groups (32 ℃ and 37 ℃ groups), respectively; "1" represents the 24th hour after the stress, and "2" represents the 72nd hour after the stress. The horizontal axis represents the ratio of the number of differential genes annotated to an entry to the total number of differential genes, and the vertical axis represents the name of the GO entry.
图 5 差异基因 KEGG 富集分析 (前20)
注:对照组 (26 ℃) 以C 表示,高温处理组 (32 ℃ 和37 ℃ 组) 分别以M 和H 表示;胁迫后的第24 小时以1 表示,第72 小时以2 表示。横坐标Rich factor 越大,富集的程度越大;纵坐标表示 KEGG 通路;点越大,通路富集的差异基因的数量越多;颜色越红,代表富集越显著。
Figure 5. KEGG enrichment analysis of differentially expressed genes (Top 20)
Note: C represents the control group (26 ℃); M and H represent the high-temperature treatment groups (32 ℃ and 37 ℃ groups), respectively; "1" represents the 24th hour after the stress, and "2" represents the 72nd hour after the stress. The larger the rich factor on the x-axis, the greater the degree of enrichment; the vertical axis represents the KEGG pathway; the larger the dot, the more differentially expressed genes are enriched in the pathway; the redder the color, the more significant the enrichment.
表 1 实验所用的引物
Table 1 Primer sequences used in this experiment
基因
Gene引物序列 (5'—3')
Primer sequence (5'–3')基因ID
Gene IDHSP71D F: CCTCCTTGATGTGGCACCTCTG LOC123759425 R: AGTGACCGCTGGTTGGTTATCTG HSP70 F: TGGCGGCTCCAGTAATTGGTATC LOC123774884 R: TGAGTCGTTCTGTGTCGGTGAAG HSP71 F: CCTACGGTGCGGCTGTTCAG LOC123759427 R: AAGACAGAGGTGCCACATCAAGG HSP90 F: GTGACCTCTCCGTGCTGTATCG LOC123770521 R: CTCGCAGTGCCTGAGCCTTC HSP7D F: CAGTCTGAGGGCGTGAAGGATG LOC123762431 R: TGAATACTTGCTGTTGCTTGGTTGG PD2 F: TGATCTGGTGAACAGGGTGCTAAC LOC123765724 R: ATGGATTCAACCGTGGCATTTGTG PN1AS F: TAGTGTGGGAGAAGGAGGAGAAGG LOC123756333 R: TGGCGGGCGAGATAGGAGTC PLP F: TGCGGTCATCGTTGGTGGTG LOC123755895 R: CATTAAGAGGCGGAGCGTTGTTC β-actin F: CACTGCCGCCTCATCCTCTTC LOC123745344 R: GAACCTCTCGTTGCCAATGGTAATG 表 2 转录组测序数据
Table 2 Data statistics of transcriptome
肠道样品
Intestinal sample测序数量
Raw reads/Gb过滤数量
Clean reads/Gb测序错误率
Error rateQ20占比
Q20 ratio/%Q30占比
Q30 ratio/%GC含量百分比
GC content/%C-1 53 149 186 51 496 446 0.03 97.33 92.50 44.37 C-2 51 693 424 49 938 434 0.03 97.20 92.21 42.72 C-3 51 340 602 49 503 334 0.03 97.18 92.22 43.79 M1-1 34 865 712 33 806 092 0.03 97.21 92.00 35.17 M1-2 43 523 082 42 563 458 0.03 97.25 92.06 35.93 M1-3 48 777 006 47 336 584 0.03 97.39 92.56 42.35 H1-1 49 274 690 48 127 640 0.03 97.43 92.57 40.25 H1-2 50 644 580 48 940 152 0.03 97.53 92.93 44.35 H1-3 51 066 558 49 289 724 0.03 97.56 93.01 44.62 M2-1 50 832 894 48 694 324 0.03 97.23 92.23 38.48 M2-2 49 202 048 46 946 612 0.03 97.26 92.31 35.95 M2-3 44 643 114 42 951 032 0.03 97.22 92.04 31.32 H2-1 48 783 458 46 894 270 0.03 97.35 92.52 39.70 H2-2 45 102 546 42 541 374 0.03 97.26 92.36 36.61 H2-3 51 312 374 48 790 394 0.03 97.30 92.45 36.66 注:对照组 (26 ℃)以C表示,高温处理组 (32 ℃和37 ℃组)分别以M和H表示;胁迫后的第24小时以1表示,第72 小时以2表示。 Note: C represents the control group (26 ℃); M and H represent the high-temperature treatment groups (32 ℃ and 37 ℃ groups), respectively; "1" represents the 24th hour after the stress, and "2" represents the 72nd hour after the stress. 表 3 Mapped Reads和Clean Reads比对情况统计
Table 3 Comparison statistics of Mapped Reads and Clean Reads
肠道样品
Intestinal sample总比对reads数
Total reads单一位点比
对率 (占比)
Unique mapped (Percentage)多位点比对率 (占比)
Multi mapped (Percentage)肠道样品
Intestinal sample总比对reads数
Total reads单一位点比
对率 (占比)
Unique mapped (Percentage)多位点比对率 (占比)
Multi mapped (Percentage)C-1 51 496 446 44 995 200 (87.38%) 2 204 550 (4.28%) H1-3 48 790 394 43 699 471 (89.57%) 2 232 553 (4.58%) C-2 49 938 434 44 427 790 (88.97%) 1 922 134 (3.85%) M2-1 33 806 092 30 711 261 (90.85%) 1 213 220 (3.59%) C-3 49 503 334 43 635 708 (88.15%) 1 943 516 (3.93%) M2-2 42 563 458 37 676 548 (88.52%) 3 365 504 (7.91%) M1-1 48 127 640 44 205 767 (91.85%) 1 418 030 (2.95%) M2-3 47 336 584 42 467 851 (89.71%) 1 908 883 (4.03%) M1-2 48 94 0152 43 851 068 (89.60%) 2 015 786 (4.12%) H2-1 48 694 324 43 573 397 (89.48%) 1 735 983 (3.57%) M1-3 49 289 724 43 822 561 (88.91%) 2 171 548 (4.41%) H2-2 46 946 612 41 788 201 (89.01%) 2 421 600 (5.16%) H1-1 46 894 270 42 198 853 (89.99%) 1 827 334 (3.90%) H2-3 42 951 032 38 155 003 (88.83%) 2 179 920 (5.08%) H1-2 42 541 374 37 729 522 (88.69%) 2 076 681 (4.88%) 注:对照组 (26 ℃)以C表示,高温处理组 (32 ℃和37 ℃组)分别以M和H表示;胁迫后的第24小时以1表示,第72 小时以2表示。 Note: C represents the control group (26 ℃); M and H represent the high-temperature treatment groups (32 ℃ and 37 ℃ groups), respectively; "1" represents the 24th hour after the stress, and "2" represents the 72nd hour after the stress. -
[1] YUE G H, LI J, BAI Z Y, et al. Genetic diversity and population structure of the invasive alien red swamp crayfish[J]. Biol Invasions, 2010, 12(8): 2697-2706. doi: 10.1007/s10530-009-9675-1
[2] XU Z Q, GAO T H, XU Y, et al. A chromosome-level reference genome of red swamp crayfish Procambarus clarkii provides insights into the gene families regarding growth or development in crustaceans[J]. Genomics, 2021, 113(5): 3274-3284. doi: 10.1016/j.ygeno.2021.07.017
[3] WAN J J, XI Q, TANG J Q, et al. Effects of pelleted and extruded feed on growth performance, intestinal histology and microbiota of juvenile red swamp crayfish (Procambarus clarkii)[J]. Animals, 2022, 12(17): 2252. doi: 10.3390/ani12172252
[4] 全国水产技术推广总站, 中国水产学会. 中国小龙虾产业发展报告(2024)[J]. 中国水产, 2024(7): 14-20. [5] 夏智宏, 刘敏, 秦鹏程, 等. 2022年长江流域高温干旱过程及其影响评估[J]. 人民长江, 2023, 54(2): 21-28. [6] LIU S Y, WANG Z Y, WANG Z, et al. Comparison of the gut microbiota and metabolism in different regions of red swamp crayfish (Procambarus clarkii)[J]. Front Microbiol, 2023, 14: 1289634. doi: 10.3389/fmicb.2023.1289634
[7] LI R H, BAI S H, YANG D C, et al. A crayfish Ras gene is involved in the defense against bacterial infection under high temperature[J]. Fish Shellfish Immunol, 2019, 86: 608-617. doi: 10.1016/j.fsi.2018.11.062
[8] 朱宏泽. 高温期间对虾病害防治[J]. 水产文摘, 2004(8): 6. [9] 凌宗欣, 唐立. 树突状细胞与肠道免疫[J]. 中国微生态学杂志, 2008, 20(3): 310-315. [10] RUAN X H, CHENG H T, SHAN J H, et al. Effect of eyestalk ablation and water temperature on the gonadal transcriptome of greasyback shrimp (Metapenaeus ensis)[J]. Aquac Rep, 2024, 39: 102374. doi: 10.1016/j.aqrep.2024.102374
[11] SHI W J, HU R H, WANG P, et al. Transcriptome analysis of acute high temperature-responsive genes and pathways in Palaemon gravieri[J]. Comp Biochem Phys D, 2022, 41: 100958.
[12] WANG Z Q, YANG L L, ZHOU F, et al. Integrated comparative transcriptome and weighted gene co-expression network analysis provide valuable insights into the response mechanisms of crayfish (Procambarus clarkii) to copper stress[J]. J Hazard Mater, 2023, 448: 130820. doi: 10.1016/j.jhazmat.2023.130820
[13] LIN W, WU J Y, LUO H M, et al. Sub-chronic ammonia exposure induces hepatopancreatic damage, oxidative stress, and immune dysfunction in red swamp crayfish (Procambarus clarkii)[J]. Ecotox Environ Safe, 2023, 254: 114724. doi: 10.1016/j.ecoenv.2023.114724
[14] LIU X, BIAN D D, JIANG J J, et al. Transcriptomic analyses provide new insights into immune response of the Procambarus clarkii intestines to Aeromonas hydrophila challenge[J]. Aquac Rep, 2024, 36: 102103. doi: 10.1016/j.aqrep.2024.102103
[15] LU X, TANG S, LIU X, et al. Comparative transcriptome analysis reveals the immune defense mechanism of the red swamp crayfish (Procambarus clarkii) in response to Vibrio parahaemolyticus challenge[J]. Aquaculture, 2024, 591: 741086. doi: 10.1016/j.aquaculture.2024.741086
[16] ZHANG L, SONG Z W, ZHONG S, et al. Acute hypoxia and reoxygenation induces oxidative stress, glycometabolism, and oxygen transport change in red swamp crayfish (Procambarus clarkii): application of transcriptome profiling in assessment of hypoxia[J]. Aquac Rep, 2022, 23: 101029. doi: 10.1016/j.aqrep.2022.101029
[17] 张龙岗, 董学飒, 朱永安, 等. 急性盐度胁迫对克氏原螯虾肝脏抗氧化酶及鳃丝Na+/K+-ATP酶活力的影响[J]. 淡水渔业, 2015, 45(4): 87-91. [18] XIAO Y, ZHANG Y M, XU W B, et al. The effects of salinities stress on histopathological changes, serum biochemical index, non-specific immune and transcriptome analysis in red swamp crayfish Procambarus clarkii[J]. Sci Total Environ, 2022, 840: 156502. doi: 10.1016/j.scitotenv.2022.156502
[19] LUO L, HUANG J H, LIU D L, et al. Transcriptome reveals the important role of metabolic imbalances, immune disorders and apoptosis in the treatment of Procambarus clarkii at super high temperature[J]. Comp Biochem Phys D, 2021, 37: 100781.
[20] 李铭, 董卫军, 邢迎春, 等. 温度对克氏原螯虾幼虾发育和存活的影响[J]. 水利渔业, 2006, 26(2): 36-37. [21] LI S X, LIN Y B, HE N J, et al. Antioxidation, immunity and hepatopancreatic histology of red swamp crayfish (Procambarus clarkii) subjected to low-temperature aerial exposure stress and re-immersion[J]. Comp Biochem Physiol A, 2023, 282: 111441. doi: 10.1016/j.cbpa.2023.111441
[22] WANG L, ZHU Q Q, HU M, et al. Toxic mechanisms of nanoplastics exposure at environmental concentrations on juvenile red swamp crayfish (Procambarus clarkii): from multiple perspectives[J]. Environ Pollut, 2024, 352: 124125. doi: 10.1016/j.envpol.2024.124125
[23] 李生瑄. 克氏原螯虾 (Procambarus clarkii) 对三种环境胁迫的生理生态响应[D]. 宜昌: 长江大学, 2023: 19-27. [24] SILVEYRA G R, SILVEYRA P, BRWN M, et al. Oxidative stress and histopathological effects by microplastic beads, in the crayfish Procambarus clarkii, and fiddler crab Leptuca pugilator[J]. Chemosphere, 2023, 343: 140260. doi: 10.1016/j.chemosphere.2023.140260
[25] WEI K Q, YANG J X. Oxidative damage induced by copper and beta-cypermethrin in gill of the freshwater crayfish Procambarus clarkii[J]. Ecotox Environ Safe, 2015, 113: 446-453. doi: 10.1016/j.ecoenv.2014.12.032
[26] CROLL S L, WATTS S A. The effect of temperature on feed consumption and nutrient absorption in Procambarus clarkii and Procambarus zonangulus[J]. J World Aquac Soc, 2004, 35(4): 478-488. doi: 10.1111/j.1749-7345.2004.tb00113.x
[27] FLOREA L, SONG L, SALZBERG S L. Thousands of exon skipping events differentiate among splicing patterns in sixteen human tissues[J]. F1000research, 2013, 2: 188. doi: 10.12688/f1000research.2-188.v1
[28] LIVAK K J, SCHMITTGEN T D. Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method[J]. Methods, 2001, 25(4): 402-408. doi: 10.1006/meth.2001.1262
[29] ZHANG X, ZHANG Y H, WU T, et al. Fish decay-accelerating factor (DAF) regulates intestinal complement pathway and immune response to bacterial challenge[J]. Fish Shellfish Immunol, 2024, 151: 109741.
[30] ZHOU Q C, BUENTELLO J A, GATLIN D M. Effects of dietary prebiotics on growth performance, immune response and intestinal morphology of red drum (Sciaenops ocellatus)[J]. Aquaculture, 2010, 309(1/2/3/4): 253-257.
[31] HUANG Y, REN Q. Research progress in innate immunity of freshwater crustaceans[J]. Dev Comp Immunol, 2020, 104: 103569. doi: 10.1016/j.dci.2019.103569
[32] 韩丰韦, 王哲奇, 刘晓静, 等. 岩藻多糖调控肠道黏膜屏障功能及其在动物生产中的应用[J]. 饲料研究, 2024(7): 151-156. [33] ZHOU J, WANG L, XIN Y, et al. Effect of temperature on antioxidant enzyme gene expression and stress protein response in white shrimp, Litopenaeus vannamei[J]. J Therm Biol, 2010, 35(6): 284-289. doi: 10.1016/j.jtherbio.2010.06.004
[34] DUAN Y F, WANG Y, ZHANG J S, et al. Elevated temperature disrupts the mucosal structure and induces an immune response in the intestine of whiteleg shrimp Litopenaeus vannamei (Boone, 1931) (Decapoda: Dendrobranchiata: Penaeidae)[J]. J Crustacean Biol, 2018, 38(5): 635-640. doi: 10.1093/jcbiol/ruy055
[35] 丁明孝, 王喜忠, 张传茂, 等. 细胞生物学[M]. 北京: 高等教育出版社, 2020: 83-84. [36] HENDERSHOT L M. The ER function BiP is a master regulator of ER function[J]. Mt Sinaij Med, 2004, 71(5): 289-297.
[37] WU Y J, PEI Y L, QIN Y H. Developmental expression of heat shock proteins 60, 70, 90, and A2 in rabbit testis[J]. Cell Tissue Res, 2011, 344(2): 355-363. doi: 10.1007/s00441-011-1151-4
[38] BALCH W E, MORIMOTO R I, DILLIN A, et al. Adapting proteostasis for disease intervention[J]. Science, 2008, 319(5865): 916-919. doi: 10.1126/science.1141448
[39] ZHANG X X, ZHU H C, YUAN J B, et al. Diversity of heat shock proteins in response to various stressors in the Pacific white shrimp Litopenaeus vannamei[J]. Aquaculture, 2024, 584: 740647. doi: 10.1016/j.aquaculture.2024.740647
[40] RUNGRASSAMEE W, LEELATANAWIT R, JIRAVANICHPAISAL P, et al. Expression and distribution of three heat shock protein genes under heat shock stress and under exposure to Vibrio harveyi in Penaeus monodon[J]. Dev Comp Immunol, 2010, 34(10): 1082-1089. doi: 10.1016/j.dci.2010.05.012
[41] HSIEH S L, LIAO W T, KUO C M. Fatty acid desaturation and stearoyl-CoA desaturase expression in the hepatic microsomes of milkfish, Chanos chanos, under cold shock[J]. 台湾水产学会刊, 2002, 29(1): 31-43. [42] GONZÁLEZ-RUÍZ R, LEYVA-CAMILLO L, PEREGÑNO-UÑARTE A B, et al. The combination of hypoxia and high temperature affects heat shock, anaerobic metabolism, and pentose phosphate pathway key components responses in the white shrimp (Litopenaeus vannamei)[J]. Cell Stress Chaperones, 2023, 28(5): 493-509. doi: 10.1007/s12192-022-01265-1
[43] OCAMPO L, PATIÑO D, RAMÍREZ C. Effect of temperature on hemolymph lactate and glucose concentrations in spiny lobster Panulirus interruptus during progressive hypoxia[J]. J Exp Mar Biol Ecol, 2003, 296(1): 71-77. doi: 10.1016/S0022-0981(03)00317-4
[44] QUEZADA-CALVILLO R, ROBAYO-TORRES C C, AO Z H, et al. Luminal substrate "brake" on mucosal maltase-glucoamylase activity regulates total rate of starch digestion to glucose[J]. J Pediatr Gastr Nutr, 2007, 45(1): 32-43. doi: 10.1097/MPG.0b013e31804216fc
[45] QUEZADA-CALVILLO R, ROBAYO-TORRES C C, OPEKUN A R, et al. Contribution of mucosal maltase-glucoamylase activities to mouse small intestinal starch α-Glucogenesis3[J]. J Nutr, 2007, 137(7): 1725-1733. doi: 10.1093/jn/137.7.1725
[46] VASUDEVAN S, RAJENDRAN S. Thermal stress induced hyperglycemia in the blue swimmer crab, Portunus pelagicus[J]. J Therm Biol, 2021, 100: 103076. doi: 10.1016/j.jtherbio.2021.103076
[47] YOHNA M A, RAY G W, YANG Q, et al. Protective effects of butyric acid during heat stress on the survival, immune response, histopathology, and gene rxpression in the hepatopancreas of juvenile Pacific shrimp (L. vannamei)[J]. Fish Shellfish Immunol, 2024, 150: 109610. doi: 10.1016/j.fsi.2024.109610