MicroRNA profiling in bovine serum according to the stage of Mycobacterium avium subsp. paratuberculosis infection

Towards a microRNA-based Johne's disease diagnostic predictive system: Preliminary results

BACKGROUND

Johne's disease, caused by Mycobacterium avium subspecies paratuberculosis (MAP), is a chronic enteritis that adversely affects welfare and productivity in cattle. Screening and subsequent removal of affected animals is a common approach for disease management, but efforts are hindered by low diagnostic sensitivity. Expression levels of small non-coding RNA molecules involved in gene regulation (microRNAs), which may be altered during mycobacterial infection, may present an alternative diagnostic method.

METHODS

The expression levels of 24 microRNAs affected by mycobacterial infection were measured in sera from MAP-positive (n = 66) and MAP-negative cattle (n = 65). They were then used within a machine learning approach to build an optimal classifier for MAP diagnosis.

RESULTS

The method provided 72% accuracy, 73% sensitivity and 71% specificity on average, with an area under the curve of 78%.

LIMITATIONS

Although control samples were collected from farms nominally MAP-free, the low sensitivity of current diagnostics means some animals may have been misclassified.

CONCLUSION

MicroRNA profiling combined with advanced predictive modelling enables rapid and accurate diagnosis of Johne's disease in cattle.

Extensive differential DNA methylation between tuberculosis skin test positive and skin test negative cattle

Bovine tuberculosis (bTB), caused by Mycobacterium bovis (M. bovis), represents a significant problem for the agriculture industry as well as posing a risk for human health. Current diagnostic tests for bTB target the cell-mediated immune (CMI) response to infection with M. bovis, primarily through screening of animals with the tuberculin skin test. Epigenetic modifications have been shown to alter the course of the immune response and differentially methylated regions (DMRs) might also influence the outcome of the skin test in cattle. Whole Genome Bisulphite Sequencing (WGBS) was used to profile DNA methylation levels from peripheral blood of a group of cattle identified as test positive for M. bovis (positive for the single intradermal comparative tuberculin test (SICTT) and/or the interferon-γ release assay compared to a test negative control group [n = 8/group, total of 16 WGBS libraries]. Although global methylation profiles were similar for both groups across the genome, 223 DMRs and 159 Differentially Promoter Methylated Genes (DPMGs) were identified between groups with an excess of hypermethylated sites in SICTT positive cattle (threshold > 15% differential methylation). Genes located within these DMRs included the Interleukin 1 receptor (IL1R1) and MHC related genes (BOLA and BOLA-DQB). KEGG pathway analysis identified enrichment of genes involved in Calcium and MAPK signalling, as well as metabolism pathways. Analysis of DMRs in a subset of SICTT negative cattle that were IFN-γ positive showed differential methylation of genes including Interleukin 10 Receptor, alpha (IL10RA), Interleukin 17 F (IL17F) and host defence peptides (DEFB and BDEF109). This study has identified a number of immune gene loci at which differential methylation is associated with SICTT test results and the degree of methylation could influence effective host immune responses.

Integrated Analysis of Transcriptome Profiles and lncRNA-miRNA-mRNA Competing Endogenous RNA Regulatory Network to Identify Biological Functional Effects of Genes and Pathways Associated with Johne's Disease in Dairy Cattle

Paratuberculosis or Johne's disease (JD), a chronic granulomatous gastroenteritis caused by Mycobacterium avium subsp. paratuberculosis (MAP), causes huge economic losses and reduces animal welfare in dairy cattle herds worldwide. At present, molecular mechanisms and biological functions involved in immune responses to MAP infection of dairy cattle are not clearly understood. Our purpose was to integrate transcriptomic profiles and competing endogenous RNA (ceRNA) network analyses to identify key messenger RNAs (mRNAs) and regulatory RNAs involved in molecular regulation of peripheral blood mononuclear cells (PBMCs) for MAP infection in dairy cattle. In total, 28 lncRNAs, 42 miRNAs, and 370 mRNAs were identified by integrating gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. In this regard, we identified 21 hub genes (CCL20, CCL5, CD40, CSF2, CXCL8, EIF2AK2, FOS, IL10, IL17A, IL1A, IL1B, IRF1, MX2, NFKB1, NFKBIA, PTGS2, SOCS3, TLR4, TNF, TNFAIP3, and VCAM1) involved in MAP infection. Furthermore, eight candidate subnets with eight lncRNAs, 29 miRNAs, and 237 mRNAs were detected through clustering analyses, whereas GO enrichment analysis of identified RNAs revealed 510, 22, and 11 significantly enriched GO terms related to MAP infection in biological process, molecular function, and cellular component categories, respectively. The main metabolic-signaling pathways related to MAP infection that were enriched included the immune system process, defense response, response to cytokine, leukocyte migration, regulation of T cell activation, defense response to bacterium, NOD-like receptor, B cell receptor, TNF, NF-kappa B, IL-17, and T cell receptor signaling pathways. Contributions of transcriptome profiles from MAP-positive and MAP-negative sample groups plus a ceRNA regulatory network underlying phenotypic differences in the intensity of pathogenicity of JD provided novel insights into molecular mechanisms associated with immune system responses to MAP infection in dairy cattle.

MicroRNAs modulate immunological and inflammatory responses in Holstein cattle naturally infected with Mycobacterium avium subsp. paratuberculosis

MicroRNAs (miRNAs) regulate the post-transcriptional expression of genes by binding to their target mRNAs. In this study, whole miRNA sequencing was used to compare the expression of miRNAs in ileocecal valve (ICV) and peripheral blood (PB) samples of cows with focal or diffuse paratuberculosis (PTB)-associated lesions in gut tissues versus (vs) control cows without lesions. Among the eight miRNAs differentially expressed in the PB samples from cows with diffuse lesions vs controls, three (miR-19a, miR-144, miR32) were also down-regulated in cows with diffuse vs focal lesions. In the ICV samples, we identified a total of 4, 5, and 18 miRNAs differentially expressed in cows with focal lesions vs controls, diffuse lesions vs controls, and diffuse vs focal lesions, respectively. The differential expression of five microRNAs (miR-19a, miR-144, miR-2425-3p, miR-139, miR-101) was confirmed by RT-qPCR. Next, mRNA target prediction was performed for each differentially expressed miRNA. A functional analysis using the predicted gene targets revealed a significant enrichment of the RNA polymerase and MAPK signaling pathways in the comparison of cows with focal vs no lesions and with diffuse vs focal lesions, respectively. The identified miRNAs could be used for the development of novel diagnostic and therapeutical tools for PTB control.

Association of plasma miRNAs with early life performance and aging in dairy cattle

Early life performance traits in dairy cattle can have important influences on lifetime productivity. Poor health and fertility are of great economical and animal welfare concern. Circulating miRNAs have been linked to several livestock traits, including resistance to infection, fertility, and muscle development. This study aimed to identify circulating miRNAs associated with early life performance traits and aging in dairy cattle. Plasma samples from female calves (n = 12) identified retrospectively as differing in health, growth, and fertility outcomes prior to first calving were analyzed using PCR arrays detecting 378 miRNAs. Levels of 6 miRNAs differed significantly in calves with poor growth/fertility relative to controls (t-test: P<0.05). Additionally, general(ized) (non)linear mixed models identified 1 miRNA associated with average daily gain until weaning, 22 with live bodyweight at one year of age, 47 with age at first service, and 19 with number of infections before first calving. Out of 85 distinct miRNAs that were associated with at least one animal trait, 9 miRNAs were validated by RT-qPCR in a larger cohort (n = 91 animals), which included longitudinal plasma samples (calf, heifer, first lactation cow). Significant associations (P<0.05) involving individual miRNAs or ratios between miRNAs and early-life performance traits were identified, but did not retain significance after multiple testing adjustment. However, levels of 8 plasma miRNAs (miR-126-3p, miR-127, miR-142-5p, miR-154b, miR-27b, miR-30c-5p, miR-34a, miR-363) changed significantly with age, most prominently during the calf-to-heifer transition. Comparative RT-qPCR analyses of these miRNAs across 19 calf tissues showed that most were ubiquitously expressed. Online database mining identified several pathways involved in metabolism and cell signaling as putative biological targets of these miRNAs. These results suggest that miR-126-3p, miR-127, miR-142-5p, miR-154b, miR-27b, miR-30c-5p, miR-34a, miR-363 are involved in regulating growth and development from birth to first lactation (~2 years old) and could provide useful biomarkers of aging in cattle.

Identification of Suitable Internal Control miRNAs in Bovine Milk Small Extracellular Vesicles for Normalization in Quantitative Real-Time Polymerase Chain Reaction

This study aimed to identify a suitable RNA extraction kit and stable internal control microRNA (miRNA) in bovine milk small extracellular vesicles (sEVs) for a quantitative polymerase chain reaction (qPCR) analysis. Two RNA extraction kits, miRNeasy Micro Kit, and Maxwell RSC miRNA Tissue Kit, were compared and evaluated using bovine milk sEVs via qPCR analysis. Five miRNAs, bta-miR-29a, bta-miR-200a, bta-miR-26b, hsa-miR-27b-3p, and hsa-miR-30b-5p, were selected by microarray analyses, and their cycle threshold (Ct) values were further evaluated mathematically using geNorm, NormFinder, BestKeeper, and ∆Ct algorithms. The results revealed that both the miRNeasy Micro Kit and Maxwell RSC miRNA Tissue Kit are useful for the efficient recovery of RNA from bovine milk sEVs. According to the final stability ranking analyzed by RefFinder, hsa-miR-27b-3p and bta-miR-29a can be used as suitable internal control miRNAs in bovine milk sEVs. The study also indicated that using a suitable internal control miRNA may improve the reliability and accuracy of the qPCR analysis for normalization in bovine milk sEVs. To the best of our knowledge, this is the first study to uncover the suitable internal control miRNAs in bovine milk sEVs.

Fetuin as a potential serum biomarker to detect subclinical shedder of bovine paratuberculosis

Paratuberculosis (PTB) is a chronic contagious granulomatous enteritis of wild and domestic ruminants caused by Mycobacterium avium subsp. paratuberculosis (MAP). PTB causes considerable economic losses to the dairy industry through decreased milk production and premature culling. PTB-affected cattle undergo a subclinical stage without clinical signs and initiate fecal shedding of MAP into the environment. Current diagnostic tools have low sensitivity for the detection of subclinical PTB infection. Therefore, alternative diagnostic tools are required to improve the diagnostic sensitivity of subclinical PTB infection. In this study, we performed ELISA for three previously identified host biomarkers (fetuin, alpha-1-acid glycoprotein, and apolipoprotein) and analyzed their diagnostic performance with conventional PTB diagnostic methods. We observed that serum fetuin levels were significantly lowered in the subclinical shedder and clinical shedder groups than in the healthy control group, indicating its potential utility as a diagnostic biomarker for bovine PTB. Also, fetuin showed an excellent discriminatory power with an AUC = 0.949, a sensitivity of 92.6%, and a specificity of 94.4% for the detection of subclinical MAP infection. In conclusion, our results demonstrated that fetuin could be used as a diagnostic biomarker for enhancing the diagnostic sensitivity for the detection of subclinical MAP infections that are difficult to detect based on current diagnostic methods.