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

Stress-induced epigenetic effects driven by maternal lactation in dairy cattle: a comethylation network approach

Epigenetic marks do not follow the Mendelian laws of inheritance. The environment can alter the epigenotype of an individual when exposed to different external stressors. In lactating cows, the first stages of gestation overlap with the lactation peak, creating a negative energy balance that is difficult to overcome with diet. This negative energy balance could affect early embryo development that must compete with the mammary tissue for nutrients. We hypothesize that the methylation profiles of calves born to nonlactating heifers are different from those of calves born to lactating cows. We found 50,277 differentially methylated cytosines and 2,281 differentially methylated regions between these two groups of animals. A comethylation network was constructed to study the correlation between the phenotypes of the mothers and the epigenome of the calves, revealing 265 regions associated with the phenotypes. Our study revealed the presence of DMCs and DMRs in calves gestated by heifers and lactating cows, which were linked to the dam's lactation and the calves' ICAP and milk EBV. Gene-specific analysis highlighted associations with vasculature and organ morphogenesis and cell communication and signalling. These finding support the hypothesis that calves gestated by nonlactating mothers have a different methylation profile than those gestated by lactating cows.

Integration of selective sweeps across the sheep genome: understanding the relationship between production and adaptation traits

BACKGROUND

Livestock populations are under constant selective pressure for higher productivity levels for different selective purposes. This pressure results in the selection of animals with unique adaptive and production traits. The study of genomic regions associated with these unique characteristics has the potential to improve biological knowledge regarding the adaptive process and how it is connected to production levels and resilience, which is the ability of an animal to adapt to stress or an imbalance in homeostasis. Sheep is a species that has been subjected to several natural and artificial selective pressures during its history, resulting in a highly specialized species for production and adaptation to challenging environments. Here, the data from multiple studies that aim at mapping selective sweeps across the sheep genome associated with production and adaptation traits were integrated to identify confirmed selective sweeps (CSS).

RESULTS

In total, 37 studies were used to identify 518 CSS across the sheep genome, which were classified as production (147 prodCSS) and adaptation (219 adapCSS) CSS based on the frequency of each type of associated study. The genes within the CSS were associated with relevant biological processes for adaptation and production. For example, for adapCSS, the associated genes were related to the control of seasonality, circadian rhythm, and thermoregulation. On the other hand, genes associated with prodCSS were related to the control of feeding behaviour, reproduction, and cellular differentiation. In addition, genes harbouring both prodCSS and adapCSS showed an interesting association with lipid metabolism, suggesting a potential role of this process in the regulation of pleiotropic effects between these classes of traits.

CONCLUSIONS

The findings of this study contribute to a deeper understanding of the genetic link between productivity and adaptability in sheep breeds. This information may provide insights into the genetic mechanisms that underlie undesirable genetic correlations between these two groups of traits and pave the way for a better understanding of resilience as a positive ability to respond to environmental stressors, where the negative effects on production level are minimized.

Transcriptome identification of differential mammary genes of Kazakh horses during early pregnancy

Kazakh mares have attracted widespread attention with their outstanding lactation traits. Lactation is a complex dynamic process regulated by multiple factors. The extensive application of transcriptome sequencing technology enables researchers to further explore this biological issue. This study selected three pregnant and three non-pregnant Kazakh mares as the research subject. Their mammary glands were taken for transcriptome sequencing. The results show that there are 9 lncRNAs and 122 mRNAs differentially expressed between the two groups. GO enrichment analysis shows that there are 175 molecular functions, 59 cellular components, and 555 biological processes, including cellular hormone metabolic process, hormone catabolic process, and I-kappaB kinase/NF-kappaB signaling. KEGG enrichment analysis exhibits that these differential genes are mainly enriched in the NF-kappa B signaling pathway, steroid hormone biosynthesis, breast cancer, ECM-receptor interaction, and MAPK signaling pathway. WNT4, DPP4, and NFKBIA are key nodes regulating breast activation. Conclusions: Through the comparative analysis of the transcriptome data of mammary tissues of pregnant and non-pregnant mares, relevant differentially expressed genes are screened and analyzed. This study provides valuable fundamental data for investigating candidate genes related to the lactation regulation and mammogenesis of Kazakh horses.

Transcriptome analysis reveals immune function-related mRNA expression in donkey mammary glands during four developmental stages

The low susceptibility to mastitis of female donkey (jenny) mammary glands and the strong immune properties of donkey milk are acknowledged, but little is known about the genes involved in mammary gland immunity in jennies. Herein, we used RNA-sequencing and bioinformatics analyses to explore jenny mammary gland transcriptomes and detect potential functional differentially expressed (DE) mRNAs related to immunity during four specific developmental stages: foetal (F), pubertal (P), adult parous nonlactation (N) and lactation (L). A total of 2497, 583 and 1820 DE mRNAs were identified in jenny mammary glands at F vs. P, P vs. N, and N vs. L, respectively. Gene Ontology (GO) and Kyoto Encyclopedia of Gene and Genomes (KEGG) analyses revealed numerous GO terms related to immune function, especially between F and P. Seven significantly enriched profiles were identified, among which 497 and 1261 DE mRNAs were upregulated in profiles 19 and 17. Eleven mRNAs were enriched in over 10 KEGG pathways. β-2-microglobulin (B2M), immunoglobulin heavy constant mu (IGHM), toll like receptor 2 (TLR2), toll like receptor 4 (TLR4) and myeloid differentiation factor 88 (MYD88) were mainly involved in phosphoinositide 3-kinase (PI3K)-Akt signalling, phagosome and nuclear factor kappa-B (NF-kappa B) signalling pathways. The findings provide insight into the molecular features underpinning the low prevalence of intramammary infections (i.e., mastitis) in donkeys.

Genomic selection pressure discovery using site-frequency spectrum and reduced local variability statistics in Pakistani Dera-Din-Panah goat

BACKGROUND

Population geneticists have long sought to comprehend various selection traces accumulated in the goat genome due to natural or human driven artificial selection through breeding practices, which led the wild animals to domestication, so understanding evolutionary process may helpful to utilize the full genetic potential of goat genome.

METHODS AND RESULTS

As a step forward to pinpoint the selection signals in Pakistani Dera-Din-Panah (DDP) goat, whole-genome pooled sequencing (n = 12) was performed, and 618,236,192 clean paired-end reads were mapped against ARS1 reference goat assembly. Five different selection signature statistics were applied using four site-frequency spectrum (SFS) methods (Tajima's D ([Formula: see text]), Fay and Wu's H ([Formula: see text]), Zeng's E ([Formula: see text]), [Formula: see text]) and one reduced local variability approach named pooled heterozygosity ([Formula: see text]). The under-selection regions were annotated with significant threshold values of [Formula: see text]≥4.7, [Formula: see text]≥6, [Formula: see text]≥2.5, Pool-HMM ≥ 12, and [Formula: see text]≥5 that resulted in accumulative 364 candidate gene hits. The highest genomic selection signals were observed on Chr. 4, 6, 10, 12, 15, 16, 18, 20, and 27 and harbor ADAMTS6, CWC27, RELN, MYCBP2, FGF14, STIM1, CFAP74, GNB1, CALML6, TMEM52, FAM149A, NADK, MMP23B, OPN3, FH, MFHAS1, KLKB1, RRM1, KMO, SPEF2, F11, KIT, KMO, ERI1, ATP8B4, and RHOG genes. Next, the validation of our captured genomic hits was also performed by more than one applied statistics which harbor meat production, immunity, and reproduction associated genes to strengthen our hypothesis of under-selection traits in this Pakistani goat breed. Furthermore, common candidate genes captured by more than one statistical method were subjected to gene ontology and KEGG pathway analysis to get insights of particular biological processes associated with this goat breed.

CONCLUSION

Current perception of genomic architecture of DDP goat provides a better understanding to improve its genetic potential and other economically important traits of medium to large body size, milk, and fiber production by updating the genomic insight driven breeding strategies to boost the livestock and agriculture-based economy of the country.

Investigation of lncRNA in Bos taurus Mammary Tissue during Dry and Lactation Periods

This study aims to collect RNA-Seq data from Bos taurus samples representing dry and lactating mammary tissue, identify lncRNA transcripts, and analyze findings for their features and functional annotation. This allows for connections to be drawn between lncRNA and the lactation process. RNA-Seq data from 103 samples of Bos taurus mammary tissue were gathered from publicly available databases (60 dry, 43 lactating). The samples were filtered to reveal 214 dry mammary lncRNA transcripts and 517 lactating mammary lncRNA transcripts. The lncRNAs met common lncRNA characteristics such as shorter length, fewer exons, lower expression levels, and less sequence conservation when compared to the genome. Interestingly, several lncRNAs showed sequence similarity to genes associated with strong hair keratin intermediate filaments. Human breast cancer research has associated strong hair keratin filaments with mammary tissue cellular resilience. The lncRNAs were also associated with several genes/proteins that linked to pregnancy using expression correlation and gene ontology. Such findings indicate that there are crucial relationships between the lncRNAs found in mammary tissue and the development of the tissue, to meet both the animal's needs and our own production needs; these relationships should be further investigated to ensure that we continue to breed the most resilient, efficient dairy cattle.

The functional and regulatory entities underlying free and peptide-bound amino acid transporters in the bovine mammary gland

Free and peptide-bound AA act as building blocks and key regulators of milk protein. To improve milk protein production, mammary epithelial cells of lactating mammals require extensive AA movement across the plasma membrane via multiple transport systems. Recent studies on bovine mammary cells/tissues have expanded the number of AA transporter systems identified and the knowledge on their contribution to AA utilization for milk protein synthesis and the regulatory machinery. However, in lactating cows, the exact intracellular location of mammary AA transporters and the extent of mammary net AA utilization for milk protein production remain unclear. This review highlights the existing knowledge on various characteristics, such as substrate specificity, kinetics, their effects on AA uptake and utilization, and regulatory mechanism, of recently examined bovine mammary free and peptide-bound AA transporters.

Haplotype-based association study of production and reproduction traits in multigenerational Vrindavani population

Haplotype-based association analysis promises to reveal important information regarding the effect of genetic variants on economic traits of interest. The present study aimed to evaluate the haplotype structure of Vrindavani cattle and explore the association of haplotypes with (re)production traits of economic interest. Genotyping array data of medium density (Bovine50KSNP BeadChip) on 96 randomly selected Vrindavani cows was used in the present study. Genotypes were called in GenomeStudio program while quality control was undertaken in PLINK using standard thresholds. The phenotypic traits used in the present study included age at first calving, dry days, lactation length, peak yield, total lactation milk yield, inter-calving period and service period. The haplotype structure of Vrindavani population was assessed, using a sliding window of 20 SNP with a shift of 5 SNPs at a time, in terms of the size of haplotype blocks regarding their length (in Kb) and frequency in chromosome-wise fashion. Haplotype blocks were assessed for possible association with important production and reproduction traits across three lactation cycles in Vrindavani cattle population. The first ten principal components were included in the model for haplotype-based association analysis to correct for stratification effects of assessed individuals. Multiple haplotypes were found to be associated with age at first calving, total lactation milk yield, peak yield, dry days, inter-calving period and service period. Various candidate genes were found to overlap haplotypes that were significantly associated with age at first calving (CDH18, MARCHF11, MYO10, FBXL7), total lactation milk yield (TGF, PDE1A, and COL8A1), peak yield (PPARGC1A, RCAN1, KCNE1, SMIM34 and MRPS6), dry days (CPNE4, ACAD11 and MRAS), inter-calving period (ABCG5, ABCG8 and COX7A2L) and service period (FOXL2 and PIK3CB). The putative candidate genes overlapping the significantly associated haplotypes revealed important pathways affecting the production and reproduction performance of animals. The identified genes and pathways may serve as good candidate markers to select animals for improved production and reproduction performance in future generations.

Identification of Key Genes and Potential Pathways Associated with Mastitis Induced by E. coli

Mastitis is one of the most important diseases of dairy cattle. It is an infectious disease leading to an inflammatory reaction in the cow's mammary gland. Escherichia coli is one of the common bacteria which induce mastitis in cows. The aim of this study was to identify key genes and potential pathways associated with mastitis induced by E. coli in dairy cattle using bioinformatics analysis. The gene expression profile of ten samples including five adjacent tissues from a quarter infected with Escherichia coli and five tissues from a healthy quarter of dairy cattle was assessed using GEO2R. Gene ontology and pathway analysis were performed using bioinformatics tools. A total of 156 differentially expressed genes were detected which 95 genes were upregulated and 61 genes were downregulated in adjacent tissue of quarter infected compared with healthy tissue. Cellular oxidant detoxification and oxidation-reduction process were the most significant biological process terms in gene ontology analysis. The most important pathways of DEGs were the biosynthesis of amino acids, p53 signaling pathway, and Metabolic pathways. Three important modules were identified and their path enrichment analysis was performed. There are 10 core genes, among which SOD2, COL1A2, COL3A1, POSTN, ALDH18A1, and CBS may be the main genes associated with mastitis, which can be considered as candidate genes in the prevention and carly diagnosis program of mastitis.

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.

Identification and functional analysis of m6A in the mammary gland tissues of dairy goats at the early and peak lactation stages

N6-methyladenosine (m⁶A) is the most common reversible epigenetic RNA modification in the mRNA of all higher eukaryotic organisms and plays an important role in the regulation of gene expression and cell function. In this study, m⁶A-modified methylated RNA immunoprecipitation sequencing (MeRIP-seq) and transcriptome sequencing (RNA-seq) were used to identify the key genes with m⁶A modification during mammary gland development and lactation in dairy goats. The results showed that m⁶A methylation occurred at 3,927 loci, which were significantly enriched in the 3′ untranslated region (3′UTR) and the termination codon region. In the early stage and peak stage of lactation, m⁶A methylation occurred extensively in mammary tissues, and a total of 725 differentially expressed m⁶A-modified genes were obtained, all negatively correlated with mRNA expression. In addition, Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed that different methylated genes were mainly involved in the growth and apoptosis of mammary epithelial cells through signaling pathways, such as the mitogen-activated protein kinase (MAPK) and phospholipase D pathways, and then affected the development and lactation of mammary gland. All in all, we identified and analyzed the methylation events related to the development and lactation regulation of mammary gland at the early and peak lactation stages, and provided a theoretical basis to reveal the physiological regulatory system of mammary gland development and lactation in dairy goats.

Maternal exposure to di-2-ethylhexyl phthalate (DEHP) depresses lactation capacity in mice

Increasing evidence shows that di-2-ethylhexyl phthalate (DEHP), mostly commonly used phthalate for the production of flexible polyvinyl chloride (PVC), has the potential to induce serious health risks in humans. However, the understanding of DEHP-induced lactation performance remains largely unknown. We sought to investigate the adverse effects of DEHP on lactation and examine the underlying mechanism linking DEHP exposure with the lactation alterations. We successfully adapted a maternal DEHP exposure model in female pregnant/lactating mice. Then we determined effects of DEHP exposure on food intake, body weight and milk production as well as the alterations in endocrine factors in lactating mice. The integrated metabonomic and transcriptomic analyses of the mammary gland were performed to measure the changed metabolites and genes related to DEHP exposure-induced lactation alterations. We observed the reduced food intake with elevated blood leptin and the decreased milk yield as well as the reduced levels of serum prolactin, growth hormone, insulin-like growth factor 1 and insulin after exposed to DEHP. Furthermore, 208 metabolites and 3452 genes were separately identified as differentially expressed features associated with DEHP exposure. Integrated metabonomic and transcriptomic analyses demonstrated that DEHP caused lactation depression mainly through impairing energy generation, inducing stress responses along with the hypoactivation of inflammation, reducing the production of antioxidants, disrupting hormone homeostasis and repressing the synthesis of milk constituents (the lower glucose availability for lactose synthesis; the disruption of milk fat globule membrane for lipid droplet formation; the ribosomal dysfunction and disruption of post-modifications for milk protein synthesis). We demonstrated that DEHP disrupted several lactation-related hormone homeostasis and multiple processes like energy insufficiency, inflammation activation, oxidative stress aggravation and disturbance of milk production in the mammary gland of female lactating mice. Our results provide valuable information for the health risk of plastic additive (DEHP) on female lactation dysfunction.

Mammary epithelial cell transcriptome reveals potential roles of lncRNAs in regulating milk synthesis pathways in Jersey and Kashmiri cattle

Background

Long noncoding RNAs (lncRNAs) are now proven as essential regulatory elements, playing diverse roles in many biological processes including mammary gland development. However, little is known about their roles in the bovine lactation process.

Results

To identify and characterize the roles of lncRNAs in bovine lactation, high throughput RNA sequencing data from Jersey (high milk yield producer), and Kashmiri cattle (low milk yield producer) were utilized. Transcriptome data from three Kashmiri and three Jersey cattle throughout their lactation stages were utilized for differential expression analysis. At each stage (early, mid and late) three samples were taken from each breed. A total of 45 differentially expressed lncRNAs were identified between the three stages of lactation. The differentially expressed lncRNAs were found co-expressed with genes involved in the milk synthesis processes such as GPAM, LPL, and ABCG2 indicating their potential regulatory effects on milk quality genes. KEGG pathways analysis of potential cis and trans target genes of differentially expressed lncRNAs indicated that 27 and 48 pathways were significantly enriched between the three stages of lactation in Kashmiri and Jersey respectively, including mTOR signaling, PI3K-Akt signaling, and RAP1 signaling pathways. These pathways are known to play key roles in lactation biology and mammary gland development.

Conclusions

Expression profiles of lncRNAs across different lactation stages in Jersey and Kashmiri cattle provide a valuable resource for the study of the regulatory mechanisms involved in the lactation process as well as facilitate understanding of the role of lncRNAs in bovine lactation biology.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08406-x.

Screening and Conjoint Analysis of Key lncRNAs for Milk Fat Metabolism in Dairy Cows

Long noncoding RNAs (lncRNAs) play an important regulatory role in various biological processes as a key regulatory factor. However, the complete expression profile of lncRNAs in dairy cows and its function in milk fat synthesis are unknown. In this study, RNA sequencing (RNA-seq) was used to research the whole genome expression of lncRNAs and mRNA transcripts in high and low milk fat percentage (MFP) bovine mammary epithelial cells (BMECs), and joint analysis was carried out. We identified a total of 47 differentially expressed genes (DEGs) and 38 differentially expressed lncRNAs (DELs, Padj <0.05), enrichment analysis screened out 11 candidate DEGs that may regulate milk fat metabolism. Downregulated differential gene ENPP2 (The expression level in BMECs of high milk fat dairy cows was lower than that of low milk fat cows) and upregulated differential gene BCAT1 are more likely to participate in the milk fat metabolism, and its function needs further experiments verification. The enrichment analysis of target genes predicted by DELs identified 7 cis (co-localization) and 10 trans (co-expression) candidate target genes related to milk lipid metabolism, corresponding to a total of 18 DELs. Among them, the targeting relationship between long intervening/intergenic noncoding RNA (lincRNA) TCONS_00082721 and FABP4 is worthy of attention. One hundred and fifty-six competing endogenous RNAs (ceRNAs) interaction regulation networks related to milk fat metabolism were constructed based on the expression information of DELs, differential microRNAs (miRNAs), and lipid metabolism-related target genes. The regulatory network centered on miR-145 will be the focus of subsequent experimental research. The ceRNAs regulatory network related to TCONS_00082721 and TCONS_00172817 are more likely to be involved in milk fat synthesis. These results will provide new ways to understand the complex biology of dairy cow milk fat synthesis and provide valuable information for breed improvement of Chinese Holstein cow.

Regulation of Milk Protein Synthesis by Free and Peptide-Bound Amino Acids in Dairy Cows

Milk protein (MP) synthesis in the mammary gland of dairy cows is a complex biological process. As the substrates for protein synthesis, amino acids (AAs) are the most important nutrients for milk synthesis. Free AAs (FAAs) are the main precursors of MP synthesis, and their supplies are supplemented by peptide-bound AAs (PBAAs) in the blood. Utilization of AAs in the mammary gland of dairy cows has attracted the great interest of researchers because of the goal of increasing MP yield. Supplying sufficient and balanced AAs is critical to improve MP concentration and yield in dairy cows. Great progress has been made in understanding limiting AAs and their requirements for MP synthesis in dairy cows. This review focuses on the effects of FAA and PBAA supply on MP synthesis and their underlying mechanisms. Advances in our knowledge in the field can help us to develop more accurate models to predict dietary protein requirements for dairy cows MP synthesis, which will ultimately improve the nitrogen utilization efficiency and lactation performance of dairy cows.

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.

In silico identification of variations in microRNAs with a potential impact on dairy traits using whole ruminant genome SNP datasets

MicroRNAs are small noncoding RNAs that have important roles in the lactation process and milk biosynthesis. Some polymorphisms have been studied in various livestock species from the perspective of pathology or production traits. To target variants that could be the causal variants of dairy traits, genetic variants of microRNAs expressed in the mammary gland or present in milk and localized in dairy quantitative trait loci (QTLs) were investigated in bovine, caprine, and ovine species. In this study, a total of 59,124 (out of 28 millions), 13,427 (out of 87 millions), and 4761 (out of 38 millions) genetic variants in microRNAs expressed in the mammary gland or present in milk were identified in bovine, caprine, and ovine species, respectively. A total of 4679 of these detected bovine genetic variants are located in dairy QTLs. In caprine species, 127 genetic variants are localized in dairy QTLs. In ovine species, no genetic variant was identified in dairy QTLs. This study leads to the detection of microRNA genetic variants of interest in the context of dairy production, taking advantage of whole genome data to identify microRNA genetic variants expressed in the mammary gland and localized in dairy QTLs.

TMT based deep proteome analysis of buffalo mammary epithelial cells and identification of novel protein signatures during lactogenic differentiation

The lactating mammary gland harbours numerous matured alveoli with their lumen surrounded by differentiated mammary epithelial cells (MECs), which are exclusively involved in milk synthesis and secretion. Buffalo (Bubalus bubalis) is the second major milk-producing animal, and its physiology is different from cattle. The complete protein machinery involved in MECs differentiation is still not defined in ruminants, in particular, buffalo. Therefore, we have studied the differential expression of regulated proteins in the in vitro grown buffalo MECs (BuMECs) at different time points (on 3, 6, 12, and 15 days) of their differentiation in the presence of lactogenic hormones. TMT-based MS analysis identified 4,934 proteins; of them, 681 were differentially expressed proteins (DEPs). The principal component analysis suggested a highly heterogeneous expression of DEPs at the four-time points of hormone treatment, with most of them (307) attained the highest expression on 12 days. Bioinformatics analysis revealed the association of DEPs with 24 KEGG pathways. We observed few new proteins, namely ABCA13, IVL, VPS37, CZIB, RFX7, Rab5, TTLL12, SMEK1, GDI2, and TMEM131 in BuMECs. The function of one of the highly upregulated proteins, namely involucrin in the differentiation of BuMECs was confirmed based on biochemical inhibition assay. The results further conclude that the proteins with higher abundance can be considered as the potential biomarkers for differentiation, and they may have a significant association with the lactation process in buffalo too. The proteome dataset obtained can be used to understand the species-specific variations among other lactating animals.

Short Time-Series Expression Transcriptome Data Reveal the Gene Expression Patterns of Dairy Cow Mammary Gland as Milk Yield Decreased Process

The existing research on dairy cow mammary gland genes is extensive, but there have been few reports about dynamic changes in dairy cow mammary gland genes as milk yield decrease. For the first time, transcriptome analysis based on short time-series expression miner (STEM) and histological observations were performed using the Holstein dairy cow mammary gland to explore gene expression patterns in this process of decrease (at peak, mid-, and late lactation). Histological observations suggested that the number of mammary acinous cells at peak/mid-lactation was significantly higher than that at mid-/late lactation, and the lipid droplets area secreted by dairy cows was almost unaltered across the three stages of lactation (p > 0.05). Totals of 882 and 1439 genes were differentially expressed at mid- and late lactation, respectively, compared to peak lactation. Function analysis showed that differentially expressed genes (DEGs) were mainly related to apoptosis and energy metabolism (fold change ≥ 2 or fold change ≤ 0.5, p-value ≤ 0.05). Transcriptome analysis based on STEM identified 16 profiles of differential gene expression patterns, including 5 significant profiles (false discovery rate, FDR ≤ 0.05). Function analysis revealed DEGs involved in milk fat synthesis were downregulated in Profile 0 and DEGs in Profile 12 associated with protein synthesis. These findings provide a foundation for future studies on the molecular mechanisms underlying mammary gland development in dairy cows.

Peptide profiling in cow urine reveals molecular signature of physiology-driven pathways and in-silico predicted bioactive properties

Peptidomics allows the identification of peptides that are derived from proteins. Urinary peptidomics has revolutionized the field of diagnostics as the samples represent complete systemic changes happening in the body. Moreover, it can be collected in a non-invasive manner. We profiled the peptides in urine collected from different physiological states (heifer, pregnancy, and lactation) of Sahiwal cows. Endogenous peptides were extracted from 30 individual cows belonging to three groups, each group comprising of ten animals (biological replicates n = 10). Nano Liquid chromatography Mass spectrometry (nLC-MS/MS) experiments revealed 5239, 4774, and 5466 peptides in the heifer, pregnant and lactating animals respectively. Urinary peptides of <10 kDa size were considered for the study. Peptides were extracted by 10 kDa MWCO filter. Sequences were identified by scanning the MS spectra ranging from 200 to 2200 m/z. The peptides exhibited diversity in sequences across different physiological states and in-silico experiments were conducted to classify the bioactive peptides into anti-microbial, anti-inflammatory, anti-hypertensive, and anti-cancerous groups. We have validated the antimicrobial effect of urinary peptides on Staphylococcus aureus and Escherichia coli under an in-vitro experimental set up. The origin of these peptides was traced back to certain proteases viz. MMPs, KLKs, CASPs, ADAMs etc. which were found responsible for the physiology-specific peptide signature of urine. Proteins involved in extracellular matrix structural constituent (GO:0005201) were found significant during pregnancy and lactation in which tissue remodeling is extensive. Collagen trimers were prominent molecules under cellular component category during lactation. Homophilic cell adhesion was found to be an important biological process involved in embryo attachment during pregnancy. The in-silico study also highlighted the enrichment of progenitor proteins on specific chromosomes and their relative expression in context to specific physiology. The urinary peptides, precursor proteins, and proteases identified in the study offers a base line information in healthy cows which can be utilized in biomarker discovery research for several pathophysiological studies.

Genome-wide association study and functional analyses for clinical and subclinical ketosis in Holstein cattle

Ketosis is one of the most frequent metabolic diseases in high-yielding dairy cows and is characterized by high concentrations of ketone bodies in blood, urine, and milk, causing high economic losses. The search for polymorphic genes, whose alleles have different effects on resistance to developing the disease, is of extreme importance to help select less susceptible animals. The aims of this study were to identify genomic regions associated with clinical and subclinical ketosis (β-hydroxybutyrate concentration) in North American Holstein dairy cattle and to investigate these regions to identify candidate genes and metabolic pathways associated with these traits. To achieve this, a GWAS was performed for 4 traits: clinical ketosis lactation 1, clinical ketosis lactation 2 to 5, subclinical ketosis lactation 1, and subclinical ketosis lactation 2 to 5. The estimated breeding values from 77,277 cows and 7,704 bulls were deregressed and used as pseudophenotypes in the GWAS. The top-20 genomic regions explaining the largest proportion of the genetic variance were investigated for putative genes associated with the traits through functional analyses. Regions of interest were identified on chromosomes 2, 5, and 6 for clinical ketosis lactation 1; 3, 6, and 7 for clinical ketosis lactation 2 to 5; 1, 2, and 12 for subclinical ketosis lactation 1; and 20, 11, and 25 for subclinical ketosis lactation 2 to 5. The highlighted genes potentially related to clinical and subclinical ketosis included ACAT2 and IGF1. Enrichment analysis of the list of candidate genes for clinical and subclinical ketosis showed molecular functions and biological processes involved in fatty acid metabolism, lipid metabolism, and inflammatory response in dairy cattle. Several genomic regions and SNPs related to susceptibility to ketosis in dairy cattle that were previously described in other studies were confirmed. The novel genomic regions identified in this study aid to characterize the most important genes and pathways that explain the susceptibility to clinical and subclinical ketosis in dairy cattle.

Amino Acid Metabolomic Profiles in Bovine Mammary Epithelial Cells under Essential Amino Acid Restriction

Simple Summary

Cells of the mammary gland obtain their necessary nutrients from the blood to produce milk components, such as casein. To achieve higher productivity, cows are excessively supplemented, thus generating a higher cost of production and affecting the environment. Therefore, this triggers the need for a reduction in the supplementation of essential amino acids without affecting the milk composition. The present in vitro study shows that, through homeostatic and homeorhetic processes, cells have the ability to maintain stable casein levels despite decreasing the percentage of essential amino acids (EAAs) supplied. These findings could contribute to the proposal of more efficient nutritional strategies at lower environmental and economic costs.

Abstract

Mammary epithelial cells (MECs) in culture are a useful model for elucidating mammary gland metabolism and changes that occur under different nutrient disponibility. MECs were exposed to different treatments: 100% EAA for 8 h and 24 h restriction (R); 2% EAA for 8 h and 24 h R; 2% EAA for 8 h and 24 h + 100% EAA for 8 h and 24 h restriction + re-feeding (R + RF). Western blotting and protein quantification was performed. The Kyoto Encyclopedia of Genes and Genomes (KEGG) software identified the amino acids (AAs) and signaling pathways. The chi-squared test, multiple classification analysis, and analysis of variance were used for the purification and identification of data. Intracellular casein levels were not affected. The KEGG analysis revealed that the important pathways of metabolism of AAs, which were involved in processes related to metabolism and biosynthesis of phenylalanine, tyrosine, and tryptophan (fumarate, acetyl-CoA, and tricarboxylic acid (TCA) cycle), were affected by both R and R + RF treatments, mainly through the glutamic-oxaloacetic transaminase-2 enzyme. Additionally, metabolic processes mediated by the mitochondrial malate dehydrogenase, S-adenosylmethionine synthetase, and asparagine synthase proteins positively regulated the carbohydrate pathway, pyruvate, and TCA cycles, as well as the metabolism of alanine, aspartate, and glutamate metabolism (carbohydrate and TCA cycle). We hypothesized that MECs have the capacity to utilize alternative pathways that ensure the availability of substrates for composing milk proteins.

Lactation Associated Genes Revealed in Holstein Dairy Cows by Weighted Gene Co-Expression Network Analysis (WGCNA)

Simple Summary

Weighted gene coexpression network analysis (WGCNA) is a novel approach that can quickly analyze the relationships between genes and traits. In the past few years, studies on the gene expression changes of dairy cow mammary glands were only based on transcriptome comparisons between two lactation stages. Few studies focused on the relationships between gene expression of the dairy mammary gland and lactation stage or milk composition in a lactation cycle. In this study, we detected milk yield and composition in a lactation cycle. For the first time, we constructed a gene coexpression network using WGCNA on the basis of 18 gene expression profiles during six stages of a lactation cycle by transcriptome sequencing, generating 10 specific modules. Genes in each module were performed with gene ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. Module–trait relationship analysis showed a series of potential candidates related to milk yield and composition. The current study provides an important theoretical basis for the further molecular breeding of dairy cows.

Abstract

Weighted gene coexpression network analysis (WGCNA) is a novel approach that can quickly analyze the relationships between genes and traits. In this study, the milk yield, lactose, fat, and protein of Holstein dairy cows were detected in a lactation cycle. Meanwhile, a total of 18 gene expression profiles were detected using mammary glands from six lactation stages (day 7 to calving, −7 d; day 30 post-calving, 30 d; day 90 post-calving, 90 d; day 180 post-calving, 180 d; day 270 post-calving, 270 d; day 315 post-calving, 315 d). On the basis of the 18 profiles, WGCNA identified for the first time 10 significant modules that may be related to lactation stage, milk yield, and the main milk composition content. Genes in the 10 significant modules were examined with gene ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. The results revealed that the galactose metabolism pathway was a potential candidate for milk yield and milk lactose synthesis. In −7 d, ion transportation was more frequent and cell proliferation related terms became active. In late lactation, the suppressor of cytokine signaling 3 (SOCS3) might play a role in apoptosis. The sphingolipid signaling pathway was a potential candidate for milk fat synthesis. Dairy cows at 315 d were in a period of cell proliferation. Another notable phenomenon was that nonlactating dairy cows had a more regular circadian rhythm after a cycle of lactation. The results provide an important theoretical basis for the further molecular breeding of dairy cows.

MicroRNA-193a-5p Regulates the Synthesis of Polyunsaturated Fatty Acids by Targeting Fatty Acid Desaturase 1 (FADS1) in Bovine Mammary Epithelial Cells

Cardiovascular diseases (CVDs) are seriously threatening to human life and health. Polyunsaturated fatty acids (PUFAs) are known for their role in preventing CVDs. It is beneficial to population health to promote the content of PUFAs in bovine milk. In recent years, limited research based on molecular mechanisms has focused on this field. The biological roles of numerous microRNAs (miRNAs) remain unknown. In this study, a promising and negatively correlated pair of the miRNA (miRNA-193a-5p) and a fatty acid desaturase 1 (FADS1) gene are identified and screened to explore whether they are potential factors of PUFAs’ synthesis in bovine milk. The targeted relationship between miRNA-193a-5p and FADS1 in bovine mammary epithelial cells (BMECs) is demonstrated by dual luciferase reporter assays. qRT-PCR and western blot assays indicate that both the expression of mRNA and the protein FADS1 show a negative correlation with miRNA-193a-5p expression in BMECs. Also, miR-193a-5p expression is positively correlated with the expression of genes associated with milk fatty acid metabolism, including ELOVL fatty acid elongase 6 (ELOVL6) and diacylglycerol O-acyltransferase 2 (DGAT2). The expression of fatty acid desaturase 2 (FADS2) is negatively correlated with miR-193a-5p expression in BMECs. The contents of triglycerides (TAG), eicosapentaenoic acid (EPA), docosapentaenoic acid (DPA) and docosahexaenoic acid (DHA) have a significant positive correlation with the expression of FADS1 and a significant negative correlation with the expression of miR-193a-5p in BMECs. For the first time, this study confirms that miRNA-193a-5p regulates PUFAs metabolism in BMECs by targeting FADS1, indicating that miRNA-193a-5p and FADS1 are underlying factors that improve PUFAs content in bovine milk.

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.

ASCT2 Is Involved in SARS-Mediated β-Casein Synthesis of Bovine Mammary Epithelial Cells with Methionine Supply

The methionine (Met) uptake into mammary cells depends upon the corresponding amino acid (AA) transporters, which play a regulatory role in the mammary protein production beyond transport. Our previous studies have identified that seryl-tRNA synthetase (SARS) could be a novel mediator to regulate essential AA-stimulated casein synthesis in primary bovine mammary epithelial cells (BMECs). However, the regulatory mechanisms of Met in milk protein production in dairy cows remain further clarified. Here, we aimed to investigate the effects of Met on milk protein synthesis in BMECs and explore the underlying mechanism. The effects of Met on the AA transporter, casein synthesis, and the related signaling pathway were evaluated in the BMECs treated with 0.6 mM Met for 6 h combined with or without the inhibition of AA transporter (ASCT2, a neutral AA transporter) activity by the corresponding inhibitor (GPNA). Besides, the effects of SARS on the cells were mainly evaluated in the BMECs treated with 0.6 mM Met for 6 h together with or without SARS knockdown by RNAi interference. The gene expression of AA transporters and pathway-related genes were analyzed by the real-time quantitative polymerase chain reaction method, and the protein expression of related proteins were determined by the western blot assay. Results showed that 0.6 mM Met remarkably enhanced cell growth and β-casein synthesis compared to the supply of other Met concentrations. Among 13 amino acid transporters, 0.6 mM Met highly increased ASCT2 expression. This Met-stimulated ASCT2 expression and the enhanced mammary intracellular Met uptake were both decreased by the addition of 500 μM GPNA, an inhibitor of ASCT2. In the presence of 0.6 mM Met, the inhibition of ASCT2 activity (by GPNA) and SARS expression (by RNAi) both reduced β-casein synthesis. Additionally, 0.6 mM Met increased the gene expression of mTOR, S6K1, 4EBP1, and Akt; in contrast, the inhibition of ASCT2 by GPNA lowered the gene expression of these four genes. Collectively, this work suggests that ASCT2 is involved in the SARS-mediated Met stimulation of β-casein synthesis through enhancing mammary Met uptake and activating the mTOR signaling pathway in BMECs.

Identification of Milk Fat Metabolism-Related Pathways of the Bovine Mammary Gland during Mid and Late Lactation and Functional Verification of the ACSL4 Gene

The concentration of bovine milk fat changes regularly with lactation stages. In particular, milk fat percentage is higher in late lactation than mid lactation. Furthermore, milk fat composition is highly subject to a few genes. Thus, transcriptome sequencing was performed to explore the expression patterns of differentially-expressed genes (DEGs) in the parenchymal mammary gland of Holstein dairy cows between mid and late lactation. The 725 DEGs were screened (fold change > 2 and p-value < 0.05), and the peroxisome proliferator-activated receptor (PPAR) signaling pathway associated with lipid synthesis had a significant variation between the two periods (p-value < 0.05). The activation of the PPAR signal pathway may a key factor in the increasing of milk fat content in late lactation compared to mid lactation. Acyl-CoA synthetase long-chain family member 4 (ACSL4), a member of the PPAR signaling pathway, was upregulated in late lactation compared to mid lactation (p < 0.05). ACSL4 catalyzes the activation of long-chain fatty acids for cellular lipid synthesis. However, it remains uncertain that the molecular mechanism of milk fat synthesis is regulated by ACSL4 in dairy cows. Subsequently, the function verification of ACSL4 was performed in bovine mammary epithelial cells (BMECs). The upregulated expression of ACSL4 was accompanied by the increase of the concentration of intracellular triglycerides, whereas knockdown of ACSL4 decreased the concentration of intracellular triglycerides, which demonstrated that ACSL4 plays an important role in modulating milk fat synthesis. In conclusion, the results displayed that ACSL4 expression regulates triglyceride metabolism in ruminant mammary cells.

Multi-omics reveals functional genomic and metabolic mechanisms of milk production and quality in dairy cows

MOTIVATION

Enhancing the utilization of human-inedible crop by-products by ruminants to produce high-quality milk for human consumption is an emerging global task. We performed a multi-omics-based study to decipher the regulatory biological processes of milk production when cows fed low-quality crop by-products with the aim to improve their utilization.

RESULTS

Seven types of different high-throughput omics data were generated across three central organs [rumen, liver and mammary gland (MG)] and biofluids (rumen fluid and blood) that involved in milk production. The integrated multi-omics analysis including metabolomics, metagenomics and transcriptomics showed altered microbiome at compositional and functional levels, microbial metabolites in the rumen, down-regulated genes and associated functions in liver and MG. These changes simultaneously contributed to down-regulated three key metabolic nodes (propionate, glucose and amino acid) across these organs and biofluids that led to lowered milk yield and quality when cows consumed corn stover (CS). Hippuric acid was identified as a biomarker that led to low milk production in CS-fed cows, suggesting a future evaluation parameter related to the metabolic mechanism of low-quality forage utilization. This study unveils the milk production-related biological mechanism across different biofluids and tissues under a low-quality forage diet, which provides a novel understanding and potential improvement strategies for future crop by-products utilization and sustainable ruminant production.

AVAILABILITY AND IMPLEMENTATION

The raw files of metagenomics, metabolomics, and transcriptomics data can be accessed at NCBI SRA (No. SRR5028206), EMBI-EBI (No. MTBLS411), and GEO (NO. GSE78524) databases respectively.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

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.

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.

Use of gene expression profile to identify potentially relevant transcripts to myofibrillar fragmentation index trait

The myofibrillar fragmentation index (MFI) is an indicative trait for meat tenderness. Longissimus thoracis muscle samples from the 20 most extreme bulls (out of 80 bulls set) for MFI (high (n = 10) and low (n = 10) groups) trait were used to perform transcriptomic analysis, using RNA Sequencing (RNA-Seq). An average of 24.616 genes was expressed in the Nellore muscle transcriptome analysis. A total of 96 genes were differentially expressed (p value ≤ 0.001) between the two groups of divergent bulls for MFI. The HEBP2 and BDH1 genes were overexpressed in animals with high MFI. The MYBPH and MYL6, myosin encoders, were identified. The differentially expressed genes were related to increase mitochondria efficiency, especially in cells under oxidative stress conditions, and these also were related to zinc and calcium binding, membrane transport, and muscle constituent proteins, such as actin and myosin. Most of those genes were involved in metabolic pathways of oxidation-reduction, transport of lactate in the plasma membrane, and muscle contraction. This is the first study applying MFI phenotypes in transcriptomic studies to identify and understand differentially expressed genes for beef tenderness. These results suggest that differences detected in gene expression between high and low MFI animals are related to reactive mechanisms and structural components of oxidative fibers under the condition of cellular stress. Some genes may be selected as positional candidate genes to beef tenderness, MYL6, MYBPH, TRIM63, TRIM55, TRIOBP, and CHRNG genes. The use of MFI phenotypes could enhance results of meat tenderness studies.

Identification and comparison of exosomal microRNAs in the milk and colostrum of two different cow breeds

Bovine milk and colostrum provide essential nutrients and immunologically active factors that are beneficial to a newborn calf. Milk-and-colostrum-derived exosomes are known as the most important for cellular communication. Exosomes also contain non-coding RNA, such as microRNA. However, there is limited information about exosomal miRNA derived from the milk and colostrum of Holstein and DAK cattle. This study aimed to identify and characterize the exosomal microRNA in the milk and colostrum of Holstein and Doğu Anadolu Kirmizisi (DAK) cows. For this purpose, total RNA isolation was carried out on the milk and colostrum samples that were collected from the Holstein and DAK cattle breeds. The RNA samples were subjected to RNA sequencing and the microRNAs were determined. Lastly, gene ontology analysis was performed for target genes. A total of 795 miRNAs that were expressed differently were identified. A total of 545 of these were known miRNAs and 260 were found to be novel miRNAs. In the functional enrichment analysis, the miRNAs expressed in Holstein milk were mostly associated with milk synthesis, and those in colostrum were mostly involved in the immunity pathways. It was also observed that the miRNAs expressed in DAK milk regulated milk fat and protein metabolism, and there were miRNAs that regulated immune pathways in the colostrum. In addition to this, many novel miRNAs were defined in DAK colostrum. When the target genes of exosomal miRNA in Holstein and DAK milk and colostrum were compared, it was suggested that the DAK breed had genes that were mostly associated with the immune system. As a result, the data obtained from this study will provide beneficial contributions to potential miRNA biomarker studies for milk yield and mastitis.

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.

Nutrigenomic Role of Acetate and β-Hydroxybutyrate in Bovine Mammary Epithelial Cells

Acetate and β-hydroxybutyrate (BHBA) are the predominant substrates for de novo fatty acid (FA) synthesis in mammary gland of dairy cow. To investigate the nutrigenomic role of acetate and BHBA in bovine mammary epithelial cells during milk fat production, RNA sequencing (RNA-seq) transcriptomic analysis was used to identify differentially expressed genes (DEGs) between acetate- and BHBA-treated cells (high-milk fat cells) and control cells. A total of 625 DEGs (358 upregulated and 267 downregulated) were identified between the high-milk fat cells and control cells. Gene ontology enrichment analysis revealed that the upregulated genes in high-milk fat cells were mainly involved in lipid biosynthetic process, steroid biosynthetic process, oxidation-reduction process, receptor binding, and vesicle and small molecule biosynthetic process. The downregulated genes were mainly associated with immune response, cytokine production, negative regulation of biological process, and peptidyl-threonine modification. Pathway analysis indicated that FA metabolism and steroid biosynthesis were significantly enriched for the upregulated genes in the high-milk fat cells, while apoptosis was enriched for the downregulated genes. This work provides a profile of gene expression changes that occur during acetate- and BHBA-induced milk fat synthesis in bovine mammary epithelial cells, which furthers our understanding of the molecular regulation of lipid metabolism.

Seryl-tRNA synthetase is involved in methionine stimulation of β-casein synthesis in bovine mammary epithelial cells

Despite the well-characterised mechanisms of amino acids (AA) regulation of milk protein synthesis in mammary glands (MG), the underlying specific AA regulatory machinery in bovine MG remains further elucidated. As methionine (Met) is one of the most important essential and limiting AA for dairy cows, it is crucial to expand how Met exerts its regulatory effects on dairy milk protein synthesis. Our previous work detected the potential regulatory role of seryl-tRNA synthetase (SARS) in essential AA (EAA)-stimulated bovine casein synthesis. Here, we investigated whether and how SARS participates in Met stimulation of casein production in bovine mammary epithelial cells (BMEC). With or without RNA interference against SARS, BMEC were treated with the medium in the absence (containing all other EAA and devoid of Met alone)/presence (containing 0·6 mm of Met in the medium devoid of Met alone) of Met. The protein abundance of β-casein and members of the mammalian target of rapamycin (mTOR) and general control nonderepressible 2 (GCN2) pathways was determined by immunoblot assay after 6 h treatment, the cell viability and cell cycle progression were determined by cell counting and propidium iodide-staining assay after 24 h treatment, and protein turnover was determined by l-[ring-3H5]phenylalanine isotope tracing assay after 48 h treatment. In the absence of Met, there was a general reduction in cell viability, total protein synthesis and β-casein production; in contrast, total protein degradation was enhanced. SARS knockdown strengthened these changes. Finally, SARS may work to promote Met-stimulated β-casein synthesis via affecting mTOR and GCN2 routes in BMEC.

The Selection of Reference Genes for Quantitative Real-Time PCR in the Ashidan Yak Mammary Gland During Lactation and Dry Period

Simple Summary

The Ashidan yak is a new cultivated breed which has polled characteristics and a mild temperament. Improving milk yield is an important aspect of a breeding program for this breed. The mammary gland undergoes dramatic physiological and metabolic changes during the transition from lactation to dry periods, which involves the expression and regulation of a great number of genes. Quantification of gene expression levels by real-time quantitative polymerase chain reaction (RT-qPCR) is important to reveal the molecular mechanisms of mammary gland development and lactation. The accuracy of RT-qPCR is strongly influenced by the expression stability of reference genes, however, a systematic approach for selecting reference genes used for analyzing gene expression of the Ashidan yak has not been developed. In this study, we selected reference genes and analyzed their expression stability at different physiological stages (lactation and dry period). We found the hydroxymethylbilane synthase gene (HMBS) and the tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein, zeta polypeptide gene (YWHAZ) were the most stable genes of the mammary gland of the Ashidan yak. These results help to improve the accuracy of gene expression analysis and provide a basis for future functional studies of target gene expression in the mammary gland of the Ashidan yak.

Abstract

Investigating the critical genes related to milk synthesis is essential for the improvement of the milk yield of the yak. Real-time quantitative polymerase chain reaction (RT-qPCR) is a reliable and widely used method to measure and evaluate gene expression levels. Selection of suitable reference genes is mandatory to acquire accurate normalization of gene expression results from RT-qPCR. To select the most stable reference genes for reliable normalization of mRNA expression by RT-qPCR in the mammary gland of the Ashidan yak, we selected 16 candidate reference genes and analyzed their expression stability at different physiological stages (lactation and dry period). The expression stability of the candidate reference genes was assessed using geNorm, NormFinder, BestKeeper, Delta Ct, and RefFinder methods. The results showed that the hydroxymethylbilane synthase gene (HMBS) and the tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein, zeta polypeptide gene (YWHAZ) were the most stable genes across all treatment samples. The reliability of selected reference genes was validated by normalizing relative expression of the lactation-related 60S ribosomal protein L35 gene (RPL35). The relative expression of RPL35 varied considerably according to the different reference genes. This work provides valuable information to further promote research in the molecular mechanisms involved in lactation and mammary gland development and provides a foundation for the improvement of the milk yield and quality of the Ashidan yak.

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

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

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

BACKGROUND

It is known that long non-coding RNAs (lncRNAs) play an important role in various biological processes, including cell proliferation, differentiation and apoptosis. However, their functions and profiles in lactation cycle of dairy cows are largely unknown. In this study, lncRNA-seq technique was employed to compare the expression profiles of lncRNAs and mRNAs from Chinese Holstein mammary gland in dry and lactation period.

RESULT

Totally 3746 differentially expressed lncRNAs (DELs) and 2890 differentially expressed genes (DEGs) were identified from the dry and lactation mammary glands of Holstein cows. Functional enrichment analysis on target genes of lncRNAs indicated that these genes were involved in lactation-related signaling pathways, including cell cycle, JAK-STAT, cell adhesion, and PI3K-Akt signaling pathways. Additionally, the interaction between lncRNAs and their potential miRNAs was predicted and partly verified. The result indicated that the lactation-associated miR-221 might interact with lncRNAs TCONS_00040268, TCONS_00137654, TCONS_00071659 and TCONS_00000352, which revealed that these lncRNAs might be important regulators for lactation cycle.

CONCLUSION

This study provides a resource for lncRNA research on lactation cycle of bovine mammary gland. Besides, the interaction between lncRNAs and the specific miRNA is revealed. It expands our knowledge about lncRNA and miRNA biology as well as contributes to clarify the regulation of lactation cycle of bovine mammary gland.

The Use of "Omics" in Lactation Research in Dairy Cows

“Omics” is the application of genomics, transcriptomics, proteomics, and metabolomics in biological research. Over the years, tremendous amounts of biological information has been gathered regarding the changes in gene, mRNA and protein expressions as well as metabolites in different physiological conditions and regulations, which has greatly advanced our understanding of the regulation of many physiological and pathophysiological processes. The aim of this review is to comprehensively describe the advances in our knowledge regarding lactation mainly in dairy cows that were obtained from the “omics” studies. The “omics” technologies have continuously been preferred as the technical tools in lactation research aiming to develop new nutritional, genetic, and management strategies to improve milk production and milk quality in dairy cows.