Functional Variants Surrounding Endothelin 2 Are Associated With Mycobacterium avium Subspecies paratuberculosis Infection

Jennifer N Kiser; Zeping Wang; Ricardo Zanella; Erik Scraggs; Mahesh Neupane; Bonnie Cantrell; Curtis P Van Tassell; Stephen N White; Jeremy F Taylor; Holly L Neibergs

doi:10.3389/fvets.2021.625323

Functional Variants Surrounding Endothelin 2 Are Associated With Mycobacterium avium Subspecies paratuberculosis Infection

Front Vet Sci. 2021 May 5:8:625323. doi: 10.3389/fvets.2021.625323. eCollection 2021.

Authors

Affiliations

¹ Department of Animal Sciences, Washington State University, Pullman, WA, United States.
² Animal Genomics and Improvement Laboratory, United States Department of Agriculture, Agricultural Research Service, Beltsville, MD, United States.
³ 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.
⁶ Division of Animal Sciences, University of Missouri, Columbia, MO, United States.

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.

Keywords: EDN2; dairy cattle; functional assays; functional variant; paratubercolosis.