Whole genome analysis of infectious bovine keratoconjunctivitis in Angus cattle using Bayesian threshold models
© Kizilkaya et al; licensee BioMed Central Ltd. 2011
Published: 3 June 2011
Infectious bovine keratoconjunctivitis (IBK), also known as pinkeye, is characterized by damage to the cornea and is an economically important, lowly heritable, categorical disease trait in beef cattle. Scores of eye damage were collected at weaning on 858 Angus cattle. SNP genotypes for each animal were obtained from BovineSNP50 Infinium-beadchips. Simultaneous associations of all SNP with IBK phenotype were determined using Bayes-C that treats SNP effects as random with equal variance for an assumed fraction (π=0.999) of SNP having no effect on IBK scores. Bayes-C threshold models were used to estimate SNP effects by classifying IBK into two, three or nine ordered categories. Magnitudes of genetic variances estimated in localized regions across the genome indicated that SNP within the most informative regions accounted for much of the genetic variance of IBK and pointed out some degree of association to IBK. There are many candidate genes in these regions which could include a gene or group of genes associated with bacterial disease in cattle.
Infectious Bovine Keratoconjunctivitis (IBK), commonly known as pinkeye, is a contagious bacterial disease caused in beef cattle by Moraxella bovis. IBK is characterized by excessive tearing, inflammation of the conjunctiva, and ulceration of the cornea. In severe cases, perforation of the cornea may occur which leads to permanent blindness. IBK is considered the most important ocular disease in cattle, due to the decreased performance of infected individuals and its subsequent economic effects. IBK is an economically important, lowly heritable (0.10 - 0.25), categorical disease trait. Breed differences in susceptibility of IBK have been established. In a number of studies, Hereford, Jersey and Holstein breeds appear to be more susceptible to infection than Bos Indicus breeds [1, 2]. Additive genetic effects within breed are also known to influence resistance to pinkeye . Putative QTL on Chromosomes 1 (66 to 110 cM) and 20 (2 to 35 cM) have been reported  to be associated with IBK.
The objective of the current study was to detect single-nucleotide polymorphism (SNP) markers in linkage disequilibrium (LD) with QTL associated with IBK in Angus cattle.
Material and Methods
Eye scores were pooled in various manners for categorical analysis as follows:
IBK classified in two categories: Incidence was scored as 0 for both unaffected eyes (63.7 %) and as 1 otherwise (37.3 %).
IBK classified in three categories: Incidence was scored as 0 for both unaffected eyes (63.7 %), 1 for single affected eye (26.4 %) and 2 for both affected eyes (10.9 %).
IBK classified in nine categories: Incidence was scored from 0 (63.7 %) to 8 (0.58 %) by adding the scores of the left and right eyes to utilize more information. For example, an animal with right eye score of 2 and a left eye score of 1 was scored as a 3.
50k SNP data
High-density SNP genotypes were obtained using the Bovine SNP50 Infinium II BeadChip (Illumina, Inc., San Diago CA). The Illumina A/B allele calls were used to compute a covariate for each locus that had values 0, 1, or 2 representing the number of B alleles. Missing genotypes represented less than 0.2% of total observations and were replaced with average covariate values. All genotypes were retained regardless of minor allele frequency.
where l c i is the underlying liability to IBK (classified into two, three or nine categories, c = 2, 3 or 9) of individual i; b is a vector of contemporary group effects defined as the interaction of year and pasture management group; x i is an incidence vector for b; K is the number of SNP marker loci; Z ij is the covariate at locus j for individual i; β j is the random substitution effect for locus j, which is conditional on and assumed normally distributed when δ j =1 but β j =0 when δ j =0, δ j is a random 0/1 variable indicating the absence (δ j =0) with probability π or presence (δ j =1) with probability 1 — π of locus j in the model; and e i is the random residual effect assumed normally distributed with µ=0 and known .
The Bayes-C threshold model was implemented in GenSel software  by assuming a fraction (π=0.999) of SNP markers have no effect on IBK scores. In the Markov Chain Monte Carlo (MCMC) implementation, a burn-in period of 5,000 MCMC cycles was used before saving samples from each of an additional 40,000 MCMC cycles. SNP marker effects β j were estimated by computing Monte-Carlo means of the posterior distribution of these effects using a Gibbs sampling strategy described in Kizilkaya et al., .
Results and Discussion
The regions on chromosomes 2 (4000, relative SNP position), 13 (30900) and 23 (48300) were found to be common across threshold model analyses and to have the highest genetic variances in the whole genome analysis of IBK scores for two, three or nine categories. Thereby, these regions could be indicative of QTL locations which are associated with IBK.
Comparison of genetic variances from analyses of IBK indicated that analysis of IBK scored into nine categories attributed smaller genetic variances to the regions on chromosomes 2, 13 and 23 than those with IBK defined in two and three categories. However, results from the nine category analysis additionally identified regions on chromosomes 1 and 20 which accounted for more genetic variance than these regions from analyses of IBK with two and three categories. SNP within the regions on chromosomes 1 and 20 were also found to reside within QTL regions identified by Casas and Stone  and within the QTL region on chromosome 20 related to incidence of pathogenic disease that include IBK, bovine respiratory disease, and infectious pododermatitis . Furthermore, SNP within the regions identified in the present study were in regions of QTL for somatic cell count .
Whole genome analysis of IBK using a marker-based threshold model identified genomic regions that were associated with IBK. Degree of associations determined in the analyses of IBK with two, three or nine classification were different for regions on chromosome 2, 13 and 23, which resulted from varying power across threshold model analyses depending on different number of observations in IBK categories with two, three or nine classification. Analysis of IBK with nine categories estimated smaller genetic variances for certain regions; however, it seemed to compensate for the loss of information by providing more indicative regions of QTL across the genome, especially on chromosomes 1 and 20. Analyses herein also provided confirmatory evidence of SNP associated with putative QTL on chromosome 1 and 20 for IBK. These findings should motivate future studies in whole genome analysis to identify the genetic basis of IBK.
Authors would like to thank Iowa State University Center for Integrated Animal Genomics for financial support of phenotype collections.
This article has been published as part of BMC Proceedings Volume 5 Supplement 4, 2011: Proceedings of the International Symposium on Animal Genomics for Animal Health (AGAH 2010). The full contents of the supplement are available online at http://www.biomedcentral.com/1753-6561/5?issue=S4.
- Frisch JE: The relative incidence and effect of bovine infectious keratoconjunctivitis in bos indicus and bos taurus cattle. Anim Prod. 1975, 21: 265-274. 10.1017/S0003356100030737.View ArticleGoogle Scholar
- Snowder GD, Van Vleck LD, Cundiff LV, Bennett GL: Genetic and environmental factors associated with incidence of infectious bovine keratoconjunctivitis in preweaned beef calves. J Anim Sci. 2005, 83: 507-518.PubMedGoogle Scholar
- Casas E, Stone RT: Putative quantitative trait loci associated with the probability of contracting infectious bovine keratoconjunctivitis. J Anim Sci. 2006, 84: 3180-3184. 10.2527/jas.2006-200.View ArticlePubMedGoogle Scholar
- Rodríguez JE: Infectious Bovine Keratoconjunctivitis in Angus cattle. Master’s thesis. 2006, Iowa State University, Genetics DepartmentGoogle Scholar
- Kizilkaya K, Fernando RL, Garrick DJ: Genomic prediction of simulated multibreed and purebred performance using observed fifty thousand single nucleotide polymorphism genotypes. J Anim Sci. 2010, 88: 544-551. 10.2527/jas.2009-2064.View ArticlePubMedGoogle Scholar
- Fernando RL, Garrick DJ: GenSel-User manual for a portfolio of genomic selection related analyses, Ver. 2.14. 2009, [http://taurus.ansci.iastate.edu/gensel]Google Scholar
- Casas E, Snowder GD: Putative quantitative trait loci associated with bovine pathogenic disease incidence. J Anim Sci. 2008, 86: 2455-2460. 10.2527/jas.2008-0933.View ArticlePubMedGoogle Scholar
- Rodriguez-Zas SL, Southey BR, Heyen DW, Lewin HA: Interval and composite interval mapping of somatic cell score, yield, and components of milk in dairy cattle. J Dairy Sci. 2002, 85: 3081-3091. 10.3168/jds.S0022-0302(02)74395-6.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.