Comparison of long noncoding RNA between muscles and adipose tissues in Hanwoo beef cattle

Screening and identification of lncRNAs in preadipocyte differentiation in sheep

Studies of preadipocyte differentiation and fat deposition in sheep have mainly focused on functional genes, and with no emphasis placed on the role that long non-coding RNAs (lncRNAs) may have on the activity of those genes. Here, the expression profile of lncRNAs in ovine preadipocyte differentiation was investigated and the differentially expressed lncRNAs were screened on day 0 (D0), day 2(D2) and day 8(D8) of ovine preadipocyte differentiation, with their target genes being predicted. The competing endogenous RNA (ceRNA) regulatory network was constructed by GO and KEGG enrichment analysis for functional annotation, and some differentially expressed lncRNAs were randomly selected to verify the RNA-Seq results by RT-qPCR. In the study, a total of 2517 novel lncRNAs and 3943 known lncRNAs were identified from ovine preadipocytes at the three stages of differentiation, with the highest proportion being intergenic lncRNAs. A total of 3455 lncRNAs were expressed at all three stages of preadipocyte differentiation, while 214, 226 and 228 lncRNAs were uniquely expressed at day 0, day 2 and day 8, respectively. By comparing the expression of the lncRNAs between the three stages of differentiation stages, a total of 405, 272 and 359 differentially expressed lncRNAs were found in D0-vs-D2, D0-vs-D8, and D2-vs-D8, respectively. Functional analysis revealed that the differentially expressed lncRNAs were enriched in signaling pathways related to ovine preadipocyte differentiation, such as mitogen-activated protein kinase (MAPK) pathway, the phosphoinositide 3-kinase protein kinase B (PI3K-Akt) pathway, and the transforming growth factor beta (TGF-β) pathway. In summary, lncRNAs from preadipocytes at different stages of differentiation in sheep were identified and screened using RNA-Seq technology, and the regulatory mechanisms of lncRNAs in preadipocyte differentiation and lipid deposition were explored. This study provides a theoretical reference for revealing the roles of lncRNAs in ovine preadipocyte differentiation and also offers a theoretical basis for further understanding the regulatory mechanisms of ovine preadipocyte differentiation.

Multi-tissue transcriptome analysis to identify candidate genes associated with weight regulation in Hanwoo cattle

While genetic markers related to meat production traits have been identified in many other cattle breeds, research on weight in Hanwoo cattle (Korean native cattle) is still insufficient. In this study, we performed expression quantitative trait loci (eQTL) analysis and differential gene expression analysis to detect candidate genes influencing the weight characteristics of 32 castrated Hanwoo cattle across 22 tissues and, we identified variants that affect gene expression levels. In total, we identified a total of 3,298 differentially expressed genes, among which we discovered key genes such as UBD, RGS2, FASN, and SCD that have functions related to adipogenesis, body weight, obesity, and lipid metabolism. Gene-set enrichment analysis revealed that candidate genes in adipose tissue are involved in metabolic pathways linked to obesity-related traits, adipose metabolism, and lipid metabolism. Additionally, we found that decreased expression of TRIM31 contributes to weight gain which can be explained by the associated candidate cis-eQTL genotypes for TRIM31 and their effect on differential gene expression between the lower and higher weight groups. Our findings revealed candidate genes associated with the weight of Hanwoo cattle and perhaps can provide comprehensive insights into the association of weight with various tissues beyond adipose tissue and muscle, indicating the potential for expanding the focus of livestock trait research.

Research progress of non-coding RNAs regulation on intramuscular adipocytes in domestic animals

Intramuscular fat (IMF) is the main determinant of the economic value of domestic animal meat, and has a vital impact on the sensory quality characteristics, while the content of IMF is mainly determined by the size and number of intramuscular adipocytes. In recent years, due to the development of sequencing technology and omics technology, a large number of non-coding RNAs have been identified in intramuscular adipocytes. Non-coding RNAs are a kind of RNA regulatory factors with biological functions but without translation function, which mainly include microRNAs (miRNAs), long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs). These non-coding RNAs regulate the key genes of intramuscular adipocyte growth and development at post-transcriptional level through a variety of regulatory mechanisms, and affect the number and size of intramuscular adipocytes, thus affecting the content of IMF. Here, the review summarizes the candidate non-coding RNAs (miRNAs, lncRNAs, circRNAs) and genes involved in the regulation of intramuscular adipocytes, the related regulation mechanism and signaling pathways, in order to provide reference for further clarifying the molecular regulation mechanism of non-coding RNAs on intramuscular adipocytes in domestic animals.

Long Non-Coding RNA BNIP3 Inhibited the Proliferation of Bovine Intramuscular Preadipocytes via Cell Cycle

The intramuscular fat (or marbling fat) content is an essential economic trait of beef cattle and improves the flavor and palatability of meat. Several studies have highlighted the correlation between long non-coding RNAs (lncRNAs) and intramuscular fat development; however, the precise molecular mechanism remains unknown. Previously, through a high-throughput sequencing analysis, we found a lncRNA and named it a long non-coding RNA BNIP3 (lncBNIP3). The 5' RACE and 3' RACE explored 1945 bp total length of lncBNIP3, including 1621 bp of 5'RACE, and 464 bp of 3'RACE. The nucleoplasmic separation and FISH results explored the nuclear localization of lncBNIP3. Moreover, the tissue expression of lncBNIP3 was higher in the longissimus dorsi muscle, followed by intramuscular fat. Furthermore, down-regulation of lncBNIP3 increased the 5-Ethynyl-2'- deoxyuridine (EdU)-EdU-positive cells. The flow cytometry results showed that the number of cells in the S phase was significantly higher in preadipocytes transfected with si-lncBNIP3 than in the control group (si-NC). Similarly, CCK8 results showed that the number of cells after transfection of si-lncBNIP3 was significantly higher than in the control group. In addition, the mRNA expressions of proliferative marker genes CyclinB1 (CCNB1) and Proliferating Cell Nuclear Antigen (PCNA) in the si-lncBNIP3 group were significantly higher than in the control group. The Western Blot (WB) results also showed that the protein expression level of PCNA transfection of si-lncBNIP3 was significantly higher than in the control group. Similarly, the enrichment of lncBNIP3 significantly decreased the EdU-positive cells in the bovine preadipocytes. The results of flow cytometry and CCK8 assay also showed that overexpression of lncBNIP3 inhibited the proliferation of bovine preadipocytes. In addition, the overexpression of lncBNIP3 significantly inhibited the mRNA expressions of CCNB1 and PCNA. The WB results showed that the overexpression of lncBNIP3 significantly inhibited the expression of the CCNB1 protein level. To further explore the mechanism of lncBNIP3 on the proliferation of intramuscular preadipocytes, RNA-seq was performed after interference with si-lncBNIP3, and 660 differentially expressed genes (DEGs) were found, including 417 up-regulated DEGs and 243 down-regulated DEGs. The KEGG pathway analysis showed that the cell cycle was the most significant pathway for the functional enrichment of DEGs, followed by the DNA replication pathway. The RT-qPCR quantified the expression of twenty DEGs in the cell cycle. Therefore, we speculated that lncBNIP3 regulated intramuscular preadipocyte proliferation through the cell cycle and DNA replication pathways. To further confirm this hypothesis, the cell cycle inhibitor Ara-C was used to inhibit DNA replication of the S phase in intramuscular preadipocytes. Herein, Ara-C and si-lncBNIP3 were simultaneously added to the preadipocytes, and the CCK8, flow cytometry, and EdU assays were performed. The results showed that the si-lncBNIP3 could rescue the inhibitory effect of Ara-C in the bovine preadipocyte proliferation. In addition, lncBNIP3 could bind to the promoter of cell division control protein 6 (CDC6), and down-regulation of lncBNIP3 promoted the transcription activity and the expression of CDC6. Therefore, the inhibitory effect of lncBNIP3 on cell proliferation might be understood through the cell cycle pathway and CDC6 expression. This study provided a valuable lncRNA with functional roles in intramuscular fat accumulation and revealed new strategies for improving beef quality.

Genome-Wide Identification and Characterization of Long Non-Coding RNAs in Longissimus dorsi Skeletal Muscle of Shandong Black Cattle and Luxi Cattle

There is an increasing understanding of the possible regulatory role of long non-coding RNAs (LncRNA). Studies on livestock have mainly focused on the regulation of cell differentiation, fat synthesis, and embryonic development. However, there has been little study of skeletal muscle of domestic animals and the potential role of lncRNA. In this study, the transcriptome numbers of longissimus muscle of different beef cattle (Shandong black catle and Luxi catlle) were used to construct muscle related lncRNAs-miRNA-mRNA interaction network through bioinformatics analysis. This is helpful to clarify the molecular mechanism of bovine muscle development, and can be used to promote animal husbandry and improve animal husbandry production. According to the screening criteria of |FC|≧2 and q < 0.05, a total of 1,415 transcripts (of which 480 were LncRNAs) were differentially expressed (q < 0.05) in the different breeds. Further, we found that the most differentially expressed LncRNAs were found on chromosome 9, in which the differentially expressed LncRNAs targeted 1,164 protein coding genes (MYORG, Wnt4, PAK1, ADCY7,etc) (upstream and downstream<50 Kb). In addition, Pearson’s correlation coefficients of co-expression levels indicated a potential trans regulatory relationship between the differentially expressed LncRNAs and 43844 mRNAs (r > 0.9). The identified co-expressed mRNAs (MYORG, Dll1, EFNB2, SOX6, MYOCD, and MYLK3) are related to the formation of muscle structure, and enriched in muscle system process, strained muscle cell differentiation, muscle cell development, striated muscle tissue development, calcium signaling, and AMPK signaling. Additionally, we also found that some LncRNAs (LOC112444238, LOC101903367, LOC104975788, LOC112441863, LOC112449549, and LOC101907194) may interact with miRNAs related to cattle muscle growth and development. Based on this, we constructed a LncRNAs-miRNA-mRNA interaction network as the putative basis for biological regulation in cattle skeletal muscle. Interestingly, a candidate differential LncRNA (LOC104975788) and a protein-coding gene (Pax7) contain miR-133a binding sites and binding was confirmed by luciferase reporter assay. LOC104975788 may combined miR-133a competitively with Pax7, thus relieving the inhibitory effect of miR-133a on Pax7 to regulate skeletal muscle development. These results will provide the theoretical basis for further study of LncRNA regulation and activity in different cattle breeds.

Mammary epithelial cell transcriptome reveals potential roles of lncRNAs in regulating milk synthesis pathways in Jersey and Kashmiri cattle

Background

Long noncoding RNAs (lncRNAs) are now proven as essential regulatory elements, playing diverse roles in many biological processes including mammary gland development. However, little is known about their roles in the bovine lactation process.

Results

To identify and characterize the roles of lncRNAs in bovine lactation, high throughput RNA sequencing data from Jersey (high milk yield producer), and Kashmiri cattle (low milk yield producer) were utilized. Transcriptome data from three Kashmiri and three Jersey cattle throughout their lactation stages were utilized for differential expression analysis. At each stage (early, mid and late) three samples were taken from each breed. A total of 45 differentially expressed lncRNAs were identified between the three stages of lactation. The differentially expressed lncRNAs were found co-expressed with genes involved in the milk synthesis processes such as GPAM, LPL, and ABCG2 indicating their potential regulatory effects on milk quality genes. KEGG pathways analysis of potential cis and trans target genes of differentially expressed lncRNAs indicated that 27 and 48 pathways were significantly enriched between the three stages of lactation in Kashmiri and Jersey respectively, including mTOR signaling, PI3K-Akt signaling, and RAP1 signaling pathways. These pathways are known to play key roles in lactation biology and mammary gland development.

Conclusions

Expression profiles of lncRNAs across different lactation stages in Jersey and Kashmiri cattle provide a valuable resource for the study of the regulatory mechanisms involved in the lactation process as well as facilitate understanding of the role of lncRNAs in bovine lactation biology.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08406-x.

Transcriptome sequencing analysis for the identification of stable lncRNAs associated with bovine Staphylococcus aureus mastitis

Background

Staphylococcus aureus (S. aureus) mastitis is one of the most difficult diseases to treat in lactating dairy cows worldwide. S. aureus with different lineages leads to different host immune responses. Long non-coding RNAs (lncRNAs) are reported to be widely involved in the progress of inflammation. However, no research has identified stable lncRNAs among different S. aureus strain infections. In addition, folic acid (FA) can effectively reduce inflammation, and whether the inflammatory response caused by S. aureus can be reduced by FA remains to be explored.

Methods

lncRNA transcripts were identified from Holstein mammary gland tissues infected with different concentrations of S. aureus (in vivo) and mammary alveolar cells (Mac-T cells, in vitro) challenged with different S. aureus strains. Differentially expressed (DE) lncRNAs were evaluated, and stable DE lncRNAs were identified in vivo and in vitro. On the basis of the gene sequence conservation and function conservation across species, key lncRNAs with the function of potentially immune regulation were retained for further analysis. The function of FA on inflammation induced by S. aureus challenge was also investigated. Then, the association analysis between these keys lncRNA transcripts and hematological parameters (HPs) was carried out. Lastly, the knockdown and overexpression of the important lncRNA were performed to validate the gene function on the regulation of cell immune response.

Results

Linear regression analysis showed a significant correlation between the expression levels of lncRNA shared by mammary tissue and Mac-T cells (P < 0.001, R² = 0.3517). lncRNAs PRANCR and TNK2–AS1 could be regarded as stable markers associated with bovine S. aureus mastitis. Several HPs could be influenced by SNPs around lncRNAs PRANCR and TNK2–AS1. The results of gene function validation showed PRANCR regulates the mRNA expression of SELPLG and ITGB2 within the S. aureus infection pathway and the Mac-T cells apoptosis. In addition, FA regulated the expression change of DE lncRNA involved in toxin metabolism and inflammation to fight against S. aureus infection.

Conclusions

The remarkable association between SNPs around these two lncRNAs and partial HP indicates the potentially important role of PRANCR and TNK2–AS1 in immune regulation. Stable DE lncRNAs PRANCR and TNK2–AS1 can be regarded as potential targets for the prevention of bovine S. aureus mastitis. FA supplementation can reduce the negative effect of S. aureus challenge by regulating the expression of lncRNAs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40104-021-00639-2.

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.

The role of long noncoding RNAs in livestock adipose tissue deposition - A review

With the development of sequencing technology, numerous, long noncoding RNAs (lncRNAs) have been discovered and annotated. Increasing evidence has shown that lncRNAs play an essential role in regulating many biological and pathological processes, especially in cancer. However, there have been few studies on the roles of lncRNAs in livestock production. In animal products, meat quality and lean percentage are vital economic traits closely related to adipose tissue deposition. However, adipose tissue accumulation is also a pivotal contributor to obesity, diabetes, atherosclerosis, and many other diseases, as demonstrated by human studies. In livestock production, the mechanism by which lncRNAs regulate adipose tissue deposition is still unclear. In addition, the phenomenon that different animal species have different adipose tissue accumulation abilities is not well understood. In this review, we summarize the characteristics of lncRNAs and their four functional archetypes and review the current knowledge about lncRNA functions in adipose tissue deposition in livestock species. This review could provide theoretical significance to explore the functional mechanisms of lncRNAs in adipose tissue accumulation in animals.

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.