Genetic parameters of infectious bovine keratoconjunctivitis and its relationship with weight and parasite infestations in Australian tropical Bos taurus cattle
© Ali et al.; licensee BioMed Central Ltd. 2012
Received: 14 January 2012
Accepted: 20 July 2012
Published: 27 July 2012
Infectious bovine keratoconjunctivitis (IBK) or ‘pinkeye’ is an economically important ocular disease that significantly impacts animal performance. Genetic parameters for IBK infection and its genetic and phenotypic correlations with cattle tick counts, number of helminth (unspecified species) eggs per gram of faeces and growth traits in Australian tropically adapted Bos taurus cattle were estimated.
Animals were clinically examined for the presence of IBK infection before and after weaning when the calves were 3 to 6 months and 15 to 18 months old, respectively and were also recorded for tick counts, helminth eggs counts as an indicator of intestinal parasites and live weights at several ages including 18 months.
Negative genetic correlations were estimated between IBK incidence and weight traits for animals in pre-weaning and post-weaning datasets. Genetic correlations among weight measurements were positive, with moderate to high values. Genetic correlations of IBK incidence with tick counts were positive for the pre-weaning and negative for the post-weaning datasets but negative with helminth eggs counts for the pre-weaning dataset and slightly positive for the post-weaning dataset. Genetic correlations between tick and helminth eggs counts were moderate and positive for both datasets. Phenotypic correlations of IBK incidence with helminth eggs per gram of faeces were moderate and positive for both datasets, but were close to zero for both datasets with tick counts.
Our results suggest that genetic selection against IBK incidence in tropical cattle is feasible and that calves genetically prone to acquire IBK infection could also be genetically prone to have a slower growth. The positive genetic correlations among weight traits and between tick and helminth eggs counts suggest that they are controlled by common genes (with pleiotropic effects). Genetic correlations between IBK incidence and tick and helminth egg counts were moderate and opposite between pre-weaning and post-weaning datasets, suggesting that the environmental and (or) maternal effects differ between these two growth phases. This preliminary study provides estimated genetic parameters for IBK incidence, which could be used to design selection and breeding programs for tropical adaptation in beef cattle.
Analysis of health records of 45 497 animals from purebred and crossbred Bos taurus and tropically adapted breeds, i.e. mainly from Brahman, Boran and Tuli cattle populations, estimated heritability for IBK incidence ranging from 0.00 to 0.28 .
A significant between-breed difference in IBK incidence has been reported [12, 13], with the purebred Hereford breed being the most susceptible compared to other pure and composite breeds. Thus, sufficient genetic variation for resistance to IBK exists in these breeds to consider genetic selection. To date, genetic parameters for resistance to IBK in Australian tropical composite breeds have not been reported. Thus, the main objectives of our study were to estimate genetic parameters and breeding values for resistance to IBK of a HS composite population at different stages of cattle growth. Results will enable northern Australian cattle producers to select animals with high genetic merit for resistance to IBK. However, before selection decisions are made, it is important to quantify possible genetic antagonisms between IBK and growth traits. For example, in dairy cattle, the increased genetic merit for milk yield is associated with a decline in average fertility and health (i.e. mastitis, claw and foot disorders) [14–16]. Therefore, another objective of our study was to estimate the genetic and phenotypic correlations between IBK and growth and parasite infestation traits by including these in a multi-trait genetic parameter estimation program. These multi-trait settings make use of economic weights more optimal in a total genetic merit index.
Genetic and phenotypic parameters for IBK, weight measurements at different ages, helminth egg count and cattle tick count were estimated in an HS composite cattle population. A detailed history of the tropically adapted Australian HS composite cattle population has been reported elsewhere . Briefly, the herd was developed by the Commonwealth Scientific and Industrial Research Organisation (CSIRO) during the 1950s at the National Cattle Breeding Station, ‘Belmont’, latitude 23.22 °S longitude 150.38 °E, near Rockhampton, Queensland, Australia, and was a nominal composition of 50% Hereford and 50% Shorthorn. Animals were born from October to December each year and weaned around March or April of the following year from 1973 to 1998. These animals were reared under similar stressful environmental conditions, which included high ambient temperatures and humidity, high IBK incidence and gastro-intestinal and external parasites .
Phenotypic data for IBK
Data on both the incidence and severity of IBK infection were collected from 1973 to 1997. Calves were observed for clinical signs of IBK at two stages of development, first, when they were 3 to 6 months old (‘pre-weaning’ time point) and second, when they were15 to 18 months old (‘post-weaning’ time point). The pre-weaning and post-weaning assessments involved 932 calves (403 males and 529 females) and 1 084 adult animals (514 males and 570 females), respectively. The pedigree of these animals was traced back to four generations and contained 3 360 animals sired by 168 sires (average of 20 progeny per sire) from 990 dams.
Animals were restrained in a head-bail located under an open-sided shelter and clinical examinations were performed by trained staff. Each eye of every animal was scored on a 1 to 6 scale, with 1 = cornea showing no clinical signs of ulceration or opacity (this score was used as a base line or reference); 2 = healing central corneal scar less than 3 mm in diameter, opaque but without active vessels; 3 = healing central cornea from 3 to 6 mm in diameter, opaque but without active vessels; 4 = healing central cornea scar greater than 6 mm diameter, opaque but without active vessels; 5 = dense corneal opacity, central opaque area less than 5 mm diameter with complete vascular ring; and 6 = complete dense corneal opacity, eye often with complete blindness. IBK score of the left and right eye were summed for each animal, ranging from a score of 2 for healthy eyes and scores 3 to12 for increasingly diseased eyes. Sizes of the lesions were qualitative estimates. The system of scoring was an attempt to classify the infection according to its severity and whether the lesion was active (and hence of direct effect at the time of scoring) or an inactive scar from a previous infection. No attempt was made to isolate the causal organism.
Phenotypic data for tick counts
Cattle tick (Rhipicephalus microplus) counts (TKC) were performed once every 21 days during the tick season, generally late autumn to early winter and/or spring to early summer for the years 1979 and 1983 through 1997. TKC counts were repeated three times per animal and consisted of counting the number of engorging female ticks 4.5 to 8.0 mm long on one side of every animal following field infestation .
Mean of these three TKC counts was used for the assessment of tick resistance. A count was considered usable if the mean count for that month was equal or greater than 20 ticks per animal side to allow for sufficient variation. The TKC counts were recorded on 875 animals (430 males and 445 females) from weaning (approximately 6 months) to 18 months of age. The animals originated from 57 sires with progeny group sizes ranging from 1 to 64, with an average of 14 animals per sire.
Phenotypic data for helminth egg counts
The number of helminth (unspecified species) eggs per gram (EPG) of faeces were counted on a fresh faecal sample using a McMaster slide under a microscope at 100 × magnification .
The EPG counts were assessed as the mean of the first three EPG counts recorded on each animal after weaning (approx. 6 to 9 months of age) during the period between 1987 and 1998. Those EPG counts were taken three weeks apart during parasite season (from June to July). A total of 727 animals (363 males and 364 females) with EPG counts were available. These animals were offspring of 53 sires (average 14 progeny per sire, ranging from 1 to 70).
Phenotypic data for growth
Four live weight measurements were taken at several ages from birth to 18 months for the period between 1978 and 1998. At birth, calves were individually identified and weighed (birth weight = BWT, kg). Weights for weaning (WWT, kg), yearling (YWT, kg) and final weight (FWT, kg) were recorded when animals were approximately 6, 12 and 18 months of age, respectively. A total 1 298 animals (636 males and 662 females) with weight measurements were available in this study. Animals were progeny of 77 sires (average 17 progeny per sire, ranging from 1 to 70).
Data provided for the present study had previously been recorded as the sum of the left and right IBK scores for each animal, thus, we could not perform ordinal regression (multiple thresholds) analyses, which would otherwise have been the most appropriate analysis. Consequently we converted the summed scores to a binary (0/1) trait, setting score 2 equal to 0 (i.e., IBK absent) and scores 3 to 12 equal to 1 (IBK present). Most health and reproduction traits are usually given as all-or-none (0/1), and it is often assumed that 0/1 traits have an underlying continuous distribution . There are several examples in the literature where ordinal multiple scores of lesions or disorders in animals have been converted to a simple binary form, e.g. for osteochondral bone disease in pigs [21, 22] and for hip dysplasia in dogs .
First, data (phenotypes and fixed effects terms) exploration was performed for each trait using numerical and graphical descriptive analyses (i.e., outliers, distributions and relationships). This has revealed that TKC and EPG counts were not normally distributed. Thus, consistent with previous studies [5, 6], counts were transformed using log10 (mean TKC + 1) for mean TKC counts and a cube-root transformation for mean EPG counts. Quantile-quantile plots indicated that these transformations produced residuals that approached normality (data not shown). All genetic parameter estimates for TKC and EPG counts presented here were based on the transformed scale. Then, significance of environmental fixed effects was assessed (using a statistical significance threshold of P < 0.05) by fitting linear regression models using the lm function for TKC, EPG and weight traits, and the glm function for IBK binary in the R software package . The effects included in the models were age in number of days as a covariate, year when traits were recorded, and sex of the calf for IBK, TKC and EPG traits. For weight traits, the previous lactation status of the dam was included as an additional fixed effect (1 = maiden, 2 = having lactated in the previous year, 3 = dry non-pregnant cow in the previous year, and 4 = dry, cow calved in the previous year but the calf did not survive).
Variance components estimation
where y = vector of phenotypes (TKC, EPG or weight traits), b = vector of fixed effects included in the model, as described previously, X = incidence matrix for fixed effects, a = vector of random animal effects, Z = incidence matrix of random animal effects and e = random residual effects.
The animal and residual terms were assumed to be where and are variances of random animal genetic effect and random residual error effect, respectively. A is a numerator relationship matrix across all animals inferred from the pedigree (n = 3 360), and I is an identity matrix with an order equal to the number of observations in y.
where π = vector of probabilities of IBK present (= 1). The fixed and random effects were the same as used in the linear mixed models. Random effects in the model were assumed to be from multivariate normal distributions as described previously.
where is the estimate of the additive genetic variance and π2/3 is the liability variance scale for the binomial distribution, where π2/3 ~ 3.28987.
Bivariate model for estimation of genetic and phenotypic correlation
In order to estimate genetic and phenotypic relationships among the traits considered, the individual datasets for TKC, EPG and weight were merged with both the IBK pre-weaning dataset and the IBK post-weaning dataset, resulting in a pre-weaning dataset with 535 animals and a post-weaning dataset with 650 animals. Although our intention was to fit a multi-trait model that analysed all traits simultaneously, this was not feasible due to computational limitations; therefore, in order to estimate genetic and phenotypic correlations, bivariate animal models were fitted for various combinations of traits.
where y i = vector of observations for the ith trait (i = 1 or 2), b i = vector of fixed effects for the ith trait (i = 1 or 2), a i = vector of random animal effects for the ith trait (i = 1 or 2), e i = vector of random residual effects for the ith trait (i = 1 or 2), X i = incidence matrices relating records of the ith trait to fixed effects (i = 1 or 2), Z i = incidence matrices relating records of the ith trait to random animal effects (i = 1 or 2).
where are additive genetic variances for traits 1 and 2, is the additive genetic covariance between traits 1 and 2, are residual variances for traits 1 and 2, is the residual covariance between traits 1 and 2,
A and I are as described previously, a and e were assumed to be normally distributed with a mean of zero and (co)variances as specified above.
Results and discussion
Summary statistics for the traits studied
Number of animals
Raw mean tick count
Transformed mean tick count 5
Raw mean EPG
Transformed mean EPG6
Genetic parameter estimates
Heritabilities and standard errors estimated with the univariate model
IBK before weaning
IBK after weaning
Our estimate of the heritability of cube-root transformed EPG counts was moderate (h2 = 0.41 ± 0.09) and similar to that reported in Brahman heifers (0.40 ± 0.12)  and in several Australian studies (0.33 – 0.44) [5, 6]. A study on Angus cattle under a natural challenge situation in New Zealand using a log10 (faecal count + 100) transformation reported a heritability of 0.32 .
In ruminants, immune response varies with the age of the host, the duration of parasite exposure, the number of parasites (dosage) and the nutrition status of the host [30, 31]. In our study, differences in heritability estimates for the parasite-resistance traits may be due to differences in these factors and/or different populations, season effects, and models used.
Heritability estimates ranged from 0.35 ± 0.10 to 0.60 ± 0.09 for weight traits (Table 2). Estimates were 0.60, 0.47, 0.38 and 0.35 for BWT, WWT, YWT and FWT, with standard errors ranging from 0.09 to 0.10. Our estimates were slightly higher than those reported in Australian Angus cattle  but since maternal genetic and permanent environmental effects were not included, the estimates may be biased upward.
Estimated breeding values for IBK
Genetic and phenotypic correlations
Genetic 1 and phenotypic 2 correlations for pre-weaning dataset
Genetic 1 and phenotypic 2 correlations for post-weaning dataset
−0. 50 (0.04)
For the pre-weaning dataset, the phenotypic correlation between IBK and FWT was −0.16 ± 0.04 but close to zero for all other weight traits (BWT, WWT and YWT) (Table 3). For the post-weaning dataset estimates of phenotypic correlations were negative between IBK and weight measurements, ranging from −0.42 ± 0.04 to -0.53 ± 0.03, except between IBK and BWT (0.10 ± 0.05) (Table 4). In general, calves with IBK appear to have a lower weight, which is consistent with previous studies .
Estimates of genetic correlations between IBK and TKC counts were 0.29 ± 0.27 and −0.14 ± 0.23 for the pre-weaning and post-weaning datasets (Table 3 and Table 4). Genetic correlations between IBK and EPG counts were negative at −0.46 ± 0.26 for the pre-weaning dataset The direction of genetic correlations of IBK with TKC and EPG counts were opposite between the two datasets. This may be due to the fact that the calves’ pre-weaning environment is different from the post-weaning environment, including with regard to maternal effects. To our knowledge, genetic relationships between IBK and TKC or EPG counts have not been previously reported.
Genetic correlations between TKC and EPG counts were positive in pre-weaning and post-weaning datasets (0.33 ± 0.23 and 0.31 ± 0.20, Tables 3 and 4). Such positive correlations suggest that pleiotropic genes may be involved in resistance to ecto-parasites and endo-parasites in this cattle population and that cattle that have developed tolerance/resistance to ticks through natural or artificial selection can also acquire tolerance/resistance to parasite infestation. The phenotypic correlations between TKC, EPG and IBK were low for both datasets.
Genetic correlations among weight traits were highly positive and significant for both datasets and ranged from 0.49 to 0.97 (Tables 3 and 4). These estimates were in line with reports from . The positive genetic correlations observed between birth weight and other weight traits is not advantageous since it could lead to calving difficulty (dystocia) and decrease calf survival and welfare of the cow.
Tables 3 and 4 also present genetic parameters for environments with varying levels of environmental stress. For instance, HS calves born in north-east Queensland move from a benign uterine environment to a relatively protected maternal pre-weaning environment and then to a tropically stressful post-weaning environment. For a breed that has evolved in a temperate environment, this represents an extreme contrast in environmental stress. In the pre-weaning environment, the genetic correlations between IBK and the growth traits were consistently small and negative (Table 3), whereas in the post-weaning environment, after removal of the maternal influence on the calves, the genetic correlations are more negative, except with BWT (uterine environment), (Table 4). Similarly, the genetic correlations between the parasite traits and the growth traits all show lower estimates for the pre-weaning dataset compared to the post-weaning dataset (Tables 3 and 4). Estimates of the heritability of growth also differ depending on the level of environmental stress. The heritability estimates for BWT and WWT (uterine and maternal environment, respectively) are higher than, and in some cases almost double the estimates for YWT and FWT (see Table 2). The finding that estimates of heritability are lower in more environmentally stressful environments has been documented previously  and studies of Herefords from the semi-arid environment of Arizona, USA, [36, 37] have also found that genetic parameters vary when estimated in “good’ versus “poor” environments. In an optimal environment, identifying sires with the optimal growth “resource allocation factor”  is straightforward but in a variable stressful environment, animal growth is more complex and must consider both inherent growth (inherent metabolic rate) and fitness (resistance to IBK and parasite infection); hence, a lower heritability for growth and modified genetic correlations in a stressful environment is expected.
Although the traditional quantitative genetic approach used here has assisted us to understand the adaptation of cattle in tropical climates and their ability to deal with environmental stressors, the true value of this information is in identifying genetic architecture affecting adaptation/fitness. The ultimate goal of such studies would be to provide information about molecular markers that could be used in marker-assisted programs and lead to the implementation of selection programs for the simultaneous improvement of growth and tropical adaptation/fitness.
Genetic parameter estimates for infectious bovine keratoconjunctivitis (IBK), also known as ‘pinkeye’, and its genetic and phenotypic correlations with growth and ecto-and endo-parasite traits were investigated. IBK infection is a heritable trait, although its heritability is low to moderate (0.17 ± 0.03 to 0.19 ± 0.06) for pre-weaning and post-weaning calves from Australian tropically adapted Bos taurus cattle. Small to moderate negative genetic correlations between IBK and different weight measurements suggest that the animals that are genetically prone to acquire IBK infection may also be genetically prone to gain weight more slowly. The genetic correlations between IBK and TKC and between IBK and EPG counts display opposite values between pre-weaning and post-weaning datasets, suggesting differences in (maternal) environment between the two stages. However, more accurate estimates are needed; hence further data collection and pedigree recording on IBK, TKC and EPG is recommended. Nevertheless, our estimates of genetic parameters for IBK indicate that there is sufficient genetic variation in beef cattle to design selection and breeding programs for tropical adaptation of these breeds. Recording IBK cases on a routine basis and at a national level would allow inclusion of IBK resistance in the multi-trait breeding goal for the Australian beef industry, as well as in other countries, and would help address animal welfare and health issues.
The authors would like to acknowledge many people involved with the development of the HS line of cattle at the CSIRO Belmont Cattle Research Station. The contribution of Dr John Frisch is particularly acknowledged for the experimental design and data collection. We thank Dr Ed Charmley, program manager and officer in charge of the CSIRO Animal, Food and Health Sciences in Townsville, for useful comments in an earlier version of the manuscript. We also thank two anonymous reviewers and editors for their valuable comments on the submitted version of the manuscript. AAA gratefully acknowledges the financial support from the CSIRO Animal, Food and Health Sciences for his PhD scholarship and publication of this paper.
- O'Neill CJ, Swain DL, Kadarmideen HN: Evolutionary process of Bos taurus cattle in favourable versus unfavourable environments and its implications for genetic selection. Evol Appl. 2010, 3: 422-433. 10.1111/j.1752-4571.2010.00151.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Jonsson NN, Vankan D, Stear MJ, Waine D, Gray CP: Selecting cattle for host resistance to tick (Rhipicephalus microplus) infestation - natural and human-directed selection. Proceedings of the 31st International Conference on Animal Genetics: 20–24 July 2008; Amsterdam. 2008, poster 2171 http://www.isag.us/Docs/2008_ISAG_Amsterdam_P2000.pdfGoogle Scholar
- Frisch JE: Changes occurring in cattle as a consequence of selection for growth-rate in a stressful environment. J Agr Sci. 1981, 96: 23-38. 10.1017/S0021859600031841.View ArticleGoogle Scholar
- Frisch JE, O'Neill CJ, Kelly MJ: Using genetics to control cattle parasites - the Rockhampton experience. Int J Parasitol. 2000, 30: 253-264. 10.1016/S0020-7519(00)00010-2.View ArticlePubMedGoogle Scholar
- Prayaga KC, Henshall JM: Adaptability in tropical beef cattle: genetic parameters of growth, adaptive and temperament traits in a crossbred population. Aust J Exp Agr. 2005, 45: 971-983. 10.1071/EA05045.View ArticleGoogle Scholar
- Henshall JM: A genetic analysis of parasite resistance traits in a tropically adapted line of Bos taurus. Aust J Agr Res. 2004, 55: 1109-1116. 10.1071/AR03085.View ArticleGoogle Scholar
- Barwick SA, Johnston DJ, Burrow HM, Holroyd RG, Fordyce G, Wolcott ML, Sim WD, Sullivan MT: Genetics of heifer performance in 'wet' and 'dry' seasons and their relationships with steer performance in two tropical beef genotypes. Anim Prod Sci. 2009, 49: 367-382. 10.1071/EA08273.View ArticleGoogle Scholar
- Prayaga KC, Corbet NJ, Johnston DJ, Wolcott ML, Fordyce G, Burrow HM: Genetics of adaptive traits in heifers and their relationship to growth, pubertal and carcass traits in two tropical beef cattle genotypes. Anim Prod Sci. 2009, 49: 413-425. 10.1071/EA08247.View ArticleGoogle Scholar
- Slatter DH, Edwards ME, Hawkins CD, Wilcox GE: A national survey of the clinical-features, treatments, and importance of infectious bovine keratoconjunctivitis. Aust Vet J. 1982, 59: 69-72. 10.1111/j.1751-0813.1982.tb02729.x.View ArticlePubMedGoogle Scholar
- Sackett D, Holmes P, Abott K, Jephcoll S, Barber M: Assessing the economic cost of endemic disease on the profitability of Australian beef cattle and sheep producer. Meat and Livestock Australia report, sydney. 2006, AHW.087: 16-18.Google Scholar
- Hansen R: New tools in the battle against pinkeye. Proceedings of Nevada Livestock Production Annual Update: 2001. 2001, University of Nevada Reno UNR Coop. Ext. SP01–01, 5-8.Google 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
- 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
- Kadarmideen HN, Thompson R, Simm G: Linear and threshold model genetic parameters for disease, fertility and milk production in dairy cattle. Anim Sci. 2000, 71: 411-419.Google Scholar
- De Haas Y, Janss LL, Kadarmideen HN: Genetic correlations between body condition scores and fertility in dairy cattle using bivariate random regression models. J Anim Breed Genet. 2007, 124: 277-285. 10.1111/j.1439-0388.2007.00675.x.View ArticlePubMedGoogle Scholar
- Koenig S, Sharifi AR, Wentrot H, Landmann D, Eise M, Simianer H: Genetic parameters of claw and foot disorders estimated with logistic models. J Dairy Sci. 2005, 88: 3316-3325. 10.3168/jds.S0022-0302(05)73015-0.View ArticlePubMedGoogle Scholar
- Turner LB, Harrison BE, Bunch RJ, Porto Neto LR, Li Y, Barendse W: A genome-wide association study of tick burden and milk composition in cattle. Anim Prod Sci. 2010, 50: 235-245. 10.1071/AN09135.View ArticleGoogle Scholar
- Wharton RH, Utech KBW, Turner HG: Resistance to cattle tick, Boophilus-microplus in a herd of Australian Illawarra Shorthorn cattle - its assessment and heritability. Aust J Agr Res. 1970, 21: 163-181. 10.1071/AR9700163.View ArticleGoogle Scholar
- Roberts FHS, O'sullivan PJ: Methods for egg counts and larval cultures for Strongyles infesting the gastro-intestinal tract of cattle. Aust J Agr Res. 1950, 1: 99-102. 10.1071/AR9500099.View ArticleGoogle Scholar
- Falconer DS, Mackey TFC: Introduction to Quantitative Genetics. 1996, Longman Group, EssexGoogle Scholar
- Kadarmideen HN, Schwörer D, Ilahi H, Malek M, Hofer A: Genetics of osteochondral disease and its relationship with meat quality and quantity, growth, and feed conversion traits in pigs. J Anim Sci. 2004, 82: 3118-3127.PubMedGoogle Scholar
- Kadarmideen HN, Janss LLG: Evidence of a major gene from Bayesian segregation analyses of liability to osteochondral diseases in pigs. Genetics. 2005, 171: 1195-1206. 10.1534/genetics.105.040956.PubMed CentralView ArticlePubMedGoogle Scholar
- Wilson BJ, Nicholas FW, James JW, Wade CM, Tammen I, Raadsma HW, Castle K, Thomson PC: Symmetry of hip dysplasia traits in the German Shepherd Dog in Australia. J Anim Breed Genet. 2011, 128: 230-243. 10.1111/j.1439-0388.2010.00903.x.View ArticlePubMedGoogle Scholar
- R Development Core Team: R: A language and environment for statistical computing. 2011, R Foundation for Statistical Computing, Vienna,http://www.R-project.org/,Google Scholar
- Gilmour AR, Gogel BJ, Cullis BR, Thompson R: ASReml User Guide Release 3.0. Hemel Hempstead: VSN International Ltd. 2009Google Scholar
- Ward JK, Nielson MK: Pinkeye (Bovine infectious keratoconjunctivitis) in beef cattle. J Anim Sci. 1979, 49: 361-366.PubMedGoogle Scholar
- Budeli MA, Nephawe KA, Norris D, Selapa NW, Bergh L, Maiwashe A: Genetic parameter estimates for tick resistance in Bonsmara cattle. South Afric J Anim Sci. 2009, 39: 321-327.Google Scholar
- Burrow HM: Variances and covariances between productive and adaptive traits and temperament in a composite breed of tropical beef cattle. Livest Prod Sci. 2001, 70: 213-233. 10.1016/S0301-6226(01)00178-6.View ArticleGoogle Scholar
- Morris CA, Green RS, Cullen NG, Hickey SM: Genetic and phenotypic relationships among faecal egg count, anti-nematode antibody level and live weight in Angus cattle. Anim Sci. 2003, 76: 167-174.Google Scholar
- Wakelin D: Nature and nurture: overcoming constraints on immunity. Parasitology. 1989, 99: S21-S35. 10.1017/S0031182000083396.View ArticlePubMedGoogle Scholar
- Gray GD, Gill HS: Host genes, parasites and parasitic infections. Int J Parasitol. 1993, 23: 485-494. 10.1016/0020-7519(93)90037-Y.View ArticlePubMedGoogle Scholar
- Jeyaruban MG, Johnston DJ, Graser HU: Estimation of genotype × environment interactions for growth, fatness and reproductive traits in Australian Angus cattle. Anim Prod Sci. 2009, 49: 1-8. 10.1071/EA08098.View ArticleGoogle Scholar
- Corbet NJ, Shepherd RK, Burrow HM, van der Westhuizen J, Strydom PE, Bosman DJ: Evaluation of Bonsmara and Belmont Red cattle breeds in South Africa. 1. Productive performance. Aust J Exper Agr. 2006, 46: 199-212. 10.1071/EA05223.View ArticleGoogle Scholar
- Hoffmann AA, Parsons PA: Evolutionary Genetics and Environmental Stress. 1991, Oxford University Press, OxfordGoogle Scholar
- Denise SK, Torabi M: Genetic parameter estimates for postweaning traits of beef - cattle in a stressful environment. J Anim Sci. 1989, 67: 2619-2626.PubMedGoogle Scholar
- Denise SK, Torabi M, Ray DE, Rice R: Genetic parameter estimates for preweaning traits of beef cattle in a stressful environment. J Anim Sci. 1988, 66: 1899-1906.PubMedGoogle Scholar
- Mirkena T, Duguma G, Haile A, Tibbo M, Okeyo AM, Wurzinger M, Soelkner J: Genetics of adaptation in domestic farm animals: A review. Livest Sci. 2010, 132: 1-12. 10.1016/j.livsci.2010.05.003.View ArticleGoogle 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.