DNA methylation haplotype block signatures responding to Staphylococcus aureus subclinical mastitis and association with production and health traits.

BACKGROUND: DNA methylation has been documented to play vital roles in diseases and biological processes. In bovine, little is known about the regulatory roles of DNA methylation alterations on production and health traits, including mastitis. RESULTS: Here, we employed whole-genome DNA methylation sequencing to profile the DNA methylation patterns of milk somatic cells from sixteen cows with naturally occurring Staphylococcus aureus (S. aureus) subclinical mastitis and ten healthy control cows. We observed abundant DNA methylation alterations, including 3,356,456 differentially methylated cytosines and 153,783 differential methylation haplotype blocks (dMHBs). The DNA methylation in regulatory regions, including promoters, first exons and first introns, showed global significant negative correlations with gene expression status. We identified 6435 dMHBs located in the regulatory regions of differentially expressed genes and significantly correlated with their corresponding genes, revealing their potential effects on transcriptional activities. Genes harboring DNA methylation alterations were significantly enriched in multiple immune- and disease-related pathways, suggesting the involvement of DNA methylation in regulating host responses to S. aureus subclinical mastitis. In addition, we found nine discriminant signatures (differentiates cows with S. aureus subclinical mastitis from healthy cows) representing the majority of the DNA methylation variations related to S. aureus subclinical mastitis. Validation of seven dMHBs in 200 cows indicated significant associations with mammary gland health (SCC and SCS) and milk production performance (milk yield). CONCLUSIONS: In conclusion, our findings revealed abundant DNA methylation alterations in milk somatic cells that may be involved in regulating mammary gland defense against S. aureus infection. Particularly noteworthy is the identification of seven dMHBs showing significant associations with mammary gland health, underscoring their potential as promising epigenetic biomarkers. Overall, our findings on DNA methylation alterations offer novel insights into the regulatory mechanisms of bovine subclinical mastitis, providing further avenues for the development of effective control measures.

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

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

Genome-Wide DNA Methylation and Transcriptome Integration Associates DNA Methylation Changes with Bovine Subclinical Mastitis Caused by Staphylococcus chromogenes.

Staphylococcus chromogenes (SC) is a common coagulase-negative staphylococcus described as an emerging mastitis pathogen and commonly found in dairy farms. This study investigated the potential involvement of DNA methylation in subclinical mastitis caused by SC. The whole-genome DNA methylation patterns and transcriptome profiles of milk somatic cells from four cows with naturally occurring SC subclinical mastitis (SCM) and four healthy cows were characterized by next-generation sequencing, bioinformatics, and integration analyses. Comparisons revealed abundant DNA methylation changes related to SCM, including differentially methylated cytosine sites (DMCs, n = 2,163,976), regions (DMRs, n = 58,965), and methylation haplotype blocks (dMHBs, n = 53,098). Integration of methylome and transcriptome data indicated a negative global association between DNA methylation at regulatory regions (promoters, first exons, and first introns) and gene expression. A total of 1486 genes with significant changes in the methylation levels of their regulatory regions and corresponding gene expression showed significant enrichment in biological processes and pathways related to immune functions. Sixteen dMHBs were identified as candidate discriminant signatures, and validation of two signatures in more samples further revealed the association of dMHBs with mammary gland health and production. This study demonstrated abundant DNA methylation changes with possible involvement in regulating host responses and potential as biomarkers for SCM.

Genome-wide association analysis identified both RNA-seq and DNA variants associated to paratuberculosis in Canadian Holstein cattle 'in vitro' experimentally infected macrophages.

BACKGROUND: Mycobacterium avium ssp. paratuberculosis (MAP) is the causative agent of paratuberculosis, or Johne's disease (JD), an incurable bovine disease. The evidence for susceptibility to MAP disease points to multiple interacting factors, including the genetic predisposition to a dysregulation of the immune system. The endemic situation in cattle populations can be in part explained by a genetic susceptibility to MAP infection. In order to identify the best genetic improvement strategy that will lead to a significant reduction of JD in the population, we need to understand the link between genetic variability and the biological systems that MAP targets in its assault to dominate macrophages. MAP survives in macrophages where it disseminates. We used next-generation RNA (RNA-Seq) sequencing to study of the transcriptome in response to MAP infection of the macrophages from cows that have been naturally infected and identified as positive for JD (JD (+); n = 22) or negative for JD (healthy/resistant, JD (-); n = 28). In addition to identifying genetic variants from RNA-seq data, SNP variants were also identified using the Bovine SNP50 DNA chip. RESULTS: The complementary strategy allowed the identification of 1,356,248 genetic variants, including 814,168 RNA-seq and 591,220 DNA chip variants. Annotation using SnpEff predicted that the 2435 RNA-seq genetic variants would produce high functional effect on known genes in comparison to the 33 DNA chip variants. Significant variants from JD(+/-) macrophages were identified by genome-wide association study and revealed two quantitative traits loci: BTA4 and 11 at (P < 5 × 10- 7). Using BovineMine, gene expression levels together with significant genomic variants revealed pathways that potentially influence JD susceptibility, notably the energy-dependent regulation of mTOR by LKB1-AMPK and the metabolism of lipids. CONCLUSION: In the present study, we succeeded in identifying genetic variants in regulatory pathways of the macrophages that may affect the susceptibility of cows that are healthy/resistant to MAP infection. RNA-seq provides an unprecedented opportunity to investigate gene expression and to link the genetic variations to biological pathways that MAP normally manipulate during the process of killing macrophages. A strategy incorporating functional markers into genetic selection may have a considerable impact in improving resistance to an incurable disease. Integrating the findings of this research into the conventional genetic selection program may allow faster and more lasting improvement in resistance to bovine paratuberculosis in dairy cattle.

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.

Association of genetic polymorphisms related to Johne's disease with estimated breeding values of Holstein sires for milk ELISA test scores.

BACKGROUND: Johne's disease (JD) is a chronic intestinal inflammatory disease caused by Mycobacterium avium subsp. paratuberculosis (MAP) infection in ruminants. Since there are currently no effective vaccine or treatment options available to control JD, genetic selection may be an alternative strategy to enhance JD resistance. Numerous Single Nucleotide Polymorphisms (SNPs) have been reported to be associated with MAP infection status based on published genome-wide association and candidate gene studies. The main objective of this study was to validate these SNPs that were previously identified to be associated with JD by testing their effect on Holstein bulls' estimated breeding values (EBVs) for milk ELISA test scores, an indirect indicator of MAP infection status in cattle. RESULTS: Three SNPs, rs41810662, rs41617133 and rs110225854, located on Bos taurus autosomes (BTA) 16, 23 and 26, respectively, were confirmed as significantly associated with Holstein bulls' EBVs for milk ELISA test score (FDR < 0.01) based on General Quasi Likelihood Scoring analysis (GQLS) analysis. Single-SNP regression analysis identified four SNPs that were associated with sire EBVs (FDR < 0.05). This includes two SNPs that were common with GQLS (rs41810662 and rs41617133), with the other two SNPs being rs110494981 and rs136182707, located on BTA9 and BTA16, respectively. CONCLUSIONS: The findings of this study validate the association of SNPs with JD MAP infection status and highlight the need to further investigate the genomic regions harboring these SNPs.

Short communication: Genetic association of variations in the osteopontin gene (SPP1) with lactation persistency in dairy cattle.

Improving lactation persistency (LP) in dairy cattle has a beneficial effect on animal health and fertility and herd productivity. A complex trait, LP not only reflects the cow's ability to maintain milk secretion activity after the lactation peak but is also a function of the postcalving development of the mammary gland and, later on, of tissue remodeling as lactation declines. This decline is a consequence of an imbalance between cell proliferation and cell removal. In a previous study, single nucleotide polymorphisms were identified in the osteopontin (OPN) gene, SPP1. Osteopontin is a multifaceted protein that plays an important role in immune regulation and tissue remodeling. Because OPN is involved in involution, it might also have an effect on LP. The objective of the present study was to evaluate whether LP could be influenced by genetic variations in the SPP1 gene. This association with LP was analyzed in the population of 578 bulls characterized in a previous study. The population mean of estimated breeding value (EBV) for LP was 100.95 ± 5.06 units. Allele and genotype association analyses were performed by comparing the frequencies of the different genotypes and alleles with EBV for LP for the respective lactation using logistic regression. The EBV for LP at the first lactation (LP1), second lactation (LP2), and third lactation (LP3) and for overall lactation (OLP) are reported for the genotypes SPP1c.-1301G>A, SPP1c.-1251C>T, SPP1c.-430G>A, and SPP1c.*40A>C. The first single nucleotide polymorphism, SPP1c.-1301G>A, affected LP1, LP2, LP3, and OLP. Analysis of the estimated average allele substitution effects also confirmed that G is a favorable allele for LP, given the gain observed over LP1, LP2, LP3, and OLP. Differences in EBV for LP were observed between animals with different haplotypes at LP1, LP2, LP3, and OLP. Contrast analysis for OLP revealed that mean EBV is greater for block H1 (101.34 ± 0.30) than for animals that do not have H1 (98.20 ± 0.77). The gain with block H1 (GCGA) suggests the presence of the favorable allele G (first position in the block: SPP1c.-1301G). The pleiotropic roles of OPN position it at the crossroads of immune regulation, tissue remodeling, and involution. From a genetic perspective, data from the present study suggest OPN as a candidate gene associated with LP for dairy cows.

Transcriptome Analysis of Long Non-Coding RNA in the Bovine Mammary Gland Following Dietary Supplementation with Linseed Oil and Safflower Oil.

This study aimed to characterize the long non-coding RNA (lncRNA) expression in the bovine mammary gland and to infer their functions in dietary response to 5% linseed oil (LSO) or 5% safflower oil (SFO). Twelve cows (six per treatment) in mid lactation were fed a control diet for 28 days followed by a treatment period (control diet supplemented with 5% LSO or 5% SFO) of 28 days. Mammary gland biopsies were collected from each animal on day-14 (D-14, control period), D+7 (early treatment period) and D+28 (late treatment period) and were subjected to RNA-Sequencing and subsequent bioinformatics analyses. Functional enrichment of lncRNA was performed via potential cis regulated target genes located within 50 kb flanking regions of lncRNAs and having expression correlation of >0.7 with mRNAs. A total of 4955 lncRNAs (325 known and 4630 novel) were identified which potentially cis targeted 59 and 494 genes in LSO and SFO treatments, respectively. Enrichments of cis target genes of lncRNAs indicated potential roles of lncRNAs in immune function, nucleic acid metabolism and cell membrane organization processes as well as involvement in Notch, cAMP and TGF-ß signaling pathways. Thirty-two and 21 lncRNAs were differentially expressed (DE) in LSO and SFO treatments, respectively. Six genes (KCNF1, STARD13, BCL6, NXPE2, HHIPL2 and MMD) were identified as potential cis target genes of six DE lncRNAs. In conclusion, this study has identified lncRNAs with potential roles in mammary gland functions and potential candidate genes and pathways via which lncRNAs might function in response to LSO and SFA.

Co-Expression Network Analysis Identifies miRNA⁻mRNA Networks Potentially Regulating Milk Traits and Blood Metabolites.

MicroRNAs (miRNA) regulate mRNA networks to coordinate cellular functions. In this study, we constructed gene co-expression networks to detect miRNA modules (clusters of miRNAs with similar expression patterns) and miRNA⁻mRNA pairs associated with blood (triacylglyceride and nonesterified fatty acids) and milk (milk yield, fat, protein, and lactose) components and milk fatty acid traits following dietary supplementation of cows' diets with 5% linseed oil (LSO) (n = 6 cows) or 5% safflower oil (SFO) (n = 6 cows) for 28 days. Using miRNA transcriptome data from mammary tissues of cows for co-expression network analysis, we identified three consensus modules: blue, brown, and turquoise, composed of 70, 34, and 86 miRNA members, respectively. The hub miRNAs (miRNAs with the most connections with other miRNAs) were miR-30d, miR-484 and miR-16b for blue, brown, and turquoise modules, respectively. Cell cycle arrest, and p53 signaling and transforming growth factor⁻beta (TGF-ß) signaling pathways were the common gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways enriched for target genes of the three modules. Protein percent (p = 0.03) correlated with the turquoise module in LSO treatment while protein yield (p = 0.003) and milk yield (p = 7 × 10-04) correlated with the turquoise model, protein and milk yields and lactose percent (p < 0.05) correlated with the blue module and fat percent (p = 0.04) correlated with the brown module in SFO treatment. Several fatty acids correlated (p < 0.05) with the blue (CLA:9,11) and brown (C4:0, C12:0, C22:0, C18:1n9c and CLA:10,12) modules in LSO treatment and with the turquoise (C14:0, C18:3n3 and CLA:9,11), blue (C14:0 and C23:0) and brown (C6:0, C16:0, C22:0, C22:6n3 and CLA:10,12) modules in SFO treatment. Correlation of miRNA and mRNA data from the same animals identified the following miRNA⁻mRNA pairs: miR-183/RHBDD2 (p = 0.003), miR-484/EIF1AD (p = 0.011) and miR-130a/SBSPON (p = 0.004) with lowest p-values for the blue, brown, and turquoise modules, respectively. Milk yield, protein yield, and protein percentage correlated (p < 0.05) with 28, 31 and 5 miRNA⁻mRNA pairs, respectively. Our results suggest that, the blue, brown, and turquoise modules miRNAs, hub miRNAs, miRNA⁻mRNA networks, cell cycle arrest GO term, p53 signaling and TGF-ß signaling pathways have considerable influence on milk and blood phenotypes following dietary supplementation of dairy cows' diets with 5% LSO or 5% SFO.

Low-depth genotyping-by-sequencing (GBS) in a bovine population: strategies to maximize the selection of high quality genotypes and the accuracy of imputation.

BACKGROUND: Genotyping-by-sequencing (GBS) has emerged as a powerful and cost-effective approach for discovering and genotyping single-nucleotide polymorphisms. The GBS technique was largely used in crop species where its low sequence coverage is not a drawback for calling genotypes because inbred lines are almost homozygous. In contrast, only a few studies used the GBS technique in animal populations (with sizeable heterozygosity rates) and many of those that have been published did not consider the quality of the genotypes produced by the bioinformatic pipelines. To improve the sequence coverage of the fragments, an alternative GBS preparation protocol that includes selective primers during the PCR amplification step has been recently proposed. In this study, we compared this modified protocol with the conventional two-enzyme GBS protocol. We also described various procedures to maximize the selection of high quality genotypes and to increase the accuracy of imputation. RESULTS: The in silico digestions of the bovine genome showed that the combination of PstI and MspI is more suitable for sequencing bovine GBS libraries than the use of single digestions with PstI or ApeKI. The sequencing output of the GBS libraries generated a total of 123,666 variants with the selective-primer approach and 272,103 variants with the conventional approach. Validating our data with genotypes obtained from mass spectrometry and Illumina's bovine SNP50 array, we found that the genotypes produced by the conventional GBS method were concordant with those produced by these alternative genotyping methods, whereas the selective-primer method failed to call heterozygotes with confidence. Our results indicate that high accuracy in genotype calling (>97%) can be obtained using low read-depth thresholds (3 to 5 reads) provided that markers are simultaneously filtered for genotype quality scores. We also show that factors such as the minimum call rate and the minor allele frequency positively influence the accuracy of imputation of missing GBS data. The highest accuracies (around 85%) of imputed GBS markers were obtained with the FIMPUTE program when GBS and SNP50 array genotypes were combined (80,190 to 100,297 markers) before imputation. CONCLUSIONS: We discovered that the conventional two-enzyme GBS protocol could produce a large number of high-quality genotypes provided that appropriate filtration criteria were used. In contrast, the selective-primer approach resulted in a substantial proportion of miscalled genotypes and should be avoided for livestock genotyping studies. Overall, our study demonstrates that carefully adjusting the different filtering parameters applied to the GBS data is critical to maximize the selection of high quality genotypes and to increase the accuracy of imputation of missing data. The strategies and results presented here provide a framework to maximize the output of the GBS technique in animal populations and qualified the PstI/MspI GBS assay as a low-cost high-density genotyping platform. The conclusions reported here regarding read-depth and genotype quality filtering could benefit many GBS applications, notably genome-wide association studies, where there is a need to increase the density of markers genotyped across the target population while preserving the quality of genotypes.

Supplements of vitamins B9 and B12 affect hepatic and mammary gland gene expression profiles in lactating dairy cows.

BACKGROUND: A combined supplement of vitamins B9 and B12 was reported to increase milk and milk component yields of dairy cows without effect on feed intake. The present study was undertaken to verify whether this supplementation positively modifies the pathways involved in milk and milk component synthesis. Thus, by studying the transcriptome activity in these tissues, the effect of supplements of both vitamins on the metabolism of both liver and mammary gland, was investigated. For this study, 24 multiparous Holstein dairy cows were assigned to 6 blocks of 4 animals each according to previous 305-day milk production. Within each block, cows were randomly assigned to weekly intramuscular injections of 5 mL of either saline 0.9 % NaCl, 320 mg of vitamin B9, 10 mg of vitamin B12 or a combination of both vitamins (B9 + B12). The experimental period began 3 weeks before the expected calving date and lasted 9 weeks of lactation. Liver and mammary biopsies were performed on lactating dairy cows 64 ± 3 days after calving. Samples from both tissues were analyzed by microarray and qPCR to identify genes differentially expressed in hepatic and mammary tissues. RESULTS: Microarray analysis identified 47 genes in hepatic tissue and 16 genes in the mammary gland whose expression was modified by the vitamin supplements. Gene ontology (GO) categorizes genes in non-overlapping domains of molecular biology. Panther is one of the online GO resources used for gene function classification. It classifies the 63 genes according to Molecular Function, Biological Process and Protein Class. Most of the biological processes modulated by the vitamin supplements were associated to developmental process, protein metabolic process, transport and response to inflammation. In the liver, most of the genes modulated by the vitamin treatments involved protein metabolic process while developmental process appeared to be more affected by the treatments in mammary gland. Out of 25 genes analysed by qPCR, 7 were validated. CONCLUSION: The results indicate that several metabolic processes were modulated by the supplementation of vitamins in early-lactating dairy cows. In addition, the results suggest that the vitamin supplements promoted liver regeneration and reduced catabolism of lipids in early lactation.

Transcriptome adaptation of the bovine mammary gland to diets rich in unsaturated fatty acids shows greater impact of linseed oil over safflower oil on gene expression and metabolic pathways.

BACKGROUND: Nutritional strategies can decrease saturated fatty acids (SFAs) and increase health beneficial fatty acids (FAs) in bovine milk. The pathways/genes involved in these processes are not properly defined. Next-generation RNA-sequencing was used to investigate the bovine mammary gland transcriptome following supplemental feeding with 5% linseed oil (LSO) or 5% safflower oil (SFO). Holstein cows in mid-lactation were fed a control diet for 28 days (control period) followed by supplementation with 5% LSO (12 cows) or 5% SFO (12 cows) for 28 days (treatment period). Milk and mammary gland biopsies were sampled on days-14 (control period), +7 and +28 (treatment period). Milk was used to measure fat(FP)/protein(PP) percentages and individual FAs while RNA was subjected to sequencing. RESULTS: Milk FP was decreased by 30.38% (LSO) or 32.42% (SFO) while PP was unaffected (LSO) or increased (SFO). Several beneficial FAs were increased by LSO (C18:1n11t, CLA:10t12c, CLA:9c11t, C20:3n3, C20:5n3, C22:5n3) and SFO (C18:1n11t, CLA:10t12c, C20:1c11, C20:2, C20:3n3) while several SFAs (C4:0, C6:0, C8:0, C14:0, C16:0, C17:0, C24:0) were decreased by both treatments (P < 0.05). 1006 (460 up- and 546 down-regulated) and 199 (127 up- and 72 down-regulated) genes were significantly differentially regulated (DE) by LSO and SFO, respectively. Top regulated genes (≥ 2 fold change) by both treatments (FBP2, UCP2, TIEG2, ANGPTL4, ALDH1L2) are potential candidate genes for milk fat traits. Involvement of SCP2, PDK4, NQO1, F2RL1, DBI, CPT1A, CNTFR, CALB1, ACADVL, SPTLC3, PIK3CG, PIGZ, ADORA2B, TRIB3, HPGD, IGFBP2 and TXN in FA/lipid metabolism in dairy cows is being reported for the first time. Functional analysis indicated similar and different top enriched functions for DE genes. DE genes were predicted to significantly decrease synthesis of FA/lipid by both treatments and FA metabolism by LSO. Top canonical pathways associated with DE genes of both treatments might be involved in lipid/cholesterol metabolism. CONCLUSION: This study shows that rich α-linolenic acid LSO has a greater impact on mammary gland transcriptome by affecting more genes, pathways and processes as compared to SFO, rich in linoleic acid. Our study suggest that decrease in milk SFAs was due to down-regulation of genes in the FA/lipid synthesis and lipid metabolism pathways while increase in PUFAs was due to increased availability of ruminal biohydrogenation metabolites that were up taken and incorporated into milk or used as substrate for the synthesis of PUFAs.

Deep sequencing shows microRNA involvement in bovine mammary gland adaptation to diets supplemented with linseed oil or safflower oil.

BACKGROUND: Bovine milk fat composition is responsive to dietary manipulation providing an avenue to modify the content of fatty acids and especially some specific unsaturated fatty acid (USFA) isomers of benefit to human health. MicroRNAs (miRNAs) regulate gene expression but their specific roles in bovine mammary gland lipogenesis are unclear. The objective of this study was to determine the expression pattern of miRNAs following mammary gland adaptation to dietary supplementation with 5 % linseed or safflower oil using next generation RNA-sequencing. METHODS: Twenty-four Canadian Holstein dairy cows (twelve per treatment) in mid lactation were fed a control diet (total mixed ration of corn:grass silages) for 28 days followed by a treatment period (control diet supplemented with 5 % linseed or safflower oil) of 28 days. Milk samples were collected weekly for fat and individual fatty acid determination. RNA from mammary gland biopsies harvested on day-14 (control period) and on days +7 and +28 (treatment period) from six randomly selected cows per treatment was subjected to small RNA sequencing. RESULTS: Milk fat percentage decreased significantly (P < 0.001) during treatment with the two diets as compared to the control period. The individual saturated fatty acids C4:0, C6:0, C8:0, C14:0 and C16:0 decreased significantly (P < 0.05) while five USFAs (C14:1, C18:1n11t, C20:3n3, C20:5n3 and CLA:t10c12) increased remarkably (P < 0.05) in response to both treatments. Analysis of 361 million sequence reads generated 321 known bovine miRNAs and 176 novel miRNAs. The expression of fourteen and twenty-two miRNAs was affected (P < 0.05) by linseed and safflower oil treatments, respectively. Seven miRNAs including six up-regulated (bta-miR-199c, miR-199a-3p, miR-98, miR-378, miR-148b and miR-21-5p) and one down-regulated (bta-miR-200a) were found to be regulated (P < 0.05) by both treatments, and thus considered core differentially expressed (DE) miRNAs. The gene targets of core DE miRNAs have functions related to gene expression and general cellular metabolism (P < 0.05) and are enriched in four pathways of lipid metabolism (3-phosphoinositide biosynthesis, 3-phosphoinositide degradation, D-myo-inisitol-5-phosphate metabolism and the superpathway of inositol phosphate compounds). CONCLUSION: Our results suggest that DE miRNAs in this study might be important regulators of bovine mammary lipogenesis and metabolism. The novel miRNAs identified in this study will further enrich the bovine miRNome repertoire and contribute to understanding mammary gland biology.

Effects of sources and routes of administration of vitamins A, D and copper on postnatal status of these micronutrients in piglets.

Twenty-six nulliparous sows were fed conventional gestation and lactation diets supplemented (N = 13) or not (N = 13) with extra daily supplements of 25-hydroxy-cholecalciferol (25-OH-D3; 4 ĸIU), ß-carotene (24 ĸIU), and copper (Cu)-proteinate (45 mg) from day 90 of gestation to 21 d of lactation (L21). In each litter, 10 piglets were divided into 5 pairs received, at 2 (L2) and 8 d (L8) of age, one of the five combinations of micronutrient sources and routes of administration (N = 260 piglets total). These neonatal treatments (N = 26 pairs or 52 piglets each) consisted of oral vitamin D3, retinol acetate and CuSO4 (T1); oral 25-OH-D3, ß-carotene, and Cu proteinate (T2); exposure to ultraviolet light (UVB), oral retinol palmitate and Cu gluconate (T3); intramuscular vitamin D3 and retinyl propionate and oral Cu acetate (T4); oral saline (CTRL). Oral or intramuscular provisions corresponded to 12 mg of Cu and 70 and 12 ĸIU of vitamins A and D, respectively. Blood samples were collected from all piglets at L2, L8, and L21 for determination of serum Cu, retinol, and 25-OH-D3. Body weight was measured at birth, L2, L8, and L21. Piglets were weaned at L21, and liver and blood samples were collected 2 d later to evaluate oxidative enzymes in blood and liver and hepatic ATP concentrations and expression of genes associated with antioxidant status. Sow treatments had marginal or no impacts on Cu, retinol, 25-OH-D3, or antioxidant status in piglet blood serum and liver. However, when supplements were given to piglets, hepatic Cu was 38% greater in for all treated piglets compared to CTRL (P < 0.01), hepatic retinol was 3 times higher in T1 than in CTRL (P < 0.01) and intermediate for other treatments whereas serum 25-OH-D3 was markedly increased with T2 and T3 at L8 and L21, respectively, compared to CTRL (Piglet treatment × Age interaction, P < 0.01). Concerning antioxidant activities, glutathione peroxidase, and superoxide dismutase were increased (P < 0.03) in plasma of T2 piglets whereas the highest values (P < 0.03) for indicators of oxidative damage to proteins were observed in T4 piglets. The study revealed that oral Cu proteinate from T2, oral retinol acetate from T1, oral 25-hydroxy-cholecalciferol from T2, and UVB light exposure from T3 were the most efficient ways of increasing the postnatal status of these micronutrients in suckling piglets and this may have some impacts on their peri-weaning antioxidant status.

Investigating in vivo Mycobacterium avium subsp. paratuberculosis microevolution and mixed strain infections.

Mycobacterium avium subsp. paratuberculosis (MAP) causes Johne's Disease (JD) in ruminants, which is responsible for significant economic loss to the global dairy industry. Mixed strain infection (MSI) refers to the concurrent infection of a susceptible host with genetically distinct strains of a pathogen, whereas within-host changes in an infecting strain leading to genetically distinguishable progeny is called microevolution. The two processes can influence host-pathogen dynamics, disease progression and outcomes, but not much is known about their prevalence and impact on JD. Therefore, we obtained up to 10 MAP isolates each from 14 high-shedding animals and subjected them to whole-genome sequencing. Twelve of the 14 animals examined showed evidence for the presence of MSIs and microevolution, while the genotypes of MAP isolates from the remaining two animals could be attributed solely to microevolution. All MAP isolates that were otherwise isogenic had differences in short sequence repeats (SSRs), of which SSR1 and SSR2 were the most diverse and homoplastic. Variations in SSR1 and SSR2, which are located in ORF1 and ORF2, respectively, affect the genetic reading frame, leading to protein products with altered sequences and computed structures. The ORF1 gene product is predicted to be a MAP surface protein with possible roles in host immune modulation, but nothing could be inferred regarding the function of ORF2. Both genes are conserved in Mycobacterium avium complex members, but SSR1-based modulation of ORF1 reading frames seems to only occur in MAP, which could have potential implications on the infectivity of this pathogen. IMPORTANCE Johne's disease (JD) is a major problem in dairy animals, and concerns have been raised regarding the association of Mycobacterium avium subsp. paratuberculosis (MAP) with Crohn's disease in humans. MAP is an extremely slow-growing bacterium with low genome evolutionary rates. Certain short sequence repeats (SSR1 and SSR2) in the MAP chromosome are highly variable and evolve at a faster rate than the rest of the chromosome. In the current study, multiple MAP isolates with genetic variations such as single-nucleotide polymorphisms, and more noticeably, diverse SSRs, could simultaneously infect animals. Variations in SSR1 and SSR2 affect the products of the respective genes containing them. Since multiple MAP isolates can infect the same animal and the possibility that the pathogen undergoes further changes within the host due to unstable SSRs, this could provide a compensative mechanism for an otherwise slow-evolving pathogen to increase phenotypic diversity for overcoming host responses.

Genomic epidemiology of Mycobacterium avium subsp. paratuberculosis isolates from Canadian dairy herds provides evidence for multiple infection events.

Mycobacterium avium subsp. paratuberculosis (MAP) is the pathogen responsible for paratuberculosis or Johne's Disease (JD) in ruminants, which is responsible for substantial economic losses worldwide. MAP transmission primarily occurs through the fecal-oral route, and the introduction of an MAP infected animal into a herd is an important transmission route. In the current study, we characterized MAP isolates from 67 cows identified in 20 herds from the provinces of Quebec and Ontario, Canada. Whole genome sequencing (WGS) was performed and an average genome coverage (relative to K-10) of ∼14.9 fold was achieved. The total number of SNPs present in each isolate varied from 51 to 132 and differed significantly between herds. Isolates with the highest genetic variability were generally present in herds from Quebec. The isolates were broadly separated into two main clades and this distinction was not influenced by the province from which they originated. Analysis of 8 MIRU-VNTR loci and 11 SSR loci was performed on the 67 isolates from the 20 dairy herds and publicly available references, notably major genetic lineages and six isolates from the province of Newfoundland and Labrador. All 67 field isolates were phylogenetically classified as Type II (C-type) and according to MIRU-VNTR, the predominant type was INMV 2 (76.1%) among four distinct patterns. Multilocus SSR typing identified 49 distinct INMV SSR patterns. The discriminatory index of the multilocus SSR typing was 0.9846, which was much higher than MIRU-VNTR typing (0.3740). Although multilocus SSR analysis provides good discriminatory power, the resolution was not informative enough to determine inter-herd transmission. In select cases, SNP-based analysis was the only approach able to document disease transmission between herds, further validated by animal movement data. The presence of SNPs in several virulence genes, notably for PE, PPE, mce and mmpL, is expected to explain differential antigenic or pathogenetic host responses. SNP-based studies will provide insight into how MAP genetic variation may impact host-pathogen interactions. Our study highlights the informative power of WGS which is now recommended for epidemiological studies and to document mixed genotypes infections.

Weaning diet supplemented with health-promoting feed additives influences microbiota and immune response in piglets challenged with Salmonella.

The aim of this study was to evaluate the potential of micronutrients and feed additives to modulate intestinal microbiota and systemic and mucosal immune responses in weaned pigs infected with Salmonella. At weaning, 32 litters of 12 piglets each were allocated to four dietary treatments: 1) control diet (CTRL), 2) CTRL supplemented with chlortetracycline (ATB), 3) CTRL supplemented with a cocktail of feed additives (CKTL); and 4) CKTL diet containing bovine colostrum in replacement of spray-dry animal plasma (CKTL+COL). The CKTL supplement included cranberry extract, encapsulated carvacrol and yeast-derived products and an enriched selenium and vitamin premix. Three weeks after weaning, four pigs per litter were orally inoculated with Salmonella Typhimurium DT104. Half of them were euthanized 3 days post-infection (dpi) and the other half, 7 dpi. The expression of IL6, TNF, IL8, monocyte chemoattractant protein 1 (MCP1), IFNG, cyclooxygenase 2 (COX2), glutathione peroxidase 2 (GPX2) and ß-defensin 2 (DEFB2) showed a peaked response at 3 dpi (P < 0.05). Results also revealed that DEFB2 expression was higher at 3 dpi in CTRL and CKTL groups than in ATB (P = 0.01 and 0.06, respectively) while GPX2 gene was markedly increased at 3 and 7 dpi in pigs fed CKTL or CKTL+COL diet compared to CTRL pigs (P < 0.05). In piglets fed CKTL or CKTL+COL diet, intestinal changes in microbial communities were less pronounced after exposure to Salmonella compared to CTRL and progressed faster toward the status before Salmonella challenge (AMOVA P < 0.01). Furthermore, the relative abundance of several families was either up- or down-regulated in pigs fed CKTL or CKTL+COL diet after Salmonella challenge. In conclusion, weaning diet enriched with bovine colostrum, vitamins and mixture of feed additives mitigated the influence of Salmonella infection on intestinal microbial populations and modulate systemic and intestinal immune defences.

Variant Calling from RNA-seq Data Using the GATK Joint Genotyping Workflow.

The Genome Analysis Toolkit (GATK) developed at the Broad Institute provides state-of-the-art pipelines for germline and somatic variant discovery and genotyping. Unfortunately, the fully validated GATK pipeline for calling variant on RNAseq data is a Per-sample workflow that does not include the recent improvements seen in modern workflows, especially the possibility to perform joint genotyping analysis. Here, we describe how modern GATK commands from distinct workflows can be combined to call variants on RNAseq samples. We provide a detailed tutorial that starts with raw RNAseq reads and ends with filtered variants, of which some were shown to be associated with bovine paratuberculosis.

Methylome and transcriptome data integration reveals potential roles of DNA methylation and candidate biomarkers of cow Streptococcus uberis subclinical mastitis.

BACKGROUND: Mastitis caused by different pathogens including Streptococcus uberis (S. uberis) is responsible for huge economic losses to the dairy industry. In order to investigate the potential genetic and epigenetic regulatory mechanisms of subclinical mastitis due to S. uberis, the DNA methylome (whole genome DNA methylation sequencing) and transcriptome (RNA sequencing) of milk somatic cells from cows with naturally occurring S. uberis subclinical mastitis and healthy control cows (n = 3/group) were studied. RESULTS: Globally, the DNA methylation levels of CpG sites were low in the promoters and first exons but high in inner exons and introns. The DNA methylation levels at the promoter, first exon and first intron regions were negatively correlated with the expression level of genes at a whole-genome-wide scale. In general, DNA methylation level was lower in S. uberis-positive group (SUG) than in the control group (CTG). A total of 174,342 differentially methylated cytosines (DMCs) (FDR < 0.05) were identified between SUG and CTG, including 132,237, 7412 and 34,693 DMCs in the context of CpG, CHG and CHH (H = A or T or C), respectively. Besides, 101,612 methylation haplotype blocks (MHBs) were identified, including 451 MHBs that were significantly different (dMHB) between the two groups. A total of 2130 differentially expressed (DE) genes (1378 with up-regulated and 752 with down-regulated expression) were found in SUG. Integration of methylome and transcriptome data with MethGET program revealed 1623 genes with significant changes in their methylation levels and/or gene expression changes (MetGDE genes, MethGET P-value < 0.001). Functional enrichment of genes harboring ≥ 15 DMCs, DE genes and MetGDE genes suggest significant involvement of DNA methylation changes in the regulation of the host immune response to S. uberis infection, especially cytokine activities. Furthermore, discriminant correlation analysis with DIABLO method identified 26 candidate biomarkers, including 6 DE genes, 15 CpG-DMCs and 5 dMHBs that discriminated between SUG and CTG. CONCLUSION: The integration of methylome and transcriptome of milk somatic cells suggests the possible involvement of DNA methylation changes in the regulation of the host immune response to subclinical mastitis due to S. uberis. The presented genetic and epigenetic biomarkers could contribute to the design of management strategies of subclinical mastitis and breeding for mastitis resistance.

Administration of probiotics influences F4 (K88)-positive enterotoxigenic Escherichia coli attachment and intestinal cytokine expression in weaned pigs.

This study evaluated the effect of the probiotics Pediococcus acidilactici and Saccharomyces cerevisiae boulardii on the intestinal colonization of O149 enterotoxigenic Escherichia coli harbouring the F4 (K88) fimbriae (ETEC F4) and on the expression of ileal cytokines in weaned pigs. At birth, different litters of pigs were randomly assigned to one of the following treatments: 1) control without antibiotics or probiotics (CTRL); 2) reference group in which chlortetracycline and tiamulin were added to weanling feed (ATB); 3) P. acidilactici; 4) S. cerevisiae boulardii; or 5) P. acidilactici + S. cerevisiae boulardii. Probiotics were administered daily (1 × 10(9) CFU per pig) during the lactation period and after weaning (day 21). At 28 days of age, all pigs were orally challenged with an ETEC F4 strain, and a necropsy was performed 24 h later. Intestinal segments were collected to evaluate bacterial colonization in the small intestine and ileal cytokine expressions. Attachment of ETEC F4 to the intestinal mucosa was significantly reduced in pigs treated with P. acidilactici or S. cerevisiae boulardii in comparison with the ATB group (P = 0.01 and P = 0.03, respectively). In addition, proinflammatory cytokines, such as IL-6, were upregulated in ETEC F4 challenged pigs treated with P. acidilactici alone or in combination with S. cerevisiae boulardii compared with the CTRL group. In conclusion, the administration of P. acidilactici or S. cerevisiae boulardii was effective in reducing ETEC F4 attachment to the ileal mucosa, whereas the presence of P. acidilactici was required to modulate the expression of intestinal inflammatory cytokines in pigs challenged with ETEC F4.