The genetic background of clinical mastitis in Holstein-Friesian cattle

Copy number variant scan in more than four thousand Holstein cows bred in Lombardy, Italy

Copy Number Variants (CNV) are modifications affecting the genome sequence of DNA, for instance, they can be duplications or deletions of a considerable number of base pairs (i.e., greater than 1000 bp and up to millions of bp). Their impact on the variation of the phenotypic traits has been widely demonstrated. In addition, CNVs are a class of markers useful to identify the genetic biodiversity among populations related to adaptation to the environment. The aim of this study was to detect CNVs in more than four thousand Holstein cows, using information derived by a genotyping done with the GGP (GeneSeek Genomic Profiler) bovine 100K SNP chip. To detect CNV the SVS 8.9 software was used, then CNV regions (CNVRs) were detected. A total of 123,814 CNVs (4,150 non redundant) were called and aggregated into 1,397 CNVRs. The PCA results obtained using the CNVs information, showed that there is some variability among animals. For many genes annotated within the CNVRs, the role in immune response is well known, as well as their association with important and economic traits object of selection in Holstein, such as milk production and quality, udder conformation and body morphology. Comparison with reference revealed unique CNVRs of the Holstein breed, and others in common with Jersey and Brown. The information regarding CNVs represents a valuable resource to understand how this class of markers may improve the accuracy in prediction of genomic value, nowadays solely based on SNPs markers.

Genome-wide CNV analysis reveals candidate genes associated with milk production traits in water buffalo (Bubalus bubalis)

Buffaloes are vital contributors to the global dairy industry. Understanding the genetic basis of milk production traits in buffalo populations is essential for breeding programs and improving productivity. In this study, we conducted whole-genome resequencing on 387 buffalo genomes from 29 diverse Asian breeds, including 132 river buffaloes, 129 swamp buffaloes, and 126 crossbred buffaloes. We identified 36,548 copy number variant (CNVs) spanning 133.29 Mb of the buffalo genome, resulting in 2,100 copy number variant regions (CNVRs), with 1,993 shared CNVRs being found within the studied buffalo types. Analyzing CNVRs highlighted distinct genetic differentiation between river and swamp buffalo subspecies, verified by evolutionary tree and principal component analyses. Admixture analysis grouped buffaloes into river and swamp categories, with crossbred buffaloes displaying mixed ancestry. To identify candidate genes associated with milk production traits, we employed 3 approaches. First, we used Vst-based population differentiation, revealing 11 genes within CNVRs that exhibited significant divergence between different buffalo breeds, including genes linked to milk production traits. Second, expression quantitative loci (eQTL) analysis revealed differential expression of CNVR-driven genes (DECGs) associated with milk production traits. Notably, known milk production-related genes were among these DECGs, validating their relevance. Last, a genome-wide association study (GWAS) identified 3 CNVRs significantly linked to peak milk yield. Our study provides comprehensive genomic insights into buffalo populations and identifies candidate genes associated with milk production traits. These findings facilitate genetic breeding programs aimed at increasing milk yield and improving quality in this economically important livestock species.

Genetic parameters for udder conformation traits derived from Cartesian coordinates generated by robotic milking systems in North American Holstein cattle

Udder conformation is directly related to milk yield, cow health, workability, and welfare. Automatic milking systems (AMS, also known as milking robots) have become popular worldwide, and the number of dairy farms adopting these systems have increased considerably over the past years. In each milking visit, AMS record the location of the 4 teats as Cartesian coordinates in a xyz plan, which can then be used to derive udder conformation traits. AMS generate a large amount of per milking visit data for individual cows, which contribute to an accurate assessment of important traits such as udder conformation without the addition of human classifier errors (in subjective scoring systems). Therefore, the primary objectives of this study were to estimate genomic-based genetic parameters for udder conformation traits derived from AMS records in North American Holstein cattle and to assess the genetic correlation between the derived traits for evaluating the feasibility of multi-trait genomic selection for breeding cows that are more suitable for milking in AMS. The Cartesian teat coordinates measured during each milking visit were collected by 36 milking robots in 4,480 Holstein cows from 2017 to 2021, resulting in 5,317,488 records. A total of 4,118 of these Holstein cows were also genotyped for 57,600 single nucleotide polymorphisms. Five udder conformation traits were derived: udder balance (UB, mm), udder depth (UD, mm), front teat distance (FTD, mm), rear teat distance (RTD, mm), and distance front-rear (DFR, mm). In addition, 2 traits directly related to cow productivity in the system were added to the study: daily milk yield (DY) and milk electroconductivity (EC; as an indicator of mastitis). Variance components and genetic parameters for UB, UD, FTD, RTD, DFR, DY, and EC were estimated based on repeatability animal models. The estimates of heritability (±standard error, SE) for UB, UD, FTD, RTD, DFR, DY, and EC were 0.41 ± 0.02, 0.79 ± 0.01, 0.53 ± 0.02, 0.40 ± 0.02, 0.65 ± 0.02, 0.20 ± 0.02, and 0.46 ± 0.02, respectively. The repeatability estimates (±SE) for UB, UD, FTD, RTD, and DFR were 0.82 ± 0.01, 0.93 ± 0.01, 0.87 ± 0.01, 0.83 ± 0.01, and 0.88 ± 0.01, respectively. The strongest genetic correlations were observed between the FTD and RTD (0.54 ± 0.03), UD and DFR (-0.47 ± 0.03), DFR and FTD (0.32 ± 0.03), and UD and FTD (-0.31 ± 0.03). These results suggest that udder conformation traits derived from Cartesian coordinates from AMS are moderately to highly heritable. Furthermore, the moderate genetic correlations between these traits should be considered when developing selection sub-indexes. The most relevant genetic correlations between traits related to cow milk productivity and udder conformation traits were between UD and EC (-0.25 ± 0.03) and between DFR and DY (0.30 ± 0.04), in which both genetic correlations are favorable. These findings will contribute to the design of genomic selection schemes for improving udder conformation in North American Holstein cattle, especially in precision dairy farms.

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.

Breed and non-genetic risk factors associated with the prevalence of subclinical mastitis in livestock systems of Arauca, Colombian orinoquia

Subclinical Mastitis (SCM) is caused by several factors associated with the interaction of the individual cow, cow management, and the environment. The aim of this study was to assess the breed and non-genetic risk factors on the prevalence of SCM in cows from the dual-purpose livestock system on the floodplain of Arauca, Colombian Orinoquia. Milk samples were taken from the individual mammary quarters of 481 cows representing 28 different farms where the electrical conductivity (EC) test applied. To determine the factors associated with SCM, a multiple logistic regression analysis was used. The response variable was the SCM presence (1), or absence (0) obtained with the EC test. Breed was included as a genetic risk factor, and as non- genetic risk factors: number of cows in production, daily milk production, lactation month, cow age, climatic period, body condition, and calving number. The factors that were significantly associated with the SCM presence were body condition, climatic period, and breed (p < 0.05). Odds ratio (OR) analysis of significant effects indicates that for each unit increase in body condition, the OR of having animals with SCM is reduced by 71%. In the dry period the OR of animals with SCM increases by 150% compared to the rainy period. Composite breeds reduce the OR of SCM animals by 73%, compared with Indicus-predominance animals. In this study, the occurrence of SCM in extensive management systems in Arauca, Colombian Orinoquia, is determined by the risk factors of breed, climatic period, and body condition.

A Genome-Wide Association Study for Resistance to Tropical Theileriosis in Two Bovine Portuguese Autochthonous Breeds

The control of Tropical Theileriosis, a tick-borne disease with a strong impact on cattle breeding, can be facilitated using marker-assisted selection in breeding programs. Genome-wide association studies (GWAS) using high-density arrays are extremely important for the ongoing process of identifying genomic variants associated with resistance to Theileria annulata infection. In this work, single-nucleotide polymorphisms (SNPs) were analyzed in the Portuguese autochthonous cattle breeds Alentejana and Mertolenga. In total, 24 SNPs suggestive of significance (p ≤ 10-4) were identified for Alentejana cattle and 20 SNPs were identified for Mertolenga cattle. The genomic regions around these SNPs were further investigated for annotated genes and quantitative trait loci (QTLs) previously described by other authors. Regarding the Alentejana breed, the MAP3K1, CMTM7, SSFA2, and ATG13 genes are located near suggestive SNPs and appear as candidate genes for resistance to Tropical Theileriosis, considering its action in the immune response and resistance to other diseases. On the other hand, in the Mertolenga breed, the UOX gene is also a candidate gene due to its apparent link to the pathogenesis of the disease. These results may represent a first step toward the possibility of including genetic markers for resistance to Tropical Theileriosis in current breed selection programs.

Environmental and Breed Risk Factors Associated with the Prevalence of Subclinical Mastitis in Dual-Purpose Livestock Systems in the Arauca Floodplain Savannah, Colombian Orinoquia

The aim of this study was to assess the environmental and breed risk factors associated with the prevalence of subclinical mastitis (SCM) in cows in the dual-purpose livestock system of Arauca, Colombian Orinoquia. Milk samples were taken from 1924 mammary quarters, corresponding to 481 cows on 28 different farms, and the California Mastitis Test (CMT) was applied. Risk factors associated with SCM were determined using multiple logistic regression analysis. The response variable was the presence (1) or absence (0) of SCM. Breed was included as a genetic risk factor, and daily milk production, number of cows in production, lactation month, calving number, cow age, climatic period, and body condition were included as environmental risk factors. The analysis of the odds ratio (OR) of significant effects indicated that the factors significantly associated with the presence of SCM were the number of cows (OR = 2.29; p = 0.005), milk production (OR = 0.88; p = 0.045), and the Taurus-Indicus breeds (OR = 1.79; p = 0.009) and composite breed (OR = 3.95; p = 0.005). In this study, the occurrence of SCM was determined by the following risk factors: number of cows, milk production, and breed. Likewise, the highest prevalence seemed to occur on farms with less technological development and sanitary management of producers from the lowest socioeconomic stratum.

Deubiquitinase OTUD3: a double-edged sword in immunity and disease

Deubiquitination is an important form of post-translational modification that regulates protein homeostasis. Ovarian tumor domain-containing proteins (OTUDs) subfamily member OTUD3 was identified as a deubiquitinating enzyme involved in the regulation of various physiological processes such as immunity and inflammation. Disturbances in these physiological processes trigger diseases in humans and animals, such as cancer, neurodegenerative diseases, diabetes, mastitis, etc. OTUD3 is aberrantly expressed in tumors and is a double-edged sword, exerting tumor-promoting or anti-tumor effects in different types of tumors affecting cancer cell proliferation, metastasis, and metabolism. OTUD3 is regulated at the transcriptional level by a number of MicroRNAs, such as miR-520h, miR-32, and miR101-3p. In addition, OTUD3 is regulated by a number of post-translational modifications, such as acetylation and ubiquitination. Therefore, understanding the regulatory mechanisms of OTUD3 expression can help provide insight into its function in human immunity and disease, offering the possibility of its use as a therapeutic target to diagnose or treat disease.

Genetical analysis of mastitis and reproductive traits in first-parity Holstein cows using standard and structural equation modelling

The present study aimed to investigate the causal relationships between clinical mastitis and some reproductive traits, including success at first insemination (SFI), the number of inseminations to pregnancy (INS), the interval from calving to first service (CTFS), first and last service interval (IFL), and open days (OD) in first-parity Holstein cows. For this purpose, the records of 58 281 first parity Holstein cows were analysed. These data sets were collected from 17 large dairy herds from 2008 to 2017. Recursive Mixed Models (RMMs) were applied and compared with the estimations under Standard Mixed Models. Then, one trivariate and three bivariate Gaussian-threshold models were used for the analyses. Recursive models were applied, considering that clinical mastitis can influence fertility traits. Mastitis is considered a covariate for the reproductive traits to determine their causal relationship. The results of this study indicated that causal effects of mastitis on SFI (on the observed scale, %), CTFS, IFL, OD, and INS were -5.7%, 3.3 days, 12.27 days, seven days, and 0.26 services, respectively. The estimated structural coefficients of the recursive models in the first parity imply that mastitis significantly lengthened the fertility interval and decreased the conception rate. In addition, genetic, residual, and phenotypic correlations between mastitis and the reproductive traits under both models were statistically significant. Results of genetic correlations between mastitis and fertility traits suggest that more incidence of mastitis during lactation is related to the delays in the heat show and pregnancy rate after insemination. In summary, considering the causal effects under RMMs may be advantageous to comprehend complicated relationships between complex traits better.

Effectiveness of Intramammary Antibiotics, Internal Teat Sealants, or Both at Dry-Off in Dairy Cows: Clinical Mastitis and Culling Outcomes

Intramammary antibiotic (AB) and internal teat sealants (TS) infusion at dry-off have been used to prevent intramammary infections (IMI) in dairy cows during the dry period and reduce the risk of mastitis during the dry period and subsequent lactation. A randomized clinal trial was completed on eight California dairy herds to estimate the effects of different dry cow therapies (AB, TS, AB + TS or None) on clinical mastitis and culling. A total of 1273 cows were randomized to one of the four treatment groups over summer and winter seasons. For each enrolled cow, microbiological testing was done on quarter milk samples collected from the first detection of clinical mastitis within the first 150 days in milk (DIM) in the subsequent lactation. Statistical analysis was done using generalized linear mixed models. There were no significant differences in the odds of clinical mastitis or culling between cows treated with AB, TS, or AB + TS compared to the controls. Dry cow therapy with AB and/or TS had no statistically significant effect on clinical mastitis and cow culling during the first 150 DIM.

Effect of Two Different Drug-Resistant Staphylococcus aureus Strains on the Physiological Properties of MAC-T Cells and Their Transcriptome Analysis

Staphylococcus aureus (S. aureus) is one of the main pathogens causing mastitis in dairy cows. The current work mainly focuses on the pathway of apoptosis induction in MAC-T cells caused by S. aureus infection or other factors. However, the physiological characteristics of S. aureus infected MAC-T cells and the resulting mRNA expression profile remain unknown particularly in the case of diverse drug resistant strains. Methicillin-resistant S. aureus (MRSA) and methicillin-susceptible S. aureus (MSSA) strains were used to infect MAC-T cells to investigate this issue. The adhesion, invasion and apoptosis ability of MRSA-infected group and MSSA-infected group was assessed over time (2, 4, 6, 8, and 12 h). After 8 h, the RNA sequencing was conducted on the MRSA-infected and the MSSA-infected with uninfected MAC-T cells as controls. The results showed that the adhesion and invasion ability of MRSA-infected and MSSA-infected to MAC-T cells increased and then decreased with infection time, peaking at 8 h. The adhesion and invasion rates of the MSSA-infected were substantially lower than those of the MRSA-infected, and the invasion rate of the MSSA-infected group was nearly non-existent. Then the apoptosis rate of MAC-T cells increased as the infection time increased. The transcriptome analysis revealed 549 differentially expressed mRNAs and 390 differentially expressed mRNAs in MRSA-infected and MSSA-infected MAC-T cells, respectively, compared to the uninfected MAC-T cells. According to GO analysis, these differentially expressed genes were involved in immune response, inflammation, apoptosis, and other processes. The Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis indicated the following pathways were linked to adhesion, invasion inflammation and apoptosis, including AMPK, FOXO, HIF-1, IL-17, JAK-STAT, MAPK, mTOR, NF-κB, p53, PI3K-Akt, TNF, Toll-like receptor, Rap1, RAS, prion disease, the bacterial invasion of epithelial cells pathway. We found 86 DEGs from 41 KEGG-enriched pathways associated with adhesion, invasion, apoptosis, and inflammation, all of which were implicated in MAC-T cells resistance to MRSA and MSSA infection. This study offers helpful data toward understanding the effect of different drug-resistant S. aureus on dairy cow mammary epithelial cells and aid in the prevention of mastitis in the dairy industry.

Analysis of mRNA and Long Non-Coding RNA Expression Profiles in Developing Yorkshire Pig Spleens

Simple Summary

Epidemic disease is a prominent problem in intensive pig production. The spleen is a blood bank and the largest immune organ, and most of the diseases in pig farms will be reflected as spleen abnormality. The results showed how the mRNA–lncRNA expression profiles in Yorkshire spleens varied with age (seven, 90, and 180 days after birth). Our study shows that 90 days after birth the gene-expression profile of pig spleen no longer changes significantly. The results are helpful for a better understanding of the transcriptome and functional genomics of spleen tissue in farm animals and could provide reference for precision pig disease research and prevention and control in pig farms.

Abstract

Epidemic diseases cause great economic loss in pig farms each year; some of these diseases are characterized mainly in the spleen, but mRNA and lncRNA (long non-coding RNA) expression networks in developing Yorkshire pig spleens remain obscure. Here, we profiled the systematic characters of mRNA and lncRNA repertoires in three groups of spleens from nine Yorkshire pigs, each three aged at seven days, 90 days, and 180 days. By using a precise mRNA and lncRNA identification pipeline, we identified 19,647 genes and 219 known and 3219 putative lncRNA transcripts; 1729 genes and 64 lncRNAs therein were found to express differentially. The gene expression characteristics of genes and lncRNAs were found to be basically fixed before 90 days after birth. Three large gene expression modules were detected. The enrichment analyses of differentially expressed genes and the potential target genes of differentially expressed lncRNAs both displayed the crucial roles of up-regulation in immune activation and hematopoiesis, and down-regulation in cell replication and division in 90 days and 180 days compared to seven days. ENSSSCT00000001325 was the only lncRNA transcript that existed in the three groups. CDK1, PCNA, and PLK were detected to be node genes that varied with age. This study contributes to a further understanding of mRNA and lncRNA expression in different developmental pig spleens.

Incidence, Etiology, and Risk Factors of Clinical Mastitis in Dairy Cows under Semi-Tropical Circumstances in Chattogram, Bangladesh

Simple Summary

Bovine clinical mastitis is an inflammatory disease of the mammary gland associated with visual changes in the milk and/or the udder. We show that the incidence of clinical mastitis in commercial dairy farms in Bangladesh is high but with large variation between farms. Streptococci and non-aureus Staphylococci were the most frequently isolated bacteria from quarter milk samples. Staphylococcus aureus and non-aureus Staphylococci were often resistant against penicillin and oxacillin. This work suggests an urgent need for improved udder health management and specifically a more prudent use of antimicrobial agents following a treatment protocol.

Abstract

Clinical mastitis (CM) is an important production disease in dairy cows, but much of the knowledge required to effectively control CM is lacking, specifically in low-income countries where most farms are small and have specific dairy management, such as regular udder cleaning and practicing hand milking. Therefore, we conducted a 6-month-long cohort study to (a) estimate the incidence rate of clinical mastitis (IRCM) at the cow and quarter level, (b) identify risk factors for the occurrence of CM, (c) describe the etiology of CM, and (d) quantify antimicrobial susceptibility (AMS) against commonly used antimicrobial agents in S. aureus and non-aureus Staphylococcus spp. (NAS) in dairy farms in the Chattogram region of Bangladesh. On 24 farms, all cows were monitored for CM during a 6-month period. Cases of CM were identified by trained farmers and milk samples were collected aseptically before administering any antimicrobial therapy. In total, 1383 lactating cows were enrolled, which totaled 446 cow-years at risk. During the study period, 196 new cases of CM occurred, resulting in an estimated crude IRCM of 43.9 cases per 100 cow-years, though this varied substantially between farms. Among the tested CM quarter samples, Streptococci (22.9%) followed by non-aureus staphylococci (20.3%) were the most frequently isolated pathogens and resistance of S. aureus and NAS against penicillin (2 out of 3 and 27 out of 39 isolates, respectively) and oxacillin (2 out of 3 and 38 out of 39 isolates, respectively) was common. The IRCM was associated with a high milk yield, 28 to 90 days in milk, and a higher body condition score. Our results show that there is substantial room for udder health improvement on most farms.

Role of the JAK-STAT Pathway in Bovine Mastitis and Milk Production

Simple Summary

The cytokine-activated Janus kinase (JAK)—signal transducer and activator of transcription (STAT) pathway has an important role in the regulation of immunity and inflammation. In addition, the signaling of this pathway has been reported to be associated with mammary gland development and milk production. Because of such important functions, the JAK-STAT pathway has been widely targeted in both human and animal diseases as a therapeutic agent. Recently, the JAK2, STATs, and inhibitors of the JAK-STAT pathway, especially cytokine signaling suppressors (SOCSs), have been reported to be associated with milk production and mastitis-resistance phenotypic traits in dairy cattle. Thus, in the current review, we attempt to overview the development of the JAK-STAT pathway role in bovine mastitis and milk production.

Abstract

The cytokine-activated Janus kinase (JAK)—signal transducer and activator of transcription (STAT) pathway is a sequence of communications between proteins in a cell, and it is associated with various processes such as cell division, apoptosis, mammary gland development, lactation, anti-inflammation, and immunity. The pathway is involved in transferring information from receptors on the cell surface to the cell nucleus, resulting in the regulation of genes through transcription. The Janus kinase 2 (JAK2), signal transducer and activator of transcription A and B (STAT5 A & B), STAT1, and cytokine signaling suppressor 3 (SOCS3) are the key members of the JAK-STAT pathway. Interestingly, prolactin (Prl) also uses the JAK-STAT pathway to regulate milk production traits in dairy cattle. The activation of JAK2 and STATs genes has a critical role in milk production and mastitis resistance. The upregulation of SOCS3 in bovine mammary epithelial cells inhibits the activation of JAK2 and STATs genes, which promotes mastitis development and reduces the lactational performance of dairy cattle. In the current review, we highlight the recent development in the knowledge of JAK-STAT, which will enhance our ability to devise therapeutic strategies for bovine mastitis control. Furthermore, the review also explores the role of the JAK-STAT pathway in the regulation of milk production in dairy cattle.

Prospecting genomic regions associated with milk production traits in Egyptian buffalo

The objectives of the current study were to detect putative genomic loci and to identify candidate genes associated with milk production traits in Egyptian buffalo. A total number of 161 479 daily milk yield (DMY) records and 60 318 monthly measures for fat and protein percentages (FP and PP, respectively), along with fat and protein yields (FY and PY, respectively) from 1670 animals were used. Genotyping was performed using Axiom® Buffalo Genotyping 90 K array. Genome-wide association study (GWAS) for each trait was performed using PLINK. After Bonferroni correction, 47 SNPs were associated with one or more milk production traits. These SNPs were distributed over 36 quantitative trait loci (QTL) and located on 20 buffalo chromosomes (BBU). For the 47 SNPs, one was overlapped for three traits (DMY, FY, and PY), six were associated with two traits (one for PP and PY and five for FY and PY) while the rest were associated with only one trait. Out of 36 identified QTL, eleven were overlapped with previously reported loci in buffalo and/or cattle populations. Some of these SNPs are placed within or close to potential candidate genes, for example: TPD52, ZBTB10, RALYL and SNX16 on BBU15, ADGRD1 on BBU17, ESRRG on BBU5 and GRIP1 on BBU4. This is the first reported study between genome-wide markers and milk components in Egyptian buffalo. Our findings provide useful information to explore the genetic mechanisms and relevant genes contributing to the variation in milk production traits. Further confirmation studies with larger population size are necessary to validate the findings and detect the causal genetic variants.

Effects of Peptidoglycan, Lipoteichoic Acid and Lipopolysaccharide on Inflammation, Proliferation and Milk Fat Synthesis in Bovine Mammary Epithelial Cells

The mammary gland of the cow is particularly susceptible to infections of a wide range of pathogenic bacteria, including both Gram-positive and Gram-negative bacteria. The endotoxins of these pathogenic bacteria include peptidoglycan (PGN), lipoteichoic acid (LTA) and lipopolysaccharide (LPS), and they are the pathogen-associated molecular patterns (PAMPs) to induce mastitis. LPS can directly inhibit proliferation and milk fat synthesis of bovine mammary epithelial cells (BMECs) while inducing mastitis, but it is unclear whether PGN and LTA also have such effects. Furthermore, since the three PAMPs usually appear simultaneously in the udder of cows with mastitis, their synergistic effects on proliferation and milk fat synthesis of BMECs are worth investigating. The immortalized BMECs (MAC-T cells) were stimulated for 24 h using various concentrations of PGN, LTA and LPS, respectively, to determine the doses that could effectively cause inflammatory responses. Next, the cells were stimulated for 24 h with no endotoxins (CON), PGN, LTA, LPS, PGN + LTA, and PGN + LTA + LPS, respectively, with the predetermined doses to analyze their effects on proliferation and milk fat synthesis of BMECs. PGN, LTA and LPS successfully induced inflammatory responses of BMECs with doses of 30, 30 and 0.1 μg/mL, respectively. Although the proliferation of BMECs was significantly inhibited in the following order: LTA < PGN + LTA < PGN + LTA + LPS, there was no change in cell morphology and cell death. LTA significantly promoted the expression of fatty acid synthesis-related genes but did not change the content of intracellular triglyceride (TG), compared with the CON group. The mRNA expression of fatty acid synthesis-related genes in the LPS group was the lowest among all the groups. Meanwhile, LPS significantly decreased the content of intracellular non-esterified fatty acids (NEFAs) and TG, compared with the CON group. PGN had no effects on milk fat synthesis. Co-stimulation with PGN, LTA and LPS significantly increased the expression of fat acid synthesis-related genes and the intracellular NEFAs, but decreased intracellular TG, compared with sole LPS stimulation. Collectively, PGN, LTA and LPS showed an additive effect on inhibiting proliferation of BMECs. The promoting role of LTA in fatty acid synthesis might offset the negative effects of LPS in this regard, but co-stimulation with PGN, LTA and LPS significantly decreased intracellular TG content.

The Synergism of PGN, LTA and LPS in Inducing Transcriptome Changes, Inflammatory Responses and a Decrease in Lactation as Well as the Associated Epigenetic Mechanisms in Bovine Mammary Epithelial Cells

Mastitis is usually caused by a variety of pathogenic bacteria that include both Gram-positive and Gram-negative bacteria. Lipopolysaccharide (LPS) is the pathogen-associated molecular pattern (PAMP) of Gram-negative bacteria, and peptidoglycan (PGN) and lipoteichoic acid (LTA) are those of Gram-positive bacteria. The effects of LPS, PGN and/or LTA on inflammatory response and lactation in bovine mammary epithelial cells (BMECs) are well studied, but the epigenetic mechanisms of their effects received less attention. Furthermore, since the three PAMPs are often simultaneously present in the udder of cows with mastitis, it has implications in practice to study their additive effects. The results show that co-stimulation of bovine mammary epithelial cells with PGN, LTA, and LPS induced a higher number of differentially expressed genes (DEGs) and greater expressions of inflammatory factors including interleukin (IL)-1β, IL-6, IL-8, tumor necrosis factor-α (TNF-α), chemokine (C-X-C motif) ligand (CXCL)1, and CXCL6. In addition, co-stimulation further increased DNA hypomethylation compared with sole LPS stimulation. Co-stimulation greatly decreased casein expression but did not further decrease histone acetylation levels and affect the activity of histone acetyltransferase (HAT) and histone deacetylase (HDAC), compared with sole LPS stimulation. Collectively, this study demonstrated that PGN, LTA, and LPS had an additive effect on inducing transcriptome changes and inflammatory responses in BMECs, probably through inducing a greater decrease in DNA methylation. Co-stimulation with PGN, LTA, and LPS decreased casein expression to a greater degree, but it might not be linked to histone acetylation and HAT and HDAC activity.

PGN and LTA from Staphylococcus aureus Induced Inflammation and Decreased Lactation through Regulating DNA Methylation and Histone H3 Acetylation in Bovine Mammary Epithelial Cells

Staphylococcus aureus (S. aureus) and Escherichia coli (E. coli) are the most common pathogens of mastitis, and S. aureus generally causes subclinical mastitis which is more persistent and resistant to treatment. Peptidoglycan (PGN) and lipoteichoic acid (LTA) are cell wall components of S. aureus. Although the roles of PGN and LTA in causing inflammation are well studied, the epigenetic mechanisms of the effects of PGN and LTA on the inflammation and lactation remain poorly understood. This study characterized the gene expression profiling by RNA sequencing and investigated DNA methylation and histone acetylation in relation to inflammation and lactation in the immortalized bovine mammary epithelial cell line (MAC-T). The cells were cultured for 24 h with neither PGN nor LTA (CON), PGN (30 μg/mL), LTA (30 μg/mL), and PGN (30 μg/mL) + LTA (30 μg/mL), respectively. The number of differentially expressed genes (DEGs) and the expression of proinflammatory factors including interleukin (IL)-1β, IL-6, IL-8, chemokine (C-X-C motif) ligand (CXCL)1, and CXCL6 of the treatments increased in the following order: CON < PGN < LTA < PGN + LTA, and the DEGs mainly enriched on the cytokine-cytokine receptor interaction and chemokine signaling pathway. LTA and PGN + LTA induced hypomethylation of global DNA by suppressing DNA methyltransferase (DNMT) activity. PGN and LTA, alone or combined, decreased the mRNA expression of casein genes (CSN1S1, CSN2, and CSN3) and the expression of two caseins (CSN2 and CSN3), and reduced histone H3 acetylation by suppressing histone acetyltransferase (HAT) activity and promoting histone deacetylase (HDAC) activity. Collectively, this study revealed that PGN and LTA induced inflammation probably due to decreasing DNA methylation through regulating DNMT activity, and decreased lactation possibly through reducing histone H3 acetylation by regulating HAT and HDAC activity in bovine mammary epithelial cells.