Genetic association analysis of paratuberculosis forms in holstein-friesian cattle

Summary-data based Mendelian randomization identifies gene expression regulatory polymorphisms associated with bovine paratuberculosis by modulation of the nuclear factor Kappa β (NF-κß)-mediated inflammatory response

Genome-wide association studies (GWAS) have identified host genetic variants associated with paratuberculosis (PTB) susceptibility. Most of the GWAS-identified SNPs are in non-coding regions. Connecting these non-coding variants and downstream affected genes is a challenge and, up to date, only a few functional mutations or expression quantitative loci (cis-eQTLs) associated with PTB susceptibility have been identified. In the current study, the associations between imputed whole-genome sequence genotypes and whole RNA-Sequencing data from peripheral blood (PB) and ileocecal valve (ICV) samples of Spanish Holstein cows (N = 16) were analyzed with TensorQTL. This approach allowed the identification of 88 and 37 cis-eQTLs regulating the expression levels of 90 and 37 genes in PB and ICV samples, respectively (False discorey rate, FDR ≤ 0.05). Next, we applied summary-based data Mendelian randomization (SMR) to integrate the cis-eQTL dataset with GWAS data obtained from a cohort of 813 culled cattle that were classified according to the presence or absence of PTB-associated histopathological lesions in gut tissues. After multiple testing corrections (FDR ≤ 0.05), we identified two novel cis-eQTLs affecting the expression of the early growth response factor 4 (EGR4) and the bovine neuroblastoma breakpoint family member 6-like protein isoform 2 (MGC134040) that showed pleiotropic associations with the presence of multifocal and diffuse lesions in gut tissues; P = 0.002 and P = 0.017, respectively. While EGR4 acts as a brake on T-cell proliferation and cytokine production through interaction with the nuclear factor Kappa β (NF-κß), MGC134040 is a target gene of NF-κß. Our findings provide a better understanding of the genetic factors influencing PTB outcomes, confirm that the multifocal lesions are localized/confined lesions that have different underlying host genetics than the diffuse lesions, and highlight regulatory SNPs and regulated-gene targets to design future functional studies.

Early-stage findings in an experimental calf model infected with Argentinean isolates of Mycobacterium avium subsp. paratuberculosis

Mycobacterium avium subsp. paratuberculosis (MAP) is an important pathogen that causes granulomatous enteritis known as Johne's disease or paratuberculosis (PTB). In this study an experimental model of calves infected with Argentinean isolates of MAP for 180 days was used to provide more data of the early PTB stages. Calves were challenged by oral route with MAP strain _IS900-RFLPA (MA; n = 3), MAP strain _IS900-RFLPC (MC; n = 2) or mock infected (MI; n = 2), and response to infection was evaluated through peripheral cytokine expression, MAP tissue distribution and histopathological early-stage findings. Specific and varied levels of IFN-γ were only detected at 80 days post-infection in infected calves. These data indicate that specific IFN-γ is not a useful indicator for early detection of MAP infection in our calf model. At 110 days post-infection, TNF-α expression was higher than IL-10 in 4 of the 5 infected animals and a significant decrease of TNF-α expression was detected in infected vs. non-infected calves. All calves challenged were identified as infected by mesenteric lymph node tissue culture and real time IS900 PCR. In addition, for lymph nodes samples, the agreement between these techniques was almost perfect (κ = 0.86). Colonization of tissues and levels of tissue infection varied between individuals. Evidence of early MAP dissemination to extraintestinal tissues such as the liver was detected by culture in one animal (MAP strain _IS900-RFLPA). In both groups microgranulomatous lesions were observed predominantly in the lymph nodes, with giant cells present only in the MA group. In summary, the findings described herein may indicate that local MAP strains induced specific immune responses with particularities that could suggest differences in their biological behavior. Further studies should be carried out in order to obtain an in-depth understanding of the influence of MAP strains in host-pathogen interactions and the outcome of disease.

The Upregulation of Cathepsin G Is Associated with Resistance to Bovine Paratuberculosis

An in silico genomic-transcriptomic combined approach allowed the identification of a polymorphism (cis-eQTL-rs41976219) in the Bos taurus genome associated with the CTSG mRNA expression in bovine blood samples, which suggests that individual genetic variation might modulate the CTSG transcriptional response. In the current study, a sandwich ELISA is used to measure the CTSG protein levels in supernatants of monocyte-derived macrophages (MDMs) isolated from cows with the AA (n = 5) and AC (n = 11) genotypes for the rs41976219 and infected ex vivo with MAP. Cows with the AC genotype have significantly higher CTSG protein levels (1.85 ng/mL) in the supernatants of enriched CD14+-MDMs after 2 h of infection when compared with infected CD14+-MDMs from cows with the AA genotype (1.68 ng/mL). Statistically significant differences in the intracellular MAP load at 7 d p.i. are observed between animals with the AA (2.16 log CFUs) and AC (1.44 log CFUs) genotypes. Finally, the association between the rs41976219 allelic variants and resistance to PTB is tested in a larger cattle population (n = 943) classified according to the presence (n = 442) or absence (n = 501) of PTB-associated lesions. The presence of the two minor alleles in the rs41976219 (CC) is more frequent among healthy cows than in cows with PTB-associated lesions in gut tissues (2.2% vs. 1.4%, OR = 0.61). In agreement with this, the CTSG levels in plasma samples of cows without lesions in gut tissues and with the CC (n = 8) genotype are significantly higher than in the plasmas of cows with the AA + AC (n = 36) genotypes.

A Genome-Wide Association Study for Tolerance to Paratuberculosis Identifies Candidate Genes Involved in DNA Packaging, DNA Damage Repair, Innate Immunity, and Pathogen Persistence

Although the genetic susceptibility to diseases has been extensively studied, the genetic loci and the primary molecular and cellular mechanisms that control disease tolerance are still largely unknown. Bovine paratuberculosis (PTB) is an enteritis caused by Mycobacterium avium subsp. paratuberculosis (MAP). PTB affects cattle worldwide and represents a major issue on animal health. In this study, the associations between host genetic and PTB tolerance were investigated using the genotypes from 277 Spanish Holstein cows with two distinct phenotypes: cases) infected animals with positive PCR and bacteriological culture results but without lesions in gut tissues (N= 24), and controls) animals with negative PCR and culture results but with PTB-associated lesions (N= 253). DNA from peripheral blood of the study population was genotyped with the Bovine EuroG MD Bead Chip, and the corresponding genotypes were imputed to whole-genome sequencing (WGS) data. A genome-wide association study was performed using the WGS data and the defined phenotypes in a case-control approach. A total of 142 single nucleotide polymorphisms (SNPs) were associated (false discovery rate ≤ 0.05, P values between 1.5 × 10^-7 and 5.7 × 10^-7) with tolerance (heritability= 0.55). The 40 SNPs with P-values < 5 × 10^-7 defined 9 QTLs and 98 candidate genes located on BTA4, BTA9, BTA16, BTA25, and BTA26. Some of the QTLs identified in this study overlap with QTLs previously associated with PTB, bovine tuberculosis, mastitis, somatic cell score, bovine diarrhea virus persistent infection, tick resistance, and length of productive life. Two candidate genes with important roles in DNA damage response (ERCC4 and RMI2) were identified on BTA25. Functional analysis using the 98 candidate genes revealed a significant enrichment of the DNA packaging process (TNP2/PRMI1/PRM2/PRM3). In addition, the TNF-signaling (bta04668; TRAF5/CREB5/CASP7/CHUK) and the toxoplasmosis (bta05145; TGFβ2/CHUK/CIITA/SOCS1) pathways were significantly enriched. Interestingly, the nuclear Factor NF-κβ Inhibitor Kinase Alpha (CHUK), a key molecule in the regulation of the NF-κB pathway, was enriched in both pathways. Taken together, our results define a distinct immunogenetic profile in the PTB-tolerant animals designed to control bacterial growth, modulate inflammation, limit tissue damage and increase repair, thus reducing the severity of the disease.

Use of ATP-Binding Cassette Subfamily A Member 13 (ABCA13) for Sensitive Detection of Focal Pathological Forms of Subclinical Bovine Paratuberculosis

Bovine paratuberculosis (PTB) is a chronic enteritis caused by Mycobacterium avium subspecies paratuberculosis (Map) that causes a heavy economic impact worldwide. Map infected animals can remain asymptomatic for years while transmitting the mycobacteria to other members of the herd. Therefore, accurate detection of subclinically infected animals is crucial for disease control. In a previous RNA-Seq study, we identified several mRNAs that were overexpressed in whole blood of cows with different PTB-associated histological lesions compared with control animals without detected lesions. The proteins encoded by two of these mRNAs, ATP binding cassette subfamily A member 13 (ABCA13) and Matrix Metallopeptidase 8 (MMP8) were significantly overexpressed in whole blood of animals with focal histological lesions, the most frequent pathological form in the subclinical stages of the disease. In the current study, the potential of sensitive early diagnostic tools of commercial ELISAs, based on the detection of these two biomarkers, was evaluated in serum samples of 704 Holstein Friesian cows (566 infected animals and 138 control animals from PTB-free farms). For this evaluation, infected animals were classified into three groups, according to the type of histological lesions present in their gut tissues: focal (n = 447), multifocal (n = 59), and diffuse (n = 60). The ELISA based on the detection of ABCA13 was successfully validated showing good discriminatory power between animals with focal lesions and control animals (sensitivity 82.99% and specificity 80.43%). Conversely, the MMP8-based ELISA showed a poor discriminatory power between the different histological groups and non-infected controls. The ABCA13-based ELISA showed a higher diagnostic value (0.822) than the IDEXX ELISA (0.517), the fecal bacterial isolation (0.523) and the real-time PCR (0.531) for the detection of animals with focal lesions. Overall, our results indicate that this ABCA13 ELISA greatly improves the identification of subclinically infected animals with focal lesions that are undetectable using current diagnostic methods.

Effect of selected single nucleotide polymorphisms in SLC11A1, ANKRA2, IFNG and PGLYRP1 genes on host susceptibility to Mycobacterium avium subspecies paratuberculosis infection in Indian cattle

Paratuberculosis (PTB) is a chronic infectious enteritis of ruminants, caused by Mycobacterium avium subspecies paratuberculosis (MAP) that brings huge economic loss to the dairy farmers. The study was conducted to explore the association of selected SNPs in IFNG, SLC11A1, ANKRA2 and PGLYRP1 genes with resistance to PTB disease in Indian cattle population. A case-control resource population was established based on the results of diagnostic tests used for detection of MAP infection status viz. ELISA, Johnin PPD test, faecal microscopy and IS900 blood PCR. The PCR-RFLP method was used for genotyping of SNPs. SNPs rs109453173 in SLC11A1, rs110853455 in IFNG and rs41933863 in ANKRA2 genes were significantly (P<0.05) associated with resistance to MAP infection. For SNP rs109453173, GG genotype and G allele was found to be associated with resistance against MAP infection than CC and CG genotypes and C allele, respectively. For SNP rs110853455, AG genotype was found to be associated with susceptibility to MAP infection than AA and GG genotype. For SNP rs41933863, the AG genotype provided three and six times more resistance against MAP infection than GG and AA genotype. The results of this study are suggestive of SNPs rs109453173, rs110853455 and rs41933863 as potential markers for screening MAP resistant cattle and a breeding programme favouring GG genotype and G allele for rs109453173, AG genotype for rs41933863 and against AG genotype for rs110853455 might confer resistance against MAP infection in Indian cattle. However, investigation of these SNPs in an independent and larger population will warrant the strength of association for resistance against MAP infection in cattle.

Identification of loci associated with pathological outcomes in Holstein cattle infected with Mycobacterium avium subsp. paratuberculosis using whole-genome sequence data

Bovine paratuberculosis (PTB), caused by Mycobacterium avium subsp. paratuberculosis (MAP), is a chronic granulomatous enteritis that affects cattle worldwide. According to their severity and extension, PTB-associated histological lesions have been classified into the following groups; focal, multifocal, and diffuse. It is unknown whether these lesions represent sequential stages or divergent outcomes. In the current study, the associations between host genetic and pathology were explored by genotyping 813 Spanish Holstein cows with no visible lesions (N = 373) and with focal (N = 371), multifocal (N = 33), and diffuse (N = 33) lesions in gut tissues and regional lymph nodes. DNA from peripheral blood samples of these animals was genotyped with the bovine EuroG MD Bead Chip, and the corresponding genotypes were imputed to whole-genome sequencing (WGS) data using the 1000 Bull genomes reference population. A genome-wide association study (GWAS) was performed using the WGS data and the presence or absence of each type of histological lesion in a case–control approach. A total of 192 and 92 single nucleotide polymorphisms (SNPs) defining 13 and 9 distinct quantitative trait loci (QTLs) were highly-associated (P ≤ 5 × 10⁻⁷) with the multifocal (heritability = 0.075) and the diffuse (heritability = 0.189) lesions, respectively. No overlap was seen in the SNPs controlling these distinct pathological outcomes. The identified QTLs overlapped with some QTLs previously associated with PTB susceptibility, bovine tuberculosis susceptibility, clinical mastitis, somatic cell score, bovine respiratory disease susceptibility, tick resistance, IgG level, and length of productive life. Pathway analysis with candidate genes overlapping the identified QTLs revealed a significant enrichment of the keratinization pathway and cholesterol metabolism in the animals with multifocal and diffuse lesions, respectively. To test whether the enrichment of SNP variants in candidate genes involved in the cholesterol metabolism was associated with the diffuse lesions; the levels of total cholesterol were measured in plasma samples of cattle with focal, multifocal, or diffuse lesions or with no visible lesions. Our results showed reduced levels of plasma cholesterol in cattle with diffuse lesions. Taken together, our findings suggested that the variation in MAP-associated pathological outcomes might be, in part, genetically determined and indicative of distinct host responses.

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

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

Genomic Association between SNP Markers and Diseases in the "Curraleiro Pé-Duro" Cattle

Susceptibility to diseases is inherited and can be transmitted between populations. Single-nucleotide polymorphism (SNPs) in genes related to immune response is associated with diseases in cattle. This study investigated SNPs in the genomic region of cytokines in 702 samples of Curraleiro Pé-Duro cattle and associated them with the occurrence of antibodies in brucellosis, leptospirosis, neosporosis, leukosis, infectious bovine rhinotracheitis (IBR), and bovine viral diarrhea (BVD) tests. DNA samples were evaluated by the kompetitive allele-specific polymerase chain reaction (KASP) method to identify polymorphisms. The gametic phase and SNP haplotypes were determined with the help of PHASE 2.1.1 software. Haplotypes were associated with serological results against Brucella abortus, Leptospira sp., Neospora caninum, leukosis, infectious rhinotracheitis, and BVD using univariate analysis followed by logistic regression. Haplotype 2 of TLR2 was present in 70% of the animals that tested positive for N. caninum infection. Haplotypes of TLR10 and TLR6 and IL10RA were more common in seronegative animals. Haplotypes related to the gene IL10RA were associated with animals negative to all infections. Curraleiro Pé-Duro cattle presented polymorphisms related to resistance to bacterial, viral, and N. caninum infections.

Bovine Intelectin 2 Expression as a Biomarker of Paratuberculosis Disease Progression

Simple Summary

The potential of the bovine intelectin 2 as a biomarker of Mycobacterium avium subsp. paratuberculosis infection was investigated using quantitative immunohistochemical analysis of ileocecal valve samples of animals with increasing degrees of lesion severity (focal, multifocal and diffuse histological lesions) and control animals without detected lesions. Significant differences were observed in the mean number of intelectin 2 immunolabelled cells between the three histopathological types and the control. Specifically, the mean number of intelectin 2 labelled cells was indicative of disease progression as the focal group had the highest number of intelectin 2 secreting cells followed by the multifocal, diffuse and control groups indicating that intelectin 2 is a good biomarker for the different stages of Mycobacterium avium subsp. paratuberculosis infection. Quantification of bovine intelectin 2 secreting cells could constitute a good post-mortem tool, complementary to histopathology, to improve detection of Mycobacterium avium subsp. Paratuberculosis infections, especially latent forms of infection.

Abstract

Paratuberculosis (PTB), a chronic granulomatous enteritis caused by Mycobacterium avium subsp. paratuberculosis (MAP), is responsible for important economic losses in the dairy industry. Our previous RNA-sequencing (RNA-Seq) analysis showed that bovine intelectin 2 (ITLN2) precursor gene was overexpressed in ileocecal valve (ICV) samples of animals with focal (log2 fold-change = 10.6) and diffuse (log2 fold-change = 6.8) PTB-associated lesions compared to animals without lesions. This study analyzes the potential use of ITLN2, a protein that has been described as fundamental in the innate immune response to infections, as a biomarker of MAP infection. The presence of ITLN2 was investigated by quantitative immunohistochemical analysis of ICV samples of 20 Holstein Friesian cows showing focal (n = 5), multifocal (n = 5), diffuse (n = 5) and no histological lesions (n = 5). Significant differences were observed in the mean number of ITLN2 immunostained goblet and Paneth cells between the three histopathological types and the control. The number of immunolabelled cells was higher in the focal histopathological type (116.9 ± 113.9) followed by the multifocal (108.7 ± 140.5), diffuse (76.5 ± 97.8) and control types (41.0 ± 81.3). These results validate ITLN2 as a post-mortem biomarker of disease progression.

Survey of Candidate Single-Nucleotide Polymorphisms in SLC11A1, TLR4, NOD2, PGLYRP1, and IFNγ in Ankole Longhorn Cattle in Central Region of Uganda to Determine Their Role in Mycobacterium avium Subspecies paratuberculosis Infection Outcome

Mycobacterium avium ssp. paratuberculosis (MAP) is the cause of Johne's disease (JD) in a wide range of domestic and wild ruminants. Single-nucleotide polymorphisms (SNPs) in several genes including solute-like carrier 11A1 (SLC11A1), interferon gamma (IFNγ), Toll-like receptor 4 (TLR4), nucleotide-binding oligomerization domain 2 gene (NOD2), and bovine peptidoglycan recognition protein 1 (PGLYRP1) have been implicated in influencing the infection outcome of MAP in cattle. We have carried out a survey in a population of Ankole cattle from three districts in the central region of Uganda including Isingiro, Lyantonde, and Rakai to determine the role played by several SNPs on the above genes in the infection outcome of local cattle in Uganda. Nine hundred fifty-five heads of cattle obtained from 93 herds were tested using ELISA. Thirty-five ELISA-positive cattle and 35 negative herd mates from a total of 955 cattle tested for MAP were genotyped using iPLEX MassARRAY genotyping systems to detect the presence of a total of 13 SNPS in five different genes (SLC11A1, IFNγ, TLR4, NOD2, and PGLYRP1). The cow-level prevalence of MAP infection in Ankole Longhorn cattle in the three districts was 3.98% (35/955), while the herd-level prevalence was 27.9% and within-herd prevalence was 12 ± 1.5% (95% CI = 9.1-14.8%). The genotypes and allele frequencies of the MAP-positive cattle were compared with those of their ELISA-negative herd mates to determine the significance of the polymorphisms. The results showed that SNPs rs109915208, rs110514940, and rs110905610 on SLC11A1, c.480G>A and c.625C>A on PGLYRP1, and c.2021C>T on TLR4 were monomorphic in both seropositive and seronegative cattle and therefore had no influence on the infection outcome. The remaining SNPs studied in the five genes [SLC11A1: rs109614179; TLR4: rs29017188 (c.226G>C), c.2021C>T; NOD2: rs110536091, rs111009394; PGLYRP1: c.102G>C, c.480G>A, c.625C>A; IFNγ: rs110853455] were polymorphic, but their allele and genotype frequencies did not show any significant difference between the seropositive and seronegative cattle. No significant difference was observed for any haplotype at the gene level.

Identification of loci associated with susceptibility to bovine paratuberculosis and with the dysregulation of the MECOM, eEF1A2, and U1 spliceosomal RNA expression

Although genome-wide association studies have identified single nucleotide polymorphisms (SNPs) associated with the susceptibility to Mycobacterium avium subsp. paratuberculosis (MAP) infection, only a few functional mutations for bovine paratuberculosis (PTB) have been characterized. Expression quantitative trait loci (eQTLs) are genetic variants typically located in gene regulatory regions that alter gene expression in an allele-specific manner. eQTLs can be considered as functional links between genomic variants, gene expression, and ultimately phenotype. In the current study, peripheral blood (PB) and ileocecal valve (ICV) gene expression was quantified by RNA-Seq from fourteen Holstein cattle with no lesions and with PTB-associated histopathological lesions in gut tissues. Genotypes were generated from the Illumina LD EuroG10K BeadChip. The associations between gene expression levels (normalized read counts) and genetic variants were analyzed by a linear regression analysis using R Matrix eQTL 2.2. This approach allowed the identification of 192 and 48 cis-eQTLs associated with the expression of 145 and 43 genes in the PB and ICV samples, respectively. To investigate potential relationships between these cis-eQTLs and MAP infection, a case–control study was performed using the genotypes for all the identified cis-eQTLs and phenotypical data (histopathology, ELISA for MAP-antibodies detection, tissue PCR, and bacteriological culture) of 986 culled cows. Our results suggested that the heterozygous genotype in the cis-eQTL-rs43744169 (T/C) was associated with the up-regulation of the MDS1 and EVI1 complex (MECOM) expression, with positive ELISA, PCR, and bacteriological culture results, and with increased risk of progression to clinical PTB. As supporting evidence, the presence of the minor allele was associated with higher MECOM levels in plasma samples from infected cows and with increased MAP survival in an ex-vivo macrophage killing assay. Moreover, the presence of the two minor alleles in the cis-eQTL-rs110345285 (C/C) was associated with the dysregulation of the eukaryotic elongation factor 1-α2 (eEF1A2) expression and with increased ELISA (OD) values. Finally, the presence of the minor allele in the cis-eQTL rs109859270 (C/T) was associated with the up-regulation of the U1 spliceosomal RNA expression and with an increased risk of progression to clinical PTB. The introduction of these novel functional variants into marker-assisted breeding programs is expected to have a relevant effect on PTB control.

Association of genetic variability in CD209 gene with bovine paratuberculosis disease: a case-control study in the Indian cattle population

This study aimed to explore the association of single nucleotide polymorphisms (SNPs) in CD209 gene with the occurrence of bovine paratuberculosis (PTB) disease caused by Mycobacterium avium subspecies paratuberculosis (MAP) in Indian cattle. A total of 213 animals were preliminarily selected on the basis of physical body condition score, which was then screened by a panel of diagnostic tests viz. Johnin, ELISA, fecal microscopy, and fecal culture, for the establishment of a case-control resource population. A total of four SNPs viz. rs208222804, rs211654540, rs208814257, and rs210748127 in CD209 gene were genotyped by PCR-RFLP. All SNPs, except rs210748127, were polymorphic in our population. Genotypic-phenotypic associations were assessed by the PROCLOGISTIC procedure of SAS 9.3. The SNP rs208814257 yielded three genotypes viz. CC, CG, and GG, which were significantly (p < 0.05) different in case as compared to the control population. The odds of CC and CG in comparison to GG genotype were 1.21 and 0.40, respectively. The CG genotype was significantly higher in control population, indicating that this genotype may provide resistance against PTB in our resource population. Upon validation in an independent, larger test population and following biological characterization, SNP rs208814257 can be incorporated in marker panel for selection of animals with greater resistance to MAP infection.

Detection of latent forms of Mycobacterium avium subsp. paratuberculosis infection using host biomarker-based ELISAs greatly improves paratuberculosis diagnostic sensitivity

Bovine paratuberculosis (PTB) is a chronic granulomatous enteritis, caused by Mycobacterium avium subsp. paratuberculosis (MAP), responsible for important economic losses in the dairy industry. Current diagnostic methods have low sensitivities for detection of latent forms of MAP infection, defined by focal granulomatous lesions and scarce humoral response or MAP presence. In contrast, patent infections correspond to multifocal and diffuse types of enteritis where there is increased antibody production, and substantial mycobacterial load. Our previous RNA-Seq analysis allowed the selection of five candidate biomarkers overexpressed in peripheral blood of MAP infected Holstein cows with focal (ABCA13 and MMP8) and diffuse (FAM84A, SPARC and DES) lesions vs. control animals with no detectable PTB-associated lesions in intestine and regional lymph nodes. The aim of the current study was to assess the PTB diagnostic potential of commercial ELISAs designed for the specific detection of these biomarkers. The ability of these ELISAs to identify animals with latent and/or patent forms of MAP infection was investigated using serum from naturally infected cattle (n = 88) and non-infected control animals (n = 67). ROC analysis revealed that the ABCA13-based ELISA showed the highest diagnostic accuracy for the detection of infected animals with focal lesions (AUC 0.837, sensitivity 79.25% and specificity 88.06%) and with any type of histological lesion (AUC 0.793, sensitivity 69.41% and specificity 86.57%) improving on the diagnostic performance of the popular IDEXX ELISA and other conventional diagnostic methods. SPARC and MMP8 showed the highest diagnostic accuracy for the detection of animals with multifocal (AUC 0.852) and diffuse lesions (AUC 0.831), respectively. In conclusion, our results suggest that quantification of ABCA13, SPARC and MMP8 by ELISA has the potential for implementation as a diagnostic tool to reliably identify MAP infection, greatly improving early detection of MAP latent infections when antibody responses and fecal shedding are undetectable using conventional diagnostic methods.

Correlations between single nucleotide polymorphisms in bovine CD209, SLC11A1, SP110 and TLR2 genes and estimated breeding values for several traits in Spanish Holstein cattle

Little is known about the correlations between the genetic susceptibility/resistance to Mycobacterium avium subsp. paratuberculosis (Map) infection and the estimated breeding values for type, production and functional traits. Previously, we identified 70 combinations of five single nucleotide polymorphisms (SNPs) in four bovine innate immune genes (SLC11A1, SP110, TLR2, CD209) that are associated with the genetic risk of paratuberculosis (PTB) or Johne's disease progression, which can be graded as low (LOWIN), latent (LATIN), or patent (PATIN) risk. Other possible combinations of these 5 SNPs were grouped in the average group (AVERIN). In the current study, differences in estimated breeding values (EBVs) for several traits were analyzed using linear models in a large cohort of Holstein cows (N = 15656) genotyped across Spain in 2016 or 2017. After the assignment of each genotyped cow to a risk group, cows within the PATIN risk group (N = 1448) had a superior combined genetic index (2797.57), type genetic index (524.62), milk yield (653.92 kg), protein yield (21.77 kg), fat yield (24.82 kg) and economic merit index (125 Euros) compared with the other three risk groups. Statistically significant differences in the longevity scores between the cows that were included in the PATIN risk group (108.85) and the LOWIN (107.82) and AVERIN (107.92) groups were also observed. The associations between the genetic risk groups and PTB diagnostic results were validated in a population of 99 cows from a Spanish farm with a high prevalence of PTB. Significant differences in ELISA readings between the PATIN (65.49 %) and the AVERIN (15.97 %), LATIN (2.11 %), and LOWIN (3.27 %) groups were observed. In addition, significant differences in Map DNA copies/gram of feces were observed between the PATIN and the other three risk groups. These results together with the substantial economic impact of PTB in dairy cattle support the selection of the animals with less susceptibility to PTB in the Spanish breeding program.

Mycobacterium tuberculosis extracellular vesicle-associated lipoprotein LpqH as a potential biomarker to distinguish paratuberculosis infection or vaccination from tuberculosis infection

BACKGROUND

Both bovine tuberculosis (TB) and paratuberculosis (PTB) are serious and widespread bacterial infections affecting many domestic and wild animal species. However, current vaccines do not confer complete protection and cause interference with other diagnostics tests, including bovine TB. Therefore, the development of "Differentiating Infected from Vaccinated Animals" (DIVA) tests are a pressing need. In this study, we have tested the feasibility of mycobacterial extracellular vesicles (EVs) as potential source of biomarkers to discriminate between Mycobacterium bovis infected, Mycobacterium avium subsp. paratuberculosis (MAP) infected and MAP-vaccinated cows. We have, initially, characterized vesicle production in the two most medically relevant species of mycobacteria for livestock, MAP and M. bovis, for being responsible for tuberculosis (TB) and paratuberculosis (PTB).

RESULTS

Our results indicate that these two species produce EVs with different kinetics, morphology and size distribution. Analysis of the immunogenicity of both type of EVs showed some cross reactivity with sera from PTB+ and TB+ cows, suggesting a limited diagnostic capacity for both EVs. Conversely, we noticed that Mycobacterium tuberculosis (Mtb) EVs showed some differential reactivity between sera from MAP-vaccinated or PTB+ cows from TB+ ones. Mass spectrometry analysis (MS) identified a 19-kDa EV-associated lipoprotein as the main source of the differential reactivity.

CONCLUSIONS

LpqH could be a good plasma biomarker with capacity to distinguish PTB+ or MAP-vaccinated cows from cows infected with TB.

Genetic association of polymorphisms in bovine TLR2 and TLR4 genes with Mycobacterium avium subspecies paratuberculosis infection in Indian cattle population

Toll like receptors (TLRs) are pattern recognition molecules involved in cellular recognition of Mycobacterium avium subspecies paratuberculosis (MAP), the infectious agent causing Paratuberculosis (PTB), a notified disease of domestic and wild ruminants. The present study was undertaken to investigate the presence of single nucleotide polymorphisms (SNPs) in TLR2 and TLR4 gene and to evaluate association of these SNPs with occurrence of PTB in Indian cattle. A total of 213 cattle, were subjected to multiple diagnostic tests viz. Johnin PPD, ELISA test (Indigenous and Parachek kit method), fecal microscopy and fecal culture for detection of MAP infection. Based on screening results 51 animals each were assigned to case and control population. Two SNPs viz. rs55617172, rs41830058 in TLR2 gene and two SNPs viz. rs8193046, rs8193060 in TLR4 gene and were genotyped by PCR-RFLP method. All SNPs were found to be polymorphic except rs41830058 in the case-control population. Both SNPs in TLR4 gene but none in TLR2 genes were significantly associated with the occurrence of PTB in our population. The genotypes in SNP rs8193046 and SNP rs8193060 were significantly (P < 0.01) different in case-control population. These findings suggest that SNPs rs8193046 and rs8193060 are likely a potential marker against MAP infection and a selection programme eliminating AG genotype for rs8193046 and CT genotype for rs8193060 might be beneficial in conferring resistance to MAP infection in Indian cattle population.

Deciphering the virulence of Mycobacterium avium subsp. paratuberculosis isolates in animal macrophages using mathematical models

Understanding why pathogenic Mycobacterium avium subsp. paratuberculosis (Map) isolates cause disparate disease outcomes with differing magnitudes of severity is important in designing and implementing new control strategies. We applied a suite of mathematical models: i) general linear, ii) and neurofuzzy logic, to explain how the host of origin of several Map isolates, Map genotype, host, macrophage-based in vitro model and time post-infection contributed to the infection. A logistic growth ordinary differential equation (ODE) model was applied to estimate within macrophage growth rates for the different Map isolates. The models revealed different susceptibilities of bovine and ovine macrophages to Map infection and confirmed distinct virulence profiles for the isolates, judged by their ability to grow within macrophages. Ovine macrophages were able to internalize Map isolates more efficiently than bovine macrophages. While bovine macrophages were able to internalize Map isolates from cattle with more efficiency, ovine macrophages were more efficient in internalizing ovine isolates. Overall, Map isolates from goat and sheep grew minimally within macrophages or did not grow but were able to persist by maintaining its initial population. In contrast, the ability of the bovine isolates and the non-domesticated animal isolates to grow to higher CFU numbers within macrophages suggests that these isolates are more virulent than the sheep and goat isolates, or that these isolates are better adapted to infect domestic ruminants. Overall, our study confirms the different virulence levels for the Map isolates and susceptibility profiles of host macrophages, which is crucial in increasing our understanding of Map infection.

Chronic regional intestinal inflammatory disease: A trans-species slow infection?

Crohn's disease and ulcerative colitis in humans and paratuberculosis in domestic and wild ruminants can be defined as chronic regional intestinal inflammatory diseases (CRIID). This review is a literature overview on these diseases in humans, non-human primates, dogs, cats, rabbits, equids and ruminants with a focus on pathological and microbiological features aimed identifying common characteristics that could lead to a unified pathological classification for a better understanding of their mechanisms and causes. The result is a framework of inflammatory forms throughout the different species indicative of common mechanisms of the slow infection type characterized by a time course varying from weeks to months or even years, and where the inflammatory component would be more prominent in the intestinal interphase between host and environment and be morphologically characterized by an infiltrate ranging from lymphoplasmacytic to histiocytic. This should provide new insights for causation demonstration and therapeutic approaches in human IBD.

Association between combinations of genetic polymorphisms and epidemiopathogenic forms of bovine paratuberculosis

Control of major mycobacterial diseases affecting livestock is a challenging issue that requires different approaches. The use of genetic markers for improving resistance to Mycobacterium avium subsp. paratuberculosis infection in cattle has been explored as a promising population strategy We performed paratuberculosis epidemiopathogenic phenotypic and genotypic characterization involving 24 SNPs in six candidate genes (NOD2, CD209, SLC11A1, SP110, TLR2 and TLR4) on 502 slaughtered Friesian cows. In the current study, we investigate whether recently proposed paratuberculosis (PTB) epidemiopathogenic (EP) forms (apparently free-AF, latent-LAT and patent-PAT) could be associated with some combination of these 24 SNPs. Best EP form grouping was obtained using a combination of 5 SNPs in four genes (CD209: rs210748127; SLC11A1: rs110090506; SP110: rs136859213 and rs110480812; and TLR2: rs41830058). These groups were defined according to the level of infection progression risk to patent epidemiopathogenic forms and showed the following distributions: LOWIN (low) with 39 (8%) cases (94.9% AF/5.1% LAT/0% PAT); LATIN (low) with 17 (3%) cases (5.9% AF/94.1% LAT/0% PAT); AVERIN (average) with 413 (82%) cases (52.1% AF/38.5% LAT/9.4% PAT) and PATIN (patent) with 33 (7%) cases (36.4% AF/24.2% LAT/39.4% PAT). Age of slaughter was significantly higher for LATIN (88.3 months) compared to AVERIN (65.3 months; p = 0.0007) and PATIN (59.1 months; p = 0.0004), and for LOWIN (73.9 months) compared to PATIN (p = 0.0233), and nearly significant compared to AVERIN (p = 0.0572) These results suggest that some selected genetic polymorphisms have a potential use as markers of PTB EP forms and thus add a new tool for the control of this widespread infection.

CARD15 Gene Polymorphisms Are Associated with Tuberculosis Susceptibility in Chinese Holstein Cows

Bovine tuberculosis (BTB) is a significant veterinary and financial problem in many parts of the world. Associations between specific host genes and susceptibility to mycobacterial infections, such as tuberculosis, have been reported in several species. The objective of this study was to identify and evaluate the relationship of single-nucleotide polymorphisms (SNPs) in the CARD15 gene with susceptibility to BTB in Chinese Holstein cows. DNA samples from 201 Chinese Holstein cows (103 cases and 98 controls) were collected from Kunming City, Yuxi City, and Dali City in China. SNPs in the CARD15 gene were assessed using polymerase chain reaction (PCR) and restriction fragment length polymorphism-polymerase chain reaction (RFLP-PCR). Case-control association testing and statistical analysis identified six SNPs associated with susceptibility to BTB in Chinese Holstein cows. The frequency of genotypes C/T, A/G, A/G, A/G, C/T, and A/G in E4 (-37), 208, 1644, 1648, 1799, and E10 (+107), respectively, was significantly higher in cases than in controls, and also the alleles C, A, A, G, T, and A, respectively, were associated with a greater relative risk in cases than in controls. The distribution of two haplotypes, TGGACA and CAGACA, was significantly different between cases and controls. Overall, this case-control study suggested that E4 (-37)(C/T), 208(A/G), 1644(A/G), 1648(A/G), 1799(C/T), and E10 (+107)(A/G) in the CARD15 gene were significantly associated with susceptibility to BTB in Chinese Holstein cows and that haplotypes TGGACA and CAGACA could be used as genetic markers in marker-assisted breeding programs for breeding cows with high resistance to BTB.

Longitudinal data collection of Mycobacterium avium subspecies Paratuberculosis infections in dairy herds: the value of precise field data

Longitudinal infection data on Mycobacterium avium subspecies paratuberculosis (MAP) was collected on three dairy farms in Northeastern United States during approximately 10 years. Precise data on animal characteristics and animal location within farm were collected on these farms. Cows were followed over time with regard to MAP status during biannual fecal and serum sampling and quarterly serum sampling. Approximately 13 000 serum samples, 6500 fecal samples and 2000 tissue samples were collected during these years. Prevalence of positive samples was 1.4% for serological samples, 2.2% in fecal samples and 16.7% in tissue samples. Infection dynamics of MAP was studied and resulted in a number of potential changes in our understanding of MAP infection dynamics. First, a high prevalence of MAP infection was observed in these herds due to lifetime follow up of cows, including slaughter. Second, two distinctly different infection patterns were observed, so called non-progressors and progressors. Non-progressors were characterized by intermittent and low shedding of MAP bacteria and a virtual absence of a humoral immune response. Progressors were characterized by continuous and progressive shedding and a clearly detectable and progressive humoral immune response. Strain typing of MAP isolates on the three farms identified on two of three farms a dominant strain type, indicating that some strains are more successful in terms of transmission and infection progression. Continuous high quality longitudinal data collection turned out to be an essential tool in our understanding of pathobiology and epidemiology of MAP infections in dairy herds.