Scientia Agricultura Sinica ›› 2020, Vol. 53 ›› Issue (3): 642-657.doi: 10.3864/j.issn.0578-1752.2020.03.015

• ANIMAL SCIENCE·VETERINARY SCIENCE·RESOURCE INSECT • Previous Articles     Next Articles

Analysis and Identification of circRNAs of Skeletal Muscle at Different Stages of Sheep Embryos Based on Whole Transcriptome Sequencing

SHI TianPei,WANG XinYue,HOU HaoBin,ZHAO ZhiDa,SHANG MingYu,ZHANG Li()   

  1. Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing 100193
  • Received:2019-03-01 Accepted:2019-05-30 Online:2020-02-01 Published:2020-02-13
  • Contact: Li ZHANG E-mail:zhangli07@caas.cn

Abstract:

【Objective】The meat production of livestock, which is closely related to the development of skeletal muscle, is an important economic trait to measure the quality of livestock. For mammals, the skeletal muscle development depends on the growth and differentiation of embryonic myocyte, which has a significant impact on the subsequent growing potential. In this study, the developmental mode of skeletal muscle, the important transformation nodes, the formation of muscle fibers and the molecular regulation mechanism of transformation were mainly explored. 【Method】 Based on the previous research, the important nodes D85, D105 and D135 related to the myotube development were used in the experiment, and the longissimus dorsi muscles were sequenced by whole transcriptome sequencing. The differentially expressed (DE) circRNAs were screened by bioinformatics analysis and verified by quantitative real-time PCR (qRT-PCR). 【Result】 1 126 DE circRNAs were obtained by conditional screening (|log2| ≥1 and P≤0.05). The 3 groups were compared and many specific expressions of circRNA were found at each stage, but in the D85 vs D135 group, the amount was the most. 374 DE circRNAs were obtained, which contained 201 up-regulated and 173 down-regulated, and 44.7% of the DE genes were differentially expressed with a difference of more than 4 times. These DE circRNAs were subjected to run GO and KEGG functional analysis and targeted prediction, and they were enriched into some pathways, such as energy metabolism and signal transduction, which involved in muscle differentiation and muscle fiber development, including MAPK, PI3K-Akt, Ras, regulation of actin cytoskeleton and other signal transduction pathways. According to the results, it was confirmed that the DE circRNAs enriched during D85 to D105 were mostly associated with cell proliferation and survival, regulation of myocyte development and cell cycle, while D105 to D135 were mainly related to energy conversion, material transport, RNA transport, and DNA repair. By drawing co-expression visualization network with the targeted prediction results used by Cytoscape, the core regulatory transcripts, such as circRNA8239, circRNA19073, circRNA2765 and circRNA1616, were identified. In the D105 period, a key factor circRNA7527 that regulated the conversion of fast and slow muscle types was found, which targets the bta-miR-135a, bta-miR-615, and chi-miR-133a-5p to regulate the MEF2C gene. According to the differential expression and functional prediction in three comparison groups, 4 circRNAs related to muscle development and 4 target miRNA were selected for qRT-PCR, and the results showed that the gene expression trend was consistent with the sequencing data. 【Conclusion】 It was verified that the stabilization of the number of muscle fibers occurred between sheep embryos at D85 and D105, and muscle fiber hypertrophy happened during the D105 to D135 period, which lead to the conclusion that D105 was probably a key time point. In this study, we firstly constructed a circRNA map in sheep embryonic skeletal muscle development based on the whole transcriptome sequencing. The transcriptome differences at key stages were revealed, and multiple circRNAs and miRNAs targeting MEF2C that involved in the MAPK signaling pathway were found, which provided reference for livestock myofiber development research and other research on non-coding RNA.

Key words: sheep, embryo, whole transcriptome, skeletal muscle, growth and development, circRNA

Table 1

The list of gene primers of qRT-PCR"

RNA 类型
RNA Type
RNA 名称
RNA ID
上游引物序列
Upstream primer sequence
下游引物序列
Downstream primer sequence
circRNA内参基因
CircRNA reference gene
Actin-β CGTGCGTGACATCAAAGAGAA AACCGCTCGTTGCCAATAGT
circRNA circRNA17939 CTGGATACTGTGTTACTTAGAC CATTCTCAGAGGATTAACGATT
circRNA circRNA1037 GGATGACTTGCCCAAGGAGAA CTTGATGGCGACTCGGACTT
circRNA circRNA1042 GGAGATTGGTTTGGAGAGAAAGC TCCTGTCCTCATCATTATCGTCTG
circRNA circRNA5212 TGCCCTCCTACAGACTTGA CGCTGGGTTGGAGATGTT
miRNA内参基因
miRNA reference gene
U6 CTCGCTTCGGCAGCACA AACGCTTCACGAATTTGCGT
miRNA oar-miR-150 ACGCGTCTCCCAACCCTT GTCGTATCCAGTGCAGGGTCCGAGGTATTC
GCACTGGATACGACCACTGGG
miRNA oar-miR-410-5p CGGCGGAGGTTGTCTGTG GTCGTATCCAGTGCAGGGTCCGAGGTATTCGCACTGGATACGACCGAACT
miRNA oar-miR-133a-5p GGCGGAGCTGGTAAAATGG GTCGTATCCAGTGCAGGGTCCGAGGTATTCGCACTGGATACGACATTTGG
miRNA bta-miR-365-3p
CGCGTAATGCCCCTAAAAAT GTCGTATCCAGTGCAGGGTCCGAGGTATTCGCACTGGATACGACATAAGG
通用下游
Universal downstream primer
CAGTGCAGGGTCCGAGGTAT

Fig. 1

Distribution of circRNAs based on FPRM in various samples"

Fig. 2

Distribution of circRNAs based on FPRM in various samples"

Fig. 3

KEGG pathway analysis of differentially expressed circRNAs A: D85 vs D105; B: D105 vs D135; C: D85 vs D135; Y-axis: Pathway name; X-axis: enrichment factor; Bubble size reflect the number of differential expressions enriched to a certain pathway/all quantities in the background; Color reflect the significance of differential expression in a certain process"

Fig. 4

PI3K-Akt Signal Pathway The white box is a generalized and detailed metabolic map based on the existing data; the greenbox is the unique gene or enzyme of the species with more detailed information; the orange box indicates that the gene differential expressed in here"

Fig. 5

The heat map of differentially expressed Horizontal axis: embryo age; Vertical axis: circRNA name; Blue: down; Red: up"

Fig. 6

miRNA-circRNA interaction network diagram Orange ellipses indicate differentially expressed miRNAs; Green arrows indicate differentially expressed circRNAs; and Edges indicate interactions between the two. The color from light to dark represent the degree of difference in expression, and the deeper the color, the more significant the difference, and vice versa"

Fig. 7

circRNA-miRNA-mRNA interaction network map Orange ellipses indicate differentially expressed mRNA; Blue arrows indicate differentially expressed circRNA; Red diamonds indicate differentially expressed miRNAs; and Edges indicate interactions between the two"

Fig. 8

GO analysis of circRNA-miRNA-mRNA Y-axis: biological process; X-axis: enrichment factor; Bubble size reflect the number of differential expressions enriched to a certain GO_Term/all quantities in the background; Color reflect the significance of differential expression in a certain process"

Fig. 9

The results of qRT-PCR The figure shows the quantitative results of RNA-seq and qRT-PCR; the line graph shows the FPKM value of the sequencing and the histogram shows the relative quantitative results; the significance test: the D105 histogram is marked with the D85 vs D105 comparison result; 2 test result in D135 histogram: left is D85 vs D135 and right is D105 vs D135; *P≤0.05, ** P≤0.01"

[1] BENTZINGER C F, XIN W Y, RUDNICKI M A . Building muscle: Molecular regulation of myogenesis. Cold Spring Harbor Perspectives in Biology, 2012,4(2):a008342.
[2] WEINTRAUB H . The MyoD family and myogenesis: Redundancy, networks, and thresholds. Cell, 1993,75(7):1241-1244.
[3] SABOURIN L A, RUDNICKI M A . The molecular regulation of myogenesis. Clinical Genetics, 2010,57(1):16-25.
[4] BUCKINGHAM M, RIGBY P W . Gene regulatory networks and transcriptional mechanisms that control myogenesis. Developmental Cell, 2014,28(3):225-238.
[5] KAWAKAMI K, SATO S, H, IKEDA K. Six family genes-structure and function as transcription factors and their roles in development. Bioessays, 2000,22(7):616-626.
[6] TESSMAR K, LOOSLI F, WITTBRODT J . A screen for co-factors of Six3. Mechanisms of Development, 2002,117(1):103-113.
[7] ZHU C C, DYER M A, UCHIKAWA M, KONDOH H, LAGUTIN O V, OLIVER G . Six3-mediated auto repression and eye development requires its interaction with members of the Groucho-related family of co-repressors. Development, 2002,129(12):2835-2849.
[8] MCKINNELL I W, JEFF I, FABIEN L G, PUNCH V G J, ADDICKS G C, GREENBLATT J F, F JEFFREY D, RUDNICKI M A. Pax7 activates myogenic genes by recruitment of a histone methyltransferase complex. Nature Cell Biology, 2008,10(1):77-84.
[9] 李雪娇, 刘晨曦, 杨开伦, 刘明军 . 德美羊与中美羊胎儿期骨骼肌组织学结构发育特征差异性研究. 草食家畜, 2017(4):1-6.
LI X J, LIU C X, YNGA K L, LIU M J . Study on differentiation of fetal skeletal muscle development characteristics between German and Chinese Merino Sheep. Grass-feeding Livestock, 2017(4):1-6.(in Chinese)
[10] 张伟 . 绵羊骨骼肌高表达miRNA靶基因的高通量获取及部分候选靶基因生物功能研究[D]. 石河子: 石河子大学, 2015.
ZHANG W . Research on acquiring target genes of miRNA expressed at a high level in sheep skeletal muscle by a high throughput way and function of partial target[D]. Shihezi: Shihezi University, 2015. ( in Chinese)
[11] 张世芳, 魏彩虹, 陆健, 张小宁, 周鑫磊, 张淑珍, 王光凯, 曹家雪, 赵福平, 张莉 , 杜立新深度测序鉴定绵羊microRNA转录组. 中国畜牧兽医, 2013,40(9):19-22.
ZHANG S F, WEI C H, LU J, ZHANG X N, ZHOU X L, ZHANG S Z, WANG G K, CAO J X, ZHAO F P, ZHANG L, DU L X . Identification of microRNAome in texel sheep by deep sequencing. China Animal Husbandry &Veterinary Medicine, 2013, 40(9):19-22. (in Chinese)
[12] 黄万龙, 张秀秀, 李嫒, 苗向阳 . 利用RNA-seq技术筛选大白猪皮下和肌内脂肪组织差异表达基因. 遗传, 2017,39(6):501-511.
HUANG W L, ZHANG X X, LI Y, MIAO X Y . Identification of differentially expressed genes between subcutaneous and intramuscular adipose tissue of Large White pig using RNA-seq. Hereditas(Beijing), 2017,39(6):501-511.(in Chinese)
[13] MARTIN M . Cutadapt removes adapter sequences from high- throughput sequencing reads. EMBnet Journal, 2011, doi: 10.14806/ej.17.1.200.
[14] LANGMEAD B, SALZBERG S L . Fast gapped-read alignment with Bowtie 2. Nature Methods, 2012,9:357-359.
[15] KIM D, PERTEA G, TRAPNELL C, PIMENTEL H, KELLEY R, SALZBERG S L . TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biology, 2013,14(4):R36.
[16] MIHAELA P, PERTEA G M, ANTONESCU C M, TSUNG-CHENG C, MENDELL J T, SALZBERG S L . StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nature Biotechnology, 2015,33(3):290-295.
[17] FRAZEE A C, GEO P, JAFFE A E, BEN L, SALZBERG S L, LEEK J T . Ballgown bridges the gap between transcriptome assembly and expression analysis. Nature Biotechnology, 2015,33(3):243.
[18] ANDERS S, HUBER W . Differential expression analysis for sequence count data. Genome Biology, 2010,11(10):R106.doi: 10.1186/gb-2010-11-10-r106.
[19] KIM D, SALZBERG S L . TopHat-Fusion: an algorithm for discovery of novel fusion transcripts. Genome Biology, 2011,12(8):R72.
[20] ZHANG X O, WANG H B, ZHANG Y, LU X, CHEN L L, YANG L . Complementary sequence-mediated exon circularization. Cell, 2014,159(1):134-147.
[21] ASHWAL-FLUSS R, MEYER M, PAMUDURTI N R, IVANOV A, BARTOK O, HANAN M, EVANTAL N, MEMCZAK S, RAJEWSKY N, KADENEr S . circRNA biogenesis competes with Pre-mRNA splicing. Molecular Cell, 2014,56(1):55-66.
[22] CAMON E, MAGRANE M, BARRELL D, LEE V, DIMMER E, MASLEN J, BINNS D, HARTE N, LOPEZ R, APWEILER R . The gene ontology annotation (GOA) database: sharing knowledge in uniprot with gene ontology. Nucleic Acids Research, 2004,32(Database issue):D262.
[23] HATTORI M, ITOH M, ARAKI M, HIRAKAWA M, KAWASHIMA S, OKUDA S, GOTO S, KATAYAMA T, TOKIMATSU T, YAMANISHI Y, KANEHISA M . KEGG for linking genomes to life and the environment. Nucleic Acids Research, 2007,36(suppl.1):D480-D484.
[24] MUKAI H, TOSHIMORI M, SHIBATA H, TAKANAGA H, KITAGAWA M, MIYAHARA M, SHIMAKAWA M, ONO Y . Interaction of PKN with alpha-actinin. Journal of Biological Chemistry, 1997,272(8):4740.
[25] MINDEN A, LIN A, MCMAHON M, LANGE-CARTER C, DÉRIJARD B, DAVIS R J, JOHNSON G L, KARIN M . Differential activation of ERK and JNK mitogen-activated protein kinases by Raf-1 and MEKK. Science, 1994,266(5191):1719-1723.
[26] LEONARDO S, LAURA P, YVONNE T, LEV K, PIER PAOLO P . A ceRNA hypothesis: the Rosetta Stone of a hidden RNA language? Cell, 2011,146(3):353-358.
[27] 李雪娇, 刘晨曦, 孙亚伟, 杨开伦, 刘明军 . 德国美利奴羊胎儿期骨骼肌组织学结构发育特征研究. 西北农林科技大学学报(自然科学版), 2018,46(5):7-13.
LI X J, LIU C X, SUN Y W, YANG K L, LIU M J . Study on structure development characteristics of German Merino sheep fetal skeletal muscle tissue. Journal of Northwest A&F University(Natural Science Edition), 2018,46(5):7-13.(in Chinese)
[28] WANG D Z, VALDEZ M R, MCANALLY J, RICHARDSON J, OLSON E N . The Mef2c gene is a direct transcriptional target of myogenic bHLH and MEF2 proteins during skeletal muscle development. Development, 2001,128(22):4623.
[29] ZHANG Y Y, VIK T A, RYDER J W, SROUR E F, JACKS T, SHANNON K, CLAPP D W . Nf1 regulates hematopoietic progenitor cell growth and ras signaling in response to multiple cytokines. Journal of Experimental Medicine, 1998,187(11):1893-1902.
[30] HEGEDUS B, DASGUPTA B, SHIN J E, EMNETT R J, HART-MAHON E K, ELGHAZI L, BERNAL-MIZRACHI E, GUTMANN D H . Neurofibromatosis-1 regulates neuronal and glial cell differentiation from neuroglial progenitorsin vivo by both cAMP- and Ras-dependent mechanisms. Cell Stem Cell, 2007,1(4):443-457.
[31] XIAOHUA W, ESTWICK S A, SHI C, MENGGANG Y, WENYU M, NEBESIO T D, YAN L, JIN Y, REUBEN K, DAVID I . Neurofibromin plays a critical role in modulating osteoblast differentiation of mesenchymal stem/progenitor cells. Human Molecular Genetics, 2006,15(19):2837-2845.
[32] NADINE K, SIGMAR S, CHRISTIAN R D, ROBINSON P N, JOHNNY K, CAROLA D, MONIKA O, JIRKO K, STEVENSON D A, THOMAS B . Neurofibromin (Nf1) is required for skeletal muscle development. Human Molecular Genetics, 2011,20(14):2697.
[33] PERRY R L, PARKER M H, RUDNICKI M A . Activated MEK1 binds the nuclear MyoD transcriptional complex to repress transactivation. Molecular Cell, 2001,8(2):291-301.
[34] RAMOCKI M B, JOHNSON S E, WHITE M A, ASHENDEL C L, KONIECZNY S F, TAPAROWSKY E J . Signaling through mitogen-activated protein kinase and Rac/Rho does not duplicate the effects of activated Ras on skeletal myogenesis. Molecular & Cellular Biology, 1997,17(7):3547-3555.
[35] PAGE J L, WANG X L, JOHNSON S E . MEKK1 signaling through p38 leads to transcriptional inactivation of E47 and repression of skeletal myogenesis. Journal of Biological Chemistry, 2004,279(30):30966-30972.
[36] GREDINGER E, GERBER A N, TAMIR Y, TAPSCOTT S J, BENGAL E . Mitogen-activated protein kinase pathway is involved in the differentiation of muscle cells. Journal of Biological Chemistry, 1998,273(17):10436-10444.
[37] BOST F, AOUADI M, CARON L, BINéTRUY B . The role of MAPKs in adipocyte differentiation and obesity. Biochimie, 2005,87(1):51-56.
[38] MYRIAM A, KATHIANE L, MATTHIEU P, YANNICK M B, BERNARD B, FRéDéRIC B . Inhibition of p38MAPK increases adipogenesis from embryonic to adult stages. Diabetes, 2006,55(2):281.
[39] FRéDéRIC B, MYRIAM A, LESLIE C, PATRICK E, NATHALIE B, MATTHIEU P, CHRISTIAN D, PAUL H, GILLES P, JACQUES P . The extracellular signal-regulated kinase isoform ERK1 is specifically required for in vitro and in vivo adipogenesis. Diabetes, 2005,54(2):402-411.
[40] LIU S, HAN W, JIANG S, ZHAO C, WU C . Integrative transcriptomics and proteomics analysis of longissimus dorsi muscles of Canadian double-muscled Large White pigs. Gene, 2016,577(1):14-23.
[41] VOILLET V, SAN CRISTOBAL M, PèRE M-C, BILLON Y, CANARIO L, LIAUBET L, LEFAUCHEUR L . Integrated analysis of proteomic and transcriptomic data highlights late fetal muscle maturation process. Molecular & Cellular Proteomics:MCP, 2018,17(4):672-693.
[42] OUYANG H, WANG Z, CHEN X, YU J, LI Z, NIE Q . Proteomic analysis of chicken skeletal muscle during embryonic development. Frontiers in Physiology, 2017,8:281-317.
[1] YE MeJin, WU Lei, MD NAHIBUZZAMAN Lohani, YIN Li, HU XinRong, LIU YaXi, JIANG YunFeng, CHEN GuoYue, PU ZhiEn, LI Yang, LI Ting, ZOU YaYa, WU JiaYi, MA Jian. Genome-Wide Association Study-Based Identification of Loci Controlling Mature Embryo Size in Chinese Wheat Landraces and Their Genetic Effects Analysis [J]. Scientia Agricultura Sinica, 2026, 59(6): 1157-1171.
[2] GAO YunJie, GUO RuoNan, CHEN ChunWen, DUAN YiFan, ZHANG ZhenJun, NIE ChanChan, LI MengYu, WANG Hui, FENG TingTing, CUI YingYing, DANG GuangHui, LIU SiGuo. Establishment and Preliminary Application of Indirect ELISA Antibody Detection Method for Mycobacterium avium subsp. paratuberculosis in Sheep [J]. Scientia Agricultura Sinica, 2026, 59(3): 655-667.
[3] YAN ChuanBo, GONG Ping, ZHENG WenXin, CHEN XinWen, GUO LeiFeng. Research on Sheep Mounting Behavior Recognition in Complex Scenes Based on an Improved YOLOv11 [J]. Scientia Agricultura Sinica, 2025, 58(18): 3783-3798.
[4] ZHOU Qi, ZHANG Liang, PAN Yu, TU Zhi, WANG Zheng, LIU HangHang, XIAN LingJin, XIA YunHong, PAN HongMei, LONG Xi. Effects of Cycloastragenol on Cellular Senescence of Pig Donor Fibroblast, Cytoskeletal Dynamic, and Early Developmental Stage of Nuclear Transfer Embryo [J]. Scientia Agricultura Sinica, 2025, 58(12): 2453-2474.
[5] JIANG Chao, ZHANG JiuPan, SONG YaPing, SONG XiaoYu, WU Hao, WEI DaWei. Study on the Role of FoxO1 in Regulating the Proliferation, Apoptosis and Differentiation of Bovine Skeletal Muscle Cells [J]. Scientia Agricultura Sinica, 2024, 57(6): 1191-1203.
[6] CHI RunQing, HAN HaiYin, WANG Peng, LI KaiYang, CHU MingXing, LIU YuFang. Estrogen Mediates CircZNF423 as a Sponge for oar-miR-541-3p to Target CALM3 for Regulating Myoblast Proliferation in Sheep [J]. Scientia Agricultura Sinica, 2024, 57(3): 597-612.
[7] ZENG XiangCui, YANG YongNian, LI RuYue, JIANG XueQian, JIANG Xu, XU YanRan, LIU ZhongKuan, LONG RuiCai, KANG JunMei, YANG QingChuan, LI MingNa. Identification of Alfalfa (Medicago sativa) MsCEP Genes and Functional Analysis of Its Regulation in Root Growth and Development [J]. Scientia Agricultura Sinica, 2024, 57(24): 4839-4853.
[8] YANG DongMei, ZHANG JiuPan, SONG YaPing, SONG XiaoYu, JIANG Chao, MA Yun, WEI DaWei. Research on the Regulatory Effects of Bovine Skeletal Muscle Cells on Adipocytes Under Co-Culture Conditions [J]. Scientia Agricultura Sinica, 2024, 57(18): 3704-3718.
[9] ZHANG YingXin, YANG Min, BAI XueBing, CHEN Chang, WU RuiZhi, YANG Ping, CHEN QiuSheng. Morphological Characteristics of Telocytes at Sheep Acupoints and Its Relationship with Surrounding Structures [J]. Scientia Agricultura Sinica, 2023, 56(7): 1417-1428.
[10] ZHAO WeiHong, HAN WenXiong, YANG Bo, MENG WeiKang, CHAI HaiLiang, MA YiMin, ZHANG ZhanSheng, WANG LiFeng, WANG Yan, WANG MingYuan, ZHANG Shan, DING YuLin, WANG JinLing, JIRINTAI Sulijid, WANG FengLong, ZHAO Li, LIU YongHong. Isolation and Genotyping of Mycobacterium avium subsp. paratuberculosis from Sheep in Inner Mongolia [J]. Scientia Agricultura Sinica, 2023, 56(6): 1204-1214.
[11] XI MengXue, SHEN Dan, SHI YiFan, LI ChunMei. Effects of TBHQ on Pyroptosis, Necroptosis and Inflammatory Damage of Chicken Embryonic Lung Tissues Induced by PM2.5 from Chicken Houses [J]. Scientia Agricultura Sinica, 2023, 56(4): 779-787.
[12] GUO YunPeng, TAN HaoYun, GUO Hong, FU MengYun, LI Xin, HU DeBao, ZHANG LinLin, DING XiangBin, GUO YiWen. LNC721 Targeted Regulation MMP9 Affects Bovineskeletal Muscle Satellite Cell Proliferation and Differentiation [J]. Scientia Agricultura Sinica, 2023, 56(24): 4944-4955.
[13] ZHU PeiPei, QIN HaoXiang, ZHANG JianXia. Changes of Endogenous Hormones and Polyamines During Ovule Development of Stenospermocarpic Seedless Grape [J]. Scientia Agricultura Sinica, 2023, 56(23): 4789-4800.
[14] YIN ZiHe, YANG ChengCheng, ZHAO YuHui, ZHAO Li, LÜ XiuRong, YANG ZhenChao, WU YongJun. Sequencing and Functional Analysis of Tomato circRNA During Flowering Stage [J]. Scientia Agricultura Sinica, 2023, 56(21): 4288-4303.
[15] SHEN LongXian, WANG LiTing, HE Ke, DU Xue, YAN FeiFei, CHEN WeiHu, LÜ YaoPing, WANG Han, ZHOU XiaoLong, ZHAO AYong. Effects of Melatonin and Nicotinamide Mononucleotides on Proliferation of Skeletal Muscle Satellite Cells in Goose [J]. Scientia Agricultura Sinica, 2023, 56(2): 391-404.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
No Suggested Reading articles found!