Genes regulating lipid and protein metabolism are highly expressed in mammary gland of lactating dairy goats

Identification of Key Proteomic Markers for Enhanced Conjugated Linoleic Acid Biosynthesis in Lactating Goats via Linseed Oil Supplementation

This study investigated the effects of linseed oil (LO) supplementation on conjugated linoleic acid (CLA) biosynthesis in lactating Saanen goats. Goats (DIM = 96 ± 14 days) were divided into control (n = 6) and LO treatment (n = 6) groups. LO supplementation was implemented in both short- and long-term experimental settings. Short-term LO supplementation significantly increased trans-11 vaccenic acid and cis-9,trans-11 CLA in milk fat without affecting lactation performance. Using two-dimensional gel electrophoresis, we identified six upregulated proteins in milk somatic cells, including translocon-associated protein (SSRD), succinyl-CoA ligase (SUCB2), ATP synthase subunit (ATPD), stress-70 protein (GRP75), NADH dehydrogenase (NDUFS2), and cytochrome complex QCR1. Long-term LO supplementation enhanced milk fat content and cis-9,trans-11 CLA levels, while significantly elevating the mRNA expression of stearoyl-CoA desaturase (SCD) and all previously identified proteins-including proteasome 20s subunit alpha 5 (PSMA5). These findings extend beyond the known SCD pathway, revealing novel protein markers and potential mechanisms associated with CLA biosynthesis in mammary tissue and milk somatic cells.

Establishment of goat mammary organoid cultures modeling the mammary gland development and lactation

BACKGROUND

Although several cell culture systems have been developed to investigate the function of the mammary gland in dairy livestock, they have potential limitations, such as the loss of alveolar structure or genetic and phenotypic differences from their native counterparts. Overcoming these challenges is crucial for lactation research. Development of protocols to establish lactating organoid of livestock represents a promising goal for the future. In this study, we developed a protocol to establish a culture system for mammary organoids in dairy goats to model the mammary gland development and lactation process.

RESULTS

The organoids cultured within an extracellular matrix gel maintained a bilayer structure that closely resembled the native architecture of mammary tissue. The expansion of mammary organoids was significantly promoted by growth factors containing epidermal growth factor and fibroblast growth factor 2 whereas the proliferative index of the organoids was significantly inhibited by the treatment with WNT inhibitors. Upon stimulation with a lactogenic medium containing prolactin, the mammary organoids exhibited efficient lactation, characterized by the accumulation of lipid droplets in the lumen space. The lactation could be sustained for more than 3 weeks. Importantly, the expression patterns of genes related to fatty acid synthesis and milk proteins in lactating organoids closely mirrored those observed in mammary tissues. These observations were confirmed by data from proteomic analysis that the bulk of milk proteins was produced in the lactating organoids.

CONCLUSION

This study is the first to establish a mammary organoid culture system modeling the mammary gland development and lactation process in ruminants. The efficient induction of lactation in ruminant mammary organoids holds promises for advancing the field of cell-based milk bio-manufacture in the food industry.

Transcriptomic analysis reveals molecular insights into lactation dynamics in Jakhrana goat mammary gland

BACKGROUND

Goat milk is gaining popularity as a superior alternative to bovine milk due to its closer resemblance to human milk. Understanding the molecular processes underlying lactation is crucial for improving milk quality and production in goats. However, the genetic mechanisms governing lactation in goats, particularly in indigenous breeds like the Jakhrana, remain largely unexplored.

RESULTS

In this study, we performed a comprehensive transcriptomic analysis of Jakhrana goat mammary glands during early and late lactation stages. We isolated milk somatic cells and conducted RNA sequencing, followed by transcript quantification and mapping against the ARS1.2 Capra hircus reference assembly. Our analysis identified differentially expressed genes (DEGs) and commonly expressed genes (CEGs) across the lactation phases. Early lactation showed enrichment of genes encoding antimicrobial peptides and lubrication proteins, while late lactation exhibited heightened expression of genes encoding major milk proteins. Additionally, DEG analysis revealed upregulation of pivotal genes, such as the ABC transporter gene MRP4, implicated in modulating milk composition and quality.

CONCLUSION

Our findings provide insights into the genetic mechanisms underlying lactation dynamics in the Jakhrana goat. Understanding these mechanisms could help in improving milk production and quality in goats, benefiting both the dairy industry and consumers.

Prolactin Modulates the Proliferation and Secretion of Goat Mammary Epithelial Cells via Regulating Sodium-Coupled Neutral Amino Acid Transporter 1 and 2

The prolactin (PRL) hormone is a major regulator of mammary gland development and lactation. However, it remains unclear whether and how PRL contributes to mammary epithelial cell proliferation and secretion. The Boer and Macheng black crossbred goats are superior in reproduction, meat, and milk, and are popular in Hubei province. To elucidate the mechanisms of PRL on mammary growth and lactation, to improve the local goat economic trade, we have performed studies on these crossbred goats during pregnancy and early lactation, and in goat mammary epithelial cells (GMECs). Here, we first found that the amino acid transporters of SNAT1 and SNAT2 expression in vivo and in vitro were closely associated with PRL levels, the proliferation and secretion of GMECs; knockdown and over-expression of SNAT1/2 demonstrated that PRL modulated the proliferation and lactation of GMECs through regulating SNAT1/2 expression. Transcriptome sequencing and qPCR assays demonstrated the effect of PRL on the transcriptional regulation of SNAT1 and SNAT2 in GMECs. Dual-luciferase reporter gene assays further verified that the binding of the potential PRL response element in the SNAT1/2 promoter regions activated SNAT1/2 transcription after PRL stimulation. Additionally, silencing of either PRLR or STAT5 nearly abolished PRL-stimulated SNAT1/2 promoter activity, suggesting PRLR-STAT5 signaling is involved in the regulation of PRL on the transcriptional activation of SNAT1/2. These results illustrated that PRL modulates the proliferation and secretion of GMECs via PRLR-STAT5-mediated regulation of the SNAT1/2 pathway. This study provides new insights into how PRL affects ruminant mammary development and lactation through regulation of amino acid transporters.

Integrated analysis of genomics and transcriptomics revealed the genetic basis for goaty flavor formation in goat milk

Goat milk exhibits a robust and distinctive "goaty" flavor. However, the underlying genetic basis of goaty flavor remains elusive and requires further elucidation at the genomic level. Through comparative genomics analysis, we identified divergent signatures of certain proteins in goat, sheep, and cow. MMUT has undergone a goat-specific mutation in the B12 binding domain. We observed the goat FASN exhibits nonsynonymous mutations in the acyltransferase domain. Structural variations in these key proteins may enhance the capacity for synthesizing goaty flavor compounds in goat. Integrated omics analysis revealed the catabolism of branched-chain amino acids contributed to the goat milk flavor. Furthermore, we uncovered a regulatory mechanism in which the transcription factor ZNF281 suppresses the expression of the ECHDC1 gene may play a pivotal role in the accumulation of flavor substances in goat milk. These findings provide insights into the genetic basis underlying the formation of goaty flavor in goat milk. STATEMENT OF SIGNIFICANCE: Branched-chain fatty acids (BCFAs) play a crucial role in generating the distinctive "goaty" flavor of goat milk. Whether there is an underlying genetic basis associated with goaty flavor is unknown. To begin deciphering mechanisms of goat milk flavor development, we collected transcriptomic data from mammary tissue of goat, sheep, cow, and buffalo at peak lactation for cross-species transcriptome analysis and downloaded nine publicly available genomes for comparative genomic analysis. Our data indicate that the catabolic pathway of branched-chain amino acids (BCAAs) is under positive selection in the goat genome, and most genes involved in this pathway exhibit significantly higher expression levels in goat mammary tissue compared to other species, which contributes to the development of flavor in goat milk. Furthermore, we have elucidated the regulatory mechanism by which the transcription factor ZNF281 suppresses ECHDC1 gene expression, thereby exerting an important influence on the accumulation of flavor compounds in goat milk. These findings provide insights into the genetic mechanisms underlying flavor formation in goat milk and suggest further research to manipulate the flavor of animal products.

Unraveling lipid metabolism reprogramming for overcoming drug resistance in melanoma

Cutaneous melanoma is the deadliest form of skin cancer, and its incidence is continuing to increase worldwide in the last decades. Traditional therapies for melanoma can easily cause drug resistance, thus the treatment of melanoma remains a challenge. Various studies have focused on reversing the drug resistance. As tumors grow and progress, cancer cells face a constantly changing microenvironment made up of different nutrients, metabolites, and cell types. Multiple studies have shown that metabolic reprogramming of cancer is not static, but a highly dynamic process. There is a growing interest in exploring the relationship between melanoma andmetabolic reprogramming, one of which may belipid metabolism. This review frames the recent research progresses on lipid metabolism in melanoma.In addition, we emphasize the dynamic ability of metabolism during tumorigenesis as a target for improving response to different therapies and for overcoming drug resistance in melanoma.

The Dietary Inclusion of Ensiled Olive Cake Increases Unsaturated Lipids in Milk and Alters the Expression of Lipogenic Genes in Mammary and Adipose Tissue in Goats

This study aimed to evaluate the effect of the dietary inclusion of ensiled OC on milk yield, composition, fatty acid (FA) profile, and the expression of selected genes involved in lipid metabolism in the udder and adipose tissue of goats. Seventy-two Damascus dairy goats in mid-lactation were assigned randomly to three iso-nitrogenous and iso-energetic diets containing 0, 10, and 20% of ensiled OC as a replacement of forage (OC0, OC10, and OC20, respectively) for 42 days. During weeks 5 and 6 of the trial, dry matter intake, milk yield, milk composition, and FA profiles were recorded, while mammary and perirenal adipose tissue samples were also collected from six animals per treatment from the OC0 and OC20 groups for gene expression analysis. No significant differences were observed among groups concerning milk yield, 4% fat-corrected milk, fat, or protein yield (kg/d). In contrast, the milk fat percentage was gradually increased with increasing OC inclusion rates in the diets, while milk protein percentages were elevated in both OC groups but significantly only in the milk of the OC20 group. The content of FA between C4:0 to C16:0 was reduced, while mono-unsaturated FA (MUFA) concentration was enhanced in the goat milk of OC groups. The OC feeding treatment was associated with the increased mammary expression of SLC2A1 (p < 0.05), VLDLR (p < 0.01), FABP3 (p < 0.01), and elevated SLC2A1 (p < 0.05) and FASN (p < 0.01) gene expression in the adipose tissue of goats fed the OC20 diet. Overall, OC can be used in goats' diets as a forage replacement, at least in the inclusion rate of 20% DM, since this could increase the milk protein and fat percentage and enrich its content with beneficial for human health lipids without adversely affecting milk production traits.

Identification of milk-related genes and regulatory networks in Bactrian camel either supplemented or under grazing

Using gene co-expression networks to understand dynamic characterizations in lactating animals becomes a common method. However, there are rarely reporters focusing on milk traits in Bactrian camel by high-throughput sequencing. We used RNA-seq to generate the camel transcriptome from the blood of 16 lactating Alxa Bactrian camel in different feeding groups. In total, we obtained 1185 milk-related genes correlated with milk yield, milk protein, milk fat, and milk lactose across the WGCNA analysis. Moreover, 364 milk-related genes were differentially expressed between supplementation and grazing feeding groups. The differential expression-camel milk-related genes CMRGs (DE-CMRGs) in supplement direct an intensive gene co-expression network to improve milk performance in lactating camels. This study provides a non-invasive method to identify the camel milk-related genes in camel blood for four primary milk traits and valuable theoretical basis and research ideas for the study of the milk performance regulation mechanism of camelid animals.

Chronic Stress Decreases Lactation Performance

The ability to provision offspring with milk is a significant adaptive feature of mammals that allows for considerable maternal regulation of offspring beyond gestation, as milk provides complete nutrition for developing neonates. For mothers, lactation is a period of marked increases in energetic and nutritive demands to support milk synthesis; because of this considerable increase in demand imposed on multiple physiological systems, lactation is particularly susceptible to the effects of chronic stress. Here, we present work that explores the impact of chronic stress during lactation on maternal lactation performance (i.e., milk quality and quantity) and the expression of key milk synthesis genes in mammary tissue using a Sprague-Dawley rat model. We induced chronic stress using a well-established, ethologically relevant novel male intruder paradigm for 10 consecutive days during the postpartum period. We hypothesized that the increased energetic burden of mounting a chronic stress response during lactation would decrease lactation performance. Specifically, we predicted that chronic exposure to this social stressor would decrease either milk quality (i.e., composition of proximate components and energy density) or quantity. We also predicted that changes in proximate composition (i.e., lipid, lactose, and protein concentrations) would be associated with changes in gene expression levels of milk synthesis genes. Our results supported our hypothesis that chronic stress impairs lactation performance. Relative to the controls, chronically stressed rats had lower milk yields. We also found that milk quality was decreased; milk from chronically stressed mothers had lower lipid concentration and lower energy density, though protein and lactose concentrations were not different between treatment groups. Although there was a change in proximate composition, chronic stress did not impact mammary gland expression of key milk synthesis genes. Together, this work demonstrates that exposure to a chronic stressor impacts lactation performance, which in turn has the potential to impact offspring development via maternal effects.

Environmental chemical TCPOBOP disrupts milk lipid homeostasis during pregnancy and lactation

Humans are exposed to different kinds of environmental contaminants or drugs throughout their lifetimes. The widespread presence of these compounds has raised concerns about the consequent adverse effects on lactating women. The constitutive androstane receptor (CAR, Nr1i3) is known as a xenobiotic sensor for environmental pollution or drugs. In this study, the model environmental chemical 1, 4-bis [2-(3, 5-dichloropyridyloxy)] benzene, TCPOBOP (TC), which is a highly specific agonist of CAR, was used to investigate the effects of exogenous exposure on lactation function and offspring health in mice. The results revealed that TC exposure decreased the proliferation of mammary epithelial cells during pregnancy. This deficiency further compromised lobular-alveolar structures, resulting in alveolar cell apoptosis, as well as premature stoppage of the lactation cycle and aberrant lactation. Furthermore, TC exposure significantly altered the size and number of milk lipid droplets, suggesting that TC exposure inhibits milk lipid synthesis. Additionally, TC exposure interfered with the milk lipid metabolism network, resulting in the inability of TC-exposed mice to efficiently secrete nutrients and feed their offspring. These findings demonstrated that restricted synthesis and secretion of milk lipids would indirectly block mammary gland form and function, which explained the possible reasons for lactation failure and retarded offspring growth.

Elucidating genes and gene networks linked to individual susceptibility to milk fat depression in dairy goats

Dietary supplementation with marine lipids modulates ruminant milk composition toward a healthier fatty acid profile for consumers, but it also causes milk fat depression (MFD). Because the dairy goat industry is mainly oriented toward cheese manufacturing, MFD can elicit economic losses. There is large individual variation in animal susceptibility with goats more (RESPO+) or less (RESPO-) responsive to diet-induced MFD. Thus, we used RNA-Seq to examine gene expression profiles in mammary cells to elucidate mechanisms underlying MFD in goats and individual variation in the extent of diet-induced MFD. Differentially expression analyses (DEA) and weighted gene co-expression network analysis (WGCNA) of RNA-Seq data were used to study milk somatic cell transcriptome changes in goats consuming a diet supplemented with marine lipids. There were 45 differentially expressed genes (DEGs) between control (no-MFD, before diet-induced MFD) and MFD, and 18 between RESPO+ and RESPO-. Biological processes and pathways such as "RNA transcription" and "Chromatin modifying enzymes" were downregulated in MFD compared with controls. Regarding susceptibility to diet-induced MFD, we identified the "Triglyceride Biosynthesis" pathway upregulated in RESPO- goats. The WGCNA approach identified 9 significant functional modules related to milk fat production and one module to the fat yield decrease in diet-induced MFD. The onset of MFD in dairy goats is influenced by the downregulation of SREBF1, other transcription factors and chromatin-modifying enzymes. A list of DEGs between RESPO+ and RESPO- goats (e.g., DBI and GPD1), and a co-related gene network linked to the decrease in milk fat (ABCD3, FABP3, and PLIN2) was uncovered. Results suggest that alterations in fatty acid transport may play an important role in determining individual variation. These candidate genes should be further investigated.

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.

Bta-miR-33a affects gene expression and lipid levels in Chinese Holstein mammary epithelial cells

MicroRNAs (miRNAs) are non-coding single-stranded RNA molecules of about 19-25 nucleotides in length that regulate different biological processes, including lipid metabolism. In this study, we explored the effect of bta-miR-33a on lipid metabolism in bovine mammary epithelial cells (BMECs) of Chinese Holstein for the first time. For this purpose, the plasmids of bta-miR-33a mimic, bta-miR-33a inhibitor and bta-miR-33a negative control were constructed to overexpress or repress bta-miR-33a in BMECs. The effects of plasmid transfection were analysed by examining the mRNA and protein expression levels of ELOVL6 and the intracellular triglycerides. The results showed that bta-miR-33a directly inhibited the expression of ELOVL6 in BMECs; decreased the mRNA levels of ELOVL5, HACD2, CPT1A and MSMO1; and increased the mRNA level of ALOX15. Sequence bta-miR-33a also increased the contents of triglycerides in the cells, presumably as a consequence of these gene expression changes. In summary, the results of the present study suggest that bta-miR-33a regulates lipid metabolism by targeting ELOVL6, which might be a potential molecular marker of milk fat composition.

Candidate Genes and Their Expressions Involved in the Regulation of Milk and Meat Production and Quality in Goats (Capra hircus)

Simple Summary

During the present decade, highly selected caprine farming has increased in popularity due to the hardiness and adaptability inherent to goats. Recent advances in genetics have enabled the improvement in goat selection efficiency. The present review explores how genetic technologies have been applied to the goat-farming sector in the last century. The main candidate genes related to economically relevant traits are reported. The major source of income in goat farming derives from the sale of milk and meat. Consequently, yield and quality must be specially considered. Meat-related traits were evaluated considering three functional groups (weight gain, carcass quality and fat profile). Milk traits were assessed in three additional functional groups (milk production, protein and fat content).

Abstract

Despite their pivotal position as relevant sources for high-quality proteins in particularly hard environmental contexts, the domestic goat has not benefited from the advances made in genomics compared to other livestock species. Genetic analysis based on the study of candidate genes is considered an appropriate approach to elucidate the physiological mechanisms involved in the regulation of the expression of functional traits. This is especially relevant when such functional traits are linked to economic interest. The knowledge of candidate genes, their location on the goat genetic map and the specific phenotypic outcomes that may arise due to the regulation of their expression act as a catalyzer for the efficiency and accuracy of goat-breeding policies, which in turn translates into a greater competitiveness and sustainable profit for goats worldwide. To this aim, this review presents a chronological comprehensive analysis of caprine genetics and genomics through the evaluation of the available literature regarding the main candidate genes involved in meat and milk production and quality in the domestic goat. Additionally, this review aims to serve as a guide for future research, given that the assessment, determination and characterization of the genes associated with desirable phenotypes may provide information that may, in turn, enhance the implementation of goat-breeding programs in future and ensure their sustainability.

Effect of Protein Nutrition Level on Protein Metabolism during Volleyball Exercise Based on Edge Computing in the Medical System

With the rapid development of the Internet of Things, 5G, and communication technologies, the growth of various types of data has shown an exponential trend. Edge computing technology provides users with almost unlimited computing power through a large number of high-performance servers in the data center. It is one of the important solutions for big data analysis and processing. Volleyball has caused a great wave in China as early as the 1960s, but people pay little attention to the physical quality of volleyball players. At the same time, in the medical field, it is difficult to give a clear value to the athlete's protein requirement. Therefore, this article aims to observe the specific values of protein metabolism in volleyball at different levels of protein nutrition. By designing controlled experiments, then these rats under three nutrient levels of protein were observed and protein metabolism was analyzed after volleyball. The results of the study show that volleyball exercise can reduce the nitrogen balance and gastrocnemius nitrogen content. The nitrogen balance of the 17% group decreased from 388 mg/day before exercise to 336 mg/day, and the gastrocnemius nitrogen content decreased by up to 5.2%; serum urea nitrogen concentration and liver nitrogen content are increased, indicating the enhancement of protein catabolism. Different protein nutrition levels have different effects on protein metabolism during volleyball. The protein intake level of 17% is more conducive to resist the protein decomposition caused by volleyball. It can be seen that, based on edge computing technology, the influence factors of protein nutrition level on protein metabolism during volleyball sports can be well explored, and the research results are also very valuable.

Multi-Omics Integration and Network Analysis Reveal Potential Hub Genes and Genetic Mechanisms Regulating Bovine Mastitis

Mastitis, inflammation of the mammary gland, is the most prevalent disease in dairy cattle that has a potential impact on profitability and animal welfare. Specifically designed multi-omics studies can be used to prioritize candidate genes and identify biomarkers and the molecular mechanisms underlying mastitis in dairy cattle. Hence, the present study aimed to explore the genetic basis of bovine mastitis by integrating microarray and RNA-Seq data containing healthy and mastitic samples in comparative transcriptome analysis with the results of published genome-wide association studies (GWAS) using a literature mining approach. The integration of different information sources resulted in the identification of 33 common and relevant genes associated with bovine mastitis. Among these, seven genes—CXCR1, HCK, IL1RN, MMP9, S100A9, GRO1, and SOCS3—were identified as the hub genes (highly connected genes) for mastitis susceptibility and resistance, and were subjected to protein-protein interaction (PPI) network and gene regulatory network construction. Gene ontology annotation and enrichment analysis revealed 23, 7, and 4 GO terms related to mastitis in the biological process, molecular function, and cellular component categories, respectively. Moreover, the main metabolic-signalling pathways responsible for the regulation of immune or inflammatory responses were significantly enriched in cytokine–cytokine-receptor interaction, the IL-17 signaling pathway, viral protein interaction with cytokines and cytokine receptors, and the chemokine signaling pathway. Consequently, the identification of these genes, pathways, and their respective functions could contribute to a better understanding of the genetics and mechanisms regulating mastitis and can be considered a starting point for future studies on bovine mastitis.

MicroRNA-432 inhibits milk fat synthesis by targeting SCD and LPL in ovine mammary epithelial cells

The microRNA (miR)-432 is differentially expressed in the mammary gland of two breeds of lactating sheep with different milk production traits, and between the non-lactating and peak-lactation periods, but there have been no reports describing the molecular mechanisms involved. In this study, the effect of miR-432 on the proliferation of ovine mammary epithelial cells (OMECs) and the target genes of miR-432 were investigated. The effects of miR-432 on the expression of the target genes and the content of triglycerides in the OMECs were also analyzed. Transfection with a miR-432 mimic was found using CCK8 and Edu assays, to inhibit the viability of OMECs and reduce the number of proliferated OMECs. In contrast, a miR-432 inhibitor had the opposite effect to the miR-432 mimic, and together these results suggest that miR-432 inhibits the proliferation of OMECs. A dual luciferase assay revealed that the genes for stearoyl-CoA desaturase (SCD) and lipoprotein lipase (LPL) are targeted by miR-432. The transfection of miR-432 mimic into OMECs resulted in decreases in the expression of SCD and LPL, and three other milk fat synthesis marker genes; FABP4, LPIN1 and ACACA. The mimic also decreased the content of triglycerides. The miR-432 inhibitor had the opposite effect to the mimic on the expression of these genes and the level of triglycerides. This is the first study to reveal the biological mechanisms by which miR-432 inhibits milk fat synthesis in sheep.

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

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

Evaluation of quality and gene expression of goat embryos produced in vivo and in vitro after cryopreservation

In the present study, we aimed to identify morphological and molecular changes of in vivo and in vitro-produced goat embryos submitted to cryopreservation. In vivo embryos were recovered by transcervical technique from superovulated goats, whereas in vitro produced embryos were produced from ovaries collected at a slaughterhouse. Embryos were frozen by two-steps slow freezing method, which is defined as freezing to -32 °C followed by transfer to liquid nitrogen. Morphological evaluation of embryos was carried out by assessing blastocoel re-expansion rate and the total number of blastomeres. The expression profile of candidate genes related to thermal and oxidative stress, apoptosis, epigenetic, and implantation control was measured using RT-qPCR based SYBR Green system. In silico analyses were performed to identify conserved genes in goat species and protein-protein interaction networks were created. In vivo-produced embryos showed greater blastocoel re-expansion and more blastomere cells (P < 0.05). The expression level of CTP2 and HSP90 genes from in vitro cryopreserved embryos was higher than their in vivo counterparts. Unlikely, no significant difference was observed in the transcription level of SOD gene between groups. The high similarity of CPT2 and HSP90 proteins to their orthologs among mammals indicates that they share conserved functions. In summary, cryopreservation negatively affects the morphology and viability of goat embryos produced in vitro and changes the CPT2 and HSP90 gene expression likely in response to the in vitro production process.

Identification and characterization of circular RNAs in mammary gland tissue from sheep at peak lactation and during the nonlactating period

Circular RNAs are a class of noncoding RNA with a widespread occurrence in eukaryote tissues, and with some having been demonstrated to have clear biological function. In sheep, little is known about the role of circular RNAs in mammary gland tissue, and therefore an RNA sequencing approach was used to compare mammary gland tissue expression of circular RNAs in 9 Small Tail Han sheep at peak lactation, and subsequently when they were not lactating. These 9 sheep had their RNA pooled for analysis into 3 libraries from peak lactation and 3 from the nonlactating period. A total of 3,278 and 1,756 circular RNAs were identified in the peak lactation and nonlactating mammary gland tissues, respectively, and the expression and identity of 9 of them was confirmed using reverse transcriptase-polymerase chain reaction analysis and DNA sequencing. The type, chromosomal location and length of the circular RNAs identified were ascertained. Forty upregulated and one downregulated circular RNAs were characterized in the mammary gland tissue at peak lactation compared with the nonlactating mammary gland tissue. Gene ontology enrichment analysis revealed that the parental genes of these differentially expressed circular RNAs were related to molecular function, binding, protein binding, ATP binding, and ion binding. Five differentially expression circular RNAs were selected for further analysis to predict their target microRNAs, and some microRNAs reportedly associated with the development of the mammary gland were found in the constructed circular RNA-microRNA network. This study reveals the expression profiles and characterization of circular RNAs at 2 key stages of mammary gland activity, thereby providing an improved understanding of the roles of circular RNAs in the mammary gland of sheep.

Comprehensive Transcriptome Profiling of Dairy Goat Mammary Gland Identifies Genes and Networks Crucial for Lactation and Fatty Acid Metabolism

Milk fatty acids secreted by the mammary gland are one of the most important determinants of the nutritional value of goat milk. Unlike cow milk, limited data are available on the transcriptome-wide changes across stages of lactation in dairy goats. In this study, goat mammary gland tissue collected at peak lactation, cessation of milking, and involution were analyzed with digital gene expression (DGE) sequencing to generate longitudinal transcript profiles. A total of 51,299 unigenes were identified and further annotated to 12,763 genes, of which 9,131 were differentially expressed across various stages of lactation. Most abundant genes and differentially expressed genes (DEGs) were functionally classified through clusters of euKaryotic Orthologous Groups (KOG), Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. A total of 16 possible expression patterns were uncovered, and 13 genes were deemed novel candidates for regulation of lactation in the goat: POLG, SPTA1, KLC, GIT2, COPS3, PDP, CD31, USP16/29/37, TLL1, NCAPH, ABI2, DNAJC4, and MAPK8IP3. In addition, PLA2, CPT1, PLD, GGA, SRPRB, and AP4S1 are proposed as novel and promising candidates regulating mammary fatty acid metabolism. “Butirosin and neomycin biosynthesis” and “Glyoxylate and dicarboxylate metabolism” were the most impacted pathways, and revealed novel metabolic alterations in lipid metabolism as lactation progressed. Overall, the present study provides new insights into the synthesis and metabolism of fatty acids and lipid species in the mammary gland along with more detailed information on molecular regulation of lactogenesis. The major findings will benefit efforts to further improve milk quality in dairy goats.

Understanding seasonal weight loss tolerance in dairy goats: a transcriptomics approach

Background

Seasonal weight loss (SWL) is a very important limitation to the production of ruminants in the Mediterranean and Tropical regions. In these areas, long dry seasons lead to poor pastures with low nutritional value. During the dry season, ruminants, particularly those raised in extensive production systems, lose around 30% of their body weight. Seasonal weight loss has important consequences on animal productive performance and health. In this study, RNA sequencing was used to characterize feed restriction effects in dairy goat of 2 breeds with different SWL tolerance: Majorera (tolerant) and Palmera (susceptible). Nine Majorera and ten Palmera goats were randomly distributed in a control and a restricted group: Majorera Control (adequately fed; MC; n = 4), Palmera Control (adequately fed; PC; n = 6), Majorera Restricted (feed restricted; ME; n = 5) and Palmera Restricted (feed restricted; PE; n = 4). On day 22 of the trial, mammary gland biopsies were collected for transcriptomics analysis.

Results

From these samples, 24,260 unique transcripts were identified. From those, 82 transcripts were differentially expressed between MC and ME, 99 between PC and PE, twelve between both control groups and twenty-nine between both restricted groups.

Conclusions

Feed restriction affected several biochemical pathways in both breeds such as: carbohydrate and lipid transport; intracellular trafficking, RNA processing and signal transduction. This research also highlights the importance or involvement of the genes in tolerance (ENPP1, S-LZ, MT2A and GPNB) and susceptibility (GPD1, CTPS1, ELOVL6 and NR4A1) to SWL with respectively higher expression in the Majorera restriced group and the Palmera restricted group in comparison to the control groups. In addition, results from the study may be extrapolated to other dairy ruminant species.

Genome-Wide Association Study towards Genomic Predictive Power for High Production and Quality of Milk in American Alpine Goats

This paper reports an exploratory study based on quantitative genomic analysis in dairy traits of American Alpine goats. The dairy traits are quality-determining components in goat milk, cheese, ice cream, etc. Alpine goat phenotypes for quality components have been routinely recorded for many years and deposited in the Council on Dairy Cattle Breeding (CDCB) repository. The data collected were used to conduct an exploratory genome-wide association study (GWAS) from 72 female Alpine goats originating from locations throughout the U.S. Genotypes were identified with the Illumina Goat 50K single-nucleotide polymorphisms (SNP) BeadChip. The analysis used a polygenic model where the dropping criterion was a call rate ≥ 0.95. The initial dataset was composed of ~60,000 rows of SNPs and 21 columns of phenotypic traits and composed of 53,384 scaffolds containing other informative data points used for genomic predictive power. Phenotypic association with the 50K BeadChip revealed 26,074 reads of candidate genes. These candidate genes segregated as separate novel SNPs and were identified as statistically significant regions for genome and chromosome level trait associations. Candidate genes associated differently for each of the following phenotypic traits: test day milk yield (13,469 candidate genes), test day protein yield (25,690 candidate genes), test day fat yield (25,690 candidate genes), percentage protein (25,690 candidate genes), percentage fat (25,690 candidate genes), and percentage lactose content (25,690 candidate genes). The outcome of this study supports elucidation of novel genes that are important for livestock species in association to key phenotypic traits. Validation towards the development of marker-based selection that provides precision breeding methods will thereby increase the breeding value.

Identification of circRNA-Associated-ceRNA Networks Involved in Milk Fat Metabolism under Heat Stress

Summer temperatures are generally high in Southern China, and cows are likely to suffer a heat stress reaction. Heat stress will have a negative impact on the performance of dairy cows; however, the mechanism by which high temperature affects lactation is not clear. CircRNA is a type of non-coding RNA discovered in recent years, which performs a crucial function in many biological activities. However, the effects of circRNA on lactation function of dairy cows under heat stress is unknown. The present study aimed to explore the expression levels of circRNA in the mammary gland tissue of cows under heat stress. Firstly, we collected blood and milk samples of summer and winter cows and evaluated lactation performance using serum indicators, milk production, and milk composition. Incorporating the calculation of the temperature and humidity index, we conformed the heat stress status of cows in summer. Heat stress increased the concentration of HSP70 and decreased the concentration of SOD and PRL. Heat stress not only reduced milk yield but also affected milk quality, with milk lactose and milk protein decreasing with increased temperature. The analysis of the fatty acid composition in summer milk found significantly reduced concentrations of unsaturated fatty acids, especially long-chain unsaturated fatty acids. Sequencing of the cow's mammary gland transcriptome revealed that compared to the appropriate temperature (ST) group, the heat stress (HS) group had a total of 2204 upregulated and 3501 downregulated transcripts. GO enrichment and KEGG pathway analysis showed that these genes were mainly related to milk fat metabolism. In addition, 19 upregulated and 19 downregulated circRNA candidates were found in response to heat stress. We used Pearson's test to establish the correlation of circRNA-mRNA and identified four pairs of circRNA-miRNA networks between four circRNAs, six miRNAs, and the CD36 gene. In this study, we revealed the possible role of circRNAs in lactation of dairy cows and identified that circRNA-miRNA-mRNA networks might exist in the cow's mammary glands, providing valuable experience for dairy lactation and milk quality.

Comparison of the Transcriptome of the Ovine Mammary Gland in Lactating and Non-lactating Small-Tailed Han Sheep

Small-Tailed Han (STH) sheep are known for their high fecundity, but the survival of lambs is compromised and influences the commercial return from farming these sheep, with this being attributed in part to starvation from insufficient milk production by the ewes. In this study, the transcriptome profiles of the mammary gland of lactating and non-lactating STH ewes were investigated using paired-end RNA sequencing (RNA-Seq). An average of 14,447 genes were found to be expressed at peak-lactation in the STH sheep, while 15,146 genes were expressed in non-lactating ewes. A total of 4,003 differentially expressed genes (DEGs) were identified. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses revealed that the DEGs were associated with a wide range of cellular components, biological processes and metabolic pathways, including binding activities, signaling pathways, cellular structures, and immune responses. The most highly expressed genes at peak-lactation included CSN2, LGB, LALBA, CSN1S1, CSN1S2, and CSN3, and the 10 most highly expressed genes accounted for 61.37% of the total Reads Per Kilobase of transcript, per Million mapped reads (RPKM). The most highly expressed genes in the mammary gland of non-lactating ewes included IgG, THYMB4X, EEF1A1, IgA, and APOE, and the 10 most highly expressed genes accounted for only 12.97% of the total gene RPKM values. This suggests that the sheep mammary gland undergoes a substantial development in milk protein synthesis infrastructure and promotion of protein transportation during lactation.

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.

Emerging Roles of Heat-Induced circRNAs Related to Lactogenesis in Lactating Sows

Heat stress negatively influences milk production and disrupts normal physiological activity of lactating sows, but the precious mechanisms by which hyperthermia adversely affects milk synthesis in sows still remain for further study. Circular RNAs are a novel class of non-coding RNAs with regulatory functions in various physiological and pathological processes. The expression profiles and functions of circRNAs of sows in lactogenesis remain largely unknown. In the present study, long-term heat stress (HS) resulted in a greater concentration of serum HSP70, LDH, and IgG, as well as decreased levels of COR, SOD, and PRL. HS reduced the total solids, fat, and lactose of sow milk, and HS significantly depressed CSNαs1, CSNαs2, and CSNκ biosynthesis. Transcriptome sequencing of lactating porcine mammary glands identified 42 upregulated and 25 downregulated transcripts in HS vs. control. Functional annotation of these differentially-expressed transcripts revealed four heat-induced genes involved in lactation. Moreover, 29 upregulated and 21 downregulated circRNA candidates were found in response to HS. Forty-two positively correlated circRNA-mRNA expression patterns were constructed between the four lactogenic genes and differentially expressed circRNAs. Five circRNA-miRNA-mRNA post-transcriptional networks were identified involving genes in the HS response of lactating sows. In this study we establish a valuable resource for circRNA biology in sow lactation. Analysis of a circRNA-miRNA-mRNA network further uncovered a novel layer of post-transcriptional regulation that could be used to improve sow milk production.

Effects of Dietary Vegetable Oils on Mammary Lipid-Related Genes in Holstein Dairy Cows

Simple Summary

This study analyzed effects of vegetable oils fed to dairy cows on abundance of genes related to lipid metabolism in milk somatic cells (MSC). During 63 days (9 weeks), 15 cows were allocated to 3 treatments: a control diet with no added lipid and the same diet supplemented with olive oil (OO, 30 g/kg DM) or hydrogenated vegetable oil (HVO, 30 g/kg DM). Dietary oil supplementation (3% DM) had a modest nutrigenomic effect on biological functions such as acetate and FA activation and intra-cellular transport, lipid droplet formation, and transcription regulation in MSC. Results suggest that long-term dietary monounsaturated and saturated lipids could alter mRNA abundance in MSC from mid-lactating cows.

Abstract

This study analyzed effects of vegetable oils fed to dairy cows on abundance of genes related to lipid metabolism in milk somatic cells (MSC). During 63 days, 15 cows were allocated to 3 treatments: a control diet with no added lipid the same diet supplemented with olive oil (OO, 30 g/kg DM) or hydrogenated vegetable oil (HVO, 30 g/kg DM). On days 21, 42 and 63, MSC were obtained from all cows. Relative abundance of genes involved in lipid metabolism in MSC from cows fed control on days 42 and 63 was compared with relative abundance at day 21 to evaluate fold-changes. Those genes without changes over the time were selected to analyze effects of OO and HVO. Compared with control, on day 42, PLIN2 and THRSP were upregulated by OO. Compared with control, on day 21, HVO up regulated ACACA, down regulated FABP3, and on day 63 THRSP and FABP4 were down regulated. Dietary oil supplementation (3% DM) had a modest nutrigenomic effect on different biological functions such as acetate and FA activation and intra-cellular transport, lipid droplet formation, and transcription regulation in MSC.

Identification of the complete coding cDNAs and expression analysis of B4GALT1, LALBA, ST3GAL5, ST6GAL1 in the colostrum and milk of the Garganica and Maltese goat breeds to reveal possible implications for oligosaccharide biosynthesis

Background

Milk sialylated oligosaccharides (SOS) play crucial roles in many biological processes. The most abundant free SOS in goat’s milk are 3’sialyllactose (3′-SL), 6’sialyllactose (6′-SL) and disialyllactose (DSL). The production of these molecules is determined genetically by the expression of glycosyltransferases and by the availability of nucleotide sugar substrates, but the precise mechanisms regulating the differential patterns of milk oligosaccharides are not known. We aimed to identify the complete cDNAs of candidate genes implicated in SOS biosynthesis (B4GALT1, LALBA, ST3GAL5, ST6GAL1) and to analyse their expression during lactation in the Garganica and Maltese goat breeds. Moreover, we analysed the colostrum and milk contents of 3′-SL, 6′-SL and disialyllactose (DSL) and the possible correlations between expressed genes and SOS.

Results

We identified the complete coding cDNAs of B4GALT1 (HQ700335.1), ST3GAL5 (KF055858.2), and ST6GAL1 (HQ709167.1), the single nucleotide polymorphism (SNPs) of these genes and 2 splicing variants of the ST6GAL1 cDNA. RT-qPCR analysis showed that LALBA and ST6GAL1 were the genes with the highest and lowest expression in both breeds, respectively. The interaction effects of the breeds and sampling times were associated with higher levels of B4GALT1 and ST3GAL5 gene expression in Garganica than in Maltese goats at kidding. B4GALT1, LALBA, and ST3GAL5 gene expression changed from kidding to 60 and 120 days in Maltese goats, while in Garganica goats, a difference was observed only for the LALBA gene. Breed and lactation effects were also found for SOS contents. Positive correlations of B4GALT1, LALBA, ST3GAL5, and ST6GAL1 with 3′-SL/6′SL and DSL were found.

Conclusions

The genetic effect on the oligosaccharide content of milk was previously highlighted in bovines, and this study is the first to investigate this effect in two goat breeds (Garganica and Maltese) during lactation. The genetic variability of candidate genes involved in SOS biosynthesis highlights their potential role in affecting gene expression and ultimately biological function. The investigation of gene regulatory regions as well as the examination of other sialyltransferase genes will be needed to identify the genetic pattern leading to a higher SOS content in the autochtonous Garganica breed and to protect it using a focused breeding strategy.

Proteomic changes in the milk of water buffaloes (Bubalus bubalis) with subclinical mastitis due to intramammary infection by Staphylococcus aureus and by non-aureus staphylococci

Subclinical mastitis by Staphylococcus aureus (SAU) and by non-aureus staphylococci (NAS) is a major issue in the water buffalo. To understand its impact on milk, 6 quarter samples with >3,000,000 cells/mL (3 SAU-positive and 3 NAS-positive) and 6 culture-negative quarter samples with <50,000 cells/mL were investigated by shotgun proteomics and label-free quantitation. A total of 1530 proteins were identified, of which 152 were significantly changed. SAU was more impacting, with 162 vs 127 differential proteins and higher abundance changes (P < 0.0005). The 119 increased proteins had mostly structural (n = 43, 28.29%) or innate immune defence functions (n = 39, 25.66%) and included vimentin, cathelicidins, histones, S100 and neutrophil granule proteins, haptoglobin, and lysozyme. The 33 decreased proteins were mainly involved in lipid metabolism (n = 13, 59.10%) and included butyrophilin, xanthine dehydrogenase/oxidase, and lipid biosynthetic enzymes. The same biological processes were significantly affected also upon STRING analysis. Cathelicidins were the most increased family, as confirmed by western immunoblotting, with a stronger reactivity in SAU mastitis. S100A8 and haptoglobin were also validated by western immunoblotting. In conclusion, we generated a detailed buffalo milk protein dataset and defined the changes occurring in SAU and NAS mastitis, with potential for improving detection (ProteomeXchange identifier PXD012355).

Fatty Acid Elongase 7 (ELOVL7) Plays a Role in the Synthesis of Long-Chain Unsaturated Fatty Acids in Goat Mammary Epithelial Cells

Simple Summary

Very-long-chain elongases are required for the synthesis of essential fatty acids in non-ruminants. Whether the fatty acid elongase 7 (ELOVL7) plays a role in ruminants is unclear. We demonstrated, in goat mammary epithelial cells, that ELOVL7 activation resulted in greater concentrations of vaccenic (C18:1n7) and linoleic (C18:2) acid, and lower concentrations of palmitoleic (C16:1n7) and oleic (C18:1n9) acid. Knockdown of ELOVL7 increased the concentration of C18:1n9. The data support a novel role of ELOVL7 in altering long-chain unsaturated fatty acids in goat mammary epithelial cells.

Abstract

In humans, fatty acid elongase 7 (ELOVL7) plays a role in synthesis of long-chain saturated fatty acids. Whether ELOVL7 protein plays a role in ruminants is unclear. The transcript abundance of ELOVL7 in goat mammary tissue was assessed at three stages of lactation. Results showed that ELOVL7 had the highest expression in the dry period compared with peak and late lactation period. Results revealed that ELOVL7 overexpression was correlated with lower expression in diacylglycerol O-acyltransferase 2 (DGAT2) and stearoyl-CoA desaturase 1 (SCD1), and had no significant effect on triacylglycerol concentration. Overexpression of ELOVL7 significantly decreased the concentration of palmitoleic (C16:1n7) and oleic (C18:1n9) acid, and increased the concentration of vaccenic (C18:1n7) and linoleic (C18:2) acid. Overexpression of ELOVL7 significantly upregulated the elongation index of C16:1 in goat epithelial mammary cells (GMEC), but had a minor effect on that of palmitate (C16:0). Knockdown of ELOVL7 decreased mRNA expression of fatty acid binding protein 3 (FABP3) and fatty acid desaturase 2 (FADS2) and had a minor effect on triacylglycerol concentration; however, it increased concentration of C18:1n9 in GMEC. The elongation indices of C16:0 and C16:1 did not differ due to knockdown of ELOVL7. The minor change for the C16:0 and stearate (C18:0) was observed after activation of ELOVL7, suggesting the two fatty acids are not the preferential substrates of ELOVL7 in cultured GMEC. However, changes in C18:1n9 and C18:2 after overexpression or knockdown of ELOVL7 indicated a biological functional role of ELOVL7. Collectively, our data highlighted a role of ELOVL7 in long-chain unsaturated fatty acid elongation in goat mammary epithelial cells.

Integrative Analysis of Transcriptome and GWAS Data to Identify the Hub Genes Associated With Milk Yield Trait in Buffalo

The mammary gland is the production organ in mammals that is of great importance for milk production and quality. However, characterization of the buffalo mammary gland transcriptome and identification of the valuable candidate genes that affect milk production is limited. Here, we performed the differential expressed genes (DEGs) analysis of mammary gland tissue on day 7, 50, 140, and 280 after calving and conducted gene-based genome-wide association studies (GWAS) of milk yield in 935 Mediterranean buffaloes. We then employed weighted gene co-expression network analysis (WGCNA) to identify specific modules and hub genes related to milk yield based on gene expression profiles and GWAS data. The results of the DEGs analysis showed that a total of 1,420 DEGs were detected across different lactation points. In the gene-based analysis, 976 genes were found to have genome-wide association (P ≤ 0.05) that could be defined as the nominally significant GWAS geneset (NSGG), 9 of which were suggestively associated with milk yield (P < 10⁻⁴). Using the WGCNA analysis, 544 and 225 genes associated with milk yield in the turquoise module were identified from DEGs and NSGG datasets, respectively. Several genes (including BNIPL, TUBA1C, C2CD4B, DCP1B, MAP3K5, PDCD11, SRGAP1, GDPD5, BARX2, SCARA3, CTU2, and RPL27A) were identified and considered as the hub genes because they were involved in multiple pathways related to milk production. Our findings provide an insight into the dynamic characterization of the buffalo mammary gland transcriptome, and these potential candidate genes may be valuable for future functional characterization of the buffalo mammary gland.

Recent progress of porcine milk components and mammary gland function

As the only nutritional source for newborn piglets, porcine colostrum and milk contain critical nutritional and immunological components including carbohydrates, lipids, and proteins (immunoglobulins). However, porcine milk composition is more complex than these three components. Recently, scientists identified additional and novel components of sow colostrum and milk, including exosomes, oligosaccharides, and bacteria, which possibly act as biological signals and modulate the intestinal environment and immune status in piglets and later in life. Evaluation of these nutritional and non-nutritional components in porcine milk will help better understand the nutritional and biological function of porcine colostrum and milk. Furthermore, some important functions of the porcine mammary gland have been reported in recent published literature. These preliminary studies hypothesized how glucose, amino acids, and fatty acids are transported from maternal blood to the porcine mammary gland for milk synthesis. Therefore, we summarized recent reports on sow milk composition and porcine mammary gland function in this review, with particular emphasis on macronutrient transfer and synthesis mechanisms, which might offer a possible approach for regulation of milk synthesis in the future.

CRISPR/Cas9-mediated Stearoyl-CoA Desaturase 1 (SCD1) Deficiency Affects Fatty Acid Metabolism in Goat Mammary Epithelial Cells

Stearoyl-CoA desaturase 1 (SCD1) is a fatty acid desaturase catalyzing cis-double-bond formation in the Δ9 position to produce monounsaturated fatty acids essential for the synthesis of milk fat. Previous studies using RNAi methods have provided support for a role of SCD1 in goat mammary epithelial cells (GMEC); however, RNAi presents several limitations that might preclude a truthful understanding of the biological function of SCD1. To explore the function of SCD1 on fatty acid metabolism in GMEC, we used CRISPR-Cas9-mediated SCD1 knockout through non-homologous end-joining (NHEJ) and homology-directed repair (HDR) pathways in GMEC. We successfully introduced nucleotide deletions and mutations in the SCD1 gene locus through the NHEJ pathway and disrupted its second exon via insertion of an EGFP-PuroR segment using the HDR pathway. In clones derived from the latter, gene- and protein-expression data indicated that we obtained a monoallelic SCD1 knockout. A T7EN1-mediated assay revealed no off-targets in the surveyed sites. The contents of triacylglycerol and cholesterol and the desaturase index were significantly decreased as a consequence of SCD1 knockout. The deletion of SCD1 decreased the expression of other genes involved in de novo fatty acid synthesis, including SREBF1 and FASN, as well the fatty acid transporters FABP3 and FABP4. The downregulation of these genes partly explains the decrease of intracellular triacylglycerols. Our results indicate a successful SCD1 knockout in goat mammary cells using CRISPR-Cas9. The demonstration of the successful use of CRISPR-Cas9 in GMEC is an important step to producing transgenic goats to study mammary biology in vivo.

Examination of the xanthosine response on gene expression of mammary epithelial cells using RNA-seq technology

BACKGROUND

Xanthosine treatment has been previously reported to increase mammary stem cell population and milk production in cattle and goats. However, the underlying molecular mechanisms associated with the increase in stem cell population and milk production remain unclear.

METHODS

Primiparous Beetal goats were assigned to the study. Five days post-partum, one mammary gland of each goat was infused with xanthosine (TRT) twice daily (2×) for 3 days consecutively, and the other gland served as a control (CON). Milk samples from the TRT and CON glands were collected on the 10th day after the last xanthosine infusion and the total RNA was isolated from milk fat globules (MEGs). Total RNA in MFGs was mainly derived from the milk epithelial cells (MECs) as evidenced by expression of milk synthesis genes. Significant differentially expressed genes (DEGs) were subjected to Gene Ontology (GO) terms using PANTHER and gene networks were generated using STRING db.

RESULTS

Preliminary analysis indicated that each individual goat responded to xanthosine treatment differently, with this trend being correlated with specific DEGs within the same animal's mammary gland. Several pathways are impacted by these DEGs, including cell communication, cell proliferation and anti-microbials.

CONCLUSIONS

This study provides valuable insights into transcriptomic changes in milk producing epithelial cells in response to xanthosine treatment. Further characterization of DEGs identified in this study is likely to delineate the molecular mechanisms of increased milk production and stem or progenitor cell population by the xanthosine treatment.

GLUT1 and lactose synthetase are critical genes for lactose synthesis in lactating sows

Background

Lactose synthesis rate is an important factor in milk production and quality in mammals. Understanding the lactose synthesis mechanism is crucial for the improvement of milk quantity and quality. However, research on the temporal gene changes regarding lactose synthesis during the whole lactation is still limited. The objective of this study was to determine gene expression profiles related to lactose synthesis in sows during lactation, and further identify the critical steps or key factors in the lactose synthesis pathway.

Methods

To determine the temporal change of factors related to lactose synthesis in sows, milk from eight multiparous Yorkshire sows (parity 3 to 6) was collected at 0 h, 2 h, 6 h, 12 h, 24 h, day 2, 3, 4, 7, 14, and 21 after birth of the first piglet. Lactose content, prolactin and progesterone concentration, and gene or protein expression related to lactose synthesis were measured.

Results

The lactose yield increased gradually from D2 to D21 and reached a maximum at D14 (3-fold from D2) during lactation (P < 0.05). A similar trend was observed in IGF-1 and insulin concentrations in milk, both of which were greatest at D3 with a subsequent decrease during middle to late lactation. Conversely, milk prolactin and progesterone concentrations moderately decreased with the progression of lactation. The mRNA or protein expressions related to glucose transportation (GLUT1), glucose-galactose interconversion (HK1 and UGP2), UDP-galactose transportation (SLC35A2), and lactose synthetase (LALBA and B4GALT1) in the lactose synthesis pathway were significantly upregulated during early to middle lactation and plateaued by late lactation (P < 0.05).

Conclusions

These novel findings suggest that the increased lactose synthesis in lactation was related to the coordinated upregulation of genes or enzymes in the lactose synthesis pathway, and glucose transportation (GLUT1) and lactose synthetase (LALBA and B4GALT1) might be the critical steps in the lactose synthesis pathway of sows during lactation.

Transcriptomic profiles of the bovine mammary gland during lactation and the dry period

The initiation and maintenance of lactation are complex phenomena governed by biochemical and endocrine processes in the mammary gland (MG). Although DNA-based approaches have been used to study the onset of lactation, more comprehensive RNA-based techniques may be critical in furthering our understanding of gene alterations that occur to support lactation in the bovine MG. To further determine how gene profiles vary during lactation compared with the dry period, RNA-seq transcriptomic analysis was used to identify differentially expressed genes (DEG) in bovine MG tissues from animals that were lactating and not lactating. A total of 881 DEG (605 upregulated and 276 downregulated) were identified in MG of 3 lactating Chinese Holstein dairy cows versus the 3 dry cows. The subcellular analysis showed that the upregulated genes were most abundantly located in "integral to membrane" and "mitochondrion," and the top number of downregulated genes existed in "nucleus" and "cytoplasm." The functional analysis indicated that the DEG were primarily associated with the support of lactation processes. The genes in higher abundance were most related to "metabolic process," "oxidation-reduction process," "transport" and "signal transduction," protein synthesis-related processes (transcription, translation, protein modifications), and some MG growth-associated processes (cell proliferation/cycle/apoptosis). The downregulated genes were mainly involved in immune-related processes (inflammatory/immune/defense responses). The KEGG analysis suggested that protein synthesis-related pathways (such as protein digestion and absorption; protein processing in endoplasmic reticulum; and glycine, serine, and threonine metabolism) were highly and significantly enriched in the bovine MG of lactating cows compared to dry cows. The results suggested that the dry cows had decreased capacity for protein synthesis, energy generation, and cell growth but enhanced immune response. Collectively, this reduced capacity in dry cows supports the physiological demands of the next lactation and the coordinated metabolic changes that occur to support these demands. A total of 51 identified DEG were validated by RT-PCR, and consistent results were found between RT-PCR and the transcriptomic analysis. This work provides a profile of gene-associated changes that occur during lactation and can be used to facilitate further investigation of the mechanisms underlying lactation in dairy cows.

Fatty acid elongase 5 (ELOVL5) alters the synthesis of long-chain unsaturated fatty acids in goat mammary epithelial cells

Increased production of long-chain unsaturated fatty acids (LCUFA) can have a positive effect on the nutritional value of ruminant milk for human consumption. In nonruminant species, fatty acid elongase 5 (ELOVL5) is a key enzyme for endogenous synthesis of long-chain unsaturated fatty acids. However, whether ELOVL5 protein plays a role (if any) in ruminant mammary tissue remains unclear. In the present study, we assessed the mRNA abundance of ELOVL5 at 3 stages of lactation in goat mammary tissue. Results revealed that ELOVL5 had the lowest expression at peak lactation compared with the nonlactating and late-lactating periods. The ELOVL5 was overexpressed or knocked down to assess its role in goat mammary epithelial cells. Results revealed that ELOVL5 overexpression increased the expression of perilipin2 (PLIN2) and decreased diacylglycerolacyltransferase 2 (DGAT2) and fatty acid desaturase 2 (FADS2) mRNA, but had no effect on the expression of DGAT1, FADS1, and stearoyl-CoA desaturase 1 (SCD1). Overexpression of ELOVL5 decreased the concentration of C16:1n-7, whereas no significant change in C18:1n-7 and C18:1n-9 was observed. Knockdown of ELOVL5 decreased the expression of PLIN2 but had no effect on DGAT1, DGAT2, FADS1, FADS2, and SCD1 mRNA expression. Knockdown of ELOVL5 increased the concentration of C16:1n-7 and decreased that of C18:1n-7. The alterations of expression of genes related to lipid metabolism after overexpression or knockdown of ELOVL5 suggested a negative feedback regulation by the products of ELOVL5 activation. However, the content of triacylglycerol was not altered by knockdown or overexpression of ELOVL5 in goat mammary epithelial cells, which might have been due to the insufficient availability of substrate in vitro. Collectively, these are the first in vitro results highlighting an important role of ELOVL5 in the elongation of 16-carbon to 18-carbon unsaturated fatty acids in ruminant mammary cells.

Akt Serine/Threonine Kinase 1 Regulates de Novo Fatty Acid Synthesis through the Mammalian Target of Rapamycin/Sterol Regulatory Element Binding Protein 1 Axis in Dairy Goat Mammary Epithelial Cells

Akt serine/threonine kinase acts as a central mediator in the phosphatidylinositol 3-kinase (PI3K)/Akt signaling pathway, regulating a series of biological processes. In lipid metabolism, Akt activation regulates a series of gene expressions, including genes related to intracellular fatty acid synthesis. However, the regulatory mechanisms of Akt in dairy goat mammary lipid metabolism have not been elaborated. In this study, the coding sequences of goat Akt1 gene were cloned and analyzed. Gene expression of Akt1 in different lactation stages was also investigated. For in vitro studies, a eukaryotic expression vector of Akt1 was constructed and transfected to goat mammary epithelial cells (GMECs), and specific inhibitors of Akt/mammalian target of rapamycin (mTOR) signaling were applied to GMECs. Results showed that Akt1 protein was highly conserved, and its mRNA was highly expressed in midlactation. In vitro studies indicated that Akt1 phosphorylation activated mTOR and subsequently enhanced sterol regulatory element binding protein 1 (SREBP1), thus increasing intracellular triacylglycerol content. Inhibition of Akt/mTOR signaling down-regulated the gene expression of lipogenic genes. Overall, Akt1 plays an important role in regulating de novo fatty acid synthesis in goat mammary epithelial cells, and this process probably is through the mTOR/SREBP1 axis.

Peroxisome proliferator-activated receptor delta regulates lipid droplet formation and transport in goat mammary epithelial cells

Even though recent evidence in goat mammary epithelial cells (GMEC) suggest a role of peroxisome proliferator-activated receptor delta (PPARD) in regulating lipid homeostasis, its role is not fully understood. Our hypothesis was that PPARD regulates lipid transport processes in GMEC and, thus, plays a crucial role in regulating fat formation. The PPARD was overexpressed using an adenovirus system (Ad-PPARD) with recombinant green fluorescent protein (Ad-GFP) as the control. Results revealed that overexpression of PPARD markedly upregulated the mRNA abundance of PPARD. Compared with the control (Ad-GFP+dimethyl sulfoxide), overexpression of PPARD alone had no effect on mRNA expression of CD36, SCD1, FABP4, ACSL1, and ADRP. The cultures overexpressing PPARD with the PPARD ligand GW0742 (GW) upregulated the expression of CD36, FABP3, FABP4, ACSL1, and ADRP. Overexpression of PPARD in GMEC plus GW increased the concentration of 16:1 and 18:1-trans and was associated with upregulation of SCD1. Compared with the control (Ad-GFP+dimethyl sulfoxide), the decrease of triacylglycerol concentration coupled with upregulation of genes related to lipid droplet secretion (e.g., ADRP and ACSL1) induced by PPARD overexpression suggests a role in lipid droplet (LD) secretion. Luciferase assay revealed that GW increased the ADRP promoter activity in a dose-dependent manner. Knockdown of PPARD impaired the increase of ADRP promoter activity induced by GW, whereas GW enhanced the activity of ADRP promoter in GMEC overexpressing PPARD. Data with the ADRP 5'-flanking truncated luciferase reporter suggest a core region (-1,444 to -990 bp) response element for the induction of GW. This core region contains a known PPARG response element (PPRE) at -1,003 to -990 bp. When the PPRE was mutated, the overexpression of PPARD had no effect on ADRP promoter activity. Collectively, these results reveal a novel role for PPARD in lipid homeostasis via promoting fatty acid transport and LD formation through a mechanism of direct binding to the promoter of key genes. Hence, PPARD activity may contribute to fatty acid transport and LD formation during lactation.

Identification of differentially expressed genes through RNA sequencing in goats (Capra hircus) at different postnatal stages

Intramuscular fat (IMF) content and fatty acid composition of longissimus dorsi muscle (LM) change with growth, which partially determines the flavor and nutritional value of goat (Capra hircus) meat. However, unlike cattle, little information is available on the transcriptome-wide changes during different postnatal stages in small ruminants, especially goats. In this study, the sequencing reads of goat LM tissues collected from kid, youth, and adult period were mapped to the goat genome. Results showed that out of total 24 689 Unigenes, 20 435 Unigenes were annotated. Based on expected number of fragments per kilobase of transcript sequence per million base pairs sequenced (FPKM), 111 annotated differentially expressed genes (DEGs) were identified among different postnatal stages, which were subsequently assigned to 16 possible expression patterns by series-cluster analysis. Functional classification by Gene Ontology (GO) analysis was used for selecting the genes showing highest expression related to lipid metabolism. Finally, we identified the node genes for lipid metabolism regulation using co-expression analysis. In conclusion, these data may uncover candidate genes having functional roles in regulation of goat muscle development and lipid metabolism during the various growth stages in goats.

MEN1/Menin regulates milk protein synthesis through mTOR signaling in mammary epithelial cells

The MEN1 gene, which encodes the protein Menin, was investigated for its regulatory role in milk protein synthesis in mammary glands. Menin responds to nutrient and hormone levels via the PI3K/Akt/mTOR pathway. Bovine mammary epithelial cells and tissues were used as experimental models in this study. The results revealed that the milk protein synthesis capacity of mammary epithelial cells could be regulated by MEN1/Menin. The overexpression of Menin caused significant suppression of factors involved in the mTOR pathway, as well as milk protein κ-casein (CSNK). In contrast, a significant increase in these factors and CSNK was observed upon MEN1/Menin knockdown. The repression of MEN1/Menin on the mTOR pathway was also observed in mammary gland tissues. Additionally, MEN1/Menin was found to elicit a negative response on prolactin (PRL) and/or insulin (INS), which caused a similar downstream impact on mTOR pathway factors and milk proteins. Collectively, our data indicate that MEN1/Menin could play a regulatory role in milk protein synthesis through mTOR signaling in the mammary gland by mediating the effects of hormones and nutrient status. The discovery of Menin's role in mammary glands suggests Menin could be potential new target for the improvement of milk performance and adjustment of lactation period of dairy cows.

Variants in the 3' untranslated region of the ovine acetyl-coenzyme A acyltransferase 2 gene are associated with dairy traits and exhibit differential allelic expression

The acetyl-CoA acyltransferase 2 (ACAA2) gene encodes an enzyme of the thiolase family that is involved in mitochondrial fatty acid elongation and degradation by catalyzing the last step of the respective β-oxidation pathway. The increased energy needs for gluconeogenesis and triglyceride synthesis during lactation are met primarily by increased fatty acid oxidation. Therefore, the ACAA2 enzyme plays an important role in the supply of energy and carbon substrates for lactation and may thus affect milk production traits. This study investigated the association of the ACAA2 gene with important sheep traits and the putative functional involvement of this gene in dairy traits. A single nucleotide substitution, a T to C transition located in the 3' untranslated region of the ACAA2 gene, was used in mixed model association analysis with milk yield, milk protein yield and percentage, milk fat yield and percentage, and litter size at birth. The single nucleotide polymorphism was significantly associated with total lactation production and milk protein percentage, with respective additive effects of 6.81 ± 2.95 kg and -0.05 ± 0.02%. Additionally, a significant dominance effect of 0.46 ± 0.21 kg was detected for milk fat yield. Homozygous TT and heterozygous CT animals exhibited higher milk yield compared with homozygous CC animals, whereas the latter exhibited increased milk protein percentage. Expression analysis from age-, lactation-, and parity-matched female sheep showed that mRNA expression of the ACAA2 gene from TT animals was 2.8- and 11.8-fold higher in liver and mammary gland, respectively. In addition, by developing an allelic expression imbalance assay, it was estimated that the T allele was expressed at an average of 18% more compared with the C allele in the udder of randomly selected ewes. We demonstrated for the first time that the variants in the 3' untranslated region of the ovine ACAA2 gene are differentially expressed in homozygous ewes of each allele and exhibit allelic expression imbalance within heterozygotes in a tissue-specific manner, supporting the existence of cis-regulatory DNA variation in the ovine ACAA2 gene. This is the first study reporting differential allelic imbalance expression of a candidate gene associated with milk production traits in dairy sheep.

Fatty acid elongase 6 plays a role in the synthesis of long-chain fatty acids in goat mammary epithelial cells

In nonruminants, it is well established that elongation of very long-chain fatty acid-like fatty acid elongase 6 (ELOVL6) catalyzes the synthesis of C18:0 from C16:0 in lipogenic tissues like adipose and liver. However, the role of ELOVL6 in regulating lipid metabolism in ruminant mammary gland remains unknown. In the present study, ELOVL6 was overexpressed or knocked down via adenoviral transfection to assess its role in goat mammary epithelial cells. Results revealed that ELOVL6 overexpression had a weak effect on the expression of genes related to triacylglycerol (TAG) synthesis and desaturation. Overexpression of ELOVL6 increased the content of C18:0 at the expense of C16:0, and increased the elongation index of C16:0. Overexpression of ELOVL6 had no significant effect on the elongation index of C16:1n-7 and the desaturation indices of C16:0 and C18:0. Knockdown of ELOVL6 had a negative effect on mRNA expression of the esterification genes GPAM and diacylglycerolacyltransferase 2 (DGAT2) and TAG concentration; however, it increased the concentration of C16:0 and decreased C18:1n-7 and C18:1n-9 in goat mammary epithelial cells. Accordingly, downregulation of ELOVL6 significantly decreased the elongation indices of C16:0 and C16:1n-7. The lack of change in the desaturation indices of C16:0 and C18:0 upon knockdown of ELOVL6 was consistent with the minor change in SCD1 expression. In conclusion, these are the first results highlighting an important role of ELOVL6 in long-chain fatty elongation and TAG synthesis in ruminant mammary cells.

miR-148a and miR-17-5p synergistically regulate milk TAG synthesis via PPARGC1A and PPARA in goat mammary epithelial cells

MicroRNA (miRNA) are a class of '18-25' nt RNA molecules which regulate gene expression and play an important role in several biologic processes including fatty acid metabolism. Here we used S-Poly (T) and high-throughput sequencing to evaluate the expression of miRNA and mRNA during early-lactation and in the non-lactating ("dry") period in goat mammary gland tissue. Results indicated that miR-148a, miR-17-5p, PPARGC1A and PPARA are highly expressed in the goat mammary gland in early-lactation and non-lactating periods. Utilizing a Luciferase reporter assay and Western Blot, PPARA, an important regulator of fatty acid oxidation, and PGC1a (PPARGC1A), a major regulator of fat metabolism, were demonstrated to be targets of miR-148a and miR-17-5p in goat mammary epithelial cells (GMECs). It was also revealed that miR-148a expression can regulate PPARA, and miR-17-5p represses PPARGC1A in GMECs. Furthermore, the overexpression of miR-148a and miR-17-5p promoted triacylglycerol (TAG) synthesis while the knockdown of miR-148a and miR-17-5p impaired TAG synthesis in GMEC. These findings underscore the importance of miR-148a and miR-17-5p as key components in the regulation of TAG synthesis. In addition, miR-148a cooperates with miR-17-5p to regulate fatty acid metabolism by repressing PPARGC1A and PPARA in GMECs. Further studies on the functional role of miRNAs in lipid metabolism of ruminant mammary cells seem warranted.

RNA-Sequencing for profiling goat milk transcriptome in colostrum and mature milk

BACKGROUND

In this work we aimed at sequencing and assembling the goat milk transcriptome corresponding at colostrum and 120 days of lactation. To reconstruct transcripts we used both the genome as reference, and a de novo assembly approach. Additionally, we aimed at identifying the differentially expressed genes (DEGs) between the two lactation stages and at analyzing the expression of genes involved in oligosaccharides metabolism.

RESULTS

A total of 44,635 different transcripts, organized in 33,757 tentative genes, were obtained using the goat genome as reference. A significant sequence similarity match was found for 40,353 transcripts (90%) against the NCBI NT and for 35,701 (80%) against the NR databases. 68% and 69% of the de novo assembled transcripts, in colostrum and 120 days of lactation samples respectively, have a significant match with the merged transcriptome obtained using Cufflinks/Cuffmerge. CSN2, PAEP, CSN1S2, CSN3, LALBA, TPT1, FTH1, M-SAA3, SPP1, GLYCAM1, EEF1A1, CTSD, FASN, RPS29, CSN1S1, KRT19 and CHEK1 were found between the top fifteen highly expressed genes. 418 loci were differentially expressed between lactation stages, among which 207 and 122 were significantly up- and down-regulated in colostrum, respectively. Functional annotation and pathway enrichment analysis showed that in goat colostrum somatic cells predominate biological processes involved in glycolysis, carbohydrate metabolism, defense response, cytokine activity, regulation of cell proliferation and cell death, vasculature development, while in mature milk, biological process associated with positive regulation of lymphocyte activation and anatomical structure morphogenesis are enriched. The analysis of 144 different oligosaccharide metabolism-related genes showed that most of these (64%) were more expressed in colostrum than in mature milk, with eight expressed at very high levels (SLCA3, GMSD, NME2, SLC2A1, B4GALT1, B3GNT2, NANS, HEXB).

CONCLUSIONS

To our knowledge, this is the first study comparing goat transcriptome of two lactation stages: colostrum and 120 days. Our findings suggest putative differences of expression between stages and can be envisioned as a base for further research in the topic. Moreover because a higher expression of genes involved in immune defense response, carbohydrate metabolism and related to oligosaccharide metabolism was identified in colostrum we here corroborate the potential of goat milk as a natural source of lactose-derived oligosaccharides and for the development of functional foods.