Identification and characterization of long intergenic noncoding RNAs in bovine mammary glands

Exploration of non-coding RNAs related to intramuscular fat deposition Xinjiang Brown cattle and Angus × Wagyu cattle

Non-coding RNAs (ncRNAs) serve as crucial regulatory elements in the process of adipogenesis in animals; however, the specific roles and interrelationships of ncRNAs in bovine fat deposition remain poorly understood. This study aims to investigate the differentially expressed ncRNAs in the longissimus dorsi muscle of Xinjiang Brown cattle (XB) and Angus × Wagyu cattle (AW), to elucidate the regulatory mechanisms underlying lipidogenesis that may involve ncRNAs. Four Xinjiang Brown cattle and four Angus × Wagyu cattle were selected, ensuring they are subjected to identical feeding conditions, in order to evaluate the intermuscular fat (IMF) of longissimus dorsi muscles. The fat content of muscle tissue was quantified using the Soxhlet extraction method, revealing that the fat levels in the AW group were significantly elevated compared to those in the XB group. Taking muscle samples for paraffin sectioning and observing their morphology, it was found that the fat richness of the AW group was significantly higher than that of the XB group. Utilizing high-throughput RNA sequencing technology, we conducted an extensive transcriptomic analysis of longissimus dorsi muscles of XB and AW to identify significant ncRNAs implicated in fat metabolism and adipogenesis. The miRNA analysis yielded between 109,343,831 117,258,570 clean reads, whereas the lncRNA and circRNA analyses produced between 81,607,756 102,917,174 clean reads. Subsequent analysis revealed the identification of 53 differentially expressed miRNAs, 176 differentially expressed lncRNAs, and 234 differentially expressed circRNAs. KEGG enrichment analysis revealed that the target genes of differentially expressed miRNAs, lncRNAs, and circRNAs are significantly enriched in 2, 17, and 22 distinct pathways, respectively. The pathways associated with the differential enrichment of miRNA target genes involve processes such as phosphorylation and protein modification. Concurrently, the pathways linked to the varying enrichment of lncRNA target genes encompass G protein-coupled receptor signaling, regulation of cell death and apoptosis, activities related to GTPase activation, and functions governing nucleotide triphosphatases, among others. The circRNA exhibiting differential expression are significantly enriched in a variety of biological processes, including signal transduction, nucleic acid synthesis, cellular architecture, GTPase activation, and phosphatase activities, among others. The analysis of the ncRNA interaction network suggests that AGBL1, THRB, and S100A13 may play pivotal roles in the formation and adipogenic differentiation of adipocytes. In conclusion, we conducted a comprehensive analysis and discussion of the complete transcriptome of intermuscular fat tissue from the longissimus dorsi muscles in Xinjiang Brown cattle and Angus × Wagyu cattle. This study provides a theoretical foundation for enhancing our understanding of the molecular mechanisms underlying fat metabolism and deposition in beef cattle.

Predictive value of cuproptosis and disulfidptosis-related lncRNA in head and neck squamous cell carcinoma prognosis and treatment

Objective

Head and neck squamous cell carcinoma (HNSCC) is a highly lethal and prevalent malignant tumor with a poor prognosis due to its high recurrence rate, This study aims to develop a prognostic index for HNSCC patients based on Cuproptosis and Disulfidptosis-related long noncoding RNA.

Methods

Gene expression and clinical data for HNSCC were obtained from The Cancer Genome Atlas (TCGA). Using Lasso regression and multivariate Cox regression, we established a risk scoring model. The predictive ability of the nomogram, based on clinical features and risk scores, was verified using receiver operating characteristics and calibration curves. We compared independent prognostic parameters, risk score distribution, and survival between high-risk and low-risk groups, followed by preliminary validity evaluations of the model.

Results

Our systematic evaluation of prognostic risk provides a new direction for improving the survival prognosis of HNSCC patients in clinical practice, The model effectively categorized patients into high- and low-risk groups with distinct outcomes, identifying numerous gene mutations in these groups, A low-risk score was associated with a better prognosis and higher survival rates.

Conclusion

The risk score prognostic prediction system developed in this study shows potential efficacy in predicting the prognosis of HNSCC patients and has practical applications in clinical settings.

Regulation of inflammatory response by LINC00346 via miR-25-3p-mediated modulation of the PTEN/PI3K/AKT/NF-κB pathway

Long intergenic non-coding RNA 346 (LINC00346) has been reported to be involved in the development of atherosclerosis and specific cancers by affecting signaling pathways. However, its function in inflammation has not been thoroughly studied. Therefore, its expression pattern and function were determined in the human macrophage-like cell line THP-1. Lipopolysaccharide (LPS) treatment induced the expression of LINC00346. LPS-induced NF-κB activation and proinflammatory cytokine expression were suppressed or enhanced by the overexpression or knockdown of LINC00346, respectively. Analyses using dual luciferase assay and decoy RNAs that could block RNA-RNA interactions indicated that LINC00346 improves phosphatase and tensin homolog (PTEN) expression by sponging miR-25-3p. Subsequently, PTEN suppresses phosphoinositide-3 kinase (PI3K)-mediated conversion of phosphatidylinositol-4,5-bisphosphate (PIP2) into phosphatidylinositol-3,4,5-trisphosphate (PIP3) as well as consequent activation of protein kinase B (AKT) and NF-κB. Interestingly, database analysis revealed that the expression levels of LINC00346 and PTEN were simultaneously decreased in breast cancer tissues. Further analyses conducted using a breast cancer cell line, MDA-MB-231, confirmed the functional relationship among LINC00346, miR-25-3p, and PTEN in LPS-induced activation of NF-κB. These results indicate that miR-25-3p-sponging activity of LINC00346 affects the balance between PTEN and PI3K as well as the downstream activation of AKT/NF-κB pathway in inflammatory conditions.

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.

Extracellular Vesicles Derived from Selenium-Deficient MAC-T Cells Aggravated Inflammation and Apoptosis by Triggering the Endoplasmic Reticulum (ER) Stress/PI3K-AKT-mTOR Pathway in Bovine Mammary Epithelial Cells

Selenium (Se) deficiency disrupts intracellular REDOX homeostasis and severely deteriorates immune and anti-inflammatory function in high-yielding periparturient dairy cattle. To investigate the damage of extracellular vesicles derived from Se-deficient MAC-T cells (SeD-EV) on normal mammary epithelial cells, an in vitro model of Se deficiency was established. Se-deficient MAC-T cells produced many ROS, promoting apoptosis and the release of inflammatory factors. Extracellular vesicles were successfully isolated by ultrahigh-speed centrifugation and identified by transmission electron microscopy, particle size analysis, and surface markers (CD63, CD81, HSP70, and TSG101). RNA sequencing was performed on exosomal RNA. A total of 9393 lncRNAs and 63,155 mRNAs transcripts were identified in the SeC and SeD groups, respectively, of which 126 lncRNAs and 955 mRNAs were differentially expressed. Furthermore, SeD-EV promoted apoptosis of normal MAC-T cells by TUNEL analysis. SeD-EV significantly inhibited Bcl-2, while Bax and Cleaved Caspase3 were greatly increased. Antioxidant capacity (CAT, T-AOC, SOD, and GSH-Px) was inhibited in SeD-EV-treated MAC-T cells. Additionally, p-PERK, p-eIF2α, ATF4, CHOP, and XBP1 were all elevated in MAC-T cells supplemented with SeD-EV. In addition, p-PI3K, p-Akt, and p-mTOR were decreased strikingly by SeD-EV. In conclusion, SeD-EV caused oxidative stress, thus triggering apoptosis and inflammation through endoplasmic reticulum stress and the PI3K-Akt-mTOR signaling pathway, which contributed to explaining the mechanism of Se deficiency causing mastitis.

High-throughput analysis of lncRNA in cows with naturally infected Staphylococcus aureus mammary gland

LncRNA (long non-coding RNA) is an RNA molecule with a length between 200 and 100,000 nt. It does not encode proteins and is involved in a variety of intracellular processes, becoming a research hotspot of genetics. To identify key lncRNAs associated with dairy mastitis, we collected mammary epithelial tissue samples of Normal disease-free Holstein cows (HCN) and unhealthy Holstein cows with Staphylococcus aureus (HCU) and performed RNA sequencing (RNA-seq) on the samples. A total of 270 differentially expressed lncRNAs and 500 differentially expressed mRNAs were identified by high-throughput sequencing and bioinformatics analysis. Furthermore, Hydrolase activity is the most enriched in GO, and ErbB signaling pathway is significantly enriched in KEGG. In addition, through qPCR validation of 5 candidate lncRNAs in HCN and HCU, four differentially expressed lncRNAs MSTRG.498, MSTRG57.1, MSTRG.41.1 and MSTRG 124.1 were confirmed to have significant differentially expressed in cow mastitis. Also, lncRNA MSTRG.498 and its target gene, SMC4, might directly or indirectly play a role in cow mastitis. The regulatory network of lncRNA-miRNA-mRNA has been inferred from a bioinformatics perspective, which may assist understand the underlying molecular mechanism of lncRNAs involved in regulating mastitis in cows. Our findings will provide meaningful resources for further research on the regulatory function of lncRNAs in cow mastitis.

CeRNA Network Reveals the Circular RNA Characterization in Goat Ear Fibroblasts Reprogramming into Mammary Epithelial Cells

Circular RNAs (circRNAs) are a type of non-coding RNA that play a crucial role in the development and lactation of mammary glands in mammals. A total of 107 differentially expressed circRNAs (DE circRNAs) were found, of which 52 were up-regulated and 55 were down-regulated. We also found that DE circRNA host genes were mainly involved in GO terms related to the development process of mammary epithelial cells and KEGG pathways were mostly related to mammary epithelial cells, lactation, and gland development. Protein network analysis found that DE circRNAs can competitively bind to miRNAs as key circRNAs by constructing a circRNA-miRNA-mRNA network. CircRNAs competitively bind to miRNAs (miR-10b-3p, miR-671-5p, chi-miR-200c, chi-miR-378-3p, and chi-miR-30e-5p) involved in goat mammary gland development, mammary epithelial cells, and lactation, affecting the expression of core genes (CDH2, MAPK1, ITGB1, CAMSAP2, and MAPKAPK5). Here, we generated CiMECs and systematically explored the differences in the transcription profile for the first time using whole-transcriptome sequencing. We also analyzed the interaction among mRNA, miRNA, and cirRNA and predicted that circRNA plays an important role in the maintenance of mammary epithelial cells.

Long-chain noncoding RNA sequencing analysis reveals the molecular profiles of chemically induced mammary epithelial cells

Long noncoding RNAs (lncRNAs) were important regulators affecting the cellular reprogramming process. Previous studies from our group have demonstrated that small molecule compounds can induce goat ear fibroblasts to reprogram into mammary epithelial cells with lactation function. In this study, we used lncRNA-Sequencing (lncRNA-seq) to analyze the lncRNA expression profile of cells before and after reprogramming (CK vs. 5i8 d). The results showed that a total of 3,970 candidate differential lncRNAs were detected, 1,170 annotated and 2,800 new lncRNAs. Compared to 0 d cells, 738 lncRNAs were significantly upregulated and 550 were significantly downregulated in 8 d cells. Heat maps of lncrnas and target genes with significant differences showed that the fate of cell lineages changed. Functional enrichment analysis revealed that these differently expressed (DE) lncRNAs target genes were mainly involved in signaling pathways related to reprogramming and mammary gland development, such as the Wnt signaling pathway, PI3K-Akt signaling pathway, arginine and proline metabolism, ECM-receptor interaction, and MAPK signaling pathway. The accuracy of sequencing was verified by real-time fluorescence quantification (RT-qPCR) of lncRNAs and key candidate genes, and it was also demonstrated that the phenotype and genes of the cells were changed. Therefore, this study offers a foundation for explaining the molecular mechanisms of lncRNAs in chemically induced mammary epithelial cells.

LncRNAs Transcriptome Analysis Revealed Potential Mechanisms of Selenium to Mastitis in Dairy Cows

The trace element selenium (Se) plays an indispensable role in the growth of humans and animals due to its antioxidant function. Mastitis is one of the most important diseases affecting the dairy industry in the world. In recent years, long non-coding RNAs (lncRNAs) have been implicated in a series of cellular processes and disease development processes. RNA-sequencing technology was used to characterize lncRNA profiles and compared transcriptomic dynamics among the control group, the LPS group, and the Se-treated group to highlight the potential roles and functions of lncRNAs in the mammary epithelial cells of dairy cows. We identified 14 specific lncRNAs related to Se and their predicted target genes. KEGG and GO functional annotation was used to elucidate their biological function and the pathways in which they may be involved. The present study provides novel insights for exploring the molecular markers for the protection of Se against mastitis in dairy cows.

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.

Insight into mechanisms of pig lncRNA FUT3-AS1 regulating E. coli F18-bacterial diarrhea

Escherichia coli F18 is a common conditional pathogen that is associated with a variety of infections in humans and animals. LncRNAs have emerged as critical players in pathogen infection, but their role in the resistance of the host to bacterial diarrhea remains unknown. Here, we used piglets as animal model and identified an antisense lncRNA termed FUT3-AS1 as a host regulator related to E. coli F18 infection by RNA sequencing. Downregulation of FUT3-AS1 expression contributed to the enhancement of E. coli F18 resistance in IPEC-J2 cells. FUT3-AS1 knockdown reduced FUT3 expression via decreasing the H4K16ac level of FUT3 promoter. Besides, the FUT3-AS1/miR-212 axis could act as a competing endogenous RNA to regulate FUT3 expression. Functional analysis demonstrated that target FUT3 plays a vital role in the resistance of IPEC-J2 cells to E. coli F18 invasion. A Fut3-knockout mice model was established and Fut3-knockout mice obviously improved the ability of resistance to bacterial diarrhea. Interestingly, FUT3 could enhance E. coli F18 susceptibility by activating glycosphingolipid biosynthesis and toll-like receptor signaling which are related to receptor formation and immune response, respectively. In summary, we have identified a novel biomarker FUT3-AS1 that modulates E. coli F18 susceptibility via histone H4 modifications or miR-212/FUT3 axis, which will provide theoretical guidance to develop novel strategies for combating bacterial diarrhea in piglets.

Preliminary study on gene regulation and its pathways in Chinese Holstein cows with clinical mastitis caused by Staphylococcus aureus

Introduction

Clinical mastitis (CM) is one of the most common diseases of dairy cows globally, has a complex aetiology and recurs easily. Staphylococcus aureus is a frequently isolated pathogen responsible for bovine mastitis and remains difficult to eradicate.

Material and Methods

To characterise the transcriptional profiles of dairy cows infected by S. aureus, we performed an RNA-seq analysis of peripheral blood leukocytes in lactating Chinese Holstein dairy cows with CM and did the same with healthy cows’ samples as controls.

Results

A total of 4,286 genes were detected in the CM cases infected with S. aureus which were differentially expressed compared to the controls, 3,085 of which were upregulated, the remainder being downregulated. Notably, we observed that some differentially expressed genes (DEGs) had strong protein–protein interaction. Of these, six downregulated DEGs (AKR1C4, PTGS2, HNMT, EPHX2, CMBL, and IDH1) were involved in the metabolic pathway, while eight upregulated DEGs (VWF, GP9, MYLK, GP6, F2RL3, ITGB3, GP5, and PRKG1) were associated with the platelet activation pathway.

Conclusion

The transcriptome dataset of CM cases would be a valuable resource for clinical guidance on anti-inflammatory medication and for deeper understanding of the biological processes of CM response to S. aureus infection, and it would enable us to identify specific genes for diagnostic markers and possibly for targeted therapy.

Progress in Expression Pattern and Molecular Regulation Mechanism of LncRNA in Bovine Mastitis

Simple Summary

Bovine mastitis is an inflammatory disease of the mammary glands that causes serious harm to cow health and huge economic losses. Susceptibility or resistance to mastitis in individual cows is mainly determined by genetic factors, including coding genes and non-coding genes. Long non-coding RNAs (lncRNAs) are non-coding RNA molecules with a length of more than 200 nucleotides (nt) that have recently been discovered. They can regulate a variety of diseases of humans and animals, especially the immune response and inflammatory disease process. This paper reviews the role of long non-coding RNA (lncRNA) in inflammatory diseases, emphasizes on the latest research progress of lncRNA expression and the molecular regulatory mechanism in bovine mastitis, and looks forward to the research and application prospect of lncRNA in bovine mastitis, intending to provide a reference for scientific researchers to systematically understand this research field.

Abstract

Bovine mastitis is an inflammatory disease caused by pathogenic microbial infection, trauma, or other factors. Its morbidity is high, and it is difficult to cure, causing great harm to the health of cows and the safety of dairy products. Susceptibility or resistance to mastitis in individual cows is mainly determined by genetic factors, including coding genes and non-coding genes. Long non-coding RNAs (lncRNAs) are a class of endogenous non-coding RNA molecules with a length of more than 200 nucleotides (nt) that have recently been discovered. They can regulate the immune response of humans and animals on three levels (transcription, epigenetic modification, and post-transcription), and are widely involved in the pathological process of inflammatory diseases. Over the past few years, extensive findings revealed basic roles of lncRNAs in inflammation, especially bovine mastitis. This paper reviews the expression pattern and mechanism of long non-coding RNA (lncRNA) in inflammatory diseases, emphasizes on the latest research progress of the lncRNA expression pattern and molecular regulatory mechanism in bovine mastitis, analyzes the molecular regulatory network of differentially expressed lncRNAs, and looks forward to the research and application prospect of lncRNA in bovine mastitis, laying a foundation for molecular breeding and the biological therapy of bovine mastitis.

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.

Exosomal RNAs: Novel Potential Biomarkers for Diseases-A Review

Exosomes are a subset of nano-sized extracellular vesicles originating from endosomes. Exosomes mediate cell-to-cell communication with their cargos, which includes mRNAs, miRNAs, lncRNAs, and circRNAs. Exosomal RNAs have cell specificity and reflect the conditions of their donor cells. Notably, their detection in biofluids can be used as a diagnostic marker for various diseases. Exosomal RNAs are ideal biomarkers because their surrounding membranes confer stability and they are detectable in almost all biofluids, which helps to reduce trauma and avoid invasive examinations. However, knowledge of exosomal biomarkers remains scarce. The present review summarizes the biogenesis, secretion, and uptake of exosomes, the current researches exploring exosomal mRNAs, miRNAs, lncRNAs, and circRNAs as potential biomarkers for the diagnosis of human diseases, as well as recent techniques of exosome isolation.

Screening and Conjoint Analysis of Key lncRNAs for Milk Fat Metabolism in Dairy Cows

Long noncoding RNAs (lncRNAs) play an important regulatory role in various biological processes as a key regulatory factor. However, the complete expression profile of lncRNAs in dairy cows and its function in milk fat synthesis are unknown. In this study, RNA sequencing (RNA-seq) was used to research the whole genome expression of lncRNAs and mRNA transcripts in high and low milk fat percentage (MFP) bovine mammary epithelial cells (BMECs), and joint analysis was carried out. We identified a total of 47 differentially expressed genes (DEGs) and 38 differentially expressed lncRNAs (DELs, Padj &lt;0.05), enrichment analysis screened out 11 candidate DEGs that may regulate milk fat metabolism. Downregulated differential gene ENPP2 (The expression level in BMECs of high milk fat dairy cows was lower than that of low milk fat cows) and upregulated differential gene BCAT1 are more likely to participate in the milk fat metabolism, and its function needs further experiments verification. The enrichment analysis of target genes predicted by DELs identified 7 cis (co-localization) and 10 trans (co-expression) candidate target genes related to milk lipid metabolism, corresponding to a total of 18 DELs. Among them, the targeting relationship between long intervening/intergenic noncoding RNA (lincRNA) TCONS_00082721 and FABP4 is worthy of attention. One hundred and fifty-six competing endogenous RNAs (ceRNAs) interaction regulation networks related to milk fat metabolism were constructed based on the expression information of DELs, differential microRNAs (miRNAs), and lipid metabolism-related target genes. The regulatory network centered on miR-145 will be the focus of subsequent experimental research. The ceRNAs regulatory network related to TCONS_00082721 and TCONS_00172817 are more likely to be involved in milk fat synthesis. These results will provide new ways to understand the complex biology of dairy cow milk fat synthesis and provide valuable information for breed improvement of Chinese Holstein cow.

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.

Regulatory network of miRNA, lncRNA, transcription factor and target immune response genes in bovine mastitis

Pre- and post-transcriptional modifications of gene expression are emerging as foci of disease studies, with some studies revealing the importance of non-coding transcripts, like long non-coding RNAs (lncRNAs) and microRNAs (miRNAs). We hypothesize that transcription factors (TFs), lncRNAs and miRNAs modulate immune response in bovine mastitis and could potentially serve as disease biomarkers and/or drug targets. With computational analyses, we identified candidate genes potentially regulated by miRNAs and lncRNAs base pair complementation and thermodynamic stability of binding regions. Remarkably, we found six miRNAs, two being bta-miR-223 and bta-miR-24-3p, to bind to several targets. LncRNAs NONBTAT027932.1 and XR_003029725.1, were identified to target several genes. Functional and pathway analyses revealed lipopolysaccharide-mediated signaling pathway, regulation of chemokine (C-X-C motif) ligand 2 production and regulation of IL-23 production among others. The overarching interactome deserves further in vitro/in vivo explication for specific molecular regulatory mechanisms during bovine mastitis immune response and could lay the foundation for development of disease markers and therapeutic intervention.

Testis transcriptome profiling identified lncRNAs involved in spermatogenic arrest of cattleyak

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

Potential Novel Biomarkers for Mastitis Diagnosis in Sheep

Simple Summary

Inflammation of the mammary gland (mastitis) is an important disease of dairy sheep. Mastitis management depends mainly on the diagnosis. Conventional diagnostic methods including somatic cell count, California Mastitis Test, and microbial culture have limitations. Therefore researchers are looking for new diagnostic biomarkers of mastitis including specific proteins produced by the liver in case of disease (acute phase proteins), unique genetic sequences (miRNAs), or antimicrobial peptides produced by immune cells during inflammation (cathelicidines).

Abstract

This review aims to characterize promising novel markers of ovine mastitis. Mastitis is considered as one of the primary factors for premature culling in dairy sheep and has noticeable financial, productional, and animal welfare-related implications. Furthermore, clinical, and subclinical mammary infections negatively affect milk yield and alter the milk composition, thereby leading to lowered quality of dairy products. It is, therefore, crucial to control and prevent mastitis through proper diagnosis, treatment or culling, and appropriate udder health management particularly at the end of the lactation period. The clinical form of mastitis is characterized by abnormalities in milk and mammary gland tissue alteration or systemic symptoms consequently causing minor diagnostic difficulties. However, to identify ewes with subclinical mastitis, laboratory diagnostics is crucial. Mastitis control is primarily dependent on determining somatic cell count (SCC) and the California Mastitis Test (CMT), which aim to detect the quantity of cells in the milk sample. The other useful diagnostic tool is microbial culture, which complements SCC and CMT. However, all mentioned diagnostic methods have their limitations and therefore novel biomarkers of ovine subclinical mastitis are highly desired. These sensitive indicators include acute-phase proteins, miRNA, and cathelicidins measurements, which could be determined in ovine serum and/or milk and in the future may become useful in early mastitis diagnostics as well as a preventive tool. This may contribute to increased detection of ovine mammary gland inflammation in sheep, especially in subclinical form, and consequently improves milk quality and quantity.

Comparative Transcriptome Analysis of Milk Somatic Cells During Lactation Between Two Intensively Reared Dairy Sheep Breeds

In dairy sheep industry, milk production dictates the value of a ewe. Milk production is directly related to the morphology and physiology of the mammary gland; both being designated targets of breeding strategies. Although within a flock breeding parameters are mutual, large differences in milk production among individual ewes are usually observed. In this work, we tested two of the most productive dairy sheep breeds reared intensively in Greece, one local the Chios breed and one foreign the Lacaune breed. We used transcriptome sequencing to reveal molecular mechanisms that render the mammary gland highly productive or not. While highly expressed genes (caseins and major whey protein genes) were common among breeds, differences were observed in differentially expressed genes. ENSOARG00000008077, as a member of ribosomal protein 14 family, together with LPCAT2, CCR3, GPSM2, ZNF131, and ASIP were among the genes significantly differentiating mammary gland's productivity in high yielding ewes. Gene ontology terms were mainly linked to the inherent transcriptional activity of the mammary gland (GO:0005524, GO:0030552, GO:0016740, GO:0004842), lipid transfer activity (GO:0005319) and innate immunity (GO:0002376, GO:0075528, GO:0002520). In addition, clusters of genes affecting zinc and iron trafficking into mitochondria were highlighted for high yielding ewes (GO:0071294, GO:0010043). Our analyses provide insights into the molecular pathways involved in lactation between ewes of different performances. Results revealed management issues that should be addressed by breeders in order to move toward increased milk yields through selection of the desired phenotypes. Our results will also contribute toward the selection of the most resilient and productive ewes, thus, will strengthen the existing breeding systems against a spectrum of environmental threats.

Profiling and integrated analysis of whole-transcriptome changes in uterine caruncles of pregnant and non-pregnant buffaloes

Improved reproductive performance in buffaloes can be achieved by understanding the basic mechanism governing the embryonic attachment and feto-maternal communication. Considering this, trascriptomic profiling and integrative analysis of long intergenic non-coding RNAs were carried out in the uterine caruncles of pregnant and non-pregnant buffaloes. Transcriptome data of pregnant and non-pregnant uterine caruncles after quality control was used to perform the analysis. Total of 86 novel lincRNAs expressed in uterine caruncular tissues were identified and characterized. Differential expression analysis revealed that 447 mRNAs and 185 mRNAs were up- and down- regulated, respectively. The number of up- and down- regulated lincRNAs were 114 and 13, respectively. Of the identified 86 novel lincRNAs, six novel lincRNAs were up-regulated in the pregnant uterine caruncles. GO terms (biological process) and PANTHER pathways associated with reproduction and embryogenesis were over-represented in differentially expressed genes. Through miRNA interaction analysis, interactions of 16 differentially expressed lincRNAs with mi-RNAs involved in reproduction were identified. This study has provided a catalogue of differentially expressed genes and novel regions previously unknown to play a significant role in buffalo reproduction. The results from the current study extends the buffalo uterine lncRNAs database and provides candidate regulators for future molecular genetic studies on buffalo uterine physiology to improve the embryo implantation and successful completion of pregnancy.

Long Intergenic Non-Coding RNAs in the Mammary Parenchyma and Fat Pad of Pre-Weaning Heifer Calves: Identification and Functional Analysis

Enhanced plane of nutrition at pre-weaning stage can promote the development of mammary gland especially heifer calves. Although several genes are involved in this process, long intergenic non-coding RNAs (lincRNAs) are regarded as key regulators in the regulated network and are still largely unknown. We identified and characterized 534 putative lincRNAs based on the published RNA-seq data, including heifer calves in two groups: fed enhanced milk replacer (EH, 1.13 kg/day, including 28% crude protein, 25% fat) group and fed restricted milk replacer (R, 0.45 kg/day, including 20% crude protein, 20% fat) group. Sub-samples from the mammary parenchyma (PAR) and mammary fat pad (MFP) were harvested from heifer calves. According to the information of these lincRNAs' quantitative trait loci (QTLs), the neighboring and co-expression genes were used to predict their function. By comparing EH vs R, 79 lincRNAs (61 upregulated, 18 downregulated) and 86 lincRNAs (54 upregulated, 32 downregulated) were differentially expressed in MFP and PAR, respectively. In MFP, some differentially expressed lincRNAs (DELs) are involved in lipid metabolism pathways, while, in PAR, among of DELs are involved in cell proliferation pathways. Taken together, this study explored the potential regulatory mechanism of lincRNAs in the mammary gland development of calves under different planes of nutrition.

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.

Transcription initiation mapping in 31 bovine tissues reveals complex promoter activity, pervasive transcription, and tissue-specific promoter usage

Characterizing transcription start sites is essential for understanding the regulatory mechanisms that control gene expression. Recently, a new bovine genome assembly (ARS-UCD1.2) with high continuity, accuracy, and completeness was released; however, the functional annotation of the bovine genome lacks precise transcription start sites and contains a low number of transcripts in comparison to human and mouse. By using the RAMPAGE approach, this study identified transcription start sites at high resolution in a large collection of bovine tissues. We found several known and novel transcription start sites attributed to promoters of protein-coding and lncRNA genes that were validated through experimental and in silico evidence. With these findings, the annotation of transcription start sites in cattle reached a level comparable to the mouse and human genome annotations. In addition, we identified and characterized transcription start sites for antisense transcripts derived from bidirectional promoters, potential lncRNAs, mRNAs, and pre-miRNAs. We also analyzed the quantitative aspects of RAMPAGE to produce a promoter activity atlas, reaching highly reproducible results comparable to traditional RNA-seq. Coexpression networks revealed considerable use of tissue-specific promoters, especially between brain and testicle, which expressed several genes in common from alternate loci. Furthermore, regions surrounding coexpressed modules were enriched in binding factor motifs representative of each tissue. The comprehensive annotation of promoters in such a large collection of tissues will substantially contribute to our understanding of gene expression in cattle and other mammalian species, shortening the gap between genotypes and phenotypes.

Epigenetics: New Insights into Mammary Gland Biology

The mammary gland undergoes important anatomical and physiological changes from embryogenesis through puberty, pregnancy, lactation and involution. These steps are under the control of a complex network of molecular factors, in which epigenetic mechanisms play a role that is increasingly well described. Recently, studies investigating epigenetic modifications and their impacts on gene expression in the mammary gland have been performed at different physiological stages and in different mammary cell types. This has led to the establishment of a role for epigenetic marks in milk component biosynthesis. This review aims to summarize the available knowledge regarding the involvement of the four main molecular mechanisms in epigenetics: DNA methylation, histone modifications, polycomb protein activity and non-coding RNA functions.

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.

Transcriptome Analysis Reveals the Potential Role of Long Non-coding RNAs in Mammary Gland of Yak During Lactation and Dry Period

The mammary gland is a remarkably dynamic organ of milk synthesis and secretion, and it experiences drastic structural and metabolic changes during the transition from dry periods to lactation, which involves the expression and regulation of numerous genes and regulatory factors. Long non-coding RNA (lncRNA) has considered as a novel type of regulatory factors involved in a variety of biological processes. However, their role in the lactation cycle of yak is still poorly understood. To reveal the involved mechanism, Ribo-zero RNA sequencing was employed to profile the lncRNA transcriptome in mammary tissue samples from yak at two physiological stages, namely lactation (LP) and dry period (DP). Notably, 1,599 lncRNA transcripts were identified through four rigorous steps and filtered through protein-coding ability. A total of 59 lncRNAs showed significantly different expression between two stages. Accordingly, the results of qRT-PCR were consistent with that of the transcriptome data. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses indicated that target genes of differentially expressed lncRNAs (DELs) were involved in pathways related to lactation, such as ECM-receptor interaction, PI3K-Akt signaling pathway, biosynthesis of amino acids and focal adhesion etc. Finally, we constructed a lncRNA-gene regulatory network containing some well known candidate genes for milk yield and quality traits. This is the first study to demonstrate a global profile of lncRNA expression in the mammary gland of yak. These results contribute to a valuable resource for future genetic and molecular studies on improving milk yield and quality, and help us to gain a better understanding of the molecular mechanisms underlying lactogenesis and mammary gland development of yak.

Long Non-Coding RNAs in Brown Adipose Tissue

Obesity has become a widespread disease that is harmful to human health. Fat homeostasis is essentially maintained by fat accumulation and energy expenditure. Studies on brown adipose tissue (BAT) represent a promising opportunity to identify a pharmaceutical intervention against obesity through increased energy expenditure. Long non-coding RNAs (lncRNAs) were thought to be critical regulators in a variety of biological processes. Recent studies have revealed that lncRNAs, including ones that are BAT-specific, conserved, and located at key protein-coding genes, function in brown adipogenesis, white adipose browning (ie, beige adipogenesis), and brown thermogenesis. In this review, we describe lncRNA properties and highlight functional lncRNAs in these biological processes, with the goal of establishing links between lncRNAs and BAT. Based on the advances of lncRNAs in the regulation of BAT, we discussed the advantages of potential lncRNA-based obesity drugs. Further BAT lncRNA-based drug development may provide new exciting approaches to defend obesity by regulation of fat homeostasis.

RNA-Seq Reveals the Expression Profiles of Long Non-Coding RNAs in Lactating Mammary Gland from Two Sheep Breeds with Divergent Milk Phenotype

Simple Summary

Long non-coding RNAs (lncRNAs) play a key role in regulating the expression level of mRNAs. The expression profiles of ovine mammary gland were investigated in two sheep breeds with divergent milk phenotype using RNA-Seq. A total of 1894 lncRNAs were found to be expressed and 68 of these were differentially expressed between the two breeds. Some important Gene Ontogeny (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways that were related to lactation and mammary gland morphogenesis were found for the target genes of differentially expressed lncRNAs. This study can improve our understanding of the functions of lncRNAs in the regulation of lactation, milk yield, and milk components in sheep.

Abstract

Long non-coding RNAs (lncRNAs) are a kind of non-coding RNA with >200 nucleotides in length. Some lncRNAs have been proven to have clear regulatory functions in many biological processes of mammals. However, there have been no reports on the roles of lncRNAs in ovine mammary gland tissues. In the study, the expression profiles of lncRNAs were studied using RNA-Seq in mammary gland tissues from lactating Small-Tailed Han (STH) ewes and Gansu Alpine Merino (GAM) ewes with different milk yield and ingredients. A total of 1894 lncRNAs were found to be expressed. Compared with the GAM ewes, the expression levels of 31 lncRNAs were significantly up-regulated in the mammary gland tissues of STH ewes, while 37 lncRNAs were remarkably down-regulated. Gene Ontogeny (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis found that the target genes of differentially expressed lncRNAs were enriched in the development and proliferation of mammary epithelial cells, morphogenesis of mammary gland, ErbB signaling pathway, and Wnt signaling pathway. Some miRNA sponges of differentially expressed lncRNAs, reported to be associated with lactation and mammary gland morphogenesis, were found in a lncRNA-miRNA network. This study reveals comprehensive lncRNAs expression profiles in ovine mammary gland tissues, thereby providing a further understanding of the functions of lncRNAs in the lactation and mammary gland development of sheep.

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.

Circ09863 Regulates Unsaturated Fatty Acid Metabolism by Adsorbing miR-27a-3p in Bovine Mammary Epithelial Cells

Fatty acid composition plays a key role in regulating flavor and quality of milk. Therefore, in order to improve milk quality, it is particularly important to investigate regulatory mechanisms of milk fatty acid metabolism. Circular RNAs (circRNAs) regulate expression genes associated with several biological processes including fatty acid metabolism. In this study, high-throughput sequencing was used to detect differentially expressed genes in bovine mammary tissue at early lactation and peak lactation. Circ09863 profiles were influenced by the lactation stage. Functional studies in bovine mammary epithelial cells (BMECs) revealed that circ09863 promotes triglyceride (TAG) synthesis together with increased content of unsaturated fatty acids (C16:1 and C18:1). These results suggested that circ09863 is partly responsible for modulating fatty acid metabolism. Additionally, software prediction identified a miR-27a-3p binding site in the circ09863 sequence. Overexpression of miR-27a-3p in BMECs led to decreased TAG synthesis. However, overexpression of circ09863 (pcDNA-circ09863) in BMECs significantly reduced expression of miR-27a-3p and enhanced gene expression of fatty acid synthase (FASN), a target of miR-27a-3p. Overall, data suggest that circ09863 relieves the inhibitory effect of miR-27a-3p on FASN expression by binding miR-27a-3p and subsequently regulating TAG synthesis and fatty acid composition. Together, these mechanisms provide new research avenues and theoretical bases to improve milk quality.

Identification and characterization of putative ovarian lincRNAs in dairy goats treated for repeated estrous synchronization

The goal of this study was to identify and characterize effects of repeated estrous synchronization (ES) treatments on the regulation of ovarian intergenic long non-coding RNAs (lincRNAs) in dairy goats. Six does were randomly assigned to a group administered three ES treatment regimens separated by 2 weeks or to a group administered only one ES treatment regimen (control) at the same time as the third ES treatment in the does administered the three hormonal regimens for ES. The paired-end RNA Sequencing procedures were used to evaluate lincRNAs of ovarian tissues. A total of 134 lincRNAs were differentially abundant between the two treatment groups. Several target genes were annotated and were related to hormone activity, cellular response to hormone stimulus, response to hormone, female pregnancy, as well as regulation of hormone secretion. These genes were noticeably enriched in MAPK, Hippo, estrogen signaling pathways, oocyte meiosis, progesterone-mediated oocyte maturation, ovarian steroidogenesis as well as GnRH signaling pathways. According to the enriched GO terms and KEGG pathways of regulated genes, 13 differentially abundant lincRNAs could be promising candidates for regulating reproductive functions of female goats. Current results indicate that repeated treatments with gonadotropins affected hormone sensitivity, estrogen synthesis, and ovarian function. The results also indicated that when there was imposing of the three hormonal treatment regimens for ES, there were several lincRNAs that could contribute to dysregulation of several genes that are important for reproduction in dairy goats. Findings provide novel insights for further investigation of lncRNAs biological functions in goats.

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 identification and characterization of long non-coding RNAs during differentiation of visceral preadipocytes in rabbit

Emerging evidence suggests that long non-coding RNAs (lncRNAs) are critical regulators of diverse biological processes, including adipogenesis. Despite being considered an ideal animal model for studying adipogenesis, little is known about the roles of lncRNAs in the regulation of rabbit preadipocyte differentiation. In the present study, visceral preadipocytes isolated from newborn rabbits were cultured in vitro and induced for differentiation, and global lncRNA expression profiles of adipocytes collected at days 0, 3, and 9 of differentiation were analyzed by RNA-seq. A total of 2066 lncRNAs were identified from nine RNA-seq libraries. Compared to protein-coding transcripts, lncRNA transcripts exhibited characteristics of a longer length and lower expression level. Furthermore, 486 and 357 differentially expressed (DE) lncRNAs were identified when comparing day 3 vs. day 0 and day 9 vs. day 3, respectively. Target genes of DE lncRNAs were predicted by the cis-regulating approach. Prediction of functions revealed that DE lncRNAs when comparing day 3 vs. day 0 were involved in gene ontology (GO) terms of developmental growth, growth, developmental cell growth, and stem cell proliferation, and involved in Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways of PI3K-Akt signaling pathway, fatty acid biosynthesis, and the insulin signaling pathway. The DE lncRNAs when comparing day 9 vs. day 3 were involved in GO terms that associated with epigenetic modification and were involved in the KEGG pathway of cAMP signaling pathway. This study provides further insight into the regulatory function of lncRNAs in rabbit visceral adipose and facilitates a better understanding of different stages of preadipocyte differentiation.

Long non-coding RNA profiling in LPS-induced intestinal inflammation model: New insight into pathogenesis

LPS can induce an inflammatory immune response in the intestine, and long non-coding RNA (lncRNA) is involved in the process of inflammatory disease. However, the biological role of lncRNA in the intestinal inflammation of piglets remains unclear. In this study, the lncRNA expression profile of the ileal mucosa of piglets challenged by LPS was analysed using lncRNA sequencing. In total, 112 novel lncRNAs were predicted, of which 58 were up-regulated and 54 down-regulated following LPS challenge. Expression of 15 selected lncRNAs was validated by quantitative PCR. We further investigated the target genes of lncRNA that were enriched in the signalling pathways involved in the inflammatory immune response by utilising Gene Ontology and Kyoto Encyclopaedia of Genes and Genomes analysis, with cell adhesion molecules and mTOR signalling pathway identified. In addition, the co-expression networks between the differentially expressed lncRNAs and the target mRNAs were constructed, with seven core lncRNAs identified, which also demonstrated that the relationship between lncRNAs and the target genes was highly correlated. Our study offers important information about the lncRNAs of the mucosal immune system in piglets and provides new insights into the inflammatory mechanism of LPS challenge, which might serve as a novel target to control intestinal inflammation.

Abnormally expressed long noncoding RNA B3GALT5-AS1 may serve as a biomarker for the diagnostic and prognostic of gastric cancer

Early diagnosis of gastric cancer (GC) is an effective method to improve prognosis. Increasing number of long noncoding RNAs (lncRNAs) have been reported as biomarkers for several cancers. We aim to detect the level of lncRNA B3GALT5-AS1 and its association with clinical parameters and to further explore its application value in GC. We measured serum B3GALT5-AS1 expression in 107 patients with GC, 40 polyp patients, and 87 normal controls to explore the significance of serum B3GALT5-AS1 in GC using the quantitative real-time polymerase chain reaction method. The result demonstrated that B3GALT5-AS1 level was markedly richer in GC patients than that in normal people (P < .001). B3GALT5-AS1 may be served as a diagnostic marker for distinguishing GC patients from healthy people, and the proportion under the receiver operating characteristics curve is 0.816 (95% confidence interval, 0.758-0.874; P = .03). Further exploration validated that high serum B3GALT5-AS1 level was related to TNM stage (P = .024), and lymph node metastasis (P = .023). Our study suggested that serum B3GALT5-AS1 may be employed as an ideal biomarker for early screening of GC.

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.

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.

Integrated analysis of long noncoding RNA and mRNA expression profiles reveals the potential role of long noncoding RNA in different bovine lactation stages

Long noncoding RNA (lncRNA) play a critical role in mammary development and breast cancer biology. Despite their important role in the mammary gland, little is known of the roles of lncRNA in bovine lactation, particularly regarding the molecular processes underlying it. To characterize the role of lncRNA in bovine lactation, 4 samples of Holstein cow mammary gland tissue at peak and late lactation stages were examined after biopsy. We then profiled the transcriptome of the mammary gland using RNA sequencing technology. Further, functional lncRNA-mRNA coexpression pairs were constructed to infer the function of lncRNA using a generalized linear model, followed by gene ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses. More than 1,000 putative lncRNA were identified, 117 of which were differentially expressed between peak and late lactation stages. Bovine lncRNA were shorter, with fewer exon numbers, and expressed at significantly lower levels than protein-coding genes. Seventy-two differentially expressed (DE) lncRNA were coexpressed with 340 different protein-coding genes. The KEGG pathway analysis showed that target mRNA for DE lncRNA were mainly related to lipid and glucose metabolism, including the peroxisome proliferator-activated receptors and 5' adenosine monophosphate-activated protein kinase signaling pathways. Further bioinformatics and integrative analyses revealed that 12 DE lncRNA potentially played important roles in bovine lactation. Our findings provide a valuable resource for future bovine transcriptome studies, facilitate the understanding of bovine lactation biology, and offer functional information for cattle lactation.

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

Background

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

Results

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

Conclusions

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

Electronic supplementary material

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

Transcriptome sequencing to detect the potential role of long non-coding RNAs in bovine mammary gland during the dry and lactation period

BACKGROUND

It is known that long non-coding RNAs (lncRNAs) play an important role in various biological processes, including cell proliferation, differentiation and apoptosis. However, their functions and profiles in lactation cycle of dairy cows are largely unknown. In this study, lncRNA-seq technique was employed to compare the expression profiles of lncRNAs and mRNAs from Chinese Holstein mammary gland in dry and lactation period.

RESULT

Totally 3746 differentially expressed lncRNAs (DELs) and 2890 differentially expressed genes (DEGs) were identified from the dry and lactation mammary glands of Holstein cows. Functional enrichment analysis on target genes of lncRNAs indicated that these genes were involved in lactation-related signaling pathways, including cell cycle, JAK-STAT, cell adhesion, and PI3K-Akt signaling pathways. Additionally, the interaction between lncRNAs and their potential miRNAs was predicted and partly verified. The result indicated that the lactation-associated miR-221 might interact with lncRNAs TCONS_00040268, TCONS_00137654, TCONS_00071659 and TCONS_00000352, which revealed that these lncRNAs might be important regulators for lactation cycle.

CONCLUSION

This study provides a resource for lncRNA research on lactation cycle of bovine mammary gland. Besides, the interaction between lncRNAs and the specific miRNA is revealed. It expands our knowledge about lncRNA and miRNA biology as well as contributes to clarify the regulation of lactation cycle of bovine mammary gland.

Genome Wide Identification of Novel Long Non-coding RNAs and Their Potential Associations With Milk Proteins in Chinese Holstein Cows

Long non-coding RNAs (lncRNAs) have emerged as a novel class of regulatory molecules involved in various biological processes. However, their role in milk performance is unknown. Here, whole transcriptome RNA sequencing was used to generate the lncRNA transcriptome profiles in mammary tissue samples from 6 Chinese Holstein cows with 3 extremely high and 3 low milk protein percentage phenotypes. In this study, 6,450 lncRNA transcripts were identified through 5 stringent steps and filtration by coding potential. In total, 31 lncRNAs and 18 novel genes were identified to be differentially expressed in high milk protein samples (HP) relative to low milk protein samples (LP), respectively. Differentially expressed lncRNAs were selected to predict target genes through bioinformatics analysis, followed by the integration of differentially expressed mRNA data, gene function, gene ontology (GO) and pathway, genome wide association study (GWAS) and quantitative trait locus (QTL) information, as well as network analysis to further characterize potential interactions. Several lncRNAs were found (such as XLOC_059976) that could be used as candidate markers for milk protein content prediction. This is the first study to perform global expression profiling of lncRNAs and mRNAs related to milk protein traits in dairy cows. These results provide important information and insights into the synthesis of milk proteins, and potential targets for the future improvement of milk quality.

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

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

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

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

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

Deciphering the roles of lncRNAs in breast development and disease

Breast cancer is the second leading cause of cancer related deaths in women. It is therefore important to understand the mechanisms underlying breast cancer development as well as raises the need for enhanced, non-invasive strategies for novel prognostic and diagnostic methods. The emergence of long non-coding RNAs (lncRNAs) as potential key players in neoplastic disease has received considerable attention over the past few years. This relatively new class of molecular regulators has been shown from ongoing research to act as critical players for key biological processes. Deregulated expression levels of lncRNAs have been observed in a number of cancers including breast cancer. Furthermore, lncRNAs have been linked to breast cancer initiation, progression, metastases and to limit sensitivity to certain targeted therapeutics. In this review we provide an update on the lncRNAs associated with breast cancer and mammary gland development and illustrate the versatility of such lncRNAs in gene control, differentiation and development both in normal physiological conditions and in diseased states. We also highlight the therapeutic and diagnostic potential of lncRNAs in cancer.