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

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.

MicroRNA-148a Targets DNMT1 and PPARGC1A to Regulate the Viability, Proliferation, and Milk Fat Synthesis of Ovine Mammary Epithelial Cells

In this study, the expression profiles of miR-148a were constructed in eight different ovine tissues, including mammary gland tissue, during six different developmental periods. The effect of miR-148a on the viability, proliferation, and milk fat synthesis of ovine mammary epithelial cells (OMECs) was investigated, and the target relationship of miR-148a with two predicted target genes was verified. The expression of miR-148a exhibited obvious tissue-specific and temporal-specific patterns. miR-148a was expressed in all eight ovine tissues investigated, with the highest expression level in mammary gland tissue (p < 0.05). Additionally, miR-148a was expressed in ovine mammary gland tissue during each of the six developmental periods studied, with its highest level at peak lactation (p < 0.05). The overexpression of miR-148a increased the viability of OMECs, the number and percentage of Edu-labeled positive OMECs, and the expression levels of two cell-proliferation marker genes. miR-148a also increased the percentage of OMECs in the S phase. In contrast, transfection with an miR-148a inhibitor produced the opposite effect compared to the miR-148a mimic. These results indicate that miR-148a promotes the viability and proliferation of OMECs in Small-tailed Han sheep. The miR-148a mimic increased the triglyceride content by 37.78% (p < 0.01) and the expression levels of three milk fat synthesis marker genes in OMECs. However, the miR-148a inhibitor reduced the triglyceride level by 87.11% (p < 0.01). These results suggest that miR-148a promotes milk fat synthesis in OMECs. The dual-luciferase reporter assay showed that miR-148a reduced the luciferase activities of DNA methyltransferase 1 (DNMT1) and peroxisome proliferator-activated receptor gamma coactivator 1-A (PPARGC1A) in wild-type vectors, suggesting that they are target genes of miR-148a. The expression of miR-148a was highly negatively correlated with PPARGC1A (r = -0.789, p < 0.001) in ovine mammary gland tissue, while it had a moderate negative correlation with DNMT1 (r = -0.515, p = 0.029). This is the first study to reveal the molecular mechanisms of miR-148a underlying the viability, proliferation, and milk fat synthesis of OMECs in sheep.

Genetic diversity of United States Rambouillet, Katahdin and Dorper sheep

BACKGROUND

Managing genetic diversity is critically important for maintaining species fitness. Excessive homozygosity caused by the loss of genetic diversity can have detrimental effects on the reproduction and production performance of a breed. Analysis of genetic diversity can facilitate the identification of signatures of selection which may contribute to the specific characteristics regarding the health, production and physical appearance of a breed or population. In this study, breeds with well-characterized traits such as fine wool production (Rambouillet, N = 745), parasite resistance (Katahdin, N = 581) and environmental hardiness (Dorper, N = 265) were evaluated for inbreeding, effective population size (Ne), runs of homozygosity (ROH) and Wright's fixation index (FST) outlier approach to identify differential signatures of selection at 36,113 autosomal single nucleotide polymorphisms (SNPs).

RESULTS

Katahdin sheep had the largest current Ne at the most recent generation estimated with both the GONe and NeEstimator software. The most highly conserved ROH Island was identified in Rambouillet with a signature of selection on chromosome 6 containing 202 SNPs called in an ROH in 50 to 94% of the individuals. This region contained the DCAF16, LCORL and NCAPG genes that have been previously reported to be under selection and have biological roles related to milk production and growth traits. The outlier regions identified through the FST comparisons of Katahdin with Rambouillet and Dorper contained genes with known roles in milk production and mastitis resistance or susceptibility, and the FST comparisons of Rambouillet with Katahdin and Dorper identified genes related to wool growth, suggesting these traits have been under natural or artificial selection pressure in these populations. Genes involved in the cytokine-cytokine receptor interaction pathways were identified in all FST breed comparisons, which indicates the presence of allelic diversity between these breeds in genomic regions controlling cytokine signaling mechanisms.

CONCLUSIONS

In this paper, we describe signatures of selection within diverse and economically important U.S. sheep breeds. The genes contained within these signatures are proposed for further study to understand their relevance to biological traits and improve understanding of breed diversity.

Transcriptome analysis to identify candidate genes related to mammary gland development of Bactrian camel (Camelus bactrianus)

Introduction

The demand for camel milk, which has unique therapeutic properties, is increasing. The mammary gland is the organ in mammals responsible for the production and quality of milk. However, few studies have investigated the genes or pathways related to mammary gland growth and development in Bactrian camels. This study aimed to compare the morphological changes in mammary gland tissue and transcriptome expression profiles between young and adult female Bactrian camels and to explore the potential candidate genes and signaling pathways related to mammary gland development.

Methods

Three 2  years-old female camels and three 5  years-old adult female camels were maintained in the same environment. The parenchyma of the mammary gland tissue was sampled from the camels using percutaneous needle biopsy. Morphological changes were observed using hematoxylin-eosin staining. High-throughput RNA sequencing was performed using the Illumina HiSeq platform to analyze changes in the transcriptome between young and adult camels. Functional enrichment, pathway enrichment, and protein-protein interaction networks were also analyzed. Gene expression was verified using quantitative real-time polymerase chain reaction (qRT-PCR).

Results

Histomorphological analysis showed that the mammary ducts and mammary epithelial cells in adult female camels were greatly developed and differentiated from those in young camels. Transcriptome analysis showed that 2,851 differentially expressed genes were obtained in the adult camel group compared to the young camel group, of which 1,420 were upregulated, 1,431 were downregulated, and 2,419 encoded proteins. Functional enrichment analysis revealed that the upregulated genes were significantly enriched for 24 pathways, including the Hedgehog signaling pathway which is closely related to mammary gland development. The downregulated genes were significantly enriched for seven pathways, among these the Wnt signaling pathway was significantly related to mammary gland development. The protein-protein interaction network sorted the nodes according to the degree of gene interaction and identified nine candidate genes: PRKAB2, PRKAG3, PLCB4, BTRC, GLI1, WIF1, DKK2, FZD3, and WNT4. The expression of fifteen genes randomly detected by qRT-PCR showed results consistent with those of the transcriptome analysis.

Discussion

Preliminary findings indicate that the Hedgehog, Wnt, oxytocin, insulin, and steroid biosynthesis signaling pathways have important effects on mammary gland development in dairy camels. Given the importance of these pathways and the interconnections of the involved genes, the genes in these pathways should be considered potential candidate genes. This study provides a theoretical basis for elucidating the molecular mechanisms associated with mammary gland development and milk production in Bactrian camels.

MicroRNA-200c Affects Milk Fat Synthesis by Targeting PANK3 in Ovine Mammary Epithelial Cells

Milk fat is the foremost nutrient of milk and a vital indicator in evaluating milk quality. Accumulating evidence suggests that microRNAs (miRNAs) are involved in the synthesis of milk fat. The miR-200c is closely related to lipid metabolism, but little is known about its effect on the synthesis of milk fat in MECs of ewes. Herein, the effect of miR-200c on the proliferation of ovine mammary epithelial cells (MECs) and its target relationship with a predicted target gene were investigated. The regulatory effects of miR-200c on the expression of the target genes and the content of triglycerides in ovine MECs were further analyzed. The results revealed that the expression level of miR-200c was differentially expressed in both eight tissues selected during lactation and in mammary gland tissues at different physiological periods. Overexpression of miR-200c inhibited the viability and proliferation of ovine MECs, while inhibition of miR-200c increased cell viability and promoted the proliferation of ovine MECs. Target gene prediction results indicated that miR-200c would bind the 3'UTR region of pantothenate kinase 3 (PANK3). Overexpression of miR-200c reduced the luciferase activity of PANK3, while inhibition of miR-200c increased its luciferase activity. These findings illustrated that miR-200c could directly interact with the target site of the PANK3. It was further found that overexpression of miR-200c reduced the expression levels of PANK3 and, thus, accelerated the synthesis of triglycerides. In contrary, the inhibitor of miR-200c promoted the expression of PANK3 that, thus, inhibited the synthesis of triglycerides in ovine MECs. Together, these findings revealed that miR-200c promotes the triglycerides synthesis in ovine MECs via increasing the lipid synthesis related genes expression by targeting PANK3.

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.

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

Background

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

Methods

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

Results

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

Conclusions

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

Graphical abstract

Supplementary Information

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

Novel Genes Associated with Dairy Traits in Sarda Sheep

The aim of the present research was to analyze the variability of 45 SNPs from different genes involved in metabolism and innate immunity to perform an association analysis with the milk yield, composition and milk coagulation traits. A population of 1112 Sarda breed sheep was sampled. Genotyping was generated by a TaqMan Open Array^TM. Thirty out of the 45 SNPs were polymorphic, and 12 displayed a minor allele frequency higher than 0.05. An association analysis showed that the variability at genes PRKAG3 and CD14 was significantly associated with the daily milk yield. The variability at PRKAG3 was also associated with the protein and casein content, somatic cell score and bacterial score. The variation at the PRKAA2 gene was associated with the milk lactose concentration. The SNPs at CD14 were also associated with the traditional milk coagulation properties, while the SNPs at GHR and GHRHR were associated with k_SR, a derived coagulation parameter related to the rate of syneresis. The information provided here is new and increases our knowledge of genotype-phenotype interactions in sheep. Our findings might be useful in appropriate breeding schemes to be set up for the Sarda sheep breed, but these should be confirmed by further studies, possibly performed on independent populations.

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.

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.

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.