Genetic mechanisms regulating the host response during mastitis

Cytokine interactions and chemokine dysregulations in mastitis immunopathogenesis: insights from transcriptomic profiling of milk somatic cells in tropical Sahiwal (Bos indicus) cows

Introduction

Bovine mastitis causes a significant loss to the dairy industry by affecting the quantity and quality of milk. Addressing this challenge, the present study will leverage advanced omics techniques for early mastitis detection in early lactating Sahiwal cows (Bos indicus). This was the first differential transcriptomic study investigating the alterations in gene expression in milk somatic cells during the progression of naturally occurring mastitis in indigenous Sahiwal cows.

Methods

Cows were grouped into healthy (H), subclinical mastitis (SCM) and clinical mastitis (CM) groups by thoroughly screening them using the California Mastitis Test (CMT) and milk somatic cell counts (SCC). This was followed by detailed milk composition analysis, differential leukocyte counts (DLC), and microbiological culture.

Results

The differential gene expression of milk SCs through transcriptome profiling identified 83 and 76, up-regulated and 157 and 192 down-regulated genes in CM vs H and SCM vs H groups (log2 fold change ≥1 and ≤-1, p < 0.05) respectively. Pathway analysis revealed that upregulated genes were enriched in pathways such as phagosome activity, IL-17 signalling, Th1 and Th2 cell differentiation, while downregulated genes were linked to RIG-I-like receptor signalling, NK cell cytotoxicity, and Toll-like receptor signalling and Cytokine-cytokine receptor interactions. Notably, the study underscores the roles of chemokines CCL8, CCL2, and CXCL10 in immune cell recruitment during mastitis, where their downregulation suggests impaired mammary immune defense that governs Chemokine signalling pathways. Further, the comparative analysis with the previously available milk SCs proteome data identified the downregulation of chemokines signalling pathways during mastitis.

Discussion

Overall, this research enhances our understanding of mastitis pathogenesis and emphasizes that these targeted chemokines may boost mammary resilience through immunomodulation, genetic selection and genome editing or by utilising adjuvants in vaccine development that restore chemokine signalling offers a potential strategy to improve mastitis resistance in dairy cattle.

LncRNA HULIB promotes LPS induced inflammatory response in bovine mammary epithelial cells via PP2AB

Bovine mastitis is regulated by genetic and environmental factors. Long non-coding RNAs (LncRNAs), which regulate various biological processes (immune system and biological development), have been found to play a role in bovine mammary inflammation responses. Here, a novel functional lncRNA, named lncRNA HULIB, was identified as a regulator during bovine mastitis. qRT-PCR and subcellular fractionation assays showed that lncRNA HULIB was significantly up-regulated in LPS-induced bMECs and was mainly localized in the cytoplasm. Gain- or loss-of-function experiments demonstrated that an increase in lncRNA HULIB expression elevated the expression of TLR4 and NF-κB1, which enhanced NF-κB activity, promoting the expression of pro-inflammatory cytokines (IL-6, IL-8, IL-1β, etc) and apoptosis-related genes (BAX, CASP9 and CASP3, etc), while the expression of proliferation-related genes (PCNA, Cyclin D1, Cyclin D2, CDK4 and CDK2) was down-regulated. Ultimately, these changes exacerbated the LPS-induced inflammatory response. Mechanistically, RNA pull-down and RNA immunoprecipitation (RIP) assays revealed that lncRNA HULIB could directly bind the PP2AB protein to regulate inflammatory responses. Overall, lncRNA HULIB is a pro-inflammatory regulator, and its silencing can alleviate the inflammatory responses of bMECs, providing a potential strategy for molecular therapy of bovine mastitis.

Unraveling Dynamic Transcriptomic Changes in Sheep's Lactating Mammary Gland Following Escherichia coli Lipopolysaccharide Exposure

Mammary gland infections constitute a significant challenge in dairy sheep, impacting productivity and welfare. Temporal RNA-Seq provide a valuable approach to evaluate the evolution of the host defensive molecular mechanisms triggered by mastitis caused by external agents or events. This study aimed to characterize the transcriptomic response of sheep mammary glands to an intramammary inflammation induced with an Escherichia coli lipopolysaccharide (LPS) inoculation based on RNA-Seq samples generated from milk somatic cells collected at 3 time points: pre-inoculation (0 h), and 6 h and 24 h post-LPS inoculation. The differential expression analyses between the analyzed time points were performed using 2 statistical approaches: one parametric (DESeq2) and one non-parametric (Wilcoxon rank sum test). The differentially expressed genes (DEGs) commonly identified by both approaches encompass 5,872 for the 0 h versus 6 h comparison, 4,063 for the 0 h versus 24 h comparison, and 1,034 for the 6 h versus 24 h comparison. At both 6 h and 24 h, transcriptomic data highlighted a significant decrease in the expression of genes linked to metabolic processes crucial for milk protein and lipid synthesis within the mammary gland. Concurrently, increased expression of genes related to the neutrophil attraction was observed for 6 and 24 h, with differences in gene expression between DEGs with the highest expression at 6 h, related to T cell activation, type I interferon-mediated signaling pathway, and 24 h, related to cell-cell neutrophil adhesion extravasation or epithelial cell proliferation. In summary, this study reveals how the sheep mammary gland transcriptome responds dynamically to an LPS inoculation, providing a comprehensive understanding of how gene expression patterns evolve over time and shedding light on the molecular mechanisms driving the initial defensive response of the mammary gland against potential inflammatory challenges.

Identification of pasteurized mare milk and powder adulteration with bovine milk using quantitative proteomics and metabolomics approaches

Adulteration in dairy products presents food safety challenges, driven by economic factors. Processing may change specific biomarkers, thus affecting their effectiveness in detection. In this study, proteomics and metabolomics approaches were to investigate the detection of bovine milk (BM) constituents adulteration in pasteurized mare milk (PMM) and mare milk powder (MMP). Several bovine proteins and metabolites were identified, with their abundances in PMM and MMP increasing upon addition of BM. Proteins like osteopontin (OPN) and serotransferrin (TF) detected adulteration down to 1 % in PMM, whereas these proteins in MMP were utilized to identify 10 % adulteration. Biotin and N6-Me-adenosine were effective in detecting adulteration in PMM as low as 10 % and 1 % respectively, while in MMP, their detection limits extend down to 0.1 %. These findings offer insights for authenticating mare milk products and underscore the influence of processing methods on biomarker levels, stressing the need to consider these effects in milk product authentication.

The role of interleukin-10 receptor alpha (IL10Rα) in Mycobacterium avium subsp. paratuberculosis infection of a mammary epithelial cell line

BACKGROUND

Johne's disease is a chronic wasting disease caused by the bacterium Mycobacterium avium subspecies paratuberculosis (MAP). Johne's disease is highly contagious and MAP infection in dairy cattle can eventually lead to death. With no available treatment for Johne's disease, genetic selection and improvements in management practices could help reduce its prevalence. In a previous study, the gene coding interleukin-10 receptor subunit alpha (IL10Rα) was associated with Johne's disease in dairy cattle. Our objective was to determine how IL10Rα affects the pathogenesis of MAP by examining the effect of a live MAP challenge on a mammary epithelial cell line (MAC-T) that had IL10Rα knocked out using CRISPR/cas9. The wild type and the IL10Rα knockout MAC-T cell lines were exposed to live MAP bacteria for 72 h. Thereafter, mRNA was extracted from infected and uninfected cells. Differentially expressed genes were compared between the wild type and the IL10Rα knockout cell lines. Gene ontology was performed based on the differentially expressed genes to determine which biological pathways were involved.

RESULTS

Immune system processes pathways were targeted to determine the effect of IL10Rα on the response to MAP infection. There was a difference in immune response between the wild type and IL10Rα knockout MAC-T cell lines, and less difference in immune response between infected and not infected IL10Rα knockout MAC-T cells, indicating IL10Rα plays an important role in the progression of MAP infection. Additionally, these comparisons allowed us to identify other genes involved in inflammation-mediated chemokine and cytokine signalling, interleukin signalling and toll-like receptor pathways.

CONCLUSIONS

Identifying differentially expressed genes in wild type and ILR10α knockout MAC-T cells infected with live MAP bacteria provided further evidence that IL10Rα contributes to mounting an immune response to MAP infection and allowed us to identify additional potential candidate genes involved in this process. We found there was a complex immune response during MAP infection that is controlled by many genes.

Identification of candidate genes associated with milk production and mastitis based on transcriptome-wide association study

Genetic research for the assessment of mastitis and milk production traits simultaneously has a long history. The main issue that arises in this context is the known existence of a positive correlation between the risk of mastitis and lactation performance due to selection. The transcriptome-wide association study (TWAS) approach endeavors to combine the expression quantitative trait loci and genome-wide association study summary statistics to decode complex traits or diseases. Accordingly, we used the farmgtex project results as a complete bovine database for mastitis and milk production. The results of colocalization and TWAS approaches were used for the detection of functional associated candidate genes with milk production and mastitis traits on multiple tissue-based transcriptome records. Also, we used the david database for gene ontology to identify significant terms and associated genes. For the identification of interaction networks, the genemania and string databases were used. Also, the available z-scores in TWAS results were used for the calculation of the correlation between tissues. Therefore, the present results confirm that LYNX1, DGAT1, C14H8orf33, and LY6E were identified as significant genes associated with milk production in eight, six, five, and five tissues, respectively. Also, FBXL6 was detected as a significant gene associated with mastitis trait. CLN3 and ZNF34 genes emerged via both the colocalization and TWAS approaches as significant genes for milk production trait. It is expected that TWAS and colocalization can improve our perception of the potential health status control mechanism in high-yielding dairy cows.

Characterization of the Hepatic Transcriptome for Divergent Immune-Responding Sheep Following Natural Exposure to Gastrointestinal Nematodes

Infections with gastrointestinal nematodes (GINs) reduce the economic efficiency of sheep operations and compromise animal welfare. Understanding the host's response to GIN infection can help producers identify animals that are naturally resistant to infection. The objective of this study was to characterize the hepatic transcriptome of sheep that had been naturally exposed to GIN parasites. The hepatic transcriptome was studied using RNA-Sequencing technology in animals characterized as high (n = 5) or medium (n = 6) based on their innate immune acute-phase (AP) response phenotype compared with uninfected controls (n = 4), and with biased antibody-mediated (AbMR, n = 5) or cell-mediated (CMR, n = 5) adaptive immune responsiveness compared to uninfected controls (n = 3). Following the assessment of sheep selected for innate responses, 0, 136, and 167 genes were differentially expressed (DE) between high- and medium-responding animals, high-responding and uninfected control animals, and medium-responding and uninfected control animals, respectively (false discovery rate (FDR) < 0.05, and fold change |FC| > 2). When adaptive immune responses were assessed, 0, 53, and 57 genes were DE between antibody- and cell-biased animals, antibody-biased and uninfected control animals, and cell-biased and uninfected control animals, respectively (FDR < 0.05, |FC| > 2). Functional analyses identified enriched gene ontology (GO) terms and metabolic pathways related to the innate immune response and energy metabolism. Six functional candidate genes were identified for further functional and validation studies to better understand the underlying biological mechanisms of host responses to GINs. These, in turn, can potentially help improve decision making and management practices to increase the overall host immune response to GIN infection.

Novel lncRNA regulatory elements in milk somatic cells of Holstein dairy cows associated with mastitis

Despite regulatory elements such as long non - coding RNAs representing most of the transcriptome, the functional understanding of long non - coding RNAs in relation to major health conditions including bovine mastitis is limited. This study examined the milk somatic cell transcriptome from udder quarters of 6 Holstein dairy cows to identify differentially expressed long non - coding RNAs using RNA - Sequencing. Ninety - four differentially expressed long non - coding RNAs are identified, 5 of which are previously annotated for gene name and length, 11 are annotated for gene name and 78 are novel, having no gene name or length previously annotated. Significant inflammatory response and regulation of immune response pathways (false discovery rate < 0.05) are associated with the differentially expressed long non - coding RNAs. QTL annotation analysis revealed 31 QTL previously annotated in the genomic regions of the 94 differentially expressed long non - coding RNAs, and the majority are associated with milk traits. This research provides a better understanding of long non - coding RNAs regulatory elements in milk somatic cells, which may enhance current breeding strategies for more adaptable or high mastitis resistant cattle.

Transmission ratio distortion detection by neutral genetic markers in the Pura Raza Española horse breed

Transmission Ratio Distortion (TRD) is a genetic phenomenon widely demonstrated in several livestock species, but barely in equine species. The TRD occurs when certain genotypes are over- or under-represented in the offspring of a particular mating and can be caused by a variety of factors during gamete formation or during embryonic development. For this study, 126 394 trios consisting of a stallion, mare, and offspring were genotyped using a panel of 17 neutral microsatellite markers recommended by the International Society for Animal Genetics for paternity tests and individual identification. The number of alleles available for each marker ranges from 13 to 18, been 268 the total number of alleles investigated. The TRDscan v.2.0 software was used with the biallelic procedure to identify regions with distorted segregation ratios. After completing the analysis, a total of 12 alleles (out of 11 microsatellites) were identified with decisive evidence for genotypic TRD; 3 and 9 with additive and heterosis patterns, respectively. In addition, 19 alleles (out of 10 microsatellites) were identified displaying allelic TRD. Among them, 14 and 5 were parent-unspecific and stallion-mare-specific TRD. Out of the TRD regions, 24 genes were identified and annotated, predominantly associated with cholesterol metabolism and homeostasis. These genes are often linked to non-specific symptoms like impaired fertility, stunted growth, and compromised overall health. The results suggest a significant impact on the inheritance of certain genetic traits in horses. Further analysis and validation are needed to better understand the TRD impact before the potential implementation in the horse breeding programme strategies.

Transcriptome profiling revealed that key rumen epithelium functions change in relation to short-chain fatty acids and rumen epithelium-attached microbiota during the weaning transition

This study aims to characterize the functional changes of the rumen epithelium associated with ruminal short-chain fatty acid (SCFA) concentration and epithelium-attached microbes during the weaning transition in dairy calves. Ruminal SCFA concentrations were determined, and transcriptome and microbiota profiling in biopsied rumen papillae were obtained from Holstein calves before and after weaning using RNA- and amplicon sequencing. Metabolic pathway analysis showed that pathways related to SCFA metabolism and cell apoptosis were up- and down-regulated postweaning, respectively. Functional analysis showed that genes related to SCFA absorption, metabolism, and protective roles against oxidative stress were positively correlated with ruminal SCFA concentrations. The relative abundance of epithelium-attached Rikenellaceae RC9 gut group and Campylobacter was positively correlated with genes involved in SCFA absorption and metabolism, suggesting that these microbes can cooperatively affect host functions. Future research should examine the contribution of attenuated apoptosis on rumen epithelial functional shifts during the weaning transition.

Reduction in mucosa thickness is associated with changes in immune function in the colon mucosa during the weaning transition in Holstein bull dairy calves

This study aims to characterize changes in the structure and the molecules related to immune function in the colon mucosa in dairy calves during the weaning transition (weaned at week 6 of age). Colon mucosa thickness, measured at week 5 to 8 and 12 of age, decreased for 2 weeks after weaning, but then recovered. Colon mucosa's transcriptome profiling at week 5, 7, and 12 of age was obtained using RNA-sequencing. Functional analysis showed that pathways related to immune function were up-regulated postweaning. A weighted gene co-expression network analysis identified 17 immune function related genes, expressed higher postweaning, which were negatively correlated with colon mucosa thickness, suggesting that these genes may be involved in colon mucosa inflammation and recovery from mucosa thickness decrement during the weaning transition. As such, it is important to determine the function of immune cells in the colon mucosa during the weaning transition in dairy calves.

Single-step genome-wide association for selected milk fatty acids in Dual-Purpose Belgian Blue cows

The aim of this study was to estimate genetic parameters and identify genomic regions associated with selected individual and groups of milk fatty acids (FA) predicted by milk mid-infrared spectrometry in Dual-Purpose Belgian Blue cows. The used data were 69,349 test-day records of milk yield, fat percentage, and protein percentage along with selected individual and groups FA of milk (g/dL milk) collected from 2007 to 2020 on 7,392 first-parity (40,903 test-day records), and 5,185 second-parity (28,446 test-day records) cows distributed in 104 herds in the Walloon Region of Belgium. Data of 28,466 SNPs, located on 29 Bos taurus autosomes (BTA), of 1,699 animals (639 males and 1,060 females) were used. Random regression test-day models were used to estimate genetic parameters through the Bayesian Gibbs sampling method. The SNP solutions were estimated using a single-step genomic best linear unbiased prediction approach. The proportion of genetic variance explained by each 25-SNP sliding window (with an average size of ~2 Mb) was calculated, and regions accounting for at least 1.0% of the total additive genetic variance were used to search for candidate genes. Average daily heritability estimated for the included milk FA traits ranged from 0.01 (C4:0) to 0.48 (C12:0) and 0.01 (C4:0) to 0.42 (C12:0) in the first and second parities, respectively. Genetic correlations found between milk yield and the studied individual milk FA, except for C18:0, C18:1 trans, C18:1 cis-9, were positive. The results showed that fat percentage and protein percentage were positively genetically correlated with all studied individual milk FA. Genome-wide association analyses identified 11 genomic regions distributed over 8 chromosomes [BTA1, BTA4, BTA10, BTA14 (4 regions), BTA19, BTA22, BTA24, and BTA26] associated with the studied FA traits, though those found on BTA14 partly overlapped. The genomic regions identified differed between parities and lactation stages. Although these differences in genomic regions detected may be due to the power of quantitative trait locus detection, it also suggests that candidate genes underlie the phenotypic expression of the studied traits may vary between parities and lactation stages. These findings increase our understanding about the genetic background of milk FA and can be used for the future implementation of genomic evaluation to improve milk FA profile in Dual-Purpose Belgian Blue cows.

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.

Effects of Serotonin on Cell Viability, Permeability of Bovine Mammary Gland Epithelial Cells and Their Transcriptome Analysis

Serotonin (5-HT) has been reported to play an important role in mammary gland involution that is defined as the process through which the gland returns to a nonlactating state. However, the overall picture of the regulatory mechanisms of 5-HT and the effects of serotonylation on mammary gland involution still need to be further investigated. The current study aimed to investigate the effects of 5-HT on global gene expression profiles of bovine mammary epithelial cells (MAC-T) and to preliminarily examine whether the serotonylation involved in the mammary gland involution by using Monodansylcadaverine (MDC), a competitive inhibitor of transglutaminase 2. Results showed that a high concentration of 5-HT decreased viability and transepithelial electrical resistance (TEER) in MAC-T cells. Transcriptome analysis indicated that 2477 genes were differentially expressed in MAC-T cells treated with 200 μg/mL of 5-HT compared with the control group, and the Notch, p53, and PI3K-Akt signaling pathways were enriched. MDC influenced 5-HT-induced MAC-T cell death, fatty acid synthesis, and the formation and disruption of tight junctions. Overall, a high concentration of 5-HT is able to accelerate mammary gland involution, which may be regulated through the Notch, p53, and PI3K-Akt signaling pathways. Serotonylation is involved in bovine mammary gland involution.

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.

Feed efficiency in dairy sheep: An insight from the milk transcriptome

Introduction

As higher feed efficiency in dairy ruminants means a higher capability to transform feed nutrients into milk and milk components, differences in feed efficiency are expected to be partly linked to changes in the physiology of the mammary glands. Therefore, this study aimed to determine the biological functions and key regulatory genes associated with feed efficiency in dairy sheep using the milk somatic cell transcriptome.

Material and methods

RNA-Seq data from high (H-FE, n = 8) and low (L-FE, n = 8) feed efficiency ewes were compared through differential expression analysis (DEA) and sparse Partial Least Square-Discriminant analysis (sPLS-DA).

Results

In the DEA, 79 genes were identified as differentially expressed between both conditions, while the sPLS-DA identified 261 predictive genes [variable importance in projection (VIP) &gt; 2] that discriminated H-FE and L-FE sheep.

Discussion

The DEA between sheep with divergent feed efficiency allowed the identification of genes associated with the immune system and stress in L-FE animals. In addition, the sPLS-DA approach revealed the importance of genes involved in cell division (e.g., KIF4A and PRC1) and cellular lipid metabolic process (e.g., LPL, SCD, GPAM, and ACOX3) for the H-FE sheep in the lactating mammary gland transcriptome. A set of discriminant genes, commonly identified by the two statistical approaches, was also detected, including some involved in cell proliferation (e.g., SESN2, KIF20A, or TOP2A) or encoding heat-shock proteins (HSPB1). These results provide novel insights into the biological basis of feed efficiency in dairy sheep, highlighting the informative potential of the mammary gland transcriptome as a target tissue and revealing the usefulness of combining univariate and multivariate analysis approaches to elucidate the molecular mechanisms controlling complex traits.

Transcriptome Profile in Dairy Cows Resistant or Sensitive to Milk Fat Depression

Feeding linseed to dairy cows results in milk fat depression (MFD), but there is a wide range of sensitivity among cows. The objectives of this study were to identify target genes containing SNP that may play a key role in the regulation of milk fat synthesis in cows resistant or sensitive to MFD. Four cows were selected from a dairy farm after a switch from a control diet to a linseed-rich diet; two were resistant to MFD with a high milk fat content in the control (4.06%) and linseed-rich (3.90%) diets; and two were sensitive to MFD with the milk fat content decreasing after the change from the control (3.87%) to linseed-rich (2.52%) diets. Transcriptome and SNP discovery analyses were performed using RNA-sequencing technology. There was a large number of differentially expressed genes in the control (n = 1316) and linseed-rich (n = 1888) diets. Of these, 15 genes were detected as key gene regulators and harboring SNP in the linseed-rich diet. The selected genes MTOR, PDPK1, EREG, NOTCH1, ZNF217 and TGFB3 may form a network with a principal axis PI3K/Akt/MTOR/SREBP1 involved in milk fat synthesis and in the response to diets that induced MFD. These 15 genes are novel candidate genes to be involved in the resistance or sensitivity of dairy cows to milk fat depression.

Identification of Novel mRNA Isoforms Associated with Acute Heat Stress Response Using RNA Sequencing Data in Sprague Dawley Rats

The molecular mechanisms underlying heat stress tolerance in animals to high temperatures remain unclear. This study identified the differentially expressed mRNA isoforms which narrowed down the most reliable DEG markers and molecular pathways that underlie the mechanisms of thermoregulation. This experiment was performed on Sprague Dawley rats housed at 22 °C (control group; CT), and three acute heat-stressed groups housed at 42 °C for 30 min (H30), 60 min (H60), and 120 min (H120). Earlier, we demonstrated that acute heat stress increased the rectal temperature of rats, caused abnormal changes in the blood biochemical parameters, as well as induced dramatic changes in the expression levels of genes through epigenetics and post-transcriptional regulation. Transcriptomic analysis using RNA-Sequencing (RNA-Seq) data obtained previously from blood (CT and H120), liver (CT, H30, H60, and H120), and adrenal glands (CT, H30, H60, and H120) was performed. The differentially expressed mRNA isoforms (DEIs) were identified and annotated by the CLC Genomics Workbench. Biological process and metabolic pathway analyses were performed using Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. A total of 225, 5764, and 4988 DEIs in the blood, liver, and adrenal glands were observed. Furthermore, the number of novel differentially expressed transcript lengths with annotated genes and novel differentially expressed transcript with non-annotated genes were 136 and 8 in blood, 3549 and 120 in the liver, as well as 3078 and 220 in adrenal glands, respectively. About 35 genes were involved in the heat stress response, out of which, Dnaja1, LOC680121, Chordc1, AABR07011951.1, Hsp90aa1, Hspa1b, Cdkn1a, Hmox1, Bag3, and Dnaja4 were commonly identified in the liver and adrenal glands, suggesting that these genes may regulate heat stress response through interactions between the liver and adrenal glands. In conclusion, this study would enhance our understanding of the complex underlying mechanisms of acute heat stress, and the identified mRNA isoforms and genes can be used as potential candidates for thermotolerance selection in mammals.

Immune mechanisms, resistance genes, and their roles in the prevention of mastitis in dairy cows

Mastitis is one of the most important diseases of the mammary gland. The increased incidence of this disease in cows is due to the breeding of dairy cattle for higher yields, which is accompanied by an increased susceptibility to mastitis. Therefore, the difficulty involved with preventing this disease has increased. An integral part of current research is the elimination of mastitis in order to reduce the consumption of antibiotic drugs, thereby reducing the resistance of microorganisms and decreasing companies' economic losses due to mastitis (i.e. decreased milk yield, increased drug costs, and reduced milk supply). Susceptibility to mastitis is based on dairy cows' immunity, health, nutrition, and welfare. Thus, it is important to understand the immune processes in the body in order to increase the resistance of animals. Recently, various studies have focused on the selection of mastitis resistance genes. An important point is also the prevention of mastitis. This publication aims to describe the physiology of the mammary gland along with its immune mechanisms and to approximate their connection with potential mastitis resistance genes. This work describes various options for mastitis elimination and focuses on genetic selection and a closer specification of resistance genes to mastitis. Among the most promising resistance genes for mastitis, we consider CD14, CXCR1, lactoferrin, and lactoglobulin.

The Analysis of the Ubiquitylomic Responses to Streptococcus agalactiae Infection in Bovine Mammary Gland Epithelial Cells

Purpose

Streptococcus agalactiae is one of the primary pathogens responsible for subclinical mastitis, a significant economic burden for dairy farms. An essential component of the immune response to infection is ubiquitination, which plays important roles in the complex interactions between the pathogen and host.

Materials and Methods

In the present study, quantitative ubiquitylomics was performed to profile changes in the global ubiquitinome of bovine mammary gland epithelial cells (BMECs) infected with S. agalactiae.

Results

The most notable changes in the BMEC ubiquitinome were related to the adherens junction, ribosome, and tight junction pathways. Ubiquitination of CTNNB1, EGFR, ITGB1, CTNNA1, CTNNA2, CDH1, YES1, and SLC9A3R1 appears to be fundamental for regulating multiple cellular processes in BMECs in response to S. agalactiae infection. In addition, broad ubiquitination of various effectors and outer membrane proteins was observed. Ubiquitinated proteins in S. agalactiae-infected BMECs were associated with regulating cell junctions in the host, with potential implications for susceptibility to infection.

Conclusion

The preliminary findings suggest that extensive ubiquitination of CTNNB1, CDH1 and SLC9A3R1 and proteins closely related to cell junctions might play an important role in mastitis progression in dairy cows. The results provide evidence that ubiquitin modification of certain proteins in S. agalactiae-infected BMECs could be a promising therapeutic strategy for reducing mammary gland injury and mastitis.

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.

Pathological insights into camel mastitis

Camel is a multipurpose animal bred to produce milk, meat, and transport and serves as a financial reserve for pastoralists by playing an important role in social prestige and prosperity. Camel milk is a good substitute for human milk because of its exceptional nutritional properties. Udder infections are considered one of the main limitations to camel farming. In recent decades, the disease has been reported by numerous camel-producing countries in Africa and Asia, such as Egypt, Somalia, Sudan, Kenya, Saudi Arabia, and Iraq. The current review provides an overview of the forms of camel mastitis, which can be clinical mastitis characterized by hardening and swelling of the breast, pain on palpation, and visible changes in the colour and texture of the milk or subclinical mastitis refers to the presence of inflammation with no obvious signs and it can be detected by indirect tests such as the California mastitis test (CMT), somatic cell count (SCC), and microbiological examination. Major pathogens of camel mastitis are Staphylococcus aureus, Streptococcus agalactiae, Escherichia coli, and Corynebacterium bovis. Regarding the risk factors for camel mastitis, this study provides an overview of the most important risk factors such as severe tick infestation, teat injuries, hygienic milking protocols, and physiological disorders causing mastitis. The use of indirect tests and bacteriological studies as diagnostic tools and their values for detecting camel mastitis will also be reviewed. Based on the above, further epidemiological studies on camel mastitis are needed to have solid scientific data on disease transmission, pathogen characterization, other possible risk factors or diagnostic methods, and the impact of the disease on public health. Proper control strategies should be adopted through early diagnosis, treatment and by avoiding potential risk factors to get good quality milk from camels.

Identification of novel alternative splicing associated with mastitis disease in Holstein dairy cows using large gap read mapping

Background

Mastitis is a very common disease in the dairy industry that producers encounter daily. Transcriptomics, using RNA-Sequencing (RNA - Seq) technology, can be used to study the functional aspect of mastitis resistance to identify animals that have a better immune response to mastitis. When the cow has mastitis, not only genes but also specific mRNA isoforms generated via alternative splicing (AS) could be differentially expressed (DE), leading to the phenotypic variation observed. Therefore, the objective of this study was to use large gap read mapping to identify mRNA isoforms DE between healthy and mastitic milk somatic cell samples (N = 12). These mRNA isoforms were then categorized based on being 1) annotated mRNA isoforms for gene name and length, 2) annotated mRNA isoforms with different transcript length and 3) novel mRNA isoforms of non - annotated genes.

Results

Analysis identified 333 DE transcripts (with at least 2 mRNA isoforms annotated, with at least one being DE) between healthy and mastitic samples corresponding to 303 unique genes. Of these 333 DE transcripts between healthy and mastitic samples, 68 mRNA isoforms are annotated in the bovine genome reference (ARS.UCD.1.2), 249 mRNA isoforms had novel transcript lengths of known genes and 16 were novel transcript lengths of non - annotated genes in the bovine genome reference (ARS.UCD.1.2). Functional analysis including gene ontology, gene network and metabolic pathway analysis was performed on the list of 288 annotated and unique DE mRNA isoforms. In total, 67 significant metabolic pathways were identified including positive regulation of cytokine secretion and immune response. Additionally, numerous DE novel mRNA isoforms showed potential involvement with the immune system or mastitis. Lastly, QTL annotation analysis was performed on coding regions of the DE mRNA isoforms, identifying overlapping QTLs associated with clinical mastitis and somatic cell score.

Conclusion

This study identified novel mRNA isoforms generated via AS that could lead to differences in the immune response of Holstein dairy cows and be potentially implemented in future breeding programs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08430-x.

Multi-Omics Integration and Network Analysis Reveal Potential Hub Genes and Genetic Mechanisms Regulating Bovine Mastitis

Mastitis, inflammation of the mammary gland, is the most prevalent disease in dairy cattle that has a potential impact on profitability and animal welfare. Specifically designed multi-omics studies can be used to prioritize candidate genes and identify biomarkers and the molecular mechanisms underlying mastitis in dairy cattle. Hence, the present study aimed to explore the genetic basis of bovine mastitis by integrating microarray and RNA-Seq data containing healthy and mastitic samples in comparative transcriptome analysis with the results of published genome-wide association studies (GWAS) using a literature mining approach. The integration of different information sources resulted in the identification of 33 common and relevant genes associated with bovine mastitis. Among these, seven genes—CXCR1, HCK, IL1RN, MMP9, S100A9, GRO1, and SOCS3—were identified as the hub genes (highly connected genes) for mastitis susceptibility and resistance, and were subjected to protein-protein interaction (PPI) network and gene regulatory network construction. Gene ontology annotation and enrichment analysis revealed 23, 7, and 4 GO terms related to mastitis in the biological process, molecular function, and cellular component categories, respectively. Moreover, the main metabolic-signalling pathways responsible for the regulation of immune or inflammatory responses were significantly enriched in cytokine–cytokine-receptor interaction, the IL-17 signaling pathway, viral protein interaction with cytokines and cytokine receptors, and the chemokine signaling pathway. Consequently, the identification of these genes, pathways, and their respective functions could contribute to a better understanding of the genetics and mechanisms regulating mastitis and can be considered a starting point for future studies on bovine mastitis.

Genome-wide association study for meat tenderness in beef cattle identifies patterns of the genetic contribution in different post-mortem stages

The beef tenderization process during the post-mortem period is one of the most important sensorial attributes and it is well-established. The aim of this study was to identify the genetic contribution pattern to meat tenderness at 7-(LMD7), 14-(LMD14), and 21-(LMD21) days post-mortem. The heritabilities for LMD7 (0.194), LMD14 (0.142) and LMD21 (0.048) are well established in the population evaluated here. However, its genetic contribution in terms of genomic candidate regions is still poorly understood. Tenderness was measured in the Longissiums thoracis using Warner-Bratzler shear force in the three post-mortem periods. A total of 4323 crossbred beef cattle were phenotyped and genotyped using the Illumina BovineSNP50K. The percentage of the total genetic variance was estimated using the weighted single-step genomic best linear unbiased prediction method. The main candidate windows for LMD7 were associated with proteolysis of myofibrillar structures and the weakening endomysium and perimysium. Candidate windows for LMD14 and LMD21 were mapped in bovine QTLs for body composition, height and growth. Results presented herein highlight, the largest contribution of proteolysis related processes before 14-days post-mortem and body composition characteristics in later stages for meat tenderness.

Transcriptome sequencing analysis for the identification of stable lncRNAs associated with bovine Staphylococcus aureus mastitis

Background

Staphylococcus aureus (S. aureus) mastitis is one of the most difficult diseases to treat in lactating dairy cows worldwide. S. aureus with different lineages leads to different host immune responses. Long non-coding RNAs (lncRNAs) are reported to be widely involved in the progress of inflammation. However, no research has identified stable lncRNAs among different S. aureus strain infections. In addition, folic acid (FA) can effectively reduce inflammation, and whether the inflammatory response caused by S. aureus can be reduced by FA remains to be explored.

Methods

lncRNA transcripts were identified from Holstein mammary gland tissues infected with different concentrations of S. aureus (in vivo) and mammary alveolar cells (Mac-T cells, in vitro) challenged with different S. aureus strains. Differentially expressed (DE) lncRNAs were evaluated, and stable DE lncRNAs were identified in vivo and in vitro. On the basis of the gene sequence conservation and function conservation across species, key lncRNAs with the function of potentially immune regulation were retained for further analysis. The function of FA on inflammation induced by S. aureus challenge was also investigated. Then, the association analysis between these keys lncRNA transcripts and hematological parameters (HPs) was carried out. Lastly, the knockdown and overexpression of the important lncRNA were performed to validate the gene function on the regulation of cell immune response.

Results

Linear regression analysis showed a significant correlation between the expression levels of lncRNA shared by mammary tissue and Mac-T cells (P < 0.001, R² = 0.3517). lncRNAs PRANCR and TNK2–AS1 could be regarded as stable markers associated with bovine S. aureus mastitis. Several HPs could be influenced by SNPs around lncRNAs PRANCR and TNK2–AS1. The results of gene function validation showed PRANCR regulates the mRNA expression of SELPLG and ITGB2 within the S. aureus infection pathway and the Mac-T cells apoptosis. In addition, FA regulated the expression change of DE lncRNA involved in toxin metabolism and inflammation to fight against S. aureus infection.

Conclusions

The remarkable association between SNPs around these two lncRNAs and partial HP indicates the potentially important role of PRANCR and TNK2–AS1 in immune regulation. Stable DE lncRNAs PRANCR and TNK2–AS1 can be regarded as potential targets for the prevention of bovine S. aureus mastitis. FA supplementation can reduce the negative effect of S. aureus challenge by regulating the expression of lncRNAs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40104-021-00639-2.

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

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

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.

Integrative Systems Biology Analysis Elucidates Mastitis Disease Underlying Functional Modules in Dairy Cattle

Background: Mastitis is the most prevalent disease in dairy cattle and one of the most significant bovine pathologies affecting milk production, animal health, and reproduction. In addition, mastitis is the most common, expensive, and contagious infection in the dairy industry. Methods: A meta-analysis of microarray and RNA-seq data was conducted to identify candidate genes and functional modules associated with mastitis disease. The results were then applied to systems biology analysis via weighted gene coexpression network analysis (WGCNA), Gene Ontology, enrichment analysis for the Kyoto Encyclopedia of Genes and Genomes (KEGG), and modeling using machine-learning algorithms. Results: Microarray and RNA-seq datasets were generated for 2,089 and 2,794 meta-genes, respectively. Between microarray and RNA-seq datasets, a total of 360 meta-genes were found that were significantly enriched as "peroxisome," "NOD-like receptor signaling pathway," "IL-17 signaling pathway," and "TNF signaling pathway" KEGG pathways. The turquoise module (n = 214 genes) and the brown module (n = 57 genes) were identified as critical functional modules associated with mastitis through WGCNA. PRDX5, RAB5C, ACTN4, SLC25A16, MAPK6, CD53, NCKAP1L, ARHGEF2, COL9A1, and PTPRC genes were detected as hub genes in identified functional modules. Finally, using attribute weighting and machine-learning methods, hub genes that are sufficiently informative in Escherichia coli mastitis were used to optimize predictive models. The constructed model proposed the optimal approach for the meta-genes and validated several high-ranked genes as biomarkers for E. coli mastitis using the decision tree (DT) method. Conclusion: The candidate genes and pathways proposed in this study may shed new light on the underlying molecular mechanisms of mastitis disease and suggest new approaches for diagnosing and treating E. coli mastitis in dairy cattle.

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.

Comprehensive RNA-Seq Profiling Reveals Temporal and Tissue-Specific Changes in Gene Expression in Sprague-Dawley Rats as Response to Heat Stress Challenges

Understanding heat stress physiology and identifying reliable biomarkers are paramount for developing effective management and mitigation strategies. However, little is known about the molecular mechanisms underlying thermal tolerance in animals. In an experimental model of Sprague–Dawley rats subjected to temperatures of 22 ± 1°C (control group; CT) and 42°C for 30 min (H30), 60 min (H60), and 120 min (H120), RNA-sequencing (RNA-Seq) assays were performed for blood (CT and H120), liver (CT, H30, H60, and H120), and adrenal glands (CT, H30, H60, and H120). A total of 53, 1,310, and 1,501 differentially expressed genes (DEGs) were significantly identified in the blood (P < 0.05 and |fold change (FC)| >2), liver (P < 0.01, false discovery rate (FDR)–adjusted P = 0.05 and |FC| >2) and adrenal glands (P < 0.01, FDR-adjusted P = 0.05 and |FC| >2), respectively. Of these, four DEGs, namely Junb, P4ha1, Chordc1, and RT1-Bb, were shared among the three tissues in CT vs. H120 comparison. Functional enrichment analyses of the DEGs identified in the blood (CT vs. H120) revealed 12 biological processes (BPs) and 25 metabolic pathways significantly enriched (FDR = 0.05). In the liver, 133 BPs and three metabolic pathways were significantly detected by comparing CT vs. H30, H60, and H120. Furthermore, 237 BPs were significantly (FDR = 0.05) enriched in the adrenal glands, and no shared metabolic pathways were detected among the different heat-stressed groups of rats. Five and four expression patterns (P < 0.05) were uncovered by 73 and 91 shared DEGs in the liver and adrenal glands, respectively, over the different comparisons. Among these, 69 and 73 genes, respectively, were proposed as candidates for regulating heat stress response in rats. Finally, together with genome-wide association study (GWAS) results in cattle and phenome-wide association studies (PheWAS) analysis in humans, five genes (Slco1b2, Clu, Arntl, Fads1, and Npas2) were considered as being associated with heat stress response across mammal species. The datasets and findings of this study will contribute to a better understanding of heat stress response in mammals and to the development of effective approaches to mitigate heat stress response in livestock through breeding.

Identification of functional candidate variants and genes for feed efficiency in Holstein and Jersey cattle breeds using RNA-sequencing

The identification of functional genetic variants and associated candidate genes linked to feed efficiency may help improve selection for feed efficiency in dairy cattle, providing economic and environmental benefits for the dairy industry. This study used RNA-sequencing data obtained from liver tissue from 9 Holstein cows [n = 5 low residual feed intake (RFI), n = 4 high RFI] and 10 Jersey cows (n = 5 low RFI, n = 5 high RFI), which were selected from a single population of 200 animals. Using RNA-sequencing, 3 analyses were performed to identify: (1) variants within low or high RFI Holstein cattle; (2) variants within low or high RFI Jersey cattle; and (3) variants within low or high RFI groups, which are common across both Holstein and Jersey cattle breeds. From each analysis, all variants were filtered for moderate, modifier, or high functional effect, and co-localized quantitative trait loci (QTL) classes, enriched biological processes, and co-localized genes related to these variants, were identified. The overlapping of the resulting genes co-localized with functional SNP from each analysis in both breeds for low or high RFI groups were compared. For the first two analyses, the total number of candidate genes associated with moderate, modifier, or high functional effect variants fixed within low or high RFI groups were 2,810 and 3,390 for Holstein and Jersey breeds, respectively. The major QTL classes co-localized with these variants included milk and reproduction QTL for the Holstein breed, and milk, production, and reproduction QTL for the Jersey breed. For the third analysis, the common variants across both Holstein and Jersey breeds, uniquely fixed within low or high RFI groups were identified, revealing a total of 86,209 and 111,126 functional variants in low and high RFI groups, respectively. Across all 3 analyses for low and high RFI cattle, 12 and 31 co-localized genes were overlapping, respectively. Among the overlapping genes across breeds, 9 were commonly detected in both the low and high RFI groups (INSRR, CSK, DYNC1H1, GAB1, KAT2B, RXRA, SHC1, TRRAP, PIK3CB), which are known to play a key role in the regulation of biological processes that have high metabolic demand and are related to cell growth and regeneration, metabolism, and immune function. The genes identified and their associated functional variants may serve as candidate genetic markers and can be implemented into breeding programs to help improve the selection for feed efficiency in dairy cattle.

Genome-wide association study to identify genomic regions and positional candidate genes associated with male fertility in beef cattle

Fertility plays a key role in the success of calf production, but there is evidence that reproductive efficiency in beef cattle has decreased during the past half-century worldwide. Therefore, identifying animals with superior fertility could significantly impact cow-calf production efficiency. The objective of this research was to identify candidate regions affecting bull fertility in beef cattle and positional candidate genes annotated within these regions. A GWAS using a weighted single-step genomic BLUP approach was performed on 265 crossbred beef bulls to identify markers associated with scrotal circumference (SC) and sperm motility (SM). Eight windows containing 32 positional candidate genes and five windows containing 28 positional candidate genes explained more than 1% of the genetic variance for SC and SM, respectively. These windows were selected to perform gene annotation, QTL enrichment, and functional analyses. Functional candidate gene prioritization analysis revealed 14 prioritized candidate genes for SC of which MAP3K1 and VIP were previously found to play roles in male fertility. A different set of 14 prioritized genes were identified for SM and five were previously identified as regulators of male fertility (SOD2, TCP1, PACRG, SPEF2, PRLR). Significant enrichment results were identified for fertility and body conformation QTLs within the candidate windows. Gene ontology enrichment analysis including biological processes, molecular functions, and cellular components revealed significant GO terms associated with male fertility. The identification of these regions contributes to a better understanding of fertility associated traits and facilitates the discovery of positional candidate genes for future investigation of causal mutations and their implications.

Relationships between novel nucleotide variants within the colony-stimulating factor 1 receptor (CSF1R) gene and mastitis indicators in sheep

Colony-stimulating factor 1 receptor (CSF1R) plays an important role in the process of innate immunity and inflammation, thus it was hypothesized that the CSF1R gene might affect the occurrence of mammalian mastitis. The purpose of this study was to investigate the association between nucleotide variations of CSF1R gene and mastitis in Australian white sheep (AUWs). Two indel variants (Intron5-27 bp and Intron5-22 bp) within the CSF1R gene have been found in AUWs. The Chi-square test for different mastitis symptoms demonstrated that individuals without symptoms of mastitis had higher 'I' allele frequencies and 'II' genotype frequencies (p < 0.01). We found strong correlation between mastitis and lactation score through Pearson correlation analysis. Therefore, we also analyzed the relationship between the two indel loci and lactation, we found that the lactation ability of individuals with type II was stronger than that of DD genotype at the Intron5-22 bp (p < 0.05). Additionally, we found that the combined genotype of the two loci was significantly associated with mastitis (p < 0.01). These findings indicated that CSF1R mutations were significantly associated with mastitis, and could affect lactation performance, suggesting that two deletion sites could be used as the effective molecular markers against mastitis in sheep breeding.