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

Research Progress on Genomic Regions and Candidate Genes Related to Milk Composition Traits of Dairy Goats Based on Functional Genomics: A Narrative Review

BACKGROUND

Goat milk has gained global attention for its unique nutritional properties and potential health benefits. Advancements in functional genomic technologies have significantly progressed genetic research on milk composition traits in dairy goats.

RESULTS

This review summarizes various research methodologies applied in this field. Genome-wide association studies (GWAS) have identified genomic regions associated with major milk components, with the diacylglycerol acyltransferase 1 (DGAT1) gene and casein gene cluster consistently linked to milk composition traits. Transcriptomics has revealed gene expression patterns in mammary tissue across lactation stages, while the role of non-coding RNAs (such as miRNAs and circRNAs) in regulating milk composition has been confirmed. Proteomic and metabolomic studies have not only helped us gain a more comprehensive understanding of goat milk composition characteristics but have also provided crucial support for the functional validation of genes related to milk components. The integration of multi-omics data has emerged as an effective strategy for elucidating complex regulatory networks from a systems biology perspective.

CONCLUSIONS

Despite progress, challenges remain, including refining reference genomes, collecting large-scale phenotypic data, and conducting functional validations. Future research should focus on improving reference genomes, expanding study populations, investigating functional milk components, exploring epigenetic regulation and non-coding RNAs, and studying microbiome-host genome interactions. These efforts will inform more precise genomic and marker-assisted selection strategies, advancing genetic improvements in milk composition traits in dairy goats.

Proteomic and Transcriptomic Profiling Revealed Vital Molecular Events in the Transition from Goat Colostrum to Mature Milk

As an important nutrient source in large areas of the world, goat milk is favored by more and more consumers; however, the composition, nutritional value, and regulation mechanism of goat milk are not fully characterized. Mammary gland development is as important as detailed annotation of protein composition to address the physiological and nutritional values of goat milk. In the present study, 4353 colostrum and mature goat milk proteins were identified. The abundance of 118 proteins was significantly different between colostrum and mature milk proteins. Our results indicate that the milk protein changes were associated with a network of mammary gene expression changes; importantly, the prime factors include enhanced mammary growth/development, decreased protein translation, attenuated protein folding, and lower lip/carbohydrate metabolism. The present study provides insights into the changes in mammary metabolisms during the transition from colostrum to mature milk, which can help deeply explore the difference and regulation mechanism of active milk protein in colostrum and mature milk and provide references for the identification and functional study of bioactive milk proteins in colostrum.

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.

Milk Lipid Regulation in Dairy Goats: A Comprehensive Review

The growth, development, and milk production traits of dairy goats, which are important sources of high-quality animal protein, are significantly influenced by a combination of genetic and environmental factors. It is imperative to identify key genetic loci that govern economically valuable traits in order to enhance breeding programs. Despite advancements in genomic technologies, there are still gaps in knowledge regarding the interplay between genetic factors and environmental influences, particularly in relation to the regulation of milk production and quality. Therefore, the aim of this paper was to synthesize advancements in the genetic and environmental factors affecting milk production and quality in dairy goats and identify key regulatory mechanisms. This review summarizes the recent progress on the identification of genes associated with milk production traits using whole-genome resequencing, the use of transcriptomic profiling to identify genes linked to milk production, the exploration of regulatory mechanisms of lipid metabolism in goat mammary epithelial cells, and the evaluation of the influence of nutritional factors on milk quality. A comprehensive understanding of these interactions is essential for enhancing breeding strategies and ensuring the sustainability of dairy goat farming. Future research should incorporate multi-omics approaches to unravel the intricate regulatory processes governing milk production and adapt practices to meet global demand while upholding economic and environmental sustainability.

Comparative expression analysis of water buffalo (Bubalus bubalis) to identify genes associated with economically important traits

The milk, meat, skins, and draft power of domestic water buffalo (Bubalus bubalis) provide substantial contributions to the global agricultural economy. The world's water buffalo population is primarily found in Asia, and the buffalo supports more people per capita than any other livestock species. For evaluating the workflow, output rate, and completeness of transcriptome assemblies within and between reference-free (RF) de novo transcriptome and reference-based (RB) datasets, abundant bioinformatics studies have been carried out to date. However, comprehensive documentation of the degree of consistency and variability of the data produced by comparing gene expression levels using these two separate techniques is lacking. In the present study, we assessed the variations in the number of differentially expressed genes (DEGs) attained with RF and RB approaches. In light of this, we conducted a study to identify, annotate, and analyze the genes associated with four economically important traits of buffalo, viz., milk volume, age at first calving, post-partum cyclicity, and feed conversion efficiency. A total of 14,201 and 279 DEGs were identified in RF and RB assemblies. Gene ontology (GO) terms associated with the identified genes were allocated to traits under study. Identified genes improve the knowledge of the underlying mechanism of trait expression in water buffalo which may support improved breeding plans for higher productivity. The empirical findings of this study using RNA-seq data-based assembly may improve the understanding of genetic diversity in relation to buffalo productivity and provide important contributions to answer biological issues regarding the transcriptome of non-model organisms.

Blood transcriptome reveals immune and metabolic-related genes involved in growth of pasteurized colostrum-fed calves

The quality of colostrum is a key factor contributing to healthy calf growth, and pasteurization of colostrum can effectively reduce the counts of pathogenic microorganisms present in the colostrum. Physiological changes in calves fed with pasteurized colostrum have been well characterized, but little is known about the underlying molecular mechanisms. In this study, key genes and functional pathways through which pasteurized colostrum affects calf growth were identified through whole blood RNA sequencing. Our results showed that calves in the pasteurized group (n = 16) had higher body height and daily weight gain than those in the unpasteurized group (n = 16) in all months tested. Importantly, significant differences in body height were observed at 3 and 4 months of age (p < 0.05), and in daily weight gain at 2, 3, and 6 months of age (p < 0.05) between the two groups. Based on whole blood transcriptome data from 6-months old calves, 630 differentially expressed genes (DEGs), of which 235 were upregulated and 395 downregulated, were identified in the pasteurized compared to the unpasteurized colostrum groups. Most of the DEGs have functions in the immune response (e.g., CCL3, CXCL3, and IL1A) and metabolism (e.g., PTX3 and EXTL1). Protein-protein interaction analyses of DEGs revealed three key subnetworks and fifteen core genes, including UBA52 and RPS28, that have roles in protein synthesis, oxidative phosphorylation, and inflammatory responses. Twelve co-expression modules were identified through weighted gene co-expression network analysis. Among them, 17 genes in the two modules that significantly associated with pasteurization were mainly involved in the tricarboxylic acid cycle, NF-kappa B signaling, and NOD-like receptor signaling pathways. Finally, DEGs that underwent alternative splicing in calves fed pasteurized colostrum have roles in the immune response (SLCO4A1, AKR1C4, and MED13L), indicative of potential roles in immune regulation. Results from multiple analytical methods used suggest that differences in calf growth between the pasteurized and unpasteurized groups may be due to differential immune activity. Our data provide new insights into the impact of pasteurization on calf immune and metabolic-related pathways through its effects on gene expression.

Feeding of palm oil fatty acids or rapeseed oil throughout lactation: Effects on mammary gene expression and milk production in Norwegian dairy goats

Lipid added as rapeseed or palm oil to the diet of dairy goats over 8 mo of one lactation alters fat secretion and milk fatty acid (FA) and protein composition. In this study, we examined the contribution of mammary gene expression to these changes and included 30 multiparous goats of Norwegian dairy goat breed for a 230-d experimental period, with indoor feeding from 1 to 120 d in milk (DIM), mountain grazing from 120 to 200 DIM, and indoor feeding from 200 to 230 DIM. After an initial period (1-60 DIM) when the control diet was given to all goats, the animals were subdivided into 3 groups of 10 goats. Treatments (60-230 DIM) were basal concentrate (control) alone or supplemented with either 8% (by weight) hydrogenated palm oil enriched with palmitic acid (POFA) or 8% (by weight) rapeseed oil (RSO). Milk was sampled individually from all animals throughout lactation, at 60, 120, 190, and 230 DIM for milk yield and composition. On d 60, 120, 190, and 230, mammary tissue was collected by biopsy to measure mRNA abundance of 19 key genes. None of the 19 genes involved in milk protein, apoptosis, lipid metabolism, transcription factors, and protein of the milk fat globule membrane, as measured by mRNA abundance, were affected by the lipid supplements, although POFA increased milk fat content, and POFA and RSO affected milk FA composition. Over the experimental period (120-230 DIM), the mRNA abundance of 13 of the 19 studied genes was affected by lactation stage. For some genes, expression either gradually increased from 120 to 230 DIM (CSN2,CASP8,CD36,GLUT4) or increased from 120 to 200 and then remained stable (XDH), or decreased (CSN3,G6PD,SREBF1,PPARG1) or increased only at 230 DIM (SCD1,SCD5,ELF3). For a second group of genes (CSN1, LALBA, FABP3, FASN, LPL, MFGE8), expression was stable over the lactation period. Our results suggest that factors other than gene expression, such as substrate availability or posttranscriptional regulation of these genes, could play an important role in the milk fat and FA responses to dietary fat composition in the goat. In conclusion, mammary gene expression in goats was more regulated by stage of lactation than by the dietary treatments applied.

Detection of candidate genes affecting milk production traits in sheep using whole-genome sequencing analysis

BACKGROUND

Artificial and natural selection for important economic traits and genetic adaptation of the populations to specific environments have led to the changes on the sheep genome. Recent advances in genome sequencing methods have made it possible to use comparative genomics tools to identify genes under selection for traits of economic interest in domestic animals.

OBJECTIVES

In this study, we compared the genomes of Assaf and Awassi sheep breeds with those of the Cambridge, Romanov and British du cher sheep breeds to explore positive selection signatures for milk traits using nucleotide diversity (Pi) and FST statistical methods.

METHODS

Genome sequences from fourteen sheep with a mean sequence depth of 9.32X per sample were analysed, and a total of 23 million single nucleotide polymorphisms (SNPs) were called and applied for this study. Genomic clustering of breeds was identified using ADMIXTURE software. The FST and Pi values for each SNP were computed between population A (Assaf and Awassi) and population B (Cambridge, British du cher, and Romanov).

RESULTS

The results of the PCA grouped two classes for these five dairy sheep breeds. The selection signatures analysis displayed 735 and 515 genes from FST and nucleotide diversity (Pi) statistical methods, respectively. Among all these, 12 genes were shared between the two approaches. The most conspicuous genes were related to milk traits, including ST3GAL1 (the synthesis of oligosacáridos), CSN1S1 (milk protein), CSN2 (milk protein), OSBPL8 (fatty acid traits), SLC35A3 (milk fat and protein percentage), VPS13B (total milk production, fat yield, and protein yield), DPY19L1 (peak yield), CCDC152 (lactation persistency and somatic cell count), NT5DC1 (lactation persistency), P4HTM (test day protein), CYTH4 (FAT Production) and METRNL (somatic cell), U1 (milk traits), U6 (milk traits) and 5S_RRNA (milk traits).

CONCLUSIONS

The findings provide new insight into the genetic basis of sheep milk properties and can play a role in designing sheep breeding programs incorporating genomic information.

Genome-wide association study on milk production and somatic cell score for Thai dairy cattle using weighted single-step approach with random regression test-day model

Genome-wide association studies are a powerful tool to identify genomic regions and variants associated with phenotypes. However, only limited mutual confirmation from different studies is available. The objectives of this study were to identify genomic regions as well as genes and pathways associated with the first-lactation milk, fat, protein, and total solid yields; fat, protein, and total solid percentage; and somatic cell score (SCS) in a Thai dairy cattle population. Effects of SNPs were estimated by a weighted single-step GWAS, which back-solved the genomic breeding values predicted using single-step genomic BLUP (ssGBLUP) fitting a single-trait random regression test-day model. Genomic regions that explained at least 0.5% of the total genetic variance were selected for further analyses of candidate genes. Despite the small number of genotyped animals, genomic predictions led to an improvement in the accuracy over the traditional BLUP. Genomic predictions using weighted ssGBLUP were slightly better than the ssGBLUP. The genomic regions associated with milk production traits contained 210 candidate genes on 19 chromosomes [Bos taurus autosome (BTA) 1 to 7, 9, 11 to 16, 20 to 21, 26 to 27 and 29], whereas 21 candidate genes on 3 chromosomes (BTA 11, 16, and 21) were associated with SCS. Many genomic regions explained a small fraction of the genetic variance, indicating polygenic inheritance of the studied traits. Several candidate genes coincided with previous reports for milk production traits in Holstein cattle, especially a large region of genes on BTA14. We identified 141 and 5 novel genes related to milk production and SCS, respectively. These novel genes were also found to be functionally related to heat tolerance (e.g., SLC45A2, IRAG1, and LOC101902172), longevity (e.g., SYT10 and LOC101903327), and fertility (e.g., PAG1). These findings may be attributed to indirect selection in our population. Identified biological networks including intracellular cell transportation and protein catabolism implicate milk production, whereas the immunological pathways such as lymphocyte activation are closely related to SCS. Further studies are required to validate our findings before exploiting them in genomic selection.

Characterization and comparative analysis of transcriptional profiles of porcine colostrum and mature milk at different parities

BACKGROUND

Porcine milk is a complex fluid, containing a myriad of immunological, biochemical, and cellular components, made to satisfy the nutritional requirements of the neonate. Whole milk contains many different cell types, including mammary epithelial cells, neutrophils, macrophages, and lymphocytes, as well nanoparticles, such as milk exosomes. To-date, only a limited number of livestock transcriptomic studies have reported sequencing of milk. Moreover, those studies focused only on sequencing somatic cells as a proxy for the mammary gland with the goal of investigating differences in the lactation process. Recent studies have indicated that RNA originating from multiple cell types present in milk can withstand harsh environments, such as the digestive system, and transmit regulatory molecules from maternal to neonate. Transcriptomic profiling of porcine whole milk, which is reflective of the combined cell populations, could help elucidate these mechanisms. To this end, total RNA from colostrum and mature milk samples were sequenced from 65 sows at differing parities. A stringent bioinformatic pipeline was used to identify and characterize 70,841 transcripts.

RESULTS

The 70,841 identified transcripts included 42,733 previously annotated transcripts and 28,108 novel transcripts. Differential gene expression analysis was conducted using a generalized linear model coupled with the Lancaster method for P-value aggregation across transcripts. In total, 1667 differentially expressed genes (DEG) were identified for the milk type main effect, and 33 DEG were identified for the milk type x parity interaction. Several gene ontology (GO) terms related to immune response were significant for the milk type main effect, supporting the well-known fact that immunoglobulins and immune cells are transferred to the neonate via colostrum.

CONCLUSIONS

This is the first study to perform global transcriptome analysis from whole milk samples in sows from different parities. Our results provide important information and insight into synthesis of milk proteins and innate immunity and potential targets for future improvement of swine lactation and piglet development.

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.

Expression patterns of major genes in fatty acid synthesis, inflammation, oxidative stress pathways from colostrum to milk in Damascus goats

The molecular regulation of milk secretion and quality in the transition period from colostrum to milk in goats is largely unknown. In the present study, mammary gland secretion of goats was collected in 0th, 4th, 7th, 14th and 28th days after parturition. In addition to composition and fatty acid profile of colostrum or milk, FASN, SCD, ACACA, COX-2, NRF2, TLR2, NF-kB, LTF and PTX3 genes expression patterns were determined from milk somatic cells. While somatic cell count (SCC), malondialdehyde (MDA), fat, fat-free dry matter, protein and lactose were highest as expression levels of the oxidative and inflammatory genes, freezing point and electrical conductivity were lowest in colostrum. With the continuation of lactation, most of the fatty acids, n3 ratio, and odour index increased but C14:0 and C16:0 decreased. While FASN was upregulated almost threefolds in 14th day, ACACA was upregulated more than fivefolds in 7th and 14th days. Separately, the major genes in fatty acid synthesis, inflammation and oxidative stress were significantly associated with each other due to being positively correlated. MDA was positively correlated with SCC and some of the genes related inflammation and oxidative stress. Furthermore, significant negative correlations were determined between SCC and fatty acid synthesis related genes. With this study, transition period of mammary secretion was particularly clarified at the molecular levels in Damascus goats.

Genome-Wide Association Study Identifies New Candidate Markers for Somatic Cells Score in a Local Dairy Sheep

In the Mediterranean basin countries, the dairy sheep production is usually based on local breeds, which are very well-adapted to their production systems and environments and can indeed guarantee income, employment, and economic viability in areas where production alternatives are scarce or non-existent. Mastitis is still one of the greatest problems affecting commercial milk production. However, genetic evaluation of mastitis is particularly difficult because of its low heritability and the categorical nature of the trait. The aim of this study was to identify genomic regions putatively associated with somatic cells count (SCC) in the local economically important Valle del Belice sheep breed using of deregressed breeding values (DEBV) as response variables. All the samples were genotyped using the Illumina OvineSNP50K BeadChip. Genome-wide association analysis was carried out based on regression of DEBV. A total of eight markers were found to be significantly associated with log-transformed SCC. Several candidate genes associated with SCC were identified related to immunity system and udder conformation. The results can help improving the competitiveness of the local Valle del Belìce breed. Further studies considering a higher sample size or independent population will be needed to confirm our results.

Genomic studies of milk-related traits in water buffalo (Bubalus bubalis) based on single-step genomic best linear unbiased prediction and random regression models

Genomic selection has been widely implemented in many livestock breeding programs, but it remains incipient in buffalo. Therefore, this study aimed to (1) estimate variance components incorporating genomic information in Murrah buffalo; (2) evaluate the performance of genomic prediction for milk-related traits using single- and multitrait random regression models (RRM) and the single-step genomic best linear unbiased prediction approach; and (3) estimate longitudinal SNP effects and candidate genes potentially associated with time-dependent variation in milk, fat, and protein yields, as well as somatic cell score (SCS) in multiple parities. The data used to estimate the genetic parameters consisted of a total of 323,140 test-day records. The average daily heritability estimates were moderate (0.35 ± 0.02 for milk yield, 0.22 ± 0.03 for fat yield, 0.42 ± 0.03 for protein yield, and 0.16 ± 0.03 for SCS). The highest heritability estimates, considering all traits studied, were observed between 20 and 280 d in milk (DIM). The genetic correlation estimates at different DIM among the evaluated traits ranged from -0.10 (156 to 185 DIM for SCS) to 0.61 (36 to 65 DIM for fat yield). In general, direct selection for any of the traits evaluated is expected to result in indirect genetic gains for milk yield, fat yield, and protein yield but also increase SCS at certain lactation stages, which is undesirable. The predicted RRM coefficients were used to derive the genomic estimated breeding values (GEBV) for each time point (from 5 to 305 DIM). In general, the tuning parameters evaluated when constructing the hybrid genomic relationship matrices had a small effect on the GEBV accuracy and a greater effect on the bias estimates. The SNP solutions were back-solved from the GEBV predicted from the Legendre random regression coefficients, which were then used to estimate the longitudinal SNP effects (from 5 to 305 DIM). The daily SNP effect for 3 different lactation stages were performed considering 3 different lactation stages for each trait and parity: from 5 to 70, from 71 to 150, and from 151 to 305 DIM. Important genomic regions related to the analyzed traits and parities that explain more than 0.50% of the total additive genetic variance were selected for further analyses of candidate genes. In general, similar potential candidate genes were found between traits, but our results suggest evidence of differential sets of candidate genes underlying the phenotypic expression of the traits across parities. These results contribute to a better understanding of the genetic architecture of milk production traits in dairy buffalo and reinforce the relevance of incorporating genomic information to genetically evaluate longitudinal traits in dairy buffalo. Furthermore, the candidate genes identified can be used as target genes in future functional genomics studies.

Transcriptome study underling difference of milk yield during peak lactation of Kazakh horse

This study was designed to provide a basis for further understanding of the mechanism of lactation based on mRNA expression differences in milk fat between different milk yields in Kazakh horses. Total RNA was extracted from the milk fat during the peak of lactation period. A total of 310 differentially expressed genes (DEGs) were identified by comparative transcriptome analysis of the high-yield and low-yield group. These DEGs regulate lactation by participated in AMPK signaling pathway, FoxO signaling pathway, ErbB signaling pathway, VEGF signaling pathway. In addition, we performed quantitative PCR to validated 5 selected DEGs and the results were in agreement with RNA-seq analysis. A new profile has been established for revealing the mechanism of equid's mammalian lactation.

Weighted gene co-expression network analysis identifies modules and functionally enriched pathways in the lactation process

The exponential growth in knowledge has resulted in a better understanding of the lactation process in a wide variety of animals. However, the underlying genetic mechanisms are not yet clearly known. In order to identify the mechanisms involved in the lactation process, various mehods, including meta-analysis, weighted gene co-express network analysis (WGCNA), hub genes identification, gene ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment at before peak (BP), peak (P), and after peak (AP) stages of the lactation processes have been employed. A total of 104, 85, and 26 differentially expressed genes were identified based on PB vs. P, BP vs. AP, and P vs. AP comparisons, respectively. GO and KEGG pathway enrichment analysis revealed that DEGs were significantly enriched in the "ubiquitin-dependent ERAD" and the "chaperone cofactor-dependent protein refolding" in BP vs. P and P vs. P, respectively. WGCNA identified five significant functional modules related to the lactation process. Moreover, GJA1, AP2A2, and NPAS3 were defined as hub genes in the identified modules, highlighting the importance of their regulatory impacts on the lactation process. The findings of this study provide new insights into the complex regulatory networks of the lactation process at three distinct stages, while suggesting several candidate genes that may be useful for future animal breeding programs. Furthermore, this study supports the notion that in combination with a meta-analysis, the WGCNA represents an opportunity to achieve a higher resolution analysis that can better predict the most important functional genes that might provide a more robust bio-signature for phenotypic traits, thus providing more suitable biomarker candidates for future studies.

Aging-Associated Alterations in Mammary Epithelia and Stroma Revealed by Single-Cell RNA Sequencing

Aging is closely associated with increased susceptibility to breast cancer, yet there have been limited systematic studies of aging-induced alterations in the mammary gland. Here, we leverage high-throughput single-cell RNA sequencing to generate a detailed transcriptomic atlas of young and aged murine mammary tissues. By analyzing epithelial, stromal, and immune cells, we identify age-dependent alterations in cell proportions and gene expression, providing evidence that suggests alveolar maturation and physiological decline. The analysis also uncovers potential pro-tumorigenic mechanisms coupled to the age-associated loss of tumor suppressor function and change in microenvironment. In addition, we identify a rare, age-dependent luminal population co-expressing hormone-sensing and secretory-alveolar lineage markers, as well as two macrophage populations expressing distinct gene signatures, underscoring the complex heterogeneity of the mammary epithelia and stroma. Collectively, this rich single-cell atlas reveals the effects of aging on mammary physiology and can serve as a useful resource for understanding aging-associated cancer risk.

Using single-cell RNA-sequencing, Li et al. compare mammary epithelia and stroma in young and aged mice. Age-dependent changes at cell and gene levels provide evidence suggesting alveolar maturation, functional deterioration, and potential pro-tumorigenic and inflammatory alterations. Additionally, identification of heterogeneous luminal and macrophage subpopulations underscores the complexity of mammary lineages.

Goat Milk Oligosaccharides: Their Diversity, Quantity, and Functional Properties in Comparison to Human Milk Oligosaccharides

Human milk is considered the golden standard in infant nutrition. Free oligosaccharides in human milk provide important health benefits. These oligosaccharides function as prebiotics, immune modulators, and pathogen inhibitors and were found to improve barrier function in the gut. Infant formulas nowadays often contain prebiotics but lack the specific functions of human milk oligosaccharides (hMOS). Milk from domesticated animals also contains milk oligosaccharides but at much lower levels and with less diversity. Goat milk contains significantly more oligosaccharides (gMOS) than bovine (bMOS) or sheep (sMOS) milk and also has a larger diversity of structures. This review summarizes structural studies, revealing a diversity of up to 77 annotated gMOS structures with almost 40 structures fully characterized. Quantitative studies of goat milk oligosaccharides range from 60 to 350 mg/L in mature milk and from 200 to 650 mg/L in colostrum. These levels are clearly lower than in human milk (5-20 g/L) but higher than in other domesticated dairy animals, e.g., bovine (30-60 mg/L) and sheep (20-40 mg/L). Finally, the review focuses on demonstrated and potential functionalities of gMOS. Some studies have shown anti-inflammatory effects of mixtures enriched in gMOS. Goat MOS also display prebiotic potential, particularly in stimulating growth of bifidobacteria preferentially. Although functional studies of gMOS are still limited, several structures are also found in human milk and have known functions as immune modulators and pathogen inhibitors. In conclusion, goat milk constitutes a promising alternative source for milk oligosaccharides, which can be used in infant formula.

Milk somatic cell derived transcriptome analysis identifies regulatory genes and pathways during lactation in Indian Sahiwal cattle (Bos indicus)

BACKGROUND

The present study is an effort to understand the genomic drivers of lactation in Sahiwal (Bos indicus), the best milch cattle breed of the tropics.

METHODS

RNA sequencing of four animals from early, mid and late lactation stages was performed using milk somatic cells as source of RNA.

RESULTS

The genes encoding the milk casein and whey proteins showed highest expression in early and mid lactation, with a declining trend towards the late stage. The enhanced expression of PLIN2, FABP5 and FABP3 genes in mid lactation suggests enrichment of the PPARα pathway which is linked to fatty acid metabolism. A gradual decline in the percentage of genes involved in metabolism of proteins, mRNA and insulin synthesis from early to late lactation reflected transition from lactogenesis to involution. Major biological pathways maintained throughout lactation were adaptive immune system, FGF signaling, EGFR signaling, activated TLR4 signaling, NFkB and MAP kinases activation mediated by TLR4 signaling repertoire. Differential expression analysis revealed 547, 1010 and 1313 differentially expressed genes (p < 0.05) between early-late, early-mid and mid-late stages, respectively. The topmost regulatory genes identified by network analysis from the differentially expressed genes, were involved in Chemokine receptor, GPCR and EGFR1 pathways.

CONCLUSION

The genes and pathways delineated in this study have regulatory implications in cell morphogenesis, lipid droplet formation and protein synthesis in the course of lactation. The study provides an insight into the expression profile of genes influencing milk properties and lactation in Sahiwal cattle.

Transcriptome signature of two lactation stages in Ghezel sheep identifies using RNA-Sequencing

The expression of genes and their regulation during lactation in Ghezel sheep breed remains less understood. To explore the underlying molecular mechanism of the lactation process in the mammary gland, transcriptome profiles of Iranian fat-tailed Ghezel sheep breed milk at two stages, before (BF) and after peak (AF) stages of lactation were investigated. Functional impacts of differentially expressed genes (DEGs) between BF and AF stages were surveyed using Gene Ontology (GO) and Protein-Protein Interaction (PPI) network analysis. Totally, 75 DEGs were identified between BF and AF stages of lactation. The RNA-Seq results were validated by Q-RT-PCR. Gene ontology of DEGs mainly enriched in metabolic process and oxidative phosphorylation. PPI network analysis also highlighted the contribution of peroxisome proliferator-activated receptors (PPAR) signaling, oxidative phosphorylation and metabolic pathways in the lactation process. Intriguingly, the genes involved in fat metabolism dominantly down-regulated at AF stage. Our results provide new insight into transcriptional changes and add to growing body of knowledge on the lactation process in fat-tailed sheep breeds.

Characterization of microRNA profiles in the mammary gland tissue of dairy goats at the late lactation, dry period and late gestation stages

MicroRNAs (miRNAs) play an important role in regulating mammary gland development and lactation. We previously analyzed miRNA expression profiles in Laoshan dairy goat mammary glands at the early (20 d postpartum), peak (90 d postpartum) and late lactation (210 d postpartum) stages. To further enrich and clarify the miRNA expression profiles during the lactation physiological cycle, we sequenced miRNAs in the mammary gland tissues of Laoshan dairy goats at three newly selected stages: the late lactation (240 d postpartum), dry period (300 d postpartum) and late gestation (140 d after mating) stages. We obtained 4038 miRNAs and 385 important miRNA families, including mir-10, let-7 and mir-9. We also identified 754 differentially expressed miRNAs in the mammary gland tissue at the 3 different stages and 6 groups of miRNA clusters that had unique expression patterns. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses showed that GO terms such as mammary gland development (GO:0030879) and mammary gland morphogenesis (GO:0060443) and important signaling pathways, including the insulin signaling pathway (chx04910), hippo signaling pathway (chx04390) and estrogen signaling pathway (chx04915), were enriched. We screened miRNAs and potential target genes that may be involved in the regulation of lactation, mammary gland growth and differentiation, cell apoptosis, and substance transport and synthesis and detected the expression patterns of important genes at the three stages. These miRNAs and critical target genes may be important factors for mammary gland development and lactation regulation and potentially valuable molecular markers, which may provide a theoretical reference for further investigation of mammary gland physiology.

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.

Genomic Identification, Evolution, and Expression Analysis of Collagen Genes Family in Water Buffalo during Lactation

Collagens, as extracellular matrix proteins, support cells for structural integrity and contribute to support mammary basic structure and development. This study aims to perform the genomic identification, evolution, and expression analyses of the collagen gene family in water buffalo (Bubalus bubalis) during lactation. A total of 128 buffalo collagen protein sequences were deduced from the 45 collagen genes identified in silico from buffalo genome, which classified into six groups based on their phylogenetic relationships, conserved motifs, and gene structure analyses. The identified collagen sequences were unequally distributed on 16 chromosomes. The tandem duplicated genes were found within three chromosomes, while only one segmental event occurred between Chr3 and Chr8. Collinearity analysis revealed that a total of 36 collagen gene pairs were orthologous between buffalo and cattle genomes despite having different chromosome numbers. Comparative transcription analyses revealed that a total of 23 orthologous collagen genes were detected in the milk samples at different lactation periods between the two species. Notably, the duplicated gene pair of COL4A1-COL4A2 during lactation had a higher mRNA expression level than that of cattle, while a higher expression level of COL6A1-COL6A2 pair was found in cattle compared with that of buffalo. The present study provides useful information for investigating the potential functions of the collagen family in buffalo during lactation and helps in the functional characterization of collagen genes in additional research.

Mammary Transcriptome Profile during Peak and Late Lactation Reveals Differentially Expression Genes Related to Inflammation and Immunity in Chinese Holstein

Simple Summary

Milk somatic cell count, referring to the total number of somatic cells per milliliter of bovine milk, changes regularly during the lactation cycle. The somatic cell count of healthy cows is usually higher in late lactation than in peak lactation. When the inflammatory response in dairy cow mammary gland becomes more intense, the milk somatic cell count increases together with the reduction of milk quality and yield. Autoimmunity was thought to play an important role in the prevention of mastitis in late lactation of dairy cattle. However, the underlying mechanisms related to the gene expression levels during the process remain unknown. In this study, transcriptome sequencing was performed to screen the differentially expressed genes related to the inflammation and immunity in healthy Chinese Holstein mammary glands. Our findings are helpful to understand the physiological functions of mammary inflammation of Chinese Holstein during late lactation.

Abstract

Somatic cell count (SCC) in milk is widely used in the dairy industry, as an indicator of the health of mammary gland. While the SCC of dairy cattle was higher in late lactation than in peak lactation, its association with gene expressions of mammary gland were largely unknown. In this study, a transcriptomic sequencing approach and bioinformatics analysis were used to investigate the differential expressed genes (DEGs) associated with inflammation and immunity between peak and late periods of lactation in Chinese Holstein. A total of 446 DEGs (padj < 0.05 and fold change >2) were identified, 50 of which belonged to seven pathways and five terms related to inflammation and immunity. Our data suggested that the activation of nuclear transcription factor-κB (NF-κB) pathway and Toll-like receptor signaling pathway caused inflammatory response, and the activation of chemokine signaling pathway and cytokine–cytokine receptor interaction signaling pathway caused a protective immune response to ensure dairy cows health during late lactation. Our findings deepen the understanding of the molecular mechanism and physiological functions of mammary inflammation in Chinese Holstein during late 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.

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.

Dynamic Temporal Variations in Bovine Lactoferrin Glycan Structures

It has been reported previously that glycosylation of bovine lactoferrin changes over time. A detailed structural overview of these changes over the whole course of lactation, including predry period milk, is lacking. In this study, a high-throughput analysis method was applied to the glycoprofile of lactoferrin isolated from colostrum, mature, and predry period mature milk, which was analyzed over two subsequent lactation cycles for 8 cows from diverse genetic backgrounds. In addition, comparisons are made with commercial bovine lactoferrin samples. During the first 72 h, dynamic changes in lactoferrin glycosylation occurred. Shifts in the oligomannose distribution and the number of sialylated and fucosylated glycans were observed. In some cows, we observed (α2,3)-linked sialic acid in the earliest colostrum samples. The glycoprofiles appeared stable from 1 month after delivery, as well as between cows. In addition, the glycosylation profiles of commercial lactoferrins isolated from pooled mature milk were stable over the year. Lactoferrin glycosylation in the predry period resembles colostrum lactoferrin. The variations in lactoferrin glycosylation profiles, lactoferrin concentrations, and other milk parameters provide detailed information that potentially assists in unraveling the functions and biosynthesis regulation of lactoferrin glycosylation.

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.

Transcriptome Profile Analysis of Mammary Gland Tissue from Two Breeds of Lactating Sheep

The mammary gland is a crucial tissue for milk synthesis and plays a critical role in the feeding and growth of mammalian offspring. The aim of this study was to use RNA-sequencing (RNA-Seq) technology to provide a transcriptome profile of the ovine mammary gland at the peak of lactation. Small-Tailed Han (STH) sheep (n = 9) and Gansu Alpine Merino (GAM) sheep (n = 9), breeds with phenotypic differences in milk production traits, were selected for the RNA-Seq analysis. This revealed 74 genes that were more highly expressed in the STHs than in the GAMs. Similarly, 143 genes that were expressed at lower levels in the STHs than in the GAMs, were identified. Gene ontogeny (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses revealed that these differentially expressed genes (DEGs) were associated with binding and catalytic activities, hematopoietic cell lineages, oxytocin signaling pathway and neuroactive ligand–receptor interaction. This is the first study of the transcriptome profile of the ovine mammary gland in these Chinese breeds at peak lactation. The results provide for a better understanding of the genetic mechanisms involved in ovine lactation.

Comparative Transcriptomic and Proteomic Analyses Identify Key Genes Associated With Milk Fat Traits in Chinese Holstein Cows

Milk fat is the most important energy substance in milk and contributes to its quality and health benefits. However, the genetic mechanisms underlying milk fat synthesis are not fully understood. The development of RNA sequencing and tandem mass tag technologies has facilitated the identification of eukaryotic genes associated with complex traits. In this study, we used these methods to obtain liver transcriptomic and proteomic profiles of Chinese Holstein cows (n = 6). Comparative analyses of cows with extremely high vs. low milk fat percentage phenotypes yielded 321 differentially expressed genes (DEGs) and 76 differentially expressed proteins (DEPs). Functional annotation of these DEGs and DEPs revealed 26 genes that were predicted to influence lipid metabolism through insulin, phosphatidylinositol 3-kinase/Akt, mitogen-activated protein kinase, 5′ AMP-activated protein kinase, mammalian target of rapamycin, and peroxisome proliferator-activated receptor signaling pathways; these genes are considered as the most promising candidate regulators of milk fat synthesis. The findings of this study enhance the understanding of the genetic basis and molecular mechanisms of milk fat synthesis, which could lead to the development of cow breeds that produce milk with higher nutritional value.

Buffalo milk transcriptome: A comparative analysis of early, mid and late lactation

The expression of genes and their regulation during lactation in buffaloes remains less understood. To understand the interplay of various genes and pathways, the milk transcriptome from three lactation stages of Murrah buffalo was analyzed by RNA sequencing. The filtered reads were mapped to the Bubalus bubalis as well as Bos taurus reference assemblies. The average mapping rate to water buffalo and Btau 4.6 reference sequence, was 75.5% and 75.7% respectively. Highly expressed genes (RPKM > 3000), throughout lactation included CSN2, CSN1S1, CSN3, LALBA, SPP1 and TPT1. A total of 12833 transcripts were common across all the stages, while 271, 205 and 418 were unique to early, mid and late lactation respectively. Majority of the genes throughout lactation were linked to biological functions like protein metabolism, transport and immune response. A discernible shift from metabolism in early stage to metabolism and immune response in mid stage, and an increase in immune response functions in late lactation was observed. The results provide information of candidate genes and pathways involved in the different stages of lactation in buffalo. The study also identified 14 differentially expressed and highly connected genes across the three lactation stages, which can be used as candidates for future research.

Gelsolin expression in sheep milk somatic cells during lactation

The identification of genes involved in phenotypes related to milk quality is important for both economic and health aspects in livestock production. The aim of this study was to assess the level of gelsolin gene expression in two breeds of dairy sheep - Sarda and Gentile - with pronounced differences in quantitative and qualitative milk traits. Gelsolin, a type of actin-modulating proteins is involved in the processes of actin remodeling during cell growth and apoptosis; therefore a role of this protein in mammary changes during lactation was here hypothesized. Individual milk samples were collected three times during lactation from 26 ewes of the two breeds. The differential gene expression of gelsolin in the two breeds and the three lactation times was estimated by quantitative PCR on RNA extracted from milk somatic cells. Correlations of gelsolin gene expression with milk yield and quality and days of lactation were also estimated. The results showed that gelsolin gene expression was significantly higher in the Sarda compared to the Gentile at each lactation stage, in agreement with the longer lactation duration and the higher daily milk yield of the first breed. Significant correlations of gelsolin gene expression were found with milk fat content in Sarda breed (-0.46, P<0.05). Gelsolin expression analysis confirmed the link between gelsolin gene function and milk fat content of sheep.

Transcriptional profiling of two different physiological states of the yak mammary gland using RNA sequencing

Yak milk is superior to common cow milk in nutrients including protein, fat and calories. However, the milk yield of the yak is very much lower compared with other dairy bovines. To understand the molecular mechanisms of lactogenesis, lactation and mammary gland development, mammary tissue samples were taken from five yaks during a dry period (DP, n = 3) and lactation period (LP, n = 2). Two types of cDNA sequence libraries that reflected the different physiological states of the mammary gland were constructed using RNA sequencing technology. After removing reads containing adapters, reads containing poly-N and low-quality reads from the raw data, 45,423,478 to 53,274,976 clean reads were obtained from these libraries. A total of 74.72% to 80.65% of the high-quality sequence reads were uniquely aligned to the BosGru v2.0 yak reference genome. Using the DESeq R package, 360 differentially expressed genes were detected between the two groups when the adjusted P value (padj < 0.05) was used as the cutoff value; this included 192 upregulated and 168 downregulated genes in the yak mammary gland tissue of the DP compared to the LP. A gene ontology analysis revealed that the most enriched GO terms were protein binding, multi-organism process, immune system and others. KEGG pathway analysis indicated that the differentially expressed genes were mostly enriched in Hippo signaling, insulin signaling, steroid biosynthesis and others. The analysis of the up- and downregulated genes provides important insights into the molecular events involved in lactogenesis, lactation and mammary gland development and will guide further research to enhance milk yield and optimize the constituents of yak milk.

Transcriptional profiling of swine mammary gland during the transition from colostrogenesis to lactogenesis using RNA sequencing

Background

Colostrum and milk are essential sources of antibodies and nutrients for the neonate, playing a key role in their survival and growth. Slight abnormalities in the timing of colostrogenesis/lactogenesis potentially threaten piglet survival. To further delineate the genes and transcription regulators implicated in the control of the transition from colostrogenesis to lactogenesis, we applied RNA-seq analysis of swine mammary gland tissue from late-gestation to farrowing. Three 2nd parity sows were used for mammary tissue biopsies on days 14, 10, 6 and 2 before (−) parturition and on day 1 after (+) parturition. A total of 15 mRNA libraries were sequenced on a HiSeq2500 (Illumina Inc.). The Dynamic Impact Approach and the Ingenuity Pathway Analysis were used for pathway analysis and gene network analysis, respectively.

Results

A large number of differentially expressed genes were detected very close to parturition (−2d) and at farrowing (+ 1d). The results reflect the extraordinary metabolic changes in the swine mammary gland once it enters into the crucial phases of lactogenesis and underscore a strong transcriptional component in the control of colostrogenesis. There was marked upregulation of genes involved in synthesis of colostrum and main milk components (i.e. proteins, fat, lactose and antimicrobial factors) with a pivotal role of CSN1S2, LALBA, WAP, SAA2, and BTN1A1. The sustained activation of transcription regulators such as SREBP1 and XBP1 suggested they help coordinate these adaptations.

Conclusions

Overall, the precise timing for the transition from colostrogenesis to lactogenesis in swine mammary gland remains uncharacterized. However, our transcriptomic data support the hypothesis that the transition occurs before parturition. This is likely attributable to upregulation of a wide array of genes including those involved in ‘Protein and Carbohydrate Metabolism’, ‘Immune System’, ‘Lipid Metabolism’, ‘PPAR signaling pathway’ and ‘Prolactin signaling pathway’ along with the activation of transcription regulators controlling lipid synthesis and endoplasmic reticulum biogenesis and stress response.

Electronic supplementary material

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

Detection of genetic variants between different Polish Landrace and Puławska pigs by means of RNA-seq analysis

Variant calling analysis based on RNA sequencing data provides information about gene variants. RNA-seq is cheaper and faster than is DNA sequencing. However, it requires individual hard filters during data processing due to post-transcriptional modifications such as splicing and RNA editing. In the present study, RNA-seq transcriptome data on two Polish pig breeds (Puławska, PUL, n = 8, and Polish Landrace, PL, n = 8) were included. The pig breeds are significantly different with regard to meat qualities such as texture, water exudation, growth traits and fat content in carcasses. A total of 2451 significant mutations were identified by a chi square tests, and functional analysis was carried out using Panther, KEGG and Kobas. Interesting missense gene variants and mutations located in regulatory regions were found in a few genes related to fatty acid metabolism and lipid storage such as ACSL5, ALDH3A2, FADS1, SCD, PLA2G12A and ATGL. A validation of mutational influences on pig traits was performed for ALDH3A2, ATGL, PLA2G12A and MYOM1 variants using association analysis including 215 pigs of the PL and PUL breeds. The ALDH3A2ENSSSCT00000019636.2:c.470T>C polymorphism was found to affect the weight of the ham and loin eye area. In turn, an ENSSSCT00000004091.2:c.2836G>A MYOM1 mutation, which could be implicated in myofibrillar network organisation, had an effect on meatiness and loin texture parameters. The study aimed to estimate the usefulness of RNA-seq results for a purpose other than differentially expressed gene analysis. The analysis performed indicated interesting gene variants that could be used in the future as markers during selection.

RNA sequencing and differential expression reveals the effects of serial oestrus synchronisation on ovarian genes in dairy goats

A total of 24 female Xinong Saanen dairy goats were used to examine differentially expressed genes (DEGs) in the ovaries of goats treated once or three times for oestrus synchronisation (ES). The goats were randomly divided into two groups: one group received three ES treatments at fortnightly intervals (repeated or triple ES group), whereas the other was only treated once on the same day as the third ES treatment for the triple group (control group) during the breeding season. Ovaries of three goats in oestrus from each group were collected for morphological examination and transcriptome sequencing, while the rest of the goats were artificially inseminated twice. Litter size and fecundity rate tended (P = 0.06) to be lower in the triple ES group. A total of 319 DEGs were identified, including carbohydrate sulphotransferase 8 (CHST8), corticosteroid-binding globulin (CBG), oestradiol 17-β-dehydrogenase 1 (DHB1), oestrogen receptor 1 (ESR1), progestin and adipoQ receptor family member 4 (PAQR4), PAQR9, prostacyclin synthase (PTGIS), contactin-associated protein (CNTNAP4), matrix metalloproteinase-2 (MMP-2), regulator of G-protein signalling 9-2 (RGS9-2) and sperm surface protein Sp17 (Sp17); these were the most promising novel candidate genes for reproductive performances in goats. Our study indicates that triple ES could cause DNA damage and alter gene expression in goat ovaries, potentially affecting ovary function, neural regulation and hormone secretion.

Effect of high-fat diet on secreted milk transcriptome in midlactation mice

High-fat diet (HFD) during lactation alters milk composition and is associated with development of metabolic diseases in the offspring. We hypothesized that HFD affects milk microRNA (miRNA) and mRNA content, which potentially impact offspring development. Our objective was to determine the effect of maternal HFD on secreted milk transcriptome. To meet this objective, 4 wk old female ICR mice were divided into two treatments: control diet containing 10% kcal fat and HFD containing 60% kcal fat. After 4 wk on CD or HFD, mice were bred while continuously fed the same diets. On postnatal day 2 (P2), litters were normalized to 10 pups, and half the pups in each litter were cross-fostered between treatments. Milk was collected from dams on P10 and P12. Total RNA was isolated from milk fat fraction of P10 samples and used for mRNA-Seq and small RNA-Seq. P12 milk was used to determine macronutrient composition. After 4 wk of prepregnancy feeding HFD mice weighed significantly more than did the control mice. Lactose and fat concentration were significantly ( P < 0.05) higher in milk of HFD dams. Pup weight was significantly greater ( P < 0.05) in groups suckled by HFD vs. control dams. There were 25 miRNA and over 1,500 mRNA differentially expressed (DE) in milk of HFD vs. control dams. DE mRNA and target genes of DE miRNA enriched categories that were primarily related to multicellular organismal development. Maternal HFD impacts mRNA and miRNA content of milk, if bioactive nucleic acids are absorbed by neonate differences may affect development.

LncRNA as ceRNAs may be involved in lactation process

The main function of the mammary gland is to secret milk for newborn growth. Milk production process is regulated by hormones, growth factors, noncoding RNAs and other factors locally. Long non-coding RNAs (lncRNAs), one type of recently discovered non-coding RNA, have been found in mammary gland and some studies suggested lncRNA may play important roles in mammary gland development. Competing endogenous RNAs (ceRNAs) are emerging to compete for miRNA binding and, in turn, regulate each other. In the current study, we sequenced mRNA, miRNA and lncRNA in goat mammary tissue at 2 points in lactation (early and mature). All data were co-expressed together from the same samples. Our data showed that the ceRNAs up-regulated during the mature lactation phase were associated with lipid, protein, carbon and amino acid synthesis and metabolism. This correlates with the function of the mature lactation phase: i.e. the continuous production of large amounts of milk, rich in proteins, lipids, amino acids and other nutrients. Alternately, the ceRNAs up-regulated during early lactation were associated with PI3K-AKT pathways and ECM-receptor interactions; these fulfil the functional role of preparing the mammary gland for full lactation. Therefore, the results suggest that ceRNAs work synergistically during different developmental stages to regulate specific functions associated with lactation control. This study suggests that ceRNAs (lncRNA-mRNA) may be involved in lactation process.

Changes in Caprine Milk Oligosaccharides at Different Lactation Stages Analyzed by High Performance Liquid Chromatography Coupled to Mass Spectrometry

Changes of the abundance of caprine milk oligosaccharides (CMO) at different lactation stages have been evaluated by hydrophilic interaction liquid chromatography coupled to mass spectrometry (HILIC-Q MS) and nanoflow liquid chromatography-quadrupole-time-of-flight mass spectrometry (nano-LC-Chip-QTOF MS). Eight major oligosaccharides (OS) were quantified at different lactation stages by HILIC-Q MS, while the use of nano-LC-Chip-QToF MS allowed expanding the study to forty-nine different OS by monitoring neutral non- and fucosylated species, as well as acidic species containing not only N-acetyl-neuraminic acid or N-glycolyl-neuraminic acid residues but also the combination of both sialic acids. Overall, the most abundant OS decreased with lactation time, whereas different trends were observed for minor OS. 6'-Sialyl-lactose was the most abundant acidic OS while galactosyl-lactose isomers were identified as the most abundant neutral OS. This is the first time that a comprehensive study regarding the changes of the abundance of CMO, both neutral and acidic, at different lactation stages is carried out.