基于高通量转录组测序的山羊睾丸和附睾头差异表达基因分析

doi:10.11843/j.issn.0366-6964.2014.03.008

摘要/Abstract

摘要：

为探明山羊睾丸和附睾差异表达基因，挖掘与精子成熟和贮存密切相关的功能基因。采用Illumina HiSeq^TM 2000高通量转录组测序RNA-seq技术对3只成年种公羊睾丸和附睾头进行测序，将组装得到的Unigene比对到Nr、GO和KEGG数据库中注释，并进行差异表达基因聚类分析。结果显示，睾丸和附睾头Unigene总数分别为57 727个和51 395个， GO功能分类注释到生物过程、分子功能和细胞组成数据库中的分别有25 201、25 456和27 318个，其中差异表达基因数分别为9 150、9 172和10 003个；注释到KEGG通路中所有基因数23 222个，其中差异表达基因数8 219个，差异基因显著富集的信号通路：mRNA监视通路和mTOR信号通路；睾丸和附睾头中差异表达的基因共有24 592个，其中上调基因有8 420个，下调基因有16 172个。在附睾头中上调幅度最大的前3位：gDEFB124、gDEFB109和gDEFB108，下调幅度最大的前3位：GTSF1L、SRRM1和FNDC8。附睾头中表达量最高的前3位基因依次为Lcn5、gDEFB124和GPx5。睾丸中表达的基因数量大于附睾头中数量；β防御素家族、Lipocalin家族和GPx5可能在附睾头精子成熟、贮存和受精中发挥重要作用。

关键词: 转录组, 差异表达基因, 睾丸, 附睾头, 公山羊

Abstract:

The objective of this study was to find the differentially expressed genes between testis and epididymis caput in adult bucks，which was necessary to better understand mechanisms involved in spermatogenesis and sperm maturation.Testis and epididymis caput of three adult bucks known fertility were isolated and their mRNA were respectively extracted.The transcriptomes were analyzed by Illumina HiSeq^TM 2000 high-throughput RNA sequencing.All of Unigenes were annotated using the BLAST search against Nr，GO and KEGG databases.The differentially expressed genes were enriched in gene ontology and KEGG pathway.The results showed that 57 727 and 51 395 Unigenes were respectively in caprine testis and epididymis caput; There were 25 201 Unigenes annotated into database of molecular function，25 456 Unigenes into cellular component，27 318 Unigenes into biological process，of which there were 9 150，9 172 and 10 003 the differentially expressed genes respectively; There were 23 222 Unigenes annotated into KEGG pathways，including 8 219 the differentially expressed genes.The signal pathways that differentially expressed genes clustered significantly were mRNA surveillance pathway and mTOR signaling pathway; There were 24 592 the differentially expressed genes in caprine testis and epididymis caput，in which 8 420 genes were up-regulated and 16 172 genes were down-regulated.The top three genes of all up-regulated genes in epididymis caput were beta-defensin 124、beta-defensin 109 and beta-defensin 108，whereas the top three genes of all down-regulated genes were gametocyte specific factor 1-like，Serine/arginine repetitive matrix 1 and fibronectin type III domain containing 8.The three of the highest expressed genes in epididymis caput were lipocalin 5，beta-defensin 124 and glutathione peroxidase 5.The transcriptomes in testis were more than those in epididymis caput; The family of bete-defensins，lipocalins and GPx5 may play important roles in sperm maturation，storage and fertilization.

Key words: transcriptomes, differentially expressed genes, testis, epididymis caput, buck

中图分类号:

S827.2

张春香，张国林，郭丽娜，赵辉，任有蛇. 基于高通量转录组测序的山羊睾丸和附睾头差异表达基因分析[J]. 畜牧兽医学报, 2014, 45(3): 391-401.

ZHANG Chun-xiang，ZHANG Guo-lin，GUO Li-na，ZHAO Hui，REN You-she. Study on Differentially Expressed Genes between Caprine Testis and Epididymis Caput Based on Transcriptomes with High-throughput RNA-seq Technology[J]. ACTA VETERINARIA ET ZOOTECHNICA SINICA, 2014, 45(3): 391-401.

参考文献

［1］GATTI J L，CASTELLA S，DACHEUX F，et al.Post-testicular sperm environment and fertility ［J］.Anim Reprod Sci，2004，82-83:321-339.
［2］DACHEUX J L，BELLEANNÉE C，JONES R，et al.Mammalian epididymal proteome ［J］.Mol Cell Endocrinol，2009，306(1-2)：45-50.
［3］RAMSKÖD，WANG E T，BURGE C B，et al.An abundance of ubiquitously expressed genes revealed by tissue transcriptome sequence data ［J］.PLoS Comput Biol，2009，5(12)：e1000598.
［4］ANAND M，PRASAD B V.The computational analysis of human testis transcriptome reveals closer ties to pluripotency ［J］.J Hum Reprod Sci，2012，5(3):266-273.
［5］ZHANG J S，LIU Q，LI Y M，et al.Genome-wide profiling of segmental-regulated transcriptomes in human epididymis using oligo microarray ［J］.Mol Cell Endocrinol，2006，250(1-2)：169-177.
［6］LI J Y，WANG H Y，LIU J，et al.Transcriptome analysis of a cDNA library from adult human epididymis ［J］.DNA Res，2008，15(3)：115-122.
［7］JOHNSON G D，LALANCETTE C，LINNEMANN A K，et al.The sperm nucleus：chromatin，RNA，and the nuclear matrix ［J］.Reproduction, 2011，141(1)：21-36.
［8］SOUMILLON M，NECSULEA A，WEIER M，et al.Cellular source and mechanisms of high transcriptome complexity in the Mammalian testis ［J］.Cell Rep，2013，3(6):2179-2190.
［9］JOHNSTON D S，JELINSKY S A，BANG H J，et al.The mouse epididymal transcriptome：transcriptional profiling of segmental gene expression in the epididymis ［J］. Biol Reprod，2005，73(3):404-413.
［10］JELINSKY S A，TURNER T T，BANG H J，et al.The rat epididymal transcriptome：comparison of segmental gene expression in the rat and mouse epididymides ［J］.Biol Reprod，2007，76(4):561-570.
［11］GUYONNET B，MAROT G，DACHEUX J L，et al.The adult boar testicular and epididymal transcriptomes ［J］.BMC Genom，2009,10：369.
［12］GRABHERR M G，HAAS B J，YASSOUR M，et al.Full-length transcriptome assembly from RNA-Seq data without a reference genome ［J］.Nat Biotechnol，2011，29(7)：644-652.
［13］LI R，YU C，LI Y R，et al.SOAP2：An improved ultrafast tool for short read alignment ［J］.Bioinformatics，2009，25(15)：1966-1967.
［14］MORTAZAVI A，WILLIAMS B A，MCCUE K，et al.Mapping and quantifying mammalian transcriptomes by RNA-Seq ［J］.Nature Methods，2008，5：621-628.
［15］YE J，FANG L，ZHENG H，et al.WEGO：a web tool for plotting GO annotations ［J］.Nucleic Acids Res，2006，34(Web Server issue):W293-297.
［16］SOUMILLON M，NECSULEA A，WEIER M，et al.Cellular source and mechanisms of high transcriptome complexity in the mammalian testis［J］.Cell Report，2013,3:2179-2190.
［17］ JHA K N，SHUMILIN I A，DIGILIO L C，et al.Biochemical and structural characterization of apolipoprotein A-I binding protein，a novel phosphoprotein with a potential role in sperm capacitation ［J］.Endocrinology，2008,149(5):2108-2120.
［18］HARAGUCHI C M，ISHIDO K，KOMINAMI E，et al.Expression of cathepsin H in differentiating rat spermatids：immunoelectron microscopic study［J］.Histochem Cell Biol, 2003,120(1):63-71.
［19］PASTÉN H K，HERNÁNDE R R，ROA A L，et al.Presence，processing，and localization of mouse ADAM15 during sperm maturation and the role of its disintegrin domain during sperm-egg binding ［J］.Reproduction,2008,136(1):41-51.
［20］TANWAR P S，KANEKO-TARUI T，ZHANG L，et al.Altered LKB1/AMPK/TSC1/TSC2/mTOR signaling causes disruption of Sertoli cell polarity and spermatogenesis ［J］.Hum Mol Genet，2012，21(20):4394-4405.
［21］SRIVASTAVA S，DESAI P，COUTINHO E，et al.Mechanism of action of L-arginine on the vitality of spermatozoa is primarily through increased biosynthesis of nitric oxide［J］.Biol Reprod，2006,74(5):954-958.
［22］PATIL A A，CAI Y，SANG Y，et al.Cross-species analysis of the mammalian β-defensin gene expression in the male reproductive tract family：presence of syntenic gene clusters and preferential expression in the male reproductive tract ［J］.Physiol Genomics，2005，23(1):5-17.
［23］HALL S H，YENUGU S，RADHAKRISHNAN Y，et al.Characterization and functions of beta defensins in the epididymis ［J］.Asian J Andro，2007，9(4)：453-462.
［24］YUDIN A I，TREECE C A，TOLLNER T L，et al.The carbohydrate structure of DEFB126，the major component of the cynomolgus macaque sperm plasma membrane glycocalyx ［J］.Membr Biol,2005，207(3):119-129.
［25］YUDIN A I，TOLLNER T L，TREECE C A，et al.Beta -defensin 22 is a major component of the mouse sperm glycocalyx ［J］.Reproduction，2008，136(6)：753-765.
［26］NARCIANDI F，LLOYD A T，CHAPWANYA A，et al.Reproductive tissue-specific expression profiling and genetic variation across a 19 gene bovine β-defensin cluster ［J］.Immungen，2011，63(10)：641-651.
［27］CHOI M K，LE M T，NGUYEN D T，et al.Genome-level identification，gene expression，andcomparative analysis of porcine β-defensin genes ［J］.BMC Genetics，2012，13：98.
［28］SALIER J P，AKERSTRÖM B，BORREGAARD N，et al.Lipocalins in bioscience：the first family gathering ［J］.Bioessays，2004，26(4):456-458.
［29］SUZUKI K，LAREYRE J J，SÁNCHEZ D，et al.Molecular evolution of epididymal lipocalin genes localized on mouse chromosome 2 ［J］.Gene，2004，339：49-59.
［30］张建新，王茜，荀文娟，等.GPx家族基因在山羊不同组织和睾丸不同发育时期的表达特性研究［J］.畜牧兽医学报，2011，42(5)：650-657.
［31］IMAI H，HAKKAKU N，IWAMOTO R，et al.Depletion of selenoprotein GPx4 in spermatocytes causes male infertility in mice［J］.J Biol Chem, 2009，284 (47)：32522-32532.
［32］CHABORY E，DAMON C，LENOIR A L，et al.Epididymis seleno-independent glutathione peroxidase 5 maintains sperm DNA integrity in mice ［J］.J Clin Invest，2009，119：2074-2085.

[1]	王璟, 滑留帅, 陈俊峰, 张家庆, 任巧玲, 白红杰, 郭红霞, 徐照学, 邢宝松, 白献晓, 曹海. 去势对淮南公猪背最长肌转录组的影响[J]. 畜牧兽医学报, 2019, 50(9): 1746-1758.
[2]	荣敏, 张然然, 涂剑峰, 王天娇, 孙丽敏, 杨颖, 徐佳萍, 邢秀梅*. 美国短毛黑水貂快速生长过程中转录组差异表达分析[J]. 畜牧兽医学报, 2019, 50(9): 1787-1801.
[3]	赵学亮, 王姝懿, 孙柯, 苏倩, 王文龙, 刘春霞. 捻转血矛线虫阿苯达唑敏感株和耐药株比较转录组学分析[J]. 畜牧兽医学报, 2019, 50(9): 1940-1944.
[4]	施会彬, 王玉琴, 张小辉, 吴至博, 杨芳. 基于转录组测序数据的山羊下丘脑组织繁殖相关信号通路及候选基因的筛选[J]. 畜牧兽医学报, 2019, 50(7): 1380-1395.
[5]	施会彬, 王玉琴, 张小辉, 吴至博, 杨芳. 基于转录组测序数据的山羊下丘脑组织繁殖相关信号通路及候选基因的筛选[J]. 畜牧兽医学报, 2019, 50(7): 1380-1395.
[6]	张石磊, 包佳鹭, 史万玉, 钟秀会. 母鼠孕期和哺乳期连续染毒双酚A对子代雄鼠睾丸发育的影响[J]. 畜牧兽医学报, 2019, 50(5): 993-1005.
[7]	张逢, 陈指龙, 易思亮, 许璐瑶, 方娟, 黎陈, 傅胜才, 杨青. 细菌脂多糖急性暴露通过炎症反应诱导小鼠睾丸损伤[J]. 畜牧兽医学报, 2019, 50(5): 1099-1105.
[8]	张逢, 陈指龙, 易思亮, 许璐瑶, 方娟, 黎陈, 傅胜才, 杨青. 细菌脂多糖急性暴露通过炎症反应诱导小鼠睾丸损伤[J]. 畜牧兽医学报, 2019, 50(5): 1099-1105.
[9]	张石磊, 包佳鹭, 史万玉, 钟秀会. 母鼠孕期和哺乳期连续染毒双酚A对子代雄鼠睾丸发育的影响[J]. 畜牧兽医学报, 2019, 50(5): 993-1005.
[10]	刘丽娟, 董然然, 王开功, 程振涛, 文明, 温贵兰, 李晨, 杨琦, 周碧君. 敲除小RNA gcvB后沙门菌的转录组分析[J]. 畜牧兽医学报, 2019, 50(4): 840-850.
[11]	袁莉刚, 宋笑, 杨婕, 梁霞霞, 陈国娟. TGF-β1及其受体在性成熟前后藏绵羊睾丸的表达比较[J]. 畜牧兽医学报, 2019, 50(4): 750-757.
[12]	袁莉刚, 宋笑, 杨婕, 梁霞霞, 陈国娟. TGF-β1及其受体在性成熟前后藏绵羊睾丸的表达比较[J]. 畜牧兽医学报, 2019, 50(4): 750-757.
[13]	刘丽娟, 董然然, 王开功, 程振涛, 文明, 温贵兰, 李晨, 杨琦, 周碧君. 敲除小RNA gcvB后沙门菌的转录组分析[J]. 畜牧兽医学报, 2019, 50(4): 840-850.
[14]	李欣, 李小俊, 陈晓丽, 赵毅强, 王栋. 转录组数据分析与功能基因挖掘[J]. 畜牧兽医学报, 2019, 50(3): 474-484.
[15]	阿果约达, 熊显荣, 王艳, 杨显英, 韩杰, 王斌, 李键. 犏牛雄性不育相关lncRNA的鉴定与分析[J]. 畜牧兽医学报, 2019, 50(3): 551-561.