A genome-wide association study reveals a locus for bilateral iridal hypopigmentation in Holstein Friesian cattle
© The Author(s). 2017
Received: 12 January 2017
Accepted: 25 March 2017
Published: 29 March 2017
Eye pigmentation abnormalities in cattle are often related to albinism, Chediak-Higashi or Tietz like syndrome. However, mutations only affecting pigmentation of coat color and eye have also been described. Herein 18 Holstein Friesian cattle affected by bicolored and hypopigmented irises have been investigated.
Affected animals did not reveal any ophthalmological or neurological abnormalities besides the specific iris color differences. Coat color of affected cattle did not differ from controls. Histological examination revealed a reduction of melanin pigment in the iridal anterior border layer and stroma in cases as cause of iris hypopigmentation. To analyze the genetics of the iris pigmentation differences, a genome-wide association study was performed using Illumina BovineSNP50 BeadChip genotypes of the 18 cases and 172 randomly chosen control animals. A significant association on bovine chromosome 8 (BTA8) was identified at position 60,990,733 with a -log10(p) = 9.17. Analysis of genotypic and allelic dependences between cases of iridal hypopigmentation and an additional set of 316 randomly selected Holstein Friesian cattle controls showed that allele A at position 60,990,733 on BTA8 (P = 4.0e–08, odds ratio = 6.3, 95% confidence interval 3.02–13.17) significantly increased the chance of iridal hypopigmentation.
The clinical appearance of the iridal hypopigmentation differed from previously reported cases of pigmentation abnormalities in syndromes like Chediak-Higashi or Tietz and seems to be mainly of cosmetic character. Iridal hypopigmentation is caused by a reduced content of melanin pigment in the anterior border layer and iridal stroma. A single genomic position on BTA8 was detected to be significantly associated with iridal hypopigmentation in examined cattle. To our knowledge this is the first report about this phenotype in Holstein Friesian cattle.
KeywordsIris hypopigmentation Cattle Heterochromia iridis Oculocutaneous hypopigmentation Albinism GWAS Holstein Friesian
For many decades, eye color and eye color genetics has been an intensively studied field of research in humans and other species. More than 100 years ago, a relatively plain dominant-recessive model of inheritance of human brown and blue eye color has been assumed . However, in recent years it became clear that iris pigmentation is under the control of a plethora of genes and is rather a quantitative than a simple Mendelian trait [2, 3].
Pigmentation of the eye, skin, and hair is the result of melanin pigment production in melanocytes. Melanin producing cells contain specialized lysosome-related organelles, the melanosomes, depositing melanin pigment in mature stage. Brighter colored eyes usually have less melanin pigment in iris stroma than darker eyes [2, 4, 5]. This and other factors as pigmentation of the posterior pigment epithelium , the type of melanin (eu- or pheomelanin) in iris melanocytes , light-scattering, and absorption processes  are believed to influence human eye color determination.
Changes in iris color have also been reported in cattle. Discolorations of the iris, either mono- or bilateral, complete or partial, are usually referred to as heterochromia iridis (HI). However, the phenotypic appearance differs remarkably between reported cases. The most prominent iris color variations were detected in cattle suffering from complete albinism, showing a pale blue iris with a white periphery [6–8]. Less distinct pigmentation anomalies were observed in non-albinotic HI cases, showing a bicolored iris with a central ring of blue and a peripheral ring of gray or brown [9, 10]. Besides albinism, severe eye color changes in cattle were also observed in syndromes like Tietz  and Chediak-Higashi [12, 13]. These syndrome related pigmentation alterations are usually accompanied by more restrictive anomalies. German Fleckvieh cattle with Tietz like syndrome exhibited bilateral deafness and colobomatous eyes . Chediak-Higashi syndrome usually manifested in bruisability and bleeding tendency [12, 13]. In albino cases with HI further clinical features related to eye development like nystagmus and blindness were detected .
Recently an alteration of iris coloration has been observed in Angus and Simmental breed. Affected cattle showed an oculocutaneous hypopigmentation (OH) with a pale blue iris and a tan periphery coupled with a change in coat color from black to chocolate. It is assumed that this aberration was introduced into the Simmental breed in the late 1950s by Angus founders and is inherited as an autosomal recessive trait. An amino acid exchange in the Ras-related Protein Rab-38 (RAB38) gene was identified as the disease causing mutation (Jon Beever, personal communication). Likewise, Rab38 cht /Rab38 cht mice with a mutation (G146T) in exon 1 of Rab38 develop a similar phenotype with chocolate coat color and ocular hypopigmentation [14, 15].
In the current study the clinical, histological and molecular examination of bilateral iris hypopigmentation in 18 Holstein Friesian cattle are described. The aim of the study was the description of the phenotype, the identification of pigmentation alterations in the eye of affected cattle and the determination of the underlying genetics of iridal hypopigmentation in cattle.
Animals, pedigree information and DNA samples
A total of 18 Holstein Friesian (HF) cattle (9 male, 9 female) with hypopigmented irises originating from eight different farms were used for this study. Pedigree data were obtained from the German livestock database service provider (VIT) and checked for shared common ancestors. Complete pedigree data were available for ten individuals, while for eight animals only paternal pedigree data were present. Pedigrees of HF cases were constructed using Pedigraph .
DNA was extracted from EDTA-blood using MagNa Pure LC DNA Isolation Kit I (Roche Diagnostics Deutschland GmbH, Mannheim, Germany). Control DNA was obtained from the depository at the Institute of Veterinary Medicine (Göttingen, Germany).
To determine the exact cause of the hypopigmentation, histological evaluation of irises of unaffected and affected animals was conducted. Eyes of each animal were completely enucleated immediately after slaughtering. The dorsal half of each freshly enucleated eyeball was opened by incision of the sclera midway between the cornea and the optic nerve. After removal of the vitreous body, eyeballs were immersed in 4% phosphate-buffered formaldehyde. After fixation for at least 48 h, lenses were removed and four cross sections were made through the anterior half of the eyeballs including dorsal, medial, ventral, and lateral aspects of the iris and adjacent structures of the anterior eye including ciliary body and cornea. Additionally, a cross section was prepared from the lens and from the caudal half of the eyeball at the level of the optic nerve. Trimmed tissue samples were paraffin-embedded, sectioned at 3 μm, and stained with hematoxylin and eosin (HE) for light microscopic examination.
Genotyping, genome-wide association study (GWAS) and statistical analysis
For GWAS 172 randomly selected HF cattle genotypes were used for GWAS. The 18 cases of iridal hypopigmentation were also genotyped using the Illumina BovineSNP50 BeadChip. Final reports were generated using GenomeStudio V2011.1 (Illumina, San Diego, USA) and imported into SNP & Variation Suite (SVS) 8.5.0 (Golden Helix, Bozeman, USA). Genotype data were screened through a series of quality control criteria, including Mendelian errors, minor allele frequency (MAF) < 1%, p-value of Fisher’s Hardy-Weinberg equilibrium (HWE) test < 0.001 (based on controls), and single nucleotide polymorphism (SNP) call rate < 98% reducing the data set from 54,610 to 44,952 SNPs. Associations were calculated under an additive, recessive, and dominant model . The additive genetic model fitted the data best. We did not detect significant evidance for population stratification. Genomic positions refer to NCBI UMD3.1.1.
Genotypic and allelic dependences were calculated between the 18 cases of iridal hypopigmentation and an additional set of randomly chosen 316 controls. Genotypes of this validation cohort were extracted from another data set generated using the Illumina BovineSNP50 BeadChip. Genotypes were compared using 3×2 or 2×2 contingency tables and Fisher’s exact or χ 2 statistics (df = 2). P < 0.005 was considered to be significant. Calculations were done using Microsoft Excel for Mac 2011 (14.7.1). HWE χ 2-values were calculated according to Rodriguez et al.  and considered significant with p > 0.05 (df = 2).
A haplotype association test was performed using a moving window of five markers. The expectation-maximization algorithm was applied using 50 iterations and a convergence tolerance of 0.0001 . Results were corrected for multiple testing according to Bonferroni.
Polymerase chain reaction (PCR) and sanger sequencing of RAB38
RAB38 primers were designed using the online software tool Primer-Blast  with the following sequences RAB38_Ex1_F: 5’-CTTCCCGGGTCCGCAG-3’, RAB38_Ex1_R: 5’-CTGGCACAGGAGATGGTCTG-3’, RAB38_Ex2_F: 5’-ACTTTGCGGAGTGATCTGCT-3’, RAB38_Ex2_R: 5’-GCTGCCTTAGCCACAAACAC-3’, RAB38_Ex3_F: 5’-CATGGGACAGGGTTTAGAAAGAGA-3’, RAB38_Ex3_R: 5’-AGGCATAGGTCTTTGGCTTG-3’. Primers were designed to amplify the complete coding region including splice sites of RAB38. All primers were synthesized by Sigma-Aldrich (Steinheim, Germany). PCR was performed in a total volume of 25 μl using FastStart Taq DNA Polymerase, dNTPack (Roche Diagnostics, Mannheim, Germany). One reaction mix (25 μl) included 1.5 U Faststart Taq DNA Polymerase, 200 μmol/L dNTP, 0.4 μmol/L of each primer, 1× PCR reaction buffer (including 20 mM MgCl2), 1× Q-Solution (Qiagen, Hilden, Germany) and 40 ng of DNA. Cycling conditions were 95 °C for 10 min, followed by 30 cycles of 95 °C for 30 s, 60 °C for 30 s and 72 °C for 30 s. Final elongation step was 72 °C for 5 min. PCR products were sequenced with the BigDye Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems, Fisher Scientific GmbH, Schwerte, Germany) on an ABI PRISM 3130xl Genetic Analyzer (Life Technologies, Foster City, USA) after purification using Rapid PCR Cleanup Enzyme Set (New England Biolabs GmbH, Frankfurt am Main, Germany). As reference sequence of RAB38 assembly UMD3.1.1 with accession number AC_000186.1 was used.
Comprehensive clinical examination of affected cattle
Detailed ophthalmological examination of the irises of affected animals showed two merging, but clearly differentiable shades of color, a bluish center and a grayish peripheral ring. Although all animals showed a clear discoloration of the iris, there was considerable variation of iris coloration between animals. The color of the central iridal parts ranged from silvery-blue to gray-blue with darker and lighter parts. In the periphery, irises were light brown to gray with occasional light gray zones (Fig. 1). In some animals, the degree of discoloration differed within the peripheral iris and showed alternating darker and lighter regions. Partial brownish corneoscleral pigmentation was visible in a few animals.
Determination of common ancestors
To identify a potential founder of the eye color phenotype, a pedigree analysis was performed. Available pedigree data of HF cattle revealed that the 18 cases descended from 11 different sires whereas dams were not closely related. None of the dams had been reported to show the typical iris hypopigmentation. In total 10 male ancestors, born between 1954 and 1983, were identified being present in the dam and sire line of every case. The available pedigree data did not reveal a single common founder for HF cases, and HF cases were not closely related.
Determination of associated chromosomal regions
Genotypic dependences of iridal hypopigmentation
Cases (n = 18)
Controls (n = 316)
Allelic dependences of iridal hypopigmentation
Cases (n = 18)
Controls (n = 316)
Haplotype association analysis for markers flanking BTB-00352779
First Marker a
0.95 CI e
Exclusion of RAB38 as candidate for bilateral iridal hypopigmention
Although the analysed cases of iridal hypopigmentation did not show the typical phenotype described for oculocutaneous hypopigmentation (OH) in Angus and Simmental cattle, i.e. eyes with pale blue irises around the pupil, a tan periphery and slightly bleached black coats looking grey or red, bovine RAB38 was screened for mutations in the 20 cases and controls. In addition, RAB38 is located on BTA29 which did not show any associated SNP in the GWAS. However, to exclude RAB38 as causative for iridal hypopigmentation coding regions including splice sites were amplified and fragments sequenced. As expected no disease causing polymorphisms or alterations from the reference genome for cases and controls were detected.
Iridal hypopigmentation in cattle is attributed to a reduced content of melanin pigment in the anterior border layer and iridal stroma. These findings clearly differ from previously reported histological observations in HI cases, showing a reduction in eye pigmentation in different uveal structures as different iris layers, retinal pigment epithelium (RPE) , and choroid [6, 25]. Leipold and Huston observed a non-albinotic case of HI with reduced iris pigmentation in anterior border layer, iris stroma, and posterior pigment epithelium in Hereford cattle. Pigmentation in the remaining pigmented eye structures was also reduced and iridal stroma was hypoplastic . Iris hypopigmentation in cases was exclusively located in the anterior border layer and the iridal stroma, and no structures were fully devoid of pigment. Pigmentation differences as seen in the present cases of iris hypopigmentation seemed to be rather comparable with naturally occuring eye color variances in humans than to HI. Human eye color variance is mainly due to differences in the amount of melanin pigment in the anterior border layer and iridal stroma. Pale colored eyes generally contain less melanin compared to brown eyes [2, 5].
In cattle with iris hypopigmentation, gradual iridal brightening can be explained by a reduced melanin content of variable intensity affecting the anterior border layer and the iridal stroma, resulting in blue to gray-brown instead of the normal black color. Differences in coloration between central and peripheral parts of the iris, namely the rather bluish color of the central iridal ring and the gray to light brown color of peripheral iridal parts, are probably attributed to the increasing thickness of the iris towards the periphery, associated with increased amounts of collagen in the peripheral iridal stroma. As the amount of collagen fibres within the iridal stroma is known to be another important determining factor of eye color, at least in humans , it may be assumed that variations in iris brightening between central and peripheral iridal parts in cattle with iris hypopigmentation result from differences in the amount of collagen and, thus, in the backscattering properties of the respective iridal regions. In this respect it is noteworthy that in the associated region on BTA8 collagen gene COL15A1 is located at position 64,437,276–64,540,739 which has been shown to be expressed in multiple ocular structures . In contrast, the amount of collagen in the iridal stroma of normal dark brown to black-eyed cattle seems to have a minor influence on iris coloration as almost all light is probably absorbed by the extensive eumelanin deposits in the anterior border layer and the iridal stroma accounting for a homogenous dark eye color [27, 28]. Variations of the degree of discoloration within the peripheral iris with alternating darker and lighter regions in affected animals of the present study are caused by regional differences in the melanin content of the anterior border layer and the iridal stroma, which were confirmed by histological analyses.
Taking all clinical findings together, there were no other abnormalities found in conjunction with the iris hypopigmentation. Although no long-term effects were examined, iris hypopigmentation seems to be mainly of cosmetic character. Under physiological conditions melanin pigment protects from ultraviolet (UV) light, and humans with lighter eye color seem to be more susceptible to age related macula degeneration  and uveal melanoma [30, 31]. This is comparable with blue-eyed horses that are at higher risk to develop ocular squamous cell carcinoma .
The bilateral iridal hypopigmentation in HF cattle described here was due to a reduction of melanin pigment in the anterior border layer and iridal stroma. The phenotype was highly associated with a chromosomal region on BTA8. Haplotype association analysis showed that the presence of the A-alleles at the associated SNPs significantly increased the chance of developing an iridal hypopigmentation.
Bos taurus chromosome
Degrees of freedom
Genome-wide association study
Holstein Friesian cattle
Minor allele frequency
Single nucleotide polymorphism
The authors thank M. Scharfenstein, S. Loos, L. Jüttner, and S. Pach for their excellent technical assistance. The German livestock database service provider (VIT) is thanked for providing pedigree data. The cattle breeders are thanked for providing samples and pictures of cases. C. Eule is thanked for critical evaluation of clinical signs. We are grateful to H. B. Urnovitz for critical review of the manuscript. We acknowledge support by the German Research Foundation and the Open Access Publication Funds of the Göttingen University.
Availability of data and materials
All relevant information supporting the results of this article are included within the article and its additional files. The raw datasets analysed during the current study are available from the corresponding author on reasonable request.
BB and ES were responsible for the study design. AKH performed the genotyping. MB performed the histological analysis. The clinical examinations were done by AT, JNN and WEW. AKH and WEW performed the pedigree analysis. BB performed the GWAS and the statistical analysis. AKH and BB wrote the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
The current study was performed according to the ethical guidelines of the University of Goettingen and the German Animal Welfare act. Genotyping data and blood sampling are based on routine diagnostic parentage control and genomic selection programs performed with written owner consent. Collection of blood samples was conducted exclusively by local veterinarians. Blood sampling by veterinarians with state examination is in accordance with the German Animal Welfare Act (§6 Abs. 1 Satz 2 TierSchG). Therefore no formal ethical approval was required, since no other samples were collected for this study. Animals were regularly slaughtered at an abattoir according to §4 of the German Animal Welfare Act.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Davenport GC, Davenport CB. Heredity of eye-color in man. Science. 1907;26(670):589–92.View ArticlePubMedGoogle Scholar
- Sturm RA, Larsson M. Genetics of human iris colour and patterns. Pigment Cell Melanoma Res. 2009;22(5):544–62.View ArticlePubMedGoogle Scholar
- Edwards M, Cha D, Krithika S, Johnson M, Cook G, Parra EJ. Iris pigmentation as a quantitative trait: variation in populations of European, East Asian and South Asian ancestry and association with candidate gene polymorphisms. Pigment Cell Melanoma Res. 2016;29(2):141–62.View ArticlePubMedGoogle Scholar
- Zhu G, Evans DM, Duffy DL, Montgomery GW, Medland SE, Gillespie NA, Ewen KR, Jewell M, Liew YW, Hayward NK, et al. A genome scan for eye color in 502 twin families: most variation is due to a QTL on chromosome 15q. Twin Res. 2004;7(2):197–210.View ArticlePubMedGoogle Scholar
- Rennie IG. Don’t it make my blue eyes brown: heterochromia and other abnormalities of the iris. Eye (Lond). 2012;26(1):29–50.View ArticleGoogle Scholar
- Leipold HW, Huston K. A Herd of glass-eyed albino Hereford cattle. J Hered. 1966;57(5):179–82.View ArticlePubMedGoogle Scholar
- Leipold HW, Huston K, Gelatt KN. Complete Albinism in a Guernsey Calf. J Hered. 1968;59(3):218–20.View ArticlePubMedGoogle Scholar
- Jayasekera U, Leipold HW. Albinism in United-States Charolais Cattle. Ann Genet Sel Anim. 1981;13(3):213–8.PubMedPubMed CentralGoogle Scholar
- Leipold HW, Huston K. Incomplete albinism and heterochromia irides in Herefords. J Hered. 1968;59(1):3–8.View ArticlePubMedGoogle Scholar
- Huston K, Leipold HW, Freeman AE. Heterochromia Irides in Dairy Cattle. J Dairy Sci. 1968;51(7):1101–2.View ArticlePubMedGoogle Scholar
- Philipp U, Lupp B, Momke S, Stein V, Tipold A, Eule JC, Rehage J, Distl O. A MITF Mutation Associated with a Dominant White Phenotype and Bilateral Deafness in German Fleckvieh Cattle. Plos One. 2011;6(12):e28857.View ArticlePubMedPubMed CentralGoogle Scholar
- Padgett GA, Leader RW, Gorham JR, O’Mary CC. The Familial Occurrence of the Chediak-Higashi Syndrome in Mink and Cattle. Genetics. 1964;49:505–12.PubMedPubMed CentralGoogle Scholar
- Shiraishi M, Ogawa H, Ikeda M, Kawashima S, Ito K. Platelet dysfunction in Chediak-Higashi syndrome-affected cattle. J Vet Med Sci. 2002;64(9):751–60.View ArticlePubMedGoogle Scholar
- Loftus SK, Larson DM, Baxter LL, Antonellis A, Chen Y, Wu X, Jiang Y, Bittner M, Hammer 3rd JA, Pavan WJ. Mutation of melanosome protein RAB38 in chocolate mice. Proc Natl Acad Sci U S A. 2002;99(7):4471–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Brooks BP, Larson DM, Chan CC, Kjellstrom S, Smith RS, Crawford MA, Lamoreux L, Huizing M, Hess R, Jiao X, et al. Analysis of ocular hypopigmentation in Rab38cht/cht mice. Invest Ophthalmol Vis Sci. 2007;48(9):3905–13.View ArticlePubMedPubMed CentralGoogle Scholar
- Garbe JR, Da Y. Pedigraph: a software tool for the graphing and analysis of large complex pedigree. User manual Version 24 Department of Animal Science, University of Minnesota 2008. Available: https://animalgene.umn.edu/pedigraph. Accessed 28 Mar 2017.
- Bae H, Perls T, Steinberg M, Sebastiani P. Bayesian polynomial regression models to fit multiple genetic models for quantitative traits. Bayesian Anal. 2015;10(1):53–74.View ArticlePubMedGoogle Scholar
- Rodriguez S, Gaunt TR, Day IN. Hardy-Weinberg equilibrium testing of biological ascertainment for Mendelian randomization studies. Am J Epidemiol. 2009;169(4):505–14.View ArticlePubMedPubMed CentralGoogle Scholar
- Excoffier L, Slatkin M. Maximum-likelihood estimation of molecular haplotype frequencies in a diploid population. Mol Biol Evol. 1995;12(5):921–7.PubMedGoogle Scholar
- Ye J, Coulouris G, Zaretskaya I, Cutcutache I, Rozen S, Madden TL. Primer-BLAST: a tool to design target-specific primers for polymerase chain reaction. BMC Bioinformatics. 2012;13(1):1–11.View ArticleGoogle Scholar
- Mayhew IGJ. Large animal neurology. Oxford: Wiley-Blackwell; 2009.Google Scholar
- Renard D, Jeanjean L, Labauge P. Heterochromia Iridis in congenital Horner’s syndrome. Eur Neurol. 2010;63(4):253.View ArticlePubMedGoogle Scholar
- Diesenhouse MC, Palay DA, Newman NJ, To K, Albert DM. Acquired heterochromia with horner syndrome in two adults. Ophthalmology. 1992;99(12):1815–7.View ArticlePubMedGoogle Scholar
- Nie C, Zhang Z, Zheng J, Sun H, Ning Z, Xu G, Yang N, Qu L. Genome-wide association study revealed genomic regions related to white/red earlobe color trait in the Rhode Island red chickens. BMC Genet. 2016;17(1):115.View ArticlePubMedPubMed CentralGoogle Scholar
- Leipold HW, Huston K. Histopathology of incomplete albinism and heterochromia irides in the Hereford. Cornell Vet. 1969;59(1):69–75.PubMedGoogle Scholar
- Wiggs JL, Howell GR, Linkroum K, Abdrabou W, Hodges E, Braine CE, Pasquale LR, Hannon GJ, Haines JL, John SW. Variations in COL15A1 and COL18A1 influence age of onset of primary open angle glaucoma. Clin Genet. 2013;84(2):167–74.View ArticlePubMedPubMed CentralGoogle Scholar
- Wilkerson CL, Syed NA, Fisher MR, Robinson NL, Wallow IL, Albert DM. Melanocytes and iris color: light microscopic findings. Arch Ophthalmol. 1996;114(4):437–42.View ArticlePubMedGoogle Scholar
- Liu Y, Hong L, Wakamatsu K, Ito S, Adhyaru BB, Cheng CY, Bowers CR, Simon JD. Comparisons of the structural and chemical properties of melanosomes isolated from retinal pigment epithelium, iris and choroid of newborn and mature bovine eyes. Photochem Photobiol. 2005;81(3):510–6.View ArticlePubMedGoogle Scholar
- Frank RN, Puklin JE, Stock C, Canter LA. Race, iris color, and age-related macular degeneration. Trans Am Ophthalmol Soc. 2000;98:109–17.PubMedPubMed CentralGoogle Scholar
- Gallagher RP, Elwood JM, Rootman J, Spinelli JJ, Hill GB, Threlfall WJ, Birdsell JM. Risk factors for ocular melanoma: Western Canada Melanoma study. J Natl Cancer Inst. 1985;74(4):775–8.PubMedGoogle Scholar
- Schmidt-Pokrzywniak A, Jockel KH, Bornfeld N, Sauerwein W, Stang A. Positive interaction between light iris color and ultraviolet radiation in relation to the risk of uveal melanoma: a case-control study. Ophthalmology. 2009;116(2):340–8.View ArticlePubMedGoogle Scholar
- Bergstrom BE, Labelle AL, Pryde ME, Hamor RE, Myrna KE. Prevalence of ophthalmic disease in blue-eyed horses. Equine Vet Educ. 2014;26(8):438–40.View ArticleGoogle Scholar