Identification of circRNA-associated ceRNA networks using longissimus thoracis of pigs of different breeds and growth stages

Background

Long-term artificial selection for growth rate and lean meat rate has eventually led to meat quality deterioration. Muscle fiber type is a key factor that markedly affects meat quality. circRNAs have been reported to participate in diverse biological activities, including myofiber growth and development; thus, we herein compared porcine circRNA transcriptome between oxidative and glycolytic muscle tissues.

Results

Longissimus thoracis muscle tissues were obtained from Lantang and Landrace pigs at birth (LT1D and LW1D, respectively) and 90 postnatal days (LT90D and LW90D, respectively). Hematoxylin and eosin staining and quantitative real-time PCR revealed that all structural traits of the muscle showed large variations between different breeds and growth stages. In total, 329 known miRNAs and 42,081 transcript candidates were identified; 6,962 differentially expressed transcripts were found to play a key role in myogenesis by gene ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses. In addition, 3,352 circRNAs were identified using five predicting algorithms, and 104 circRNA candidates were differentially expressed. Integrated analysis of differentially expressed miRNAs, mRNAs, and circRNAs led to the identification of 777, 855, and 22 convincing ceRNA interactions in LT1D vs. LT90D, LW1D vs. LW90D, and LT90D vs. LW90D, respectively. Finally, we identified a circRNA candidate circKANSL1L, which showed high homology between mice and pigs, and it was found to inhibit the proliferation of C₂C₁₂ cells but promote their differentiation.

Conclusions

We identified genome-wide circRNAs in 0- and 90-day-old Lantang and Landrace pigs by RNA-seq and found that circRNAs were abundant, differentially expressed, and associated with myogenesis. Our results should serve as a reference for future studies on pork quality.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08515-7.

Identification and Characterization of Circular RNAs in Mammary Tissue from Holstein Cows at Early Lactation and Non-Lactation

In this study, circular RNAs (circRNAs) from Holstein cow mammary tissues were identified and compared between early lactation and non-lactation. After analysis, 10,684 circRNAs were identified, ranging from 48 to 99,406 bp, and the average size was 882 bp. The circRNAs were mainly distributed on chromosomes 1 to 11, and 89.89% of the circRNAs belonged to sense-overlapping circRNA. The exons contained with circRNAs ranged from 1 to 47 and were concentrated from 1 to 5. Compared with the non-lactating cows, 87 circRNAs were significantly differentially expressed in the peak lactation cows. There were 68 upregulated circRNAs and 19 downregulated circRNAs. Enrichment analysis of circRNAs showed that GO analysis mainly focused on immune response, triglyceride transport, T cell receptor signaling pathway, etc. Pathway analysis mainly focused on cytokine-cytokine receptor interaction, T helper 17 cell differentiation, fatty acid biosynthesis, the JAK-STAT signaling pathway, etc. Specific primers were designed for two proximal ends of the circRNA junction sites to allow for PCR validation of four randomly selected circRNAs and carry out circRNA-miRNA interaction research. This study revealed the expression profile and characteristics of circRNAs in mammary tissue from Holstein cows at early lactation and non-lactation, thus providing rich information for the study of circRNA functions and mechanisms, as well as potential candidate miRNA genes for studying lactation in Holstein cows.

The Competitive Endogenous RNA (ceRNA) Regulation in Porcine Alveolar Macrophages (3D4/21) Infected by Swine Influenza Virus (H1N1 and H3N2)

H1N1 and H3N2 are the two most common subtypes of swine influenza virus (SIV). They not only endanger the pig industry, but are also a huge risk of zoonotic diseases. However, the molecular mechanism and regulatory network of pigs (hosts) against influenza virus infection are still unclear. In this study, porcine alveolar macrophage cell (3D4/21) models infected by swine influenza virus (H1N1 and H3N2) were constructed. The expression profiles of miRNAs, mRNAs, lncRNAs and circRNAs after H1N1 and H3N2 infected 3D4/21 cells were revealed in this study. Then, two ceRNAs (TCONS_00166432-miR10391-MAN2A1 and novel_circ_0004733-miR10391-MAN2A1) that regulated H1N1 and H3N2 infection in 3D4/21 cells were verified by the methods of bioinformatics analysis, gene overexpression, gene interference, real-time quantitative PCR (qPCR), dual luciferase activity assay and RNA immunoprecipitation (RIP). In addition, the important candidate molecules (miR-10391, TCONS_00166432, and novel_circ_0004733) were identified by qPCR and enzyme linked immunosorbent assay (ELISA). Finally, the regulatory effect and possible molecular mechanism of the target gene MAN2A1 were identified by the methods of gene interference, qPCR, Western blot and ELISA. The results of this study suggested that TCONS_00166432 and novel_circ_0004733 could competitively bind miR-10391 to target the MAN2A1 gene to regulate swine influenza virus infecting 3D4/21 cells. This study reported for the first time the ceRNA networks involved in the regulation of the swine influenza virus infecting 3D4/21 cells, which provided a new insight into the molecular mechanism of 3D4/21 cells against swine influenza virus infection.

Circular RNA, microRNA and Protein Profiles of the Longissimus Dorsi of Germany ZIKA and Sichuan White Rabbits

Due to the dietetic properties and remarkable nutritive value of rabbit meat, its industry is increasing rapidly. However, the association between circular RNAs, microRNAs, and proteins and muscle fiber type, and meat quality of rabbit is still unknown. Here, using deep sequencing and iTRAQ proteomics technologies we first identified 3159 circRNAs, 356 miRNAs, and 755 proteins in the longissimus dorsi tissues from Sichuan white (SCWrabs) and Germany great line ZIKA rabbits (ZIKArabs). Next, we identified 267 circRNAs, 3 miRNAs, and 29 proteins differentially expressed in the muscle tissues of SCWrabs and ZIKArabs. Interaction network analysis revealed some key regulation relationships between noncoding RNAs and proteins that might be associated with the muscle fiber type and meat quality of rabbit. Further, miRNA isoforms and gene variants identified in SCWrabs and ZIKArabs revealed some pathways and biological processes related to the muscle development. This is the first study of noncoding RNA and protein profiles for the two rabbit breeds. It provides a valuable resource for future studies in rabbits and will improve our understanding of the molecular regulation mechanisms in the muscle development of livestock. More importantly, the output of our study will benefit the researchers and producers in the rabbit breeding program.

Circular RNA screening identifies circMYLK4 as a regulator of fast/slow myofibers in porcine skeletal muscles

The type of myofiber is related to the quality of meat. The slow oxidized myofiber helps to increase the tenderness and juiciness of muscle. Numerous studies have shown that circRNA plays a key role in skeletal muscle development. However, the role of circRNA in porcine skeletal myofiber types is unclear. In this study, we performed high-throughput RNA sequencing to study the differential expression of circRNA in the longissimus dorsi and the soleus muscle. A total of 40,757 circRNAs were identified, of which 181 were significantly different. Interestingly, some circRNAs were involved in metabolism pathways, AMPK, FoxO, and PI3K-Akt signaling pathways. Besides, we focused on a novel circRNA-circMYLK4. By injecting circMYLK4-AAV into piglets, we found that circMYLK4 significantly increased the mRNA and protein levels of the slow muscle marker genes. In summary, our study laid an essential foundation for further research of circRNA in myofiber type conversion and higher meat quality.

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.

Integrated Analysis Reveals a lncRNA-miRNA-mRNA Network Associated with Pigeon Skeletal Muscle Development

Growing evidence has demonstrated the emerging role of long non-coding RNA as competitive endogenous RNA (ceRNA) in regulating skeletal muscle development. However, the mechanism of ceRNA regulated by lncRNA in pigeon skeletal muscle development remains unclear. To reveal the function and regulatory mechanisms of lncRNA, we first analyzed the expression profiles of lncRNA, microRNA (miRNA), and mRNA during the development of pigeon skeletal muscle using high-throughput sequencing. We then constructed a lncRNA-miRNA-mRNA ceRNA network based on differentially expressed (DE) lncRNAs, miRNAs, and mRNAs according to the ceRNA hypothesis. Functional enrichment and short time-series expression miner (STEM) analysis were performed to explore the function of the ceRNA network. Hub lncRNA-miRNA-mRNA interactions were identified by connectivity degree and validated using dual-luciferase activity assay. The results showed that a total of 1625 DE lncRNAs, 11,311 DE mRNAs, and 573 DE miRNAs were identified. A ceRNA network containing 9120 lncRNA-miRNA-mRNA interactions was constructed. STEM analysis indicated that the function of the lncRNA-associated ceRNA network might be developmental specific. Functional enrichment analysis identified potential pathways regulating pigeon skeletal muscle development, such as cell cycle and MAPK signaling. Based on the connectivity degree, lncRNAs TCONS_00066712, TCONS_00026594, TCONS_00001557, TCONS_00001553, and TCONS_00003307 were identified as hub genes in the ceRNA network. lncRNA TCONS_00026594 might regulate the FSHD region gene 1 (FRG1)/ SRC proto-oncogene, non-receptor tyrosine kinase (SRC) by sponge adsorption of cli-miR-1a-3p to affect the development of pigeon skeletal muscle. Our findings provide a data basis for in-depth elucidation of the lncRNA-associated ceRNA mechanism underlying pigeon skeletal muscle development.

Spatiotemporal Regulation and Functional Analysis of Circular RNAs in Skeletal Muscle and Subcutaneous Fat during Pig Growth

Recently, thousands of circular RNAs have been reported in different pig breeds. However, researches on the temporal and spatial expression patterns of circRNA over the period of animal growth are limited. Here, we systematically analyzed circRNAs in skeletal muscle and subcutaneous fat in four growth time points (30 days, 90 days, 150 days and 210 days after birth) of a Chinese native pig breed, Ningxiang pigs. A total of 1171 differentially expressed (DE) circRNAs between muscle and fat were identified, including 562 upregulated and 609 downregulated circRNAs. KEGG pathway enrichment analysis of these DE circRNAs revealed that host genes were mainly involved in glycogen metabolism signaling pathways, muscle development signaling pathways such as ErbB pathway and adipocytokine signaling pathways and AMPK signaling pathways and fatty acid biosynthesis. The circRNAs have striking spatiotemporal specificity in the form of dynamic expression at 90 d. Short Time-Series Expression Miner analysis showed multiple model spectra that are significantly enriched with time changes in muscle and fat. Our findings provide new ideas and perspectives about the role of circular RNAs and their targeting relations with mRNA and miRNA in skeletal muscle and fat tissue during pig growth.

Comparison of long non-coding RNAs in adipose and muscle tissues between seven indigenous Chinese and the Yorkshire pig breeds

lncRNAs play crucial roles in fat metabolism in animals. Previously, we have compared the mRNA transcriptome profiles between seven fat-type Chinese pig breeds and one lean-type Western breed (Yorkshire, YY). The associations between differentially expressed (DE) genes and phenotypical traits were investigated. In the present study, to further explore the underlying regulatory mechanisms, lncRNAs were sequenced and compared between YY and Chinese indigenous breeds. The results showed 9114 and 7538 DE lncRNAs between at least one Chinese breed and the YY breed in the adipose and muscle tissue respectively. KEGG enrichment analysis revealed that the target genes of these DE lncRNAs mainly influenced the glucolipid metabolism, which is an important process affecting meat quality. Correlation analyses between the DE lncRNA and DE mRNA genes related to meat quality and growth traits were performed. The results showed that LTCONS_00073280 was associated with intramuscular fat content. Four lncRNAs (LTCONS_00101781, LTCONS_00037879, LTCONS_00088260 and LTCONS-00128343) might mediate backfat thickness. Overall, this study provides candidate lncRNAs that potentially affect meat quality, which might be useful for molecular breeding of pig breeds in future.

Biological mechanisms of growth performance and meat quality in porcine muscle tissue

Growth performance and meat quality are important traits for pig production. The aim of the present study was to investigate the molecular mechanisms underlying growth performance and meat quality, and to identify novel target molecules for predicting the growth performance and meat quality. The differentially expressed genes (DEGs) in Diannan small ears pigs (DSP) and Landrace pigs (LP) were assessed by RNA-sequencing analyzing technology. A total of 339 DEGs were obtained between DSP and LP. 146 DEGs were upregulated in LP compared with DSP and 193 DEGs were upregulated in DSP compared with LP. The DEGs were significantly enriched in 26 GO and 3 KEGG pathways. The protein-protein interaction (PPI) network with 201 nodes and 382 edges was constructed and 5 modules were extracted from the entire network. The identified upregulated expression of genes involved in glycolysis and myogenesis as well as extracellular matrix may be associated with fast body and muscle deposition rates in LP. Increased expression of genes involved in PPAR signaling pathway and fatty acid metabolism as well as oxidative phosphate processes could be related to the intramuscular fat deposition and meat quality in DSP. The present study may provide an improved understanding of the growth performance and meat quality.

The Roles of CircRNAs in Regulating Muscle Development of Livestock Animals

The muscle growth and development of livestock animals is a complex, multistage process, which is regulated by many factors, especially the genes related to muscle development. In recent years, it has been reported frequently that circular RNAs (circRNAs) are involved widely in cell proliferation, cell differentiation, and body development (including muscle development). However, the research on circRNAs in muscle growth and development of livestock animals is still in its infancy. In this paper, we briefly introduce the discovery, classification, biogenesis, biological function, and degradation of circRNAs and focus on the molecular mechanism and mode of action of circRNAs as competitive endogenous RNAs in the muscle development of livestock and poultry. In addition, we also discuss the regulatory mechanism of circRNAs on muscle development in livestock in terms of transcription, translation, and mRNAs. The purpose of this article is to discuss the multiple regulatory roles of circRNAs in the process of muscle development in livestock, to provide new ideas for the development of a new co-expression regulation network, and to lay a foundation for enriching livestock breeding and improving livestock economic traits.

Functional Role of circRNAs in the Regulation of Fetal Development, Muscle Development, and Lactation in Livestock

circRNAs are a class of endogenous noncoding RNA molecules with closed loop structures. They are mainly responsible for regulating gene expression in eukaryotic cells. With the emergence of high-throughput RNA sequencing (RNA-Seq) and new types of bioinformatics tools, thousands of circRNAs have been discovered, making circRNA one of the research hotspots. Recent studies have shown that circRNAs play an important regulatory role in the growth, reproduction, and formation of livestock products. They can not only regulate mammalian fetal growth and development but also have important regulatory effects on livestock muscle development and lactation. In this review, we briefly introduce the putative biogenic pathways and regulatory functions of circRNA and highlight our understanding of circRNA and its latest advances in fetal development, muscle development, and lactation biogenesis as well as expression in livestock. This review will provide a theoretical basis for the research and development of related industries.

Circular RNA in cardiovascular disease: Expression, mechanisms and clinical prospects

Circular RNAs (circRNAs) are a group of covalently closed, endogenous, non‐coding RNAs, which exist widely in human tissues including the heart. Increasing evidence has shown that cardiac circRNAs play crucial regulatory roles in cardiovascular diseases (CVDs). In this review, we aimed to provide a systemic understanding of circRNAs with a special emphasis on the cardiovascular system. We have summarized the current research on the classification, biogenesis and properties of circRNAs as well as their participation in the pathogenesis of CVDs. CircRNAs are conserved, stable and have specific spatiotemporal expression; thus, they have been accepted as a potential diagnostic marker or an incremental prognostic biomarker for CVDs.

Identification of Candidate Circular RNAs Underlying Intramuscular Fat Content in the Donkey

Intramuscular fat (IMF) content is a crucial indicator of meat quality. Circular RNAs (circRNAs) are a large class of endogenous RNAs that are involved in many physiological processes. However, the expression and function of circRNA in IMF in the donkey remains unresolved. Here we performed an expression profiling of circRNAs in the donkey longissimus dorsi muscle and identified 12,727 candidate circRNAs. Among these, 70% were derived from the exons of protein genes. Furthermore, a total of 127 differentially expressed (DE) circRNAs were identified in high (H) and low (L) IMF content groups, including 63 upregulated and 64 downregulated circRNAs. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of the host genes of the DE circRNAs showed that the host genes were enriched in lipid metabolism related GO terms (e.g., fatty acid beta-oxidation using acyl-CoA dehydrogenase and MLL3/4 complex), and signaling pathways (e.g., TGF-beta and lysine degradation signaling pathway). Further analyses indicated that 127 DE circRNAs were predicted to potentially interact with miRNAs, leading to the construction of circRNA-miRNA regulatory network. Multiple circRNAs can potentially function as sponges of miRNAs that regulate the differentiation of adipocytes. Our results provide valuable expression profile information for circRNA in the donkey and new insight into the regulatory mechanisms of circRNAs in the regulation of IMF content.

Oxidized low-density lipoprotein accelerates the injury of endothelial cells via circ-USP36/miR-98-5p/VCAM1 axis

Circular RNAs (circRNAs) are a group of RNAs featured by a covalently closed continuous loop structure. This study aimed to uncover the function and mechanism of circ-ubiquitin specific peptidase 36 (USP36) in endothelial cells treated with oxidized low-density lipoprotein (ox-LDL). The levels of circ-USP36, microRNA-98-5p (miR-98-5p) and vascular cell adhesion molecule 1 (VCAM1) were examined by a quantitative real-time polymerase chain reaction (qRT-PCR). The viability, apoptosis and inflammation were detected by (4,5-Dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) assay, flow cytometry and enzyme-linked immunosorbent assay (ELISA), respectively. Western blot assay was performed to detect the expression of apoptosis and proliferation-related markers and VCAM1 protein level. The targets of circ-USP36 and miR-98-5p were searched using starBase website, and dual-luciferase reporter assay and RNA immunoprecipitation (RIP) assay were applied to validate the above predictions. Ox-LDL exposure induced the upregulation of circ-USP36 in HUVEC cells. Circ-USP36 accelerated ox-LDL-induced apoptosis, inflammatory and viability inhibition of HUVEC cells. MiR-98-5p was a direct downstream gene of circ-USP36. Circ-USP36 promoted the injury of ox-LDL-induced HUVEC cells through targeting miR-98-5p. VCAM1 could bind to miR-98-5p, and the protective effects of miR-98-5p accumulation on ox-LDL-induced HUVEC cells were reversed by the transfection of VCAM1. VCAM1 was regulated by circ-USP36/miR-98-5p signaling in HUVEC cells. Ox-LDL promoted the apoptosis and inflammation but suppressed the viability of HUVEC cells through upregulating circ-USP36, thus elevating the expression of VCAM1 via miR-98-5p.

LncRNA IMFlnc1 promotes porcine intramuscular adipocyte adipogenesis by sponging miR-199a-5p to up-regulate CAV-1

Background

Local Chinese local pig breeds have thinner muscle fiber and higher intramuscular-fat (IMF) content. But its regulation mechanism has not been discussed in-depth. Studies indicated that long non coding RNAs (lncRNAs) play important role in muscle and fat development.

Results

The lncRNAs expressional differences in the longissimus dorsi (LD) muscle were identified between Huainan pigs (local Chinese pigs, fat-type, HN) and Large White pigs (lean-type, LW) at 38, 58, and 78 days post conception (dpc). In total, 2131 novel lncRNAs were identified in 18 samples, and 291, 305, and 683 differentially expressed lncRNAs (DELs) were found between these two breeds at three stages, respectively. The mRNAs that co-expressed with these DELs were used for GO and KEGG analysis, and the results showed that muscle development and energy metabolism were more active at 58 dpc in HN, but at 78 dpc in LW pigs. Muscle cell differentiation and myofibril assembly might associated with earlier myogenesis and primary-muscle-fiber assembly in HN, and cell proliferation, insulin, and the MAPK pathway might be contribute to longer proliferation and elevated energy metabolism in LW pigs at 78 dpc. The PI3K/Akt and cAMP pathways were associated with higher IMF deposition in HN. Intramuscular fat deposition-associated long noncoding RNA 1 (IMFlnc1) was selected for functional verification, and results indicated that it regulated the expressional level of caveolin-1 (CAV-1) by acting as competing endogenous RNA (ceRNA) to sponge miR-199a-5p.

Conclusions

Our data contributed to understanding the role of lncRNAs in porcine-muscle development and IMF deposition, and provided valuable information for improving pig-meat quality.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12860-020-00324-8.

Genomewide analysis of circular RNA in pituitaries of normal and heat-stressed sows

BACKGROUND

As a newly characterized type of noncoding RNA, circular RNA (circRNA) has been shown to have functions in diverse biological processes of animals. It has been reported that several noncoding RNAs may regulate animals' response to heat stress which can be easily induced by hyperthermia in summer. However, the expression and functions of circRNAs in the pituitary of sows and whether they participate in heat stress adaption are still unclear.

RESULTS

In this study, we found that high temperature over the thermoneutral zone of sows during the summer increased the serum heat shock protein 70 (HSP70) level, decreased the superoxide dismutase (SOD) vitality and prolactin (PRL) concentration, and induced heat stress in sows. Then, we explored circRNA in the pituitary of heat-stressed and normal sows using RNA sequencing and bioinformatics analysis. In total, 12,035 circRNAs were detected, with 59 circRNAs differentially expressed, including 42 up-regulated and 17 down-regulated circRNAs in pituitaries of the heat-stressed sows. Six randomly selected circRNAs were identified through reverse transcription PCR followed by DNA sequencing and other 7 randomly selected differentially expressed circRNAs were verified by quantitative real-time PCR analysis. The predicted target genes regulated by circRNAs through sponging microRNAs (miRNAs) were enriched in metabolic pathway. Furthermore, the predicted circRNA-miRNA-mRNA interactions showed that some circRNAs might sponge miRNAs to regulate pituitary-specific genes and heat shock protein family members, indicating circRNA's roles in pituitary hormone secretion and heat stress response.

CONCLUSIONS

Our results provided a meaningful reference to understand the functions of circRNA in the porcine pituitary and the mechanisms by which circRNA may participate in animals' response to heat stress.

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.

Comparison of microRNAs in adipose and muscle tissue from seven indigenous Chinese breeds and Yorkshire pigs

Elucidation of the pig microRNAome is essential for interpreting functional elements of the genome and understanding the genetic architecture of complex traits. Here, we extracted small RNAs from skeletal muscle and adipose tissue, and we compared their expression levels between one Western breed (Yorkshire) and seven indigenous Chinese breeds. We detected the expression of 172 known porcine microRNAs (miRNAs) and 181 novel miRNAs. Differential expression analysis found 92 and 12 differentially expressed miRNAs in adipose and muscle tissue respectively. We found that different Chinese breeds shared common directional miRNA expression changes compared to Yorkshire pigs. Some miRNAs differentially expressed across multiple Chinese breeds, including ssc-miR-129-5p, ssc-miR-30 and ssc-miR-150, are involved in adipose tissue function. Functional enrichment analysis revealed that the target genes of the differentially expressed miRNAs are associated mainly with signaling pathways rather than metabolic and biosynthetic processes. The miRNA-target gene and miRNA-phenotypic traits networks identified many hub miRNAs that regulate a large number of target genes or phenotypic traits. Specifically, we found that intramuscular fat content is regulated by the greatest number of miRNAs in muscle tissue. This study provides valuable new candidate miRNAs that will aid in the improvement of meat quality and production.

Characterization of Circular RNAs in Chinese Buffalo (Bubalus bubalis) Adipose Tissue: A Focus on Circular RNAs Involved in Fat Deposition

Circular RNAs (circRNAs) have been identified as a novel type of regulators involved in multiple biological processes. However, circRNAs with a potential function in fat deposition in buffalo are poorly understood. In this study, six RNA libraries of adipose tissue were constructed for three young and three adult Chinese buffaloes with paired-ends RNA sequencing using the Illumina HiSeq 3000 platform. A total of 5141 circRNAs were computationally identified. Among them, 252 circRNAs were differentially expressed (DE) between the young and adult buffaloes. Of these, 54 were upregulated and 198 were downregulated in the adult group. Eight DE circRNAs were further identified by quantitative real-time-PCR (qRT-PCR) and Sanger sequencing. Co-expression analysis revealed that 34 circRNAs demonstrated a strong correlation with fat deposition-associated genes (|r| > 0.980). Among these, expressional correlation between two circRNAs (19:45387150|45389986 and 21:6969877|69753491) and PR/SET domain 16 was further verified using qRT-PCR, and a strong correlation was revealed (1 > |r| > 0.8). These results strongly suggest that circRNAs 19:45387150|45389986 and 21:6969877|69753491 are potential regulators of buffalo fat deposition. In summary, this study characterized the circRNA profiles of adipose tissues at different stages for the first time and revealed two circRNAs strongly correlated with fat deposition-associated genes, which provided new candidate regulators for fat deposition in buffalo.

Genome-wide identification and comparison of mRNAs, lncRNAs and circRNAs in porcine intramuscular, subcutaneous, retroperitoneal and mesenteric adipose tissues

Many types of RNAs, including messenger RNAs (mRNAs), long noncoding RNAs (lncRNAs) and circular RNAs (circRNAs), play crucial roles in regulating fat cell differentiation and tissue development. However, the expression profiles of these RNAs in different adipose tissues are still largely unknown. To shed light on this issue, we performed a transcriptome analysis of mRNAs, lncRNAs and circRNAs obtained from intramuscular adipose tissue, subcutaneous adipose tissue, retroperitoneal adipose tissue and mesenteric adipose tissue of Chinese Erhualian pigs. A number of differentially expressed mRNAs, lncRNAs and circRNAs were identified among the four adipose tissues. Tissue-specific analysis indicated that circRNAs exhibited the highest tissue specificity among mRNAs, lncRNAs and circRNAs, whereas intramuscular adipose tissue had the most tissue-specific genes among the four adipose tissues. Gene Ontology analysis showed that differentially expressed mRNAs among groups were involved mainly in lipid metabolism and immune inflammatory response processes. Furthermore, the co-expression network construction of mRNAs-lncRNAs revealed that several important lncRNAs, such as MSTRG.426159 and MSTRG.604206, might associate with lipid metabolic process. Taken together, these data provide a genome-wide resource of mRNAs, lncRNAs and circRNAs potentially involved in porcine fat metabolism, thus improving understanding of their function in diverse adipose tissues.

Up-regulation of circular RNA hsa_circ_0037909 promotes essential hypertension

Aims

Essential hypertension (EH) is a high prevalence disease facing a public health challenge. People were little known about the genetics of diagnosing the cause of EH. Circular RNAs that have a continuous cycle of covalent closure, without affected by RNA exonuclease, and are more stable and hard to degrade may involve into the molecule regulation mechanism of EH as an important biomedical.

Methods

qRT‐PCR was used to analyze circRNAs in total volume of human blood and the induced human aortic endothelial cells (HAECs) and human umbilical vein endothelial cells (HUVECs). Our case‐control study was involved with 48 pairs of case controls with sex and age (±3 years) match. We conducted t test, Pearson's χ² test, and receiver operating characteristics (ROC) curve analysis for the corresponding analysis.

Results

The expression level of hsa_circ_0037909 in EH patients was significantly higher than that in the healthy controls (P = 0.007), and the expression level of hsa‐miR‐637 in EH patients was significantly lower in than that in the healthy controls (P = 0.039); the same result appears in the HAECs and HUVECs. Hsa‐miR‐637 (adjusted P = 0.018), hsa_circ_0037909 (adjusted P = 0.005), HDL (adjusted P = 0.024), and serum creatinine (adjusted P = 0.014) were brought into the model which performed logistic regression analysis. The combination of two RNAs was excellent (P < 0.001) through ROC curve analysis. Hsa_circ_0037909 was significantly positively correlated with serum creatinine (P < 0.001) and low‐density lipoprotein (LDL) (P = 0.017).

Conclusions

Our findings suggested that the combination of hsa_circ_0037911 and hsa‐miR‐637 may be a significant important biomarker for early diagnosis of EH. Hsa_circ_0037909 may affect serum creatinine or LDL leading to the formation of EH.

Expression profiles of circular RNAs in sheep skeletal muscle

Objective

Circular RNAs (circRNAs) are a newfound class of non-coding RNA in animals and plants. Recent studies have revealed that circRNAs play important roles in cell proliferation, differentiation, autophagy and apoptosis during development. However, there are few reports about muscle development-related circRNAs in livestock.

Methods

RNA sequencing analysis was employed to identify and annotate circRNAs from longissimus dorsi of sheep. Reverse transcription followed by real-time quantitative (q) polymerase chain reaction (PCR) analysis verified the presence of these circRNAs. Targetscan7.0 and miRanda were used to analyse the interaction of circRNA-microRNA (miRNA). To investigate the function of circRNAs, an experiment was conducted to perform enrichment analysis hosting genes of circRNAs using gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) pathways.

Results

About 75.5 million sequences were obtained from RNA libraries of sheep skeletal muscle. These sequences were mapped to 729 genes in the sheep reference genome. We identified 886 circRNAs, including numerous circular intronic RNAs and exonic circRNAs. Reverse transcription PCR (RT-PCR) and DNA sequencing analysis confirmed the presence of several circRNAs. Real-Time RT-PCR analysis exhibited resistance of sheep circRNAs to RNase R digestion. We found that many circRNAs interacted with muscle-specific miRNAs involved in growth and development of muscle, especially circ776. The GO and KEGG enrichment analysis showed that hosting genes of circRNAs was involved in muscle cell development and signaling pathway.

Conclusion

The study provides comprehensive expression profiles of circRNAs in sheep skeletal muscle. Our study offers a large number of circRNAs to facilitate a better understanding of their roles in muscle growth. Meanwhile, we suggested that circ776 could be analyzed in future study.

Genome-wide analysis of circular RNAs in prenatal and postnatal muscle of sheep

Circular RNAs (circRNAs), a type of non-coding RNA with circular structure, were generated by back splicing and widely expressed in animals and plants. Recent studies have shown that circRNAs extensively participate in cell proliferation, cell differentiation, cell autophagy and other biological processes. However, the role and expression of circRNAs in the development and growth of muscle have not been studied in sheep. In our study, we first used RNA-seq to study the circRNAs in prenatal and postnatal longissimus dorsi muscle of sheep. A total of 6113 circRNAs were detected from the RNA-seq data. Several circRNAs were identified using reverse transcription PCR, DNA sequencing and RNase R digestion experiments. The expression levels of several circRNAs in prenatal and postnatal muscle were confirmed by Real-Time RT-PCR. The gene ontology (GO) and KEGG enrichment analysis of the host gene of the circRNAs showed that these circRNAs were mainly involved in the growth and development of muscle related signaling pathways. These circRNAs might sponge microRNAs (miRNAs) in predicted circRNA-miRNA-mRNA networks. The circRNAs expression profiles in muscle provided an important reference for the study of circRNAs in sheep.