A catalogue of novel bovine long noncoding RNA across 18 tissues

Novel lncRNA regulatory elements in milk somatic cells of Holstein dairy cows associated with mastitis

Despite regulatory elements such as long non - coding RNAs representing most of the transcriptome, the functional understanding of long non - coding RNAs in relation to major health conditions including bovine mastitis is limited. This study examined the milk somatic cell transcriptome from udder quarters of 6 Holstein dairy cows to identify differentially expressed long non - coding RNAs using RNA - Sequencing. Ninety - four differentially expressed long non - coding RNAs are identified, 5 of which are previously annotated for gene name and length, 11 are annotated for gene name and 78 are novel, having no gene name or length previously annotated. Significant inflammatory response and regulation of immune response pathways (false discovery rate < 0.05) are associated with the differentially expressed long non - coding RNAs. QTL annotation analysis revealed 31 QTL previously annotated in the genomic regions of the 94 differentially expressed long non - coding RNAs, and the majority are associated with milk traits. This research provides a better understanding of long non - coding RNAs regulatory elements in milk somatic cells, which may enhance current breeding strategies for more adaptable or high mastitis resistant cattle.

Comprehensive transcriptomic analysis of long non-coding RNAs in bovine ovarian follicles and early embryos

Long non-coding RNAs (lncRNAs) have been the subject of numerous studies over the past decade. First thought to come from aberrant transcriptional events, lncRNAs are now considered a crucial component of the genome with roles in multiple cellular functions. However, the functional annotation and characterization of bovine lncRNAs during early development remain limited. In this comprehensive analysis, we review lncRNAs expression in bovine ovarian follicles and early embryos, based on a unique database comprising 468 microarray hybridizations from a single platform designed to target 7,724 lncRNA transcripts, of which 5,272 are intergenic (lincRNA), 958 are intronic, and 1,524 are antisense (lncNAT). Compared to translated mRNA, lncRNAs have been shown to be more tissue-specific and expressed in low copy numbers. This analysis revealed that protein-coding genes and lncRNAs are both expressed more in oocytes. Differences between the oocyte and the 2-cell embryo are also more apparent in terms of lncRNAs than mRNAs. Co-expression network analysis using WGCNA generated 25 modules with differing proportions of lncRNAs. The modules exhibiting a higher proportion of lncRNAs were found to be associated with fewer annotated mRNAs and housekeeping functions. Functional annotation of co-expressed mRNAs allowed attribution of lncRNAs to a wide array of key cellular events such as meiosis, translation initiation, immune response, and mitochondrial related functions. We thus provide evidence that lncRNAs play diverse physiological roles that are tissue-specific and associated with key cellular functions alongside mRNAs in bovine ovarian follicles and early embryos. This contributes to add lncRNAs as active molecules in the complex regulatory networks driving folliculogenesis, oogenesis and early embryogenesis all of which are necessary for reproductive success.

Investigation of lncRNA in Bos taurus Mammary Tissue during Dry and Lactation Periods

This study aims to collect RNA-Seq data from Bos taurus samples representing dry and lactating mammary tissue, identify lncRNA transcripts, and analyze findings for their features and functional annotation. This allows for connections to be drawn between lncRNA and the lactation process. RNA-Seq data from 103 samples of Bos taurus mammary tissue were gathered from publicly available databases (60 dry, 43 lactating). The samples were filtered to reveal 214 dry mammary lncRNA transcripts and 517 lactating mammary lncRNA transcripts. The lncRNAs met common lncRNA characteristics such as shorter length, fewer exons, lower expression levels, and less sequence conservation when compared to the genome. Interestingly, several lncRNAs showed sequence similarity to genes associated with strong hair keratin intermediate filaments. Human breast cancer research has associated strong hair keratin filaments with mammary tissue cellular resilience. The lncRNAs were also associated with several genes/proteins that linked to pregnancy using expression correlation and gene ontology. Such findings indicate that there are crucial relationships between the lncRNAs found in mammary tissue and the development of the tissue, to meet both the animal's needs and our own production needs; these relationships should be further investigated to ensure that we continue to breed the most resilient, efficient dairy cattle.

Karrow N, Mao Y. Identification and Classification of Long Non-Coding RNAs in the Mammary Gland of the Holstein Cow

The mammary glands, responsible for milk secretion, are regulated at a local level by various hormones, growth factors, non-coding RNAs, and other elements. Recent research has discovered the presence of lncRNAs in these glands, with suggestions that they may be essential for the maintenance and function of mammary glands. Besides directly controlling the gene and protein expression, lncRNAs are believed to play a significant part in numerous physiological and pathological processes. This study focused on examining the mammary gland tissues of Chinese Holstein cows, to identify and categorize long non-coding RNAs (lncRNAs). The research intended to distinguish lncRNAs in the mammary tissues of Holstein cows and contrast them between lactation and non-lactation periods. In this study, mammary gland tissues were sampled from three Holstein cows in early lactation (n = 3, 30 days postpartum) and non-lactation (n = 3, 315 days postpartum) on a large dairy farm in Jiangsu province. Mammary tissue samples were collected during early lactation and again during non-lactation. In total, we detected 1905 lncRNAs, with 57.3% being 500 bp and 612 intronic lncRNAs. The exon count for lncRNAs varied from 2 to 10. It was observed that 96 lncRNA expressions markedly differed between the two stages, with 83 genes being upregulated and 53 downregulated. Enrichment analysis results revealed that Gene Ontology (GO) analysis was primarily abundant in cellular processes. The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis indicated that target genes were predominantly abundant in metabolic pathways, fatty acid biosynthesis, the immune system, and glycosphingolipid biosynthesis. This study analyzed the expression profile and characteristics of lncRNAs in the mammary gland tissues of Holstein cows during both lactation and non-lactation stages, forming a foundation for further investigation into the functional roles of lncRNAs in Holstein cows throughout lactation.

Long Non-Coding RNAs and Their "Discrete" Contribution to IBD and Johne's Disease-What Stands out in the Current Picture? A Comprehensive Review

Non-coding RNAs (ncRNA) have paved the way to new perspectives on the regulation of gene expression, not only in biology and medicine, but also in associated fields and technologies, ensuring advances in diagnostic means and therapeutic modalities. Critical in this multistep approach are the associations of long non-coding RNA (lncRNA) with diseases and their causal genes in their networks of interactions, gene enrichment and expression analysis, associated pathways, the monitoring of the involved genes and their functional roles during disease progression from one stage to another. Studies have shown that Johne's Disease (JD), caused by Mycobacterium avium subspecies partuberculosis (MAP), shares common lncRNAs, clinical findings, and other molecular entities with Crohn's Disease (CD). This has been a subject of vigorous investigation owing to the zoonotic nature of this condition, although results are still inconclusive. In this review, on one hand, the current knowledge of lncRNAs in cells is presented, focusing on the pathogenesis of gastrointestinal-related pathologies and MAP-related infections and, on the other hand, we attempt to dissect the associated genes and pathways involved. Furthermore, the recently characterized and novel lncRNAs share common pathologies with IBD and JD, including the expression, molecular networks, and dataset analysis results. These are also presented in an attempt to identify potential biomarkers pertinent to cattle and human disease phenotypes.

Investigation of rumen long noncoding RNA before and after weaning in cattle

Background

This study aimed to identify long non-coding RNA (lncRNA) from the rumen tissue in dairy cattle, explore their features including expression and conservation levels, and reveal potential links between lncRNA and complex traits that may indicate important functional impacts of rumen lncRNA during the transition to the weaning period.

Results

A total of six cattle rumen samples were taken with three replicates from before and after weaning periods, respectively. Total RNAs were extracted and sequenced with lncRNA discovered based on size, coding potential, sequence homology, and known protein domains. As a result, 404 and 234 rumen lncRNAs were identified before and after weaning, respectively. However, only nine of them were shared under two conditions, with 395 lncRNAs found only in pre-weaning tissues and 225 only in post-weaning samples. Interestingly, none of the nine common lncRNAs were differentially expressed between the two weaning conditions. LncRNA averaged shorter length, lower expression, and lower conservation scores than the genome overall, which is consistent with general lncRNA characteristics. By integrating rumen lncRNA before and after weaning with large-scale GWAS results in cattle, we reported significant enrichment of both pre- and after-weaning lncRNA with traits of economic importance including production, reproduction, health, and body conformation phenotypes.

Conclusions

The majority of rumen lncRNAs are uniquely expressed in one of the two weaning conditions, indicating a functional role of lncRNA in rumen development and transition of weaning. Notably, both pre- and post-weaning lncRNA showed significant enrichment with a variety of complex traits in dairy cattle, suggesting the importance of rumen lncRNA for cattle performance in the adult stage. These relationships should be further investigated to better understand the specific roles lncRNAs are playing in rumen development and cow performance.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08758-4.

Emerging Roles of Noncoding RNAs in Bovine Mastitis Diseases

Non-coding RNAs (ncRNAs) are an abundant class of RNA with varying nucleotide lengths. They have been shown to have great potential in eutherians/human disease diagnosis and treatments and are now gaining more importance for the improvement of diseases in livestock. To date, thousands of ncRNAs have been discovered in the bovine genome and the continuous advancement in deep sequencing technologies and various bioinformatics tools has enabled the elucidation of their roles in bovine health. Among farm animals' diseases, mastitis, a common inflammatory disease in cattle, has caused devastating economic losses to dairy farmers over the last few decades. Here, we summarize the biology of bovine mastitis and comprehensively discuss the roles of ncRNAs in different types of mastitis infection. Based on our findings and relevant literature, we highlighted various evidence of ncRNA roles in mastitis. Different approaches (in vivo versus in vitro) for exploring ncRNA roles in mastitis are emphasized. More particularly, the potential applications of emerging genome editing technologies, as well as integrated omics platforms for ncRNA studies and implications for mastitis are presented.

Effect of Diet Enriched With Hemp Seeds on Goat Milk Fatty Acids, Transcriptome, and miRNAs

In dairy ruminants, a diet supplemented with feed rich in unsaturated fatty acids can be an effective medium to increase the health-promoting properties of milk, although their effect on the pathways/genes involved in these processes has not been properly and completely defined to date. To improve our knowledge of the cell’s activity in specific conditions, next-generation RNA-sequencing technology was used to allow whole transcriptome characterization under given conditions. In addition to this, microRNAs (miRNAs) have recently been known as post-transcriptional regulators in fatty acid and cholesterol metabolism by targeting lipid metabolism genes. In this study, to analyze the transcriptome and miRNAs in goat milk after a supplemental diet enriched with linoleic acid (hemp seeds), next-generation RNA-sequencing was used in order to point out the general biological mechanisms underlying the effects related to milk fat metabolism. Ten pluriparous Alpine goats were fed with the same pretreatment diet for 40 days; then, they were arranged to two dietary treatments consisting of control (C) and hemp seed (H)-supplemented diets. Milk samples were collected at 40 (time point = T0) and 140 days of lactation (time point = T1). Milk fatty acid (FA) profiles revealed a significant effect of hemp seeds that determined a strong increment in the preformed FA, causing a reduction in the concentration of de-novo FA. Monounsaturated and polyunsaturated n−3 FAs were increased by hemp treatment, determining a reduction in the n−6/n−3 ratio. After removing milk fats and proteins, RNA was extracted from the milk cells and transcriptomic analysis was conducted using Illumina RNA-sequencing. A total of 3,835 genes were highly differentially expressed (p-value &lt; 0.05, fold change &gt; 1.5, and FDR &lt; 0.05) in the H group. Functional analyses evidenced changes in metabolism, immune, and inflammatory responses. Furthermore, modifications in feeding strategies affected also key transcription factors regulating the expression of several genes involved in milk fat metabolism, such as peroxisome proliferator-activated receptors (PPARs). Moreover, 38 (15 known and 23 novel) differentially expressed miRNAs were uncovered in the H group and their potential functions were also predicted. This study gives the possibility to improve our knowledge of the molecular changes occurring after a hemp seed supplementation in the goat diet and increase our understanding of the relationship between nutrient variation and phenotypic effects.

Transcriptome Sequencing to Identify Important Genes and lncRNAs Regulating Abdominal Fat Deposition in Ducks

Simple Summary

Abdominal fat deposition affects the quality of duck meat and the feed conversion ratio. Here, we performed transcriptome sequencing of the abdominal adipose tissue of ducks with high and low abdominal fat rate by RNA sequencing, exploring the key regulatory genes and lncRNAs related to abdominal fat deposition. As a result, several candidate genes, lncRNAs, and pathways related to abdominal fat deposition in ducks were retrieved. This study lays the foundations for exploring molecular mechanisms underlying the regulation of abdominal fat deposition in ducks, providing a theoretical reference for breeding high-quality meat-producing ducks.

Abstract

Abdominal fat deposition is an important trait in meat-producing ducks. F2 generations of 304 Cherry Valley and Runzhou Crested White ducks were studied to identify genes and lncRNAs affecting abdominal fat deposition. RNA sequencing was used to study abdominal fat tissue of four ducks each with high or low abdominal fat rates. In all, 336 upregulated and 297 downregulated mRNAs, and 95 upregulated and 119 downregulated lncRNAs were identified. Target gene prediction of differentially expressed lncRNAs identified 602 genes that were further subjected to Gene Ontology and KEGG pathway analysis. The target genes were enriched in pathways associated with fat synthesis and metabolism and participated in biological processes, including Linoleic acid metabolism, lipid storage, and fat cell differentiation, indicating that these lncRNAs play an important role in abdominal fat deposition. PPAPA, FOXO3, FASN, PNPLA2, FKBP5, TCF7L2, BMP2, FGF2, LIFR, ZBTB16, SIRT, GYG2, NCOR1, and NR3C1 were involved in the regulation of abdominal fat deposition. PNPLA2, TCF7L2, FGF2, LIFR, BMP2, FKBP5, GYG2, and ZBTB16 were regulated by the lncRNAs TCONS_00038080, TCONS_0033547, TCONS_00066773, XR_001190174.3, XR_003492471.1, XR_003493494.1, XR_001192142.3, XR_002405656.2, XR_002401822.2, XR_003497063.1, and so on. This study lays foundations for exploring molecular mechanisms underlying the regulation of abdominal fat deposition in ducks and provides a theoretical basis for breeding high-quality meat-producing ducks.

Identification of Long Noncoding RNAs Involved in Eyelid Pigmentation of Hereford Cattle

Several ocular pathologies in cattle, such as ocular squamous cell carcinoma and infectious keratoconjunctivitis, have been associated with low pigmentation of the eyelids. The main objective of this study was to analyze the transcriptome of eyelid skin in Hereford cattle using strand-specific RNA sequencing technology to characterize and identify long noncoding RNAs (lncRNAs). We compared the expression of lncRNAs between pigmented and unpigmented eyelids and analyzed the interaction of lncRNAs and putative target genes to reveal the genetic basis underlying eyelid pigmentation in cattle. We predicted 4,937 putative lncRNAs mapped to the bovine reference genome, enriching the catalog of lncRNAs in Bos taurus. We found 27 differentially expressed lncRNAs between pigmented and unpigmented eyelids, suggesting their involvement in eyelid pigmentation. In addition, we revealed potential links between some significant differentially expressed lncRNAs and target mRNAs involved in the immune response and pigmentation. Overall, this study expands the catalog of lncRNAs in cattle and contributes to a better understanding of the biology of eyelid pigmentation.

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.

Differential expression of long non-coding RNAs under Peste des petits ruminants virus (PPRV) infection in goats

Peste des petits ruminants (PPR) characterized by fever, sore mouth, conjunctivitis, gastroenteritis, and pneumonia, is an acute, highly contagious viral disease of sheep and goats. The role of long non-coding RNAs (lncRNAs) in PPRV infection has not been explored to date. In this study, the transcriptome profiles of virulent Peste des petits ruminants virus (PPRV) infected goat tissues – lung and spleen were analyzed to identify the role of lncRNAs in PPRV infection. A total of 13,928 lncRNA transcripts were identified, out of which 170 were known lncRNAs. Intergenic lncRNAs (7625) formed the major chunk of the novel lncRNA transcripts. Differential expression analysis revealed that 15 lncRNAs (11 downregulated and 4 upregulated) in the PPRV infected spleen samples and 16 lncRNAs (13 downregulated and 3 upregulated) in PPRV infected lung samples were differentially expressed as compared to control. The differentially expressed lncRNAs (DElncRNAs) possibly regulate various immunological processes related to natural killer cell activation, antigen processing and presentation, and B cell activity, by regulating the expression of mRNAs through the cis- or trans-regulatory mechanism. Functional enrichment analysis of differentially expressed mRNAs (DEmRNAs) revealed enrichment of immune pathways and biological processes in concordance with the pathways in which correlated lncRNA-neighboring genes were enriched. The results suggest that a coordinated immune response is raised in both lung and spleen tissues of the goat through mRNA-lncRNA crosstalk.

Expression of lncRNAs in response to bacterial infections of goat mammary epithelial cells reveals insights into mammary gland diseases

Mastitis or inflammation of the mammary gland is a highly economic and deadly alarming disease for the dairy sector as well as policymakers caused by microbial infection. Transcriptomic and proteomic approaches have been widely employed to identify the underlying molecular mechanisms of bacterial infections in the mammary gland. Numerous differentially expressed mRNAs, miRNAs, and proteins together with their associated signaling pathways have been identified during bacterial infection, paving the way for analysis of their biological functions. Long noncoding RNAs (lncRNAs) are important regulators of multiple biological processes. However, little is known regarding their role in bacterial infection in mammary epithelial cells. Hence, RNA-sequencing was performed by infecting primary mammary epithelial cells (pMECs) with both gram-negative (E. coli) and gram-positive bacteria (S. aureus). Using stringent pipeline, a set of 1957 known and 1175 novel lncRNAs were identified, among which, 112 lncRNAs were found differentially expressed in bacteria challenged PMECs compared with the control. Additionally, potential targets of the lncRNAs were predicted in cis- and trans-configuration. KEGG analysis revealed that DE lncRNAs were associated with at least 15 immune-related pathways. Therefore, our study revealed that bacterial challenge triggers the expression of lncRNAs associated with immune response and defense mechanisms in goat mammary epithelial cells.

Identification of Potential Long Non-Coding RNA Candidates that Contribute to Triple-Negative Breast Cancer in Humans through Computational Approach

Breast cancer (BC) is the most frequent malignancy identified in adult females, resulting in enormous financial losses worldwide. Owing to the heterogeneity as well as various molecular subtypes, the molecular pathways underlying carcinogenesis in various forms of BC are distinct. Therefore, the advancement of alternative therapy is required to combat the ailment. Recent analyses propose that long non-coding RNAs (lncRNAs) perform an essential function in controlling immune response, and therefore, may provide essential information about the disorder. However, their function in patients with triple-negative BC (TNBC) has not been explored in detail. Here, we analyzed the changes in the genomic expression of messenger RNA (mRNA) and lncRNA in standard control in response to cancer metastasis using publicly available single-cell RNA-Seq data. We identified a total of 197 potentially novel lncRNAs in TNBC patients of which 86 were differentially upregulated and 111 were differentially downregulated. In addition, among the 909 candidate lncRNA transcripts, 19 were significantly differentially expressed (DE) of which three were upregulated and 16 were downregulated. On the other hand, 1901 mRNA transcripts were significantly DE of which 1110 were upregulated and 791 were downregulated by TNBCs subtypes. The Gene Ontology (GO) analyses showed that some of the host genes were enriched in various biological, molecular, and cellular functions. The Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis showed that some of the genes were involved in only one pathway of prostate cancer. The lncRNA-miRNA-gene network analysis showed that the lncRNAs TCONS_00076394 and TCONS_00051377 interacted with breast cancer-related micro RNAs (miRNAs) and the host genes of these lncRNAs were also functionally related to breast cancer. Thus, this study provides novel lncRNAs as potential biomarkers for the therapeutic intervention of this cancer subtype.

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.

PacBio Single-Molecule Long-Read Sequencing Provides New Light on the Complexity of Full-Length Transcripts in Cattle

Cattle (Bos taurus) is one of the most widely distributed livestock species in the world, and provides us with high-quality milk and meat which have a huge impact on the quality of human life. Therefore, accurate and complete transcriptome and genome annotation are of great value to the research of cattle breeding. In this study, we used error-corrected PacBio single-molecule real-time (SMRT) data to perform whole-transcriptome profiling in cattle. Then, 22.5 Gb of subreads was generated, including 381,423 circular consensus sequences (CCSs), among which 276,295 full-length non-chimeric (FLNC) sequences were identified. After correction by Illumina short reads, we obtained 22,353 error-corrected isoforms. A total of 305 alternative splicing (AS) events and 3,795 alternative polyadenylation (APA) sites were detected by transcriptome structural analysis. Furthermore, we identified 457 novel genes, 120 putative transcription factors (TFs), and 569 novel long non-coding RNAs (lncRNAs). Taken together, this research improves our understanding and provides new insights into the complexity of full-length transcripts in cattle.

Emerging Role of MicroRNA-200 Family in Dentistry

MicroRNAs (miRNAs) are endogenous non-coding RNAs ~22 nucleotides in length, which have been shown to participate in various biological processes. As one of the most researched miRNAs, the miR-200 family has been found to regulate several factors that are associated with the epithelial to mesenchymal transition (EMT) and cancer stem cells (CSCs) behavior. In this review, we briefly summarize the background of the miR-200 family and their implication in various dental diseases. We focus on the expression changes, biological functions, and clinical significance of the miR-200 family in oral cancer; periodontitis; oral potentially malignant disorder; gingival overgrowth; and other periodontal diseases. Additionally, we discuss the use of the miR-200 family as molecular biomarkers for diagnosis, prognostic, and therapeutic application.

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.

Identification of Long Non-coding RNA Isolated From Naturally Infected Macrophages and Associated With Bovine Johne's Disease in Canadian Holstein Using a Combination of Neural Networks and Logistic Regression

Mycobacterium avium ssp. paratuberculosis (MAP) causes chronic enteritis in most ruminants. The pathogen MAP causes Johne's disease (JD), a chronic, incurable, wasting disease. Weight loss, diarrhea, and a gradual drop in milk production characterize the disease's clinical phase, culminating in death. Several studies have characterized long non-coding RNA (lncRNA) in bovine tissues, and a previous study characterizes (lncRNA) in macrophages infected with MAP in vitro. In this study, we aim to characterize the lncRNA in macrophages from cows naturally infected with MAP. From 15 herds, feces and blood samples were collected for each cow older than 24 months, twice yearly over 3–5 years. Paired samples were analyzed by fecal PCR and blood ELISA. We used RNA-seq data to study lncRNA in macrophages from 33 JD(+) and 33 JD(–) dairy cows. We performed RNA-seq analysis using the “new Tuxedo” suite. We characterized lncRNA using logistic regression and multilayered neural networks and used DESeq2 for differential expression analysis and Panther and Reactome classification systems for gene ontology (GO) analysis. The study identified 13,301 lncRNA, 605 of which were novel lncRNA. We found seven genes close to differentially expressed lncRNA, including CCDC174, ERI1, FZD1, TWSG1, ZBTB38, ZNF814, and ZSCAN4. None of the genes associated with susceptibility to JD have been cited in the literature. LncRNA target genes were significantly enriched for biological process GO terms involved in immunity and nucleic acid regulation. These include the MyD88 pathway (TLR5), GO:0043312 (neutrophil degranulation), GO:0002446 (neutrophil-mediated immunity), and GO:0042119 (neutrophil activation). These results identified lncRNA with potential roles in host immunity and potential candidate genes and pathways through which lncRNA might function in response to MAP infection.

Identification and co-expression analysis of long noncoding RNAs and mRNAs involved in the deposition of intramuscular fat in Aohan fine-wool sheep

BACKGROUND

Intramuscular fat (IMF) content has become one of the most important indicators for measuring meat quality, and levels of IMF are affected by various genes. Long non-coding RNAs (lncRNAs) are widely expressed non-coding RNAs that play an important regulatory role in a variety of biological processes; however, research on the lncRNAs involved in sheep IMF deposition is still in its infancy. Aohan fine-wool sheep (AFWS), one of China's most important meat-hair, dual-purpose sheep breed, provides a great model for studying the role of lncRNAs in the regulation of IMF deposition. We identified lncRNAs by RNA sequencing in Longissimus thoracis et lumborum (LTL) samples of sheep at two ages: 2 months (Mth-2) and 12 months (Mth-12).

RESULTS

We identified a total of 26,247 genes and 6935 novel lncRNAs in LTL samples of sheep. Among these, 199 mRNAs and 61 lncRNAs were differentially expressed. We then compared the structural characteristics of lncRNAs and mRNAs. We obtained target genes of differentially expressed lncRNAs (DELs) and performed enrichment analyses using Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG). We found that target mRNAs were enriched in metabolic processes and developmental pathways. One pathway was significantly enriched, namely tight junction. Based on the analysis of critical target genes, we obtained seven candidate lncRNAs that potentially regulated lipid deposition and constructed a lncRNA-mRNA co-expression network that included MSTRG.4051.3-FZD4, MSTRG.16157.3-ULK1, MSTRG.21053.3-PAQR3, MSTRG.19941.2-TPI1, MSTRG.12864.1-FHL1, MSTRG.2469.2-EXOC6 and MSTRG.21381.1-NCOA1. We speculated that these candidate lncRNAs might play a role by regulating the expression of target genes. We randomly selected five mRNAs and five lncRNAs to verify the accuracy of the sequencing data by qRT-PCR.

CONCLUSIONS

Our study identified the differentially expressed mRNAs and lncRNAs during intramuscular lipid deposition in Aohan fine-wool sheep. The work may widen the knowledge about the annotation of the sheep genome and provide a working basis for investigating intramuscular fat deposition in sheep.

Differential expression and prediction of function of lncRNAs in the ovaries of low and high fecundity Hanper sheep

Litter size is an important trait that determines the production efficiency of sheep bred for meat. Its detailed investigation can reveal the molecular mechanisms that control the fecundity of sheep and possibly accelerate the breeding process of new varieties of sheep that have high prolificacy. Long non-coding RNAs (lncRNAs) have proven to be an important factor in the regulation of follicular development. However, the mechanisms by which lncRNAs regulate litter size in sheep remain unclear. In the present study, ovarian tissues from the follicular (F) or luteal phase (L) of Hanper sheep that were either monotocous (M) or polytocous (P; FM, FP, LM and LP groups) were collected and sequenced to identify differentially expressed lncRNAs and predict their function. The results indicate that the number of up- and down-regulated lncRNAs in the follicular phase (FM vs. FP) was 95 and 111 and 109 and 49, respectively, in the luteal phase (LM vs. LP). The functional enrichment of the different lncRNAs coexpressed with mRNA was analysed. The results demonstrated that the KISS1-GnRH-LH/FSH-E2 and EGF-EGFR-RAS-PI3K signalling pathways promoted the initiation of the primordial period, follicular development and ovulation in the follicular phase (FM vs. FP). During the luteal phase (LM vs. LP), the production and development of the corpus luteum in ewes was influenced by the KITLG-KIT/FGF-FGFR/HGF-MET-RAS-ERK signalling pathway. STEM clustering functional enrichment analysis of the differentially expressed lncRNAs indicated that profile11 was principally enriched in the Cytokine-Jak-STAT, PDGF-PDGFR-PI3K and KITLG-KIT-RAS-ERK signalling pathways. By analysis of the differential expression of the lncRNAs and their expression in each group, lncRNAs Xist (loc101112291) and Gtl2 (loc101123329) were found to be highly expressed, suggesting that regulation of follicular development was mediated through methylation processes.

Dual Effects of Non-Coding RNAs (ncRNAs) in Cancer Stem Cell Biology

The identification of cancer stem cells (CSCs) as initiators of carcinogenesis has revolutionized the era of cancer research and our perception for the disease treatment options. Additional CSC features, including self-renewal and migratory and invasive capabilities, have further justified these cells as putative diagnostic, prognostic, and therapeutic targets. Given the CSC plasticity, the identification of CSC-related biomarkers has been a serious burden in CSC characterization and therapeutic targeting. Over the past decades, a compelling amount of evidence has demonstrated critical regulatory functions of non-coding RNAs (ncRNAs) on the exclusive features of CSCs. We now know that ncRNAs may interfere with signaling pathways, vital for CSC phenotype maintenance, such as Notch, Wnt, and Hedgehog. Here, we discuss the multifaceted contribution of microRNAs (miRNAs), long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs), as representative ncRNA classes, in sustaining the CSC-like traits, as well as the underlying molecular mechanisms of their action in various CSC types. We further discuss the use of CSC-related ncRNAs as putative biomarkers of high diagnostic, prognostic, and therapeutic value.

Exploring the Regulatory Potential of Long Non-Coding RNA in Feed Efficiency of Indicine Cattle

Long non-coding RNA (lncRNA) can regulate several aspects of gene expression, being associated with complex phenotypes in humans and livestock species. In taurine beef cattle, recent evidence points to the involvement of lncRNA in feed efficiency (FE), a proxy for increased productivity and sustainability. Here, we hypothesized specific regulatory roles of lncRNA in FE of indicine cattle. Using RNA-Seq data from the liver, muscle, hypothalamus, pituitary gland and adrenal gland from Nellore bulls with divergent FE, we submitted new transcripts to a series of filters to confidently predict lncRNA. Then, we identified lncRNA that were differentially expressed (DE) and/or key regulators of FE. Finally, we explored lncRNA genomic location and interactions with miRNA and mRNA to infer potential function. We were able to identify 126 relevant lncRNA for FE in Bos indicus, some with high homology to previously identified lncRNA in Bos taurus and some possible specific regulators of FE in indicine cattle. Moreover, lncRNA identified here were linked to previously described mechanisms related to FE in hypothalamus-pituitary-adrenal axis and are expected to help elucidate this complex phenotype. This study contributes to expanding the catalogue of lncRNA, particularly in indicine cattle, and identifies candidates for further studies in animal selection and management.

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.

Genome-Wide Analysis Reveals Changes in Polled Yak Long Non-coding RNAs in Skeletal Muscle Development

Long non-coding RNAs (lncRNAs) have been extensively studied in recent years. Numerous lncRNAs have been identified in mice, rats, and humans, some of which play important roles in muscle formation and development. However, little is known about lncRNA regulators that affect muscle development in yak (Bos grunniens). LncRNA expression during skeletal muscle development in yak was analyzed by RNA sequencing at three development stages: 3 years (group A), 6 months (group M), and 90-day-old fetuses (group E). A total of 1180 lncRNAs were identified in the three development stages. Compared with group E, 154 were upregulated and 130 were downregulated in group A. Compared with group A, 31 were upregulated and 29 were downregulated in group M. Compared with group E, 147 were upregulated and 149 were downregulated in group M (padj < 0.001, |log2FC| > 1.2). In addition, functional annotation analysis based on gene ontology (GO) and the Kyoto protocol encyclopedia of genes and genomes (KEGG) database showed that differentially expressed lncRNAs (DElncRNAs) were cis–trans target genes. The results showed that DElncRNAs were mainly involved in PI3K-Akt signaling pathway, focal adhesion, MAPK signaling pathway, apoptosis, and p53 signaling pathway. Furthermore, RTL1, IGF2, MEF2C, Pax7, and other well-known muscle development regulators were included in a co-expression network of differentially expressed target genes and lncRNAs. These data will help to further clarify the function of lncRNAs in the different stages of skeletal muscle developmental in yak.

Systematic identification of long intergenic non-coding RNAs expressed in bovine oocytes

BACKGROUND

Long non-coding RNAs (lncRNAs) are key regulators of diverse cellular processes. Although a number of studies have reported the identification of bovine lncRNAs across many tissues, very little is known about the identity and characteristics of lncRNAs in bovine oocytes.

METHODS

A bovine oocyte cDNA library was constructed and sequenced using the Illumina HiSeq 2000 sequencing system. The oocyte transcriptome was constructed using the ab initio assembly software Scripture and Cufflinks. The assembled transcripts were categorized to identify the novel intergenic transcripts, and the coding potential of these novel transcripts was assessed using CPAT and PhyloCSF. The resulting candidate long intergenic non-coding RNAs (lincRNAs) transcripts were further evaluated to determine if any of them contain any known protein coding domains in the Pfam database. RT-PCR was used to analyze the expression of oocyte-expressed lincRNAs in various bovine tissues.

RESULTS

A total of 85 million raw reads were generated from sequencing of the bovine oocyte library. Transcriptome reconstruction resulted in the assembly of a total of 42,396 transcripts from 37,678 genomic loci. Analysis of the assembled transcripts using the step-wide pipeline resulted in the identification of 1535 oocyte lincRNAs corresponding to 1183 putative non-coding genes. A comparison of the oocyte lincRNAs with the lncRNAs reported in other bovine tissues indicated that 970 of the 1535 oocyte lincRNAs appear to be unique to bovine oocytes. RT-PCR analysis of 5 selected lincRNAs showed either specific or predominant expression of 4 lincRNAs in the fetal ovary. Functional prediction of the oocyte-expressed lincRNAs suggested their involvement in oogenesis through regulating their neighboring protein-coding genes.

CONCLUSIONS

This study provides a starting point for future research aimed at understanding the roles of lncRNAs in controlling oocyte development and early embryogenesis in cattle.

Analysis of Long Non-Coding RNAs and mRNAs Associated with Lactation in the Crop of Pigeons (Columba livia)

Pigeons have the ability to produce milk and feed their squabs. The genetic mechanisms underlying milk production in the crops of 'lactating' pigeons are not fully understood. In this study, RNA sequencing was employed to profile the transcriptome of lncRNA and mRNA in lactating and non-'lactating' pigeon crops. We identified 7066 known and 17,085 novel lncRNAs. Of these lncRNAs, 6166 were differentially expressed. Among the 15,138 mRNAs detected, 6483 were differentially expressed, including many predominant genes with known functions in the milk production of mammals. A GO annotation analysis revealed that these genes were significantly enriched in 55, 65, and 30 pathways of biological processes, cellular components, and molecular functions, respectively. A KEGG pathway enrichment analysis revealed that 12 pathways (involving 544 genes), including the biosynthesis of amino acids, the propanoate metabolism, the carbon metabolism and the cell cycle, were significantly enriched. The results provide fundamental evidence for the better understanding of lncRNAs' and differentially expressed genes' (DEGs) regulatory role in the molecular pathways governing milk production in pigeon crops. To our knowledge, this is the first genome-wide investigation of the lncRNAs in pigeon crop associated with milk production. This study provided valuable resources for differentially expressed lncRNAs and mRNAs, improving our understanding of the molecular mechanism of pigeon milk production.

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.

Genome-wide identification and characterization of novel long non-coding RNA in Ruminal tissue affected with sub-acute Ruminal acidosis from Holstein cattle

Sub-acute ruminal acidosis is a type of metabolic disorder in which affected cattle show a considerable depression of rumen pH. This leads to a dramatic decline in productivity and consequent loss of income for many dairy farms. The objective of the present study is to identify and characterize novel long non-coding RNAs (lncRNAs) in Holstein cattle affected by sub-acute ruminal acidosis. Two replicates from six animals were sequenced that bioinformatically analyzed. Results showed 6679 novel lncRNAs among which 12 intergenic lncRNAs showed differential expression (p value ≤0.05). GO and KEGG analysis revealed that calcium signaling and G protein couple-receptor pathways may be involved in regulating metabolic processes during sub-acute ruminal acidosis. Furthermore, other biological processes including transmembrane transport, adult behavior, neuroactive ligand-receptor interaction, GABAergic synapse, cholinergic synapse were significantly enriched. The present data suggest that these differentially expressed lncRNAs may play regulatory roles in modulating biological processes associated with sub-acute ruminal acidosis in cattle rumen.

MIR221HG Is a Novel Long Noncoding RNA that Inhibits Bovine Adipocyte Differentiation

Adipogenesis is a complicated but precisely orchestrated process mediated by a series of transcription factors. Our previous study has identified a novel long noncoding RNA (lncRNA) that was differentially expressed during bovine adipocyte differentiation. Because this lncRNA overlaps with miR-221 in the genome, it was named miR-221 host gene (MIR221HG). The purpose of this study was to clone the full length of MIR221HG, detect its subcellular localization, and determine the effects of MIR221HG on bovine adipocyte differentiation. The 5′ rapid amplification of cDNA ends (RACE) and 3′ RACE analyses demonstrated that MIR221HG is a transcript of 1064 nucleotides, is located on the bovine X chromosome, and contains a single exon. Bioinformatics analyses suggested that MIR221HG is an lncRNA and the promoter of MIR221HG includes the binding consensus sequences of the forkhead box C1 (FOXC1) and krüppel-like factor5 (KLF5). The semi-quantitative PCR and quantitative real-time PCR (qRT-PCR) of nuclear and cytoplasmic fractions revealed that MIR221HG mainly resides in the nucleus. Inhibition of MIR221HG significantly increased adipocyte differentiation, as indicated by a dramatic increment in the number of mature adipocytes and in the expression of the respective adipogenic markers, peroxisome proliferator-activated receptor γ (PPARγ), CCAAT/enhancer-binding protein α (C/EBPα), and fatty acid binding protein 4 (FABP4). Our results provide a basis for elucidating the mechanism by which MIR221HG regulates adipocyte differentiation.

Transcriptome profiling of the liver among the prenatal and postnatal stages in chickens

The liver is an important organ that has pivotal functions in the synthesis of several vital proteins, the metabolism of various biologically useful materials, the detoxification of toxic substances, and immune defense. Most liver functions are not mature at a young age and many changes happen during postnatal liver development, which lead to differential functions of the liver at different developmental stages. However, the transcriptome details of what changes occur in the liver after birth and the molecular mechanisms for the regulation of the developmental process are not clearly known in chickens. Here, we used RNA-sequencing to analyze the transcriptome of chicken liver from the prenatal (at an embryonic day of 13) to the postnatal stages (at 5 wk and 42 wk of age). A total of approximately 161.17 Gb of raw data were obtained, with 4,127 putative and 539 differentially expressed lncRNAs, and with 13,949 putative and 6,370 differentially expressed mRNAs. Coexpression of lncRNAs-mRNAs in hepatic transcriptome analysis showed that the liver plays important roles in providing energy for organisms through the mitochondrial respiratory chain in chickens, meanwhile, acting as a crucial part of antioxidant stress. The developmental transcriptome date revealed that antioxidant defenses are likely to act on chicken embryo development and that significant functional changes during postnatal liver development are associated with the liver maturation of chickens. These results provide a timeline for the functional transcriptome transition from the prenatal to adult stages in chickens and will be helpful to reveal the underlying molecular mechanisms of liver development.

Genome-wide association analysis reveals QTL and candidate mutations involved in white spotting in cattle

BACKGROUND

White spotting of the coat is a characteristic trait of various domestic species including cattle and other mammals. It is a hallmark of Holstein-Friesian cattle, and several previous studies have detected genetic loci with major effects for white spotting in animals with Holstein-Friesian ancestry. Here, our aim was to better understand the underlying genetic and molecular mechanisms of white spotting, by conducting the largest mapping study for this trait in cattle, to date.

RESULTS

Using imputed whole-genome sequence data, we conducted a genome-wide association analysis in 2973 mixed-breed cows and bulls. Highly significant quantitative trait loci (QTL) were found on chromosomes 6 and 22, highlighting the well-established coat color genes KIT and MITF as likely responsible for these effects. These results are in broad agreement with previous studies, although we also report a third significant QTL on chromosome 2 that appears to be novel. This signal maps immediately adjacent to the PAX3 gene, which encodes a known transcription factor that controls MITF expression and is the causal locus for white spotting in horses. More detailed examination of these loci revealed a candidate causal mutation in PAX3 (p.Thr424Met), and another candidate mutation (rs209784468) within a conserved element in intron 2 of MITF transcripts expressed in the skin. These analyses also revealed a mechanistic ambiguity at the chromosome 6 locus, where highly dispersed association signals suggested multiple or multiallelic QTL involving KIT and/or other genes in this region.

CONCLUSIONS

Our findings extend those of previous studies that reported KIT as a likely causal gene for white spotting, and report novel associations between candidate causal mutations in both the MITF and PAX3 genes. The sizes of the effects of these QTL are substantial, and could be used to select animals with darker, or conversely whiter, coats depending on the desired characteristics.

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.

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 Gene Expression Atlas of the Domestic Water Buffalo (Bubalus bubalis)

The domestic water buffalo (Bubalus bubalis) makes a major contribution to the global agricultural economy in the form of milk, meat, hides, and draught power. The global water buffalo population is predominantly found in Asia, and per head of population more people depend upon the buffalo than on any other livestock species. Despite its agricultural importance, there are comparatively fewer genomic and transcriptomic resources available for buffalo than for other livestock species. We have generated a large-scale gene expression atlas covering multiple tissue and cell types from all major organ systems collected from three breeds of riverine water buffalo (Mediterranean, Pandharpuri and Bhadawari) and used the network analysis tool Graphia Professional to identify clusters of genes with similar expression profiles. Alongside similar data, we and others have generated for ruminants as part of the Functional Annotation of Animal Genomes Consortium; this comprehensive transcriptome supports functional annotation and comparative analysis of the water buffalo genome.

Analysis of Long Non-Coding RNA and mRNA Expression Profiling in Immature and Mature Bovine (Bos taurus) Testes

Testis development and spermatogenesis are strictly regulated by numbers of genes and non-coding genes. However, long non-coding RNAs (lncRNAs) as key regulators in multitudinous biological processes have not been systematically identified in bovine testes during sexual maturation. In this study, we comprehensively analyzed lncRNA and mRNA expression profiling of six bovine testes at 3 days after birth and 13 months by RNA sequencing. 23,735 lncRNAs and 22,118 mRNAs were identified, in which 540 lncRNAs (P-value < 0.05) and 3,525 mRNAs (P-adjust < 0.05) were significantly differentially expressed (DE) between two stages. Correspondingly, the results of RT-qPCR analysis showed well correlation with the transcriptome data. Moreover, GO and KEGG enrichment analyses showed that DE genes and target genes of DE lncRNAs were enriched in spermatogenesis. Furthermore, we constructed lncRNA–gene interaction networks; consequently, 15 DE lncRNAs and 12 cis-target genes were involved. The target genes (SPATA16, TCF21, ZPBP, PACRG, ATP8B3, COMP, ACE, and OSBP2) were found associated with bovine sexual maturation. In addition, the expression of lncRNAs and cis-target genes was detected in bovine Leydig cells, Sertoli cells, and spermatogonia. Our study identified and analyzed lncRNAs and mRNAs in testis tissues, suggesting that lncRNAs may regulate testis development and spermatogenesis. Our findings provided new insights for further investigation of biological function in bovine lncRNA.

DNA methylation and functional characterization of the XIST gene during in vitro early embryo development in cattle

XIST, in association with the shorter ncRNA RepA, are essential for the initiation of X chromosome inactivation (XCI) in mice. The molecular mechanisms controlling XIST and RepA expression are well characterized in that specie. However, little is known in livestock. We aimed to characterize the DNA methylation status along the 5' portion of XIST and to characterize its transcriptional profile during early development in cattle. Three genomic regions of XIST named here as promoter, RepA and DMR1 had their DNA methylation status characterized in gametes and embryos. Expression profile of XIST was evaluated, including sense and antisense transcription. Oocytes showed higher levels of methylation than spermatozoa that was demethylated. DMR1 was hypermethylated throughout oogenesis. At the 8-16-cell embryo stage DMR1 was completed demethylated. Interestingly, RepA gain methylation during oocyte maturation and was demethylated at the blastocyst stage, later than DMR1. These results suggest that DMR1 and RepA are transient differentially methylated regions in cattle. XIST RNA was detected in matured oocytes and in single cells from the 2-cell to the morula stage, confirming the presence of maternal and embryonic transcripts. Sense and antisense transcripts were detected along the XIST in blastocyst. In silico analysis identified 63 novel transcript candidates at bovine XIST locus from both the plus and minus strands. Taking together these results improve our understanding of the molecular mechanisms involved in XCI initiation in cattle. This information may be useful for the improvement of assisted reproductive technologies in livestock considering that in vitro conditions may impair epigenetic reprogramming.

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

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

The genetic background of clinical mastitis in Holstein-Friesian cattle

Mastitis is an inflammatory disease of the mammary gland, which has a significant economic impact and is an animal welfare concern. This work examined the association between single nucleotide polymorphisms (SNPs) and copy number variations (CNVs) with the incidence of clinical mastitis (CM). Using information from 16 half-sib pairs of Holstein-Friesian cows (32 animals in total) we searched for genomic regions that differed between a healthy (no incidence of CM) and a mastitis-prone (multiple incidences of CM) half-sib. Three cows with average sequence depth of coverage below 10 were excluded, which left 13 half-sib pairs available for comparisons. In total, 191 CNV regions were identified, which were deleted in a mastitis-prone cow, but present in its healthy half-sib and overlapped in at least nine half-sib pairs. These regions overlapped with exons of 46 genes, among which APP (BTA1), FOXL2 (BTA1), SSFA2 (BTA2), OTUD3 (BTA2), ADORA2A (BTA17), TXNRD2 (BTA17) and NDUFS6 (BTA20) have been reported to influence CM. Moreover, two duplicated CNV regions present in nine healthy individuals and absent in their mastitis-affected half-sibs overlapped with exons of a cholinergic receptor nicotinic α 10 subunit on BTA15 and a novel gene (ENSBTAG00000008519) on BTA27. One CNV region deleted in nine mastitis-affected sibs overlapped with two neighbouring long non-coding RNA sequences located on BTA12. Single nucleotide polymorphisms with differential genotypes between a healthy and a mastitis-affected sib included 17 polymorphisms with alternate alleles in eight affected and healthy half-sib families. Three of these SNPs were located introns of genes: MET (BTA04), RNF122 (BTA27) and WRN (BTA27). In summary, structural polymorphisms in form of CNVs, putatively play a role in susceptibility to CM. Specifically, sequence deletions have a greater effect on reducing resistance against mastitis, than sequence duplications have on increasing resistance against the disease.

Characterization of lncRNA-miRNA-mRNA Network to Reveal Potential Functional ceRNAs in Bovine Skeletal Muscle

There is growing evidence that non-coding RNAs are emerging as critical regulators of skeletal muscle development. In order to reveal their functional roles and regulatory mechanisms, we constructed a lncRNA–miRNA–mRNA network according to the ceRNA (competitive endogenous RNA) theory, using our high-throughput sequencing data. Subsequently, the network analysis, GO (Gene Ontology) analysis, and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway analysis were performed for functional annotation and exploration of lncRNA ceRNAs. The results uncovered a scale-free characteristics network which exhibited high functional specificity for bovine skeletal muscle development: co-expression lncRNAs were significantly enriched in muscle development related biological processes and the Wnt signaling pathway. Furthermore, GSEA (Gene Set Enrichment Analysis) indicated that the risk score has a tendency to associate with myogenesis, and differentially expressed RNAs were validated by qPCR, further confirming the credibility of our network. In summary, this study provides insights into lncRNA-mediated ceRNA function and mechanisms in bovine skeletal muscle development and will expand our understanding of lncRNA biology in mammals.

Analysis of long non-coding RNA and mRNA expression in bovine macrophages brings up novel aspects of Mycobacterium avium subspecies paratuberculosis infections

Paratuberculosis is a major disease in cattle that severely affects animal welfare and causes huge economic losses worldwide. Development of alternative diagnostic methods is of urgent need to control the disease. Recent studies suggest that long non-coding RNAs (lncRNAs) play a crucial role in regulating immune function and may confer valuable information about the disease. However, their role has not yet been investigated in cattle with respect to infection towards Paratuberculosis. Therefore, we investigated the alteration in genomic expression profiles of mRNA and lncRNA in bovine macrophages in response to Paratuberculosis infection using RNA-Seq. We identified 397 potentially novel lncRNA candidates in macrophages of which 38 were differentially regulated by the infection. A total of 820 coding genes were also significantly altered by the infection. Co-expression analysis of lncRNAs and their neighbouring coding genes suggest regulatory functions of lncRNAs in pathways related to immune response. For example, this included protein coding genes such as TNIP3, TNFAIP3 and NF-κB2 that play a role in NF-κB2 signalling, a pathway associated with immune response. This study advances our understanding of lncRNA roles during Paratuberculosis infection.

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

Long noncoding RNAs (lncRNAs) regulate the expression of mRNA and can affect various biological processes and phenotypes. Currently, studies of lncRNAs in cattle are under way, but their exact function for several tissues has not yet been established. Hanwoo cattle (Bos taurus coreanae) have inhabited the Korean peninsula for about 6000 years and are one of the representative domesticated animals in Korea. As a result of intensive breeding, the meat of Hanwoo cattle is high in marbling content and is preferred by Koreans and other East Asian people. In this study, the expression of lncRNAs was identified in 36 samples from skeletal muscle and three adipose tissues (intramuscular, subcutaneous, and omental) of nine Hanwoo individuals. We identified 76 tissue-specific lncRNAs for each of the four tissues using the differences in expression levels. Through QTL information, we could identify 12 lncRNAs associated with shear force and six lncRNAs associated with body weight, which are two important traits in the Hanwoo population breeding strategy. By the physical position comparison of lncRNA and Bovine transcripts information, we could identify 11 lncRNAs that were in bovine transcripts, and four of the 11 genes related to transcripts of lncRNAs were biologically associated with muscle function. We believe this Hanwoo lncRNAs study will help reveal the lncRNA role in the physiological mechanisms of these four tissues.

Sequencing the mosaic genome of Brahman cattle identifies historic and recent introgression including polled

Brahman cattle have a Bos indicus and Bos taurus mosaic genome, as a result of the process used to create the breed (repeat backcrossing of Bos taurus females to Bos indicus bulls). With the aim of identifying Bos taurus segments in the Brahman genome at sequence level resolution, we sequenced the genomes of 46 influential Brahman bulls. Using 36 million variants identified in the sequences, we searched for regions close to fixation for Bos indicus or Bos taurus segments that were longer than expected by chance (from simulation of the breed formation history of Brahman cattle). Regions close to fixation for Bos indicus content were enriched for protein synthesis genes, while regions of higher Bos taurus content included genes of the G-protein coupled receptor family (including genes implicated in puberty, such as THRS). The region with the most extreme Bos taurus enrichment was on chromosome 14 surrounding PLAG1. The introgressed Bos taurus allele at PLAG1 increases stature and the high frequency of the allele likely reflects strong selection for the trait. Finally, we provide evidence that the polled mutation in Brahmans, a desirable trait under very strong recent selection, is of Celtic origin and is introgressed from Bos taurus.

Genome-wide identification and characterization of long non-coding RNAs during postnatal development of rabbit adipose tissue

Background

The rabbit is widely used as an important experimental model for biomedical research, and shows low adipose tissue deposition during growth. Long non-coding RNAs (lncRNAs) are associated with adipose growth, but little is known about the function of lncRNAs in the rabbit adipose tissue.

Methods

Deep RNA-sequencing and comprehensive bioinformatics analyses were used to characterize the lncRNAs of rabbit visceral adipose tissue (VAT) at 35, 85 and 120 days after birth. Differentially expressed (DE) lncRNAs were identified at the three growth stages by DESeq. The cis and trans prediction ways predicted the target genes of the DE lncRNAs. To explore the function of lncRNAs, Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed on the candidate genes.

Results

A total of 991,157,544 clean reads were generated after RNA-Seq of the three growth stages, of which, 30,353 and 107 differentially expressed (DE) lncRNAs were identified. Compared to the protein-coding transcripts, the rabbit lncRNAs shared some characteristics such as shorter length and fewer exons. Cis and trans target gene prediction revealed, 43 and 64 DE lncRNAs respectively, corresponding to 72 and 20 protein-coding genes. GO enrichment and KEGG pathway analyses revealed that the candidate DE lncRNA target genes were involved in oxidative phosphorylation, glyoxylate and dicarboxylate metabolism, and other adipose growth-related pathways. Six DE lncRNAs were randomly selected and validated by q-PCR.

Conclusions

This study is the first to profile the potentially functional lncRNAs in the adipose tissue growth in rabbits, and contributes to our understanding of mammalian adipogenesis.

Electronic supplementary material

The online version of this article (10.1186/s12944-018-0915-1) contains supplementary material, which is available to authorized users.

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.

Coordinate regulation of long non-coding RNAs and protein-coding genes in germ-free mice

BACKGROUND

Long non-coding RNAs (lncRNAs) are increasingly recognized as regulators of tissue-specific cellular functions and have been shown to regulate transcriptional and translational processes, acting as signals, decoys, guides, and scaffolds. It has been suggested that some lncRNAs act in cis to regulate the expression of neighboring protein-coding genes (PCGs) in a mechanism that fine-tunes gene expression. Gut microbiome is increasingly recognized as a regulator of development, inflammation, host metabolic processes, and xenobiotic metabolism. However, there is little known regarding whether the gut microbiome modulates lncRNA gene expression in various host metabolic organs. The goals of this study were to 1) characterize the tissue-specific expression of lncRNAs and 2) identify and annotate lncRNAs differentially regulated in the absence of gut microbiome.

RESULTS

Total RNA was isolated from various tissues (liver, duodenum, jejunum, ileum, colon, brown adipose tissue, white adipose tissue, and skeletal muscle) from adult male conventional and germ-free mice (n = 3 per group). RNA-Seq was conducted and reads were mapped to the mouse reference genome (mm10) using HISAT. Transcript abundance and differential expression was determined with Cufflinks using the reference databases NONCODE 2016 for lncRNAs and UCSC mm10 for PCGs. Although the constitutive expression of lncRNAs was ubiquitous within the enterohepatic (liver and intestine) and the peripheral metabolic tissues (fat and muscle) in conventional mice, differential expression of lncRNAs by lack of gut microbiota was highly tissue specific. Interestingly, the majority of gut microbiota-regulated lncRNAs were in jejunum. Most lncRNAs were co-regulated with neighboring PCGs. STRING analysis showed that differentially expressed PCGs in proximity to lncRNAs form tissue-specific networks, suggesting that lncRNAs may interact with gut microbiota/microbial metabolites to regulate tissue-specific functions.

CONCLUSIONS

This study is among the first to demonstrate that gut microbiota critically regulates the expression of lncRNAs not only locally in intestine but also remotely in other metabolic organs, suggesting that common transcriptional machinery may be shared to transcribe lncRNA-PCG pairs, and lncRNAs may interact with PCGs to regulate tissue-specific pathways.

Transcriptome Analysis of Long Non-Coding RNA in the Bovine Mammary Gland Following Dietary Supplementation with Linseed Oil and Safflower Oil

This study aimed to characterize the long non-coding RNA (lncRNA) expression in the bovine mammary gland and to infer their functions in dietary response to 5% linseed oil (LSO) or 5% safflower oil (SFO). Twelve cows (six per treatment) in mid lactation were fed a control diet for 28 days followed by a treatment period (control diet supplemented with 5% LSO or 5% SFO) of 28 days. Mammary gland biopsies were collected from each animal on day-14 (D-14, control period), D+7 (early treatment period) and D+28 (late treatment period) and were subjected to RNA-Sequencing and subsequent bioinformatics analyses. Functional enrichment of lncRNA was performed via potential cis regulated target genes located within 50 kb flanking regions of lncRNAs and having expression correlation of >0.7 with mRNAs. A total of 4955 lncRNAs (325 known and 4630 novel) were identified which potentially cis targeted 59 and 494 genes in LSO and SFO treatments, respectively. Enrichments of cis target genes of lncRNAs indicated potential roles of lncRNAs in immune function, nucleic acid metabolism and cell membrane organization processes as well as involvement in Notch, cAMP and TGF-β signaling pathways. Thirty-two and 21 lncRNAs were differentially expressed (DE) in LSO and SFO treatments, respectively. Six genes (KCNF1, STARD13, BCL6, NXPE2, HHIPL2 and MMD) were identified as potential cis target genes of six DE lncRNAs. In conclusion, this study has identified lncRNAs with potential roles in mammary gland functions and potential candidate genes and pathways via which lncRNAs might function in response to LSO and SFA.

Functional Annotation of Animal Genomes (FAANG): Current Achievements and Roadmap

Functional annotation of genomes is a prerequisite for contemporary basic and applied genomic research, yet farmed animal genomics is deficient in such annotation. To address this, the FAANG (Functional Annotation of Animal Genomes) Consortium is producing genome-wide data sets on RNA expression, DNA methylation, and chromatin modification, as well as chromatin accessibility and interactions. In addition to informing our understanding of genome function, including comparative approaches to elucidate constrained sequence or epigenetic elements, these annotation maps will improve the precision and sensitivity of genomic selection for animal improvement. A scientific community-driven effort has already created a coordinated data collection and analysis enterprise crucial for the success of this global effort. Although it is early in this continuing process, functional data have already been produced and application to genetic improvement reported. The functional annotation delivered by the FAANG initiative will add value and utility to the greatly improved genome sequences being established for domesticated animal species.