基于单细胞转录组数据探究不同海拔牦牛心肌细胞通讯和基因表达差异

doi:10.16829/j.slxb.150817

摘要/Abstract

摘要： 心脏作为血液循环最重要的器官在动物环境适应中发挥着重要的作用。心肌细胞是心脏主要的功能细胞，然而高海拔地区牦牛心肌细胞基因表达变化及细胞间通讯关系尚不清楚。为此，本研究以高海拔地区牦牛 (青海祁连，海拔 4 000 m) 和低海拔地区牦牛 (青海循化，海拔 2 600 m) 为研究对象，基于心脏组织 10 × 单细胞转录组测序数据，对心肌细胞和心脏其他细胞进行细胞通讯及配体-靶基因调控预测，同时对高、低海拔牦牛心肌细胞差异表达基因进行功能注释分析，以期探究心肌细胞在牦牛适应高原环境过程中的作用。结果显示，牦牛心肌细胞与内皮细胞、上皮细胞的关联最强，心肌细胞-树突状细胞及心肌细胞-巨噬细胞两个“细胞对”中CD74_COAP 和 CD74_APP 配受体表达量相对较高；配受体活力在 TNF 信号通路中相对较强；在免疫相关通路发挥重要作用的配体 PTPRC、PECAM1、ITGB2、ANXA1、BDNF 等对所有“细胞对”影响明显，且 PI3K-Akt 信号通路在配体-靶基因调控中潜力分值最高；高、低海拔牦牛心肌细胞中差异表达基因功能主要富集在代谢途径，其中氧化磷酸化、糖酵解代谢通路相关的基因表达明显增强。本研究结果提示，牦牛可通过加强心肌细胞与其他细胞之间的联系及调节免疫作用相关过程维持自身稳态；同时高海拔牦牛可通过增强心肌细胞有氧代谢和无氧代谢进而获取更多能量以适应高海拔严酷环境。

关键词: 牦牛, 心肌细胞, 单细胞转录组测序, 高原环境适应

Abstract: As the most important organ for blood circulation, the heart plays an important role in the environmental adaptation of animals. Cardiomyocyte is the main functional cell in the animal heart, however, gene expression changes and cellular communication in cardiomyocytes of yaks living at high altitudes are still unclear. To study the effect of yak cardiomyocytes in the process of adaptation to the plateau environment, yaks in high altitude areas (Qilian, Qinghai, altitude 4000 m) and low altitude areas (Xunhua, Qinghai, altitude 2 600 m) were taken as research objects, and the cellular communications and ligand-target gene regulations in cardiomyocytes and other cells were predicted based on the 10× singlecell RNA sequencing data. In addition, the functions of differentially expressed genes (DEGs) between different altitude yak cardiomyocytes were annotated. The results showed that yak cardiomyocytes were most highly correlated with the endothelial cells and epithelial cells. The expression level of ligand and receptor of CD74_COAP and CD74_APP was relatively higher in the cardiomyocyte-dendritic cell and cardiomyocyte-macrophage, the vitality of ligands and receptors was relatively higher in TNF signal pathway; The immune related pathway ligands, such as PTPRC, PECAM1, ITGB2, ANXA1, BDNF had a significant influence on all ‘cell pairs’and the PI3K-Akt signal pathway showed the highest potential score in ligand-target gene regulation; The functions of DEGs between high and low altitude yak cardiomyocytes were mainly enriched in metabolic pathways as the oxidative phosphorylation and glycolysis pathways. These results suggest that yaks may maintain their homeostasis by strengthening the connection between cardiomyocytes and other cells and regulating the processes related to immune action. Meanwhile, high-altitude yaks enhance their aerobic and anaerobic metabolism of cardiomyocytes to acquire more energy to adapt to the harsh environment of high altitude.

Key words: Yak, Cardiomyocytes, scRNA-seq, Plateau environment adaptation

中图分类号:

Q494

黄娜婷, 陈一博, 宋娜, 李静怡, 魏青, 陈家瑞. 基于单细胞转录组数据探究不同海拔牦牛心肌细胞通讯和基因表达差异[J]. 兽类学报, 2024, 44(2): 171-182.

HUANG Nating, CHEN Yibo, SONG Na, LI Jingyi, WEI Qing, CHEN Jiarui. Cardiomyocyte communication and its gene expression difference between different altitudes yaks based on single cell RNA-seq data[J]. ACTA THERIOLOGICA SINICA, 2024, 44(2): 171-182.

参考文献

Averdunk L, Rückbeil M V, Zarbock A, Martin L, Marx G, Jalaie H, Ja cobs M J, Stoppe C, Gombert A. 2020. SLPI-A biomarker of acute kidney injury after open and endovascular thoracoabdominal aortic aneurysm (TAAA) repair. Scientific Reports, 10 (1):3453.
Barashdi M A A, Ali A, Mcmullin M F, Mills K. 2021. Protein tyro sine phosphatase receptor type C (PTPRC or CD45). Journal of Clinical Pathology, 74 (9):548-552.
Browaeys R, Saelens W, Saeys Y. 2020. NicheNet:modeling intercel lular communication by linking ligands to target genes. Nature Methods, 17 (2):159-162.
Chen X, Chakravarty T, Zhang Y, Li X, Zhong J F, Wang C. 2016.Single-cell transcriptome and epigenomic reprogramming of cardiomyocyte-derived cardiac progenitor cells. Scientific Data, 3:160079.
Chen Y, Meng X, Wan R, Cheng R, Zhang G, Zhang Q, Yu H, Wei Q. 2023. Single-cell transcriptomic survey of cell diversity and functional changes in yak hearts at different altitude. Proteomics, 23(11):e2200345.
Cheng W L, Kao Y H, Chen Y C, Lin Y K, Chen S A, Chen Y J. 2020.Macrophage migration inhibitory factor increases atrial arrhyth mogenesis through CD74 signaling. Translation Research, 216:43-56.
Cornwell W K, Baggish A L, Bhatta Y K D, Brosnan M J, Dehnert C, Guseh J S, Hammer D, Levine B D, Parati G, Wolfel E E. 2021.
Clinical implications for exercise at altitude among individuals with cardiovascular disease:A scientific statement from the American heart association. Journal of the American Heart Asso ciation, 10 (19):e023225.
de Jong R, Leoni G, Drechsler M, Soehnlein O. 2017. The advantageous role of annexin A1 in cardiovascular disease. Cell Adhe sion & Migration, 11 (3):261-274.
Ding X Z, Liang C N, Guo X, Wu X Y, Wang H B, Johnson K A, Yan P. 2014. Physiological insight into the high-altitude adaptations in domesticated yaks (Bosgrunniens) along the Qinghai-Tibetan Plateau altitudinal gradient. Livestock Science, 162:233-239.
Efremova M, Vento-Tormo M, Teichmann S A, Vento-Tormo R. 2020. CellPhoneDB:inferring cell-cell communication from com bined expression of multi-subunit ligand-receptor complexes. Nature Protocols, 15 (4):1484-1506.
Falcon S, Gentleman R. 2007. Using GOstats to test gene lists for GO term association. Bioinformatics, 23 (2):257-258.
Guan J, Long K, Ma J, Zhang J, He D, Jin L, Tang Q, Jiang A, Wang X, Hu Y, Tian S, Jiang Z, Li M, Luo X. 2017. Comparative analysis of the microRNA transcriptome between yak and cattle provides insight into high-altitude adaptation. PeerJ, 5:e3959.
Guengerich F P. 2019. Cytochrome P450 research and The Journal of Biological Chemistry. Journal of Biological Chemistry, 294 (5):1671-1680.
Hang P Z, Zhu H, Li P F, Liu J, Ge F Q, Zhao J, Du Z M. 2021. The emerging role of BDNF/TrkB signaling in cardiovascular dis eases. Life (Basel), 11 (1). DOI:10. 3390/life11010070.
Jha P K, Vijay A, Sahu A, Ashraf M Z. 2016. Comprehensive gene expression meta-analysis and integrated bioinformatic approaches reveal shared signatures between thrombosis and myeloproliferative disorders. Scientific Reports, 6:37099.
Krishnan G, Paul V, Hanah S S, Bam J, Das P J. 2016. Effects of climate change on yak production at high altitude. The Indian Journal of Animal Sciences, 86 (6):621-626.
Kuang L, Zheng Y, Lin Y, Xu Y, Jin S, Li Y, Dong F, Jiang Z. 2010.
High-altitude adaptation of yak based on genetic variants and activity of lactate dehydrogenase-1. Biochemical Genetics, 48 (5-6):418-427.
Leng L, Metz C N, Fang Y, Xu J, Donnelly S, Baugh J, Delohery T, Chen Y, Mitchell R A, Bucala R. 2003. MIF signal transduction initiated by binding to CD74. Journal of Experimental Medicine, 197 (11):1467-1476.
Leone M, Engel F B. 2021. Isolation, culture, and live-cell imaging of primary rat cardiomyocytes. Methods in Molecular Biology, 2158:109-124.
Lertkiatmongkol P, Liao D, Mei H, Hu Y, Newman P J. 2016. Endothelial functions of platelet/endothelial cell adhesion molecule-1(CD31). Current Opinion in Hematology, 23 (3):253-259.
Li Q L, Tang J, Zhao L, Ruze A, Shan X F, Gao X M. 2023. The role of CD74 in cardiovascular disease. Frontiers in Cardiovascular Medicine, 9:1049143.
Li R, Chen S, Li C, Xiao H, Costa V, Bhuiyan M S A, Baig M, BejaPereira A. 2022. Whole-genome analysis deciphers population structure and genetic introgression among bovine species. Fron tiers in Genetics, 13:847492.
Li X, Gao Y H, Peng Z L, Guo C H, Wang Q, Liang, Wang D. 2022.Effect of altitude changes on blood cell physiological parameters and cardiopulmonary histomorphological characteristics of yaks.Journal of Domestic Animal Ecology, 43 (4):47-54. (in Chinese)
Li X, Wei L, Wang Y, Xu L N, Wei L N, Wei D B. 2015. The expression of the sperm-specific lactate dehydrogenase gene Ldh-c in plateau pika (Ochotona curzoniae) cardiac muscle and its effect on the anaerobic glycolysis. Sheng Li Xue Bao, 67 (3):312-318.(in Chinese)
Long K R. 2020. Comparative microRNA transcriptomic analysis provides insights into high-altitude acclimation across different species. Ph. D thesis. Ya'an:Sichuan Agricultural University. (in Chinese)
Lu Y, Wu Q, Liao J, Zhang S, Lu K, Yang S, Wu Y, Dong Q, Yuan J, Zhao N, Du Y. 2021. Identification of the distinctive role of DPT in dilated cardiomyopathy:a study based on bulk and single-cell transcriptomic analysis. Annals of Translational Medicine, 9 (18):1401.
Marsh L M, Cakarova L, Kwapiszewska G, von Wulffen W, Herold S, Seeger W, Lohmeyer J. 2009. Surface expression of CD74 by type II alveolar epithelial cells:a potential mechanism for macro phage migration inhibitory factor-induced epithelial repair.American Journal of Physiology-Lung Cellular and Molecular Physiology, 296 (3):L442-L452.
Mason F E, Pronto J R D, Alhussini K, Maack C, Voigt N. 2020. Cellular and mitochondrial mechanisms of atrial fibrillation. Basic Research in Cardiology, 115 (6):72.
Monaco C, Nanchahal J, Taylor P, Feldmann M. 2015. Anti-TNF therapy:past, present and future. International Immunology, 27(1):55-62.
Mu D, Wu X, Feijó A, Wu W, Wen Z, Cheng J, Xia L, Yang Q, Shan W, Ge D. 2022. Transcriptome analysis of pika heart tissue re veals mechanisms underlying the adaptation of a keystone species on the roof of the world. Frontiers Genetics, 13:1020789.
Naeije R. 2018. Cœur droit et altitude[Altitude and the right heart]. Revuedes Maladies Respiratoires, 35 (4):441-451.
Niu J, Huang Y, Liu X, Zhang Z, Tang J, Wang B, Lu Y, Cai J, Jian J. 2020. Single-cell RNA-seq reveals different subsets of nonspecific cytotoxic cells in teleost. Genomics, 112 (6):5170-5179.
Ogata H, Goto S, Sato K, Fujibuchi W, Bono H, Kanehisa M. 1999.KEGG:Kyoto encyclopedia of genes and genomes. Nucleic Ac ids Research, 27 (1):29-34.
Ovchinnikova E, Hoes M, Ustyantsev K, Bomer N, de Jong T V, van der Mei H, Berezikov E, vander Meer P. 2018. Modeling human cardiac hypertrophy in stem cell-derived cardiomyocytes. Stem Cell Reports, 10 (3):794-807.
Podyma B, Parekh K, Güler A D, Deppmann C D. 2021. Metabolic homeostasis via BDNF and its receptors. Trends in Endocrinology and Metabolism, 32 (7):488-499.
Rigoulet M, Bouchez C L, Paumard P, Ransac S, Cuvellier S, DuvezinCaubet S, Mazat J P, Devin A. 2020. Cell energy metabolism:An update. Biochim Biophys Acta Bioenerg, 1861 (11):148276.
Sheikh M H, Solito E. 2018. Annexin A1:Uncovering the many talents of an old protein. International Journal of Molecular Sciences, 19 (4):1045.
Sogkas G, Klose K, Baerlecken N, Schweikhard E, Matthias T, Kniesch K, Schmidt R E, Witte T. 2020. CD74 is a T cell antigen in spondyloarthritis. Clinical and Experimental Rheumatology, 38 (2):195-202.
Tseng W W, Wei A C. 2022. Kinetic mathematical modeling of oxidative phosphorylation in cardiomyocyte mitochondria. Cells, 11(24):4020.
Wang D P, Li H G, Li Y J, Guo S C, Yang J, Qi D L, Jin C, Zhao X Q. 2006. Hypoxia-inducible factor 1alpha cDNA cloning and its mRNA and protein tissue specific expression in domestic yak(Bosgrunniens) from Qinghai-Tibetan Plateau. Biochemical and Biophysical Research Communications, 348 (1):310-319.
Wang G, Zhao X, Zhong J, Cao M, He Q, Liu Z, Lin Y, Xu Y, Zheng Y. 2013. Cloning and polymorphisms of yak lactate dehydrogenase B gene. International Journal of Molecular Sciences, 14 (6):11994-12003.
Wang J Y, Zhu S G, Xu C F. 2007. Biochemistry. 3rd ed. Beijing:Higher Education Press, 341-343. (in Chinese)
Wang K, Yang Y, Wang L, Ma T, Shang H, Ding L, Han J, Qiu Q. 2016. Different gene expressions between cattle and yak provide insights into high-altitude adaptation. Animal Genetics, 47 (1):28-35.
Wei Q, Yu H X. 2008. Comparison of histological structure of pulmonary alveoli between 180 days old yak and plain cattle. Journal of Qinghai University, 98 (4):36-39. (in Chinese)
Xin J W, Chai Z X, Zhang C F, Zhang Q, Zhu Y, Cao H W, YangJi C, Chen X Y, Jiang H, Zhong J C, Ji Q M. 2022. Differences in proteomic profiles between yak and three cattle strains provide in sights into molecular mechanisms underlying high-altitude adapta tion. Journal of Animal Physiology and Animal Nutrition, 106(3):485-493.
Xiong X, Fu M, Lan D, Li J, Zi X, Zhong J. 2015. Yak response to high-altitude hypoxic stress by altering mRNA expression and DNA methylation of hypoxia-inducible factors. Animal Biotechnology, 26 (3):222-229.
Yang M S, Zhou N, Wang Z G, Hao H F. 2016. Research progress on the role of transcription factor HIF-1α and its signal pathway in the pathogenesis. Biotechnology Bulletin, 32 (8):8-13. (in Chinese)
Yoon S M, Haritunians T, Chhina S, Liu Z, Yang S, Landers C, Li D, Ye B D, Shih D, Vasiliauskas E A, Ippoliti A, Rabizadeh S, Targan S R, Melmed G Y, McGovern D P B. 2017. Colonic phenotypes are associated with poorer response to anti-TNF therapies in pa tients with IBD. Inflammatory Bowel Diseases, 23 (8):1382-1393.
Zhang Q W, Yu H X, Li L, Jing H X, Wei Q, Liang L, Wang Z Q, Zhang J J. 2018. Histological adaptation of Datongyak cardiac muscle to hypoxia from different altitude areas. Chinese Journal of Veterinary Science, 38 (1):230-234. (in Chinese)
Zhu S H, Zhang G S, Yuan Q. 2013. Comparison of physiological characteristics of two ecotypes of yak in Qinghai Province for en vironmental adaptation. Heilongjiang Animal Science and Veterinary Medicine, (1):152-153. (in Chinese)
王镜岩, 朱圣庚, 徐长法. 2007. 生物化学. 第三版. 北京:高等教育出版社, 341-343.
朱世海, 张国顺, 袁琦. 2013. 青海两种生态类型牦牛对环境适应的生理特性比较. 黑龙江畜牧兽医, (1):152-153.
李筱, 魏莲, 汪洋, 许利娜, 魏琳娜, 魏登邦. 2015. 高原鼠兔心脏中Ldh-c 基因的表达及其对无氧糖酵解水平的影响. 生理学报, 67 (3):312-318.
李潇, 高彦华, 彭忠利, 郭春华, 王琦, 梁红, 王丹. 2022. 海拔高度变化对牦牛血细胞生理指标及心、肺组织形态特征的影响. 家畜生态学报, 43 (4):47-54.
杨梦思, 周娜, 王志钢, 郝慧芳. 2016. 转录因子 HIF-1α 及其信号通路在疾病发生中的作用研究进展. 生物技术通报, 32 (8):8-13.
张勤文, 俞红贤, 李莉, 荆海霞, 魏青, 梁林, 王志强, 张晶晶. 2018.不同海拔梯度大通牦牛心肌组织低氧适应的组织学特点. 中国兽医学报, 38 (1):230-234.
魏青, 俞红贤. 2008. 180 日龄高原牦牛和平原黄牛肺泡组织结构的比较研究. 青海大学学报 (自然科学版), 98 (4):36-39.

[1]	王强龙, 张同享, 高泽川, 王立斌, 余四九, 潘阳阳. 牦牛DBI基因的分子特征及其在乳腺组织中的表达与定位[J]. 兽类学报, 2022, 42(3): 295-303.
[2]	傅芳, 王利, 字向东. 麦洼牦牛不同生长阶段肝脏差异表达基因分析[J]. 兽类学报, 2022, 42(1): 85-94.
[3]	李娟, 王利, 罗晓林, 官久强, 张翔飞. 牦牛水通道蛋白AQP1和AQP3基因克隆及在不同组织中表达和定位[J]. 兽类学报, 2021, 41(1): 51-58.
[4]	傅芳官久强曲秀龙王利罗晓林安添午张翔飞. 不同海拔饲育的麦洼牦牛肺脏组织转录组学分析[J]. 兽类学报, 2020, 40(5): 475-484.
[5]	陈健楠张全伟王琪王媛媛马友记张勇赵兴绪. MTNR1a和MTNR1b在牦牛不同组织中的表达及其在睾丸中的定位[J]. 兽类学报, 2020, 40(1): 64-70.
[6]	王靖雷潘阳阳马睿王萌崔燕余四九. 牦牛Atg5基因克隆及其在输卵管、卵巢和子宫中的表达定位[J]. 兽类学报, 2019, 39(5): 546-555.
[7]	刘传发张良志付海波李文靖张贺李吉叶皮立张堰铭. 野牦牛和家牦牛粪便菌群与短链脂肪酸关系的研究[J]. 兽类学报, 2019, 39(1): 1-7.
[8]	米玛旺堆. 植被覆盖度和土地利用对高原鼠兔觅食与警戒行为的影响[J]. 兽类学报, 2019, 39(1): 69-76.
[9]	王琪张全伟张勇马友记赵兴绪. 天祝白牦牛OXGR1基因多态性与体尺性状的关联性分析[J]. 兽类学报, 2017, 37(4): 379-388.
[10]	汪迁昊刘健锋丁艳平王建林邵宝平. 牦牛Endothelin-1基因CDS序列克隆及其真核表达质粒构建[J]. 兽类学报, 2017, 37(3): 293-299.
[11]	田奎丽潘阳阳李月琴崔燕余四九. 牦牛SND1基因特征序列分析及其蛋白在乳腺的表达[J]. 兽类学报, 2017, 37(2): 207-214.
[12]	孙娟崔燕段德勇何俊峰张倩贺延玉黄玉风. 初生与成年牦牛心肌的组织学结构比较[J]. 兽类学报, 2016, 36(1): 104-111.
[13]	周芸芸张于光卢慧刘芳李迪强冯金朝. 西藏金丝野牦牛的遗传分类地位初步分析[J]. 兽类学报, 2015, 35(1): 48-54.
[14]	景慕娴崔燕杨琨温泽星刘鹏刚申青峰余四九. 急性冷应激对牦牛乳腺上皮细胞 HSP70 mRNA 表达的影响[J]. 兽类学报, 2015, 35(1): 95-101.
[15]	宋乔乔钟金城张成福信金伟姬秋梅柴志欣. 西藏牦牛线粒体DNA 的遗传多样性及系统进化分析[J]. 兽类学报, 2014, 34(4): 356-.