Skip to main content
  • Research article
  • Open access
  • Published:

Detecting the effects of selection at the population level in six bovine immune genes



The capacity of a species or population to respond to and survive novel infectious disease challenge is one of the most significant selective forces shaping genetic diversity and the period following animal domestication was likely one of the most important in terms of newly emerging diseases. Inter-specific genome-wide comparison has suggested that genes, including cluster of differentiation 2 (CD2), ADP-ribosyltransferase 4 (ART4), tyrosine kinase binding protein (TYROBP) and interleukins IL2, IL5, IL13, may have undergone positive selection during the evolution of the bovine lineage. Past adaptive change implies that more recent variation may have also been subject to selective forces.


In this paper, we re-sequence each of these genes in cattle cohorts from Europe, Africa and Asia to investigate patterns of polymorphism at the population level. Patterns of diversity are higher within Bos indicus suggesting different demographic history to that of Bos taurus. Significant coding polymorphism was observed within each of the cell-surface receptors. In particular, CD2 shows two divergent haplotypes defined by a series of six derived nonsynonymous substitutions that are significantly clustered on the extracellular surface of the protein and give significant values for Fay and Wu's H, strongly suggesting a recent adaptive history. In contrast, the signaling molecules (especially IL13) display outlying allele frequency spectra which are consistent with the effects of selection, but display negligible coding polymorphism.


We present evidence suggestive of recent adaptive history in bovine immune genes; implying some correspondence between intra- and inter-specific signals of selection. Interestingly, three signaling molecules have negligible nonsynonymous variation but show outlying test statistics in contrast to three receptors, where it is protein sequence diversity that suggests selective history.


Domestication has involved severe and novel selective pressures on cattle and other domesticates, undoubtedly leaving behind population genetic signatures of adaptation [1]. In particular, the epidemiology of infectious diseases would have changed during domestication due to sharp increases in population densities and the new proximity of previously separated species. The effect of the domestication process in terms of emerging diseases is often discussed with respect to human populations [2], however, similar processes are as likely to have shaped the genomes of domesticated animals. By screening the bovine genome for selective signatures associated with immunity or disease susceptibility, we may be able to identify those genes that have been of critical importance to the development of disease resistance [3].

We have previously reported a comparative genomics study identifying several bovine genes, including cluster of differentiation 2 (CD2) and ADP-ribosyltransferase 4 (ART4), as having significant evidence of adaptive evolution on the evolutionary lineage leading to modern cattle from the bovine-pig common ancestor [4]. Briefly, we investigated evidence of variable selective pressure on the bovine lineage in a dataset of approximately 3,000 orthologs from human, mouse, cow and pig. A gene was inferred to be subject to adaptive evolution on the bovine lineage when a model of variable selective pressure specifically on that lineage was significantly favoured over the alternative model and the ratio of nonsynonymous substitutions per nonsynonymous site to synonymous substitutions per synonymous site, ω, was significantly greater than 1. CD2 (ω = 3.858) and ART4 (ω = 1.356) were ranked first and third in a list of significant scores ranked by this estimated ratio.

Emerging data from the bovine genome project has allowed us to similarly investigate a larger number of bovine genes. Interleukins 2, 5, 13 (IL2, IL5 and IL13), and tyrosine kinase binding protein (TYROBP) each gave preliminary evidence suggestive of adaptive evolution on the bovine lineage (Lynn et al, unpublished).

Several studies of human genes, which have signatures of adaptive evolution between species, found that many of these genes have also been subject to more recent population selection [57]. To assess the level of within-species polymorphism and to identify potential population-specific selective signatures in these six bovine immune genes we have re-sequenced exonic, intron and intergenic regions of each in population samples from Africa, Asia and Europe, plus an outgroup species, bison. Thirty-nine individual cattle, representing 16 different breeds were included in these samples. The three continental groups, which probably reflect the products of separate domestication events, certainly have had different post-domestic histories and have endured markedly different infectious challenges [8, 9].


Patterns of sequence diversity

The three samples chosen for re-sequencing were multi-breed collections of European and West African Bos taurus plus B. indicus sampled across four breeds of South Asian origin. Each of these three continental populations is thought to have separate histories since domestication and here the legacy of migration was minimized by careful selection from populations known from prior data to have experienced minimal introgression [10]. Summary statistics and tests of neutrality are given for each gene within each population sample in Table 1. Each locus showed polymorphism, within and between the three continental samples (alignments are given in Additional file 1). Interestingly neither IL2 nor IL5 showed any nonsynonymous sequence changes and IL13 displayed only a single nonsynonymous change present in one African and four B. indicus chromosomes. There is a marked tendency for the B. indicus sample to give the highest sequence diversity (five from six loci) and also the highest haplotype diversity (five from six loci).

Table 1 Neutrality statistics for six genes re-sequenced in European and African B.taurus and B.indicus cattle

Neutrality tests at each locus

Interleukin 13

Interleukin 13 (IL13) is a cytokine which plays a critical role in allergic inflammation and in susceptibility to parasites, such as, helminths, schistosomes and nematodes [1114]. The pattern of diversity at the IL13 locus is consistent with non-neutral evolution in the African population sample, with significant results in four tests of neutrality (Table 1). Significantly negative values of Tajima's D, Fu and Li's, D and F statistics, and Fu's Fs were obtained for this population, suggesting positive selection; nucleotide diversity was also somewhat reduced relative to the other samples [1517]. These significance values are calculated on coalescent simulations assuming a constant effective population size. In order to consider the demography of domestication, which probably included a significant post-domestic population expansion, we also calculated Tajima's D and Fay and Wu's H conditioned on a significant population expansion. The IL13 African D value remained below the 5th percentile. Aside from the relatively robust signature of positive selection in the African population, Fu and Li's D and Fu's Fs also showed somewhat less significant departures from neutrality in European samples (Table 1).

IL13 displays almost no coding sequence polymorphism with only a single (nonsynonymous) SNP in some B. indicus samples and one African heterozygote (possibly due to B. indicus introgression). The nonsynonymous SNP (A7T) is not predicted to have an affect on protein structure by the SIFT or PolyPhen algorithms [18, 19]. Apart from one synonymous substitution compared to plains bison, there are no differences between the bovine IL13 coding sequence and this outgroup species, although there are several intergenic and intronic SNPs.

Interleukin 5

This gene is part of a cytokine cluster on Bta7 along with IL13 and IL4, all with known involvement in the Th2-type response. Co-expression at this cluster is controlled by elements over a 120 kb range in humans [20]. The gene product is a main regulator of eosinophil maturation and activation and has been implicated as influencing Shistosoma infection in humans [21]. The gene showed moderate nucleotide and haplotypic diversity but no nonsynonymous polymorphisms within the sampled populations. Fay and Wu's H tests were significant in both the African (P = 0.001) and European (P = 0.002) samples and remained outside the 95% interval in comparison to values generated under a simulated domestication bottleneck and expansion.

Interleukin 2

IL2 is a central regulator of the immune system involved in processes including production of cytotoxic T-cells, proliferation of immunoglobulins by B-lymphocytes and induction and maintenance of natural killer (NK) cells [22]. Three neutrality tests, based on the allele frequency spectrum (Tajima's D, Fay and Wu's H and Fu's Fs), show significant deviation from neutrality at the IL2 locus in European populations (Table 1) [1517, 23]. On comparison to the coalescent model, modeled on a population expansion akin to domestication, Fay and Wu's H retained an outlying value (P < 0.025) but Tajima's D did not retain significance.

The African population contains only a single SNP over the 3461 bp re-sequenced at this locus; the lowest level of sequence diversity (0.0001) within any locus-sample combination (Table 1). No nonsynonymous polymorphism was observed at this locus.

Tyrosine Kinase Binding Protein

Tyrosine kinase binding protein (TYROBP; also known as DAP12) is involved in the activation of Natural Killer (NK) cell anti-viral responses, via 'non-inhibitory' members of killer-cell inhibitory receptor (KIR) family [2427]. NK cell KIRs have been shown to be rapidly evolving in a wide range of species, including cattle, and have undergone lineage-specific expansions or contractions of the gene family [2830]. It is tempting to speculate that adaptive evolution of bovine TYROBP, may be driven by co-evolution with rapidly evolving members of the KIRs, with which it interacts.

The results of the standard neutrality tests at this locus are not as clear as for the cytokines discussed previously (Table 1). Fu's Fs test [17] and Tajima's D are significantly negative in the European population, and Fay and Wu's H test [23] suggests an excess of derived SNPs in the African and European samples. Significance at the H test proved robust to a population expansion model but the D value fell within the 95% confidence interval.

TYROBP has two nonsynonymous SNPs within cattle, one only within the African sample (S12C) and one derived variant which is almost fixed within both the African and European samples (Q24L). Both of these occur in the signal peptide and are predicted by SIFT to have an impact on protein structure and/or function [18]. An additional nonsynonymous substitution (P106A), predicted by PolyPhen as 'possibly damaging', occurs between cattle and plains bison, and is located in the cytoplasmic domain of TYROBP.

ADP-ribosyltransferase 4

Human ADP-ribosyltransferase 4 (ART4) is an erythrocyte glycoprotein that has been identified as the Dombrock blood group [31]. The two major antigens of this blood group (Doa and Dob) have significantly different frequencies in human populations and are associated with severe transfusion reactions [32]. Three other high incidence antigens (Gya, Hy and Joa) are also carried on the Dombrock blood group molecule.

Bovine ART4 was found to have five nonsynonymous and six synonymous SNPs within cattle populations and one nonsynonymous and synonymous substitution compared to plains bison. More specifically, four of the five nonsynonymous changes (G5R, R41M, I130V, S213A) and all of the within-cattle synonymous changes occur only in the B. indicus sample. Indeed, the B. indicus population exhibits much higher nucleotide diversity with 33 segregating SNPs compared to just 6 within each of the African and European populations (Table 1).

Cattle do not appear to be polymorphic at the orthologous position to the human mutation (N265D) which is responsible for the Doa/Dob antigenicity. Bovine ART4 is, however, polymorphic at the position associated with the Joa antigen in humans (human T117I; bovine I130V) [32]. The antigenicity associated with this polymorphism in humans suggests that mutations at this position may result in structural alteration of the molecule. Two of the other nonsynonymous changes (G5R, G126W) are also predicted to alter protein structure by SIFT [18]. The I130V and G5R mutations are found only in B. indicus cattle while G126W is found only in African samples.

Cluster of Differentiation 2

Antigen recognition by T cells is one of the crucial steps in controlling adaptive immune responses. Cluster of Differentiation 2 (CD2) is a T cell and natural killer cell surface protein that optimizes T cell activation through interaction with its counter receptor; CD58 in non-rodent mammals [3335] and CD48 in rodents [36]. CD2 and other similar cell surface proteins have low affinity for their counter receptors [37], a fact that may be exploited by viruses which can evolve a higher affinity for these receptors in order to invade the host cell [38]. We have previously identified the specific codon sites that have undergone adaptive changes in CD2 in the bovine lineage and revealed that they are almost exclusively limited to the extracellular portion of the molecular structure [4]. We speculated that a pathogen or more specifically, a virus interaction could have driven the rapid evolution observed in a region of the extracellular domain not involved in CD58/48 binding [4].

Here we identified a significant Fay and Wu's H value within the European sample, indicating an excess of derived variants (P < 0.05; Table 1). This significance was robust to imposition of a domestication/post domestic expansion demographic model on test simulations; and using these as a null the African sample attained significance also. More strikingly, bovine CD2 displayed high coding polymorphism. Eleven nonsynonymous SNPs were present within the sampled populations, and there were an additional three nonsynonymous substitutions compared to the outgroup species, plains bison. Seven of the eleven within-cattle nonsynonymous SNPs and all of the outgroup nonsynonymous substitutions are located in the extracellular domain of CD2, a region which we have previously shown to be under adaptive evolution in cattle [4]. This genomic region (exon 3) displayed a clear peak in sequence diversity (Figure 1). This is reinforced by a visual inspection of haplotypes; six of these nonsynonymous changes lie within 130 bp and segregate in two strikingly divergent exon 3 haplotypes which are common in all three populations.

Figure 1
figure 1

Sliding window analysis of per-site sequence diversity (π) within each population for resequenced regions of CD2. A window size of 150 bp and a slide of 25 bp were used. Regions sequenced for each locus are denoted as filled and exons as thicker rectangles in the scaled locus representation underneath. Plots are colour coded with reference to population (blue, African Bos taurus; red, European Bos taurus; yellow, Bos indicus). Note the peak in diversity within each population at exon 3.

The locations of the seven extracellular domain nonsynonymous SNPs occurring within cattle populations were plotted on the structure of human CD2 based on their orthologous positions. It was noted that these sites appeared to cluster together in a particular region of the extracellular domain in the C2-set immunoglobulin domain (Figure 2). The median Cα distance between residues where nonsynonymous changes are observed within bovine populations was calculated as 11.12 Å. In comparison, the median distance between all residues in the structure human CD2 extracellular region was 24.47 Å. In simulations, there were only 16/1000 cases where the median distance of seven randomly selected residues was lower than that of the median distance between the cattle nonsynonymous changes, indicating that these residues are non-randomly clustered on the structure of CD2 (P = 0.016).

Figure 2
figure 2

Two alternative perspectives of the 3D structure of human CD2 extracellular domain. Sites involved in CD58/CD48 binding are shown in blue. Codon sites (amino acids) predicted to be subject to positive selection from inter-specific comparison are shown in red. Yellow sites indicate the orthologous positions of nonsynonymous SNPs occurring within cattle populations. The structure was displayed using RasMol v2.

PolyPhen [19] predicts two substitutions, K180M and S189K, as "probably" and "possibly" damaging, respectively. The SIFT program [18] also predicts K180M, plus T156M, and V197M as likely to affect protein function. These variants each form part of the cluster haplotype of coding changes in exon 3 that are situated together on the surface of the extracellular domain of the protein.


The close proximity and subsequent migrations of humans and animals, following domestication, caused fundamental changes in the epidemiology of many infectious diseases [2]. Different cattle populations have been exposed differently to a variety of pathogens, such as rinderpest in early Asian populations [39] or trypanosomasis in sub-Saharan African populations [40]. Screening the genomes of domesticated animals for signatures of selection, particularly in relevant immune genes, may help pinpoint those genes that have been most important in the evolution of differential disease susceptibilities.

We queried whether six genes which have signatures suggestive of positive selection between species may have experienced continued selective pressure at the population level [6, 7, 41]. Notably, these genes were also chosen on the basis of having been previously implicated in susceptibility to bacterial, parasitic or viral infections. We re-sequenced each in 39 individuals comprising three distinct continental cattle populations. Several test statistics within the compilation of locus-population combinations yield evidence that recent adaptive change may have imprinted patterns of diversity around these genes. There is varying (but incomplete) concordance across tests. This is unsurprising as the statistics examine different aspects of locus diversity and are sensitive to selective processes that vary, for example in intensity or time depth [42]. The results fall into two interesting categories: first, outlying test statistics of allele frequency spectra in cytokines that imply the imprint of adaptive change but which are not associated with appreciable protein sequence polymorphism; and secondly, receptors which display coding variation suggestive of allelic functional variation, but which give less significant summary statistic values.

Of the latter loci, TYROBP gives the strongest statistical departure from neutrality, and has a derived nonsynonymous mutation almost fixed within European and African B. taurus cattle with predicted impact. In humans, ART4 encodes the Dombrock blood group glycoprotein [31]. It has been speculated that the polymorphic nature of this molecule in humans may be driven by an evolutionary advantage due to mutation of the RGD motif that may influence adhesion of the malaria parasite to erythrocytes [31]. Cattle do not appear to be polymorphic at the position responsible for the Doa/Dob antigenicity in humans, but are polymorphic at the position associated with the Joa antigen, suggesting a structural or functional impact [32]. However, of these three, the most interesting evidence arises with CD2.

CD2 is a natural killer-cell and T-cell surface molecule that stabilizes interactions with antigen-presenting cells [34]. CD2 diversity gives outlying scores for Fay and Wu's H in both European and African B. taurus under a domestication-modeled simulation and within Europe under a standard coalescent model. Taken alone, these scores are not secure indications of recent adaptive evolution at this locus, but a closer look at the pattern of derived polymorphism underlying them gives several layers of additional support.

Firstly, CD2 was found to have a high number (eleven) of nonsynonymous changes within cattle (three between cattle and the outgroup species, plains bison). The majority of these are localized in the extracellular domain of the molecule, a region with a strong signal of positive selection (from inter-specific comparison), which may be driven by a host-pathogen genetic conflict [4]. Secondly, six of these changes are derived, located within a narrow window within exon 3 and define a divergent haplotype that appears at appreciable frequency in each of the three sampled populations. This region yields a signal of selection in sliding window analysis and has a peak of sequence diversity. Thirdly this genetically divergent haplotype is likely also to be functionally divergent. Three of these clustered substitutions are predicted to have phenotypic impact from biochemical and comparative information. Lastly, examination of a predicted 3D structure of the extracellular domain (Figure 2) showed that the seven substitutions clustered spatially – a configuration confirmed to be unlikely to be random by bootstrapping.

Taken together, these observations provide population genetic inference of adaptive evolution at the CD2 locus, possibly driven by receptor-pathogen interactions. Note that several similar receptors to CD2 are known to be exploited by viruses in order to invade the host cell; for example, the CD2 family member, SLAM, is bound by the measles virus [38].

Within the three cytokines sequenced, the strongest population genetic signatures of selection occur at the closely linked IL5 and IL13 loci within the African B. taurus population sample. Fay and Wu's H (IL5) plus Tajima's D, Fu and Li's D and F (IL13) values all show highly significant departures from the neutral model. Interestingly, these two loci comprise part of a cytokine cluster on chromosome 7 (syntenic with human chromosome region 5q31–33) which has key involvement in the Th2 mediated immune response and consequently, susceptibility to extracellular parasitic infections. There are well described and considerable parasitic challenges particular to African livestock; for example, non-indigenous breeds continue to be excluded from parts of West Africa by endemic trypanosomiasis [10]. Both of these loci show some reduction in nucleotide diversity within Africa (as does IL2) relative to the other populations but a focus of positive selection within each region is not easily localized from the data. Also, each interleukin sequenced showed a dearth of nonsynonymous change. IL5 and IL13 (along with IL4) are known to be influenced by cis factors over a 120 Kbp region and the focus of selection may be a remote regulatory sequence.


This study alludes to two generalities. First, we show some correspondence between inference of selection from intra-species variation and that from the million year-order comparison between species. This hints that the diversity that was a matter of life and death for bovine antecedents has also been important in more recent eras, and probably continues to matter today. Second, we see contrasting patterns of diversity between the three signaling molecules (IL2, IL5 and IL13), which have outlying allele frequency spectra but insubstantial nonsynonymous polymorphism and three cell-surface receptors (TYROBP, ART4 and CD2), which interact directly with the extracellular environment and were observed to have significant coding polymorphism, giving some suggestion of adaptive history. This difference accords with prior work that records an increase of protein coding sequence adaptive events at the cellular periphery in humans, with more constraint on change at those loci which are internal and more central in interacting networks [43].


Re-sequencing strategy

We have re-sequenced almost all coding sequence, plus extensive intronic and intergenic regions from the IL13 (NM_174089), IL2 (NM_180997), IL5 (NM_173922), TYROBP (NM_174627), ART4 (XM_583864), and CD2 (NM_001011676) genes, in a panel of 18 European Bos taurus, 10 African B. taurus and 11 B. indicus cattle. These three population samples were drawn from 16 breeds of cattle (see Additional file 2). Breeds with known admixture were not included, as admixture of the scale in these breeds is a confounding factor in tests of selection. The orthologous regions were also re-sequenced in the outgroup species, plains bison (Bison bison). Primer details and PCR conditions are available for each gene sequenced (see Additional file 3). Samples were sequenced on the forward and reverse strands by Agowa, Germany The regions re-sequenced for each gene are indicated figuratively (see Additional file 4). Sequence data from this article have been deposited with the EMBL/GenBank Data Libraries under accession no's. [GenBank:EF118301–EF118739; EU914955–EU914988; EU914991–EU915048].

SNP detection and sequence analysis

Sequence assembly and SNP analysis was carried out using the Phred/Phrap/Consed/Polyphred pipeline Bases were called and sequences assembled into contigs using Phred v0.020425.c [44, 45] and Phrap v0.990319. Default parameters were used except for the forcelevel flag, which was set to a maximum value of 10, reducing the stringency of matching required for final contig assembly, and allowing outgroup sequences to be incorporated in subsequent SNP detection. Polyphred v.5.0 was used for polymorphism detection [46, 47]. Reads were compared to the UCSC genomic reference sequence, and the source flag was used to optimize heterozygote identification. SNPs were checked manually and verified by examining trace files in Consed [48]. Alignments of polymorphic positions for each gene are available (see Additional file 1). PHASE v.2.1 was used to reconstruct haplotypes [49]. Nucleotide positions with lower confidence in sequence quality, i.e. fewer than 20 individuals with high quality sequence were masked and excluded from further analysis.

Neutrality tests

We carried out several tests of neutrality based on allele frequency spectra derived from the re-sequencing data. These statistics were calculated over the entire re-sequenced regions (Table 1) using DnaSP [50, 51]. The Tajima's D test [15] compared the difference between the number of segregating sites, S, and the average number of pairwise differences, π. Negative values of Tajima's D are due to an excess of low frequency variants, which may be due to positive selection or population expansion.

Fu and Li's D and F statistics incorporate information from a genealogy constructed with the outgroup sequence, plains bison. External branches of this genealogy tend to be populated by more recent mutations. Fu and Li's D is based on the differences between Se, the total number of mutations in external branches of the genealogy, and S, the total number of mutations, while, their F test statistic is based on the differences between Se, the total number of mutations in external branches of the genealogy, and π [16].

The H statistic, developed by Fay and Wu, requires an outgroup sequence, and is based on the difference between two estimators of θ : π and θ H , which is based on the frequency of the derived variants [23]. Fu's Fs [17] is based on the probability of observing a greater number of alleles in the population than expected under neutrality. Either genetic hitch-hiking or recent population expansion can cause a significant result in this test. For all of these tests, negative values result from an excess of low frequency and/or derived polymorphisms, while a deficiency of low frequency alleles causes a shift towards positive values.

One problem with neutrality test statistics can be their sensitivity to different demographies. Consequently, to test the robustness of our findings, we examined an alternative demographic history for whole gene coalescent simulations. These were used to construct empirical null distributions of two complementary test statistics (Tajima's D and Fay and Wu's H) presuming a constant population size of 104 until 10,000 years ago, followed by an instant reduction to 103 and an exponential increase to a present day size of 105, a scenario consistent with the domestication and subsequent expansion of livestock. Simulations were carried out using MS [52].

Predictive polymorphism phenotyping

Predictive polymorphism phenotyping of nonsynonymous SNPs was carried out using the PolyPhen and Sorting Intolerant From Tolerant (SIFT) programs [18, 19]. PolyPhen assesses a number of parameters including, amino acid physiochemical properties, the level of sequence conservation in homologous proteins, and the proximity to structural or functional domains, to predict functionally significant nonsynonymous SNPs. Similarly, the SIFT program uses an alignment of orthologous or paralogous sequences to assess the usual level of sequence conservation at a polymorphic site and thus the potential functional importance of a substitution.

Non-random Clustering of CD2 nonsynonymous SNPs

The distance between the alpha carbon atoms of each residue (Cα distance) in the 3 dimensional structure of the human CD2 extracellular domain (1 HNF.pdb) was calculated using the BioShell program





ADP-ribosyltransferase 4


base pair

Cα distance:

the distance between the alpha carbon atoms of each residue


cluster of differentiation 2 (CD2)






killer-cell inhibitory receptor


natural killer


nucleotide diversity


ratio of nonsynonymous substitutions per nonsynonymous site to synonymous substitutions per synonymous site


segregating sites


the total number of mutations in external branches of the genealogy


Sorting Intolerant From Tolerant


signaling lymphocytic activation molecule


tyrosine kinase binding protein


University of California Santa Cruz


  1. Zeder MA, Emshwiller E, Smith BD, Bradley DG: Documenting domestication: the intersection of genetics and archaeology. Trends Genet. 2006, 22 (3): 139-155. 10.1016/j.tig.2006.01.007.

    Article  CAS  PubMed  Google Scholar 

  2. Pearce-Duvet JM: The origin of human pathogens: evaluating the role of agriculture and domestic animals in the evolution of human disease. Biol Rev Camb Philos Soc. 2006, 81 (3): 369-382. 10.1017/S1464793106007020.

    Article  PubMed  Google Scholar 

  3. White SN, Taylor KH, Abbey CA, Gill CA, Womack JE: Haplotype variation in bovine Toll-like receptor 4 and computational prediction of a positively selected ligand-binding domain. Proceedings of the National Academy of Sciences of the United States of America. 2003, 100 (18): 10364-10369. 10.1073/pnas.1333957100.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  4. Lynn DJ, Freeman AR, Murray C, Bradley DG: A genomics approach to the detection of positive selection in cattle: adaptive evolution of the T-cell and natural killer cell-surface protein CD2. Genetics. 2005, 170 (3): 1189-1196. 10.1534/genetics.104.039040.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  5. Hughes AL, Packer B, Welch R, Chanock SJ, Yeager M: High level of functional polymorphism indicates a unique role of natural selection at human immune system loci. Immunogenetics. 2005, 57 (11): 821-827. 10.1007/s00251-005-0052-7.

    Article  CAS  PubMed  Google Scholar 

  6. Akey JM, Eberle MA, Rieder MJ, Carlson CS, Shriver MD, Nickerson DA, Kruglyak L: Population history and natural selection shape patterns of genetic variation in 132 genes. PLoS Biol. 2004, 2 (10): e286-10.1371/journal.pbio.0020286.

    Article  PubMed Central  PubMed  Google Scholar 

  7. Nielsen R, Bustamante C, Clark AG, Glanowski S, Sackton TB, Hubisz MJ, Fledel-Alon A, Tanenbaum DM, Civello D, White TJ, et al: A scan for positively selected genes in the genomes of humans and chimpanzees. PLoS Biol. 2005, 3 (6): e170-10.1371/journal.pbio.0030170.

    Article  PubMed Central  PubMed  Google Scholar 

  8. Bradley DG, MacHugh DE, Cunningham P, Loftus RT: Mitochondrial diversity and the origins of African and European cattle. Proc Natl Acad Sci USA. 1996, 3 (10): 5131-5135. 10.1073/pnas.93.10.5131.

    Article  Google Scholar 

  9. Troy CS, MacHugh DE, Bailey JF, Magee DA, Loftus RT, Cunningham P, Chamberlain AT, Sykes BC, Bradley DG: Genetic evidence for Near-Eastern origins of European cattle. Nature. 2001, 410 (6832): 1088-1091. 10.1038/35074088.

    Article  CAS  PubMed  Google Scholar 

  10. Freeman AR, Hoggart CJ, Hanotte O, Bradley DG: Assessing the relative ages of admixture in the bovine hybrid zones of Africa and the Near East using X chromosome haplotype mosaicism. Genetics. 2006, 173 (3): 1503-1510. 10.1534/genetics.105.053280.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  11. McKenzie GJ, Bancroft A, Grencis RK, McKenzie AN: A distinct role for interleukin-13 in Th2-cell-mediated immune responses. Curr Biol. 1998, 8 (6): 339-342. 10.1016/S0960-9822(98)70134-4.

    Article  CAS  PubMed  Google Scholar 

  12. Gibbs BF: Human basophils as effectors and immunomodulators of allergic inflammation and innate immunity. Clin Exp Med. 2005, 5 (2): 43-49. 10.1007/s10238-005-0064-5.

    Article  CAS  PubMed  Google Scholar 

  13. Bancroft AJ, McKenzie AN, Grencis RK: A critical role for IL-13 in resistance to intestinal nematode infection. J Immunol. 1998, 160 (7): 3453-3461.

    CAS  PubMed  Google Scholar 

  14. Fallon PG, Richardson EJ, McKenzie GJ, McKenzie AN: Schistosome infection of transgenic mice defines distinct and contrasting pathogenic roles for IL-4 and IL-13: IL-13 is a profibrotic agent. J Immunol. 2000, 164 (5): 2585-2591.

    Article  CAS  PubMed  Google Scholar 

  15. Tajima F: Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989, 123 (3): 585-595.

    PubMed Central  CAS  PubMed  Google Scholar 

  16. Fu YX, Li WH: Statistical tests of neutrality of mutations. Genetics. 1993, 33 (3): 693-709.

    Google Scholar 

  17. Fu YX: Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics. 1997, 147 (2): 915-925.

    PubMed Central  CAS  PubMed  Google Scholar 

  18. Ng PC, Henikoff S: Predicting deleterious amino acid substitutions. Genome Res. 2001, 11 (5): 863-874. 10.1101/gr.176601.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  19. Ramensky V, Bork P, Sunyaev S: Human non-synonymous SNPs: server and survey. Nucleic Acids Res. 2002, 30 (17): 3894-3900. 10.1093/nar/gkf493.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  20. Loots GG, Locksley RM, Blankespoor CM, Wang ZE, Miller W, Rubin EM, Frazer KA: Identification of a coordinate regulator of interleukins 4, 13, and 5 by cross-species sequence comparisons. Science. 2000, 288 (5463): 136-140. 10.1126/science.288.5463.136.

    Article  CAS  PubMed  Google Scholar 

  21. Ellis MK, Zhao ZZ, Chen HG, Montgomery GW, Li YS, McManus DP: Analysis of the 5q31 33 locus shows an association between single nucleotide polymorphism variants in the IL-5 gene and symptomatic infection with the human blood fluke, Schistosoma japonicum. J Immunol. 2007, 179 (12): 8366-8371.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  22. Gaffen SL, Liu KD: Overview of interleukin-2 function, production and clinical applications. Cytokine. 2004, 28 (3): 109-123. 10.1016/j.cyto.2004.06.010.

    Article  CAS  PubMed  Google Scholar 

  23. Fay JC, Wu CI: Hitchhiking under positive Darwinian selection. Genetics. 2000, 155 (3): 1405-1413.

    PubMed Central  CAS  PubMed  Google Scholar 

  24. Diefenbach A, Tomasello E, Lucas M, Jamieson AM, Hsia JK, Vivier E, Raulet DH: Selective associations with signaling proteins determine stimulatory versus costimulatory activity of NKG2D. Nat Immunol. 2002, 3 (12): 1142-1149. 10.1038/ni858.

    Article  CAS  PubMed  Google Scholar 

  25. Lanier LL, Corliss BC, Wu J, Leong C, Phillips JH: Immunoreceptor DAP12 bearing a tyrosine-based activation motif is involved in activating NK cells. Nature. 1998, 391 (6668): 703-707. 10.1038/35642.

    Article  CAS  PubMed  Google Scholar 

  26. Sjolin H, Tomasello E, Mousavi-Jazi M, Bartolazzi A, Karre K, Vivier E, Cerboni C: Pivotal role of KARAP/DAP12 adaptor molecule in the natural killer cell-mediated resistance to murine cytomegalovirus infection. J Exp Med. 2002, 195 (7): 825-834. 10.1084/jem.20011427.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  27. Mohamadzadeh M, Coberley SS, Olinger GG, Kalina WV, Ruthel G, Fuller CL, Swenson DL, Pratt WD, Kuhns DB, Schmaljohn AL: Activation of triggering receptor expressed on myeloid cells-1 on human neutrophils by marburg and ebola viruses. J Virol. 2006, 80 (14): 7235-7244. 10.1128/JVI.00543-06.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  28. Hao L, Nei M: Genomic organization and evolutionary analysis of Ly49 genes encoding the rodent natural killer cell receptors: rapid evolution by repeated gene duplication. Immunogenetics. 2004, 56 (5): 343-354. 10.1007/s00251-004-0703-0.

    Article  CAS  PubMed  Google Scholar 

  29. Sambrook JG, Bashirova A, Palmer S, Sims S, Trowsdale J, Abi-Rached L, Parham P, Carrington M, Beck S: Single haplotype analysis demonstrates rapid evolution of the killer immunoglobulin-like receptor (KIR) loci in primates. Genome Res. 2005, 15 (1): 25-35. 10.1101/gr.2381205.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  30. Storset AK, Slettedal IO, Williams JL, Law A, Dissen E: Natural killer cell receptors in cattle: a bovine killer cell immunoglobulin-like receptor multigene family contains members with divergent signaling motifs. Eur J Immunol. 2003, 33 (4): 980-990. 10.1002/eji.200323710.

    Article  CAS  PubMed  Google Scholar 

  31. Gubin AN, Njoroge JM, Wojda U, Pack SD, Rios M, Reid ME, Miller JL: Identification of the dombrock blood group glycoprotein as a polymorphic member of the ADP-ribosyltransferase gene family. Blood. 2000, 96 (7): 2621-2627.

    CAS  PubMed  Google Scholar 

  32. Reid ME: The Dombrock blood group system: a review. Transfusion. 2003, 43 (1): 107-114. 10.1046/j.1537-2995.2003.00283.x.

    Article  CAS  PubMed  Google Scholar 

  33. Brossay A, Hube F, Moreau T, Bardos P, Watier H: Porcine CD58: cDNA cloning and molecular dissection of the porcine CD58-human CD2 interface. Biochem Biophys Res Commun. 2003, 309 (4): 992-998. 10.1016/j.bbrc.2003.08.099.

    Article  CAS  PubMed  Google Scholar 

  34. Selvaraj P, Plunkett ML, Dustin M, Sanders ME, Shaw S, Springer TA: The T lymphocyte glycoprotein CD2 binds the cell surface ligand LFA-3. Nature. 1987, 326 (6111): 400-403. 10.1038/326400a0.

    Article  CAS  PubMed  Google Scholar 

  35. Shimojima M, Nishimura Y, Miyazawa T, Kato K, Nakamura K, Izumiya Y, Akashi H, Tohya Y: A feline CD2 homologue interacts with human red blood cells. Immunology. 2002, 105 (3): 360-366. 10.1046/j.0019-2805.2001.01371.x.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  36. Kato K, Koyanagi M, Okada H, Takanashi T, Wong YW, Williams AF, Okumura K, Yagita H: CD48 is a counter-receptor for mouse CD2 and is involved in T cell activation. J Exp Med. 1992, 176 (5): 1241-1249. 10.1084/jem.176.5.1241.

    Article  CAS  PubMed  Google Scholar 

  37. Dustin ML, Golan DE, Zhu DM, Miller JM, Meier W, Davies EA, Merwe van der PA: Low affinity interaction of human or rat T cell adhesion molecule CD2 with its ligand aligns adhering membranes to achieve high physiological affinity. J Biol Chem. 1997, 272 (49): 30889-30898. 10.1074/jbc.272.49.30889.

    Article  CAS  PubMed  Google Scholar 

  38. Wang J: Protein recognition by cell surface receptors: physiological receptors versus virus interactions. Trends Biochem Sci. 2002, 27 (3): 122-126. 10.1016/S0968-0004(01)02038-2.

    Article  CAS  PubMed  Google Scholar 

  39. Spinage CA: Cattle plague: a history. 2003, New York: Kluwer Academic/Plenum Publishers

    Chapter  Google Scholar 

  40. Gifford-Gonzalez D: Animal Disease Challenges to the Emergence of Pastoralism in Sub-Saharan Africa. African Archaeological Review. 2000, 17 (3): 95-139. 10.1023/A:1006601020217.

    Article  Google Scholar 

  41. Gilbert SL, Dobyns WB, Lahn BT: Genetic links between brain development and brain evolution. Nat Rev Genet. 2005, 6 (7): 581-590. 10.1038/nrg1634.

    Article  CAS  PubMed  Google Scholar 

  42. Nielsen R, Hellmann I, Hubisz M, Bustamante C, Clark AG: Recent and ongoing selection in the human genome. Nat Rev Genet. 2007, 8 (11): 857-868. 10.1038/nrg2187.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  43. Kim PM, Korbel JO, Gerstein MB: Positive selection at the protein network periphery: evaluation in terms of structural constraints and cellular context. Proc Natl Acad Sci USA. 2007, 104 (51): 20274-20279. 10.1073/pnas.0710183104.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  44. Ewing B, Green P: Base-calling of automated sequencer traces using phred. II. Error probabilities. Genome Res. 1998, 8 (3): 186-194.

    Article  CAS  PubMed  Google Scholar 

  45. Ewing B, Hillier L, Wendl MC, Green P: Base-calling of automated sequencer traces using phred. I. Accuracy assessment. Genome Res. 1998, 8 (3): 175-185.

    Article  CAS  PubMed  Google Scholar 

  46. Stephens M, Sloan JS, Robertson PD, Scheet P, Nickerson DA: Automating sequence-based detection and genotyping of SNPs from diploid samples. Nat Genet. 2006, 38 (3): 375-381. 10.1038/ng1746.

    Article  CAS  PubMed  Google Scholar 

  47. Nickerson DA, Tobe VO, Taylor SL: PolyPhred: automating the detection and genotyping of single nucleotide substitutions using fluorescence-based resequencing. Nucleic Acids Res. 1997, 25 (14): 2745-2751. 10.1093/nar/25.14.2745.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  48. Gordon D, Abajian C, Green P: Consed: a graphical tool for sequence finishing. Genome Res. 1998, 8 (3): 195-202.

    Article  CAS  PubMed  Google Scholar 

  49. Stephens M, Smith NJ, Donnelly P: A new statistical method for haplotype reconstruction from population data. Am J Hum Genet. 2001, 68 (4): 978-989. 10.1086/319501.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  50. Rozas J, Rozas R: DnaSP version 3: an integrated program for molecular population genetics and molecular evolution analysis. Bioinformatics. 1999, 15 (2): 174-175. 10.1093/bioinformatics/15.2.174.

    Article  CAS  PubMed  Google Scholar 

  51. Rozas J, Sanchez-DelBarrio JC, Messeguer X, Rozas R: DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics. 2003, 19 (18): 2496-2497. 10.1093/bioinformatics/btg359.

    Article  CAS  PubMed  Google Scholar 

  52. Hudson RR: Generating samples under a Wright-Fisher neutral model of genetic variation. Bioinformatics. 2002, 18 (2): 337-338. 10.1093/bioinformatics/18.2.337.

    Article  CAS  PubMed  Google Scholar 

Download references


Thanks to Sarah Connell and Lilian Lau for helping with manuscript preparation. This material is based on work supported by the Science Foundation Ireland under Grant No. 02-IN.1-B256.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Daniel G Bradley.

Additional information

Authors' contributions

ARF and DJL contributed equally to this work, both carried out data analysis and co-wrote the manuscript. ARF re-sequenced the genes studied in this project. CM assisted in SNP calling and in data analysis. DGB conceived of and oversaw the study and assisted in writing the manuscript.

Abigail R Freeman, David J Lynn contributed equally to this work.

Electronic supplementary material


Additional File 1:Alignments of polymorphic positions in re-sequenced bovine IL13, IL2, IL5, TYROBP, ART4 and CD2 genes.(XLS 292 KB)

Additional File 2:Breed and country of origin information for each bovine sample.(XLS 18 KB)

Additional File 3:Primer information for re-sequenced bovine genes.(XLS 44 KB)


Additional File 4: Genomic structures of the genes encoding bovine A) Interleukin 13 (IL13); B) IL2; C) IL5; D) TYROBP; E) ART4 and F) CD2. Exons are represented by shaded grey rectangles. Untranslated regions (UTRs) are similarly shown but are represented by a smaller height. Ns indicate regions not re-sequenced. Nonsynonymous SNPs ▲ (red), synonymous SNPs ■, and non-coding (intronic, upstream and downstream) SNPs | are shown. (TIFF 13 MB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Freeman, A.R., Lynn, D.J., Murray, C. et al. Detecting the effects of selection at the population level in six bovine immune genes. BMC Genet 9, 62 (2008).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: