研究论文

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

  • 周婧晶 ,
  • 丘银彬 ,
  • 万韬 ,
  • 何水旺 ,
  • 蒋学龙 ,
  • 潘星华 ,
  • 何锴
展开
  • 1 南方医科大学基础医学院, 广东省单细胞技术与应用重点实验室, 广州 510515;
    2 陕西师范大学生命科学学院, 西安 710119;
    3 中国科学院昆明动物研究所, 昆明 650223
周婧晶(1995-),女,硕士研究生,主要从事进化生物学研究.

收稿日期: 2019-12-30

  网络出版日期: 2021-04-02

基金资助

国家自然科学基金面上基金(31970389);中国博士后科学基金面上项目(2009M652952);广东省珠江人才计划本土创新科研团队(2017BT01S131)

Transcriptome analyses for three Soricid shrews

  • ZHOU Jingjing ,
  • QIU Yinbin ,
  • WAN Tao ,
  • HE Shuiwang ,
  • JIANG Xuelong ,
  • PAN Xinghua ,
  • HE Kai
Expand
  • 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 date: 2019-12-30

  Online published: 2021-04-02

摘要

中国是鼩鼱科物种十分丰富的国家,已知有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基因,该基因是缺氧诱导因子信号通路中的关键基因,调节下游与缺氧应激相关基因的表达,可能与水生适应有关。

本文引用格式

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

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.

参考文献

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.中国兽类图鉴.福州市:海峡书局.
文章导航

/

〈 〉