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

Transcriptome identification of differential mammary genes of Kazakh horses during early pregnancy

Kazakh mares have attracted widespread attention with their outstanding lactation traits. Lactation is a complex dynamic process regulated by multiple factors. The extensive application of transcriptome sequencing technology enables researchers to further explore this biological issue. This study selected three pregnant and three non-pregnant Kazakh mares as the research subject. Their mammary glands were taken for transcriptome sequencing. The results show that there are 9 lncRNAs and 122 mRNAs differentially expressed between the two groups. GO enrichment analysis shows that there are 175 molecular functions, 59 cellular components, and 555 biological processes, including cellular hormone metabolic process, hormone catabolic process, and I-kappaB kinase/NF-kappaB signaling. KEGG enrichment analysis exhibits that these differential genes are mainly enriched in the NF-kappa B signaling pathway, steroid hormone biosynthesis, breast cancer, ECM-receptor interaction, and MAPK signaling pathway. WNT4, DPP4, and NFKBIA are key nodes regulating breast activation. Conclusions: Through the comparative analysis of the transcriptome data of mammary tissues of pregnant and non-pregnant mares, relevant differentially expressed genes are screened and analyzed. This study provides valuable fundamental data for investigating candidate genes related to the lactation regulation and mammogenesis of Kazakh horses.

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.

Preliminary Transcriptome Analysis of Long Noncoding RNA in Hypothalamic-Pituitary-Mammary Gland Axis of Dairy Cows under Heat Stress

Heat stress (HS) is directly correlated to mammary gland dysfunction in dairy cows, especially in summer. The hypothalamic-pituitary-mammary gland axis (HPM axis) plays an important role in the regulation of stress response and lactation physiology in heat-stressed dairy cows. The aim of this study was to explore the lncRNA profile, and the competitive endogenous RNA (ceRNA) regulatory network in hypothalamus, pituitary, and mammary gland tissues of heat-stressed and normal dairy cows. We performed RNA sequencing (RNA-seq) to identify differentially expressed (DE) lncRNAs, and the ceRNA regulatory network was established in HPM-axis-related tissues. Our results showed that a total of 13, 702 and 202 DE lncRNAs were identified in hypothalamus, pituitary, and mammary glands, respectively. Of lncRNAs, 8, 209 and 45 were up-regulated, and 5, 493 and 157 lncRNAs were down-regulated. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses indicated that DE lncRNAs target genes that might play a role in hormone synthesis, secretion and action, apoptosis, mitogen-activated protein kinase (MAPK), AMP-activated protein kinase (AMPK), and mechanistic target of rapamycin (mTOR) signaling pathway. Moreover, the ceRNA regulatory network associated with the MAPK signaling pathway in HPM-axis-related tissues contains 3286 lncRNA-mRNA pairs. Furthermore, the ceRNA regulatory network associated with apoptosis, prolactin, AMPK, and mTOR signaling pathway in the mammary gland contains 772 lncRNA-mRNA pairs. Thus, some lncRNAs may be involved in the regulation of stress response and the physiological process of lactation. The changes in lncRNA expression profiles and ceRNAs (lncRNA-miRNA-mRNA) in HPM-axis-related tissues are the key to affect the stress response and lactation physiology of dairy cows under HS, which provide a theoretical basis for the molecular mechanism in the stress response of HPM-axis-related tissues in dairy cows under HS.

Lnc-TRTMFS promotes milk fat synthesis via the miR-132x/RAI14/mTOR pathway in BMECs

As an important index to evaluate the quality of milk, milk fat content directly determines the nutrition and flavor of milk. Recently, growing evidence has suggested that long noncoding RNAs (lncRNAs) play important roles in bovine lactation, but little is known about the roles of lncRNAs in milk fat synthesis, particularly the underlying molecular processes. Therefore, the purpose of this study was to explore the regulatory mechanism of lncRNAs in milk fat synthesis. Based on our previous lncRNA-seq data and bioinformatics analysis, we found that Lnc-TRTMFS (transcripts related to milk fat synthesis) was upregulated in the lactation period compared to the dry period. In this study, we found that knockdown of Lnc-TRTMFS significantly inhibited milk fat synthesis, resulting in a smaller amount of lipid droplets and lower cellular triacylglycerol levels, and significantly decreased the expression of genes related to adipogenesis. In contrast, overexpression of Lnc-TRTMFS significantly promoted milk fat synthesis in bovine mammary epithelial cells (BMECs). In addition, Bibiserv2 analysis showed that Lnc-TRTMFS could act as a molecular sponge for miR-132x, and retinoic acid induced protein 14 (RAI14) was a potential target of miR-132x, which was further confirmed by dual-luciferase reporter assays, quantitative reverse transcription PCR, and western blots. We also found that miR-132x significantly inhibited milk fat synthesis. Finally, rescue experiments showed that Lnc-TRTMFS could weaken the inhibitory effect of miR-132x on milk fat synthesis and rescue the expression of RAI14. Taken together, these results revealed that Lnc-TRTMFS regulated milk fat synthesis in BMECs via the miR-132x/RAI14/mTOR pathway.

Transcriptome Analysis of Goat Mammary Gland Tissue Reveals the Adaptive Strategies and Molecular Mechanisms of Lactation and Involution

To understand how genes precisely regulate lactation physiological activity and the molecular genetic mechanisms underlying mammary gland involution, this study investigated the transcriptome characteristics of goat mammary gland tissues at the late gestation (LG), early lactation (EL), peak lactation (PL), late lactation (LL), dry period (DP), and involution (IN) stages. A total of 13,083 differentially expressed transcripts were identified by mutual comparison of mammary gland tissues at six developmental stages. Genes related to cell growth, apoptosis, immunity, nutrient transport, synthesis, and metabolism make adaptive transcriptional changes to meet the needs of mammary lactation. Notably, platelet derived growth factor receptor beta (PDGFRB) was screened as a hub gene of the mammary gland developmental network, which is highly expressed during the DP and IN. Overexpression of PDGFRB in vitro could slow down the G1/S phase arrest of goat mammary epithelial cell cycle and promote cell proliferation by regulating the PI3K/Akt signaling pathway. In addition, PDGFRB overexpression can also affect the expression of genes related to apoptosis, matrix metalloproteinase family, and vascular development, which is beneficial to the remodeling of mammary gland tissue during involution. These findings provide new insights into the molecular mechanisms involved in lactation and mammary gland involution.

lncRNA MPFAST Promotes Proliferation and Fatty Acid Synthesis of Bovine Mammary Epithelial Cell by Sponging miR-103 Regulating PI3K-AKT Pathway

Long noncoding RNAs (lncRNAs) have an essential role in mammary gland development and lactation. Our earlier study showed that the lncRNA mammary proliferation and fatty acid synthesis-associated transcript (MPFAST) is highly expressed in the Holstein cow mammary gland during the middle lactation period compared to the dry period, which indicates its potential role in lactation. Therefore, gain- and loss-of-function experiments were performed on bovine mammary epithelial cells (BMECs) by cell counting kit 8 (CCK8), 5-ethynyl-2'-deoxyuridine (EdU), real-time quantitative polymerase chain reaction (RT-qPCR), and western blot. The results indicated that MPFAST promoted the viability and proliferation of BMECs. The oil red O staining and cellular triglyceride assay further showed that MPFAST promoted the number of lipid droplets and cellular triglyceride synthesis in BMECs. Bioinformatics analysis showed that MPFAST could act as a molecular sponge for miR-103, and PIK3R1 was a potential target of miR-103, which was further confirmed by the dual-luciferase reporter assay, RT-qPCR, and western blot. The overexpression of MPFAST promoted the expression of PIK3R1 at mRNA and protein levels. It also significantly increased the mRNA relative expression levels of AKT, mTOR, and SREBP1, and the protein relative expression levels of AKT and p-AKT in the PI3K-AKT signaling pathway. In contrast, the inhibition of MPFAST resulted in the downregulation of the PI3K-AKT signaling pathway genes. These results indicated that MPFAST regulates the expression of the genes in the PI3K-AKT signaling pathway through sponging miR-103 and promotes the proliferation and synthesis of fatty acids of BMECs. Our results would provide a new direction for further exploring the regulatory mechanism of lncRNA in the mammary gland.

Tissue-Specific Expression of Circ_015343 and Its Inhibitory Effect on Mammary Epithelial Cells in Sheep

Circular RNAs (circRNAs) are a kind of non-coding RNA that have an important molecular function in mammary gland development and lactation of mammals. In our previous study, circ_015343 was found to be highly expressed in the ovine mammary gland tissue at the peak-lactation period by using RNA sequencing (RNA-seq). In the present study, the authenticity of circ_015343 was confirmed by using reverse transcriptase-polymerase chain reaction (RT-PCR) analysis and Sanger sequencing. The circ_015343 was derived from the complete 10 exons of aminoadipic semialdehyde synthase (AASS), ranging from exon 2 to exon 11 and mainly located in cytoplasm of ovine mammary epithelial cells. The circRNA was found to be expressed in eight ovine tissues, with the highest expression level in the mammary gland and the least expression in Longissimus dorsi muscle. The circ_015343 had a lower level of expression in a sheep breed with higher milk yield and milk fat content. The disturbed circ_015343 increased the viability and proliferation of the ovine mammary epithelial cells. The inhibition of circ_015343 also increased the expression levels of three milk fat synthesis marker genes: acetyl-coenzyme A carboxylase alpha (ACACA), fatty acid-binding protein 4 (FABP4), and sterol regulatory element-binding protein 1 (SREBP1), as well as three proliferation-related genes: cyclin dependent kinase 2 (CDK2), cyclin dependent kinase 4 (CDK4) and proliferating cell nuclear antigen (PCNA), but decreased the expression level of its parent gene AASS. A circRNA-miRNA-mRNA interaction network showed that circ_015343 would bind some microRNAs (miRNAs) to regulate the expression of functional genes related to the development of mammary gland and lactation. This study contributes to a better understanding of the roles of circ_015343 in the mammary gland of sheep.

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.

Comparative Transcriptomic Analysis of Hu Sheep Pituitary Gland Prolificacy at the Follicular and Luteal Phases

The pituitary gland directly regulates the reproduction of domestic animals. Research has increasingly focused on the potential regulatory mechanism of non-coding RNA in pituitary development. Little is known about the differential expression pattern of lncRNAs in Hu sheep, a famous sheep breed with high fecundity, and its role in the pituitary gland between the follicular phase and luteal phase. Herein, to identify the transcriptomic differences of the sheep pituitary gland during the estrus cycle, RNA sequencing (RNA-Seq) was performed. The results showed that 3529 lncRNAs and 16,651 mRNAs were identified in the pituitary gland. Among of them, 144 differentially expressed (DE) lncRNA transcripts and 557 DE mRNA transcripts were screened in the follicular and luteal phases. Moreover, GO and KEGG analyses demonstrated that 39 downregulated and 22 upregulated genes interacted with pituitary functions and reproduction. Lastly, the interaction of the candidate lncRNA XR_001039544.4 and its targeted gene LHB were validated in sheep pituitary cells in vitro. LncRNA XR_001039544.4 and LHB showed high expression levels in the luteal phase in Hu sheep. LncRNA XR_001039544.4 is mainly located in the cytoplasm, as determined by FISH analysis, indicating that XR_001039544.4 might act as competing endogenous RNAs for miRNAs to regulate LHB. LncRNA XR_001039544.4 knockdown significantly inhibited LH secretion and cell proliferation. LncRNA XR_001039544.4 may regulate the secretion of LH in the luteal-phase pituitary gland via affecting cell proliferation. Taken together, these findings provided genome-wide lncRNA- and mRNA-expression profiles for the sheep pituitary gland between the follicular and luteal phases, thereby contributing to the elucidation of the molecular mechanisms of pituitary function.

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 <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.

Metabogenomic analysis to functionally annotate the regulatory role of long non-coding RNAs in the liver of cows with different nutrient partitioning phenotype

Long non-coding RNAs (lncRNAs) hold gene regulatory potential, but require substantial further functional annotation in livestock. Applying two metabogenomic approaches by combining transcriptomic and metabolomic analyses, we aimed to identify lncRNAs with potential regulatory function for divergent nutrient partitioning of lactating crossbred cows and to establish metabogenomic interaction networks comprising metabolites, genes and lncRNAs. Through correlation analysis of lncRNA expression with transcriptomic and metabolomic data, we unraveled lncRNAs that have a putative regulatory role in energy and lipid metabolism, the urea and tricarboxylic acid cycles, and gluconeogenesis. Especially FGF21, which correlated with a plentitude of differentially expressed genes, differentially abundant metabolites, as well as lncRNAs, suggested itself as a key metabolic regulator. Notably, lncRNAs in close physical proximity to coding-genes as well as lncRNAs with natural antisense transcripts appear to perform a fine-tuning function in gene expression involved in metabolic pathways associated with different nutrient partitioning phenotypes.

Whole transcriptome analysis of RNA expression profiles reveals the potential regulating action of long noncoding RNA in lactating cows fed a high concentrate diet

Subacute ruminal acidosis (SARA) is a common metabolic disease in the dairy farming industry which is usually caused by an excessive amount of high concentrate diet. SARA not only threatens animal welfare but also leads to economic losses in the farming industry. The liver plays an important role in the distribution of nutritional substances and metabolism; however, a high concentrate diet can cause hepatic metabolic disorders and liver injury. Recently, noncoding RNA has been considered as a critical regulator of hepatic disease, however, its role in the bovine liver is limited. In this study, 12 mid-lactating dairy cows were randomly assigned to a control (CON) group (40% concentrate of dry matter, n = 6) and a SARA group (60% concentrate of dry matter, n = 6). After 21 d of treatment, all cows were sacrificed, and liver tissue samples were collected. Three dairy cows were randomly selected from the CON and SARA groups respectively to perform whole transcriptome analysis. More than 20,000 messenger RNA (mRNA), 10,000 long noncoding RNA (lncRNA), 3,500 circular RNA (circRNA) and 1,000 micro RNA (miRNA) were identified. Furthermore, 43 mRNA, 121 lncRNA and 3 miRNA were differentially expressed, whereas no obvious differentially expressed circRNA were detected between the 2 groups. Gene Ontology (GO) annotation revealed that the differentially expressed genes were mainly enriched in oxidoreductase activity, stress, metabolism, the immune response, cell apoptosis, and cell proliferation. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis showed that the deferentially expressed genes were highly enriched in the phosphatidylinositol 3 kinase (PI3K)-serine/threonine kinase (AKT) signaling pathway (P < 0.05). According to KEGG pathway analysis, the differentially expressed lncRNA (DElncRNA) target genes were mainly related to proteasomes, peroxisomes, and the hypoxia-inducible factor-1 signaling pathway (P < 0.005). Further bioinformatics and integrative analyses revealed that the lncRNA were strongly correlated with mRNA; therefore, it is reasonable to speculate that lncRNA potentially play important roles in the liver dysfunction induced by SARA. Our study provides a valuable resource for future investigations on the mechanisms of SARA to facilitate an understanding of the importance of lncRNA, and offer functional RNA information.

Identification of Regulatory Functions of LncRNAs Associated With T. circumcincta Infection in Adult Sheep

Several recent studies have demonstrated the role of long non-coding RNAs (lncRNAs) in regulating the defense mechanism against parasite infections, but no studies are available that investigated their relevance for immune response to nematode infection in sheep. Thus, the aim of the current study was to (i) detect putative lncRNAs that are expressed in the abomasal lymph node of adult sheep after an experimental infection with the gastrointestinal nematode (GIN) Teladorsagia circumcincta and (ii) to elucidate their potential functional role associated with the differential host immune response. We hypothesized that putative lncRNAs differentially expressed (DE) between samples from animals that differ in resistance to infection may play a significant regulatory role in response to nematode infection in adult sheep. To obtain further support for our hypothesis, we performed co-expression and functional gene enrichment analyses with the differentially expressed lncRNAs (DE lncRNAs). In a conservative approach, we included for this predictive analysis only those lncRNAs that are confirmed and supported by documentation of expression in gastrointestinal tissues in the current sheep gene atlas. We identified 9,105 putative lncRNA transcripts corresponding to 7,124 gene loci. Of these, 457 were differentially expressed lncRNA loci (DELs) with 683 lncRNA transcripts. Based on a gene co-expression analysis via weighted gene co-expression network analysis, 12 gene network modules (GNMs) were found significantly correlated with at least one of 10 selected target DE lncRNAs. Based on the principle of “guilt-by-association,” the DE genes from each of the three most significantly correlated GNMs were subjected to a gene enrichment analysis. The significant pathways associated with DE lncRNAs included ERK5 Signaling, SAPK/JNK Signaling, RhoGDI Signaling, EIF2 Signaling, Regulation of eIF4 and p70S6K Signaling and Oxidative Phosphorylation pathways. They belong to signaling pathway categories like Cellular Growth, Proliferation and Development, Cellular Stress and Injury, Intracellular and Second Messenger Signaling and Apoptosis. Overall, this lncRNA study conducted in adult sheep after GIN infection provided first insights into the potential functional role of lncRNAs in the differential host response to nematode infection.

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.

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.

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.

lncRNA MT1JP Suppresses Biological Activities of Breast Cancer Cells in vitro and in vivo by Regulating the miRNA-214/RUNX3 Axis

Introduction

The purpose of our research was to evaluate MT1JP in breast cancer.

Material and Methods

For clinical purpose, tissues were collected, and a correlation analysis ofMT1JP and miRNA-214 gene expressions was conducted. Using an in vitro study, MDA-MB-231 and MCF-7 cell lines were used as research objects in our research. Colony, flow cytometry, TUNEL, transwell, adhesion and wound healing assay were used to discuss the biological activities of the cells. In an in vivo study, tumor weight and volume were measured, and cell apoptosis was measured by TUNEL assay. The relative mechanism's proteins were evaluated by Western blotting or immunohistochemistry assay.

Results

Compared with adjacent tissues, MT1JP and miRNA-214 gene expressions were significantly different (P<0.001, respectively). By in vitro and in vivo studies, the biological activities of the cells were significantly decreased in MDA-MB-231 and MCF-7 cell lines with MT1JP overexpression. The relative mechanism was correlated with miRNA-214/RUNX3 axis.

Conclusion

The overexpression of MT1JP suppresses the biological activities of breast cancer cells by regulation miRNA-214/RUNX3 axis in vitro and vivo study.

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.

The Downregulation of lncRNA pgm5-as1 Inhibits the Proliferation and Metastasis Via Increasing miR-484 Expression in Colorectal Cancer

Background: Bioinformatics showed that long non-coding RNA (lncRNA) pgm5-as1 was regulated in patients with colorectal cancer (CRC), and miR-484 was also regulated in CRC. We aimed at determining the modulatory pathway of lncRNA pgm5-as1 in CRC cells, and whether miR-484 was involved in the pathway. Materials and Methods: The target gene of pgm5-as1 was predicted by bioinformatics and verified by dual luciferase assay. Transcription levels of pgm5-as1 and miR-484 were determined by quantitative real-time polymerase chain reaction. Viability, migration rate, invasion, and growth of SW480 and HCT116 cells were determined by Cell Counting Kit-8 (CCK-8), wound healing assay, transwell, and colony formation assay, respectively. Results: pgm5-as1 was upregulated in CRC tissues and cell lines; however, its downregulation contributed to the decreasing of cell viability, growth, migration, and invasion of SW480 and HCT116 cells. Moreover, miR-484 was predicted as the target of pgm5-as1, and the downregulation of pgm5-as1 partially restored the elevated cell viability, growth, migration, and invasion that were induced by the inhibition of miR-484 expression in SW480 and HCT116 cells. Conclusions: The loss of miR-484 expression in CRC might be involved in the promotion and metastasis of CRC, which may be caused by the overexpression of pgm5-as1. Hence, the downregulation of pgm5-as1 could be a therapeutic target in the prevention or intervention of CRC.

Genome‑wide integrated analysis demonstrates widespread functions of lncRNAs in mammary gland development and lactation in dairy goats

Background

The mammary gland is a unique organ for milk synthesis, secretion and storage, and it undergoes cyclical processes of development, differentiation, lactation and degeneration. At different developmental periods, the biological processes governing mammary gland physiology and internal environmental homeostasis depend on a complex network of genes and regulatory factors. Emerging evidence indicates that lncRNAs have arbitrarily critical functions in regulating gene expression in many organisms; however, the systematic characteristics, expression, and regulatory roles of lncRNAs in the mammary gland tissues of dairy goats have not been determined.

Result

In the present study, we profiled long noncoding RNA (lncRNA) expression in the mammary gland tissues of Laoshan dairy goats (Capra hircus) from different lactation periods at the whole-genome level, to identify, characterize and explore the regulatory functions of lncRNAs. A total of 37,249 transcripts were obtained, of which 2381 lncRNAs and 37,249 mRNAs were identified, 22,488 transcripts, including 800 noncoding transcripts and 21,688 coding transcripts, differed significantly (p ≤ 0.01) among the different lactation stages. The results of lncRNA-RNA interaction analysis showed that six known lncRNAs belonging to four families were identified as the precursors of 67 known microRNAs; 1478 and 573 mRNAs were predicted as hypothetical cis-regulation elements and antisense mRNAs, respectively. GO annotation and KEGG analysis indicated that the coexpressed mRNAs were largely enriched in biological processes related to such activities as metabolism, immune activation, and stress,., and most genes were involved in pathways related to such phenomena as inflammation, cancer, signal transduction, and metabolism.

Conclusions

Our results clearly indicated that lncRNAs involved in responses to stimuli, multiorganism processes, development, reproductive processes and growth, are closely related to mammary gland development and lactation.

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.

Genome-wide analysis of mRNAs and lncRNAs in Mycoplasma bovis infected and non-infected bovine mammary gland tissues

Mycoplasma bovis (M. bovis) causes diseases such as arthritis, pneumonia, abortion, and mastitis, leading to great losses in the bovine dairy industries. RNA types such as messenger RNAs (mRNAs) and long non-coding RNAs (lncRNAs) play significant roles in regulating the immune responses triggered by bacteria. The expression profiles of mRNA and lncRNA as they occur in bovine mammary gland tissues infected with M. bovis are still not well understood. To illuminate this issue, transcription analysis of mRNA and LncRNAs were conducted on the mammary gland tissues belonging to Holstein cattle infected and not infected with M. bovis. The analysis revealed 1310 differentially expressed mRNAs and 57 differentially expressed lncRNAs in the bovine mammary gland tissues infected and not infected with M. bovis. In addition, 392 novel lncRNAs were detected, 19 of which were differentially expressed. Gene ontology analysis reveals that differentially expressed mRNAs and lncRNAs play significant roles in such vital biological pathways as metabolic pathways, T-cell receptor signaling, TGF-beta signaling, pathways in cancer, PI3K-Akt signaling, NF-kappa B signaling, mTOR signaling, and apoptosis, including in the immune response to cancer. Based on our literature review, this study is the first genome-wide lncRNA research conducted on bovine mammary gland tissues infected with M. bovis. Our results provide bovine mammary gland lncRNA and mRNA resources to understand their roles in the regulation of the immune response against the agent M. bovis in bovine mammary gland tissues.

The Selection of Reference Genes for Quantitative Real-Time PCR in the Ashidan Yak Mammary Gland During Lactation and Dry Period

Simple Summary

The Ashidan yak is a new cultivated breed which has polled characteristics and a mild temperament. Improving milk yield is an important aspect of a breeding program for this breed. The mammary gland undergoes dramatic physiological and metabolic changes during the transition from lactation to dry periods, which involves the expression and regulation of a great number of genes. Quantification of gene expression levels by real-time quantitative polymerase chain reaction (RT-qPCR) is important to reveal the molecular mechanisms of mammary gland development and lactation. The accuracy of RT-qPCR is strongly influenced by the expression stability of reference genes, however, a systematic approach for selecting reference genes used for analyzing gene expression of the Ashidan yak has not been developed. In this study, we selected reference genes and analyzed their expression stability at different physiological stages (lactation and dry period). We found the hydroxymethylbilane synthase gene (HMBS) and the tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein, zeta polypeptide gene (YWHAZ) were the most stable genes of the mammary gland of the Ashidan yak. These results help to improve the accuracy of gene expression analysis and provide a basis for future functional studies of target gene expression in the mammary gland of the Ashidan yak.

Abstract

Investigating the critical genes related to milk synthesis is essential for the improvement of the milk yield of the yak. Real-time quantitative polymerase chain reaction (RT-qPCR) is a reliable and widely used method to measure and evaluate gene expression levels. Selection of suitable reference genes is mandatory to acquire accurate normalization of gene expression results from RT-qPCR. To select the most stable reference genes for reliable normalization of mRNA expression by RT-qPCR in the mammary gland of the Ashidan yak, we selected 16 candidate reference genes and analyzed their expression stability at different physiological stages (lactation and dry period). The expression stability of the candidate reference genes was assessed using geNorm, NormFinder, BestKeeper, Delta Ct, and RefFinder methods. The results showed that the hydroxymethylbilane synthase gene (HMBS) and the tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein, zeta polypeptide gene (YWHAZ) were the most stable genes across all treatment samples. The reliability of selected reference genes was validated by normalizing relative expression of the lactation-related 60S ribosomal protein L35 gene (RPL35). The relative expression of RPL35 varied considerably according to the different reference genes. This work provides valuable information to further promote research in the molecular mechanisms involved in lactation and mammary gland development and provides a foundation for the improvement of the milk yield and quality of the Ashidan yak.

Effect of Lactoferrin on the Expression Profiles of Long Non-coding RNA during Osteogenic Differentiation of Bone Marrow Mesenchymal Stem Cells

Lactoferrin (LF) has demonstrated stimulation of osteogenic differentiation of mesenchymal stem cells (MSCs). Long non-coding RNAs (lncRNAs) participate in regulating the osteogenic differentiation processes. However, the impact of LF on lncRNA expression in MSC osteogenic differentiation is poorly understood. Our aim was to investigate the effects of LF on lncRNAs expression profiles, during osteogenic differentiation of rat bone marrow mesenchymal stem cells (rBMSCs), by RNA sequencing. A total number of 1331 putative lncRNAs were identified in rBMSCs during osteogenic differentiation in the study. LF influenced the expression of 120 lncRNAs (differentially expressed lncRNAs [DELs], Fold change > 1.5 or < −1.5; p < 0.05) in rBMSCs on day 14 of osteogenic differentiation, consisted of 60 upregulated and 60 down-regulated. Furthermore, the potential functions of DELs were of prediction by searching their target cis- and trans-regulated protein-coding genes. The bioinformatic analysis of DELs target gene revealed that LF led to the disfunction of transforming growth factor beta stimulus (TGF-β) and positive regulation of I-κappa B kinase/NF-κappa B signaling pathway, which may relate to osteogenic differentiation of rBMSCs. Our work is the first profiling of lncRNA in osteogenic differentiation of rBMSCs induced by LF, and provides valuable insights into the potential mechanisms for LF promoting osteogenic activity.

A novel mutation in the promoter region of RPL8 regulates milk fat traits in dairy cattle by binding transcription factor Pax6

Ribosomal protein L8 (RPL8) was considered as a promising candidate gene for the milk fat percentage trait in dairy cattle in our previous genome-wide association studies, but the mechanism remains to be determined. Here we investigated the molecular mechanism underlying the effect of bovine RPL8 on milk fat percentage. We demonstrated that RPL8 silencing in bovine mammary epithelial cells affected the expression of genes encoding fat-related enzymes (ACACA, FASN, ACSS1, FABP3, SREBP-1, DGAT1, GPAM, PLIN2, PLIN5 and CIDEA). Furthermore, we showed here that a single nucleotide polymorphism, g.-931G > T (chr14:1508300, UMD3.1) in the putative RPL8 promoter region significantly reduced its promoter activity. Interestingly, this decrease in activity was paralleled by lower RPL8 expression in mammary gland tissues of dairy cattle with the homozygous TT genotype compared to that of cattle with the wild-type homozygous GG genotype. Importantly, we found g.-931G > T added a paired box 6 (Pax6)-binding site and this mutation located in the presumed Pax6-binding site. EMSA and co-immunoprecipitation (Co-IP) assays confirmed the interaction between RPL8 and Pax6 and the T allele exhibited a higher affinity of DNA/protein interactions than G allele, suggesting that Pax6 is an important transcription factor for RPL8 expression. In addition, lactating cows with the GG and GT genotypes presented a significant decrease in milk fat percentage compared to cows with TT genotypes. Altogether, our study indicated that g.-931G > T at RPL8 promoter altered its expression by affecting the interplay between Pax6 and RPL8, which may account for the association with milk fat traits. Findings herein first elucidated the biological function of RPL8 gene in milk fat and the identified SNP g.-931G > T may be considered as genetic makers for breeding in dairy cattle.

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.