An atlas and analysis of bovine skeletal muscle long noncoding RNAs

Understanding lncRNAs: key regulators of myogenesis and lipogenesis in farm animals

Long non-coding RNAs (lncRNAs) are RNA molecules exceeding 200 nucleotides in length. Recent studies have demonstrated their involvement in regulating gene expression and various biological processes. Among these, myogenesis and lipogenesis are particularly important because of their direct effects on muscle development and fat deposition in farm animals. These processes are crucial for determining meat quality, growth rates, and overall economic value in animal husbandry. Although the specific mechanisms through which lncRNAs influence these pathways are still under investigation, further research into their roles in muscle and fat development is crucial for optimizing farm animal breeding strategies. Here, we review the characteristics of lncRNAs, including their biogenesis, localization, and structures, with a particular focus on their association with myogenesis and adipogenesis. This review seeks to establish a theoretical foundation for enhancing farm animal production. In particular, focusing on lncRNAs may reveal how these molecules can enhance the economic traits of farm animals, thereby contributing to the optimization of farm animal breeding processes.

Transcriptome analysis of norepinephrine-induced lipolysis in differentiated adipocytes of Bama pig

Sympathetic innervation of white adipose tissue (WAT) plays a key role in the regulation of lipid metabolism. Sympathetic activation promotes release of norepinephrine (NE), which binds to adrenergic receptors on adipocytes, promoting adipocyte lipolysis and enhanced oxidative metabolism. However, the mechanism by which sympathetic nerves regulate lipid metabolism in pig adipose tissue remains unclear. We used NE to simulate the process of sympathetic driving in pig adipocytes. RNA sequencing (RNA-seq) was used to determine the gene expression profile of pig adipocytes responding to NE stimulation. Our data suggests that the lipolytic signaling pathway is activated in pig adipocytes upon acute stimulation of NE, resulting in enhanced lipid metabolism and lipolysis, consistent with the phenomena found in humans and mice. Specifically, differentially expressed protein coding genes (PCGs) (SIRT4, SLC27A1) are mainly associated with functions that inhibit fatty acid oxidation and promote lipid synthesis. Similarly, we investigated the changes in regulatory transcripts such as long non-coding RNAs (lncRNAs) and transcripts of uncertain coding potential (TUCP) in response to NE and found that differentially expressed lncRNAs (lncG47338, lncG30660, lncG29516, lncG3790) and TUCP (TUCP_G38001) were co-expressed with target genes related to the promotion of fatty acid β-oxidation, lipolysis and oxidative metabolism, thus acting as regulators. These results indicate a broad suite of gene expression alterations in response to NE stimulation and promote the understanding of the molecular mechanisms by which NE regulates lipid metabolism in pigs.

Buffalo long non-coding RNA gene11007 promotes myoblasts proliferation

Buffalo meat is of good quality because it is lean and tender, and could bring significant cardiovascular benefits. The underlying difference in muscle development and meat quality is a complex and precisely orchestrated process which has been demonstrated to be regulated by long non-coding RNAs (lncRNAs). However, the regulatory role of lncRNAs in the growth and development of buffalo skeletal muscle is still unclear. In this study, the Ribo-Zero RNA-Seq method was used to explore the lncRNA expression profiles of buffalo myoblasts during the proliferation and differentiation phases. A specific set of 9,978 lncRNAs was found. By comparing the expression profiles of lncRNAs, it was found that there were 1,576 differentially expressed lncRNAs (DELs) during buffalo myoblast differentiation. Twelve DELs were chosen and subsequently verified in eight different buffalo tissues during fetal and adult stages by using qPCR. Gene11007 was found to be one of the most down-regulated lncRNAs during buffalo myoblasts differentiation and it was subsequently characterized. EdU, CCK-8, qPCR and western blotting assays showed that gene11007 promoted the proliferation of buffalo myoblasts but it had no effect on cell differentiation. Our research may enrich the genome annotations of buffalo and provide a new molecular target for the in-depth understanding of the regulation of lncRNAs in skeletal muscle.

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.

Transcriptional atlas analysis from multiple tissues reveals the expression specificity patterns in beef cattle

BACKGROUND

A comprehensive analysis of gene expression profiling across tissues can provide necessary information for an in-depth understanding of their biological functions. We performed a large-scale gene expression analysis and generated a high-resolution atlas of the transcriptome in beef cattle.

RESULTS

Our transcriptome atlas was generated from 135 bovine tissues in adult beef cattle, covering 51 tissue types of major organ systems (e.g., muscular system, digestive system, immune system, reproductive system). Approximately 94.76% of sequencing reads were successfully mapped to the reference genome assembly ARS-UCD1.2. We detected a total of 60,488 transcripts, and 32% of them were not reported before. We identified 2654 housekeeping genes (HKGs) and 477 tissue-specific genes (TSGs) across tissues. Using weighted gene co-expression network analysis, we obtained 24 modules with 237 hub genes (HUBGs). Functional enrichment analysis showed that HKGs mainly maintain the basic biological activities of cells, while TSGs were involved in tissue differentiation and specific physiological processes. HKGs in bovine tissues were more conserved in terms of expression pattern as compared to TSGs and HUBGs among multiple species. Finally, we obtained a subset of tissue-specific differentially expressed genes (DEGs) between beef and dairy cattle and several functional pathways, which may be involved in production and health traits.

CONCLUSIONS

We generated a large-scale gene expression atlas across the major tissues in beef cattle, providing valuable information for enhancing genome assembly and annotation. HKGs, TSGs, and HUBGs further contribute to better understanding the biology and evolution of multiple tissues in cattle. DEGs between beef and dairy cattle also fill in the knowledge gaps about differential transcriptome regulation of bovine tissues underlying economically important traits.

Differential Expression Profiles of lncRNA Following LPS-Induced Inflammation in Bovine Mammary Epithelial Cells

Bovine mastitis is an inflammatory response of mammary glands caused by pathogenic microorganisms such as Escherichia coli (E. coli). As a key virulence factor of E. coli, lipopolysaccharide (LPS) triggers innate immune responses via activation of the toll-like-receptor 4 (TLR4) signaling pathway. However, the molecular regulatory network of LPS-induced bovine mastitis has yet to be fully mapped. In this study, bovine mammary epithelial cell lines MAC-T were exposed to LPS for 0, 6 and 12 h to assess the expression profiles of long non-coding RNAs (lncRNAs) using RNA-seq. Differentially expressed lncRNAs (DElncRNAs) were filtered out of the raw data for subsequent analyses. A total of 2,257 lncRNAs, including 210 annotated and 2047 novel lncRNAs were detected in all samples. A large proportion of lncRNAs were present in a high abundance, and 112 DElncRNAs were screened out at different time points. Compared with 0 h, there were 22 up- and 25 down-regulated lncRNAs in the 6 h of post-infection (hpi) group, and 27 up- and 22 down-regulated lncRNAs in the 12 hpi group. Compared with the 6 hpi group, 32 lncRNAs were up-regulated and 25 lncRNAs were down-regulated in the 12 hpi group. These DElncRNAs are involved in the regulation of a variety of immune-related processes including inflammatory responses bMECs exposed to LPS. Furthermore, lncRNA TCONS_00039271 and TCONS_00139850 were respectively significance down- and up-regulated, and their target genes involve in regulating inflammation-related signaling pathways (i.e.,Notch, NF-κB, MAPK, PI3K-Akt and mTOR signaling pathway), thereby regulating the occurrence and development of E. coli mastitis. This study provides a resource for lncRNA research on the molecular regulation of bovine mastitis.

Testis transcriptome profiling identified lncRNAs involved in spermatogenic arrest of cattleyak

Cattleyaks are the crossbred offspring between cattle and yaks, exhibiting the prominent adaptability to the harsh environment as yaks and much higher growth performances than yaks around Qinghai-Tibet plateau. Unfortunately, cattleyak cannot be effectively used in yak breeding due to its male infertility resulted from spermatogenic arrest. In this study, we performed RNA sequencing (RNA-seq) and bioinformatics analysis to determine the expression profiles of long noncoding RNA (lncRNA) from cattleyak and yak testis. A total of 604 differentially expressed (DE) lncRNAs (135 upregulated and 469 downregulated) were identified in cattleyak with respect to yak. Through gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses, we identified several DE lncRNAs regulating the mitotic cell cycle processes by targeting the genes significantly associated with the mitotic cell cycle checkpoint and DNA damage checkpoint term and also significantly involved in p53 signaling pathway, mismatch repair and homologous recombination pathway (P < 0.05). The reverse transcription PCR (RT-PCR) and quantitative Real-Time PCR (qRT-PCR) analysis of the randomly selected fourteen DE lncRNAs and the seven target genes validated the RNA-seq data and their true expressions during spermatogenesis in vivo. Molecular cloning and sequencing indicated that the testis lncRNAs NONBTAT012170 and NONBTAT010258 presented higher similarity among different cattleyak and yak individuals. The downregulation of these target genes in cattleyak contributed to the abnormal DNA replication and spermatogenic arrest during the S phase of mitotic cell cycle. This study provided a novel insight into lncRNA expression profile changes associated with spermatogenic arrest of cattleyak.

A novel lncRNA promotes myogenesis of bovine skeletal muscle satellite cells via PFN1-RhoA/Rac1

Myogenesis, the process of skeletal muscle formation, is a highly coordinated multistep biological process. Accumulating evidence suggests that long non-coding RNAs (lncRNAs) are emerging as a gatekeeper in myogenesis. Up to now, most studies on muscle development-related lncRNAs are mainly focussed on humans and mice. In this study, a novel muscle highly expressed lncRNA, named lnc23, localized in nucleus, was found differentially expressed in different stages of embryonic development and myogenic differentiation. The knockdown and over-expression experiments showed that lnc23 positively regulated the myogenic differentiation of bovine skeletal muscle satellite cells. Then, TMT 10-plex labelling quantitative proteomics was performed to screen the potentially regulatory proteins of lnc23. Results indicated that lnc23 was involved in the key processes of myogenic differentiation such as cell fusion, further demonstrated that down-regulation of lnc23 may inhibit myogenic differentiation by reducing signal transduction and cell fusion among cells. Furthermore, RNA pulldown/LC-MS and RIP experiment illustrated that PFN1 was a binding protein of lnc23. Further, we also found that lnc23 positively regulated the protein expression of RhoA and Rac1, and PFN1 may negatively regulate myogenic differentiation and the expression of its interacting proteins RhoA and Rac1. Hence, we support that lnc23 may reduce the inhibiting effect of PFN1 on RhoA and Rac1 by binding to PFN1, thereby promoting myogenic differentiation. In short, the novel identified lnc23 promotes myogenesis of bovine skeletal muscle satellite cells via PFN1-RhoA/Rac1.

Transcriptome analysis of long noncoding RNAs ribonucleic acids from the livers of Hu sheep with different residual feed intake

Long noncoding RNAs (LncRNAs), as key regulators, have vital functions in various biological activities. However, in sheep, little has been reported concerning the genetic mechanism of LncRNA regulation of feed efficiency. In the present study, we explored the genome-wide expression of LncRNAs and transcripts of uncertain coding potential (TUCPs) in the livers of sheep with extreme residual feed intake (RFI) using RNA sequencing. We identified 1 523 TUCPs and 1 996 LncRNAs, among which 10 LncRNAs and 16 TUCPs were identified as being differentially expressed between the High-RFI and Low-RFI groups. Co-expression and co-localization methods were used to search for LncRNA and TUCP target genes, which identified 970/1 538 and 23/27 genes, respectively. Ontology and pathways analysis revealed that the LncRNAs/TUCPs that were highly expressed in the Low-RFI group are mostly concentrated in energy metabolism pathways. For example, LNC_000890 and TUCP_000582 might regulate liver tissue metabolic efficiency. The LncRNAs/TUCPs that were highly expressed in the High-RFI group are mostly enriched in immune function pathways. For example, TUCP_000832 might regulate animal health, thereby affecting feed efficiency. Subsequently, a co-expression network was established by applying the expression information of both the differentially expressed LncRNAs and TUCPs and their target mRNAs. The network indicated that differentially expressed genes targeted by the upregulated LncRNAs and TUCPs were mainly related to energy metabolism, while those genes targeted by the downregulated LncRNAs and TUCPs were mainly related to immune response. These results provide the basis for further study of LncRNA/TUCP-mediated regulation of feed efficiency.

Identification and Comparative Analysis of Long Non-Coding RNA in the Skeletal Muscle of Two Dezhou Donkey Strains

Long non-coding RNA (lncRNA) has been extensively studied in many livestock. However, compared with other livestock breeds, there is less research regarding donkey lncRNA function. It has been reported that lncRNA plays an important role in the timing control of development, aging, and death of livestock. Therefore, the study of donkey skeletal muscle lncRNA is of great significance for exploring donkey meat production performance. In this study, RNA-Seq was used to perform high-throughput sequencing of skeletal muscle (longissimus dorsi and gluteus) of two Dezhou donkey strains (SanFen and WuTou). The differentially expressed lncRNAs were screened between different strains and tissues. Then candidate genes for conjoint analysis were screened based on Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. Finally, the accuracy of the sequencing data was verified by real-time quantitative polymerase chain reaction (RT-qPCR). Herein, we identified 3869 novel lncRNAs and 73 differentially expressed lncRNAs. Through the comparison between groups, the specific expression of lncRNAs were found in different strains and muscle tissues. Importantly, we constructed the lncRNA-miRNA-mRNA interaction network and found three important candidate lncRNAs (MSTRG.9787.1, MSTRG.3144.1, and MSTRG.9886.1) and four candidate genes (ACTN1, CDON, FMOD, and BMPR1B). Analysis of the KEGG pathway indicates that the TGF-β signaling pathway plays a pivotal role in the growth and development of skeletal muscle in Dezhou donkey strains.

Comparison of characteristics of long noncoding RNA in Hanwoo according to sex

Objective

Cattle were some of the first animals domesticated by humans for the production of milk, meat, etc. Long noncoding RNA (lncRNA) is defined as longer than 200 bp in non-protein coding transcripts. lncRNA is known to function in regulating gene expression and is currently being studied in a variety of livestock including cattle. The purpose of this study is to analyze the characteristics of lncRNA according to sex in Hanwoo cattle.

Methods

This study was conducted using the skeletal muscles of 9 Hanwoo cattle include bulls, steers and cows. RNA was extracted from skeletal muscle of Hanwoo. Sequencing was conducted using Illumina HiSeq2000 and mapped to the Bovine Taurus genome. The expression levels of lncRNAs were measured by DEGseq and quantitative trait loci (QTL) data base was used to identify QTLs associated with lncRNA. The python script was used to match the nearby genes

Results

In this study, the expression patterns of transcripts of bulls, steers and cows were identified. And we identified significantly differentially expressed lncRNAs in bulls, steers and cows. In addition, characteristics of lncRNA which express differentially in muscles according to the sex of Hanwoo were identified. As a result, we found differentially expressed lncRNAs according to sex were related to shear force and body weight.

Conclusion

This study was classified and characterized lncRNA which differentially expressed by sex in Hanwoo cattle. We believe that the characterization of lncRNA by sex of Hanwoo will be helpful for future studies of the physiological mechanisms of Hanwoo cattle.

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.

Novel lncRNA lncFAM200B: Molecular Characteristics and Effects of Genetic Variants on Promoter Activity and Cattle Body Measurement Traits

Skeletal muscle is one of the three major muscle types in an organism and has key roles in the motor system, metabolism, and homeostasis. RNA-Seq analysis showed that novel lncRNA, lncFAM200B, was differentially expressed in embryonic, neonatal, and adult cattle skeletal muscles. The main aim of this study was to investigate the molecular and expression characteristics of lncFAM200B along with its crucial genetic variations. Our results showed that bovine lncFAM200B was a 472 nucleotide (nt) non-coding RNA containing two exons. The transcription factor binding site prediction analysis found that lncFAM200B promoter region was enriched with SP1 transcription factor, which promotes the binding of myogenic regulatory factor MyoD and DNA sequence. The mRNA expression analysis showed that lncFAM200B was differentially expressed in embryonic, neonatal, adult bovine muscle tissues, and the lncFAM200B expression trend positively correlated with that of MyoG and Myf5 in myoblast proliferation and differential stages. To identify the promoter active region of lncFAM200B, we constructed promoter luciferase reporter gene vector pGL3-Basic plasmids containing lncFAM200B promoter sequences and transfected them into 293T, C2C12, and 3T3-L1 cells. Our results suggested that lncFAM200B promoter active region was from −403 to −139 (264 nt) of its transcription start site, covering 6 SP1 potential binding sites. Furthermore, we found a novel C-T variation, named as SNP2 (ERZ990081 in European Variation Archive) in the promoter active region, which was linked to the nearby SNP1 (rs456951291 in Ensembl database). The genotypes of SNP1 and combined genotypes of SNP1 and SNP2 were significantly associated with Jinnan cattle hip height. The luciferase activity analysis found that the SNP1-SNP2 haplotype CC had the highest luciferase activity, which was consistent with the association analysis result that the combined genotype CC-CC carriers had the highest hip height in Jinnan cattle. In conclusion, our data showed that lncFAM200B is a positive regulator of muscle development and that SNP1 and SNP2 could be used as genetic markers for marker-assisted selection (MAS) breeding of beef cattle.

lnc9141-a and -b Play a Different Role in Bovine Myoblast Proliferation, Apoptosis, and Differentiation

Previously, our transcriptome sequencing revealed that lnc9141 was differentially expressed in muscles of fetal bovine, calf, and adult bovine, which is considered to provide the basis for raising the muscle mass. In this study, we identified lnc9141 characters. lnc9141 has different transcription start sites and 3′ alternative splicing sites of exon 1, producing lnc9141-a and lnc9141-b transcripts that were highly expressed in the heart and lung. Moreover, neither lnc9141-a nor lnc9141-b had the ability to encode proteins. The functions of lnc9141-a and lnc9141-b were explored by cell cycle, 5-ethynyl-2'-deoxyuridine (EdU), and 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) assay. The results showed that lnc9141-a or lnc9141-b overexpression decreased the number of myoblasts in the S phase and increased the proportion of cells in the G₀/G₁ phase. Furthermore, overexpressing lnc9141-a and lnc9141-b respectively downregulated the expression of Cyclin D1. However, lnc9141-a or lnc9141-b interference was found to increase the number of S-phase myoblasts, and upregulate Cyclin D1 and Cyclin E expression. Through Annexin V-FITC/propidium iodide (PI) double staining and the expression of apoptosis marker genes (Bax, Bcl2, and Caspase-3), it was found that lnc9141-b could regulate the expression of Bax gene. Meantime, high expression of lnc9141-b could decrease MyHC expression. In addition, the intergenic region between lnc9141 and IRX5 was 2.3 kb, with a head-to-head orientation. The study also revealed the core regions of the lnc9141 and IRX5 promoter. Our study demonstrated that both lnc9141-a and -b expression inhibited bovine myoblast proliferation. However, lnc9141-b regulated Bax and MyHC expression. The regulatory mechanism of lnc9141-a and lnc9141-b needs to be further explored.

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.

Gene Expression Maps in Plants: Current State and Prospects

For many years, progress in the identification of gene functions has been based on classical genetic approaches. However, considerable recent omics developments have brought to the fore indirect but high-resolution methods of gene function identification such as transcriptomics, proteomics, and metabolomics. A transcriptome map is a powerful source of functional information and the result of the genome-wide expression analysis of a broad sampling of tissues and/or organs from different developmental stages and/or environmental conditions. In plant science, the application of transcriptome maps extends from the inference of gene regulatory networks to evolutionary studies. However, only some of these data have been integrated into databases, thus enabling analyses to be conducted without raw data; without this integration, extensive data preprocessing is required, which limits data usability. In this review, we summarize the state of plant transcriptome maps, analyze the problems associated with the combined analysis of large-scale data from various studies, and outline possible solutions to these problems.

Transcriptomic analysis to affirm the regulatory role of long non-coding RNA in horn cancer of Indian zebu cattle breed Kankrej (Bos indicus)

Long non-coding RNA (lncRNA) was previously considered as a non-functional transcript, which now established as part of regulatory elements of biological events such as chromosome structure, remodeling, and regulation of gene expression. The study presented here showed the role of lncRNA through differential expression analysis on cancer-related coding genes in horn squamous cell carcinoma of Indian zebu cattle. A total of 10,360 candidate lncRNAs were identified and further analyzed for its coding potential ability using three tools (CPC, CPAT, and PLEK) that provide 8862 common lncRNAs. Pfam analysis of these common lncRNAs gave 8612 potential candidates for lncRNA differential expression analysis. Differential expression analysis showed a total of 59 significantly differentially expressed genes and 19 lncRNAs. Pearson's correlation analysis was used to identify co-expressed mRNA-lncRNAs to established relation of the regulatory role of lncRNAs in horn cancer. We established a positive relation of seven upregulated (XLOC_000016, XLOC_002198, XLOC_002851, XLOC_ 007383, XLOC_010701, XLOC_010272, and XLOC_011517) and one downregulated (XLOC_011302) lncRNAs with eleven genes that are related to keratin family protein, keratin-associated protein family, cornifelin, corneodesmosin, serpin family protein, and metallothionein that have well-established role in squamous cell carcinoma through cellular communication, cell growth, cell invasion, and cell migration. These biological events were found to be related to the MAPK pathway of cell cycle regulation indicating the role of lncRNAs in manipulating cell cycle regulation during horn squamous cell carcinomas that will be useful in identifying molecular portraits related to the development of horn cancer.

A novel long non-coding RNA, lncKBTBD10, involved in bovine skeletal muscle myogenesis

Accumulating evidence suggests that long non-coding RNAs (lncRNAs) play a crucial role in regulating skeletal muscle myogenesis, a highly coordinated multistep biological process. However, most studies of lncRNAs have focused on humans, mouse, and other model animals. In this study, we identified a novel lncRNA, named lncKBTBD10, located in the nucleus and involved in the proliferation and differentiation of bovine skeletal muscle satellite cells. Prediction of coding potential and in vitro translation system illustrated that lncKBTBD10 has no encoding capability. With the process of myogenic differentiation, the expression of lncKBTBD10 gradually increased. To elucidate the functions of lncKBTBD10 during myogenesis, the gain/loss-of-function experiments were performed. Results showed that the proliferation and differentiation of bovine skeletal muscle satellite cells were all suppressed whether lncKBTBD10 was knocked down or over-expressed. Furthermore, we found that lncKBTBD10 may affect its proximity gene KBTBD10 to involve in myogenesis. Results indicated that the protein level of KBTBD10 was all diminished after induced differentiation for 2 d in differentiation medium (DM2) whether lncKBTBD10 was knocked down or over-expressed. It may support why the altering of lncKBTBD10 can suppress the proliferation and differentiation of bovine skeletal muscle satellite cells. In short, our study elucidated that lncKBTBD10 could induce a decrease of KBTBD10 protein and further to affect bovine skeletal muscle myogenesis. The novel identified lncKBTBD10 may provide a reference for lncRNAs involved in myogenesis of bovine skeletal muscle.

Genome-wide identification of tissue-specific long non-coding RNA in three farm animal species

Background

Numerous long non-coding RNAs (lncRNAs) have been identified and their roles in gene regulation in humans, mice, and other model organisms studied; however, far less research has been focused on lncRNAs in farm animal species. While previous studies in chickens, cattle, and pigs identified lncRNAs in specific developmental stages or differentially expressed under specific conditions in a limited number of tissues, more comprehensive identification of lncRNAs in these species is needed. The goal of the FAANG Consortium (Functional Annotation of Animal Genomes) is to functionally annotate animal genomes, including the annotation of lncRNAs. As one of the FAANG pilot projects, lncRNAs were identified across eight tissues in two adult male biological replicates from chickens, cattle, and pigs.

Results

Comprehensive lncRNA annotations for the chicken, cattle, and pig genomes were generated by utilizing RNA-seq from eight tissue types from two biological replicates per species at the adult developmental stage. A total of 9393 lncRNAs in chickens, 7235 lncRNAs in cattle, and 14,429 lncRNAs in pigs were identified. Including novel isoforms and lncRNAs from novel loci, 5288 novel lncRNAs were identified in chickens, 3732 in cattle, and 4870 in pigs. These transcripts match previously known patterns of lncRNAs, such as generally lower expression levels than mRNAs and higher tissue specificity. An analysis of lncRNA conservation across species identified a set of conserved lncRNAs with potential functions associated with chromatin structure and gene regulation. Tissue-specific lncRNAs were identified. Genes proximal to tissue-specific lncRNAs were enriched for GO terms associated with the tissue of origin, such as leukocyte activation in spleen.

Conclusions

LncRNAs were identified in three important farm animal species using eight tissues from adult individuals. About half of the identified lncRNAs were not previously reported in the NCBI annotations for these species. While lncRNAs are less conserved than protein-coding genes, a set of positionally conserved lncRNAs were identified among chickens, cattle, and pigs with potential functions related to chromatin structure and gene regulation. Tissue-specific lncRNAs have potential regulatory functions on genes enriched for tissue-specific GO terms. Future work will include epigenetic data from ChIP-seq experiments to further refine these annotations.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5037-7) contains supplementary material, which is available to authorized users.

Cross-species inference of long non-coding RNAs greatly expands the ruminant transcriptome

Background

mRNA-like long non-coding RNAs (lncRNAs) are a significant component of mammalian transcriptomes, although most are expressed only at low levels, with high tissue-specificity and/or at specific developmental stages. Thus, in many cases lncRNA detection by RNA-sequencing (RNA-seq) is compromised by stochastic sampling. To account for this and create a catalogue of ruminant lncRNAs, we compared de novo assembled lncRNAs derived from large RNA-seq datasets in transcriptional atlas projects for sheep and goats with previous lncRNAs assembled in cattle and human. We then combined the novel lncRNAs with the sheep transcriptional atlas to identify co-regulated sets of protein-coding and non-coding loci.

Results

Few lncRNAs could be reproducibly assembled from a single dataset, even with deep sequencing of the same tissues from multiple animals. Furthermore, there was little sequence overlap between lncRNAs that were assembled from pooled RNA-seq data. We combined positional conservation (synteny) with cross-species mapping of candidate lncRNAs to identify a consensus set of ruminant lncRNAs and then used the RNA-seq data to demonstrate detectable and reproducible expression in each species. In sheep, 20 to 30% of lncRNAs were located close to protein-coding genes with which they are strongly co-expressed, which is consistent with the evolutionary origin of some ncRNAs in enhancer sequences. Nevertheless, most of the lncRNAs are not co-expressed with neighbouring protein-coding genes.

Conclusions

Alongside substantially expanding the ruminant lncRNA repertoire, the outcomes of our analysis demonstrate that stochastic sampling can be partly overcome by combining RNA-seq datasets from related species. This has practical implications for the future discovery of lncRNAs in other species.

Electronic supplementary material

The online version of this article (10.1186/s12711-018-0391-0) contains supplementary material, which is available to authorized users.

Long noncoding RNAs are associated with metabolic and cellular processes in the jejunum mucosa of pre-weaning calves in response to different diets

Long noncoding RNAs (lncRNAs) emerged as important regulatory component of mechanisms involved in gene expression, chromatin modification and epigenetic processes, but they are rarely annotated in the bovine genome.

Our study monitored the jejunum transcriptome of German Holstein calves fed two different milk diets using transcriptome sequencing (RNA-seq). To identify potential lncRNAs within the pool of unknown transcripts, four bioinformatic lncRNA prediction tools were applied. The intersection of the alignment-free lncRNA prediction tools (CNCI, PLEK and FEELnc) predicted 1,812 lncRNA transcripts concordantly comprising a catalogue of 1,042 putative lncRNA loci expressed in the calves’ intestinal mucosa.

Nine lncRNA loci were differentially expressed (DE lncRNAs) between both calf groups. To elucidate their biological function, we applied a systems biology approach that combines weighted gene co-expression network analysis with functional enrichment and biological pathway analysis. Four DE lncRNAs were found to be strongly correlated with a gene network module (GNM) enriched for genes from canonical pathways of remodeling of epithelial adherens junction, tight junction and integrin signaling. Another DE lncRNA was strongly correlated with a GNM enriched for genes associated with energy metabolism and maintaining of cellular homeostasis with a focus on mitochondrial processes.

Our data suggest that these DE lncRNAs may play potential regulatory roles in modulating biological processes associated with energy metabolism pathways and cellular signaling processes affecting the barrier function of intestinal epithelial cells of calves in response to different feeding regimens in the pre-weaning period.

Integration of lncRNA and mRNA Transcriptome Analyses Reveals Genes and Pathways Potentially Involved in Calf Intestinal Growth and Development during the Early Weeks of Life

A better understanding of the factors that regulate growth and immune response of the gastrointestinal tract (GIT) of calves will promote informed management practices in calf rearing. This study aimed to explore genomics (messenger RNA (mRNA)) and epigenomics (long non-coding RNA (lncRNA)) mechanisms regulating the development of the rumen and ileum in calves. Thirty-two calves (≈5-days-old) were reared for 96 days following standard procedures. Sixteen calves were humanely euthanized on experiment day 33 (D33) (pre-weaning) and another 16 on D96 (post-weaning) for collection of ileum and rumen tissues. RNA from tissues was subjected to next generation sequencing and 3310 and 4217 mRNAs were differentially expressed (DE) between D33 and D96 in ileum and rumen tissues, respectively. Gene ontology and pathways enrichment of DE genes confirmed their roles in developmental processes, immunity and lipid metabolism. A total of 1568 (63 known and 1505 novel) and 4243 (88 known and 4155 novel) lncRNAs were detected in ileum and rumen tissues, respectively. Cis target gene analysis identified BMPR1A, an important gene for a GIT disease (juvenile polyposis syndrome) in humans, as a candidate cis target gene for lncRNAs in both tissues. LncRNA cis target gene enrichment suggested that lncRNAs might regulate growth and development in both tissues as well as posttranscriptional gene silencing by RNA or microRNA processing in rumen, or disease resistance mechanisms in ileum. This study provides a catalog of bovine lncRNAs and set a baseline for exploring their functions in calf GIT development.

Identification and characterization of long non-coding RNA in prenatal and postnatal skeletal muscle of sheep

lncRNAs are a class of transcriptional RNA molecules of >200 nucleotides in length. However, the overall expression pattern and function of lncRNAs in sheep muscle is not clear. Here, we identified 1566 lncRNAs and 404 differentially expressed lncRNAs in sheep muscle from prenatal (110 days of fetus) and postnatal (2 to 3 years old of adult sheep) developmental stages by using RNA-seq technology. Several lncRNAs were identified by using RT-PCR and DNA sequencing. The expression levels of several lncRNAs were confirmed by qRT-PCR. We analyzed the effect of lncRNAs that act cis to the target genes. lncRNA targeting genes were involved in signaling pathways associated with growth and development of muscle by GO and KEGG enrichment analysis. Through our study, we provide a comprehensive expression profile of muscle lncRNAs in sheep, which provides valuable resources for further understanding genetic regulation of muscle growth and development from the perspective of lncRNA.