Functional development of the adult ovine mammary gland--insights from gene expression profiling

Transcriptome analysis of mammary epithelial cell between Sewa sheep and East FriEsian sheep from different localities

Mammary epithelial cells, the only milk-producing cell type in the mammary gland, undergo dynamic proliferation and differentiation during pregnancy, culminating in lactation postpartum. The East FriEsian sheep ranks among the world's most prolific dairy breeds, while the Sewa sheep, a unique dual-purpose breed autochthonous to the Qinghai-Tibet Plateau, exhibits significantly lower milk production. Employing tissue culture methods, we successfully established mammary epithelial cell lines from both breeds. Morphological assessment of mammary epithelial cells and immunofluorescence identification of Cytokeratin 7 and Cytokeratin 8 confirmed the epithelial identity of the isolated cells. Subsequent RNA-seq analysis of these in vitro epithelial cell lines revealed 1813 differentially expressed genes (DEGs). Among these, 1108 were significantly up-regulated and 705 were down-regulated in Sewa epithelial sheep cells compared to East FriEsian epithelial cells. KEGG enrichment analysis identified cellular processes, environmental information processing, human diseases, metabolism, and organismal systems as the primary functional categories associated with DEGs. Gene ontology (GO) terms annotation, categorized into molecular function, biological processes, and cellular component, yielded "binding and catalytic activity," "molecular function regulator activity," and "cellular process," "biological regulation," and "regulation of biological process" as the top three terms within each domain, respectively. Clusters of Orthologous Groups of proteins (KOG) classification further revealed that "signal transduction mechanisms" accounted for the largest proportion of DEGs among all KOG categories. Finally, based on these analyses, ATF3 and MPP7 were identified as promising candidate genes for regulating lactation.

Melatonin promotes the proliferation of dermal papilla cells in cashmere goats via activation of chi-let-7d-5p/WNT2 axis

Exogenous melatonin promotes the differentiation of secondary hair follicles in Cashmere goats, thereby improving cashmere production. MicroRNAs (miRNAs) play a crucial role in regulating post-transcriptional gene expression and influence hair follicle growth. However, the mechanism through which melatonin regulates hair follicle development via miRNA mediation remains unclear. In this study, we used RNA-seq to identify differentially expressed (DE) miRNAs during melatonin-induced growth of secondary hair follicles in inner Mongolian Cashmere goats. In total, 170 DE miRNAs were identified. Enrichment analysis revealed that the target genes of these DE miRNAs were related to biological processes such as protein modification; cytoskeletal components; and the Notch, Wnt, and MAPK signaling pathways. The miRNA-mRNA regulatory network suggested that the DE miRNA chi-let-7d-5p negatively regulates WNT2 expression. Mechanistic studies revealed that melatonin promotes the proliferation of DP cells in Cashmere goats via the chi-let-7d-5p/WNT2 axis.

Molecular Characteristics of JAK2 and Its Effect on the Milk Fat and Casein Synthesis of Ovine Mammary Epithelial Cells

In addition to its association with milk protein synthesis via the Janus kinase-signal transducer and activator of transcription (JAK-STAT) pathway, JAK2 also affects milk fat synthesis. However, to date, there have been no reports on the effect of JAK2 on ovine mammary epithelial cells (OMECs), which directly determine milk yield and milk contents. In this study, the coding sequence (CDS) region of ovine JAK2 was cloned and identified and its tissue expression and localization in ovine mammary glands, as well as its effects on the viability, proliferation, and milk fat and casein levels of OMECs, were also investigated. The CDS region of ovine JAK2, 3399 bp in length, was cloned and its authenticity was validated by analyzing its sequence similarity with JAK2 sequences from other animal species using a phylogenetic tree. JAK2 was found to be expressed in six ovine tissues, with the highest expression being in the mammary gland. Over-expressed JAK2 and three groups of JAK2 interference sequences were successfully transfected into OMECs identified by immunofluorescence staining. When compared with the negative control (NC) group, the viability of OMECs was increased by 90.1% in the pcDNA3.1-JAK2 group. The over-expression of JAK2 also increased the number and ratio of EdU-labeled positive OMECs, as well as the expression levels of three cell proliferation marker genes. These findings show that JAK2 promotes the viability and proliferation of OMECs. Meanwhile, the triglyceride content in the over-expressed JAK2 group was 2.9-fold higher than the controls and the expression levels of four milk fat synthesis marker genes were also increased. These results indicate that JAK2 promotes milk fat synthesis. Over-expressed JAK2 significantly up-regulated the expression levels of casein alpha s2 (CSN1S2), casein beta (CSN2), and casein kappa (CSN3) but down-regulated casein alpha s1 (CSN1S1) expression. In contrast, small interfered JAK2 had the opposite effect to JAK2 over-expression on the viability, proliferation, and milk fat and milk protein synthesis of OMECs. In summary, these results demonstrate that JAK2 promotes the viability, proliferation, and milk fat synthesis of OMECs in addition to regulating casein expression in these cells. This study contributes to a better comprehension of the role of JAK2 in the lactation performance of sheep.

Mitochondrial Calcium Uniporter Regulator 1 (MCUR1) Relieves Mitochondrial Damage Induced by Lipopolysaccharide by Mediating Mitochondrial Ca2+ Homeostasis in Bovine Mammary Epithelial Cells

The metabolic stress triggered by negative energy balance after calving induces mitochondrial damage of bovine mammary epithelial cells. Mitochondrial calcium uniporter regulator 1 (MCUR1) is a key protein-coding gene that mediates mitochondrial calcium ion (Ca²⁺) uptake and plays an important role in mediating homeostasis of mitochondria. The aim of the present study was to elucidate the effects of MCUR1-mediated Ca²⁺ homeostasis on mitochondria of bovine mammary epithelial cells in response to an inflammatory challenge with lipopolysaccharide (LPS). Exogenous LPS resulted in upregulation of the MCUR1 mRNA and protein abundance, mitochondrial Ca²⁺ content, and mitochondrial reactive oxygen species (Mito-ROS) content while decreasing mitochondrial membrane potential, causing mitochondrial damage, and increasing the rate of apoptosis. Ryanodine pretreatment attenuated the upregulation of the mitochondrial Ca²⁺ content and Mito-ROS content induced by LPS. Overexpression of MCUR1 increased the mitochondrial Ca²⁺ content and Mito-ROS content, while it decreased mitochondrial membrane potential, damaged mitochondria, and induced cell apoptosis. In addition, knockdown of MCUR1 by small interfering RNA attenuated LPS-induced mitochondrial dysfunction by inhibiting mitochondrial Ca²⁺ uptake. Our results revealed that exogenous LPS induces MCUR1-mediated mitochondrial Ca²⁺ overload in bovine mammary epithelial cells, which leads to mitochondrial injury. Thus, MCUR1-mediated Ca²⁺ homeostasis may be a potential therapeutic target against mitochondrial damage induced by metabolic challenges in bovine mammary epithelial cells.

Effect of parturition time and photoperiod on milk production, quality, and somatic cell count traits of pure and crossbred goats in a different production system

This is the first study to compare the effects of parturition time and photoperiod on milk yield, quality, and somatic cell count (SCC) traits in pure and crossbred doe genotypes reared with supplementary feeding at various physiological stages in a mountain-pasture grazing system. Data were collected from Hair, Alpine × Hair F₁ (AHF₁), and Saanen × Hair F₁ (SHF₁) crossbred doe genotypes with an average live weight of 49.60±0.40 kg. Hair doe had significantly lower (P < 0.001) milk yields and lactation length, but higher (P < 0.001) fat, protein, and lactose contents and electrical conductivity (EC) than AHF₁ and SHF₁ crossbred doe genotypes. Does giving parturition in the hours of darkness had higher (P < 0.05) milk volume and lower (P < 0.05) SCC compared to does giving parturition in the hours of daylight. Daily milk yield was positively correlated with daylight (P < 0.05; r = 0.50 to 0.53), while milk fat, protein, and lactose contents were negatively correlated with daylight (P < 0.05 to P < 0.001; r = -0.51 to -0.82, respectively) in Hair, AHF₁, and SHF₁ crossbred doe genotypes. Lactation stages and daily milk yield levels significantly affected (P < 0.05) the SCC, pH, total solids, and freezing point of milk. It was concluded that the milk yield traits of goats can be rapidly improved by better modelling of physiological and environmental variations such as photoperiod, reproduction, and hormonal effects in sustainable production systems.

Oar-miR-16b and oar-miR-27a: negatively correlated with milk yield and milk protein in sheep

The characterization of miRNAs from sheep milk and their effect on milk yield and composition in sheep are remain unclear. Therefore, the aim of this study was to determine the expression pattern of several important miRNAs, which are associated with lactation in the sheep milk between high and low lactating-yield ewe groups. In addition to experimentally obtained miRNA expression results, the miRNA target genes were determined by bioinformatics analysis to identify biological pathways involved. miRNAs found to differ significantly in the expression level between the groups were oar-miR-181a, oar-miR-23a, oar-miR-27a, oar-miR-16b and oar-miR-374. Also, oar-miR-27a was shown negative correlation with milk protein and lactose contents while oar-miR-16b was shown negative correlation with milk yield in the high milk yield group. The highest connected hub genes for miR-27a target genes were determined as MAPK14 and PPARG. Also, six genes (HSPA4L, DNAJA2, ATP6V1B2, PPP2R1A, PPP2R1B, and PRKAR2A) were detected as hub genes for miR-16b. In this study, the relationship between expression profiles of several important miRNAs in sheep milk and milk yield and milk composition were investigated for the first time in high and low lactating yield groups.

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.

Non-coding transcriptomic profiles in the sheep mammary gland during different lactation periods

Sheep milk production is a dynamic and multifactorial trait regulated by diverse biological mechanisms. To improve the quality and production of sheep milk, it is necessary to understand the underlying non-coding transcriptomic mechanisms. In this study, ribonucleic acid-sequencing (RNA-seq) was used to profile the expression of microRNAs (miRNAs) and circular RNAs (circRNAs) in the sheep mammary gland at three key lactation time points (perinatal period, PP; early lactation, EL; and peak lactation, PL). A total of 2,369 novel circRNAs and 272 miRNAs were profiled, of which 348, 373, and 36 differentially expressed (DE) circRNAs and 30, 34, and 7 DE miRNAs were detected in the comparison of EL vs. PP, PL vs. PP, and PL vs. EL, respectively. A series of bioinformatics analyses including functional enrichment, machine learning prediction, and competing endogenous RNA (ceRNA) network analyses were conducted to identify subsets of the potential candidate miRNAs (e.g., oar_miR_148a, oar_miR_362, and oar_miR_432) and circRNAs (e.g., novel_circ_0011066, novel_circ_0010460, and novel_circ_0006589) involved in sheep mammary gland development. Taken together, this study offers a window into the dynamics of non-coding transcriptomes that occur during sheep lactation and may provide further insights into miRNA and circRNA that influence sheep mammary gland development.

Transcriptome Adaptation of the Ovine Mammary Gland to Dietary Supplementation of Extruded Linseed

Several dietary strategies were adopted to reduce saturated fatty acids and increase beneficial fatty acids (FA) for human health. Few studies are available about the pathways/genes involved in these processes. Illumina RNA-sequencing was used to investigate changes in the ovine mammary gland transcriptome following supplemental feeding with 20% extruded linseed. Comisana ewes in mid-lactation were fed a control diet for 28 days (control period) followed by supplementation with 20% DM of linseed panel for 28 days (treatment period). Milk production was decreased by 30.46% with linseed supplementation. Moreover, a significant reduction in fat, protein and lactose secretion was also observed. Several unsaturated FAs were increased while short and medium chain saturated FAs were decreased by linseed treatment. Around four thousand (1795 up- and 2133 down-regulated) genes were significantly differentially regulated by linseed supplementation. The main pathways affected by linseed supplementation were those involved in the energy balance of the mammary gland. Principally, the mammary gland of fed linseed sheep showed a reduced abundance of transcripts related to the synthesis of lipids and carbohydrates and oxidative phosphorylation. Our study suggests that the observed decrease in milk saturated FA was correlated to down-regulation of genes in the lipid synthesis and lipid metabolism pathways.

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.

Functional characterization of Copy Number Variations regions in Djallonké sheep

A total of 184 Djallonké (West African Dwarf) sheep of Burkina Faso were analysed for Copy Number Variations (CNV) using Ovine 50 K SNP BeadChip genotyping data and two different CNV calling platforms: PennCNV and QuantiSNP. Analyses allowed to identify a total of 63 candidate Copy Number Variations Regions (CNVR) on 11 different ovine chromosomes covering about 82.5 Mb of the sheep genome. Gene-annotation enrichment analysis allowed to identify a total of 751 potential candidate ovine genes located in the candidate CNVR bounds. Functional annotation allowed to identify five statistically significant Functional Clusters (FC; enrichment factor > 1.3) involving 61 candidate genes. All genes forming significantly enriched FC were located on ovine chromosome (OAR) 21. FC1 (22 genes including PAG4 and PAG6) and FC5 (three genes: CTSC, CTSW and CTSF), coding proteases (peptidases and cathepsins, respectively), were involved in reproductive performance and modulation of gestation. Both FC3 and FC4 were involved in inflammatory and immunologic response through coding serum amyloid A and B-box-type zinc finger proteins, respectively. Finally, FC2 consisted of 27 genes (including OR10G6 and OR8B8) involved in olfactory receptor activity, key for animals adapting to new food resources. CNVR identified on at least 15% of individuals were considered CNVR hotspots and further overlapped with previously reported quantitative trait loci (QTL). CNVR hotspots spanning genes putatively involved with lipid metabolism (SKP1, TCF7, JADE2, UBE2B and SAR1B) and differential expression in mammary gland (SEC24A and CDKN2AIPNL) on OAR5 and dairy traits (CCDC198 and SLC35F4) on OAR7 overlapped with QTL associated with lipid metabolism, milk protein yield and milk fat percentage. Information obtained from local sheep populations naturally adapted to harsh environments contributes to increase our understanding of the genomic importance of CNV.

Deep small RNA-Seq reveals microRNAs expression profiles in lactating mammary gland of 2 sheep breeds with different milk performance

Micro ribonucleic acid (miRNA) is a type of noncoding RNA, and it has been revealed to play important roles in the activity of the mammary gland (MG) in some species. However, the function of miRNAs in MG of sheep is poorly understood. In the study, Gansu Alpine Merino (GAM; n = 9) and Small-tailed Han sheep (STH; n = 9) with different milk production traits were investigated. Microstructures and the expression profile of miRNAs of MG tissues at peak lactation were studied. Mature alveolar lumens of MG in appearance were larger in STH than GAM. The expression levels of CSN2 and the content of rough endoplasmic reticulum were also higher in STH ewes than GAM ewes. A total of 124 mature miRNAs were expressed, and 18 of these were differentially expressed between the 2 breeds. The KEGG analysis results showed that the targeted genes of differentially expressed miRNAs were mainly involved in some metabolic pathways and signaling pathways related to MG development, milk protein, and fat synthesis. The findings in the study can improve our understanding of the roles of miRNAs in the development and lactation of MG in sheep.

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.

Comparative mammary gland postnatal development and tumourigenesis in the sheep, cow, cat and rabbit: Exploring the menagerie

Sheep, cows, cats, and rabbits are kept by humans for agricultural purposes and as companion animals. Much of the mammary research in these species has focussed on mastitis in the case of ruminants and rabbits, and mammary tumourigenesis in cats and rabbits. However, similarities with the human breast suggest that these species may be currently underutilised as valuable comparative models of breast development and disease. The mammary gland undergoes cyclical postnatal development that will be considered here in the context of these non-traditional model species, with a focus on the mammary microenvironment at different postnatal developmental stages. The second part of this review will consider mammary tumour development. Ruminants are thought to be relatively 'resistant' to mammary tumourigenesis, likely due to multiple factors including functional properties of ruminant mammary stem/progenitor cells, diet, and/or the fact that production animals undergo a first parity soon after puberty. By contrast, unneutered female cats and rabbits have a propensity to develop mammary neoplasms, and subsets of these may constitute valuable comparative models of breast cancer.

Small RNA deep sequencing reveals the expressions of microRNAs in ovine mammary gland development at peak-lactation and during the non-lactating period

MicroRNAs (miRNAs) are small non-coding RNAs that are involved in mammary gland development and lactation in livestock. Little is known about the roles of miRNAs in ovine mammary gland development, hence in this study the expression profiles of miRNAs of the mammary gland tissues of ewes at peak-lactation and during the non-lactating period were investigated using RNA sequencing. A total of 147 mature miRNAs were expressed in the two periods. Compared with peak-lactation, eight miRNAs in the non-lactating ewe mammary gland were significantly up-regulated, whereas fifteen miRNAs were down-regulated. A KEGG analysis revealed that the target genes of the up-regulated miRNAs were significantly enriched in lysosome, Wnt and MAPK signaling pathways, while the target genes of down-regulated miRNAs were significantly enriched in the PI3K-Akt signaling pathway, protein processing in endoplasmic reticulum and axon guidance. These results suggest that further study of the differentially expressed miRNAs could provide a better understanding of the molecular mechanisms of mammary development and lactation in sheep.

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

Simple Summary

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

Abstract

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

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.

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.

Effect of Diet and Type of Pregnancy on Transcriptional Expression of Selected Genes in Sheep Mammary Gland

Simple Summary

An experiment was designed to determine the effect of diet and type of pregnancy on the mammary gland development, measured by the transcriptional expression of genes involved in angiogenesis and cell turnover/lactogenesis. To that end, twin and single-bearing ewes were fed naturalized pasture or red clover from day −45 pre-partum until day +60 post-partum, taking samples of mammary tissue at day −10, +30 and +60 post-partum. The results showed that the group of twin-bearing ewes fed red clover was the best combination to increase the expression of genes associated to angiogenesis and cell turnover/lactogenesis in the mammary gland.

Abstract

These trials were carried out to determine firstly the effect of diet and type of pregnancy on the transcriptional expression of genes involved in angiogenesis and cell turnover/lactogenesis inside the sheep mammary gland from late gestation to late lactation. Eighteen Ile de France sheep, 8 twin- and 10 single-bearing ewes were alloted into two groups according to their diet, either based on ad libitum naturalized pasture or red clover hay plus lupine from day −45 pre-partum until day +60 post-partum. Samples from diets and mammary glands were collected at day −10 pre partum (time 1), day +30 (time 2) and day +60 post-partum (time 3) and analyzed by qRT-PCR. Additionally, samples from longissimus dorsi muscle were taken from lambs twice, at weaning and 45 days later, to determine the effect of the maternal treatment with regard to diet and type of pregnancy, on the mRNA expression of genes involved in lipid metabolism. The data was processed using the lme4 package for R, and SPSS Statistics 23.0 for Windows^®. The results showed that the group of twin-bearing ewes fed red clover showed a higher expression of genes involved in angiogenesis before lambing and in cell turnover/lactogenesis during late lactation, explained by a lamb survival mechanism to delay apoptosis as a way to keep a secretory cells population and boosted by the diet quality, assuring a longer milk production potential during late lactation. Regarding lambs, apparently the maternal diet would influence the transcriptional expression of lipogenic enzymes in the longissimus dorsi muscle after weaning, but further studies are necessary to validate these results. In summary, Twin-bearing ewes fed red clover performed best at increasing the expression of genes associated with angiogenesis and cell turnover/lactogenesis in the mammary gland.

Development and evaluation of a milk protein transcript depletion method for differential transcriptome analysis in mammary gland tissue

Background

In the mammary gland transcriptome of lactating dairy cows genes encoding milk proteins are highly abundant, which can impair the detection of lowly expressed transcripts and can bias the outcome in global transcriptome analyses. Therefore, the aim of this study was to develop and evaluate a method to deplete extremely highly expressed transcripts in mRNA from lactating mammary gland tissue.

Results

Selective RNA depletion was performed by hybridization of antisense oligonucleotides targeting genes encoding the caseins (CSN1S1, CSN1S2, CSN2 and CSN3) and whey proteins (LALBA and PAEP) within total RNA followed by RNase H-mediated elimination of the respective transcripts. The effect of the RNA depletion procedure was monitored by RNA sequencing analysis comparing depleted and non-depleted RNA samples from Escherichia coli (E. coli) challenged and non-challenged udder tissue of lactating cows in a proof of principle experiment. Using RNase H-mediated RNA depletion, the ratio of highly abundant milk protein gene transcripts was reduced in all depleted samples by an average of more than 50% compared to the non-depleted samples. Furthermore, the sensitivity for discovering transcripts with marginal expression levels and transcripts not yet annotated was improved. Finally, the sensitivity to detect significantly differentially expressed transcripts between non-challenged and challenged udder tissue was increased without leading to an inadvertent bias in the pathogen challenge-associated biological signaling pathway patterns.

Conclusions

The implementation of selective RNase H-mediated RNA depletion of milk protein gene transcripts from the mammary gland transcriptome of lactating cows will be highly beneficial to establish comprehensive transcript catalogues of the tissue that better reflects its transcriptome complexity.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5781-3) contains supplementary material, which is available to authorized users.

Milk-derived mammary epithelial cells as non-invasive source to define stage-specific abundance of milk protein and fat synthesis transcripts in native Sahiwal cows and Murrah buffaloes

The molecular physiology of milk production of two important dairy species; Sahiwal cows (Bos indicus) and Murrah buffaloes (Bubalus bubalis) are not fully understood due to constraints in obtaining mammary tissue samples because of sacred and ethical reasons. The present study suggests the use of milk-derived mammary epithelial cells (MECs) as a non-invasive method to understand molecular aspects of lactation biology in dairy animals. A total of 76 MECs were collected from five different lactation periods viz. colostrum (0-2), early (5-20), peak (30-50), mid (90-140) and late lactation (> 215 days) stages from Sahiwal cows and Murrah buffaloes to study the transcription kinetics of milk protein, fat synthesis, and their regulatory genes. Significant changes were observed in milk composition of both dairy species with lactation stages. High mRNA abundance of all milk protein and fat synthesis genes was observed in MECs of Murrah buffaloes as compared to Sahiwal cows. The mRNA abundance of caseins (CSN1S1, CSN1S2, CSN2, and CSN3) and whey protein (LALBA, LF) were higher in early lactation stage. Similarly, the expression of milk fat synthesis genes (SCD, BTN1A1, ACACA, GPAM, FAPB3, FASN) was also high in early lactation stage. The relative abundance of 4 regulatory genes (JAK2, STAT5, SREBF1 and EIF4BP41) remained high during early lactation indicating their regulatory roles in lactogenesis process. Overall, results suggested a significant effect of lactation stages on milk composition and transcription abundance of milk protein and fat synthesis genes. The present study establishes the fact that milk-derived MECs could be utilized as a valuable source to understand mammary gland functioning of native cows and buffaloes.

Genetic Analyses Confirm SNPs in HSPA8 and ERBB2 are Associated with Milk Protein Concentration in Chinese Holstein Cattle

Heat shock 70 kDa protein 8 (HSPA8) and erb-b2 receptor tyrosine kinase 2 (ERBB2) were the promising candidates for milk protein concentration in dairy cattle revealed through previous RNA sequencing (RNA-Seq) study. The objective of this post-RNA-Seq study was to confirm genetic effects of HSPA8 and ERBB2 on milk protein concentration in a large Chinese Holstein population and to evaluate the genetic effects of both genes on other milk production traits. There were 2 single-nucleotide polymorphisms (SNPs) identified for HSPA8 and 11 SNPs for ERBB2 by sequencing 17 unrelated Chinese Holstein sires. The SNP-rs136632043 in HSPA8 had significant associations with all five milk production traits (p = 0.0086 to p < 0.0001), whereas SNP-rs132976221 was remarkably associated with three yield traits (p < 0.0001). Nine (ss1996900615, rs109017161, rs109122971, ss1996900614, rs110133654, rs109941438, rs110552983, rs133031530, and rs109763505) of 11 SNPs in ERBB2 were significantly associated with milk protein percentage (p = 0.0177 to p < 0.0001). A 12 Kb haplotype block was formed in ERBB2 and haplotype associations revealed similar effects on milk protein traits. Our findings confirmed the significant genetic effects of HSPA8 and ERBB2 on milk protein concentration and other milk production traits and SNP phenotypic variances above 1% may serve as genetic markers in dairy cattle breeding programs.

Increasing metabolisable energy and protein supplementation to stimulate the subsequent milk production during late gestation by increasing proliferation and reducing apoptosis in goat mammary gland prepartum

In total, 32 pregnant goats were assigned randomly to four diets fed from Day 100 of pregnancy to Day 30 after parturition, to determine the effects of metabolisable energy (ME) and metabolisable protein (MP) supplementation levels on feed intake, subsequent colostrum and milk production and expression of genes regulating mammary-cell proliferation and apoptosis. Diets were as follows: (1) diet with ME and MP provided according to NRC recommendations (control), (2) diet with extra 10% ME, (3) diet with extra 10% MP, and (4) diet 1 with 10% extra of both ME and MP. Mammary biopsies were obtained from each udder half 24 h after parturition. Feed intake (g/day), and colostrum (kg/day) and milk (kg/month) production increased when the extra ME and MP were provided together prepartum and in early lactation (P &lt; 0.05). Relative mRNA expressions significantly increased in the mammary gland of insulin-like growth factor 1 (IGF-1, 4.3-fold), IGF-1 receptor (IGF-1R, 3.6-fold) and B-cell lymphoma 2 (Bcl-2, 4.6-fold), whereas insulin-like growth factor binding protein 3 (IGFBP-3, 3.2-fold), Bcl-2-associated X protein (Bax, 16.7-fold) and the ratio of Bax:Bcl-2 expressions significantly decreased (69.8-fold) with increased ME and MP levels fed in late gestation. In conclusion, colostrum production and milk yield in the early lactation period are sensitive to nutrient supply during gestation, where increased dietary ME as well as MP supplementation levels during late gestation will favour mammary development, by increasing expression of genes stimulating cellular proliferation (IGF-1, IGF-1R, Bcl-2) and reduced those stimulating apoptosis (IGFBP-3, Bax).

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

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

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.

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.

Electronic supplementary material

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

Transcriptomic Impacts of Rumen Epithelium Induced by Butyrate Infusion in Dairy Cattle in Dry Period

The purpose of this study was to evaluate the effects of butyrate infusion on rumen epithelial transcriptome. Next-generation sequencing (NGS) and bioinformatics are used to accelerate our understanding of regulation in rumen epithelial transcriptome of cattle in the dry period induced by butyrate infusion at the level of the whole transcriptome. Butyrate, as an essential element of nutrients, is a histone deacetylase (HDAC) inhibitor that can alter histone acetylation and methylation, and plays a prominent role in regulating genomic activities influencing rumen nutrition utilization and function. Ruminal infusion of butyrate was following 0-hour sampling (baseline controls) and continued for 168 hours at a rate of 5.0 L/day of a 2.5 M solution as a continuous infusion. Following the 168-hour infusion, the infusion was stopped, and cows were maintained on the basal lactation ration for an additional 168 hours for sampling. Rumen epithelial samples were serially collected via biopsy through rumen fistulae at 0-, 24-, 72-, and 168-hour (D1, D3, D7) and 168-hour post-infusion (D14). In comparison with pre-infusion at 0 hours, a total of 3513 genes were identified to be impacted in the rumen epithelium by butyrate infusion at least once at different sampling time points at a stringent cutoff of false discovery rate (FDR) < 0.01. The maximal effect of butyrate was observed at day 7. Among these impacted genes, 117 genes were responsive consistently from day 1 to day 14, and another 42 genes were lasting through day 7. Temporal effects induced by butyrate infusion indicate that the transcriptomic alterations are very dynamic. Gene ontology (GO) enrichment analysis revealed that in the early stage of rumen butyrate infusion (on day 1 and day 3 of butyrate infusion), the transcriptomic effects in the rumen epithelium were involved with mitotic cell cycle process, cell cycle process, and regulation of cell cycle. Bioinformatic analysis of cellular functions, canonical pathways, and upstream regulator of impacted genes underlie the potential mechanisms of butyrate-induced gene expression regulation in rumen epithelium. The introduction of transcriptomic and bioinformatic technologies to study nutrigenomics in the farm animal presented a new prospect to study multiple levels of biological information to better apprehend the whole animal response to nutrition, physiological state, and their interactions. The nutrigenomics approach may eventually lead to more precise management of utilization of feed resources in a more effective approach.

SNP co-association and network analyses identify E2F3, KDM5A and BACH2 as key regulators of the bovine milk fatty acid profile

The fatty acid (FA) profile has a considerable impact on the nutritional and technological quality of milk and dairy products. The molecular mechanism underlying the regulation of fat metabolism in bovine mammary gland have been not completely elucidated. We conducted genome-wide association studies (GWAS) across 65 milk FAs and fat percentage in 1,152 Brown Swiss cows. In total, we identified 175 significant single nucleotide polymorphism (SNPs) spanning all chromosomes. Pathway analyses revealed that 12:0 was associated with the greatest number of overrepresented categories/pathways (e.g. mitogen-activated protein kinase (MAPK) activity and protein phosphorylation), suggesting that it might play an important biological role in controlling milk fat composition. An Associated Weight Matrix approach based on SNP co-associations predicted a network of 791 genes related to the milk FA profile, which were involved in several connected molecular pathways (e.g., MAPK, lipid metabolism and hormone signalling) and undetectable through standard GWAS. Analysis of transcription factors and their putative target genes within the network identified BACH2, E2F3 and KDM5A as key regulators of milk FA metabolism. These findings contribute to increasing knowledge of FA metabolism and mammary gland functionality in dairy cows and may be useful in developing targeted breeding practices to improve milk quality.

Quantitative proteome analysis of bovine mammary gland reveals protein dynamic changes involved in peak and late lactation stages

Mammary gland is an important organ for milk synthesis and secretion. It undergoes dramatic physiological changes to adapt the shift from peak to late lactation stage. Protein plays a final very vital role in many life functions, and the protein changes during different lactation stages potentially reflect the biology of lactation and the functions of mammary gland in cows. In current study, we adopted tandem mass tags label-based quantitative analysis technique and to investigate proteome changes occurring in bovine mammary gland from peak to late lactation stages. A total of 3753 proteins from mammary tissues taken at two lactation points from four individual cows by biopsy were quantified, out of which 179 proteins were expressed differentially between two stages. We observed five new DEPs (AACS, DHCR7, GSTM3, SFRP1 and SFRP4) and nine functional well-studies known proteins (PLIN2, LPIN1, PLIN3, GSN, CD74, MMP2, SOD1, SOD3 and GPX3) related to milk performance and mammary morphology. Bioinformatics analyses of the DEPs showed a majority of the up-regulated proteins during late lactation stage were related to apoptosis and immune process, while the downregulated proteins were mainly involved in localization, lipid metabolic and transport process. This suggests that the mammary gland can adapt to different molecular functions according to the biological need of the animal. From the integrated analysis of the differentially expressed proteins with known quantitative trait loci and genome-wide association study data, we identified 95 proteins may potentially affect milking performance. We expect findings in this study could be a valuable resource for future studies investigating the bovine proteome and functional studies.

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.

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.

Comprehensive RNA-Seq profiling to evaluate lactating sheep mammary gland transcriptome

RNA-Seq enables the generation of extensive transcriptome information providing the capability to characterize transcripts (including alternative isoforms and polymorphism), to quantify expression and to identify differential regulation in a single experiment. Our aim in this study was to take advantage of using RNA-Seq high-throughput technology to provide a comprehensive transcriptome profiling of the sheep lactating mammary gland. Eight ewes of two dairy sheep breeds with differences in milk production traits were used in this experiment (four Churra and four Assaf ewes). Milk samples from these animals were collected on days 10, 50, 120 and 150 after lambing to cover the various physiological stages of the mammary gland across the complete lactation. RNA samples were extracted from milk somatic cells. The RNA-Seq dataset was generated using an Illumina HiSeq 2000 sequencer. The information reported here will be useful to understand the biology of lactation in sheep, providing also an opportunity to characterize their different patterns on milk production aptitude.

Zinc Finger Homeodomain Factor Zfhx3 Is Essential for Mammary Lactogenic Differentiation by Maintaining Prolactin Signaling Activity

The zinc finger homeobox 3 (ZFHX3, also named ATBF1 for AT motif binding factor 1) is a transcription factor that suppresses prostatic carcinogenesis and induces neuronal differentiation. It also interacts with estrogen receptor α to inhibit cell proliferation and regulate pubertal mammary gland development in mice. In the present study, we examined whether and how Zfhx3 regulates lactogenic differentiation in mouse mammary glands. At different stages of mammary gland development, Zfhx3 protein was expressed at varying levels, with the highest level at lactation. In the HC11 mouse mammary epithelial cell line, an in vitro model of lactogenesis, knockdown of Zfhx3 attenuated prolactin-induced β-casein expression and morphological changes, indicators of lactogenic differentiation. In mouse mammary tissue, knock-out of Zfhx3 interrupted lactogenesis, resulting in underdeveloped glands with much smaller and fewer alveoli, reduced β-casein expression, accumulation of large cytoplasmic lipid droplets in luminal cells after parturition, and failure in lactation. Mechanistically, Zfhx3 maintained the expression of Prlr (prolactin receptor) and Prlr-Jak2-Stat5 signaling activity, whereas knockdown and knock-out of Zfhx3 in HC11 cells and mammary tissues, respectively, decreased Prlr expression, Stat5 phosphorylation, and the expression of Prlr-Jak2-Stat5 target genes. These findings indicate that Zfhx3 plays an essential role in proper lactogenic development in mammary glands, at least in part by maintaining Prlr expression and Prlr-Jak2-Stat5 signaling activity.