鼩鼱科三物种转录组初步分析

doi:10.16829/j.slxb.150403

兽类学报 ›› 2021, Vol. 41 ›› Issue (2): 173-181.DOI: 10.16829/j.slxb.150403

鼩鼱科三物种转录组初步分析

周婧晶¹, 丘银彬¹, 万韬², 何水旺³, 蒋学龙³, 潘星华¹, 何锴^1,3

1 南方医科大学基础医学院, 广东省单细胞技术与应用重点实验室, 广州 510515;
2 陕西师范大学生命科学学院, 西安 710119;
3 中国科学院昆明动物研究所, 昆明 650223

收稿日期:2019-12-30 出版日期:2021-03-30 发布日期:2021-04-02
通讯作者: 潘星华,E-mail:panvictor@qq.com;何锴,E-mail:hekai2018@smu.edu.cn
作者简介:周婧晶(1995-),女,硕士研究生,主要从事进化生物学研究.
基金资助:
国家自然科学基金面上基金（31970389）；中国博士后科学基金面上项目（2009M652952）；广东省珠江人才计划本土创新科研团队（2017BT01S131）

Transcriptome analyses for three Soricid shrews

ZHOU Jingjing¹, QIU Yinbin¹, WAN Tao², HE Shuiwang³, JIANG Xuelong³, PAN Xinghua¹, HE Kai^1,3

1 Southern Medical University, School of Basic Medical Science, Guangdong Provincial Key Laboratory of Single Cell Technology and Application, Guangzhou 510515, China;
2 College of Life Sciences, Shaanxi Normal University, Xi'an 710119, China;
3 Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming 650223, China

Received:2019-12-30 Online:2021-03-30 Published:2021-04-02

摘要/Abstract

摘要： 中国是鼩鼱科物种十分丰富的国家，已知有12属61个物种在中国分布，并演化出多种生态适应型。本研究采用高通量测序技术对小纹背鼩鼱、云南长尾鼩鼱和蹼足鼩3个物种的心和肺组织样本进行转录组测序。其中小纹背鼩鼱和云南长尾鼩营地表生活，蹼足鼩适应半水生生活。3个物种的转录组测序一共获得20.5 Gb的Clean Data，拼接后分别获得131 045条、153 621条和135 838条Unigenes，BUSCO分析显示其中完整的保守基因最高为82.2%，共有同源基因家族超过8 000个。此外，我们发现通过优化分析方案，可以将单拷贝基因的比例提高18%~25%，同时将多拷贝基因比例从28.4%~31.8%降为0.7%~1.1%。同物种不同组织间的基因差异表达分析找到各组织的标志性基因。转录组差异表达分析发现，蹼足鼩心脏和肺均高表达HIF1A基因，该基因是缺氧诱导因子信号通路中的关键基因，调节下游与缺氧应激相关基因的表达，可能与水生适应有关。

关键词: 劳亚食虫类, 鼩鼱科, 转录组从头组装, 低氧适应

Abstract: China harbors high species diversity of Soricid shrews.Sixty-one species representing 12 genera are distributed in China,inhabiting various habitats.In this study,we sequenced the transcriptome of heart and lung tissue from three species,namely, Episoriculus umbrinus,Nectogale elegans,and Sorex bedfordiae,and obtained 20.5 Gb clean data.After de novo assembly,we obtained 131 045,153 621,and 135 838 unigenes for E.umbrinus,N.elegans,and S.bedfordiae,respectively.BUSCO analyses suggested we recovered as many as 82.2% full-length transcripts of genes.Ortho Finder recognized more than 8 000 homologous gene families as shared by soricid and talpid species.We compared two post-Trinity pipelines to optimize the quality of unigenes.When combining Trinity scripts,CD-Hit and TransRate,the proportions of single-copy genes increased 18%-25% compared with using CD-Hit alone,and the proportions of multi-copy genes decreased from 28.4%-31.8% to 0.7%-1.1%.Analysis of differentially expressed genes (DEGs) identified tissue enriched genes,validating our analysis pipeline.Our analyses HIF 1 A highly expressed in both heart and lung of N.elegans.HIF 1 A is a key gene in the hypoxia-inducible factor signaling pathway and regulates the expression of downstream genes related to hypoxia,probably associated with aquatic adaptation.This is the first report of shrew transcriptome,providing important data for understanding the evolution and physiology in Soricidae.Our analysis pipeline could also inspire transcriptomic study in other non-model organisms.

Key words: Eulipotyphla, Soricidae, Transcriptome de novo assembly, Adaptation to hypoxia

中图分类号:

周婧晶, 丘银彬, 万韬, 何水旺, 蒋学龙, 潘星华, 何锴. 鼩鼱科三物种转录组初步分析[J]. 兽类学报, 2021, 41(2): 173-181.

ZHOU Jingjing, QIU Yinbin, WAN Tao, HE Shuiwang, JIANG Xuelong, PAN Xinghua, HE Kai. Transcriptome analyses for three Soricid shrews[J]. ACTA Theriologica Sinica, 2021, 41(2): 173-181.

参考文献

Andrews S.2010.FastQC:A quality control tool for high throughput sequence data.In:Babraham Bioinformatics,Babraham Institute.Cambridge,United Kingdom.
Bolger A M,Lohse M,Usadel B.2014.Trimmomatic:A flexible trimmer for Illumina sequence data.Bioinformatics,30:2114-2120.
Burgin C,He K.2018.Family Soricidae (Shrews).Barcelona:Lynx Edicions,332-551.
Bushnell B.2014.BBMap:A fast,accurate,splice-aware aligner.In:Lawrence Berkeley National Lab.(LBNL),Berkeley,CA (United States).
Camacho C,Coulouris G,Avagyan V,Ma N,Papadopoulos J,Bealer K,Madden T L.2009.BLAST+:Architecture and applications.BMC Bioinformatics,10:421-421.
Dress A W M,Flamm C,Fritzsch G,Grünewald S,Kruspe M,Prohaska S J,Stadler P F.2008.Noisy:Identification of problematic columns in multiple sequence alignments.Algorithms for Molecular Biology,3:7.
Emms D M,Kelly S.2015.OrthoFinder:solving fundamental biases in whole genome comparisons dramatically improves orthogroup inference accuracy.Genome Biology,16:157.
Fan H,Xiao Y,Yang Y,Xia W,Mason A S,Xia Z,Qiao F,Zhao S,Tang H.2013.RNA-Seq analysis of Cocos nucifera:transcriptome sequencing and de novo assembly for subsequent functional genomics approaches.PLoS ONE,8:e59997.
Grabherr M G,Haas B J,Yassour M,Levin J Z,Thompson D A,Amit I,Adiconis X,Fan L,Raychowdhury R,Zeng Q,Chen Z,Mauceli E,Hacohen N,Gnirke A,Rhind N,di Palma F,Birren B W,Nusbaum C,Lindblad-Toh K,Friedman N,Regev A.2011.Full-length transcriptome assembly from RNA-Seq data without a reference genome.Nature Biotechnology,29:644-652.
Haas B J,Papanicolaou A,Yassour M,Grabherr M,Blood P D,Bowden J,Couger M B,Eccles D,Li B,Lieber M,MacManes M D,Ott M,Orvis J,Pochet N,Strozzi F,Weeks N,Westerman R,William T,Dewey C N,Henschel R,LeDuc R D,Friedman N,Regev A.2013.De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis.Nature Protocols,8:1494-1512.
Hutterer R.1985.Anatomical adaptations of shrews.Mammal Review,15:43-55.
Jain A,Tuteja G.2019.TissueEnrich:Tissue-specific gene enrichment analysis.Bioinformatics,35:1966-1967.
Joshi S,Singh A R,Durden D L.2014.MDM2 regulates hypoxic hypoxia-inducible factor 1α stability in an E3 ligase,proteasome,and PTEN-phosphatidylinositol 3-kinase-AKT-dependent manner.The Journal of Biological Chemistry,289:22785-22797.
Katoh K,Misawa K,Kuma K,Miyata T.2002.MAFFT:A novel method for rapid multiple sequence alignment based on fast Fourier transform.Nucleic Acids Research,30:3059-3066.
Langmead B,Salzberg S L.2012.Fast gapped-read alignment with Bowtie2.Nature Methods,9:357-359.
Li B,Dewey C N.2011.RSEM:Accurate transcript quantification from RNA-Seq data with or without a reference genome.BMC Bioinformatics,12:323-323.
Li W,Godzik A.2006.Cd-hit:a fast program for clustering and comparing large sets of protein or nucleotide sequences.Bioinformatics,22:1658-1659.
Liang H,He S,Yang J,Jia X,Wang P,Chen X,Zhang Z,Zou X,McNutt M A,Shen W H,Yin Y.2014.PTENα,a PTEN isoform translated through alternative initiation,regulates mitochondrial function and energy metabolism.Cell Metabolism,19:836-848.
Moreno-Santillan D D,Machain-Williams C,Hernandez-Montes G,Ortega J.2019.De Novo Transcriptome Assembly and Functional Annotation in Five Species of Bats.Scientific Reports,9.
Nishimura O,Hara Y,Kuraku S.2017.gVolante for standardizing completeness assessment of genome and transcriptome assemblies.Bioinformatics,33:3635-3637.
Nowak R M.1999.Walker's Mammals of the World.Baltimore:Johns Hopkins University Press
Prakash A,Jeffryes M,Bateman A,Finn R.2017.The HMMER web server for protein sequence similarity search.Current Protocols in Bioinformatics,60:3.15.11-13.15.23.
Robinson M D,McCarthy D J,Smyth G K.2010.edgeR:A Bioconductor package for differential expression analysis of digital gene expression data.Bioinformatics,26:139-140.
Semenza G L.2012.Hypoxia-inducible factors in physiology and medicine.Cell,148:399-408.
Simão F A,Waterhouse R M,Ioannidis P,Kriventseva E V,Zdobnov E M.2015.BUSCO:Assessing genome assembly and annotation completeness with single-copy orthologs.Bioinformatics,31:3210-3212.
Smith-Unna R,Boursnell C,Patro R,Hibberd J M,Kelly S.2016.TransRate:Reference-free quality assessment of de novo transcriptome assemblies.Genome Research,26:1134-1144.
Smith A T,Xie Y.2008.A Guide to the Mammals of China.Princeton:Princeton University Press
Stamatakis A.2014.RAxML version 8:a tool for phylogenetic analysis and post-analysis of large phylogenies.Bioinformatics,30:1312-1313.
Takagi H,King G L,Ferrara N,Aiello L P.1996.Hypoxia regulates vascular endothelial growth factor receptor KDR/Flk gene expression through adenosine A2 receptors in retinal capillary endothelial cells.Investigative Ophthalmology Visual Science,37:1311-1321.
Warnes G R,Bolker B,Bonebakker L,Gentleman R,Huber W,Liaw A,Lumley T,Maechler M,Magnusson A,Moeller S.2009.gplots:Various R programming tools for plotting data.R Package Version,2:1.
Ye J,Fang L,Zheng H,Zhang Y,Chen J,Zhang Z,Wang J,Li S,Li R,Bolund L.2006.WEGO:a web tool for plotting GO annotations.Nucleic Acids Research,34:W293-W297.
Zhao Q Y,Wang Y,Kong Y M,Luo D,Li X,Hao P.2011.Optimizing de novo transcriptome assembly from short-read RNA-Seq data:a comparative study.BMC Bioinformatics,12:S2-S2.
刘少英,吴毅.2019.中国兽类图鉴.福州市:海峡书局.

[1]	傅芳官久强曲秀龙王利罗晓林安添午张翔飞. 不同海拔饲育的麦洼牦牛肺脏组织转录组学分析[J]. 兽类学报, 2020, 40(5): 475-484.
[2]	刘铸张隽晟白薇刘欢解瑞雪杨茜金志民. 中国东北地区鼩鼱科动物分类与分布[J]. 兽类学报, 2019, 39(1): 8-26.
[3]	王琼王露权梦一汤燕燕韦海雪付长坤张修月陈贵英陈顺德宗浩. 鼩鼱科六个物种背部针毛的微观结构及其分类学意义[J]. 兽类学报, 2014, 34(4): 399-.
[4]	邓小弓王堃张守栋苏建平张同作林恭华. 利用转录组数据分析高原鼢鼠和裸鼹鼠基因的趋同进化[J]. 兽类学报, 2014, 34(2): 129-137.
[5]	刘芳, 格日力. 藏羚羊组蛋白去乙酰化酶1 基因编码区的克隆与序列分析[J]. , 2011, 31(4): 396-403.
[6]	涂飞云, 刘少英, 刘洋, 孙治宇, 朱天辉. 鼩鼱科内6 个种的阴茎形态学[J]. , 2010, 30(3): 278-282.

鼩鼱科三物种转录组初步分析

Transcriptome analyses for three Soricid shrews

PDF

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 6

编辑推荐

Metrics