Identification of the ABCC4, IER3, and CBFA2T2 candidate genes for resistance to paratuberculosis from sequence-based GWAS in Holstein and Normande dairy cattle.

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.

Longitudinal study of Mycobacterium avium ssp. paratuberculosis fecal shedding patterns and concurrent serological patterns in naturally infected dairy cattle.

Mycobacterium avium ssp. paratuberculosis (MAP) is the etiological agent of paratuberculosis, a disease that affects ruminants worldwide. Despite global interest in the control of this disease, gaps exist in our knowledge of fecal shedding patterns and concurrent serological patterns. This longitudinal study in dairy cattle herds with high MAP seroprevalence in France aimed at accurately describing fecal shedding patterns over 1 year; relating those shedding patterns to individual animal characteristics (age, breed, parity); and exploring the association between fecal shedding patterns and serological patterns. To describe temporal fecal shedding patterns and continuity of shedding, along with the standard quantitative PCR (qPCR) threshold cycle we used a cutoff value that related to low or nonculturable fecal shedding. We also defined a threshold cycle indicative of shedding in high quantities to describe infection progression patterns. Twenty-one herds completed the study, and 782 cows were tested 4 times each. We obtained 4 sets of paired fecal qPCR and serum ELISA results from 757 cows. Although we targeted highly likely infectious animals, we found a large diversity of shedding patterns, as well as high variability between herds in the proportion of animals showing a given pattern. The fecal qPCR results of almost 20% of the final study sample were positioned at least once in the range that indicated low or nonculturable fecal shedding (between the adjusted and the standard cutoff value). Although these animals would typically be classified as non-shedders, they could be important to infection dynamics on the farm. Animals that shed at least twice consecutively and animals that shed in high quantities rarely reverted to negativity. Repeated fecal qPCR can be used to detect temporal fecal shedding traits, and the decision to cull an animal could practically be based on temporal, semiquantitative results. Overall, we found a mismatch between fecal shedding and ELISA seropositivity (637 animals were ELISA-negative 4 times, but only 13% of those animals were qPCR-negative 4 times). We found that having more than 2 ELISA-positive samples was strongly related to persistent and continuous shedding. We suggest that although serological testing is much less sensitive than qPCR, it can also be used, particularly over the course of multiple testing events, to identify animals that are most likely to contribute to the contamination of the farm environment.