Identification of novel molecular markers of mastitis caused by Staphylococcus aureus using gene expression profiling in two consecutive generations of Chinese Holstein dairy cattle

Decoding epigenetic markers: implications of traits and genes through DNA methylation in resilience and susceptibility to mastitis in dairy cows

Cattle farming faces challenges linked to intensive exploitation and climate change, requiring the reinforcement of animal resilience in response to these dynamic environments. Currently, genetic selection is used to enhance resilience by identifying animals resistant to specific diseases; however, certain diseases, such as mastitis, pose difficulties in genetic prediction. This study introduced the utilization of enzymatic methyl sequencing (EM-seq) of the blood genomic DNA from twelve dairy cows to identify DNA methylation biomarkers, with the aim of predicting resilience and susceptibility to mastitis. The analysis uncovered significant differences between cows resilient and susceptible to mastitis, with 196,275 differentially methylated cytosines (DMCs) and 1,227 Differentially Methylated Regions (DMRs). Key genes associated with the immune response and morphological traits, including ENOPH1, MYL10 and KIR2DL5A, were identified by our analysis. Quantitative trait loci (QTL) were also highlighted and the body weight trait was the most targeted by DMCs and DMRs. Based on our results, the risk of developing mastitis can potentially be estimated with as few as fifty methylation biomarkers, paving the way for early animal selection. This research sets the stage for improved animal health management and economic yields within the framework of agricultural sustainability through early selection based on the epigenetic status of animals.

UID-Dual Transcriptome Sequencing Analysis of the Molecular Interactions between Streptococcus agalactiae ATCC 27956 and Mammary Epithelial Cells

Streptococcus agalactiae ATCC 27956 is a highly contagious Gram-positive bacterium that causes mastitis, has a high infectivity for mammary epithelial cells, and becomes challenging to treat. However, the molecular interactions between it and mammary epithelial cells remain poorly understood. This study analyzed differential gene expression in mammary epithelial cells with varying levels of S. agalactiae infection using UID-Dual transcriptome sequencing and bioinformatics tools. This study identified 211 differentially expressed mRNAs (DEmRNAs) and 452 differentially expressed lncRNAs (DElncRNAs) in host cells, primarily enriched in anti-inflammatory responses, immune responses, and cancer-related processes. Additionally, 854 pathogen differentially expressed mRNAs (pDEmRNAs) were identified, mainly enriched in protein metabolism, gene expression, and biosynthesis processes. Mammary epithelial cells activate pathways, such as the ERK1/2 pathway, to produce reactive oxygen species (ROS) to eliminate bacteria. The bacteria disrupt the host's innate immune mechanisms by interfering with the alternative splicing processes of mammary epithelial cells. Specifically, the bacterial genes of tsf, prfB, and infC can interfere with lncRNAs targeting RUNX1 and BCL2L11 in mammary epithelial cells, affecting the alternative splicing of target genes and altering normal molecular regulation.

Identification and characterization of whole blood gene expression and splicing quantitative trait loci during early to mid-lactation of dairy cattle

BACKGROUND

Characterization of regulatory variants (e.g., gene expression quantitative trait loci, eQTL; gene splicing QTL, sQTL) is crucial for biologically interpreting molecular mechanisms underlying loci associated with complex traits. However, regulatory variants in dairy cattle, particularly in specific biological contexts (e.g., distinct lactation stages), remain largely unknown. In this study, we explored regulatory variants in whole blood samples collected during early to mid-lactation (22-150 days after calving) of 101 Holstein cows and analyzed them to decipher the regulatory mechanisms underlying complex traits in dairy cattle.

RESULTS

We identified 14,303 genes and 227,705 intron clusters expressed in the white blood cells of 101 cattle. The average heritability of gene expression and intron excision ratio explained by cis-SNPs is 0.28 ± 0.13 and 0.25 ± 0.13, respectively. We identified 23,485 SNP-gene expression pairs and 18,166 SNP-intron cluster pairs in dairy cattle during early to mid-lactation. Compared with the 2,380,457 cis-eQTLs reported to be present in blood in the Cattle Genotype-Tissue Expression atlas (CattleGTEx), only 6,114 cis-eQTLs (P < 0.05) were detected in the present study. By conducting colocalization analysis between cis-e/sQTL and the results of genome-wide association studies (GWAS) from four traits, we identified a cis-e/sQTL (rs109421300) of the DGAT1 gene that might be a key marker in early to mid-lactation for milk yield, fat yield, protein yield, and somatic cell score (PP4 > 0.6). Finally, transcriptome-wide association studies (TWAS) revealed certain genes (e.g., FAM83H and TBC1D17) whose expression in white blood cells was significantly (P < 0.05) associated with complex traits.

CONCLUSIONS

This study investigated the genetic regulation of gene expression and alternative splicing in dairy cows during early to mid-lactation and provided new insights into the regulatory mechanisms underlying complex traits of economic importance.

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.

Chlorogenic acid enhances PPARγ-mediated lipogenesis through preventing Lipin 1 nuclear translocation in Staphylococcus aureus-exposed bovine mammary epithelial cells

Chlorogenic acid (CGA) as one of the most ubiquitously dietary polyphenolic compounds, has been reported to have various antimicrobial effects and exhibit strong anti-inflammatory ability. Staphylococcus aureus is a gram-positive bacterium that can induce mastitis. However, the mechanism through which S. aureus infection affects lipid synthesis and whether CGA have protective effect on S. aureus reduced lipid synthesis is not fully understood. In this study, the internalization of S. aureus reduced intracellular lipid droplet formation, decreased the levels of intracellular triacylglycerol, total cholesterol and 7 types of fatty acid and downregulated the expression of lipogenic genes FAS, ACC, and DGAT1 in bovine mammary epithelial cells (BMECs). In addition, we found that S. aureus intracellular infection attenuated mTORC1 activation resulting in Lipin 1 nuclear localization. Remarkablely, S. aureus infection-mediated repression of lipid synthesis related to the mTORC1 signaling and Lipin 1 nuclear localization can be alleviated by CGA. Thus, our findings provide a novel mechanism by which lipid synthesis is regulated under S. aureus infection and the protective effects of CGA on lipid synthesis in BMECs.

Identification of hub genes associated with somatic cell score in dairy cow

CONTEXT

Somatic cell count (SCC) is used as an indicator of udder health. The log transformation of SCC is called somatic cell score (SCS).

AIM

Several QTL and genes have been identified that are associated with SCS. This study aimed to identify the most important genes associated with SCS.

METHODS

This study compiled 168 genes that were reported to be significantly linked to SCS. Pathway analysis and network analysis were used to identify hub genes.

KEY RESULTS

Pathway analysis of these genes identified 73 gene ontology (GO) terms associated with SCS. These GO terms are associated with molecular function, biological processes, and cellular components, and the identified pathways are directly or indirectly linked with the immune system. In this study, a gene network was constructed, and from this network, the 17 hub genes (CD4, CXCL8, TLR4, STAT1, TLR2, CXCL9, CCR2, IGF1, LEP, SPP1, GH1, GHR, VWF, TNFSF11, IL10RA, NOD2, and PDGFRB) associated to SCS were identified. The subnetwork analysis yielded 10 clusters, with cluster 1 containing all identified hub genes (except for the VWF gene).

CONCLUSION

Most hub genes and pathways identified in our study were mainly involved in inflammatory and cytokine responses.

IMPLICATIONS

Result obtained in current study provides knowledge of the genetic basis and biological mechanisms controlling SCS. Therefore, the identified hub genes may be regarded as the main gene for the genomic selection of mastitis resistance.

Integrated analysis of inflammatory mRNAs, miRNAs, and lncRNAs elucidates the molecular interactome behind bovine mastitis

Mastitis is known as intramammary inflammation, which has a multifactorial complex phenotype. However, the underlying molecular pathogenesis of mastitis remains poorly understood. In this study, we utilized a combination of RNA-seq and miRNA-seq techniques, along with computational systems biology approaches, to gain a deeper understanding of the molecular interactome involved in mastitis. We retrieved and processed one hundred transcriptomic libraries, consisting of 50 RNA-seq and 50 matched miRNA-seq data, obtained from milk-isolated monocytes of Holstein-Friesian cows, both infected with Streptococcus uberis and non-infected controls. Using the weighted gene co-expression network analysis (WGCNA) approach, we constructed co-expressed RNA-seq-based and miRNA-seq-based modules separately. Module-trait relationship analysis was then performed on the RNA-seq-based modules to identify highly-correlated modules associated with clinical traits of mastitis. Functional enrichment analysis was conducted to understand the functional behavior of these modules. Additionally, we assigned the RNA-seq-based modules to the miRNA-seq-based modules and constructed an integrated regulatory network based on the modules of interest. To enhance the reliability of our findings, we conducted further analyses, including hub RNA detection, protein-protein interaction (PPI) network construction, screening of hub-hub RNAs, and target prediction analysis on the detected modules. We identified a total of 17 RNA-seq-based modules and 3 miRNA-seq-based modules. Among the significant highly-correlated RNA-seq-based modules, six modules showed strong associations with clinical characteristics of mastitis. Functional enrichment analysis revealed that the turquoise module was directly related to inflammation persistence and mastitis development. Furthermore, module assignment analysis demonstrated that the blue miRNA-seq-based module post-transcriptionally regulates the turquoise RNA-seq-based module. We also identified a set of different RNAs, including hub-hub genes, hub-hub TFs (transcription factors), hub-hub lncRNAs (long non-coding RNAs), and hub miRNAs within the modules of interest, indicating their central role in the molecular interactome underlying the pathogenic mechanisms of S. uberis infection. This study provides a comprehensive insight into the molecular crosstalk between immunoregulatory mRNAs, miRNAs, and lncRNAs during S. uberis infection. These findings offer valuable directions for the development of molecular diagnosis and biological therapies for mastitis.

Gene co-expression in response to Staphylococcus aureus infection reveals networks of genes with specific functions during bovine subclinical mastitis

Staphylococcus aureus is one of the most prevalent contagious bacterial pathogen of bovine mastitis. The subclinical mastitis it causes has long-term economic implications and it is difficult to control. To further understanding of the genetic basis of mammary gland defense against S. aureus infection, the transcriptomes of milk somatic cells from 15 cows with persistent natural S. aureus infection (S. aureus-positive, SAP) and 10 healthy control cows (HC) were studied by deep RNA-sequencing technology. Comparing the transcriptomes of SAP to HC group revealed 4,077 differentially expressed genes (DEG; 1,616 up- and 2,461 downregulated). Functional annotation indicated enrichment of DEG in 94 Gene Ontology (GO) and 47 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Terms related to the immune response and disease processes were mostly enriched for by upregulated DEG, whereas biological process terms related to cell adhesion, cell movement and localization, and tissue development were mostly enriched for by downregulated DEG. Weighted gene co-expression network analysis grouped DEG into 7 modules, the most important module (colored turquoise by software and here referred to as Turquoise module) was positively significantly correlated with S. aureus subclinical mastitis. The 1,546 genes in the Turquoise module were significantly enriched in 48 GO terms and 72 KEGG pathways, with 80% of them being disease- and immune-related terms [e.g., immune system process (GO:0002376), cytokine-cytokine receptor interaction (bta04060) and S. aureus infection (bta05150)]. Some DEG such as IFNG, IL18, IL1B, NFKB1, CXCL8, and IL12B were enriched in immune and disease pathways suggesting their possible involvement in the regulation of the host response to S. aureus infection. Four modules (Yellow, Brown, Blue, and Red) were negatively correlated (significantly) with S. aureus subclinical mastitis, and were enriched in functional annotations involved in the regulation of cell migration, cell communication, metabolic process, and blood circulatory system development, respectively. Application of sparse partial least squares discriminant analysis to genes of the Turquoise module identified 5 genes (NR2F6, PDLIM5, RAB11FIP5, ACOT4, and TMEM53) capable of explaining the majority of the differences in the expression patterns between SAP and HC cows. In conclusion, this study has furthered understanding of the genetic changes in the mammary gland and the molecular mechanisms underlying S. aureus mastitis, as well as revealed a list of candidate discriminant genes with potential regulatory roles in response to S. aureus infection.

Transcriptome-wide mapping of milk somatic cells upon subclinical mastitis infection in dairy cattle

BACKGROUND

Subclinical intramammary infection (IMI) represents a significant problem in maintaining dairy cows' health. Disease severity and extent depend on the interaction between the causative agent, environment, and host. To investigate the molecular mechanisms behind the host immune response, we used RNA-Seq for the milk somatic cells (SC) transcriptome profiling in healthy cows (n = 9), and cows naturally affected by subclinical IMI from Prototheca spp. (n = 11) and Streptococcus agalactiae (S. agalactiae; n = 11). Data Integration Analysis for Biomarker discovery using Latent Components (DIABLO) was used to integrate transcriptomic data and host phenotypic traits related to milk composition, SC composition, and udder health to identify hub variables for subclinical IMI detection.

RESULTS

A total of 1,682 and 2,427 differentially expressed genes (DEGs) were identified when comparing Prototheca spp. and S. agalactiae to healthy animals, respectively. Pathogen-specific pathway analyses evidenced that Prototheca's infection upregulated antigen processing and lymphocyte proliferation pathways while S. agalactiae induced a reduction of energy-related pathways like the tricarboxylic acid cycle, and carbohydrate and lipid metabolism. The integrative analysis of commonly shared DEGs between the two pathogens (n = 681) referred to the core-mastitis response genes, and phenotypic data evidenced a strong covariation between those genes and the flow cytometry immune cells (r² = 0.72), followed by the udder health (r² = 0.64) and milk quality parameters (r² = 0.64). Variables with r ≥ 0.90 were used to build a network in which the top 20 hub variables were identified with the Cytoscape cytohubba plug-in. The genes in common between DIABLO and cytohubba (n = 10) were submitted to a ROC analysis which showed they had excellent predictive performances in terms of discriminating healthy and mastitis-affected animals (sensitivity > 0.89, specificity > 0.81, accuracy > 0.87, and precision > 0.69). Among these genes, CIITA could play a key role in regulating the animals' response to subclinical IMI.

CONCLUSIONS

Despite some differences in the enriched pathways, the two mastitis-causing pathogens seemed to induce a shared host immune-transcriptomic response. The hub variables identified with the integrative approach might be included in screening and diagnostic tools for subclinical IMI detection.

Targeted Transcriptional Analysis of IgA Vasculitis, IgA Nephropathy, and IgA-Dominant Infection-Related Glomerulonephritis Reveals Both Distinct and Overlapping Immune Signatures

Key Points

Skin IL-9, calprotectin, and KIR gene expression may be predictive of subsequent kidney involvement in patients with IgAV. Histologically similar patients with IgAN, IgAV, and IgA-IRGN can be distinguished by their immune transcriptomes. Kidney biopsies from patients with IgA-IRGN are enriched for transcripts involved in granulocyte chemotaxis.

Background

IgA vasculitis (IgAV), IgA nephropathy (IgAN), and IgA-dominant infection-related glomerulonephritis (IgA-IRGN) have shared histopathologic features, but differences in clinical management and prognosis. The most serious IgAV organ involvement is in the kidneys (IgAV nephritis). In this study, we hypothesized that targeted immune transcript profiling could aid in (1 ) predicting the development of IgAV nephritis in patients with cutaneous IgAV and (2 ) differentiating IgAN, IgAV, and IgA-IRGN.

Methods

RNA was extracted from 24 formalin-fixed paraffin-embedded tissue specimens (16 kidney, 8 skin) from 21 patients with IgAV nephritis (n=7), IgAN (n=5), and IgA-IRGN (n=4), and IgAV skin biopsies from patients with (n=3) and without (n=5) IgAV nephritis. Differential gene expression and gene set enrichment analysis were performed on a total of 594 transcripts (Nanostring immunology panel) profiled using the nCounter system.

Results

Skin biopsies in patients with IgAV who develop kidney involvement exhibit reduced S100A8/S100A9 , IL9 , and killer cell immunoglobulin-like receptor expression. The kidney tissue immune transcriptomes of IgAN, IgAV, and IgA-IRGN are largely overlapping. IgA-IRGN kidney biopsies are, however, uniquely enriched for transcripts involved in granulocyte chemotaxis.

Conclusion

This study identifies immune transcript signatures that may predict IgAV nephritis in skin biopsies and distinguish IgA-IRGN from IgAN and IgAV in kidney biopsies.

Genome-Wide Signal Selection Analysis Revealing Genes Potentially Related to Sheep-Milk-Production Traits

Natural selection and domestication have shaped modern sheep populations into a vast range of phenotypically diverse breeds. Among these breeds, dairy sheep have a smaller population than meat sheep and wool sheep, and less research is performed on them, but the lactation mechanism in dairy sheep is critically important for improving animal-production methods. In this study, whole-genome sequences were generated from 10 sheep breeds, including 57 high-milk-yield sheep and 44 low-milk-yield sheep, to investigate the genetic signatures of milk production in dairy sheep, and 59,864,820 valid SNPs (Single Nucleotide Polymorphisms) were kept after quality control to perform population-genetic-structure analyses, gene-detection analyses, and gene-function-validation analyses. For the population-genetic-structure analyses, we carried out PCA (Principal Component Analysis), as well as neighbor-joining tree and structure analyses to classify different sheep populations. The sheep used in our study were well distributed in ten groups, with the high-milk-yield-group populations close to each other and the low-milk-yield-group populations showing similar classifications. To perform an exact signal-selection analysis, we used three different methods to find SNPs to perform gene-annotation analyses within the 995 common regions derived from the fixation index (F_ST), nucleotide diversity (Ɵπ), and heterozygosity rate (ZHp) results. In total, we found 553 genes that were located in these regions. These genes mainly participate in the protein-binding pathway and the nucleoplasm-interaction pathway, as revealed by the GO- and KEGG-function-enrichment analyses. After the gene selection and function analyses, we found that FCGR3A, CTSK, CTSS, ARNT, GHR, SLC29A4, ROR1, and TNRC18 were potentially related to sheep-milk-production traits. We chose the strongly selected genes, FCGR3A, CTSK, CTSS, and ARNT during the signal-selection analysis to perform a RT-qPCR (Reale time Quantitative Polymerase Chain Reaction) experiment to validate their expression-level relationship with milk production, and the results showed that FCGR3A has a significant negative relationship with sheep-milk production, while other three genes did not show any positive or negative relations. In this study, it was discovered and proven that the candidate gene FCGR3A potentially contributes to the milk production of dairy sheep and a basis was laid for the further study of the genetic mechanism underlying the strong milk-production traits of sheep.

Genetic polymorphisms in immune- and inflammation-associated genes and their association with bovine mastitis resistance/susceptibility

Bovine mastitis, the inflammation of the mammary gland, is a contagious disease characterized by chemical and physical changes in milk and pathological changes in udder tissues. Depressed immunity and higher expression of inflammatory cytokines with an elevated milk somatic cell count can be observed during mastitis in dairy cattle. The use of somatic cell count (SCC) and somatic cell score (SCS) as correlated traits in the indirect selection of animals against mastitis resistance is in progress globally. Traditional breeding for mastitis resistance seems difficult because of the low heritability (0.10-0.16) of SCC/SCS and clinical mastitis. Thus, genetic-marker-selective breeding to improve host genetics has attracted considerable attention worldwide. Moreover, genomic selection has been found to be an effective and fast method of screening for dairy cattle that are genetically resistant and susceptible to mastitis at a very early age. The current review discusses and summarizes the candidate gene approach using polymorphisms in immune- and inflammation-linked genes (CD4, CD14, CD46, TRAPPC9, JAK2, Tf, Lf, TLRs, CXCL8, CXCR1, CXCR2, C4A, C5, MASP2, MBL1, MBL2, LBP, NCF1, NCF4, MASP2, A2M, and CLU, etc.) and their related signaling pathways (Staphylococcus aureus infection signaling, Toll-like receptor signaling, NF-kappa B signaling pathway, Cytokine-cytokine receptor, and Complement and coagulation cascades, etc.) associated with mastitis resistance and susceptibility phenotypic traits (IL-6, interferon-gamma (IFN-γ), IL17, IL8, SCS, and SCC) in dairy cattle.

Milk yield variation partially attributed to blood oxygen-mediated neutrophil activation in lactating dairy goats

Blood oxygen is an essential component for numerous biological processes of mammalian animals. Milk production of ruminants largely relies on the supply of nutrients, such as glucose, amino acids and fatty acids. To define the regulatory role of blood oxygen availability in regard to milk production, seventy-five healthy Guanzhong dairy goats with similar body weight, days in milk and parities were selected. For each animal, milk yield was recorded and milk sample was collected to determine compositions. Milk vein blood was collected to determine parameters including blood gas, physio-biochemistry and haematology. Another blood sample was prepared for transcriptome and RT-qPCR. Results showed that both pressure of oxygen (pO₂) in the milk vein (positively) and numbers of neutrophils in mammary vein (negatively) were associated with milk yield of the animals. To learn the role of pO₂ in blood cell functionality, twelve animals (six with higher yield (H-group) and six with lower yield (L-group)) from seventy-five goats were selected. Compared with animals in L-group, goats in H-group were higher in pO₂ but lower in pCO₂, lactate, lactate dehydrogenase activity and neutrophil abundance in milk vein, compared with L-group. The blood transcriptome analysis suggested that compared with L-group, animals in H-group were depressed in functionality including neutrophil activation and metabolic pathways including glycolysis, NF-κB and HIF-1. Our result revealed that lower milk production could be associated with neutrophil activation responding to low pO₂ in the mammary vein. In the meantime, we highlighted the potential importance of blood oxygen as a milk yield regulator.

Genetic polymorphisms of TRAPPC9 and CD4 genes and their association with milk production and mastitis resistance phenotypic traits in Chinese Holstein

The present study was designed to evaluate the association of polymorphisms in bovine trafficking protein particle complex subunit 9 (TRAPPC9) and cluster of differentiation 4 (CD4) genes with milk production and mastitis resistance phenotypic traits in a different cattle population. Three single nucleotide polymorphisms (SNPs) (SNP1 Position: Chr14:2484891, SNP2 (rs110017379), SNP3 Position: Chr14:2525852) in bovine TRAPPC9 and one SNP (Position: Chr5:104010752) in CD4 were screened through Chinese Cow's SNPs Chip-I (CCSC-I) and genotyped in a population of 312 Chinese Holsteins (156: Mastitis, 156: Healthy). The results were analyzed using the general linear model in SAS 9.4. Our analysis revealed that milk protein percentage, somatic cell count (SCC), somatic cell score (SCS), serum cytokines interleukin 6 (IL-6) and interferon-gamma (IFN-γ) were significantly (P < 0.05) associated with at least one or more identified SNPs of TRAPPC9 and CD4 genes. Furthermore, the expression status of SNPs in CD4 and TRAPPC9 genes were verified through RT-qPCR. The expression analysis showed that genotypes GG in SNP3 of TRAPPC9 and TT genotype in SNP4 of CD4 showed higher expression level compared to other genotypes. The GG genotype in SNP2 and TT genotype in SNP3 of TRAPPC9 were associated with higher bovine milk SCC and lower IL6. Altogether, our findings suggested that the SNPs of TRAPPC9 and CD4 genes could be useful genetic markers in selection for milk protein improvement and mastitis resistance phenotypic traits in dairy cattle. The CCSC-I used in current study is proposed to be validate in different and large population of dairy cattle not only in China but also in other countries. Moreover, our analyses recommended that besides SCC and SCS, the association of genetic markers could also be considered with the serum cytokines (IL-6, IFN-γ) while selecting genetically mastitis resistance dairy cattle.

First report on delineation of differentially expressed genes and pathways in milk somatic cells of mastitic and healthy Murrah buffaloes

Despite immense contribution of buffaloes as dairy species, limited studies have addressed the bubaline mastitis as compared to cattle. This was the first differential transcriptomic study investigating the alterations induced by clinical mastitis in buffalo milk relative to healthy controls. Comparative gene expression profiling of three biological replicates of each group identified 1014 up-regulated and 999 down-regulated genes in the diseased buffaloes (Fold change > 2, FDR < 0.05). Activation of immune and inflammatory responses were the most enriched GO terms in the mastitic animals, with higher transcript abundance of many genes coding for anti-microbial proteins such as β-defensins, perforin, granzymes, granulysin, cathelicidins etc. Analysis of the gene regulatory interactions of the up-regulated DEGs identified many hub genes that govern the cellular and macromolecular metabolic processes (E2F4, E2F1, RBL2, FOXM1, IRF1 and MYB). This study contributes to an insightful understanding of molecular mechanisms governing immune response of buffaloes to mastitis.

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.

Preliminary study on gene regulation and its pathways in Chinese Holstein cows with clinical mastitis caused by Staphylococcus aureus

Introduction

Clinical mastitis (CM) is one of the most common diseases of dairy cows globally, has a complex aetiology and recurs easily. Staphylococcus aureus is a frequently isolated pathogen responsible for bovine mastitis and remains difficult to eradicate.

Material and Methods

To characterise the transcriptional profiles of dairy cows infected by S. aureus, we performed an RNA-seq analysis of peripheral blood leukocytes in lactating Chinese Holstein dairy cows with CM and did the same with healthy cows’ samples as controls.

Results

A total of 4,286 genes were detected in the CM cases infected with S. aureus which were differentially expressed compared to the controls, 3,085 of which were upregulated, the remainder being downregulated. Notably, we observed that some differentially expressed genes (DEGs) had strong protein–protein interaction. Of these, six downregulated DEGs (AKR1C4, PTGS2, HNMT, EPHX2, CMBL, and IDH1) were involved in the metabolic pathway, while eight upregulated DEGs (VWF, GP9, MYLK, GP6, F2RL3, ITGB3, GP5, and PRKG1) were associated with the platelet activation pathway.

Conclusion

The transcriptome dataset of CM cases would be a valuable resource for clinical guidance on anti-inflammatory medication and for deeper understanding of the biological processes of CM response to S. aureus infection, and it would enable us to identify specific genes for diagnostic markers and possibly for targeted therapy.

Weighted Single-Step Genome-Wide Association Study Uncovers Known and Novel Candidate Genomic Regions for Milk Production Traits and Somatic Cell Score in Valle del Belice Dairy Sheep

Simple Summary

Milk production is the most economically crucial dairy sheep trait and constitutes the major genetic enhancement purpose via selective breeding. Also, mastitis is one of the most frequently encountered diseases, having a significant impact on animal welfare, milk yield, and quality. The aim of this study was to identify genomic region(s) associated with the milk production traits and somatic cell score (SCS) in Valle del Belice sheep using single-step genome-wide association (ssGWA) and genotyping data from medium density SNP panels. We identified several genomic regions (OAR1, OAR2, OAR3, OAR4, OAR6, OAR9, and OAR25) and candidate genes implicated in milk production traits and SCS. Our findings offer new insights into the genetic basis of milk production traits and SCS in dairy sheep.

Abstract

The objective of this study was to uncover genomic regions explaining a substantial proportion of the genetic variance in milk production traits and somatic cell score in a Valle del Belice dairy sheep. Weighted single-step genome-wide association studies (WssGWAS) were conducted for milk yield (MY), fat yield (FY), fat percentage (FAT%), protein yield (PY), protein percentage (PROT%), and somatic cell score (SCS). In addition, our aim was also to identify candidate genes within genomic regions that explained the highest proportions of genetic variance. Overall, the full pedigree consists of 5534 animals, of which 1813 ewes had milk data (15,008 records), and 481 ewes were genotyped with a 50 K single nucleotide polymorphism (SNP) array. The effects of markers and the genomic estimated breeding values (GEBV) of the animals were obtained by five iterations of WssGBLUP. We considered the top 10 genomic regions in terms of their explained genomic variants as candidate window regions for each trait. The results showed that top ranked genomic windows (1 Mb windows) explained 3.49, 4.04, 5.37, 4.09, 3.80, and 5.24% of the genetic variances for MY, FY, FAT%, PY, PROT%, and total SCS, respectively. Among the candidate genes found, some known associations were confirmed, while several novel candidate genes were also revealed, including PPARGC1A, LYPLA1, LEP, and MYH9 for MY; CACNA1C, PTPN1, ROBO2, CHRM3, and ERCC6 for FY and FAT%; PCSK5 and ANGPT1 for PY and PROT%; and IL26, IFNG, PEX26, NEGR1, LAP3, and MED28 for SCS. These findings increase our understanding of the genetic architecture of six examined traits and provide guidance for subsequent genetic improvement through genome selection.

Protective Roles of Folic Acid in the Responses of Bovine Mammary Epithelial Cells to Different Virulent Staphylococcus aureus Strains

Mastitis caused by Staphylococcus aureus (S. aureus) infection is one of the most difficult diseases to treat in dairy cattle. Exploring the biological progression of S. aureus mastitis via the interaction between host, pathogen, and environment is the key to an effective and sustainable improvement of animal health. Here, two strains of S. aureus and a strain of MRSA (Methicillin-resistant Staphylococcus aureus) isolated from cows with different inflammation phenotypes were used to challenge Mac-T cells and to investigate their effects on the global transcriptome of the cells, then to explore the potential regulatory mechanisms of folic acid on S. aureus mastitis prevention. Differential gene expression or splicing analysis showed that different strains of S. aureus led to distinct transcriptional responses from the host immune system. Folic acid could protect host defense against the challenge of S. aureus and MRSA partially through activating cytoplasmic DNA sensing and tight junction pathway. ZBP1 at the upstream of cytoplasmic DNA sensing pathway was verified and related to anti-pathogen through RNA interference. Further enrichment analysis using these transcriptome data with cattle large-scale genome-wide association study (GWAS) data confirmed that ZBP1 gene is highly associated with bovine somatic cell score (SCS) trait. Our data shed light on the potential effect of FA through regulating key gene and then protect host cells' defense against S. aureus and MRSA.

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.

The Recombinant Expression Proteins FnBP and ClfA From Staphylococcus aureus in Addition to GapC and Sip From Streptococcus agalactiae Can Protect BALB/c Mice From Bacterial Infection

Dairy cow mastitis is a serious disease that is mainly caused by intramammary infection with Staphylococcus aureus and Streptococcus agalactiae [group B streptococcus (GBS)]. FnBP and ClfA are the virulence factors of S. aureus, while GapC is the respective factor for S. agalactiae. Sip is a highly immunogenic protein, and it is conserved in all GBS serotypes. In this study, we analyzed the abovementioned four genes prepared a FnBP+ClfA chimeric protein (FC), a GapC+Sip chimeric protein (GS), and a FnBP+ClfA+GapC+Sip chimeric protein (FCGS) based on the antigenic sites to evaluate their use in vaccine development. After expression and purification of the recombinant proteins in Escherichia coli, BALB/c mice were immunized with them to examine resistance effects. The total lethal and half lethal doses of S. aureus and S. agalactiae were then measured, and the immunoprotective effects of the fusion proteins were evaluated. The FC and FCGS chimeric proteins could induce mice to produce high levels of antibodies, and bacterial loads were significantly reduced in the spleens and livers after challenge. After immunization with FCGS, the recipients resisted the attacks of both S. aureus and S. agalactiae, indicating the potential of the fusion protein as a mastitis vaccine.

Global transcriptomic profiles of circulating leucocytes in early lactation cows with clinical or subclinical mastitis

Bovine mastitis, an inflammatory disease of the mammary gland, is classified as subclinical or clinical. Circulating neutrophils are recruited to the udder to combat infection. We compared the transcriptomic profiles in circulating leukocytes between healthy cows and those with naturally occurring subclinical or clinical mastitis. Holstein Friesian dairy cows from six farms in EU countries were recruited. Based on milk somatic cell count and clinical records, cows were classified as healthy (n = 147), subclinically (n = 45) or clinically mastitic (n = 22). Circulating leukocyte RNA was sequenced with Illumina NextSeq single end reads (30 M). Differentially expressed genes (DEGs) between the groups were identified using CLC Genomics Workbench V21, followed by GO enrichment analysis. Both subclinical and clinical mastitis caused significant changes in the leukocyte transcriptome, with more intensive changes attributed to clinical mastitis. We detected 769 DEGs between clinical and healthy groups, 258 DEGs between subclinical and healthy groups and 193 DEGs between clinical and subclinical groups. Most DEGs were associated with cell killing and immune processes. Many upregulated DEGs in clinical mastitis encoded antimicrobial peptides (AZU1, BCL3, CAMP, CATHL1, CATHL2, CATHL4,CATHL5, CATHL6, CCL1, CXCL2, CXCL13, DEFB1, DEFB10, DEFB4A, DEFB7, LCN2, PGLYRP1, PRTN3, PTX3, S100A8, S100A9, S100A12, SLC11A1, TF and LTF) which were not upregulated in subclinical mastitis. The use of transcriptomic profiles has identified a much greater up-regulation of genes encoding antimicrobial peptides in circulating leukocytes of cows with naturally occurring clinical compared with subclinical mastitis. These could play a key role in combatting disease organisms.

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

Staphylococcus aureus- induced mastitis is one of the most intractable problems for the dairy industry, which causes loss of milk yield and early slaughter of cows worldwide. Few studies have used a comprehensive approach based on the integrative analysis of miRNA and mRNA expression profiles to explore molecular mechanism in bovine mastitis caused by S. aureus. In this study, S. aureus (A1, B1 and C1) and sterile phosphate buffered saline (PBS) (A2, B2 and C2) were introduced to different udder quarters of three individual cows, and transcriptome sequencing and microarrays were utilized to detected miRNA and gene expression in mammary glands from the challenged and control groups. A total of 77 differentially expressed microRNAs (DE miRNAs) and 1625 differentially expressed genes (DEGs) were identified. Gene Ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed that multiple DEGs were enriched in significant terms and pathways associated with immunity and inflammation. Integrative analysis between DE miRNAs and DEGs proved that miR-664b, miR-23b-3p, miR-331-5p, miR-19b and miR-2431-3p were potential factors regulating the expression levels of CD14 Molecule (CD14), G protein subunit gamma 2 (GNG2), interleukin 17A (IL17A), collagen type IV alpha 1 chain (COL4A1), microtubule associated protein RP/EB family member 2 (MAPRE2), member of RAS oncogene family (RAP1B), LDOC1 regulator of NFKB signaling (LDOC1), low-density lipoprotein receptor (LDLR) and S100 calcium binding protein A9 (S100A9) in bovine mastitis caused by S. aureus. These findings could enhance the understanding of the underlying immune response in bovine mammary glands against S. aureus infection and provide a useful foundation for future application of the miRNA–mRNA-based genetic regulatory network in the breeding cows resistant to S. aureus.

Transcriptome analysis reveals potential immune function-related regulatory genes/pathways of female Lubo goat submandibular glands at different developmental stages

Background

The submandibular glands, as major salivary glands, participate in rumen digestion in goats. Sialic acid, lysozyme, immunoglobulin A (IgA), lactoferrin and other biologically active substances secreted in the submandibular glands were reported in succession, which suggests that the submandibular gland may have immune functions in addition to participating in digestion. The aim of this study was to map the expression profile of differentially expressed genes (DEGs) at three different stages by transcriptome sequencing, screen immune-related genes and pathways by bioinformatics methods, and predict the immune function of submandibular glands at different developmental stages.

Methods

Nine submandibular gland tissue samples were collected from groups of 1-month-old kids, 12-month-old adolescent goats and 24-month-old adult goats (3 samples from each group), and high-throughput transcriptome sequencing was conducted on these samples. The DEGs among the three stages were screened and analysed. Key genes and signalling pathways were selected via protein-protein interaction (PPI) network analysis.

Results

The results revealed 2,706, 2,525 and 52 DEGs between 1-month-old and 12-month-old goats, between 1-month-old and 24-month-old goats, and between 12-month-old and 24-month-old goats, respectively. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses indicated that most of the DEGs were enriched in immune- related GO terms and pathways. Based on functional enrichment analysis and network analysis, 10 genes (PTPRC, CD28, SELL, LCP2, MYC, LCK, ZAP70, ITGB2, SYK and CCR7), two signalling pathways (the T cell receptor signalling pathway and the NF-κβ signalling pathway) and eight GO terms (T cell receptor signalling pathway, neutrophil mediated immunity, B cell mediated immunity, regulation of alpha-beta T cell activation, positive regulation of T cell proliferation, regulation of leukocyte differentiation, positive regulation of antigen receptor-mediated signalling pathway, positive regulation of lymphocyte proliferation) that may play key roles in the immune functions of the goat submandibular glands at different developmental stages were identified. Moreover, we found that eight antibacterial peptide-encoding genes were downregulated in the tuberculosis and salivary secretion pathways, while all immunoglobulins were upregulated in 10 immune system pathways. These findings indicate that the submandibular glands may be important immunological organs during the growth process of goats and that the immune function of these glands gradually weakens with age up to 12 months but remains relatively stable after 12 months of age. Overall, this study will improve our understanding of transcriptional regulation related to goat submandibular gland immune function.