Genome-wide association study of infectious bovine keratoconjunctivitis in Angus cattle
© Kizilkaya et al.; licensee BioMed Central Ltd. 2013
Received: 18 October 2012
Accepted: 22 February 2013
Published: 26 March 2013
Skip to main content
© Kizilkaya et al.; licensee BioMed Central Ltd. 2013
Received: 18 October 2012
Accepted: 22 February 2013
Published: 26 March 2013
Infectious Bovine Keratoconjunctivitis (IBK) in beef cattle, commonly known as pinkeye, is a bacterial disease caused by M o r a x e l l a bovis. IBK is characterized by excessive tearing and ulceration of the cornea. Perforation of the cornea may also occur in severe cases. IBK is considered the most important ocular disease in cattle production, due to the decreased growth performance of infected individuals and its subsequent economic effects. IBK is an economically important, lowly heritable categorical disease trait. Mass selection of unaffected animals has not been successful at reducing disease incidence. Genome-wide studies can determine chromosomal regions associated with IBK susceptibility. The objective of the study was to detect single-nucleotide polymorphism (SNP) markers in linkage disequilibrium (LD) with genetic variants associated with IBK in American Angus cattle.
The proportion of phenotypic variance explained by markers was 0.06 in the whole genome analysis of IBK incidence classified as two, three or nine categories. Whole-genome analysis using any categorisation of (two, three or nine) IBK scores showed that locations on chromosomes 2, 12, 13 and 21 were associated with IBK disease. The genomic locations on chromosomes 13 and 21 overlap with QTLs associated with Bovine spongiform encephalopathy, clinical mastitis or somatic cell count.
Results of these genome-wide analyses indicated that if the underlying genetic factors confer not only IBK susceptibility but also IBK severity, treating IBK phenotypes as a two-categorical trait can cause information loss in the genome-wide analysis. These results help our overall understanding of the genetics of IBK and have the potential to provide information for future use in breeding schemes.
Infectious Bovine Keratoconjunctivitis (IBK), commonly known as pinkeye, is a highly contagious ocular bacterial cattle disease which occurs in cattle populations throughout the world. IBK is caused by the Gram negative bacterium Moraxella bovis and is characterized by excessive tearing, inflammation of the conjunctiva, and ulceration of the cornea in one or both eyes. As the disease progresses the cornea becomes cloudy or white. In severe cases, perforation of the cornea may occur which can lead to permanent blindness.
IBK is non-fatal; however, it is considered the most important ocular disease in cattle, due to the decreased growth performance of infected individuals. The United State National Animal Health Monitoring System survey  and Australian postal survey  reported IBK as an economically important disease in animal production. Major economic losses are the result of inappetence and poor weight gain in affected animals suffering from ocular pain and visual impairment . It has been estimated that IBK costs beef producers 150 million US$ in the United States  and 22 million AUD$ in Australia  per annum.
There has been little research to quantify genetic variation in susceptibility to IBK. Breed differences in susceptibility to IBK have been demonstrated and Hereford cattle were found to be more susceptible compared with all other purebreds (Angus, Braunvieh, Charolais, Gelbvieh, Limousin, Pinzgauer, Red Poll and Simmental) and Bos Indicus breeds [5, 6]. Heritability estimates (0.10 - 0.25) different from zero have been reported, highlighting the presence of within-breed genetic variation in susceptibility to IBK.
Putative QTL on Chromosomes 1 (66 to 110 cM) and 20 (2 to 35 cM) have been reported  to be associated with resistance to IBK. Reducing IBK through selection would be advantageous, because genetic gain is cumulative and permanent. Once the loci associated with susceptibility or resistance to IBK are identified, relevant markers can be used for selection in breeding schemes.
In recent years, high-density single nucleotide polymorphism (SNP) genotyping assays for cattle have become available . In addition, statistical methodologies have been developed to make use of these SNP panels to detect association between SNP and economically important traits. The use of genome-wide studies helps identify chromosomal regions associated with disease incidence in immunity-related diseases such as IBK. The objective of this study was to detect SNP markers in linkage disequilibrium (LD) with genetic variants associated with IBK in Angus cattle.
All animal procedures were approved by the Iowa State University Animal Care and Use Committee.
Number of animals diagnosed with different severities of Infectious Bovine Keratoconjunctivitis for left, right or both eyes
Number of animals diagnosed with different severities of Infectious Bovine Keratoconjunctivitis classified within two, three or nine categories
High-density (53,367) SNP genotypes of purebred American Angus cattle 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 to represent 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 in the analysis regardless of minor allele frequency.
The threshold model for ordered categorical data assumes the existence of a set of ordered thresholds τ 0=−∞<τ 1<τ 2<…<τ c−1<τ c =∞ for the trait and an underlying or latent variable l i for each animal i. When the latent variable for animal i is between thresholds τ j−1 and τ j the categorical phenotype y i for animal i takes value j[9, 10].
where β is a px1 vector of fixed effects, is a known incidence row vector corresponding to fixed effects in β, K is the number of marker loci, z i k is the value of the covariate for marker k in individual i, α k is the random substitution effect for locus k, with with probability 1- π or α k =0 with probability π. For a locus included in the model, the locus-specific variance, , is assumed to have a scaled inverse chi-square distribution with degrees of freedom ν α =4 and scale parameter , where p k is the B allele frequency of marker k and is the additive genetic variance explained by markers [13, 14]. For a locus not in the model, is set to zero. A flat prior was assumed for the fixed effects and thresholds, and the residual, e i , was assumed to be distributed .
where and l are the vectors of marker effects, marker variances, thresholds and liabilities.
where I(.) is an indicator function taking the value 1 when expression (.) is true and 0 otherwise.
A Gibbs sampler was used to construct an irreducible, Markov chain with the joint posterior as its stationary distribution [10, 16, 17]. Inferences were made from the Markov chain. In the Gibbs sampler, samples are obtained from the full conditional posteriors rather than the joint posterior [10, 16, 17]. The full conditional posteriors used in this paper are described below.
X −m is the matrix X with the column associated with m deleted, and β −m is β with the mth element deleted.
At locus k, following Meuwissen et al, the locus-specific marker variance, , and the SNP effect, α k , were sampled from their joint full conditional posterior distribution. The strategy they used was to sample from the marginal of the full conditional posterior and then sample α k from the conditional posterior given the sampled value for .
where and .
as described by Devroye .
The arithmetic means of these samples were used to estimate their posterior means.
Samples of window variances were obtained using (15) but with g (t) computed with only the sampled SNP effects for window w. These samples were used to compute posterior means and posterior probabilities for percent window variances and were used for inferences on locations of causal variants.
SNP effects from genome-wide analysis were obtained using the threshold model in GenSel software  and the estimated effects were uploaded into SNPLOTz (http://www.animalgenome.org/tools/snplotz). SNPLOTz visualized the estimated SNP effects and their genomic location and was dynamically linked to Gbrowse (http://www.animalgenome.org/gbrowse), which allowed visualization of SNP with other types of genome features (i.e., annotated genes, curated QTL, transcripts, etc.) . The SNP positions within a chromosome were based on Bos taurus genome assembly (UMD 3.1).
The distribution of IBK records by left and right eyes are presented in Table 1. The rates for severity of IBK that was scored into five categories from cornea with no apparent lesions to perforation of the cornea were 75, 12, 7, 3 and 3% for left eye and 77, 11, 6, 3 and 3% for right eye. These results indicated similarity between left and right eyes for IBK incidence rates. Rodriguez  studied the effects of IBK incidence severity on production traits in Angus breed and found similar incidence rates between left and right eyes.
The distribution of IBK scores classified in two, three or nine categories of infection can be seen in Table 2. The rates of unaffected, affected, only one-eye affected and two-eye affected animals were 63, 37, 26 and 11%, respectively. The severity score of IBK within the nine category classification system varied from 12.8% to 0.58% for infected animals. These results were found to be similar with those from Rodriguez ; however, higher than those (3.7%) from Snowder et al., that studied environmental effects and genetic factors influencing the incidence of IBK among nine breeds including Angus. A significant between-breed difference in IBK incidence has been reported , with purebred Herefords (22.4%) being more susceptible compared to other pure (Simmental 7.6%, Charolais 6.5%, Angus 3.7%, Limousin 3.4%) and composite breeds. Thus, sufficient genetic variation for resistance to IBK exists in these breeds to consider selection.
The posterior distribution of heritabilities seemed similar across the different classification of IBK scores. Posterior means of heritabilities for two- (0.064) and nine-category (0.064) IBK scores were the same; however, they were slightly lower than 0.066 for three-category IBK scores. The 95% PPI of heritabilities overlapped, which indicates that there is no significant difference among heritability estimates of IBK incidence based on different classifications.
There are very few studies reporting estimates of heritability for IBK incidence in the literature. Our heritability estimates were lower, but in good agreement with the estimates (0.06 - 0.10) obtained from Rodriguez . However, our heritability estimates are lower than the heritability estimates of 0.10-0.25 obtained by Snowder et al. based on the models including only animal effect, and animal and maternal effects for Angus cattle. In the same study, Snowder et al. also found high heritability estimates (0.20-0.28) for Hereford and low heritability (0.0 - 0.13) estimates for Red Poll, Charolais, Simmental, Limousin, Gelbvieh, Pinzgauer, Braunvieh breeds. Ali et al. indicated that estimates of heritability for IBK were small to moderate (0.17 - 0.19) for both pre-weaning and post-weaning calves in a Hereford and Shorthorn composite population. These differences between heritability estimates could be attributed partly by the differences in phenotypes measured, failure to accurately distinguish between phenotypes of healthy and sick animals; false assumption that disease observed is the primary infection, lack of knowledge of the influence of passive immunity on disease incidence, susceptibility of disease biased by time, age, or season dependency and/or differences in the models of analyses used in these studies ,.
Two-category classification of IBK was performed by considering all animals whose IBK status was at least as severe as the threshold level to be ‘affected’, while considering all animals whose IBK status was less severe than the threshold level to be ‘unaffected’. It was also assumed that animals with IBK had the same levels of exposure to the genetic risk factors. In the two-category classification, however, some animals considered ‘affected’ could be in more severe states than other ‘affected’ animals . This additional information could be used in the multiple-category classification of IBK, whereas it was lost upon two-category classification. In the linkage analyses of two-, three- and multiple-categorical traits, Corbet el al. indicated that two-category classification of multiple-categorical traits could lead to crippling power loss, especially in the case of many loci of small effects. In particular, if the underlying genetic factors for any trait confer not only IBK susceptibility but also IBK severity, treating IBK phenotypes as multi-categorical could provide additional information in genome-wide analysis.
Figure 2 and Tables in the Additional file 1 indicated many genomic regions with different degree of association with IBK incidence in American Angus breed. Casas and Stone  had previously reported presence of QTL for IBK tolerance on chromosome 1 and 20. Casas and Snowder  further indicated that the region on chromosome 20 might be associated with general resistance to bacterial diseases (including IBK). In this present study, several SNP windows associated with IBK disease status were identified on chromosome 1 (Additional file 1). However, three SNP windows ([109.01 Mb] - [109.87 Mb], [110.04 Mb] - [110.98 Mb] and [112.04 Mb] - [113.00 Mb]) on chromosome 1 overlapped previously described QTL region for IBK susceptibility . There are two annotated genes within these 1-Mb SNP windows on chromosome 1; PTX3 (pentraxin 3) and IL12A (interleukin 12A, natural killer cell stimulatory factor 1, cytotoxic lymphocyte maturation factor 1, p35).
PTX3 has been reported to be rapidly produced and released by several cell types, in particular by mononuclear phagocytes, dendritic cells, fibroblasts and endothelial cells , and recognizes microbial moieties, opsonizes fungi and selected Gram-positive and Gram-negative bacteria and activates complement. Opsonization resulted in facilitated pathogen recognition and innate immune cell activation; moreover, opsonization by PTX3 is likely to be involved in the activation of an appropriate adaptive immune response . IL12A, an immunomodulatory cytokine secreted by antigen presenting cells, is critical for differentiation of T helper (Th)1 and Th2 lymphocytes . IL12 has been shown to augment the growth of activated T- and natural killer (NK)-cells , stimulate interferon gamma (IFN- γ) production by T-cells and NK cells, and suppress the expansion of Th2 cell clones . This information about PTX3 and IL12A genes indicates that they could be strong biological as well as positional candidates involved in response to IBK infections.
It is important to note that in the present study there were no common 1-Mb SNP windows on chromosome 20 for IBK incidence. While this result was not consistent with the results from the previous studies , about IBK incidence, the difference may be in part due to design differences between this study (a single breed of purebred Bos taurus animals) and the previous studies (a single Bos taurus x Bos indicus sire bred to multiple breeds of cows). These differences in results may indicate different metabolic pathways or different segregating QTL (within breed vs. across breed) for IBK incidence.
Common 1-Mb SNP windows from the analyses of IBK incidence across two, three and nine categories
rs number of
rs number of
Proportion of genetic variance in IBK classified in
SNP in window
Positional candidate genes (PCG) within 1-Mb SNP window (SNPW) and their putative functions
19.11 - 19.97
It encodes a protein that catalyzes the step of lipid biosynthesis
and the removal of long chain acid anion
It encodes a transcription factor which is a member of a small
family of basic Leucine (bZIP) proteins
22.09 - 22.77
It encodes a protein that plays an important role in the organization
of the actin cytoskeleton
It hydrolyzes ATP and can hydrolyze GTP with lower efficiency
regulates transcription by binding consensus GC and GT- box
elements in target genes
53.01 - 53.99
It encodes a tissue-specific basic helix-loop-helix (bHLH) protein
with a pivotal role in hemopoises
it mediates actin by association with G-Protein that activates
a phosphaidylinositol-calcium second messenger system
14.36 - 14.92
It encodes a transmembrane protein that is a member of soluble
N-ethylmeleimide-sensitive factor attachment protein receptor
19.00 - 19.98
It encodes a member of a neurotrophic tyrosine receptor kinase family
It involved in post transcriptional regulation of gene expression
It involved in post transcriptional regulation of gene expression
The second SNP window [22.09 - 22.77 Mb] harbors 3 genes; WIPF1 (WAS/WASL interacting protein family, member 1), OLA1 (Obg-like ATPase 1), SP3 (Sp3 transcription factor). WIPF1 encodes a protein that plays an important role in the organization of the actin cytoskeleton. Wickramarachchi et al. indicated that deficiency of actin regulators can result in defects limited to the immune system or sometimes to a single immune cell type and lack of an actin cytoskeletal regulator can cause immunodeficiency, autoimmunity, autoinflammatory disease, or a combination of these manifestations. In addition, Dustin and Cooper  pointed out that the actin cytoskeleton seems to play two critical roles in the activation of T cells: T cell shape development and movement, including formation of the immunological synapse and the formation of a scaffold for signaling components. OLA1 is a negative regulator of the antioxidative process. Knockdown of OLA1 in human cells elicited an increased resistance to oxidizing agents. Zhanga et al. also reported that knockdown of OLA1, a newly discovered regulatory protein of oxidative stress response, inhibits cell migration and invasion ability in breast cancer cells. Sp factors are able to stimulate transcription from proximal promoters or from distal enhancers. They can also physically interact with other transcription factors. van Loo et al. studied the impact of the absence of the widely expressed transcription factor SP3 on the developing hematopoietic system in the mouse. They showed that the absence of SP3 results in cell-autonomous differentiation defects in the erythroid and myeloid cell lineage. WIPF1, OLA1and SP3 genes are found to be associated with defects in immune system or cell lineage making them biological candidate for IBK disease.
Two genes were identified within the SNP window [53.01 - 53.99 Mb] on chromosome 12; SCEL (sciellin), EDNRB (endothelin receptor type B). Sciellin is a protein that is encoded by the SCEL gene. it is a precursor of the cornified enveloped. The cornified envelope is an insoluble protein complex formed under the plasma membrane in the uppermost layers of stratified squamous epithelium and plays a major role in the barrier properties of the stratum corneum . Champliaud et al. reported proteins interacting with sciellin and identified vitamin D-upregulated protein 1 (VDUP1), which plays multiple roles in a wide range of cellular processes such as proliferation or apoptosis. Recently, it has been reported that VDUP1 is also involved in the immune system via positive regulation of natural killer development . The endothelin receptor B (EDNRB) gene encodes a g-protein-coupled receptor-mediated endothelin, inducing development and transformation of the neural crest cell-specific lineage. Hypermethylation of the EDNRB promotor has been shown in multiple tumor types, and EDNRB has been proposed as a putative tumor suppressor gene .
One gene is located within SNP window [14.36 - 14.92 Mb] on chromosome 13; VAMP7 (vesicle-associated membrane protein 7). VAMP7 inhibits the release of lytic granules and severely impairs natural killer (NK) cell cytotoxic activity which is used to eliminate cancer and virus-infected cells. Furthermore, VAMP7 is involved in IFN secretion in NK cells, which indicates that VAMP7 is involved in many fusion processes and thus plays a general function in NK cell activity . In addition, this SNP window lies in the region where putative QTLs that affect Bovine spongiform encephalopathy at 22.997-61.84 cM resides . These findings indicate that these genes could be candidates for IBK susceptibility based on their functions.
The SNP window [19.00 - 19.98 Mb] on chromosome 21 harbors three genes that encode proteins; NTRK3 (neurotrophic tyrosine kinase, receptor, type 3), MIR1179 (microRNA mir-1179), MIR7-1 (microRNA mir-7-1). NTRK3 gene encodes a member of the NTRK family. These neurotrophins (NTs) receptors are best known for their role in the differentiation and survival of various types of neurons . In the study of latitude-driven adaptation for both schizophrenia and vitamin D related genes, Amato et al. found 9 genes including NTRK3 in common among those ones related to latitude, vitamin D and schizophrenia. MicroRNAs (miRNAs) are small 21-23 nucleotide-long noncoding RNAs involved in several biological process including development, differentiation, apoptosis, survival, motility, invasion and proliferation. Many miRNAs are implicated as proto-oncogenes or as tumor suppressors and aberrantly expressed in various cancer types. MIR7 has been characterized as a tumor suppressor in several human cancers. This SNP window also lies in the region where putative QTLs that affect clinical mastitis at 12.60-35.89 cM and somatic cell scores at 12.601-29.77 cM, have been identified by Schulman et al. and Schnabel et al., respectively. These results indicate that genetic variant(s) that reside in this chromosomal region could be associated with bacterial disease such as IBK in cattle.
The present analyses, using Bovine SNP50 Infinium II BeadChips, identified several 1-Mb SNP regions and genes within these regions that were associated with IBK. They may be of interest for breeding schemes as they could be used to identify animals susceptible to IBK.
IBK disease is a lowly heritable complex trait that is polygenic in nature, where many loci with small effects are expected and disease incidence is affected by the environment. Such traits require a large number of individuals to enhance the power of the genome-wide analysis and obtain good statistical support for the detection of the causal loci.
Authors would like to thank Iowa State University Center for Integrated Animal Genomics for financial support of phenotype collections. The authors acknowledge Scientific and Technological Research Council of Turkey (TUBITAK) for financial support. Authors would like to thank Mary Sue Mayes for assistance with phenotype collections and processing DNA samples for genotyping.
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.