Sequencing and bioinformatic analysis for transcriptome of Micropterus salmoides based on RNA-seq
-
摘要:
文章以大口黑鲈(Micropterus salmoides)组织作为研究对象,利用RNA-seq技术进行转录本测序和数据分析,经拼接组装,最终获得35 659条unigenes,序列平均长度738 bp,序列长度中位数 (N50)为1 052 bp。另外从长度分布与GC含量等方面对unigenes进行评估,数据显示测序质量好、可信度高。使用6大数据库 (KOG、Nr、Pfam、Swiss-Prot、GO和KEGG) 注释大口黑鲈转录组unigenes,分别对应有15 832、21 279、14 524、16 973、15 024和11 185条unigenes获得注释。其中,5 617条unigenes在以上所有数据库中同时注释成功,17 253条unigenes至少被一个数据库注释。KEGG分析结果显示,获得注释的11 185条unigenes被划分到267个代谢通路中,参与信号转导通路的unigenes数量最多,共有1 349条(12.06%)。另外还检测到4 030个微卫星 (SSR)位点。通过对大口黑鲈转录组测序,获得了大量的转录组信息,为大口黑鲈的功能基因克隆、基因组学、遗传多样性分析、分子标记开发及遗传改良等研究奠定了基础。
Abstract:The transcripts of Micropterus salmoides were obtained by RNA-seq technology to conduct a transcriptomic analysis. A total of 35 659 unigenes were generated by de novo assembly with an average length of 738 bp and N50 of 1 052 bp. The unigenes were assessed for length distribution and GC content. The sequencing data were of high quality and reliability. A total of 15 832, 21 279, 14 524, 16 973, 15 024 and 11 185 unigenes were annotated from the KOG, Nr, Pfam, Swiss-Prot, GO and KEGG databases, respectively, among which 5 617 were annotated in all the databases and 17 253 were annotated in at least one database. Blasted with KEGG pathway, 11 185 unigenes were annotated, belonging to 267 categories. Of all the pathways, the number of unigenes joining in the signal transduction was the most (1 349, 12.06%). Finally, 4 030 SSRs were identified. The results provide rich data to understand transcriptome information of M. salmoides and lay the foundation for further research on functional gene cloning, genomics, genetic diversity analysis, molecular marker exploitation and genetic improvement in M. salmoides.
-
Keywords:
- Micropterus salmoides /
- transcriptome /
- unigenes /
- gene annotation
-
三角帆蚌 (Hyriopsis cumingii) 所产珍珠占中国珍珠产量的80%以上[1],且品质良好、经济价值较高。但三角帆蚌培育珍珠的大小、光泽度等均与性别有关,在产珠性能方面雄性优于雌性[2],而人工雌雄分养情况下能促进雌性三角帆蚌的个体发育及产珠性能[3]。三角帆蚌属于雌雄异体,不存在异型染色体,目前性别决定无法从染色体层面进行研究,所以性别分化与决定相关基因的研究成为性别调控的热点,对揭示复杂的性别决定与分化机制研究及创造经济效益具有重大意义[4]。
丝裂原活化蛋白激酶 (MAPK) 也称为细胞外信号调节激酶,受各种细胞外生长因子及受体的相互作用而被激活,会影响许多组织特异性生物活性,如细胞增殖、存活和分化[5-6]。然而,MAPK在雌性生殖细胞中不遵循这种传统的模式[7]。在卵母细胞中,MAPK独立于生长因子和酪氨酸激酶受体而被激活,独立于转录调控发挥作用[8],受特殊上游调节剂激酶 MOS的控制,在卵母细胞减数分裂中起关键作用[9-10]。Mos/MAPK通路在脊椎动物中启动卵母细胞的成熟发育[10],C-MOS基因在Mos/MAPK通路中位于MAP2K1基因的上游。MAP2K1基因作为MAPK家族的一分子,参与细胞增殖、分化、基因转录调控和发育过程[11-13]。在脊椎动物中,MAP2K1基因参与小鼠 (Mus musculus) 的卵巢发育[14],促进水牛 (Bubalus) 卵母细胞的表达与发育[15],对鹅 (Anser cygnoides orientalis) 的卵泡发育起到支持与促进作用[5]。在无脊椎动物中,MAP2K1基因对斑节对虾 (Penaeus monodon) 卵母细胞的生长发育起促进作用[14];能促进果蝇 (Drosophilid) 的卵母细胞成熟[16]。
在贝类中,紫贻贝 (Mytilus galloprovincialis)、虾夷扇贝 (Mizuhopecten yessoensis) 的MAP2K1基因均有报道[17-18]。本文对三角帆蚌MAP2K1基因进行了克隆,分析了该基因在性腺不同阶段和雌雄性腺组织中的表达模式,探讨了该基因在三角帆蚌性腺发育中的作用。
1. 材料与方法
1.1 实验材料
从浙江省金华市武义实验基地选取健康的1—8月龄和1、2、3龄三角帆蚌,暂养1周后对1—8月龄个体的性腺组织、1龄和2龄个体各组织 (性腺、闭壳肌、肝胰腺、鳃、外套膜、斧足)、3龄个体的性腺组织进行取样,样本立刻放入液氮,速冻后用镊子放入样本盒,−80 ℃保存。
1.2 总RNA提取与cDNA合成
按TRIzol法提取总RNA,利用NanoDrop 2000分度计进行核酸检测,通过1%琼脂糖凝胶电泳[19] (电压180 V,电流200 mA) 检查RNA完整性情况 (3条条带5.8S、18S、28S)。检查合格后将RNA保存于−80 ℃。将RNA进行反转录 (Prime ScriptTM RT reagent Kit with gDNA Eraser试剂盒,TaKaRa, 日本) 得到cDNA模板,cDNA保存于−40 ℃。
1.3 基因序列的获取
通过转录组库[20]中得到的序列设计3'RACE的inner和outer引物 (表1),并参照SMART 3'RACE试剂盒 (Clontech, America) 说明书进行MAP2K1基因的3'端PCR扩增。电泳 (电压180 V,电流200 mA) 检测产物,割胶回收纯化。用PMD19-T (TaKaRa,大连) 载体连接纯化产物,16 ℃下连接4 h。准备冰盒,在冰上将连接产物转化入大肠杆菌DH5α (TaKaRa,大连) 中 30 min。准备含氨苄的固体培养基 (4 ℃冰箱保存),将所得菌液均匀地涂在上面,37 ℃培养10 h。挑选白色菌株送至生工 (上海) 测序。
表 1 本实验所需引物Table 1. Primers for this study引物名
Primer's name引物序列 (5'—3')
Primer sequence (5'−3')用途
PurposeK1-O TATTCCTACAGCAGACTCGTGATTG 3'RACE outer K1-I GACCTGGATTTGATAATAGATGTGGG 3'RACE inner K1-YG-F AACAAGCTGAATCTGACGCTG 荧光定量正向引物 K1-YG-R GCTCCTTCAGTTTCTTGGTCAGT 荧光定量反向引物 EFl-αF GGAACTTCCCAGGCAGACTGTGC 内参正向引物 EFl-αR TCAAAACGGGCCGCAGAGAAT 内参反向引物 K1-YW-F AACAAGCTGAATCTGACGCTG 原位杂交正向引物 K1-YW-R TAATACGACTCACTATAGGGGCTCCTTCAGTTTCTTGGTCAGT 原位杂交反向引物 M-GR-1F TTGGGTATGGTTCGCAGTT 干扰链1正向引物 M-GR-1R GCCCTTGTAAGCTCATGTTT 干扰链1正向引物 T7+M-GR-1F TAATACGACTCACTATAGGGTTGGGTATGGTTCGCAGTT T7+干扰链1正向引物 T7+M-GR-1R TAATACGACTCACTATAGGGGCCCTTGTAAGCTCATGTTT T7+干扰链1反向引物 M-GR-2F TATGAACGAACATCGGAGGC 干扰链2正向引物 M-GR-2R CCCATACGCAACAACCC 干扰链2反向引物 T7+M-GR-2F TAATACGACTCACTATAGGGTATGAACGAACATCGGAGGC T7+干扰链2正向引物 T7+M-GR-2R TAATACGACTCACTATAGGGCCCATACGCAACAACCC T7+干扰链2反向引物 M-GR-3F TAACGCCTCACGACGAC 干扰链3正向引物 M-GR-3R CCGGCTCTTGCATTCC 干扰链3反向引物 T7+M-GR-3F TAATACGACTCACTATAGGGTAACGCCTCACGACGAC T7+干扰链3正向引物 T7+M-GR-3R TAATACGACTCACTATAGGGCCGGCTCTTGCATTCC T7+干扰链3反向引物 1.4 基因生物信息学分析
使用NCBI (https://www.ncbi.nlm.nih.gov/) 中的ORF finder预测MAP2K1基因的开放阅读框 (ORF)、Blast得到不同物种MAP2K1基因序列、Premier设计引物;蛋白基本理化性质利用ProtParam tool (https://web.expasy.org/aprotparam/) 分析;二级结构利用SOPMA (https://npsa-prabi.ibcp.fr/NPSA/npsa_sopma.html) 预测;三级结构利用I-TASSER (https://zhanglab.ccmb.med.umich.edu/I-TASSER/) 预测;跨膜结构域利用TMHMM Serverv 2.0程序 (http://www.cbs.dtu.dk/services/TMHMM/) 预测并用 (http://www.genedenovo.com/news/582.html) 解读结果;信号肽利用SignalP-5.0 (http://www.cbs.dtu.dk/services/SignalP/) 预测并用 (http://www.cbs.dtu.dk/services/SignalP-3.0/output.php) 解读结果,基因表达量差别利用SPSS 18.0软件分析;表达量差别条形图利用SigmaPlot 12.5绘制;NJ (邻接法) 系统进化树利用MEGA 7.0构建,Bootstrap分析其可靠性,重复1 000次。
1.5 荧光定量分析
选取EF-1α为内参基因,使用CFX96 TouchTM Real-Time PCR Detection System (Bio-Rad) 仪器进行定量分析。定量体系 (20 μL) 为cDNA (使用前混合均匀) 1.6 μL,2×TB Green Premix ExTaq 10 μL,F (上游引物) 0.8 μL,R (下游引物) 0.8 μL,RNase Free Water 6.8 μL。定量程序为95 ℃ 3 min;95 ℃ 5 s,60 ℃ 30 s,循环40次。
1.6 原位杂交
原位杂交的引物具有特异性,上游引物保持不变,将T7启动子序列 (TAATACGACTCACTATAGGG) 加入到下游引物5'端前。三角帆蚌性腺cDNA作为模板进行扩增,得到目的条带割胶回收,纯化后作为模板,使用T7 (T7 High Efficiency Transcription) 试剂盒和DIG RNA Labeling Mix得到标记探针,cDNA 1 µL,4 μL 5×T7 Transcription Reaction、8 μL 10 mmol·L−1 NTP Mix,2 μL T7 Transcription Enzyme Mix,RNase-free water 5 μL。反应程序为37 ℃,2 h。将所得探针纯化后保存于−80 ℃。
取2龄健康雌、雄三角帆蚌的性腺组织 (组织上连带少许表皮),将其在4%多聚甲醛与DEPC中固定2 h,转移至70%乙醇中,在乙醇中可保存约10 d。进而进行组织梯度脱水、组织透明、浸蜡、包埋、切片 (厚度控制在6~8 μm)、粘片。利用DIG nucleic acid detection kit进行原位杂交,并进行封片处理,在显微镜下观察杂交信号并拍照。
1.7 RNA干扰
1.7.1 dsRNA的合成
使用Primer Premier 5.0软件在C-MOS基因的ORF区设计3对引物,在上下游引物前加上T7序列,将每组中正常引物与加有T7序列的引物交叉重组,共形成6对引物,分别进行PCR扩增,扩增后回收纯化cDNA序列 (质量浓度≥167 ng·μL−1)。利用T7 High Efficiency Transcription 试剂盒 (全式金,北京) 将cDNA逆转录成单链RNA,混匀,37 ℃ PCR孵育2 h。将同一片段的两管单链RNA在冰上等量混匀,将混匀物放置在70 ℃ PCR仪中反应10 min,反应结束立刻取出后室温静置20 min。向每20 μL的反应液加入1 μL RNase A (用无酶水稀释200倍) 和1 μL DNase I,放置在37 ℃ PCR仪中反应30 min。这一步骤是为了去除DNA,使之退火合成dsRNA。将所得溶液中加入1/10体积的3 mol·L−1醋酸铵和等量异丙醇,冰中静置5 min;放置离心机中12 000 r·min−1、4 ℃、10 min。弃上清,70%乙醇洗涤沉淀;12 000 r·min−1、4 ℃离心10 min。超净工作台内风干15 min;加入100 μL RNase-free water溶解。检测dsRNA的质量和浓度,−80 ℃保存。
1.7.2 dsRNA干扰链注射
取1龄三角帆蚌按照每组15只分为4组 (干扰链1、干扰链2、干扰链3和对照组),实验组每只注射50 μL (200 ng·μL−1) 干扰链,阴性对照组每只注射50 μL (200 ng·μL−1) 生理盐水 (Normal saline, NS)。在第12、第24和第48小时提取性腺组织并立刻放置于液氮中。利用TRIzol法提取总RNA,试剂盒将其反转为cDNA,将反转后的cDNA作为模板进行荧光定量实验,检测直接干扰C-MOS基因后在性腺中的表达量情况及下游基因MAP2K1的表达量情况。
2. 结果
2.1 MAP2K1基因cDNA的克隆及序列特征
本实验克隆得到三角帆蚌MAP2K1基因cDNA的序列,其中5' 非编码区 (UTR) 为145 bp,3'UTR为2 070 bp,开放阅读框 (ORF) 1 194 bp,编码397个氨基酸 (图1)。预测相对分子质量44.24 kD,等电点为6.81,平均亲水系数为−0.379,推测为亲水蛋白。根据MAP2K1蛋白既不具备跨膜结构、也不具备信号肽的特点,推测其不属于膜蛋白,而属于胞内蛋白。通过SOPMA分析得到MAP2K1基因具有S-TKC结构域 (72—372 aa,图2-a),通过I-TASSER预测MAP2K1基因的蛋白质三级结构 (图2-b)。
2.2 MAP2K1基因系统进化
利用GeneDoc进行不同物种之间的序列比对,结果显示来自多种物种的MAP2K1基因结构域重合度很高 (图3),说明在多种物种之间,MAP2K1基因进化和结构具有保守性,推断其具有功能相似性。
根据同源性比对结果分析,三角帆蚌MAP2K1基因与紫贻贝、虾夷扇贝、泥蚶 (Tegillarca granosa) 的同源性较高,同源性分别为85.52%、82.65%、80.59%。与智人 (Homo sapiens)、北美属兔 (Ochotona princeps) 等同源性较低,同源性分别为72.27%、72.53%。系统进化树显示三角帆蚌与紫贻贝、虾夷扇贝、泥蚶聚为一支,亲缘关系较近 (图4)。
2.3 MAP2K1基因在各组织和各时期的表达量差异
在雌雄三角帆蚌的不同组织中,MAP2K1基因均有一定的表达量,但是MAP2K1基因在雌性性腺中的相对表达量高于另5个组织,雌雄差异也比较显著 (P<0.05,图5)。
在三角帆蚌幼龄 (1—8月) 性腺发育阶段过程中,MAP2K1基因均有一定程度的表达 (图6)。说明MAP2K1基因在三角帆蚌幼龄性腺发育阶段起到一定作用。
三角帆蚌1龄性腺中,雌性相对表达量显著高于雄性 (P<0.01),但在雄性中也有一定表达 (图7);等到2龄左右,性腺相对表达量雌雄差异极显著 (P<0.01),雄性中几乎不表达;在3龄性腺中,雌雄表达量均处于非常低的水平。结果表明,MAP2K1基因在三角帆蚌2龄卵巢发育过程中起到非常重要的作用。
2.4 性腺组织原位杂交
结果显示,实验组中三角帆蚌雌性性腺的卵母细胞、卵子上,蓝紫色杂交信号非常显著。实验组中三角帆蚌雄性性腺的精原细胞、精母细胞上,蓝紫色杂交信号存在极少,雌性及雄性阴性对照组中均未出现杂交信号 (图8)。根据结果推测,MAP2K1基因可能参与三角帆蚌的卵巢发育过程。
2.5 RNA干扰结果
三条干扰链的干扰结果为,干扰链1在雌性和雄性中的干扰率分别为88.43%和77.90%;干扰链2在雌性和雄性中的干扰率分别为93.10%和78.10%;干扰链3在雌性和雄性中的干扰率分别为96.39%和86.26% (图9)。干扰率结果显示干扰链3对C-MOS基因的干扰效果最佳,为此选择干扰链3进行实验,检测C-MOS基因在12、24和48 h干扰后在性腺中的表达情况,结果显示干扰链3在第24小时的干扰效率最佳 (图10)。
干扰C-MOS基因后,荧光定量检测其下游基因MAP2K1的表达量情况 (图11)。结果显示,MAP2K1基因在雌性中的表达量下降了82.31%,在雄性中的表达量下降了73.60%。说明在三角帆蚌中,C-MOS基因对MAP2K1基因有显著的调控作用。
3. 讨论
本研究克隆了三角帆蚌MAP2K1基因的cDNA序列,编码一个长为1 194 bp的ORF。蛋白质结构预测分析显示,MAP2K1基因中含有S-TKC催化结构域。蛋白苏氨酸/丝氨酸激酶是十分常见的一类蛋白激酶[21],主要作用是蛋白磷酸化、去磷酸化,有利于细胞内活性[22]。S-TKC是这类激酶中一个位于核心部位重要的催化结构域[23]。在S-TKC中有一个极其保守的区域,在催化结构域它已被证明是涉及ATP的结合,即赖氨酸残基附近N-末端的一个富含甘氨酸的残基。在S-TKC的中间区域,存在一个天冬氨酸残基,在催化结构域它已被证明对酶的催化活性起显著作用[24-26]。拥有S-TKC催化结构域的基因通常参与许多生理生化过程[27-30],若其构造折叠或变化,会对丝氨酸/苏氨酸蛋白激酶产生影响,从而对个体的生殖发育过程造成极大损害[31],这说明MAP2K1基因可能参与了三角帆蚌个体的生殖发育过程。
不同物种之间的序列比对结果显示,MAP2K1基因在三角帆蚌中的氨基酸序列与其他物种结构域重合度较高,说明与其他物种相似度较高,尤其与虾夷扇贝等软体动物相似度更高,由此推断在进化过程中,MAP2K1基因呈现高度保守状态。各组织、早期发育阶段、不同发育时期的定量结果说明,MAP2K1基因在卵巢中特异性高表达,在1—8月龄时期均有不同程度的表达,在1、2、3龄卵巢中的表达量均高于睾丸,对小鼠[14]、水牛[15]、鹅[5]、斑节对虾[14]及果蝇[16]MAP2K1基因的研究中,均发现雌性表达量高于雄性,与本实验结果相近,由此推测MAP2K1基因是三角帆蚌中的偏雌性基因,既参与了三角帆蚌的早期发育阶段,也参与了卵巢发育过程。根据原位杂交的探针定位显示,MAP2K1基因在卵母细胞及卵子上存在明显的信号,推测MAP2K1基因可能参与了三角帆蚌的卵巢发育过程。
RNA干扰结果显示,3条干扰链中干扰链2和干扰链3对C-MOS基因起到抑制作用,但是干扰链3的干扰效果最佳,在雌性和雄性中的干扰率分别为84.21%和96.87%。在选取干扰链3观察干扰C-MOS基因后MAP2K1基因的表达情况中,MAP2K1基因在雌性中的表达量下降了82.31%,在雄性中的表达量下降了73.60%。结果显示在三角帆蚌中,C-MOS基因对MAP2K1基因有显著的调控作用。在脊椎动物卵母细胞MAPK信号通路中,MAPK独立于生长因子和酪氨酸激酶受体而被激活,受特殊上游调节剂激酶MOS的控制,在卵母细胞减数分裂中起关键作用[9]。综上,推测在三角帆蚌中,MAP2K1基因可能参与了三角帆蚌的卵巢发育过程,是偏雌性基因。
-
图 3 Unigenes GO功能注释
1. 转录、DNA依赖性;2. 转录调控、DNA依赖性;3. 蛋白质转运;4. 凋亡;5. 细胞分化;6. 细胞黏附;7. 蛋白质水解;8. 多细胞器官发育;9. 细胞周期;10. 细胞分化;11. 胞内信号转导;12. mRNA加工;13. 信号转导;14. 小G蛋白介导信号转导;15. 转运;16. 转录正调控;17. 有丝分裂;18. 染色质修饰;19. 转录负调控;20. RNA剪切;21. 胞内蛋白质运输;22. 翻译;23. DNA修复;24. 转录;25. 内噬作用;26. 细胞核;27. 细胞浆;28. 必须膜;29. 胞液;30. 细胞质膜;31. 内质网膜;32. 核仁;33. 线粒体;34. 胞外区;35. 核质;36. 隔膜;37. 高尔基体膜;38. 细胞骨架;39. 高尔基体;40. 细胞核周区;41. ATP结合;42. 锌离子结合;43. 蛋白结合;44. DNA结合;45. 金属离子结合;46. RNA结合;47. 钙离子结合;48. 特异序列DNA结合;49. 结合;50. 蛋白质丝氨酸特异酶结合
Figure 3. GO functional annotation of unigenes
1. transcription, DNA-dependent; 2. regulation of transcription, DNA-dependent; 3. protein transport; 4. apoptosis; 5. cell division; 6. cell adhesion; 7. proteolysis; 8. multicellular organismal development; 9. cell cycle; 10. cell differentiation; 11. intracellular signal transduction; 12. mRNA procession; 13. signal transduction; 14. small GTPase mediated signal transduction; 15. transport; 16. negative regulation of transcription; 17. mitosis; 18. chromatin modification; 19. positive regulation of transcription; 20. RNA splicing; 21. intracellular protein; 22. translation; 23. DNA repair; 24. transcription; 25. endocytosis; 26. nucleus; 27. cytoplasm; 28. integral to membrane; 29. cytosol; 30. plasma membrane; 31. endoplasmic reticulum membrane; 32. nuclelous; 33. mitochondrion; 34. extracellular region; 35. nucleoplasm; 36. membrane; 37. golgi membrane; 38. cytoskeleton; 39. golgi apparatus; 40. perinuclear region of cytoplasm; 41. ATP binding; 42. zinc ion binding; 43. protein binding; 44. DNA binding; 45. metal ion binding; 46. RNA binding; 47. calcium ion binding; 48. specific DNA sequence binding; 49. binding; 50. protein serine/threonine-specific kinase
图 4 Unigenes的KOG注释
A. RNA加工和修饰;B. 染色体结构和动力学;C. 能源生产和转换;D. 细胞周期调控-细胞分裂-染色体分离;E. 氨基酸转运和代谢;F. 氨基酸转运和代谢;G. 碳水化合物转运和代谢;H. 辅酶转运和代谢;I. 脂质转运和代谢;J. 翻译-核糖体结构-生物合成;K. 转录;L. 复制-重组-修复;M. 细胞壁-细胞膜合成;N. 细胞运动;O. 翻译后修饰-蛋白质周转-分子伴侣;P. 无机离子转运与代谢;Q. 次生代谢产物生物合成、转运和分解代谢;R. 一般功能预测;S. 功能未知;T. 信号转导机制;U. 细胞内分泌和囊泡运输;V. 防御机制;W. 胞外结构;Y. 核结构;Z. 细胞骨架
Figure 4. KOG annotation of unigenes
A. RNA processing and modification; B. chromatin structure and dynamics; C. energy production and conversion; D. cell cycle control, cell division, chromosome partitioning; E. amino acid transport and metabolism; F. nucleotide transport and metabolism; G. carbohydrate transport and metabolism; H. coenzyme transport and metabolism; I. lipid transport and metabolism; J. translation, ribosomal structure and biogenesis; K. transcription; L. replication, recombination and repair; M. cell wall/membrane/envelope biogenesis; N. cell motility; O. posttranslational modification, protein turnover, chaperones; P. inorganic ion transport and metabolism; Q. secondary metabolites biosynthesis, transport and catabolism; R. general function prediction only; S. function unknown; T. signal transduction mechanisms; U. intracellular trafficking, secretion, and vesicular transport; V. defense mechanisms; W. extracellular structure; Y. nuclear structure; Z. cytoskeleton
表 1 转录组数据拼接结果统计
Table 1 Assembly result of transcriptome data
序列种类
category总数
total numberN50长度/bp
N50 length总长度/bp
total length最大长度/bp
maximum length中等长度/bp
median length最小长度/bp
minimum length平均长度/bp
average lengthtranscript 42 264 1 331 36 549 208 35 516 522 201 864 unigene 35 659 1 052 26 340 493 35 516 472 201 738 注:N50表示将转录本从长到短排序,依次累加碱基数,当累计碱基数达到转录本总碱基数的50%时的转录本的长度 Note: N50 of transcript or unigenes was calculated by ordering all sequences, then adding the lengths from longest to shortest until the summed length exceeded 50% of the total length of all sequences. -
[1] LI Y, ZHOU Z, TIAN M, et al. Exploring single nucleotide polymorphism (SNP), microsatellite (SSR) and differentially expressed genes in the jellyfish (Rhopilema esculentum) by transcriptome sequencing[J]. Mar Genom, 2017, 34: 31-37. doi: 10.1016/j.margen.2017.01.007
[2] CARRUTHERS M, YURCHENKO A A, AUGLEY J J, et al. De novo transcriptome assembly, annotation and comparison of four ecological and evolutionary model salmonid fish species[J]. BMC Genom, 2018, 19(1): 32. doi: 10.1186/s12864-017-4379-x
[3] JORGE P H, MASTROCHIRCO F A, HATA M E, et al. Genetic characterization of the fish Piaractus brachypomus by microsatellites derived from transcriptome sequencing[J]. Front Genet, 2018, 9: 46. doi: 10.3389/fgene.2018.00046
[4] 罗辉, 叶华, 肖世俊, 等. 转录组学技术在水产动物研究中的运用[J]. 水产学报, 2015, 39(4): 598-607. [5] CUI J, XU J, ZHANG S, et al. Transcriptional profiling reveals differential gene expression of Amur ide (Leuciscus waleckii) during spawning migration[J]. Int J Mol Sci, 2015, 16(6): 13959-13972.
[6] 邱忠营, 刘菲, 张克瑶, 等. 疣蝗转录组分析[J]. 基因组学与应用生物学, 2016, 35(8): 1989-1998. [7] SUDHAGAR A, KUMAR G, EL-MATBOULI M. Transcriptome analysis based on RNA-Seq in understanding pathogenic mechanisms of diseases and the immune system of fish: a comprehensive review[J]. Int J Mol Sci, 2018, 19(1): 1-19.
[8] CAO S, ZHU L, NIE H, et al. De novo assembly, gene annotation, and marker development using Illumina paired-end transcriptome sequencing in the Crassadoma gigantea[J]. Gene, 2018, 658: 54-62. doi: 10.1016/j.gene.2018.03.019
[9] JIA Z, WANG Q, WU K, et al. De novo transcriptome sequencing and comparative analysis to discover genes involved in ovarian maturity in Strongylocentrotus nudus[J]. Comp Biochem Physiol D, 2017, 23: 27-38.
[10] LU T, SUN Y, MA Q, et al. De novo transcriptomic analysis and development of EST-SSR markers in the Siberian tiger (Panthera tigris altaica)[J]. Mol Genet Genom, 2016, 291(6): 2145-2157. doi: 10.1007/s00438-016-1246-4
[11] RAI R, CHAUHAN S K, SINGH V V, et al. RNA-seq analysis reveals unique transcriptome signatures in systemic lupus erythematosus patients with distinct autoantibody specificities[J]. PLoS One, 2016, 11(11): e0166312. doi: 10.1371/journal.pone.0166312
[12] LI S, SHEN L, SUN L, et al. Small RNA-Seq analysis reveals microRNA-regulation of the Imd pathway during Escherichia coli infection in Drosophila[J]. Dev Comp Immunol, 2017, 70: 80-87. doi: 10.1016/j.dci.2017.01.008
[13] SALEM M, REXROAD C E, WANG J N, et al. Characterization of the rainbow trout transcriptome using Sanger and 454-pyrosequencing approaches[J]. BMC Genomics, 2010, 11: 564. doi: 10.1186/1471-2164-11-564
[14] COPPE A, PUJOLAR J M, MAES G E, et al. Sequencing, de novo annotation and analysis of the first Anguilla anguilla transcriptome: EeelBase opens new perspectives for the study of the critically endangered european eel[J]. BMC Genomics, 2010, 11: 635. doi: 10.1186/1471-2164-11-635
[15] PEREIRO P, BALSEIRO P, ROMERO A, et al. High-throughput sequence analysis of turbot (Scophthalmus maximus) transcriptome using 454-pyrosequencing for the discovery of antiviral immune genes[J]. PLoS One, 2012, 7(5): e35369. doi: 10.1371/journal.pone.0035369
[16] 赵文, 高峰英, 石振广. 达氏鳇鱼肌肉组织转录组测序与功能分析[J]. 水产学报, 2014, 38(9): 1255-1262. [17] MU Y N, DING F, CUI P, et al. Transcriptome and expression profiling analysis revealed changes of multiple signaling pathways involved in immunity in the large yellow croaker during Aeromonas hydrophila infection[J]. BMC Genom, 2010, 11: 506. doi: 10.1186/1471-2164-11-506
[18] WANG R J, SUN L Y, BAO L S, et al. Bulk segregant RNA-seq reveals expression and positional candidate genes and allele-specific expression for disease resistance against enteric septicemia of catfish[J]. BMC Genom, 2013, 14: 929. doi: 10.1186/1471-2164-14-929
[19] SUTHERLAND B J, KOCZKA K W, YASUIKE M, et al. Comparative transcriptomics of Atlantic Salmo salar, chum Oncorhynchus keta and pink salmon O. gorbuscha during infections with salmon lice Lepeophtheirus salmonis[J]. BMC Genom, 2014, 15(1): 200. doi: 10.1186/1471-2164-15-200
[20] DANG Y, XU X, SHEN Y, et al. Transcriptome analysis of the innate immunity-related complement system in spleen tissue of Ctenopharyngodon idella infected with Aeromonas hydrophila[J]. PLoS One, 2016, 11(7): e0157413. doi: 10.1371/journal.pone.0157413
[21] 赵刚, 龚全, 刘亚, 等. 基于Illumina高通量测序的岩原鲤转录组分析[J]. 西南农业学报, 2016, 29(7): 1743-1749. [22] 许建, 赵建, 徐礼鸣, 等. 基于RNA-Seq技术的鲮转录组分析[J]. 大连海洋大学学报, 2014, 29(6): 556-560. doi: 10.3969/J.ISSN.2095-1388.2014.06.003 [23] 全迎春, 马冬梅, 白俊杰, 等. 大口黑鲈转录组SNPs筛选及其与生长的关联分析[J]. 水生生物学报, 2016, 11(40): 294-300. -
期刊类型引用(1)
1. 路俊怡,姜立妍,龙凯,王韬,吴正理,李艳红. HcTLR1通过MyD88-NF-κB信号通路参与三角帆蚌抗菌免疫应答. 南方水产科学. 2024(06): 19-30 . 本站查看
其他类型引用(0)