Genetic effects of single nucleotide polymorphisms in JAK2 and STAT5A genes on susceptibility of Chinese Holsteins to mastitis

Genome-wide association and functional genomic analyses for various hoof health traits in North American Holstein cattle

Hoof diseases are a major welfare and economic issue in the global dairy cattle production industry, which can be minimized through improved management and breeding practices. Optimal genetic improvement of hoof health could benefit from a deep understanding of the genetic background and biological underpinning of indicators of hoof health. Therefore, the primary objectives of this study were to perform genome-wide association studies, using imputed high-density genetic markers data from North American Holstein cattle, for 8 hoof-related traits: digital dermatitis, sole ulcer, sole hemorrhage, white line lesion, heel horn erosion, interdigital dermatitis, interdigital hyperplasia, and toe ulcer, and a hoof health index. De-regressed estimated breeding values from 25,580 Holstein animals were used as pseudo-phenotypes for the association analyses. The genomic quality control, genotype phasing, and genotype imputation were performed using the PLINK (version 1.9), Eagle (version 2.4.1), and Minimac4 software, respectively. The functional genomic analyses were performed using the GALLO R package and the DAVID platform. We identified 22, 34, 14, 22, 28, 33, 24, 43, and 15 significant markers for digital dermatitis, heel horn erosion, interdigital dermatitis, interdigital hyperplasia, sole hemorrhage, sole ulcer, toe ulcer, white line lesion disease, and the hoof health index, respectively. The significant markers were located across all autosomes, except BTA10, BTA12, BTA20, BTA26, BTA27, and BTA28. Moreover, the genomic regions identified overlap with various previously reported quantitative trait loci for exterior, health, meat and carcass, milk, production, and reproduction traits. The enrichment analyses identified 44 significant gene ontology terms. These enriched genomic regions harbor various candidate genes previously associated with bone development, metabolism, and infectious and immunological diseases. These findings indicate that hoof health traits are highly polygenic and influenced by a wide range of biological processes.

Characterization of runs of homozygosity islands in American mink using whole-genome sequencing data

The genome-wide analysis of runs of homozygosity (ROH) islands can be an effective strategy for identifying shared variants within a population and uncovering important genomic regions related to complex traits. The current study performed ROH analysis to characterize the genome-wide patterns of homozygosity, identify ROH islands and annotated genes within these candidate regions using whole-genome sequencing data from 100 American mink (Neogale vison). After sequence processing, variants were called using GATK and Samtools pipelines. Subsequent to quality control, 8,373,854 bi-allelic variants identified by both pipelines remained for further analysis. A total of 34,652 ROH segments were identified in all individuals, among which shorter segments (0.3-1 Mb) were abundant throughout the genome, approximately accounting for 84.39% of all ROH. Within these segments, we identified 63 ROH islands housing 156 annotated genes. The genes located in ROH islands were associated with fur quality (EDNRA, FGF2, FOXA2 and SLC24A4), body size/weight (MYLK4, PRIM2, FABP2, EYS and PHF3), immune capacity (IL2, IL21, PTP4A1, SEMA4C, JAK2, CCNA2 and TNIP3) and reproduction (ADAD1, KHDRBS2, INSL6, PGRMC2 and HSPA4L). Furthermore, Gene Ontology and KEGG pathway enrichment analyses revealed 56 and 9 significant terms (FDR-corrected p-value < 0.05), respectively, among which cGMP-PKG signalling pathway, regulation of actin cytoskeleton, and calcium signalling pathway were highlighted due to their functional roles in growth and fur characteristics. This is the first study to present ROH islands in American mink. The candidate genes from ROH islands and functional enrichment analysis suggest possible signatures of selection in response to the mink breeding targets, such as increased body length, reproductive performance and fur quality. These findings contribute to our understanding of genetic characteristics, and provide complementary information to assist with implementation of breeding strategies for genetic improvement in American mink.

The emerging role of miRNA-122 in infectious diseases: Mechanisms and potential biomarkers

microRNAs (miRNAs) are small, non-coding RNA molecules that play crucial regulatory roles in numerous cellular processes. Recent investigations have highlighted the significant involvement of miRNA-122 (miR-122) in the pathogenesis of infectious diseases caused by diverse pathogens, encompassing viral, bacterial, and parasitic infections. In the context of viral infections, miR-122 exerts regulatory control over viral replication by binding to the viral genome and modulating the host's antiviral response. For instance, in hepatitis B virus (HBV) infection, miR-122 restricts viral replication, while HBV, in turn, suppresses miR-122 expression. Conversely, miR-122 interacts with the hepatitis C virus (HCV) genome, facilitating viral replication. Regarding bacterial infections, miR-122 has been found to regulate host immune responses by influencing inflammatory cytokine production and phagocytosis. In Vibrio anguillarum infections, there is a significant reduction in miR-122 expression, contributing to the pathophysiology of bacterial infections. Toll-like receptor 14 (TLR14) has been identified as a novel target gene of miR-122, affecting inflammatory and immune responses. In the context of parasitic infections, miR-122 plays a crucial role in regulating host lipid metabolism and immune responses. For example, during Leishmania infection, miR-122-containing extracellular vesicles from liver cells are unable to enter infected macrophages, leading to a suppression of the inflammatory response. Furthermore, miR-122 exhibits promise as a potential biomarker for various infectious diseases. Its expression level in body fluids, particularly in serum and plasma, correlates with disease severity and treatment response in patients affected by HCV, HBV, and tuberculosis. This paper also discusses the potential of miR-122 as a biomarker in infectious diseases. In summary, this review provides a comprehensive and insightful overview of the emerging role of miR-122 in infectious diseases, detailing its mechanism of action and potential implications for the development of novel therapeutic strategies.

Polymorphism in neutrophil cytosolic factor 4 (NCF4) of dairy cows had mastitis in previous lactations, and the relationship with the respiratory burst

Nicotinamide adenine dinucleotide phosphate oxidase (NADPH oxidase), as a key factor in innate immunity, consists of several components, one of them is p40phox which is encoded by neutrophil cytosolic factor 4 (NCF4). Respiratory burst and reactive oxygen species (ROS) production are antimicrobial mechanisms associated with NADPH oxidase. This study evaluated the effects of g.18174 A > G and g.18270C > T single-nucleotide polymorphisms (SNP) in NCF4 on bovine mastitis and the respiratory burst capacity of neutrophils. SNPs of 160 dairy cattle were determined using a novel PCR-RFLP protocol by employing restriction enzymes, MboI and FokI. Also, the flow cytometry measured respiratory burst in 82 blood samples. Our results indicated that only g.18174 A > G SNP reduced the respiratory burst capacity. However, both SNPs were not significantly correlated with clinical mastitis. We concluded that g.18174 A > G decreases the function of NADPH oxidase. However, both SNPs were not significantly correlated with clinical mastitis.

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.

Candidate genes for mastitis resistance in dairy cattle: a data integration approach

BACKGROUND

Inflammation of the mammary tissue (mastitis) is one of the most detrimental health conditions in dairy ruminants and is considered the most economically important infectious disease of the dairy sector. Improving mastitis resistance is becoming an important goal in dairy ruminant breeding programmes. However, mastitis resistance is a complex trait and identification of mastitis-associated alleles in livestock is difficult. Currently, the only applicable approach to identify candidate loci for complex traits in large farm animals is to combine different information that supports the functionality of the identified genomic regions with respect to a complex trait.

METHODS

To identify the most promising candidate loci for mastitis resistance we integrated heterogeneous data from multiple sources and compiled the information into a comprehensive database of mastitis-associated candidate loci. Mastitis-associated candidate genes reported in association, expression, and mouse model studies were collected by searching the relevant literature and databases. The collected data were integrated into a single database, screened for overlaps, and used for gene set enrichment analysis.

RESULTS

The database contains candidate genes from association and expression studies and relevant transgenic mouse models. The 2448 collected candidate loci are evenly distributed across bovine chromosomes. Data integration and analysis revealed overlaps between different studies and/or with mastitis-associated QTL, revealing promising candidate genes for mastitis resistance.

CONCLUSION

Mastitis resistance is a complex trait influenced by numerous alleles. Based on the number of independent studies, we were able to prioritise candidate genes and propose a list of the 22 most promising. To our knowledge this is the most comprehensive database of mastitis associated candidate genes and could be helpful in selecting genes for functional validation studies.

Characterization of peripheral white blood cells transcriptome to unravel the regulatory signatures of bovine subclinical mastitis resistance

Subclinical bovine mastitis is a pathogenic infection of the breast characterized by a marked decrease in milk production and quality. As it has no obvious clinical symptoms, diagnosis and treatment are challenging. Therefore, searching for biomarkers in cows’ peripheral white blood cells is valuable for preventing and treating subclinical mastitis. Thus, in this study, the transcriptome of peripheral blood from healthy and subclinical mastitis cows was characterized to find the regulatory signatures of bovine subclinical mastitis using RNA-seq. A total of 287 differentially expressed genes (DEGs) and 70 differentially expressed lncRNAs (DELs) were detected, and 37 DELs were documented near known Quantitative Trait Loci (QTL) associated with the mastitis of cows. Bioinformatic analysis indicated that lncRNAs MSTRG25101.2, MSTRG.56327.1, and MSTRG.18968.1, which are adjacent to the SCS QTL and SCC QTL, may be candidate lncRNAs that influence the pathogenesis of mastitis in cows by up-regulating the expression of genes TLR4, NOD2, CXCL8, and OAS2. Moreover, the alternative splicing (AS) pattern of transcriptional sequence differences between healthy cows and subclinical mastitis cows suggested a molecular mechanism of mastitis resistance and susceptibility. A total of 2,212 differential alternative splicing (DAS) events, corresponding to 1,621 unique DAS genes, were identified in both groups and significantly enriched in immune and inflammatory pathways. Of these, 29 DAS genes were subject to regulation by 32 alternative splicing SNPs, showing diverse and specific splicing patterns and events. It is hypothesized that the PIK3C2B and PPRPF8 splice variants associated with AS SNPs (rs42705933 and rs133847062) may be risk factors for susceptibility to bovine subclinical mastitis. Altogether, these key blood markers associated with resistance to subclinical mastitis and SNPs associated with alternative splicing of genes provide the basis for genetic breeding for resistance to subclinical mastitis in cows.

Induction of mastitis by cow-to-mouse fecal and milk microbiota transplantation causes microbiome dysbiosis and genomic functional perturbation in mice

BACKGROUND

Mastitis pathogenesis involves a wide range of opportunistic and apparently resident microorganims including bacteria, viruses and archaea. In dairy animals, microbes reside in the host, interact with environment and evade the host immune system, providing a potential for host-tropism to favor mastitis pathogenesis. To understand the host-tropism phenomena of bovine-tropic mastitis microbiomes, we developed a cow-to-mouse mastitis model.

METHODS

A cow-to-mouse mastitis model was established by fecal microbiota transplantation (FMT) and milk microbiota transplantation (MMT) to pregnant mice to assess microbiome dysbiosis and genomic functional perturbations through shotgun whole metagenome sequencing (WMS) along with histopathological changes in mice mammary gland and colon tissues.

RESULTS

The cow-to-mouse FMT and MMT from clinical mastitis (CM) cows induced mastitis syndromes in mice as evidenced by histopathological changes in mammary gland and colon tissues. The WMS of 24 samples including six milk (CM = 3, healthy; H = 3), six fecal (CM = 4, H = 2) samples from cows, and six fecal (CM = 4, H = 2) and six mammary tissue (CM = 3, H = 3) samples from mice generating 517.14 million reads (average: 21.55 million reads/sample) mapped to 2191 bacterial, 94 viral and 54 archaeal genomes. The Kruskal-Wallis test revealed significant differences (p = 0.009) in diversity, composition, and relative abundances in microbiomes between CM- and H-metagenomes. These differences in microbiome composition were mostly represented by Pseudomonas aeruginosa, Lactobacillus crispatus, Klebsiella oxytoca, Enterococcus faecalis, Pantoea dispersa in CM-cows (feces and milk), and Muribaculum spp., Duncaniella spp., Muribaculum intestinale, Bifidobacterium animalis, Escherichia coli, Staphylococcus aureus, Massilia oculi, Ralstonia pickettii in CM-mice (feces and mammary tissues). Different species of Clostridia, Bacteroida, Actinobacteria, Flavobacteriia and Betaproteobacteria had a strong co-occurrence and positive correlation as the indicator species of murine mastitis. However, both CM cows and mice shared few mastitis-associated microbial taxa (1.14%) and functional pathways regardless of conservation of mastitis syndromes, indicating the higher discrepancy in mastitis-associated microbiomes among lactating mammals.

CONCLUSIONS

We successfully induced mastitis by FMT and MMT that resulted in microbiome dysbiosis and genomic functional perturbations in mice. This study induced mastitis in a mouse model through FMT and MMT, which might be useful for further studies- focused on pathogen(s) involved in mastitis, their cross-talk among themselves and the host.

Association of Aberrant DNA Methylation Level in the CD4 and JAK-STAT-Pathway-Related Genes with Mastitis Indicator Traits in Chinese Holstein Dairy Cattle

The present study was designed to evaluate the gene expression and DNA methylation level in the promoter region of the CD4 and the JAK-STAT-pathway-related genes. A total of 24 samples were deployed in the gene expression and 118 samples were used in the DNA methylation study. Student's t-tests were used to analyze the gene expression and DNA methylation. The evaluation of DNA methylation in promoter regions of JAK2 and STAT5A revealed hypo-methylation levels of CpG sites and higher gene expression in cows diagnosed with mastitis as compared to the healthy control, and vice versa in those with CD4. DNA methylation was negatively correlated with gene expression in JAK2, STAT5A, and CD4 genes. Six, two, and four active transcription factors were identified on the CpG sites in the promoter regions of JAK2, STAT5A, and CD4 genes, respectively. Regarding correlation analysis, the DNA methylation levels of CD4 showed significantly higher positive correlations with somatic cell counts (p < 0.05). Findings of the current study inferred that aberrant DNA methylation in the CpG sites at the 1 kb promoter region in JAK2, STAT5A, and CD4 genes due to mastitis in cows can be used as potential epigenetic markers to estimate bovine mastitis susceptibility in dairy cattle.

Elucidating the Role of Serum tRF-31-U5YKFN8DYDZDD as a Novel Diagnostic Biomarker in Gastric Cancer (GC)

Background

Gastric cancer (GC) is one of the malignant tumors with the highest morbidity and mortality in the world. Early diagnosis combined with surgical treatment can significantly improve the prognosis of patients. Therefore, it is urgent to seek higher sensitivity and specificity biomarkers in GC. tRNA-derived small RNAs are a new non-coding small RNA that widely exists in tumor cells and body fluids. In this study, we explore the expression and biological significance of tRNA-derived small RNAs in GC.

Materials and Methods

First of all, we screened the differentially expressed tRNA-derived small RNAs in tumor tissues by high-throughput sequencing. Agarose gel electrophoresis (AGE), Sanger sequencing, and Nuclear and Cytoplasmic RNA Separation Assay were used to screen tRF-31-U5YKFN8DYDZDD as a potential tumor biomarker for the diagnosis of GC. Then, we detected the different expressions of tRF-31-U5YKFN8DYDZDD in 24 pairs of GC and paracancerous tissues, the serum of 111 GC patients at first diagnosis, 89 normal subjects, 48 superficial gastritis patients, and 28 postoperative GC patients by quantitative real-time PCR (qRT-PCR). Finally, we used the receiver operating characteristic (ROC) curve to analyze its diagnostic efficacy.

Results

The expression of tRF-31-U5YKFN8DYDZDD has good stability and easy detection. tRF-31-U5YKFN8DYDZDD was highly expressed in tumor tissue, serum, and cell lines of GC, and the expression was significantly related to TNM stage, depth of tumor invasion, lymph node metastasis, and vascular invasion. The expression of serum tRF-31-U5YKFN8DYDZDD in the GC patients decreased after the operation (P = 0.0003). Combined with ROC curve analysis, tRF-31-U5YKFN8DYDZDD has better detection efficiency than conventional markers.

Conclusions

The expressions of tRF-31-U5YKFN8DYDZDD in the tumor and paracancerous tissues, the serum of GC patients and healthy people, and the serum of GC patients before and after operation were different. tRF-31-U5YKFN8DYDZDD is not only a diagnostic biomarker of GC but also a predictor of poor prognosis.

Genome-Wide Association Studies of Somatic Cell Count in the Assaf Breed

Simple Summary

Mastitis causes economic loss due to discarded milk and reduced milk production and quality, increased medical care costs and somatic cell count (SCC) penalties. The use of genetic markers associated with the variability of this trait through marker-assisted selection (MAS) could help traditional methods. Our objectives were to identify new single nucleotide polymorphisms (SNPs) and genes associated with mastitis resistance in Assaf sheep by using the Illumina Ovine Infinium® HD SNP BeadChip (680K). Firstly, corrected phenotype estimates for somatic cell score (SCS) were calculated using 6173 records from 1894 multiparous Assaf ewes, and were used to select 192 extreme animals (low SCS group: n = 96; and high SCS group: n = 96) for the genome-wide association study (GWAS). Four SNPs (rs419096188, rs415580501, rs410336647, and rs424642424), three of them totally linked, were found to be significant at the chromosome level (FDR 10%) in two different regions of OAR19 close to genes related to the immune system response. Validation studies of two SNPs (rs419096188 and rs424642424) by Kompetitive Allele-Specific PCR (KASP) genotyping in the total population (n = 1894) confirmed previous GWAS association results for the SCS trait. Finally, the SNP rs419096188 was also associated with lactose content trait.

Abstract

A genome-wide association study (GWAS) was performed to identify new single nucleotide polymorphisms (SNPs) and genes associated with mastitis resistance in Assaf sheep by using the Illumina Ovine Infinium^® HD SNP BeadChip (680K). In total, 6173 records from 1894 multiparous Assaf ewes with at least three test day records and aged between 2 and 7 years old were used to estimate a corrected phenotype for somatic cell score (SCS). Then, 192 ewes were selected from the top (n = 96) and bottom (n = 96) tails of the corrected SCS phenotype distribution to be used in a GWAS. Although no significant SNPs were found at the genome level, four SNPs (rs419096188, rs415580501, rs410336647, and rs424642424) were significant at the chromosome level (FDR 10%) in two different regions of OAR19. The SNP rs419096188 was located in intron 1 of the NUP210 and close to the HDAC11 genes (61 kb apart), while the other three SNPs were totally linked and located 171 kb apart from the ARPP21 gene. These three genes were related to the immune system response. These results were validated in two SNPs (rs419096188 and rs424642424) in the total population (n = 1894) by Kompetitive Allele-Specific PCR (KASP) genotyping. Furthermore, rs419096188 was also associated with lactose content.

Mechanisms by which mastitis affects reproduction in dairy cow: A review

Reproductive performance is a key factor in determining the profitability of dairy farm, which is affected by many factors such as environment and diseases. Mastitis is a common and important disease, which has caused huge economic losses to the dairy industries worldwide. Mammary gland infection causes immune responses, resulting in the abnormal secretion of cytokines and hormones and abnormal function of the reproductive system such as the ovary, corpus luteum, uterus and embryo. Cows with mastitis have delayed oestrus, decreased pregnancy rate and increased risk of abortion. The adverse effects of mastitis on reproductive performance are affected by many factors, such as occurrence time, pathogen and cow factors. This paper primarily reviews the progress in the effects and mechanisms of mastitis on reproductive performance, with emphasis on maternal transcriptome, genomic analysis, epigenetic modification, microbiota, inflammatory regulation and immune evasion mechanism of mastitis, aiming to provide directions for the prevention and control of mastitis in the future.

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.

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

Background

Mastitis in dairy cows caused by Staphylococcus aureus is a major problem hindering economic growth in dairy farms worldwide. It is difficult to prevent or eliminate due to its asymptomatic nature and long persistence of infection. Although transcriptomic responses of bovine mammary gland cells to pathogens that cause mastitis have been studied, the common responses of peripheral blood leukocytes to S. aureus infection across two consecutive generations of dairy cattle have not been investigated.

Methods

In the current study, RNA-Seq was used to profile the transcriptomes of peripheral blood leukocytes sampled from S. aureus-infected mothers and their S. aureus-infected daughters, and also healthy non-infected mothers and their healthy daughters. Differential gene expression was evaluated as follows: 1) S. aureus-infected cows versus healthy non-infected cows (S vs. H, which include all the mothers and daughters), 2) S. aureus-infected mothers versus healthy non-infected mothers (SM vs. HM), and 3) S. aureus-infected daughters versus healthy non-infected daughters (SMD vs. HMD).

Results

Analysis of all identified expressed genes in the four groups (SM, SMD, HM, and HMD) showed that EPOR, IL9, IFNL3, CCL26, IL26 were exclusively expressed in both the HM and HMD groups, and that they were significantly (P <  0.05) enriched for the cytokine-cytokine receptor interaction pathway. A total of 17, 13 and 10 differentially expressed genes (DEGs) (FDR P_adj. < 0.1 and |FC| > 1.2) were detected in the three comparisons, respectively. DEGs with P <  0.05 and |FC| > 2 were used for functional enrichment analyses. For the S vs. H comparison, DEGs detected included CCL20, IL13 and MMP3, which are associated with the IL-17 signaling pathway. In the SM vs. HM and SMD vs. HMD comparisons, five (BLA-DQB, C1R, C2, FCGR1A, and KRT10) and six (BLA-DQB, C3AR1, CFI, FCAR, FCGR3A, and LOC10498484) genes, respectively, were involved in the S. aureus infection pathway.

Conclusions

Our study provides insights into the transcriptomic responses of bovine peripheral blood leukocytes across two generations of cattle naturally infected with S. aureus. The genes highlighted in this study could serve as expression biomarkers for mastitis and may also contain sequence variation that can be used for genetic improvement of dairy cattle for resilience to mastitis.

Association of GH, STAT5A, MYF5 gene polymorphisms with milk somatic cell count, EC and pH levels of Holstein dairy cattle

This study was conducted to ivnestigate the associations of GH-AluI, STAT5A-AvaI and MYF5-TaqI gene polymorphisms with milk somatic cell count (SCC), electrical conductivity (EC) and pH levels in Holstein dairy cows. For this purpose, 167 blood and 1670 milk samples of 167 Holstein cows in their 2nd lactation were used. There were significant relationships between GH-AluI genotypes and milk EC (p < 0.001) and between STAT5A-AvaI genotypes and milk EC (p = 0.007), but there were not any significant relationships between MYF5 gene polymorphism and the investigated traits (p > 0.05). The greatest EC values were observed in GH-AluI-LV and STAT5A-AvaI-TT-genotyped individuals. Just because of association of EC with mastitis, it was concluded that present GH-AluI and STAT5A-AvaI polymorphisms could be used in further studies to be conducted to improve mastitis resistance and milk quality traits of Holstein dairy cows.

Polymorphisms in JAK2 Gene are Associated with Production Traits and Mastitis Resistance in Dairy Cattle

The present study was designed to investigate the effects of single nucleotide polymorphisms (SNPs) in the JAK2 gene on the production and mastitis related traits in dairy cattle. Blood and milk samples were collected from 201 lactating dairy cattle of three breeds, i.e. Holstein Friesian (HF), Jersey (J) and Achai (A) and their crosses maintained at well-established dairy farms in Khyber Pakhtunkhwa, Pakistan. Generalized linear model was used to evaluate the association between genotypes and the studied traits. A DNA pool was made from randomly selected 30 samples which revealed three SNPs, i.e. SNP 1 in 5’ upstream region (G>A, rs379754157), SNP 2 in intron 15 (A>G, rs134192265), and SNP 3 in exon 20 (A>G, rs110298451) that were further validated in the population under study using SNaPshot technique. Of the three SNPs, SNP 1 did not obey Hardy-Weinberg equilibrium (P<0.05). SNP 2 and SNP 3 were found to be in strong linkage disequilibrium and allele G was highly prevalent compared to allele A in these SNPs. in SNP 1, the GG genotype was associated with significantly (P<0.01) higher SCC, whereas SNP 2 and SNP 3 were significantly (P<0.01) associated with higher lactose percentage compared to the other geno-types. The haplogroups association analysis revealed that H1H2 (GG GG AG) has significantly lower SCC than H2H2 (GG GG GG). The results infer that JAK2 could be an important candidate gene and the studied SNPs might be useful genetic markers for production and mastitis related traits.

Significant genetic effects of JAK2 and DGAT1 mutations on milk fat content and mastitis resistance in Holsteins

Improving the production traits and resistance against mastitis in dairy cattle is a challenge for animal scientists across the globe. The present study was designed to investigate the genetic effects of single nucleotide polymorphisms (SNPs) in Janus kinase 2 (JAK2) and diacylglycerol acyltransferase (DGAT1) genes with production and mastitis-related traits. Four SNPs in JAK2 and one in DGAT1 were analyzed through Chinese Cow's SNPs Chip-I (CCSC-I) and genotyped in a population of 312 Chinese Holsteins. Our findings demonstrated that milk fat percentage, somatic cell count (SCC), somatic cell score (SCS), serum cytokines interleukin 6 (IL-6) and interferon gamma (IFN-γ) showed significant associations (P < 0.05) with at least one or more identified SNPs. Consequently, the analysis based on haplotypes amongst the SNPs in JAK2 revealed noteworthy (P < 0.05) association with SCC and IL-6. Collectively, our results verified the pleiotropic ability of detected SNPs in bovine JAK2 and DGAT1 for milk fat percentage as well as mastitis-related traits. The significant SNPs in both the genes could serve as powerful genetic markers to minimize mastitis risk. In addition, besides SCC and SCS, the IFN-γ and IL-6 could also be used as indicators of improved genetic resistance against mastitis.

Single-step genome-wide association for longitudinal traits of Canadian Ayrshire, Holstein, and Jersey dairy cattle

Estimating single nucleotide polymorphism (SNP) effects over time is essential to identify and validate candidate genes (or quantitative trait loci) associated with time-dependent variation of economically important traits and to better understand the underlying mechanisms of lactation biology. Therefore, in this study, we aimed to estimate time-dependent effects of SNP and identifying candidate genes associated with milk (MY), fat (FY), and protein (PY) yields, and somatic cell score (SCS) in the first 3 lactations of Canadian Ayrshire, Holstein, and Jersey breeds, as well as suggest their potential pattern of phenotypic effect over time. Random regression coefficients for the additive direct genetic effect were estimated for each animal using single-step genomic BLUP, based on 2 random regression models: one considering MY, FY, and PY in the first 3 lactations and the other considering SCS in the first 3 lactations. Thereafter, SNP solutions were obtained for random regression coefficients, which were used to estimate the SNP effects over time (from 5 to 305 d in lactation). The top 1% of SNP that showed a high magnitude of SNP effect in at least 1 d in lactation were selected as relevant SNP for further analyses of candidate genes, and clustered according to the trajectory of their SNP effects over time. The majority of SNP selected for MY, FY, and PY increased the magnitude of their effects over time, for all breeds. In contrast, for SCS, most selected SNP decreased the magnitude of their effects over time, especially for the Holstein and Jersey breeds. In general, we identified a different set of candidate genes for each breed, and similar genes were found across different lactations for the same trait in the same breed. For some of the candidate genes, the suggested pattern of phenotypic effect changed among lactations. Among the lactations, candidate genes (and their suggested phenotypic effect over time) identified for the second and third lactations were more similar to each other than for the first lactation. Well-known candidate genes with major effects on milk production traits presented different suggested patterns of phenotypic effect across breeds, traits, and lactations in which they were identified. The candidate genes identified in this study can be used as target genes in studies of gene expression.

Overexpression of lncRNA H19 changes basic characteristics and affects immune response of bovine mammary epithelial cells

The function of long non-coding RNA H19 (H19) on cell proliferation has been observed in various cell types, and the increased expression of H19 was also found in the lipopolysaccharide (LPS)-induced inflammatory bovine mammary epithelial cells (MAC-T). However, the roles of H19 in the inflammatory response and physiological functions of bovine mammary epithelial cell are not clear. In the present study, we found that overexpression of H19 in MAC-T cells significantly promoted cell proliferation, increased the protein and mRNA level of β-casein, and enhanced the expression of tight junction (TJ)-related proteins while inhibited staphylococcus aureus adhesion to cells. In addition, results demonstrated that overexpression of H19 affected the LPS-induced immune response of MAC-T cells by promoting expressions of inflammatory factors, including TNF-α, IL-6, CXCL2 and CCL5, and activating the NF-κB signal pathway. Our findings indicate that H19 is likely to play an important role in maintaining normal functions and regulating immune response of bovine mammary epithelial cells.

Prioritizing candidate genes post-GWAS using multiple sources of data for mastitis resistance in dairy cattle

Background

Improving resistance to mastitis, one of the costliest diseases in dairy production, has become an important objective in dairy cattle breeding. However, mastitis resistance is influenced by many genes involved in multiple processes, including the response to infection, inflammation, and post-infection healing. Low genetic heritability, environmental variations, and farm management differences further complicate the identification of links between genetic variants and mastitis resistance. Consequently, studies of the genetics of variation in mastitis resistance in dairy cattle lack agreement about the responsible genes.

Results

We associated 15,552,968 imputed whole-genome sequencing markers for 5147 Nordic Holstein cattle with mastitis resistance in a genome-wide association study (GWAS). Next, we augmented P-values for markers in genes in the associated regions using Gene Ontology terms, Kyoto Encyclopedia of Genes and Genomes pathway analysis, and mammalian phenotype database. To confirm results of gene-based analyses, we used gene expression data from E. coli-challenged cow udders. We identified 22 independent quantitative trait loci (QTL) that collectively explained 14% of the variance in breeding values for resistance to clinical mastitis (CM). Using association test statistics with multiple pieces of independent information on gene function and differential expression during bacterial infection, we suggested putative causal genes with biological relevance for 12 QTL affecting resistance to CM in dairy cattle.

Conclusion

Combining information on the nearest positional genes, gene-based analyses, and differential gene expression data from RNA-seq, we identified putative causal genes (candidate genes with biological evidence) in QTL for mastitis resistance in Nordic Holstein cattle. The same strategy can be applied for other traits.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5050-x) contains supplementary material, which is available to authorized users.

Genome-wide association mapping for type and mammary health traits in French dairy goats identifies a pleiotropic region on chromosome 19 in the Saanen breed

Type traits and mammary health traits are important to dairy ruminant breeding because they influence animal health, milking ability, and longevity, as well as the economic sustainability of farms. The availability of the genomic sequence and a single nucleotide polymorphism chip in goats has opened up new fields of investigation to better understand the genes and mechanisms that underlie such complex traits and to be able to select them. Our objective was to perform a genome-wide association study in dairy goats for 11 type traits and somatic cell count (SCC) as proxies for mastitis resistance. A genome-wide association study was implemented using a daughter design composed of 1,941 Alpine and Saanen goats sired by 20 artificial insemination bucks, genotyped with the Illumina GoatSNP50 BeadChip (Illumina Inc., San Diego, CA). This association study was based on both linkage analyses and linkage disequilibrium using QTLmap software (http://dga7.jouy.inra.fr/qtlmap/) interval mapping was performed with the likelihood ratio test using linear regressions. Breeds were analyzed together and separately. The study highlighted 37 chromosome-wide significant quantitative trait loci (QTL) with linkage analyses and 222 genome-wide significant QTL for linkage disequilibrium, for type and SCC traits in dairy goats. Genomic control of those traits was mostly polygenic and breed-specific, suggesting that within-breed selection would be favored for those traits. Of note, Capra hircus autosome (CHI) 19 appeared to be highly enriched in single nucleotide polymorphisms associated with type and SCC, with 2 highly significant regions in the Saanen breed. One region (33-42 Mb) was significantly associated with SCC and includes candidate genes associated with response to intramammary infections (RARA, STAT3, STAT5A, and STAT5B). Another region of the CHI 19 (24.5-27 Mb) exhibited an adverse pleiotropic effect on milk production (milk, fat yield, and protein yield) and udder traits (udder floor position and rear udder attachment) that agreed with the negative genetic correlations that exist between those 2 groups of traits. These QTL were not found in the Alpine breed. In Alpine, the 2 most significant regions were associated with chest depth on CHI 6 (45.8-46.0 Mb) and CHI 8 (80.7-81.1 Mb). These results will be helpful for goat selection in the future and could lead to identification of causal mutations.

A Point Mutation in Suppressor of Cytokine Signalling 2 (Socs2) Increases the Susceptibility to Inflammation of the Mammary Gland while Associated with Higher Body Weight and Size and Higher Milk Production in a Sheep Model

Mastitis is an infectious disease mainly caused by bacteria invading the mammary gland. Genetic control of susceptibility to mastitis has been widely evidenced in dairy ruminants, but the genetic basis and underlying mechanisms are still largely unknown. We describe the discovery, fine mapping and functional characterization of a genetic variant associated with elevated milk leukocytes count, or SCC, as a proxy for mastitis. After implementing genome-wide association studies, we identified a major QTL associated with SCC on ovine chromosome 3. Fine mapping of the region, using full sequencing with 12X coverage in three animals, provided one strong candidate SNP that mapped to the coding sequence of a highly conserved gene, suppressor of cytokine signalling 2 (Socs2). The frequency of the SNP associated with increased SCC was 21.7% and the Socs2 genotype explained 12% of the variance of the trait. The point mutation induces the p.R96C substitution in the SH2 functional domain of SOCS2 i.e. the binding site of the protein to various ligands, as well-established for the growth hormone receptor GHR. Using surface plasmon resonance we showed that the p.R96C point mutation completely abrogates SOCS2 binding affinity for the phosphopeptide of GHR. Additionally, the size, weight and milk production in p.R96C homozygote sheep, were significantly increased by 24%, 18%, and 4.4%, respectively, when compared to wild type sheep, supporting the view that the point mutation causes a loss of SOCS2 functional activity. Altogether these results provide strong evidence for a causal mutation controlling SCC in sheep and highlight the major role of SOCS2 as a tradeoff between the host’s inflammatory response to mammary infections, and body growth and milk production, which are all mediated by the JAK/STAT signaling pathway.

Mastitis is an inflammation of the mammary gland mainly caused by invading bacteria. Ruminants show natural variability in their predisposition to mastitis, and therefore provide unique models for study of the genetics and physiology of host response to bacterial infection. A genome-wide association study was conducted in a dairy sheep population for milk somatic cell counts as a proxy for mastitis. Fine mapping, using whole genome sequencing, led to the identification of a mutation in the Suppressor of Cytokine Signaling 2 gene (socs2). This mutation was shown to cause a loss of functional activity of the SOCS2 protein, which suggested impairment of feedback control of the JAK/STAT signaling pathways in susceptible animals. Additionally, size, weight and milk production were increased in animals carrying the susceptible variant suggesting a pleiotropic effect of the gene on production versus health traits. Results gave strong evidence of the role of SOCS2 in the host’s inflammation of the udder and provided new insights into the key mechanisms underlying the genetic control of mastitis.

Association study of single nucleotide polymorphisms in JAK2 and STAT5B genes and their differential mRNA expression with mastitis susceptibility in Chinese Holstein cattle

The JAK-STAT pathway plays a key role in mediating immune responses. The genetic effects of single nucleotide polymorphisms (SNPs) in JAK2 and STAT5B were investigated for serum cytokines, mastitis indicators and productions traits in a population of 468 Chinese Holstein cattle. Pooled DNA sequencing revealed one SNP (BTA8:g.39645396A>G) in JAK2 and two SNPs (BTA19:g.43673888A>G and BTA19:g.43660093T>C) in STAT5B. A fixed effect model considering the effects of SNPs, parity, herd, season and year of calving was used by way of the general linear model procedure of sas. Genotype frequencies of these SNPs in the population were in Hardy-Weinberg equilibrium (P > 0.05). A novel SNP (g.39645396A>G) in JAK2 was predicted to change the amino acid from lysine to asparagine and was significantly associated with the somatic cell count (SCC) and somatic cell score (SCS), whereas g.43673888A>G in STAT5B was significantly associated with SCC, SCS and interleukin-4 (IL-4) (P < 0.05). The dominant effect of g.39645396A>G in JAK2 was significant for SCS, and its additive effect was significant for SCC, whereas the dominant effect of g.43673888A>G in STAT5B was significant for SCS and IL-4 (P < 0.05). The combination of g.39645396A>G in JAK2 and g.43673888A>G in STAT5B showed a significant effect on SCC, SCS, IL-4 and TNF-α (P < 0.05). As for mRNA expression analysis, the AA genotype g.39645396A>G and GG genotype g.43673888A>G indicated higher mRNA expression level and were significantly different from other genotypes (P < 0.05). The results imply that JAK2 and STAT5B genes could be useful candidate genes, and the identified polymorphisms might potentially be strong genetic markers for selection of dairy cattle against mastitis development.