Mammary tissue transcriptomic analysis for construction of integrated regulatory networks involved in lactogenesis of Ovis aries

Transcriptomic analysis of mammary gland tissues in lactating and non-lactating dairy goats reveals miRNA-mediated regulation of lactation, involution, and remodeling

Background

Dynamic changes in the mammary gland during lactation and the dry period involve proliferation, secretion, apoptosis, and remodeling of mammary epithelial cells. MicroRNAs (miRNAs) are recognized as critical regulators of mammary gland development and lactation. However, their expression patterns and regulatory mechanisms at different lactation stages-particularly during mammary involution and remodeling-remain poorly understood in dairy goats.

Methods

In this study, high-throughput sequencing was employed to analyze miRNA expression profiles in goat mammary tissues at five key stages: late gestation (LG), early lactation (EL), peak lactation (PL), late lactation (LL), and the dry period (DP). Differential expression analysis, miRNA clustering, Gene Ontology (GO) annotation, and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment were performed to predict the functions of target genes. A miRNA-mRNA regulatory network associated with mammary gland development was constructed, and functional validation experiments were conducted to confirm key regulatory relationships.

Results

A total of 1,120 miRNAs were identified, including 408 known and 712 newly predicted miRNAs. Among them, 383 were significantly differentially expressed, with the largest number observed between the dry period and late gestation. Six expression-specific miRNA clusters were identified. Functional enrichment analysis indicated that these miRNAs may regulate epithelial cell proliferation, apoptosis, and tissue remodeling by targeting pathways such as energy metabolism, cell adhesion, and the PI3K/Akt signaling pathway. IGF1R was identified as a key regulatory gene in the miRNA-mRNA network related to mammary gland development. Experimental validation showed that chi-miR-423-3p inhibited mammary epithelial cell proliferation, induced G1/S cell cycle arrest, and promoted apoptosis by targeting IGF1R and suppressing the PI3K/Akt pathway.

Conclusion

This study highlights the dynamic regulatory roles of miRNAs in the goat mammary gland across lactation stages. Notably, the miR-423-3p/IGF1R axis is a key regulator of mammary remodeling during the dry period, offering new insights into the molecular basis of mammary gland functional transitions.

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.

Transcriptome analysis of the Bactrian camel (Camelus bactrianus) reveals candidate genes affecting milk production traits

BACKGROUND

Milk production traits are complex traits with vital economic importance in the camel industry. However, the genetic mechanisms regulating milk production traits in camels remain poorly understood. Therefore, we aimed to identify candidate genes and metabolic pathways that affect milk production traits in Bactrian camels.

METHODS

We classified camels (fourth parity) as low- or high-yield, examined pregnant camels using B-mode ultrasonography, observed the microscopic changes in the mammary gland using hematoxylin and eosin (HE) staining, and used RNA sequencing to identify differentially expressed genes (DEGs) and pathways.

RESULTS

The average standard milk yield over the 300 days during parity was recorded as 470.18 ± 9.75 and 978.34 ± 3.80 kg in low- and high-performance camels, respectively. Nine female Junggar Bactrian camels were subjected to transcriptome sequencing, and 609 and 393 DEGs were identified in the low-yield vs. high-yield (WDL vs. WGH) and pregnancy versus colostrum period (RSQ vs. CRQ) comparison groups, respectively. The DEGs were compared with genes associated with milk production traits in the Animal Quantitative Trait Loci database and in Alashan Bactrian camels, and 65 and 46 overlapping candidate genes were obtained, respectively. Functional enrichment and protein-protein interaction network analyses of the DEGs and candidate genes were conducted. After comparing our results with those of other livestock studies, we identified 16 signaling pathways and 27 core candidate genes associated with maternal parturition, estrogen regulation, initiation of lactation, and milk production traits. The pathways suggest that emerged milk production involves the regulation of multiple complex metabolic and cellular developmental processes in camels. Finally, the RNA sequencing results were validated using quantitative real-time PCR; the 15 selected genes exhibited consistent expression changes.

CONCLUSIONS

This study identified DEGs and metabolic pathways affecting maternal parturition and milk production traits. The results provides a theoretical foundation for further research on the molecular mechanism of genes related to milk production traits in camels. Furthermore, these findings will help improve breeding strategies to achieve the desired milk yield in camels.

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

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

Identification of Unique Key miRNAs, TFs, and mRNAs in Virulent MTB Infection Macrophages by Network Analysis

Although Mycobacterium tuberculosis (MTB) has existed for thousands of years, its immune escape mechanism remains obscure. Increasing evidence signifies that microRNAs (miRNAs) play pivotal roles in the progression of tuberculosis (TB). RNA sequencing was used to sequence miRNAs in human acute monocytic leukemia cells (THP-1) infected by the virulent MTB-1458 strain and the avirulent vaccine strain Mycobacterium bovis Bacillus Calmette-Guérin (BCG). Sets of differentially expressed miRNAs (DE-miRNAs) between MTB-1458/BCG-infected groups and uninfected groups were identified, among which 18 were differentially expressed only in the MTB-1458-infected THP-1 group. Then, 13 transcription factors (TFs) and 81 target genes of these 18 DE-miRNAs were matched. Gene Ontology classification as well as Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis showed that the candidate targets were predominantly involved in apoptotic-associated and interferon-γ-mediated signaling pathways. A TF-miRNA-mRNA interaction network was constructed to analyze the relationships among these 18 DE-miRNAs and their targets and TFs, as well as display the hub miRNAs, TFs, and target genes. Considering the degrees from network analysis and the reported functions, this study focused on the BHLHE40-miR-378d-BHLHE40 regulation axis and confirmed that BHLHE40 was a target of miR-378d. This cross-talk among DE-miRNAs, mRNAs, and TFs might be an important feature in TB, and the findings merited further study and provided new insights into immune defense and evasion underlying host-pathogen interactions.