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

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.

miR-30a-3p Regulates Autophagy in the Involution of Mice Mammary Glands

The mammary gland undergoes intensive remodeling during the lactation cycle, and the involution process of mammary gland contains extensive epithelial cells involved in the process of autophagy. Our studies of mice mammary glands suggest that miR-30a-3p expression was low during involution compared with its high expression in the mammary glands of lactating mice. Then, we revealed that miR-30a-3p negatively regulated autophagy by autophagy related 12 (Atg12) in mouse mammary gland epithelial cells (MMECs). Restoring ATG12, knocking down autophagy related 5 (Atg5), starvation, and Rapamycin were used to further confirm this conclusion. Overexpression of miR-30a-3p inhibited autophagy and altered mammary structure in the involution of the mammary glands of mice, which was indicative of alteration in mammary remodeling. Taken together, these results elucidated the molecular mechanisms of miR-30a-3p as a key induction mediator of autophagy by targeting Atg12 within the transition period between lactation and involution in mammary glands.

miR-145 Modulates Fatty Acid Metabolism by Targeting FOXO1 to Affect SERBP1 Activity in Bovine Mammary Epithelial Cells

MicroRNA-mediated gene regulation is important for the regulation of fatty acid metabolism and synthesis. Our previous study uncovered that the miR-145 expression is higher in the lactating mammary gland of dairy cows than in the dry-period, but the underlying molecular mechanism is incompletely understood. In this study, we have investigated the potential role of miR-145 in bovine mammary epithelial cells (BMECs). We found that the expression of miR-145 gradually increased during lactation. CRISPR/Cas9-mediated knockout (KO) of miR-145 in BMECs results in the downregulated expression of fatty acid metabolism-associated genes. Further results revealed that miR-145 KO reduced total triacylglycerol (TAG) and cholesterol (TC) accumulation and altered the composition of intracellular fatty acids (C16:0, C18:0, and C18:1). Conversely, miR-145 overexpression had the opposite effect. Bioinformatics online program predicted that miR-145 targets the 3'-UTR of the Forkhead Box O1 (FOXO1) gene. Subsequently, FOXO1 was identified as a direct target of miR-145 by qRT-PCR, Western blot analysis, and luciferase reporter assay. Furthermore, siRNA-mediated silencing of FOXO1 promoted fatty acid metabolism and TAG synthesis in BMECs. Additionally, we observed the involvement of FOXO1 in the transcriptional activity of the sterol regulatory element-binding protein 1 (SREBP1) gene promoter. Overall, our findings indicated that miR-145 relieves the inhibitory effect of FOXO1 on SREBP1 expression by targeting FOXO1 and subsequently regulating fatty acid metabolism. Thus, our results provide valuable information on the molecular mechanisms for improving milk yield and quality from the perspective of miRNA-mRNA networks.

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.

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.

Genome-Wide Association Study for Udder Conformation Traits in Chinese Holstein Cattle

Udder conformation traits are one of the most economic traits in dairy cows, greatly affecting animal health, milk production, and producer profitability in the dairy industry. Genetic analysis of udder structure and scores have been developed in Holstein cattle. In our research, we conducted a genome-wide association study for five udder traits, including anterior udder attachment (AUA), central suspensory ligament (CSL), posterior udder attachment height (PUAH), posterior udder attachment width (PUAW), and udder depth (UD), in which the fixed and random model circulating probability unification (FarmCPU) model was applied for the association analysis. The heritability and the standard errors of these five udder traits ranged from 0.04 ± 0.00 to 0.49 ± 0.03. Phenotype data were measured from 1000 Holstein cows, and the GeneSeek Genomic Profiler (GGP) Bovine 100 K SNP chip was used to analyze genotypic data in Holstein cattle. For GWAS analysis, 984 individual cows and 84,407 single-nucleotide polymorphisms (SNPs) remained after quality control; a total of 18 SNPs were found at the GW significant threshold (p < 5.90 × 10−7). Many candidate genes were identified within 200kb upstream or downstream of the significant SNPs, which include MGST1, MGST2, MTUS1, PRKN, STXBP6, GRID2, E2F8, CDH11, FOXP1, SLF1, TMEM117, SBF2, GC, ADGRB3, and GCLC. Pathway analysis revealed that 58 Gene Ontology (GO) terms and 18 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were enriched with adjusted p values, and these GO terms and the KEGG pathway analysis were associated with biological information, metabolism, hormonal growth, and development processes. These results could give valuable biological information for the genetic architecture of udder conformation traits in dairy Holstein cattle.

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

Genome-Wide Association Study Candidate Genes on Mammary System-Related Teat-Shape Conformation Traits in Chinese Holstein Cattle

In the dairy industry, mammary system traits are economically important for dairy animals, and it is important to explain their fundamental genetic architecture in Holstein cattle. Good and stable mammary system-related teat traits are essential for producer profitability in animal fitness and in the safety of dairy production. In this study, we conducted a genome-wide association study on three traits—anterior teat position (ATP), posterior teat position (PTP), and front teat length (FTL)—in which the FarmCPU method was used for association analyses. Phenotypic data were collected from 1000 Chinese Holstein cattle, and the GeneSeek Genomic Profiler Bovine 100K single-nucleotide polymorphisms (SNP) chip was used for cattle genotyping data. After the quality control process, 984 individual cattle and 84,406 SNPs remained for GWAS work analysis. Nine SNPs were detected significantly associated with mammary-system-related teat traits after a Bonferroni correction (p < 5.92 × 10⁻⁷), and genes within a region of 200 kb upstream or downstream of these SNPs were performed bioinformatics analysis. A total of 36 gene ontology (GO) terms and 3 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were significantly enriched (p < 0.05), and these terms and pathways are mainly related to metabolic processes, immune response, and cellular and amino acid catabolic processes. Eleven genes including MMS22L, E2F8, CSRP3, CDH11, PEX26, HAL, TAMM41, HIVEP3, SBF2, MYO16 and STXBP6 were selected as candidate genes that might play roles in the teat traits of cows. These results identify SNPs and candidate genes that give helpful biological information for the genetic architecture of these teat traits, thus contributing to the dairy production, health, and genetic selection of Chinese Holstein cattle.

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.

Progress on the Regulation of Ruminant Milk Fat by Noncoding RNAs and ceRNAs

Milk fat is not only a key factor affecting the quality of fresh milk but also a major target trait forbreeding. The regulation of milk fat involves multiple genes, network regulation and signal transduction. To explore recent discoveries of pathway regulation, we reviewed the published literature with a focus on functional noncoding RNAs and epigenetic regulation in ruminants. Results indicate that miRNAs play key roles in the regulation of milk fat synthesis and catabolism in ruminants. Although few data are available, merging evidence indicates that lncRNAs and circRNAs act on milk fat related genes through indirect action with microRNAs or RNAs in the ceRNA network to elicit positive effects on transcription. Although precise regulatory mechanisms remain unclear, most studies have focused on the regulation of the function of target genes through functional noncoding RNAs. Data to help identify factors that can regulate their own expression and function or to determine whether self-regulation involves positive and/or negative feedback are needed. Despite the growing body of research on the role of functional noncoding RNA in the control of ruminant milk fat, most data are still not translatable for field applications. Overall, the understanding of mechanisms whereby miRNA, lncRNA, circRNA, and ceRNA regulate ruminant milk fat remains an exciting area of research.

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.

Transcriptomic in silico analysis of bovine Escherichia coli mastitis highlights its immune-related expressed genes as an effective biomarker

Background

Mastitis is one of the major diseases causing economic loss to the dairy industry by reducing the quantity and quality of milk. Thus, the objective of this scientific study was to find new biomarkers based on genes for the early prediction before its severity.

Methods

In the present study, advanced bioinformatics including hierarchical clustering, enrichment analysis, active site prediction, epigenetic analysis, functional domain identification, and protein docking were used to analyze the important genes that could be utilized as biomarkers and therapeutic targets for mastitis.

Results

Four differentially expressed genes (DEGs) were identified in different regions of the mammary gland (teat cistern, gland cistern, lobuloalveolar, and Furstenberg’s rosette) that resulted in 453, 597, 577, and 636 DEG, respectively. Also, 101 overlapped genes were found by comparing 27 different expressed genes. These genes were associated with eight immune response pathways including NOD-like receptor signaling pathway (IL8, IL18, IL1B, PYDC1) and chemokine signaling pathway (PTK2, IL8, NCF1, CCR1, HCK). Meanwhile, 241 protein-protein interaction networks were developed among overlapped genes. Fifty-seven regulatory events were found between miRNAs, expressed genes, and the transcription factors (TFs) through micro-RNA and transcription factors (miRNA-DEG-TF) regulatory network. The 3D structure docking model of the expressed genes proteins identified their active sites and the binding ligands that could help in choosing the appropriate feed or treatment for affected animals.

Conclusions

The novelty of the distinguished DEG and their pathways in this study is that they can precisely improve the detection biomarkers and treatments techniques of cows’ Escherichia coli mastitis disease due to their high affinity with the target site of the mammary gland before appearing the symptoms.

Graphical abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s43141-021-00235-x.

Bta-miR-34b controls milk fat biosynthesis via the Akt/mTOR signaling pathway by targeting RAI14 in bovine mammary epithelial cells

BACKGROUND

The biosynthesis of milk fat affects both the technological properties and organoleptic quality of milk and dairy products. MicroRNAs (miRNAs) are endogenous small non-coding RNAs that inhibit the expression of their mRNA targets and are involved in downstream signaling pathways that control several biological processes, including milk fat synthesis. miR-34b is a member of the miR-34 miRNA cluster, which is differentially expressed in the mammary gland tissue of dairy cows during lactation and dry periods. Previous studies have indicated miR-34b is a potential candidate gene that plays a decisive role in regulating milk fat synthesis; therefore, it is important to focus on miR-34b and investigate its regulatory effect on the biosynthesis of milk fat in bovine mammary epithelial cells (BMECs).

RESULTS

In this study, elevated miR-34b levels reduced milk fat synthesis, upregulated 1,999 genes, and downregulated 2,009 genes in BMECs. Moreover, Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of differentially expressed genes suggested that miR-34b may play an inhibitory role in milk fat synthesis via the protein kinase B (Akt)/mammalian target of rapamycin (mTOR) signaling pathway by reducing phosphorylation levels. Notably, the mTOR activator MHY1485 rescued the inhibitory effect of miR-34b. Furthermore, we demonstrated that retinoic acid-induced protein 14 (RAI14) is a target of miR-34b via TargetScan and immunofluorescence assays. RAI14 mRNA and protein levels were significantly decreased by the miR-34b mimic and increased by the miR-34b inhibitor. Moreover, the reduction in RAI14 levels led to the inhibition of the Akt/mTOR signaling pathway.

CONCLUSIONS

Overall, our results identified a miR-34b-RAI14-Akt/mTOR regulatory network, while also providing a theoretical basis for the molecular breeding of dairy cows.

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.

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.

Comprehensive MicroRNA Expression Profile of the Mammary Gland in Lactating Dairy Cows With Extremely Different Milk Protein and Fat Percentages

A total of 31 differentially expressed genes in the mammary glands were identified in our previous study using RNA sequencing (RNA-Seq), for lactating cows with extremely high and low milk protein and fat percentages. To determine the regulation of milk composition traits, we herein investigated the expression profiles of microRNA (miRNA) using small RNA sequencing based on the same samples as in the previous RNA-Seq experiment. A total of 497 known miRNAs (miRBase, release 22.1) and 49 novel miRNAs among the reads were identified. Among these miRNAs, 71 were found differentially expressed between the high and low groups (p < 0.05, q < 0.05). Furthermore, 21 of the differentially expressed genes reported in our previous RNA-Seq study were predicted as target genes for some of the 71 miRNAs. Gene ontology and KEGG pathway analyses showed that these targets were enriched for functions such as metabolism of protein and fat, and development of mammary gland, which indicating the critical role of these miRNAs in regulating the formation of milk protein and fat. With dual luciferase report assay, we further validated the regulatory role of 7 differentially expressed miRNAs through interaction with the specific sequences in 3'UTR of the targets. In conclusion, the current study investigated the complexity of the mammary gland transcriptome in dairy cattle using small RNA-seq. Comprehensive analysis of differential miRNAs expression and the data from previous study RNA-seq provided the opportunity to identify the key candidate genes for milk composition traits.

Circular RNAs Are Regulators of Diverse Animal Transcriptomes: One Health Perspective

Derived from linear (parental) precursor mRNA, circRNA are recycled exons and introns whose ends are ligated. By titrating microRNAs and RNA binding proteins, circRNA interconnect networks of competing endogenous RNAs. Without altering chromosomal DNA, circRNA regulates skeletal muscle development and proliferation, lactation, ovulation, brain development, and responses to infections and metabolic stress. This review integrates emerging knowledge of circRNA activity coming from genome-wide characterizations in many clades of animals. circRNA research addresses one of the main pillars of the One Health vision – to improve the health and productivity of food animals and generate translational knowledge in animal species.

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.

Identification of Gene Modules and Hub Genes Involved in Mastitis Development Using a Systems Biology Approach

Objective

Mastitis is defined as the inflammation of the mammary gland, which impact directly on the production performance and welfare of dairy cattle. Since, mastitis is a multifactorial complex disease and the molecular pathways underlying this disorder have not been clearly understood yet, a system biology approach was used in this study to a better understanding of the molecular mechanisms behind mastitis.

Methods

Publicly available RNA-Seq data containing samples from milk of five infected and five healthy Holstein cows at five time points were retrieved. Gene Co-expression network analysis (WGCNA) approach and functional enrichment analysis were then applied with the aim to find the non-preserved module of genes that their connectivity were altered under infected condition. Hub genes were identified in the non-preserved modules and were subjected to protein-protein interactions (PPI) network construction.

Results

Among the 25 modules identified, eight modules were non-preserved and were also biologically associated with inflammation, immune response and mastitis development. Interestingly most of the hub genes in the eight modules were also densely connected in the PPI network. Of the hub genes, 250 genes were hubs in both co-expression and PPI networks and most of them were reported to play important roles in immune response or inflammatory pathways. The blue module was highly enriched in inflammatory responses and STAT1 was suggested to play an important role in mastitis development by regulating the immune related genes in this module. Moreover, a set of highly connected genes were identified such as BIRC3, PSMA6, FYN, F11R, NFKBIZ, NFKBIA, GRO1, PHB, CD3E, IL16, GSN, SOCS2, HCK, VAV1 and TLR6, which have been established to be critical for mastitis pathogenesis.

Conclusion

This study improved the understanding of the mechanisms underlying bovine mastitis and suggested eight non-preserved modules along with several most important genes with promising potential in etiology of mastitis.

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.

Long noncoding RNA and mRNA profiling of hypothalamic-pituitary-mammary gland axis in lactating sows under heat stress

Heat stress (HS) seriously affects sow lactation performance and Long non-coding RNAs (lncRNAs) play vital roles in the regulation of transcription and post transcription. However, the mechanism of lncRNAs expression affecting lactation performance on the hypothalamus-pituitary-mammary axis of sows is still unclear. In this study, we performed RNA sequencing and bioinformatics analysis of the hypothalamus, pituitary, and mammary gland tissues of lactating sows under HS and thermal comfort. In total, the analysis identified 658, 6021, and 6745 differently expressed (DE) mRNAs, 26, 126, and 169 DE lncRNAs between comparison groups in the hypothalamus, pituitary, and mammary glands, respectively. The hormone genes and most DE mRNAs encoding heat shock protein were differently expressed in the HS group. In addition, 2, 60, and 86 pairs of DE lncRNAs and mRNAs correlation were observed in those tissues, respectively. Some lncRNAs may be involved in the regulation of lactation performance in the HS sows.

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.

Analyzing the genomic and transcriptomic architecture of milk traits in Murciano-Granadina goats

Background

In this study, we aimed to investigate the molecular basis of lactation as well as to identify the genetic factors that influence milk yield and composition in goats. To achieve these two goals, we have analyzed how the mRNA profile of the mammary gland changes in seven Murciano-Granadina goats at each of three different time points, i.e. 78 d (T1, early lactation), 216 d (T2, late lactation) and 285 d (T3, dry period) after parturition. Moreover, we have performed a genome-wide association study (GWAS) for seven dairy traits recorded in the 1st lactation of 822 Murciano-Granadina goats.

Results

The expression profiles of the mammary gland in the early (T1) and late (T2) lactation were quite similar (42 differentially expressed genes), while strong transcriptomic differences (more than one thousand differentially expressed genes) were observed between the lactating (T1/T2) and non-lactating (T3) mammary glands. A large number of differentially expressed genes were involved in pathways related with the biosynthesis of amino acids, cholesterol, triglycerides and steroids as well as with glycerophospholipid metabolism, adipocytokine signaling, lipid binding, regulation of ion transmembrane transport, calcium ion binding, metalloendopeptidase activity and complement and coagulation cascades. With regard to the second goal of the study, the performance of the GWAS allowed us to detect 24 quantitative trait loci (QTLs), including three genome-wide significant associations: QTL1 (chromosome 2, 130.72-131.01 Mb) for lactose percentage, QTL6 (chromosome 6, 78.90-93.48 Mb) for protein percentage and QTL17 (chromosome 17, 11.20 Mb) for both protein and dry matter percentages. Interestingly, QTL6 shows positional coincidence with the casein genes, which encode 80% of milk proteins.

Conclusions

The abrogation of lactation involves dramatic changes in the expression of genes participating in a broad array of physiological processes such as protein, lipid and carbohydrate metabolism, calcium homeostasis, cell death and tissue remodeling, as well as immunity. We also conclude that genetic variation at the casein genes has a major impact on the milk protein content of Murciano-Granadina goats.

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

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

Bta-miR-34b regulates milk fat biosynthesis by targeting mRNA decapping enzyme 1A (DCP1A) in cultured bovine mammary epithelial cells1

Milk fat is a main nutritional component of milk, and it has become one of the important traits of dairy cow breeding. Recently, there is increasing evidence that microRNAs (miRNA) play significant roles in the process of milk fat synthesis in the mammary gland. Primary bovine mammary epithelial cells (BMEC) were harvested from midlactation cows and cultured in DMEM/F-12 medium with 10% fetal bovine serum, 100 units/mL penicillin, 100 µg/mL streptomycin, 5 µg/mL bovine insulin, 1 µg/mL hydrocortisone, and 2 µg/mL bovine prolactin. We found that miR-34b mimic transfection in BMEC reduced the content of intracellular triacylglycerol (TAG) and lipid droplet accumulation via triacylglycerol assay and Oil Red O staining; meanwhile, overexpression of miR-34b inhibited mRNA expression of lipid metabolism-related genes such as peroxisome proliferator-activated receptor gamma (PPARγ), fatty acid synthase (FASN), fatty acid binding protein 4 (FABP4), and CCAAT enhancer binding protein alpha (C/EBPα). Whereas miR-34b inhibitor resulted in completely opposite results. Furthermore, q-PCR and western blot analysis revealed the mRNA and protein expression levels of DCP1A were downregulated in miR-34b mimic transfection group and upregulated in miR-34b inhibitor group. Moreover, luciferase reporter assays verified that DCP1A was the direct target of miR-34b and DCP1A gene silencing in BMEC-inhibited TAG accumulation and suppressed lipid droplet formation. In conclusion, these findings revealed a novel miR-34b-DCP1A axis that has a significant role in regulating milk fat synthesis and suggested that miR-34b may be used to improve the beneficial ingredients in milk.

MiR-21-3p Centric Regulatory Network in Dairy Cow Mammary Epithelial Cell Proliferation

MicroRNA-mediated gene regulation is important for the development of the mammary gland and the lactating process. A previous study has shown that the expression of microRNA-21 (miR-21) is different in the dry and early lactation period of the dairy cow mammary gland, but the molecular mechanisms underlying the lactation cycle are not fully understood. Here, the function of miR-21-3p on bovine mammary gland epithelial cells (BMECs) was detected by MTT assay and flow cytometry analysis, which showed that miR-21-3p significantly promoted the cell viability and proliferation. Then, the regulating mechanism of miR-21-3p on cell viability and proliferation was elucidated. Dual luciferase assay, RT-qPCR, and Western blot results revealed that IGFBP5 was a target gene of miR-21-3p. It was known that lncRNA could act as a competing endogenous RNA to sequester miRNAs and reduce the regulatory effect of miRNA-targeted genes. Based on our previous lncRNA-seq data and bioinformatics analysis, lncRNA NONBTAT017009.2 was potentially associated with miR-21-3p, and its expression was specifically inhibited with the transfection of miR-21-3p mimic into BMECs. Inversely, the overexpression of NONBTAT017009.2 significantly decreased the expression level of miR-21-3p in BMECs, while the expression of IGFBP5, the target gene of miR-21-3p, was significantly upregulated. In addition, the promoter region of miR-21 contained two STAT3 binding sites, and the dual luciferase reporter assays revealed that the overexpression of STAT3 significantly reduced the promoter activity of miR-21, implying that the transcription factor STAT3 may act as an upstream regulator affecting the regulation process of miR-21-3p. The overexpression of STAT3 significantly inhibited the expression of miR-21-3p, while the mRNA expression of IGFBP5 was significantly increased compared with the control group. Besides, there are no STAT3 binding sites in the promoter region of IGFBP5 as we predicted by gene-regulation and JASPAR software. Therefore, it could infer that STAT3 might regulate the expression of IGFBP5 by miR-21-3p. Taken together, these results established a regulatory network of miR-21-3p to illustrate the regulating mechanism on promoting cow mammary epithelial cell proliferation.

MiR-16a Regulates Milk Fat Metabolism by Targeting Large Tumor Suppressor Kinase 1 (LATS1) in Bovine Mammary Epithelial Cells

Milk contains a number of beneficial fatty acids including short and medium chain and unsaturated conjugated and nonconjugated fatty acids. In this study, microRNA sequencing of mammary tissue collected in early-, peak-, mid-, and late-lactation periods was performed to determine the miRNA expression profiles. miR-16a was one of the differentially expressed miRNA and was selected for in-depth functional studies pertaining to fatty acid metabolism. The mimic of miR-16a impaired fat metabolism [triacylglycerol (TAG) and cholesterol] while knock-down of miR-16a promoted fat metabolism in vitro in bovine mammary epithelial cells (BMECs). In addition, the in vitro work with BMECs also revealed that miR-16a had a negative effect on the cellular concentration of cis 9-C18:1, total C18:1, C20:1, and C22:1 and long-chain polyunsaturated fatty acids. Therefore, these data suggesting a negative effect on fatty acid metabolism extend the discovery of the key role of miR-16a in mediating adipocyte differentiation. Through a combination of bioinformatics analysis, target gene 3' UTR luciferase reporter assays, and western blotting, we identified large tumor suppressor kinase 1 (LATS1) as a target of miR-16a. Transfection of siRNA-LATS1 into BMECs led to increases in TAG, cholesterol, and cellular fatty acid concentrations, suggesting a positive role of LATS1 in mammary cell fatty acid metabolism. In summary, data suggest that miR-16a regulates biological processes associated with intracellular TAG, cholesterol, and unsaturated fatty acid synthesis through LATS1. These data provide a theoretical and experimental framework for further clarifying the regulation of lipid metabolism in mammary cells of dairy cows.

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.

Exploration of long noncoding RNA in bovine milk exosomes and their stability during digestion in vitro

Previous studies have demonstrated that bovine milk contains mRNA and microRNA that are largely encapsulated in milk-derived exosomes. However, little information is available about long noncoding RNAs (lncRNA) in bovine milk. Increasing evidence suggests that lncRNA are of particular interest given their key role in gene expression and development. We performed a comprehensive analysis of lncRNA in bovine milk exosomes by RNA sequencing. We used a validated human in vitro digestion model to investigate the stability of lncRNA encapsulated in bovine milk exosomes during the digestion process. We identified 3,475 novel lncRNA and 6 annotated lncRNA. The lncRNA shared characteristics with those of other mammals in terms of length, exon number, and open reading frames. However, lncRNA showed higher expression than mRNAs. We selected 12 lncRNA of high-expression abundance and identified them by PCR. Gene ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses showed that lncRNA regulate immune function, osteoblastogenesis, neurodevelopment, reproduction, cell proliferation, and cell-cell communication. We also investigated the 12 lncRNA using quantitative real-time PCR to reveal their expression profiles in milk exosomes during different stages of lactation (colostrum 2 d, 30 d, 150 d, and 270 d); their resulting expression levels in milk exosomes showed variations across the stages. A digestion experiment showed that bovine milk exosome lncRNA was resistant to in vitro digestion with different digestive juices, including saliva, gastric juice, pancreatic juice, and bile juice. Taken together, these results show for the first time that cow milk contains lncRNA, and that their abundance varied at different stages of lactation. As expected, bovine milk exosomal lncRNA were stable during in vitro digestion. These findings provide a basis for further understanding of the physiological role of milk lncRNA.

Invited review: Advances and challenges in application of feedomics to improve dairy cow production and health

Dairy cattle science has evolved greatly over the past century, contributing significantly to the improvement in milk production achieved today. However, a new approach is needed to meet the increasing demand for milk production and address the increased concerns about animal health and welfare. It is now easy to collect and access large and complex data sets consisting of molecular, physiological, and metabolic data as well as animal-level data (such as behavior). This provides new opportunities to better understand the mechanisms regulating cow performance. The recently proposed concept of feedomics could help achieve this goal by increasing our understanding of interactions between the different components or levels and their impact on animal production. Feedomics is an emerging field that integrates a range of omics technologies (e.g., genomics, epigenomics, transcriptomics, proteomics, metabolomics, metagenomics, and metatranscriptomics) to provide these insights. In this way, we can identify the best strategies to improve overall animal productivity, product quality, welfare, and health. This approach can help research communities elucidate the complex interactions among nutrition, environment, management, animal genetics, metabolism, physiology, and the symbiotic microbiota. In this review, we summarize the outcomes of the most recent research on omics in dairy cows and highlight how an integrated feedomics approach could be applied in the future to improve dairy cow production and health. Specifically, we focus on 2 topics: (1) improving milk yield and milk quality, and (2) understanding metabolic physiology in transition dairy cows, which are 2 important challenges faced by the dairy industry worldwide.