Three Novel Players: PTK2B, SYK, and TNFRSF21 Were Identified to Be Involved in the Regulation of Bovine Mastitis Susceptibility via GWAS and Post-transcriptional Analysis

A Single-Step Genome-Wide Association Study for Semen Traits of Egyptian Buffalo Bulls

The present study aimed to contribute to the limited research on buffalo (Bubalus bubalis) semen traits by incorporating genomic data. A total of 8465 ejaculates were collected. The genotyping procedure was conducted using the Axiom® Buffalo Genotyping 90 K array designed by the Affymetrix Expert Design Program. After conducting a quality assessment, we utilized 67,282 SNPs genotyped in 192 animals. We identified several genomic loci explaining high genetic variance by employing single-step genomic evaluation. The aforementioned regions were located on buffalo chromosomes no. 3, 4, 6, 7, 14, 16, 20, 22, and the X-chromosome. The X-chromosome exhibited substantial influence, accounting for 4.18, 4.59, 5.16, 5.19, and 4.31% of the genomic variance for ejaculate volume, mass motility, livability, abnormality, and concentration, respectively. In the examined genomic regions, we identified five novel candidate genes linked to male fertility and spermatogenesis, four in the X-chromosome and one in chromosome no. 16. Additional extensive research with larger sample sizes and datasets is imperative to validate these findings and evaluate their applicability for genomic selection.

Retinoic acid ameliorates low-grade endotoxemia-induced mastitis by limiting inflammatory responses in mice

Mastitis is a serious disease for humans and animals, which causes huge economic losses in the dairy industry and is hard to prevent due to the complex and unclear pathogenesis. Subacute ruminal acidosis (SARA) has contributed to the development of mastitis by inducing ruminal dysbiosis and subsequent low-grade endotoxemia (LGE), however, how ruminal metabolic changes regulate this progress is still unclear. Our previous study revealed that cows with SARA had increased ruminal retinoic acid (RA) levels, a metabolic intermediate of vitamin A that plays an essential role in mucosal immune responses. Hence, the aim of this study was to investigate the protective effect of RA on LGE-induced mastitis and the underlying mechanisms in mice. The results showed that RA alleviated LGE-induced mastitis, as evidenced by RA significantly reduced the increase in mammary proinflammatory cytokines and improved blood-milk barrier injury caused by LGE. In addition, RA increased the expression of tight junction proteins, including ZO-1, occludin and claudin-3. Furthermore, we found that RA limited the mammary inflammatory responses by inhibiting the activation of NF-κB and NLRP3 signaling pathways. These findings suggest that RA effectively alleviates LGE-induced mastitis and implies a potential strategy for the treatment and prevention of mastitis and other diseases.

A Novel TLR4-SYK Interaction Axis Plays an Essential Role in the Innate Immunity Response in Bovine Mammary Epithelial Cells

Mammary gland epithelium, as the first line of defense for bovine mammary gland immunity, is crucial in the process of mammary glands’ innate immunity, especially that of bovine mammary epithelial cells (bMECs). Our previous studies successfully marked SYK as an important candidate gene for mastitis traits via GWAS and preliminarily confirmed that SYK expression is down-regulated in bMECs with LPS (E. coli) stimulation, but its work mechanism is still unclear. In this study, for the first time, in vivo, TLR4 and SYK were colocalized and had a high correlation in mastitis mammary epithelium; protein−protein interaction results also confirmed that there was a direct interaction between them in mastitis tissue, suggesting that SYK participates in the immune regulation of the TLR4 cascade for bovine mastitis. In vitro, TLR4 also interacts with SYK in LPS (E. coli)-stimulated or GBS (S. agalactiae)-infected bMECs, respectively. Moreover, TLR4 mRNA expression and protein levels were little affected in bMECsSYK- with LPS stimulation or GBS infection, indicating that SYK is an important downstream element of the TLR4 cascade in bMECs. Interestingly, IL-1β, IL-8, NF-κB and NLRP3 expression in LPS-stimulated or GBS-infected bMECsSYK- were significantly higher than in the control group, while AKT1 expression was down-regulated, implying that SYK could inhibit the IL-1β, IL-8, NF-κB and NLRP3 expression and alleviate inflammation in bMECs with LPS and GBS. Taken together, our solid evidence supports that TLR4/SYK/NF-κB signal axis in bMECs regulates the innate immunity response to LPS or GBS.

Emerging Roles of Noncoding RNAs in Bovine Mastitis Diseases

Non-coding RNAs (ncRNAs) are an abundant class of RNA with varying nucleotide lengths. They have been shown to have great potential in eutherians/human disease diagnosis and treatments and are now gaining more importance for the improvement of diseases in livestock. To date, thousands of ncRNAs have been discovered in the bovine genome and the continuous advancement in deep sequencing technologies and various bioinformatics tools has enabled the elucidation of their roles in bovine health. Among farm animals' diseases, mastitis, a common inflammatory disease in cattle, has caused devastating economic losses to dairy farmers over the last few decades. Here, we summarize the biology of bovine mastitis and comprehensively discuss the roles of ncRNAs in different types of mastitis infection. Based on our findings and relevant literature, we highlighted various evidence of ncRNA roles in mastitis. Different approaches (in vivo versus in vitro) for exploring ncRNA roles in mastitis are emphasized. More particularly, the potential applications of emerging genome editing technologies, as well as integrated omics platforms for ncRNA studies and implications for mastitis are presented.

Within- and between-Breed Selection Signatures in the Original and Improved Valachian Sheep

This study explored the genomic diversity and selection signatures in two Slovakian national breeds, the Original Valachian and the Improved Valachian sheep. As they are an important animal genetic resource within the country, but with decreasing population size, our aim is to identify potentially valuable genomic regions. A total of 97 sheep (18 male and 79 female) from the Original Valachian, and 69 sheep (25 male and 44 female) from the Improved Valachian populations were genotyped using the GeneSeek GGP Ovine 50 K chip. The inbreeding levels were assessed with runs of homozygosity (ROH). The selection signatures within breeds were identified based on the top 1% of most homozygous regions within the breed, the so-called ROH islands. The selection signatures between breeds were assessed based on variance in linkage disequilibrium. Overall, we have identified selection signatures with quantitative trait loci (QTL) and genes pointing towards all three production purposes of the Valachian sheep, milk, meat, and wool, including their quality characteristics. Another group with apparent large importance was the various traits related to health and resistance to parasites, which is well in line with the sturdy nature of this breed.

Genetic Parameter Estimation and Genome-Wide Association Study-Based Loci Identification of Milk-Related Traits in Chinese Holstein

Accurately estimating the genetic parameters and revealing more genetic variants underlying milk production and quality are conducive to the genetic improvement of dairy cows. In this study, we estimate the genetic parameters of five milk-related traits of cows-namely, milk yield (MY), milk fat percentage (MFP), milk fat yield (MFY), milk protein percentage (MPP), and milk protein yield (MPY)-based on a random regression test-day model. A total of 95,375 test-day records of 9,834 cows in the lower reaches of the Yangtze River were used for the estimation. In addition, genome-wide association studies (GWASs) for these traits were conducted, based on adjusted phenotypes. The heritability, as well as the standard errors, of MY, MFP, MFY, MPP, and MPY during lactation ranged from 0.22 ± 0.02 to 0.31 ± 0.04, 0.06 ± 0.02 to 0.15 ± 0.03, 0.09 ± 0.02 to 0.28 ± 0.04, 0.07 ± 0.01 to 0.16 ± 0.03, and 0.14 ± 0.02 to 0.27 ± 0.03, respectively, and the genetic correlations between different days in milk (DIM) within lactations decreased as the time interval increased. Two, six, four, six, and three single nucleotide polymorphisms (SNPs) were detected, which explained 5.44, 12.39, 8.89, 10.65, and 7.09% of the phenotypic variation in MY, MFP, MFY, MPP, and MPY, respectively. Ten Kyoto Encyclopedia of Genes and Genomes pathways and 25 Gene Ontology terms were enriched by analyzing the nearest genes and genes within 200 kb of the detected SNPs. Moreover, 17 genes in the enrichment results that may play roles in milk production and quality were selected as candidates, including CAMK2G, WNT3A, WNT9A, PLCB4, SMAD9, PLA2G4A, ARF1, OPLAH, MGST1, CLIP1, DGAT1, PRMT6, VPS28, HSF1, MAF1, TMEM98, and F7. We hope that this study will provide useful information for in-depth understanding of the genetic architecture of milk production and quality traits, as well as contribute to the genomic selection work of dairy cows in the lower reaches of the Yangtze River.

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.

Exploring the potential mechanism of Rhodomyrtus tomentosa (Ait.) Hassk fruit phenolic rich extract on ameliorating nonalcoholic fatty liver disease by integration of transcriptomics and metabolomics profiling

Nonalcoholic fatty liver disease (NAFLD), as the commonest form of chronic liver disease, is accompanied by liver oxidative stress and inflammatory responses. Rhodomyrtus tomentosa (Ait.) Hassk fruit phenolic rich extract (RTE) possesses multiple pharmacological effects in management of chronic diseases. In this study, the liver-protective effect of RTE on mice with high-fat-diet (HFD)-induced NAFLD was investigated for the first time, and the underlying molecular mechanism was explored via integration of transcriptomics and metabolomics. The results showed that RTE mitigated liver damage, which was evidenced by declined inflammatory cell infiltration in liver, decreased liver function markers, oxidative stress indexes, lipid profile levels and inflammatory cytokines levels. The differential metabolites by metabonomics illustrated supplementation of RTE affected metabolomics pathways including tryptophan metabolism, alanine, aspartate and glutamate metabolism, D-glutamine and D-glutamate metabolism, cysteine and methionine metabolism, arginine and proline metabolism, which are all involved in oxidative stress and inflammation. Furthermore, the five differential expression genes (DEGs) through liver transcriptomics were screened and recognized, namely Tnfrsf21, Ifit1, Inhbb, Mapk15 and Gadd45g, which revealed that HFD induced Cytokine-cytokine receptor interaction pathway, NF-κB signaling pathway NOD-like receptor pathway, TNF signaling pathway. Integrated analysis of transcriptomics and metabolomics confirmed the supplementation of RTE had significantly regulatory effects on the metabolic pathways involved in inflammatory responses. Additionally, RT-PCR and western blot authenticated RTE intervention regulated the mRNA levels of liver genes involved in inflammation response and inhibited the liver endotoxin-TLR4-NF-κB pathway triggered by HFD, thus alleviating NAFLD. Our findings strongly support the possibility that RTE can be regarded as a potential therapeutic method for obesity-associated NAFLD.

Identification of loci associated with susceptibility to Mycobacterium avium subsp. paratuberculosis infection in Holstein cattle using combinations of diagnostic tests and imputed whole-genome sequence data

Bovine paratuberculosis (PTB) is a chronic inflammatory disease caused by Mycobacterium avium susbp. paratuberculosis (MAP). Genome-wide association studies (GWAS) have identified single nucleotide polymorphisms (SNPs) significantly associated with susceptibility to bovine PTB. The main objective of this study was to identify quantitative trait loci (QTLs) associated with MAP infection in Spanish Holstein cows (N = 983) using combinations of diagnostic tests and imputed whole-genome sequence (WGS) data. The infection status of these animals was defined by three diagnostic methods including ELISA for MAP-antibodies detection, and tissue culture and PCR for MAP detection. The 983 cows included in this study were genotyped with the Bovine MD SNP50 Bead Chip, and the corresponding genotypes were imputed to WGS using the 1,000 Bull genomes reference population. In total, 33.77 million SNP variants per animal were identified across the genome. Linear mixed models were used to calculate the heritability (h²) estimates for each diagnostic test and test combinations. Next, we performed a case-control GWAS using the imputed WGS datasets and the phenotypes and combinations of phenotypes with h² estimates > 0.080. After performing the GWAS, the test combinations that showed SNPs with a significant association (P_FDR ≤ 0.05), were the ELISA-tissue PCR-tissue culture, ELISA-tissue culture, and ELISA-tissue PCR. A total of twelve quantitative trait loci (QTLs) highly associated with MAP infection status were identified on the Bos taurus autosomes (BTA) 4, BTA5, BTA11, BTA12, BTA14, BTA23, BTA24, and BTA28, and some of these QTLs were linked to immune-modulating genes. The identified QTLs on BTA23 spanning from 18.81 to 22.95 Mb of the Bos taurus genome overlapped with several QTLs previously found to be associated with PTB susceptibility, bovine tuberculosis susceptibility, and clinical mastitis. The results from this study provide more clues regarding the molecular mechanisms underlying susceptibility to PTB infection in cattle and might be used to develop national genetic evaluations for PTB in Spain.

Chitosan functionalized graphene oxide nanocomposites for fluorescence imaging of apoptotic processes and targeted anti-inflammation study

This work reports novel chitosan functionalized graphene oxide (GO) nanocomposites combined fluorescence imaging and therapeutic functions in one agent, which can serve as a promising alternative to alleviate related diseases caused hyperinflammation. Briefly, GO was designed to be conjugated with chitosan, fluorescein-labeled peptide, toll-like receptor 4 antibody and hydroxycamptothecin/aloe emodin. We have demonstrated that such nanocomposites could effectively achieve active targeted delivery of pro-apoptotic and anti-inflammatory drugs into inflammatory cells and cause cells apoptosis by acid-responsive drug release. Moreover, confocal fluorescence imaging confirms that the drug-induced inflammatory cells apoptosis could be visualized the light-up fluorescence of fluorescein activated by caspase-3. Meanwhile, inflammatory-related biomarkers have down-regulated after the nanocomposites' treatment in both vitro and vivo experiments consistent with the results in histological sections. In summary, the bifunctional nanocomposites that possess anti-inflammation and fluorescence imaging could serve as a promising therapeutic agent for reducing hyperinflammation caused by numerous diseases.

Genomic studies of milk-related traits in water buffalo (Bubalus bubalis) based on single-step genomic best linear unbiased prediction and random regression models

Genomic selection has been widely implemented in many livestock breeding programs, but it remains incipient in buffalo. Therefore, this study aimed to (1) estimate variance components incorporating genomic information in Murrah buffalo; (2) evaluate the performance of genomic prediction for milk-related traits using single- and multitrait random regression models (RRM) and the single-step genomic best linear unbiased prediction approach; and (3) estimate longitudinal SNP effects and candidate genes potentially associated with time-dependent variation in milk, fat, and protein yields, as well as somatic cell score (SCS) in multiple parities. The data used to estimate the genetic parameters consisted of a total of 323,140 test-day records. The average daily heritability estimates were moderate (0.35 ± 0.02 for milk yield, 0.22 ± 0.03 for fat yield, 0.42 ± 0.03 for protein yield, and 0.16 ± 0.03 for SCS). The highest heritability estimates, considering all traits studied, were observed between 20 and 280 d in milk (DIM). The genetic correlation estimates at different DIM among the evaluated traits ranged from -0.10 (156 to 185 DIM for SCS) to 0.61 (36 to 65 DIM for fat yield). In general, direct selection for any of the traits evaluated is expected to result in indirect genetic gains for milk yield, fat yield, and protein yield but also increase SCS at certain lactation stages, which is undesirable. The predicted RRM coefficients were used to derive the genomic estimated breeding values (GEBV) for each time point (from 5 to 305 DIM). In general, the tuning parameters evaluated when constructing the hybrid genomic relationship matrices had a small effect on the GEBV accuracy and a greater effect on the bias estimates. The SNP solutions were back-solved from the GEBV predicted from the Legendre random regression coefficients, which were then used to estimate the longitudinal SNP effects (from 5 to 305 DIM). The daily SNP effect for 3 different lactation stages were performed considering 3 different lactation stages for each trait and parity: from 5 to 70, from 71 to 150, and from 151 to 305 DIM. Important genomic regions related to the analyzed traits and parities that explain more than 0.50% of the total additive genetic variance were selected for further analyses of candidate genes. In general, similar potential candidate genes were found between traits, but our results suggest evidence of differential sets of candidate genes underlying the phenotypic expression of the traits across parities. These results contribute to a better understanding of the genetic architecture of milk production traits in dairy buffalo and reinforce the relevance of incorporating genomic information to genetically evaluate longitudinal traits in dairy buffalo. Furthermore, the candidate genes identified can be used as target genes in future functional genomics studies.

Identification of SNPs Associated with Somatic Cell Score in Candidate Genes in Italian Holstein Friesian Bulls

Mastitis is an infectious disease affecting the mammary gland, leading to inflammatory reactions and to heavy economic losses due to milk production decrease. One possible way to tackle the antimicrobial resistance issue stemming from antimicrobial therapy is to select animals with a genetic resistance to this disease. Therefore, aim of this study was to analyze the genetic variability of the SNPs found in candidate genes related to mastitis resistance in Holstein Friesian bulls. Target regions were amplified, sequenced by Next-Generation Sequencing technology on the Illumina® MiSeq, and then analyzed to find correlation with mastitis related phenotypes in 95 Italian Holstein bulls chosen with the aid of a selective genotyping approach. On a total of 557 detected mutations, 61 showed different genotype distribution in the tails of the deregressed EBVs for SCS and 15 were identified as significantly associated with the phenotype using two different approaches. The significant SNPs were identified in intergenic or intronic regions of six genes, known to be key components in the immune system (namely CXCR1, DCK, NOD2, MBL2, MBL1 and M-SAA3.2). These SNPs could be considered as candidates for a future genetic selection for mastitis resistance, although further studies are required to assess their presence in other dairy cattle breeds and their possible negative correlation with other traits.

Functional analysis of bovine interleukin-10 receptor alpha in response to Mycobacterium avium subsp. paratuberculosis lysate using CRISPR/Cas9

Background

The interleukin-10 receptor alpha (IL10RA) gene codes for the alpha chain of the IL-10 receptor which binds the cytokine IL-10. IL-10 is an anti-inflammatory cytokine with immunoregulatory function during the pathogenesis of many inflammatory disorders in livestock, including Johne’s disease (JD). JD is a chronic enteritis in cattle caused by Mycobacterium avium subsp. paratuberculosis (MAP) and is responsible for significant economic losses to the dairy industry. Several candidate genes including IL10RA have been found to be associated with JD. The aim of this study was to better understand the functional significance of IL10RA in the context of immune stimulation with MAP cell wall lysate.

Results

An IL10RA knock out (KO) bovine mammary epithelial cell (MAC-T) line was generated using the CRISPR/cas9 (Clustered Regularly Interspaced Short Palindromic Repeats/CRISPR-associated protein 9) gene editing system. These IL10RA KO cells were stimulated with the immune stimulant MAP lysate +/− IL-10, or with LPS as a positive control. In comparison to unedited cells, relative quantification of immune-related genes after stimulation revealed that knocking out IL10RA resulted in upregulation of pro-inflammatory cytokine gene expression (TNFA, IL1A, IL1B and IL6) and downregulation of suppressor of cytokine signaling 3 (SOCS3), a negative regulator of pro-inflammatory cytokine signaling. At the protein level knocking out IL10RA also resulted in upregulation of inflammatory cytokines - TNF-α and IL-6 and chemokines - IL-8, CCL2 and CCL4, relative to unedited cells.

Conclusions

The findings of this study illustrate the broad and significant effects of knocking out the IL10RA gene in enhancing pro-inflammatory cytokine expression and further support the immunoregulatory role of IL10RA in eliciting an anti-inflammatory response as well as its potential functional involvement during the immune response associated with JD.

A Fixed Cohort Field Study of Gene Expression in Circulating Leukocytes From Dairy Cows With and Without Mastitis

Specifically designed gene expression studies can be used to prioritize candidate genes and identify novel biomarkers affecting resilience against mastitis and other diseases in dairy cattle. The primary goal of this study was to assess whether specific peripheral leukocyte genes expressed differentially in a previous study of dairy cattle with postpartum disease, also would be expressed differentially in peripheral leukocytes from a diverse set of different dairy cattle with moderate to severe clinical mastitis. Four genes were selected for this study due to their differential expression in a previous transcriptomic analysis of circulating leukocytes from dairy cows with and without evidence of early postpartum disease. An additional 15 genes were included based on their cellular, immunologic, and inflammatory functions associated with resistance and tolerance to mastitis. This fixed cohort study was conducted on a conventional dairy in Washington state. Cows >50 days in milk (DIM) with mastitis (n = 12) were enrolled along with healthy cows (n = 8) selected to match the DIM and lactation numbers of mastitic cows. Blood was collected for a complete blood count (CBC), serum biochemistry, leukocyte isolation, and RNA extraction on the day of enrollment and twice more at 6 to 8-days intervals. Latent class analysis was performed to discriminate healthy vs. mastitic cows and to describe disease resolution. RNA samples were processed by the Primate Diagnostic Services Laboratory (University of Washington, Seattle, WA). Gene expression analysis was performed using the Nanostring System (Nanostring Technologies, Seattle, Washington, USA). Of the four genes (C5AR1, CATHL6, LCN2, and PGLYRP1) with evidence of upregulation in cows with mastitis, three of those genes (CATHL6, LCN2, and PGLYRP1) were investigated due to their previously identified association with postpartum disease. These genes are responsible for immunomodulatory molecules that selectively enhance or alter host innate immune defense mechanisms and modulate pathogen-induced inflammatory responses. Although further research is warranted to explain their functional mechanisms and bioactivity in cattle, our findings suggest that these conserved elements of innate immunity have the potential to bridge disease states and target tissues in diverse dairy populations.

Comparative Transcriptomic Analysis of the Pituitary Gland between Cattle Breeds Differing in Growth: Yunling Cattle and Leiqiong Cattle

The hypothalamic-pituitary-thyroid (HPT) axis hormones regulate the growth and development of ruminants, and the pituitary gland plays a decisive role in this process. In order to identify pivotal genes in the pituitary gland that could affect the growth of cattle by regulating the secretion of hormones, we detected the content of six HPT hormones related to growth in the plasma of two cattle breeds (Yunling and Leiqiong cattle, both also known as the zebu cattle) with great differences in growth and compared the transcriptome data of their pituitary glands. Our study found that the contents of GH, IGF, TSH, thyroxine, triiodothyronine, and insulin were significantly different between the two breeds, which was the main cause of the difference in growth; 175 genes were identified as differentially expressed genes (DEGs). Functional association analyses revealed that DEGs were mainly involved in the process of transcription and signal transduction. Combining the enrichment analysis and protein interaction analysis, eight DEGs were predicted to control the growth of cattle by affecting the expression of growth-related hormones in the pituitary gland. In summary, our results suggested that SLC38A1, SLC38A3, DGKH, GNB4, GNAQ, ESR1, NPY, and GAL are candidates in the pituitary gland for regulating the growth of Yunling and Leiqiong cattle by regulating the secretion of growth-related hormones. This study may help researchers further understand the growth mechanisms and improve the artificial selection of zebu cattle.