Open Access

Genetic susceptibility to infectious disease in East African Shorthorn Zebu: a genome-wide analysis of the effect of heterozygosity and exotic introgression

  • Gemma GR Murray1, 2Email author,
  • Mark EJ Woolhouse1,
  • Miika Tapio3,
  • Mary N Mbole-Kariuki4,
  • Tad S Sonstegard5,
  • Samuel M Thumbi1, 9,
  • Amy E Jennings6,
  • Ilana Conradie van Wyk7,
  • Margo Chase-Topping1,
  • Henry Kiara8,
  • Phil Toye8,
  • Koos Coetzer7,
  • Barend M deC Bronsvoort6 and
  • Olivier Hanotte4
BMC Evolutionary Biology201313:246

https://doi.org/10.1186/1471-2148-13-246

Received: 2 May 2013

Accepted: 4 November 2013

Published: 9 November 2013

Abstract

Background

Positive multi-locus heterozygosity-fitness correlations have been observed in a number of natural populations. They have been explained by the correlation between heterozygosity and inbreeding, and the negative effect of inbreeding on fitness (inbreeding depression). Exotic introgression in a locally adapted population has also been found to reduce fitness (outbreeding depression) through the breaking-up of co-adapted genes, or the introduction of non-locally adapted gene variants.

In this study we examined the inter-relationships between genome-wide heterozygosity, introgression, and death or illness as a result of infectious disease in a sample of calves from an indigenous population of East African Shorthorn Zebu (crossbred Bos taurus x Bos indicus) in western Kenya. These calves were observed from birth to one year of age as part of the Infectious Disease in East African Livestock (IDEAL) project. Some of the calves were found to be genetic hybrids, resulting from the recent introgression of European cattle breed(s) into the indigenous population. European cattle are known to be less well adapted to the infectious diseases present in East Africa. If death and illness as a result of infectious disease have a genetic basis within the population, we would expect both a negative association of these outcomes with introgression and a positive association with heterozygosity.

Results

In this indigenous livestock population we observed negative associations between heterozygosity and both death and illness as a result of infectious disease and a positive association between European taurine introgression and episodes of clinical illness.

Conclusion

We observe the effects of both inbreeding and outbreeding depression in the East African Shorthorn Zebu, and therefore find evidence of a genetic component to vulnerability to infectious disease. These results indicate that the significant burden of infectious disease in this population could, in principle, be reduced by altered breeding practices.

Keywords

Inbreeding depression Heterozygosity Fitness Infectious disease Exotic introgression Local adaptation Cattle

Background

Positive multi-locus heterozygosity-fitness correlations have been observed in a number of natural populations [1, 2]. Genome-wide heterozygosity is correlated with inbreeding within a population [3] and heterozygosity-fitness correlations have therefore been associated with inbreeding depression [4]. Inbreeding may be detrimental to fitness if: (i) fitness related traits are over-dominant, such that the heterozygous genotype results in greater fitness than either homozygote, or (ii) there are recessive alleles present in the population that result in lower fitness [5]. The extent to which the heterozygosity of genetic markers reflects genome-wide heterozygosity and inbreeding is a contentious issue, and heterozygosity-fitness correlations may result from a few loci of major effect in linkage disequilibrium with the genetic markers rather than the effect of inbreeding across the genome [2, 3, 6, 7]. However, when large numbers of markers are used, estimates of genome-wide heterozygosity are greatly improved, are more likely to correlate with inbreeding, and are less likely to correlate with the fitness effects of a few loci [2, 3].

Inbreeding has been observed to have a negative effect on fitness-related traits in many populations [8], e.g. [9, 10], including some European cattle populations [11, 12]. In a number of natural populations susceptibility to infectious disease has been found to be associated with heterozygosity e.g. [1315] and inbreeding e.g. [16]. Infectious disease susceptibility in humans has been found to have a strong genetic component in numerous studies reviewed in [17], with many genes and genetic pathways associated with susceptibility to different diseases. The major histocompatibility complex (MHC) is considered to play a significant role in determining susceptibility to various pathogens [18], and MHC-targeted breeding programmes have been suggested for captive populations [19, 20]. In one natural population it was found that the effect of a single locus linked to the MHC exceeded the effect of genome-wide heterozygosity on fitness [21].

Exotic introgression in a locally adapted population may also reduce fitness. A number of studies have observed a negative effect of introgression on fitness (outbreeding depression) [22, 23]. These may result from either (i) the break-up of co-adapted epistatic interaction between genes, or (ii) the introduction of non-locally adapted alleles [24].

In this study, we tested for associations of heterozygosity and exotic European taurine breed introgression with death and illness due to infectious disease during the first year of life, in a population of East African shorthorn zebu (EASZ) cattle (crossbred Bos taurus x Bos indicus) from western Kenya. More than 500 calves were followed for the first year of their life as part of the Infectious Disease in East African Livestock (IDEAL) project [25]. The calves were closely monitored throughout; deaths and episodes of clinical illness were reported and followed up by trained veterinary staff. During the study, more than 15% of the calves died, with the majority of deaths due to infectious diseases, most commonly East Coast Fever, helminth infections, heartwater and trypanosomiasis [26]. Given the heavy infectious disease burden on these cattle coupled with minimal disease control or treatment, an ability to survive infection is likely to be significantly associated with fitness. Furthermore, more than half the surviving calves were reported to have at least one episode of clinical illness during the study period, and almost all these episodes were attributed to infectious disease (see Methods). We propose that clinical illness is also indicative of a vulnerability to infectious disease that could affect fitness in this environment.

All calves in the IDEAL population were genotyped using the 50K Illumina® BovineSNP50 beadchip v. 1 (55,777 SNPs before quality control). These SNPs are relatively evenly distributed across the genome, with an average of 1,895 on each autosome (ranging from 1,009 to 3,553) and 1,362 on the X chromosome (102 remain unassigned).

The study population was found to be partly introgressed by one or possibly more European taurine breeds [27, 28] (determined from analysis of genome-wide SNP data using STRUCTURE [29]). This reflects the effects of a substantial breed improvement programme, involving the crossing of local cattle with European breeds, that took place in the mid-1990s with intermittent cross-breeding thereafter. As a result, approximately 20% of the calves show levels of introgression consistent with crossing with European breeds ≤ 5 generations ago (fraction European taurine > 2-6) [27].

We would also expect some inbreeding to occur in this population, since mating is largely unmanaged, there is a relatively low number of breeding bulls, and cattle are generally only transported over short distances (only a small fraction of cattle are traded and these over distances typically just 20-30 km). Also, historically, this region has experienced major rinderpest epidemics, which have reduced the effective population size of the EASZ [28]. In this study we calculate the genome-wide heterozygosity for the calves in the IDEAL population using the same SNP data set as was used to calculate exotic introgression [28].

In this study we aim to determine whether inbreeding (as estimated by SNP heterozygosity) and outbreeding (as estimated by the estimated proportion of genetic markers that descend from a European taurine breed) affect the likelihood of death or illness due to infectious disease and therefore whether there is evidence of a genetic basis to vulnerability to infectious disease in EASZ. We use results from earlier studies that have determined the degree of European taurine introgression [28] and describe and analyse the causes of instances of death and illness for all the calves in the IDEAL project [25, 26]. Through analysis of these data we find associations of both heterozygosity and European taurine introgression with vulnerability to infectious disease during the first year of life in this population of EASZ.

Results

The sample size available for analysis was 518 calves (see Methods for inclusion/exclusion criteria), of which 68 (13%) died of infection-related disease within one year of birth (see Methods for more details). Of the survivors, episodes of clinical illness (almost always related to infections) were reported in 243 (54%) (Table 1, Additional file 1: Table S1).
Table 1

Numbers of calves that fall into the different categories used for the statistical analysis of the effects of heterozygosity and introgression on likelihood of death and clinical illness

 

Low heterozygosity

High heterozygosity

 

Introgressed

Non-introgressed

Introgressed

Non-introgressed

Died

2

11

11

44

Survived with clinical episode

6

13

64

160

Survived without clinical episode

1

4

29

173

Cut-offs for high/low heterozygosity and introgressed/non-introgressed are defined in Methods.

The distributions of observed heterozygosity values for non-introgressed (defined as < 2-6 European taurine) and introgressed calves (> 2-6 European taurine) are shown in Figure 1. The mean heterozygosity of non-introgressed calves is 0.28 (with a standard deviation (SD) of 0.01). This is not an unusual level of heterozygosity for a cattle population. Holstein-Friesians (European taurine) have been found to have average heterozygosity of 0.33 (SD = 0.01), Jersey (European taurine) have an average heterozygosity of 0.25 (SD = 0.03), N’dama (African taurine) have an average heterozygosity of 0.17 (SD = 0.03) and Ethiopian Sheko (admixed African taurine and Asian zebu) have an average heterozygosity of 0.26 (SD = 0.0003) [27].
Figure 1

Histogram of frequency of heterozygosity values for European taurine introgressed (red) and non-introgressed (blue) East African Shorthorn Zebu calves.

In order to test whether or not our estimates of heterozygosity are likely to reflect true genome-wide heterozygosity, we randomly sampled our SNPs without replacement 100 times and determined the degree of correlation in the heterozygosity estimates for calves based on one half of the markers and the estimates based on the other half [2, 6]. We found a consistently strong correlation between the estimates (mean Pearson’s r = 0.98, SD = 0.0012, range = 0.977-0.982). Since some SNPs may be in linkage disequilibrium with one another, we repeated the analysis with two smaller subsets (each representing a different random sample, without replacement, of one 10th of the total SNPs). Even with this much smaller set of markers, a strong correlation between the estimates was found (mean Pearson’s r = 0.91, SD = 0.0053, range = 0.89-0.92).

The non-introgressed calves show a skewed distribution of heterozygosity values, with the majority (93%) contained within a near-symmetrical distribution ± 0.010 from a mode of 0.284 (Figure 1). However, there is a long tail of 28 calves with lower values. We categorise calves with values < 0.274 as “low” heterozygosity, assumed to reflect a higher degree of inbreeding than the population norm. Introgressed calves show a different distribution of heterozygosity values that are generally higher than those for non-introgressed calves (Figure 1) (Mann Whitney U test, P < 0.001).

The relationship between heterozygosity and introgression is shown in Figure 2. As expected, there is a clear increase in heterozygosity with increasing introgression. However, in addition to the 28 non-introgressed calves with low heterozygosity, there are several introgressed calves with heterozygosity values that appear low relative to the amount of introgression. In order to investigate the effect of variation in heterozygosity on the fitness of introgressed calves independently of the effect of introgression we define a group of calves with low heterozygosity relative to their degree of introgression. We define this group by finding the best fit linear regression and iteratively excluding calves with heterozygosity values > 0.01 below their expected value. This procedure results in the fitted regression line shown in Figure 2 (with the -0.01 threshold also indicated). Using this criterion, 9 of the 113 introgressed calves are categorised as low heterozygosity. Two of these 9 died and 6 others experienced at least one clinical episode (see Table 1).
Figure 2

Plot of heterozygosity against proportion European introgression for all calves in our analysis. Calves that died are highlighted (red symbols). The best fit relationship between heterozygosity (H) and introgression (I) for non-inbred calves is shown (thick line: regression equation H = 0.284 (95% confidence intervals ±0.0003) + I*0.312 (±0.006), R2 = 0.972). The cut-off to define ‘low’ heterozygosity calves varies with introgression, since the expected or mean heterozygosity varies with introgression. Calves that are more that 0.01 below the expected heterozygosity (blue dashed line) for their amount of introgression are defined as ‘low’ heterozygosity (see Results and Methods). The orange dashed line indicates the cut-off for introgressed calves (2-6%).

We compared both calves that died and those that survived but experienced any episode of clinical illness to a reference category of healthy calves (those that survived and had no reported episode of clinical illness). We considered two categorical predictors of these outcomes, low heterozygosity and being introgressed, using a multinomial logistic regression (described in Methods; results in Table 2). We used these categorical predictors, rather than analysing the estimates directly, due to the non-normality of the data (the distribution of values of both heterozygosity and introgression are highly skewed) and to allow us to robustly account for the relationship between heterozygosity and introgression (Figures 1 and 2, Additional file 1: Table S1).
Table 2

Summary of results of a multinomial logistic regression of the effects of genome-wide heterozygosity and introgression (as two level factors) on the likelihood of death and the likelihood of clinical illness

Phenotype

Dead calves (n = 68)

Ill calves (n = 243)

Factors

Estimate (SE)

OR (95% CI)

P-value

Estimate (SE)

OR (95% CI)

P-value

Low heterozygosity

1.13 (0.27)

9.5 (3.5-27.9)

< 0.001

0.61 (0.26)

3.4 (1.2-9.3)

0.018

Introgressed

0.16 (0.19)

1.4 (0.7-2.9)

0.39

0.43 (0.12)

2.4 (1.5-3.8)

< 0.001

Cut-offs for high/low heterozygosity and introgressed/non-introgressed are defined in Methods. Healthy calves were the reference category. CI = confidence interval. SE = standard error. OR = odds ratio. P-values in bold type are statistically significant.

We first consider the association between predictors and calf death due to infectious disease. For non-introgressed calves (fraction European taurine < 2-6) only there is strong evidence of a relationship between death and low heterozygosity: odds ratio (OR) = 10.8 (95% CIs 3.0-48.1); P < 0.001 (Table 3). For all calves, this relationship remains when both low heterozygosity and introgression are included as predictors: OR = 9.5 (3.5-27.9); P < 0.001 (Table 2). There is no evidence of an interaction between low heterozygosity and introgression (P > 0.05) and there is no evidence of an association between death and degree of introgression (OR = 1.4 [0.7-2.9]; P = 0.39) (Table 2). No significant differences were found in the nature of infectious disease-related deaths among the low heterozygosity calves. These results are robust to variation in the cut-off for introgressed calves (see Methods).
Table 3

Summary of results of Fisher’s exact tests of the effects of low genome-wide heterozygosity (as a two-level factor) on the likelihood of death and the likelihood of clinical illness for the sub-populations of non-introgressed and introgressed calves

Phenotypes

Dead vs healthy calves

Ill vs healthy calves

Subpopulation

OR

95% CI

P-value

OR

95% CI

P-value

Non-introgressed

10.8

3.0-48.1

< 0.001

3.5

1.05-15.1

0.026

Introgressed

5.3

0.24-321

0.213

2.72

0.31-129

0.671

Cut-offs for high/low heterozygosity and introgressed/non-introgressed are defined in Methods. OR = odds ratio. CI = confidence interval. P-values in bold type are statistically significant.

Clinical illness is a less well defined phenotype than death, since it includes a wide range of disease signs and it cannot be unequivocally determined whether or not an instance of clinical illness was the result of infectious disease or not (although it is likely that the overwhelming majority were). Despite this, clinical illness and death are clearly related outcomes. Moreover, because many more calves experienced clinical illness during their first year than died, this outcome may provide increased power to detect small effects. For non-introgressed calves there is evidence of a relationship between illness and low heterozygosity: OR = 3.5 (95% CIs 1.05-15.1); P = 0.026 (Table 3). For all calves, this relationship remains when both low heterozygosity and introgression are included as predictors: OR = 3.4 [1.3-9.3]; P = 0.018 (Table 2). However, for this phenotype there is also strong evidence of an association with introgression, with introgressed calves being more likely to experience illness (OR = 2.4 [1.5-3.8]; P < 0.001). Again there is no evidence of any interaction between low heterozygosity and introgression (P > 0.05). These results are robust to variation in the cut-off for introgressed calves (see Methods).

Discussion

Our results indicate that there is a genetic basis to vulnerability to infectious disease in East African Shorthorn Zebu. We find a positive association between low heterozygosity and both death and illness. Since we use a reasonable proxy for genome-wide heterozygosity (40,457 SNPs across all autosomes), and the estimates of genome-wide heterozygosity are highly correlated for different subsets of SNPs, we consider it likely that our results are indicative of the presence of inbreeding depression in the study population, rather than being a consequence of localised associations between fitness and loci of large effect [6, 7, 17].

Our results indicate the presence of fitness-related genotypes that are either over-dominant or have recessive deleterious alleles in the population. Although the proportion of calves that we classified as “low” heterozygosity was small (approximately 7%), the effect on infectious disease-related mortality was large (OR ~ 9) and there was also an effect on clinical illness (OR ~ 3). Our results give population attributable fractions (for the non-introgressed population) due to low heterozygosity of 14% of mortality and 4% of clinical illness [30]. In principle, therefore, cases of illness and deaths could be reduced by better management practices (e.g. bull rotation).

We found evidence that European taurine introgression affected health as it was a significant predictor of clinical illness (Table 2). As there were only 13 deaths among those with fractional introgression greater than 2-6 we had very little power to detect any effect on death alone. Nonetheless, our results for these two phenotypes are not inconsistent, as the confidence intervals for the odds ratios overlap (see Table 2). We suggest that clinical illness is a marker for vulnerability in this population and consequently that European taurine introgressed calves are more vulnerable in this sense.

We could find no non-genetic confounding factors that could be driving the observed associations between introgression and either death or illness. The introgressed calves (even those with high introgression) could not be distinguished from non-introgressed calves by the field team [25]; exotic introgression was only discovered upon genetic analysis of the calves. And, although the introgressed calves were not found uniformly distributed among the 20 sub-locations of the study [28], there is no evidence of geographic variation in the rate of mortality or of clinical illness.

Since the EASZ ancestors of the study population have been present in East Africa for many generations [28], it might be expected that they may have in some way adapted to the specific challenges posed by that environment, in particular to the infectious diseases present in that region. Infectious disease was the primary cause of first year mortality and illness in this population, and thus the ability to withstand infection is likely to be of significant selective advantage. On the other hand, the introgressed animals descend, in part, from European cattle, which have not had the same history of selective pressures. It has been found that European breeds are more susceptible than local breeds to some infectious diseases that are prevalent in Africa, such as East Coast Fever [31, 32], bovine tuberculosis [33, 34], and trypanosomiasis [35, 36]. Thus, the European breeds are, in this important respect, less well adapted to the East African environment, and introgressed animals, while more heterozygous due to outbreeding, are also less locally adapted. While increased heterozygosity is likely to be beneficial if fitness related traits are over-dominant or there are recessive deleterious alleles in the population, the breaking up of or loss of locally adapted genotypes is likely to be detrimental to fitness. Both these competing forces appear to be at work in this population of East African Shorthorn Zebu calves in western Kenya.

Conclusions

In conclusion, in this indigenous cattle population from the tropics we find strong associations between genome-wide heterozygosity and both death and clinical illness due to infectious diseases. There is evidence that introgression of European taurine breeds also increases vulnerability to infectious diseases. We therefore believe that we have observed the effect of both inbreeding and outbreeding depression, and have found evidence of a genetic component to vulnerability to infectious disease and of adaptation in the East African Shorthorned Zebu population in western Kenya. This implies that the significant burden of infectious disease in this population could, in principle, be reduced by better breeding practices.

Methods

Data set

The IDEAL study population consisted of 548 calves, which were monitored closely throughout their first year of life [25]. Three of these calves were excluded from this analysis because they received treatment for disease. The calves were sampled from a region in western Kenya (of approximately 45 × 90 km) (covering some of the Busia, Teso, Siaya, Butere/Mumias and Bungoma districts). This region traversed four different types of agro-ecological zones. The region was divided into sub-locations of approximately 10 × 10 km, of which 20 were randomly selected (within each agro-ecological zone) for study [25]. Small-scale cattle farming is common in this region and the sub-locations contain approximately 80 to 90 households, of which approximately 60% own cattle (generally ranging from 1-10 cattle per household). New born calves were recruited into the study from these 20 sub-locations over a period of 2 years (2007-2009). Calves that showed visible signs of exotic breed introgression were excluded, and only calves from small farms, where the cattle spent at least part of their time outside, exposed to infectious disease, were recruited [25].

Genetic data

All calves in the IDEAL population were genotyped using the 50K Illumina® BovineSNP50 beadchip v. 1 including 55,777 SNPs before quality control. Standard quality control procedures were applied to the data prior to analysis reviewed in [37]. This was implemented using GenABEL. A minor allele frequency for SNPs cut-off of 1% was applied and a SNP call-rate cut-off of 90%. A cut-off of call rate for an animal of 90% was applied. And an identity by state (IBS) threshold of 90% and a cut-off odds of 1000-1 was applied to decide whether a calf should be excluded based on sex/X-linked marker data inconsistency (which both eliminate errors due to sample mis-identification). These quality control thresholds left a total of 40,457 autosomal SNPs for subsequent analysis. A Hardy-Weinberg equilibrium (HWE) cut-off was not applied, since a divergence from HWE can result from factors other than genotyping error (in this case introgression).

As a result of quality control, 9 calves and 13,856 of the 54,313 initial autosomal SNPs had to be excluded from the 50K SNP chip data set due to unsuccessful genotyping or mis-identification of individuals. This reduced the genotyping data set for this study to 536 calves and 40,457 autosomal SNPs.

Fitness phenotype – death

The best proxy for fitness in this study is survival to first year. Eighty-eight of the original 548 calves died during the study, and all but nine of these were subject to a detailed post-mortem examination [25, 26]. In most cases the cause of death was found to be infectious disease. These results are described in detail elsewhere [25, 26], but the most frequent causes of death were East Coast Fever (caused by the tick-borne protozoan Theileria parva), heartwater (caused by the tick-borne rickettsia Ehrlichia ruminantium) and helminthosis (especially associated with the nematode Haemonchus placei). Of the 88 deaths, 9 were found to have causes unlikely to be associated with any heritable fitness, such as trauma and starvation, and 10 had an undetermined cause of death. These calves were excluded from the analysis, so as to make the phenotype more likely to have a consistent (although likely complex) genetic basis. This resulted in a data set, after quality control, of 68 calves that died due to infectious disease, and 450 calves that survived their first year (Table 1, Additional file 1: Table S1). Follow-up studies indicated that the death rate over the next 1-2 years was much lower than during the first year.

Fitness phenotype – clinical illness

All calves were monitored for signs of clinical illness throughout the study. Episodes of clinical illness were reported by the calf’s owner, visiting animal health workers or the project staff and were subject to detailed investigation using a standardised protocol and, where necessary, followed up by the project’s veterinary staff [25]. Of the surviving calves in the study, 243 experienced at least one reported episode of clinical illness. Clinical signs ranged from mild to severe. Treatment was provided if the calf was in distress and its keeper agreed; treated calves (n = 3) were subsequently excluded from the study. In a few instances a severely ill calf was euthanized; this was regarded as equivalent to a death. Almost all episodes of clinical illness could be attributed to infections, however causal links are difficult to definitively establish and so episodes cannot be unequivocally divided into those that were caused by infectious disease and those that were not. For this reason all episodes of clinical illness were used to distinguish between healthy and unhealthy calves (where this difference in health can generally be attributed to variation in susceptibility to infectious disease). This division results in 207 calves that neither died nor experienced any reported episodes of clinical illness up to one year old.

European taurine introgression

Using a cut-off of fraction 2-6 introgression we identified 113 introgressed calves, with a maximum of 36.1% European taurine introgression. The introgressed calves were not found uniformly distributed among the 20 sub-locations of the study: they are found mainly in 12 sub-locations in the northern and central regions of the study area [28]. The most likely origin of this European taurine introgression in the study area is a number of recent and continuing breed improvement programmes that use exotic animals and semen, and cattle markets in the study area continue to sell crossbred cattle [28]. The most active programme took place in the mid 1990s, equivalent to approximately 5 EASZ generations before our study took place. This is reflected in our cut-off of fraction 2-6 introgression, corresponding to the expected level of introgression after 5 generations of back-crossing with indigenous cattle since the F1 generation. However, to confirm the robustness to the precise value of the cut-off, we repeated all analyses using alternative cut-off values but obtained very similar results.

Heterozygosity

Individual mean genome-wide heterozygosity was computed as the proportion of all autosomal genotyped loci that passed quality control (40,457 SNPs) that were heterozygous. This was implemented using GenABEL [38]. To test for the whether this is likely to be a good estimate of genome-wide heterozygosity, we randomly sampled SNPs without replacement 100 times and calculated the Spearman correlation coefficient of the heterozygosity estimates for calves based on one half of the markers and the estimates based on the other half [2, 6]. Since SNPs may in linkage disequilibrium with one another we repeated the analysis comparing smaller subsets.

Statistical methods

The comparisons between healthy calves (n = 217) and those that died as a result of infectious disease (n = 68) or survived but experienced clinical illness (n = 243) was investigated using a 3-level multinomial logistic regression analysis using Proc Logistic SAS version 9.3 (SAS Institute Inc., Cary, NC). Both heterozygosity and introgression were included as predictors and were categorised for the purposes of analysis. Using a cut-off of fraction 2-6 introgression we identified 113 introgressed calves. We repeated all analyses using alternative cut-offs of 2-7, 2-5 and 2-4 (6, 4 and 3 generations of crossing respectively). We categorised non-introgressed calves with heterozygosity < 0.274 (more than 0.10 below the mean for this group) as low heterozygosity, assumed to reflect a higher degree of inbreeding than the population norm (28 calves). We defined low heterozygosity calves in the introgressed group by finding the best fit linear regression of heterozygosity against introgression and iteratively excluding calves with heterozygosity values > 0.01 below their expected value (resulting in 9 of the 113 introgressed calves being categorised as low heterozygosity).

We conducted confirmatory univariate analyses using Fisher’s exact test (StatXact version 8, Cytel Software Corp, Cambridge, MA, USA) for the non-introgressed and introgressed categories, allowing us to investigate the effect of heterozygosity for these separately.

Availability of supporting data

Genotyping data can be found at Dryad (http://datadryad.org/resource/doi:10.5061/dryad.bc598).

Declarations

Acknowledgements

The IDEAL project was funded by the Wellcome Trust (grant reference WT079445MF). USDA-ARS bovine functional laboratory and Geneseek veterinary diagnostics provided invaluable technical assistance for the genotyping of the samples.

This project was also partially supported by USDA Agricultural Research Service project 1265-31000-104-00D. Mention of trade names or commercial products in this article is solely for the purpose of providing specific information and does not imply recommendation or endorsement by the U.S. Department of Agriculture.

Authors’ Affiliations

(1)
Institute of Evolutionary Biology, and Centre for Immunity, Infection and Evolution, University of Edinburgh
(2)
Department of Genetics, University of Cambridge
(3)
Biotechnology and Food Research, MTT Agrifood Research Finland
(4)
School of Life Sciences, University of Nottingham, University Park
(5)
United States Department of Agriculture, Bovine Functional Genomics Laboratory, Agricultural Research Service
(6)
Roslin Institute, Easter Bush, University of Edinburgh, Roslin
(7)
Department of Veterinary Tropical Diseases, Faculty of Veterinary Science, University of Pretoria
(8)
International Livestock Research
(9)
Paul G Allen School for Global Animal Health, Washington State University

References

  1. Coltman DW, Slate J: Microsatellite measures of inbreeding: a meta-analysis. Evolution. 2003, 57: 971-983.PubMedView ArticleGoogle Scholar
  2. Chapman JR, Nakagawa S, Coltman DW, Slate J, Sheldon BC: A quantitative review of heterozygosity-fitness correlations in animal populations. Mol Ecol. 2009, 18: 2746-2765. 10.1111/j.1365-294X.2009.04247.x.PubMedView ArticleGoogle Scholar
  3. Slate J, David P, Dodds KG, Veenvliet BA, Glass BC, Broad TE, McEwan JC: Understanding the relationship between the inbreeding coefficient and multilocus heterozygosity: theoretical expectations and empirical data. Heredity. 2004, 93: 255-265. 10.1038/sj.hdy.6800485.PubMedView ArticleGoogle Scholar
  4. Keller LF, Waller DM: Inbreeding effects in wild populations. Trends Ecol Evol. 2002, 17: 230-241. 10.1016/S0169-5347(02)02489-8.View ArticleGoogle Scholar
  5. Charlesworth D, Charlesworth B: Inbreeding depression and its evolutionary consequences. Annu Rev Ecol Syst. 1987, 18: 237-268. 10.1146/annurev.es.18.110187.001321.View ArticleGoogle Scholar
  6. Balloux F, Amos W, Coulson T: Does heterozygosity estimate inbreeding in real populations?. Mol Ecol. 2004, 13: 3021-3031. 10.1111/j.1365-294X.2004.02318.x.PubMedView ArticleGoogle Scholar
  7. Pemberton J: Measuring inbreeding depression in the wild: the old ways are the best. Trends Ecol Evol. 2004, 19: 613-615. 10.1016/j.tree.2004.09.010.PubMedView ArticleGoogle Scholar
  8. Crnokrak P, Roff DA: Inbreeding depression in the wild. Heredity. 1999, 83: 260-270. 10.1038/sj.hdy.6885530.PubMedView ArticleGoogle Scholar
  9. Walling CA, Nussey DH, Morris A, Clutton-Brock TH, Kruuk LEB, Pemberton JM: Inbreeding depression in red deer calves. BMC Evol Biol. 2011, 11: 318-10.1186/1471-2148-11-318.PubMed CentralPubMedView ArticleGoogle Scholar
  10. Coltman DW, Pilkington JG, Smith JA, Pemberton JM: Parasite-mediated selection against inbred soay sheep in a free-living, island population. Evolution. 1999, 53: 1259-1267. 10.2307/2640828.View ArticleGoogle Scholar
  11. Sewalem A, Kistemaker GJ, Miglior F, Van Doormaal BJ: Analysis of inbreeding and its relationship with functional longevity in canadian dairy cattle. J Dairy Sci. 2006, 89: 2210-2216. 10.3168/jds.S0022-0302(06)72291-3.PubMedView ArticleGoogle Scholar
  12. Thompson JR, Everett RW, Hammerschmidt NL: Effects of inbreeding on production and survival in holsteins. J Dairy Sci. 2000, 83: 1856-1864. 10.3168/jds.S0022-0302(00)75057-0.PubMedView ArticleGoogle Scholar
  13. Luikart G, Pilgrim K, Visty J, Ezenwa VO, Schwartz MK: Candidate gene microsatellite variation is associated with parasitism in wild bighorn sheep. Biol Lett. 2008, 4: 228-231. 10.1098/rsbl.2007.0633.PubMed CentralPubMedView ArticleGoogle Scholar
  14. Lyons EJ, Amos W, Berkley JA, Mwangi I, Shafi M, Williams TN, Newton CR, Peshu N, Marsh K, Scott JAG, Hill AV: Homozygosity and risk of childhood death due to invasive bacterial disease. BMC Med Genet. 2009, 10: 55-PubMed CentralPubMedView ArticleGoogle Scholar
  15. Rijks JM, Hoffman JI, Kuiken T, Osterhaus ADME, Amos W: Heterozygosity and lungworm burden in harbour seals (Phoca vitulina). Heredity. 2008, 100: 587-593. 10.1038/hdy.2008.18.PubMedView ArticleGoogle Scholar
  16. Butcher AR, Palethorpe HM, Grove DI: The susceptibility of inbred mice to infection with Brachylaima cribbi (Digenea: Brachylaimidae). Parasitol Int. 2002, 51: 109-115. 10.1016/S1383-5769(01)00102-7.PubMedView ArticleGoogle Scholar
  17. Chapman SJ, Hill AVS: Human genetic susceptibility to infectious disease. Nat Rev Genet. 2012, 13: 175-188.PubMedGoogle Scholar
  18. Horton R, Wilming L, Rand V, Lovering RC, Bruford EA, Khodiyar VK, Lush MJ, Povey S, Talbot CC, Wright MW, Wain HM, Trowsdale J, Ziegler A, Beck S: Gene map of the extended human MHC. Nat Rev Genet. 2004, 5: 889-899. 10.1038/nrg1489.PubMedView ArticleGoogle Scholar
  19. Hughes AL: MHC polymorphism and the design of captive breeding programs. Conserv Biol. 1991, 5: 249-251. 10.1111/j.1523-1739.1991.tb00130.x.View ArticleGoogle Scholar
  20. Wedekind C: Sexual selection and life-history decisions: Implications for supportive breeding and the management of captive populations. Conserv Biol. 2002, 16: 1204-1211. 10.1046/j.1523-1739.2002.01217.x.View ArticleGoogle Scholar
  21. Lenz TL, Mueller B, Trillmich F, Wolf JBW: Divergent allele advantage at MHC-DRB through direct and maternal genotypic effects and its consequences for allele pool composition and mating. Proc R Soc B. 2013, 280: 176-View ArticleGoogle Scholar
  22. Fournier-Level A, Korte A, Cooper MD, Nordborg M, Schmitt J, Wilczek AM: A Map of Local Adaptation in Arabidopsis thaliana. Science. 2011, 334: 86-89. 10.1126/science.1209271.PubMedView ArticleGoogle Scholar
  23. Olano-Marin J, Mueller JC, Kempenaers B: Heterozygosity and survival in blue tits (Cyanistes caeruleus): contrasting effects of presumably functional and neutral loci. Mol Ecol. 2011, 20: 4028-4041. 10.1111/j.1365-294X.2011.05177.x.PubMedView ArticleGoogle Scholar
  24. Szulkin M, David P: Negative heterozygosity-fitness correlations observed with microsatellites located in functional areas of the genome. Mol Ecol. 2011, 20: 3949-3952. 10.1111/j.1365-294X.2011.05277.x.PubMedView ArticleGoogle Scholar
  25. de C Bronsvoort BM, Thumbi SM, Poole EJ, Kiara H, Tosas Auguet O, Handel IG, Jennings A, Conradie I, Mbole-Kariuki MN, Toye PG, Hanotte O, Coetzer J, Woolhouse ME: Design and descriptive epidemiology of the Infectious Diseases of East African Livestock (IDEAL) project, a longitudinal calf cohort study in western Kenya. BMC Vet Res. 2013, 9: 171-10.1186/1746-6148-9-171.View ArticleGoogle Scholar
  26. Thumbi S, Bronsvoort B, Kiara H, Toye P, Poole J, Mbole-Kariuki MN, van Wyk Conradie I, Jennings AE, Handel I, Coetzer JAW, Steyl JCA, Hanotte O, Woolhouse MEJ: Parasite co-infections and their impact on survival of indigenous cattle. PLoS One. In PressGoogle Scholar
  27. Mbole-Kariuki MN: PhD Thesis. Genomic diversity of East African Shorthorn zebu in Western Kenya. 2012Google Scholar
  28. Mbole-Kariuki MN, Sonstegard T, Orth A, Thumbi S, Bronsvoort B, Kiara H, Toye P, Conradie I, Jennings AE, Coetzer K, Woolhouse MEJ, Hanotte O, Tapio M: Genome–wide analysis reveal the ancient and recent genetic admixture of East African Shorthorn Zebu. Under reviewGoogle Scholar
  29. Pritchard JK, Stephens M, Donnelly P: Inference of Population Structure Using Multilocus Genotype Data. Genetics. 2000, 155: 945-959.PubMed CentralPubMedGoogle Scholar
  30. Dohoo I, Martin W, Stryhn H: Veterinary Epidemiologic Research. 2009, Canada: AVC, 2Google Scholar
  31. Jensen K, Paxton E, Waddington D, Talbot R, Darghouth MA, Glass EJ: Differences in the transcriptional responses induced by Theileria annulata infection in bovine monocytes derived from resistant and susceptible cattle breeds. Int J Parasitol. 2008, 38: 313-325. 10.1016/j.ijpara.2007.08.007.PubMedView ArticleGoogle Scholar
  32. Glass EJ, Preston PM, Springbett A, Craigmile S, Kirvar E, Wilkie G, Brown CGD: Bos taurus and Bos indicus (Sahiwal) calves respond differently to infection with Theileria annulata and produce markedly different levels of acute phase proteins. Int J Parasitol. 2005, 35: 337-347. 10.1016/j.ijpara.2004.12.006.PubMedView ArticleGoogle Scholar
  33. Vordermeier M, Ameni G, Berg S, Bishop R, Robertson BD, Aseffa A, Hewinson RG, Young DB: The influence of cattle breed on susceptibility to bovine tuberculosis in Ethiopia. Comp Immunol Microbiol Infect Dis. 2012, 35: 227-232. 10.1016/j.cimid.2012.01.003.PubMed CentralPubMedView ArticleGoogle Scholar
  34. Ameni G, Aseffa A, Engers H, Young D, Gordon S, Hewinson G, Vordermeier M: High prevalence and increased severity of pathology of bovine tuberculosis in Holsteins compared to zebu breeds under field cattle husbandry in central Ethiopia. Clin Vaccine Immunol. 2007, 14: 1356-1361. 10.1128/CVI.00205-07.PubMed CentralPubMedView ArticleGoogle Scholar
  35. Naessens J, Teale AJ, Sileghem M: Identification of mechanisms of natural resistance to African trypanosomiasis in cattle. Vet Immunol Immunopathol. 2002, 87: 187-194. 10.1016/S0165-2427(02)00070-3.PubMedView ArticleGoogle Scholar
  36. Murray M, Morrison WI, Whitelaw DD: Host susceptibility to African trypanosomiasis: trypanotolerance. Adv Parasitol. 1982, 21: 1-68.PubMedView ArticleGoogle Scholar
  37. Anderson CA, Pettersson FH, Clarke GM, Cardon LR, Morris AP, Zondervan KT: Data quality control in genetic case-control association studies. Nat Protocols. 2010, 5: 1564-1573. 10.1038/nprot.2010.116.View ArticleGoogle Scholar
  38. Aulchenko YS, Ripke S, Isaacs A, van Duijn CM: GenABEL: an R library for genome-wide association analysis. Bioinformatics. 2007, 23: 1294-1296. 10.1093/bioinformatics/btm108.PubMedView ArticleGoogle Scholar

Copyright

© Murray et al.; licensee BioMed Central Ltd. 2013

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.