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.