Transcriptome analysis reveals long intergenic non-coding RNAs involved in skeletal muscle growth and development in pig

Identification and functional prediction of lncRNAs associated with intramuscular lipid deposition in Guangling donkeys

Many studies have shown that long non-coding RNAs (lncRNAs) play key regulatory roles in various biological processes. However, the importance and molecular regulatory mechanisms of lncRNAs in donkey intramuscular fat deposition remain to be further investigated. In this study, we used published transcriptomic data from the longissimus dorsi muscle of Guangling donkeys to identify lncRNAs and obtained 196 novel lncRNAs. Compared with the coding genes, the novel lncRNAs and the known lncRNAs exhibited some typical features, such as shorter transcript length and smaller exons. A total of 272 coding genes and 52 lncRNAs were differentially expressed between the longissimus dorsi muscles of the low-fat and high-fat groups. The differentially expressed genes were found to be involved in various biological processes related to lipid metabolism. The potential target genes of differentially expressed lncRNAs were predicted by cis and trans. Functional analysis of lncRNA targets showed that some lncRNAs may act on potential target genes involved in lipid metabolism processes and regulate lipid deposition in the longissimus dorsi muscle. This study provides valuable information for further investigation of the molecular mechanisms of lipid deposition traits in donkeys, which may improve meat traits and facilitate the selection process of donkeys in future breeding.

LncR-133a Suppresses Myoblast Differentiation by Sponging miR-133a-3p to Activate the FGFR1/ERK1/2 Signaling Pathway in Goats

Long noncoding RNAs (lncRNAs) are involved in a variety of biological processes and illnesses. While a considerable number of lncRNAs have been discovered in skeletal muscle to far, their role and underlying processes during myogenesis remain mostly unclear. In this study, we described a new functional lncRNA named lncR-133a. Gene overexpression and interference studies in goat skeletal muscle satellite cells (MuSCs) were used to establish its function. The molecular mechanism by which lncR-133a governs muscle differentiation was elucidated primarily using quantitative real-time PCR (qRT-PCR), Western blotting, dual-luciferase activity assays, RNA immunoprecipitation, biotin-labeled probe, and RNA fluorescence in situ hybridization analyses. LncR-133a was found to be substantially expressed in longissimus thoracis et lumborum muscle, and its expression levels changed during MuSC differentiation in goats. We validated that lncR-133a suppresses MuSC differentiation in vitro. Dual-luciferase reporter screening, Argonaute 2 (AGO2) RNA immunoprecipitation assays, biotin-labeled lncR-133a capture, and fluorescence in situ hybridization showed that lncR-133a interacted with miR-133a-3p. Additionally, miR-133a-3p facilitated MuSC differentiation, but lncR-133a reversed this effect. The luciferase reporter assay and functional analyses established that miR-133a-3p directly targets fibroblast growth factor receptor 1 (FGFR1). Moreover, lncR-133a directly reduced miR-133a-3p’s capacity to suppress FGFR1 expression, and positively regulated the phosphorylation of extracellular signal-regulated kinase 1/2 (ERK1/2). In summary, our results suggested that lncR-133a suppresses goat muscle differentiation by targeting miR-133a-3p and activating FGFR1/ERK1/2 signaling pathway.

Comprehensive Analysis of Long Noncoding RNA Modified by m6A Methylation in Oxidative and Glycolytic Skeletal Muscles

N⁶-methyladenosine (m⁶A) is the most common modification in eukaryotic RNAs. Accumulating evidence shows m⁶A methylation plays vital roles in various biological processes, including muscle and fat differentiation. However, there is a lack of research on lncRNAs’ m⁶A modification in regulating pig muscle-fiber-type conversion. In this study, we identified novel and differentially expressed lncRNAs in oxidative and glycolytic skeletal muscles through RNA-seq, and further reported the m⁶A-methylation patterns of lncRNAs via MeRIP-seq. We found that most lncRNAs have one m⁶A peak, and the m⁶A peaks were preferentially enriched in the last exon of the lncRNAs. Interestingly, we found that lncRNAs’ m⁶A levels were positively correlated with their expression homeostasis and levels. Furthermore, we performed conjoint analysis of MeRIP-seq and RNA-seq data and obtained 305 differentially expressed and differentially m⁶A-modified lncRNAs (dme-lncRNAs). Through QTL enrichment analysis of dme-lncRNAs and PPI analysis for their cis-genes, we finally identified seven key m⁶A-modified lncRNAs that may play a potential role in muscle-fiber-type conversion. Notably, inhibition of one of the key lncRNAs, MSTRG.14200.1, delayed satellite cell differentiation and stimulated fast-to-slow muscle-fiber conversion. Our study comprehensively analyzed m⁶A modifications on lncRNAs in oxidative and glycolytic skeletal muscles and provided new targets for the study of pig muscle-fiber-type conversion.

Knockdown of CDR1as Decreases Differentiation of Goat Skeletal Muscle Satellite Cells via Upregulating miR-27a-3p to Inhibit ANGPT1

Myogenesis is a complex process controlled by several coding and non-coding RNAs (ncRNAs), such as circular RNAs (circRNAs) that are known to function as endogenous microRNAs (miRNAs) sponges. Cerebellar Degeneration-Related protein 1 antisense (CDR1as) is the most spotlighted circRNA that is known as an miR-7 sponge, which has bloomed circRNAs’ research in animal disease and physiology. Here, we screened for miRNAs and mRNA associated with CDR1as and further characterized their regulatory function during muscle differentiation. We found that a total of 43 miRNAs (including miR-107-3p, miR-125b-5p, miR-140-5p, miR-29a-3p, and miR-27a-3p upregulated) and 789 mRNAs (including ANGPT1, E2F2, CCN1, FGFR1, and MEF2C downregulated) were differentially expressed in goat skeletal muscle satellite cells (SMSCs). Further, knockdown of CDR1as and ANGPT1 inhibited SMSCs differentiation. miR-27a-3p was differentially upregulated after the knockdown of CDR1as in SMSCs. Overexpressed miR-27a-3p decreased SMSCs differentiation. Via RNAhybrid and luciferase, miR-27a-3p was identified to regulate ANGPT1. We discovered that miR-27a-3p has an inverse relationship with CDR1as and decreases the expression level of ANGPT1 during SMSCs differentiation. In summary, our study demonstrates that siCDR1as inhibits myoblast differentiation by downregulating ANGPT1 mRNA via miR-27a-3p in SMSCs.

LncRNAs in domesticated animals: from dog to livestock species

Animal genomes are pervasively transcribed into multiple RNA molecules, of which many will not be translated into proteins. One major component of this transcribed non-coding genome is the long non-coding RNAs (lncRNAs), which are defined as transcripts longer than 200 nucleotides with low coding-potential capabilities. Domestic animals constitute a unique resource for studying the genetic and epigenetic basis of phenotypic variations involving protein-coding and non-coding RNAs, such as lncRNAs. This review presents the current knowledge regarding transcriptome-based catalogues of lncRNAs in major domesticated animals (pets and livestock species), covering a broad phylogenetic scale (from dogs to chicken), and in comparison with human and mouse lncRNA catalogues. Furthermore, we describe different methods to extract known or discover novel lncRNAs and explore comparative genomics approaches to strengthen the annotation of lncRNAs. We then detail different strategies contributing to a better understanding of lncRNA functions, from genetic studies such as GWAS to molecular biology experiments and give some case examples in domestic animals. Finally, we discuss the limitations of current lncRNA annotations and suggest research directions to improve them and their functional characterisation.

Analysis of long intergenic non-coding RNAs transcriptomic profiling in skeletal muscle growth during porcine embryonic development

Skeletal muscle growth plays a critical role during porcine muscle development stages. Genome-wide transcriptome analysis reveals that long intergenic non-coding RNAs (lincRNAs) are implicated as crucial regulator involving in epigenetic regulation. However, comprehensive analysis of lincRNAs in embryonic muscle development stages remain still elusive. Here, we investigated the transcriptome profiles of Duroc embryonic muscle tissues from days 33, 65, and 90 of gestation using RNA-seq, and 228 putative lincRNAs were identified. Moreover, these lincRNAs exhibit the characteristics of shorter transcripts length, longer exons, less exon numbers and lower expression level compared with protein-coding transcripts. Expression profile analysis showed that a total of 120 lincRNAs and 2638 mRNAs were differentially expressed. In addition, we also performed quantitative trait locus (QTL) mapping analysis for differentially expressed lincRNAs (DE lincRNAs), 113 of 120 DE lincRNAs were localized on 2200 QTLs, we observed many QTLs involved in growth and meat quality traits. Furthermore, we predicted potential target genes of DE lincRNAs in cis or trans regulation. Gene ontology and pathway analysis reveals that potential targets of DE lincRNAs mostly were enriched in the processes and pathways related to tissue development, MAPK signaling pathway, Wnt signaling pathway, TGF-beta signaling pathway and insulin signaling pathway, which involved in skeletal muscle physiological functions. Based on cluster analysis, co-expression network analysis of DE lincRNAs and their potential target genes indicated that DE lincRNAs highly regulated protein-coding genes associated with skeletal muscle development. In this study, many of the DE lincRNAs may play essential roles in pig muscle growth and muscle mass. Our study provides crucial information for further exploring the molecular mechanisms of lincRNAs during skeletal muscle development.

Analysis of long non-coding RNAs in skeletal muscle of Bama Xiang pigs in response to heat stress

The aim of this study was to identify differentially expressed long non-coding RNAs (lncRNAs) molecules and predict their target genes related to muscle development and lipid metabolism in longissimus dorsi (LD) muscles of Bama Xiang pigs under constant heat stress. Ten male Bama Xiang pigs with an average initial body weight of 14 kg were randomly divided into control group (22°C) and heat stress (35 °C) group. The experiment lasted for 28 days. All the pigs were slaughtered at the end of the experiment, and LD muscles were collected for muscle quality analysis and transcriptome sequencing. Heat stress reduced meat quality of Bama Xiang pigs. lncRNAs in LD were identified systematically by deep RNA sequencing between the two groups. The results showed that 365 lncRNAs from the LD were identified, including 128 intergenic lncRNAs, 82 intronic lncRNAs, and 155 anti-sense lncRNAs. The differences lie in transcript of length, number of exons and wider size distribution, and expression level per KB fragment in three subtypes of lncRNAs. The three types of transposable elements coverage, including Line/L1, SINE/tRNA, and LTR/ERVL-MaLR, are the highest in mRNA and the three subtypes of lncRNAs in pigs. lncRNAs and mRNAs were different in comparison of features. The results predicted the target genes of the significant differentially expressed lncRNAs related to muscle development and lipid metabolism. This is the first study to expand the knowledge about muscle-related lncRNAs biology in Bama Xiang pigs under heat stress and will contribute to the development of alleviating the adverse effects of heat stress on pork quality targeting lncRNAs.

Comparative study of weaning pigs' muscle proteins using two-dimensional electrophoresis

The proteostasis system of animals, including various types of protein modification during the growth stage, leads to an almost incomprehensible number of possible forms of protein, and each can regulate numerous functions. In the presented work, the composition of muscle tissue protein from different portions of piglets was studied to understand the main muscle protein formation. Comparative analysis of weaned piglets' main muscle protein from l. dorsi, biceps femoris, and brachiocephalicus were analyzed using two-dimensional electrophoresis. Changes in the staining intensity of protein fractions inherent in different muscles were revealed. As part of this work, candidate groups of pig muscle proteins have been selected. Eleven protein spots were revealed for the longest muscle of the back, and seven for the biceps; the muscles of the neck are characterized by indicators of low protein fraction volume. Among the proteins found, myosin light chains, phosphoglycerate mutase, troponins, and adenylate kinase is most likely present. The obtained results of protein identification in muscle tissues, obtained during the intensive growth period, will allow a more detailed understanding of protein regulation, function, and interactions in complex biological systems, which will subsequently be significantly important for biomonitoring health and predicting farm animals productivity.

Profiling and Functional Analysis of Long Noncoding RNAs and mRNAs during Porcine Skeletal Muscle Development

Gene transcripts or mRNAs and long noncoding RNAs (lncRNAs) are differentially expressed during porcine skeletal muscle development. However, only a few studies have been conducted on skeletal muscle transcriptome in pigs based on timepoints according to the growth curve for porcine. Here, we investigated gene expression in Qingyu pigs at three different growth stages: the inflection point with the maximum growth rate (MGI), the inflection point of the gradually increasing stage to the rapidly increasing stage (GRI), and the inflection point of the rapidly increasing stage to the slowly increasing stage (RSI). Subsequently, we explored gene expression profiles during muscle development at the MGI, GRI and RSI stages by Ribo-Zero RNA sequencing. Qingyu pigs reached the MGI, GRI and RSI stages at 156.40, 23.82 and 288.97 days of age with 51.73, 3.14 and 107.03 kg body weight, respectively. A total of 14,530 mRNAs and 11,970 lncRNAs were identified at the three stages, and 645, 323 differentially expressed genes (DEGs) and 696, 760 differentially expressed lncRNAs (DELs) were identified in the GRI vs. MGI, and RSI vs. MGI, comparisons. Functional enrichment analysis revealed that genes involved in immune system development and energy metabolism (mainly relate to amino acid, carbohydrate and lipid) were enriched at the GRI and MGI stages, respectively, whereas genes involved in lipid metabolism were enriched at the RSI stage. We further characterized G1430, an abundant lncRNA. The full-length sequence (316 nt) of lncRNA G1430 was determined by rapid amplification of cDNA ends (RACE). Subcellular distribution analysis by quantitative real-time PCR (qRT-PCR) revealed that G1430 is a cytoplasmic lncRNA. Binding site prediction and dual luciferase assay showed that lncRNA G1430 directly binds to microRNA 133a (miR-133a). Our findings provide the basis for further investigation of the regulatory mechanisms and molecular genetics of muscle development in pigs.

Profiling of skeletal muscle tissue for long non-coding RNAs related to muscle metabolism in the QingYu pig at the growth inflection point

OBJECTIVE

Investigation of muscle growth at different developmental stages is an appropriate strategy for studying the mechanisms underlying muscle development and differences in phenotypes. In particular, the muscle development mechanisms and the difference between the fastest and slowest growth rates.

METHODS

In this study, we used a growth curve model to fit the growth inflection point (IP) of QingYu pigs and compared differences in the long non-coding RNA (lncRNA) transcriptome of muscle both at the growth IP and plateau phase (PP).

RESULTS

The growth curve of the QingYu pig had a good fit (R2 = 0.974) relative to a typical S-curve and reached the IP at day 177.96. At the PP, marbling, intramuscular fat, and monounsaturated fatty acids had increased significantly and the percentage of lean muscle and polyunsaturated fatty acids had decreased. A total of 1,199 mRNAs and 62 lncRNAs were differentially expressed at the IP compared with the PP. Additional to gene ontology and Kyoto encyclopedia of genes and genomes pathway analyses, these differentially expressed protein coding genes were principally related to muscle growth and lipid metabolism.

CONCLUSION

Our results suggest that the identified differentially expressed lncRNAs, could play roles in muscle growth, fat deposition and regulation of fatty acid composition at the IP and PP.

Review: Long non-coding RNA in livestock

Less than 2% of mammalian genomes code for proteins, but 'the majority of its bases can be found in primary transcripts' - a phenomenon termed the pervasive transcription, which was first reported in 2007. Even though most of the transcripts do not code for proteins, they play a variety of biological functions, with regulation of gene expression appearing as the most common one. Those transcripts are divided into two groups based on their length: small non-coding RNAs, which are maximally 200 bp long, and long non-coding RNAs (lncRNAs), which are longer than 200 nucleotides. The advances in next-generation sequencing methods provided a new possibility of investigating the full set of RNA molecules in the cell. In this review, we summarized the current state of knowledge on lncRNAs in three major livestock species - Sus scrofa, Bos taurus and Gallus gallus, based on the literature and the content of biological databases. In the NONCODE database, the largest number of identified lncRNA transcripts is available for pigs, but cattle have the largest number of lncRNA genes. Poultry is represented by less than a half of records. Genomic annotation of lncRNAs showed that the majority of them are assigned to introns (pig, poultry) or intergenic (cattle). The comparison with well-annotated human and mouse genomes indicates that such annotation is a result of lack of proper lncRNA annotation data. Since lncRNAs play an important role in genomic studies, their characterization in farm animals' genomes is critical in bridging the gap between genotype and phenotype.

Long Noncoding Ribonucleic Acid MSTRG.59589 Promotes Porcine Skeletal Muscle Satellite Cells Differentiation by Enhancing the Function of PALLD

Skeletal muscle satellite cells are a class of undifferentiated mononuclear myogenic stem cells distributed between the myofibroblast and membrane basement. Since their development determines the development of skeletal muscles, knowledge of their proliferation, differentiation, and fate is vital for understanding skeletal muscle development. Increasing evidence have shown that long noncoding RNA (lncRNA) plays an important role in regulating the development process of satellite cells. Based on the results of our previous studies, we screened lncRNA MSTRG.59589, which is highly expressed in skeletal muscle tissue. In the present study, knockdown of MSTRG.59589 significantly inhibited satellite cell differentiation at various time points, whereas overexpression of MSTRG.59589 demonstrated opposite effects. An MSTRG.59589 knockdown cell model was constructed for transcriptome sequencing, and RNA sequencing analysis screened out a large number of differentially expressed genes. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses of these differentially expressed genes revealed that they are mainly enriched in actin cytoskeleton, muscle contraction, and other pathways related to muscle development. Mechanistic analyses showed that MSTRG.59589 could promote the differentiation process of skeletal muscle satellite cells by positively regulating the expression level of the target gene PALLD. This experiment lays a theoretical foundation for deeper studies on the mechanism of MSTRG.59589 in the differentiation of porcine skeletal muscle satellite cells.

A Novel Long Noncoding RNA, lncR-125b, Promotes the Differentiation of Goat Skeletal Muscle Satellite Cells by Sponging miR-125b

Long noncoding RNAs (lncRNAs) have emerged as essential regulators of skeletal myogenesis, but few myogenesis-associated lncRNAs have been identified and our understanding of their regulatory mechanisms remains limited, particularly in goat. Here, we identified a novel lncRNA, TCONS_00006810 (named lncR-125b), from our previous lncRNA sequencing data on fetal (45, 60, and 105 days of gestation, three biological replicates for each point) and postnatal (3 days after birth, n = 3) goat skeletal muscle, and found that it is highly expressed in skeletal muscle and gradually upregulated during skeletal muscle satellite cell (SMSC) differentiation in goat. Notably, overexpression of lncR-125b accelerated the expression of myogenic differentiation 1 (MyoD 1) and myogenin (MyoG), and the formation of myotubes, and knockdown of lncR-125b showed opposite effects in SMSCs. Results of dual-luciferase assay and quantitative real-time polymerase chain reaction revealed that lncR-125b acts as a molecular sponge for miR-125b and that insulin-like growth factor 2 (IGF2), a critical regulator of skeletal myogenesis, is a direct target gene of miR-125b. Further analyses showed that lncR-125b negatively regulates miR-125b expression and positively regulates IGF2 expression in SMSCs. Mechanistically, lncR-125b promotes SMSC differentiation by functioning as a competing endogenous RNA (ceRNA) for miR-125b to control IGF2 expression. These findings identify lncR-125b as a novel noncoding regulator of muscle cell differentiation and skeletal muscle development in goat.

Functions and Regulatory Mechanisms of lncRNAs in Skeletal Myogenesis, Muscle Disease and Meat Production

Myogenesis is a complex biological process, and understanding the regulatory network of skeletal myogenesis will contribute to the treatment of human muscle related diseases and improvement of agricultural animal meat production. Long noncoding RNAs (lncRNAs) serve as regulators in gene expression networks, and participate in various biological processes. Recent studies have identified functional lncRNAs involved in skeletal muscle development and disease. These lncRNAs regulate the proliferation, differentiation, and fusion of myoblasts through multiple mechanisms, such as chromatin modification, transcription regulation, and microRNA sponge activity. In this review, we presented the latest advances regarding the functions and regulatory activities of lncRNAs involved in muscle development, muscle disease, and meat production. Moreover, challenges and future perspectives related to the identification of functional lncRNAs were also discussed.

Transcriptome Analysis Reveals the Effect of Long Intergenic Noncoding RNAs on Pig Muscle Growth and Fat Deposition

Muscle growth and fat deposition are the two important biological processes in the development of pigs which are closely related to the pig production performance. Long intergenic noncoding RNAs (lincRNAs), with lack of coding potential and the length of at least 200nt, have been extensively studied to play important roles in many biological processes. However, the importance and molecular regulation mechanism of lincRNAs in the process of muscle growth and fat deposition in pigs are still to be further studied comprehensively. In our study, we used the data, including liver, abdominal fat, and longissimus dorsi muscle of 240 days' age of two F2 full-sib female individuals from the white Duroc and Erhualian crossbreed, to identify 581 putative lincRNAs associated with pig muscle growth and fat deposition. The 581 putative lincRNAs shared many common features with other mammalian lincRNAs, such as fewer exons, lower expression levels, and shorter transcript lengths. Cross-tissue comparisons showed that many transcripts were tissue-specific and were involved in the important biological processes in their corresponding tissues. Gene ontology and pathway analysis revealed that many potential target genes (PTGs) of putative lincRNAs were involved in pig muscle growth and fat deposition-related processes, including muscle cell proliferation, lipid metabolism, and fatty acid degradation. In Quantitative Trait Locus (QTLs) analysis, some PTGs were screened from putative lincRNAs, MRPL12 is associated with muscle growth, GCGR and SLC25A10 were associated with fat deposition, and PPP3CA, DPYD, and FGGY were related not only to muscle growth but also to fat deposition. Therefore, it implied that these lincRNAs might participate in the biological processes related to muscle growth or fat deposition through homeostatic regulation of PTGs, but the detailed molecular regulatory mechanisms still needed to be further explored. This study lays the molecular foundation for the in-depth study of the role of lincRNAs in the pig muscle growth and fat deposition and further provides the new molecular markers for understanding the complex biological mechanisms of pig muscle growth and fat deposition.

Transcriptome Analysis Suggests the Roles of Long Intergenic Non-coding RNAs in the Growth Performance of Weaned Piglets

Long intergenic non-coding RNAs (lincRNAs) have been considered to play a key regulatory role in various biological processes. An increasing number of studies have utilized transcriptome analysis to obtain lincRNAs with functions related to cancer, but lincRNAs affecting growth rates in weaned piglets are rarely described. Although lincRNAs have been systematically identified in various mouse tissues and cell lines, studies of lincRNA in pigs remain rare. Therefore, identifying and characterizing novel lincRNAs affecting the growth performance of weaned piglets is of great importance. Here, we reconstructed 101,988 lincRNA transcripts and identified 1,078 lincRNAs in two groups of longissimus dorsi muscle (LDM) and subcutaneous fat (SF) based on published RNA-seq datasets. These lincRNAs exhibit typical characteristics, such as shorter lengths and lower expression relative to protein-encoding genes. Gene ontology analysis revealed that some lincRNAs could be involved in weaned piglet related processes, such as insulin resistance and the AMPK signaling pathway. We also compared the positional relationship between differentially expressed lincRNAs (DELs) and quantitative trait loci (QTL) and found that some of DELs may play an important role in piglet growth and development. Our work details part of the lincRNAs that may affect the growth performance of weaned piglets and promotes future studies of lincRNAs for molecular-assisted development in weaned piglets.

Identification and Functional Prediction of Long Intergenic Non-coding RNAs Related to Subcutaneous Adipose Development in Pigs

An increasing number of studies have shown that long intergenic non-coding RNAs (lincRNAs) are a very important class of non-coding RNAs that plays a vital role in many biological processes. Adipose tissue is an important place for storing energy, but few studies on lincRNAs were related to pig subcutaneous fat development. Here, we used published RNA-seq data from subcutaneous adipose tissue of Italian Large White pigs and identified 252 putative lincRNAs, wherein 34 were unannotated. These lincRNAs had relatively shorter length, lower number of exons, and lower expression level compared with protein-coding transcripts. Gene ontology and pathway analysis indicated that the adjacent genes of lincRNAs were involved in lipid metabolism. In addition, differentially expressed lincRNAs (DELs) between low and high backfat thickness pigs were identified. Through the detection of quantitative trait locus (QTL), DELs were mainly located in QTLs related to adipose development. Based on the expression correlation of DEL genes and their differentially expressed potential target genes, we constructed a co-expression network and a potential pathway of DEL's effect on lipid metabolism. Our study identified and analyzed lincRNAs in subcutaneous adipose tissue, and results suggested that lincRNAs may be involved in the regulation of subcutaneous fat development. Our findings provided new insights into the biological function of porcine lincRNAs.

Identification and Functional Analysis of Long Intergenic Non-coding RNAs Underlying Intramuscular Fat Content in Pigs

Intramuscular fat (IMF) content is an important trait that can affect pork quality. Previous studies have identified many genes that can regulate IMF. Long intergenic non-coding RNAs (lincRNAs) are emerging as key regulators in various biological processes. However, lincRNAs related to IMF in pig are largely unknown, and the mechanisms by which they regulate IMF are yet to be elucidated. Here we reconstructed 105,687 transcripts and identified 1,032 lincRNAs in pig longissimus dorsi muscle (LDM) of four stages with different IMF contents based on published RNA-seq. These lincRNAs show typical characteristics such as shorter length and lower expression compared with protein-coding genes. Combined with methylation data, we found that both the promoter and genebody methylation of lincRNAs can negatively regulate lincRNA expression. We found that lincRNAs exhibit high correlation with their protein-coding neighbors in expression. Co-expression network analysis resulted in eight stage-specific modules, gene ontology and pathway analysis of them suggested that some lincRNAs were involved in IMF-related processes, such as fatty acid metabolism and peroxisome proliferator-activated receptor signaling pathway. Furthermore, we identified hub lincRNAs and found six of them may play important roles in IMF development. This work detailed some lincRNAs which may affect of IMF development in pig, and facilitated future research on these lincRNAs and molecular assisted breeding for pig.