1
|
Kravitz A, Liao M, Morota G, Tyler R, Cockrum R, Manohar BM, Ronald BSM, Collins MT, Sriranganathan N. Retrospective Single Nucleotide Polymorphism Analysis of Host Resistance and Susceptibility to Ovine Johne's Disease Using Restored FFPE DNA. Int J Mol Sci 2024; 25:7748. [PMID: 39062990 PMCID: PMC11276633 DOI: 10.3390/ijms25147748] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/12/2024] [Revised: 06/26/2024] [Accepted: 06/28/2024] [Indexed: 07/28/2024] Open
Abstract
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.
Collapse
Affiliation(s)
- Amanda Kravitz
- Center for One Health Research, Department of Biomedical Sciences and Pathobiology, Virginia-Maryland College of Veterinary Medicine, Virginia Polytechnic Institute and State University, Blacksburg, VA 24061, USA
| | - Mingsi Liao
- Department of Animal and Poultry Sciences, Virginia Polytechnic Institute and State University, Blacksburg, VA 24061, USA
| | - Gota Morota
- Department of Animal and Poultry Sciences, Virginia Polytechnic Institute and State University, Blacksburg, VA 24061, USA
| | - Ron Tyler
- Center for One Health Research, Department of Biomedical Sciences and Pathobiology, Virginia-Maryland College of Veterinary Medicine, Virginia Polytechnic Institute and State University, Blacksburg, VA 24061, USA
| | - Rebecca Cockrum
- Department of Animal and Poultry Sciences, Virginia Polytechnic Institute and State University, Blacksburg, VA 24061, USA
| | - B. Murali Manohar
- Department of Veterinary Pathology, Tamilnadu Veterinary and Animal Sciences University, Madhavaram Milk Colony, Chennai 600051, Tamil Nadu India, India
| | - B. Samuel Masilamoni Ronald
- Department of Veterinary Pathology, Tamilnadu Veterinary and Animal Sciences University, Madhavaram Milk Colony, Chennai 600051, Tamil Nadu India, India
| | - Michael T. Collins
- Department of Pathobiological Sciences, School of Veterinary Medicine, University of Wisconsin-Madison, Madison, WI 53706, USA
| | - Nammalwar Sriranganathan
- Center for One Health Research, Department of Biomedical Sciences and Pathobiology, Virginia-Maryland College of Veterinary Medicine, Virginia Polytechnic Institute and State University, Blacksburg, VA 24061, USA
| |
Collapse
|
2
|
Somenzi E, Partel E, Barbato M, Chero Osorio AM, Colli L, Franceschi N, Mantovani R, Pilla F, Komjanc M, Achilli A, Hauffe HC, Ajmone Marsan P. Genetic legacy and adaptive signatures: investigating the history, diversity, and selection signatures in Rendena cattle resilient to eighteenth century rinderpest epidemics. Genet Sel Evol 2024; 56:32. [PMID: 38698323 PMCID: PMC11064358 DOI: 10.1186/s12711-024-00900-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/20/2023] [Accepted: 04/05/2024] [Indexed: 05/05/2024] Open
Abstract
BACKGROUND Rendena is a dual-purpose cattle breed, which is primarily found in the Italian Alps and the eastern areas of the Po valley, and recognized for its longevity, fertility, disease resistance and adaptability to steep Alpine pastures. It is categorized as 'vulnerable to extinction' with only 6057 registered animals in 2022, yet no comprehensive analyses of its molecular diversity have been performed to date. The aim of this study was to analyse the origin, genetic diversity, and genomic signatures of selection in Rendena cattle using data from samples collected in 2000 and 2018, and shed light on the breed's evolution and conservation needs. RESULTS Genetic analysis revealed that the Rendena breed shares genetic components with various Alpine and Po valley breeds, with a marked genetic proximity to the Original Braunvieh breed, reflecting historical restocking efforts across the region. The breed shows signatures of selection related to both milk and meat production, environmental adaptation and immune response, the latter being possibly the result of multiple rinderpest epidemics that swept across the Alps in the eighteenth century. An analysis of the Rendena cattle population spanning 18 years showed an increase in the mean level of inbreeding over time, which is confirmed by the mean number of runs of homozygosity per individual, which was larger in the 2018 sample. CONCLUSIONS The Rendena breed, while sharing a common origin with Brown Swiss, has developed distinct traits that enable it to thrive in the Alpine environment and make it highly valued by local farmers. Preserving these adaptive features is essential, not only for maintaining genetic diversity and enhancing the ability of this traditional animal husbandry to adapt to changing environments, but also for guaranteeing the resilience and sustainability of both this livestock system and the livelihoods within the Rendena valley.
Collapse
Affiliation(s)
- Elisa Somenzi
- DIANA Dipartimento di Scienze Animali, della Nutrizione e degli Alimenti, Università Cattolica del Sacro Cuore, Piacenza, Italy.
| | - Erika Partel
- Unità risorse foraggere e produzioni zootecniche, Centro Trasferimento Tecnologico, Fondazione Edmund Mach, S. Michele all'Adige, Trento, Italy
| | - Mario Barbato
- DIANA Dipartimento di Scienze Animali, della Nutrizione e degli Alimenti, Università Cattolica del Sacro Cuore, Piacenza, Italy
- Department of Veterinary Science, Università degli Studi di Messina, Messina, Italy
| | - Ana María Chero Osorio
- Dipartimento di Biologia e Biotecnologie "L. Spallanzani", University of Pavia, Pavia, Italy
| | - Licia Colli
- DIANA Dipartimento di Scienze Animali, della Nutrizione e degli Alimenti, Università Cattolica del Sacro Cuore, Piacenza, Italy
- Centro di Ricerca Sulla Biodiversità e sul DNA Antico, BioDNA, Università Cattolica del Sacro Cuore, Piacenza, Italy
| | - Niccolò Franceschi
- DIANA Dipartimento di Scienze Animali, della Nutrizione e degli Alimenti, Università Cattolica del Sacro Cuore, Piacenza, Italy
| | - Roberto Mantovani
- Department of Agronomy, Food, Natural Resources, Animals, and Environment, University of Padua, Padua, Italy
| | - Fabio Pilla
- Department of Agriculture Environment and Food Science, University of Molise, Campobasso, Italy
| | - Matteo Komjanc
- Unità risorse foraggere e produzioni zootecniche, Centro Trasferimento Tecnologico, Fondazione Edmund Mach, S. Michele all'Adige, Trento, Italy
| | - Alessandro Achilli
- Dipartimento di Biologia e Biotecnologie "L. Spallanzani", University of Pavia, Pavia, Italy
| | - Heidi Christine Hauffe
- Conservation Genomics Research Unit, Research and Innovation Centre, Fondazione Edmund Mach, S. Michele all'Adige, Trento, Italy
| | - Paolo Ajmone Marsan
- DIANA Dipartimento di Scienze Animali, della Nutrizione e degli Alimenti, Università Cattolica del Sacro Cuore, Piacenza, Italy
- Centro di Ricerca Nutrigenomica e Proteomica-PRONUTRIGEN, Università Cattolica del Sacro Cuore, Piacenza, Italy
| |
Collapse
|
3
|
Usai MG, Casu S, Sechi T, Salaris SL, Miari S, Mulas G, Cancedda MG, Ligios C, Carta A. 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. Genet Sel Evol 2024; 56:5. [PMID: 38200416 PMCID: PMC10777618 DOI: 10.1186/s12711-023-00873-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2023] [Accepted: 12/20/2023] [Indexed: 01/12/2024] Open
Abstract
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.
Collapse
Affiliation(s)
- Mario Graziano Usai
- Research Unit Genetics and Biotechnology - Agris Sardegna, 07100, Sassari, Italy
| | - Sara Casu
- Research Unit Genetics and Biotechnology - Agris Sardegna, 07100, Sassari, Italy.
| | - Tiziana Sechi
- Research Unit Genetics and Biotechnology - Agris Sardegna, 07100, Sassari, Italy
| | - Sotero L Salaris
- Research Unit Genetics and Biotechnology - Agris Sardegna, 07100, Sassari, Italy
| | - Sabrina Miari
- Research Unit Genetics and Biotechnology - Agris Sardegna, 07100, Sassari, Italy
| | - Giuliana Mulas
- Research Unit Genetics and Biotechnology - Agris Sardegna, 07100, Sassari, Italy
| | | | - Ciriaco Ligios
- Istituto Zooprofilattico Sperimentale Della Sardegna G. Pegreffi, 07100, Sassari, Italy
| | - Antonello Carta
- Research Unit Genetics and Biotechnology - Agris Sardegna, 07100, Sassari, Italy
| |
Collapse
|
4
|
Zhang X, Li Y, Han X, Xu Y, Wang H, Wang T, Zhang T. Role of the Killer Immunoglobulin-like Receptor and Human Leukocyte Antigen I Complex Polymorphisms in Kaposi Sarcoma-Associated Herpesvirus Infection. Open Forum Infect Dis 2023; 10:ofad435. [PMID: 37636520 PMCID: PMC10456215 DOI: 10.1093/ofid/ofad435] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/09/2023] [Accepted: 08/15/2023] [Indexed: 08/29/2023] Open
Abstract
Background Kaposi sarcoma, caused by the pathogen Kaposi sarcoma-associated herpesvirus (KSHV), is the most common neoplasm for patients with AIDS. Susceptibility to KSHV has been associated with several different genetic risk variants. The purpose of this study was to test whether variants of killer cell immunoglobulin-like receptors (KIRs) and their human leukocyte antigen (HLA-I) ligands influence the risk of KSHV infection. Methods A case-control study was performed in Xinjiang, a KSHV-endemic region of China. We recruited 299 individuals with HIV, including 123 KSHV-seropositive persons and 176 KSHV-seronegative controls. We used logistic regression and the MiDAS package to evaluate the association between KIR/HLA-I polymorphisms and KSHV infection. Results HLA-A*31:01, HLA-C*03:04, and HLA-C*12:03 were found to be associated with KSHV infection, with A*31:01 showing a protective effect under 3 different models (dominant: 0.30 [95% confidence interval {CI}, .08-.82], P = .031; additive: 0.30 [95% CI, .09-.80], P = .030; overdominant: 0.31 [95% CI, .09-.88], P = .042). The effect of A*31:01 might cause the variants of amino acid at HLA-A position 56, with individuals carrying an arginine having a lower KSHV infection risk. The increased homozygous KIR2DL3 was associated with a relatively high KSHV viral load (16.30% vs 41.94%, P = .010). Conclusions This study provides further insight into the link between HLA-I alleles and KIR genes and KSHV infection, highlighting KSHV-susceptible variants of HLA-I and KSHV replication caused by specific KIR genotype, and revealing a potential role of KIR-mediated natural killer cell activation in anti-KSHV infection.
Collapse
Affiliation(s)
- Xin Zhang
- Department of Epidemiology, School of Public Health, Fudan University, Shanghai, China
- Key Laboratory of Public Health Safety, Fudan University, Ministry of Education, Shanghai, China
| | - Yi Li
- Department of Epidemiology, School of Public Health, Fudan University, Shanghai, China
- Key Laboratory of Public Health Safety, Fudan University, Ministry of Education, Shanghai, China
| | - Xinyu Han
- Department of Epidemiology, School of Public Health, Fudan University, Shanghai, China
- Key Laboratory of Public Health Safety, Fudan University, Ministry of Education, Shanghai, China
| | - Yiyun Xu
- Department of Epidemiology, School of Public Health, Fudan University, Shanghai, China
- Key Laboratory of Public Health Safety, Fudan University, Ministry of Education, Shanghai, China
| | - Haili Wang
- Department of Epidemiology, School of Public Health, Fudan University, Shanghai, China
- Key Laboratory of Public Health Safety, Fudan University, Ministry of Education, Shanghai, China
| | - Tianye Wang
- Department of Epidemiology, School of Public Health, Fudan University, Shanghai, China
- Key Laboratory of Public Health Safety, Fudan University, Ministry of Education, Shanghai, China
| | - Tiejun Zhang
- Department of Epidemiology, School of Public Health, Fudan University, Shanghai, China
- Key Laboratory of Public Health Safety, Fudan University, Ministry of Education, Shanghai, China
- Shanghai Institute of Infectious Disease and Biosecurity, School of Public Health, Fudan University, Shanghai, China
- Yiwu Research Institute, Fudan University, Yiwu, China
| |
Collapse
|
5
|
Sanchez MP, Tribout T, Fritz S, Guatteo R, Fourichon C, Schibler L, Delafosse A, Boichard D. New insights into the genetic resistance to paratuberculosis in Holstein cattle via single-step genomic evaluation. Genet Sel Evol 2022; 54:67. [PMID: 36243688 PMCID: PMC9569073 DOI: 10.1186/s12711-022-00757-z] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2022] [Accepted: 09/24/2022] [Indexed: 11/10/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Marie-Pierre Sanchez
- Université Paris Saclay, INRAE, AgroParisTech, GABI, 78350, Jouy-en-Josas, France.
| | - Thierry Tribout
- Université Paris Saclay, INRAE, AgroParisTech, GABI, 78350, Jouy-en-Josas, France
| | - Sébastien Fritz
- Université Paris Saclay, INRAE, AgroParisTech, GABI, 78350, Jouy-en-Josas, France.,Eliance, 149 Rue de Bercy, 75012, Paris, France
| | | | | | | | | | - Didier Boichard
- Université Paris Saclay, INRAE, AgroParisTech, GABI, 78350, Jouy-en-Josas, France
| |
Collapse
|
6
|
Alonso-Hearn M, Badia-Bringué G, Canive M. Genome-wide association studies for the identification of cattle susceptible and resilient to paratuberculosis. Front Vet Sci 2022; 9:935133. [PMID: 36172612 PMCID: PMC9510743 DOI: 10.3389/fvets.2022.935133] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2022] [Accepted: 08/22/2022] [Indexed: 12/18/2022] Open
Abstract
Mycobacterium avium subsp. paratuberculosis (MAP) causes Johne's disease or paratuberculosis (PTB), with important animal health and economic implications. There are no therapeutic strategies to control this disease, and vaccination with inactivated vaccines is limited in many countries because it can interfere with the intradermal test used for bovine tuberculosis detection. Thus, infected animals either get culled after a positive ELISA or fecal PCR result or die due to clinical disease. In this study, we review recent studies aimed to discover genetic markers which could help to identify and select cattle less susceptible and more resilient to PTB. In recent years, the genotyping and subsequent imputation to whole-genome sequence (WGS) has allowed the identification of single-nucleotide polymorphisms (SNPs), quantitative trait loci (QTL), and candidate genes in the Bos taurus genome associated with susceptibility to MAP infection. In most of these genome-wide association studies (GWAS), phenotypes were based on ante-mortem test results including serum ELISA, milk ELISA, and detection of MAP by fecal PCR and bacteriological culture. Cattle infected with MAP display lesions with distinct severity but the associations between host genetics and PTB-associated pathology had not been explored until very recently. On the contrary, the understanding of the mechanisms and genetic loci influencing pathogen resistance, and disease tolerance in asymptomatic individuals is currently very limited. The identification of long-time asymptomatic cattle that is able to resist the infection and/or tolerate the disease without having their health and milk production compromised is important for disease control and breeding purposes.
Collapse
|
7
|
Heidari M, Pakdel A, Bakhtiarizadeh MR, Dehghanian F. A framework for non-preserved consensus gene module detection in Johne's disease. Front Vet Sci 2022; 9:974444. [PMID: 35968017 PMCID: PMC9363878 DOI: 10.3389/fvets.2022.974444] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2022] [Accepted: 07/08/2022] [Indexed: 11/29/2022] Open
Abstract
Johne's disease caused by Mycobacterium avium subsp. paratuberculosis (MAP) is a major concern in dairy industry. Since, the pathogenesis of the disease is not clearly known, it is necessary to develop an approach to discover molecular mechanisms behind this disease with high confidence. Biological studies often suffer from issues with reproducibility. Lack of a method to find stable modules in co-expression networks from different datasets related to Johne's disease motivated us to present a computational pipeline to identify non-preserved consensus modules. Two RNA-Seq datasets related to MAP infection were analyzed, and consensus modules were detected and were subjected to the preservation analysis. The non-preserved consensus modules in both datasets were determined as they are modules whose connectivity and density are affected by the disease. Long non-coding RNAs (lncRNAs) and TF genes in the non-preserved consensus modules were identified to construct integrated networks of lncRNA-mRNA-TF. These networks were confirmed by protein-protein interactions (PPIs) networks. Also, the overlapped hub genes between two datasets were considered hub genes of the consensus modules. Out of 66 consensus modules, 21 modules were non-preserved consensus modules, which were common in both datasets and 619 hub genes were members of these modules. Moreover, 34 lncRNA and 152 TF genes were identified in 12 and 19 non-preserved consensus modules, respectively. The predicted PPIs in 17 non-preserved consensus modules were significant, and 283 hub genes were commonly identified in both co-expression and PPIs networks. Functional enrichment analysis revealed that eight out of 21 modules were significantly enriched for biological processes associated with Johne's disease including “inflammatory response,” “interleukin-1-mediated signaling pathway”, “type I interferon signaling pathway,” “cytokine-mediated signaling pathway,” “regulation of interferon-beta production,” and “response to interferon-gamma.” Moreover, some genes (hub mRNA, TF, and lncRNA) were introduced as potential candidates for Johne's disease pathogenesis such as TLR2, NFKB1, IRF1, ATF3, TREM1, CDH26, HMGB1, STAT1, ISG15, CASP3. This study expanded our knowledge of molecular mechanisms involved in Johne's disease, and the presented pipeline enabled us to achieve more valid results.
Collapse
Affiliation(s)
- Maryam Heidari
- Department of Animal Science, College of Agriculture, Isfahan University of Technology, Isfahan, Iran
| | - Abbas Pakdel
- Department of Animal Science, College of Agriculture, Isfahan University of Technology, Isfahan, Iran
- *Correspondence: Abbas Pakdel
| | - Mohammad Reza Bakhtiarizadeh
- Department of Animal and Poultry Science, College of Aburaihan, University of Tehran, Tehran, Iran
- Mohammad Reza Bakhtiarizadeh
| | | |
Collapse
|
8
|
Kirkpatrick BW, Cooke ME, Frie M, Sporer KRB, Lett B, Wells SJ, Coussens PM. Genome-wide association analysis for susceptibility to infection by Mycobacterium avium ssp. paratuberculosis in US Holsteins. J Dairy Sci 2022; 105:4301-4313. [PMID: 35307176 DOI: 10.3168/jds.2021-21276] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2021] [Accepted: 01/27/2022] [Indexed: 12/19/2022]
Abstract
Paratuberculosis, or Johne's disease, is a chronic, granulomatous, gastrointestinal tract disease of cattle and other ruminants caused by the bacterium Mycobacterium avium subspecies paratuberculosis (MAP). Control of Johne's disease is based on programs of testing and culling animals positive for infection with MAP and concurrently modifying management to reduce the likelihood of infection. The current study was motivated by the hypothesis that genetic variation in host susceptibility to MAP infection can be dissected and quantifiable associations with genetic markers identified. Two separate GWAS analyses were conducted, the first using 897 genotyped Holstein artificial insemination sires with phenotypes derived from incidence of MAP infection among daughters based on milk ELISA testing records. The second GWAS analysis was a case-control design using US Holstein cows phenotyped for MAP infection by serum ELISA or fecal culture tests. Cases included cows positive for either serum ELISA, fecal culture, or both. Controls consisted of animals negative for all tests conducted. A total of 376 samples (70 cases and 306 controls) from a University of Minnesota Johne's management demonstration project and 184 samples (76 cases and 108 controls) from a Michigan State University study were used. Medium-density (sires) and high-density (cows) genotype data were imputed to full genome sequence for the analyses. Marker-trait associations were analyzed using the single-step (ss)GWAS procedure implemented in the BLUPF90 suite of programs. Evidence of significant genomic contributions for susceptibility to MAP infection were observed on multiple chromosomes. Results were combined across studies in a meta-analysis, and increased support for genomic regions on BTA7 and BTA21 were observed. Gene set enrichment analysis suggested pathways for antigen processing and presentation, antimicrobial peptides and natural killer cell-mediated cytotoxicity are relevant to variation in host susceptibility to MAP infection, among others. Genomic prediction was evaluated using a 5-fold cross-validation, and moderate correlations were observed between genomic breeding value predictions and daughter averages (∼0.43 to 0.53) for MAP infection in testing data sets. These results suggest that genomic selection against susceptibility to MAP infection is feasible in Holstein cattle.
Collapse
Affiliation(s)
- B W Kirkpatrick
- Department of Animal and Dairy Sciences, University of Wisconsin-Madison, 1675 Observatory Drive, Madison 53706.
| | - M E Cooke
- Department of Animal and Dairy Sciences, University of Wisconsin-Madison, 1675 Observatory Drive, Madison 53706
| | - M Frie
- Department of Animal Science, Michigan State University, 474 S Shaw Ln, East Lansing 48824
| | - K R B Sporer
- Department of Animal Science, Michigan State University, 474 S Shaw Ln, East Lansing 48824
| | - B Lett
- Department of Animal and Dairy Sciences, University of Wisconsin-Madison, 1675 Observatory Drive, Madison 53706
| | - S J Wells
- Department of Veterinary Population Medicine, University of Minnesota, 1365 Gortner Avenue, St. Paul 55108
| | - P M Coussens
- Department of Animal Science, Michigan State University, 474 S Shaw Ln, East Lansing 48824
| |
Collapse
|
9
|
Li H, Huang J, Zhang J, Gao Y, Han B, Sun D. Identification of Alternative Splicing Events Associated with Paratuberculosis in Dairy Cattle Using Multi-Tissue RNA Sequencing Data. Genes (Basel) 2022; 13:genes13030497. [PMID: 35328051 PMCID: PMC8948961 DOI: 10.3390/genes13030497] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2022] [Revised: 03/01/2022] [Accepted: 03/03/2022] [Indexed: 01/18/2023] Open
Abstract
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.
Collapse
|
10
|
Zhang X, Fang Q, Zhu S, Wu X, HuangboYuan, Liu Z, Xu Y, Chen T, Zeng Y, Zhang T. Environmental risk factors and genetic markers of Kaposi's sarcoma associated herpesvirus infection among Uygur population in Xinjiang, China. J Med Virol 2022; 94:2755-2765. [PMID: 35043408 DOI: 10.1002/jmv.27600] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2021] [Revised: 01/07/2022] [Accepted: 01/16/2022] [Indexed: 11/10/2022]
Abstract
OBJECTIVES Kaposi sarcoma-associated herpesvirus (KSHV) is endemic in Xinjiang, China. Determinants of KSHV seropositivity among high-risk groups are not well understood. We seek to identify genetic and environmental predisposing factors for KSHV infection among Uygurs in this endemic region. METHODS A cross-sectional study was performed among Uygur population in Xinjiang, China. KSHV-antibodies were detected using immunofluorescence assay (IFA) and human leukocyte antigen (HLA) alleles were genotyped. Univariate and multivariate logistic regression analysis were applied to explore the environmental and genetic risk factors of KSHV seropositivity. RESULTS A total of 721 participants were included. The seroprevalence of KSHV was 24.1% among this population. Sweet-food preference (OR 1.85, 95%CI 1.03-3.34), and coronary heart disease (OR 1.91, 95 %CI 1.24-2.94) were statistically correlated with KSHV infection. HLA-DQB1*06:09 were found to significantly increase the risk of KSHV infection under all 3 models (ORAllelic =4.06; ORDominant =3.27; ORRecessive =8.06). Six SNPs (SNP0260, SNP0361, SNP0797, SNP0852, SNP1159, SNP1375) in the DQB1 and DRB1 region and haploid type GTCTAACTAATC in block 17 were statistically associated with KSHV infection. CONCLUSION We demonstrated that genetic variations in HLA-DQB1/DRB1 and environmental risk factors were strongly associated with KSHV infection among this population. This article is protected by copyright. All rights reserved.
Collapse
Affiliation(s)
- Xin Zhang
- Department of Epidemiology, School of Public Health, Fudan University, Shanghai, China. Key Laboratory of Public Health Safety (Fudan University), Ministry of Education, China.,Shanghai Institute of Infectious Disease and Biosecurity, School of Public Health, Fudan University, Shanghai, China
| | - Qiwen Fang
- Department of Epidemiology, School of Public Health, Fudan University, Shanghai, China. Key Laboratory of Public Health Safety (Fudan University), Ministry of Education, China.,Shanghai Institute of Infectious Disease and Biosecurity, School of Public Health, Fudan University, Shanghai, China
| | - Sibo Zhu
- State Key Laboratory of Genetic Engineering and Collaborative Innovation Center for Genetics and Development, School of Life Sciences, Fudan University, Shanghai, 200438, China
| | - Xuefu Wu
- Department of Epidemiology, School of Public Health, Fudan University, Shanghai, China. Key Laboratory of Public Health Safety (Fudan University), Ministry of Education, China
| | - HuangboYuan
- Department of Epidemiology, School of Public Health, Fudan University, Shanghai, China. Key Laboratory of Public Health Safety (Fudan University), Ministry of Education, China
| | - Zhenqiu Liu
- State Key Laboratory of Genetic Engineering and Collaborative Innovation Center for Genetics and Development, School of Life Sciences, Fudan University, Shanghai, 200438, China
| | - Yiyun Xu
- Department of Epidemiology, School of Public Health, Fudan University, Shanghai, China. Key Laboratory of Public Health Safety (Fudan University), Ministry of Education, China
| | - Tao Chen
- Xinjiang Ili center for diseases control and prevention, Xinjiang, China
| | - Yan Zeng
- Key Laboratory of Xinjiang Endemic and Ethnic Disease and Department of Biochemistry, School of Medicine, Shihezi University, 832000, Shihezi, Xinjiang, China
| | - Tiejun Zhang
- Department of Epidemiology, School of Public Health, Fudan University, Shanghai, China. Key Laboratory of Public Health Safety (Fudan University), Ministry of Education, China.,Shanghai Institute of Infectious Disease and Biosecurity, School of Public Health, Fudan University, Shanghai, China.,Yiwu Research Institue, Fudan University, Yiwu, China
| |
Collapse
|
11
|
Ibeagha-Awemu EM, Bissonnette N, Bhattarai S, Wang M, Dudemaine PL, McKay S, Zhao X. Whole Genome Methylation Analysis Reveals Role of DNA Methylation in Cow's Ileal and Ileal Lymph Node Responses to Mycobacterium avium subsp. paratuberculosis Infection. Front Genet 2021; 12:797490. [PMID: 34992636 PMCID: PMC8724574 DOI: 10.3389/fgene.2021.797490] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/18/2021] [Accepted: 11/22/2021] [Indexed: 12/27/2022] Open
Abstract
Johne's Disease (JD), caused by Mycobacterium avium subsp paratuberculosis (MAP), is an incurable disease of ruminants and other animal species and is characterized by an imbalance of gut immunity. The role of MAP infection on the epigenetic modeling of gut immunity during the progression of JD is still unknown. This study investigated the DNA methylation patterns in ileal (IL) and ileal lymph node (ILLN) tissues from cows diagnosed with persistent subclinical MAP infection over a one to 4 years period. DNA samples from IL and ILLN tissues from cows negative (MAPneg) (n = 3) or positive for MAP infection (MAPinf) (n = 4) were subjected to whole genome bisulfite sequencing. A total of 11,263 and 62,459 differentially methylated cytosines (DMCs), and 1259 and 8086 differentially methylated regions (DMRs) (FDR<0.1) were found between MAPinf and MAPneg IL and ILLN tissues, respectively. The DMRs were found on 394 genes (denoted DMR genes) in the IL and on 1305 genes in the ILLN. DMR genes with hypermethylated promoters/5'UTR [3 (IL) and 88 (ILLN)] or hypomethylated promoters/5'UTR [10 (IL) and 25 (ILLN)] and having multiple functions including response to stimulus/immune response (BLK, BTC, CCL21, AVPR1A, CHRNG, GABRA4, TDGF1), cellular processes (H2AC20, TEX101, GLA, NCKAP5L, RBM27, SLC18A1, H2AC20BARHL2, NLGN3, SUV39H1, GABRA4, PPA1, UBE2D2) and metabolic processes (GSTO2, H2AC20, SUV39H1, PPA1, UBE2D2) are potential DNA methylation candidate genes of MAP infection. The ILLN DMR genes were enriched for more biological process (BP) gene ontology (GO) terms (n = 374), most of which were related to cellular processes (27.6%), biological regulation (16.6%), metabolic processes (15.4%) and response to stimulus/immune response (8.2%) compared to 75 BP GO terms (related to cellular processes, metabolic processes and transport, and system development) enriched for IL DMR genes. ILLN DMR genes were enriched for more pathways (n = 47) including 13 disease pathways compared with 36 enriched pathways, including 7 disease/immune pathways for IL DMR genes. In conclusion, the results show tissue specific responses to MAP infection with more epigenetic changes (DMCs and DMRs) in the ILLN than in the IL tissue, suggesting that the ILLN and immune processes were more responsive to regulation by methylation of DNA relative to IL tissue. Our data is the first to demonstrate a potential role for DNA methylation in the pathogenesis of MAP infection in dairy cattle.
Collapse
Affiliation(s)
- Eveline M. Ibeagha-Awemu
- Sherbrooke Research and Development Centre, Agriculture and Agri-Food Canada, Sherbrooke, QC, Canada
| | - Nathalie Bissonnette
- Sherbrooke Research and Development Centre, Agriculture and Agri-Food Canada, Sherbrooke, QC, Canada
| | - Suraj Bhattarai
- Department of Animal and Veterinary Sciences, University of Vermont, Burlington, VT, United States
| | - Mengqi Wang
- Sherbrooke Research and Development Centre, Agriculture and Agri-Food Canada, Sherbrooke, QC, Canada
| | - Pier-Luc Dudemaine
- Sherbrooke Research and Development Centre, Agriculture and Agri-Food Canada, Sherbrooke, QC, Canada
| | - Stephanie McKay
- Department of Animal and Veterinary Sciences, University of Vermont, Burlington, VT, United States
| | - Xin Zhao
- Department of Animal Science, McGill University, Ste-Anne-Be-Bellevue, QC, Canada
| |
Collapse
|
12
|
Zhang J, Han B, Zheng W, Lin S, Li H, Gao Y, Sun D. Genome-Wide DNA Methylation Profile in Jejunum Reveals the Potential Genes Associated With Paratuberculosis in Dairy Cattle. Front Genet 2021; 12:735147. [PMID: 34721525 PMCID: PMC8554095 DOI: 10.3389/fgene.2021.735147] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/02/2021] [Accepted: 09/23/2021] [Indexed: 12/04/2022] Open
Abstract
Paratuberculosis in cattle causes substantial economic losses to the dairy industry. Exploring functional genes and corresponding regulatory pathways related to resistance or susceptibility to paratuberculosis is essential to the breeding of disease resistance in cattle. Co-analysis of genome-wide DNA methylation and transcriptome profiles is a critically important approach to understand potential regulatory mechanism underlying the development of diseases. In this study, we characterized the profiles of DNA methylation of jejunum from nine Holstein cows in clinical, subclinical, and healthy groups using whole-genome bisulfite sequencing (WGBS). The average methylation level in functional regions was 29.95% in the promoter, 29.65% in the 5’ untranslated region (UTR), 68.24% in exons, 71.55% in introns, and 72.81% in the 3’ UTR. A total of 3,911, 4,336, and 4,094 differentially methylated genes (DMGs) were detected in clinical vs. subclinical, clinical vs. healthy, and subclinical vs. healthy comparative group, respectively. Gene ontology (GO) and analysis based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) showed that these DMGs were significantly enriched in specific biological processes related to immune response, such as Th1 and Th2 cell differentiation, wnt, TNF, MAPK, ECM-receptor interaction, cellular senescence, calcium, and chemokine signaling pathways (q value <0.05). The integration of information about DMGs, differentially expressed genes (DEGs), and biological functions suggested nine genes CALCRL, TNC, GATA4, CD44, TGM3, CXCL9, CXCL10, PPARG, and NFATC1 as promising candidates related to resistance/susceptibility to Mycobacterium avium subspecies paratuberculosis (MAP). This study reports on the high-resolution DNA methylation landscapes of the jejunum methylome across three conditions (clinical, subclinical, and healthy) in dairy cows. Our investigations integrated different sources of information about DMGs, DEGs, and pathways, enabling us to find nine functional genes that might have potential application in resisting paratuberculosis in dairy cattle.
Collapse
Affiliation(s)
- Junnan Zhang
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Ministry of Agriculture and Rural Affairs, National Engineering Laboratory for Animal Breeding, Department of Animal Genetics, Breeding and Reproduction, College of Animal Science and Technology, China Agricultural University, Beijing, China
| | - Bo Han
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Ministry of Agriculture and Rural Affairs, National Engineering Laboratory for Animal Breeding, Department of Animal Genetics, Breeding and Reproduction, College of Animal Science and Technology, China Agricultural University, Beijing, China
| | - Weijie Zheng
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Ministry of Agriculture and Rural Affairs, National Engineering Laboratory for Animal Breeding, Department of Animal Genetics, Breeding and Reproduction, College of Animal Science and Technology, China Agricultural University, Beijing, China
| | - Shan Lin
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Ministry of Agriculture and Rural Affairs, National Engineering Laboratory for Animal Breeding, Department of Animal Genetics, Breeding and Reproduction, College of Animal Science and Technology, China Agricultural University, Beijing, China
| | - Houcheng Li
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Ministry of Agriculture and Rural Affairs, National Engineering Laboratory for Animal Breeding, Department of Animal Genetics, Breeding and Reproduction, College of Animal Science and Technology, China Agricultural University, Beijing, China
| | - Yahui Gao
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Ministry of Agriculture and Rural Affairs, National Engineering Laboratory for Animal Breeding, Department of Animal Genetics, Breeding and Reproduction, College of Animal Science and Technology, China Agricultural University, Beijing, China
| | - Dongxiao Sun
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Ministry of Agriculture and Rural Affairs, National Engineering Laboratory for Animal Breeding, Department of Animal Genetics, Breeding and Reproduction, College of Animal Science and Technology, China Agricultural University, Beijing, China
| |
Collapse
|
13
|
Liu D, Chen Z, Zhao W, Guo L, Sun H, Zhu K, Liu G, Shen X, Zhao X, Wang Q, Ma P, Pan Y. Genome-wide selection signatures detection in Shanghai Holstein cattle population identified genes related to adaption, health and reproduction traits. BMC Genomics 2021; 22:747. [PMID: 34654366 PMCID: PMC8520274 DOI: 10.1186/s12864-021-08042-x] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2021] [Accepted: 09/28/2021] [Indexed: 12/26/2022] Open
Abstract
Background Over several decades, a wide range of natural and artificial selection events in response to subtropical environments, intensive pasture and intensive feedlot systems have greatly changed the customary behaviour, appearance, and important economic traits of Shanghai Holstein cattle. In particular, the longevity of the Shanghai Holstein cattle population is generally short, approximately the 2nd to 3rd lactation. In this study, two complementary approaches, integrated haplotype score (iHS) and runs of homozygosity (ROH), were applied for the detection of selection signatures within the genome using genotyping by genome-reduced sequence data from 1092 cows. Results In total, 101 significant iHS genomic regions containing selection signatures encompassing a total of 256 candidate genes were detected. There were 27 significant |iHS| genomic regions with a mean |iHS| score > 2. The average number of ROH per individual was 42.15 ± 25.47, with an average size of 2.95 Mb. The length of 78 % of the detected ROH was within the range of 1–2 MB and 2–4 MB, and 99 % were shorter than 8 Mb. A total of 168 genes were detected in 18 ROH islands (top 1 %) across 16 autosomes, in which each SNP showed a percentage of occurrence > 30 %. There were 160 and 167 genes associated with the 52 candidate regions within health-related QTL intervals and 59 candidate regions within reproduction-related QTL intervals, respectively. Annotation of the regions harbouring clustered |iHS| signals and candidate regions for ROH revealed a panel of interesting candidate genes associated with adaptation and economic traits, such as IL22RA1, CALHM3, ITGA9, NDUFB3, RGS3, SOD2, SNRPA1, ST3GAL4, ALAD, EXOSC10, and MASP2. In a further step, a total of 1472 SNPs in 256 genes were matched with 352 cis-eQTLs in 21 tissues and 27 trans-eQTLs in 6 tissues. For SNPs located in candidate regions for ROH, a total of 108 cis-eQTLs in 13 tissues and 4 trans-eQTLs were found for 1092 SNPs. Eighty-one eGenes were significantly expressed in at least one tissue relevant to a trait (P value < 0.05) and matched the 256 genes detected by iHS. For the 168 significant genes detected by ROH, 47 gene-tissue pairs were significantly associated with at least one of the 37 traits. Conclusions We provide a comprehensive overview of selection signatures in Shanghai Holstein cattle genomes by combining iHS and ROH. Our study provides a list of genes associated with immunity, reproduction and adaptation. For functional annotation, the cGTEx resource was used to interpret SNP-trait associations. The results may facilitate the identification of genes relevant to important economic traits and can help us better understand the biological processes and mechanisms affected by strong ongoing natural or artificial selection in livestock populations. Supplementary Information The online version contains supplementary material available at 10.1186/s12864-021-08042-x.
Collapse
Affiliation(s)
- Dengying Liu
- Department of Animal Science, School of Agriculture and Biology, Shanghai Jiao Tong University, 200240, Shanghai, PR China
| | - Zhenliang Chen
- Department of Animal Science, School of Agriculture and Biology, Shanghai Jiao Tong University, 200240, Shanghai, PR China
| | - Wei Zhao
- Department of Animal Science, School of Agriculture and Biology, Shanghai Jiao Tong University, 200240, Shanghai, PR China
| | - Longyu Guo
- Department of Animal Science, School of Agriculture and Biology, Shanghai Jiao Tong University, 200240, Shanghai, PR China
| | - Hao Sun
- Department of Animal Science, School of Agriculture and Biology, Shanghai Jiao Tong University, 200240, Shanghai, PR China
| | - Kai Zhu
- Shanghai Dairy Cattle Breeding Centre Co., Ltd, 201901, Shanghai, P.R. China
| | - Guanglei Liu
- Shanghai Dairy Cattle Breeding Centre Co., Ltd, 201901, Shanghai, P.R. China
| | - Xiuping Shen
- Shanghai Agricultural Development Promotion Center, 200335, Shanghai, PR China
| | - Xiaoduo Zhao
- Shanghai Dairy Cattle Breeding Centre Co., Ltd, 201901, Shanghai, P.R. China
| | - Qishan Wang
- Department of Animal Breeding and Reproduction, College of Animal Science, Zhejiang University, 310058, Hangzhou, PR China
| | - Peipei Ma
- Department of Animal Science, School of Agriculture and Biology, Shanghai Jiao Tong University, 200240, Shanghai, PR China.
| | - Yuchun Pan
- Department of Animal Breeding and Reproduction, College of Animal Science, Zhejiang University, 310058, Hangzhou, PR China.
| |
Collapse
|
14
|
Canive M, Badia-Bringué G, Vázquez P, González-Recio O, Fernández A, Garrido JM, Juste RA, Alonso-Hearn M. Identification of loci associated with pathological outcomes in Holstein cattle infected with Mycobacterium avium subsp. paratuberculosis using whole-genome sequence data. Sci Rep 2021; 11:20177. [PMID: 34635747 PMCID: PMC8505495 DOI: 10.1038/s41598-021-99672-4] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2021] [Accepted: 09/27/2021] [Indexed: 12/03/2022] Open
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.
Collapse
Affiliation(s)
- Maria Canive
- Department of Animal Health, NEIKER-Basque Institute for Agricultural Research and Development, Basque Research and Technology Alliance (BRTA), Derio, Bizkaia, Spain.,Doctoral Program in Immunology, Microbiology and Parasitology, Universidad del País Vasco/Euskal Herriko Unibertsitatea (UPV/EHU), Leioa, Bizkaia, Spain
| | - Gerard Badia-Bringué
- Department of Animal Health, NEIKER-Basque Institute for Agricultural Research and Development, Basque Research and Technology Alliance (BRTA), Derio, Bizkaia, Spain.,Doctoral Program in Molecular Biology and Biomedicine, Universidad del País Vasco/Euskal Herriko Unibertsitatea (UPV/EHU), Leioa, Bizkaia, Spain
| | - Patricia Vázquez
- Department of Animal Health, NEIKER-Basque Institute for Agricultural Research and Development, Basque Research and Technology Alliance (BRTA), Derio, Bizkaia, Spain
| | - Oscar González-Recio
- Departamento de Mejora Genética Animal, Instituto Nacional de Investigación y Tecnología Agraria y Alimentaria, CSIC, Madrid, Spain.,Departamento de Producción Agraria, Escuela Técnica Superior de Ingeniería Agronómica, Alimentaria y de Biosistemas, Universidad Politécnica de Madrid, Ciudad Universitaria, Madrid, Spain
| | - Almudena Fernández
- Departamento de Mejora Genética Animal, Instituto Nacional de Investigación y Tecnología Agraria y Alimentaria, CSIC, Madrid, Spain
| | - Joseba M Garrido
- Department of Animal Health, NEIKER-Basque Institute for Agricultural Research and Development, Basque Research and Technology Alliance (BRTA), Derio, Bizkaia, Spain
| | - Ramón A Juste
- Department of Animal Health, NEIKER-Basque Institute for Agricultural Research and Development, Basque Research and Technology Alliance (BRTA), Derio, Bizkaia, Spain
| | - Marta Alonso-Hearn
- Department of Animal Health, NEIKER-Basque Institute for Agricultural Research and Development, Basque Research and Technology Alliance (BRTA), Derio, Bizkaia, Spain.
| |
Collapse
|
15
|
Mallikarjunappa S, Brito LF, Pant SD, Schenkel FS, Meade KG, Karrow NA. Johne's Disease in Dairy Cattle: An Immunogenetic Perspective. Front Vet Sci 2021; 8:718987. [PMID: 34513975 PMCID: PMC8426623 DOI: 10.3389/fvets.2021.718987] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/01/2021] [Accepted: 07/19/2021] [Indexed: 12/25/2022] Open
Abstract
Johne's disease (JD), also known as paratuberculosis, is a severe production-limiting disease with significant economic and welfare implications for the global cattle industry. Caused by infection with Mycobacterium avium subspecies paratuberculosis (MAP), JD manifests as chronic enteritis in infected cattle. In addition to the economic losses and animal welfare issues associated with JD, MAP has attracted public health concerns with potential association with Crohn's disease, a human inflammatory bowel disease. The lack of effective treatment options, such as a vaccine, has hampered JD control resulting in its increasing global prevalence. The disease was first reported in 1895, but in recognition of its growing economic impact, extensive recent research facilitated by a revolution in technological approaches has led to significantly enhanced understanding of the immunological, genetic, and pathogen factors influencing disease pathogenesis. This knowledge has been derived from a variety of diverse models to elucidate host-pathogen interactions including in vivo and in vitro experimental infection models, studies measuring immune parameters in naturally-infected animals, and by studies conducted at the population level to enable the estimation of genetic parameters, and the identification of genetic markers and quantitative trait loci (QTL) putatively associated with susceptibility or resistance to JD. The main objectives of this review are to summarize these recent developments from an immunogenetics perspective and attempt to extract the principal and common findings emerging from this wealth of recent information. Based on these analyses, and in light of emerging technologies such as gene-editing, we conclude by discussing potential future avenues for effectively mitigating JD in cattle.
Collapse
Affiliation(s)
- Sanjay Mallikarjunappa
- Department of Animal Biosciences, Centre for Genetic Improvement of Livestock, University of Guelph, Guelph, ON, Canada
| | - Luiz F Brito
- Department of Animal Sciences, Purdue University, West Lafayette, IN, United States
| | - Sameer D Pant
- Graham Centre for Agricultural Innovation, Charles Sturt University, Wagga Wagga, NSW, Australia
| | - Flavio S Schenkel
- Department of Animal Biosciences, Centre for Genetic Improvement of Livestock, University of Guelph, Guelph, ON, Canada
| | - Kieran G Meade
- School of Agriculture and Food Science, University College Dublin, Dublin, Ireland
| | - Niel A Karrow
- Department of Animal Biosciences, Centre for Genetic Improvement of Livestock, University of Guelph, Guelph, ON, Canada
| |
Collapse
|
16
|
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. PLoS One 2021; 16:e0256091. [PMID: 34449805 PMCID: PMC8396740 DOI: 10.1371/journal.pone.0256091] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/28/2020] [Accepted: 07/29/2021] [Indexed: 12/16/2022] Open
Abstract
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.
Collapse
|
17
|
Kravitz A, Pelzer K, Sriranganathan N. The Paratuberculosis Paradigm Examined: A Review of Host Genetic Resistance and Innate Immune Fitness in Mycobacterium avium subsp. Paratuberculosis Infection. Front Vet Sci 2021; 8:721706. [PMID: 34485444 PMCID: PMC8414637 DOI: 10.3389/fvets.2021.721706] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2021] [Accepted: 07/23/2021] [Indexed: 11/13/2022] Open
Abstract
Paratuberculosis, or Johne's Disease (JD) is a debilitating chronic enteritis mainly affecting ruminants caused by Mycobacterium avium subsp. paratuberculosis (MAP). This organism causes worldwide economic losses to the livestock industry, and is of public health importance due to the potential zoonotic risk between MAP and Crohn's disease (CD) in humans. Without economical treatments, or a vaccine capable of preventing infection without causing cross-reactions with bovine tuberculosis, test-and-cull methods for disease control are imperative. Unfortunately, difficulties in diagnostics and long subclinical stage hinder adequate control and is further complicated by variation in MAP exposure outcome. Interestingly, the majority of infections result in asymptomatic presentation and never progress to clinical disease. One contributing factor is host genetics, where polymorphisms in innate immune genes have been found to influence resistance and susceptibility to disease. Candidate genes identified across studies overlap with those found in CD and tuberculosis including; Solute carrier family 11 member 1 gene (SLC11A1), Nucleotide-binding-oligomerization domain containing gene 2 (NOD2), Major histocompatibility complex type II (MHC-II), and Toll-like receptor (TLR) genes. This review will highlight evidence supporting the vital role of these genes in MAP infection outcome, associated challenges, and implications for the future of JD research.
Collapse
Affiliation(s)
- Amanda Kravitz
- Department of Biomedical Sciences and Pathobiology, Center for One Health Research, Virginia-Maryland College of Veterinary Medicine, Virginia Polytechnic Institute and State University, Blacksburg, VA, United States
| | - Kevin Pelzer
- Department of Large Animal Clinical Sciences, Virginia-Maryland College of Veterinary Medicine, Virginia Polytechnic Institute and State University, Blacksburg, VA, United States
| | - Nammalwar Sriranganathan
- Department of Biomedical Sciences and Pathobiology, Center for One Health Research, Virginia-Maryland College of Veterinary Medicine, Virginia Polytechnic Institute and State University, Blacksburg, VA, United States
| |
Collapse
|
18
|
Kiser JN, Wang Z, Zanella R, Scraggs E, Neupane M, Cantrell B, Van Tassell CP, White SN, Taylor JF, Neibergs HL. Functional Variants Surrounding Endothelin 2 Are Associated With Mycobacterium avium Subspecies paratuberculosis Infection. Front Vet Sci 2021; 8:625323. [PMID: 34026885 PMCID: PMC8131860 DOI: 10.3389/fvets.2021.625323] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2020] [Accepted: 03/04/2021] [Indexed: 02/04/2023] Open
Abstract
Bovine paratuberculosis, caused by Mycobacterium avium subspecies paratuberculosis (MAP), continues to impact the dairy industry through increased morbidity, mortality, and lost production. Although genome-wide association analyses (GWAAs) have identified loci associated with susceptibility to MAP, limited progress has been made in identifying mutations that cause disease susceptibility. A 235-kb region on Bos taurus chromosome 3 (BTA3), containing a 70-kb haplotype block surrounding endothelin 2 (EDN2), has previously been associated with the risk of MAP infection. EDN2 is highly expressed in the gut and is involved in intracellular calcium signaling and a wide array of biological processes. The objective of this study was to identify putative causal mutations for disease susceptibility in the region surrounding EDN2 in Holstein and Jersey cattle. Using sequence data from 10 Holstein and 10 Jersey cattle, common variants within the 70-kb region containing EDN2 were identified. A custom SNP genotyping array fine-mapped the region using 221 Holstein and 51 Jersey cattle and identified 17 putative causal variants (P < 0.01) located in the 5′ region of EDN2 and a SNP in the 3′ UTR (P = 0.00009) associated with MAP infection. MicroRNA interference assays, mRNA stability assays, and electrophoretic mobility shift assays were performed to determine if allelic changes at each SNP resulted in differences in EDN2 stability or expression. Two SNPs [rs109651404 (G/A) and rs110287192 (G/T)] located within the promoter region of EDN2 displayed differential binding affinity for transcription factors in binding sequences harboring the alternate SNP alleles. The luciferase reporter assay revealed that the transcriptional activity of the EDN2 promoter was increased (P < 0.05) with the A allele for rs109651404 and the G allele for rs110287192. These results suggest that the variants rs109651404 and rs110287192 are mutations that alter transcription and thus may alter susceptibility to MAP infection in Holstein and Jersey cattle.
Collapse
Affiliation(s)
- Jennifer N Kiser
- Department of Animal Sciences, Washington State University, Pullman, WA, United States
| | - Zeping Wang
- Department of Animal Sciences, Washington State University, Pullman, WA, United States
| | - Ricardo Zanella
- Department of Animal Sciences, Washington State University, Pullman, WA, United States
| | - Erik Scraggs
- Department of Animal Sciences, Washington State University, Pullman, WA, United States
| | - Mahesh Neupane
- Department of Animal Sciences, Washington State University, Pullman, WA, United States
| | - Bonnie Cantrell
- Department of Animal Sciences, Washington State University, Pullman, WA, United States
| | - Curtis P Van Tassell
- Animal Genomics and Improvement Laboratory, United States Department of Agriculture, Agricultural Research Service, Beltsville, MD, United States
| | - Stephen N White
- Animal Disease Research, United States Department of Agriculture, Agricultural Research Service, Pullman, WA, United States.,Department of Veterinary Microbiology and Pathology, Washington State University, Pullman, WA, United States.,Center for Reproductive Biology, Washington State University, Pullman, WA, United States
| | - Jeremy F Taylor
- Division of Animal Sciences, University of Missouri, Columbia, MO, United States
| | - Holly L Neibergs
- Department of Animal Sciences, Washington State University, Pullman, WA, United States
| |
Collapse
|
19
|
Identification of loci associated with susceptibility to bovine paratuberculosis and with the dysregulation of the MECOM, eEF1A2, and U1 spliceosomal RNA expression. Sci Rep 2021; 11:313. [PMID: 33432064 PMCID: PMC7801378 DOI: 10.1038/s41598-020-79619-x] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2020] [Accepted: 12/07/2020] [Indexed: 12/13/2022] Open
Abstract
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.
Collapse
|
20
|
Hasenauer FC, Rossi UA, Caffaro ME, Raschia MA, Maurizio E, Poli MA, Rossetti CA. Association of TNF rs668920841 and INRA111 polymorphisms with caprine brucellosis: A case-control study of candidate genes involved in innate immunity. Genomics 2020; 112:3925-3932. [PMID: 32629097 DOI: 10.1016/j.ygeno.2020.06.050] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/12/2020] [Revised: 05/25/2020] [Accepted: 06/28/2020] [Indexed: 01/24/2023]
Abstract
Caprine brucellosis is an infectious, contagious zoonotic disease caused by Brucella melitensis. Multiple factors, including host genetics, can influence the outcome of the exposure to Brucella; and it is expected that genetic variants that affect the host innate immune response could have a key role in Brucella infection and pathogenesis. In this study, we evaluated if polymorphisms in innate immunity-related genes are associated with results of Brucella infection in goats. Nine polymorphisms within interferon gamma (IFNG), tumor necrosis factor (TNF), MyD88 innate immune signal transduction adaptor (MYD88), interleukin 10 (IL10) and IL-10 receptor subunit alpha (IL10RA) genes and two molecular markers (BMS2753 and INRA111) were resolved by PCR-capillary electrophoresis in samples from 81 seronegative and 61 seropositive goats for brucellosis. A heterozygous genotype at INRA111, a microsatellite near the VRK serine/threonine kinase 2 (VRK2) gene, was associated with absence of Brucella-specific antibodies in goats naturally exposed to the pathogen (P = .004). Conversely, variants in the TNF gene (rs668920841) and near the IFN gamma receptor 1 (IFNGR1) gene (microsatellite BMS2753) were significantly associated with presence of Brucella-specific antibodies at allelic (P = .042 and P = .046) and genotypic level (P = .012 and P = .041, respectively). Moreover, an in silico analysis predicted a functional role of the insertion-deletion polymorphism rs668920841 on the transcriptional regulation of the caprine TNF gene. Altogether, these results contribute to the identification of genetic factors that have a putative effect on the resistance / susceptibility phenotype of goats to Brucella infection.
Collapse
Affiliation(s)
- F C Hasenauer
- Instituto de Patobiología, CICVyA, INTA, Nicolás Repetto y de Los Reseros s/n, Hurlingham, B1686, Buenos Aires, Argentina; CONICET, Godoy Cruz 2290, C1425 CABA, Argentina
| | - U A Rossi
- Instituto de Patobiología, CICVyA, INTA, Nicolás Repetto y de Los Reseros s/n, Hurlingham, B1686, Buenos Aires, Argentina; CONICET, Godoy Cruz 2290, C1425 CABA, Argentina
| | - M E Caffaro
- Instituto de Genética "Ewald A. Favret", CICVyA, INTA, Nicolás Repetto y de Los Reseros s/n, Hurlingham, B1686, Buenos Aires, Argentina
| | - M A Raschia
- Instituto de Genética "Ewald A. Favret", CICVyA, INTA, Nicolás Repetto y de Los Reseros s/n, Hurlingham, B1686, Buenos Aires, Argentina
| | - E Maurizio
- Instituto de Patobiología, CICVyA, INTA, Nicolás Repetto y de Los Reseros s/n, Hurlingham, B1686, Buenos Aires, Argentina; CONICET, Godoy Cruz 2290, C1425 CABA, Argentina
| | - M A Poli
- Instituto de Genética "Ewald A. Favret", CICVyA, INTA, Nicolás Repetto y de Los Reseros s/n, Hurlingham, B1686, Buenos Aires, Argentina
| | - C A Rossetti
- Instituto de Patobiología, CICVyA, INTA, Nicolás Repetto y de Los Reseros s/n, Hurlingham, B1686, Buenos Aires, Argentina.
| |
Collapse
|
21
|
Sanchez MP, Guatteo R, Davergne A, Saout J, Grohs C, Deloche MC, Taussat S, Fritz S, Boussaha M, Blanquefort P, Delafosse A, Joly A, Schibler L, Fourichon C, Boichard D. Identification of the ABCC4, IER3, and CBFA2T2 candidate genes for resistance to paratuberculosis from sequence-based GWAS in Holstein and Normande dairy cattle. Genet Sel Evol 2020; 52:14. [PMID: 32183688 PMCID: PMC7077142 DOI: 10.1186/s12711-020-00535-9] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/05/2019] [Accepted: 03/11/2020] [Indexed: 12/29/2022] Open
Abstract
Background Bovine paratuberculosis is a contagious disease, caused by Mycobacterium avium subsp. paratuberculosis (MAP), with adverse effects on animal welfare and serious economic consequences. Published results on host genetic resistance to MAP are inconsistent, mainly because of difficulties in characterizing the infection status of cows. The objectives of this study were to identify quantitative trait loci (QTL) for resistance to MAP in Holstein and Normande cows with an accurately defined status for MAP. Results From MAP-infected herds, cows without clinical signs of disease were subjected to at least four repeated serum ELISA and fecal PCR tests over time to determine both infected and non-infected statuses. Clinical cases were confirmed using PCR. Only cows that had concordant results for all tests were included in further analyses. Positive and control cows were matched within herd according to their birth date to ensure a same level of exposure to MAP. Cows with accurate phenotypes, i.e. unaffected (control) or affected (clinical or non-clinical cases), were genotyped with the Illumina BovineSNP50 BeadChip. Genotypes were imputed to whole-genome sequences using the 1000 Bull Genomes reference population (run6). A genome-wide association study (GWAS) of MAP status of 1644 Holstein and 649 Normande cows, using either two (controls versus cases) or three classes of phenotype (controls, non-clinical and clinical cases), revealed three regions, on Bos taurus (BTA) chromosomes 12, 13, and 23, presenting significant effects in Holstein cows, while only one of those was identified in Normande cows (BTA23). The most significant effect was found on BTA13, in a short 8.5-kb region. Conditional analyses revealed that only one causal variant may be responsible for the effects observed on each chromosome with the ABCC4 (BTA12), CBFA2T2 (BTA13), and IER3 (BTA23) genes as good functional candidates. Conclusions A sequence-based GWAS on cows for which resistance to MAP was accurately defined, was able to identify candidate variants located in genes that were functionally related to resistance to MAP; these explained up to 28% of the genetic variance of the trait. These results are very encouraging for efforts towards implementation of a breeding strategy aimed at improving resistance to paratuberculosis in Holstein cows.
Collapse
Affiliation(s)
- Marie-Pierre Sanchez
- Université Paris Saclay, INRAE, AgroParisTech, GABI, 78350, Jouy-en-Josas, France.
| | | | | | - Judikael Saout
- Université Paris Saclay, INRAE, AgroParisTech, GABI, 78350, Jouy-en-Josas, France
| | - Cécile Grohs
- Université Paris Saclay, INRAE, AgroParisTech, GABI, 78350, Jouy-en-Josas, France
| | - Marie-Christine Deloche
- Université Paris Saclay, INRAE, AgroParisTech, GABI, 78350, Jouy-en-Josas, France.,Allice, 149 Rue de Bercy, 75012, Paris, France
| | - Sébastien Taussat
- Université Paris Saclay, INRAE, AgroParisTech, GABI, 78350, Jouy-en-Josas, France.,Allice, 149 Rue de Bercy, 75012, Paris, France
| | - Sébastien Fritz
- Université Paris Saclay, INRAE, AgroParisTech, GABI, 78350, Jouy-en-Josas, France.,Allice, 149 Rue de Bercy, 75012, Paris, France
| | - Mekki Boussaha
- Université Paris Saclay, INRAE, AgroParisTech, GABI, 78350, Jouy-en-Josas, France
| | | | | | | | | | | | - Didier Boichard
- Université Paris Saclay, INRAE, AgroParisTech, GABI, 78350, Jouy-en-Josas, France
| |
Collapse
|
22
|
Single-nucleotide polymorphisms in CLEC7A, CD209 and TLR4 gene and their association with susceptibility to paratuberculosis in Indian cattle. J Genet 2020. [DOI: 10.1007/s12041-019-1172-4] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
|