nav emailalert searchbtn searchbox tablepage yinyongbenwen piczone journalimg journalInfo journalinfonormal searchdiv searchzone qikanlogo popupnotification paper paperNew
2024, 06, v.43 1039-1053
基于RNA-seq技术的不同发育期中华大蟾蜍耳后腺和腹皮转录组比较分析
基金项目(Foundation): 山东省重点研发(乡村振兴提振计划)项目(2022TZXD0035); 日照市重点研发计划项目(2021ZDYF010118)共同资助
邮箱(Email): hujinghong97@126.com;zhangjin78527@126.com;
DOI: 10.13417/j.gab.043.001039
摘要:

通过对中华大蟾蜍(Bufo bufo gargarizans)5个不同发育时期的耳后腺和腹皮进行转录组测序分析,挖掘参与毒素合成的关键因子,为深入研究毒素合成调控机理提供理论依据。分别取3月龄、 6月龄、 12月龄、 24月龄、 36月龄等5个发育时期的中华大蟾蜍耳后腺和腹皮进行转录组测序,利用GO和KEGG对差异表达基因功能进行分析,利用加权共表达网络分析筛选关键基因。将处于相同发育时期腹皮与耳后腺的转录组数据进行比较,分别筛选出3 053、 3 026、 1 516、 1 028和2 061个差异表达基因。由GO分类统计结果可知,差异表达基因在生物过程类别主要集中于细胞过程、单一生物体过程、代谢过程和生物调节;在分子功能类别主要富集于结合、催化活性等方面;在细胞组分类别主要聚集于细胞、细胞组分和细胞器。KEGG通路分析结果表明,差异表达基因在蟾蜍发育过程中参与多条通路,如酪氨酸代谢、过氧化物酶体和初级胆汁酸代谢等。将5个发育时期的共同差异表达基因进行比较,筛选出与类固醇代谢、初级胆汁酸代谢和过氧化物酶体相关的21个基因,利用加权共表达网络分析筛选出4个模块的6个关键基因。在不同发育时期中华大蟾蜍的腹皮与耳后腺差异表达基因中,大量代谢途径基因发生了富集,这些基因可能在调控耳后腺毒素合成中具有重要作用,为解析蟾蜍产毒分子机制奠定了基础。

Abstract:

To explore the key factors involved in toxin synthesis by the transcriptome sequencing analysis of the retroauricular glands and abdominal skin of Bufo bufo gargarizans at five different developmental stages, it provides a theoretical basis for further study on the regulation mechanism of toxin synthesis. Transcriptome sequencing was performed on the abdominal skin and retroauricular gland of Bufo bufo gargarizans at 5 developmental stages, including 3-month-old, 6-month-old, 12-month-old, 24-month-old, and 36-month-old. The function of differentially expressed genes was analyzed by GO and KEGG. The key genes were screened by weighted co-expression network analysis. The transcriptome data of the abdominal skin and retroauricular gland at the same developmental stage were compared, and 3053, 3 026, 1 516, 1 028, and 2 061 differentially expressed genes were screened, respectively. According to the GO classification statistics, the differentially expressed genes were mainly concentrated in biological processes categories such as cellular processes, monomer processes, metabolic processes, and biological regulations. In the molecular function categories, it is mainly concentrated in the binding, catalytic activity, and so on. In cellular component categories, it is mainly concentrated in the cell, component and organelle. The analysis of the KEGG pathway showed that the differentially expressed genes were involved in many pathways during toad development, such as tyrosine metabolism, peroxisome metabolism, and primary bile acid metabolism. Comparing common differentially expressed genes from 5 developmental stages, twenty-one genes related to steroid metabolism, primary bile acid metabolism, and the peroxisome were identified. Five hub genes and one key transcription factor from four modules were screened by weighted co-expression network analysis. A large number of genes of the metabolic pathway were enriched in the differentially expressed genes between the skin and retroauricular glands of Bufo bufo gargarizans at different developmental stages. These genes may play an important role in regulating the synthesis of retroauricular adenotoxins. This study lays the foundation for understanding the molecular mechanism of toad toxin production.

参考文献

陈文杰,胡耀廷,王如锋,等,2021.蟾蜍二烯酸内酯类成分的生物转化和生物合成研究进展.中草药,52(15):4741-4751.[CHEN W J,HU Y T,WANG R F,et al.,2021.Research progress on biotransformation and biosynthesis of bufadienolides.Chinese Traditional and Herbal Drugs,52(15):4741-4751.]

陈瀛澜,郝艳艳,郭夫江,等,2017.蟾酥化学成分及药理活性研究进展.中草药,48(12):2579-2588.[CHEN Y L,HAO Y Y,GUO F J,et al.,2017.Research progress on chemical constituents and pharmacological activities of Bufonis Venenum.Chinese Traditional and Herbal Drugs,48(12):2579-2588.]

胥文才,蔺杼阳,肖雅丹,等,2020.中华蟾蜍多组织差异转录组分析.四川动物,39(6):601-615.[XU W C,LIN Z Y,XIAO Y D,et al.,2020.Differential transcriptome analysis in multi-tissue of Bufo gargarizans.Sichuan Journal of Zoology,39(6):601-615.]

于海滨,2020.线粒体甘油-3-磷酸乙酰转移酶(GPAM)基因沉默对牛乳脂代谢的影响,博士学位论文.长春:吉林大学.[YU H B,2020.The effect of Mitochondrial glycerol-3-phosphate acyltransferase(GPAM) gene silencing on milk fat metabolism in cows,Dissertation for Ph.D.Changchun:Jilin University.]

ASHBURNER M,BALL C A,BLAKE J A,et al.,2020.Gene ontology:tool for the unification of biology.Nat.Genet.,25(1):25-29.

BLACKBURN A C,COGGAN M,SHIELD A J,et al.,2011.Glutathione transferase kappa deficiency causes glomerular nephropathy without overt oxidative stress.Lab.Investig.,91(11):1572-1583.

BRUZZONE A,BUSCHIAZZO J,ALONSO T S,2003.Lipids during Bufo arenarum oogenesis.Zygote,11(2):95-100.

DING G F,LI Y,GAO J R,et al.,2022.Associations between AGT,MTHFR,and VEGF gene polymorphisms and pree-clampsia in the Chinese population.Placenta,118:38-45.

FAN J J,LI X L,ISSOP L,et al.,2016.ACBD2/ECI2-mediated peroxisome-mitochondria interactions in leydig cell ste-roid biosynthesis.Mol.Endocrinol.,30(7):763-782.

FERDINANDUSSE S,ZOMER A W M,KOMEN J C,et al.,2008.Ataxia with loss of Purkinje cells in a mouse model for Refsum disease.Proc.Natl.Acad.Sci.USA,105(46):17712-17717.

FIRNENO T J,RAMESH B,MALDONADO J A,et al.,2022.Transcriptomic analysis reveals potential candidate pathways and genes involved in toxin biosynthesis in true toads.J.Hered.,113(3):311-324.

FRACCASCIA P,CASTEELS M,DE SCHRYVER E,et al.,2011.Role of thiamine pyrophosphate in oligomerisation,functioning and import of peroxisomal 2-hydroxyacyl-CoA lyase.Biochim.Biophys.Acta (BBA) Proteins Proteom.,1814(10):1226-1233.

FURUSAWA Y,KUBO T,FUKAZAWA T,2016.Phyhd1,an XPhyH-like homologue,is induced in mouse T cells upon T cell stimulation.Biochem.Biophys.Res.Commun.,472(3):551-556.

GUO P P,JIN X,ZHANG J F,et al.,2023.Overexpression of DGAT2 regulates the differentiation of bovine preadipocytes.Animals,13(7):1195.

HOUTEN S M,DENIS S,ARGMANN C A,et al.,2012.Peroxisomal L-bifunctional enzyme (Ehhadh) is essential for the production of medium-chain dicarboxylic acids.J.Lipid Res.,53(7):1296-1303.

KANEHISA M,GOTO S,2000.KEGG:Kyoto encyclopedia of genes and genomes.Nucleic.Acids Res.,28(1):27-30.

LAN Y,HE L W,DONG X,et al.,2022.Comparative transcriptomes of three different skin sites for the Asiatic toad (Bufo gargarizans).PeerJ,10:e12993.

LANGFELDER P,HORVATH S,2008.WGCNA:an R package for weighted correlation network analysis.BMC Bioinforma-tics,9:559.

LI B,DEWEY C N,2011.RSEM:accurate transcript quantification from RNA-seq data with or without a reference genome.BMC Bioinformatics,12:323.

LI F J,HU J H,REN X,et al.,2021.Toad venom:a comprehensive review of chemical constituents,anticancer activities,and mechanisms.Arch.der Pharm.,354(7):e2100060.

LIANG Z J,CHEN Y B,WANG L P,et al.,2019.CYP27A1 inhibits bladder cancer cells proliferation by regulating cholesterol homeostasis.Cell Cycle,18(1):34-45.

LIU J S,DENG L J,TIAN H Y,et al.,2019.Anti-tumor effects and 3D-quantitative structure-activity relationship analysis of bufadienolides from toad venom.Fitoterapia,134:362-371.

LOVE M I,HUBER W,ANDERS S,2014.Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2.Genome Biol.,15(12):1-21.

MA H Y,KOU J P,WANG J R,et al.,2007.Evaluation of the anti-inflammatory and analgesic activities of Liu-Shen-Wan and its individual fractions.J.Ethnopharmacol.,112(1):108-114.

MAILHO-FONTANA P L,ANTONIAZZI M M,SCIANI J M,et al.,2018.Morphological and biochemical characterization of the cutaneous poison glands in toads (Rhinella marina group) from different environments.Front.Zool.,15(1):1-15.

MARTISKAINEN H,PALDANIUS K M A,NATUNEN T,et al.,2017.DHCR24 exerts neuroprotection upon inflammation-induced neuronal death.J.Neuroinflammation,14(1):215.

NIKOLAOU N,GATHERCOLE L L,KIRKWOOD L,et al.,2019.AKR1D1 regulates glucocorticoid availability and glucocorticoid receptor activation in human hepatoma cells.J.Steroid Biochem.Mol.Biol.,189:218-227.

OKAMURA S,YOSHINO H,KUROSHIMA K,et al.,2021.EHHADH contributes to cisplatin resistance through regulation by tumor-suppressive microRNAs in bladder cancer.BMC Cancer,21(1):48.

SHANNON P,MARKIEL A,OZIER O,et al.,2003.Cytoscape:a software environment for integrated models of biomolecular interaction networks.Genome Res.,13(11):2498-2504.

SIPERSTEIN M D,MURRAY A W,TITUS E,1957.Biosynthesis of cardiotonic sterols from cholesterol in the toad,Bufo marinus.Arch.Biochem.Biophys.,67(1):154-160.

SOUMOY L,GHANEM G E,SAUSSEZ S,et al.,2022.Bufalin for an innovative therapeutic approach against cancer.Pharmacol.Res.,184:106442.

SUI J S Y,MARTIN P,KEOGH A,et al.,2022.Altered expression of ACOX2 in non-small cell lung cancer.BMC Pulm.Med.,22(1):1-21.

TóTH Z,KURALI A,MóRICZ á M,et al.,2019.Changes in toxin quantities following experimental manipulation of toxin reserves in Bufo bufo tadpoles.J.Chem.Ecol.,45(3):253-263.

WANG T,LIU C,XIONG Y Z,et al.,2007.Isolation and cloning of porcine SLC27A2 gene and detection of its polymorphism associated with growth and carcass traits.Asian Australas.J.Anim.Sci,20(8):1169-1173.

ZHENG Y R,DENG L J,CAO H H,et al.,2022.Screening of bufadienolides from toad venom identifies gammabufotalin as a potential anti-inflammatory agent.Planta Med.,88(1):43-52.

基本信息:

DOI:10.13417/j.gab.043.001039

中图分类号:S865.49

引用信息:

[1]陈文潇,陈晶,郭媛媛,等.基于RNA-seq技术的不同发育期中华大蟾蜍耳后腺和腹皮转录组比较分析[J].基因组学与应用生物学,2024,43(06):1039-1053.DOI:10.13417/j.gab.043.001039.

基金信息:

山东省重点研发(乡村振兴提振计划)项目(2022TZXD0035); 日照市重点研发计划项目(2021ZDYF010118)共同资助

发布时间:

2024-01-22

出版时间:

2024-01-22

网络发布时间:

2024-01-22

检 索 高级检索