Analysis of Hypoxia-adaptation Related circRNA in Heart Tissues of Tibetan Cattle (Bos taurus) and Sanjiang Cattle (Bos taurus)
BAI Jia-Ling1, WANG Hui1, ZHONG Jin-Cheng1,*, CHEN Zhi-Hua1, CHAI Zhi-Xin1, WANG Ji-Kun1, WANG Jia-Bo1, WU Zhi-Juan1, XIN Jin-Wei2
1 Key Laboratory of Qinghai-Tibetan Plateau Animal Genetic Resource Reservation and Utilization of the Sichuan Province/Key Laboratory of Qinghai-Tibetan Plateau Animal Genetic Resource Reservation and Utilization, Ministry of Education, Southwest Minzu University, Chengdu 610041, China; 2 State Key Laboratory of Hulless Barley and Yak Germplasm Resources and Genetic Improvement, Lhasa 850000, China
Abstract：Circular RNA (circle RNA, circRNA) is a covalently enclosed endogenous biomolecule, which plays an important role in the hypoxia adaptation of plateau animals. High-throughput sequencing technology has been widely used in mammalian and human disease research due to its advantages of low cost, high throughput, and fast speed. However, there are few reports on the circRNAs related hypoxia adaptation in Tibetan cattle (Bos taurus). The heart tissue is sensitive to the high-altitude hypoxia environment. In this study, the circRNAs in the heart tissues of Tibetan cattle and Sanjiang cattle were sequenced using the Illumina HiSeqTM 4000 platform. Then analyzed the differential expression profiling of circRNAs and predicted the targeting relationship between circRNA-miRNA-mRNA, constructed a circRNA-miRNA-mRNA visualization regulatory network. The results showed that a total of 283 differentially expressed circRNAs were filtered between the two species, which included 278 up-regulated and 5 down-regulated. Gene Ontology functional enrichment analysis showed that the source genes of differentially expressed circRNA were annotated to 45 functional terms, and the KEGG pathway analysis showed that they mainly significantly enriched in ubiquitin mediated proteolysis, cGMP-PKG signaling pathway and ect.. In addition, novel_circ_018959 was found had a targeting relationship with miR-142-x, miR-144-y, miR-1908-x, miR- 2302-y, and miR-1273-y, and the targets of these miRNAs including ECE1, NFATC1, CAST, and COX7C were reported as hypoxia adaptation candidate genes, indicating that novel_circ_018959 was likely to be involved in the hypoxia adaptive regulation in Tibetan cattle, Its function and regulation mechanism need to be further studied. The results provide the data support for further revealing the molecular mechanism of circRNA in Tibetan cattle hypoxia adaptation process.
[1] 邓磊, 柴志欣, 王嘉博, 等. 2020. 西藏牛和三江牛大脑组织中低氧适应相关circRNAs的比较分析[J]. 畜牧兽医学报, 51(03): 490-504. (Deng L, Chai Z X, Wang J B, et al.2020. Comparative analysis of hypoxia-adapted circRNAs in brain tissues of Zang cattle and Sanjiang cattle[J]. Acta Veterinaria et Zootechnica Sinica, 51(03): 490-504.) [2] 国家畜禽遗传资源委员会. 2011.中国畜禽遗传资源志: 牛志[M]. 中国农业出版社, 北京. pp. 192-194. (National Animal and Poultry Genetic Resources Committee. 2011. Animal Genetic Resources in China Bovines[M]. China Agriculture Press, Beijing. pp. 192-194.) [3] 侯孟典, 王会, 钟金城, 等. 2020. 基于转录组测序筛选黄牛低氧适应性相关差异基因[J]. 华北农学报, 35(02): 217-227. (Hou M D, Wang H, Zhong J C, et al.2020. Screening differential genes related to hypoxia adaptation in cattle by using transcriptome analysis[J]. Acta Agriculturae Boreali-Sinica, 35(02): 217-227.) [4] 孙福勇, 刘君. 2009. 三江黄牛的生态分布及其品种特点[J]. 草业与畜牧, 2009(09): 51-52. (Sun F Y, Liu J.2009. The ecological distribution and breed characteristics of Sanjiang yellow cattle[J]. Journal of Grassland and Forage Science, 2009(09): 51-52.) [5] 王延东. 2016. 西藏藏猪耐低氧基因ECE1的分子克隆及其表达分析[D]. 硕士学位论文, 山东农业大学, 导师: 曾勇庆. pp. 47-48. (Wang Y D.2016. Molecular cloning and expression analysis of the Xizang Tibetan pig anti-hypoxia gene ECE1[D]. Thesis for M.S., Shandong Agricultural University, Supervisor: Zeng Y Q. pp. 47-48.) [6] 袁晓龙, 潘金春, 蒋瑶, 等. 2020. miR-144、miR-21-3p、miR-142-5p及miR-27b-3p对心肌梗塞后血管再生的作用[J]. 中国实验动物学报, 28(03): 297-306. (Yuan X L, Pan J C, Jiang Y, et al.2020. Effects of miR-144, miR-21-3p, miR-142-5p and miR-27b-3p on angiogenesis after myocardial infarction[J]. Acta Laboratorium Animalis Scientia Sinica, 28(03): 297-306.) [7] Abdelmohsen K, Panda A C, De S, et al.2015. Circular RNAs in monkey muscle: Age-dependent changes[J]. Aging (Albany NY). 7(11): 903-910. [8] Adelman K, Egan E.2017. Non-coding RNA: More uses for genomic junk[J]. Nature. 543(7644): 183-185. [9] Balasubramaniam B, VenkataKrishna L M, Vinitha T, et al.2020. Salmonella enterica Serovar Typhi exposure elicits deliberate physiological alterations and triggers the involvement of ubiquitin mediated proteolysis pathway in Caenorhabditis elegans[J]. International Journal of Biological Macromolecules. 149: 215-233. [10] Cai L D, Qi B Z, Wu X Y, et al.2019. Circular RNA Ttc3 regulates cardiac function after myocardial infarction by sponging miR-15b[J]. Journal of Molecular and Cellular Cardiology. 130: 10-22. [11] Cao Y, You S, Yao Y, et al.2018. Expression profiles of circular RNAs in sheep skeletal muscle[J]. Asian-Australasian Journal of Animal Sciences. 31(10): 1550-1557. [12] Chen Y, Li T, Gao Q, et al.2018. MiR-1908 improves cardiac fibrosis after myocardial infarction by targeting TGF-β1[J]. European Review for Medical and Pharmacological Sciences. 22(7): 2061-2069. [13] Geng H H, Li R, Su Y M, et al.2016. The circular RNA Cdr1as promotes myocardial infarction by mediating the regulation of miR-7a on its target genes expression[J]. PLOS ONE. 11(3): e0151753. [14] Guo M W, Qiu J, Shen F, et al.2020. Comprehensive analysis of circular RNA profiles in skeletal muscles of aging mice and after aerobic exercise intervention[J]. Aging (Albany NY). 12(6): 5071-5090. [15] Hansen T B, Jensen T I, Clausen B H, et al.2013. Natural RNA circles function as efficient microRNA sponges[J]. Nature. 495(7441): 384-388. [16] Huang A Q, Zheng H X, Wu Z Y, et al.2020. Circular RNA-protein interactions: Functions, mechanisms, and identification[J]. Theranostics. 10(8): 3503-3517. [17] Huang M J, Shen Y F, Mao H G, et al.2018. Circular RNA expression profiles in the porcine liver of two distinct phenotype pig breeds[J]. Asian-Australasian Journal of Animal Sciences. 31(6): 812-819. [18] Iguchi K, Saotome M, Yamashita K, et al.2019. Pinacidil, a KATP channel opener, stimulates cardiac Na+/Ca2+ exchanger function through the NO/cGMP/PKG signaling pathway in guinea pig cardiac ventricular myocytes[J]. Naunyn Schmiedebergs Arch Pharmacol. 392(8): 949-959. [19] Jiang R, Li H, Yang J M, et al.2020. CircRNA profiling reveals an abundant circFUT10 that promotes adipocyte proliferation and inhibits adipocyte differentiation via sponging let-7[J]. Molecular Therapy-Nucleic Acids. 20: 491-501. [20] Jiang S Y, Guo C F, Zhang W, et al.2019. The integrative regulatory network of circRNA, microRNA, and mRNA in atrial fibrillation[J]. Frontiers in Genetics. 10: 526. [21] Li M Z, Tian S L, Jin L, et al.2013. Genomic analyses identify distinct patterns of selection in domesticated pigs and Tibetan wild boars[J]. Nature Genetics. 45(12): 1431-1438. [22] Li X L, Wang J, Bai Y, et al.2019. Effect of hypoxia on the muscle fiber switching signal pathways CnA/NFATc1 and myostatin in mouse myocytes[J]. Acta Histochem. 121(5): 539-545. [23] Li Z Y, Huang C, Bao C, et al.2017. Correction:Corrigendum: Exon-intron circular RNAs regulate transcription in the nucleus[J]. Nature Structural & Molecular Biology. 24(2): 194. [24] Liang G M, Yang Y L, Niu G L, et al.2017. Genome-wide profiling of Sus scrofa circular RNAs across nine organs and three developmental stages[J]. DNA Research. 24(5): 523-535. [25] Lin F, Chen H W, Zhao G A, et al.2020. Advances in research on the circRNA-miRNA-mRNA network in coronary heart disease treated with traditional Chinese medicine[J]. Evidence-Based Complementary and Alternative Medicine. 2020(3): 1-10. [26] Persoon S, Paulus M, Hirt S, et al.2018. Cardiac unloading by LVAD support differentially influences components of the cGMP-PKG signaling pathway in ischemic and dilated cardiomyopathy[J]. Heart Vessels. 33(8): 948-957. [27] Shen M M, Li T T, Zhang G X, et al.2019. Dynamic expression and functional analysis of circRNA in granulosa cells during follicular development in chicken[J]. BMC Genomics. 20(1): 96. [28] Wang J, Ren Q L, Hua L S, et al.2019. Comprehensive analysis of differentially expressed mRNA, lncRNA and circRNA and their ceRNA networks in the longissimus dorsi muscle of two different pig breeds[J]. International Journal of Molecular Sciences. 20(5): 1107. [29] Wang K, Yang Y, Wang L, et al.2016. Different gene expressions between cattle and yak provide insights into high-altitude adaptation[J]. Animal Genetics. 47(1): 28-35. [30] Wang K, Long B, Liu F, et al.2016. A circular RNA protects the heart from pathological hypertrophy and heart failure by targeting miR-223[J]. European Heart Journal. 37(33): 2602-2611. [31] Wang Y D, Chen S R, Zhang T, et al.2020. Bioinformatic analysis of peripheral blood miRNA of breast cancer patients in relation with anthracycline cardiotoxicity[J]. BMC Cardiovascular Disorders. 20(1): 43. [32] Wei X F, Li H, Yang J M, et al.2017. Circular RNA profiling reveals an abundant circLMO7 that regulates myoblasts differentiation and survival by sponging miR-378a-3p[J]. Cell Death & Disease. 8(10): e3153. [33] Xu T Y, Wu J, Han P, et al.2017. Circular RNA expression profiles and features in human tissues: a study using RNA-seq data[J]. BMC Genomics. 18(6): 131-142. [34] Yan X M, Zhang Z, Meng Y, et al.2020. Genome-wide identification and analysis of circular RNAs differentially expressed in the longissimus dorsi between Kazakh cattle and Xinjiang brown cattle[J]. PeerJ. 8(3): e8646. [35] Yang M, Wan Q, Hu X, et al.2019. Coexpression modules constructed by weighted gene co-expression network analysis indicate ubiquitin-mediated proteolysis as a potential biomarker of uveal melanoma[J]. Experimental and Therapeutic Medicine. 17(1): 237-243. [36] Yang Y B, Pandurangan M, Hwang I.2012. Changes in proteolytic enzymes mRNAs and proteins relevant for meat quality during myogenesis and hypoxia of primary bovine satellite cells[J]. In Vitro Cellular & Developmental Biology-Animal. 48(6): 359-368. [37] Yuan X L, Pan J C, Wen L J, et al.2019. MiR-144-3p enhances cardiac fibrosis after myocardial infarction by targeting PTEN[J]. Frontiers in Cell and Developmental Biology. 7: 249. [38] Zhang M, Wang Z, Cheng Q S, et al.2020. Circular RNA (circRNA) CDYL induces myocardial regeneration by ceRNA after myocardial infarction[J]. Medical Science Monitor. 26: e923188. [39] Zhang Q, Sun W X, Han J, et al.2020a. The circular RNA hsa_circ_0007623 acts as a sponge of microRNA-297 and promotes cardiac repair[J]. Biochemical and Biophysical Research Communications. 523(4): 993-1000. [40] Zhang Y W, Gou W Y, Zhang Y, et al.2019. Insights into hypoxic adaptation in Tibetan chicken embryos from comparative proteomics[J]. Comparative Biochemistry and Physiology D-Genomics & Proteomics. 31: 100602. [41] Zhang Y F, Guo X, Pei J, et al.2020b. CircRNA expression profile during yak adipocyte differentiation and screen potential circRNAs for adipocyte differentiation[J]. Genes. 11(4): 414. [42] Zhou B, Yu J W.2017. A novel identified circular RNA, circRNA_010567, promotes myocardial fibrosis via suppressing miR-141 by targeting TGF-β1[J]. Biochemical and Biophysical Research Communications. 487(4): 769-775.
