Integrative Analysis of miRNA and mRNA Expression Profiles in Mammary Glands of Holstein Cows Artificially Infected with Staphylococcus aureus

Transcriptome analysis reveals the regulation of miR-19b on inflammation in bovine mammary epithelial cells

MicroRNAs (miRNAs) are involved in various biological processes where they regulate the expression of mRNAs. Bovine mammary epithelial cells (bMECs) are functional cells that mediate mammary inflammatory immunity. Although numerous miRNAs regulate the function of bMECs, the role of miR-19b in bMECs has not been reported. In this study, the transcriptome of miR-19b overexpressed bMECs was analyzed by RNA-seq. Additionally, the differentially expressed genes (DEGs) were analyzed to establish the role of miR-19b in bMECs. The results revealed 269 DEGs between the miR-19b overexpression group and the negative control, including 199 up-regulated and 70 down-regulated genes. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses revealed that the DEGs regulated immune and inflammatory responses through Staphylococcus aureus (S. aureus) infection and phosphatidylinositol 3-kinase (PI3K)-Akt signaling pathway. In addition, the expression of miR-19b was significantly upregulated in lipophosphoric acid (LTA)-induced bMECs, and overexpression of miR-19b negatively regulated the expression of inflammatory cytokines IL-1β and IL-6, thereby alleviating the inflammatory response of LTA-induced bMECs. Based on the above results, we speculate that miR-19b may inhibit in dairy cow mammary inflammation caused by S. aureus, and this process may be mediated through the regulation of relevant gene expression and signaling pathways. The findings from this study provide a new reference for analyzing the molecular regulation of miR-19b in bMECs.

Multi-omics integration identifies regulatory factors underlying bovine subclinical mastitis

BACKGROUND

Mastitis caused by multiple factors remains one of the most common and costly disease of the dairy industry. Multi-omics approaches enable the comprehensive investigation of the complex interactions between multiple layers of information to provide a more holistic view of disease pathogenesis. Therefore, this study investigated the genomic and epigenomic signatures and the possible regulatory mechanisms underlying subclinical mastitis by integrating RNA sequencing data (mRNA and lncRNA), small RNA sequencing data (miRNA) and DNA methylation sequencing data of milk somatic cells from 10 healthy cows and 20 cows with naturally occurring subclinical mastitis caused by Staphylococcus aureus or Staphylococcus chromogenes.

RESULTS

Functional investigation of the data sets through gene set analysis uncovered 3458 biological process GO terms and 170 KEGG pathways with altered activities during subclinical mastitis, provided further insights into subclinical mastitis and revealed the involvement of multi-omics signatures in the altered immune responses and impaired mammary gland productivity during subclinical mastitis. The abundant genomic and epigenomic signatures with significant alterations related to subclinical mastitis were observed, including 30,846, 2552, 1276 and 57 differential methylation haplotype blocks (dMHBs), differentially expressed genes (DEGs), lncRNAs (DELs) and miRNAs (DEMs), respectively. Next, 5 factors presenting the principal variation of differential multi-omics signatures were identified. The important roles of Factor 1 (DEG, DEM and DEL) and Factor 2 (dMHB and DEM), in the regulation of immune defense and impaired mammary gland functions during subclinical mastitis were revealed. Each of the omics within Factors 1 and 2 explained about 20% of the source of variation in subclinical mastitis. Also, networks of important functional gene sets with the involvement of multi-omics signatures were demonstrated, which contributed to a comprehensive view of the possible regulatory mechanisms underlying subclinical mastitis. Furthermore, multi-omics integration enabled the association of the epigenomic regulatory factors (dMHBs, DELs and DEMs) of altered genes in important pathways, such as 'Staphylococcus aureus infection pathway' and 'natural killer cell mediated cytotoxicity pathway', etc., which provides further insights into mastitis regulatory mechanisms. Moreover, few multi-omics signatures (14 dMHBs, 25 DEGs, 18 DELs and 5 DEMs) were identified as candidate discriminant signatures with capacity of distinguishing subclinical mastitis cows from healthy cows.

CONCLUSION

The integration of genomic and epigenomic data by multi-omics approaches in this study provided a better understanding of the molecular mechanisms underlying subclinical mastitis and identified multi-omics candidate discriminant signatures for subclinical mastitis, which may ultimately lead to the development of more effective mastitis control and management strategies.

Genomic clues of association between clinical mastitis and SNPs identified by ddRAD sequencing in Murrah buffaloes

The total milk production of India is 209.96 MT out of which 45% is contributed by the indigenous buffalo and due to their high producing virtue, the prevalence of mastitis is 5-20%. Despite the increasing level of technological advancement, mastitis is still an issue of concern for dairy industry in India as well as across the world. Therefore, the present study aimed to identify the SNPs and associate them with the incidence of clinical mastitis in Murrah buffalo using the ddRAD sequencing approach taking mastitis incidence data of 96 Murrah buffaloes. A total of 246 million quality controlled reads were obtained with an average alignment rate of 99.01% and at a read depth of 10, quality controlled SNPs obtained were 18,056. The logistic regression model was used and a total of seven SNPs were found significantly associated (p < 0.001) with mastitis incidence and seven genes were identified viz., NCBP1, FOXN3, TPK1, XYLT2, CPXM2, HERC1, and OPCML. The majority of them were having tumor suppressing action, related to immunogenetics or glycolytic and energy production. Conclusively, the SNPs identified in this study may be useful for future studies on mastitis incidence in Murrah buffalo and the SNP associations can be further validated.

Genome-wide diversity and admixture of five indigenous cattle populations from the Tigray region of northern Ethiopia

The Tigray region, where we found around eight per cent of the indigenous cattle population of Ethiopia, is considered as the historic centre of the country, with the ancient pre-Aksumite and Aksumite civilisations in contact with the civilisations of the Fertile Crescent and the Indian subcontinent. Here, we used whole genome sequencing data to characterise the genomic diversity, relatedness, and admixture of five cattle populations (Abergelle, Arado, Begait, Erob, and Raya) indigenous to the Tigray region of Ethiopia. We detected 28 to 29 million SNPs and 2.7 to 2.9 million indels in each population, of which 7% of SNPs and 34% of indels were novel. Functional annotation of the variants showed around 0.01% SNPs and 0.22%-0.27% indels in coding regions. Enrichment analysis of genes overlapping missense private SNPs revealed 20 significant GO terms and KEGG pathways that were shared by or specific to breeds. They included important genes associated with morphology (SCN4A, TAS1R2 and KCNG4), milk yield (GABRG1), meat quality (MMRN2, VWC2), feed efficiency (PCDH8 and SLC26A3), immune response (LAMC1, PCDH18, CELSR1, TLR6 and ITGA5), heat resistance (NPFFR1 and HTR7) and genes belonging to the olfactory gene family, which may be related to adaptation to harsh environments. Tigray indigenous cattle are very diverse. Their genome-wide average nucleotide diversity ranged from 0.0035 to 0.0036. The number of heterozygous SNPs was about 0.6-0.7 times higher than homozygous ones. The within-breed average number of ROHs ranged from 777.82 to 1000.45, with the average sum of the length of ROHs ranging from 122.01 Mbp to 163.88 Mbp. The genomic inbreeding coefficients differed among animals and breeds, reaching up to 10% in some Begait and Raya animals. Tigray indigenous cattle shared a common ancestry with Asian indicine (85.6%-88.7%) and African taurine (11.3%-14.1%) cattle, with very small, if any, European taurine introgression. This study identified high within-breed genetic diversity representing an opportunity for breeding improvement programs and, also, significant novel variants that could increase the number of known cattle variants, an important contribution to the knowledge of domestic cattle genetic diversity.

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.

Emerging roles of noncoding micro RNAs and circular RNAs in bovine mastitis: Regulation, breeding, diagnosis, and therapy

Bovine mastitis is one of the most troublesome and costly problems in the modern dairy industry, which is not only difficult to monitor, but can also cause economic losses while having significant implications on public health. However, efficacious preventative methods and therapy are still lacking. Moreover, new drugs and therapeutic targets are in increasing demand due to antibiotic restrictions. In recent years, noncoding RNAs have gained popularity as a topic in pathological and genetic studies. Meanwhile, there is growing evidence that they play a role in regulating various biological processes and developing novel treatment platforms. In light of this, this review focuses on two types of noncoding RNAs, micro RNAs and circular RNAs, and summarizes their characterizations, relationships, potential applications as selection markers, diagnostic or treatment targets and potential applications in RNA-based therapy, in order to shed new light on further research.

Analyzing the interactions of mRNAs, miRNAs and lncRNAs to predict ceRNA networks in bovine cystic follicular granulosa cells

Cross-talk between competitive endogenous RNAs (ceRNAs) may play a critical role in revealing potential mechanism of bovine follicular cysts. Ovarian cyst has always been an intractable scientific problem and has led to considerable economic losses to bovine breeding industry. However, its pathogenesis and molecular mechanisms are still not well understood. Here, this study aimed to investigate the role of non-coding RNAs (ncRNAs) and the ceRNA networks in bovine follicular cyst. Whole transcriptome sequencing of bovine follicular granulosa cells (GCs) was conducted to obtain the expression profiles of mRNAs, lncRNAs and miRNAs. The results for the identified expressions of 8,003 mRNAs, 579 lncRNAs and 205 miRNAs were often altered between cystic and normal follicular GCs. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses were performed on these differentially expressed mRNAs. Furthermore, the ceRNA network combining mRNAs, miRNAs, and lncRNAs using several bioinformatics methods based on co-expression analysis between the differentially expressed RNAs was conducted. Finally, the lncRNA NONBTAT027373.1-miR-664b-HSD17B7 pathway was verified by dual-luciferase reporting assay and RNA binding protein immunoprecipitation (RIP) assay. LncRNA NONBTAT027373.1 sponged miR-664b in GCs and prevented miR-664b from binding to the HSD17B7 3′-UTR. These results indicated that genes and lncRNAs related to steroid hormone synthesis and energy metabolism could play important roles in the formation of bovine cystic follicles through the ceRNA mechanism and represent candidate targets for further research. This can be used as a practical guideline for promoting healthy and highly efficient development in the bovine industry.

Integrated Analysis of Transcriptome mRNA and miRNA Profiles Reveals Self-Protective Mechanism of Bovine MECs Induced by LPS

Many studies have investigated the molecular crosstalk between mastitis-pathogens and cows by either miRNA or mRNA profiles. Here, we employed both miRNA and mRNA profiles to understand the mechanisms of the response of bovine mammary epithelial cells (bMECs) to lipopolysaccharide (LPS) by RNA-Seq. The total expression level of miRNAs increased while mRNAs reduced after LPS treatment. About 41 differentially expressed mRNAs and 45 differentially expressed miRNAs involved in inflammation were screened out. We found the NFκB-dependent chemokine, CXCL1, CXCL3, CXCL6, IL8, and CX3CL1 to be strongly induced. The anti-apoptosis was active because BCL2A1 and BIRC3 significantly increased with a higher expression. The effects of anti-microbe and inflammation were weakly activated because TNF, IL1, CCL20, CFB, S100A, MMP9, and NOS2A significantly increased but with a low expression, IL6 and β-defensin decreased. These activities were supervised by the NFKBIA to avoid excessive damage to bMECs. The bta-let-7a-5p, bta-miR-30a-5p, bta-miR-125b, and bta-miR-100 were essential to regulate infection process in bMECs after LPS induction. Moreover, the lactation potential of bMECs was undermined due to significantly downregulated SOSTDC1, WNT7B, MSX1, and bta-miR-2425-5p. In summary, bMECs may not be good at going head-to-head with the pathogens; they seem to be mainly charged with sending out signals for help and anti-apoptosis for maintaining lives after LPS induction.

Identification of Bovine miRNAs with the Potential to Affect Human Gene Expression

Milk and other products from large mammals have emerged during human evolution as an important source of nutrition. Recently, it has been recognized that exogenous miRNAs (mRNA inhibited RNA) contained in milk and other tissues of the mammalian body can enter the human body, which in turn have the ability to potentially regulate human metabolism by affecting gene expression. We studied for exogenous miRNAs from Bos taurus that are potentially contain miRNAs from milk and that could act postprandially as regulators of human gene expression. The interaction of 17,508 human genes with 1025 bta-miRNAs, including 245 raw milk miRNAs was studied. The milk bta-miR-151-5p, bta-miR-151-3p, bta-miRNA-320 each have 11 BSs (binding sites), and bta-miRNA-345-5p, bta-miRNA-614, bta-miRNA-1296b and bta-miRNA-149 has 12, 14, 15 and 26 BSs, respectively. The bta-miR-574-5p from cow’s milk had 209 human genes in mRNAs from one to 25 repeating BSs. We found 15 bta-miRNAs that have 100% complementarity to the mRNA of 13 human target genes. Another 12 miRNAs have BSs in the mRNA of 19 human genes with 98% complementarity. The bta-miR-11975, bta-miR-11976, and bta-miR-2885 BSs are located with the overlap of nucleotide sequences in the mRNA of human genes. Nucleotide sequences of BSs of these miRNAs in 5′UTR mRNA of human genes consisted of GCC repeats with a total length of 18 nucleotides (nt) in 18 genes, 21 nt in 11 genes, 24 nt in 14 genes, and 27–48 nt in nine genes. Nucleotide sequences of BSs of bta-miR-11975, bta-miR-11976, and bta-miR-2885 in CDS mRNA of human genes consisted of GCC repeats with a total length of 18 nt in 33 genes, 21 nt in 13 genes, 24 nt in nine genes, and 27–36 nt in 11 genes. These BSs encoded polyA or polyP peptides. In only one case, the polyR (SLC24A3 gene) was encoded. The possibility of regulating the expression of human genes by exogenous bovine miRNAs is discussed.

The Role of microRNAs in the Mammary Gland Development, Health, and Function of Cattle, Goats, and Sheep

Milk is an integral and therefore complex structural element of mammalian nutrition. Therefore, it is simple to conclude that lactation, the process of producing milk, is as complex as the mammary gland, the organ responsible for this biochemical activity. Nutrition, genetics, epigenetics, disease pathogens, climatic conditions, and other environmental variables all impact breast productivity. In the last decade, the number of studies devoted to epigenetics has increased dramatically. Reports are increasingly describing the direct participation of microRNAs (miRNAs), small noncoding RNAs that regulate gene expression post-transcriptionally, in the regulation of mammary gland development and function. This paper presents a summary of the current state of knowledge about the roles of miRNAs in mammary gland development, health, and functions, particularly during lactation. The significance of miRNAs in signaling pathways, cellular proliferation, and the lipid metabolism in agricultural ruminants, which are crucial in light of their role in the nutrition of humans as consumers of dairy products, is discussed.

Bovine milk somatic cell transcriptomic response to Staphylococcus aureus is dependent on strain genotype

Background

Mastitis is an economically important disease of dairy cows with Staphylococcus aureus a major cause worldwide. Challenge of Holstein-Friesian cows demonstrated that S. aureus strain MOK124, which belongs to Clonal Complex (CC)151, caused clinical mastitis, while strain MOK023, belonging to CC97, caused mild or subclinical mastitis. The aim of this study was to elucidate the molecular mechanisms of the host immune response utilising a transcriptomic approach. Milk somatic cells were collected from cows infected with either S. aureus MOK023 or MOK124 at 0, 24, 48, 72 and 168 h post-infection (hpi) and analysed for differentially expressed (DE) genes in response to each strain.

Results

In response to MOK023, 1278, 2278, 1986 and 1750 DE genes were found at 24, 48, 72 and 168 hpi, respectively, while 2293, 1979, 1428 and 1544 DE genes were found in response to MOK124 at those time points. Genes involved in milk production (CSN1, CSN10, CSN1S2, CSN2, a-LACTA and PRLR) were downregulated in response to both strains, with a more pronounced decrease in the MOK124 group. Immune response pathways such as NF-κB and TNF signalling were overrepresented in response to both strains at 24 hpi. These immune pathways continued to be overrepresented in the MOK023 group at 48 and 72 hpi, while the Hippo signalling, extracellular matrix interaction (ECM) and tight junction pathways were overrepresented in the MOK124 group between 48 and 168 hpi. Cellular composition analysis demonstrated that a neutrophil response was predominant in response to MOK124, while M1 macrophages were the main milk cell type post-infection in the MOK023 group.

Conclusions

A switch from immune response pathways to pathways involved in maintaining the integrity of the epithelial cell layer was observed in the MOK124 group from 48 hpi, which coincided with the occurrence of clinical signs in the infected animals. The higher proportion of M1 macrophages in the MOK023 group and lack of substantial neutrophil recruitment in response to MOK023 may indicate immune evasion by this strain. The results of this study highlight that the somatic cell transcriptomic response to S. aureus is dependent on the genotype of the infecting strain.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08135-7.