Dataset Information

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

ABSTRACT: 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^-7) 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.

SUBMITTER: Canive M

PROVIDER: S-EPMC8505495 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:The sensitivity of current diagnostics for Johne's disease, a slow, progressing enteritis in ruminants caused by Mycobacterium avium subsp. paratuberculosis (MAP), is too low to reliably detect all infected animals in the subclinical stage. The objective was to identify individual metabolites or metabolite profiles that could be used as biomarkers of early MAP infection in ruminants. In a monthly follow-up for 17 months, calves infected at 2 weeks of age were compared with aged-matched controls. Sera from all animals were analyzed by 1H nuclear magnetic resonance spectrometry. Spectra were acquired, processed, and quantified for analysis. The concentration of many metabolites changed over time in all calves, but some metabolites only changed over time in either infected or non-infected groups and the change in others was impacted by the infection. Hierarchical multivariate statistical analysis achieved best separation between groups between 300 and 400 days after infection. Therefore, a cross-sectional comparison between 1-year-old calves experimentally infected at various ages with either a high- or a low-dose and age-matched non-infected controls was performed. Orthogonal Projection to Latent Structures Discriminant Analysis (OPLS DA) yielded distinct separation of non-infected from infected cattle, regardless of dose and time (3, 6, 9 or 12 months) after infection. Receiver Operating Curves demonstrated that constructed models were high quality. Increased isobutyrate in the infected cattle was the most important agreement between the longitudinal and cross-sectional analysis. In general, high- and low-dose cattle responded similarly to infection. Differences in acetone, citrate, glycerol and iso-butyrate concentrations indicated energy shortages and increased fat metabolism in infected cattle, whereas changes in urea and several amino acids (AA), including the branched chain AA, indicated increased protein turnover. In conclusion, metabolomics was a sensitive method for detecting MAP infection much sooner than with current diagnostic methods, with individual metabolites significantly distinguishing infected from non-infected cattle.

Project description: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 (h2) 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 h2 estimates > 0.080. After performing the GWAS, the test combinations that showed SNPs with a significant association (PFDR ≤ 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.

Project description:Purpose: Since the peripheral blood (PB) is the most available physiological fluid for the detection of biomarkers, in the current we examine whether it recapitulates, at least in part, the transcriptome of the ileocecal valve (ICV), the primary site of MAP colonization. For this purpose, RNA-Seq was used to identify host genes differentially expressed (DE) in ICV and PB samples collected from PTB-infected animals with focal or diffuse lesions in gut tissues versus control animals. Methods: RNA-Seq libraries were single-end sequenced in a 1x75 format using an Illumina NextSeq 500 sequencer. The raw reads were filtered by their length (minimum size 60 bp long) and percentage of ambiguous base N less than 5 % using Prinseq-lite. Trimmed reads were subsequently mapped to the Bos Taurus reference genome (Bos_taurus.UMD3.1. version 87) with TopHat mapper. The resulting alignment files were provided to Cufflinks to generate a transcriptome assembly for each condition. These assemblies were then merged together using Cuffmerge, which is included in the Cufflinks package. This merged assembly provided a uniform basis for calculating gene and transcript expression levels in each condition. The reads and merged assembly were fed to Cuffdiff which calculates expression levels and tested the statistical significance of each observed change in expression. The fold change (in log2 scale), P-values and false discovery rates (FDR) for each gene were obtained. The genes with a FDR-adjusted threshold ˂ 0.05 were considered differentially expressed (DE) when compared to the control group. Results: Our results demonstrated both shared, and PB and ICV-specific gene expression in response to a natural Map infection. As expected, the number of differentially expressed (DE) genes was larger in the ICV than in the PB samples. Among the DE genes in the PB and ICV samples, there were some common genes irrespective of the type of lesion including the C-X-C motif chemokine ligand 8 (CXCL8/IL8), apolipoprotein L (APOLD1), and the interferon inducible protein 27 (IFI27). The biological processes (BP) enriched in the ICV gene expression profiles from the cows with focal lesions included the killing of cells of other organism, defense response, immune response and the regulation of neutrophil chemotaxis. Two of these BP, the defense and immune response, were also enriched in the ICV and PB from the cows with diffuse lesions. Metabolic analysis of the DE genes revealed that the N-glycan biosynthesis, bile secretion, one-carbon pool by folate and purine metabolism were significantly enriched in the ICV from the cows with focal lesions. In the ICV from cows with diffuse lesions; the valine, leucine and isoleucine degradation route, purine metabolism, vitamin digestion and absorption and the cholesterol routes were enriched Conclusions: This study is the first to simultaneously describe transcriptomic changes in PB and ICV samples of cattle naturally infected with MAP. Several important observations were made. First, by comparing the transcriptomic profiles of two MAP-targeted tissues we described both unique and overlapping changes in the transcriptome of the infected cows versus the control group. Second, our study highlighted the ability of the RNA-Seq technology to reveal roles for genes that have not been previously implicated in the host response to MAP infection. Third, our transcriptomic analysis provided potential biomarkers for the development of future diagnostic tools, vaccines and therapeutics.

Dataset Information

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

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets