Identification and functional prediction of long non-coding RNAs related to skeletal muscle development in Duroc pigs

다국어 초록 (Multilingual Abstract)

Objective: The growth of pigs involves multiple regulatory mechanisms, and modern molecular breeding techniques can be used to understand the skeletal muscle growth and development to promote the selection process of pigs. This study aims to explore candidate lncRNAs and mRNAs related to skeletal muscle growth and development among Duroc pigs with different average daily gain (ADG).
Methods: A total of 8 pigs were selected and divided into two groups: H group (high-ADG) and L group (low-ADG). And followed by whole transcriptome sequencing to identify differentially expressed (DE) lncRNAs and mRNAs.
Results: In RNA-seq, 703 DE mRNAs (263 up-regulated and 440 down-regulated) and 74 DE lncRNAs (45 up-regulated and 29 down-regulated) were identified. In addition, 1,418 Transcription factors (TFs) were found. Compared with mRNAs, lncRNAs had fewer exons, shorter transcript length and open reading frame length. DE mRNAs and DE lncRNAs can form 417 lncRNA-mRNA pairs (antisense, cis and trans). DE mRNAs and target genes of lncRNAs were enriched in cellular processes, biological regulation, and regulation of biological processes. In addition, quantitative trait locus (QTL) analysis was used to detect the functions of DE mRNAs and lncRNAs, the most of DE mRNAs and target genes of lncRNAs were enriched in QTLs related to growth traits and skeletal muscle development.
In single-nucleotide polymorphism/insertion-deletion (SNP/INDEL) analysis, 1,081,182 SNP and 131,721 INDEL were found, and transition was more than transversion. Over 60% of percentage were skipped exon events among alternative splicing events.
Conclusion: The results showed that different ADG among Duroc pigs with the same diet maybe due to the DE mRNAs and DE lncRNAs related to skeletal muscle growth and development.

다국어 초록 (Multilingual Abstract)

Objective: The growth of pigs involves multiple regulatory mechanisms, and modern molecular breeding techniques can be used to understand the skeletal muscle growth and development to promote the selection process of pigs. This study aims to explore candidate lncRNAs and mRNAs related to skeletal muscle growth and development among Duroc pigs with different average daily gain (ADG).Methods: A total of 8 pigs were selected and divided into two groups: H group (high-ADG) and L group (low-ADG). And followed by whole transcriptome sequencing to identify differentially expressed (DE) lncRNAs and mRNAs.Results: In RNA-seq, 703 DE mRNAs (263 up-regulated and 440 down-regulated) and 74 DE lncRNAs (45 up-regulated and 29 down-regulated) were identified. In addition, 1,418 Transcription factors (TFs) were found. Compared with mRNAs, lncRNAs had fewer exons, shorter transcript length and open reading frame length. DE mRNAs and DE lncRNAs can form 417 lncRNA-mRNA pairs (antisense, cis and trans). DE mRNAs and target genes of lncRNAs were enriched in cellular processes, biological regulation, and regulation of biological processes. In addition, quantitative trait locus (QTL) analysis was used to detect the functions of DE mRNAs and lncRNAs, the most of DE mRNAs and target genes of lncRNAs were enriched in QTLs related to growth traits and skeletal muscle development. In single-nucleotide polymorphism/insertion-deletion (SNP/INDEL) analysis, 1,081,182 SNP and 131,721 INDEL were found, and transition was more than transversion. Over 60% of percentage were skipped exon events among alternative splicing events.Conclusion: The results showed that different ADG among Duroc pigs with the same diet maybe due to the DE mRNAs and DE lncRNAs related to skeletal muscle growth and development.

참고문헌 (Reference)

1 Lorenz R, "Vienna RNA Package 2. 0" 6 : 26-, 2011

2 Hao Y, "Transcriptome analysis reveals that constant heat stress modifies the metabolism and structure of the porcine longissimus dorsi skeletal muscle" 291 : 2101-2115, 2016

3 Sebastian S, "Tissue-specific splicing of a ubiquitously expressed transcription factor is essential for muscle differentiation" 27 : 1247-1259, 2013

4 Zhou M, "The significance of serum S100A9 and TNC levels as biomarkers in colorectal cancer" 10 : 5315-5323, 2019

5 Midwood KS, "Tenascin-C at a glance" 129 : 4321-4327, 2016

6 Yoshida T, "Tenascin-C and integrins in cancer" 9 : 96-104, 2015

7 Kumar A, "SWELL1 regulates skeletal muscle cell size, intracellular signaling, adiposity and glucose metabolism" 9 : e58941-, 2020

8 Tomida T, "Roles of p38 MAPK signaling in the skeletal muscle formation, regeneration, and pathology" 155 : 241-247, 2020

9 Tafer H, "RNAplex : A fast tool for RNA-RNA interaction search" 24 : 2657-2663, 2008

10 Bunch H, "RNA polymerase II promoter-proximal pausing in mammalian long non-coding genes" 108 : 64-77, 2016

11 Aphasizhev R, "RNA and DNA editing: methods and protocols(Chapter 11)" Humana Press 171-184, 2011

12 Shen S, "RMATS : Robust and flexible detection of differential alternative splicing from replicate RNA-Seq data" 111 (111): E5593-E5601, 2014

13 Simionescu-Bankston A, "Noncoding RNAs in the regulation of skeletal muscle biology in health and disease" 94 : 853-866, 2016

14 Wu Y, "MiRPara : A SVM-based software tool for prediction of most probable microRNA coding regions in genome scale sequences" 12 : 107-, 2011

15 Li Y, "Long noncoding RNA SAM promotes myoblast proliferation through stabilizing Sugt1 and facilitating kinetochore assembly" 11 : 2725-, 2020

16 Li Y, "Long non-coding RNAs in the regulation of skeletal myogenesis and muscle diseases" 417 : 58-64, 2018

17 Li Y, "Long non-coding RNAs in the regulation of skeletal myogenesis and muscle diseases" 417 : 58-64, 2018

18 Ma M, "LncRNA-Six1 is a target of miR-1611that functions as a ceRNA to regulate six1 protein expression and fiber type switching in chicken myogenesis" 7 : 243-, 2018

19 Cai B, "LncRNA-Six1 encodes a micropeptide to activate six1 in cis and is involved in cell proliferation and muscle growth" 8 : 230-, 2017

20 Wang L, "LncRNA Dum interacts with Dnmts to regulate Dppa2 expression during myogenic differentiation and muscle regeneration" 25 : 335-350, 2015

21 Jin CF, "Lnc133b, a novel, long non-coding RNA, regulates bovine skeletal muscle satellite cell proliferation and differentiation by mediating miR-133b" 630 : 35-43, 2017

22 Chen LL, "Linking long noncoding RNA localization and function" 41 : 761-772, 2016

23 Zillikens MC, "Large meta-analysis of genome-wide association studies identifies five loci for lean body mass" 8 : 80-, 2017

24 Sarah, Djebali, "Landscape of transcription in human cells" 489 : 101-108, 2012

25 Sun JY, "Knockdown of UACA inhibitsproliferation and invasion and promotes senescence of hepatocellular carcinoma cells" 11 : 4666-4675, 2018

26 Kanehisa M, "KEGG for linking genomes to life and the environment" 36 : D480-D484, 2008

27 Cabili MN, "Integrative annotation of human large intergenic noncoding RNAs reveals global properties and specific subclasses" 25 : 1915-1927, 2011

28 Ramaswami G, "Identifying RNA editing sites using RNA sequencing data alone" 10 : 128-132, 2013

29 Kim D, "HISAT : A fast spliced aligner with low memory requirements" 12 : 357-360, 2015

30 Li D, "Genomic analyses revealed the genetic difference and potential selection genes of growth traits in two duroc lines" 8 : 725367-, 2021

31 Yu J, "Genomewide detection of selection signatures in duroc revealed candidate genes relating to growth and meat quality" 10 : 3765-3773, 2020

32 Shanshan W, "Functions and regulatory mechanisms of lncRNAs in skeletal myogenesis, muscle disease and meat production" 8 : 107-, 2019

33 Zammit PS, "Function of the myogenic regulatory factors Myf5, MyoD, Myogenin and MRF4 in skeletal muscle, satellite cells and regenerative myogenesis" 72 : 19-32, 2017

34 Jin W, "Function analysis of Mef2c promoter in muscle differentiation" 64 : 647-656, 2017

35 Der-Auwera GAV, "From FastQ data to high confidence variant calls : the genome analysis toolkit best practices pipeline" 43 : 11-, 2013

36 Chen S, "Fastp : An ultra-fast all-inone FASTQ preprocessor" 34 : i884-i890, 2018

37 Langmead B, "Fast gapped-read alignment with Bowtie 2" 9 : 357-359, 2012

38 Li R, "Exploring the lncRNAs related to skeletal muscle fiber types and meat quality traits in pigs" 11 : 883-, 2020

39 Lam MT, "Enhancer RNAs and regulated transcriptional programs" 39 : 170-182, 2014

40 Gennarelli M, "Different expression of the myotonin protein kinase gene in discrete areas of human brain" 216 : 489-494, 1995

41 Yotsukura S, "Computational recognition for long noncoding RNA(lncRNA) : Software and databases" 18 : 9-27, 2017

42 Zhou Y, "Comparative whole genome DNA methylation profiling across cattle tissues reveals global and tissue-specific methylation patterns" 18 : 85-, 2020

43 Liu H, "Comparative transcriptomic analysis of skeletal muscle tissue during prenatal stages in Tongcheng and Yorkshire pig using RNA-seq" 18 : 195-209, 2018

44 Luo W, "C-Myc inhibits myoblast differentiation and promotes myoblast proliferation and muscle fibre hypertrophy by regulating the expression of its target genes, miRNAs and lincRNAs" 26 : 426-442, 2019

45 Gu H, "Analysis of differential gene expression of the transgenic pig with overexpression of PGC1alpha in muscle" 46 : 3427-3435, 2019

46 Zhang M, "Alternative splicing of MEF2C pre-mRNA controls its activity in normal myogenesis and promotes tumorigenicity in rhabdomyosarcoma cells" 290 : 310-324, 2015

47 Bahn JH, "Accurate identification of A-to-I RNA editing in human by transcriptome sequencing" 22 : 142-150, 2012

48 Wang K, "ANNOVAR : Functional annotation of genetic variants from high-throughput sequencing data" 38 : e164-, 2010

49 Maas S, "A-to-I RNA editing and human disease" 3 : 1-9, 2006