Retrospective Single Nucleotide Polymorphism Analysis of Host Resistance and Susceptibility to Ovine Johne's Disease Using Restored FFPE DNA.

Johne's disease (JD), also known as paratuberculosis, is a chronic, untreatable gastroenteritis of ruminants caused by Mycobacterium avium subsp. paratuberculosis (MAP) infection. Evidence for host genetic resistance to disease progression exists, although it is limited due to the extended incubation period (years) and diagnostic challenges. To overcome this, previously restored formalin-fixed paraffin embedded tissue (FFPE) DNA from archived FFPE tissue cassettes was utilized for a novel retrospective case-control genome-wide association study (GWAS) on ovine JD. Samples from known MAP-infected flocks with ante- and postmortem diagnostic data were used. Cases (N = 9) had evidence of tissue infection, compared to controls (N = 25) without evidence of tissue infection despite positive antemortem diagnostics. A genome-wide efficient mixed model analysis (GEMMA) to conduct a GWAS using restored FFPE DNA SNP results from the Illumina Ovine SNP50 Bead Chip, identified 10 SNPs reaching genome-wide significance of p < 1 × 10-6 on chromosomes 1, 3, 4, 24, and 26. Pathway analysis using PANTHER and the Kyoto Encyclopedia of Genes and Genomes (KEGG) was completed on 45 genes found within 1 Mb of significant SNPs. Our work provides a framework for the novel use of archived FFPE tissues for animal genetic studies in complex diseases and further evidence for a genetic association in JD.

The role of interleukin-10 receptor alpha (IL10Rα) in Mycobacterium avium subsp. paratuberculosis infection of a mammary epithelial cell line.

BACKGROUND: Johne's disease is a chronic wasting disease caused by the bacterium Mycobacterium avium subspecies paratuberculosis (MAP). Johne's disease is highly contagious and MAP infection in dairy cattle can eventually lead to death. With no available treatment for Johne's disease, genetic selection and improvements in management practices could help reduce its prevalence. In a previous study, the gene coding interleukin-10 receptor subunit alpha (IL10Rα) was associated with Johne's disease in dairy cattle. Our objective was to determine how IL10Rα affects the pathogenesis of MAP by examining the effect of a live MAP challenge on a mammary epithelial cell line (MAC-T) that had IL10Rα knocked out using CRISPR/cas9. The wild type and the IL10Rα knockout MAC-T cell lines were exposed to live MAP bacteria for 72 h. Thereafter, mRNA was extracted from infected and uninfected cells. Differentially expressed genes were compared between the wild type and the IL10Rα knockout cell lines. Gene ontology was performed based on the differentially expressed genes to determine which biological pathways were involved. RESULTS: Immune system processes pathways were targeted to determine the effect of IL10Rα on the response to MAP infection. There was a difference in immune response between the wild type and IL10Rα knockout MAC-T cell lines, and less difference in immune response between infected and not infected IL10Rα knockout MAC-T cells, indicating IL10Rα plays an important role in the progression of MAP infection. Additionally, these comparisons allowed us to identify other genes involved in inflammation-mediated chemokine and cytokine signalling, interleukin signalling and toll-like receptor pathways. CONCLUSIONS: Identifying differentially expressed genes in wild type and ILR10α knockout MAC-T cells infected with live MAP bacteria provided further evidence that IL10Rα contributes to mounting an immune response to MAP infection and allowed us to identify additional potential candidate genes involved in this process. We found there was a complex immune response during MAP infection that is controlled by many genes.

The Fur-like regulatory protein MAP3773c modulates key metabolic pathways in Mycobacterium avium subsp. paratuberculosis under in-vitro iron starvation.

Johne's disease (JD) is a chronic enteric infection of dairy cattle worldwide. Mycobacterium avium subsp. paratuberculosis (MAP), the causative agent of JD, is fastidious often requiring eight to sixteen weeks to produce colonies in culture-a major hurdle in the diagnosis and therefore in implementation of optimal JD control measures. A significant gap in knowledge is the comprehensive understanding of the metabolic networks deployed by MAP to regulate iron both in-vitro and in-vivo. The genome of MAP carries MAP3773c, a putative metal regulator, which is absent in all other mycobacteria. The role of MAP3773c in intracellular iron regulation is poorly understood. In the current study, a field isolate (K-10) and an in-frame MAP3773c deletion mutant (ΔMAP3773c) derived from K-10, were exposed to iron starvation for 5, 30, 60, and 90 min and RNA-Seq was performed. A comparison of transcriptional profiles between K-10 and ΔMAP3773c showed 425 differentially expressed genes (DEGs) at 30 min time post-iron restriction. Functional analysis of DEGs in ΔMAP3773c revealed that pantothenate (Pan) biosynthesis, polysaccharide biosynthesis and sugar metabolism genes were downregulated at 30 min post-iron starvation whereas ATP-binding cassette (ABC) type metal transporters, putative siderophore biosynthesis, PPE and PE family genes were upregulated. Pathway analysis revealed that the MAP3773c knockout has an impairment in Pan and Coenzyme A (CoA) biosynthesis pathways suggesting that the absence of those pathways likely affect overall metabolic processes and cellular functions, which have consequences on MAP survival and pathogenesis.

Genetic parameters and trends for Johne's disease in US Holsteins: An updated study.

Johne's disease (JD) is an infectious enteric disease in ruminants, causing substantial economic loss annually worldwide. This work aimed to estimate JD's genetic parameters and the phenotypic and genetic trends by incorporating recent data. It also explores the feasibility of a national genetic evaluation for JD susceptibility in Holstein cattle in the United States. The data were extracted from a JD data repository, maintained at the Council on Dairy Cattle Breeding, and initially supplied by 2 dairy record processing centers. The data comprised 365,980 Holstein cows from 1,048 herds participating in a voluntary control program for JD. Two protocol kits, IDEXX Paratuberculosis Screening Ab Test (IDX) and Parachek 2 (PCK), were used to analyze milk samples with the ELISA technique. Test results from the first 5 parities were considered. An animal was considered infected if it had at least one positive outcome. The overall average of JD incidence was 4.72% in these US Holstein cattle. Genotypes of 78,964 SNP markers were used for 25,000 animals randomly selected from the phenotyped population. Variance components and genetic parameters were estimated based on 3 models, namely, a pedigree-only threshold model (THR), a single-step threshold model (ssTHR), and a single-step linear model (ssLR). The posterior heritability estimates of JD susceptibility were low to moderate: 0.11 to 0.16 based on the 2 threshold models and 0.05 to 0.09 based on the linear model. The average reliability of EBVs of JD susceptibility using single-step analysis for animals with or without phenotypes varied from 0.18 (THR) to 0.22 (ssLR) for IDX and from 0.14 (THR) to 0.18 (ssTHR and ssLR) for PCK. Despite no prior direct genetic selection against JD, the estimated genetic trends of JD susceptibility were negative and highly significant. The correlations of bulls' PTA with economically important traits such as milk yield, milk protein, milk fat, somatic cell score, and mastitis were low, indicating a nonoverlapping genetic selection process with traits in current genetic evaluations. Our results suggest the feasibility of reducing the JD incidence rate by incorporating it into the national genetic evaluation programs.

Alternative splicing of pre-mRNA modulates the immune response in Holstein cattle naturally infected with Mycobacterium avium subsp. paratuberculosis.

Little is known about the role of alternative splicing (AS) in regulating gene expression in Mycobacteria-infected individuals in distinct stages of infection. Pre-mRNA AS consists of the removal of introns and the assembly of exons contained in eukaryotic genes. AS events can influence transcript stability or structure with important physiological consequences. Using RNA-Seq data from peripheral blood (PB) and ileocecal valve (ICV) samples collected from Holstein cattle with focal and diffuse paratuberculosis (PTB)-associated histopathological lesions in gut tissues and without lesions (controls), we detected differential AS profiles between the infected and control groups. Four of the identified AS events were experimentally validated by reverse transcription-digital droplet PCR (RT-ddPCR). AS events in several genes correlated with changes in gene expression. In the ICV of animals with diffuse lesions, for instance, alternatively spliced genes correlated with changes in the expression of genes involved in endocytosis, antigen processing and presentation, complement activation, and several inflammatory and autoimmune diseases in humans. Taken together, our results identified common mechanisms of AS involvement in the pathogenesis of PTB and human diseases and shed light on novel diagnostic and therapeutic interventions to control these diseases.

Advances in understanding the genetic architecture of antibody response to paratuberculosis in sheep by heritability estimate and LDLA mapping analyses and investigation of candidate regions using sequence-based data.

BACKGROUND: Paratuberculosis is a contagious and incurable disease that is caused by Mycobacterium avium subsp. paratuberculosis (MAP) with significant negative effects on animal welfare and farm profitability. Based on a large naturally infected flock over 12 years, we analyzed repeated enzyme-linked immunosorbent assay tests (ELISA), OvineSNP50 BeadChip genotypes and whole-genome sequences imputed from 56 influential animals. The main goals were to estimate the genetic parameters of proxy traits for resistance to MAP, identify genomic regions associated with the host's immune response against MAP and search for candidate genes and causative mutations through association and functional annotation analyses of polymorphisms identified by sequencing. RESULTS: Two variables were derived from ELISA tests. The first, a binary variable, assessed the infection status of each animal over the entire productive life, while the second considered the level of antibody recorded over time. Very similar results were obtained for both variables. Heritability estimates of about 0.20 were found and a significant region capturing 18% and 13% of the genetic variance was detected on ovine chromosome 20 by linkage disequilibrium and linkage analysis on OvineSNP50 positions. Functional annotation and association analyses on the imputed sequence polymorphisms that were identified in this region were carried out. No significant variants showed a functional effect on the genes that mapped to this region, most of which belong to the major histocompatibility complex class II (MHC II). However, the conditional analysis led to the identification of two significant polymorphisms that can explain the genetic variance associated with the investigated genomic region. CONCLUSIONS: Our results confirm the involvement of the host's genetics in susceptibility to MAP in sheep and suggest that selective breeding may be an option to limit the infection. The estimated heritability is moderate with a relevant portion being due to a highly significant region on ovine chromosome 20. The results of the combined use of sequence-based data and functional analyses suggest several genes belonging to the MHC II as the most likely candidates, although no mutations in their coding regions showed a significant association. Nevertheless, information from genotypes of two highly significant polymorphisms in the region can enhance the efficiency of selective breeding programs.

The pan-genome of Mycobacterium avium subsp. paratuberculosis (Map) confirms ancestral lineage and reveals gene rearrangements within Map Type S.

BACKGROUND: To date genomic studies on Map have concentrated on Type C strains with only a few Type S strains included for comparison. In this study the entire pan-genome of 261 Map genomes (205 Type C, 52 Type S and 4 Type B) and 7 Mycobacterium avium complex (Mac) genomes were analysed to identify genomic similarities and differences between the strains and provide more insight into the evolutionary relationship within this Mycobacterial species. RESULTS: Our analysis of the core genome of all the Map isolates identified two distinct lineages, Type S and Type C Map that is consistent with previous phylogenetic studies of Map. Pan-genome analysis revealed that Map has a larger accessory genome than Mycobacterium avium subsp. avium (Maa) and Type C Map has a larger accessory genome than Type S Map. In addition, we found large rearrangements within Type S strains of Map and little to none in Type C and Type B strains. There were 50 core genes identified that were unique to Type S Map and there were no unique core genes identified between Type B and Type C Map strains. In Type C Map we identified an additional CE10 CAZyme class which was identified as an alpha/beta hydrolase and an additional polyketide and non-ribosomal peptide synthetase cluster. Consistent with previous analysis no plasmids and only incomplete prophages were identified in the genomes of Map. There were 45 hypothetical CRISPR elements identified with no associated cas genes. CONCLUSION: This is the most comprehensive comparison of the genomic content of Map isolates to date and included the closing of eight Map genomes. The analysis revealed that there is greater variation in gene synteny within Type S strains when compared to Type C indicating that the Type C Map strain emerged after Type S. Further analysis of Type C and Type B genomes revealed that they are structurally similar with little to no genetic variation and that Type B Map may be a distinct clade within Type C Map and not a different strain type of Map. The evolutionary lineage of Maa and Map was confirmed as emerging after M. hominissuis.

Association between major histocompatibility complex haplotypes and susceptibility of unvaccinated and vaccinated cattle to paratuberculosis.

Bovine Johne's disease (BJD) or paratuberculosis is caused by Mycobacterium avium spp. paratuberculosis (MAP) and is a worldwide problem among domestic and wild ruminants. While vaccines are available, natural differences in background immunity between breeds within species and between individuals within herds suggest that genetic differences may be able to be exploited in marker-assisted selection as an aid to disease control. The major histocompatibility complex (MHC) is an important component in immune recognition with considerable genetic variability. In this study, associations between the MHC and resistance to BJD were explored in dairy cattle across two herds in which some of the cattle had been vaccinated with Silirum® (n = 540 cows). A BJD susceptible animal was exposed to MAP and became infected, while a resistant animal was exposed but did not become infected. There are different ways to define both exposure and infection, with different levels of stringency, therefore many classifications of the same set of animals are possible and were included in the analysis. The polymorphic regions of major histocompatibility complex class I (MHC I) and class II (MHC II) genes were amplified from the genomic DNA by PCR and sequenced, targeting exons 2 and 3 of the classical and non-classical MHC I genes and exon 2 from the DRB3, DQA1, DQA2 + 3 and DQB MHC II genes. The frequencies of MHC I and MHC II haplotypes and alleles were determined in susceptible and resistant populations. In unvaccinated animals, seven MHC I haplotypes and seven MHC II haplotypes were associated with susceptibility while two MHC I and six MHC II haplotypes were associated with resistance (P < 0.05). In vaccinated animals, two MHC I and three MHC II haplotypes were associated with susceptibility, while one MHC I and two MHC II haplotypes were associated with resistance (P < 0.05). The alleles in significant haplotypes were also identified. Case definitions with higher stringency resulted in fewer animals being included in the analyses, but the power to detect an association was not reduced and there was an increase in strength and consistency of associations. Consistent use of stringent case definitions is likely to improve agreement in future association studies.

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.

Long Non-Coding RNAs and Their "Discrete" Contribution to IBD and Johne's Disease-What Stands out in the Current Picture? A Comprehensive Review.

Non-coding RNAs (ncRNA) have paved the way to new perspectives on the regulation of gene expression, not only in biology and medicine, but also in associated fields and technologies, ensuring advances in diagnostic means and therapeutic modalities. Critical in this multistep approach are the associations of long non-coding RNA (lncRNA) with diseases and their causal genes in their networks of interactions, gene enrichment and expression analysis, associated pathways, the monitoring of the involved genes and their functional roles during disease progression from one stage to another. Studies have shown that Johne's Disease (JD), caused by Mycobacterium avium subspecies partuberculosis (MAP), shares common lncRNAs, clinical findings, and other molecular entities with Crohn's Disease (CD). This has been a subject of vigorous investigation owing to the zoonotic nature of this condition, although results are still inconclusive. In this review, on one hand, the current knowledge of lncRNAs in cells is presented, focusing on the pathogenesis of gastrointestinal-related pathologies and MAP-related infections and, on the other hand, we attempt to dissect the associated genes and pathways involved. Furthermore, the recently characterized and novel lncRNAs share common pathologies with IBD and JD, including the expression, molecular networks, and dataset analysis results. These are also presented in an attempt to identify potential biomarkers pertinent to cattle and human disease phenotypes.

Impact of Mycobacterium avium subsp. paratuberculosis infection on bovine IL10RA knockout mammary epithelial (MAC-T) cells.

Mycobacterium avium subsp. Paratuberculosis (MAP) is an intracellular pathogen that causes Johne's disease (JD) in cattle and other ruminants. IL10RA encodes the alpha chain of the IL-10 receptor that binds the cytokine IL-10, and is one of the candidate genes that have been found to be associated with JD infection status. In this study, a previously developed IL10RA knockout (IL10RAKO) bovine mammary epithelial (MAC-T) cell line and wild-type (WT) MAC-T cells were infected with live MAP for 72 h to identify potential immunoregulatory miRNAs, inflammatory genes, and cytokines/chemokines impacted by MAP infection in the presence/absence of IL10RA. Cytokine and chemokine concentrations in culture supernatants were measured by multiplexing immunoassay. Total RNA was extracted from the MAC-T cells, and qPCR was performed to determine the expression of inflammatory genes and selected bovine miRNAs. Results showed that the levels of TNF-α, IL-6, CXCL8, CXCL10, CCL2, and CCL3 were significantly induced in WT MAC-T cells and IL-10 was significantly inhibited post-MAP infection. However, IL10RAKO MAC-T cells had greater secretion of TNF-α, IL-6, IFN-γ, CCL3, CCL4, CXCL8, and CXCL10, and lower secretion of VEGF-α. Moreover, the expression of inflammatory genes (TNF-α, IL-1α, IL-6) was also more significantly induced in IL10RAKO cells than in WT MAC-T cells post-MAP-infection, and unlike the WT cells, anti-inflammatory cytokines IL-10 and SOCS3 and chemokines CCL2 were not significantly induced. In addition, the expression of miRNAs (miR133b, miR-92a, and miR-184) was increased in WT MAC-T cells post-MAP-infection; however, there was no significant induction of these miRNAs in the IL10RAKO cells, which suggests IL10 receptor is somehow involved in regulating the miRNA response to MAP infection. Target gene function analysis further suggests that miR-92a may be involved in interleukin signaling, and miR-133b and miR-184 may be involved in other signaling pathways. These findings support the involvement of IL10RA in the regulation of innate immune response to MAP.

Improved DNA Amplification of the Hallmark IS900 Element in Mycobacterium avium subsp. paratuberculosis: a Reexamination Based on Whole-Genome Sequence Analysis.

Amplification of the IS900 multicopy element is a hallmark nucleic acid-based diagnostic test for Mycobacterium avium subsp. paratuberculosis, which causes Johne's disease in ruminants. This assay is frequently used to determine the presence of the bacterium in feces of infected cattle and sheep. Two IS900 primer sets developed in the 1990s were widely used for decades, and their use has continued in current studies. However, these primers were developed prior to the availability of complete genome sequences. Recent sequence analysis of the binding locations for one primer pair (P90/P91) identified errors and binding inefficiencies that can be easily corrected to further increase detection sensitivity. The P90 primer is missing two nucleotides that should be present near the 3' end, and it does not bind all copies of IS900 due to 5' deletions at some IS900 loci. These IS900 primer pairs, along with newly developed primers, were tested by real-time PCR on purified genomic DNA to determine which primer set performed the best and how primer design errors affect amplification efficiencies. The newly designed PCR primer set (JB5) showed increased sensitivity by two to three quantification cycles using purified genomic DNA and was similar in efficiency to 150C/921. These tests were extended using DNA from feces and tissues of infected cows, which showed similar results. Finally, a 167-bp partial duplication of IS900 was found in type I strains. Although P90 and P91 primers successfully amplify M. avium subsp. paratuberculosis DNA, their use should be discontinued in favor of more efficient primer pairs in future studies. IMPORTANCE This study is an example of how applied genomic analysis can aid diagnostic test improvements. Detection of Mycobacterium avium subsp. paratuberculosis infection of livestock prior to the appearance of clinical disease signs is very difficult but essential for identifying animals shedding the bacterium to prevent transmission of Johne's disease. Total M. avium subsp. paratuberculosis quantity in the feces as determined by real-time PCR (qPCR) using the IS900 target indicates bacterial shedding status and potential for transmission of the pathogen. However, legacy primers designed prior to the availability of complete genome sequences that are used in these tests to detect M. avium subsp. paratuberculosis were based on data from only a single copy of IS900 and not considering all copies collectively as a group. This approach resulted in primer design errors which can be easily corrected to improve test sensitivities. We tested original primers that contain these errors and their corrected versions by qPCR and showed improved sensitivity on purified genomic DNA as well as fecal and tissue samples. These findings may help detect the organism from environmental samples on farms where sensitivity is currently lacking.

Genome-wide association study for antibody response to Mycobacterium avium ssp. paratubeculosis in goats.

Mycobacterium avium ssp. paratuberculosis (MAP), causes Johne's disease (JD), or paratuberculosis, a chronic enteritis of ruminants, which in goats is characterized by ileal lesions. The work described here is a case-control association study using the Illumina Caprine SNP50 BeadChip to unravel the genes involved in susceptibility of goats to JD. Goats in herds with a high occurrence of Johne's disease were classified as healthy or infected based on the level of serum antibodies against MAP, and 331 animals were selected for the association study. Goats belonged to the Jonica (157) and Siriana breeds (174). Whole-genome association analysis identified one region suggestive of significance associated with an antibody response to MAP on chromosome 7 (p-value = 1.23 × 10-5 ). These results provide evidence for genetic loci involved in the antibody response to MAP in goats.

Analysis of long non-coding RNA expression profile of bovine monocyte-macrophage infected by Mycobacterium avium subsp. paratuberculosis.

Mycobacterium avium subsp. paratuberculosis (MAP) is the causative agent of paratuberculosis. As a potential zoonotic pathogen, MAP also seriously threatens human health and social security. At present, long non-coding RNA (lncRNA) has attracted wide attention as an useful biomarker in various diseases. Therefore, our study analyzed the lncRNA expression profiles and lncRNA-mRNA regulatory network of MAP infected bovine monocytes-macrophages and uninfected bovine cells by high-throughput sequencing. A total of 4641 differentially expressed lncRNAs genes were identified, including 3111 up-regulated genes and 1530 down-regulated genes. In addition, lncRNA-mRNA interaction analysis was performed to predict the target genes of lncRNA. Among them, after MAP infection, 86 lncRNAs targeted to mRNA, of which only 6 genes were significantly different. The results of Gene Ontology (GO) enrichment analysis showed that the differentially expressed genes significantly enriched in functional groups were related to immune regulation. Multiple signal pathways including NF-κB, NOD-like receptor, Cytokine-cytokine receptor, Toll-like receptor signaling pathway, Chemokine signaling pathway, and other important biochemical, metabolic and signal transduction pathways were enriched in Kyoto Encyclopedia of Genes and Genomes (KEGG). In this study, analysis of macrophage transcriptomes in response to MAP infection is expected to provide key information to deeply understand role of the pathogen in initiating an inappropriate and persistent infection in susceptible hosts and molecular mechanisms that might underlie the early phases of paratuberculosis.

New insights into the genetic resistance to paratuberculosis in Holstein cattle via single-step genomic evaluation.

BACKGROUND: Bovine paratuberculosis, or Johne's disease (JD), is a contagious and incurable disease caused by Mycobacterium avium subsp. paratuberculosis (MAP). It has adverse effects on animal welfare and is very difficult to control, leading to serious economic consequences. An important line of defense to this disease is host genetic resistance to MAP, which, when it will be more fully understood, could be improved through selective breeding. Using a large dataset of Holstein cows (161,253 animals including 56,766 cows with ELISA serological phenotypes and 12,431 animals with genotypes), we applied a single-step single nucleotide polymorphism (SNP) best linear unbiased prediction approach to investigate the genetic determinism underlying resistance to this disease (heritability estimate and identification of relevant genomic regions) and estimated genetic trends, reliability, and relative risk factors associated with genomic predictions. RESULTS: Resistance to JD was moderately heritable (0.14) and 16 genomic regions were detected that accounted for at least 0.05% of the breeding values variance (GV) in resistance to JD, and were located on chromosomes 1, 3, 5, 6, 7, 19, 20, 21, 23, 25, and 27, with the highest percentage of variance explained by regions on chromosomes 23 (0.36% GV), 5 (0.22% GV), 1 (0.14% GV), and 3 (0.13% GV). When estimated for the whole chromosomes, the autosomes with the largest overall contributions were chromosomes 3 (5.3% GV), 10 (4.8%), 23 (4.7%), 1 (3.6%), 7 (3.4%), 5 (2.9%), 12 (2.5%), 11 (2.2%), and 13 (2%). We estimated a slightly favorable genetic trend in resistance to JD over the last two decades, which can be explained by a low positive genetic correlation between resistance to JD and total merit index (+ 0.06). Finally, in a validation population of 907 cows, relatively reliable genomic predictions (reliability = 0.55) were obtained, which allowed the identification of cows at high risk of infection. CONCLUSIONS: This study provides new insights into the genetic determinism of resistance to JD and shows that this trait can be predicted from SNP genotypes. It has led to the implementation of a single-step genomic evaluation that should rapidly become an effective tool for controlling paratuberculosis on French Holstein farms.

Identification of Alternative Splicing Events Associated with Paratuberculosis in Dairy Cattle Using Multi-Tissue RNA Sequencing Data.

Paratuberculosis is a major endemic disease caused by Mycobacterium avium subspecies paratuberculosis (MAP) infection and leads to huge economic loss in the dairy sector worldwide. Alternative splicing (AS) events, playing indispensable regulatory roles in many protein functions and biological pathways, are shown to be associated with complex traits and diseases. In this study, by integrating the RNA sequencing (RNA-seq) data of 24 samples from three tissues (peripheral blood, jejunum and salivary gland) of Holstein cows, we obtained 2,706,541,696 uniquely mapped reads in total that represented 12,870 expressed genes, and detected 4285 differentially expressed genes (DEGs) between MAP-infected and healthy cows (p < 0.05). Of them, 92 differentially expressed splicing factors (DESFs) were included. Further, 119, 150 and 68 differential alternative splicing (DAS) events between MAP-infected and healthy cows were identified in peripheral blood, jejunum and salivary glands, respectively. Of note, six DAS events were highly and significantly correlated with the DESFs (R2 > 0.9; p < 0.01), and their corresponding genes (COPI coat complex subunit gamma 2gene (COPG2), kinesin family member 2C gene (KIF2C), exocyst complex component 7 (EXOC7), Rab9 effector protein with kelch motifs gene (RABEPK), deoxyribonuclease 1 gene (DNASE1) and early endosome antigen 1gene (EEA1)) were significantly enriched in immune response such as vesicle-mediated transport, regulation of acute inflammatory response and tuberculosis through gene ontology (GO) and KEGG analysis. KS test showed that the DAS events in the EXOC7 and KIF2C genes indeed displayed differences between MAP-infected cows and healthy cows. The DAS in EXOC7 might produce a new protein sequence with lack of 23 amino acids, and the DAS in KIF2C induced a stop codon of premature occurrence and resulted in a lack of functional domain. In summary, this study identified the DAS events and corresponding genes related to MAP-infection base on the RNA-seq data from multiple tissues of Holstein cows, providing novel insights into the regulatory mechanisms underpinning paratuberculosis in dairy cattle.

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.

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-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.

Validation of suitable reference genes for quantitative expression analysis by qPCR in bovine terminal ileum and ileocecal valve.

The use of reference genes is commonly accepted as the most reliable approach to normalize qRT-PCR and to reduce possible errors in the quantification of gene expression. The aim of this study was to identify a set of reference genes suitable for gene expression analysis in the distal portion of small intestine and ileocecal valve in cattle. These sites of intestine are of interest in veterinary science as they are the main sites of inflammation caused by Mycobacterium avium subsp. paratuberculosis, agent of paratuberculosis. We employed ten PCR assays for commonly used reference genes belonging to various functional classes and then determined their expression stability. The most stable genes were RPL13A and HMBS, followed by TFRC and B-ACT. NormFinder analysis provided similar results with B-ACT as the best reference gene, followed by RLP13A and TFRC. This validated gene panel may be useful for studies on paratuberculosis aiming to identify genes differentially expressed by qRT-PCR.