Molecular evolution of bovine Toll-like receptor 2 suggests substitutions of functional relevance

Abstract

Background: There is accumulating evidence that polymorphism in Toll-like receptor (TLR) genes might be associated with disease resistance or susceptibility traits in livestock. Polymorphic sites affecting TLR function should exhibit signatures of positive selection, identified as a high ratio of non-synonymous to synonymous nucleotide substitutions (omega). Phylogeny based models of codon substitution based on estimates of omega for each amino acid position can therefore offer a valuable tool to predict sites of functional relevance. We have used this approach to identify such polymorphic sites within the bovine TLR2 genes from ten Bos indicus and Bos taurus cattle breeds. By analysing TLR2 gene phylogeny in a set of mammalian species and a subset of ruminant species we have estimated the selective pressure on individual sites and domains and identified polymorphisms at sites of putative functional importance.

Results: The omega were highest in the mammalian TLR2 domains thought to be responsible for ligand binding and lowest in regions responsible for heterodimerisation with other TLR-related molecules. Several positively-selected sites were detected in or around ligand-binding domains. However a comparison of the ruminant subset of TLR2 sequences with the whole mammalian set of sequences revealed that there has been less selective pressure among ruminants than in mammals as a whole. This suggests that there have been functional changes during ruminant evolution. Twenty newly-discovered non-synonymous polymorphic sites were identified in cattle. Three of them were localised at positions shaped by positive selection in the ruminant dataset (Leu227Phe, His305Pro, His326Gln) and in domains involved in the recognition of ligands. His326Gln is of particular interest as it consists of an exchange of differentially-charged amino acids at a position which has previously been shown to be crucial for ligand binding in human TLR2.

Conclusion: Within bovine TLR2, polymorphisms at amino acid positions 227, 305 and 326 map to functionally important sites of TLR2 and should be considered as candidate SNPs for immune related traits in cattle. A final proof of their functional relevance requires further studies to determine their functional effect on the immune response after stimulation with relevant ligands and/or their association with immune related traits in animals.

Figure 1

Average dN/dS ratio of a 20 AA sliding window along the protein sequence. The analysis is based on the complete dataset of 22 mammalian sequences. The predicted domain structure, position of polymorphic and positively selected (P > 0.95) sites is indicated. Sites displayed as "functional" refer to positions which were found to be crucial for ligand binding or heterodimerisation in human (additional file 7: Functional relevant TLR2 sites identified in former studies in human or mouse).

Figure 2

Schematic phylogenetic tree of the TLR2sequences. Numbers represent nodes and are discussed in the main text.