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.

Multi-omics integration identifies regulatory factors underlying bovine subclinical mastitis

BACKGROUND

Mastitis caused by multiple factors remains one of the most common and costly disease of the dairy industry. Multi-omics approaches enable the comprehensive investigation of the complex interactions between multiple layers of information to provide a more holistic view of disease pathogenesis. Therefore, this study investigated the genomic and epigenomic signatures and the possible regulatory mechanisms underlying subclinical mastitis by integrating RNA sequencing data (mRNA and lncRNA), small RNA sequencing data (miRNA) and DNA methylation sequencing data of milk somatic cells from 10 healthy cows and 20 cows with naturally occurring subclinical mastitis caused by Staphylococcus aureus or Staphylococcus chromogenes.

RESULTS

Functional investigation of the data sets through gene set analysis uncovered 3458 biological process GO terms and 170 KEGG pathways with altered activities during subclinical mastitis, provided further insights into subclinical mastitis and revealed the involvement of multi-omics signatures in the altered immune responses and impaired mammary gland productivity during subclinical mastitis. The abundant genomic and epigenomic signatures with significant alterations related to subclinical mastitis were observed, including 30,846, 2552, 1276 and 57 differential methylation haplotype blocks (dMHBs), differentially expressed genes (DEGs), lncRNAs (DELs) and miRNAs (DEMs), respectively. Next, 5 factors presenting the principal variation of differential multi-omics signatures were identified. The important roles of Factor 1 (DEG, DEM and DEL) and Factor 2 (dMHB and DEM), in the regulation of immune defense and impaired mammary gland functions during subclinical mastitis were revealed. Each of the omics within Factors 1 and 2 explained about 20% of the source of variation in subclinical mastitis. Also, networks of important functional gene sets with the involvement of multi-omics signatures were demonstrated, which contributed to a comprehensive view of the possible regulatory mechanisms underlying subclinical mastitis. Furthermore, multi-omics integration enabled the association of the epigenomic regulatory factors (dMHBs, DELs and DEMs) of altered genes in important pathways, such as 'Staphylococcus aureus infection pathway' and 'natural killer cell mediated cytotoxicity pathway', etc., which provides further insights into mastitis regulatory mechanisms. Moreover, few multi-omics signatures (14 dMHBs, 25 DEGs, 18 DELs and 5 DEMs) were identified as candidate discriminant signatures with capacity of distinguishing subclinical mastitis cows from healthy cows.

CONCLUSION

The integration of genomic and epigenomic data by multi-omics approaches in this study provided a better understanding of the molecular mechanisms underlying subclinical mastitis and identified multi-omics candidate discriminant signatures for subclinical mastitis, which may ultimately lead to the development of more effective mastitis control and management strategies.

Comparative study of HA and HNB staining RT-LAMP assays for peste des petits ruminants virus detection in West African Dwarf goats

Peste des petits ruminants (PPR) cause severe economic losses to many countries of the world where the disease is endemic. It has been targeted for global eradication by 2030 following the successful eradication of rinderpest in 2011. The proposed eradication program would benefit from efficient and relatively reliable diagnostic tools for early PPR virus (PPRV) detection. A total of 33 eight to 12 months old West African Dwarf (WAD) goats were used. Nineteen goats infected by commingling with two PPR virus-positive animals formed the infected group (PPRV-infected goats) while 14 non-infected goats formed the control group (CTG). The suitability of hydroxyl naphthol blue (HNB) staining of reverse transcription loop-mediated isothermal amplification (RT-LAMP) and haemagglutination (HA) assays was compared for their sensitivity to detect the PPRV in PPRV-infected goats and non-infected CTG. PPR disease severity in WAD goats at different days post infection (dpi) was evaluated by clinical scoring and haemagglutination titre (HAT). HNB staining RT-LAMP reaction and HA showed sensitivities of 100% and 73.68%, respectively, for PPRV detection. Expression of PPR clinical signs began from 3 dpi, attained peak at 5 dpi, thereafter showed irregular patterns till 24 dpi. Evaluation of HAT in PPRV-infected goats at 12 dpi ranged from 2 to 64 haemagglutination units (HAU), while CTG goats had 0 HAU. In conclusion, HA could be a good tool for rapid diagnosis of PPRV in a developing country setting. However, HNB staining RT-LAMP assay demonstrated high sensitivity for accurate diagnoses of PPRV and as an important diagnostic tool when precise phenotyping is desired.

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.

A sheep pangenome reveals the spectrum of structural variations and their effects on tail phenotypes

Structural variations (SVs) are a major contributor to genetic diversity and phenotypic variations, but their prevalence and functions in domestic animals are largely unexplored. Here we generated high-quality genome assemblies for 15 individuals from genetically diverse sheep breeds using Pacific Biosciences (PacBio) high-fidelity sequencing, discovering 130.3 Mb nonreference sequences, from which 588 genes were annotated. A total of 149,158 biallelic insertions/deletions, 6531 divergent alleles, and 14,707 multiallelic variations with precise breakpoints were discovered. The SV spectrum is characterized by an excess of derived insertions compared to deletions (94,422 vs. 33,571), suggesting recent active LINE expansions in sheep. Nearly half of the SVs display low to moderate linkage disequilibrium with surrounding single-nucleotide polymorphisms (SNPs) and most SVs cannot be tagged by SNP probes from the widely used ovine 50K SNP chip. We identified 865 population-stratified SVs including 122 SVs possibly derived in the domestication process among 690 individuals from sheep breeds worldwide. A novel 168-bp insertion in the 5′ untranslated region (5′ UTR) ofHOXB13is found at high frequency in long-tailed sheep. Further genome-wide association study and gene expression analyses suggest that this mutation is causative for the long-tail trait. In summary, we have developed a panel of high-quality de novo assemblies and present a catalog of structural variations in sheep. Our data capture abundant candidate functional variations that were previously unexplored and provide a fundamental resource for understanding trait biology in sheep.

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.

PSVIII-B-15 Potential Involvement of DNA Methylation of First Intron in Transcriptional Regulation During Bovine Subclinical Mastitis Caused by Staphylococcus Aureus

DNA methylation involvement in the regulation of mammary gland inflammatory defense against pathogens, including Staphylococcus aureus(SA) has been documented. However, the mechanisms are not clear. The possible influence of DNA methylation of first introns (DNA-MeFI) on transcriptional activities during subclinical mastitis was studied by profiling the methylome (whole genome methylation sequencing) and transcriptome (RNA-sequencing) of milk somatic cells from cows(n=15) with SA subclinical mastitis (SACs) and healthy cows (HCs)(n=13), followed by bioinformatics processing with standard tools. The DNA-MeFI was calculated, for each of 19411 genes found as expressed in the samples, by averaging the methylation levels of all CpG sites in the first intron. At genome-wide, the DNA-MeFIs was inversely correlated with gene expression (GE) (Pearson’s r=-0.173, p=1.60×10-122)(Figure 1-A). Similarly, the difference in the DNA-MeFI and GE between SACs and HCs was also inversely correlated (Pearson’s r=-0.149, p=4.25×10-97)(Figure 1-B). Application of machine learning by Gaussian Mixed model (GMM) using scikit-learn in python to the changes of DNA-MeFI and GE between SACs and HCs revealed 2866 outliers (GMM’s p&lt; 0.005) that showed significant changes in the DNA-MeFI and/or GE. 644 genes with &gt;10% difference in DNA-MeFI and |log2FC| &gt;1 in GE were selected and referred to as differentially methylated and expressed genes (DME-genes). 644 DME-genes were significantly enriched in 5 GO terms and 11 KEGG pathways related to diseases and immune functions, including Staphylococcus aureus infection and Natural killer cell mediated cytotoxicity (Table 1). Furthermore, 74.84% DME-genes showed inverse changes in DNA-MeFI and GE. Besides, the DNA-MeFI of 264 DME-genes were found to correlate significantly with their GE (|rho| &gt;0.3, FDR&lt;0.05), including 225 DME-genes(85.23%) with inverse correlations(rho&lt; -0.3, FDR&lt; 0.05). 39 DME-genes with positive correlations were significantly enriched in 21 GO terms mainly related to metabolic processes. In conclusion, the DNA-MeFI possibly participate in the regulation of gene expression during bovine subclinical mastitis caused by SA.

Mammary epithelial cell transcriptome reveals potential roles of lncRNAs in regulating milk synthesis pathways in Jersey and Kashmiri cattle

Background

Long noncoding RNAs (lncRNAs) are now proven as essential regulatory elements, playing diverse roles in many biological processes including mammary gland development. However, little is known about their roles in the bovine lactation process.

Results

To identify and characterize the roles of lncRNAs in bovine lactation, high throughput RNA sequencing data from Jersey (high milk yield producer), and Kashmiri cattle (low milk yield producer) were utilized. Transcriptome data from three Kashmiri and three Jersey cattle throughout their lactation stages were utilized for differential expression analysis. At each stage (early, mid and late) three samples were taken from each breed. A total of 45 differentially expressed lncRNAs were identified between the three stages of lactation. The differentially expressed lncRNAs were found co-expressed with genes involved in the milk synthesis processes such as GPAM, LPL, and ABCG2 indicating their potential regulatory effects on milk quality genes. KEGG pathways analysis of potential cis and trans target genes of differentially expressed lncRNAs indicated that 27 and 48 pathways were significantly enriched between the three stages of lactation in Kashmiri and Jersey respectively, including mTOR signaling, PI3K-Akt signaling, and RAP1 signaling pathways. These pathways are known to play key roles in lactation biology and mammary gland development.

Conclusions

Expression profiles of lncRNAs across different lactation stages in Jersey and Kashmiri cattle provide a valuable resource for the study of the regulatory mechanisms involved in the lactation process as well as facilitate understanding of the role of lncRNAs in bovine lactation biology.

Supplementary Information

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

Genome wide expression analysis of circular RNAs in mammary epithelial cells of cattle revealed difference in milk synthesis

Milk is an excellent source of nutrients for humans. Therefore, in order to enhance the quality and production of milk in cattle, it is interesting to examine the underlying mechanisms. A number of new investigations and research have found that, circRNA; a specific class of non-coding RNAs, is linked with the development of mammary gland and lactation. In the present study, genome wide identification and expression of the circRNAs in mammary epithelial cells of two distinct cattle breeds viz Jersey and Kashmiri at peak lactation was conducted. We reported 1554 and 1286 circRNA in Jersey and Kashmiri cattle, respectively, with 21 circRNAs being differentially expressed in the two breeds. The developmental genes of the established differentially expressed circRNAs were found to be largely enriched in antioxidant activity, progesterone, estradiol, lipid, growth hormone, and drug response. Certain pathways like MAPK, IP3K and immune response pathways were found significantly enriched in KEGG analysis. These results add to our understanding of the controlling mechanisms connected with the lactation process, as well as the function of circRNAs in bovine milk synthesis. Additionally, the comparative analysis of differentially expressed circRNAs showed significant conservation across different species.

Regionally Distinct Immune and Metabolic Transcriptional Responses in the Bovine Small Intestine and Draining Lymph Nodes During a Subclinical Mycobacterium avium subsp. paratuberculosis Infection

Mycobacterium avium subsp. paratuberculosis (MAP) is the causative infectious agent of Johne’s disease (JD), an incurable granulomatous enteritis affecting domestic livestock and other ruminants around the world. Chronic MAP infections usually begin in calves with MAP uptake by Peyer’s patches (PP) located in the jejunum (JE) and ileum (IL). Determining host responses at these intestinal sites can provide a more complete understanding of how MAP manipulates the local microenvironment to support its long-term survival. We selected naturally infected (MAPinf, n=4) and naive (MAPneg, n=3) cows and transcriptionally profiled the JE and IL regions of the small intestine and draining mesenteric lymph nodes (LN). Differentially expressed (DE) genes associated with MAP infection were identified in the IL (585), JE (218), jejunum lymph node (JELN) (205), and ileum lymph node (ILLN) (117). Three DE genes (CD14, LOC616364 and ENSBTAG00000027033) were common to all MAPinf versus MAPneg tissues. Functional enrichment analysis revealed immune/disease related biological processes gene ontology (GO) terms and pathways predominated in IL tissue, indicative of an activated immune response state. Enriched GO terms and pathways in JE revealed a distinct set of host responses from those detected in IL. Regional differences were also identified between the mesenteric LNs draining each intestinal site. More down-regulated genes (52%) and fewer immune/disease pathways (n=5) were found in the ILLN compared to a higher number of up-regulated DE genes (56%) and enriched immune/disease pathways (n=13) in the JELN. Immunohistochemical staining validated myeloid cell transcriptional changes with increased CD172-positive myeloid cells in IL and JE tissues and draining LNs of MAPinf versus MAPneg cows. Several genes, GO terms, and pathways related to metabolism were significantly DE in IL and JE, but to a lesser extent (comparatively fewer enriched metabolic GO terms and pathways) in JELN suggesting distinct regional metabolic changes in IL compared to JE and JELN in response to MAP infection. These unique tissue- and regional-specific differences provides novel insight into the dichotomy in host responses to MAP infection that occur throughout the small intestine and mesenteric LN of chronically MAP infected cows.

Whole Genome Methylation Analysis Reveals Role of DNA Methylation in Cow's Ileal and Ileal Lymph Node Responses to Mycobacterium avium subsp. paratuberculosis Infection

Johne's Disease (JD), caused by Mycobacterium avium subsp paratuberculosis (MAP), is an incurable disease of ruminants and other animal species and is characterized by an imbalance of gut immunity. The role of MAP infection on the epigenetic modeling of gut immunity during the progression of JD is still unknown. This study investigated the DNA methylation patterns in ileal (IL) and ileal lymph node (ILLN) tissues from cows diagnosed with persistent subclinical MAP infection over a one to 4 years period. DNA samples from IL and ILLN tissues from cows negative (MAPneg) (n = 3) or positive for MAP infection (MAPinf) (n = 4) were subjected to whole genome bisulfite sequencing. A total of 11,263 and 62,459 differentially methylated cytosines (DMCs), and 1259 and 8086 differentially methylated regions (DMRs) (FDR<0.1) were found between MAPinf and MAPneg IL and ILLN tissues, respectively. The DMRs were found on 394 genes (denoted DMR genes) in the IL and on 1305 genes in the ILLN. DMR genes with hypermethylated promoters/5'UTR [3 (IL) and 88 (ILLN)] or hypomethylated promoters/5'UTR [10 (IL) and 25 (ILLN)] and having multiple functions including response to stimulus/immune response (BLK, BTC, CCL21, AVPR1A, CHRNG, GABRA4, TDGF1), cellular processes (H2AC20, TEX101, GLA, NCKAP5L, RBM27, SLC18A1, H2AC20BARHL2, NLGN3, SUV39H1, GABRA4, PPA1, UBE2D2) and metabolic processes (GSTO2, H2AC20, SUV39H1, PPA1, UBE2D2) are potential DNA methylation candidate genes of MAP infection. The ILLN DMR genes were enriched for more biological process (BP) gene ontology (GO) terms (n = 374), most of which were related to cellular processes (27.6%), biological regulation (16.6%), metabolic processes (15.4%) and response to stimulus/immune response (8.2%) compared to 75 BP GO terms (related to cellular processes, metabolic processes and transport, and system development) enriched for IL DMR genes. ILLN DMR genes were enriched for more pathways (n = 47) including 13 disease pathways compared with 36 enriched pathways, including 7 disease/immune pathways for IL DMR genes. In conclusion, the results show tissue specific responses to MAP infection with more epigenetic changes (DMCs and DMRs) in the ILLN than in the IL tissue, suggesting that the ILLN and immune processes were more responsive to regulation by methylation of DNA relative to IL tissue. Our data is the first to demonstrate a potential role for DNA methylation in the pathogenesis of MAP infection in dairy cattle.

Whole Genome DNA Methylation Variations in Mammary Gland Tissues from Holstein Cattle Producing Milk with Various Fat and Protein Contents

Milk fat and protein contents are among key elements of milk quality, and they are attracting more attention in response to consumers' demand for high-quality dairy products. To investigate the potential regulatory roles of DNA methylation underlying milk component yield, whole genome bisulfite sequencing was employed to profile the global DNA methylation patterns of mammary gland tissues from 17 Canada Holstein cows with various milk fat and protein contents. A total of 706, 2420 and 1645 differentially methylated CpG sites (DMCs) were found between high vs. low milk fat (HMF vs. LMF), high vs. low milk protein (HMP vs. LMP), and high vs. low milk fat and protein (HMFP vs. LMFP) groups, respectively (q value < 0.1). Twenty-seven, 56 and 67 genes harboring DMCs in gene regions (denoted DMC genes) were identified for HMF vs. LMF, HMP vs. LMP and HMFP vs. LMFP, respectively. DMC genes from HMP vs. LMP and HMFP vs. LMFP comparisons were significantly overrepresented in GO terms related to aerobic electron transport chain and/or mitochondrial ATP (adenosine triphosphate) synthesis coupled electron transport. A total of 83 (HMF vs. LMF), 708 (HMP vs. LMP) and 408 (HMFP vs. LMFP) DMCs were co-located with 87, 147 and 158 quantitative trait loci (QTL) for milk component and yield traits, respectively. In conclusion, the identified methylation changes are potentially involved in the regulation of milk fat and protein yields, as well as the variation in reported co-located QTLs.

156 Genome-wide DNA Methylation Profile of Mammary Gland Tissues from Holstein Cattle Producing Various Milk Protein Yields

Milk protein content is attracting considerable attention due to increasing demand for high-quality products. Epigenetic marks including DNA methylation have demonstrated considerable potential to contribute to the achievement of increased gain in bovine milk traits. To investigate potential roles of DNA methylation in regulating milk protein productivity, whole genome bisulfite sequencing was used to profile the global DNA methylation of mammary gland (MG) tissues from Canadian Holstein cows producing milk with high protein (HMP) (&gt;4.0%, n = 5) or low protein (LMP) (&lt; 3.0%, n = 7) contents. Bioinformatics processing using standard tools and construction of the methylation landscape of MG was based on methylation sites with ≥10× coverage. A total of 2420 differentially methylated CpG (DMC) sites, including 1842 with higher and 578 with lower methylation levels (FDR&lt; 0.1) were identified when HMP group was compared with LMP group. DMCs were most abundant on chromosomes 2 (n=347) and 27 (N=535), respectively. The most DMCs were annotated to mRNAs, followed by lncRNAs, pseudogenes and rRNAs. Eight genes harbored five or more DMCs each, including LOC112443250 (5S ribosomal RNA) and LOC112443149 (liprin-alpha-1-like, LIP1) (harbored 20 hypermethylated and 1 hypomethylated sites, each), followed by miR-2287-1 (12 hypomethylated sites), etc. LIP1, a liprin family member interacts with LAR family members to play important roles in MG development. MiR-2287-1 and other family members are frequently identified in MG tissues and milk and are associated with many livestock diseases. Interestingly, among the DMCs that overlap with QTLs related to milk production were clusters on chromosomes 27 (N=534) and 14 (N=91) that associated with milk protein yield. In conclusion, DNA methylation showed potential involvement in the genetic variation underlying milk protein content and MG biological processes. More analysis with a higher sample size is however necessary to adequately quantify the involvement of DNA methylation in the regulation of milk protein content and mammary gland development.

72 A One-health Perspective on S. aureus Induced Mastitis

Bovine mastitis, commonly caused by bacteria, is one of the most devastating diseases for the dairy industry worldwide, with high economic, social, and public health impacts. Among the pathogenic bacteria, Staphylococcus aureus is one of the leading pathogens in most countries. Success rate of antibiotic treatment is low such that effective control of S. aureus induced mastitis is currently only possible through prevention of new infections and culling of infected animals. The infection is usually subclinical, causing elevated somatic cell counts but no detectable changes in milk or the udder. S. aureus persists in mammary glands, teat canals, and teat lesions of infected cows and is transmitted to other cows during milking. The direct evidence of cross-species transmission of S. aureus is still scarce. Nevertheless, accumulating indirect evidence supports the zoonotic potential for S. aureus from bovine mastitis to humans. In addition to direct contact, aerosol exposures may also be an important mechanism for the transmission. While animal-adapted livestock-associated methicillin-resistant S. aureus has been known for many years, most S. aureus isolated from mastitic milk in Canada are methicillin-sensitive. S. aureus induced mastitis is a common reason for therapeutic and/or prophylactic use of antibiotics on dairy farms. Occurrence of multidrug resistant, especially methicillin-resistant, S. aureus has been a major public health concern. In this talk, we will address the pathogenesis and strain spectrum of S. aureus induced bovine mastitis, review existing evidence for inter-species transmission, and discuss possible transmission of methicillin-resistance determinants. This talk aims to emphasize the need for the one-health approach for prevention and treatment of S. aureus induced mastitis, in order to provide safe and nutritious milk and milk products to consumers.

Impacts of Epigenetic Processes on the Health and Productivity of Livestock

The dynamic changes in the epigenome resulting from the intricate interactions of genetic and environmental factors play crucial roles in individual growth and development. Numerous studies in plants, rodents, and humans have provided evidence of the regulatory roles of epigenetic processes in health and disease. There is increasing pressure to increase livestock production in light of increasing food needs of an expanding human population and environment challenges, but there is limited related epigenetic data on livestock to complement genomic information and support advances in improvement breeding and health management. This review examines the recent discoveries on epigenetic processes due to DNA methylation, histone modification, and chromatin remodeling and their impacts on health and production traits in farm animals, including bovine, swine, sheep, goat, and poultry species. Most of the reports focused on epigenome profiling at the genome-wide or specific genic regions in response to developmental processes, environmental stressors, nutrition, and disease pathogens. The bulk of available data mainly characterized the epigenetic markers in tissues/organs or in relation to traits and detection of epigenetic regulatory mechanisms underlying livestock phenotype diversity. However, available data is inadequate to support gainful exploitation of epigenetic processes for improved animal health and productivity management. Increased research effort, which is vital to elucidate how epigenetic mechanisms affect the health and productivity of livestock, is currently limited due to several factors including lack of adequate analytical tools. In this review, we (1) summarize available evidence of the impacts of epigenetic processes on livestock production and health traits, (2) discuss the application of epigenetics data in livestock production, and (3) present gaps in livestock epigenetics research. Knowledge of the epigenetic factors influencing livestock health and productivity is vital for the management and improvement of livestock productivity.

Comparative accuracies of genetic values predicted for economically important milk traits, genome-wide association, and linkage disequilibrium patterns of Canadian Holstein cows

Genomic selection methodologies and genome-wide association studies use powerful statistical procedures that correlate large amounts of high-density SNP genotypes and phenotypic data. Actual 305-d milk (MY), fat (FY), and protein (PY) yield data on 695 cows and 76,355 genotyping-by-sequencing-generated SNP marker genotypes from Canadian Holstein dairy cows were used to characterize linkage disequilibrium (LD) structure of Canadian Holstein cows. Also, the comparison of pedigree-based BLUP, genomic BLUP (GBLUP), and Bayesian (BayesB) statistical methods in the genomic selection methodologies and the comparison of Bayesian ridge regression and BayesB statistical methods in the genome-wide association studies were carried out for MY, FY, and PY. Results from LD analysis revealed that as marker distance decreases, LD increases through chromosomes. However, unexpected high peaks in LD were observed between marker pairs with larger marker distances on all chromosomes. The GBLUP and BayesB models resulted in similar heritability estimates through 10-fold cross-validation for MY and PY; however, the GBLUP model resulted in higher heritability estimates than BayesB model for FY. The predictive ability of GBLUP model was significantly lower than that of BayesB for MY, FY, and PY. Association analyses indicated that 28 high-effect markers and markers on Bos taurus autosome 14 located within 6 genes (DOP1B, TONSL, CPSF1, ADCK5, PARP10, and GRINA) associated significantly with FY.

Cholesterol deficiency haplotype frequency and its impact on milk production and milk cholesterol content in Canadian Holstein cows

An autosomal recessive defect caused by a loss of function mutation in the bovine apolipoprotein B (APOB) gene causes cholesterol deficiency with strong effects on calf survival and development. This study examined the distribution of cholesterol deficiency haplotype (CDH) in the APOB gene in Canadian Holstein cows and its impacts on milk production and milk cholesterol content. The CDH in 917 cows from 29 herds was determined via polymerase chain reaction (PCR). A longer insertion of approximately 7 kb potentially characterizing the mutation was amplified. The average carrier frequency of herds in Quebec based on PCR results was 14.6%. Cholesterol deficiency haplotype status had non-significant (P > 0.05) effects on milk cholesterol content, other milk components, and milk yield.

PSVIII-15 Genome wide DNA methylation analysis reveals role of DNA methylation in cow’s ileal response to Mycobacterium avium subsp. paratuberculosis

Several investigations on disease progression of Mycobacterium avium subsp. paratuberculosis (MAP) infection in dairy cows have revealed molecular mechanisms including genes and pathways implicated in MAP pathogenesis. Epigenetic processes including DNA methylation are known to regulate the expression of genes and many biological processes. However, limited studies have examined the role of DNA methylation in the pathogenesis of Johne’s disease (JD). This study examined the impact of subclinical MAP infection on DNA methylation profile in the ileum of cows, the site of initial interaction between MAP and host. DNA from ileum tissues from five cows positive for MAP (direct fecal qPCR and blood ELISA test; MAP+/+) and 5 cows negative for MAP (MAP-/-) were subjected to whole genome bisulfite sequencing and bioinformatics analysis. 2000 differentially methylated cytosines (DMCs; FDR&lt; 0.05) and 205 differentially methylated regions (DMRs; P &lt; 0.01) were detected. Majority of DMCs and DMRs are located in intergenic regions (87.2% and 57.1%) followed by intronic regions (12.8% and 30.7%) of genes, respectively. Some DMCs are located on 250 genes including genes that were previously identified to be associated with JD (e.g. IL-12RB2, CD38). Interestingly, CD38, known to play roles in the effective containment of mycobacteria within granulomata in cows and genetic polymorphisms in IL-12RB2 are associated with JD and human Crohn’s disease. Also, several genes of the solute carrier family including SLC13A3, SLC15A1, SLC17A7, SLC25A21, SLC25A38 and SLC9A9 harbored DMCs. Some members of this gene family participate in pathogen clearance and have associations with JD. Our data suggest that DNA methylation changes may have regulatory roles in host (ileum) response to MAP infection. Our data therefore suggest that DNA methylation changes contribute to the regulation of host response to MAP pathogenesis and may be one of the mechanisms that MAP uses to subvert host immune responses for its survival.

29 Conditional GWAS using sequence-based genotypes for susceptibility to Mycobacterium avium subsp paratuberculosis infection in Canadian Holstein

Paratuberculosis is an incurable disease caused by Mycobacterium avium spp. paratuberculosis (MAP). Variability in the definition of a true infection could explain the lack of consensus in published results on host genetic susceptibility to MAP. The objective of this study was to identify quantitative trait loci (QTL) for susceptibility to paratuberculosis in Canadian Holstein with an accurately defined status. Cows with known parentage were recruited from 23 Canadian dairy herds and enrolled in a longitudinal study. Blood and feces were collected semi-annually over a 3-5-yr period. Serum samples were analyzed by ELISA (IDEXX) and feces by qPCR. Cows were labelled based on positive ELISA (&gt; 45%) and on at least two consecutive fecal excretion of MAP (n = 382). Control cows (n = 669) did not excrete MAP nor have MAP antibodies (&lt; 10%) over the 3–5 year period. Their genotypes (GGP Bovine 50K chip) were imputed to HD-SNP (within-breed), then to WGS using FImpute software and Run6 of 1,000 Bulls’ genome project as a multi-breed reference. Using GCTA software, heritability (h2) was estimated and conditional GWAS was performed to define multiple true causal mutations in QTL regions from those caused by linkage disequilibrium. The h2(se) was 0.52 (0.01) and after FDR correction, 1,357 variants presented significant effects (-log10 P &gt; 6), grouped within 7 QTL located on chromosome BTA7, BTA16, BTA20, and BTA23. The most significant variant was on BTA20 (-log10 P =12.39). Functional annotation showed that 76% of all the variants of the QTL were intergenic, with the rest being mostly intronic. Of interest was variant located on BTA7 with effect size (se) of -0.77 (0.02). In conclusion, using accurate phenotypes, chromosomal regions and candidate genes that could be functionally related to MAP infection were uncovered. The significant variants found in this study could contribute to the enhancement of genomic evaluation for MAP resistance.

63 DNA methylome wide profile associates differentially methylated loci and regions with cow’s ileal lymph node response to Mycobacterium avium subsp. paratuberculosis

Johne’s Disease (JD), caused by Mycobacterium avim spp paratuberculosis (MAP), is a chronic and incurable disease of ruminants with devastating consequences to the dairy industry. MAP can alter the expression of genes and biological processes during the progression of JD. While some studies have examined the role of gene expression regulators like microRNA in the pathogenesis of JD, no study has explored the role of DNA methylation. This study therefore examined the effect of MAP on DNA methylation profile in the ileum lymph node (ILLN) of cows with subclinical MAP infection. DNA from ILLN tissues from five cows positive for MAP (MAP positive) and five negative cows (MAP negative) were extracted and subjected to whole genome bisulfite sequencing and bioinformatics analysis. A total of 6,394 differentially methylated cytosines (DMCs) and 3,946 differentially methylated regions (DMRs) (FDR &lt; 0.05) were identified between MAP positive and negative cows. DMRs were annotated to 2,488 genes, including the promoters of 238 genes. Some genes with hypermethylated promoters like GRB10, EIF4E, SLC30A3, SOX30 or hypomethylated promoters like SLC11A1, HOXA4, SLC18A1 have been associated with JD or mycobacterial infections in cattle and/or humans. Functional annotation of DMR genes indicated enrichment in pathways previously associated with JD or human diseases with similar pathological conditions as JD, such as T/B cell receptor signaling pathway, Th17 cell differentiation, Cell adhesion molecules, Leukocyte transendothelial migration, HIF-1 signaling pathway and Chagas disease. Furthermore, enriched gene ontology terms like negative regulation of immune system process, negative regulation of cytokine secretion/production and negative regulation of inflammatory response suggest that MAP prevented or reduced the host immune response. Our data demonstrate that DNA methylation changes contribute to regulation of host immune responses to MAP infection and may be one of the mechanisms that MAP uses to subvert host immune responses for its long-term survival.

Genome-Wide DNA Methylation Analysis of Mammary Gland Tissues From Chinese Holstein Cows With Staphylococcus aureus Induced Mastitis

Staphylococcus aureus intramammary infection is one of the most common causes of chronic mastitis in dairy cows, whose development may be associated with epigenetic changes in the expression of important host defense genes. This study aimed to construct a genome-wide DNA methylation profile of the mammary gland of Chinese Holstein cows (n = 3) following experimentally induced S. aureus mastitis, and to explore the potential gene regulatory mechanisms affected by DNA methylation during S. aureus mastitis. DNA was extracted from S. aureus-positive (n = 3) and S. aureus-negative (n = 3) mammary gland quarters and subjected to methylation-dependent restriction-site associated DNA sequencing (Methyl-RAD Seq). Results showed that CmCGG/CmCWGG DNA methylation sites were unevenly distributed and concentrated on chromosomes 5, 11, and 19, and within intergenic regions and intron regions of genes. Compared with healthy control quarters, 9,181 significantly differentially methylated (DM) CmCGG sites and 1,790 DM CmCWGG sites were found in the S. aureus-positive quarters (P < 0.05, |log2FC| > 1). Furthermore, 363 CmCGG differently methylated genes (DMGs) and 301 CmCWGG DMGs (adjusted P < 0.05, |log2FC| > 1) were identified. Gene ontology and KEGG enrichment analysis indicated that CmCGG DMGs are involved in immune response pathways, while the CmCWGG DMGs were mainly enriched in gene ontology terms related to metabolism. The mRNAs of 526 differentially methylated CmCGG genes and 124 differentially methylated CmCWGG genes were also significantly differentially expressed (RNA-Seq data) in the same samples, herein denoted differentially methylated and expressed genes (DMEGs) (P < 0.05). Functional enrichment analysis of DMEGs revealed roles related to biological processes, especially the regulation of immune response to diseases. CmCGG DMEGs like IL6R, TNF, BTK, IL1R2, and TNFSF8 enriched in several immune-related GO terms and pathways indicated their important roles in host immune response and their potential as candidate genes for S. aureus mastitis. These results suggest potential regulatory roles for DNA methylation in bovine mammary gland processes during S. aureus mastitis and serves as a reference for future epigenetic regulation and mechanistic studies.

Distinct miRNA Profile of Cellular and Extracellular Vesicles Released from Chicken Tracheal Cells Following Avian Influenza Virus Infection

Innate responses provide the first line of defense against viral infections, including the influenza virus at mucosal surfaces. Communication and interaction between different host cells at the early stage of viral infections determine the quality and magnitude of immune responses against the invading virus. The release of membrane-encapsulated extracellular vesicles (EVs), from host cells, is defined as a refined system of cell-to-cell communication. EVs contain a diverse array of biomolecules, including microRNAs (miRNAs). We hypothesized that the activation of the tracheal cells with different stimuli impacts the cellular and EV miRNA profiles. Chicken tracheal rings were stimulated with polyI:C and LPS from Escherichia coli 026:B6 or infected with low pathogenic avian influenza virus H4N6. Subsequently, miRNAs were isolated from chicken tracheal cells or from EVs released from chicken tracheal cells. Differentially expressed (DE) miRNAs were identified in treated groups when compared to the control group. Our results demonstrated that there were 67 up-regulated miRNAs, 157 down-regulated miRNAs across all cellular and EV samples. In the next step, several genes or pathways targeted by DE miRNAs were predicted. Overall, this study presented a global miRNA expression profile in chicken tracheas in response to avian influenza viruses (AIV) and toll-like receptor (TLR) ligands. The results presented predicted the possible roles of some DE miRNAs in the induction of antiviral responses. The DE candidate miRNAs, including miR-146a, miR-146b, miR-205a, miR-205b and miR-449, can be investigated further for functional validation studies and to be used as novel prophylactic and therapeutic targets in tailoring or enhancing antiviral responses against AIV.

Comparative milk proteome analysis of Kashmiri and Jersey cattle identifies differential expression of key proteins involved in immune system regulation and milk quality

Background

Exploration of the bioactive components of bovine milk has gained global interest due to their potential applications in human nutrition and health promotion. Despite advances in proteomics profiling, limited studies have been carried out to fully characterize the bovine milk proteome. This study explored the milk proteome of Jersey and Kashmiri cattle at day 90 of lactation using high-resolution mass spectrometry based quantitative proteomics nano-scale LC-MS/Q-TOF technique. Data are available via ProteomeXchange with identifier PXD017412.

Results

Proteins from whey were fractionated by precipitation into high and low abundant proteins. A total of 81 high-abundant and 99 low-abundant proteins were significantly differentially expressed between Kashmiri and Jersey cattle, clearly differentiating the two breeds at the proteome level. Among the top differentiating proteins, the Kashmiri cattle milk proteome was characterised by increased concentrations of immune-related proteins (apelin, acid glycoprotein, CD14 antigen), neonatal developmental protein (probetacellulin), xenobiotic metabolising enzyme (flavin monooxygenase 3 (FMO3), GLYCAM1 and HSP90AA1 (chaperone) while the Jersey milk proteome presented higher concentrations of enzyme modulators (SERPINA1, RAC1, serine peptidase inhibitor) and hydrolases (LTF, LPL, CYM, PNLIPRP2). Pathway analysis in Kashmiri cattle revealed enrichment of key pathways involved in the regulation of mammary gland development like Wnt signalling pathway, EGF receptor signalling pathway and FGF signalling pathway while a pathway (T-cell activation pathway) associated with immune system regulation was significantly enriched in Jersey cattle. Most importantly, the high-abundant FMO3 enzyme with an observed 17-fold higher expression in Kashmiri cattle milk seems to be a characteristic feature of the breed. The presence of this (FMO3) bioactive peptide/enzyme in Kashmiri cattle could be economically advantageous for milk products from Kashmiri cattle.

Conclusion

In conclusion, this is the first study to provide insights not only into the milk proteome differences between Kashmiri and Jersey cattle but also provides potential directions for application of specific milk proteins from Kashmiri cattle in special milk preparations like infant formula.

Transcriptome Profiling of Bovine Macrophages Infected by Mycobacterium avium spp. paratuberculosis Depicts Foam Cell and Innate Immune Tolerance Phenotypes

Mycobacterium avium spp. paratuberculosis (MAP) is the causative agent of Johne's disease (JD), also known as paratuberculosis, in ruminants. The mechanisms of JD pathogenesis are not fully understood, but it is known that MAP subverts the host immune system by using macrophages as its primary reservoir. MAP infection in macrophages is often studied in healthy cows or experimentally infected calves, but reports on macrophages from naturally infected cows are lacking. In our study, primary monocyte-derived macrophages (MDMs) from cows diagnosed as positive (+) or negative (–) for JD were challenged in vitro with live MAP. Analysis using next-generation RNA sequencing revealed that macrophages from JD(+) cows did not present a definite pattern of response to MAP infection. Interestingly, a considerable number of genes, up to 1436, were differentially expressed in JD(–) macrophages. The signatures of the infection time course of 1, 4, 8, and 24 h revealed differential expression of ARG2, COL1A1, CCL2, CSF3, IL1A, IL6, IL10, PTGS2, PTX3, SOCS3, TNF, and TNFAIP6 among other genes, with major effects on host signaling pathways. While several immune pathways were affected by MAP, other pathways related to hepatic fibrosis/hepatic stellate cell activation, lipid homeostasis, such as LXR/RXR (liver X receptor/retinoid X receptor) activation pathways, and autoimmune diseases (rheumatoid arthritis or atherosclerosis) also responded to the presence of live MAP. Comparison of the profiles of the unchallenged MDMs from JD(+) vs. JD(–) cows showed that 868 genes were differentially expressed, suggesting that these genes were already affected before monocytes differentiated into macrophages. The downregulated genes predominantly modified the general cell metabolism by downregulating amino acid synthesis and affecting cholesterol biosynthesis and other energy production pathways while introducing a pro-fibrotic pattern associated with foam cells. The upregulated genes indicated that lipid homeostasis was already supporting fat storage in uninfected JD(+) MDMs. For JD(+) MDMs, differential gene expression expounds long-term mechanisms established during disease progression of paratuberculosis. Therefore, MAP could further promote disease persistence by influencing long-term macrophage behavior by using both tolerance and fat-storage states. This report contributes to a better understanding of MAP's controls over the immune cell response and mechanisms of MAP survival.

Targeted genotyping to identify potential functional variants associated with cholesterol content in bovine milk

High blood cholesterol concentration, mainly caused by high dietary cholesterol, is a potential risk factor for human health. Dairy products are important sources of human dietary cholesterol intake. Therefore, monitoring bovine milk cholesterol concentration is important for human health benefit. Genetic selection for improvement of cow milk cholesterol content requires understanding of the genetics of milk cholesterol. For this purpose, we performed analyses of additive and dominance effects of 126 potentially functional SNPs within 43 candidate genes with milk cholesterol content [expressed as mg of cholesterol in 100 g of fat (CHL_fat) or in 100 mg of milk (CHL_milk)]. The additive and dominance effects of SNPs rs380643365 in AGPAT1 (P = 0.04) and rs134357240 in SOAT1 (P = 0.035) genes associated significantly with CHL_fat. Moreover, five (rs109326954 and rs523413537 in DGAT1, rs109376747 in LDLR, rs42781651 in FAM198B and rs109967779 in ACAT2) and four (rs137347384 in RBM19, rs109376747 in LDLR, rs42016945 in PPARG and rs110862179 in SCAP) SNPs were significantly associated with CHL_milk (P < 0.05) based on additive and dominance effect analyses respectively. Rs109326954 and rs523413537 in DGAT1 explained a considerable portion of the phenotypic variance of CHL_milk (7.54 and 6.84% respectively), and might be useful in selection programs for reduced milk cholesterol content. Several significantly associated SNPs were in genes (such as ACAT2 and LDLR) involved in cholesterol metabolism in the liver or cholesterol transport, suggesting multiple mechanisms regulating milk cholesterol content. Nine and seven SNPs identified by additive or dominance effect analyses associated significantly with milk yield and fat yield respectively. Further analyses are required to better understand the consequences of these variants and their potential use in genomic selection of the studied traits.

A targeted genotyping approach to enhance the identification of variants for lactation persistency in dairy cows

Lactation persistency (LP), defined as the ability of a cow to maintain milk production at a high level after milk peak, is an important phenotype for the dairy industry. In this study, we used a targeted genotyping approach to scan for potentially functional single nucleotide polymorphisms (SNPs) within 57 potential candidate genes derived from our previous genome wide association study on LP and from the literature. A total of 175,490 SNPs were annotated within 10-kb flanking regions of the selected candidate genes. After applying several filtering steps, a total of 105 SNPs were retained for genotyping using target genotyping arrays. SNP association analyses were performed in 1,231 Holstein cows with 69 polymorphic SNPs using the univariate liner mixed model with polygenic effects using DMU package. Six SNPs including rs43770847, rs208794152, and rs208332214 in ADRM1; rs209443540 in C5orf34; rs378943586 in DDX11; and rs385640152 in GHR were suggestively significantly associated with LP based on additive effects and associations with 4 of them (rs43770847, rs208794152, rs208332214, and rs209443540) were based on dominance effects at P < 0.05. However, none of the associations remained significant at false discovery rate adjusted P (FDR) < 0.05. The additive variances explained by each suggestively significantly associated SNP ranged from 0.15% (rs43770847 in ADRM1) to 5.69% (rs209443540 in C5orf34), suggesting that these SNPs might be used in genetic selection for enhanced LP. The percentage of phenotypic variance explained by dominance effect ranged from 0.24% to 1.35% which suggests that genetic selection for enhanced LP might be more efficient by inclusion of dominance effects. Overall, this study identified several potentially functional variants that might be useful for selection programs for higher LP. Finally, a combination of identification of potentially functional variants followed by targeted genotyping and association analysis is a cost-effective approach for increasing the power of genetic association studies.

PSVI-14 Differentially expressed microRNAs with potential regulatory roles in ileum of Holstein cows with subclinical Johne’s disease

Paratuberculosis or Johne’s disease (JD), caused by Mycobacterium avium subsp. paratuberculosis (MAP), is a chronic inflammatory disorder of the intestinal tract of ruminants. The mechanisms of JD pathogenesis remain unclear but MAP can cause chronic inflammatory lesions in the terminal small intestine of cattle. Using microRNA (miRNA) transcriptome analysis, this study investigated the impact of MAP infection on miRNA expression and their potential regulatory roles in JD progression. Ileal tissues were collected from MAP infected Holstein cows (n = 5) and healthy cows (n = 5). Purified total RNA were used to construct libraries followed by next generation sequencing, bioinformatic processing using a standard pipeline, differential miRNA expression analysis with DESeq2 software, and functional enrichment with ClueGo app in Cytoscape. A total of 269 miRNAs were expressed in ileal tissue including three abundantly expressed miRNAs (Bta-miR-143; 54.05% of total read counts), bta-miR-145 (11.54%) and bta-miR-192 (4.31%)). Twenty miRNAs were significantly (P &lt; 0.05) differently expressed (DE) when comparing ileal tissues from MAP infected and uninfected cows. Bta-miR-383 (fold change = -12.93) and bta-miR-1842 (fold-change = 8.17) were the most down- and up-regulated miRNAs, respectively. Amongst the DE miRNAs, bta-miR-146a, bta-miR-146b and bta-miR-393 are inflammation-related miRNAs and considered potential biomarkers of bovine mastitis and other diseases, while bta-miR-370 and bta-miR-383 are associated with lipid metabolism. The predicted target genes of DE miRNAs were enriched (P &lt; 0.05) for gene ontology terms related to regulation of metabolic, cell and developmental processes, and immune system development. KEGG pathway analysis revealed potential involvement in immune regulation (e.g. T-cell receptor signaling and C-type lectin receptor signaling pathways) and pathways in cancer were also enriched (P &lt; 0.05). Our results indicate that some DE miRNAs play important roles by regulating metabolism and immune processes during JD and provide DE miRNAs that might form targets for future development of intervention measures.

23 Differential microRNA expression in jejunal tissue and jejunal lymph nodes following naturally occurring Mycobacterium avium subspecies paratuberculosis infection in Holstein cows

Mycobacterium avium subsp. Paratuberculosis (MAP) is the causal agent of Johne’s disease (JD), a chronic intestinal disease affecting ruminants worldwide. This study investigated miRNA expression in jejunal intestine (JE) and jejunal lymph nodes (JELN), and the potential regulatory roles of miRNAs during JD progression. JE and JELN tissues were collected from 5 MAP positive (JD subclinical stage) Holstein cows and 5 MAP negative cows. Following miRNA sequencing, bioinformatic processing with a standard pipeline and functional analysis with ClueGo, 272 and 333 miRNAs were identified in JE and JELN, respectively. Compared with MAP negative cows, 13 and 71 miRNAs were differently expressed (DE) (P &lt; 0.05) in MAP infected JE and JELN, respectively. The most up-regulated and down-regulated miRNAs were bta-miR-485 (fold change = 6.18) and bta-miR-451 (fold change = -6.81), and bta-miR-331-5p (fold change = 35.56) and bta-miR-2285bk (fold change = -61.25) in JE and JELN, respectively. In JE, some DE miRNAs (miR-147 and miR-199a-5p) have been associated with glucose metabolism while several DE miRNAs in JELN have associations with immune response to disease (e.g. miR-146a and miR-146b have been associated with bovine mastitis and miR331-5p and miR-184 are important for human cancers). Target genes of JELN DE miRNAs were enriched (P &lt; 0.05) for more gene ontology terms (n = 180) and KEGG pathways (n = 123) as compared with JE (n = 90 and 25, respectively). Furthermore, JELN DE miRNAs were uniquely enriched in several immune related pathways including immune response regulating/activating signal transduction, T-cell/B-cell receptor signaling, Toll-like receptor signaling, TNF signaling pathways, etc., suggesting potential regulatory roles for JELN DE miRNAs in the host immune response to JD. In addition, interactions between DE miRNAs and DE mRNAs (expressed in the same samples) demonstrated important functions for miR-2284/2285 family members in JELN response to MAP infection. In conclusion, these results suggest site specific regulatory roles for miRNAs during MAP infection.

PSVI-15 Transcriptome analysis of ileal lymph nodes identifies key microRNAs affecting disease progression in Holstein cows with subclinical Johne’s disease

Johne’s disease (JD) is a chronic enteric disease of ruminants caused by Mycobacterium avium subsp. Paratuberculosis (MAP) but immune and metabolic mechanisms affecting JD progression have not been clearly elucidated. MicroRNAs (miRNA) have been reported to play important roles in numerous biological processes by regulating gene expression. In this study, miRNA expression in ileal lymph nodes (ILLN) from MAP positive (JD subclinical stage) Holstein cows (n = 5) and MAP negative cows (n = 5) were characterized by high throughput sequencing. MiRNA-sequence data was processed using a standard pipeline followed by functional enrichment with ClueGo app in Cytoscape. Among 324 expressed miRNAs, 37 were differently expressed (DE) (18 up-regulated and 19 down-regulated) in MAP infected cows relative to MAP negative cows. Bta-miR-100 (fold change = -2.29), a down-regulated miRNA, has previously been associated with Streptococcus uberis and MAP infections in cows, while bta-miR-330 (fold change = 63.42), the most up-regulated miRNA, was previously reported to be present in the sera of MAP infected calves. Predicted target genes of 37 DE miRNAs were enriched (P &lt; 0.05) in gene ontology terms related to metabolic, cellular and developmental processes, and immune system development, as well as pathways related to cancer and Th17 cell differentiation, and lysosome compartmentalization. This suggests that miRNAs DE following MAP infection may be involved in regulating a wide variety of processes. In addition, 5 DE miRNAs were negatively correlated (P &lt; 0.05) with 39 DE mRNAs (mRNAs expressed in the same tissue), comprising 42 negatively correlated DE miRNA-mRNA pairs. Bta-miR-2447 was negatively correlated with 29 DE mRNAs, suggesting it may be an important hub miRNA in the regulation of MAP infection in ILLN tissue. Our data thus suggest that bta-miR-100, bta-miR-330 and bta-miR-2447, which are implicated in the Th17 cell differentiation pathway, may be key regulators of the host response during MAP infection in ILLN.

Identification and characterization of differentially expressed exosomal microRNAs in bovine milk infected with Staphylococcus aureus

Background

MicroRNAs (miRNAs) in milk-derived exosomes may reflect pathophysiological changes caused by mastitis. This study profiled miRNAs in exosomes from both normal milk and mastitic milk infected by Staphylococcus aureus (S. aureus). The potential targets for differentially expressed (DE) miRNAs were predicted and the target genes for bta-miR-378 and bta-miR-185 were also validated.

Results

Total RNA from milk exosomes was collected from healthy cows (n = 3, the control group) and S. aureus infected cows (n = 6, the SA group). Two hundred ninety miRNAs (221 known and 69 novel ones) were identified. Among them, 22 known and 15 novel miRNAs were differentially expressed. Target genes of DE miRNAs were significantly enriched in intracellular protein transport, endoplasmic reticulum and identical protein binding. The expression of two miRNAs (bta-miR-378 and bta-miR-185) with high read counts and log₂ fold changes (> 3.5) was significantly higher in mastitic milk infected with S. aureus. One target gene (VAT1L) of bta-miR-378 and five target genes (DYRK1B, MLLT3, HP1BP3, NPR2 and PGM1) of bta-miR-185 were validated.

Conclusion

DE miRNAs in exosomes from normal and S. aureus infected milk were identified. The predicted targets for two DE miRNAs (bta-miR-378 and bta-miR-185) were further validated. The linkage between the validated target genes and diseases suggested that we should pay particular attention to exosome miRNAs from mastitic milk in terms of milk safety.

Characterization of gonadotropin releasing hormone receptor gene in Sokoto and Kalahari Red goats

The objective of this study was to characterize GnRHR gene in Sokoto (n = 70) and Kalahari Red (n = 70) goats. Three SNPs, (g.-29T > G, g.48 G > A and g.209 T > G), were detected in Sokoto Red (SR) and one (g.48 G > A) in Kalahari Red (KR) goats. All the mutations occurred within the 5'UTR and Exon one of the gene and the g.209 T > G was non-synonymous and, therefore, resulted in an amino acid change from methionine to arginine at Position 70 of the GnRHR polypeptide. The homozygous mutant genotypes at the three SNP loci were not detected in both breeds but minor allele frequencies were ≥ 0.1 for the three SNP loci in SR goats. Frequency of the T allele, however, was 0.93 at the only SNP locus detected in KR goats. There was a strong linkage disequilibrium (LD; r²>0.98) among the detected mutations in SR goats resulting in two haplotypes (T-G-T and G-A-G) with a frequency of 86% and 13%, respectively. There was no significant association between genotypes at the polymorphic loci and litter size (P > 0.05) in the two breeds. The non-synonymous mutation (g.209T>G) appears to have changed the nucleotide binding region and area spanning exposed/buried regions on the predicted secondary structure of the two variants of the receptor. This change led to variation in the tertiary structure between the two variants of the receptor and can influence the function of the transmembrane receptor. Comparison of the GnRH receptor domains for goats, sheep, cattle and swine confirmed that the seven transmembrane domains of the receptor are conserved in all the farm animals considered.

Leveraging Available Resources and Stakeholder Involvement for Improved Productivity of African Livestock in the Era of Genomic Breeding

The African continent is home to diverse populations of livestock breeds adapted to harsh environmental conditions with more than 70% under traditional systems of management. Animal productivity is less than optimal in most cases and is faced with numerous challenges including limited access to adequate nutrition and disease management, poor institutional capacities and lack of adequate government policies and funding to develop the livestock sector. Africa is home to about 1.3 billion people and with increasing demand for animal proteins by an ever growing human population, the current state of livestock productivity creates a significant yield gap for animal products. Although a greater section of the population, especially those living in rural areas depend largely on livestock for their livelihoods; the potential of the sector remains underutilized and therefore unable to contribute significantly to economic development and social wellbeing of the people. With current advances in livestock management practices, breeding technologies and health management, and with inclusion of all stakeholders, African livestock populations can be sustainably developed to close the animal protein gap that exists in the continent. In particular, advances in gene technologies, and application of genomic breeding in many Western countries has resulted in tremendous gains in traits like milk production with the potential that, implementation of genomic selection and other improved practices (nutrition, healthcare, etc.) can lead to rapid improvement in traits of economic importance in African livestock populations. The African livestock populations in the context of this review are limited to cattle, goat, pig, poultry, and sheep, which are mainly exploited for meat, milk, and eggs. This review examines the current state of livestock productivity in Africa, the main challenges faced by the sector, the role of various stakeholders and discusses in-depth strategies that can enable the application of genomic technologies for rapid improvement of livestock traits of economic importance.

Comparative transcriptome analysis of mammary epithelial cells at different stages of lactation reveals wide differences in gene expression and pathways regulating milk synthesis between Jersey and Kashmiri cattle

Jersey and Kashmiri cattle are important dairy breeds that contribute significantly to the total milk production of the Indian northern state of Jammu and Kashmir. The Kashmiri cattle germplasm has been extensively diluted through crossbreeding with Jersey cattle with the goal of enhancing its milk production ability. However, crossbred animals are prone to diseases resulting to unsustainable milk production. This study aimed to provide a comprehensive transcriptome profile of mammary gland epithelial cells at different stages of lactation and to find key differences in genes and pathways regulating milk traits between Jersey and Kashmiri cattle. Mammary epithelial cells (MEC) isolated from milk obtained from six lactating cows (three Jersey and three Kashmiri cattle) on day 15 (D15), D90 and D250 in milk, representing early, mid and late lactation, respectively were used. RNA isolated from MEC was subjected to next-generation RNA sequencing and bioinformatics processing. Casein and whey protein genes were found to be highly expressed throughout the lactation stages in both breeds. Largest differences in differentially expressed genes (DEG) were between D15 vs D90 (1,805 genes) in Kashmiri cattle and, D15 vs D250 (3,392 genes) in Jersey cattle. A total of 1,103, 1,356 and 1,397 genes were differentially expressed between Kashmiri and Jersey cattle on D15, D90 and D250, respectively. Antioxidant genes like RPLPO and RPS28 were highly expressed in Kashmiri cattle. Differentially expressed genes in both Kashmiri and Jersey were enriched for multicellular organismal process, receptor activity, catalytic activity, signal transducer activity, macromolecular complex and developmental process gene ontology terms. Whereas, biological regulation, endopeptidase activity and response to stimulus were enriched in Kashmiri cattle and, reproduction and immune system process were enriched in Jersey cattle. Most of the pathways responsible for regulation of milk production like JAK-STAT, p38 MAPK pathway, PI3 kinase pathway were enriched by DEG in Jersey cattle only. Although Kashmiri has poor milk production efficiency, the present study suggests possible physicochemical and antioxidant properties of Kashmiri cattle milk that needs to be further explored.

Genetic parameters of milk cholesterol content in Holstein cattle

This study aimed to estimate heritability for milk cholesterol (CHL) and genetic correlations between milk CHL and other production traits (test-day milk, fat, and protein yields, fat and protein percentages, and somatic cell score). Milk CHL content was determined by gas chromatography and expressed as mg of CHL in 100 g of fat (CHL_fat) or in 100 mg of milk (CHL_milk). Univariate models were used to estimate variances and heritability, whereas bivariate models were used to compute correlations using data from 1793 cows. The average concentrations (standard deviation) of CHL_fat and CHL_milk were 275.63 (75) mg and 11.16 (3.63) mg, respectively. Milk CHL content was significantly affected by days in milk and herd (P < 0.05), but not by parity, regardless of the scale of expression. Heritability estimates for CHL_fat and CHL_milk were 0.06 ± 0.04 and 0.17 ± 0.06, respectively. Phenotypic and genetic correlations between CHL_fat and CHL_milk were 0.82 and 0.44 ± 0.24, respectively. CHL_fat had nonsignificant genetic correlations with all production traits, whereas CHL_milk had significant (P < 0.05) genetic correlations with milk yield (−0.47), fat yield (0.51), protein percentage (0.56), and fat percentage (0.88). This is the first study to estimate genetic parameters for milk CHL content. Further studies are required to assess the possibility of genetically selecting cows with lower milk CHL content.

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.

Integration of miRNA and mRNA Co-Expression Reveals Potential Regulatory Roles of miRNAs in Developmental and Immunological Processes in Calf Ileum during Early Growth

This study aimed to investigate the potential regulatory roles of miRNAs in calf ileum developmental transition from the pre- to the post-weaning period. For this purpose, ileum tissues were collected from eight calves at the pre-weaning period and another eight calves at the post-weaning period and miRNA expression characterized by miRNA sequencing, followed by functional analyses. A total of 388 miRNAs, including 81 novel miRNAs, were identified. A total of 220 miRNAs were differentially expressed (DE) between the two periods. The potential functions of DE miRNAs in ileum development were supported by significant enrichment of their target genes in gene ontology terms related to metabolic processes and transcription factor activities or pathways related to metabolism (peroxisomes), vitamin digestion and absorption, lipid and protein metabolism, as well as intracellular signaling. Integration of DE miRNAs and DE mRNAs revealed several DE miRNA-mRNA pairs with crucial roles in ileum development (bta-miR-374a-FBXO18, bta-miR-374a-GTPBP3, bta-miR-374a-GNB2) and immune function (bta-miR-15b-IKBKB). This is the first integrated miRNA-mRNA analysis exploring the potential roles of miRNAs in calf ileum growth and development during early life.

Genome wide association study identifies novel potential candidate genes for bovine milk cholesterol content

This study aimed to identify single nucleotide polymorphisms (SNPs) associated with milk cholesterol (CHL) content via a genome wide association study (GWAS). Milk CHL content was determined by gas chromatography and expressed as mg of CHL in 100 g of fat (CHL_fat) or in 100 mg of milk (CHL_milk). GWAS was performed with 1,183 cows and 40,196 SNPs using a univariate linear mixed model. Two and 20 SNPs were significantly associated with CHL_fat and CHL_milk, respectively. The important regions for CHL_fat and CHL_milk were at 41.9 Mb on chromosome (BTA) 17 and 1.6-3.2 Mb on BTA 14, respectively. DGAT1, PTPN1, INSIG1, HEXIM1, SDS, and HTR5A genes, also known to be associated with human plasma CHL phenotypes, were identified as potential candidate genes for bovine milk CHL. Additional new potential candidate genes for milk CHL were RXFP1, FAM198B, TMEM144, CXXC4, MAML2 and CDH13. Enrichment analyses suggested that identified candidate genes participated in cell-cell signaling processes and are key members in tight junction, focal adhesion, Notch signaling and glycerolipid metabolism pathways. Furthermore, identified transcription factors such as PPARD, LXR, and NOTCH1 might be important in the regulation of bovine milk CHL content. The expression of several positional candidate genes (such as DGAT1, INSIG1 and FAM198B) and their correlation with milk CHL content were further confirmed with RNA sequence data from mammary gland tissues. This is the first GWAS on bovine milk CHL. The identified markers and candidate genes need further validation in a larger cohort for use in the selection of cows with desired milk CHL content.

Saccharomyces cerevisiae var. boulardii CNCM I-1079 and Lactobacillus acidophilus BT1386 influence innate immune response and serum levels of acute-phase proteins during weaning in Holstein calves

The aims of this study were to investigate the effect of Saccharomyces cerevisiae var. boulardii CNCM I-1079 (SCB) or Lactobacillus acidophilus BT1386 (LA) on (1) innate immune response, (2) markers of acute-phase reaction, and (3) immune gene expression of rumen and ileum tissues of Holstein calves. Forty eight calves (∼5 d old) were randomly allocated to four treatments as follows: (1) control (CTRL) fed milk replacer followed by starter feed, (2) CTRL supplemented with SCB in milk and feed, (3) CTRL supplemented with LA in milk and feed, and (4) CTRL supplemented with antibiotics (ATB; chlortetracycline and neomycin in milk, and chlortetracycline in feed). Tumor necrosis factor α (TNF-α) decreased (P < 0.05) on day 66 (post-weaning) for the ATB-treated calves. There were no treatment effects on production of interferon γ (IFN-γ) and interleukin 6 (IL-6) proteins and on expression of TLR4, TLR6, TLR9, TLR10, CLDN3, MUC1, and MUC20 genes. Calves fed SCB or LA had a greater (P < 0.05) oxidative burst at weaning (day 53) compared with CTRL. Oxidative burst was also greater (P < 0.05) after weaning (day 59 and day 87) for SCB-fed calves. Calves fed SCB and ATB had higher (P < 0.05) phagocytosis activity during weaning (day 47) compared with CTRL. The concentration of serum amyloid A2 (SAA2) increased (P < 0.05) in SCB- and LA-fed calves (day 53), whereas the concentration of C-reactive protein (CRP) increased (P < 0.05) in SCB-fed calves during weaning as compared with CTRL. Our results suggest that SCB could improve innate immune response (oxidative burst and phagocytosis) and markers of acute-phase reaction (CRP and SAA2), especially during critical periods like weaning.

Integration of miRNA weighted gene co-expression network and miRNA-mRNA co-expression analyses reveals potential regulatory functions of miRNAs in calf rumen development

This study aimed to explore the roles of microRNAs (miRNAs) in calf rumen development during early life. Rumen tissues were collected from 16 calves (8 at pre-weaning and 8 at post-weaning) for miRNA-sequencing, differential expression (DE), miRNA weighted gene co-expression network (WGCNA) and miRNA-mRNA co-expression analyses. 295 miRNAs were identified. Bta-miR-143, miR-26a, miR-145 and miR-27b were the most abundantly expressed. 122 miRNAs were significantly DE between the pre- and post-weaning periods and the most up- and down-regulated miRNAs were bta-miR-29b and bta-miR-493, respectively. Enrichment analyses of the target genes of DE miRNAs revealed important roles for miRNA in rumen developmental processes, immune system development, protein digestion and processes related to the extracellular matrix. WGCNA indicated that bta-miR-145 and bta-miR-199a-3p are important hub miRNAs in the regulation of these processes. Therefore, bta-miR-143, miR-29b, miR-145, miR-493, miR-26a and miR-199 family members might be key regulators of calf rumen development during early life.

Polymorphisms of caprine GnRHR gene and their association with litter size in West African Dwarf goats

Gonadotropin-releasing hormone receptor (GnRHR) gene is considered a candidate gene for litter size due to its critical role in regulating the activities of hypothalamo-pituitary-gonadal axis which synthesizes and releases gonadotropins. This study was designed to identify mutations within the caprine GnRHR gene and investigate their association with litter size at various parities. Polymorphisms scanning and genotyping of GnRHR gene in West African Dwarf (WAD) goats (n = 226) revealed three single nucleotide polymorphisms (SNPs), one mutation (g.-29T > G) was detected within 5'UTR region while two others (g.48G > A and g.209T > G) were identified in exon 1. Mutation at g.209T > G locus resulted in amino acid change from Methionine to Arginine at position 70 on the polypeptide residue. Based on heterozygosity and polymorphism information content (PIC), WAD goat population diversity at the SNP loci was moderate. Strong linkage disequilibrium (LD) (r² > 0.98) existed among the detected mutations resulting in three observed haplotypes, two (T-G-T and G-A-G) had cumulative frequency of > 97%. The mutation within 5'UTR region of GnRHR gene (g.-29T > G) is novel, being reported in goats for the first time. Association analysis revealed a significant (p < 0.05) association between allele G at g.-29T > G with higher mean litter size for homozygous (GG) mutant does compared with heterozygotes (GT) or homozygotes (TT), while the relationship between SNPs at the two loci detected in exon 1 and litter size was not significant.

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.

MicroRNA roles in signalling during lactation: an insight from differential expression, time course and pathway analyses of deep sequence data

The study examined microRNA (miRNA) expression and regulatory patterns during an entire bovine lactation cycle. Total RNA from milk fat samples collected at the lactogenesis (LAC, day1 [D1] and D7), galactopoiesis (GAL, D30, D70, D130, D170 and D230) and involution (INV, D290 and when milk production dropped to 5 kg/day) stages from 9 cows was used for miRNA sequencing. A total of 475 known and 238 novel miRNAs were identified. Fifteen abundantly expressed miRNAs across lactation stages play regulatory roles in basic metabolic, cellular and immunological functions. About 344, 366 and 209 miRNAs were significantly differentially expressed (DE) between GAL and LAC, INV and GAL, and INV and LAC stages, respectively. MiR-29b/miR-363 and miR-874/miR-6254 are important mediators for transition signals from LAC to GAL and from GAL to INV, respectively. Moreover, 58 miRNAs were dynamically DE in all lactation stages and 19 miRNAs were significantly time-dependently DE throughout lactation. Relevant signalling pathways for transition between lactation stages are involved in apoptosis (PTEN and SAPK/JNK), intracellular signalling (protein kinase A, TGF-β and ERK5), cell cycle regulation (STAT3), cytokines, hormones and growth factors (prolactin, growth hormone and glucocorticoid receptor). Overall, our data suggest diverse, temporal and physiological signal-dependent regulatory and mediator functions for miRNAs during 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.