Paratuberculosis (Johne’s disease), an enteric disorder in ruminants caused by Mycobacterium avium subspecies paratuberculosis (MAP), causes economic losses in excess of $200 million annually to the US dairy industry. To identify genomic regions underlying susceptibility to MAP infection in Jersey cattle, a case-control genome-wide association study (GWAS) was performed. Blood and fecal samples were collected from ∼5,000 mature cows in 30 commercial Jersey herds from across the US. Discovery data consisted of 450 cases and 439 controls genotyped with the Illumina BovineSNP50 BeadChip. Cases were animals with positive ELISA and fecal culture (FC) results. Controls were animals negative to both ELISA and FC tests that matched cases on birth date and herd. Validation data consisted of 180 animals including 90 cases (positive to FC) and 90 controls (negative to ELISA and FC), selected from discovery herds and genotyped by Illumina BovineLD BeadChip (∼7K SNPs). Two analytical approaches were used: single-marker GWAS using the GRAMMAR-GC method and Bayesian variable selection (Bayes C) using GenSel software. GRAMMAR-GC identified one SNP on BTA7 at 68 megabases (Mb) surpassing a significance threshold of 5×10⁻⁵. ARS-BFGL-NGS-11887 on BTA23 (27.7 Mb) accounted for the highest percentage of genetic variance (3.3%) in the Bayes C analysis. SNPs identified in common by GRAMMAR-GC and Bayes C in both discovery and combined data were mapped to BTA23 (27, 29 and 44 Mb), 3 (100, 101, 106 and 107 Mb) and 17 (57 Mb). Correspondence between results of GRAMMAR-GC and Bayes C was high (70–80% of most significant SNPs in common). These SNPs could potentially be associated with causal variants underlying susceptibility to MAP infection in Jersey cattle. Predictive performance of the model developed by Bayes C for prediction of infection status of animals in validation set was low (55% probability of correct ranking of paired case and control samples).

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 while concurrently modifying management to reduce the likelihood of infection. The current study is motivated by the hypothesis that genetic variation in host susceptibility to MAP infection can be dissected and quantifiable associations with genetic markers identified. For this purpose, a case-control, genome-wide association study was conducted using US Holstein cattle phenotyped for MAP infection using a serum ELISA and/or fecal culture test. Cases included cows positive for either serum ELISA, fecal culture or both. Controls consisted of animals negative for the serum ELISA test or both serum ELISA and fecal culture when both were available. Controls were matched by herd and proximal birth date with cases. A total of 856 cows (451 cases and 405 controls) were used in initial discovery analyses, and an additional 263 cows (159 cases and 104 controls) from the same herds were used as a validation data set. Data were analyzed in a single marker analysis controlling for relatedness of individuals (GRAMMAR-GC) and also in a Bayesian analysis in which multiple marker effects were estimated simultaneously (GenSel). For the latter, effects of non-overlapping 1 Mb marker windows across the genome were estimated. Results from the two discovery analyses were generally concordant; however, discovery results were generally not well supported in analysis of the validation data set. A combined analysis of discovery and validation data sets provided strongest support for SNPs and 1 Mb windows on chromosomes 1, 2, 6, 7, 17 and 29.

Half-sib daughters sired by a bull believed to be a carrier of a major gene for high ovulation rate were evaluated for ovulation rate and genotyped in an effort to both test the hypothesis of segregation of a major gene and to map the gene’s location. A total of 131 daughters were produced over four consecutive years at a University of Wisconsin-Madison research farm. All were evaluated for ovulation rate over an average of four estrous cycles using transrectal ultrasonography. The sire and all daughters were genotyped using a 3K SNP chip and the genotype and phenotype data were used in a linkage analysis. Subsequently, daughters recombinant within the QTL region and the sire were genotyped successively with 50K and 777K SNP chips to refine the location of the causative polymorphism. Positional candidate genes within the fine-mapped region were examined for polymorphism by Sanger sequencing of PCR amplicons encompassing coding and 5’ and 3’ flanking regions of the genes. Sire DNA was used as template in the PCR reactions. Strong evidence of a major gene for ovulation rate was observed (p<1x10^-28) with the gene localized to bovine chromosome 10. Fine-mapping subsequently reduced the location to a 1.2 Mb region between 13.6 and 14.8 Mb on chromosome 10. The location identified does not correspond to that for any previously identified major gene for ovulation rate. This region contains three candidate genes, SMAD3, SMAD6 and IQCH. While candidate gene screening failed to identify the causative polymorphism, three polymorphisms were identified that can be used as a haplotype to track inheritance of the high ovulation rate allele in descendants of the carrier sire.

Twinning in cattle is a complex trait that is associated with economic loss and health issues such as abortion, dystocia, and reduced calf survival. Twinning-rate QTL have been detected previously on BTA5 in the North American Holstein and Norwegian dairy cattle populations and in a USDA herd selected for high twinning rate. In previous work with the North American Holstein population, the strongest evidence for a QTL was obtained from analysis of an extended, multiple-generation family. Using additional animals, an increased density of SNP marker association tests, and a combined linkage and linkage disequilibrium mapping method, we refined the position of this QTL in the North American Holstein population. Two sets of twinning-rate predicted transmitting abilities estimated during 2 different time periods in the North American dairy cattle population were used to provide validation of results. A total of 106 SNP and 3 microsatellites were used to scan the genomic region between 5 and 80 Mb on BTA5. Combined linkage-linkage disequilibrium analysis identified significant evidence for QTL within the 25- to 35-Mb and 64- to 70-Mb regions of BTA5. The IGF-1 gene (IGF1) was examined as a positional candidate gene and an SNP in intron 2 of IGF1 was significantly associated with twinning rate by using both data sets (P = 0.003 and P = 1.05 × 10⁻⁶). Replication of this association in other cattle populations will be required to examine the extent of linkage disequilibrium with the underlying quantitative trait nucleotide across breeds.